[发明专利]一种连续流交替式污水处理工艺的两相循环模拟方法

权利要求书

1.一种连续流交替式污水处理工艺的两相循环模拟方法，所述连续流交替式污水处理工艺由中间池和两侧的边池串联组成，其特征在于，将边池在一个连续流交替式污水处理工艺运行周期内的状态分成两个阶段：完全混合反应阶段和分层沉淀阶段，所述完全混合反应阶段将边池视为处于完全混合的均相状态，仅发生生化反应；所述分层沉淀阶段将边池沿水深方向分成若干层进行污泥沉淀，同时发生重力沉淀和生化反应；中间池始终处于完全混合反应阶段；

对完全混合反应阶段，采用生化反应动力学模型进行建模，得到各物质浓度的净变化率的表达式；

对分层沉淀阶段，将边池沿水深方向分成10层，采用沉淀模型和生化反应动力学模型进行建模，得到每一层上各物质浓度的净变化率的表达式；

模拟所涉及到污水中的物质有：溶解氧、溶解性有机物、颗粒性有机物、氨氮、硝态氮、异养菌、自养菌；

1)当边池处于完全混合反应阶段时及中间池中，各物质浓度的净变化率的表达式为：

溶解氧：

溶解性有机物：

颗粒性有机物：

氨氮：

硝态氮：

异养菌：

自养菌：

式中：

S_O：溶解氧的浓度，单位mgO₂/L；

S_COD：溶解性有机物的浓度，单位mgCOD/L；

S_NH：氨态氮的浓度，单位mgN/L；

S_NO：硝态氮的浓度，单位mgN/L；

X_COD：颗粒性有机物的浓度，单位mgCOD/L；

X_H：异养菌的浓度，单位mgCOD/L；

X_A：自养菌的浓度，单位mgCOD/L；

Q_in：边池或中间池的进水流量，单位m³/d；

Q_out：边池或中间池的出水流量，单位m³/d；

V：边池或中间池的总容积，单位m³；

S_o,in：边池或中间池进水中溶解氧的浓度，单位mgO₂/L；

S_COD,in：边池或中间池进水中溶解性有机物的浓度，单位mgCOD/L；

S_NH,in：边池或中间池进水中氨态氮的浓度，单位mgN/L；

S_NO,in：边池或中间池进水中硝态氮的浓度，单位mgN/L；

X_COD,in：边池或中间池进水中颗粒性有机物的浓度，单位mgCOD/L；

X_H,in：边池或中间池进水中异养菌的浓度，单位mgCOD/L；

X_A,in：边池或中间池进水中自养菌的浓度，单位mgCOD/L；

μ_H：异养菌最大比生长速率，单位h^-1；

μ_A：自养菌最大比生长速率，单位h^-1；

K_S：异养菌生长半速率常数，单位mgCOD/L；

b_H：异养菌比衰减速率，单位h^-1；

K_OH：异养菌氧呼吸半速率常数，单位mgO2/L；

K_NO：异养菌硝态氮呼吸半速率常数，单位mgN/L；

η_g：异养菌缺氧状态生长修正系数；

k_h：水解速率，单位mgCOD/(mg生物量(COD)·h)；

K_X：水解半速率常数，单位mgCOD/(mg生物量(COD)；

η_h：缺氧状态水解修正系数；

K_OA：自养菌氧呼吸半速率常数，单位mgO₂/L；

K_NH：自养菌生长半速率常数，单位mgN/L；

b_A：自养菌比衰减速率，单位h^-1；

Y_H：异养菌产率系数，单位mg生物量(COD)/mg COD去除；

Y_A：自养菌产率系数，单位mg生物量(COD)/mg被氧化的N；

对上述各项物质浓度随时间的净变化率进行联立得到微分方程组，进行数值求解，得中间池和处于完全混合反应阶段的边池中，各物质浓度随时间的变化；

2)当边池处于分层沉淀阶段时，每一层上各固体物质浓度的净变化率的表达式为：

颗粒性有机物：

异养菌：

自养菌：

式中：

Q_out是边池的出水流量，单位m³/d；

S是边池的表面积，单位m²；

h_m是当前层的高度，数值上等于池深和层数的比值，单位m；

X_COD,m，X_H,m和X_A,m分别是当前层上的颗粒性有机物、异养菌和自养菌浓度，单位g/m³；

X_COD,m+1，X_H,m+1和X_A,m+1分别是当前层上一层的颗粒性有机物、异养菌和自养菌浓度，单位g/m³；

J_XCOD,m-1和J_XCOD,m分别是当前层下面一层和当前层颗粒性有机物的通量；

J_XH,m-1和J_XH,m分别是当前层下面一层和当前层异养菌的通量；

J_XA,m-1和J_XA,m分别是当前层下面一层和当前层自养菌的通量。

当前层的通量算法：

J_XCOD,m＝V_mX_COD,m

J_XH,m＝V_mX_H,m

J_XA,m＝V_mX_A,m

V_m是当前层的固体沉降速度，通过下式计算：

式中，V₀是最大理论沉降速度，单位m/d；

r_h为沉降参数，单位为m³/g固体；

f_ns是边池进水中不可沉降物质所占的最大比例；

X_in是边池进水中的总悬浮固体浓度，g/m³；

X_in是边池进水中的总悬浮固体浓度，通过下式计算：

其中，X_i,in代表边池进水中的各种固体物质，有颗粒性有机物X_COD,in、异养菌X_H,in、自养菌X_A,,in，它们的加和作为进水中的总悬浮固体浓度；

X_m是当前层的总悬浮固体浓度，通过下式计算：

其中X_i,m代表当前层中的各种固体物质，有颗粒性有机物X_COD、异养菌X_H、自养菌X_A，它们浓度的加和作为当前层的总悬浮固体浓度；

μ_H：异养菌最大比生长速率，单位h^-1；

μ_A：自养菌最大比生长速率，单位h^-1；

K_S：异养菌生长半速率常数，单位mgCOD/L；

b_H：异养菌比衰减速率，单位h^-1；

K_OH：异养菌氧呼吸半速率常数，单位mgO2/L；

K_NO：异养菌硝态氮呼吸半速率常数，单位mgN/L；

η_g：异养菌缺氧状态生长修正系数；

k_h：水解速率，单位mgCOD/(mg生物量(COD)·h)；

K_X：水解半速率常数，单位mgCOD/(mg生物量(COD)；

η_h：缺氧状态水解修正系数；η_g

K_OA：自养菌氧呼吸半速率常数，单位mgO₂/L；

K_NH：自养菌生长半速率常数，单位mgN/L；

b_A：自养菌比衰减速率，单位h^-1；

对上述各项物质浓度随时间的净变化率进行联立得到微分方程组，进行数值求解，得处于分层沉淀阶段的边池中，每一层上各固体物质浓度随时间的变化；

3)当边池处于分层沉淀阶段时，每一层上各溶解性物质浓度的净变化率的表达式为：

溶解氧：

溶解性有机物：

氨氮：

硝态氮：

式中：

Q_out：边池的出水流量，单位m³/d；

V_layer：边池每层的容积，单位m³；

S_o,m、S_o,m+1：分别是当前层和当前层上面一层溶解氧的浓度，单位mgO₂/L；

S_COD,m、S_COD,m+1：分别是当前层和当前层上面一层溶解性有机物的浓度，单位mgCOD/L；

S_NH,m、S_NH,m+1：分别是当前层和当前层上面一层氨态氮的浓度，单位mgN/L；

S_NO,m、S_NO,m+1：分别是当前层和当前层上面一层硝态氮的浓度，单位mgN/L；

μ_H：异养菌最大比生长速率，单位h^-1；

μ_A：自养菌最大比生长速率，单位h^-1；

K_S：异养菌生长半速率常数，单位mgCOD/L；

K_OH：异养菌氧呼吸半速率常数，单位mgO2/L；

K_NO：异养菌硝态氮呼吸半速率常数，单位mgN/L；

K_OA：自养菌氧呼吸半速率常数，单位mgO₂/L；

K_NH：自养菌生长半速率常数，单位mgN/L；

Y_H：异养菌产率系数，单位mg生物量(COD)/mg COD去除；

Y_A：自养菌产率系数，单位mg生物量(COD)/mg被氧化的N；

X_H,m：当前层异养菌的浓度，单位mgCOD/L；

X_A,m：当前层自养菌的浓度，单位mgCOD/L；

对上述各项物质浓度随时间的净变化率进行联立得到微分方程组，进行数值求解，得处于分层沉淀阶段的边池中，每一层上各溶解性物质浓度随时间的变化。