Subsurface maps

Series Development Geology Reference Manual Methods in Exploration Geological methods Subsurface maps K. W. Weissenburger Web page AAPG Store

Reservoir properties are mapped to promote optimal field development. Subsurface maps dictate well placement and enable engineers to calculate reserves and monitor trends in reservoir performance. Geologists play a key role in subsurface mapping by using interpretations of depositional environments and diagenetic events to project reservoir data away from relatively few well control points. In this sense, subsurface mapping is in great contrast to geological mapping of the earth's surface. Whether using traditional concepts[1] or “high technology” computer contouring hardware/software systems,[2] mapping interwell areas places a premium on interpretation rather than straightforward plotting of precise data. “Mapping” is here limited to projections in plan view.

Mapping surfaces

A number of surfaces are typically mapped during reservoir development to show closure and other limits to reservoir production. Maps of top of pay and bottom of pay can also be “subtracted” to determine pay thickness.

Figure 1 Structure map of the top of the T5 marker, Frlo Formation, Brazoria County, Texas. (After Bebout el al.)[3]

Structure

Structure maps show lines of equal elevation or depth for a selected marker horizon (Figure 1) (see Evaluating structurally complex reservoirs). Mean sea level is a useful reference datum. Commonly contoured horizons are top of zone or top of net pay. Control points are provided by surveyed wells and can be supplemented by seismic interpretations, especially offshore. In highly developed fields, typically onshore, sufficient well control might exist to allow geostatistical interpolation between control points.

Figure 2 Structure of the base of the Humber unconformity (top of the Brent Group), Dunlin field, U.K. Northern North Sea mapped with 1979 and 1989 vintage data. Contours are marked in ft subsea × 100; contour interval is 100 ft. (From Braithwaite et al.)[4]

Fault planes

Faults are special surfaces whose traces will show on structure contour maps (Figures 1 and 2). Faults form bounding surfaces for some reservoirs, and sufficient well control might exist to contour map the fault surface itself. Projections of subsurface data into the plane of the fault are also useful “maps” for reservoir development, but are more appropriately described as cross sections. (For details of construction of fault plane maps, see Conversion of well log data to subsurface stratigraphic and structural information.)

Figure 3 Map of pressure response to pattern flood, Judy Creek field, western Canada, 1974 and 1975. Contour interval is 2750 kPa. (After Jardine et al.)[5]

Unconformities and subcrops

Surfaces of unconformity can be especially useful marker horizons for structure contour mapping (Figure 2). In many fields, unconformities are the location of sealing shales and/or source rocks above reservoir pay. Subcrop maps, traces of productive zones, barriers, or marker horizons mapped on the unconformity surface are invaluable for planning well placement and for reservoir development.

Pressure

Maps of reservoir pressure are useful throughout reservoir life (Figure 3). Pressures should be converted to a common depth datum, such as mid-reservoir, prior to contouring. (For information on obtaining pressure data, see Production testing and Pressure transient testing, Wireline formation testers, and Drill stem testing.)

Mapping thicknesses

Figure 4 (a) Cross section and (b) net pay Isopach map of the Strachan gas field, western Canada. Contour Interval is 100 ft. (From Hriskevich et al.)[6]

Interpretations of depositional trends, pre- and syndepositional structural development, and reservoir storage capacity are based in large part on thickness information. An accurate meaning of thickness is critical in these and other analyses (see Conversion of well log data to subsurface stratigraphic and structural information).

Isopach

A contour map of equal values of true stratigraphic thickness is an isopach map (Figure 4). Except for vertical wells in horizontal beds, corrections for wellbore deviation and formation dip are needed to make isopach maps.

Isochore

A contour map of equal values of true vertical thickness is an isochore map.[7] Note that in common practice, isochore maps are informally referred to as “isopach” maps, a term that properly should be restricted to true stratigraphic thickness.

Isochron

An isochron map is a contour map of equal values of seismic traveltime between selected events.[7] Isochron maps are the seismic analog of isochore maps and, as such, are intended to derive thickness information from seismic data. Isochroning between events above and below a pay horizon, for example, would estimate pay thickness. Renick and Gunn[8] present a good case history of using isochron and time-structure maps to generate “isopach” and elevation-structure maps. Their isochron-isopach approach delineated reef trends for further development drilling and used well penetrations through a shallow horizon for depth control on a deeper horizon. Phipps[9] documents the pros and cons of using isochron thins and structural highs as exploration drilling criteria for dolomitized Devonian limestones.

Mapping to calculate reserves

When few production performance data are available, typically early in the life of a reservoir, reserves can be calculated by a volumetric analysis (see Reserves estimation). For an oil reservoir, the basic volumetric equation is as follows:

$N=7758\times A\times H\times \phi \times (1-S_{{{\rm {w}}}})/{\mbox{Boi}}$

where

• N = original oil in place, stock tank barrels (STBstock tank barrels)
• 7758 = conversion factor (acre-feet to barrels)
• A = area of reservoir (acres)
• H = height or thickness of pay zone (feet)
• ϕ = porosity (fraction of bulk volume)
• Sw = connate water saturation (fraction of pore volume)
• Boi = oil formation volume factor (dimensionless, reservoir bbls/ stock tank bbls)

A similar basic equation applies to gas reservoirs.

Net pay

The product A × H is the reservoir bulk volume, and the product A × H × ϕ is the reservoir pore volume. The general determination of bulk reservoir volume involves mapping reservoir area in plan view and mapping net pay in terms of true vertical thickness to provide a common presentation of dipping beds or deviated wells. An isochore map of net pay should be contoured using well control points and interpolated or extrapolated using available seismic and well test data and the geologist's interpretation of depositional and diagenetic history.

“Net” pay (see Effective pay determination) implies that some formation thickness has been excluded from consideration by either (1) occurring below an oil-water contact (or above a gas-water contact), or (2) having porosity and/or permeability values below a “cutoff” limit for productivity. Not all net pay is necessarily productive at a given well spacing. Discontinuous productive horizons between wells might be described, for example, by the concept of net pay to net connected pay ratio.[10]

Porosity

Figure 5 Porosity-weighted average water saturation map for Layer 2 of a Middle Eastern carbonate reservoir.

The porosity (ϕ) in a reservoir zone can be determined from log and/or core data. The data in an individual well within the net pay interval can be averaged arithmetically and posted on a map for contouring. The averages should be weighted by thickness.

Water saturation

The water saturation (Sw) within the net pay interval is typically estimated from well logs. Water saturations can also be derived from capillary pressure testing of cores to determine the relationship of water saturation versus height above the oil-water contact. Like porosity, the water saturation data in an individual well within the net pay interval can be averaged arithmetically and posted on a map for contouring (Figure 5). The averages should be weighted by porosity.

Oil saturation

In an oil-water system, the water saturation and oil saturation (So) sum to 1. Therefore, once Sw has been determined, oil saturation can be calculated and mapped as So = 1–Sw.

Mapping for reservoir management

A variety of maps are used to predict or monitor reservoir performance.

Figure 6 Porosity thickness (ϕH) maps for the B and C zones from the San Andres Formation reservoir, Jordan field, Ector and Crane Counties, Texas. Contours in PV fraction-feet. (After Major et al.)[11]

Permeability

Permeability (k) can also be mapped and contoured (see Core-log transformations and porosity-permeability relationships). As for saturation values, some care must be exercised in mapping permeability because values must be derived from indirect measurements. Typically, permeabilities are derived from wireline log porosities transformed on the basis of core permeability versus porosity cross plots. Permeabilities can be reported at ambient laboratory conditions of pressure or adjusted to reservoir conditions of confining pressure. Similarly, permeabilities can be absolute permeabilities to air (nitrogen) or liquid or effective permeabilities]] to oil in the presence of irreducible water. Permeability values in an individual well are thickness weighted and typically averaged harmonically, arithmetically, or geometrically, depending on flow geometry. Alternatively, flow capacity (kH) values derived from pressure transient testing can be divided by net pay thickness (H) to yield a liquid permeability value for a well.

Porosity thickness

Reservoir storage capacity or porosity thickness (ϕH) is the product of porosity and net pay (Figure 6).

Productivity index

To avoid coning, sand production, pipe collapse, or other harmful effects, wells might not be produced at their maximum wide-open flow rates. Therefore, the ability of a well to produce is usually determined by a productivity index (PI).[12] The PI is a measure of the stock tank barrels (STBstock tank barrels) of oil produced per day per psi drawdown under steady state or pseudosteady state flow conditions. Changes will show on periodic maps of PI during reservoir life indicating trends in reservoir depletion or formation damage.

Solution gas to oil ratio

Engineers forecast ultimate recoverable reserves by applying material balance equations or decline curve analysis to production history records. For example, in a depletion-type reservoir, the solution gas to oil ratio is sometimes plotted versus cumulative oil production on semilog paper.[13] If such a curve shows a good straight-line relationship, the curve can be used to predict the trend of a cumulative gas or cumulative oil plot to estimate ultimate recovery.

The solution gas to oil ratio (GORgas-oil ratio) is the amount of dissolved gas that will evolve from the oil as the pressure is reduced to atmospheric from some higher pressure. GORgas-oil ratio is usually expressed in units of SCFstandard cubic feet gas/STBstock tank barrels oil. A barrel of oil and its solution gas at reservoir conditions of temperature and pressure will usually “shrink” as the fluid is produced and brought to stock tank conditions (normally reported at 60°F288.706 K
15.556 °C
519.67 °R
and 14.7 psia). As GORgas-oil ratio changes during reservoir life, GORs for individual wells can be mapped periodically to identify areas of the reservoir receiving or not receiving pressure support and serving as indicators for reservoir management action.

Figure 7 Cumulative oil production map for the A, B, C, and D zones from the San Andres Formation reservoir, Jordan field, Ector and Crane counties, Texas. Contours in MSTB/year/acre. (After Major et al.)[11]

Water cut

Water cut is the fraction of a liquid production stream that is water, where oil cut = 1 - water cut. Like GORgas-oil ratio, water cut will change during the life of a reservoir, and periodic mapping can serve as a performance indicator for reservoir management. A variety of performance features can be indicated by water cut maps, including water coning, directional permeability or channeling, and formation damage.

Cumulative production

Cumulative oil or gas production is a parameter useful for ultimate reserves forecasts. Cumulative production can also be mapped periodically as a performance indicator signaling areas of the reservoir that may be responding in a manner seemingly unrelated to initial potential. Figure 7 shows an example of cumulative production that was concluded to be only poorly correlated to storage capacity (Figure 6) in individual and summed zones of a carbonate reservoir.[11] In this case, porosity did not necessarily indicate effective porosity.

Other maps

A variety of other maps can come into play during the development of a specific reservoir. Maps of facies, facies architecture for more information, see Lithofacies and environmental analysis of clastic depositional systems and Carbonate reservoir models: facies, diagenesis, and flow characterization), paleoenvironment, and isolithology might be particularly important in selecting stepout well locations and planning reservoir development strategy. Other reservoir properties such as temperature can have value for specific reservoir engineering applications, particularly where potentially temperature-sensitive chemical stimulation, production, or recovery technology might be involved.

References

1. Landes, K. L., 1951, Subsurface maps and sections, in Petroleum Geology: New York, John Wiley, 660 p.
2. Jones, T. A., Hamilton, D. E., Johnson, C. R., 1986, Contouring geologic surfaces with the computer: New York, Van Nostrand Reinhold, 320 p.
3. Bebout, D. G., Loucks, R. G., Gregory, A. R., 1978, Frio Sandstone reservoirs in the deep subsurface along the Texas Gulf coast—their potential for production of geopressured geothermal energy: The Univ. of Texas Bureau of Economic Geology Report of Investigations, n. 91, 93 p.
4. Braithwaite, C. I. M., Marshall, J. D., Holland, T. C., 1989, Improving recovery from the Dunlin field, U., K. Northern North Sea, in Formation Evaluation and Reservoir Geology: 64th Annual Technical Conference of the Society of Petroleum Engineers, San Antonio, TX, Oct. 8–11, SPE 19878, 18 p.
5. Jardine, D., Wilshart, J. W., 1987, Carbonate reservoir description, in Tillman, R. W., Weber, K. J., eds. Reservoir Sedimentology: SEPM Special Publication 40, p. 129–152.
6. Hriskevich, M. E., Faber, J. M., Langton, J. R., 1980, Strachan and Ricinus West gas fields, Alberta, Canada, in Halbouty, M. T., ed., Giant Oil and Gas Fields of the Decade 1968–1978: AAPG Memoir 30, p. 315–327.
7. Tucker, P. M., 1988, Seismic contouring—a unique skill: Geophysics, v. 53, n. 6, p. 741–749, DOI: 10.1190/1.1442509.
8. Renick, H. Jr., Gunn, R. D., 1989, Triangle Ranch Headquarters field development using shallow core holes and high-resolution seismic data: Geophysics, v. 54, n. 11, p. 1384–1396, DOI: 10.1190/1.1442602.
9. Phipps, G. G., 1989, Exploring for dolomitized Slave Point carbonates in northeastern British Columbia: Geophysics, v. 54, n. 7, p. 806–814, DOI: 10.1190/1.1442709.
10. Poston, S. W., 1987, Development plan for oil and gas reservoirs, in Bradley, H. B., ed., Petroleum Engineering Handbook: Richardson, TX, Society of Petroleum Engineers, p. 36-1–36-11.
11. Major, R. P., Holtz, M. H., 1989, Depositionally and Diagenetically Controlled Reservoir Heterogeneity at Jordan Field: 64th Annual Technology Conference of the Society of Petroleum Engineers, San Antonio, TX, Oct., SPE 19874, p. 633–640.
12. Kimmel, J. D., Dalati, R. N., 1987, Potential tests of oil wells, in Bradley, H. B., ed., Petroleum Engineering Handbook: Richardson, TX, Society of Petroleum Engineers, p. 32-1–32-16.
13. Garb, F. A., Smith, G. L., 1987, Estimation of oil and gas reserves, in Bradley, H. B., ed., Petroleum Engineering Handbook: Richardson, TX, Society of Petroleum Engineers, p. 4-1–40-38.