首页 / 专利库 / 地球科学 / 岩石 / 沉积岩 / 砂岩 / 一种通过测井数据获取页岩储层含气量的方法

一种通过测井数据获取页岩储层含气量的方法

阅读:888发布:2021-12-31

专利汇可以提供一种通过测井数据获取页岩储层含气量的方法专利检索,专利查询,专利分析的服务。并且本 发明 涉及一种通过 测井 数据获取 页岩 储层含气量的方法,其包括以下步骤:步骤一:在研究区中选择页岩气井进行取心,然后对同一个页岩样品分别测量孔隙度有机 碳 含量TOCcore、含 水 饱和度 Sw-core、胶结指数m、饱和度指数n、兰格缪尔体积VL、兰格缪尔压 力 PL、 地层 骨架 密度 ρm和地层 流体 密度ρfl;步骤二:通过测井数据计算获取页岩储层含气量的关键参数;步骤三:利用上述关键参数分别计算 吸附 气含量和游离气含量,并用吸附气和游离气之和来表征页岩储层总含气量。本发明通过测井数据计算页岩储层含气量,减少了钻井取心的成本,缩短了储层含气性评价的周期。同时,本发明分别从吸附气、游离气两个相态计算储层的页岩储层中 天然气 的总量,符合页岩储层含气性特点,计算 精度 更高。,下面是一种通过测井数据获取页岩储层含气量的方法专利的具体信息内容。

1.一种通过测井数据获取页岩储层含气量的方法,其包括以下步骤:
步骤一:在研究区中选择页岩气井进行取心,然后对同一个页岩样品分别测量孔隙度有机含量TOCcore、含饱和度Sw-core、胶结指数m、饱和度指数n、兰格缪尔体积VL、兰格缪尔压PL、地层骨架密度ρm和地层流体密度ρfl;
步骤二:通过测井数据计算获取页岩储层含气量的关键参数:
1)利用多元线性回归法获取测井有机碳含量TOClog:
①相关性检测:使用SPSS软件将步骤一中岩心实测的页岩样品TOCcore数据与页岩样品采样点相同深度的各测井曲线值进行双变量相关性分析,得到TOCcore与各测井曲线的相关性表,并从相关性表选择在0.01水平上显著相关的测井曲线,包括自然伽GR、体积密度DEN和中子CNL;
②多元线性回归:使用步骤一中岩心实测的页岩样品TOCcore与显著相关的测井曲线使用SPSS软件进行多元线性回归,得到复相关系数R、统计量F和多元线性回归方程(1):
TOCcore=a×GR+b×DEN+c×CNL+d (1)
式中,TOCcore为步骤一中岩心实测的页岩样品有机碳含量,wt%;a、b、c、d为多元回归系数,通过SPSS软件进行线性回归后得到;
③F统计检验:给定显著性水平0.01、多元回归样品容量M和多元线性回归方程元数N,通过查F分布临界值表得到F0.01(N,M-N-1),如果F>F0.01(N,M-N-1),多元线性关系显著,则多元线性回归方程(1)成立;如果F≤F0.01(N,M-N-1),多元线性关系不显著,则略去相关性分析中Pearson相关系数绝对值较小的测井曲线,重新进行多元线性回归,直到F>F0.01(n,M-n-1)成立为止,则通过测井计算的有机碳含量的公式为:
TOClog=a×GR+b×DEN+c×CNL+d (2)
2)计算孔隙度:
首先建立岩石体积模型,并通过体积密度得到体积方程:
式中, 为测井计算的孔隙度,%;ρm为地层骨架密度,g/cm3;ρb为地层体积密度,g/
3
cm;VTOC为有机碳含量体积百分比,vol%;wTOC为测井计算的有机碳含量TOClog,wt%;ρfl
3
为地层流体密度,g/cm;k为校正因子;
根据以上体积方程得到孔隙度的计算公式为:
将公式(5)计算得到的孔隙度 与步骤一中岩心实测的页岩样品孔隙度 进行线性回归:
式中,e、f为线性回归系数,通过SPSS软件进行线性回归后得到; 为校正后的测井计算孔隙度,%;
3)计算含气饱和度:
页岩含水饱和度根据阿尔奇公式计算:
式中,Sw为含水饱和度,%;Rw为地层水电阻率,ohmm; 为测井计算的孔隙度,%;Rt为地层电阻率,ohmm;m为胶结指数;n为饱和度指数;
将公式(8)计算得到的含水饱和度Sw与步骤一中岩心实测的页岩样品含水饱和度Sw-core进行线性回归:
Sw-core=g×Sw+h (9)
Sw-log=g×Sw+h (10)
式中,Sw-log为经过岩心校正的测井计算的含水饱和度,%;g、h为线性回归系数,通过SPSS软件进行线性回归后得到;
页岩含气饱和度为:
Sg=1-Sw-log (11)
步骤三:计算页岩含气量:
由于在计算页岩储层原地资源量时使用吸附气和游离气之和来表征总含气量,故需要对吸附气含量和游离气含量分别进行计算:
1)吸附气含量:
根据地层实际温度、压力和有机碳含量对页岩样品进行校正,具体包括利用公式(12)进行兰格缪尔压力校正,利用公式(13)进行兰格缪尔体积校正,以及利用公式(14)进行有机碳含量校正:
Vlt=33.3357Vl×1.0062Ti-T (12)
Plt=145.1379Pl×1.0116T-Ti (13)
3
式中,Vlt为储层温度校正的兰格缪尔体积,m /t;Plt为储层温度校正的兰格缪尔压力,
3
MPa;T为储层温度,℃;Ti为等温吸附实验温度,℃;Vl为页岩样品兰格缪尔体积,m /t;Pl
3
为页岩样品兰格缪尔压力,MPa;Vlc为经过储层温度和有机碳含量校正的兰格缪尔体积,m /t;
根据兰格缪尔方程计算页岩储层吸附气量:
3
式中,Ga为吸附气含量,m /t;P为储层压力,MPa;H为储层深度,m;
2)游离气含量:
游离气量是页岩孔隙空间中的含气量,通过体积模型计算:
3 3 3
式中,Gf为游离气量,m /t;Bg为气体压缩系数,m /m; 为测井计算的孔隙度,%;Sg为页岩含气饱和度,%;ρb为地层体积密度,g/cm3。
由于吸附气的存在会影响游离气的容纳空间,故需减去吸附气所占的体积空间,因此游离气实际含量为:
3
式中, 为天然气的视分子重量,g/mol;ρs为吸附态甲烷密度,g/cm ;
最终得到页岩储层总含气量为:
Gt=Ga+Gf (19)
3 3 3
式中,Gt为页岩储层总含气量,m /t;Ga为页岩吸附气量,m /t;Gf为页岩游离气量,m /t。
2.如权利要求1所述的一种通过测井数据获取页岩储层含气量的方法,其特征在于,在公式(8)中,地层电阻率Rt在页岩储层中常取与页岩相邻或相近的灰岩或砂岩地层水电阻率。
3.如权利要求1所述的一种通过测井数据获取页岩储层含气量的方法,其特征在于,在公式(15)中,储层压力P根据地层压力梯度计算,正常压力梯度情况下P=
9.8×H/1000。
4.如权利要求1所述的一种通过测井数据获取页岩储层含气量的方法,其特征在于,在公式(18)中,天然气的视分子重量 根据解吸气组分计算,甲烷为16;吸附态甲烷密度ρs=0.375-0.4233,页岩储层根据经验取0.4。

说明书全文

一种通过测井数据获取页岩储层含气量的方法

技术领域

[0001] 本发明涉及一种获取页岩储层含气量的方法,具体涉及一种通过测井数据获取页岩储层含气量的方法。

背景技术

[0002] 含气量是页岩储层评价的关键参数,对页岩含气性评价、储量预测具有重要意义。含气量确定包括现场解吸法、等温吸附法和测井解释法。现场解吸法反映了样品的实际含气量,但受取心方式和提钻时间影响,损失气估算误差较大。等温吸附法获得的是页岩的最大吸附气量,未考虑游离气量部分,与实际含气量有较大差别。测井解释法是通过测井数据结合岩心实验建立含气量的测井解释模型,可以预测页岩的含气量,减少实验成本。Decker研究发现,页岩实测含气量与TOC(总有机含量)存在很好的正线性相关关系,TOC与体积密度存在很好的负线性相关关系,因此可以通过建立体积密度与含气量的计算模型来预测页岩含气量。Cluff等以等温吸附和体积模型为基础,利用测井解释参数计算了Delaware盆地中Barnett页岩和Woodford页岩原地资源量。Utley等使用神经网络计算关键参数基于等温吸附和体积模型计算了Fayetteville页岩含气量。
[0003] 此外,申请号为201410249240.5,发明名称为“一种利用测井曲线获得页岩气含气量的方法”的中国发明专利申请提出通过解析实验所得页岩气含气量与页岩测井数据线性回归,得到页岩气含气量与声波时差、电阻率的一元线性回归方程及二元线性回归方程计算页岩含气量。郭旭升等通过现场实测总含气量与实验分析TOC的线性关系建立一元线性模型计算涪陵气田页岩储层含气量。但是,通过建立实测含气量与测井参数之间的线性模型预测含气量在层气井中应用效果较好,而线性回归法对含气量较低的井误差较高,页岩含气量总体远小于煤层含气量。另外页岩气井在钻井取心过程中,受取心方式影响气体损失严重,损失气量占总含气量的40-70%,损失气估算结果误差较大,因此使用测井数据与实测含气量之间的线性或非线性模型预测页岩含气量可能带来很大的误差。

发明内容

[0004] 针对上述问题,本发明的目的是提供一种通过测井数据获取页岩储层含气量的方法。
[0005] 为实现上述目的,本发明采取以下技术方案:一种通过测井数据获取页岩储层含气量的方法,其包括以下步骤:
[0006] 步骤一:在研究区中选择页岩气井进行取心,然后对同一个页岩样品分别测量孔隙度 有机碳含量TOCcore、含饱和度Sw-core、胶结指数m、饱和度指数n、兰格缪尔体积VL、兰格缪尔压PL、地层骨架密度ρm和地层流体密度ρfl;
[0007] 步骤二:通过测井数据计算获取页岩储层含气量的关键参数:
[0008] 1)利用多元线性回归法获取测井有机碳含量TOClog:
[0009] ①相关性检测:使用SPSS软件将步骤一中岩心实测的页岩样品TOCcore数据与页岩样品采样点相同深度的各测井曲线值进行双变量相关性分析,得到TOCcore与各测井曲线的相关性表,并从相关性表选择在0.01水平上显著相关的测井曲线,包括自然伽GR、体积密度DEN和中子CNL;
[0010] ②多元线性回归:使用步骤一中岩心实测的页岩样品TOCcore与显著相关的测井曲线使用SPSS软件进行多元线性回归,得到复相关系数R、统计量F和多元线性回归方程(1):
[0011] TOCcore=a×GR+b×DEN+c×CNL+d (1)
[0012] 式中,TOCcore为步骤一中岩心实测的页岩样品有机碳含量,wt%;a、b、c、d为多元回归系数,通过SPSS软件进行线性回归后得到;
[0013] ③F统计检验:给定显著性水平0.01、多元回归样品容量M和多元线性回归方程元数N,通过查F分布临界值表得到F0.01(N,M-N-1),如果F>F0.01(N,M-N-1),多元线性关系显著,则多元线性回归方程(1)成立;如果F≤F0.01(N,M-N-1),多元线性关系不显著,则略去相关性分析中Pearson相关系数绝对值较小的测井曲线,重新进行多元线性回归,直到F>F0.01(n,M-n-1)成立为止,则通过测井计算的有机碳含量的公式为:
[0014] TOClog=a×GR+b×DEN+c×CNL+d (2)
[0015] 2)计算孔隙度:
[0016] 首先建立岩石体积模型,并通过体积密度得到体积方程:
[0017]
[0018]
[0019] 式中,为测井计算的孔隙度,%;ρm为地层骨架密度,g/cm3;ρb为地层体积密3
度,g/cm;VTOC为有机碳含量体积百分比,vol%;wTOC为测井计算的有机碳含量TOClog,wt%;
3
ρfl为地层流体密度,g/cm ;k为校正因子;
[0020] 根据以上体积方程得到孔隙度的计算公式为:
[0021]
[0022] 将公式(5)计算得到的孔隙度 与步骤一中岩心实测的页岩样品孔隙度 进行线性回归:
[0023]
[0024]
[0025] 式中,e、f为线性回归系数,通过SPSS软件进行线性回归后得到; 为校正后的测井计算孔隙度,%;
[0026] 3)计算含气饱和度:
[0027] 页岩含水饱和度根据阿尔奇公式计算:
[0028]
[0029] 式中,Sw为含水饱和度,%;Rw为地层水电阻率,ohmm; 为测井计算的孔隙度,%;Rt为地层电阻率,ohmm;m为胶结指数;n为饱和度指数;
[0030] 将公式(8)计算得到的含水饱和度Sw与步骤一中岩心实测的页岩样品含水饱和度Sw-core进行线性回归:
[0031] Sw-core=g×Sw+h (9)
[0032] Sw-log=g×Sw+h (10)
[0033] 式中,Sw-log为经过岩心校正的测井计算的含水饱和度,%;g、h为线性回归系数,通过SPSS软件进行线性回归后得到;
[0034] 页岩含气饱和度为:
[0035] Sg=1-Sw-log (11)
[0036] 步骤三:计算页岩含气量:
[0037] 由于在计算页岩储层原地资源量时使用吸附气和游离气之和来表征总含气量,故需要对吸附气含量和游离气含量分别进行计算:
[0038] 1)吸附气含量:
[0039] 根据地层实际温度、压力和有机碳含量对页岩样品进行校正,具体包括利用公式(12)进行兰格缪尔压力校正,利用公式(13)进行兰格缪尔体积校正,以及利用公式(14)进行有机碳含量校正:Ti-T
[0040] Vlt=33.3357Vl×1.0062 (12)T-Ti
[0041] Plt=145.1379Pl×1.0116 (13)
[0042]3
[0043] 式中,Vlt为储层温度校正的兰格缪尔体积,m /t;Plt为储层温度校正的兰格缪尔压力,MPa;T为储层温度,℃;Ti为等温吸附实验温度,℃;Vl为页岩样品兰格缪尔体积,3
m/t;Pl为页岩样品兰格缪尔压力,MPa;Vlc为经过储层温度和有机碳含量校正的兰格缪尔
3
体积,m/t;
[0044] 根据兰格缪尔方程计算页岩储层吸附气量:
[0045]
[0046] 式中,Ga为吸附气含量,m3/t;P为储层压力,MPa;H为储层深度,m;
[0047] 2)游离气含量:
[0048] 游离气量是页岩孔隙空间中的含气量,通过体积模型计算:
[0049]
[0050] 式中,Gf为游离气量,m3/t;Bg为气体压缩系数,m3/m3;为测井计算的孔隙度,%;3
Sg为页岩含气饱和度,%;ρb为地层体积密度,g/cm 。
[0051] 由于吸附气的存在会影响游离气的容纳空间,故需减去吸附气所占的体积空间,因此游离气实际含量为:
[0052]
[0053]
[0054] 式中, 为天然气的视分子重量,g/mol;ρs为吸附态甲烷密度,g/cm3;
[0055] 最终得到页岩储层总含气量为:
[0056] Gt=Ga+Gf (19)
[0057] 式中,Gt为页岩储层总含气量,m3/t;Ga为页岩吸附气量,m3/t;Gf为页岩游离气量,3
m/t。
[0058] 在公式(8)中,地层电阻率Rt在页岩储层中常取与页岩相邻或相近的灰岩或砂岩地层水电阻率。
[0059] 在公式(15)中,储层压力P根据地层压力梯度计算,正常压力梯度情况下P=9.8×H/1000。
[0060] 在公式(18)中,天然气的视分子重量 根据解吸气组分计算,甲烷为16;吸附态甲烷密度ρs=0.375-0.4233,页岩储层根据经验取0.4。
[0061] 本发明由于采取以上技术方案,其具有以下优点:1、本发明通过测井数据计算页岩储层含气量,减少了钻井取心的成本,缩短了储层含气性评价的周期。2、本发明分别从吸附气、游离气两个相态计算储层的页岩储层中天然气的总量,符合页岩储层含气性特点,计算精度更高。3、本发明基于页岩储层的岩电特征和吸附特征,具有普遍使用性,区域推广应用效果更好。

具体实施方式

[0062] 本发明提供的通过测井数据获取页岩储层含气量的方法包括以下步骤:
[0063] 步骤一:在研究区中选择页岩气井进行取心,然后对同一个页岩样品分别测量孔隙度 有机碳含量TOCcore、含水饱和度Sw-core、胶结指数m、饱和度指数n、兰格缪尔体积VL、兰格缪尔压力PL、地层骨架密度ρm和地层流体密度ρfl。取心的井位和层位以能够反映研究区页岩基本特征为原则。
[0064] 步骤二:通过测井数据计算获取页岩储层含气量的关键参数:
[0065] 1)利用多元线性回归法获取测井有机碳含量TOClog:
[0066] ①相关性检测:使用SPSS(Statistical Product and Service Solutions)软件将步骤一中岩心实测的页岩样品TOCcore数据与页岩样品采样点相同深度的各测井曲线值进行双变量相关性分析,得到TOCcore与各测井曲线的相关性表,并从相关性表选择在0.01水平(双侧)上显著相关的测井曲线。多数情况下,与TOCcore显著相关的测井曲线为自然伽马GR、体积密度DEN和中子CNL。
[0067] ②多元线性回归:使用步骤一中岩心实测的页岩样品TOCcore与显著相关的测井曲线(GR、DEN、CNL)使用SPSS软件进行多元线性回归,得到复相关系数R(R只是作为一个相关性参考)、统计量F和多元线性回归方程(1):
[0068] TOCcore=a×GR+b×DEN+c×CNL+d (1)
[0069] 式中,TOCcore为步骤一中岩心实测的页岩样品有机碳含量,wt%;a、b、c、d为多元回归系数,可以通过SPSS软件进行线性回归后得到。
[0070] ③F统计检验:给定显著性水平0.01、多元回归样品容量M和多元线性回归方程元数N,通过查F分布临界值表得到F0.01(N,M-N-1),如果F>F0.01(N,M-N-1),多元线性关系显著,则多元线性回归方程(1)成立;如果F≤F0.01(N,M-N-1),多元线性关系不显著,则略去相关性分析中Pearson相关系数绝对值较小的测井曲线,重新进行多元线性回归,直到F>F0.01(n,M-n-1)成立为止,则通过测井计算的有机碳含量的公式为:
[0071] TOClog=a×GR+b×DEN+c×CNL+d (2)
[0072] 2)计算孔隙度:
[0073] 首先建立岩石体积模型,通过体积密度DEN得到体积方程:
[0074]
[0075]3
[0076] 式中,为测井计算的孔隙度,%;ρm为地层骨架密度,g/cm ;ρb为地层体积密度,即体积密度DEN,g/cm3;VTOC为有机碳含量体积百分比,vol%;wTOC为测井计算的有机碳含量TOClog,通过公式(2)获得,wt%;ρfl为地层流体密度,g/cm3;k为校正因子,一般取1。
[0077] 根据以上体积方程可以得到孔隙度的计算公式为:
[0078]
[0079] 将公式(5)计算得到的孔隙度 与步骤一中岩心实测的页岩样品孔隙度 进行线性回归:
[0080]
[0081]
[0082] 式中,e、f为线性回归系数,可以通过SPSS软件进行线性回归后得到; 为校正后的测井计算孔隙度,%。
[0083] 3)计算含气饱和度:
[0084] 成熟页岩不含油,因此含气饱和度和含水饱和度之和等于100%,故一旦含水饱和度确定,就可以通过含水饱和度计算含气饱和度。
[0085] 页岩含水饱和度可以根据阿尔奇公式计算:
[0086]
[0087] 式中,Sw为含水饱和度,%;Rw为地层水电阻率,ohmm; 为测井计算的孔隙度,通过公式(5)获得,%;Rt为地层电阻率,ohmm,页岩储层中常取与页岩相邻或相近的灰岩或砂岩地层水电阻率;m为胶结指数;n为饱和度指数。
[0088] 将公式(8)计算得到的含水饱和度Sw与步骤一中岩心实测的页岩样品含水饱和度Sw-core进行线性回归:
[0089] Sw-core=g×Sw+h (9)
[0090] Sw-log=g×Sw+h (10)
[0091] 式中,Sw-log为经过岩心校正的测井计算的含水饱和度,%;g、h为线性回归系数,可以通过SPSS软件进行线性回归后得到。
[0092] 页岩含气饱和度为:
[0093] Sg=1-Sw-log (11)
[0094] 步骤三:计算页岩含气量:
[0095] 在计算页岩储层原地资源量时通常使用吸附气和游离气之和来表征总含气量,故需要对吸附气含量和游离气含量分别进行计算:
[0096] 1)吸附气含量:
[0097] 根据地层实际温度、压力和有机碳含量对页岩样品进行校正,具体包括利用公式(12)进行兰格缪尔压力校正,利用公式(13)进行兰格缪尔体积校正,以及利用公式(14)进行有机碳含量校正:
[0098] Vlt=33.3357Vl×1.0062Ti-T (12)
[0099] Plt=145.1379Pl×1.0116T-Ti (13)
[0100]
[0101] 式中,Vlt为储层温度校正的兰格缪尔体积,m3/t;Plt为储层温度校正的兰格缪尔压力,MPa;T为储层温度,℃;Ti为等温吸附实验温度,℃;Vl为页岩样品兰格缪尔体积,m3/t;Pl为页岩样品兰格缪尔压力,MPa;Vlc为经过储层温度和有机碳含量校正的兰格缪尔体积,m3/t。
[0102] 根据兰格缪尔方程计算页岩储层吸附气量:
[0103]3
[0104] 式中,Ga为吸附气含量,m /t;P为储层压力,MPa,根据地层压力梯度计算,正常压力梯度情况下P=9.8×H/1000;H为储层深度,m。
[0105] 2)游离气含量:
[0106] 游离气量是页岩孔隙空间中的含气量,可以通过体积模型计算:
[0107]
[0108] 式中,Gf为游离气量,m3/t;Bg为气体压缩系数,m3/m3,根据地层温压条件确定;3
为测井计算的孔隙度,%;Sg为页岩含气饱和度,%;ρb为地层体积密度,g/cm 。
[0109] 由于吸附气的存在会影响游离气的容纳空间,故需减去吸附气所占的体积空间,因此游离气实际含量为:
相关专利内容
标题 发布/更新时间 阅读量
一种砂岩制砂工艺 2020-05-13 152
砂岩制砂工艺 2020-05-11 890
菱镁砂岩浮雕板 2020-05-11 648
一种平稳可靠型砂岩机 2020-05-12 241
一种砂岩储层酸化方法 2020-05-12 9
菱镁砂岩保温外墙板 2020-05-12 986
砂岩板材 2020-05-11 249
砂岩装饰板的制备方法 2020-05-13 191
砂岩环保吸声板 2020-05-11 982
砂岩磨盘 2020-05-11 694
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈