Abstract

CO2 sequestration is a kind of technology that can effectively mitigate the greenhouse effect, and the wellbore integrity of the storage system is the key issue to ensure successful CCUS. Failure to timely diagnose tubing leakage in offshore gas wells with high CO2 will lead to annular pressure risk and wellbore corrosion problems. The annular pressure caused by tubing leakage is characterized by high pressure and rapid rise rate, resulting in wellhead jacking, gas leakage, and wellbore structure corrosion, which has become the main cause of wellbore integrity failure in offshore gas wells. Therefore, a model of heat and mass transfer and distribution in offshore gas wells was established firstly, and the physical process of gas leakage and accumulation was described based on the pinhole model and the principle of gas PVT characteristics and volume compatibility. The results show that the depth and equivalent size of the leakage point are two important factors affecting the pressure rise process and leakage rate. Taking case well parameters as an example, after inversion calculation, the annulus pressure in case well was 23.0 MPa, the leakage point equivalent diameter was 2.3 mm, the maximum leakage rate was 0.30 m3/min, and the wellbore safety risk is relatively high. Mechanical repair and chemical plugging agent are recommended to seal the tubing leakage. These risk assessment technologies provide reference for wellbore integrity design and management to reduce CO2 leakage risk caused by wellbore integrity failure.

1. Introduction

With the development of offshore oil and gas industry, most of the offshore gas well production faces high temperature, high pressure, and high corrosive gas containing downhole environment, and the phenomenon of abnormal pressure carrying in the annulus column was gradually increasing [16]. According to the U.S. Bureau of Mines, up to 60% of gas wells have annular pressure buildup phenomenon, and the failure of oil casing column leakage accounts for more than 50%, becoming the primary factor. Once a tubing leak occurs, the tubing-casing annular pressure rises sharply, accompanied by corrosive media such as hydrogen sulfide gas and carbon dioxide gas in the natural gas of the marine well fleeing into the annulus, the degree of leakage expands, and the wellbore barrier fails. If the annulus gas with high CO2 leaks to the surface, it not only causes the CO2 pollution and makes environmental problems but also leads to wellbore scrapping and endangers the lives and property of personnel. At present, continuous annular pressure analysis mainly focuses on the casing annular pressure caused by cement sheath failure [710], usually using formation gas as the source of leakage, using the cement sheath as the transport channel, using the PVT conservation relationship to calculate the pressure value, and controlling the annular pressure by optimizing the cement sheath system. Some scholars have also conducted studies on the characteristics of pressure changes under gas-driven recovery tubular column puncture [1114] and completion tubular column leakage [1518]. However, the annular pressure caused by tubing leakage is different from the pressure caused by cement sheath integrity failure, which is mainly reflected in the difference of gas leakage source, leakage, and gathering process, which is not conducive to carry out relevant risk prediction and control measure research. Therefore, how to use production data such as wellbore pressure, temperature, flow rate, gas components, and liquid level to make accurate judgments on the extent of tubing-casing annulus leakage through surface diagnostic techniques has become a research challenge. This paper predicts the depth of the leak point based on the pressure balance principle, describes the pressurization process of gas leakage using a small-hole leakage model, analyzes the influence of leak point parameters on the change pattern of annular pressure, and summarizes four types of typical pressure initiation patterns of tubing-casing annular pressure by combining the characteristics of annular pressure relief recovery test curves. It provides technical support for the integrity evaluation of offshore production wells, well repair, oil recovery operations, and the management of CO2 leakage risk.

2. Wellbore Heat Mass Transfer and Temperature-Pressure Coupling Model

The difference in temperature and pressure inside and outside the tubular column of an offshore production well is the initial driver of gas leakage. Due to the relatively high temperature of the fluid extracted from the formation, during the upward transport along the wellbore, the high-temperature formation fluid will continue to dissipate heat to the low-temperature medium around the wellbore, resulting in a constant change in temperature within the production tubular column with well depth [19]. As shown in Figure 1, a microelement body of length is taken within the tubing column and the conservation of momentum and compliance with the energy conservation relationship within the microelement, and the two conservation equations are shown in where is the pressure (Pa), is the gas density in the tubing (kg/m3), is the gravitational acceleration (m/s2), is the well slope angle (°), is the friction coefficient (factorless), is the gas flow rate (m/s), is the gas-specific heat capacity (J/(kg-°C)), is the gas temperature (°C), is the gas mass flow rate (kg/s), and is the wellbore radial heat flow rate (J/s).

The friction coefficient in Equation (1) is calculated as shown in the following equation: where is the oil pipe roughness (m), is the oil pipe inner diameter (m), is the Reynolds number (dimensionless), and is the gas viscosity (Pa-s).

The heat transfer within the wellbore and between the wellbore and the formation is in accordance with the principle of radial heat flow conservation [20, 21], as shown in where is the radial heat flow rate from the center of the tubing to the edge of the wellbore (J/s), is the radial heat flow rate from the edge of the cement sheath to the formation (J/s), is the temperature at the edge of the wellbore (°C), is the radial heat transfer thermal resistance in the wellbore (m-°C/W), is the formation thermal conductivity (W/(m-°C)), is the formation temperature (°C), and is the dimensionless stratigraphic temperature (dimensionless).

Combining Equations (3) and (4) yields the radial heat flow rate as shown in

The transient factorless temperature in Equation (3) is shown in where is the dimensionless time (dimensionless), is the ratio of borehole-specific heat capacity to formation-specific heat capacity (dimensionless), is the time (s), is the formation heat spreading coefficient (m2/s), and is the borehole radius (m).

Substitute Equation (5) into Equation (1), and substitute the conservation of momentum relationship into the conservation of energy relationship; at this time, the pressure and temperature along the direction of the tube column are shown in

There is a coupling relationship between the temperature and pressure fields inside the oil pipe column. Therefore, a spatial segmentation method is used to obtain the temperature and pressure distribution within the tubular column. Firstly, the tubing is divided into segments of length along the direction from the bottom of the well to the wellhead, and the temperature and pressure are kept constant within the segments, and Equation (9) is obtained from Equation (8). where is the calculated parameter, as shown in where is the segment number, (dimensionless), and is the segment length (m).

According to the PVT nature of the gas and the law of mass conservation, the velocity change and the gas density within each segment are shown in

The compression factor is temperature-pressure dependent and can be expressed by Ehsan’s equation as shown in where , , , , , , , , , and are the constants (dimensionless); is the proposed comparison pressure (dimensionless); is the proposed comparison temperature (dimensionless); and is the relative density of gas (dimensionless).

The gas viscosity is closely related to the temperature and can be expressed by the Sartland formular as shown in where is the measured viscosity of the gas (Pa-s), is the temperature corresponding to the measured viscosity (°C), and is the constant, taken as 110.4.

3. Physical Model of Gas Leakage

After the failure of the tubular column of the offshore gas well under the effect of its own gravity, high temperature and pressure, and fluid corrosion, the gas enters the tubing-casing annulus through the leakage point of the tubular column, leaks from the inside of the tubular column to the outside annulus under the effect of the temperature and pressure difference between the inside and outside, and then floats and transports to the top of the annulus through the protective fluid of the annulus, forming a gas cavity (as shown in Figure 2). The whole physical process can be divided into three stages: tube column failure leakage, gas floating transport, and cavitation gathering and pressurization, as follows: (1)After completion of a natural gas well, the annular wellhead pressure Pco should be zero, excluding the effect of wellbore temperature(2)When a leak occurs in the production tubing column, driven by pressure, high-pressure oil and gas will flow into the annulus from the inside of the production tubing column. At the same time, this fluid will compress the existing fluid in the annulus, causing the annulus wellhead pressure Pc1 to rise(3)And when the production pipe column and the pressure of the annulus on both sides of the leak point are balanced, the leak stops and the annular pressure Pc2 rises to a steady state(4)The wellbore pressure-depth distribution diagram shows that as fluid continues to leak into the annulus, the annular pressure profile gradually shifts to the right until the annular pressure PCL at the leak point is equal to the production tubing column pressure PTL

Therefore, the following modeling assumptions are made before the calculation: (1)There is only one-dimensional stable fluid flow in the tubing(2)Heat transfer in wellbore and between wellbore and formation conforms to the principle of conservation of radial heat flow(3)The fluid thermodynamic properties of tubing are consistent at the same depth(4)The leak point below the liquid column in the annulus and the height could be measured by surface equipment(5)The model involves only one leak point

The temperature and pressure inside the tubing column are higher than the outer sleeve annulus, so the gas inside the column will leak into the sleeve annulus upon failure of the column integrity. The failure of the tubing column that causes sustained annular pressure is small, so the small-hole leakage model can be used to describe the leakage process of gas through the failure location. The small-hole model treats the flow as noncritical and critical, where the noncritical flow is related to the pressure inside and outside the leak point, while the critical flow is related to the internal pressure at the leak point only, as shown in

The critical pressure is shown in where is the leak rate (nominal) (m3/s), is the flow coefficient (dimensionless), is the pressure inside the leak point (Pa), is the equivalent area of the leak point (m2), is the nominal density of the gas (kg/m3), is the adiabatic index of the gas (dimensionless), is the molar mass of the gas (kg/mol), is the compression factor (dimensionless), is the pressure outside the leak (Pa), is the temperature at the leak (K), and is the critical pressure ratio (dimensionless).

According to Equation (15), the leak rate is related to the pressure inside and outside the leak point of the tubing column. The pressure inside the leak point is the internal pressure of the tubing column, while the external pressure consists of the pressure in the annulus and the pressure generated by the fluid inside the annulus, as shown in where is the annular pressure (Pa), is the density of gas in the annulus (kg/m3), and is the density of liquid in the annulus (kg/m3).

The gas density is lower than the annulus liquid, so it will float up from the leak point to the wellhead and finally gather in the annulus, so only the floating of the gas in the tubing-casing annulus liquid belongs to the intermediate process, and the final distribution state of the gas in the annulus is the factor that determines the size of the annular pressure. After the gas enters the annulus, it floats upward from the leak point to the wellhead and finally collects in the upper part of the annulus, in which the tubing-casing annulus conforms to the volume compatibility principle. Therefore, the volume relationship after the leak is shown in where is the volume of gas column in the annulus (m3), is the initial volume of liquid column in the annulus (m3), is the volume change of liquid column in the annulus (m3), is the volume of annulus (m3), and is the volume change of annulus (rigid in the annulus (m) and space is zero (m3)).

In the above equation, the volume of the annulus gas column is affected by temperature and pressure. Since the temperature of the annulus varies along the depth of the well, the gas volume is calculated by using the method of segment superposition, and the gas temperature within the segment is the same when the length is small enough. According to the PVT law for gases, the volume of the entire gas column is shown in where is the cumulative gas leakage (standard condition) (m3), is the number of segments (dimensionless), is the annulus segment number (dimensionless), is the standard condition pressure (Pa), is the standard condition temperature (K), is the annular temperature (K), and is the pressure in the annulus (Pa).

The volume change of the liquid occurs under the combined effect of temperature and pressure, as shown in where is the average temperature change of the annular liquid (°C), is the liquid isobaric expansion coefficient (°C-1), is the liquid isothermal compression coefficient (MPa-1), and is the annular pressure (Pa).

Equations (19) and (20) in the gas volume and liquid volume are affected by the change in the temperature of the annulus. According to the principle of radial heat flow conservation, the average temperature of the inner and outer sides of the annulus is taken as the annulus temperature, and the annulus temperature is shown in

Correspondingly, for the average segmentation of the annular fluid, the average temperature change is shown in where is the annulus temperature (°C), is the annulus outer to tubing thermal resistance (m-°C/W), is the annulus inner to tubing thermal resistance (m-°C/W), is the annular initial temperature (°C), is the annular liquid segment number (the no factor), and is the total number of annulus liquid segments (no factor).

After obtaining the temperature, pressure distribution and variation, and gas leakage, substituting Equations (19) and (20) into Equation (18) constitutes an equation for the annular pressure. In this case, there is a coupling between the annular gas leakage rate and the annular pressure in the relationship. For this reason, the total amount of gas entering the annular, when the pressure is unchanged for a shorter period of time, is shown in where is the initial gas volume in the annulus (standard condition) (m3), is the number of time segments (dimensionless), is the time step number (dimensionless), and is the time step (s).

The analysis shows that there is a coupling relationship between the temperature and pressure field inside the oil pipe column, and the total pressure and gas in the tubing-casing annulus also changes with time. Therefore, a spatial segmentation and time iteration method is used for the solution. Before a leak occurs, the liquid within the annulus undergoes thermal expansion to generate the confined annular pressure. Therefore, the initial conditions are shown in where is the time when the leak occurred (s), is the temperature at the bottom of the well (°C), is the pressure at the bottom of the well (Pa), and is the pressure generated by thermal expansion of the annular fluid before the leak (Pa).

After obtaining the initial conditions, the gas physical parameters under the temperature-pressure conditions of segment are used to calculate the pressure temperature of each subsegment of segment , and then, the gas physical parameters of segment are recalculated, and the difference between the two calculated rates is used as the rate change value, and then, the temperature pressure of the new segment is solved, and the error analysis is performed with the previous one, and the calculation of the next subsegment up to the wellhead is continued after meeting the error. Then, calculate the leak rate according to the location of the leak point, obtain the volume of gas inside the annulus, and then solve for the annular pressure, after which make time continue and cycle the calculation until the ratio of pressure difference inside and outside the leak point reaches a predetermined value.

4. Case Study

The case well is a high-temperature and high-pressure well with gas reservoir temperature of 150°C, CO2 content of 45%, and CO2 partial pressure of 20.5 MPa. There is obvious sustained annular pressure in tubing-casing annulus. In order to avoid the damage to people and property which caused by the CO2 leakage from the annulus, it is necessary to find an accurate diagnosis method for downhole leakage. The relevant calculated parameters are shown in Table 1.

As shown in Figure 3, the annular pressure and gas volume rose gradually after the leakage occurred, but the rising rate slowed down gradually. The tubing-casing annular pressure rose to 21.24 MPa within 1400 min, of which the rate of pressure rise was faster in the first 500 min, and the volume of annulus gas grew faster. This shows that the continuous annular pressure can lead to the accumulation of a large amount of gas in the annulus, and the risks caused by gas accumulation should be paid attention to in the field operation. At the same time, the gas leakage changes with time from critical leakage to noncritical leakage, and the leakage rate remains high and basically unchanged during the critical leakage stage and gradually decreases during the noncritical leakage stage.

As shown in Figure 4, the depth of the leak point influences both the process of pressure rise in the annulus and the rate of leakage. At the early stage of pressure generation, the rate of rise of the annular pressure caused by different depth leak points is basically the same, after which the sustained annular pressure caused by shallow leak points has the characteristics of fast rise rate, long pressure rise period, and high-pressure value. The shallow leakage point has the characteristics of rapid annular pressure rise, long pressure rise period, and high-pressure value, which is caused by the long duration of noncritical leakage phase of shallow leak point. Taking the four curves of 1000 m, 2000 m, 3000 m, and 4000 m depth leak points in Figure 4 as an example, the highest-pressure value at the leak point is 23.3 MPa and the lowest is 15.2 MPa, respectively. The above analysis shows that the sustained annular pressure caused by the shallow leak point is more hazardous than the deep leak point.

As shown in Figure 5, the highest values of the sustained annular pressure under different leak point aperture conditions did not change, but the sustained annular pressure in the case of larger apertures rose quickly and had a short rise period, and the initial value of its leakage rate was large and fell quickly. The above analysis shows that the increase in leak point orifice diameter will not change the final pressure value but will lead to a significant increase in the leak rate and pressure rise. The above analysis shows that the increase in the diameter of the leak point will not change the final pressure value but will lead to a significant increase in the leak rate and pressure rise rate.

According to the above analysis, the depth of the leak point and the hole size are two important factors that affect the pressure rise process and the leak rate. According to the depth of the leak point and the size of the leak point pressure, the pressure mode of the annulus can be divided into four categories: small-hole shallow leakage, small-hole deep leakage, large-hole shallow leakage, and large-hole deep leakage. As shown in Table 2, the four modes in the pressure value, pressure rise speed, pressure rise cycle, and amount of gas accumulation in the annulus have differences, taking into account the high-pressure value that affects the reliability of the wellbore safety barrier and pressure rise speed. Considering that a high-pressure value affects the reliability of the wellbore safety barrier, a fast pressure rise rate and a short rise period are unfavorable for early warning and control, and a larger volume of accumulated gas makes it more difficult to release the spray. The risks of the four models were evaluated qualitatively. Among them, the shallow leakage of large holes has the characteristics of high-pressure value, fast speed, and high gas volume, and its risk is the highest. In contrast, the small-hole deep leakage has the lowest risk. Therefore, extra attention should be paid to the integrity of the upper production tubular column of offshore gas wells to enhance its corrosion, tensile and compressive properties, and gas tightness to reduce the probability of shallow leakage from large boreholes.

Since the case well was put into production, the annulus sleeve pressure continued to rise. Testing was carried out using the annulus with pressure surface diagnostic device, and the components of annulus sample testing were basically the same as the gas phase components of the reservoir, so it was initially determined that the tubing was leaking. According to the construction requirements of B-B test, the annular pressure relief process was operated, as shown in Figure 6. The well annular pressure start-up pattern was consistent with the “small-hole shallow leakage” type characteristics, and the annular pressure was stable at about 23.0 MPa, with relatively high wellbore safety risk. Through fitting calculations, it was diagnosed that the pore size at the leak point of the well was 2.3 mm and the maximum leak rate was 0.30 m3/min, which was lower than the critical standard of 0.42 m3/min [16].

Due to the high difficulty of deep gas well annular pressure management and expensive well repair costs, the well-adopted ultrafine calcium carbonate leak plugging agent was used to reduce the degree of annular pressure fracture and regular injection of annular protection fluid to prevent stress corrosion and other methods to comprehensively manage annular pressure. The well has been able to ensure the smooth and safe production of gas.

5. Conclusions

(1)A model was developed to calculate the sustained annular pressure in the annular space of tubing-casing of an offshore high-temperature and high-pressure gas well with high CO2. The analysis shows that the temperature and pressure differences between the inner and outer annuli of the tubular column provide the driving force for gas leakage generation. After the integrity failure of the tubular column, the annular pressure and gas volume gradually increase, but the rate of increase gradually slows down. As the pressure difference between the inside and outside of the leak point decreases, the leak rate gradually decreases, and the risk caused by gas accumulation should be considered in field operations(2)The depth of the leak point and the hole diameter are two important factors that affect the pressure rise process and the rate of leakage. Circumferential pressure rise mode is divided into four categories: small-hole shallow leakage, small-hole deep leakage, large-hole shallow leakage, and large-hole deep leakage. Among them, large-hole shallow leakage has the characteristics of high-pressure value, fast speed, and high gas volume, and its risk is the highest, while small-hole deep leakage has the lowest risk. Therefore, extra attention should be paid to the integrity of the upper production tubular column of the offshore gas wells to enhance its corrosion resistance, tensile and compressive properties, and gas tightness to reduce the probability of shallow leakage from large boreholes(3)The model involves only one leak point calculation method; in the actual working conditions, there may be multiple leak points in the oil pipe column, which will lead to more complex calculation of leak point depth, size, and leak rate; this model can be used as the superposition effect of multiple leak points and provide the judgment basis for the diagnosis of annular pressure on a macrolevel

Data Availability

The authors confirm that the data supporting the findings of this study are available within the article and its supplementary materials.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors at the China University of Petroleum, SINOPEC Research Institute of Petroleum Engineering Co., Ltd., and State Key Laboratory of Coal Resources and Safe Mining, China University of Mining and Technology would like to acknowledge the support under the grant by the National Key Research and Development Program of China (No. 2022YFC28064001).