CN110320139B - 缝洞型储层裂缝孔隙度定量评价方法及*** - Google Patents
缝洞型储层裂缝孔隙度定量评价方法及*** Download PDFInfo
- Publication number
- CN110320139B CN110320139B CN201810286736.8A CN201810286736A CN110320139B CN 110320139 B CN110320139 B CN 110320139B CN 201810286736 A CN201810286736 A CN 201810286736A CN 110320139 B CN110320139 B CN 110320139B
- Authority
- CN
- China
- Prior art keywords
- fracture
- porosity
- resistivity
- lateral
- pure
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000011158 quantitative evaluation Methods 0.000 title claims abstract description 39
- 238000000034 method Methods 0.000 title claims abstract description 38
- 239000011435 rock Substances 0.000 claims abstract description 45
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 42
- 239000011159 matrix material Substances 0.000 claims abstract description 28
- 238000012545 processing Methods 0.000 claims abstract description 9
- 230000006870 function Effects 0.000 claims description 42
- 238000004458 analytical method Methods 0.000 claims description 23
- 230000002146 bilateral effect Effects 0.000 claims description 13
- 230000035699 permeability Effects 0.000 claims description 9
- 238000009736 wetting Methods 0.000 claims description 6
- 239000008398 formation water Substances 0.000 claims description 4
- 206010017076 Fracture Diseases 0.000 description 123
- 208000010392 Bone Fractures Diseases 0.000 description 118
- 238000004364 calculation method Methods 0.000 description 13
- 230000009977 dual effect Effects 0.000 description 11
- 230000004044 response Effects 0.000 description 9
- 238000011160 research Methods 0.000 description 7
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 6
- 230000008901 benefit Effects 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 4
- 239000007789 gas Substances 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- 239000011148 porous material Substances 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 235000010363 Vitex negundo Nutrition 0.000 description 1
- 244000248021 Vitex negundo Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000012625 in-situ measurement Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 238000004457 water analysis Methods 0.000 description 1
Images
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
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开了一种缝洞型储层裂缝孔隙度定量评价方法及***，包括：利用J函数对连续深度的岩心毛管压力数据进行处理，获得连续深度J函数值；基于J函数与含水饱和度的关系式，获得连续深度含水饱和度值；基于连续深度含水饱和度值，计算基质部分电阻率；基于串并联导电模型，计算纯裂缝岩石双侧向电阻率；利用裂缝角度判别式和双侧向裂缝孔隙度公式，计算裂缝孔隙度。本发明提供的缝洞型储层裂缝孔隙度定量评价方法及***，能够对裂缝孔隙度进行快速准确的定量评价。
Description
技术领域
本发明属于石油与地质领域，更具体地，涉及一种缝洞型储层裂缝孔隙度定量评价方法及***。
背景技术
在油气资源勘探中，碳酸盐岩缝洞型储层占有重要的地位，裂缝对储层产能有着积极重要的控制作用，作为评价该类储层的重要参数之一，裂缝孔隙度的求取和标定一直是世界性难题。
Sibbit A.和Faiver O.(Sibbit A,Faivre O.The dual laterolog response infractured rocks[C].SPWLA 26th Annual Logging Symposium.Dallas,Texas,1985,17-20.)利用三维有限元数值计算方法，建立了垂直裂缝和水平裂缝开度计算模型；罗贞耀(罗贞耀.用侧向资料计算裂缝张开度的初步研究[J].地球物理测井,1990,14(2):83-92.)基于双侧向测井资料建立了任意视角度裂缝开度模型；随后Philippe A.P.和Roger N.A.(Philippe A.P.,Roger N.A.In situ measurements of electrical resistivity,formation anisotropy and tectonic context[C].SPWLA 31st Annual LoggingSymposium,Lafayette,Louisiana,1990,24-27.)正演导出了双侧向在任意倾角情况下裂缝的测井响应，并得到了准水平缝与准立缝条件下双侧向响应值与裂缝张开度的关系；通过双侧向深浅视电阻率差异来判断裂缝角度，但该模型的建立是在较为特殊的条件下进行的，对于中间角度裂缝的响应规律没有更进一步总结；在Sibbit A.和Philippe A.P.的研究基础上，李善军(李善军,肖承文,汪涵明,等.裂缝的双侧向测井响应的数学模型及裂缝孔隙度的定量解释[J].地球物理学报,1996,39(6):845-852.)等人利用三维有限元数值计算方法以及电导率张量的方法分别模拟了不同倾角情况下，单一裂缝、等间距平行裂缝、交叉缝及宏观各向异性裂缝地层的双侧向测井响应，并且归纳了各种情况裂缝的测井响应值与裂缝倾角的关系，将裂缝进一步划分为三类，即高角度缝，倾斜裂缝和低角度缝，并具体给出了各种情况下反演裂缝孔隙度与裂缝倾角的方法；邓少贵等(邓少贵,卢春利.裂缝-孔隙双孔介质的双侧向测井解释方法初探[J].测井技术,2010,34(6),537-541.)建立了双侧向测井响应与基块电导率、裂缝孔隙度、倾角的关系；秦启荣(秦启荣，黄平辉，周志远，等.全直径样品分析在测井解释裂缝孔隙度中的应用[J].天然气地球科学，2005,16(5):637-640)等利用全直径样品进行裂缝孔隙度计算；吕洪志(吕洪志，陆云龙，崔云江，等.改进的孔隙度模型评价流体性质与裂缝孔隙度[J].应用声学，2016,35(4):351-356)依据Biot理论，利用阵列声波测井数据进行反演得到裂缝孔隙度；张莹(张莹，潘保芝.基于FMI的3种火山岩储层裂缝孔隙度求取方法[J].测井技术，2012,36(4):365-369)利用FMI资料分析了三种火山岩裂缝孔隙度求取方法；尚根华(尚根华，候晓春，刘学伟，等.利用X-CT研究大庆油田双重介质裂缝和砂岩孔隙度[J].CT理论与应用研究，1997,6(1):26-31)利用X-CT研究了双重介质裂缝和砂岩孔隙度；刘利波(刘利波.三孔隙度测井方法在识别裂缝中的应用[内蒙古石油化工].2014,11:155-156)通过对三孔隙度进行一致化处理，引入裂缝指示特征系数，达到对裂缝识别及孔隙度计算；秦海旭(秦海旭，吴国枕.裂缝储层孔隙度、饱和度计算方法[J].石油地球物理勘探，2015,50(1):103-110.)运用裂缝介质各向异性梯度计算裂缝介质孔隙度和饱和度的新方法。
综合分析前人已有的成果，针对裂缝孔隙度的评价主要集中在五类方法上，1.基于双侧向测井资料的计算方法，该方法缺点在于待定参数众多，且双侧向的响应处理受裂缝影响外，还受孔隙中流体影响，多解性强；2.基于电阻率成像测井资料方法，该方法局限性在于费用昂贵，并不是所有井都进行成像测量，同时电成像计算的裂缝仅是图像上见到的裂缝，图像上看不到的裂缝并不参与计算；3.基于三孔隙度方法，该方法难点在于岩石骨架参数的准确选取，且裂缝高角度时，并不适用；4.基于全直径岩心分析法，该方法缺点在于，裂缝段取芯成功率难以保证，且需要对岩心进行压实、归位等校正；5.岩心X-CT扫描法，该方法计算量大，且代表性岩性样品获取困难。
因此，有必要提供一种缝洞型储层裂缝孔隙度定量评价方法及***，解决油气勘探中缝洞型碳酸盐岩储层裂缝孔隙度精细定量评价问题。
发明内容
本发明通过提供一种缝洞型储层裂缝孔隙度定量评价方法及***，目的是通过对岩心样品实验分析，校正并标定双侧向测井资料，剔除基质部分对双侧向的影响，得到纯裂缝部分电阻率值，进而获取精度较高的裂缝孔隙度，最终实现了利用双侧向测井资料对裂缝孔隙度进行快速准确的定量评价。
根据本发明的一方面，提出了一种缝洞型储层裂缝孔隙度定量评价方法，该方法包括：
1)利用J函数对连续深度的岩心毛管压力数据进行处理，获得连续深度J函数值；
2)基于所述J函数与含水饱和度的关系式，获得连续深度含水饱和度值；
3)基于所述连续深度含水饱和度值，计算基质部分电阻率；
4)基于串并联导电模型，计算纯裂缝岩石双侧向电阻率；
5)利用裂缝角度判别式和双侧向裂缝孔隙度公式，计算裂缝孔隙度。
优选地，所述J函数表示为：
其中，Pc为毛管曲线毛管压力值；σ为界面张力；θ为湿润角；K为渗透率；
为孔隙度。
优选地，所述J函数与含水饱和度的关系式表示为：
Sw＝cedJ (2)
其中，Sw为含水饱和度；c、d为常数；J为J函数。
优选地，步骤3)中，利用阿尔奇公式计算基质部分电阻率，所述阿尔奇公式表示为：
其中，Rtb为基质部分电阻率；a为岩性系数；b为岩性相关常数；Rw为地层水电阻率；Sw为含水饱和度；n为饱和度指数；
为孔隙度；m为岩性指数。
优选地，所述串并联导电模型表示为：
其中，Rd为深侧向电阻率测井值；Rdf为纯裂缝岩石深侧向电阻率；Rs为浅侧向电阻率测井值；Rsf为纯裂缝岩石浅侧向电阻率；Rtb为基质部分电阻率；
为孔隙度。
优选地，步骤5)包括：
5.1)基于步骤4)获得的纯裂缝岩石双侧向电阻率，利用裂缝角度判别式获得裂缝的类型；
5.2)基于步骤5.1)判断出的裂缝的类型，确认所述双侧向裂缝孔隙度公式的分析系数；
5.3)基于步骤5.2)获得的分析系数，利用所述双侧向裂缝孔隙度公式获得裂缝孔隙度。
优选地，所述裂缝角度判别式表示为：
其中，Y为判别指数；Rdf为纯裂缝岩石深侧向电阻率；Rsf为纯裂缝岩石浅侧向电阻率。
优选地，所述双侧向裂缝孔隙度公式表示为：
其中，
为裂缝孔隙度；a1、a2和a3为分析系数；Rd为深侧向电阻率；Rs为浅侧向电阻率；Rmf为泥浆电阻率。
根据本发明的另一方面，提出了一种缝洞型储层裂缝孔隙度定量评价***，所述***包括：存储器，存储有计算机可执行指令；处理器，运行所述存储器上的计算可执行指令时，所述处理器实现以下步骤：
1)利用J函数对连续深度的岩心毛管压力数据进行处理，获得连续深度J函数值；
2)基于所述J函数与含水饱和度的关系式，获得连续深度含水饱和度值；
3)基于所述连续深度含水饱和度值，计算基质部分电阻率；
4)基于串并联导电模型，计算纯裂缝岩石双侧向电阻率；
5)利用裂缝角度判别式和双侧向裂缝孔隙度公式，计算裂缝孔隙度。
优选地，所述J函数表示为：
其中，Pc为毛管曲线毛管压力值；σ为界面张力；θ为湿润角；K为渗透率；
为孔隙度。
本发明的有益效果在于：通过提供一种缝洞型储层裂缝孔隙度定量评价方法及***，对岩心样品实验分析，校正并标定双侧向测井资料，剔除基质部分对双侧向的影响，得到纯裂缝部分电阻率值，进而获取精度较高的裂缝孔隙度，进而建立裂缝孔隙度定量计算模型，最终实现了利用双侧向测井资料对裂缝孔隙度进行快速准确的定量评价，解决了油气勘探中缝洞型碳酸盐岩储层裂缝孔隙度精细定量评价问题。
本发明的其它特征和优点将在随后具体实施方式部分予以详细说明。
附图说明
通过结合附图对本发明示例性实施方式进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施方式中，相同的参考标号通常代表相同部件。
图1示出了根据本发明的缝洞型储层裂缝孔隙度定量评价方法的流程图。
图2示出了根据本发明的缝洞型储层的示意图。
图3示出了根据本发明的一个实施例的J函数与含水饱和度关系图。
图4示出了根据本发明的一个实施例的成果图。
具体实施方式
下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。
实施例1
在该实施例中，根据本发明的缝洞型储层裂缝孔隙度定量评价方法可以包括：
1)利用J函数对连续深度的岩心毛管压力数据进行处理，获得连续深度J函数值；
2)基于J函数与含水饱和度的关系式，获得连续深度含水饱和度值；
3)基于连续深度含水饱和度值，计算基质部分电阻率；
4)基于串并联导电模型，计算纯裂缝岩石双侧向电阻率；
5)利用裂缝角度判别式和双侧向裂缝孔隙度公式，计算裂缝孔隙度。
图2示出了根据本发明的缝洞型储层的示意图。如图2所示，缝洞型碳酸盐岩储层由基质部分和裂缝构成，因此，基质部分孔隙与裂缝都对双侧向有分流作用，都会影响双侧向电阻率。该实施例通过提供一种缝洞型储层裂缝孔隙度定量评价方法及***，通过对岩心样品实验分析，校正并标定双侧向测井资料，剔除基质部分对双侧向的影响，得到纯裂缝部分电阻率值，进而获取精度较高的裂缝孔隙度，最终实现了利用双侧向测井资料对裂缝孔隙度进行快速准确的定量评价。
图1示出了根据本发明的缝洞型储层裂缝孔隙度定量评价方法的流程图。下面参考图1详细说明根据本发明的缝洞型储层裂缝孔隙度定量评价方法的具体步骤。
步骤101，利用J函数对连续深度的岩心毛管压力数据进行处理，获得连续深度J函数值。
具体地，毛管压力数据是实验室根据某几个深度点测量的，J函数处理的目的就是找到能够表征整个井段的一个关系式。
在一个示例中，J函数表示为：
其中，Pc为毛管曲线毛管压力值；σ为界面张力；θ为湿润角；K为渗透率；
为孔隙度。
具体地，孔隙度
可根据中子密度交会求取；渗透率K可通过岩心孔隙度与渗透率回归得到计算模型；σ与θ根据地层流体情况确定，当为油水接触时，σcosθ取26，当为气水接触时，σcosθ取50。
另外，毛管曲线毛管压力值Pc表示为：
Pc＝0.01(ρw-ρ0)H (8)
其中，ρw为油的密度；ρ0为水的密度；H为油藏高度。
步骤102，基于J函数与含水饱和度的关系式，获得连续深度含水饱和度值。
在一个示例中，J函数与含水饱和度的关系式表示为：
Sw＝cedJ (2)
其中，Sw为含水饱和度；c、d为常数；J为J函数。
具体地，c、d为常数，是无因次量，具体数值根据研究区实验数据获取，结合研究区实验数据，利用公式(1)、公式(2)和公式(8)，即可获得连续深度含水饱和度值。
步骤103，基于连续深度含水饱和度值，计算基质部分电阻率。
在一个示例中，步骤103中，利用阿尔奇公式计算基质部分电阻率，阿尔奇公式表示为：
其中，Rtb为基质部分电阻率；a为岩性系数；b为岩性相关常数；Rw为地层水电阻率；Sw为含水饱和度；n为饱和度指数；
为孔隙度；m为岩性指数。
具体地，a、b、m和n均为无因次量，可通过岩电实验获取；再通过水分析资料获取地层水电阻率Rw，利用公式(3)即可计算出基质部分电阻率。
步骤104，基于串并联导电模型，计算纯裂缝岩石双侧向电阻率。
在一个示例中，串并联导电模型表示为：
其中，Rd为深侧向电阻率测井值；Rdf为纯裂缝岩石深侧向电阻率；Rs为浅侧向电阻率测井值；Rsf为纯裂缝岩石浅侧向电阻率；Rtb为基质部分电阻率；
为孔隙度。
具体地，根据双侧向电阻率测井值及步骤103计算的基质部分电阻率，利用如公式(4)和公式(5)所示的串并联导电模型，即可得到由纯裂缝引起的双侧向电阻率。
步骤105，利用裂缝角度判别式和双侧向裂缝孔隙度公式，计算裂缝孔隙度。
在一个示例中，步骤105包括：
基于纯裂缝岩石双侧向电阻率，利用裂缝角度判别式获得裂缝的类型；
基于判断出的裂缝的类型，确认双侧向裂缝孔隙度公式的分析系数；
基于获得的分析系数，利用双侧向裂缝孔隙度公式获得裂缝孔隙度。
在一个示例中，裂缝角度判别式表示为：
其中，Y为判别指数；Rdf为纯裂缝岩石深侧向电阻率；Rsf为纯裂缝岩石浅侧向电阻率。
在一个示例中，双侧向裂缝孔隙度公式表示为：
其中，
为裂缝孔隙度；a1、a2和a3为分析系数；Rd为深侧向电阻率；Rs为浅侧向电阻率；Rmf为泥浆电阻率。
具体地，利用如公式(6)所示的裂缝角度判别式和如公式(7)所示的双侧向裂缝孔隙度公式，计算裂缝孔隙度包括：
首先利用如公式(6)所示的裂缝角度判别式判断裂缝类型，当Y＞0.1时，为高角度裂缝；当0.1≥Y＞0时，为倾斜裂缝；当Y＜0时，为低角度裂缝；
其次，利用判断出的裂缝类型，通过数值模拟，确认如公式(7)所示的双侧向裂缝孔隙度公式中的分析系数，该分析系数与裂缝发育状态有关；
最后，基于获得的分析系数和实测的双侧向电阻率测井值以及泥浆电阻率，利用公式(7)计算裂缝孔隙度。
本实施例通过提供一种缝洞型储层裂缝孔隙度定量评价方法及***，对岩心样品实验分析，校正并标定双侧向测井资料，剔除基质部分对双侧向的影响，得到纯裂缝部分电阻率值，进而获取精度较高的裂缝孔隙度，进而建立裂缝孔隙度定量计算模型，基于上述步骤建立的裂缝孔隙度计算模型，理论依据强、精度高并且操作简单，能用于定量评价裂缝孔隙度，实用性强，解决了油气勘探中缝洞型碳酸盐岩储层裂缝孔隙度精细定量评价问题。
应用示例
为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
图3示出了根据本发明的一个实施例的J函数与含水饱和度关系图。图4示出了根据本发明的一个实施例的成果图。
在本应用示例中利用本申请提供的缝洞型储层裂缝孔隙度定量评价方法对塔河油田缝洞型储层进行处理，对该地区的裂缝孔隙度进行定量评价。
第一步，根据表1中岩心样品的实验数据(由于一个岩心样品的实验数据有16行，若所有岩心样品都展示比较大，因此表1只选了两个岩心样品的实验数据展示)，利用公式(1)将岩心分析毛管曲线转换成J函数，此步骤中，孔隙度
和渗透率K为表1中实验数据，Pc为毛管压力分析实验数据，σcosθ取26。
表1
第二步，根据第一步计算的J函数，绘制J函数与含水饱和度关系图，如图3所示，从图3中可以很明显看出，可以将J函数分为两类进行拟合，利用公式(2)建立J函数与含水饱和度关系表示为：
第三步，利用公式(1)计算连续深度J函数，此时公式(1)中孔隙度
在本实例中利用中子密度交会方法计算，获得全井段的孔隙度数据，K根据K与
的回归方程求取，得到连续深度J函数后，然后利用公式(9)计算连续深度含水饱和度。
第四步，利用如公式(3)所示的阿尔奇公式反算基质部分电阻率Rtb；根据该地区岩电实验，式中a、b、m、n分别取2.1472、1.10137、1.3452、3.3506，Rw取0.02Ω·m。
第五步，利用如公式(4)和式(5)所示的串并联导电模型，计算纯裂缝岩石双侧向电阻率Rdf与Rsf。
第六步，利用第5步获得的纯裂缝岩石双侧向电阻率Rdf与Rsf，利用公式(6)和表2确认公式(7)中的分析系数，表2中的参数通过三维有限元数值模拟获得，具体如下所示：
表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 |
利用公式(7)即可计算出该地区的裂缝孔隙度
结果如图4所示，图中倒数第二道(Jfunction)为利用公式(1)计算的J函数，倒数第一道虚线为利用成像测井资料计算的裂缝孔隙度，实线为利用本发明计算的裂缝孔隙度，整体上可以看到二者数值大小和趋势基本一致。
本应用示例通过提供一种缝洞型储层裂缝孔隙度定量评价方法及***，对岩心样品实验分析，校正并标定双侧向测井资料，剔除基质部分对双侧向的影响，得到纯裂缝部分电阻率值，进而获取精度较高的裂缝孔隙度，进而建立裂缝孔隙度定量计算模型，基于上述步骤建立的裂缝孔隙度计算模型，理论依据强、精度高并且操作简单，能用于定量评价裂缝孔隙度，实用性强，解决了油气勘探中缝洞型碳酸盐岩储层裂缝孔隙度精细定量评价问题。
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
实施例2
根据本发明的实施例，提供了一种缝洞型储层裂缝孔隙度定量评价***，***包括：存储器，存储有计算机可执行指令；处理器，运行存储器上的计算可执行指令时，处理器实现以下步骤：
1)利用J函数对连续深度的岩心毛管压力数据进行处理，获得连续深度J函数值；
2)基于J函数与含水饱和度的关系式，获得连续深度含水饱和度值；
3)基于连续深度含水饱和度值，计算基质部分电阻率；
4)基于串并联导电模型，计算纯裂缝岩石双侧向电阻率；
5)利用裂缝角度判别式和双侧向裂缝孔隙度公式，计算裂缝孔隙度。
该实施例通过提供一种缝洞型储层裂缝孔隙度定量评价***，目的是通过对岩心样品实验分析，校正并标定双侧向测井资料，剔除基质部分对双侧向的影响，得到纯裂缝部分电阻率值，进而获取精度较高的裂缝孔隙度，最终实现了利用双侧向测井资料对裂缝孔隙度进行快速准确的定量评价。
在一个示例中，J函数表示为：
其中，Pc为毛管曲线毛管压力值；σ为界面张力；θ为湿润角；K为渗透率；
为孔隙度。
本实施例通过提供一种缝洞型储层裂缝孔隙度定量评价***，对岩心样品实验分析，校正并标定双侧向测井资料，剔除基质部分对双侧向的影响，得到纯裂缝部分电阻率值，进而获取精度较高的裂缝孔隙度，进而建立裂缝孔隙度定量计算模型，最终实现了利用双侧向测井资料对裂缝孔隙度进行快速准确的定量评价，解决了油气勘探中缝洞型碳酸盐岩储层裂缝孔隙度精细定量评价问题。
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (8)
1.一种缝洞型储层裂缝孔隙度定量评价方法，其特征在于，该方法包括：
1)利用J函数对连续深度的岩心毛管压力数据进行处理，获得连续深度J函数值；
2)基于所述J函数与含水饱和度的关系式，获得连续深度含水饱和度值；
3)基于所述连续深度含水饱和度值，计算基质部分电阻率；
4)基于串并联导电模型，计算纯裂缝岩石双侧向电阻率；
5)利用裂缝角度判别式和双侧向裂缝孔隙度公式，计算裂缝孔隙度；
其中，所述串并联导电模型表示为：
其中，Rd为深侧向电阻率测井值；Rdf为纯裂缝岩石深侧向电阻率；Rs为浅侧向电阻率测井值；Rsf为纯裂缝岩石浅侧向电阻率；Rtb为基质部分电阻率；
为孔隙度；
其中，所述裂缝角度判别式表示为：
其中，Y为判别指数；Rdf为纯裂缝岩石深侧向电阻率；Rsf为纯裂缝岩石浅侧向电阻率。
2.根据权利要求1所述的缝洞型储层裂缝孔隙度定量评价方法，其中，所述J函数表示为：
其中，Pc为毛管曲线毛管压力值；σ为界面张力；θ为湿润角；K为渗透率；
为孔隙度。
3.根据权利要求1所述的缝洞型储层裂缝孔隙度定量评价方法，其中，所述J函数与含水饱和度的关系式表示为：
Sw＝cedJ (2)
其中，Sw为含水饱和度；c、d为常数；J为J函数。
4.根据权利要求1所述的缝洞型储层裂缝孔隙度定量评价方法，其中，步骤3)中，利用阿尔奇公式计算基质部分电阻率，所述阿尔奇公式表示为：
其中，Rtb为基质部分电阻率；a为岩性系数；b为岩性相关常数；Rw为地层水电阻率；Sw为含水饱和度；n为饱和度指数；
为孔隙度；m为岩性指数。
5.根据权利要求1所述的缝洞型储层裂缝孔隙度定量评价方法，其中，步骤5)包括：
5.1)基于步骤4)获得的纯裂缝岩石双侧向电阻率，利用裂缝角度判别式获得裂缝的类型；
5.2)基于步骤5.1)判断出的裂缝的类型，确认所述双侧向裂缝孔隙度公式的分析系数；
5.3)基于步骤5.2)获得的分析系数，利用所述双侧向裂缝孔隙度公式获得裂缝孔隙度。
6.根据权利要求1所述的缝洞型储层裂缝孔隙度定量评价方法，其中，所述双侧向裂缝孔隙度公式表示为：
其中，
为裂缝孔隙度；a1、a2和a3为分析系数；Rd为深侧向电阻率；Rs为浅侧向电阻率；Rmf为泥浆电阻率。
7.一种缝洞型储层裂缝孔隙度定量评价***，其特征在于，所述***包括：
存储器，存储有计算机可执行指令；
处理器，运行所述存储器上的计算可执行指令时，所述处理器实现以下步骤：
1)利用J函数对连续深度的岩心毛管压力数据进行处理，获得连续深度J函数值；
2)基于所述J函数与含水饱和度的关系式，获得连续深度含水饱和度值；
3)基于所述连续深度含水饱和度值，计算基质部分电阻率；
4)基于串并联导电模型，计算纯裂缝岩石双侧向电阻率；
5)利用裂缝角度判别式和双侧向裂缝孔隙度公式，计算裂缝孔隙度；
其中，所述串并联导电模型表示为：
其中，Rd为深侧向电阻率测井值；Rdf为纯裂缝岩石深侧向电阻率；Rs为浅侧向电阻率测井值；Rsf为纯裂缝岩石浅侧向电阻率；Rtb为基质部分电阻率；
为孔隙度；
其中，所述裂缝角度判别式表示为：
其中，Y为判别指数；Rdf为纯裂缝岩石深侧向电阻率；Rsf为纯裂缝岩石浅侧向电阻率。
8.根据权利要求7所述的缝洞型储层裂缝孔隙度定量评价***，其中，所述J函数表示为：
其中，Pc为毛管曲线毛管压力值；σ为界面张力；θ为湿润角；K为渗透率；
为孔隙度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810286736.8A CN110320139B (zh) | 2018-03-30 | 2018-03-30 | 缝洞型储层裂缝孔隙度定量评价方法及*** |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810286736.8A CN110320139B (zh) | 2018-03-30 | 2018-03-30 | 缝洞型储层裂缝孔隙度定量评价方法及*** |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110320139A CN110320139A (zh) | 2019-10-11 |
| CN110320139B true CN110320139B (zh) | 2021-12-24 |
Family
ID=68112121
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810286736.8A Active CN110320139B (zh) | 2018-03-30 | 2018-03-30 | 缝洞型储层裂缝孔隙度定量评价方法及*** |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110320139B (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110967771B (zh) * | 2019-12-19 | 2021-04-09 | 西南石油大学 | 一种利用常规测井计算裂缝、溶蚀孔洞孔隙度的方法 |
| CN111751905A (zh) * | 2020-05-11 | 2020-10-09 | 中国石油天然气股份有限公司 | 一种基于岩石三孔隙度模型的碳酸盐岩储层有效性识别方法 |
| CN114185108B (zh) * | 2020-09-14 | 2024-04-12 | 中国石油化工股份有限公司 | 裂缝孔隙度确定方法、存储介质和计算机设备 |
| CN112485174B (zh) * | 2020-10-19 | 2021-09-14 | 中国地质大学（北京） | 基于堆叠立方体模型计算含水合物储层渗透率的方法 |
| CN113029908B (zh) * | 2021-03-16 | 2021-11-26 | 中国石油大学(华东) | 一种致密储层饱和度指数的实验室测量方法 |
| CN113759440B (zh) * | 2021-08-18 | 2022-08-19 | 西南石油大学 | 一种地层岩石Biot系数的校正方法 |
| CN113989433B (zh) * | 2021-10-26 | 2022-06-07 | 重庆科技学院 | 一种基于孔隙类型细分的缝洞型储层饱和度模型建立方法 |
| CN114076727B (zh) * | 2022-01-10 | 2022-05-13 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种基于电阻率的冰的孔隙率测量方法 |
| CN114429044B (zh) * | 2022-01-21 | 2023-07-28 | 成都理工大学 | 一种建立缝洞型储层电导率模型的方法及应用 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103422853A (zh) * | 2012-05-24 | 2013-12-04 | 中国石油天然气集团公司 | 一种确定地层含水饱和度的方法及装置 |
| CN104847340A (zh) * | 2015-03-26 | 2015-08-19 | 中国海洋石油总公司 | 一种水淹层测井定量评价方法 |
| CN105781539A (zh) * | 2016-03-15 | 2016-07-20 | 中国石油大学(华东) | 一种致密油气储层饱和度测井计算方法 |
| CN106680172A (zh) * | 2016-11-16 | 2017-05-17 | 中国石油大学(北京) | 评价致密油储层的裂缝的方法 |
| CN107301255A (zh) * | 2016-04-15 | 2017-10-27 | 中国石油化工股份有限公司 | 一种计算储层裂缝孔隙度的方法 |
-
2018
- 2018-03-30 CN CN201810286736.8A patent/CN110320139B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103422853A (zh) * | 2012-05-24 | 2013-12-04 | 中国石油天然气集团公司 | 一种确定地层含水饱和度的方法及装置 |
| CN104847340A (zh) * | 2015-03-26 | 2015-08-19 | 中国海洋石油总公司 | 一种水淹层测井定量评价方法 |
| CN105781539A (zh) * | 2016-03-15 | 2016-07-20 | 中国石油大学(华东) | 一种致密油气储层饱和度测井计算方法 |
| CN107301255A (zh) * | 2016-04-15 | 2017-10-27 | 中国石油化工股份有限公司 | 一种计算储层裂缝孔隙度的方法 |
| CN106680172A (zh) * | 2016-11-16 | 2017-05-17 | 中国石油大学(北京) | 评价致密油储层的裂缝的方法 |
Non-Patent Citations (2)
| Title |
|---|
| 用J函数法求取碳酸盐岩储层饱和度方法探讨;司马立强 等;《岩性油气藏》;20141231;第26卷(第6期);第106-110页 * |
| 碳酸盐岩地层中裂缝孔隙度的定量解释;李善军 等;《测井技术》;19971231;第21卷(第3期);第208页第2段至第212页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110320139A (zh) | 2019-10-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110320139B (zh) | 缝洞型储层裂缝孔隙度定量评价方法及*** | |
| Xiao et al. | A new methodology of constructing pseudo capillary pressure (Pc) curves from nuclear magnetic resonance (NMR) logs | |
| McKinley et al. | How porosity and permeability vary spatially with grain size, sorting, cement volume, and mineral dissolution in fluvial Triassic sandstones: the value of geostatistics and local regression | |
| CN105659074B (zh) | 三轴核磁共振测试仪器 | |
| Verwer et al. | Effect of pore structure on electrical resistivity in carbonates | |
| US7532983B2 (en) | Method and apparatus for measuring the wettability of geological formations | |
| CN105089663B (zh) | 一种大斜度/水平井地层电阻率各向异性校正方法 | |
| Minh et al. | Determination of wettability from magnetic resonance relaxation and diffusion measurements on fresh-state cores | |
| BRPI0903061B1 (pt) | Método para identificação de fraturas induzidas por perfuração e instrumento para identificação de fraturas induzidas por perfuração em uma formação que circunda um furo de poço | |
| MX2007007060A (es) | Metodo para la determinacion de la saturacion de agua de una formacion subterranea. | |
| CN107829731B (zh) | 一种黏土蚀变的火山岩孔隙度校正方法 | |
| US8005619B2 (en) | Method of determining reservoir parameters | |
| WO2017197059A1 (en) | Nuclear magnetic resonance and saturation well logs for determining free water level and reservoir type | |
| CN107524437B (zh) | 确定储层裂缝开度的方法及*** | |
| CN110927794A (zh) | 一种致密储层裂缝识别与孔隙度定量计算的方法 | |
| Xiao et al. | Tight-gas-sand permeability estimation from nuclear-magnetic-resonance (NMR) logs based on the hydraulic-flow-unit (HFU) approach | |
| Kumar et al. | Mapping 3D pore scale fluid distributions: how rock resistivity is influenced by wettability and saturation history | |
| Al-Sulami et al. | The unconventional shale reservoirs of jafurah basin: An integrated petrophysical evaluation using cores and advanced well logs | |
| CN107861162A (zh) | 基于微电极测井资料的天然裂缝识别方法及*** | |
| CN112115592A (zh) | 一种基于岩石物理实验的大斜度/水平井电阻率校正方法 | |
| CN108627878B (zh) | 致密砂岩地层的裂缝识别方法及*** | |
| CN106353813B (zh) | 基于阵列声波测井的流体性质识别方法 | |
| Jin et al. | Capillary pressure prediction from rock models reconstructed using well log data | |
| Kumar et al. | Mapping fluid distributions in 3D at the pore scale: Quantifying the influence of wettability and saturation history on rock resistivity | |
| Al-Bulushi et al. | Capillary Pressure Corrections, Quality Control and Curve Fitting Workflow |
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 |