Gas lift drainage assisting optimization method for horizontal well with high water yield

Abstract

The invention discloses a gas lift drainage-aiding optimization method for a high-water-yield horizontal well, belongs to the technical field of oil-gas field development, and aims to solve the problems that in the prior art, oil pipe annulus wellbore pressure distribution of a high-water-yield gas well is unknown, collaborative optimization of gas lift and annulus drainage is little, and a liquid drainage channel production system lacks system optimization. The method comprises the following steps: firstly, collecting well bore structure and production data of a target gas well, introducing dimensionless times in a Froude number form to characterize the working condition of the high-pressure gas well, respectively establishing oil pipe and annulus wellbore pressure drop models for the high-water-yield gas well based on experimental test data, introducing a node analysis method, comparing two production modes of oil pipe and annulus production, and establishing an oil pipe and annulus pressure drop model for the high-water-yield gas well. A production channel capable of effectively reducing the flowing bottomhole pressure is selected, and then the reasonable gas injection amount is determined according to the TPR curve characteristics under the oil pipe or annulus production condition. The method does not need to replace a tubular column, is simple and convenient to operate, accurate in pressure calculation and high in adaptability, and provides theoretical and technical support for process optimization of a gas well with high water yield.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of oil and gas field development, and particularly relates to a gas lifting and cleanup optimization method for a high-yield water horizontal well. BACKGROUND

[0002] After entering the middle and late production stage, high-yield water gas wells generally face the problems of formation energy attenuation and serious wellbore liquid loading, which leads to an increase in bottom hole back pressure, blocked gas flow, and a significant reduction in production, seriously affecting the efficient production of gas wells; as a mainstream cleanup method, gas lifting can effectively lift the liquid loading by injecting gas into the wellbore to supplement the energy, and is suitable for gas wells of different well depths and liquid production scales. The determination of the gas injection rate is the key to improving the efficiency of gas lifting, and insufficient gas injection rate will lead to insufficient lifting energy and difficulty in discharging liquid loading, while excessive gas injection rate will increase the gas-liquid mixed flow friction, which has a negative impact on the production of gas wells and causes waste of gas lifting resources.

[0003] In the method of nodal system analysis, the IPR curve is the formation inflow dynamic curve, which describes the relationship between the gas supply capacity of the gas layer to the wellbore and the bottom hole flowing pressure, and reflects the productivity characteristics of the formation; the TPR curve is the tubing flow dynamic curve, which reflects the relationship between the bottom hole flowing pressure and the gas rate when the gas flows in the wellbore, and embodies the resistance characteristics of the tubing flow; by analyzing the two curves and their intersection, the production performance can be explained and the gas injection rate can be optimized. The drawing of the TPR curve is closely related to the calculation of the wellbore pressure drop, and for high-yield water gas wells, the conventional wellbore pressure drop model often cannot meet the precision requirements of the planned gas lifting system, and the low-pressure experimental data is needed to represent the high-pressure flow state of the gas well, and through experimental simulation, a wellbore pressure drop model suitable for high-yield water conditions is established. In addition, high-yield water gas wells usually use tubing production during the production process, and due to the smaller flow area, the flow friction is larger at high gas-liquid flow rate, which easily increases the bottom hole flowing pressure and reduces the gas well production; while the annular production provides a new idea for the production optimization of high-yield water gas wells due to its larger flow area and lower friction.

[0004] The invention patent (CN120145937B) divides the wellbore into two parts of liquid accumulation section and non-liquid accumulation section, and establishes a wellbore pressure drop model considering the influence of liquid accumulation, but does not consider the influence of pipe diameter on liquid holdup, and the method of testing liquid holdup in the field is not clear; the invention patent (CN117738624A) determines the optimal gas lift injection rate by coupling the optimal gas injection rate of the minimum wellbore energy loss, the minimum gas injection rate of critical liquid carrying and the minimum gas injection rate of reservoir-well coupling, but lacks pertinence for pressure drop calculation, and is difficult to adapt to high-yield water conditions; the invention patent (CN120211694A) determines the injection working gas volume suitable for the current working condition by judging the relationship between liquid volume, gas volume and injection time and judging the relationship between liquid volume, gas volume and starting pressure through field working parameters, which focuses on calculating the gas injection rate at the start of the gas well, and does not involve the optimization calculation of the gas injection rate of the gas well; the invention patent (CN120373534A) adopts the Gray model to calculate the wellbore pressure gradient, obtains the outflow dynamic curve of the bottom hole flowing pressure and the liquid production rate, and dynamically optimizes the oil nozzle size through the node system analysis method; therefore, in the gas lift production process of the high-yield water gas well, the annular wellbore pressure distribution of the tubing is not clear, the gas lift and annular drainage optimization method is less, and the production system of different drainage channels lacks a systematic optimization method.

[0005] Therefore, the present application provides a gas lift assisted drainage optimization method for high-yield water horizontal wells, which uses the dimensionless number form of the Froude number to characterize the high-pressure gas well working condition, establishes a tubing and annular wellbore pressure drop model for high-yield water gas wells based on experimental test data, introduces the node analysis method, compares the two production modes of tubing and annulus production, selects the production channel that can effectively reduce the bottom hole flowing pressure, and then determines the reasonable gas injection rate according to the TPR curve characteristics under the tubing or annular production conditions, thereby providing theoretical and technical support for the process optimization of high-yield water gas wells. SUMMARY

[0006] The present application provides a gas lift assisted drainage optimization method for high-yield water horizontal wells, which uses the dimensionless number form of the Froude number to characterize the high-pressure gas well working condition, establishes a tubing and annular wellbore pressure drop model based on experimental test data, draws IPR and TPR curves, compares the two production modes of tubing and annulus production, selects the production channel that can effectively reduce the bottom hole flowing pressure, and then determines the reasonable gas injection rate according to the TPR curve characteristics under the tubing or annular production conditions, thereby providing theoretical and technical support for the process optimization of high-yield water gas wells.

[0007] To achieve the above purpose, the specific steps of the gas lift assisted drainage optimization method for high-yield water horizontal wells according to the present application are as follows: Step 1: Collect the wellbore structure and production data of the target gas well, including well inclination angle, tubing inner diameter, tubing outer diameter, casing inner diameter, wellhead oil pressure, casing pressure, temperature, gas production rate, liquid production rate, gas phase density, liquid phase density, pressure measurement data, etc.

[0008] Step 2: Based on the production data collected in Step 1, preliminarily plot the IPR curve and TPR curve; The calculation of the IPR curve uses the empirical formula for gas well productivity commonly used in engineering: In the formula, Q SC For gas production, m 3 / d; J The gas production index, m 3 / (d×MPa 2 ); p r The mean formation pressure is expressed in MPa. p wf The bottom hole flowing pressure is in MPa. Specifically, based on the pressure measurement data collected in step 1, the formation pressure is estimated. p r and production capacity index J By sensing the bottom hole flowing pressure under different gas production rates, a preliminary formation inflow curve relating bottom hole flowing pressure to gas production rate was obtained. The plotting of the TPR curve requires the use of a wellbore pressure drop model, and the key to calculating the wellbore pressure drop is the calculation of the wellbore liquid holdup. Specifically, a tubing holdup model for high-yield water-gas wells was established. This model is based on the Mukehejee-Brill holdup model and uses the Froude number as a dimensionless criterion to more accurately achieve flow similarity between the experiment and high-pressure gas wells. It also considers the range of gas production and liquid production parameters for high-yield water-gas wells. Experimental tests were conducted, and the tubing holdup formula was obtained after fitting and correcting the experimental results. The formula is as follows: In the formula, H L tu This represents the tubing holdup, which is dimensionless. θ The inclination angle is °; N LV The apparent flow rate of the dimensionless liquid phase is dimensionless. N GV The apparent velocity of the gas phase is dimensionless and dimensionless. t 1. t 2. t 3. t 4. t 5 represents the fitting coefficient for the simulation experiment of high-water-yield gas well tubing production; Furthermore, experimental data show a clear logarithmic relationship between annular flow liquid holdup and dimensionless liquid holdup. The formula for calculating annular liquid holdup, obtained from experimental test results, is as follows: In the formula, H L an The annular liquid holdup is dimensionless. N λ , where is the dimensionless, slip-free liquid holdup, and is dimensionless; a1, a2, a3, a4, a5, a7, a8, and a9 are the fitting coefficients of the simulation experiment of annular production in high-water-yield gas wells; Among them, the dimensionless gas and liquid apparent velocities and the dimensionless no-slip holdup take into account the effects of density, flow velocity, and pipe diameter, so that the holdup calculation results can characterize high-pressure flow conditions and show good applicability in pipes of different sizes. The calculation formula is as follows: In the formula, p L The density is in the gas phase, kg / m³ 3 ; p G The density of the liquid phase is kg / m³. 3 ; v SG The apparent airflow velocity is in m / s; v SL The apparent fluid velocity is given in m / s. D Pipe diameter, in meters (m); g The acceleration due to gravity is expressed in m / s². After determining the tubing and annular holdup models, substituting them into the gas-liquid two-phase pipe flow pressure drop formula yields the tubing and annular pressure drop models, allowing for the further plotting of the TPR curve. The gas-liquid two-phase pipe flow pressure drop formula is as follows: In the formula, dp Pressure, Pa; p m For mixed density, kg / m³ 3 ; f m The coefficient of friction is a mixed friction factor, dimensionless; v m The velocity is the gas-liquid mixing velocity, in m / s; dz H is the distance along the central axis of the wellbore, in meters (m). L This is the liquid holdup, dimensionless; After clarifying the pressure drop model, based on the well structure and oil pressure, casing pressure, gas production, and fluid production data collected in step 1, the pressure drop model of the tubing or annulus production wellbore can be used to calculate the pressure value at a distance dz from the wellhead, starting from the wellhead oil pressure or casing pressure. Then, based on this pressure value, the pressure value at the next dz distance position can be calculated. The calculation is iterated step by step until the bottom hole flowing pressure is obtained. By sensing the gas production, the bottom hole flowing pressure under different gas production conditions can be obtained, and finally the TPR curve is completed.

[0009] Step 3: Adjust the formation pressure based on the intersection of the IPR curve and the TPR curve from Step 2. p r and production capacity index J This makes the gas production at the intersection of the IPR curve and the TPR curve close to the gas production of the target gas well. Specifically, by adjusting formation pressure p r and production capacity index J This ensures that the gas production at the intersection of the IPR curve and the TPR curve drawn in step 2 deviates from the actual gas production within ±2000 m³ / d, thus making the obtained IPR curve more representative of the formation's supply capacity.

[0010] Step 4: Based on the data collected in Step 1 and the annular production wellbore pressure drop model established in Step 2, plot the TPR curve under annular production conditions. Compare the gas production corresponding to the intersection of the tubing production TPR curve and the annular production TPR curve with the IPR curve obtained in Step 3, and select the production channel that can effectively reduce bottom hole flowing pressure as the recommended production channel.

[0011] Step 5: Determine the recommended gas injection rate for the target well based on the characteristics of the TPR curve of the tubing or annulus production. Step 5.1: If the recommended production method is tubing production, then calculate the gas production Q based on the lowest point of the TPR curve for tubing production. min Gas production Q at the intersection with IPR-TPR tu The difference determines the injection volume, and the recommended injection volume is Q. min -Q tu ; Step 5.2: If the recommended production method is annular production, then the gas production rate Q at the recommended injection point is determined based on the annular production TPR curve. T Gas production Q at the intersection with IPR-TPR an The difference determines the injection volume, and the recommended injection volume is Q. T -Q an ; Among them, the recommended gas injection point for annular production is selected as the point on the annular production TPR curve where the gas production increases by 10,000 cubic meters per day when the flowing pressure decreases by 1 MPa. This point can be adaptively adjusted according to the gas well production characteristics and equipment process conditions. This is because the lowest point of the annular production TPR curve corresponds to a large gas production volume, and the lowest point of the flowing pressure can no longer be used as the basis for determining the gas injection volume of the annular production. Therefore, it is stipulated that the point where the production increases by 10,000 cubic meters per 1 MPa decrease in flowing pressure is the recommended gas injection point, which is used as the gas injection judgment boundary. When the production is greater than this point, it is considered that increasing production by gas injection is no longer economical. Attached Figure Description

[0012] Figure 1 A technical roadmap for an optimized gas lift-assisted drainage method for high-yield horizontal wells; Figure 2 A schematic diagram illustrating the determination of gas injection volume during annular production; Figure 3 A schematic diagram for optimizing production channels and determining the gas injection volume during oil pipe production. Detailed Implementation

[0013] To make the objectives and calculation process of this invention clearer, the invention will be further described in detail below with reference to the accompanying drawings, so as to highlight the advantages of this invention.

[0014] like Figure 1 As shown, Figure 1 This is a technical roadmap for the present invention. First, collect the target gas well's wellbore structure and production data, including well inclination angle, tubing inner diameter, tubing outer diameter, casing inner diameter, wellhead oil pressure, casing pressure, temperature, gas production, liquid production, gas phase density, liquid phase density, pressure measurement data, etc.; based on the collected production data, initially plot the IPR curve and TPR curve; The calculation of the IPR curve uses the empirical formula for gas well productivity commonly used in engineering: In the formula, Q SC For gas production, m 3 / d; J The gas production index, m 3 / (d×MPa 2 ); p r The mean formation pressure is expressed in MPa. p wf The bottom hole flowing pressure is in MPa. Estimate formation pressure based on pressure measurement data. p r and production capacity index J By sensing the bottom hole flowing pressure under different gas production rates, a preliminary formation inflow curve relating bottom hole flowing pressure to gas production rate was obtained. The plotting of the TPR curve requires the use of a wellbore pressure drop model, and the key to calculating the wellbore pressure drop is the calculation of the wellbore liquid holdup. The formula for calculating tubing holdup is: In the formula, H L tu This represents the tubing holdup, which is dimensionless. θ The inclination angle is °; N LV The apparent flow rate of the dimensionless liquid phase is dimensionless. N GV The apparent velocity of the gas phase is dimensionless and dimensionless. t 1 , t 2 , t 3 , t 4 , t 5 ; The formula for calculating the annular liquid holdup is: In the formula, H L an The annular liquid holdup is dimensionless. N λ , where is the dimensionless, slip-free liquid holdup, and is dimensionless; a1, a2, a3, a4, a5, a7, a8, and a9 are the fitting coefficients of the simulation experiment of annular production in high-water-yield gas wells; Among them, the dimensionless gas and liquid apparent velocities and the dimensionless no-slip holdup take into account the effects of density, flow velocity, and pipe diameter, so that the holdup calculation results can characterize high-pressure flow conditions and show good applicability in pipes of different sizes. The calculation formula is as follows: , , In the formula, p L The density is in the gas phase, kg / m³ 3 ; p G The density of the liquid phase is kg / m³. 3 ; v SG The apparent airflow velocity is in m / s; v SL The apparent fluid velocity is given in m / s. D Pipe diameter, in meters (m); g The acceleration due to gravity is expressed in m / s². After determining the tubing and annular holdup models, substituting them into the gas-liquid two-phase pipe flow pressure drop formula yields the tubing and annular pressure drop models, allowing for the further plotting of the TPR curve. The gas-liquid two-phase pipe flow pressure drop formula is as follows: , In the formula, dp Pressure, Pa; p m For mixed density, kg / m³ 3 ; f m The coefficient of friction is a mixed friction factor, dimensionless; v m The velocity is the gas-liquid mixing velocity, in m / s; dz The distance along the central axis of the wellbore is in meters (m). H L This is the liquid holdup, dimensionless; Furthermore, by adjusting formation pressure p r and production capacity index J This ensures that the gas production at the intersection of the IPR curve and the TPR curve deviates from the actual gas production within ±2000 m³ / d, thus making the obtained IPR curve more representative of the formation's supply capacity.

[0015] Plot the TPR curve under annular production conditions, compare the gas production corresponding to the intersection of the tubing production TPR curve and the annular production TPR curve with the IPR curve obtained in step 3, and select the production channel that can effectively reduce bottom hole flowing pressure as the recommended production channel; a lower bottom hole flowing pressure corresponds to a larger production pressure difference from the formation to the bottom of the well, which is more conducive to the flow of formation fluid into the wellbore.

[0016] Furthermore, if the recommended production method is pipeline production, then the gas production Q at the lowest point of the pipeline production TPR curve should be used. min Gas production Q at the intersection of IPR and TPR tu The difference determines the injection volume, and the recommended injection volume is Q. min -Q tu If the recommended production method is annular production, then the recommended gas production rate Q at the injection point is based on the annular production TPR curve. T Gas production Q at the intersection with IPR-TPR an The difference determines the injection volume, and the recommended injection volume is Q. T -Q an ; like Figure 2 As shown, since the lowest point of the annular TPR curve corresponds to a large gas production volume, the lowest point of the flowing pressure can no longer be used as the basis for determining the gas injection volume of the annular production. Therefore, it is stipulated that the point where the production increases by 10,000 cubic meters per 1 MPa decrease in flowing pressure is the recommended gas injection point, which is used as the gas injection judgment boundary. When the production is greater than this point, it is considered that increasing production by gas injection is no longer economical.

[0017] Example 1: Well X1 is a high-yield water-gas well in Block X of the Sichuan Basin. The well has a casing inner diameter of 139.7 mm, an outer diameter of 60 mm, and an inner diameter of 50 mm. It uses tubing for production. As of June 2024, its average daily gas production was 16,000 cubic meters / day, daily water production was 30 cubic meters / day, oil pressure was 2.3 MPa, and casing pressure was 10.1 MPa. Based on the production data of the block corresponding to the well, the Froude number was introduced as a similarity criterion to calculate the production conditions of the block. A gas-water two-phase pipe flow experiment was designed and carried out in a targeted manner. The liquid holdup formula was fitted with the experimental data to establish a pressure drop model for the tubing and annulus wellbore. Two-phase flow experiments were conducted on the production data of high-yield water-gas wells in Block X. Based on the experimental data, coefficients in the formulas for tubing and annular liquid holdup were fitted to obtain a wellbore pressure drop calculation model for tubing and annularity. The values ​​of the fitting coefficients are as follows: t 1 = -2.05 t 2 = 0.15 t 3 = -0.14, t 4 = 0.10, t 5 = 0.15; a 1 = 0.33 a 2 = 0.04 a 3 = 0.5 a 4 = -0.04 a 5 = 0.4 a 6 = 0.51 a 7 = 0.74 a 8 = 0.29 a 9 = 1.06; like Figure 3 As shown, the IPR curve corresponding to this block is plotted according to steps 2-3. Simultaneously, the TPR curves for tubing and annulus production are calculated and plotted using the tubing and annulus wellbore pressure drop calculation model. The IPR curve intersects the tubing and annulus production TPR curves at two points. Comparing the bottomhole flowing pressure values ​​corresponding to the two intersection points, the bottomhole flowing pressure for tubing production is lower than that for annulus production; therefore, tubing production is selected, meaning there is no need to change the production method. The lowest point of the tubing production TPR curve in the figure corresponds to the gas production rate Q. min =2.2×10 4 m 3 / d, the gas production Q corresponds to the intersection of the TPR curve and the IPR curve in tubing production. tu =1.6×10 4 m 3 / d, therefore the recommended gas injection rate for this well at this time is Q. min -Q tu =6000m 3 / d.

[0018] Compared with the shortcomings and deficiencies of existing technologies, the present invention has the following beneficial effects: (1) The gas lift-assisted exhaust optimization method proposed in this invention includes the selection of production channels and determination of gas injection volume, the optimization system is more complete, and the process does not require replacement of tubing, and the process is simple and easy to operate.

[0019] (2) Establish a calculation model for tubing and annulus liquid holdup for high-yield water and gas wells, and provide a more accurate and adaptable calculation method for pressure calculation of high-yield water and gas wells.

[0020] (3) A method for determining the gas injection volume was designed for tubular production and annular production, providing specific implementation standards for the selection of production channels.

[0021] Obviously, the above description is only an implementation idea of ​​the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A method for optimizing production and determining gas injection volume in high water-cut gas wells, characterized in that, Includes the following steps: Step 1: Collect the target gas well's wellbore structure and production data, including well inclination angle, tubing inner diameter, tubing outer diameter, casing inner diameter, wellhead oil pressure, casing pressure, temperature, gas production, liquid production, gas phase density, liquid phase density, pressure measurement data, etc. Step 2: Based on the production data collected in Step 1, preliminarily plot the IPR curve and TPR curve; The plotting of the TPR curve requires a wellbore pressure drop model, and the key to calculating the wellbore pressure drop is the calculation of the wellbore fluid holdup; the formula for calculating the tubing fluid holdup is: In the formula, H L tu This represents the tubing holdup, which is dimensionless. θ The inclination angle is °; N LV The apparent flow rate of the dimensionless liquid phase is dimensionless. N GV The apparent velocity of the gas phase is dimensionless and dimensionless. t 1. t 2. t 3. t 4. t 5 represents the fitting coefficient for the simulation experiment of high-water-yield gas well tubing production; The formula for calculating the annular liquid holdup is: In the formula, H L an The annular liquid holdup is dimensionless. N λ Dimensionless, slip-free liquid holding capacity, dimensionless; a 1. a 2. a 3. a 4. a 5. a 7. a 8. a 9 represents the fitting coefficient for the simulation experiment of high-yield water-producing gas well annulus production; The formulas for calculating the dimensionless apparent flow velocities of the gas and liquid phases and the dimensionless non-slip liquid holdup are as follows: , , In the formula, ρ L The density is in the gas phase, kg / m³ 3 ; ρ G The density of the liquid phase is kg / m³. 3 ; v SG The apparent airflow velocity is in m / s; v SL The apparent fluid velocity is given in m / s. D Pipe diameter, in meters (m); g The acceleration due to gravity is expressed in m / s². After determining the tubing and annular holdup models, substituting them into the gas-liquid two-phase pipe flow pressure drop formula yields the tubing and annular pressure drop models, allowing for the further plotting of the TPR curve. The gas-liquid two-phase pipe flow pressure drop formula is as follows: , In the formula, dp Pressure, Pa; ρ m For mixed density, kg / m³ 3 ; f m The coefficient of friction is a mixed friction factor, dimensionless; v m The velocity is the gas-liquid mixing velocity, in m / s; dz The distance along the central axis of the wellbore is in meters (m). H L This is the liquid holdup, dimensionless; Step 3: Adjust the formation pressure based on the intersection of the IPR curve and the TPR curve from Step 2. p r and production capacity index J This makes the gas production at the intersection of the IPR curve and the TPR curve close to the gas production of the target gas well. Step 4: Based on the data collected in Step 1 and the annular production wellbore pressure drop model established in Step 2, plot the TPR curve under annular production conditions. Compare the gas production corresponding to the intersection of the tubing production TPR curve and the annular production TPR curve with the IPR curve obtained in Step 3, and select the production channel that can effectively reduce bottom hole flowing pressure as the recommended production channel. Step 5: Determine the recommended gas injection rate for the target well based on the characteristics of the TPR curve of the tubing or annulus production. Step 5.1: If the recommended production method is tubing production, then calculate the gas production Q based on the lowest point of the TPR curve for tubing production. min Gas production Q at the intersection with IPR-TPR tu The difference determines the injection volume, and the recommended injection volume is Q. min -Q tu ; Step 5.2: If the recommended production method is annular production, then the gas production rate Q at the recommended injection point is determined based on the annular production TPR curve. T Gas production Q at the intersection with IPR-TPR an The difference determines the injection volume, and the recommended injection volume is Q. T -Q an ; Since the lowest point of the annular production TPR curve corresponds to a large gas production volume, the lowest point of the flowing pressure can no longer be used as the basis for determining the gas injection volume of the annular production. Therefore, it is stipulated that the point where the production increases by 10,000 cubic meters per 1 MPa decrease in flowing pressure is the recommended gas injection point, which is used as the gas injection judgment boundary. When the production is greater than this point, it is considered that increasing production by gas injection is no longer economical.

Images