基于电阻率法的人工湿地堵塞区域探测方法

丁彦礼; 肖敬瑞; 白少元; 熊彬; 孙琳

doi:10.12030/j.cjee.202109076

基于电阻率法的人工湿地堵塞区域探测方法

1.
桂林理工大学地球科学学院，桂林 541004
2.
广西环境污染控制理论与技术重点实验室，桂林 541004
3.
桂林理工大学环境科学与工程学院，桂林 541004

作者简介: 丁彦礼(1979—)，男，硕士，副教授，dyl@glut.edu.cn

通讯作者: 白少元(1981—)，女，博士，教授，baisy@glut.edu.cn;

基金项目:
广西科技重大专项项目(桂科AA20161001)；国家自然科学基金项目(51768013，41404116)
中图分类号: X505

Detection method of constructed wetland clogging based on resistivity method

1.
College of Earth Sciences, Guilin University of Technology, Guilin 541004, China
2.
Guangxi Key Laboratory of Environmental Pollution Control Theory and Technology, Guilin 541004, China
3.
College of Environmental Science and Engineering, Guilin University of Technology, Guilin 541004, China

Corresponding author: BAI Shaoyuan, baisy@glut.edu.cn ;

摘要: 人工湿地污水处理技术已经得到广泛的应用，但目前对人工湿地堵塞问题，无法对堵塞区域进行精确定位。为了解决这个问题，根据潜流人工湿地结构和堵塞区域的特点，建立了长、宽、高分别为146、119和102 cm水槽物理实验模型，利用改进的高密度电阻率测井法开展了探测和定位人工湿地堵塞区域的实验研究，并结合人工湿地堵塞模型的Visual MODFLOW(VMOD)流场模拟特征，探讨了基于电阻率法探测和定位人工湿地堵塞区域的可行性。结果表明，改进的高密度电阻率测井法在电阻率探杆的电极间距不大于5 cm，供电极距D_AB/2不超过12.5 cm的情况下，通过绘制电阻率数据的二维和三维视电阻率等值线图，能够更好地定位高阻体的空间位置和形态。可见，利用改进的高密度电阻率测井方法可以探测和定位人工湿地高阻堵塞区域，实现堵塞区域的精确定位。
- 潜流人工湿地 /
- 堵塞 /
- 改进的高密度电阻率测井法 /
- VMOD /
- 堵塞定位
Abstract: Constructed wetland for sewage treatment has been widely used, but at present, the blocked area in constructed wetland can not be accurately located. To solve the problem, according to the characteristics of subsurface flow constructed wetland structure and blocked area, a physical experimental model of a flume with length, width and height of 146, 119 and 102 cm, respectively, was established. Experimental research on detecting and locating the blocked area of constructed wetland was carried out by using an improved high-density resistivity logging method. Combined with Visual MODFLOW (VMOD) flow field simulation characteristics of the constructed wetland blocking model, the feasibility of detecting and locating the blocked area of the constructed wetland based on the resistivity method was discussed. The experimental results show that the improved high-density resistivity logging method could well locate the spatial position and shape of high resistivity body by drawing two-dimensional and three-dimensional apparent resistivity contours of resistivity data under the condition that the electrode spacing of the resistivity probe rod and the power supply distance of D_AB/2 were less than 5 cm and 12.5 cm, respectively. Therefore, the improved high-density resistivity logging method can be used to detect and locate the high resistance blocked area of constructed wetland, which can realize the accurate location of the blocked area.
- subsurface flow constructed wetland /
- clogging /
- improved high density resistivity logging method /
- VMOD /
- blockage positioning

随着工业化和城镇化的加速推进，对废水的集中处理备受关注^[1]。1932年开始应用的Wuhrmann工艺是最早的脱氮工艺，称之为O/A工艺，遵循硝化、反硝化的流程顺序而设置^[2]。然而，在硝化过程中需要供氧，反硝化过程中需要外加碳源，这造成了能耗和碳源的双重浪费。对此，将生物单元的顺序进行倒置，便产生了A/O工艺，A/O工艺成为最早使用的生物脱氮技术。这是工艺单元不同排列顺序构成组合工艺的开端，后续发展的废水生物处理工艺几乎均为厌氧、缺氧/水解、好氧单元的组装(图1)。典型的工艺有A/A/O和O/A/O，组合工艺中的不同单元反应器排序会影响碳源利用和脱氮效果，因此，需要根据废水组成与处理目标选择合适的工艺技术。

图 1 前置厌氧与前置好氧工艺的演化配置

Figure 1. Evolutionary configuration of pre-anaerobic and pre-aerobic processes.

下载: 全尺寸图片幻灯片

厌氧置前的工艺可以控制碳源转化为小分子有机物或者甲烷，提高废水的可生化性，为后续反硝化反应提供碳源。HAO等^[3]采用A/A/O工艺处理制革废水，考察了沿程溶解性有机物的浓度变化，发现A1的厌氧水解单元能优先去除小分子量的物质和蛋白质，后续的A/O工艺可更彻底地去除残余有机物。O/A/O工艺可在O1单元反应器中好氧降解部分有机物，实现含氮有机物的氨化，有助于硝化反应的实现。李国令等^[4]对比了O/A/O和A/O工艺处理同一城镇污水的结果，在O1单元反应器中降解了大部分有机物，可为O2提供良好的硝化环境，因此，O/A/O脱氮效果优于A/O工艺。A/A/O工艺对高毒性工业废水的处理不具有优势，这是因为A1中的微生物增殖速度慢，难以消除毒性抑制作用。兼顾脱氮和除磷是A/A/O工艺的特征，脱氮效率受回流比的影响，无法实现完全脱除总氮，也存在着与除磷菌在碳源利用分配之间的矛盾。然而，前置好氧的O/A/O工艺因大幅度削减了毒性物质而有利于后续单元硝化菌的生长。与A/A/O工艺不同的是，该工艺不能利用废水中存在的易降解有机物作为碳源进行反硝化脱氮，造成一定程度的碳源浪费。由此可见，前置厌氧或者前置好氧对后续的脱氮工艺有着不同的影响机制，A/A/O工艺多用于生活污水^[5-6]，而O/A/O工艺可能更适合于工业废水^[7]。

焦化废水是典型的高碳氮比工业废水，含有多种高浓度有毒物质。其中的有机污染物主要包括酚类^[8]、苯系物、杂环芳烃和多环芳烃等物质^[9]；其无机物中，S^2-、SCN⁻、CN⁻等均为典型的毒性物质，并且对废水的COD值有较大的贡献^[10]。LI等^[11]研究了在相同水力停留时间下A/A/O与A/O工艺分别对焦化废水中COD和NH₄⁺-N的去除效果，发现两者的去除率几乎相同，但A/A/O比A/O工艺对总氮的去除效果更好。汤清泉等^[12]比较了A/A/O与O/A/O工艺对焦化废水的处理效果，认为碳氮比是决定二者对总氮去除效果的关键因素。当碳氮比为15~20时可以选择A/A/O工艺，当碳氮比为20~35时则O/A/O工艺效果更好。其原因是：前置好氧单元可以去除高碳氮废水中的有机物而降低后续处理的负荷。本课题组在长期实践的基础上开发了针对焦化废水处理三污泥系统的好氧-水解-好氧流化床脱氮工艺(命名为O/H/O工艺，其中，O1为除碳氨化单元，H为水解脱氮单元，O2为完全硝化单元)^[13-15]，已有 5个实际工程应用案例，最长运行时间达到12年。O/H/O工艺具有独特的三相分离器，可以保证在不需要污泥回流的情况下实现各个单元反应器独立的污泥特征和生物量，节省了能耗，并促进了污泥生态与水质环境的相容性^[16]。新型结构生物三相流化床作为O1反应器，在进水有机负荷达到2.4 kg ·(m³_·d)⁻¹ 的运行情况下，其耗氧有机物的去除率可以达到93.0%以上，反应器中氧的利用率为50%~60%。面对高毒性、高浓度的焦化废水，A/A/O工艺需要1~2倍稀释后才能进入生物系统，而O/A/O或O/H/O工艺则不需要稀释。

厌氧、水解、好氧单元不同顺序的排列组合构成了不同的废水生物处理工艺技术。在废水性质转化方面，厌氧单元可提高B/C值^[17]，而好氧单元可降低B/C值，分别有利于异养反硝化与硝化反应；在脱氮模式中，要考虑硝化反硝化^[18]、短程硝化反硝化^[19]、厌氧氨氧化^[20]、自养反硝化^[21]、好氧反硝化^[22]等原理的选用、协同及条件控制。A/A/O工艺和O/A/O工艺都需要回流才能保持反应器内的污泥浓度，A/A/O工艺的运行属于单污泥系统，O/A/O工艺中设置了2个二沉池，属于双污泥系统，而O/H/O工艺属于三污泥系统。根据废水的性质选择合适的工艺，可以在达标排放的基础上实现能耗与物耗的减量化。由于目前缺乏不同工艺特征的比较，为此，本文分析了不同工艺的碳源利用模式和脱氮模式，提出了一种代表性的焦化废水组成并通过研究A/A/O、O/A/O、O/H/O的组合工艺对焦化废水中核心污染物的去除及其能耗分配关系，阐明了工艺技术选择的原则，为复杂工业废水生物处理技术的工艺优选提供参考。

1. 研究方法

1.1 数据来源

本课题组对国内38个焦化厂进行了实地调查和数据采集，分析了焦化废水的水质特征与地域差异的关系，发现华北、华中、华东地区废水中的COD值略高，华中和西南地区废水的氨氮浓度略低^[23]。焦化废水中的含氮物质主要由氨氮、有机氮、SCN⁻、CN⁻等组分构成，由于蒸氨工艺的差异，含氮物质的比例各有不同。综合国内外的焦化废水原水水质^[24-26]，结合我们的调查，为了消除差异性和增强可比性，本文定义代表性的焦化废水组成为： COD为4 000 mg·L⁻¹，苯酚、 ${\rm{NH}}_4^{+}$ -N、SCN⁻、CN⁻、S²⁻以及总氮的质量浓度分别为800、 100、 500、 50、 50 和280 mg·L⁻¹。

A/A/O工艺借鉴宝武韶钢公司的运行数据，水量为60 m³·h⁻¹，3个单元反应器的水力停留时间分别为34、22和52 h，COD负荷分别为1.22、1.46和0.47 kg·(m³·d)⁻¹；O/H/O工艺参考实验室和焦化厂的运行数据^[27-28]，废水处理量为60 m³·h⁻¹，3个单元反应器的水力停留时间分别为36、40和24 h，COD负荷分别为2.30、0.38和0.55 kg·(m³·d)⁻¹；选取韩国某厂实验室数据作为O/A/O工艺的案例^[29]，实验规模为0.03 L·h⁻¹，3个反应器的水力停留时间分别为28.8、12和19.2 h，进水中添加KH₂PO₄和Na₂CO₃以维持碱度，在缺氧池中加入3倍总氮浓度的甲醇作为碳源，工艺装置总水力停留时间为2.5 d。通过实际与假设相结合的方法进行分析，以3个焦化厂的实际废水数据（见表1）来剖析不同工艺的碳源利用和脱氮模式。O/A/O和O/H/O工艺的反应器排列顺序相同，反应器的性能和运行模式不同。因此，在分析碳源利用和脱氮模式时只考虑A/A/O与O/A/O的对比，而在能耗分析时，再考虑O/A/O与O/H/O的差异性。

表 1 3种工艺实际运行水质

Table 1. Actual operating water quality in three processes mg·L⁻¹

工艺	COD	挥发酚	${\rm{NH}}_4^{+}$ -N	SCN⁻	CN⁻
A/A/O	1 727±60	742±69	173±12	175±18	26.2±4.5
O/A/O	2 300±100	635±15	235±15	375±25	-
O/H/O	3 451±215	973±74	245±15	450±17	25±3
注：以集水调池的水质作为生物上水。

| Show Table

DownLoad: CSV

1.2 过程分析

根据污染物的降解途径计算了污染物的COD当量和TN当量，结果见表2，在生物系统里，SCN⁻和CN⁻中的氮转化为氨氮^[30-31]。

表 2 不同污染物对COD和总氮的贡献

Table 2. Contribution of various pollutants to COD and nitrogen mg·mg^-1

当量	挥发酚	SCN⁻	CN⁻	S²⁻	${\rm{NO}}_3^{-}$	${\rm{NO}}_2^{-}$
COD当量	2.380	1.100	0.615	2.000	-	0.348
N当量	-	0.241	0.538	-	0.226	0.304

| Show Table

DownLoad: CSV

通过分析不同污染物对COD和总氮的贡献，检验废水组成的合理性。如式(1)所示，废水中的含氮量主要由 ${\rm{NO}}_3^{-}$ 、 ${\rm{NO}}_2^{-}$ 、 ${\rm{NH}}_4^{+}$ 、SCN⁻、CN⁻以及其他有机氮提供。如式(2)所示，废水中的COD主要由有机物和还原性无机离子构成，其中，挥发酚、苯系物、SCN⁻、S²⁻的贡献比例比较大，部分难降解的有机物也导致生物出水中检出较高的COD值。

$C_{\text{TN}}\text{= 0.226}{C}_{{\rm{NO}}_3^{-}}\text{+0.538}{C}_{{\text{CN}}^{-}}\text{+0.241}{C}_{{\text{SCN}}^{-}}\text+{C}_{{\rm{NH}}_4^{+}{-{\rm{N}}}}\text{+0.304}{C}_{{\text{NO}}_{\text{2}}^{-}}\text+{{C}}_{\text{其他含氮物质}}$

$(1)$

$C_{\text{COD}}\text{= 2.380}{C}_{\text{phenol}}\text{+1.100}{C}_{{\text{SCN}}^{-}}\text{+0.615}{C}_{{\text{CN}}^{-}}\text{+2.000}{C}_{{\text{S}}^{2-}}\text+{\text{0.348}C}_{{\text{NO}}_{\text{2}}^{-}}\text+{C}_{\text{其他有机物}}$

$(2)$

式中：C_TN、 $C_{{\rm{NO}}_3^{-}}$ 、 $C_{{\rm{NO}}_2^{-}}$ 、 $C_{CN^{-}}$ 、 $C_{SCN^{-}}$ 、 $C_{{\rm{NH}}_4^{+}}$ _-N、C_{其他含氮物质}分别表示总氮、硝酸根、亚硝酸根、氰化物、硫氰化物、氨氮以及其他含氮物质的质量浓度，mg·L⁻¹；C_COD为废水中耗氧有机物(以COD计)的质量浓度，mg·L⁻¹；C_phenol、C_S²⁻、C_{其他有机物}分别表示废水中苯酚、硫离子以及其他有机物的质量浓度，mg·L⁻¹。

在每一个单元反应器的出水中，都通过以上的方法进行检验，以确定废水组成的合理关系。

1.3 数据处理

根据式(3)~式(7)计算A/A/O工艺中每个反应器对污染物i总体去除的贡献率，分别以 $P_{{\rm{A1}}}^{\rm{i}}$ 、 $P_{{\rm{A2}}}^{\rm{i}}$ 、 $P_{{\rm{O}}}^{\rm{i}}$ 表示。根据式(8)~式(13)计算O/A/O每个反应器对污染物i总体去除率，分别以 $P_{{\rm{O1}}}^{\rm{i}}$ 、 $P_{{\rm{A}}}^{\rm{i}}$ 、 $P_{{\rm{O2}}}^{\rm{i}}$ 表示。

${P}_{\text{A1}}^{\text{i}}\text=\left(\text{1+}{R}_{\text{1}}\right)\times \frac{{C}_{\text{A1-I}}^{\text{i}}\text-{C}_{\text{A1-E}}^{\text{i}}}{{C}_{\text{0}}^{\text{i}}}\times 100\%$

$(3)$

${P}_{\text{A2}}^{\text{i}}\text=\left(\text{1+}{R}_{\text{1}}\text+{R}_{\text{2}}\right)\times \frac{{C}_{\text{A2-I}}^{\text{i}}\text-{C}_{\text{A2-E}}^{\text{i}}}{{C}_{\text{0}}^{\text{i}}}\times 100\%$

$(4)$

${P}_{\text{O}}^{\text{i}}\text=\left(\text{1+}{R}_{\text{1}}\text+{R}_{\text{2}}\right)\times \frac{{C}_{\text{O-I}}^{\text{i}}\text-{C}_{\text{O-E}}^{\text{i}}}{{C}_{\text{0}}^{\text{i}}}\times 100\%$

$(5)$

${\text{C}}_{\text{A2-I}}^{\text{i}}\text=\frac{{\text{(1+}{\text{R}}_{\text{1}}\text{)}\times \text{C}}_{\text{A1-E}}^{\text{i}}\text+{\text{R}}_{\text{1}}\times {\text{C}}_{\text{O-E}}^{\text{i}}}{\text{1+}{{\text{R}}_{\text{1}}\text+\text{R}}_{\text{2}}}$

$(6)$

${C}_{\text{O-I}}^{\text{i}}\text={C}_{\text{A2-E}}^{\text{i}}$

$(7)$

${P}_{\text{O1}}^{\text{i}}\text{=(1+}{R}_{\text{3}}\text{)}\frac{{C}_{\text{O1-I}}^{\text{i}}\text-{C}_{\text{O1-E}}^{\text{i}}}{{C}_{\text{0}}^{\text{i}}}\times 100\%$

$(8)$

${P}_{\text{A}}^{\text{i}}\text=\left(\text{1+}{R}_{\text{4}}\text+{R}_{\text{5}}\right)\times \frac{{C}_{\text{A-I}}^{\text{i}}\text-{C}_{\text{A-E}}^{\text{i}}}{{C}_{\text{0}}^{\text{i}}}\times 100\%$

$(9)$

$P_{\text{O2}}^{\text{i}}\text=\left(\text{1+}{R}_{\text{4}}\text+{R}_{\text{5}}\right)\times \frac{{C}_{\text{O2-I}}^{\text{i}}\text-{C}_{\text{O2-E}}^{\text{i}}}{{C}_{\text{0}}^{\text{i}}}\times 100\%$

$(10)$

${C}_{\text{O1-I}}^{\text{i}}\text=\frac{{C}_{\text{0}}^{\text{i}}\text+{R}_{\text{3}}{C}_{\text{O1-E}}^{\text{i}}}{\text{1+}{R}_{\text{3}}}$

$(11)$

$C_{\text{A-I}}^{\text{i}}\text=\frac{C_{\text{O1-E}}^{\text{i}}\text+{\text{(}R}_{\text{4}}\text+{R}_{\text{5}}\text{)}\times {C}_{\text{O2-E}}^{\text{i}}}{\text{1+}{{R}_{\text{4}}\text+R}_{\text{5}}}$

$(12)$

$C_{\text{O2-I}}^{\text{i}}\text={C}_{\text{A-E}}^{\text{i}}$

$(13)$

式中：i为各种污染物(COD、苯酚、硫氰化物、氰化物、氨氮、亚硝酸根、硝酸根和总氮)。R₁和R₂分别为A/A/O工艺中污泥回流比和硝化液回流比，污泥回流比取值1，硝化液回流比取值3；R₃、R₄、R₅分别为O/A/O工艺中初沉池回流至O1的污泥回流比、二沉池回流至A的污泥回流比以及硝化液回流比，均取值为1。C₀ⁱ为未处理废水中污染物i的质量浓度，mg·L⁻¹； $C_{{\rm{A1-I}}}^{\rm{i}}$ 和 $C_{{\rm{A1-E}}}^{\rm{i}}$ 、 $C_{{\rm{A2-I}}}^{\rm{i}}$ 和 $C_{{\rm{A2-E}}}^{\rm{i}}$ 、 $C_{{\rm{O-I}}}^{\rm{i}}$ 和 $C_{{\rm{O-E}}}^{\rm{i}}$ 、 $C_{{\rm{O1-I}}}^{\rm{i}}$ 和 $C_{{\rm{O1-E}}}^{\rm{i}}$ 、 $C_{{\rm{A-I}}}^{\rm{i}}$ 和 $C_{{\rm{A-E}}}^{\rm{i}}$ 、 $C_{{\rm{O2-I}}}^{\rm{i}}$ 和 $C_{{\rm{O2-E}}}^{\rm{i}}$ 分别表示A1、A2、O反应器和O1、A、O2反应器中污染物i的进水和出水的质量浓度，mg·L⁻¹。

排除水力停留时间对工艺对比造成差异，假设A/A/O与O/A/O工艺具有相同的总水力停留时间，结合文献调研和实际考虑，每个工艺各个反应器的体积比为1:1:2，处理水量为60 m³·h⁻¹。

污染物在反应器中会进行到氨化碳氧化、亚硝化氮氧化或硝化氮氧化3种不同的处理阶段，不同阶段的耗氧量分别根据式(14)~式(16)进行计算。

${O}_{\text{S}}\text=\left[a\cdot{K}_{\mathrm{C}\mathrm{O}\mathrm{D}}\cdot{C}_{\text{COD}}\text+{C}_{\text{DO}}\right]\frac{Q}{\text{24 000}}$

$(14)$

${O}_{\text{S}}\text=\left[a\cdot{K}_{\text{COD}}\cdot{C}_{\text{COD}}\text+b\left(\text{1}\text-{K}_{\text{d}}\right)\text{×}\left({C}_{\text{N}}\text+\frac{{C}_{\text{CN}}}{\text{1.86}}\text+\frac{{C}_{\text{SCN}}}{\text{4.14}}\right)\text+\left(\text{1+}{R}_{\text{S}}\text+{R}_{\mathrm{d}}\right)\cdot{C}_{\text{DO}}\right]\frac{Q}{\text{24 000}}$

$(15)$

$O_{\text{S}}\text=\left[a\cdot{K}_{\text{COD}}\cdot{C}_{\text{COD}}\text+c\left(\text{1}\text-{K}_{\text{d}}\right)\text{×}\left({C}_{\text{N}}\text+\frac{{C}_{\text{CN}}}{\text{1.86}}\text+\frac{{C}_{\text{SCN}}}{\text{4.14}}\right)\text+\left(\text{1+}{R}_{\mathrm{S}}\text+{R}_{\text{d}}\right)\cdot{C}_{\text{DO}}\right]\frac{Q}{\text{24 000}}$

$(16)$

${K}_{\text{d}}\text=\left(\text{1-}\frac{{N}_{\text{o}}}{{N}_{\mathrm{i}}}\right)\text{×}100\%$

$(17)$

式中：Q为生物系统进水量，m³·d⁻¹；a、b、c分别为氧化COD、氨氮到亚硝氮、氨氮到硝态氮的有关的耗氧系数，在本研究中为1.4、3.43、4.57；O_s为好氧单元的理论需氧量，kg·h⁻¹；C_COD为耗氧有机物(以COD计)的质量浓度，mg·L⁻¹；C_DO为好氧单元溶解氧的质量浓度，mg·L⁻¹；K_COD为COD去除率，%；R_s、R_d分别为活性污泥和硝化液回流比；C_N、C_CN、C_SCN分别为以脱氮为目标的好氧池中含氨氮、氰化物、硫氰化物的质量浓度，mg·L⁻¹；K_d为反硝化率，%；N_i、N_o分别为脱氮系统进、出水总氮的质量浓度，mg·L⁻¹。

A/A/O中的好氧单元主要发挥硝化作用，通过式(16)和式(17)计算可知其耗氧量；O/A/O工艺中，O1易氧化降解耗氧有机物(以COD计)，不考虑硝化作用，耗氧量通过式(14)计算可知；在O2中进行硝化作用，耗氧量通过式(16)和式(17)计算可知；O/H/O工艺与O/A/O工艺相似，但不需要污泥回流，因此，在计算O/H/O工艺中O2的曝气能耗时，式(16)的污泥回流比Rs为0。

污泥回流的能耗是A/A/O与O/A/O工艺所必不可少的，只有通过污泥回流才能保证生物池活性污泥的浓度，回流泵的能耗通过式(18)进行计算。

${W}_{\mathrm{S}}\text=K\cdot Q\cdot H$

$(18)$

$K\text=\frac{k\text{}}{\text{183.5}}$

$(19)$

式中：Ws为污泥回流泵的能耗，kW·h；K是安全系数，由式(19)计算，当水泵功率和污泥回流泵功率超过5 kW时，式(19)中的k取值1.15^[32]；Q为回流的流量，m³·d⁻¹；H为水泵总水头损失，m。

由于回流污泥含水率高达99.5%~99.9%，所以，污泥回流与废水回流的能耗以相同方法计算。A/A/O与O/A/O工艺的回流比已经明确，O/H/O工艺仅存在硝化液回流，回流比为1，能耗估算值可由泵能耗的公式给出。A/A/O工艺中污泥回流至厌氧池的水头损失为1.5~2 m，硝化液回流至缺氧池的水头损失为1~1.2 m。O/A/O工艺有2个污泥回流系统，二沉池至O1的水头损失为0.5~0.8 m，另一个二沉池至A的水头损失为1~1.5 m。O/H/O不存在污泥回流，硝化液回流的水头损失为1~1.6 m。

2. 结果与讨论

2.1 碳源利用

首先考察了2种工艺中COD的沿程变化，分析2种工艺的碳源利用模式差异。由图2可以看出，在A/A/O工艺中，O单元对耗氧有机物(以COD计)的去除效果最好，A1的水解作用使难降解有机物断链、开环，转化为小分子有机酸，为后续的反硝化脱氮所利用；而在O/A/O工艺中，O1对COD的去除率高达90.0%以上，使后续单元工艺主要为脱氮服务。两者不同的是，A/A/O工艺通过微生物反硝化作用去除了废水中的耗氧有机物，而O/A/O工艺则通过生物耗氧直接氧化废水中的耗氧有机物。

图 2 各工艺沿程COD的变化

Figure 2. COD changes along each process.

下载: 全尺寸图片幻灯片

LI等^[11]对比了A/A/O与A/O工艺的处理效果，指出2个工艺对于有机物和氨氮的去除效果几乎相同，但A/A/O工艺更有利于总氮脱除，这是因为A/A/O工艺设置了产酸阶段。CHAKRABORTY等^[33]发现，在A1中COD的去除率为5%~11%，CN⁻降解率为35.0%，没有发现苯酚降解的中间产物和甲烷的生成。王子兴等^[34]指出，在A/A/O-MBR工艺处理焦化废水的过程中，单个反应器COD去除率分别为9.2%、73.5%、14.7%；经过GC/MS检测分析，苯酚在A1中的降解率为26.7%，而含氮杂环化合物以及苯系物的去除率分别为49.5%和65.8%。此外，有研究^[35]表明，在A/A/O工艺中，A1单元去除污染物效果不明显，COD去除率低于10%；A2单元的COD去除率最高，尤其是易降解有机物在此阶段几乎全部被利用；在O单元中，利用异养微生物好氧氧化残留的有机物，CN⁻和SCN⁻在O2中也被彻底去除。SHARMA^[36]研究了厌氧、缺氧、好氧单个单元的处理效果时发现，好氧单元可去除83.3%的CN⁻ 和62.0%的COD；当加入氰化物后，好氧单元中COD的去除率下降到52.0%。由此可见，废水组成的复杂性会影响单组分的去除效果。马昕等^[37]采用O/A/O工艺处理焦化废水时发现，在O1停留时间为16 h时对COD的去除率达到75.0%，这与我们调查的工艺结果相似。由图2(b)可见，在O/A/O工艺中，O1对COD的去除率很高，浪费了部分有机碳源，而添加的外部有机碳源是造成A单元COD去除率降低的原因之一^[38]；另一方面，O1中的氨化过程可为O2提供良好的硝化环境。以上研究结果表明，2种工艺对废水中碳源的利用在原理上存在非常大的差异。

2.2 脱氮模式

脱氮的效果可通过协调碳源、电子供受体以及DO等因素来实现，故根据2种工艺中氨氮浓度沿流程变化来分析不同脱氮模式的有效性。由图3可见，虽然O/A/O工艺进水氨氮偏高，但出水氨氮却很低，在O2单元中已经彻底硝化。可见，前置好氧工艺可以为后续O2创造良好的硝化条件。A1去除了27.0%的氨氮，而O1去除了87.5%的氨氮，即在A1中仍然保留着较高浓度的氨氮，而在O1中氨氮几乎完全硝化，这与在进水中是否添加磷盐有关^[39]。O1、A1中氨氮浓度的变化以及微生物同化、有机氮氨化、氰化物及硫氰化物氨化等可以同时发生。在工程研究中发现，O1中还存在亚硝化和硝化的可能性^[17]。

图 3 各工艺沿程铵离子浓度变化

Figure 3. Change in NH₄⁺-N concentration along each process.

下载: 全尺寸图片幻灯片

焦化废水中的含氮物质除了铵离子/氨分子外，还有SCN⁻、CN⁻以及含氮有机物。ZHANG等^[40]发现，A/A/O中各个单元对氨氮的去除率分别为-2.5%、3%、97%，A1出水中氨氮升高的原因是其他含氮物质氨化作用所致。吕鹏飞等^[41]的研究表明，2种流化床工艺的前置厌氧单元对氨氮有少量的降解，氨氮去除率分别为18.1%和35.6%，体现出反应器对于处理效果的影响不同，流化床反应器面对复杂毒性废水比传统的沸腾床反应器表现出更好的耐毒性抑制作用。经过缺氧反应器A2后，氨氮浓度的变化主要有回流导致的直接稀释以及微生物降解的共同作用。GUI等研究了2个A/A/O系统，在硝化液回流比为200%的情况下，氨氮的质量浓度由250 mg·L⁻¹降低至80 mg·L^−1[42]。易欣怡等^[28]考察了O/H/O工艺的焦化废水处理，发现O1单元能够把氰化物、硫氰化物氧化为氨氮，有机氮全部氨化，从而造成O1出水氨氮浓度的升高；而在H单元中，环状含氮化合物通过水解作用可实现分子开环转变为氨氮，回流液中的硝态氮实现反硝化转变为氮气；接下来的O2单元能够将残余低价状态的含氮化合物转变为硝态氮，所以对氨氮的去除非常彻底。由于多种含氮物质之间具有不同价态转化机制，工艺中合理安排碳源进行脱氮，以及通过回流/超越或微生物功能调控实现总氮的彻底去除将是工艺理论中具有挑战性的研究方向。

2.3 能耗分析

1)各单元反应器的去除效率。能耗分配受工艺的单元反应器组合的影响。单元反应器的不同组合顺序可构成多样的生物处理工艺，前置好氧与厌氧工艺对同一种废水会产生不同的污染物去除效率，较优的工艺应该是在达标排放(即核心污染物去除)的基础上实现时间和空间上的减量化，还要降低二次污染。图4反映了A/A/O和O/A/O工艺污染物浓度的沿程变化。沿流程图中的百分比数据代表反应单元出水污染物浓度占进水中污染物浓度的比例。除了内部降解外，还要考虑因回流引起的反应器内污染物浓度的稀释作用。结合文献调查，综合实际情况，总结出代表性焦化废水典型污染物在单元反应器中的去除效率，如图5所示。其中，假设SCN⁻和CN⁻在O/A/O工艺的O1中完全氨化。

图 4 典型污染物含量沿各工艺流程的变化

Figure 4. The variation of typical pollutant content along each process.

下载: 全尺寸图片幻灯片

图 5 典型污染物在各工艺单元反应器中的降解率

Figure 5. The degradation rate of typical pollutants in each process unit reactor.

下载: 全尺寸图片幻灯片

2)不同工艺的能耗分配。废水中的污染物在不同工艺各单元反应器中的总体去除率如图6所示。A/A/O工艺对污染物的降解主要集中在O单元中，O/A/O工艺的降解则集中在O1单元中。这两者的差异反映了前置好氧工艺与前置厌氧工艺在曝气能耗上的差别。通过式(14)~式(16)计算，各工艺需氧单元的曝气量如图7所示。A/A/O工艺中O单元的需氧量为102.7 kg·h⁻¹，O/A/O中O1和O2的需氧量分别为260.8 kg·h⁻¹和35.1 kg·h⁻¹。由图7可看出，O/A/O工艺的O1大部分的曝气量是用来去除易降解有机物，因此，需氧量较高。但当废水中有机物的浓度很低时(当不考虑有机物耗氧时)，A/A/O工艺氧化含氮物质需氧量为100.4 kg·h⁻¹，O/A/O工艺氧化含氮类物质的需氧量为83.9 kg·h⁻¹。因此，对于脱氮性能，O/A/O工艺比A/A/O工艺能耗更高。这归因于：在O1中解除了SCN⁻、CN⁻等有毒物质对A反应器微生物的抑制作用，使得在A中降解的含氮物质相对较多，可以实现O2单元的低能耗硝化反应。因此，当废水中的耗氧有机物的预处理较为彻底时，前置好氧工艺可以实现低耗能高效率脱氮。O/H/O工艺在保留了O/A/O工艺优点的基础上，实现了反应器内部流态化的颗粒污泥特征，氧传质系数是一般活性污泥的2倍左右^[43]，因此，与O/A/O工艺相比，O/H/O工艺在耗氧量的节能方面更能体现出优势。本课题组根据多年的O/H/O运行经验数据统计得出，在仅考虑脱氮目标时，O/H/O工艺的需氧量约为53.26 kg·h⁻¹。

图 6 不同工艺单元反应器对各污染物的总体去除率

Figure 6. Overall removal rate of various pollutants in the unit reactor of different process.

下载: 全尺寸图片幻灯片

图 7 各工艺段的曝气量分配

Figure 7. Aeration distribution of each process section.

下载: 全尺寸图片幻灯片

由图4所示的计算可得出，在A/A/O工艺中，进入A2的废水COD为1 140.0 mg·L⁻¹，硝化液回流的硝酸根为84.4 mg·L⁻¹，在A2中主要去除总氮中的硝酸根，其余的氨氮、SCN⁻、CN⁻等含氮物质只是发生了少量的生物降解，经过A2可去除80.0 mg·L⁻¹左右的硝态氮，满足微生物生长的碳源需求量为723.2 g·m⁻³ (缺氧条件下C∶N∶P = 200∶5∶1)，因为废水中含有一定量的有机物，故实际可以供微生物利用的量约为540.0 g·m⁻³，需要外加碳源122.1 g·m⁻³ (以甲醇计)。在O/A/O工艺中，进入A单元的废水COD值为633.3 mg·L⁻¹，其总氮类型为硝酸根和氨氮，浓度分别为93.9 mg·L⁻¹和50.0 mg·L⁻¹，在A中降解90.0 mg·L⁻¹的硝态氮，满足微生物正常生长的碳源需求量约为813.6 g·m⁻³，进入A的废水中可降解有机物的含量约为83.3 g·m⁻³，不足的碳源需要从外部添加486.9 g·m⁻³(以甲醇计)。以上的讨论是在不考虑O/A/O工艺中有超越进水的情况，但在实际工程中，往往会使部分集水调节池中的出水以超越O1池的方式进入A池，这样既可以降低O1的曝气能耗，又可减少A单元的外部碳源的需求量。当超越1/3处理量的废水进入A单元时，O1的曝气量变为174.0 kg·h⁻¹ ，超越之后A单元进水的有机物浓度达到977.8 mg·L⁻¹，可供微生物利用的量约为427.8 mg·L⁻¹，因此，折合计算1 m³废水仅需要257.2 g的外加碳源，节省了229.7 g的外部碳源(以甲醇计)。可以看出，O/A/O系统的模式多样性，可以实现总氮的低能耗高效率去除。在实际运行的O/H/O工艺中，由于不需要污泥回流，每个反应器可以灵活调控，因此，O/H/O工艺比O/A/O工艺更容易实现厌氧氨氧化反应，并且可以利用FeS进行自养反硝化脱氮而节省能耗，故实际的O/H/O工艺的外部碳源需求约0~220 g·m⁻³，具体的需求量取决于厌氧氨氧化与自养反硝化的耦合性能^[44]。

污泥回流可以保证生物单元中的污泥浓度即生物量。通过式(18)和式(19)的计算，A/A/O工艺的污泥回流和硝化液回流的总能耗约为42.37 kW·h；O/A/O系统污泥回流与硝化液回流的总能耗约为23.55 kW·h；O/H/O系统只存在硝化液回流，回流能耗约为9.42 kW·h。除了曝气和回流的能耗外，考虑综合因素，3种工艺归纳为2大类：厌氧-缺氧-好氧以及好氧-水解/缺氧-好氧。由于反应器的设置不同，好氧-水解/缺氧-好氧工艺又可以分类为O/A/O和O/H/O，分化出二污泥法和三污泥法，反应器的类型决定了工艺的耗能。若只考虑生物阶段的处理，废水COD在3 000~4 000 mg·L⁻¹、铵离子质量浓度在100~200 mg·L⁻¹时，A/A/O的处理费用为6~8 元·t^−1[45,46]，O/A/O的处理费用为7~9 元·t^{−1 [47-49]}，而O/H/O流化床工艺的处理费用仅为4~5元·t⁻¹，体现了不同技术的成本差异。

2.4 性能比较

单元工艺的摆放顺序不仅决定了整体工艺运行的能耗，还会对冲击负荷、系统中微生物菌落和处理效果产生很大的影响。李国令等^[4]指出，热单胞菌属、脱氯单胞菌属是O/A/O工艺好氧池中的优势菌属；热单胞菌属、脱氯单胞菌属、球形红假单胞菌属是O/A/O工艺缺氧池中的优势菌属。WANG等也发现^[50]，热单胞菌属与硝酸盐还原酶基因呈正相关，对同时厌氧氨氧化-反硝化系统中的硝酸盐还原起重要作用。 WEI等^[15]指出，丛毛单胞菌属在反应器O1中对COD去除起到了关键作用，有助于去除O1反应器中的NH₄⁺-N；硫杆菌则在H反应器中起着主要的反硝化作用，AOB和NOB(亚硝化单胞菌和硝化螺菌)对反应器中硝化作用的贡献最大。三污泥法的O/H/O工艺各单元在污染物组成、去除、功能和微生物群落等方面存在显著进步，有望实现厌氧氨氧化脱氮与深度脱氮的结合，也表明废水水质和反应器的组合对微生物功能分布具有调控功能。

根据污泥回流的设置与否，A/A/O、O/A/O、O/H/O工艺可以分为单污泥系统、双污泥系统及三污泥系统，3个工艺的主要区别见表3。据报道，A/A/O工艺中A1单元对COD去除效率小于10%，检测不到甲烷的产生^[51]。因此，A/A/O工艺仅仅在缺氧和有氧反应器中实现了对COD的去除。由于回流的存在，A/A/O工艺表现为单污泥特征，异养细菌具有较高的比生长速率，因污泥排放量高而导致其在处理高COD/TN废水时，大量自养硝化细菌被排洗。前置好氧工艺对高浓度毒性废水有很好的抗负荷冲击能力，并且O/H/O工艺中的新型结构流化床反应器的强化传质功能与污泥原位分离原理加强了各单元反应器中的微生物能力^[22]。在H单元中，根据投加的电子供体不同而具有多种反硝化模式：如利用O1池的剩余COD作为碳源及其他电子供体进行异养反硝化脱氮；通过投加无机还原性电子供体以利用其作为营养源进行自养反硝化脱氮^[21,52]，还可以避免二次碳源的污染。另外，有研究表明，控制O1反应器在短程硝化水平，可使亚硝酸盐直接得到富集和积累，然后实现厌氧氨氧化模式脱氮，从而使工艺过程节能效果更好^[17,19]。可见，复杂废水的脱氮模式多种多样，需要根据实际情况合理选择或耦合新原理，从而进一步实现低能耗、低物耗目标下的总氮去除。

表 3 不同工艺系统的特点

Table 3. Characteristics of different process systems

工艺	污泥系统	毒性物质的去除	COD/TN	脱氮途径	能耗影响因素	平均运行单价/(元·m^-3)	优点	缺点
A/A/O	单污泥系统	A1对大分子有机物的去除	11.4	异养反硝化	一次回流、一次曝气	7	有利于含氮有机物的水解；反硝化可利用废水中有机物作为碳源	不耐冲击负荷，受毒性抑制，需要稀释进水
O/A/O	双污泥系统	O1对SCN⁻、CN⁻的去除及氨化	12.5	异养反硝化、自养反硝化	二次回流、二次曝气	8	耐冲击负荷，进水不需要稀释；硝化效果好	耗氧量大，污泥回流频繁，耗能多
O/H/O	三污泥系统	O1对SCN⁻、CN⁻的去除及氨化	13.8	异养反硝化、自养反硝化、厌氧氨氧化及其耦合脱氮	二次曝气	4.5	耐冲击负荷，颗粒污泥耐毒性抑制，硝化效果好，不需要沉淀池；不需要回流	耗氧量大

| Show Table

DownLoad: CSV

3. 结论与展望

1)每处理1 m³设定浓度的焦化废水(不考虑O/A/O的超越进水)，A/A/O和O/A/O工艺分别需要122.1 g 与486.9 g的外部碳源(以甲醇计)。当废水中的易降解有机物较少且只考虑脱氮目标时，O/A/O工艺的曝气需氧量为83.9 kg·h⁻¹，A/A/O工艺的曝气需氧量为100.4 kg·h⁻¹；当O/A/O工艺中有1/3的进水流量超越至A单元时，其碳源需求量由486.9 g·m⁻³减至257.2 g·m⁻³(以甲醇计)，曝气量也将显著降低。

2)由于废水组成的复杂性，污染物的降解效率除了受到彼此的相互制约外，工艺条件和反应器的设计也至关重要。具有高毒性、高碳氮含量的焦化废水，更适合于选择前置好氧的工艺。O/H/O工艺由于其独特的三相分离器的设置而节省了污泥回流部分的能耗，反应器中的颗粒污泥更加耐毒性抑制和抗冲击负荷，并且传氧速率高，工艺耗氧量仅为53.26 kg·h⁻¹，外部碳源的消耗可以由486.9 g·m⁻³降至0~220 g·m⁻³。

3)反应器的高效性和可控性，使O/H/O工艺比O/A/O工艺更容易实现自养反硝化与异养反硝化协同脱氮、自养型短程反硝化与厌氧氨氧化的协同脱氮等其他脱氮途径，进而使O/H/O工艺成为一种更具潜力的低能耗、低物耗的生物脱氮技术工艺。针对不同的废水水质与物质组成特征，O/H/O工艺能够对不同功能的单元进行组合和编辑，从时间与空间、药剂与能耗、处理效率等方面追求更加丰富的优化模式，以满足各种不同的出水需求，特别是满足总氮浓度趋零的要求。

图 1 普通电阻率测井测量原理图

Figure 1. Schematic diagram of common resistivity logging

下载: 全尺寸图片幻灯片

图 2 改进的高密度电阻率电极杆与实物图

Figure 2. Improved high density resistivity electrode rod and physical diagram

下载: 全尺寸图片幻灯片

图 3 VMOD模拟堵塞三维立体图

Figure 3. 3D stereogram of constructed wetland clogging

下载: 全尺寸图片幻灯片

图 4 人工湿地堵塞流场拟合图

Figure 4. Constructed wetland clogging flow field fitting

下载: 全尺寸图片幻灯片

图 5 供电极距D_AB/2=7.5 cm的视电阻率断面等值线图

Figure 5. Apparent resistivity cross-section contour map at the supply pole distance D_AB/2=7.5 cm

下载: 全尺寸图片幻灯片

图 6 供电极距D_AB/2=12.5 cm的视电阻率断面等值线图

Figure 6. Apparent resistivity cross-section contour map at the supply pole distance D_AB/2=12.5 cm

下载: 全尺寸图片幻灯片

图 7 供电极距D_AB/2=22.5 cm的视电阻率断面等值线图

Figure 7. Apparent resistivity cross-section contour map at the supply pole distance D_AB/2=22.5 cm

下载: 全尺寸图片幻灯片

图 8 供电极距D_AB/2=7.5 cm的视电阻率三维立体图

Figure 8. 3D stereogram of apparent resistivity at the supply pole distance D_AB/2=7.5 cm

下载: 全尺寸图片幻灯片

[1]	CORBELLA C, GARCÍA J, PUIGAGUT J. Microbial fuel cells for clogging assessment in constructed wetlands[J]. Science of the Total Environment, 2016, 569: 1060-1063.
[2]	MATOS M P, SPERLING M V, MATOS A T, et al. Clogging in constructed wetlands: Estimation of medium porosity by ground-penetrating radar[C]// 15^th International conference on wetland systems for water Pollution Control. London: IWA Publishing, 2016: 4-9.
[3]	COOPER D, GRIFFIN P, COOPER P. Factors affecting the longevity of sub-surface horizontal flow systems operating as tertiary treatment for sewage effluent[J]. Water Science and Technology, 2005, 51(9): 127-135. doi: 10.2166/wst.2005.0303
[4]	NIVALA J, KNOWLES P, DOTRO G, et al. Clogging in subsurface-flow treatment wetlands: Measurement, modeling and management[J]. Water Research, 2012, 46(6): 1625-1640. doi: 10.1016/j.watres.2011.12.051
[5]	MATOS M P, SPERLING M V, MATOS A T, et al. Clogging in constructed wetlands: Indirect estimation of medium porosity by analysis of ground-penetrating radar images[J]. Science of the Total Environment, 2019, 676: 333-342. doi: 10.1016/j.scitotenv.2019.04.168
[6]	TAPIAS J C, HIMI M, LOVERA R, et al. Assessing clogging development in infiltration-percolation systems for wastewater treatment by electrical resistivity and induced polarisation methods[J]. //EGU General Assembly 2013. EGU General Assembly Conference Abstracts. Austria, 2013: EGU2013-8050.
[7]	张明珍, 徐栋, 武俊梅, 等. 人工湿地堵塞监测方法的研究进展[J]. 水生态学杂志, 2021, 42(3): 121-126.
[8]	MORVANNOU A. Dynamic modelling of nitrification in vertical flow constructed wetlands[D]. Brussels: Université Catholique de Louvain, 2012.
[9]	MORRIS R H, NEWTON M I, KNOWLES P R, et al. Analysis of clogging in constructed wetlands using magnetic resonance[J]. The Analyst, 2011, 136(11): 2283-2286. doi: 10.1039/c0an00986e
[10]	LAURENT J, BOIS P, NUEL M, et al. Systemic models of full-scale surface flow treatment wetlands: Determination by application of fluorescent tracers[J]. Chemical Engineering Journal, 2015, 264: 389-398. doi: 10.1016/j.cej.2014.11.073
[11]	李曼, 华国芬, 姬雨雨. 垂直流人工湿地堵塞和轮休过程中水流变化规律[J]. 水处理技术, 2017, 43(1): 57-61.
[12]	MARTÍN I, SALAS J J, SARDÓN N, et al. Research on clogging and effective surface in subsurface flow constructed wetlands[C]// 12^th International conference on wetland systems for water Pollution Control. London: IWA Publishing, 2010: 4-8.
[13]	SMITH E, GORDON R, MADANI A, et al. Cold climate hydrological flow characteristics of constructed wetlands[J]. Canadian Biosystems Engineering, 2005, 47: 1.1-1.7.
[14]	WHITMER S, BAKER L, WASS R. Loss of bromide in a wetland tracer experiment[J]. Journal of Environmental Quality, 2000, 29(6): 2043-2045.
[15]	张亮, 曹营渠, 易畅, 等. 人工湿地堵塞的研究进展[J]. 四川环境, 2017, 36(S1): 174-178.
[16]	AIELLO R, BAGARELLO V, BARBAGALLO S, et al. Evaluation of clogging in full-scale subsurface flow constructed wetlands[J]. Ecological Engineering, 2016, 95: 505-513. doi: 10.1016/j.ecoleng.2016.06.113
[17]	楚泽涵. 地球物理测井方法与原理[M]. 北京: 石油工业出版社出版, 2007: 26-42.
[18]	程志平. 电法勘探教程[M]. 北京: 冶金工业出版社, 2007: 13-20.
[19]	ARCHIE G E. The electrical resistivity log as an aid in determining some reservoir characteristics[J]. Transactions of the AIME, 1942, 146(3): 54-62.
[20]	丁彦礼, 李为一, 李文龙, 等. 一种小型可监测堵塞的人工湿地: 201820973071.3[P]. 2019-04-12.
[21]	王帅兵, 张雪蕾, 李常斌, 等. 基于Visual MODFLOW的武威盆地地下水动态模拟及预测[J]. 兰州大学学报(自然科学版), 2015, 51(2): 159-165.
[22]	郝治福, 康绍忠. 地下水系统数值模拟的研究现状和发展趋势[J]. 水利水电科技进展, 2006, 26(1): 77-81. doi: 10.3880/j.issn.1006-7647.2006.01.023
[23]	王庆永, 贾忠华, 刘晓峰, 等. Visual MODFLOW及其在地下水模拟中的应用[J]. 水资源与水工程学报, 2007, 18(5): 90-92. doi: 10.3969/j.issn.1672-643X.2007.05.023
[24]	何杉. Processing Modflow软件在地下水污染防治中的应用[J]. 水资源保护, 1999, 57(3): 16-18.
[25]	白少元, 宋志鑫, 丁彦礼, 等. 潜流人工湿地基质结构与水力特性相关性研究[J]. 环境科学, 2014, 35(2): 592-596.
[26]	贺映全, 曹红军, 胡武林, 等. 垂直流人工湿地基质堵塞分析与处理措施[J]. 山西建筑, 2019, 45(10): 175-176. doi: 10.3969/j.issn.1009-6825.2019.10.093
[27]	丁彦礼, 李为一, 韦益华, 等. 一种强化电阻率差异的人工湿地堵塞区域探测方法: 201810656942.3[P]. 2018-11-16.

点击查看大图

图( 8)

计量

文章访问数: 3673
HTML全文浏览数: 3673
PDF下载数: 87
施引文献: 0

全文HTML

利用人工湿地进行污水处理的技术目前已经得到广泛应用，但人工湿地运行过程中频繁出现的堵塞问题，已严重影响到人工湿地的持久和高效运行。而目前对人工湿地发生堵塞的判断及堵塞程度仅能进行定性评价，无法对堵塞区域进行精确定位。因此，在治理堵塞的时候没有具体针对性，在治理的时候只能针对于人工湿地整体进行，时间成本高、经济效果差。因此，针对人工湿地堵塞区域的定位探测是亟待解决的问题。

目前在实际针对人工湿地堵塞问题上，较为常用的有生物电池法^[1]、探地雷达法^[2-5] 、电阻率法^[6-8]、渗透系数法^[9]、示踪剂法^[10-11]、分析堵塞物质性质^[12]等方法。对于示踪剂法，投入的示踪剂会被湿地植物和其他生物吸收^[13-14]，且单一的示踪实验结果可能不具有代表性^[4]。通过水力传导率法测得的值仅表示轴流方向上横截面内的平均水力传导率，并不能表示该横截面内某特定垂直和横向位置处的堵塞严重程度^[15]，实际操作过程中一般将水力传导率法和示踪剂法相结合。AIELLO等^[16]分析了水平潜流人工湿地的水力特性，通过现场测量砾石层的水力传导率、量化累积堵塞物质，最后通过示踪实验可视化流动路径来研究堵塞现象。对于利用地球物理方法探测人工湿地堵塞问题方面，目前的文献主要还是集中在电阻率法和探地雷达方法，探地雷达精度较高、速度快，适合中、小尺度监测，可得到湿地内部能量衰减图像^[3-4]。特别是MATOS等^[5]使用探地雷达探测了种植香蒲与未种植植被的2个全尺寸水平潜流人工湿地的堵塞特性，结果表明，探地雷达可以探测出堵塞区域，但地质雷达方法也存在纵向尺度定位不精确的问题，并且由于人工湿地一般纵向尺度较小和水饱和状态，存在基底和边界反射较强、堵塞区域反射不够清晰的问题，因而探测精度不够理想。电阻率法对湿地内部造成的干扰较小，可以快速、无损、有效地反映湿地堵塞状况，在湿地堵塞物定性定量方面具有良好的发展前景^[7-8]。但目前使用的电阻率探测方法多采用的是地面电阻率法，在纵向尺度探测效果不佳，并且干扰因素过多。

基于上述情况，本研究以较为广泛的潜流人工湿地为研究对象，根据潜流人工湿地结构和堵塞区域的特点，基于电阻率法提出了利用改进的高密度电阻率法来探测和定位人工湿地堵塞区域，并使用Visual MODFLOW(VMOD)建立了堵塞模型，主要从人工湿地堵塞探测的电阻率方法讨论、电阻率法的水槽模拟实验和流场模拟几个方面，探讨了基于电阻率法探测和定位人工湿地堵塞区域的可行性。

3. 结论

1)只要电阻率差异存在并且差异够大，改进的高密度电阻率测井方法能够较好地探测和定位人工湿地高阻堵塞区域。

2)改进的高密度测井方法采用的数据供电极距不宜过大，如果供电极距过大的话，不但有效测量范围会减小，而且异常效果由于受边界和电极距范围内介质综合效应的影响而变差。根据实验结果，建议高密度电阻率探杆的电极间距不大于5 cm，供电极距D_AB/2不超过12.5 cm。

3)对于视电阻率的计算方面，采用全空间的视电阻率计算公式计算得到的视电阻率，除了边界附近存在较小的误差外，其他区域误差相对较小。因此，在方法探索阶段，可以采用全空间的视电阻率计算公式计算视电阻率。但后续还需进一步研究边界处全空间条件和半空间条件的转换问题，在视电阻率的计算上需要研究全空间和半空间的混合算法，以期消除边界处的计算误差，增强边界处的分辨率。

4)对于人工湿地堵塞探测的数据处理来说，根据探测数据绘制三维立体图，可以更好的展示高阻堵塞体的空间位置和形态，更精确的实现堵塞区域的定位。

5)对于人工湿地工程实例来说，人工湿地未堵塞区域和堵塞区域的电阻率差异并不是很大，并不能保证有较好的探测效果，可以通过在人工湿地中加入电解质溶液的方法来加大人工湿地未堵塞区域和堵塞区域的电阻率差异，提升人工湿地堵塞区域探测精度和探测效果。

参考文献 (27)

基于电阻率法的人工湿地堵塞区域探测方法

作者简介: 丁彦礼(1979—)，男，硕士，副教授，dyl@glut.edu.cn

通讯作者: 白少元(1981—)，女，博士，教授，baisy@glut.edu.cn;