摘要
利用X射线计算机断层成像及图像三维重构技术获取了陶瓷透水路面砖的孔隙特征参数与渗透系数.基于沿程水头损失及局部水头损失影响因素,选取合适孔隙特征参数,并利用量纲分析法建立了渗透系数与孔隙特征参数的定量关系模型.结果表明:通过统计学检验方法对建立的定量关系模型进行检验,发现所建立模型是合理可靠的;渗透系数受有效孔隙迂曲度、有效孔隙半径均值的影响最大,受有效孔隙比表面积、有效孔隙等效喉道半径及有效孔隙等效喉道长度的影响较小.
随着海绵城市战略的提出,透水路面砖受到越来越多的关注,已被广泛应用于人行道、公园及停车场等区域.其优良的抗滑和降噪特性为人们出行提供了安全与舒适保障条件.为此国内外学者对透水路面材料的性能展开了大量研究.
崔新壮
本文拟利用X射线计算机断层成像(X‑ray computed tomography,XCT)与图像三维重构技术,获得孔隙三维模型,提取孔隙特征参数.基于沿程水头损失及局部水头损失影响机制,对孔隙特征参数进行相关性分析.然后,结合量纲分析法建立渗透系数与孔隙特征参数之间的定量关系模型,并基于孔隙模型的渗流数值模拟结果进行验证,希望为透水路面材料的应用及优化选用提供理论参考.
以陶瓷透水路面砖为研究对象,规格为200 mm(长)×100 mm(宽)×55 mm(高)的长方体,抗压强度为40 MPa,抗折强度大于6 MPa,透水系数大于0.2 mm/s,耐磨性能满足磨坑长度小于35 mm,抗滑性能满足防滑指标(BPN值)大于60,其他性能指标符合GB/T 25993—2010《透水路面砖和透水路面板》要求.透水砖烧制流程为:在透水砖高度方向进行两次布料和一次烧结成型,上层为细集料,下层为粗集料,从而形成上层小而下层大的孔隙结构,其中,粗细集料均为废旧陶瓷经过高温煅烧之后的陶瓷颗粒.当淤堵物随雨水从表面进入透水砖内部时,结构的反滤特性可满足透水砖透水功能的长期有效性.
试样制样方法:为获得较高精度孔隙结构信息,对整砖试件进行钻芯取样,得到保留原整砖试样高度方向结构分布的圆柱芯样,其直径为50 mm,高度为55 mm,如

图1 透水砖圆柱芯样(截取自实际试样)
Fig.1 Cylindrical core sample of pervious brick(actual specimens intercepted)
透水路面材料中的孔隙主要分为封闭孔隙、半封闭孔隙和连通孔隙3类.其中连通孔隙又称为有效孔隙,对透水路面材料渗透系数的影响最大.另外两类孔隙对于透水路面材料的渗透系数影响较小.故本文主要讨论陶瓷透水路面砖中有效孔隙及其特征参数,如无特别说明,下文中孔隙特征参数均指有效孔隙特征参数.
将扫描得到的图像进行三维重构,并对其进行滤波处理,得到圆柱芯样的灰度切片图像.将灰度图像导入Avizo软件进行阈值分割,将孔隙相与集料相分离.最终得到透水砖圆柱芯样的三维孔隙模型,如

图2 透水砖圆柱芯样的三维孔隙模型
Fig.2 3D pore model of pervious brick cylindrical core sample
基于三维孔隙模型,选取能够代表孔隙模型结构特征的表征单元体(representative volume element,RVE)代替原模型.表征单元体是多孔介质渗流中重要的概念,是宏观尺度上保持介质性质基本稳定的最小单

图3 RVE尺寸与孔隙率变化曲线
Fig.3 RVE size and porosity variation curve
从

图4 RVE形心位置
Fig.4 RVE shape center locations
达西定律中渗透系数的计算式见
(1) |
式中:K为渗透系数,mm/s;v为渗流流速,mm/s;i为渗流水力坡降.
分析
将圆柱芯样孔隙模型作分割并剔除非有效孔隙后,对所得24个RVE孔隙模型进行特征参数的全面获取与筛选.综合分析比较有效孔隙率、有效孔隙比表面积、有效孔隙迂曲度、有效孔隙半径均值、有效孔隙三维分形维数、有效孔隙平均孔喉比、有效孔隙平均配位数、有效孔隙团簇离心率、有效孔隙等效喉道半径、有效孔隙等效喉道长度、有效孔隙费雷特直径和欧拉示性数等孔隙特征参数,根据“能够独立控制”、“对渗透性能影响大”和“相近代替取更优”原则,将以下3类参数筛去:(1)描述孔隙综合信息的孔隙三维分形维数(难以独立控制);(2)经数据分析对渗透系数影响较小的有效孔隙率、有效孔隙平均配位数、孔隙费雷特直径和欧拉示性数;(3)能够由更优参数替代的有效孔隙平均孔喉比(由喉道长度、喉道半径代替)和团簇离心率(由比表面积代替).利用SPSS软件得到剩余各孔隙特征参数间的线性相关程度,结果见
Parameter | Tortuosity | Specific surface | Mean pore radius | Equivalent throat length of pore | Equivalent throat radius of pore |
---|---|---|---|---|---|
Tortuosity | 1.000 | 0.145 | -0.455 | -0.467 | -0.439 |
Specific surface | 1.000 | -0.626 | -0.430 | -0.692 | |
Mean pore radius | 1.000 | 0.897 | 0.930 | ||
Equivalent throat length of pore | 1.000 | 0.875 | |||
Equivalent throat radius of pore | 1.000 |
由
基于上述分析,确定了影响沿程水头损失的有效孔隙迂曲度、有效孔隙半径均值,以及影响局部水头损失的有效孔隙比表面积(表征孔隙形状)、有效孔隙等效喉道半径与有效孔隙等效喉道长度作为本文透水路面砖孔隙特征参数.
有效孔隙迂曲度τ表征整个渗流路径沿程方向的信息,其计算式见
(2) |
式中:L为渗流试样长度,mm;Le为实际渗流路径长度,mm.
Le的计算方法为:将孔隙结构模型沿Z轴(高度方向)切割为多个微元平面,获取每个微元平面上的有效孔隙形心点,连接各形心点,计算各连接微元段长度的总和,即Le.
有效孔隙比表面积S表征有效孔隙形状信息,其计算式为:
(3) |
式中:Sp为有效孔隙表面积,μ
有效孔隙等效喉道半径req表征渗流路径中有效孔隙喉道的横截面信息,其计算使用最大内切球算法,即有效孔隙内最窄处等效体积最大球的半径.
有效孔隙等效喉道长度leq表征有效孔隙喉道的沿程方向信息.leq为喉道两端所连接2个有效孔隙等效体积球体心之间距离减去2个等效体积球半径之后的长度.
有效孔隙半径均值r指有效孔隙以渗径为圆柱母线,以体积不变为原则等效为圆柱型孔隙的半径.
通过计算共获得24组数据,计算孔隙特征参数与其均值的比值,得到如

图5 孔隙特征参数波动曲线
Fig.5 Variation range curves of pore characteristic parameters
利用质量法对经XCT扫描的透水砖圆柱芯样进行有效孔隙率测定,并与计算值对比,以验证软件仿真的可信度.试验步骤如下:
(1)测算无水干燥条件下透水砖圆柱芯样的直径d、高度h和质量md,精度分别为0.02 mm和0.01 g;
(2)将透水砖圆柱芯样置入饱和缸内,对饱和缸进行抽真空处理之后静置24 h,使透水砖圆柱芯样充分饱水;
(3)利用保鲜膜在饱和缸水面下将圆柱芯样紧密贴合包裹完全,之后小心取出,测量其质量m0;
(4)将含水保鲜膜快速小心取下,测定其质量mp;
(5)按下式计算实测透水砖圆柱芯样有效孔隙率.
(4) |
式中:ρw为水的密度.
试验共测量7次,剔除最大值与最小值后,取其均值作为最终的实测透水砖圆柱芯样有效孔隙率,为21.93%,其与仿真计算得到的有效孔隙率15.23%相差较大,主要原因有以下几点:
(1)有部分微孔隙的尺寸小于XCT扫描分辨率,在图像处理时被包含在圆柱芯样集料相中,无法被划分为孔隙相,导致软件计算值较真实值偏低.
(2)存在回孔效
(3)宏观有效孔隙率试验中,抽真空会导致半封闭孔隙被水充满,导致实测有效孔隙率较真实值偏大.
综上,有效孔隙率误差主要由XCT扫描及有效孔隙率实测试验的系统误差导致,因此单纯由软件本身引起的有效孔隙率计算误差可以忽略不计,即认为软件的计算结果是可靠的.
根据选定的对孔隙渗透性能影响较大的几个参数(τ、r、S、req、leq)的物理意义,建立绝对渗透率k与孔隙特征参数之间的关系,见
(5) |
绝对渗透率k与渗透系数K存在
(6) |
式中:γ为渗流流体的容重,N/
结合式(
(7) |
在Avizo软件中,对选取的24个RVE孔隙模型使用Axis Connectivity命令获得其连通孔隙模型,该命令可以将孔隙模型中不与Z轴上下表面连通的孔隙进行剔除,获得的连通孔隙模型见

图6 RVE连通孔隙模型
Fig.6 RVE connectivity pore model

图7 渗流模拟流速云图
Fig.7 Simulated flow velocity clouds for seepage
RVE | K/(mm∙ | RVE | K/(mm∙ |
---|---|---|---|
1 | 1.91 | 13 | 1.91 |
2 | 1.87 | 14 | 1.89 |
3 | 1.92 | 15 | 1.96 |
4 | 1.93 | 16 | 1.90 |
5 | 1.95 | 17 | 1.85 |
6 | 1.94 | 18 | 1.92 |
7 | 2.03 | 19 | 1.91 |
8 | 1.88 | 20 | 1.88 |
9 | 1.92 | 21 | 1.89 |
10 | 1.97 | 22 | 1.89 |
11 | 2.06 | 23 | 1.93 |
12 | 2.12 | 24 | 1.88 |
将由

图8 计算值与仿真值对比
Fig.8 Comparison of calculated values with simulation values
利用SPSS软件对渗透系数计算值与仿真值进行回归分析,得到方差分析结果和回归系数,分别见
Model | Sum of squares | Degrees of freedom | Mean square | F | Conspicuousness |
---|---|---|---|---|---|
Regression | 0.074 | 1 | 0.074 | 108.006 | <0.001 |
Residual | 0.015 | 22 | 0.001 | ||
Total | 0.089 | 23 |
Model | Unstandardized coefficient | Standardized coefficient | t | Conspicuous‑ness | |
---|---|---|---|---|---|
B | Standard error | Beta | |||
Constant | 0.008 | 0.185 | 0.043 | 0.966 | |
Calculated value | 0.997 | 0.096 | 0.911 | 10.393 | <0.001 |
从
经计算,渗透系数计算值均值与渗透系数仿真值均值的相对误差为0.16%,远小于5%,说明
(1)建立了陶瓷透水路面砖渗透系数与孔隙特征参数的定量关系模型,经过检验确定了该模型的可靠性与合理性.
(2)陶瓷透水路面砖渗透系数受有效孔隙迂曲度和有效孔隙半径均值的影响最大,受有效孔隙比表面积、有效孔隙等效喉道半径及有效孔隙等效喉道长度的影响次之.
参考文献
崔新壮, 张炯, 黄丹, 等. 暴雨作用下透水混凝土路面快速堵塞试验模拟[J]. 中国公路学报, 2016, 29(10):1‑11,19. [百度学术]
CUI Xinzhuang, ZHANG Jiong, HUANG Dan, et al. Experimental simulation of rapid clogging of pervious concrete pavement under effects of rainstorm[J]. China Journal of Highway and Transport, 2016, 29(10):1‑11,19. (in Chinese) [百度学术]
刘冠志, 侯泽彪, 廖霖, 等. 透水混凝土的细观孔隙结构分析[J]. 混凝土与水泥制品, 2022(7):40‑44. [百度学术]
LIU Guanzhi, HOU Zebiao, LIAO Lin, et al. Analysis of microscopic pore structure of pervious concrete[J]. China Concrete and Cement Products, 2022(7):40‑44. (in Chinese) [百度学术]
蒋正武, 孙振平, 王培铭. 若干因素对多孔透水混凝土性能的影响[J]. 建筑材料学报, 2005, 8(5):513‑519. [百度学术]
JIANG Zhengwu, SUN Zhenping, WANG Peiming. Effects of some factors on properties of porous pervious concrete[J]. Journal of Building Materials, 2005, 8(5):513‑519. (in Chinese) [百度学术]
郭远臣, 叶青, 向凯, 等. 净水功能型粉煤灰透水混凝土性能研究[J]. 建筑材料学报, 2023,26(9):988‑995. [百度学术]
GUO Yuanchen, YE Qing, XIANG Kai, et al. Performance of fly ash pervious concrete with water purification function[J]. Journal of Building Materials, 2023, 26(9):988‑995. (in Chinese) [百度学术]
丁崧, 陈潇, 夏飞跃, 等. 净水型赤泥-矿渣基地聚合物透水混凝土的研究[J]. 建筑材料学报, 2020, 23(1):48‑55. [百度学术]
DING Song, CHEN Xiao, XIA Feiyue, et al. Study on red mud‑slag based geopolymer pervious concrete with function of water purification [J]. Journal of Building Materials, 2020, 23(1):48‑55. (in Chinese) [百度学术]
贺拿, 陈宁生, 朱云华, 等. 泥石流源区砾石土分形特征及其与渗透系数关系试验研究[J]. 岩土力学, 2014, 35(9):2543‑2548. [百度学术]
HE Na, CHEN Ningsheng, ZHU Yunhua, et al. Experimental study on fractal feature and relationship between fractal dimension and permeability coefficient of gravelly soil in debris flow source area[J]. Rock and Soil Mechanics, 2014, 35(9):2543‑2548. (in Chinese) [百度学术]
ZHENG J J, ZHOU X Z. Analytical method for prediction of water permeability of cement paste[J]. ACI Materials Journal, 2008, 105(2):200‑206. [百度学术]
BAGEL L, ZIVICA V. Relationship between pore structure and permeability of hardened cement mortars:On the choice of effective pore structure parameter[J]. Cement and Concrete Research, 1997, 27(8):1225‑1235. [百度学术]
YANG Z J, REN W Y, SHARMA R, et al. In‑situ X‑ray computed tomography characterisation of 3D fracture evolution and image‑based numerical homogenisation of concrete[J]. Cement and Concrete Composites, 2017, 75:74‑83 [百度学术]
任文渊, 杨贞军, 黄宇劼. 基于X射线计算断层扫描图像的混凝土细观断裂模拟[J]. 水利学报, 2015, 46(4):452‑459. [百度学术]
REN Wenyuan, YANG Zhenjun, HUANG Yujie. Meso‑scale fracture modelling of concrete based on X‑ray computed tomography images[J]. Journal of Hydraulic Engineering, 2015, 46(4):452‑459. (in Chinese) [百度学术]
魏汉林, 任文渊, 张爱军, 等. 基于XCT扫描图像的混凝土粗骨料二维细观形态[J]. 科学技术与工程, 2021, 21(19):8183‑8189. [百度学术]
WEI Hanlin, REN Wenyuan, ZHANG Aijun, et al. Two‑dimensional mesoscale morphology of coarse aggregates in concrete based on XCT scanning images[J]. Science Technology and Engineering, 2021, 21(19):8183‑8189. (in Chinese) [百度学术]
叶强, 李喜龙, 梁莎, 等. 混凝土透水砖性能影响因素试验研究及孔结构定量计算[J]. 混凝土, 2022, 388(2):167‑171. [百度学术]
YE Qiang, LI Xilong, LIANG Sha, et al. Performance optimization of permeable concrete and quantitative calculation of pore structure[J]. Concrete, 2022, 388(2):167‑171. (in Chinese) [百度学术]
张振, 饶烽瑞, 叶观宝, 等. 基于X‑CT技术的气泡轻质土孔隙结构研究[J]. 建筑材料学报, 2020, 23(5):1104‑1112. [百度学术]
ZHANG Zhen, RAO Fengrui, YE Guanbao, et al. Investigation on viod structure of foamed light‑weight soil with X‑CT scanning technique[J]. Journal of Building Materials, 2020, 23(5):1104‑1112. (in Chinese) [百度学术]
高志涵, 陈波, 陈家林, 等. 基于X‑CT的泡沫混凝土孔隙结构与导热性能[J]. 建筑材料学报, 2023, 26(7):723‑730. [百度学术]
GAO Zhihan, CHEN Bo, CHEN Jialin, et al. Pore structure and thermal conductivity of foam concrete based on X‑CT[J]. Journal of Building Materials, 2023, 26(7):723‑730. (in Chinese) [百度学术]
郭玉柱, 陈徐东, 宁英杰, 等. 基于X‑CT的蒸养大掺量矿物掺合料砂浆孔结构[J]. 建筑材料学报, 2022, 25(9):885‑892. [百度学术]
GUO Yuzhu, CHEN Xudong, NING Yingjie, et al. Pore structure of steam cured high volume mineral admixture mortar based on X‑CT technology[J]. Journal of Building Materials, 2022, 25(9):885‑892. (in Chinese) [百度学术]
孙卫, 史成恩, 赵惊蛰, 等. X‑CT扫描成像技术在特低渗透储层微观孔隙结构及渗流机理研究中的应用——以西峰油田庄19井区长82储层为例[J]. 地质学报, 2006(5):775‑779,789‑790. [百度学术]
SUN Wei, SHI Cheng’en, ZHAO Jingzhe, et al. Application of X‑CT scanned image technique in the research of micro‑pore texture and percolation mechanism in ultra‑permeable oil field — Taking an example from the chang 82 formation in the Xifeng oil field[J]. Acta Geologica Sinica, 2006(5):775‑779,789‑790. (in Chinese) [百度学术]
杨鹏飞, 李显, 阙云, 等. RVE孔隙模型细观结构特征分析与对比[J]. 长江科学院院报, 2022, 39(8):99‑104. [百度学术]
YANG Pengfei, LI Xian, QUE Yun, et al. RVE‑based quantitative characterization of 3D mesostructure characteristics of granite residual soil with macropores[J]. Journal of Changjiang River Scientific Research Institute, 2022, 39(8):99‑104. (in Chinese) [百度学术]
延克军, 甄树聪. 流体力学[M]. 北京:中国建筑工业出版社, 2017:75‑77. [百度学术]
YAN Kejun, ZHEN Shucong. Fluid mechanics[M]. Beijing:China Architecture & Building Press, 2017:75‑77. (in Chinese) [百度学术]
SALZER M, SPETTL A, STENZEL O, et al. A two‑stage approach to the segmentation of FIB‑SEM images of highly porous materials[J]. Materials Characterization, 2012, 69:115‑126. [百度学术]