Click to view PDF version with images, as originally presented in AAPG
Explorer.
Stratton Field, South Texas: Example of Integration in Reservoir Characterization
By
Bob A. Hardage
Search and Discovery Article #40020, (2001).
Senior research scientist, Bureau of Economic Geology at the University of Texas at Austin ([email protected] ).
Adapted for online presentation from two articles by same author, entitled “Integration Aids Reservoir Effort” in Geophysical Corner, AAPG Explorer, December, 1997, and “Combining Data Aids Interpretation” in Geophysical Corner, AAPG Explorer, January, 1998. Appreciation is expressed to the author and to M. Ray Thomasson, former Chairman of the AAPG Geophysical Integration Committee, and Larry Nation, AAPG Communications Director, for their support of this online version. For a somewhat expanded version of this B.E.G. study, the viewers are referred to “3-D Seismic Imaging and Interpretation of Fluvially Deposited Thin-Bed Reservoirs,” by Bob A. Hardage, R.A. Levey, V. Pendleton, J. Simmons, and R. Edson, in AAPG Studies in Geology No. 42, p. 27-33. This publication is available from Bookstore at www.aapg.org.
This article summaries a study of the Stratton Field, a large, Frio gas-producing property in Kleberg and Nueces counties in South Texas. The stratigraphic interval involved was the Oligocene Frio Formation – a thick, fluvially deposited sand-shale sequence that has been a prolific gas producer in Stratton Field and in several other fields along the FR-4 depositional trend (Figure 1). This reservoir-characterization effort is an example of integrating geophysics, geology, and reservoir engineering technologies to detect thin-bed compartmented reservoirs in a fluvially deposited reservoir system. An integrated interpretation of a complex reservoir system is made using 3-D seismic, well log, and bottom-hole pressure data. Techniques for thin-bed interpretations and recognition of reservoir heterogeneity are essential elements in the interpretation.
Two examples of seismic thin-bed interpretation in a fluvially deposited gas reservoirs are shown, and these interpretations are supported with geologic and reservoir engineering data. In these examples, the 3-D seismic data reveal stratigraphic variations where reservoir pressure information implies that a compartment boundary should exist. These examples illustrate that, although fluvial deposition creates numerous compartment boundaries, determining which seismically imaged stratigraphic changes are compartment boundaries requires that geologic and reservoir-engineering data (particularly reservoir-pressure data) be incorporated into the seismic interpretation.
In this study it was particularly important to have an accurate and reliable way to translate thin-bed stratigraphy (known in depth) into precisely defined seismic time windows. VSP data, when properly recorded and processed, are the best information to establish the detailed depth-versus-time calibration required to seismically distinguish closely spaced thin-beds. The VSP calibration procedure used was able seismically to distinguish thin-beds that are vertically separated by as little as 4 ms.
Figure
1--Index map showing FR-4 depositional trend in which Stratton Field is located.
Figure
2--Map of Stratton field area, showing well locations and 3-D seismic grid.
Figure
3--Zero-offset image is spliced into a north-south vertical slice from the 3-D
data volume passing through the VSP well (Figure 2). Also shown is a graphic
representation of the stratigraphic column penetrated by the VSP well. The
F37 reservoir is about 20 feet above the F39 reservoir.
Figure
5--Reflection amplitude behavior on the F39 depositional surface shown in Figure
4.
Click here for sequence of Figures 5 and 9 (for comparison of F37 and F39).
Figure
6--Magnified view of this F39 surface in the vicinity of four key wells.
Figure
9-- Reflection amplitude behavior on the F37 surface.
Click here for sequence of Figures 5 and 9 (for comparison of F37 and F39).
Figure
10--Enlargement of the meander features of F37.
Click here for sequence of Figures 10 and 13 (seismic reflection event vs. model).
Figure
11--Log-based stratigraphic cross-section of the F37 reservoir.
Figure
12--Rapid F37 pressure decay in well 189.
Figure
13--Multi-discipline reservoir model, illustrating the proposed model for F37
reservoir.
Click here for sequence of Figures 10 and 13 (seismic reflection event vs. model).
Thin-Bed Interpretation Procedure
Defining Chronostratigraphic Depositional Surfaces
The study covered a 7.6-square-mile area (Figure 2) where 3-D seismic data were acquired and where a large number of wells were used in making a geologic analyses of the Frio reservoirs. Additional well data (with well locations corresponding to the circled dots) were used to supplement the historic well log, production and reservoir pressure databases; they consist of modern well logs, cores, and various pressure tests.
Vertical seismic profile (VSP) data were recorded in two closely spaced wells inside the triangle shown near the center of the 3-D grid.
Thin-Bed Interpretation Procedure
The seismic interpretation at Stratton Field was particularly challenging because most of the Frio reservoirs are thin (<15 feet, or 5 meters), and they were closely stacked--in some areas individually separated vertically by only 10-15 feet (3-5 meters). These conditions required precise calibration of stratigraphic depth-versus-seismic travel-time to extract a depositional stratal surface from the 3-D data volume that would reliably depict the areal distribution of a particular Frio thin-bed reservoir. Zero-offset VSP data recorded in one of the wells shown in Figure 2 were used to establish the precise depth-versus-time control needed for the thin-bed interpretation. Figure 3 shows the zero-offset image spliced into a north-south vertical slice from the 3-D data volume passing through the VSP well. Also shown in Figure 3 is a graphic representation of the stratigraphic column penetrated by the VSP well. Only producing or potentially-producing Frio reservoirs are shown in this diagram, and not all of the reservoirs are labeled by name. The top and base of each reservoir are accurately positioned in terms of two-way VSP travel-time, and because there is no difference in the VSP and 3-D time datum in this case, the reservoirs are also correctly positioned vertically inside the 3-D seismic data volume at the VSP well.
Using these VSP travel-time control data, each thin-bed reservoir can be placed in the correct reflection phase position in the 3-D seismic reflection wavefield at the VSP well. This thin-bed calibration was extended away from the VSP well and across the entire 7.6-square-mile area imaged by the 3-D data using the interpretation principles of seismic stratigraphy.
Defining Chronostratigraphic Depositional Surfaces
The fundamental assumption made in the seismic interpretation was that seismic reflections follow chronostratigraphic depositional surfaces (Vail and Mitchum, 1977). Thus, a continuous seismic reflection even imaged over the entire 7.6-square-mile area by the 3-D seismic data defines a geologic surface that corresponds to a fixed, constant depositional time. Two such areally continuous reflection events were found in the Frio interval. These two surfaces are shown on the east-west vertical section crossing the VSP well (Figure 4).
At the VSP control well, the apex of the peak associated with the shallower stratal surface (the orange surface in Figure 4) corresponds to the thick C38 reservoir (Figure 3), and the apex of the peak at the deeper stratal surface (the green surface in Figure 4) correlates with the F11 reservoir. Thus, the seismic time surface following the apexes of all of the peaks of the orange event is assumed to define the ancient topographic Frio surface at the time when the C38 reservoir sediments were deposited.
Similarly the seismic time surface following the apexes of the peaks of the deeper green event define the ancient depositional surface associated with the F11 reservoir. Once the 3-D data volume was flattened relative to one of these two reference stratal surfaces, it follows that any horizontal time slice in this flattened data volume also followed an ancient Frio depositional surface – as long as the seismic reflection character in the immediate neighborhood of the time slice was time-conformable with the reflection character in the immediate vicinity of the reference surface used to flatten the data volume.
This interpretation is based on the assumption that the entire Frio section inside the 7.6-square-mile grid was seismically conformable to one of the two seismic reference surfaces. In this specific interpretation problem, with many closely spaced (vertically) thin-beds, the VSP-defined position of a particular thin-bed reservoir was rarely at the apex of a reflection peak or trough. Invariably, each thin-bed of interest was positioned at some intermediate, commonly non-descript phase point in the reflection waveform at the VSP control well.
To create a seismic image that emphasizes the internal complex architecture of a given thin-bed reservoir system, the migrated 3-D data volume was:
-
First time shifted so the proper pre-defined reference stratal surface was flat.
-
Then a horizontal time slice was made through this flattened data volume at the exact VSP-defined time for the targeted thin-bed, regardless of where that time slice is positioned in the reflection waveform at the VSP control well.
By prior assumption, the seismic time surface contained in this horizontal slice was the fixed depositional stratal surface where that thin-bed unit was deposited, and any seismic anomalies seen on this surface would be related directly to stratigraphic heterogeneities within the targeted thin-bed and, to a lesser degree, would be related to stratigraphic variations in thin-beds positioned immediately above and below the target thin-bed.
The F39 reservoir was the deepest Frio reservoir studied. (The depositional surface for the F39 reservoir is shown by the yellow horizon in Figure 4; the reflection amplitude behavior on the F39 depositional surface is shown in Figure 5.) The linear north-south trends near the center of the image are assumed to be residual effects from the deeper Vicksburg faults (a magnified view of this F39 surface in the vicinity of four key wells is shown in Figure 6).
F39 reservoir pressure measurements were acquired in all four wells (Figure 7), and the differences in these static pressures indicate that each well is in a different F39 compartment. The 3-D seismic image and the available geologic control gave clues as to where the boundaries are that segregate the F39 reservoir into these distinct compartments.
Figure 8 displays the available geologic control. The log curves infer that the F39 reservoir in each well was deposited in a channel environment that shows some evidence of splay deposition. The seismic image defines some possible compartment boundaries. For example, the most likely cause of the compartment boundary that separates well 197 from the other wells is the depositional variation that created the red/blue (positive/negative) amplitude changes, which trend north-south between crossline coordinates 130 and 140 (Figure 6).
Similarly, a probable seismic indication of the compartment boundary that segregates well 75 from the other wells is the positive-to-negative (red-to-blue) amplitude variations trending north-south between crossline coordinates 110 and 120. By analyzing the seismic, geologic, and engineering data associated with the F39 reservoir, it is possible to detect F39 reservoir compartments seismically--at least in the vicinity of wells 75, 175, and 197.
To create this reservoir compartment model, it is essential that the seismic image be interpreted with the assistance of reservoir-pressure data to infer which of the many stratigraphic changes revealed in the seismic image are most likely to be the compartment boundaries.
The F37 reservoir is approximately 20 feet (6 meters) above the F39 reservoir in the VSP calibration well (Figure 3). The two-way travel-time difference between F37 and F39 is only four milliseconds (4 ms). Using the thin-bed interpretation procedure above, a time slice was made through the flattened 3-D data volume 4ms above the F39 stratal surface. This F37 surface is displayed in Figure 9. Comparing this image with the F39 surface (Figure 5), red, linear north-south apparent channels in the central part of the F37 image are similar to those observed in the F39 image, implying that Vicksburg faulting was still controlling sedimentation in this part of the field.
However, there is a significant difference in the southeast quadrant of the F39 and F37 images. Specifically, meander channel features occur at the F37 level but are not present at the deeper F39 surface. (An enlarged plot of the meander features in Figure 9 is shown on Figure 10.) A log-based stratigraphic cross-section of the F37 reservoir across the meander features and extending southward beyond the seismic grid was constructed (Figure 11). The depositional environment (either channel or splay) at each well is an interpretation based on log-curve shape and was made before the 3-D seismic data were recorded.
This geologically-based interpretation of the F37 depositional environments indicates that the meander feature seen in the F37 seismic surface is indeed a depositional channel. Specifically, the log interpretation (Figure 11) implies the F37 reservoirs found in wells 189 and 185 were deposited as channel fill, and the seismic image shows these wells to be coincident with a meander feature.
The interpretation of depositional environment for the extremely thin F37 reservoir in well 211 is a splay (Figure 11). The 211 wellhead is approximately 300 feet (91 meters) north of the meander feature (Figure 10). The log-based interpretation of the F37 depositional environment at the 211 well is thus supported by seismic evidence.
Pressure histories recorded in several F37 reservoirs near these seismic meander features were analyzed to determine if reservoir compartmentalization exists. These pressure histories (Figure 12) show there are at least three, and perhaps four, individual F37 reservoir compartments in this area of the field.
A proposed reservoir model that honors all three databases--seismic, geologic, and reservoir engineering--is shown in Figure 13. This model assumes that the F37 reservoir in the southeast quadrant of the 3-D grid is composed of three intermeshed channels, labeled A, B and C, and a grid overlay of seismic inline and crossline coordinates is provided so these channels can be correlated with features in the 3-D seismic image.
The location of the F37 stratigraphic cross-section (Figure 11) is shown, but this geologic information defines channel locations along only a single 2-D profile of the model. The important information is the reservoir-pressure data, because without this engineering data there would be no reason to conclude that a three-channel model would be appropriate. Thus, the reservoir-compartment model places well 129 in Channel A and well 185 in Channel B; this allows these two wells to be in different F37 pressure regimes; i.e., in different compartments (channels).
Wells 127 and 161 are proposed to be in channel C, south of the 3-D seismic coverage. Only one meander loop of this hypothesized channel C extends into the 3-D seismic grid. The rapid F37 pressure decay observed in well 189 (Figure 12) implies that this well is not in pressure communication with well 185, even though both wells are in channel B. There may be an intrachannel compartment boundary in channel B.
The reservoir model in Figure 13 is hypothetical and may not yet be the correct picture of the compartmentalized nature of these F37 reservoirs. However, the F37 reservoir in this portion of Stratton Field is segregated into distinct compartments, and this compartmentalization must be caused by the variable elements of fluvial deposition, because the seismic data show no evidence of faulting in these particular reservoirs.
The proposed reservoir model honors all existing data that provide any information about the F37 reservoir system. The seismic image in Figure 10 reveals not just one meander channel system but at least three intermeshed thin-bed channels. By using pressure histories it was possible to use 3-D seismic images to define where compartment boundaries most likely exist in the interwell space.
Raymond A. Levey, Virginia Pendleton, James Simmons and Rick Edson, all at the Bureau of Economic Geology at the time this work was done, assisted in the writing.