skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Exploring the effects of data quality, data worth, and redundancy of CO2 gas pressure and saturation data on reservoir characterization through PEST Inversion

Abstract

This study examined the impacts of reservoir properties on CO2 migration after subsurface injection and evaluated the possibility of characterizing reservoir properties using CO2 monitoring data such as saturation distribution. The injection reservoir was assumed to be located 1400-1500 m below the ground surface such that CO2 remained in the supercritical state. The reservoir was assumed to contain layers with alternating conductive and resistive properties, which is analogous to actual geological formations such as the Mount Simon Sandstone unit. The CO2 injection simulation used a cylindrical grid setting in which the injection well was situated at the center of the domain, which extended up to 8000 m from the injection well. The CO2 migration was simulated using the PNNL-developed simulator STOMP-CO2e (the water-salt-CO2 module). We adopted a nonlinear parameter estimation and optimization modeling software package, PEST, for automated reservoir parameter estimation. We explored the effects of data quality, data worth, and data redundancy on the detectability of reservoir parameters using CO2 saturation monitoring data, by comparing PEST inversion results using data with different levels of noises, various numbers of monitoring wells and locations, and different data collection spacing and temporal sampling intervals. This study yielded insight into the use ofmore » CO2 saturation monitoring data for reservoir characterization and how to design the monitoring system to optimize data worth and reduce data redundancy.« less

Authors:
; ; ; ; ;
Publication Date:
Research Org.:
Pacific Northwest National Lab. (PNNL), Richland, WA (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1124030
Report Number(s):
PNNL-SA-88296
DOE Contract Number:
AC05-76RL01830
Resource Type:
Journal Article
Resource Relation:
Journal Name: Environmental Earth Sciences, 71(7):3025-3037
Country of Publication:
United States
Language:
English
Subject:
Carbon sequestration; data worth; inversion; monitoring network

Citation Formats

Fang, Zhufeng, Hou, Zhangshuan, Lin, Guang, Engel, David W., Fang, Yilin, and Eslinger, Paul W. Exploring the effects of data quality, data worth, and redundancy of CO2 gas pressure and saturation data on reservoir characterization through PEST Inversion. United States: N. p., 2014. Web. doi:10.1007/s12665-013-2680-9.
Fang, Zhufeng, Hou, Zhangshuan, Lin, Guang, Engel, David W., Fang, Yilin, & Eslinger, Paul W. Exploring the effects of data quality, data worth, and redundancy of CO2 gas pressure and saturation data on reservoir characterization through PEST Inversion. United States. doi:10.1007/s12665-013-2680-9.
Fang, Zhufeng, Hou, Zhangshuan, Lin, Guang, Engel, David W., Fang, Yilin, and Eslinger, Paul W. Tue . "Exploring the effects of data quality, data worth, and redundancy of CO2 gas pressure and saturation data on reservoir characterization through PEST Inversion". United States. doi:10.1007/s12665-013-2680-9.
@article{osti_1124030,
title = {Exploring the effects of data quality, data worth, and redundancy of CO2 gas pressure and saturation data on reservoir characterization through PEST Inversion},
author = {Fang, Zhufeng and Hou, Zhangshuan and Lin, Guang and Engel, David W. and Fang, Yilin and Eslinger, Paul W.},
abstractNote = {This study examined the impacts of reservoir properties on CO2 migration after subsurface injection and evaluated the possibility of characterizing reservoir properties using CO2 monitoring data such as saturation distribution. The injection reservoir was assumed to be located 1400-1500 m below the ground surface such that CO2 remained in the supercritical state. The reservoir was assumed to contain layers with alternating conductive and resistive properties, which is analogous to actual geological formations such as the Mount Simon Sandstone unit. The CO2 injection simulation used a cylindrical grid setting in which the injection well was situated at the center of the domain, which extended up to 8000 m from the injection well. The CO2 migration was simulated using the PNNL-developed simulator STOMP-CO2e (the water-salt-CO2 module). We adopted a nonlinear parameter estimation and optimization modeling software package, PEST, for automated reservoir parameter estimation. We explored the effects of data quality, data worth, and data redundancy on the detectability of reservoir parameters using CO2 saturation monitoring data, by comparing PEST inversion results using data with different levels of noises, various numbers of monitoring wells and locations, and different data collection spacing and temporal sampling intervals. This study yielded insight into the use of CO2 saturation monitoring data for reservoir characterization and how to design the monitoring system to optimize data worth and reduce data redundancy.},
doi = {10.1007/s12665-013-2680-9},
journal = {Environmental Earth Sciences, 71(7):3025-3037},
number = ,
volume = ,
place = {United States},
year = {Tue Apr 01 00:00:00 EDT 2014},
month = {Tue Apr 01 00:00:00 EDT 2014}
}
  • In this study we developed an efficient Bayesian inversion framework for interpreting marine seismic amplitude versus angle (AVA) and controlled source electromagnetic (CSEM) data for marine reservoir characterization. The framework uses a multi-chain Markov-chain Monte Carlo (MCMC) sampler, which is a hybrid of DiffeRential Evolution Adaptive Metropolis (DREAM) and Adaptive Metropolis (AM) samplers. The inversion framework is tested by estimating reservoir-fluid saturations and porosity based on marine seismic and CSEM data. The multi-chain MCMC is scalable in terms of the number of chains, and is useful for computationally demanding Bayesian model calibration in scientific and engineering problems. As a demonstration,more » the approach is used to efficiently and accurately estimate the porosity and saturations in a representative layered synthetic reservoir. The results indicate that the seismic AVA and CSEM joint inversion provides better estimation of reservoir saturations than the seismic AVA-only inversion, especially for the parameters in deep layers. The performance of the inversion approach for various levels of noise in observational data was evaluated – reasonable estimates can be obtained with noise levels up to 25%. Sampling efficiency due to the use of multiple chains was also checked and was found to have almost linear scalability.« less
  • In this paper we developed an efficient Bayesian inversion framework for interpreting marine seismic Amplitude Versus Angle and Controlled-Source Electromagnetic data for marine reservoir characterization. The framework uses a multi-chain Markov-chain Monte Carlo sampler, which is a hybrid of DiffeRential Evolution Adaptive Metropolis and Adaptive Metropolis samplers. The inversion framework is tested by estimating reservoir-fluid saturations and porosity based on marine seismic and Controlled-Source Electromagnetic data. The multi-chain Markov-chain Monte Carlo is scalable in terms of the number of chains, and is useful for computationally demanding Bayesian model calibration in scientific and engineering problems. As a demonstration, the approach ismore » used to efficiently and accurately estimate the porosity and saturations in a representative layered synthetic reservoir. The results indicate that the seismic Amplitude Versus Angle and Controlled-Source Electromagnetic joint inversion provides better estimation of reservoir saturations than the seismic Amplitude Versus Angle only inversion, especially for the parameters in deep layers. The performance of the inversion approach for various levels of noise in observational data was evaluated — reasonable estimates can be obtained with noise levels up to 25%. Sampling efficiency due to the use of multiple chains was also checked and was found to have almost linear scalability.« less
  • Ground deformation, commonly seen in storage projects, carries useful information about processes occurring in the injection formation. The Krechba gas field at In Salah (Algeria) is one of the best-known sites for studying ground surface deformation during geological carbon storage. At this first industrial-scale on-shore CO 2 demonstration project, satellite-based ground-deformation monitoring data of high quality are available and used to study the large-scale hydrological and geomechanical response of the system to injection. In this work, we carry out coupled fluid flow and geomechanical simulations to understand the uplift at three different CO 2 injection wells (KB-501, KB-502, KB-503). Previousmore » numerical studies focused on the KB-502 injection well, where a double-lobe uplift pattern has been observed in the ground-deformation data. The observed uplift patterns at KB-501 and KB-503 have single-lobe patterns, but they can also indicate a deep fracture zone mechanical response to the injection.The current study improves the previous modeling approach by introducing an injection reservoir and a fracture zone, both responding to a Mohr-Coulomb failure criterion. In addition, we model a stress-dependent permeability and bulk modulus, according to a dual continuum model. Mechanical and hydraulic properties are determined through inverse modeling by matching the simulated spatial and temporal evolution of uplift to InSAR observations as well as by matching simulated and measured pressures. The numerical simulations are in agreement with both spatial and temporal observations. The estimated values for the parameterized mechanical and hydraulic properties are in good agreement with previous numerical results. In addition, the formal joint inversion of hydrogeological and geomechanical data provides measures of the estimation uncertainty.« less
  • We report seismic full-waveform inversion (FWI) methods hold strong potential to recover multiple subsurface elastic properties for hydrocarbon reservoir characterization. Simultaneously updating multiple physical parameters introduces the problem of interparameter tradeoff, arising from the covariance between different physical parameters, which increases nonlinearity and uncertainty of multiparameter FWI. The coupling effects of different physical parameters are significantly influenced by model parameterization and acquisition arrangement. An appropriate choice of model parameterization is critical to successful field data applications of multiparameter FWI. The objective of this paper is to examine the performance of various model parameterizations in isotropic-elastic FWI with walk-away vertical seismicmore » profile (W-VSP) dataset for unconventional heavy oil reservoir characterization. Six model parameterizations are considered: velocity-density (α, β and ρ'), modulus-density (κ, μ and ρ), Lamé-density (λ, μ' and ρ'''), impedance-density (IP, IS and ρ''), velocity-impedance-I (α', β' and I' P), and velocity-impedance-II (α'', β'' and I'S). We begin analyzing the interparameter tradeoff by making use of scattering radiation patterns, which is a common strategy for qualitative parameter resolution analysis. In this paper, we discuss the advantages and limitations of the scattering radiation patterns and recommend that interparameter tradeoffs be evaluated using interparameter contamination kernels, which provide quantitative, second-order measurements of the interparameter contaminations and can be constructed efficiently with an adjoint-state approach. Synthetic W-VSP isotropic-elastic FWI experiments in the time domain verify our conclusions about interparameter tradeoffs for various model parameterizations. Density profiles are most strongly influenced by the interparameter contaminations; depending on model parameterization, the inverted density profile can be over-estimated, under-estimated or spatially distorted. Among the six cases, only the velocity-density parameterization provides stable and informative density features not included in the starting model. Field data applications of multicomponent W-VSP isotropic-elastic FWI in the time domain were also carried out. Finally, the heavy oil reservoir target zone, characterized by low α-to-β ratios and low Poisson’s ratios, can be identified clearly with the inverted isotropic-elastic parameters.« less
  • The porous plate method was used to determine capillary pressure curves and saturation exponents at ambient conditions and at full reservoir conditions (reservoir pressure, temperature, and fluid) on sandstone core plugs from a North Sea reservoir. This paper reports that the main objective of this study is to determine whether any difference exists between the capillary pressure curves and the saturation exponents obtained at ambient conditions and those obtained at reservoir conditions for a North Sea reservoir. Results show that the capillary pressure and electrical resistivity measurements performed at ambient conditions with nonreservoir fluid displacing water can give results differentmore » from measurements performed at conditions prevailing in the reservoir.« less