CN114906906A - 一种高效处理有机废水的梯级电催化氧化工艺设计方法 - Google Patents

一种高效处理有机废水的梯级电催化氧化工艺设计方法 Download PDF

Info

Publication number
CN114906906A
CN114906906A CN202210529802.6A CN202210529802A CN114906906A CN 114906906 A CN114906906 A CN 114906906A CN 202210529802 A CN202210529802 A CN 202210529802A CN 114906906 A CN114906906 A CN 114906906A
Authority
CN
China
Prior art keywords
reactor
electrode
value
organic wastewater
oxidation process
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202210529802.6A
Other languages
English (en)
Inventor
王立章
侯静
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Mining and Technology CUMT
Original Assignee
China University of Mining and Technology CUMT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Mining and Technology CUMT filed Critical China University of Mining and Technology CUMT
Priority to CN202210529802.6A priority Critical patent/CN114906906A/zh
Publication of CN114906906A publication Critical patent/CN114906906A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • C—CHEMISTRY; METALLURGY
    • C02—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F1/00—Treatment of water, waste water, or sewage
    • C02F1/46—Treatment of water, waste water, or sewage by electrochemical methods
    • C02F1/461—Treatment of water, waste water, or sewage by electrochemical methods by electrolysis
    • C02F1/467—Treatment of water, waste water, or sewage by electrochemical methods by electrolysis by electrochemical disinfection; by electrooxydation or by electroreduction
    • C02F1/4672—Treatment of water, waste water, or sewage by electrochemical methods by electrolysis by electrochemical disinfection; by electrooxydation or by electroreduction by electrooxydation
    • C—CHEMISTRY; METALLURGY
    • C02—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F1/00—Treatment of water, waste water, or sewage
    • C02F1/46—Treatment of water, waste water, or sewage by electrochemical methods
    • C02F1/461—Treatment of water, waste water, or sewage by electrochemical methods by electrolysis
    • C02F1/46104—Devices therefor; Their operating or servicing
    • C02F1/46109—Electrodes
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • C—CHEMISTRY; METALLURGY
    • C02—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F2101/00—Nature of the contaminant
    • C02F2101/30—Organic compounds
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02W—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO WASTEWATER TREATMENT OR WASTE MANAGEMENT
    • Y02W10/00—Technologies for wastewater treatment
    • Y02W10/10—Biological treatment of water, waste water, or sewage

Landscapes

  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Water Supply & Treatment (AREA)
  • Organic Chemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Hydrology & Water Resources (AREA)
  • Electrochemistry (AREA)
  • Environmental & Geological Engineering (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • General Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Water Treatment By Electricity Or Magnetism (AREA)

Abstract

本发明公开了一种高效处理有机废水的梯级电催化氧化工艺设计方法，包括：S1、以DSA电极作为阳极，Ti作为阴极，填充颗粒作为粒子电极，组成多个梯级反应器串联的填充床电极反应器；S2、将有机废水以动态进水方式通入填充床电极反应器，接通电源使阳极、阴极通电，通过各梯级反应器对有机废水进行有机污染物处理；S3、有机废水有机污染物处理期间，对梯级电催化氧化工艺进行修正及调整。本发明通过有机物降解效能强化模型能够为梯级反应器的精确设计和优化操作提供科学而实用的计算工具，也能为电化学水处理工艺的工程应用奠定重要的理论和技术基础。

Description

一种高效处理有机废水的梯级电催化氧化工艺设计方法
技术领域
本发明涉及有机废水处理技术领域，具体是涉及一种高效处理有机废水的梯级电催化氧化工艺设计方法。
背景技术
梯级电催化氧化技术，是一种随降解过程中体系内有机质浓度的降低而逐级递减操作电流密度，以保障体系内氧化进程始终处于反应控制状态的一种新型电催化氧化技术。其核心在于电催化氧化反应不再以恒定的电流密度，在同一反应器内完成，而是多个串联运行的梯级反应器在极限电流密度下实现对有机质的逐级降解。
在恒电流填充床反应体系中，随电催化氧化反应的进行，体系内有机质逐渐被分解，COD逐渐降低。而梯级电催化氧化技术中，采用多个串联运行的梯级反应器，故每一级反应器都可以依据COD的变化情况而调整操作电流密度，因为操作电流密度与极限电流密度之间的关系可以控制有机物氧化进程控制步骤，进而保证每级反应器均处于反应控制阶段，进而达到低耗、高效处理有机质的目的。研究结果表明，相同操作条件下，梯级电催化氧化技术与传统的恒电流填充床反应器相比电耗可降低至少1/3以上，且反应器分级越多节能效果越好，但过多的梯级数量将会大幅增加基建费用和操作难度，导致效能强化目标难以实现。
而目前国内对梯级电催化氧化技术研究较少，缺乏系统的搭建方式。因此，急切需要建立同时具备普遍性和精确性的梯级操作电流密度的确定方法，明确梯级电催化氧化技术的操作流程，为梯级电催化氧化技术实际应用中的稳定运行及能耗预估奠定理论和实践基础。
发明内容
为解决上述技术问题，本发明提供了一种高效处理有机废水的梯级电催化氧化工艺设计方法。
本发明的技术方案是：一种高效处理有机废水的梯级电催化氧化工艺设计方法，包括以下步骤：
S1、以DSA电极作为阳极，Ti作为阴极，填充颗粒作为粒子电极，组成多个梯级反应器串联的填充床电极反应器；
S2、将有机废水以动态进水方式通入填充床电极反应器，接通电源使阳极、阴极通电，通过各梯级反应器对有机废水进行有机污染物处理；
S3、有机废水有机污染物处理期间，对梯级电催化氧化工艺进行修正及调整；
S3-1、通过修正电极电流密度的定量函数模型，以及反映多项体系耦合传质机制的实时传质模型，对梯级电催化氧化工艺施加的操作电流密度、传质系数进行修正与调整；
S3-2、通过有机物降解效能强化模型调整填充床电极反应器的梯级反应器数量，并对填充床电极反应器的影响参数进行效能预测与优化。
进一步地，所述DSA电极采用刷涂热氧化法或电沉积法制备，如：PbO2/Ti电极。
进一步地，所述填充颗粒为活性炭、涂膜活性炭中的任意一种或两种任意比的混合物。
更进一步地，所述活性炭、涂膜活性炭的粒径为5mm，涂膜活性炭采用浸泡、低温烘干定型方式制取。
进一步地，步骤S3-1中用于修正电极电流密度的所述定量函数模型，具体为：
利用CV曲线计算伏安电量，建立的q、qin和qout之间的关系拟合数值，用式(1)、(2)、(3)计算。以式(4)计算粒子电极相对阳极的活化面积拓展系数(λ)；构建等效电路，拟合EIS得到有机物氧化电阻(Rct)和析氧电阻(R0)，用式(5)计算交换电流密度i0并与极化曲线所得数值进行比对以校核EIS等效电路模型，而后据式(6)计算参与有机废水氧化的阳极和粒子电极有效电流占比系数γ、β，再应用填充床电极反应器在反应控制(零级反应)状态下所得的有机废水中有机物降解实时数据并根据式(7)进行γ、β的修正；
q＝qin+qout (1)
q(v)＝qout+k1v-1/2 (2)
q(v)-1＝q-1+k2v1/2 (3)
λ＝qout,P/qout,A (4)
i0＝RT/(NFRct) (5)
γ或(λ+β)＝R0/(R0+Rct) (6)
Figure BDA0003645722630000031
其中，q为体系的总伏安电量；qin、qout分别为体系的内层、外层伏安电量；q(v)为CV曲线的积分面积；v为扫速；k1、k2为常数；qout,A、qout,P分别为阳极、粒子电极的外层伏安电量；R为气体常数；T为绝对温度；c0、c(t)分别为填充床电极反应器进水、处理出水有机物浓度；ε为床层空隙率；x0为极板间距；N为电子转移数量；F为法拉第常数；α为电流密度操作系数；km为传质系数；t为反应时间；c(t)为t时的填充床电极反应器处理出水有机物浓度。
进一步地，步骤S3-1中用于反映多项体系耦合传质机制的所述实时传质模型，具体为：
结合雷诺数Re和施密特数Sc之间的关系并考虑梯级反应器几何构型、周向返混、液相特征等因素的影响，引入操作电流密度iR(以阳极几何面积计算)描述反应物实时传质进程：
Figure BDA0003645722630000041
其中，D为扩散系数；c0为反应物起始浓度；dp为填料平均粒径；RH为梯级反应器水力半径；y0为主反应区高度；χ为梯级反应器形状因子校正系数；ε为床层空隙率；N为电子转移数量；F为法拉第常数；
同时，利用在扩散控制(一级反应)状态下所得的有机废水有机物降解实时数据进行传质系数的修正，关联式如下所示：
Figure BDA0003645722630000042
其中，c(t)为t时的废水污染物浓度；c0为原水污染物浓度；λ为床层拓展系数；km为传质系数；ε为床层空隙率；x0为极板间距；t为反应时间；
利用式(9)所得的传质系数数值代入式(8)，校正常数项可获得反映多项体系耦合传质机制的实时传质模型。
在填充床电极反应器工作时，其内部传质路径、有效电流数值、电极活化面积与水力学及电化学条件、废水物理属性、电极类型、反应器构型的耦合机制密切关联，它们之间的关系很难利用准确、完整的数学公式进行表达，而这也关系到电催化氧化进程控制步骤的判定，而通过本发明定量函数模型以及实时传质模型的建立，既能反映多维场复杂耦合实际情况又能精确表述有机污染物实时传质系数和操作电流密度数值的定量函数关系。
进一步地，步骤S3-2中所述有机物降解效能强化模型具体为：
填充床电极反应器在反应控制状态时经济因子之间的关系采用式(10)表示：
Figure BDA0003645722630000051
其中，p1、p2、p3表示经济指标因子，在电极和填充床电极反应器主体材料已定时为常数；Q为废水流量；b为直流电源经济因子；I为电流；Esp为电耗；gE为电价；
梯级操作流程的效能强化模型f(αn)为各梯级单元fn(αn,Xn)的总和；X已定时，对于基建费用(式10前两项)和处理成本(式10最后一项)有相反的影响，因此f(αn)存在极小值；由于∑Xn＝X(n＝1,2,3,......)，且Xn与操作电流密度αn、fn(αn,Xn)密切相关；f(αn)对第n级梯级反应器进水或出水口处(命名为节点，其数量为n+1)的有机物无量纲浓度Cj(Cj＝cj/c0)求导并代入有关的关联式可得以下方程：
Figure BDA0003645722630000052
其中，
Figure BDA0003645722630000053
CFn为第n级梯级反应器进水有机物无量纲浓度；γ为阳极对有机物的氧化选择系数；
Figure BDA0003645722630000056
为节点虚去除率，表示处理流程的年费用折算值随节点归一化浓度增加而增加的比率；
Figure BDA0003645722630000057
为第n级梯级反应器的年费用折算值随节点归一化浓度增加而增加的比率。
式(11)所示的关系为标准的遗传算法目标函数表达式，可采用以下步骤求解各梯级反应器的有机污染物去除率
Figure BDA0003645722630000054
在假定n值的前提下，拟定各节点有机物无量纲浓度初值
Figure BDA0003645722630000055
其中C1＝1、Cn+1＝1-X，其与节点无量纲浓度均在1和(1-X)之间取值；在该初值点处采用Taylor公式将
Figure BDA0003645722630000058
展开并舍却高阶项可得到遗传算子：
Figure BDA0003645722630000061
其中，θ＝γ2(p1x0n+p2)/[(γ+β)km,n]；x0n为第n级梯级反应器的极板间距；km,n为第n级梯级反应器的传质系数；ΔCj为各节点的归一化浓度初值微量；
Figure BDA0003645722630000062
为节点虚去除率，表示处理流程的初始年费用折算值随节点归一化浓度增加而增加的比率；
Figure BDA0003645722630000063
为第n级梯级反应器的初始年费用折算值随节点归一化浓度增加而增加的比率；
Figure BDA0003645722630000064
为第n级梯级反应器进水有机物初始无量纲浓度；对于各节点Cj (0)反复施加应用式(12)计算的微量ΔCj，则梯级操作流程的f(αn)数值将不断减小，各节点满足式(11)的条件时取得极小值，使得f(αn)为最小值时的n值即为最佳梯级数量，进而可据Xn确定各梯级反应器几何尺寸与运行条件。
反应控制状态时，传质系数和操作电流密度在设定有机废水有机物去除率情况下可用于各梯级反应器的电耗、极板耗量、反应器大小等经济因子的确定，在由若干梯级反应器串联组成的梯级电催化氧化工艺中，各梯级反应器去除率相互约束，但其总值与所需有机物去除率相等；各梯级反应器有机物去除率组合情况不同时，梯级电催化氧化工艺的经济因子数值也随之发生变化，通过本发明构建的效能强化模型可对梯级电催化氧化工艺流程优化，从而指导梯级电催化氧化工艺流程的优化构建。
更进一步地，采用不同的电极材料，在极板间距数值允许范围之内分别计算f(αn)的数值，f(αn)的数值越小则所选用的电极材料和极板间距更优。
进一步地，所述步骤S3-2中影响参数包括：阳极材质、填料类型、极板间距、床层厚度与高度以及粒径与床层高度的比例、进水流速以及填料粒径。
本发明的有益效果是：
(1)本发明引入操作电流密度的影响构建多相耦合体系有机污染物实时传质模型，通过定量函数模型以及实时传质模型的建立，既能反映多维场复杂耦合实际情况又能精确表述有机污染物实时传质系数和操作电流密度数值的定量函数关系。
(2)本发明提出了一种以伏安电量、EIS等效电路模型数值分别定量表述电极活化面积与有效电流的新方法。
(3)本发明构建了一种采用极化性能所得的极限电流密度校核EIS等效电路准确性的新方法。
(4)本发明建立了一种用于梯级电催化氧化工艺的有机物降解效能强化模型，可用于经济因子及有机废水处理效能的精确预测与优化。
(5)本发明为梯级电催化氧化工艺的优化提供了新的方式，通过有机物降解效能强化模型能够为梯级反应器的精确设计和优化操作提供科学而实用的计算工具，也能为电化学水处理工艺的工程应用奠定重要的理论和技术基础。
附图说明
图1梯级填充床电催化氧化工艺实验装置(A)以及反应机理图(B)，1-8分别代表进水、计量泵、取样阀、流量计、储液槽、搅拌器、双电极体系、出水。
图2恒电流与梯级体系在不同苯酚初始浓度下CODt/COD0与反应时间的关系图。
图3恒电流与梯级体系在不同苯酚初始浓度下Yst与反应时间的关系图。
图4恒电流与梯级体系在不同苯酚初始浓度下Esp与反应时间的关系图。
图5恒电流与梯级体系在不同苯酚初始浓度下ACE与反应时间的关系图。
图6梯级电催化氧化模式示意图，阴影区描绘了实际的能耗。
图7梯级电催化氧化处理过程中所需电流步数的优化结果(A)以及梯级模式下COD、η、Esp变化(B)图。
具体实施方式
下面结合具体实施方式来对本发明进行更进一步详细的说明，以更好地体现本发明的优势。
实施例
一种高效处理有机废水的梯级电催化氧化工艺设计方法，包括以下步骤：
S1、以PbO2/Ti电极作为阳极，Ti作为阴极，粒径为5mm活性炭作为填充颗粒，组成多个梯级反应器串联的填充床电极反应器；
S2、将有机废水以动态进水方式通入填充床电极反应器，接通电源使阳极、阴极通电，通过各梯级反应器对有机废水进行有机污染物处理；
S3、有机废水有机污染物处理期间，对梯级电催化氧化工艺进行修正及调整；
S3-1、通过修正电极电流密度的定量函数模型，以及反映多项体系耦合传质机制的实时传质模型，对梯级电催化氧化工艺施加的操作电流密度、传质系数进行修正与调整；
1)定量函数模型具体为：利用CV曲线计算伏安电量，建立的q、qin和qout之间的关系拟合数值，用式(1)、(2)、(3)计算。以式(4)计算粒子电极相对阳极的活化面积拓展系数(λ)；构建等效电路，拟合EIS得到有机物氧化电阻(Rct)和析氧电阻(R0)，用式(5)计算交换电流密度i0并与极化曲线所得数值进行比对以校核EIS等效电路模型，而后据式(6)计算参与有机废水氧化的阳极和粒子电极有效电流占比系数γ、β，再应用填充床电极反应器在反应控制(零级反应)状态下所得的有机废水中有机物降解实时数据并根据式(7)进行γ、β的修正；
q＝qin+qout (1)
q(v)＝qout+k1v-1/2 (2)
q(v)-1＝q-1+k2v1/2 (3)
λ＝qout,P/qout,A (4)
i0＝RT/(NFRct) (5)
γ或(λ+β)＝R0/(R0+Rct) (6)
Figure BDA0003645722630000091
其中，q为体系的总伏安电量；qin、qout分别为体系的内层、外层伏安电量；q(v)为CV曲线的积分面积；v为扫速；k1、k2为常数；qout,A、qout,P分别为阳极、粒子电极的外层伏安电量；R为气体常数；T为绝对温度；c0、c(t)分别为填充床电极反应器进水、处理出水有机物浓度；ε为床层空隙率；x0为极板间距；N为电子转移数量；F为法拉第常数；α为电流密度操作系数；km为传质系数；t为反应时间；c(t)为t时的填充床电极反应器处理出水有机物浓度。
2)实时传质模型具体为：结合雷诺数Re和施密特数Sc之间的关系并考虑梯级反应器几何构型、周向返混、液相特征等因素的影响，引入操作电流密度iR(以阳极几何面积计算)描述反应物实时传质进程：
Figure BDA0003645722630000092
其中，D为扩散系数；c0为反应物起始浓度；dp为填料平均粒径；RH为梯级反应器水力半径；y0为主反应区高度；χ为梯级反应器形状因子校正系数；ε为床层空隙率；N为电子转移数量；F为法拉第常数；
同时，利用在扩散控制(一级反应)状态下所得的有机废水有机物降解实时数据进行传质系数的修正，关联式如下所示：
Figure BDA0003645722630000101
其中，c(t)为t时的废水污染物浓度；c0为原水污染物浓度；λ为床层拓展系数；km为传质系数；ε为床层空隙率；x0为极板间距；t为反应时间；
利用式(9)所得的传质系数数值代入式(8)，校正常数项可获得反映多项体系耦合传质机制的实时传质模型；在填充床电极反应器工作时，其内部传质路径、有效电流数值、电极活化面积与水力学及电化学条件、废水物理属性、电极类型、反应器构型的耦合机制密切关联，它们之间的关系很难利用准确、完整的数学公式进行表达，而这也关系到电催化氧化进程控制步骤的判定，而通过本发明定量函数模型以及实时传质模型的建立，既能反映多维场复杂耦合实际情况又能精确表述有机污染物实时传质系数和操作电流密度数值的定量函数关系；
S3-2、通过有机物降解效能强化模型调整填充床电极反应器的梯级反应器数量，并对填充床电极反应器的影响参数进行效能预测与优化，影响参数包括：阳极材质、填料类型、极板间距、床层厚度与高度以及粒径与床层高度的比例、进水流速以及填料粒径等；
有机物降解效能强化模型具体为：填充床电极反应器在反应控制状态时经济因子之间的关系采用式(10)表示：
Figure BDA0003645722630000102
其中，p1、p2、p3表示经济指标因子，在电极和填充床电极反应器主体材料已定时为常数；Q为废水流量；b为直流电源经济因子；I为电流；Esp为电耗；gE为电价；
梯级操作流程的效能强化模型f(αn)为各梯级单元fn(αn,Xn)的总和；X已定时，对于基建费用(式10前两项)和处理成本(式10最后一项)有相反的影响，因此f(αn)存在极小值；由于∑Xn＝X(n＝1,2,3,......)，且Xn与操作电流密度αn、fn(αn,Xn)密切相关；f(αn)对第n级梯级反应器进水或出水口处(命名为节点，其数量为n+1)的有机物无量纲浓度Cj(Cj＝cj/c0)求导并代入有关的关联式可得以下方程：
Figure BDA0003645722630000111
其中，
Figure BDA0003645722630000112
为CFn第n级梯级反应器进水有机物无量纲浓度；γ为阳极对有机物的氧化选择系数；
Figure BDA0003645722630000116
为节点虚去除率，表示处理流程的年费用折算值随节点归一化浓度增加而增加的比率；
Figure BDA0003645722630000117
为第n级梯级反应器的年费用折算值随节点归一化浓度增加而增加的比率。
式(11)所示的关系为标准的遗传算法目标函数表达式，可采用以下步骤求解各梯级反应器的有机污染物去除率
Figure BDA0003645722630000113
在假定n值的前提下，拟定各节点有机物无量纲浓度初值Cj (0)，其中C1＝1、Cn+1＝1-X，其与节点无量纲浓度均在1和(1-X)之间取值；在该初值点处采用Taylor公式将
Figure BDA0003645722630000114
展开并舍却高阶项可得到遗传算子：
Figure BDA0003645722630000115
其中，θ＝γ2(p1x0n+p2)/[(γ+β)km,n]；x0n为第n级梯级反应器的极板间距；km,n为第n级梯级反应器的传质系数；ΔCj为各节点的归一化浓度初值微量；
Figure BDA0003645722630000121
为节点虚去除率，表示处理流程的初始年费用折算值随节点归一化浓度增加而增加的比率；
Figure BDA0003645722630000122
为第n级梯级反应器的初始年费用折算值随节点归一化浓度增加而增加的比率；CFn (0)为第n级梯级反应器进水有机物初始无量纲浓度；对于各节点Cj (0)反复施加应用式(12)计算的微量ΔCj，则梯级操作流程的f(αn)数值将不断减小，各节点满足式(11)的条件时取得极小值，使得f(αn)为最小值时的n值即为最佳梯级数量，进而可据Xn确定各梯级反应器几何尺寸与运行条件；采用不同的电极材料，在极板间距数值允许范围之内分别计算f(αn)的数值，f(αn)的数值越小则所选用的电极材料和极板间距更优；反应控制状态时，传质系数和操作电流密度在设定有机废水有机物去除率情况下可用于各梯级反应器的电耗、极板耗量、反应器大小等经济因子的确定，在由若干梯级反应器串联组成的梯级电催化氧化工艺中，各梯级反应器去除率相互约束，但其总值与所需有机物去除率相等；各梯级反应器有机物去除率组合情况不同时，梯级电催化氧化工艺的经济因子数值也随之发生变化，通过本发明构建的效能强化模型可对梯级电催化氧化工艺流程优化，从而指导梯级电催化氧化工艺流程的优化构建。
应用例1
以PbO2/Ti为阳极，Ti为阴极，粒径为5mm的活性炭为填充颗粒，以0.21mol/L的Na2SO4溶液为电解液。在室温条件下，采用动态进水方式，控制进水流量Q为0.50L/h，设置初始操作电流密度i为100A/m2，分别配置浓度为0.6、1.0、1.5g/L的苯酚模拟废水，以对应苯酚初始浓度下恒电流填充床体系出水口处的COD值为基础，考察不同苯酚初始浓度下，梯级体系以及恒电流体系对有机质的电催化氧化能力，实验结果如图2所示。由图2可知，各梯级体系中COD的下降速度均快于同等苯酚初始浓度下恒电流填充床体系中的COD下降速度。梯级体系在输入总电量较小的情况下，对COD的处理效果强于恒电流填充床体系，体现出梯级体系的高效性。
为更直观得对比梯级与恒电流填充床体系对苯酚的降解效果，引入氧化时空产率(Yst)，该值表示单位体积施加单位电流在单位时间内去除的污染物量，计算公式如下：
Figure BDA0003645722630000131
从图3总体来看，在不同苯酚初始浓度下，梯级体系的Yst均高于恒电流填充床体系。正是梯级体系的反应控制状态保证了体系在输入电流降低的情况下，仍有较强的有机物去除效果，高电流利用效率保证了梯级体系对苯酚的高处理效果以及低能耗特性。
由图4可知，相同苯酚初始浓度条件下，梯级体系的能耗明显低于恒电流填充床体系，表明相较于恒电流填充床体系，梯级体系在节能方面的性能较高，体现出梯级体系成为解决填充床体系高能耗问题的一种行之有效的方法。在梯级系统中，每一次操作电流密度的下降，会造成体系电位的显著降低，而经过前期的计算，已保证体系恒定处于反应控制阶段，即电流效率恒定为100％，高效的输入电能使用效率在保证有机质去除效果的同时，有效解决输入电能用于析氧副反应、转换为热能等形式的浪费问题，从而使得系统的电耗明显低于恒电流填充床体系，达成了对填充床反应器低能耗的改进目标。
不同苯酚初始浓度条件下，梯级体系与恒电流填充床体系的电流效率(ACE)如附图5所示。梯级体系在不同苯酚初始浓度下，在反应的全过程中均可以保持100％的电流效率，而在恒电流填充床体系中，在各苯酚初始浓度下，均呈现出ACE随反应时间的延长而逐渐降低的现象。
应用例2
如图1B所示反应机理，以IrO2–Ta2O5/Ti为阳极，Ti为阴极，粒径为5mm的活性炭为填充颗粒，实验装置如图1A所示，在连续流动模式下，将噻吩-2，5-二羧酸废水直接泵入装置，并在室温下通过系统进行单程处理。图6为本应用例梯级填充床电催化氧化工艺节能模式示意图。从图6中可以看出，电流密度越小，功耗越低。根据噻吩-2，5-二羧酸废水处理的模拟结果，图7A表明处理废水理论最低成本为9.86$/m3(n＝8)，其中基建投资和电力成本分别为3.80$/m3和6.06$/m3。图7B提供了连续流动电流控制氧化过程中的COD演变，从图7B中可以看出，梯级填充床电催化氧化工艺可以实现高效氧化，其电流效率(η)约为100％，功耗更低。

Claims (10)

1.一种高效处理有机废水的梯级电催化氧化工艺设计方法，其特征在于，包括以下步骤：
S1、以DSA电极作为阳极，Ti作为阴极，填充颗粒作为粒子电极，组成多个梯级反应器串联的填充床电极反应器；
S2、将有机废水以动态进水方式通入填充床电极反应器，接通电源使阳极、阴极通电，通过各梯级反应器对有机废水进行有机污染物处理；
S3、有机废水有机污染物处理期间，对梯级电催化氧化工艺进行修正及调整；
S3-1、通过修正电极电流密度的定量函数模型，以及反映多项体系耦合传质机制的实时传质模型，对梯级电催化氧化工艺施加的操作电流密度、传质系数进行修正与调整；
S3-2、通过有机物降解效能强化模型调整填充床电极反应器的梯级反应器数量，并对填充床电极反应器的影响参数进行效能预测与优化。
2.根据权利要求1所述的一种高效处理有机废水的梯级电催化氧化工艺设计方法，其特征在于，所述DSA电极采用刷涂热氧化法或电沉积法制备，如：PbO2/Ti电极。
3.根据权利要求1所述的一种高效处理有机废水的梯级电催化氧化工艺设计方法，其特征在于，所述填充颗粒为活性炭、涂膜活性炭中的任意一种或两种任意比的混合物。
4.根据权利要求3所述的一种高效处理有机废水的梯级电催化氧化工艺设计方法，其特征在于，所述活性炭、涂膜活性炭的粒径为5mm，涂膜活性炭采用浸泡、低温烘干定型方式制取。
5.根据权利要求1所述的一种高效处理有机废水的梯级电催化氧化工艺设计方法，其特征在于，步骤S3-1中用于修正电极电流密度的所述定量函数模型，具体为：
利用CV曲线计算伏安电量，建立的q、qin和qout之间的关系拟合数值，用式(1)、(2)、(3)计算；以式(4)计算粒子电极相对阳极的活化面积拓展系数(λ)；构建等效电路，拟合EIS得到有机物氧化电阻(Rct)和析氧电阻(R0)，用式(5)计算交换电流密度i0并与极化曲线所得数值进行比对以校核EIS等效电路模型，而后据式(6)计算参与有机废水氧化的阳极和粒子电极有效电流占比系数γ、β，再应用填充床电极反应器在反应控制状态下所得的有机废水中有机物降解实时数据并根据式(7)进行γ、β的修正；
q＝qin+qout (1)
q(v)＝qout+k1v-1/2 (2)
q(v)-1＝q-1+k2v1/2 (3)
λ＝qout,P/qout,A (4)
i0＝RT/(NFRct) (5)
γ或(λ+β)＝R0/(R0+Rct) (6)
Figure FDA0003645722620000021
其中，q为体系的总伏安电量；qin、qout分别为体系的内层、外层伏安电量；q(v)为CV曲线的积分面积；v为扫速；k1、k2为常数；qout,A、qout,P分别为阳极、粒子电极的外层伏安电量；R为气体常数；T为绝对温度；c0、c(t)分别为填充床电极反应器进水、处理出水有机物浓度；ε为床层空隙率；x0为极板间距；N为电子转移数量；F为法拉第常数；α为电流密度操作系数；km为传质系数；t为反应时间；c(t)为t时的填充床电极反应器处理出水有机物浓度。
6.根据权利要求1所述的一种高效处理有机废水的梯级电催化氧化工艺设计方法，其特征在于，步骤S3-1中用于反映多项体系耦合传质机制的所述实时传质模型，具体为：
结合雷诺数Re和施密特数Sc之间的关系并考虑梯级反应器几何构型、周向返混、液相特征等因素的影响，以阳极几何面积计算，引入操作电流密度iR描述反应物实时传质进程：
Figure FDA0003645722620000031
其中，D为扩散系数；c0为反应物起始浓度；dp为填料平均粒径；RH为梯级反应器水力半径；y0为主反应区高度；χ为梯级反应器形状因子校正系数；ε为床层空隙率；N为电子转移数量；F为法拉第常数；
同时，利用在扩散控制状态下所得的有机废水有机物降解实时数据进行传质系数的修正，关联式如下所示：
Figure FDA0003645722620000032
其中，c(t)为t时的废水污染物浓度；c0为原水污染物浓度；λ为床层拓展系数；km为传质系数；ε为床层空隙率；x0为极板间距；t为反应时间；
利用式(9)所得的传质系数数值代入式(8)，校正常数项可获得反映多项体系耦合传质机制的实时传质模型。
7.根据权利要求1所述的一种高效处理有机废水的梯级电催化氧化工艺设计方法，其特征在于，步骤S3-2中所述有机物降解效能强化模型具体为：
填充床电极反应器在反应控制状态时经济因子之间的关系采用式(10)表示：
Figure FDA0003645722620000041
其中，p1、p2、p3表示经济指标因子，在电极和填充床电极反应器主体材料已定时为常数；Q为废水流量；b为直流电源经济因子；I为电流；Esp为电耗；gE为电价；
梯级操作流程的效能强化模型f(αn)为各梯级单元fn(αn,Xn)的总和；X已定时，对于式(10)前两项的基建费用和式(10)最后一项的处理成本有相反的影响，因此f(αn)存在极小值；由于∑Xn＝X(n＝1,2,3,......)，且Xn与操作电流密度αn、fn(αn,Xn)密切相关；对第n级梯级反应器进水或出水口处命名为节点，其数量为n+1，f(αn)对第n级梯级反应器进水或出水口处的有机物无量纲浓度Cj(Cj＝cj/c0)求导并代入有关的关联式可得以下方程：
Figure FDA0003645722620000042
其中，
Figure FDA0003645722620000043
CFn为第n级梯级反应器进水有机物无量纲浓度；γ为阳极对有机物的氧化选择系数；
Figure FDA0003645722620000044
为节点虚去除率，表示处理流程的年费用折算值随节点归一化浓度增加而增加的比率；
Figure FDA0003645722620000045
为第n级梯级反应器的年费用折算值随节点归一化浓度增加而增加的比率。
式(11)所示的关系为标准的遗传算法目标函数表达式，可采用以下步骤求解各梯级反应器的有机污染物去除率
Figure FDA0003645722620000046
在假定n值的前提下，拟定各节点有机物无量纲浓度初值
Figure FDA0003645722620000047
其中C1＝1、Cn+1＝1-X，其与节点无量纲浓度均在1和(1-X)之间取值；在该初值点处采用Taylor公式将
Figure FDA0003645722620000048
展开并舍却高阶项可得到遗传算子：
Figure FDA0003645722620000051
其中，θ＝γ2(p1x0n+p2)/[(γ+β)km,n]；x0n为第n级梯级反应器的极板间距；km,n为第n级梯级反应器的传质系数；ΔCj为各节点的归一化浓度初值微量；
Figure FDA0003645722620000052
为节点虚去除率，表示处理流程的初始年费用折算值随节点归一化浓度增加而增加的比率；
Figure FDA0003645722620000053
为第n级梯级反应器的初始年费用折算值随节点归一化浓度增加而增加的比率；CFn (0)为第n级梯级反应器进水有机物初始无量纲浓度；对于各节点Cj (0)反复施加应用式(12)计算的微量ΔCj，则梯级操作流程的f(αn)数值将不断减小，各节点满足式(11)的条件时取得极小值，使得f(αn)为最小值时的n值即为最佳梯级数量，进而可据Xn确定各梯级反应器几何尺寸与运行条件。
8.根据权利要求7所述的一种高效处理有机废水的梯级电催化氧化工艺设计方法，其特征在于，采用不同的电极材料，在极板间距数值允许范围之内分别计算f(αn)的数值，f(αn)的数值越小则所选用的电极材料和极板间距更优。
9.根据权利要求1所述的一种高效处理有机废水的梯级电催化氧化工艺设计方法，其特征在于，所述DSA电极为PbO2/Ti电极。
10.根据权利要求1所述的一种高效处理有机废水的梯级电催化氧化工艺设计方法，其特征在于，所述步骤S3-2中影响参数包括：阳极材质、填料类型、极板间距、床层厚度与高度以及粒径与床层高度的比例、进水流速以及填料粒径。
CN202210529802.6A 2022-05-16 2022-05-16 一种高效处理有机废水的梯级电催化氧化工艺设计方法 Pending CN114906906A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210529802.6A CN114906906A (zh) 2022-05-16 2022-05-16 一种高效处理有机废水的梯级电催化氧化工艺设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210529802.6A CN114906906A (zh) 2022-05-16 2022-05-16 一种高效处理有机废水的梯级电催化氧化工艺设计方法

Publications (1)

Publication Number Publication Date
CN114906906A true CN114906906A (zh) 2022-08-16

Family

ID=82766588

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210529802.6A Pending CN114906906A (zh) 2022-05-16 2022-05-16 一种高效处理有机废水的梯级电催化氧化工艺设计方法

Country Status (1)

Country Link
CN (1) CN114906906A (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4012319A (en) * 1974-02-01 1977-03-15 Swift And Company Waste water treatment
CN101833314A (zh) * 2010-03-30 2010-09-15 深圳达实智能股份有限公司 污水处理控制系统及污水处理控制方法
CN106673139A (zh) * 2016-12-15 2017-05-17 中国矿业大学 电化学分段脱氮装置
CN106986426A (zh) * 2017-05-09 2017-07-28 中国矿业大学 一种用于废水处理的床电极反应器设计方法
CN113087092A (zh) * 2021-04-12 2021-07-09 浙江树人学院（浙江树人大学） 一种基于电流密度梯级控制的废水电氧化处理方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4012319A (en) * 1974-02-01 1977-03-15 Swift And Company Waste water treatment
CN101833314A (zh) * 2010-03-30 2010-09-15 深圳达实智能股份有限公司 污水处理控制系统及污水处理控制方法
CN106673139A (zh) * 2016-12-15 2017-05-17 中国矿业大学 电化学分段脱氮装置
CN106986426A (zh) * 2017-05-09 2017-07-28 中国矿业大学 一种用于废水处理的床电极反应器设计方法
CN113087092A (zh) * 2021-04-12 2021-07-09 浙江树人学院（浙江树人大学） 一种基于电流密度梯级控制的废水电氧化处理方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张一弛: ""电催化氧化苯酚节能策略研究"", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 *
张一弛等: "电催化氧化有机污染物的电极性能评估策略", 《中国环境科学》 *

Similar Documents

Publication Publication Date Title
CN108193225A (zh) 一种膜电极构型co2电还原电解池
CN108439553B (zh) 一种用于去除水中氯离子的电化学氧化装置
CN110194508A (zh) 一种连续式电催化反应系统及其电流间接调节方法
CN114069661B (zh) 一种适于波动性输入的电解水制氢系统性能优化方法
CN106006866B (zh) 基于变电位控制的废水电氧化处理方法
CN113800606B (zh) 一种循环冷却水处理用涂层钛阳极、制备方法及应用
CN102583662B (zh) 一种网板柱塞流电解装置及用于处理有机废水的方法
Mirzazadeh et al. Optimization of caustic current efficiency in a zero-gap advanced chlor-alkali cell with application of genetic algorithm assisted by artificial neural networks
Ying et al. Performance of electrochemical cell with various flow channels for Bunsen reaction in the sulfur–iodine hydrogen production process
Fonseca et al. Improving microbial electrolysis stability using flow-through brush electrodes and monitoring anode potentials relative to thermodynamic minima
CN114906906A (zh) 一种高效处理有机废水的梯级电催化氧化工艺设计方法
Shojaei et al. 1-D transient microbial fuel cell simulation considering biofilm growth and temperature variation
Yao et al. Modeling and characterization of the biochar electrodes for vanadium redox flow battery
WO2020133611A1 (zh) 一种基于操作模式动态匹配的重金属废水净化控制方法
Adams et al. Electrochemical oxidation of ferrous iron in very dilute solutions
CN114606523B (zh) 一种变功率制氢系统及其自动控制系统和方法
Xu et al. Enhanced electricity generation and H2O2 production in a photocatalytic fuel cell and Fenton hybrid system assisted with reverse electrodialysis
CN106521610B (zh) 一种六价铬电镀用组合钛阳极及其制备方法
CN210122499U (zh) 一种连续式电催化反应系统
Deng et al. Spatiotemporal distribution model for zinc electrowinning process and its parameter estimation
Azwara et al. Optimal production of biohydrogen gas via Microbial Electrolysis cells (MEC) in a controlled batch reactor system
CN114781188A (zh) 一种面向波动电能消纳的固体氧化物电解制氢系统动态建模方法
Grau et al. Electrochemical removal of cadmium from dilute aqueous solutions using a rotating cylinder electrode of wedge wire screens
CN107576913B (zh) 主动配电网中氢负荷运行分析方法
CN113314723A (zh) 一种碳基催化剂的电化学性能及其在微生物燃料电池中的应用

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination