CN109408972A - The method that computer programing based on Delft3D model calculates environmental capacity of water - Google Patents

The method that computer programing based on Delft3D model calculates environmental capacity of water Download PDF

Info

Publication number
CN109408972A
CN109408972A CN201811274280.XA CN201811274280A CN109408972A CN 109408972 A CN109408972 A CN 109408972A CN 201811274280 A CN201811274280 A CN 201811274280A CN 109408972 A CN109408972 A CN 109408972A
Authority
CN
China
Prior art keywords
water
model
concentration
environmental capacity
river
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.)
Granted
Application number
CN201811274280.XA
Other languages
Chinese (zh)
Other versions
CN109408972B (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.)
Henan Research Academy Of Environment Protection Sciences
Original Assignee
Henan Research Academy Of Environment Protection Sciences
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 Henan Research Academy Of Environment Protection Sciences filed Critical Henan Research Academy Of Environment Protection Sciences
Priority to CN201811274280.XA priority Critical patent/CN109408972B/en
Publication of CN109408972A publication Critical patent/CN109408972A/en
Application granted granted Critical
Publication of CN109408972B publication Critical patent/CN109408972B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C13/00Surveying specially adapted to open water, e.g. sea, lake, river or canal
    • G01C13/002Measuring the movement of open water
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • 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

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Hydrology & Water Resources (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The method that the computer programing that the present invention relates to a kind of based on Delft3D model calculates environmental capacity of water, effective solution of the present invention early stage can only estimate roughly environmental capacity of water, time and spatial accuracy are poor simultaneously, the larger problem of experience dependence;The technical solution of solution includes: to choose certain river to be research range and target section Classification of water Qualities is arranged;Setting model boundary condition includes upstream incoming flow and tributary boundary condition on the way;The input of Delft3D model system is used to export to obtain pollutant concentration value when target section water quality reaches requirement with upper boundary conditions parameter;The calculation of water environmental capacity that related data carries out pollutant is extracted from model, since Delft3D model can utilize specific discrete hydrology-water quality data, be computed output the time be accurate to the second, space is accurate to Hydrodynamic field and the water quality field of rice, it realizes that automation calculates by preparing a computer program, the calculation of water environmental capacity of any water body any time period finally can be achieved.

Description

The method that computer programing based on Delft3D model calculates environmental capacity of water
Technical field
The present invention relates to the environmental managements such as pollution emission reduction and pollutant emission factor fields, and in particular to one kind is based on The method that the computer programing of Delft3D model calculates environmental capacity of water.
Background technique
Environmental carrying capacity is one of Issues of Basic Theory of environmental science, is one of important actual application problem of environmental management Environmental capacity of water is the important component of environmental carrying capacity, is one of the core content of capacity total amount technical system.With China Water environment management system realizes quality in watershed target tube from concentration control, transformation of the target overall control to Capacity-Overall Control Reason receives dirty red line management with Water Functional Zone limitation, and environmental capacity of water is theoretical and the importance of Research on Calculation more highlights.In The initial calculation of water environmental capacity method of state first is that directly establish its calculation formula from definition, this meter can be claimed Calculation method is equation;With going deep into for research, and Water Environment Mathematical Model formula is combined, i.e., is defined based on environmental capacity of water And Water Environment Mathematical Model, the calculation of water environmental capacity formula under certain condition is derived, hydrodynamic model and water quality model are based on Required parameters in calculation of water environmental capacity formula are calculated, and then substitutes into formula and calculates environmental capacity of water.By what is derived Calculation of water environmental capacity formula is also referred to as environmental capacity of water model, and (careful wait has used " environmental capacity of water mathematical modulo earliest Type ", the concept of " environmental capacity of water model "), therefore equation also known as environmental capacity of water modelling.
Very much, but its citation form is equal for the computation model (calculation formula) of environmental capacity of water are as follows: environmental capacity of water=dilution Capacity+from net capacity+migration capacity.With gradually deep, the calculation of water environmental capacity formula gradual perfection of research, and according to not With pollutant, different water bodys and establish different calculation formula.Equation may be considered most basic in all kinds of methods Method, the calculating of other kinds method is also based on calculation of water environmental capacity formula, due to the general side of change of formula under different operating conditions There are difference for method, and as illustrated, equation has different expression ways, this also can make the calculated result of environmental capacity of water It is different;Especially early stage, due to environmental capacity of water definition, i.e. difference itself is larger, then calculated result certainly exists difference, Therefore early stage can only estimate roughly environmental capacity of water, while time and spatial accuracy are poor, experience dependence is larger.
Recently as the development of numerical model and computer technology, Hydrodynamic field, the simulation of water quality field and prediction model It makes great progress, by taking Delft3D model as an example, inputs data (flow, water in particular point (upstream, downstream, tributary etc.) Position, water quality etc.), model can simulate the Hydrodynamic field of entire water body and water quality field by these discrete data, when Between resolution ratio can be as accurate as the second, spatial resolution can be as accurate as 1m,
Then it is calculated by the realization automation for preparing a computer program, finally realizes the water environment of any water body any time Calculation of capacity, therefore, the present invention provide a kind of side of computer programing calculating environmental capacity of water based on Delft3D model Method.
Summary of the invention
For above situation, for the defect for overcoming the prior art, the present invention provides a kind of calculating based on Delft3D model The method of machine sequencing calculating environmental capacity of water.Effective solution early stage can only estimate roughly environmental capacity of water It calculates, while time and spatial accuracy are poor, the larger problem of experience dependence.
The method that computer programing based on Delft3D model calculates environmental capacity of water, which is characterized in that calculate step It is as follows:
Step 1: choosing certain river is research range, and by taking the environmental carrying capacity of river X calculates as an example, the river target side Face function water quality is F class water quality, when calculating the section and meet F class water standard, the X environmental capacity of water in the river;
Step 2: carrying out the drafting of orthonormal curvilinear grid to the river according to the research range of model, obtaining one and contain Grid number is the grid of A, and the mean depth in the river is set as N;
Step 3: setting model boundary condition includes upstream incoming flow boundary condition, on the way tributary boundary condition;
Step 4: a point of release, flow, the concentration of emission of point of release is additionally arranged in the upstream of the river research range Be respectively set to M, L, use Delft3D model system by input upstream incoming flow boundary condition, on the way tributary boundary condition with And the flow M and concentration of emission L of point of release X, output obtain the concentration of target assessment section X, point of release flow M are determined, no Disconnected adjustment concentration of emission L numerical value, so that the X concentration in target section part reaches exactly to F class standard value;
Step 6: the calculating that related data carries out the environmental capacity of water of X is extracted from model after model calculates, it is public Formula is as follows:
Pass through calculation formula:
ABS is the symbol of an absolute value in formula, and expression takes absolute value to subsequent calculated value;Assuming that being provided with n+1 A data, time are denoted as respectively: t0,t1…tn, unit: min, corresponding flow are respectively as follows: q0,q1…qn, unit: m3/ s is right The concentration answered is respectively as follows: c0,c1…cn, unit: kg/m3
Wherein, first data is the data of model time started, and nth data is the data of model end time.Meter The value calculated, divided by 1000, obtains the X environmental carrying capacity that unit is t/d divided by simulated time number of days.
Preferably, boundary condition set in the step 3 is the flow of pollutant X and the prison of pollutant X concentration Measured value.
Preferably, the related data extracted from model in the step 6 refers to: t0,t1…tn, unit: min, often The flow and concentration of pollutant X in a period.
(1) since Delft3D model can be using specific discrete hydrology-water quality data, being computed the output time is accurate to Second, space is accurate to Hydrodynamic field and the water quality field of rice, realizes that automation calculates by preparing a computer program, final to can be achieved The calculation of water environmental capacity of any water body any time period;
(2) present invention can fast and accurately calculate the environmental capacity of water of water body of interest, be the disposal of pollutants of water body periphery Formulation with water environment relevant policies provides scientific basis.
Detailed description of the invention
Fig. 1 is the river model grid schematic diagram;
Fig. 2 is the river topography schematic diagram;
Fig. 3 is model modeling flow chart;
Fig. 4 is the model boundary distribution map;
Fig. 5 is downstream targets monitoring point and the modal position relation schematic diagram.
Specific embodiment
For the present invention aforementioned and other technology contents, feature and effect, in following cooperation with reference to figures 1 through Fig. 5 pairs In the detailed description of embodiment, can clearly it present.The structure content being previously mentioned in following embodiment is attached with specification Figure is reference.
Each exemplary embodiment of the invention is described below with reference to accompanying drawings.
Embodiment 1, the method that the computer programing based on Delft3D model calculates environmental capacity of water, which is characterized in that Steps are as follows for calculating:
Step 1: choosing certain river is research range, by taking the environmental carrying capacity of the river ammonia nitrogen calculates as an example, it is assumed that the river Target end face function water matter is V class water quality, when calculating the section and meet V class water standard, the ammonia nitrogen water environment in the river Capacity, the V class water quality are the standard in " water environment quality standard " GB3838-2002, " surface water environment quality mark It is quasi- " surface water quality is classified as I class, II class, Group III, IV class, V class in GB3838-2002;
Step 2: carrying out the drafting of orthonormal curvilinear grid to the river according to the research range of model, obtaining one and contain The grid that grid number is 5971, the mean depth in the river are set as 1 meter;
Step 3: setting model boundary condition includes upstream incoming flow boundary condition, on the way tributary boundary condition, the perimeter strip Part be different upstream incoming flows and on the way in tributary flow condition and condition of water quality monitor value;
Step 4: a point of release, flow, the concentration of emission of point of release is additionally arranged in the upstream of the river research range Be respectively set to M, L, use Delft3D model system by input upstream incoming flow boundary condition, on the way tributary boundary condition with And the flow M and concentration of emission L of point of release ammonia nitrogen, output obtains the concentration of target assessment section ammonia nitrogen, and point of release flow M is true It is set to 5m3/ s constantly adjusts concentration of emission L numerical value, so that the ammonia nitrogen concentration in target section part reaches exactly to V class standard value;
Step 5: the calculating that related data carries out the environmental capacity of water of ammonia nitrogen is extracted from model after model calculates, Pass through calculation formula:
ABS is the symbol of an absolute value in formula, and expression takes absolute value to subsequent calculated value, the flow of ammonia nitrogen and dense The product adduction of degree is exactly the environmental carrying capacity of the section ammonia nitrogen;
Assuming that being provided with n+1 data, the time is denoted as respectively: t0,t1…tn, unit: min, corresponding flow are respectively as follows: q0,q1…qn, unit: m3/ s, corresponding concentration are respectively as follows: c0,c1…cn, unit: kg/m3
Wherein, first data is the data of model time started, and nth data is the data of model end time.Meter The value calculated, divided by 1000, obtains the ammonia nitrogen environmental carrying capacity that unit is t/d divided by simulated time number of days;
The related data extracted from model in step 5 refers to: t0,t1…tn, unit: min, it is dirty in each period The flow and concentration of object ammonia nitrogen are contaminated, is then calculated by formula by t0,t1…tn, unit: min, the ammonia nitrogen ring in all periods Border capacity sums up when show that target section meets V class water standard, the environmental capacity of water of the section ammonia nitrogen.
S=(2 × 0.002+6 × 0.3 × 0.0015+5 × 0.0017) × 1440 × 60/1000=1.3
In formula 2 × 0.002 be step 3 middle and upper reaches inlet flow rate and water quality product;6 × 0.3 × 0.0015 is six The product of tributary flow and water quality on the way;5 × 0.0017 for the flow of point of release and when target section ammonia nitrogen concentration value reaches V class The product of the concentration of emission of point of release ammonia nitrogen when standard.
Since Delft3D model can be using specific discrete hydrology-water quality data, being computed the output time is accurate to the second, Space is accurate to Hydrodynamic field and the water quality field of rice, realizes that automation calculates by preparing a computer program, final to can be achieved to appoint The calculation of water environmental capacity of meaning water body any time period.The water environment that the present invention can fast and accurately calculate water body of interest is held Amount, provides scientific basis for the disposal of pollutants of water body periphery and the formulation of water environment relevant policies.
Embodiment described above is not intended to limit the scope of the present invention, before not departing from design concept of the present invention It puts, the various changes and improvements that technical staff described in this field makes technical solution of the present invention should all be included in the present invention Claims determine protection scope in.

Claims (3)

1. the method that the computer programing based on Delft3D model calculates environmental capacity of water, which is characterized in that calculate step such as Under:
Step 1: choosing certain river is research range, and by taking the environmental carrying capacity of river X calculates as an example, the river target end face function Can water quality be F class water quality, when calculating the section and meet F class water standard, the X environmental capacity of water in the river;
Step 2: carrying out the drafting of orthonormal curvilinear grid to the river according to the research range of model, obtaining one and contain grid Number is the grid of A, and the mean depth in the river is set as N;
Step 3: setting model boundary condition includes upstream incoming flow boundary condition, on the way tributary boundary condition;
Step 4: a point of release, the flow of point of release, concentration of emission difference is additionally arranged in the upstream of the river research range It is set as M, L, using Delft3D model system by inputting upstream incoming flow boundary condition, on the way tributary boundary condition and releasing The flow M and concentration of emission L of point X are put, output obtains the concentration of target assessment section X, and point of release flow M is determined, constantly adjusts Whole concentration of emission L numerical value, so that the X concentration in target section part reaches exactly to F class standard value;
Step 5: extracting the calculating that related data carries out the environmental capacity of water of X from model after model calculates, passing through meter Calculate formula:
ABS is the symbol of an absolute value in formula, and expression takes absolute value to subsequent calculated value;
Assuming that being provided with n+1 data, the time is denoted as respectively: t0,t1…tn, unit: min, corresponding flow are respectively as follows: q0, q1…qn, unit: m3/ s, corresponding concentration are respectively as follows: c0,c1…cn, unit: kg/m3
Wherein, first data is the data of model time started, and nth data is the data of model end time.It calculates The value come, divided by 1000, obtains the X environmental carrying capacity that unit is t/d divided by simulated time number of days.
2. the method that the computer programing according to claim 1 based on Delft3D model calculates environmental capacity of water, It is characterized in that, set boundary condition is the flow of pollutant X and the monitor value of pollutant X concentration in the step 3.
3. special according to the method that the computer programing described in claim 1 based on Delft3D model calculates environmental capacity of water Sign is that the related data extracted from model in the step 6 refers to: t0,t1…tn, unit: min, each period The flow and concentration of interior pollutant X.
CN201811274280.XA 2018-10-30 2018-10-30 Method for calculating water environment capacity through computer programming based on Delft3D model Active CN109408972B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811274280.XA CN109408972B (en) 2018-10-30 2018-10-30 Method for calculating water environment capacity through computer programming based on Delft3D model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811274280.XA CN109408972B (en) 2018-10-30 2018-10-30 Method for calculating water environment capacity through computer programming based on Delft3D model

Publications (2)

Publication Number Publication Date
CN109408972A true CN109408972A (en) 2019-03-01
CN109408972B CN109408972B (en) 2021-03-16

Family

ID=65469720

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811274280.XA Active CN109408972B (en) 2018-10-30 2018-10-30 Method for calculating water environment capacity through computer programming based on Delft3D model

Country Status (1)

Country Link
CN (1) CN109408972B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024051131A1 (en) * 2022-09-09 2024-03-14 河北先河环保科技股份有限公司 Water environmental capacity calculation method and apparatus, electronic device, and storage medium

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103218485A (en) * 2013-04-03 2013-07-24 郑州大学 Method for calculating small watershed environmental capacity under support of GIS (Geographic Information System) technology
CN104199887A (en) * 2014-08-25 2014-12-10 北京三易思创科技有限公司 Water environment analysis method and system based on geographic information system and Delft3D-FLOW calculation engine
CN105320828A (en) * 2014-07-04 2016-02-10 中国环境科学研究院 Pollutant total-amount control method under large-scale planar grid condition
CN107657912A (en) * 2017-09-26 2018-02-02 天津大学 A kind of river water environmental capacity determining methods based on Mike11
CN108549744A (en) * 2018-03-14 2018-09-18 西安理工大学 A kind of river water environmental capacity allocation method

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103218485A (en) * 2013-04-03 2013-07-24 郑州大学 Method for calculating small watershed environmental capacity under support of GIS (Geographic Information System) technology
CN105320828A (en) * 2014-07-04 2016-02-10 中国环境科学研究院 Pollutant total-amount control method under large-scale planar grid condition
CN104199887A (en) * 2014-08-25 2014-12-10 北京三易思创科技有限公司 Water environment analysis method and system based on geographic information system and Delft3D-FLOW calculation engine
CN107657912A (en) * 2017-09-26 2018-02-02 天津大学 A kind of river water environmental capacity determining methods based on Mike11
CN108549744A (en) * 2018-03-14 2018-09-18 西安理工大学 A kind of river water environmental capacity allocation method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
周吉: "无锡市区水环境容量估算与排污总量分配", 《水利科技与经济》 *
王菊翠等: "泾河流域陕西段水环境容量估算和分配模型研究", 《水土保持研究》 *
邬彬等: "深圳水环境容量及其承载力评价", 《环境科学研究》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024051131A1 (en) * 2022-09-09 2024-03-14 河北先河环保科技股份有限公司 Water environmental capacity calculation method and apparatus, electronic device, and storage medium

Also Published As

Publication number Publication date
CN109408972B (en) 2021-03-16

Similar Documents

Publication Publication Date Title
CN109657418B (en) Lake water environment capacity calculation method based on MIKE21
CN109815608B (en) Ecological coupling scheduling analysis method for water quality and water volume of shallow lake group
CN107091911B (en) river water quality prediction method
Shen et al. An overview of research on agricultural non-point source pollution modelling in China
CN108256177B (en) Parameter optimization method and system for river water sand model
CN106202746B (en) Simulate the Yeh- multi-level finite element modeling method of Water in Porous Medium stream Darcy velocity
CN109101781B (en) The calculation method of pollution sources contribution proportion in a kind of Complex River
Keupers et al. Development and testing of a fast conceptual river water quality model
CN112182866A (en) Water quality rapid simulation method and system based on water environment coupling model
CN111259607B (en) River and lake transition region hydrological boundary defining method
Wagener et al. Dynamic identifiability analysis of the transient storage model for solute transport in rivers
Kremer et al. Simulating property exchange in estuarine ecosystem models at ecologically appropriate scales
Butenschön et al. Sensitivity of a marine coupled physical biogeochemical model to time resolution, integration scheme and time splitting method
CN109460605B (en) A method of predicting large-scale low-lift pump flow
CN109711107A (en) Water resource burst water pollutant Diffusion Simulation system
CN109185211A (en) A kind of pumping plant unit pressure pulsation prediction method instantaneously calculated based on Fluid Mechanics Computation
Rauschenbach Modeling, control, and optimization of water systems
Jia et al. The development of a multi-species algal ecodynamic model for urban surface water systems and its application
Gumiere et al. Evaluating the impact of the spatial distribution of land management practices on water erosion: case study of a Mediterranean catchment
CN109408972A (en) The method that computer programing based on Delft3D model calculates environmental capacity of water
CN110847112B (en) River flood discharge early warning method based on hydraulics simulation
Zhao et al. Development of a zoning-based environmental–ecological coupled model for lakes: a case study of Baiyangdian Lake in northern China
CN116432439A (en) Urban river sewage receiving capacity planning method and system based on numerical simulation
CN110322140A (en) A kind of tidal reach calculation of water environmental capacity method based on MIKE21
Zhang et al. Discretization approach in integrated hydrologic model for surface and groundwater interaction

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