CN100397405C - Land source sewage discharging quantity inversion method based on variation algorithm - Google Patents
Land source sewage discharging quantity inversion method based on variation algorithm Download PDFInfo
- Publication number
- CN100397405C CN100397405C CNB2005101238253A CN200510123825A CN100397405C CN 100397405 C CN100397405 C CN 100397405C CN B2005101238253 A CNB2005101238253 A CN B2005101238253A CN 200510123825 A CN200510123825 A CN 200510123825A CN 100397405 C CN100397405 C CN 100397405C
- Authority
- CN
- China
- Prior art keywords
- yuan
- flow rate
- blowdown flow
- sewage discharging
- sea
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 21
- 239000010865 sewage Substances 0.000 title abstract description 15
- 238000007599 discharging Methods 0.000 title abstract description 7
- 239000003344 environmental pollutant Substances 0.000 claims abstract description 13
- 231100000719 pollutant Toxicity 0.000 claims abstract description 13
- 230000008569 process Effects 0.000 claims abstract description 6
- 230000006870 function Effects 0.000 claims description 12
- 238000009792 diffusion process Methods 0.000 claims description 5
- 238000004088 simulation Methods 0.000 claims description 3
- 230000032258 transport Effects 0.000 claims description 3
- 238000010230 functional analysis Methods 0.000 claims description 2
- 230000010354 integration Effects 0.000 claims description 2
- 238000013507 mapping Methods 0.000 claims description 2
- 238000012544 monitoring process Methods 0.000 abstract description 23
- 238000011160 research Methods 0.000 abstract description 3
- 230000008901 benefit Effects 0.000 abstract description 2
- 238000005070 sampling Methods 0.000 description 16
- 238000005457 optimization Methods 0.000 description 3
- 238000003912 environmental pollution Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 239000013535 sea water Substances 0.000 description 1
- 238000003902 seawater pollution Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000004800 variational method Methods 0.000 description 1
- 239000003643 water by type Substances 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Sampling And Sample Adjustment (AREA)
Abstract
The present invention belongs to the technical field of monitoring and managing land source sewage discharging quantity and provides a monitoring method for sea environment pollution, and the monitoring method is characterized in that the land source sewage discharging quantity can be obtained inverting the pollutant concentration monitored by sea samples. First of all, an adjacent sea dynamics numerical value mode is established; afterwards, a variation method is introduced into the research on sea environment problems and the inversion problem from sea sample parameters to sewage discharging parameters is converted to a variational problem under a restricted condition (numerical model equation); then each parameter in the computation process is decided to realize the optimum estimation of the algorithm. Therefore, the corresponding land source sewage discharging quantity is obtained by inverting the pollutant concentration monitored by sea samples. The present invention has the advantages of strong practicability and high versatility, and is mainly used for monitoring and managing the land source sewage discharging quantity.
Description
Technical field
The inversion method that pollutant levels and genetic algorithm by marine sampling monitoring obtains the Lu Yuan blowdown flow rate.
Background technology
In marine environmental management and pollution control work, obtaining the pollution from land-based sources source data is core link.Yet, because the extra large pollution source monitoring of row difficulty, the investigation of Lu Yuan blowdown flow rate itself is one and has that quite big probabilistic process etc. is all multifactor in addition, cause Lu Yuan blowdown flow rate data scarcity, the low of inferior quality problem of data quality, greatly hampered monitoring, improvement and the research of marine environmental pollution.
At this present situation, the seawater that the present invention proposes the sewage draining exit surrounding waters carries out this new monitoring method of sample analysis (hereinafter to be referred as " marine sampling monitoring "), it has following major advantage: (1) is from largely going up, marine sampling can less consider not participate in the influence of the solid shape polluter of ocean diffusion, obtain more accurately marine environment is produced the dissolving of considerable influence and is suspended in marine pollutant discharge amount, thereby more can reflect Lu Yuan blowdown information.(2) the sea pollution parameter is the integral mean to lasting blowdown, and the sewage draining exit blowdown flow rate is risen and fallen certain filter action, helps obtaining average emission information more reliably by small sample.(3) marine sampling monitoring can be avoided the monitoring difficulty that sewage draining exit position, landform etc. are brought, and reduces the sampling risk, improves observed efficiency.(4) comparatively speaking, marine sampling monitoring all is easier at aspects such as instrument and equipment and operation measures implement than the monitoring of sewage draining exit.Therefore, marine sampling monitoring is the effective means of monitoring land-sourced pollutant discharging.
Yet,,, directly infer that from the sampling analysis result actual blowdown flow rate can produce very mistake if do not consider the dynamic process of ocean though the seawater pollution substrate concentration of marine sampling monitoring is closely related with the Lu Yuan blowdown flow rate.Have only to have solved, could realize really that marine sampling monitoring replaces the monitoring of Lu Yuan sewage draining exit from the anti-problem that pushes away parameters sewage of marine sampling parameters.Rely on the ocean dynamics relation also not have the precedent of application from the anti-problem that pushes away parameters sewage of marine sampling parameters the ocean, but tentatively developed similar method in air monitoring, it is generalized in the ocean is fully feasible.
Summary of the invention
Purpose of the present invention is exactly by variational technique, develop a kind of practical, versatility is high, can obtain the marine environmental pollution monitoring method of Lu Yuan blowdown flow rate by the pollutant levels inverting of marine sampling monitoring.
The present invention at first sets up coastal ocean dynamics numerical model, then the variational method is incorporated in the research of marine environment problem, the variational problem that is converted into (numerical model equation) under certain constraint condition by marine sampling parameters to the inverse problem of parameters sewage, determine each parameter in the computation process then, realize the optimal estimation of this algorithm, thereby obtain corresponding Lu Yuan blowdown flow rate by the pollutant levels inverting of marine sampling monitoring.
Description of drawings
Concrete implementing procedure figure of the present invention
Embodiment
1. the dynamic process that Pollutants Diffusion is transported is regarded the function of Lu Yuan blowdown flow rate S as, language with functional analysis is exactly as operator the diffusion transport equation of pollutant, the distribution of sea pollution thing is a kind of mapping of this operator to the Lu Yuan blowdown flow rate, to observe concentration field and utilize dynamic mode to compare from the concentration field that certain Lu Yuan blowdown flow rate distribution field calculates, definition distance function J (" distance " of marine concentration field (Cm) and actual observation concentration (Co) that generally is simulation is to integration in finite time or summation), then J is the functional of Lu Yuan blowdown flow rate S, ask the minimum Lu Yuan blowdown flow rate of functional is distributed, obtain variational problem.Dynamic mode as constraint condition, by the size of continuous adjustment Lu Yuan blowdown flow rate, is made distance function J reach minimum, thereby reaches the purpose of inverting.
2. adopt accurate Newton's algorithm to realize the solution procedure of above-mentioned optimization problem: the governing equation of employing logarithm value pattern is directly write the method for adjoint mode program, by to the adjoint mode reverse integral, can get the gradient of objective function to S; Then with the negative gradient direction as the direction of search, try to achieve and make J reach minimum Lu Yuan blowdown flow rate S.
The present invention wish the Lu Yuan blowdown flow rate that obtains by inverting by numerical simulation calculate marine observation station pollutant levels (Cik) the most approaching with the pollutant levels (Cok) of actual marine sampling monitoring, therefore get
As objective function (wherein k is the sewage draining exit number).Because in fact a certain regional blowdown flow rate can only be a nonnegative value, it is S 〉=0, situation according to inversion result, can with due regard to construct penalty function about blowdown flow rate S, the optimization problem of former cause inequality constrain is converted into the unconstrained optimization problem of finding the solution a series of penalty functions, utilizes accurate newton's descent algorithm to find the solution as objective function with it.
Claims (1)
1. inversion method based on the Lu Yuan blowdown flow rate of variational algorithm, it is characterized in that the dynamic process that Pollutants Diffusion transports is regarded as the function of Lu Yuan blowdown flow rate S, language with functional analysis is expressed, be exactly as operator the diffusion transport equation of pollutant, the distribution of sea pollution thing as of a kind of mapping of this operator to the Lu Yuan blowdown flow rate, to observe concentration field and utilize dynamic mode to compare from the concentration field that certain Lu Yuan blowdown flow rate distribution field calculates, definition distance function J for " distance " of the marine concentration field Cm of simulation and actual observation concentration C o to integration in finite time or summation, then J is the functional of Lu Yuan blowdown flow rate S, so just can distribute in the hope of making the minimum Lu Yuan blowdown flow rate of functional, obtain variational problem, again with dynamic mode as constraint condition, size by continuous adjustment Lu Yuan blowdown flow rate, make distance function J reach minimum, and the governing equation that wherein adopts the logarithm value pattern is directly write the method for adjoint mode program, by to the adjoint mode reverse integral, obtain the gradient of objective function to S, then with the negative gradient direction as the direction of search, try to achieve and make J reach minimum Lu Yuan blowdown flow rate S, thereby reach the purpose of inverting.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005101238253A CN100397405C (en) | 2005-11-23 | 2005-11-23 | Land source sewage discharging quantity inversion method based on variation algorithm |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005101238253A CN100397405C (en) | 2005-11-23 | 2005-11-23 | Land source sewage discharging quantity inversion method based on variation algorithm |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1794251A CN1794251A (en) | 2006-06-28 |
CN100397405C true CN100397405C (en) | 2008-06-25 |
Family
ID=36805683
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2005101238253A Expired - Fee Related CN100397405C (en) | 2005-11-23 | 2005-11-23 | Land source sewage discharging quantity inversion method based on variation algorithm |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100397405C (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101604019B (en) * | 2009-07-13 | 2012-07-04 | 中国船舶重工集团公司第七一五研究所 | Method for quickly calculating characterization and transfer of uncertainties of marine environment and sound fields |
CN102254081B (en) * | 2011-04-06 | 2014-11-05 | 北方工业大学 | Statistical location method of sudden river water pollution based on stochastic probability model |
CN102156817B (en) * | 2011-04-11 | 2013-04-03 | 国家海洋信息中心 | Rapid prediction method for drifting path and diffusion of marine spilled oil |
CN107526880B (en) * | 2017-08-03 | 2021-02-23 | 中国科学院南京地理与湖泊研究所 | Method for calculating upstream river-entering pollution flux based on river water quality observation |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5813798A (en) * | 1997-03-28 | 1998-09-29 | Whiffen; Greg | Piecewise continuous control of groundwater remediation |
CN1527119A (en) * | 2003-09-24 | 2004-09-08 | 成 程 | Genetic algorithm design of Er-doping fiber and Er-doping fiber amplifier |
-
2005
- 2005-11-23 CN CNB2005101238253A patent/CN100397405C/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5813798A (en) * | 1997-03-28 | 1998-09-29 | Whiffen; Greg | Piecewise continuous control of groundwater remediation |
CN1527119A (en) * | 2003-09-24 | 2004-09-08 | 成 程 | Genetic algorithm design of Er-doping fiber and Er-doping fiber amplifier |
Non-Patent Citations (6)
Title |
---|
地下水污染强度反演的数值方法. 李功胜,王孝勤,高希报.地下水,第26卷第2期. 2004 |
地下水污染强度反演的数值方法. 李功胜,王孝勤,高希报.地下水,第26卷第2期. 2004 * |
污染源浓度分布分数阶模型及其解. 蒋晓芸,徐明瑜.山东大学学报(理学版),第39卷第3期. 2004 |
污染源浓度分布分数阶模型及其解. 蒋晓芸,徐明瑜.山东大学学报(理学版),第39卷第3期. 2004 * |
遗传算法在点源扩散浓度反演排放源强中的应用. 陈军明,徐大海,朱蓉.气象,第28卷第9期. 2002 |
遗传算法在点源扩散浓度反演排放源强中的应用. 陈军明,徐大海,朱蓉.气象,第28卷第9期. 2002 * |
Also Published As
Publication number | Publication date |
---|---|
CN1794251A (en) | 2006-06-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109657418B (en) | Lake water environment capacity calculation method based on MIKE21 | |
Sadeghian et al. | Improving in-lake water quality modeling using variable chlorophyll a/algal biomass ratios | |
Gordon et al. | LOICZ biogeochemical modelling guidelines | |
CN114881544B (en) | Reservoir flow regulation and control method and device, electronic equipment and storage medium | |
Chen et al. | Assimilating multi-source data into a three-dimensional hydro-ecological dynamics model using ensemble Kalman filter | |
CN118657091B (en) | Ocean pollution degree assessment method and system based on pollutant diffusion model | |
CN109598428A (en) | One kind being based on administration cell and water-based pollutant fluxes distribution method | |
CN115422850B (en) | Method for judging pollution contribution rate of river inflow river drain outlet based on EFDC and edge calculation | |
Ascione Kenov et al. | Advances in modeling of water quality in estuaries | |
CN100397405C (en) | Land source sewage discharging quantity inversion method based on variation algorithm | |
Zhao et al. | Development of a zoning-based environmental–ecological coupled model for lakes: a case study of Baiyangdian Lake in northern China | |
CN119720731A (en) | Water environment management method and platform based on water environment mathematical model and artificial intelligence | |
Sheng et al. | Skill assessment of an integrated modeling system for shallow coastal and estuarine ecosystems | |
CN110322140A (en) | MIKE 21-based tidal river reach water environment capacity calculation method | |
Huggett et al. | Modelling the impact of river flow, macronutrients and solar radiation on the eutrophication status of small shallow estuaries | |
Ekdal et al. | Application of WASP and SWAT models for a Mediterranean coastal lagoon with limited seawater exchange | |
CN100401308C (en) | Method for inverting calculating land source draining away sewage quantity based on genetic calculating | |
Mannina et al. | A hydrodynamic water quality model for propagation of pollutants in rivers | |
Håkanson | Factors and criteria to quantify coastal area sensitivity/vulnerability to eutrophication: presentation of a sensitivity index based on morphometrical parameters | |
Moore et al. | Oxygen dynamics in a well mixed estuary, the lower Delaware Bay, USA | |
Kim et al. | Three-dimensional water quality modelling for tidal lake and coastal waters with ROMS-ICM | |
Jang et al. | Spatiotemporal algal bloom prediction of geum river, Korea using the deep learning models in company with the EFDC model | |
Alam et al. | Application of a mathematical model for nutrient pollution assessment in a small catchment of the Latrobe River Basin, Australia | |
Cheng | Assessment system of water conservancy project impact on environmental pollution based on HSPF model | |
Plew et al. | CLUES Estuary-a tool for predicting estuary water quality |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C19 | Lapse of patent right due to non-payment of the annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |