Embedding a Petroelastic Model in a Multipurpose Flow Simulator To Enhance the Value of 4D Seismic

2010 ◽  
Vol 13 (01) ◽  
pp. 37-43 ◽  
Author(s):  
John R. Fanchi

Summary Time-lapse (4D) seismic can be effectively integrated into the reservoir-management process by embedding the calculation of seismic attributes in a flow simulator. This paper describes a petroelastic model (PEM) embedded in a multipurpose flow simulator. The flow simulator may be used to model gas, black-oil, compositional, and thermal systems. The PEM can calculate reservoir geophysical attributes such as compressional-wave (P-wave) and shear-wave (S-wave) velocities and impedances, dynamic and static Young's moduli, and dynamic and static Poisson's ratios. Examples illustrate how to use the PEM to facilitate the integration of 4D seismic and reservoir flow modeling.

Geophysics ◽  
2012 ◽  
Vol 77 (6) ◽  
pp. T221-T237 ◽  
Author(s):  
Steven Shawn Smith ◽  
Ilya Tsvankin

Modeling of time shifts associated with time-lapse (4D) seismic surveys is helpful in evaluating reservoir depressurization and inverting for subsurface stress. Using coupled geomechanical and full-waveform seismic modeling, we study the influence of compaction-induced stress and strain around a simplified reservoir on compressional (P), shear (S), and mode-converted (PS) waves. We estimate compaction-induced time shifts and analyze their dependence on reflector depth and pressure drop inside the reservoir. Time shifts between synthetic baseline and monitor surveys are obtained by processing techniques that are potentially applicable to field data. Although P-wave time-shift lags for reflectors in the overburden are indicative of induced anisotropy, they are two to three times smaller than S-wave time-shift leads for reflectors beneath the reservoir. We also investigate the contributions of the deviatoric and volumetric stains to the time shifts for all three modes. Time shifts for S- and PS-waves are strongly influenced by elevated volumetric and deviatoric strains inside the reservoir. Almost constant S-wave time shifts for a range of offsets and source locations indicate that the contribution of stress-induced velocity anisotropy to shear-wave signatures is weak because the symmetry is close to elliptical. Our modeling also shows that mild tilt of a rectangular reservoir, or its replacement with an elliptically shaped reservoir of the same aspect ratio, has little influence on time shifts. Potentially, the developed methodology can be applied to estimate compaction-induced stress fields using simple compartmentalized reservoir models.


1990 ◽  
Vol 80 (2) ◽  
pp. 430-449 ◽  
Author(s):  
William Menke ◽  
Arthur L. Lerner-Lam ◽  
Bruce Dubendorff ◽  
Javier Pacheco

Abstract Except for its very onset, the P wave of earthquakes and chemical explosions observed at two narrow-aperture arrays on hard-rock sites in the Adirondack Mountains have a nearly random polarization. The amount of energy on the vertical, radial, and transverse components is about equal over the frequency range 5 to 30 Hz, for the entire seismogram. The spatial coherence of the seismograms is approximately exp(−cfΔx), where c is in the range 0.4 to 0.7 km−1Hz−1, f is frequency and Δx is the distance between array elements. Vertical, radial, and transverse components were quite coherent over the aperture of the array, indicating that the transverse motion of the compressional wave is a property of relatively large (106 m3) volumes of rock, and not just an anomaly caused by a malfunctioning instrument, poor instrument-rock coupling, or out-crop-scale effects. The spatial coherence is approximately independent of component, epicentral azimuth and range, and whether P- or S-wave coda is being considered, at least for propagation distances between 5 and 170 km. These results imply a strongly and three-dimensionally heterogeneous crust, with near-receiver scattering in the uppermost crust controlling the coherence properties of the waves.


Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. C107-C123
Author(s):  
Ivan Pšenčík ◽  
Bohuslav Růžek ◽  
Petr Jílek

We have developed a practical concept of compressional wave (P-wave) traveltime inversion in weakly to moderately anisotropic media of arbitrary symmetry and orientation. The concept provides sufficient freedom to explain and reproduce observed anisotropic seismic signatures to a high degree of accuracy. The key to this concept is the proposed P-wave anisotropy parameterization (A-parameters) that, together with the use of the weak-anisotropy approximation, leads to a significantly simplified theory. Here, as an example, we use a simple and transparent formula relating P-wave traveltimes to 15 P-wave A-parameters describing anisotropy of arbitrary symmetry. The formula is used in the inversion scheme, which does not require any a priori information about anisotropy symmetry and its orientation, and it is applicable to weak and moderate anisotropy. As the first step, we test applicability of the proposed scheme on a blind inversion of synthetic P-wave traveltimes generated in vertical seismic profile experiments in homogeneous models. Three models of varying anisotropy are used: tilted orthorhombic and triclinic models of moderate anisotropy (approximately 10%) and an orthorhombic model of strong anisotropy (>25%) with a horizontal plane of symmetry. In all cases, the inversion yields the complete set of 15 P-wave A-parameters, which make reconstruction of corresponding phase-velocity surfaces possible with high accuracy. The inversion scheme is robust with respect to noise and the source distribution pattern. Its quality depends on the angular illumination of the medium; we determine how the absence of nearly horizontal propagation directions affects inversion accuracy. The results of the inversion are applicable, for example, in migration or as a starting model for inversion methods, such as full-waveform inversion, if a model refinement is desired. A similar procedure could be designed for the inversion of S-wave traveltimes in anisotropic media of arbitrary symmetry.


2020 ◽  
Author(s):  
Hanneke Paulssen ◽  
Wen Zhou

<p>Between 2013 and 2017, the Groningen gas field was monitored by several deployments of an array of geophones in a deep borehole at reservoir level (3 km). Zhou & Paulssen (2017) showed that the P- and S-velocity structure of the reservoir could be retrieved from noise interferometry by cross-correlation. Here we show that deconvolution interferometry of high-frequency train signals from a nearby railroad not only allows determination of the velocity structure with higher accuracy, but also enables time-lapse measurements. We found that the travel times within the reservoir decrease by a few tens of microseconds for two 5-month periods. The observed travel time decreases are associated to velocity increases caused by compaction of the reservoir. However, the uncertainties are relatively large. <br>Striking is the large P-wave travel time anomaly (-0.8 ms) during a distinct period of time (17 Jul - 2 Sep 2015). It is only observed for inter-geophone paths that cross the gas-water contact (GWC) of the reservoir. The anomaly started 4 days after drilling into the reservoir of a new well at 4.5 km distance and ended 4 days after the drilling operations stopped. We did not find an associated S-wave travel time anomaly. This suggests that the anomaly is caused by a temporary elevation of the GWC (water replacing gas) of approximately 20 m. We suggest that the GWC is elevated due to pore-pressure variations during drilling. The 4-day delay corresponds to a pore-pressure diffusivity of ~5m<sup>2</sup>/s, which is in good agreement with the value found from material parameters and the diffusivity of (induced) seismicity for various regions in the world. </p>


Geophysics ◽  
2020 ◽  
Vol 85 (5) ◽  
pp. A31-A36
Author(s):  
Qizhen Du ◽  
Qiang Zhao ◽  
Qingqing Li ◽  
Liyun Fu ◽  
Qifeng Sun

Methods to decompose the elastic wavefield into compressional wave (P-wave) and shear wave (S-wave) components in heterogeneous media without wavefield distortions or energy leakage are the key issues in elastic imaging and inversion. We have introduced a decoupled P- and S-wave propagator to form an efficient elastic reverse time migration (RTM) framework, without assuming homogeneous Lamé parameters. Also, no wave-mode conversions occur using the proposed propagator in the presence of strong heterogeneities, which avoids the potential imaging artifacts caused by wave-mode conversions in the receiver-side backward extrapolation. In the proposed elastic RTM framework, the source-side forward wavefield is simulated with a P-wave propagator. The receiver-side wavefield is back extrapolated with the proposed propagator, using the recorded multicomponent seismic data as input. Compared to the conventional elastic RTM, the proposed framework reduces the computational complexity while preserving the imaging accuracy. We have determined its accuracy and efficiency using two synthetic examples.


2000 ◽  
Vol 3 (01) ◽  
pp. 88-97 ◽  
Author(s):  
R.D. Benson ◽  
T.L. Davis

Summary This article presents the results of a multidisciplinary, four-dimensional (4D) (time-lapse), three-component (3C) (multicomponent) seismic study of a CO2 injection project in vacuum field, New Mexico. The ability to sense bulk rock/fluid properties with 4D, 3C seismology enables characterization of the most important transport property of a reservoir, namely, permeability. Because of the high volume resolution of the 4D, 3C seismology, we can monitor the sweep efficiency of a production process to see if reserves are bypassed by channeling around lower permeability parts of the reservoir and the rate at which the channeling occurs. In doing so, we can change production processes to sweep the reservoir more efficiently. Introduction Improving reservoir performance and enhancing hydrocarbon recovery while reducing environmental impact are critical to the future of the petroleum industry. To do this, it must be possible to characterize reservoir parameters including fluid properties and their changes with time, i.e., dynamic reservoir characterization. The objectives of our research arerepeated acquisition of a three-dimensional (3D), three-component (3C) seismic survey;demonstrate the ability of 3D, 3C, and four-dimensional (4D), 3C seismology to detect and monitor rock/fluid property change associated with a production process;incorporate geological, petrophysical, petroleum engineering, and other geophysical studies;refine the reservoir model and recommend procedures for scaling up from a pilot injection program to partial field flood to achieve maximum sweep efficiency and minimize bypassed reservoir zones;link bulk rock/fluid property variation monitored by time-lapse multicomponent (4D, 3C) seismic surveying to dynamic attributes of the reservoir including permeability, fluids, and flow characterization. Three-dimensional, 3C seismology involves seismic data acquisition in three orientations at each receiver location—two orthogonal horizontal and one vertical. When three source components are used, nine times the amount of data of a conventional P-wave 3D survey can be recorded. Horizontal components of source and receiver displacements enable shear- (S-) wave recording; this is a powerful complement to vertical P-wave recording. Three-dimensional, 3C seismic surveys provide significantly more information about the rock/fluid properties of a reservoir than can be achieved from conventional P-wave seismic surveys alone. By combining P- and S-wave recording, the seismic ability to determine rock/fluid property changes in the subsurface is increased. Seismic wave propagation includes travel time, reflectivity, and the effects of anisotropy and attenuation. In-situ stress orientation and relative magnitudes can be derived from seismic anisotropy. Rock/fluid properties, including lithology and porosity, may be obtained from comparative travel times or velocities of P and S waves. Other rock/fluid properties, including permeability, may be determined from comparative P and S anisotropy, travel time, reflectivity, and attenuation measurements. By combining P- and S-wave recording, seismic wave propagation characteristics can be transformed into reservoir parameters. Introducing time as the "fourth dimension," new time-lapse (4D), 3C seismology is a tool to monitor production processes and to determine reservoir property variations under changing conditions. Using 4D, 3C seismic monitoring as an integral part of dynamic reservoir characterization, refinements can be made to production processes to improve reservoir hydrocarbon recovery. VP/VS ratios for both the fast S1 shear component and slow S2 shear component may provide a tool for separating bulk rock changes due to fluid property variations from bulk rock changes due to effective stress variations. Changes in shear wave anisotropy may reflect varying concentrations of open fractures and low aspect ratio pore structure in both a spatial and temporal sense across the reservoir. The permeability of a formation, or the connectivity of the pore space, will be the target in 4D, 3C seismology. Refinements made to reservoir characterization techniques and their applications, now extending into the fourth dimension, are an important new area of research. Benefits of this research will include improved reservoir characterization and correlative increased hydrocarbon recovery and reduction in operating costs through improved reservoir management. Geologic Setting Since early Permian time, the general evolution of the portion of the Permian Basin which includes vacuum field is that of a progressively shallowing-upward carbonate platform. Aggrading and prograding cycles represent, respectively, the results of high stand and still stand sea levels. At the shelf edge these platform carbonates typically grade into reef-type deposits such as the Abo, Goat Seep, and Capitan formations. The San Andres is an exception; no reef-like rocks have been detected. Beyond the shelf edge, in the Delaware basin, clastic rocks, especially siliciclastics, were deposited during a lowstand sea level. Vacuum field is located on a large anticlinal structure that plunges slightly to the east-northeast. The San Andres and Grayburg formations correspond to the rim of a broad carbonate shelf province to the north and northwest, northwest shelf, and of a deeper intracratonic basin, Delaware basin, on the southeast and east.1 The overall area including the Midland basin, northern and eastern shelves, and central basin platform are part of a major restricted intracratonic basin which existed during Permian time. West Texas and southeast New Mexico were in the low latitudes throughout the late Paleozoic period, making them an ideal location for carbonate sedimentation. As a consequence of this tropical environment, broad carbonate shelves were established on the margins of the Delaware and Midland basins.2


2021 ◽  
Vol 22 (3) ◽  
pp. 1-9
Author(s):  
Qahtan Abdul Aziz ◽  
Hassan Abdul Hussein

The Compressional-wave (Vp) data are useful for reservoir exploration, drilling operations, stimulation, hydraulic fracturing employment, and development plans for a specific reservoir. Due to the different nature and behavior of the influencing parameters, more complex nonlinearity exists for Vp modeling purposes. In this study, a statistical relationship between compressional wave velocity and petrophysical parameters was developed from wireline log data for Jeribe formation in Fauqi oil field south Est Iraq, which is studied using single and multiple linear regressions. The model concentrated on predicting compressional wave velocity from petrophysical parameters and any pair of shear waves velocity, porosity, density, and fluid saturation in carbonate rocks. A strong linear correlation between P-wave velocity and S-wave velocity and between P-wave velocity and density rock was found. The resulting linear equations can be used to estimate P-wave velocity from the S-wave velocity in the case of both. The results of multiple regression analysis indicated that the density, porosity, water-saturated, and shear wave velocity (VS) are strongly related to Vp.


Geophysics ◽  
1995 ◽  
Vol 60 (3) ◽  
pp. 631-650 ◽  
Author(s):  
Mark E. Mathisen ◽  
Anthony A. Vasiliou ◽  
Paul Cunningham ◽  
J. Shaw ◽  
J. H. Justice ◽  
...  

Time‐lapse crosswell seismic data acquired with a cemented receiver cable have been processed into P‐ and S‐wave tomograms which image heavy oil sand lithofacies and changes as a result of steam injection. Twenty‐seven crosswell surveys were acquired between two wells over a 3.5 month period before, during, and after a 34‐day, 30 MBBL [Formula: see text] steam injection cycle. Interpretation was based on correlations with reservoir data and models, observation well data, and engineering documentation of the production history and steam cycle. Baseline S‐ and P‐wave tomograms image reservoir sand flow units and areas affected by past cyclic steam injection. S‐wave tomograms define lithology and porosity contrasts between the excellent reservoir quality, “high flow” turbidite channel facies and the interbedded “low to moderate flow” bioturbated levee facies. The reservoir dip of approximately 20° is defined by the velocity contrast between lithofacies. P‐wave baseline tomograms image lithology, porosity, structure, and several low velocity zones caused by past steam injection. Previous steam‐heat injection caused the formation of gas which reduced velocities as much as several thousand ft/s (600 m/s), an amount which obscures the velocity contrast between lithofacies and smaller velocity reductions as a result of temperature alone. Time‐lapse and difference P‐wave tomograms document several areas with small decreases in velocity during steam injection and larger decreases after cyclic steam injection. Velocity reductions range from 300 to 900 ft/s (90 to 270 m/s) adjacent to and above injectors located 20 to 50 feet (6 to 15 m) from the tomogram cross‐section. Poisson’s ratio tomograms show a significant decrease (.10) in the same area, and include low values indicative of gas saturation. Continuous injectors located 50 to 350 feet (15 to 100 m) from the survey area also caused a progressive decrease in velocity of the “high flow” channel sands during the time‐lapse survey. Interdisciplinary interpretation indicates that tomograms not only complement other borehole‐derived reservoir characterization and temperature monitoring data but can be used to quantitatively characterize interwell reservoir properties and monitor changes as a result of the thermal recovery process. Monitoring results over 3.5 months confirms that stratification has controlled the flow of steam, in contrast to gravity override. This suggests that tomographic images of reservoir flow‐units and gas‐bearing high temperature zones should be useful for positioning wells and optimizing injection intervals, steam volumes, and producing well completions.


2011 ◽  
Vol 2011 ◽  
pp. 1-16 ◽  
Author(s):  
Paritosh Singh ◽  
Thomas Davis

The Upper Morrow sandstones in the western Anadarko Basin have been prolific oil producers for more than five decades. Detection of Morrow sandstones is a major problem in the exploration of new fields and the characterization of existing fields because they are often very thin and laterally discontinuous. Until recently compressional wave data have been the primary resource for mapping the lateral extent of Morrow sandstones. The success with compressional wave datasets is limited because the acoustic impedance contrast between the reservoir sandstones and the encasing shales is small. Here, we have performed full waveform modeling study to understand the Morrow sandstone signatures on compressional wave (P-wave), converted-wave (PS-wave) and pure shear wave (S-wave) gathers. The contrast in rigidity between the Morrow sandstone and surrounding shale causes a strong seismic expression on the S-wave data. Morrow sandstone shows a distinct high amplitude event in pure S-wave modeled gathers as compared to the weaker P- and PS-wave events. Modeling also helps in understanding the adverse effect of interbed multiples (due to shallow high velocity anhydrite layers) and side lobe interference effects at the Morrow level. Modeling tied with the field data demonstrates that S-waves are more robust than P-waves in detecting the Morrow sandstone reservoirs.


SPE Journal ◽  
2021 ◽  
pp. 1-10
Author(s):  
Jing Fu ◽  
Carl Sondergeld ◽  
Chandra Rai

Summary Elastic wave velocities are commonly used to predict porosity, mineralogy, and lithology from formation properties. When only P-wave sonics are available in historical wells, systematics for predicting shear velocities are useful for developing elastic models. Although much research has been done on conventional reservoir velocity systematics, the equivalency for unconventional formations is still a work in progress. There has also been a limited number of research studies with laboratory measures published. Using laboratory pulse transmission ultrasonic data, we created a Vp-Vs systematic for the Meramec Formation in this study. The effects of porosity and mineralogy on velocities are explored, as well as a comparison of Meramec velocity systematics with well-established literature systematics. Vp and Vs measurements were taken on 385 dodecane-saturated core samples from seven Meramec wells (106 vertical and 279 horizontal plugs). S-wave and P-wave anisotropy in Meramec Formation samples used in this study are typically less than 10%. Each sample was also tested for porosity and mineralogy. We find that velocities are more sensitive to porosity than mineralogy by a factor of 10. Below are our equations for predicting Vp and Vs (in km/s), when only clay content and porosity are known. In these equations, φ is the volume fraction pores, and Clays is the weight fraction of clay. These equations are for those samples in which there is low P-wave and S-wave anisotropies:(1)Vp=6.4−1.2*Clays−15.4*φ(R2=0.5),(2)Vs=3.6−0.5*Clays−5.2*φ(R2=0.4). We suggest two methods for calculating Vs from Vp: Ignoring anisotropy, we combined both Vp and Vs measurements from all vertical plugs and low anisotropy horizontal plugs to create a single shear wave predictor; and considering anisotropy, Vp measurements from horizontal plugs were corrected using Thomsen’s compressional wave anisotropy parameter, after which a shear velocity predictor was generated. The shear wave predictors for dodecane-saturated measurements are as follows (all velocities are km/s):(3)Method 1: Vs= 0.90 + 0.42*Vp (R2=0.7),(4)Method 2: Vs= 0.80 + 0.45*Vp (R2=0.6). The residual and estimated error in Eq. 3 is slightly less than in Eq. 4. Even though there is a significant variance in measurement frequency, the Meramec velocity systematic shows good agreement with dipole wireline measurements using the first equation. The Meramec velocity systematics differ significantly from previously published systematics, such as the trend line by Greenberg and Castagna (1992) and the shale trend line by Vernik et al. (2018). Using the correlations by Greenberg and Castagna (1992) for limestone or dolomite, the shear velocities of the samples in this study cannot be predicted. These data have yielded shear wave systematics, which can be used in wireline and seismic investigations. The results suggest that the method of ignoring anisotropy yields a better Vs estimate than the one that takes anisotropy into account. Using well-established shear wave velocity systematics from the published literature can result in an estimated inaccuracy of greater than 16%. It is important to calibrate velocity systematics to the target formation.


Sign in / Sign up

Export Citation Format

Share Document