CN110728035B - Pollutant total amount control method based on control of section water quality reaching standard - Google Patents
Pollutant total amount control method based on control of section water quality reaching standard Download PDFInfo
- Publication number
- CN110728035B CN110728035B CN201910903473.5A CN201910903473A CN110728035B CN 110728035 B CN110728035 B CN 110728035B CN 201910903473 A CN201910903473 A CN 201910903473A CN 110728035 B CN110728035 B CN 110728035B
- Authority
- CN
- China
- Prior art keywords
- water quality
- pollution
- river
- amount
- water
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 155
- 239000003344 environmental pollutant Substances 0.000 title claims abstract description 66
- 231100000719 pollutant Toxicity 0.000 title claims abstract description 66
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000011160 research Methods 0.000 claims abstract description 43
- 238000004364 calculation method Methods 0.000 claims abstract description 31
- 238000013178 mathematical model Methods 0.000 claims abstract description 28
- 238000012544 monitoring process Methods 0.000 claims abstract description 20
- 238000011835 investigation Methods 0.000 claims abstract description 17
- 230000001360 synchronised effect Effects 0.000 claims abstract description 6
- 239000010865 sewage Substances 0.000 claims description 52
- 230000015556 catabolic process Effects 0.000 claims description 12
- 238000006731 degradation reaction Methods 0.000 claims description 12
- 238000004088 simulation Methods 0.000 claims description 12
- 238000011144 upstream manufacturing Methods 0.000 claims description 12
- 244000144977 poultry Species 0.000 claims description 10
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 6
- 238000009360 aquaculture Methods 0.000 claims description 6
- 244000144974 aquaculture Species 0.000 claims description 6
- 238000004519 manufacturing process Methods 0.000 claims description 5
- 238000009395 breeding Methods 0.000 claims description 4
- 230000001488 breeding effect Effects 0.000 claims description 4
- 238000007599 discharging Methods 0.000 claims description 4
- 244000144972 livestock Species 0.000 claims description 4
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims description 3
- 239000003905 agrochemical Substances 0.000 claims description 3
- XKMRRTOUMJRJIA-UHFFFAOYSA-N ammonia nh3 Chemical compound N.N XKMRRTOUMJRJIA-UHFFFAOYSA-N 0.000 claims description 3
- 229910052757 nitrogen Inorganic materials 0.000 claims description 3
- 229910052698 phosphorus Inorganic materials 0.000 claims description 3
- 239000011574 phosphorus Substances 0.000 claims description 3
- 238000003908 quality control method Methods 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 3
- 239000002351 wastewater Substances 0.000 claims description 3
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 2
- 229910052760 oxygen Inorganic materials 0.000 claims description 2
- 239000001301 oxygen Substances 0.000 claims description 2
- -1 permanganate index Chemical compound 0.000 claims 1
- 230000006872 improvement Effects 0.000 abstract description 12
- 238000004458 analytical method Methods 0.000 abstract description 10
- 230000009467 reduction Effects 0.000 abstract description 8
- 238000013480 data collection Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 8
- 230000008859 change Effects 0.000 description 7
- 238000004590 computer program Methods 0.000 description 7
- 238000009792 diffusion process Methods 0.000 description 6
- 238000010276 construction Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 238000013508 migration Methods 0.000 description 4
- 230000005012 migration Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 229910000069 nitrogen hydride Inorganic materials 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 239000013256 coordination polymer Substances 0.000 description 1
- 230000000994 depressogenic effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000010842 industrial wastewater Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000013441 quality evaluation Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000005067 remediation Methods 0.000 description 1
- 230000001932 seasonal effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A20/00—Water conservation; Efficient water supply; Efficient water use
- Y02A20/152—Water filtration
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Educational Administration (AREA)
- Economics (AREA)
- Development Economics (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Feedback Control In General (AREA)
- Sewage (AREA)
Abstract
The invention discloses a method for controlling the total amount of pollutants based on controlling the water quality of a section to reach the standard, which comprises data collection, on-site investigation and monitoring; carrying out investigation and analysis on the current situation of a pollution source and the current situation of water environment quality in a research area, and identifying a main pollution source and a main overproof factor; constructing a nested river network hydrodynamic mathematical model and a nested water quality mathematical model; calibrating according to the synchronous monitoring data of water quantity and water quality; calculating allowable discharge according to the pollution source investigation result and the model calibration result; and a water quality improvement scheme is provided. According to the invention, by dividing the research area range and the calculation unit, a nested area hydrodynamic water quality mathematical model is established, on the basis, the allowable discharge calculation research of the area is developed, a water quality improvement scheme is provided, and a technical support and a theoretical reference basis are provided for the pollutant reduction and environment improvement scheme.
Description
Technical Field
The invention relates to the field of water environment numerical simulation, in particular to a pollutant total amount control method based on section water quality control reaching the standard.
Background
Along with the development of economy and the increase of population, the water consumption is increased day by day, the shortage of water resource and the water environment pollution gradually become the bottleneck restricting the development of economy, the discharge amount of domestic sewage and industrial wastewater is increased continuously, the problem that a series of urban domestic sewage is directly discharged due to the lag of the construction of sewage collecting pipe networks in partial urban areas and the like is caused, and in addition, the influence of rural area source pollution and industrial point source pollution causes the water quality in partial areas to be in a stage and seasonal standard exceeding state, which needs to draw high attention to us, the calculation and research of the allowable discharge amount of regional pollutants needs to be carried out to provide technical support and theoretical reference basis for pollutant reduction and environment improvement scheme establishment, and a harmonious, healthy and perfect ecological water system combining manpower and nature is created in an effort.
The existing river treatment technology has the defects that simple pollution control is the traditional river treatment idea in China, the method is subjective, and the allowable discharge amount of pollutants is not quantified based on the actual conditions of migration, diffusion and water quality change of the pollutants in a water body, so that the aim of improving the water quality cannot be achieved, and the recovery of a water body ecosystem is more difficult.
Disclosure of Invention
The purpose of the invention is as follows: in order to solve the defects that the prior art only has simple pollution control, is subjective and does not quantify the allowable pollutant discharge amount based on the actual conditions of migration, diffusion and water quality change of pollutants in a water body, the invention provides a method for controlling the total quantity of pollutants on the basis of controlling the section water quality to reach the standard.
In order to achieve the purpose, the invention adopts the following technical scheme:
a method for controlling the total amount of pollutants reaching the water quality of a section comprises the following steps:
generalizing river channels in the selected investigation and research range and dividing a calculation unit of river channel sections, and determining a pollution source needing total amount control based on the pollution source in the calculation unit; constructing a river network hydrodynamic mathematical model and a water quality mathematical model;
determining boundary conditions and initial conditions required by a river network hydrodynamic mathematical model and a water quality mathematical model according to the collected and sorted hydrological water quality time data of the boundary of the research area and the upper monitoring point;
carrying out model parameter calibration on the river network hydrodynamic mathematical model and the water quality mathematical model based on water quantity and water quality synchronous monitoring data to obtain a water quality degradation coefficient;
and establishing a response relation between the water quality of the control section and the discharge capacity of the sewage discharge outlet according to the determined main pollution source, the pollution overproof factor and the water quality degradation coefficient obtained by the model calibration through the established mathematical model of the river network hydrodynamic force and the mathematical model of the water quality of the research area, and obtaining the allowable discharge capacity of pollutants from each sewage discharge outlet to the control section through a calculation method based on the standard-reaching allowable discharge capacity of the control section.
Compared with the prior art, the invention has the following beneficial technical effects:
according to the invention, by dividing the research area range and the calculation unit, a nested area hydrodynamic water quality mathematical model is established, on the basis, the allowable emission calculation research of the area is developed, and technical support and theoretical reference basis are provided for pollutant reduction and environmental remediation schemes. The research of the allowable emission amount has operability and foresight, and can provide reliable reference for further planning work. The numerical model simulation can analyze the water quantity and water quality change, can very conveniently show the migration and diffusion process of pollutants in the water body, can quickly predict the water quality change under specific hydrological conditions and specific working conditions, and can calculate important reference indexes such as water environment capacity, maximum pollutant discharge and the like according to the pollutant degradation coefficient determined by the model so as to provide a scientific and reasonable research area water environment improvement scheme.
Drawings
FIG. 1 is a flow chart of a method for controlling total amount of pollutants based on cross-section water quality control;
FIG. 2 is a comparison graph of a calculated value and an actual value of a point location model of a large watershed according to an embodiment of the present invention;
FIG. 3 is a comparison graph of the calculated value and the measured value of the point model in the study area according to the embodiment of the present invention.
Detailed Description
The present invention is further illustrated by the following examples, which are intended to be purely exemplary and are not intended to limit the scope of the invention, as various equivalent modifications of the invention will occur to those skilled in the art upon reading the present disclosure and fall within the scope of the appended claims.
As shown in fig. 1, a method for controlling the total amount of pollutants based on the control of the water quality of a section to reach the standard comprises the following steps:
(1) collecting data, surveying and monitoring on site;
(2) carrying out investigation and analysis on the current situation of a pollution source and the current situation of water environment quality in a research area, and identifying a main pollution source and a main overproof factor;
(3) constructing a nested river network hydrodynamic mathematical model and a nested water quality mathematical model;
(4) calibrating according to the water quantity and water quality synchronous monitoring data to obtain a water quality degradation coefficient;
(5) calculating allowable discharge according to the pollution source investigation result and the model calibration result;
the data collection in the step (1) generally comprises water level, water quantity, meteorological data and DEM of a research area, and the on-site investigation and monitoring are mainly the distribution positions of a gate, a pump and a sewage outlet.
The investigation and analysis of the current situation of the pollution source and the current situation of the water environment quality in the research area in the step (2) are mainly to determine a research investigation range by combining river basin DEM information and an actual catchment range, and further to perform refined calculation unit division on the research range, namely to divide the research investigation range based on factors such as administrative boundaries, key control nodes (such as important hydrological stations, gate dams, important branch inflow rivers, conventional monitoring sections, important pollution source sewage outlets), water function divisions and the like of the research investigation range, so as to facilitate the calculation of subsequent pollution sources. The calculated pollution sources comprise a living pollution source, an industrial pollution source and an agricultural non-point source, a required current situation investigation result of the pollution sources is obtained, and a main pollution source is identified; the current situation analysis and investigation of the water environment quality is based on the data of the section water quality and aims to identify main overproof factors. The water quality evaluation indexes comprise COD, ammonia nitrogen, total phosphorus, total nitrogen and the like, each control section has a corresponding water quality standard, and the exceeding indexes are regarded as exceeding factors.
The method for determining the main source of pollution is as follows: calculating the river inflow amount of each pollution source, and comparing the river inflow amount of each pollution source, wherein the highest pollution source is the main pollution source; the method for calculating the river inflow of the pollution source comprises the following steps:
in the specific embodiment, the calculation formulas of the domestic pollution source, the industrial pollution source and the agricultural non-point source are as follows:
the calculation method of the life pollution source comprises the following steps:
production of domestic pollutants:
wherein: n is a radical ofCity (a city)Is the town population; alpha is alpha1The urban domestic pollution discharge coefficient;
produce pollution for urban life
Wherein: n is a radical ofAgricultural chemicalThe number of rural population; alpha is alpha2Pollution discharge coefficient for rural life;
produce the sewage quantity for the rural area
② domestic sewage treatment rate:
the centralized treatment rate of the urban domestic sewage is equal to the amount of urban domestic sewage taken over by a sewage plant/the amount of urban domestic sewage generated; the rural domestic sewage treatment rate is about 10%, and the pollutant removal rate is about 40%.
③ the river inflow of pollutants in urban and rural life:
wherein: beta is a3Taking 0.4 as the river coefficient; theta2The discharge amount of the domestic pollutant part of the sewage treatment plant; wRaw rThe total river inflow of pollutants for urban life and rural life.
The method for calculating the industrial pollution source comprises the following steps:
industrial pollutant river inflow:
wherein: beta is a3Taking 0.4 as the river coefficient; theta3The discharge amount of the industrial pollution source part of the sewage treatment plant; wI am o 1The pollution amount is industrial; wWorker rThe river inflow amount of industrial pollutants.
The agricultural non-point source calculation method comprises the following steps:
firstly, calculating the pollutant discharge amount of livestock and poultry breeding:
Wlivestock and poultry p＝NLivestock and poultry×α3
Wherein: n is a radical ofLivestock and poultryFor breeding number, alpha3The pollution discharge coefficient of livestock and poultry.
Secondly, calculating pollutant discharge amount of aquaculture:
Waquatic product p＝MAquatic product×α4
Wherein: mAquatic productIs the aquaculture area, alpha4Is the pollution discharge coefficient of aquaculture.
Thirdly, calculating the discharge amount of the farmland pollutants:
Wagricultural p＝M×α5
Wherein: m is the area of cultivated land, alpha5The pollution discharge coefficient of the farmland.
Fourthly, the river inflow of agricultural non-point source pollutants:
Wagricultural non-point sourcer＝WDischarging×β3
WDischarging＝WLivestock and poultry p+WAquatic product p+WAgricultural p
Wherein: beta is a3For the number of entries, 0.4 was used.
And (4) constructing a nested river network hydrodynamic mathematical model in the step (3) which comprises a large watershed hydrodynamic model where the research area is located and a hydrodynamic model of the research area. The construction of the hydrodynamic model of the large watershed comprises the following steps: river network generalization, section generalization, setting of boundary conditions and initial conditions and rating verification of model parameters.
In the boundary conditions, an upstream boundary is set as a flow boundary, and a downstream boundary is set as a water level boundary (if the flow boundary is lacked, digital watershed water system extraction can be adopted, then sub-watersheds are divided, finally, the flow of the afflux river is calculated by a method of producing confluence simulation for each sub-watershed, and the flow value is used as the flow boundary). The initial water level is set above the river bed elevation and the initial flow rate is given a value close to 0. The regional hydrodynamic model construction comprises the following steps: and (4) river network outlines, setting of boundary conditions and rating verification of model parameters. Wherein the water volume calculation boundary of the research area river network model is provided by a large watershed hydrodynamic model.
The step (3) of establishing a water quality mathematical model of the research area comprises the following steps: selecting a pollutant diffusion coefficient, and calibrating initial conditions, boundary conditions and model parameters.
Setting the boundary conditions of the hydrodynamic model in the step (3), if a flow boundary is lacked, extracting a digital watershed water system, dividing sub-watersheds, and finally calculating the flow of the afflux river by a method for producing confluence simulation for each sub-watershed, wherein the flow value is used as the flow boundary, and the specific process is as follows:
i, carrying out depression and heightening treatment on the DEM in the drainage basin;
ii, calculating the water flow direction of the non-depressed area;
iii confluence accumulation amount calculation;
iv setting a threshold value to generate a grid river network;
v grid river network vectorization.
And (4) calibrating according to the synchronous water quantity and water quality monitoring data, wherein the water level and the flow of the hydrologic monitoring data are used for calibrating parameters of a hydrodynamic model, and the water temperature, the pH value, the dissolved oxygen, the permanganate index, the ammonia nitrogen, the total nitrogen and the total phosphorus of the water quality monitoring data are used for calibrating parameters of a water quality model.
Calculating allowable discharge according to a pollution source investigation result and a model calibration result in the step (5), establishing a response relation between the water quality of the control section and the discharge capacity of the sewage discharge outlets according to boundary water quality conditions and design hydrologic conditions through an established mathematical model of the water environment of a river network area of the research area, and then obtaining the allowable discharge capacity of pollutants of each sewage discharge outlet through an allowable discharge capacity calculation method based on the standard reaching of the control section, wherein the calculation method comprises the following steps:
basically, the water quality concentration at the downstream x of the sewage draining outlet is as follows:
in the formula: c' -the water quality concentration after mixing,
C1q is the sewage outlet wastewater concentration, mg/L; amount of waste water, m3/s；
C0、Q0-upstream incoming water concentration, mg/L; flow rate, m3/s；
k is the water quality degradation coefficient, 1/d (obtained by model calibration);
x is the distance from the sewage draining exit, m;
u-flow velocity, m/s.
Allowable pollutant holding amount W considering non-point source pollutionC，
The water quality concentration at a certain position of the downstream of the sewage draining outlet is as follows:
and because:
QP、CPupstream amount of incoming water, m3S; the upstream water concentration, mg/L;
QE、CE-discharge of sewage at sewage outlet m3S; emission concentration, mg/L;
CS、C0-side inflow concentration, mg/L; controlling the river water concentration at the section x, wherein the concentration is mg/L so that:
wherein:
when the water quality at the position of the control section x meets C0When S is the standard of water quality of control section or the requirement of concentration control,
WCi.e. the allowable discharge from the sewage draining outlet to the control section.
Reducing the discharge amount of pollutants according to the obtained allowable discharge amount (namely the allowable discharge amount of the sewage discharge outlet), such as perfecting a town sewage collecting pipe network system and sewage lifting pump station facilities; and determining a water diversion drainage scheme according to the actual condition of the area.
The determination of the designed hydrological conditions is established on the annual daily rainfall data of the regional rainfall station, and the P-III frequency curve is adopted for analysis; and determining the boundary water quality condition according to the water quality target of the functional area where each boundary section is positioned in the model.
The following takes the study area a as an example to illustrate the specific implementation steps of the present invention:
(1) data collection, site survey and monitoring
The water coming from the upstream of the whole watershed of the area A is basically the rainfall confluence of mountains and hills, the hydrological station a and the hydrological station b are two outlets of the whole watershed and are respectively positioned at the northwest corner and the northeast corner of the watershed, the water system in the watershed is developed, and 32 main watercourses are arranged.
(2) Investigating and analyzing the current situation of the regional pollution source and the current situation of the water environment quality, and identifying a main pollution source and a main overproof factor;
and determining a research investigation range by combining the information of the drainage basin DEM and the actual catchment range, and performing refined calculation unit division on the research range, wherein the division condition of the calculation unit of the area A part is shown in table 1.
TABLE 1 part A of the area calculation Unit division
The discharge conditions of the in-zone in-line enterprises are as follows:
TABLE 2 partial in-line enterprise emissions in the area
The total amount of emissions from inline enterprises in the area of study was not large by analysis. Direct point source emissions are not a major source of pollution in the area of interest.
Calculated, the total COD river-entering amount generated by the surface source in the region is 8757.87 tons/year and NH3The river volume of N entering is 974.40 tons/year, and the river volume of TP entering is 127.78 tons/year. The amount of pollutants entering the river is far larger than that of the pollutants entering a straight-line point source, wherein the amount of COD entering the river generated by a surface source is about 106 times that of the pollutants entering the river, and NH is larger than that of the pollutants entering the point source3The river volume N being about a point source159 times, the amount of TP entering the river is about 36 times of the point source. Non-point source fouling is therefore the main source of contamination in the area of investigation.
The construction of sewage pipe networks in the research range is gradually accelerated in recent years, but because residential districts and industrial enterprises in the areas are different in construction times, the pipe networks of old districts are connected disorderly, and the pipe networks are not completely matched in place, so that the domestic sewage cannot be effectively taken over and treated, and the domestic sewage entering the river also becomes a source of pollutants. For example, in the area A-1, COD and NH are generated because the sewage cannot be effectively taken over3The direct river inflow amount in N, TP years is 2307.9 tons, 261.94 tons and 24.09 tons respectively, which respectively account for 89.2 percent, 90.3 percent and 79.0 percent of the river inflow amount in the region.
(3) Constructing nested river network hydrodynamic mathematical model, water quality mathematical model and rating
And constructing a river network hydrodynamic model nested in a basin where the research area is located and the research area, wherein the large network model is used for hydrodynamic simulation of the whole basin and providing a hydrodynamic boundary under a designed hydrological condition for the nested river network of the research area.
a. 15 rivers are generalized in a large river area, and the basic geometric shapes of the generalized cross sections are flat slopes and trapezoidal cross sections. In the boundary condition, the upstream boundary is set as a flow boundary, the downstream boundary is set as a water level boundary, and the actually measured water level value is adopted. The initial water level is set above the river bed elevation and the initial flow rate is given a value close to 0. The roughness of the river bed is rated from n equal to 0.03, and the roughness of the river bed obtained by the rating is 0.025-0.032. The comparison graph of the calculated value and the measured value of a certain point location model is shown in FIG. 2, and the analysis of the large net hydrodynamic simulation calibration error is shown in Table 3.
TABLE 3 analysis of large net hydrodynamic simulation calibration error
| Monitoring point location | Item | Mean measured value | Mean calculated value | Error rate | Qualified (No) |
| Q-1 | Flow rate | 12.98m3/s | 15.63m3/s | 20.42％ | Qualified |
| Q-2 | Flow rate | 23.23m3/s | 24.63m3/s | 6.03％ | Qualified |
| Q-3 | Flow rate | 28.86m3/s | 32.56m3/s | 12.82％ | Qualified |
| Q-4 | Water level | 7.70m | 7.76m | 6cm | Qualified |
b. 10 river channels are generalized in the research area, the water quantity calculation boundary of the river network model of the research area is provided by the large network model, and the river channel roughness rating result is 0.025-0.032. The comparison graph of the calculated value and the measured value of a certain point model is shown in FIG. 3, and the analysis of hydrodynamic force simulation calibration errors in a research area is shown in Table 4.
TABLE 4 analysis of hydrodynamic simulation calibration error in research area
| Monitoring point location | Item | Mean measured value | Mean calculated value | Error rate | Qualified (No) |
| Q-1 | Flow rate | 12.98m3/s | 14.23m3/s | 9.63％ | Qualified |
| Q-2 | Flow rate | 23.23m3/s | 25.15m3/s | 8.27％ | Qualified |
| Q-3 | Flow rate | 28.86m3/s | 33.78m3/s | 17.05％ | Qualified |
| Q-4 | Water level | 7.70m | 7.78m | 8cm | Qualified |
c. Water quality model
The longitudinal diffusion coefficient is selected according to different water flow conditions, and the initial condition is the average value of measured values. The outer boundary adopts the measured value, and the inner boundary refers to the pollution discharge condition of the direct discharge enterprise. Calibration to obtain NH3the-N degradation coefficient is 0.05 to 0.09d-1The TP degradation coefficient is 0.05-0.08 d-1。
(4) Allowable emissions calculation
The result shows that the COD in the research area is calculated by the allowable discharge amount of the research areaMnThe river inflow amount of TP pollutants is less than the allowable discharge amount, but NH3-N river intake is above the allowable emission limit. The allowable pollutant discharge amount based on the qualified cross section water quality is respectively as follows: CODMn10120.5 ton/year NH3-N342.4 ton/year, TP 180.4 ton/year.
(5) In order to further provide a scientific and reasonable research area water environment improvement scheme based on the calculated allowable discharge amount of pollutants from each sewage discharge outlet to the control section on the basis of the above embodiment. A general water quality improvement scheme includes:
(1) the water diversion and drainage scheme can be determined according to the actual condition of the area.
(2) And (4) treating the sewage outlet from the source.
The specific embodiment comprises the following water quality improvement scheme:
the first water quality improvement scheme is to study the water NH of the upstream section of the area under the condition of W1 water diversion3-N water quality extraction functional zone target value, CODMnAnd the TP water quality is taken as a current state monitoring value and simultaneously completely reaches a preset engineering reduction amount.
The second water quality improvement scheme is to study the water NH of the upstream section of the area under the condition that W1 does not introduce water3-N water quality extraction functional zone target value, CODMnAnd the TP water quality is taken as a current state monitoring value and simultaneously completely reaches a preset engineering reduction amount.
The third water quality improvement scheme is to study the water NH of the upstream section of the area under the condition that W1 does not introduce water3-N、CODMnAnd the TP water quality takes the current monitoring value and completely reaches the preset engineering reduction amount at the same time, and the three schemes are compared and selected.
According to the invention, through investigating the current situation of a pollution source, a main pollution source and a main overproof factor are identified, a water environment numerical model is utilized, the water quality change rule of a water body is fully known, the water power and water quality degradation coefficient are obtained through calibration, the allowable pollutant emission amount is calculated by selecting a calculation method based on the standard reaching of a control section, the reduction amount of pollutants is obtained, the conversion from simple pollution control to pollution control and emission reduction is realized, the emission total amount is controlled, and a water ecological system is recovered.
The research of the allowable emission amount has operability and foresight, and can provide reliable reference for further planning work. The water environment is scientifically planned and managed, and social and economic rapid development services are better served.
The numerical model simulation can analyze the water quantity and water quality change and can very conveniently show the migration and diffusion processes of pollutants in the water body, so the numerical model is constructed as a main component in the water environment pollution prevention and control work. The numerical model is used for simulation, the water quality change under specific hydrological conditions and specific working conditions can be conveniently and rapidly predicted, and meanwhile, important reference indexes such as water environment capacity, maximum pollutant discharge and the like can be calculated according to the pollutant degradation coefficient determined by the model, so that a scientific and reasonable research area water environment improvement scheme can be provided.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
While the present invention has been described with reference to the embodiments shown in the drawings, the present invention is not limited to the embodiments, which are illustrative and not restrictive, and it will be apparent to those skilled in the art that various changes and modifications can be made therein without departing from the spirit and scope of the invention as defined in the appended claims.
Claims (7)
1. A total pollutant amount control method based on section water quality control standard is characterized by comprising the following steps:
generalizing river channels in the selected investigation and research range and dividing a calculation unit of river channel sections, and determining a pollution source needing total amount control based on the pollution source in the calculation unit; constructing a river network hydrodynamic mathematical model and a water quality mathematical model;
determining boundary conditions and initial conditions required by a river network hydrodynamic mathematical model and a water quality mathematical model according to the collected and sorted hydrological water quality time data of the boundary of the research area and the upper monitoring point;
carrying out model parameter calibration on the river network hydrodynamic mathematical model and the water quality mathematical model based on water quantity and water quality synchronous monitoring data to obtain a water quality degradation coefficient;
according to the determined main pollution source, pollution overproof factors and a water quality degradation coefficient obtained by model calibration, establishing a response relation between the water quality of the control section and the discharge capacity of the sewage discharge outlet through an established mathematical model of river network hydrodynamic force and a mathematical model of water quality of the research area, and obtaining the allowable discharge capacity of pollutants from each sewage discharge outlet to the control section through a calculation method based on the standard-reaching allowable discharge capacity of the control section; the calculation formula of the allowable pollutant discharge amount WC from each sewage discharge outlet to the control section is as follows:
wherein:
s is a standard for controlling the water quality of the cross section; x is the distance from the downstream of the sewage draining exit;
CSis the side inflow concentration;
a is the area of the cross section and q is the amount of wastewater;
QPis the upstream water flow; cPIs the upstream incoming water concentration;
QEthe sewage discharge amount of a sewage discharge outlet is obtained; cEIs the emission concentration.
2. The method for controlling the total amount of pollutants reaching the water quality standard based on the control section of claim 1, wherein the pollution sources comprise a living pollution source, an industrial pollution source and an agricultural non-point source, and the method for determining the main pollution sources comprises the following steps: calculating the river inflow amount of each pollution source, and comparing the river inflow amount of each pollution source, wherein the highest pollution source is the main pollution source; the method for calculating the river inflow of the pollution source comprises the following steps:
the calculation method of the life pollution source comprises the following steps:
(1) production of domestic pollutants:
urban domestic sewage yield WGeneration of 1The expression of (a) is as follows:
Wgeneration of 1＝NCity (a city)×α1
Wherein: n is a radical ofCity (a city)Is the town population; alpha is alpha1The urban domestic pollution discharge coefficient; wGeneration of 1The pollution amount of the urban life is obtained;
rural sewage yield WGeneration 2The expression of (a) is as follows:
Wgeneration 2＝NAgricultural chemical×α2
Wherein: n is a radical ofAgricultural chemicalThe number of rural population; alpha is alpha2Pollution discharge coefficient for rural life; wGeneration 2The pollution yield in rural areas is high;
(2) domestic sewage treatment rate:
the centralized treatment rate of the urban domestic sewage is equal to the ratio of the amount of the urban domestic sewage taken over by a sewage plant to the amount of the urban domestic sewage;
(3) the calculation expression of the river inflow amount of pollutants in urban life and rural life is as follows:
Wraw r＝WGeneration of 1X (1-centralized treatment rate of urban domestic sewage). times.beta3+θ2+WGeneration 2X (1-rural domestic sewage treatment rate pollutant removal rate) × beta3，
Wherein: beta is a3Is the river entering coefficient; theta2The discharge amount of the domestic pollutant part of the sewage treatment plant; wRaw rThe total river inflow of pollutants for urban life and rural life;
the method for calculating the industrial pollution source comprises the following steps:
(1) river inflow of industrial pollutants:
Wworker r＝WI am o 1X (1-Sewage centralized treatment efficiency). times.beta3+θ3
Wherein: beta is a3Is the river entering coefficient; theta3The discharge amount of the industrial pollution source part of the sewage treatment plant; wI am o 1The pollution amount is industrial; wWorker rThe river inflow amount of industrial pollutants;
the agricultural non-point source calculation method comprises the following steps:
(1) and (3) calculating the discharge amount of livestock and poultry breeding pollutants:
Wlivestock and poultry p＝NLivestock and poultry×α3
Wherein: n is a radical ofLivestock and poultryFor breeding number, alpha3The pollution discharge coefficient of the livestock and poultry is shown;
(2) calculating the pollutant discharge amount of aquaculture:
Waquatic product p＝MAquatic product×α4
Wherein: mAquatic productIs the aquaculture area, alpha4The pollution discharge coefficient of aquaculture;
(3) calculating the discharge amount of farmland pollutants:
Wagricultural p＝M×α5
Wherein: m is the area of cultivated land, alpha5The pollution discharge coefficient of the farmland;
(4) the river inflow amount of agricultural non-point source pollutants:
Wagricultural non-point sourcer＝WDischarging×β3
WDischarging＝WLivestock and poultry p+WAquatic product p+WAgricultural p
Wherein: beta is a3The river entering coefficient.
3. The method of claim 2, wherein the river coefficient β is a factor of total pollutant production based on the control of the water quality of the profile3Take 0.4.
4. The method of claim 1, wherein the river network hydrodynamic mathematical model comprises a large basin hydrodynamic model of a research area and a hydrodynamic model of the research area.
5. The method of claim 1, wherein determining the boundary conditions of the model comprises setting an upstream boundary as a flow boundary and a downstream boundary as a water level boundary; and if the flow boundary is lacked, adopting a digital watershed water system to extract, dividing sub-watersheds, finally calculating the flow of the afflux river by a method of producing confluence simulation for each sub-watershed, and taking the flow value as the flow boundary.
6. The method as claimed in claim 1, wherein the synchronous monitoring data of water quantity and water quality comprises water temperature, pH, dissolved oxygen, permanganate index, ammonia nitrogen, total nitrogen and total phosphorus.
7. The method for controlling the total amount of pollutants reaching the water quality standard based on the control section of claim 1, wherein the method for determining the pollution source needing the total amount control comprises the following steps:
determining a pollution source with the highest emission based on the pollution source in the calculation unit and determining a pollution overproof factor based on the section water quality data;
and determining pollutants needing total amount control according to the pollution source with the highest emission and the pollution overproof factor.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910903473.5A CN110728035B (en) | 2019-09-24 | 2019-09-24 | Pollutant total amount control method based on control of section water quality reaching standard |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910903473.5A CN110728035B (en) | 2019-09-24 | 2019-09-24 | Pollutant total amount control method based on control of section water quality reaching standard |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110728035A CN110728035A (en) | 2020-01-24 |
| CN110728035B true CN110728035B (en) | 2021-09-17 |
Family
ID=69218316
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910903473.5A Active CN110728035B (en) | 2019-09-24 | 2019-09-24 | Pollutant total amount control method based on control of section water quality reaching standard |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110728035B (en) |
Families Citing this family (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111241490A (en) * | 2020-02-13 | 2020-06-05 | 北京市环境保护科学研究院 | Method for processing emission reduction data of sewage treatment plant |
| CN112101693B (en) * | 2020-05-12 | 2023-11-24 | 天津市政工程设计研究总院有限公司 | Method for analyzing current urban river water quality standard based on orthogonal analysis |
| CN111680845B (en) * | 2020-06-15 | 2021-03-19 | 生态环境部南京环境科学研究所 | Water pollution discharge and water environment management system |
| CN111680929B (en) * | 2020-06-16 | 2023-09-01 | 北京嘉博文生物科技有限公司 | Agricultural non-point source pollution evaluation method and application of chemical fertilizer for farmland |
| CN111599418B (en) * | 2020-06-16 | 2023-06-13 | 北京嘉博文生物科技有限公司 | Evaluation method and application of agricultural non-point source pollution of organic waste resources |
| CN113807626A (en) * | 2020-06-17 | 2021-12-17 | 天津大学 | Pollution census-based atmospheric pollutant emission reduction grading method for power and heating power industries |
| CN111914204B (en) * | 2020-07-03 | 2023-08-08 | 河海大学 | Water quality treatment method based on target difference of water quality of cross section examined at upstream and downstream of river network |
| CN111912947B (en) * | 2020-08-04 | 2022-07-29 | 中国农业科学院农业资源与农业区划研究所 | Monitoring method and system for determining river basin agricultural non-point source pollution load |
| CN112116229A (en) * | 2020-09-08 | 2020-12-22 | 深圳市广汇源环境水务有限公司 | Drainage basin water quality scheduling management method, system and platform |
| CN112287613B (en) * | 2020-10-20 | 2022-11-04 | 河海大学 | Pollutant reduction method for watershed water environment control section |
| CN112766669A (en) * | 2021-01-06 | 2021-05-07 | 中国科学院东北地理与农业生态研究所 | Drainage basin pollution limit value emission verification method based on water quality model |
| CN112883644B (en) * | 2021-02-10 | 2022-03-01 | 中国环境科学研究院 | Dynamic water environment management method |
| CN113284005B (en) * | 2021-05-12 | 2022-08-26 | 河海大学 | Sewage treatment system classification method and system |
| CN113361114B (en) * | 2021-06-11 | 2022-05-17 | 中国科学院精密测量科学与技术创新研究院 | Multi-scale non-point source pollutant river entering coefficient measuring and calculating method based on runoff path |
| CN113327056B (en) * | 2021-06-23 | 2022-02-08 | 中国环境科学研究院 | Method for defining water pollution type |
| CN114818211B (en) * | 2022-04-12 | 2024-05-10 | 长江生态环保集团有限公司 | Quantitative simulation technical method for constructing sewage collection system based on MIKE + |
| CN115952246B (en) * | 2022-11-29 | 2023-07-04 | 中国环境科学研究院 | Grid-based factory boundary pollutant measuring and displaying method and system |
| CN116090846A (en) * | 2023-02-16 | 2023-05-09 | 生态环境部华南环境科学研究所(生态环境部生态环境应急研究所) | Water environment control partition dividing method based on pollutant allowable discharge amount |
| CN116431976B (en) * | 2023-06-14 | 2023-08-29 | 江苏省环境工程技术有限公司 | Sewage treatment plant and pipe network efficiency calculation method and device considering external water |
| CN117332959B (en) * | 2023-09-28 | 2024-05-28 | 中建生态环境集团有限公司 | River basin river entering pollutant regulation and control method based on section water quality reaching standard |
| CN117892983B (en) * | 2024-03-14 | 2024-05-28 | 生态环境部华南环境科学研究所(生态环境部生态环境应急研究所) | Method and system for determining offshore area land pollution load distribution |
| CN118134264A (en) * | 2024-05-06 | 2024-06-04 | 生态环境部华南环境科学研究所(生态环境部生态环境应急研究所) | Risk division method and system for black and odorous water body in underdeveloped water system city |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20140015812A (en) * | 2012-07-25 | 2014-02-07 | 인하대학교 산학협력단 | Methodology of gis based automatic schematic diagram creation system for supporting input data generation of 1-dimensional water quality model |
| CN105590024A (en) * | 2015-12-16 | 2016-05-18 | 河海大学 | Method for compiling port atmospheric pollutant discharge inventory based on activities |
| CN107563562A (en) * | 2017-09-08 | 2018-01-09 | 新奥泛能网络科技股份有限公司 | The appraisal procedure and device of gridding urban atmospheric pollution thing discharge |
| CN109242367A (en) * | 2018-11-07 | 2019-01-18 | 山东建筑大学 | A kind of city wastewater treatment rate assessment calculation method |
| CN109389251A (en) * | 2018-10-30 | 2019-02-26 | 江苏省环境科学研究院 | Total amount of pollutant optimizing distribution method based on control section water quality reaching standard |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103218485B (en) * | 2013-04-03 | 2015-10-07 | 郑州大学 | Small watershed environmental capacity computing method under a kind of GIS technology support |
-
2019
- 2019-09-24 CN CN201910903473.5A patent/CN110728035B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20140015812A (en) * | 2012-07-25 | 2014-02-07 | 인하대학교 산학협력단 | Methodology of gis based automatic schematic diagram creation system for supporting input data generation of 1-dimensional water quality model |
| CN105590024A (en) * | 2015-12-16 | 2016-05-18 | 河海大学 | Method for compiling port atmospheric pollutant discharge inventory based on activities |
| CN107563562A (en) * | 2017-09-08 | 2018-01-09 | 新奥泛能网络科技股份有限公司 | The appraisal procedure and device of gridding urban atmospheric pollution thing discharge |
| CN109389251A (en) * | 2018-10-30 | 2019-02-26 | 江苏省环境科学研究院 | Total amount of pollutant optimizing distribution method based on control section water quality reaching standard |
| CN109242367A (en) * | 2018-11-07 | 2019-01-18 | 山东建筑大学 | A kind of city wastewater treatment rate assessment calculation method |
Non-Patent Citations (3)
| Title |
|---|
| 基于控制断面水质达标的允许排放量研究;王烨等;《水电能源科学》;20150731;第33卷(第7期);48-54 * |
| 基尼系数法在黄河中上游流域水污染物总量分配中的应用研究;邱俊永等;《安全与环境学报》;20110630;第11卷(第3期);69-73 * |
| 王烨等.基于控制断面水质达标的允许排放量研究.《水电能源科学》.2015,第33卷(第7期),48-54. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110728035A (en) | 2020-01-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110728035B (en) | Pollutant total amount control method based on control of section water quality reaching standard | |
| CN108664647B (en) | Basin fine management system of integrated water environment model | |
| CN108287950B (en) | Water quality simulation method based on water environment quality target management of control unit | |
| CN108764573B (en) | Inland arid region-oriented multidimensional equilibrium configuration system for water resources | |
| US20230054713A1 (en) | Method for determining contribution rate of pollution load in water quality assessment section of annular river network system based on water quantity constitute | |
| CN107563139A (en) | A kind of point source sudden water pollution source accident source percentage contribution computational methods | |
| CN111259530B (en) | Dynamic water supplementing quantity prediction method for water quality reaching standard of medium-small river basin with multiple river ports | |
| CN109636670B (en) | River type division method based on small-watershed water ecological environment protection stage | |
| Xin et al. | Estimation of non-point source pollution loads with flux method in Danjiangkou Reservoir area, China | |
| CN106599562B (en) | River ecological water demand computational methods based on probability weight FDC methods | |
| CN109598428B (en) | Pollutant reduction and distribution method based on administrative units and water system | |
| CN111310327A (en) | Water environment standard-reaching quantitative design method based on reduction model | |
| CN111444595B (en) | Sewage interception engineering environmental effect evaluation method and device based on SWMM and EFDC coupling model | |
| CN107622363B (en) | Urban water environment quality evaluation method | |
| CN114169266B (en) | River pollution source analysis and pollutant migration calculation method and system | |
| CN111428914B (en) | Method and device for determining water environment treatment engineering scheme and storage medium | |
| Bai et al. | Contribution rates analysis for sources apportionment to special river sections in Yangtze River Basin | |
| CN110728062A (en) | SWMM-based rural non-point source pollution simulation method | |
| CN115758944A (en) | River water power and water quality coupling model and construction method | |
| CN113343604B (en) | Method for defining protection range of river and lake ecological space | |
| CN113284005B (en) | Sewage treatment system classification method and system | |
| CN112287613B (en) | Pollutant reduction method for watershed water environment control section | |
| CN111353718B (en) | Wetland and water replenishing engineering environmental effect evaluation method and device based on SWMM and EFDC | |
| CN113158591B (en) | Method for determining bearing capacity of basin land development and utilization | |
| CN114219371A (en) | Moon-scale basin artificial water pollution load accounting method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |