A method and device for determining water invasion wave propagation coefficient, electronic equipment and storage medium

Abstract

The application discloses a water invasion wave and coefficient determination method and device, electronic equipment and a storage medium. The method comprises the following steps: obtaining reservoir basic parameters, geological layering data and well logging curves of a water drive gas reservoir; calculating the permeability average and thickness of each reservoir of the water drive gas reservoir according to the geological layering data and well logging curves; for each reservoir of the water drive gas reservoir, determining the water production gas ratio when water breakthrough of the reservoir according to the reservoir basic parameters, the permeability average and thickness of the reservoir; determining the water breakthrough layer when abandonment and the water non-breakthrough layer when abandonment from all the reservoirs of the water drive gas reservoir according to the water production gas ratio when water breakthrough; and calculating the water invasion wave and coefficient when abandonment according to the reservoir basic parameters, the permeability average of the water non-breakthrough layer when abandonment and the thickness of each reservoir of the water drive gas reservoir. The scheme can accurately calculate the water invasion wave and coefficient when abandonment of the water drive gas reservoir, and provides technical support for the recovery efficiency calibration and recovery efficiency improvement of the water drive gas reservoir.

Description

Technical Field

[0001] The present invention relates to the technical field of oil and gas reservoir development, and in particular to a method, device, electronic equipment and storage medium for determining a water invasion sweep coefficient. Background Art

[0002] The water invasion sweep coefficient is a key parameter for calibrating the recovery factor of water-flooded gas reservoirs. It is dimensionless and ranges from 0 to 1. To a certain extent, the water invasion sweep coefficient directly determines the development effect and ultimate recovery factor of water-flooded gas reservoirs. Therefore, accurately predicting the water invasion sweep coefficient of water-flooded gas reservoirs is crucial for calibrating the recovery factors of edge- and bottom-water gas reservoirs and formulating recovery enhancement strategies.

[0003] Currently, the water invasion sweep coefficient for water-driven gas reservoirs can only be calculated based on water invasion and water production. Related technologies provide a formula for calculating the water invasion sweep coefficient based on dynamic water invasion and water production, but this primarily evaluates the current sweep coefficient. Without knowing the water invasion and water production at the time of abandonment, it is impossible to predict the sweep coefficient of a water-driven gas reservoir at the time of abandonment, directly impacting the calibration of the recovery factor and the evaluation of the enhanced oil recovery effect of the water-driven gas reservoir. Therefore, there is an urgent need to develop a method for predicting the sweep coefficient of a water-driven gas reservoir at the time of abandonment. Summary of the Invention

[0004] The present invention provides a method, device, electronic equipment and storage medium for determining a water invasion sweep coefficient, which can accurately predict the water invasion sweep coefficient when a water drive gas reservoir is abandoned.

[0005] According to one aspect of the present invention, a method for determining a water intrusion sweep coefficient is provided, comprising:

[0006] Obtain basic reservoir parameters, geological stratification data and well logging curves of water-driven gas reservoirs;

[0007] Calculating the mean permeability and thickness of each reservoir in the water-flooded gas reservoir based on the geological stratification data and the well logging curve;

[0008] For each reservoir of the water-drive gas reservoir, determining the production water-gas ratio of the reservoir when water invades water according to the basic parameters of the reservoir, the average permeability and the thickness of the reservoir;

[0009] Determining, from all reservoirs of the water-flooded gas reservoir, layers where water has been invaded during abandonment and layers where water has not been invaded during abandonment, based on the water-gas ratio produced during water-flooding.

[0010] The water invasion sweep coefficient at abandonment is calculated according to the basic reservoir parameters, the average permeability of the water-invasion-unseen water layer at abandonment, and the thickness of each reservoir of the water-flooding gas reservoir.

[0011] According to another aspect of the present invention, there is provided a device for determining a water invasion sweep coefficient, comprising:

[0012] Reservoir parameter acquisition module, used to obtain basic reservoir parameters, geological stratification data and logging curves of water-driven gas reservoirs;

[0013] a permeability mean calculation module, configured to calculate the permeability mean and thickness of each reservoir of the water-flooded gas reservoir based on the geological stratification data and the well logging curve;

[0014] a production water-gas ratio determination module, configured to determine, for each reservoir of the water-drive gas reservoir, the production water-gas ratio of the reservoir when water invades water according to the basic parameters of the reservoir, the average permeability value and the thickness of the reservoir;

[0015] a reservoir classification module for determining, from all reservoirs of the water-flooding gas reservoir, layers where water was seen during abandonment and layers where water was not seen during abandonment, based on the water-gas ratio of production during water-flooding;

[0016] The water invasion sweep coefficient calculation module is used to calculate the water invasion sweep coefficient at abandonment based on the reservoir basic parameters, the average permeability of the water invasion unseen water layer at abandonment, and the thickness of each reservoir of the water drive gas reservoir.

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

[0018] at least one processor; and

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

[0020] The memory stores a computer program executable by the at least one processor. The computer program is executed by the at least one processor so that the at least one processor can perform the method for determining the water invasion sweep coefficient according to any embodiment of the present invention.

[0021] 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 enable a processor to implement the method for determining the water intrusion sweep coefficient according to any embodiment of the present invention when executed.

[0022] The scheme for determining the water invasion sweep coefficient of the embodiment of the present invention includes: obtaining basic reservoir parameters, geological stratification data and well logging curves of a water-driven gas reservoir; calculating the mean permeability and thickness of each reservoir of the water-driven gas reservoir based on the geological stratification data and the well logging curves; for each reservoir of the water-driven gas reservoir, determining the water-to-gas ratio of production when water invades water in the reservoir based on the basic reservoir parameters, the mean permeability and thickness of the reservoir; determining the water-invaded water layer and the water-uninvaded water layer from all reservoirs of the water-driven gas reservoir based on the water-invaded water ratio; calculating the water invasion sweep coefficient at abandonment based on the basic reservoir parameters, the mean permeability of the water-uninvaded water layer and the thickness of each reservoir of the water-driven gas reservoir. The technical scheme provided by the embodiment of the present invention can accurately calculate the water invasion sweep coefficient of a water-driven gas reservoir at abandonment, provide technical support for the calibration and improvement of the recovery rate of water-driven gas reservoirs, and effectively improve the development effect of water-driven gas reservoirs.

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

[0024] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. 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.

[0025] Figure 1 This is a flow chart of a method for determining a water intrusion sweep coefficient according to the first embodiment of the present invention;

[0026] Figure 2 is a schematic diagram of the well logging curve before reordering;

[0027] Figure 3 1 is a schematic diagram of the re-arranged reservoir stratification and the water-gas ratio of production when water invades water, provided by an embodiment of the present invention;

[0028] Figure 4 Schematic diagram of the corresponding relationship between the natural gas compression factor Z and the formation pressure P provided by an embodiment of the present invention;

[0029] Figure 5 is the natural gas viscosity μ provided by the embodiment of the present invention g Schematic diagram of the corresponding relationship between it and the formation pressure P;

[0030] Figure 6Schematic diagram of water invasion during abandonment of a water-driven gas reservoir according to an embodiment of the present invention;

[0031] Figure 7 2 is a schematic structural diagram of a device for determining a water intrusion sweep coefficient according to a second embodiment of the present invention;

[0032] Figure 8 It is a structural diagram of an electronic device for implementing the method for determining the water intrusion sweep coefficient according to an embodiment of the present invention. DETAILED DESCRIPTION

[0033] In order to enable those skilled in the art to better understand the solutions of the present invention, the technical solutions 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 embodiments described 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 making creative efforts should fall within the scope of protection of the present invention.

[0034] It should be noted that the terms "first", "second", etc. in the description 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 numbers 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.

[0035] Example 1

[0036] Figure 1 A flow chart of a method for determining a water intrusion sweep coefficient is provided for the first embodiment of the present invention. This embodiment is applicable to situations where a water intrusion sweep coefficient is to be determined. The method can be executed by a device for determining a water intrusion sweep coefficient. The device can be implemented in the form of hardware and / or software and can be configured in an electronic device. Figure 1 As shown, the method includes:

[0037] S110. Obtain basic reservoir parameters, geological stratification data, and well logging curves of the water-drive gas reservoir.

[0038] In the embodiment of the present invention, the reservoir basic parameters, geological stratification data and well logging curves of the water drive gas reservoir are obtained, wherein the reservoir basic parameters can also be called reservoir basic physical property parameters. The reservoir basic parameters can include the original gas saturation S of the reservoir gi , the gas phase relative permeability of the reservoir at the original gas saturation Krgi, the residual gas saturation of the reservoir after water flooding S gr , relative permeability of reservoir water phase after water flooding K rwc , reservoir temperature T, original formation pressure P i , reservoir porosity Formation water viscosity μ w , formation water volume coefficient B w , natural gas PVT data. Natural gas PVT data includes natural gas compressibility, formation pressure, and natural gas viscosity.

[0039] S120. Calculate the mean permeability and thickness of each reservoir in the water-drive gas reservoir based on the geological stratification data and the well logging curve.

[0040] In the embodiment of the present invention, the permeability average of each reservoir of the porous water-drive gas reservoir is calculated from top to bottom based on the well logging curve and geological stratification data by using the arithmetic mean method, and the thickness of each reservoir is calculated based on the top depth and bottom depth of each reservoir. The permeability average and thickness of each reservoir can be respectively recorded as k j and h j , where j is a serial number, taking values ​​of 1, 2, …, N, where N represents the number of geological layers.

[0041] S130. For each reservoir of the water-drive gas reservoir, determine the production water-gas ratio of the reservoir when water invades water according to the basic parameters of the reservoir, the average permeability of the reservoir, and the thickness of the reservoir.

[0042] In an embodiment of the present invention, the water-gas ratio produced when water invades water in the reservoir is determined based on the basic reservoir parameters, the average permeability and thickness of each reservoir of the water-drive gas reservoir. Optionally, the water-gas ratio produced when water invades water in each reservoir of the water-drive gas reservoir is determined based on the basic reservoir parameters, the average permeability and thickness of the reservoir, including: sorting each reservoir of the water-drive gas reservoir in the order of the average permeability from large to small or from small to large; determining the water-gas ratio produced when water invades water in each reservoir of the water-drive gas reservoir according to the basic reservoir parameters, the average permeability and thickness of the reservoir. Exemplarily, the reservoirs are re-sorted in the order of the average permeability of each reservoir from large to small, and the corresponding average permeability and thickness of the re-sorted reservoirs are respectively recorded as K j and H jDue to the difference in vertical permeability of each reservoir, the layer with the highest permeability in the porous water-drive gas reservoir is the first to be invaded by water, followed by the layer with the second highest permeability, and so on. According to the order of water invasion of each reservoir, based on the basic parameters of the reservoir, the thickness of each reservoir, and the average permeability of each reservoir, the production water-gas ratio of each reservoir when water invades is calculated. For example, Figure 2 This is a schematic diagram of the logging curve before reordering. Figure 3 A schematic diagram of the reordered reservoir layers and the water-gas ratio of production when water invades water provided by an embodiment of the present invention.

[0043] Optionally, the basic reservoir parameters include the gas phase relative permeability of the reservoir at the original gas saturation, the water phase relative permeability of the reservoir after water flooding, the reservoir temperature, the original formation pressure, the formation water viscosity, the formation water volume coefficient, the natural gas volume coefficient and the natural gas PVT data; the determination of the production water-gas ratio of the reservoir when water invades water according to the order of water invasion of each reservoir of the water-flooded gas reservoir, based on the basic reservoir parameters, the average permeability and the thickness of the reservoir, includes: establishing a first corresponding relationship between the natural gas compressibility factor and the formation pressure and a second corresponding relationship between the natural gas viscosity and the formation pressure based on the natural gas PVT data; determining the natural gas compressibility factor according to the first corresponding relationship and the original formation pressure, and determining the natural gas viscosity according to the second corresponding relationship and the original formation pressure; and determining the production water-gas ratio of the reservoir when water invades water according to the following formula:

[0044]

[0045] Among them, H i represents the thickness of the i-th reservoir after sorting, μ w Indicates the formation water viscosity, K rwc represents the relative permeability of water phase in the reservoir after water flooding, B w represents the formation water volume coefficient, K rgi represents the gas phase relative permeability of the reservoir at the original gas saturation, μ g Represents the viscosity of reservoir natural gas, P i represents the original formation pressure, Z represents the reservoir natural gas compression factor, P sc Indicates standard atmospheric pressure, T sc represents the standard temperature of the oil industry, T represents the reservoir temperature, B g represents the natural gas volume coefficient, α j represents the water-gas ratio of production when water invades the jth reservoir, and N represents the total number of reservoirs.

[0046] In the embodiment of the present invention, natural gas PVT data includes natural gas compressibility factor Z, natural gas viscosity μ gand formation pressure P, therefore, the first correspondence between natural gas compressibility factor Z and formation pressure P and the natural gas viscosity μ can be established based on natural gas PVT data. g The second corresponding relationship between and the formation pressure P. For example, Figure 4 A schematic diagram of the corresponding relationship between the natural gas compression factor Z and the formation pressure P provided in an embodiment of the present invention, Figure 5 The natural gas viscosity μ provided in the embodiment of the present invention g Schematic diagram of the corresponding relationship between the natural gas compressibility factor Z and the formation pressure P. Among them, the binomial fitting algorithm can be used to establish the first corresponding relationship between the natural gas compressibility factor Z and the formation pressure P and the natural gas viscosity μ g The second corresponding relationship between the natural gas compression factor Z and the formation pressure P has a fitting correlation coefficient greater than 0.90. For example, the first corresponding relationship between the natural gas compression factor Z and the formation pressure P can be expressed as: Z = 0.8084 + 0.0048P + 0.00054P 2 ; Natural gas viscosity μ g The second corresponding relationship between μ and formation pressure P can be expressed as: g =0.01063+0.000349P; where Z represents the natural gas compression factor, f; μ g represents the viscosity of natural gas, mPa.s; P represents the formation pressure, MPa; a0, a1, a2, b0, b1, and b2 are the corresponding fitting parameters.

[0047] In this embodiment of the present invention, the natural gas compressibility factor is determined based on the first correspondence and the original formation pressure, and the natural gas viscosity is determined based on the second correspondence and the original formation pressure. The water-to-gas ratio produced when water intrudes into the reservoir is then determined based on the natural gas compressibility factor and the natural gas viscosity according to the following formula: Among them, H i represents the thickness of the i-th reservoir after sorting, μ w Indicates the formation water viscosity, K rwc represents the relative permeability of water phase in the reservoir after water flooding, B w represents the formation water volume coefficient, K rgi represents the gas phase relative permeability of the reservoir at the original gas saturation, μ g Represents the viscosity of reservoir natural gas, P i represents the original formation pressure, Z represents the reservoir natural gas compression factor, P sc Indicates standard atmospheric pressure, T sc represents the standard temperature of the oil industry, T represents the reservoir temperature, B g represents the natural gas volume coefficient, α j represents the water-gas ratio of production when water invades the jth reservoir, and N represents the total number of reservoirs.

[0048] S140. Determine, from all reservoirs of the water drive gas reservoir, layers where water has been invaded during abandonment and layers where water has not been invaded during abandonment based on the water-to-gas ratio produced during water invasion.

[0049] For example, the water-gas ratio of each reservoir produced when water invades water is compared with a preset threshold value, and the reservoir whose water-gas ratio is greater than the preset threshold value when water invades water is regarded as a water-invaded layer when abandoned, and the reservoir whose water-gas ratio is less than the preset threshold value when water invades water is regarded as a water-uninvaded layer when abandoned.

[0050] Optionally, the method of determining the abandoned water-invasion layer and the abandoned water-uninvasion layer from all reservoirs of the water-drive gas reservoir based on the water-invasion water-gas ratio, includes: obtaining the abandoned water-gas ratio; comparing the water-invasion water-gas ratio of each reservoir of the water-drive gas reservoir with the abandoned water-gas ratio, and treating the reservoir with the water-invasion water-gas ratio greater than the abandoned water-gas ratio as the abandoned water-invasion layer, and treating the reservoir with the water-invasion water-gas ratio less than the abandoned water-gas ratio as the abandoned water-uninvasion layer. Exemplarily, the abandoned water-gas ratio WGR input by the user is obtained, wherein the abandoned water-gas ratio WGR can be understood as the abandonment condition of the water-drive gas reservoir. For example, the abandoned water-gas ratio WGR can be 5.0m 3 / 10 4 m 3 For each reservoir of the water-drive gas reservoir, the water-gas production ratio when water invades water and the abandoned water-gas ratio of the reservoir are compared. When the water-gas production ratio when water invades water is greater than the abandoned water-gas ratio, the reservoir is regarded as a water-invaded layer when water is abandoned. When the water-gas production ratio when water invades water is less than the abandoned water-gas ratio, the reservoir is regarded as a water-uninvaded layer when water is abandoned.

[0051] S150, calculating the water invasion sweep coefficient at abandonment based on the reservoir basic parameters, the average permeability of the water-invasion-unseen water layer at abandonment, and the thickness of each reservoir of the water-flooding gas reservoir.

[0052] Optionally, the basic reservoir parameters include the original gas saturation of the reservoir and the residual gas saturation of the reservoir after water flooding; before calculating the water invasion sweep coefficient at abandonment based on the basic reservoir parameters, the average permeability of the water-invaded unseen water layer at abandonment, and the thickness of each reservoir of the water-flooded gas reservoir, the following further comprises: obtaining the minimum bottom hole pressure allowed for production and the water body parameters of the water-flooded gas reservoir; wherein the water body parameters of the water-flooded gas reservoir include the water body multiple, the pore compressibility coefficient, and the formation water compressibility coefficient; according to the water body parameters of the water-flooded gas reservoir, the original gas saturation of the reservoir, the residual gas saturation of the reservoir after water flooding ... calculating the ultimate water invasion sweep coefficient of the water-drive gas reservoir based on the gas saturation, the residual gas saturation of the reservoir after water flooding, and the minimum bottom hole pressure allowed for production; and calculating the water invasion sweep coefficient at the time of abandonment based on the basic reservoir parameters, the average permeability of the unseen water layers after water flooding at the time of abandonment, and the thickness of each reservoir layer of the water-drive gas reservoir, comprising: calculating the water invasion sweep coefficient at the time of abandonment based on the ultimate water invasion sweep coefficient of the water-drive gas reservoir, the basic reservoir parameters, the average permeability of the unseen water layers after water flooding at the time of abandonment, and the thickness of each reservoir layer of the water-drive gas reservoir.

[0053] In the embodiment of the present invention, the minimum bottom hole pressure allowed for production and water body parameters of the water drive gas reservoir input by the user are obtained, wherein the water body parameters of the water drive gas reservoir include the water body multiple N, the pore compression coefficient C p and formation water compressibility C w The minimum bottom hole pressure allowed for production is another abandonment condition for water drive gas reservoirs. According to the water parameters of porous water drive gas reservoirs (water body multiples, pore compressibility coefficient, formation water compressibility coefficient), the original water saturation of the reservoir, the residual gas saturation after water drive, the minimum bottom hole pressure allowed for production P a , the ultimate water invasion sweep coefficient of porous water drive gas reservoir is calculated from the perspective of water body energy.

[0054] Optionally, the calculating the ultimate water invasion sweep coefficient of the water-flooding gas reservoir according to the water parameters of the water-flooding gas reservoir, the original gas saturation of the reservoir, the residual gas saturation of the reservoir after water flooding, and the minimum bottom hole pressure allowed for production includes: calculating the ultimate water invasion sweep coefficient of the water-flooding gas reservoir according to the following formula:

[0055]

[0056] Among them, N represents the water volume multiple, C w represents the formation water compressibility coefficient, C p represents the pore compressibility coefficient, P i Indicates the original formation pressure, P a Indicates the minimum bottom hole pressure allowed for production, S gi Indicates the original gas saturation of the reservoir, S gr It represents the residual gas saturation of the reservoir after water flooding, and β1 represents the ultimate water invasion coefficient of the water flooding gas reservoir.

[0057] Optional. Calculating the water invasion sweep coefficient at abandonment based on the water flooding gas reservoir's ultimate water invasion sweep coefficient, the reservoir basic parameters, the average permeability of the water-invasion-unseen water layer at abandonment, and the thickness of each reservoir of the water flooding gas reservoir includes: calculating the water invasion sweep coefficient at abandonment based on the following formula:

[0058]

[0059] Among them, β1 represents the ultimate water invasion sweep coefficient of the water drive gas reservoir, β2 represents the water invasion sweep coefficient when abandoned, the 1st reservoir to the nth reservoir are the water invasion visible layers when abandoned, the n+1th reservoir to the Nth reservoir are the water invasion invisible layers when abandoned, H i represents the thickness of the i-th reservoir, K i represents the mean permeability of the i-th reservoir, K n+1 , represents the mean permeability of the n+1th reservoir, and N represents the total number of reservoirs in the water-flooded gas reservoir.

[0060] For example, Figure 6 A schematic diagram of water invasion during abandonment of a water-driven gas reservoir according to an embodiment of the present invention.

[0061] Example 2

[0062] Figure 7 This is a schematic diagram of the structure of a device for determining a water invasion sweep coefficient provided in the second embodiment of the present invention. Figure 7 As shown, the device includes:

[0063] Reservoir parameter acquisition module 710, used to obtain basic reservoir parameters, geological layer data and well logging curves of water-driven gas reservoirs;

[0064] a permeability mean calculation module 720, configured to calculate the permeability mean and thickness of each reservoir of the water-flooded gas reservoir based on the geological stratification data and the well logging curve;

[0065] The production water-gas ratio determination module 730 is configured to determine, for each reservoir of the water-flooding gas reservoir, the production water-gas ratio of the reservoir when water invades the reservoir based on the reservoir basic parameters, the average permeability and the thickness of the reservoir;

[0066] A reservoir classification module 740 is configured to determine, from all reservoirs of the water flooding gas reservoir, layers with water invasion at the time of abandonment and layers without water invasion at the time of abandonment based on the water-gas ratio of production at the time of water invasion;

[0067] The water invasion sweep coefficient calculation module 750 is used to calculate the water invasion sweep coefficient at abandonment based on the reservoir basic parameters, the average permeability of the water invasion unseen water layer at abandonment, and the thickness of each reservoir of the water drive gas reservoir.

[0068] Optionally, the reservoir classification module is used to:

[0069] Get the waste water-gas ratio;

[0070] The produced water-gas ratio of each reservoir of the water-drive gas reservoir when water invades water is compared with the abandoned water-gas ratio respectively, and the reservoir whose produced water-gas ratio when water invades water is greater than the abandoned water-gas ratio is regarded as the water-invaded layer when abandoned, and the reservoir whose produced water-gas ratio when water invades water is less than the abandoned water-gas ratio is regarded as the water-uninvaded layer when abandoned.

[0071] Optionally, the production water-gas ratio determination module includes:

[0072] A reservoir sorting unit, configured to sort each reservoir of the water drive gas reservoir in order of the permeability average from large to small or from small to large;

[0073] The production water-gas ratio determination unit is used to determine the production water-gas ratio of the reservoir when water invades water according to the order of water invasion of each reservoir in the water-drive gas reservoir, based on the basic parameters of the reservoir, the average permeability and thickness of the reservoir.

[0074] Optionally, the basic reservoir parameters include gas phase relative permeability of the reservoir at original gas saturation, water phase relative permeability of the reservoir after water flooding, reservoir temperature, original formation pressure, formation water viscosity, formation water volume coefficient, natural gas volume coefficient and natural gas PVT data;

[0075] The production water-gas ratio determining unit is used to:

[0076] establishing a first correspondence between the natural gas compressibility factor and the formation pressure and a second correspondence between the natural gas viscosity and the formation pressure based on the natural gas PVT data;

[0077] determining a natural gas compressibility factor based on the first corresponding relationship and an original formation pressure, and determining a natural gas viscosity based on the second corresponding relationship and the original formation pressure;

[0078] The water-gas ratio of the reservoir produced when water invades water is determined according to the following formula:

[0079]

[0080] Among them, H i represents the thickness of the i-th reservoir after sorting, μ w Indicates the formation water viscosity, K rwc represents the relative permeability of water phase in the reservoir after water flooding, B w represents the formation water volume coefficient, K rgi represents the gas phase relative permeability of the reservoir at the original gas saturation, μ gRepresents the viscosity of reservoir natural gas, P i represents the original formation pressure, Z represents the reservoir natural gas compression factor, P sc Indicates standard atmospheric pressure, T sc represents the standard temperature of the oil industry, T represents the reservoir temperature, B g represents the natural gas volume coefficient, α j represents the water-gas ratio of production when water invades the jth reservoir, and N represents the total number of reservoirs.

[0081] Optionally, the basic reservoir parameters include the original gas saturation of the reservoir and the residual gas saturation of the reservoir after water flooding;

[0082] The device further comprises:

[0083] a water parameter acquisition module, configured to obtain the minimum bottom hole pressure allowed for production and water parameters of the water-drive gas reservoir before calculating the water invasion sweep coefficient at abandonment based on the reservoir basic parameters, the average permeability of the water-invasion-unseen water layer at abandonment, and the thickness of each reservoir layer of the water-drive gas reservoir; wherein the water parameters of the water-drive gas reservoir include a water multiplier, a pore compressibility coefficient, and a formation water compressibility coefficient;

[0084] a limit water invasion sweep coefficient calculation module, configured to calculate the limit water invasion sweep coefficient of the water drive gas reservoir based on the water parameters of the water drive gas reservoir, the original gas saturation of the reservoir, the residual gas saturation of the reservoir after water drive, and the minimum bottom hole pressure allowed for production;

[0085] The water invasion sweep coefficient calculation module includes:

[0086] The water invasion sweep coefficient calculation unit is used to calculate the water invasion sweep coefficient at abandonment based on the water flooding gas reservoir's ultimate water invasion sweep coefficient, the reservoir basic parameters, the average permeability of the water-invasion-unseen water layer at abandonment, and the thickness of each reservoir of the water flooding gas reservoir.

[0087] Optionally, the extreme water intrusion sweep coefficient calculation module is used to:

[0088] The ultimate water invasion sweep coefficient of a water-flooded gas reservoir is calculated according to the following formula:

[0089]

[0090] Among them, N represents the water volume multiple, C w represents the formation water compressibility coefficient, C p represents the pore compressibility coefficient, P i Indicates the original formation pressure, P a Indicates the minimum bottom hole pressure allowed for production, S gi Indicates the original gas saturation of the reservoir, S grIt represents the residual gas saturation of the reservoir after water flooding, and β1 represents the ultimate water invasion coefficient of the water flooding gas reservoir.

[0091] Optionally, the water invasion sweep coefficient calculation unit is configured to:

[0092] The water intrusion coefficient during abandonment is calculated according to the following formula:

[0093]

[0094] Among them, β1 represents the ultimate water invasion sweep coefficient of the water drive gas reservoir, β2 represents the water invasion sweep coefficient when abandoned, the 1st reservoir to the nth reservoir are the water invasion visible layers when abandoned, the n+1th reservoir to the Nth reservoir are the water invasion invisible layers when abandoned, H i represents the thickness of the i-th reservoir, K i represents the mean permeability of the i-th reservoir, K n+1 , represents the mean permeability of the n+1th reservoir, and N represents the total number of reservoirs in the water-flooded gas reservoir.

[0095] The device for determining the water intrusion sweep coefficient provided in the embodiment of the present invention can execute the method for determining the water intrusion sweep coefficient provided in any embodiment of the present invention, and has the corresponding functional modules and beneficial effects of the execution method.

[0096] Example 3

[0097] Figure 8 A 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 claimed herein.

[0098] like Figure 8As shown, the electronic device 10 includes at least one processor 11 and a memory, such as a read-only memory (ROM) 12, a random access memory (RAM) 13, etc., which is communicatively connected to the at least one processor 11. The memory stores a computer program that can be executed by the at least one processor. 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 into the random access memory (RAM) 13. Various programs and data required for the operation of the electronic device 10 can also be stored in the RAM 13. The processor 11, ROM 12, and RAM 13 are connected to each other via a bus 14. An input / output (I / O) interface 15 is also connected to the bus 14.

[0099] Multiple 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 magnetic 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 via a computer network such as the Internet and / or various telecommunication networks.

[0100] The processor 11 can be any general-purpose and / or specialized processing component 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 specialized artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any other suitable processor, controller, microcontroller, etc. The processor 11 executes the various methods and processes described above, such as the method for determining the water intrusion sweep coefficient.

[0101] In some embodiments, the method for determining the water intrusion conformity coefficient may be implemented as a computer program tangibly embodied in a computer-readable storage medium, such as storage unit 18. In some embodiments, part or all of the computer program may be loaded and / or installed on electronic device 10 via ROM 12 and / or communication unit 19. When the computer program is loaded into RAM 13 and executed by processor 11, one or more steps of the method for determining the water intrusion conformity coefficient described above may be performed. Alternatively, in other embodiments, processor 11 may be configured to perform the method for determining the water intrusion conformity coefficient in any other appropriate manner (e.g., via firmware).

[0102] Various embodiments of the systems and techniques described 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), system-on-chip systems (SOCs), programmable logic devices (CPLDs), 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 interpreted on a programmable system that includes 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.

[0103] 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, such that when the computer program is executed by the processor, the functions / operations specified in the flowcharts and / or block diagrams are implemented. The computer program may 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.

[0104] In the context of the present invention, computer-readable storage media can be tangible media that can contain or store a computer program for use with an instruction execution system, device or equipment or used in combination with an instruction execution system, device or equipment. Computer-readable storage media can include but are not limited to electronic, magnetic, optical, electromagnetic, infrared or semiconductor systems, devices or equipment, or any suitable combination of the foregoing. Alternatively, computer-readable storage media can be machine-readable signal media. More specific examples of machine-readable storage media can include electrical connections based on one or more lines, portable computer disks, hard disks, random access memories (RAM), read-only memories (ROM), erasable programmable read-only memories (EPROM or flash memory), optical fibers, portable compact disk read-only memories (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination of the foregoing.

[0105] To provide interaction with a user, the systems and techniques described herein can 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 pointing device (e.g., a mouse or trackball) through which the user can provide input to the electronic device. Other types 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).

[0106] The systems and techniques described herein can be implemented in a computing system that includes back-end 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 front-end components (e.g., a user computer with a graphical user interface or 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 back-end components, middleware components, or front-end components. The components of the system can 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.

[0107] A computing system may include clients and servers. The clients and servers are typically remote from each other and typically interact via a communication network. This client-server relationship arises through computer programs running on the respective computers, creating a client-server relationship. The server may be a cloud server, also known as a cloud computing server or cloud host. This server is a hosting product within a cloud computing service ecosystem that addresses the management difficulties and limited scalability of traditional physical hosting and VPS services.

[0108] It should be understood that the various forms of the processes shown above can be used to reorder, add, or delete steps. For example, the steps described in the present invention can be performed in parallel, sequentially, or in a different order, as long as the desired results of the technical solution of the present invention can be achieved. This is not limited herein.

[0109] The above specific embodiments do not limit the scope of protection of the present invention. Those skilled in the art will appreciate that various modifications, combinations, sub-combinations, and substitutions may be made based on design requirements and other factors. Any modifications, equivalent substitutions, and improvements made within the spirit and principles of the present invention are intended to be included within the scope of protection of the present invention.

Claims

1. A method for determining a water intrusion sweep coefficient, characterized in that: include: Obtain basic reservoir parameters, geological stratification data and well logging curves of water-driven gas reservoirs; Calculating the mean permeability and thickness of each reservoir in the water-flooded gas reservoir based on the geological stratification data and the well logging curve; For each reservoir of the water-drive gas reservoir, determining the production water-gas ratio of the reservoir when water invades water according to the basic parameters of the reservoir, the average permeability and the thickness of the reservoir; Determining, from all reservoirs of the water-flooded gas reservoir, layers where water has been invaded during abandonment and layers where water has not been invaded during abandonment, based on the water-gas ratio produced during water-flooding. Calculating the water invasion sweep coefficient at abandonment based on the reservoir basic parameters, the average permeability of the water-invasion-unseen water layer at abandonment, and the thickness of each reservoir of the water-flooding gas reservoir; The basic reservoir parameters include the original gas saturation of the reservoir and the residual gas saturation of the reservoir after water flooding; Before calculating the water invasion sweep coefficient at abandonment based on the reservoir basic parameters, the average permeability of the water invasion-unseen water layer at abandonment, and the thickness of each reservoir of the water drive gas reservoir, the method further includes: Obtaining the minimum bottom hole pressure allowed for production and water parameters of the water-driven gas reservoir; wherein the water parameters of the water-driven gas reservoir include water multiple, pore compressibility coefficient, and formation water compressibility coefficient; Calculating the ultimate water invasion sweep coefficient of the water-flooding gas reservoir according to the water parameters of the water-flooding gas reservoir, the original gas saturation of the reservoir, the residual gas saturation of the reservoir after water flooding, and the minimum bottom hole pressure allowed for production; Calculating the water invasion sweep coefficient at abandonment based on the reservoir basic parameters, the average permeability of the water-invasion-unseen water layer at abandonment, and the thickness of each reservoir of the water-flooding gas reservoir includes: The water invasion sweep coefficient at abandonment is calculated based on the water flooding gas reservoir's ultimate water invasion sweep coefficient, the reservoir basic parameters, the average permeability of the unseen water layer at abandonment, and the thickness of each reservoir of the water flooding gas reservoir.

2. The method according to claim 1, characterized in that The step of determining, from all reservoirs of the water drive gas reservoir, layers where water has been invaded during abandonment and layers where water has not been invaded during abandonment based on the water-to-gas ratio of production during water invasion, comprises: Get the waste water-gas ratio; The produced water-gas ratio of each reservoir of the water-drive gas reservoir when water invades water is compared with the abandoned water-gas ratio respectively, and the reservoir whose produced water-gas ratio when water invades water is greater than the abandoned water-gas ratio is regarded as the water-invaded layer when abandoned, and the reservoir whose produced water-gas ratio when water invades water is less than the abandoned water-gas ratio is regarded as the water-uninvaded layer when abandoned.

3. The method according to claim 1, characterized in that The step of determining the production water-gas ratio of each reservoir layer of the water-drive gas reservoir when water invades water based on the reservoir basic parameters, the average permeability value, and the thickness of the reservoir layer comprises: sorting each reservoir of the water drive gas reservoir in descending or ascending order according to the permeability average; According to the order of water invasion of each reservoir in the water drive gas reservoir, the production water-gas ratio of the reservoir when water invades the reservoir is determined according to the basic parameters of the reservoir, the average permeability and the thickness of the reservoir.

4. The method according to claim 3, characterized in that The basic reservoir parameters include gas phase relative permeability of the reservoir at original gas saturation, water phase relative permeability of the reservoir after water flooding, reservoir temperature, original formation pressure, formation water viscosity, formation water volume coefficient, natural gas volume coefficient and natural gas PVT data; The method of determining the production water-gas ratio of each reservoir when water invades water according to the order of water invasion in each reservoir of the water drive gas reservoir, based on the reservoir basic parameters, the average permeability and the thickness of the reservoir, includes: establishing a first correspondence between the natural gas compressibility factor and the formation pressure and a second correspondence between the natural gas viscosity and the formation pressure based on the natural gas PVT data; determining a natural gas compressibility factor based on the first corresponding relationship and an original formation pressure, and determining a natural gas viscosity based on the second corresponding relationship and the original formation pressure; The water-gas ratio of the reservoir produced when water invades water is determined according to the following formula: Among them, H i represents the thickness of the i-th reservoir after sorting, μ w Indicates the formation water viscosity, K rwc represents the relative permeability of water phase in the reservoir after water flooding, B w represents the formation water volume coefficient, K rgi represents the gas phase relative permeability of the reservoir at the original gas saturation, μ g Represents the viscosity of reservoir natural gas, P i represents the original formation pressure, Z represents the reservoir natural gas compression factor, P sc Indicates standard atmospheric pressure, T sc represents the standard temperature of the oil industry, T represents the reservoir temperature, B g represents the natural gas volume coefficient, α j represents the water-gas ratio of production when water invades the jth reservoir, and N represents the total number of reservoirs.

5. The method according to claim 1, wherein Calculating the ultimate water invasion sweep coefficient of the water drive gas reservoir based on the water body parameters of the water drive gas reservoir, the original gas saturation of the reservoir, the residual gas saturation of the reservoir after water drive, and the minimum bottom hole pressure allowed for production includes: The ultimate water invasion sweep coefficient of a water-flooded gas reservoir is calculated according to the following formula: Among them, N represents the water volume multiple, C w represents the formation water compressibility coefficient, C p represents the pore compressibility coefficient, P i Indicates the original formation pressure, P a Indicates the minimum bottom hole pressure allowed for production, S gi Indicates the original gas saturation of the reservoir, S gr It represents the residual gas saturation of the reservoir after water flooding, and β1 represents the ultimate water invasion coefficient of the water flooding gas reservoir.

6. The method according to claim 1, characterized in that The water invasion sweep coefficient at abandonment is calculated based on the extreme water invasion sweep coefficient of the water-flooding gas reservoir, the basic reservoir parameters, the average permeability of the water-invasion-unseen water layer at abandonment, and the thickness of each reservoir of the water-flooding gas reservoir, including: The water intrusion coefficient during abandonment is calculated according to the following formula: Among them, β1 represents the ultimate water invasion sweep coefficient of the water drive gas reservoir, β2 represents the water invasion sweep coefficient when abandoned, the 1st reservoir to the nth reservoir are the water invasion visible layers when abandoned, the n+1th reservoir to the Nth reservoir are the water invasion invisible layers when abandoned, H i represents the thickness of the i-th reservoir, K i represents the mean permeability of the i-th reservoir, K n+1 , represents the mean permeability of the n+1th reservoir, and N represents the total number of reservoirs in the water-flooded gas reservoir.

7. A device for determining a water invasion sweep coefficient, characterized in that: include: Reservoir parameter acquisition module, used to obtain basic reservoir parameters, geological stratification data and logging curves of water-driven gas reservoirs; a permeability mean calculation module, configured to calculate the permeability mean and thickness of each reservoir of the water-flooded gas reservoir based on the geological stratification data and the well logging curve; a production water-gas ratio determination module, configured to determine, for each reservoir of the water-drive gas reservoir, the production water-gas ratio of the reservoir when water invades water according to the basic parameters of the reservoir, the average permeability value and the thickness of the reservoir; a reservoir classification module for determining, from all reservoirs of the water-flooding gas reservoir, layers where water was seen during abandonment and layers where water was not seen during abandonment, based on the water-gas ratio of production during water-flooding; a water invasion sweep coefficient calculation module, configured to calculate the water invasion sweep coefficient at abandonment based on the reservoir basic parameters, the average permeability of the water-invasion-unseen water layer at abandonment, and the thickness of each reservoir of the water-flooding gas reservoir; The basic reservoir parameters include the original gas saturation of the reservoir and the residual gas saturation of the reservoir after water flooding; The device further comprises: a water parameter acquisition module, configured to obtain the minimum bottom hole pressure allowed for production and water parameters of the water-drive gas reservoir before calculating the water invasion sweep coefficient at abandonment based on the reservoir basic parameters, the average permeability of the water-invasion-unseen water layer at abandonment, and the thickness of each reservoir layer of the water-drive gas reservoir; wherein the water parameters of the water-drive gas reservoir include a water multiplier, a pore compressibility coefficient, and a formation water compressibility coefficient; a limit water invasion sweep coefficient calculation module, configured to calculate the limit water invasion sweep coefficient of the water drive gas reservoir based on the water parameters of the water drive gas reservoir, the original gas saturation of the reservoir, the residual gas saturation of the reservoir after water drive, and the minimum bottom hole pressure allowed for production; The water invasion sweep coefficient calculation module includes: The water invasion sweep coefficient calculation unit is used to calculate the water invasion sweep coefficient at abandonment based on the water flooding gas reservoir's ultimate water invasion sweep coefficient, the reservoir basic parameters, the average permeability of the water-invasion-unseen water layer at abandonment, and the thickness of each reservoir of the water flooding gas reservoir.

8. 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. The computer program is executed by the at least one processor so that the at least one processor can perform the method for determining the water invasion sweep coefficient according to any one of claims 1 to 6.

9. 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 the water intrusion sweep coefficient according to any one of claims 1 to 6 when executed.

Images