CN108287972B - Method for calculating deep sea optimized discharge amount of sewage based on tidal current change - Google Patents
Method for calculating deep sea optimized discharge amount of sewage based on tidal current change Download PDFInfo
- Publication number
- CN108287972B CN108287972B CN201810101666.4A CN201810101666A CN108287972B CN 108287972 B CN108287972 B CN 108287972B CN 201810101666 A CN201810101666 A CN 201810101666A CN 108287972 B CN108287972 B CN 108287972B
- Authority
- CN
- China
- Prior art keywords
- sewage
- optimal
- discharge amount
- calculation
- sewage discharge
- 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
- 239000010865 sewage Substances 0.000 title claims abstract description 76
- 230000008859 change Effects 0.000 title claims abstract description 14
- 238000000034 method Methods 0.000 title claims abstract description 14
- 238000004364 calculation method Methods 0.000 claims abstract description 34
- 239000003344 environmental pollutant Substances 0.000 claims abstract description 25
- 231100000719 pollutant Toxicity 0.000 claims abstract description 25
- 210000000349 chromosome Anatomy 0.000 claims abstract description 23
- 230000002068 genetic effect Effects 0.000 claims abstract description 23
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 19
- 239000012895 dilution Substances 0.000 claims abstract description 16
- 238000010790 dilution Methods 0.000 claims abstract description 16
- 238000005265 energy consumption Methods 0.000 claims abstract description 14
- 238000004088 simulation Methods 0.000 claims abstract description 14
- 230000000694 effects Effects 0.000 claims abstract description 13
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 13
- 238000009792 diffusion process Methods 0.000 claims abstract description 11
- 238000013178 mathematical model Methods 0.000 claims abstract description 11
- 230000009191 jumping Effects 0.000 claims abstract description 5
- 230000007613 environmental effect Effects 0.000 claims description 4
- 238000011156 evaluation Methods 0.000 claims description 3
- 230000035772 mutation Effects 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 239000000243 solution Substances 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 abstract description 2
- 108090000623 proteins and genes Proteins 0.000 description 4
- 238000011161 development Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
Abstract
The invention discloses a method for calculating deep sea optimized discharge amount of sewage based on tidal current change, which is characterized in that a simulation inverse algorithm is adopted, a sewage outlet pollutant diffusion mathematical model and an operation energy consumption formula are embedded into a genetic algorithm model, a coupling method model for performing simulation inverse calculation on optimal sewage discharge amount is used as targets of optimal sewage dilution effect, lowest energy consumption and optimal discharge efficiency, the sewage outlet pollutant diffusion mathematical model is adopted to calculate each chromosome and output pollutant concentration data of an observation point in a calculated water area; calculating and analyzing the dilution, sewage conveying energy consumption and discharge efficiency of each chromosome by adopting a multi-objective function; and (4) taking the maximum generation number set by the model as a judgment condition, jumping out of a calculation program when the condition is met, and entering genetic calculation when the condition is not met. The method determines the sewage discharge amount, and can highly dilute the sewage within a space-time range as small as possible so as to improve the initial dilution effect and finally discharge the sewage up to the standard.
Description
Technical Field
The invention relates to a sewage discharge method, in particular to a deep sea optimized discharge amount calculation method of sewage based on tidal current change.
Background
In recent years, with the rapid development of coastal economy in China, a series of major maritime engineering is frequently built, and newly added pollution discharge inevitably causes greater ecological environment pressure on engineering sea areas. By means of the extremely strong mixing and transporting capacity of waves, tides and currents, the offshore deep sea discharge of sewage becomes an engineering measure of priority for the sewage discharge in coastal areas due to the low engineering cost and the low operation cost. However, the current deep sea discharge engineering cannot realize a mechanism for optimizing dynamic regulation in real time along with the motion of complex tidal waves, so that the sewage entering the environmental water body cannot obtain the optimal real-time dilution effect. Therefore, research and development of a deep-sea discharge port sewage optimal discharge amount calculation method based on tidal current change are urgently needed, and the method is used for guiding sewage discharge, so that sewage is highly diluted in a space-time range as small as possible, the initial dilution effect of the sewage is improved, and the sewage finally reaches the standard for discharge.
Disclosure of Invention
The invention provides a method for calculating deep sea optimized discharge amount of sewage based on tidal current change, which is used for determining the discharge amount of sewage, so that the sewage can be highly diluted in a space-time range as small as possible, the initial dilution effect is improved, and the sewage finally reaches the standard and is discharged.
The technical scheme adopted by the invention for solving the technical problems in the prior art is as follows: a method for calculating deep sea optimized discharge of sewage based on tidal current change adopts a simulation inverse algorithm and comprises the following specific steps:
1) establishing a simulation back calculation method model
1.1) adopting a sewage outlet pollutant diffusion mathematical model;
1.2) establishing a Multi-objective function
MAX.DD＝C0/C1 (1)
MIN.EC＝0.0414×Q+13.828 (2)
MAX.OE＝OE (3)
In the formula (1), DD is the dilution of pollutants; c0The concentration value of the pollutants in the jet flow water body is obtained; c1The concentration value of the pollutant at a certain point in the environmental water body is obtained; EC is energy consumption value; q is the sewage discharge amount of the pump, m is m multiplied by v multiplied by S, m is the number of the nozzles, v is the discharge rate, and S is the sectional area of the nozzles; OE is the sewage discharge efficiency;
1.3) Using a genetic Algorithm model
1.3.1) calculation operator by evolution
1.3.2) using cross-computing operators
1.3.3) calculation operator by mutation
2) Initial setting
2.1) taking the sewage discharge as a control variable of an optimization process, and adopting a real number coding form;
2.2) setting a value range of the control variable, and taking the value range as a constraint condition of the genetic algorithm population value;
3) aiming at the tidal current change interval (u) of the engineering sea areamin,u1,…,ui,…,un,umax) According to each tidal flow velocityThe simulation inverse calculation of the optimal sewage discharge amount is carried out according to the following steps, and the n +2 optimal sewage discharge amounts (Q) aiming at the n +2 tidal current flow velocities can be obtainedmin,Q1,…,Qi,…,Qn,Qmax)
3.1) randomly selecting r numerical values in the control variable value range according to each tidal current flow rate to obtain r chromosomes, and forming a population (r) of the possible solution of the optimal sewage discharge amount under the tidal current flow rate condition;
3.2) taking each chromosome as a numerical value example, adopting a pollutant diffusion mathematical model to calculate r times, and simultaneously outputting and storing r pollutant concentration results of observation points in the water area;
3.3) carrying out chromosome evaluation by adopting the multi-target function established in the step 1.2) with the target of optimal sewage dilution effect, lowest energy consumption and optimal sewage discharge efficiency;
3.4) when the maximum generation number is met, jumping out of the program and outputting the optimal chromosome, namely the optimal sewage discharge amount; if the maximum generation number target is not met, adopting the genetic algorithm model in the step 1.3) to perform genetic calculation;
3.5) updating the chromosomes selected in the step 3.1) through genetic calculation, obtaining a brand new population, and repeating the calculation according to the steps (3.2) to (3.4) until the best chromosomes, namely the optimal sewage discharge amount, is output.
The invention has the advantages and positive effects that: the optimal sewage discharge amount corresponding to the flow speed in the tidal current change interval of the engineering sea area is calculated by taking the optimal sewage dilution effect, the lowest energy consumption and the optimal sewage discharge efficiency as targets, the optimal sewage discharge amount is used for optimally adjusting the deep sea discharge amount of the sewage, the sewage is highly diluted in a space-time range as small as possible, the initial dilution effect is improved, the sewage is finally discharged up to the standard, the management level of the deep sea discharge of the sewage can be improved, the ecological quality level of the sea area environment is maintained, and a certain protection effect and practical significance are achieved for the ecological quality of the engineering and the adjacent sea area environment.
Detailed Description
To further understand the contents, features and effects of the present invention, the following examples are illustrated as follows:
a method for calculating deep sea optimized discharge of sewage based on tidal current change adopts a simulation inverse algorithm and comprises the following specific steps:
1) establishing a simulation back calculation method model
1.1) the mathematical model of the diffusion of pollutants at sewage outlets is adopted, please refer to the following two documents, namely, first) the numerical simulation of the influence of the Jet angle on the Jet characteristics in the environment of Flowing water [ J ]. university of continental Engineering, 2007,47(4): 583-.
1.2) establishing a Multi-objective function
MAX.DD＝C0/C1 (1)
MIN.EC＝0.0414×Q+13.828 (2)
MAX.OE＝OE (3)
In the formula (1), DD is the dilution of pollutants; c0The concentration value of the pollutants in the jet flow water body is obtained; c1The concentration value of the pollutant at a certain point in the environmental water body is obtained; EC is energy consumption value; q is the sewage discharge amount of the pump, m is m multiplied by v multiplied by S, m is the number of the nozzles, v is the discharge rate, and S is the sectional area of the nozzles; OE is sewage discharge efficiency, and for practical engineers, the higher the efficiency, the better; the formula (2) is obtained by adopting least square fitting according to the flow and the energy consumption of the Germany KSB company K400-500 type horizontal sewage pump under different operation modes.
1.3) Using a genetic Algorithm model
The genetic algorithm applies the prior art, comprises three calculation parts, namely evolution, intersection and variation, and specifically comprises a genetic algorithm model:
1.3.1) calculation operator by evolution
1.3.2) using cross-computing operators
1.3.3) calculation operator by mutation
The genetic algorithm comprises the following steps: evolutionary calculations (roulette selection), cross calculations (as shown in equations (4) and (5)), and variant calculations (as shown in equations (6) to (8)).
Os1＝ω×Pa1+(1-ω)×Pa2 (4)
Os2＝ω×Pa2+(1-ω)×Pa1 (5)
In the formulas (4) and (5), omega is a random parameter between-0.25 and 1.25; os and Pa represent offspring and father, respectively; 1 and 2 are labels for children and parents;
and (3) calculating variation:
X＝X'+Δ(t,y) (6)
in the formulas (6) and (7), X is a mutated gene, X' is an initial gene, T is the current generation number, T is the maximum generation number, r is a random number between 0 and 1, and b is a system parameter. The parameter y takes the value according to the formula (8), LD is the minimum value of the gene, UD is the maximum value of the gene, and i is a random number which is 0 or 1 respectively.
2) Initial setting
2.1) taking the sewage discharge amount as a control variable of an optimization process, and adopting a real number coding form.
And 2.2) setting the value range of the control variable, and taking the value range as a constraint condition of the population value of the genetic algorithm.
3) Aiming at the tidal current change interval (u) of the engineering sea areamin,u1,…,ui,…,un,umax) The simulation reverse calculation of the optimal sewage discharge amount is carried out according to the following steps for each tidal current flow velocity, and then n +2 optimal sewage discharge amounts (Q) aiming at n +2 tidal current flow velocities can be obtainedmin,Q1,…,Qi,…,Qn,Qmax)
3.1) randomly selecting r numerical values in the control variable value range according to each tidal current flow velocity to obtain r chromosomes, and forming a population (r) of the possible solution of the optimal sewage discharge amount under the tidal current flow velocity condition.
And 3.2) taking each chromosome as a numerical example, adopting a pollutant diffusion mathematical model to perform r times of calculation, and simultaneously outputting and storing r pollutant concentration results of observation points in the water area.
3.3) carrying out chromosome evaluation by adopting the multi-objective function established in the step 1.2) with the objectives of optimal sewage dilution effect, lowest energy consumption and optimal sewage discharge efficiency.
3.4) when the maximum generation number is met, jumping out of the program and outputting the optimal chromosome, namely the optimal sewage discharge amount; if the maximum generation number target is not met, adopting the genetic algorithm model in the step 1.3) to perform genetic calculation.
3.5) updating the chromosomes selected in the step 3.1) through genetic calculation, obtaining a brand new population, and repeating the calculation according to the steps (3.2) to (3.4) until the best chromosomes, namely the optimal sewage discharge amount, is output.
The simulation inverse algorithm embeds a sewage discharge port pollutant diffusion mathematical model (also called a jet flow model) and an operation energy consumption formula into a genetic algorithm model, takes the optimal sewage dilution effect, the lowest energy consumption and the optimal discharge efficiency as targets, carries out a coupling method model of optimal sewage discharge amount simulation inverse calculation, adopts the sewage discharge port pollutant diffusion mathematical model to calculate each chromosome and outputs pollutant concentration data of an observation point in a calculated water area; calculating and analyzing the dilution, sewage conveying energy consumption and discharge efficiency of each chromosome by adopting a multi-objective function; taking the maximum generation number set by the model as a judgment condition, jumping out of a calculation program when the condition is met, and entering genetic calculation when the condition is not met; and (4) adjusting chromosomes by genetic calculation operators of evolution, intersection and variation to form a brand new population, and substituting the brand new population into a pollutant diffusion mathematical model of the sewage outlet to perform calculation of the next generation.
(umin,u1,…,ui,…,un,umax) The flow velocity in the engineering sea area tidal current change interval is obtained by collecting historical actual measurement data of engineering and adjacent sea areas tidal current.
Although the preferred embodiments of the present invention have been described, the present invention is not limited to the above-mentioned embodiments, which are only illustrative and not restrictive, and those skilled in the art can make many modifications without departing from the spirit and scope of the present invention as claimed in the appended claims.
Claims (1)
1. A method for calculating deep sea optimized discharge amount of sewage based on tidal current change is characterized in that a simulation inverse algorithm is adopted, and the method comprises the following specific steps:
1) establishing a simulation back calculation method model
1.1) adopting a sewage outlet pollutant diffusion mathematical model;
1.2) establishing a Multi-objective function
MAX.DD＝C0/C1 (1)
MIN.EC＝0.0414×Q+13.828 (2)
MAX.OE＝OE (3)
In the formula (1), DD is the dilution of pollutants; c0The concentration value of the pollutants in the jet flow water body is obtained; c1The concentration value of the pollutant at a certain point in the environmental water body is obtained; EC is energy consumption value; q is the sewage discharge amount of the pump, m is m multiplied by v multiplied by S, m is the number of the nozzles, v is the discharge rate, and S is the sectional area of the nozzles; OE is the sewage discharge efficiency;
1.3) Using a genetic Algorithm model
1.3.1) calculation operator by evolution
1.3.2) using cross-computing operators
1.3.3) calculation operator by mutation
2) Initial setting
2.1) taking the sewage discharge as a control variable of an optimization process, and adopting a real number coding form;
2.2) setting a value range of the control variable, and taking the value range as a constraint condition of the genetic algorithm population value;
3) aiming at engineering sea areaTidal current change interval (u)min,u1,…,ui,…,un,umax) The simulation reverse calculation of the optimal sewage discharge amount is carried out according to the following steps for each tidal current flow velocity, and then n +2 optimal sewage discharge amounts (Q) aiming at n +2 tidal current flow velocities can be obtainedmin,Q1,…,Qi,…,Qn,Qmax)
3.1) randomly selecting r numerical values in the control variable value range according to each tidal current flow rate to obtain r chromosomes, and forming a population (r) of the possible solution of the optimal sewage discharge amount under the tidal current flow rate condition;
3.2) taking each chromosome as a numerical value example, adopting a pollutant diffusion mathematical model to calculate r times, and simultaneously outputting and storing r pollutant concentration results of observation points in the water area;
3.3) carrying out chromosome evaluation by adopting the multi-target function established in the step 1.2) with the target of optimal sewage dilution effect, lowest energy consumption and optimal sewage discharge efficiency;
3.4) when the maximum generation number is met, jumping out of the program and outputting the optimal chromosome, namely the optimal sewage discharge amount; if the maximum generation number target is not met, adopting the genetic algorithm model in the step 1.3) to perform genetic calculation;
3.5) updating the chromosomes selected in the step 3.1) through genetic calculation, obtaining a brand new population, and repeating the calculation according to the steps (3.2) to (3.4) until the best chromosomes, namely the optimal sewage discharge amount, is output.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810101666.4A CN108287972B (en) | 2018-02-01 | 2018-02-01 | Method for calculating deep sea optimized discharge amount of sewage based on tidal current change |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810101666.4A CN108287972B (en) | 2018-02-01 | 2018-02-01 | Method for calculating deep sea optimized discharge amount of sewage based on tidal current change |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108287972A CN108287972A (en) | 2018-07-17 |
| CN108287972B true CN108287972B (en) | 2021-07-27 |
Family
ID=62836319
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810101666.4A Expired - Fee Related CN108287972B (en) | 2018-02-01 | 2018-02-01 | Method for calculating deep sea optimized discharge amount of sewage based on tidal current change |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108287972B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116341192A (en) * | 2023-01-17 | 2023-06-27 | 山东省海洋资源与环境研究院(山东省海洋环境监测中心、山东省水产品质量检验中心) | Sea area pollutant transport flux calculation method, system, equipment and medium |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105224772A (en) * | 2015-11-04 | 2016-01-06 | 松辽流域水资源保护局松辽流域水环境监测中心 | A kind of Nenjiang representative region Ecology methods of risk assessment based on system dynamics decision model |
| CN105843072A (en) * | 2016-03-22 | 2016-08-10 | 华南理工大学 | Wastewater treatment energy conservation optimization method based on improved local search immune genetic algorithm |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7500437B2 (en) * | 2004-08-27 | 2009-03-10 | Neuco, Inc. | Method and system for SCR optimization |
-
2018
- 2018-02-01 CN CN201810101666.4A patent/CN108287972B/en not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105224772A (en) * | 2015-11-04 | 2016-01-06 | 松辽流域水资源保护局松辽流域水环境监测中心 | A kind of Nenjiang representative region Ecology methods of risk assessment based on system dynamics decision model |
| CN105843072A (en) * | 2016-03-22 | 2016-08-10 | 华南理工大学 | Wastewater treatment energy conservation optimization method based on improved local search immune genetic algorithm |
Non-Patent Citations (1)
| Title |
|---|
| 动水射流多参数多目标非线性遗传优化耦合反演方法研究;李明昌等;《水力发电学报》;20140831;第33卷(第4期);第92-96页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108287972A (en) | 2018-07-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110160526B (en) | Air route planning method based on genetic algorithm | |
| CN108805434B (en) | Improved NSGA-III-based multi-objective optimization scheduling method for cascade hydropower station | |
| CN107506909B (en) | Cascade reservoir hydropower station group scheduling control system and method for fish habitat protection | |
| CN111126847B (en) | Cascade reservoir short-term optimization scheduling method and system coupled with riverway water power process | |
| CN105911867B (en) | Ship thrust distribution method based on NSGA-II algorithm | |
| CN111005346B (en) | Reservoir group multi-objective action mechanism and optimization scheduling scheme analysis method | |
| CN109948235B (en) | Water resource scheduling and accurate configuration method | |
| CN111476676B (en) | Classification, classification and classification intelligent scheduling method for small river basin in town | |
| CN110363343B (en) | Hybrid self-adaptive hydropower station group intelligent optimization scheduling method and system | |
| CN115375198B (en) | Method and system for communication joint scheduling and water quality safety guarantee of regional river and lake water systems | |
| CN105045091B (en) | Dredging technique intelligent decision analysis method based on Fuzzy Neural Control system | |
| CN103761583A (en) | Reservoir sediment ejection power generation multi-target optimized dispatching method | |
| Song et al. | Basin-scale multi-objective simulation-optimization modeling for conjunctive use of surface water and groundwater in northwest China | |
| Azamathulla et al. | Prediction of scour depth at culvert outlets using gene-expression programming | |
| CN103886141A (en) | Method for computing fluvial cross section information in automatic-batch mode | |
| CN108287972B (en) | Method for calculating deep sea optimized discharge amount of sewage based on tidal current change | |
| CN115758886A (en) | Optimal distribution method of storage tanks based on rain flood pipe network model and decision tree algorithm | |
| CN108304654B (en) | Sewage deep sea optimized discharge method based on tidal current change | |
| CN105512383A (en) | Dredging process regulation and control parameter screening method based on BP neural network | |
| CN105334824A (en) | Aluminum electrolysis production optimization method based on NSGA-II algorithm | |
| CN108280791B (en) | Sewage deep sea optimized discharge response relation determination method based on tidal current change | |
| CN115618480A (en) | Dredger energy consumption optimization method and system based on genetic algorithm | |
| CN107818220A (en) | Evaluation method based on dynamics of ecosystem collective model to estuarine environment capacity | |
| CN115600527A (en) | Reservoir operation state prediction analysis method based on reservoir environment data | |
| CN115730739A (en) | Urban waterlogging prediction method based on LSTM neural network |
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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20220316 Address after: No.37, Xingang No.2 Road, Binhai New Area, Tianjin, 300456 Patentee after: Tianjin Research Institute for Water Transport Engineering, M.O.T Patentee after: TIWTE ENVIRONMENTAL TECHNOLOGY DEVELOPMENT (TIANJIN) CO.,LTD. Address before: No.2618, No.2 Road, Tanggu, Binhai New Area, Tianjin Patentee before: Tianjin Research Institute for Water Transport Engineering, M.O.T |
|
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210727 |