By: Rocky Roden and Mike Forrest, Rose & Associates, Rodger Holeywell, Marathon Oil Company
Published with permission: The Leading Edge
In the process of quantifying resources/reserves, geoscientists attempt to employ all the available pertinent information to produce the most accurate results. The presence of direct hydrocarbon indicators (DHI) on seismic data can have a significant impact on the reserve/resource calculations not only for volumes, but also uncertainty levels. In 2001 a consortium of oil companies was organized in an attempt to understand seismic amplitude anomalies interpreted as DHIs and their impact on prospect risking and resource calculations (Roden et al., 2005; Forrest et al., 2010). The geologic setting, seismic and rock physics data quality, DHI characteristics, and calibration of drilling results are all incorporated into a database in a consistent and systematic process. From this process, the evaluation of 217 prospects and associated well results has enabled an in-depth understanding of the relevant key aspects of seismic amplitude anomalies and how they relate to drilling results.
DHI Consortium background
In the DHI Consortium, 41 oil companies have contributed to developing a seismic amplitude analysis process that systematically provides a probability of geologic success (Pg) values and associated implications for resource evaluations (Figure 1). The goals of the DHI consortium over the past 11 years have been and continue to be the following:
- Gain a better understanding of how seismic amplitudes impact prospect risking (i.e., predrill chance of geological success or Pg).
- Objectively characterize DHI observations using documented occurrences of recoverable hydrocarbons in the subsurface via prospect reviews and risk analysis discussions.
- Archive a statistically significant library of drilled prospect results.
- Use the prospect database to improve the predictability of Pg and range of uncertainty of hydrocarbon volumes.
- Improve the DHI work analysis process to help risk seismic amplitudes and provide an educational tool for interpreters.
- Discussion and review of pertinent technologies for amplitude interpretation.
These companies have supplied 217 prospects from around the world with approximately an equal number of successes, dry holes, and objective reservoirs ranging in age from Triassic to Pleistocene. The prospects range in size from less than 100 acres to greater than 10,000 acres and depths from 2000 to more than 20,000 ft. This database is predominantly made up of exploratory wells with only 14% coming from “development” or “extension to known reservoir” wells. The closure types include structural (59%), stratigraphic (26%), and combination structural/ stratigraphic (15%) traps.
Critical in the evaluation of DHI characteristics is an understanding of the geologic setting. Because expected DHI responses vary depending on the geologic setting, the prospects in the database are categorized by AVO classes 1–4 as displayed in Figure 4 (Rutherford and Williams, 1989; Castagna et al., 1998). However, 76% of the prospects are class 3 and 22% are class 2, with only a few from class 1 and 4 settings.
A general description of the input information and quantified output for this database as is described in Table 1.
Table 1. General input and quantified output from DHI analysis process. Note that the input parameters are not necessarily linked directly to the corresponding quantified output row.
Figure 2. Locations around the world of the 217 wells in DHI Consortium database
There are several ways the most important DHI characteristics can be determined from the well results. The following list is a result of when a DHI characteristic was rated to be a 4 or 5 on a scale of 1–5, with 1 being the worst and 5 the best. These results were then correlated with whether the wells were successful or dry holes. This approach eliminates any scaling or weighting effects and compares only whether the DHI characteristic is strongly correlated with the actual postdrill well results. The top five DHI characteristics have been determined for both class 2 and 3 AVO categories, which relate to different geologic settings, and for which there are a sufficient number of example prospects to be meaningful.
Figure 3. Prospect types of wells in the DHI Consortium database.
A class 3 AVO setting relates to an unconsolidated gas sand with a gross interval velocity usually less than 8500 ft/s (2650 m/s). These reservoirs typically are low-impedance sands encased in relatively higher-impedance shales. Gas sands in this setting exhibit high amplitudes on the stack section as well as all offsets and angles. The seismic response from the top of these reservoirs produces high amplitudes with respect to background and increase (at times only slightly) with increasing offset/angle (Figure 4).
Top five class 3 DHI characteristics
- Amplitude downdip conformance (fit to closure) on stacked or far offset seismic data. This is the top-rated DHI characteristic in the database and describes how well the seismic amplitudes of the anomaly conform to the downdip structural contours of the reservoir. Key issues are whether the structural contours are in time or depth, potential velocity variations that can affect imaging, the accuracy of how the amplitudes were picked (top, bottom, or window) at the anomaly, and consideration of any stratigraphic variations. It is important to remember that this characteristic relates only to the downdip edge of the amplitude anomaly and how it compares to the interpreted structure. It does not relate to an amplitude’s conformance on the flanks of stratigraphic traps. This lateral conformance characteristic in stratigraphic traps depends on the geologic model and the strength of the independent evidence that supports that model.
- Phase or character change at the downdip edge of the anomaly. This DHI characteristic relates to a change in character such as phase or frequency at the downdip edge of the anomaly. As the model in Figure 6 displays, the transition from the hydrocarbon leg to the water leg in a sand can produce a seismic response indicative of the edge of a reservoir.
- Amplitude consistency in the mapped target area. The internal consistency of the seismic amplitude in a reservoir was found to be a significant DHI characteristic. This relates to the uniformity of the amplitude response within the mapped target area as interpreted from the stacked seismic data (Figure 7). When evaluating this characteristic, consideration should be given to possible faulting and stratigraphic changes that may modify internal consistency.
Figure 4. (top) AVO classes based on amplitude change with offset from the top of gas sands (Rutherford and Williams, 1989; Castagna et al., 1998). (bottom) Angle gathers, zero-offset, and stack responses for typical class 2 and 3 AVO classes.
- Flat spots. This characteristic represents the seismic response at a hydrocarbon contact that presumably is relatively flat. This contact can be at the gas/oil contact, oil/ water contact, or gas/water contact. The gross reservoir thickness of the hydrocarbon unit must be greater than the tuning thickness (vertical resolution) to image a flat spot. Extreme care must be taken because the base or edge of channels, low-angle faults, diagenetic boundaries, or even processing artifacts are often misinterpreted as flat spots. Low-saturation gas in a reservoir can also produce a flat spot and this possibility is usually assessed by regional geologic studies.
- Excluding possible stacked pays, the AVO response is anomalous compared to events above and below. On prestack data, usually gathers, this characteristic refers to whether there are high-amplitude events above and below the targeted anomaly that look similar (Figure 8). The reasoning is that the anomaly is relatively unique suggesting a hydrocarbon-bearing reservoir.
A class 2 setting contains gas sands more consolidated than class 3 sands with gross interval velocities usually between 8500 ft/s and 12,000 ft/s (2650 m/s and 3650 m/s). The acoustic impedances of these gas sands and encasing shales are about equal. The intercept or near-offset amplitude can vary from a weak positive to a weak negative. The AVO effect can be strongly more negative with offset (large gradient) in these settings (Figure 4).
Top five class 2 DHI characteristics
- Amplitude downdip conformance (fit to closure) on far-offset seismic data. The principles behind this characteristic are the same as the top characteristic for class 3 settings, but for class 2 prospects this is usually evaluated on the far offset seismic data.
- Consistency in mapped target area (typically on gathers, faroffset/angle stacks, or windowed attributes). Internal consistency of the seismic amplitude identified on far-offset data was found to be an important DHI characteristic for class 2 rocks just as it is for class 3 rocks.
- AVO observations using gathers, far-offset/angle stacks, or windowed attributes. This characteristic relates to an interpreter’s confidence that the AVO response is proper for a class 2 setting. In other words, the near-offset is low in amplitude (small peak or trough) and the amplitude increases in negative amplitude with offset. (Figure 4) Noisy gathers, incorrect NMO corrections, multiples, insufficient offset during acquisition, and processing artifacts often complicate evaluation of this characteristic.
- The AVO event is anomalous compared to the same event outside the closure. This characteristic describes whether the class 2 AVO response is unique compared to the correlative event outside the closure. This characteristic is often interpreted from far-offset/angle stacks, as well as, intercept, gradient, intercept x-gradient, far-near, and (farnear) x-far displays.
- Change in AVO compared to model (wet versus hydrocarbon- filled). Modeling of the AVO response, usually applying Gassman’s equation, typically involves substitution modeling of gas, oil, and water responses. A comparison of the in-situ and modeled responses to actual gathers provides confidence that hydrocarbons are present or not.
Figure 5. Examples of amplitude conformance to the downdip structural contours from grade 1 (worst) to grade 5 (best).
Reasons for failure
Figure 6. Model of hydrocarbon and water legs of a reservoir and location of possible phase or frequency change at edge of fluid contact (courtesy Quantum Earth Corp).
The consortium defines successful wells as geologic successes (i.e., wells with flowable hydrocarbons). However, in this database all successes were determined to also be commercially successful except for a few wells. These exceptions contained from 50 to more than 100 ft of gas, but were determined to be noneconomic because of their location in the world and lack of infrastructure. Wells containing only low-saturation gas are considered by the consortium to be dry holes even though they may represent successful predictions of acoustic response.
Wet sands in the target interval accounted for 49% of the dry holes in the database. Nearly half of these wet sands were thick. In fact, these thick wet sands were greater than the predrill P10 estimates for reservoir thickness. Hard shale on top of a wet sand, blocky sands with higher than expected porosities, and tuning effects of wet sands all produced amplitude anomalies misinterpreted to be direct hydrocarbon effects.
Low-saturation gas (LSG) accounted for 23% for the failures in the database. This well known phenomenon is caused by a small percentage of residual gas in a reservoir (5–10%) producing an acoustic effect similar to commercial saturation and is usually associated with a breach or break in reservoir seal. All but one of the LSG dry holes were in class 3 rocks, with one in a class 2 setting. All but one of the LSG dry holes were in normally pressured sediment columns. Four of the LSG wells exhibited acoustic flat spots indicative of a paleo hydrocarbon contact. Low vertical effective stress caused by undercompaction accounted for several LSG results in deepwater settings, but was also found where a relatively young sediment column overlaid the amplitude anomaly. Even though not dry holes, several wells exhibited hydrocarbons, usually gas, overlying a LSG column in the sand.
No reservoir present was responsible for 17% of the dry holes. Amplitudes misinterpreted to be hydrocarbons included low-density shale, ash, coal beds, top of hard overpressures, seismic processing artifacts, mud volcano, and diagenetic boundary.
Tight reservoirs accounted for 11% of the dry holes in the database. Examples of these include oil-charged marl, gas-charged condensed section, low-impedance siltstone/ mudstone and low-permeability reservoir.
Figure 7. Examples of amplitude consistency within a defined DHI anomaly from grade 1 (worst) to grade 5 (best).
Implications for resource calculations in exploration
The presence of direct hydrocarbon indicators can have a significant impact on methods to compute prospective resources, especially in exploration settings. How do DHIs change area determination for resources? How do DHIs impact thickness determinations in resource computations? How does a seismic-amplitude (DHI) defined area and thickness relate to a geologically defined area and thickness? The presence of true direct hydrocarbon indicators indicates that all geologic chance factors (e.g., source, migration, seal, reservoir, and trap) are working to form a petroleum system; therefore, the critical issue is determining the confidence level that the seismic is truly displaying a direct hydrocarbon indicator because not all anomalous amplitude events are DHIs.
Figure 8. An unmuted CDP gather displaying an AVO anomaly that appears prominent compared to the events above and below (Roden et al., 2005).The near vertical lines represent offset angles (purple = 10°, red = 20°, green = 30°).
Figure 9. Illustration of DHI Index method to define and weight geologic and amplitude areas for resource calculations
For area determinations, one approach is to consider two separate distribution estimates. One distribution is based on geologic evidence, independent of the amplitude as a DHI, while the other distribution is based on the amplitude-defined area (Figure 9). Each distribution can be defined by a P10 (largest reasonable) and a P90 (smallest reasonable) value. The variance or P10/P90 ratio of the area is often large for the geologic distribution, but much smaller for the amplitude-defined distribution because the anomaly usually has a sharp cutoff. From statistics based on the DHI consortium results, a “DHI Index” has been calculated that indicates the likelihood that the interpreted amplitude anomaly is truly a direct hydrocarbon indicator. This “DHI Index” can be used as a relative weighting factor for the two distribution estimates.
Figure 10. DHI Index values color-coded by drilling outcomes. Note that above 20% almost all wells are successful. There are negative DHI Index values, usually associated with dry holes, where DHI characteristics should be present, but were absent.
The histograms of Figure 10 show that, when the DHI Index is more than 20%, essentially all wells are successful. To get a DHI Index of 20% or higher requires having numerous positive DHI characteristics. Therefore, using this logic, prospects with a DHI Index of 20% or larger employ the DHI-defined area, whereas, a DHI Index of 0% or less employ the geologically defined area (Figure 9). A DHI Index of 0% or less indicates there is no DHI element at all in the area determination. When a prospect has a DHI Index between 0% and 20%, the area is defined by a simple linear weighting (e.g., a DHI Index of 10% would have a 50/50 weighting of geology and amplitude-defined areas) between the two distribution estimates.
For thickness determinations, similar to the DHI Index approach for area, two thickness distributions can be estimated. One is defined geologically and one by the amplitude-de- fined thickness. The geologically defined distribution should be based on well control, sand studies, depositional trends, and usually involves constructing regional and local sand isopachs consistent with the depositional model in terms of reservoir shape and thickness. From the combination of sand isopachs, depositional environment determinations, models and stratigraphy studies, an estimate of the P10 and P90 average net pay thicknesses can be calculated.
For the amplitude-defined thickness, the interpreter first needs to determine whether the anomaly is above tuning, below tuning, or both. Above tuning, the thickness is a function of the time separation between the top and bottom of the anomaly. Gross pay maps (isochrons) can be generated for the P10 and P90 cases including any geometric factors. An estimated net-to-gross and interval velocity are applied to determine the final thickness distributions. Below tuning, in the ideal situation, the composite amplitude from the top and bottom of the bed decreases linearly with thickness starting from tuning thickness. The amplitude should be calibrated with well control and stratigraphy studies. Below tuning, a starting point P10 could be the tuning thickness and a P90 could be the detection limit (approximately 1/30 wavelength of the dominant frequency). Appropriate interval velocities may need to be applied for thickness conversions. Below tuning, the composite amplitude usually incorporates any net-to-gross effects. Once a geologically defined thickness distribution and a seismic amplitude-defined thickness are determined, the DHI Index can be used to determine the weighting of these two distributions. It can be more difficult to separate and differentiate between a geologically and amplitude-defined thickness as compared to interpreting distributions for the area. In this situation, it may be more appropriate to combine geologic and amplitude thickness information and use one distribution without weighting to compute resources. Spectral decomposition has been found to be helpful in certain situations to help determine thickness trends.
Employing the DHI Index method to determine area and thickness for resource calculations is one approach that gives credit to both geologically and amplitude defined distributions. However, it is not uncommon in strong amplitude driven plays to employ the amplitude-defined distributions only, especially for area. The logic here is that the presence of amplitudes is defining the reservoir and the geologic distributions do not apply. Care should be taken in that all amplitudes are not DHIs and not all geologic settings exhibit DHIs equally. In addition, even if there is evidence that the area can be narrowly estimated, there will often be significant remaining uncertainty in resource size because of the variance in average net pay thickness and recovery factor estimates.
Implications for reserve calculations
The U. S. Securities and Exchange Commission (SEC) defines “proved oil and gas reserves” in part as “those quantities of oil and gas, which, by analysis of geosciences and engineering data, can be estimated with reasonable certainty to be economically producible— from a given date forward, from known reservoirs, and under existing economic conditions, operating methods, and government regulations — prior to the time at which contracts providing the right to operate expire, unless evidence indicates that renewal is reasonably certain, regardless of whether deterministic or probabilistic methods are used for the estimation.” Unlike resource calculations in exploration, reserve calculations have had a well or wells drilled to confirm geologic conditions and to some degree seismic amplitude features. Therefore the interpretation of DHI features relates to the determination of “reliable technology” as defined by the SEC. As defined by the SEC, “reliable technology is a grouping of one or more technologies (including computational methods) that has been field tested and has been demonstrated to provide reasonably certain results with consistency and repeatability in the formation being evaluated or in an analogous formation.” Therefore for DHI characteristics to qualify as reliable technology requires the correlation of amplitude features with well control, reservoir engineering tests, known geologic trends and stratigraphy.
Even with well control and substantial geologic control of thickness and area extents, care should be taken in incorporating DHI information from seismic data. Considerations for amplitude anomalies and reserve determinations include:
- Strong amplitude is primarily at wedge where tuning occurs, weaker amplitude updip
- Stratigraphic trap where conventional downdip and updip limits do not apply
- Reservoir thins updip, bald structure
- Velocity tilt can affect conformance to closure interpretations
- Faults affecting amplitudes for area determinations
- Stratigraphic changes across field
- Masking of amplitudes from bright spots above prospect
- Amplitude effects from low-saturation gas downdip of proven pay in the same reservoir
- Beware of large column heights for gas reservoirs
- Data issues, such as prospect at edge of survey, complicate area estimates
In development, seismic technologies such as spectral decomposition, neural networks, statistical approaches, and advanced seismic inversion techniques are often applied to help determine reserves. The transition from exploration to development and the use of seismic inversion requires the incorporation of all available geologic, seismic, and engineering data in an attempt to discriminate specifically between lithology, porosity, and fluid effects and how these parameters relate to reserve calculations.
Over the past 11 years, 41 oil companies have contributed data for the evaluation of more than 200 drilled prospects (including assessment of the associated technical information and drilling results). The subsequent database developed from these assessments has enabled the identification of the pertinent aspects of direct hydrocarbon indicators and their impact on risking and resource/reserve determination. This database is composed primarily of exploratory wells (86%) and has direct relevance for computing resources in exploration. However, the lessons learned from this database are also pertinent for reserve assessment in development because an understanding of the most important DHI characteristics, reasons for failure, and potential pitfalls in DHI interpretation provide valuable information in the reserve assessment process. Critical in the evaluation of DHI characteristics is the pertinent technical approaches employed and their calibration to known drilling results. However, even in development settings with a certain amount of well control for calibration, DHI interpretation requires a consistent and systematic evaluation procedure, including recognition of potential pitfalls and an understanding of the cause of the DHIs and their significance in reserve assessment.
Castagna, J., H. Swan, and D. Foster, 1998, Framework for AVO gradient and intercept interpretation: Geophysics, 63, no. 3, 948– 956, http://dx.doi.org/10.1190/1.1444406.
Forrest, M., R. Roden, and R. Holeywell, 2010, Risking seismic amplitude anomaly prospects based on database trends: The Leading Edge, 29, 936–940, http://dx.doi.org/ 10.1190/ 1.3422455.
Roden, R., M. Forrest, and R. Holeywell, 2005, The impact of seismic amplitudes on prospect risk analysis: The Leading Edge, 24, no. 7, 706–711, http://dx.doi.org/10.1190/1.1993262.
Roden, R., J. Castagna, and G. Jones, 2005, The impact of prestack data phase on the AVO interpretation workflow-A case study: The Leading Edge, 24, 890–895. http://dx.doi.org/10.1190/1.2056369.
Rutherford, S. E. and R. H. Williams, 1989, Amplitude versus off- set variation in gas sands: Geophysics, 54, no. 6, 680–688, http:// dx.doi.org/10.1190/1.1442696.
U. S. Securities and Exchange Commission, 2008, Modernization of oil and gas reporting, December 31.
Acknowledgments: The authors thank the member companies of the DHI Consortium for providing invaluable information necessary to develop the resulting interpretation process and prospect database.