CN107578343A - 一种碳酸盐岩缝洞型储层含水饱和度计算方法及其装置 - Google Patents
一种碳酸盐岩缝洞型储层含水饱和度计算方法及其装置 Download PDFInfo
- Publication number
- CN107578343A CN107578343A CN201710615928.4A CN201710615928A CN107578343A CN 107578343 A CN107578343 A CN 107578343A CN 201710615928 A CN201710615928 A CN 201710615928A CN 107578343 A CN107578343 A CN 107578343A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- porosity
- mfrac
- water saturation
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明提供了一种碳酸盐岩缝洞型储层含水饱和度计算方法及其装置。该碳酸盐岩缝洞型储层含水饱和度计算方法是在获取碳酸盐岩缝洞型储层相关参数的基础上,基于一种新型胶结指数计算公式获取高精度胶结指数m,然后根据得到的胶结指数m通过阿尔奇公式计算碳酸盐岩的含水饱和度。该新型胶结指数计算公式不仅考虑了孔、洞、缝各部分对储层电性特征的影响,更重要的是可以反映裂缝部分对电性的影响,尤其是裂缝倾角的变化。因此,使得本发明提供的碳酸盐岩储层含水饱和度的计算方法,可以更全面的反映碳酸盐岩储层中复杂的孔隙类型和孔隙结构,从而使得利用该方法计算的含水饱和度更加精确,更加符合实际地质情况。
Description
技术领域
本发明属于油田开发技术领域,具体涉及一种碳酸盐岩缝洞型储层含水饱和度计算方法及其装置。
背景技术
含水饱和度评价是油气储集层定量评价的核心,如何提高储层含水饱和度的计算精度一直以来是储层测井评价的难点,尤其是碳酸盐岩储层。碳酸盐岩由于其孔隙结构复杂、储集空间类型多样且不均匀分布的特点,导致储层具有极强的非均质性,使得传统意义上的测井解释方法存在明显的不适用性(曾文冲,刘学锋,碳酸盐岩非阿尔奇特性的诠释[J],测井技术,2013,37(4):341-351)。储层强烈的非均质性导致储层饱和度解释模型中的胶结指数m和饱和度指数n值的变化范围大,尤其是胶结指数m,而准确的参数值对饱和度的计算至关重要。
阿尔奇公式作为连接电阻率和孔隙度、饱和度之间的桥梁,其是测井计算饱和度的基础,而要想能够准确计算出储层含水饱和度情况,获取准确的胶结指数是关键。
目前研究方法和思路主要分为以下几步:
1)岩电关系分析:首先对取心井的岩心钻孔、获取岩塞样品,然后在实验室中对岩样进行洗油、洗盐、烘干等处理后,测量每块岩心样品孔隙度,最后利用配制的模拟地层水饱和岩心,测量岩心的电阻以及不同含水饱和度状态下的电阻。
2)技术方法:针对不同的层系,分别对同一层系的岩石样品所测出的地层因素、孔隙度和各种含水饱和度状态下的电阻率指数利用数学方法(回归分析)来确定岩石电阻率与孔隙度关系、岩石电阻率与含水饱和度的关系,既而确定所需的饱和度模型参数。
3)技术应用:利用岩电实验得到的饱和度模型参数(a、b、m、n)代入阿尔奇公式中,在获知工区地层水的情况下,可以直接应用到工区其它井的含水饱和度计算中。
目前国内关于碳酸盐岩储层含水饱和度的研究很多,对于碳酸盐岩储层含水饱和度的计算方法主要有以下六大类:基于经典阿尔奇公式的经验性扩展、基于孔隙结构的饱和度模型、基于孔隙尺寸的多重孔隙结构的饱和度模型、基于岩石导电效率理论的饱和度模型、基于网络导电理论的饱和度模型和基于J函数的饱和度模型(田瀚,李昌,贾鹏,碳酸盐岩储层含水饱和度解释模型研究[J].地球物理学进展,2017,32(1):0279-0286.),这些方法要不就是过于简单不能很好的反映出碳酸盐岩储层强烈的非均质性特征或者就是过于复杂无法有效的推广仍处于理论阶段,要不就是只针对某一种类型的碳酸盐岩储层不具有通用性或者就是考虑的因素不够全面。
国外针对碳酸盐岩含水饱和度的研究较多也较早,其中要以Aguilera等提出的三孔隙度模型最为主流(Aguilera R F,A triple porosity model for petrophysicalanalysis of naturally fractured reservoirs[J].Petrophysics,2004,45(2):157-166)。Aguilera根据碳酸盐岩储层的特征,按储集空间类型将储层三分类,即基质孔隙、连通缝洞和非连通缝洞,且认为基质孔隙和连通缝洞为并联导电,后再与非连通缝洞串联导电。在此理论基础之上,经严格公式推导给出了计算胶结指数m的三孔隙度模型,并于2011年再次给出了最新的三孔隙度模型(AI-Ghamdi A,Chen Bo,Behmanesh H,et al.Animproved triple porosity model for evaluation of naturally fracturedreservoirs[J].SPE Reservoir Evaluation&Engineers,2011,14(4):377-384),这种方法在一定程度上改变了利用岩电实验经回归分析得到固定胶结指数的方法,比传统的饱和度解释方法应用方便,而且精度有所提高,但仍待提高。
发明内容
为解决上述问题,本发明的目的是提供一种碳酸盐岩缝洞型储层含水饱和度计算方法。
本发明的另一目的是提供一种碳酸盐岩缝洞型储层含水饱和度计算装置。
为达到上述目的,本发明提供了一种碳酸盐岩缝洞型储层含水饱和度计算方法,该方法是在获取碳酸盐岩缝洞型储层相关参数的基础上,基于以下胶结指数计算公式获取胶结指数m,然后根据得到的胶结指数m通过阿尔奇公式计算碳酸盐岩的含水饱和度;
所述胶结指数计算公式为:
其中,Fθ=0=1/[φ2+(1-φ2)·φb mb];Fθ=90=φ2+(1-φ2)·φb -mb;
上述公式中,m为胶结指数;Φ为总孔隙度,小数;θ为裂缝倾角,°;Φnc为孤立孔洞孔隙度,小数;Fθ=0为水平裂缝情况下混合体的地层因素值;Fθ=90为垂直裂缝情况下混合体的地层因素值;Φ2为裂缝孔隙度,小数;Φb为基质孔隙度,小数;mb为基质部分的胶结指数。
上述胶结指数计算公式不仅考虑了孔、洞、缝各部分对储层电性特征的影响,更重要的是可以反映裂缝部分对电性的影响,尤其是裂缝倾角的变化。可见,该计算公式更加符合实际地质情况,可以获得更可靠的胶结指数m。因此,研究人员可以在此基础上,通过阿尔奇公式计算出更准确的碳酸盐岩缝洞型储层含水饱和度。
在上述碳酸盐岩缝洞型储层含水饱和度计算方法中,优选地,所述胶结指数计算公式的推导过程为:
1.先考虑基质与裂缝混合部分
该部分先假设基质部分的电导率为εe,混合部分中的包含体(裂缝部分)的电导率为εi,整个混合体的电导率为εeff,根据Maxwell Garnett方程可知,沿着任意x方向混合体的电导率为:
其中:Nx为去极化因子;f为包含体(裂缝部分)的孔隙度,小数;εeff,x为沿x方向的混合体电导率。
在此基础上,我们考虑两种极端情况,即在水平裂缝和垂直裂缝的情况下,混合部分的电导率的变化。
当是水平裂缝时,此时去极化因子Nx=0,混合体的电导率处于最大值状态,即:
εeff,max=fεi+(1-f)εe
当是垂直裂缝时,此时去极化因子Nx=1,混合体的电导率处于最大值状态,即:
从而我们就可以得到在水平裂缝和垂直裂缝情况下,混合体的地层因素值,即:
水平裂缝情况下:Fθ=0=Rθ=0/Rw=1/[φ2+(1-φ2)·φb mb]
垂直裂缝情况下:Fθ=90=Rθ=90/Rw=φ2+(1-φ2)·φb -mb
其中:Φ2为裂缝孔隙度,小数;Φb为基质孔隙度,小数;mb为基质部分的胶结指数,可通过岩电分析得到,若无岩电分析资料,可取地区经验值或取2;Fθ=0为水平裂缝情况下混合体的地层因素值;Fθ=90为垂直裂缝情况下混合体的地层因素值。
而对于任意裂缝倾角θ的情况下,混合体的地层因素Ffo可有Fθ=0、Fθ=90和裂缝倾角θ三者来表示,即:
2.再考虑混合部分与孤立孔洞部分
认为基质与裂缝并联导电后再与孤立孔洞部分串联,同理认为混合部分是与孤立孔洞部分是串联导电。假设整体的电阻率为R,基质与裂缝混合部分的电阻率为Rfo,孤立孔洞孔隙度为Φnc,按照串联导电原理,则有:
R=(1-φnc)·Rfo+φnc·Rw
根据阿尔奇公式F=RO/RW,可将上式转化为:
F=(1-φnc)·Ffo+φnc
结合1所推导的混合部分的地层因素表达式,就可以得到整体的地层因素表达式,即:
再根据阿尔奇公式F=RO/RW=Φ-m,就可以得到最终的胶结指数m计算方法,即:
该胶结指数m的计算公式综合考虑了基质孔隙、连通缝洞孔隙、孤立孔洞和裂缝倾角的影响(示意图见图1),相对Aguilera等提出的改进后的三孔隙度模型更加精确,且适用于任何类型的储层。
在上述碳酸盐岩缝洞型储层含水饱和度计算方法中,优选地,该方法具体包括以下步骤:
步骤1,获取胶结指数计算公式中的相关参数;所述相关参数包括:总孔隙度Φ、裂缝倾角θ、孤立孔洞孔隙度Φnc、裂缝孔隙度Φ2、基质孔隙度Φb以及基质部分的胶结指数mb;
步骤2,将获取的相关参数带入胶结指数计算公式中,计算胶结指数m;
步骤3,将获得的胶结指数m代入阿尔奇公式中计算含水饱和度。
在上述碳酸盐岩缝洞型储层含水饱和度计算方法中,优选地,在上述步骤1中,所述总孔隙度Φ是利用变骨架的中子-密度交会获取的,具体计算公式为:
DENma=v1·DENillite+v2·DENquartz+v3·DENcalcite+v4·DENdolomite
CNLma=v1·CNLillite+v2·CNLquartz+v3·CNLcalcite+v4·CNLdolomite
DTma=v1·DTillite+v2·DTquartz+v3·DTcalcite+v4·DTdolomite
上述公式中,DENma、CNLma、DTma分别为变骨架的岩石骨架密度值、补偿中子值和声波时差值,单位分别为g/cm3、V/V、us/ft;DENillite、DENquartz、DENcalcite、DENdolomite分别为泥质、硅质、方解石和白云石的密度值,g/cm3;CNLillite、CNLquartz、CNLcalcite、CNLdolomite分别为泥质、硅质、方解石和白云石的中子值,小数;DTillite、DTquartz、DTcalcite、DTdolomite分别为泥质、硅质、方解石和白云石的声波时差值,us/ft;DEN、CNL分别为测井仪器测量的密度值和补偿中子值,单位分别为g/cm3、小数;DENf、CNLf分别为孔隙流体的密度值和补偿中子值,采用理论值1.0g/cm3和1;ΦD、ΦC分别为利用变骨架密度和变骨架中子计算的密度孔隙度和中子孔隙度,小数;Φ为总孔隙度,小数。
在上述碳酸盐岩缝洞型储层含水饱和度计算方法中,优选地,在上述步骤1中,所述基质孔隙度Φb是利用变骨架声波时差结合Wylie公式获取的,具体计算公式为:
上述公式中,Φb为基质孔隙度,小数;DT为测井仪器测量的声波时差值,us/ft;DTf为孔隙流体的声波时差值,采用理论值189us/ft;DTma为岩石骨架的骨架值,us/ft,通过多矿物模型求得。
在上述碳酸盐岩缝洞型储层含水饱和度计算方法中,优选地,在上述步骤1中,所述孤立孔洞孔隙度Φnc是通过以下公式获取的:
上述公式中,Φ、Φb、Φ2分别为总孔隙度、基质孔隙度和裂缝孔隙度,小数。
在上述碳酸盐岩缝洞型储层含水饱和度计算方法中,优选地,在上述步骤1中,所述基质部分的胶结指数mb是通过以下方式获取的:
通过岩电实验分析得到,即,取无裂缝和溶蚀孔发育的岩心柱塞做岩电实验,从而获取基质部分的胶结指数mb;
当无法进行岩电实验获取数据时,则取研究地区经验值或理论值2。
在上述碳酸盐岩缝洞型储层含水饱和度计算方法中,优选地,在上述步骤1中,获取所述裂缝孔隙度Φ2以及裂缝倾角θ的方式为:
当有成像资料时,通过商业软件对成像资料进行处理和解释,从而获取裂缝孔隙度Φ2和裂缝倾角θ;优选的,所述商业软件为中国石油勘探开发研究院测井与遥感技术研究所研制的Ciflog软件;或者,
当无成像资料时,通过利用深浅侧向电阻率计算裂缝孔隙度Φ2,对于裂缝倾角,通过了解邻井相同层段的情况,取平均值代替。
目前成像测井由于其高分辨率和高覆盖率的特点,是识别裂缝最为有效的手段。而深浅电阻率计算的裂缝孔隙度可能偏大,这时需要利用成像测井解释结果对其进行校正。
在上述碳酸盐岩缝洞型储层含水饱和度计算方法中,优选地,在上述步骤3中,计算地层含水饱和度的阿尔奇公式为:
上述公式中:Sw为含水饱和度,小数;Φe用地层有效孔隙度Φ代替,小数;Rw为地层水电阻率,Ω.m;Rt为地层电阻率,Ω.m;a、b为阿尔奇参数;m为孔隙度胶结指数;n为饱和度指数。
在上述碳酸盐岩缝洞型储层含水饱和度计算方法中,优选地,在上述地层含水饱和度的计算公式中,
地层水电阻率Rw的获取方式为:通过井下地层水取样,分析其矿化度,然后根据以下公式将地层水矿化度转换地层水电阻率;
地层水矿化度转换地层水电阻率的计算公式为:
上述公式中,Cw为地层水矿化度,ppm;T为地层水温度,°F;
阿尔奇参数a、b和饱和度指数n的获取方式为:利用岩-电实验获取;或者,根取采用经验参数:a=b=1、n=2;
地层电阻率Rt的获取方式为:通过CLS-5700测井仪器进行测量获取。
本发明还提供了一种碳酸盐岩缝洞型储层含水饱和度计算装置,该装置包括:
第一单元,所述第一单元用于获取胶结指数计算公式中的相关参数;所述相关参数包括:总孔隙度Φ、裂缝倾角θ、孤立孔洞孔隙度Φnc、裂缝孔隙度Φ2、基质孔隙度Φb以及基质部分的胶结指数mb;
所述胶结指数计算公式为:
其中,Fθ=0=1/[φ2+(1-φ2)·φb mb];Fθ=90=φ2+(1-φ2)·φb -mb;
上述公式中,m为胶结指数;Φ为总孔隙度,小数;θ为裂缝倾角,°;Φnc为孤立孔洞孔隙度,小数;Fθ=0为水平裂缝情况下混合体的地层因素值;Fθ=90为垂直裂缝情况下混合体的地层因素值;Φ2为裂缝孔隙度,小数;Φb为基质孔隙度,小数;mb为基质部分的胶结指数;
第二单元,所述第二单元将获取的相关参数带入胶结指数计算公式中,计算胶结指数m;
第三单元,所述第三单元用于将获得的胶结指数m代入阿尔奇公式中计算含水饱和度。
本次发明提供的方案中,最重要之处是给出了一种胶结指数m的高精度计算公式。在此基础上,使得本发明提供的碳酸盐岩储层含水饱和度的计算方法,可以更全面的反映碳酸盐岩储层中复杂的孔隙类型和孔隙结构,从而使得利用该方法计算的含水饱和度更加精确,更加符合实际地质情况。
附图说明
图1为考虑了基质部分、连通缝洞、孤立孔洞和裂缝倾角影响的等效示意图;
图2为利用ciflog软件处理磨溪204井成像测井计算裂缝孔隙度和裂缝倾角;
图3为利用本发明方法综合解释磨溪204井,计算各种孔隙度、胶结指数m及含水饱和度,并与传统的利用岩电分析得到固定胶结指数m计算的结果进行对比;
图4为利用本发明方法和传统的阿尔奇方法解释的含水饱和度与岩心分析含水饱和度对比图(高石28井)。
具体实施方式
为了对本发明的技术特征、目的和有益效果有更加清楚的理解,现对本发明的技术方案进行以下详细说明,但不能理解为对本发明的可实施范围的限定。
实施例1
本实施例提供了一种碳酸盐岩缝洞型储层含水饱和度计算方法,该计算方法包括以下步骤:
步骤1,获取胶结指数计算公式中的相关参数;所述相关参数包括:总孔隙度Φ、裂缝倾角θ、孤立孔洞孔隙度Φnc、裂缝孔隙度Φ2、基质孔隙度Φb以及基质部分的胶结指数mb;
步骤2,将获取的相关参数带入胶结指数计算公式中,计算胶结指数m;
所述胶结指数计算公式为:
其中,Fθ=0=1/[φ2+(1-φ2)·φb mb];Fθ=90=φ2+(1-φ2)·φb -mb;
上述公式中,m为胶结指数;Φ为总孔隙度,小数;θ为裂缝倾角,°;Φnc为孤立孔洞孔隙度,小数;Fθ=0为水平裂缝情况下混合体的地层因素值;Fθ=90为垂直裂缝情况下混合体的地层因素值;Φ2为裂缝孔隙度,小数;Φb为基质孔隙度,小数;mb为基质部分的胶结指数;
步骤3,将获得的胶结指数m代入阿尔奇公式中计算含水饱和度。
在本实施例中,将上述方法具体应用到了磨溪204井和高石28井的含水饱和度计算中,具体的:
磨溪204井获取相关参数的过程为:
(1)利用多矿物模型对磨溪204井进行处理,其中根据岩屑录井和邻井的岩性情况,在处理过程中,选择的矿物类型有方解石、白云石、石英和伊利石;而所需的地层水、泥浆滤液的电阻率和密度、井温、钻头尺寸等参数可以通过测井参数卡获取,处理后可获得各矿物的含量,具体岩性剖面见图3第三道;
(2)利用多矿物模型得到各矿物含量后,利用发明内容中步骤1的计算方法就可以得到岩石的变骨架密度值、补偿中子值和声波时差值,然后在利用中子-密度交会法计算岩石的总孔隙度和基质孔隙度,在图3第七道;
(3)利用商业软件ciflog可以对磨溪204井成像数据进行处理,处理后可以得到裂缝孔隙度、裂缝倾角等信息,见图2;
(4)在获取各种所需的参数后,将这些参数代入新给出的胶结指数计算公式中,就可以得到随深度变化的变胶结指数m值,具体见图3最后一道;
(5)将得到的变胶结指数m值代入阿尔奇公式中,其中公式中地层水电阻率可通过地层水矿化度获得,最终可以得到含水饱和度,见图3倒数第二道。
图3中显示了本实施例提供的方法和传统的利用岩电分析方法解释磨溪204井的含水饱和度与岩心分析含水饱和度对比结果。通过图3可知,利用变骨架值的中子-密度交会法计算的孔隙度值与岩心分析的孔隙度值吻合很好,且可以发现,随深度变化,胶结指数m存在很明显的变化,这符合碳酸盐岩地层强烈非均质性的特点,而且计算的含水饱和度与岩心分析的很一致,对于下面产水层,利用变胶结指数m计算的含水饱和度在40%-70%,而传统方法计算的含水饱和度在20%-40%,利用变胶结指数m计算的含水饱和度比传统解释模型精度更高,更符合实际试油情况。高石28井获取相关参数的过程为:
高石28井的各参数的获取和胶结指数m的计算方法与上述磨溪204井的相一致,不再重复叙述。
图4中显示了本实施例提供的方法和传统解释方法解释高石28井的含水饱和度与岩心分析含水饱和度对比结果。通过图4可知,本实施提供的方法计算的含水饱和度比传统解释方法计算的含水饱和度精度很高,与岩心分析结果更加吻合。
Claims (10)
1.一种碳酸盐岩缝洞型储层含水饱和度计算方法,其特征在于,该方法是在获取碳酸盐岩缝洞型储层相关参数的基础上,基于以下胶结指数计算公式获取胶结指数m,然后根据得到的胶结指数m通过阿尔奇公式计算碳酸盐岩的含水饱和度;
所述胶结指数计算公式为:
<mrow>
<mi>m</mi>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>n</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&theta;</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>&theta;</mi>
<mo>=</mo>
<mn>90</mn>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&theta;</mi>
<mo>=</mo>
<mn>90</mn>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>+</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>&theta;</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>n</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
<mi>&phi;</mi>
</mrow>
</mfrac>
</mrow>
其中,Fθ=0=1/[φ2+(1-φ2)·φb mb];Fθ=90=φ2+(1-φ2)·φb -mb;
上述公式中,m为胶结指数;Φ为总孔隙度,小数;θ为裂缝倾角,°;Φnc为孤立孔洞孔隙度,小数;Fθ=0为水平裂缝情况下混合体的地层因素值;Fθ=90为垂直裂缝情况下混合体的地层因素值;Φ2为裂缝孔隙度,小数;Φb为基质孔隙度,小数;mb为基质部分的胶结指数。
2.根据权利要求1所述的碳酸盐岩缝洞型储层含水饱和度计算方法,其特征在于,该方法具体包括以下步骤:
步骤1,获取胶结指数计算公式中的相关参数;所述相关参数包括:总孔隙度Φ、裂缝倾角θ、孤立孔洞孔隙度Φnc、裂缝孔隙度Φ2、基质孔隙度Φb以及基质部分的胶结指数mb;
步骤2,将获取的相关参数带入胶结指数计算公式中,计算胶结指数m;
步骤3,将获得的胶结指数m代入阿尔奇公式中计算含水饱和度。
3.根据权利要求2所述的碳酸盐岩缝洞型储层含水饱和度计算方法,其特征在于,在上述步骤1中,所述总孔隙度Φ是利用变骨架值的中子-密度交会获取的,具体计算公式为:
DENma=v1·DENillite+v2·DENquartz+v3·DENcalcite+v4·DENdolomite
CNLma=v1·CNLillite+v2·CNLquartz+v3·CNLcalcite+v4·CNLdolomite
DTma=v1·DTillite+v2·DTquartz+v3·DTcalcite+v4·DTdolomite
<mrow>
<msub>
<mi>&phi;</mi>
<mi>D</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>D</mi>
<mi>E</mi>
<mi>N</mi>
<mo>-</mo>
<msub>
<mi>DEN</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>DEN</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>DEN</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>&phi;</mi>
<mi>C</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>C</mi>
<mi>N</mi>
<mi>L</mi>
<mo>-</mo>
<msub>
<mi>CNL</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>CNL</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>CNL</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
<mrow>
<mi>&phi;</mi>
<mo>=</mo>
<msqrt>
<mfrac>
<mrow>
<msup>
<msub>
<mi>&phi;</mi>
<mi>D</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<msub>
<mi>&phi;</mi>
<mi>C</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
<mn>2</mn>
</mfrac>
</msqrt>
</mrow>
上述公式中,DENma、CNLma、DTma分别为变骨架的岩石骨架密度值、补偿中子值和声波时差值,单位分别为g/cm3、V/V、us/ft;DENillite、DENquartz、DENcalcite、DENdolomite分别为泥质、硅质、方解石和白云石的密度值,g/cm3;CNLillite、CNLquartz、CNLcalcite、CNLdolomite分别为泥质、硅质、方解石和白云石的中子值,小数;DTillite、DTquartz、DTcalcite、DTdolomite分别为泥质、硅质、方解石和白云石的声波时差值,us/ft;DEN、CNL分别为测井仪器测量的密度值和补偿中子值,单位分别为g/cm3、小数;DENf、CNLf分别为孔隙流体的密度值和补偿中子值,采用理论值1.0g/cm3和1;ΦD、ΦC分别为利用变骨架密度和变骨架中子计算的密度孔隙度和中子孔隙度,小数;Φ为总孔隙度,小数。
4.根据权利要求2所述的碳酸盐岩缝洞型储层含水饱和度计算方法,其特征在于,在上述步骤1中,所述基质孔隙度Φb是利用变骨架声波时差结合Wylie公式获取的,具体计算公式为:
<mrow>
<msub>
<mi>&Phi;</mi>
<mi>b</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>D</mi>
<mi>T</mi>
<mo>-</mo>
<msub>
<mi>DT</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>DT</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>DT</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
上述公式中,Φb为基质孔隙度,小数;DT为测井仪器测量的声波时差值,us/ft;DTf为孔隙流体的声波时差值,采用理论值189us/ft;DTma为岩石骨架的骨架值,us/ft,通过多矿物模型求得。
5.根据权利要求2所述的碳酸盐岩缝洞型储层含水饱和度计算方法,其特征在于,在上述步骤1中,所述孤立孔洞孔隙度Φnc是通过以下公式获取的:
<mrow>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>n</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&phi;</mi>
<mo>-</mo>
<msub>
<mi>&phi;</mi>
<mi>b</mi>
</msub>
</mrow>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>&phi;</mi>
<mi>b</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<msub>
<mi>&phi;</mi>
<mn>2</mn>
</msub>
</mrow>
上述公式中,Φ、Φb、Φ2分别为总孔隙度、基质孔隙度和裂缝孔隙度,小数。
6.根据权利要求2所述的碳酸盐岩缝洞型储层含水饱和度计算方法,其特征在于,在上述步骤1中,所述基质部分的胶结指数mb是通过以下方式获取的:
通过岩电实验分析得到,即,取无裂缝和溶蚀孔发育的岩心柱塞做岩电实验,从而获取基质部分的胶结指数mb;
当无法进行岩电实验获取数据时,则取研究地区经验值或理论值2。
7.根据权利要求2所述的碳酸盐岩缝洞型储层含水饱和度计算方法,其特征在于,在上述步骤1中,获取所述裂缝孔隙度Φ2以及裂缝倾角θ的方式为:
当有成像资料时,通过商业软件对成像资料进行处理和解释,从而获取裂缝孔隙度Φ2和裂缝倾角θ;优选的,所述商业软件为中国石油勘探开发研究院测井与遥感技术研究所研制的Ciflog软件;或者,
当无成像资料时,通过利用深浅侧向电阻率计算裂缝孔隙度Φ2,对于裂缝倾角,通过了解邻井相同层段的情况,取平均值代替。
8.根据权利要求2所述的碳酸盐岩缝洞型储层含水饱和度计算方法,其特征在于,在上述步骤3中,计算地层含水饱和度的阿尔奇公式为:
<mrow>
<msub>
<mi>S</mi>
<mi>w</mi>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>a</mi>
<mo>&CenterDot;</mo>
<mi>b</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>R</mi>
<mi>w</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mi>t</mi>
</msub>
<mo>&CenterDot;</mo>
<msup>
<msub>
<mi>&phi;</mi>
<mi>e</mi>
</msub>
<mi>m</mi>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mi>n</mi>
</mfrac>
</msup>
</mrow>
上述公式中:Sw为含水饱和度,小数;Φe用地层有效孔隙度Φ代替,小数;Rw为地层水电阻率,Ω.m;Rt为地层电阻率,Ω.m;a、b为阿尔奇参数;m为孔隙度胶结指数;n为饱和度指数。
9.根据权利要求8所述的碳酸盐岩缝洞型储层含水饱和度计算方法,其特征在于,在上述地层含水饱和度的计算公式中,
地层水电阻率Rw的获取方式为:通过井下地层水取样,分析其矿化度,然后根据以下公式将地层水矿化度转换地层水电阻率;
地层水矿化度转换地层水电阻率的计算公式为:
<mrow>
<msub>
<mi>R</mi>
<mi>w</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>0.0123</mn>
<mo>+</mo>
<mfrac>
<mn>3647.5</mn>
<mrow>
<msup>
<msub>
<mi>C</mi>
<mi>w</mi>
</msub>
<mn>0.955</mn>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>*</mo>
<mfrac>
<mn>45.5</mn>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>+</mo>
<mn>21.5</mn>
<mo>)</mo>
</mrow>
</mfrac>
</mrow>
上述公式中,Cw为地层水矿化度,ppm;T为地层水温度,°F;
阿尔奇参数a、b和饱和度指数n的获取方式为:利用岩-电实验获取,或者,根据地区经验给定;
地层电阻率Rt的获取方式为:通过CLS-5700测井仪器进行测量获取。
10.一种碳酸盐岩缝洞型储层含水饱和度计算装置,其特征在于,该装置包括:
第一单元,所述第一单元用于获取胶结指数计算公式中的相关参数;所述相关参数包括:总孔隙度Φ、裂缝倾角θ、孤立孔洞孔隙度Φnc、裂缝孔隙度Φ2、基质孔隙度Φb以及基质部分的胶结指数mb;
所述胶结指数计算公式为:
<mrow>
<mi>m</mi>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>n</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&theta;</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>&theta;</mi>
<mo>=</mo>
<mn>90</mn>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&theta;</mi>
<mo>=</mo>
<mn>90</mn>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>+</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>&theta;</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>n</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
<mi>&phi;</mi>
</mrow>
</mfrac>
</mrow>
其中,Fθ=0=1/[φ2+(1-φ2)·φb mb];Fθ=90=φ2+(1-φ2)·φb -mb;
上述公式中,m为胶结指数;Φ为总孔隙度,小数;θ为裂缝倾角,°;Φnc为孤立孔洞孔隙度,小数;Fθ=0为水平裂缝情况下混合体的地层因素值;Fθ=90为垂直裂缝情况下混合体的地层因素值;Φ2为裂缝孔隙度,小数;Φb为基质孔隙度,小数;mb为基质部分的胶结指数;
第二单元,所述第二单元将获取的相关参数带入胶结指数计算公式中,计算胶结指数m;
第三单元,所述第三单元用于将获得的胶结指数m代入阿尔奇公式中计算含水饱和度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710615928.4A CN107578343B (zh) | 2017-07-26 | 2017-07-26 | 一种碳酸盐岩缝洞型储层含水饱和度计算方法及其装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710615928.4A CN107578343B (zh) | 2017-07-26 | 2017-07-26 | 一种碳酸盐岩缝洞型储层含水饱和度计算方法及其装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107578343A true CN107578343A (zh) | 2018-01-12 |
CN107578343B CN107578343B (zh) | 2020-07-10 |
Family
ID=61033684
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710615928.4A Active CN107578343B (zh) | 2017-07-26 | 2017-07-26 | 一种碳酸盐岩缝洞型储层含水饱和度计算方法及其装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107578343B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109582922A (zh) * | 2018-11-06 | 2019-04-05 | 广州海洋地质调查局 | 一种水合物饱和度的现场快速判别方法及处理终端 |
CN109632597A (zh) * | 2018-11-06 | 2019-04-16 | 广州海洋地质调查局 | 一种水合物与游离气存在状态的显示方法及处理终端 |
CN110344822A (zh) * | 2018-04-03 | 2019-10-18 | 中国石油天然气股份有限公司 | 含钙质致密储层含水饱和度的确定方法 |
CN110967771A (zh) * | 2019-12-19 | 2020-04-07 | 西南石油大学 | 一种利用常规测井计算裂缝、溶蚀孔洞孔隙度的方法 |
CN110968963A (zh) * | 2019-12-20 | 2020-04-07 | 西南石油大学 | 一种基于电磁理论计算m值的改进三孔隙度模型及方法 |
CN111894568A (zh) * | 2020-08-04 | 2020-11-06 | 中国地质大学(北京) | 裂缝性碳酸盐岩油藏饱和度模型数字岩心分析方法 |
CN111911142A (zh) * | 2020-08-04 | 2020-11-10 | 中国地质大学(北京) | 裂缝性致密砂岩气藏含水饱和度模型数字岩心构建方法 |
CN112627810A (zh) * | 2019-09-24 | 2021-04-09 | 中国石油天然气股份有限公司 | 碳酸盐岩储层中溶洞占比的检测方法 |
CN112696197A (zh) * | 2020-12-30 | 2021-04-23 | 中国石油天然气集团有限公司 | 一种油田储层指数曲线构建方法、系统、设备及存储介质 |
CN113420441A (zh) * | 2021-06-22 | 2021-09-21 | 成都理工大学 | 考虑裂缝倾角的缝洞型储层新三孔隙度模型及构建方法 |
CN113989433A (zh) * | 2021-10-26 | 2022-01-28 | 重庆科技学院 | 一种基于孔隙类型细分的缝洞型储层饱和度模型建立方法 |
CN114429044A (zh) * | 2022-01-21 | 2022-05-03 | 成都理工大学 | 一种建立缝洞型储层电导率模型的方法及应用 |
CN117386359A (zh) * | 2023-12-13 | 2024-01-12 | 中国石油集团川庆钻探工程有限公司 | 一种缝洞型碳酸盐岩储层孔隙度随钻定量评价方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012031264A2 (en) * | 2010-09-02 | 2012-03-08 | Schlumberger Canada Limited | Method to predict dense hydrocarbon saturations for high pressure high temperature |
CN103116192A (zh) * | 2012-07-30 | 2013-05-22 | 中国石油化工股份有限公司 | 一种缝洞型碳酸盐岩油藏储集体建模方法 |
-
2017
- 2017-07-26 CN CN201710615928.4A patent/CN107578343B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012031264A2 (en) * | 2010-09-02 | 2012-03-08 | Schlumberger Canada Limited | Method to predict dense hydrocarbon saturations for high pressure high temperature |
CN103116192A (zh) * | 2012-07-30 | 2013-05-22 | 中国石油化工股份有限公司 | 一种缝洞型碳酸盐岩油藏储集体建模方法 |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110344822A (zh) * | 2018-04-03 | 2019-10-18 | 中国石油天然气股份有限公司 | 含钙质致密储层含水饱和度的确定方法 |
CN109632597A (zh) * | 2018-11-06 | 2019-04-16 | 广州海洋地质调查局 | 一种水合物与游离气存在状态的显示方法及处理终端 |
CN109632597B (zh) * | 2018-11-06 | 2019-08-16 | 广州海洋地质调查局 | 一种水合物与游离气存在状态的显示方法及处理终端 |
CN109582922B (zh) * | 2018-11-06 | 2019-10-18 | 广州海洋地质调查局 | 一种水合物饱和度的现场快速判别方法及处理终端 |
CN109582922A (zh) * | 2018-11-06 | 2019-04-05 | 广州海洋地质调查局 | 一种水合物饱和度的现场快速判别方法及处理终端 |
CN112627810A (zh) * | 2019-09-24 | 2021-04-09 | 中国石油天然气股份有限公司 | 碳酸盐岩储层中溶洞占比的检测方法 |
CN112627810B (zh) * | 2019-09-24 | 2023-08-22 | 中国石油天然气股份有限公司 | 碳酸盐岩储层中溶洞占比的检测方法 |
CN110967771A (zh) * | 2019-12-19 | 2020-04-07 | 西南石油大学 | 一种利用常规测井计算裂缝、溶蚀孔洞孔隙度的方法 |
CN110967771B (zh) * | 2019-12-19 | 2021-04-09 | 西南石油大学 | 一种利用常规测井计算裂缝、溶蚀孔洞孔隙度的方法 |
CN110968963B (zh) * | 2019-12-20 | 2022-05-06 | 西南石油大学 | 一种基于电磁理论计算m值的改进三孔隙度模型及方法 |
CN110968963A (zh) * | 2019-12-20 | 2020-04-07 | 西南石油大学 | 一种基于电磁理论计算m值的改进三孔隙度模型及方法 |
CN111894568A (zh) * | 2020-08-04 | 2020-11-06 | 中国地质大学(北京) | 裂缝性碳酸盐岩油藏饱和度模型数字岩心分析方法 |
CN111911142A (zh) * | 2020-08-04 | 2020-11-10 | 中国地质大学(北京) | 裂缝性致密砂岩气藏含水饱和度模型数字岩心构建方法 |
CN111911142B (zh) * | 2020-08-04 | 2022-05-13 | 中国地质大学(北京) | 裂缝性致密砂岩气藏含水饱和度模型数字岩心构建方法 |
CN112696197A (zh) * | 2020-12-30 | 2021-04-23 | 中国石油天然气集团有限公司 | 一种油田储层指数曲线构建方法、系统、设备及存储介质 |
CN112696197B (zh) * | 2020-12-30 | 2024-05-14 | 中国石油天然气集团有限公司 | 一种油田储层指数曲线构建方法、系统、设备及存储介质 |
CN113420441B (zh) * | 2021-06-22 | 2022-03-11 | 成都理工大学 | 考虑裂缝倾角的缝洞型储层新三孔隙度模型及构建方法 |
CN113420441A (zh) * | 2021-06-22 | 2021-09-21 | 成都理工大学 | 考虑裂缝倾角的缝洞型储层新三孔隙度模型及构建方法 |
CN113989433A (zh) * | 2021-10-26 | 2022-01-28 | 重庆科技学院 | 一种基于孔隙类型细分的缝洞型储层饱和度模型建立方法 |
CN114429044A (zh) * | 2022-01-21 | 2022-05-03 | 成都理工大学 | 一种建立缝洞型储层电导率模型的方法及应用 |
CN117386359A (zh) * | 2023-12-13 | 2024-01-12 | 中国石油集团川庆钻探工程有限公司 | 一种缝洞型碳酸盐岩储层孔隙度随钻定量评价方法 |
CN117386359B (zh) * | 2023-12-13 | 2024-03-08 | 中国石油集团川庆钻探工程有限公司 | 一种缝洞型碳酸盐岩储层孔隙度随钻定量评价方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107578343B (zh) | 2020-07-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107578343A (zh) | 一种碳酸盐岩缝洞型储层含水饱和度计算方法及其装置 | |
CN101487390B (zh) | 一种确定油层原始含油饱和度的阿尔奇模式方法 | |
Wu et al. | Source discrimination of mine water inrush using multiple methods: a case study from the Beiyangzhuang Mine, Northern China | |
CN102175832B (zh) | 一种确定典型储层最佳饱和度计算模型的方法 | |
Ritzi Jr et al. | Relation between anisotropic transmissivity and azimuthal resistivity surveys in shallow, fractured, carbonate flow systems | |
CN107703561A (zh) | 一种页岩气水平井中计算有机碳含量的测井方法 | |
CN108107478B (zh) | 大地电磁同步探测与实时反演方法及系统 | |
CN106875471A (zh) | 煤系含或隔水层三维可视化建模方法 | |
Voeckler et al. | Estimating regional-scale fractured bedrock hydraulic conductivity using discrete fracture network (DFN) modeling | |
CN104863574A (zh) | 一种适用于致密砂岩储层的流体识别方法 | |
CN111058837A (zh) | 一种基于多元逐步回归的页岩油岩性评价方法 | |
CN107515957A (zh) | 泥页岩层序地层划分方法 | |
CN105317429A (zh) | 一种确定页岩地层干酪根含量的方法 | |
CN110501387B (zh) | 一种含天然气水合物沉积物的电阻率测井物理模拟方法 | |
CN105205296A (zh) | 一种求取页岩气储层孔隙度的方法 | |
Kuchuk et al. | Determination of in situ two-phase flow properties through downhole fluid movement monitoring | |
CN103604836A (zh) | 一种测定天然气水合物储层饱和度的方法及设备 | |
Feng et al. | Accurate determination of water saturation in tight sandstone gas reservoirs based on optimized Gaussian process regression | |
Yuan et al. | Fracture characterization in oil-based mud boreholes using image logs: example form tight sandstones of Lower Cretaceous Bashijiqike Formation of KS5 well area, Kuqa Depression, Tarim Basin, China | |
CN105275459A (zh) | 一种确定页岩地层可动水体积含量的方法 | |
Aghli et al. | Structural and fracture analysis using EMI and FMI image Log in the carbonate Asmari reservoir (Oligo-Miocene), SW Iran | |
Holden et al. | Integration of production logs helps to understand heterogeneity of Mishrif reservoir in Rumaila | |
Teikeu et al. | Hydrogeophysical parameters estimation for aquifer characterisation in hard rock environment: a case study from Yaounde, Cameroon | |
Luo et al. | Three‐Dimensional Steady‐State Hydraulic Tomography Analysis With Integration of Cross‐Hole Flowmeter Data at a Highly Heterogeneous Site | |
Li et al. | In situ estimation of relative permeability from resistivity measurements |
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 |