Richard R. Hillis

The Australian Stress Map

Knowledge of the in situ stress field of the Australian continent has increased greatly since compilation of the World Stress Map in 1992, principally by analysis of borehole breakouts and drilling-induced tensile fractures in petroleum wells. Stress orientations are variable across the Australian continent as a whole. However, within 15 of 16 individual stress provinces defined in the Australian continent (of one to a few hundred kilometres scale), mean stress orientations are statistically significant. The stress provinces, and stress trajectory mapping, reveal that there are systematic, continental-scale rotations of stress orientation within Australia. Unlike many other continental areas, stress orientations do not parallel the direction of absolute plate motion. Nonetheless, the regional pattern of stress orientation is consistent with control by plate boundary forces, if the complex nature of the convergent northeastern boundary of the Indo-Australian plate, and stress focusing by collisional segments of the boundary, is recognized.

Topography, boundary forces, and the Indo-Australian intraplate stress field

The relative contribution of topographic (e.g., ridge push, continental margins, and elevated continental crust) and plate boundary (e.g., subduction and collisional) forces to the intraplate stress field in the Indo-Australian plate (IAP) is evaluated through a finite element analysis. Two important aspects of the IAP intraplate stress field are highlighted in the present study: (1) if substantial focusing of the ridge push torque occurs along the collisional boundaries (i.e., Himalaya, New Guinea, and New Zealand), many of the first-order features of the observed stress field can be explained without appealing to either subduction or basal drag forces; and (2) it is possible to fit the observed SHmax, (maximum horizontal stress orientation) and stress regime information with a set of boundary conditions that results in low tectonic stress magnitudes (e.g., tens of megapascals, averaged over the thickness of the lithosphere) throughout the plate. This study therefore presents a plausible alternative to previous studies of the IAP intraplate stress field, which predicted very large tectonic stress magnitudes (hundreds of megapascals) in some parts of the plate. In addition, topographic forces due to continental margins and elevated continental material were found to play an important role in the predicted stress fields of continental India and Australia, and the inclusion of these forces in the modeling produced a significant improvement in the fit of the predicted intraplate stresses to the available observed stress information in these continental regions. A central focus of this study is the relative importance of the boundary conditions used to represent forces acting along the northern plate margin. We note that a wide range of boundary conditions can be configured to match the large portion of the observed intraplate stress field, and this nonuniqueness continues to make modeling the IAP stress field problematic. While our study is an important step forward in understanding the sources of the IAP intraplate stress field, a more complete understanding awaits a better understanding of the relative magnitude of the boundary forces acting along the northern plate margin.

The origins of the intraplate stress field in continental Australia

The ridge push force acting on the Indo-Australian plate exerts a significant torque (8.5 × 1025N m) about a pole at 30.3°N, 34.5°E. The angular difference between this torque pole and the observed pole of rotation for the plate (19.2°N, 35.6°E) is less than 12° and suggests that the ridge push force plays an important role in the dynamics of the Indo-Australian plate. We have used an elastic finite-element analysis to study the predicted intraplate stress field in continental Australia for four models which employ different boundary conditions to balance the ridge push torque acting on the plate. The modeling indicates that a number of important features of the observed stress field within the Australian continent can be explained in terms of balancing the ridge push torque with resistance imposed along the Himalaya, Papua New Guinea, and New Zealand collisional boundaries segments. These features include NS-to NE-SW-oriented compression in the northern Australia and E-W-oriented compression in southern Australia. Our analysis also shows that subduction processes along the northern and eastern boundaries provide only second-order controls on the intraplate stress field in continental Australia.

Origin of overpressure and pore-pressure prediction in the Baram province, Brunei

Accurate pore-pressure prediction is critical in hydrocarbon exploration and is especially important in the rapidly deposited Tertiary Baram Delta province where all economic fields exhibit overpressures, commonly of high magnitude and with narrow transition zones. A pore-pressure database was compiled using wireline formation interval tests, drillstem tests, and mud weights from 157 wells in 61 fields throughout Brunei. Overpressures are observed in 54 fields both in the inner-shelf deltaic sequences and in the underlying prodelta shales. Porosity vs. vertical effective stress plots from 31 fields reveal that overpressures are primarily generated by disequilibrium compaction in the prodelta shales but have been generated by fluid expansion in the inner-shelf deltaic sequences. However, the geology of Brunei precludes overpressures in the inner-shelf deltaics being generated by any conventional fluid expansion mechanism (e.g., kerogen-to-gas maturation), and we propose that these overpressures have been vertically transferred into reservoir units, via faults, from the prodelta shales. Sediments overpressured by disequilibrium compaction exhibit different physical properties to those overpressured by vertical transfer, and hence, different pore-pressure prediction strategies need to be applied in the prodelta shales and inner-shelf deltaic sequences. Sonic and density log data detect overpressures generated by disequilibrium compaction, and pore pressures are accurately predicted using an Eaton exponent of 3.0. Sonic log data detect vertically transferred overpressures even in the absence of a porosity anomaly, and pore pressures are reasonably predicted using an Eaton exponent of 6.5.

Present-day stresses, seismicity and Neogene-to-Recent tectonics of Australia's 'passive' margins: intraplate deformation controlled by plate boundary forces

Abstract Neogene-to-Recent deformation is widespread on and adjacent to Australias ‘passive’ margins. Elevated historical seismic activity and relatively high levels of Neogene-to-Recent tectonic activity are recognized in the Flinders and Mount Lofty Ranges, the SE Australian Passive Margin, SW Western Australia and the North West Shelf. In all cases the orientation of palaeostresses inferred from Neogene-to-Recent structures is consistent with independent determinations of the orientation of the present-day stress field. Present-day stress orientations (and neotectonic palaeostress trends) vary across the Australian continent. Plate-scale stress modelling that incorporates the complex nature of the convergent plate boundary of the Indo-Australian Plate (with segments of continent–continent collision, continent–arc collision and subduction) indicates that present-day stress orientations in the Australian continent are consistent with a first-order control by plate-boundary forces. The consistency between the present-day, plate-boundary-sourced stress orientations and the record of deformation deduced from neotectonic structures implicates plate boundary forces in the ongoing intraplate deformation of the Australian continent. Deformation rates inferred from seismicity and neotectonics (as high as 10−16 s−1) are faster than seismic strain rates in many other ‘stable’ intraplate regions, suggestive of unusually high stress levels imposed on the Australian intraplate environment from plate boundary interactions many thousands of kilometres distant. The spatial overlap of neotectonic structures and zones of concentrated historical seismicity with ancient fault zones and/or regions of enhanced crustal heat flow indicates that patterns of active deformation in Australia are in part, governed, by prior tectonic structuring and are also related to structural and thermal weakening of continental crust. Neogene-to-Recent intraplate deformation within the Australian continent has had profound and under-recognized effects on hydrocarbon occurrence, both by amplifying some hydrocarbon-hosting structures and by inducing leakage from pre-existing traps due to fault reactivation or tilting.

Coupled changes in pore pressure and stress in oil fields and sedimentary basins

Repeated pressure measurements undertaken throughout the depletion of oil fields demonstrate that reduction in pore pressure is accompanied by a reduction in total minimum horizontal stress (σh), a phenomenon described herein as oil field-scale pore pressure/stress (Pp/σh) coupling. Virgin pressure measurements (i.e. those unaffected by depletion) through normally and overpressured sequences in sedimentary basins demonstrate that overpressure development is accompanied by an increase in σh, described herein as sedimentary basin-scale Pp/σh coupling. With depletion of the Ekofisk Field, North Sea, minimum horizontal stress decreased at approximately 80% of the rate of reduction of reservoir pore pressure (i.e. Δσh/ΔPp≈0.8). Virgin pressures measured in exploration wells surrounding the Ekofisk Field (Norwegian quadrants 1 and 2) indicate that with overpressure development Δσh/δPp≈0.73 (assuming shallow, normally pressured sequences are representative of overpressured sequences prior to overpressure development). Hence, despite the different temporal and spatial scales, the rate of decrease of minimum horizontal stress with pore pressure due to depletion of the Ekofisk Field is similar to the rate of increase of minimum horizontal stress with pore pressure due to overpressure development in the surrounding region. Basin-scale exploration pressure data in the Ekofisk region may thus provide an indication of the reservoir stress changes associated with depletion. Knowledge of such stress changes is critical because they can lead to the collapse of uncased wellbores, sand production and to faulting/fracturing and seismicity with field development.

Cenozoic exhumation of the southern British Isles

Rocks that crop out across southern Britain were exhumed from depths of as much as 2.5 km during Cenozoic time. This has been widely attributed to Paleocene regional uplift resulting from igneous underplating related to the Iceland mantle plume. Our compilation of paleothermal and compaction data reveals spatial and temporal patterns of exhumation showing little correspondence with the postulated influence of underplating, instead being dominated by kilometer-scale variations across Cenozoic compressional structures, which in several basins are demonstrably of Neogene age. We propose that crustal compression, due to plate boundary forces transmitted into the plate interior, was the major cause of Cenozoic uplift in southern Britain, witnessing a high strength crust in western Europe.

Present-day stress orientation in Brunei: a snapshot of ‘prograding tectonics’ in a Tertiary delta

The Baram Delta province of NW Borneo is unusual when compared with most other Tertiary deltas, as it has built up upon an active margin. Hence, structures observed in the Baram Delta province are the result of both margin-parallel gravity-driven deltaic tectonics and approximately margin-normal transpressive tectonics associated with the active margin. Image and dipmeter logs have been examined for breakouts and drilling-induced tensile fractures (DITFs) in 47 wells throughout Brunei. Breakouts and DITFs observed in 19 wells suggest that the maximum horizontal stress is oriented margin-normal (NW–SE) in the proximal parts of the basin and margin-parallel (NE–SW) in the outer shelf region. The margin-parallel outer shelf stress field is interpreted as a local ‘deltaic’ stress field caused by the shape of the clastic wedge. The margin-normal maximum horizontal stress in the inner shelf is interpreted to reflect basement stresses associated with the active margin. However, the maximum horizontal stress in the inner shelf is approximately perpendicular to the strike of Miocene–Pliocene normal growth faults, suggesting that maximum horizontal stress in the inner shelf has rotated from margin-parallel (‘deltaic’) to margin-normal (‘basement-associated’) over time. Hence, approximately the same stress rotation has occurred over time in the inner shelf as is currently observed spatially from the outer to inner shelf. The spatial and temporal stress rotations in Brunei are thus interpreted to be the result of ‘deltaic’ and ‘basement-associated’ tectonic regimes that are ‘prograding’ basin-wards. The proximity of the active margin has resulted in progressive uplift and inversion of the hinterland that has ‘forced’ the delta system to prograde rapidly. The zone of active deltaic growth faulting (and margin-parallel maximum horizontal stress) has shifted basin-wards (‘prograded’) as the delta system has rapidly prograded across the shelf. After uplift and delta progradation, the old growth faults of the inner shelf ceased being active and have then been successively reactivated by a similarly ‘prograding’ margin-normal inversion front.

Present-day stress and neotectonics of Brunei: Implications for petroleum exploration and production

The present-day state of stress in Tertiary deltas is poorly understood but vital for a range of applications such as wellbore stability and fracture stimulation. The Tertiary Baram Delta province, Brunei, exhibits a range of contemporary stress values that reflect the competing influence of the northwest Borneo active margin (situated underneath the basin) and local stresses generated within the delta. Vertical stress (v) gradients at 1500-m (4921-ft) depth range from 18.3 MPa/km (0.81 psi/ft) at the shelf edge to 24.3 MPa/km (1.07 psi/ft) in the hinterland, indicating a range in the shallow bulk density across the delta of 2.07–2.48 g/cm3. The maximum horizontal stress (Hmax) orientation rotates from margin parallel (northeast–southwest; deltaic) in the outer shelf to margin normal (northwest–southeast; basement associated) in the inner shelf. Minimum horizontal stress (hmin) gradients in normally pressured sequences range from 13.8 to 17.0 MPa/km (0.61–0.75 psi/ft) with higher gradients observed in older parts of the basin. The variation in contemporary stress across the basin reveals a delta system that is inverting and self-cannibalizing as the delta system rapidly progrades across the margin. The present-day stress in the delta system has implications for a range of exploration and production issues affecting Brunei. Underbalanced wells are more stable if deviated toward the hmin direction, whereas fracture stimulation in mature fields and tight reservoirs can be more easily conducted in wells deviated toward Hmax. Finally, faults near the shelf edge are optimally oriented for reactivation, and hence exploration targets in this region are at a high risk of fault seal breach.

Quantification of Tertiary Exhumation in the United Kingdom Southern North Sea Using Sonic Velocity Data

Sonic velocities from the Upper and Middle Chalk (Upper Cretaceous), the Bunter Sandstone, and the Bunter Shale (both Lower Triassic) were used to independently quantify apparent exhumation (height above maximum burial depth) in the United Kingdom (UK) southern North Sea. Apparent exhumation is the displacement, on the depth axis, of a given velocity/depth trend from the normal (unaffected by exhumation) trend. Apparent exhumation results derived from the Upper and Middle Chalk, the Bunter Sandstone, and the Bunter Shale are statistically similar. The consistency of results from carbonate and clastic units suggests that, at a formational and regional scale, overcompaction (i.e., anomalously high sonic velocity) in all three units analyzed reflects previously greater buria depth, rather than sedimentological and/or diagenetic processes, and validates the use of lithologies other than shale in maximum burial depth studies. The consistency of results from units of Early Triassic to Late Cretaceous age suggests that Tertiary exhumation was of sufficiently great magnitude to mask any earlier Mesozoic periods of exhumation, and that maximum Mesozoic-Cenozoic burial depth in the southern North Sea was attained prior to Tertiary exhumation. The proposed magnitudes of exhumation are generally greater than those previously published for the southern North Sea, but they are consistent with recent estimates from apatite fission track analysis. The amount of exhumation of the rock column is greatest on the recognized inversion axes (i.e., Sole Pit-Cleveland), where it reaches approximately 2.5 km. However, there was a regional component of 1.0-1.5 km exhumation during the Tertiary that affected structurally uninverted areas, and on which the more localized inversion-related exhumation was superimposed. Cretaceous-Tertiary burial prior to exhumation must have been of greater magnitude and more rapid than suggested by the preserved stratigraphy. The effect of this extra burial and subsequent exhumation on sedimentary rock decompaction procedure and thermal maturation modeling is illustrated for the Cleethorpes-1 and 44/7-1 wells, and must also be incorporated in modeling reservoir diagenesis. The regional, Tertiary tectonic uplift associated with exhumation must have had a thick-skinned origin.