CN104819923A - 基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法 - Google Patents
基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法 Download PDFInfo
- Publication number
- CN104819923A CN104819923A CN201510251449.XA CN201510251449A CN104819923A CN 104819923 A CN104819923 A CN 104819923A CN 201510251449 A CN201510251449 A CN 201510251449A CN 104819923 A CN104819923 A CN 104819923A
- Authority
- CN
- China
- Prior art keywords
- nuclear
- magnetism
- type
- throat distribution
- free fluid
- 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)
Abstract
本发明公开了一种基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法,包括以下步骤:第一、挑选基础样品;第二、划分孔隙结构类型;第三、构建识别孔隙结构的测井图版;第四、分孔隙结构类型,建立孔喉分布反演公式;第五、验证公式的正确性;第六、利用测井图版识别井筒剖面的孔隙结构类型;第七、利用核磁共振测井T2谱定量反演自由流体孔喉分布。本发明的有益效果如下:提供定量评价孔隙结构的一种可靠方法和技术,打破现行方法的研究思想。沿井眼剖面,通过核磁共振测井反演出连续的自由流体孔喉分布,为低渗透砂岩储层有效性评价提供了直接依据,也是测井信息用于定量反演储层微观孔隙结构的有益探索,推动了核磁测井技术的应用与发展。
Description
技术领域
本发明涉及岩心核磁共振、压汞实验及核磁共振测井领域,特别涉及一种基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法。
背景技术
核磁共振的基础是原子核的磁性及其外加磁场的相互作用。在外加磁场作用下使得地层流体中的质子极化,核自旋空间取向从无序变为有序。极化后,关闭外加磁场,测量自旋***从不平衡状态向平衡状态过程的时间(弛豫时间),即测量纵向弛豫时间T1、横向弛豫时间T2。
磁共振测井的原始数据是随时间衰减的自旋回波串,自旋回波串包含了储层物性、孔隙类型、孔径大小、流体类型及其分布等十分丰富的信息。自旋一回波串可以表示为多个衰减指数的迭加:
式中:Echo(t),核磁测井测量得到的随时间t变化的回波幅度;
T2i,第i个组分的横向弛豫时间,i=1,…,N;
φi,对应于T2i组分的核磁共振孔隙度,i=1,…,N;
通过对核磁测井仪器测量得到的自旋回波串(即式(1))进行多指数拟合,可以求出每一测量点的T2i-φi的分布曲线,即所谓的核磁共振测井T2分布谱。该T2分布谱的横坐标为T2值,纵坐标为信号相对幅度,总的信号幅度与岩石孔隙度有关。
在岩石孔隙中,流体分子弛豫时间不仅与本身弛豫特点有关,还与岩石孔隙结构、成分关系密切,弛豫过程较复杂。目前普遍认为有三种机制控制孔隙流体的横向弛豫过程,即体积弛豫、表面弛豫和扩散弛豫,分别用T2B、T2S、T2D表示相应的弛豫时间常数。进而横向弛豫时间表示为:
式中:T2B为流体的体积(自由)弛豫时间,ms;
T2S为横向表面弛豫时间,ms;
T2D为横向扩散弛豫时间,ms;
D为扩散系数,μm2/ms;
G为磁场梯度,G/cm;
TE为回波间隔,ms;
S为孔隙表面积,cm2;
V为孔隙体积,cm3;
ρ2为岩石的横向表面弛豫强度,μm/ms;
T2B的数值通常为2~3s,比T2大得多。即T2B>>T2。因此式3中右边的第一项可忽略。当磁场均匀时(相对G很小),且TE足够短时,式2右边第三项也可忽略,于是T2为:
得到T2与孔径rc的关系式为:
式中:FS为形状几何因子。对于球形孔隙,FS=3,;对柱状管道,FS=2。由式5可知,孔隙内的流体的弛豫时间与孔隙空间的大小及形状有关,孔隙越小,比表面积越大,表面相互作用的影响就越强烈,T2时间也越短。观测的弛豫时间T2和平均孔径rc是一一对应的,因此,可以利用T2谱来评价岩石孔隙大小及其孔喉分布。
目前现有技术提出横向驰豫时间和毛管压力之间的转换关系,如式5,该方法假设孔隙结构可以简化为球形和柱状管道。其比表面积与孔径成线性关系,Pc=C/T2线性关系进行转换获得毛细管压力曲线,其转换系数C由相似原理确定,该方法的一个关键假设条件是要求核磁T2谱曲线与毛细管压力曲线的形态一致性高,而从大量的岩心核磁共振测量与压汞测量结果对比发现,T2谱与压汞微分曲线形态存在一定差异。并且,未考虑T2谱与压汞曲线反应的孔隙空间不同这一关键。故而,这种方法计算误差就会较大。
在上述的研究基础上,基于最大相似性原理的方法确定T2与Pc之间的转换系数C。首先假设一个C,将与Pc~SHg两条曲线重合在一张图上(图1)。通过调节C值来移动曲线。明显存在唯一值C,是使与Pc~SHg的相关性最大。去其相关性最大时的C值为T2与Pc之间的转换系数,进而将核磁T2谱转换为压汞微分曲线。这一方法的计算精度不高。原因在于忽略了孔喉空间中薄膜束缚水对小孔核磁信号的影响。即核磁T2谱测量过程中,大孔喉空间薄膜束缚水信号叠加到小孔核磁信号上,使得小孔信号增强。因此,核磁T2谱与压汞分布曲线在小孔后部分不成线性关系。故而,在小孔吼部分转换系数并不适用。
现有技术还提出了2种核磁共振分布评价毛细管压力曲线的改进方法。
第1种方法:岩心核磁T2分布与相应的压汞毛管半径分布进行的相关分析。发现在大孔、小孔处两者呈现明显不同两种幂函数关系。所以,在T2分布为双峰时需要用不同幂函数分段构造毛管压力曲线对T2分布呈单峰分布时,采用单一幂函数来构建伪毛管压力曲线。
第2种方法:何雨丹从压汞、NMR核磁所反映的岩石孔隙空间入手,提出薄膜束缚水对核磁T2谱的影响,并提出与压汞信息对应的是核磁自由流体谱而非核磁饱和谱。故而,消除薄膜束缚水影响,利用自由流体分布构造毛细管压力曲线。
该技术中方法2中利用转换系数C将自由流体谱转换为伪毛管压力曲线,未考虑大小孔喉的转换关系不同这一特征。
方法1是对方法2的改进,利用不同幂函数构建伪毛管压力曲线,与压汞得到的毛细管压力曲线对某些岩样较以前吻合度有了一定的提高,但它还是通过对饱和水T2谱进行积分所得,其最大饱和度为100%,没有考虑到最大进汞饱和度问题,存在一定的局限性。
为了解决核磁测井定量评价孔隙结构的难题,在现有的研究基础上,即回到核磁原理,直接研究核磁饱和T2谱与压汞孔喉分布曲线的对应关系,及核磁自由流体T2谱与压汞孔喉分布≥0.0248μm部分的曲线(<0.0248μm的压汞孔喉分布曲线主要包含未连通孔喉的信息)的对应关系,并在分孔隙结构、分孔喉尺度的情况下构建反演核磁孔喉分布、自由流体孔喉分布公式。从而定量反演核磁孔喉分布和自由流体孔喉分布。提供定量评价孔隙结构的一种可靠方法和技术,打破现行方法的研究思想。推动核磁测技术的应用与发展。
发明内容
本发明针对现有技术的缺陷,提供了一种基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法,能有效的解决上述现有技术存在的问题。
为了实现以上发明目的,本发明采取的技术方案如下:
1.一种基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法,以下步骤:
第一步、挑选基础样品;
第二步、以压汞数据为基础,划分低渗砂岩储层孔隙结构类型,将低渗砂岩储层孔隙结构类型分为五种类型;
第三步、提取基本样品点对应的常规测井信息,常规测井中表现出AC、HRID值越高,而GR值越低的特征,以此特征为依据,构建AC-HRID-GR气泡图图版、雷达图图版识别孔隙结构;
第四步、挑选标准样品;
第五步、针对每个标准样品,提取自由流体T2谱与压汞孔喉分布的对应数据点,在分孔喉尺度<0.1μm、0.1-1μm、1-4μm的基础上,建立自由流体孔喉分布、饱和T2谱与压汞孔喉分布的对应关系,从而计算核磁孔喉分布、核磁自由流体孔喉分布;
第六步、依据第五部中所构建的对应关系,反演自由流体孔喉分布、核磁孔喉分布,并建立核磁反演孔喉分布直方图,将其与压汞孔喉分布直方图比较、验证;
第七步、在分孔隙结构的基础上,建立核磁饱和谱与核磁自由流体谱的对应关系,将核磁共振测井T2谱转换为核磁共振自由流体T2谱,进而计算地层自由流体孔喉分布;
第八步、沿井眼剖面定量反演低渗砂岩储层的孔隙结构类型。
作为优选,所述步骤一中基础样品的选取依据如下:
a、样品需用配套核磁实验和压汞实验的数据;
b、样品的常规测井信息能真实反映地层信息;
c、样品范围涵盖超低孔、特低孔、低孔和中孔,非渗、超低渗、特低渗和低渗孔隙结构品类型样。
作为优选,所述第二步中孔隙结构具体分为I1型、I2型、II1型、II2型和III型;
所述I1型排驱压力为Pd<0.2,最大进汞饱和度为SHgmax≥80;
所述I2型排驱压力为0.2≤Pd<0.4,最大进汞饱和度为70≤SHgmax<80;
所述II1型排驱压力为0.4≤Pd<0.7,最大进汞饱和度为60≤SHgmax<70;
所述II2型排驱压力为0.7≤Pd<1.6,最大进汞饱和度为50≤SHgmax<60;
所述III型排驱压力为Pd≥1.6,最大进汞饱和度为SHgmax<50。
作为优选,所述第四步中标准样品的选取依据如下:
a、从孔隙结构分类结果中筛选出标准样品;
b、标准样品能代表该种类型孔隙结构样品的特征。
作为优选,所述第五步中饱和T2谱与压汞孔喉分布的对应关系,如下式:
rc=a1×T2+b1 1≤rc<4μm
rc:压汞孔喉分布,μm;
T2:与rc对应的自由流体谱横向弛豫时间,ms。
作为优选,所述步骤七中建立核磁饱和谱与核磁自由流体谱的对应关系,如下式:
I1、I2、II1型孔隙结构:
y=a1×x+b1 T2饱和谱≥40ms
y=a2×x2+b2×x+c2 4.5≤T2饱和谱<40ms
y=a3×x2+b3×x+c3 T2饱和谱<4.5ms
II2、III型孔隙结构:
y=a1×ln(x)+b1 T2饱和谱≥2.25ms
y=a2×x2+b2×x+c2 0.7≤T2饱和谱<2.25ms
y=a3×x2+b3×x+c3 T2饱和谱<0.7ms
y:为核磁自由流体谱分量;
x:核磁饱和谱分量。
与现有技术相比本发明的优点在于:在前人研究基础上,采用新的研究思路与方法,即回到核磁原理,直接研究核磁饱和T2谱与压汞孔喉分布曲线的对应关系,及核磁自由流体T2谱与压汞孔喉分布≥0.0248μm部分的曲线(<0.0248μm的压汞孔喉分布曲线主要包含未连通孔喉的信息)的对应关系,并在分孔隙结构、分孔喉尺度的情况下构建反演核磁孔喉分布、自由流体孔喉分布公式。从而定量反演核磁孔喉分布和自由流体孔喉分布。提供定量评价孔隙结构的一种可靠方法和技术,打破现行方法的研究思想。自由流体T2与压汞孔喉分布对应关系更好,以此建立了不同孔隙结构类型二者之间不同孔喉尺度对应的关系式可在井筒剖面上通过识别孔隙结构类型,进而利用NMR定量反演核磁孔喉分布、自由流体孔喉分布,省去了构建伪毛管曲线环节,为低渗透砂岩储层有效性评价提供了直接依据,也是测井信息用于定量反演储层微观孔隙结构的有益探索,推动核磁测技术的应用与发展。
附图说明
图1为东营凹陷南斜坡沙河街组沙四段低渗致密砂岩压汞数据孔隙结构分类图;
图2为基本样品压汞特征参数排驱压力-物***会图;
图3为基本样品压汞特征参数SHgmax-物***会图;
图4为基本样品压汞特征参数平均孔喉半径-物***会图;
图5为常规测井“三性”雷达图识别图版;
图6为常规测井AC-GR-HRID气泡图识别图版;
图7为FX井AC-GR-HRID气泡图;
图8为I1型标准样品核磁T2谱特征;
图9为I1型标准样品压汞孔喉分布特征;
图10为I2型标准样品核磁T2谱特征;
图11为I2型标准样品压汞孔喉分布特征;
图12为II1型标准样品核磁T2谱特征;
图13为II1型标准样品压汞孔喉分布特征;
图14为II2型标准样品核磁T2谱特征;
图15为II2型标准样品压汞孔喉分布特征;
图16为III型标准样品核磁T2谱特征;
图17为III型标准样品压汞孔喉分布特征;
图18为I1型标准样品核磁自由流体谱与压汞孔喉分布对应图;
图19为I2型标准样品核磁自由流体谱与压汞孔喉分布对应图;
图20为II1型标准样品核磁自由流体谱与压汞孔喉分布对应图;
图21为II2型标准样品核磁自由流体谱与压汞孔喉分布对应图;
图22为III型标准样品核磁自由流体谱与压汞孔喉分布对应图;
图23为I1型标准样品自由流体孔喉分布反演公式图;
图24为I2型标准样品自由流体孔喉分布反演公式图;
图25为II1型标准样品自由流体孔喉分布反演公式图;
图26为II2型标准样品自由流体孔喉分布反演公式图;
图27为III型标准样品自由流体孔喉分布反演公式图;
图28为I1型标准样品核磁孔喉分布反演公式图;
图29为I2型标准样品核磁孔喉分布反演公式图;
图30为II1型标准样品核磁孔喉分布反演公式图;
图31为II2型标准样品核磁孔喉分布反演公式图;
图32为III型标准样品核磁孔喉分布反演公式图;
图33为I1型标准样品岩石孔喉分布反演图;
图34为I2型标准样品岩石孔喉分布反演图;
图35为II1型标准样品岩石孔喉分布反演图;
图36为II2型标准样品岩石孔喉分布反演图;
图37为III型标准样品岩石孔喉分布反演图;
图38为低渗透砂岩储层孔隙结构定量反演实例分析图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下参照附图并举实施例,对本发明做进一步详细说明。
一种基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法,包括以下步骤:
第一步、挑选基础样品,基本样品的选取依据:
a、样品需用配套核磁实验和压汞实验的数据,
b、样品必须是常规测井信息,
c、样品包括不同孔隙结构类型;
第二步、以压汞数据为基础,划分低渗砂岩储层孔隙结构类型:
由于低渗砂岩孔隙结构复杂,不同孔隙结构样品的压汞曲线特征差异较大。经对比分析,将低渗砂岩储层分为三大类五小类,如图1,分别为I1型(排驱压力低Pd<0.2,最大进汞饱和度高SHgmax≥80);I2型(排驱压力较低,0.2≤Pd<0.4,最大进汞饱和度较高,70≤SHgmax<80);II1型(排驱压力集中在0.4≤Pd<0.7,最大进汞饱和度分布在60≤SHgmax<70);II2型(排驱压力偏高,0.7≤Pd<1.6,最大进汞饱和度偏低,50≤SHgmax<60);III型(排驱压力高,Pd≥1.6,最大进汞饱和度低,SHgmax<50)。并结合物性资料,建立压汞特征参数-物性参数交会图,样品孔隙度、渗透率与排驱压力、最大进汞饱和度等参数存在良好的对应关系,如图2、3、4。当孔隙结构越好,即孔隙度、渗透率越大时,排驱压力越小,最大进汞饱和度越大。结合核磁T2谱分析发现,孔隙结构越好的样品束缚水饱和度越低。以此为依据,建立不同孔隙结构物性参数划分范围,见表1。
表1 ES4层段5种孔隙结构分类表
第三步、建立常规测井识别孔隙结构的方法:
在划分孔隙结构的情况下,提取基本样品点对应的常规测井信息,发现孔隙结构越好,常规测井中表现出AC、HRID值越高,而GR值越低的特征。以此特征为依据,构建AC-HRID-GR气泡图图版、雷达图图版识别孔隙结构,如图5、6,沿井筒剖面识别低渗砂岩孔隙结构类型。
以FX井为例,如图7,3241.25-3243.2m层段为II1孔隙结构,3243.2-3244.8m层段为在II2孔隙结构,3246.13-3248.7m层段为II1孔隙结构,3249.63-3251.8m层段为I2孔隙结构。评价准则见表2。
表2 常规测井识别孔隙结构准则表
第四步、挑选标准样品:
(1)基本样品的选取依据:
a、每种标准样品可代表该类孔隙结构类型;
b、孔隙结构特征明显;
5种标准样品特征描述:
a、I1型孔隙结构
连通性较好,排驱压力值低,一般<0.2MPa,最大进汞饱和度值高,往往>80%;核磁自由流体T2谱主要集中在10-300ms之间(图8);孔喉分布集中在0.73-1.495μm范围内(图9)。其CPOR=17.1%,CPERM=10.991mD,物性级别属于中孔低渗;此孔隙结构层段具有良好的渗透性。
b、I2型孔隙结构
连通性较差,Pd较低,在0.2-0.4MPa,SHgmax分布于70-80%;;核磁自由流体T2谱主要集中在10-300ms之间,但其含量低于I1型(图10);孔喉分布集中在0.1374-0.3721μm范围内(图11);CPOR=17.12%,CPERM=8.233mD,属于低孔特低渗结构,此孔隙结构层段的渗透性较好。
c、II1型孔隙结构
Pd分布于0.4-0.7MPa,SHgmax分布于60-70%;核磁自由流体T2谱主要集中在5-300ms之间(见图12);孔喉分布集中在0.1497-0.3801μm范围内,但含量低于I2型(图13)。CPOR=11.46%,CPERM=0.3627mD,属于特低孔超低渗结构。
d、II2型孔隙结构
Pd分布于0.7-1.6MPa,SHgmax分布于50-60%;核磁自由流体T2谱分布广,但自由孔隙含量低(见图14);孔喉分布在0.0936-0.1477μm范围内有一单峰,但其微孔含量高(图15)。CPOR=10%,CPERM=0.168mD,属属于特低孔超低渗结构,此孔隙结构层段的渗透性较差,为差渗透层。
e、III型孔隙结构
Pd值高,往往>1.6MPa,SHgmax值低,一般<50%;T2谱可见束缚孔隙度值大,而有效孔隙度值小,SBVI则>75%。(见图16),孔喉分布以微孔(<0.0248μm)分布为主(见图17);CPOR=5.87%,CPERM=0.0141mD,属于超低孔非渗结构,此孔隙结构层段的渗透性差,为典型非渗层。
第五步、分孔隙结构类型、孔喉尺度构建反演公式
针对每个标准样品,提取自由流体T2谱与压汞孔喉分布的对应数据点,如图18、19、20、21、22,在分孔喉尺度(<0.1μm、0.1-1μm、1-4μm)的基础上,如图23、24、25、26、27,自由流体孔喉分布反演公式见表3,同时建立饱和T2谱与压汞孔喉分布的对应关系,如图28、29、30、31、32,从而建立核磁孔喉分布反演公式见表4。
表3 5种孔隙结构自由流体谱与压汞孔喉分布的对应关系表
表4 5种孔隙结构核磁饱和谱与压汞孔喉分布的对应关系
第六步、模型验证
依据所构建的反演公式,反演自由流体孔喉分布、核磁孔喉分布,并建立核磁反演孔喉分布直方图,将其与压汞孔喉分布直方图比较、验证,如图33、34、35、36、37。
第七步、核磁自由流体谱反演
实际核磁测井不能直接提取核磁共振自由流体谱。故而,在分孔隙结构的基础上,建立核磁饱和谱与核磁自由流体谱的对应关系(如表5)。进而可利用核磁共振测井T2谱计算地层自由流体孔喉分布。
表5 5种孔隙结构核磁饱和谱与自由流体谱的对应关系
(注:y为自由流体谱分量;x为饱和谱分量)
第八步、沿井眼剖面定量反演低渗砂岩储层的自由流体孔喉分布。
以FX井为例,在常规测井信息沿井筒剖面识别孔隙结构的基础上,运用相应孔隙结构的反演公式将核磁共振测井T2谱沿井筒剖面方向定量反演为岩石孔隙空间分布,如图38。分析FX井核磁共振测井T2谱,发现3238.125-3241.2m、3244.8-3246.13m、3248.7-3249.625m层段T2谱位于最左边,明显孔隙结构差,为III型。3241.25-3243.2m(87号)、3246.13-3248.7m(89号)T2谱分布范围广,明显双峰特征为II1型,解释结论为差油层,且在核磁T2谱中高值分布含量高。3243.2-3244.8m(88号)T2谱单峰分布,为II2型,解释结论为油干层。3249.63-3251.8m(90号)T2谱分布范围广,并且分布在T2高值范围,为I2型,解释结论为油层。与图6AC-GR-HRID识别图结果一致,说明核磁T2谱识别孔隙结构和有效储层的适用性良好。进而并沿井筒剖面,通过核磁共振T2谱反演的核磁自由流体谱来计定量反演FX井核磁孔喉分布和自由流体孔喉分布。
需要注意的是,针对本技术方案,可以通过别的替代方案同样能完成发明目的,
如研究对象转换,可将低渗砂岩换为其他岩性,如砂岩、泥页岩、碳酸盐岩;
孔隙结构类型划分可利用薄片、电镜、CT扫面等其他资料建立;
井筒剖面孔隙结构分类可采用时频分析方法或直接用核磁测井进行分类。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的实施方法,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (6)
1.一种基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法,其特征在于包括以下步骤:
第一步、挑选基础样品;
第二步、以压汞数据为基础,划分低渗砂岩储层孔隙结构类型,将低渗砂岩储层孔隙结构类型分为五种类型;
第三步、提取基本样品点对应的常规测井信息,常规测井中表现出AC、HRID值越高,而GR值越低的特征,以此特征为依据,构建AC-HRID-GR气泡图图版、雷达图图版识别孔隙结构;
第四步、挑选标准样品;
第五步、针对每个标准样品,提取自由流体T2谱与压汞孔喉分布的对应数据点,在分孔喉尺度<0.1μm、0.1-1μm、1-4μm的基础上,建立自由流体孔喉分布、饱和T2谱与压汞孔喉分布的对应关系,从而计算核磁孔喉分布、核磁自由流体孔喉分布;
第六步、依据第五部中所构建的对应关系,反演自由流体孔喉分布、核磁孔喉分布,并建立核磁反演孔喉分布直方图,将其与压汞孔喉分布直方图比较、验证;
第七步、在分孔隙结构的基础上,建立核磁饱和谱与核磁自由流体谱的对应关系,将核磁共振测井T2谱转换为核磁共振自由流体T2谱,进而计算地层自由流体孔喉分布;
第八步、沿井眼剖面定量反演低渗砂岩储层的自由流体孔喉分布。
2.根据权利要求1的一种基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法,其特征在于:所述步骤一中基础样品的选取依据如下:
a、样品需用配套核磁实验和压汞实验的数据;
b、样品的常规测井信息能真实反映地层信息;
c、样品范围涵盖超低孔、特低孔、低孔和中孔,非渗、超低渗、特低渗和低渗孔隙结构品类型样。
3.根据权利要求1或2的一种基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法,其特征在于:所述第二步中孔隙结构具体分为I1型、I2型、II1型、II2型和III型;
所述I1型排驱压力为Pd<0.2,最大进汞饱和度为SHgmax≥80;
所述I2型排驱压力为0.2≤Pd<0.4,最大进汞饱和度为70≤SHgmax<80;
所述II1型排驱压力为0.4≤Pd<0.7,最大进汞饱和度为60≤SHgmax<70;
所述II2型排驱压力为0.7≤Pd<1.6,最大进汞饱和度为50≤SHgmax<60;
所述III型排驱压力为Pd≥1.6,最大进汞饱和度为SHgmax<50。
4.根据权利要求3的一种基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法,其特征在于:所述第四步中标准样品的选取依据如下:
a、从孔隙结构分类结果中筛选出标准样品;
b、标准样品能代表该种类型孔隙结构样品的特征。
5.根据权利要求4的一种基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法,其特征在于:所述第五步中饱和T2谱与压汞孔喉分布的对应关系,如下式:
rc=a1×T2+b1 1≤rc<4μm
rc:压汞孔喉分布,μm;
T2:与rc对应的自由流体谱横向弛豫时间,ms。
6.根据权利要求5的一种基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法,其特征在于:所述步骤七中建立核磁饱和谱与核磁自由流体谱的对应关系,如下式:
I1、I2、II1型孔隙结构:
y=a1×x+b1 T2饱和谱≥40ms
y=a2×x2+b2×x+c2 4.5≤T2饱和谱<40ms
y=a3×x2+b3×x+c3 T2饱和谱<4.5ms
II2、III型孔隙结构:
y=a1×ln(x)+b1 T2饱和谱≥2.25ms
y=a2×x2+b2×x+c2 0.7≤T2饱和谱<2.25ms
y=a3×x2+b3×x+c3 T2饱和谱<0.7ms
y:为核磁自由流体谱分量;
x:核磁饱和谱分量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510251449.XA CN104819923B (zh) | 2015-05-17 | 2015-05-17 | 基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510251449.XA CN104819923B (zh) | 2015-05-17 | 2015-05-17 | 基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104819923A true CN104819923A (zh) | 2015-08-05 |
CN104819923B CN104819923B (zh) | 2017-10-10 |
Family
ID=53730280
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510251449.XA Expired - Fee Related CN104819923B (zh) | 2015-05-17 | 2015-05-17 | 基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104819923B (zh) |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105241798A (zh) * | 2015-09-23 | 2016-01-13 | 中国海洋石油总公司 | 一种复杂碳酸盐岩储层渗透率的定量表征方法 |
CN106053315A (zh) * | 2016-06-06 | 2016-10-26 | 中国石油大学(华东) | 一种生物碎屑灰岩储层孔隙结构分类方法 |
CN106351652A (zh) * | 2016-11-25 | 2017-01-25 | 中国地质大学(北京) | 一种含烃储集层核磁共振测井t2谱形态校正方法 |
CN106770415A (zh) * | 2016-12-13 | 2017-05-31 | 中国石油天然气股份有限公司 | 一种高粘度低流度致密油在不同温压条件下的微观可流动性的评价方法 |
CN106908370A (zh) * | 2017-03-22 | 2017-06-30 | 中国石油天然气股份有限公司 | 砂岩孔渗参数的确定方法和装置 |
CN107817199A (zh) * | 2016-09-14 | 2018-03-20 | 中国石油化工股份有限公司 | 一种致密砂岩多尺度孔隙模型的构建方法和应用 |
CN107816349A (zh) * | 2016-09-12 | 2018-03-20 | 中国石油化工股份有限公司 | 一种分析致密砂岩孔隙结构状态的方法 |
CN107894386A (zh) * | 2017-11-14 | 2018-04-10 | 西安石油大学 | 超临界二氧化碳注入对低渗透砂岩油藏孔喉结构影响的定量评价方法 |
CN108169099A (zh) * | 2018-01-17 | 2018-06-15 | 西南石油大学 | 一种基于核磁共振的页岩气储层孔隙结构定量计算方法 |
CN108181219A (zh) * | 2017-11-15 | 2018-06-19 | 中国石油天然气股份有限公司 | 一种基于碳酸盐岩孔隙结构预测渗透率的方法及其装置 |
CN109883889A (zh) * | 2019-01-24 | 2019-06-14 | 西南石油大学 | 模拟co2在致密基质-裂缝扩散的实验装置及前缘预测方法 |
CN110160934A (zh) * | 2019-06-06 | 2019-08-23 | 西安石油大学 | 基于核磁共振技术评价致密储层孔隙结构分类标准的方法 |
CN110987761A (zh) * | 2019-12-13 | 2020-04-10 | 重庆科技学院 | 一种定量计算岩心孔隙连通性的方法 |
CN111379560A (zh) * | 2018-12-29 | 2020-07-07 | 中国石油化工股份有限公司 | 一种构建储层油水性质判识图版的方法 |
CN112505085A (zh) * | 2021-02-05 | 2021-03-16 | 西南石油大学 | 基于核磁共振的孔隙度有效应力系数测定方法 |
CN114112841A (zh) * | 2021-11-19 | 2022-03-01 | 中国石油大学(华东) | 一种致密砂岩储层束缚水饱和度的计算方法 |
CN114233284A (zh) * | 2021-09-28 | 2022-03-25 | 重庆科技学院 | 基于核磁视自由水孔隙度反演的致密储层高阻水层识别方法 |
WO2023184214A1 (en) * | 2022-03-30 | 2023-10-05 | Saudi Arabian Oil Company | Parameterization of nuclear magnetic resonance transverse relaxation time distribution |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102141637A (zh) * | 2010-01-28 | 2011-08-03 | 中国石油天然气股份有限公司 | 一种利用核磁共振测井资料连续定量评价储集层孔隙结构的方法 |
CN102998322A (zh) * | 2011-09-14 | 2013-03-27 | 中国石油天然气股份有限公司 | 恒定梯度场核磁共振岩样分析方法及仪器 |
CN103256046A (zh) * | 2013-04-28 | 2013-08-21 | 北京大学 | 非常规油气藏水平井全缝长压裂参数模拟的方法及装置 |
US20140055134A1 (en) * | 2012-08-23 | 2014-02-27 | Schlumberger Technology Corporation | Magnetic resonance examination of porous samples |
CN103778328A (zh) * | 2014-01-10 | 2014-05-07 | 中国石油大学(华东) | 基于标准偏差分析的储层敏感孔喉提取方法 |
CN103884633A (zh) * | 2014-03-05 | 2014-06-25 | 中国石油天然气股份有限公司 | 一种确定岩石渗透率的方法及装置 |
CN104278989A (zh) * | 2013-07-02 | 2015-01-14 | 中国石油天然气股份有限公司 | 一种获取低孔低渗储层饱和度指数的方法 |
-
2015
- 2015-05-17 CN CN201510251449.XA patent/CN104819923B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102141637A (zh) * | 2010-01-28 | 2011-08-03 | 中国石油天然气股份有限公司 | 一种利用核磁共振测井资料连续定量评价储集层孔隙结构的方法 |
CN102998322A (zh) * | 2011-09-14 | 2013-03-27 | 中国石油天然气股份有限公司 | 恒定梯度场核磁共振岩样分析方法及仪器 |
US20140055134A1 (en) * | 2012-08-23 | 2014-02-27 | Schlumberger Technology Corporation | Magnetic resonance examination of porous samples |
CN103256046A (zh) * | 2013-04-28 | 2013-08-21 | 北京大学 | 非常规油气藏水平井全缝长压裂参数模拟的方法及装置 |
CN104278989A (zh) * | 2013-07-02 | 2015-01-14 | 中国石油天然气股份有限公司 | 一种获取低孔低渗储层饱和度指数的方法 |
CN103778328A (zh) * | 2014-01-10 | 2014-05-07 | 中国石油大学(华东) | 基于标准偏差分析的储层敏感孔喉提取方法 |
CN103884633A (zh) * | 2014-03-05 | 2014-06-25 | 中国石油天然气股份有限公司 | 一种确定岩石渗透率的方法及装置 |
Cited By (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105241798A (zh) * | 2015-09-23 | 2016-01-13 | 中国海洋石油总公司 | 一种复杂碳酸盐岩储层渗透率的定量表征方法 |
CN106053315A (zh) * | 2016-06-06 | 2016-10-26 | 中国石油大学(华东) | 一种生物碎屑灰岩储层孔隙结构分类方法 |
CN106053315B (zh) * | 2016-06-06 | 2019-02-05 | 中国石油大学(华东) | 一种生物碎屑灰岩储层孔隙结构分类方法 |
CN107816349A (zh) * | 2016-09-12 | 2018-03-20 | 中国石油化工股份有限公司 | 一种分析致密砂岩孔隙结构状态的方法 |
CN107817199A (zh) * | 2016-09-14 | 2018-03-20 | 中国石油化工股份有限公司 | 一种致密砂岩多尺度孔隙模型的构建方法和应用 |
CN107817199B (zh) * | 2016-09-14 | 2019-09-03 | 中国石油化工股份有限公司 | 一种致密砂岩多尺度孔隙模型的构建方法和应用 |
CN106351652B (zh) * | 2016-11-25 | 2019-03-29 | 中国地质大学(北京) | 一种含烃储集层核磁共振测井t2谱形态校正方法 |
CN106351652A (zh) * | 2016-11-25 | 2017-01-25 | 中国地质大学(北京) | 一种含烃储集层核磁共振测井t2谱形态校正方法 |
CN106770415A (zh) * | 2016-12-13 | 2017-05-31 | 中国石油天然气股份有限公司 | 一种高粘度低流度致密油在不同温压条件下的微观可流动性的评价方法 |
CN106770415B (zh) * | 2016-12-13 | 2019-03-15 | 中国石油天然气股份有限公司 | 一种高粘度低流度致密油在不同温压条件下的微观可流动性的评价方法 |
CN106908370A (zh) * | 2017-03-22 | 2017-06-30 | 中国石油天然气股份有限公司 | 砂岩孔渗参数的确定方法和装置 |
CN106908370B (zh) * | 2017-03-22 | 2019-08-06 | 中国石油天然气股份有限公司 | 砂岩孔渗参数的确定方法和装置 |
CN107894386A (zh) * | 2017-11-14 | 2018-04-10 | 西安石油大学 | 超临界二氧化碳注入对低渗透砂岩油藏孔喉结构影响的定量评价方法 |
CN108181219A (zh) * | 2017-11-15 | 2018-06-19 | 中国石油天然气股份有限公司 | 一种基于碳酸盐岩孔隙结构预测渗透率的方法及其装置 |
CN108181219B (zh) * | 2017-11-15 | 2020-05-08 | 中国石油天然气股份有限公司 | 一种基于碳酸盐岩孔隙结构预测渗透率的方法及其装置 |
CN108169099A (zh) * | 2018-01-17 | 2018-06-15 | 西南石油大学 | 一种基于核磁共振的页岩气储层孔隙结构定量计算方法 |
CN111379560A (zh) * | 2018-12-29 | 2020-07-07 | 中国石油化工股份有限公司 | 一种构建储层油水性质判识图版的方法 |
CN111379560B (zh) * | 2018-12-29 | 2023-04-07 | 中国石油化工股份有限公司 | 一种构建储层油水性质判识图版的方法 |
CN109883889A (zh) * | 2019-01-24 | 2019-06-14 | 西南石油大学 | 模拟co2在致密基质-裂缝扩散的实验装置及前缘预测方法 |
CN109883889B (zh) * | 2019-01-24 | 2021-09-28 | 西南石油大学 | 模拟co2在致密基质-裂缝扩散的实验装置及前缘预测方法 |
CN110160934A (zh) * | 2019-06-06 | 2019-08-23 | 西安石油大学 | 基于核磁共振技术评价致密储层孔隙结构分类标准的方法 |
CN110987761B (zh) * | 2019-12-13 | 2022-03-18 | 重庆科技学院 | 一种定量计算岩心孔隙连通性的方法 |
CN110987761A (zh) * | 2019-12-13 | 2020-04-10 | 重庆科技学院 | 一种定量计算岩心孔隙连通性的方法 |
CN112505085B (zh) * | 2021-02-05 | 2021-04-09 | 西南石油大学 | 基于核磁共振的孔隙度有效应力系数测定方法 |
CN112505085A (zh) * | 2021-02-05 | 2021-03-16 | 西南石油大学 | 基于核磁共振的孔隙度有效应力系数测定方法 |
CN114233284A (zh) * | 2021-09-28 | 2022-03-25 | 重庆科技学院 | 基于核磁视自由水孔隙度反演的致密储层高阻水层识别方法 |
CN114233284B (zh) * | 2021-09-28 | 2022-07-15 | 重庆科技学院 | 基于核磁视自由水孔隙度反演的致密储层高阻水层识别方法 |
CN114112841A (zh) * | 2021-11-19 | 2022-03-01 | 中国石油大学(华东) | 一种致密砂岩储层束缚水饱和度的计算方法 |
WO2023184214A1 (en) * | 2022-03-30 | 2023-10-05 | Saudi Arabian Oil Company | Parameterization of nuclear magnetic resonance transverse relaxation time distribution |
Also Published As
Publication number | Publication date |
---|---|
CN104819923B (zh) | 2017-10-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104819923A (zh) | 基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法 | |
Zheng et al. | Nuclear magnetic resonance T2 cutoffs of coals: A novel method by multifractal analysis theory | |
CN102175832B (zh) | 一种确定典型储层最佳饱和度计算模型的方法 | |
Xiao et al. | Combining rate-controlled porosimetry and NMR to probe full-range pore throat structures and their evolution features in tight sands: A case study in the Songliao Basin, China | |
CN103267721B (zh) | 一种致密砂岩储层孔隙含水特征与赋存状态的评价方法 | |
CN105927218B (zh) | 一种陆相泥页岩储层含气量预测方法及装置 | |
CN102141637B (zh) | 一种利用核磁共振测井资料连续定量评价储集层孔隙结构的方法 | |
CN103278436B (zh) | 特低渗透双重介质砂岩油藏微观孔隙结构的定量表征方法 | |
CN103884633B (zh) | 一种确定岩石渗透率的方法及装置 | |
CN102200008B (zh) | 一种基于电成像测井的储层有效性识别方法 | |
Jian-Ping et al. | A QUANTITATIVE EVALUATION METHOD OF LOW PERMEABLE SANDSTONE PORE STRUCTURE BASED ON NUCLEAR MAGNETIC RESONANCE (NMR) LOGGING: A CASE STUDY OF Es4 FORMATION IN THE SOUTH SLOPE OF DONGYING SAG | |
CN102434152A (zh) | 一种储层含油饱和度的计算方法 | |
CN104278989B (zh) | 一种获取低孔低渗储层饱和度指数的方法 | |
CN109030311A (zh) | 基于核磁共振t2谱敏感参数的孔隙结构分类与识别方法 | |
CN104806232B (zh) | 一种确定碳酸盐岩储层孔隙度下限的方法 | |
CN104634718A (zh) | 应用核磁共振表征致密砂岩孔径分布的标定方法 | |
CN105525917B (zh) | 识别储层的方法以及识别储层中流体类型的方法 | |
CN105781539A (zh) | 一种致密油气储层饱和度测井计算方法 | |
CN105447762B (zh) | 一种流体替换的低渗透油藏水淹信息的计算方法 | |
Li et al. | The pore-fracture system properties of coalbed methane reservoirs in the Panguan Syncline, Guizhou, China | |
CN108444881B (zh) | 一种适用于陆相泥页岩微纳米尺度储集空间的表征方法 | |
Hu et al. | A method to determine nuclear magnetic resonance T2 cutoff value of tight sandstone reservoir based on multifractal analysis | |
CN105445441A (zh) | 储层含气饱和度确定方法及装置 | |
CN106443800A (zh) | 一种中基性火山岩的核磁孔隙度校正方法 | |
CN103412335B (zh) | 一种利用地震物相体预测储层的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171010 Termination date: 20190517 |
|
CF01 | Termination of patent right due to non-payment of annual fee |