Christini Fandridi, "Optimization of petroleum production under industrial constraints using alternative objective functions and adjoint gradient-based techniques.", Master Thesis, School of Mineral Resources Engineering, Technical University of Crete, Chania, Greece, 2017
https://doi.org/10.26233/heallink.tuc.67321
The optimization of oil production is a tedious and computationally intensive pro-cess that requires the solution of time dependent nonlinear set of partial differ-ential equations describing the flow of hydrocarbons in anisotropic porous me-dia. Optimization of production is usually performed using either gradient freetechniques like genetic algorithms, particle swarm algorithms, or gradient-basedtechniques where the gradients are computed through the solution of the adjointproblem. A gradient-based optimization method, in which the gradient is com-puted using an adjoint formulation, is often the method of choice since in con-trast to numerical perturbation techniques that require as many objective functionevaluations as the number of control parameters, the gradient using adjoint-basedtechniques is obtained only at a small fraction of the time spent for the evaluationof the objective function. It is well known that for non-convex optimisation prob-lems, gradient-based techniques are likely to get trapped in poor local optima. Acommon practise is to lunch several independent optimisation runs from differentinitial guesses or to combine ideas from gradient-free algorithms with gradient-based to benefit from the merits of both. An adequate sampling of the search space would require an intractable number of simulations and it is thus impossible.The aim of this work is to exploit an observation in homogeneous reservoirs,where the global optimum, when optimising cumulative oil recovery, is usuallyachieved from practically any initial guess. This observation suggest to optimizecumulative oil by adopting a “geology continuation” method. In this novel ap-proach the porosity and permeability fields, gradually switch from some averagehomogeneous values chosen heuristically for the particular benchmark, to the in-homogeneous geological properties characterizing the reservoir. The optimal con-trols from each step become the initial controls to the next step.In addition instead of maximizing the cumulative oil we suggest to minimize mod-ified versions of the residual oil function which are likely to be more convex andthus less likely to lead in poor local optima.