Encrypted well location determination method and device based on old well water ridge waves and model

Abstract

The invention provides an infilled well location determination method and device based on old well water ridge waves and a model. The infilled well location determination method comprises the steps that reservoir physical property parameters and fluid parameters of a production horizon of a target horizontal well and production dynamic data of the target horizontal well are obtained; obtaining a relative permeability curve and oil displacement efficiency through a rock core relative permeability experiment; the height of a water ridge during sidetrack drilling is calculated; the water ridge height, the reservoir physical property parameters, the fluid parameters and the production dynamic data are substituted into a horizontal well oil-water interface description formula, and the formula is nondimensionalized through normalization deformation; the corresponding relation between the dimensionless cumulative oil production and the dimensionless water ridge width and the dimensionless water ridge height is calculated through the oil displacement efficiency according to the water-oil displacement theory; establishing a data map comparison table of the dimensionless cumulative oil production, the dimensionless water ridge width and the dimensionless water ridge height, and judging the limit well spacing of the infilled well according to the data map comparison table; and according to the limit well spacing, carrying out encryption or branch reconstruction on the target horizontal well.

Description

Technical Field

[0001] The invention relates to the technical field of oil and natural gas exploitation, and in particular to a method and a device for determining infill well locations based on an old well water ridge wave model. Background Art

[0002] At present, most marine sandstones in the eastern South China Sea are in the middle and late stages of development of producing oil fields and are in the "high-ultra-high" water-cut stage, making it difficult to continuously tap their potential. Most old wells are horizontal wells, and their large-scale use in the eastern South China Sea has ensured the rapid and efficient development of oil fields. The advantages of large-scale mobilization and efficient development of horizontal wells are particularly significant. However, the resulting difficulty is that it is difficult to predict the affected range of horizontal wells and the remaining oil around the wells. Especially for high-water-cut reservoirs, how to design infill wells or side-drilling branches of old wells to achieve the best development effect and economy is of great significance to the economic potential of remaining oil in offshore oil fields. In order to avoid rapid flooding of infill wells or side-drilling branches of old wells after they are put into production, it is crucial to accurately and quantitatively describe the morphology and range of water ridges in old wells during design.

[0003] However, existing methods for analyzing the range of water ridges in horizontal wells are mostly based on empirical formulas calibrated with numerical models and analysis of experimental data, which are often based on the understanding of reservoir-level laws. Given the heterogeneity of actual reservoirs and the large production differences between different wells, there are limitations in the accurate quantitative evaluation of water ridges in actual single wells, and there is still a lack of targeted inversion of water ridges in a specific old well. Summary of the invention

[0004] The purpose of the present invention is to provide a method and device for determining the location of encrypted wells based on the old well water ridge wave model, so as to solve the problem that the existing horizontal well water ridge wave and range analysis methods mentioned in the above background technology are mostly based on digital model calibration empirical formulas and experimental data analysis, which are often the understanding of reservoir-level laws and have limitations.

[0005] To achieve the above-mentioned purpose, the present invention provides the following technical solutions: a method for determining the infill well location based on the old well water ridge sweep model, the steps comprising: obtaining the reservoir physical property parameters and fluid parameters of the production layer of the target horizontal well, and the production dynamic data of the target horizontal well; obtaining the phase permeability curve and oil recovery efficiency through the core phase permeability experiment; determining the reservoir physical property parameters corresponding to the side drilling according to the phase permeability curve, and calculating the water ridge height during side drilling in combination with the fluid parameters and the production dynamic data; and calculating the water ridge height, the reservoir physical property parameters, the fluid parameters and the production dynamic data. The volume parameters and the production dynamic data are substituted into the horizontal well oil-water interface description formula, the formula is dimensionless through normalization and deformation, and the corresponding relationship between the dimensionless cumulative oil production and the dimensionless water ridge width and the dimensionless water ridge height is obtained by combining the water-to-oil theory with the oil recovery efficiency calculation; a data map comparison table of the dimensionless cumulative oil production and the dimensionless water ridge width and the 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 property parameters include permeability, porosity, water saturation, and volume coefficient; the fluid parameters include fluid density and viscosity; and the production dynamic data include daily oil production, daily liquid production, water content, bottom hole flow pressure, and initial oil column height.

[0007] Optionally, the step of obtaining the phase permeability curve and oil recovery efficiency through the core phase permeability test includes: the oil recovery efficiency is calculated by the initial water saturation and residual oil saturation obtained through the core test, and the calculation formula is: Among them, E d is the oil displacement efficiency, S wi is the initial water saturation, S or is the residual oil saturation.

[0008] Optionally, the step of obtaining the phase permeability curve and oil recovery efficiency through the core phase permeability experiment further includes: using the approximate theoretical water drive curve to correct the oil recovery efficiency, and the correction formula is: in is the average water saturation, S wi is the initial water saturation, S or is the residual oil saturation, w and b are parameters calculated from the approximate theoretical water drive curve.

[0009] Optionally, the step of determining the reservoir physical property parameters corresponding to the sidetracking according to the phase permeability curve, and calculating the water ridge height during the sidetracking in combination with the fluid parameters and the production dynamic data includes: determining the reservoir physical property parameters corresponding to the sidetracking according to the phase permeability curve, and calculating the dimensionless fluid production index during the sidetracking, the calculation formula is: Among them JDL is the dimensionless liquid extraction index, K ro is the oil phase permeability, K rw is the water phase permeability, μ o is the crude oil viscosity, μ w is the viscosity of water, B o is the crude oil volume coefficient, B w is the water volume coefficient; the liquid production index during side drilling is calculated based on the production dynamic data combined with the dimensionless liquid production index, and the calculation formula is ΔP initial =P res-initial -bhp initial , J=J initial ×J DL , where ΔP initial is the initial reservoir pressure difference, P res-initial is the initial reservoir pressure, bhp initial is the bottom hole flowing pressure at the initial stage of production, J initial is the production index at the initial stage of production, q is the flow rate data, J is the production index during side drilling, and J DL is a dimensionless liquid production index; the decrease in reservoir pressure compared to the time of production is calculated based on the production dynamic data combined with the liquid production index, and the calculation formula is: P res = bhp + ΔP, where ΔP is the drop in reservoir pressure compared to the time of production, q is the flow rate data, J is the liquid production index during side drilling, and P res is the reservoir pressure during sidetracking, bhp is the bottom hole flow pressure; the water ridge height during sidetracking is calculated based on the production dynamic data, the fluid parameters, and the drop in the reservoir pressure compared to the time of production. The calculation formula is: Among them, y G is the water ridge height during sidetracking, y G-initial is the initial oil column height, P res-initial is the initial reservoir pressure, P res is the reservoir pressure during sidetracking, ρ is the density of formation water, and g is the gravitational acceleration.

[0010] Optionally, the step of substituting 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 includes: in the horizontal well oil-water interface description formula, is defined as the far-well region, then the oil-water interface in the far-well region can be expressed as Among them, x is the horizontal coordinate of any point on the oil-water interface, y is the vertical coordinate of any point on the oil-water interface, and y G is the water ridge height, D is the control parameter, D coneis the water ridge width, μ0 is the crude oil viscosity, q is the flow data, K H is the plane permeability, ρ is the density of formation water, and g is the gravitational acceleration; in the horizontal well oil-water interface description formula, is defined as the near-wellbore area, then the oil-water interface in the near-wellbore area can be expressed as Among them, x is the horizontal coordinate of any point on the oil-water interface, y is the vertical coordinate of any point on the oil-water interface, and y G is the water ridge height, D is the control parameter, μ0 is the crude oil viscosity, q is the flow data, K H is the plane permeability, ρ is the density of formation water, and g is the gravitational acceleration.

[0011] Optionally, the step of transforming the formula into a dimensionless form by normalization includes: normalizing and transforming the formula into a dimensionless form, Among them, x′ is the dimensionless horizontal coordinate of any point on the oil-water interface, y′ is the dimensionless vertical coordinate of any point on the oil-water interface, and y G ′ is the dimensionless water ridge height, D′ is the dimensionless water ridge width, x is the horizontal coordinate of any point on the oil-water interface, y is the vertical coordinate of any point on the oil-water interface, K H is the plane permeability, ρ is the density of formation water, g is the gravitational acceleration, μ o is the viscosity of crude oil, q is the flow rate data; the relationship between the dimensionless horizontal coordinate x′ and the dimensionless vertical coordinate y′ of any point on the oil-water interface in the far well area is: The relationship between the dimensionless horizontal coordinate x′ and the dimensionless vertical coordinate y′ of any point on the oil-water interface near the wellbore is: Among them, x′ is the dimensionless horizontal coordinate of any point on the oil-water interface, y′ is the dimensionless vertical coordinate of any point on the oil-water interface, and y G ′ is the dimensionless water ridge height, and D′ is the dimensionless water ridge width.

[0012] Optionally, the step of calculating the corresponding relationship between the dimensionless cumulative oil production and the dimensionless water ridge width and the dimensionless water ridge height by combining the water-to-oil theory with the oil displacement efficiency comprises: expressing the cumulative oil production as N by the water-to-oil theory p =VE d , where N P is the cumulative oil production of the old well at a certain moment, V is the water drive volume, E d is the oil displacement efficiency; its dimensionless form is N p ′=V′E d , Among them, N' p is the cumulative dimensionless oil production of the old well at a certain moment, V' is the dimensionless water drive volume, E dis the oil displacement efficiency, K H is the plane permeability, ρ is the density of formation water, g is the gravitational acceleration, μ o is the viscosity of crude oil, q is the flow rate data; the horizontal section of the water ridge of the horizontal well in the bottom water reservoir is reasonably split, and the section can be regarded 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 formula of the rotating body and cylinder

[0013] Where L' is the dimensionless well length of the old well, V' is the dimensionless water drive volume, y G ′ is the dimensionless water ridge height, D′ is the dimensionless water ridge width, x′ is the dimensionless horizontal coordinate of any point on the oil-water interface, dx is the derivative, K H is the plane permeability, ρ is the density of formation water, g is the acceleration of gravity, μ0 is the viscosity of crude oil, q is the flow data, and L is the length of the old well; the dimensionless water drive volume and oil recovery efficiency calculation formula is substituted into the cumulative dimensionless oil production calculation formula of the old well at a certain moment, and the relationship between the dimensionless cumulative oil production and the dimensionless water ridge width and dimensionless water ridge height can be obtained.

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

[0015] On the other hand, the present invention also provides an infill well location determination device based on the old well water ridge sweep model, comprising: an acquisition module, used to obtain reservoir physical property parameters, fluid parameters, and production dynamic data of the target horizontal well's production layer; a core phase permeability experiment module, used to obtain a phase permeability curve and oil recovery efficiency through a core phase permeability experiment; a water ridge height acquisition module, used to determine the reservoir physical property parameters corresponding to side drilling according to the phase permeability curve, and calculate the water ridge height during side drilling in combination with the fluid parameters and the production dynamic data; an oil-water interface description formula normalization module, used to convert 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 dimensionless through normalization and deformation, and the corresponding relationship between the dimensionless cumulative oil production and the dimensionless water ridge width and the dimensionless water ridge height is obtained through the oil recovery efficiency calculation in combination with the water-to-oil theory; a data map module is used to establish a data map comparison table of the dimensionless cumulative oil production and the dimensionless water ridge width and the dimensionless water ridge height, and determine the limit well spacing of the infill well according to the data map comparison table; an infill module is used to infill or re-branch the target horizontal well according to the limit well spacing.

[0016] Compared with the prior art, the present invention has the following beneficial effects:

[0017] Through the acquired parameter information and dynamic data, the phase permeability curve and oil recovery efficiency are obtained through the core phase permeability experiment, and then the water ridge height during side drilling is calculated. Based on the horizontal well oil-water interface description formula, the relationship between the dimensionless cumulative oil production and the two variables of dimensionless water ridge height and dimensionless water ridge width is derived through normalized deformation combined with water-drive theory. Finally, the established data map is used to determine the extreme well spacing of the infill wells, and the target horizontal wells are infilled or re-branched according to the extreme well spacing. This overcomes the limitations of the conventional mainstream water ridge analysis method, such as the narrow application range and the inability to give accurate quantitative evaluation results for specific production wells. It is a more reasonable and effective old well sweep range analysis and re-branch design evaluation method. BRIEF DESCRIPTION OF THE DRAWINGS

[0018] Figure 1 The figure is a schematic flow chart of the steps of the method of the present invention.

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

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

[0021] Figure 4 It is a schematic diagram of splitting the bottom surface of the water ridge of the present invention.

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

[0023] Figure 6 It is the calculation result of the water ridge width of the old well and the schematic diagram of the optimization of the well location of the new well in the example of the present invention. DETAILED DESCRIPTION

[0024] The scheme of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, rather than all the embodiments.

[0025] It should be noted that the terms "first", "second", etc. in the specification and claims of the present application and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequential order. It should be understood that the data used in this way can be interchanged where appropriate, so that the embodiments of the present application described here. In addition, the terms "including" and "having" and any of their variations are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device comprising 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.

[0026] It will be understood by those skilled in the art that, unless expressly stated, the singular forms "one", "said", and "the" used herein may also include plural forms. It should be further understood that the term "comprising" used in the specification of the present application refers to 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 refer to an element as being "connected" or "coupled" to another element, it may be directly connected or coupled to the other element, or there may be an intermediate element. In addition, the "connection" or "coupling" used herein may include wireless connection or wireless coupling. The term "and / or" used herein includes all or any unit 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 (including technical and scientific terms) used herein have the same meaning as those generally understood by those skilled in the art to which this application belongs. It should also be understood that terms such as those defined in general dictionaries should be understood to have meanings consistent with those in the context of the prior art, and will not be interpreted with idealized or overly formal meanings unless specifically defined as here.

[0028] It should be understood that the sequence numbers and sizes of the steps in this embodiment do not mean 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 the embodiment of the present application.

[0029] It should be noted that, in the absence of conflict, the embodiments and features in the embodiments of the present application can be combined with each other. The present application will be described in detail below with reference to the accompanying drawings and in combination with the embodiments.

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

[0031] S100, obtaining reservoir physical property parameters and fluid parameters of a production layer of a target horizontal well, as well as production dynamic data of the target horizontal well.

[0032] Specifically, the reservoir physical parameters include oil layer thickness, water avoidance height, that is, 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; the production dynamic data includes the development date and the corresponding daily liquid production, daily water production, and water content, which are derived from the oil well production dynamic data acquisition system.

[0033] S200, obtain the phase permeability curve and oil recovery efficiency through core phase permeability experiment.

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

[0035] S300, determining the reservoir physical property parameters corresponding to the sidetracking according to the phase permeability curve, and calculating the water ridge height during the sidetracking in combination with the fluid parameters and the production dynamic data.

[0036] Specifically, taking the marine sandstone horizontal well in the eastern South China Sea as an example, the reservoir physical properties such as permeability, porosity, water saturation, fluid parameters such as fluid density and viscosity of the production layer of the target well, as well as the production dynamic data such as daily oil production, daily liquid production, water content, bottom hole pressure and so on of the old well are obtained; the phase permeability curve is obtained through the core phase permeability experiment, and the dimensionless liquid production index under the corresponding water content during MRC side drilling is calculated; the liquid production index in the initial stage of production is calculated by the initial reservoir pressure and the bottom hole pressure and flow rate data of the old well in the early stage of production, and the liquid production index during side drilling is calculated in combination with the dimensionless liquid production index; the reservoir pressure during side drilling is calculated according to the flow rate, bottom hole pressure and liquid production index of the old well before side drilling, and then the decrease in reservoir pressure compared with the time of production is calculated; the water ridge height during MRC side drilling is calculated according to the initial oil column height of the old well before side drilling, the decrease in reservoir pressure, the density of formation water, and the gravity acceleration.

[0037] S400, substituting 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, making the formula dimensionless through normalization and deformation, and combining the water drive oil theory with the oil drive efficiency to obtain the corresponding relationship between the dimensionless cumulative oil production and the dimensionless water ridge width and the dimensionless water ridge height.

[0038] Specifically, combined with the cumulative oil production of the old well during MRC sidetracking, the dimensionless cumulative oil production and the dimensionless oil column height are calculated, where the oil recovery efficiency can be calculated by obtaining the initial water saturation and residual oil saturation through core tests, and then substituted into the above calculation formula to obtain the relationship between the dimensionless cumulative oil production and the dimensionless water ridge width and dimensionless water ridge height.

[0039] S500, establishing a data graph comparison table of dimensionless cumulative oil production, dimensionless water ridge width, and dimensionless water ridge height, and determining the limit well spacing of infill wells according to the data graph comparison table.

[0040] S600, infilling or re-branching the target horizontal well according to the limit well spacing.

[0041] Specifically, taking a test well 20H2 in a certain reservoir as an example, the initial water saturation is 0.2, the residual oil saturation is 0.26, and the oil recovery efficiency of the reservoir is calculated to be 0.56; according to the dynamic data, the water ridge height of the 20H2 well during MRC side drilling is 10 meters, and the cumulative oil production is 146,000 cubic meters. The corresponding column of the initial oil column height is found in the comparison table, and the dimensionless water ridge width is determined according to the dimensionless cumulative oil production. Finally, the dimensionless water ridge width is converted into the water ridge width, and the water ridge width is calculated to be about 125 meters; the 1 / 2 range of the water ridge width is the area greatly affected by the old well water ridge, and the new branch is 1 / 2 of the old branch water ridge. / 2 width is less than the total length * 0.25 as the limit, and the length of the reconstructed branch in the large affected area is calculated according to the length of the reconstructed branch and the side drilling angle, and this is used as the basis for evaluating whether the design of the reconstructed branch is reasonable; in this case, after the final calculation, the side drilling branch lengths of the MRC design of the 20H2 well are 548 meters and 743 meters respectively, and the corresponding side drilling branch angle is greater than 15° to obtain a better side drilling effect; if the well network needs to be intensified, the reasonable well spacing range should be greater than 1 / 2 width of the old branch water ridge, that is, 62.5 meters, so as to achieve the optimization evaluation of the well location of the intensified well or MRC reconstructed branch.

[0042] It can be understood that, through the obtained parameter information and dynamic data, the phase permeability curve and oil displacement efficiency are obtained through the core phase permeability experiment, and then the water ridge height during side drilling is calculated. Based on the horizontal well oil-water interface description formula, the relationship between the dimensionless cumulative oil production and the two variables of dimensionless water ridge height and dimensionless water ridge width is derived through normalized deformation combined with water displacement theory. Finally, the established data map is used to determine the extreme well spacing of the infill well, and the target horizontal well is infilled or re-branched according to the extreme well spacing. The limitations of the conventional mainstream water ridge analysis method, such as the narrow application range and the inability to give accurate quantitative evaluation results for specific production wells, are overcome. This is a more reasonable and effective old well sweep range analysis and re-branch design evaluation method. This method can realize the accurate evaluation of the old well water ridge sweep range in the inter-well infill and old well side drilling MRC design, and at the same time guide the optimization of the infill well location, and continuously tap the potential of old oil fields in the high water content stage, which has a good application prospect.

[0043] Optionally, the reservoir physical property parameters include permeability, porosity, water saturation, and volume coefficient; the fluid parameters include fluid density and viscosity; and the production dynamic data include daily oil production, daily liquid production, water content, bottom hole flow pressure, and initial oil column height.

[0044] Optionally, the step of obtaining the phase permeability curve and oil recovery efficiency through the core phase permeability test includes: the oil recovery efficiency is calculated by the initial water saturation and residual oil saturation obtained through the core test, and the calculation formula is: Among them, E d is the oil displacement efficiency, S wi is the initial water saturation, S or is the residual oil saturation.

[0045] Optionally, the step of obtaining the phase permeability curve and oil recovery efficiency through the core phase permeability experiment further includes: using the approximate theoretical water drive curve to correct the oil recovery efficiency, and the correction formula is: in is the average water saturation, S wi is the initial water saturation, S or is the residual oil saturation, w and b are parameters calculated from the approximate theoretical water drive curve.

[0046] Please refer to Figure 3 Optionally, the step of determining the reservoir physical property parameters corresponding to the sidetracking according to the phase permeability curve, and calculating the water ridge height during the sidetracking in combination with the fluid parameters and the production dynamic data includes: determining the reservoir physical property parameters corresponding to the sidetracking according to the phase permeability curve, and calculating the dimensionless fluid production index during the sidetracking, the calculation formula is: Among them J DL is the dimensionless liquid extraction index, K ro is the oil phase permeability, K rw is the water phase permeability, μ o is the crude oil viscosity, μ w is the viscosity of water, B o is the crude oil volume coefficient, B w is the water volume coefficient; the liquid production index during side drilling is calculated based on the production dynamic data combined with the dimensionless liquid production index, and the calculation formula is ΔP initial =P res-initial -bhp initial , J=J initial ×J DL , where ΔP initial is the initial reservoir pressure difference, P res-initial is the initial reservoir pressure, bhp initial is the bottom hole flowing pressure at the initial stage of production, J initial is the production index at the initial stage of production, q is the flow rate data, J is the production index during side drilling, and J DL is a dimensionless liquid production index; the decrease in reservoir pressure compared to the time of production is calculated based on the production dynamic data combined with the liquid production index, and the calculation formula is: P res = bhp + ΔP, where ΔP is the drop in reservoir pressure compared to the time of production, q is the flow rate data, J is the liquid production index during side drilling, and P resis the reservoir pressure during sidetracking, bhp is the bottom hole flow pressure; the water ridge height during sidetracking is calculated based on the production dynamic data, the fluid parameters, and the drop in the reservoir pressure compared to the time of production. The calculation formula is: Among them, y G is the water ridge height during sidetracking, y G-initial is the initial oil column height, P res-initial is the initial reservoir pressure, P res is the reservoir pressure during sidetracking, ρ is the density of formation water, and g is the gravitational acceleration.

[0047] Optionally, the step of substituting 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 includes: in the horizontal well oil-water interface description formula, is defined as the far-well region, then the oil-water interface in the far-well region can be expressed as Among them, x is the horizontal coordinate of any point on the oil-water interface, y is the vertical coordinate of any point on the oil-water interface, and y G is the water ridge height, D is the control parameter, D cone is the water ridge width, μ0 is the crude oil viscosity, q is the flow data, K H is the plane permeability, ρ is the density of formation water, and g is the gravitational acceleration; in the horizontal well oil-water interface description formula, is defined as the near-wellbore area, then the oil-water interface in the near-wellbore area can be expressed as Among them, x is the horizontal coordinate of any point on the oil-water interface, y is the vertical coordinate of any point on the oil-water interface, and y G is the water ridge height, D is the control parameter, μ0 is the crude oil viscosity, q is the flow data, K H is the plane permeability, ρ is the density of formation water, and g is the gravitational acceleration.

[0048] Optionally, the step of transforming the formula into a dimensionless form by normalization includes: normalizing and transforming the formula into a dimensionless form, Among them, x′ is the dimensionless horizontal coordinate of any point on the oil-water interface, y′ is the dimensionless vertical coordinate of any point on the oil-water interface, and y G ′ is the dimensionless water ridge height, D′ is the dimensionless water ridge width, x is the horizontal coordinate of any point on the oil-water interface, y is the vertical coordinate of any point on the oil-water interface, K H is the plane permeability, ρ is the density of formation water, g is the gravitational acceleration, μ o is the viscosity of crude oil, q is the flow rate data; the relationship between the dimensionless horizontal coordinate x′ and the dimensionless vertical coordinate y′ of any point on the oil-water interface in the far well area is: The relationship between the dimensionless horizontal coordinate x′ and the dimensionless vertical coordinate y′ of any point on the oil-water interface near the wellbore is: Among them, x′ is the dimensionless horizontal coordinate of any point on the oil-water interface, y′ is the dimensionless vertical coordinate of any point on the oil-water interface, and y G ′ is the dimensionless water ridge height, and D′ is the dimensionless water ridge width.

[0049] Please refer to Figure 4 and Figure 5 Optionally, the step of calculating the corresponding relationship between the dimensionless cumulative oil production and the dimensionless water ridge width and the dimensionless water ridge height by combining the water-to-oil theory with the oil displacement efficiency includes: expressing the cumulative oil production as N by the water-to-oil theory p =VE d , where N P is the cumulative oil production of the old well at a certain moment, V is the water drive volume, E d is the oil displacement efficiency; its dimensionless form is N p ′=V′E d , Among them, N' p is the cumulative dimensionless oil production of the old well at a certain moment, V' is the dimensionless water drive volume, E d is the oil displacement efficiency, K H is the plane permeability, ρ is the density of formation water, g is the gravitational acceleration, μ o is the viscosity of crude oil, q is the flow rate data; the horizontal section of the water ridge of the horizontal well in the bottom water reservoir is reasonably split, and the section can be regarded 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 formula of the rotating body and cylinder Where L' is the dimensionless well length of the old well, V' is the dimensionless water drive volume, y G ′ is the dimensionless water ridge height, D′ is the dimensionless water ridge width, x′ is the dimensionless horizontal coordinate of any point on the oil-water interface, dx is the derivative, K H is the plane permeability, ρ is the density of formation water, g is the acceleration of gravity, μ0 is the viscosity of crude oil, q is the flow data, and L is the length of the old well; the dimensionless water drive volume and oil recovery efficiency calculation formula is substituted into the cumulative dimensionless oil production calculation formula of the old well at a certain moment, and the relationship between the dimensionless cumulative oil production and the dimensionless water ridge width and dimensionless water ridge height can be obtained.

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

[0051] Specifically, a data table comparing dimensionless cumulative oil production, dimensionless water ridge width, and dimensionless water ridge height is established as shown in Table 1 below:

[0052] Table 1: Data map comparison table.

[0053]

[0054]

[0055] On the other hand, the present invention also provides an infill well location determination device based on the old well water ridge sweep model, comprising: an acquisition module, used to obtain reservoir physical property parameters, fluid parameters, and production dynamic data of the target horizontal well's production layer; a core phase permeability experiment module, used to obtain a phase permeability curve and oil recovery efficiency through a core phase permeability experiment; a water ridge height acquisition module, used to determine the reservoir physical property parameters corresponding to side drilling according to the phase permeability curve, and calculate the water ridge height during side drilling in combination with the fluid parameters and the production dynamic data; an oil-water interface description formula normalization module, used to convert 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 dimensionless through normalization and deformation, and the corresponding relationship between the dimensionless cumulative oil production and the dimensionless water ridge width and the dimensionless water ridge height is obtained through the oil recovery efficiency calculation in combination with the water-to-oil theory; a data map module is used to establish a data map comparison table of the dimensionless cumulative oil production and the dimensionless water ridge width and the dimensionless water ridge height, and determine the limit well spacing of the infill well according to the data map comparison table; an infill module is used to infill or re-branch the target horizontal well according to the limit well spacing.

[0056] On the other hand, the present invention further provides a computer device, including a memory and a processor, wherein the memory stores a computer program, and the processor implements the steps of the above control method when executing the computer program.

[0057] On the other hand, the present invention further provides a computer-readable storage medium having a computer program stored thereon, and the steps of the above control method are implemented when the computer program is executed by a processor.

[0058] If the integrated unit is implemented in the form of 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, which is stored in a storage medium and includes several instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the methods of each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM), random access memory (RAM), disk or optical disk, etc., various media that can store program codes.

[0059] A person skilled in the art can understand that all or part of the processes in the above-mentioned embodiments can be completed by instructing related hardware through a computer program. 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-mentioned methods; wherein, any reference 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; as an illustration and not limitation, RAM is available in various forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).

[0060] The above are only embodiments of the present invention, and are not intended to limit the patent scope of the present invention. Any equivalent transformations made using the contents of the specification and drawings of the present invention, or directly or indirectly applied in related technical fields, are also included in the patent protection scope of the present invention.

Claims

1. A method for determining the location of infill wells based on the old well water ridge wave model, characterized in that the steps include: Acquiring reservoir physical property parameters and fluid parameters of the production layer of the target horizontal well, as well as production performance data of the target horizontal well; Obtain phase permeability curve and oil displacement efficiency through core phase permeability test; Determine the reservoir physical property parameters corresponding to the sidetracking according to the phase permeability curve, and calculate the water ridge height during the sidetracking in combination with the fluid parameters and the production dynamic data; Substituting 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, transforming the formula into dimensionless through normalization deformation, and combining the water drive oil theory with the oil drive efficiency to obtain the corresponding relationship between the dimensionless cumulative oil production and the dimensionless water ridge width and the dimensionless water ridge height; Establishing a data map comparison table of the dimensionless cumulative oil production, the dimensionless water ridge width, and the dimensionless water ridge height, and determining the limit well spacing of the infill wells according to the data map comparison table; The target horizontal well is infilled or branched according to the limit well spacing.

2. The method for determining the infilled well locations based on the old well water ridge and wave model according to claim 1 is characterized in that: The reservoir physical property parameters include permeability, porosity, water saturation, and volume coefficient; the fluid parameters include fluid density and viscosity; and the production dynamic data include daily oil production, daily liquid production, water content, bottom hole flow pressure, and initial oil column height.

3. The method for determining the infilled well locations based on the old well water ridge wave model according to claim 2 is characterized in that: The step of obtaining the phase permeability curve and oil displacement efficiency through the core phase permeability experiment comprises: The oil displacement efficiency is calculated by the initial water saturation and residual oil saturation obtained from the core test, and the calculation formula is: Among them, E d is the oil displacement efficiency, S wi is the initial water saturation, S or is the residual oil saturation.

4. The method for determining the infilled well locations based on the old well water ridge wave model according to claim 3 is characterized in that: The step of obtaining the phase permeability curve and oil displacement efficiency through the core phase permeability experiment also includes: The oil displacement efficiency is corrected using the approximate theoretical water displacement curve. The correction formula is: in is the average water saturation, S wi is the initial water saturation, S or is the residual oil saturation, w and b are parameters calculated from the approximate theoretical water drive curve.

5. The method for determining the infilled well locations based on the old well water ridge and wave model according to claim 2 is characterized in that: The step of determining the reservoir physical property parameters corresponding to the sidetracking according to the phase permeability curve, and calculating the water ridge height during the sidetracking in combination with the fluid parameters and the production dynamic data comprises: The reservoir physical property parameters corresponding to the sidetracking are determined according to the phase permeability curve, and the dimensionless liquid production index during sidetracking is calculated. The calculation formula is: Among them J DL is the dimensionless liquid extraction index, K ro is the oil phase permeability, K rw is the water phase permeability, μ o is the crude oil viscosity, μ w is the viscosity of water, B o is the crude oil volume coefficient, B w is the water volume coefficient; The production index during sidetracking is calculated based on the production dynamic data combined with the dimensionless production index. The calculation formula is ΔP initial =P res-initial -bhp initial , J=J initial ×J DL , where ΔP initial is the initial reservoir pressure difference, P res-initial is the initial reservoir pressure, bhp initial is the bottom hole flowing pressure at the initial stage of production, J initial is the production index at the initial stage of production, q is the flow rate data, J is the production index during side drilling, and J DL is the dimensionless liquid extraction index; The decrease in reservoir pressure compared to the time of production is calculated based on the production dynamic data combined with the liquid production index. The calculation formula is: P res = bhp + ΔP, where ΔP is the drop in reservoir pressure compared to the time of production, q is the flow rate data, J is the liquid production index during side drilling, and P res is the reservoir pressure during sidetracking, bhp is the bottom hole flowing pressure; The water ridge height during sidetracking is calculated based on the production dynamic data, the fluid parameters, and the drop in reservoir pressure compared to the time of production. The calculation formula is: Among them, y G is the water ridge height during sidetracking, y G-initial is the initial oil column height, P res-initial is the initial reservoir pressure, P res is the reservoir pressure during sidetracking, ρ is the density of formation water, and g is the gravitational acceleration.

6. The method for determining the infilled well locations based on the old well water ridge and wave model according to claim 2 is characterized in that: The step of substituting 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 includes: In the horizontal well oil-water interface description formula, is defined as the far-well region, then the oil-water interface in the far-well region can be expressed as Among them, x is the horizontal coordinate of any point on the oil-water interface, y is the vertical coordinate of any point on the oil-water interface, and y G is the water ridge height, D is the control parameter, D cone is the water ridge width, μ0 is the crude oil viscosity, q is the flow data, K H is the plane permeability, ρ is the density of formation water, and g is the gravitational acceleration; In the horizontal well oil-water interface description formula, is defined as the near-wellbore area, then the oil-water interface in the near-wellbore area can be expressed as Among them, x is the horizontal coordinate of any point on the oil-water interface, y is the vertical coordinate of any point on the oil-water interface, and y G is the water ridge height, D is the control parameter, μ0 is the crude oil viscosity, q is the flow data, K H is the plane permeability, ρ is the density of formation water, and g is the gravitational acceleration.

7. The method for determining the infilled well locations based on the old well water ridge and wave model according to claim 6 is characterized in that: The step of transforming the formula dimensionlessly by normalization includes: Normalize the formula to a dimensionless form, Among them, x′ is the dimensionless horizontal coordinate of any point on the oil-water interface, y′ is the dimensionless vertical coordinate of any point on the oil-water interface, and y G ′ is the dimensionless water ridge height, D′ is the dimensionless water ridge width, x is the horizontal coordinate of any point on the oil-water interface, y is the vertical coordinate of any point on the oil-water interface, K H is the plane permeability, ρ is the density of formation water, g is the gravitational acceleration, μ o is the viscosity of crude oil, q is the flow rate data; The relationship between the dimensionless horizontal coordinate x′ and the dimensionless vertical coordinate y′ of any point on the oil-water interface in the far-well area is: The relationship between the dimensionless horizontal coordinate x′ and the dimensionless vertical coordinate y′ of any point on the oil-water interface near the wellbore is: Among them, x′ is the dimensionless horizontal coordinate of any point on the oil-water interface, y′ is the dimensionless vertical coordinate of any point on the oil-water interface, and y G ′ is the dimensionless water ridge height, and D′ is the dimensionless water ridge width.

8. The method for determining the infilled well locations based on the old well water ridge and wave model according to claim 7 is characterized in that: The step of calculating the corresponding relationship between the dimensionless cumulative oil production and the dimensionless water ridge width and the dimensionless water ridge height by combining the water flooding theory with the oil flooding efficiency comprises: The cumulative oil production is expressed as N by water displacement theory. p =VE d , where N P is the cumulative oil production of the old well at a certain moment, V is the water drive volume, E d for oil recovery efficiency; Its dimensionless form is N p ′=V′E d , Among them, N' p is the cumulative dimensionless oil production of the old well at a certain moment, V' is the dimensionless water drive volume, E d is the oil displacement efficiency, K H is the plane permeability, ρ is the density of formation water, g is the gravitational acceleration, μ o is the viscosity of crude oil, q is the flow rate data; The horizontal section of the water ridge of the horizontal well in the bottom water reservoir is reasonably split, and the section can be regarded 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 formula of the rotation body and cylinder. Where L' is the dimensionless well length of the old well, V' is the dimensionless water drive volume, y G ′ is the dimensionless water ridge height, D′ is the dimensionless water ridge width, x′ is the dimensionless horizontal coordinate of any point on the oil-water interface, dx is the derivative, K H is the plane permeability, ρ is the density of formation water, g is the acceleration of gravity, μ0 is the viscosity of crude oil, q is the flow data, and L is the length of the old well; Substituting the dimensionless water drive volume and oil recovery efficiency calculation formula into the cumulative dimensionless oil production calculation formula of the old well at a certain moment, 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 infilled well locations based on the old well water ridge and wave model according to claim 1 is characterized in that: The step of determining the limit well spacing of the infill wells according to the data spectrum comparison table comprises: The maximum well spacing for adding new wells is 1 / 2 of the water ridge width.

10. A device for determining the location of encrypted wells based on the old well water ridge wave model, characterized in that: include: An acquisition module, used to acquire reservoir physical property parameters and fluid parameters of the production layer of the target horizontal well, as well as production dynamic data of the target horizontal well; Core phase permeability experiment module, used to obtain phase permeability curve and oil recovery efficiency through core phase permeability experiment; A water ridge height acquisition module, used to determine the reservoir physical property parameters corresponding to the sidetracking according to the phase permeability curve, and calculate the water ridge height during the sidetracking in combination with the fluid parameters and the production dynamic data; The oil-water interface description formula normalization module 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, transform the formula into dimensionless through normalization deformation, and calculate the corresponding relationship between the dimensionless cumulative oil production and the dimensionless water ridge width and the dimensionless water ridge height through the oil recovery efficiency in combination with the water recovery theory; A data map module is used to 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 according to the data map comparison table; The encryption module is used to infill or re-branch the target horizontal well according to the limit well spacing.

Images