自由水面湿地系统耗氧系数K1几种实验室求法计算和讨论.docVIP

自由水面湿地系统耗氧系数K1几种实验室求法计算和讨论 　　摘要 　　 BOD（生化需氧量）是水质污染研究中一项重要的氧平稳指标，在描述和评价水质有机污染时常采用。当有机污染物质进入水体后，如果向水中供给溶解氧，并随着耗氧速度的变化，水体的复氧速度亦发生改变。此时，如果耗氧速度大于复氧速度，则水质将变成嫌气性，在嫌气微生物作用下，有机物可发生腐烂，发酵，变成有害恶臭物质，而使水质受到污染。如当复氧速度大于耗氧速度时，则可向相反方向发展，水质可以逐渐恢复，即自解作用，因此对受有机污染形成的水域，研究其耗氧和复氧的速度变化是水质保护工作中的一个重要问题。本文就计算耗氧系数（K1）的几种方法进行讨论。 　　关键词：自由水面湿地 BOD 好氧系数 计算 　　一、Phelps的耗氧作用定律 　　 BOD的耗氧过程分为两个阶段，第一阶段为碳化阶段，在温度20℃下耗氧主要在10天以内，第二阶段为硝化阶段，是一个缓慢的过程，在温度20℃下耗氧主要在7－20天以内。河流水质影响主要是碳化阶段的耗氧，所以实验室K1是以10天BOD进行计算的，用来描述第一阶段或碳BOD的一级反应动力学以公式为： 　　 　　式中：K1――是碳BOD的反应速度常数； 　　t――是反应时间； 　　 这里L0是在反应时间为O时BOD的初始值。在实际工作中，我们不是通过随着时间BOD的减少来测量BOD，而是通过随着时间累积的氧耗来测量BOD，设这种累积的氧耗是Y。 　　二、BOD原始资料的收集 　　 本文用于计算的BOD来自某基地漫流地段三个断面采集的水样，在实验室利用BOD测定的标准方法，在20℃下做1－10天的序列培养，分别测定1－10天BOD值而得。 　　三、计算K1值的方法及实例 　　1最小二乘法 　　 把最小二乘法用于一组BOD数据，以取得在某种意义上应是最精确的降解速率（K1）和最终BODL0，这时BOD基本上已不随时间变化，使用最小二乘法的基本思想是：首先找到近似于线性函数的关系式，使离开均值的平方总和为小。依据此原则，本文选用二个线性关系式然后使用最小二乘法。 　　1）h最小二乘法 　　 在方程中，，令是一个可以确定的近似值，K1是真实的降解速度，而h是一个很小的变量，将上式代入BOD方程便有： 　　 （a） 　　因为h值极小，所以有： 　　 　　即式（a）变成： 　　 （b） 　　这里： 　　方程（b）是一线性方程，可以用最小二乘法来定出a、b。 　　设：为n个实验数据，Y1为公式所算得数据： 　　 　　必须使真实值和计算值的平方差和为最小。所以的一阶导数应为0： 　　 （c） 　　同理： （d） 　　其中： （e） 　　 （f） 　　使导数为零，（c），（d）为： 　　 　　解得： 　　 　　 　　 　　 为了求解a和b，首先需要使用函数的值；即使用K1的偿试值计算每一个时间t的，计算出的加和值，并将这些加和值代入上面两式求出a、b。 　　 为了得到K1和L0值，在推导最小二乘法的过程中必须满足h极小，所以选择的K1的尝试值必须足够地接近于K1的真实值，在实际运算中，当h值小于0.0001时，将会满足这一条件。 　　 由a和b确定L0与h，如果h0.0001，结果的误差较大，说明假设的这个不太适合，我们用作为新的起始值；重新开始这个计算过程，为了使h足够小，必须重复此过程。 　　以进水的计算机叠带过程为例： 　　 　　 　　 　　最终选K1=0.2136 　　 同样，我们得到5－1的K1=0.22615－6的K1=0.1053 　　（2）取对数最小二乘法 　　由BOD反应动力学推导公式： 　　 　　式中：L－剩余BOD（mg/L） 　　L0－t=0时的BOD（mg/L） 　　t－时间（天） 　　设： 　　则有：Y=a+bX 　　计算 　　 　　 　　 此方法认为总BOD用第10天的BOD代替，用第10天的BOD分别减去前9天的BOD，得到前9天剩余的BOD，然后取对数，表1 为3个断面对数剩余BOD值。 　　下面分别计算进水，5－1,5－6的K1值 　　进： 　　 　　表1取对数最小二乘法计算表 　　 　　 　　 　　5－1： 　　 　　同样可算出5－6的b=-0.3932K1=0.3932 　　（3）两点法 　　Rhm叙述了两段时间内BOD值之间的关系，第二段时间恰好等于第一段时间的2倍。他所采用的以下两个表达式是根据实数推导出来的。 　　 　　式中： 　　X－在t时间（天）测出的BOD，即BODt 　　Z－在T时间（天）测出的BOD，即BODT 　　T、t－分别为第二次，第一次测定BOD时间（天） 　　T=2t，一般t取5天，结果见表2。 　　表2两点法计算表 　　 　　 　　 　　表2可证明，尽可能使用高的t值是正确的。在t=5，T＝10（天）时所