AADE-10-DF-HO-15 Implementing and Monitoring a Geomechanical Model for Wellbore Stability in an Area of High Geological Complexity J. G. Vargas, Halliburton; J. Mateus, Halliburton; A. Jaramillo, Halliburton; M. Serrano, Halliburton; M. Rivera, Ecopetrol; O. Mercado, Ecopetrol. Copyright 2010, AADE This paper was prepared for presentation at the 2010 AADE Fluids Conference and Exhibition held at the Hilton Houston North, Houston, Texas, April 6-7, 2010. This conference was sponsored by the Houston Chapter of the American Association of Drilling Engineers. The information presented in this paper does not reflect any position, claim or endorsement made or implied by the American Association of Drilling Engineers, their officers or members. Questions concerning the content of this paper should be directed to the individuals listed as authors of this work. Abstract The drilling of an oil well in an area of geological complexity presents many challenges in the planning and drilling phases. In the case presented here from Colombia, the complexities included a highly-deviated well profile, dipping formations, and naturally fractured formations. Moreover, the wellbore challenges were compounded by the high stresses in a tectonically-active area. Hence, to help ensure success of the operation, it was necessary to incorporate new technologies to define the risk that the drilling case presented and to prepare for the potential problems predicted in the pre-well planning. For this project, the best trajectory to drill the well from a geomechanical and operational point of view could not be used because of the formation dips. A different azimuth was used, which required a new wellbore trajectory. With the new wellbore profile, the safe drilling window for wellbore stability was redefined. During the successful drilling of the difficult interval, the wellbore stability issues were monitored closely with frequent updates in the modeling. From the various issues described in this paper, the reader can better understand the necessary planning and wellbore monitoring that is required for drilling difficult formations in tectonically-active areas. Introduction In order to increase well production, the operator company selected an alternative trajectory to drill the wells in the area, positioning the well at the highest point of the structure and drilling in the down-dip direction, thus allowing this well to cross a thicker drainage area. This highly deviated wellpath would drill through an interbedded formation of shale and naturally fractured limestone lithologies. The axis would be almost parallel to the dip formation and main faulting plane. These aspects reflect the complex drilling environment of the project. Hence, it was crucial to apply geomechanical concepts to decrease the risk associated with drilling the well. Geomechanical modeling for wellbore stability can address the different conditions that constitute this type of well. However, in order to quantify the effect of the variables mentioned above (drilling parallel to bedding planes in a shale formation and drilling through a naturally fractured formation), it is necessary to have fracture porosity and permeability data, fluid viscosity inside the fractures, poroelastic parameters like undrained Poisson, Biot coeffient to model the natural fracture effects, and parameters like cohesion and the friction angle of the bedding planes and the bedding plane orientation that account for the angle of attack analysis (angle between the well and the bedding plane). Most of this information was not available, so the wellbore stability model did not include the effect of the natural fractures and angle of attack in a quantitative way but in a qualitative way. Therefore, the collapse and fracture pressures estimated during the analysis were not influenced by these aspects, but the interpretation of the results and the recommendations about mud weight window, drilling fluid aspects related to wellbore stability and some operational aspects took into account the effect of the issues mentioned. The results of the geomechanical and wellbore stability modeling and how the various aspects were used are discussed below, including the following up during drilling. Project Development The geomechanical contributions to the project are divided into three principal phases: 1. Pre-drilling where the geomechanical model (pore pressure, mechanical properties and stresses) and wellbore stability model (collapse and fracture pressures) were initially defined based on the off-set wells. 2. While-drilling where a real time following up was done at the rig site. 3. Post-drilling where with electrical logs and the events documented during while-drilling stage the final wellbore stability model was defined. This paper addresses the Pre-drilling and While Drilling phases. It is important to mention that during most of the planning and execution phases, interdisciplinary work was performed by reservoir and operation personnel from the operator company, as well as directional drilling and geomechanics personnel from the services company. This synergy allowed the project to benefit from dynamic planning even during the execution.
2 J. G. Vargas, J. Mateus, A. Jaramillo, M. Serrano, M. Rivera, O. Mercado AADE-10-DF-HO-15 Pre-drilling Phase The field is located in the upper Magdalena valley (Colombia) where only one offset wells (PPX-1) had the complete information needed to make the geomechanical model. Information from the other eight wells was used to calibrate the models and to define the principal risks associated with drilling operations in the field. Information like electrical logs, image logs, rock mechanic tests, final drilling reports, leak off and formation integrity tests, mini frac tests, petrophysical model defined by the operator company, injection data and production data were used to perform the analysis. In Figure 1 a seismic section of the area shows the locations of the prospect well (PPX-2) and PPX-1 well. The differences are clear as to how the wells cross the structure with the subsequent uncertainty for the data extrapolation from well to well. Figure 1 also shows the relationship between the well trajectory and the bedding plane dips as well as the faulting zone interpreted by the operator using seismic data. A geomechanical model is defined by pore pressure, mechanical properties and the in situ stress state (including stress orientation). With the geomechanical model as input data, a wellbore stability model can be estimated by calculating the formation s collapse and fracture pressures. This methodology was applied to the PPX-1 well and the wellbore stability model was compared with drilling data to validate that the model could reproduce by itself what was observed during drilling. The model obtained for PPX-1 well confirmed most of the events reported. Once this validation was done, the geomechanical model was extrapolated to the PPX-2 well and evaluated with the prognosis formation tops, and trajectory. To perform a geomechanical evaluation, the basic log suite includes sonic, density and a lithology log such as gamma ray. The logs available for the project are shown in Figure 2; all these logs are from PPX-1 well, which had more logs of better quantity and quality. However, as can be seen in Figure 2, there is an interval without density information. To complete the density log, the Gardner et al [1] (1974) relation was used and the result is shown in Figure 3 where the synthetic density log is shown in green color. Once the electric logs were complete, the geomechanical modeling began with mechanical properties calculations, followed by pore pressure estimations and ending with stresses magnitude and orientation. There are different methods for obtaining the mechanical rock properties to be used in a geomechanical modeling; in this study two different sources were used. One of them was calculating the mechanical properties from different mathematical correlations from electric logs. The profiles obtained with the methodology mentioned before were calibrated by using results from rock mechanics laboratory tests performed on a core from the reservoir in PPX-1. Table 1 presents the lab results. Figure 4 presents the mechanical properties results estimated with PPX-1 well logs and the comparison with the laboratory results. To estimate the pore pressure values, a relation between compressive sonic log and depth was defined by using normal compaction trend lines [2]. The pore pressure profile was defined by correlating the trend line and the sonic log measures. The pore pressure profile was compared and calibrated with mud weight used and events during drilling operations. Figure 5 shows the final pore pressure profile used during the pre-drilling phase. Knowledge of stress magnitudes and orientations at great depth is of appreciable interest in both the geologic sciences and engineering. One of the most important uses of in-situ stress data in the petroleum industry is associated with problems of wellbore Stability [3]. The stress magnitudes were estimated using some correlations as a function of some mechanical properties and pore pressure. The magnitude of this theoretical horizontal stress needs to be calibrated with field data to include the effects of lateral tectonic strains and thermal rock deformations [4]. The predominant stress regime observed according with the results is strike slip (according to Anderson s classification [5]) as can be seen in Figure 6. The orientation of breakouts corresponds to the direction of the minimum horizontal stress whereas the orientation of induced fractures corresponds to the direction of the maximum horizontal stress in vertical wells [4]. The orientation of the principal horizontal stresses was calculated using image logs, where some open natural fractures and hydraulic fractures were identified. The average orientation of the maximum horizontal stress is E-W. Figure 7 shows diagrams where the orientation of the fractures can been seen. Generally, a wellbore fails either by exceeding the tensile strength or the shear strength of the formation. In addition, due to the laminated texture of shale, wellbore stability may be governed by failure of bedding planes [6]. Once the geomechanical model for the PPX-1 well was defined, the wellbore stability model was calculated using an elastic analysis with Mohr-Coulomb failure criteria. As a result, collapse and fracture pressures were defined and used to determine the possible mechanical status and the mud weight window. The trajectory sensibility of the model was determined with polar charts where the required mud weight for a specific depth is a function of wellbore azimuth and inclination. As an example, Figures 8-9 demonstrate the effect of the stress magnitudes and mechanical properties on the most stable trajectory at different depths, even with the same horizontal stress orientation. According to multiple polar charts created during analysis, the optimum well orientation should be NE- SW. However, it is important to remember that variables like attack angle could not be introduced in the quantitative analysis, so the polar charts did not reveal the effect of bedding dip. Figure 10 shows the wellbore stability model for the PPX- 1 well. The model was calibrated using the final borehole shape according to the PPX-1 caliper log. Some examples of
AADE-10-DF-HO-15 Implementing and Monitoring a Geomechanical Model for Wellbore Stability in an Area of High Geological Complexity 3 the comparison between the model results and caliper log can be observed in Figures 11-12. Since, a good correlation between wellbore stability model and the caliper log was observed, the geomechanical model was extrapolated to prognosis formation tops and well trajectory for the PPX-2 well. The wellbore stability model for the PPX-2 well is shown in Figure 13 with the safe mud weight window proposed. With the observations made from geomechanical and wellbore stability modeling, some risks and mitigation plans were proposed including maintaining concentrations of a sealing, bridging agent and shale stabilizers in the drilling fluid. This would help control stability related to natural fracture because friction angle reduction of fractures due to mud infiltration significantly affects wellbore stability during drilling [7]. Once this is accomplished, the mud weight can be increased to control shear failure in the borehole. A particle size distribution (PSD) analysis was recommended to determine the correct concentrations of the drilling fluid additives. While-drilling Phase Due to some operational issues, a sidetrack was performed in PPX-2 well. Due to this side track, it was necessary to calibrate and monitor wellbore stability model in real time. The main activities during while-drilling phase included caving monitoring (volume, shape and sizes), review of pressure-while-drilling (PWD) measurements and whiledrilling events related to stability and hole cleaning. On the other hand, the actualization of the model with observations made at well site and according to the new well trajectory was done. The results during follow up showed that the previous wellbore stability modeling was very close to the mechanical behavior of the well. Figure 14 shows the previous modeling (collapse and fracture pressures) and the mud weight used in the original well before the first sidetrack (ST1). This figure shows how the model predicted zones where caving rates increase (cyan line) and it became necessary to increase the mud weight. It is clear that as the collapse pressure increases, the caving volume increases too. The arrows in Figure 14 indicate the behavior of the caving rate and the correspondence between it and the collapse pressure predicted in the planning phase. Taking into account that the predicted potentiality wellbore stability tendencies of the formations were observed by means of different factors like caving rates, increases in annular pressures and difficult trips, and that modeling does not include critical factors like the angle of attack and natural fractures influences, it is considered that the modeling represents the geomechanical response of the formations drilled. The fracture pressure calibration can be done by means of several sources: step rate test (SRT), extended leak-off test (XLOT) or hydraulic fracturing. The formation integrity test (FIT) indicates only the minimum value that the fracture gradient may be. Due to the fact that there was no more information, the FIT values were used in the way mentioned. Once the adjustment of the modeling in the original hole was completed, the data extrapolation was done according to the tops and trajectory of the ST1. Figure 15 shows the wellbore stability model after the extrapolation to the ST1. The monitoring mentioned previously was the base line for the changes made to the mud weight window to drill the ST1. The main change was to increase the mud weight according to the collapse pressure estimated for ST1. As a result of the mentioned changes, a lower caving production rate was seen as shown in Figure 16. Figure 17 shows the predominant morphology registered. Figure 18 shows an image log acquired in the 11133 ft to 11772 ft (MD) interval. In this image a difference of less than 1.5 between bit size and borehole diameter can be observed. This is evidence of the stability of the interval as predicted by the wellbore stability model. The follow up of the model at the rig site explains this widening as a combined effect of tabular caving production and drill string eccentricity due to the high inclination of the hole. Several more similar examples were done during project development. After about 43 days of drilling the ST1, the well reached total depth achieving the principal objectives and indicating that these types of well geometries can be drilled with a daily planning work during execution, and the use of tools during the project development that can help drilling engineers reduce risk by identifying it. Conclusions 1. The estimated mechanical properties from logs match with the test results from the laboratory. 2. The predominant stress regime in the area corresponds to a strike slip regime. 3. The most stable trajectory to drill wells in the area corresponds with E-W orientation. 4. The wellbore stability model developed during planning was accurate based on events observed during drilling, especially in the interval prior to the first sidetrack. 5. To handle instability from natural fractured zones and shear failure at the same time it is recommended to increase the concentration of sealing materials such as graphite and asphalts prior to increase mud weight. 6. The detailed caving monitoring during drilling is a very helpful tool for making decisions about increasing mud weight to manage wellbore instability. Acknowledgments The authors would like to thank both Halliburton and Ecopetrol for their cooperation and support in developing this project. Nomenclature PWD = Pressure While Drilling ST1 = Side Track 1 UCS = Unconfined Compressive Strength
4 J. G. Vargas, J. Mateus, A. Jaramillo, M. Serrano, M. Rivera, O. Mercado AADE-10-DF-HO-15 To = Tensile Strength P Velocity = Compressive wave velocity S Velocity = Shear wave velocity References 1. Gardner, G.H.F., L.W. Gardner, and A.R. Gregory. Formation Velocity and Density ; The Diagnostic for Stratigraphic traps. Geophysics, vol. 39, 770 780, 1974. 2. Mouchet, J.P. Mitchell, A. Abnormal Pressures While Drilling. Elf Aquitaine, Manuels Techniques 2. 1989. 3. Zoback, M.D et all. Determination of Stress Orientation and Magnitude in Deep Wells. International Journal fo Rock Mechanics and Miniing Sciences 40 (2003) 1049-1076. 4. Anderson, E.M. The Dynamics of faulting and dyke formation with applications to Britain. Edinburgh: Oliver and Boyd; 1951. 5. Franquet, J.A et all. Critically-Stressed Fracture Analysis Contibutes to Determining the Optimal Drilling Trajectory in Naturally Fractured Reservoirs. International PetroleumTechnology Conference, 2008. 6. Cheng, X. Tan, C.P. Haberfield, C.M. Guidelines for Efficient Wellbore Stability Analisis. Int. J. Rock Mechanics and Min Sci. Vol 34, No 3-4, 1997. 7. Cheng, X. Tan, C.P. Detournay, C. A Study on Wellbor Stability in Fractured Rock Masses With Impact of Mud Infiltration. Journal of Petroleum Science and Engineering 38. (2003) 145-154. Tables Table 1. Result for mechanical properties measured using a core samples from reservoir formation in PPX-1 well. SAMPLE DEPTH (ft) 12113.63 12115.17 12119 UCS (psi) 10134.24 11947.4 12488.39 To (psi) 572.28 N/D 1396.6 P Velocity (m/seg) S Velocity (m/seg) Min. Min. 2624.61 1361.35 2667.7 1647.63 2728.12 1696.071 Max. Max. 3969.92 2069.01 4185.16 2099.74 3394.22 2075.91 Average Poisson Dinamic 0.352 0.279 0.174 Relation Pseudo-static 0.313 0.252 0.155 Average Young Dinamic 2.826 3.598 3.213 Modulus (1E06 psi) Pseudo-static 2.669 2.301 2.08
Figures Principal offset well. Less complex area Offset well TD: 12640 (MD) Drilling parallel to faulting zone (too many open natural fractures) Down dip drilling direction Figure 1. Seismic section of the area under analysis. Figure 2. Logs inventory in PPX-1 well. An interval without density log is observed.
6 J. G. Vargas, J. Mateus, A. Jaramillo, M. Serrano, M. Rivera, O. Mercado AADE-10-DF-HO-15 Figure 3. Synthetic density log (green curve), modeled to complete the input data to develop the geomechanical analysis. Poisson Young UCS Tensil 11800 0.00 0.20 0.40 0.60 11800 1.50E+05 3.15E+06 6.15E+06 9.15E+06 11800 0 5000 10000 15000 20000 25000 30000 11800 0 800 1600 2400 3200 4000 4800 5600 11900 11900 11900 11900 12000 12000 12000 12000 Depth TVD Depth TVD Depth TVD Depth TVD 12100 12100 12100 12100 12200 12200 12200 12200 12300 12300 12300 12300 Figure 4. Comparison of the mechanical properties estimated during the analysis (yellow profiles) with lab test (red points).
AADE-10-DF-HO-15 Implementing and Monitoring a Geomechanical Model for Wellbore Stability in an Area of High Geological Complexity 7 Pore pressure PPX-1 well PP (ppg) 4 6 8 10 12 14 0 TD Seccion 26 2000 4000 TVD (ft) 6000 TD Seccion 17.5 8000 10000 TD Seccion 12.25 12000 TD Seccion 8.5 Figure 5. Pore pressure profile (blue line) estimated from sonic log in PPX-1 well. Red line corresponds to MW used during drilling. Horizontal dashed lines correspond to end of hole sections and horizontal continuous lines correspond to formation tops. 500 Stress Magnitude Gradient (psi/ft) 0.30 0.55 0.80 1.05 1.30 1500 2500 3500 4500 5500 Depth (ft) 6500 7500 8500 9500 10500 11500 Minimum horizontal stress Maximum Horizontal stress Vertical stress Minifract Test Figure 6. Stresses magnitudes estimated for the project. A predominant strike slip regime is observed.
8 J. G. Vargas, J. Mateus, A. Jaramillo, M. Serrano, M. Rivera, O. Mercado AADE-10-DF-HO-15 Open Natural Fractures Induced Fractures Figure 7. Open natural fractures and induced fractures orientation determined from image log from PPX-1 well. Model: Isotropic; Elastic; Impermeable; Vertical Stress = 2149.8 PSI (0. 940 PSI /feet) Max Hor Stress = 2492. 8 PSI (1.090 PSI/feet ) Min Hor Stress = 1646.6 PSI (0.720 PSI /feet) Pore Pressure = 960.5 PSI (0.420 PSI/feet) Critical Mudweight Polar Charts -- Shear Failure -- Collapse 330 0 Shmin 30 Distance into formation (r/r) = 1.05 True Vertical Dept h = 2287 feet Cohesion = 700.00 PSI; Friction Angle = 35.00 Failure Criterion = Mohr- Coulomb No BreakOut Angle Always Stable (M W < 0. 00) Hole Azimuth (deg) 300 60 Always Fail (MW > 20.983) (l b/g al ) 8. 605 8. 505 8. 405 8. 305 8. 205 Hole In cl inatio n Angle 270 SHmax 90 SHmax 8. 105 8. 006 7. 906 7. 806 7. 706 7. 606 240 120 7. 506 7. 406 N 7. 306 7. 206 210 150 W E 7. 106 180 S 7. 006 Shmin PBORE- 3D 7. 10, 2008 Figure 8. Polar chart at 2287ft TVD.
AADE-10-DF-HO-15 Implementing and Monitoring a Geomechanical Model for Wellbore Stability in an Area of High Geological Complexity 9 Model: Isotropic; Elastic; Impermeable; Vertical Stress = 4968.6 PSI (0. 980 PSI /feet) Max Hor Stress = 5424. 9 PSI (1.070 PSI/feet ) Min Hor Stress = 3650.4 PSI (0.720 PSI /feet) Pore Pressure = 2129.4 PSI (0. 420 PSI/ feet) Critical Mudweight Polar Charts -- Shear Failure -- Collapse 330 0 Shmin 30 Distance into formation (r/r) = 1.05 True Vertical Depth = 5070 feet Cohesion = 1171.00 PSI ; Friction Angle = 28.00 Failure Criterion = Mohr- Coulomb No BreakOut Angle Always Stable (MW < 0. 00) Hole Azimuth (deg) 300 60 Always F ail (MW > 20.598) (l b/g al ) 11.148 11.050 10.952 10.853 10.755 Hole In cl inatio n Angle 270 SHmax 90 SHmax 10.657 10.558 10.460 10.362 10.264 10.165 240 120 10.067 9. 969 N 9. 871 9. 772 210 150 W E 9. 674 180 S 9. 576 Shmin PBORE- 3D 7. 10, 2008 Figure 9. Polar chart at 5070ft TVD. Wellbore Stability Model for PPX-1 RÉGIMEN RUMBO-DESLIZANTE well. MW (ppg) 0 5 10 15 20 0 2000 4000 TVD (ft) 6000 8000 10000 12000 Collapse Mudw eight Fracture Mudw eight MW (ppg) Figure 10. Wellbore Stability Model estimated for PPX-1 well.
10 J. G. Vargas, J. Mateus, A. Jaramillo, M. Serrano, M. Rivera, O. Mercado AADE-10-DF-HO-15 RÉGIMEN RUMBO-DESLIZANTE Wellbore Stability Model MW (ppg) 3 8 13 1000 Bit Size - CAL (in) 10 15 20 1000 1200 1200 1400 1400 TVD (ft) TVD (ft) 1600 1600 1800 1800 2000 2000 Collapse Mudw eight Fracture Mudw eight MW (ppg) BS CALIPER Figure 11. Comparison between wellbore stability model and caliper log in PPX-1 well. Example from 1000 to 2000 ft TVD RÉGIMEN RUMBO-DESLIZANTE Wellbore Stability Model MW (ppg) 3 8 13 18 6100 Bit Size - CAL (in) 7 12 17 6100 TVD (ft) TVD (ft) 6400 6400 Collapse Mudw eight Fracture Mudw eight MW (ppg) BS CALIPER Figure 12. Comparison between wellbore stability model and caliper log in PPX-1 well. Example from 6100 to 6550 ft TVD
AADE-10-DF-HO-15 Implementing and Monitoring a Geomechanical Model for Wellbore Stability in an Area of High Geological Complexity 11 Wellbore Stability Window, PPX-2 Well PPG 0 5 10 15 20 25 1000 3000 5000 Fm. 1 Fm. 2 Profundidad (ft) 7000 Fm. 3 Fm. 4 9000 Fm. 5 Fm. 6 11000 Fm. 7 Collapse Fracture Pore Pressure Figure 13. Wellbore stability model for the prospect well (PPX-2), and mud weight window proposed (Purple lines). Wellbore Stability Model (following up) 25 5.8 20 4.8 15 3.8 PPG 2.8 10 5 1.8 0.8 0-0.2 6500 7000 7500 8000 8500 9000 9500 Fm. 3 Fm. 4 Fm. 5 Fm. 6 Profundidad TVD (ft) Collapse Pressure Fracture Pressure FIT MW Caving Figure 14. Comparison between wellbore stability model and caving rate during PPX2 well drilling. Observe how where MW is too close to collapse pressure the caving rate increase.
12 J. G. Vargas, J. Mateus, A. Jaramillo, M. Serrano, M. Rivera, O. Mercado AADE-10-DF-HO-15 25 Wellbore Stability Model (Updated) 20 15 PPG 10 5 0 5000 5500 6000 6500 7000 7500 8000 8500 9000 9500 10000 10500 11000 11500 12000 12500 Fm. 2 Fm. 3 Fm. 4 Fm. 5 Fm. 6 Fm. 7 Profundidad TVD (ft) Extrapolated Collapse Extrapolated Fracture FIT Planning Fracture Planning Collapse Figure 15. Comparison between wellbore stability model before updating (planning) and model after following up (extrapolated). Daily Caving Volume 10.0 9.0 8.0 7.0 Volume (bbl/hr) 6.0 5.0 4.0 This high caving volumes are not product of an increasing instability, this correspond to an enhance hole cleaning. 3.0 2.0 1.0 0.0 9300 9600 9900 10200 10500 10800 11100 11400 11700 12000 12300 12600 12900 13200 13500 13800 Depth MD (ft)
AADE-10-DF-HO-15 Implementing and Monitoring a Geomechanical Model for Wellbore Stability in an Area of High Geological Complexity 13 Figure 16. Caving volume during drilling ST1. Figure 17. Example of the predominant caving morphology produced during drilling ST1. It can be observed the tabular shape.
14 J. G. Vargas, J. Mateus, A. Jaramillo, M. Serrano, M. Rivera, O. Mercado AADE-10-DF-HO-15 Figure 18. Example of the image log took in PPX-2 well. No more than 10 diameter is observed in an interval drilled with 8.5 bit, in a formation drilled with high deviation angle and producing tabular cavings.