A method and apparatus for determining the location of intensified wells based on old well water ridge waves and models.

Abstract

This invention provides a method and apparatus for determining infill well locations based on an old well water ridge sweep tracing model. The method includes acquiring reservoir physical parameters, fluid parameters, and production dynamic data of the target horizontal well's production layer; obtaining relative permeability curves and oil displacement efficiency through core reactivity experiments; calculating the water ridge height during sidetracking; substituting the water ridge height, reservoir physical parameters, fluid parameters, and production dynamic data into the horizontal well oil-water interface description formula, normalizing the formula to make it dimensionless, and calculating the correspondence between dimensionless cumulative oil production, dimensionless water ridge width, and dimensionless water ridge height using the oil displacement efficiency in conjunction with water-drive theory; establishing a data map comparison table of dimensionless cumulative oil production, dimensionless water ridge width, and dimensionless water ridge height, and determining the limiting well spacing of the infill wells based on the data map comparison table; and performing infill or re-branching of the target horizontal well based on the limiting well spacing.

Description

Technical Field

[0001] This invention relates to the field of oil and gas extraction technology, and in particular to a method and apparatus for determining in-depth well locations based on old well water ridge wave pattern. Background Technology

[0002] Currently, most marine sandstone formations in the eastern South China Sea are in the mid-to-late stages of oilfield development, at a "high to ultra-high" water cut level, making continuous tapping of potential difficult. Many older wells are horizontal wells, and their large-scale use in the eastern South China Sea has ensured the rapid and efficient development of oilfields, highlighting the significant advantages of their extensive deployment and efficient development. However, this also presents challenges in predicting the impact range of horizontal wells and the remaining oil around them, especially for high water-cut reservoirs. Therefore, designing infill wells or sidetracking branches of older wells to achieve optimal development results and economic efficiency is crucial for the economic tapping of remaining oil in offshore oilfields. To avoid rapid water flooding after the commissioning of infill wells or sidetracking branches of older wells, accurate quantitative description of the water ridge morphology and extent of older wells is essential during the design phase.

[0003] However, existing methods for analyzing the sweep range of water ridges in horizontal wells are mostly based on empirical formulas for numerical model calibration and analysis of experimental data. They are often based on reservoir-level regularities. Given the heterogeneity of actual reservoirs and the large differences in production between different wells, there are limitations in accurately and quantitatively evaluating water ridges in actual single wells. There is still a lack of targeted inversion of water ridges for specific old wells. Summary of the Invention

[0004] The purpose of this invention is to provide a method and apparatus for determining the location of in-depth wells based on the old well water ridge wave pattern, in order to solve the problem mentioned in the background art that the existing horizontal well water ridge wave pattern range analysis methods are mostly based on numerical model calibration empirical formulas and experimental data analysis, which are often based on reservoir level regularity and have limitations.

[0005] To achieve the above objectives, the present invention provides the following technical solution: a method for determining the location of infill wells based on an old well water ridge sweep tracing model, comprising the following steps: acquiring reservoir physical parameters, fluid parameters, and production dynamic data of the production layer of a target horizontal well; obtaining the relative permeability curve and oil displacement efficiency through core relative permeability experiments; determining the reservoir physical parameters corresponding to sidetracking based on the relative permeability curve, and calculating the water ridge height during sidetracking by combining the fluid parameters and the production dynamic data; and combining the water ridge height, the reservoir physical parameters, and the fluid parameters with the production dynamic data. The body parameters and the production dynamic data are substituted into the horizontal well oil-water interface description formula. The formula is normalized to be dimensionless, and the correspondence between dimensionless cumulative oil production, dimensionless water ridge width, and dimensionless water ridge height is calculated by combining the water drive theory and the oil displacement efficiency. A data map comparison table of the dimensionless cumulative oil production, dimensionless water ridge width, and dimensionless water ridge height is established, and the limit well spacing of the infill well is determined according to the data map comparison table. The target horizontal well is infilled or re-branched according to the limit well spacing.

[0006] Optionally, the reservoir physical properties include permeability, porosity, water saturation, and volume factor; the fluid parameters include fluid density and viscosity; and the production dynamic data include daily oil production, daily fluid production, water cut, bottom hole flowing pressure, and initial oil column height.

[0007] Optionally, the step of obtaining the relative permeability curve and oil displacement efficiency through core relative permeability experiments includes: the oil displacement efficiency is calculated using the initial water saturation and residual oil saturation obtained from the core experiment, and the calculation formula is as follows: Among them, E d For oil displacement efficiency, S wi S represents the initial water saturation. or This represents the residual oil saturation.

[0008] Optionally, the step of obtaining the relative permeability curve and oil displacement efficiency through core relative permeability experiments further includes: correcting the oil displacement efficiency using an approximate theoretical water drive curve, wherein the correction formula is as follows: in S represents the average water saturation. wi S represents the initial water saturation. or denoted as residual oil saturation, and w and b are parameters calculated from approximate theoretical water drive curves.

[0009] Optionally, the step of determining the reservoir physical property parameters corresponding to sidetracking based on the relative permeability curve, and calculating the water ridge height during sidetracking in conjunction with the fluid parameters and the production dynamic data, includes: determining the reservoir physical property parameters corresponding to sidetracking based on the relative permeability curve, and calculating the dimensionless fluid production index during sidetracking, using the following formula: JDL K is the dimensionless fluid collection index. ro K represents the oil phase permeability. rw For water phase permeability, μ o The viscosity of crude oil is μ. w B represents the viscosity of water. o B is the crude oil volume factor. w The water volume factor is used; the fluid production index during side-drilling is calculated based on the aforementioned production dynamic data and the dimensionless fluid production index, using the formula ΔP. initial =P res-initial -bhp initial , J = J initial ×J DL , where ΔP initial P represents the initial pressure difference of the reservoir. res-initial bhp is the initial pressure of the reservoir. initial For the initial wellbore flowing pressure during the initial production phase, J initial The production index is q, where q is the flow rate and J is the production index during sidetracking. DL The fluid recovery index is a dimensionless index; the decrease in reservoir pressure compared to the initial production level is calculated based on the production dynamic data and the fluid recovery index, using the following formula: P res = bhp + ΔP, where ΔP is the decrease in reservoir pressure compared to the initial production level, q is the flow rate data, J is the fluid production index during sidetracking, and P... res Here, is the reservoir pressure during sidetracking, and bhp is the bottomhole flowing pressure. The water ridge height during sidetracking is calculated based on the production dynamics data, the fluid parameters, and the decrease in reservoir pressure compared to the initial production level. The calculation formula is as follows: Among them, y G The water ridge height during side drilling, y G-initial P is the initial oil column height. res-initial P is the initial pressure of the reservoir. res ρ is the reservoir pressure during side-drilling, ρ is the density of formation water, and g is the acceleration due to gravity.

[0010] Optionally, the step of substituting the water ridge height, the reservoir physical properties, the fluid parameters, and the production dynamic data into the horizontal well oil-water interface description formula includes: in the horizontal well oil-water interface description formula, the... Defined as the far-well region, the oil-water interface within the far-well region can be represented as: Where x is the x-coordinate of any point at the oil-water interface, and y is the y-coordinate of any point at the oil-water interface. G D is the water ridge height, and D is a control parameter. coneWhere is the width of the water ridge, μ0 is the viscosity of the crude oil, q is the flow rate, and K is the flow rate. H Let ρ be the planar permeability, ρ be the density of formation water, and g be the acceleration due to gravity; in the formula describing the oil-water interface in the horizontal well, Defined as the near-wellbore region, the oil-water interface within the near-wellbore region can be represented as: Where x is the x-coordinate of any point at the oil-water interface, and y is the y-coordinate of any point at the oil-water interface. G Where is the water ridge height, D is the control parameter, μ0 is the crude oil viscosity, q is the flow rate data, and K is the flow rate. H ρ is the planar permeability, ρ is the density of the formation water, and g is the gravitational acceleration.

[0011] Optionally, the step of normalizing the formula to make it dimensionless includes: normalizing the formula to a dimensionless form. Where x′ is the dimensionless abscissa of any point at the oil-water interface, and y′ is the dimensionless ordinate of any point at the oil-water interface. G ′ is the dimensionless water ridge height, D′ is the dimensionless water ridge width, x is the x-coordinate of any point at the oil-water interface, y is the y-coordinate of any point at the oil-water interface, and K H ρ is the planar permeability, ρ is the density of formation water, g is the acceleration due to gravity, and μ is the groundwater density. o Let q represent crude oil viscosity and q represent flow rate data; the relationship between the dimensionless x′ and dimensionless y′ of any point on the oil-water interface in the far-well region is as follows: The relationship between the dimensionless x′ and the dimensionless y′ of any point on the oil-water interface in the near-wellbore region is as follows: Where x′ is the dimensionless abscissa of any point at the oil-water interface, and y′ is the dimensionless ordinate of any point at the oil-water interface. G ′ represents the dimensionless height of the water ridge, and D′ represents the dimensionless width of the water ridge.

[0012] Optionally, the step of calculating the correspondence between the dimensionless cumulative oil production and the dimensionless water ridge width and dimensionless water ridge height by combining water-drive oil displacement theory includes: expressing the cumulative oil production as N using water-drive oil displacement theory. p =VE d , where N P V represents the cumulative oil production of an old well at a certain moment, where V is the water drive volume and E is the total oil production. d The oil displacement efficiency; its dimensionless form is N. p ′=V′E d , Among them, N' p V' represents the cumulative dimensionless oil production of an old well at a certain moment, and E represents the dimensionless water drive volume. dFor oil displacement efficiency, K H ρ is the planar permeability, ρ is the density of formation water, g is the acceleration due to gravity, and μ is the groundwater density. o Let q represent the crude oil viscosity and q represent the flow rate. The horizontal cross-section of the water ridge in the bottom water reservoir is rationally divided; the cross-section can be considered as a combination of two semicircles and a rectangle. The final volume can be obtained by combining the integral of the dimensionless oil-water interface morphology equation and the formulas for the volumes of solids of revolution and cylinders.

[0013] Where L′ is the dimensionless well length of the old well, V′ is the dimensionless water drive volume, and y G ′ is the dimensionless water ridge height, D′ is the dimensionless water ridge width, x′ is the dimensionless abscissa of any point at the oil-water interface, dx is the derivative, and K H Let ρ be the planar permeability, g be the density of formation water, μ0 be the acceleration due to gravity, q be the flow rate, and L be the length of the old well. By substituting the formulas for calculating the dimensionless water drive volume and oil displacement efficiency into the formula for calculating the cumulative dimensionless oil production of the old well at a certain moment, we can obtain the relationship between the dimensionless cumulative oil production and the dimensionless water ridge width and dimensionless water ridge height.

[0014] Optionally, the step of determining the limit well spacing of the densified wells according to the data map comparison table includes: taking 1 / 2 of the width of the water ridge as the limit well spacing of the densified new wells.

[0015] On the other hand, the present invention also provides a device for determining the location of infiltrated wells based on an old well water ridge sweep tracing model, comprising: an acquisition module for acquiring reservoir physical property parameters, fluid parameters, and production dynamic data of the production layer of a target horizontal well; a core relative permeability test module for obtaining relative permeability curves and oil displacement efficiency through core relative permeability tests; a water ridge height acquisition module for determining the reservoir physical property parameters corresponding to sidetracking based on the relative permeability curves, and calculating the water ridge height during sidetracking by combining the fluid parameters and the production dynamic data; and an oil-water interface description formula normalization module for normalizing the water ridge height, The reservoir physical parameters, the fluid parameters, and the production dynamic data are substituted into the horizontal well oil-water interface description formula. The formula is normalized to be dimensionless, and the correspondence between dimensionless cumulative oil production, dimensionless water ridge width, and dimensionless water ridge height is calculated using the oil displacement efficiency in conjunction with waterflooding theory. The data mapping module is used to establish a data mapping comparison table of dimensionless cumulative oil production, dimensionless water ridge width, and dimensionless water ridge height, and to determine the limit well spacing of the infill wells based on the data mapping comparison table. The infill module is used to infill or regenerate branches of the target horizontal well based on the limit well spacing.

[0016] Compared with the prior art, the beneficial effects of the present invention are:

[0017] By acquiring parameter information and dynamic data, the relative permeability curve and oil displacement efficiency are obtained through core relative permeability experiments. Then, the water ridge height during sidetracking is calculated. Based on the horizontal well oil-water interface description formula, and through normalization deformation combined with water-drive theory, the relationship between dimensionless cumulative oil production and two variables, dimensionless water ridge height and dimensionless water ridge width, is derived. Finally, the limit well spacing of the infill well is determined using the established data map, and the target horizontal well is infilled or re-branched according to the limit well spacing. This method overcomes the limitations of conventional mainstream water ridge analysis methods, such as narrow applicability and inability to provide accurate quantitative evaluation results for specific production wells. It is a more reasonable and effective method for analyzing the affected area of ​​old wells and designing and evaluating re-branching. Attached Figure Description

[0018] Figure 1 This is a schematic diagram of the method steps of the present invention.

[0019] Figure 2 This is a schematic diagram of the technical process of the present invention.

[0020] Figure 3 This is a schematic diagram of the dimensionless oil-water interface morphology of the present invention.

[0021] Figure 4 This is a schematic diagram of the splitting of the bottom surface of the water ridge in this invention.

[0022] Figure 5 This is a schematic diagram of the dimensionless water drive volume of the present invention.

[0023] Figure 6 This is a schematic diagram showing the calculation results of the water ridge width of the old well and the optimized location of the new well in an example of the present invention. Detailed Implementation

[0024] The present invention will now be clearly and completely described in conjunction with the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments.

[0025] It should be noted that the terms "first," "second," etc., in the specification, claims, and accompanying drawings of this application are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate for the embodiments of this application described herein. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover non-exclusive inclusion; for example, a process, method, system, product, or apparatus that comprises a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, methods, products, or apparatus.

[0026] Those skilled in the art will understand that, unless explicitly stated otherwise, the singular forms “a,” “an,” “the,” and “the” used herein may also include the plural forms. It should be further understood that the term “comprising” as used in the specification of this application means the presence of features, integers, steps, operations, elements, and / or components, but does not exclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and / or groups thereof. It should be understood that when we say an element is “connected” or “coupled” to another element, it can be directly connected or coupled to the other element, or there may be intermediate elements. Furthermore, “connected” or “coupled” as used herein can include wireless connections or wireless coupling. The term “and / or” as used herein includes all or any units and all combinations of one or more associated listed items.

[0027] It will be understood by those skilled in the art that, unless otherwise defined, all terms used herein (including technical and scientific terms) have the same meaning as commonly understood by one of ordinary skill in the art to which this application pertains. It should also be understood that terms such as those defined in general dictionaries should be understood to have the same meaning as in the context of the prior art, and should not be interpreted in an idealized or overly formal sense unless specifically defined as herein.

[0028] It should be understood that the sequence number and size of each step in this embodiment do not imply the order of execution. The execution order of each process is determined by its function and internal logic, and should not constitute any limitation on the implementation process of this application embodiment.

[0029] It should be noted that, unless otherwise specified, the embodiments and features described in this application can be combined with each other. This application will now be described in detail with reference to the accompanying drawings and embodiments.

[0030] Please refer to Figure 1 and Figure 2 The present invention provides a method for determining the location of infiltrated wells based on an old well water ridge wavelet model, comprising the following steps:

[0031] S100: Obtain reservoir physical property parameters, fluid parameters, and production dynamic data of the target horizontal well's production layer.

[0032] Specifically, the reservoir physical parameters include oil layer thickness, water-avoidance height (i.e., the distance between the oil-water interface and the horizontal well), porosity, and permeability, which are obtained from development data such as oil well reservoir design documents; the fluid parameters include oil / water density, oil / water viscosity, and crude oil volume coefficient, which are obtained from development data such as fluid physical property analysis reports; and the production dynamic data includes the development date and the corresponding daily fluid production, daily water production, and water cut, which are exported from the oil well production dynamic data acquisition system.

[0033] S200, the relative permeability curve and oil displacement efficiency were obtained through core relative permeability experiments.

[0034] Specifically, the relative permeability curve is the relationship curve between the relative permeability of the oil phase and the water phase and the water saturation, and is called the oil-water two-phase relative permeability curve. As the water saturation increases, the relative permeability of the oil phase decreases and the relative permeability of the water phase increases. The oil displacement efficiency can be calculated by obtaining the initial water saturation and residual oil saturation through core tests.

[0035] S300. Determine the reservoir physical property parameters corresponding to the side-drilling based on the relative permeability curve, and calculate the water ridge height during side-drilling by combining the fluid parameters and the production dynamic data.

[0036] Specifically, taking a horizontal well in marine sandstone in the eastern South China Sea as an example, the following parameters are obtained: reservoir physical properties such as permeability, porosity, and water saturation; fluid parameters such as fluid density and viscosity; and production dynamic data such as daily oil production, daily fluid production, water cut, and bottom hole flowing pressure of the old well. Relative permeability curves are obtained through core permeability experiments, and the dimensionless fluid production index at the corresponding water cut is calculated during MRC sidetracking. The fluid production index at the initial production stage is calculated using the initial reservoir pressure and the bottom hole flowing pressure and flow rate data of the old well during the initial production stage. The fluid production index during sidetracking is then calculated in conjunction with the dimensionless fluid production index. Based on the flow rate, bottom hole flowing pressure, and fluid production index of the old well before sidetracking, the reservoir pressure during sidetracking is calculated, and then the decrease in reservoir pressure is compared to the decrease during production. The water ridge height during MRC sidetracking is calculated based on the initial oil column height of the old well before sidetracking, the decrease in reservoir pressure, and the density and gravity acceleration of the formation water.

[0037] S400. Substitute the water ridge height, the reservoir physical property parameters, the fluid parameters, and the production dynamic data into the horizontal well oil-water interface description formula. Normalize the formula to make it dimensionless. Combine the water drive theory with the oil displacement efficiency to calculate the correspondence between the dimensionless cumulative oil production and the dimensionless water ridge width and dimensionless water ridge height.

[0038] Specifically, by combining the cumulative oil production of old wells during MRC sidetracking, the dimensionless cumulative oil production and dimensionless oil column height are calculated. The oil displacement efficiency can be calculated by obtaining the initial water saturation and residual oil saturation through core tests. Substituting these values ​​into the aforementioned calculation formula, the relationship between dimensionless cumulative oil production and dimensionless water ridge width and dimensionless water ridge height can be obtained.

[0039] S500. Establish a data map comparison table of dimensionless cumulative oil production, dimensionless water ridge width, and dimensionless water ridge height, and determine the limit well spacing of the infill wells based on the data map comparison table.

[0040] S600. Based on the extreme well spacing, densify or regenerate branches of the target horizontal well.

[0041] Specifically, taking well 20H2 in a certain oil reservoir as an example, the initial water saturation was 0.2, the residual oil saturation was 0.26, and the calculated oil displacement efficiency of the reservoir was 0.56. Based on dynamic data, the water ridge height of well 20H2 during MRC sidetracking was calculated to be 10 meters, with a cumulative oil production of 146,000 cubic meters. The column corresponding to the initial oil column height was found in the reference table. The dimensionless water ridge width was determined based on the dimensionless cumulative oil production, and finally converted to the actual water ridge width, which was calculated to be approximately 125 meters. The area within half the water ridge width is the region most affected by the old well's water ridge. Taking the new branch at the old branch water ridge... The length within the width of / 2 is less than the total length * 0.25 as the limit. The length of the reconstructed branch within the large influence area is calculated based on the reconstructed branch length and the sidetracking angle, and this is used as the basis for evaluating whether the reconstructed branch design is reasonable. In this example, after final calculation, the design sidetracking branch lengths for the 20H2 well MRC are 548 meters and 743 meters, respectively. A sidetracking angle greater than 15° can achieve a better sidetracking effect. If it is necessary to densify the well network, the reasonable well spacing range should be greater than 1 / 2 the width of the old tributary ridge, i.e., 62.5 meters, so as to achieve the optimization evaluation of the location of the densified well or the MRC reconstructed branch well.

[0042] Understandably, by acquiring parameter information and dynamic data, obtaining the relative permeability curve and oil displacement efficiency through core relative permeability experiments, and then calculating the water ridge height during sidetracking, based on the horizontal well oil-water interface description formula, and through normalization deformation combined with water-drive theory, the relationship between dimensionless cumulative oil production and two variables—dimensional water ridge height and dimensionless water ridge width—is derived. Finally, the established data map is used to determine the limit well spacing of infill wells, and the target horizontal well is infilled or re-branched according to the limit well spacing. This method overcomes the limitations of conventional mainstream water ridge analysis methods, such as their narrow applicability and inability to provide accurate quantitative evaluation results for specific production wells. It is a more reasonable and effective method for old well sweep range analysis and re-branch design evaluation. This method can achieve accurate assessment of the water ridge sweep range of old wells in inter-well infilling and old well sidetracking MRC design, while guiding the optimization of infill well locations. It has good application prospects for the continuous tapping of potential in old oilfields in the high water-cut stage.

[0043] Optionally, the reservoir physical properties include permeability, porosity, water saturation, and volume factor; the fluid parameters include fluid density and viscosity; and the production dynamic data include daily oil production, daily fluid production, water cut, bottom hole flowing pressure, and initial oil column height.

[0044] Optionally, the step of obtaining the relative permeability curve and oil displacement efficiency through core relative permeability experiments includes: the oil displacement efficiency is calculated using the initial water saturation and residual oil saturation obtained from the core experiment, and the calculation formula is as follows: Among them, E d For oil displacement efficiency, S wi S represents the initial water saturation. or This represents the residual oil saturation.

[0045] Optionally, the step of obtaining the relative permeability curve and oil displacement efficiency through core relative permeability experiments further includes: correcting the oil displacement efficiency using an approximate theoretical water drive curve, wherein the correction formula is as follows: in S represents the average water saturation. wi S represents the initial water saturation. or denoted as residual oil saturation, and w and b are parameters calculated from approximate theoretical water drive curves.

[0046] Please refer to Figure 3 Optionally, the step of determining the reservoir physical property parameters corresponding to sidetracking based on the relative permeability curve, and calculating the water ridge height during sidetracking in conjunction with the fluid parameters and the production dynamic data, includes: determining the reservoir physical property parameters corresponding to sidetracking based on the relative permeability curve, and calculating the dimensionless fluid production index during sidetracking, using the following formula: J DL K is the dimensionless fluid collection index. ro K represents the oil phase permeability. rw For water phase permeability, μ o The viscosity of crude oil is μ. w B represents the viscosity of water. o B is the crude oil volume factor. w The water volume factor is used; the fluid production index during side-drilling is calculated based on the aforementioned production dynamic data and the dimensionless fluid production index, using the formula ΔP. initial =P res-initial -bhp initial , J = J initial ×J DL , where ΔP initial P represents the initial pressure difference of the reservoir. res-initial bhp is the initial pressure of the reservoir. initial For the initial wellbore flowing pressure during the initial production phase, J initial The production index is q, where q is the flow rate and J is the production index during sidetracking. DL The fluid recovery index is a dimensionless index; the decrease in reservoir pressure compared to the initial production level is calculated based on the production dynamic data and the fluid recovery index, using the following formula: P res = bhp + ΔP, where ΔP is the decrease in reservoir pressure compared to the initial production level, q is the flow rate data, J is the fluid production index during sidetracking, and P... resHere, is the reservoir pressure during sidetracking, and bhp is the bottomhole flowing pressure. The water ridge height during sidetracking is calculated based on the production dynamics data, the fluid parameters, and the decrease in reservoir pressure compared to the initial production level. The calculation formula is as follows: Among them, y G The water ridge height during side drilling, y G-initial P is the initial oil column height. res-initial P is the initial pressure of the reservoir. res ρ is the reservoir pressure during side-drilling, ρ is the density of formation water, and g is the acceleration due to gravity.

[0047] Optionally, the step of substituting the water ridge height, the reservoir physical properties, the fluid parameters, and the production dynamic data into the horizontal well oil-water interface description formula includes: in the horizontal well oil-water interface description formula, the... Defined as the far-well region, the oil-water interface within the far-well region can be represented as: Where x is the x-coordinate of any point at the oil-water interface, and y is the y-coordinate of any point at the oil-water interface. G D is the water ridge height, and D is a control parameter. cone Where is the width of the water ridge, μ0 is the viscosity of the crude oil, q is the flow rate, and K is the flow rate. H Let ρ be the planar permeability, ρ be the density of formation water, and g be the acceleration due to gravity; in the formula describing the oil-water interface in the horizontal well, Defined as the near-wellbore region, the oil-water interface within the near-wellbore region can be represented as: Where x is the x-coordinate of any point at the oil-water interface, and y is the y-coordinate of any point at the oil-water interface. G Where is the water ridge height, D is the control parameter, μ0 is the crude oil viscosity, q is the flow rate data, and K is the flow rate. H ρ is the planar permeability, ρ is the density of the formation water, and g is the gravitational acceleration.

[0048] Optionally, the step of normalizing the formula to make it dimensionless includes: normalizing the formula to a dimensionless form. Where x′ is the dimensionless abscissa of any point at the oil-water interface, and y′ is the dimensionless ordinate of any point at the oil-water interface. G ′ is the dimensionless water ridge height, D′ is the dimensionless water ridge width, x is the x-coordinate of any point at the oil-water interface, y is the y-coordinate of any point at the oil-water interface, and K H ρ is the planar permeability, ρ is the density of formation water, g is the acceleration due to gravity, and μ is the groundwater density. o Let q represent crude oil viscosity and q represent flow rate data; the relationship between the dimensionless x′ and dimensionless y′ of any point on the oil-water interface in the far-well region is as follows: The relationship between the dimensionless x′ and the dimensionless y′ of any point on the oil-water interface in the near-wellbore region is as follows: Where x′ is the dimensionless abscissa of any point at the oil-water interface, and y′ is the dimensionless ordinate of any point at the oil-water interface. G ′ represents the dimensionless height of the water ridge, and D′ represents the dimensionless width of the water ridge.

[0049] Please refer to Figure 4 and Figure 5 Optionally, the step of calculating the correspondence between dimensionless cumulative oil production and dimensionless water ridge width and dimensionless water ridge height by combining water-drive oil displacement theory includes: expressing the cumulative oil production as N using water-drive oil displacement theory. p =VE d , where N P V represents the cumulative oil production of an old well at a certain moment, where V is the water drive volume and E is the total oil production. d The oil displacement efficiency; its dimensionless form is N. p ′=V′E d , Among them, N' p V' represents the cumulative dimensionless oil production of an old well at a certain moment, and E represents the dimensionless water drive volume. d For oil displacement efficiency, K H ρ is the planar permeability, ρ is the density of formation water, g is the acceleration due to gravity, and μ is the groundwater density. o Let q represent the crude oil viscosity and q represent the flow rate. The horizontal cross-section of the water ridge in the bottom water reservoir is rationally divided; the cross-section can be considered as a combination of two semicircles and a rectangle. The final volume can be obtained by combining the integral of the dimensionless oil-water interface morphology equation and the formulas for the volumes of solids of revolution and cylinders. Where L′ is the dimensionless well length of the old well, V′ is the dimensionless water drive volume, and y G ′ is the dimensionless water ridge height, D′ is the dimensionless water ridge width, x′ is the dimensionless abscissa of any point at the oil-water interface, dx is the derivative, and K H Let ρ be the planar permeability, g be the density of formation water, μ0 be the acceleration due to gravity, q be the flow rate, and L be the length of the old well. By substituting the formulas for calculating the dimensionless water drive volume and oil displacement efficiency into the formula for calculating the cumulative dimensionless oil production of the old well at a certain moment, we can obtain the relationship between the dimensionless cumulative oil production and the dimensionless water ridge width and dimensionless water ridge height.

[0050] Please refer to Figure 6Optionally, the step of determining the limit well spacing of the densification wells according to the data map comparison table includes: taking 1 / 2 of the width of the water ridge as the limit well spacing of the densification wells.

[0051] Specifically, the data map comparison table between dimensionless cumulative oil production and dimensionless water ridge width and dimensionless water ridge height is shown in Table 1 below:

[0052] Table 1: Data graph comparison table.

[0053]

[0054]

[0055] On the other hand, the present invention also provides a device for determining the location of infiltrated wells based on an old well water ridge sweep tracing model, comprising: an acquisition module for acquiring reservoir physical property parameters, fluid parameters, and production dynamic data of the production layer of a target horizontal well; a core relative permeability test module for obtaining relative permeability curves and oil displacement efficiency through core relative permeability tests; a water ridge height acquisition module for determining the reservoir physical property parameters corresponding to sidetracking based on the relative permeability curves, and calculating the water ridge height during sidetracking by combining the fluid parameters and the production dynamic data; and an oil-water interface description formula normalization module for normalizing the water ridge height, The reservoir physical parameters, the fluid parameters, and the production dynamic data are substituted into the horizontal well oil-water interface description formula. The formula is normalized to be dimensionless, and the correspondence between dimensionless cumulative oil production, dimensionless water ridge width, and dimensionless water ridge height is calculated using the oil displacement efficiency in conjunction with waterflooding theory. The data mapping module is used to establish a data mapping comparison table of dimensionless cumulative oil production, dimensionless water ridge width, and dimensionless water ridge height, and to determine the limit well spacing of the infill wells based on the data mapping comparison table. The infill module is used to infill or regenerate branches of the target horizontal well based on the limit well spacing.

[0056] On the other hand, the present invention also provides a computer device, including a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to implement the steps of the above-described control method.

[0057] On the other hand, the present invention also provides a computer-readable storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the steps of the above-described control method.

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

[0059] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above methods. Any references to memory, storage, database or other media used in the embodiments provided by the present invention can include non-volatile and / or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM) or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in various forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), dual data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), RAMbus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and RAMbus dynamic RAM (RDRAM), etc.

[0060] The above are merely embodiments of the present invention and do not limit the patent scope of the present invention. Any equivalent modifications made based on the content of the present invention's specification and drawings, or direct or indirect applications in related technical fields, are similarly included within the patent protection scope of the present invention.

Claims

1. A method for determining the location of intensified wells based on an old well water ridge wave pattern, characterized by the following steps: include: Acquire reservoir physical properties and fluid parameters of the production layer of the target horizontal well, as well as the production dynamic data of the target horizontal well; The relative permeability curves and oil displacement efficiency were obtained through core relative permeability experiments. The reservoir physical parameters corresponding to sidetracking are determined based on the relative permeability curve, and the water ridge height during sidetracking is calculated by combining the fluid parameters and the production dynamic data. The water ridge height, reservoir physical parameters, fluid parameters, and production dynamic data are substituted into the horizontal well oil-water interface description formula. The formula is normalized and transformed to be dimensionless. The correspondence between dimensionless cumulative oil production and dimensionless water ridge width and dimensionless water ridge height is obtained by calculating the oil displacement efficiency in combination with water drive theory. Establish a data map comparison table of the dimensionless cumulative oil production, the dimensionless water ridge width, and the dimensionless water ridge height, and determine the limit well spacing of the infill wells based on the data map comparison table; The target horizontal well is densified or re-branched based on the specified limit well spacing.

2. The method for determining the location of infiltrated wells based on the old well water ridge wave pattern according to claim 1, characterized in that, The reservoir physical properties include permeability, porosity, water saturation, and volume factor; the fluid parameters include fluid density and viscosity; and the production dynamic data include daily oil production, daily liquid production, water cut, bottom hole flowing pressure, and initial oil column height.

3. The method for determining the location of infiltrated wells based on the old well water ridge wavelet model according to claim 2, characterized in that, The steps for obtaining the relative permeability curve and oil displacement efficiency through core relative permeability experiments include: The oil displacement efficiency was calculated using the initial water saturation and residual oil saturation obtained from core tests, and the calculation formula is as follows: Among them, E d For oil displacement efficiency, S wi S represents the initial water saturation. or This represents the residual oil saturation.

4. The method for determining the location of infiltrated wells based on the old well water ridge wave pattern according to claim 3, characterized in that, The step of obtaining the relative permeability curve and oil displacement efficiency through core relative permeability experiments also includes: The oil displacement efficiency is corrected using an approximate theoretical water drive curve, and the correction formula is as follows: in S represents the average water saturation. wi S represents the initial water saturation. or denoted as residual oil saturation, and w and b are parameters calculated from approximate theoretical water drive curves.

5. The method for determining the location of infiltrated wells based on the old well water ridge wave pattern according to claim 2, characterized in that, The step of determining the reservoir physical property parameters corresponding to sidetracking based on the relative permeability curve, and calculating the water ridge height during sidetracking in conjunction with the fluid parameters and the production dynamic data includes: The reservoir physical parameters corresponding to sidetracking are determined based on the relative permeability curve, and the dimensionless fluid production index during sidetracking is calculated using the following formula: J DL K is the dimensionless fluid collection index. ro K represents the oil phase permeability. rw For water phase permeability, μ o The viscosity of crude oil is μ. w B represents the viscosity of water. o B is the crude oil volume factor. w The volume coefficient is the water volume factor. The fluid production index during sidetracking is calculated based on the aforementioned production dynamic data and the dimensionless fluid production index. The calculation formula is ΔP. initial =P res-initial -bhp initial , J = J initial ×J DL , where ΔP initial P represents the initial pressure difference of the reservoir. res-initial bhp is the initial pressure of the reservoir. initial For the initial wellbore flowing pressure during the initial production phase, J initial The production index is q, where q is the flow rate and J is the production index during sidetracking. DL The dimensionless fluid collection index; The decrease in reservoir pressure compared to the initial production level was calculated based on the aforementioned production dynamic data and the fluid recovery index. The calculation formula is as follows: P res = bhp + ΔP, where ΔP is the decrease in reservoir pressure compared to the initial production level, q is the flow rate data, J is the fluid production index during sidetracking, and P... res is the reservoir pressure during sidetracking, and bhp is the bottom hole flowing pressure. Based on the production dynamic data, the fluid parameters, and the reservoir pressure drop compared to the initial production level, the water ridge height during side-drilling is calculated using the following formula: Among them, y G The water ridge height during side drilling, y G-initial P is the initial oil column height. res-initial P is the initial pressure of the reservoir. res ρ is the reservoir pressure during side-drilling, ρ is the density of formation water, and g is the acceleration due to gravity.

6. The method for determining the location of infiltrated wells based on the old well water ridge wave pattern according to claim 2, characterized in that, The step of substituting the water ridge height, the reservoir physical properties, the fluid parameters, and the production dynamic data into the horizontal well oil-water interface description formula includes: In the formula describing the oil-water interface in the horizontal well, Defined as the far-well region, the oil-water interface within the far-well region can be represented as: Where x is the x-coordinate of any point at the oil-water interface, and y is the y-coordinate of any point at the oil-water interface. G D is the water ridge height, and D is a control parameter. cone Where is the width of the water ridge, μ0 is the viscosity of the crude oil, q is the flow rate, and K is the flow rate. H ρ is the planar permeability, ρ is the density of the formation water, and g is the gravitational acceleration. In the formula describing the oil-water interface in the horizontal well, Defined as the near-wellbore region, the oil-water interface within the near-wellbore region can be represented as: Where x is the x-coordinate of any point at the oil-water interface, and y is the y-coordinate of any point at the oil-water interface. G Where is the water ridge height, D is the control parameter, μ0 is the crude oil viscosity, q is the flow rate data, and K is the flow rate. H ρ is the planar permeability, ρ is the density of the formation water, and g is the gravitational acceleration.

7. The method for determining the location of infiltrated wells based on the old well water ridge wave pattern according to claim 6, characterized in that, The step of making the formula dimensionless through normalization transformation includes: Normalize the formula to a dimensionless form. Where x′ is the dimensionless abscissa of any point at the oil-water interface, and y′ is the dimensionless ordinate of any point at the oil-water interface. G ′ is the dimensionless water ridge height, D′ is the dimensionless water ridge width, x is the x-coordinate of any point at the oil-water interface, y is the y-coordinate of any point at the oil-water interface, and K H ρ is the planar permeability, ρ is the density of formation water, g is the acceleration due to gravity, and μ is the groundwater density. o Here, q represents the viscosity of the crude oil, and q represents the flow rate. The relationship between the dimensionless x-coordinate x′ and the dimensionless y-coordinate y′ of any point on the oil-water interface in the far-well region is as follows: The relationship between the dimensionless x′ and the dimensionless y′ of any point on the oil-water interface in the near-wellbore region is as follows: Where x′ is the dimensionless abscissa of any point at the oil-water interface, and y′ is the dimensionless ordinate of any point at the oil-water interface. G ′ represents the dimensionless height of the water ridge, and D′ represents the dimensionless width of the water ridge.

8. The method for determining the location of infiltrated wells based on the old well water ridge wave pattern according to claim 7, characterized in that, The steps for calculating the correspondence between dimensionless cumulative oil production and dimensionless water ridge width and dimensionless water ridge height using the water displacement efficiency, based on the water displacement theory, include: The cumulative oil production is expressed as N using the water-drive oil recovery theory. p =VE d , where N P V represents the cumulative oil production of an old well at a certain moment, where V is the water drive volume and E is the total oil production. d For oil displacement efficiency; Its dimensionless form is N p ′=V′E d , Among them, N' p V' represents the cumulative dimensionless oil production of an old well at a certain moment, and E represents the dimensionless water drive volume. d For oil displacement efficiency, K H ρ is the planar permeability, ρ is the density of formation water, g is the acceleration due to gravity, and μ is the groundwater density. o Here, q represents the viscosity of the crude oil, and q represents the flow rate. By rationally dividing the horizontal cross-section of the water ridge in a bottom-water reservoir horizontal well, the cross-section can be considered as a combination of two semicircles and a rectangle. The final volume can be obtained by combining the integral of the dimensionless oil-water interface morphology equation and the volume formulas for solids of revolution and cylinders. Where L′ is the dimensionless well length of the old well, V′ is the dimensionless water drive volume, and y G ′ is the dimensionless water ridge height, D′ is the dimensionless water ridge width, x′ is the dimensionless abscissa of any point at the oil-water interface, dx is the derivative, and K H ρ is the planar permeability, g is the density of formation water, μ0 is the viscosity of crude oil, q is the flow rate, and L is the length of the old well. By substituting the formulas for calculating dimensionless water drive volume and oil displacement efficiency into the formula for calculating the cumulative dimensionless oil production at a certain moment in an old well, the relationship between the dimensionless cumulative oil production and the dimensionless water ridge width and dimensionless water ridge height can be obtained.

9. The method for determining the location of infiltrated wells based on the old well water ridge wavelet model according to claim 1, characterized in that, The step of determining the limit well spacing of the infill wells based on the data map comparison table includes: The maximum well spacing for new densification wells is set at half the width of the water ridge.

10. A device for determining the location of densely packed wells based on an old well water ridge wavelet model, characterized in that, include: The acquisition module is used to acquire reservoir physical property parameters, fluid parameters, and production dynamic data of the target horizontal well's production layer. The core relative permeability test module is used to obtain relative permeability curves and oil displacement efficiency through core relative permeability tests. The water ridge height acquisition module is used to determine the reservoir physical property parameters corresponding to side-drilling based on the relative permeability curve, and to calculate the water ridge height during side-drilling by combining the fluid parameters and the production dynamic data. The normalization module for the oil-water interface description formula is used to substitute the water ridge height, the reservoir physical property parameters, the fluid parameters, and the production dynamic data into the horizontal well oil-water interface description formula. The formula is made dimensionless through normalization deformation, and the correspondence between the dimensionless cumulative oil production and the dimensionless water ridge width and dimensionless water ridge height is calculated by combining the water drive theory and the oil displacement efficiency. The data mapping module is used to establish a data mapping comparison table of dimensionless cumulative oil production, dimensionless water ridge width, and dimensionless water ridge height, and to determine the limit well spacing of infill wells based on the data mapping comparison table. An encryption module is used to encrypt or regenerate branches of the target horizontal well based on the limit well spacing.

Images