CN108491361A - heterogeneous reservoir permeability distribution inversion method and device - Google Patents
heterogeneous reservoir permeability distribution inversion method and device Download PDFInfo
- Publication number
- CN108491361A CN108491361A CN201810187730.5A CN201810187730A CN108491361A CN 108491361 A CN108491361 A CN 108491361A CN 201810187730 A CN201810187730 A CN 201810187730A CN 108491361 A CN108491361 A CN 108491361A
- Authority
- CN
- China
- Prior art keywords
- reservoir
- permeability
- lorentz curve
- coefficient
- oil 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.)
- Granted
Links
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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Algebra (AREA)
- Evolutionary Biology (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Investigating Or Analysing Biological Materials (AREA)
Abstract
The invention discloses a kind of heterogeneous reservoir Permeability Distribution inversion method and device, wherein this method includes：Determine the corresponding Lorentz curve inverting expression formula of heterogeneous reservoir difference coefficient of permeability variation；According to the oil reservoir number of plies of heterogeneous reservoir and each layer reservoir thickness, the abscissa of each layer oil reservoir respective coordinates point on Lorentz curve is determined；Using the corresponding Lorentz curve inverting expression formula of different coefficient of permeability variation, the corresponding ordinate of different abscissas is determined；According to the corresponding ordinate of different abscissas, the permeability of the respective coordinates point of each layer oil reservoir is determined.The present invention solves the problems, such as that current heterogeneous reservoir Permeability Distribution inverting exists, establish a set of simplicity it is easy-to-use and with unified standard heterogeneous reservoir Permeability Distribution inverting solution, the understanding influenced on oil reservoir development on the anisotropism of heterogeneous reservoir is improved, preferably heterogeneous reservoir can be instructed to develop.
Description
Technical field
The present invention relates to the heterogeneous reservoir development technique field of petroleum industry more particularly to heterogeneous reservoir permeabilities point
Cloth inversion method and device.
Background technology
Oil reservoir generally existing interlayer heterogeneity, the common parameter for weighing interlayer heterogeneity includes differential, advancing by leaps and bounds is
Number, coefficient of variation etc., wherein the coefficient of variation is the most commonly used, it reflects that mono-layer osmotic rate deviates the degree of ensemble average value, is to comment
The most important parameter of valence oil reservoir macroscopic heterogeneity, coefficient of permeability variation VKBigger expression Reservoir Heterogeneity is more serious.Mesh
Before be usually used in calculate the coefficient of variation method have statistical method, Dykstra H methods and Lorentz method.Statistical method is suitable
For the reservoir of any Permeability Distribution type, but the value for coefficient of variation that this method calculates is between 0~∞, to non-homogeneous degree
Measurement lack relativity；Dykstra H methods are a kind of empirical methods, the value for coefficient of variation that this method calculates 0~1 it
Between, but it is only applicable to the reservoir that permeability is logarithm normal distribution；Lorentz method then has good applicability, it is not only fitted
For the reservoir of any Permeability Distribution type and the value for coefficient of variation of calculating is between 0~1.
However, although Lorentz method can calculate the coefficient of variation of a certain determining oil reservoir, its anisotropism is analyzed,
Since research needs, it is sometimes desirable to the heterogeneous reservoir for building some or certain several specific coefficient of variation, come study it is non-
Influence of the matter to effect of reservoir development.But lacks more efficiently inversion method now and go to obtain corresponding permeability point
Cloth.It is directed to the anisotropism research of the oil reservoir of some or certain several specific coefficient of variation at present, is more by tentative calculation
Method obtains, and lacks unified standard using the Permeability Distribution that this method obtains.
Invention content
The embodiment of the present invention provides a kind of heterogeneous reservoir Permeability Distribution inversion method, easy-to-use to establish a set of simplicity
And there is the heterogeneous reservoir Permeability Distribution inverting solution of unified standard, this method includes：
Determine the corresponding Lorentz curve inverting expression formula of heterogeneous reservoir difference coefficient of permeability variation；
According to the oil reservoir number of plies of heterogeneous reservoir and each layer reservoir thickness, determine that each layer oil reservoir corresponds to seat on Lorentz curve
The abscissa of punctuate；
Using the corresponding Lorentz curve inverting expression formula of different coefficient of permeability variation, determine that different abscissas are corresponding
Ordinate；
According to the corresponding ordinate of different abscissas, the permeability of the respective coordinates point of each layer oil reservoir is determined.
The embodiment of the present invention also provides a kind of heterogeneous reservoir Permeability Distribution inverting device, a set of easy to be easy to establish
With and with unified standard heterogeneous reservoir Permeability Distribution inverting solution, which includes：
Expression formula determining module, for determining that the corresponding Lorentz curve of heterogeneous reservoir difference coefficient of permeability variation is anti-
Drill expression formula；
Abscissa determining module is used for the oil reservoir number of plies according to heterogeneous reservoir and each layer reservoir thickness, determines Lorentz
The abscissa of each layer oil reservoir respective coordinates point on curve；
Ordinate determining module, for utilizing the corresponding Lorentz curve inverting expression formula of different coefficient of permeability variation,
Determine the corresponding ordinate of different abscissas；
Permeability determining module, for according to the corresponding ordinate of different abscissas, determining the respective coordinates of each layer oil reservoir
The permeability of point.
The embodiment of the present invention also provides a kind of computer equipment, including memory, processor and storage are on a memory simultaneously
The computer program that can be run on a processor, the processor realize above-mentioned heterogeneous reservoir when executing the computer program
Permeability Distribution inversion method.
The embodiment of the present invention also provides a kind of computer readable storage medium, and the computer-readable recording medium storage has
Execute the computer program of above-mentioned heterogeneous reservoir Permeability Distribution inversion method.
The embodiment of the present invention solves the problems, such as that current heterogeneous reservoir Permeability Distribution inverting exists, and establishes a set of letter
It is easy to use and with the heterogeneous reservoir Permeability Distribution inverting solution of unified standard, improve to the non-of heterogeneous reservoir
The understanding that homogenieity influences oil reservoir development preferably can instruct heterogeneous reservoir to develop.
Description of the drawings
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technology description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with
Obtain other attached drawings according to these attached drawings.In the accompanying drawings：
Fig. 1 is the schematic diagram of heterogeneous reservoir Permeability Distribution inversion method in the embodiment of the present invention；
Fig. 2 is the schematic diagram of Lorentz curve in the embodiment of the present invention；
Fig. 3 is the corresponding Lorentz curve schematic diagram of difference coefficient of permeability variation in the embodiment of the present invention；
Fig. 4 is the schematic diagram of heterogeneous reservoir Permeability Distribution inverting device in the embodiment of the present invention.
Specific implementation mode
Understand in order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the accompanying drawings to this hair
Bright embodiment is described in further details.Here, the illustrative embodiments of the present invention and their descriptions are used to explain the present invention, but simultaneously
It is not as a limitation of the invention.
In order to solve the problems, such as that current heterogeneous reservoir Permeability Distribution inversion method exists, provided in the embodiment of the present invention
A kind of heterogeneous reservoir Permeability Distribution inversion method, as shown in Figure 1, this method may include：
Step 101 determines the corresponding Lorentz curve inverting expression formula of heterogeneous reservoir difference coefficient of permeability variation；
Step 102, the oil reservoir number of plies according to heterogeneous reservoir and each layer reservoir thickness determine each layer oil on Lorentz curve
Hide the abscissa of respective coordinates point；
Step 103, using the corresponding Lorentz curve inverting expression formula of different coefficient of permeability variation, determine different horizontal seats
Mark corresponding ordinate；
Step 104, according to the corresponding ordinate of different abscissas, determine the permeability of the respective coordinates point of each layer oil reservoir.
Flow is it is known that the embodiment of the present invention solves current heterogeneous reservoir Permeability Distribution inverting as shown in Figure 1
There are the problem of, establish a set of simplicity it is easy-to-use and with unified standard heterogeneous reservoir Permeability Distribution inverting solution party
Case improves the understanding influenced on oil reservoir development on the anisotropism of heterogeneous reservoir, can preferably instruct heterogeneous reservoir
Exploitation.
When it is implemented, determining the corresponding Lorentz curve inverting expression of heterogeneous reservoir difference coefficient of permeability variation
Formula.The method that heterogeneous reservoir coefficient of permeability variation can be determined based on Lorentz curve in embodiment, determines different infiltrations
The corresponding Lorentz curve inverting expression formula of the rate coefficient of variation.
Lorentz curve is proposed by U.S.'s statistician's Lorentz, is initially used to study national income its people's
Between assignment problem.As shown in Fig. 2, if distribution is uniform, which is the straight line that a slope is 45 °, i.e. in Fig. 2
Diagonal line AC, referred to as " substantially uniformity line "；When the distribution that its people all take in converges, which is broken line, i.e. Fig. 2
Middle polyline ADC, referred to as " complete uneven line "；Lorentz curve under normal circumstances between both of these case, for " on
It is convex " form of curve, i.e. curve AEC in Fig. 2, therefore, Gini coefficient may be defined as：
Wherein, SACEAFor the area in the regions ACEA in Fig. 2；SACDAFor the area in the regions ACDA in Fig. 2；Y, which is that curve is vertical, to be sat
Mark；X is curvilinear abscissa.
Lorentz curve is commented for evaluating oil reservoir entirety anisotropism with coefficient of permeability variation (i.e. Gini coefficient)
The size of valence non-homogeneous degree.
In embodiment, the corresponding Lorentz curve inverting expression formula of heterogeneous reservoir difference coefficient of permeability variation is determined,
May include：
By each layer Reservoir Permeability KiBy being ranked sequentially from small to large, K is determinediRespective layer reservoir thickness hi, obtain n groups
(Ki, hi), (x is calculated as followsm, ym), obtain the corresponding Lorentz curve of different coefficient of permeability variation：
Wherein, xmFor the 1st group of Lorentz curve to m group abscissas；ymIt is sat to m groups are vertical for the 1st group for Lorentz curve
Mark；N is the oil reservoir number of plies；m≤n；
And function inversion is carried out to the corresponding Lorentz curve of different coefficient of permeability variation as follows, it obtains
The corresponding Lorentz curve inverting expression formula of different coefficient of permeability variation：
Fig. 3 shows the corresponding Lorentz curve of difference coefficient of permeability variation in an embodiment.With Lorentz shown in Fig. 2
For curve, for the relationship of quantitative study Reservoir Heterogeneity and development features, calculated according to above-mentioned coefficient of permeability variation
Method, by carrying out function inversion to the corresponding curve of different coefficient of permeability variation, to obtain in Lorentz curve plate
The corresponding permeability coefficient of variation is respectively 0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0 Lorentz curve
Inverting expression formula, as shown in table 1：
1 Lorentz curve inverting expression formula of table
In embodiment, work as VKValue is in each VKWhen between value, arbitrary V is determined with linear interpolation methodKCorresponding Lorentz
Log inversion expression formula, interpolating function are as follows：
Similarly, corresponding inflection point x is determined with linear interpolation methodcValue, interpolating function is as follows：
Wherein, VkFor arbitrary VKValue；Vk0For VkPrevious VKValue；Vk1For VkLatter VKValue.
For example, taking Vk0=0.1, Vk1=0.2, then VkKnee value when=0.15 is shown below：
As x ＜ xcWhen, corresponding Lorentz curve inverting expression formula is：
As x ＞ xcWhen, corresponding Lorentz curve inverting expression formula is：
When it is implemented, according to the oil reservoir number of plies of heterogeneous reservoir and each layer reservoir thickness, determine each on Lorentz curve
The abscissa of layer oil reservoir respective coordinates point.In embodiment, it can determine that each layer oil reservoir corresponds on Lorentz curve as follows
The abscissa of coordinate points：
Wherein, xmFor the 1st group of Lorentz curve to m group abscissas；hiFor each layer Reservoir Permeability KiBy from small to large
It is ranked sequentially, KiRespective layer reservoir thickness；N is the oil reservoir number of plies；m≤n.In embodiment, the desired oil of research can be provided first
Hide number of plies n and each layer reservoir thickness h1, h2, h3……hn, require to determine that the arrangement of each layer thickness is advised according to experiment or research
Then, abscissa x is determined with formula (10)m。
When it is implemented, using the corresponding Lorentz curve inverting expression formula of different coefficient of permeability variation, determine different
The corresponding ordinate of abscissa.It, can be based on given coefficient of permeability variation V in embodimentK, with VKCorresponding Lorentz
Log inversion expression formula calculates different abscissa xmUnder corresponding ordinate ymValue.
When it is implemented, according to the corresponding ordinate of different abscissas, the infiltration of the respective coordinates point of each layer oil reservoir is determined
Rate.In embodiment, the permeability of the respective coordinates point of each layer oil reservoir can be determined as follows：
Wherein, ymFor the 1st group of Lorentz curve to m group ordinates；N is the oil reservoir number of plies；M≤n, KiIt is oozed for each layer oil reservoir
Saturating rate, by being ranked sequentially from small to large.
As described above, in embodiment, the corresponding Lorentz curve inverting expression of different coefficient of permeability variation can be based on
Formula gives the oil reservoir number of plies and each layer reservoir thickness needed for research, determines the corresponding coordinate points (x of each layer on Lorentz curve1,
y1), (x2, y2), (x3, y3)……(xn, yn), thus it is finally inversed by each layer permeability value K of studied oil reservoir1, K2, K3……Kn。
It is given below specific embodiment the present invention will be described in detail again the heterogeneous reservoir Permeability Distribution inverting side of embodiment
Method.
Embodiment 1：
In this example, calculating coefficient of permeability variation with the Lorentz curve inverting expression formula of the embodiment of the present invention is respectively
0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9 when Permeability Distribution, the total number of plies of oil reservoir be 5 layers, each layer thickness is equal
For 5m, mean permeability KIt is flatFor 100md.Table 2 is the calculating data of Permeability Distribution inversion result in this example, can be obtained by table 2
Know, the embodiment of the present invention can carry out inverting to Permeability Distribution under the conditions of different heterogeneity, heterogeneous to further study
Influence of the power to oil reservoir development.
Permeability Distribution inversion result under the conditions of 2 different heterogeneity of table
| VK | KIt is flat, 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 |
Embodiment 2：
In this example, coefficient of permeability variation difference is calculated with the Lorentz curve inverting expression formula of the embodiment of the present invention
0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9 when Permeability Distribution, the total number of plies of oil reservoir be 10 layers, each layer thickness is equal
For 5m, mean permeability KIt is flatFor 300md when Permeability Distribution.Table 3 is the calculating number of Permeability Distribution inversion result in this example
According to, by table 3 it is known that the embodiment of the present invention can carry out inverting to Permeability Distribution under the conditions of different heterogeneity, thus into
One step studies the heterogeneous strong and weak influence to oil reservoir development.
Permeability Distribution inversion result under the conditions of 3 different heterogeneity of table
Embodiment 3：
In this example, it is 0.15 to calculate coefficient of permeability variation with the Lorentz curve inverting expression formula of the embodiment of the present invention
When Permeability Distribution, the total number of plies of oil reservoir be 6 layers, each layer thickness be 5m, mean permeability KIt is flatFor 200md.Table 4 is in this example
The calculating data of Permeability Distribution inversion result, by table 4 it is known that the embodiment of the present invention can be to different heterogeneity under the conditions of
Permeability Distribution carries out inverting, to lay the foundation to study the heterogeneous strong and weak influence to oil reservoir development.
4 Permeability Distribution inversion result of table
| VK | KIt is flat, 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 same inventive concept, a kind of heterogeneous reservoir Permeability Distribution inverting dress is additionally provided in the embodiment of the present invention
It sets, as described in the following examples.The principle and heterogeneous reservoir Permeability Distribution inversion method solved the problems, such as due to the device
It is similar, therefore the implementation of the device may refer to the implementation of heterogeneous reservoir Permeability Distribution inversion method, repeat place no longer
It repeats.
Fig. 4 is the schematic diagram of heterogeneous reservoir Permeability Distribution inverting device in the embodiment of the present invention, as shown in figure 4, should
Device may include：
Expression formula determining module 401, for determining that the corresponding Lorentz of heterogeneous reservoir difference coefficient of permeability variation is bent
Line inversion expression formula；
Abscissa determining module 402 is used for the oil reservoir number of plies according to heterogeneous reservoir and each layer reservoir thickness, determines long-range navigation
Hereby on curve each layer oil reservoir respective coordinates point abscissa；
Ordinate determining module 403, for being expressed using the corresponding Lorentz curve inverting of different coefficient of permeability variation
Formula determines the corresponding ordinate of different abscissas；
Permeability determining module 404, for according to the corresponding ordinate of different abscissas, determining that the corresponding of each layer oil reservoir is sat
The permeability of punctuate.
In one embodiment, expression formula determining module 401 can be further used for：
By each layer Reservoir Permeability KiBy being ranked sequentially from small to large, K is determinediRespective layer reservoir thickness hi, obtain n groups
(Ki, hi), (x is calculated as followsm, ym), obtain the corresponding Lorentz curve of different coefficient of permeability variation：
Function inversion is carried out to the corresponding Lorentz curve of different coefficient of permeability variation as follows, difference is obtained and oozes
The corresponding Lorentz curve inverting expression formula of the rate coefficient of variation thoroughly：
Wherein, xmFor the 1st group of Lorentz curve to m group abscissas；ymIt is sat to m groups are vertical for the 1st group for Lorentz curve
Mark；N is the oil reservoir number of plies；m≤n；VKFor coefficient of permeability variation.
In one embodiment, expression formula determining module 401 can be further used for：
Work as VKValue is in each VKWhen between value, arbitrary V is determined with linear interpolation methodKCorresponding Lorentz curve inverting
Expression formula, interpolating function are as follows：
And determine corresponding inflection point x with linear interpolation methodcValue, interpolating function is as follows：
Wherein, VkFor arbitrary VKValue；Vk0For VkPrevious VKValue；Vk1For VkLatter VKValue.
In one embodiment, abscissa determining module 402 can be further used for：
The abscissa of each layer oil reservoir respective coordinates point on Lorentz curve is determined as follows：
In one embodiment, permeability determining module 404 can be further used for：
The permeability of the respective coordinates point of each layer oil reservoir is determined as follows：
The embodiment of the present invention also provides a kind of computer equipment, including memory, processor and storage are on a memory simultaneously
The computer program that can be run on a processor, the processor realize above-mentioned heterogeneous reservoir when executing the computer program
Permeability Distribution inversion method.
The embodiment of the present invention also provides a kind of computer readable storage medium, and the computer-readable recording medium storage has
Execute the computer program of above-mentioned heterogeneous reservoir Permeability Distribution inversion method.
In conclusion the embodiment of the present invention solves the problems, such as that current heterogeneous reservoir Permeability Distribution inverting exists, build
Stood a set of simplicity it is easy-to-use and with unified standard heterogeneous reservoir Permeability Distribution inverting solution, improve to it is non-
The understanding that the anisotropism of matter oil reservoir influences oil reservoir development preferably can instruct heterogeneous reservoir to develop.
It should be understood by those skilled in the art that, the embodiment of the present invention can be provided as method, system or computer program
Product.Therefore, complete hardware embodiment, complete software embodiment or reality combining software and hardware aspects can be used in the present invention
Apply the form of example.Moreover, the present invention can be used in one or more wherein include computer usable program code computer
The computer program production implemented in usable storage medium (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.)
The form of product.
The present invention be with reference to according to the method for the embodiment of the present invention, the flow of equipment (system) and computer program product
Figure and/or block diagram describe.It should be understood that can be realized by computer program instructions every first-class in flowchart and/or the block diagram
The combination of flow and/or box in journey and/or box and flowchart and/or the block diagram.These computer programs can be provided
Instruct the processor of all-purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine so that the instruction executed by computer or the processor of other programmable data processing devices is generated for real
The device for the function of being specified in present one flow of flow chart or one box of multiple flows and/or block diagram or multiple boxes.
These computer program instructions, which may also be stored in, can guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works so that instruction generation stored in the computer readable memory includes referring to
Enable the manufacture of device, the command device realize in one flow of flow chart or multiple flows and/or one box of block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device so that count
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, in computer or
The instruction executed on other programmable devices is provided for realizing in one flow of flow chart or multiple flows and/or block diagram one
The step of function of being specified in a box or multiple boxes.
Particular embodiments described above has carried out further in detail the purpose of the present invention, technical solution and advantageous effect
Describe in detail it is bright, it should be understood that the above is only a specific embodiment of the present invention, the guarantor being not intended to limit the present invention
Range is protected, all within the spirits and principles of the present invention, any modification, equivalent substitution, improvement and etc. done should be included in this
Within the protection domain of invention.
Claims (10)
1. a kind of heterogeneous reservoir Permeability Distribution inversion method, which is characterized in that including：
Determine the corresponding Lorentz curve inverting expression formula of heterogeneous reservoir difference coefficient of permeability variation；
According to the oil reservoir number of plies of heterogeneous reservoir and each layer reservoir thickness, each layer oil reservoir respective coordinates point on Lorentz curve is determined
Abscissa；
Using the corresponding Lorentz curve inverting expression formula of different coefficient of permeability variation, the corresponding vertical seat of different abscissas is determined
Mark；
According to the corresponding ordinate of different abscissas, the permeability of the respective coordinates point of each layer oil reservoir is determined.
2. the method as described in claim 1, which is characterized in that the determining heterogeneous reservoir difference coefficient of permeability variation pair
The Lorentz curve inverting expression formula answered, including：
By each layer Reservoir Permeability KiBy being ranked sequentially from small to large, K is determinediRespective layer reservoir thickness hi, obtain n groups (Ki,
hi), (x is calculated as followsm, ym), obtain the corresponding Lorentz curve of different coefficient of permeability variation：
Function inversion is carried out to the corresponding Lorentz curve of different coefficient of permeability variation as follows, obtains different permeabilities
The corresponding Lorentz curve inverting expression formula of the coefficient of variation：
Wherein, xmFor the 1st group of Lorentz curve to m group abscissas；ymFor the 1st group of Lorentz curve to m group ordinates；N is
The oil reservoir number of plies；m≤n；VKFor coefficient of permeability variation.
3. method as claimed in claim 2, which is characterized in that the determining heterogeneous reservoir difference coefficient of permeability variation pair
The Lorentz curve inverting expression formula answered, including：
Work as VKValue is in each VKWhen between value, arbitrary V is determined with linear interpolation methodKCorresponding Lorentz curve inverting expression
Formula, interpolating function are as follows：
And determine corresponding inflection point x with linear interpolation methodcValue, interpolating function is as follows：
Wherein, VkFor arbitrary VKValue；Vk0For VkPrevious VKValue；Vk1For VkLatter VKValue.
4. method as claimed in claim 2 or claim 3, which is characterized in that the oil reservoir number of plies according to heterogeneous reservoir and each layer
Reservoir thickness determines the abscissa of each layer oil reservoir respective coordinates point on Lorentz curve, including determines Lorentz as follows
The abscissa of each layer oil reservoir respective coordinates point on curve：
5. method as claimed in claim 2 or claim 3, which is characterized in that it is described according to the corresponding ordinate of different abscissas, it determines
The permeability of the respective coordinates point of each layer oil reservoir, including the permeability of the respective coordinates point of each layer oil reservoir is determined as follows：
6. a kind of heterogeneous reservoir Permeability Distribution inverting device, which is characterized in that including：
Expression formula determining module, for determining the corresponding Lorentz curve inverting table of heterogeneous reservoir difference coefficient of permeability variation
Up to formula；
Abscissa determining module is used for the oil reservoir number of plies according to heterogeneous reservoir and each layer reservoir thickness, determines Lorentz curve
The abscissa of upper each layer oil reservoir respective coordinates point；
Ordinate determining module, for using the corresponding Lorentz curve inverting expression formula of different coefficient of permeability variation, determining
The corresponding ordinate of different abscissas；
Permeability determining module, for according to the corresponding ordinate of different abscissas, determining the respective coordinates point of each layer oil reservoir
Permeability.
7. device as claimed in claim 6, which is characterized in that the expression formula determining module is further used for：
By each layer Reservoir Permeability KiBy being ranked sequentially from small to large, K is determinediRespective layer reservoir thickness hi, obtain n groups (Ki,
hi), (x is calculated as followsm, ym), obtain the corresponding Lorentz curve of different coefficient of permeability variation：
Function inversion is carried out to the corresponding Lorentz curve of different coefficient of permeability variation as follows, obtains different permeabilities
The corresponding Lorentz curve inverting expression formula of the coefficient of variation：
Wherein, xmFor the 1st group of Lorentz curve to m group abscissas；ymFor the 1st group of Lorentz curve to m group ordinates；N is
The oil reservoir number of plies；m≤n；VKFor coefficient of permeability variation.
8. device as claimed in claim 7, which is characterized in that the expression formula determining module is further used for：
Work as VKValue is in each VKWhen between value, arbitrary V is determined with linear interpolation methodKCorresponding Lorentz curve inverting expression
Formula, interpolating function are as follows：
And determine corresponding inflection point x with linear interpolation methodcValue, interpolating function is as follows：
Wherein, VkFor arbitrary VKValue；Vk0For VkPrevious VKValue；Vk1For VkLatter VKValue.
9. a kind of computer equipment, including memory, processor and storage are on a memory and the meter that can run on a processor
Calculation machine program, which is characterized in that the processor realizes any side of claim 1 to 5 when executing the computer program
Method.
10. a kind of computer readable storage medium, which is characterized in that the computer-readable recording medium storage has perform claim
It is required that the computer program of 1 to 5 any the method.
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 true CN108491361A (en) | 2018-09-04 |
| CN108491361B 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) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110242291A (en) * | 2019-07-02 | 2019-09-17 | 中国石油化工股份有限公司 | A kind of superpower heterogeneous oil and gas reservoir heterogeneity characterizing method |
| CN111456710A (en) * | 2020-04-27 | 2020-07-28 | 中国石油天然气股份有限公司 | Profile control well selection method |
| CN112347707A (en) * | 2019-08-07 | 2021-02-09 | 中国石油天然气股份有限公司 | Characterization method and device for flow heterogeneity of oil reservoir |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105525909A (en) * | 2015-11-19 | 2016-04-27 | 薛云飞 | Method for analyzing heterogeneous property of oil reservoir |
| US20160275044A1 (en) * | 2012-08-31 | 2016-09-22 | The Government Of The United States Of America, As Represented By The Secretary Of The Navy | Estimation of Hidden Variance Distribution Parameters |
| CN107219159A (en) * | 2017-04-19 | 2017-09-29 | 中国石油天然气股份有限公司 | A kind of method and apparatus for determining reservoir heterogeneity |
-
2018
- 2018-03-07 CN CN201810187730.5A patent/CN108491361B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20160275044A1 (en) * | 2012-08-31 | 2016-09-22 | The Government Of The United States Of America, As Represented By The Secretary Of The Navy | Estimation of Hidden Variance Distribution Parameters |
| 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)
| Title |
|---|
| 康晓东等: "洛伦茨曲线在油藏工程中的应用", 《新疆石油地质》 * |
| 朱小影等: "渗透率变异系数的几种计算方法", 《海洋石油》 * |
| 林伶等: "洛伦兹曲线在微观非均质性研究中的应用", 《复杂油气藏》 * |
| 王增林等: "非均质油藏开发规律研究", 《油气地质与采收率》 * |
| 邵先杰等: "储层渗透率非均质性表征新参数", 《石油实验地质》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110242291A (en) * | 2019-07-02 | 2019-09-17 | 中国石油化工股份有限公司 | A kind of superpower heterogeneous oil and gas reservoir heterogeneity characterizing method |
| CN112347707A (en) * | 2019-08-07 | 2021-02-09 | 中国石油天然气股份有限公司 | Characterization method and device for flow heterogeneity of oil reservoir |
| CN111456710A (en) * | 2020-04-27 | 2020-07-28 | 中国石油天然气股份有限公司 | Profile control well selection method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108491361B (en) | 2021-02-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108491361A (en) | heterogeneous reservoir permeability distribution inversion method and device | |
| Usuda et al. | All-electron GW calculation based on the LAPW method: Application to wurtzite ZnO | |
| Hussain et al. | Environmental Impact of Sectoral Energy Consumption on Economic Growth in Malaysia: Evidence from ARDL Bound Testing Approach. | |
| US8725478B2 (en) | Reservoir upscaling method with preserved transmissibility | |
| Zhang et al. | Preliminary discussion regarding SVM kernel function selection in the twofold rock slope prediction model | |
| CN114819696A (en) | Method and device for determining pollutant concentration influence factors and electronic equipment | |
| CN108182531A (en) | Shale gas development evaluation method, apparatus and terminal device | |
| CN110633529A (en) | Data processing method, device and system for determining oil reservoir water drive development effect | |
| Wu et al. | Local and regional flood frequency analysis based on hierarchical Bayesian model: application to annual maximum streamflow for the Huaihe River basin | |
| CN111507591A (en) | Power system state determination method, device, computer medium and storage medium | |
| CN108563889A (en) | A kind of sampled analog method of stochastic variable | |
| CN112419493B (en) | Shale reservoir three-dimensional attribute model building method and device | |
| CN113392372B (en) | Gas well full life cycle non-blocking flow determination method and system | |
| Rausch et al. | Capturing natural resource heterogeneity in top-down energy-economic equilibrium models | |
| CN110929999B (en) | Voltage sag severity calculation method considering tolerance characteristics of different devices | |
| CN106441032A (en) | Detection method of ultra-thin float glass quantification flexibility value | |
| Luo et al. | Evaluation and improvement of routing procedure for nonlinear Muskingum models | |
| CN108510778A (en) | A kind of implementation method of the vehicle polymerization display based on recorder management platform | |
| CN104574516B (en) | Point cloud smoothing system and method | |
| LÓPEZ‐BARRÓN et al. | Measurement of geometrical parameters in cocontinuous polymer blends: 3D versus 2D image analysis | |
| CN110287545A (en) | A kind of method and system of the calculating oil reservoir inside original water saturation of any part | |
| CN109614691B (en) | Random fracture network generation method and device | |
| CN103559402A (en) | Statistical analyzing method applying linear regression method | |
| CN105513984B (en) | The test method and device of the actual channel length of metal-oxide-semiconductor | |
| Dong et al. | A study on the relationship between economic production performance and pollution in chinese coastal cities: an analysis of the undesired output |
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 |