CN110320139A - Fractured-vuggy reservoir fracture porosity quantitative evaluation method and system - Google Patents
Fractured-vuggy reservoir fracture porosity quantitative evaluation method and system Download PDFInfo
- Publication number
- CN110320139A CN110320139A CN201810286736.8A CN201810286736A CN110320139A CN 110320139 A CN110320139 A CN 110320139A CN 201810286736 A CN201810286736 A CN 201810286736A CN 110320139 A CN110320139 A CN 110320139A
- Authority
- CN
- China
- Prior art keywords
- porosity
- fracture porosity
- resistivity
- crack
- fractured
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- 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
Abstract
The invention discloses a kind of fractured-vuggy reservoir fracture porosity quantitative evaluation method and systems, comprising: is handled using rock core capillary pressure data of the J function to continuous depth, obtains continuous depth J functional value；Relational expression based on J function and water saturation obtains continuous depth water saturation value；Based on continuous depth water saturation value, base portion resistivity is calculated；Based on series-parallel electrical conduction model, pure crack rock al-lateral resistivity is calculated；Using crack angle discriminate and bilaterally fracture porosity formula, fracture porosity is calculated.Fractured-vuggy reservoir fracture porosity quantitative evaluation method and system provided by the invention, being capable of fracture porosity progress fast and accurately quantitative assessment.
Description
Technical field
The invention belongs to petroleum and geology field, quantitatively comment more particularly, to a kind of fractured-vuggy reservoir fracture porosity
Valence method and system.
Background technique
In Exploration of Oil And Gas, carbonate rock fractured cave type reservoir occupies an important position, and there is reservoir productivity in crack
The control action of positive important, as one of the important parameter for evaluating such reservoir, fracture porosity being sought and demarcating always
It is global problem.
Sibbit A. and Faiver O. (Sibbit A, Faivre O.The dual laterolog response in
fractured rocks[C].SPWLA 26th Annual Logging Symposium.Dallas,Texas,1985,17-
20.) three-dimensional finite element numerical computation method is utilized, vertical fracture and horizontal fracture aperture computation model are established；Luo Zhenyao (sieve
Primary Study [J] geophysical log of the chastity credit lateral material computation fracture opening of, 1990,14 (2): 83-92.) base
Visual angle degree fracture aperture model is established in dual laterolog data；Subsequent Philippe A.P. and Roger N.A.
(Philippe A.P.,Roger N.A.In situ measurements of electrical resistivity,
formation anisotropy and tectonic context[C].SPWLA 31st Annual Logging
Symposium, Lafayette, Louisiana, 1990,24-27.) forward modeling is derived and bilaterally splits in arbitrary angle
The log response of seam, and obtained under the conditions of quasi- horizontal joint and quasi- vertical seam the bilaterally relationship of response and fracture opening；It is logical
It crosses bilaterally depth apparent resistivity difference and judges crack angle, but the foundation of the model is carried out under conditions of more special
, the response pattern in intermediate angle crack is not summarized further；In grinding for Sibbit A. and Philippe A.P.
On the basis of studying carefully, (Li Shanjun, Xiao Chengwen, Wang Hanming wait mathematical model and the crack of the dual laterolog response in the crack to Li Shanjun
Quantitative interpretation [J] Chinese Journal of Geophysics of porosity, 1996,39 (6): 845-852.) et al. utilize three-dimensional finite element numerical value meter
In the case of calculation method and the method for conductivity tensor simulate different angle respectively, single crack, equidistant parallel fracture, friendship
The dual laterolog response of fork seam and macroscopic anisotropy fracture formation, and summarize the log response value in various situation cracks
With the relationship of fracture dip, crack is further divided into three classes, i.e. high angle frature, dipping fracture and low angle frature, and specific
Give the method for inverting fracture porosity and fracture dip in various situations；(the crack Deng Shaogui, Lu Chunli-hole such as Deng Shaogui
The pre-test of dual laterolog means of interpretation [J] logging technique of gap double porosity media, 2010,34 (6), 537-541.) establish bilateral
To log response and matrix conductivity, fracture porosity, the relationship at inclination angle；(Qin Qirong, Huangping brightness, Zhou Zhiyuan wait to Qin Qirong
Application [J] natural gas geoscience of the full diameter sample analysis in well log interpretation fracture porosity, 2005,16 (5): 637-
Etc. 640) full diameter sample is utilized to carry out fracture porosity calculating；(Lv Hongzhi, Lu Yunlong, Cui Yunjiang wait improved to Lv Hongzhi
Porosity model evaluates fluid properties and fracture porosity [J] applied acoustics, 2016,35 (4): 351-356) it is managed according to Biot
By, using array sonic log data carry out inverting obtain fracture porosity；Zhang Ying (3 kind fire of Zhang Ying, the Pan Baozhi based on FMI
Mountain rock Reservoir Fracture porosity acquiring method [J] logging technique, 2012,36 (4): 365-369) utilize FMI analysis three
Kind Volcanic Rock porosity acquiring method；(Shang Genhua, Hou Xiaochun, Liu Xuewei wait to study Daqing oil using X-CT to Shang Genhua
Field dual media crack and sandstone porosity [J] .CT theory and application research, 1997,6 (1): 26-31) it is had studied using X-CT
Dual media crack and sandstone porosity；(application of the Liu Libo tri-porosity logging method in identification crack is [interior by Liu Libo
Mongolian petrochemical industry] .2014,11:155-156) by carrying out unification processing to three porosity, introduce crack indicative character system
Number reaches fracture identification and porosity calculation；(fracture reservoir porosity, saturation computation are rested the head on to Qin Haixu by Qin Haixu, Wu state
Method [J] geophysical prospecting for oil, 2015,50 (1): 103-110.) with fracture medium anisotropic gradient calculating crack
The new method of medium porosity and saturation degree.
The existing achievement of comprehensive analysis forefathers, the evaluation for fracture porosity are concentrated mainly in five class methods, 1. bases
In the calculation method of dual laterolog data, this method disadvantage is that undetermined parameter is numerous, and response processing bilaterally is split
Seam is also influenced by Fluid in Pore, multi-solution is strong outside influencing；2. being based on resistivity imaging logging document method, this method limitation
Property be somewhat expensive, not all well all carries out imaging measurement, while the crack that Electrical imaging calculates is seen on image
Crack, the crack that can't see on image is simultaneously not involved in calculating；3. being based on three porosity method, this method difficult point is rock bone
The accurate selection of frame parameter, and when the high angle of crack, and be not suitable for；4. being based on bulk method, this method disadvantage exists
In corrections such as crack section coring success rate are difficult to ensure, and need to be compacted rock core, playback；5. rock core X-CT scanning method,
This method is computationally intensive, and representative lithology sample obtains difficulty.
Therefore, it is necessary to provide a kind of fractured-vuggy reservoir fracture porosity quantitative evaluation method and system, solves oil gas and survey
The fine quantitative assessment problem of fracture-cavity type carbonate Reservoir Fracture porosity in spy.
Summary of the invention
The present invention is by providing a kind of fractured-vuggy reservoir fracture porosity quantitative evaluation method and system, it is therefore an objective to by right
Core sample experimental analysis corrects and demarcates dual laterolog data, rejects base portion to influence bilaterally, obtains pure split
Partial ohmic rate value is stitched, and then obtains the higher fracture porosity of precision, finally realizes and utilizes dual laterolog data counterincision
Slot apertures porosity carries out fast and accurately quantitative assessment.
According to an aspect of the invention, it is proposed that a kind of fractured-vuggy reservoir fracture porosity quantitative evaluation method, this method
Include:
1) it is handled using rock core capillary pressure data of the J function to continuous depth, obtains continuous depth J functional value；
2) relational expression based on the J function and water saturation, obtains continuous depth water saturation value；
3) it is based on the continuous depth water saturation value, calculates base portion resistivity；
4) it is based on series-parallel electrical conduction model, calculates pure crack rock al-lateral resistivity；
5) crack angle discriminate and bilaterally fracture porosity formula are utilized, fracture porosity is calculated.
Preferably, the J function representation are as follows:
Wherein, PcFor hollow billet curve capillary pressure value；σ is interfacial tension；θ is angle of wetting；K is permeability；For hole
Degree.
Preferably, the relational expression of the J function and water saturation indicates are as follows:
Sw=cedJ (2)
Wherein, SwFor water saturation；C, d is constant；J is J function.
Preferably, in step 3), base portion resistivity is calculated using Archie formula, the Archie formula indicates
Are as follows:
Wherein, RtbFor base portion resistivity；A is lithology factor；B is lithology dependent constant；RwFor formation water resistivity；
SwFor water saturation；N is saturation exponent；For porosity；M is lithologic index.
Preferably, the series-parallel electrical conduction model indicates are as follows:
Wherein, RdFor deep lateral resistivity log value；RdfFor pure crack rock deep lateral resistivity；RsFor shallow lateral resistance
Rate log value；RsfFor pure crack rock shallowly lateral resistivity；RtbFor base portion resistivity；For porosity.
Preferably, step 5) includes:
5.1) the pure crack rock al-lateral resistivity obtained based on step 4) obtains crack using crack angle discriminate
Type；
5.2) type in the crack judged based on step 5.1), the analysis of the confirmation bilaterally fracture porosity formula
Coefficient；
5.3) coefficient of analysis obtained based on step 5.2) obtains crack hole using the bilaterally fracture porosity formula
Porosity.
Preferably, the crack angle discriminate indicates are as follows:
Wherein, Y is discriminant index；RdfFor pure crack rock deep lateral resistivity；RsfFor pure crack rock shallowly lateral resistance
Rate.
Preferably, the bilaterally fracture porosity formula indicates are as follows:
Wherein,For fracture porosity；a1、a2And a3For coefficient of analysis；RdFor deep lateral resistivity；RsFor shallow lateral resistance
Rate；RmfFor mud resistivity.
According to another aspect of the invention, it is proposed that a kind of fractured-vuggy reservoir fracture porosity Quantitative Evaluation System, described
System includes: memory, is stored with computer executable instructions；Processor runs the executable finger of calculating on the memory
When enabling, the processor is performed the steps of
1) it is handled using rock core capillary pressure data of the J function to continuous depth, obtains continuous depth J functional value；
2) relational expression based on the J function and water saturation, obtains continuous depth water saturation value；
3) it is based on the continuous depth water saturation value, calculates base portion resistivity；
4) it is based on series-parallel electrical conduction model, calculates pure crack rock al-lateral resistivity；
5) crack angle discriminate and bilaterally fracture porosity formula are utilized, fracture porosity is calculated.
Preferably, the J function representation are as follows:
Wherein, PcFor hollow billet curve capillary pressure value；σ is interfacial tension；θ is angle of wetting；K is permeability；For hole
Degree.
The beneficial effects of the present invention are: by providing a kind of fractured-vuggy reservoir fracture porosity quantitative evaluation method and being
System corrects and demarcates dual laterolog data to core sample experimental analysis, rejects base portion to influence bilaterally, obtains
To pure rift portion resistivity value, and then the higher fracture porosity of precision is obtained, and then establish fracture porosity and quantitatively calculate
Model finally realizes and carries out fast and accurately quantitative assessment using dual laterolog data fracture porosity, solves oil
The fine quantitative assessment problem of fracture-cavity type carbonate Reservoir Fracture porosity in gas exploration.
Other features and advantages of the present invention will then part of the detailed description can be specified.
Detailed description of the invention
Exemplary embodiment of the invention is described in more detail in conjunction with the accompanying drawings, it is of the invention above-mentioned and its
Its purpose, feature and advantage will be apparent, wherein in exemplary embodiment of the invention, identical reference label
Typically represent same parts.
Fig. 1 shows the flow chart of fractured-vuggy reservoir fracture porosity quantitative evaluation method according to the present invention.
Fig. 2 shows the schematic diagrames of fractured-vuggy reservoir according to the present invention.
Fig. 3 shows J function according to an embodiment of the invention and water saturation relational graph.
Fig. 4 shows result map according to an embodiment of the invention.
Specific embodiment
The preferred embodiment of the present invention is described in more detail below.Although the following describe preferred implementations of the invention
Mode, however, it is to be appreciated that may be realized in various forms the present invention without that should be limited by the embodiments set forth herein.Phase
Instead, these embodiments are provided so that the present invention is more thorough and complete, and can be by the scope of the present invention completely
It is communicated to those skilled in the art.
Embodiment 1
In this embodiment, fractured-vuggy reservoir fracture porosity quantitative evaluation method according to the present invention may include:
1) it is handled using rock core capillary pressure data of the J function to continuous depth, obtains continuous depth J functional value；
2) relational expression based on J function and water saturation obtains continuous depth water saturation value；
3) it is based on continuous depth water saturation value, calculates base portion resistivity；
4) it is based on series-parallel electrical conduction model, calculates pure crack rock al-lateral resistivity；
5) crack angle discriminate and bilaterally fracture porosity formula are utilized, fracture porosity is calculated.
Fig. 2 shows the schematic diagrames of fractured-vuggy reservoir according to the present invention.As shown in Fig. 2, fracture-cavity type carbonate reservoir
It is made of base portion and crack, therefore, base portion hole and crack can all influence bilateral all to bilaterally there is shunting function
To resistivity.The embodiment is by providing a kind of fractured-vuggy reservoir fracture porosity quantitative evaluation method and system, by rock
The analysis of heart sample experiments, corrects and demarcates dual laterolog data, rejects base portion to influence bilaterally, obtains pure crack
Partial ohmic rate value, and then the higher fracture porosity of precision is obtained, it finally realizes and utilizes dual laterolog data fracture
Porosity carries out fast and accurately quantitative assessment.
Fig. 1 shows the flow chart of fractured-vuggy reservoir fracture porosity quantitative evaluation method according to the present invention.Join below
Examine the specific steps that fractured-vuggy reservoir fracture porosity quantitative evaluation method according to the present invention is described in detail in Fig. 1.
Step 101, it is handled using rock core capillary pressure data of the J function to continuous depth, obtains continuous depth J letter
Numerical value.
Specifically, capillary pressure data are that laboratory is measured according to certain several depth point, and the purpose of J function processing is exactly
Find the relational expression that can characterize entire well section.
In one example, J function representation are as follows:
Wherein, PcFor hollow billet curve capillary pressure value；σ is interfacial tension；θ is angle of wetting；K is permeability；For hole
Degree.
Specifically, porosityIt can be sought according to neutron density crossplot；Permeability K can pass through core porosity and permeability
Recurrence obtains computation model；σ and θ determines that when for oil water contact, σ cos θ takes 26, connects when for air water according to formation fluid situation
When touching, σ cos θ takes 50.
In addition, hollow billet curve capillary pressure value PcIt indicates are as follows:
Pc=0.01 (ρw-ρ0)H (8)
Wherein, ρwFor the density of oil；ρ0For the density of water；H is reservoir height.
Step 102, the relational expression based on J function and water saturation obtains continuous depth water saturation value.
In one example, the relational expression of J function and water saturation indicates are as follows:
Sw=cedJ (2)
Wherein, SwFor water saturation；C, d is constant；J is J function.
Specifically, c, d are constant, are nondimensional numbers, and specific value is obtained according to research area's experimental data, binding area
Experimental data can be obtained continuous depth water saturation value using formula (1), formula (2) and formula (8).
Step 103, it is based on continuous depth water saturation value, calculates base portion resistivity.
In one example, in step 103, base portion resistivity, Archie formula table are calculated using Archie formula
It is shown as:
Wherein, RtbFor base portion resistivity；A is lithology factor；B is lithology dependent constant；RwFor formation water resistivity；
SwFor water saturation；N is saturation exponent；For porosity；M is lithologic index.
Specifically, a, b, m and n are nondimensional number, can be obtained by rock-electric test；Obtain ground by water analysis of data again
Layer water resistance rate Rw, base portion resistivity can be calculated using formula (3).
Step 104, it is based on series-parallel electrical conduction model, calculates pure crack rock al-lateral resistivity.
In one example, series-parallel electrical conduction model indicates are as follows:
Wherein, RdFor deep lateral resistivity log value；RdfFor pure crack rock deep lateral resistivity；RsFor shallow lateral resistance
Rate log value；RsfFor pure crack rock shallowly lateral resistivity；RtbFor base portion resistivity；For porosity.
Specifically, the base portion resistivity calculated according to al-lateral resistivity log value and step 103, utilizes such as formula
(4) and series-parallel electrical conduction model shown in formula (5) al-lateral resistivity as caused by pure crack, can be obtained.
Step 105, using crack angle discriminate and bilaterally fracture porosity formula, fracture porosity is calculated.
In one example, step 105 includes:
Based on pure crack rock al-lateral resistivity, the type in crack is obtained using crack angle discriminate；
Based on the type in the crack judged, the coefficient of analysis of bilaterally fracture porosity formula is confirmed；
Coefficient of analysis based on acquisition obtains fracture porosity using bilaterally fracture porosity formula.
In one example, angle discriminate in crack indicates are as follows:
Wherein, Y is discriminant index；RdfFor pure crack rock deep lateral resistivity；RsfFor pure crack rock shallowly lateral resistance
Rate.
In one example, bilaterally fracture porosity formula indicates are as follows:
Wherein,For fracture porosity；a1、a2And a3For coefficient of analysis；RdFor deep lateral resistivity；RsFor shallow lateral resistance
Rate；RmfFor mud resistivity.
Specifically, the bilaterally crack hole using the crack angle discriminate as shown in formula (6) and as shown in formula (7)
Porosity formula, calculating fracture porosity includes:
Types of fractures is judged first with the crack angle discriminate as shown in formula (6), is high angle as Y > 0.1
Crack；It is dipping fracture as 0.1 >=Y > 0；It is low angle crack as Y < 0；
Secondly, using the types of fractures judged, by numerical simulation, confirmation bilaterally crack as shown in formula (7)
Coefficient of analysis in porosity formula, the coefficient of analysis are related with fracture development state；
Finally, the al-lateral resistivity log value and mud resistivity of coefficient of analysis and actual measurement based on acquisition, utilize
Formula (7) calculates fracture porosity.
The present embodiment is by providing a kind of fractured-vuggy reservoir fracture porosity quantitative evaluation method and system, to core sample
Experimental analysis corrects and demarcates dual laterolog data, rejects base portion to influence bilaterally, obtains pure rift portion electricity
Values of resistivity, and then the higher fracture porosity of precision is obtained, and then establish fracture porosity quantitative calculation, it is based on above-mentioned step
Suddenly the fracture porosity computation model established, theoretical foundation is strong, precision is high and easy to operate, can be used for quantitative assessment crack hole
Porosity, it is practical, solve the problems, such as the fine quantitative assessment of fracture-cavity type carbonate Reservoir Fracture porosity in oil-gas exploration.
Using example
A concrete application example is given below in the scheme and its effect of the embodiment of the present invention for ease of understanding.This field
It should be understood to the one skilled in the art that the example is only for the purposes of understanding the present invention, any detail is not intended to be limited in any way
The system present invention.
Fig. 3 shows J function according to an embodiment of the invention and water saturation relational graph.Fig. 4 shows root
According to the result map of one embodiment of the present of invention.
Using fractured-vuggy reservoir fracture porosity quantitative evaluation method provided by the present application to Tahe in this application example
Oil field fractured-vuggy reservoir is handled, and carries out quantitative assessment to the fracture porosity of this area.
The first step, according to the experimental data of core sample in table 1 (since the experimental data of a core sample has 16 rows,
If all core samples all show it is bigger, therefore 1, table has been selected the experimental data of two core samples to show), utilize formula
(1) by core analysis hollow billet Curve transform at J function, in this step, porosityIt is experimental data in table 1, P with permeability Kc
Experimental data is analyzed for capillary pressure, σ cos θ takes 26.
Table 1
Second step draws J function and water saturation relational graph according to the J function that the first step calculates, as shown in figure 3, from
In Fig. 3 J function can be established using formula (2) and satisfied with aqueous to be fitted it is clear to see that J function can be divided into two classes
It is expressed as with degree relationship:
Third step calculates continuous depth J function using formula (1), at this time formula (1) porosityIt is sharp in this example
Calculated with neutron density crossplot method, obtain the porosity data of full well section, K according to K withRegression equation seek, connected
After continuous depth J function, continuous depth water saturation then is calculated using formula (9).
4th step utilizes the Archie formula inverse base portion resistivity R as shown in formula (3)tb；According to this area's rock
Electricity experiment, in formula a, b, m, n take 2.1472 respectively, 1.10137,1.3452,3.3506, Rw take 0.02 Ω m.
5th step calculates pure crack rock bilaterally using the series-parallel electrical conduction model as shown in formula (4) and formula (5)
Resistivity RdfWith Rsf。
6th step, the pure crack rock al-lateral resistivity R obtained using step 5dfWith Rsf, utilize formula (6) and table 2
Confirming the coefficient of analysis in formula (7), the parameter in table 2 is obtained by three-dimensional finite element simulation, specific as follows shown:
Table 2
| a1 | a2 | a3 | |
| Y>0.1 | -1711.5 | 1625.4 | 0.1153 |
| 0.1≥Y>0 | 2597.1 | -2348.4 | -0.3135 |
| Y<0 | 430.4707 | -327.3297 | -0.0131 |
The fracture porosity of this area can be calculated using formula (7)As a result as shown in figure 4, it is second from the bottom in figure
Road (Jfunction) is the J function calculated using formula (1), and road dotted line last is to be calculated using Image Logging Data
Fracture porosity, solid line is the fracture porosity calculated using the present invention, on the whole it can be seen that the two numerical values recited and trend
It is almost the same.
This application example is by providing a kind of fractured-vuggy reservoir fracture porosity quantitative evaluation method and system, to core sample
Product experimental analysis corrects and demarcates dual laterolog data, rejects base portion to influence bilaterally, obtains pure rift portion
Resistivity value, and then the higher fracture porosity of precision is obtained, and then establish fracture porosity quantitative calculation, based on above-mentioned
The fracture porosity computation model that step is established, theoretical foundation is strong, precision is high and easy to operate, can be used for quantitative assessment crack
Porosity, it is practical, solve the problems, such as the fine quantitative assessment of fracture-cavity type carbonate Reservoir Fracture porosity in oil-gas exploration.
It will be understood by those skilled in the art that above to the purpose of the description of the embodiment of the present invention only for illustratively saying
The beneficial effect of bright the embodiment of the present invention is not intended to limit embodiments of the invention to given any example.
Embodiment 2
According to an embodiment of the invention, providing a kind of fractured-vuggy reservoir fracture porosity Quantitative Evaluation System, system packet
Include: memory is stored with computer executable instructions；Processor, when calculating executable instruction on run memory, processor
It performs the steps of
1) it is handled using rock core capillary pressure data of the J function to continuous depth, obtains continuous depth J functional value；
2) relational expression based on J function and water saturation obtains continuous depth water saturation value；
3) it is based on continuous depth water saturation value, calculates base portion resistivity；
4) it is based on series-parallel electrical conduction model, calculates pure crack rock al-lateral resistivity；
5) crack angle discriminate and bilaterally fracture porosity formula are utilized, fracture porosity is calculated.
The embodiment is by providing a kind of fractured-vuggy reservoir fracture porosity Quantitative Evaluation System, it is therefore an objective to by rock core
Sample experiments analysis, corrects and demarcates dual laterolog data, rejects base portion to influence bilaterally, obtains pure crack portion
Sub-resistance rate value, and then the higher fracture porosity of precision is obtained, it finally realizes and utilizes dual laterolog data counterincision slot apertures
Porosity carries out fast and accurately quantitative assessment.
In one example, J function representation are as follows:
Wherein, PcFor hollow billet curve capillary pressure value；σ is interfacial tension；θ is angle of wetting；K is permeability；For hole
Degree.
The present embodiment is by providing a kind of fractured-vuggy reservoir fracture porosity Quantitative Evaluation System, to core sample experiment point
Analysis corrects and demarcates dual laterolog data, rejects base portion to influence bilaterally, obtains pure rift portion resistivity
Value, and then the higher fracture porosity of precision is obtained, and then establish fracture porosity quantitative calculation, finally realize utilization
Dual laterolog data fracture porosity carries out fast and accurately quantitative assessment, solves fracture hole type carbonate in oil-gas exploration
The fine quantitative assessment problem of rock Reservoir Fracture porosity.
It will be understood by those skilled in the art that above to the purpose of the description of the embodiment of the present invention only for illustratively saying
The beneficial effect of bright the embodiment of the present invention is not intended to limit embodiments of the invention to given any example.
Various embodiments of the present invention are described above, above description is exemplary, and non-exclusive, and
It is not limited to disclosed each embodiment.Without departing from the scope and spirit of illustrated each embodiment, for this skill
Many modifications and changes are obvious for the those of ordinary skill in art field.
Claims (10)
1. a kind of fractured-vuggy reservoir fracture porosity quantitative evaluation method, which is characterized in that this method comprises:
1) it is handled using rock core capillary pressure data of the J function to continuous depth, obtains continuous depth J functional value；
2) relational expression based on the J function and water saturation, obtains continuous depth water saturation value；
3) it is based on the continuous depth water saturation value, calculates base portion resistivity；
4) it is based on series-parallel electrical conduction model, calculates pure crack rock al-lateral resistivity；
5) crack angle discriminate and bilaterally fracture porosity formula are utilized, fracture porosity is calculated.
2. fractured-vuggy reservoir fracture porosity quantitative evaluation method according to claim 1, wherein the J function representation
Are as follows:
Wherein, PcFor hollow billet curve capillary pressure value；σ is interfacial tension；θ is angle of wetting；K is permeability；For porosity.
3. fractured-vuggy reservoir fracture porosity quantitative evaluation method according to claim 1, wherein the J function with contain
The relational expression of water saturation indicates are as follows:
Sw=cedJ (2)
Wherein, SwFor water saturation；C, d is constant；J is J function.
4. fractured-vuggy reservoir fracture porosity quantitative evaluation method according to claim 1, wherein in step 3), utilize
Archie formula calculates base portion resistivity, and the Archie formula indicates are as follows:
Wherein, RtbFor base portion resistivity；A is lithology factor；B is lithology dependent constant；RwFor formation water resistivity；SwFor
Water saturation；N is saturation exponent；For porosity；M is lithologic index.
5. fractured-vuggy reservoir fracture porosity quantitative evaluation method according to claim 1, wherein the series-parallel conduction
Model is expressed as:
Wherein, RdFor deep lateral resistivity log value；RdfFor pure crack rock deep lateral resistivity；RsIt is surveyed for shallow lateral resistivity
Well value；RsfFor pure crack rock shallowly lateral resistivity；RtbFor base portion resistivity；For porosity.
6. fractured-vuggy reservoir fracture porosity quantitative evaluation method according to claim 1, wherein step 5) includes:
5.1) the pure crack rock al-lateral resistivity obtained based on step 4) obtains the class in crack using crack angle discriminate
Type；
5.2) type in the crack judged based on step 5.1), the analysis system of the confirmation bilaterally fracture porosity formula
Number；
5.3) coefficient of analysis obtained based on step 5.2) obtains fracture pore using the bilaterally fracture porosity formula
Degree.
7. fractured-vuggy reservoir fracture porosity quantitative evaluation method according to claim 1, wherein the crack angle is sentenced
Other formula indicates are as follows:
Wherein, Y is discriminant index；RdfFor pure crack rock deep lateral resistivity；RsfFor pure crack rock shallowly lateral resistivity.
8. fractured-vuggy reservoir fracture porosity quantitative evaluation method according to claim 1, wherein the bilaterally crack
Porosity formula indicates are as follows:
Wherein,For fracture porosity；a1、a2And a3For coefficient of analysis；RdFor deep lateral resistivity；RsFor shallow lateral resistivity；
RmfFor mud resistivity.
9. a kind of fractured-vuggy reservoir fracture porosity Quantitative Evaluation System, which is characterized in that the system comprises:
Memory is stored with computer executable instructions；
Processor, when running the calculating executable instruction on the memory, the processor is performed the steps of
1) it is handled using rock core capillary pressure data of the J function to continuous depth, obtains continuous depth J functional value；
2) relational expression based on the J function and water saturation, obtains continuous depth water saturation value；
3) it is based on the continuous depth water saturation value, calculates base portion resistivity；
4) it is based on series-parallel electrical conduction model, calculates pure crack rock al-lateral resistivity；
5) crack angle discriminate and bilaterally fracture porosity formula are utilized, fracture porosity is calculated.
10. fractured-vuggy reservoir fracture porosity Quantitative Evaluation System according to claim 9, wherein the J function representation
Are as follows:
Wherein, PcFor hollow billet curve capillary pressure value；σ is interfacial tension；θ is angle of wetting；K is permeability；For porosity.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810286736.8A CN110320139B (en) | 2018-03-30 | 2018-03-30 | Fracture-cavity reservoir fracture porosity quantitative evaluation method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810286736.8A CN110320139B (en) | 2018-03-30 | 2018-03-30 | Fracture-cavity reservoir fracture porosity quantitative evaluation method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110320139A true CN110320139A (en) | 2019-10-11 |
| CN110320139B CN110320139B (en) | 2021-12-24 |
Family
ID=68112121
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810286736.8A Active CN110320139B (en) | 2018-03-30 | 2018-03-30 | Fracture-cavity reservoir fracture porosity quantitative evaluation method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110320139B (en) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110967771A (en) * | 2019-12-19 | 2020-04-07 | 西南石油大学 | Method for calculating porosity of cracks and corrosion holes by utilizing conventional logging |
| CN111751905A (en) * | 2020-05-11 | 2020-10-09 | 中国石油天然气股份有限公司 | Carbonate reservoir effectiveness identification method based on rock three-porosity model |
| CN112485174A (en) * | 2020-10-19 | 2021-03-12 | 中国地质大学（北京） | Method for calculating permeability of reservoir containing hydrate based on stacked cube model |
| CN113029908A (en) * | 2021-03-16 | 2021-06-25 | 中国石油大学(华东) | Laboratory measurement method for compact reservoir saturation index |
| CN113759440A (en) * | 2021-08-18 | 2021-12-07 | 西南石油大学 | Method for correcting Biot coefficient of formation rock |
| CN113989433A (en) * | 2021-10-26 | 2022-01-28 | 重庆科技学院 | Fracture-cavity reservoir saturation model building method based on pore type subdivision |
| CN114076727A (en) * | 2022-01-10 | 2022-02-22 | 中国空气动力研究与发展中心低速空气动力研究所 | Resistivity-based ice porosity measurement method |
| CN114185108A (en) * | 2020-09-14 | 2022-03-15 | 中国石油化工股份有限公司 | Crack porosity determination method, storage medium, and computer device |
| CN114429044A (en) * | 2022-01-21 | 2022-05-03 | 成都理工大学 | Method for establishing fracture-cavity reservoir conductivity model and application |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103422853A (en) * | 2012-05-24 | 2013-12-04 | 中国石油天然气集团公司 | Method and device for determining water saturation of stratum |
| CN104847340A (en) * | 2015-03-26 | 2015-08-19 | 中国海洋石油总公司 | Flooded-layer well-logging quantitative evaluation method |
| CN105781539A (en) * | 2016-03-15 | 2016-07-20 | 中国石油大学(华东) | Saturability well logging calculation method of tight oil and gas reservoir |
| CN106680172A (en) * | 2016-11-16 | 2017-05-17 | 中国石油大学(北京) | Method for evaluating fractures of compact oil reservoirs |
| CN107301255A (en) * | 2016-04-15 | 2017-10-27 | 中国石油化工股份有限公司 | A kind of method for calculating Reservoir Fracture porosity |
-
2018
- 2018-03-30 CN CN201810286736.8A patent/CN110320139B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103422853A (en) * | 2012-05-24 | 2013-12-04 | 中国石油天然气集团公司 | Method and device for determining water saturation of stratum |
| CN104847340A (en) * | 2015-03-26 | 2015-08-19 | 中国海洋石油总公司 | Flooded-layer well-logging quantitative evaluation method |
| CN105781539A (en) * | 2016-03-15 | 2016-07-20 | 中国石油大学(华东) | Saturability well logging calculation method of tight oil and gas reservoir |
| CN107301255A (en) * | 2016-04-15 | 2017-10-27 | 中国石油化工股份有限公司 | A kind of method for calculating Reservoir Fracture porosity |
| CN106680172A (en) * | 2016-11-16 | 2017-05-17 | 中国石油大学(北京) | Method for evaluating fractures of compact oil reservoirs |
Non-Patent Citations (2)
| Title |
|---|
| 司马立强 等: "用J函数法求取碳酸盐岩储层饱和度方法探讨", 《岩性油气藏》 * |
| 李善军 等: "碳酸盐岩地层中裂缝孔隙度的定量解释", 《测井技术》 * |
Cited By (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110967771B (en) * | 2019-12-19 | 2021-04-09 | 西南石油大学 | Method for calculating porosity of cracks and corrosion holes by utilizing conventional logging |
| CN110967771A (en) * | 2019-12-19 | 2020-04-07 | 西南石油大学 | Method for calculating porosity of cracks and corrosion holes by utilizing conventional logging |
| CN111751905A (en) * | 2020-05-11 | 2020-10-09 | 中国石油天然气股份有限公司 | Carbonate reservoir effectiveness identification method based on rock three-porosity model |
| CN114185108A (en) * | 2020-09-14 | 2022-03-15 | 中国石油化工股份有限公司 | Crack porosity determination method, storage medium, and computer device |
| CN114185108B (en) * | 2020-09-14 | 2024-04-12 | 中国石油化工股份有限公司 | Crack porosity determination method, storage medium, and computer device |
| CN112485174A (en) * | 2020-10-19 | 2021-03-12 | 中国地质大学（北京） | Method for calculating permeability of reservoir containing hydrate based on stacked cube model |
| CN112485174B (en) * | 2020-10-19 | 2021-09-14 | 中国地质大学（北京） | Method for calculating permeability of reservoir containing hydrate based on stacked cube model |
| CN113029908B (en) * | 2021-03-16 | 2021-11-26 | 中国石油大学(华东) | Laboratory measurement method for compact reservoir saturation index |
| CN113029908A (en) * | 2021-03-16 | 2021-06-25 | 中国石油大学(华东) | Laboratory measurement method for compact reservoir saturation index |
| CN113759440A (en) * | 2021-08-18 | 2021-12-07 | 西南石油大学 | Method for correcting Biot coefficient of formation rock |
| CN113759440B (en) * | 2021-08-18 | 2022-08-19 | 西南石油大学 | Method for correcting Biot coefficient of formation rock |
| CN113989433A (en) * | 2021-10-26 | 2022-01-28 | 重庆科技学院 | Fracture-cavity reservoir saturation model building method based on pore type subdivision |
| CN114076727A (en) * | 2022-01-10 | 2022-02-22 | 中国空气动力研究与发展中心低速空气动力研究所 | Resistivity-based ice porosity measurement method |
| CN114076727B (en) * | 2022-01-10 | 2022-05-13 | 中国空气动力研究与发展中心低速空气动力研究所 | Resistivity-based ice porosity measurement method |
| CN114429044A (en) * | 2022-01-21 | 2022-05-03 | 成都理工大学 | Method for establishing fracture-cavity reservoir conductivity model and application |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110320139B (en) | 2021-12-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110320139A (en) | Fractured-vuggy reservoir fracture porosity quantitative evaluation method and system | |
| EP2965126B1 (en) | Determining continuous capillary pressure curves for subsurface earth formations using saturation and nmr log data | |
| US7532983B2 (en) | Method and apparatus for measuring the wettability of geological formations | |
| Nazemi et al. | The effect of carbonate reservoir heterogeneity on Archie's exponents (a and m), an example from Kangan and Dalan gas formations in the central Persian Gulf | |
| CN106154343B (en) | The method for calculating the oil saturation of fine and close oily reservoir | |
| CN103630939A (en) | Air layer identification and evaluation method | |
| Liang et al. | Tight gas sandstone reservoirs evaluation from nuclear magnetic resonance (NMR) logs: case studies | |
| EP2758810A1 (en) | Integrated dip from 3d resistivity tool and borehole imaging tool | |
| CN111381292B (en) | Logging interpretation method and device for predicting sandstone hydrocarbon-bearing reservoir | |
| Ding et al. | Reliability analysis of seismic attribute in the detection of fault-karst | |
| Al-Sulami et al. | The unconventional shale reservoirs of jafurah basin: An integrated petrophysical evaluation using cores and advanced well logs | |
| Lopez et al. | Petrophysical quantification of multiple porosities in shale petroleum reservoirs | |
| Xiao et al. | Tight-gas-sand permeability estimation from nuclear-magnetic-resonance (NMR) logs based on the hydraulic-flow-unit (HFU) approach | |
| Garcia et al. | Integrated workflow to estimate permeability through quantification of rock fabric using joint interpretation of nuclear magnetic resonance and electric measurements | |
| CN110658555A (en) | Reservoir fluid factor construction method and reservoir fluid identification method | |
| CN108956410A (en) | The recognition methods of rock sample pore structure fast qualitative and device | |
| Spain et al. | Petrophysical rock typing of tight gas sands-beyond porosity and saturation: example from the cotton valley formation, East Texas | |
| CN106353813B (en) | Flow net model method based on array sonic log | |
| Deng et al. | A new index used to characterize the near-wellbore fracture network in naturally fractured gas reservoirs | |
| CN109826623A (en) | Knowledge method is sentenced in a kind of geophysical log of tight sandstone reservoir stratification seam | |
| CN109709610A (en) | A kind of rock fracture detection method and system | |
| Shen et al. | Study on fluid mobility in sandwich-type shale oil reservoir using two-dimensional nuclear magnetic resonance approaches | |
| Nooh et al. | Comparison of quantitative analysis of image logs for shale volume and net to gross calculation of a thinly laminated reservoir between VNG-NERGE and LAGIA-EGYPT | |
| Newgord et al. | Wettability Assessment in Complex Formations Using NMR Measurements: Workflow Development and Experimental Core-Scale Verificiation | |
| McCarney et al. | Core plug nuclear magnetic resonance (NMR) analysis as a method to estimate permeability anisotropy |
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 |