Numerical modeling of influences of discontinuous centrifugation on contaminant migration in cutoff wall
-
摘要: 土-膨润土竖向阻隔墙广泛应用于场地管控和修复工程。正确评估阻隔墙的服役寿命，成为阻隔墙设计的关键。使用现场试验评估阻隔墙的长期服役性能往往需要数十年，而使用超重力模型试验能够将试验周期大大缩短。在高离心加速度下，土工离心机无法长时间工作，为解决该问题本文提出采用非连续转机的方法来模拟污染物长时间的运移。在非连续转机情况下，污染物的运移发生了变化，为正确评估阻隔墙的服役寿命，建立了一个数值模型来模拟污染物在离心加速度周期性变化下的运移行为。研究结果表明：①采用有效转机时长作为非连续转机情况下污染物运移的时间，则阻隔墙的服役寿命将低于连续转机。②离心机停机时，阻隔墙内污染物运移速度降低有限，污染物运移以扩散为主。③阻隔墙的服役寿命随离心机转机间隔的增大而减少。Abstract: The soil-bentonite vertical cutoff walls are widely used in the remediation and control of contaminated sites. The correct evaluation of the service life of the cutoff walls is the key to their design. It often takes decades to evaluate the long-term performance of the cutoff walls by field tests, and the centrifugal model tests can obviously shorten the test period. In the case of high centrifugal acceleration, the geotechnical centrifuge cannot work for a long time. To solve this problem, a method of discontinuous centrifugation is proposed to simulate the long-term migration of contaminants. In the case of discontinuous centrifugation, the migration behavior of contaminants is changed. To correctly evaluate the service life of the cutoff walls, a numerical model is established to simulate the migration of the contaminants under the periodic change of centrifugal acceleration. The results show that: (1) The service life of the cutoff walls in the case of discontinuous centrifugation will be lower than that of continuous centrifugation if the effective centrifugation time is used as the time of contaminant transport. (2) When the centrifuge shuts down, the decreased migration velocity of contaminants in the cutoff wall is limited, and diffusion has become the main migration mode of contaminants. (3) The service life of the cutoff walls decreases with the increase of the work interval of the centrifuge.
-
Keywords:
- cutoff wall /
- centrifugal model test /
- discontinuous centrifugation /
- numerical model /
- service life
-
0. 引言
北疆供水渠道在进行渠道设计时，不仅仅要考虑经济最优和保证流量最优的问题，还需要考虑渠道的冬季冻胀破坏。梯形渠道占中国输水渠道的比例较大，应用较广。因此，梯形渠道的冻胀问题受到了北方寒区广泛关注。王正中等[1-2]对梯形渠道冻胀模型进行了理论分析，将梯形渠道考虑为在切向冻结力、法向冻胀力的以及衬砌板约束下的简支梁，王正中等[3]、刘旭东等[4]、李爽等[5]用数值模拟的方法模拟梯形断面形式下冻胀的规律，得到了温度场，位移场沿着渠道分布情况。
本文主要研究不同梯形断面情况下，以北疆某供水工程总干渠退水渠段梯形渠道为例，对该渠道用comsol进行数值模拟，对不同的断面渠道进行参数化模拟，得到最佳参数范围，并结合水力最佳断面，得到双优断面形式，为寒区渠道设计提供参考。
1. 力学模型的建立
1.1 基本假设
渠基土冻结时，土体、水和冰之间的相互作用的微观结构及动态过程相当复杂，目前很难准确模拟整个冻胀过程。为了便于分析，对其进行适当的简化，以便抓住影响冻结过程中及冻胀变形的主要特征，主要假设如下[6]：①根据现场及室内试验研究，假设冻土是均匀连续各向同性体；②尽管土的冻胀与其温度、水分、土质密切相关，对具体工程当水分及土壤条件确定时土体最终冻胀主要取决于温度，将水分迁移对冻土的体积影响，以线膨胀系数表达；③根据试验研究假定相变温度在同一种土中和同种外力条件下为常值，即暂取相变温度为0℃[7]；④由于渠道为细长结构，不考虑长方向上土颗粒对温度的影响，选取平面应变问题进行模拟。
1.2 热传导方程
根据以上假设，将温度热传递视为二维瞬态热传递的过程，建立二维热传导控制方程：
(λx∂2T∂x2)+(λy∂2T∂y2)=ρc∂T∂t。 (1) 式中T为温度；λx，λy分别是冻土沿x，y方向的导热系数；ρ为土体的密度；c为土体的比热容；t为时间。
1.3 本构方程
冻土属于冷胀热缩材料，冻土在冻结过程中水冻结成冰，除了原位水冻结体积膨胀，还有从未冻结区向冻结区迁移的水分冻结成冰，本文将原位水以及迁移水冻结成冰的体积膨胀以关于温度T的函数的线膨胀系数表示。因此，土体的本构方程可以表示为
ϵx=σxE−μσyE+α(T−Tref) ，ϵy=σyE−μσxE+α(T−Tref) ，γxy=2(1+μ)Eτxy 。} (2) 式中εx，εy为正应变；γxy为剪应变；σx、σy为正应力；τxy为剪应力；E为弹性模量；μ为泊松比，α为混凝土或者冻土自由冻胀时的热膨胀系数；T为计算温度；Tref为参考温度。
2. 计算模型
2.1 原型概况
本文选择北疆地区某输水渠道退水渠段的梯形渠道（图 1）为研究对象，原型渠道的基本情况如表 1，2所示，并对该模型渠道进行了冻胀模拟分析计算，因原型渠道各部位坡向，土质，水分不同，在模型边界设置时对阴阳破和渠堤分别对温度和土体冻胀率赋值。
表 1 渠道各部位的表面温度以及冻结期Table 1. Surface temperatures and freezing periods of various parts of canal部位 月平均表面温度/(℃) 冻结期、（月-日） 12月 1月 2月 阴坡 -14.92 -18.85 -10.72 11-27—02-27 渠底 -14.56 -16.22 -9.15 11-27—02-26 阳坡 -12.55 -14.75 -10.54 11-27—02-27 表 2 原型渠道基本情况Table 2. Basic information of prototype canal部位 渠床
土质基土干密度/(g·cm-3) 冻深h/cm 冻胀量Δh/cm 冻胀率η/% 阴坡 砂砾石 1.80 171 5.0 2.92 渠底 159 4.4 2.77 阳坡 146 3.7 2.53 2.2 有限元模型
有限元计算模型为：①原型梯形断面渠道冻胀数值模拟；②梯形断面渠道冻胀“参数化”断面渠道数值模拟。
有限元模型如图 2所示，选取各表面温度接近原型渠道，渠道两边设置为热绝缘边界，下边界为固定温度边界取10℃，阴坡边界温度取冻结期平均值-14.8℃，渠底温度取-13.3℃，阳坡温度为-12.6℃。几何模型如图 2所示，渠底长2.0 m，渠坡的横向长度为3.75 m，坡比为1.5，基础底向下取2.5 m，左右边界取0.75 m。模型将衬砌板与冻土作为一个整体进行数值模拟，力学的边界设置时，左右边界设置为辊支撑，底部设置指定位移0。
2.3 参数取值
参考文献[5]中线膨胀系数按照η/Tmin取值，η为冻胀率，Tmin为相应部位月平均表面温度的最小值。由式（2），将T-Tref看成1℃，那么可以将线膨胀系数α考虑跟温度Tmin相关的函数，数值上与冻胀率相等。未冻土与混凝土的导热系数参考文献[3]，冻土的导热系数与土体的未冻水含量相关，冻土融土导热系数比与未冻水含量相关[9]，模型渠道的未冻水含量选取19.4%，因此冻土的导热系数为未冻土导热系数的0.9倍。
3. 计算结果分析
3.1 温度场
图 3可以看出：接近于渠底的温度分布接近于一组平行的直线，0℃等温线以上，由于未冻水冻结导致剧烈相变，该区域导热系数小，因此等温线较密集，温度梯度比较小；0℃等温线以下，导热系数较大，等温线较稀疏。阴坡冻深为162 cm，渠底冻深为146 cm，阳坡的冻深为141 cm。与表 3模型试验较符合，最大相对误差为8.2%。
表 3 材料热力学参数Table 3. Thermodynamic parameters of materials介质 弹性模量E/Pa 泊松比ν 导热系数λ/(W·m-1·℃) 线膨胀系数α/℃-1 混凝土 2.4×1010 0.167 1.580 1.1×10-5 冻土 4.6×107 0.330 1.188[9] 阴坡2.92%、渠底2.77%、阳坡2.53% 未冻土 1.5×107 0.375 1.320 0 3.2 位移场
如图 4所示的是模拟得到的结果与试验结果进行对比，模型试验结果冻胀量边坡较大，渠底较小。模型试验的冻胀量整体比模拟冻胀量大，因为模型试验所用的衬砌板受到边坡的约束较大，且衬砌板厚度小于模拟设定厚度。模拟得到渠底冻胀量比较平缓而试验得到的渠底中间部位冻胀量大而两边小，模拟设定的衬砌板是整体现浇型板，在冻胀力的作用下容易产生整体变形；而模型试验的衬砌板分布较为离散。
3.3 “参数化”位移场模拟结果
根据梯形断面的设计规范要求，应满足：①渠道边坡系数应不小于允许最小边坡系数，应不大于允许最大边坡系数，1≤m≤2。②渠道的宽深比应满足规范给定的要求，0.8 < b/h < 3.5。
根据要求，“参数化”分析模拟了宽深比β一定时，不同边坡比的渠道断面冻胀量情况；和边坡比一定时，不同宽深比的渠道断面冻胀量情况。
由图 5，宽深比一定时，边坡比越小，渠坡的法向冻胀量越大；m=2时，左侧阴坡法向冻胀量最大为5.2 cm，阳坡为1.7 cm。m=1时，阴坡法向冻胀量最大值为3.3 cm。由图 6，边坡比一定时，宽深比对法向冻胀量影响不大。
![]()
图 5 不同坡比法向冻胀量比较Figure 5. Comparison of normal frost-heave quantity under different slopes![]()
图 6 不同宽深比法向冻胀量比较Figure 6. Comparison of normal frost-heave quantity under different breadth depth ratios4. 水力最佳断面
从设计角度考虑水力最佳断面通常是指过水流量一定，所要求的过水断面面积最小，也就是所需要的材料最少，施工且最方便[12]。或者是过水断面面积一定，所通过的流量最大。
考虑最小冻胀量，所以采用边坡比为2的冻胀断面。水力学中考虑流量一定，过水断面面积最小可采用以下公式计算宽深比[13]：
β=2（√1+m2−1)。 (3) 将相关参数代入式（3），可得到考虑水力冻胀最优宽深比为2.472。
5. 结语
本文利用comsol有限元软件，对北疆供水渠道以及渠基土冻胀破坏规律进行有限元分析，与试验数据对比，基本与试验数据相似。对不同坡比以及不同宽深比的渠道进行有限元模拟分析，发现坡比影响冻胀量的大小，坡比越大，渠道冻胀量越小。理论分析受力情况易得坡比越大越能释放变形，符合实际。再结合水力最佳断面的公式，可以得到最佳宽深比。最终的断面形式能有效的防止冻胀，且满足水力断面要求。为北疆地区渠道防冻胀断面设计提供了理论支撑。
-
![]()
图 1 离心模型中污染物二维运移示意图
Figure 1. Schematic diagram of two-dimensional transport of contaminants in centrifuge model
![]()
图 3 连续与非连续转机计算结果对比(te=122.5 h）
Figure 3. Comparison of calculated results between continuous and discontinuous centrifugation (te=122.5 h)
![]()
图 4 连续转机与非连续转机污染物浓度沿深度的变化
Figure 4. Variation of concentration of contaminants along depth under continuous and discontinuous centrifugations
![]()
图 5 未改性阻隔墙-A连续转机与非连续转机污染物通量随时间的变化
Figure 5. Variation of fluxes of contaminants of unmodified cutoff wall A with time under continuous and discontinuous centrifugations
![]()
图 6 改性阻隔墙-B连续转机与非连续转机污染物通量随时间的变化
Figure 6. Variation of fluxes of contaminants of modified cutoff wall B with time under continuous and discontinuous centrifugations
![]()
图 7 污染物最大相对浓度随时间的变化
Figure 7. Variation of maximum relative concentration of contaminants with time
![]()
图 8 不同转机间隔下污染物最大相对浓度随时间的变化
Figure 8. Variation of the maximum relative concentration of contaminants with time at different work intervals
表 1 模型计算输入参数[8]
Table 1 Input parameters for proposed model
位置 层厚L/cm 层高H/cm 有效扩散系数D∗/
(m2·s-1)干密度ρd/(g·cm-3) 分配系数Kd/(mg·L-1) 孔隙率n 弥散度αL/m 渗透系数k/(m·s-1) 浅部 中部 底部 浅部 中部 底部 浅部 中部 底部 浅部 中部 底部 上游含水层 7 23.6 4.82×10-10 1.51 1.51 1.51 0 0 0 0.43 0.43 0.43 6 1×10-3 1×10-3 1×10-3 下游含水层 8.5 21.6 4.82×10-10 1.51 1.51 1.51 0 0 0 0.43 0.43 0.43 6 1×10-3 1×10-3 1×10-3 未改性阻隔墙 4 28 4.5×10-10 1.486 1.471 1.48 1.50 2.00 3.01 0.443 0.449 0.446 0.0028 9×10-8 9.74×10-8 9.67×10-8 改性阻隔墙 4 27 4.5×10-10 1.484 1.491 1.487 3.00 5.70 8.66 0.445 0.443 0.444 0.0028 6.2×10-8 8.03×10-8 6.89×10-8 透水石 2 30 4.16×10-10 1.66 1.66 1.66 0 0 0 0.38 0.38 0.38 0.5 1×10-3 1×10-3 1×10-3 
下载: 导出CSV
-
[1] YEO S S, SHACKELFORD C D, EVANS J C. Consolidation and hydraulic conductivity of nine model soil-bentonite backfills[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2005, 131(10): 1189–1198. doi: 10.1061/(ASCE)1090-0241(2005)131:10(1189)
[2] EVANS J C. Geotechnics of hazardous waste control systems[M]//Foundation Engineering Handbook. Boston, MA: Springer US, 1991: 750–777.
[3] D'APPOLONIA D J. Soil-bentonite slurry trench cutoffs[J]. Journal of the Geotechnical Engineering Division, 1980, 106(4): 399–417. doi: 10.1061/AJGEB6.0000945
[4] BAXTER D Y. Mechanical behavior of soil-bentonite cutoff walls[D]. Virginia Tech, 2000.
[5] TAYLOR R N. Geotechnical Centrifuge Technology[M]. London: Blackie Academic & Professional, 1995.
[6] ARULANANDAN K, THOMPSON P Y, KUTTER B L, et al. Centrifuge modeling of transport processes for pollutants in soils[J]. Journal of Geotechnical Engineering, 1988, 114(2): 185–205. doi: 10.1061/(ASCE)0733-9410(1988)114:2(185)
[7] LO I M C, HU J Z & L. Centrifuge modeling of cadmium migration in saturated and unsaturated soils[J]. Soil and Sediment Contamination: An International Journal, 2005.
[8] ZHAN L T, YOU Y Q, ZHAO R, et al. Centrifuge modelling of lead retardation in soil–bentonite cut-off walls. International Journal of Physical Modelling in Geotechnics, 2022: 1–14. doi: 10.1061/(ASCE)1090-0241(2005)131:4(418)
[9] YEO S S, SHACKELFORD C D, EVANS J C. Membrane behavior of model soil–bentonite backfills[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2005, 131(4): 418–429. doi: 10.1061/(ASCE)1090-0241(2005)131:4(418)
[10] HENNING J T, EVANS J C, SHACKELFORD C D. Membrane behavior of two backfills from field-constructed soil-bentonite cutoff walls[J]. Journal of geotechnical and geoenvironmental engineering, 2006, 132(10): 1243–1249. doi: 10.1061/(ASCE)1090-0241(2006)132:10(1243)
[11] EVANS J C, SHACKELFORD C D, YEO S S, et al. Membrane behavior of soil-bentonite slurry-trench cutoff walls[J]. Soil & Sediment Contamination, 2008, 17(4): 316–322.
[12] ROTH K, JURY W A, FLÜHLER H, et al. Transport of chloride through an unsaturated field soil[J]. Water Resources Research, 1991, 27(10): 2533–2541. doi: 10.1029/91WR01771
[13] BEAR J, VERRUIJT A. Modeling Groundwater Flow and Pollution[M]. Springer: Science & Business Media, 1987.
[14] SMITH D W. One‐dimensional contaminant transport through a deforming porous medium: theory and a solution for a quasi‐steady‐state problem. International journal for numerical and analytical methods in geomechanics, 2000, 24(8): 693–722. doi: 10.1002/1096-9853(200007)24:8<693::AID-NAG91>3.0.CO;2-E
[15] DING X H, FENG S J, ZHENG Q T, et al. A two-dimensional analytical model for organic contaminants transport in a transition layer-cutoff wall-aquifer system[J]. Computers and Geotechnics, 2020, 128: 103816. doi: 10.1016/j.compgeo.2020.103816
[16] WHITMORE M K B, TROTT D W, PEERCY B E, et al. COMSOL modeling of groundwater flow and contaminant transport in two-dimensional geometries with heterogeneities[J]. UMBC Geography and Environmental Systems Department, 2011.
[17] FAISAL A A H, AL-WAKEL S F A, ASSI H A, et al. Waterworks sludge-filter sand permeable reactive barrier for removal of toxic lead ions from contaminated groundwater[J]. Journal of Water Process Engineering, 2020, 33: 101112. doi: 10.1016/j.jwpe.2019.101112
[18] 丁祥鸿, 冯世进. 固结作用下非均质隔离墙中污染物二维运移规律研究[J]. 岩土工程学报, 2022, 44(3): 584–590. doi: 10.11779/CJGE202203021 DING Xiang-hong, FENG Shi-jin. Two-dimensional migration of contaminants in non-homogeneous cutoff wall considering consolidation[J]. Chinese Journal of Geotechnical Engineering, 2022, 44(3): 584–590. (in Chinese) doi: 10.11779/CJGE202203021
计量
- 文章访问数: 113
- HTML全文浏览量: 20
- PDF下载量: 18
