The deepwater shallow formation has low fracture pressure and narrow safe window of mud density, which results in a high risk of wellbore instability in this kind of formation.
Without considering the plastic state of shallow formation around the borehole or the influence of in-situ stress difference on wellbore stability, the prediction accuracy of the traditional wellbore stability analysis models are relatively low. This paper can provide a reliable method to accurately predict the safe window of drilling fluid density.
In this paper, the shallow formation around the borehole is divided into plastic zone and elastic zone considering it under non-uniform in-situ stress. The collapse pressure formula of shallow formation is derived by taking the shrinkage rate of the borehole as the instability criterion. The fracturing pressure calculation model of shallow sediment under non-uniform in-situ stress is derived by combining the theory of excess pore pressure and hydraulic fracturing.
The calculated results indicate that the horizontal in-situ stress difference has a significant effect on the shape of the plastic zone, the shrinkage rate of borehole, collapse pressure and fracturing pressure. The calculated results are in good agreement with the field test results, and the prediction accuracy of this model is higher than that of other traditional models.
Marine oil and gas resources are very rich in the world, which is much more difficult to develop than in onshore oil field [1J. Shelton, J.R. Smith, and A. Gupta, "Experimental evaluation of separation methods for a riser dilution approach to dual density drilling", Journal of Energy Resources and Technology, vol. 133, no. 3, pp. 611-620, 2011.
[http://dx.doi.org/10.1115/1.4004968] , 2B.H. Yu, C.L. Yan, J.G. Deng, S.J. Liu, Q. Tan, and K. Xiao, "Evaluation and application of wellbore stability in deep water", Oil Drilling & Production Technology, vol. 33, no. 6, pp. 1-4, 2011.]. Deepwater shallow sediment is less-consolidated and has high permeability and low strength. Because of its unique mechanical properties, the shallow formation around the borehole often has higher collapse pressure and lower fracture pressure. As a result, the narrow safe window of drilling fluid density usually causes the occurrence of drilling accidents [3J. Dodson, T. Dodson, and V. Schmidt, "Gulf of Mexico ‘trouble time’ creates major drilling expenses: use of cost-effective technologies needed", Offshore, vol. 64, no. 1, pp. 46-48, 2004.-5L.A.S. Rocha, P. Junqueira, and J.L. Roque, "Overcoming deep and ultra deepwater drilling challenges", In: Offshore Technology Conference, 2003, p. OTC-15233-MS.
[http://dx.doi.org/10.4043/15233-MS] ].
Many studies have been conducted on wellbore stability of deepwater shallow sediments, for example, Wojtanowicz et al. established an elastic model which can reflect the mechanical characteristics of deepwater shallow sediment, and studied its fracturing pressure based on both theoretical analysis and numerical simulation [6A. Paknejad, J. Schubert, and M. Amani, "A new method to evaluate leak-off tests in shallow marine sediments (SMS)", In: SPE Saudi Arabia Technical Symposium, 2007, p. SPE-110953-MS.
[http://dx.doi.org/10.2118/110953-MS] ]; R. Chhajlani and Z. Q. Zheng analyzed the stability of borehole wall in Mexico Bay deepwater oil field based on the linear elastic constitutive model and Mohr-Coulomb strength criterion and used the research results to guide the field work [7R. Chhajlani, and Z.Q. Zheng, "Utilization of geomechanics for Medusa field development, deepwater gulf of Mexico", SPE Annual Technical Conference and Exhibition, 2002 SPE-77779-MS.
[http://dx.doi.org/10.2118/77779-MS] ]. However, the traditional models of borehole wall stability analysis are based on elastic mechanics [8L. Cui, Y. Abousleiman, A.H.D. Cheng, and J.C. Roegiers, "Time dependent failure analysis of inclined boreholes in fluid-saturated formation", Journal of Energy Resources Technology, vol. 121, no. 1, pp. 31-39, 1999.
[http://dx.doi.org/10.1115/1.2795057] , 9A.M. Al-Ajmi, and M.H. Al-Harthy, "Probabilistic wellbore collapse analysis", Journal of Petroleum Science & Engineering, vol. 74, no. 3, pp. 171-177, 2010.
[http://dx.doi.org/10.1016/j.petrol.2010.09.001] ]. Since the shallow sediment is similar to saturated soil, even if the stress around the borehole exceeds the rock strength, the borehole will remain stable and not collapse or fracture immediately. Therefore, the traditional analytical model cannot reveal the real borehole instability mechanism of deepwater shallow sediment. D. G. Ritter and B. Grollimund studied borehole stability in Brazil deepwater oil field based on the method of quantitative risk assessment, and the results showed that the prediction accuracy of in-situ stress and formation strength had the greatest influence on the analysis results [10D.G. Ritter, and B. Grollimund, "Wellbore stability in the deepwater Bijupira & Salema fields, offshore Brazil-a probabilistic approach", In: Offshore Technology Conference, 2003, p. OTC-15205-MS.
[http://dx.doi.org/10.4043/15205-MS] ]. E. Kaarstad and B. S. Aadnoy established an elastoplastic model for the fracturing pressure calculation of shallow formation based on the field data of deepwater drilling in Norway. But the empirical model based on field practice cannot truly reflect the actual stress state around the borehole [11E. Kaarstad, and B.S. Aadnoy, "Fracture model for general off-shore applications", In: SPE Asia Pacific Oil & Gas Conference and Exhibition, 2006, p. SPE-101178-MS.
[http://dx.doi.org/10.2118/101178-MS] , 12E. Kaarstad, and B.S. Aadnoy, "Improved prediction of shallow sediment fracturing for offshore applications", Spe Drilling & Completion, vol. 23, no. 2, pp. 88-92, 2008.
[http://dx.doi.org/10.2118/101178-PA] ]. L. Yan et al. divided the shallow formation around the borehole into plastic zone and elastic zone and deduced the collapse pressure and fracturing pressure calculation formula based on the assumptions that the formation is in a uniform stress state and the plastic zone around the borehole is circular [13C.L. Yan, J.G. Deng, and X.R. Lai, "Borehole stability analysis in deepwater shallow sediments", Journal of Energy Resources Technology, vol. 137, pp. 1-6, 2015.
[http://dx.doi.org/10.1115/1.4027564] ]. However, many researches and field tests showed that the horizontal stress difference has a significant impact on the collapse pressure and fracturing pressure although the horizontal stress difference in deepwater shallow formation is smaller than that in terrestrial formation. Therefore, the prediction accuracy of the calculation model under the assumption of uniform in-situ stress is relatively low.
In this paper, considering the shallow formation is in the condition of non-uniform in-situ stress and combining the rock mechanics and some theories of soil mechanics, the shallow formation around the borehole is divided into plastic zone and elastic zone. The formulas for calculating the plastic zone boundary and the shrinkage rate of the borehole under non-uniform in-situ stress are deduced by assuming that the stress-strain relation in the plastic zone is in accordance with the theory of large deformation. The shrinkage rate is chosen as the criterion of borehole collapse, and the collapse pressure formula of shallow formation is obtained. Then the fracturing pressure calculation model of shallow sediment under non-uniform in-situ stress is derived by combining the theory of excess pore pressure and hydraulic fracturing theory.
Due to the lack of sufficient geological and logging data in deepwater drilling, the errors of calculation results of safety drilling fluid density window by the traditional models are always large. It is very important to establish a model to accurately predict the collapse pressure and fracture pressure of deepwater shallow formation. Although there are great differences in the physical and mechanical properties between the deepwater shallow sediment and the terrestrial formation, their research methods of borehole stability are same.
However, it is more difficult to analyze the wellbore stability in deepwater shallow formation. Firstly, overburden pressure of deepwater shallow sediment is lower than that of the terrestrial formation, which leads to a low fracturing pressure of shallow formation and a narrow safe window of drilling fluid density. Therefore, it is very important for deepwater drilling to improve the prediction accuracy of the safe window of drilling fluid density. Secondly, deepwater shallow formation is characterized by high porosity, high permeability, weak cementation and low strength, which is prone to transform into plastic state and appear as large deformation under the action of in-situ stress. In the process of drilling, some serious accidents often occur in shallow formation, such as wellbore shrinkage, drill pipe sticking and wellbore collapse. Thirdly, shallow formation may have shallow gas and shallow water, and well blowout or gas channeling is easy to occur. Because the shallow gas and shallow water can only be formed under special conditions, the first two problems are more important to evaluate the well bore stability in deepwater shallow formation.
Deepwater shallow sediment is similar to saturated soil and can be considered as homogeneous, isotropic and ideally elastoplastic material. The traditional elastic model is relatively conservative in calculating the collapse pressure and fracturing pressure of shallow sediment. Many researches have proved that even if there is a certain plastic area appeared around the borehole in the shallow formation, the borehole wall can still maintain stability. Considering the effect of non-uniform in-situ stress, the mechanical model is established, as shown in Fig. (1). The shallow formation around the borehole is divided into elastic zone (E) and plastic zone (P). When borehole pressure (P_{w}) approaches the lower limit of the safe window of drilling fluid density, the formation around the borehole is in a plastic state. With the increase of borehole pressure, the formation is in gradual transition to the elastic state. But when borehole pressure continues increasing, the formation will turn back into the plastic state again. In Fig. (1), r_{w0} is bit radius during the drilling, mm; r_{w} is the borehole radius after drilling, mm; r_{b} is the distance elastic-plastic interface to the borehole center, mm; θ is angle from the orientation of minimum horizontal principal stress.
Fig. (1) Borehole mechanics model of shallow sediment. |
Deepwater shallow sediment has obvious plastic characteristics. When a borehole is drilled into the deepwater shallow formation, the stresses are redistributed around the wellbore. If the stresses exceed a yield criterion, a plastic zone can appear around the wellbore, which means that the shallow formation around wellbore experiences certainly a plastic stage before instability. Many studies have shown that the plastic zone boundary is not a standard shape of circle due to the non-uniformity of in-situ stresses [14L.S. Zhang, X.Z. Yan, X.J. Yang, Z.L. Tian, and H.L. Yang, "An elastoplastic model of collapse pressure for deep coal seam drilling based on Hoek–Brown criterion related to drilling fluid loss to reservoir", Journal of Petroleum Science and Engineering, vol. 134, pp. 205-213, 2015.
[http://dx.doi.org/10.1016/j.petrol.2015.07.003] , 15Y. Jin, M. Chen, and C. Wu, "Borehole stability prediction before drilling for formations under the second section in an exploration well", Petroleum Exploration and Development, vol. 35, no. 6, pp. 742-745, 2008.]. Fig. (1) shows the wellbore stress elastoplastic analysis model, where the approximate boundaries of the elastic zone and plastic zone are depicted. Based on the linear elastic theory, the stress distribution in the elastic zone around wellbore for shallow formation can be determined as follows [16W.B. Bradley, "Mathematical concept—stress cloud can predict borehole failure", Oil & Gas Journal, vol. 77, no. 8, pp. 92-102, 1979.-18F.L. Ning, "Research on Wellbore Stability in Gas Hydrae Formation", Ph.D thesis, China University of Geosciences, Beijing, China, 2005.]:
(1) |
Where σ_{r} and σ_{θ} are the radial and tangential stress.
Q. H. Sun et al. believe that the deepwater shallow formation in elastic zone is the ideal elastic material and obeys the Mohr Coulomb strength criterion [19Q.H. Sun, J.G. Deng, C.L. Yan, B.H. Wei, S.J. Liu, Z.L. Liu, and K. Xiao, "Calculation method for fracture pressure of deep sea shallow formation", Journal of Central South University (Science and Technology), vol. 46, no. 4, pp. 1402-1408, 2015.]. It shows that Eq.(1) can be applied to the calculation of the stress in the elastic zone around wellbore for shallow formation. According to the study of B. H. Yu et al. [20B.H. Yu, C.L. Yan, J.G. Deng, J.L. Zhou, S.J. Liu, and X.J. Xing, "Study and application of calculation model of safe drilling fluid density window—A case study of AKPO Oil field, West Africa", China Offshore Oil and Gas, vol. 24, no. 2, pp. 58-60, 2012.], the shallow sediment around the borehole will turn into the plastic state after the formation of wellbore. There must be a transition point from the elastic zone to the plastic zone in the formation. Mohr-Coulomb criterion can be used to determine the boundary between the two zones [21B.S. Jiang, L. Yang, and L.P. Shi, "Stress analysis of cracked surrounding rock based on Hoek-Brown criterion", Chinese Journal of Solid Mechanics, vol. 32, pp. 300-305, 2011.]. According to the Mohr-Coulomb strength criterion, the yield function of deepwater shallow sediments is:
(2) |
Where σ_{θ} and σ_{r} represent the minimum and maximum principal stresses, respectively N = (1 + sinθ) / (1 - sinθ), σ_{ 0} = (2C cosθ) / (1- sinθ). C and φ are the cohesion and internal friction angle, respectively. σ_{0 } is the uniaxial compressive strength (UCS), MPa.
Combining Eqs. (1) and (2), the distance of elastic-plastic interface from the borehole center(r_{b}) can be written as:
(3) |
Where
Deepwater shallow sediment is similar to saturated soil and has strong plasticity. When stress around the borehole exceeds the rock strength, the formation will turn into plastic state. In this case, the borehole will remain stable and not collapse immediately. The collapse pressure calculation model based on the elastic is relatively conservative. Therefore, considering the effects of borehole pressure, the balance equation can be established as:
(4) |
According to the Mohr-Coulomb strength criterion, the stress field in the plastic zone can be obtained by combining Eqs. (2) and (4). When σ_{r}_{|}_{r = rw}=P_{w}, the radial stress and tangential stress are:
(5) |
(6) |
Substituting Eq. (3) into (5) and (6), the stress on the elastic–plastic interface are:
(7) |
(8) |
Many studies have shown that deepwater shallow sediment has a short diagenetic time. It is less-consolidated and has all physical properties of saturated soft soil [22C.W. Wang, R.H. Wang, J.Y. Pu, R.D. Chen, S.Q. Zhong, and F.D. Li, "Problems existent in deepwater cementing and its solution", Drilling & Production Technology, vol. 29, no. 3, pp. 11-14, 2006., 23Y.R. Zhou, J.H. Wang, and S.Z. Li, "Static and cyclic behaviors of remoulded deepwater sediments of Liwan in South China Sea", Marine Science Bulletin, vol. 32, no. 5, pp. 521-526, 2013.]. Therefore, the plastic zone around the borehole for shallow formation is controlled by the consolidation theory of saturated soil [24Q. Tan, B.H. Wei, J.G. Deng, C.L. Yan, and K. Zhao, "Calculation method of drilling safety density window in deepwater oil and gas fields", Journal of Oil and Gas Technology, vol. 34, no. 10, pp. 98-104, 2012.]. It is not suitable to analyze the stability of borehole wall with the traditional porous elastic theory because of the strong plasticity of shallow formation near the borehole wall. Y. Jin et al. have studied the borehole shrinkage deformation of salt gypsum formation and point out that the small deformation theory cannot describe the deformation law of the elastic-plastic formation, and the elastic-plastic deformation law of soft rock must be solved by large deformation theory [25Y. Jin, and M. Chen, "A new approach for controlling tight hole in salt formation", Chinese Journal of Rock Mechanics and Engineering, vol. 19, pp. 1111-1114, 2000.]. The study of L. Yan et al. also suggests that the deformation of the plastic zone around the borehole for deepwater shallow sediment satisfies the large deformation theory. According to the large deformation theory, we can get the following formulas:
(9) |
Where ε_{r} is the radial strain, ε_{θ} is the tangential strain; u is the formation displacement. Deepwater shallow sediment is similar to saturated soil which is undrained during the drilling, so the volumetric strain of the plastic zone is 0.
(10) |
Assuming the displacement on the elastic-plastic interface is μ_{b} .
(11) |
Zhou(2002) deduced the displacement field around borehole under uniform in-situ stress [26D. Zhou, and A.K. Wojtanowicz, "Analysis of leak-off tests in shallow marine sediments", Journal of Energy Resources Technology, vol. 124, no. 4, pp. 231-238, 2002.
[http://dx.doi.org/10.1115/1.1506322] ]. Considering the difference between maximum and minimum horizontal principal stress, the displacement on the elastic-plastic interface in non-uniform in-situ stress field can be written as
(12) |
The displacement in the plastic zone can be obtained by substituting Eq. (12) into Eq. (11).
(13) |
By substituting Eqs. (3) and (7) into (10), when r=r_{w}, the borehole shrinkage rate of shallow sediments is obtained.
(14) |
The collapse pressure calculation formula under the control of borehole shrinkage rate can be written as:
(15) |
In deepwater shallow sediments, the borehole does not directly collapse. However, if the density of drilling fluid is too low, it will cause serious borehole shrinkage, which will influence wellbore stability and subsequent cementing operation. By (14) and (15), the shrinkage rate of borehole and the collapse pressure under the allowable shrinkage rate can be calculated. Many scholars have studied the collapse pressure of elastic-plastic formation on the basis of borehole shrinkage [27Z. Nie, B.R. Xia, L.Z. Zou, and J.G. Deng, "Modeling of wellbore stability for gas drilling", Natural Gas Industry, vol. 31, no. 6, pp. 71-76, 2011.-29L.S. Zhang, X.Z. Yan, X.J. Yang, Z.L. Tian, and H.L. Yang, "Elasto-plastic analysis of collapse pressure for deep coal seam drilling based on Hoek-Brown criterion", Journal of China Coal Society, vol. 38, no. 1, pp. 85-90, 2013.]. J.Q. Liu believe that when the stress state of the borehole wall exceeds the limit stress of the elastic-plastic formation, the formation of the surrounding rock will turn into the plastic state, which makes the radius of the borehole decrease continuously. When the borehole shrinkage rate reaches a critical value, the borehole wall will collapse [30J.Q. Liu, "Research on the Mechanism of Wellbore Instability and Technical Methods in Baonan Area", Ph.D thesis, Xi'an Shiyou University, Xian, China, 2013.]. Y. Jin et al. and Y. Yang et al. suggest that borehole shrinkage rate can be used as the borehole instability criterion for plastic formation [15Y. Jin, M. Chen, and C. Wu, "Borehole stability prediction before drilling for formations under the second section in an exploration well", Petroleum Exploration and Development, vol. 35, no. 6, pp. 742-745, 2008., 31Y. Yang, G.H. Wu, L.Q. Tang, and J.Y. Yang, "Analysis of wellbore stability by the elasto-plastic finite deformation theory", Journal of Harbin Engineering University, vol. 31, no. 11, pp. 1460-1464, 2010.]. The above studies indicate that the collapse pressure calculation formula of deepwater shallow sediments under the control of different shrinkage rates is reasonable and appropriate.
Different from terrestrial formation, the deepwater shallow formation around the wellbore will turn into plastic state before its fracturing. Shallow sediments are under the control of saturated soil consolidation theory. When the pore pressure is equal to the external pressure, the formation is under critical condition. The effective stress around the borehole will turn into a tensile stress if the pore pressure exceeds the critical value, and the formation will fracture if the tensile stress exceeds its tensile strength. At that time, the borehole pressure is the fracturing pressure. Therefore, the theory of excess pore pressure is introduced to analyze the changes of pore pressure.
When the borehole pressure exceeds the horizontal in-situ stress, radial stress on the borehole will be greater than the tangential stress and the borehole will turn into plastic state. When the vertical fracture appears on the borehole wall, the stress state on the borehole wall should be expressed as σ_{r}>σ_{z}>σ_{θ}. Thus, the yield function is:
(16) |
(17) |
Substituting Eq. (16) into Eq. (4), and when σ_{r}_{|}_{r=rw}=P_{w}, the tangential stress and axial stress in the plastic zone can be written as
(18) |
(19) |
Shallow sediment is composed of solid particle skeleton and pores which are full of water. Drilling will result in stress concentration around the borehole, and the external pressure will cause a balance between effective stress and pore pressure. The pore pressure generated by external pressure is called excess pore pressure [32A.A. Al-Karni, "Evaluation of shear strength of cohesionless soil due to excess pore water pressure", Arabian Journal of Geosciences, vol. 4, no. 7-8, pp. 1095-1101, 2011.
[http://dx.doi.org/10.1007/s12517-009-0112-7] , 33J. Lovisa, W.W. Read, and N. Sivakugan, "Calculating c(v) based on non-uniform initial excess pore pressure", Géotechnique, vol. 62, no. 8, pp. 741-748, 2012.
[http://dx.doi.org/10.1680/geot.11.P.055] ]. The solution of excess pore pressure is mainly based on the researches of two scholars: Skempton and Henkel [34A.W. Skempton, "The pore-pressure coefficients A and B", Géotechnique, vol. 4, no. 4, pp. 143-147, 1954.
[http://dx.doi.org/10.1680/geot.1954.4.4.143] , 35D.J. Henkel, "The relationship between the strength, pore waterpressure and volume change characteristics of saturated clays", Géotechnique, vol. 9, pp. 119-135, 1959.
[http://dx.doi.org/10.1680/geot.1959.9.3.119] ]. Yan et al. combined their calculation model and rewrote the excess pore pressure formula as
(20) |
(21) |
(22) |
Where ΔP is the excess pore pressure. β is Henkel’s pore pressure parameters, for saturated soil, β=1. A is Skempton’s pore pressure parameter under the effect of deviatoric stress, for shallow sediments, A=0.8 [36J.H. Qian, and Z.Z. Yin, Geotechnical Principles and Calculation, China Water Resource and Hydropower Press: Beijing, China, 1994., 37X.G. Wen, X. Zhang, and J. Zhou, "Influencing factors about pore pressure generated for clay under complex stress path", Chinese Journal of Geotechnical Engineering., vol. 32, no. 11, pp. 1709-1716, 2010.]. Δσ_{1} and Δσ_{3} are the major and minor principal stress increments.
According to Eq. (18) and (19), the stress increment in plastic zone induced by drilling can be written as
(23) |
(24) |
On the basis of the plane strain theory, the axial stress increment in plastic zone around the borehole can be written as [38B.S. Aadnoy, and M. Belayneh, "Elasto-plastic fracturing model for wellbore stability using non-penetrating fluids", Journal of Petroleum Science & Engineering, vol. 45, no. 3-4, pp. 179-192, 2004.
[http://dx.doi.org/10.1016/j.petrol.2004.07.006] ]:
(25) |
Substituting Eq. (23) and (24) into Eq. (25), the increment of axial stress in the plastic zone is obtained.
(26) |
By substituting Δσ_{OCT} and Δτ_{OCT} into (20), the excess pore pressure in the plastic zone can be written as:
(27) |
According to the effective stress principle, the effective stress in the plastic zone can be expressed as
(28) |
(29) |
(30) |
Where P_{p0} means the in-situ pore pressure.
The above formulas indicate that the tangential and axial effective stress will increase with the increase of the radial distance. Therefore, the fracture appears on the borehole wall initially, which conforms to the general rule of borehole fracturing [39E. Fjær, R.M. Holt, P. Horsrud, and A.M. Raaen, Petroleum Related Rock Mechanics, 2^{nd} ed. Elsevier, Amsterdam: The Netherlands, 2008.]. The tangential effective stress on the borehole wall is:
(31) |
When the vertical fracture appears, the calculation formula of fracturing pressure can be written as:
(32) |
When the horizontal fracture appears, the stress state in plastic zone satisfies σ_{r}>σ_{θ}>σ_{z}. Therefore the yield function can be written as:
(33) |
According to Eq. (25), the axial stress on the borehole wall can be obtained.
(34) |
Considering Eqs. (32) and (33), when σ_{r}=P_{w}, the tangential and axial stress of the wellhole can be expressed as
(35) |
(36) |
Then the excess pore pressure is obtained:
(37) |
When the fracture appears on the borehole wall, the tangential effective stress should meet Eq. (39).
(38) |
When the horizontal fracture appears, the calculation formula of fracturing pressure can be written as:
(39) |
The minor value of these two above fracturing pressures ought to be regarded as the fracturing pressure of shallow sediments (P_{f}):
(40) |
In order to verify the theoretical model proposed in this paper. The data of wells in a block from the South China Sea are used to case study. The physical and mechanical parameters of the formations are shown in Table 1.
The boundary of plastic zone and borehole shrinkage rate are calculated and analyzed based on the parameters of Well-1. σ_{h} is kept constant and σ_{H} is gradually increased, then the variation of the plastic zone boundary and the borehole shrinkage rate under the condition of different horizontal stress difference is obtained. From Fig. (2) we can see that the horizontal stress difference has a significant impact on the shape of the plastic zone. The plastic zone is a circular under the uniform in-situ stress, and with increase of horizontal stress difference, the plastic zone gradually changed into an ellipse. Fig. (3) shows the variation rule of borehole shrinkage rate. Many scholars believe that the borehole shrinkage rate is a fixed value, that is, the degrees of shrinkage in different directions are different. However, the calculation results in this paper show that the borehole shrinkage rate changes with θ under the non-uniform in-situ stress. With the increase of horizontal stress difference, the anisotropy of borehole shrinkage rate increases. In Fig. (3), when the horizontal stress difference reaches 2.5MPa, the borehole shrinkage rate in the direction of σ_{H} is 0.87, and that in the direction of σ_{h} is 2.13. Based on Ewy’s criterion [40R.T. Ewy, "Yield and closure of directional and horizontal wells", International journal of rock mechanics and mining sciences, vol. 30, no. 7, pp. 1061-1067, 1993.
[http://dx.doi.org/10.1016/0148-9062(93)90072-L] ], when the borehole shrinkage rate exceeds 2%, it will influence the safety of drilling. According to this criterion, the borehole wall in the direction of σ_{h} has collapsed, but the borehole wall in the direction of σ_{H} remains stable. It is proved that the anisotropy of in-situ stress has a significant influence on the stability of borehole wall.
Fig. (2) Plastic zone boundaries under different horizontal stress difference (Well-1). |
Fig. (3) Shrinkage rates of borehole under different horizontal stress difference (Well-1). |
First, the collapse pressures under different horizontal stress difference are calculated based on the based on the parameters of Well-1. It can be seen from Fig. (4) that there is a significant difference in the collapse pressure in different directions under non-uniform in-situ stress. With the increase of horizontal stress difference, the difference of collapse pressure in different directions increases gradually. Therefore, it is not appropriate to consider the formation in a uniform in-situ stress state in the evaluation of wellbore stability for deepwater shallow formation. Then, the collapse pressures of three wells are calculated. In Fig. (5), curves a, b and c represent the calculated collapse pressures of three wells, respectively. Under a 2% hole shrinking rate, the collapse pressures of three wells change periodically with θ. Take curve a as an example, collapse pressures of Well-1 are 18.4MPa and 21.1MPa in the direction σ_{h} and σ_{H}, respectively. Therefore, the collapse pressure of shallow sediments is 18.4MPa. At this depth, the pore pressure after drilling is 20.07 MPa, which is greater than the calculated collapse pressure. Through further calculation, the borehole shrinkage rate is about 1.4%, which there was no drilling accident. The field data and test result also indicated that the drilling process is safe. The calculation results of another two wells are also in agreement with the actual measured values, which verifies the accuracy of the model established in this paper.
Fig. (4) Collapse pressures under different horizontal stress difference (Well-1). |
Fig. (5) Variation of calculated collapse pressure with the angle θ. |
(a: collapse pressure of well 1,b: collapse pressure of well 2,c: collapse pressure of well 3)
Fig. (6) shows the calculation results of the fracturing pressure under different horizontal stress difference (Well-1). P_{fv} exhibits periodic variation with θ, which indicates that the horizontal stress difference has a significant impact on P_{fv}. The difference of P_{fv} in different directions increases with the increase of the horizontal stress difference. However, P_{fh} is independent of θ and not affected by the horizontal stress difference. For Well-1, P_{fh} is always less than the minimum of pressure P_{fv}, which indicates that the borehole wall of Well-1 at this depth should generate vertical fracture. In order to further verify the reasonability of the proposed model for calculating the fracturing pressure of shallow formation, the fracturing pressures of three wells are calculated. As shown in Table 2, the results are compared with the calculation results of the traditional model: elasto-plastic fracture model of Aadnoy (2004) and shallow sediment model of Kaarstad (2008). The actual measured fracture pressure was obtained from leak off tests (LOT). As shown in Table 2, the calculation results by Kaarstad’s model have a large deviation from the actual measured values, which is due to the fact that Kaarstad’s model belongs to the empirical model and does not consider the plastic state of shallow formation around the borehole. The Aadnoy’s model is based on the assumption that the fracture occurs first on the elastic–plastic interface. The calculated fracturing pressure of this model is much higher than the result of the elastic model, which is not consistent with the fact that the deepwater shallow formation has a low fracturing pressure. The error between the results of this model and measured values is also large. Besides, Aadnoy’s model indicates that when the borehole in shallow formation is broken, only horizontal fracture can be formed on the borehole wall, but no vertical fracture can be formed. But the field data show that the shallow sediment of these three wells all generate vertical fracture, which also proves the limitations of Aadnoy’s model. The error of fracture pressure calculated by the model established in this paper is the lowest, and it is mostly in agreement with the actual measured values, which verifies the accuracy of the model established in this paper. It also proves that the consideration of the influence of non-uniform in-situ stress and the introduction of excess pore pressure theory in the analysis of fracture mechanism in deepwater shallow sediments is feasible.
Fig. (6) Fracturing pressures under different horizontal stress difference(Well-1). |
In order to further analyze the particularity of wellbore stability in deepwater drilling, the borehole stability of Well-2 in the shallow formation is evaluated according to the field data. Field data show that the density of shallow formation is very low, and the mud density is about 1.1 g/cm^{3}. Together with the impact of the upper sea water, the overburden pressure of deepwater shallow formation is very low, which is 0.7 g/cm^{3} smaller than that of terrestrial formation at the same depth. As shown in Fig. (7), the pressure system of the well is normal, and the water depth has no effect on the formation pressure. The horizontal stress of this well is low, and there are some differences between the two horizontal principal stresses. The safe window of drilling fluid density of shallow formation is narrow and the fracturing pressure is very low. If the minimum horizontal stress is selected as leakage pressure, the density window between leakage pressure and collapse pressure is approximately 0.3 g/cm^{3}, which indicates that the leakage of mud is easy to occur during drilling, so it is necessary to enhance the sealing property of drilling fluid. In order to avoid instability of borehole wall caused by surge pressure, the drilling technology should be optimized based on the reasonable selection of drilling fluid density. The calculated value of collapse pressure is consistent with the practical drilling fluid density, and the fracturing pressure is consistent with the measured value, which shows the reliability of the calculation model proposed in this paper.
Fig. (7) Safety density window of drilling fluid in deepwater shallow sediment (Well-2). |
The total depth and water depth of the Well-2 are 3453m and 1367m, respectively. The lithology of the target formation is mainly the interbeds of gray fine sandstone and gray mudstone in different thickness. The actual drilling fluid density in the shallow formation is 1.03~1.15g/cm^{3}. Caliper logging data show that the borehole shrinkage rate of Well-2 in the shallow formation is 0.68~1.35%. During the process of drilling, there is no borehole wall collapses. At the drilling depth of 2079m and 2382m, leak off tests were carried out. The error analysis of the measured value of fracture pressure and the results of this paper are shown in Table 3. The prediction accuracy of each test point is more than 90%, which can meet the requirements of engineering construction. The results indicate that the safe drilling fluid density window based on the model in this paper is reliable. The collapse pressure and fracture pressure of the Well-2 are also calculated by using seismic velocity data, as shown in Fig. (7) (dotted lines). Compared with the calculated results based on the logging data, the drilling fluid density window calculated by the model in this paper is narrower and more close to the measured values, which proves that wellbore stability analysis in this paper is accurate and has an important guiding significance for the subsequent deepwater drilling.
σ_{H} | = Maximum horizontal principal stress, MPa |
σ_{h} | = Minimum horizontal principal stress, MPa |
P_{w} | = The borehole pressure, MPa |
r_{w0} | = Bit radius during the drilling, mm |
r_{w} | = Borehole radius after drilling, mm |
r_{b} | = Distance elastic-plastic interface to the borehole center, mm |
θ | = Angle from the orientation of minimum horizontal principal stress, deg |
σ_{r} | = Radial stress on the wall of borehole |
σ_{θ} | = Tangential stress on the wall of borehole |
C | = Cohesion, MPa; |
φ | = Internal friction angle, deg |
σ_{0 } | = Uniaxial compressive strength (UCS), MPa |
ε_{r} | = The radial strain |
ε_{θ} | = The tangential strain |
u | = Is the formation displacement |
μ_{b} | = Displacement on the elastic-plastic interface |
μ_{rb} | = Displacement on the elastic-plastic interface |
ΔP | = Excess pore pressure, MPa |
β | = Henkel’s pore pressure parameters, for saturated soil, β=1 |
A | = Skempton’s pore pressure parameter under the effect of deviatoric stress |
Δσ_{1} | = Major principal stress increments |
Δσ_{3} | = Minor principal stress increments |
P_{p0} | = In-situ pore pressure, MPa |
P_{c} | = Collapse pressure of shallow sedimentary formations, MPa |
P_{fv} | = Fracture pressure of vertical fracture, MPa |
P_{fh} | = Fracture pressure of horizontal fracture, MPa |
P_{f} | = Fracture pressure of shallow sedimentary formations, MPa |
The authors confirm that this article content has no conflict of interest.
This work was supported by the National Science and Technology Major Project of China (2016ZX05046004-004) and National Natural Science Foundation of China (No.51490650).
[1] | J. Shelton, J.R. Smith, and A. Gupta, "Experimental evaluation of separation methods for a riser dilution approach to dual density drilling", Journal of Energy Resources and Technology, vol. 133, no. 3, pp. 611-620, 2011. [http://dx.doi.org/10.1115/1.4004968] |
[2] | B.H. Yu, C.L. Yan, J.G. Deng, S.J. Liu, Q. Tan, and K. Xiao, "Evaluation and application of wellbore stability in deep water", Oil Drilling & Production Technology, vol. 33, no. 6, pp. 1-4, 2011. |
[3] | J. Dodson, T. Dodson, and V. Schmidt, "Gulf of Mexico ‘trouble time’ creates major drilling expenses: use of cost-effective technologies needed", Offshore, vol. 64, no. 1, pp. 46-48, 2004. |
[4] | W. Jeong-Yun, and S. Hadi, "An empirical correlation for compressibility of shallow depth clays in deepwater Gulf of Mexico", In: ASME 2014 33^{rd} International Conference on Ocean, Offshore and Arctic Engineering, 2014, p. V003T10A021. [http://dx.doi.org/10.1115/OMAE2014-24223] |
[5] | L.A.S. Rocha, P. Junqueira, and J.L. Roque, "Overcoming deep and ultra deepwater drilling challenges", In: Offshore Technology Conference, 2003, p. OTC-15233-MS. [http://dx.doi.org/10.4043/15233-MS] |
[6] | A. Paknejad, J. Schubert, and M. Amani, "A new method to evaluate leak-off tests in shallow marine sediments (SMS)", In: SPE Saudi Arabia Technical Symposium, 2007, p. SPE-110953-MS. [http://dx.doi.org/10.2118/110953-MS] |
[7] | R. Chhajlani, and Z.Q. Zheng, "Utilization of geomechanics for Medusa field development, deepwater gulf of Mexico", SPE Annual Technical Conference and Exhibition, 2002 SPE-77779-MS. [http://dx.doi.org/10.2118/77779-MS] |
[8] | L. Cui, Y. Abousleiman, A.H.D. Cheng, and J.C. Roegiers, "Time dependent failure analysis of inclined boreholes in fluid-saturated formation", Journal of Energy Resources Technology, vol. 121, no. 1, pp. 31-39, 1999. [http://dx.doi.org/10.1115/1.2795057] |
[9] | A.M. Al-Ajmi, and M.H. Al-Harthy, "Probabilistic wellbore collapse analysis", Journal of Petroleum Science & Engineering, vol. 74, no. 3, pp. 171-177, 2010. [http://dx.doi.org/10.1016/j.petrol.2010.09.001] |
[10] | D.G. Ritter, and B. Grollimund, "Wellbore stability in the deepwater Bijupira & Salema fields, offshore Brazil-a probabilistic approach", In: Offshore Technology Conference, 2003, p. OTC-15205-MS. [http://dx.doi.org/10.4043/15205-MS] |
[11] | E. Kaarstad, and B.S. Aadnoy, "Fracture model for general off-shore applications", In: SPE Asia Pacific Oil & Gas Conference and Exhibition, 2006, p. SPE-101178-MS. [http://dx.doi.org/10.2118/101178-MS] |
[12] | E. Kaarstad, and B.S. Aadnoy, "Improved prediction of shallow sediment fracturing for offshore applications", Spe Drilling & Completion, vol. 23, no. 2, pp. 88-92, 2008. [http://dx.doi.org/10.2118/101178-PA] |
[13] | C.L. Yan, J.G. Deng, and X.R. Lai, "Borehole stability analysis in deepwater shallow sediments", Journal of Energy Resources Technology, vol. 137, pp. 1-6, 2015. [http://dx.doi.org/10.1115/1.4027564] |
[14] | L.S. Zhang, X.Z. Yan, X.J. Yang, Z.L. Tian, and H.L. Yang, "An elastoplastic model of collapse pressure for deep coal seam drilling based on Hoek–Brown criterion related to drilling fluid loss to reservoir", Journal of Petroleum Science and Engineering, vol. 134, pp. 205-213, 2015. [http://dx.doi.org/10.1016/j.petrol.2015.07.003] |
[15] | Y. Jin, M. Chen, and C. Wu, "Borehole stability prediction before drilling for formations under the second section in an exploration well", Petroleum Exploration and Development, vol. 35, no. 6, pp. 742-745, 2008. |
[16] | W.B. Bradley, "Mathematical concept—stress cloud can predict borehole failure", Oil & Gas Journal, vol. 77, no. 8, pp. 92-102, 1979. |
[17] | W.B. Bradley, "Failure of inclined boreholes", Journal of Energy Resources Technology, vol. 101, no. 4, pp. 232-240, 1979. [http://dx.doi.org/10.1115/1.3446925] |
[18] | F.L. Ning, "Research on Wellbore Stability in Gas Hydrae Formation", Ph.D thesis, China University of Geosciences, Beijing, China, 2005. |
[19] | Q.H. Sun, J.G. Deng, C.L. Yan, B.H. Wei, S.J. Liu, Z.L. Liu, and K. Xiao, "Calculation method for fracture pressure of deep sea shallow formation", Journal of Central South University (Science and Technology), vol. 46, no. 4, pp. 1402-1408, 2015. |
[20] | B.H. Yu, C.L. Yan, J.G. Deng, J.L. Zhou, S.J. Liu, and X.J. Xing, "Study and application of calculation model of safe drilling fluid density window—A case study of AKPO Oil field, West Africa", China Offshore Oil and Gas, vol. 24, no. 2, pp. 58-60, 2012. |
[21] | B.S. Jiang, L. Yang, and L.P. Shi, "Stress analysis of cracked surrounding rock based on Hoek-Brown criterion", Chinese Journal of Solid Mechanics, vol. 32, pp. 300-305, 2011. |
[22] | C.W. Wang, R.H. Wang, J.Y. Pu, R.D. Chen, S.Q. Zhong, and F.D. Li, "Problems existent in deepwater cementing and its solution", Drilling & Production Technology, vol. 29, no. 3, pp. 11-14, 2006. |
[23] | Y.R. Zhou, J.H. Wang, and S.Z. Li, "Static and cyclic behaviors of remoulded deepwater sediments of Liwan in South China Sea", Marine Science Bulletin, vol. 32, no. 5, pp. 521-526, 2013. |
[24] | Q. Tan, B.H. Wei, J.G. Deng, C.L. Yan, and K. Zhao, "Calculation method of drilling safety density window in deepwater oil and gas fields", Journal of Oil and Gas Technology, vol. 34, no. 10, pp. 98-104, 2012. |
[25] | Y. Jin, and M. Chen, "A new approach for controlling tight hole in salt formation", Chinese Journal of Rock Mechanics and Engineering, vol. 19, pp. 1111-1114, 2000. |
[26] | D. Zhou, and A.K. Wojtanowicz, "Analysis of leak-off tests in shallow marine sediments", Journal of Energy Resources Technology, vol. 124, no. 4, pp. 231-238, 2002. [http://dx.doi.org/10.1115/1.1506322] |
[27] | Z. Nie, B.R. Xia, L.Z. Zou, and J.G. Deng, "Modeling of wellbore stability for gas drilling", Natural Gas Industry, vol. 31, no. 6, pp. 71-76, 2011. |
[28] | Q.H. Feng, Y.F. Cheng, and J.G. Zhang, "Elastic& plastic model for borehole stability and its apllication", Petroleum Drilling Techniques, vol. 28, no. 4, pp. 8-10, 2000. |
[29] | L.S. Zhang, X.Z. Yan, X.J. Yang, Z.L. Tian, and H.L. Yang, "Elasto-plastic analysis of collapse pressure for deep coal seam drilling based on Hoek-Brown criterion", Journal of China Coal Society, vol. 38, no. 1, pp. 85-90, 2013. |
[30] | J.Q. Liu, "Research on the Mechanism of Wellbore Instability and Technical Methods in Baonan Area", Ph.D thesis, Xi'an Shiyou University, Xian, China, 2013. |
[31] | Y. Yang, G.H. Wu, L.Q. Tang, and J.Y. Yang, "Analysis of wellbore stability by the elasto-plastic finite deformation theory", Journal of Harbin Engineering University, vol. 31, no. 11, pp. 1460-1464, 2010. |
[32] | A.A. Al-Karni, "Evaluation of shear strength of cohesionless soil due to excess pore water pressure", Arabian Journal of Geosciences, vol. 4, no. 7-8, pp. 1095-1101, 2011. [http://dx.doi.org/10.1007/s12517-009-0112-7] |
[33] | J. Lovisa, W.W. Read, and N. Sivakugan, "Calculating c(v) based on non-uniform initial excess pore pressure", Géotechnique, vol. 62, no. 8, pp. 741-748, 2012. [http://dx.doi.org/10.1680/geot.11.P.055] |
[34] | A.W. Skempton, "The pore-pressure coefficients A and B", Géotechnique, vol. 4, no. 4, pp. 143-147, 1954. [http://dx.doi.org/10.1680/geot.1954.4.4.143] |
[35] | D.J. Henkel, "The relationship between the strength, pore waterpressure and volume change characteristics of saturated clays", Géotechnique, vol. 9, pp. 119-135, 1959. [http://dx.doi.org/10.1680/geot.1959.9.3.119] |
[36] | J.H. Qian, and Z.Z. Yin, Geotechnical Principles and Calculation, China Water Resource and Hydropower Press: Beijing, China, 1994. |
[37] | X.G. Wen, X. Zhang, and J. Zhou, "Influencing factors about pore pressure generated for clay under complex stress path", Chinese Journal of Geotechnical Engineering., vol. 32, no. 11, pp. 1709-1716, 2010. |
[38] | B.S. Aadnoy, and M. Belayneh, "Elasto-plastic fracturing model for wellbore stability using non-penetrating fluids", Journal of Petroleum Science & Engineering, vol. 45, no. 3-4, pp. 179-192, 2004. [http://dx.doi.org/10.1016/j.petrol.2004.07.006] |
[39] | E. Fjær, R.M. Holt, P. Horsrud, and A.M. Raaen, Petroleum Related Rock Mechanics, 2^{nd} ed. Elsevier, Amsterdam: The Netherlands, 2008. |
[40] | R.T. Ewy, "Yield and closure of directional and horizontal wells", International journal of rock mechanics and mining sciences, vol. 30, no. 7, pp. 1061-1067, 1993. [http://dx.doi.org/10.1016/0148-9062(93)90072-L] |