Integration of Advanced Borehole Sonic and Resistivity Image Analysis for Fracture and Stress Characterisation - Implications to Carbon Sequestration Feasibility

Borehole resistivity images and dipole sonic data analysis helps a great deal to identify fractured zones and obtain reasonable estimates of the in-situ stress conditions of geologic formations. Especially when assessing geologic formations for carbon sequestration feasibility, borehole resistivity image and borehole sonic assisted analysis provides answers on presence of fractured zones and stress-state of these fractures. While in deeper formations open fractures would favour carbon storage, in shallower formations, on the other hand, storage integrity would be potentially compromised if these fractures get reactivated, thereby causing induced seismicity due to fluid injection. This paper discusses a methodology adopted to assess the carbon dioxide sequestration feasibility of a formation in the Newark Basin in the United States, using borehole resistivity image(FMI™ Schlumberger) and borehole sonic data (SonicScaner™ Schlumberger). The borehole image was interpreted for the presence of natural and drilling-induced fractures, and also to find the direction of the horizontal stress azimuth from the identified induced fractures. Cross-dipole sonic anisotropy analysis was done to evaluate the presence of intrinsic or stress-based anisotropy in the formation and also to obtain the horizontal stress azimuth. The open or closed nature of natural fractures was deduced from both FMI fracture filling electrical character and the Stoneley reflection wave attenuation from SonicScanner monopole low frequency waveform. The magnitudes of the maximum and minimum horizontal stresses obtained from a 1-Dimensional Mechanical Earth Model were calibrated with stress magnitudes derived from the ‘Integrated Stress Analysis’ approach which takes into account the shear wave radial variation profiles in zones with visible crossover indications of dipole flexural waves. This was followed by a fracture stability analysis in order to identify critically stressed fractures. The borehole resistivity image analysis revealed the presence of abundant natural fractures and microfaults throughout the interval which was also supported by the considerable sonic slowness anisotropy present in those intervals. Stoneley reflected wave attenuation confirmed the openness of some natural fractures identified in the resistivity image. The strike of the natural fractures and microfaults showed an almost NE-SW trend, albeit with considerable variability. The azimuth of maximum horizontal stress obtained in intervals with crossover of dipole flexural waves was also found to be NE-SW in the middle part of the interval, thus coinciding with the overall trend of natural fractures. This might indicate that the stresses in those intervals are also driven by the natural fracture network. However, towards the bottom of the interval, especially from 1255ft-1380ft, where there were indications of drilling induced fractures but no stress-based sonic anisotropy, it was found that that maximum horizontal stress azimuth rotated almost about 30 degrees in orientation to an ESE-WNW trend. The stress magnitudes obtained from the 1D-Mechanical Earth Model and Integrated Stress Analysis approach point to a normal fault stress regime in that interval. The fracture stability analysis indicated some critically stressed open fractures and microfaults, mostly towards the lower intervals of the well section. These critically stressed open fractures and microfaults present at these comparatively shallower depths of the basin point to risks associated with carbon dioxide(CO2) leakage and also to induced seismicity that might result from the injection of CO2 anywhere in or immediately below this interval.

Identifying Producing Horizons in Fractured Reservoirs a Far Field Acoustic Approach

The present study attempts to use 3D slowness time coherence (STC) technique to characterize the far-field fractures based on the reflector locations and attributes such as the dip and azimuth of fractures. These, in integration with the rest of the available data are used to accurately characterize the producing horizons in fractured basement reservoirs. The first step of the workflow involves the generation of 2D image to see if there are evidences of near and far wellbore reflectors. Since this is subjective in nature and does not directly provide quantitative results for discrete reflections, a new automated sonic imaging technique – 3D slowness time coherence (STC), has been incorporated to address this challenge. This method complements the image by providing the dip and azimuth for each event. The 2D and 3D maps of the reflectors can be readily available to integrate with the interpretations provided by other measurements, to better correlate and map the producing horizons. A field example is presented from the western offshore, India in which a fractured basement reservoir was examined using 3D STC technique to provide insight to the near and far field fracture network around the borehole. Few of the interpreted fractures from the resistivity image and conventional sonic fracture analysis coincide with the far field 3D STC reflectors, indicated by significant acoustic impedance. Further, the zones where the near and far field events coincide, represent a producing horizon. Comparing the near wellbore structures from the borehole images with the reflectors identified through the far field sonic imaging workflow provides necessary information to confirm the structural setting and characteristics of fractures away from the borehole. For the present case, it indicates the continuity of the fracture network away from the wellbore and explains the possibility of high production from the reservoir horizon. This study opens new perspective for identifying and evaluating fractured basement reservoirs using the sonic imaging technique. As more wells are drilled, it will be possible to better correlate and map the producing horizons in the field. This will allow better planning of location of future wells and help in optimizing field economics. A robust, automated and synergistic approach is used to locate and characterize individual arrival events which allows a more reliable understanding of the fracture extent and geologic structures. The 2D and 3D visualizations/maps can be readily integrated with the interpretations provided by other measurements.

Anisotropy analysis for optimized gas production in Coalbed Methane reservoir, Bokaro field, India

The efficient production of Coalbed Methane (CBM) gas is facing challenges due to the larger dewatering period from fracture connectivity to the aquifer zone. Also, commingled production from well makes it more difficult to identify the coal seam-wise problem. Therefore, prior knowledge of sub-surface fractures in coal seams is necessary to execute an accurate simulation model for planning hydraulic fracturing treatment. This paper highlights the studies in Bokaro CBM reservoir to mitigate challenges in few wells by characterizing anisotropy, determining fast shear wave polarization angle, maximum horizontal stress direction, fracture orientation, and analysis of low resistivity signature. Both the fast shear wave polarization angle and fracture orientation in resistivity image are observed in the same direction (N26°-35°E) in coal. The fast and slow shear slowness versus frequency plot concludes stress-induced anisotropy resulting from fractures that are supported by resistivity image and drilling core. Processing of the resistivity image log shows the maximum horizontal stress is along NE-SW direction, as identified from drilling-induced fractures. The observation of low resistivity signature with resistivity ranging from 0.4 to 0.8 ohm-m in few wells confirms the presence of conducting minerals such as siderite and pyrite from the x-ray diffraction studies of sidewall core. The present work guides in making production, drilling, and hydraulic fracturing design strategies to better understand the fluid propagation for optimized CBM production and will also help in future geomechanical studies.

APPLICATION OF IMAGE LOGS FOR ENHANCED RESISTIVITY-BASED WATER SATURATION ASSESSMENT IN ORGANIC-RICH MUDROCKS

Organic-rich mudrocks are complex in terms of rock fabric (i.e., the spatial distribution of rock components), which impacts electrical resistivity measurements and, therefore, estimates of hydrocarbon reserves. Conventional resistivity-saturation-porosity methods for assessment of water/hydrocarbon saturation do not reliably incorporate the spatial distribution of rock components and pores in the assessment of fluid saturation. Extensive calibration efforts are required for indirectly projecting the impact of rock fabric on resistivity models. For instance, none of the existing shaly-sand models incorporate a realistic distribution of clay network. This might be acceptable in conventional reservoirs. However, oversimplifying assumptions can cause significant uncertainty in reserves evaluation in organic-rich mudrocks. It should be noted that even the methods which incorporate the realistic distribution of rock components are difficult to calibrate. To address the aforementioned challenge, we introduce a joint interpretation of conventional resistivity and resistivity image logs to improve water saturation assessment by honoring the type of rock component, the spatial distribution of the conductive and non-conductive rock components, and the volumetric concentration of fluids and minerals in the rock. Borehole image logs are a source of high-resolution continuous rock sequence records and can provide detailed rock-fabric-related features. In this paper, we propose a method for the estimation of lamination density and mean resistivity value from image logs within each rock type. These fabric-related features are used to quantify the geometric model parameters for each conductive component of the rock. We use these geometric model parameters as inputs to a new resistivity model that considers volumetric concentration and spatial distribution of rock components for a depth-by-depth assessment of water saturation. The other inputs to the workflow are the volumetric concentration of conductive and non-conductive rock components, electrical conductivity of rock components, and porosity estimates from the joint interpretation of well logs. We successfully applied the proposed workflow to a dataset from the Wolfcamp formation in the Permian Basin in which resistivity image logs were available. We observed a measurable variation in estimated image-log-based geometric model parameters, which were in agreement with the visual content of the images. Incorporation of the estimated rock-class-based geometric model parameters in the resistivity model improved water saturation assessment. Results demonstrated a relative improvement in water saturation estimates of 44.2% and 59.1% against Waxman-Smits and Archie's models, respectively. We then used the estimated geometric model parameters for each rock type for a depth-by-depth assessment of water saturation in one additional well without image logs. This led to a faster and more reliable assessment of water saturation within a certain distance from the well with image logs, where the rock types remain comparable. This distance can be evaluated using variogram analysis. We demonstrated that using the estimated geometric model parameters could improve estimates of hydrocarbon reserves in the Permian Basin by approximately 34%. It should be noted that the proposed method for assessment of geometric model parameters is completely based on the actual spatial distribution of rock components and does not require core-based calibration efforts.

Magma Reservoir Beneath Azumayama Volcano, NE Japan as Inferred from Three-dimensional Electrical Resistivity Image by Magnetotellurics

An electrical resistivity image beneath Azumayama Volcano, NE Japan is modeled using magnetotellurics to probe the magma/hydrothermal fluid distribution. The 3-D inversion modeling images the conductive magma reservoir beneath Oana crater at depths of 3–15 km. The resolution scale for the conductor is 5 km by checkerboard resolution tests and the 67 % and 90 % confidential intervals of resistivity are 0.2–5 Ωm and 0.02–70 Ωm, respectively, for the region of less than 3 Ωm resistivity. The shallower part of the conductor is not explained by a water-saturated (5.5 wt %) dacitic melt, and the more probable interpretation is that it consists of a water-saturated, dacitic melt-silicic rock-hydrothermal fluid complex. The deeper part of the conductor is interpreted as a water-saturated (8 wt %) andesitic melt-mafic rock complex. The Mogi inflation source determined from GNSS and tilt data is located near the top boundary of the conductor at a depth of 2.7–3.7 km, which suggests that the ascent of hydrothermal fluids exsolved from the dacitic melt is interrupted by the impermeable wall and conduit. Assuming two phases of hydrothermal fluid and silicic rock, the resistivity at the inflation source, regarded as the upper bound resistivity of the conductor, is realized by the hydrothermal fluid fraction below the percolation threshold porosity in an effusive eruption. This indicates that the percolation threshold porosity in an effusive eruption characterizes the impermeable wall and conduit associated with the Mogi inflation source.

Understanding of subsurface conditions controlling flow liquefaction occurrence during the 2018 Palu earthquake based on resistivity profiles

The 7.4 Mw earthquake on 28th September 2018 in Palu City triggered a flow liquefaction phenomenon in the Balaroa and Petobo areas, contributing to significant casualties and building damage. This paper presents the results of a liquefaction study to map subsurface conditions in these areas using the multi-electrode resistivity method with the dipole-dipole configuration. The objective of this study is to understand factors controlling the flow liquefaction phenomenon. Based on the interpretation of 2-D resistivity images, the liquefied soil layers are characterized by lower resistivity values than the non-liquified layers. These contrasts of resistivity values form a gently sloping boundary between the liquefied and non-liquefied soil layers. The resistivity image perpendicular to the flow direction indicates the presence of a subsurface basinal morphology in the Balaroa area, suggesting that a shallow groundwater zone is present within the liquefiable soil layer. Thus, the subsurface topographical condition is the main governing factor of flow liquefaction phenomena during the 2018 Palu earthquake.