CN111855537A - 一种基于hrtem的煤中微小孔径的测定方法 - Google Patents

一种基于hrtem的煤中微小孔径的测定方法 Download PDF

Info

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
Application number
CN202010852224.0A
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.)
Taiyuan University of Technology
Original Assignee
Taiyuan University of Technology
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 Taiyuan University of Technology filed Critical Taiyuan University of Technology
Priority to CN202010852224.0A priority Critical patent/CN111855537A/zh
Publication of CN111855537A publication Critical patent/CN111855537A/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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/62Analysis of geometric attributes of area, perimeter, diameter or volume
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10056Microscopic image
    • G06T2207/10061Microscopic image from scanning electron microscope
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20048Transform domain processing
    • G06T2207/20056Discrete and fast Fourier transform, [DFT, FFT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20112Image segmentation details
    • G06T2207/20152Watershed 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的煤中微小孔径的测定方法
技术领域
本发明涉及煤孔隙结构测定技术领域,具体涉及一种基于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中反转之后的黑色区块近似为煤中微小孔。
CN202010852224.0A 2020-08-21 2020-08-21 一种基于hrtem的煤中微小孔径的测定方法 Pending CN111855537A (zh)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
王小令 等: ""基于HRTEM的煤中不同聚集态结构表征"", 《煤炭学报》 *
谢婷 等: ""贵州遵义黔XY 1井龙马溪组页岩孔隙特征及主控因素"", 《石油实验地质》 *
金智敏 等: ""煤岩CT图像的孔隙度和比表面积测量方法"", 《西安科技大学学报》 *

Cited By (5)

* Cited by examiner, † Cited by third party
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