Time-Lapse Seismic Monitoring and Dynamic Reservoir Characterization, Central Vacuum Unit, Lea County, New Mexico

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

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.


2022 ◽  
Vol 41 (1) ◽  
pp. 47-53
Author(s):  
Zhiwen Deng ◽  
Rui Zhang ◽  
Liang Gou ◽  
Shaohua Zhang ◽  
Yuanyuan Yue ◽  
...  

The formation containing shallow gas clouds poses a major challenge for conventional P-wave seismic surveys in the Sanhu area, Qaidam Basin, west China, as it dramatically attenuates seismic P-waves, resulting in high uncertainty in the subsurface structure and complexity in reservoir characterization. To address this issue, we proposed a workflow of direct shear-wave seismic (S-S) surveys. This is because the shear wave is not significantly affected by the pore fluid. Our workflow includes acquisition, processing, and interpretation in calibration with conventional P-wave seismic data to obtain improved subsurface structure images and reservoir characterization. To procure a good S-wave seismic image, several key techniques were applied: (1) a newly developed S-wave vibrator, one of the most powerful such vibrators in the world, was used to send a strong S-wave into the subsurface; (2) the acquired 9C S-S data sets initially were rotated into SH-SH and SV-SV components and subsequently were rotated into fast and slow S-wave components; and (3) a surface-wave inversion technique was applied to obtain the near-surface shear-wave velocity, used for static correction. As expected, the S-wave data were not affected by the gas clouds. This allowed us to map the subsurface structures with stronger confidence than with the P-wave data. Such S-wave data materialize into similar frequency spectra as P-wave data with a better signal-to-noise ratio. Seismic attributes were also applied to the S-wave data sets. This resulted in clearly visible geologic features that were invisible in the P-wave data.


1974 ◽  
Vol 64 (4) ◽  
pp. 1159-1180 ◽  
Author(s):  
F. A. Dahlen

abstract We construct a theoretical three-dimensional kinematical model of shallow-focus earthquake faulting in order to investigate the ratio of the P- and S-wave corner frequencies of the far-field elastic radiation. We attempt to incorporate in this model all of the important gross kinematical features which would arise if ordinary mechanical friction should be the dominant traction resisting fault motion. These features include a self-similar nucleation at a single point, a subsonic spreading of rupture away from that point, and a termination of faulting by smooth deceleration. We show that the ratio of the P-wave corner frequency to the S-wave corner frequency for any model which has these features will be less than unity at all points on the focal sphere.


2019 ◽  
Vol 23 (3) ◽  
pp. 209-223 ◽  
Author(s):  
Caglar Ozer ◽  
Mehmet Ozyazicioglu

Erzurum and its surroundings are one of the seismically active and hydrothermal areas in the Eastern part of Turkey. This study is the first approach to characterize the crust by seismic features by using the local earthquake tomography method. The earthquake source location and the three dimensional seismic velocity structures are solved simultaneously by an iterative tomographic algorithm, LOTOS-12. Data from a combined permanent network comprising comprises of 59 seismometers which was installed by Ataturk University-Earthquake Research Center and Earthquake Department of the Disaster and Emergency Management Authority  to monitor the seismic activity in the Eastern Anatolia, In this paper, three-dimensional Vp and Vp/Vs characteristics of Erzurum geothermal area were investigated down to 30 km by using 1685 well-located earthquakes with 29.894 arrival times, consisting of 17.298 P- wave and 12.596 S- wave arrivals. We develop new high-resolution depth-cross sections through Erzurum and its surroundings to provide the subsurface geological structure of seismogenic layers and geothermal areas. We applied various size horizontal and vertical checkerboard resolution tests to determine the quality of our inversion process. The basin models are traceable down to 3 km depth, in terms of P-wave velocity models. The higher P-wave velocity areas in surface layers are related to the metamorphic and magmatic compact materials. We report that the low Vp and high Vp/Vs values are observed in Yedisu, Kaynarpinar, Askale, Cimenozu, Kaplica, Ovacik, Yigitler, E part of Icmeler, Koprukoy, Uzunahmet, Budakli, Soylemez, Koprukoy, Gunduzu, Karayazi, Icmesu, E part of Horasan and Kaynak regions indicated geothermal reservoir.


Geophysics ◽  
2010 ◽  
Vol 75 (5) ◽  
pp. 75A15-75A29 ◽  
Author(s):  
Ilya Tsvankin ◽  
James Gaiser ◽  
Vladimir Grechka ◽  
Mirko van der Baan ◽  
Leon Thomsen

Recent advances in parameter estimation and seismic processing have allowed incorporation of anisotropic models into a wide range of seismic methods. In particular, vertical and tilted transverse isotropy are currently treated as an integral part of velocity fields employed in prestack depth migration algorithms, especially those based on the wave equation. We briefly review the state of the art in modeling, processing, and inversion of seismic data for anisotropic media. Topics include optimal parameterization, body-wave modeling methods, P-wave velocity analysis and imaging, processing in the [Formula: see text] domain, anisotropy estimation from vertical-seismic-profiling (VSP) surveys, moveout inversion of wide-azimuth data, amplitude-variation-with-offset (AVO) analysis, processing and applications of shear and mode-converted waves, and fracture characterization. When outlining future trends in anisotropy studies, we emphasize that continued progress in data-acquisition technology is likely to spur transition from transverse isotropy to lower anisotropic symmetries (e.g., orthorhombic). Further development of inversion and processing methods for such realistic anisotropic models should facilitate effective application of anisotropy parameters in lithology discrimination, fracture detection, and time-lapse seismology.


2000 ◽  
Vol 31 ◽  
pp. 85-90 ◽  
Author(s):  
N. Deichmann ◽  
J. Ansorge ◽  
F. Scherbaum ◽  
A. Aschwanden ◽  
F. Bernard ◽  
...  

AbstractTo obtain more reliable information about the focal-depth distribution of icequakes, in April 1997 we operated an array of seven portable digital seismographs on Unteraargletscher, central Swiss Alps. Over 5000 events were detected by at least two instruments during the 9 day recording period. P-wave velocities (3770 m f) were determined from several calibration shots detonated at the glacier surface as well as in a 49 m deep borehole, whereas S-wave velocities (1860 ms–1) were derived from a simultaneous inversion for Vp/Vs6 applied to 169 icequakes. So far, hypocentral locations have been calculated for over 300 icequakes. Besides confirming the occurrence of shallow events associated with the opening of crevasses, our results show that a small but significant fraction of the hypocenters are located at or near the glacier bed. One event was found at an intermediate depth of about 120 m. Three-dimensional particle-motion diagrams of both explosions and icequakes clearly demonstrate that all vertical component seismograms from shallow sources are dominated by the Rayleigh wave. On the other hand, for events occurring at depths greater than about 40 m, the Rayleigh wave disappears almost entirely. Therefore, a qualitative analysis of the signal character provides direct information on the focal depth of an event and was used as an independent check of the locations obtained from traditional arrival-time inversions. Thus, our results demonstrate that deep icequakes do occur and that simple rheological models, according to which brittle deformation is restricted to the uppermost part of a glacier, may need revision.


1974 ◽  
Vol 64 (6) ◽  
pp. 1789-1808 ◽  
Author(s):  
Robert Burridge ◽  
Cesar Levy

abstract It has recently been shown (Burridge, 1973) that in two dimensions plane-strain shear cracks lacking cohesion may run at speeds ranging from the Rayleigh-wave to the S-wave speed or possibly at the P-wave speed. On the other hand, it has long been known that in antiplane strain, cracks lacking cohesion must run at least at the S-wave speed. Since locally at the edge of a three-dimensional crack there is a combination of plane and antiplane strain, we have sought and found solutions for circular shear cracks expanding at the S-wave speed and at the P-wave speed. These have finite shear tractions ahead of the crack and so may correspond to frictional sliding in the absence of cohesion. The analysis combines the method of Kostrov (1964b) with the results of Burridge (1973). We carry out a complete evaluation for the displacements and tractions on the fault plane, and far-field radiation for the S-wave-speed crack. The corresponding evaluations for the P-wave speed are not elementary and are not attempted here. As far as the authors are aware, these are the first analytic solutions of three-dimensional crack problems which satisfy a physically plausible fracture criterion for failure under shear.


Geophysics ◽  
1992 ◽  
Vol 57 (4) ◽  
pp. 643-646
Author(s):  
Hans A. K. Edelmann

If shear waves are to be recorded, all other types of waves (including P waves) have to be regarded as noise. All data processing applied later is limited in its success, not so much by the character of the signal, but by the character of the noise superimposed on the signal. Therefore an optimum method for simultaneous P‐ and S‐wave recording does not exist per se. All efforts made in the field that help to enhance the relatively weak S‐wave signal enhance the possibility of a more detailed interpretation such as polarization analysis. In the course of shear‐wave investigations over a period of more than ten years, simultaneous P‐ and SV‐wave recording has yielded fairly good results for velocity ratio determination, but has never produced satisfying results for polarization analysis because of the interfering P‐wave events. When generating pure SH‐waves, however, P‐wave arrival amplitudes in a shot record can, under favorable conditions, be kept well below the SH‐wave amplitude (−40 dB). Through additional processing, a ratio of P‐ to SH‐signal amplitude of −60 dB can be reached. The improvement achieved by making separate shear‐wave recordings, obviously, must be weighed against the additional costs caused by these recordings.


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.


Sign in / Sign up

Export Citation Format

Share Document