CN107063961A - 一种水泥水化三维孔隙结构特征的表征方法 - Google Patents

一种水泥水化三维孔隙结构特征的表征方法 Download PDF

Info

Publication number
CN107063961A
CN107063961A CN201611160205.1A CN201611160205A CN107063961A CN 107063961 A CN107063961 A CN 107063961A CN 201611160205 A CN201611160205 A CN 201611160205A CN 107063961 A CN107063961 A CN 107063961A
Authority
CN
China
Prior art keywords
hole
hydrated cementitious
sample
characterizing method
structures feature
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
CN201611160205.1A
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.)
China University of Mining and Technology CUMT
Original Assignee
China University of Mining and Technology CUMT
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 China University of Mining and Technology CUMT filed Critical China University of Mining and Technology CUMT
Priority to CN201611160205.1A priority Critical patent/CN107063961A/zh
Publication of CN107063961A publication Critical patent/CN107063961A/zh
Pending legal-status Critical Current

Links

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

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

本发明提供了一种水泥水化三维孔隙结构特征的表征方法,先取水泥基材料制样,获取试样的背散射和能谱图片,通过“Bentz分相规则”得到试样各矿物相参数,基于CEMHYD3D水化模型软件,建立试样的水化模型;然后分别通过“内侵蚀法”、“燃烧算法”和“随机走法”三种算法实现对水泥水化三维孔隙结构特征的表征。本发明针对现有实验手段无法测试出重要的特征参数,无法进行连续监测,实验周期长以及对试件产生破坏等问题,提供了一种可连续监测、可视化、无破坏、省时的水泥水化三维孔隙结构特征的表征方法。

Description

一种水泥水化三维孔隙结构特征的表征方法
技术领域
本发明属于微观数值试验技术领域,具体涉及一种水泥水化三维孔隙结构特征的表征方法。
背景技术
水泥基材料在水化过程中形成的孔隙对其诸多性能具有重要影响,尤其是其作为水分及恶劣环境中各种离子传输的通道,腐蚀材料内部造成破坏。因此,充分了解水泥基材料微结构中的各种孔结构特征,有助于分析水泥基材料的传输性能与耐久性。
目前,随着各种精密仪器和显微技术的成熟,已经出现了各种用于表征孔隙结构的新技术,对于水化过程和微结构的演变的研究则主要是通过实验方法进行。现在较多使用的测试手段包括小角射线散射法、压汞法、磁共振测孔法、同步加速 X 射线层析扫描测孔法等。但由于各种条件以及技术的限制,现有的各种测试手段都具有一定的局限性,比如样品制备过程要求高、难度大;测试的周期较长、费用高;无法进行连续可视化监测;实验过程中可能对样品产生部分破坏以及部分参数无法表征等缺陷。随着计算机技术的发展以及对水泥水化演变机理认识的不断深入,利用计算机技术模拟水泥水化演变过程的水泥水化模型得到不断优化发展。其中的 CEMHYD3D 模型、HYMOSTRUC3D 模型、DuCOM模型、 μic 模型以及 THAMES 模型是比较成功的水泥水化模型。
CEMHYD3D水化模型充分考虑了各水泥物相的反应,利用元胞自动机技术控制水泥像素的溶解、扩散及成核反应过程,考虑其代码开源的优点,本发明基于数字图像的CEMHYD3D水化模型进行三维孔隙结构特征的表征。
发明内容
解决的技术问题:本发明的目的是针对现有实验手段无法测试出重要的特征参数,无法进行连续监测,实验周期长以及对试件产生破坏等问题,提供了一种可连续监测、无破坏、省时的水泥水化三维孔隙结构特征的表征方法。
技术方案:一种水泥水化三维孔隙结构特征的表征方法,包括以下步骤:
步骤1,取水泥基材料制样,获取试样的背散射和能谱图片,通过“Bentz分相规则”得到试样中各矿物相参数;
步骤2,根据步骤1所得试样的各矿物相参数,基于CEMHYD3D水化模型软件,建立试样的水化三维可视化模型;
步骤3,根据步骤2建立的水化模型,对微结构像素点进行扫描以统计孔隙结构的孔隙率,通过“内侵蚀法”得到孔隙结构的孔径累积分布曲线;
步骤4,根据步骤2建立的水化模型,通过“燃烧算法”得到孔隙结构的封闭孔、半连通和连通孔孔隙率;
步骤5,根据步骤2建立的水化模型,通过“随机走法”计算孔隙结构的曲折度。
进一步地,步骤1中制样是先取水泥粉体和低粘度环氧树脂混合成粘稠状浆体,压入模具养护后脱模,经切割成厚度1cm圆柱体;然后经打磨、抛光和表面喷碳,即得。
进一步地,步骤1中获取试样的背散射和能谱图片是在加速电压12kV,管电流2nA获取背散射图片,管电流10nA 获取能谱图片。
进一步地,步骤1中“Bentz分相规则”是先通过能谱图片的灰度直方图选取灰度直方图右端肩状部位置为阈值,灰度值大于此值的像素点表示存在该元素,否则不存在该元素;然后,对背散射图片采用中值滤波进行降噪处理,划分出图片中的固相、孔隙,结合分相树进行四种矿物相与石膏的划分。
进一步地,步骤1中所述各矿物相参数是指各矿物相的体积分数与表面积分数。
进一步地,步骤3中“内侵蚀法”是先将图片中的孔相和固相进行区分;再对微结构依次逐层进行扫描,当该像素为孔相且其至少有一个面与固相相接处,则将此孔孔像素点亮,标记为1,当该孔像素点周边都为孔相,且已进行标记,则该孔像素标记成已标记数字的周边孔相最小值加1,以此类推,直到所有孔像素被标记;最后搜索标记为最大值的数字,当其邻近位置存在相同标记的像素时,则其直径为标记数字的2倍,否则为2倍减1。
进一步地,步骤4中“燃烧算法”是在微结构第一层进行搜索,当检索到孔位置时,将其标记并检查六个临近像素点,当有孔隙存在时再次进行标记,以此类推,直至某一像素点的六个临近位置为固相,此时可以统计得到连通孔和半连通孔的体积,由总孔隙的体积减去连通孔和半连通孔的体积即为微结构中封闭孔的体积。
进一步地,步骤5中“随机走法”是随机产生一定数量的数组,在微结构中随机进行坐标的行走,直到一定的步数之后,统计均方位移和步数,随后通过计算均方位移和步数的比值作为曲折度。
有益效果:相比于实验技术手段的测试,本发明的水泥水化三维孔隙结构特征的表征方法具有如下优点:首先,建立了可连续化监测水泥水化三维孔结构演变的过程;其次,相对于实验手段,本发明不会产生有害影响;最后,本发明建立了三维孔结构特征的参数表征手段,实现了参数的物理模型表示并得到参数结果。
附图说明
图1为实施例1中试样未水化的CEMHYD3D水化模型微结构示意图,及试样水化后的CEMHYD3D水化模型微结构示意图;
图2为实施例1中试样在水化进行过程中孔隙率的变化曲线;
图3为实施例1中试样在水化7天后微结构的孔径积分与微分分布曲线;
图4为实施例1中试样在水化进行过程中连通度的变化曲线;
图5为实施例1中试样在水化进行过程中总孔隙率、连通孔、半连通孔及封闭孔孔隙率的变化曲线;
图6为实施例1中试样在水化进行过程中曲折度的变化曲线。
具体实施方式
下面结合实施例和附图对本发明做更进一步地解释。下列实施例中仅用于说明本发明,但并不用来限定本发明的实施范围。
孔隙率,即微结构中孔隙占总体积的比例。
封闭孔、连通孔、半连通孔的孔隙率为各类型孔占总孔隙的比例。
孔径积分分布曲线为相应直径及小于此直径的各孔隙体积之和占总孔隙体积的比例。
孔径微分分布曲线为相应直径的孔隙体积占总孔隙体积的比例。
曲折度,即表征内部孔隙结构曲折程度的参数,使用均方位移与步数的比值来表示。
实施例1
水灰比为0.35时,水泥水化三维孔隙结构特征的表征:
步骤1,选用的水泥信息为:PI 525.普通硅酸盐水泥,密度3.15 g/cm3,比表面积 370m2/kg,初凝时间2小时12分,终凝时间3小时7分。称取约 25g水泥粉体,将其与低粘度环氧树脂混合成粘稠状浆体,压入圆柱形模具,在60℃条件下养护24h 脱模,切割成厚度1cm 圆柱体,然后在磨样机打磨 2 min,使用0.1μm的金刚石悬浮液抛光60min,放到乙醇中,随后进行表面喷碳,得到试样。设定加速电压12 kV,管电流2 nA获取背散射图片,设定管电流10nA 获取各元素能谱图片。
通过能谱图片的灰度直方图选取合适的阈值(一般选取灰度直方图右端肩状部位置),灰度值大于此值的像素点表示存在该元素,否则不存在该元素;然后,对背散射图片采用中值滤波进行降噪处理,划分出图片中的固相、孔隙,结合分相树进行四种矿物相与石膏的划分,得到试样中各矿物相参数。统计得到的C3S、C2S、C3A、 C4AF的体积分数分别为52.36%、29.75%、4.77%、13.12%,表面积分数分别为48.82%、34.98%、4.96%、11.24%。分相后未水化的微结构如图1中Cycle=0对应的图。
步骤2,将水泥颗粒的相关信息和设定条件,连同步骤1中计算出来的体积分数与表面积分数作为参数,输入CEMHYD3D水化模型软件,执行水化程序,得到水灰比为0.35 时的水泥水化模型三维可视化微结构示意图。分相后水化的微结构如图2中Cycle=500对应的图。
步骤3,将建立的水化微结构图输入到开发的孔径统计程序中,程序对微结构像素点进行扫描的同时对各孔相进行标记,随后进行孔径孔相统计与孔径累计分数的计算,统计的孔隙率如图3,计算的水化7天之后孔径累积分布与微分分布结果如图4。
步骤4,将建立的水化微结构图输入到开发的孔分类程序中,程序首先应用“燃烧算法”的思想将孔进行标记,随后进行连通孔、封闭孔的统计,并计算连通度。计算的连通度如图5,连通孔和封闭孔的孔隙率如图6。
步骤5,将建立的水化微结构图输入到开发的曲折度计算程序中,控制生成的数组随机进行移动,当移动的步数达到设定数量时,统计行走之后的最终位置与最初位置之间的均方位移,通过计算均方位移与步数之间的比值来表示曲折度的大小。计算的曲折度大小如图6。
对图2进行分析,可以看出随着水化过程的进行,水泥水化微结构内部的孔隙随水化程度逐渐减少,且具有线性关系;
对图3进行分析,可以看出水化7天之后,微结构内部孔隙的孔径的范围在 1um 到14um之间,其中,孔径多集中在 3—8um,孔径为 5um 的孔隙占的体积最多;
对图4进行分析,可以看出随着水化时间的进行,微结构内部孔隙的连通度逐渐在下降并直到内部不存在连通的孔;
对图5进行分析,可以看出随着水化时间的进行,总孔隙率和连通孔孔隙率逐渐下降,封闭孔孔隙率则上升。半连通孔孔隙率先上升后下降,可能原因为前期水化产物填充孔隙阻断孔的连通形成半连通孔,而随着水化的进行,水化产物逐渐增多,占据的孔隙越来越多导致半连通孔的数量也在下降;
对图6进行分析,可以看出随着水化进行曲折度在上升,但在水化18天时,其曲折度出现减小,因为在水化18天之后,内部水化产物阻断了与连通孔相连的周边延伸孔隙,由此而造成连通孔弯曲程度减小。

Claims (8)

1.一种水泥水化三维孔隙结构特征的表征方法,其特征在于:包括以下步骤:
步骤1,取水泥基材料制样,获取试样的背散射和能谱图片,通过“Bentz分相规则”得到试样各矿物相参数;
步骤2,根据步骤1所得试样的各矿物相参数,基于CEMHYD3D水化模型软件,建立试样的水化模型;
步骤3,根据步骤2建立的水化模型,对微结构像素点进行扫描以统计孔隙结构的孔隙率,通过“内侵蚀法”得到孔隙结构的孔径累积分布曲线;
步骤4,根据步骤2建立的水化模型,通过“燃烧算法”得到孔隙结构的封闭孔、半连通和连通孔孔隙率;
步骤5,根据步骤2建立的水化模型,通过“随机走法”计算孔隙结构的曲折度。
2.根据权利要求1所述的水泥水化三维孔隙结构特征的表征方法,其特征在于:步骤1中制样是先取水泥粉体和低粘度环氧树脂混合成粘稠状浆体,压入模具养护后脱模,经切割成厚度1cm圆柱体;然后经打磨、抛光和表面喷碳,即得。
3.根据权利要求1所述的水泥水化三维孔隙结构特征的表征方法,其特征在于:步骤1中获取试样的背散射和能谱图片是在加速电压12kV,管电流2nA 获取背散射图片,管电流10nA 获取能谱图片。
4.根据权利要求1所述的水泥水化三维孔隙结构特征的表征方法,其特征在于:步骤1中“Bentz分相规则”是先通过能谱图片的灰度直方图选取灰度直方图右端肩状部位置为阈值,灰度值大于此值的像素点表示存在该元素,否则不存在该元素;然后,对背散射图片采用中值滤波进行降噪处理,划分出图片中的固相、孔隙,结合分相树进行四种矿物相与石膏的划分。
5.根据权利要求1所述的水泥水化三维孔隙结构特征的表征方法,其特征在于:步骤1中所述各矿物相参数是指各矿物相的体积分数与表面积分数。
6.根据权利要求1所述的水泥水化三维孔隙结构特征的表征方法,其特征在于:步骤3中“内侵蚀法”是先将图片中的孔相和固相进行区分;再对微结构依次逐层进行扫描,当该像素为孔相且其至少有一个面与固相相接处,则将此孔孔像素点亮,标记为1,当该孔像素点周边都为孔相,且已进行标记,则该孔像素标记成已标记数字的周边孔相最小值加1,以此类推,直到所有孔像素被标记;最后搜索标记为最大值的数字,当其邻近位置存在相同标记的像素时,则其直径为标记数字的2倍,否则为2倍减1。
7.根据权利要求1所述的水泥水化三维孔隙结构特征的表征方法,其特征在于:步骤4中“燃烧算法”是在微结构第一层进行搜索,当检索到孔位置时,将其标记并检查六个临近像素点,当有孔隙存在时再次进行标记,以此类推,直至某一像素点的六个临近位置为固相,此时可以统计得到连通孔和半连通孔的体积,由总孔隙的体积减去连通孔和半连通孔的体积即为微结构中封闭孔的体积。
8.根据权利要求1所述的水泥水化三维孔隙结构特征的表征方法,其特征在于:步骤5中“随机走法”是随机产生数组,在微结构中随机进行坐标的行走,统计均方位移和步数,随后通过计算均方位移和步数的比值作为曲折度。
CN201611160205.1A 2016-12-15 2016-12-15 一种水泥水化三维孔隙结构特征的表征方法 Pending CN107063961A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611160205.1A CN107063961A (zh) 2016-12-15 2016-12-15 一种水泥水化三维孔隙结构特征的表征方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611160205.1A CN107063961A (zh) 2016-12-15 2016-12-15 一种水泥水化三维孔隙结构特征的表征方法

Publications (1)

Publication Number Publication Date
CN107063961A true CN107063961A (zh) 2017-08-18

Family

ID=59619067

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611160205.1A Pending CN107063961A (zh) 2016-12-15 2016-12-15 一种水泥水化三维孔隙结构特征的表征方法

Country Status (1)

Country Link
CN (1) CN107063961A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112150614A (zh) * 2020-09-24 2020-12-29 济南大学 基于微观结构图像的水泥三维仿真模型建立方法及系统
CN113237809A (zh) * 2021-04-16 2021-08-10 贵州电网有限责任公司 一种复合绝缘子芯棒孔隙率评估方法
CN113405965A (zh) * 2021-06-08 2021-09-17 浙江广天构件集团股份有限公司 一种水泥基材料颗粒堆积体系孔隙连通性分析方法
CN116297581A (zh) * 2023-03-06 2023-06-23 中国长江三峡集团有限公司 一种混凝土中方镁石反应程度的定量分析方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103940721A (zh) * 2014-04-17 2014-07-23 山东大学 一种新拌多孔混凝土材料孔隙率快速测定方法
CN104596911A (zh) * 2015-02-05 2015-05-06 山东大学 一种基于图像处理的透水混凝土孔隙特征的检测方法
CN105259335A (zh) * 2015-10-26 2016-01-20 中国建筑材料科学研究总院 水泥熟料性能的评价方法及煅烧工艺优化方法
CN105352873A (zh) * 2015-11-26 2016-02-24 中国石油大学(北京) 页岩孔隙结构的表征方法
CN105372168A (zh) * 2015-12-08 2016-03-02 清华大学 一种混凝土孔隙三维结构的构建方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103940721A (zh) * 2014-04-17 2014-07-23 山东大学 一种新拌多孔混凝土材料孔隙率快速测定方法
CN104596911A (zh) * 2015-02-05 2015-05-06 山东大学 一种基于图像处理的透水混凝土孔隙特征的检测方法
CN105259335A (zh) * 2015-10-26 2016-01-20 中国建筑材料科学研究总院 水泥熟料性能的评价方法及煅烧工艺优化方法
CN105352873A (zh) * 2015-11-26 2016-02-24 中国石油大学(北京) 页岩孔隙结构的表征方法
CN105372168A (zh) * 2015-12-08 2016-03-02 清华大学 一种混凝土孔隙三维结构的构建方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
M.A.B. PROMENTILLA ET AL.: "qantification of tortuosity in hardened cement pastes using synchrotron-based X-ray computed microtomography", 《CEMENT AND CONCRETE RESEARCH》 *
MINGZHONG ZHANG ET AL.: "Computational investigation on mass diffusivity in Portland cement paste based on X-ray computed microtomography (μCT) image", 《CONSTRUCTION AND BUILDING MATERIALS》 *
SHASHANLC BISHNOI ET AL.: "μic: A new platform for modelling the hydration of cements", 《CEMENT AND CONCRETE RESEARCH》 *
谢德擎: "基于不规则颗粒的水泥基材料水化进程及其传输性能的数值模拟", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 *
魏小胜: "《用电阻率表征水泥混凝土结构形成动力学性能》", 30 September 2016, 武汉理工大学出版社 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112150614A (zh) * 2020-09-24 2020-12-29 济南大学 基于微观结构图像的水泥三维仿真模型建立方法及系统
CN113237809A (zh) * 2021-04-16 2021-08-10 贵州电网有限责任公司 一种复合绝缘子芯棒孔隙率评估方法
CN113405965A (zh) * 2021-06-08 2021-09-17 浙江广天构件集团股份有限公司 一种水泥基材料颗粒堆积体系孔隙连通性分析方法
CN116297581A (zh) * 2023-03-06 2023-06-23 中国长江三峡集团有限公司 一种混凝土中方镁石反应程度的定量分析方法
CN116297581B (zh) * 2023-03-06 2024-04-26 中国长江三峡集团有限公司 一种混凝土中方镁石反应程度的定量分析方法

Similar Documents

Publication Publication Date Title
CN107063961A (zh) 一种水泥水化三维孔隙结构特征的表征方法
Shan et al. Pore characteristics of pervious concrete and their influence on permeability attributes
CN103364241B (zh) 重塑细粒土试样制备装置
CN106447776A (zh) 基于3d打印制作的复杂裂隙岩体物理模型及建模方法
Fu et al. Experimental study on pore characteristics and fractal dimension calculation of pore structure of aerated concrete block
Liu et al. Concrete damage evolution and three-dimensional reconstruction by integrating CT test and fractal theory
CN104751473A (zh) 碳酸岩多尺度孔洞特征的确定方法及装置
CN113063810B (zh) 获得砂岩冻融作用下宏细观损伤演化规律的方法
CN109632429B (zh) 一种土石混合体双轴压缩试验的试样制备方法
CN103076347A (zh) 基于原位x射线断层照相的脆性材料力学损伤的测量方法
CN104359763A (zh) 一种水泥基材料在荷载作用下内部裂缝开展的检测方法
Zhang et al. Effects of initial defects within mortar cover on corrosion of steel and cracking of cover using X-ray computed tomography
CN116631536A (zh) 图像处理视域下透水混凝土集料均匀性指标计算方法
CN112270093A (zh) 一种基于微地震矩张量反演的储层天然裂缝建模方法
CN116050917A (zh) 一种基于网格分布特征的路面三维纹理均匀性评价方法
Liu et al. Digital concrete modelling: an alternative approach to microstructural pore analysis of cement hydrates
Tian et al. Corrosion damages of reinforced concrete characterized by X-ray CT and DVC techniques
CN116754587A (zh) 三维层面的再生骨料表面残余砂浆层分布的量化表征方法
CN110599477A (zh) 一种植草混凝土有效孔隙率的检测方法及系统
CN113405962B (zh) 碳酸盐岩储层的储集空间评价方法和装置
CN110988004B (zh) 一种粒间压溶成因石英胶结物含量的评价方法
WANG et al. Mechanical characteristics and fracture mechanism of 3D printed rock samples with cross joints
Zhang et al. Study on the non-uniformity of rust layer and residual rebar caused by pitting corrosion in mortar
Ma et al. Micro-constructive Damage Model of Dry-Wet Cyclic Red Clay Based on Weibull Distribution
Xie et al. Study of the Mesodamage Analysis Method for Frozen–Thawed Concrete Based on CT Image Recognition

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

RJ01 Rejection of invention patent application after publication