Drainage basin space living source pollution discharge amount calculation method based on POI distribution

Abstract

The invention discloses a drainage basin space living source pollution discharge amount calculation method based on POI distribution, and relates to the technical field of pollution discharge amount calculation. Obtaining research data; obtaining the number of town POIs and the number of rural POIs contained in different administrative regions in each small watershed; calculating to obtain the urban living source pollutant discharge amount and the rural living source pollutant discharge amount in different administrative regions in each small watershed; and the total space living pollution emission amount of a certain small watershed is obtained through the urban living source pollutant emission amount and the rural living source pollutant emission amount in different administrative regions in each small watershed. According to the drainage basin space living source pollution discharge amount calculation method based on POI distribution, through space distribution of POI data, the distribution characteristics of population can be reflected more accurately, and therefore the pollution source calculation precision is improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of pollution emission calculation, in particular to a watershed spatial life source pollution emission calculation method based on POI distribution. BACKGROUND

[0002] With the acceleration of urbanization, life source pollution has become one of the important factors of water environmental quality deterioration in the watershed. The traditional pollution source calculation method is mostly based on administrative division or data in statistical yearbook, which is difficult to accurately reflect the actual distribution of life sewage emission, especially at the watershed scale. This statistical method is often inconsistent with the water flow collection path, resulting in large deviation of the calculation result. Therefore, a pollution source calculation method based on spatial distribution characteristics is urgently needed to improve the accuracy and scientificity of pollution emission estimation. SUMMARY

[0003] The purpose of the present application is to provide a watershed spatial life source pollution emission calculation method based on POI distribution, to solve the problems raised in the background art, and to realize the fine calculation of life source pollution emission in the watershed through high-precision spatial data acquisition and analysis.

[0004] To achieve the above purpose, the present application provides a watershed spatial life source pollution emission calculation method based on POI distribution, comprising the following steps:

[0005] Step S1, using ArcSWAT to divide the study area into watersheds;

[0006] Step S2, obtaining research data according to Baidu POI big data interface service;

[0007] Step S3, obtaining the number of urban POIs and the number of rural POIs in each small watershed in different administrative areas;

[0008] Step S4, according to the number of urban POIs and the number of rural POIs, combining the life source pollution emission coefficient manual to calculate the urban life source pollution emission and the rural life source pollution emission in each small watershed in different administrative areas;

[0009] Step S5, obtaining the total spatial life pollution emission of a small watershed by the urban life source pollution emission and the rural life source pollution emission in different administrative areas in each small watershed.

[0010] Preferably, the research data in step S2 is POI data related to restaurant, shopping mall, school, hospital and population activity in the study area.

[0011] Preferably, the number of urban POIs and the number of rural POIs in step S3 are obtained by histogram in ArcGIS.

[0012] Preferably, the amount of urban life source pollutant emissions in a certain administrative region in step S4 is calculated by the following formula:

[0013] The amount of urban life pollutant emissions in a certain administrative region = the number of POIs in a certain administrative region × the amount of comprehensive life water consumption of each POI in a certain administrative region × the pollution conversion coefficient of a certain administrative region × 365 / 1000 - the amount of life sewage reuse of all sewage treatment plants in a certain administrative region.

[0014] Preferably, the amount of rural life source pollutant emissions in a certain administrative region in step S4 is calculated by the following formula:

[0015] The amount of rural life sewage emissions in a certain administrative region = the number of POIs in a certain administrative region × the sewage emission coefficient of each POI in a certain administrative region × 365 / 1000.

[0016] Preferably, the total amount of life source pollutant emissions in a certain small watershed in step S5 is obtained by summing the total amount of life source pollutant emissions in a plurality of administrative regions in the small watershed, and is calculated by the following formula:

[0017] Q L =∑Q i ;

[0018] In the formula, Q L represents the total amount of life source pollutant emissions in a certain small watershed, ∑ represents a summation symbol, and Q i represents the total amount of life source pollutant emissions in a certain administrative region.

[0019] Preferably, the total amount of life source pollutant emissions Q i in a certain administrative region is calculated by the following formula:

[0020] Q i =A cp ×e c +A np ×e n ;

[0021] In the formula, Q i represents the total amount of life source pollutant emissions in a certain region, with the unit being tons / year; A cp represents the total number of POIs in a certain administrative region, with the unit being people; e c represents the life source pollutant coefficient of the population in a certain administrative region, with the unit being tons / person; A np represents the total number of POIs in a certain administrative region, with the unit being people; and e n represents the life source pollutant coefficient of the population in a certain administrative region, with the unit being tons / person.

[0022] Therefore, the application adopts the above-mentioned basin spatial life source pollution discharge calculation method based on POI distribution, and has the following beneficial effects:

[0023] (1) High-precision spatial distribution analysis: through the spatial distribution of POI data, the distribution characteristics of the population can be more accurately reflected, thereby improving the precision of pollution source calculation.

[0024] (2) Consistent with the water flow collection path: the traditional method is to count the population in administrative areas, and the present method makes the population distribution consistent with the water flow collection path through spatial distribution, thereby improving the rationality of the calculation results.

[0025] (3) Suitable for large-scale basin analysis: the present method is not only suitable for small basin pollution source calculation, but also suitable for large-scale basin pollution source calculation, and has good expansibility.

[0026] The technical solutions of the present application will be further described in detail below with reference to the drawings and examples. BRIEF DESCRIPTION OF DRAWINGS

[0027] Figure 1 The POI data distribution diagram of the POI distribution-based basin spatial life source pollution discharge calculation method embodiment of the present application. DETAILED DESCRIPTION

[0028] The technical solutions of the present application will be further described in detail below with reference to the drawings and examples.

[0029] Unless otherwise defined, the technical terms or scientific terms used in the present application should be understood as the usual meanings understood by those skilled in the art to which the present application belongs. The "first", "second" and similar words used in the present application do not represent any order, quantity or importance, but are only used to distinguish different components. "Include" or "contain" and similar words mean that the elements or objects before the word cover the elements or objects listed after the word and their equivalents, without excluding other elements or objects. "Connected" or "connected" and similar words are not limited to physical or mechanical connection, but can include electrical connection, whether direct or indirect. "Up", "down", "left", "right" and the like are only used to represent relative positional relationship, when the absolute position of the described object changes, the relative positional relationship may also change accordingly.

[0030] EMBODIMENT

[0031] Please refer to Figure 1 The present application provides a POI distribution-based basin spatial life source pollution discharge calculation method, comprising the following steps:

[0032] Step S1, using ArcSWAT to divide the study area into basins.

[0033] Step S2, according to the Baidu POI big data interface service, the POI data of the research area restaurant, shopping mall, school, hospital and population activity closely related are obtained.

[0034] Step S3, the histogram in ArcGIS is used to obtain the number of urban POIs and the number of rural POIs in each small watershed in different administrative regions.

[0035] Step S4, according to the number of urban POIs and the number of rural POIs, the urban life source pollutant discharge and the rural life source pollutant discharge in each small watershed in different administrative regions are calculated respectively by combining the life source and discharge coefficient manual.

[0036] The urban life source pollutant discharge in a certain administrative region is calculated by the following formula:

[0037] The urban life pollutant discharge of a certain administrative region = the number of urban POIs in a certain administrative region × the urban POI comprehensive life water consumption coefficient of a certain administrative region × the pollution coefficient of a certain administrative region × 365 / 1000 - the life sewage reuse amount of all sewage treatment plants in a certain administrative region.

[0038] The rural life source pollutant discharge in a certain administrative region is calculated by the following formula:

[0039] The rural life pollutant discharge of a certain administrative region = the number of rural POIs in a certain administrative region × the rural POI sewage discharge coefficient of a certain administrative region × 365 / 1000.

[0040] Step S5, the total space life pollution discharge of a certain small watershed is obtained by the urban life source pollutant discharge and the rural life source pollutant discharge in different administrative regions in each small watershed. The total space life pollution discharge of a certain small watershed is obtained by summing the life source pollutant discharge in a plurality of administrative regions in a certain small watershed, which is calculated by the following formula:

[0041] Q L =∑Q i ;

[0042] In the formula, Q L represents the total space life source pollution discharge of a certain small watershed, ∑ represents the summation symbol, and Q i represents the total life source pollutant discharge in a certain administrative region.

[0043] The total life source pollutant discharge Q i in a certain administrative region is calculated by the following formula:

[0044] Q i =A cp ×e c +Anp x e n ;

[0045] In the formula, Q i represents the total amount of life source pollutant emissions in a certain area, in tons / year; A cp represents the total number of POIs in a certain administrative district, in units of people; e c represents the life source pollutant coefficient of the population in a certain administrative district, in tons / person; A np represents the total number of POIs in a certain administrative district, in units of people; e n represents the life source pollutant coefficient of the population in a certain administrative district, in tons / person.

[0046] The technical solutions of the present application are further illustrated by the following examples.

[0047] The administrative districts a and b are taken as the research area, and ArcSWAT is used for watershed division. The sub-watersheds a-001 and b-001 in the shadowed part are taken as the research objects. The POI data closely related to population activities are obtained through the Baidu POI big data interface service, and the number of urban POIs and the number of rural POIs in each small watershed in different administrative districts are obtained by using the histogram in ArcGIS.

[0048] Among them, the POI density in the urban area is usually high and the distribution is relatively concentrated; the POI density in the rural area is low and the distribution is relatively dispersed. The POI density can be used to intuitively see the intensity and spatial distribution pattern of POI. As shown in Figure 1 The POIs with a distance of more than two and less than three are considered as rural POIs, and the rest are urban POIs. In the a-001 research area, there are 6 rural POIs and 10 urban POIs, and in the b-001 research area, there are 5 rural POIs and 6 urban POIs. In this example, COD, NH3-H, TN and TP are selected for calculation according to the calculation formula of urban and rural life pollutant emissions, because they not only have clear coefficients for calculation in the life source pollutant emission coefficient manual, but also can more comprehensively reflect the pollution status of domestic sewage.

[0049] The amount of rural life pollutant COD, NH3-H, TN and TP in the sub-basin a-001 is 53.13, 4.12, 6.53 and 0.48 t / a respectively, and the amount of urban life pollutant COD, NH3-H, TN and TP is 214.13, 20.53, 28.22 and 2.69 t / a respectively; the amount of rural life pollutant COD, NH3-H, TN and TP in the sub-basin b-001 is 44.27, 3.43, 5.44 and 0.40 t / a respectively, and the amount of urban life pollutant COD, NH3-H, TN and TP is 128.48, 12.32, 16.93 and 1.61 t / a respectively.

[0050] The total amount of COD, NH3-H, TN and TP of the life source pollution discharge of the basin is 440.01, 40.40, 57.12 and 5.18 t / a respectively.

[0051] The above results can improve the accuracy of the calculation of the amount of pollution source discharge of the basin, and provide a reference for formulating a scientific control scheme.

[0052] Therefore, the above-mentioned calculation method of the amount of life source pollution discharge of the basin based on the POI distribution can more accurately reflect the distribution characteristics of the population by the spatial distribution of the POI data, so as to improve the accuracy of the calculation of the pollution source; by the spatial distribution, the population distribution is consistent with the water flow collection path, so as to improve the rationality of the calculation result; the method is not only suitable for small basins, but also suitable for the calculation of the pollution source of large-scale basins, and has good expansibility.

[0053] Finally, it should be noted that: the above examples are only used to illustrate the technical solutions of the present application but not to limit it, although the present application has been described in detail with reference to the preferred embodiments, those skilled in the art should understand that the technical solutions of the present application can still be modified or replaced by equivalents, and these modifications or equivalent replacements cannot make the modified technical solutions deviate from the spirit and scope of the technical solutions of the present application.

Claims

1.A method for calculating spatial life source pollution discharge of a watershed based on POI distribution, characterized in that, The method comprises the following steps: Step S1, using ArcSWAT to divide the study area into basins; Step S2, obtaining research data according to the Baidu POI big data interface service; Step S3, obtaining the number of urban POIs and the number of rural POIs in each small watershed in different administrative regions; Step S4, according to the number of urban POIs and the number of rural POIs, combined with the life source pollution coefficient manual, the urban life source pollutant emissions and the rural life source pollutant emissions in each small watershed in different administrative regions are calculated respectively; Step S5, obtaining the total spatial life pollution emission of a small watershed through the urban life source pollutant emissions and the rural life source pollutant emissions in different administrative regions in the small watershed. 2.The POI distribution-based calculation method of spatial life source pollution emission of a river basin according to claim 1, characterized in that: The research data in step S2 is the POI data related to the activities of restaurants, shopping malls, schools, hospitals and population in the study area. 3.The POI distribution-based calculation method of spatial life source pollution emission of a river basin according to claim 2, characterized in that: The number of urban POIs and the number of rural POIs in step S3 are obtained by histogram in ArcGIS. 4.The method of claim 3, wherein, The urban life source pollutant emissions in a certain administrative region in step S4 are calculated by the following formula: Urban life pollutant emissions in a certain administrative region = urban POI number in a certain administrative region × urban POI comprehensive water consumption coefficient in a certain administrative region × a certain administrative region pollution coefficient × 365 / 1000 - all the life sewage reuse amount of sewage treatment plants in a certain administrative region. 5.The POI distribution-based calculation method of spatial life source pollution emission of a river basin according to claim 4, characterized in that, The rural life source pollutant emissions in a certain administrative region in step S4 are calculated by the following formula: Rural life sewage emissions in a certain administrative region = rural POI number in a certain administrative region × rural POI sewage emission coefficient in a certain administrative region × 365 / 1000. 6.The POI distribution-based calculation method of spatial life source pollution emission of a river basin according to claim 5, characterized in that: The total spatial life pollution emission of a small watershed in step S5 is obtained by summing the total life source pollutant emissions in multiple administrative regions in the small watershed, which is calculated by the following formula: Q L =∑Q i ; In the formula, Q L represents the total amount of non-point source pollution in a small watershed, and Σ represents the summation symbol. i represents the total amount of non-point source pollution in a certain administrative region. 7.The method of claim 6, wherein, The total amount of life source pollutant discharge Q in the certain administrative region i The total amount of life source pollutant discharge Q in the certain administrative region i The total amount of life source pollutant discharge Q in the certain administrative region Q i = A cp x e c + A np x e n ; In the formula, Q i represents the total amount of life source pollutant emissions in a certain area, in tons / year; A cp represents the total number of POIs in a certain administrative district, in people; e c represents the coefficient of source pollution of urban population in a certain administrative region, with unit of ton / person; A np represents the total number of POIs in rural areas of a certain administrative region, with unit of person; e n The coefficient of rural population life source pollutants in a certain administrative region is represented, with the unit of tons / person.

Images