基于图像法的多孔介质双分子反应溶质运移模拟

doi:10.13745/j.esf.sf.2022.1.28

地学前缘 ›› 2022, Vol. 29 ›› Issue (3): 248-255.DOI: 10.13745/j.esf.sf.2022.1.28

基于图像法的多孔介质双分子反应溶质运移模拟

刘咏¹(), 张琪², 钱家忠²^,^*(), 吴盾³, 张文永³

1.合肥工业大学食品与生物工程学院, 安徽合肥 230601
2.合肥工业大学资源与环境工程学院, 安徽合肥 230009
3.安徽省煤田地质局勘探公司, 安徽合肥 230088

收稿日期:2021-07-01 修回日期:2022-02-18 出版日期:2022-05-25 发布日期:2022-04-28
通讯作者: 钱家忠
作者简介:刘咏(1968—),女,教授,硕士生导师,主要从事地下水环境污染防治与修复研究工作。E-mail: liuyong@hfut.edu.cn
基金资助:
国家自然科学基金项目(41877191);国家自然科学基金项目(41831289);安徽省重点研究和开发计划项目(201904a07020071)

Simulation of bimolecular reactive solute transport in porous media via image analysis

LIU Yong¹(), ZHANG Qi², QIAN Jiazhong²^,^*(), WU Dun³, ZHANG Wenyong³

1. School of Food and Biological Engineering, Hefei University of Technology, Hefei 230601, China
2. School of Resources and Environmental Engineering, Hefei University of Technology, Hefei 230009, China
3. Exploration Research Institute, Anhui Provincial Bureau of Coal Geology, Hefei 230088, China

Received:2021-07-01 Revised:2022-02-18 Online:2022-05-25 Published:2022-04-28
Contact: QIAN Jiazhong

摘要/Abstract

摘要：

采用图像分析法确定溶质浓度,以硫酸铜和EDTA二钠作为双分子反应物,在多孔介质模型中开展了不同粒径(1.52.0,2.53.0,3.54.0 mm)和流量(1.0,1.5,2.0 mL/s)下反应性溶质运移实验,探讨了应用不完全混合的对流弥散模型(IM-ADRE)对双分子反应溶质运移的模拟和预测,并进行了参数敏感性分析。结果表明:图像分析法可准确获取多孔介质中显色反应性溶质的浓度,灰度值与浓度的决定系数R²大于0.96;用IM-ADRE模型能够准确预测双分子反应性溶质硫酸铜和EDTA二钠在3种不同多孔介质中的运移过程,误差低于3.71%;实验条件的改变对IM-ADRE模型参数D、m和β₀的影响显著,说明模型参数依赖于环境条件,其变化规律需要根据实际环境条件进一步率定,便于IM-ADRE模型的进一步推广应用。

关键词: 双分子反应性溶质运移, 多孔介质, 模拟, IM-ADRE模型, 图像法

Abstract:

In this study, reactive solute transport experiments in porous media were carried out under different particle sizes (ϕ1.5-2.0, 2.5-3.0, and 3.5-4.0 mm) and flow rates (1.0, 1.5, and 2.0 mL/s). We used copper sulfate and disodium EDTA as bimolecular reactors and simulated their transport via image analysis. The incomplete mixing advection-dispersion reaction equation (IM-ADRE) based simulation and prediction of bimolecular reactive solute transport was investigated, and parameter sensitivity analysis was performed. The results showed that the concentration of color reactive solutes in porous media can be accurately obtained via image analysis. The determination coefficient between the gray value and concentration was greater than 0.96. The IM-ADRE model gave high accuracy predictions for bimolecular reactive solute transport in porous media under different conditions, with less than 3.71% relative errors. The change of experimental conditions had significant effects on parameters D, M and β₀ of IM-ADRE model, indicating that the model parameters are dependent on the environmental conditions, and the dependency rules need to be further calibrated according to actual environmental conditions, which can help promoting IM-ADRE model applications.

Key words: bimolecular reactive solute transport, porous media, simulating, IM-ADRE model, image method

中图分类号:

P641.2
X143

刘咏, 张琪, 钱家忠, 吴盾, 张文永. 基于图像法的多孔介质双分子反应溶质运移模拟[J]. 地学前缘, 2022, 29(3): 248-255.

LIU Yong, ZHANG Qi, QIAN Jiazhong, WU Dun, ZHANG Wenyong. Simulation of bimolecular reactive solute transport in porous media via image analysis[J]. Earth Science Frontiers, 2022, 29(3): 248-255.

图/表 6

图1 实验装置

Fig.1 Experimental setup

图2 50 cm和100 cm处溶质的灰度值-浓度标准曲线

Fig.2 The standard curve of gray value-concentration for each solute at 50 cm and 100 cm

图3 不同粒径和流量条件下双分子反应溶质运移穿透曲线(100 cm处) 图中A-S,B-S,AB-S与A-E,B-E,AB-E分别为Na2EDTA2-,CuSO4,CuEDTA2-的实验值与拟合值;(a),(b),(c)分别代表流量为1.0 mL/s的3组粒径1.5~2.0,2.5~3.0,3.5~4.0 mm; (d),(e),(f)分别代表流量为1.5 mL/s的3组粒径1.5~2.0,2.5~3.0,3.5~4.0 mm; (g),(h),(i)分别代表流量为2.0 mL/s的3组粒径1.5~2.0,2.5~3.0,3.5~4.0 mm。

Fig.3 BTCs of bimolecular reactive solute transport under different flow rates and grain sizes (at 100 cm)

表1 流量与粒径对IM-ADRE模型拟合参数的影响(100 cm 处)(据文献[10])

Table 1 Effect of flow rate and grain size on the fitting parameters of IM-ADRE model (at 100 cm). Adapted from [10].

流量/(mL·s^-1)	粒径ϕ /mm	弥散系数D/(10^-5 m²·s^-1)	β₀/(m³·mol^-1·s^-1)	m	峰值拟合误差/%
1.0	1.52.0	1.27	698	1.26	1.30
	2.53.0	1.43	652	0.90	3.31
	3.54.0	2.30	535	0.41	3.71
1.5	1.52.0	1.44	699	1.19	1.41
	2.53.0	1.86	665	0.82	0.61
	3.54.0	2.79	598	0.36	1.39
2.0	1.52.0	1.89	712	1.18	2.52
	2.53.0	2.02	683	0.75	2.42
	3.54.0	2.79	604	0.28	2.09

图4 不同粒径和流量条件下反应性溶质运移穿透曲线(50 cm处) 图中A-S,B-S,AB-S与A-E,B-E,AB-E分别为Na2EDTA2-,CuSO4,CuEDTA2-的实验值与拟合值;(a),(b),(c)分别代表粒径为1.5~2.0 mm的3组流量1.0,1.5,2.0 mL/s;(d),(e),(f)分别代表粒径为2.5~3.0 mm的3组流量1.0,1.5,2.0 mL/s;(g),(h),(i)分别代表粒径为3.5~4.0 mm的3组流量1.0,1.5,2.0 mL/s。

Fig.4 BTCs of reactive solutes migration under different flow rates and grain sizes (at 50 cm)

表2 粒径和流量对IM-ADRE模型拟合参数的影响(50 cm处)(据文献[10])

Table 2 Effects of grain size and flow rate on the fitting parameters of IM-ADRE model (at 50 cm). Adapted from [10].

粒径ϕ /mm	流量/(mL·s^-1)	弥散系数D/(10^-5 m²·s^-1)	β₀/(m³·mol^-1·s^-1)	m	峰值拟合误差/%
1.52.0	1.0	0.97	672	1.23	3.33
	1.5	1.16	684	1.16	0.22
	2	1.49	692	1.14	0.89
2.53.0	1.0	1.17	628	0.86	0.96
	1.5	1.53	645	0.78	0.18
	2	1.69	657	0.75	1.04
3.54.0	1.0	1.86	535	0.41	3.71
	1.5	2.24	579	0.30	0.26
	2	2.79	604	0.28	0.22

参考文献 14

[1]	LIU Y J, LIU Y, LUO Q K, et al. Experimental and numerical study of bimolecular reactive transport in a single rough-wall fracture[J]. Journal of Hydrology, 2021, 594: 125944.
[2]	GRAMLING C M, HARVEY C F, MEIGS L C. Reactive transport in porous media:a comparison of model prediction with laboratory visualization[J]. Environmental Science and Technology, 2002, 36(11): 2508-2514. DOI URL
[3]	HAGGERTY R, HARVEY C F, VON SCHWERIN C F, et al. What controls the apparent timescale of solute mass transfer in aquifers and soils? A comparison of experimental results[J]. Water Resources Research, 2004, 40(1): W01510.
[4]	SANCHEZ-VILA X, FERNÀNDEZ-GARCIA D, GUADAGNINI A. Interpretation of column experiments of transport of solutes undergoing an irreversible bimolecular reaction using a continuum approximation[J]. Water Resources Research, 2010, 46(12): W12510.
[5]	QIAN J Z, ZHAN H B, ZHANG Y, et al. Numerical simulation and experimental study of bimolecular reactive transport in porous media[J]. Transport in Porous Media, 2015, 109 (3): 727-746. DOI URL
[6]	LIU Y, QIAN J Z, MA L, et al. Improvement of numerical simulation of bimolecular reactive transport using time-dependent parameters[J]. Journal of Hydrodynamics, 2015, 27(1): 62-67. DOI URL
[7]	RAJE D S, KAPOOR V. Experimental study of bimolecular reaction kinetics in porous media[J]. Environmental Science and Technology, 2000, 34(7): 1234-1239. DOI URL
[8]	QIAN J Z, WANG Z K, GARRARD R M, et al. Non-invasive image processing method to map the spatiotemporal evolution of solute concentration in two-dimensional porous media[J]. Journal of Hydrodynamics, 2018, 30(4): 758-761. DOI URL
[9]	钱家忠, 刘咏, 张勇, 等. 多孔介质双分子反应溶质运移实验与模拟研究[J]. 中国科学技术大学学报, 2013, 43(2): 104-109.
[10]	张琪. 多孔介质中双分子反应物质运移及尺度依赖性研究[D]. 合肥: 合肥工业大学, 2020.
[11]	郑志成. 裂隙中双分子反应性溶质运移及其尺度依赖性研究[D]. 合肥: 合肥工业大学, 2018.
[12]	高晓梅. 多孔介质中反应性溶质运移机理研究: 以CuSO₄和EDTA-Na为例[D]. 合肥: 合肥工业大学, 2015.
[13]	徐玉璐. 多孔介质中污染物运移及弥散系数对流速依赖性实验研究[D]. 合肥: 合肥工业大学, 2017.
[14]	PORTA G M, RIVA M, GUADAGNINI A. Upscaling solute transport in porous media in the presence of an irreversible bimolecular reaction[J]. Advances in Water Resources, 2012, 35: 151-162. DOI URL

[1]	彭光荣, 杜晓东, 姜素华, 李三忠, 关雪婷, 魏董, 姜衍, 陆蕾蕾. 基于古地貌和源-汇特征分析的烃源岩评价:以珠江口盆地阳江凹陷为例[J]. 地学前缘, 2022, 29(5): 188-202.
[2]	孙辉, 刘晓东. 青藏高原隆升气候效应的数值模拟研究进展概述[J]. 地学前缘, 2022, 29(5): 300-309.
[3]	孙炜毅, 刘健, 严蜜, 宁亮. 全新世亚洲季风百年-千年尺度变化的模拟研究进展[J]. 地学前缘, 2022, 29(5): 342-354.
[4]	燕青. 亚洲高山区全新世中期气候及其影响下的冰川模拟[J]. 地学前缘, 2022, 29(5): 372-381.
[5]	何雨江, 刘肖, 邢林啸, 谭现锋, 卜宪标. 河北保定岩溶地热结垢过程模拟及防垢对策[J]. 地学前缘, 2022, 29(4): 430-437.
[6]	李栋, 赵敏, 刘再华, 陈波. 普定岩溶水-碳循环模拟试验场水体双碳同位素特征(δ¹³C-Δ¹⁴C)与碳足迹[J]. 地学前缘, 2022, 29(3): 155-166.
[7]	李东东, 高丹丹, 边绍菊, 李武, 董亚萍. Na⁺,K⁺,Mg²⁺,Ca²⁺//Cl^-,SO₄^2--H₂O六元体系中杂卤石形成条件的再认识[J]. 地学前缘, 2021, 28(6): 46-55.
[8]	沈晓芳, 万玉玉, 王利刚, 苏小四, 董维红. 基于多相流数值模拟的某石油污染场地地下水中VOCs自然衰减过程识别及能力评估[J]. 地学前缘, 2021, 28(5): 90-103.
[9]	安文通, 陈建平, 朱鹏飞. 基于成矿过程数值模拟的隐伏矿双向预测研究[J]. 地学前缘, 2021, 28(3): 97-111.
[10]	肖凡, 王恺其. 德兴斑岩铜矿床断裂与侵入体产状对成矿的控制作用:从力-热-流三场耦合数值模拟结果分析[J]. 地学前缘, 2021, 28(3): 190-207.
[11]	张华添, 李江海, 陶春辉. 西南印度洋中脊斜向扩张分段特征及构造成因探讨[J]. 地学前缘, 2021, 28(2): 271-283.
[12]	罗琪, 黄时卓, 史德锋, 汪锐. 利用能量半衰属性识别涠西南凹陷砂泥薄互层组合中浊积砂体[J]. 地学前缘, 2021, 28(1): 273-281.
[13]	刘成林, 张禹, 杨晟颢, 李宗星, 田继先, 彭博, 马寅生, 杨元元, 孔骅. 海相低勘探程度地区油气资源评价:以柴达木盆地德令哈坳陷石炭系为例[J]. 地学前缘, 2021, 28(1): 295-307.
[14]	王殿举, 李江海, 李一赫. 下地壳流变学性质对南大西洋两岸盆地结构不对称性的控制[J]. 地学前缘, 2020, 27(3): 254-261.
[15]	孔彦龙, 黄永辉, 郑天元, 陆仁超, 潘晟, 邵亥冰, 庞忠和. 地热能可持续开发利用的数值模拟软件OpenGeoSys:原理与应用[J]. 地学前缘, 2020, 27(1): 170-177.

基于图像法的多孔介质双分子反应溶质运移模拟

Simulation of bimolecular reactive solute transport in porous media via image analysis

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

图/表 6

参考文献 14

相关文章 15

编辑推荐

Metrics

本文评价