CN119827552B - 一种考虑有机碳含量的煤岩核磁共振孔隙度校正方法 - Google Patents
一种考虑有机碳含量的煤岩核磁共振孔隙度校正方法Info
- Publication number
- CN119827552B CN119827552B CN202510030184.4A CN202510030184A CN119827552B CN 119827552 B CN119827552 B CN 119827552B CN 202510030184 A CN202510030184 A CN 202510030184A CN 119827552 B CN119827552 B CN 119827552B
- Authority
- CN
- China
- Prior art keywords
- nuclear magnetic
- coal
- rock
- porosity
- magnetic resonance
- 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
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种考虑有机碳含量的煤岩核磁共振孔隙度校正方法,包括对煤岩样品进行饱和水处理,并获取煤岩样品的饱和水核磁孔隙度;选取相同煤岩样品开展总有机碳测定,获取煤岩样品的总有机碳含量;开展煤岩镜质体反射率测定,获取煤岩样品的镜质体反射率;根据饱和水核磁孔隙度、总有机碳含量、镜质体反射率进行煤岩核磁共振饱水孔隙度校正,得到校正孔隙度。本发明引入煤岩基础数据测试所得的总有机碳含量和代表煤岩演化阶段的镜质体反射率数据,进行校正以消除煤岩岩石基质中的核磁信号,校正后的数据更为真实可靠且方法方便实惠,可以大面积展开,对于煤岩气的大规模勘探与开发具有重要的指导意义。
Description
技术领域
本发明涉及一种考虑有机碳含量的煤岩核磁共振孔隙度校正方法,属于地质勘探技术领域。
背景技术
煤岩气作为天然气行业新兴的重要领域,煤岩孔隙度分析对于评价煤岩气的储集能力及煤岩气的勘探开发至关重要。孔隙度是指岩石中未被固态物质填满的空洞体积与岩石整体体积的比值。它是一个重要的地质学参数,因为它直接影响到岩石的渗透性和储存能力(SY/T 5336-2006岩心分析方法)。然而,由于煤岩岩石特征的特殊性,用真空干燥器中干燥的煤层岩心进行煤的孔隙度测定。煤的风化研究表明这种干燥方法改变了煤的结构,得出了错误的结果。用未干燥岩心测定层理孔隙度的混相驱替法中,示踪剂发生吸附就会产生错误的结果。
基于气体在岩石孔隙中的扩散和吸附行为,通过测量气体在岩石中的平衡压力来计算孔隙度,测量过程可能需要较长时间,主要用于常规储层的孔隙度测量,尤其是在孔隙结构简单的岩石中,主要测量连通孔隙度,难以直接测量不连通孔隙度,数据解释相对直观,但可能需要考虑气体饱和度和岩石类型的影响。
因此,为避免传统测试方法的实验误差,现引入核磁共振(NMR)技术进行孔隙测试,核磁共振(NMR)技术基于原子核自旋在外加磁场中的磁性响应,通过测量核磁共振信号的弛豫时间来计算孔隙度,是一种快速、可重复、无损的检测方法,能够区分不同类型的孔隙度,如总孔隙度、毛细孔隙度、有效孔隙度等,并能区分自由流体和束缚流体,广泛应用于岩心和储层(尤其是非常规油气藏,如煤岩气)的孔隙度测试。通过测量岩心样品中水分子的弛豫时间,可以间接得到孔隙结构的信息,从而计算出孔隙度。在测试孔隙度时,核磁共振的原理是对岩体孔隙结构中流体氢核1H的检测,由于岩体骨架中的碳元素在测量中不会产生核磁信号,从而不会影响测量结果。但在针对煤岩的应用过程中,核磁共振技术存在部分缺陷。煤岩由古代植物大量堆积在特定的物理化学及地质条件下形成,主要是有机物质组成,含有丰富的有机物质。煤岩的核磁共振孔隙度测试,特别是总有机碳含量高(普遍高于70%)的煤岩,有机质分子中的大量氢核也会被核磁共振测试出,这会导致测量所得的核磁孔隙度偏高,影响核磁结果的可靠性,如图2。杨文宽于1987年提出的煤分子式C102H96O26可以看出,煤岩含有与C含量相当的H元素,可以利用其总有机碳含量代表其H核含量。煤岩总有机碳含量主要受煤岩形成时的有机物质含量,以及煤岩形成后的埋藏演化控制,且总有机碳含量和代表煤岩演化程度的镜质体反射率可以快速实惠测得。
发明内容
为了克服现有技术中存在的缺陷,本发明旨在提供一种考虑有机碳含量的煤岩核磁共振孔隙度校正方法。
本发明解决上述技术问题所提供的技术方案是:一种考虑有机碳含量的煤岩核磁共振孔隙度校正方法,包括以下步骤:
S1、获取煤岩样品;
S2、对煤岩样品进行饱和水处理,并获取煤岩样品的饱和水核磁孔隙度φw;
S3、对煤岩样品进行烘干处理,并获取干燥状态下煤岩样品的岩石骨架氢核数据Rm;
S4、选取相同煤岩样品开展总有机碳测定,获取煤岩样品的总有机碳含量TOC;
S5、开展煤岩镜质体反射率测定,获取煤岩样品的镜质体反射率Ro;
S6、分别建立总有机碳与镜质体反射率和岩石骨架核磁数据的一元回归方程;
S7、建立总有机碳与镜质体反射率和煤岩岩石骨架核磁数据的二元回归方程;
S8、建立饱和水核磁孔隙度φw与总有机碳含量TOC、镜质体反射率Ro和校正后核磁共振孔隙度的计算公式;
φR=φw-2.9768ln(TOC)+0.0543ln(Ro)+10.2738
式中:φR为校正孔隙度;φw为饱和水核磁孔隙度;TOC为总有机碳含量;Ro为镜质体反射率;
S9、根据饱和水核磁孔隙度φw、总有机碳含量TOC、镜质体反射率Ro进行煤岩核磁共振饱水孔隙度校正,得到校正孔隙度φR。
进一步的技术方案是,所述步骤S1中煤岩样品为目标区域的煤岩岩心或剖面露头样品。
进一步的技术方案是,所述步骤S2中在真空条件下加压15MPa饱和水48h进行饱和水处理。
进一步的技术方案是,所述步骤S2中饱水后开展样品饱和水核磁共振分析测试,得到饱和水样品核磁共振信号和弛豫时间,并基于饱和水样品核磁共振信号计算得到煤岩样品的饱和水核磁孔隙度φw。
进一步的技术方案是,所述步骤S3中对煤岩样品在120℃条件下干燥24h进行烘干处理。
进一步的技术方案是,所述步骤S3中对不饱和流体进行核磁共振分析,得到干燥样品核磁共振信号和弛豫时间,并计算得出干燥状态下煤岩样品岩石骨架氢核数据。
进一步的技术方案是,所述步骤S6中的一元回归方程包括:
Rm=3.0164ln(TOC)-10.475
Rm=-4.006ln(Ro)+4.7922
式中:Rm为岩石骨架氢核数据;TOC为总有机碳含量;Ro为镜质体反射率。
进一步的技术方案是,所述步骤S7中的二元回归方程为:
Rm=2.9768ln(TOC)-0.0543ln(Ro)-10.2738
式中:Rm为岩石骨架氢核数据;TOC为总有机碳含量;Ro为镜质体反射率。
本发明具有以下有益效果:本发明针对目前煤岩核磁共振孔隙度测试存在的误差和弊端,引入煤岩基础数据测试所得的总有机碳含量和代表煤岩演化阶段的镜质体反射率数据,进行校正以消除煤岩岩石基质中的核磁信号,校正后的数据更为真实可靠且方法方便实惠,可以大面积展开,对于煤岩气的大规模勘探与开发具有重要的指导意义。
附图说明
图1为流程图;
图2为核磁孔隙测试概念图;
图3为不同地区样品饱水核磁共振T2谱图;
图4为不同地区样品干燥后核磁共振T2谱图;
图5为TOC与Rm(a)、Ro与Rm(b)线性回归图;
图6为二元回归线性拟合图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的一种深层页岩天然裂缝水化后力学强度的测试方法,具体包括以下步骤:
S1、采集样品,获取盆地不同地区不同深度段不同物理性质的煤岩岩心和剖面露头样品;
采集鄂尔多斯盆地扒楼沟剖面二叠、神木地区石炭系井下岩心、佳县地区石炭系井下岩心和横山地区二叠系与石炭系井下岩心共15块煤岩样品,对其开展总有机碳和镜质体反射率测定。
表1煤岩样品分析测试结果表
S2、对所有煤岩样品在真空条件下加压15MPa饱和水48h进行饱和水处理,饱水后开展样品饱和水核磁共振分析测试,得到饱和水样品核磁共振信号和弛豫时间;并基于样品饱水核磁共振信号计算得出所有样品的饱和水核磁孔隙度φw;
S3、对煤岩样品在120℃条件下干燥24h进行烘干处理,然后不饱和流体进行核磁共振分析,得到干燥样品核磁共振信号和弛豫时间,并计算得出干燥状态下煤岩样品的岩石骨架氢核数据Rm;
S4、选取相同煤岩样品开展总有机碳测定,获取煤岩样品的总有机碳含量TOC;
S5、开展煤岩镜质体反射率测定,获取煤岩样品的镜质体反射率Ro;
S6、分别建立总有机碳与镜质体反射率Ro和岩石骨架核磁数据Rm的一元回归方程;
其中总有机碳含量TOC与煤岩样品的岩石骨架核磁数据Rm的一元回归方程为:
Rm=3.0164ln(TOC)-10.475,(R2=0.91)
式中:Rm为岩石骨架氢核数据;TOC为总有机碳含量;
镜质体反射率Ro与煤岩岩石骨架核磁数据Rm的一元回归方程为:
Rm=-4.006ln(Ro)+4.7922,(R2=0.92)
式中:Rm为岩石骨架氢核数据;Ro为镜质体反射率;
S7、建立总有机碳与镜质体反射率和煤岩岩石骨架核磁数据的二元回归方程;
y=-39.1202x1+174.5889x2-84.0028
即,
Rm=2.9768ln(TOC)-0.0543ln(Ro)-10.2738,(R2=0.94)
式中:Rm为岩石骨架氢核数据;TOC为总有机碳含量;Ro为镜质体反射率;
S8、建立饱和水核磁孔隙度φw与总有机碳含量TOC、镜质体反射率Ro和校正后核磁共振孔隙度的计算公式;
φR=φw-2.9768ln(TOC)+0.0543ln(Ro)+10.2738
式中:φR为校正孔隙度;φw为饱和水核磁孔隙度;TOC为总有机碳含量;Ro为镜质体反射率;
S9、根据饱和水核磁孔隙度φw、总有机碳含量TOC、镜质体反射率Ro进行煤岩核磁共振饱水孔隙度校正,得到校正孔隙度φR。
实施例
采集三个煤岩样品,开展饱水核磁共振测试,并利用本发明的校正公式进行校正,计算得到校正后煤岩孔隙度数据,结果如下。
样本1实测孔隙度3.12%,常规饱水核磁测试孔隙度5.54%,校正后孔隙度3.10%,校正前误差43.68%,校正后误差0.854%;
样本2实测孔隙度2.47%,常规饱水核磁测试孔隙度4.43%,校正后孔隙度2.49%,校正前误差44.24%,校正后误差0.977%;
样本3实测孔隙度2.46%,常规饱水核磁测试孔隙度4.59%,校正后孔隙度2.44%,校正前误差46.41%,校正后误差0.677%。
通过本发明进行核磁共振孔隙度校正结果与实测数据误差均小于1%。本发明技术准确可靠,计算效果好。
以上所述,并非对本发明作任何形式上的限制,虽然本发明已通过上述实施例揭示,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,可利用上述揭示的技术内容作出一些变动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均属于本发明技术方案的范围内。
Claims (8)
1.一种考虑有机碳含量的煤岩核磁共振孔隙度校正方法,其特征在于,包括以下步骤:
S1、获取煤岩样品;
S2、对煤岩样品进行饱和水处理,并获取煤岩样品的饱和水核磁孔隙度φw;
S3、对煤岩样品进行烘干处理,并获取干燥状态下煤岩样品的岩石骨架氢核数据Rm;
S4、选取相同煤岩样品开展总有机碳测定,获取煤岩样品的总有机碳含量TOC;
S5、开展煤岩镜质体反射率测定,获取煤岩样品的镜质体反射率Ro;
S6、分别建立总有机碳与镜质体反射率和岩石骨架核磁数据的一元回归方程;
S7、建立总有机碳与镜质体反射率和煤岩岩石骨架核磁数据的二元回归方程;
S8、建立饱和水核磁孔隙度φw与总有机碳含量TOC、镜质体反射率Ro和校正后核磁共振孔隙度的计算公式;
φR=φw-2.9768ln(TOC×100)+0.0543ln(Ro×100)+10.2738
式中:φR为校正孔隙度;φw为饱和水核磁孔隙度;TOC为总有机碳含量;Ro为镜质体反射率;
S9、根据饱和水核磁孔隙度φw、总有机碳含量TOC、镜质体反射率Ro进行煤岩核磁共振饱水孔隙度校正,得到校正孔隙度φR。
2.根据权利要求1所述的一种考虑有机碳含量的煤岩核磁共振孔隙度校正方法,其特征在于,所述步骤S1中煤岩样品为目标区域的煤岩岩心或剖面露头样品。
3.根据权利要求1所述的一种考虑有机碳含量的煤岩核磁共振孔隙度校正方法,其特征在于,所述步骤S2中在真空条件下加压15MPa饱和水48h进行饱和水处理。
4.根据权利要求3所述的一种考虑有机碳含量的煤岩核磁共振孔隙度校正方法,其特征在于,所述步骤S2中饱水后开展样品饱和水核磁共振分析测试,得到饱和水样品核磁共振信号和弛豫时间,并基于饱和水样品核磁共振信号计算得到煤岩样品的饱和水核磁孔隙度φw。
5.根据权利要求1所述的一种考虑有机碳含量的煤岩核磁共振孔隙度校正方法,其特征在于,所述步骤S3中对煤岩样品在120℃条件下干燥24h进行烘干处理。
6.根据权利要求5所述的一种考虑有机碳含量的煤岩核磁共振孔隙度校正方法,其特征在于,所述步骤S3中对不饱和流体进行核磁共振分析,得到干燥样品核磁共振信号和弛豫时间,并计算得出干燥状态下煤岩样品岩石骨架氢核数据。
7.根据权利要求1所述的一种考虑有机碳含量的煤岩核磁共振孔隙度校正方法,其特征在于,所述步骤S6中的一元回归方程包括:
Rm=3.0164ln(TOC×100)-10.475
Rm=-4.006ln(Ro×100)+4.7922
式中:Rm为岩石骨架氢核数据;TOC为总有机碳含量;Ro为镜质体反射率。
8.根据权利要求1所述的一种考虑有机碳含量的煤岩核磁共振孔隙度校正方法,其特征在于,所述步骤S7中的二元回归方程为:
Rm=2.9768ln(TOC×100)-0.0543ln(Ro×100)-10.2738
式中:Rm为岩石骨架氢核数据;TOC为总有机碳含量;Ro为镜质体反射率。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202510030184.4A CN119827552B (zh) | 2025-01-08 | 2025-01-08 | 一种考虑有机碳含量的煤岩核磁共振孔隙度校正方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202510030184.4A CN119827552B (zh) | 2025-01-08 | 2025-01-08 | 一种考虑有机碳含量的煤岩核磁共振孔隙度校正方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN119827552A CN119827552A (zh) | 2025-04-15 |
| CN119827552B true CN119827552B (zh) | 2025-07-15 |
Family
ID=95291958
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202510030184.4A Active CN119827552B (zh) | 2025-01-08 | 2025-01-08 | 一种考虑有机碳含量的煤岩核磁共振孔隙度校正方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN119827552B (zh) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115728337A (zh) * | 2021-08-30 | 2023-03-03 | 中国石油天然气集团有限公司 | 差值核磁共振法测量页岩气储层孔隙度的方法 |
| CN116413184A (zh) * | 2021-12-31 | 2023-07-11 | 中国石油天然气股份有限公司 | 一种快速评价页岩总孔隙度、含烃孔隙度和渗透率的方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116415699A (zh) * | 2021-12-24 | 2023-07-11 | 中国石油天然气股份有限公司 | 一种计算页岩初始有机碳含量的方法 |
| CN118008277A (zh) * | 2022-11-08 | 2024-05-10 | 中国石油天然气集团有限公司 | 一种泥页岩储层的可动油含量评价方法、系统、设备与存储介质 |
-
2025
- 2025-01-08 CN CN202510030184.4A patent/CN119827552B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115728337A (zh) * | 2021-08-30 | 2023-03-03 | 中国石油天然气集团有限公司 | 差值核磁共振法测量页岩气储层孔隙度的方法 |
| CN116413184A (zh) * | 2021-12-31 | 2023-07-11 | 中国石油天然气股份有限公司 | 一种快速评价页岩总孔隙度、含烃孔隙度和渗透率的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN119827552A (zh) | 2025-04-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109060623B (zh) | 致密砂岩孔隙大小分布联合表征方法 | |
| CN105910971B (zh) | 富有机质致密岩心气体渗透率和扩散系数的联测方法 | |
| US5265462A (en) | Method and apparatus for determining permeability, diffusivity, porosity, and gas storage in gas-containing substrates | |
| CN104502224B (zh) | 饱和水条件下煤岩等温解吸曲线测定装置及方法 | |
| CN110501272B (zh) | 三轴应力和孔隙压力条件下同时测试多孔岩石孔隙率和渗透率的方法 | |
| CN110261274B (zh) | 自发渗吸作用对水驱驱油效率静态贡献率的评价方法 | |
| CN110296931B (zh) | 一种致密砂岩油水相对渗透率信息的表征方法及系统 | |
| CN102879306A (zh) | 一种利用磁共振成像技术检测气液扩散过程的装置及其方法 | |
| CN110160933B (zh) | 定量评价致密砂岩油藏自发渗吸驱油速度的方法 | |
| CN116660306B (zh) | 一种基于柱塞样和碎样测量页岩含气量的实验方法 | |
| CN115420766B (zh) | 基于t2谱分析的砂岩注水采油剩余油分布测定方法 | |
| WO2025087342A1 (zh) | 深层煤岩气储层含气量确定方法、装置及可读存储介质 | |
| CN115876659B (zh) | 一种获取页岩孔隙和裂缝体积的方法及系统 | |
| CN115755195B (zh) | 一种储层岩石模拟生产条件下的润湿性表征方法 | |
| CN111426617B (zh) | 一种地质条件下陆相湖盆致密储层孔隙度的定量预测方法 | |
| CN119827552B (zh) | 一种考虑有机碳含量的煤岩核磁共振孔隙度校正方法 | |
| CN110672495A (zh) | 一种基于低场磁共振技术的水泥基材料水分渗透率的预测方法 | |
| CN117630079B (zh) | 一种煤层气吸附态和游离态含量的动态评估方法及系统 | |
| CN119395070A (zh) | 一种页岩油储层含油量的评价方法 | |
| CN119780142B (zh) | 多孔介质内混溶溶液中水的核磁信号分离方法 | |
| CN119064509A (zh) | 一种页岩多尺度孔隙中二氧化碳吞吐动用凝析气的流体动态含量量化实验装置及方法 | |
| CN114894671A (zh) | 基于扩散-渗流解耦实验的煤体甲烷吸附时间测定方法 | |
| CN116773578B (zh) | 对煤储层原位含水饱和度进行测试与计算的方法 | |
| CN116223213B (zh) | 一种基于声波响应的地层应力敏感性评价方法 | |
| CN119936096B (zh) | 一种页岩油储层原位孔隙油水分布核磁共振测试方法 |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |