岩土力学 ›› 2022, Vol. 43 ›› Issue (2): 528-538.doi: 10.16285/j.rsm.2021.1255

• 数值分析 • 上一篇    下一篇

污染场地六价铬迁移转化机制与数值模拟研究

贺勇1, 2,胡广1, 2,张召1, 2,娄伟3,邹艳红1, 2,李星3,张可能1, 2   

  1. 1. 中南大学 有色金属成矿预测与地质环境监测教育部重点实验室,湖南 长沙 410083; 2. 中南大学 地球科学与信息物理学院,湖南 长沙 410083;3. 湖南省和清环境科技有限公司,湖南 长沙 410221
  • 收稿日期:2021-08-08 修回日期:2021-12-01 出版日期:2022-02-11 发布日期:2022-02-22
  • 作者简介:贺勇,男,1987年生,博士,副教授,博士生导师,主要从事环境工程地质和非饱和土力学方面的研究。
  • 基金资助:
    国家重点研发计划(No. 2019YFC1805905);国家自然科学基金(No. 42072318,No. 41972282,No. 41807253);环境地球化学国家重点实验室开放课题(No. SKLEG2021208);湖南省自然科学基金(No. 2019JJ50763);中南大学中央高校基本科研业务费专项资金(No. 2021zzts0254)。

Numerical simulation on the migration and transformation mechanism of hexavalent chromium in contaminated site

HE Yong1, 2, HU Guang1, 2, ZHANG Zhao1, 2, LOU Wei3, ZOU Yan-hong1, 2, LI Xing3, ZHANG Ke-neng1, 2   

  1. 1. Key Laboratory of Metallogenic Prediction of Nonferrous Metals and Geological Environment Monitoring, Ministry of Education, Central South University, Changsha, Hunan 410083, China; 2. School of Geosciences and Info-Physics, Central South University, Changsha, Hunan 410083, China; 3. Hunan HIKEE Environmental Technology Co., Ltd., Changsha, Hunan 410221, China
  • Received:2021-08-08 Revised:2021-12-01 Online:2022-02-11 Published:2022-02-22
  • Supported by:
    This work was supported by the National Key Research and Development Program of China (2019YFC1805905), the National Natural Science Foundation of China (42072318, 41972282, 41807253), the Opening Fund of the State Key Laboratory of Environmental Geochemistry (SKLEG2021208), the Natural Science Foundation of Hunan Province (2019JJ50763) and the Fundamental Research Funds for the Central Universities of Central South University (2021zzts0254).

摘要: 随着全球工业化迅猛发展,土壤和地下水六价铬污染日益严重。基于某铁合金厂铬渣场地现场调查与采样分析,开展铬渣场地土样吸附、渗透和弥散试验,研究六价铬在粉质黏土土样中的吸附特性和迁移规律,建立考虑对流-弥散-吸附的六价铬迁移三维动力学模型,结合数值软件获取污染源位于场地上、下游时地下水中六价铬迁移分布特征,并揭示弥散度?和分配系数 对六价铬时空分布的影响。试验结果表明,粉质黏土对六价铬吸附符合Langmuir等温吸附模型,最大吸附量为466.6 mg/kg;蒸馏水和160 mg/L 六价铬溶液入渗下粉质黏土渗透系数约为6.5×10–7~6.7×10–7 cm/s,1 000 mg/L六价铬溶液的渗透系数增大至4.4×10–6 cm/s;粉质黏土水动力弥散系数D为1.4×10–4 m2/d,计算得到阻滞因子 为4.2~10。数值模拟结果表明,场地下游受到六价铬污染时,即使不考虑分子扩散作用,上游仍有被污染的风险,污染程度取决于含水层的弥散度;考虑含水层对六价铬吸附时,土体分配系数越大,六价铬污染羽分布范围越小,在预测地下水中六价铬浓度分布时应重点考虑六价铬吸附等转化过程。

关键词: 污染场地, 六价铬, 迁移转化, 弥散, 数值模拟

Abstract: With the rapid development of global industrialization, the pollution of hexavalent chromium (Cr(Ⅵ)) in soil and groundwater has become increasingly serious. Field investigation and laboratory tests were carried out for the soil polluted by the chromium slag of a ferroalloy plant. The adsorption, infiltration and dispersion experiments were conducted to study the adsorption characteristics and migration mechanism of Cr(Ⅵ) in silty clay. A three-dimensional kinetic mathematical model of Cr(Ⅵ) migration considering convection-dispersion-adsorption was established. The migration and distribution characteristics of Cr(Ⅵ) in groundwater with the pollution source located upstream or downstream of the contaminated site were obtained using the numerical approach. Meanwhile, the effects of dispersity (?) and distribution coefficient ( ) on the spatial and temporal distribution of Cr(Ⅵ) were revealed. The experimental results show that the Langmuir isotherm model well fits the adsorption data of silty clay. The maximum adsorption capacity of silty clay for Cr(Ⅵ) was 466.6 mg/kg. The hydraulic conductivity of silty clay under the infiltration of distilled water and 160 mg/L Cr(Ⅵ) solution was 6.5×10–7–6.7×10–7 cm/s, while it increased to 4.4×10–6 cm/s under infiltration of Cr(Ⅵ) solution with a concentration of 1 000 mg/L. The hydrodynamic dispersion coefficient (D) of silty clay was 1.4×10–4 m2/d. The value of the retardation factor ( ) was found to be 4.2–10. The results of the numerical simulation indicated that when the downstream was contaminated by Cr(Ⅵ), there was still a risk of pollution in the upstream even if molecular diffusion was not considered. The degree of pollution depended on the dispersity of the aquifer. Considering the adsorption of Cr(Ⅵ) by the aquifer, the higher the soil distribution coefficient, the smaller was the distribution range of the Cr(Ⅵ) pollution plume. Therefore, the transformation processes such as Cr(Ⅵ) adsorption should be focused on when predicting the distribution of Cr(Ⅵ) in contaminated sites.

Key words: contaminated site, hexavalent chromium Cr(Ⅵ), migration and transformation, dispersion, numerical simulation

中图分类号: 

  • X 53
[1] 吕茂淋, 朱珍德, 周露明, 葛鑫梁, . 基于相场法的预制双裂隙岩体水力压裂扩展数值模拟研究[J]. 岩土力学, 2024, 45(6): 1850-1862.
[2] 马登辉, 韩迅, 蔡正银, 关云飞, . 静压桩的桩侧土压力分布规律数值分析[J]. 岩土力学, 2024, 45(6): 1863-1872.
[3] 陈磊, 张强, 贾朝军, 雷明锋, 黄娟, 胡晶, . 强降雨对库岸堆积体边坡稳定性影响的离心模型试验和数值模拟研究[J]. 岩土力学, 2024, 45(5): 1423-1434.
[4] 周慧颖, 李树忱, 王曼灵, 袁超, 张俊艇, 冯健鹏, . 气孔杏仁状玄武岩随机三维模型构建与 单轴压缩力学行为研究[J]. 岩土力学, 2024, 45(4): 1214-1232.
[5] 文磊, 刘 钟, 马晓华, 张振. 粉土地基中劲性复合桩抗压承载特性与荷载传递机制研究[J]. 岩土力学, 2024, 45(2): 511-524.
[6] 刘新荣, 王浩, 郭雪岩, 罗新飏, 周小涵, 许彬, . 考虑消落带岩体劣化影响的典型危岩岸坡稳定性研究[J]. 岩土力学, 2024, 45(2): 563-576.
[7] 朱寅斌, 李长冬, 周佳庆, 项林语, 姜茜慧, 朱文宇, . 考虑基质渗透性的粗糙单裂隙非达西流动特性研究[J]. 岩土力学, 2024, 45(2): 601-611.
[8] 华涛, 申林方, 王志良, 李泽, 徐则民. 基于近场动力学的岩石水力压裂数值模拟及裂隙网络定量分析[J]. 岩土力学, 2024, 45(2): 612-622.
[9] 黄生根, 张义, 霍昊, 陈常青. 软土地区深基坑支护工程格构柱变形规律研究[J]. 岩土力学, 2023, 44(S1): 533-538.
[10] 王凯, 付强, 徐超, 艾子博, 李丹, 王磊, 舒龙勇, . 原生煤岩组合体界面力学效应数值模拟研究[J]. 岩土力学, 2023, 44(S1): 623-633.
[11] 乔亚飞, 闫凯, 赵腾腾, 丁文其, . 软土地区超深圆形竖井的坑底隆起特性与机制[J]. 岩土力学, 2023, 44(9): 2707-2716.
[12] 张坤勇, 张梦, 孙斌, 李福东, 简永洲, . 考虑时空效应的软土狭长型深基坑地连墙变形计算方法[J]. 岩土力学, 2023, 44(8): 2389-2399.
[13] 尹鑫晟, 舒营, 梁禄钜, 张世民, . 考虑渗流的饱和粉土地层盾构开挖面稳定分析[J]. 岩土力学, 2023, 44(7): 2005-2016.
[14] 季雨坤, 王钦科, 赵国良, 张健, 马建林, . 斜坡上嵌岩抗拔桩竖向承载变形特性模型试验及数值模拟[J]. 岩土力学, 2023, 44(6): 1604-1614.
[15] 孙彦晓, 刘松玉, 童立元, 王峻, 崔佳, 李世龙, 李敏, . 长江漫滩区明挖隧道基坑降承压水优化分析[J]. 岩土力学, 2023, 44(6): 1800-1810.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 徐金明,羌培,张鹏飞. 粉质黏土图像的纹理特征分析[J]. , 2009, 30(10): 2903 -2907 .
[2] 梁桂兰,徐卫亚,谈小龙. 基于熵权的可拓理论在岩体质量评价中的应用[J]. , 2010, 31(2): 535 -540 .
[3] 于琳琳,徐学燕,邱明国,闫自利,李鹏飞. 冻融作用对饱和粉质黏土抗剪性能的影响[J]. , 2010, 31(8): 2448 -2452 .
[4] 王 伟,刘必灯,周正华,王玉石,赵纪生. 刚度和阻尼频率相关的等效线性化方法[J]. , 2010, 31(12): 3928 -3933 .
[5] 王海波,徐 明,宋二祥. 基于硬化土模型的小应变本构模型研究[J]. , 2011, 32(1): 39 -43 .
[6] 曹光栩,宋二祥,徐 明. 山区机场高填方地基工后沉降变形简化算法[J]. , 2011, 32(S1): 1 -5 .
[7] 刘华丽 ,朱大勇 ,钱七虎 ,李宏伟. 边坡三维端部效应分析[J]. , 2011, 32(6): 1905 -1909 .
[8] 刘年平 ,王宏图 ,袁志刚 ,刘竟成. 砂土液化预测的Fisher判别模型及应用[J]. , 2012, 33(2): 554 -557 .
[9] 王卫东 ,李永辉 ,吴江斌 . 超长灌注桩桩-土界面剪切模型及其有限元模拟[J]. , 2012, 33(12): 3818 -3824 .
[10] 卫振海 ,王梦恕 ,张顶立 . 土结构强度模型研究[J]. , 2013, 34(1): 40 -46 .