CN108491361B - Heterogeneous reservoir permeability distribution inversion method and device - Google Patents

Heterogeneous reservoir permeability distribution inversion method and device Download PDF

Info

Publication number
CN108491361B
CN108491361B CN201810187730.5A CN201810187730A CN108491361B CN 108491361 B CN108491361 B CN 108491361B CN 201810187730 A CN201810187730 A CN 201810187730A CN 108491361 B CN108491361 B CN 108491361B
Authority
CN
China
Prior art keywords
permeability
determining
inversion
layer
reservoir
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810187730.5A
Other languages
Chinese (zh)
Other versions
CN108491361A (en
Inventor
王敬
刘慧卿
刘化普
冯其红
钱其豪
单高军
由春梅
郭军辉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum Beijing
China University of Petroleum East China
Original Assignee
China University of Petroleum Beijing
China University of Petroleum East China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Petroleum Beijing, China University of Petroleum East China filed Critical China University of Petroleum Beijing
Priority to CN201810187730.5A priority Critical patent/CN108491361B/en
Publication of CN108491361A publication Critical patent/CN108491361A/en
Application granted granted Critical
Publication of CN108491361B publication Critical patent/CN108491361B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Abstract

The invention discloses a heterogeneous reservoir permeability distribution inversion method and a heterogeneous reservoir permeability distribution inversion device, wherein the method comprises the following steps: determining Lorentz curve inversion expressions corresponding to different permeability variation coefficients of the heterogeneous oil reservoir; determining the abscissa of a coordinate point corresponding to each layer of oil deposit on the Lorentz curve according to the number of oil deposit layers and the thickness of each layer of oil deposit; determining the vertical coordinates corresponding to different horizontal coordinates by using the Lorentz curve inversion expression corresponding to different permeability variation coefficients; and determining the permeability of the corresponding coordinate points of each layer of oil deposit according to the vertical coordinates corresponding to different horizontal coordinates. The invention solves the problems existing in the current heterogeneous oil reservoir permeability distribution inversion, establishes a set of simple and easy-to-use heterogeneous oil reservoir permeability distribution inversion solution with uniform standard, improves the understanding of the influence of the heterogeneity of the heterogeneous oil reservoir on the oil reservoir development, and can better guide the heterogeneous oil reservoir development.

Description

Heterogeneous reservoir permeability distribution inversion method and device
Technical Field
The invention relates to the technical field of heterogeneous reservoir development in the petroleum industry, in particular to a heterogeneous reservoir permeability distribution inversion method and device.
Background
Interlayer heterogeneity commonly exists in the oil reservoir, and commonly used parameters for measuring the interlayer heterogeneity comprise level difference, breakthrough coefficient, variation coefficient and the like, wherein the variation coefficient is most commonly used, reflects the degree of deviation of single-layer permeability from the overall average value, is the most important parameter for evaluating the macroscopic heterogeneity of the oil reservoir, and has a permeability variation coefficient VKLarger indicates more heterogeneous reservoir. The methods commonly used at present for calculating the coefficient of variation are the statistical method, the Dykstra H method and the lorentz method. The statistical method is suitable for reservoirs with any permeability distribution type, but the coefficient of variation value calculated by the method is between 0 and infinity, and the measurement of the heterogeneous degree lacks relativity; the Dykstra H method is an empirical method, the coefficient of variation value calculated by the method is between 0 and 1, and the method is only suitable for reservoirs with lognormal distribution of permeability; the Lorentz method has good applicability, is suitable for reservoirs with any permeability distribution type, and has a calculated coefficient of variation value of 0-1.
However, although the lorentz method can calculate the variation coefficient of a certain reservoir and analyze the heterogeneity, it is sometimes necessary to construct a heterogeneous reservoir with a certain or several specific variation coefficients to study the influence of the heterogeneity on the reservoir development effect. But now a more efficient inversion method to obtain the corresponding permeability distribution is lacking. At present, the heterogeneity research of a certain oil reservoir or oil reservoirs with certain specific variation coefficients is obtained by a trial calculation method, and the permeability distribution obtained by the method lacks a uniform standard.
Disclosure of Invention
The embodiment of the invention provides an inversion method of permeability distribution of a heterogeneous oil reservoir, which is used for establishing a set of simple and easy-to-use inversion solution of permeability distribution of the heterogeneous oil reservoir with unified standard, and comprises the following steps:
determining Lorentz curve inversion expressions corresponding to different permeability variation coefficients of the heterogeneous oil reservoir;
determining the abscissa of a coordinate point corresponding to each layer of oil deposit on the Lorentz curve according to the number of oil deposit layers and the thickness of each layer of oil deposit;
determining the vertical coordinates corresponding to different horizontal coordinates by using the Lorentz curve inversion expression corresponding to different permeability variation coefficients;
determining the permeability of the corresponding coordinate points of each layer of oil deposit according to the vertical coordinates corresponding to different horizontal coordinates;
the determining the Lorentz curve inversion expression corresponding to different permeability variation coefficients of the heterogeneous oil reservoir comprises the following steps:
and performing function inversion on the Lorentz curves corresponding to different permeability variation coefficients according to the following formula to obtain Lorentz curve inversion expressions corresponding to different permeability variation coefficients:
Figure GDA0002580992260000021
wherein x ismThe abscissa of the Lorentz curve from the 1 st group to the m th group is shown; y ismThe longitudinal coordinates of the Lorentz curve from the 1 st group to the m th group are shown; n is the number of oil reservoir layers; m is less than or equal to n; vKIs the coefficient of permeability variation; kiThe permeability of the reservoir at the ith layer; h isiThe thickness of the ith layer of reservoir;
the method comprises the following steps of determining the abscissa of the coordinate point corresponding to each layer of oil deposit on the Lorentz curve according to the number of oil deposit layers and the thickness of each layer of oil deposit, wherein the abscissa of the coordinate point corresponding to each layer of oil deposit on the Lorentz curve is determined according to the following formula:
Figure GDA0002580992260000022
and determining the permeability of the corresponding coordinate points of each layer of oil deposit according to the vertical coordinates corresponding to different horizontal coordinates, wherein the permeability of the corresponding coordinate points of each layer of oil deposit is determined according to the following formula:
Figure GDA0002580992260000023
the embodiment of the invention also provides an inversion device for permeability distribution of a heterogeneous oil reservoir, which is used for establishing a set of simple and easy-to-use inversion solution for permeability distribution of the heterogeneous oil reservoir with unified standard, and the device comprises:
the expression determining module is used for determining Lorentz curve inversion expressions corresponding to different permeability variation coefficients of the heterogeneous oil reservoir;
the abscissa determining module is used for determining the abscissas of the coordinate points corresponding to each layer of oil deposit on the Lorentz curve according to the number of oil deposit layers and the thickness of each layer of oil deposit;
the ordinate determining module is used for determining the ordinates corresponding to different abscissas by utilizing the Lorentz curve inversion expressions corresponding to different permeability variation coefficients;
the permeability determining module is used for determining the permeability of the corresponding coordinate points of each layer of oil deposit according to the vertical coordinates corresponding to different horizontal coordinates;
the expression determination module is further to:
and performing function inversion on the Lorentz curves corresponding to different permeability variation coefficients according to the following formula to obtain Lorentz curve inversion expressions corresponding to different permeability variation coefficients:
Figure GDA0002580992260000031
wherein x ismThe abscissa of the Lorentz curve from the 1 st group to the m th group is shown; y ismThe longitudinal coordinates of the Lorentz curve from the 1 st group to the m th group are shown; n is the number of oil reservoir layers; m is less than or equal to n; vKIs the coefficient of permeability variation; kiThe permeability of the reservoir at the ith layer; h isiThe thickness of the ith layer of reservoir;
the abscissa determination module is further to:
determining the abscissa of the coordinate point corresponding to each layer of oil deposit on the Lorentz curve according to the following formula:
Figure GDA0002580992260000032
the permeability determination module is further to:
and determining the permeability of the corresponding coordinate points of each layer of oil deposit according to the following formula:
Figure GDA0002580992260000033
the embodiment of the invention also provides computer equipment which comprises a memory, a processor and a computer program which is stored on the memory and can be operated on the processor, wherein the processor realizes the heterogeneous reservoir permeability distribution inversion method when executing the computer program.
An embodiment of the present invention further provides a computer-readable storage medium, where a computer program for executing the heterogeneous reservoir permeability distribution inversion method is stored in the computer-readable storage medium.
The embodiment of the invention solves the problems existing in the current heterogeneous reservoir permeability distribution inversion, establishes a set of simple and easy-to-use heterogeneous reservoir permeability distribution inversion solution with unified standard, improves the understanding of the influence of the heterogeneity of the heterogeneous reservoir on the reservoir development, and can better guide the heterogeneous reservoir development.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts. In the drawings:
FIG. 1 is a schematic diagram of an inversion method of heterogeneous reservoir permeability distribution according to an embodiment of the present disclosure;
FIG. 2 is a schematic diagram of Lorentzian curves in an embodiment of the present invention;
FIG. 3 is a Lorentzian curve diagram illustrating different permeability coefficients of variation according to an embodiment of the present invention;
fig. 4 is a schematic diagram of an inversion apparatus for heterogeneous reservoir permeability distribution in an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the embodiments of the present invention are further described in detail below with reference to the accompanying drawings. The exemplary embodiments and descriptions of the present invention are provided to explain the present invention, but not to limit the present invention.
In order to solve the problems of the current heterogeneous reservoir permeability distribution inversion method, an embodiment of the present invention provides an heterogeneous reservoir permeability distribution inversion method, as shown in fig. 1, where the method may include:
step 101, determining Lorentz curve inversion expressions corresponding to different permeability variation coefficients of the heterogeneous oil reservoir;
102, determining the abscissa of a coordinate point corresponding to each layer of oil deposit on a Lorentz curve according to the number of oil deposit layers and the thickness of each layer of oil deposit;
103, determining vertical coordinates corresponding to different horizontal coordinates by using Lorentz curve inversion expressions corresponding to different permeability coefficient of variation;
and 104, determining the permeability of the corresponding coordinate points of each layer of oil deposit according to the vertical coordinates corresponding to different horizontal coordinates.
The process shown in fig. 1 shows that the embodiment of the invention solves the problems existing in the current heterogeneous reservoir permeability distribution inversion, establishes a set of simple and easy-to-use heterogeneous reservoir permeability distribution inversion solution with uniform standard, improves the understanding of the influence of the heterogeneity of the heterogeneous reservoir on reservoir development, and can better guide the development of the heterogeneous reservoir.
In specific implementation, Lorentz curve inversion expressions corresponding to different permeability variation coefficients of the heterogeneous oil reservoir are determined. In the embodiment, a Lorentz curve inversion expression corresponding to different permeability variation coefficients can be determined by a method for determining the permeability variation coefficient of the heterogeneous oil reservoir based on the Lorentz curve.
Lorentz curves were developed by the American statistician Lorentz and were originally used to study the problem of allocation of national income among the nations. As shown in fig. 2, if the distribution is uniform, the curve is a straight line with a slope of 45 °, i.e., the diagonal line AC in fig. 2, which is called "completely uniform line"; when the distribution of the national total income is concentrated at one point, the curve is a broken line, namely a broken line ADC in FIG. 2, which is called a completely uneven line; the lorentz curve in the general case is intermediate between these two cases, in the form of a "convex up" curve, i.e. curve AEC in fig. 2, and therefore the kini coefficient can be defined as:
Figure GDA0002580992260000051
wherein S isACEAArea of the area of ACEA in FIG. 2; sACDAIs the area of the ACDA region in fig. 2; y is a curve ordinate; x is the abscissa of the curve.
The Lorentz curve is used for evaluating the overall heterogeneity of the oil reservoir, and the permeability variation coefficient (namely the Gini coefficient) is used for evaluating the magnitude of the heterogeneity degree.
In an embodiment, determining the lorentz curve inversion expression corresponding to different permeability coefficient of variation of the heterogeneous reservoir may include:
permeability K of each layer of oil depositiDetermining K according to the sequence from small to largeiCorresponding layer reservoir thickness hiObtaining n groups (K)i，hi) Is calculated as followsm，ym) Obtaining Lorentz curves corresponding to different permeability variation coefficients:
Figure GDA0002580992260000052
Figure GDA0002580992260000053
wherein x ismThe abscissa of the Lorentz curve from the 1 st group to the m th group is shown; y ismThe longitudinal coordinates of the Lorentz curve from the 1 st group to the m th group are shown; n is the number of oil reservoir layers; m is less than or equal to n;
and performing function inversion on the Lorentz curves corresponding to different permeability variation coefficients according to the following formula to obtain Lorentz curve inversion expressions corresponding to different permeability variation coefficients:
Figure GDA0002580992260000054
FIG. 3 illustrates Lorentzian curves for different permeability coefficients of variation in one embodiment. Taking the lorentz curve shown in fig. 2 as an example, in order to quantitatively research the relationship between the reservoir heterogeneity and the development characteristics, according to the permeability variation coefficient calculation method, the curves corresponding to different permeability variation coefficients are subjected to function inversion, so that lorentz curve inversion expressions with corresponding permeability variation coefficients of 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9 and 1.0 in a lorentz curve plate are obtained, and the lorentz curve inversion expressions are shown in table 1:
TABLE 1 Lorentz curve inversion expression
VK xi＜xc xi＞xc xc
0 yi＝xi yi＝xi 0.50
0.1 yi＝-0.37xi 2+1.3409xi yi＝-0.2132xi 2+1.1657xi+0.0475 0.46073
0.2 yi＝-0.88xi 2+1.7423xi yi＝-0.3000xi 2+1.1556xi+0.1444 0.42174
0.3 yi＝-1.66xi 2+2.2619xi yi＝-0.3264xi 2+1.0632xi+0.2632 0.38009
0.4 yi＝-2.88xi 2+2.9395xi yi＝-0.3191xi 2+0.9337xi+0.3854 0.336726
0.5 yi＝-4.95xi 2+3.8785xi yi＝-0.3000xi 2+0.7972xi+0.5028 0.29076
0.6 yi＝-8.85xi 2+5.2769xi yi＝-0.2651xi 2+0.6478xi+0.6173 0.241635
0.7 yi＝-17.36xi 2+7.5697xi yi＝-0.21xi 2+0.4798xi+0.7302 0.18922
0.8 yi＝-41.90xi 2+12.1xi yi＝-0.1438xi 2+0.3156xi+0.8282 0.132454
0.9 yi＝-175.03xi 2+25.535xi yi＝-0.0606xi 2+0.1359xi+0.9247 0.070089
1.0 yi＝0 yi＝1 0
In the examples, when VKTake value in each VKWhen the values are in the middle, any V is determined by using a linear interpolation methodKThe corresponding Lorentzian curve inversion expression has the following interpolation function:
Figure GDA0002580992260000061
similarly, the corresponding inflection point x is determined by using a linear interpolation methodcThe interpolation function is as follows:
Figure GDA0002580992260000062
wherein, VkIs any VKA value; vk0Is a VkIs a previous VKA value; vk1Is a VkLast V ofKThe value is obtained.
For example, take Vk0＝0.1，Vk1When the value is 0.2, then VkThe inflection point value at 0.15 is shown by the following equation:
Figure GDA0002580992260000071
when x < xcAnd then, the corresponding Lorentz curve inversion expression is as follows:
Figure GDA0002580992260000072
when x > xcAnd then, the corresponding Lorentz curve inversion expression is as follows:
Figure GDA0002580992260000073
and during specific implementation, determining the abscissa of the coordinate point corresponding to each layer of oil deposit on the Lorentz curve according to the number of the oil deposit layers of the heterogeneous oil deposit and the thickness of each layer of oil deposit. In an embodiment, the abscissa of the coordinate point corresponding to each layer of the reservoir on the lorentz curve can be determined according to the following formula:
Figure GDA0002580992260000074
wherein x ismThe abscissa of the Lorentz curve from the 1 st group to the m th group is shown; h isiFor each layer of reservoir permeability KiArranged in order from small to large, KiThe thickness of the corresponding layer oil deposit; n is the number of oil reservoir layers; m is less than or equal to n. In the embodiment, the number n of reservoir layers and the thickness h of each layer of reservoir required for research can be given first1，h2，h3……hnDetermining the arrangement rule of the thicknesses of the layers according to experimental or research requirements, and determining the abscissa x by using the formula (10)m。
In specific implementation, the Lorentz curve inversion expressions corresponding to different permeability coefficient of variation are utilized to determine the vertical coordinates corresponding to different horizontal coordinates. In an embodiment, the permeability coefficient of variation V may be given based onKApplication of VKCorresponding Lorentz curve inversion expression, calculating different abscissa xmLower corresponding ordinate ymThe value of (c).
And during specific implementation, determining the permeability of the corresponding coordinate points of each layer of oil deposit according to the vertical coordinates corresponding to different horizontal coordinates. In an embodiment, the permeability of the corresponding coordinate point of each layer of oil reservoir may be determined according to the following formula:
Figure GDA0002580992260000075
wherein, ymThe longitudinal coordinates of the Lorentz curve from the 1 st group to the m th group are shown; n is the number of oil reservoir layers; m is less than or equal to n, KiThe permeability of each layer of oil deposit is arranged from small to large.
As described aboveIn an embodiment, based on the lorentz curve inversion expression corresponding to different permeability coefficient of variation, given the number of reservoir layers and the thickness of each layer of reservoir required for research, the coordinate point (x) corresponding to each layer on the lorentz curve can be determined1，y1)，(x2，y2)，(x3，y3)……(xn，yn) Thereby reversing the permeability value K of each layer of the reservoir to be researched1，K2，K3……Kn。
Specific examples are given below to further explain the heterogeneous reservoir permeability distribution inversion method of the embodiment of the invention in detail.
Example 1:
in this example, permeability distribution with permeability variation coefficients of 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8 and 0.9 was calculated using the lorentz curve inversion expression of the embodiment of the present invention, the total number of oil reservoirs was 5, the thicknesses of the layers were 5m, and the average permeability K wasFlat plateIs 100 md. Table 2 is calculation data of the inversion result of permeability distribution in this example, and it can be known from table 2 that the embodiment of the present invention can invert permeability distribution under different heterogeneous conditions, thereby further studying the influence of heterogeneous intensity on reservoir development.
TABLE 2 inversion results of permeability distributions under different heterogeneous conditions
VK KFlat plate，md K1，md K2，md K3，md K4，md K5，md
0.1 100 126.69 111.89 96.504 86.722 78.194
0.2 100 156.63 121.43 86.82 73.56 61.56
0.3 100 192.99 125.138 73.68 60.624 47.568
0.4 100 236.35 117.562 61.46 48.696 35.932
0.5 100 288.85 97.99 49.72 37.72 25.72
0.6 100 350.69 66.312 38.27 27.666 17.062
0.7 100 408.88 35.38 26.98 18.58 10.18
0.8 100 442.784 22.932 17.18 11.428 5.676
0.9 100 474.728 9.954 7.53 5.106 2.682
Example 2:
in this example, permeability distributions when permeability variation coefficients are 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8 and 0.9 respectively are calculated by using a lorentz curve inversion expression of the embodiment of the invention, the total number of oil reservoirs is 10, the thicknesses of all layers are 5m, and the average permeability K isFlat platePermeability distribution at 300 md. Table 3 is calculation data of the inversion result of permeability distribution in this example, and it can be known from table 3 that the embodiment of the present invention can invert permeability distribution under different heterogeneous conditions, thereby further studying the influence of heterogeneous intensity on reservoir development.
TABLE 3 inversion results of permeability distributions under different heterogeneous conditions
Figure GDA0002580992260000081
Figure GDA0002580992260000091
Example 3:
in this example, the permeability distribution when the permeability variation coefficient is 0.15 is calculated by using the lorentz curve inversion expression of the embodiment of the invention, the total number of layers of the oil reservoir is 6, the thickness of each layer is 5m, and the average permeability K isFlat plateIs 200 md. Table 4 is the calculation data of the inversion result of permeability distribution in this example, and it can be known from table 4 that the embodiment of the present invention can invert permeability distribution under different heterogeneous conditions, thereby laying a foundation for studying the influence of heterogeneous intensity on reservoir development.
TABLE 4 permeability distribution inversion results
VK KFlat plate，md K1，md K2，md K3，md K4，md K5，md K6，md
0.15 200 287.64 246.24 200.64 172.32 155.11 138.05
Based on the same inventive concept, the embodiment of the invention also provides an inversion device for permeability distribution of a heterogeneous reservoir, as described in the following embodiments. Because the principle of the device for solving the problems is similar to the heterogeneous reservoir permeability distribution inversion method, the implementation of the device can refer to the implementation of the heterogeneous reservoir permeability distribution inversion method, and repeated parts are not repeated.
Fig. 4 is a schematic diagram of an apparatus for inverting permeability distribution of a heterogeneous reservoir according to an embodiment of the present invention, and as shown in fig. 4, the apparatus may include:
the expression determining module 401 is configured to determine lorentz curve inversion expressions corresponding to different permeability variation coefficients of the heterogeneous oil reservoir;
the abscissa determining module 402 is configured to determine an abscissa of a coordinate point corresponding to each layer of the oil deposit on the lorentz curve according to the number of oil deposit layers and the thickness of each layer of the oil deposit;
a ordinate determining module 403, configured to determine, by using the lorentz curve inversion expressions corresponding to different permeability variation coefficients, ordinates corresponding to different abscissas;
and a permeability determining module 404, configured to determine, according to the ordinate corresponding to the different abscissa, the permeability of the corresponding coordinate point of each layer of the oil reservoir.
In one embodiment, the expression determination module 401 may be further configured to:
permeability K of each layer of oil depositiDetermining K according to the sequence from small to largeiCorresponding layer reservoir thickness hiObtaining n groups (K)i，hi) Is calculated as followsm，ym) Obtaining Lorentz curves corresponding to different permeability variation coefficients:
Figure GDA0002580992260000101
Figure GDA0002580992260000102
and performing function inversion on the Lorentz curves corresponding to different permeability variation coefficients according to the following formula to obtain Lorentz curve inversion expressions corresponding to different permeability variation coefficients:
Figure GDA0002580992260000103
wherein x ismThe abscissa of the Lorentz curve from the 1 st group to the m th group is shown; y ismThe longitudinal coordinates of the Lorentz curve from the 1 st group to the m th group are shown; n is the number of oil reservoir layers; m is less than or equal to n; vKThe permeability coefficient of variation.
In one embodiment, the expression determination module 401 may be further configured to:
when V isKTake value in each VKWhen the values are in the middle, any V is determined by using a linear interpolation methodKThe corresponding Lorentzian curve inversion expression has the following interpolation function:
Figure GDA0002580992260000111
and determining the corresponding inflection point x by using a linear interpolation methodcThe interpolation function is as follows:
Figure GDA0002580992260000112
wherein, VkIs any VKA value; vk0Is a VkIs a previous VKA value; vk1Is a VkLast V ofKThe value is obtained.
In one embodiment, the abscissa determination module 402 may be further configured to:
determining the abscissa of the coordinate point corresponding to each layer of oil deposit on the Lorentz curve according to the following formula:
Figure GDA0002580992260000113
in one embodiment, the permeability determination module 404 may be further configured to:
and determining the permeability of the corresponding coordinate points of each layer of oil deposit according to the following formula:
Figure GDA0002580992260000114
the embodiment of the invention also provides computer equipment which comprises a memory, a processor and a computer program which is stored on the memory and can be operated on the processor, wherein the processor realizes the heterogeneous reservoir permeability distribution inversion method when executing the computer program.
An embodiment of the present invention further provides a computer-readable storage medium, where a computer program for executing the heterogeneous reservoir permeability distribution inversion method is stored in the computer-readable storage medium.
In conclusion, the embodiment of the invention solves the problems existing in the current heterogeneous reservoir permeability distribution inversion, establishes a set of simple and easy-to-use heterogeneous reservoir permeability distribution inversion solution with unified standard, improves the understanding of the influence of the heterogeneity of the heterogeneous reservoir on reservoir development, and can better guide the heterogeneous reservoir development.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.

Claims (6)

1. An inversion method for permeability distribution of a heterogeneous reservoir is characterized by comprising the following steps:
determining Lorentz curve inversion expressions corresponding to different permeability variation coefficients of the heterogeneous oil reservoir;
determining the abscissa of a coordinate point corresponding to each layer of oil deposit on the Lorentz curve according to the number of oil deposit layers and the thickness of each layer of oil deposit;
determining the vertical coordinates corresponding to different horizontal coordinates by using the Lorentz curve inversion expression corresponding to different permeability variation coefficients;
determining the permeability of the corresponding coordinate points of each layer of oil deposit according to the vertical coordinates corresponding to different horizontal coordinates;
the determining the Lorentz curve inversion expression corresponding to different permeability variation coefficients of the heterogeneous oil reservoir comprises the following steps:
and performing function inversion on the Lorentz curves corresponding to different permeability variation coefficients according to the following formula to obtain Lorentz curve inversion expressions corresponding to different permeability variation coefficients:
Figure FDA0002779786110000011
wherein x ismThe abscissa of the Lorentz curve from the 1 st group to the m th group is shown; y ismThe longitudinal coordinates of the Lorentz curve from the 1 st group to the m th group are shown; n is the number of oil reservoir layers; m is less than or equal to n; vKIs the coefficient of permeability variation; kiThe permeability of the reservoir at the ith layer; h isiThe thickness of the ith layer of reservoir;
the method comprises the following steps of determining the abscissa of the coordinate point corresponding to each layer of oil deposit on the Lorentz curve according to the number of oil deposit layers and the thickness of each layer of oil deposit, wherein the abscissa of the coordinate point corresponding to each layer of oil deposit on the Lorentz curve is determined according to the following formula:
Figure FDA0002779786110000012
and determining the permeability of the corresponding coordinate points of each layer of oil deposit according to the vertical coordinates corresponding to different horizontal coordinates, wherein the permeability of the corresponding coordinate points of each layer of oil deposit is determined according to the following formula:
Figure FDA0002779786110000013
2. the method of claim 1, wherein determining the Lorentzian curve inversion expression corresponding to different permeability coefficients of variation for the heterogeneous reservoir comprises:
when V iskTake value in each VKApplying linearity between valuesInterpolation method for determining arbitrary VkThe corresponding Lorentzian curve inversion expression has the following interpolation function:
Figure FDA0002779786110000021
and determining the corresponding inflection point x by using a linear interpolation methodcThe interpolation function is as follows:
Figure FDA0002779786110000022
wherein, V is more than 0k< 1 and Vk≠0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9；VK∈{0,0.1,0.2,…,1.0}；Vk0Is a VkIs a previous VKA value; vk1Is a VkLast V ofKThe value is obtained.
3. An inversion device for permeability distribution of heterogeneous reservoirs, comprising:
the expression determining module is used for determining Lorentz curve inversion expressions corresponding to different permeability variation coefficients of the heterogeneous oil reservoir;
the abscissa determining module is used for determining the abscissas of the coordinate points corresponding to each layer of oil deposit on the Lorentz curve according to the number of oil deposit layers and the thickness of each layer of oil deposit;
the ordinate determining module is used for determining the ordinates corresponding to different abscissas by utilizing the Lorentz curve inversion expressions corresponding to different permeability variation coefficients;
the permeability determining module is used for determining the permeability of the corresponding coordinate points of each layer of oil deposit according to the vertical coordinates corresponding to different horizontal coordinates;
the expression determination module is further to:
and performing function inversion on the Lorentz curves corresponding to different permeability variation coefficients according to the following formula to obtain Lorentz curve inversion expressions corresponding to different permeability variation coefficients:
Figure FDA0002779786110000023
wherein x ismThe abscissa of the Lorentz curve from the 1 st group to the m th group is shown; y ismThe longitudinal coordinates of the Lorentz curve from the 1 st group to the m th group are shown; n is the number of oil reservoir layers; m is less than or equal to n; vKIs the coefficient of permeability variation; kiThe permeability of the reservoir at the ith layer; h isiThe thickness of the ith layer of reservoir;
the abscissa determination module is further to:
determining the abscissa of the coordinate point corresponding to each layer of oil deposit on the Lorentz curve according to the following formula:
Figure FDA0002779786110000031
the permeability determination module is further to:
and determining the permeability of the corresponding coordinate points of each layer of oil deposit according to the following formula:
Figure FDA0002779786110000032
4. the apparatus of claim 3, wherein the expression determination module is further to:
when V iskTake value in each VKWhen the values are in the middle, any V is determined by using a linear interpolation methodkThe corresponding Lorentzian curve inversion expression has the following interpolation function:
Figure FDA0002779786110000033
and determining the corresponding inflection point x by using a linear interpolation methodcThe interpolation function is as follows:
Figure FDA0002779786110000034
wherein, V is more than 0k< 1 and Vk≠0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9；VK∈{0,0.1,0.2,…,1.0}；Vk0Is a VkIs a previous VKA value; vk1Is a VkLast V ofKThe value is obtained.
5. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method of any of claims 1 to 2 when executing the computer program.
6. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program for executing the method of any one of claims 1 to 2.
CN201810187730.5A 2018-03-07 2018-03-07 Heterogeneous reservoir permeability distribution inversion method and device Active CN108491361B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810187730.5A CN108491361B (en) 2018-03-07 2018-03-07 Heterogeneous reservoir permeability distribution inversion method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810187730.5A CN108491361B (en) 2018-03-07 2018-03-07 Heterogeneous reservoir permeability distribution inversion method and device

Publications (2)

Publication Number Publication Date
CN108491361A CN108491361A (en) 2018-09-04
CN108491361B true CN108491361B (en) 2021-02-23

Family

ID=63341819

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810187730.5A Active CN108491361B (en) 2018-03-07 2018-03-07 Heterogeneous reservoir permeability distribution inversion method and device

Country Status (1)

Country Link
CN (1) CN108491361B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110242291B (en) * 2019-07-02 2022-03-01 中国石油化工股份有限公司 Method for representing heterogeneity of super-strong heterogeneous oil and gas reservoir
CN112347707B (en) * 2019-08-07 2022-10-04 中国石油天然气股份有限公司 Method and device for characterizing flow heterogeneity of oil reservoir
CN111456710A (en) * 2020-04-27 2020-07-28 中国石油天然气股份有限公司 Profile control well selection method

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105525909A (en) * 2015-11-19 2016-04-27 薛云飞 Method for analyzing heterogeneous property of oil reservoir
CN107219159A (en) * 2017-04-19 2017-09-29 中国石油天然气股份有限公司 A kind of method and apparatus for determining reservoir heterogeneity

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140068353A1 (en) * 2012-08-31 2014-03-06 The Government Of The United States Of America, As Represented By The Secretary Of The Navy Estimation of Hidden Variance Distribution Parameters

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105525909A (en) * 2015-11-19 2016-04-27 薛云飞 Method for analyzing heterogeneous property of oil reservoir
CN107219159A (en) * 2017-04-19 2017-09-29 中国石油天然气股份有限公司 A kind of method and apparatus for determining reservoir heterogeneity

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
储层渗透率非均质性表征新参数;邵先杰等;《石油实验地质》;20100831;第32卷(第4期);全文 *
洛伦兹曲线在微观非均质性研究中的应用;林伶等;《复杂油气藏》;20160630;第9卷(第2期);全文 *
洛伦茨曲线在油藏工程中的应用;康晓东等;《新疆石油地质》;20020228;第23卷(第1期);全文 *
渗透率变异系数的几种计算方法;朱小影等;《海洋石油》;20090622;第29卷(第2期);全文 *
非均质油藏开发规律研究;王增林等;《油气地质与采收率》;20110930;第18卷(第5期);全文 *

Also Published As

Publication number Publication date
CN108491361A (en) 2018-09-04

Similar Documents

Publication Publication Date Title
CN108491361B (en) Heterogeneous reservoir permeability distribution inversion method and device
CN111681162B (en) Defect sample generation method and device, electronic equipment and storage medium
CN112101530A (en) Neural network training method, device, equipment and storage medium
CN104809311A (en) Structural part remaining life predicting method based on multi-factor fusion correction
CN112182910B (en) Method and device for controlling preparation effect of cross-linked material
CN114757587B (en) Product quality control system and method based on big data
CN112506752A (en) Operation and maintenance method, device and equipment based on time series data trend prediction
CN114444815B (en) Defect-plate stock layout method and system based on branch cutting
CN112631881A (en) Operation and maintenance method, device and equipment based on time series data trend prediction
CN111783209B (en) Self-adaptive structure reliability analysis method combining learning function and kriging model
CN107346350B (en) Distribution method, device and cluster system for integrated circuit layout data processing tasks
CN105224558A (en) The evaluation disposal route of speech business and device
CN110633529A (en) Data processing method, device and system for determining oil reservoir water drive development effect
CN112289379B (en) Method and device for determining cell type, storage medium and electronic device
CN110751400B (en) Risk assessment method and device
CN107291843A (en) Hierarchical clustering improved method based on Distributed Computing Platform
CN110543691B (en) Planar fracture grouting quasi-three-dimensional simulation method and system based on finite volume method
CN110795255B (en) Data batch value adjusting method and device, readable storage medium and equipment
CN114706303A (en) Control method and device for strip steel width allowance, medium and electronic equipment
CN112419493B (en) Shale reservoir three-dimensional attribute model building method and device
CN113568733A (en) Resource allocation method, device, electronic equipment and storage medium
CN113900800B (en) Distribution method of edge computing system
CN110570058B (en) Method and device for determining sensitivity
KR102265361B1 (en) Method and apparatus for deep-learning using smoothing
CN109885981A (en) A kind of fatigue analysis method based on design margin evaluation

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant