网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

陶瓷透水路面砖渗透系数与孔隙特征的定量关系模型  PDF

  • 王仪良 1
  • 谭建平 1
  • 李薇 1
  • 任文渊 1
  • 袁克阔 2
1. 西北农林科技大学 水利与建筑工程学院,陕西 杨凌 712100; 2. 西京学院 土木工程学院,陕西 西安 710123

中图分类号: TU997

最近更新:2024-02-29

DOI:10.3969/j.issn.1007-9629.2024.02.011

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
EN
目录contents

摘要

利用X射线计算机断层成像及图像三维重构技术获取了陶瓷透水路面砖的孔隙特征参数与渗透系数.基于沿程水头损失及局部水头损失影响因素,选取合适孔隙特征参数,并利用量纲分析法建立了渗透系数与孔隙特征参数的定量关系模型.结果表明:通过统计学检验方法对建立的定量关系模型进行检验,发现所建立模型是合理可靠的;渗透系数受有效孔隙迂曲度、有效孔隙半径均值的影响最大,受有效孔隙比表面积、有效孔隙等效喉道半径及有效孔隙等效喉道长度的影响较小.

随着海绵城市战略的提出,透水路面砖受到越来越多的关注,已被广泛应用于人行道、公园及停车场等区域.其优良的抗滑和降噪特性为人们出行提供了安全与舒适保障条件.为此国内外学者对透水路面材料的性能展开了大量研究.

崔新壮

1和刘冠志2分别研究了暴雨作用下淤塞以及骨料粒径的改变对于透水路面材料孔隙结构的影响.透水混凝土渗透系数的影响因素也受到了众多学者关3‑6.国外众多学者对孔隙特征参数与渗透系数的关系进行了大量研7‑8,提出了基于孔隙体积和孔隙半径预测渗透系数的经验公式,获得了较好效果.然而,这些研究主要基于孔隙半径和孔隙率等基本参数,而未涉及孔隙形状及分布等参数与渗透特性之间的关系.

本文拟利用X射线计算机断层成像(X‑ray computed tomography,XCT)与图像三维重构技术,获得孔隙三维模型,提取孔隙特征参数.基于沿程水头损失及局部水头损失影响机制,对孔隙特征参数进行相关性分析.然后,结合量纲分析法建立渗透系数与孔隙特征参数之间的定量关系模型,并基于孔隙模型的渗流数值模拟结果进行验证,希望为透水路面材料的应用及优化选用提供理论参考.

1 试验

1.1 陶瓷透水路面砖

以陶瓷透水路面砖为研究对象,规格为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所示.

图1  透水砖圆柱芯样(截取自实际试样)

Fig.1  Cylindrical core sample of pervious brick(actual specimens intercepted)

1.2 XCT扫描

XCT技术具有分辨率高、无损测试和三维可视化等优点,可用来无损捕捉混凝土试件中的真实空间分布和形态特

9‑10,还可为骨料与孔隙三维形态表征以及更接近真实形态随机骨料模型的生成提供参考和依11,在材料孔隙结构研究中已经得到了广泛应12‑16.

设备使用Multiscale‑Voxel 450 kV工业XCT扫描仪(CD‑300BX),扫描圆柱芯样获得图像模型大小为1 200 像素×1 200 像素,共采集到1 300幅图像,XCT扫描空间分辨率为43.63 μm.

2 孔隙特征参数

透水路面材料中的孔隙主要分为封闭孔隙、半封闭孔隙和连通孔隙3类.其中连通孔隙又称为有效孔隙,对透水路面材料渗透系数的影响最大.另外两类孔隙对于透水路面材料的渗透系数影响较小.故本文主要讨论陶瓷透水路面砖中有效孔隙及其特征参数,如无特别说明,下文中孔隙特征参数均指有效孔隙特征参数.

2.1 孔隙三维模型建立

将扫描得到的图像进行三维重构,并对其进行滤波处理,得到圆柱芯样的灰度切片图像.将灰度图像导入Avizo软件进行阈值分割,将孔隙相与集料相分离.最终得到透水砖圆柱芯样的三维孔隙模型,如图2所示.

图2  透水砖圆柱芯样的三维孔隙模型

Fig.2  3D pore model of pervious brick cylindrical core sample

基于三维孔隙模型,选取能够代表孔隙模型结构特征的表征单元体(representative volume element,RVE)代替原模型.表征单元体是多孔介质渗流中重要的概念,是宏观尺度上保持介质性质基本稳定的最小单

17.以结构有效孔隙率维持稳定为选取原则确定RVE尺寸.考虑到圆柱芯样原模型是分层建立的,且下层高度较高,故确定下层孔隙模型RVE尺寸的具体步骤如下:随机在三维孔隙模型中选取4个不同形心的RVE;逐步增加这4个RVE的尺寸,并计算其有效孔隙率的变化;当RVE的有效孔隙率不随尺寸大小变化时,即为目标RVE尺寸.分析结果见图3,其中软件计算得到的圆柱芯样整体结构有效孔隙率为15.23%.

图3  RVE尺寸与孔隙率变化曲线

Fig.3  RVE size and porosity variation curve

图3可以看出,当选取的RVE尺寸逐步增加到13.1 mm时,其有效孔隙率基本维持稳定.考虑设备算力,应尽可能缩小模型以减少计算量,故最终取下层孔隙结构模型的RVE尺寸为13.1 mm×13.1 mm×13.1 mm.同时,鉴于上层孔隙层高度较小,因此直接选取上层孔隙结构模型厚度6.5 mm作为上层RVE高度尺寸,故上下层贯通后RVE尺寸为13.1 mm×13.1 mm×19.6 mm.在圆柱芯样孔隙结构模型中随机选取24个RVE,考虑到其与圆柱芯样边界的位置关系,RVE形心位置如图4所示.

图4  RVE形心位置

Fig.4  RVE shape center locations

2.2 孔隙特征参数

2.2.1 孔隙特征参数的选取

达西定律中渗透系数的计算式见式(1).

K=vi (1)

式中:K为渗透系数,mm/s;v为渗流流速,mm/s;i为渗流水力坡降.

分析式(1)中流速v与水力坡降i的影响因素,以选取合适的孔隙特征参数.在给定进出口压力和渗流试样长度的前提下,其水力坡降i是一定的,而流速v的大小受到沿程水头损失与局部水头损失的影响.其中沿程水头损失为渗流过程中流体与孔隙内壁间的摩擦力于渗流路径长度内所做的负功,受到孔隙迂曲度(直接控制渗流路径长度)及孔隙半径(影响渗流速度梯度大小,从而间接影响黏滞力大小

18的影响;局部水头损失主要发生在孔隙截面突变处,受到孔隙喉道及截面突变程度(由孔隙形状表征)的影响.

将圆柱芯样孔隙模型作分割并剔除非有效孔隙后,对所得24个RVE孔隙模型进行特征参数的全面获取与筛选.综合分析比较有效孔隙率、有效孔隙比表面积、有效孔隙迂曲度、有效孔隙半径均值、有效孔隙三维分形维数、有效孔隙平均孔喉比、有效孔隙平均配位数、有效孔隙团簇离心率、有效孔隙等效喉道半径、有效孔隙等效喉道长度、有效孔隙费雷特直径和欧拉示性数等孔隙特征参数,根据“能够独立控制”、“对渗透性能影响大”和“相近代替取更优”原则,将以下3类参数筛去:(1)描述孔隙综合信息的孔隙三维分形维数(难以独立控制);(2)经数据分析对渗透系数影响较小的有效孔隙率、有效孔隙平均配位数、孔隙费雷特直径和欧拉示性数;(3)能够由更优参数替代的有效孔隙平均孔喉比(由喉道长度、喉道半径代替)和团簇离心率(由比表面积代替).利用SPSS软件得到剩余各孔隙特征参数间的线性相关程度,结果见表1,表中数字为皮尔逊相关性(R).

表1  各孔隙特征参数间的线性相关程度
Table 1  Degree of linear correlation between pore characteristic parameters
ParameterTortuositySpecific surfaceMean pore radiusEquivalent throat length of poreEquivalent 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

表1可见,孔隙半径均值与孔隙等效喉道长度、喉道半径存在较为显著的线性相关关系,这是因为喉道作为孔隙的一部分,必然会对孔隙半径均值产生影响.孔隙喉道作为孔隙渗流过程的关键制约段(类比于“木桶效应”中的“短板”),其重要性不言而喻.想要将孔隙结构特征进行完整表征,必然需要孔隙喉道的参与,且孔隙喉道作为孔隙中形状较为规则的部分,由沿程方向的喉道长度及横截面方向的喉道半径来表征.孔隙半径均值作为孔隙横截面方向信息的综合,还受到除喉道之外的其他因素影响,其在极大程度上影响渗流流速梯度,进而影响黏滞力大小,最终影响沿程水头损失,故必须选取其作为孔隙特征参数.

基于上述分析,确定了影响沿程水头损失的有效孔隙迂曲度、有效孔隙半径均值,以及影响局部水头损失的有效孔隙比表面积(表征孔隙形状)、有效孔隙等效喉道半径与有效孔隙等效喉道长度作为本文透水路面砖孔隙特征参数.

2.2.2 孔隙特征参数的获取

有效孔隙迂曲度τ表征整个渗流路径沿程方向的信息,其计算式见式(2).

τ=LeL (2)

式中:L为渗流试样长度,mm;Le为实际渗流路径长度,mm.

Le的计算方法为:将孔隙结构模型沿Z轴(高度方向)切割为多个微元平面,获取每个微元平面上的有效孔隙形心点,连接各形心点,计算各连接微元段长度的总和,即Le.

有效孔隙比表面积S表征有效孔隙形状信息,其计算式为:

S=SpVp (3)

式中:Sp为有效孔隙表面积,μm2Vp为有效孔隙体积,μm3.

有效孔隙等效喉道半径req表征渗流路径中有效孔隙喉道的横截面信息,其计算使用最大内切球算法,即有效孔隙内最窄处等效体积最大球的半径.

有效孔隙等效喉道长度leq表征有效孔隙喉道的沿程方向信息.leq为喉道两端所连接2个有效孔隙等效体积球体心之间距离减去2个等效体积球半径之后的长度.

有效孔隙半径均值r指有效孔隙以渗径为圆柱母线,以体积不变为原则等效为圆柱型孔隙的半径.

通过计算共获得24组数据,计算孔隙特征参数与其均值的比值,得到如图5所示的孔隙特征参数波动曲线.由图5可知,陶瓷透水砖内有效孔隙迂曲度、有效孔隙半径均值、有效孔隙等效喉道半径和有效孔隙等效喉道长度围绕各自均值在±10%范围内波动,有效孔隙比表面积则基本稳定,无明显变化.

图5  孔隙特征参数波动曲线

Fig.5  Variation range curves of pore characteristic parameters

2.3 宏观有效孔隙率试验

利用质量法对经XCT扫描的透水砖圆柱芯样进行有效孔隙率测定,并与计算值对比,以验证软件仿真的可信度.试验步骤如下:

(1)测算无水干燥条件下透水砖圆柱芯样的直径d、高度h和质量md,精度分别为0.02 mm和0.01 g;

(2)将透水砖圆柱芯样置入饱和缸内,对饱和缸进行抽真空处理之后静置24 h,使透水砖圆柱芯样充分饱水;

(3)利用保鲜膜在饱和缸水面下将圆柱芯样紧密贴合包裹完全,之后小心取出,测量其质量m0

(4)将含水保鲜膜快速小心取下,测定其质量mp

(5)按下式计算实测透水砖圆柱芯样有效孔隙率ϕ.

ϕ=4(m0-mp-md)/(ρwπd2h)×100% (4)

式中:ρw为水的密度.

试验共测量7次,剔除最大值与最小值后,取其均值作为最终的实测透水砖圆柱芯样有效孔隙率,为21.93%,其与仿真计算得到的有效孔隙率15.23%相差较大,主要原因有以下几点:

(1)有部分微孔隙的尺寸小于XCT扫描分辨率,在图像处理时被包含在圆柱芯样集料相中,无法被划分为孔隙相,导致软件计算值较真实值偏低.

(2)存在回孔效

19.因为XCT扫描无法精准识别深度,故一部分孔隙在XCT扫描时会被该扫描面下层集料“穿透”,造成一部分孔隙被误判成集料,导致软件计算的有效孔隙率较真实值偏低.

(3)宏观有效孔隙率试验中,抽真空会导致半封闭孔隙被水充满,导致实测有效孔隙率较真实值偏大.

综上,有效孔隙率误差主要由XCT扫描及有效孔隙率实测试验的系统误差导致,因此单纯由软件本身引起的有效孔隙率计算误差可以忽略不计,即认为软件的计算结果是可靠的.

3 理论公式

3.1 理论公式设计

根据选定的对孔隙渗透性能影响较大的几个参数(τrSreqleq)的物理意义,建立绝对渗透率k与孔隙特征参数之间的关系,见式(5).

k=9.964τ-0.518S0.363req-0.120r0.887leq-0.014 (5)

绝对渗透率k与渗透系数K存在式(6)所示关系.

K=k×γμ (6)

式中:γ为渗流流体的容重,N/m3μ为渗流流体的动力黏滞系数,Pa∙s.

结合式(5)、(6),得到拟定的渗透系数与孔隙特征参数之间的关系,如式(7)所示.其中,式(7)中各孔隙特征参数应该严格以μm或μm-1为单位.由式(7)中各孔隙特征参数的幂次大小可知:孔隙迂曲度τ及孔隙半径均值r对于透水材料渗透系数影响最大,其余参数影响次之.

K=9.964τ-0.518S0.363req-0.120r0.887leq-0.014γμ (7)

3.2 试验验证

在Avizo软件中,对选取的24个RVE孔隙模型使用Axis Connectivity命令获得其连通孔隙模型,该命令可以将孔隙模型中不与Z轴上下表面连通的孔隙进行剔除,获得的连通孔隙模型见图6.而后对其进行渗流模拟,在给定的入口压力、出口压力、渗流液体动力黏滞系数和主要渗流方向等参数条件下,对连通孔隙模型进行渗流模拟,渗流流速云图见图7.同时,计算得到各RVE孔隙模型的渗透系数仿真值,见表2.

图6  RVE连通孔隙模型

Fig.6  RVE connectivity pore model

图7  渗流模拟流速云图

Fig.7  Simulated flow velocity clouds for seepage

表2  RVE孔隙模型渗透系数仿真值
Table 2  Simulated permeability coefficient of RVE connectivity pore model
RVEK/(mm∙s-1)RVEK/(mm∙s-1)
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

将由式(7)得到的渗透系数公式计算值Kcv与数值软件得到的渗透系数仿真值Ksv进行对比,如图8所示,发现两者的线性相关关系较好.

图8  计算值与仿真值对比

Fig.8  Comparison of calculated values with simulation values

利用SPSS软件对渗透系数计算值与仿真值进行回归分析,得到方差分析结果和回归系数,分别见表3表4.其中:FF检验值;B为标准化前每个自变量在回归方程中的系数;Beta为标准化后自变量在回归方程中的系数;tt检验值.

表3  渗透系数方差分析结果
Table 3  Analysis of variance for permeability coefficient
ModelSum of squaresDegrees of freedomMean squareFConspicuousness
Regression 0.074 1 0.074 108.006 <0.001
Residual 0.015 22 0.001
Total 0.089 23
表4  回归系数
Table 4  Regression coefficients
ModelUnstandardized coefficientStandardized coefficienttConspicuousness
BStandard errorBeta
Constant 0.008 0.185 0.043 0.966
Calculated value 0.997 0.096 0.911 10.393 <0.001

表3表4可以看出,其F检验及t检验的显著性均小于0.001,远小于0.05,说明所建线性回归模型具有较为显著的可靠性.因此,可以认为式(7)模型是合理可靠的.

经计算,渗透系数计算值均值与渗透系数仿真值均值的相对误差为0.16%,远小于5%,说明式(7)能够在较大程度上反映孔隙结构特征对于透水路面材料渗透性能的影响规律.

4 结论

(1)建立了陶瓷透水路面砖渗透系数与孔隙特征参数的定量关系模型,经过检验确定了该模型的可靠性与合理性.

(2)陶瓷透水路面砖渗透系数受有效孔隙迂曲度和有效孔隙半径均值的影响最大,受有效孔隙比表面积、有效孔隙等效喉道半径及有效孔隙等效喉道长度的影响次之.

参考文献

1

崔新壮张炯黄丹. 暴雨作用下透水混凝土路面快速堵塞试验模拟[J]. 中国公路学报20162910):111,19. [百度学术] 

CUI XinzhuangZHANG JiongHUANG Danet al. Experimental simulation of rapid clogging of pervious concrete pavement under effects of rainstorm[J]. China Journal of Highway and Transport20162910):111,19. (in Chinese) [百度学术] 

2

刘冠志侯泽彪廖霖. 透水混凝土的细观孔隙结构分析[J]. 混凝土与水泥制品20227):4044. [百度学术] 

LIU GuanzhiHOU ZebiaoLIAO Linet al. Analysis of microscopic pore structure of pervious concrete[J]. China Concrete and Cement Products20227):4044. (in Chinese) [百度学术] 

3

蒋正武孙振平王培铭. 若干因素对多孔透水混凝土性能的影响[J]. 建筑材料学报200585):513519. [百度学术] 

JIANG ZhengwuSUN ZhenpingWANG Peiming. Effects of some factors on properties of porous pervious concrete[J]. Journal of Building Materials200585):513519. (in Chinese) [百度学术] 

4

郭远臣叶青向凯. 净水功能型粉煤灰透水混凝土性能研究[J]. 建筑材料学报2023269):988995. [百度学术] 

GUO YuanchenYE QingXIANG Kaiet al. Performance of fly ash pervious concrete with water purification function[J]. Journal of Building Materials2023269):988995. (in Chinese) [百度学术] 

5

丁崧陈潇夏飞跃. 净水型赤泥-矿渣基地聚合物透水混凝土的研究[J]. 建筑材料学报2020231):4855. [百度学术] 

DING SongCHEN XiaoXIA Feiyueet al. Study on red mudslag based geopolymer pervious concrete with function of water purification [J]. Journal of Building Materials2020231):4855. (in Chinese) [百度学术] 

6

贺拿陈宁生朱云华. 泥石流源区砾石土分形特征及其与渗透系数关系试验研究[J]. 岩土力学2014359):25432548. [百度学术] 

HE NaCHEN NingshengZHU Yunhuaet 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 Mechanics2014359):25432548. (in Chinese) [百度学术] 

7

ZHENG J JZHOU X Z. Analytical method for prediction of water permeability of cement paste[J]. ACI Materials Journal20081052):200206. [百度学术] 

8

BAGEL LZIVICA V. Relationship between pore structure and permeability of hardened cement mortars:On the choice of effective pore structure parameter[J]. Cement and Concrete Research1997278):12251235. [百度学术] 

9

YANG Z JREN W YSHARMA Ret al. Insitu Xray computed tomography characterisation of 3D fracture evolution and imagebased numerical homogenisation of concrete[J]. Cement and Concrete Composites2017757483 [百度学术] 

10

任文渊杨贞军黄宇劼. 基于X射线计算断层扫描图像的混凝土细观断裂模拟[J]. 水利学报2015464):452459. [百度学术] 

REN WenyuanYANG ZhenjunHUANG Yujie. Mesoscale fracture modelling of concrete based on Xray computed tomography images[J]. Journal of Hydraulic Engineering2015464):452459. (in Chinese) [百度学术] 

11

魏汉林任文渊张爱军. 基于XCT扫描图像的混凝土粗骨料二维细观形态[J]. 科学技术与工程20212119):81838189. [百度学术] 

WEI HanlinREN WenyuanZHANG Aijunet al. Twodimensional mesoscale morphology of coarse aggregates in concrete based on XCT scanning images[J]. Science Technology and Engineering20212119):81838189. (in Chinese) [百度学术] 

12

叶强李喜龙梁莎. 混凝土透水砖性能影响因素试验研究及孔结构定量计算[J]. 混凝土20223882):167171. [百度学术] 

YE QiangLI XilongLIANG Shaet al. Performance optimization of permeable concrete and quantitative calculation of pore structure[J]. Concrete20223882):167171. (in Chinese) [百度学术] 

13

张振饶烽瑞叶观宝. 基于XCT技术的气泡轻质土孔隙结构研究[J]. 建筑材料学报2020235):11041112. [百度学术] 

ZHANG ZhenRAO FengruiYE Guanbaoet al. Investigation on viod structure of foamed lightweight soil with XCT scanning technique[J]. Journal of Building Materials2020235):11041112. (in Chinese) [百度学术] 

14

高志涵陈波陈家林. 基于XCT的泡沫混凝土孔隙结构与导热性能[J]. 建筑材料学报2023267):723730. [百度学术] 

GAO ZhihanCHEN BoCHEN Jialinet al. Pore structure and thermal conductivity of foam concrete based on XCT[J]. Journal of Building Materials2023267):723730. (in Chinese) [百度学术] 

15

郭玉柱陈徐东宁英杰. 基于XCT的蒸养大掺量矿物掺合料砂浆孔结构[J]. 建筑材料学报2022259):885892. [百度学术] 

GUO YuzhuCHEN XudongNING Yingjieet al. Pore structure of steam cured high volume mineral admixture mortar based on XCT technology[J]. Journal of Building Materials2022259):885892. (in Chinese) [百度学术] 

16

孙卫史成恩赵惊蛰. XCT扫描成像技术在特低渗透储层微观孔隙结构及渗流机理研究中的应用——以西峰油田庄19井区长82储层为例[J]. 地质学报20065):775779789790. [百度学术] 

SUN WeiSHI Cheng’enZHAO Jingzheet al. Application of XCT scanned image technique in the research of micropore texture and percolation mechanism in ultrapermeable oil field — Taking an example from the chang 82 formation in the Xifeng oil field[J]. Acta Geologica Sinica20065):775779789790. (in Chinese) [百度学术] 

17

杨鹏飞李显阙云. RVE孔隙模型细观结构特征分析与对比[J]. 长江科学院院报2022398):99104. [百度学术] 

YANG PengfeiLI XianQUE Yunet al. RVEbased quantitative characterization of 3D mesostructure characteristics of granite residual soil with macropores[J]. Journal of Changjiang River Scientific Research Institute2022398):99104. (in Chinese) [百度学术] 

18

延克军甄树聪. 流体力学[M]. 北京中国建筑工业出版社20177577. [百度学术] 

YAN KejunZHEN Shucong. Fluid mechanics[M]. BeijingChina Architecture & Building Press20177577. (in Chinese) [百度学术] 

19

SALZER MSPETTL ASTENZEL Oet al. A twostage approach to the segmentation of FIBSEM images of highly porous materials[J]. Materials Characterization201269115126. [百度学术]