CHAPTER 7

Interpretation

7.1 Introduction

Interpretation of the multicomponent data volumes addresses different aspects of the RCP Phase VI project. The most general aspect of the interpretation is to document the basic reflectivity, signal quality and polarization of the compressional and shear-wave data. Another level of interpretation is integrating the seismic, geological and reservoir engineering data to form a static, multidisciplinary model of the reservoir. The third aspect is manipulating the differences between the initial and repeat surveys to reveal characteristics of the dynamic response of the bulk properties to the changing reservoir conditions. Integration with reservoir engineering data, fluid composition, volumes and pressures forms a dynamic model of the permeability structure of the reservoir. I will describe my efforts toward these aspects, particularly the third aspect, dynamic response.

7.2 Static Data Interpretation

Examples of the north - south vertical migrated section Inline 69 are shown in Figures 6.38 through 6.42 for the P-wave, S-wave S1 and S-wave S2 data. Inline 69 traverses both the CO2 injection well, CVU 97, and well CVU 200 which is the location of the downhole 3-C geophone and 9-C VSP survey. CVU 97 is at crossline 67, while CVU 200 is located at crossline 56. Figure 6.38 are examples of the p-wave data for the initial and repeat surveys. The reflected events include both the shallow reflections, the San Andres reservoir interval (0.68 seconds at CVU 200) and the deeper Paleozoic horizons below 1.3 seconds. Both data are filtered 10-72 Hz. and displayed with identical gain applied. Data within the migration halo, approximately 24 CDP bins or 1320 feet, have been excluded from the data volumes used for interpretation.

Figure 6.39 and 6.40 show the four shear data volumes processed through migration for the initial or pre-CO2 survey. These volumes represent the S1 (118o) and S2 (28o) polarizations and the two off-diagonals, D1 and D2. Data are filtered 10-24 Hz. and displayed with the same gain. The delay of approximately 10 ms due to shear-wave birefringence can be seen and most of the reflected energy is contained on the S1 and S2 sections. There is more energy evident on the off diagonal sections in the deeper section below 2.0 seconds which may be due to lower S/N ratio but may also be an indication of variations in local anisotropy in the deeper Paleozoic section.

Figures 6.41 and 6.42 are the four shear data volumes processed through migration for the repeat or post-CO2 survey. These four volumes, S1, S2, D1 and D2 can be compared to those from the initial survey. Again, the energy in the off-diagonals is minimal indicating the Alford rotation of the four horizontal components, within the offset window used in processing, is appropriate. Time to depth relationships for P-wave and S-wave data are from analysis of the VSP data described in Voorhies, 1996. The P-wave and S-wave data examples are displayed at a 2:1 time scale to facilitate comparison of the reflectivity between compressional and shear data.

One of the objectives in RCP Phase VI is to demonstrate the applications of repeated acquisition and processing of 3-C, 3-D seismic data. The primary application is to show the ability of time-lapse 3-C, 3-D to detect and monitor rock and fluid property changes associated with a production process. A requirement of the data is that it is repeatable. Examination of Figures 6.38 through 6.42 provides a visual, qualitative examination of the repeatability of the data acquired and processed in RCP Phase VI.

Vertical seismic sections from the migrated data volumes are displayed in Figures 7.1, 7.2, 7.4, 7.5, 7.7 and 7.8 at an expanded time scale. The vertical time sections are north -south profiles through wells CVU 97 and CVU 200. Figure 7.1 is a P-wave profile from the initial or pre-CO2 injection survey. Significant horizons are annotated. Reflection time of the top of San Andres is 0.68 seconds at CVU 200 based on the VSP data. Five hundred feet of reservoir, encompassing 160 feet of Upper San Andres and 340 feet of Lower San Andres, represent approximately 52 ms of two way travel time. The same north - south traverse from the repeat or post-CO2 injection P-wave survey is included for comparison (Figure 7.2). A bandpass filter, 10-72 Hz., is applied to both data volumes.

FIGURE 1. P-wave migrated data example for north-south profile Inline 68 from initial survey. Time scale is in ms.
FIGURE 2. P-wave migrated data example for north-south profile Inline 68 from repeat survey. Time scale is in ms.

The P-wave data are visually very similar between the initial and repeat surveys. Lithologic boundaries with higher reflectivity, Rustler, Yates and Glorietta, are essentially identical. One measure of repeatability is to form a simple difference stack by subtracting one volume from another. Figure 7.3 is a difference ensemble computed by subtracting the repeat survey data from the initial survey. This two-panel display shows data from the repeat survey in the upper panel with the difference section in the lower panel. The same display gain is applied to both panels.

FIGURE 3. P-wave difference ensemble. Upper panel are data from the repeat survey; lower panel are the difference traces. Both data sets are plotted at the same gain.

Figures 7.4 and 7.5 show the S1 and S2 polarizations for the S-wave data from the initial survey. The Rustler, Yates, Queen and San Andres are annotated. Comparison with the initial survey P-wave data (Figure 7.1), yields observations about bandwidth, reflectivity and reflection continuity.

FIGURE 4. S-wave (S1) migrated data example for north-south profile Inline 68 from the initial survey. Time scale is in ms.
FIGURE 5. S-wave (S2) migrated data example for north-south profile Inline 68 from the initial survey. Time scale is in ms.

The bandwidth of the P-wave migrated data is approximately 10-72 hz. while the S-wave bandwidth is much lower, approximately 10-24 Hz. Bandwidth, in this context, is the portion of the signal spectrum that is repeatable between the initial and repeat surveys using the processing flow outlined in Section 6. Please note that more aggressive processing techniques could be applied to the data to increase the apparent frequency content, but results in a greater amount of noise in the repeatability measurements such as difference stacks.

The overall reflectivity can be correlated between the P-wave and S-wave data but differences can be seen, particularly in the shallow section above the Rustler. In the shallow well control there is an interface with an abrupt increase in density. This coincides with the strong reflectivity at 0.80 seconds on the S-wave data. This reflection is the dominate energy on the S-wave data but not on the P-wave data. For both the P-wave and S-wave volumes, the reflection strength within the reservoir zone is weak due to low reflection impedance. The low reflectivity was verified using dipole sonic log information from two wells in the survey area (Roche, RCP Sponsors Meeting Notes, September 25, 1995; DeVault, RCP Sponsors Meeting Notes, April 15, 1996). Reflection continuity is good on both the P-wave and S-wave data for the section above the San Andres horizon (Queen to Rustler). Continuity of the S1 and S2 data is comparable. The interface of high reflectivity at 0.800 seconds on the S-wave data corresponds to the base of an interval of high Vp/Vs ratio identified in the VSP data (RCP Sponsors Meeting Notes, October 14, 1996).

Over the reservoir interval on the S-wave data, approximately 1600 to 1700 ms, the reflection continuity is less continuous. I ascribe the discontinuous nature of the reflected energy to be caused by fracturing and local stress variations in the more rigid dolomite as opposed to the less rigid lithology above the San Andres. The anhydrite and evaporite content of the section above the San Andres increases (see Figure 2.5) as the basin desiccates. The Salado is primarily salt and anhydrite. Lower S/N ratio could also be a factor, with the high S-wave reflectivity at 0.80 seconds acting as an elastic wave barrier to energy transmission.

Analysis of reflection continuity by Talley, 1997, associates faulting and fracturing with discontinuities in the reflection character. His work identifies faults and fracture systems propagating upward from basement faults such as the north - south trending regional fault system identified by Hill, 1984, and Shumaker, 1990 (Figures 2.2 and 2.3). Fracture systems in the San Andres are influenced by the deeper faulting and modified by the east - west trending Permian shelf edge margin (Figures 2.8 and 2.9). Talley, 1997, notes the intensity of the reflection discontinuity, related to faults and fractures, diminishes above the San Andres. The increased salt and anhydrite content may reduce the rigidity of the shallower lithology, lessening the fracture intensity and spatial variations in the local stress field.

Within the reservoir interval, the S2 reflection character is more continuous than that of the S1 data. Core analysis (Scuta, in RCP Sponsors Meeting Notes, April 15, 1996) show fractures are present, both open and healed with anhydrite. Fractures are more numerous along the shelf edge, CVU 100, as opposed to north of the shelf edge, CVU 345. Analysis of fracture orientation from a single FMI log in the Warn State 26-1 well, approximately 2000 feet south - southwest of CVU 97, measures high-angle fracture sets, parallel and perpendicular to the southwest to northeast trending shelf margin (Scuta, in RCP Sponsors Meeting Notes, April 10, 1997). The more discontinuous nature of the S1 reflections in the reservoir zone may be related to more numerous fractures trending along the shelf edge. Fault and fracture trends delineated by seismic continuity analysis also parallel the San Andres shelf edge (Talley, 1997).

Evidence for southwest to northeast trends can be seen in the P-wave data and reservoir production statistics. Interpretation of the P-wave data provided a structural image of reservoir geometry. Ties to well control are excellent. Time to depth conversion was accomplished using the time : depth relationships (P-wave and S-wave) obtained from the VSP survey at CVU 200. Data were datumed at the Queen horizon, then a single spatially invariant velocity function was applied to convert the deeper horizons to depth. Ties to well control were in good agreement through the Grayburg to the top of the upper San Andres. The seismic response beneath the top of San Andres is primarily a "tuning" phenomena, thus a horizon pick on a waveform associated with porosity changes within the reservoir did not produce an accurate depth tie to a specific marker picked from well control. Figure 7.6 is a perspective view of the time structure of the top of San Andres from the P-wave data volume. Perspective view is from the southwest toward the northeast. The shelf crest is visible, trending away from the viewer. The blacked-out contour highlights a southwest to northeast lineation through the center of the survey area. This feature coincides with reflection character discontinuities (Talley, 1997) and with the general demarcation of the reservoir production data dividing the reservoir into two general regions. The reservoir produces more fluid at a higher water cut with less gas to the north and northwest of this lineation (Figures 3.2 to 3.7).

FIGURE 6. Perspective view of P-wave San Andres time structure, southwest to northeast view.

These data provide several fundamental insights concerning time-lapse, multicomponent seismology. The first insight is not related to time-lapse issues, rather the different images of the subsurface reflectivity provided by the S1 and S2 volumes. The S1 and S2 reflectivity are providing different and complementary information of the subsurface bulk properties. Note that reflectivity is more similar in the shallow section, Rustler to Yates interval, which is dominated by thick salt and anhydrite formations. Below the Yates marker the section is predominately carbonates. With respect to stress, salt is a more isotropic material that may absorb local stress variations while the more rigid carbonate rock would preserve local stress variations. It is possible that the variations in S1 and S2 reflectivity are delineating variations in lithology, porosity, stress conditions and fluid characteristics (pressure and properties). An alternate explanation is a decrease in signal / noise ratio beneath the high shallow reflectivity. Note that examination of the deeper S-wave data (figures 6.39, 6.40, 6.41 and 6.42) shows good S/N but poorer repeatability at the deeper Paleozoic reflections.

The compressional and shear wave reflectivity are similar but different in the Phase VI data volume. This reaffirms the observation that p-wave and shear-wave data provide more information about the subsurface reflectivity when used together. Compressional data provide a measure of the bulk rock compressibility, rigidity and density while shear-wave data is sensitive to rigidity and density. Compressibility and rigidity are bulk properties and are likely to be anisotropic. In porous media, pore fluid pressure and fluid properties affect the bulk properties. The combined use of p-wave, S1 shear and S2 shear seismic data allow different views of the bulk properties in the subsurface. Interpretation strives to use these volumes to delineate spatial variations in the subsurface related to lithology, porosity, variations in pore structure related to preferred permeability directions, and variations in pore fluid pressure and properties.

Figures 7.7 and 7.8 show the S1 and S2 polarizations for the S-wave migrated data for Inline 68 from the repeat survey. Reflectivity above the reservoir zone is very similar between the initial and repeat surveys for both the S1 and S2 data. Subtle waveform variations between the Rustler and Yates interval compare very closely ( for example 0.93 seconds or 1.08 seconds). Difference ensembles for the S1 and S2 data for the Inline 68 traverse are shown in Figures 7.9 and 7.10. The high amplitude reflectivity at 0.80 is virtually absent on the difference stacks indicating the acquisition and processing have produced the same image, in amplitude and phase, of the shallow reflectivity that is assumed to be constant between the initial and repeat surveys. There are differences between the surveys in the reservoir zone and even above the reservoir zone that will be discussed in the dynamic interpretation section.

FIGURE 7. S-wave (S1) migrated data example for north-south profile Inline 68 from the repeat survey. Time scale is in ms.
FIGURE 8. S-wave (S2) migrated data example for north-south profile Inline 68 from the repeat survey. Time scale is in ms.
FIGURE 9. S-wave (S1) difference ensemble. Upper panel are data from the repeat survey; lower panel are the difference traces. Both datasets are plotted at the same gain.
FIGURE 10. S-wave (S2) difference ensemble. Upper panel are data from the repeat survey; lower panel are the difference traces. Both datasets are plotted at the same gain.

The north - south traverse, inline 68, extracted from the migrated data volume through wells CVU 97 and 200 was selected since significant reservoir processes occurred at these locations. There are other traverses which could be shown where the reflectivity between the initial and repeat surveys is virtually identical and the difference stacks show very little energy. This forms a second significant insight, the repeatability of the S-wave data volumes, especially above the San Andres formation. Although the S1 and S2 volumes exhibit less continuity than the P-wave volumes, shear-wave data provide very similar images for the initial and repeat surveys. I interpret this observation to indicate shear-wave data provide real and valid information about subsurface bulk properties. The lack of continuity of S-wave data when compared to P-wave data should not be construed as "poorer" data quality. It is likely to be actual information regarding the spatial distribution of rigidity, density and velocity anisotropy in the subsurface.

7.3 Vp/Vs Ratio

A measure of Vp/Vs ratio was derived from the multicomponent data volume. To improve the S/N ratio of the P, S1 and S2 data, the initial and repeat survey volumes were summed. Using the depth control from the VSP, horizons were picked on the P, S1 and S2 migrated and summed volumes that correspond to the Queen and a reflection in the lower San Andres called the "cycle 1" marker. This depth interval is approximately 1000 feet thick and spans approximately 400 feet of upper and lower San Andres reservoir. The resulting Vp/Vs1 and Vp/Vs2 ratios are shown in Figures 7.11 and 712. Due to the probability that the "cycle 1" horizon pick may be influenced by tuning, and is not a precise interface pick, I present these Ts/Tp measures as estimates of Vp/Vs. The values shown have been scaled by 0.82 to calibrate to those values measured from the dipole sonic at CVU 200. The dipole sonic shows the San Andres reservoir zone, 4300 feet to 4800 feet, has a Vp/Vs ratio of 1.8 to 1.9. The application of a single scalar correction, 0.82, to all Vp/Vs ratio measurements obtained from the surface data indicates the Vp/Vs ratio to be approximately 1.85 at the CVU 345 location. CVU 345 is at the northern limits of the migrated data volume and is one of the three wells with core available. A Vp/Vs ratio of 1.85 is in excellent agreement with the Vp/Vs measurements obtained from core samples from CVU 345 in the San Andres reservoir zone (Capello de Passalacqua, 1995).

FIGURE 11. Map of Vp/Vs1 ratio from summed 4-D, 3-C data volumes. The S1 data was rotated to a constant 118 and does not reflect local variation in the fast shear wave polarization.
FIGURE 12. Map of Vp/Vs2 ratio from summed 4-D, 3-C data volumes. The S2 data was rotated to a constant 28 and does not reflect local variation in the slow shear wave polarization.

There is a clarification required in the use of the terminology "S1" and "S2" to describe S-wave data volumes. In this thesis, S1 And S2 data are the results of 4-C, Alford, rotation to constant azimuthal directions 118o and 28o. These represent regional S1 and S2 polarization directions based on literature, VSP analysis and borehole breakout study. However, local changes in the S1, or fast, S-wave polarization may exist due to local variations in stress direction. These local variations in S1 polarization are imbedded in the regional "S1" and "S2" volumes. A possible example can be seen in Figures 7.11 and 7.12. In the region of wells CVU 96, 244 and 144 (east-northeast of CVU 97) there are high Vp/Vs1 anomalies oriented in the northeast direction. Within the context of the regional Vp/Vs1 map, this anomalous area may represent a change in the local S1 polarization direction to the northeast. The high Vp/Vs1 ratio is due to regional "S1" slowing in response to fractures and pore structure open in the northwest direction within this local northeast trend. In a sense, this local anomaly on the Vp/Vs1 map is a slow S-wave (or S2) feature.

The Vp/Vs1 and Vp/Vs2 measures show regions of lower Vp/Vs ratio in the southwest portion of the image area. This portion of the reservoir produces less fluid with a lower watercut than do areas of the field exhibiting a higher Vp/Vs ratio. The Vp/Vs measures in Figures 7.11 and 7.12 are compared to reservoir production characteristics shown in Figures 7.13 through 7.20. These production data represent the daily average over the three-month period before the CO2 injection project (see Figures 3.2 through 3.7).

Figures 7.13 and 7.14 are Vp/Vs ratios with contours of the average barrels of oil produced per day (BOPD). Figures 7.15 and 7.16 are Vp/Vs ratios with contours of the average gas - oil ratio (GOR). In the regions of lower Vp/Vs ratio, these figures indicate portions of the field, near wells CVU 204 and 206, produce more oil per day with a higher gas - oil ratio (GOR). A higher saturation of oil with gas would result in a less dense and more compressible fluid, reducing the P-wave velocity relative to the S-wave velocity (which would increase) thus lowering the Vp/Vs ratio. Analysis of Vp, Vs and Vp/Vs ratios from core samples for brine-saturated, oil-saturated and dry conditions show the same trend towards lower Vp/Vs ratios (Capello de Passalacqua, 1995). Even subtle variations such as the increase in Vp/Vs across well CVU 203 then decrease at well location CVU 201 correlate to the higher BOPD and GOR features in the production data.

FIGURE 13. Map of Vp/Vs1 ratio from summed 4-D, 3-C data volumes (Figure 7.11) with BOPD contour overlay (Figure 3.3).
FIGURE 14. Map of Vp/Vs2 ratio from summed 4-D, 3-C data volumes (Figure 7.12) with BOPD contour overlay (Figure 3.3).
FIGURE 15. Map of Vp/Vs1 ratio from summed 4-D, 3-C data volumes (Figure 7.11) with GOR contour overlay (Figure 3.7).
FIGURE 16. Map of Vp/Vs2 ratio from summed 4-D, 3-C data volumes (Figure 7.12) with GOR contour overlay (Figure 3.7).

Note that contours connect high GOR on the east-west row of production wells with the next row of production wells to the north (through CVU 87 and 186). Specifically, CVU 196 is grouped with CVU 186. These contours cross the east-west row of water-injection wells (CVU 93 through CVU 244). Since the GOR contouring is controlled by production wells only, I propose the Vp/Vs information from the seismic data is indicating higher water saturation in the interwell region between CVU 196 and CVU 186.

The production data shows that well CVU 97 is similar to the wells located to the northwest, that is, along the direction of maximum horizontal stress. Aspects such as watercut contours overlayed on Vp/Vs (Figures 7.17 and 7.18) show a northwest to southeast trend from wells CVU 88 to 93 and then to 97. The Vp/Vs measures correspond to these watercut trends, also suggesting a different spatial distribution than that of the production well contours.

FIGURE 17. Map of Vp/Vs1 ratio from summed 4-D, 3-C data volumes (Figure 7.11) with percent watercut contour overlay (Figure 3.4).
FIGURE 18. Map of Vp/Vs2 ratio from summed 4-D, 3-C data volumes (Figure 7.12) with percent watercut contour overlay (Figure 3.4).

Another indication of a general permeability trend of northwest to southeast through well CVU 97 can be inferred from the water injection volumes and surface tubing pressures. The area of higher Vp/Vs ratios greater than 1.9 in figures 7.11 and 7.12 delineates a "fairway" approximately 2500 feet wide in the southwest to northeast dimension, trending west - northwest to east - southeast through the center of the survey area. Figures 7.19 and 7.20 are Vp/Vs ratios with contours of surface tubing pressure overlayed. This trend of Vp/Vs higher than 1.9 generally coincides with the S1 direction. The preferential permeability direction will be along the S1 direction due to open fractures systems, microfractures and low-aspect-ratio pore structure. Although there is an operational overprint to the water injection pressures and volumes, this "fairway" allows the injection of a greater volume of water at a lower surface tubing pressure than outside the fairway. Note the anomalous region of a possible change in the local S1 direction near CVU 244 is in agreement with the surface tubing pressure data.

FIGURE 19. Map of Vp/Vs1 ratio from summed 4-D, 3-C data volumes (Figure 7.11) with surface tubing pressure contour overlay (Figure 3.6).
FIGURE 20. Map of Vp/Vs2 ratio from summed 4-D, 3-C data volumes (Figure 7.12) with surface tubing pressure contour overlay (Figure 3.6).

My interpretation of the higher Vp/Vs ratios in the fairway is indicative of more fractures and crack-like or low-aspect-ratio pore structure with better permeability; i.e., the reservoir is more compliant with less filling of porosity by anhydrite. S-wave velocity is reduced relative to P-wave velocity, increasing the Vp/Vs ratio. Arestad, 1995, presented similar results from the Joffre Field study in RCP Phase V. At Joffre, the reservoir is also a dolomite with the reservoir production characteristics governed by the anhydrite plugging of porosity and fractures. The RCP Phase V report included the observation that better reservoir production coincided with areas of higher Vp/Vs ratio.

There are different spatial patterns to the Vp/Vs1 and Vp/Vs2 ratios shown in Figures 7.11 through 7.20 that can be explained in terms of local changes in anisotropy. The data volumes were rotated to a single S1 and S2 azimuth; local variations in anisotropy will cause the fast shear-wave polarization to change. Analysis of anisotropy by Talley, 1997, showed evidence for spatial variation of shear-wave splitting in the San Andres reservoir interval and that areas of high anisotropy measures coincided with the better producing wells. The same spatial variation of S1 and S2 velocity associated with shear-wave birefringence is incorporated in the Vp/Vs1 and Vp/Vs2 measurements. For example, the narrow region of high Vp/Vs1 between wells CVU 196 and CVU 244 may indicate an area with open fractures or pore structure aligned in the northeast direction. A second example, the area near CVU 97 shows a higher Vp/Vs2 ratio, indicating open fractures or pore structure aligned to the northwest. Areas with higher Vp/Vs ratios for both the S1 and S2 data may suggest multiple azimuths of open fractures.

Note that the area of higher Vp/Vs2 ratio appears to cover more area within the fairway as opposed to the area of higher Vp/Vs1 ratio. The regional horizontal stress or S1 direction is dominate, opening more pores and fractures aligned in this direction thus increasing the area of higher Vp/Vs2 ratios. Areas with higher Vp/Vs1 ratios are smaller, more localized features, indicating changes in open fracture direction, micropore structure or stress conditions. Remember that these anomalous areas on "S1" maps (regional S1) may represent local S2 regions which are affected by fracture and pore structure.

During the infill drilling down to 10 acre spacing, two wells produced approximately four times the expected rate (Wehner, 1997). These wells are CVU 195, located 3 locations east of CVU 97 at the eastern edge of the migrated data volume, and CVU 203, located south - southeast of CVU 97. Both of these wells are located near regions of higher Vp/Vs ratios.

7.4 Repeatability

The design and acquisition of the 3-C, 3-D surveys provided nine possible data volumes to be processed and interpreted. With repeat survey or time-lapse element, there are 18 data volumes available for analysis. In this section I present the results concerning five of the nine possible volumes, the conventional p-wave and four shear-wave volumes. The four shear-wave datasets are the products of applying an Alford rotation to produce a survey invariant S1 volume, or fast shear-wave polarization of 118o, an S2 volume, slow shear-wave polarization of 28o, and the two off-diagonal datasets.

It is important to provide quantitative measures of the repeatability of the seismic data when interpreting time-lapse data volumes. One measure is to compute the RMS amplitude level of the difference section and compare this RMS amplitude difference value to the RMS amplitude level of the input sections. This analysis was done on the migrated data volumes, excluding data from the migration aperture along the edges of the survey. For the shallow data, all three data volumes show differences with a mean lower than -20 db, meaning the data are repeatable to less than one part in ten. This measure of computing the ratio of the RMS amplitude of the difference to the RMS amplitude of the input comprehends both amplitude and phase variations between the initial and repeat surveys. Figures 7.21, 7.22 and 7.23 are map views of the repeatability measure over the survey area for the p-wave, S1 and S2 data volumes. The analysis was derived from the reflectivity in the shallow portion of the section above the San Andres reservoir zone. This level of repeatability was attained using 3-C geophones planted at the earth's surface, vibratory sources at the surface and the processing techniques outlined in this report. I use this repeatability measure to state the background repeatability of the P-wave data to be approximately -20 to -24 db, as measured over the reflectivity in the shallow section above the San Andres reservoir zone.

FIGURE 21. P-wave difference measure between initial and repeat surveys over the shallow Yates interval. Measure is expressed in db as the ratio of the RMS amplitude difference trace to the RMS amplitude level of the input traces.
FIGURE 22. S-wave (S1) difference measure between initial and repeat surveys over the shallow Rustler interval. Measure is expressed in db as the ratio of the RMS amplitude difference trace to the RMS amplitude level of the input traces.
FIGURE 23. S-wave (S2) difference measure between initial and repeat surveys over the shallow Rustler interval. Measure is expressed in db as the ratio of the RMS amplitude difference trace to the RMS amplitude level of the input traces.

The measure of repeatability is also performed for analysis windows centered on the San Andres reservoir zone and deeper Paleozoic reflectivity (Figures 7.24 through 7.29). Repeatability diminishes with increased time. A probable explanation is a general decrease in the S/N ratio. Note that the P-wave repeatability is approximately -24 db at the shallow horizon (Figure 7.21) then decreases to -18 db for the San Andres interval (Figure 7.24) and decreases further to -12 db for the deeper Paleozoic reflectivity.

FIGURE 24. P-wave difference measure between initial and repeat surveys over the San Andres interval. Measure is expressed in db as the ratio of the RMS amplitude difference trace to the RMS amplitude level of the input traces.
FIGURE 25. S-wave (S1) difference measure between initial and repeat surveys over the San Andres interval. Measure is expressed in db as the ratio of the RMS amplitude difference trace to the RMS amplitude level of the input traces.
FIGURE 26. S-wave (S2) difference measure between initial and repeat surveys over the San Andres interval. Measure is expressed in db as the ratio of the RMS amplitude difference trace to the RMS amplitude level of the input traces.
FIGURE 27. P-wave difference measure between initial and repeat surveys over the deeper Paleozoic horizons. Measure is expressed in db as the ratio of the RMS amplitude difference trace to the RMS amplitude level of the input traces.
FIGURE 28. S-wave (S1) difference measure between initial and repeat surveys over the deeper Paleozoic horizons. Measure is expressed in db as the ratio of the RMS amplitude difference trace to the RMS amplitude level of the input traces.
FIGURE 29. S-wave (S2) difference measure between initial and repeat surveys over the deeper Paleozoic horizons. Measure is expressed in db as the ratio of the RMS amplitude difference trace to the RMS amplitude level of the input traces.

The decrease in repeatability is more rapid with time for the S-wave data. Reflectivity above and containing the predominately salt section is repeatable to approximately -24 to -30 db. However, beneath the salt and evaporite section the repeatability decreases abruptly to -6 to -12 db. It is probable that the higher reflectivity in the shallow section for the S-wave data allows less transmission of S-wave energy through the salt and evaporite section, thus decreasing the S/N ratio.

An alternate explanation is the more rigid carbonate strata beneath the Rustler and Salado are actually exhibiting a significant time-lapse signature associated with reservoir processes, tidal effects or other processes affecting the rigidity of the subsurface.

7.5 Dynamic Interpretation

The time-lapse, multi-component seismic data volumes have been analyzed for a temporal response related to the changes in bulk properties due to the changing reservoir conditions. There are observed variations in p-wave reflection amplitude and changes in shear-wave polarization. I believe these time-lapse seismic anomalies observed between the initial and repeat surveys are associated with changes in pore fluid properties and changes in effective stress due to pore pressure variations. The pore fluid property changes include density, compressibility and viscosity associated with the CO2 Huff-n-Puff project.

Figure 7.30 is a map view of the difference in p-wave RMS amplitude measured over a 40 ms window centered at the Yates reflection. This measure is normalized to be percent change relative to the reflection strength at the Yates horizon. Computation of the amplitude difference is similar to the repeatability measure described in the previous section. This measure indicates the initial and repeat surveys are repeatable to -20 Db, which corresponds to a +/- 10% reflection strength change. Figure 7.30 shows the repeatability at the Yates horizon to be approximately +/- 10%. The analysis at the Queen horizon (Figure 7.31), located above the San Andres reservoir zone, also indicates the majority of the data exhibits +/- 10% variation.

FIGURE 30. Time-lapse seismic attribute measurement: P-wave amplitude difference from a 40 ms window at the shallow Yates horizon.
FIGURE 31. Time-lapse seismic attribute measurement: P-wave amplitude difference from a 40 ms window at the Queen horizon.

Continue Chapter 7