理想凝聚沉淀试验数据的整理与分析_secret

1、理想凝聚沉淀试验数据的整理与分析摘 要：对凝聚性沉淀试验数据的整理是设计沉淀池的一个重要步骤。试验首先要测出不同深度、不同沉淀时间下悬浮物的去除率，然后采用插值法绘出去除率等值线，选定某一深度，计算出不同停留时间下得到的总的去除率，最后将结果点绘成停留时间 t与去除率 C的关系曲线。本文对matlab编程整理凝聚性沉淀试验数据进行了研究。关键字：凝聚沉淀试验 沉淀池 去除率 Sorting out and Analyses of ideal Coagulation Sedimentation Experimental Data Abstract：In sorting out the exper

2、imental data of coagulating sedimentatation，manual calculation and curving are usually used to find out the relationship among the removal rate C，the surface load 0 and the residence time t.If Matlab program is used for the sorting out of the data，artificial effects can be minimized，full use can be

3、made of the data already available to get the results of calculations quickly，efficiently and accurately，and the results can be adjustedanalyzed and evaluated directly，which will optimize the design dataKeywords：coagulation；sedimentation； coagulation calculation；surface load：residence time对凝聚性沉淀试验数据

4、的整理是设计沉淀池的一个重要步骤。试验首先要测出不同深度、不同沉淀时间下悬浮物的去除率，然后采用插值法绘出去除率等值线，选定某一深度，计算出不同停留时间下得到的总的去除率，最后将结果点绘成停留时间 t与去除率 C的关系曲线。本文对matlab编程整理凝聚性沉淀试验数据进行了研究。1 理想凝聚沉淀的数学模型从实验数据可以看出：悬浮物去除率C是沉淀深度ho、停留时间L的函数，记作C=f（h0，t）。同时，可以根据公式(1)计算截流沉速。0hot (1)式中：0-理想沉淀池的截流沉速；ho-沉淀深度；t-颗粒在沉淀区中的停留时间。由于表面负荷在数值上等于截流沉速，所以又可以换算出表面负荷与去除率的函

5、数关系，点绘出二者的关系曲线。因此，C又可以看作0和t的函数C=g（0，t）。当0，t，ho中的任意两个参数确定后，C的值也随之确定。根据C与t，ho的测定结果，利用数值插值的有关计算方法进行三次插值，可以得到任意时间、任意深度的去除率计算结果。并且结果可用三维或二维坐标图表示出来。2 数据处理21 Matlab方法Matlab是集数学计算、图形处理和程序语言设计于一体的著名数学软件。与其它高级语言相比，Matlab提供了一个人机交互的数学系统环境，并以矩阵作为基本的数据结构，可以大大节省编程时间。 22 对凝聚沉淀的分析为了分析问题的方便，把最初的凝聚性悬浮物浓度当作100，以后所测得的去除

6、结果都用相对去除率表示。221 计算某浓度下的去除率假定池深为ho，选定不同的停留时间，计算出某一时刻水中悬浮物去除浓度。总的去除百分率为：C=(1-C0)0co（Ui0）d Ci （2）式中：C0-所有沉速小于0的颗粒重量占原水中全部颗粒重量的百分率；0-理想沉淀池的截流沉速；i-小于截流沉速的颗粒沉速；Ci-所有沉速小于U；的颗粒重量占原水中全部颗粒重量的百分率；dCi-具有沉速为认的颗粒重量占原水中全部颗粒重量的百分率。对于上式中积分的处理，可以近似地认为：式中：h0-沉淀池深度；hi-去除率等值线图中相邻两条等值线的中点所对应的深度；Ci-去除率等值线图中相邻两条等值线所代表的去除率之

7、差。Ci取值越小划分得越密，积分与累积求和的差值越小，精度就越高。利用Matlab编程可以很方便地实现这一点。已知池深为h0，需要画出0t时段上的总去除率曲线，方法如下：根据前述插值拟合的结果，可以找出在不同时刻、在h0深度处沉速大于或等于0的已全部去除的颗粒的去除率CO。对于某一时刻t，根据插值拟合的结果，找出不同深度处的去除率，代人（3）式，可以求出在h0深处，未能全部去除悬浮物的去除率。总去除率由以上两部分相加即可求得。假定不同的沉淀深度，重复上述步，可得出不同深度、不同时间下的总的去除率。这一计算结果与根据实验结果得出的二维值拟合结果不会相同。前述插值计算的结果反映的是某一时刻、某一深

8、度该点的去除率。而后者是考虑了沉淀池的出水实际上是某时刻不同深度处的出水相混合的结果。假定不同的沉淀深度，会有不同的总的去除率。或者说，将不同深度的去除率作了平均化计算，取其平均数。 3 计算实例例：凝聚性悬浮物的质量浓度为400mgL。静置沉淀试验所得的沉淀时间和取样深度以及相应的悬浮物去除率如表1；表1悬浮物去除率取样深度m不同沉淀时间（min）的去除率/%510204060901200.61415060677273761.22193345586270741.8315313854596371根据上表的试验数据，试确定采用平流沉淀池的悬浮物去除率、表面负荷和停留时间之间的关系。此时的初始条件

9、为：t0时，任意深度处的去除率 C0；h0（而t0）时，去除率C100，插值的结果用三维图形表示如图1所示。从图1可以看出，池深越浅、沉淀时间越长，去除率越高。例如h0183m时，计算结果如图2。由图2可以看出：当沉淀时间在040min以内时，去除率随时间增长很快；超出该范围时，去除率随时间增长缓慢。因此将时间控制在40min左右时，既经济又高效。超出该范围延长时间，对去除率的影响不大。由于0=h0/t1.83t，所以利用Matlab进行简单地换算即可计算出表面负荷0与去除率C之间的关系，绘出关系曲线，如图3所示。至此我们得出了在池深为1.83m时，在不同停留时间或不同表面负荷时所对应的去除率

10、。图3Oc关系曲线如果改变有效池深，可以得出另外一条t-C曲线和0一C曲线。例如，池深h015m时，t一c曲线和0一c曲线如图4、图5。比较图2和图4，可以发现：相同停留时间，池深浅时，去除率高，并可以在曲线图上准确地找出某一深度时t-C的准确对应关系。31 计算结果的验证假定有效池深1.83m，停留时间35min，则相应的表面负荷为：075.3m3（m3d）。利用Matlab计算出的总去除百分数为64.4，手工计算的总去除百分数为60.9。如果假定池深和停留时间不同，还可以计算出不同的去除率。 利用Matlab计算与手工计算有一定的出入，其主要原因分析如下：计算机计算的去除率略高于手算的去除率，主要是由于在水深浅时，计算机的程序中对深度作了较为精细的划分。且在理想条件下认为h00时，C100。而手算时，并没有将这些因素全部考虑进去。原始数据中所取的相邻深度相差较大，而利用Matlab计算是为了有较高的精度，深度划分、时间划分都较细，而手算很难做到