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.

Joint Application of Diagenetic, Petrophysical and Geomechanical Data for Selecting Hydraulic Fracturing Candidate Zone

The more comprehensive information on the reservoir properties will help to better plan drilling and design production. Herein, diagenetic processes and geomechanical properties are notable parameters that determine reservoir quality. Recognizing the geomechanical properties of the reservoir as well as building a mechanical earth model play a strong role in the hydrocarbon reservoir life cycle and are key factors in analyzing wellbore instability, drilling operation optimization, and hydraulic fracturing designing operation. Therefore, the present study focuses on selecting the candidate zone for hydraulic fracturing through a novel approach that simultaneously considers the diagenetic, petrophysical, and geomechanical properties. The diagenetic processes were analyzed to determine the porosity types in the reservoir. After that, based on the laboratory test results for estimating reservoir petrophysical parameters, the zones with suitable reservoir properties were selected. Moreover, based on the reservoir geomechanical parameters and the constructed mechanical earth model, the best zones were selected for hydraulic fracturing operation in one of the Iranian fractured carbonate reservoirs. Finally, a new empirical equation for estimating pore pressure in nine zones of the studied well was developed. This equation provides a more precise estimation of stress profiles and thus leads to more accurate decision-making for candidate zone selection. Based on the results, vuggy porosity was the best porosity type, and zones C2, E2 and G2, having suitable values of porosity, permeability, and water saturation, showed good reservoir properties. Therefore, zone E2 and G2 were chosen as the candidate for hydraulic fracturing simulation based on their E (Young’s modulus) and ν (Poisson’s ratio) values. Based on the mechanical earth model and changes in the acoustic data versus depth, a new equation is introduced for calculating the pore pressure in the studied reservoir. According to the new equation, the dominant stress regime in the whole well, especially in the candidate zones, is SigHmax>SigV>Sighmin, while according to the pore pressure equation presented in the literature, the dominant stress regime in the studied well turns out to be SigHmax>Sighmin>SigV.

A Workflow to Derive Rock Mechanics Correlations and Stress Profile by the Integration of Dynamic and Static Formation Evaluation Data

Rock mechanics utilizes empirical formulas which are based on studies of certain environments. The shortcoming of such criteria is having estimations of rock physical properties with high uncertainty and not field/formation specific. The objective of this paper is to apply a core-log integration to convert dynamic mechanical properties captured from formation evaluation logs and calibrate them with core static data to generate a continuous profile of data with low uncertainty and generate correlations applicable to the specific physical environment. To obtain proper rock mechanical correlations, building a mechanical earth model (MEM) calibrated with core data and stimulation data is essential. Multiple wells drilled in a certain sandstone field with rock mechanical physical tests are analyzed. Multi-arm caliber data is also put in use to establish knowledge about in-situ stress directions. The procedure starts with gathering and filtering acoustic slowness & shear, formation pressure, density, and oriented multi-arm caliper logs. Next, calibration of dynamic to core static mechanical data collected in the lab is established. The geomechanical analysis includes an understanding of the state of stresses in a chosen reservoir along with rock elastic and failure properties. The complied data is then integrated using different workflows to develop Mechanical Earth Model (MEM). The intended rock mechanics correlations include elastic constants (Young's Modulus and Poisson's ratio), and rock failure parameters. Once Mechanical Earth Model (MEM) is established, dynamic logging data and core static data are correlated to produce key rock mechanics elements that are field and formation specific. The correlations include Young's Modulus, Poisson's Ratio, Unconfined Compressive Strength (UCS) correlation, and Friction Angle (FANG) correlation. A range of each rock mechanic element is also highlighted for the specific environment showcasing the limits expected for collapse and fracture. Ultimately, stress profile is generated with low uncertainty highlighting magnitudes of maximum and minimum horizontal stresses along with the given interval.

Perforation location optimization through 1-D mechanical earth model for high-pressure deep formations

Optimum perforation location selection is an important study to improve well production and hence in the reservoir development process, especially for unconventional high-pressure formations such as the formations under study. Reservoir geomechanics is one of the key factors to find optimal perforation location. This study aims to detect optimum perforation location by investigating the changes in geomechanical properties and wellbore stress for high-pressure formations and studying the difference in different stress type behaviors between normal and abnormal formations. The calculations are achieved by building one-dimensional mechanical earth model using the data of four deep abnormal wells located in Southern Iraqi oil fields. The magnitude of different stress types and geomechanical properties was estimated from well-log data using the Techlog software. The directions of the horizontal stresses are determined in the current wells utilizing image-log formation micro-imager (FMI) and caliper logs. The results in terms of rock mechanical properties showed a reduction in Poisson’s ratio, Young modulus, and bulk modulus near the high-pressure zones as compared to normal pressure zones because of the presence of anhydrite, salt cycles, and shales. Low maximum and minimum horizontal stress values are also observed in high-pressure zones as compared to normal pressure zones indicating the effects of geomechanical properties on horizontal stress estimation. Around the wellbore of the studied wells, formation breakouts are the most expected situation according to the results of the wellbore stress state (effective vertical stress (σzz) > effective tangential stress (σθθ) > effective radial stress (σrr)).

3D mechanical earth model for optimized wellbore stability, a case study from South of Iraq

Wellbore instability issues represent the most critical problems in Iraq Southern fields. These problems, such as hole collapse, tight hole and stuck pipe result in tremendous increasing in the nonproductive time (NPT) and well costs. The present study introduced a calibrated three-dimensional mechanical earth model (3DMEM) for the X-field in the South of Iraq. This post-drill model can be used to conduct a comprehensive geomechanical analysis of the trouble zones from Sadi Formation to Zubair Reservoir. A one-dimensional mechanical earth model (1DMEM) was constructed using Well logs, mechanical core tests, pressure measurements, drilling reports, and mud logs. Mohr–Coulomb and Mogi–Coulomb failure criteria determined the possibility of wellbore deformation. Then, the 1DMEMs were interpolated to construct a three-dimensional mechanical earth model (3DMEM). 3DMEM indicated relative heterogeneity in rock properties and field stresses between the southern and northern of the studied field. The shale intervals revealed prone to failure more than others, with a relatively high Poisson's ratio, low Young's modulus, low friction angle, and low rock strength. The best orientation for directional Wells is 140° clockwise from the North. Vertical and slightly inclined Wells (less than 40°) are more stable than the high angle directional Wells. This integration between 1 and 3DMEM enables anticipating the subsurface conditions for the proactive design and drilling of new Wells. However, the geomechanics investigations still have uncertainty due to unavailability of enough calibrating data, especially which related with maximum horizontal stresses magnitudes.

Geomechanics Insights for Successful Well Delivery in Complex Kutch - Saurashtra Offshore Region

Subsurface lithofacies sequences encountered in the Kutch & Saurashtra Basin has its own set of challenges brought about due to its complex geological settings. These challenges are related to drilling, logging and completion and demand rigorous planning for the upcoming wells with detailed analysis of hazards associated with the overburden and reservoir rocks. In the study, these challenges are found to be linked with three prime geological sequences. Detailed integrated geomechanical analysis with inputs from drilling parameters, real-time formation experience, geophysical and geological are conducted for the improvement in borehole condition and improvising the effective drilling rate. A customized geomechanical workflow has been adopted to construct Mechanical Earth Model (MEM, Plumb et al., 2000) for strategic wells across the basin. Wellbore stability events related to geomechanics were reproduced and analyzed. The cause of the events was established and mitigatory methods were proposed. In addition, stress orientation along the wellbore trajectory and across the basin was estimated using breakouts identified on images and multi-arm calipers. Fast shear azimuth from Dipole Shear Sonic anisotropy analysis was also integrated to provide more robust and accurate estimates. Wells in the region are characterized by slow ROP, high torque and drag, wellbore instabilities (severe held ups, cavings, stuck pipes, string stalling etc.) and challenges while logging and running casing. The study has characterized these challenges and identified required solutions linked to the three geological sequences - weak Tertiary, Late Cretaceous Deccan Trap and Early Cretaceous to Jurassic clastic formations. The Tertiary formations are relatively weak (UCS∼300 to 1500psi) and prone to sanding and cavings due to breakouts. MEM based mud weight window estimation predicts that shear/failure hole collapse can be prevented using 10ppg to 11ppg mud weight. The formations below the Deccan Trap are locally categorized under Mesozoic sequence. The Deccan Trap and Mesozoic formations are extremely hard, tight, extremely stressed, heavily fractured and in some areas are also of HPHT nature. Rock strength shows a wide variation (UCS ∼5,000psi to 25,000psi) making bit selection a difficult task. Borehole failure is complex and cuttings analysis shows the signature of both shear and weak plane failure. Fractures on the image logs, rotation of breakouts, and fast shear azimuth support this theory. Mixing fracture sealing agents along with the use of optimal mud weights is found to be the most likely drilling solution. The understanding developed in the region and implementation of recommended steps assisted in successful drilling of two recent wells wherein gun-barrel shape borehole condition in both Tertiary and the Mesozoic sequence was achieved. The non-productive time was reduced by nearly 40 days increasing the effective ROP by 40%. In addition, smooth borehole prevented any major issues while carrying out casing and cementing operations.

Building 1D Mechanical Earth Model for Zubair Oilfield in Iraq

Many problems were encountered during the drilling operations in Zubair oilfield. Stuckpipe, wellbore instability, breakouts and washouts, which increased the critical limits problems, were observed in many wells in this field, therefore an extra non-productive time added to the total drilling time, which will lead to an extra cost spent. A 1D Mechanical Earth Model (1D MEM) was built to suggest many solutions to such types of problems. An overpressured zone is noticed and an alternative mud weigh window is predicted depending on the results of the 1D MEM. Results of this study are diagnosed and wellbore instability problems are predicted in an efficient way using the 1D MEM. Suitable alternative solutions are presented ahead to the drilling process commences in the future operations.