Multi‐scale stratigraphic forward modelling of the Surat Basin for geological storage of CO2

Underground geological storage of CO₂ (GSC) requires a high level of subsurface understanding that is often hindered by a lack of data. This study demonstrates the use of stratigraphic forward modelling (SFM) in generating and characterising a static reservoir model using limited well data, with multiple potential applications within the GSC workflow. Sedsim SFM software was used to create a static model of the Surat Basin, including a high‐resolution nested model of the EPQ‐7 GSC tenement within the basin. Deposition and burial of the Jurassic Precipice Sandstone, Evergreen Formation and Hutton Sandstone were simulated. Modelling results show a close match with gamma‐ray well logs in the tenement area, and the model can be considered a credible model of the subsurface. The Sedsim‐predicted formation thicknesses and porosity and permeability distributions meet criteria set for GSC, suggesting that the EPQ‐7 tenement may be a prospective GSC location.     

Multigrain sedimentation/erosion model based on cross-shore equilibrium sediment distribution: Application to nourishment design

In the light of global warming and sea level rise there are many coastal beaches that suffer from erosion. Beach nourishment has become a common practice to maintain the sediment balance on a shore-face. In this paper, a three-dimensional numerical model for evaluating long-term impact of beach nourishment projects has been developed. The model addresses the longstanding complex issue of coastal morphology and sediment grain size distribution from an unconventional angle, which exploits the strong links between grain size distribution and the prevailing transport direction of each sediment constituent under ‘average’ wave and storm action. The present model predicts the redistribution of nourished sediment according to the subtle clues implied by equilibrium distribution curves and latest coastal wave transformation theories. After verification against recent field observations in Terschelling, The Netherlands, the model was used to predict long-term effects of different beach nourishment strategies. It was found that: (a) given the source sediment available in Terschelling the tactics of large volume and less frequent implementation are better than otherwise; and (b) from a pure engineering point of view, waterline nourishment outperforms offshore trough nourishment.The model offers an additional tool for coastal engineers to evaluate the feasibility, effectiveness and the optimization of dumping locations for beach nourishment projects. It is also a useful tool for stratigraphic modelling of shallow-marine sedimentation in conjunction with sea level changes.     

Method for the determination of residual carbon dioxide saturation using reactive ester tracers

The mechanisms for storage of CO₂ in rock formations include structural/stratigraphic, mineral, solubility and residual trapping. Residual trapping is very important in terms of both containment security and storage capacity. However, to date, the contribution from residual trapping (i.e. immobilisation of supercritical fluid via capillarity in pore spaces) is still relatively difficult to quantify accurately. Using a laboratory-based testing program, this study demonstrates the feasibility of using reactive ester tracers (i.e. triacetin, propylene glycol diacetate and tripropionin), which partition between a mobile water phase and a stationary supercritical CO₂ phase, to quantify the residual CO₂ saturation, S_gr, of a rock formation. The proposed single-well test involves injecting these tracers into the subsurface, followed by CO₂ saturated water, where the ester tracers slowly hydrolyse to form products with differing partition coefficients. After a suitable period of time, allowing for partial hydrolysis, water containing the tracer mixture is produced from the subsurface and analysed using gas chromatography mass spectrometry (GCMS). A numerical simulator of the tracer behaviour in a reservoir is used to explain the differential breakthrough of these tracer compounds during water production to estimate S_gr. Computer modelling suggests that the use of esters tracers to determine CO₂ residual saturation is a potentially robust method. The supercritical CO₂/water partition coefficients directly dictate the amount of time that each tracer spends in the CO₂ and water phases. As such for modelling of tracer behaviour and estimating S_gr, knowing the tracer partition coefficient is essential; in this paper, the first laboratory study to determine the partition coefficients of these reactive ester tracers is described.     

A numerical simulation of supersonic thermal convection

A computational scheme has been constructed for solving the equations that describe strong thermal convection in a two-dimensional gas layer that is heated from below and is stratified across many scaleheights by a uniform gravitational field. The purpose of this scheme is to mimic the physical conditions that may have existed in a section of the proto-solar cloud from which the planetary system formed. The vertical temperature gradient of the initial quiescent layer of diatomic gas is strongly superadiabatic and matches that of a polytrope of index m  = 1. The temperature at the upper boundary is kept fixed during the computation. Because of the highly compressible nature of the gas and the steep spatial gradients, a modified version of a flux-corrected transport scheme due to Zalesak is devised. The computations show that after the convection adopts a steady-state configuration, the flow consists of horizontal pairs of giant convective cells of opposing circulation. At the cell boundaries, the downflows are rapid and spatially concentrated while the upflows are broad and sluggish. Supersonic speeds are easily achieved in the downflows. Contrary to the expectations of the mixing length theory of convection, there is a net downward flux of kinetic energy at each level in the layer. The convecting layer is cooler on average compared with the initial temperature profile, and there is a net shift of mass towards the lower boundary. The implications of these results for the modern Laplacian theory of Solar system origin are briefly discussed.     

Automating conceptual models to easily assess trap integrity and oil preservation risks associated with fault reactivation

Extensive Neogene fault reactivation and leakage in the Timor Sea, Australian North-West Shelf, has long been identified as the likely cause for traps in the region being dry or underfilled, despite widespread evidence for hydrocarbon charge and larger palaeo-hydrocarbon columns than those preserved today. A structural model has been previously proposed (Gartrell et al., 2005) to explain the distribution of palaeo- and live hydrocarbon columns. This basic geometrical model is used to define the volume of closure that is protected from fault reactivation, thereby providing the potential for preserving hydrocarbons. Although the model is retrospectively successful in identifying which traps contain hydrocarbons, application to other areas is subjective. To assess the application of the Gartrell et al. (2005) concept a computational model has been generated and applied to a regional 3250 km² seismic interpretation, allowing sensitivity testing of the critical input parameters.The automation has allowed us to demonstrate that the model successfully identifies dry traps and derives reasonable approximations of hydrocarbon-water contacts in commercial accumulations. Therefore, the basic geometrical model is deemed valid to apply to other hydrocarbon provinces, both frontier and mature, that are adversely impacted by fault reactivation after hydrocarbon charge.