CN113989433A - 一种基于孔隙类型细分的缝洞型储层饱和度模型建立方法 - Google Patents
一种基于孔隙类型细分的缝洞型储层饱和度模型建立方法 Download PDFInfo
- Publication number
- CN113989433A CN113989433A CN202111249858.8A CN202111249858A CN113989433A CN 113989433 A CN113989433 A CN 113989433A CN 202111249858 A CN202111249858 A CN 202111249858A CN 113989433 A CN113989433 A CN 113989433A
- Authority
- CN
- China
- Prior art keywords
- pore
- fracture
- porosity
- model
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种基于孔隙类型细分的缝洞型储层饱和度模型建立方法,包括以下步骤,步骤1:确定缝洞型储层孔隙类型;步骤2:基于缝洞型储层中的孔隙类型,进行饱和度模型分析;步骤3:建立缝洞型储层饱和度模型;步骤4:确定饱和度模型参数;步骤5:依据所述饱和度模型和饱和度参数进行饱和度计算。本发明解决了现有技术在各向异性储层中的饱和度模型中孔隙指数m和饱和度指数n变化范围大(不适用性)的问题。
Description
技术领域
本发明属于地质勘探领域,具体涉及一种基于孔隙类型细分的缝洞型储层饱和度模型建立方法。
背景技术
饱和度评价是油气储集层定量评价的核心,如何提高储层含水饱和度的计算精度一直以来是储层测井评价的难点。缝洞型储层由于其孔隙结构的复杂性和储集空间的多元性,以及极不均匀的随机分布,导致储层强烈的非均质性,使得建立在均质、各向同性地层基础上的传统测井解释方法出现了明显的不适用性,其原因就在于强烈的非均质性导致储层饱和度模型中孔隙度指数m_ex和饱和度指数n_ex值的变化范围加大。
缝洞型储层的孔隙度指数、饱和度指数以及基质、裂缝、溶蚀孔洞的孔隙度参数是储层测井评价中最重要的参数,也是测井界一直努力攻关的难题,前人一直在这方便展开大量的研究工作,取得了一些重要的理论成果并且发现相应的问题。
基于孔隙类型的多重孔隙结构的饱和度模型中,由于孔隙类型多样,不同孔隙类型,其导电机理存在较大差异。一些现有技术重点考虑碳酸盐岩中基质孔隙、裂缝和溶蚀孔隙之间的串、并联方式,从而建立等效岩石导电网络。但是采用串、并联方式来等效地层真实导电网络是缺乏真实导电机理的研究与证实。
发明内容
本发明要解决的技术问题是:
(1)现有技术缝洞型储层不同孔隙类型导电机理不明确;
(2)现有技术的经验关系中饱和度指数与孔隙度指数缺乏理论依据和难以确定合适参数值的问题;
(3)现有技术缝洞型储层饱和度评价的困难。
本发明提供一种基于孔隙类型细分的缝洞型储层饱和度模型建立方法,包括以下步骤,
步骤1:确定缝洞型储层孔隙类型;
步骤2:基于缝洞型储层中的孔隙类型,进行饱和度模型分析;
步骤3:建立缝洞型储层饱和度模型;
步骤4:确定饱和度模型参数;
步骤5:依据所述饱和度模型和饱和度参数进行饱和度计算。
进一步的,所述步骤1包括
步骤1.1、储层柱塞岩心NMR实验获取T2谱转换为伪毛管压力曲线或压汞实验获取毛管压力曲线;
步骤1.2、依据毛管压力曲线、岩心二维图像与岩心层段电成像测井图像特征确定缝洞型储层中发育的孔隙类型。
进一步的,所述步骤2包括,
步骤2.1基于缝洞型储层中发育的孔隙类型,采用基质岩心嵌套平板状裂缝模型构建基质孔隙与裂缝发育的复合孔隙类型的三维数字岩心模型;
步骤2.2基于构建的三维数字岩心模型,采用有限元法模拟岩心电阻率,模拟结果同基于基质孔隙与裂缝并联导电的理论模型进行对比,验证基质孔隙与裂缝孔隙的串并联耦合关系。
进一步的,所述步骤3包括,
步骤3.1.基于孔隙类型划分和不同类型孔隙串并联导电耦合关系,建立岩心完全含水条件下的孔隙度模型;
采用以下公式将基质孔隙与裂缝构成复合系统电阻率Rfo表示为基质孔隙电阻率R0与裂缝孔隙电阻率并联结果
式中,γ2=φ2/φ,表明裂缝孔隙度φ2占总孔隙度φ的比值,Rw表示地层水电阻率,R0表示100%饱含地层水岩石电阻率,mb表示基质孔隙部分对应孔隙度指数,φfracture为裂缝孔隙度,φvug为溶蚀孔洞孔隙度,φ为总孔隙度
采用以下公式建立裂缝与基质孔隙发育岩心的孔隙度模型
Φm表示基质孔隙的孔隙度,FF表示地层因素,m_ex为孔隙度指数。
步骤3.2.基于孔隙类型划分和不同类型孔隙串并联导电耦合关系,采用以下公式建立基质孔隙与裂缝发育的部分含水岩心的饱和度导电模型:
式中,Rt表示储层电阻率,Sw为含水饱和度,λp为基质孔隙的孔隙分布参数;λf为裂缝孔隙的孔隙分布参数;
步骤3.3.基于基质孔隙与裂缝发育的孔隙度模型与饱和度导电模型,采用以下公式建立基质孔隙与裂缝发育的饱和度模型;
该模型的参数,包括基质孔隙度φm,裂缝孔隙度φ2,总孔隙度φ,以及基质孔隙和裂缝孔隙的孔隙分布参数λP和λf。
进一步的,所述步骤4包括
步骤4.1依据三孔隙、侧向测井曲线计算基质孔隙度φm,裂缝孔隙度φ2,总孔隙度φ;
步骤4.2依据研究层段岩心压汞曲线在划分孔隙类型后,分段计算基质孔隙与裂缝孔隙的孔隙分布参数λP和λf。
本发明的有益效果是,实现了缝洞型储层中不同类型孔隙导电机理的研究,同时实现了孔隙度与饱和度模型的理论推导,并最终提出了基于孔隙类型细分的测井饱和度模型,从而解决已有饱和度模型缺乏导电机理的研究,解决现有技术在各向异性储层中的饱和度模型中孔隙指数m_ex和饱和度指数n_ex变化范围大(不适用性)的问题。本发明饱和度模型的提出有助于提高缝洞型储层的测井饱和度评价,从而提高流体识别能力和测井解释准确度。同时该技术还可以推广到孔隙类型发育的致密砂岩储层中,从而进一步提高缝洞孔隙类型发育的碳酸盐岩、致密砂岩等非常规性储层的测井饱和度评价。
附图说明
图1.基于孔隙类型细分的缝洞型储层饱和度模型建立流程图。
图2.依据(a)电成像资料,(b)岩心压汞数据,(c)岩心二维图片识别研究层段主要发育的孔隙类型。
图3.依据岩心嵌套技术将裂缝(a)嵌入基质孔隙岩心模型(b)中建立基质孔隙-裂缝发育的双孔隙岩心模型(c)。
图4.依据基质孔隙与裂缝孔隙并联导电(a)计算的理论数值与基于三维数字岩心电阻率模拟结果进行对比(b)。
图5.(a)依据毛管压力曲线划分孔隙类型后获得不同孔隙类型的分布参数,(b)依据该方法获取不同深度岩心的孔隙分布参数和计算的饱和度指数。
图6.针对研究层段依据三孔隙度曲线(第2道)计算饱和度模型参数(孔隙度,第6道),而后采用基于孔隙度类型细分的饱和度模型计算含水饱和度SWRT与现有技术饱和度SWAR对比(第七道),发现新饱和度模型计算饱和度与试油试采结果相吻合。
具体实施方式
下面对本发明推导过程进行说明。
本发明包括以下步骤
第一.储层孔隙类型细分;缝洞型储层中孔隙类型为基质孔隙和次生孔隙(如裂缝、溶蚀孔洞等)。
1、储层柱塞岩心NMR实验获取T2谱转换为伪毛管压力曲线或压汞实验获取毛管压力曲线;
2、依据毛管压力曲线、岩心二维图像与岩心层段电成像测井图像特征确定缝洞型储层中发育的孔隙类型;
(1)毛管压力识别孔隙类型依据:由于一般情况下,从裂缝,连通孔洞到连通孔隙的大小依次降低。压汞实验中,随着压力增大,汞依次进入裂缝,孔洞和基质孔隙,压汞曲线上表现出分段特点:裂缝段-毛管阻力小,毛管压力曲线几乎贴近横坐标分布;孔洞段与裂缝段-存在明显台阶,表明孔洞与裂缝尺寸的突变,而非渐变;孔隙段与孔洞段-存在明显的台阶,表明二者孔喉直径的突变。
所以,可依据毛管压力曲线的特点确定缝洞型储层中发育的孔隙类型:裂缝型-毛管压力曲线紧贴横坐标,汞充满裂缝后,压力增大不再进入,出现垂直段;孔(洞)缝型-裂缝段为较长垂直段,而后汞进入孔隙(孔洞);孔洞缝型-裂缝段,孔洞段,孔隙段三者发育,呈阶梯演化。
(2)电成像资料识别孔隙类型依据:
一般情况下,有溶孔溶洞发育的地层,由于井壁的孔洞里充满了低电阻的泥浆或钻井液,在电成像测井图像上,孔洞常常为深色或黑色特征;在形态特征方面,溶孔在图像上一般表现为类圆状、斑点状,溶洞表现为不规则的块状,延伸比较短的片状和条带分布。
裂缝在电成像测井图像上呈现正弦曲线的形态:高导缝在图像上表现为连续的黑色正弦曲线;高阻缝在图像上往往表现为晕圈状的亮黄色-白色正弦曲线,反映裂缝被方解石填充,属于无效缝;诱导缝是钻井过程中产生的,对储层原始孔渗空间没有贡献,在图像上往往分布在相距180度图像的两侧。
第二.缝洞型储层饱和度导电模型研究;
1、岩心二维图像确定储层孔隙类型以及基质孔隙、裂缝与溶蚀孔洞的接触关系;
2、岩心三维重建技术构建孔缝、孔洞或孔洞缝发育的三维数字岩心;
3、有限元计算储层三维数字岩心电阻率;
4、不同孔隙类型串并联耦合关系的研究与验证;依据有限元计算的孔缝、孔洞或孔洞缝等复合孔隙类型的三维数字岩心电阻率,与单一孔隙类型岩心电阻率依据串并联公式推导出的复合孔隙类型三维数字岩心电阻率对比,研究孔缝、孔洞或孔洞缝等复合孔隙类型的岩心中各个孔隙类型的串并联耦合关系。
第三.基于不同孔隙串并联耦合关系,理论推导地层水电阻率、完全含水岩石电阻率与孔隙度、孔隙度指数m_ex关系的孔隙度模型,以及部分含水岩石电阻率、完全含水岩石电阻率与饱和度、饱和度指数n_ex关系的饱和度模型的建立。
1、推导缝洞型孔隙度、孔隙度指数、完全含水岩心电阻率、地层水电阻率关系的孔隙度模型;
首先:依据孔隙类型以及孔隙类型串并联耦合关系,推导出复合岩心完全含水岩石电阻率;
(1)当存在两种孔隙类型串联导电时,岩心完全饱和地层水电阻率Rfnco可表示为孔隙类型I与孔隙类型II电阻率的串联结果(公式1):
式中,φI和φII分别表示两种孔隙类型的孔隙度;RI0和RII0分别表示完全含水岩心中孔隙类型I和孔隙类型II电阻率;
(2)当存在两种孔隙类型并联导电时,岩心完全饱和地层水电阻率Rfnco可表示为孔隙类型I与孔隙类型II电阻率的并联结果(公式2):
式中,φI和φII分别表示两种孔隙类型的孔隙度;RI0和RII0分别表示完全含水岩心中孔隙类型I和孔隙类型II电阻率;
(3)当存在三种孔隙类型串联导电时,岩心完全饱和地层水电阻率Rfnco可表示为孔隙类型I、孔隙类型II与孔隙类型III电阻率的串联结果(公式3):
式中,φI、φII和φIII分别表示三种孔隙类型的孔隙度;RI0、RII0和RIII0分别表示完全含水岩心中孔隙类型I、孔隙类型II和孔隙类型III的电阻率;
(4)当存在三种孔隙类型并联导电时,岩心完全饱和地层水电阻率Rfnco可表示为孔隙类型I、孔隙类型II与孔隙类型III电阻率的并联结果(公式4):
式中,φI、φII和φIII分别表示三种孔隙类型的孔隙度;RI0、RII0和RIII0分别表示完全含水岩心中孔隙类型I、孔隙类型II和孔隙类型III的电阻率;
(5)当存在三种孔隙类型中两种首先并联而后与第三种串联导电时,岩心完全饱和地层水电阻率Rfnco可表示为孔隙类型I、孔隙类型II首先并联而后与孔隙类型III串联的结果(公式5):
式中,φI、φII和φIII分别表示三种孔隙类型的孔隙度;RI0、RII0和RIII0分别表示完全含水岩心中孔隙类型I、孔隙类型II和孔隙类型III的电阻率;
(6)当存在三种孔隙类型中两种串联而后与第三种并联导电时,岩心完全饱和地层水电阻率Rfnco可表示为孔隙类型I、孔隙类型II首先串联而后与孔隙类型III并联的结果(公式6):
式中,φI、φII和φIII分别表示三种孔隙类型的孔隙度;RI0、RII0和RIII0分别表示完全含水岩心中孔隙类型I、孔隙类型II和孔隙类型III的电阻率;
其次,依据存在的孔隙类型,确定完全含水岩心时,各个孔隙类型的电阻率;
(1)岩心完全含水时,基质孔隙电阻率(公式7)
式中,mb,基质孔隙部分对应孔隙度指数,一般取2.0;φ总孔隙度,(小数);φfracture,为裂缝孔隙度,小数;φvug,为溶蚀孔洞孔隙度,(小数);
(2)岩心完全含水时,裂缝孔隙电阻率(公式8)
R0=φfracture, (8)
式中,φfracture,为裂缝孔隙度,小数;
(3)岩心完全含水时,岩心溶蚀孔隙电阻率(公式9)
R0=φvug, (9)
式中,φvug,为溶蚀孔隙的孔隙度,小数;
最后,定义地层因素FF为Rw(地层水电阻率)与R0(100%饱含地层水岩石电阻率)的比值,建立出地层孔隙度模型:即孔隙度φ,孔隙度指数m_ex与复合孔隙类型岩心电阻率Rfnco之间的定量关系。
2、缝洞型储层中部分含水岩石电阻率、完全含水岩石电阻率与饱和度的关系即饱和度模型的推导。
首先:依据电流与渗流相似性原理结合渗流方程,推导出电阻率增大指数(部分含水岩石电阻率与完全含水岩石电阻率的比值)与含水饱和度、孔隙结构(孔隙分布参数)的关系;
(1)对于纯净的均匀介质,不考虑粘土导电的情况,获取电阻率增大指数与水相相对渗透率的关系如下(公式11)
式中,RI电阻率增大指数,Rt储层电阻率;R0完全饱含地层水储层电阻率;krw,地层水相对渗透率。
(2)依据渗流模型方程(公式12),
式中λ取决于孔隙大小与孔隙分布,称为孔隙分布参数
和毛管压力曲线确定λ的方法(公式13),
Pc=Pe(Sw)-1/λ, (13)
式中,Pc毛管压力,Pe毛管入口压力;
(3)给出电阻率增大指数的表达式(14)。
其次,依据验证后的各孔隙类型串并联关系,推导出孔洞、孔缝双孔隙类型或孔洞缝三孔隙类型发育的复合岩心部分含水岩心电阻率;
(1)当存在两种孔隙类型串联导电时,岩心部分饱和地层水电阻率Rfnct可表示为孔隙类型I与孔隙类型II电阻率的串联结果(公式15):
式中,φI和φII分别表示两种孔隙类型的孔隙度;RIt和RIIt分别表示部分含水岩心中孔隙类型I和孔隙类型II电阻率;
(2)当存在两种孔隙类型并联导电时,岩心部分饱和地层水电阻率Rfnct可表示为孔隙类型I与孔隙类型II电阻率的并联结果(公式16):
式中,φI和φII分别表示两种孔隙类型的孔隙度;RIt和RIIt分别表示部分含水岩心中孔隙类型I和孔隙类型II电阻率;
(3)当存在三种孔隙类型串联导电时,岩心部分饱和地层水电阻率Rfnct可表示为孔隙类型I、孔隙类型II与孔隙类型III电阻率的串联结果(公式17):
式中,φI、φII和φIII分别表示三种孔隙类型的孔隙度;RIt、RIIt和RIIIt分别表示部分含水岩心中孔隙类型I、孔隙类型II和孔隙类型III的电阻率;
(4)当存在三种孔隙类型并联导电时,岩心部分饱和地层水电阻率Rfnct可表示为孔隙类型I、孔隙类型II与孔隙类型III电阻率的并联结果(公式18):
式中,φI、φII和φIII分别表示三种孔隙类型的孔隙度;RIt、RIIt和RIIIt分别表示部分含水岩心中孔隙类型I、孔隙类型II和孔隙类型III的电阻率;
(5)当存在三种孔隙类型中两种首先并联而后与第三种串联导电时,岩心部分饱和地层水电阻率Rfnct可表示为孔隙类型I、孔隙类型II首先并联而后与孔隙类型III串联的结果(公式19):
式中,φI、φII和φIII分别表示三种孔隙类型的孔隙度;RIt、RIIt和RIIIt分别表示部分含水岩心中孔隙类型I、孔隙类型II和孔隙类型III的电阻率;
(6)当存在三种孔隙类型中两种首先串联而后与第三种并联时,岩心部分饱和地层水电阻率Rfnct可表示为孔隙类型I、孔隙类型II首先串联而后与孔隙类型III并联的结果(公式20):
式中,φI、φII和φIII分别表示三种孔隙类型的孔隙度;RIt、RIIt和RIIIt分别表示部分含水岩心中孔隙类型I、孔隙类型II和孔隙类型III的电阻率;
再次,依据存在的孔隙类型,确定岩石部分含水时孔隙类型电阻率;
(1)岩心部分含水时,基质孔隙电阻率(公式21)
式中,φm,基质孔隙的孔隙度(小数);φ总孔隙度,(小数);φfracture,为裂缝孔隙度,小数;φvug,为溶蚀孔洞孔隙度,(小数);λp为基质孔隙的孔隙分布参数;Sw为含水饱和度(小数).
(2)岩心部分含水时,裂缝孔隙电阻率(公式22)
式中,φfracture,为裂缝孔隙度,小数;λf为裂缝孔隙的孔隙分布参数;Sw为含水饱和度(小数).
(3)岩心部分含水时,溶蚀孔隙电阻率(公式23)
式中,φvug,为溶蚀孔隙的孔隙度,小数;λv为裂缝孔隙的孔隙分布参数;Sw为含水饱和度(小数).
最后,定义电阻率增大指数RI为Rt(部分含水岩石电阻率)与R0(100%饱含地层水岩石电阻率)的比值,建立出地层饱和度模型:即饱和度Sw,饱和度指数n_ex与复合孔隙类型岩心部分含水岩石电阻率Rfnct之间的定量关系。
3、基于孔隙度模型和饱和度模型,建立新的饱和度模型。
第四.缝洞型储层饱和度模型参数确定;新建立的饱和度模型中,主要的模型参数为:各个孔隙类型的孔隙度和孔隙分布参数。
1、在缝洞型储层孔隙类型划分后的基础上,依据三孔隙度、侧向测井和电成像测井资料计算出储层基质孔隙度φm,裂缝孔隙度φfracture,溶蚀孔洞孔隙度φvug和总孔隙度φ。
2、缝洞型储层中孔隙分布参数λ的确定;
基于毛管压力曲线,首先对毛管压力曲线进行分段识别孔隙类型;
其次,在每个孔隙类型对应压力曲线段内,依据毛管入口压力Pe与幂函数拟合(公式25);
Pc=Pe(Sw)-1/λ, (25)
式中,Pc毛管压力,Pe毛管入口压力。
最后,确定各个孔隙类型的孔隙分布参数;
第五.新的饱和度模型应用;
针对研究层段建立的饱和度模型和模型参数,进行饱和度模型解释工作,验证模型准确性。
下面结合附图和具体实施例对本发明的具体实施方式做进一步说明:
实施例1
以中国某油田例井Z碳酸盐岩缝洞型储层3260m-3300m基于孔隙类型细分的饱和度模型提出为例,图1为基于孔隙类型细分的饱和度模型提出流程图,一种基于孔隙类型细分的缝洞型储层饱和度模型建立方法,具体包括以下步骤:
步骤1:缝洞型储层孔隙类型细分,具体包括以下步骤:
依据研究层段的电成像资料,岩心压汞数据和岩心二维图片识别出研究层段主要发育的孔隙类型(如图2)。
从图2(a)中可以见黑色正弦曲线,表明研究层段裂缝发育;
图2(b)的压汞曲线明显存在一个台阶,首先在较低压力下,汞迅速进行缝洞型储层的孔隙,这是压汞曲线基本与横坐标重合,表明孔隙类型为裂缝;当汞充满裂缝后,压力增大,汞饱和度不变,直到压力增大到40Mpa,汞再次进入缝洞型储层的孔隙,从压汞曲线可以看出该层段存在两种孔隙类型;
图2(c)可以看出,岩心中发育(黑色)细长延伸的裂缝和(黑色)基质孔隙。
通过上述数据,可以确定研究层段缝洞型储层中的孔隙类型为基质孔隙-裂缝型。
步骤2:基于缝洞型储层中的孔隙类型,进行饱和度模型分析;
步骤2.1.基于缝洞型储层中发育的孔隙类型,采用基质岩心(b)嵌套平板状裂缝模型构(a)建基质孔隙与裂缝发育的复合孔隙类型的三维数字岩心模型(c),如图3所示。
步骤2.2.基于构建的三维数字岩心模型,采用有限元法模拟岩心电阻率,模拟结果同基于基质孔隙与裂缝并联导电的理论模型进行对比,验证基质孔隙与裂缝孔隙的串并联耦合关系,如图4。
图4(a)给出当基质孔隙与裂缝并联导电时,岩心总的电阻率Rfo倒数等于裂缝电阻率Rw倒数与裂缝孔隙度在总孔隙度中比例(γ2=φ2/φ)的乘积与基质孔隙电阻率R0的倒数乘以基质孔隙度在总孔隙度中的比例(1-γ2)之和,公式(26)
图4(b)对比基于三维数字岩心采用有限元模拟计算的岩心电阻率,与基质孔隙与裂缝并联导电理论结果,理论结果与模拟结果基本吻合,表明在基质孔隙与裂缝孔隙发育的缝洞型储层中,两种孔隙类型并联导电。
步骤3:缝洞型储层饱和度模型建立;
步骤3.1.基于孔隙类型划分和不同类型孔隙串并联导电耦合关系的研究,采用基质孔隙与裂缝并联导电的思路建立岩心完全含水条件下的孔隙度模型。
基质孔隙与裂缝构成复合系统电阻率Rfo可以表示为基质孔隙电阻率R0与裂缝孔隙电阻率并联结果(公式27)。
式中,γ2=φ2/φ,表明裂缝孔隙度,φ2占总孔隙度φ的比值,小数。
建立裂缝与基质孔隙发育岩心的孔隙度模型(公式28)。
步骤3.2.基于孔隙类型划分和不同类型孔隙串并联导电耦合关系,建立基质孔隙与裂缝发育的部分含水岩心的饱和度导电模型(公式29)。
步骤3.3.基于基质孔隙与裂缝发育的孔隙度模型与饱和度导电模型,建立基质孔隙与裂缝发育的饱和度模型(公式30)。
该模型的关键参数是:基质孔隙度φm,裂缝孔隙度φ2,总孔隙度φ,以及基质孔隙和裂缝孔隙的孔隙分布参数λP和λf。
步骤4:确定饱和度模型参数;
步骤4.1.依据三孔隙、侧向测井曲线计算基质孔隙度φm,裂缝孔隙度φ2,总孔隙度φ。
步骤4.2.依据研究层段岩心压汞曲线在划分孔隙类型后,分段计算基质孔隙与裂缝孔隙的孔隙分布参数λP和λf。如图5所示。
步骤5:依据上述饱和度模型和饱和度参数对研究层段进行饱和度计算,如图6所示。基于三孔隙度(第2道)计算基质与裂缝孔隙度(第6道)的饱和度模型参数,计算含水饱和度(第7道),新模型的含水饱和度SWRT比现有技术饱和度SWAR更符合测井解释结论和实际试油试采结果。
Claims (5)
1.一种基于孔隙类型细分的缝洞型储层饱和度模型建立方法,其特征在于,包括以下步骤,
步骤1:确定缝洞型储层孔隙类型;
步骤2:基于缝洞型储层中的孔隙类型,进行饱和度模型分析;
步骤3:建立缝洞型储层饱和度模型;
步骤4:确定饱和度模型参数;
步骤5:依据所述饱和度模型和饱和度参数进行饱和度计算。
2.如权利要求1所述的一种基于孔隙类型细分的缝洞型储层饱和度模型建立方法,其特征在于,所述步骤1包括
步骤1.1、储层柱塞岩心NMR实验获取T2谱转换为伪毛管压力曲线或压汞实验获取毛管压力曲线;
步骤1.2、依据毛管压力曲线、岩心二维图像与岩心层段电成像测井图像特征确定缝洞型储层中发育的孔隙类型。
3.如权利要求1所述的一种基于孔隙类型细分的缝洞型储层饱和度模型建立方法,其特征在于,所述步骤2包括,
步骤2.1基于缝洞型储层中发育的孔隙类型,采用基质岩心嵌套平板状裂缝模型构建基质孔隙与裂缝发育的复合孔隙类型的三维数字岩心模型;
步骤2.2基于构建的三维数字岩心模型,采用有限元法模拟岩心电阻率,模拟结果同基于基质孔隙与裂缝并联导电的理论模型进行对比,验证基质孔隙与裂缝孔隙的串并联耦合关系。
4.如权利要求1所述的一种基于孔隙类型细分的缝洞型储层饱和度模型建立方法,其特征在于,所述步骤3包括,
步骤3.1.基于孔隙类型划分和不同类型孔隙串并联导电耦合关系,建立岩心完全含水条件下的孔隙度模型;
采用以下公式将基质孔隙与裂缝构成复合系统电阻率Rfo表示为基质孔隙电阻率R0与裂缝孔隙电阻率并联结果
式中,γ2=φ2/φ,表明裂缝孔隙度φ2占总孔隙度φ的比值,Rw表示地层水电阻率,R0表示100%饱含地层水岩石电阻率,mb表示基质孔隙部分对应孔隙度指数,φfracture为裂缝孔隙度,φvug为溶蚀孔洞孔隙度,φ为总孔隙度;
采用以下公式建立裂缝与基质孔隙发育岩心的孔隙度模型
Φm表示基质孔隙的孔隙度,FF表示地层因素,m_ex为孔隙度指数。
步骤3.2.基于孔隙类型划分和不同类型孔隙串并联导电耦合关系,采用以下公式建立基质孔隙与裂缝发育的部分含水岩心的饱和度导电模型:
式中,Rt表示部分含水岩石电阻率,Sw为含水饱和度,λp为基质孔隙的孔隙分布参数;λf为裂缝孔隙的孔隙分布参数;
步骤3.3.基于基质孔隙与裂缝发育的孔隙度模型与饱和度导电模型,采用以下公式建立基质孔隙与裂缝发育的饱和度模型;
该模型的参数,包括基质孔隙度φm,裂缝孔隙度φ2,总孔隙度φ,以及基质孔隙和裂缝孔隙的孔隙分布参数λP和λf。
5.如权利要求1所述的一种基于孔隙类型细分的缝洞型储层饱和度模型建立方法,其特征在于,所述步骤4包括
步骤4.1依据三孔隙、侧向测井曲线计算基质孔隙度φm,裂缝孔隙度φ2,总孔隙度φ;
步骤4.2依据研究层段岩心压汞曲线在划分孔隙类型后,分段计算基质孔隙与裂缝孔隙的孔隙分布参数λP和λf。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111249858.8A CN113989433B (zh) | 2021-10-26 | 2021-10-26 | 一种基于孔隙类型细分的缝洞型储层饱和度模型建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111249858.8A CN113989433B (zh) | 2021-10-26 | 2021-10-26 | 一种基于孔隙类型细分的缝洞型储层饱和度模型建立方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113989433A true CN113989433A (zh) | 2022-01-28 |
CN113989433B CN113989433B (zh) | 2022-06-07 |
Family
ID=79741785
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111249858.8A Active CN113989433B (zh) | 2021-10-26 | 2021-10-26 | 一种基于孔隙类型细分的缝洞型储层饱和度模型建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113989433B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117421890A (zh) * | 2023-10-19 | 2024-01-19 | 重庆科技学院 | 一种基于矿物填充的缝洞型储层导电模型构建方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103543474A (zh) * | 2012-07-16 | 2014-01-29 | 中国石油化工股份有限公司 | 一种基于非导电孔隙模型的含油气饱和度评价方法 |
CN107578343A (zh) * | 2017-07-26 | 2018-01-12 | 中国石油天然气股份有限公司 | 一种碳酸盐岩缝洞型储层含水饱和度计算方法及其装置 |
CN110320139A (zh) * | 2018-03-30 | 2019-10-11 | 中国石油化工股份有限公司 | 缝洞型储层裂缝孔隙度定量评价方法及系统 |
WO2020224539A1 (zh) * | 2019-05-06 | 2020-11-12 | 西安华线石油科技有限公司 | 一种基于广义管流渗流耦合的流动模拟及瞬变井分析方法 |
US20210255349A1 (en) * | 2020-02-17 | 2021-08-19 | China Petroleum & Chemical Corporation | Computer-implemented method and system for small cave recognition using seismic reflection data |
CN113420441A (zh) * | 2021-06-22 | 2021-09-21 | 成都理工大学 | 考虑裂缝倾角的缝洞型储层新三孔隙度模型及构建方法 |
-
2021
- 2021-10-26 CN CN202111249858.8A patent/CN113989433B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103543474A (zh) * | 2012-07-16 | 2014-01-29 | 中国石油化工股份有限公司 | 一种基于非导电孔隙模型的含油气饱和度评价方法 |
CN107578343A (zh) * | 2017-07-26 | 2018-01-12 | 中国石油天然气股份有限公司 | 一种碳酸盐岩缝洞型储层含水饱和度计算方法及其装置 |
CN110320139A (zh) * | 2018-03-30 | 2019-10-11 | 中国石油化工股份有限公司 | 缝洞型储层裂缝孔隙度定量评价方法及系统 |
WO2020224539A1 (zh) * | 2019-05-06 | 2020-11-12 | 西安华线石油科技有限公司 | 一种基于广义管流渗流耦合的流动模拟及瞬变井分析方法 |
US20210255349A1 (en) * | 2020-02-17 | 2021-08-19 | China Petroleum & Chemical Corporation | Computer-implemented method and system for small cave recognition using seismic reflection data |
CN113420441A (zh) * | 2021-06-22 | 2021-09-21 | 成都理工大学 | 考虑裂缝倾角的缝洞型储层新三孔隙度模型及构建方法 |
Non-Patent Citations (2)
Title |
---|
沈金松等: "缝洞型储层中次生孔隙类型对胶结指数和饱和度指数影响的机理分析", 《测井技术》 * |
赵辉等: "裂缝性储层孔隙指数、饱和度及裂缝孔隙度计算研究", 《地球物理学进展》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117421890A (zh) * | 2023-10-19 | 2024-01-19 | 重庆科技学院 | 一种基于矿物填充的缝洞型储层导电模型构建方法 |
CN117421890B (zh) * | 2023-10-19 | 2024-05-10 | 重庆科技学院 | 一种基于矿物填充的缝洞型储层导电模型构建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113989433B (zh) | 2022-06-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107701180B (zh) | 一种基于密闭取心的原始油藏含水饱和度计算方法 | |
RU2009135608A (ru) | Способ генерирования численных псевдокернов с использованием изображений скважины, цифровых образов породы и многоточечной статистики | |
CN104616350A (zh) | 缝洞型碳酸盐油藏三维物理模型建立方法 | |
CN106285664B (zh) | 基于逾渗网络模拟的双重介质储层岩石含水饱和度计算法 | |
CN102540265B (zh) | 一种基于网络模拟的多孔介质含水饱和度计算方法 | |
CN104712329A (zh) | 一种泥页岩油气饱和度的计算模型 | |
CN107045147A (zh) | 一种探查河道溶洞的多维探测装置及方法 | |
CN106401574B (zh) | 一种钻前高温地热井地层压力的预测方法 | |
CN106570262B (zh) | 一种储层构型结构的描述方法 | |
CN111911142B (zh) | 裂缝性致密砂岩气藏含水饱和度模型数字岩心构建方法 | |
CN206960673U (zh) | 一种探查河道溶洞的多维探测装置 | |
CN108694264A (zh) | 一种确定页岩气储层渗透率的方法 | |
CN106503295B (zh) | 一种利用状态空间模型解释油田水淹层的方法及装置 | |
CN113989433B (zh) | 一种基于孔隙类型细分的缝洞型储层饱和度模型建立方法 | |
CN106202763B (zh) | 双重介质储层岩石含水饱和度的计算方法 | |
CN112394404A (zh) | 一种渐进式储层精细表征方法 | |
CN109211745B (zh) | 一种富含有机质泥页岩孔隙度的演化过程的恢复方法 | |
CN112145165B (zh) | 一种微裂缝-孔隙型储层动静态渗透率转换方法 | |
CN111963159A (zh) | 一种砾岩致密油储层流体性质的识别方法 | |
CN109763812B (zh) | 基于摩尔-库伦理论的水合物开采力学稳定性定量分析方法 | |
Bordbar et al. | A pseudo-TOF based streamline tracing for streamline simulation method in heterogeneous hydrocarbon reservoirs | |
CN103615239B (zh) | 一种测井地层成像系统及其方法 | |
CN114037805B (zh) | 一种非常规储层饱和度模型构建方法 | |
CN109582923A (zh) | 一种基于双参数计算水合物饱和度的方法及处理终端 | |
CN111862306A (zh) | 一种基于ct图像的孔隙填充型水合物沉积物三维建模方法 |
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 |