The facies distribution of a reservoir is one of the biggest concerns for geologists,geophysicists,reservoir modelers,and reservoir engineers due to its high importance in the setting of any reliable decisionmaking/op...The facies distribution of a reservoir is one of the biggest concerns for geologists,geophysicists,reservoir modelers,and reservoir engineers due to its high importance in the setting of any reliable decisionmaking/optimization of field development planning.The approach for parameterizing the facies distribution as a random variable comes naturally through using the probability fields.Since the prior probability fields of facies come either from a seismic inversion or from other sources of geologic information,they are not conditioned to the data observed from the cores extracted from the wells.This paper presents a regularized element-free Galerkin(R-EFG)method for conditioning facies probability fields to facies observation.The conditioned probability fields respect all the conditions of the probability theory(i.e.all the values are between 0 and 1,and the sum of all fields is a uniform field of 1).This property achieves by an optimization procedure under equality and inequality constraints with the gradient projection method.The conditioned probability fields are further used as the input in the adaptive pluri-Gaussian simulation(APS)methodology and coupled with the ensemble smoother with multiple data assimilation(ES-MDA)for estimation and uncertainty quantification of the facies distribution.The history-matching of the facies models shows a good estimation and uncertainty quantification of facies distribution,a good data match and prediction capabilities.展开更多
The Ensemble Kalman Filter(EnKF),as the most popular sequential data assimilation algorithm for history matching,has the intrinsic problem of high computational cost and the potential inconsistency of state variables ...The Ensemble Kalman Filter(EnKF),as the most popular sequential data assimilation algorithm for history matching,has the intrinsic problem of high computational cost and the potential inconsistency of state variables updated at each loop of data assimilation and its corresponding reservoir simulated result.This problem forbids the reservoir engineers to make the best use of the 4D seismic data,which provides valuable information about the fluid change inside the reservoir.Moreover,only matching the production data in the past is not enough to accurately forecast the future,and the development plan based on the false forecast is very likely to be suboptimal.To solve this problem,we developed a workflow for geophysical and production data history matching by modifying ensemble smoother with multiple data assimilation(ESMDA).In this work,we derived the mathematical expressions of ESMDA and discussed its scope of applications.The geophysical data we used is P-wave impedance,which is typically included in a basic seismic interpretation,and it directly reflects the saturation change in the reservoir.Full resolution of the seismic data is not necessary,we subsampled the P-wave impedance data to further reduce the computational cost.With our case studies on a benchmark synthetic reservoir model,we also showed the supremacy of matching both geophysical and production data,than the traditional reservoir history matching merely on the production data:the overall percentage error of the observed data is halved,and the variances of the updated forecasts are reduced by two orders of the magnitude.展开更多
文摘The facies distribution of a reservoir is one of the biggest concerns for geologists,geophysicists,reservoir modelers,and reservoir engineers due to its high importance in the setting of any reliable decisionmaking/optimization of field development planning.The approach for parameterizing the facies distribution as a random variable comes naturally through using the probability fields.Since the prior probability fields of facies come either from a seismic inversion or from other sources of geologic information,they are not conditioned to the data observed from the cores extracted from the wells.This paper presents a regularized element-free Galerkin(R-EFG)method for conditioning facies probability fields to facies observation.The conditioned probability fields respect all the conditions of the probability theory(i.e.all the values are between 0 and 1,and the sum of all fields is a uniform field of 1).This property achieves by an optimization procedure under equality and inequality constraints with the gradient projection method.The conditioned probability fields are further used as the input in the adaptive pluri-Gaussian simulation(APS)methodology and coupled with the ensemble smoother with multiple data assimilation(ES-MDA)for estimation and uncertainty quantification of the facies distribution.The history-matching of the facies models shows a good estimation and uncertainty quantification of facies distribution,a good data match and prediction capabilities.
基金supported by Chinese National Science and Technology Major Project(2016ZX05015-005).
文摘The Ensemble Kalman Filter(EnKF),as the most popular sequential data assimilation algorithm for history matching,has the intrinsic problem of high computational cost and the potential inconsistency of state variables updated at each loop of data assimilation and its corresponding reservoir simulated result.This problem forbids the reservoir engineers to make the best use of the 4D seismic data,which provides valuable information about the fluid change inside the reservoir.Moreover,only matching the production data in the past is not enough to accurately forecast the future,and the development plan based on the false forecast is very likely to be suboptimal.To solve this problem,we developed a workflow for geophysical and production data history matching by modifying ensemble smoother with multiple data assimilation(ESMDA).In this work,we derived the mathematical expressions of ESMDA and discussed its scope of applications.The geophysical data we used is P-wave impedance,which is typically included in a basic seismic interpretation,and it directly reflects the saturation change in the reservoir.Full resolution of the seismic data is not necessary,we subsampled the P-wave impedance data to further reduce the computational cost.With our case studies on a benchmark synthetic reservoir model,we also showed the supremacy of matching both geophysical and production data,than the traditional reservoir history matching merely on the production data:the overall percentage error of the observed data is halved,and the variances of the updated forecasts are reduced by two orders of the magnitude.
