Wellbore stability is important to deep-sea drilling (DSD) for exploiting offshore oil and gas resources. Accurate prediction of the mechanical behaviors of rock formation around the wellbore can effectively prevent wellbore failure via adjusting mud weight. In this paper, a fully coupled thermo-poro-elastic (THM) model is developed by combining fluid circulation in the wellbore and rock deformation in the reservoir. Transient heat transfer in the wellbore/reservoir and heat exchange between the tubing, annulus, riser, casing, cement, seawater and rock formation are considered. By decomposing the problem into axisymmetric and antisymmetric sub-problems, the analytical solutions for transient temperature, pressure and stress in the rock formation and fluid temperature changes in the drill pipe, annulus and seawater near the riser are obtained and validated against with the numerical solutions, existing fluid circulation model, field measurement data and the existing model without considering fluid circulation. After that, a typical case study is carried out to show the mechanical behavior around the wellbore. It is found that the proposed fully coupled THM model has a high accuracy with a relative error of around 1.2%. There is a large difference in the predicted bottomhole temperature (BHT) and thus the temperature in the reservoir between the models with and without considering fluid circulation in the wellbore. An assumption of constant temperature boundary condition may lead to overestimated or underestimated results depending on the selection of the assigned wellbore temperature. In the DSD case studied, the rock formation at the upper part of the reservoir (Z < 0.40) is continuously heated and the rock formation at the lower part (Z > 0.6) of the reservoir is continuously cooled during fluid circulation, which thus results in a complex heat transfer and thus mechanical behavior in the wellbore/reservoir system, especially near the wellbore. The bottomhole temperature (BHT) decreases by 38.9 K due to cold fluid circulation from 3.68 min to 25.68 days, which can lead to 14.59 MPa tensile stress that cannot be ignored for wellbore stability analysis. The proposed analytical model provides a useful tool to predict the mechanical behaviors such as the transient temperature, pressure, displacement and stresses in the wellbore/reservoir system during fluid circulation in DSD. It runs efficiently in terms of seconds per case study in a portable laptop computer and can also work a benchmark for other numerical methods. However, the fully coupled THM model and its analytical solution are based on several assumptions, such as the heat generated by the drill bit rotation is ignored and only the liquid phase is considered, which requires further investigation.
Read full abstract