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 PDF

Info

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
Application number
CN201910903473.5A
Other languages
Chinese (zh)
Other versions
CN110728035A (en
Inventor
徐芸蔚
李一平
施媛媛
黄亚男
程一鑫
程月
魏蓥蓥
赵晓磊
王凯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hohai University HHU
Original Assignee
Hohai University HHU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hohai University HHU filed Critical Hohai University HHU
Priority to CN201910903473.5A priority Critical patent/CN110728035B/en
Publication of CN110728035A publication Critical patent/CN110728035A/en
Application granted granted Critical
Publication of CN110728035B publication Critical patent/CN110728035B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393Score-carding, benchmarking or key performance indicator [KPI] analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10Services
    • G06Q50/26Government or public services
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A20/00Water conservation; Efficient water supply; Efficient water use
    • Y02A20/152Water 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)
  • Sewage (AREA)
  • Feedback Control In General (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

Pollutant total amount control method based on control of section water quality reaching standard
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:
Figure BDA0002212560130000054
wherein: n is a radical ofCity (a city)Is the town population; alpha is alpha1The urban domestic pollution discharge coefficient;
Figure BDA0002212560130000056
produce pollution for urban life
Figure BDA0002212560130000055
Wherein: n is a radical ofAgricultural chemicalThe number of rural population; alpha is alpha2Pollution discharge coefficient for rural life;
Figure BDA0002212560130000051
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:
Figure BDA0002212560130000052
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:
Figure BDA0002212560130000053
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:
Figure BDA0002212560130000071
in the formula: c' -the water quality concentration after mixing,
Figure BDA0002212560130000072
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:
Figure BDA0002212560130000081
and because:
Figure BDA0002212560130000082
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:
Figure BDA0002212560130000083
wherein:
Figure BDA0002212560130000084
Figure BDA0002212560130000085
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,
Figure BDA0002212560130000091
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
Figure BDA0002212560130000101
The discharge conditions of the in-zone in-line enterprises are as follows:
TABLE 2 partial in-line enterprise emissions in the area
Figure BDA0002212560130000102
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:
Figure FDA0003036692750000011
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;
Figure FDA0003036692750000021
k is a water quality degradation coefficient obtained by calibration;
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.beta32+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.beta33
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.
CN201910903473.5A 2019-09-24 2019-09-24 Pollutant total amount control method based on control of section water quality reaching standard Active CN110728035B (en)

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 (21)

* Cited by examiner, † Cited by third party
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

Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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
CN103276686A (en) Method of determining downstream design flood of cascade reservoirs
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
CN111428914B (en) Method and device for determining water environment treatment engineering scheme and storage medium
CN110765213A (en) Method for compiling emission list (dynamic list) of pollution sources in surface water basin
CN104679993A (en) Assimilative capacity calculating method based on binary water circulation
Bai et al. Contribution rates analysis for sources apportionment to special river sections in Yangtze River Basin
CN114169266A (en) River pollution source analysis and pollutant migration calculation method and system
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

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