1. 引言
煤矸石作为工业固体废渣之一,将其应用于新建或改扩建高速公路路堤层位,不仅具有良好的经济价值和社会效益,而且可以实现废旧材料的循环利用 [1],最终达到保护环境和节约土地的目的。但在煤矸石填筑路基工程中,煤矸石中存在的化学成分中包含重金属及有害物质 [2],导致路基运维期的土壤、地表水以及地下水污染。降雨可使煤矸石中的污染物淋溶并渗入到地下,造成周围的地下水污染。研究煤矸石中污染物释放传输的动力学行为对污染物的减排监控具有重要的科学意义。煤矸石释放污染组分向环境释放传输是一个复杂物理和化学变化过程 [3]。
在煤矸石与水反应的初始阶段,矿物的分解速率要快一些,反应液中有害离子的浓度会不断增加。随着反应的进行,水中溶解氧的消耗,矿物的分解速率将逐渐减慢,溶液中污染组分的浓度将趋于平衡。污染组分进入溶液中主要的过程有:反应物通过孔隙向煤矸石中扩散;在扩散途中,溶解氧与相遇到的矿物在其表面发生氧化还原反应导致矿物分解;溶解下来的污染物向溶液中扩散;在离开固相之前,其中的一部分污染物会被骨架基质吸附;在溶液中有一些污染物可能因为酸度的改变形成沉淀或被胶体和残余骨架中的有机物所吸附。其溶解反应过程不仅受到煤矸石的物理化学性质所制约,同时也与煤矸石的岩相组成、水分含量、煤矸石的比表面积、孔隙率以及堆体内部温度等有关 [4]。因此,要弄清溶解作用、沉淀作用、离子交换和吸附作用对污染物浓度分布贡献程度,必须建立煤矸石污染组分的水化反应动力学模型,定量描述其反应的程度,同时从细观角度分析温度水流耦合作用煤矸石内污染物演化的动力学特征。
煤矸石内污染物对于环境条件的变化是敏感的,特别是温度、水分、pH值以及煤矸石在表生环境条件下矿物发生氧化反应、吸附作用、沉淀作用、离子交换作用等一系列反应控制着污染物释放传输过程。开展煤矸石中污染物质释放传输规律的研究不仅能为煤矸石污染控制及资源化利用奠定理论基础,而且对预防煤矸石中污染物对生态环境的影响以及环境危害性评价也具有理论指导意义和现实意义 [5]。很多学者对煤矸石淋溶污染物传输动力学进行了研究,本文从不同学者的不同角度对近年来的研究进展和成果进行综述。
2. 淋溶污染物迁移的耦合动力学模型
通常情况下,煤矸石可被看成非饱和的土壤。土壤中水的流动形式属于多孔介质中流体运动。雨水通过煤矸石后会含有一定浓度的污染物,这算是淋溶的结果。煤矸石路基淋溶所产生的污染物在土壤–水系统中存在扩散、吸附、分配以及生物降解等过程。当这些过程不能完全降解煤矸石淋溶所产生的污染物时,污染物就会随着雨水的下渗到土壤和地下水中,从而会给道路沿线的环境造成污染。根据污染物随水在非饱和土壤中的运动扩散规律,建立煤矸石中污染物迁移转化的耦合动力学模型 [6]。
2.1. 土壤水分运移数学模型的建立
淋溶煤矸石路堤所产生的污染物在土壤–水系统中进行的扩散是三维流体运动 [7],假设土壤为各向同性、均质的多孔介质,根据淋溶煤矸石所产生的污染物随水分迁移的规律,建立了煤矸石路堤下土壤水分运移的一维控制数学模型为:
(1)
式中:,但在实际计算中常取(L为某一时刻或距离);z为污染物垂向运移距离,取向下为正,cm;t为时间,min;T为某一时刻,min;θ为土壤含水量,g/cm3;θm为初始含水量,g/cm3;θb为上边界含水量,g/cm3;K(θ)、D(θ)分别为非饱和土壤导水率和扩散率;i(t)为降雨强度,mm/min。
取;。其中θa为土壤饱和含水量,g/cm3;D0、n、a、b为经验参数;s为土壤水吸力,kPa。由于比水容量,进而根据可得:
(2)
2.2. 煤矸石中污染物迁移数学模型的建立
假设稳定水流条件下,同时考虑深度变化下的一价降解和随深度变化的线性平衡吸附时,一维的煤矸石中污染物迁移的数学模型如下:
(3)
式中:为土壤密度,g/cm3;Dz为水动力纵向弥散系数,cm;c为土壤中污染物的浓度,mg/L;S为污染物吸附浓度,mg/L;、分别为水在z和y方向的速度,cm/s;和分别为随深度变化的降解系数和吸附系数,d−1;和分别为降解系数和吸附系数的常数,代表和的最大值,d−1;为经验参数;ch为双曲余弦函数;y为污染物水平方向运移距离,取向右为正,m。
第一类初始条件:
(4)
(5)
第二类初始条件:
(6)
(7)
式中:为土壤中污染物的初始质量浓度,mg/L;D为水动力弥散系数,cm;c为土壤中污染物的浓度,mg/L。第一类初始条件考虑污染物浓度的z向迁移,第二类初始条件考虑的是水的动力扩散和弥散。
2.3. 数学模型的求解
土壤水分运移模型和污染物迁移模型是耦合存在的,需要联立求解。先求解土壤水分运移模型的解定流速场,然后将其带入污染物的迁移模型中求解,即可得到土壤中污染物的浓度分布。土壤水分运移模型和煤矸石中污染物的迁移模型都是偏微分方程,一般条件下无法得到解析解,只有数值解。
利用建立的模型,根据经验选取相应的计算参数,可以模拟计算大气降雨淋溶煤矸石产生的污染物在土壤中的扩散迁移分布情况 [8]。所建立的模型可以为评价煤矸石作为路基填料的使用而对道路沿线带来的环境污染问题提供理论依据。
3. 通过耗氧量耦合的动力学模型
在煤矸石反应过程中,控制着化学反应发生的黄铁矿在氧化反应中反应生成Fe3+,同时还有H+,其量取决于地下水的化学成分,同时也取决于耗氧量,此类动力学模型 [9] [10] [11] 就是通过耗氧量将有关式子耦合。其反应方程为:
(8)
有关的式子为:
(9)
式中,为含水率;t为时间;为负压水头,为空间坐标;为非饱和水压传导系数,其计算公式为
(10)
式中,为饱和水压传导系数;为有效饱和度,其计算公式为
(11)
其中为容器中多孔介质的饱和含水率;为容器中的残留含水率;含水率和其它参数之间的关系为
(12)
式中,及均为经验系数。
在传输过程中,忽略水蒸汽的弥散作用,在平衡吸附的条件下,二维溶质传输通过对流弥散方程表示,即
(13)
式中,C为微量元素在水中的浓度;为水动力学弥散系数;为达西流速;R为源汇项。
(14)
其中为弥散度;为多孔介质中有效分子弥散系数。饱和、非饱和流动模型方程和微量元素传输方程之间的关系通过孔隙水流速耦合,其控制方程为
(15)
式中,v为孔隙水流速。
氧在气相中的弥散传输一维基本控制方程为
(16)
式中,为含气量;为氧在气相中的弥散系数;Q为氧的消耗量。
通过以上模型可以预测微量元素在地下环境中的分布及扩散趋势,再结合研究区域的水文地质资料可以有效地提出污染防治的控制措施。通过编制的计算程序可以定量分析微量元素在不同质地土壤和地下水中分布的规律,从而来优化污染控制的技术方案和工艺参数,这可以为污染控制提供科学的技术支持;同时,可以结合各种实验参数,模拟不同污染强度条件下土壤–水环境系统中微量元素的动态变化规律,预测污染的范围与程度,根据预测结果,提出防治煤矸石淋溶液污染的具体措施与治理方案。
4. 水分运移和重金属运移模型耦合
煤矸石淋溶的重金属在包气带土壤中的迁移、转化主要是由土壤水运移及重金属与土壤间的各种物理、化学吸附引起的。采用数学模型来描述重金属污染物对地下水潜在性污染过程,包括两类控制方程:一类包气带土壤水分运移;另一类是重金属污染物运移模型。为此,建立如下耦合数学模型 [12] [13]。
4.1. 土壤水分迁移的数学模型
对于土壤中水分运动可用如下数学模型
(17)
其中为比水容量、分别为横向和纵向水力传导系数;h为压力水头。
4.2. 重金属污染物运移的数学模型
对于重金属污染物运移采用如下数学模型:
(18)
其中θ为土壤水分含量;C为污染物在溶解相中的质量浓度;S为污染物在吸附相中的质量浓度;土壤体积密度;γ为最大沉淀量:和分别为横向和纵向水动力弥散系数;q为达西流速;为沉淀速率;为解沉淀速率;为初始溶解相污染物质量浓度。
水分迁移的达西流速为
(19)
4.3. 耦合数学模型的耦合求解
由土壤水分运移控制方程和污染物运移控制方程构成了重金属在地下环境系统中运移的数值模型。
对土壤水分运移的控制方程,令
(20)
则水分控制方程的差分格式为
(21)
令(22)
(23)
(24)
(25)
将上式联立可得,
(26)
对应矩阵表达式为:
(27)
其中[T]为系数矩阵;[F]为常数项列阵;{h}为待求的水头列阵。
对于溶质运移的控制方程对应的矩阵表达时为
(28)
(29)
其中、、、分别为弥散系数张量、孔隙土壤水流速系数矩阵;、分别为污染物结点浓度列矩。
通过本文水分运移和重金属迁移转化耦合数值模型的研究,可以再现煤矸石淋溶重金属在地下环境系统浓度分布以及演化与发展,并可对不同质地土壤进行定量化分析,其研究结果不仅可为煤矸石所引起的地下水中污染物的运移及其对环境影响的预测、预报提供科学手段。
5. 其他动力学模型
5.1. 淋滤液在地下水中迁移的动力学数学模型
淋滤液在地下水中迁移的动力学行为是由地下水流和污染物运移组成。为了进行地下水污染预测,须同时求解地下水流控制方程和污染物在地下水中迁移的控制方程 [14]。
第四纪含水层中z方向受到底板隔水层的限制,所以建立了地下水在含水层中的二维运移方程,
(30)
式中,Ss——贮水系数(L−1);h——压力水头(L);Kxx,Kyy——分别为横向和纵向渗透系数(LT−1);t——时间(T);qs——单位体积的源的体积流率(T−1)。
在综合考虑淋滤液在地下水环境体系中发生对流、弥散、吸附解吸以及存在源汇项情况下,建立淋滤液迁移转化的控制方程:
(31)
式中,——有效孔隙度;C——溶质在地下水中的浓度(ML−3);DL,DT ——纵向、横向弥散系数(L2T−1);vx,vy——x,y方向平均孔隙流速(LT−1);Cs——源的污染物质浓度(ML−3);Rd——迟滞因子。
地下水运移的定解初始条件:
,(32)
边界条件:
,(33)
,,(34)
淋滤液在地下水中迁移的定解初始条件:
,(35)
边界条件:
,(36)
,(37)
式中,——初始压力水头值(L);——边界压力水头值(L);n——边界Γ2的外法线方向;q——单位面积上流出的流量(LT−1);——溶质初始浓度;C1——污染物在边界处的浓度。
煤矸石淋滤液在地下水中迁移的动力学数学模型,可以采用有限差分的方法对模型进行了数值求解,为定量化预测、预报污染物浓度变化趋势及采取有效及时的污染防治措施提供了理论依据。
5.2. 孔隙–裂隙双重介质模型
煤矸石的淋溶液淋溶后先后经过土壤的表土层、下包气带而进人含水层,对地下水造成一定程度的污染。在含有孔隙裂隙的地下含水层中,裂隙发育的孔隙裂隙岩层中同时存在着两种不同孔隙率的渗流系统:孔隙总体积较大但渗透性相对弱的多孔岩块和总体积较小但渗透性却相对较强的分割多孔块的裂隙.岩层中的水流大部分是由体积相对较小的那部分裂隙传输的,地下水主要贮存在孔隙中,水的运动主要在裂隙中进行,为此,建立如下孔隙–裂隙双重介质模型 [15]。
1) 地下水水流模型
裂隙域:
(I) (38)
岩块域:
(II) (39)
2) 淋溶液污染物迁移模型
裂隙域:
(III) (40)
岩块域:
(IV) (41)
(42)
(43)
式中G、G'分别表示裂隙域和孔隙岩块域;表示裂隙域的边界;f1和f2分别表示垂向补给裂隙和岩块域的补给强度;h和h'分别表示裂隙和岩块域的水头分布;B表示孔隙岩块G'边界;表示单位时间单位面积裂隙一岩块界面上的变换量;M为含水层厚度,n、b分别为岩块厚度和裂隙宽度之半;C和C分别为裂隙域和孔隙岩块域的浓度分布;n为孔隙度;R和R'分别为裂隙域和孔隙岩块域的延迟系数;dc/dt价表示化学动力学溶解速度;fc为含水层补给量或排泄量中的溶质浓度;D和D'分别为裂隙域的弥散系数和岩块域分子扩散系数;为边界B的单位外法向量;π为特定参数。
淋溶液污染物在含水层中迁移的双重介质模型,可为淋溶液污染机理的分析提供理论依据。
5. 结论
在释放传输特性定量化研究方面,现在文献主要从矿物学、地球化学及模拟试验方面对煤矸石研究颇多,使得煤矸石矿环境污染综合治理上取得了较大进展,但是应考虑多场耦合条件下煤矸石内污染物释放特性,尤其需要从定量角度对煤矸石衍生环境效应的过程和机制进行研究,这是研究的一个发展方向。在溶解反应动力学研究方面,对水–岩相互作用的动力学的理论研究和数值模拟已比较透彻,但对于水–岩化学作用的非平衡态非线性动力学理论的研究还处于初始阶段,这是另一个研究的发展方向。
基金项目
山东省交通运输厅科技计划(2019B24)。