CN106501149B - 一种基于煤岩工业成分分析的孔隙度计算方法 - Google Patents
一种基于煤岩工业成分分析的孔隙度计算方法 Download PDFInfo
- Publication number
- CN106501149B CN106501149B CN201610930124.9A CN201610930124A CN106501149B CN 106501149 B CN106501149 B CN 106501149B CN 201610930124 A CN201610930124 A CN 201610930124A CN 106501149 B CN106501149 B CN 106501149B
- Authority
- CN
- China
- Prior art keywords
- density
- coal petrography
- porosity
- coal
- mineral
- 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.)
- Expired - Fee Related
Links
- 239000003245 coal Substances 0.000 title claims abstract description 85
- 238000004458 analytical method Methods 0.000 title claims abstract description 20
- 238000004364 calculation method Methods 0.000 title claims abstract description 18
- 229910052500 inorganic mineral Inorganic materials 0.000 claims abstract description 38
- 239000011707 mineral Substances 0.000 claims abstract description 38
- 239000005416 organic matter Substances 0.000 claims abstract description 33
- 238000000034 method Methods 0.000 claims abstract description 17
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims abstract description 15
- 229910052799 carbon Inorganic materials 0.000 claims abstract description 15
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 claims abstract description 11
- 229910052717 sulfur Inorganic materials 0.000 claims abstract description 11
- 239000011593 sulfur Substances 0.000 claims abstract description 11
- 239000000203 mixture Substances 0.000 claims abstract description 6
- 239000011435 rock Substances 0.000 claims description 12
- 238000002441 X-ray diffraction Methods 0.000 claims description 5
- 238000012360 testing method Methods 0.000 claims description 5
- 238000013480 data collection Methods 0.000 claims description 4
- 239000000463 material Substances 0.000 claims description 2
- 238000011160 research Methods 0.000 abstract description 9
- 239000000470 constituent Substances 0.000 abstract description 3
- 239000000523 sample Substances 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 238000005259 measurement 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
- 239000011148 porous material Substances 0.000 description 2
- 239000010453 quartz Substances 0.000 description 2
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- PXGOKWXKJXAPGV-UHFFFAOYSA-N Fluorine Chemical compound FF PXGOKWXKJXAPGV-UHFFFAOYSA-N 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 229910052570 clay Inorganic materials 0.000 description 1
- 239000002734 clay mineral Substances 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000010433 feldspar Substances 0.000 description 1
- 239000011737 fluorine Substances 0.000 description 1
- 229910052731 fluorine Inorganic materials 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 239000011028 pyrite Substances 0.000 description 1
- 229910052683 pyrite Inorganic materials 0.000 description 1
- NIFIFKQPDTWWGU-UHFFFAOYSA-N pyrite Chemical compound [Fe+2].[S-][S-] NIFIFKQPDTWWGU-UHFFFAOYSA-N 0.000 description 1
- 239000000243 solution Substances 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
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)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明涉及一种计算煤岩孔隙度的方法。原则上所计算的孔隙度要求研究区具有煤质分析数据,包括煤的工业成份(灰分、挥发分、固定碳,含硫量)、煤岩视密度。其主要计算步骤为:获取煤样的体密度,获取煤样的工业成份数据(灰分、挥发分、固定碳,含硫量),结合已有的有机质密度和矿物密度数据,通过本发明专利的计算公式进行孔隙度计算。利用本发明的孔隙度计算方法,可以简单、大量的计算煤岩孔隙度,对地震岩石物理的研究有很大帮助。
Description
所属领域
本实用发明涉及一种基于煤岩工业成分分析的孔隙度计算方法。
技术背景
岩石物理学是架起孔隙介质的储层特性和地震特性之间的桥梁。在岩石物理分析中,根据储层特性建立精确的模型有助于提高正演的准确性,从而更好的模拟介质的地震特性以及提高反演解释的精度。储层特性中孔隙度、岩性、流体类型是建立岩石物理模型的基础,其中孔隙度是一个重要的参数。在以往的研究中,孔隙度通常采用压汞法和波尔法测定,不过这两种方法实验过程复杂,压汞法在测量过程中由于煤岩极软,容易破坏内部结构而造成结果偏差,波尔法在测量过程中由于煤本身具有吸附的性质,会造成结果偏大。
在煤炭领域,作为煤质关键指标的工业成分分析、体密度等数据则大量存在。在目前还未见到通过煤工业成分数据求取孔隙度的方法,这里建立一种快速计算煤岩孔隙的方法,利用工业成分获取孔隙度参数,即方便快捷又能提高岩石物理正演模拟的精度。
发明内容
本技术方法主要目的是计算煤岩岩石中的孔隙度,利用研究区煤质分析成果报告获取的煤岩视密度、煤工业成份(灰分、挥发分、固定碳、硫含量),通过本发明计算公式获得煤岩孔隙度。它给出了快速计算煤岩孔隙度的方法,具有简单、快速的特点。
本发明解决其技术问题所采取的技术方案:
(1)煤工业分析数据获取:通过资料收集获或者实验测试的方式取研究区内煤工业分析数据。
(2)有机质密度:参考已有煤岩中的有机质密度参数。
(3)矿物密度计算:通过本发明中的公式计算矿物平均密度,其中需要用到的某种矿物质量分数可以根据X-ray或者电子探针分析获得,也可以根据研究区内已有的资料获得。
(4)孔隙度计算:通过本发明中的公式计算煤岩孔隙度。
优选地,所述煤岩孔隙度计算方法具体包括以下步骤:
(1)煤工业分析数据获取:收集研究区内已有的煤质分析报告,或者根据国标GB/T212-2008实验测试获取煤工业分析数据,为煤岩孔隙度的计算提供必要的参数:空气干燥基灰分(Aad)——灰分质量分数;空气干燥基挥发分(Vad)——挥发分质量分数、空气干燥基固定碳(FCad)——固定碳质量分数;黄铁矿硫含量(Sp,d)——黄铁矿硫质量分数。
(2)煤岩体密度:根据国标GB/T6949-2010煤的相对视密度测定方法确定,或者根据研究区内已有资料获得。
(3)有机质密度:参考已有煤岩中的有机质密度参数。
(4)矿物密度确定:通过X-ray分析或者电子探针,或者已有资料,获得研究区煤岩矿物质量含量数据,通过本发明中的公式计算矿物平均密度。
(5)孔隙度计算:输入上述步骤中得到的参数,通过本发明公式计算煤岩孔隙度。
附图说明
下面结合附图和实施例对本发明进一步说明。
图1是孔隙度计算流程图
图2是煤岩组成模型
具体实施方式
以下将结合附图对本发明进行详细说明,但本发明并不限于此。
根据大量的岩石样品成分分析,岩石的组成主要包括矿物(碳酸盐岩矿物、粘土、石英、黄铁矿等)、孔隙流体(水、油或气)和有机质。根据组成成分,建立简单的岩石模型,如图2所示。
煤工业成分分析中的四大组分为水分、灰分、挥发分和固定碳。其中水分是吸附在煤孔隙表面的内在水,属于孔隙空间的一部分;灰分是是赋存于煤中的矿物质在815o±10℃条件下完全燃烧后所得残留物但并不完全等同于矿物,可以利用Parr校正公式进行校正从而表征矿物含量;挥发分是煤中有机质高温隔绝空气状态下热解的产物;固定碳是煤中有机质的主要部分。因此,利用工业成分分析中的灰分表征矿物含量,挥发分和固定碳表征有机质含量就可以计算出煤岩的孔隙度。
根据图2所示,孔隙度、矿物体积含量以及有机质含量应该满足如下关系:
其中:
Φ-孔隙度,%;
F矿-全部矿物体积分量,%;
F有机质-煤岩有机质体积分量,%;
根据如上公式,要确定煤岩孔隙度,关键是求取煤岩中矿物以及有机质的体积分量。下面分以下三个步骤实现煤岩孔隙度的计算。
①煤岩中矿物体积含量的确定:
式中
F矿=矿物占煤岩体积分量,%;
WM=煤中矿物占煤基质的质量分数,%;
ρ矿=全部矿物平均体密度,g/cm3;
ρ体=煤岩的体密度,g/cm3;
全部矿物平均密度ρ矿,通过公式(3)计算:
其中Wi和ρi分别为岩石各组成矿物的质量分数和体密度。通过电子探针或者X-ray分析、或者已有资料,获得各种矿物含量数据。其中各种矿物密度值可以参考已有的测试成果,比如岩石物理手册。
WM由Parr公式(1975年)给定
WM=1.08Aad+0.55Sp,d (4)
式中
Aad=灰分质量分量,%;
Sp,d=煤岩中黄铁矿硫质量分数,%;
②对于岩石的有机质,通过煤的工业成份分析,可以用挥发分与固定碳之和来代替。
式中
F有机质=有机质占煤体积分量,%;
Vad=挥发分质量分数,%;
FCad=固定碳质量分数,%;
ρ体=煤岩体密度,g/cm3;
ρ有机质=有机质密度,g/cm3;
③煤岩孔隙度求取
式中
φ预测=煤岩预测孔隙度值,%;
Vad=挥发分质量分数,%;
Aad=灰分质量分数,%;
FCad=固定碳质量分数,%;
Sp,d=煤岩中黄铁矿硫质量分数,%;
ρ体=煤岩体密度,g/cm3;
ρ有机质=有机质密度,g/cm3;
公式(6)即为煤岩孔隙度的计算公式。从中可以看出,孔隙度Φ计算时的输入参数是矿物平均体密度ρ矿、煤岩体密度ρ体、挥发分质量分数、灰分质量分数Aad、固定碳质量分数FCad、黄铁矿硫质量分数Sp,d。
实施例1:山西晋城赵庄矿煤岩孔隙度计算,于2016年7月15日进行。
利用公式(6)计算煤岩孔隙度需要的输入参数是体密度(ρ体)、矿物平均密度(ρ矿)、有机质平均密度(ρ有机质)、空气干燥基灰分(Aad)、空气干燥基挥发分(Vad)、空气干燥基固定碳(FCad)、黄铁矿硫含量(Sp,d)。收集研究区内煤质分析报告,其中ρ体、Aad、Vad、FCad、Sp,d可以直接从煤质分析报告中获取。通过查阅已有文献资料,得知有机质密度为1.37g/cc。通过X-ray衍射得到工区煤岩矿物含量为:粘土矿物占8.59%、石英占11.4%、长石占1.7%、黄铁矿为1%,根据此结果,由本专利公式(3)计算得到矿物密度ρm为2.7g/cc。根据本专利中的公式(6)计算获得研究区内煤岩孔隙度,结果如表1所示。
表1研究区煤岩孔隙度计算表
注:表中密度单位均为(g/cc)。
Claims (5)
1.一种基于煤岩工业成分分析的孔隙度计算方法,包括四个步骤:(1)煤工业分析数据获取,(2)确定有机质密度,(3)计算矿物密度,(4)计算孔隙度;其特征是:(1)通过收集资料获取工区内煤工业分析数据,(2)通过已有文献确定有机质密度,(3)计算矿物密度,(4)根据所述有机质密度和所述矿物密度计算煤岩孔隙度;
步骤(4)中,通过公式(2),计算获得煤岩孔隙度φ预测;
其中,φ预测是煤岩预测孔隙度值;
F矿是矿物占煤岩体积分量;
F有机质是煤岩有机质体积分量;
Vad是挥发分质量分数;
Aad是灰分质量分数;
FCad是固定碳质量分数;
Sp,d是煤岩中黄铁矿硫质量分数;
ρ体是煤岩体密度;
ρ矿是矿物平均密度;
ρ有机质是有机质密度。
2.根据权利要求1所述的方法,其特征是,(1)通过资料收集获取煤岩工业成分数据,(2)通过已有文献资料确定煤岩中有机质密度,(3)结合X-ray衍射结果,利用矿物密度计算矿物平均密度,(4)根据所述有机质密度和所述矿物平均密度计算孔隙度。
3.根据权利要求2中所述的方法,其特征在于,步骤(1)中,通过工区资料收集获取煤岩工业成分数据,为煤岩孔隙度计算提供5个输入参数,即空气干燥基灰分Aad、空气干燥基挥发分Vad、空气干燥基固定碳FCad、黄铁矿硫含量Sp,d、ρ体。
4.根据权利要求2中所述的方法,其特征在于,步骤(2)中,通过已有测试成果获取有机质密度,为煤岩孔隙度计算提供一个输入参数ρ有机质。
5.根据权利要求2中所述的方法,其特征在于,步骤(3)中,利用X-ray射线衍射得到工区矿物含量,利用公式(1)计算矿物平均密度,为煤岩孔隙度计算提供一个输入参数ρ矿;其中Wi和ρi分别为岩石各组成矿物的质量分数和体密度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610930124.9A CN106501149B (zh) | 2016-10-31 | 2016-10-31 | 一种基于煤岩工业成分分析的孔隙度计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610930124.9A CN106501149B (zh) | 2016-10-31 | 2016-10-31 | 一种基于煤岩工业成分分析的孔隙度计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106501149A CN106501149A (zh) | 2017-03-15 |
CN106501149B true CN106501149B (zh) | 2019-08-30 |
Family
ID=58319645
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610930124.9A Expired - Fee Related CN106501149B (zh) | 2016-10-31 | 2016-10-31 | 一种基于煤岩工业成分分析的孔隙度计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106501149B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108536994A (zh) * | 2018-02-13 | 2018-09-14 | 武汉科技大学 | 一种煤基因库的建立方法及其应用方法 |
CN108412488B (zh) * | 2018-03-13 | 2021-09-10 | 中石化石油工程技术服务有限公司 | 快速确定页岩气储层有机孔隙度的测井方法 |
CN110715888B (zh) * | 2019-05-30 | 2022-04-22 | 西南石油大学 | 一种高含硫气藏岩心硫沉积孔隙尺寸分布变化的测量方法 |
CN113533160B (zh) * | 2021-07-06 | 2023-03-24 | 长江大学 | 一种岩石孔隙度测量方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013148632A1 (en) * | 2012-03-29 | 2013-10-03 | Ingrain, Inc. | A method and system for estimating properties of porous media such as fine pore or tight rocks |
CN103790581A (zh) * | 2014-01-16 | 2014-05-14 | 燕山大学 | 一种确定煤岩工业组分物理参数的方法 |
CN104806232A (zh) * | 2014-07-31 | 2015-07-29 | 吴见萌 | 一种确定碳酸盐岩储层孔隙度下限的方法 |
CN105628576A (zh) * | 2014-10-31 | 2016-06-01 | 中国石油天然气股份有限公司 | 灰质白云岩孔隙度测量方法及装置 |
CN106019403A (zh) * | 2016-06-08 | 2016-10-12 | 西北大学 | 自生自储油气层孔隙度测量方法 |
-
2016
- 2016-10-31 CN CN201610930124.9A patent/CN106501149B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013148632A1 (en) * | 2012-03-29 | 2013-10-03 | Ingrain, Inc. | A method and system for estimating properties of porous media such as fine pore or tight rocks |
CN103790581A (zh) * | 2014-01-16 | 2014-05-14 | 燕山大学 | 一种确定煤岩工业组分物理参数的方法 |
CN104806232A (zh) * | 2014-07-31 | 2015-07-29 | 吴见萌 | 一种确定碳酸盐岩储层孔隙度下限的方法 |
CN105628576A (zh) * | 2014-10-31 | 2016-06-01 | 中国石油天然气股份有限公司 | 灰质白云岩孔隙度测量方法及装置 |
CN106019403A (zh) * | 2016-06-08 | 2016-10-12 | 西北大学 | 自生自储油气层孔隙度测量方法 |
Non-Patent Citations (2)
Title |
---|
吐哈盆地巴喀气田八道湾组致密砂岩储层分析及孔隙度演化定量模拟;王国亭 等;《地质学报》;20121231;第86卷(第11期);第1847-1854页 * |
砂砾岩储层孔隙度和渗透率预测方法;张丽艳;《测井技术》;20050630;第29卷(第3期);第212-215页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106501149A (zh) | 2017-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106501149B (zh) | 一种基于煤岩工业成分分析的孔隙度计算方法 | |
CN103792338B (zh) | 一种烃源岩有机碳含量的确定方法和装置 | |
Li et al. | Experimental research on coal permeability: the roles of effective stress and gas slippage | |
El Nady et al. | Evaluation of organic matters, hydrocarbon potential and thermal maturity of source rocks based on geochemical and statistical methods: Case study of source rocks in Ras Gharib oilfield, central Gulf of Suez, Egypt | |
Cai et al. | Petrophysical characterization of Chinese coal cores with heat treatment by nuclear magnetic resonance | |
WO2017024530A1 (zh) | 一种计算烃源岩中有机碳含量的方法 | |
CN104278991B (zh) | 盐湖相烃源岩有机碳和生烃潜量的多元测井计算方法 | |
Lecompte et al. | Quantifying source rock maturity from logs: how to get more than TOC from delta log R | |
CN103344541A (zh) | 一种泥页岩总孔隙度测定方法 | |
Zhou et al. | Evolution of coal permeability with cleat deformation and variable Klinkenberg effect | |
CN104656162A (zh) | 一种富铀烃源岩铀含量及有机碳含量确定方法 | |
CN104007482A (zh) | 一种基于各向异性有效场的泥页岩岩石物理模型方法 | |
Xu et al. | Brittleness and rock strength of the Bakken formation, Williston basin, North Dakota | |
Yang et al. | Desorption-diffusion model and lost gas quantity estimation of coalbed methane from coal core under drilling fluid medium | |
Deng et al. | Characterization of mineral composition and its influence on microstructure and sorption capacity of coal | |
Zhong et al. | Experimental study of the impact on methane adsorption capacity of continental shales with thermal evolution | |
CN110909487A (zh) | 一种精细各向异性页岩岩石物理模型的建立方法 | |
CN106054279A (zh) | 一种煤岩脆性指数的确定方法 | |
Wang et al. | Dynamic-change laws of the porosity and permeability of low-to medium-rank coals under heating and pressurization treatments in the eastern junggar basin, China | |
Yao et al. | Methods for petrological and petrophysical characterization of gas shales | |
Wang et al. | Adsorption effect of overlying strata on carbon dioxide in coalfield fire area | |
Cao et al. | A new gas-content-evaluation method for organic-rich shale using the fractionation of carbon isotopes of methane | |
Chen et al. | Shale gas reservoir characterization: A typical case in the Southeast Chongqing of Sichuan Basin, China | |
Zhu et al. | New parameters for charactering the gas-bearing properties of shale gas | |
Wu et al. | Uncertainty analysis of shale gas simulation: consideration of basic petrophysical properties |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190830 Termination date: 20211031 |