CN101799846A - Multi-objective groundwater remediation optimization method - Google Patents
Multi-objective groundwater remediation optimization method Download PDFInfo
- Publication number
- CN101799846A CN101799846A CN201010018361A CN201010018361A CN101799846A CN 101799846 A CN101799846 A CN 101799846A CN 201010018361 A CN201010018361 A CN 201010018361A CN 201010018361 A CN201010018361 A CN 201010018361A CN 101799846 A CN101799846 A CN 101799846A
- Authority
- CN
- China
- Prior art keywords
- neighborhood
- subsolution
- separated
- pareto
- solution
- 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
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 37
- 239000003673 groundwater Substances 0.000 title abstract description 10
- 238000005067 remediation Methods 0.000 title abstract description 3
- 230000008569 process Effects 0.000 claims abstract description 12
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 33
- 238000004422 calculation algorithm Methods 0.000 claims description 25
- 230000006978 adaptation Effects 0.000 claims description 12
- 238000011160 research Methods 0.000 claims description 4
- 238000004088 simulation Methods 0.000 claims description 4
- 238000000151 deposition Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 4
- 238000003895 groundwater pollution Methods 0.000 abstract 1
- 239000000243 solution Substances 0.000 description 44
- HNBKNDFZJLEZLZ-UHFFFAOYSA-N 2-[methyl-[3-[4-(4-methylbenzoyl)phenoxy]-3-phenylpropyl]amino]acetic acid Chemical compound C=1C=CC=CC=1C(CCN(C)CC(O)=O)OC(C=C1)=CC=C1C(=O)C1=CC=C(C)C=C1 HNBKNDFZJLEZLZ-UHFFFAOYSA-N 0.000 description 10
- 238000011109 contamination Methods 0.000 description 10
- 238000013461 design Methods 0.000 description 10
- 230000002068 genetic effect Effects 0.000 description 9
- 230000008439 repair process Effects 0.000 description 6
- 238000012546 transfer Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 239000003344 environmental pollutant Substances 0.000 description 4
- 231100000719 pollutant Toxicity 0.000 description 4
- 238000011084 recovery Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 229910052757 nitrogen Inorganic materials 0.000 description 2
- 229910002651 NO3 Inorganic materials 0.000 description 1
- NHNBFGGVMKEFGY-UHFFFAOYSA-N Nitrate Chemical compound [O-][N+]([O-])=O NHNBFGGVMKEFGY-UHFFFAOYSA-N 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a multi-objective groundwater remediation optimization method. The individual fitness function value archiving strategy, the Pareto dominance concept and a niche fitness sharing method are combined with tabu search (TS) to systematically form a new multi-objective optimization method. The method integrates the advantages of various optimization technology and is high in computing efficiency and strong in global search capability. Compared with the similar methods, the method of the invention mainly avoids repeated search in the optimization process and improves the diversity and the global property of seed solutions. The technology is coupled with the groundwater flow program (MODFLOW) and the solute transport program (MT3DMS) and has strong applicability to solving the multi-objective management model for groundwater pollution control.
Description
Technical field
The present invention relates to the Hydrology and Water Resources field, specifically is that underground water is repaired Multipurpose Optimal Method.
Background technology
Groundwater contamination is repaired in the optimum management system, often relate to a plurality of conflicting management objectives, have some characteristics between these targets, the linear module that comprises each target can not commensurability, conflicting between each objective function, this makes separating of multi-objective optimization question not be unique, and the task of multiobjective management is the pros and cons of each target of balance, seeks to satisfy the compromise of each target interests and separates.
Simulative optimization model (simulation-optimization model) is to find the solution the groundwater contamination repair system to manage effective method, also is that purposes is found the solution instrument the most widely.Analogy model is used for bringing in constant renewal in state variable, comprises groundwater flow motion model and solute transfer model; Optimization Model is used for selecting the optimizing decision variable, finds the solution the balance of multi-objective problem and separates, the foundation of making a strategic decision as the supvr.Repair in the multiobjective management problem at underground water, often relate to a plurality of conflicting objective functions, wherein of paramount importance is the minimum economic cost of realizing the recovery well optimal design on the basis of various constraint conditions of sewerage contentedly, and the recovery well optimal design comprises that well location puts the optimization with well yield.In addition the quality of residual contamination thing has become an important indicator differentiating the optimal design success or not after the groundwater contamination reparation.
In recent years based on population and in finding the solution multi-objective optimization question, obtained using widely based on the Multiobjective Intelligent algorithm of Pareto, in underground water remediation management model, also obtained development rapidly.For example, compared Pareto genetic algorithm (Pareto GA), vector is estimated genetic algorithm (Vector Evaluated GA, VEGA), microhabitat multi-objective genetic algorithm (Niched Pareto Genetic Alogrithm, NPGA), improved non-domination ordering genetic algorithm (Improved Non-dominated Sorting Genetic Algorithm, ε-NSGAII), at random multi-objective genetic algorithm (Probabilistic multi-objective genetic algorithm, PMOGA).More than finding the solution underground water multiobjective management model all is the Multiobjective Intelligent algorithm of utilization based on GA.Yet the Multiobjective Intelligent algorithm based on tabu search does not use as yet in finding the solution underground water multiobjective management problem.
Summary of the invention
The invention provides a kind of underground water and repair the new method of multi-objective optimization question, this method has higher counting yield, has guaranteed the ability of searching optimum of algorithm simultaneously.
Underground water of the present invention is repaired Multipurpose Optimal Method, and it may further comprise the steps:
Step 3 coupled simulation model and Optimization Model;
The result is optimized in step 5 output, analyzes the of overall importance and counting yield that balance is separated, and adjusts parameters optimization.
Wherein, step 4 comprises following detailed process:
1) producing initial solution and corresponding neighborhood separates
At first, from initial solution, adopt random fashion to produce an initial solution to the search in feasible solution territory, it then is to be basic point with the initial solution that neighborhood is separated, and carries out neighborhood according to the neighborhood move function and moves the back and produce, and the neighborhood move function is defined as follows:
s
i=(2×random()-1)×stepsize,i=1,2,…,i
max????????????????????????(1)
S wherein
i(i=1 ..., i
Max) represent that neighborhood moves rule function, i
MaxBe the number of times that neighborhood moves, random () is the random function that produces real number between the 0-1, and stepsize is the neighborhood moving step length.Initial solution and neighborhood are separated to produce and are had randomness in this research;
2) calculate target function value, Pareto rank and the individual crowding that neighborhood is separated
After the generation neighborhood is separated, calculate the target function value that all neighborhoods are separated,, calculate Pareto rank and individual crowding value that each neighborhood is separated according to target function value; The Pareto rank number of separating equal neighborhood separate in all Pareto control the skill that this separates, the Pareto rank is 0 separates and be nondominated solution;
3) select kind of a subsolution
Kind of subsolution is the initial solution of searching in per generation, and the selection of planting subsolution is derived from the previous generation neighborhood and separates, and promptly finds to separate separating as candidate seed that Pareto is dominant with respect to other neighborhood and separate from neighborhood is separated, and must satisfy the uncontrolled property of Pareto;
4) upgrade Pareto table, candidate list and taboo table
Pareto table and candidate list are respectively applied for the kind subsolution of depositing each generation and do not have selected other candidate solution for kind of subsolution; To no longer have separating from Pareto table and candidate list of Pareto optimality in per generation search and remove, the individuality that is chosen as kind of subsolution in the candidate list will be removed from table; Upgrade the taboo table with kind of subsolution as the taboo object, selected kind subsolution is avoided in the taboo algebraically of regulation, and seed is lifted a ban the new taboo kind subsolution in back replacement is deposited in the taboo table;
5) judge whether to satisfy stopping criterion
In the optimizing process, if reached predefined maximum iteration time, seed selection subsolution collection also is empty for empty candidate list simultaneously when perhaps searching certain generation, then can't find the next subsolution of planting, algorithm can't enter the search in next stage, and algorithm stops, output Pareto optimal solution set.
Said process 2) adopt in ideal adaptation degree functional value archival strategy note first search to separate the information of separating in the direct call function value storehouse of separating then of repeat search.
The selection of planting subsolution said process 3) has following several situation: if a) only exist a candidate seed to separate, then this is separated as follow-on kind of subsolution; B) if existence is separated more than one candidate seed, contrast the individual crowding that candidate seed is separated, the individuality of choosing crowding value minimum is as follow-on kind of subsolution; C) if there is no candidate seed is separated, and then selects to enter separating as new kind subsolution in the table the earliest from candidate list, and this is separated from candidate list delete.
Optimization method of the present invention can be coupled with groundwater simulation model (MODFLOW+MT3DMS), is applied to optimal design underground water and repairs the multiple-objection optimization problem of management.The notion that the present invention is dominant according to Pareto is found the solution the balance of multi-objective optimization question and is separated, and has only the nondominated solution of neighborhood in separating just might win in competition in each iteration of search optimization solution.The mode of NPTS search is based on the neighborhood search of certain initial solution, and progressively to the process of global search, nondominated solution and be positioned at the neighborhood in loose zone and separate and be chosen as follow-on initial solution in each neighborhood search has strengthened the ability of algorithm global search greatly.
This method has been inherited tabu search search capability efficiently, guarantees the diversity of understanding simultaneously by taboo strategy and microhabitat fitness technology of sharing, has strengthened the ability of searching optimum of algorithm.NPTS repairs in the multiobjective management problem at underground water and has broad application prospects.
Description of drawings
Fig. 1 is the NPTS process flow diagram;
Fig. 2 is the Pareto ordering chart, and the Pareto diagram of sorting: the rank number equals to separate other number of separating that comprises in the square frame of lower left;
Fig. 3 is individual crowding comparison diagram, " △ " expression separate in the microhabitat radius is the round territory of radius population density greater than
Represented separates;
Fig. 4 administers place hydrogeological characteristics planimetric map;
Fig. 5 is that the Pareto that NPTS optimization obtains separates along the distribution plan of balance curve, and wherein horizontal ordinate is first objective function a---treatment cost, and ordinate is second objective function---residual contamination amount ratio." zero " is for satisfying the Pareto optimum solution of optimal design;
Fig. 6 is to use the step figure of ideal adaptation degree when strategy filing calculating target function number of times, the relation when ideal adaptation degree filing strategy is not used in straight line "-" expression that the figure top is parallel to transverse axis between objective function calculation times and the iterations.
Embodiment
1, the present invention is used for finding the solution underground water to repair the concrete steps of multiobjective management problem as follows:
Minimize J
2=MR=(mass
End/ mass
0) * 100% (3)
J wherein
1Expression treatment cost (RC), N
wFor administering the number of well, y
iBinary number (y during employing whether expression well i adopts
i=1, y when not adopting
i=0), Q
i tBe taking out in the t management phase/water injection rate (negative value is represented to draw water, on the occasion of the expression water filling) of well i, N
tBe total management issue, Δ t
tBe the duration of t stress phase, M
i tBe the quality of in the t management phase, removing pollutant of well i, α
i(i=1,2,3) expression well is installed, is taken out/water filling, the cost coefficient of pollutant control.J
2End of term residual contamination amount number percent (MR), mass are administered in expression
0And mass
EndExpression is administered under the original state and the pollutant quality in the improvement end of term respectively.
Constraint condition comprises administers the constraint of well sum, head constraint, hydraulic gradient constraint, pollutant levels constraint, individual well draw water traffic constraints and overall balance constraint.Objective function and constraint condition have constituted the mathematical model of underground water reparation optimum management problem.
Step 3 coupled simulation model and Optimization Model.The target of optimum management is the function (formula (2)) of decision variable (well yield) and state variable (head and solute concentration), need bring in constant renewal in state variable by analogy model, calculating target function value and judge whether to satisfy constraint condition; By Optimization Model trade-off decision variable, turn back to update mode variable in the analogy model simultaneously.Decision variable and state variable satisfy analogy model and state model simultaneously in the optimum management model, and upgrade synchronously.
Its Chinese style (4) and (5) expression system be t state variable h at a time
tAnd C
tBe this period decision variable Q
tAn and last period t-1 state variable h
T-1And C
T-1Function, state transition function trans
hAnd trans
CExpression.
The concrete steps of optimization method are (see figure 1):
1) producing initial solution and corresponding neighborhood separates
NPTS from initial solution, adopts random fashion to produce an initial solution to the search in feasible solution territory.It then is to be basic point with the initial solution that neighborhood is separated, and carries out neighborhood according to the neighborhood move function and moves that the back produces, and the neighborhood move function is defined as follows:
s
i=(2×random()-1)×stepsize,i=1,2,…,i
max????????????????????????(1)
S wherein
i(i=1 ..., i
Max) represent that neighborhood moves rule function, i
MaxBe the number of times that neighborhood moves, random () is the random function that produces real number between the 0-1, and stepsize is the neighborhood moving step length.Initial solution and neighborhood are separated to produce and are had randomness in this research.
2) calculate target function value, Pareto rank and the individual crowding that neighborhood is separated
After producing neighborhood and separating, calculate the target function value that all neighborhoods are separated, this method introduce ideal adaptation degree functional value archival strategy note first search to separate the information of separating in the direct call function value storehouse of separating then of repeat search; According to target function value, calculate Pareto rank and individual crowding value that each neighborhood is separated.The Pareto rank number of separating equal neighborhood separate in all Pareto control the skill that this separates, the Pareto rank is 0 separates and be nondominated solution (seeing Fig. 2,3).
3) select kind of a subsolution
Planting subsolution is the initial solution of searching in per generation, and the selection of planting subsolution is derived from the previous generation neighborhood and separates, and must satisfy the uncontrolled property of Pareto.Before determining kind of subsolution, need from neighborhood is separated, find and separate separating that Pareto is dominant with respect to other neighborhood and separate as candidate seed, the selection of planting subsolution can be summarized as following several situation: if a) only exist a candidate seed to separate, then this is separated as follow-on kind of subsolution; B) if existence is separated more than one candidate seed, contrast the individual crowding that candidate seed is separated, the individuality of choosing crowding value minimum is as follow-on kind of subsolution; C) if there is no candidate seed is separated, and then selects to enter separating as new kind subsolution in the table the earliest from candidate list, and this is separated from candidate list delete.
4) upgrade Pareto table, candidate list and taboo table
Pareto table and candidate list are respectively applied for the kind subsolution of depositing each generation and do not have selected other candidate solution for kind of subsolution.To no longer have separating from Pareto table and candidate list of Pareto optimality in per generation search and remove, the individuality that is chosen as kind of subsolution in the candidate list will be removed from table.Upgrade the taboo table with kind of subsolution as the taboo object, for fear of being absorbed in cyclic search, selected kind subsolution is avoided in the taboo algebraically of regulation, and seed is lifted a ban the new taboo kind subsolution in back replacement is deposited in the taboo table.
5) judge whether to satisfy stopping criterion
If reached predefined maximum iteration time in the optimizing process, seed selection subsolution collection then can't find the next subsolution of planting for empty candidate list simultaneously also is empty when perhaps searching certain generation, and algorithm can't enter the search in next stage, algorithm stops, output Pareto optimal solution set.
The result is optimized in step 5 output, analyzes the of overall importance and counting yield that balance is separated, and adjusts parameters optimization, makes optimization method can reach optimum optimization effect in different problems.
2, underground water of following design is repaired the multiobjective management problem, utilizes the present invention to find the solution to satisfy the balance of management objectives and constraint condition to separate.
2.1 problem overview
Study area is average, isotropic two dimension confined aquifer, has existed nitrate (in nitrogen) to pollute, and pollutes the distribution of plumage and sees Fig. 4.The long 3300m in water-bearing zone, wide 2400m, the square net that with the length of side is 150m is the finite difference grid of 17 row, 23 row with the survey region subdivision.The current boundary condition is: east and western part are respectively 25m and 35m for to decide head boundary, and northern and south is the water proof border.The solute transfer boundary condition is: western, northern and south is solute zero flux border, and the east is the given flux of solute convection border.Other relevant hydrogeological parameter, current and solute transfer parameter are seen Fig. 4.4 mouthfuls of preliminary elections are administered well and are distributed as shown in Figure 4.4 decision variables are set in the Optimization Model, are the flow that draws water of 4 mouthfuls of wells, and limit the traffic constraints (0≤Q that draws water of every mouthful of well
i≤ 10000m
3/ d, i=1,2 ..., 4).The management objectives of Optimization Model satisfy concentration constraint condition simultaneously for minimizing treatment cost (formula (2)) and minimizing the residual contamination amount than (formula (3)), and the concentration that promptly satisfies concentration constraint (shadow region among Fig. 4) internal contamination thing is less than 3ppb.The optimization improvement cycle is assumed to be 5 years.
2.2 application the present invention
Being provided with of correlation parameter of the present invention is as follows: calculating algebraically is 400; The number of times that neighborhood moves is i
Max=2; Neighborhood moving step length stepzize=3; The discretize interval number of each parameter is 32; The microhabitat radius is 0.05; The length of taboo table is 10.Often adopt penalty function that the administrative model of constraint is converted to the unconfinement administrative model when using intelligent algorithm to be optimized design, promptly to not satisfying separating of constraint condition, the form of employing penalty function is added in the corresponding objective function, optimizes to the results are shown in Figure 5.
2.3 optimize result's comparative analysis
The calculating of objective function needs to call analogy model in a large number repeatedly, has occupied more than 90% of total computation optimization time its computing time, and for the simulative optimization problem of actual large scale, the ratio that the time occupies is infinitely near 100%.So the calculation times of objective function, promptly the call number of analogy model is a reflection algorithm computation efficient important indicator.Fig. 6 has compared NPTS and has used and do not using each generation under the ideal adaptation degree functional value archival strategy to call the number of times of analogy model.Under the situation of not using ideal adaptation degree functional value archival strategy, each generation the number of times that calls analogy model equal number of times product (n * i that decision variable number and neighborhood move
Max).As can be seen from Figure 6 along with the increase of iterations, the individuality of a large amount of double countings occurred, table 1 has been added up the number of times of the actual analogy model that calls of convergence algebraic sum, calculates the double counting rate of separating.The number of individuals of double counting reaches 831, and repetition rate reaches 38.0%, and this has proved absolutely that also ideal adaptation degree functional value archival strategy has improved algorithm computation efficient widely.
Table 1. uses and does not use the Pareto optimum solution statistical form of ideal adaptation degree functional value archival strategy
The number of times that neighborhood is separated in per generation | Optimize iterations when reaching convergence | Objective function calculation times when not using ideal adaptation degree functional value archival strategy | Objective function calculation times when using ideal adaptation degree functional value archival strategy | Objective function does not calculate repetition rate (%) when not using ideal adaptation degree functional value file plan |
??8 | ??273 | ??2184 | ??1353 | ??38.0 |
With the present invention with based on the multi-target evolution algorithm of genetic algorithm, comprise NPGA and VEGA, it is as follows that NPGA and VEGA are provided with identical parameter: population size, 100; Crossing-over rate P
c, 0.95; Aberration rate P
m, 0.05, all the other settings are identical with NPTS.The Pareto optimum solution that pairs of three kinds multiple goal algorithm optimizations of table 2 obtain is added up, NPTS optimizes that to obtain number and NPGA that Pareto separates suitable with the VEGA number, but the calculation times of objective function then is 1/9 and 1/10 of NPGA and VEGA, the result shows NPTS under the situation that guarantees global search, has improved the counting yield of optimizing greatly.
Table 2.NPTS and NPGA and the VEGA Pareto optimum solution contrast table under abundant convergence situation.
Claims (3)
1. a underground water is repaired Multipurpose Optimal Method, and it may further comprise the steps:
Step 1 is set up analogy model, is used to portray reparation place underground water head and the distribution of solute concentration on time and space;
Step 2 is determined management objectives, sets up Optimization Model;
Step 3 coupled simulation model and Optimization Model;
Step 4 is selected for use optimization method to find the solution the balance of multiobjective management problem and is separated;
The result is optimized in step 5 output, analyzes the of overall importance and counting yield that balance is separated, and adjusts parameters optimization;
It is characterized in that step 4 comprises following detailed process:
1) producing initial solution and corresponding neighborhood separates
At first, from initial solution, adopt random fashion to produce an initial solution to the search in feasible solution territory, it then is to be basic point with the initial solution that neighborhood is separated, and carries out neighborhood according to the neighborhood move function and moves the back and produce, and the neighborhood move function is defined as follows:
s
i=(2×random()-1)×stepsize,i=1,2,…,i
max????????????????(1)
S wherein
i(i=1 ..., i
Max) represent that neighborhood moves rule function, i
MaxBe the number of times that neighborhood moves, random () is the random function that produces real number between the 0-1, and stepsize is the neighborhood moving step length.Initial solution and neighborhood are separated to produce and are had randomness in this research;
2) calculate target function value, Pareto rank and the individual crowding that neighborhood is separated
After the generation neighborhood is separated, calculate the target function value that all neighborhoods are separated,, calculate Pareto rank and individual crowding value that each neighborhood is separated according to target function value; The Pareto rank number of separating equal neighborhood separate in all Pareto control the skill that this separates, the Pareto rank is 0 separates and be nondominated solution;
3) select kind of a subsolution
Kind of subsolution is the initial solution of searching in per generation, and the selection of planting subsolution is derived from the previous generation neighborhood and separates, and promptly finds to separate separating as candidate seed that Pareto is dominant with respect to other neighborhood and separate from neighborhood is separated, and must satisfy the uncontrolled property of Pareto;
4) upgrade Pareto table, candidate list and taboo table
Pareto table and candidate list are respectively applied for the kind subsolution of depositing each generation and do not have selected other candidate solution for kind of subsolution; To no longer have separating from Pareto table and candidate list of Pareto optimality in per generation search and remove, the individuality that is chosen as kind of subsolution in the candidate list will be removed from table; Upgrade the taboo table with kind of subsolution as the taboo object, selected kind subsolution is avoided in the taboo algebraically of regulation, and seed is lifted a ban the new taboo kind subsolution in back replacement is deposited in the taboo table;
5) judge whether to satisfy stopping criterion
In the optimizing process, if reached predefined maximum iteration time, seed selection subsolution collection also is empty for empty candidate list simultaneously when perhaps searching certain generation, then can't find the next subsolution of planting, algorithm can't enter the search in next stage, and algorithm stops, output Pareto optimal solution set.
2. underground water according to claim 1 is repaired Multipurpose Optimal Method, it is characterized in that process 2) in adopt ideal adaptation degree functional value archival strategy note first search to separate the information of separating in the direct call function value storehouse of separating then of repeat search.
3. underground water according to claim 1 and 2 is repaired Multipurpose Optimal Method, it is characterized in that process 3) in plant subsolution selection following several situation is arranged: if a) only exist a candidate seed to separate, then this is separated as follow-on kind of subsolution; B) if existence is separated more than one candidate seed, contrast the individual crowding that candidate seed is separated, the individuality of choosing crowding value minimum is as follow-on kind of subsolution; C) if there is no candidate seed is separated, and then selects to enter separating as new kind subsolution in the table the earliest from candidate list, and this is separated from candidate list delete.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010100183610A CN101799846B (en) | 2010-01-14 | 2010-01-14 | Multi-objective groundwater remediation optimization method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010100183610A CN101799846B (en) | 2010-01-14 | 2010-01-14 | Multi-objective groundwater remediation optimization method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101799846A true CN101799846A (en) | 2010-08-11 |
CN101799846B CN101799846B (en) | 2011-08-31 |
Family
ID=42595519
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010100183610A Expired - Fee Related CN101799846B (en) | 2010-01-14 | 2010-01-14 | Multi-objective groundwater remediation optimization method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101799846B (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102542076A (en) * | 2010-12-09 | 2012-07-04 | 东北电力科学研究院有限公司 | Correction method of simulation track of power grid frequency dynamic process |
CN104912025A (en) * | 2015-04-23 | 2015-09-16 | 河海大学 | Determination method for low-water ecological restoration project layout of straight stream channel |
CN105718632A (en) * | 2016-01-14 | 2016-06-29 | 河海大学 | Multi-target optimal management method for remediation of underground water with uncertainty |
CN105894132A (en) * | 2016-04-27 | 2016-08-24 | 河海大学 | Multi-target optimization management method for underground water in land subsidence region |
CN107423564A (en) * | 2017-07-26 | 2017-12-01 | 中国环境科学研究院 | The method of decision analysis of river basin ecological correcting strategy |
CN109034479A (en) * | 2018-07-28 | 2018-12-18 | 河南工业大学 | A kind of Multiobjective Scheduling method and device based on differential evolution algorithm |
CN109117517A (en) * | 2018-07-19 | 2019-01-01 | 江苏地基工程有限公司 | A kind of base pit dewatering optimum design method based on genetic algorithm and groundwater model |
CN109635510A (en) * | 2019-01-15 | 2019-04-16 | 西南交通大学 | The parallel repair line setting method of railway freight-car |
CN112836459A (en) * | 2021-02-02 | 2021-05-25 | 山东大学 | Circuit analysis model assisted analog circuit automatic optimization method, device and storage medium |
-
2010
- 2010-01-14 CN CN2010100183610A patent/CN101799846B/en not_active Expired - Fee Related
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102542076B (en) * | 2010-12-09 | 2013-10-02 | 东北电力科学研究院有限公司 | Correction method of simulation track of power grid frequency dynamic process |
CN102542076A (en) * | 2010-12-09 | 2012-07-04 | 东北电力科学研究院有限公司 | Correction method of simulation track of power grid frequency dynamic process |
CN104912025A (en) * | 2015-04-23 | 2015-09-16 | 河海大学 | Determination method for low-water ecological restoration project layout of straight stream channel |
CN104912025B (en) * | 2015-04-23 | 2016-09-28 | 河海大学 | A kind of determination method of straight river course low Water Ecological Recovery engineering layout |
CN105718632A (en) * | 2016-01-14 | 2016-06-29 | 河海大学 | Multi-target optimal management method for remediation of underground water with uncertainty |
CN105718632B (en) * | 2016-01-14 | 2018-10-19 | 河海大学 | A kind of uncertainty underground water repairs multiple-objection optimization management method |
CN105894132B (en) * | 2016-04-27 | 2019-05-21 | 河海大学 | A kind of surface subsidence area underground water multiple-objection optimization management method |
CN105894132A (en) * | 2016-04-27 | 2016-08-24 | 河海大学 | Multi-target optimization management method for underground water in land subsidence region |
CN107423564A (en) * | 2017-07-26 | 2017-12-01 | 中国环境科学研究院 | The method of decision analysis of river basin ecological correcting strategy |
CN109117517A (en) * | 2018-07-19 | 2019-01-01 | 江苏地基工程有限公司 | A kind of base pit dewatering optimum design method based on genetic algorithm and groundwater model |
CN109034479A (en) * | 2018-07-28 | 2018-12-18 | 河南工业大学 | A kind of Multiobjective Scheduling method and device based on differential evolution algorithm |
CN109635510A (en) * | 2019-01-15 | 2019-04-16 | 西南交通大学 | The parallel repair line setting method of railway freight-car |
CN112836459A (en) * | 2021-02-02 | 2021-05-25 | 山东大学 | Circuit analysis model assisted analog circuit automatic optimization method, device and storage medium |
CN112836459B (en) * | 2021-02-02 | 2022-09-20 | 山东大学 | Circuit analysis model assisted analog circuit automatic optimization method, device and storage medium |
Also Published As
Publication number | Publication date |
---|---|
CN101799846B (en) | 2011-08-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101799846B (en) | Multi-objective groundwater remediation optimization method | |
CN109670650B (en) | Multi-objective optimization algorithm-based solving method for cascade reservoir group scheduling model | |
CN105354646B (en) | Power load forecasting method for hybrid particle swarm optimization and extreme learning machine | |
CN105162151B (en) | A kind of grid-connected real-time control method of intelligent energy-storage system based on artificial fish-swarm algorithm | |
CN101996242B (en) | Three-dimensional R-tree index expansion structure-based three-dimensional city model adaptive method | |
CN105809279A (en) | Multi-objective quantum Shuffled Frog Leaping Algorithm (SFLA) based water resource optimization and diversion method | |
CN103902783B (en) | A kind of drainage pipeline networks optimization method dividing algorithm based on the reverse poor learning of broad sense | |
CN106682355A (en) | Hydrological model parameter calibration method based on PSO (particle swarm optimization)-GA (genetic algorithm) mixed algorithm | |
CN105718632B (en) | A kind of uncertainty underground water repairs multiple-objection optimization management method | |
Tung et al. | Application of the genetic algorithm for optimizing operation rules of the liyutan reservoir in taiwan 1 | |
CN107704965B (en) | A kind of reservoir usable regulation figure optimization method for taking into account ecological protection demand | |
CN104331749A (en) | AGV optimization scheduling method based on simulated annealing particle swarm | |
CN103150614A (en) | Automatic configuration method for land utilization space | |
CN106651011B (en) | Canal system optimized water distribution method based on particle swarm optimization | |
CN104732091A (en) | Cellular automaton river bed evolution prediction method based on natural selection ant colony algorithm | |
CN107939372B (en) | Optimal well position deployment method and device for small fault block oil reservoir | |
CN103473465A (en) | Method for optimizing spatial layout of land resources based on multi-target artificial immunization system | |
CN104993966A (en) | Power integrated service network flow prediction method | |
CN105894132B (en) | A kind of surface subsidence area underground water multiple-objection optimization management method | |
CN108915927B (en) | Method, device and equipment for determining parameters of wave power generation device | |
CN108133286B (en) | Underground water multi-target calculation method based on ground settlement substitution model | |
CN106780747A (en) | A kind of method that Fast Segmentation CFD calculates grid | |
CN113592318B (en) | Urban green space optimization layout model building method and system | |
CN109787247A (en) | A kind of reactive compensation planing method based on multi-parametric programming | |
CN109117517A (en) | A kind of base pit dewatering optimum design method based on genetic algorithm and groundwater model |
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 | ||
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: 20110831 Termination date: 20130114 |