岩土力学 ›› 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]. 岩土力学, 2025, 46(S1): 507-518.
[2] 张奇, 王驹, 刘江峰, 曹胜飞, 谢敬礼, 成建峰, . 热-水-力多场耦合下高放废物处置库核心处置单元间距设计研究[J]. 岩土力学, 2025, 46(8): 2626-2638.
[3] 梁庆国, 李景, 张崇辉, 刘彤彤, 孙志涛, . 基底均匀膨胀作用下黄土−泥岩复合地层隧道衬砌力学响应研究[J]. 岩土力学, 2025, 46(6): 1811-1824.
[4] 朱先祥, 张琦, 马俊鹏, 王永军, 孟凡贞, . 浆−水置换效应下含水砂层渗透注浆扩散机制[J]. 岩土力学, 2025, 46(6): 1957-1966.
[5] 杨明云, 陈川, 赖莹, 陈云敏. 串联锚在黏土中的三向受荷承载力分析[J]. 岩土力学, 2025, 46(2): 582-590.
[6] 张凌博, 孙宜松, 程星磊, 郭群录, 赵川, 刘京红. 基于损伤能量耗散的三维土体切削破坏面表征方法研究[J]. 岩土力学, 2025, 46(11): 3626-3636.
[7] 张昕晔, 刘志伟, 翁效林, 李铉聪, 赵建崇, 刘小光. 上砂下黏复合地层隧道开挖面稳定性及破坏模式研究[J]. 岩土力学, 2025, 46(11): 3637-3648.
[8] 吴迪, 陈嵘, 孔纲强, 牛庚, 缪玉松, 王振兴. 冷-热循环温度下桥梁能量排桩热-力响应特性现场试验与数值模拟[J]. 岩土力学, 2025, 46(11): 3649-3660.
[9] 许国庆, 黄高翔, 王协康, 罗登泽, 李洪涛, 姚强, . 新型弧形聚能爆破作用下的岩石破裂演化机制研究[J]. 岩土力学, 2025, 46(10): 3267-3279.
[10] 王帅, 王豫徽, 王玲, 李佳祺, 赵梓皓, 庞凯旋, . 基于晶体模型的岩石孔隙结构与矿物组成对裂纹扩展影响机制研究[J]. 岩土力学, 2025, 46(10): 3289-3301.
[11] 杨立. 平板载荷试验数值分析及承载力判定标准研究[J]. 岩土力学, 2024, 45(S1): 723-730.
[12] 薛秀丽, 谢伟睿, 廖欢, 曾超峰, 陈宏波, 徐长节, 韩磊, . 邻近深埋地铁车站水−土阻隔效应及其对基坑抽水致沉的影响[J]. 岩土力学, 2024, 45(9): 2786-2796.
[13] 吕茂淋, 朱珍德, 周露明, 葛鑫梁, . 基于相场法的预制双裂隙岩体水力压裂扩展数值模拟研究[J]. 岩土力学, 2024, 45(6): 1850-1862.
[14] 马登辉, 韩迅, 蔡正银, 关云飞, . 静压桩的桩侧土压力分布规律数值分析[J]. 岩土力学, 2024, 45(6): 1863-1872.
[15] 陈磊, 张强, 贾朝军, 雷明锋, 黄娟, 胡晶, . 强降雨对库岸堆积体边坡稳定性影响的离心模型试验和数值模拟研究[J]. 岩土力学, 2024, 45(5): 1423-1434.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!