Rock-physics templates based on seismic Q

Geophysics ◽  
2019 ◽  
Vol 84 (1) ◽  
pp. MR13-MR23 ◽  
Author(s):  
Stefano Picotti ◽  
José M. Carcione ◽  
Jing Ba

We build rock-physics templates (RPTs) for reservoir rocks based on seismic quality factors. In these templates, the effects of partial saturation, porosity, and permeability on the seismic properties are described by generalizing the Johnson mesoscopic-loss model to a distribution of gas-patch sizes in brine- and oil-saturated rocks. This model addresses the wave-induced fluid flow attenuation mechanism, by which part of the energy of the fast P-wave is converted into the slow P (Biot) diffusive mode. We consider patch sizes, whose probability density function is defined by a normal (Gaussian) distribution. The complex bulk modulus of the composite medium is obtained with the Voigt-Reuss-Hill average, and we show that the results are close to those obtained with the Hashin-Shtrikman average. The templates represent the seismic dissipation factor (reciprocal of seismic quality factor) as a function of the P-wave velocity, acoustic impedance, and [Formula: see text] (P to S velocity ratio), for isolines of saturation, porosity, and permeability. They differentiate between oil and brine on the basis of the quality factor, with the gas-brine case showing more dissipation than the gas-oil case. We obtain sensitivity maps of the seismic properties to gas saturation and porosity for brine and oil. Unlike the gas-brine case, which shows higher sensitivity of attenuation to gas saturation, the gas-oil case shows higher sensitivity to porosity, and higher acoustic impedance and [Formula: see text] sensitivity values versus saturation. The RPTs can be used for a robust sensitivity analysis, which provides insights on seismic attributes for hydrocarbon detection and reservoir delineation. The templates are also relevant for studies related to [Formula: see text]-storage monitoring.

2021 ◽  
Vol 40 (10) ◽  
pp. 716-722
Author(s):  
Yangjun (Kevin) Liu ◽  
Michelle Ellis ◽  
Mohamed El-Toukhy ◽  
Jonathan Hernandez

We present a basin-wide rock-physics analysis of reservoir rocks and fluid properties in Campeche Basin. Reservoir data from discovery wells are analyzed in terms of their relationship between P-wave velocity, density, porosity, clay content, Poisson's ratio (PR), and P-impedance (IP). The fluid properties are computed by using in-situ pressure, temperature, American Petroleum Institute gravity, gas-oil ratio, and volume of gas, oil, and water. Oil- and gas-saturated reservoir sands show strong PR anomalies compared to modeled water sand at equivalent depth. This suggests that PR anomalies can be used as a direct hydrocarbon indicator in the Tertiary sands in Campeche Basin. However, false PR anomalies due to residual gas or oil exist and compose about 30% of the total anomalies. The impact of fluid properties on IP and PR is calibrated using more than 30 discovery wells. These calibrated relationships between fluid properties and PR can be used to guide or constrain amplitude variation with offset inversion for better pore fluid discrimination.


Energies ◽  
2021 ◽  
Vol 14 (21) ◽  
pp. 7225
Author(s):  
Chuantong Ruan ◽  
Jing Ba ◽  
José M. Carcione ◽  
Tiansheng Chen ◽  
Runfa He

Low porosity-permeability structures and microcracks, where gas is produced, are the main characteristics of tight sandstone gas reservoirs in the Sichuan Basin, China. In this work, an analysis of amplitude variation with offset (AVO) is performed. Based on the experimental and log data, sensitivity analysis is performed to sort out the rock physics attributes sensitive to microcrack and total porosities. The Biot–Rayleigh poroelasticity theory describes the complexity of the rock and yields the seismic properties, such as Poisson’s ratio and P-wave impedance, which are used to build rock-physics templates calibrated with ultrasonic data at varying effective pressures. The templates are then applied to seismic data of the Xujiahe formation to estimate the total and microcrack porosities, indicating that the results are consistent with actual gas production reports.


Geophysics ◽  
2016 ◽  
Vol 81 (1) ◽  
pp. WA233-WA245 ◽  
Author(s):  
Shan Dou ◽  
Seiji Nakagawa ◽  
Douglas Dreger ◽  
Jonathan Ajo-Franklin

Saline permafrost is sensitive to thermal disturbances and is prone to subsidence, which renders it a major source of geohazard in Arctic coastal environments. Seismic methods could be used to map and monitor saline permafrost at scales of geotechnical interests because of the ice-content dependencies of seismic properties. We have developed a comprehensive study of the ultrasonic P-wave properties (i.e., velocity and attenuation) of synthetic saline permafrost samples for a range of salinities and temperatures, and measurements conducted on a fine-grained permafrost core obtained from Barrow, Alaska. The resulting data consist of P-wave properties presented as functions of temperature and salinity. Notable observations include the following: P-wave velocities showed marked reductions in the presence of dissolved salts and complex variations resulting from the water-to-ice phase transitions; strong P-wave attenuation was present in the temperature intervals in which the samples were partially frozen. When presented as functions of ice saturation, the data sets lead us to two key findings: (1) neither a purely cementing nor a purely pore-filling model of the pore-scale distributions of ice could adequately fit the observed velocity data and (2) although the velocities increase monotonically with increasing ice saturations, P-wave attenuation reaches a maximum at intermediate ice saturations—contrary to the ordinary expectation of decreasing attenuation with increasing velocities. The observed ice-content dependencies of P-wave properties, along with the implications on the probable pore-scale distributions of ice, provide a valuable basis for rock-physics modeling, which in turn could facilitate seismic characterizations of saline permafrost.


2021 ◽  
Vol 9 ◽  
Author(s):  
Han Jin ◽  
Cai Liu ◽  
Zhiqi Guo ◽  
Yiming Zhang ◽  
Cong Niu ◽  
...  

Gas identification using seismic data is challenging for tight gas reservoirs with low porosity and permeability due to the complicated poroelastic behaviors of tight sandstone. In this study, the Chapman theory was used to simulate the dispersion and attenuation caused by the squirt flow of fluids in the complex pore spaces, which are assumed to consist of high aspect-ratio pores (stiff pores) and low aspect-ratio microcracks (soft pores). The rock physics modeling revealed that as the gas saturation varies, P-wave velocity dispersion and attenuation occurs at seismic frequencies, and it tends to move to high frequencies as the gas saturation increases. The velocity dispersion of the tight gas sandstone causes a frequency-dependent contrast in the P-wave impedance between the tight sandstone and the overlying mudstone, which consequently leads to frequency-dependent incidence reflection coefficients across the interface. In the synthetic seismic AVO modeling conducted by integrating the rock physics model and the propagator matrix method, the variations in the amplitudes and phases of the PP reflections can be observed for various gas saturations. The tests of the frequency-dependent AVO inversion of these synthetic data revealed that the magnitude of the inverted P-wave dispersion attribute can be used to indicate gas saturation in tight sandstone reservoirs. The applications of the frequency-dependent AVO inversion to the field pre-stacked seismic data revealed that the obtained P-wave dispersion attribute is positively correlated with the gas production from the pay zone at the well locations. Thus, the methods of the rock physics modeling and the frequency-dependent AVO inversion conducted in this study have good potential for the evaluation of the gas saturation in tight gas sandstone reservoirs.


Geophysics ◽  
2017 ◽  
Vol 82 (1) ◽  
pp. O11-O22 ◽  
Author(s):  
James Beckwith ◽  
Roger Clark ◽  
Linda Hodgson

The intrinsic seismic quality factor [Formula: see text] is known from poroelastic rock-physics theory to be frequency dependent, even within typical bandwidths of individual surface- and borehole-based surveys in which measurement methods usually deliver frequency-independent [Formula: see text]. Thus, measuring frequency-dependent [Formula: see text] instead offers better characterization of seismic properties and moreover a potential step toward estimating permeability directly from seismic data. Therefore, we have introduced a method to measure frequency-dependent [Formula: see text] from pairs of reflections in prestack [Formula: see text]-[Formula: see text] domain surface seismic data — a data type that, unlike a vertical seismic profile, offers useful areal coverage. Although, in principle, any analytic form with a manageable number of parameters could be prescribed, the frequency dependence of [Formula: see text] is modeled as a power law, [Formula: see text]. Inversion is done with a simple grid search over coefficient ([Formula: see text]) and exponent [Formula: see text], seeking a minimum [Formula: see text]-norm. We have found, using a numerical experiment and a synthetic data set, that it is robust and also accurate down to a signal-to-noise ratio of approximately 0.65. Then, [Formula: see text] is estimated for some 955 [Formula: see text] superbins of a 3D prestack ocean bottom cable data set over the Kinnoull field, central North Sea. All combinations of eight prominent reflections between Top Miocene and Base Cretaceous were treated together to give some 21,000 frequency-dependent and (for comparison) constant-[Formula: see text] results. The median coefficient ([Formula: see text]) and exponent [Formula: see text] were 0.0074 and 0.06, respectively, with sharply peaked distributions (excess kurtosis [Formula: see text]). Outlier, strongly frequency-dependent results, given by [Formula: see text], coincide with low-frequency “shadows” under amplitude anomalies, adversely affecting the spectra of reflections. The inferred frequency-dependent [Formula: see text] at 32.5 Hz, the center of the available bandwidth, is not statistically different from the frequency-independent [Formula: see text], 181 with a standard error from the distribution of one, derived from the same data. Hence for this data set, a constant-[Formula: see text] assumption would in fact be adequate. However, our method has the ability to measure stable estimates of [Formula: see text].


Geophysics ◽  
2016 ◽  
Vol 81 (4) ◽  
pp. M55-M69 ◽  
Author(s):  
Bastien Dupuy ◽  
Stéphane Garambois ◽  
Amir Asnaashari ◽  
Hadi M. Balhareth ◽  
Martin Landrø ◽  
...  

The estimation of quantitative rock physics properties is of great importance for reservoir characterization and monitoring in [Formula: see text] storage or enhanced oil recovery as an example. We have combined the high-resolution results of full-waveform inversion (FWI) methods with rock physics inversion. Because we consider a generic and dynamic rock physics model, our method is applicable to most kinds of rocks for a wide range of frequencies. The first step allows determination of viscoelastic effective properties, i.e., quantitative seismic attributes, whereas the rock physics inversion estimates rock physics properties (porosity, solid frame moduli, fluid phase properties, or saturation). This two-step workflow is applied to time-lapse synthetic and field cases. The sensitivity tests that we had previously carried out showed that it can be crucial to use multiparameter inputs to accurately recover fluid saturations and fluid properties. However, due to the limited data availability and difficulties in getting reliable multiparameter FWI results, we are limited to acoustic FWI results. The synthetic tests are conclusive even if they are favorable cases. For the first time-lapse fluid substitution synthetic case, we first characterize the rock frame parameters on the baseline model using P-wave velocity estimations obtained by acoustic FWI. Then, we obtain an accurate estimation of fluid bulk modulus from the time-lapse P-wave velocity. In the Marmousi synthetic case, the rock frame properties are accurately recovered for the baseline model, whereas the gas saturation change in the monitor model is not estimated correctly. On the field data example (time-lapse monitoring of an underground blowout in the North Sea), the estimation of rock frame properties gives results on a relatively narrow range, and we use this estimation as a starting model for the gas saturation inversion. We have found that the estimation of the gas saturation is not accurate enough, and the use of attenuation data is then required. However, the uncertainty on the estimation of baseline rock frame properties is not critical to monitor gas saturation changes.


2021 ◽  
Vol 114 (1) ◽  
Author(s):  
Alba Zappone ◽  
Eduard Kissling

AbstractThe Swiss Atlas of Physical Properties of Rocks (SAPHYR) project aims at centralize, uniform, and digitize dispersed and often hardly accessible laboratory data on physical properties of rocks from Switzerland and surrounding regions. The goal of SAPHYR is to make the quality-controlled and homogenized data digitally accessible to an open public, including industrial, engineering, land and resource planning companies as well as governmental and academic institutions, or simply common people interested in rock physics. The physical properties, derived from pre-existing literature or newly measured, are density, porosity and permeability as well as seismic, magnetic, thermal and electrical properties. The data were collected on samples either from outcrops or from tunnels and boreholes. At present, data from literature have been collected extensively for density, porosity, seismic and thermal properties. In the past years, effort has been placed especially on collecting samples and measuring the physical properties of rock types that were poorly documented in literature. A workflow for quality control on reliability and completeness of the data was established. We made the attempt to quantify the variability and the uncertainty of the data. The database has been recently transferred to the Federal Office of Topography swisstopo with the aim to develop the necessary tools to query the database and open it to the public. Laboratory measurements are continuously collected, therefore the database is ongoing and in continuous development. The spatial distribution of the physical properties can be visualized as maps using simple GIS tools. Here the distribution of bulk density and velocity at room conditions are presented as examples of data representation; the methodology to produce these maps is described in detail. Moreover we also present an exemplification of the use of specific datasets, for which pressure and temperatures derivatives are available, to develop crustal models.


Sign in / Sign up

Export Citation Format

Share Document