DERIVATIVE-FREE OPTIMIZATION METHODS FOR HANDLING FIXED COSTS IN OPTIMAL GROUNDWATER REMEDIATION DESIGN

We consider a hydraulic capture application for water resources management that includes a fixed installation cost in addition to operating costs. The result is a simulationbased, nonlinear, mixed-integer optimization problem. The motivation is that our preliminary studies have shown that convergence to an unsatisfactory, local minimum with many wells operating at low pumping rates is common when the fixed cost is ignored. Such optimization tasks are not unique to subsurface management, rather efficient simulation-based methods are needed in the whole field of computational engineering. All the approaches used below do not need the gradient of the objective function, only function values for minimization. In one approach, we bypass including the number of wells as a decision variable by defining an inactive-well threshold. In another approach, we use penalty coefficients proposed in the literature to transform the discontinuous problem into a continuous one. For the two above formulations, we use the implicit filtering algorithm. In the third approach, we introduce a mixed-integer problem formulation and use an iterative stochastic modeling technique to build surrogate functions that approximate the objective function. With this new procedure the use of a branch-and-bound technique becomes possible to solve the mixed-integer problem in contrast to methods working directly on the simulation results, which impedes relaxation of integer variables. We present promising numerical results on the benchmarking problem and point the way towards improvement and future work.

[1]  Arthur H. Delmege,et al.  Design-by-Computer , 1971 .

[2]  Philip E. Gill,et al.  Practical optimization , 1981 .

[3]  Arlen W. Harbaugh,et al.  A modular three-dimensional finite-difference ground-water flow model , 1984 .

[4]  P. Gilmore,et al.  IFFCO (Implicit filtering for constrained optimization), users' guide , 1993 .

[5]  C. T. Kelley,et al.  An Implicit Filtering Algorithm for Optimization of Functions with Many Local Minima , 1995, SIAM J. Optim..

[6]  D. McKinney,et al.  Approximate Mixed‐Integer Nonlinear Programming Methods for Optimal Aquifer Remediation Design , 1995 .

[7]  George F. Pinder,et al.  Optimal Ground‐Water Remediation Methods Applied to a Superfund Site: From Formulation to Implementation , 1995 .

[8]  Christodoulos A. Floudas,et al.  Deterministic global optimization - theory, methods and applications , 2010, Nonconvex optimization and its applications.

[9]  Cass T. Miller,et al.  Solution of a Groundwater Control Problem with Implicit Filtering , 2002 .

[10]  Michael James Sasena,et al.  Flexibility and efficiency enhancements for constrained global design optimization with kriging approximations. , 2002 .

[11]  Søren Nymand Lophaven,et al.  DACE - A Matlab Kriging Toolbox , 2002 .

[12]  I. Grossmann Review of Nonlinear Mixed-Integer and Disjunctive Programming Techniques , 2002 .

[13]  Michael A. Saunders,et al.  USER’S GUIDE FOR SNOPT 5.3: A FORTRAN PACKAGE FOR LARGE-SCALE NONLINEAR PROGRAMMING , 2002 .

[14]  Cass T. Miller,et al.  Solution of a Well-Field Design Problem with Implicit Filtering , 2004 .

[15]  Christopher E. Kees,et al.  A hydraulic capture application for optimal remediation design , 2004 .

[16]  Barbara S. Minsker,et al.  Which Groundwater Remediation Objective is Better: A Realistic One or a Simple One? , 2005 .

[17]  Herbert De Gersem,et al.  Mixed-integer simulation-based optimization for a superconductive magnet design , 2006 .

[18]  James R. Craig,et al.  Pump-and-treat optimization using analytic element method flow models , 2006 .