3.2.2. Drilling Events Review
Drilling Events Review is a process that involves a comprehensive evaluation of all well drilling-related events [
16]. A correct understanding of the formation’s pore pressure is essential not only for the safe and economical drilling of wells but also for assessing exploration risk factors, such as fluid migration and sediment integrity [
17]. This process encompasses the analysis of all events, from preparation to well completion, including safety, technical, and environmental incidents, to enhance safety, efficiency, and overall outcomes of well drilling operations while minimizing potential risks.
Overburden Stress. The overburden stress
is the pressure exerted on a formation at a specific depth due to the total weight of rocks and fluids above that depth [
18]. Overburden stress is calculated by integrating the total formation density from the surface to the total depth (TD), following the equation provided by Oloruntobi [
19]. Total density can also be obtained from wireline or LWD logs or from core density measurements. In this case, measured density data is available for most intervals in wells S-3 and S-4.
Figure 6 illustrates the calculation of overburden stress in two wells, S-3 and S-4. The analysis presents several curves representing essential information related to these wells’ density, vertical pressure, and overburden stress gradients. In the figure, the red curve on track 4 represents extrapolated density, while the curve on track 5 shows the vertical pressure calculated in psi. Track 6 displays the overburden stress gradient in ppg.
Observations from
Figure 6 indicate that the overburden stress across the research area exhibits good consistency and trend. The similarity in the shape and pattern of the curves in both wells, particularly the vertical pressure curve (track 5) showing a linear increase with depth, indicates that overburden stress increases proportionally with depth [
20]. The overburden stress gradient (track 6) shows relatively constant values throughout the area, indicating a uniform stress distribution. The consistency and trend observed in the overburden stress data confirm the accuracy and reliability of the data obtained in this study, which is crucial for ensuring valid geological and geotechnical interpretations.
Pore Pressure. Pore pressure refers to the fluid pressure generated within the pore spaces of a formation. It is estimated using the Eaton method, with compressional velocity as the primary input [
21]. Two trendlines are applied to account for the higher organic content in the Talang Akar and PRT formations. The overpressure zone extends through the Telisa, BRF Equivalent, Pendopo, and Talang Akar formations. Additionally, three pressure inversions are identified in the DTCO model. Pore pressure in the PRT formation is assumed to follow a gas gradient of 0.1 psi/ft. This result is then validated using pressure measurements from well S-3 and a DST result of 3794 psi at a depth of 2490 mMD in well S-4. The pore pressure models for wells S-3 and S-4 are illustrated in
Figure 7 and
Figure 8.
Rock Parameters. Elasticity parameters of rocks such as Young’s modulus, Poisson’s ratio, shear modulus, and bulk modulus represent the elastic deformation behavior of rocks [
22] These parameters serve as inputs for calculating rock strength and in-situ stress. Dynamic elastic properties are obtained from compressional sonic logs, shear sonic logs, and bulk density logs, which are constructed using equations following the method by Bozorgi [
23].
For the geomechanical modeling of this field, static elastic properties measured from rock mechanics tests are required. In this well analysis, rock mechanics test data is unavailable, so typical values from the nearest field have been utilized. The static Young’s modulus is estimated using the Horsrud correlation [
24], leveraging input from compressional velocity and assuming that the static Poisson’s ratio is 0.9 of the dynamic Poisson’s ratio.
Unconfined Compressive Strength (UCS) is used to determine rock strength parameters. UCS calculations for the Laman Lului III to LTAF formations utilize the Horsrud equation with compressional velocity as input, while the Coates-Denoo correlation uses static Young’s modulus, bulk modulus, and shale volume as inputs to calculate UCS in the PRT formation. The friction angle is calculated using the Lal equation with computed velocity input. Meanwhile, the formation’s tensile strength (TSTR) is used to evaluate borehole tensile failure due to stress relaxation. Tensile strength is typically low, ranging from 1/12 to 1/8 of UCS, and in this case, it is estimated to be 1/10 of UCS by Al-Zubaidy and Al-Jawad [
25].
The graphs displayed in
Figure 9 and
Figure 10 illustrate elastic properties and rock strength in wells S-3 and S-4, respectively. The graphs show, in sequence, the measured depth, formation tops, inclination & azimuth, mechanical stratigraphy, bit size and caliper, gamma ray, compressional and shear velocity, bulk density and neutron porosity, static Young’s modulus & Poisson’s ratio, TSTR and UCS, friction angle, well scheme, and true vertical depth.
Horizontal Stress Profile. No existing research exists on the orientation of horizontal stresses for the S field or the nearby area. Therefore, the direction of horizontal stresses is assumed based on the regional stress map [
26], as shown in
Figure 11. The direction of maximum horizontal stress is NE-SW, approximately 135 degrees NW-SE.
With the values of vertical stress, minimum horizontal stress, and maximum horizontal stress, the stress regime of the formation can be identified [
10]. In this context, the study employs a poroelastic model to determine horizontal stress, incorporating the effects of geothermal and tectonic strains into the calculations. (Marelis et al. [
27]). By neglecting thermal or temperature effects, the poroelastic model can be simplified into an equation that encompasses both minimum and maximum horizontal stresses.
This approach provides a mathematical framework that enables the calculation of horizontal stress magnitudes without the need for direct measurements, accounting for the strain effects from geothermal sources and tectonic activity [
28].
WBS Analysis. Wellbore stability analysis was conducted to validate the MEM. The MEM is a numerical representation of the state of stress and rock mechanical properties for a specific stratigraphic section in a field or basin [
29]. The shear failure model for breakout (Track 11) was validated using caliper data (borehole enlargement) (Track 12). The Breakout model from the S-3 well, shown in
Figure 12, demonstrates good consistency with the measured caliper log and drilling events related to wellbore instability (narrow hole). Similarly, the breakout model from the S-4 well, shown in
Figure 13, exhibits good consistency with the measured caliper log and drilling events related to wellbore instability (point narrow), implying that the breakout model can accurately predict the location and severity of breakouts occurring in the wellbore [
30]. This alignment indicates that the breakout model can be effectively used to understand and analyze wellbore. Although wellbore degradation can be observed in image and caliper logs, damage to the wellbore wall may go undetected during the drilling process if the damaged material remains in place [
31].
The MEM from S-3 is then projected onto S-8, while the MEM from S-4 is propagated to S-9. This projection helps understand and anticipate potential wellbore instability at the location based on previous experience and relevant data.
In
Figure 14, the results of the wellbore prediction model for the curtain section of the S-8 well are shown. This section is also referred to as the Curtain section. It should be noted that no log data was directly collected from the S-8 well and its lateral tracks. Therefore, log information from the 1D MEM, including DTCO, DTSH, and RHOB, along with elastic properties and rock strength derived from the S-3 well, was applied to S-8.
This process involves propagating log data from S-3 to S-8 to perform WBS calculations at the S-8 location. In other words, information obtained from the S-3 well is used to model the wellbore conditions and characteristics in S-8 [
32]. This step is crucial due to the lack of direct data from S-8, and by leveraging available data from the nearest well, the prediction model can provide a more accurate representation of the geological and mechanical conditions at the desired location.
Furthermore, in
Figure 15, the limitations of calibration data for S-8 are shown. All parameters are kept consistent across all support wells, and the breakout model is adjusted based on caliper data and/or drilling events. The rock strength in the PRT formation in S-8 is reduced to model the presence of natural fractures. This is reflected in the decrease of UCS to 30% of its original value, with an average value of 3708 psi. The TSTR is taken as 0.1 of the reduced UCS, and the shear angle is reduced to 60% of its original value, with an average value of 24.6 degrees. The breakout model from S-8 shows good consistency with drilling events related to WBS, including a pack-off at a depth of 2138 mMD. The validated parameters for rock strength reduction are adopted for S-9 as a low-case scenario for WBS analysis.
WBS is an essential factor influencing the safety and efficiency of operations. (Aslannezhar et al. [
33]). This analysis underscores the importance of understanding the rock and geological characteristics of the well to prevent potential failures and accidents during the drilling process. The simulated reduction in rock strength, accounting for natural fractures in the PRT formation, emphasizes the need to understand the local geological environment (Zhao et al. [
34]). Moreover, the consistency of the breakout model with drilling events in S-8 demonstrates the model’s accuracy of the model used in predicting well conditions. Using validated parameters in WBS analysis for subsequent wells, such as S-9, reflects a systematic and evidence-based approach to managing risks and improving operational efficiency.
Figure 16 depicts the Curtain section and wellbore prediction model of S-9, which geometrically evolves from S-4 due to its proximity to the planned well. This modeling is conducted by considering a set of 1D MEM logs, including DTCO, DTSH, RHOB, as well as elastic and rock strength parameters. Additionally, in the development of S-9, pore pressure, horizontal stress, and pressure load are recalculated.
The analysis of the selection of S-4 as a reference point for S-9 underscores a strategic approach to new well placement. By selecting a nearby location, it is anticipated that the geological and rock characteristics in the vicinity will exhibit significant similarities, enabling more accurate propagation of MEM parameters. The use of 1D MEM logs, including DTCO, DTSH, and RHOB, provides detailed insights into rock properties and elasticity, while rock strength parameters are essential for assessing formation stability.
The recalculation of pore pressure, horizontal stress, and pressure load in S-9 highlights the importance of adjusting critical parameters to account for the well’s unique geological environment [
35]. These recalculations provide a more accurate representation of subsurface conditions, aiding in the identification of potential risks and ensuring well stability during the drilling and production phases.
Figure 17 provides an overview of the wellbore stability analysis and mud weight recommendations for well S-9. The pore pressure estimation in S-9 accounts for a pressure drop of 1600 psi in the PRT formation. Two scenarios were considered in the wellbore stability analysis: the first scenario uses the original rock strength from S-4 for the entire interval, while the second scenario incorporates a reduction in rock strength in the PRT formation to capture the potential for well collapse, as observed in S-8, while other formations retain the rock strength properties of S-4.
One of the drilling challenges arises from suboptimal mud weight planning, which can lead to borehole instability, characterized by the formation of a cave along the borehole wall [
36]. Mud weight (MW) recommendations are determined by considering the planned well cases and ensuring a low risk of shear failure. These recommendations are taken above the general trend of the 0% Depth of Damage, although in some layers, the recommended mud weight is lower than the shear failure gradient. This reflects a cautious approach to maintaining well stability without compromising safety. A sensitivity analysis is conducted to evaluate whether changes in inclination and azimuth can expand the safe MW window [
37]. This analysis highlights the critical role of pore pressure estimation, proper rock strength selection, and MW recommendations in maintaining the stability of well S-9. The two-case strategy in sensitivity analysis demonstrates a prudent approach to managing the potential risk of well collapse.
Sensitivity Analysis. The results of the study in
Figure 18 demonstrate that sensitivity analysis plays a critical role in evaluating wellbore stability under various formation and operational conditions. In the MW formation, the deviation has a significant impact, with a high sensitivity to increasing deviation. As deviation increases, the MW window narrows due to a rise in the shear failure gradient. Recommendations for the PRT formation suggest drilling with a maximum deviation of 30 degrees. As above this value, there is no safe MW window in Case 2. Therefore, maintaining deviation below 30 degrees ensures a safe MW window for both cases.
Sensitivity analysis also highlights significant risk factors related to rock strength and stress values. In the PRT formation, wellbore stability analysis reveals a high risk, particularly in the use of Case 2. Practical recommendations include drilling with mud loss to prevent pack-off or wellbore collapse. The importance of accurately monitoring cuttings and signs of collapse is emphasized, with the recommendation to consider a gradual increase in MW if collapse occurs. To improve the accuracy of wellbore stability analysis, obtaining complete logging data, including shear velocity measurements using sonic technology, is advised. This step helps address uncertainties in shear velocity measurements. The results of the wellbore stability analysis are also presented in Schmidt contour plots, showing critical MW or ECD limits for breakout and breakdown at various azimuths and deviations. These sensitivity plots offer valuable visual information for designing MW and determining casing points for future well trajectories. In the MW analysis, the left stereo net provides a planned view as the “target board” for stability values, highlighting the sensitivity of MW breakout to well deviation and azimuth. Emphasis on breakout risk and crack initiation is represented in the MW window plot, which accounts for deviation and azimuth. Yellow and dark blue shading indicates the minimum and maximum MW required to prevent breakout and crack initiation. The white area between breakout and mud loss represents a stable MW window relative to deviation and azimuth, providing critical guidance for well-operational planning.