一种A2O生物池工艺建模方法

摘要

本发明涉及污水处理领域，提供了一种A2O生物池工艺建模方法，该方法包括：设定输入参数及输出参数；根据所述输入参数、输出参数、ASM2D模型以及物料平衡方程建立A2O生物池工艺模型，所述A2O生物池工艺模型包括厌氧池工艺模型、缺氧池工艺模型以及好氧池工艺模型，所述输出参数为组分浓度变化率；该A2O生物池工艺模型简单，且通过Matlab仿真软件即可实现，用以模拟仿真A2O生物池工艺，从而为后期的控制策略研究提供了基础，同时避免了购买国外的类似于ASMs系列模拟软件高昂费用，节约了A2O生物池工艺处理污水的成本。

主权项

一种A2O生物池工艺建模方法，该方法基于仿真软件实现，其特征在于，所述方法包括以下步骤：设定输入参数及输出参数；根据所述输入参数、输出参数、ASM2D模型以及物料平衡方程建立A2O生物池工艺模型，所述A2O生物池工艺模型包括厌氧池工艺模型、缺氧池工艺模型以及好氧池工艺模型；其中，所述物料平衡方程为：$\frac{{\mathrm{\Delta Z}}_i}{\Delta t}=\frac{1}{V}(Q_{in}{Z}_{in}+r_{i}V-Q_{out}{Z}_{out})---\left(1\right)$其中，△Z_i为组分浓度变化量，△t为时间变化量，V为该池的体积，Q_in为进入水量，Z_in为进入组分浓度，r_i为影响组分浓度变化的过程速率，Q_out为流出水量，Z_out为流出组分浓度；所述厌氧池工艺模型包括以下微分方程：$\begin{array}{c}\frac{{\mathrm{dS}}_{O2,1}}{dt}=\frac{Q_e}{V_y}{S}_{O2,0}+\frac{Q_r}{V_y}{S}_{O2,4}-\frac{Q_1}{V_y}{S}_{O2,1}-0.6{\rho }_4-0.6{\rho }_5-0.2{\rho }_{11}\\ -0.6{\rho }_{13}-18{\rho }_{18}\end{array}---\left(2\right)$$\begin{array}{c}\frac{{\mathrm{dS}}_{A,1}}{dt}=\frac{Q_e}{V_y}{S}_{A,0}+\frac{Q_r}{V_y}{S}_{A,4}-\frac{Q_1}{V_y}{S}_{A,1}-1.6{\rho }_5-1.6{\rho }_7+{\rho }_8\\ -{\rho }_{10}+{\rho }_{17}\end{array}---\left(3\right)$$\begin{array}{c}\frac{{\mathrm{dS}}_{F,1}}{dt}=\frac{Q_e}{V_y}{S}_{F,0}+\frac{Q_r}{V_y}{S}_{F,4}-\frac{Q_1}{V_y}{S}_{F,1}+{\rho }_1+{\rho }_2+{\rho }_3-1.6{\rho }_4\\ -1.6{\rho }_6-{\rho }_8\end{array}---\left(4\right)$$\begin{array}{c}\frac{{\mathrm{dS}}_{{\mathrm{NH}}_4,1}}{dt}=\frac{Q_e}{V_y}{S}_{{\mathrm{NH}}_4,0}+\frac{Q_r}{V_y}{S}_{{\mathrm{NH}}_4,4}-\frac{Q_1}{V_y}{S}_{{\mathrm{NH}}_4,1}+0.01{\rho }_1+0.01{\rho }_2\\ +0.01{\rho }_3-0.022{\rho }_4-0.07{\rho }_5-0.022{\rho }_6-0.07{\rho }_7\\ +0.03{\rho }_8+0.03{\rho }_9-0.07{\rho }_{13}-0.07{\rho }_{14}+0.031{\rho }_{15}\\ -4.24{\rho }_{18}+0.031{\rho }_{19}\end{array}---\left(5\right)$$\begin{array}{c}\frac{{\mathrm{dS}}_{{\mathrm{NO}}_3,1}}{dt}=\frac{Q_e}{V_y}{S}_{{\mathrm{NO}}_3,0}+\frac{Q_r}{V_y}{S}_{{\mathrm{NO}}_3,4}-\frac{Q_1}{V_y}{S}_{{\mathrm{NO}}_3,1}-0.21{\rho }_6-0.21{\rho }_7\\ -0.07{\rho }_{12}-0.21{\rho }_{14}+4.17{\rho }_{18}\end{array}---\left(6\right)$$\begin{array}{c}\frac{{\mathrm{dS}}_{{\mathrm{PO}}_4,1}}{dt}=\frac{Q_e}{V_y}{S}_{{\mathrm{PO}}_4,0}+\frac{Q_r}{V_y}{S}_{{\mathrm{PO}}_4,4}-\frac{Q_1}{V_y}{S}_{{\mathrm{PO}}_4,1}-0.004{\rho }_4-0.02{\rho }_5\\ -0.004{\rho }_6-0.02{\rho }_7+0.01{\rho }_8+0.01{\rho }_9+0.04{\rho }_{10}\\ -{\rho }_{11}-{\rho }_{12}-0.02{\rho }_{13}+0.01{\rho }_{15}+{\rho }_{16}-0.02{\rho }_{19}+0.01{\rho }_{19}\end{array}---\left(7\right)$$\frac{{\mathrm{dX}}_{I,1}}{dt}=\frac{Q_e}{V_y}{X}_{I,0}+\frac{Q_r}{V_y}{X}_{I,4}-\frac{Q_1}{V_y}{X}_{I,1}+0.1{\rho }_9+0.1{\rho }_{15}+0.1{\rho }_{19}---\left(8\right)$$\begin{array}{c}\frac{{\mathrm{dX}}_{S,1}}{dt}=\frac{Q_e}{V_y}{X}_{S,0}+\frac{Q_r}{V_y}{X}_{S,4}-\frac{Q_1}{V_y}{X}_{S,1}-{\rho }_1-{\rho }_2-{\rho }_3+0.9{\rho }_9\\ +0.9{\rho }_{15}+0.9{\rho }_{19}\end{array}---\left(9\right)$$\frac{{\mathrm{dX}}_{H,1}}{dt}=\frac{Q_e}{V_y}{X}_{H,0}+\frac{Q_r}{V_y}{X}_{H,4}-\frac{Q_1}{V_y}{X}_{H,1}+{\rho }_4+{\rho }_5+{\rho }_6+{\rho }_7-{\rho }_9---\left(10\right)$$\frac{{\mathrm{dX}}_{PAO,1}}{dt}=\frac{Q_e}{V_y}{X}_{PAO,0}+\frac{Q_r}{V_y}{X}_{PAO,4}-\frac{Q_1}{V_y}{X}_{PAO,1}+{\rho }_{13}+{\rho }_{14}-{\rho }_{15}---\left(11\right)$$\frac{{\mathrm{dX}}_{PP,1}}{dt}=\frac{Q_e}{V_y}{X}_{PP,0}+\frac{Q_r}{V_y}{X}_{PP,4}-\frac{Q_1}{V_y}{X}_{PP,1}-0.4{\rho }_{10}+{\rho }_{11}+{\rho }_{12}-{\rho }_{16}---\left(12\right)$$\frac{{\mathrm{dX}}_{PHA,1}}{dt}=\frac{Q_e}{V_y}{X}_{PHA,0}+\frac{Q_r}{V_y}{X}_{PHA,4}-\frac{Q_1}{V_y}{X}_{PHA,1}+{\rho }_{10}-{\rho }_{17}---\left(13\right)$$\frac{{\mathrm{dX}}_{AUT,1}}{dt}=\frac{Q_e}{V_y}{X}_{AUT,0}+\frac{Q_r}{V_y}{X}_{AUT,4}-\frac{Q_1}{V_y}{X}_{AUT,1}-{\rho }_{19}---\left(14\right)$所述缺氧池工艺模型包括以下微分方程：$\begin{array}{c}\frac{{\mathrm{dS}}_{O_2,2}}{dt}=\frac{Q_1}{V_q}{S}_{O2,1}+\frac{Q_t}{V_q}{S}_{O2,3}-\frac{Q_2}{V_q}{S}_{O2,2}-0.6{\rho }_4-0.6{\rho }_5-0.2{\rho }_{11}\\ -0.6{\rho }_{13}-18{\rho }_{18}\end{array}---\left(15\right)$$\begin{array}{c}\frac{{\mathrm{dS}}_{A,2}}{dt}=\frac{Q_1}{V_q}{S}_{A,1}+\frac{Q_t}{V_q}{S}_{A,3}-\frac{Q_2}{V_q}{S}_{A,2}-1.6{\rho }_5-1.6{\rho }_7+{\rho }_8\\ -{\rho }_{10}+{\rho }_{17}\end{array}---\left(16\right)$$\begin{array}{c}\frac{{\mathrm{dS}}_{F,2}}{dt}=\frac{Q_1}{V_q}{S}_{F,1}+\frac{Q_t}{V_q}{S}_{F,3}-\frac{Q_2}{V_q}{S}_{F,2}+{\rho }_1+{\rho }_2+{\rho }_3-1.6{\rho }_4\\ -1.6{\rho }_6-{\rho }_8\end{array}---\left(17\right)$$\begin{array}{c}\frac{{\mathrm{dS}}_{{\mathrm{NH}}_4,2}}{dt}=\frac{Q_1}{V_q}{S}_{{\mathrm{NH}}_4,1}+\frac{Q_t}{V_q}{S}_{{\mathrm{NH}}_4,3}-\frac{Q_2}{V_q}{S}_{{\mathrm{NH}}_4,2}+0.01{\rho }_1+0.01{\rho }_2\\ +0.01{\rho }_3-0.022{\rho }_4-0.07{\rho }_5-0.022{\rho }_6-0.07{\rho }_7\\ +0.03{\rho }_8+0.03{\rho }_9-0.07{\rho }_{13}-0.07{\rho }_{14}+0.031{\rho }_{15}\\ -4.24{\rho }_{18}+0.031{\rho }_{19}\end{array}---\left(18\right)$$\begin{array}{c}\frac{{\mathrm{dS}}_{{\mathrm{NO}}_3,2}}{dt}=\frac{Q_1}{V_q}{S}_{{\mathrm{NO}}_3,1}+\frac{Q_t}{V_q}{S}_{{\mathrm{NO}}_3,3}-\frac{Q_2}{V_q}{S}_{{\mathrm{NO}}_3,2}-0.21{\rho }_6-0.21{\rho }_7-0.07{\rho }_{12}\\ -0.21{\rho }_{14}+4.17{\rho }_{18}\end{array}---\left(19\right)$$\begin{array}{c}\frac{{\mathrm{dS}}_{{\mathrm{PO}}_4,2}}{dt}=\frac{Q_1}{V_q}{S}_{{\mathrm{PO}}_4,1}+\frac{Q_t}{V_q}{S}_{{\mathrm{PO}}_4,3}-\frac{Q_2}{V_q}{S}_{{\mathrm{PO}}_4,2}-0.004{\rho }_4-0.02{\rho }_5-0.004{\rho }_6\\ -0.02{\rho }_7+0.01{\rho }_8+0.01{\rho }_9+0.4{\rho }_{10}-{\rho }_{11}-{\rho }_{12}-0.02{\rho }_{13}\\ +0.01{\rho }_{15}+{\rho }_{16}-0.02{\rho }_{18}+0.01{\rho }_{19}\end{array}---\left(20\right)$$\frac{{\mathrm{dX}}_{I,2}}{dt}=\frac{Q_1}{V_q}{X}_{I,1}+\frac{Q_t}{V_q}{X}_{I,3}-\frac{Q_2}{V_q}{X}_{I,2}+0.1{\rho }_9+0.1{\rho }_{15}+0.1{\rho }_{19}---\left(21\right)$$\begin{array}{c}\frac{{\mathrm{dX}}_{S,2}}{dt}=\frac{Q_1}{V_y}{X}_{S,1}+\frac{Q_t}{V_q}{X}_{S,3}-\frac{Q_2}{V_q}{X}_{S,2}-{\rho }_1-{\rho }_2-{\rho }_3+0.9{\rho }_9\\ +0.9{\rho }_{15}+0.9{\rho }_{19}\end{array}---\left(22\right)$$\frac{{\mathrm{dX}}_{H,2}}{dt}=\frac{Q_1}{V_q}{X}_{H,1}+\frac{Q_t}{V_q}{X}_{H,3}-\frac{Q_2}{V_q}{X}_{H,2}+{\rho }_4+{\rho }_5+{\rho }_6+{\rho }_7-{\rho }_9---\left(23\right)$$\frac{{\mathrm{dX}}_{PAO,2}}{dt}=\frac{Q_t}{V_q}{X}_{PAO,1}+\frac{Q_t}{V_q}{X}_{PAO,3}-\frac{Q_2}{V_q}{X}_{PAO,2}+{\rho }_{13}+{\rho }_{14}-{\rho }_{15}---\left(24\right)$$\frac{{\mathrm{dX}}_{PP,2}}{dt}=\frac{Q_t}{V_q}{X}_{PP,1}+\frac{Q_t}{V_q}{X}_{PP,3}-\frac{Q_2}{V_q}{X}_{PP,2}-0.4{\rho }_{10}+{\rho }_{11}+{\rho }_{12}-{\rho }_{16}---\left(25\right)$$\frac{{\mathrm{dX}}_{PHA,2}}{dt}=\frac{Q_t}{V_q}{X}_{PHA,1}+\frac{Q_t}{V_q}{X}_{PHA,3}-\frac{Q_2}{V_q}{X}_{PHA,2}+{\rho }_{10}-{\rho }_{17}---\left(26\right)$$\frac{{\mathrm{dX}}_{AUT,2}}{dt}=\frac{Q_1}{V_q}{X}_{AUT,1}+\frac{Q_t}{V_q}{X}_{AUT,3}-\frac{Q_2}{V_q}{X}_{AUT,2}-{\rho }_{19}---\left(27\right)$所述好氧池工艺模型包括以下微分方程：$\begin{array}{c}\frac{{\mathrm{dS}}_{O_2,3}}{dt}=\frac{Q_2}{V_h}{S}_{O2,2}-\frac{Q_t+Q_3+Q_r+Q_w}{V_h}{S}_{O2,3}-0.6{\rho }_4-0.6{\rho }_5\\ -0.2{\rho }_{11}-0.6{\rho }_{13}-18{\rho }_{18}+KLa\left(8.56-S_{O2,3}\right)\end{array}---\left(28\right)$$\begin{array}{c}\frac{{\mathrm{dS}}_{A,3}}{dt}=\frac{Q_2}{V_h}{S}_{A,2}-\frac{Q_t+Q_3+Q_r+Q_w}{V_h}{S}_{A,3}-1.6{\rho }_5-1.6{\rho }_7\\ +{\rho }_8-{\rho }_{10}+{\rho }_{17}\end{array}---\left(29\right)$$\begin{array}{c}\frac{{\mathrm{dS}}_{F,3}}{dt}=\frac{Q_2}{V_h}{S}_{F,2}-\frac{Q_t+Q_3+Q_r+Q_w}{V_h}{S}_{F,3}+{\rho }_1+{\rho }_2+{\rho }_3\\ -1.6{\rho }_4-1.6{\rho }_6-{\rho }_8\end{array}---\left(30\right)$$\begin{array}{c}\frac{{\mathrm{dS}}_{{\mathrm{NH}}_4,3}}{dt}=\frac{Q_2}{V_h}{S}_{{\mathrm{NH}}_4,2}-\frac{Q_t+Q_3+Q_r+Q_w}{V_h}{S}_{{\mathrm{NH}}_4,3}+0.01{\rho }_1+0.01{\rho }_2\\ +0.01{\rho }_3-0.022{\rho }_4-0.07\rho 5-0.022\rho 6-0.07{\rho }_7\\ +0.03{\rho }_8+0.03{\rho }_9-0.07{\rho }_{13}-0.07{\rho }_{14}+0.031{\rho }_{15}\\ -4.24{\rho }_{18}+0.031{\rho }_{19}\end{array}---\left(31\right)$$\begin{array}{c}\frac{{\mathrm{dS}}_{{\mathrm{NO}}_3,3}}{dt}=\frac{Q_2}{V_h}{S}_{{\mathrm{NO}}_3,2}-\frac{Q_t+Q_3+Q_r+Q_w}{V_h}{S}_{{\mathrm{NO}}_3,3}-0.21{\rho }_6-0.21{\rho }_7\\ -0.07{\rho }_{12}-0.21{\rho }_{14}+4.17{\rho }_{18}\end{array}---\left(32\right)$$\begin{array}{c}\frac{{\mathrm{dS}}_{{\mathrm{PO}}_4,3}}{dt}=\frac{Q_2}{V_h}{S}_{{\mathrm{PO}}_4,2}-\frac{Q_t+Q_3+Q_r+Q_w}{V_h}{S}_{{\mathrm{PO}}_4,3}-0.004{\rho }_4-0.02{\rho }_5\\ -0.004{\rho }_6-0.02{\rho }_7++0.01{\rho }_8+0.01{\rho }_9+0.4{\rho }_{10}-{\rho }_{11}\\ -{\rho }_{12}-0.02{\rho }_{13}+0.01{\rho }_{15}++{\rho }_{16}-0.02{\rho }_{18}+0.01{\rho }_{19}\end{array}---\left(33\right)$$\begin{array}{c}\frac{{\mathrm{dX}}_{I,3}}{dt}=\frac{Q_2}{V_h}{X}_{I,2}-\frac{Q_t+Q_3+Q_r+Q_w}{V_h}{X}_{I,3}+0.1{\rho }_9\\ +0.1{\rho }_{15}+0.1{\rho }_{19}\end{array}---\left(34\right)$$\begin{array}{c}\frac{{\mathrm{dX}}_{S,3}}{dt}=\frac{Q_2}{V_h}{X}_{S,2}-\frac{Q_t+Q_3+Q_r+Q_w}{V_h}{X}_{S,3}-{\rho }_1-{\rho }_2-{\rho }_3\\ +0.9{\rho }_9+0.9{\rho }_{15}+0.9{\rho }_{19}\end{array}---\left(35\right)$$\begin{array}{c}\frac{{\mathrm{dX}}_{H,3}}{dt}=\frac{Q_2}{V_h}{X}_{H,2}-\frac{Q_t+Q_3+Q_r+Q_w}{V_h}{X}_{H,3}+{\rho }_4+{\rho }_5\\ +{\rho }_6+{\rho }_7-{\rho }_9\end{array}---\left(36\right)$$\begin{array}{c}\frac{{\mathrm{dX}}_{PAO,3}}{dt}=\frac{Q_2}{V_h}{X}_{PAO,2}-\frac{Q_t+Q_3+Q_r+Q_w}{V_h}{X}_{PAO,3}+{\rho }_{13}\\ +{\rho }_{14}-{\rho }_{15}\end{array}---\left(37\right)$$\begin{array}{c}\frac{{\mathrm{dX}}_{PP,3}}{dt}=\frac{Q_2}{V_h}{X}_{PP,2}-\frac{Q_t+Q_3+Q_r+Q_w}{V_h}{X}_{PP,3}-0.4{\rho }_{10}+{\rho }_{11}\\ +{\rho }_{12}-{\rho }_{16}\end{array}---\left(38\right)$$\frac{{\mathrm{dX}}_{PHA,3}}{dt}=\frac{Q_2}{V_h}{X}_{PHA,2}-\frac{Q_t+Q_3+Q_r+Q_w}{V_h}{X}_{PHA,3}+{\rho }_{10}-{\rho }_{17}---\left(39\right)$$\frac{{\mathrm{dX}}_{AUT,3}}{dt}=\frac{Q_2}{V_h}{X}_{AUT,2}-\frac{Q_t+Q_3+Q_r+Q_w}{V_h}{X}_{AUT,3}-{\rho }_{19}---\left(40\right)$公式2‑40中，Q_e为生物池进水流量，Q_r为污泥回流量，Q₁为厌氧池出水量，Q₂为缺氧池出水量，Q₃为好氧池出水量，Q_t混合液回流量，Q_w为排泥量，其中Q₁＝Q_e+Q_r,Q₂＝Q₁+Q_t,Q₃＝Q₂‑Q_t‑Q_r‑Q_w,V_y为厌氧池的容积，V_q为缺氧池的容积，V_h为好氧池的容积，t为反应时间，KLa是氧传质系数，S_O2、S_F、S_A、S_NH4、S_NO3和S_PO4均为溶解性组分，S_O2、S_F、S_A、S_NH4、S_NO3和S_PO4分别表示溶解氧、可发酵的易生物降解有机物、发酵产物、氨氮、硝酸盐氮与亚硝酸盐氮、溶解态无机磷，ρ₁，ρ₂...，ρ₁₉为ASM2D模型中的过程速率，X_I、X_s、X_H、X_PAO、X_PP、X_PHA和X_AUT均为颗粒性组分，X_I、X_s、X_H、X_PAO、X_PP、X_PHA和X_AUT分别表示惰性颗粒性有机物、慢速可降解基质、异养微生物、聚磷菌PAO、聚磷酸盐、聚磷菌PAO的胞内贮存物、硝化菌，式2‑40中S_O2、S_F、S_A、S_NH4、S_NO3、S_PO4、S_O2、S_F、S_A、S_NH4、S_NO3和S_PO4的下标中逗号后的0表示生物池，下标中逗号后的1、2、3、4分别代表厌氧池、缺氧池、好氧池和二沉池。