CN108226002A - The method that secondary porosity is obtained based on scale Electrical imaging secondary porosity spectrum cutoff value - Google Patents
The method that secondary porosity is obtained based on scale Electrical imaging secondary porosity spectrum cutoff value Download PDFInfo
- Publication number
- CN108226002A CN108226002A CN201711271137.0A CN201711271137A CN108226002A CN 108226002 A CN108226002 A CN 108226002A CN 201711271137 A CN201711271137 A CN 201711271137A CN 108226002 A CN108226002 A CN 108226002A
- Authority
- CN
- China
- Prior art keywords
- porosity
- spectrum
- secondary porosity
- cutoff value
- rock core
- 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
- 238000001228 spectrum Methods 0.000 title claims abstract description 142
- 238000003384 imaging method Methods 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 35
- 239000011435 rock Substances 0.000 claims abstract description 56
- 238000004364 calculation method Methods 0.000 claims abstract description 47
- 238000005481 NMR spectroscopy Methods 0.000 claims abstract description 30
- 239000011148 porous material Substances 0.000 claims description 55
- 238000005266 casting Methods 0.000 claims description 45
- 238000005119 centrifugation Methods 0.000 claims description 24
- 238000004458 analytical method Methods 0.000 claims description 9
- 230000003595 spectral effect Effects 0.000 claims description 9
- 238000012795 verification Methods 0.000 claims description 8
- 239000012530 fluid Substances 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 4
- 238000005260 corrosion Methods 0.000 claims description 4
- 230000007797 corrosion Effects 0.000 claims description 4
- 238000004090 dissolution Methods 0.000 claims description 4
- 238000012797 qualification Methods 0.000 claims description 3
- 238000013459 approach Methods 0.000 abstract description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 10
- 238000009825 accumulation Methods 0.000 description 5
- 238000000205 computational method Methods 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 235000019738 Limestone Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002902 bimodal effect Effects 0.000 description 1
- 230000004071 biological effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research 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
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10088—Magnetic resonance imaging [MRI]
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Biochemistry (AREA)
- Dispersion Chemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Pathology (AREA)
- Geometry (AREA)
- Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses it is a kind of based on scale Electrical imaging secondary porosity spectrum cutoff value obtain secondary porosity method, including：S1, secondary porosity spectrum cutoff value calculation formula is established：Secondary porosity composes cutoff value=a × δ+porosity intermediate value, and a is fixed coefficient, and δ is the porosity variance in selected window in minimal amount of porosity and porosity median value range on electric imaging logging image；S2, rock core secondary porosity is calculated using rock core nuclear magnetic resonance T 2 spectrum；S3, the numerical value using fixed coefficient a in the rock core secondary porosity scale secondary porosity spectrum cutoff value calculation formula of rock core nuclear magnetic resonance T 2 spectrum calculating gained, to obtain secondary porosity after scale.The present invention uses rock core nuclear magnetic resonance T 2 spectrum scale, the fixed coefficient a in secondary porosity spectrum cutoff value calculation formula is obtained using the method for Approach by inchmeal, the reliability of secondary porosity spectrum cutoff value result of calculation is increased, and then obtains the high secondary porosity acquisition methods of accuracy.
Description
Technical field
The present invention relates to carbonate reservoir secondary porosity acquiring technology field more particularly to it is a kind of based on scale electricity into
The method for obtaining secondary porosity as secondary porosity spectrum cutoff value.
Background technology
Carbonate reservoir has matrix pores (primary pore) and secondary pore is developed simultaneously, pore throat character is complicated, storage
Collect the reservoir characteristics such as spatial diversity, reservoir heterogeneity be strong.Due to influences such as diagenesis, primary pore is often difficult to effectively protect
It deposits, and the aperture of secondary pore is usually bigger than the aperture of primary pore, permeability is also better than primary pore, is carbonate reservoir
The major reservoir spaces and transport channel of middle oil gas, therefore carry out carbonate reservoir secondary porosity Logging Evaluation Method and grind
Study carefully, especially accurately calculate secondary porosity, have great importance for carbonate reservoir efficiency evaluation, also can be
Oilfield development program is formulated and is adjusted, effectively improves recovery ratio offer theoretical foundation.
Since the formation of secondary pore is by series of factors such as corrosion, cementation, diagenesis, biological effects
Control, the response relation of single log parameter and secondary porosity is generally poor, can not effectively be commented using Conventional Logs
Valency secondary porosity, therefore secondary porosity logging evaluation research is carried out using electric imaging logging data both at home and abroad at present,
The middle method for calculating secondary porosity is broadly divided into following 3 steps：
(1) it is composed using electric imaging logging material computation reservoir porosity：An electric imaging logging image window is chosen first,
Borehole wall flushed zone conductivity value in window is converted into borehole wall flushed zone resistivity value, is rinsed the borehole wall using Archie formula
Strip resistance rate matrix is converted to the porosity size of each imaging logging pixel, then the hole to reservoir in a certain depth bounds
Angle value carries out Frequency statistics in different porosity value sections, and it is straight that the result after statistics is depicted as frequency in coordinate system
Fang Tu, the Electrical imaging porosity spectrum on the frequency histogram, that is, stratum.
(2) cutoff value that primary pore and secondary pore are distinguished in porosity spectrum is calculated：Schlumberger
B.M.Newberry etc. 1996 propose based on porosity intermediate value Newberry porositys spectrum cut-off value calculating method and
The discriminate porosity based on normal linearity discriminant analysis that T.S.Ramakrishnan etc. was proposed in 1997 composes cutoff value
Computational methods are current mainstream computational methods, both porositys spectrum cut-off value calculating method is the mathematics based on porosity
Computational methods can either reflect reservoir pore space feature, and not influenced by porosity spectrum morphological feature, but porosity spectrum is cut
Only value computational accuracy depends on the value of regional experience coefficient, and computational methods applicability and result of calculation reliability be not high, directly
Influence the computational accuracy of secondary porosity in next step.Domestic Shi Feizhou et al. is in order to solve these problems, it is proposed that
Gaussian function fitting porosity spectrum cutoff value method computational methods based on porosity spectral shape, but the computational accuracy of this method and hole
Porosity spectral shape is in close relations, is not suitable for the situation that bimodal spacing is larger in porosity spectrum, answering in carbonate reservoir
With being limited to.
(3) the reservoir secondary porosity based on secondary porosity spectrum cutoff value calculates：It is calculated according to porosity spectrum cutoff value
As a result the line of demarcation of primary pore and secondary pore is determined in porosity spectrum, to the secondary pore spectrum envelope on rear side of the line of demarcation
The area that line is surrounded is integrated, and integral and calculating result is secondary porosity.
According to the above method and steps for calculating carbonate reservoir secondary porosity, by electric imaging logging data, have
Effect solves the problems, such as conventional logging can not be used to calculate carbonate rock secondary porosity.But link is calculated in secondary porosity,
Or computational accuracy depends on regional experience parameter unduly, parameter selection lacks rationale.In addition, computational accuracy also can be by hole
Porosity composes morphology influence, and method applicability is restricted.
Therefore, it is necessary to need to establish new secondary porosity acquisition methods, secondary porosity precision is improved to reach
Demand.
Invention content
The technical problem to be solved in the present invention is, provide it is a kind of improve reliability and accuracy based on scale Electrical imaging
The method that secondary porosity spectrum cutoff value obtains secondary porosity.
The technical solution adopted by the present invention to solve the technical problems is：It provides a kind of based on scale Electrical imaging secondary pore
The method that degree spectrum cutoff value obtains secondary porosity, includes the following steps：
S1, secondary porosity spectrum cutoff value calculation formula is established：
Secondary porosity composes cutoff value=a × δ+porosity intermediate value
Wherein, a is fixed coefficient, and δ is in minimal amount of porosity and porosity in selected window on electric imaging logging image
Porosity variance in the range of value, porosity intermediate value are the porosity intermediate value in selected window；
S2, rock core secondary porosity is calculated using rock core nuclear magnetic resonance T 2 spectrum；
S3, the rock core secondary porosity scale secondary porosity spectrum cutoff value that gained is calculated using rock core nuclear magnetic resonance T 2 spectrum
The numerical value of fixed coefficient a in calculation formula：
Cutoff value calculation formula is composed using the secondary porosity in step S1, adjusts the number of fixed coefficient a between 1 and 3
Value, control secondary porosity spectrum cutoff value result of calculation, to the area that the porosity spectral envelope line after cutoff value is surrounded into
Row integral and calculating secondary porosity obtains porosity spectrum secondary porosity；By the porosity gradually calculated spectrum secondary porosity with
The secondary porosity that step S2 is obtained carries out error analysis, and corresponding a values is determine value, corresponding hole after scale during error minimum
Porosity spectrum secondary porosity result of calculation is secondary porosity after scale.
Preferably, in step S1, the selected window on electric imaging logging image, minimal amount of porosity and hole in calculation window
Spend the sum of Porosimeter digit rate in median value range；The counting rate for asking for each porosity in the range of this accounts for the percentage of gross-count rate
Than calculating minimal amount of porosity and the porosity variance δ in porosity median value range accordingly；One fixed coefficient a of setting is multiplied by this
Variance obtains secondary porosity spectrum cutoff value, obtains secondary porosity spectrum cutoff value calculation formula.
Preferably, in step S2, the large aperture hole stream discharged after centrifugation is equal to according to the size of rock core secondary porosity
Body volume according to sharp peaks characteristic in rock core nuclear magnetic resonance T 2 spectrum distribution map, obtains before centrifuging in porosity component the relaxation time most
The long peak time started, which corresponds to secondary porosity spectrum T2 cutoff values, before obtaining rock core secondary porosity=centrifugation
Porosity is accumulated before accumulation total porosity-corresponding centrifugation of secondary porosity spectrum T2 cutoff values.
Preferably, it is further comprising the steps of：
S4, casting body flake secondary porosity is calculated；
S5, the precision of porosity spectrum secondary porosity obtained in casting body flake secondary porosity verification step S3 is utilized.
Preferably, in step S4, the corresponding reservoir of selected the window's position, system on electric imaging logging image are taken in step S1
Standby casting body flake；
According to casting body flake qualification result, casting body flake type of secondary pore and corresponding relative amount, fitting casting are obtained
The relationship of body thin slice total pore surface rate and total porosity obtains the conversion relational expression of casting body flake total pore surface rate and total porosity, will
Secondary pore Areal porosity is scaled casting body flake secondary porosity.
Preferably, it is molten to include intergranular dissolved pore, intragranular dissolved pore, matrix dissolution pore, construction for the casting body flake type of secondary pore
Erosion seam.
Preferably, in step S5, by between casting body flake secondary porosity and porosity spectrum calculating secondary porosity
Mean absolute error verifies the accuracy of the porosity spectrum secondary porosity of acquisition.
Preferably, in step S5, the mean absolute error is less than 5%
The present invention provides another method for obtaining secondary porosity based on scale Electrical imaging secondary porosity spectrum cutoff value,
Include the following steps：
S1, secondary porosity spectrum cutoff value calculation formula is established：
Secondary porosity composes cutoff value=a+ average pores
Wherein, a is fixed coefficient, and average pore is all porositys in selected window on electric imaging logging image
Arithmetic average；
S2, rock core secondary porosity is calculated using rock core nuclear magnetic resonance T 2 spectrum；
S3, the rock core secondary porosity scale secondary porosity spectrum cutoff value that gained is calculated using rock core nuclear magnetic resonance T 2 spectrum
The numerical value of fixed coefficient a in calculation formula：
Cutoff value calculation formula is composed using the secondary porosity in step S1, adjusts the number of fixed coefficient a between 1 and 3
Value, control secondary porosity spectrum cutoff value result of calculation, to the area that the porosity spectral envelope line after cutoff value is surrounded into
Row integral and calculating secondary porosity obtains porosity spectrum secondary porosity；By the porosity gradually calculated spectrum secondary porosity with
The secondary porosity that step S2 is obtained carries out error analysis, and corresponding a values is determine value, corresponding hole after scale during error minimum
Porosity spectrum secondary porosity result of calculation is secondary porosity after scale.
Preferably, in step S2, the large aperture hole stream discharged after centrifugation is equal to according to the size of rock core secondary porosity
Body volume according to sharp peaks characteristic in rock core nuclear magnetic resonance T 2 spectrum distribution map, obtains before centrifuging in porosity component the relaxation time most
The long peak time started, which corresponds to secondary porosity spectrum T2 cutoff values, before obtaining rock core secondary porosity=centrifugation
Porosity is accumulated before accumulation total porosity-corresponding centrifugation of secondary porosity spectrum T2 cutoff values.
Beneficial effects of the present invention：Using rock core nuclear magnetic resonance T 2 spectrum scale, obtained using the method for Approach by inchmeal secondary
Porosity spectrum cutoff value calculation formula in fixed coefficient a, be significantly greatly increased secondary porosity spectrum cutoff value result of calculation can
By property, and then obtain the high secondary porosity acquisition methods of accuracy.
Using casting body flake secondary porosity verification porosity spectrum secondary porosity obtain as a result, effective guarantee is secondary
Porosity obtains the accuracy of result.
Description of the drawings
Below in conjunction with accompanying drawings and embodiments, the invention will be further described, in attached drawing：
Fig. 1 is first embodiment of the invention based on scale Electrical imaging secondary porosity spectrum cutoff value acquisition secondary porosity
Method flow diagram；
Fig. 2 is the carbonate reservoir rock core nuclear magnetic resonance T 2 spectrum distribution map in the present invention；
Fig. 3 is the fit correlation figure of casting body flake total pore surface (seam) rate and total porosity in the present invention.
Specific embodiment
With reference to figure 1, first embodiment of the invention obtains secondary hole based on scale Electrical imaging secondary porosity spectrum cutoff value
The method of porosity, includes the following steps：
S1, secondary porosity spectrum cutoff value calculation formula is established：
Secondary porosity composes cutoff value=a × δ+porosity intermediate value
Wherein, a is fixed coefficient；δ is in minimal amount of porosity and porosity in selected window on electric imaging logging image
Porosity variance in the range of value；Porosity intermediate value be selected window porosity intermediate value, cumulative porosity degree percentile
Reach porosity corresponding when 50%, physical significance is greater than its porosity and accounts for 50%, is also accounted for less than its porosity
50%.
Specifically, in step S1, the selected window on electric imaging logging image, minimal amount of porosity and hole in calculation window
Spend the sum of Porosimeter digit rate in median value range；Ask in window each hole in minimal amount of porosity and porosity median value range
The counting rate of degree accounts for the percentage of gross-count rate, calculates minimal amount of porosity and the porosity variance in porosity median value range accordingly
δ；One fixed coefficient a of setting is multiplied by the variance, obtains secondary porosity spectrum cutoff value, obtains secondary porosity spectrum cutoff value meter
Calculate formula.
S2, rock core secondary porosity is calculated using rock core nuclear magnetic resonance T 2 spectrum.
Specifically, the large aperture hole discharged after (or being approximately equal to) centrifugation is equal to according to the size of rock core secondary porosity
Fluid volume according to sharp peaks characteristic in rock core nuclear magnetic resonance T 2 spectrum distribution map, obtains centrifuging in preceding porosity component the relaxation time
Longest peak time started, the time started correspond to secondary porosity spectrum T2 cutoff values, obtain rock core secondary porosity=centrifugation
Porosity is accumulated before preceding accumulation total porosity-corresponding centrifugation of secondary porosity spectrum T2 cutoff values.
S3, fixed coefficient in cutoff value calculation formula is composed using rock core nuclear magnetic resonance secondary porosity scale secondary porosity
The numerical value of a：
Cutoff value calculation formula is composed using the secondary porosity in step S1, adjusts the number of fixed coefficient a between 1 and 3
Value, control secondary porosity spectrum cutoff value result of calculation, to the area that the porosity spectral envelope line after cutoff value is surrounded into
Row integral and calculating secondary porosity obtains porosity spectrum secondary porosity；By the porosity gradually calculated spectrum secondary porosity with
The secondary porosity that step S2 is obtained carries out error analysis, and corresponding a values is determine value, corresponding hole after scale during error minimum
Porosity spectrum secondary porosity result of calculation is secondary porosity after scale.
Further, first embodiment of the invention based on scale Electrical imaging secondary porosity spectrum cutoff value obtains secondary hole
The method of porosity is further comprising the steps of：
S4, casting body flake secondary porosity is calculated.
Wherein, selected window on electric imaging logging image in corresponding step S1, prepares casting body flake.
According to casting body flake qualification result, casting body flake type of secondary pore and relative amount are obtained, is fitted casting body flake
The relationship of total pore surface rate and total porosity obtains the conversion relational expression of casting body flake total pore surface rate and total porosity, by secondary hole
Gap Areal porosity is scaled casting body flake secondary porosity.Casting body flake type of secondary pore includes intergranular dissolved pore, intragranular dissolved pore, base
Matter dissolution pore, construction corrosion seam.Above-mentioned casting body flake total pore surface rate includes poroid in casting body flake and gap-like hole total pore surface rate.
S5, the precision of porosity spectrum secondary porosity obtained in casting body flake secondary porosity verification step S3 is utilized：
Mean absolute error between secondary porosity, verification are calculated by casting body flake secondary porosity and porosity spectrum
The accuracy (i.e. precision) of the porosity spectrum secondary porosity of acquisition.
If verification result meets the precision of related request, export secondary porosity and obtain result；If it is not, then return to step
S3 determines secondary porosity spectrum cutoff value again, until meeting acquisition precision.
The side that secondary porosity is obtained based on scale Electrical imaging secondary porosity spectrum cutoff value of second embodiment of the invention
Method includes the following steps：
S1, secondary porosity spectrum cutoff value calculation formula is established：
Secondary porosity composes cutoff value=a+ average pores
Wherein, a is fixed coefficient, and average pore is all porositys in selected window on electric imaging logging image
Arithmetic average.
S2, rock core secondary porosity is calculated using rock core nuclear magnetic resonance T 2 spectrum.
Specifically, the large aperture hole discharged after (or being approximately equal to) centrifugation is equal to according to the size of rock core secondary porosity
Fluid volume according to sharp peaks characteristic in rock core nuclear magnetic resonance T 2 spectrum distribution map, obtains centrifuging in preceding porosity component the relaxation time
Longest peak time started, the time started correspond to secondary porosity spectrum T2 cutoff values, obtain rock core secondary porosity=centrifugation
Porosity is accumulated before preceding accumulation total porosity-corresponding centrifugation of secondary porosity spectrum T2 cutoff values.
S3, the rock core secondary porosity scale secondary porosity spectrum cutoff value that gained is calculated using rock core nuclear magnetic resonance T 2 spectrum
The numerical value of fixed coefficient a in calculation formula：
Cutoff value calculation formula is composed using the secondary porosity in step S1, adjusts the number of fixed coefficient a between 1 and 3
Value, control secondary porosity spectrum cutoff value result of calculation, to the area that the porosity spectral envelope line after cutoff value is surrounded into
Row integral and calculating secondary porosity obtains porosity spectrum secondary porosity；By the porosity gradually calculated spectrum secondary porosity with
The secondary porosity that step S2 is obtained carries out error analysis, and corresponding a values is determine value, corresponding hole after scale during error minimum
Porosity spectrum secondary porosity result of calculation is secondary porosity after scale.
Further, second embodiment of the invention based on scale Electrical imaging secondary porosity spectrum cutoff value obtains secondary hole
The method of porosity is further comprising the steps of：
S4, casting body flake secondary porosity is calculated.
S5, the precision of porosity spectrum secondary porosity obtained in casting body flake secondary porosity verification step S3 is utilized.
Step S4 and S5 concrete operations can refer to above-mentioned first embodiment it is related described in.
Below by taking Zhujiangkou Basin stream flower 11-1 biohermal limestones oil field as an example, target zone is carried out with LH11-1-D5P1 wells
Section secondary porosity elaborates the present invention for obtaining：
Step 1 establishes secondary porosity spectrum cutoff value calculation formula：The selected window on electric imaging logging image, window
Size is determined according to actual conditions.Porosimeter digit rate in calculation window in minimal amount of porosity and porosity median value range it
With the counting rate for asking for each porosity in the range of this accounts for the percentage of gross-count rate, calculates minimal amount of porosity and hole accordingly
The porosity variance in median value range is spent, one fixed coefficient a of setting is multiplied by the variance, obtains secondary porosity spectrum cutoff value.
Calculation formula is：
Secondary porosity composes cutoff value=a × δ+porosity intermediate value
In formula：δ-minimal amount of porosity and the porosity variance in porosity median value range；
A-fixed coefficient, the coefficient determine secondary porosity spectrum cutoff value computational accuracy.
Step 2 calculates rock core secondary porosity using rock core nuclear magnetic resonance T 2 spectrum：As shown in Fig. 2, it is stored up for carbonate rock
Layer rock core nuclear magnetic resonance T 2 spectrum distribution map, curve 1 are porosity component before rock sample centrifugation, show apparent three peak feature, one
As think that relaxation time longer third peak represents large aperture hole, i.e. secondary pore；Curve 2 is porosity after rock sample centrifugation
Apparent double-peak feature is presented in component, and relaxation time longest third peak disappears in porosity component before centrifugation, illustrates macropore
Fluid in diameter hole is because centrifugal action is discharged.Curve 3 and curve 4 are respectively to be accumulated after accumulating and centrifuge before centrifuging.
Secondary porosity size is approximately equal to the large aperture pore-fluid volume discharged after centrifugation.Find porosity before centrifuging
Longest third peak time started in relaxation time in component, the time correspond to secondary porosity spectrum T2 cutoff values, tire out before centrifugation
Accumulation porosity is secondary porosity before product total porosity subtracts the corresponding centrifugation of secondary porosity spectrum T2 cutoff values.
Step 3 uses the fixed coefficient a in rock core nuclear magnetic resonance secondary porosity scale secondary porosity calculation formula
Value：Cutoff value calculation formula is composed using the secondary porosity in step 1, by adjusting the number of the fixed coefficient a in calculation formula
Value controls secondary porosity spectrum cutoff value result of calculation (between 1-3), the porosity spectral envelope line after cutoff value is surrounded
Area carry out integral and calculating secondary porosity.
Porosity spectrum secondary porosity and the nuclear magnetic resonance T 2 spectrum that LH11-1-D5P1 well intervals of interest are gradually calculated calculate
Secondary porosity carries out error analysis, and when a takes 1.01, secondary porosity and the nuclear magnetic resonance T 2 spectrum calculating being calculated are secondary
Porosity error is minimum, thereby determines that the fixed coefficient a in secondary porosity spectrum cutoff value formula is 1.01, completes to secondary hole
The scale of porosity.
Step 4 calculates casting body flake secondary porosity：It is identified according to LH11-1-D5P1 well 1399.5m positions casting body flake
As a result, show that thin slice intergranular dissolved pore, intragranular dissolved pore, matrix dissolution pore, construction corrosion seam are respectively 2%, 5%, 2%, 3%, i.e., it is secondary
Dough hole (seam) rate is 12%.
Fig. 3 is the fit correlation figure of casting body flake total pore surface (seam) rate and total porosity, it follows that casting body flake face
The conversion relational expression y=0.61x+6.75 of (seam) rate and one embodiment of total porosity changes secondary pore face (seam) rate
It calculates as casting body flake secondary porosity.
Step 5, the precision using casting body flake secondary porosity verification porosity spectrum secondary porosity：It is as shown in table 1 below
It calculates casting body flake secondary porosity and porosity spectrum in LH11-1-D5P1 well 1389-1470m intervals of interest and calculates secondary pore
Mean absolute error and average relative error between degree are respectively 0.86% and 7.7%, it was demonstrated that the secondary porosity spectrum of acquisition
Precision it is higher.
1 secondary porosity precision analysis of table
It is compared by the way that two kinds of scale methods of casting body flake and nuclear magnetic resonance are obtained secondary porosity, porosity spectrum obtains
Secondary porosity is taken to coincide and be on speaking terms with casting body flake secondary porosity, nuclear magnetic resonance T 2 spectrum calculating secondary porosity.
The foregoing is merely the embodiment of the present invention, are not intended to limit the scope of the invention, every to utilize this hair
The equivalent structure or equivalent flow shift that bright specification and accompanying drawing content are made directly or indirectly is used in other relevant skills
Art field, is included within the scope of the present invention.
Claims (10)
- A kind of 1. method that secondary porosity is obtained based on scale Electrical imaging secondary porosity spectrum cutoff value, which is characterized in that packet Include following steps：S1, secondary porosity spectrum cutoff value calculation formula is established：Secondary porosity composes cutoff value=a × δ+porosity intermediate valueWherein, a is fixed coefficient, and δ is minimal amount of porosity in selected window on electric imaging logging image and porosity intermediate value model Interior porosity variance is enclosed, porosity intermediate value is the porosity intermediate value in selected window；S2, rock core secondary porosity is calculated using rock core nuclear magnetic resonance T 2 spectrum；S3, the rock core secondary porosity scale secondary porosity spectrum cutoff value calculating that gained is calculated using rock core nuclear magnetic resonance T 2 spectrum The numerical value of fixed coefficient a in formula：Cutoff value calculation formula is composed using the secondary porosity in step S1, adjusts the numerical value of fixed coefficient a, control between 1 and 3 Secondary porosity processed composes cutoff value result of calculation, and the area surrounded to the porosity spectral envelope line after cutoff value integrates Secondary porosity is calculated, obtains porosity spectrum secondary porosity；By the porosity gradually calculated spectrum secondary porosity and step S2 The secondary porosity of acquisition carries out error analysis, and corresponding a values is determine value, corresponding porosity spectrum after scale during error minimum Secondary porosity result of calculation is secondary porosity after scale.
- 2. the side according to claim 1 that secondary porosity is obtained based on scale Electrical imaging secondary porosity spectrum cutoff value Method, which is characterized in that in step S1, the selected window on electric imaging logging image, minimal amount of porosity and hole in calculation window Spend the sum of Porosimeter digit rate in median value range；The counting rate for asking for each porosity in the range of this accounts for the percentage of gross-count rate Than calculating minimal amount of porosity and the porosity variance δ in porosity median value range accordingly；One fixed coefficient a of setting is multiplied by this Variance obtains secondary porosity spectrum cutoff value, obtains secondary porosity spectrum cutoff value calculation formula.
- 3. the side according to claim 1 that secondary porosity is obtained based on scale Electrical imaging secondary porosity spectrum cutoff value Method, which is characterized in that in step S2, the large aperture pore-fluid discharged after centrifugation is equal to according to the size of rock core secondary porosity Volume according to sharp peaks characteristic in rock core nuclear magnetic resonance T 2 spectrum distribution map, obtains centrifuging relaxation time longest in preceding porosity component The peak time started, which corresponds to secondary porosity spectrum T2 cutoff values, is accumulated before obtaining rock core secondary porosity=centrifugation Porosity is accumulated before total porosity-corresponding centrifugation of secondary porosity spectrum T2 cutoff values.
- 4. being based on scale Electrical imaging secondary porosity spectrum cutoff value according to claim 1-3 any one of them obtains secondary pore The method of degree, which is characterized in that further comprising the steps of：S4, casting body flake secondary porosity is calculated；S5, the precision of porosity spectrum secondary porosity obtained in casting body flake secondary porosity verification step S3 is utilized.
- 5. the side according to claim 4 that secondary porosity is obtained based on scale Electrical imaging secondary porosity spectrum cutoff value Method, which is characterized in that in step S4, take in step S1 the corresponding reservoir of selected the window's position, system on electric imaging logging image Standby casting body flake；According to casting body flake qualification result, casting body flake type of secondary pore and corresponding relative amount are obtained, the fitting body of casting is thin The relationship of piece total pore surface rate and total porosity obtains the conversion relational expression of casting body flake total pore surface rate and total porosity, will be secondary Hole Areal porosity is scaled casting body flake secondary porosity.
- 6. the side according to claim 5 that secondary porosity is obtained based on scale Electrical imaging secondary porosity spectrum cutoff value Method, which is characterized in that the casting body flake type of secondary pore includes intergranular dissolved pore, intragranular dissolved pore, matrix dissolution pore, construction corrosion Seam.
- 7. the side according to claim 4 that secondary porosity is obtained based on scale Electrical imaging secondary porosity spectrum cutoff value Method, which is characterized in that in step S5, by flat between casting body flake secondary porosity and porosity spectrum calculating secondary porosity Equal absolute error verifies the accuracy of the porosity spectrum secondary porosity of acquisition.
- 8. the side according to claim 7 that secondary porosity is obtained based on scale Electrical imaging secondary porosity spectrum cutoff value Method, which is characterized in that in step S5, the mean absolute error is less than 5%.
- A kind of 9. method that secondary porosity is obtained based on scale Electrical imaging secondary porosity spectrum cutoff value, which is characterized in that packet Include following steps：S1, secondary porosity spectrum cutoff value calculation formula is established：Secondary porosity composes cutoff value=a+ average poresWherein, a is fixed coefficient, and average pore is that all porositys count in selected window on electric imaging logging image Average value；S2, rock core secondary porosity is calculated using rock core nuclear magnetic resonance T 2 spectrum；S3, the rock core secondary porosity scale secondary porosity spectrum cutoff value calculating that gained is calculated using rock core nuclear magnetic resonance T 2 spectrum The numerical value of fixed coefficient a in formula：Cutoff value calculation formula is composed using the secondary porosity in step S1, adjusts the numerical value of fixed coefficient a, control between 1 and 3 Secondary porosity processed composes cutoff value result of calculation, and the area surrounded to the porosity spectral envelope line after cutoff value integrates Secondary porosity is calculated, obtains porosity spectrum secondary porosity；By the porosity gradually calculated spectrum secondary porosity and step S2 The secondary porosity of acquisition carries out error analysis, and corresponding a values is determine value, corresponding porosity spectrum after scale during error minimum Secondary porosity result of calculation is secondary porosity after scale.
- 10. the side according to claim 9 that secondary porosity is obtained based on scale Electrical imaging secondary porosity spectrum cutoff value Method, which is characterized in that in step S2, the large aperture pore-fluid discharged after centrifugation is equal to according to the size of rock core secondary porosity Volume according to sharp peaks characteristic in rock core nuclear magnetic resonance T 2 spectrum distribution map, obtains centrifuging relaxation time longest in preceding porosity component The peak time started, which corresponds to secondary porosity spectrum T2 cutoff values, is accumulated before obtaining rock core secondary porosity=centrifugation Porosity is accumulated before total porosity-corresponding centrifugation of secondary porosity spectrum T2 cutoff values.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711271137.0A CN108226002B (en) | 2017-12-05 | 2017-12-05 | Method for obtaining secondary porosity based on cutoff value of secondary porosity spectrum of scale electrography |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711271137.0A CN108226002B (en) | 2017-12-05 | 2017-12-05 | Method for obtaining secondary porosity based on cutoff value of secondary porosity spectrum of scale electrography |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108226002A true CN108226002A (en) | 2018-06-29 |
| CN108226002B CN108226002B (en) | 2020-08-18 |
Family
ID=62653839
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711271137.0A Active CN108226002B (en) | 2017-12-05 | 2017-12-05 | Method for obtaining secondary porosity based on cutoff value of secondary porosity spectrum of scale electrography |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108226002B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020236784A1 (en) * | 2019-05-20 | 2020-11-26 | Saudi Arabian Oil Company | Methods and systems for measuring pore volume compressibility with low field nuclear magnetic resonance techniques |
| US11614417B2 (en) | 2021-07-06 | 2023-03-28 | Saudi Arabian Oil Company | Determining saturation in low resistivity pay zones |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5359194A (en) * | 1992-05-01 | 1994-10-25 | Texaco Inc. | X-ray CT measurement of secondary (vugular) porosity in reservoir core material |
| CN102141637A (en) * | 2010-01-28 | 2011-08-03 | 中国石油天然气股份有限公司 | Method for continuously and quantitatively evaluating pore structure of reservoir by utilizing nuclear magnetic resonance logging data |
| CN105866009A (en) * | 2016-05-30 | 2016-08-17 | 中国石油大学(北京) | Method and device for calculating effective porosity of compact oil storage layer |
-
2017
- 2017-12-05 CN CN201711271137.0A patent/CN108226002B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5359194A (en) * | 1992-05-01 | 1994-10-25 | Texaco Inc. | X-ray CT measurement of secondary (vugular) porosity in reservoir core material |
| CN102141637A (en) * | 2010-01-28 | 2011-08-03 | 中国石油天然气股份有限公司 | Method for continuously and quantitatively evaluating pore structure of reservoir by utilizing nuclear magnetic resonance logging data |
| CN105866009A (en) * | 2016-05-30 | 2016-08-17 | 中国石油大学(北京) | Method and device for calculating effective porosity of compact oil storage layer |
Non-Patent Citations (2)
| Title |
|---|
| 史飞洲: "碳酸盐岩地层电成像测井孔隙度谱截止值计算方法", 《测井技术》 * |
| 谢蕾蕾: "基于成像测井孔隙度谱的储层次生孔隙度评价", 《能源与环保》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020236784A1 (en) * | 2019-05-20 | 2020-11-26 | Saudi Arabian Oil Company | Methods and systems for measuring pore volume compressibility with low field nuclear magnetic resonance techniques |
| US10948439B2 (en) | 2019-05-20 | 2021-03-16 | Saudi Arabian Oil Company | Methods and systems for measuring pore volume compressibility with low field nuclear magnetic resonance techniques |
| US11293886B2 (en) | 2019-05-20 | 2022-04-05 | Saudi Arabian Oil Company | Methods and systems for measuring pore volume compressibility with low field nuclear magnetic resonance techniques |
| US11614417B2 (en) | 2021-07-06 | 2023-03-28 | Saudi Arabian Oil Company | Determining saturation in low resistivity pay zones |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108226002B (en) | 2020-08-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2447469B1 (en) | Quantitative calculation method for hydrocarbon saturation of fractured reservoir | |
| CN106066901B (en) | A kind of datum mark method for analyzing stability of GNSS automation deformation monitoring | |
| CN104278989B (en) | Method for obtaining saturation index of low-porosity and low-permeability reservoir | |
| CN109100793B (en) | Method for quantitatively analyzing influence of fracture factors on reservoir | |
| CN102434152A (en) | Method for calculating oil saturation of reservoir | |
| CN105525917B (en) | The method for identifying fluid type in the method and identification reservoir of reservoir | |
| CN110955982B (en) | Metamorphic rock reservoir permeability calculation method and device and computer storage medium | |
| CN106285664A (en) | The dual media reservoir rock water saturation of percolation-based network analog calculates method | |
| CN107944156B (en) | The choosing method of wave height threshold value | |
| CN112363242B (en) | Reservoir fluid identification method and device based on logging fusion | |
| CN103437759B (en) | Method for non-experimentally measuring cut-off value of natural gas layer T2 | |
| CN108226002A (en) | The method that secondary porosity is obtained based on scale Electrical imaging secondary porosity spectrum cutoff value | |
| CN105447762A (en) | Calculation method for low permeability reservoir flooding information of fluid replacement | |
| CN109386281A (en) | A method of obtaining Fractured low porosity and low permeability reservoir well logging saturation degree | |
| CN106204302B (en) | Method for calculating original water saturation and application | |
| CN103995297A (en) | Method for correcting resistivity of crack carbonatite stratum | |
| CN105401937B (en) | Saturation index prediction method based on pore structure | |
| CN106202763A (en) | The computational methods of dual media reservoir rock water saturation | |
| CN106802432A (en) | The method for surveying and device of soil layer construction | |
| CN109031436A (en) | A kind of hypotonic loose sand reservoir water saturation evaluation method of low-resistance based under Reservoir Classification | |
| CN109814174B (en) | Comprehensive well logging identification method for clastic rock unconformity structure body | |
| CN107725044A (en) | The method of Sandstone Gas Bearing reservoir producing water ratio prediction based on array induction, laterolog | |
| CN106320255B (en) | A kind of computational methods of Flood Plain Flow swale inflow-rate of water turbine | |
| CN105043390A (en) | Gravity field interpolation method based on universal kriging | |
| CN104374679A (en) | Method for determining permeability contribution value of reservoir |
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 |