Estimation method for high-resolution large-scale domestic sewage and domestic sewage pollutant output

Abstract

The invention provides a high-resolution large-scale domestic sewage and domestic sewage pollutant output estimation method. The method comprises the following steps: acquiring open source data of a target area; dividing a target area into a plurality of independent calculation areas based on the spatial resolution of the open source data, and processing the open source data to obtain a parameter value of each independent calculation area; and based on the parameter value of each independent calculation area and a preset calculation formula, calculating the domestic sewage output and the domestic sewage pollutant output of each independent calculation area. According to the method, high-resolution and large-scale estimation of the domestic sewage output and the domestic sewage pollutant output can be realized, and a solid data support is provided for subsequent sewage treatment resource optimization.

Description

Technical Field

[0001] The present invention relates to the technical field of domestic sewage estimation statistics, and more specifically, to a method for estimating high-resolution, large-scale domestic sewage and the amount of domestic sewage pollutants produced. Background Art

[0003] Estimating the amount of domestic sewage and pollutants generated is the basis for the allocation of sewage treatment resources. High-resolution and high-precision estimation of domestic sewage and pollutants is conducive to more refined resource allocation optimization. Common sewage estimation methods include per capita comprehensive index method, unit construction land comprehensive index method, unit sub-item construction land index method, and industrial output value prediction method. Some scholars also construct the relationship between POI density and sewage volume based on POI (Point of Information) data; estimate grid sewage discharge distribution based on night light images; estimate BOD (Biochemical oxygen demand) load of grids in the basin based on hydrological models; and estimate the production, collection, treatment and reuse of global grid sewage based on national data downscaling. Summary of the invention

[0004] In order to overcome the problem that the grid resolution estimated by the prior art is difficult to meet the demand for fine allocation of sewage treatment resources, the present invention provides a high-resolution, large-scale method for estimating the amount of domestic sewage and domestic sewage pollutants generated.

[0005] In order to solve the above technical problems, the technical solution of the present invention is as follows:

[0006] The present invention proposes a high-resolution and large-scale method for estimating the amount of domestic sewage and domestic sewage pollutants generated, comprising the following steps:

[0007] Obtain open source data on the target area;

[0008] Dividing the target area into multiple independent calculation areas based on the spatial resolution of the open source data, and processing the open source data to obtain parameter values ​​of each independent calculation area;

[0009] Based on the parameter values ​​of each independent calculation area and the preset calculation formula, the amount of domestic sewage generated and the amount of domestic sewage pollutants generated in each independent calculation area are calculated.

[0010] Preferably, the open source data of the target area include the land use type, population, provincial administrative region, municipal administrative region, pollution reduction coefficient, pollution concentration coefficient, and rural sewage discharge coefficient of the target area.

[0011] Preferably, dividing the target area into a plurality of independent calculation areas based on the spatial resolution of the open source data, and processing the open source data to obtain the parameter value of each independent calculation area includes:

[0012] The spatial resolution of the population in the open source data is used as the spatial resolution of each independent calculation area. The study area is divided into multiple square grids according to the spatial resolution of each independent calculation area, and each square grid is used as an independent calculation area.

[0013] The open source data of each independent computing area is unified in coordinate system, resampled, cropped and outlier processed to realize the division of open source data into each independent computing area as the parameter value of each independent computing area.

[0014] Preferably, the parameter values ​​of each independent calculation area include population parameter values, space parameter values, and water consumption parameter values.

[0015] Preferably, the population parameter value includes the population size of the independent calculation area.

[0016] Preferably, the spatial parameter values ​​include the land use type, province and city of the independent calculation area.

[0017] Preferably, the water consumption parameter value includes the pollution concentration coefficient, pollution reduction coefficient or rural sewage discharge coefficient of the independent calculation area.

[0018] Preferably, based on the parameter value of each independent calculation area and a preset calculation formula, calculating the amount of domestic sewage generated and the amount of domestic sewage pollutants generated in each independent calculation area includes:

[0019] When the land use type of the independent calculation area is urban, the amount of urban domestic sewage and the amount of urban domestic sewage pollutants generated are calculated based on the population, pollution reduction coefficient, pollution concentration coefficient and the preset calculation formula for the amount of urban domestic sewage and the amount of urban domestic sewage pollutants generated;

[0020] When the land use type of the independent calculation area is rural, the rural domestic sewage generation and rural domestic sewage pollutant generation are calculated based on the population size, rural sewage discharge coefficient, pollution concentration coefficient and the preset rural domestic sewage generation and rural domestic sewage pollutant generation calculation formula.

[0021] Preferably, the calculation formula for the amount of urban domestic sewage generated and the amount of pollutants generated by urban domestic sewage includes:

[0022] WW=POP×PUW p ×α p

[0023] POLi =WW×β p,i

[0024] Wherein, WW is the sewage generation in the independent calculation area (liters / day); POP is the population in the independent calculation area (people); p represents the province; PUW p is the per capita comprehensive water consumption coefficient of the province where the independent calculation area is located (liter / (person / day)); α p is the pollution reduction coefficient of the province where the independent calculation area is located; i represents the type of pollutant, including chemical oxygen demand, ammonia nitrogen, total nitrogen, and total phosphorus; POL i is the amount of pollutant i produced in the independent calculation area (mg / day); β p,i It is the pollution concentration coefficient (mg / L) of pollutant i in the province where the independent calculation area is located.

[0025] Preferably, the calculation formula for the amount of rural domestic sewage generated and the amount of pollutants generated in rural domestic sewage includes:

[0026] WW=POP×γ c

[0027] POL i =WW×β c,i

[0028] Where c represents a city; γ c β is the rural sewage discharge coefficient of the city where the independent calculation area is located (liter / (person / day)); c,i It is the pollution concentration coefficient (g / (person·L)) of pollutant i in the city where the independent calculation area is located.

[0029] Compared with the prior art, the technical solution of the present invention has the following beneficial effects:

[0030] The present invention proposes a high-resolution and large-scale estimation method for the amount of domestic sewage and domestic sewage pollutants generated. First, the open source data of the target area is obtained; the target area is divided into multiple independent calculation areas based on the spatial resolution of the open source data, and the open source data is processed to obtain the parameter value of each independent calculation area; based on the parameter value of each independent calculation area and the preset calculation formula, the amount of domestic sewage generated and the amount of domestic sewage pollutants generated in each independent calculation area are calculated. The present invention can achieve high-resolution and large-scale estimation of the amount of domestic sewage generated and the amount of domestic sewage pollutants generated, and provide solid data support for the subsequent optimization of sewage treatment resources. BRIEF DESCRIPTION OF THE DRAWINGS

[0031] Figure 1 This is a schematic diagram of the process of the estimation method described in Example 1;

[0032] Figure 2This is a schematic diagram of the overall scheme of the estimation method described in Example 2;

[0033] Figure 3 This is a flow chart of the estimation method described in Example 2;

[0034] Figure 4 This is a schematic diagram of the overall scheme of the estimation method described in Example 3;

[0035] Figure 5 This is a diagram showing the sewage volume estimation result described in Example 3;

[0036] Figure 6 This is a diagram showing the estimation results of the four types of sewage pollutants described in Example 3. DETAILED DESCRIPTION

[0037] The drawings are for illustrative purposes only and should not be construed as limiting the present patent;

[0038] In order to better illustrate the present embodiment, some parts in the drawings may be omitted, enlarged or reduced, and do not represent the size of the actual product;

[0039] It is understandable to those skilled in the art that some well-known structures and their descriptions may be omitted in the drawings.

[0040] The technical solution of the present invention is further described below in conjunction with the accompanying drawings and embodiments.

[0041] Example 1

[0042] This embodiment provides a high-resolution and large-scale method for estimating the amount of domestic sewage and domestic sewage pollutants generated. Figure 1 As shown, the following steps are included:

[0043] Obtain open source data on the target area;

[0044] Dividing the target area into multiple independent calculation areas based on the spatial resolution of the open source data, and processing the open source data to obtain parameter values ​​of each independent calculation area;

[0045] Based on the parameter values ​​of each independent calculation area and the preset calculation formula, the amount of domestic sewage generated and the amount of domestic sewage pollutants generated in each independent calculation area are calculated.

[0046] In the specific implementation process, firstly, the open source data of the target area is obtained, then the target area is divided into multiple independent calculation areas based on the spatial resolution of the open source data, and the open source data is processed to obtain the parameter value of each independent calculation area. Finally, based on the parameter value of each independent calculation area and the preset calculation formula, the amount of domestic sewage and the amount of domestic sewage pollutants generated in each independent calculation area are calculated.

[0047] Example 2

[0048] The present invention proposes a high-resolution and large-scale method for estimating the amount of domestic sewage and domestic sewage pollutants generated, comprising the following steps:

[0049] Obtain open source data of the target area, including land use type, population, provincial administrative region, municipal administrative region, pollution reduction coefficient, pollution concentration coefficient, and rural sewage discharge coefficient of the target area, and store the data on the GIS platform.

[0050] like Figure 2 As shown in FIG. 1 , it is a general scheme diagram of the method for estimating the generation of domestic sewage and domestic sewage pollutants with high resolution and large scale according to the present embodiment. The target area is divided into multiple independent calculation areas based on the spatial resolution of the open source data, and the open source data is processed to obtain the parameter value of each independent calculation area, specifically:

[0051] The spatial resolution of the population in the open source data is used as the spatial resolution of each independent calculation area. The study area is divided into multiple square grids according to the spatial resolution of each independent calculation area, and each square grid is used as an independent calculation area.

[0052] like Figure 3 As shown, this is a flowchart for estimating the amount of domestic sewage and domestic sewage pollutants generated with high resolution and large scale as described in this embodiment. The open source data of each independent calculation area are subjected to a unified coordinate system, resampled, cropped and outlier processing to realize the division of the open source data into each independent calculation area as the parameter value of each independent calculation area.

[0053] The parameter values ​​of each independent calculation area include population parameter values, spatial parameter values, and water consumption parameter values. The population parameter values ​​include the population size of the independent calculation area, the spatial parameter values ​​include the land use type, province, and city of the independent calculation area, and the water consumption parameter values ​​include the pollution concentration coefficient, pollution reduction coefficient, or rural sewage discharge coefficient of the independent calculation area.

[0054] Based on the parameter values ​​of each independent calculation area and the preset calculation formula, the amount of domestic sewage and domestic sewage pollutants generated in each independent calculation area are calculated, specifically:

[0055] When the land use type of the independent calculation area is urban, the amount of urban domestic sewage and the amount of urban domestic sewage pollutants generated are calculated based on the population, pollution reduction coefficient, pollution concentration coefficient and the preset calculation formula for the amount of urban domestic sewage and the amount of urban domestic sewage pollutants generated;

[0056] WW=POP×PUW p ×α p

[0057] POL i =WW×β p,i

[0058] Wherein, WW is the sewage generation in the independent calculation area (liters / day); POP is the population in the independent calculation area (people); p represents the province; PUW p is the per capita comprehensive water consumption coefficient of the province where the independent calculation area is located (liter / (person / day)); α p is the pollution reduction coefficient of the province where the independent calculation area is located; i represents the type of pollutant, including chemical oxygen demand, ammonia nitrogen, total nitrogen, and total phosphorus; POL i is the amount of pollutant i produced in the independent calculation area (mg / day); β p,i It is the pollution concentration coefficient (mg / L) of pollutant i in the province where the independent calculation area is located.

[0059] When the land use type of the independent calculation area is rural, the rural domestic sewage generation and rural domestic sewage pollutant generation are calculated based on the population size, rural sewage discharge coefficient, pollution concentration coefficient and the preset rural domestic sewage generation and rural domestic sewage pollutant generation calculation formula.

[0060] WW=POP×γ c

[0061] POL i =WW×β c,i

[0062] Where c represents a city; γ c β is the rural sewage discharge coefficient of the city where the independent calculation area is located (liter / (person / day)); c,i It is the pollution concentration coefficient (g / (person·L)) of pollutant i in the city where the independent calculation area is located.

[0063] This embodiment also visualizes the amount of domestic sewage and domestic sewage pollutants generated in each independent calculation area by redefining the projection to make it consistent with the projection of the original data, and through geographic alignment and resampling, making it consistent with the position of the original data. Finally, the color expression of the numerical value is adjusted and legends, compasses, scales and other elements are added to fully display the results.

[0064] Example 3

[0065] Based on Example 2, this example takes the Dongjiang River Basin as the target area to calculate the amount of domestic sewage and domestic sewage pollutants generated in the Dongjiang River Basin. The Dongjiang River is one of the main water systems of the Pearl River. It originates from Yajibo, Xunwu County, Jiangxi Province. After entering Guangdong, it flows through Heyuan and Huizhou. It merges into the Pearl River Delta network river area in Shilong Town, Dongguan City, and flows into the Shiziyang in the southeast of Yudong Lianwei, Huangpu District, Guangzhou City. It flows into the sea through Humen. The catchment area is 35,340 square kilometers, accounting for 5.96% of the total area of ​​the Pearl River Basin and about 24.3% of the Pearl River Basin area in Guangdong. There are two large reservoirs in the Dongjiang River Basin, the Xinfengjiang Reservoir and the Fengshuba Reservoir, as well as the Dongshen Water Supply Project that supplies water to Hong Kong. As a water transport channel between Jiangxi and Guangdong and connecting the Pearl River Delta in eastern Guangdong, the Dongjiang River has promoted the prosperity of towns along the line. The overall process of implementation can be seen. Figure 4 .

[0066] The specific process is described as follows:

[0067] Step 1: Data download and data preprocessing

[0068] The 100-m population distribution data comes from the document "A100-m gridded population dataset of China's seventh census using ensemble learning and geospatial big data", the 30-m land use type comes from the Resource and Environmental Science and Data Platform of the Chinese Academy of Sciences, the provincial and municipal administrative boundaries come from Tiandi Map, and the pollution reduction coefficient, pollution concentration coefficient, and rural sewage discharge coefficient come from the "Handbook of Pollution Production and Discharge Coefficients for Domestic Sources".

[0069] First, the population distribution data, land use data, and provincial and municipal administrative boundaries are uniformly projected into the projection coordinate system WGS_1984_Web_Mercator_Auxiliary_Sphere. Then, the land use 30m is resampled to 100m through NEAREST to make it the same resolution as the population. Then, the province and city are converted to a 100m resolution raster so that they can participate in the calculation. While converting to a 100m resolution, the capture grid in the environment settings is set to the population distribution data to unify the grid position. Finally, the above data is clipped according to the shape of the Dongjiang River Basin.

[0070] Step 2: Outlier processing

[0071] Due to the different data sources, 100% accuracy and fit cannot be guaranteed. Although the population of some grids is greater than 0, the land use type is neither urban nor rural, and subsequent calculations cannot be performed. It has been observed that most of these grids are scattered, so they are classified as rural areas. Since the shape of the Dongjiang River Basin does not completely fit the provincial and municipal administrative boundaries, for the very small number of boundary parts that are difficult to fit, they will be merged into the nearest province / city.

[0072] Step 3: Estimate the amount of domestic sewage and pollutants generated by the grid

[0073] Use Python to calculate the amount of domestic sewage and pollutants generated in the grid.

[0074] The estimation method for the amount of urban domestic sewage and pollutants (chemical oxygen demand, ammonia nitrogen, total nitrogen, total phosphorus) is as follows:

[0075] WW=POP×PUW p ×α p

[0076] POL i =WW×β p,i

[0077] Among them, WW is the sewage generation in the grid (liters / day); POP is the population in the grid (people); p represents the province, including Guangdong Province and Jiangxi Province; PUW p is the per capita comprehensive water consumption coefficient of the province p where the grid is located (liter / (person / day)); p is the pollution reduction coefficient of the province p where the grid is located; i represents the type of pollutant, including chemical oxygen demand, ammonia nitrogen, total nitrogen, and total phosphorus; POL i is the amount of pollutant i produced in the grid (mg / day); β p,i is the pollution concentration coefficient of pollutant i (mg / L). p , β p,i From the Manual of Pollution Emission Coefficients from Domestic Sources.

[0078] The amount of rural domestic sewage and pollutants (chemical oxygen demand, ammonia nitrogen, total nitrogen, total phosphorus) generated is estimated as follows:

[0079] WW=POP×γ c

[0080] POL i =WW×β c,i

[0081] Among them, c represents a city, including Guangzhou, Shaoguan, Shenzhen, Huizhou, Meizhou, Heyuan, Dongguan, and Ganzhou; γ cis the rural sewage discharge coefficient of city c where the grid is located (liter / (person / day)); c,i is the pollution intensity coefficient of pollutant i in city c where the grid is located (g / (person·L)). c and β c,i From the Manual of Pollution Emission Coefficients from Domestic Sources.

[0082] Step 4: Visualization

[0083] Since the python output is only in .tif format, without other geographic information, it is not conducive to subsequent processing and data display. Therefore, the results are imported into ArcGIS. First, the projection consistent with the original data is defined, that is, WGS_1984_Web_Mercator_Auxiliary_Sphere, and then georeferencing and resampling are performed to make the position of the result consistent with the original data. Finally, the color expression of the numerical value is adjusted and elements such as legend, compass, scale, etc. are added to fully display the results.

[0084] In the output result, the resolution is 100m, the image pixels are 3329×2726, and the number of grids in the Dongjiang River Basin is 4216795. Figure 5 As shown in the figure, it is the sewage volume estimation result diagram, such as Figure 6 The figure shows the estimation results of four types of sewage pollutants, namely COD (chemical oxygen demand), ammonia nitrogen, TN (total nitrogen) and TP (total phosphorus). Figure 5 and Figure 6 It can be seen that the amount of sewage and pollutants generated are concentrated in the lower reaches of the Dongjiang River Basin and the center of the Dongjiang River Basin, especially in cities in the Pearl River Delta such as Guangzhou and Shenzhen, where the amount of domestic sewage and pollutants generated is large.

[0085] The same or similar reference numerals correspond to the same or similar components;

[0086] The terms used in the drawings to describe positional relationships are only used for illustrative purposes and should not be construed as limiting this patent;

[0087] Obviously, the above embodiments of the present invention are merely examples for clearly illustrating the present invention, and are not intended to limit the embodiments of the present invention. For those skilled in the art, other different forms of changes or modifications can be made based on the above description. It is not necessary and impossible to list all the embodiments here. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should be included in the protection scope of the claims of the present invention.

Claims

1. A high-resolution, large-scale method for estimating the amount of domestic sewage and domestic sewage pollutants produced, characterized in that: The following steps are involved: Obtain open source data on the target area; Dividing the target area into multiple independent calculation areas based on the spatial resolution of the open source data, and processing the open source data to obtain parameter values ​​of each independent calculation area; Based on the parameter values ​​of each independent calculation area and the preset calculation formula, the amount of domestic sewage generated and the amount of domestic sewage pollutants generated in each independent calculation area are calculated.

2. The high-resolution, large-scale method for estimating the amount of domestic sewage and domestic sewage pollutants according to claim 1 is characterized in that: The open source data of the target area include the land use type, population, provincial administrative region, municipal administrative region, pollution reduction coefficient, pollution concentration coefficient, and rural sewage discharge coefficient of the target area.

3. The high-resolution, large-scale method for estimating the amount of domestic sewage and domestic sewage pollutants according to claim 2 is characterized in that: The target area is divided into multiple independent calculation areas based on the spatial resolution of the open source data, and the open source data is processed to obtain parameter values ​​of each independent calculation area, including: The spatial resolution of the population in the open source data is used as the spatial resolution of each independent calculation area. The study area is divided into multiple square grids according to the spatial resolution of each independent calculation area, and each square grid is used as an independent calculation area. The open source data of each independent computing area is unified in coordinate system, resampled, cropped and outlier processed to realize the division of open source data into each independent computing area as the parameter value of each independent computing area.

4. The high-resolution, large-scale method for estimating the amount of domestic sewage and domestic sewage pollutants according to claim 3 is characterized in that: The parameter values ​​of each independent calculation area include population parameter values, space parameter values, and water consumption parameter values.

5. The high-resolution, large-scale method for estimating the amount of domestic sewage and domestic sewage pollutants according to claim 4 is characterized in that: The population parameter value includes the population size of the independent calculation area.

6. The high-resolution, large-scale method for estimating the amount of domestic sewage and domestic sewage pollutants according to claim 4 is characterized in that: The spatial parameter values ​​include the land use type, province and city of the independent calculation area.

7. The high-resolution, large-scale method for estimating the amount of domestic sewage and domestic sewage pollutants according to claim 4, characterized in that: The water consumption parameter value includes the pollution concentration coefficient, pollution reduction coefficient or rural sewage discharge coefficient of the independent calculation area.

8. The method for estimating the amount of domestic sewage and domestic sewage pollutants produced with high resolution and large scale according to claims 5-7 is characterized in that: Based on the parameter values ​​of each independent calculation area and the preset calculation formula, the amount of domestic sewage generated and the amount of domestic sewage pollutants generated in each independent calculation area are calculated including: When the land use type of the independent calculation area is urban, the amount of urban domestic sewage and the amount of urban domestic sewage pollutants generated are calculated based on the population, pollution reduction coefficient, pollution concentration coefficient and the preset calculation formula for the amount of urban domestic sewage and the amount of urban domestic sewage pollutants generated; When the land use type of the independent calculation area is rural, the rural domestic sewage generation and rural domestic sewage pollutant generation are calculated based on the population size, rural sewage discharge coefficient, pollution concentration coefficient and the preset rural domestic sewage generation and rural domestic sewage pollutant generation calculation formula.

9. The high-resolution, large-scale method for estimating the amount of domestic sewage and domestic sewage pollutants according to claim 8, characterized in that: The calculation formulas for the amount of urban domestic sewage generated and the amount of pollutants generated by urban domestic sewage include: WW=POP×PUW p ×α p POL i =WW×β p,i Wherein, WW is the sewage generation in the independent calculation area (liters / day); POP is the population in the independent calculation area (people); p represents the province; PUW p is the per capita comprehensive water consumption coefficient of the province where the independent calculation area is located (liter / (person / day)); α p is the pollution reduction coefficient of the province where the independent calculation area is located; i represents the type of pollutant, including chemical oxygen demand, ammonia nitrogen, total nitrogen, and total phosphorus; POL i is the amount of pollutant i produced in the independent calculation area (mg / day); β p,i It is the pollution concentration coefficient (mg / L) of pollutant i in the province where the independent calculation area is located.

10. The high-resolution, large-scale method for estimating the amount of domestic sewage and domestic sewage pollutants according to claim 8, characterized in that: The calculation formulas for the amount of rural domestic sewage generated and the amount of pollutants generated by rural domestic sewage include: WW=POP×γ c POL i =WW×β c,i Where c represents a city; γ c β is the rural sewage discharge coefficient of the city where the independent calculation area is located (liter / (person / day)); c,i It is the pollution concentration coefficient (g / (person·L)) of pollutant i in the city where the independent calculation area is located.

Images