Method, device and equipment for determining water breakthrough time of bottom water gas reservoir and medium

Abstract

The invention discloses a water breakthrough time determination method, device and equipment of a bottom water gas reservoir and a medium. The method comprises the following steps: determining the bottom daily gas production rate of a perforation section of a target bottom water gas reservoir; the water invasion process of the target bottom water gas reservoir is divided into three stages, and the water breakthrough time of each stage is determined based on the daily gas production rate of the bottom of the perforation interval; wherein the three stages comprise the first stage, the second stage and the third stage, the first stage is used for representing that the gas-water interface rises to the interlayer from the position below the interlayer, the second stage is used for representing that the gas-water interface moves to the center of the interlayer from the edge of the interlayer, and the third stage is used for representing that the gas-water interface rises to the bottom of the perforation section from the center of the interlayer; and determining the water breakthrough time of the target bottom water gas reservoir according to the sum of the water breakthrough time of each stage. According to the scheme, the bottom water gas reservoir is divided into three stages, the water breakthrough time is calculated in stages in combination with the daily gas production rate change in the actual production process, and the accuracy and applicability of prediction of the water breakthrough time of the bottom water gas reservoir are effectively improved.

Description

Technical Field

[0001] The present invention relates to the technical field of bottom water gas reservoir development, and in particular to a method, device, equipment and medium for determining water breakthrough time of a bottom water gas reservoir. Background Art

[0002] In recent years, my country has made remarkable achievements in natural gas exploration and development. A series of bottom water gas reservoirs have been developed in the main domestic natural gas producing areas, showing great development potential. Most bottom water gas reservoirs are affected by water invasion to varying degrees as the development stage changes. Water in the gas reservoir will reduce the production capacity of the gas well, seriously affecting the economic benefits and development process of the gas field. Some bottom water gas reservoirs are seriously heterogeneous, and there are low-permeability or impermeable interlayers in the reservoir. The interlayers will cause the bottom water to bypass and extend the route of the bottom water cone, thereby greatly extending the water breakthrough time and increasing the gas production of the gas well.

[0003] In the related art, a steady-state model is used to calculate the water breakthrough time of bottom water gas reservoirs. In the steady-state model, certain parameters (such as daily gas production) are regarded as quantitative values ​​that do not change with time. However, in the actual production process, the above parameters will change with time. Therefore, it is difficult for the steady-state model to accurately predict the water breakthrough time, and the applicability of the steady-state model is poor, and it cannot provide an effective basis for reasonable development. Summary of the invention

[0004] The present invention provides a method, device, equipment and medium for determining the water breakthrough time of a bottom water gas reservoir, which divides the bottom water gas reservoir into three stages, and calculates the water breakthrough time in stages based on the change of daily gas production in the actual production process, thereby effectively improving the accuracy and applicability of the water breakthrough time prediction of the bottom water gas reservoir.

[0005] According to one aspect of the present invention, a method for determining water breakthrough time of a bottom water gas reservoir is provided, the method comprising:

[0006] Determine the daily gas production at the bottom of the perforation section of the target bottom water gas reservoir; wherein the target bottom water gas reservoir includes interlayers;

[0007] The water invasion process of the target bottom water gas reservoir is divided into three stages, and the water breakthrough time of each stage is determined based on the daily gas production at the bottom of the perforation section; wherein the three stages include the first stage, the second stage and the third stage, the first stage is used to characterize the rise of the gas-water interface from below the interlayer to the interlayer, the second stage is used to characterize the movement of the gas-water interface from the edge of the interlayer to the center of the interlayer, and the third stage is used to characterize the rise of the gas-water interface from the center of the interlayer to the bottom of the perforation section;

[0008] The water breakthrough time of the target bottom water gas reservoir is determined according to the sum of the water breakthrough times of each stage.

[0009] According to another aspect of the present invention, a device for determining water breakthrough time of a bottom water gas reservoir is provided, the device comprising:

[0010] A perforation section bottom production determination module is used to determine the daily gas production at the bottom of the perforation section of a target bottom water gas reservoir; wherein the target bottom water gas reservoir includes an interlayer;

[0011] A three-stage water breakthrough time determination module is used to divide the water invasion process of the target bottom water gas reservoir into three stages, and determine the water breakthrough time of each stage based on the daily gas production at the bottom of the perforation section; wherein the three stages include a first stage, a second stage and a third stage, the first stage is used to characterize the rise of the gas-water interface from below the interlayer to the interlayer, the second stage is used to characterize the movement of the gas-water interface from the edge of the interlayer to the center of the interlayer, and the third stage is used to characterize the rise of the gas-water interface from the center of the interlayer to the bottom of the perforation section;

[0012] The bottom water gas reservoir water breakthrough time determination module is used to determine the water breakthrough time of the target bottom water gas reservoir according to the sum of the water breakthrough time of each stage.

[0013] According to another aspect of the present invention, an electronic device is provided, the electronic device comprising:

[0014] at least one processor; and,

[0015] a memory communicatively connected to the at least one processor; wherein,

[0016] The memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can execute the method for determining the water breakthrough time of a bottom water gas reservoir described in any embodiment of the present invention.

[0017] According to another aspect of the present invention, a computer-readable storage medium is provided, wherein the computer-readable storage medium stores computer instructions, and the computer instructions are used to implement the method for determining water breakthrough time of a bottom water gas reservoir described in any embodiment of the present invention when executed by a processor.

[0018] The technical solution of the embodiment of the present invention determines the daily gas production at the bottom of the perforation section of the target bottom water gas reservoir; wherein the target bottom water gas reservoir includes an interlayer; the water invasion process of the target bottom water gas reservoir is divided into three stages, and the water breakthrough time of each stage is determined based on the daily gas production at the bottom of the perforation section; wherein the three stages include the first stage, the second stage and the third stage, the first stage is used to characterize the rise of the gas-water interface from the bottom of the interlayer to the interlayer, the second stage is used to characterize the movement of the gas-water interface from the edge of the interlayer to the center of the interlayer, and the third stage is used to characterize the rise of the gas-water interface from the center of the interlayer to the bottom of the perforation section; the water breakthrough time of the target bottom water gas reservoir is determined according to the sum of the water breakthrough time of each stage. This technical solution, by dividing the bottom water gas reservoir into three stages and calculating the water breakthrough time in stages in combination with the daily gas production changes in the actual production process, can effectively improve the accuracy and applicability of the prediction of the water breakthrough time of the bottom water gas reservoir, so as to provide theoretical support for the development of the bottom water gas reservoir, help optimize the development plan, improve the development efficiency, and promote technological progress and industrial development in the field of gas reservoir development.

[0019] It should be understood that the contents described in this section are not intended to identify the key or important features of the embodiments of the present invention, nor are they intended to limit the scope of the present invention. Other features of the present invention will become easily understood through the following description. BRIEF DESCRIPTION OF THE DRAWINGS

[0020] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0021] Figure 1 This is a flow chart of a method for determining water breakthrough time of a bottom water gas reservoir provided according to Embodiment 1 of the present invention;

[0022] Figure 2 is a schematic diagram of water intrusion stage division provided according to Embodiment 1 of the present invention;

[0023] Figure 3 is a schematic diagram of a pseudo-stabilization principle provided according to the first embodiment of the present invention;

[0024] Figure 4 This is a flow chart of a method for determining water breakthrough time of a bottom water gas reservoir provided in accordance with Embodiment 2 of the present invention;

[0025] Figure 5 2 is a schematic structural diagram of a device for determining water breakthrough time of a bottom water gas reservoir provided according to a third embodiment of the present invention;

[0026] Figure 6The present invention is a schematic diagram of the structure of an electronic device for implementing a method for determining water breakthrough time of a bottom water gas reservoir according to an embodiment of the present invention. DETAILED DESCRIPTION

[0027] In order to enable those skilled in the art to better understand the scheme of the present invention, the technical scheme in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of the present invention.

[0028] It should be noted that the terms "first", "second", "target", etc. in the specification and claims of the present invention and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchanged where appropriate, so that the embodiments of the present invention described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0029] Embodiment 1

[0030] Figure 1 This is a flow chart of a method for determining the breakthrough time of a bottom water gas reservoir provided in the first embodiment of the present invention. This embodiment can be applied to the case of improving the prediction accuracy of the breakthrough time of a bottom water gas reservoir. The method can be executed by a device for determining the breakthrough time of a bottom water gas reservoir. The device for determining the breakthrough time of a bottom water gas reservoir can be implemented in the form of hardware and / or software. The device for determining the breakthrough time of a bottom water gas reservoir can be configured in an electronic device with data processing capabilities. Figure 1 As shown, the method includes:

[0031] S110, determining the daily gas production at the bottom of the perforation section of the target bottom water gas reservoir; wherein the target bottom water gas reservoir includes interlayers.

[0032] In this embodiment, first, according to the different flow states of the target bottom water gas reservoir, the gas production section of the target bottom water gas reservoir is divided into two parts: the perforation section and the bottom of the perforation section. Among them, the perforation section presents a non-Darcy plane radial flow, and the bottom of the perforation section includes a near-wellbore zone at the bottom of the perforation section and a far-wellbore zone at the bottom of the perforation section. The near-wellbore zone at the bottom of the perforation section presents a non-Darcy hemispherical center-facing flow, and the far-wellbore zone at the bottom of the perforation section presents a Darcy plane radial flow. Then, the daily gas production at the bottom of the perforation section of the target bottom water gas reservoir can be determined based on the following formula:

[0033]

[0034] in,

[0035] In the above formula, q is the daily gas production of the gas well; q2 is the production at the bottom of the perforated section; D1 is the inertia coefficient or turbulence coefficient of the non-Darcy plane radial flow; D2 is the inertia or turbulence coefficient of the non-Darcy hemispherical center-facing flow; h p is the length of the perforation section; r1 is the range of the non-Darcy hemisphere facing the center flow; r e is the supply radius; r w is the wellbore radius; S is the skin factor, dimensionless; β is the velocity coefficient, which describes the influencing parameter of turbulence and is related to permeability, β = 7.644 × 10 10 / K 1.5 , unit is m -3 ; γ g is the relative density of gas; K is the permeability; For viscosity.

[0036] S120, dividing the water invasion process of the target bottom water gas reservoir into three stages, and determining the water breakthrough time of each stage based on the daily gas production at the bottom of the perforated section.

[0037] Among them, the three stages include the first stage, the second stage and the third stage. The first stage is used to characterize the rise of the gas-water interface from the bottom of the interlayer to the interlayer, the second stage is used to characterize the movement of the gas-water interface from the edge of the interlayer to the center of the interlayer, and the third stage is used to characterize the rise of the gas-water interface from the center of the interlayer to the bottom of the perforation section. Figure 2 FIG. 1 is a schematic diagram of water intrusion stage division provided in Embodiment 1 of the present invention. Figure 2 As shown, stage I (AB) represents the first stage, stage II (BC) represents the second stage, and stage III (CD) represents the third stage; h b is the height from the interlayer to the initial air-water interface, h c is the height from the interlayer to the bottom of the perforated section, h p is the length of the perforation section, h is the thickness of the gas layer, h = h b +h c +h p , rb is the distance from the edge to the center of the interlayer, r e is the supply radius.

[0038] In this embodiment, the water invasion process of the target bottom water gas reservoir is divided into three stages, and the water invasion process in each stage is subjected to quasi-steady-state processing to determine the water breakthrough time of each stage. The quasi-steady-state processing can be understood as dividing the water invasion process of each stage into multiple time periods, and each time period is regarded as a steady state, that is, the daily gas production at the bottom of the perforated section in each time period is regarded as a fixed value, which is determined based on the average daily gas production at the bottom of the perforated section at each time in the corresponding time period. Figure 3 A schematic diagram of a pseudo-steady-state principle provided in the first embodiment of the present invention. n Indicates the water invasion height (i.e., the height of the air-water interface rise) at different time periods, and the dotted curves indicate the air-water interface corresponding to different time periods. Above and below the air-water interface are gas and water, respectively. Figure 3 As shown, the target bottom water gas reservoir includes a directional well and an interlayer. After the target bottom water gas reservoir produces gas, the gas will enter the ground from underground through the directional well, resulting in a reduction of underground gas, causing water to move upward, thereby causing the gas-water interface to change.

[0039] In this embodiment, optionally, the water breakthrough time of each stage is determined based on the daily gas production at the bottom of the perforation section, including: dividing the water invasion process of the first stage into multiple time periods based on a preset time interval; wherein the daily gas production at the bottom of the perforation section in each time period is regarded as a fixed value; determining the rising heights of the gas-water interface corresponding to adjacent time periods in the first stage based on the daily gas production at the bottom of the perforation section; judging whether the gas-water interface reaches the boundary of the first stage according to the rising heights of the gas-water interface corresponding to adjacent time periods and the height from the interlayer to the initial gas-water interface; if so, determining the water breakthrough time of the first stage based on the daily gas production at the bottom of the perforation section.

[0040] In this embodiment, the water invasion process of the first stage is first divided into multiple time periods based on a preset time interval, and the average daily gas production at the bottom of the perforation section at each moment in each time period is calculated as the daily gas production at the bottom of the perforation section in the corresponding time period (fixed value). Among them, the preset time interval may refer to a time interval pre-set according to actual needs, and the time interval may be the same time interval or a different time interval. Exemplarily, according to the actual production law, the time interval in the early stage of production may be set to be smaller (such as 5 days), and the time interval may be continuously increased as production continues (such as 15-30 days). Then, based on the daily gas production at the bottom of the perforation section, the gas-water interface rise height corresponding to the adjacent time periods in the first stage is determined. It should be noted that in the first stage, since the gas-water interface is far away from the bottom of the perforation section, this process can be regarded as a Darcy radial flow process (corresponding to the Darcy plane radial flow in the far well area at the bottom of the perforation section). Optionally, determining the gas-water interface rise height corresponding to the adjacent time periods in the first stage based on the daily gas production at the bottom of the perforation section includes: determining the gas-water interface rise height corresponding to the adjacent time periods in the first stage based on the first iterative formula, and the first iterative formula is expressed as follows:

[0041]

[0042] Where l = h c +h b -z, z is the height of the air-water interface, h c is the height from the interlayer to the bottom of the perforated section, h b is the height from the interlayer to the initial air-water interface, i and i+1 represent the previous time period and the next time period in the adjacent time period, respectively, r b is the distance from the edge of the interlayer to the center, q2 is the daily gas production at the bottom of the perforated section, φ is the porosity, S wi is the bound water saturation, S gr is the residual gas saturation, and t is the time.

[0043] After determining the rising height of the gas-water interface corresponding to the adjacent time periods in the first stage, it is possible to determine whether the gas-water interface has reached the boundary of the first stage based on the rising height of the gas-water interface corresponding to the adjacent time periods and the height from the interlayer to the initial gas-water interface. If so, it indicates that the first stage has been completed, and the water breakthrough time of the first stage can be determined based on the daily gas production at the bottom of the perforated section; otherwise, it indicates that the first stage has not been completed, and the above-mentioned calculation of the rising height of the gas-water interface and the first stage boundary judgment process need to be repeated until the gas-water interface reaches the first stage boundary.

[0044] In this embodiment, optionally, judging whether the gas-water interface has reached the boundary of the first stage according to the rising height of the gas-water interface corresponding to the adjacent time periods and the height from the interlayer to the initial gas-water interface includes: if the rising height of the gas-water interface corresponding to the previous time period in the adjacent time periods is less than the height from the interlayer to the initial gas-water interface, and the rising height of the gas-water interface corresponding to the subsequent time period in the adjacent time periods is greater than the height from the interlayer to the initial gas-water interface, then determining that the gas-water interface has reached the boundary of the first stage; correspondingly, determining the water breakthrough time of the first stage based on the daily gas production at the bottom of the perforation section includes: determining the water breakthrough time of the first stage based on the first time formula, and the first time formula is expressed as follows:

[0045]

[0046] Among them, t AB Indicates the water breakthrough time of the first stage.

[0047] For example, refer to Figure 2 , assuming that the adjacent time periods are i and i+1, i is the earlier time period and i+1 is the later time period. i <h b And z i+1 >h b It indicates that the previous time period i is in the first stage, and the subsequent time period i+1 is in the second stage. Therefore, it can be determined that the air-water interface reaches the boundary of the first stage, and the moment when the air-water interface reaches the boundary of the first stage is between i and i+1. In this case, the water breakthrough time of the first stage can be calculated based on the above-mentioned first time formula using the relevant parameters of the previous time period i.

[0048] In this embodiment, optionally, the water breakthrough time of each stage is determined based on the daily gas production at the bottom of the perforation section, including: dividing the water invasion process of the second stage into multiple time periods based on a preset time interval; determining the distance from the gas-water interface to the edge of the interlayer corresponding to adjacent time periods in the second stage based on the daily gas production at the bottom of the perforation section; judging whether the gas-water interface reaches the boundary of the second stage according to the distance from the gas-water interface to the edge of the interlayer and the distance from the edge of the interlayer to the center corresponding to the adjacent time periods; if so, determining the water breakthrough time of the second stage based on the daily gas production at the bottom of the perforation section.

[0049] It should be noted that in the second stage, after the gas-water interface rises to the interlayer, it needs to flow radially to the center of the interlayer and then vertically to the bottom of the perforation section. Since the gas-water interface is also far away from the bottom of the perforation section in this stage, it can also be regarded as a Darcy hemispheric flow process, but the flow direction is different from that in the first stage.

[0050] In this embodiment, optionally, determining the distance from the gas-water interface to the edge of the interlayer corresponding to the adjacent time periods in the second stage based on the daily gas production at the bottom of the perforation section includes: determining the distance from the gas-water interface to the edge of the interlayer corresponding to the adjacent time periods in the second stage based on a second iterative formula, and the second iterative formula is expressed as follows:

[0051]

[0052] Among them, h c is the height from the interlayer to the bottom of the perforation section, i and i+1 represent the previous time period and the next time period in the adjacent time period, q2 is the daily gas production at the bottom of the perforation section, φ is the porosity, S wi is the bound water saturation, S gr is the residual gas saturation, t is the time, and x is the distance from the gas-water interface to the edge of the interlayer in the second stage.

[0053] After determining the distance to the edge of the interlayer corresponding to the adjacent time period in the second stage, it is possible to determine whether the gas-water interface has reached the boundary of the second stage based on the distance from the gas-water interface to the edge of the interlayer and the distance from the edge of the interlayer to the center corresponding to the adjacent time period. If so, it indicates that the second stage has been completed, and the water breakthrough time of the second stage can be determined based on the daily gas production at the bottom of the perforated section; otherwise, it indicates that the second stage has not been completed, and the above-mentioned distance calculation from the gas-water interface to the edge of the interlayer and the second stage boundary judgment process need to be repeated until the gas-water interface reaches the boundary of the second stage.

[0054] In this embodiment, optionally, judging whether the gas-water interface has reached the boundary of the second stage according to the distance from the gas-water interface to the edge of the interlayer and the distance from the edge of the interlayer to the center corresponding to the adjacent time periods, including: if the distance from the gas-water interface to the edge of the interlayer corresponding to the previous time period in the adjacent time periods is less than the distance from the edge of the interlayer to the center, and the distance from the gas-water interface to the edge of the interlayer corresponding to the subsequent time period in the adjacent time periods is greater than the distance from the edge of the interlayer to the center, then determining that the gas-water interface has reached the boundary of the second stage; correspondingly, determining the water breakthrough time of the second stage based on the daily gas production at the bottom of the perforation section, including: determining the water breakthrough time of the second stage based on the second time formula, the second time formula is expressed as follows:

[0055]

[0056] Among them, t BC represents the water breakthrough time of the second stage, r b Indicates the distance from the edge to the center of the interlayer.

[0057] For example, refer to Figure 2 , assuming that the adjacent time periods are i and i+1, i is the earlier time period and i+1 is the later time period. i<r b And x i+1 >r b It indicates that the previous time period i is in the second stage, and the subsequent time period i+1 is in the third stage. Therefore, it can be determined that the air-water interface reaches the boundary of the second stage, and the moment when the air-water interface reaches the boundary of the second stage is between i and i+1. In this case, the water breakthrough time of the second stage can be calculated based on the above-mentioned second time formula using the relevant parameters of the previous time period i.

[0058] In this embodiment, optionally, the water breakthrough time of each stage is determined based on the daily gas production at the bottom of the perforation section, including: dividing the water invasion process of the third stage into multiple time periods based on a preset time interval; wherein the daily gas production at the bottom of the perforation section in each time period is regarded as a fixed value; determining the distance from the gas-water interface corresponding to the adjacent time periods in the third stage to the bottom of the perforation section based on the daily gas production at the bottom of the perforation section; judging whether the gas-water interface reaches the boundary of the third stage according to the distance from the gas-water interface corresponding to the adjacent time periods to the bottom of the perforation section; if so, determining the water breakthrough time of the third stage based on the daily gas production at the bottom of the perforation section.

[0059] It should be noted that in the third stage, after the gas-water interface rises above the center of the interlayer, the gas-water interface will rise further and move toward the bottom of the perforation section, which is closer to the bottom of the perforation section. Therefore, the flow state in this stage can be regarded as high-speed non-Darcy flow (corresponding to the non-Darcy hemispherical center-facing flow in the near-wellbore area at the bottom of the perforation section).

[0060] In this embodiment, optionally, determining the distance from the gas-water interface corresponding to the adjacent time periods in the third stage to the bottom of the perforation section based on the daily gas production at the bottom of the perforation section includes: determining the distance from the gas-water interface corresponding to the adjacent time periods in the third stage to the bottom of the perforation section based on a third iterative formula, and the third iterative formula is expressed as follows:

[0061]

[0062] in, β is the velocity coefficient, ρ g is the gas density, q2 is the daily gas production at the bottom of the perforation section, φ is the porosity, r is the distance from the gas-water interface to the bottom of the perforation section in the third stage, i-1 and i represent the previous time period and the next time period in the adjacent time period, respectively, S wi is the bound water saturation, S gr is the residual gas saturation, and t is the time.

[0063] After determining the distance from the gas-water interface to the bottom of the perforated section corresponding to the adjacent time periods in the third stage, it is possible to determine whether the gas-water interface has reached the boundary of the third stage based on the distance from the gas-water interface to the bottom of the perforated section corresponding to the adjacent time periods. If so, it indicates that the third stage has been completed, and the water breakthrough time of the third stage can be determined based on the daily gas production at the bottom of the perforated section; otherwise, it indicates that the third stage has not been completed, and the above-mentioned calculation of the distance from the gas-water interface to the bottom of the perforated section and the third stage boundary judgment process need to be repeated until the gas-water interface reaches the boundary of the third stage.

[0064] In this embodiment, optionally, judging whether the gas-water interface has reached the third stage boundary according to the distance from the gas-water interface to the bottom of the perforation section corresponding to the adjacent time periods includes: if the distance from the gas-water interface to the bottom of the perforation section corresponding to the previous time period in the adjacent time periods is greater than zero, and the distance from the gas-water interface to the bottom of the perforation section corresponding to the subsequent time period in the adjacent time periods is less than zero, then determining that the gas-water interface has reached the third stage boundary; correspondingly, determining the water breakthrough time of the third stage based on the daily gas production at the bottom of the perforation section includes: determining the water breakthrough time of the third stage based on the third time formula, and the third time formula is expressed as follows:

[0065]

[0066] Among them, t CD Indicates the water breakthrough time of the third stage.

[0067] For example, refer to Figure 2 , assuming that the adjacent time periods are i-1 and i, i-1 is the earlier time period and i is the later time period. i-1 >0 and r i When <0, it indicates that the previous time period i-1 is in the third stage, and the gas-water interface corresponding to the subsequent time period i is already at a position above the bottom of the perforation section, so it can be determined that the gas-water interface reaches the boundary of the third stage, and the moment when the gas-water interface reaches the boundary of the third stage is between i-1 and i. In this case, the water breakthrough time of the third stage can be calculated based on the above third time formula using the relevant parameters of the previous time period i-1. It should be noted that the distance from the gas-water interface to the bottom of the perforation section described here is less than zero only for the convenience of position representation, that is, it means that the gas-water interface is at a position above the bottom of the perforation section.

[0068] S130: determining the water breakthrough time of the target bottom water gas reservoir according to the sum of the water breakthrough times in each stage.

[0069] In this embodiment, after obtaining the water breakthrough time of each stage, the water breakthrough time of each stage can be added up, and the sum is used as the water breakthrough time of the target bottom water gas reservoir. For example, assuming that the water breakthrough time of the three stages is t AB ,tBC and t CD , then the water breakthrough time of the target bottom water gas reservoir can be expressed as t = t AB +t BC +t CD .

[0070] The technical solution of the embodiment of the present invention determines the daily gas production at the bottom of the perforation section of the target bottom water gas reservoir; wherein the target bottom water gas reservoir includes an interlayer; the water invasion process of the target bottom water gas reservoir is divided into three stages, and the water breakthrough time of each stage is determined based on the daily gas production at the bottom of the perforation section; wherein the three stages include the first stage, the second stage and the third stage, the first stage is used to characterize the rise of the gas-water interface from the bottom of the interlayer to the interlayer, the second stage is used to characterize the movement of the gas-water interface from the edge of the interlayer to the center of the interlayer, and the third stage is used to characterize the rise of the gas-water interface from the center of the interlayer to the bottom of the perforation section; the water breakthrough time of the target bottom water gas reservoir is determined according to the sum of the water breakthrough time of each stage. This technical solution, by dividing the bottom water gas reservoir into three stages and calculating the water breakthrough time in stages in combination with the daily gas production changes in the actual production process, can effectively improve the accuracy and applicability of the prediction of the water breakthrough time of the bottom water gas reservoir, so as to provide theoretical support for the development of the bottom water gas reservoir, help optimize the development plan, improve the development efficiency, and promote technological progress and industrial development in the field of gas reservoir development.

[0071] In this embodiment, optionally, whether the gas-water interface has reached the boundary of the first stage is judged based on the rising height of the gas-water interface corresponding to adjacent time periods and the height from the interlayer to the initial gas-water interface, including: if the rising height of the gas-water interface corresponding to the previous time period or the subsequent time period in the adjacent time periods is equal to the height from the interlayer to the initial gas-water interface, then it is determined that the gas-water interface has reached the boundary of the first stage; accordingly, the water breakthrough time of the first stage is determined based on the daily gas production at the bottom of the perforation section, including: the time corresponding to the time when the gas-water interface reaches the boundary of the first stage is determined as the water breakthrough time of the first stage.

[0072] For example, refer to Figure 2 , assuming that the adjacent time periods are i and i+1, i is the earlier time period and i+1 is the later time period. i =h b (At this time z i+1 >h b ) or z i+1 =h b (At this time z i <h b ), it indicates that the first stage boundary is just reached in the previous time period i or the subsequent time period i+1. In this case, the time corresponding to when the gas-water interface reaches the first stage boundary (i.e. the end time of i or i+1) can be directly determined as the water breakthrough time of the first stage.

[0073] In this embodiment, optionally, whether the gas-water interface has reached the second stage boundary is judged based on the distance from the gas-water interface to the edge of the interlayer and the distance from the edge of the interlayer to the center corresponding to adjacent time periods, including: if the distance from the gas-water interface to the edge of the interlayer corresponding to the previous time period or the subsequent time period in the adjacent time periods is equal to the distance from the edge of the interlayer to the center, then it is determined that the gas-water interface has reached the second stage boundary; accordingly, the water breakthrough time of the second stage is determined based on the daily gas production at the bottom of the perforation section, including: determining the time corresponding to the time when the gas-water interface reaches the second stage boundary as the water breakthrough time of the second stage.

[0074] For example, refer to Figure 2 , assuming that the adjacent time periods are i and i+1, i is the earlier time period and i+1 is the later time period. i =r b (At this time x i+1 >r b ) or x i+1 =r b (At this time x i <r b ), it indicates that the first time period i or the second time period i+1 has just reached the boundary of the second stage. In this case, the time corresponding to when the gas-water interface reaches the boundary of the second stage (i.e. the end time of i or i+1) can be directly determined as the water breakthrough time of the second stage.

[0075] In this embodiment, optionally, whether the gas-water interface has reached the boundary of the third stage is judged based on the distance from the gas-water interface to the bottom of the perforation section corresponding to adjacent time periods, including: if the distance from the gas-water interface to the bottom of the perforation section corresponding to the preceding time period or the succeeding time period in the adjacent time periods is equal to zero, then it is determined that the gas-water interface has reached the boundary of the third stage; correspondingly, the water breakthrough time of the third stage is determined based on the daily gas production at the bottom of the perforation section, including: determining the time corresponding to when the gas-water interface reaches the boundary of the third stage as the water breakthrough time of the third stage.

[0076] For example, refer to Figure 2 , assuming that the adjacent time periods are i-1 and i, i-1 is the earlier time period and i is the later time period. i-1 =0 (at this time z i <0) or r i =0 (at this time z i-1 >0), it indicates that the first time period i-1 or the second time period i has just reached the boundary of the third stage. In this case, the time corresponding to when the gas-water interface reaches the boundary of the third stage (i.e., the end time of i-1 or i) can be directly determined as the water breakthrough time of the third stage.

[0077] Embodiment 2

[0078] Figure 4This is a flow chart of a method for determining water breakthrough time of a bottom water gas reservoir provided in the second embodiment of the present invention. This embodiment is optimized based on the above embodiment.

[0079] like Figure 4 As shown, the method of this embodiment specifically includes the following steps:

[0080] S210, determining the daily gas production at the bottom of the perforation section of the target bottom water gas reservoir; wherein the target bottom water gas reservoir includes interlayers.

[0081] S220, dividing the water invasion process of the target bottom water gas reservoir into three stages, and dividing the water invasion process of each stage into a plurality of time periods based on a preset time interval.

[0082] Among them, the three stages include the first stage, the second stage and the third stage. The first stage is used to characterize the rise of the gas-water interface from the bottom of the interlayer to the interlayer, the second stage is used to characterize the movement of the gas-water interface from the edge of the interlayer to the center of the interlayer, and the third stage is used to characterize the rise of the gas-water interface from the center of the interlayer to the bottom of the perforation section.

[0083] S230, determining the rising height of the air-water interface corresponding to adjacent time periods in the first stage based on the first iterative formula.

[0084] The first iteration formula is as follows:

[0085]

[0086] Where l = h c +h b -z, z is the height of the air-water interface, h c is the height from the interlayer to the bottom of the perforated section, h b is the height from the interlayer to the initial air-water interface, i and i+1 represent the previous time period and the next time period in the adjacent time period, respectively, r b is the distance from the edge of the interlayer to the center, q2 is the daily gas production at the bottom of the perforated section, φ is the porosity, S wi is the bound water saturation, S gr is the residual gas saturation, and t is the time.

[0087] S240: If it is determined that the gas-water interface reaches the boundary of the first stage according to the gas-water interface rising height corresponding to the adjacent time periods and the height from the interlayer to the initial gas-water interface, the water breakthrough time of the first stage is determined based on the daily gas production at the bottom of the perforation section.

[0088] S250, determining the distance from the air-water interface to the edge of the interlayer corresponding to the adjacent time periods in the second stage based on the second iterative formula.

[0089] The second iteration formula is as follows:

[0090]

[0091] Among them, h c is the height from the interlayer to the bottom of the perforation section, i and i+1 represent the previous time period and the next time period in the adjacent time period, q2 is the daily gas production at the bottom of the perforation section, φ is the porosity, S wi is the bound water saturation, S gr is the residual gas saturation, t is the time, and x is the distance from the gas-water interface to the edge of the interlayer in the second stage.

[0092] S260: If the gas-water interface reaches the second stage boundary according to the distance from the gas-water interface to the edge of the interlayer and the distance from the edge of the interlayer to the center corresponding to the adjacent time periods, the water breakthrough time of the second stage is determined based on the daily gas production at the bottom of the perforation section.

[0093] S270: Determine the distance from the gas-water interface to the bottom of the perforation section corresponding to the adjacent time periods in the third stage based on the third iterative formula.

[0094] The third iteration formula is as follows:

[0095]

[0096] in, β is the velocity coefficient, ρ g is the gas density, q2 is the daily gas production at the bottom of the perforation section, φ is the porosity, r is the distance from the gas-water interface to the bottom of the perforation section in the third stage, i-1 and i represent the previous time period and the next time period in the adjacent time period, respectively, S wi is the bound water saturation, S gr is the residual gas saturation, and t is the time.

[0097] S280: If it is determined that the gas-water interface reaches the boundary of the third stage according to the distance from the gas-water interface to the bottom of the perforation section corresponding to the adjacent time periods, the water breakthrough time of the third stage is determined based on the daily gas production at the bottom of the perforation section.

[0098] S290: Determine the water breakthrough time of the target bottom water gas reservoir according to the sum of the water breakthrough times in each stage.

[0099] The technical solution of the embodiment of the present invention can effectively improve the accuracy and applicability of the prediction of the water breakthrough time of the bottom water gas reservoir by dividing the bottom water gas reservoir into three stages, and performing quasi-steady-state processing on the three stages respectively, and calculating the water breakthrough time in stages in combination with the daily gas production changes in the actual production process, so as to provide theoretical support for the development of the bottom water gas reservoir, help optimize the development plan, improve the development efficiency, and promote technological progress and industrial development in the field of gas reservoir development.

[0100] Embodiment 3

[0101] Figure 5 This is a schematic diagram of the structure of a device for determining the breakthrough time of a bottom water gas reservoir provided in the third embodiment of the present invention. The device can execute the method for determining the breakthrough time of a bottom water gas reservoir provided in any embodiment of the present invention, and has the corresponding functional modules and beneficial effects of the execution method. Figure 5 As shown, the device comprises:

[0102] The perforation section bottom production determination module 310 is used to determine the daily gas production at the bottom of the perforation section of the target bottom water gas reservoir; wherein the target bottom water gas reservoir includes interlayers;

[0103] The three-stage water breakthrough time determination module 320 is used to divide the water invasion process of the target bottom water gas reservoir into three stages, and determine the water breakthrough time of each stage based on the daily gas production at the bottom of the perforation section; wherein the three stages include the first stage, the second stage and the third stage, the first stage is used to characterize the rise of the gas-water interface from below the interlayer to the interlayer, the second stage is used to characterize the movement of the gas-water interface from the edge of the interlayer to the center of the interlayer, and the third stage is used to characterize the rise of the gas-water interface from the center of the interlayer to the bottom of the perforation section;

[0104] The bottom water gas reservoir water breakthrough time determination module 330 is used to determine the water breakthrough time of the target bottom water gas reservoir according to the sum of the water breakthrough time of each stage.

[0105] Optionally, the three-stage water breakthrough time determination module 320 includes:

[0106] A first stage division unit, used to divide the water invasion process of the first stage into a plurality of time periods based on a preset time interval; wherein the daily gas production at the bottom of the perforation section in each time period is regarded as a fixed value;

[0107] A gas-water interface rise height determination unit, configured to determine the gas-water interface rise heights corresponding to adjacent time periods in the first stage based on the daily gas production at the bottom of the perforation section;

[0108] A first stage boundary judgment unit, used to judge whether the air-water interface reaches the first stage boundary according to the air-water interface rising height corresponding to the adjacent time periods and the height from the interlayer to the initial air-water interface;

[0109] The first stage water breakthrough time determining unit is used to determine the first stage water breakthrough time based on the daily gas production at the bottom of the perforation section if yes.

[0110] Optionally, the air-water interface rise height determination unit is used to:

[0111] The rising height of the air-water interface corresponding to adjacent time periods in the first stage is determined based on a first iterative formula, and the first iterative formula is expressed as follows:

[0112]

[0113] Where l = h c +h b -z, z is the height of the air-water interface, h c is the height from the interlayer to the bottom of the perforated section, h b is the height from the interlayer to the initial air-water interface, i and i+1 represent the previous time period and the next time period in the adjacent time period, respectively, r b is the distance from the edge of the interlayer to the center, q2 is the daily gas production at the bottom of the perforated section, φ is the porosity, S wi is the bound water saturation, S gr is the residual gas saturation, and t is the time.

[0114] Optionally, the first-stage boundary judgment unit is used to:

[0115] If the rising height of the air-water interface corresponding to the previous time period in the adjacent time periods is less than the height from the interlayer to the initial air-water interface, and the rising height of the air-water interface corresponding to the subsequent time period in the adjacent time periods is greater than the height from the interlayer to the initial air-water interface, it is determined that the air-water interface has reached the first stage boundary;

[0116] Accordingly, the first-stage water breakthrough time determination unit is used to:

[0117] The water breakthrough time of the first stage is determined based on the first time formula, and the first time formula is expressed as follows:

[0118]

[0119] Among them, t AB Indicates the water breakthrough time of the first stage.

[0120] Optionally, the first-stage boundary judgment unit is further used to:

[0121] Judging whether the air-water interface reaches the first stage boundary according to the air-water interface rising height corresponding to the adjacent time periods and the height from the interlayer to the initial air-water interface comprises:

[0122] If the rising height of the air-water interface corresponding to the previous time period or the next time period in the adjacent time periods is equal to the height from the interlayer to the initial air-water interface, it is determined that the air-water interface has reached the first stage boundary;

[0123] Accordingly, the first-stage water breakthrough time determination unit is further used for:

[0124] The time corresponding to when the air-water interface reaches the boundary of the first stage is determined as the water breakthrough time of the first stage.

[0125] Optionally, the three-stage water breakthrough time determination module 320 further includes:

[0126] A second stage division unit, for dividing the water invasion process of the second stage into a plurality of time periods based on a preset time interval, wherein the daily gas production at the bottom of the perforation section in each time period is regarded as a fixed value;

[0127] A unit for determining the distance from the gas-water interface to the edge of the interlayer, used to determine the distance from the gas-water interface to the edge of the interlayer corresponding to adjacent time periods in the second stage based on the daily gas production at the bottom of the perforation section;

[0128] A second-stage boundary judgment unit, for judging whether the air-water interface reaches the second-stage boundary according to the distance from the air-water interface to the edge of the interlayer and the distance from the edge of the interlayer to the center corresponding to the adjacent time periods;

[0129] The second stage water breakthrough time determining unit is used to determine the second stage water breakthrough time based on the daily gas production at the bottom of the perforation section if yes.

[0130] Optionally, the unit for determining the distance to the edge of the interlayer is used to:

[0131] The distances to the edge of the interlayer corresponding to adjacent time periods in the second stage are determined based on a second iterative formula, wherein the second iterative formula is expressed as follows:

[0132]

[0133] Among them, h c is the height from the interlayer to the bottom of the perforation section, i and i+1 represent the previous time period and the next time period in the adjacent time period, q2 is the daily gas production at the bottom of the perforation section, φ is the porosity, S wi is the bound water saturation, S gr is the residual gas saturation, t is the time, and x is the distance from the gas-water interface to the edge of the interlayer in the second stage.

[0134] Optionally, the second-stage boundary judgment unit is used to:

[0135] If the distance from the air-water interface to the edge of the interlayer corresponding to the previous time period in the adjacent time periods is less than the distance from the edge of the interlayer to the center, and the distance from the air-water interface to the edge of the interlayer corresponding to the subsequent time period in the adjacent time periods is greater than the distance from the edge of the interlayer to the center, it is determined that the air-water interface reaches the second stage boundary;

[0136] Accordingly, the second-stage water breakthrough time determination unit is used to:

[0137] The water breakthrough time of the second stage is determined based on the second time formula, and the second time formula is expressed as follows:

[0138]

[0139] Among them, t BC represents the water breakthrough time of the second stage, r b Indicates the distance from the edge to the center of the interlayer.

[0140] Optionally, the second-stage boundary judgment unit is further used to:

[0141] If the distance from the air-water interface to the edge of the interlayer corresponding to the previous time period or the next time period in the adjacent time periods is equal to the distance from the edge of the interlayer to the center, it is determined that the air-water interface reaches the boundary of the second stage;

[0142] Accordingly, the second-stage water breakthrough time determination unit is further used for:

[0143] The time corresponding to when the air-water interface reaches the boundary of the second stage is determined as the water breakthrough time of the second stage.

[0144] Optionally, the three-stage water breakthrough time determination module 320 further includes:

[0145] A third stage division unit is used to divide the water invasion process of the third stage into a plurality of time periods based on a preset time interval; wherein the daily gas production at the bottom of the perforation section in each time period is regarded as a fixed value;

[0146] A unit for determining the distance from the gas-water interface to the bottom of the perforation section, used to determine the distance from the gas-water interface to the bottom of the perforation section corresponding to the adjacent time periods in the third stage based on the daily gas production at the bottom of the perforation section;

[0147] A third stage boundary judgment unit, used to judge whether the gas-water interface reaches the third stage boundary according to the distance from the gas-water interface to the bottom of the perforation section corresponding to the adjacent time periods;

[0148] The third stage water breakthrough time determination unit is used to determine the water breakthrough time of the third stage based on the daily gas production at the bottom of the perforation section if yes.

[0149] Optionally, the unit for determining the distance to the bottom of the perforation section is used to:

[0150] The distance from the gas-water interface to the bottom of the perforation section corresponding to the adjacent time periods in the third stage is determined based on the third iterative formula, and the third iterative formula is expressed as follows:

[0151]

[0152] in, β is the velocity coefficient, ρ g is the gas density, q2 is the daily gas production at the bottom of the perforation section, φ is the porosity, r is the distance from the gas-water interface to the bottom of the perforation section in the third stage, i-1 and i represent the previous time period and the next time period in the adjacent time period, respectively, S wi is the bound water saturation, S gr is the residual gas saturation, and t is the time.

[0153] Optionally, the third-stage boundary judgment unit is used to:

[0154] If the distance from the gas-water interface corresponding to the previous time period in the adjacent time periods to the bottom of the perforation section is greater than zero, and the distance from the gas-water interface corresponding to the next time period in the adjacent time periods to the bottom of the perforation section is less than zero, it is determined that the gas-water interface reaches the boundary of the third stage;

[0155] Accordingly, the third-stage water breakthrough time determination unit is used to:

[0156] The water breakthrough time of the third stage is determined based on the third time formula, and the third time formula is expressed as follows:

[0157]

[0158] Among them, t CD Indicates the water breakthrough time of the third stage.

[0159] Optionally, the third-stage boundary judgment unit is further used to:

[0160] If the distance from the gas-water interface corresponding to the previous time period or the next time period in the adjacent time periods to the bottom of the perforation section is zero, it is determined that the gas-water interface has reached the boundary of the third stage;

[0161] Accordingly, the third-stage water breakthrough time determination unit is further used for:

[0162] The time corresponding to when the air-water interface reaches the boundary of the third stage is determined as the water breakthrough time of the third stage.

[0163] A device for determining water breakthrough time of a bottom water gas reservoir provided in an embodiment of the present invention can execute a method for determining water breakthrough time of a bottom water gas reservoir provided in any embodiment of the present invention, and has functional modules and beneficial effects corresponding to the execution method.

[0164] Embodiment 4

[0165] Figure 6A schematic diagram of the structure of an electronic device 10 that can be used to implement an embodiment of the present invention is shown. The electronic device is intended to represent various forms of digital computers, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device can also represent various forms of mobile devices, such as personal digital processing, cellular phones, smart phones, wearable devices (such as helmets, glasses, watches, etc.) and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely examples and are not intended to limit the implementation of the present invention described and / or required herein.

[0166] like Figure 6 As shown, the electronic device 10 includes at least one processor 11, and a memory connected to the at least one processor 11, such as a read-only memory (ROM) 12, a random access memory (RAM) 13, etc., wherein the memory stores a computer program that can be executed by at least one processor, and the processor 11 can perform various appropriate actions and processes according to the computer program stored in the read-only memory (ROM) 12 or the computer program loaded from the storage unit 18 to the random access memory (RAM) 13. In the RAM 13, various programs and data required for the operation of the electronic device 10 can also be stored. The processor 11, the ROM 12, and the RAM 13 are connected to each other through a bus 14. An input / output (I / O) interface 15 is also connected to the bus 14.

[0167] A number of components in the electronic device 10 are connected to the I / O interface 15, including: an input unit 16, such as a keyboard, a mouse, etc.; an output unit 17, such as various types of displays, speakers, etc.; a storage unit 18, such as a disk, an optical disk, etc.; and a communication unit 19, such as a network card, a modem, a wireless communication transceiver, etc. The communication unit 19 allows the electronic device 10 to exchange information / data with other devices through a computer network such as the Internet and / or various telecommunication networks.

[0168] The processor 11 may be a variety of general and / or special processing components with processing and computing capabilities. Some examples of the processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various special artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any appropriate processor, controller, microcontroller, etc. The processor 11 executes the various methods and processes described above, such as a method for determining the breakthrough time of a bottom water gas reservoir.

[0169] In some embodiments, the method for determining the breakthrough time of a bottom water gas reservoir can be implemented as a computer program, which is tangibly contained in a computer-readable storage medium, such as a storage unit 18. In some embodiments, part or all of the computer program can be loaded and / or installed on the electronic device 10 via the ROM 12 and / or the communication unit 19. When the computer program is loaded into the RAM 13 and executed by the processor 11, one or more steps of the method for determining the breakthrough time of a bottom water gas reservoir described above can be performed. Alternatively, in other embodiments, the processor 11 can be configured to execute the method for determining the breakthrough time of a bottom water gas reservoir by any other appropriate means (for example, by means of firmware).

[0170] Various implementations of the systems and techniques described above herein can be implemented in digital electronic circuit systems, integrated circuit systems, field programmable gate arrays (FPGAs), application specific integrated circuits (ASICs), application specific standard products (ASSPs), systems on chips (SOCs), load programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various implementations can include: being implemented in one or more computer programs that can be executed and / or interpreted on a programmable system including at least one programmable processor, which can be a special purpose or general purpose programmable processor that can receive data and instructions from a storage system, at least one input device, and at least one output device, and transmit data and instructions to the storage system, the at least one input device, and the at least one output device.

[0171] Computer programs for implementing the methods of the present invention may be written in any combination of one or more programming languages. These computer programs may be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, so that when the computer program is executed by the processor, the functions / operations specified in the flow chart and / or block diagram are implemented. The computer program may be executed entirely on the machine, partially on the machine, partially on the machine and partially on a remote machine as a stand-alone software package, or entirely on a remote machine or server.

[0172] In the context of the present invention, a computer-readable storage medium may be a tangible medium that may contain or store a computer program for use by or in combination with an instruction execution system, device or equipment. A computer-readable storage medium may include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, devices or equipment, or any suitable combination of the foregoing. Alternatively, a computer-readable storage medium may be a machine-readable signal medium. A more specific example of a machine-readable storage medium may include an electrical connection based on one or more lines, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

[0173] To provide interaction with a user, the systems and techniques described herein may be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and a pointing device (e.g., a mouse or trackball) through which the user can provide input to the electronic device. Other types of devices may also be used to provide interaction with the user; for example, the feedback provided to the user may be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user may be received in any form (including acoustic input, voice input, or tactile input).

[0174] The systems and techniques described herein may be implemented in a computing system that includes backend components (e.g., as a data server), or a computing system that includes middleware components (e.g., an application server), or a computing system that includes frontend components (e.g., a user computer with a graphical user interface or a web browser through which a user can interact with implementations of the systems and techniques described herein), or a computing system that includes any combination of such backend components, middleware components, or frontend components. The components of the system may be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include: a local area network (LAN), a wide area network (WAN), a blockchain network, and the Internet.

[0175] A computing system may include a client and a server. The client and the server are generally remote from each other and usually interact through a communication network. The client and server relationship is generated by computer programs running on the corresponding computers and having a client-server relationship with each other. The server may be a cloud server, also known as a cloud computing server or cloud host, which is a host product in the cloud computing service system to solve the defects of difficult management and weak business scalability in traditional physical hosts and VPS services.

[0176] It should be understood that the various forms of processes shown above can be used to reorder, add or delete steps. For example, the steps described in the present invention can be executed in parallel, sequentially or in different orders, as long as the desired results of the technical solution of the present invention can be achieved, and this document does not limit this.

[0177] The above specific implementations do not constitute a limitation on the protection scope of the present invention. It should be understood by those skilled in the art that various modifications, combinations, sub-combinations and substitutions can be made according to design requirements and other factors. Any modification, equivalent substitution and improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims

1. A method for determining water breakthrough time of a bottom water gas reservoir, characterized in that: The method comprises: Determine the daily gas production at the bottom of the perforation section of the target bottom water gas reservoir; wherein the target bottom water gas reservoir includes interlayers; The water invasion process of the target bottom water gas reservoir is divided into three stages, and the water breakthrough time of each stage is determined based on the daily gas production at the bottom of the perforation section; wherein the three stages include the first stage, the second stage and the third stage, the first stage is used to characterize the rise of the gas-water interface from below the interlayer to the interlayer, the second stage is used to characterize the movement of the gas-water interface from the edge of the interlayer to the center of the interlayer, and the third stage is used to characterize the rise of the gas-water interface from the center of the interlayer to the bottom of the perforation section; The water breakthrough time of the target bottom water gas reservoir is determined according to the sum of the water breakthrough times of each stage.

2. The method according to claim 1, characterized in that The water breakthrough time of each stage is determined based on the daily gas production at the bottom of the perforation section, including: The water invasion process of the first stage is divided into a plurality of time periods based on a preset time interval; wherein the daily gas production at the bottom of the perforation section in each time period is regarded as a fixed value; Determining the rising height of the gas-water interface corresponding to adjacent time periods in the first stage based on the daily gas production at the bottom of the perforation section; According to the rising height of the air-water interface corresponding to the adjacent time periods and the height from the interlayer to the initial air-water interface, it is judged whether the air-water interface reaches the first stage boundary; If yes, the water breakthrough time of the first stage is determined based on the daily gas production at the bottom of the perforation section.

3. The method according to claim 2, characterized in that Determining the rising height of the gas-water interface corresponding to adjacent time periods in the first stage based on the daily gas production at the bottom of the perforation section includes: The rising height of the air-water interface corresponding to adjacent time periods in the first stage is determined based on a first iterative formula, and the first iterative formula is expressed as follows: Where l = h c +h b -z, z is the height of the air-water interface, h c is the height from the interlayer to the bottom of the perforated section, h b is the height from the interlayer to the initial air-water interface, i and i+1 represent the previous time period and the next time period in the adjacent time period, respectively, r b is the distance from the edge of the interlayer to the center, q2 is the daily gas production at the bottom of the perforated section, φ is the porosity, S wi is the bound water saturation, S gr is the residual gas saturation, and t is the time.

4. The method according to claim 3, characterized in that: Judging whether the air-water interface reaches the first stage boundary according to the air-water interface rising height corresponding to the adjacent time periods and the height from the interlayer to the initial air-water interface comprises: If the rising height of the air-water interface corresponding to the previous time period in the adjacent time periods is less than the height from the interlayer to the initial air-water interface, and the rising height of the air-water interface corresponding to the subsequent time period in the adjacent time periods is greater than the height from the interlayer to the initial air-water interface, it is determined that the air-water interface has reached the first stage boundary; Accordingly, determining the water breakthrough time of the first stage based on the daily gas production at the bottom of the perforation section includes: The water breakthrough time of the first stage is determined based on the first time formula, and the first time formula is expressed as follows: Among them, t AB Indicates the water breakthrough time of the first stage.

5. The method according to claim 3, characterized in that: Judging whether the air-water interface reaches the first stage boundary according to the air-water interface rising height corresponding to the adjacent time periods and the height from the interlayer to the initial air-water interface comprises: If the rising height of the air-water interface corresponding to the previous time period or the next time period in the adjacent time periods is equal to the height from the interlayer to the initial air-water interface, it is determined that the air-water interface has reached the first stage boundary; Accordingly, determining the water breakthrough time of the first stage based on the daily gas production at the bottom of the perforation section includes: The time corresponding to when the air-water interface reaches the boundary of the first stage is determined as the water breakthrough time of the first stage.

6. The method according to claim 1, characterized in that The water breakthrough time of each stage is determined based on the daily gas production at the bottom of the perforation section, including: The water invasion process of the second stage is divided into a plurality of time periods based on a preset time interval, wherein the daily gas production at the bottom of the perforation section in each time period is regarded as a fixed value; Determine the distance from the gas-water interface to the edge of the interlayer corresponding to the adjacent time periods in the second stage based on the daily gas production at the bottom of the perforation section; According to the distance from the air-water interface to the edge of the interlayer and the distance from the edge of the interlayer to the center corresponding to the adjacent time periods, it is determined whether the air-water interface reaches the boundary of the second stage; If yes, the water breakthrough time of the second stage is determined based on the daily gas production at the bottom of the perforation section.

7. The method according to claim 6, characterized in that Determining the distance from the gas-water interface to the edge of the interlayer corresponding to the adjacent time periods in the second stage based on the daily gas production at the bottom of the perforation section includes: The distance from the air-water interface to the edge of the interlayer corresponding to the adjacent time periods in the second stage is determined based on the second iterative formula, and the second iterative formula is expressed as follows: Among them, h c is the height from the interlayer to the bottom of the perforation section, i and i+1 represent the previous time period and the next time period in the adjacent time period, q2 is the daily gas production at the bottom of the perforation section, φ is the porosity, S wi is the bound water saturation, S gr is the residual gas saturation, t is the time, and x is the distance from the gas-water interface to the edge of the interlayer in the second stage.

8. The method according to claim 7, characterized in that Judging whether the air-water interface reaches the second stage boundary according to the distance from the air-water interface to the edge of the interlayer and the distance from the edge of the interlayer to the center corresponding to the adjacent time periods includes: If the distance from the air-water interface to the edge of the interlayer corresponding to the previous time period in the adjacent time periods is less than the distance from the edge of the interlayer to the center, and the distance from the air-water interface to the edge of the interlayer corresponding to the subsequent time period in the adjacent time periods is greater than the distance from the edge of the interlayer to the center, it is determined that the air-water interface reaches the second stage boundary; Accordingly, determining the water breakthrough time of the second stage based on the daily gas production at the bottom of the perforation section includes: The water breakthrough time of the second stage is determined based on the second time formula, and the second time formula is expressed as follows: Among them, t BC represents the water breakthrough time of the second stage, r b Indicates the distance from the edge to the center of the interlayer.

9. The method according to claim 7, characterized in that: Judging whether the air-water interface reaches the second stage boundary according to the distance from the air-water interface to the edge of the interlayer and the distance from the edge of the interlayer to the center corresponding to the adjacent time periods includes: If the distance from the air-water interface to the edge of the interlayer corresponding to the previous time period or the next time period in the adjacent time periods is equal to the distance from the edge of the interlayer to the center, it is determined that the air-water interface reaches the boundary of the second stage; Accordingly, determining the water breakthrough time of the second stage based on the daily gas production at the bottom of the perforation section includes: The time corresponding to when the air-water interface reaches the boundary of the second stage is determined as the water breakthrough time of the second stage.

10. The method according to claim 1, characterized in that The water breakthrough time of each stage is determined based on the daily gas production at the bottom of the perforation section, including: The water invasion process of the third stage is divided into a plurality of time periods based on a preset time interval; wherein the daily gas production at the bottom of the perforation section in each time period is regarded as a fixed value; Determine the distance from the gas-water interface to the bottom of the perforation section corresponding to the adjacent time periods in the third stage based on the daily gas production at the bottom of the perforation section; Judging whether the gas-water interface has reached the third stage boundary according to the distance from the gas-water interface corresponding to the adjacent time periods to the bottom of the perforation section; If yes, the water breakthrough time of the third stage is determined based on the daily gas production at the bottom of the perforation section.

11. The method according to claim 10, characterized in that Determining the distance from the gas-water interface to the bottom of the perforation section corresponding to the adjacent time periods in the third stage based on the daily gas production at the bottom of the perforation section includes: The distance from the gas-water interface to the bottom of the perforation section corresponding to the adjacent time periods in the third stage is determined based on the third iterative formula, and the third iterative formula is expressed as follows: in, β is the velocity coefficient, ρ g is the gas density, q2 is the daily gas production at the bottom of the perforation section, φ is the porosity, r is the distance from the gas-water interface to the bottom of the perforation section in the third stage, i-1 and i represent the previous time period and the next time period in the adjacent time period, respectively, S wi is the bound water saturation, S gr is the residual gas saturation, and t is the time.

12. The method according to claim 11, characterized in that Judging whether the gas-water interface reaches the third stage boundary according to the distance from the gas-water interface corresponding to the adjacent time periods to the bottom of the perforation section includes: If the distance from the gas-water interface corresponding to the previous time period in the adjacent time periods to the bottom of the perforation section is greater than zero, and the distance from the gas-water interface corresponding to the next time period in the adjacent time periods to the bottom of the perforation section is less than zero, it is determined that the gas-water interface reaches the boundary of the third stage; Accordingly, determining the water breakthrough time of the third stage based on the daily gas production at the bottom of the perforation section includes: The water breakthrough time of the third stage is determined based on the third time formula, and the third time formula is expressed as follows: Among them, t CD Indicates the water breakthrough time of the third stage.

13. The method according to claim 11, characterized in that Judging whether the gas-water interface reaches the third stage boundary according to the distance from the gas-water interface corresponding to the adjacent time periods to the bottom of the perforation section includes: If the distance from the gas-water interface corresponding to the previous time period or the next time period in the adjacent time periods to the bottom of the perforation section is zero, it is determined that the gas-water interface has reached the boundary of the third stage; Accordingly, determining the water breakthrough time of the third stage based on the daily gas production at the bottom of the perforation section includes: The time corresponding to when the air-water interface reaches the boundary of the third stage is determined as the water breakthrough time of the third stage.

14. A device for determining water breakthrough time of a bottom water gas reservoir, characterized in that: The device comprises: A perforation section bottom production determination module is used to determine the daily gas production at the bottom of the perforation section of a target bottom water gas reservoir; wherein the target bottom water gas reservoir includes an interlayer; A three-stage water breakthrough time determination module is used to divide the water invasion process of the target bottom water gas reservoir into three stages, and determine the water breakthrough time of each stage based on the daily gas production at the bottom of the perforation section; wherein the three stages include a first stage, a second stage and a third stage, the first stage is used to characterize the rise of the gas-water interface from below the interlayer to the interlayer, the second stage is used to characterize the movement of the gas-water interface from the edge of the interlayer to the center of the interlayer, and the third stage is used to characterize the rise of the gas-water interface from the center of the interlayer to the bottom of the perforation section; The bottom water gas reservoir water breakthrough time determination module is used to determine the water breakthrough time of the target bottom water gas reservoir according to the sum of the water breakthrough time of each stage.

15. An electronic device, characterized in that: The electronic device comprises: at least one processor; and, a memory communicatively connected to the at least one processor; wherein, The memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can execute the method for determining the water breakthrough time of a bottom water gas reservoir according to any one of claims 1-13.

16. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a processor to implement the method for determining water breakthrough time of a bottom water gas reservoir according to any one of claims 1-13 when executed.

Images