黄河中下游污染物扩散模型的研究-3rd-Edition.docVIP

黄河中下游污染物扩散模型研究 Xxx1,2，Xxx 3，Xxx 2，Xxx 1 （1.西安交通大学电工程学院，西 西安 7100049 摘要：黄河流域的突发性水污染事件是应急工作的重要内容,本文以圣维南方程组为基础,建立了黄河中下游流域一维污染物扩散模型，结合黄河中下游流域的特点进行关键参数率定，使用有限差分格式进行离散化，快速有效地实现了水体演进与污染物扩散的预测。实验结果表明，该模型的可用性和可靠性高，可应用到实际的工程中。 关键词：污染物扩散，圣维南方程组，有限差分法，参数率定 中图分类号：TP 212.9 文献标识码：A ： 52m3.s-1到最大值4240m3.s-1 。这些检测记录是一维的，也就是说，在给定的截面和时间，对于类似水位或者流速这样的水文数据，监测结果经平均后只有一个数值。每天对水位与流量进行一次测量，每4天再加一次包含截面面积的测量。如下所示。 表1 监测数据的格式 河名 Channel 站名 Station 时间 Time 水位 H (m) Water level 流量 Q(m3.s-1) Water flow 测流面积 A(m2) Cross-sectional area 黄河 花园口 2005-01-01 08:00 91,27 445 黄河 花园口 2005-01-01 14:18 91,21 404 600 … … … … … … 污染物扩散模型 污染物扩散模型是对污染物扩散的模拟，是对水体中污染物随时间和空间迁移转化规律的描述，即用于描述水体中污染物质与时间、空间定量关系的数学方程。目前应用广泛的是确定性数学模型[5]，可以根据一组确定的输入条件，求解出对应的值。 污染物扩散模型采用一维平流扩散方程表示，其根据是悬浮物与溶解物的一维质量守恒定律，其方程如下： 模型求解方法 数值解法 数值计算方法将给定的偏微分方程按照一定的离散格式离散化，获得对此方程的近似的解决方案。在这里将河流离散成一个个较小的线性计算单元。广泛使用的标准离散格式有有限差分法、有限体积法和有限元法。 有限元法是将整个方程视为一系列相互连接的局部形状函数的总和，局部形状函数用差值函数进行局部模拟单元体内的特征点，从而将偏微分方程组离散化为代数方程组。该方法适用于模拟复杂的几何学问题以及变化剧烈的结构，而在机械工程学中应用广泛，相应的它也需要耗费更多的时间用于计算，且编程实现较难。在本例中，不变的河流环境使我们没有必要在此耗费大量的资源。 有限差分法是在离散的网格节点上把微分方程的各阶偏导数用差商取代，使微分方程变为代数方程。遵循基本差分规律，局部区域特征被附加到每个离散点上。 有限容积法将计算域划分为许多互不重叠的控制容积，使每个控制容积内部都包含有一个网格节点，在每个控制容积内对微分方程进行积分求解，就可以得到一组包含有网格节点上待求变量值的离散化方程。结合有限元的有限差分，有限容积法就称为差分方程的自然扩展。然而对于一维水质模型，水体密度恒定使得按照有限差分法与按照有限容积法的离散方程格式一致。 于是本文使用有限差分法作为计算模型。水动力学模型是一个双曲型方程组，用有限差分法对方程进行离散，将河段离散成N个点，离散点越多，计算精度越高，但模型需要更多的计算时间。 如图1所示，横坐标为空间步长，纵坐标为时间步长，表示河段上第j点在n时刻的流量，相应的有其他因变量如表示水位等等。 图1 离散点关系图 将连续方程 离散为如下格式方程： 其中 ;;; ;； 同理对于动量方程 得到如下格式方程： 其中 ；； ；； ； 污染物扩散模型离散方法与上述的过程相同，不再赘述。 加入边界条件，与(4)、(5)联立建立方程组并求解，就可以得到河段上任意点在任意时刻的流量和水位值。水动力学模型求解的结果作为水质方程的输入，进行模型的污染物扩散模拟。 特殊处理 上述模型在水资源领域应用广泛，为了能在黄河中下游河段获得更好的计算结果，需要对模型进行一些特殊处理。 过水断面面积 在模型计算过程中，需要过水断面面积作为输入，然而它是随着时间、地点不断变化的，尤其受水位的影响显著。在黄河这样的多泥沙河流中，河床每年都因泥沙沉积而发生明显变化，这对过水断面面积的计算产生了较大的困难。问题的关键就在于如何在模拟时，确定每一个过水断面的面积值。 参考黄河中下游几个监测站点过水断面面积的数据，引入一个梯形模型来模拟黄河的河道形状，如图2所示。 图2 河道截面模型图 阴影部分为过水断面，O点处高程记为为，由梯形面积公式得： 记然而希望在此式中仅有H一个变量即得到，这就需要我们率定系数和的值或者是将其变化成关于H的表达式。进一步地，将E与之间的距离代入公式消除，将此差值记为。在excel中对数据进行分析，不断调