水环境数学模型课件.ppt

对河口的BOD分布，即可以写出矩阵方程： 式中的 G 为 n 阶系数矩阵，其中的第 i 行,第 j 列的矩阵元素 Gij 可按下式计算： 当j=i 当j=i-1 当j=i+1 其余矩阵元素为零。 河口做模型 同样对于河口的氧亏，也可以写出矩阵方程。 矩阵 H 为 n 阶三对角矩阵，其元素可通过下式计算： 当j=i 当j=i-1 当j=i+1 F 为 n 阶对角矩阵，其对角线上的元素为： 将BOD的矩阵表达式代入氧亏的矩阵表达式中，得： 边界条件处理： i=1 时，出现了 Q0,L0 和 D0 ，在计算时可将这些已知项计入源项中。 i=n时，需要直到第n+1个河段的BOD和DO浓度及参数，对此可以将下游的浓度梯度视为 0 ，即令Ln+1= Ln ， Dn+1= Dn 该模型的解析解为： 多河段水质模型 BOD-DO耦合矩阵模型 0 i-1 i i+1 n Qi k1i li Li k2i ui Oi ksi ti Qi-1 Li-1 Oi-1 Qn Ln On Qo,i Lo,i Oo,i Qin,i Li n,i Oin,i 断面划分的原则： 断面形状变化处 支流或污水汇入处 取水口处 其它 多河段矩阵模型 根据连续性原理，写出每一个断面的流量 Q 平衡关系： 根据 S-P 模型写出由 i-1 断面至 i 断面之间的BOD衰减关系： 令 同时根据连续性方程，可得到： 令 可以得到： 该方程组可以用一个矩阵方程表达： 式中： 为 n 维向量， g1 ＝α 0 L0 BOD-DO耦合矩阵模型 根据 S-P 模型可以写出第 i 断面的溶解氧浓度： 同样，根据质量平衡原理可得到： 即 令 再令 最后得到 该方程用矩阵形式表达： 式中 C 和 D 为 n 维矩阵。 对于每个断面的溶解氧，可表达为： 将BOD的表达式代入，得到耦合方程 令 可以得到多河段BOD-DO的耦合矩阵模型。 上式中的 U 是BOD的响应矩阵， V 是DO的响应矩阵。 已知一维河流的输入,输出数据如下图所示设河流的饱和溶解氧值操作系统＝ 10mg/L 。试用多河段模型模拟河流的BOD和DO。 单位： Q － m3/s,L －毫克/L ， O －毫克/L ， k1 － 1/d ， k2 － 1/d ， t － d 0 I II III IV Qin,1= 0.5 L1 = 200 Oin,1 = 1 Qin,2= 0.3 L2 = 200 Oin,2 = 1 Qin,3= 0.4 L1 = 200 Qin,3 = 1 Qin,4= 0.5 L4 = 200 Qin,4 = 1 Q0 = 10 L0 = 2 O0 = 8 k1,0= 0.3 k2,0= 0.6 t0 = 1 k1,1= 0.3 k2,1= 0.6 t1 = 1 k1,2= 0.3 k2,2= 0.6 t2= 1 k1,3= 0.3 k2,3= 0.6 t3= 1 Qo,1 = 0.2 Qo,2 = 1 Qo,3 = 0 Qo,4 = 1 解： 第一步，计算 i=1 ～ 4 的α i ，β i ，γ i, δ i 第二步，计算矩阵 A ， B ， C,D 的元素 第三步，计算逆矩阵，求出BOD和DO的响应矩阵及向量 第四步，利用 U,V 计算各断面的BOD和DO浓度。 含支流的河流矩阵模型 假设主流含有 1,2 ，… i … n 个断面，支流含有 1,2,第一… m 个断面，在主流断面汇入主流对支流写出矩阵方程，计算支流最下游断面 m 的水质，将支流作为污染源计入主流的矩阵方程。。 其它水质模型 QUALL-2 模型 模型描述的主要成分之间的作用 1、复氧作用； 2、河底生物 (包括底泥) 耗氧； 3、碳化阶段降解耗氧； 4、光合作用产氧； 5、氨氮氧化耗氧； 6、亚硝酸盐氮氧化耗氧： 7、碳化合物的沉淀； 8 、浮游植物对硝酸盐氮的吸收； 9 、浮游植物对磷 (磷酸盐) 的吸收； 10 、浮游植物呼吸产生磷（碳