In mature oil fields undergoing enhanced oil recovery methods, such as CO2 injection, monitoring the reservoir changes becomes important. To understand how reservoir changes influence compressional wave (P) and shear wave (S) velocities, we conducted laboratory core experiments on five core samples taken from the Morrow A sandstone at Postle Field, Oklahoma. The laboratory experiments measured P- and S-wave velocities as a function of confining pressure, pore pressure, and fluid type (which included CO2 in the gas and supercritical phase). P-wave velocity shows a response that is sensitive to both pore pressure and fluid saturation. However, S-wave velocity is primarily sensitive to changes in pore pressure. We use the fluid and pore pressure response measured from the core samples to modify velocity well logs through a log facies model correlation. The modified well logs simulate the brine- and CO2-saturated cases at minimum and maximum reservoir pressure and are inputs for full waveform seismic modeling. Modeling shows how P- and S-waves have a different time-lapse amplitude response with offset. The results from the laboratory experiments and modeling show the advantages of combining P- and S-wave attributes in recognizing the mechanism responsible for time-lapse changes due to CO2 injection.

1. Introduction

Currently, there are 128 enhanced oil recovery (EOR) projects worldwide using CO2 injection as a tertiary recovery method, with 113 taking place in the United States [1]. Often the criterion for successful seismic monitoring of CO2 projects is the ability to collect high-resolution seismic data, since seismic methods have shown to be very effective in detecting time-lapse changes [26]. The effectiveness of seismic monitoring depends upon the changes in acoustic and elastic impedance caused by injection and production. These changes are then observed in the amplitude and travel time of seismic waves.

Time-lapse seismic methods utilizing compressional waves (P-wave) and shear waves (S-wave) can ideally improve the identification of flooded or pressurized compartments within the reservoir during a CO2 EOR operation. The use of P-wave data is not a definitive method for pressure detection, but the addition of S-wave data may provide a distinction between a gas-bearing zone and a zone of high pore pressure [7], this is a motivation for using S-waves as a time-lapse tool. An example provided by Xue et al. [8] explains if an observed time-lapse P-wave anomaly is caused by a decrease in velocity without an associated S-wave anomaly, the zone has most likely been flushed with CO2 without any changes in pore pressure. However, if the S-wave data show a similar anomaly, due to a reduction in velocity, the zone may have undergone a pore pressure buildup in addition to being flushed with CO2. Identifying the impact of changing stresses within the reservoir is as important as recognizing the effects of saturation.

Saturation effects are typically estimated from fluid substitution modeling (e.g., Gassmann [9] and Kuster and Toksoz [10]). Since fluid substitution methods do not directly account for changing rock properties with respect to varying pore pressure, additional information must be used to account for the effects of pore pressure. This additional information may be obtained from, for example, reservoir flow simulation or experimental core studies. The purpose of this paper is to show how the results from rock physics experiments that include saturation and pore pressure variation can be used to estimate amplitude changes in a reservoir undergoing CO2 injection.

The contribution of this study applies directly to a major ongoing research project at Colorado School of Mines. The Reservoir Characterization Project (RCP) is a consortium at Colorado School of Mines whose research goal is to monitor CO2 EOR projects to effectively manage and enhance the recovery process. The RCP has used P-, S-, and converted-wave (PS) time-lapse seismic acquisition to monitor CO2 floods for EOR. The RCP conducted the first ever time-lapse, multicomponent seismic survey in conjunction with a CO2 huff-n-puff project at Vacuum Field, New Mexico. Multicomponent time-lapse seismic data was used to monitor the EOR project at Weyburn Field, Saskatchewan, Canada. The current research project is a time-lapse study of an incised valley fill system undergoing a water alternating gas (i.e., CO2) injection scheme for EOR at Postle Field, Texas County, Oklahoma, located in the Anadarko Basin. The results from the experiment and analysis reflect the use of reservoir core samples and actual reservoir fluids. The aim of this work is to provide insight related directly to this specific type of environment.

1.1. Previous Work

Studies on the effects of saturation, pressure, and pore fluid have been documented as early as 1966 [12]. King [12] reported on the change in velocities of five sandstone samples saturated with air, an NaCl solution, and kerosene as a function of hydrostatic confining pressure and internal pore pressure. Furthermore, Domenico [13] showed that S-wave velocity is more sensitive to increasing differential pressure than P-wave velocity in three sandstone cores of low, intermediate, and high porosity.

A series of papers have been published on the results of the P-wave velocity response to the migration of CO2 in water-saturated sandstones using an experimental setup that simulates cross-well seismic profiling. Xue et al. [8] mapped the velocity changes caused by CO2 injection in water-saturated Tako and Shirahama sandstones at constant pore pressure. In a similar experiment, Xue and Ohsumi [14] mapped the velocity changes caused by the injection of gas, liquid, and supercritical phase CO2 into a water-saturated Tako sandstone. They investigated the changes in velocity with increasing pore pressure and concluded that CO2 injection causes larger velocity changes than the change associated with increasing pore pressure in their core samples. In addition to the P-wave velocity, strain was measured as a function of hydrostatic pressure in a dry sample. Their observations showed how strain normal and parallel to the bedding plane are not equal and more deformation occurred normal to the bedding plane. To improve upon their cross-well seismic profiling experiment and the accuracy of the P-wave velocity changes relative to the CO2 injection, Xue and Lei [15] applied a difference tomography method to the measured slowness. This technique correlated areas of higher porosity with greater velocity reduction. Velocity reduction when CO2 was in liquid or supercritical phase was twice as much as when in the gas phase. From the same dataset as Xue and Lei [15], Shi et al. [16] focused on P-wave velocity reductions due to supercritical CO2 displacement of pore water, where the velocity reductions deviated significantly from the Gassmann-predicted velocities for both patchy and uniform saturations. Lei and Xue [17] expanded on the results of Xue and Lei [15] and Shi et al. [16] and showed that P-wave velocity decreased, and the attenuation coefficient increased, due to gaseous, liquid, and supercritical CO2 partially replacing pore water. Kim et al. [18] also investigated the use of P-wave velocity and resistivity as a way to monitor saturation of a sandstone during a CO2 injection process. This was accomplished by injecting supercritical CO2 into a water-saturated Berea sandstone sample at simulated reservoir conditions.

The effects of saturation, pore pressure, and effective stress, studied on sandstone core samples from the Northwest Shelf of Australia, show that the / ratio is an indicator for the saturation state between oil-saturated and the dry state [19]. Siggins [20] also studied the response of CO2-saturated sandstones at in-situ reservoir conditions for two synthetic sandstones and one reservoir sandstone/mudstone from the Otway Basin, C seam, Waarre Formation. His findings were both P-wave and S-wave velocities decreased when the samples were saturated with liquid CO2, compared to dry rock. Also, the synthetic sandstones showed better agreement with Gassman’s fluid substitution theory than did the reservoir sandstone. Similarly, Siggins et al. [21] compared synthetic sandstones and Waarre C formation reservoir sandstones when dry and saturated with gaseous and liquid-phase CO2. Larger changes were observed in both P-wave and S-wave velocity when CO2 was in the liquid phase.

The effects of saturation and pore pressure have also been studied with respect to the effective pressure law and the effective pressure coefficient [22, 23] and how fractures change the constituency of the rock and the effective pressure coefficient [24]. The work by Xu et al. [24] on effective pressure shows the importance of why assuming the effective pressure coefficient to be unity could lead to significant errors when predicting time-lapse changes.

The RCP has used experimental rock physics to understand how the changes in confining pressure, pore pressure, and fluid saturation affect P-wave and S-wave velocities. This knowledge has then been applied to the interpretation of multicomponent time-lapse surface seismic monitoring. Capello [25] used rock physics as a feasibility study for seismic monitoring of a CO2 huff-n-puff flood in the Permian San Andres Formation, a dolomitized carbonate deposit at Vacuum Field. Also using samples from the San Andres Formation, Duranti [26] investigated the coexistence of equant pores and fractures and found that a dual porosity rock physics model was more successful in predicting the changes in differential pressure and saturation observed in the experimental data.

Brown [27] integrated experimental rock physics measurements from the carbonate Marly beds of the Mississippian Charles Formation to improve the interpretation of time-lapse seismic data for the monitoring of a CO2 flood at Weyburn Field. Yamamoto [28] included the rock physics results obtained by Brown [27] to improve the flow modeling of CO2 by minimizing the difference between calculated and observed acoustic impedance while simultaneously matching the production history.

Lastly, Rojas [29] measured dry rock velocities for the purpose of analyzing the different responses from Gassmann fluid substitution. This study on tight gas sands linked / sensitivity to pore pressure changes, lithology, and fluid content to aid in identifying zones with a higher reservoir potential of the Late Cretaceous Williams Fork Formation at Rulison Field, CO.

Perhaps the earliest and most cited laboratory experiment related to monitoring a CO2 EOR operation is that by Wang and Nur [30]. They measured ultrasonic P-wave and S-wave velocities on seven sandstones and one unconsolidated sand saturated with n-hexadecane and then flushed with CO2. Their measurements were done at constant confining pressure over a range of pore pressures and temperatures above and below the critical point of CO2. Their findings show a significant decrease in P-wave velocity which was attributed to CO2 saturation. The S-wave velocity exhibited a greater effect to pore pressure than saturation and was shown to be less sensitive to saturation.

Outside of research done on the Tako, Shirahama, Waarre C Formation, Williams Fork Formation sandstones and the carbonates of San Andres and Charles Formation, most of the experimental studies are not done on reservoir samples obtained from well core. This is evident in the statement made by Siggins et al. [21], “literature on laboratory measurements of the effects of CO2 on the seismic properties of reservoir rocks is relatively sparse.”

2. Postle Field Overview

The geologic environment of the producing reservoir at Postle Field is the Upper Morrow sandstone Formation, which is a Pennsylvanian-age siliciclastic incised valley fill system. Typical Morrow sediments consist of shales punctuated by valley fill sand deposits [31]. The Morrow Formation is encased between two high-stand limestones [32]. The upper boundary is a conformable surface at the base of the Atoka Thirteen Fingers Limestone. The lower boundary of the Morrow Formation is an unconformable surface above the Mississippian Chester Limestone. The deposition of Postle Field during Pennsylvanian time was a relative sea-level low-stand shoreline where the Morrow Formation is much sandier [33]. The primary producing unit in the Morrow Formation is the Upper Morrow A, A1, and A2 sands. In Postle Field, the producing formation that is undergoing the CO2 flood is the Morrow A sandstone.

Postle Field is a large mature oil field undergoing a water alternating gas (WAG) flood for enhanced oil recovery of the Morrow A sandstone. The volume of water to CO2 injected into the reservoir is referred to as the WAG ratio and is expressed in reservoir volumes at reservoir conditions. The average WAG ratio implemented at Postle Field is 0.35 to 1. The reservoir pressure is maintained above the minimum miscibility pressure of 14.48 MPa by injecting water and CO2 at 29.65 MPa and by keeping the bottom hole pressure of a producing well higher than 6.89 MPa [34].

3. Sample Selection

We retrieved five core plug samples from three different lithological zones within the Morrow A sandstone interval. The three zones were a low-porosity and low-permeability cemented zone, a higher-porosity and high-permeability zone, and a higher-porosity and lower-permeability zone. Each sample was cut parallel to the bedding plane. The cylindrical dimensions of the samples were roughly 3.81 cm in length and 2.54 cm in diameter. The depths from where the core samples were retrieved, as well as density, porosity, and permeability are listed in Table 1.

3.1. Mineralogy

We used QEMSCAN (quantitative evaluation of minerals by scanning electron microscopy) to quantitatively evaluate mineral content on each of the three lithological zones. The cemented zone (Figure 1) consists of subrounded quartz grains cemented by ankerite and kaolinite clay. This cemented zone shows a lower amount of grain-to-grain contact which suggests the possibility of early cementation. The lower amount of grain-to-grain contact is likely related to the presence of ankerite which may indicate that the pore water was rich in iron. In addition, there are occurrences of angular overgrowths over subrounded quartz grains. Even though the QEMSCAN measured 14% porosity, most of the porosity for this cemented zone is interstitial between the fine-grained kaolinite clay.

The higher-porosity and high-permeability zone (Figure 2) consists mostly of subrounded quartz grains with lesser amounts of kaolinite clay cement, where most of the porosity (25%) occurs between quartz grains. Unlike the other zones, only a trace amount of Ankerite occurs as cement.

The higher-porosity and lower-permeability zone (Figure 3) consists mostly of subrounded quartz grains with lesser amounts of ankerite and kaolinite clay cement. This zone is more similar to the cemented zone with respect to the presence of ankerite and interstitial porosity. Much like the cemented zone, most of the porosity (20%) is interstitial between the fine-grained kaolinite and chlorite.

Table 2 shows the mineral volume percentage related to each of the zones as described above.

4. Methodology

4.1. Experimental Description

For the laboratory core experiments, we used the direct pulse transmission technique to record compressional and shear-wave velocities at ultrasonic frequencies (1 MHz). The experimental setup included a pressure vessel to encase and apply an isostatic confining pressure to the core sample; heat tapes wrapped around the pressure vessel to control the temperature; independent pumps for regulating confining pressure and pore pressure; a pulse generator; a digital oscilloscope for recording wave forms. A transducer was placed at each flat end of the core plug sample and was capable of recording a compressional wave and two shear waves in orthogonal directions. The two shear waves were oriented along the fast- and slow-velocity directions relative to the core sample at bench-top conditions. Each transducer casing had a fluid line for injecting the different fluids through the sample. In addition to recording the different wave modes, average axial strain was measured along the major axis of the sample. Uncertainties from the experimental setup correspond to a velocity error of 0.3% to 0.4%. However, the actual error is higher due to hand-picking first arrivals.

4.2. Testing Sequence

We started the testing sequence with dry rock measurements where we incrementally increased the confining pressure from 3.45 MPa to 41.37 MPa at the estimated reservoir temperature of 63.9°C.

Following the dry rock measurements, we vacuum saturated the core sample with brine then recorded measurements at different pore pressures. The fluids flushed through the samples following brine were live oil, live oil with a 0.334 mol fraction of CO2, and pure CO2. The pore pressure range for the brine-saturated samples was 29.65 MPa to 3.45 MPa, to simulate the stress path from an injection well to a production well. For the other fluids, live oil, live oil with a 0.334 mol fraction of CO2, and pure CO2, the pore pressure range was 29.65 MPa to 6.89 MPa. For these fluids, we did not record data below 6.89 MPa since these pressures would be below the bubble point. For the fluid-saturated and fluid-flushed cases, a constant confining pressure of 44.82 MPa is maintained to stimulate the stress of the overburden, and all the measurements were recorded at 63.9°C.

4.3. Fluid Properties

The fluids used to saturate and flush through the core samples were reproduced to simulate the actual reservoir fluids at Postle Field. These fluids consist of brine with an NaCl concentration of 142,000 ppm, live oil with a gas oil ratio (GOR) of 103 L/L, live oil with a 0.334 mole fraction of CO2, and pure CO2. The live oil has an API gravity of 38.6 and a gas gravity of 0.786. This information was obtained from a reservoir fluid study provided by the current field operator. Using the Fluid Acoustics for Geophysics (FLAG) calculator [11], the live oil bubble point is estimated to be 13.33 MPa at 63.9°C. The mixture consisting of live oil and a mole fraction of 0.334 CO2 has an estimated bubble point of 17.14 MPa at 63.9°C. The GOR and volume of gas for the live oil with a 0.334 mole fraction of CO2 was calculated assuming a gas gravity of 1.53 and a molecular weight of 44.01 gm/mole for CO2. This approach was used by Wang et al. [35] based on experimental measurements of Postle Field oil and mixtures of approximately 10 mole% and 50 mole% CO2. Figures 4 and 5 show the properties of the fluids over the pore pressure range used in the experiment calculated with the FLAG calculator [11]. The phase of CO2 for these pressures is supercritical fluid with the exception of the lowest pressure, 6.89 MPa, where CO2 is in the gas phase.

5. Results and Discussion

First, we measured compressional wave (P-wave) and shear wave (S-wave) velocities in air-saturated samples (dry rock) while loading and unloading an isostatic confining pressure. Following the dry rock testing, we measured P- and S-wave velocities, as a function of pore pressure, on brine-saturated core samples which were then flooded with live oil, live oil with a 0.334 mol fraction of CO2, and pure CO2, with a constant confining pressure of 44.82 MPa.

5.1. Core-Measured Velocities

The measured velocities, P- and S-wave, are shown in Figures 610 for the dry rock and fluid-saturated core samples. The S-wave velocity shown in these figures is the average of the two orthogonally measured S-waves. Differential pressure, , is defined as confining pressure, , minus pore pressure, : (1);

In general, the high-permeability samples show a greater sensitivity to both confining pressure and pore pressure. Both P-wave and S-wave velocities decrease with decreasing differential pressure (i.e., increasing pore pressure). The P-wave velocities show a larger decrease in velocities due to pressure when flooded with CO2. Whereas pore fluid has little bearing on S-wave velocities. Tables 3 and 4 show the percentage change in velocity for dry rock and fluid-saturated core samples over the differential pressure range of 15.17 MPa to 37.92 MPa.

In terms of different lithological characteristics of the core samples with respect to velocity, the higher the permeability of the sample, the lower the P-wave and S-wave velocities. This relationship cannot be extended to porosity since the porosities of the lower-permeability samples were very to the high-permeability samples.

The P-wave velocities show a response to both pressure and fluid type. The P-wave velocities systematically decrease, where brine saturated is the fastest and CO2 flushed is usually the slowest. The S-wave velocities’ response is mainly dominated by pressure and generally the fluid effect is negligible.

Some specific observations with respect to each of the core samples are as follows.

5.1.1. Cemented Sample

The S-wave velocity for the brine-saturated sample has a slightly higher velocity than when CO2 flushed (Figure 6(b)). More commonly the opposite behavior is observed.

5.1.2. High Permeability Number 1

The dry rock P-wave velocities (Figure 7(a)) show an unusual pattern where the velocity due to loading is higher than the velocity while unloading. This loading and unloading cycle may have caused permanent damage to the core leading to this response. Additionally, the P-wave velocities of the mixture consisting of live oil with a mole fraction of 0.334 CO2 and pure CO2 become nearly equal at a differential pressures of 31.03 MPa and 37.92 MPa (equivalent to pore pressures of 13.79 MPa and 6.89 MPa). The estimated bubble point for the oil-CO2 mixture is 17.14 MPa, suggesting the dominant phase of the mixture at these pressures could be gas.

In regards to the S-wave velocities (Figure 7(b)), the dry rock velocities do not show the same pattern as the P-wave velocities during the loading and unloading cycle. The dry rock S-wave velocities exhibit a response consistent with the other samples, where the velocity slightly increases during the unloading cycle. This sample showed the most S-wave velocity response to fluid; where the brine-saturated velocity is slower than when flushed with the oil-CO2 mixture and CO2.

5.1.3. High Permeability Number 2

Much like the high-permeability number 1 sample, the high-permeability number 2 sample (Figure 8(a)) shows a very similar P-wave velocity response to live oil with a mole fraction of 0.334 CO2 and pure CO2. These velocities nearly converge at differential pressures of 31.03 MPa and 37.92 MPa. Also, the velocity of the sample when flushed with live oil begins to converge with the oil-CO2 mixture and CO2 at differential pressures of 31.03 MPa and above. Again, this is most likely due to the effect of gas coming out of solution (live oil has a bubble point estimated to be 13.33 MPa at 63.9°C).

5.1.4. Low Permeability Number 1

The most obvious difference between the low-permeability number 1 sample and the other core samples is the much higher dry rock velocity for both P- and S-waves when compared to the fluid-saturated and fluid-flushed velocities (Figure 9). These higher dry rock velocities are the result of being recorded at 20°C, while the fluid measurements were recorded at the estimated reservoir temperature of 63.9°C.

5.1.5. Low Permeability Number 2

The low-permeability number 2 sample S-wave velocities (Figure 10(b)) show a similar trend as the cemented sample where the brine-saturated velocity is at some pressures marginally higher than when the sample was flushed with live oil and CO2. Again, typically the opposite effect is observed.

5.2. Average Axial Strain

Average axial strains measured during the dry rock loading and unloading cycle for three samples, each from a different lithological zone, are shown in Figure 11. The isostatic confining pressure has a greater effect on the high-permeability sample and least affects the cemented sample. We used these results to correct for sample length, porosity, and density at each differential pressure.

5.3. Log Facies Model

One method to link the core sample information to well log response is through log facies modeling. We generated a log facies model using 12 wells logged with the same technology within Postle Field. These well logs were corrected for environmental effects related to pressure, temperature, and borehole diameter. The well logs used as inputs for the log facies model were gamma ray, neutron porosity, bulk density, and photoelectric logs. The reservoir zones from each well log were spliced together in series to capture the variability of the reservoir. To identify the different log responses within the reservoir zone, principal component analysis was applied to maximize variability followed by cluster analysis to group data by similar log response. Figure 12 shows the box and whisker plot for each of the four clusters and the log response associated with each cluster. For each box, the central green line is the median of the data, the edges of the box are the 25th and 75th percentiles, and vertical lines (whiskers) extend to the extreme data points not considered outliers, the outliers are plotted as red diamonds. Based on the trends of the four box and whisker plots for each log response, log facies number 1 is interpreted as high-quality reservoir, log facies number 2 as intermediate-quality reservoir, log facies number 3 as low-quality reservoir, and log facies number 4 represents an interbedded shale.

Once the log facies are established, a supervised classification method is then applied to wells with additional data such as sonic or core data. This is shown in Figure 13 where the log facies model has been applied to a well with core-measured porosity and permeability. From Figure 13, an inference can be made that the dark blue log facies (number 1) relates to high-permeability zones, the light blue log facies (number 2) correlates to intermediate-permeability zones, and the green log facies (number 3) corresponds to zones of low permeability.

Unfortunately, the well which the core samples used for the laboratory velocity measurements did not have the necessary input logs for supervised classification. To relate the core samples to the log facies model, we compared the permeability of the core samples to the permeability associated with each log facies taken from the well in Figure 13. Figure 14 shows this comparison where the core samples used in the laboratory velocity measurements are shown as stars and each box and whisker plot represents the variability in permeability for each of the reservoir log facies. This assumes each core sample used in the laboratory velocity measurements corresponds to a log facies. Hence, the cemented sample is associated with the low-quality reservoir log facies, the high-permeability number 1 and number 2 samples are related to the high-quality log facies and the lower-permeability number 1 and number 2 samples correspond to the intermediate-quality log facies.

5.3.1. Integrating the Log Facies Model and Core-Measured Velocity

We achieved only limited success using Gassmann’s relation [9] and the effective medium theory of Mori and Tanaka [36, 37] to predict velocity changes related to fluid saturation and differential pressure in the core samples. Generally, Gassmann’s relation is not expected to model fluid effects at ultrasonic frequencies. This result only means Gassmann’s relation is not appropriate for this laboratory data but may be suitable at the field scale [38]. Therefore, we apply an alternative approach which assumes a direct relationship between the log facies model and the core-measured velocities. This direct relationship is based on estimating the reservoir condition in the laboratory velocity measurements, then applying the changes associated with fluid type and pressure observed in the laboratory measurements directly to the only cross-dipole velocity well logs recorded in the Postle Field study area.

The reservoir condition from petrophysical analysis estimates a water saturation of 45% [34]. Figure 15 shows the properties of live oil, brine, and the reservoir condition with a water saturation of 45%. To estimate the reservoir condition the difference in velocity between oil and brine is taken and is 76% less than the velocity of brine.

This estimated velocity at reservoir conditions now allows for changes to be applied to the log response based on laboratory velocity measurements made on the core samples. For example, to estimate the log response at complete brine saturation, the log would be multiplied by the percent change between brine-saturated velocity and the estimated reservoir velocity from the experimental core measurement.

Since very little fluid effect is observed in the S-wave velocity core measurements, it is assumed that the fluid effect is negligible and the only changes in S-wave velocity are due to pressure. For estimating changes in the S-wave log response, the S-wave velocity is multiplied by the percent change from the reservoir condition to the desired pore pressure.

The core samples we selected to estimate the changes in P- and S-wave velocity log response were the cemented sample which is associated to the low-quality reservoir log facies, the high-permeability sample number 2 which is associated to the high-quality reservoir log facies, and the lower-permeability sample number 2 which is associated to the intermediate-quality reservoir log facies. Changes in bulk density were calculated using the density transformation equation (2) where is bulk density saturated with the new fluid, is the initial bulk density, is the density of the initial fluid, and is the density of the new fluid:

Figure 16 shows the changes in log response for the conditions of complete brine saturation at pore pressures of 6.89 MPa and 29.65 MPa. Figure 17 shows the changes in log response for the conditions of complete CO2 saturation for pore pressures of 6.89 MPa and 29.65 MPa.

5.4. Time-Lapse Amplitude Difference

We generated synthetic seismic gathers using the full waveform modeling method described in Singh and Davis [39]. To understand the P- and S-wave response at the Morrow A sandstone level, we created 1D models for the brine- and CO2-saturated cases using the modified well logs (Figures 16 and 17). The logs were blocked at 1.5 m and simplified to have three main reflectors: the Atoka limestone above the reservoir, the Morrow A sandstone (reservoir), and a limestone layer below the reservoir. Above the Atoka limestone, the log response was smoothed to remove the effect of multiples.

From the 1D models, we generated 2D multicomponent gathers simulating vertical and horizontal sources, and receivers with vertical and horizontal components. We extracted statistical wavelets from the processed P- and S-wave seismic volumes recorded at Postle Field for use in the full waveform modeling. The P-wave wavelet was zero phase with a peak frequency of 18 Hz and the S-wave wavelet was also zero phase with a peak frequency of 13 Hz. The gathers (Figure 18) are modeled with a maximum offset of 1829 m, which is consistent with the usable offset from the seismic data recorded at Postle Field.

Seismic acquisition parameters and examples of P- and S-wave seismic data recorded at Postle Field can be found in Singh and Davis [39]. Using nearly identical wavelets and the same cross-dipole sonic logs used in this study, Singh and Davis [39] also address the tuning and interference effect at the Morrow A sandstone level.

To evaluate the changes in seismic response, we analyzed the amplitude difference of the P- and S-wave gathers for the brine- and CO2-saturated cases. The amplitude difference for the brine-saturated case represents a pore pressure increase from 6.89 MPa to 29.65 MPa in the Morrow A sandstone. Similarly, the amplitude difference for the CO2-saturated case represents a pore pressure increase from 6.89 MPa to 29.65 MPa in the Morrow A sandstone. The time shift for the brine and CO2 cases was minor, calculated to be less than 0.5 milliseconds. The P-wave amplitude differences are shown in Figure 19 and the S-wave amplitude differences are shown in Figure 20.

For both the brine and CO2 cases, the P-wave modeling shows a negative amplitude difference for an increase in pressure (Figure 19), this occurs because an increase in pore pressure results in a lower acoustic impedance. For both cases, the amplitude difference is greatest at zero offset and becomes less with increasing offset. The brine-saturated case shows a large amplitude difference starting from the reflection at the Morrow A sandstone level. The CO2 case shows a very weak amplitude decrease at the Morrow A sandstone reflection and a larger amplitude difference from the reflectors below the Morrow A sandstone. A change in waveform around the Morrow A sandstone for the CO2 case affects the signature of the amplitude difference. These observations suggest near-angle P-wave stacks would be better suited for time-lapse analysis with larger amplitude differences occurring below the reservoir.

Increasing pore pressure lowers the elastic impedance and hence the S-wave amplitude difference is negative with increasing pore pressure for both the brine and CO2 cases (Figure 20). The most notable observation is the amplitude difference with increasing offset between the two cases. The amplitude difference is highest at zero offset for the brine-saturated case and decreases with increasing offset. However, the CO2-saturated case shows the amplitude difference to increase with increasing offset. Analysis on S-wave angle stacks, rather than full offset stacks, may be more effective in detecting changes in saturation and, therefore, useful for CO2 monitoring.

The amplitude difference of the full waveform seismic modeling shows the effects fluid saturation and pressure have on P- and S-wave reflections. P-wave reflections show the near offsets to have the most amplitude difference. While S-wave reflections show a maximum amplitude difference at near offsets for the brine case and at far offsets for the CO2 case. Also, the S-wave amplitude difference is greater than the P-wave amplitude difference for these two cases.

To determine if the time-lapse changes are mainly due to time shift or only an amplitude difference, we compared the root mean square (RMS) amplitude difference for the two cases. Since we did not take into account any time shift when we differenced the modeled synthetic seismic gathers, representing the minimum and maximum pressures for the brine and CO2 cases (Figures 19 and 20), we calculated the RMS amplitude difference at the Morrow A sandstone top. We examined a 20 ms window centered at the Morrow A sandstone top for the P-wave gathers and a 40 ms window for the S-wave gathers. These RMS amplitude differences are shown in Figure 21 for both the brine and CO2 cases. We observe that the RMS amplitude difference at the Morrow A top is comparable to the amplitude difference of the synthetic gathers shown in Figures 19 and 20. The RMS amplitude difference with offset also agrees with the amplitude difference of the synthetic gathers. A comparison of the RMS amplitude difference and the differenced synthetic seismic gathers demonstrate that the effect of time shift is almost negligible at the Morrow A top. The time-lapse changes at the Morrow A top (observed in Figures 19 and 20) are mainly due to amplitude differences. Figure 21 also confirms that S-wave gathers show a larger time-lapse change as compared to the P-wave gathers.

6. Conclusions

Laboratory core measurements show the effects of differential pressure and fluid saturation through the P- and S-wave velocity response of five Morrow A sandstone samples. The samples were saturated with brine and flushed with live oil, an oil-CO2 mixture, and pure CO2. The P-wave velocity shows a sensitivity to both fluid saturation and pressure. Whilst the S-wave velocity shows a dependence to changes in differential and a nearly negligible sensitivity to fluid saturation. Core samples taken from the high-permeability zone show the most velocity change to confining pressure and pore pressure. The core sample from the cemented zone shows the least amount of velocity change to confining pressure and pore pressure. The overall velocity decrease was greater for S-waves as shown in Tables 3 and 4. Average axial strain also showed the high-permeability sample to be the most stress-sensitive and the cemented sample the least sensitive to confining pressure.

We attempted modeling the core-measured velocities using Gassmann’s relation [9] and the effective medium theory of Mori and Tanaka [36, 37]. Neither of these models successfully predicted the measured P-wave or S-wave velocities. The P- and S-wave velocity well logs were modified to estimate changes due to pressure and fluid saturation based on the laboratory core-measured velocities. This approximation was implemented through a log facies model where each log facies was associated with a core sample.

The P-wave amplitude difference for the brine- and CO2-saturated cases, modeling a pore pressure increase from 6.89 MPa to 29.65 MPa, displayed a maximum amplitude difference at the near offsets. Whereas the S-wave shows a maximum amplitude difference at the near offsets for the brine case and a maximum amplitude difference at far offsets for the CO2 case.

The full waveform modeling results show how fluid and pressure changes in the Morrow A sandstone reservoir affect the characteristics of P- and S-wave amplitude differences. In addition, the full waveform modeling shows an advantage to using angle-limited stacks as compared to full offset stacks for time-lapse analysis. These results demonstrate the benefit and establish motivation for using S-wave seismic data in time-lapse studies involving changes in pressure and fluid saturation.


The authors would like to thank the sponsors of The Reservoir Characterization Project at Colorado School of Mines, specifically Whiting Petroleum: the current operator of Postle Field; the Center for Rock Abuse at Colorado School of Mines for the use of their Rock Squeezing Laboratory, especially Professor Michael Batzle and Weiping Wang; Tom R. Bratton for his help and guidance generating the log facies model; BP America for the use of the seismic modeling code. A. V. Wondler would like to thank BP for the BP Colorado School of Mines Ph.D. Fellowship.