CN107219159A - Method and device for determining reservoir heterogeneity - Google Patents
Method and device for determining reservoir heterogeneity Download PDFInfo
- Publication number
- CN107219159A CN107219159A CN201710256713.8A CN201710256713A CN107219159A CN 107219159 A CN107219159 A CN 107219159A CN 201710256713 A CN201710256713 A CN 201710256713A CN 107219159 A CN107219159 A CN 107219159A
- Authority
- CN
- China
- Prior art keywords
- depth
- reservoir
- measured
- mrow
- permeability
- 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
- 238000000034 method Methods 0.000 title claims abstract description 38
- 230000035699 permeability Effects 0.000 claims abstract description 150
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000001764 infiltration Methods 0.000 claims description 4
- 230000008595 infiltration Effects 0.000 claims description 4
- 230000035515 penetration Effects 0.000 claims description 4
- 238000009738 saturating Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 239000011148 porous material Substances 0.000 description 4
- 239000000047 product Substances 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000010410 layer Substances 0.000 description 3
- 230000008878 coupling Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000011229 interlayer Substances 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 238000011017 operating method Methods 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 239000004568 cement Substances 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 238000000151 deposition Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The invention provides a method and a device for determining reservoir heterogeneity, wherein the method comprises the following steps: averagely grouping reservoirs to be detected according to preset depth intervals to obtain a plurality of depth groups; obtaining the average permeability value of each depth group in a plurality of depth groups; calculating to obtain a longitudinal heterogeneity index of the reservoir to be tested according to the average permeability value of each depth group; and determining the heterogeneity degree of the reservoir to be detected according to the longitudinal heterogeneity index of the reservoir to be detected. In the embodiment of the invention, the reservoir to be detected is divided into a plurality of depth groups, and the average permeability value is calculated for each depth group, so that the heterogeneity of the whole reservoir is finally determined.
Description
Technical field
The present invention relates to technical field of geological exploration, more particularly to a kind of method and apparatus for determining reservoir heterogeneity.
Background technology
Reservoir heterogeneity also refers to the various properties of reservoir with the attribute of spatial position change, can mainly show
The heterogeneous and interstitial space constituted in rock substance it is heterogeneous.For example：Clastic reservoir rock is due to depositing and the work after diagenesis
With difference, its rock forming mineral composition, matrix content, cement content are differed, so as to influence pore shape, size with
And the change of reservoir properties, and form the anisotropism of interior reservoir layer, plane and interlayer.
The existing method for determining reservoir heterogeneity, is usually the core sample for first gathering the multiple depth of reservoir, then
The permeability of each core sample is tested, recycles statistical method to determine the anisotropism of whole reservoir to be measured.That is, it is determined that
What the process of anisotropism was relied on is the core sample taken out from reservoir to be measured.
However, due in reservoir geologic structure distribution there is polytropy, reservoir is only determined with the permeability at core sample
Anisotropism easily lead to the inaccurate of final result, it is difficult to react the real anisotropism situation of reservoir to be measured.
In view of the above-mentioned problems, not yet proposing effective solution at present.
The content of the invention
The invention provides a kind of method and apparatus for determining reservoir heterogeneity, reservoir difference to be measured is portrayed deeply to reach
The purpose of anisotropism under degree.
The embodiments of the invention provide a kind of method for determining reservoir heterogeneity, it can include：By reservoir to be measured according to
Default depth interval carries out average packet, obtains multiple depth groups；Obtain oozing for each depth groups in the multiple depth groups
Saturating rate average value；According to the permeability average value of each depth groups, the vertical heterogeneity for obtaining the reservoir to be measured is calculated
Sex index；According to the vertical heterogeneity sex index of the reservoir to be measured, the anisotropism degree of the reservoir to be measured is determined.
In one embodiment, it can be calculated according to permeability average value of the below equation according to each depth groups
Obtain the vertical heterogeneity sex index of the reservoir to be measured：
Wherein, H represents the vertical heterogeneity sex index of the reservoir to be measured, ZlRepresent described l-th of depth of reservoir to be measured
The permeability average value of group, l=1,2 ..., m-1, m represents the total number of depth groups in the reservoir to be measured.
In one embodiment, reservoir to be measured is subjected to average packet according to default depth interval, obtains multiple depths
Before degree group, it can also include：According to the well-log information of the reservoir to be measured, each is obtained under the multiple depth of reservoir to be measured
The permeability of depth, wherein, include at least one depth in multiple depth groups in each depth groups；It is many to the reservoir to be measured
The permeability of each depth is standardized under individual depth.
In one embodiment, the permeability average value of each depth groups in the multiple depth groups is obtained, can be included：
Obtain the permeability of each depth under the multiple depth of current depth group；The average value of the permeability of each depth is calculated, will
The average value of the permeability of each depth as the current depth group permeability average value.
In one embodiment, the average value of the permeability of each depth can be calculated according to below equation：
Wherein, ZlThe permeability average value of described l-th of depth groups of reservoir to be measured is represented, l=1,2 ..., m, m represents institute
The total number of depth groups in reservoir to be measured is stated, t represents the depth number that the current depth group is divided, yiRepresent described to be measured
Permeability under i-th of depth of reservoir, i=t × (l-1)+1, t × (l-1)+2 ..., t × (l-1)+t.
In one embodiment, infiltration that can be according to below equation to each depth under the multiple depth of reservoir to be measured
Rate is standardized：
Wherein, yiThe standardization permeability under described i-th of depth of reservoir to be measured is represented, i=1,2 ..., n, n represents institute
State the total number of depth of reservoirs to be measured, KiRepresent the permeability under described i-th of depth of reservoir to be measured, KminRepresent described to be measured
The minimum value of the permeability of each depth, K under the multiple depth of reservoirmaxRepresent each depth under the multiple depth of reservoir to be measured
Permeability maximum.
The embodiment of the present invention additionally provides a kind of device for determining reservoir heterogeneity, can include：It is spaced grouping module,
It can be used for reservoir to be measured carrying out average packet according to default depth interval, obtain multiple depth groups；Average value obtains mould
Block, can be used for the permeability average value for obtaining each depth groups in the multiple depth groups；Anisotropism computing module, can be with
For the permeability average value according to each depth groups, the vertical heterogeneity sex index for obtaining the reservoir to be measured is calculated；
Anisotropism determining module, can be used for the vertical heterogeneity sex index according to the reservoir to be measured, determine the reservoir to be measured
Anisotropism degree.
In one embodiment, the anisotropism computing module specifically can be used for according to below equation according to described each
The permeability average value of individual depth groups, calculates the vertical heterogeneity sex index for obtaining the reservoir to be measured：
Wherein, H represents the vertical heterogeneity sex index of the reservoir to be measured, ZlRepresent described l-th of depth of reservoir to be measured
The permeability average value of group, l=1,2 ..., m-1, m represents the total number of depth groups in the reservoir to be measured.
In one embodiment, the interval grouping module can include：Permeability acquiring unit, can be used for treating
Survey reservoir and carry out average packet according to default depth interval, before obtaining multiple depth groups, according to the survey of the reservoir to be measured
Well data, obtains the permeability of each depth under the multiple depth of reservoir to be measured, wherein, each depth groups in multiple depth groups
In include at least one depth；Standardization unit, can be used for each depth under the multiple depth of reservoir to be measured
Permeability be standardized.
In one embodiment, the average value acquisition module can include：Depth of penetration rate acquiring unit, can be used for
Obtain the permeability of each depth under the multiple depth of current depth group；Permeability average calculation unit, can be used for calculating institute
The average value of the permeability of each depth is stated, the average value of the permeability of each depth is regard as the current depth group
Permeability average value.
In embodiments of the present invention, by reservoir to be measured according to default depth interval carry out average packet after, according to many
The permeability average value of each depth groups in individual depth groups, calculates the vertical heterogeneity sex index for obtaining reservoir to be measured, and determine
The anisotropism degree of reservoir to be measured.In this application, it is multiple depth groups by reservoir division to be measured, for each depth groups
The permeability average value of this depth groups is all calculated, so that the anisotropism of whole reservoir is finally determined, because it is determined that non-
It is based on stratum contents all in reservoir to be measured each depth groups during homogenieity, and not exclusively core sample,
For prior art, the anisotropism result determined is more accurate, based on more accurate anisotropism knot
Really, it is possibility to have effect improves the efficiency of oil-gas mining.
Brief description of the drawings
, below will be to embodiment or existing in order to illustrate more clearly of the embodiment of the present application or technical scheme of the prior art
There is the accompanying drawing used required in technology description to be briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments described in application, for those of ordinary skill in the art, are not paying the premise of creative labor
Under, other accompanying drawings can also be obtained according to these accompanying drawings.
Fig. 1 is a kind of method flow diagram for determination reservoir heterogeneity that the application is provided；
Fig. 2 is a kind of structured flowchart of the device for determination reservoir heterogeneity that the application is provided.
Embodiment
In order that those skilled in the art more fully understand the technical scheme in the application, it is real below in conjunction with the application
The accompanying drawing in example is applied, the technical scheme in the embodiment of the present application is clearly and completely described, it is clear that described implementation
Example only some embodiments of the present application, rather than whole embodiments.Based on the embodiment in the application, this area is common
The every other embodiment that technical staff is obtained under the premise of creative work is not made, should all belong to the application protection
Scope.
During in view of calculating reservoir heterogeneity in the prior art, due to referring to according only to core data or according only to sandstone
When number determines reservoir heterogeneity, the problem of parameters precision of resulting sign reservoir heterogeneity is relatively low, inventor proposes
Reservoir to be measured is first subjected to average packet according to default depth interval, then multiple there is default depth interval to what is obtained
Depth groups carry out anisotropism index and ask for, so that it is determined that the anisotropism degree of reservoir to be measured.Specifically, proposing one
The method for determining reservoir heterogeneity is planted, as shown in figure 1, may comprise steps of：
S101：Reservoir to be measured is subjected to average packet according to default depth interval, multiple depth groups are obtained.
Above-mentioned reservoir to be measured can be carbonate reservoir, pore type biocalcarenite reservoir.It is of course also possible to be it
His Reservoir type, the application is not construed as limiting to this.
The specific size of above-mentioned default depth interval, can be according to adjacent permeability during each depth of penetration rate of determination
Between depth difference determine, if for example, the depth difference between adjacent permeability is 0.1 meter, then above-mentioned default depth
Interval can select 0.1*n meters, i.e. default depth interval is the multiple of the depth difference between adjacent permeability.
Further, by adjusting the value of default depth interval, it can obtain under different default depth intervals
Reservoir to be measured anisotropism degree, also can be obtained by non-homogeneous degree of the reservoir to be measured under a variety of yardsticks.For example,
When default depth interval takes 0.1*n meters, the anisotropism under decimeter yardstick can be obtained, is taken in default depth interval
When 1 meter, it is possible to obtain the anisotropism under metrical scale.
In one embodiment of the application, reservoir to be measured is being subjected to average packet according to default depth interval, obtained
To before multiple depth groups, it can also include：According to the well-log information of reservoir to be measured, each is obtained under the multiple depth of reservoir to be measured
The permeability of depth, wherein, include at least one depth in multiple depth groups in each depth groups；To the multiple depths of reservoir to be measured
The permeability of each lower depth of degree is standardized.
In one embodiment, can be according to permeability of the following steps to each depth under the multiple depth of reservoir to be measured
It is standardized：
S1：Maximum, the minimum value in multiple permeabilities are obtained, and calculates the maximum of permeability and the logarithm of minimum value
Difference, obtains the permeability logarithmic difference of reservoir to be measured；
S2：The standardization permeability for obtaining each depth under the multiple depth of reservoir to be measured is calculated according to below equation：
Wherein, yiThe standardization permeability under described i-th of depth of reservoir to be measured is represented, i=1,2 ..., n, n represents institute
State the total number of depth of reservoirs to be measured, KiRepresent the permeability under described i-th of depth of reservoir to be measured, KminRepresent described to be measured
The minimum value of the permeability of each depth, K under the multiple depth of reservoirmaxRepresent each depth under the multiple depth of reservoir to be measured
Permeability maximum.
In one embodiment of the application, the multiple depths of reservoir to be measured can be obtained according to the well-log information of reservoir to be measured
The permeability of each lower depth of degree.According still further to default depth interval, such as：1m etc., it is pair corresponding with each depth multiple to ooze
Saturating rate carries out average packet, so as to obtain multiple depth groups, that is to say multiple permeability groups.I.e., it is possible to by n permeability
M permeability group is equally divided into, wherein including t permeability in each permeability group.
S102：Obtain the permeability average value of each depth groups in multiple depth groups.
Specifically, the permeability average value of current depth group can be obtained in such a way：
S2-1：Obtain the permeability of each depth under the multiple depth of current depth group.
S2-2：The average value of the permeability of each depth is calculated, using the average value of the permeability of each depth as current
The permeability average value of depth groups.
During method in S101 is got under multiple depth after the permeability of each depth, current depth group is chosen
Under corresponding multiple depth after the permeability of each depth, can by the average value of the permeability for calculating each depth,
Using the average value of the permeability of each depth as current depth group permeability average value.Specifically, can be according to following public affairs
Formula calculates the permeability average value for determining current depth group：
Wherein, ZlThe permeability average value of described l-th of depth groups of reservoir to be measured is represented, l=1,2 ..., m, m represents institute
The total number of depth groups in reservoir to be measured is stated, t represents the depth number that the current depth group is divided, yiRepresent described to be measured
Permeability under i-th of depth of reservoir, i=t × (l-1)+1, t × (l-1)+2 ..., t × (l-1)+t.
S103：According to the permeability average value of each depth groups, the vertical heterogeneity sex index for obtaining reservoir to be measured is calculated.
, can be according to below equation according to corresponding with each standardization permeability group in one embodiment of the application
Average magnitude angle value, calculates the vertical heterogeneity sex index for obtaining reservoir to be measured：
Wherein, H represents the vertical heterogeneity sex index of the reservoir to be measured, ZlRepresent described l-th of depth of reservoir to be measured
The permeability average value of group, l=1,2 ..., m-1, m represents the total number of depth groups in the reservoir to be measured.
By the permeability average value of each depth groups of reservoir to be measured resulting in S102, above-mentioned vertical heterogeneity is brought into
In the calculation formula of index, the vertical heterogeneity sex index for obtaining reservoir to be measured can be calculated.
S104：According to the vertical heterogeneity sex index of reservoir to be measured, the anisotropism degree of reservoir to be measured is determined.
Reservoir heterogeneity also refers to the various properties of reservoir with the attribute of spatial position change, can mainly show
The heterogeneous and interstitial space constituted in rock substance it is heterogeneous.Anisotropism index H is bigger, shows the non-equal of reservoir to be measured
Matter is stronger, and the physical property characteristic in terms of permeability and porosity is poorer；Anisotropism index H is smaller, shows reservoir to be measured
More homogeneous, physical property characteristic is better.
Different reservoir study objects are directed to, the depth that can be divided according to oil reservoir actual conditions selected depth group
Number, i.e. the number t of the permeability included in permeability group.By selecting suitable t, it is possible to achieve different depth interval is treated
The anisotropism for surveying reservoir is characterized.
The distribution of vertical heterogeneity sex index is 0~1.Specifically, table 1 is given according to vertical heterogeneity sex index
Judge the standard of reservoir heterogeneity degree.
Relation between the anisotropism index of table 1 and anisotropism degree
In one embodiment of the application, above-mentioned S101 calculates the side of the anisotropism degree of reservoir to be measured into S105
Method is applied to individual well., can be by calculating the institute in interval to be measured when calculating the vertical heterogeneity sex index in interval to be measured
There is the heterogeneous index H of individual well, the longitudinal direction determined in interval to be measured is integrated according to resulting multiple heterogeneous exponential quantities non-
Matter sex index.
With reference to a specific embodiment to the corresponding mark of each depth under the multiple depth of above-mentioned acquisition reservoir to be measured
The method of standardization permeability is specifically described, it is important to note, however, that the specific embodiment is merely to be better described
The present invention, does not constitute inappropriate limitation of the present invention.
Somewhere, which is one, has 5 intervals, and each interval thickness in more than 10m.This area is pore type
Biocalcarenite reservoir, average pore is 19.3%, and mean permeability is 38mD, and crack agensis, reservoir inside is developed
Every interlayer and hypertonic band, the anisotropism of reservoir is very strong.
For above-mentioned zone, the anisotropism degree of the individual well of certain in this area can be obtained in such a way：
To somewhere each standardization permeability be grouped when, it is contemplated that each well section thickness 10m with
On, meanwhile, the difference between the corresponding depth of adjacent permeability is 0.1m, therefore, it can the well section of this area being equally divided into
10 depth groups, each depth groups include 10 depth, that is, obtain 100 permeabilities, and are classified as 10 permeability groups, often
There are 10 permeabilities in group.
Obtain a series of with the discrete data point that depth (0.1m) is variable according to interpretation of logging data result, wherein, obtain
A series of discrete data points using depth as variable taken can be：X1, X2, X3…XnEtc. a series of depth points, and, it is and upper
State a series of corresponding permeability in depth points：K1, K2, K3…Kn.The characteristics of in view of permeability, permeability can be carried out
Standardization, obtains standardization permeability y1, y2, y3……yn.Specifically, can be expressed as follows：
Wherein, yiThe standardization permeability under described i-th of depth of reservoir to be measured is represented, i=1,2 ..., n, n represents institute
State the total number of depth of reservoirs to be measured, KiRepresent the permeability under described i-th of depth of reservoir to be measured, KminRepresent described to be measured
The minimum value of the permeability of each depth, K under the multiple depth of reservoirmaxRepresent each depth under the multiple depth of reservoir to be measured
Permeability maximum.
Understand：yi∈[0,1].That is, the standardization permeability range of reservoir to be measured is between zero and one.
Respectively into permeability mean value calculation formula is brought to the corresponding permeability of 10 depth in 10 depth groups, can be with
10 permeability average value corresponding with each depth groups is obtained, and is obtained according to the permeability mean value calculation of 10 depth groups
The vertical heterogeneity sex index of somewhere individual well well section, according to vertical heterogeneity sex index H, determines somewhere individual well well
The anisotropism degree of section.
Further, the vertical heterogeneity sex index H of other individual wells of somewhere can be repeatedly calculated in the manner described above, and
The multiple vertical heterogeneity sex index H obtained according to calculating distribution, and table 1 is combined, each well of somewhere can be obtained
The anisotropism degree of section, basis and reference are provided for later stage oil field development.
Based on same inventive concept, a kind of device for determining reservoir heterogeneity is additionally provided in the embodiment of the present invention, such as
Described in the following examples.Due to determining that the device of reservoir heterogeneity solves the principle of problem with determining reservoir heterogeneity
Method is similar, it is thus determined that the implementation of the device of reservoir heterogeneity may refer to determine the reality of the method for reservoir heterogeneity
Apply, repeat part and repeat no more.Used below, term " unit " or " module " can realize the software of predetermined function
And/or the combination of hardware.Although the device described by following examples is preferably realized with software, hardware, or soft
The realization of the combination of part and hardware is also that may and be contemplated.Fig. 2 is the determination reservoir heterogeneity of the embodiment of the present invention
A kind of structured flowchart of device, as shown in Fig. 2 can include：Be spaced grouping module 201, average value acquisition module 202, it is non-
Matter computing module 203, anisotropism determining module 204, are illustrated to the structure below.
Grouping module 201 is spaced, can be used for reservoir to be measured carrying out average packet according to default depth interval, obtain
Multiple depth groups；
Average value acquisition module 202, the permeability that can be used for obtaining each depth groups in the multiple depth groups is averaged
Value；
Anisotropism computing module 203, can be used for the permeability average value according to each depth groups, and calculating is obtained
The vertical heterogeneity sex index of the reservoir to be measured；
Anisotropism determining module 204, can be used for the vertical heterogeneity sex index according to the reservoir to be measured, determines institute
State the anisotropism degree of reservoir to be measured.
In one embodiment, the anisotropism computing module specifically can be used for according to below equation according to described each
The permeability average value of individual depth groups, calculates the vertical heterogeneity sex index for obtaining the reservoir to be measured：
Wherein, H represents the vertical heterogeneity sex index of the reservoir to be measured, ZlRepresent described l-th of depth of reservoir to be measured
The permeability average value of group, l=1,2 ..., m-1, m represents the total number of depth groups in the reservoir to be measured.
In one embodiment, the interval grouping module can include：Permeability acquiring unit, for by storage to be measured
Layer carries out average packet according to default depth interval, obtains before multiple depth groups, is provided according to the well logging of the reservoir to be measured
Material, obtains the permeability of each depth under the multiple depth of reservoir to be measured, wherein, wrapped in multiple depth groups in each depth groups
Contain at least one depth；Standardization unit, for the permeability to each depth under the multiple depth of reservoir to be measured
It is standardized.
In one embodiment, the average value acquisition module can include：Depth of penetration rate acquiring unit, for obtaining
The permeability of each depth under the multiple depth of current depth group；Permeability average calculation unit, for calculating each described depth
The average value of the permeability of degree, puts down the average value of the permeability of each depth as the permeability of the current depth group
Average.
In one embodiment, the permeability average calculation unit specifically can be used for calculating institute according to below equation
State the average value of the permeability of each depth：
Wherein, ZlThe permeability average value of described l-th of depth groups of reservoir to be measured is represented, l=1,2 ..., m, m represents institute
The total number of depth groups in reservoir to be measured is stated, t represents the depth number that the current depth group is divided, yiRepresent described to be measured
Permeability under i-th of depth of reservoir, i=t × (l-1)+1, t × (l-1)+2 ..., t × (l-1)+t.
In one embodiment, the standardization unit specifically can be used for according to below equation to the storage to be measured
The permeability of each depth is standardized under the multiple depth of layer：
Wherein, yiThe standardization permeability under described i-th of depth of reservoir to be measured is represented, i=1,2 ..., n, n represents institute
State the total number of depth of reservoirs to be measured, KiRepresent the permeability under described i-th of depth of reservoir to be measured, KminRepresent described to be measured
The minimum value of the permeability of each depth, K under the multiple depth of reservoirmaxRepresent each depth under the multiple depth of reservoir to be measured
Permeability maximum.
As can be seen from the above description, the embodiment of the present invention realizes following technique effect：By reservoir to be measured according to
Default depth interval is carried out after average packet, according to the permeability average value of each depth groups in multiple depth groups, is calculated
The vertical heterogeneity sex index of reservoir to be measured is obtained, and determines the anisotropism degree of reservoir to be measured.In this application, will be to be measured
Reservoir division is multiple depth groups, the permeability average value of this depth groups is calculated for each depth groups, so that finally
The anisotropism of whole reservoir is determined, because it is determined that being based on reservoir to be measured each depth groups during anisotropism
In all stratum content, and not exclusively core sample, for prior art, the anisotropism determined
As a result it is more accurate, based on more accurate anisotropism result, it is possibility to have effect improves the efficiency of oil-gas mining.
Although mentioning the average value-acquiring method of permeability, vertical heterogeneity index calculation method, infiltration in teachings herein
Rate standardization processing method etc. is described, and still, the application is not limited to be the situation described by the embodiment of the present application.Certain
Embodiment amended slightly also may be used in a little professional standards or the practice processes described using self-defined mode or embodiment
With realize above-described embodiment it is identical, equivalent or close or deformation after it is anticipated that implementation result.After these modifications or deformation
The acquisition such as permeability average value-acquiring method, vertical heterogeneity index calculation method, permeability standardization processing method reality
Example is applied, within the scope of the optional embodiment that still may belong to the application.
Although this application provides the method operating procedure as described in embodiment or flow chart, based on conventional or noninvasive
The means for the property made can include more or less operating procedures.The step of being enumerated in embodiment order is only numerous steps
A kind of mode in execution sequence, unique execution sequence is not represented., can be with when device in practice or end product execution
Performed according to embodiment or method shown in the drawings order or parallel execution (such as parallel processor or multiple threads
Environment, even distributed data processing environment).Term " comprising ", "comprising" or its any other variant are intended to
Nonexcludability is included, so that process, method, product or equipment including a series of key elements not only will including those
Element, but also other key elements including being not expressly set out, or also include being this process, method, product or equipment
Intrinsic key element.In the absence of more restrictions, be not precluded from the process including the key element, method, product or
Also there are other identical or equivalent elements in person's equipment.
Unit, device or module that above-described embodiment is illustrated etc., can specifically be realized by computer chip or entity, or
Realized by the product with certain function.For convenience of description, describe to be divided into various modules point during apparatus above with function
Do not describe.Certainly, when implementing the application can the function of each module in same or multiple softwares and/or hardware it is real
It is existing, the module for realizing same function can also be realized by the combination of multiple submodule or subelement etc..Dress described above
It is only schematical to put embodiment, for example, the division of the unit, only a kind of division of logic function, when actually realizing
There can be other dividing mode, such as multiple units or component can combine or be desirably integrated into another system, or one
A little features can be ignored, or not perform.It is another, shown or discussed coupling or direct-coupling or communication link each other
Connect can be can be electrical, machinery or other shapes by some interfaces, the INDIRECT COUPLING or communication connection of device or unit
Formula.
It is also known in the art that in addition to realizing controller in pure computer readable program code mode, it is complete
Controller can be caused with gate, switch, application specific integrated circuit, programmable by the way that method and step is carried out into programming in logic entirely
Logic controller realizes identical function with the form of embedded microcontroller etc..Therefore this controller is considered one kind
Hardware component, and the device for realizing various functions included to its inside can also be considered as the structure in hardware component.Or
Person even, not only can be able to will be the software module of implementation method but also can be hardware for realizing that the device of various functions is considered as
Structure in part.
The application can be described in the general context of computer executable instructions, such as program
Module.Usually, program module includes performing particular task or realizes routine, program, object, the group of particular abstract data type
Part, data structure, class etc..The application can also be put into practice in a distributed computing environment, in these DCEs,
Task is performed by the remote processing devices connected by communication network.In a distributed computing environment, program module can
With positioned at including in the local and remote computer-readable storage medium including storage device.
As seen through the above description of the embodiments, those skilled in the art can be understood that the application can
Realized by the mode of software plus required general hardware platform.Understood based on such, the technical scheme essence of the application
On the part that is contributed in other words to prior art can be embodied in the form of software product, the computer software product
It can be stored in storage medium, such as ROM/RAM, magnetic disc, CD, including some instructions are to cause a computer equipment
(can be personal computer, mobile terminal, server, or network equipment etc.) performs each embodiment of the application or implementation
Method described in some parts of example.
Each embodiment in this specification is described by the way of progressive, same or analogous portion between each embodiment
Divide mutually referring to what each embodiment was stressed is the difference with other embodiment.The application can be used for crowd
In more general or special purpose computing system environments or configuration.For example：Personal computer, server computer, handheld device or
Portable set, laptop device, multicomputer system, the system based on microprocessor, set top box, programmable electronics are set
Standby, network PC, minicom, DCE of mainframe computer including any of the above system or equipment etc..
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application have it is many deformation and
Change is without departing from spirit herein, it is desirable to which appended claim includes these deformations and changed without departing from the application's
Spirit.
Claims (10)
1. a kind of method for determining reservoir heterogeneity, it is characterised in that including：
Reservoir to be measured is subjected to average packet according to default depth interval, multiple depth groups are obtained；
Obtain the permeability average value of each depth groups in the multiple depth groups；
According to the permeability average value of each depth groups, the vertical heterogeneity sex index for obtaining the reservoir to be measured is calculated；
According to the vertical heterogeneity sex index of the reservoir to be measured, the anisotropism degree of the reservoir to be measured is determined.
2. the method as described in claim 1, it is characterised in that according to permeability of the below equation according to each depth groups
Average value, calculates the vertical heterogeneity sex index for obtaining the reservoir to be measured：
<mrow>
<mi>H</mi>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<mo>|</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>Z</mi>
<mi>l</mi>
</msub>
<mo>|</mo>
</mrow>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
</mrow>
Wherein, H represents the vertical heterogeneity sex index of the reservoir to be measured, ZlRepresent oozing for described l-th of depth groups of reservoir to be measured
Saturating rate average value, l=1,2 ..., m-1, m represents the total number of depth groups in the reservoir to be measured.
3. the method as described in claim 1, it is characterised in that be averaged by reservoir to be measured according to default depth interval
Packet, is obtained before multiple depth groups, methods described also includes：
According to the well-log information of the reservoir to be measured, the permeability of each depth under the multiple depth of reservoir to be measured is obtained, its
In, include at least one depth in multiple depth groups in each depth groups；
Permeability to each depth under the multiple depth of reservoir to be measured is standardized.
4. method as claimed in claim 3, it is characterised in that obtain the permeability of each depth groups in the multiple depth groups
Average value, including：
Obtain the permeability of each depth under the multiple depth of current depth group；
The average value of the permeability of each depth is calculated, the average value of the permeability of each depth is worked as described
The permeability average value of preceding depth groups.
5. method as claimed in claim 4, it is characterised in that the permeability of each depth is calculated according to below equation
Average value：
<mrow>
<msub>
<mi>Z</mi>
<mi>l</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>t</mi>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>t</mi>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>t</mi>
</mrow>
</munderover>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
</mrow>
<mi>t</mi>
</mfrac>
</mrow>
Wherein, ZlThe permeability average value of described l-th of depth groups of reservoir to be measured is represented, l=1,2 ..., m, m represents described to be measured
The total number of depth groups in reservoir, t represents the depth number that the current depth group is divided, yiRepresent the reservoir to be measured
Permeability under i depth, i=t × (l-1)+1, t × (l-1)+2 ..., t × (l-1)+t.
6. method as claimed in claim 3, it is characterised in that according to below equation to each under the multiple depth of reservoir to be measured
The permeability of individual depth is standardized：
<mrow>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>lg</mi>
<mi> </mi>
<msub>
<mi>K</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mi>lg</mi>
<mi> </mi>
<msub>
<mi>K</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mi>lg</mi>
<mi> </mi>
<msub>
<mi>K</mi>
<mi>max</mi>
</msub>
<mo>-</mo>
<mi>lg</mi>
<mi> </mi>
<msub>
<mi>K</mi>
<mi>min</mi>
</msub>
</mrow>
</mfrac>
</mrow>
Wherein, yiThe standardization permeability under described i-th of depth of reservoir to be measured is represented, i=1,2 ..., n, n represents described to be measured
The total number of depth of reservoirs, KiRepresent the permeability under described i-th of depth of reservoir to be measured, KminRepresent that the reservoir to be measured is more
The minimum value of the permeability of each depth, K under individual depthmaxRepresent the infiltration of each depth under the multiple depth of reservoir to be measured
The maximum of rate.
7. a kind of device for determining reservoir heterogeneity, it is characterised in that including：
Grouping module is spaced, for reservoir to be measured to be carried out into average packet according to default depth interval, multiple depth groups are obtained；
Average value acquisition module, the permeability average value for obtaining each depth groups in the multiple depth groups；
Anisotropism computing module, for the permeability average value according to each depth groups, calculating obtains the storage to be measured
The vertical heterogeneity sex index of layer；
Anisotropism determining module, for the vertical heterogeneity sex index according to the reservoir to be measured, determines the reservoir to be measured
Anisotropism degree.
8. device as claimed in claim 7, it is characterised in that the anisotropism computing module is specifically for according to following public affairs
Formula calculates the vertical heterogeneity sex index for obtaining the reservoir to be measured according to the permeability average value of each depth groups：
<mrow>
<mi>H</mi>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<mo>|</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>Z</mi>
<mi>l</mi>
</msub>
<mo>|</mo>
</mrow>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
</mrow>
Wherein, H represents the vertical heterogeneity sex index of the reservoir to be measured, ZlRepresent oozing for described l-th of depth groups of reservoir to be measured
Saturating rate average value, l=1,2 ..., m-1, m represents the total number of depth groups in the reservoir to be measured.
9. device as claimed in claim 7, it is characterised in that the interval grouping module includes：
Permeability acquiring unit, for reservoir to be measured to be carried out into average packet according to default depth interval, obtains multiple depths
Before degree group, according to the well-log information of the reservoir to be measured, the infiltration of each depth under the multiple depth of reservoir to be measured is obtained
Rate, wherein, include at least one depth in multiple depth groups in each depth groups；
Standardization unit, place is standardized for the permeability to each depth under the multiple depth of reservoir to be measured
Reason.
10. device as claimed in claim 9, it is characterised in that the average value acquisition module includes：
Depth of penetration rate acquiring unit, the permeability for obtaining each depth under the multiple depth of current depth group；
Permeability average calculation unit, the average value of the permeability for calculating each depth will each described depth
Permeability average value as the current depth group permeability average value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710256713.8A CN107219159B (en) | 2017-04-19 | 2017-04-19 | Method and device for determining reservoir heterogeneity |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710256713.8A CN107219159B (en) | 2017-04-19 | 2017-04-19 | Method and device for determining reservoir heterogeneity |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107219159A true CN107219159A (en) | 2017-09-29 |
| CN107219159B CN107219159B (en) | 2019-07-09 |
Family
ID=59927659
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710256713.8A Active CN107219159B (en) | 2017-04-19 | 2017-04-19 | Method and device for determining reservoir heterogeneity |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107219159B (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108414415A (en) * | 2018-01-27 | 2018-08-17 | 东北石油大学 | A kind of device and manufacturing method using Centimeter Level rock core simulated formation anisotropism |
| CN108491361A (en) * | 2018-03-07 | 2018-09-04 | 中国石油大学(北京) | heterogeneous reservoir permeability distribution inversion method and device |
| CN109113732A (en) * | 2018-08-09 | 2019-01-01 | 中国石油天然气股份有限公司 | Reservoir heterogeneity determination method and device |
| CN110320278A (en) * | 2018-03-30 | 2019-10-11 | 中国石油化工股份有限公司 | Full-hole core heterogeneity measurement device and method |
| CN110568501A (en) * | 2018-06-05 | 2019-12-13 | 中国石油化工股份有限公司 | Quantitative characterization method for longitudinal heterogeneous characteristics of hydrocarbon source rock in stratum |
| CN110966000A (en) * | 2020-01-02 | 2020-04-07 | 中国石油集团测井有限公司华北分公司 | Method and system for representing glutenite reservoir stratiform index |
| CN113107463A (en) * | 2021-04-22 | 2021-07-13 | 中国石油天然气集团有限公司 | Reservoir heterogeneity evaluation method |
| CN114089421A (en) * | 2021-12-06 | 2022-02-25 | 中国矿业大学 | Oil and gas reservoir heterogeneity analysis method |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2015107499A1 (en) * | 2014-01-20 | 2015-07-23 | National Centre For Cell Science | A tumor deconstruction platform for the analysis of intra-tumor heterogeneity |
| WO2015153506A1 (en) * | 2014-03-31 | 2015-10-08 | Ingrain, Inc. | Digital rock physics-based trend determination and usage for upscaling |
| CN105092820A (en) * | 2015-08-06 | 2015-11-25 | 中国石油集团长城钻探工程有限公司地质研究院 | Bioclastic limestone reservoir heterogeneity research method |
| CN105092446A (en) * | 2014-05-08 | 2015-11-25 | 中国海洋石油总公司 | Two-dimensional visual percolation experiment method capable of simulating in-layer heterogeneity |
| CN105651676A (en) * | 2016-02-17 | 2016-06-08 | 中国石油大学(华东) | Reservoir heterogeneity characterization method under regular development well pattern of horizontal wells |
| CN106355571A (en) * | 2016-08-15 | 2017-01-25 | 中国石油天然气股份有限公司 | Method and device for determining quality of dolomite reservoir |
-
2017
- 2017-04-19 CN CN201710256713.8A patent/CN107219159B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2015107499A1 (en) * | 2014-01-20 | 2015-07-23 | National Centre For Cell Science | A tumor deconstruction platform for the analysis of intra-tumor heterogeneity |
| WO2015153506A1 (en) * | 2014-03-31 | 2015-10-08 | Ingrain, Inc. | Digital rock physics-based trend determination and usage for upscaling |
| CN105092446A (en) * | 2014-05-08 | 2015-11-25 | 中国海洋石油总公司 | Two-dimensional visual percolation experiment method capable of simulating in-layer heterogeneity |
| CN105092820A (en) * | 2015-08-06 | 2015-11-25 | 中国石油集团长城钻探工程有限公司地质研究院 | Bioclastic limestone reservoir heterogeneity research method |
| CN105651676A (en) * | 2016-02-17 | 2016-06-08 | 中国石油大学(华东) | Reservoir heterogeneity characterization method under regular development well pattern of horizontal wells |
| CN106355571A (en) * | 2016-08-15 | 2017-01-25 | 中国石油天然气股份有限公司 | Method and device for determining quality of dolomite reservoir |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108414415A (en) * | 2018-01-27 | 2018-08-17 | 东北石油大学 | A kind of device and manufacturing method using Centimeter Level rock core simulated formation anisotropism |
| CN108414415B (en) * | 2018-01-27 | 2020-04-28 | 东北石油大学 | Device for simulating formation heterogeneity by centimeter-level core and manufacturing method |
| CN108491361A (en) * | 2018-03-07 | 2018-09-04 | 中国石油大学(北京) | heterogeneous reservoir permeability distribution inversion method and device |
| CN108491361B (en) * | 2018-03-07 | 2021-02-23 | 中国石油大学(北京) | Heterogeneous reservoir permeability distribution inversion method and device |
| CN110320278A (en) * | 2018-03-30 | 2019-10-11 | 中国石油化工股份有限公司 | Full-hole core heterogeneity measurement device and method |
| CN110320278B (en) * | 2018-03-30 | 2021-09-14 | 中国石油化工股份有限公司 | Device and method for measuring heterogeneity of full-diameter core |
| CN110568501A (en) * | 2018-06-05 | 2019-12-13 | 中国石油化工股份有限公司 | Quantitative characterization method for longitudinal heterogeneous characteristics of hydrocarbon source rock in stratum |
| CN109113732A (en) * | 2018-08-09 | 2019-01-01 | 中国石油天然气股份有限公司 | Reservoir heterogeneity determination method and device |
| CN110966000A (en) * | 2020-01-02 | 2020-04-07 | 中国石油集团测井有限公司华北分公司 | Method and system for representing glutenite reservoir stratiform index |
| CN113107463A (en) * | 2021-04-22 | 2021-07-13 | 中国石油天然气集团有限公司 | Reservoir heterogeneity evaluation method |
| CN114089421A (en) * | 2021-12-06 | 2022-02-25 | 中国矿业大学 | Oil and gas reservoir heterogeneity analysis method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107219159B (en) | 2019-07-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107219159A (en) | Method and device for determining reservoir heterogeneity | |
| CN106526693A (en) | crack identification method and device | |
| CN104712330B (en) | One kind well logging permeability means of interpretation | |
| CN107402176B (en) | method and device for determining porosity of crack | |
| CN105445441B (en) | Reservoir gas saturation determination method and device | |
| Santos et al. | Characterization of natural fracture systems: Analysis of uncertainty effects in linear scanline results | |
| CN105388525A (en) | Reservoir prediction method and device | |
| CN104914470A (en) | Carbonate rock fracture-cavity reservoir reserve correction method | |
| CN107366534A (en) | Method and device for determining coarsening permeability | |
| US11008835B2 (en) | Horizontal reservoir description systems | |
| CN106546525A (en) | The method and apparatus for setting up three-dimensional penetration rate model | |
| CN106355571A (en) | Method and device for determining quality of dolomite reservoir | |
| CN106323835B (en) | Method for determining cementing index of heterogeneous carbonate reservoir | |
| Li et al. | A novel approach of two-dimensional representation of rock fracture network characterization and connectivity analysis | |
| CN105257284A (en) | Method and device for determining content of coagulated ash by utilizing element capture energy spectrum logging | |
| CN104374679A (en) | Method for determining permeability contribution value of reservoir | |
| CN110232208A (en) | The Reservoir Parameter Models construction method and equipment of advantage seepage flow section for identification | |
| US11249218B2 (en) | Horizontal reservoir description systems | |
| CN107701177A (en) | Geological exploration method and device | |
| CN105844011A (en) | Capillary model-based permeability computing method | |
| CN107780926A (en) | Borehole section shape measuring method, device and terminal | |
| CN103901475B (en) | Attribute contour map drawing method and device | |
| CN105240003A (en) | Ideal formation factor determination method, hydrocarbon saturation determination method and porosity effectiveness evaluation method | |
| CN111155980B (en) | Water flow dominant channel identification method and device | |
| CN106934725A (en) | Rock formation median radius forecast model method for building up, apparatus and system |
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 |