Abstract

Morrow sandstones constitute important oil-producing reservoirs in the Anadarko Basin in the Mid-Continent Region of the USA. Characterization of the Morrow A sandstone reservoir in Postle Field, Texas County, Oklahoma, is challenging due to its small thickness, low acoustic impedance contrast with the surrounding Morrow shale, and lithological heterogeneity. Shear wave data have been documented as a promising solution for imaging the Morrow A sandstone. Vertical seismic profiling (VSP) offers the potential to enhance shear wave imaging of the thin heterogeneous Morrow A sandstone at Postle Field. The zero-offset VSP results confirm the advantages of shear wave over compressional wave in imaging the Morrow A sandstone. Also, the final shear wave VSP image shows that, applying the proposed processing flow, we are able to image the Morrow A sandstone where the thickness is as thin as 8.5 m.

1. Introduction

Pennsylvanian upper Morrow sandstones (Figure 1) constitute major oil-producing reservoirs across southeastern Colorado, southwestern Kansas, western Oklahoma, and northern Texas. They consist of multiple-stacked lenticular sandstone bodies formed within valley-fill complexes. The Morrow A is an oil-producing sandstone mainly in Oklahoma [13]. In this study, we mainly focus on the Morrow A sandstone in Postle Field, Texas County, Oklahoma (Figure 2).

Compressional wave studies have been mostly used in the past for characterizing the Morrow sandstones [47]. Briefly, there are two main challenges with characterizing the Morrow A sandstone based on compressional data: first, the thickness of the Morrow A sandstone is below the tuning thickness [6, 8] and secondly, the acoustic impedance contrast between the Morrow A and the surrounding Morrow shale layer is extremely low to the point that these sandstones are acoustically transparent [9, 10]. Dipole sonic well logs reveal the low acoustic impedance contrast between the Morrow A and overlying shale (Table 1 and Figure 3). Singh and Davis [10] showed that the low acoustic contrast between the Morrow A sandstone and surrounding Morrow shale layer and the presence of strong multiples from shallow anhydrite layers are important challenges associated with the detection of the Morrow A sandstone in compressional data (Figure 4). Full waveform modeling [10, 11] suggested that shear wave data can improve the Morrow A imaging and characterization significantly in Postle Field. Hardage et al. [12] drew the same conclusion after their study on thicker section of the Morrow sandstone using three zero-offset 9C VSP surveys. Also, a 3D VSP study using shear data addressed the advantages of shear wave in imaging of sandstone layers which have low acoustic impedance is contrast with surrounding layers [13].

Shear VSP data have characteristics which make this method a unique source of information. The advantages of VSP methods used in this study involve higher frequency, better wavefield separation, and higher signal-to-noise ratio. In this study, we expect to map the Morrow A sandstone and its variations around the VSP well. To achieve this purpose, certain processing steps and enhancements were applied that will be discussed.

2. Field VSP Data

As a part of Reservoir Characterization Project Phase XII at Colorado School of Mines, two 3D multicomponent VSP surveys were acquired from well HMU 13-2 in March and December 2008. Figures 5(a) and 5(b) demonstrate source locations and 3C geophone arrays in the well, respectively. For imaging purposes, we used only the March 2008 VSP survey in this paper.

The March 2008 survey was mainly acquired with horizontal vibrators, with the exception of two vertical shots: one 91 m to the north of the VSP well (called zero-offset data) and one 1829 m to the south of the VSP well (for the purpose of receiver reorientation). The horizontal vibrator truck generated two mutually orthogonal horizontal (E-W and N-S) ground motions. The geophones in the borehole were located from 933 m to 1818 m MD (their depth spacing is 15 m) from well head (Figure 5(b)). The last geophone was 60 m above the top of the Morrow A sandstone. Table 2 summarizes source and receiver parameters of the March 2008 VSP survey.

In addition to VSP data, well log and surface seismic data were used in this study. The locations of studied wells and surface seismic area are depicted in Figure 6.

3. VSP Data Processing

VSP processing was applied to three datasets from the March 2008 survey: zero-offset shear dataset, zero-offset compressional dataset, and multioffset (or 3D) shear dataset. The main goal of the proposed processing flow (Table 3) is to preserve the original amplitude and its variations due to lithology changes.

One of the early steps in shear wave data processing is to reorient the horizontal components of all geophones toward the constant geographical coordinate system, commonly called “Receiver Reorientation.” Hodogram analysis on direct compressional wave arrivals is one of the common approaches for finding the azimuth of each different geophone components in the borehole. Modeling showed that the hodogram analysis for a source with 488 m offset introduces up to 5 deg error in estimation of receiver azimuth. The error is expected to increase for more complex geological structures and noisy data. To reduce the uncertainty of the hodogram analysis, usually a number of shot gathers are analyzed. Since we had only one far-offset vertical source in our dataset, the error in the hodogram analysis and the receiver reorientation could have been considerable. Amplitude and time variations of the first shear wave arrivals in reoriented data (Figure 8(a)), using hodogram analysis of only the vertical source, confirmed the uncertainty in estimation of receivers azimuth. Uncertainty in receiver reorientation can diminish weak reflections from the reservoir or make invalid amplitude variations in the final VSP image. Because of this reason, we modified the receiver reorientation process by utilizing two more analyses: (1) hodogram analysis of compressional wave arrivals from a few horizontal vibrators and (2) 3C amplitude analysis. Hodogram analysis of the compressional wave generated by horizontal sources improved receiver reorientation results for some shallow geophones. Besides the mentioned methods, 3C amplitude analysis was done on the data which improved the results. In this approach, we rotated all three components of each geophone in a 3D half space and the RMS amplitude maps of the first arrivals of compressional and shear waves were projected onto two separate stereonets (Figure 7). The amplitude maps from 3C amplitude analysis are supposed to show the direction of maximum direct compressional and shear wave energy in 3D space around the receivers. Therefore, accurate receiver reorientation should lead to the similar amplitude pattern on all geophones. Using this approach, some residual reorientation angles were calculated. Interpretation of the projected amplitude maps, in addition to previous hodogram analyses, provided more accurate rotation angle for aligning receivers (Figure 8(b)).

After the receiver reorientation, based on the polarity of the shear wave direct arrivals, the data polarity for different shots was fixed. Different receivers were rotated to the local coordinate system (radial-transverse), and then they were tilted in the direction of direct shear wave to maximize shear wave energy on the horizontal components. Many of the compressional wave arrivals and reflections disappeared on the horizontal components after this step.

4. Shear Wave Splitting

When a shear wave enters an anisotropic medium, it splits into two polarized shear waves that travel with different velocities [14, 15]. Shear wave splitting is an important phenomenon that requires attention during shear wave data processing and interpretation. Typically, when the medium is azimuthally anisotropic, shear wave data will be processed and presented in two modes: fast shear wave data (S1) and slow shear wave data (S2). Finding the fast and slow shear wave polarization directions (or the natural coordinate system) and their spatial and temporal variations has a significant impact on the quality of the final shear wave image. Therefore, information from different tools gives us better insight about the reason for shear wave splitting, the polarization directions of fast and slow shear waves, and their temporal and spatial variations. Here, to show the importance of contributing different sources of information, we analyzed four different tools: full waveform sonic log (sonic scanner), image logs (FMI), multicomponent zero-offset VSP data, and multi-offset VSP data.

Slowness dispersion analysis (Figure 9(a)) on sonic scanner data from well HMU 24-4 (Figure 6) shows that the anisotropy in the study area, especially at the reservoir level, is mostly stress-induced anisotropy [16, 17]. In addition, the sonic scanner data suggest that the fast shear wave polarization direction is 126 deg (all the angles in this analysis are respect to north) at the reservoir level and 108 deg for the entire logged section (Figure 9(b)).

Since the origin of anisotropy in the field is presumed to be stress-induced anisotropy, fast and slow shear wave polarization directions should match maximum and minimum horizontal stress directions at the field. Drilling-induced fractures usually form in the direction of maximum horizontal stress [18], and they are detectable on micro-imaging logs (e.g., FMI logs). Figure 10 depicts the direction of some drilling-induced fractures in well HMU 24-4 and JULIAN 6-15 (Figure 6). Most of the drilling-induced fractures formed around 125 deg which is very close to the sonic scanner results. Although both image and sonic scanner logs are useful tools in determining the anisotropy direction, both datasets are limited to the deep part of the well. In addition, the radius of investigation of these sources of information is limited to the vicinity of the VSP well. To find the azimuthal anisotropy direction at shallower layers and also away from the borehole, we did a comprehensive study on the shear wave VSP data.

The other source of information for determining the azimuthal anisotropy direction is zero-offset shear VSP data. There are two available zero-offset shear VSP datasets in well HMU 13-2: the March 2008 and December 2008 data. The geophones were located in the borehole from 933 m to 1818 m MD in March 2008 survey and from 453 m to 1128 m MD in the December 2008 survey (Figure 5). The significant depth spread helps to follow the azimuthal anisotropy direction from shallow layers to deep layers. The zero-offset VSP data (March 2008 data) were acquired using a horizontal vibrator that made shear waves polarized in only the E-W direction. Because of the single source polarization, 2C rotation analysis [19] was used to estimate the polarization directions of fast and slow shear waves (Figure 11(c)). The second zero-offset VSP data, December 2008, were acquired using one vibrator that made mutual shear wave motion in the E-W and N-S directions. These crossed-dipole sources enabled us to apply Alford rotation analysis [20] to this data. The direct application of Alford rotation analysis is to find the natural polarization direction through a mathematical operation, assuming that the two split shear waves are orthogonally polarized and the anisotropy direction in the subsurface remains constant with depth [21]. The data were rotated from −90 deg to +90 deg (10 deg angle increment), and the RMS amplitude was calculated using a 200 ms time window around the first arrivals. Using (1), the final map (Figure 11(b)) was constructed to point out where the diagonal components (i.e., and ) had more energy than off-diagonal components (i.e., and ):

On average, the March 2008 and December 2008 surveys suggested 134 deg as the fast shear wave polarization direction.

To understand spatial variation of the S11 polarization direction in the field, multioffset VSP data, using Alford rotation analysis, were analyzed. Alford rotation analysis on multioffset VSP data (i.e., near-offset, mid-offset, and far-offset data) depicts that the fast shear wave polarization direction is approximately 138 deg (Figure 12). Similar to the zero-offset VSP results, multioffset VSP data show that the fast shear wave polarization direction does not change significantly with depth.

Analyzing the anisotropy direction using different tools (Figure 13) suggests a 30 deg change (from 108 deg to 138 deg) that may happen because of factors such as different radius of investigation, difference in measurement frequency, complex heterogeneity of the reservoir and overlying layers, among others. A discussion of these factors is outside the scope of this paper, but the data exhibited in Figure 13 indicated that analysis of different sources of information is required to understand azimuthal anisotropy.

5. Wavefield Separation and Imaging

Wavefield separation was the next step in the VSP processing flow. Generally, different wave types can be recognized in VSP data better and more accurately than in surface seismic data. Because of this ease of recognition, we are usually able to effectively filter undesirable waves from VSP data. Because VSP data in Postle Field will be used to image details inside the sand channels (e.g., braid bars and channel boundaries) around the borehole, we need to preserve discontinuities in the data. Median filtering offers the unique ability to preserve discontinuities in data without smearing discontinuities [22]. Although the median filter preserves discontinuities in the data, it diminishes some weak upgoing waves and introduces high frequency noise where upgoing waves interfere with strong downgoing waves. Because some geophones are close to the Morrow A sandstone, reflections (upgoing waves) from this horizon are covered by strong direct arrivals and their multiples (downgoing waves). As expected, the quality of the separated upgoing waves in the interference area is affected by the median filter and the wavefield separation process. To overcome this issue, synthetic zero-offset VSP data were built, and different parameters and approaches for wavefield separation were examined using different median filters.

In one approach, upgoing waves were aligned and they were separated using a median filter. The quality of these separated upgoing waves was not good, but after stacking this result with the separated upgoing waves using the common approach (i.e., removing downgoing waves from data), we obtained a higher signal-to-noise ratio and preserved weak reflections better (Figure 14). We called this method “Dual Median Filtering.” Because the method is adapted for zero-offset VSP data, we applied it to only the zero-offset VSP datasets. The multioffset VSP dataset was filtered by using a common median filtering. F-K filtering is another popular method for separating different wavefields but, due to the high level of noise in the data, discriminating different wavefields in the F-K domain is challenging. Also, F-K filtering cannot preserve sharp discontinuities in the data. After wavefield separation, the downgoing waves were used for making a deconvolution filter. Corridor stacks were made of the zero-offset compressional and shear wave datasets after applying the deconvolution filter (Figure 15).

VSP migration and VSP-CDP transformation are two common approaches for final VSP imaging. Lack of data fold in our VSP survey makes VSP migration method ineffective. An alternative for VSP migration in this case is the VSP-CDP transformation. A VSP-CDP transformation transforms VSP wave data from the depth-time domain to the offset-time domain [23]. This technique uses 2D ray tracing to map reflected signals into the offset-time domain [24]. After imaging the data using VSP-CDP transformation method, we calculated the fast shear wave RMS amplitude for a 30 ms window centered around the Morrow A sandstone (Figure 16).

6. Results and Discussion

Figures 15(c) and 15(d) compare the prestack zero-offset compressional and shear wave images as well as their corridor stacks. Although the Morrow A sandstone is visible on both datasets, it has higher amplitude when shear waves are used. It seems that the deconvolution filter has a considerable effect on compressional wave data and the Morrow A sandstone is not visible before applying this filter. Deconvolution filtering also assists us in resolving the Morrow A sandstone from the underlying layer in the shear wave data. Also, qualitatively, the zero-offset compressional and shear wave data have approximately the same temporal resolution.

Figure 16 displays the shear wave RMS amplitude map for a 30 ms window centered around the Morrow A sandstone from some shots on the northern part of the VSP well. Based on the Postle Field geology and the size of the high-amplitude areas, one hypothesis is that these features represent braid bars which are good quality and thick reservoir rocks. Also, the importance of the VSP data at this area is to map the edge of the channel (mainly NW-SE at this part), and some linear features could show the evolution of the channel edge during time too. At the far end of the map (data were recorded by geophones from 1 to 16, located from 933 m to 1158 m MD in the well), very high amplitude events could result from braid bars as well as strong noise (pay attention to arc shape of the high amplitude area). One reason for the high level of noise at this interval is a weak bond between the casing and the borehole wall. In addition, the presence of shale layers in the Morrow A sandstone can affect the RMS amplitude of the Morrow A. Some low-amplitude areas correspond to the presence of the mentioned shale layers and indicate low quality (or shaly sandstone) reservoir rocks. Unfortunately, as stated before, the data fold at this area is low and is the main limitation for the VSP image in this area. More 3D VSP coverage could image the Morrow A better and more reliably away from the VSP well.

Comparison of surface seismic and VSP RMS amplitude maps reveals advantages of the VSP data. Figures 17(a) and 17(b) represent the Morrow A sandstone S11 RMS amplitude map of the surface seismic (acquired on March 2008) and VSP data. Three amplitude anomalies are highlighted by circles showing similarities and differences between these two datasets. The most interesting anomaly is highlighted by blue circles on both images. The VSP data appear to map details better. The ability of the shear VSP data to delineate the extension of this anomaly reconfirms the advantages of the S-wave VSP technique. As mentioned before, due to high heterogeneity in the reservoir, preserving original amplitudes and amplitudes variations is crucial. The ability of VSP data, after applying the proposed processing techniques in this paper, to delineate high-amplitude anomalies can reduce drilling risk.

7. Conclusions

Shear wave VSP data were processed using a modified processing flow. The main purposes of the proposed processing flow and some introduced techniques were to enhance the weak seismic response from the Morrow A sandstone in the presence of noise as well as to preserve amplitude variations caused by lithology changes (especially between the Morrow A sandstone and the Morrow shale). 3C amplitude analysis made VSP receiver reorientation more robust and avoided some processing artifacts. Also, estimating fast shear wave polarization direction using different tools addressed the scaling issue. The results of this analysis demonstrate that reliance on one source of information, either well logs or seismic data, can introduce significant uncertainty in final results. Dual median filtering was another method which was introduced to enhance VSP weak reflections where upgoing waves interfere with the strong downgoing waves.

Zero-offset VSP results reconfirmed the fact that the Morrow A sandstone signature is weak in compressional wave data. Also, it confirmed the ability of shear wave data to image the thin Morrow A sandstone. The Morrow A sandstone was mapped using the VSP-CDP transformation. High-amplitude anomalies in the shear wave RMS amplitude map could correspond to geologic features like braid bars, which fits the sediment deposition history of Postle Field. Comparison between surface seismic and VSP data showed the potential of VSP data in imaging detailed changes in subsurface rock properties and delineating anomalies with higher spatial resolution. Therefore, shear wave VSP data are able to map detailed geologic features where the heterogeneity in the Morrow A sandstone plays a major role.

Acknowledgments

The authors would like to acknowledge Tom Bratton from Schlumberger, Bruce Mattocks from CGGVeritas, and Richard Van Dok from HiPoint Reservoir Imaging for their advice. Also, they would like to thank their colleagues Dr. Robert Benson, Paritosh Singh, and Mohsen Minaei from RCP for their technical support. Special thanks to Whiting Petroleum for providing the multicomponent data. The authors would like to sincerely thank the Reservoir Characterization Project sponsors and students too.