CN104729972B - 确定碳酸盐岩气水相对渗透率的方法及装置 - Google Patents
确定碳酸盐岩气水相对渗透率的方法及装置 Download PDFInfo
- Publication number
- CN104729972B CN104729972B CN201510167334.2A CN201510167334A CN104729972B CN 104729972 B CN104729972 B CN 104729972B CN 201510167334 A CN201510167334 A CN 201510167334A CN 104729972 B CN104729972 B CN 104729972B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msup
- msubsup
- mfrac
- 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
Landscapes
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种确定碳酸盐岩气水相对渗透率的方法及装置,该方法包括:确定碳酸盐岩中孔隙的分布;确定各孔隙的气相及水相流速;确定各孔隙的气相及水相流量;根据各孔隙的气相及水相流速、气相及水相流量,确定碳酸盐岩气水相对渗透率。本发明无需人工参与,也无需采用水或者汞等的注入过程,因此能够解决碳酸盐岩气水相对渗透率确定时出现的时间长、误差大的不足。
Description
技术领域
本发明涉及油田开采技术领域,尤其涉及一种确定碳酸盐岩气水相对渗透率的方法及装置。
背景技术
对于碳酸盐岩储层的孔洞、裂缝发育的研究中,目前对碳酸盐岩储层的孔、洞、缝微观流动研究的方法和认识较少,主要采用常规的驱替实验、压汞法等手段,很难精确,快速求取岩心的相对渗透率。下面对驱替实验、压汞法进行说明。
一、流动规律研究中的微观分析技术之一:驱替实验。
驱替实验是将岩心饱和盐水,然后装入岩芯夹持器,设置好流量,在所设置的流量下测定岩芯入口端的压力,并用量杯和秒表测其出口端流量;待压力表数据趋于稳定,出口流速稳定时,记录此时的压力和流量,然后带入达西公式,求得岩心渗透率。此方法简单,成本低廉,是目前储层渗流规律研究中的常用方法。该分析方法的不足是:时间长,实验过程容易出现误差。具体如下:
1、利用驱替实验测定岩心渗透率,需要专业技术人士去判读,对操作人员的技术水平经验依赖性大,仪器的精度会直接影响驱替结果的好坏;实际工作中,因操作人员、摄像参数设置等的不同,分析结果可能存在差别;
2、测定时间长。由于碳酸盐岩连通性差、非均质性很强,做一组驱替实验需要花费比较长的时间,如出现差错,需要重新测定,同时,由于渗透率低,流量小,计量不准也会造成测试的误差;
3、驱替过程中,只能记录一些实时数据,时间和流量,不能实时动态的观察驱替的过程,反映不出真实的流体的流动规律。
二、储层研究中的微观分析技术之一:压汞法。
常规压汞法以毛细管束模型为基础,假设多孔介质由直径大小不同的毛细管束组成。将非润湿相的汞注入被抽空的岩样空间中去时,一定要克服岩石孔隙系统对汞的毛细管阻力。汞作为非润湿相,汞的注入过程中汞首先进入较大的孔隙喉道中;随着注入压力不断增加,汞即不断进入较小的孔隙。岩样的毛细管压力与孔径分布、孔隙体积分布、孔喉连通关系、孔隙度、渗透率及流体饱和度等都有关。因此,利用毛管曲线可以分析不同渗透率岩心的孔喉分布及渗透率贡献率分布特征,探讨渗透率、平均喉道半径、分选系数等孔隙结构参数。
压汞法利用毛管曲线可以分析不同渗透率岩心的孔喉分布及渗透率贡献率分布特征,但有其不足在于:
1、压汞实验过程中,使用的汞有毒,会对人体和环境造成危害;
2、压汞法得到的岩石的孔隙结构特征,是一个平均参数,不能得到不同大小的孔喉的定量分布特征;
3、压汞实验时,汞的注入过程中会对原来的孔隙及其内部的粘土矿物等产生破坏,反映不出真实的岩石孔隙结构特征。注入压力越高,这种影响越大;
4、压汞实验时,只能记录一些实时数据,不能实时动态的观察驱替的过程,反映不出真实的流体的流动规律。
发明内容
本发明实施例提供一种确定碳酸盐岩气水相对渗透率的方法,用以解决现有技术在确定碳酸盐岩气水相对渗透率时时间长、误差大的不足,该方法包括:
确定碳酸盐岩中孔隙的分布;
确定各孔隙的气相及水相流速;
确定各孔隙的气相及水相流量;
根据各孔隙的气相及水相流速、气相及水相流量,确定碳酸盐岩气水相对渗透率;
所述碳酸盐岩中孔隙的分布根据多孔介质分形几何原理确定,包括:
取特征长度为L0的正方形岩样,则毛细管束直线长度为L0,单元截面面积为A,单元截面上孔隙半径在r和r+dr之间的孔隙个数为:
孔隙大小分布的概率密度函数为:
多孔介质毛细管实际弯曲长度为:
所述确定各孔隙的气相及水相流速,包括:根据以下方式确定各孔隙的气相及水相流速:
其中:Df为孔隙分形维数;DT为迂曲度分形维数;下标:max表示最大;min表示最小;r为毛细管中心到任意位置的径向距离,μm;dr为径向距离增量,μm;r1为气、水交界面相半径,μm;r0为毛细管半径,μm;μw为水相黏度,mPa·s;μg为气相黏度,mPa·s;p1为入口端压力,MPa;p2为出口端压力,MPa;δ为束缚水膜厚度,μm;下标:w为水相;g为气相;
所述确定各孔隙的气相及水相流量,包括:根据达西公式按以下方式确定各孔隙的气相及水相流量为:
其中,Qw为水相总流量,cm3/s;Qg为气相总流量,cm3/s;A为单元截面面积,cm2,K为岩石绝对渗透率,μm2;Krw为水的相对渗透率,百分数;Krg为气的相对渗透率,百分数。
一个实施例中,所述根据各孔隙的气相及水相流速、气相及水相流量,确定碳酸盐岩气水相对渗透率,包括:根据以下方式确定碳酸盐岩气水相对渗透率为:
其中:Sg为气相饱和度,百分数。
一个实施例中,所述确定碳酸盐岩中孔隙的分布,包括:通过CT扫描获取的数据确定碳酸盐岩中孔隙的分布。
本发明实施例还提供一种确定碳酸盐岩气水相对渗透率的装置,用以解决现有技术在确定碳酸盐岩气水相对渗透率时时间长、误差大的不足,该装置包括:
孔隙确定模块,用于确定碳酸盐岩中孔隙的分布;
流速确定模块,用于确定各孔隙的气相及水相流速;
流量确定模块,用于确定各孔隙的气相及水相流量;
渗透率确定模块,用于根据各孔隙的气相及水相流速、气相及水相流量,确定碳酸盐岩气水相对渗透率;
所述孔隙确定模块进一步用于按如下方式根据多孔介质分形几何原理确定碳酸盐岩中孔隙的分布:
取特征长度为L0的正方形岩样,则毛细管束直线长度为L0,单元截面面积为A,单元截面上孔隙半径在r和r+dr之间的孔隙个数为:
孔隙大小分布的概率密度函数为:
多孔介质毛细管实际弯曲长度为:
所述流速确定模块进一步用于根据以下方式确定各孔隙的气相及水相流速:
其中:Df为孔隙分形维数;DT为迂曲度分形维数;下标:max表示最大;min表示最小;r为毛细管中心到任意位置的径向距离,μm;dr为径向距离增量,μm;r1为气、水交界面相半径,μm;r0为毛细管半径,μm;μw为水相黏度,mPa·s;μg为气相黏度,mPa·s;p1为入口端压力,MPa;p2为出口端压力,MPa;δ为束缚水膜厚度,μm;下标:w为水相;g为气相;
所述流量确定模块进一步用于根据达西公式按以下方式确定各孔隙的气相及水相流量为:
其中,Qw为水相总流量,cm3/s;Qg为气相总流量,cm3/s;A为单元截面面积,cm2,K为岩石绝对渗透率,μm2;Krw为水的相对渗透率,百分数;Krg为气的相对渗透率,百分数。
一个实施例中,所述渗透率确定模块进一步用于根据以下方式确定碳酸盐岩气水相对渗透率为:
其中:Sg为气相饱和度,百分数。
一个实施例中,所述孔隙确定模块进一步用于通过CT扫描获取的数据确定碳酸盐岩中孔隙的分布。
在本发明实施例提供的技术方案中,能够在确定碳酸盐岩中孔隙的分布后,确定各孔隙的气相及水相流速,确定各孔隙的气相及水相流量,进而确定碳酸盐岩气水相对渗透率,无需人工参与,也无需采用水或者汞等的注入过程,因此能够解决碳酸盐岩气水相对渗透率确定时出现的时间长、误差大的不足。
进一步的,由于确定碳酸盐岩中孔隙的分布是通过CT扫描进行确定的,因此可以在岩心无损的情况下,可以准确的分析出岩心的孔喉分布,进一步的减少了确定碳酸盐岩气水相对渗透率过程中的时间,减小了误差。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本发明实施例的确定碳酸盐岩气水相对渗透率的方法实施流程示意图;
图2为本发明实施例的确定碳酸盐岩气水相对渗透率的方法实例图;
图3为本发明实施例的确定碳酸盐岩气水相对渗透率的装置结构示意图;
图4为本发明实施例的新疆某碳酸盐岩气藏岩心CT扫描得到的岩心孔隙半径频率分布示意图;
图5为本发明实施例的岩心A的孔喉直径分布特征示意图;
图6为本发明实施例的岩心A的气水相渗曲线示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白,下面结合附图对本发明实施例做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。
针对目前常规岩心及非常规岩心流动模拟时间长,误差大的不足,在本发明实施例中提供的技术方案中,通过图像处理技术,结合程序设计,利用CT扫描图像,以碳酸盐岩CT扫描图像为实例,通过分析岩心孔喉分布,推导了碳酸盐岩渗透率计算的方案,计算了气水的相对渗透率,为碳酸盐岩油藏开发提供了可靠的方法和技术。具体实施可以如下:
图1为确定碳酸盐岩气水相对渗透率的方法实施流程示意图,如图1所示,可以包括如下步骤:
步骤101、确定碳酸盐岩中孔隙的分布;
步骤102、确定各孔隙的气相及水相流速;
步骤103、确定各孔隙的气相及水相流量;
步骤104、根据各孔隙的气相及水相流速、气相及水相流量,确定碳酸盐岩气水相对渗透率。
下面对各步骤的具体实施方式进行说明。
步骤101中,在确定碳酸盐岩中孔隙的分布时,首先确定分形维数和孔隙密度分布函数。
为了研究碳酸盐岩岩心孔隙分布,可以通过CT扫描,对碳酸盐岩岩心进行孔隙分布测定。将测试所得的岩心孔隙半径按从小到大顺序进行排列,绘制岩心孔隙半径r和大于该孔隙半径的孔隙数目N(R>r)关系曲线图。
根据多孔介质分形几何原理,取特征长度为L0的正方形岩样,则毛细管束直线长度为L0,单元截面面积为A,单元截面上孔隙半径在r和r+dr之间的孔隙个数为:
孔隙大小分布的概率密度函数为:
多孔介质毛细管实际弯曲长度为:
式中:Df为孔隙分形维数;DT为迂曲度分形维数;下标:max表示最大;min表示最小;r为毛细管中心到任意位置的径向距离,μm;dr为径向距离增量,μm。
步骤102中,在确定各孔隙的气相及水相流速时,具体可以如下:
气、水两相在毛细管内发生层流流动,可以把流体流动看成一组长度相同、速度不同的同心液筒,可以得到水相、气相速度方程分别为:
vw=0,r0-δ≤r≤r0 (2)
式中:v为流体在毛细管内流动速度,cm/s;L为毛细管长度,cm;r为毛细管中心到任意位置的径向距离,μm;r1为气、水交界面相半径,μm;r0为毛细管半径,μm;μw为水相黏度,mPa·s;μg为气相黏度,mPa·s;p1为入口端压力,MPa;p2为出口端压力,MPa;c为常数;δ为束缚水膜厚度,μm;下标:w为水相;g为气相。
在气、水交界面处:
vw=vg,r=r1 (4)
将式(2)、式(4)代入式(1)和式(3),可得:
步骤103中,在确定各孔隙的气相及水相流量时,可以如下:
通过半径为r0毛细管的气相及水相流量分别为:
式中:q为流体在毛细管内流量,cm3/s。
对式(7)和式(8)进行积分可以得到岩石气、水总流量分别为:
式中:Qw为水相总流量,cm3/s;Qg为气相总流量,cm3/s。
岩石含水饱和度、束缚水饱和度分别为:
Sw=1-Sg (14)
式中:Sw为岩石含水饱和度,百分数;Sg为气相饱和度,百分数;Swc为岩石束缚水饱和度,百分数。
由达西公式可知通过岩石水相、气相流量分别为:
式中:A为单元截面面积,cm2,K为岩石绝对渗透率,μm2;Krw为水的相对渗透率,百分数;Krg为气的相对渗透率,百分数。
步骤104中,在确定碳酸盐岩气水相对渗透率时,气水相渗曲线表达式为:
图2为确定碳酸盐岩气水相对渗透率的方法实例图。如图2所示,本例中,可以先读取多孔介质二值化图像,例如可以通过CT扫描获取的数据获得多孔介质二值化图像;然后分析岩心孔喉结构,求取孔径大小,再计算出孔隙密度分布函数和分形维数;接着计算水相速度和气相速度,计算水相流量和气相流量;最后计算出气水相对渗透率。
基于同一发明构思,本发明实施例中还提供了一种确定碳酸盐岩气水相对渗透率的装置,由于该装置解决问题的原理与一种确定碳酸盐岩气水相对渗透率的方法相似,因此该装置的实施可以参见方法的实施,重复之处不再赘述。
图3为确定碳酸盐岩气水相对渗透率的装置结构示意图,如图3所示,装置中包括:
孔隙确定模块301,用于确定碳酸盐岩中孔隙的分布;
流速确定模块302,用于确定各孔隙的气相及水相流速;
流量确定模块303,用于确定各孔隙的气相及水相流量;
渗透率确定模块304,用于根据各孔隙的气相及水相流速、气相及水相流量,确定碳酸盐岩气水相对渗透率。
实施中,孔隙确定模块还可以进一步用于根据多孔介质分形几何原理确定碳酸盐岩中孔隙的分布。
实施中,孔隙确定模块还可以进一步用于按如下方式根据多孔介质分形几何原理确定碳酸盐岩中孔隙的分布:
取特征长度为L0的正方形岩样,则毛细管束直线长度为L0,单元截面面积为A,单元截面上孔隙半径在r和r+dr之间的孔隙个数为:
孔隙大小分布的概率密度函数为:
多孔介质毛细管实际弯曲长度为:
其中:Df为孔隙分形维数;DT为迂曲度分形维数;下标:max表示最大;min表示最小;r为毛细管中心到任意位置的径向距离,μm;dr为径向距离增量,μm。
实施中,流速确定模块还可以进一步用于根据以下方式确定各孔隙的气相及水相流速:
其中:v为流体在毛细管内流动速度,cm/s;r1为气、水交界面相半径,μm;r0为毛细管半径,μm;μw为水相黏度,mPa·s;μg为气相黏度,mPa·s;p1为入口端压力,MPa;p2为出口端压力,MPa;δ为束缚水膜厚度,μm;下标:w为水相;g为气相。
实施中,流量确定模块还可以进一步用于根据达西公式按以下方式确定各孔隙的气相、水相流量为:
其中,Qw为水相总流量,cm3/s;Qg为气相总流量,cm3/s;A为单元截面面积,cm2,K为岩石绝对渗透率,μm2;Krw为水的相对渗透率,百分数;Krg为气的相对渗透率,百分数。
实施中,渗透率确定模块还可以进一步用于根据以下方式确定碳酸盐岩气水相对渗透率为:
其中:Sg为气相饱和度,百分数。
实施中,孔隙确定模块还可以进一步用于通过CT扫描获取的数据确定碳酸盐岩中孔隙的分布。
下面以实例说明。
图4为新疆某碳酸盐岩气藏岩心CT扫描得到的岩心孔隙半径频率分布示意图,如图4所示,以新疆某碳酸盐岩气藏岩心为例,通过电镜扫描技术得到岩心孔隙半径频率分布。
图5为岩心A的孔喉直径分布特征示意图,如图5所示,根据图4的扫描图获得孔喉直径分布。
图6为岩心A的气水相渗曲线示意图,如图6所示,根据本发明实施例中提供的方案,可得到如图6所示的岩心气水相渗曲线。
由上述实施例可见,采用本发明实施例提供的技术方案,可以在岩心无损的情况下,基于微纳米CT扫描可以准确的分析出岩心的孔喉分布,进而推导了压汞曲线计算方案,求取了岩心的相对渗透率,解决像碳酸盐岩油藏这样非均质性强、非常规储层岩石的微观流动问题,为非常规储层物性参数计算、有效驱替和调堵提供明确的孔隙特征和动态效果描述,有非常规储层的高效、合理开发提供有力的条件。显然,采用本发明实施例提供的技术方案,可以在岩心无损的情况下,快速、准确的分析出岩心的渗透率,可视化流动模拟的过程,有油田开发提供有力的条件。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种确定碳酸盐岩气水相对渗透率的方法,其特征在于,该方法包括:
确定碳酸盐岩中孔隙的分布;
确定各孔隙的气相及水相流速;
确定各孔隙的气相及水相流量;
根据各孔隙的气相及水相流速、气相及水相流量,确定碳酸盐岩气水相对渗透率;
所述碳酸盐岩中孔隙的分布根据多孔介质分形几何原理确定,包括:
取特征长度为L0的正方形岩样,则毛细管束直线长度为L0,单元截面面积为A,单元截面上孔隙半径在r和r+dr之间的孔隙个数为:
<mrow>
<mi>N</mi>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>/</mo>
<msub>
<mi>r</mi>
<mi>min</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>D</mi>
<mi>f</mi>
</msub>
</msup>
<mo>;</mo>
</mrow>
孔隙大小分布的概率密度函数为:
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>D</mi>
<mi>f</mi>
</msub>
<msubsup>
<mi>r</mi>
<mi>min</mi>
<msub>
<mi>D</mi>
<mi>f</mi>
</msub>
</msubsup>
<msup>
<mi>r</mi>
<mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>D</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mo>;</mo>
</mrow>
多孔介质毛细管实际弯曲长度为:
<mrow>
<mi>L</mi>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
</mrow>
</msup>
<msubsup>
<mi>L</mi>
<mn>0</mn>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
</msubsup>
<mo>;</mo>
</mrow>
所述确定各孔隙的气相及水相流速,包括:根据以下方式确定各孔隙的气相及水相流速:
<mrow>
<msub>
<mi>v</mi>
<mi>w</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>4</mn>
<msub>
<mi>&mu;</mi>
<mi>w</mi>
</msub>
<mi>L</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>&delta;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>4</mn>
<msub>
<mi>&mu;</mi>
<mi>w</mi>
</msub>
<mi>L</mi>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>v</mi>
<mi>g</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>4</mn>
<msub>
<mi>&mu;</mi>
<mi>g</mi>
</msub>
<mi>L</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<msubsup>
<mi>r</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
</mrow>
<mrow>
<mn>4</mn>
<msub>
<mi>&mu;</mi>
<mi>g</mi>
</msub>
<mi>L</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>&delta;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msubsup>
<mi>r</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>4</mn>
<msub>
<mi>&mu;</mi>
<mi>w</mi>
</msub>
<mi>L</mi>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
其中:Df为孔隙分形维数;DT为迂曲度分形维数;下标:max表示最大;min表示最小;r为毛细管中心到任意位置的径向距离,μm;dr为径向距离增量,μm;r1为气、水交界面相半径,μm;r0为毛细管半径,μm;μw为水相黏度,mPa·s;μg为气相黏度,mPa·s;p1为入口端压力,MPa;p2为出口端压力,MPa;δ为束缚水膜厚度,μm;下标:w为水相;g为气相;
所述确定各孔隙的气相及水相流量,包括:根据达西公式按以下方式确定各孔隙的气相及水相流量为:
<mrow>
<msub>
<mi>Q</mi>
<mi>w</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>AKK</mi>
<mrow>
<mi>r</mi>
<mi>w</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&mu;</mi>
<mi>w</mi>
</msub>
<msub>
<mi>L</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>Q</mi>
<mi>g</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>AKK</mi>
<mrow>
<mi>r</mi>
<mi>g</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&mu;</mi>
<mi>g</mi>
</msub>
<msub>
<mi>L</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
其中,Qw为水相总流量,cm3/s;Qg为气相总流量,cm3/s;A为单元截面面积,cm2,K为岩石绝对渗透率,μm2;Krw为水的相对渗透率,百分数;Krg为气的相对渗透率,百分数。
2.如权利要求1所述的方法,其特征在于,所述根据各孔隙的气相及水相流速、气相及水相流量,确定碳酸盐岩气水相对渗透率,包括:根据以下方式确定碳酸盐岩气水相对渗透率为:
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>r</mi>
<mi>w</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msup>
<mn>2</mn>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<msubsup>
<mi>L</mi>
<mn>0</mn>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>r</mi>
<mi>min</mi>
</msub>
<msub>
<mi>r</mi>
<mi>max</mi>
</msub>
</msubsup>
<msup>
<mi>r</mi>
<mrow>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msup>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mi>r</mi>
<mo>-</mo>
<mi>&delta;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
<msub>
<mi>S</mi>
<mi>g</mi>
</msub>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
<mi>d</mi>
<mi>r</mi>
</mrow>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>r</mi>
<mi>min</mi>
</msub>
<msub>
<mi>r</mi>
<mi>max</mi>
</msub>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mi>r</mi>
<mo>-</mo>
<mi>&delta;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
<msup>
<mi>r</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>f</mi>
</msub>
</mrow>
</msup>
<mi>d</mi>
<mi>r</mi>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>r</mi>
<mi>g</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msubsup>
<mi>S</mi>
<mi>g</mi>
<mn>2</mn>
</msubsup>
<mrow>
<msup>
<mn>2</mn>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<msubsup>
<mi>L</mi>
<mn>0</mn>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>r</mi>
<mi>min</mi>
</msub>
<msub>
<mi>r</mi>
<mi>max</mi>
</msub>
</msubsup>
<msup>
<mi>r</mi>
<mrow>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msup>
<mi>d</mi>
<mi>r</mi>
</mrow>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>r</mi>
<mi>min</mi>
</msub>
<msub>
<mi>r</mi>
<mi>max</mi>
</msub>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mi>r</mi>
<mo>-</mo>
<mi>&delta;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
<msup>
<mi>r</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
</mrow>
</msup>
<mi>d</mi>
<mi>r</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>S</mi>
<mi>g</mi>
</msub>
<msub>
<mi>&mu;</mi>
<mi>g</mi>
</msub>
</mrow>
<mrow>
<msup>
<mn>2</mn>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<msubsup>
<mi>L</mi>
<mn>0</mn>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<msub>
<mi>&mu;</mi>
<mi>w</mi>
</msub>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>r</mi>
<mi>min</mi>
</msub>
<msub>
<mi>r</mi>
<mi>max</mi>
</msub>
</msubsup>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mi>r</mi>
<mo>-</mo>
<mi>&delta;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msub>
<mi>S</mi>
<mi>g</mi>
</msub>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<msup>
<mi>r</mi>
<mrow>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>f</mi>
</msub>
</mrow>
</msup>
<mi>d</mi>
<mi>r</mi>
</mrow>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>r</mi>
<mi>min</mi>
</msub>
<msub>
<mi>r</mi>
<mi>max</mi>
</msub>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mi>r</mi>
<mo>-</mo>
<mi>&delta;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
<msup>
<mi>r</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
</mrow>
</msup>
<mi>d</mi>
<mi>r</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
其中:Sg为气相饱和度,百分数。
3.如权利要求1或2所述的方法,其特征在于,所述确定碳酸盐岩中孔隙的分布,包括:通过CT扫描获取的数据确定碳酸盐岩中孔隙的分布。
4.一种确定碳酸盐岩气水相对渗透率的装置,其特征在于,包括:
孔隙确定模块,用于确定碳酸盐岩中孔隙的分布;
流速确定模块,用于确定各孔隙的气相及水相流速;
流量确定模块,用于确定各孔隙的气相及水相流量;
渗透率确定模块,用于根据各孔隙的气相及水相流速、气相及水相流量,确定碳酸盐岩气水相对渗透率;
所述孔隙确定模块进一步用于按如下方式根据多孔介质分形几何原理确定碳酸盐岩中孔隙的分布:
取特征长度为L0的正方形岩样,则毛细管束直线长度为L0,单元截面面积为A,单元截面上孔隙半径在r和r+dr之间的孔隙个数为:
<mrow>
<mi>N</mi>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>/</mo>
<msub>
<mi>r</mi>
<mi>min</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>D</mi>
<mi>f</mi>
</msub>
</msup>
<mo>;</mo>
</mrow>
孔隙大小分布的概率密度函数为:
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>D</mi>
<mi>f</mi>
</msub>
<msubsup>
<mi>r</mi>
<mi>min</mi>
<msub>
<mi>D</mi>
<mi>f</mi>
</msub>
</msubsup>
<msup>
<mi>r</mi>
<mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>D</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mo>;</mo>
</mrow>
多孔介质毛细管实际弯曲长度为:
<mrow>
<mi>L</mi>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
</mrow>
</msup>
<msubsup>
<mi>L</mi>
<mn>0</mn>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
</msubsup>
<mo>;</mo>
</mrow>
所述流速确定模块进一步用于根据以下方式确定各孔隙的气相及水相流速:
<mrow>
<msub>
<mi>v</mi>
<mi>w</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>4</mn>
<msub>
<mi>&mu;</mi>
<mi>w</mi>
</msub>
<mi>L</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>&delta;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>4</mn>
<msub>
<mi>&mu;</mi>
<mi>w</mi>
</msub>
<mi>L</mi>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>v</mi>
<mi>g</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>4</mn>
<msub>
<mi>&mu;</mi>
<mi>g</mi>
</msub>
<mi>L</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<msubsup>
<mi>r</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
</mrow>
<mrow>
<mn>4</mn>
<msub>
<mi>&mu;</mi>
<mi>g</mi>
</msub>
<mi>L</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>&delta;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msubsup>
<mi>r</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mrow>
<mrow>
<mn>4</mn>
<msub>
<mi>&mu;</mi>
<mi>w</mi>
</msub>
<mi>L</mi>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
其中:Df为孔隙分形维数;DT为迂曲度分形维数;下标:max表示最大;min表示最小;r为毛细管中心到任意位置的径向距离,μm;dr为径向距离增量,μm;r1为气、水交界面相半径,μm;r0为毛细管半径,μm;μw为水相黏度,mPa·s;μg为气相黏度,mPa·s;p1为入口端压力,MPa;p2为出口端压力,MPa;δ为束缚水膜厚度,μm;下标:w为水相;g为气相;
所述流量确定模块进一步用于根据达西公式按以下方式确定各孔隙的气相及水相流量为:
<mrow>
<msub>
<mi>Q</mi>
<mi>w</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>AKK</mi>
<mrow>
<mi>r</mi>
<mi>w</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&mu;</mi>
<mi>w</mi>
</msub>
<msub>
<mi>L</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>Q</mi>
<mi>g</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>AKK</mi>
<mrow>
<mi>r</mi>
<mi>g</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&mu;</mi>
<mi>g</mi>
</msub>
<msub>
<mi>L</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
其中,Qw为水相总流量,cm3/s;Qg为气相总流量,cm3/s;A为单元截面面积,cm2,K为岩石绝对渗透率,μm2;Krw为水的相对渗透率,百分数;Krg为气的相对渗透率,百分数。
5.如权利要求4所述的装置,其特征在于,所述渗透率确定模块进一步用于根据以下方式确定碳酸盐岩气水相对渗透率为:
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>r</mi>
<mi>w</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msup>
<mn>2</mn>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<msubsup>
<mi>L</mi>
<mn>0</mn>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>r</mi>
<mi>min</mi>
</msub>
<msub>
<mi>r</mi>
<mi>max</mi>
</msub>
</msubsup>
<msup>
<mi>r</mi>
<mrow>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msup>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mi>r</mi>
<mo>-</mo>
<mi>&delta;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
<msub>
<mi>S</mi>
<mi>g</mi>
</msub>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
<mi>d</mi>
<mi>r</mi>
</mrow>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>r</mi>
<mi>min</mi>
</msub>
<msub>
<mi>r</mi>
<mi>max</mi>
</msub>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mi>r</mi>
<mo>-</mo>
<mi>&delta;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
<msup>
<mi>r</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>f</mi>
</msub>
</mrow>
</msup>
<mi>d</mi>
<mi>r</mi>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>r</mi>
<mi>g</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msubsup>
<mi>S</mi>
<mi>g</mi>
<mn>2</mn>
</msubsup>
<mrow>
<msup>
<mn>2</mn>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<msubsup>
<mi>L</mi>
<mn>0</mn>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>r</mi>
<mi>min</mi>
</msub>
<msub>
<mi>r</mi>
<mi>max</mi>
</msub>
</msubsup>
<msup>
<mi>r</mi>
<mrow>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msup>
<mi>d</mi>
<mi>r</mi>
</mrow>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>r</mi>
<mi>min</mi>
</msub>
<msub>
<mi>r</mi>
<mi>max</mi>
</msub>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mi>r</mi>
<mo>-</mo>
<mi>&delta;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
<msup>
<mi>r</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
</mrow>
</msup>
<mi>d</mi>
<mi>r</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>S</mi>
<mi>g</mi>
</msub>
<msub>
<mi>&mu;</mi>
<mi>g</mi>
</msub>
</mrow>
<mrow>
<msup>
<mn>2</mn>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<msubsup>
<mi>L</mi>
<mn>0</mn>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<msub>
<mi>&mu;</mi>
<mi>w</mi>
</msub>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>r</mi>
<mi>min</mi>
</msub>
<msub>
<mi>r</mi>
<mi>max</mi>
</msub>
</msubsup>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mi>r</mi>
<mo>-</mo>
<mi>&delta;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msub>
<mi>S</mi>
<mi>g</mi>
</msub>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<msup>
<mi>r</mi>
<mrow>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>f</mi>
</msub>
</mrow>
</msup>
<mi>d</mi>
<mi>r</mi>
</mrow>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>r</mi>
<mi>min</mi>
</msub>
<msub>
<mi>r</mi>
<mi>max</mi>
</msub>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mi>r</mi>
<mo>-</mo>
<mi>&delta;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
<msup>
<mi>r</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>T</mi>
</msub>
</mrow>
</msup>
<mi>d</mi>
<mi>r</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
其中:Sg为气相饱和度,百分数。
6.如权利要求4或5所述的装置,其特征在于,所述孔隙确定模块进一步用于通过CT扫描获取的数据确定碳酸盐岩中孔隙的分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510167334.2A CN104729972B (zh) | 2015-04-10 | 2015-04-10 | 确定碳酸盐岩气水相对渗透率的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510167334.2A CN104729972B (zh) | 2015-04-10 | 2015-04-10 | 确定碳酸盐岩气水相对渗透率的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104729972A CN104729972A (zh) | 2015-06-24 |
CN104729972B true CN104729972B (zh) | 2017-10-27 |
Family
ID=53454069
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510167334.2A Active CN104729972B (zh) | 2015-04-10 | 2015-04-10 | 确定碳酸盐岩气水相对渗透率的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104729972B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105424575B (zh) * | 2015-11-10 | 2019-06-11 | 中国石油天然气股份有限公司 | 一种确定渗透率的方法和装置 |
CN105717025A (zh) * | 2016-02-24 | 2016-06-29 | 西南石油大学 | 一种岩石有效渗透率的测定装置及其使用方法 |
CN105699268A (zh) * | 2016-02-24 | 2016-06-22 | 西南石油大学 | 一种岩石绝对渗透率的测定装置及其使用方法 |
CN105628587A (zh) * | 2016-03-09 | 2016-06-01 | 西南石油大学 | 一种气水两相有效渗透率精确测定新装置及方法 |
CN106680172B (zh) * | 2016-11-16 | 2019-09-27 | 中国石油大学(北京) | 评价致密油储层的裂缝的方法 |
CN106979918B (zh) * | 2017-04-20 | 2020-02-07 | 中国石油大学(北京) | 一种获取致密油藏岩心的液体渗透率的方法及装置 |
CN108986627B (zh) * | 2018-06-13 | 2021-01-01 | 中国石油天然气股份有限公司 | 一种人造岩心微观可视化模型及其制备方法和应用 |
CN111521532B (zh) * | 2019-02-01 | 2023-02-10 | 中国石油天然气股份有限公司 | 一种致密储层含气饱和度的测定方法 |
CN110146940B (zh) * | 2019-04-04 | 2021-03-26 | 中国石油大学(华东) | 一种气藏束缚水饱和度、气藏可动水范围确定方法和装置 |
CN112284991B (zh) * | 2019-07-22 | 2022-09-13 | 中海石油(中国)有限公司 | 一种基于分形理论的孔隙型多孔介质渗透率预测方法 |
CN110865011B (zh) * | 2019-11-14 | 2021-03-09 | 西南石油大学 | 一种基于数字成像技术计算致密岩心相对渗透率的方法 |
CN112345429B (zh) * | 2020-11-02 | 2021-08-06 | 北京科技大学 | 非常规油气藏非均匀地质参数的获取方法及装置 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1015661B (zh) * | 1988-12-19 | 1992-02-26 | 大庆石油管理局勘探开发研究院 | 用于测定油藏岩心油、气、水三相相对渗透率实验方法及装置 |
CN103573234B (zh) * | 2012-08-06 | 2016-06-08 | 中国石油化工股份有限公司 | 一种确定完整油水相对渗透率曲线的方法 |
CN104316449A (zh) * | 2014-11-05 | 2015-01-28 | 辽宁石油化工大学 | 一种用于测定火山岩气、水相对渗透率的实验方法及实验装置 |
CN104359819A (zh) * | 2014-11-10 | 2015-02-18 | 中国石油天然气股份有限公司 | 低渗致密岩心气水相对渗透率测定装置及测定方法 |
-
2015
- 2015-04-10 CN CN201510167334.2A patent/CN104729972B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104729972A (zh) | 2015-06-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104729972B (zh) | 确定碳酸盐岩气水相对渗透率的方法及装置 | |
US10161891B1 (en) | Method for characterizing rock physical characteristics of deeply buried carbonate rocks | |
CN108896599B (zh) | 一种测试气水相对渗透率曲线的系统及方法 | |
CN111624147B (zh) | 岩心的相对渗透率测定方法及装置 | |
CN104568694B (zh) | 一种致密岩心气‑水相对渗透率的测试方法 | |
EP3593159B1 (en) | Absolute porosity and pore size determination of pore types in media with varying pore sizes using nmr | |
EP2765409B1 (en) | Nuclear magnetic resonance rock sample analysis method and instrument with constant gradient field | |
CN103954544B (zh) | 一种聚合物控水增气效果评价的实验装置和实验方法 | |
CN104237099B (zh) | 测定致密岩心径向渗透率的装置及方法 | |
CN104634804B (zh) | 一种利用核磁共振t2谱确定储层相对渗透率的方法 | |
CN106525690A (zh) | 致密砂岩稳态法气水相对渗透率曲线测定方法 | |
CN103513270B (zh) | 一种基于岩石声学特性的气层识别评价方法及装置 | |
CN110927035A (zh) | 一种低渗致密砂岩束缚水饱和度计算方法 | |
CN106979918A (zh) | 一种获取致密油藏岩心的液体渗透率的方法及装置 | |
CN105891089B (zh) | 储层渗透率确定方法及装置 | |
CN104359819A (zh) | 低渗致密岩心气水相对渗透率测定装置及测定方法 | |
CN104316449A (zh) | 一种用于测定火山岩气、水相对渗透率的实验方法及实验装置 | |
CN103498669A (zh) | 一种非均质岩心模型层间窜流量的定量测定方法 | |
CN112098155B (zh) | 油藏油水岩反应实验装置、方法以及取样位置的确定方法 | |
CN106777515B (zh) | 一种基于岩心实验资料分析致密气井产能的方法 | |
CN109100812A (zh) | 基于核磁共振的岩石孔隙分形维数评价方法及装置 | |
CN206161492U (zh) | 一种可实现变水压力作用的渗透装置 | |
CN104374679B (zh) | 一种确定储层渗透率贡献值的方法 | |
CN105676309A (zh) | 裂缝性油藏逆向渗吸采收率预测方法 | |
CN110441209A (zh) | 一种基于致密储层数字岩心计算岩石渗透率的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |