发明名称 城市污水处理厂出水有机物浓度的随机过程预测方法
摘要 本发明公开了一种城市污水处理厂出水有机物浓度的随机过程预测方法,即以污水处理厂出水有机物为研究对象,建立进出水有机物的质量守恒微分方程,并引入随机过程,推导得到城市污水处理厂出水有机物浓度的预测结果。本发明需要实测异养菌的产率系数、最大比生长速率及其半饱和常数,同时利用污水处理厂曝气池混合液挥发性悬浮固体浓度、水力停留时间、实际进出水总BOD5及日常实测数据,即可预测当实际进水水质水量等波动时出水有机物浓度的动态变化过程。本发明为新型的活性污泥动态模型方法,具有操作简便可行、实测参数少、预测值与实测值拟合较好等优点。
申请公布号 CN101786721A 申请公布日期 2010.07.28
申请号 CN201019161004.7 申请日期 2010.02.05
申请人 重庆大学 发明人 许劲;涂茂;张海霞;谈涛;齐龙;高溢璟;洪国强;赵绪光
分类号 C02F3/12(2006.01)I;G06F17/15(2006.01)I 主分类号 C02F3/12(2006.01)I
代理机构 重庆博凯知识产权代理有限公司 50212 代理人 张先芸
主权项 1.城市污水处理厂出水有机物浓度的随机过程预测方法,其特征在于包括下述步骤:1)利用城市污水处理厂的现有设备实测异养菌产率系数Y<sub>H</sub>、异养菌最大比生长速率<img file="FSA00000024924700011.GIF" wi="57" he="54" />与半饱和常数K<sub>S</sub>、曝气池混合液挥发性悬浮固体浓度X<sub>V</sub>、水力停留时间HRT以及实际进水有机物浓度S<sub>0</sub>;其中,<maths num="0001"><![CDATA[<math><mrow><mi>HRT</mi><mo>=</mo><mfrac><mi>V</mi><mi>Q</mi></mfrac><mo>,</mo></mrow></math>]]></maths>Q为进水流量,V为处理城市污水的反应器体积;2)利用城市污水处理厂日常实测的进水有机物浓度、进水质量和污水动力学参数进行数理统计,分别得到进水有机物浓度的标准差<img file="FSA00000024924700013.GIF" wi="129" he="52" />进水水量的标准差σ<sub>Q(t)</sub>以及污水动力学参数的标准差σ<sub>P</sub>;3)将所述实测值输入计算机,利用Matlab拟合一个维纳过程W(t)作为简化误差的补偿系数,并利用Matlab通过以下函数预测出水有机物浓度S(t):<maths num="0002"><![CDATA[<math><mrow><mi>S</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>&alpha;</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></msup><mi>S</mi><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>+</mo><mfrac><mi>&beta;</mi><mi>&alpha;</mi></mfrac><mo>-</mo><mfrac><mi>&beta;</mi><mi>&alpha;</mi></mfrac><msup><mi>e</mi><mrow><mo>-</mo><mi>&alpha;</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></msup><mo>+</mo><msub><mi>&sigma;</mi><mrow><mi>S</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></msub><msubsup><mo>&Integral;</mo><msub><mi>t</mi><mn>0</mn></msub><mi>t</mi></msubsup><msup><mi>e</mi><mrow><mi>&alpha;</mi><mrow><mo>(</mo><mi>x</mi><mo>-</mo><mi>t</mi><mo>)</mo></mrow></mrow></msup><mi>dB</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow></math>]]></maths>其中,t为时间,t<sub>0</sub>表示预测起始时间,e为自然常数;<maths num="0003"><![CDATA[<math><mrow><mi>&alpha;</mi><mo>=</mo><mfrac><mi>Q</mi><mi>V</mi></mfrac><mo>+</mo><mfrac><msub><mover><mi>&mu;</mi><mo>^</mo></mover><mi>H</mi></msub><msub><mi>K</mi><mi>S</mi></msub></mfrac><mo>&CenterDot;</mo><mfrac><msub><mi>X</mi><mi>V</mi></msub><msub><mi>Y</mi><mi>H</mi></msub></mfrac><mo>,</mo></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><mi>&beta;</mi><mo>=</mo><mfrac><mi>Q</mi><mi>V</mi></mfrac><msub><mi>S</mi><mn>0</mn></msub><mo>;</mo></mrow></math>]]></maths>σ<sub>S(t)</sub>为扰动系数且<maths num="0005"><![CDATA[<math><mrow><msub><mi>&sigma;</mi><mrow><mi>S</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></msub><mo>=</mo><msqrt><msub><mi>w</mi><mn>1</mn></msub><msubsup><mi>&sigma;</mi><mrow><msub><mi>S</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mn>2</mn></msubsup><mo>+</mo><msub><mi>w</mi><mn>2</mn></msub><msubsup><mi>&sigma;</mi><mrow><mi>Q</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mn>2</mn></msubsup><mo>+</mo><msub><mi>w</mi><mn>3</mn></msub><msubsup><mi>&sigma;</mi><mi>P</mi><mn>2</mn></msubsup></msqrt><mo>,</mo></mrow></math>]]></maths>其中w<sub>1</sub>、w<sub>2</sub>、w<sub>3</sub>分别表示<img file="FSA00000024924700018.GIF" wi="127" he="51" />σ<sub>Q(t)</sub>、σ<sub>P</sub>的权重,且w<sub>1</sub>+w<sub>2</sub>+w<sub>3</sub>=1;B(t)表示布朗运动,且dB(t)=W(t)dt。
地址 400030 重庆市沙坪坝区沙正街174号