CN104991280B - A kind of perforated zone system of selection and device based on micro throat structure index - Google Patents
A kind of perforated zone system of selection and device based on micro throat structure index Download PDFInfo
- Publication number
- CN104991280B CN104991280B CN201510382782.4A CN201510382782A CN104991280B CN 104991280 B CN104991280 B CN 104991280B CN 201510382782 A CN201510382782 A CN 201510382782A CN 104991280 B CN104991280 B CN 104991280B
- Authority
- CN
- China
- Prior art keywords
- pore throat
- distributed data
- throat radius
- mrow
- magnetic resonance
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention relates to a kind of perforated zone system of selection based on micro throat structure index and device, wherein, method includes：The relational expression set up between nuclear magnetic resonance T2 distributed data and rock core sample pore throat radius distributed data；Choose the nuclear magnetic resonance T2 distributed data of target interval；Formation porosity is obtained according to the nuclear magnetic resonance T2 distributed data of target interval；The pore throat radius distributed data of target interval is obtained using the nuclear magnetic resonance T2 distributed data of target interval；The maximum pore throat radius and pore throat radius intermediate value of target interval are determined using the pore throat radius distributed data of target interval；Pore throat qualitative index is set up according to formation porosity, maximum pore throat radius and pore throat radius intermediate value；Target interval is classified according to pore throat qualitative index value, the pore throat qualitative index difference of different intervals is contrasted, the processing interval of selective value greatly is used as perforated zone.
Description
Technical field
The present invention relates to oil-gas exploration technical field, more particularly to a kind of perforated zone based on micro throat structure index
System of selection and device.
Background technology
With the increase of world energy sources demand, unconventional petroleum resources is increasingly taken seriously.Fine and close oil gas is unconventional
Compare real part in petroleum resources, be also the important making ＆ breaking produced in following oil gas increasing the storage.China's densification petroleum resources
It is widely distributed, have a high potential, be the petroleum resources that most realistic meaning and exploration and development are worth at present.Ooze in fine and close oil and gas reservoir hole
Low, anisotropism is strong, and the micro throat structure and production capacity of fine and close oil and gas reservoir are closely related, reservoir quality Quantitative Evaluation with Well Logging with
And reservoir well logging sophisticated category has become the key for improving exploration and development benefit.Therefore, compact reservoir micro throat knot is carried out
Structure is studied, and proposes rationally effective compact reservoir micro throat structure evaluation and sorting technique, for finding favorable exploration object,
Effectively carry out reservoir quality evaluation, realize that fine and close oil gas scale and benefit exploitation is significant.
Compact reservoir complex pore structure, reservoir property difference is big, is difficult to have by traditional rock physicses sorting technique
Effect ground carries out Reservoir Classification.At present generally according to the physical property characteristics such as porosity, permeability and casting body flake, pressure mercury capillary pressure
Porosity type and pore structural information that curve is obtained, such as Areal porosity, average pore size, sorting coefficient, aperture combination type, row drive
Pressure-actuatable fluid saturation etc. sets up Reservoir Classification standard, but when not having core analyzing data, then is difficult to apply the party
Method is effectively classified to reservoir.In actual applications, it usually needs Reservoir Classification is carried out using well-log information.In view of rock core
Press mercury experiment and casting body flake cycle length and may not apply to stratum and continuously evaluate, the conventional logging evaluation based on porosity
Method application effect in the classification of fine and close oil and gas reservoir is undesirable, current to lack micro- using the fine and close oil and gas reservoir of well-log information development
See pore throat character quantitative assessment and the on this basis effective ways of preferred perforated zone.Nuclear magnetic resonance log can continuously reflect storage
How layer pore-size distribution, extract the characterization parameter of reservoir micro throat information by nuclear magnetic resonance log and build the microcosmic hole of reservoir
Adam's apple structure quantitative assessment parameter is the emphasis of fine and close oil and gas reservoir quality logging evaluation, based on this, realizes fine and close oil gas storage
Layer classification and perforated zone are preferred.
The content of the invention
To solve problem of the prior art, the present invention proposes a kind of perforated zone selection based on micro throat structure index
Method and device, calculates reservoir pore throat radius using NMR logging data and is distributed, extract micro throat information, build reservoir
Micro throat qualitative index PTI, is realized to the quantitative assessment of fine and close oil and gas reservoir micro throat structure, Reservoir Classification and perforation layer
Duan Youxuan.
To achieve the above object, the invention provides a kind of perforated zone selecting party based on micro throat structure index
Method, including：
Determine the nuclear magnetic resonance T2 distributed data T of rock core sample2With pore throat radius distributed data, and set up
Nuclear magnetic resonance T2 distributed data T2Relational expression between rock core sample pore throat radius distributed data；
Choose the nuclear magnetic resonance T2 distributed data T of target interval2；
According to the nuclear magnetic resonance T2 distributed data T of target interval2Obtain formation porosity φ；
According to nuclear magnetic resonance T2 distributed data T2Pass between rock core sample pore throat radius distributed data
It is formula, utilizes the nuclear magnetic resonance T2 distributed data T of target interval2Obtain the pore throat radius distribution number of target interval
According to；
In maximum pore throat radius and pore throat radius that target interval is determined using the pore throat radius distributed data of target interval
Value；
Pore throat qualitative index PTI is set up according to formation porosity φ, maximum pore throat radius and pore throat radius intermediate value；
Target interval is classified according to pore throat qualitative index PTI values, the pore throat qualitative index PTI of different intervals is contrasted
Difference, the big processing interval of selection PTI values is used as perforated zone.
Preferably, the nuclear magnetic resonance T2 distributed data T2With rock core sample pore throat radius distributed data it
Between relational expression be：
Wherein, A, B, C, D, E determine for laboratory experiments supporting data scale；T2For nuclear magnetic resonance T2 point
Cloth data, RptFor pore throat radius distributed data.
Preferably, the expression formula of the formation porosity φ is：
Wherein, BIN is porosity component number, i.e. nuclear magnetic resonance T2The cloth points of distribution, AmiFor T2It is distributed i-th of component
Range value.
Preferably, the expression formula of the pore throat qualitative index PTI is：
PTI=ω 1*f1(Rp max)+ω2*f2(Rp50)+ω3*f3(φ)
Wherein, w1、w2、w3It is weight coefficient, is determined by analytic hierarchy process (AHP)；RpmaxMaximum pore throat radius is represented, is non-
The venturi value charged at first during wetting phase fluid-entering aperture gap network；Rp50Pore throat radius intermediate value is represented, is that pore throat radius distribution is tired
Product frequency is the corresponding pore throat radius value in 50% place；f1、f2、f3It is the normalized function of pore throat parameter；
min(Rp max) represent to study the minimum of the maximum pore throat radius distribution in area；max(Rp max) represent the maximum pore throat in research area
The maximum of radius distribution；min(Rp50) represent to study the minimum of Distribution value in area's pore throat radius；max(Rp50) represent research
The maximum of Distribution value in area's pore throat radius；Min (φ) represents the minimum of research area's porosity distribution；Max (φ) represents to grind
Study carefully the maximum of area's porosity distribution.
To achieve the above object, dress is selected present invention also offers a kind of perforated zone based on micro throat structure index
Put, including：
Relation reality unit, the nuclear magnetic resonance T2 distributed data T for determining rock core sample2And pore throat
Radius distribution data, and set up nuclear magnetic resonance T2 distributed data T2With rock core sample pore throat radius distributed data it
Between relational expression；
Nuclear magnetic resonance T2 distributed data chooses unit, and the nuclear magnetic resonance for choosing target interval laterally relaxes
Henan Annual distribution data T2；
Formation porosity acquiring unit, for the nuclear magnetic resonance T2 distributed data T according to target interval2Obtain
Obtain formation porosity φ；
Pore throat radius distributed data acquiring unit, for according to nuclear magnetic resonance T2 distributed data T2And rock core
Relational expression between sample pore throat radius distribution data, utilizes the nuclear magnetic resonance T2 distributed data T of target interval2
Obtain the pore throat radius distributed data of target interval；
Maximum pore throat radius and pore throat radius intermediate value acquiring unit, for the pore throat radius distributed data using target interval
Determine the maximum pore throat radius and pore throat radius intermediate value of target interval；
Pore throat qualitative index sets up unit, for according to formation porosity φ, maximum pore throat radius and pore throat radius intermediate value
Set up pore throat qualitative index PTI；
Perforated zone selecting unit, for being classified according to pore throat qualitative index PTI values to target interval, contrast is different
The pore throat qualitative index difference of interval, the maximum processing interval of selection pore throat qualitative index difference is used as perforated zone.
Preferably, the nuclear magnetic resonance T2 distributed data T that the relation reality unit is set up2And core sample
Relational expression between this pore throat radius distributed data is：
Wherein, A, B, C, D, E determine for laboratory experiments supporting data scale；T2For nuclear magnetic resonance T2 point
Cloth data, RptFor pore throat radius distributed data.
Preferably, the expression formula for the formation porosity φ that the formation porosity acquiring unit is obtained is：
Wherein, BIN is porosity component number, i.e. nuclear magnetic resonance T2The cloth points of distribution, AmiFor T2It is distributed i-th of component
Range value.
Preferably, the pore throat qualitative index PTI expression formula that the pore throat qualitative index sets up unit acquisition is：
PTI=ω 1*f1(Rp max)+ω2*f2(Rp50)+ω3*f3(φ)
Wherein, w1、w2、w3It is weight coefficient, is determined by analytic hierarchy process (AHP)；RpmaxMaximum pore throat radius is represented, is non-
The venturi value charged at first during wetting phase fluid-entering aperture gap network；Rp50Pore throat radius intermediate value is represented, is that pore throat radius distribution is tired
Product frequency is the corresponding pore throat radius value in 50% place；f1、f2、f3It is the normalized function of pore throat parameter；
min(Rp max) represent to study the minimum of the maximum pore throat radius distribution in area；max(Rp max) represent that research area is maximum
The maximum of pore throat radius distribution；min(Rp50) represent to study the minimum of Distribution value in area's pore throat radius；max(Rp50) represent
Study the maximum of Distribution value in area's pore throat radius；Min (φ) represents the minimum of research area's porosity distribution；Max (φ) table
Show the maximum of research area's porosity distribution.
Above-mentioned technical proposal has the advantages that：The technical program is according to nuclear magnetic resonance T2 distribution number
According to pore throat radius distribution and formation porosity is calculated, the pore throat qualitative index for characterizing reservoir micro throat structure is built, is utilized
Pore throat qualitative index realizes that the quantitative calculating of reservoir micro throat structural parameters, reservoir pore throat qualitative index are quantitatively calculated and perforation
Interval preferably, perforation formation testing is carried out to preferred interval.Have the advantages that simple, directly perceived, discrimination is high, good reliability, with bright
Aobvious practical application effect.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
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 of invention, for those of ordinary skill in the art, on the premise of not paying creative work, can be with
Other accompanying drawings are obtained according to these accompanying drawings.
Fig. 1 is a kind of perforated zone system of selection flow chart based on micro throat structure index proposed by the present invention；
Fig. 2 be using the present embodiment 1 respectively to the reservoir micro throat qualitative index result of calculation comparison diagrams of two mouthfuls of wells with
Perforated zone selects example；
Fig. 3 is a kind of perforated zone selection device frame based on micro throat structure index that the embodiment of the present invention 1 is provided
Figure；
Fig. 4 is the reservoir micro throat qualitative index computing module structural representation that the embodiment of the present invention 2 is provided.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete
Site preparation is described, it is clear that described embodiment is only a part of embodiment of the invention, rather than whole embodiments.It is based on
Embodiment in the present invention, it is every other that those of ordinary skill in the art are obtained under the premise of creative work is not made
Embodiment, belongs to the scope of protection of the invention.
As shown in figure 1, being a kind of perforated zone system of selection stream based on micro throat structure index proposed by the present invention
Cheng Tu.Including：
Step 101)：Determine the nuclear magnetic resonance T2 distributed data T of rock core sample2With pore throat radius distribution number
According to, and set up nuclear magnetic resonance T2 distributed data T2Relation between rock core sample pore throat radius distributed data
Formula；
Step 102)：Choose the nuclear magnetic resonance T2 distributed data T of target interval2；
Step 103)：According to the nuclear magnetic resonance T2 distributed data T of target interval2Obtain formation porosity
φ；
In step 103, the expression formula of the formation porosity φ is：
Wherein, BIN is porosity component number, i.e. nuclear magnetic resonance T2The cloth points of distribution, AmiFor T2It is distributed i-th of component
Range value.
Step 104)：According to nuclear magnetic resonance T2 distributed data T2With rock core sample pore throat radius distributed data
Between relational expression, utilize the nuclear magnetic resonance T2 distributed data T of target interval2Obtain the pore throat half of target interval
Footpath distributed data；
At step 104, the nuclear magnetic resonance T2 distributed data T2With the distribution of rock core sample pore throat radius
Relational expression between data is：
Wherein, A, B, C, D, E determine for laboratory experiments supporting data scale；T2For nuclear magnetic resonance T2 point
Cloth data, RptFor pore throat radius distributed data.
Step 105)：Maximum pore throat radius and the hole of target interval are determined using the pore throat radius distributed data of target interval
Larynx median radius；
Step 106)：Pore throat qualitative index is set up according to formation porosity φ, maximum pore throat radius and pore throat radius intermediate value
PTI；
In step 106, the expression formula of the pore throat qualitative index PTI is：
PTI=ω 1*f1(Rp max)+ω2*f2(Rp50)+ω3*f3(φ)
Wherein, w1、w2、w3It is weight coefficient, is determined by analytic hierarchy process (AHP)；RpmaxMaximum pore throat radius is represented, is non-
The venturi value charged at first during wetting phase fluid-entering aperture gap network；Rp50Pore throat radius intermediate value is represented, is that pore throat radius distribution is tired
Product frequency is the corresponding pore throat radius value in 50% place；f1、f2、f3It is the normalized function of pore throat parameter；
min(Rp max) represent to study the minimum of the maximum pore throat radius distribution in area；max(Rp max) represent that research area is maximum
The maximum of pore throat radius distribution；min(Rp50) represent to study the minimum of Distribution value in area's pore throat radius；max(Rp50) represent
Study the maximum of Distribution value in area's pore throat radius；Min (φ) represents the minimum of research area's porosity distribution；Max (φ) table
Show the maximum of research area's porosity distribution.
For the technical program, step 106 is specifically included：
S1：Calculate weight coefficient w1、w2、w3。
Using analytic hierarchy process (AHP), weight coefficient w is determined1、w2、w3。
Analytic hierarchy process (AHP) is also known as AHP structures power method (Analytic hierarchy process, be abbreviated as AHP), and being will be multiple
Miscellaneous evaluation object is arranged as the entirety of an orderly recursive hierarchy structure, is then carried out two-by-two between each assessment item
Comparison, judge, calculate the relatively important property coefficient of each assessment item, i.e. weight.AHP structures power method is divided into single criterion structure power again
Method and multiple criteria structure power method, introduce single criterion structure power method and specific steps herein.
(1) quantitative criteria of agriculture products
The key problem of analytic hierarchy process (AHP) is to set up a reasonable and consistent judgment matrix of construction, judgment matrix it is reasonable
Rational influence of the property by scale.So-called scale refers to estimator to each evaluation index (or project) importance rate
The quantization concept of difference.The conventional method of the quantitative criteria of agriculture products importance has：Proportion quotiety method and exponential scale method.Than
Example scaling law is based on the judgment criteria to things qualitative difference, typically to differentiate that grade represents things qualitative difference with 5 kinds.
When evaluation analysis needs higher accuracy, 9 kinds can be used to differentiate grade to evaluate, 1 is shown in Table.
Proportion quotiety value diagram of system (the importance scores x of table 1ij)
(2) initial flexible strategy are determined
The initial flexible strategy permanent method being combined frequently with qualitative analysis and quantitative analysis really.Concrete operations are：According to
The quantitative criteria of the purpose of analysis and research, the assessment indicator system having built up and the index importance primarily determined that provides each
The weight of index, calculates the average and standard deviation of each index weights, initial flexible strategy is obtained by continuing to optimize respectively.
(3) initial flexible strategy are handled
The first step, sets up judgment matrix A.By the evaluation to evaluation index, compared two-by-two, its initial flexible strategy is formed
Judgment matrix A, see the table below.The element x of ith row and jth column in judgment matrix AijRepresent index xiWith xjCompare the mark of rear gained
Spend coefficient.
The judgment matrix A of the parameter of table 2
| Index | Rpmax | Rp50 | φ |
| Rpmax | 1 | 1/5 | 1/3 |
| Rp50 | 5 | 1 | 3 |
| φ | 3 | 1/3 | 1 |
Second step, calculates the geometric mean of each calibration data of every a line in judgment matrix A, is denoted as wi。
3rd step, is normalized.Normalized is to utilize formulaCalculate, according to result of calculation
Determine the weight coefficient of each index.
Σwi=w1+w2+w3=3.8717
w1=0.1047
w2=0.6370
w3=0.2583
S2：Calculate the normalized function f of pore throat parameter1、f2、f3。
Choose the R of work area reservoir minimum and maximumpmax、Rp50, φ, i.e.,：min(Rp max)、max(Rp max), min (Rp50)、
max(Rp50), min (φ), max (φ) obtain normalized function f using following formula1、f2、f3：
Wherein：max(Rp max)=0.75, min (Rp max)=0.3, max (Rp50)=0.25, min (Rp50)=0.1, max
(φ)=12, min (φ)=4.
S3：Calculate pore throat qualitative index PTI.
According to the maximum pore throat radius Rpmax, pore throat radius intermediate value Rp50With porosity φ weight coefficient and normalizing
Change function and set up pore throat qualitative index PTI, its expression formula is：
PTI=ω 1*f1(Rp max)+ω2*f2(Rp50)+ω3*f3(φ)
Step 107)：Target interval is classified according to pore throat qualitative index PTI values, the pore throat product of different intervals are contrasted
Matter indices P TI differences, the big processing interval of selection PTI values is used as perforated zone.
As shown in figure 4, being a kind of perforated zone selection device frame based on micro throat structure index proposed by the present invention
Figure.Including：
Relation reality unit 401, the nuclear magnetic resonance T2 distributed data T for determining rock core sample2With
Pore throat radius distributed data, and set up nuclear magnetic resonance T2 distributed data T2With rock core sample pore throat radius distribution number
Relational expression between；
Nuclear magnetic resonance T2 distributed data chooses unit 402, and the nuclear magnetic resonance for choosing target interval is horizontal
To relaxation time distributed data T2；
Formation porosity acquiring unit 403, for the nuclear magnetic resonance T2 distributed data according to target interval
T2Obtain formation porosity φ；
Pore throat radius distributed data acquiring unit 404, for according to nuclear magnetic resonance T2 distributed data T2With
Relational expression between rock core sample pore throat radius distributed data, utilizes the nuclear magnetic resonance T2 distribution number of target interval
According to T2Obtain the pore throat radius distributed data of target interval；
Maximum pore throat radius and pore throat radius intermediate value acquiring unit 405, for being distributed using the pore throat radius of target interval
Data determine the maximum pore throat radius and pore throat radius intermediate value of target interval；
Pore throat qualitative index sets up unit 406, for according in formation porosity φ, maximum pore throat radius and pore throat radius
Value sets up pore throat qualitative index PTI；
Perforated zone selecting unit 407, for being classified according to pore throat qualitative index PTI values to target interval, contrast
The pore throat qualitative index difference of different intervals, the maximum processing interval of selection pore throat qualitative index difference is used as perforated zone.
For the present embodiment, the perforated zone obtained using the technical program carries out perforation formation testing to verify this technology
The superiority of scheme.
Under normal circumstances, the PTI values that selected interval is calculated are bigger, illustrate that reservoir pore throat quality is better, reservoir quality is better.
Therefore, according to PTI values, reservoir is classified (table 3), and preferably I, II class reservoir is developed as perforated zone formation testing.
Reservoir Classification table of the table 3 based on reservoir micro throat qualitative index PTI
As shown in Fig. 2 to utilize the pore throat qualitative index result figure of the different intervals of the A1 wells of the technical program acquisition.Such as
Shown in Fig. 3, to utilize the pore throat qualitative index result figure of the different intervals of the A2 wells of the technical program acquisition.Tied according to calculating
Really, No. 45 floor of A1 wells, the PTI values of No. 46 floor are respectively 0.24~0.57 (average value 0.39), 0.14~0.46 (average value
0.32), logging evaluation result is poor reservoir (IV class).No. 104 floor of A2 wells, No. 105 floor, the PTI Distribution values of No. 106 floor are 0.57
~0.87 (average value 0.73), 0.29~0.66 (average value 0.51), 0.51~1.00 (average value 0.75), logging evaluation 104
Number floor, No. 106 floor preferably reservoir (II class), illustrate No. 104 floor of A2 wells, No. 106 floor reservoir qualities preferably, therefore it is preferred that this two layers
Duan Jinhang perforation formation testings, formation testing result is 13.35 tons of oil-producing of day, and production capacity is preferable.
Realized using the inventive method and effectively classification and perforated zone are carried out to reservoir preferably, with simple, directly perceived, area
The advantages of indexing height, good reliability, with obvious practical application effect.
Above-described embodiment, has been carried out further to the purpose of the present invention, technical scheme and beneficial effect
Describe in detail, should be understood that the embodiment that the foregoing is only the present invention, be not intended to limit the present invention
Protection domain, within the spirit and principles of the invention, any modification, equivalent substitution and improvements done etc. all should be included
Within protection scope of the present invention.
Claims (6)
1. a kind of perforated zone system of selection based on micro throat structure index, it is characterised in that including：
Determine the nuclear magnetic resonance T2 distributed data T of rock core sample2With pore throat radius distributed data, and nuclear-magnetism is set up
Resonate T2 distributed data T2Relational expression between rock core sample pore throat radius distributed data；
Choose the nuclear magnetic resonance T2 distributed data T of target interval2；
Formation porosity φ is obtained according to the nuclear magnetic resonance T2 distributed data T2 of target interval；
According to nuclear magnetic resonance T2 distributed data T2Relational expression between rock core sample pore throat radius distributed data,
Utilize the nuclear magnetic resonance T2 distributed data T of target interval2Obtain the pore throat radius distributed data of target interval；
The maximum pore throat radius and pore throat radius intermediate value of target interval are determined using the pore throat radius distributed data of target interval；
Pore throat qualitative index PTI is set up according to formation porosity φ, maximum pore throat radius and pore throat radius intermediate value；
Target interval is classified according to pore throat qualitative index PTI values, the pore throat qualitative index PTI for contrasting different intervals is poor
Different, the big processing interval of selection PTI values is used as perforated zone；
Wherein, the expression formula of the pore throat qualitative index PTI is：
PTI=w1*f1(Rpmax)+w2*f2(Rp50)+w3*f3(φ)
Wherein, w1、w2、w3It is weight coefficient, is determined by analytic hierarchy process (AHP)；RpmaxMaximum pore throat radius is represented, is non-wetted
Phase fluid enters the venturi value charged at first during pore network；Rp50Pore throat radius intermediate value is represented, is pore throat radius profile accumulation frequency
Rate is the corresponding pore throat radius value in 50% place；f1、f2、f3It is the normalized function of pore throat parameter；
min(Rpmax) represent to study the minimum of the maximum pore throat radius distribution in area；max(Rpmax) represent the maximum pore throat half in research area
The maximum of footpath distribution；min(Rp50) represent to study the minimum of Distribution value in area's pore throat radius；max(Rp50) represent research area
The maximum of Distribution value in pore throat radius；Min (φ) represents the minimum of research area's porosity distribution；Max (φ) represents research
The maximum of area's porosity distribution.
2. the method as described in claim 1, it is characterised in that the nuclear magnetic resonance T2 distributed data T2And rock
Relational expression between heart sample pore throat radius distribution data is：
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>p</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>ET</mi>
<mn>2</mn>
<mi>D</mi>
</msubsup>
</mrow>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mi>A</mi>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>BT</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>C</mi>
</msup>
</mfrac>
<mo>)</mo>
</mrow>
</mfrac>
</mrow>
Wherein, A, B, C, D, E determine for laboratory experiments supporting data scale；T2For nuclear magnetic resonance T2 distribution number
According to RptFor pore throat radius distributed data.
3. the method as described in claim 1, it is characterised in that the expression formula of the formation porosity φ is：
<mrow>
<mi>&phi;</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>B</mi>
<mi>I</mi>
<mi>N</mi>
</mrow>
</munderover>
<msub>
<mi>Am</mi>
<mi>i</mi>
</msub>
</mrow>
Wherein, BIN is porosity component number, i.e. nuclear magnetic resonance T2The cloth points of distribution, AmiFor T2It is distributed the width of i-th of component
Angle value.
4. a kind of perforated zone selection device based on micro throat structure index, it is characterised in that including：
Relation reality unit, the nuclear magnetic resonance T2 distributed data T for determining rock core sample2And pore throat radius
Distributed data, and set up nuclear magnetic resonance T2 distributed data T2Between rock core sample pore throat radius distributed data
Relational expression；
When nuclear magnetic resonance T2 distributed data chooses unit, nuclear magnetic resonance transverse relaxation for choosing target interval
Between distributed data T2；
Formation porosity acquiring unit, for the nuclear magnetic resonance T2 distributed data T according to target interval2Obtain ground
Layer porosity φ；
Pore throat radius distributed data acquiring unit, for according to nuclear magnetic resonance T2 distributed data T2And rock core sample
Relational expression between pore throat radius distributed data, utilizes the nuclear magnetic resonance T2 distributed data T of target interval2Obtain
The pore throat radius distributed data of target interval；
Maximum pore throat radius and pore throat radius intermediate value acquiring unit, for being determined using the pore throat radius distributed data of target interval
The maximum pore throat radius and pore throat radius intermediate value of target interval；
Pore throat qualitative index sets up unit, for being set up according to formation porosity φ, maximum pore throat radius and pore throat radius intermediate value
Pore throat qualitative index PTI；
Perforated zone selecting unit, for classifying according to pore throat qualitative index PTI values to target interval, contrasts different intervals
Pore throat qualitative index difference, the maximum processing interval of selection pore throat qualitative index difference is used as perforated zone；
Wherein, the pore throat qualitative index PTI expression formula that the pore throat qualitative index sets up unit acquisition is：
PTI=w1*f1(Rpmax)+w2*f2(Rp50)+w3*f3(φ)
w1、w2、w3It is weight coefficient, is determined by analytic hierarchy process (AHP)；RpmaxMaximum pore throat radius is represented, is non-wetted phase fluid
The venturi value charged at first during into pore network；Rp50Pore throat radius intermediate value is represented, is that pore throat radius profile accumulation frequency is
The corresponding pore throat radius value in 50% place；f1、f2、f3It is the normalized function of pore throat parameter；
min(Rpmax) represent to study the minimum of the maximum pore throat radius distribution in area；max(Rpmax) represent the maximum pore throat half in research area
The maximum of footpath distribution；min(Rp50) represent to study the minimum of Distribution value in area's pore throat radius；max(Rp50) represent research area
The maximum of Distribution value in pore throat radius；Min (φ) represents the minimum of research area's porosity distribution；Max (φ) represents research
The maximum of area's porosity distribution.
5. device as claimed in claim 4, it is characterised in that the nuclear magnetic resonance that the relation reality unit is set up laterally relaxes
Henan Annual distribution data T2Relational expression between rock core sample pore throat radius distributed data is：
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>p</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>ET</mi>
<mn>2</mn>
<mi>D</mi>
</msubsup>
</mrow>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mi>A</mi>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>BT</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>C</mi>
</msup>
</mfrac>
<mo>)</mo>
</mrow>
</mfrac>
</mrow>
Wherein, A, B, C, D, E determine for laboratory experiments supporting data scale；T2For nuclear magnetic resonance T2 distribution number
According to RptFor pore throat radius distributed data.
6. device as claimed in claim 4, it is characterised in that the formation porosity that the formation porosity acquiring unit is obtained
φ expression formula is：
<mrow>
<mi>&phi;</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>B</mi>
<mi>I</mi>
<mi>N</mi>
</mrow>
</munderover>
<msub>
<mi>Am</mi>
<mi>i</mi>
</msub>
</mrow>
2
Wherein, BIN is porosity component number, i.e. nuclear magnetic resonance T2The cloth points of distribution, AmiFor T2It is distributed the width of i-th of component
Angle value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510382782.4A CN104991280B (en) | 2015-07-02 | 2015-07-02 | A kind of perforated zone system of selection and device based on micro throat structure index |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510382782.4A CN104991280B (en) | 2015-07-02 | 2015-07-02 | A kind of perforated zone system of selection and device based on micro throat structure index |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104991280A CN104991280A (en) | 2015-10-21 |
| CN104991280B true CN104991280B (en) | 2017-09-01 |
Family
ID=54303116
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510382782.4A Active CN104991280B (en) | 2015-07-02 | 2015-07-02 | A kind of perforated zone system of selection and device based on micro throat structure index |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104991280B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105785457B (en) * | 2016-05-13 | 2018-04-06 | 吉林大学 | The measuring method of ground nuclear magnetic resonance T2 based on bipolar pulse |
| CN106526079B (en) * | 2016-10-27 | 2019-04-16 | 中国石油大学(北京) | A method of research tight sand pore throat character dynamic change |
| CN110244369B (en) * | 2019-06-28 | 2020-10-13 | 中国石油大学(北京) | Reservoir constraint and movable fluid distribution determination method, device and system |
| CN110987751B (en) * | 2019-11-15 | 2021-05-18 | 东北石油大学 | Quantitative grading evaluation method for pore throat of compact reservoir in three-dimensional space |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102141637A (en) * | 2010-01-28 | 2011-08-03 | 中国石油天然气股份有限公司 | Method for continuously quantitative evaluation of pore structures of reservoir strata by utilizing nuclear magnetic resonance well logging data |
| CN104634718A (en) * | 2015-03-05 | 2015-05-20 | 中国石油大学(华东) | Calibration method for representing dense sandstone pore size distribution by adopting nuclear magnetic resonance |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6714009B2 (en) * | 2002-05-16 | 2004-03-30 | Schlumberger Technology Corporation | Method for the inversion of CPMG measurements enhanced by often repeated short wait time measurements |
| US20140132259A1 (en) * | 2012-11-13 | 2014-05-15 | Schlumberger Technology Corporation | Nmr method to determine grain size distribution in mixed saturation |
-
2015
- 2015-07-02 CN CN201510382782.4A patent/CN104991280B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102141637A (en) * | 2010-01-28 | 2011-08-03 | 中国石油天然气股份有限公司 | Method for continuously quantitative evaluation of pore structures of reservoir strata by utilizing nuclear magnetic resonance well logging data |
| CN104634718A (en) * | 2015-03-05 | 2015-05-20 | 中国石油大学(华东) | Calibration method for representing dense sandstone pore size distribution by adopting nuclear magnetic resonance |
Non-Patent Citations (1)
| Title |
|---|
| 利用核磁共振测井资料计算平均孔喉半径;赵文杰;《油气地质与采收率》;20090331;第16卷(第2期);第43-45页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104991280A (en) | 2015-10-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104278991B (en) | Saline Lake Facies hydrocarbon source rock organic carbon and the polynary well logging computational methods of hydrocarbon potential | |
| CN103592690B (en) | Method based on electric imaging logging porosity spectrum information automatic Reservoir Recognition crack | |
| CN104819923A (en) | Low-permeability sandstone reservoir pore structure quantitative inversion method based on nuclear magnetic resonance | |
| CN108222925A (en) | Shale gas reservoir grading integrated evaluating method | |
| CN104991280B (en) | A kind of perforated zone system of selection and device based on micro throat structure index | |
| CN104278989B (en) | Method for obtaining saturability index of low porosity and low permeability reservoir | |
| CN106526693A (en) | Crack recognition method and device | |
| CN105134195A (en) | Shale gas reservoir quality evaluation method based on logging information | |
| CN105447762B (en) | A kind of calculation method of the low-permeability oil deposit water logging information of fluid replacement | |
| Zhang et al. | Experimental study on the pore structure characteristics of tight sandstone reservoirs in Upper Triassic Ordos Basin China | |
| CN104806232B (en) | A kind of method for determining porosity lower limit of fracture | |
| CN104564042B (en) | A kind of brittle evaluation method of shale reservoir | |
| CN106050225A (en) | Method for determining 100% pure water spectrum through nuclear magnetic resonance logging spectrum T2 | |
| CN105781539A (en) | Saturability well logging calculation method of tight oil and gas reservoir | |
| CN104863574B (en) | A kind of Fluid Identification Method suitable for tight sandstone reservoir | |
| CN104564041A (en) | Low-permeability clastic rock reservoir effectiveness evaluation method based on development permeability lower limit | |
| CN103122762A (en) | Detection method and device for effective fractured intervals in unconventional shale oil and gas reservoir | |
| CN109025982A (en) | Classification evaluation method, device and the terminal device of shale gas exploitation interval | |
| CN105525917B (en) | The method for identifying fluid type in the method and identification reservoir of reservoir | |
| CN108661630A (en) | One kind being based on the preferred geology dessert quantitative evaluation method of parameter | |
| CN104343445B (en) | The method for recognizing gas-bearing formation and water layer, high gas-oil ratio (HGOR) reservoir and gas saturation | |
| CN106443800B (en) | A kind of nuclear-magnetism porosity correction method of mafic-intermediate volcanic rocks | |
| CN104712330A (en) | Well logging permeability interpretation method | |
| CN105986813A (en) | Quasi tight reservoir rapid evaluation method and quasi tight reservoir multi-index evaluation method | |
| CN106597543A (en) | Stratigraphic sedimentary facies division method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |