Well Optimisation With Goal-Based Sensitivity Maps Using Time Windows And Ensemble Perturbations
C.E. Heaney, P. Salinas, C.C. Pain, F. Fang and I.M. Navon
Event name: ECMOR XVI - 16th European Conference on the Mathematics of Oil Recovery
Session: ECMOR XVI Poster Session 2
Publication date: 03 September 2018
Info: Extended abstract, PDF ( 781.42Kb )
Price: € 20
Knowledge of the sensitivity of a solution to small changes in the model parameters is exploited in many areas in computational physics and used to perform mesh adaptivity, or to correct errors based on discretisation and sub-grid-scale modelling errors, to perform the assimilation of data based on adjusting the most sensitive parameters to the model-observation misfit, and similarly to form optimised sub-grid-scale models. We present a goal-based approach for forming sensitivity (or importance) maps using ensembles. These maps are defined as regions in space and time of high relevance for a given goal, for example, the solution at an observation point within the domain. The presented approach relies solely on ensembles obtained from the forward model and thus can be used with complex models for which calculating an adjoint is not a practical option. This provides a simple approach for optimisation of sensor placement, goal based mesh adaptivity, assessment of goals and data assimilation. We investigate methods which reduce the number of ensembles used to construct the maps yet which retain reasonable fidelity of the maps. The fidelity comes from an integrated method including a goal-based approach, in which the most up-to-date importance maps are fed back into the perturbations to focus the algorithm on the key variables and domain areas. Also within the method smoothing is applied to the perturbations to obtain a multi-scale, global picture of the sensitivities; the perturbations are orthogonalised in order to generate a well-posed system which can be inverted; and time windows are applied (for time dependent problems) where we work backwards in time to obtain greater accuracy of the sensitivity maps. The approach is demonstrated on a multi-phase flow problem.