Production Uncertainty Quantification for a Naturally Fractured Reservoir using a New Data-Space Inversion Procedure

A new method for production forecasting and uncertainty quantification, applicable for realistic naturally fractured reservoirs (NFRs) represented as general discrete-fracture-matrix (DFM) models, is developed and applied. The forecasting procedure extends a recently developed data-space inversion (DSI) technique that generates production predictions using only prior-model simulation results and observed data. The method does not provide posterior (history-matched) geological models. Rather, the DSI method treats production data as random variables. The prior distribution is estimated from the flow simulations performed on prior geological models, and the posterior data-variable distribution is sampled using a data-space randomized maximum likelihood method. The DSI treatment requires the parameterization of data variables to render them approximately multivariate Gaussian. The complex production data considered here (resulting from frequent well shut-ins) is treated using a new reparameterization that involves principal component analysis combined with histogram transformation. The DSI procedure is applied to a real NFR that has undergone years of primary production and is now under consideration for waterflooding. To construct the DSI representation, 400 prior DFM models, which correspond to different geologic concepts and properties, are simulated. Two different reference `true’ models, along with different data-assimilation durations, are considered to evaluate the performance of the DSI procedure. In all cases, the DSI predictions are shown to be consistent with the forecasts from the `true’ model, and to provide reasonable quantification of the forecast uncertainty.