Pollution load river entering quantity calculation method

Abstract

The invention provides a pollution load in-river quantity calculation method, which comprises the following steps of: carrying out unit decomposition on a research range, distinguishing and classifying pollution loads entering a river, determining the contribution degree of each pollution load type to river water body pollution, carrying out statistics on different types of pollution load quantities in each unit, and calculating the pollution load quantity of each unit. The corresponding rivers are converged according to the analysis units, the pollution loads entering the rivers are summarized respectively, the pollution load river entering amount of the rivers in the research range is obtained, the pollution loads are classified scientifically and logically clearly, pollution types are positioned accurately and are independent of one another, the pollution loads cover comprehensively, the statistical method is proper, the calculation result is good in reference, and the calculation efficiency is high. Technical support and scientific basis can be provided for water quality simulation, river water environment management and the like of each river in a research area.

Description

Technical Field

[0001] The present invention belongs to the technical field of hydrological water resources management and supervision, and particularly relates to a method for calculating the amount of pollution load entering a river. Background Art

[0002] In the research of environmental impact assessment, demonstration of sewage outfall setting, and issues related to water environmental pollution, calculating the pollution load entering the target water body can not only master the current situation of water environmental pollution in the research area, but also provide data support for water resources protection and water pollution prevention and control work. At present, the amount of pollution load entering the river mostly focuses on the research of non-point source pollution entering the river. For example, in the Chinese patent "A method, system, medium and electronic device for regional non-point source pollution allocation and accounting" with the application number CN202410435750.5, the accounting method in this patent only considers the calculation of the amount of pollution load entering the river from non-point source pollution, without comprehensively considering other pollution sources in the river water body, and does not include point source pollution such as industrial enterprises and endogenous pollution. The statistical results are not scientific and comprehensive enough, and it is easy to occur the situation of logical confusion between the data of the amount of pollution load entering the river in the investigation and analysis method and the water quality status monitoring data and the pollution carrying capacity data. Therefore, there is an urgent need for a method for calculating the amount of pollution load entering the river that comprehensively considers various pollution emission sources, accurately divides the pollution load entering the river, classifies different pollution types, and conducts accurate quantitative statistics. Summary of the Invention

[0003] In view of this, the purpose of the present invention is to provide a method for calculating the amount of pollution load entering the river, which systematically considers various pollution sources and types in the analysis unit, and has the advantages of being scientific, comprehensive, and accurately positioned.

[0004] A method for calculating the amount of pollution load entering the river disclosed by the present invention includes the following steps:

[0005] Determine the research scope, which includes at least one river within the research scope. Divide the research scope into several analysis units, and the pollution load in the analysis unit includes at least one of point source pollution load, endogenous pollution load, non-point source pollution load, and runoff pollution load;

[0006] The point source pollution load includes at least one of the effluent load of the sewage treatment plant, domestic sewage load, and industrial enterprise sewage load;

[0007] The endogenous pollution load includes at least one of the river bottom sediment release load and the lake bottom sediment release load;

[0008] The non-point source pollution load includes at least one of the landfill leachate load, livestock and poultry breeding load, aquaculture load, urban runoff load, and agricultural non-point source load;

[0009] According to the rivers corresponding to the confluence of each analysis unit, the pollution loads entering each river are respectively summarized and counted to obtain the amount of pollution load entering the river for each river within the research scope.

[0010] Further, the formula for the effluent load of the sewage treatment plant in the point source pollution load is: WS ij = 3.65 * AP i * E j ,

[0011] In the formula: WS ij is the effluent load of the jth pollutant factor of the ith sewage treatment plant, AP i is the actual treatment capacity of the ith sewage treatment plant, and E j is the effluent concentration of the jth pollutant factor.

[0012] Further, the formula for the domestic sewage load in the point source pollution load is: DS ij = 3.65 * S i *(1 - C) * E j ,

[0013] In the formula: DS ij is the direct discharge sewage load of the jth pollutant factor of the ith analysis unit, S i is the sewage volume considering groundwater infiltration of the ith analysis unit, C is the sewage collection rate; E j is the effluent concentration of the jth pollutant factor.

[0014] Further, the formula for the industrial enterprise sewage load in the point source pollution load is:

[0015] In the formula, Qout i is the sewage discharge of the ith enterprise, COC ij is the concentration of the jth pollutant factor in the sewage discharge of the ith enterprise, is the sewage into-river coefficient of the ith enterprise, and Le ij is the into-river discharge of the jth pollutant factor of the ith enterprise.

[0016] Further, the formula for the sediment release load at the river bottom or lake bottom in the endogenous pollution load is: IL ij = 3.65 * 10 -4 * S i * B i * R ij ,

[0017] In the formula, IL ij is the sediment release load at the river bottom or lake bottom in the endogenous pollution load of the jth pollutant factor of the ith analysis unit, and Si is the length of the river or lake for the endogenous pollution load of the i-th analysis unit, B i is the average river width or lake width for the endogenous pollution load of the i-th analysis unit, R ij is the release rate of the j-th pollutant factor in the i-th analysis unit.

[0018] Furthermore, the calculation formula for the landfill leachate load in the non-point source pollution load is:

[0019] In the formula, C ij is the pollutant emission amount of the j-th pollutant factor in the i-th analysis unit, A i is the landfill area of the i-th analysis unit, P i is the rainfall infiltration runoff depth of the i-th analysis unit, is the runoff coefficient of the i-th analysis unit, Lp i is the landfill leachate load amount of the i-th analysis unit, E ij is the leachate concentration of the j-th pollutant factor in the i-th analysis unit.

[0020] Furthermore, the calculation formula for the aquaculture pollution load amount in the non-point source pollution load is:

[0021] In the formula: LY ij is the generation amount of the j-th pollutant factor in the i-th analysis unit, is the aquaculture area of the i-th analysis unit, A i y is the average aquaculture amount per mu of the i-th analysis unit, δ ij is the pollution discharge coefficient of the j-th pollutant factor in the i-th analysis unit.

[0022] Furthermore, the calculation formula for the livestock and poultry breeding load amount in the non-point source pollution load is:

[0023] In the formula: LS i is the generation amount of the j-th pollutant factor in the i-th analysis unit, are the numbers of poultry, pigs, cows, and sheep in the i-th analysis unit respectively, are the unit pollutant generation amounts of the j-th pollutant factor in the pollutants generated by poultry, pigs, cows, and sheep respectively, and K is the coefficient of entry into the river.

[0024] Furthermore, the calculation formula for the urban runoff load amount in the non-point source pollution load is:

[0025] Where: Lr ij is the river input load of the j-th pollutant factor in the i-th analysis unit, St i is the town area of the i-th analysis unit, P i is the annual rainfall of the i-th analysis unit, α i is the comprehensive runoff coefficient of the i-th analysis unit; is the average concentration of the j-th pollutant factor in the i-th analysis unit during a single rainfall, δ i is the river input coefficient of the i-th analysis unit.

[0026] Furthermore, the pollutant factors include at least one of COD, NH3-N, and TP.

[0027] The present invention has the following beneficial effects:

[0028] The method for calculating the river input amount of pollution load in the present invention decomposes the research scope into units, classifies and differentiates the pollution loads entering the river, selects various types of collectable data or existing reliable data, establishes relevant formulas between each pollution load and the above data, quantitatively calculates and classifies and statistically analyzes the different types of pollution load amounts in each unit, determines the contribution degree of each pollution load type to the river water body pollution, and according to the rivers corresponding to the confluence of each analysis unit, respectively summarizes each pollution load entering each river to obtain the river input amount of pollution load of each river within the research scope. Each pollution load is scientifically classified, logically clear, each pollution type is accurately positioned and independent of each other, the pollution load is comprehensively covered, the statistical method is appropriate, the calculation result has good reference value, and it can provide technical support and scientific basis for water quality simulation of rivers and river water environment management in the research area. Description of the Drawings

[0029] Figure 1 is a schematic flow chart of the method for calculating the river input amount of pollution load provided by some embodiments of the present invention,

[0030] Figure 2 is a schematic diagram of the research area and the division of analysis units provided by an embodiment of the present invention,

[0031] Figure 3 is a schematic flow chart of the load calculation of industrial enterprises provided by an embodiment of the present invention. Detailed Embodiments

[0032] In order to more clearly and completely describe the technical solution of the present invention, the present invention will be further described in detail below through specific embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention, and various changes can be made within the scope defined by the rights of the present invention.

[0033] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as those commonly understood by those skilled in the art of the present invention. The terms used in the specification of the present invention herein are only for the purpose of describing specific embodiments and are not intended to limit the present invention.

[0034] The present invention divides the pollution sources that cause the deterioration of river water environment into point source pollution load, internal source pollution load and non-point source pollution load.

[0035] Point source pollution loads are caused by the rapid development of society and economy. Many cities are facing development and construction, which has led to insufficient urban industrial and domestic sewage collection capacity, resulting in industrial and domestic sewage being discharged directly into the river without entering the sewage interception network; some cities have insufficient sewage treatment capacity, and there are problems with direct sewage discharge or overloaded operation of sewage treatment facilities, resulting in substandard tailwater entering the river; the combined rainwater and sewage pipe network in old urban areas causes overflow pollution into the river during the rainy season, mixed and wrong connections of rainwater and sewage pipes cause rainwater outlets to discharge sewage into the river during dry days, and the initial rainwater in the diversion rainwater pipe is not effectively controlled and discharged into the river. Point sources have many pollutants and complex components, and their changing patterns are affected by the discharge patterns of industrial wastewater and domestic sewage.

[0036] Endogenous pollution load refers to the organic matter and nitrogen and phosphorus pollutants in polluted water bodies that are deposited in the sediment through precipitation or adsorption of particulate matter. Under appropriate conditions, a large amount of pollutants are released from the sediment into the water body.

[0037] Non-point source pollution load refers to the pollutants in dissolved or particulate form that flow from non-specific locations into receiving water bodies such as rivers through the runoff process after being washed by precipitation. The sources of non-point source pollution are relatively complex, including: the random accumulation of domestic and construction waste on both sides of the river and the leachate of garbage discharged into the river with rainwater, causing pollution; the rainwater runoff pollution caused by the erosion of the surface, soil and road surfaces along the river by rainfall; the loss of fertilizers from agricultural cultivation on both sides of the river and the discharge of wastewater from livestock, poultry and aquaculture by individual households, which causes a large amount of nitrogen, phosphorus and other pollutants to enter the water body with rainwater, causing pollution, etc. Non-point source pollution is seasonal and random.

[0038] The above-mentioned point source pollution load, endogenous pollution load and non-point source pollution load are measured using pollutant factors. Pollutant factors refer to any substance that enters and acts on the river water environment system in inappropriate concentration, quantity, speed, form and pathway, and pollutes or destroys the river water environment ecosystem, including but not limited to COD (Chemical Oxygen Demand), NH3-N (ammonia nitrogen content in water) and TP (Total Phosphorus, the total amount of phosphorus in the water sample).

[0039] like Figure 1As shown in the figure, a method for calculating the amount of pollution load entering the river provided by the present invention specifically includes the following steps:

[0040] Step 1: Determine the research scope and divide the analysis units. When determining the research scope and dividing the analysis units, the basis for dividing the analysis units includes hydrometeorological data, topographic and geomorphic information, and administrative division information. According to the above division basis, the analysis units are divided into different types. Among them, there are analysis units with clear water use characteristics and water supply sources, such as village-townships divided according to administrative regions, or agricultural and forestry lands divided according to usage, such as livestock plots and aquaculture lands. It can also include other water bodies entering the river within the research scope, such as river runoff, lakes, or reservoirs. According to water resource utilization and underlying surface conditions, the research area can also be divided into multiple analysis units, and each analysis unit has clear water use characteristics and water supply sources. The above is only an example and does not limit the analysis units, and other division types can also be included.

[0041] Step 2: Statistically calculate the water consumption data of different types in each analysis unit divided in Step 1, including industrial water use, urban domestic and commercial water use, and rural domestic and commercial water use.

[0042] Step 3: According to the water consumption data in Step 2, calculate the sewage volume of each analysis unit. The formula is as follows:

[0043] S ij =WC iJ *PC j

[0044] In the formula: S ij is the sewage volume of the jth water volume data of the ith analysis unit, unit: 10,000 m 3 / d; WC ij is the jth water volume data of the ith analysis unit; PC j is the pollution generation coefficient corresponding to the jth water volume data.

[0045] Step 4: Multiply the sewage volume of each analysis unit obtained in Step 3 by the corresponding groundwater infiltration coefficient of each analysis unit to obtain the sewage volume considering groundwater infiltration of each analysis unit.

[0046] Step 5: Read the sewage collection rate of each analysis unit, multiply it by the sewage volume considering groundwater infiltration of each analysis unit obtained in Step 4 to obtain the sewage treatment volume of each sewage treatment plant, with the unit of 10,000 m 3 / d; Read the concentration of each pollutant factor of the influent concentration and effluent concentration of each sewage treatment plant or sewage treatment station, with the unit of mg / L; Calculate the pollutant factor load of each sewage treatment plant or sewage treatment station. The formula is as follows:

[0047] WS ij =3.65*APi *E j

[0048] Where: WS ij is the effluent load of the jth pollutant factor of the ith sewage treatment plant / station, unit: t / a; AP i is the actual treatment capacity of the ith sewage treatment plant / station (considering the sewage volume infiltrated from groundwater * sewage collection rate), unit: 10,000 m 3 / d; E j is the effluent concentration of the jth pollutant factor, unit: mg / L.

[0049] Step 6, according to the sewage volume infiltrated from groundwater obtained in Step 4 for each analysis unit, calculate the industrial enterprise load and its river inflow. The industrial enterprise load algorithm mainly uses the enterprise outlet data collected to calculate the industrial sewage concentration of each analysis unit in the industrial enterprise load analysis.

[0050] Read the sewage discharge volume, pollutant concentration and river inflow coefficient of each enterprise; read the data of the analysis unit and the basin to which each enterprise belongs; calculate the pollutant discharge volume and river inflow of each enterprise, and the formula is as follows:

[0051] LO ij = Qout i * COC ij * 3.65

[0052]

[0053] Where: LO ij is the generation volume of the jth pollutant factor of the ith enterprise, unit: t / a; Qout i is the sewage discharge volume of the ith enterprise, unit: m 3 / d; COC ij is the concentration of the jth pollutant factor in the sewage discharge volume of the ith enterprise, unit: mg / L; is the river inflow coefficient of the sewage of the ith enterprise; Le ij is the river inflow discharge volume of the jth pollutant factor of the ith enterprise, unit: t / a.

[0054] Perform interpolation calculation on the industrial water use of each analysis unit. The interpolation calculation method can adopt the Chinese patent "A Multi-attribute Data Interpolation Method" with the application number "202011260411.6". According to the enterprise-analysis unit matching relationship, enterprise outlet-analysis unit matching relationship and the industrial water use of each analysis unit, count the industrial sewage concentration, total generation load and total direct discharge load of each analysis unit. The schematic diagram of the industrial enterprise load calculation process is as Figure 2 shown.

[0055] Step 7: Based on the sewage volume considering groundwater infiltration obtained in Step 4, calculate the domestic sewage load and its river discharge. The calculation of domestic sewage load includes two statistical methods: the survey and statistics method and the population quota method. If relevant information of domestic sewage outlets is collected, the survey and statistics method is adopted; if population data is collected, the population quota method is adopted. The survey and statistics method reads the information of domestic sewage outlets, including sewage discharge type, town of origin, river of origin, sewage discharge volume, concentration of each pollutant factor in the discharged sewage, etc.; among them, for the survey and statistics method - urban areas, read the urban domestic and commercial water consumption (unit: 10,000 m 3 / d), urban sewage collection rate information, concentration data of each pollutant in domestic sewage (unit: mg / L), and calculate the total load generated in urban areas and the total river discharge load according to the data information, with the unit of both being t / a; for the survey and statistics method - rural areas, read the rural domestic and commercial water consumption (unit: 10,000 m 3 / d), rural sewage collection rate information, concentration data of each pollutant in domestic sewage (unit: mg / L), and calculate the total load generated in rural areas and the total river discharge load according to the data information, with the unit of both being t / a; the population quota method reads the urban population and rural population information of each analysis unit, calculates the urban sewage generation amount of each analysis unit according to the urban population and the urban per capita sewage generation quota, calculates the rural sewage generation amount of each analysis unit according to the rural population and the rural per capita sewage generation quota, and calculates the river discharge amount of urban and rural domestic sewage of each analysis unit according to the river discharge coefficient respectively;

[0056] DS ij = 3.65 * S i *(1 - C) * E j

[0057] In the formula: DS ij is the direct discharge sewage load of the jth pollutant factor in the ith analysis unit, unit: t / a; S i is the sewage volume considering groundwater infiltration in the ith analysis unit, unit: 10,000 m 3 / d; C is the sewage collection rate; E j is the effluent concentration of the jth pollutant factor, unit: mg / L.

[0058] Step 8: Calculate the endogenous load. The endogenous load is mainly the release load from river / lake bottom mud. Read the endogenous load river length, average river width of the endogenous load, and the pollutant release rate of each pollutant factor of all analysis units; statistically calculate the release load of each pollutant factor in the river / lake. The formula is as follows:

[0059] IL ij = 3.65 * 10 -4 * S i * B i * R ij

[0060] Where: $I_{L}$ ij is the endogenous load release load of the $j$-th pollutant factor in the $i$-th analysis unit, unit: t / a; $S$ i is the endogenous load river length of the $i$-th analysis unit, unit: km; $B$ i is the average river width of the endogenous load of the $i$-th analysis unit, unit: m; $R$ ij is the release rate of the $j$-th pollutant factor in the $i$-th analysis unit, unit: mg / (m 2 *d).

[0061] Step 9, landfill leachate load calculation. Under the long-term scouring of rainfall, the pollutants generated by the landfill may eventually flow into the river through groundwater. Read the data of each analysis unit, including the landfill area, rainfall infiltration runoff depth and runoff coefficient; calculate the landfill leachate generation amount in each analysis unit, and the formula is as follows:

[0062]

[0063] Where: $L_{p}$ i is the landfill leachate load of the $i$-th analysis unit, unit: m 3 / a; $A$ i is the landfill area of the $i$-th analysis unit, unit: m 2 ; $P$ i is the rainfall infiltration runoff depth of the $i$-th analysis unit, unit: mm; is the runoff coefficient of the $i$-th analysis unit. Read the pollutant concentration data of the leachate of each analysis unit; calculate the emission amount of each pollutant factor in each analysis unit, and the formula is as follows:

[0064]

[0065] Where: $C$ ij is the pollutant emission amount of the $j$-th pollutant factor in the $i$-th analysis unit, unit: t / a; $L_{p}$ i is the landfill leachate load of the $i$-th analysis unit, unit: m 3 / a; $E$ ij is the leachate concentration of the $j$-th pollutant factor in the $i$-th analysis unit, unit: mg / L.

[0066] Step 10, Calculation of livestock and poultry breeding load. The livestock and poultry breeding load is mainly used to calculate the load of each pollutant factor entering the river according to the quota table of livestock and poultry load and the number of livestock and poultry breeding in each analysis unit. Read the pollutant generation amount per unit of different livestock and poultry species. The livestock and poultry species can be divided into four types, including: poultry, pigs, cows, and sheep; obtain the coefficient of the livestock and poultry breeding load entering the river; read the number of different livestock and poultry species in each analysis unit; calculate the load generation amount of different pollutant factors. The formula is as follows:

[0067]

[0068] In the formula: LS i is the generation amount of the jth pollutant factor in the ith analysis unit, unit: t / a; are the numbers of poultry, pigs, cows, and sheep in the ith analysis unit respectively, unit: heads; are the unit pollutant generation amounts of the jth pollutant factor in the pollutants generated by poultry, pigs, cows, and sheep respectively, unit: kg / a; K is the coefficient entering the river.

[0069] Step 11, Calculation of aquaculture load. This algorithm mainly calculates the pollution load brought by aquaculture. Before calculation, it is necessary to obtain the aquaculture information of each analysis unit, including the aquaculture area, sewage discharge coefficient, etc. Read the aquaculture area and the average aquaculture amount per mu of each analysis unit; read the discharge coefficient of each pollutant factor of each analysis unit; calculate the pollutant generation amount. The calculation formula is as follows:

[0070]

[0071] LY ij =3.65*10 -4 *Y i *δ ij

[0072] In the formula: Y i is the aquaculture amount of the ith analysis unit, unit: kg; is the aquaculture area of the ith analysis unit, unit: mu; is the average aquaculture amount per mu of the ith analysis unit, unit: kg / mu; LY ij is the generation amount of the jth pollutant factor in the ith analysis unit, unit: t / a; δ ij is the sewage discharge coefficient of the jth pollutant factor in the ith analysis unit.

[0073] Step 12, Calculation of urban runoff load. This algorithm mainly calculates the load generated by rainfall runoff based on annual rainfall, comprehensive runoff coefficient, average concentration of pollutants in a single rainfall event, and river entry coefficient. Read the data of the town area, annual rainfall, and comprehensive runoff coefficient of each analysis unit; obtain the average concentration of pollutants in a single rainfall event (EMC) data of each pollutant factor in each analysis unit; calculate the discharge load, and the formula is as follows:

[0074]

[0075] In the formula: LC ij is the discharge load of the jth pollutant in the ith analysis unit, unit: t / a; St i is the town area of the ith analysis unit, unit: km 2 ; P i is the annual rainfall of the ith analysis unit, unit: mm; α i is the comprehensive runoff coefficient of the ith analysis unit; is the average concentration of pollutants in a single rainfall event of the jth pollutant factor in the ith analysis unit, unit: mg / L.

[0076] Read the river entry coefficient; calculate the river entry load, and the formula is as follows:

[0077] Lr ij = LC ij *δ i

[0078] In the formula: Lr ij is the river entry load of the jth pollutant in the ith analysis unit, unit: t / a; LC ij is the discharge load of the jth pollutant in the ith analysis unit, unit: t / a; δ i is the river entry coefficient of the ith analysis unit.

[0079] Step 13, Calculation of agricultural non-point source load. The algorithm for agricultural non-point source load is similar to that for livestock and poultry breeding load. The river entry amount of each analysis unit is calculated through the area of different agricultural land types and the loss amount of pollutants.

[0080] Obtain the data of the agricultural load quota table, mainly the loss amount of pollutants per mu of different agricultural land types, unit: kg / (mu·a). The agricultural land types include 6 categories: paddy fields, vegetable fields, dry land, orchards, tea gardens, and forest lands; read the area data of 6 types of agricultural land in each analysis unit, unit: mu; based on the above data, calculate the river entry load of each pollutant index, unit: t / a.

[0081] Step 14: Statistically analyze the pollution load calculation results of all the analysis units in Steps 5 to 13 above. According to the topological relationship of the river channels within the research scope, including the flow direction, confluence points, branching conditions, and connectivity of the rivers, determine which river each analysis unit flows into respectively, and incorporate the pollution load quantity into the corresponding river to complete the analysis and calculation of the pollution load entering the river.

[0082] Taking a certain basin as an example, the pollution load entering the river in this basin is calculated in combination with the attached drawings. As Figure 1 shown, based on the water supply data and population data of a certain basin, the calculation of the pollution load entering the river in this basin includes the following steps:

[0083] Step 1: Determine the research scope and divide the analysis units. The basis for dividing the analysis units includes hydrometeorological data, topographic and geomorphic information, and administrative division information; as Figure 2 shown, there are three rivers, namely River A, River B, and River C, and two sewage treatment plants, namely W1 and W2, within the basin scope. According to the water resource utilization and the underlying surface conditions, the research area is divided into 8 analysis units, and each analysis unit has clear water use characteristics and water supply sources.

[0084] Step 2: Statistically analyze the water consumption data of different types for each of the analysis units divided in Step 1, including industrial water use, urban domestic and commercial water use, and rural domestic and commercial water use. The results are shown in Table 1.

[0085] Table 1 Water Consumption Statistics (Unit: 10,000 m 3 / d)

[0086] Analysis unit Industrial water Urban domestic and commercial water Rural domestic and commercial water 1 8.634 1.09 1.02 2 13.61 19.87 5.32 3 219.127 4.95 0.468 4 57.704 1.25 1.16 5 112.03 4.52 0.89 6 25.32 10.11 2.15 7 72.739 1.39 1.3 8 404.444 7.62 0.95

[0087] Step 3: Multiply the various types of water consumption obtained in Step 2 by their corresponding pollution production coefficients to obtain the sewage volume of each analysis unit. The formula is as follows:

[0088] S ij = WC iJ * PC j (1-1)

[0089] In the formula: S ij is the sewage volume of the jth water volume data of the ith analysis unit, with the unit of 10,000 m 3 / d; WC ij is the jth water volume data of the ith analysis unit; PC j is the pollution production coefficient corresponding to the jth water volume data.

[0090] For convenience of calculation, the pollutant production coefficients of industrial water use for all 8 analysis units are set to 0.5, those of urban domestic and commercial water use are set to 0.85, and those of rural domestic and commercial water use are set to 0.5. The calculation results obtained by multiplying the water consumption in Table 1 by their corresponding pollutant production coefficients are shown in Table 2.

[0091] Table 2 Statistical Analysis Table of Sewage Volume (Unit: 10,000 m 3 / d)

[0092]

[0093]

[0094] Step 4: Multiply the sewage volume of each analysis unit obtained in Step 3 by the corresponding groundwater infiltration coefficient of each analysis unit to obtain the sewage volume considering groundwater infiltration for each analysis unit. The groundwater infiltration coefficient is taken as 1.1, and the calculation results are shown in Table 3.

[0095] Table 3 Statistical Analysis Table of Sewage Volume Considering Leakage (Unit: 10,000 m 3 / d)

[0096] Analysis unit Industrial water Urban domestic and commercial water Rural domestic and commercial water 1 4.7487 1.01915 0.561 2 7.4855 18.57845 2.926 3 120.51985 4.62825 0.2574 4 31.7372 1.16875 0.638 5 61.6165 4.2262 0.4895 6 13.926 9.45285 1.1825 7 40.00645 1.29965 0.715 8 222.4442 7.1247 0.5225

[0097] Step 5: Multiply the sewage collection rate of each analysis unit by the sewage volume considering groundwater infiltration of each analysis unit obtained in Step 4 to obtain the sewage treatment volume of each sewage treatment plant, with the unit of 10,000 m 3 / d; Read the concentration of each pollutant factor of the influent and effluent of each sewage treatment plant or sewage treatment station, with the unit of mg / L; Calculate the load of each pollutant factor of each sewage treatment plant or sewage treatment station, and the formula is as follows:

[0098] WS ij =3.65*AP i *E j (1 - 2)

[0099] In the formula: WS ij is the effluent load of the jth pollutant factor of the ith sewage treatment plant / station, with the unit of t / a; AP i is the actual treatment capacity of the ith sewage treatment plant / station (sewage volume considering groundwater infiltration * sewage collection rate), with the unit of 10,000 m 3 / d; E j is the effluent concentration of the jth pollutant factor, with the unit of mg / L.

[0100] The sewage collection rate of industrial sewage is taken as 1, that of urban sewage is taken as 0.65, and that of rural sewage is taken as 0.35. The water collection range of W1 is Analysis Units 2, 3, and 4, and the water collection range of W2 is Analysis Unit 8. The calculation results are shown in Table 4.

[0101] Table 4 Effluent Load of Sewage Treatment Plant (Unit: t / a)

[0102]

[0103] Step 6. The sewage volume except that entering the sewage treatment plant is the river inflow directly discharged into the river, which is divided into two categories: one is industrial enterprises and the other is domestic sewage. For the load calculation method of industrial enterprises, mainly using the enterprise outlet data collected, calculate the industrial sewage concentration of each analysis unit in the industrial enterprise load analysis. Read the sewage discharge volume, pollutant concentration and river inflow coefficient of each enterprise; read the data of the analysis unit and the basin to which each enterprise belongs; calculate the pollutant discharge volume and river inflow volume of each enterprise, and the formula is as follows:

[0104] LO i = Qout i * COC ij * 3.65 (1 - 3)

[0105]

[0106] In the formula: LO i is the generation amount of the jth pollutant factor of the ith enterprise, unit: t / a; Qout i is the sewage discharge volume of the ith enterprise, unit: m 3 / d; COC ij is the concentration of the jth pollutant factor in the sewage discharge volume of the ith enterprise, unit: mg / L; is the river inflow coefficient of the sewage of the ith enterprise; Le ij is the river inflow discharge volume of the jth pollutant factor of the ith enterprise, unit: t / a.

[0107] Interpolate the industrial water use in Table 3. The interpolation calculation method follows the Chinese patent "A Multi-Attribute Data Interpolation Method" with the application number 202011260411.6. According to the enterprise-analysis unit matching relationship, enterprise outlet-analysis unit matching relationship and the industrial water use of each analysis unit, count the industrial sewage concentration, total generation load and total direct discharge load of each analysis unit. The schematic diagram of the industrial enterprise load calculation process is as Figure 3 shown, and the calculation results are shown in Table 5.

[0108] Table 5 River Inflow Load of Industrial Enterprises

[0109]

[0110] Step 7. Calculate the domestic sewage load. Using the survey and statistics method, read the urban domestic and commercial water use (unit: 10,000 m 3 / d), the urban sewage collection rate is taken as 0.65, and the concentration data of each pollutant in domestic sewage (unit: mg / L). According to the data information, calculate the total load generated in the urban area and the total load entering the river, both in t / a; Survey and statistics method - rural area, read the domestic and commercial water consumption in step 4 (unit: 10,000 m 3 / d), the rural sewage collection rate is taken as 0.35, and the concentration data of each pollutant in domestic sewage (unit: mg / L). According to the data information, calculate the total load generated in the rural area and the total load entering the river, both in t / a;

[0111] DS ij = 3.65 * S i *(1 - C)*E j (1 - 5)

[0112] In the formula: DS ij is the direct discharge sewage load of the jth pollutant factor in the ith analysis unit, unit: t / a; S i is the sewage volume considering groundwater infiltration in the ith analysis unit, unit: 10,000 m 3 / d;, C is the sewage collection rate; E j is the effluent concentration of the jth pollutant factor, unit: mg / L.

[0113] Table 6 Domestic sewage load entering the river

[0114]

[0115]

[0116] Step 8, endogenous load calculation. The endogenous load is mainly the load released from the river bottom sediment. Read the endogenous load river length, endogenous load average river width, and the pollutant release rate of each pollutant factor of all analysis units; Statistically calculate the release load of each pollutant factor in the river, and the formula is as follows:

[0117] IL ij = 3.65 * 10 -4 *S i *B i *R ij (1 - 6)

[0118] In the formula: IL ij is the endogenous load release load of the jth pollutant factor in the ith analysis unit, unit: t / a; S i is the endogenous load river length of the ith analysis unit, unit: km; B i is the endogenous load average river width of the ith analysis unit, unit: m; R ij is the release rate of the jth pollutant factor in the ith analysis unit, unit: mg / (m2 *d).

[0119] Substitute into the above formula for calculation, and the calculation results are shown in Table 7 below.

[0120] Table 7 Endogenous Load (Unit: t / a)

[0121]

[0122] Step 9: Calculate the landfill leachate load. Under the long-term scouring of rainfall, the pollutants generated by the landfill may eventually flow into the river through groundwater. Read the data of each analysis unit, including the landfill area, rainfall infiltration runoff depth, and runoff coefficient; calculate the landfill leachate generation amount in each analysis unit, and the formula is as follows:

[0123]

[0124] In the formula: Lp i is the landfill leachate load of the i-th analysis unit, unit: m 3 / a; A i is the landfill area of the i-th analysis unit, unit: m 2 ; P i is the rainfall infiltration runoff depth of the i-th analysis unit, unit: mm; is the runoff coefficient of the i-th analysis unit. Read the pollutant concentration data of the leachate in each analysis unit; calculate the emission amount of each pollutant factor in each analysis unit, and the formula is as follows:

[0125]

[0126] In the formula: C ij is the pollutant emission amount of the j-th pollutant factor in the i-th analysis unit, unit: t / a; Lp i is the landfill leachate load of the i-th analysis unit, unit: m 3 / a; E ij is the leachate concentration of the j-th pollutant factor in the i-th analysis unit, unit: mg / L.

[0127] Substitute into formula (1-6) for calculation, and the calculation results are shown in Table 8 below.

[0128] Table 8 Analysis Table of Leachate Generation Amount (Unit: 10,000 m 3 / a)

[0129]

[0130] Substitute into formula (1-7) for calculation, and the calculation results are shown in Table 9 below.

[0131] Table 9 Landfill Leachate Load (Unit: t / a)

[0132]

[0133] Step 10: Calculation of livestock and poultry breeding load. The livestock and poultry breeding load is mainly obtained by solving the river input load of each pollutant factor according to the quota table of livestock and poultry load and the number of livestock and poultry breeding in each analysis unit. Read the unit pollutant generation amount of different livestock and poultry species, which can be divided into four types: poultry, pigs, cows, and sheep; obtain the river input coefficient of livestock and poultry breeding load; read the number of different livestock and poultry species in each analysis unit; calculate the load generation amount of different pollutant factors, and the formula is as follows:

[0134]

[0135] In the formula: LS i is the generation amount of the jth pollutant factor in the ith analysis unit, unit: t / a; are the numbers of poultry, pigs, cows, and sheep in the ith analysis unit respectively, unit: head; are the unit pollutant generation amounts of the jth pollutant factor in the pollutants generated by poultry, pigs, cows, and sheep respectively, unit: kg / a; K is the river input coefficient.

[0136] Table 10 Quota Table of Livestock and Poultry Load

[0137]

[0138] The river input coefficient of livestock and poultry breeding load is taken as 0.23. The calculated livestock and poultry breeding load is shown in the following table.

[0139] Table 11 Livestock and Poultry Breeding Load Table

[0140]

[0141] Step 11: Calculation of aquaculture load. This algorithm mainly calculates the pollution load brought by aquaculture. Before calculation, it is necessary to obtain the aquaculture information of each analysis unit, including aquaculture area, sewage discharge coefficient, etc. Read the aquaculture area and average aquaculture amount per mu of each analysis unit; read the emission coefficient of each pollutant factor of each analysis unit; calculate the pollutant generation amount, and the calculation formula is as follows:

[0142]

[0143] LY ij =3.65*10 -4 *Y i *δ ij (1 - 11)

[0144] In the formula: Yi is the aquaculture volume of the i-th analysis unit, unit: kg; is the aquaculture area of the i-th analysis unit, unit: mu; is the average aquaculture volume per mu of the i-th analysis unit, unit: kg / mu; LY ij is the generation amount of the j-th pollutant factor in the i-th analysis unit, unit: t / a; δ ij is the pollution discharge coefficient of the j-th pollutant factor in the i-th analysis unit.

[0145] Table 12 Aquaculture volume

[0146]

[0147] Table 13 Generation amount of aquaculture pollutants Unit: t / a

[0148]

[0149] Step 12, calculation of urban runoff load. This algorithm mainly calculates the load generated by rainfall runoff according to annual rainfall, comprehensive runoff coefficient, average concentration of pollutants in a single rainfall, and river entry coefficient. Read the data of the town area, annual rainfall, and comprehensive runoff coefficient of each analysis unit; obtain the average concentration of pollutants in a single rainfall EMC data of each pollutant factor in each analysis unit; calculate the discharge load, and the formula is as follows:

[0150]

[0151] In the formula: LC ij is the discharge load of the j-th pollutant in the i-th analysis unit, unit: t / a; St i is the town area of the i-th analysis unit, unit: km 2 ; P i is the annual rainfall of the i-th analysis unit, unit: mm; α i is the comprehensive runoff coefficient of the i-th analysis unit; is the average concentration of the j-th pollutant factor in a single rainfall of the i-th analysis unit, unit: mg / L.

[0152] Read the river entry coefficient; calculate the river entry load, and the formula is as follows:

[0153] Lr ij = LC ij *δ i (1-13)

[0154] In the formula: Lr ij is the river entry load of the j-th pollutant in the i-th analysis unit, unit: t / a; LC ijis the pollutant emission load of the jth pollutant in the ith analysis unit, unit: t / a; δ i is the river entry coefficient of the ith analysis unit.

[0155] Table 14 Pollutant Emission Load of Urban Rainfall Runoff (unit: t / a)

[0156]

[0157]

[0158] Table 15 Pollutant River Entry Load of Urban Rainfall Runoff (unit: t / a)

[0159]

[0160] Step 13, calculation of agricultural non-point source load. The algorithm for agricultural non-point source load is similar to that for livestock and poultry breeding load. The river entry amount of each analysis unit is calculated through the area of different agricultural land types and the pollutant loss amount.

[0161] Obtain the data of the agricultural load quota table, mainly the pollutant loss amount per mu of different agricultural land types, unit: kg / (mu·a). The agricultural land types include 6 categories: paddy fields, vegetable fields, dry land, orchards, tea gardens, and forest lands; read the area data of the 6 types of agricultural land in each analysis unit, unit: mu; calculate the river entry load of each pollutant index according to the above data, unit: t / a.

[0162] Table 16 Agricultural Load Quota Table

[0163]

[0164]

[0165] Table 17 Agricultural Non-point Source Load Table

[0166]

[0167] Step 14, for the pollution load calculation results of all analysis units in the above steps 5 to 13, respectively judge which river each analysis unit enters. The corresponding relationship is shown in Table 18. Record the pollution load into the corresponding river to complete the analysis and calculation of the pollution load river entry amount. The final calculation results are shown in Tables 19 - 21.

[0168] Table 18 Corresponding Relationship Table of River, Sewage Treatment Plant, and Analysis Unit

[0169]

[0170] Table 19 Total Load Statistics Table of River A

[0171]

[0172]

[0173] Table 20 Statistical Table of Total Load of River B

[0174]

[0175] Table 21 Statistical Table of Total Load of River C

[0176]

[0177] The above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that they can still modify the technical solutions recorded in the foregoing embodiments, or perform equivalent replacements on some of the technical features; and these modifications or replacements all fall within the protection scope of the present invention. Therefore, the protection scope of this invention patent shall be subject to the appended claims.

Claims

1. A method for calculating the amount of pollution load entering a river, characterized in that, It includes the following steps: Determine the research scope, which includes at least one river. Divide the research scope into several analysis units, and the pollution load in each analysis unit includes at least one of point source pollution load, endogenous pollution load, non-point source pollution load, and runoff pollution load; The point source pollution load includes at least one of the effluent load of sewage treatment plants, domestic sewage load, and industrial enterprise sewage load; The endogenous pollution load includes at least one of the sediment release load from the river bottom and the sediment release load from the lake bottom; The non-point source pollution load includes at least one of the landfill leachate load, livestock and poultry breeding load, aquaculture load, urban runoff load, and agricultural non-point source load; According to the rivers corresponding to the confluence of each analysis unit, respectively summarize and statistically calculate the pollution loads entering each river to obtain the pollution load inflow into each river within the research scope.

2. The pollution load inflow calculation method according to claim 1, characterized in that The calculation formula for the effluent load of the sewage treatment plant in the point source pollution load is: WS ij = 3.65 * AP i * E j , Where: WS ij is the effluent load of the j-th pollutant factor of the i-th sewage treatment plant, AP i is the actual treatment capacity of the i-th sewage treatment plant, E j is the effluent concentration of the j-th pollutant factor.

3. The method for calculating the amount of pollution load entering the river according to claim 1, characterized in that The calculation formula for the domestic sewage load in the point source pollution load is: DS ij = 3.65 * S i *(1 - C) * E j , Where: DS ij is the direct discharge sewage load of the j-th pollutant factor in the i-th analysis unit, S i is the sewage volume considering groundwater infiltration in the i-th analysis unit, C is the sewage collection rate; E j is the effluent concentration of the j-th pollutant factor.

4. The method for calculating the amount of pollution load entering the river according to claim 1, characterized in that The calculation formula for the sewage load of industrial enterprises in the point source pollution load is as follows: Where Qout i is the sewage discharge of the i-th enterprise, and COC ij is the concentration of the j-th pollutant factor in the sewage discharge of the i-th enterprise, is the sewage into-river coefficient of the i-th enterprise, and Le ij is the into-river discharge of the j-th pollutant factor of the i-th enterprise.

5. The method for calculating the amount of pollution load entering the river according to claim 1, characterized in that The calculation formula for the release load of sediment at the bottom of a river or lake in the endogenous pollution load is: IL ij = 3.65 * 10 -4 * S i * B i * R ij , where, IL ij is the release load of the sediment at the bottom of the river or lake in the endogenous pollution load of the j-th pollutant factor in the i-th analysis unit, S i is the length of the river or lake of the endogenous pollution load in the i-th analysis unit, B i is the average river width or lake width of the endogenous pollution load in the i-th analysis unit, R ij is the release rate of the j-th pollutant factor in the i-th analysis unit.

6. The method for calculating the amount of pollution load entering the river according to claim 1, characterized in that, The calculation formula for the landfill leachate load in the non-point source pollution load is as follows: Where C ij is the pollutant emission of the j-th pollutant factor in the i-th analysis unit, A i is the landfill area of the i-th analysis unit, P i is the rainfall infiltration runoff depth of the i-th analysis unit, is the runoff coefficient of the i-th analysis unit, Lp i is the leachate load of the i-th analysis unit, E ij is the leachate concentration of the j-th pollutant factor in the i-th analysis unit.

7. The method for calculating the amount of pollution load entering the river according to claim 1, characterized in that The calculation formula for the aquaculture pollution load in the non-point source pollution load is as follows: Where: LY ij is the generation amount of the jth pollutant factor in the ith analysis unit, is the aquaculture area of the ith analysis unit, is the average aquaculture amount per mu of the ith analysis unit, δ ij is the pollutant discharge coefficient of the jth pollutant factor in the ith analysis unit.

8. The method for calculating the amount of pollution load entering the river according to claim 1, characterized in that, The calculation formula for the livestock and poultry breeding load in the non-point source pollution load is as follows: Where: LS i is the generation amount of the jth pollutant factor in the ith analysis unit, are respectively the numbers of poultry, pigs, cattle, and sheep in the ith analysis unit, are respectively the unit pollutant generation amounts of the jth pollutant factor in the pollutants generated by poultry, pigs, cattle, and sheep, and K is the river entry coefficient.

9. The method for calculating the amount of pollution load entering the river according to claim 1, characterized in that, The formula for calculating the urban runoff load in the non-point source pollution load is as follows: Where: Lr ij is the river input load of the jth pollutant factor in the ith analysis unit, St i is the town area of the ith analysis unit, P i is the annual rainfall of the ith analysis unit, α i is the comprehensive runoff coefficient of the ith analysis unit; C j i is the average concentration of the jth pollutant factor in the ith analysis unit during a single rainfall, δ i is the river input coefficient of the ith analysis unit.

10. The method for calculating the amount of pollution load entering the river according to any one of claims 2-9, characterized in that, The pollutant factors include at least one of COD, NH3-N, and TP.

Images