CN111855537A - 一种基于hrtem的煤中微小孔径的测定方法 - Google Patents
一种基于hrtem的煤中微小孔径的测定方法 Download PDFInfo
- Publication number
- CN111855537A CN111855537A CN202010852224.0A CN202010852224A CN111855537A CN 111855537 A CN111855537 A CN 111855537A CN 202010852224 A CN202010852224 A CN 202010852224A CN 111855537 A CN111855537 A CN 111855537A
- Authority
- CN
- China
- Prior art keywords
- coal
- image
- hrtem
- pores
- block
- 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
- 239000003245 coal Substances 0.000 title claims abstract description 49
- 239000011148 porous material Substances 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 title claims abstract description 18
- 238000002173 high-resolution transmission electron microscopy Methods 0.000 title claims abstract description 15
- 238000004458 analytical method Methods 0.000 claims abstract description 7
- 238000005115 demineralization Methods 0.000 claims abstract description 7
- 230000002328 demineralizing effect Effects 0.000 claims abstract description 7
- 230000005540 biological transmission Effects 0.000 claims abstract description 4
- 230000011218 segmentation Effects 0.000 claims abstract description 4
- 229920002521 macromolecule Polymers 0.000 abstract description 4
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 239000007789 gas Substances 0.000 description 3
- 238000001179 sorption measurement Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 2
- 239000007788 liquid Substances 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
- 229910052757 nitrogen Inorganic materials 0.000 description 2
- 238000000235 small-angle X-ray scattering Methods 0.000 description 2
- 241001436672 Bhatia Species 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000003795 desorption Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000001000 micrograph Methods 0.000 description 1
- 238000012900 molecular simulation Methods 0.000 description 1
- 238000001956 neutron scattering Methods 0.000 description 1
- 238000003860 storage 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10056—Microscopic image
- G06T2207/10061—Microscopic image from scanning electron microscope
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20056—Discrete and fast Fourier transform, [DFT, FFT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20112—Image segmentation details
- G06T2207/20152—Watershed segmentation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Dispersion Chemistry (AREA)
- Geometry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Length-Measuring Devices Using Wave Or Particle Radiation (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明涉及煤孔隙结构测定技术领域,具体涉及一种基于HRTEM的煤中微小孔径的测定方法;包括以下步骤:1),对煤样进行脱矿物质处理,然后进行高分辨透射电镜分析,对原始图像进行去噪,傅里叶‑反傅里叶变换,图像的二值化处理;2)利用ImageJ软件进行图像反转(Invert);3),对反转后的图像进行分水岭分割算法;4),对步骤3)得到的图像计算出每个分割后区块的面积,即为单个微小孔的面积;5),将得到的每个区块近似拟合为正方形,计算其边长,即可得到孔径分布;本发明方法较为直观的再现了煤大分子基本结构单元中孔的分布情况,尤其是对于微小孔的具***置,对于研究煤中微小孔隙分布和煤大分子结构基本单元之间的关系起到了良好的指导作用。
Description
技术领域
本发明涉及煤孔隙结构测定技术领域,具体涉及一种基于HRTEM的煤中微小孔径的测定方法。
背景技术
煤中微小孔已被证明是煤层气储存的主要空间,也是煤层气运输的重要隧道,有充分证据表明,煤中微小孔提供了大部分的总孔体积和总比表面积,直接影响了游离气体储存、吸收、解吸和扩散。因此,煤中微小孔已成为当下研究的热点,诸多学者展开了研究,并取得显著成果,煤化学结构的变化对煤中的纳米级孔有显着影响,特别是对于微孔(<2nm)。由于实验条件的限制,微小孔很难直观的体现,因此,只能通过构建煤的3D分子结构模型来显示不同煤中微小孔的分布情况。关于这项工作,最早的学者Faulon等使用三维(3D)煤分子结构来研究微孔尺寸分布,证实了在3D煤分子结构内存在<1nm的超微孔,其超微孔体积与低压CO2吸附数据得到的一致。在随后的研究中,纳米孔隙尺度渐渐变小,已经达到0.6nm的水平,Feng和Bhatia认为小微孔(约0.6nm)可以代表两个微晶的基面之间的区域,较大的微孔(1.3nm)应该代表两个微晶的边缘位置之间的空间。
此前,对于微小孔的研究大部分采用的分子模拟的办法,却未能更真实的反应煤中孔隙分布。近年来,多种方法应运而生,诸如压汞法、液氮吸附法、小角X-射线散射(SAXS)分析、中子散射(SANS)分析和扫描电镜(SEM)分析等。由于煤是由煤大分子基本结构单元(BSU)构成,很难通过普通的技术手段量化纳米级孔,例如,压汞实验主要用来测定100 nm以上的大孔,液氮吸附实验主要测定孔径2 nm以上的孔隙,原子力显微镜,扫描电镜等主要用来观察孔隙结构和形态。而煤中小于2nm的纳米级孔隙报道较少。
高分辨率透射电子显微镜(HRTEM)的出现开启了研究物质的微观世界的大门,已经能检测到煤的大分子结构基本单元,并且定量研究,虽然对煤中晶格条纹取得了定量化效果,但是对于煤纳米孔隙还研究较少。
发明内容
本发明目的是提供一种基于HRTEM的煤中微小孔径的测定方法,通过HRTEM直观量化煤中微小孔的结构信息。
为解决上述技术问题,本发明所采用的技术方案为:
一种基于HRTEM的煤中微小孔径的测定方法,包括以下步骤:
步骤1),对煤样进行脱矿物质处理,然后对脱矿物质煤样进行高分辨透射电镜分析,对原始图像进行去噪,傅里叶-反傅里叶变换,图像的二值化处理等过程,在此之前,去除显微图像的边界,标记和比例尺。
步骤2),对得到的二值化图导入ImageJ软件中,截取含有晶格条纹信息的图像,去除区域之外的部分进行反转(Invert)。
步骤3),对反转后的图像进行分水岭分割算法(WaterShed)。
步骤4),对得到的图像按照原有比例尺进行设置,划线定义图像的长度,最终计算出每个分割后区块的面积,即为单个微小孔的面积。
步骤5),将得到的每个区块近似拟合为正方形,计算其边长,即可得到孔径分布。
进一步的,所述步骤3中,设置的反转区域为黑白二值化图所覆盖区域,不是电镜图像所有区域。
进一步的,所述步骤4中,反转之后的黑色区块可近似为煤中微小孔。
与现有技术相比本发明具有以下有益效果:
本发明在对煤中微小孔的研究中,摒弃了传统的方法,本发明方法简单,且本发明方法较为直观的再现了煤大分子基本结构单元中孔的分布情况,尤其是对于微小孔的具***置,对于研究煤中微小孔隙分布和煤大分子结构基本单元之间的关系起到了良好的指导作用。
附图说明
图1是本发明实施例所借助的杜儿坪2号煤的电镜原始图像及二值化图像。
图2是本发明实施例中图像处理过程(a是对二值化图进行分割,b是对分割后的二值化图像反转,c是对反转图进行分水岭算法得到若干个独立的黑色区块)。
图3是本发明实施例最终得到的微小孔径分布模型图。
图4为本发明实施例最终得到的微小孔孔径分布图。
具体实施方式
以下结合具体实施例对本发明作进一步说明。
实施例
本实施例以山西杜儿坪2号煤的高分辨率透射电镜图像处理过程进行说明。
取5g杜儿坪2号煤,基于HRTEM测定煤中微小孔径,进行以下步骤:
步骤1),对煤样进行脱矿物质处理,然后对脱矿物质煤样进行高分辨透射电镜分析,得到杜儿坪3号煤的电镜原始图像,如图1a所示。对原始图像进行去噪,傅里叶-反傅里叶变换,图像的二值化处理等过程,在此之前,去除显微图像的边界,标记和比例尺。
步骤2),对得到的二值化图导入ImageJ软件中,截取含有晶格条纹信息的图像,如图2a所示,去除区域之外的部分进行反转(Invert),得到图2b。
步骤3),对反转后的图像进行分水岭分割算法(WaterShed),得到图2c,可以清楚的看到分割后的孔的分布情况。
步骤4),对得到的图像按照原有比例尺进行设置,划线定义图像的长度,最终计算出每个分割后区块的面积,即为单个微小孔的面积。最终得到杜儿坪2号煤的孔隙分布模型图,如图3所示。
步骤5),将得到的每个区块近似拟合为正方形,计算其边长,即可得到孔径分布,如图4所示。
Claims (3)
1.一种基于HRTEM的煤中微小孔径的测定方法,其特征在于,包括以下步骤:
步骤1),对煤样进行脱矿物质处理,然后对脱矿物质煤样进行高分辨透射电镜分析,去除显微图像的边界,标记和比例尺后,对原始图像进行去噪,傅里叶-反傅里叶变换,图像的二值化处理;
步骤2),对得到的二值化图导入ImageJ软件中,截取含有晶格条纹信息的图像进行反转;
步骤3),对反转后的图像进行分水岭分割算法;
步骤4),对步骤3)得到的图像按照原有比例尺进行设置,划线定义图像的长度,最终计算出每个分割后区块的面积,即为单个微小孔的面积;
步骤5),将得到的每个区块近似拟合为正方形,计算其边长,即可得到孔径分布。
2.根据权利要求1所述的一种基于HRTEM的煤中微小孔径的测定方法,其特征在于,设置的反转区域为黑白二值化图所覆盖区域。
3.根据权利要求1所述的一种基于HRTEM的煤中微小孔径的测定方法,其特征在于,所述步骤4中反转之后的黑色区块近似为煤中微小孔。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010852224.0A CN111855537A (zh) | 2020-08-21 | 2020-08-21 | 一种基于hrtem的煤中微小孔径的测定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010852224.0A CN111855537A (zh) | 2020-08-21 | 2020-08-21 | 一种基于hrtem的煤中微小孔径的测定方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111855537A true CN111855537A (zh) | 2020-10-30 |
Family
ID=72970428
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010852224.0A Pending CN111855537A (zh) | 2020-08-21 | 2020-08-21 | 一种基于hrtem的煤中微小孔径的测定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111855537A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112461870A (zh) * | 2020-12-09 | 2021-03-09 | 中国矿业大学(北京) | 一种基于hrtem量化煤晶格条纹长度的方法 |
CN112862816A (zh) * | 2021-03-15 | 2021-05-28 | 太原理工大学 | 一种hrtem图像中煤芳香烃晶格条纹的智能提取方法 |
CN117409408A (zh) * | 2023-12-15 | 2024-01-16 | 北京大学 | 层理缝参数获取方法、装置、设备及可读存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102626628A (zh) * | 2012-05-17 | 2012-08-08 | 太原理工大学 | 一种用于合成气甲烷化催化剂的制备方法 |
WO2013144865A1 (en) * | 2012-03-27 | 2013-10-03 | University Of The Western Cape | Synthesis of zeolite x with hierarchical morphology from fly ash |
CN106525691A (zh) * | 2016-12-09 | 2017-03-22 | 河南理工大学 | 一种煤全孔径孔隙结构多数据融合的测定方法 |
CN109839401A (zh) * | 2019-01-29 | 2019-06-04 | 太原理工大学 | 一种采空区裂隙发育区的判定和处理方法 |
US20200005013A1 (en) * | 2018-06-29 | 2020-01-02 | Saudi Arabian Oil Company | Identifying geometrical properties of rock structure through digital imaging |
-
2020
- 2020-08-21 CN CN202010852224.0A patent/CN111855537A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013144865A1 (en) * | 2012-03-27 | 2013-10-03 | University Of The Western Cape | Synthesis of zeolite x with hierarchical morphology from fly ash |
CN102626628A (zh) * | 2012-05-17 | 2012-08-08 | 太原理工大学 | 一种用于合成气甲烷化催化剂的制备方法 |
CN106525691A (zh) * | 2016-12-09 | 2017-03-22 | 河南理工大学 | 一种煤全孔径孔隙结构多数据融合的测定方法 |
US20200005013A1 (en) * | 2018-06-29 | 2020-01-02 | Saudi Arabian Oil Company | Identifying geometrical properties of rock structure through digital imaging |
CN109839401A (zh) * | 2019-01-29 | 2019-06-04 | 太原理工大学 | 一种采空区裂隙发育区的判定和处理方法 |
Non-Patent Citations (3)
Title |
---|
王小令 等: ""基于HRTEM的煤中不同聚集态结构表征"", 《煤炭学报》 * |
谢婷 等: ""贵州遵义黔XY 1井龙马溪组页岩孔隙特征及主控因素"", 《石油实验地质》 * |
金智敏 等: ""煤岩CT图像的孔隙度和比表面积测量方法"", 《西安科技大学学报》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112461870A (zh) * | 2020-12-09 | 2021-03-09 | 中国矿业大学(北京) | 一种基于hrtem量化煤晶格条纹长度的方法 |
CN112862816A (zh) * | 2021-03-15 | 2021-05-28 | 太原理工大学 | 一种hrtem图像中煤芳香烃晶格条纹的智能提取方法 |
CN112862816B (zh) * | 2021-03-15 | 2024-03-15 | 太原理工大学 | 一种hrtem图像中煤芳香烃晶格条纹的智能提取方法 |
CN117409408A (zh) * | 2023-12-15 | 2024-01-16 | 北京大学 | 层理缝参数获取方法、装置、设备及可读存储介质 |
CN117409408B (zh) * | 2023-12-15 | 2024-03-08 | 北京大学 | 层理缝参数获取方法、装置、设备及可读存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111855537A (zh) | 一种基于hrtem的煤中微小孔径的测定方法 | |
Ziehmer et al. | A principle curvatures analysis of the isothermal evolution of nanoporous gold: Quantifying the characteristic length-scales | |
ElHadidy et al. | Development of a pore construction data analysis technique for investigating pore size distribution of ultrafiltration membranes by atomic force microscopy | |
Wierzbicki et al. | Mapping the complex morphology of cell interactions with nanowire substrates using FIB-SEM | |
Paredes et al. | Application of scanning tunneling and atomic force microscopies to the characterization of microporous and mesoporous materials | |
Wyart et al. | Membrane characterization by microscopic methods: multiscale structure | |
US20210074387A1 (en) | Simulation method for analyzing diffusion property of water-soluble monomer in hydrogel membrane | |
CN105352873A (zh) | 页岩孔隙结构的表征方法 | |
CN104574420A (zh) | 一种纳米级泥页岩数字岩心构建方法 | |
CN104751473B (zh) | 碳酸岩多尺度孔洞特征的确定方法及装置 | |
US11360037B1 (en) | Classified characterization method for connectivity of organic matter (OM)-hosted pores in shale | |
Babel et al. | Characterizing mineral wettabilities on a microscale by colloidal probe atomic force microscopy | |
CN110927194B (zh) | 确定泥页岩有机孔含量和孔径分布的方法 | |
Feng et al. | Influences of planarization modification and morphological filtering by AFM probe-tip on the evaluation accuracy of fractal dimension | |
CN112461870B (zh) | 一种基于hrtem量化煤晶格条纹长度的方法 | |
Laskaris et al. | AFM and SIMS surface and cation profile investigation of archaeological obsidians: New data | |
Blomqvist et al. | Interconnectivity imaged in three dimensions: Nano-particulate silica-hydrogel structure revealed using electron tomography | |
Katsiaounis et al. | Graphene nano-sieves by femtosecond laser irradiation | |
Fang et al. | Mechanics of bio-sediment transport | |
Palash et al. | An approach for quantitative analysis of pore size distribution of silica gel using atomic force microscopy | |
CN113570505B (zh) | 一种页岩三维超分辨率数字岩心分级重构方法及*** | |
CN114428040B (zh) | 一种页岩油储层储渗空间定量表征及参数获取方法 | |
CN112116583B (zh) | 一种基于sem图像处理的绝缘纸板老化判别检验方法 | |
CN114818542A (zh) | 一种毛细凝聚量关键参数孔隙结构凝聚比例的确定方法 | |
Dou et al. | Reducing molecular simulation time for AFM images based on super-resolution methods |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20201030 |
|
RJ01 | Rejection of invention patent application after publication |