CN111122408A - 一种孔径介于0.3-10nm的页岩有机孔隙度估算方法 - Google Patents

一种孔径介于0.3-10nm的页岩有机孔隙度估算方法 Download PDF

Info

Publication number
CN111122408A
CN111122408A CN201911335718.5A CN201911335718A CN111122408A CN 111122408 A CN111122408 A CN 111122408A CN 201911335718 A CN201911335718 A CN 201911335718A CN 111122408 A CN111122408 A CN 111122408A
Authority
CN
China
Prior art keywords
shale
pore volume
pore
pore diameter
organic
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.)
Pending
Application number
CN201911335718.5A
Other languages
English (en)
Inventor
何陈诚
何生
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Geosciences
Original Assignee
China University of Geosciences
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Geosciences filed Critical China University of Geosciences
Priority to CN201911335718.5A priority Critical patent/CN111122408A/zh
Publication of CN111122408A publication Critical patent/CN111122408A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/088Investigating volume, surface area, size or distribution of pores; Porosimetry
    • G01N15/0893Investigating 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

本发明公开了一种孔径介于0.3‑10nm的页岩有机孔隙度估算方法,包括以下步骤:通过页岩样品二氧化碳和氮气吸附测试,利用密度函数理论DFT模型计算页岩样品内孔径为0.3‑10nm的孔隙体积;建立页岩样品孔径介于0.3‑10nm的孔隙体积与TOC含量的线性正相关关系;估算0.3‑10nm有机孔隙体积以及占总孔体积的比值,结合页岩总孔隙度估算页岩孔径介于0.3‑10nm的页岩有机孔隙度,该方法思路清晰、技术可靠、涉及的参数物理意义明确、估算方法和计算公式应用方便。

Description

一种孔径介于0.3-10nm的页岩有机孔隙度估算方法
技术领域
本发明涉及页岩气储层研究及页岩气勘探,具体涉及一种孔径介于0.3-10nm的页岩有机孔隙度估算方法。
背景技术
页岩气储层总有机孔隙度是指页岩气储层中全部有机质孔隙体积与页岩气储层总体积的比值,可用百分比表示。页岩气储层中有机质孔隙或称有机孔隙是页岩有机质在热演化裂解生气过程中形成的孔隙,有机孔隙是页岩气储层中广泛存在且是最重要的孔隙类型,有机孔隙发育在页岩气储层中的有机质颗粒内,以纳米级孔隙为主,有时也发育少量的微米级孔隙,是页岩气赋存的主要孔隙空间。页岩气储层有机孔隙度是页岩储气能力和含气性评价的重要参数,在页岩气的勘探开发过程中,页岩气储层有机孔隙度也是计算页岩气储量和制定开发方案的关键参数。页岩气储层中有机孔隙度计算是页岩气储层评价的重要内容,也是优选页岩气有利钻探目标的依据。
页岩气储层有机孔隙度无法通过仪器测定,而目前国内也还没有标准的页岩气储层有机孔隙度的计算公式。现在常用的是利用扫描电镜镜下观察有机孔面孔率,通过有机孔面孔率计算页岩有机孔隙度,但此方法仍然受限于场发射扫描电镜(FE-SEM)的分辨率,具有极高的分辨率FE-SEM可基本解决孔径>10nm有机孔隙的观察统计,但也不能完全观察统计孔径2-10nm的有机孔隙,更不可能观察统计孔径<2nm的有机微孔;纳米CT和FIB-SEM三维重构也仅能解决大于20nm有机孔隙的估算。页岩中0.3-10nm的有机孔隙数量众多、比表面积大、占总孔隙体积的比例约在10-30%,占全有机孔隙体积的比例约在20-50%,是页岩气吸附气和游离气富集的孔隙空间。
目前孔径介于0.3-10nm页岩有机孔隙度的估算是有机孔隙度评价的难点问题,仍未见到有效解决该问题的相关报道,故本专利提出一种孔径介于0.3-10nm页岩有机孔隙度的估算方法,该方法可解决页岩气储层中孔径介于0.3-10nm有机孔隙度的评价,该方法思路清晰、技术可靠、涉及的参数物理意义明确、估算方法和计算公式应用方便。
孔径介于0.3-10nm的页岩孔隙体积与TOC有非常好的正相关关系,且微孔-中孔主要发育在高成熟-过成熟页岩中,从而可以在沉积环境、岩石学、地球化学和成岩作用等地质条件相似的页岩层系,可根据页岩样品孔隙体积与TOC的关系获得孔径介于0.3-10nm的无机孔体积的变化范围,一般接近并取某一固定值,进而可根据TOC含量估算该部分有机孔体积。该方法计算误差主要来自孔径介于0.3-10nm的有机孔体积的计算。
发明内容
有鉴于此,本发明提供了一种孔径介于0.3-10nm的页岩有机孔隙度估算方法。
本发明提供一种孔径介于0.3-10nm的页岩有机孔隙度估算方法,包括以下步骤:
步骤1:利用二氧化碳吸附测试的密度函数理论DFT模型计算出页岩孔径介于0.3-1.4nm范围内的孔隙体积及其孔径分布;利用氮气吸附测试中的QSDFT模型计算页岩孔径介于1.4-30nm的孔隙体积及其孔径分布,再运用BJH模型计算出大于30nm以上的孔隙体积及其孔径分布,结合三类模型计算得到页岩气储层总孔隙体积V总孔以及孔径介于0.3-10nm的页岩孔隙体积V10
步骤2:建立孔径介于0.3-10nm的页岩孔隙体积V10与其页岩样品TOC值之间的线性正相关关系,得到介于0.3-10nm的页岩孔隙体积V10与其页岩样品TOC值的线性关系式,其中页岩样品TOC值由实验室实验仪器测试页岩样品得到;
步骤3:在地质条件包括沉积环境、岩石学、地球化学和成岩作用相似的页岩层系,根据页岩孔隙体积V10与其页岩样品TOC值的关系获得孔径介于0.3-10nm的无机孔隙体积,即步骤2线性关系式中TOC值为零时的孔隙体积,然后根据TOC值计算得到页岩孔径0.3-10nm范围的有机孔隙体积V10oc,以及V10oc与页岩气储层总孔隙体积V总孔的比值K10oc
步骤4:页岩孔径介于0.3-10nm的有机孔隙体积与总孔隙体积的比值K10oc,乘以页岩气储层总孔隙度
Figure BDA0002330875580000031
即为孔径介于0.3-10nm的页岩有机孔隙度
Figure BDA0002330875580000032
进一步地,步骤2中孔径介于0.3-10nm的页岩孔隙体积V10与其页岩样品TOC值的线性关系式如下所示:
V10=0.0018*TOC+0.0067cc/g
进一步地,步骤3中,有机孔隙体积V10oc以及有机孔隙体积与总孔隙体积的比值K10oc的计算公式分别如下:
V10oc=V10–0.0067cc/g
K10oc=V10oc/V总孔
其中,参数0.0067cc/g为无机孔隙体积,是依据孔径介于0.3-10nm的页岩孔隙体积V10与其页岩样品TOC值之间线性正相关关系得出的,即TOC值为零时得到的值;V10oc为有机孔隙体积,V总孔为页岩气储层总孔隙体积,V10为孔径介于0.3-10nm的页岩孔隙体积。
进一步地,步骤4中,孔径介于0.3-10nm的页岩有机孔隙度
Figure BDA0002330875580000041
计算公式为:
Figure BDA0002330875580000042
其中,K10oc为页岩孔径介于0.3-10nm的有机孔隙体积与页岩气储层总孔隙体积V总孔的比值,
Figure BDA0002330875580000043
为页岩气储层总孔隙度,由实验室实验仪器测试页岩样品得到。
本发明提供的技术方案带来的有益效果是:解决页岩气储层中孔径介于0.3-10nm有机孔隙度的评价,该方法思路清晰、技术可靠、涉及的参数物理意义明确、估算方法和计算公式应用方便。
附图说明
图1是本发明一种孔径介于0.3-10nm的页岩有机孔隙度估算方法流程图;
图2是本发明一种孔径介于0.3-10nm的页岩有机孔隙度估算方法中DFT模型结合BJH模型计算的页岩纳米孔隙的孔径分布图;
图3是本发明一种孔径介于0.3-10nm的页岩有机孔隙度估算方法DFT模型计算的孔径介于0.3-10nm孔隙体积与样品TOC关系图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地描述。
请参考图1,本发明的实施例提供了一种孔径介于0.3-10nm的页岩有机孔隙度估算方法。
步骤1:利用二氧化碳吸附测试的密度函数理论DFT模型计算出页岩孔径介于0.3-1.4nm范围内的孔隙体积及其孔径分布;利用氮气吸附测试中的QSDFT模型计算页岩孔径介于1.4-30nm的孔隙体积及其孔径分布,再运用BJH模型计算出大于30nm以上的孔隙体积及其孔径分布,结合三类模型计算得到页岩气储层总孔隙体积V总孔以及孔径介于0.3-10nm的页岩孔隙体积V10,见图2和表1:
表1DFT模型结合BJH模型计算的孔径介于0.3-10nm孔隙体积和有机孔体积计算表
Figure BDA0002330875580000051
步骤2:建立孔径介于0.3-10nm的页岩孔隙体积V10与其页岩样品TOC值之间的线性正相关关系,得到介于0.3-10nm的页岩孔隙体积V10与其页岩样品TOC值的线性关系式,其中页岩样品TOC值由实验室实验仪器测试岩石样品得到,参考图3,关系式如下所示:
V10=0.0018*TOC+0.0067cc/g
步骤3:在地质条件包括沉积环境、岩石学、地球化学和成岩作用相似的页岩层系,根据页岩孔隙体积V10与其页岩样品TOC值的关系获得孔径介于0.3-10nm的无机孔隙体积,即步骤2线性关系式中TOC值为零时的孔隙体积,如图1关系式中的0.0067cc/g,然后根据TOC值计算得到页岩孔径0.3-10nm范围的有机孔隙体积V10oc,以及V10oc与页岩气储层总孔隙体积V总孔的比值K10oc,计算公式分别如下:
V10oc=V10–0.0067cc/g
K10oc=V10oc/V总孔
其中,参数0.0067cc/g为无机孔体积,是依据孔径介于0.3-10nm内的页岩孔隙体积V10与其页岩样品TOC值之间线性正相关关系得出的,即TOC值为零时得到的值;V10oc为有机孔隙体积,V总孔为页岩气储层孔隙总体积,V10为孔径介于0.3-10nm内的页岩孔隙体积。
步骤4:页岩孔径介于0.3-10nm的有机孔隙体积与总孔隙体积的比值K10oc,乘以页岩气储层总孔隙度
Figure BDA0002330875580000061
即为孔径介于0.3-10nm的页岩有机孔隙度
Figure BDA0002330875580000062
如表1中的估算结果,计算公式为:
Figure BDA0002330875580000063
其中,K10oc为页岩孔径介于0.3-10nm的有机孔隙体积与页岩气储层总孔隙体积V总孔的比值,
Figure BDA0002330875580000064
为页岩气储层总孔隙度,由实验室实验仪器测试岩石样品得到。
由于二氧化碳和和氮气吸附实验的限制,页岩气储层孔隙总体积可能不完全,可结合压汞实验获得页岩气储层全孔径的孔隙体积,本实施例中的页岩气储层孔隙总体积是由二氧化碳和氮气吸附实验计算获得,此类计算孔径范围限于350nm以内的孔隙体积,与全孔体积有所差异。在观察同页岩层系的页岩样品的压汞实验数据中,孔径在350nm以上的孔隙体积基本集中在0.001-0.002cc/g,增加此类孔隙体积后,计算的孔径介于0.3-10nm的有机孔隙度会有10%以内的误差。
在不冲突的情况下,本文中上述实施例及实施例中的特征可以相互结合。以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (4)

1.一种孔径介于0.3-10nm的页岩有机孔隙度估算方法,其特征在于,包括以下步骤:
步骤1:利用二氧化碳吸附测试的密度函数理论DFT模型计算出页岩孔径介于0.3-1.4nm范围内的孔隙体积及其孔径分布;利用氮气吸附测试中的QSDFT模型计算页岩孔径介于1.4-30nm的孔隙体积及其孔径分布,再运用BJH模型计算出大于30nm以上的孔隙体积及其孔径分布,结合三类模型计算得到页岩气储层总孔隙体积V总孔以及孔径介于0.3-10nm的页岩孔隙体积V10
步骤2:建立孔径介于0.3-10nm的页岩孔隙体积V10与其页岩样品TOC值之间的线性正相关关系,得到介于0.3-10nm的页岩孔隙体积V10与其页岩样品TOC值的线性关系式,其中页岩样品TOC值由实验室实验仪器测试页岩样品得到;
步骤3:在地质条件包括沉积环境、岩石学、地球化学和成岩作用相似的页岩层系,根据页岩孔隙体积V10与其页岩样品TOC值的关系获得孔径介于0.3-10nm的无机孔隙体积,即步骤2线性关系式中TOC值为零时的孔隙体积,然后根据TOC值计算得到页岩孔径0.3-10nm范围的有机孔隙体积V10oc,以及V10oc与页岩气储层总孔隙体积V总孔的比值K10oc
步骤4:页岩孔径介于0.3-10nm的有机孔隙体积与总孔隙体积的比值K10oc,乘以页岩气储层总孔隙度
Figure FDA0002330875570000011
即为孔径介于0.3-10nm的页岩有机孔隙度
Figure FDA0002330875570000012
2.根据权利要求1所述的一种孔径介于0.3-10nm的页岩有机孔隙度估算方法,其特征在于,步骤2中孔径介于0.3-10nm的页岩孔隙体积V10与其页岩样品TOC值的线性关系式如下所示:
V10=0.0018*TOC+0.0067cc/g 。
3.根据权利要求1所述的一种孔径介于0.3-10nm的页岩有机孔隙度估算方法,其特征在于,步骤3中,有机孔隙体积V10oc以及有机孔隙体积与总孔隙体积的比值K10oc的计算公式分别如下:
V10oc=V10–0.0067cc/g
K10oc=V10oc/V总孔
其中,参数0.0067cc/g为无机孔隙体积,是依据孔径介于0.3-10nm的页岩孔隙体积V10与其页岩样品TOC值之间线性正相关关系得出的,即TOC值为零时得到的值;V10oc为有机孔隙体积,V总孔为页岩气储层总孔隙体积,V10为孔径介于0.3-10nm的页岩孔隙体积。
4.根据权利要求3所述的一种孔径介于0.3-10nm的页岩有机孔隙度估算方法,其特征在于,步骤4中,孔径介于0.3-10nm的页岩有机孔隙度
Figure FDA0002330875570000021
计算公式为:
Figure FDA0002330875570000022
其中,K10oc为页岩孔径介于0.3-10nm的有机孔隙体积与页岩气储层总孔隙体积V总孔的比值,
Figure FDA0002330875570000023
为页岩气储层总孔隙度,由实验室实验仪器测试页岩样品得到。
CN201911335718.5A 2019-12-23 2019-12-23 一种孔径介于0.3-10nm的页岩有机孔隙度估算方法 Pending CN111122408A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911335718.5A CN111122408A (zh) 2019-12-23 2019-12-23 一种孔径介于0.3-10nm的页岩有机孔隙度估算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911335718.5A CN111122408A (zh) 2019-12-23 2019-12-23 一种孔径介于0.3-10nm的页岩有机孔隙度估算方法

Publications (1)

Publication Number Publication Date
CN111122408A true CN111122408A (zh) 2020-05-08

Family

ID=70501096

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911335718.5A Pending CN111122408A (zh) 2019-12-23 2019-12-23 一种孔径介于0.3-10nm的页岩有机孔隙度估算方法

Country Status (1)

Country Link
CN (1) CN111122408A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111827969A (zh) * 2020-07-16 2020-10-27 长江大学 一种页岩气储量计算方法、设备及可读存储介质
CN112858136A (zh) * 2021-01-29 2021-05-28 中国地质大学(武汉) 一种页岩有机质孔隙结构的定量评价方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103344541A (zh) * 2013-07-08 2013-10-09 中国石油大学(华东) 一种泥页岩总孔隙度测定方法
CN103512838A (zh) * 2013-09-18 2014-01-15 中国石油大学(华东) 一种测定泥页岩储层不同孔径孔隙对孔隙度贡献的方法
CN105675468A (zh) * 2016-01-15 2016-06-15 山东科技大学 一种泥页岩有机质孔隙度的评价方法
CN106223941A (zh) * 2016-07-20 2016-12-14 中石化石油工程技术服务有限公司 基于测井资料的页岩气储层有机孔隙度确定方法
CN107194104A (zh) * 2017-06-07 2017-09-22 成都理工大学 一种采用测井曲线计算页岩储层中无机质孔隙度、总孔隙度以及组分的方法
CN108561126A (zh) * 2018-03-13 2018-09-21 中石化石油工程技术服务有限公司 一种确定页岩气储层有机孔隙度的简易方法
CN110487693A (zh) * 2018-05-14 2019-11-22 中国石油化工股份有限公司 一种确定泥页岩不同类型孔隙度的方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103344541A (zh) * 2013-07-08 2013-10-09 中国石油大学(华东) 一种泥页岩总孔隙度测定方法
CN103512838A (zh) * 2013-09-18 2014-01-15 中国石油大学(华东) 一种测定泥页岩储层不同孔径孔隙对孔隙度贡献的方法
CN105675468A (zh) * 2016-01-15 2016-06-15 山东科技大学 一种泥页岩有机质孔隙度的评价方法
CN106223941A (zh) * 2016-07-20 2016-12-14 中石化石油工程技术服务有限公司 基于测井资料的页岩气储层有机孔隙度确定方法
CN107194104A (zh) * 2017-06-07 2017-09-22 成都理工大学 一种采用测井曲线计算页岩储层中无机质孔隙度、总孔隙度以及组分的方法
CN108561126A (zh) * 2018-03-13 2018-09-21 中石化石油工程技术服务有限公司 一种确定页岩气储层有机孔隙度的简易方法
CN110487693A (zh) * 2018-05-14 2019-11-22 中国石油化工股份有限公司 一种确定泥页岩不同类型孔隙度的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
李军 等: "页岩气储层"四孔隙"模型建立及测井定量表征方法", 《石油与天然气地质》 *
陈曼霏等: "涪陵页岩气田平桥区块页岩气储层有机质孔发育特征", 《石油学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111827969A (zh) * 2020-07-16 2020-10-27 长江大学 一种页岩气储量计算方法、设备及可读存储介质
CN112858136A (zh) * 2021-01-29 2021-05-28 中国地质大学(武汉) 一种页岩有机质孔隙结构的定量评价方法

Similar Documents

Publication Publication Date Title
Wei et al. Nanopore structure characterization for organic-rich shale using the non-local-density functional theory by a combination of N2 and CO2 adsorption
Liu et al. Nanoscale pore structure characterization of the Bakken shale in the USA
Sun et al. Nanoscale pore characteristics of the Lower Cambrian Niutitang Formation Shale: a case study from Well Yuke# 1 in the Southeast of Chongqing, China
CN103512838B (zh) 一种测定泥页岩储层不同孔径孔隙对孔隙度贡献的方法
CN104794351A (zh) 一种页岩气储层总有机孔隙度的计算方法及装置
CN111122408A (zh) 一种孔径介于0.3-10nm的页岩有机孔隙度估算方法
Wang et al. Exploration potential of different lithofacies of deep marine shale gas systems: Insight into organic matter accumulation and pore formation mechanisms
CN110487693A (zh) 一种确定泥页岩不同类型孔隙度的方法
CN107132171B (zh) 一种基于压汞-氮吸附联测数据确定致密储层孔径分布的方法
Sun et al. Pore connectivity and water accessibility in Upper Permian transitional shales, southern China
CN111007233A (zh) 一种分析页岩微观孔隙中甲烷-二氧化碳运动行为的方法
CN110470584B (zh) 一种评价渗吸和水锁综合效应的方法
GB2607836A (en) Method for determining change in nanoscale pore structure, and use thereof
Xie et al. Nano-pore structure and fractal characteristics of shale gas reservoirs: a case study of Longmaxi Formation in southeastern Chongqing, China
CN106285652B (zh) 确定页岩游离气体饱和度的方法
Wei et al. Main factors influencing the development of nanopores in over-mature, organic-rich shales
Tao et al. Comparative analysis of soil-water characteristic curve in fractal and empirical models
Liu et al. Reservoir characteristics and methane adsorption capacity of the Upper Triassic continental shale in Western Sichuan Depression, China
Zhang et al. Mechanism of methane adsorption on groove space in organic matter surface
Li et al. Pore‐Throat Combination Types and Gas‐Water Relative Permeability Responses of Tight Gas Sandstone Reservoirs in the Zizhou Area of East Ordos Basin, China
CN111487176B (zh) 一种页岩油系统中液态烃所占孔隙度的计算方法
CN110632677B (zh) 一种碳酸盐岩地热储层分类方法
CN111664825B (zh) 一种提高孔结构参数测定准确性的方法
Ji et al. Investigation into the apparent permeability and gas-bearing property in typical organic pores in shale rocks
Tang et al. Investigation of shale nano-pore characteristics by scanning electron microscope and low-pressure nitrogen adsorption

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
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20200508

WD01 Invention patent application deemed withdrawn after publication