David Kulikowski

Geomechanical modelling of fault reactivation in the Cooper Basin, Australia

ABSTRACT The Australian Cooper Basin is a structurally complex intra-cratonic basin with large unconventional hydrocarbon potential. Fracture stimulation treatments are used extensively in this basin to improve the economic feasibility; however, such treatments may induce fault activity and risk the integrity of hydrocarbon accumulations. Fault reactivation may not only encourage tertiary fluid migration but also decrease porosity through cataclasis and potentially compartmentalise the reservoir. Relatively new depth-converted three-dimensional seismic surveys covering the Dullingari and Swan Lake 3D seismic surveys were structurally interpreted and geomechanically modelled to constrain the slip tendency, dilation tendency and fracture stability of faults under the present-day stress. A field-scale pore pressure study found a maximum pressure gradient of 11.31 kPa/m within the Dullingari 3D seismic survey, and 11.14 kPa/m within the Swan Lake 3D seismic survey. The present-day stress tensor was taken from previously published work, and combined with local pore pressure gradients and depth-converted field-scale fault geometries, to conclude that SE–NW-striking strike-slip faults are optimally oriented to reactivate and dilate. High-angle faults striking approximately E–W appear most likely to dilate, and act as fluid conduits irrespective of being modelled under a strike-slip or compressional stress regime. Near-vertical SE–NW and NE–SW-striking faults were modelled to be preferentially oriented to slip and reactivate under a strike-slip stress regime. Considering that SE–NW-striking strike-slip faults have only recently been interpreted in the literature, it is possible that many reservoir simulations and development plans have overlooked or underestimated the effect that fault reactivation may have on reservoir properties. Future work investigating the likelihood that fracture stimulation treatments may be interacting, and reactivating, pre-existing faults and fractures would benefit field development programs utilising high-pressure hydraulic fracture stimulation treatments.

Combining geophysical data and calcite twin stress inversion to refine the tectonic history of subsurface and offshore provinces: A case study on the Cooper‐Eromanga Basin, Australia

The use of core samples, borehole image logs, and seismic data is common practice for obtaining valuable structural data; however, these data are often obtained in isolation from other methods and not usually used for paleostress inversion processes. Therefore, for the first time, we present a new approach for constraining paleoprincipal stress orientations and regimes by integrating geophysical data (seismic and well data) with calcite twin principal stress orientation inversion analysis to refine the evolution of entirely subsurface or offshore basins; a case study on the subsurface Cooper-Eromanga Basin, Australia. Calcite samples were collected from oriented core, natural fracture data collected from borehole image logs, and fault data interpreted from three-dimensional seismic surveys. The analysis of microscale, mesoscale and macroscale data constrained the paleoprincipal stress orientations and regimes of six successive tectonic events: (1) NNW-SSE oriented strike-slip Carboniferous Alice Springs event; (2) SE-NW oriented compressional Mid-Permian event; (3) NE-SW oriented strike-slip Late Permian Daralingie event; (4) E-W compressional Late Triassic Hunter-Bowen event; (5) E-W compressional Late Cretaceous event; and (6) N-S compressional Paleogene event. This study shows the applicability of integrating geophysics with calcite twin stress inversion to decipher the tectonic evolution of entirely subsurface and offshore provinces.

Basement structural architecture and hydrocarbon conduit potential of polygonal faults in the Cooper-Eromanga Basin, Australia

ABSTRACT A thorough and complete understanding of the structural geology and evolution of the Cooper‐Eromanga Basin has been hampered by low‐resolution seismic data that becomes particularly difficult to interpret below the thick Permian coal measures. As a result, researchers are tentative to interpret the basement fault architecture within the basin, which is largely undefined. To provide a better understanding of the basement fault geometry, all available two‐dimensional seismic lines together with 12 three‐dimensional seismic surveys were structurally interpreted with assistance from seismic attribute analysis. The Upper Cretaceous Cadna‐owie Formation and top Permian reflectors were analysed using a common seismic attribute technique (incoherency) that was used to infer the presence of faults that may have otherwise been overlooked. Detailed basement fault maps for each seismic survey were constructed and used in conjunction with two‐dimensional seismic data interpretation to produce a regional basement fault map. Large north‐northeast–south‐southwest‐striking sinistral strike–slip faults were identified within the Patchawarra Trough appearing to splay from the main northeast–southwest‐striking ridge. These sinistral north‐northeast–south‐southwest‐striking faults, together with field‐scale southeast–northwest‐striking dextral strike–slip faults, are optimally oriented to have potentially developed as a conjugated fault set under a south‐southeast–north‐northwest‐oriented strike–slip stress regime. Geomechanical modelling for a regionally extensive system of Cretaceous polygonal faults was performed to calculate the Leakage Factor and Dilation Tendency of individual faults. Faults that extend into Lower Cretaceous oil‐rich reservoirs with strikes of between 060°N and 140°N and a high to near‐vertical dip angle were identified to most likely be acting as conduits for the tertiary migration of hydrocarbons from known Lower Cretaceous hydrocarbon reservoirs into shallow Cretaceous sediments. This research provides valuable information on the regional basement fault architecture and a more detailed exploration target for the Cooper‐Eromanga Basin, which were previously not available in literature.

Development of a new approach for hydraulic fracturing in tight sand with pre-existing natural fractures

Hydraulic fracturing in tight gas reservoirs has been performed in the Cooper Basin for decades in reservoirs containing high stress and pre-existing natural fractures, especially near faults. The hydraulic fracture is affected by factors such as tortuosity, high entry pressures, and the rock fabric including natural fractures. These factors cause fracture plane rotation and complexities, leading to fracture disconnection or reduced proppant placement during the treatment. In this paper, rock properties are estimated for a targeted formation using well logs to create a geomechanical model. Natural fracture and stress azimuths within the interval were interpreted from borehole image logs. The image log interpretations inferred that fissures are oriented 30–60° relative to the maximum horizontal stress. Next, diagnostic fracture injection test (DFIT) data was used with the poro-elastic stress equations to predict tectonic strains. Finally, the geomechanical model was history-matched with a planar 3D hydraulic fracturing simulator, and gave more insight into fracture propagation in an environment of pre-existing natural fractures. The natural fracture azimuths and calibrated geomechanical model are input into a framework to evaluate varying scenarios that might result based on a vertical or inclined well design. A well design is proposed based on the natural fracture orientation relative to the hydraulic fracture that minimises complexity to optimise proppant placement. In addition, further models and diagnostics are proposed to aid predicting the hydraulically induced fracture geometry, its impact on gas production, and optimising wellbore trajectory to positively interact with pre-existing natural fractures.

Composition changes of hydrocarbons during secondary petroleum migration

We investigate secondary migration of hydrocarbons with significant composition difference between the source and oil pools in the Cooper-Eromanga Basin, Australia. The secondary migration period is significantly shorter than the time of the hydrocarbon pulse generation, so neither adsorption nor dispersion of components can explain the concentration difference. The filtration coefficients, obtained from oil compositions in source rock (Patchawarra Formation) and in the reservoir (Poolowanna Formation and Hutton Sandstone), monotonically increase as carbon number increases. The monotonicity takes place for heavy hydrocarbons (n > 10). Loss of monotonicity for light and intermediate hydrocarbons can be explained by their evaporation into the gas phase. The evaporation of light and intermediate hydrocarbons into the gas phase is supported by their concentrations in oil, which are higher in source rock than in trapped reservoir oil. The paper proposes deep bed filtration of hydrocarbons with component kinetic retention by the rock. Introduction of the component capture rate into the mass balance transport equation allows matching the concentration difference, and the tuned filtration coefficients are in the common range. The results suggest that deep bed filtration controls the final reservoir oil composition during secondary migration in the Cooper-Eromanga Basin petroleum system, which was not previously considered.

Determining paleo-structural environments through natural fracture and calcite twin analyses: a case study in the Otway Basin, Australia

The structural history of the Otway Basin has been widely studied; however, previous works have focussed on large kilometre scale, basin and seismic structures, or have over-simplified natural fracture analysis with an excessive focus on fracture strike direction and a disregard for 3D geometry, a crucial characteristic when considering states of stress responsible for natural fracture formation. In this paper, we combine techniques of natural fracture analysis and calcite twin stress inversion to investigate the meso (outcrop and borehole) and micro (crystal) scale evidence for structural environments that have contributed to basin evolution. Our results indicate that basin evolution during the post-Albian may be markedly more complex than the previously thought stages of late Cretaceous inversion, renewed rifting and long-lived mid-Eocene to recent compression, with evidence for up to six structural environments detected across the basin, including; NE–SW and NW–SE extension, NW–SE compression, a previously undetected regime of NE–SW compression, and two regimes of strike-slip activity. By constraining structural environments on the meso- and micro-scale we can deliver higher levels of detail into structural evolution, which in turn, provides better-quality insights into multiple petroleum system elements, including secondary migration pathways and trap formation. Our research also shows that the Otway Basin presents a suitable environment for additional micro-scale structural investigations through calcite twin analyses.

An automated cross-validation method to assess seismic time-to-depth conversion accuracy: a case study on the Cooper and Eromanga basins, Australia

Selecting a seismic time-to-depth conversion method can be a subjective choice that is made by geophysicists, and is particularly difficult if the accuracy of these methods is unknown. This study presents an automated statistical approach for assessing seismic time-to-depth conversion accuracy by integrating the cross-validation method with four commonly used seismic time-to-depth conversion methods. To showcase this automated approach, we use a regional dataset from the Cooper and Eromanga basins, Australia, consisting of 13 three-dimensional (3D) seismic surveys, 73 two-way-time surface grids and 729 wells. Approximately 10,000 error values (predicted depth vs. measured well depth) and associated variables were calculated. The average velocity method was the most accurate overall (7.6 m mean error); however, the most accurate method and the expected error changed by several metres depending on the combination and value of the most significant variables. Cluster analysis tested the significance of the associated variables to find that the seismic survey location (potentially related to local geology (i.e. sedimentology, structural geology, cementation, pore pressure, etc.), processing workflow, or seismic vintage), formation (potentially associated with reduced signal-to-noise with increasing depth or the changes in lithology), distance to the nearest well control, and the spatial location of the predicted well relative to the existing well data envelope had the largest impact on accuracy. Importantly, the effect of these significant variables on accuracy were found to be more important than choosing between the four methods, highlighting the importance of better understanding seismic time-to-depth conversions, which can be achieved by applying this automated cross-validation method.

4D modelling of fault reactivation using complete paleostress tensors from the Cooper–Eromanga Basin, Australia

ABSTRACT Determining fault activity through time has typically utilised high-resolution seismic data to identify stratigraphic thickness changes or displacement vs distance plots; however, this approach is not possible in regions with low-resolution seismic data. We present a new approach for determining fault reactivation (tensile and shear) through time by integrating three-dimensional seismic data, geomechanical modelling and complete paleostress tensors from calcite twin stress inversion. The Cooper–Eromanga Basin is used as a case study to model the stress conditions present during six tectonic events that have affected the basin and, in doing so, constrain the effective paleostress magnitudes through time. Results show that the likelihood of dilation and shear reactivation of individual fault sets varies through time, with N–S- and E–W-striking faults likely to have been open to fluid flow after the critical moment in the hydrocarbon system. These results have substantial implications for hydrocarbon migration pathway models and structural and stratigraphic models for the Cooper–Eromanga Basin. This approach would benefit other provinces with low-resolution seismic data preventing fault growth analysis, or in regions where hydrocarbon migration pathways are poorly defined.