A Novel Workflow to Model Geomechanical Properties Using Fracture Coherence Interpretation and Density Inversion in the Rangal Coal Measures of the Bowen Basin in Australia

Abstract

Abstract The study focused on the Permian aged Rangal coal measures within a 3D seismic volume in middle of Bowen basin. The Rangal Coal Measures (RCM) are a middle-rank coal seam gas (CSG) target with low permeability. Approximately 150 to 200 wells, with geophysical log data, are unevenly distributed in the seismic area of 51 km2. The effect of stress and related geomechanics on producing wells has been poorly understood. Petrophysical data, 3D seismic cube, geographic information system (GIS), borehole imager and well stress data were used to conduct the interpretative workflow. A novel four-step workflow was created to provide geomechanical insights in order to optimize future well designs for improved productivity and well integrity. The first step in the workflow was to finely delineate fracture distribution by using coherence and max curvature methods in the seismic volume. This was followed by determining the most appropriate number of coal-sensitive attributes using step-wise regression methods. The multi-attribute emerged density inversion was conducted utilizing probabilistic neural network training and predicting, which was derived from the density cube and predicted coal thickness map. The resulting density cubes included the overburden zone which was calculated by well data interpolation in consideration of GIS overlays. This was followed by 2D stress simulation that was created by analyzing a combination of maximum/ minimum stress, principal stress, Young's Modulus, and Possion's ratio in measured well data. Inputs such as fracture networks, coal thickness, density and velocity cubes were combined to process the simulation using designed fracture reservoir characterization software (FRS). The workflow was completed by using GeoPressure Analysis (GPA) software and data such as: stress, density and the velocity cubes to calculate the overburden pressure (Pov) and pore pressure (Pp). The overburden pressure was computed by integrating the density volume. The pore pressure volume was analyzed using the Fillippone formula. Vertical effective stress (VES) and the pore pressure coefficient were calculated using the principles of rock theory. The minimum / maximum horizontal principal stress (SHmin / SHmax) were calculated by utilizing the Zoback and Healy formulas and using the following input data: calculated Pp, Pov, measured effective stress ratio of the minimum/maximum horizontal principal stress and measured horizontal stress. Direct measured properties from geophysical logs such as principal stress showed good correlation with adjacent extracted values from the modeled principal stress. The well production properties indicated good correlation with geomechanical properties such as small fractures and stress. Single well average peak gas production and average peak gas production both decreased with the diversity factor of horizontal stress (DHSR) which was extrapolated from SHmax and SHmin. The originality of this study lies in the comprehensive integration of geological, geophysical and petrophysical data to create a functional geomechanical model. The results indicate that this geomechanical study will be useful for well planning and predicting production properties in each future study area.