Gas production data processing method for deviated well in bottom water gas reservoir

Abstract

This application discloses a method for processing gas production data of deviated wells in bottom-water gas reservoirs. The method includes: acquiring the mechanical parameters of the target bottom-water gas reservoir, the mechanical parameters including at least rock permeability, liquid viscosity, gas viscosity, liquid pressure, gas pressure, liquid density, and gas density; determining a second Green's formula based on the mechanical parameters, the second Green's formula describing the volume of oil and gas space in the target bottom-water gas reservoir; performing dimensionality reduction processing on the second Green's formula to obtain a first gas production model, and based on the first gas production model, inferring a first reference gas production of the deviated well in the target bottom-water gas reservoir, the first reference gas production describing the critical gas production of the deviated well under ideal production conditions. The technical solution provided by this application can improve the accuracy of inferring the critical gas production of deviated wells in bottom-water gas reservoirs.

Description

Technical Field

[0001] This application belongs to the field of production data processing technology for inclined wells in bottom water gas reservoirs, and particularly relates to a method for processing gas production data of inclined wells in bottom water gas reservoirs. Background Technology

[0002] In bottom-water gas reservoirs, the bottom contains a large amount of bottom water (i.e., groundwater). During oil and gas extraction from these reservoirs, the pressure in the wellbore of the inclined well decreases, causing the bottom water to rise and leading to water intrusion. This not only reduces the oil and gas production of the inclined well but may also damage the equipment inside the well. When the daily gas production of the inclined well in a bottom-water gas reservoir is lower than the critical production rate, the rate of bottom water rise slows down, reducing the probability of water intrusion. Conversely, when the daily gas production of the inclined well in a bottom-water gas reservoir is higher than the critical production rate, the rate of bottom water rise accelerates, increasing the probability of water intrusion. Therefore, accurately predicting the critical production rate of the inclined well in a bottom-water gas reservoir is a technical problem that urgently needs to be solved. Summary of the Invention

[0003] The embodiments of this application provide a method for processing gas production data of inclined wells in bottom-water gas reservoirs, which can improve the accuracy of inferring the critical gas production of inclined wells in bottom-water gas reservoirs.

[0004] Other features and advantages of this application will become apparent from the following detailed description, or may be learned in part from practice of this application.

[0005] According to a first aspect of the embodiments of this application, a method for processing gas production data of a deviated well in a bottom-water gas reservoir is provided. The method comprises: acquiring mechanical parameters of a target bottom-water gas reservoir, the mechanical parameters including at least rock permeability, liquid viscosity, gas viscosity, liquid pressure, gas pressure, liquid density, and gas density; determining a second Green's formula based on the mechanical parameters, the second Green's formula describing the volume of oil and gas space in the target bottom-water gas reservoir; performing dimensionality reduction processing on the second Green's formula to obtain a first gas production model; and, based on the first gas production model, inferring a first reference gas production of the deviated well in the target bottom-water gas reservoir, the first reference gas production describing the critical gas production of the deviated well under ideal production conditions.

[0006] In some embodiments of this application, based on the foregoing scheme, the second Green's formula is:

[0007]

[0008] In the formula:

[0009] U = Φ g

[0010] Φ g =Φw -Δρ wg gz1

[0011] Φ w =p w +ρ w gz

[0012] Δρ wg =ρ w -ρ g

[0013] W = ln(r / r) w )

[0014] Where, Φ g Φ represents the gas phase potential of a gas. w The liquid phase potential of a liquid is represented by Δρ. wg ρ represents the density difference between the density of a liquid and the density of a gas. w ρ represents the density of a liquid. g Let g represent the gas density, g represent the gravitational acceleration, z1 represent the height of the gas-liquid interface in the target bottom water gas reservoir, and p represent the gas density. w r represents the liquid pressure. w V represents the completion radius of the deviated well, and V represents the volume of the oil and gas space in the target bottom water and gas reservoir.

[0015] In some embodiments of this application, based on the foregoing scheme, the step of reducing the dimensionality of the second Green's formula to obtain the first gas production model includes: performing a dimensionality reduction process on the second Green's formula to obtain a first dimensionality-reduced model.

[0016]

[0017] The first dimensionality reduction model is used to describe the surface area of ​​the oil and gas space in the target bottom water gas reservoir; a second dimensionality reduction process is performed on the first dimensionality reduction model to obtain the first gas production model:

[0018]

[0019] Where, q c1 The first reference gas production rate of the deviated well is represented by K, which represents the rock permeability, and Δρ is the rock permeability. wg The density difference between the liquid and gas densities is represented by g, g represents gravitational acceleration, h represents the effective thickness of the target bottom water gas reservoir, b represents the perforation thickness of the deviated well, and μ represents the density difference between the liquid and gas densities. g r represents the viscosity of the gas. e The radius of the inclined well gas supply is represented by r. w Indicates the completion radius of the inclined shaft.

[0020] In some embodiments of this application, based on the foregoing scheme, the method further includes: adding a skin factor and / or anisotropy correction coefficient to the first gas production model to determine a second gas production model, wherein the skin factor is used to characterize the wellbore contamination degree of the deviated well, and the anisotropy correction coefficient is used to correct the anisotropy of the target bottom water gas reservoir; converting the units of each physical parameter in the second gas production model to international standard units to obtain a target gas production model; and determining the critical gas production of the bottom water gas reservoir deviated well under actual production conditions based on the target gas production model, as a second reference gas production.

[0021] In some embodiments of this application, based on the aforementioned scheme, a skin factor is added to the first gas production model, and the determined second gas production model is as follows:

[0022]

[0023] Where, q c2 The value represents the critical gas production rate of the deviated well under wellbore contamination conditions, K represents the rock permeability, and Δρ represents the critical gas production rate. wg The density difference between the liquid and gas densities is represented by g, g represents gravitational acceleration, h represents the effective thickness of the target bottom water gas reservoir, b represents the perforation thickness of the deviated well, and μ represents the density difference between the liquid and gas densities. g r represents the viscosity of the gas. e The radius of the inclined well gas supply is represented by r. w The radius of the deviated well completion is represented by S, and the skin factor is represented by S.

[0024] In some embodiments of this application, based on the aforementioned scheme, a skin factor and anisotropy correction coefficient are added to the first gas production model, and the determined second gas production model is as follows:

[0025]

[0026] In the formula,

[0027] β=1 / η

[0028]

[0029] Where, q c3 This indicates the critical gas production rate considering the anisotropy of the target bottom water gas reservoir and the impact of wellbore contamination on the deviated well. K represents the rock permeability, and Δρ represents the rock permeability. wg The density difference between the liquid and gas densities is represented by g, g represents gravitational acceleration, h represents the effective thickness of the target bottom water gas reservoir, b represents the perforation thickness of the deviated well, and μ represents the density difference between the liquid and gas densities. g r represents the viscosity of the gas. e The radius of the inclined well gas supply is represented by r. w Let S represent the completion radius of the deviated well, S represent the skin factor, η represent the anisotropy correction coefficient, and K represent the surface factor.v K represents the vertical permeability of rock. h This indicates the horizontal permeability of the rock.

[0030] In some embodiments of this application, based on the foregoing scheme, the target gas production model is as follows:

[0031]

[0032] In the formula,

[0033] β=1 / η

[0034]

[0035] Where, q c4 This represents the second reference gas production rate of the deviated well, K represents the rock permeability, and Δρ wg The density difference between the liquid and gas densities is represented by h, which represents the effective thickness of the target bottom water gas reservoir, and b represents the perforation thickness of the deviated well. g μ represents the gas volume coefficient. g r represents the viscosity of the gas. e The radius of the inclined well gas supply is represented by r. w Let S represent the completion radius of the deviated well, S represent the skin factor, η represent the anisotropy correction coefficient, and K represent the surface factor. v K represents the vertical permeability of rock. h This indicates the horizontal permeability of the rock.

[0036] In some embodiments of this application, based on the foregoing scheme, the method further includes: obtaining a first reference gas production rate or a second reference gas production rate of the deviated well in the target bottom water gas reservoir; determining the daily gas production rate of the deviated well based on the first reference gas production rate or the second reference gas production rate, wherein the daily gas production rate is less than the first reference gas production rate or the second reference gas production rate.

[0037] According to a second aspect of the embodiments of this application, a computer program product is provided, the computer program product including computer instructions stored in a computer-readable storage medium and adapted to be read and executed by a processor to cause a computer device having the processor to perform an operation as described in any of the embodiments of the first aspect above.

[0038] According to a third aspect of the embodiments of this application, a computer-readable storage medium is provided, the computer-readable storage medium storing at least one computer program instruction, the at least one computer program instruction being loaded and executed by a processor to perform the operation performed by the method described in any of the embodiments of the first aspect above.

[0039] According to a fourth aspect of the present application, an electronic device is provided, the electronic device including one or more processors and one or more memories, the one or more memories storing at least one computer program instruction, the at least one computer program instruction being loaded and executed by the one or more processors to perform the operation performed by the method as described in any of the embodiments of the first aspect above.

[0040] Based on the technical solution proposed in this application, by acquiring various mechanical parameters of the target bottom-water gas reservoir (such as rock permeability, viscosity, density, and pressure of liquids and gases), a data foundation is provided for the establishment of the subsequent first gas production model and the estimation of the first reference gas production. In this way, by comprehensively collecting various mechanical parameters, the accuracy and reliability of the subsequent model establishment can be improved, thereby improving the accuracy of the estimation of the critical gas production of bottom-water gas reservoir deviated wells. In addition, by adopting the second Green's formula in the gas production data processing process, the rigor of the gas production data processing process of bottom-water gas reservoir deviated wells can be improved, thereby further improving the accuracy of the estimation of the critical gas production of bottom-water gas reservoir deviated wells.

[0041] It should be understood that the above general description and the following detailed description are exemplary and explanatory only, and do not limit this application. Attached Figure Description

[0042] The accompanying drawings, which are incorporated in and form part of this specification, illustrate embodiments consistent with this application and, together with the description, serve to explain the principles of this application. It is obvious that the drawings described below are merely some embodiments of this application, and those skilled in the art can obtain other drawings based on these drawings without any inventive effort. In the drawings:

[0043] Figure 1 A flowchart of a gas production data processing method for a bottom water gas reservoir deviated well in one embodiment of this application is shown;

[0044] Figure 2 A schematic cross-sectional view of the target bottom water gas reservoir in the vertical direction is shown in one embodiment of this application;

[0045] Figure 3 This invention provides a schematic diagram of the surface units of the gas-liquid interface of a target bottom water gas reservoir in one embodiment of the present application.

[0046] Figure 4 A schematic diagram of the structure of an electronic device according to one embodiment of this application is shown. Detailed Implementation

[0047] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, and not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0048] Furthermore, the described features, structures, or characteristics can be combined in any suitable manner in one or more embodiments. Numerous specific details are provided in the following description to give a thorough understanding of embodiments of this application. However, those skilled in the art will recognize that the technical solutions of this application can be practiced without one or more of the specific details, or other methods, components, apparatuses, steps, etc., can be employed. In other instances, well-known methods, apparatuses, implementations, or operations are not shown or described in detail to avoid obscuring various aspects of this application.

[0049] The block diagrams shown in the accompanying drawings are merely functional entities and do not necessarily correspond to physically independent entities. That is, these functional entities can be implemented in software, in one or more hardware modules or integrated circuits, or in different network and / or processor devices and / or microcontroller devices.

[0050] The flowcharts shown in the accompanying drawings are merely illustrative and do not necessarily include all content and operations / steps, nor do they necessarily have to be performed in the described order. For example, some operations / steps can be broken down, while others can be combined or partially combined; therefore, the actual execution order may change depending on the specific circumstances.

[0051] It should also be noted that the terms "first," "second," etc., in the specification, claims, and accompanying drawings of this application are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such uses of these terms can be interchanged where appropriate so that the embodiments of this application described herein can be implemented in orders other than those illustrated or described.

[0052] To enable those skilled in the art to better understand this application, the critical gas production rate of the bottom water gas reservoir deviated well proposed in this application will be briefly described below.

[0053] In bottom-water gas reservoirs, the bottom contains a large amount of bottom water (i.e., groundwater). During oil and gas extraction from these reservoirs, the pressure in the wellbore of the inclined well decreases, causing the bottom water to rise and leading to water intrusion. This not only reduces the oil and gas production of the inclined well but may also damage the equipment inside the well. When the daily gas production of the inclined well in a bottom-water gas reservoir is lower than the critical production rate, the rate of bottom water rise slows down, reducing the probability of water intrusion. Conversely, when the daily gas production of the inclined well in a bottom-water gas reservoir is higher than the critical production rate, the rate of bottom water rise accelerates, increasing the probability of water intrusion. Therefore, accurately predicting the critical production rate of an inclined well in a bottom-water gas reservoir is a technical problem that urgently needs to be solved. Based on this, the inventors of this application propose a data processing method for gas production of inclined wells in bottom-water gas reservoirs to improve the accuracy of predicting the critical production rate of inclined wells in bottom-water gas reservoirs.

[0054] Next, we will combine Figure 1 The gas production data processing method of the bottom water gas reservoir deviated well proposed in this application is described in detail.

[0055] See Figure 1 This document illustrates a flowchart of a gas production data processing method for a bottom-water gas reservoir deviated well according to one embodiment of this application. The method can be executed by a device with speculative processing capabilities, such as... Figure 1 As shown, the method may include at least steps 110 to 130:

[0056] Step 110: Obtain the mechanical parameters of the target bottom water gas reservoir. The mechanical parameters include at least rock permeability, liquid viscosity, gas viscosity, liquid pressure, gas pressure, liquid density, and gas density.

[0057] Step 120: Determine the second Green's formula based on the mechanical parameters. The second Green's formula is used to describe the volume of oil and gas space in the target bottom water gas reservoir.

[0058] Step 130: Dimensionality reduction is performed on the second Green's formula to obtain the first gas production model. Based on the first gas production model, the first reference gas production of the deviated well in the target bottom water gas reservoir is inferred. The first reference gas production is used to describe the critical gas production of the deviated well under ideal production conditions.

[0059] In this application, it should be noted that the ideal production conditions can specifically be that both gas and liquid flow stably without undergoing physicochemical changes; the rock permeability in the bottom water gas reservoir is equal everywhere; capillary pressure is ignored, and there is a clear interface between the gas and liquid; the density and viscosity of the gas and liquid are constant, that is, the temperature and pressure remain unchanged. According to actual needs, the ideal production conditions can also be other conditions, which are not specifically limited in this application.

[0060] In this application, by acquiring various mechanical parameters of the target bottom-water gas reservoir (such as rock permeability, viscosity, density, and pressure of liquids and gases), a data foundation is provided for the establishment of the subsequent first gas production model and the estimation of the first reference gas production. In this way, by comprehensively collecting various mechanical parameters, the accuracy and reliability of the subsequent model establishment can be improved, thereby improving the accuracy of the estimation of the critical gas production of the bottom-water gas reservoir deviated well. In addition, by adopting the second Green's formula in the gas production data processing process, the rigor of the gas production data processing process can be improved, thereby further improving the accuracy of the estimation of the critical gas production of the bottom-water gas reservoir deviated well.

[0061] In step 120 above, the second Green's formula can specifically be formula (1):

[0062]

[0063] In the formula:

[0064] U = Φ g (2)

[0065] Φ g =Φ w -Δρ wg gz1 (3)

[0066] Φ w =p w +ρ w gz (4)

[0067] Δρ wg =ρ w -ρ g (5)

[0068] W = ln(r / r) w (6)

[0069] Where, Φ g Φ represents the gas phase potential of a gas. w The liquid phase potential of a liquid is represented by Δρ. wg ρ represents the density difference between the density of a liquid and the density of a gas. w ρ represents the density of a liquid. g Let g represent the gas density, g represent the gravitational acceleration, z1 represent the height of the gas-liquid interface in the target bottom water gas reservoir, and p represent the gas density. w r represents the liquid pressure. w V represents the completion radius of the deviated well, and V represents the volume of the oil and gas space in the target bottom water and gas reservoir.

[0070] In this application, the second Green's formula can be specifically derived as follows: Under ideal production conditions, if the deviated well in the target bottom water gas reservoir produces gas at the critical production rate, and the pressure of the stagnant bottom water is the hydrostatic pressure, then the liquid phase potential of the bottom water is constant, which leads to the conclusion that:

[0071] Φ w =p w +ρ w gz (4)

[0072] Under ideal production conditions, since capillary pressure is negligible, the gas pressure and liquid pressure on both sides of the gas-liquid interface of the bottom water gas reservoir are equal. Therefore, at the height z1 of the gas-liquid interface in the target bottom water gas reservoir, we can conclude that:

[0073] p g (z1)=p w (z1)=Φ w -ρ w gz1 (7)

[0074] At the height z1 of the gas-liquid interface in the target bottom water gas reservoir, the following can be derived based on the gas phase potential:

[0075] φ g =φ w -Δρ wg gz1 (3)

[0076] The density difference between the liquid density and the gas density is:

[0077] Δρ wg =ρ w -ρ g (5)

[0078] In the second Green's formula, if functions U and W have continuous first-order partial derivatives with respect to V (volume of the oil and gas space in the target bottom-water gas reservoir) and S1 (surface area of ​​the oil and gas space in the target bottom-water gas reservoir), and continuous second-order partial derivatives with respect to V (volume of the oil and gas space in the target bottom-water gas reservoir), then the following equation holds:

[0079]

[0080] Formula (1) is the second Green's formula.

[0081] Among them, functions U and W can be:

[0082] U = φ g (2)

[0083] W = ln(r / r) w (6)

[0084] In this application, by applying the second Green's formula for data processing, the rigor and accuracy of gas production data processing can be improved to a certain extent, thereby improving the accuracy of critical gas production prediction for bottom water gas reservoir deviated wells.

[0085] In step 130 above, the dimensionality reduction of the second Green's formula to obtain the first gas production model can be specifically performed according to steps 131 to 132 as follows:

[0086] Step 131: Perform a dimensionality reduction on the second Green's formula to obtain the first dimensionality-reduced model:

[0087]

[0088] The first dimensionality reduction model is used to describe the surface area of ​​the oil and gas space in the target bottom water gas reservoir.

[0089] Step 132: Perform a second dimensionality reduction on the first dimensionality reduction model to obtain the first gas production model:

[0090]

[0091] Where, q c1 The first reference gas production rate of the deviated well is represented by K, which represents the rock permeability, and Δρ is the rock permeability. wg The density difference between the liquid and gas densities is represented by g, g represents gravitational acceleration, h represents the effective thickness of the target bottom water gas reservoir, b represents the perforation thickness of the deviated well, and μ represents the density difference between the liquid and gas densities. g r represents the viscosity of the gas. e The radius of the inclined well gas supply is represented by r. w Indicates the completion radius of the inclined shaft.

[0092] In this application, the dimensionality reduction of the second Green's formula can be specifically derived as follows: since both functions U and W satisfy the Laplace equation, substituting functions U and W into the second Green's formula for dimensionality reduction yields the first dimensionality reduction model as formula (8):

[0093]

[0094] In this application, the first dimensionality reduction model is subjected to a second dimensionality reduction process to obtain the first gas production model. Specifically, this can be combined with... Figure 2 The following method is used to derive the result, referring to... Figure 2 The diagram shows a cross-sectional view of the target bottom water gas reservoir in the vertical direction in one embodiment of this application.

[0095] like Figure 2As shown in the figure, segment AD represents the upper boundary of the target bottom water gas reservoir in the vertical direction, segment AB represents the right boundary of the target bottom water gas reservoir in the horizontal direction, segment BC represents the gas-liquid interface (i.e., gas-liquid interface) in the target bottom water gas reservoir, and segment CD represents the wellbore of the deviated well in the target bottom water gas reservoir.

[0096] In segment AB, the following formulas (10) to (13) can be derived:

[0097] ln(r / r w )=ln(r′ e / r w (10)

[0098]

[0099] dS1=2πr′ e dz (13)

[0100] In segment BC, such as Figure 3 As shown, a schematic diagram of the surface unit of the gas-liquid interface of the target bottom water gas reservoir in one embodiment of this application is shown, from which the following formulas (14) to (16) can be derived:

[0101]

[0102] In section CD, since water intrusion in the wellbore is mainly concentrated in the middle and lower sections, the middle and upper sections of the wellbore can be considered as a vertical well, and the following formulas (17) to (19) can be derived:

[0103] ln(r / r w )=0 (17)

[0104]

[0105] dS=2πr w dz (19)

[0106] In segment AD, the following formulas (20) to (21) can be derived:

[0107]

[0108] Substituting the above formulas (10) to (21) into the first dimensionality reduction model, i.e., formula (8), and simplifying and rearranging, we can obtain formula (22):

[0109]

[0110] At the well point, formulas (23) to (24) can be derived:

[0111] r = r w,z1=h′-b cosθ (23)

[0112] Φ g =Φ w -Δρ wg gz1 (24)

[0113] At the right boundary of the target bottom water gas reservoir, formulas (25) to (27) can be derived:

[0114] r = r' e ,z1=0,Φ g =Φ w (25)

[0115] h′=hb cosθ / 2 (26)

[0116] r′ e =r e +b sinθ (27)

[0117] Substituting formulas (23) to (27) into formula (22) and simplifying, we obtain formula (9), which is the first gas production model.

[0118]

[0119] In this application, by performing dimensionality reduction on the second Green's formula, the advantage is that the first gas production model can be accurately derived. The first gas production model can accurately predict the critical gas production of the deviated well in the target bottom water gas reservoir under ideal production conditions, i.e., the first reference gas production, thereby further improving the accuracy of the prediction of the critical gas production of the deviated well in the bottom water gas reservoir.

[0120] In the gas production data processing method for bottom water gas reservoir deviated wells proposed in this application, the method can also be specifically executed according to the following steps 140 to 160:

[0121] Step 140: Add a skin factor and / or anisotropy correction coefficient to the first gas production model to determine the second gas production model. The skin factor is used to characterize the wellbore contamination level of the deviated well, and the anisotropy correction coefficient is used to correct the anisotropy of the target bottom water gas reservoir.

[0122] Step 150: Convert the units of each physical parameter in the second gas production model into international standard units to obtain the target gas production model.

[0123] Step 160: Determine the critical gas production rate of the bottom water gas reservoir deviated well under actual production conditions based on the target gas production rate model, and use it as the second reference gas production rate.

[0124] In this application, a second gas production model is obtained by adding a skin factor and / or anisotropy correction coefficient to the first gas production model. In this way, by taking into account the wellbore contamination level of the deviated well and the anisotropy of the target bottom water gas reservoir, the rationality of the second gas production model can be improved, making the inferred critical gas production closer to the critical gas production under actual production conditions. In addition, by converting the units of each physical parameter in the second gas production model to international standard units, the units in the inference process can be unified, which can improve the convenience of the inference process and also improve the accuracy of the inference of the critical gas production of the deviated well in the bottom water gas reservoir.

[0125] In step 140 above, the epidermal factor is added to the first gas production model, and the determined second gas production model is formula (28):

[0126]

[0127] Where, q c2 The value represents the critical gas production rate of the deviated well under wellbore contamination conditions, K represents the rock permeability, and Δρ represents the critical gas production rate. wg The density difference between the liquid and gas densities is represented by g, g represents gravitational acceleration, h represents the effective thickness of the target bottom water gas reservoir, b represents the perforation thickness of the deviated well, and μ represents the density difference between the liquid and gas densities. g r represents the viscosity of the gas. e The radius of the inclined well gas supply is represented by r. w The radius of the deviated well completion is represented by S, and the skin factor is represented by S.

[0128] In this application, a skin factor is added to the first gas production model, which is specifically derived using the following method.

[0129] like Figure 2 As shown, applying the second Green's theorem to V2 yields formula (29):

[0130]

[0131] Applying the second Green's theorem to V1, we obtain formula (30):

[0132]

[0133] Formula (29) can be simplified to:

[0134]

[0135] Formula (30) can be simplified to:

[0136]

[0137] Adding formulas (31) and (32) together, we get formula (33):

[0138]

[0139] Among them, epidermal factors

[0140] After processing, a second gas production model (28) can be obtained by adding the epidermal factor to the first gas production model:

[0141]

[0142] In step 140 above, the skin factor and anisotropy correction coefficient are added to the first gas production model, and the determined second gas production model is formula (34):

[0143]

[0144] In the formula,

[0145] β=1 / η (35)

[0146]

[0147] Where, q c3 This indicates the critical gas production rate considering the anisotropy of the target bottom water gas reservoir and the impact of wellbore contamination on the deviated well. K represents the rock permeability, and Δρ represents the rock permeability. wg The density difference between the liquid and gas densities is represented by g, g represents gravitational acceleration, h represents the effective thickness of the target bottom water gas reservoir, b represents the perforation thickness of the deviated well, and μ represents the density difference between the liquid and gas densities. g r represents the viscosity of the gas. e The radius of the inclined well gas supply is represented by r. w Let S represent the completion radius of the deviated well, S represent the skin factor, η represent the anisotropy correction coefficient, and K represent the surface factor. v K represents the vertical permeability of rock. h This indicates the horizontal permeability of the rock.

[0148] In this application, by multiplying the anisotropy correction coefficient η in formula (28) and rearranging it, the second gas production model (34) determined by adding the skin factor and the anisotropy correction coefficient to the first gas production model can be obtained:

[0149]

[0150] In this application, by simultaneously considering the wellbore contamination level of the deviated well and the anisotropy of the target bottom water gas reservoir, the rationality of the second gas production model can be improved, making the predicted critical gas production closer to the critical gas production under actual production conditions, thereby improving the accuracy of the prediction of the critical gas production of the deviated well in the bottom water gas reservoir.

[0151] In step 150 above, the target gas production model can specifically be formula (37):

[0152]

[0153] In the formula,

[0154] β=1 / η (35)

[0155]

[0156] Where, q c4 This represents the second reference gas production rate of the deviated well, K represents the rock permeability, and Δρ wg The density difference between the liquid and gas densities is represented by h, which represents the effective thickness of the target bottom water gas reservoir, and b represents the perforation thickness of the deviated well. g μ represents the gas volume coefficient. g r represents the viscosity of the gas. e The radius of the inclined well gas supply is represented by r. w Let S represent the completion radius of the deviated well, S represent the skin factor, η represent the anisotropy correction coefficient, and K represent the surface factor. v K represents the vertical permeability of rock. h This indicates the horizontal permeability of the rock.

[0157] In this application, by converting the units of each physical parameter in the second gas production model into international standard units, the units in the prediction process can be unified. This not only improves the convenience of the prediction process, but also improves the accuracy of the prediction of the critical gas production of the bottom water gas reservoir deviated well.

[0158] In the gas production data processing method for bottom water gas reservoir deviated wells proposed in this application, the method can also be specifically executed according to the following steps 170 to 180:

[0159] Step 170: Obtain the first or second reference gas production of the deviated well in the target bottom water gas reservoir.

[0160] Step 180: Determine the daily gas production of the deviated well based on the first reference gas production or the second reference gas production, wherein the daily gas production is less than the first reference gas production or the second reference gas production.

[0161] In this application, it is understood that if the daily gas production of the deviated well is higher than the first reference gas production or the second reference gas production, it will accelerate the rise of bottom water in the target bottom water gas reservoir and increase the probability of water invasion of the deviated well.

[0162] In this application, the daily gas production is determined by the first reference gas production or the second reference gas production, and the production of the deviated well is guided according to the daily gas production. This can ensure that the deviated well of the target bottom water gas reservoir can operate normally, slow down the rate of bottom water rise, and reduce the probability of water invasion.

[0163] Based on the technical solution proposed in this application, by acquiring various mechanical parameters of the target bottom-water gas reservoir (such as rock permeability, viscosity, density, and pressure of liquids and gases), a data foundation is provided for the establishment of the subsequent first gas production model and the estimation of the first reference gas production. In this way, by comprehensively collecting various mechanical parameters, the accuracy and reliability of the subsequent model establishment can be improved, thereby improving the accuracy of the estimation of the critical gas production of the bottom-water gas reservoir deviated well. In addition, by adopting the second Green's formula in the gas production data processing process, the rigor of the gas production data processing process can be improved, thereby further improving the accuracy of the estimation of the critical gas production of the bottom-water gas reservoir deviated well.

[0164] Based on the same inventive concept, embodiments of this application provide a computer program product, the computer program product including computer instructions stored in a computer-readable storage medium and adapted to be read and executed by a processor to cause a computer device having the processor to perform the operations performed as described above.

[0165] Based on the same inventive concept, embodiments of this application provide a computer-readable storage medium storing at least one computer program instruction, which is loaded and executed by a processor to perform the operations described above.

[0166] Figure 4 A schematic diagram of the structure of an electronic device according to one embodiment of this application is shown.

[0167] Based on the same inventive concept, embodiments of this application also provide an electronic device. (Reference) Figure 4 The diagram shows a schematic of the structure of an electronic device according to an embodiment of this application. The electronic device includes one or more memories 404, one or more processors 402, and at least one computer program (program code) stored in the memory 404 and executable on the processor 402. When the processor 402 executes the computer program, it implements the method described above.

[0168] Among them, Figure 4In this document, a bus architecture (represented by bus 400) is used. Bus 400 may include any number of interconnected buses and bridges, linking various circuits including one or more processors represented by processor 402 and memory represented by memory 404. Bus 400 may also link various other circuits such as peripheral devices, voltage regulators, and power management circuits, which are well known in the art and therefore will not be described further herein. Bus interface 405 provides an interface between bus 400 and receiver 401 and transmitter 403. Receiver 401 and transmitter 403 may be the same element, i.e., a transceiver, providing a unit for communicating with various other devices over a transmission medium. Processor 402 is responsible for managing bus 400 and general processing, while memory 404 can be used to store data used by processor 402 during operation.

[0169] The functions described herein may be implemented in hardware, software executed by a processor, firmware, or any combination thereof. If implemented in software executed by a processor, the functions may be stored as one or more instructions or codes on or transmitted via a computer-readable medium. Other examples and embodiments are within the scope and spirit of this application and the appended claims. For example, due to the nature of software, the functions described above may be implemented using software executed by a processor, hardware, firmware, hardwired, or any combination thereof. Furthermore, the functional units may be integrated into a single processing unit, or each unit may exist physically separately, or two or more units may be integrated into a single unit.

[0170] In the several embodiments provided in this application, it should be understood that the disclosed technical content can be implemented in other ways. The device embodiments described above are merely illustrative; for example, the division of units can be a logical functional division, and in actual implementation, there may be other division methods. For instance, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Furthermore, the displayed or discussed mutual coupling, direct coupling, or communication connection may be through some interfaces; the indirect coupling or communication connection between units or modules may be electrical or other forms.

[0171] The units described as separate components may or may not be physically separate. Similarly, the components of the control device may or may not be physical units; they may be located in one place or distributed across multiple units. Some or all of the units can be selected to achieve the purpose of this embodiment, depending on actual needs.

[0172] If the integrated unit is implemented as a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, or all or part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this application. The aforementioned storage medium includes various media capable of storing program code, such as a USB flash drive, read-only memory (ROM), random access memory (RAM), portable hard drive, magnetic disk, or optical disk.

[0173] The above description is merely an embodiment of this application and is not intended to limit this application. Various modifications and variations can be made to this application by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of this application should be included within the scope of the claims of this application.

Claims

1. A method for processing gas production data of a deviated well in a gas reservoir with bottom water, characterized in that, The method includes: Obtain the mechanical parameters of the target bottom water gas reservoir, wherein the mechanical parameters include at least rock permeability, liquid viscosity, gas viscosity, liquid pressure, gas pressure, liquid density, and gas density; The second Green's formula is determined based on the mechanical parameters. The second Green's formula is used to describe the volume of oil and gas space in the target bottom water gas reservoir. The second Green's formula is reduced in dimensionality to obtain the first gas production model. Based on the first gas production model, the first reference gas production of the deviated well in the target bottom water gas reservoir is inferred. The first reference gas production is used to describe the critical gas production of the deviated well under ideal production conditions.

2. The method of claim 1, wherein, The second Green's formula is: In the formula: U = Φ g Φ g = Φ w - Δρ wg gz1 Φ w = p w + ρ w gz Δρ wg = ρ w - ρ g W = ln(r / r w ) wherein Φ g represents the gas phase potential of the gas, Φ w represents the liquid phase potential of the liquid, Δρ wg represents the density difference between the liquid density and the gas density, ρ w represents the liquid density, ρ g represents the gas density, g represents the acceleration of gravity, z1 represents the height of the interface between the gas and the liquid in the target bottom water gas reservoir, p w represents the liquid pressure, r w represents the completion radius of the deviated well, V represents the volume of the oil and gas space in the target bottom water gas reservoir.

3. The method of claim 2, wherein, The dimensionality reduction of the second Green's formula to obtain the first gas production model includes: Performing a dimensionality reduction operation on the second Green's formula yields the first dimensionality-reduced model: The first dimensionality reduction model is used to describe the surface area of ​​the oil and gas space in the target bottom water gas reservoir; A second dimensionality reduction process is performed on the first dimensionality reduction model to obtain the first gas production model: wherein q c1 represents the first reference gas production of the inclined well, K represents the rock permeability, Δρ wg represents the density difference between the liquid density and the gas density, g represents the gravitational acceleration, h represents the effective thickness of the target bottom water gas reservoir, b represents the inclined well shooting thickness, μ g represents the gas viscosity, r e represents the gas supply radius of the inclined well, r w represents the completion radius of the inclined well.

4. The method according to claim 1, characterized in that, The method further includes: A skin factor and / or anisotropy correction coefficient are added to the first gas production model to determine the second gas production model. The skin factor is used to characterize the wellbore contamination degree of the deviated well, and the anisotropy correction coefficient is used to correct the anisotropy of the target bottom water gas reservoir. The units of each physical parameter in the second gas production model are converted into international standard units to obtain the target gas production model. Based on the target gas production model, the critical gas production of the bottom water gas reservoir deviated well under actual production conditions is determined and used as the second reference gas production.

5. The method according to claim 4, characterized in that, By incorporating the epidermal factor into the first gas production model, the determined second gas production model is as follows: wherein q c2 represents the critical gas production rate of the inclined well under the wellbore pollution condition, K represents the rock permeability, Δρ wg represents the density difference between the liquid density and the gas density, g represents the gravitational acceleration, h represents the effective thickness of the target bottom water gas reservoir, b represents the inclined well shooting thickness, μg represents the gas viscosity, r e represents the gas supply radius of the inclined well, r w represents the completion radius of the inclined well, and S represents the skin factor.

6. The method of claim 4, wherein, By incorporating the skin factor and anisotropy correction coefficient into the first gas production model, the second gas production model is determined as follows: In the formula, β=1 / η Where, q c3 This indicates the critical gas production rate considering the anisotropy of the target bottom water gas reservoir and the impact of wellbore contamination on the deviated well. K represents the rock permeability, and Δρ represents the rock permeability. wg The density difference between the liquid and gas densities is represented by g, g represents gravitational acceleration, h represents the effective thickness of the target bottom water gas reservoir, b represents the perforation thickness of the deviated well, and μ represents the density difference between the liquid and gas densities. g r represents the viscosity of the gas. e The radius of the inclined well gas supply is represented by r. w Let S represent the completion radius of the deviated well, S represent the skin factor, η represent the anisotropy correction coefficient, and K represent the surface factor. v K represents the vertical permeability of rock. h This indicates the horizontal permeability of the rock.

7. The method of claim 6, wherein, The target gas production model is as follows: In the formula, β=1 / η Where, q c4 This represents the second reference gas production rate of the deviated well, K represents the rock permeability, and Δρ wg The density difference between the liquid and gas densities is represented by h, which represents the effective thickness of the target bottom water gas reservoir, and b represents the perforation thickness of the deviated well. g μ represents the gas volume coefficient. g r represents the viscosity of the gas. e The radius of the inclined well gas supply is represented by r. w Let S represent the completion radius of the deviated well, S represent the skin factor, η represent the anisotropy correction coefficient, and K represent the surface factor. v K represents the vertical permeability of rock. h This indicates the horizontal permeability of the rock.

8. The method of claim 4, wherein, The method further includes: Obtain the first or second reference gas production rate of the deviated well in the target bottom water gas reservoir; Based on the first reference gas production rate or the second reference gas production rate, the daily gas production rate of the deviated well is determined, wherein the daily gas production rate is less than the first reference gas production rate or the second reference gas production rate.

9. A computer program product, characterised in that, The computer program product includes computer instructions stored in a computer-readable storage medium and adapted to be read and executed by a processor to cause a computer device having the processor to perform the method of any one of claims 1 to 8.

10. An electronic device, comprising: The electronic device includes one or more processors and one or more memories, wherein at least one piece of program code is stored in the one or more memories, and the at least one piece of program code is loaded and executed by the one or more processors to perform the operation performed by the method as described in any one of claims 1 to 8.

Images