Invention content
The present invention provides a kind of sewage deep-sea based on power flow changing to solve technical problem present in known technology
Optimize discharge method, real-time optimization adjustment can be carried out with complicated Tidal Movement using the quantity of wastewater effluent of this method, make sewage
It is diluted in space-time unique inner height as small as possible, to promote its Initial dilution effect and final qualified discharge.
The present invention is adopted the technical scheme that solve technical problem present in known technology:One kind is become based on trend
The sewage deep-sea of change optimizes discharge method, and in Marine Outfall Diffuser or close position installs power flow monitor equipment;In blowdown
The response in the power flow changing section of engineering marine site between flow velocity and corresponding optimal quantity of wastewater effluent is established in control system to close
System:
(Qmin,Q1,…,Qi,…,Qn,Qmax)=F (umin,u1,…,ui,…,un,umax)
In formula, (umin,u1,…,ui,…,un,umax) be engineering marine site power flow changing section in flow velocity, pass through collect work
Journey and Adjacent Sea Area trend historic measured data obtain;(Qmin,Q1,…,Qi,…,Qn,Qmax) be and tidal current speed pair in section
The optimal quantity of wastewater effluent answered;F is response relation between the two;
Blowdown control system receives the trend reality that power flow monitor equipment monitors obtained Sewage outlet near zone in real time
Velocity measurement u ' data, and be input in response relation, the real-time optimal quantity of wastewater effluent for trend actual measurement flow velocity u ' is obtained, and
According to this, optimal quantity of wastewater effluent sends quantity of wastewater effluent instruction to sewage pump in real time;
Optimal quantity of wastewater effluent in response relation is calculated using simulation anti-inference method to be obtained, and is as follows:
1) simulation reverse calculation algorithms model is established
1.1) sewage draining exit pollutant is used to spread mathematical model;
1.2) multiple objective function is established
MAX.DD=C0/C1 (1)
MIN.EC=0.0414 × Q+13.828 (2)
MAX.OE=OE (3)
In formula (1), DD is pollutant dilution degree;C0For pollutant concentration value in jet stream water body;C1For certain point in environment water
Locate pollutant concentration value;EC is energy consumption value;Q is pump quantity of wastewater effluent, and Q=m × v × S, m are spout number, and v is discharge
Rate, S are spout cross-sectional area;OE is sewage discharge efficiency;
1.3) Genetic Algorithm Model is used
1.3.1 evolutionary computation operator) is used
1.3.2 calculated crosswise operator) is used
1.3.3 operator) is calculated using variation
2) initial setting up
2.1) using quantity of wastewater effluent as the control variable of optimization process, and real coding form is used;
2.2) setting control variable-value range, and in this, as the constraints of Population in Genetic Algorithms value;
3) each tidal current speed being directed in n+2 tidal current speed follows the steps below optimal quantity of wastewater effluent
Simulate inverse, you can obtain the n+2 optimal quantity of wastewater effluent for n+2 tidal current speed
3.1) it is directed to each tidal current speed, r numerical value is randomly selected within the scope of control variable-value, obtains r items dye
Colour solid constitutes the population (r) that the optimal quantity of wastewater effluent under the conditions of the tidal current speed may solve;
3.2) using each chromosome as the numerical example, r calculating is carried out using sewage draining exit pollutant diffusion mathematical model,
And r pollutant concentration result for calculating waters internal observation point is exported and stored simultaneously;
3.3) using best dilution of sewage effect, lowest energy consumption and optimal sewage discharge efficiency as target, using step
1.2) establish multiple objective function carry out evaluated chromosome;
3.4) it when meeting maximum generation number target, jumps out program and exports the i.e. optimal quantity of wastewater effluent of optimum dyeing body;
If being unsatisfactory for maximum generation number target, genetic computation is carried out using the Genetic Algorithm Model in step 1.3);
3.5) by genetic computation, selected chromosome in step 3.1) is updated, and obtain completely new population, according to step
(3.2)-(3.4) are computed repeatedly, until exporting optimum dyeing body, that is, optimal quantity of wastewater effluent.
The method for building up of the response relation uses following steps:
1) the BP Algorithm model of data-driven is established, it is non-linear in BP Algorithm model
Transfer function uses hyperbolic type Sigmoid functions, as shown in formula (9):
F (x)=1/1+e-x+θ (9)
θ indicates threshold values in formula;
Using the predictive ability of root-mean-square error (RMSE) evaluation network in BP Algorithm model, such as formula
(10) shown in:
A is data amount check in formula;YiFor measured value;Y′iFor predicted value;
The normalized of all data uses formula (11) in BP Algorithm model:
Wherein:Indicate the value after normalization;YiIndicate the value before normalization;YmaxIndicate the maximum value in all data;
YminIndicate the minimum value in all data;
2) network weight is initialized, in conjunction with n+2 tidal current speed and corresponding optimal quantity of wastewater effluent setting network
Framework, and the basic parameter of BP Algorithm model is set;
3) using n+2 tidal current speed data as the input of BP Algorithm model, accordingly most by n+2
Output of the excellent quantity of wastewater effluent as BP Algorithm model, is updated in artificial nerve network model and is trained
Study calculates root-mean-square error;
4) it when meeting root-mean-square error (RMSE) error requirements, jumps out program and establishes response relation;If being unsatisfactory for square
Root error (RMSE) error requirements then enter follow-up calculate;
5) using root-mean-square error as foundation, after reverse adjustment network weight, according still further to step 2) -4) it is computed repeatedly,
Until establishing response relation.
The invention has the advantages and positive effects that:By being consumed with best dilution of sewage effect, lowest energy and most
Excellent sewage discharge efficiency is target, establishes flow velocity and corresponding optimal quantity of wastewater effluent in the power flow changing section of engineering marine site
Between response relation, using the response relation, using Observed current data flow velocity as foundation, real-time optimization adjust sewage deep sea emission amount,
Sewage is set to be diluted in space-time unique inner height as small as possible, it, can to promote its Initial dilution effect and final qualified discharge
The management level of lifting sewage deep sea emission safeguards that marine environment ecological quality is horizontal, engineering and Adjacent Sea Area environment is given birth to
State quality has certain protective effect and practical significance.And the present invention also has at low cost, good economy performance, does not need largely
Engineering construction cost, existing engineered fusion degree height and environmental ecological benefit outstanding feature.
Specific implementation mode
In order to further understand the content, features and effects of the present invention, the following examples are hereby given, and coordinate attached drawing
Detailed description are as follows:
Referring to Fig. 1, a kind of sewage deep sea emission system, is equipped with blowdown control system 2 and sewage pump in sewage pumping station 1
3, blow-off line 4 is connected on sewage pump 3, the discharge section of blow-off line 4 is located at sea level 7m hereinafter, in blow-off line 4
Discharge section is equipped with multiple Marine Outfall Diffusers 5, and spout is equipped in Marine Outfall Diffuser 5, and spout is horizontally disposed.
Above-mentioned exhaust system optimizes discharge method using the sewage deep-sea based on power flow changing, based on the trend of region,
Using sewage deep-sea row's Haikou region power flow monitor as foundation, using the blowdown control system in land sewage pumping station as main tool,
Blowdown control system 2 calculates optimal quantity of wastewater effluent in real time, and sends work order to sewage pump 3, controls the work of sewage pump
Efficiency, concrete measure are:In Marine Outfall Diffuser 5 or close position installs power flow monitor equipment 6;In blowdown control system
The response relation between flow velocity and corresponding optimal quantity of wastewater effluent in the power flow changing section of engineering marine site is established in 2:
(Qmin,Q1,…,Qi,…,Qn,Qmax)=F (umin,u1,…,ui,…,un,umax)
In formula, (umin,u1,…,ui,…,un,umax) be engineering marine site power flow changing section in flow velocity, pass through collect work
Journey and Adjacent Sea Area trend historic measured data obtain;(Qmin,Q1,…,Qi,…,Qn,Qmax) be and tidal current speed pair in section
The optimal quantity of wastewater effluent answered;F is response relation between the two.
Blowdown control system 2 receives the trend that power flow monitor equipment 6 monitors obtained Sewage outlet near zone in real time
Flow velocity u ' data are surveyed, and are input in above-mentioned response relation, the real-time optimal sewage row for trend actual measurement flow velocity u ' is obtained
High-volume, and according to this optimal quantity of wastewater effluent sends quantity of wastewater effluent instruction to sewage pump 3 in real time.
Before power flow monitor equipment 6 is installed, need to collect engineering and Adjacent Sea Area trend historic measured data, in conjunction with work
Journey marine site tidal current dynamics numerical simulation technology, the essential characteristics such as engineering marine site tidal currenttype, Flow Field Distribution are grasped in analysis, and are extracted
Tidal current speed maxima and minima determines engineering marine site power flow changing section and relevant location information.According to engineering marine site tide
The essential characteristics such as stream type, Flow Field Distribution and tidal current speed maxima and minima relevant location information, according to regional representativeness
By force, it constructs the basic principles such as simple, reliable and stable, in Marine Outfall Diffuser or close position, chooses suitable point, as
Power flow monitor equipment 6 is mounted in power flow monitor point by power flow monitor point.
Optimal quantity of wastewater effluent in above-mentioned response relation is calculated using simulation anti-inference method to be obtained, and is as follows:
1) simulation reverse calculation algorithms model is established
1.1) it uses sewage draining exit pollutant to spread mathematical model, please refers to following two documents, one) Zhou Feng, Liang Shuxiu, grandson
Jet stream spray angle influences numerical simulation [J] Journal of Dalian University of Technology Total, 2007,47 (4) to jet characteristics in clear morning circumstance of flowing water:
583-588. bis-) Li M C, Si Q, Liang S X, et al.Multiple Objectives for Genetically
Optimized Coupled Inversion Method for Jet Models in Flowing Ambient Fluid
[J].Engineering Applications of Computational Fluid Mechanics,2014,8(1):82-
90。
1.2) multiple objective function is established
MAX.DD=C0/C1 (1)
MIN.EC=0.0414 × Q+13.828 (2)
MAX.OE=OE (3)
In formula (1), DD is pollutant dilution degree;C0For pollutant concentration value in jet stream water body;C1For certain point in environment water
Locate pollutant concentration value;EC is energy consumption value;Q is pump quantity of wastewater effluent, and Q=m × v × S, m are spout number, and v is discharge
Rate, S are spout cross-sectional area;OE is sewage discharge efficiency, and for Practical Project person, the higher the better for efficiency;Formula (2) root
According to flow and energy consumption of the horizontal sewage pump of German KSB companies K400-500 types under different operational modes, using least square method
Fitting obtains.
1.3) Genetic Algorithm Model is used
The genetic algorithm application prior art, including three calculating sections respectively evolve, intersect and make a variation, specifically lose
Propagation algorithm model includes:
1.3.1 evolutionary computation operator) is used
1.3.2 calculated crosswise operator) is used
1.3.3 operator) is calculated using variation
Genetic algorithm includes:Evolutionary computation (roulette selection), calculated crosswise (as shown in formula (4) and formula (5)) and variation
It calculates (as shown in formula (6)-formula (8)).
Os1=ω × Pa1+ (1- ω) × Pa2 (4)
Os2=ω × Pa2+ (1- ω) × Pa1 (5)
In formula (4) and formula (5), ω is the random parameter between -0.25 to 1.25;Os and Pa respectively represent filial generation and parent;1
It is the label of filial generation and parent with 2;
Variation calculates:
X=X'+ Δs (t, y) (6)
In formula (6) and formula (7), X is mutant gene, and X' is initial gene, and t is current generation number, and T is maximum generation number, r
Random number between 0-1, b are systematic parameter.For parameter y according to formula (8) value, LD is gene minimum value, and UD is that gene is maximum
Value, i is the random number for taking 0 or 1 respectively.
2) initial setting up
2.1) using quantity of wastewater effluent as the control variable of optimization process, and real coding form is used.
2.2) setting control variable-value range, and in this, as the constraints of Population in Genetic Algorithms value.
3) each tidal current speed being directed in n+2 tidal current speed follows the steps below optimal quantity of wastewater effluent
Simulate inverse, you can the n+2 optimal quantity of wastewater effluent for n+2 tidal current speed are obtained, the specific steps are:
3.1) it is directed to each tidal current speed, r numerical value is randomly selected within the scope of control variable-value, obtains r items dye
Colour solid constitutes the population (r) that the optimal quantity of wastewater effluent under the conditions of the tidal current speed may solve.
3.2) using each chromosome as the numerical example, r calculating is carried out using sewage draining exit pollutant diffusion mathematical model,
And r pollutant concentration result for calculating waters internal observation point is exported and stored simultaneously.
3.3) using best dilution of sewage effect, lowest energy consumption and optimal sewage discharge efficiency as target, using step
1.2) multiple objective function established carries out evaluated chromosome.
3.4) it when meeting maximum generation number target, jumps out program and exports the i.e. optimal quantity of wastewater effluent of optimum dyeing body;
If being unsatisfactory for maximum generation number target, genetic computation is carried out using the Genetic Algorithm Model in step 1.3).
3.5) by genetic computation, selected chromosome in step 3.1) is updated, and obtain completely new population, according to step
(3.2)-(3.4) are computed repeatedly, until exporting optimum dyeing body, that is, optimal quantity of wastewater effluent.
Sewage draining exit pollutant is spread mathematical model (being also known as Jet model) to above-mentioned simulation anti-inference method and the operation energy disappears
Consumption formula is embedded into Genetic Algorithm Model, using best dilution of sewage effect, lowest energy consumption and optimum discharge efficiency as mesh
Mark carries out the coupling process model of optimal quantity of wastewater effluent simulation inverse, is carried out using sewage draining exit pollutant diffusion mathematical model
The calculating of each chromosome simultaneously exports calculating waters internal observation point pollution object concentration data;It is each that analysis is calculated using multiple objective function
Dilution, sewage transport energy consumption and the discharge efficiency of a chromosome;Using the maximum generation number of model specification as Rule of judgment,
Meet condition and jump out calculation procedure, is unsatisfactory for then entering genetic computation;Genetic computation is with three heredity of evolving, intersect and make a variation
It calculates operator and adjusts chromosome, form completely new population, then substitute into sewage draining exit pollutant diffusion mathematical model, carry out next generation
Calculating.
The method for building up of above-mentioned response relation uses following steps:
1) the BP Algorithm model of data-driven is established, it is non-linear in BP Algorithm model
Transfer function uses hyperbolic type Sigmoid functions, as shown in formula (9):
F (x)=1/1+e-x+θ (9)
θ indicates threshold values in formula;
Using the predictive ability of root-mean-square error (RMSE) evaluation network in BP Algorithm model, such as formula
(10) shown in:
A is data amount check in formula;YiFor measured value;Yi' it is predicted value.
The normalized of all data uses formula (11) in BP Algorithm model:
Wherein:Indicate the value after normalization;YiIndicate the value before normalization;YmaxIndicate the maximum value in all data;
YminIndicate the minimum value in all data.
2) network weight is initialized, in conjunction with n+2 tidal current speed and corresponding optimal quantity of wastewater effluent setting network
Framework, and the basic parameter of BP Algorithm model is set.The determination of basic parameter can refer to document:Li Mingchang,
Zhang Guangyu, department's fine jade wait to couple optimization with the marine site assembled unit water quality model multi-parameter substep of genetic computation based on data-driven
Inversion method studies the practice and understanding of [J] mathematics, 2015,45 (12):167-175.
3) using n+2 tidal current speed data as the input of BP Algorithm model, accordingly most by n+2
Output of the excellent quantity of wastewater effluent as BP Algorithm model, is updated in artificial nerve network model and is trained
Study calculates root-mean-square error.
4) it when meeting root-mean-square error (RMSE) error requirements, jumps out program and establishes response relation;If being unsatisfactory for square
Root error (RMSE) error requirements then enter follow-up calculate.
5) using root-mean-square error as foundation, after reverse adjustment network weight, according still further to step 2) -4) it is computed repeatedly,
Until establishing response relation.
The advantage of above-mentioned response relation method for building up is, can using data-driven model-BP Algorithm
With on the basis of knowing little about it to system physical knowledge, only using system state variables and control variable as mode input, defeated
Go out, and the characteristics of analysis system data, you can establish the response relation between system variable.
Although the preferred embodiment of the present invention is described above in conjunction with attached drawing, the invention is not limited in upper
The specific implementation mode stated, the above mentioned embodiment is only schematical, be not it is restrictive, this field it is common
Technical staff under the inspiration of the present invention, in the case where not departing from present inventive concept and scope of the claimed protection, goes back
Many forms can be made, within these are all belonged to the scope of protection of the present invention.