Method and device for calculating gas cut flow during drilling of high-pressure gas layer

Abstract

The invention discloses a method and device for calculating gas cut flow during drilling of a high-pressure gas reservoir, and relates to the technical field of data processing.The method comprises the steps that the physical property of a rock sample is detected, the rock permeability and the rock porosity corresponding to the rock sample are obtained, and the vertical depth, the borehole diameter, the drill rod outer diameter and the annulus inner wall roughness of a target well are obtained; the rock sample is obtained by sampling rocks of a target gas layer in the target well; the non-Darcy flow coefficient of the target gas reservoir is calculated based on the rock permeability and the rock porosity, and the bottom hole pressure of the target well during normal drilling is calculated based on the vertical depth, the borehole diameter, the drill rod outer diameter and the annulus inner wall roughness; and calculating the gas cut flow in the drilling process of the target well based on a preset gas cut flow calculation model in combination with the non-Darcy flow coefficient and the bottom hole pressure. By comprehensively considering a plurality of key factors including physical characteristics of rocks and basic parameters of a target well, calculation of the gas cut flow during drilling of a high-pressure gas layer is realized.

Description

Technical Field

[0001] The present disclosure relates to the technical field of data processing, and particularly to a method and device for calculating gas invasion flow rate when drilling through a high-pressure gas layer. Background Art

[0002] The development of deep wells and ultra-deep wells is an inevitable choice for the modern energy industry to cope with the increasingly depleted resources and meet the growing energy demand. Against the background of the gradual depletion of traditional shallow oil and gas resources, the development of deep wells and ultra-deep wells has become a key means to ensure energy supply, promote the optimization of the energy structure, and enhance energy security. However, during the drilling process of deep wells and ultra-deep wells, the harm of gas invasion is very significant. Since deep wells and ultra-deep wells are usually located in complex geological structures with high underground pressure and temperature, the risk of gas invading the wellbore increases significantly. Gas invasion not only causes pressure imbalance between the wellhead and the downhole, but also may cause serious accidents such as wellbore instability, blowout, and explosion. These problems are far more complex and dangerous than gas invasion accidents in ordinary wells. If effective preventive and control measures are not taken in time, the harm caused by gas invasion may rapidly amplify in a short time, leading to serious safety accidents, and even causing large-scale environmental pollution and resource waste.

[0003] During the drilling process of deep wells and ultra-deep wells, great attention must be paid to the problem of gas invasion. The key to solving the problem of gas invasion lies in the determination of the gas invasion flow rate. However, there is currently no clear method to reasonably calculate the gas invasion flow rate when drilling through a high-pressure gas layer. Therefore, providing a method for calculating the gas invasion flow rate when drilling through a high-pressure gas layer is an urgent problem to be solved at present. Summary of the Invention

[0004] The present disclosure provides a method and device for calculating the gas invasion flow rate when drilling through a high-pressure gas layer, and the main purpose is to provide a method for calculating the gas invasion flow rate when drilling through a high-pressure gas layer to realize the calculation of the gas invasion flow rate when drilling through a high-pressure gas layer.

[0005] According to the first aspect of the present disclosure, a method for calculating the gas invasion flow rate when drilling through a high-pressure gas layer is provided, which includes:

[0006] Detect the physical properties of the rock sample to obtain the rock permeability and rock porosity corresponding to the rock sample, and obtain the vertical depth, wellbore diameter, drill pipe outer diameter, and annulus inner wall roughness of the target well. The rock sample is obtained by sampling the rock of the target gas layer in the target well;

[0007] Calculate the non-Darcy flow coefficient of the target gas layer based on the rock permeability and the rock porosity, and calculate the bottom hole pressure of the target well during normal drilling based on the vertical depth, the wellbore diameter, the drill pipe outer diameter, and the annulus inner wall roughness;

[0008] Calculate the gas invasion flow rate during the drilling of the target well based on a preset gas invasion flow rate calculation model in combination with the non-Darcy flow coefficient and the bottom hole pressure.

[0009] Optionally, the formula of the preset gas invasion flow rate calculation model is expressed as:

[0010]

[0011] Among them, P f is the gas layer pressure of the target gas layer; P w is the bottom hole pressure of the target well during normal drilling; μ g is the bottom hole gas viscosity of the target well; K is the rock permeability; q g is the gas invasion flow rate; r w is the wellbore radius of the target well; h is the gas layer thickness of the target gas layer; δ is the non-Darcy flow coefficient; r f is the gas layer radius of the target gas layer; ρ g is the bottom hole gas density.

[0012] Optionally, the formula for calculating the bottom hole pressure of the target well during normal drilling based on the vertical depth, the wellbore diameter, the drill pipe outer diameter, and the annulus inner wall roughness is expressed as:

[0013]

[0014] Among them, ρ l is the drilling fluid density of the target well; g is the acceleration due to gravity in the area where the target well is located; H is the vertical depth; ξ is the friction factor; Q is the drilling fluid displacement of the target well; D is the wellbore diameter; d is the drill pipe outer diameter; ε is the annulus inner wall roughness; μ l is the drilling fluid viscosity of the target well; L is the annulus length; a is the first intermediate calculation quantity; b is the second intermediate calculation quantity.

[0015] Optionally, the formula for calculating the drilling fluid density is expressed as:

[0016] ρ l = ρ 0 + c 0 P w + c 1 T + c 2 P w T + c 3 T 2 + c 4 T 3

[0017] Among them, ρ 0 is the density of the drilling fluid under atmospheric pressure; c 0For the first correlation parameter value related to ρ 0 c 1 For the second correlation parameter value related to ρ 0 c 2 For the second correlation parameter value related to ρ 0 c 3 For the second correlation parameter value related to ρ 0 c 4 For the second correlation parameter value related to ρ 0 c, and 0 c 1 c 2 c 3 c 4 The values of c, c, c, c, c, c are all determined based on the composition and type of the drilling fluid.

[0018] Optionally, the formula for calculating the bottom-hole gas viscosity is expressed as:

[0019]

[0020] μ′ g = 1.55×10 -6 (1.8T + 32)+9.596×10 -3

[0021] Where: μ g is the bottom-hole gas viscosity; μ′ g is the viscosity of the gas under standard conditions that satisfies the corrected gas temperature; P r is the reduced pressure of the gas; T r is the reduced temperature of the gas; T is the formation temperature of the target gas reservoir; a 0 ~a 15 are different intermediate calculation quantities.

[0022] Optionally, the formula for calculating the bottom-hole gas density is expressed as:

[0023]

[0024] Where, Z sc is the gas compressibility factor under standard conditions; T sc is the formation temperature of the target gas reservoir under standard conditions; P sc is the formation pressure under standard conditions; Z is the gas compressibility factor at the bottom-hole state; b 1 ~b 8 are intermediate calculation quantities.

[0025] Optionally, the formula for calculating the non-Darcy flow coefficient of the target gas reservoir based on the rock permeability and the rock porosity is expressed as:

[0026]

[0027] Wherein, K is the permeability of the rock, and φ is the porosity of the rock.

[0028] According to a second aspect of the present disclosure, there is provided a calculating device for gas invasion flow rate when drilling through a high-pressure gas layer, including:

[0029] An acquisition unit, configured to detect the physical properties of the rock sample, obtain the rock permeability and rock porosity corresponding to the rock sample, and acquire the vertical depth, wellbore diameter, drill pipe outer diameter, and annulus inner wall roughness of the target well, wherein the rock sample is obtained by sampling the rock of the target gas layer in the target well;

[0030] A first calculation unit, configured to calculate the non-Darcy flow coefficient of the target gas layer based on the rock permeability and the rock porosity, and calculate the bottom hole pressure of the target well during normal drilling based on the vertical depth, the wellbore diameter, the drill pipe outer diameter, and the annulus inner wall roughness;

[0031] A second calculation unit, configured to calculate the gas invasion flow rate during the drilling process of the target well based on a preset gas invasion flow rate calculation model in combination with the non-Darcy flow coefficient and the bottom hole pressure.

[0032] According to a third aspect of the present disclosure, there is provided an electronic device, including:

[0033] At least one processor; and

[0034] A memory communicatively connected to the at least one processor; wherein,

[0035] The memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor so that the at least one processor can execute the method described in the foregoing first aspect.

[0036] According to a fourth aspect of the present disclosure, there is provided a non-transitory computer-readable storage medium storing computer instructions, wherein the computer instructions are used to cause the computer to execute the method described in the foregoing first aspect.

[0037] According to a fifth aspect of the present disclosure, there is provided a computer program product, including a computer program, where the computer program, when executed by a processor, implements the method described in the foregoing first aspect.

[0038] The calculation method and device for gas invasion flow rate when drilling through a high-pressure gas layer provided by the present disclosure detect the physical properties of a rock sample to obtain the rock permeability and rock porosity corresponding to the rock sample, and acquire the vertical depth, wellbore diameter, drill pipe outer diameter, and annulus inner wall roughness of the target well. The rock sample is obtained by sampling the rock of the target gas layer in the target well; calculate the non-Darcy flow coefficient of the target gas layer based on the rock permeability and the rock porosity, and calculate the bottom hole pressure of the target well during normal drilling based on the vertical depth, the wellbore diameter, the drill pipe outer diameter, and the annulus inner wall roughness; calculate the gas invasion flow rate during the drilling process of the target well based on a preset gas invasion flow rate calculation model in combination with the non-Darcy flow coefficient and the bottom hole pressure. Compared with the related art, by comprehensively considering multiple key factors, including the physical properties of the rock (such as rock permeability, rock porosity), and the vertical depth, wellbore diameter, drill pipe outer diameter, and annulus inner wall roughness of the target well, the calculation of the gas invasion flow rate when drilling through a high-pressure gas layer is realized. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] The specification drawings forming a part of the present invention are used to provide a further understanding of the present invention. The schematic embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute an improper limitation to the present invention. In the drawings:

[0040] Figure 1 It is a schematic flow chart of a calculation method for gas invasion flow rate when drilling through a high-pressure gas layer provided by an embodiment of the present disclosure;

[0041] Figure 2 It is a schematic structural diagram of a calculation device for gas invasion flow rate when drilling through a high-pressure gas layer provided by an embodiment of the present disclosure;

[0042] Figure 3 It is a schematic block diagram of an exemplary electronic device 300 provided by an embodiment of the present disclosure. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0043] The present invention will be described in detail below with reference to the drawings and in combination with embodiments. It should be noted that, without conflict, the embodiments in the present invention and the features in the embodiments can be combined with each other.

[0044] The following detailed descriptions are all exemplary descriptions, aiming to provide a further detailed description of the present invention. Unless otherwise specified, all technical terms adopted by the present invention have the same meaning as commonly understood by those of ordinary skill in the technical field to which the present invention belongs. The terms used in the present invention are only for describing specific embodiments and are not intended to limit the exemplary embodiments according to the present invention.

[0045] In addition, the terms "first", "second", etc. in the description, claims, and the above-mentioned drawings of the present disclosure are used to distinguish similar objects and do not necessarily describe a specific order or sequence. It should be understood that the data used in this way can be interchanged under appropriate circumstances so that the embodiments of the present disclosure described herein can be implemented in an order different from those illustrated or described herein.

[0046] The following describes a method and device for calculating gas invasion flow rate when drilling through a high-pressure gas layer according to an embodiment of the present disclosure with reference to the drawings.

[0047] In order to at least provide a method for calculating gas invasion flow rate when drilling through a high-pressure gas layer to achieve the calculation of gas invasion flow rate when drilling through a high-pressure gas layer. This embodiment provides a method for calculating gas invasion flow rate when drilling through a high-pressure gas layer.

[0048] Figure 1 It is a schematic flow chart of a method for calculating gas invasion flow rate when drilling through a high-pressure gas layer provided by an embodiment of the present disclosure. As Figure 1 shown, the method includes the following steps:

[0049] Step 101: Detect the physical properties of the rock sample to obtain the rock permeability and rock porosity corresponding to the rock sample, and obtain the vertical depth, wellbore diameter, drill pipe outer diameter, and annulus inner wall roughness of the target well. The rock sample is obtained by sampling the rock of the target gas layer in the target well;

[0050] Step 102: Calculate the non-Darcy flow coefficient of the target gas layer based on the rock permeability and the rock porosity, and calculate the bottom hole pressure of the target well during normal drilling based on the vertical depth, the wellbore diameter, the drill pipe outer diameter, and the annulus inner wall roughness;

[0051] As a refinement of the above step 102, the annulus inner wall roughness generally takes a value of 5×10 -5 m.

[0052] Step 103: Calculate the gas invasion flow rate during the drilling process of the target well based on a preset gas invasion flow rate calculation model in combination with the non-Darcy flow coefficient and the bottom hole pressure.

[0053] The calculation method of gas invasion flow rate when drilling through a high-pressure gas layer provided by the present disclosure detects the physical properties of a rock sample to obtain the rock permeability and rock porosity corresponding to the rock sample, and acquires the vertical depth, wellbore diameter, drill pipe outer diameter, and annulus inner wall roughness of the target well. The rock sample is obtained by sampling the rock of the target gas layer in the target well; calculates the non-Darcy flow coefficient of the target gas layer based on the rock permeability and the rock porosity, and calculates the bottom hole pressure of the target well during normal drilling based on the vertical depth, the wellbore diameter, the drill pipe outer diameter, and the annulus inner wall roughness; calculates the gas invasion flow rate during the drilling process of the target well based on a preset gas invasion flow rate calculation model in combination with the non-Darcy flow coefficient and the bottom hole pressure. Compared with the related art, by comprehensively considering multiple key factors, including the physical properties of the rock (such as rock permeability, rock porosity) and the vertical depth, wellbore diameter, drill pipe outer diameter, and annulus inner wall roughness of the target well, the calculation of the gas invasion flow rate when drilling through a high-pressure gas layer is realized.

[0054] As a refinement of the above embodiment, the formula of the preset gas invasion flow rate calculation model is expressed as:

[0055]

[0056] Wherein, P f is the gas layer pressure of the target gas layer, with the dimension of Pa; P w is the bottom hole pressure of the target well during normal drilling, with the dimension of Pa; μ g is the bottom hole gas viscosity of the target well, with the dimension of Pa·s; K is the rock permeability, with the dimension of m 2 ; q g is the gas invasion flow rate, with the dimension of m 3 / s; r w is the wellbore radius of the target well, with the dimension of m; h is the gas layer thickness of the target gas layer, with the dimension of m; δ is the non-Darcy flow coefficient, with the dimension of 1 / m; r f is the gas layer radius of the target gas layer, with the dimension of m; ρ g is the bottom hole gas density. When the preset gas invasion flow rate calculation model calculates the gas invasion flow rate, the scale of the required data parameters is small and the structure is simple, so that the time consumption for calculating the gas invasion flow rate is significantly reduced, which not only ensures the high accuracy of the result, but also greatly improves the operation speed and efficiency. This design provides strong support for the efficient decision-making of subsequent overflow disposal, provides a more scientific basis for the optimization decision-making and operation in actual engineering, and greatly promotes the progress of related technologies and the improvement of application value.

[0057] In some embodiments, the gas layer radius generally takes a value of five times the wellbore radius, that is, r f= 5r w 。

[0058] As a refinement of the above embodiment, the formula for calculating the bottom hole pressure of the target well during normal drilling based on the vertical depth, the wellbore diameter, the outer diameter of the drill pipe, and the roughness of the inner wall of the annulus is expressed as:

[0059]

[0060]

[0061] where ρ l is the density of the drilling fluid of the target well, with the dimension of kg / m3; g is the acceleration due to gravity in the area where the target well is located, with the dimension of m / s2; H is the vertical depth, with the dimension of m; ξ is the friction factor; Q is the displacement of the drilling fluid of the target well, with the dimension of m3 / s; D is the wellbore diameter, with the dimension of m; d is the outer diameter of the drill pipe, with the dimension of m; ε is the roughness of the inner wall of the annulus, with the dimension of m; μ l is the viscosity of the drilling fluid of the target well, with the dimension of Pa·s; L is the length of the annulus, with the dimension of m; a is the first intermediate calculation quantity; b is the second intermediate calculation quantity.

[0062] As a refinement of the above embodiment, the formula for calculating the density of the drilling fluid is expressed as:

[0063] ρ l = ρ 0 + c 0 P w + c 1 T + c 2 P w T + c 3 T 2 + c 4 T 3

[0064] where ρ 0 is the density of the drilling fluid under normal pressure, and ρ 0 is the density of the drilling fluid at a temperature of 338.15 K (K represents the Kelvin temperature unit) and a pressure of 101.325 KPa, with the dimension of kg / cm3; c 0 is the first correlation parameter value related to ρ 0 ; c 1 is the second correlation parameter value related to ρ 0 ; c 2 is the second correlation parameter value related to ρ 0 ; c 3 is the second correlation parameter value related to ρ 0 ; c 4 is the second correlation parameter value related to ρ 0The relevant second relevant parameter value, and c 0 , c 1 , c 2 , c 3 , c 4 The values of are all determined based on the composition and type of the drilling fluid and are dimensionless. At different drilling fluid densities, c 0 ~c 4 The values are shown in Table 1. At other drilling fluid densities, c 0 ~c 4 The values are determined by interpolation.

[0065] Table 1 shows the values of c 0 ~c 4 The values

[0066]

[0067] As a refinement of the above embodiment, the formula for calculating the bottom-hole gas viscosity is expressed as:

[0068]

[0069] μ′ g =1.55×10 -6 (1.8T + 32)+9.596×10 -3

[0070] Where: μ g is the bottom-hole gas viscosity, with the dimension of Pa·s; μ′ g is the viscosity of the gas under the standard state of the gas temperature that satisfies the correction, with the dimension of Pa·s. The standard state is the environmental condition with a temperature of 273.15K and a pressure of 101.325KPa; P r is the reduced pressure of the gas, dimensionless; T r is the reduced temperature of the gas, dimensionless; T is the gas layer temperature of the target gas layer, with the dimension of K; a 0 ~a 15 are different intermediate calculation quantities, dimensionless; a 0 ~a 15 The values of are respectively a 0 = - 2.46; a 1 =2.97; a 2 = - 0.28; a 3 =0.008; a 4 =0.008; a 5 = - 3.50; a 6 =0.36; a 7 = - 0.01; a 8 = - 0.79; a 9 =1.39; a 10= -0.15; a 11 = 0.004; a 12 = 0.08; a 13 = -0.19; a 14 = 0.02; a 15 = -0.0006.

[0071] As a refinement of the above embodiment, the formula for calculating the bottom-hole gas density is expressed as:

[0072]

[0073] where Z sc is the gas compressibility factor under the standard state, dimensionless; T sc is the formation temperature of the target gas reservoir under the standard state, with the dimension of K; P sc is the formation pressure under the standard state, with the dimension of Pa; Z is the compressibility factor of the gas under the bottom-hole state, dimensionless; b 1 ~b 8 are intermediate calculation quantities. It should be understood that there is no gas in the target well originally. However, since the temperature and pressure of the gas in the target gas reservoir, i.e., the high-pressure gas reservoir, are greater than those in the target well, there is a pressure difference between the two. Therefore, the gas in the target gas reservoir will enter the target well, making there be bottom-hole gas in the target well. Further, the gas components of the gas under the standard state and the bottom-hole gas mentioned above are essentially the same as those of the gas in the target gas reservoir.

[0074] As a refinement of the above embodiment, the formula for calculating the non-Darcy flow coefficient of the target gas reservoir based on the rock permeability and the rock porosity is expressed as:

[0075]

[0076] where K is the rock permeability and φ is the rock porosity.

[0077] This embodiment provides the recorded data of the application of the method proposed in the present disclosure to different target wells. The recorded data more intuitively shows the technical effects of the present disclosure, and the specific content is as follows:

[0078] In this embodiment, when drilling the fourth open hole to 8500 meters, gas invasion occurred in the first target well. The basic parameters of the first target well are shown in the following table:

[0079] Table 2 shows the basic parameters of the first target well

[0080]

[0081] Taking into comprehensive consideration factors including rock physical properties, changes in gas density and viscosity under high temperature and pressure, and fluid properties of the drilling fluid, etc., the gas invasion flow rate results when encountering a high-pressure gas layer during the drilling of the oil and gas well are calculated by the method disclosed in the present disclosure as shown in the table:

[0082] Table 3 shows the gas invasion flow rate results when the target well encounters a high-pressure gas layer obtained by the method disclosed in the present disclosure

[0083]

[0084] Table 4 shows the gas invasion flow rate results when the first target well encounters a high-pressure gas layer obtained by using other models

[0085]

[0086] By comparing the calculation results in the two tables, it can be clearly seen that when calculating the gas invasion flow rate when encountering a high-pressure gas layer, the error between the result of the preset gas invasion flow rate calculation model proposed by the present invention and the actual value is 2.81%. This error value has a significant advantage compared with the calculation results of other existing models. Specifically, the error between the gas invasion flow rate calculated by other models when encountering a high-pressure gas layer and the actual value is as high as 8.62%, which is significantly larger. Generally speaking, the calculation method provided by the present invention is significantly superior to other models in terms of accuracy, can more accurately reflect the change of gas invasion flow rate in the actual working condition, and has higher practical value and application prospect. Therefore, it can be concluded that the calculation result of the present invention not only has effectively improved accuracy, but also is more in line with the actual working condition and has strong engineering application significance.

[0087] In this embodiment, when the second target well is drilled to 8025 meters in the third section, gas invasion occurs in the second target well. The basic parameters of the second target well are shown in the table:

[0088] Table 5 shows the basic parameters of the second target well

[0089]

[0090] Taking into comprehensive consideration factors including rock physical properties, changes in gas density and viscosity under high temperature and pressure, and fluid properties of the drilling fluid, etc., the gas invasion flow rate results when encountering a high-pressure gas layer during the drilling of the oil and gas well are calculated by the method disclosed in the present disclosure as shown in the table:

[0091] Table 6 shows the gas invasion flow rate results when the target well encounters a high-pressure gas layer obtained by the method disclosed in the present disclosure

[0092]

[0093]

[0094] Table 7 shows the gas invasion flow rate results when the target well encounters a high-pressure gas layer obtained using other models.

[0095] Gas invasion flow rate L / min Actual gas invasion flow rate L / min Error analysis % 12.25 14.97 18.16

[0096] By comparing the calculation results in the tables, it can be clearly seen that when calculating the gas invasion flow rate when the well encounters a high-pressure gas layer using the preset gas invasion flow rate calculation model proposed in the present invention, the error between the result and the actual value is 4.21%. This error value has a significant advantage compared to the calculation results of other existing models. Specifically, the error between the gas invasion flow rate calculated by other models when the well encounters a high-pressure gas layer and the actual value is as high as 18.16%, which is significantly larger. Overall, the calculation method provided by the present invention is significantly superior to other models in terms of accuracy, can more accurately reflect the change of gas invasion flow rate in the actual working conditions, and has higher practical value and application prospects. Therefore, it can be concluded that the calculation result of the present invention not only effectively improves the accuracy, but also more conforms to the actual working conditions and has strong engineering application significance.

[0097] In this embodiment, when the fourth section of drilling reached 4375 meters, gas invasion occurred in the third target well. The basic parameters of the third target well are shown in the table:

[0098] Table 8 shows the basic parameters of the third target well

[0099]

[0100] Taking into account factors such as rock physical properties, changes in gas density and viscosity under high temperature and high pressure, and fluid properties of drilling fluid, the gas invasion flow rate results when the oil and gas well encounters a high-pressure gas layer during the drilling process are calculated using the method disclosed in the present invention as shown in the table:

[0101] Table 9 shows the gas invasion flow rate results when the target well encounters a high-pressure gas layer obtained using the method disclosed in the present invention

[0102]

[0103]

[0104] Table 10 shows the gas invasion flow rate results when the target well encounters a high-pressure gas layer obtained using other models

[0105] Gas invasion flow rate L / min Actual gas invasion flow rate L / min Error analysis % 15.54 18.37 15.41

[0106] By comparing the table and the calculation results in the table, it can be clearly seen that when calculating the gas invasion flow rate when drilling through a high-pressure gas layer using the preset gas invasion flow rate calculation model proposed in the present invention, the error between the result and the actual value is 4.57%. This error value has significant advantages compared with the calculation results of other existing models. Specifically, the error between the gas invasion flow rate calculated by other models when drilling through a high-pressure gas layer and the actual value is as high as 15.41%, which is significantly larger. Generally speaking, the calculation method provided by the present invention is significantly superior to other models in terms of accuracy, can more accurately reflect the change of gas invasion flow rate in actual working conditions, and has higher practical value and application prospects. Therefore, it can be concluded that the calculation results of the present invention not only have effectively improved accuracy, but also are more in line with actual working conditions and have strong engineering application significance.

[0107] In summary, the embodiments of the present disclosure can achieve the following effects:

[0108] 1. By comprehensively considering multiple key factors, including the physical properties of the rock (such as rock permeability, rock porosity) and the vertical depth, wellbore diameter, drill pipe outer diameter, and annulus inner wall roughness of the target well, the calculation of the gas invasion flow rate when drilling through a high-pressure gas layer is realized.

[0109] 2. When calculating the gas invasion flow rate using the preset gas invasion flow rate calculation model proposed in the present disclosure, the scale of the required data parameters is small and the structure is simple, thus significantly reducing the time consumption for calculating the gas invasion flow rate. This not only ensures the high accuracy of the result, but also greatly improves the operation speed and efficiency. This design provides strong support for the efficient decision-making of subsequent overflow disposal, provides a more scientific basis for the optimization decision-making and operation in actual engineering, and greatly promotes the progress of related technologies and the improvement of application value.

[0110] Corresponding to the above calculation method of the gas invasion flow rate when drilling through a high-pressure gas layer, the present invention also proposes a calculation device for the gas invasion flow rate when drilling through a high-pressure gas layer. Since the device embodiment of the present invention corresponds to the above method embodiment, for the details not disclosed in the device embodiment, reference can be made to the above method embodiment, and the present invention will not be elaborated herein.

[0111] Figure 2 The structural schematic diagram of a calculation device for the gas invasion flow rate when drilling through a high-pressure gas layer provided by the embodiments of the present disclosure is as Figure 2 shown, including:

[0112] An acquisition unit 21, configured to detect the physical properties of the rock sample, obtain the rock permeability and rock porosity corresponding to the rock sample, and acquire the vertical depth, wellbore diameter, drill pipe outer diameter, and annulus inner wall roughness of the target well, where the rock sample is obtained by sampling the rock of the target gas layer in the target well;

[0113] The first calculation unit 22 is configured to calculate the non-Darcy flow coefficient of the target gas reservoir based on the rock permeability and the rock porosity, and calculate the bottom-hole pressure of the target well during normal drilling based on the vertical depth, the wellbore diameter, the outer diameter of the drill pipe, and the roughness of the inner wall of the annulus.

[0114] The second calculation unit 23 is configured to calculate the gas invasion flow rate during the drilling of the target well based on a preset gas invasion flow rate calculation model in combination with the non-Darcy flow coefficient and the bottom-hole pressure.

[0115] The calculation method of gas invasion flow rate when encountering a high-pressure gas reservoir provided by the present disclosure detects the physical properties of a rock sample to obtain the rock permeability and rock porosity corresponding to the rock sample, and obtains the vertical depth, wellbore diameter, outer diameter of the drill pipe, and roughness of the inner wall of the annulus of the target well. The rock sample is obtained by sampling the rock of the target gas reservoir in the target well; calculates the non-Darcy flow coefficient of the target gas reservoir based on the rock permeability and the rock porosity, and calculates the bottom-hole pressure of the target well during normal drilling based on the vertical depth, the wellbore diameter, the outer diameter of the drill pipe, and the roughness of the inner wall of the annulus; calculates the gas invasion flow rate during the drilling of the target well based on a preset gas invasion flow rate calculation model in combination with the non-Darcy flow coefficient and the bottom-hole pressure. Compared with the related art, by comprehensively considering multiple key factors, including the physical properties of the rock (such as rock permeability, rock porosity) and the vertical depth, wellbore diameter, outer diameter of the drill pipe, and roughness of the inner wall of the annulus of the target well, the calculation of the gas invasion flow rate when encountering a high-pressure gas reservoir is realized.

[0116] It should be noted that the foregoing explanation of the method embodiment also applies to the device of this embodiment, and the principle is the same, so it is not limited in this embodiment.

[0117] According to an embodiment of the present disclosure, the present disclosure also provides an electronic device, a readable storage medium, and a computer program product.

[0118] Figure 3 FIG. shows a schematic block diagram of an exemplary electronic device 300 that can be used to implement the embodiments of the present disclosure. The electronic device is intended to represent various forms of digital computers, such as, for example, a laptop computer, a desktop computer, a workbench, a personal digital assistant, a server, a blade server, a mainframe computer, and other suitable computers. The electronic device can also represent various forms of mobile devices, such as, for example, a personal digital processor, a cellular phone, a smart phone, a wearable device, and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely exemplary and are not intended to limit the implementation of the present disclosure described and / or claimed herein.

[0119] AsFigure 3 As shown, device 300 includes a computing unit 301, which can execute various appropriate actions and processes according to computer programs stored in a ROM (Read-Only Memory) 302 or computer programs loaded from a storage unit 308 into a RAM (Random Access Memory) 303. In the RAM 303, various programs and data required for the operation of device 300 can also be stored. The computing unit 301, the ROM 302, and the RAM 303 are connected to each other via a bus 304. An I / O (Input / Output) interface 305 is also connected to the bus 304.

[0120] Multiple components in device 300 are connected to the I / O interface 305, including: an input unit 306, such as a keyboard, a mouse, etc.; an output unit 307, such as various types of displays, speakers, etc.; a storage unit 308, such as a magnetic disk, an optical disc, etc.; and a communication unit 309, such as a network card, a modem, a wireless communication transceiver, etc. The communication unit 309 allows device 300 to exchange information / data with other devices via a computer network such as the Internet and / or various telecommunication networks.

[0121] The computing unit 301 can be various general-purpose and / or special-purpose processing components with processing and computing capabilities. Some examples of the computing unit 301 include but are not limited to a CPU (Central Processing Unit), a GPU (Graphic Processing Units), various dedicated AI (Artificial Intelligence) computing chips, various computing units running machine learning model algorithms, a DSP (Digital Signal Processor), and any appropriate processor, controller, microcontroller, etc. The computing unit 301 executes the various methods and processes described above, such as the method for calculating the gas invasion flow rate when drilling through a high-pressure gas layer. For example, in some embodiments, the method for calculating the gas invasion flow rate when drilling through a high-pressure gas layer can be implemented as a computer software program, which is tangibly contained in a machine-readable medium, such as the storage unit 308. In some embodiments, part or all of the computer program can be loaded and / or installed onto device 300 via the ROM 302 and / or the communication unit 309. When the computer program is loaded into the RAM 303 and executed by the computing unit 301, one or more steps of the method described above can be executed. Alternatively, in other embodiments, the computing unit 301 can be configured to execute the aforementioned method for calculating the gas invasion flow rate when drilling through a high-pressure gas layer in any other appropriate manner (e.g., by means of firmware).

[0122] The various embodiments of the systems and techniques described above in this specification can be implemented in digital electronic circuitry, integrated circuit systems, FPGAs (Field Programmable Gate Arrays), ASICs (Application-Specific Integrated Circuits), ASSPs (Application Specific Standard Products), SoCs (System On Chip), CPLDs (Complex Programmable Logic Devices), computer hardware, firmware, software, and / or combinations thereof. These various embodiments can include: being implemented in one or more computer programs that are executable and / or interpretable on a programmable system including at least one programmable processor, which can be a special-purpose or general-purpose programmable processor that receives data and instructions from a storage system, at least one input device, and at least one output device, and transmits the data and instructions to the storage system, the at least one input device, and the at least one output device.

[0123] The program code for implementing the methods of the present disclosure can be written in any combination of one or more programming languages. These program codes can be provided to a processor or controller of a general-purpose computer, a special-purpose computer, or other programmable data processing device, such that when the program codes are executed by the processor or controller, the functions / operations specified in the flowchart and / or block diagram are implemented. The program code can be executed entirely on the machine, partially on the machine, as a stand-alone software package partially on the machine and partially on a remote machine, or entirely on a remote machine or server.

[0124] In the context of this disclosure, a machine-readable medium can be a tangible medium that can contain or store a program for use by or in connection with an instruction execution system, apparatus, or device. A machine-readable medium can be a machine-readable signal medium or a machine-readable storage medium. A machine-readable medium can include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, apparatus, or devices, or any suitable combination of the foregoing. More specific examples of a machine-readable storage medium would include an electrical connection based on one or more wires, a portable computer disk, a hard disk, a RAM, a ROM, an EPROM (Electrically Programmable Read-Only-Memory), or a flash memory, an optical fiber, a CD-ROM (Compact Disc Read-Only Memory), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

[0125] In order to provide interaction with a user, the systems and techniques described herein can be implemented on a computer having: a display device for displaying information to the user (e.g., a CRT (Cathode-Ray Tube) or an LCD (Liquid Crystal Display) monitor); and a keyboard and a pointing device (e.g., a mouse or a trackball) through which the user can provide input to the computer. Other kinds of devices can also be used to provide interaction with the user; for example, the feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form (including acoustic input, voice input, or tactile input).

[0126] The systems and techniques described herein can be implemented in a computing system including backend components (e.g., as a data server), or a computing system including middleware components (e.g., an application server), or a computing system including frontend components (e.g., a user computer having a graphical user interface or a web browser through which the user can interact with an implementation of the systems and techniques described herein), or a computing system including any combination of such backend components, middleware components, or frontend components. The components of the system can be interconnected to each other by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include: a LAN (Local Area Network), a WAN (Wide Area Network), the Internet, and a blockchain network.

[0127] A computer system may include a client and a server. The client and the server are generally far from each other and usually interact through a communication network. The client-server relationship is generated by computer programs running on the respective computers and having a client-server relationship with each other. The server can be a cloud server, also known as a cloud computing server or a cloud host, which is a host product in the cloud computing service system, solving the defects of difficult management and weak business scalability existing in traditional physical hosts and VPS services ("Virtual Private Server", or simply "VPS" for short). The server can also be a server of a distributed system or a server combined with blockchain.

[0128] Among them, it should be noted that artificial intelligence is a discipline that studies how to make a computer simulate certain thinking processes and intelligent behaviors of humans (such as learning, reasoning, thinking, planning, etc.), and there are both hardware-level technologies and software-level technologies. Artificial intelligence hardware technologies generally include technologies such as sensors, dedicated artificial intelligence chips, cloud computing, distributed storage, and big data processing; artificial intelligence software technologies mainly include several major directions such as computer vision technology, speech recognition technology, natural language processing technology, and machine learning / deep learning, big data processing technology, and knowledge graph technology.

Claims

1. A method for calculating gas intrusion flow when drilling into a high-pressure gas layer, characterized in that: include: Detecting the physical properties of rock samples, obtaining the rock permeability and rock porosity corresponding to the rock samples, and acquiring the vertical depth, borehole diameter, drill pipe outer diameter, and annulus inner wall roughness of the target well, wherein the rock samples are obtained by sampling the rocks of the target gas layer in the target well; Calculating the non-Darcy flow coefficient of the target gas layer based on the rock permeability and the rock porosity and calculating the bottom hole pressure of the target well during normal drilling based on the vertical depth, the wellbore diameter, the drill pipe outer diameter and the annulus inner wall roughness; The gas invasion flow rate during the drilling of the target well is calculated based on a preset gas invasion flow rate calculation model combined with the non-Darcy flow coefficient and the bottom hole pressure.

2. The method according to claim 1, characterized in that The formula of the preset gas intrusion flow calculation model is expressed as: Among them, P f is the gas layer pressure of the target gas layer; P w is the bottom hole pressure of the target well during normal drilling; μ g is the bottom hole gas viscosity of the target well; K is the rock permeability; q g is the gas intrusion flow rate; r w is the wellbore radius of the target well; h is the gas layer thickness of the target gas layer; δ is the non-Darcy flow coefficient; r f is the air layer radius of the target air layer; ρ g is the bottom hole gas density.

3. The method according to claim 2, characterized in that The formula for calculating the bottom hole pressure of the target well during normal drilling based on the vertical depth, the wellbore diameter, the outer diameter of the drill pipe and the roughness of the inner wall of the annulus is expressed as: Among them, ρ l is the drilling fluid density of the target well; g is the gravity acceleration of the target well; H is the vertical depth; ξ is the friction factor; Q is the drilling fluid displacement of the target well; D is the wellbore diameter; d is the outer diameter of the drill pipe; ε is the roughness of the inner wall of the annulus; μ l is the drilling fluid viscosity of the target well; L is the annulus length, a is the first intermediate calculation quantity; b is the second intermediate calculation quantity.

4. The method according to claim 3, characterized in that The formula for calculating the drilling fluid density is expressed as: r l =ρ0+c0P w +c1T+c2P w T+c3T 2 +c4T 3 Among them, ρ0 is the density of the drilling fluid under normal pressure; c0 is the first relevant parameter value related to ρ0, c1 is the second relevant parameter value related to ρ0, c2 is the second relevant parameter value related to ρ0, c3 is the second relevant parameter value related to ρ0, c4 is the second relevant parameter value related to ρ0, and the values ​​of c0, c1, c2, c3, and c4 are all determined based on the composition and type of the drilling fluid.

5. The method according to claim 4, characterized in that The formula for calculating the bottom hole gas viscosity is expressed as: m′ g =1.55×10 -6 (1.8T+32)+9.596×10 -3 Where: μ g is the bottom hole gas viscosity; μ′ g The viscosity of the gas under standard conditions that meet the calibration gas temperature; P r is the relative pressure of the gas; T r is the relative temperature of the gas; T is the gas layer temperature of the target gas layer; a0~a 15 For different intermediate calculation quantities.

6. The method according to claim 5, characterized in that The formula for calculating the bottom hole gas density is expressed as: Among them, Z sc is the gas compression factor under the standard state; T sc is the gas layer temperature of the target gas layer under the standard state; P sc is the gas layer pressure under the standard state; Z is the gas compression factor under the bottom hole state; b1~b8 are intermediate calculation quantities.

7. The method according to claim 6, characterized in that The formula for calculating the non-Darcy flow coefficient of the target gas layer based on the rock permeability and the rock porosity is expressed as: Wherein, K is the rock permeability, and φ is the rock porosity.

8. A device for calculating gas intrusion flow when drilling into a high-pressure gas layer, characterized in that: include: an acquisition unit, used to detect the physical properties of the rock sample, obtain the rock permeability and rock porosity corresponding to the rock sample, and obtain the vertical depth, wellbore diameter, drill pipe outer diameter and annulus inner wall roughness of the target well, wherein the rock sample is obtained by sampling the rock of the target gas layer in the target well; A first calculation unit is used to calculate the non-Darcy flow coefficient of the target gas layer based on the rock permeability and the rock porosity, and calculate the bottom hole pressure of the target well during normal drilling based on the vertical depth, the wellbore diameter, the drill pipe outer diameter and the annulus inner wall roughness; The second calculation unit is used to calculate the gas invasion flow rate during the drilling of the target well based on a preset gas invasion flow rate calculation model combined with the non-Darcy flow coefficient and the bottom hole pressure.

9. An electronic device, characterized in that: include: at least one processor; as well as a memory communicatively connected to the at least one processor; wherein, The memory stores instructions that can be executed by the at least one processor, and the instructions are executed by the at least one processor to enable the at least one processor to perform the method according to any one of claims 1 to 7.

10. A non-transitory computer-readable storage medium storing computer instructions, characterized in that: The computer instructions are used to cause the computer to execute the method according to any one of claims 1-7.

Images