CN109001236A - 一种基于纳米ct技术表征水泥浆体中孔隙率分布的方法 - Google Patents

一种基于纳米ct技术表征水泥浆体中孔隙率分布的方法 Download PDF

Info

Publication number
CN109001236A
CN109001236A CN201810719363.9A CN201810719363A CN109001236A CN 109001236 A CN109001236 A CN 109001236A CN 201810719363 A CN201810719363 A CN 201810719363A CN 109001236 A CN109001236 A CN 109001236A
Authority
CN
China
Prior art keywords
porosity distribution
cement slurry
porosity
nanometer
techniques
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
CN201810719363.9A
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.)
Southeast University
Original Assignee
Southeast University
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 Southeast University filed Critical Southeast University
Priority to CN201810719363.9A priority Critical patent/CN109001236A/zh
Publication of CN109001236A publication Critical patent/CN109001236A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/02Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
    • G01N23/04Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
    • G01N23/046Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
    • 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

Landscapes

  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • Pathology (AREA)
  • Immunology (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Pulmonology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Radiology & Medical Imaging (AREA)
  • Dispersion Chemistry (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本发明公开了一种基于纳米CT技术表征水泥浆体中孔隙率分布的方法,包括以下步骤:1)获取三维纳米CT灰度图像;2)对原始CT图像进行处理,将各像素的灰度值转换成孔隙率;3)计算孔隙率分布的配分函数;4)获取表征孔隙率分布的谱函数。本发明解决了现有技术中对孔结构造成测试破坏、需要预先假定孔隙具有规则几何形貌的问题,从而达到无损、精确表征水泥浆体中的孔隙率分布信息。

Description

一种基于纳米CT技术表征水泥浆体中孔隙率分布的方法
技术领域
本发明涉及一种数字图像信息分析技术,尤其涉及一种针对材料微观结构具有非均质分布特征的数字图像分析的方法。
背景技术
水泥浆体作为混凝土强度形成的关键,具有复杂的非均质微观结构。水泥浆体中的孔隙率分布对混凝土的强度及其它性能如介质传输均能产生重要的影响。水泥浆体中孔隙率分布的均匀程度越好,表明水泥浆体及混凝土具有更优良的性能,如更高的强度、更低的传输系数。如何准确定量地表征水泥浆体中的孔隙率分布对水泥基材料研究而言具有重要的理论与应用价值。
当前研究水泥浆体中孔隙率分布有多种方法,常用的测试手段包括压汞测试、扫描电镜等。压汞测试通过在不同静水压力下,将液态汞压入到水泥浆体内,记录汞压入时的压力与体积数据,基于Laplace方程获取水泥浆体中孔隙率分布的信息。扫描电镜通过观察具有理想平面的水泥浆体,对扫描电镜图片进行数字图像分析,获取水泥浆体中的孔隙率分布信息。值得注意的是,无论是压汞测试还是扫描电镜方法都存在一些显著的缺点。首先,这两种测试方法本身可能对样品造成一定的损害,带来测试数据的误差,如压汞测试在高压范围时、扫描电镜在样品制样时都可能破坏孔结构;其次,这两种测试方法需要假定水泥浆体中的孔隙具有规则的几何形貌,而实际上孔隙具有不规则的复杂几何形貌,规则几何形貌的假定导致基于水泥浆体孔结构预测其强度、介质传输等性能上与实验结果存在较大的偏差。
发明内容
针对上述存在的问题,本发明提供一种方法,解决了现有技术中孔结构造成测试破坏、需要预先假定孔隙具有规则几何形貌的问题,从而达到无损、精确表征水泥浆体中的孔隙率分布信息。
为了达到上述目的,本发明采用的技术方案如下:
一种基于纳米CT技术表征水泥浆体中孔隙率分布的方法,包括如下步骤:
1)获取三维纳米CT灰度图像,以h(j)表示第j像素的灰度值;
2)对原始CT图像进行处理,将各像素的灰度值转换成孔隙率P(j),转换关系式为P(j)=1-h(j)/255,其中常数255为CT灰度图像的最大像素值;
3)将三维CT图像划分大小相等的小立方体集合,以i表示小立方体的编号,以δ表示小立方体的大小,以Ωi(δ)表示小立方体所占据的空间,定义孔隙率分布的尺度变量mi(δ)为mi(δ)=∑Pj,j∈Ωi(δ),∑表示求和符号;
4)基于孔隙率分布的尺度变量mi(δ)定义孔隙率分布的概率测度
5)基于孔隙率分布的概率测度μi(δ)计算孔隙率分布的配分函数,其中q为取定指数,通常取为-10至10之间的整数值;
6)基于配分函数χ(q,δ)计算配分指数,其中符号≈表示当δ趋近于0时的取值,ln表示自然对数符号;
7)基于配分指数τ(q)计算孔隙率分布的奇异指数
8)基于配分指数τ(q)和奇异指数α(q)获取孔隙率分布的谱函数f(α)=q·α(q)-τ(q)。
本发明的有益效果是:
相较于常用的压汞测试及扫描电镜实验,本发明所述的纳米CT技术在样品制备和测试过程不对水泥浆体中的孔结构产生破坏性影响。同时,基于纳米CT数字图像分析获取的反映孔隙率分布的谱函数不需要假定水泥浆体中的孔隙具有规则的几何形貌,避免与孔隙真实的复杂不规则几何形貌产生偏差,达到无损、精确表征水泥浆体中的孔隙率分布信息的效果。
附图说明
图1为本发明实例纯水泥浆体三维纳米CT图像;
图2为本发明实例矿渣水泥浆体三维纳米CT图像;
图3为本发明实例纯水泥浆体孔隙率分布图像;
图4为本发明实例矿渣水泥浆体孔隙率分布图像;
图5为本发明实例纯水泥浆体孔隙率分布的配分函数;
图6为本发明实例矿渣水泥浆体孔隙率分布的配分函数;
图7为本发明实例纯水泥浆体和矿渣水泥浆体的配分指数;
图8为本发明实例纯水泥浆体和矿渣水泥浆体的奇异指数;
图9为本发明实例纯水泥浆体和矿渣水泥浆体的孔隙率分布谱函数。
具体实施方式
下面结合附图和具体实施方式,进一步阐明本发明,应理解下述具体实施方式仅用于说明本发明而不用于限制本发明的范围。
本发明提及的水泥浆体包括纯水泥浆体和矿渣水泥浆体。矿渣作为常用的矿物掺合料,其相对较低的水化速率对于降低水泥浆体的早期放热具有重要的作用。此外,大量研究表明,相较于纯水泥浆体,矿渣的掺入可以显著改善水泥浆体的孔结构,使孔隙分布更均匀。本实例将通过谱函数表征水泥浆体的孔隙率分布特征。本实例中的水泥浆体包含两种,一种是纯水泥浆体,另一种是60%水泥和40%矿渣的混合浆体。水灰比固定为0.4,水泥浆体样品在标准条件下养护56天。
实施例1:如图1、2、3、4、5、6、7、8和9所示的一种基于纳米CT技术表征水泥浆体中孔隙率分布的方法,包括如下:
(1)获取三维纳米CT灰度图像,以h(j)表示第j像素的灰度值:
如图1、2所示分别为纯水泥浆体和矿渣水泥浆体的三维纳米CT灰度图像,大小为128×128×128像素,每个像素的空间分辨率为1微米。
(2)对原始CT图像进行处理,将各像素的灰度值转换成孔隙率P(j):
利用MATLAB软件读取三维纳米CT灰度图像,基于转换关系式为P(j)=1-h(j)/255,其中常数255为CT灰度图像的最大像素值,将各像素的灰度值h(j)转换成孔隙率P(j),如图3、4所示分别为纯水泥浆体和矿渣水泥浆体的孔隙率分布图;
(3)将三维CT图像划分大小相等的小立方体集合,以i表示小立方体的编号,以δ表示小立方体的大小,以Ωi(δ)表示小立方体所占据的空间,定义孔隙率分布的尺度变量mi(δ)=∑Pj,j∈Ωi(δ),∑表示求和符号;
(4)基于孔隙率分布的尺度变量mi(δ)定义孔隙率分布的概率测度
(5)基于孔隙率分布的概率测度μi(δ)计算孔隙率分布的配分函数,如图5、6所示(针对两种不同的材料,图5和图6是有差别的,只是差别比较小),
(6)基于配分函数χ(q,δ)计算配分指数,其中符号≈表示当δ趋近于0时的取值,ln表示自然对数符号,如图7所示;
(7)如图8所示,基于配分指数τ(q)计算孔隙率分布的奇异指数
(8)基于配分指数τ(q)和奇异指数α(q)获取孔隙率分布的谱函数f(α)=q·α(q)-τ(q)。
如图8、9所示,矿渣水泥浆体孔隙率分布的谱函数宽度小于纯水泥浆体,说明矿渣水泥浆体中的孔隙率分布均匀性更好,与事实相一致,证明了该方法的有效性。
本实施例表明,本发明的一种基于纳米CT技术表征水泥浆体中孔隙率分布的方法,纳米CT技术不会对水泥浆体样品产生测试损害,更好地反映样品的原始微观结构;同时,采用谱函数表征孔隙率分布不需要假定水泥浆体中的孔隙具有规则几何形貌,避免不当假定带来的误差。
需要说明的是,本发明方案所公开的技术手段不仅限于上述实施方式所公开的技术手段,还包括由以上技术特征任意组合所组成的技术方案。

Claims (8)

1.一种基于纳米CT技术表征水泥浆体中孔隙率分布的方法,其特征在于,包括以下步骤:
1)获取三维纳米CT灰度图像,以h(j)表示第j像素的灰度值;
2)对原始CT图像进行处理,将各像素的灰度值转换成孔隙率P(j);
3)计算孔隙率分布的配分函数;
4)获取表征孔隙率分布的谱函数。
2.如权利要求1所述的基于纳米CT技术表征水泥浆体中孔隙率分布的方法,其特征在于,所述步骤2)的操作过程是对原始CT图像进行处理,将各像素的灰度值转换成孔隙率P(j),计算方法为:P(j)=1-h(j)/255,其中常数255为CT灰度图像的最大像素值。
3.如权利要求1所述的基于纳米CT技术表征水泥浆体中孔隙率分布的方法,其特征在于,所述步骤3)的操作过程中,将三维CT图像划分大小相等的小立方体集合,以i表示小立方体的编号,以δ表示小立方体的大小,以Ωi(δ)表示小立方体所占据的空间,定义孔隙率分布的尺度变量如下:mi(δ)=∑Pj,j∈Ωi(δ),其中∑表示求和符号。
4.如权利要求3所述的基于纳米CT技术表征水泥浆体中孔隙率分布的方法,其特征在于,所述步骤3)的操作过程中,基于孔隙率分布的尺度变量mi(δ)定义孔隙率分布的概率测度μi(δ)如下:
5.如权利要求4所述的基于纳米CT技术表征水泥浆体中孔隙率分布的方法,其特征在于,所述步骤3)的操作过程中,基于孔隙率分布的概率测度μi(δ)计算孔隙率分布的配分函数χ(q,δ)如下:
其中q为取定指数,通常取为-10至10之间的整数值。
6.如权利要求1所述的基于纳米CT技术表征水泥浆体中孔隙率分布的方法,其特征在于,所述步骤4)的操作过程中,基于配分函数χ(q,δ)计算配分指数τ(q)如下:
其中符号≈表示当δ趋近于0时的取值,ln表示自然对数符号。
7.如权利要求6述的基于纳米CT技术表征水泥浆体中孔隙率分布的方法,其特征在于,所述步骤4)的操作过程中,基于配分指数τ(q)计算孔隙率分布的奇异指数α(q)如下:
8.如权利要求7述的基于纳米CT技术表征水泥浆体中孔隙率分布的方法,其特征在于,所述步骤4)的操作过程中,基于配分指数τ(q)和奇异指数α(q)
获取表征孔隙率分布的谱函数f(α),计算方法如下:
f(α)=q·α(q)-τ(q)。
CN201810719363.9A 2018-09-03 2018-09-03 一种基于纳米ct技术表征水泥浆体中孔隙率分布的方法 Pending CN109001236A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810719363.9A CN109001236A (zh) 2018-09-03 2018-09-03 一种基于纳米ct技术表征水泥浆体中孔隙率分布的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810719363.9A CN109001236A (zh) 2018-09-03 2018-09-03 一种基于纳米ct技术表征水泥浆体中孔隙率分布的方法

Publications (1)

Publication Number Publication Date
CN109001236A true CN109001236A (zh) 2018-12-14

Family

ID=64598978

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810719363.9A Pending CN109001236A (zh) 2018-09-03 2018-09-03 一种基于纳米ct技术表征水泥浆体中孔隙率分布的方法

Country Status (1)

Country Link
CN (1) CN109001236A (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111272624A (zh) * 2019-11-29 2020-06-12 中国商用飞机有限责任公司 一种孔隙率检测方法
CN112129676A (zh) * 2019-06-24 2020-12-25 中国航发商用航空发动机有限责任公司 孔隙率试块的制作方法及孔隙率快速检测方法
CN112557254A (zh) * 2020-11-25 2021-03-26 东南大学 一种预测氯离子在水泥浆中有效扩散系数的方法
CN112881265A (zh) * 2021-03-01 2021-06-01 西南石油大学 固井水泥浆凝固过程孔隙连通性的定量原位评价方法
CN113567654A (zh) * 2021-06-18 2021-10-29 长江大学 一种用于评价气藏固井水泥石自愈合性能的实验方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101556703A (zh) * 2009-05-16 2009-10-14 中国石油大学(华东) 基于连续切片图像的网络模型建立方法
US8068579B1 (en) * 2008-04-09 2011-11-29 Xradia, Inc. Process for examining mineral samples with X-ray microscope and projection systems
CN102944568A (zh) * 2012-12-03 2013-02-27 东南大学 工业x-ct用水泥基试件原位加载仪及其使用方法
WO2014142976A1 (en) * 2013-03-15 2014-09-18 Bp Corporation North America Inc. Systems and methods for improving direct numerical simulation of material properties from rock samples and determining uncertainty in the material properties
CN105487121A (zh) * 2015-12-03 2016-04-13 长江大学 基于ct扫描图像与电成像图像融合构建多尺度数字岩心方法
CN105806765A (zh) * 2016-04-13 2016-07-27 南京大学(苏州)高新技术研究院 一种显微ct扫描土体空间孔隙结构的精细化表征方法
CN107655908A (zh) * 2017-11-07 2018-02-02 中国石油天然气股份有限公司 一种构建数字岩心的方法及装置
CN107817199A (zh) * 2016-09-14 2018-03-20 中国石油化工股份有限公司 一种致密砂岩多尺度孔隙模型的构建方法和应用
CN108311409A (zh) * 2017-12-26 2018-07-24 江苏仲博敬陈信息科技有限公司 一种基于人工智能图像识别的印刷电路板检测系统
CN108387495A (zh) * 2018-01-22 2018-08-10 青岛理工大学 一种多孔混凝土孔隙率计算和孔隙参数表征方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8068579B1 (en) * 2008-04-09 2011-11-29 Xradia, Inc. Process for examining mineral samples with X-ray microscope and projection systems
CN101556703A (zh) * 2009-05-16 2009-10-14 中国石油大学(华东) 基于连续切片图像的网络模型建立方法
CN102944568A (zh) * 2012-12-03 2013-02-27 东南大学 工业x-ct用水泥基试件原位加载仪及其使用方法
WO2014142976A1 (en) * 2013-03-15 2014-09-18 Bp Corporation North America Inc. Systems and methods for improving direct numerical simulation of material properties from rock samples and determining uncertainty in the material properties
CN105487121A (zh) * 2015-12-03 2016-04-13 长江大学 基于ct扫描图像与电成像图像融合构建多尺度数字岩心方法
CN105806765A (zh) * 2016-04-13 2016-07-27 南京大学(苏州)高新技术研究院 一种显微ct扫描土体空间孔隙结构的精细化表征方法
CN107817199A (zh) * 2016-09-14 2018-03-20 中国石油化工股份有限公司 一种致密砂岩多尺度孔隙模型的构建方法和应用
CN107655908A (zh) * 2017-11-07 2018-02-02 中国石油天然气股份有限公司 一种构建数字岩心的方法及装置
CN108311409A (zh) * 2017-12-26 2018-07-24 江苏仲博敬陈信息科技有限公司 一种基于人工智能图像识别的印刷电路板检测系统
CN108387495A (zh) * 2018-01-22 2018-08-10 青岛理工大学 一种多孔混凝土孔隙率计算和孔隙参数表征方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
刘数华 等: "《混凝土辅助胶凝材料》", 30 June 2010, 中国建材工业出版社 *
张云强 等: "铁谱磨粒多重分形特征研究", 《润滑与密封》 *
方海平 等: "五大连池玄武岩三维孔隙组构的多重分形特征", 《地质科技情报》 *
王毅: "《块状褐煤高温蒸汽热解的宏细观特性分析及应用》", 31 March 2012, 中国矿业大学出版社 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112129676A (zh) * 2019-06-24 2020-12-25 中国航发商用航空发动机有限责任公司 孔隙率试块的制作方法及孔隙率快速检测方法
CN112129676B (zh) * 2019-06-24 2023-09-22 中国航发商用航空发动机有限责任公司 孔隙率试块的制作方法及孔隙率快速检测方法
CN111272624A (zh) * 2019-11-29 2020-06-12 中国商用飞机有限责任公司 一种孔隙率检测方法
CN112557254A (zh) * 2020-11-25 2021-03-26 东南大学 一种预测氯离子在水泥浆中有效扩散系数的方法
CN112557254B (zh) * 2020-11-25 2022-04-15 东南大学 一种预测氯离子在水泥浆中有效扩散系数的方法
CN112881265A (zh) * 2021-03-01 2021-06-01 西南石油大学 固井水泥浆凝固过程孔隙连通性的定量原位评价方法
CN112881265B (zh) * 2021-03-01 2022-04-12 西南石油大学 固井水泥浆凝固过程孔隙连通性的定量原位评价方法
CN113567654A (zh) * 2021-06-18 2021-10-29 长江大学 一种用于评价气藏固井水泥石自愈合性能的实验方法
CN113567654B (zh) * 2021-06-18 2024-02-27 长江大学 一种用于评价气藏固井水泥石自愈合性能的实验方法

Similar Documents

Publication Publication Date Title
CN109001236A (zh) 一种基于纳米ct技术表征水泥浆体中孔隙率分布的方法
Lü et al. Fractal dimension of concrete incorporating silica fume and its correlations to pore structure, strength and permeability
Wang et al. X-ray computed tomography for pore-related characterization and simulation of cement mortar matrix
Brisard et al. Morphological quantification of hierarchical geomaterials by X-ray nano-CT bridges the gap from nano to micro length scales
CN104685382B (zh) 可靠地预测孔隙渗透性趋势的数字岩石分析系统及方法
CN108763711A (zh) 一种基于岩心扫描图像分块数值模拟的渗透率预测方法
CN110415167A (zh) 一种基于数字图像技术的粗糙面裂隙生成方法及试验系统
CN110133035B (zh) 核磁共振耦合恒速压汞定量表征砂岩储层孔喉结构的方法
Xu et al. Microstructure and water absorption of ancient concrete from Pompeii: An integrated synchrotron microtomography and neutron radiography characterization
CN110210176B (zh) 一种花岗岩微观裂缝结构识别和微观渗流分析方法
Gao et al. The multifractal property of heterogeneous microstructure in cement paste
Zhao et al. Digital measurement of 2D and 3D cracks in sandstones through improved pseudo color image enhancement and 3D reconstruction method
CN113075102B (zh) 一种建立多孔介质自发渗吸量与时间关系数学模型的方法
CN108535162A (zh) 一种基于x射线显微ct表征木材孔隙结构的方法
Zhang et al. A testing assembly for combination measurements on gas hydrate-bearing sediments using x-ray computed tomography and low-field nuclear magnetic resonance
CN111537544B (zh) 一种提高核磁共振t2谱表征致密储层孔径分布精度的转换方法
CN108061697B (zh) 土体三维孔隙率计算方法
Zhang et al. Porosity models for determining the pore-size distribution of rocks and soils and their applications
Arnaud et al. Imaging of firn and bubbly ice in coaxial reflected light: a new technique for the characterization of these porous media
Li et al. Three-dimensional multiscale fusion for porous media on microtomography images of different resolutions
CN111504875A (zh) 一种致密砂岩孔喉参数的提取及计算方法及提取装置
CN108827849B (zh) 一种基于超临界碳化的水泥砂浆孔隙率标定的方法
CN111006985B (zh) 地质条件下陆相湖盆致密储层孔喉有效性定量评价方法
Cerepi et al. Petrophysical properties of porous medium from petrographic image analysis data
Yio Characterising the microstructure of cement-based materials using laser scanning confocal microscopy

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: 20181214

RJ01 Rejection of invention patent application after publication