CN111122408A - 一种孔径介于0.3-10nm的页岩有机孔隙度估算方法 - Google Patents
一种孔径介于0.3-10nm的页岩有机孔隙度估算方法 Download PDFInfo
- 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
Links
- 239000011148 porous material Substances 0.000 title claims abstract description 188
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 238000012360 testing method Methods 0.000 claims abstract description 12
- 238000003775 Density Functional Theory Methods 0.000 claims abstract description 11
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims abstract description 10
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims abstract description 10
- 238000001179 sorption measurement Methods 0.000 claims abstract description 9
- 229910002092 carbon dioxide Inorganic materials 0.000 claims abstract description 5
- 239000001569 carbon dioxide Substances 0.000 claims abstract description 5
- 229910052757 nitrogen Inorganic materials 0.000 claims abstract description 5
- 239000007789 gas Substances 0.000 claims description 44
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000011156 evaluation Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 2
- 229910052753 mercury Inorganic materials 0.000 description 2
- 239000005416 organic matter Substances 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000000197 pyrolysis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
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
本发明公开了一种孔径介于0.3‑10nm的页岩有机孔隙度估算方法,包括以下步骤:通过页岩样品二氧化碳和氮气吸附测试,利用密度函数理论DFT模型计算页岩样品内孔径为0.3‑10nm的孔隙体积;建立页岩样品孔径介于0.3‑10nm的孔隙体积与TOC含量的线性正相关关系;估算0.3‑10nm有机孔隙体积以及占总孔体积的比值,结合页岩总孔隙度估算页岩孔径介于0.3‑10nm的页岩有机孔隙度,该方法思路清晰、技术可靠、涉及的参数物理意义明确、估算方法和计算公式应用方便。
Description
技术领域
本发明涉及页岩气储层研究及页岩气勘探,具体涉及一种孔径介于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;
进一步地,步骤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的页岩孔隙体积。
本发明提供的技术方案带来的有益效果是:解决页岩气储层中孔径介于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孔隙体积和有机孔体积计算表
步骤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内的页岩孔隙体积。
由于二氧化碳和和氮气吸附实验的限制,页岩气储层孔隙总体积可能不完全,可结合压汞实验获得页岩气储层全孔径的孔隙体积,本实施例中的页岩气储层孔隙总体积是由二氧化碳和氮气吸附实验计算获得,此类计算孔径范围限于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;
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的页岩孔隙体积。
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)
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)
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 | 中国石油化工股份有限公司 | 一种确定泥页岩不同类型孔隙度的方法 |
-
2019
- 2019-12-23 CN CN201911335718.5A patent/CN111122408A/zh active Pending
Patent Citations (7)
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)
Title |
---|
李军 等: "页岩气储层"四孔隙"模型建立及测井定量表征方法", 《石油与天然气地质》 * |
陈曼霏等: "涪陵页岩气田平桥区块页岩气储层有机质孔发育特征", 《石油学报》 * |
Cited By (2)
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 |