Data Assimilation of Coupled Fluid Flow and Geomechanics Using the Ensemble Kalman Filter

Summary In reservoir history matching or data assimilation, dynamic data, such as production rates and pressures, are used to constrain reservoir models and to update model parameters. As such, even if under certain conceptualization the model parameters do not vary with time, the estimate of such parameters may change with the available observations and, thus, with time. In reality, the production process may lead to changes in both the flow and geomechanics fields, which are dynamically coupled. For example, the variations in the stress/strain field lead to changes in porosity and permeability of the reservoir and, hence, in the flow field. In weak formations, such as the Lost Hills oil field, fluid extraction may cause a large compaction to the reservoir rock and a significant subsidence at the land surface, resulting in huge economic losses and detrimental environmental consequences. The strong nonlinear coupling between reservoir flow and geomechanics poses a challenge to constructing a reliable model for predicting oil recovery in such reservoirs. On the other hand, the subsidence and other geomechanics observations can provide additional insight into the nature of the reservoir rock and help constrain the reservoir model if used wisely. In this study, the ensemble-Kalman-filter (EnKF) approach is used to estimate reservoir flow and material properties by jointly assimilating dynamic flow and geomechanics observations. The resulting model can be used for managing and optimizing production operations and for mitigating the land subsidence. The use of surface displacement observations improves the match to both production and displacement data. Localization is used to facilitate the assimilation of a large amount of data and to mitigate the effect of spurious correlations resulting from small ensembles. Because the stress, strain, and displacement fields are updated together with the material properties in the EnKF, the issue of consistency at the analysis step of the EnKF is investigated. A 3D problem with reservoir fluid-flow and mechanical parameters close to those of the Lost Hills oil field is used to test the applicability.