CN106442268B - 一种页岩介孔孔径分布的检测方法 - Google Patents
一种页岩介孔孔径分布的检测方法 Download PDFInfo
- Publication number
- CN106442268B CN106442268B CN201610941737.2A CN201610941737A CN106442268B CN 106442268 B CN106442268 B CN 106442268B CN 201610941737 A CN201610941737 A CN 201610941737A CN 106442268 B CN106442268 B CN 106442268B
- Authority
- CN
- China
- Prior art keywords
- shale
- pore
- nitrogen adsorption
- pressure
- micropore
- 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
- 239000011148 porous material Substances 0.000 title claims abstract description 112
- 238000001514 detection method Methods 0.000 title claims abstract description 25
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims abstract description 254
- 229910052757 nitrogen Inorganic materials 0.000 claims abstract description 129
- 238000001179 sorption measurement Methods 0.000 claims abstract description 71
- 238000002336 sorption--desorption measurement Methods 0.000 claims abstract description 32
- 238000002474 experimental method Methods 0.000 claims abstract description 20
- 238000003795 desorption Methods 0.000 claims description 36
- 238000000034 method Methods 0.000 claims description 26
- 239000007788 liquid Substances 0.000 claims description 23
- 239000007789 gas Substances 0.000 claims description 16
- 238000010521 absorption reaction Methods 0.000 claims description 15
- 239000003463 adsorbent Substances 0.000 claims description 13
- 238000002203 pretreatment Methods 0.000 claims description 7
- 238000007872 degassing Methods 0.000 claims description 6
- 230000018044 dehydration Effects 0.000 claims description 6
- 238000006297 dehydration reaction Methods 0.000 claims description 6
- 239000002156 adsorbate Substances 0.000 claims description 5
- 230000011218 segmentation Effects 0.000 claims description 3
- 239000011435 rock Substances 0.000 claims 1
- 230000007547 defect Effects 0.000 abstract description 4
- 238000012360 testing method Methods 0.000 abstract description 4
- 238000004458 analytical method Methods 0.000 description 9
- 230000000052 comparative effect Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 230000018109 developmental process Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- QJGQUHMNIGDVPM-UHFFFAOYSA-N nitrogen group Chemical group [N] QJGQUHMNIGDVPM-UHFFFAOYSA-N 0.000 description 4
- 238000005354 coacervation Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000003079 shale oil Substances 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005345 coagulation Methods 0.000 description 1
- 230000015271 coagulation Effects 0.000 description 1
- 238000009833 condensation Methods 0.000 description 1
- 230000005494 condensation Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 239000003921 oil Substances 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- ZSDSQXJSNMTJDA-UHFFFAOYSA-N trifluralin Chemical compound CCCN(CCC)C1=C([N+]([O-])=O)C=C(C(F)(F)F)C=C1[N+]([O-])=O ZSDSQXJSNMTJDA-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
- G01N15/0893—Investigating volume, surface area, size or distribution of pores; Porosimetry by measuring weight or volume of sorbed fluid, e.g. B.E.T. method
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Solid-Sorbent Or Filter-Aiding Compositions (AREA)
Abstract
本发明提供了一种页岩介孔孔径分布的检测方法,包括以下步骤:a)将页岩待测样品进行低温低压氮气等温吸附实验,得到氮气等温吸附/脱附曲线;b)以微孔隙和微裂缝为孔隙形态,对得到的氮气等温吸附/脱附曲线进行分析,得到页岩介孔孔径分布。与现有技术相比,本发明提供的检测方法以微孔隙和微裂缝作为页岩孔隙形态,进一步通过对其氮气等温吸附/脱附曲线进行分析,得到相应孔隙形态下的页岩介孔孔径分布,克服了现有技术中仅将页岩中介孔视为单一孔形造成孔隙体积计算误差大的缺陷,检测结果对于实际中的页岩而言更加准确,从而实现科学合理地对页岩内部孔隙结构和形态进行正确评价。
Description
技术领域
本发明涉及非常规油气勘探开发技术领域,更具体地说,是涉及一种页岩介孔孔径分布的检测方法。
背景技术
随着非常规油气的大规模勘探开发,页岩油气成为当今最令人关注的非常规油气能源之一。在对页岩油气进行开发之前,为了正确评价页岩储层的结构与性质,需对页岩储层中的孔隙形态以及孔隙体积进行定量分析评价,尤其是作为页岩油气藏研究重点的介孔孔隙。准确的介孔孔隙信息可为页岩气赋存状态以及流动状态的研究提供关键的帮助,从而为提高采收率、降低开发成本及有效避免勘探开发过程中的各种损失提供基础数据。
页岩内存在大量微纳量级的孔隙及裂缝,为页岩气藏的开发提供了可能性,根据国际理论和应用化学联合会(IUPAC)的分类,将孔隙按照直径大小分为微孔(<2nm)、介孔(2nm~50nm)和宏孔(>50nm)。其中,介孔由于其介于微观尺度与宏观尺度之间的特点,成为页岩储层中研究赋存状态和流动状态的重点内容。
目前,现有的检测页岩介孔孔径分布的方法,通常将页岩中的孔隙单一地视为圆柱形孔隙进行分析;而对于实际中的页岩而言,其孔隙的形状是多变的,不仅包括圆柱形孔隙,还包括平板形孔隙、椭圆柱形孔隙等,甚至还可以是更为复杂的孔隙形状。因此,在现有的检测方法中,仅将孔隙视为圆柱形孔隙进行页岩介孔孔径分布的检测,会造成孔隙体积的较大偏差,无法合理准确地认识页岩内部的孔隙结构及形态,极大地影响了页岩储层的客观正确评价。
发明内容
有鉴于此,本发明提供了一种页岩介孔孔径分布的检测方法,本发明提供的检测方法能够得到更准确的页岩介孔孔径分布,从而实现科学合理地对页岩内部孔隙结构和形态进行正确评价。
本发明提供了一种页岩介孔孔径分布的检测方法,包括以下步骤:
a)将页岩待测样品进行低温低压氮气等温吸附实验,得到氮气等温吸附/脱附曲线;
b)以微孔隙和微裂缝为孔隙形态,对得到的氮气等温吸附/脱附曲线进行分析,得到页岩介孔孔径分布。
优选的,步骤a)中所述页岩待测样品的质量为0.5g~5g,粒径小于等于150μm。
优选的,步骤a)中所述低温低压氮气等温吸附实验的温度为-200℃~-190℃,压力为0.005MPa~0.1MPa。
优选的,所述步骤b)具体包括以下步骤:
b1)采用修正Kelvin方程对氮气等温吸附/脱附曲线中的数据进行处理,得到处理后的数据;
b2)基于BJH模型分别对氮气等温吸附曲线和脱附曲线进行分析,得到孔隙形态为微孔隙和微裂缝的基本方程;
b3)采用步骤b2)得到的基本方程并运用中值定理和积分区间分割,对步骤b1)得到的处理后的数据进行分析,得到页岩介孔孔径分布。
优选的,步骤b1)中所述修正Kelvin方程为:
r=rk+t
式中,rk为Kelvin半径,c为孔隙类型常数,γ为表面张力,νm为吸附质的摩尔体积,φ为弯曲界面与孔壁之间的接触角,R为普适摩尔气体常数,T为液氮温度,x为氮气饱和压力下的相对压力,r为孔隙半径,t为吸附层厚度;
所述孔隙形态为微孔隙的模型中,氮气等温吸附曲线的c=2,氮气等温脱附曲线的c=1;
所述孔隙形态为微裂缝的模型中,氮气等温吸附曲线的c=1,氮气等温脱附曲线的c=0。
优选的,所述修正Kelvin方程中,吸附层厚度t采用Halsey方程得到;
所述Halsey方程为:
式中,
优选的,步骤b1)中所述进行处理的过程采用线性插值的方法得到吸附曲线和脱附曲线之间相对压力数据及其吸附量数据。
优选的,步骤b2)中所述孔隙形态为微孔隙和微裂缝的基本方程为:
式中,xi为相对压力值,xi-1为xi相邻压力点的相对压力值,rk1为微孔隙在吸附时相对压力对应的Kelvin半径,rk2为微孔隙在脱附时以及微裂缝在吸附时相对压力对应的Kelvin半径/裂缝宽度,Vc为微孔隙的孔径分布,Vs为微裂缝的孔径分布,Δνad为吸附曲线上相邻压力点之间对应的吸附量差值,Δνde为脱附曲线上相邻压力点之间对应的吸附量差值。
优选的,步骤a)中所述进行低温低压氮气等温吸附实验前,还包括:
对页岩待测样品进行前处理;
所述前处理具体为:
将页岩待测样品进行脱水和脱气。
优选的,步骤a)中所述脱水的温度为105℃~120℃,时间为12h~24h;所述脱气的真空度为小于等于1.01325×10-4MPa,时间为8h~12h。
本发明提供了一种页岩介孔孔径分布的检测方法,包括以下步骤:a)将页岩待测样品进行低温低压氮气等温吸附实验,得到氮气等温吸附/脱附曲线;b)以微孔隙和微裂缝为孔隙形态,对得到的氮气等温吸附/脱附曲线进行分析,得到页岩介孔孔径分布。与现有技术相比,本发明提供的检测方法以微孔隙和微裂缝作为页岩孔隙形态,进一步通过对其氮气等温吸附/脱附曲线进行分析,得到相应孔隙形态下的页岩介孔孔径分布,克服了现有技术中仅将页岩中介孔视为单一孔形造成孔隙体积计算误差大的缺陷,检测结果对于实际中的页岩而言更加准确,从而实现科学合理地对页岩内部孔隙结构和形态进行正确评价。
附图说明
图1为本发明得到页岩介孔孔径分布的分析过程;
图2为本发明实施例1得到的页岩待测样品的氮气等温吸附/脱附曲线;
图3为本发明实施例1和对比例得到的页岩介孔孔径分布图。
具体实施方式
下面将结合本发明实施例,对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明提供了一种页岩介孔孔径分布的检测方法,包括以下步骤:
a)将页岩待测样品进行低温低压氮气等温吸附实验,得到氮气等温吸附/脱附曲线;
b)以微孔隙和微裂缝为孔隙形态,对得到的氮气等温吸附/脱附曲线进行分析,得到页岩介孔孔径分布。
本发明首先将页岩待测样品进行低温低压氮气等温吸附实验,得到氮气等温吸附/脱附曲线。在本发明中,所述页岩待测样品的质量优选为0.5g~5g,更优选为1g~3g;所述页岩待测样品的粒径优选小于等于150μm,更优选为50μm~150μm。在本发明中,所述页岩待测样品均出自同一页岩,并且质量和粒径范围一致,能够消除非均质性所造成的离散误差,保证后续实验的准确性。
在本发明中,所述进行低温低压氮气等温吸附实验前,优选还包括:
对页岩待测样品进行前处理。在本发明中,所述前处理具体优选为:
将页岩待测样品进行脱水和脱气。本发明对所述脱水和脱气的方式没有特殊限制,采用本领域技术人员熟知的脱水和脱气的技术方案即可。在本发明中,所述脱水的温度优选为105℃~120℃,更优选为110℃~115℃;所述脱水的时间优选为12h~24h,更优选为16h~20h。在本发明中,所述脱气的真空度优选为小于等于1.01325×10-4MPa;所述脱气的时间优选为8h~12h,更优选为10h~11h。
完成所述前处理过程后,本发明将前处理后得到的页岩待测样品进行低温低压氮气等温吸附实验,得到氮气等温吸附/脱附曲线。在本发明中,所述低温低压氮气等温吸附实验能够得到压力增加或是减少的过程中页岩中孔隙氮气吸附量的具体变化,吸附/脱附过程中不同压力下对应的不同氮气吸附量构成氮气等温吸附/脱附曲线,即横坐标为相对压力x=p/pm,纵坐标为氮气吸附量ν(单位:cm3/g)。本发明对所述低温低压氮气等温吸附实验的设备没有特殊限制,采用本领域技术人员熟知的全自动比表面积和孔隙分析仪即可。
在本发明中,所述低温低压氮气等温吸附实验的温度优选为-200℃~-190℃,更优选为-196℃~-194℃。在本发明中,所述低温低压氮气等温吸附实验的压力优选为0.005MPa~0.1MPa。在本发明一个优选的实施例中,所述低温低压氮气等温吸附实验具体实验步骤为:将页岩待测样品放入全自动比表面积和孔隙分析仪中,在-196℃下,使压力从0.005MPa开始每增加一定压力记录其压力以及待体系平衡后该压力下的氮气吸附量,直至压力达到0.1MPa,得到氮气等温吸附曲线;再使压力从0.1MPa开始每降低一定压力记录其压力及待体系平衡后该压力下的氮气吸附量,直至压力达到0.005MPa,得到氮气等温脱附曲线。
得到页岩待测样品的氮气等温吸附/脱附曲线后,本发明以微孔隙和微裂缝为孔隙形态,对得到的氮气等温吸附/脱附曲线进行分析,得到页岩介孔孔径分布。在本发明中,上述分析过程优选具体包括以下步骤:
b1)采用修正Kelvin方程对氮气等温吸附/脱附曲线中的数据进行处理,得到处理后的数据;
b2)基于BJH模型分别对氮气等温吸附曲线和脱附曲线进行分析,得到孔隙形态为微孔隙和微裂缝的基本方程;
b3)采用步骤b2)得到的基本方程并运用中值定理和积分区间分割,对步骤b1)得到的处理后的数据进行分析,得到页岩介孔孔径分布。
本发明首先采用修正Kelvin方程对氮气等温吸附/脱附曲线中的数据进行处理,得到处理后的数据。在本发明中,所述修正Kelvin方程为:
r=rk+t
式中,rk为Kelvin半径,c为孔隙类型常数,γ为表面张力,νm为吸附质的摩尔体积,φ为弯曲界面与孔壁之间的接触角,R为普适摩尔气体常数,T为液氮温度,x为氮气饱和压力下的相对压力,r为孔隙半径,t为吸附层厚度。在本发明中,所述孔隙形态为微孔隙的模型中,氮气等温吸附曲线的c=2,氮气等温脱附曲线的c=1;所述孔隙形态为微裂缝的模型中,氮气等温吸附曲线的c=1,氮气等温脱附曲线的c=0。在本发明中,所述γ与吸附质性质相关。在本发明一个优选的实施例中,T=77.3K,所述液氮温度下γ=8.85×10-3N/m,所述νm在液氮平衡状态下为34.677cm3/mol,所述φ在计算临界半径时为0。
在本发明中,所述修正Kelvin方程中,吸附层厚度t采用Halsey方程得到;所述Halsey方程为:
式中,
在本发明中,所述修正Kelvin方程是关于压力-半径的方程,能够将氮气等温吸附/脱附曲线数据中的相对压力读数转化为对应的孔隙半径值。在本发明中,所述rk为Kelvin半径,并分为rk1(微孔隙在吸附时相对压力对应的Kelvin半径)和rk2(微孔隙在脱附时以及微裂缝在吸附时相对压力对应的Kelvin半径/裂缝宽度);所述r为孔隙半径,同理分为r1(微孔隙在吸附时相对压力对应的孔隙半径)和r2(微孔隙在脱附时以及微裂缝在吸附时相对压力对应的孔隙半径/裂缝宽度)。
在本发明中,氮气在毛细孔中发生吸附以及凝聚现象后,由气态转变为液态,可以将氮气吸附量通过状态方程及质量守恒计算,转化为填充进孔隙中的液氮体积,即孔隙中被液氮填充占用的体积,进而将氮气吸附量转化为孔隙体积分布,具体孔径范围为所述相对压力测量值对应的孔径范围。
在本发明中,所述采用修正Kelvin方程对氮气等温吸附/脱附曲线中的数据进行处理的过程,优选采用线性插值的方法得到吸附曲线和脱附曲线之间相对压力数据及其吸附量数据。在本发明中,由于同一孔径的孔隙根据修正Kelvin方程在吸附/脱附过程中对应不同的相对压力读数,故需利用线性插值的方法得到吸附曲线和脱附曲线两者之间对应的相对压力读数及其氮气吸附量读数,以方便后续分析。
本发明基于BJH模型分别对氮气等温吸附曲线和脱附曲线进行分析,得到孔隙形态为微孔隙和微裂缝的基本方程。在本发明中,基于BJH模型的毛细凝聚现象以及体积等效代换的原理,分别对吸附曲线和脱附曲线进行分析。吸附过程中,随着压力增大一定量时,氮气填充进页岩储集空间中所占体积分为三部分,其一为较大的微孔隙壁面上吸附层厚度增加所占的体积,其二为较小的微孔隙发生凝聚后凝聚液填满孔隙所占的体积,其三为微裂缝壁面上吸附层厚度增加所占的体积,值得说明的是,微裂缝在吸附过程中不会出现毛细凝聚现象;同理,脱附过程中,随着压力减小一定量时,氮气释放页岩储集空间所空出的体积分为四部分,其一为较大的微孔隙的壁面上吸附层厚度减小所占的体积,其二为较小的微孔隙发生解凝后释放出的体积,其三为较大的微裂缝壁面上吸附层厚度减小所占的体积,其四为较小的微裂缝发生解凝后释放出的体积。由上可见,利用吸附过程以及脱附过程中氮气吸附量的不同即可计算出不同孔隙形态的孔径分布。
在本发明中,所述孔隙形态为微孔隙和微裂缝的基本方程为:
式中,xi为相对压力值,xi-1为xi相邻压力点的相对压力值,rk1为微孔隙在吸附时相对压力对应的Kelvin半径,rk2为微孔隙在脱附时以及微裂缝在吸附时相对压力对应的Kelvin半径/宽度,Vc为微孔隙的孔径分布,Vs为微裂缝的孔径分布,Δνad为吸附曲线上相邻压力点之间对应的吸附量差值,Δνde为脱附曲线上相邻压力点之间对应的吸附量差值。
在本发明中,所述Δνad、Δνde、xi和xi-1由氮气等温吸附/脱附曲线数据直接得到;Vc和Vs为需要得到的页岩介孔孔径分布。同时,rk1、rk2和t均在步骤b1)中得到。在此基础上,本发明采用步骤b2)得到的基本方程并运用中值定理和积分区间分割,对步骤b1)得到的处理后的数据进行分析,得到页岩介孔孔径分布。
在本发明中,所述分析过程如图1所示,图1为本发明得到页岩介孔孔径分布的分析过程。
本发明提供了一种页岩介孔孔径分布的检测方法,包括以下步骤:a)将页岩待测样品进行低温低压氮气等温吸附实验,得到氮气等温吸附/脱附曲线;b)以微孔隙和微裂缝为孔隙形态,对得到的氮气等温吸附/脱附曲线进行分析,得到页岩介孔孔径分布。与现有技术相比,本发明提供的检测方法以微孔隙和微裂缝作为页岩孔隙形态,进一步通过对其氮气等温吸附/脱附曲线进行分析,得到相应孔隙形态下的页岩介孔孔径分布,克服了现有技术中仅将页岩中介孔视为单一孔形造成孔隙体积计算误差大的缺陷,检测结果对于实际中的页岩而言更加准确,从而实现科学合理地对页岩内部孔隙结构和形态进行正确评价。
为了进一步说明本发明,下面通过以下实施例进行详细说明。以下实施例和对比例中所用的页岩均为清江盆地曲江向斜地区曲页1井625.42米处泥页岩。
实施例1
(1)将页岩进行加工制成粒径为50μm~150μm的粉末样,然后将上述粉末样在110℃下脱水18h,再在真空度为1.01325×10-4MPa下脱气10h,得到页岩待测样品。
(2)将页岩待测样品进行低温低压氮气等温吸附实验,得到氮气等温吸附/脱附曲线:将页岩待测样品放入Tristar II 3020M型全自动比表面积和孔隙分析仪中,在-196℃下,使压力从0.005MPa开始每增加一定压力记录其压力以及待体系平衡后该压力下的氮气吸附量,直至压力达到0.1MPa,得到氮气等温吸附曲线;再使压力从0.1MPa开始每降低一定压力记录其压力及待体系平衡后该压力下的氮气吸附量,直至压力达到0.005MPa,得到氮气等温脱附曲线。参见图2所示,图2为本发明实施例1得到的页岩待测样品的氮气等温吸附/脱附曲线。
(3)采用修正Kelvin方程对氮气等温吸附/脱附曲线中的数据进行处理,得到处理后的数据,方便后续分析:
所述修正Kelvin方程为:
r=rk+t
式中,rk为Kelvin半径,c为孔隙类型常数,γ为表面张力,νm为吸附质的摩尔体积,φ为弯曲界面与孔壁之间的接触角,R为普适摩尔气体常数,T为液氮温度,x为氮气饱和压力下的相对压力,r为孔隙半径,t为吸附层厚度。所述孔隙形态为微孔隙的模型中,氮气等温吸附曲线的c=2,氮气等温脱附曲线的c=1;所述孔隙形态为微裂缝的模型中,氮气等温吸附曲线的c=1,氮气等温脱附曲线的c=0;T=77.3K,所述液氮温度下γ=8.85×10-3N/m,所述νm在液氮平衡状态下为34.677cm3/mol,所述φ在计算临界半径时为0;吸附层厚度t采用Halsey方程得到;所述Halsey方程为:
式中,
而氮气在毛细孔中发生吸附以及凝聚现象后,由气态转变为液态,故可将氮气吸附量通过状态方程转化为填充进孔隙中的液氮体积,即孔隙中被液氮填充占用的体积:将氮气等温吸附/脱附曲线中的氮气吸附量代入质量守恒方程及理想气体和液体状态方程可得ρgVg=ρlVl,其中ρg=1.25g/L,为氮气在77.3K下的密度,ρl=808.3kg/m3,为液氮体积。
具体为:脱附曲线上相对压力为0.97563,通过修正Kelvin方程得到相对应的微孔隙的孔隙半径值为29.42nm,在吸附曲线上对应的相对压力为0.96605,位于测量压力点区间(0.94490,0.97297)之间,利用线性插值的方法,从相对应的测量氮气吸附量区间(14.333cm3/g,17.603cm3/g)得到相对压力0.96605下对应的氮气吸附量为16.797cm3/g。
(4)基于BJH模型分别对氮气等温吸附曲线和脱附曲线进行分析,得到孔隙形态为微孔隙和微裂缝的基本方程。所述孔隙形态为微孔隙和微裂缝的基本方程为:
式中,xi为相对压力值,xi-1为xi相邻压力点的相对压力值,rk1为微孔隙在吸附时相对压力对应的Kelvin半径,rk2为微孔隙在脱附时以及微裂缝在吸附时相对压力对应的Kelvin半径/宽度,Vc为微孔隙的孔径分布,Vs为微裂缝的孔径分布,Δνad为吸附曲线上相邻压力点之间对应的吸附量差值,Δνde为脱附曲线上相邻压力点之间对应的吸附量差值。
(5)采用上述基本方程并运用中值定理和积分区间分割,对步骤(3)得到的处理后的数据进行分析,得到页岩介孔孔径分布。运用中值定理和积分区间分割,可得:
将式(IV)和式(V)代入到式(II),通过递推,所需要的所有ΔVi可由递推法逐级分析出来,具体分析过程如图1所示,得到页岩介孔孔径分布,即两种介孔孔形——微孔隙和微裂缝的孔径分布。
本发明实施例1得到页岩介孔孔径分布的数据结果参见表1。
表1本发明实施例1和对比例得到的页岩介孔孔径分布数据
本发明实施例1得到的页岩介孔孔径分布图参见图3。
值得说明的是,由于拉伸强度效应(Tensile strength effect/TSE),计算出的微裂缝在约2nm处会出现假峰,假峰的出现是不符合实际情况的,且假峰会使该孔径及以下的孔隙体积计算出现误差。因此,最终得到的孔径分布图中截取了孔径大于2nm结果可靠的部分。
对比例
以圆柱形孔隙为孔隙形态,对将实施例1得到的氮气等温吸附/脱附曲线进行分析,得到页岩介孔孔径分布:
(1)将氮气等温吸附/脱附曲线数据中的相对压力读数代入修正Kelvin方程及Halsey方程中,获得其对应的孔隙半径的值,如等温吸附曲线中对应相对压力为0.95893,将其代入Kelvin方程及Halsey方程中,得到对应的Kelvin半径rk(x)=22.73nm,对应的吸附层厚度t(x)=1.74nm,而对应的孔隙半径为两者之和,即r=rk+t=24.47nm。
(2)将等温吸附实验数据中的氮气吸附量代入质量守恒方程及理想气体和液体的状态方程可得ρgVg=ρlVl,其中ρg=1.25g/L,为氮气在77.3K下的密度,ρl=808.3kg/m3为液氮密度,则可将氮气吸附量转化为相对应的液氮体积,如等温吸附曲线中对应氮气吸附量为Vg=15.968cm3/g,转化后的液氮体积即为Vl=0.0246cm3/g。
(3)采用孔隙形态为圆柱形孔隙的基本方程:
将式(IV)和式(V)代入到式(VI),通过递推,所需要的所有ΔVi可由递推法逐级分析出来,具体为:最高压力区间(x1,x0,)为(0.97297,0.98137),更高的压力区间由于吸附层的厚度变化极小,可忽略。而对应的氮气吸附量区间为(17.603cm3/g,19.434cm3/g),已知转化为液氮体积的氮气吸附变化量ΔV1=2.832×10-3cm3/g,得到对应孔隙半径区间为(36.77nm,50.30nm),则进一步得到该半径区间内的孔隙体积为0.00372cm3/g;而次高压力区间(x2,x1)为(0.94490,0.97297),对应的氮气吸附量区间为(14.333cm3/g,17.603cm3/g),转化为液氮体积的氮气吸附变化量ΔV2=5.058×10-3cm3/g,对应的孔隙半径区间为(18.39nm,36.77nm),则分析得到该半径区间内的孔隙体积为0.00653cm3/g;以此类推可以得到孔径分布的所有ΔVi。
对比例得到页岩介孔孔径分布的数据结果参见表1,页岩介孔孔径分布图参见图3。
通过比较可知,本发明提供的检测方法以微孔隙和微裂缝作为页岩孔隙形态,进一步通过对其氮气等温吸附/脱附曲线进行分析,得到相应孔隙形态下的页岩介孔孔径分布,克服了现有技术中仅将页岩中介孔视为单一孔形造成孔隙体积计算误差大的缺陷,检测结果对于实际中的页岩而言更加准确,从而实现科学合理地对页岩内部孔隙结构和形态进行正确评价。
所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (8)
1.一种页岩介孔孔径分布的检测方法,包括以下步骤:
a)将页岩待测样品进行低温低压氮气等温吸附实验,得到氮气等温吸附/脱附曲线;
b)以微孔隙和微裂缝为孔隙形态,对得到的氮气等温吸附/脱附曲线进行分析,得到页岩介孔孔径分布;
所述步骤b)具体包括以下步骤:
b1)采用修正Kelvin方程对氮气等温吸附/脱附曲线中的数据进行处理,得到处理后的数据;
所述修正Kelvin方程为:
式中,rk为Kelvin半径,c为孔隙类型常数,γ为表面张力,νm为吸附质的摩尔体积,φ为弯曲界面与孔壁之间的接触角,R为普适摩尔气体常数,T为液氮温度,x为氮气饱和压力下的相对压力,r为孔隙半径,t为吸附层厚度;
所述孔隙形态为微孔隙的模型中,氮气等温吸附曲线的c=2,氮气等温脱附曲线的c=1;
所述孔隙形态为微裂缝的模型中,氮气等温吸附曲线的c=1,氮气等温脱附曲线的c=0;
b2)基于BJH模型分别对氮气等温吸附曲线和脱附曲线进行分析,得到孔隙形态为微孔隙和微裂缝的基本方程;
b3)采用步骤b2)得到的基本方程并运用中值定理和积分区间分割,对步骤b1)得到的处理后的数据进行分析,得到页岩介孔孔径分布。
2.根据权利要求1所述的检测方法,其特征在于,步骤a)中所述页岩待测样品的质量为0.5g~5g,粒径小于等于150μm。
3.根据权利要求1所述的检测方法,其特征在于,步骤a)中所述低温低压氮气等温吸附实验的温度为-200℃~-190℃,压力为0.005MPa~0.1MPa。
4.根据权利要求1所述的检测方法,其特征在于,所述修正Kelvin方程中,吸附层厚度t采用Halsey方程得到;
所述Halsey方程为:
式中,
5.根据权利要求1所述的检测方法,其特征在于,步骤b1)中所述进行处理的过程采用线性插值的方法得到吸附曲线和脱附曲线之间相对压力数据及其吸附量数据。
6.根据权利要求1所述的检测方法,其特征在于,步骤b2)中所述孔隙形态为微孔隙和微裂缝的基本方程为:
式中,xi为相对压力值,xi-1为xi相邻压力点的相对压力值,t(xi)为相对压力值为xi的吸附层厚度,t(xi-1)为相对压力值为xi-1的吸附层厚度,rk1为微孔隙在吸附时相对压力对应的Kelvin半径,rk2为微孔隙在脱附时以及微裂缝在吸附时相对压力对应的Kelvin半径/裂缝宽度,Vc为微孔隙的孔径分布,Vs为微裂缝的孔径分布,Δνad为吸附曲线上相邻压力点之间对应的吸附量差值,Δνde为脱附曲线上相邻压力点之间对应的吸附量差值。
7.根据权利要求1~6任一项所述的检测方法,其特征在于,步骤a)中所述进行低温低压氮气等温吸附实验前,还包括:
对页岩待测样品进行前处理;
所述前处理具体为:
将页岩待测样品进行脱水和脱气。
8.根据权利要求7所述的检测方法,其特征在于,步骤a)中所述脱水的温度为105℃~120℃,时间为12h~24h;所述脱气的真空度为小于等于1.01325×10-4MPa,时间为8h~12h。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610941737.2A CN106442268B (zh) | 2016-10-31 | 2016-10-31 | 一种页岩介孔孔径分布的检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610941737.2A CN106442268B (zh) | 2016-10-31 | 2016-10-31 | 一种页岩介孔孔径分布的检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106442268A CN106442268A (zh) | 2017-02-22 |
CN106442268B true CN106442268B (zh) | 2018-11-27 |
Family
ID=58178275
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610941737.2A Active CN106442268B (zh) | 2016-10-31 | 2016-10-31 | 一种页岩介孔孔径分布的检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106442268B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108802073B (zh) * | 2017-05-05 | 2021-03-12 | 中国石油化工股份有限公司 | 一种基于数字岩心的岩石电性参数获取方法及装置 |
CN107991215B (zh) * | 2017-11-23 | 2020-02-07 | 中南大学 | 大尺寸低渗岩样天然孔径与比表面积的测试方法 |
CN111537416B (zh) * | 2020-04-10 | 2021-10-22 | 中国地质大学(北京) | 一种纳米级孔隙结构变化的判断方法及应用 |
CN111664825B (zh) * | 2020-06-10 | 2021-07-09 | 西安石油大学 | 一种提高孔结构参数测定准确性的方法 |
CN111965067B (zh) * | 2020-08-19 | 2023-11-03 | 东莞东阳光科研发有限公司 | 一种循环后三元材料微裂纹的表征方法 |
CN114136855B (zh) * | 2020-09-03 | 2024-05-28 | 中国石油化工股份有限公司 | 一种判别页岩孔隙连通性的方法、存储介质和计算机设备 |
CN114428043B (zh) * | 2020-09-18 | 2024-06-21 | 中国石油化工股份有限公司 | 多孔介质孔径分布表征方法及电子设备 |
CN113484220B (zh) * | 2021-07-26 | 2022-07-22 | 中国石油大学(北京) | 确定页岩有机质和无机质分布的方法、装置及电子设备 |
CN117805000B (zh) * | 2024-03-01 | 2024-06-18 | 江苏龙城精锻集团有限公司 | 针对氢化环境下材料内部微孔隙的检测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002042744A1 (en) * | 2000-11-21 | 2002-05-30 | Akzo Nobel N.V. | Method of analysing microporous material |
CN105043957A (zh) * | 2015-07-06 | 2015-11-11 | 成都理工大学 | 通过泥页岩等温吸附曲线分类判断孔隙结构的方法 |
CN105974092A (zh) * | 2016-07-08 | 2016-09-28 | 重庆科技学院 | 一种致密储层孔喉全尺度表征及分析方法 |
-
2016
- 2016-10-31 CN CN201610941737.2A patent/CN106442268B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002042744A1 (en) * | 2000-11-21 | 2002-05-30 | Akzo Nobel N.V. | Method of analysing microporous material |
CN105043957A (zh) * | 2015-07-06 | 2015-11-11 | 成都理工大学 | 通过泥页岩等温吸附曲线分类判断孔隙结构的方法 |
CN105974092A (zh) * | 2016-07-08 | 2016-09-28 | 重庆科技学院 | 一种致密储层孔喉全尺度表征及分析方法 |
Non-Patent Citations (5)
Title |
---|
Pore Size Distribution Analysis of Microporous Carbons:A Density Functional Theory Approach;Christian Lastoskie et al.;《The Journal of Physical Chemistry》;19931231;第97卷(第18期);第4786-4796页 * |
Pore structure characterization of different rank coals using gas adsorption and scanning electron microscopy;Nie Bsisheng et al.;《Fuel》;20150623(第158期);摘要以及第2节 * |
The effects of pore geometry on adsorption equilibrium in shale formations and coal-beds: Lattice density functional theory study;Ali Qajar et al.;《Fuel》;20151003(第163期);第206页第2节第1-2段,表2,第212页左栏第2段 * |
压汞法和气体吸附法研究富有机质页岩孔隙特征;田华 等;《石油学报》;20120531;第33卷(第3期);第1.2、2.3节 * |
多孔活性炭孔径分布的表征;张超 等;《多孔活性炭孔径分布的表征》;20060430;第22卷(第2期);第187-192页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106442268A (zh) | 2017-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106442268B (zh) | 一种页岩介孔孔径分布的检测方法 | |
Schlumberger et al. | Characterization of hierarchically ordered porous materials by physisorption and mercury porosimetry—a tutorial review | |
Liu et al. | Determination of the absolute adsorption/desorption isotherms of CH4 and n-C4H10 on shale from a nano-scale perspective | |
CN106979917B (zh) | 一种页岩气储层孔隙结构的表征方法以及页岩气储层的评价方法 | |
Kong et al. | Adsorption/desorption isotherms of CH4 and C2H6 on typical shale samples | |
CN104677771B (zh) | 一种基于磁悬浮重量法确定页岩孔隙度的方法 | |
Wang et al. | Investigation of methane desorption and its effect on the gas production process from shale: experimental and mathematical study | |
CN110672813B (zh) | 一种页岩含气量计算方法 | |
CN111337410B (zh) | 一种确定页岩孔径分布的方法、装置、设备及系统 | |
Hu et al. | Influence of reservoir primary water on shale gas occurrence and flow capacity | |
CN105486621A (zh) | 一种页岩孔径分布测试方法 | |
CN103115844A (zh) | 一种煤页岩等温吸附/解吸曲线的测定方法 | |
CN108460219A (zh) | 基于超临界单分子层吸附模型计算页岩吸附气含气量方法 | |
CN108801879B (zh) | 一种页岩基质颗粒孔隙率和渗透率一体化测量系统及方法 | |
CN110702585B (zh) | 一种岩石压缩系数计算方法 | |
CN102980843A (zh) | 一种焦炭气孔特征的检测方法 | |
CN109342297A (zh) | 基于压汞实验的煤中孔隙校正方法 | |
Zhang et al. | The effect of water vapor on methane adsorption in the nanopores of shale | |
CN110032744A (zh) | 一种非常规油气储量计算方法及系统 | |
CN109696380A (zh) | 一种评价微孔材料气体分离性能的方法及装置 | |
CN111664825B (zh) | 一种提高孔结构参数测定准确性的方法 | |
CN203241324U (zh) | 一种页岩气体渗透率测定仪 | |
CN105445159A (zh) | 一种得到孔径分布曲线及样品比表面积的方法 | |
CN110274863A (zh) | 一种页岩油富集孔径下限的确定方法 | |
Harrison et al. | Characterization of superactivated carbons |
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 |