一种电动辅助微生物修复过程的数值模拟方法技术

本发明专利技术公开了一种电动辅助微生物修复过程的数值模拟方法，属于数值模拟方法体系，首先通过初始化研究区设置，确定研究区的概念模型、化学物质设置以及初始条件和边界条件。随后，基于网格的自适应剖分，设置模拟节点数量。基于迁移

【技术实现步骤摘要】
一种电动辅助微生物修复过程的数值模拟方法


[0001]本专利技术涉及一种数值模拟方法体系，具体涉及一种可用于精细刻画电动辅助微生物修复过程的数值模拟方法。

技术介绍

[0002]近年来由于经济社会的发展，部分地块受到了一定程度的土壤和地下水污染，而低渗透介质的存在给水
‑
土修复相关技术的应用带来了巨大的挑战。电动辅助微生物修复为这一领域的修复提供了一个新的思路，而数值模拟修复过程可以预测、评估物质迁移和修复效率，为修复过程提供指导。
[0003]目前已有的地下水数值模拟方法大多限于水力梯度控制流或一维条件，且模拟反应较为简单，对电动辅助微生物修复涉及的多场耦合过程，微生物生长反应过程刻画不全面。进行电动辅助微生物修复过程数值模拟需要进行多场
‑
多反应耦合，其技术门槛较高，且应用时操作复杂，建立通用的电动辅助微生物修复数值模拟方法体系可实现数值模拟效率和精度的提升。COMSOL和PHREEQC作为多场耦合和反应领域通用的模拟软件可作为方法构建基础，进一步拓展并应用于电动辅助微生物修复数值模拟过程。

技术实现思路

[0004]专利技术针对现有地下水数值模拟程序不能较好的模拟电动辅助微生物修复过程，以及模拟过程中出现难以精准评估微生物生长过程等的不足，构建基于COMSOL和PHREEQC的数值模拟方法体系，利用COMSOL来计算物质迁移的对流、弥散、电迁移项等；利用PHREEQC来处理水相物质分配和微生物生长及污染物质消耗的多元Monod反应，考虑微生物修复过程中涉及的微生物生长、降解、介质阻塞等过程，并用微生物生长阻塞模型，GBM模型来评估微生物生长对渗透系数的影响。最终拓展建立的EK
‑
BIO耦合模型，实现对电动辅助微生物修复的精确刻画。
[0005]本专利技术是这样实现的：
[0006]一种电动辅助微生物修复过程的数值模拟方法，其特征在于，所述的方法步骤如下：
[0007]步骤一、对COMSOL和PHREEQC软件内的设置进行初始化，确定研究区的概念模型、涉及研究区结构定义、化学物质设置，初始和边界条件设置，并基于网格的自适应剖分，设置模拟节点数量；
[0008]步骤二、利用运移
‑
反应算子分离法，分别计算迁移与反应过程；其中，迁移项由COMSOL进行多场耦合环境刻画，具体包括水力梯度下的对流，浓度梯度下的弥散，电势梯度下的电迁移、电渗流及可能发生的电极反应；反应项则由PHREEQC进行水相物质分配、以及对微生物修复多元抑制Monod方程的刻画；
[0009]步骤三、在中定义微生物生长对渗透系数影响评估的模型，即微生物
生长阻塞模型(GBM)的参数设置；本专利技术提出的用于刻画微生物生长对渗透系数影响的微生物生长阻塞模型(GBM)，此模型建立了微生物生长同渗透系数之间的变化关系，使得电动辅助微生物修复模型更为精确。
[0010]步骤四、迁移与反应模拟耦合过程是逐步长进行的：本步长内，COMSOL计算得到的物质空间分布结果被导入至PHREEQC内计算各类反应；之后PHREEQC中计算得到的新浓度数据被视为初始条件带入下一个步长进行运移模拟；由MATLAB内编译的子程序控制模拟步长及模拟时长；基于自编译得到耦合计算模块，组合形成CPM
‑
BIO数值模拟流程，实现精细模拟电动辅助微生物修复过程；利用自编译算子分离的逐步长耦合过程。
[0011]步骤五、在完成指定数目步长的数值模拟后，将所有结果导入至COMSOL或在外部程序中进行后处理。
[0012]进一步，步骤二中基于能斯特
‑
普兰克
‑
泊松的多场环境刻画：
[0013]宏观尺度下，使用能斯特
‑
普朗克方程来描述某种物质i移动行为，其通量表达式主要由弥散项活度项和电迁移项组成：
[0014][0015]其中是水相物质扩散系数，c
i
是水相物质浓度，γ
i
是活度系数，z
i
离子价态，F是法拉第常数，R理想气体常数，T是温度，Φ是电势；
[0016]在地下水含水介质中，由于固相介质存在，使得物质i在水相里的扩散受到限制；因此，式1中的水相扩散系数被替换为D
i
＝D
iaq
τ，其中τ为迂曲率；同样，介质孔隙度n的大小对物质的扩散也又有类似影响；除此之外，总的通量考虑对流的影响
[0017][0018]依据驱动力的不同，对流项，q
*
主要分为两类：一、外加水力梯度条件下，会产生压力水流，即传统意义上的地下水流；二、外加电场的条件下，会在介质体间产生电渗流；因此，q
*
表示为q
*
＝q+u
e
，其中q是水力梯度下形成的通量，宏观条件下依据达西定律表示：
[0019][0020]其中k
h
为渗透率，ρ为流体密度，g为重力加速度，μ为动力粘滞系数，
▽
h为水头差；
[0021]根据霍姆赫兹
‑
斯莫鲁霍夫斯基理论，电渗流表示为u
e
＝
‑
k
e
▽
Φ，其中k
e
为电渗系数，对应于达西定律中的渗透系数；在介质体内部溶液pH接近中性时，该值为10
‑8‑
10
‑9m2/(V
·
s)；
[0022]由此，带入质量守恒方程，物质i的迁移本构方程主要由下式表示：
[0023][0024]其中，R
i
为源汇项；
[0025]外加电场作用下或者带电离子本身在外力作用下发生定向移动，会影响周围的电势分布，泊松方程被用来定量这一过程：
[0026][0027]其中ρ
e
为电荷密度，表示为ρ
e
＝F∑z
i
c
i
，ε0为真空的介电常数，ε
r
为相对介电常数；
[0028]联合式(2)、(4)、(5)即为能斯特
‑
普兰克
‑
泊松方程，该组方程即描述物质在介质中运移的本构方程。
[0029]进一步，步骤二中是基于多元Monod的微生物修复过程刻画：
[0030]在对微生物降解速率的评估中，使用基于多元抑制的Monod方程来对其中的微生物生长和有机物降解进行刻画；Monod动力学方程可以刻画物质在链式降解过程中发生的竞争性抑制现象：
[0031][0032]其中，C
EA,i
与C
ED
分别为第i种电子受体EA和电子供体ED的浓度；t为时间(s)；最大表观降解速率(1/s)；K
s,ED
和K
s,EA,i
分别为ED和EA的半饱和系数；K
inh,j
为除了本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种电动辅助微生物修复过程的数值模拟方法，其特征在于，所述的方法步骤如下：步骤一、对COMSOL和PHREEQC软件内的设置进行初始化，确定研究区的概念模型、涉及研究区结构定义、化学物质设置，初始和边界条件设置，并基于网格的自适应剖分，设置模拟节点数量；步骤二、利用运移
‑
反应算子分离法，分别计算迁移与反应过程；其中，迁移项由COMSOL进行多场耦合环境刻画，具体包括水力梯度下的对流，浓度梯度下的弥散，电势梯度下的电迁移、电渗流及可能发生的电极反应；反应项则由PHREEQC进行水相物质分配、以及对微生物修复多元抑制Monod方程的刻画；步骤三、在中定义微生物生长对渗透系数影响评估的模型，即微生物生长阻塞模型(GBM)的参数设置；步骤四、迁移与反应模拟耦合过程是逐步长进行的：本步长内，COMSOL计算得到的物质空间分布结果被导入至PHREEQC内计算各类反应；之后PHREEQC中计算得到的新浓度数据被视为初始条件带入下一个步长进行运移模拟；由MATLAB内编译的子程序控制模拟步长及模拟时长；基于自编译得到耦合计算模块，组合形成CPM
‑
BIO数值模拟流程，实现精细模拟电动辅助微生物修复过程；步骤五、在完成指定数目步长的数值模拟后，将所有结果导入至COMSOL或在外部程序中进行后处理。2.根据权利要求1所述的一种电动辅助微生物修复过程的数值模拟方法，其特征在于，步骤二中基于能斯特
‑
普兰克
‑
泊松的多场环境刻画：宏观尺度下，使用能斯特
‑
普朗克方程来描述某种物质i移动行为，其通量表达式主要由弥散项活度项和电迁移项组成：其中是水相物质扩散系数，c
i
是水相物质浓度，γ
i
是活度系数，z
i
离子价态，F是法拉第常数，R理想气体常数，T是温度，Φ是电势；在地下水含水介质中，由于固相介质存在，使得物质i在水相里的扩散受到限制；因此，式1中的水相扩散系数被替换为D
i
＝D
iaq
τ，其中τ为迂曲率；同样，介质孔隙度n的大小对物质的扩散也又有类似影响；除此之外，总的通量考虑对流的影响考虑对流的影响依据驱动力的不同，对流项，q
*
主要分为两类：一、外加水力梯度条件下，会产生压力水流，即传统意义上的地下水流；二、外加电场的条件下，会在介质体间产生电渗流；因此，q
*
表示为q
*
＝q+u
e
，其中q是水力梯度下形成的通量，宏观条件下依据达西定律表示：其中k
h
为渗透率，ρ为流体密度，g为重力加速度，μ为动力粘滞系数，为水头差；
根据霍姆赫兹
‑
斯莫鲁霍夫斯基理论，电渗流表示为其中k
e
为电渗系数，对应于达西定律中的渗透系数；在介质体内部溶液pH接近中性时，该值为10
‑8‑
10
‑9m2/(V
·
s)；由此，带入质量守恒方程，物质i的迁移本构方程主要由下式表示：其中，R
i
为源汇项；外加电场作用下或者带电离子本身在外力作用下发生定向移动，会影响周围的电势分布，泊松方程被用来定量这一过程：其中ρ
e
为电荷密度，表示为ρ
e
＝F∑z
i
c
i
，ε0为真空的介电常数，ε
r
为相对介电常数；联合式(2)、(4)、(5)即为能斯特
‑
普兰克
‑
泊松方程，该组方程即描述物质在介质中运移的本构方程。3.根据权利要求1所述的一种电动辅助微生物修复过程的数值模拟方法，其特征在于，步骤二中是基于多元Monod方程的微生物修复过程刻画；在对微生物降解速率的评估中，使用基于多元抑制的Monod方程来对其中的微生物生长和有机物降解进行刻画；Monod动力学方程可以刻画物质降解过程中发生的竞争性抑制现象：其中，C
EA,i
与C
ED
分别为第i种电子受体EA和电子供体ED的浓度；t为时间(s)；κ
max,EA,Xk
最大表观降解速率(1/s)；K
s,ED
和K
s,EA,i
分别为ED和EA的半饱和系数；K
inh,j
为除了第i种EA之外的其他j种EA对目标EA的抑制系数；X
k
为微生物数量；对降解过程引起的微生物数量随时间的变化，则由下式动力学反应来刻画：其中，Y为微生物生长系数，b
k
为微生物消亡系数。4.根据权利要求1所述的一种电动辅助微生物修复过程的数值模拟方法，其特征在于，步骤三中微生物生长阻塞模型GBM中：实际修复过程中，须评估生物生长过程对水动力参数的影响，GBM模型在进行迁移
‑
反应模拟时考虑微生物存在的两种相态，即吸附在固体表面上和游离于液体中，假定吸附相微生物是造成介质阻塞的原因；两种相态的微生物进一步考虑活性和死亡两种形态，能进行微生物降解的活性生物及不能进行微生物降解的死亡生物在两相上均有分配；模型的相关设定如下：(1)设定存在最大吸附生物量(X
surface
_
max
)和最大渗透系数降低倍率(P
max

【专利技术属性】
技术研发人员：叶淑君，杨平，吴吉春，吴剑锋，
申请(专利权)人：南京大学，
类型：发明
国别省市：