河流水质数学模型分析.pptVIP

令： 可以得到： 这是一个递推方程，可以用矩阵方程表达： 式中，C和D是n×n维矩阵，分别为： 对于每一个断面的溶解氧，可以表达为： O2=C-1BO-C-1DL2+C-1(f+h) 式中： f=(f0 f1 … fn-1)T h=(h1 0 … 0)T h1=C0O20-d0L20 如果将BOD的表达式代入，形成偶合方程： O2=C-1BO-C-1DA-1BL+C-1(f+h)-C-1DA-1g 若令： U=A-1B V=-C-1DA-1B m=A-1g n=C-1BO+C-1(f+h)-C-1DA-1g 则： L2=UL+m O2=VL+n 这是描述多河段BOD-ＤＯ偶合关系的矩阵模型。Ｕ和Ｖ是根据给定数据计算的下三角矩阵。Ｕ是BOD（对BOD的）响应矩阵，V是溶解氧（对BOD的）响应矩阵。 4、正交曲线坐标系统 在一个给定的河段中，沿水流方向将河段分成m个流带，同时在垂直水流方向将河段分成n个子河段，组成一个包含m×n个有限单元的平面网格系统。 流管划分的原则是保持每一个流管流过相同的流量。流管的流量可以通过断面的单宽流量计算。单宽流量是局部水深h与断面平均水深H的函数： 式中，Q表示通过河流断面的总流量；a和b表示根据断面流量分布估计的参数，休姆根据观察数据给出了估值范围。 在平直河道中： 若50≤B/H≤70，a=1.0，b=5/3 若70≤B/H，a=0.92，b=7/4 在弯曲河道中，当50≤B/H≤100时，0.095≥a≥0.08,2.48≥b≥1.78。 5、断面累积流量曲线 根据单宽流量的计算结果作出断面累积流量曲线(如下图)。如果以纵坐标表示流量，将总流量分成m等分，对应的横坐标则表示流管的宽度。 由流线和断面线组成一个正交曲线坐标系统。假设第(i，j)个单元的长度是△xij,宽度为△yij，平均水深hij。 6、BOD模型 对任意一个单元（i，j），可以写出质量平衡方程。 由纵向水流输入、输出该单元的BOD总量为： 由纵向弥散作用输入、输出该单元的BOD总量为： 由横向弥散作用输入、输出该单元的BOD总量为： 式中，qi——第j个流带的流量； Lij——第ij个单元中的BOD浓度； 　　　Dij,kl——单元ij和kl间的弥撒系数； ——上下游相邻单元间的平均距离； ——横向相邻单元间的平均距离； Aij ,kl ——相邻单元ij和kl间的界面面积。 上式中 等项相当于相邻两单元间的BOD 浓度梯度， 即为单元i-1向单元i的弥 撒通量。 令： 并设在ij单位内的BOD衰减量为vijKd(ij)Lij、由系统外输入ij单位的BOD总量为 ,对ij单元写出质量平衡关系： 如果问题可以简化为稳态，则 如果将所有单元中的ＢＯＤ值写成ｍ×ｎ维向量： 将所有系统外输入也写成ｍ×ｎ维向量： 可以对这个河段写出矩阵方程： 　　　　　　　 GL=WL G是一个（m×n）×（m×n）维矩阵，称为变换矩 阵。G的各个元素g kl计算如下： 对 对 对 对 对 对其余l, 矩阵G的元素是流带流量、弥散系数、单元几何 尺寸及BOD衰减速度的函数，如果已知上述参数， 在给定外部输入的BOD（污染源）时，每个单元的 BOD值可以计算如下： L=G-1WL 式中，G-1是（m×n）×（m×n）维矩阵，称BOD 响应矩阵。 7、DO有限单元模型 与BOD模型相似，可以写出一个单元的DO平衡： 式中，Oij为DO单元的