黄河中下游污染物扩散模型研究rdEdition.docVIP

黄河中下游污染物扩散模型研究 Xxx1,2，Xxx 3，Xxx 2，Xxx 1 （1.西安交通大学电信工程学院，陕西 西安 7100049） 摘要：黄河流域的突发性水污染事件是应急工作的重要内容,本文以圣维南方程组为基础,建立了黄河中下游流域一维污染物扩散模型，结合黄河中下游流域的特点进行关键参数率定，使用有限差分格式进行离散化，快速有效地实现了水体演进与污染物扩散的预测。实验结果表明，该模型的可用性和可靠性高，可应用到实际的工程中。 关键词：污染物扩散，圣维南方程组，有限差分法，参数率定 中图分类号：TP 212.9 文献标识码：A 文章编号： 近年来，黄河流域的突发性水污染事件造成的严重后果引起了人们的重视，传统的监测管理方法由于其局限性已经无法满足对突发性水污染事件应急工作的需要[1]。如何有效、科学地对黄河的突发性水污染事件进行监督和管理，降低突发性事件对流域人民的生产、生活的影响，成为目前急需解决的热点问题[2]。 在污染物扩散模型及其预报方面，国内外学者做了大量的研究工作，在模型建立与求解方面取得了众多有意义的成果[3]。在实际应用中，要考虑到河床变化等流域特点对水体演进、污染物扩散预测所产生的影响，以求获得更高精度的计算结果。 本文在国内外河道水体演进以及污染物扩散预测方面的研究成果的基础上，结合黄河预警预报系统中对计算的高效、准确等要求以及黄河中下游特点，对黄河污染物扩散模型的建立和求解进行了研究，并进行了实验数据验证，取得了较好的效果。 1. 模型建立 1.1. 研究区段 黄河是我国第二大河，是西北、华北地区最重要的水资源。黄河因其高含沙量而闻名于世，是公认的水动力条件复杂的河流之一。其流程长、汛期复杂多变，尤其是泥沙的沉淀不断改变着河床的形状，难以建立统一模型。 黄河中下游段多在平原，比降低,约为0.013%,流速变化相对平稳,使得建立一个精确的模型成为可能。本文所研究的区段起自于小浪底坝，至利津站入海区段，总长764公里，平均河面宽度为5百米。区段内有8个主要监测站提供常规的日常断面监测数据。这些数据是必需的，是对该河流进行研究的主要元素，关联着污染物出现时间，浓度，类型，河流水文信息。检测数据由YRCC提供，包含05与06年间的大约12,800个记录，平均每个站点有1,600个记录。数据涵盖了实际中的旱季与雨季情况：流量变化范围较大，从最小的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 … … … … … … 这些监测数据可用于几个途径。一是作为边界条件与初始条件；二是作为检验数据；三是用于重要参数的校准。 1.2. 选取模型 对河流建模包含两个部分：概念模型（数学模型）与计算模型（求解方法）。概念模型由描述该河流的观察资料、建模数据与数学方程组成，它同样包含初始条件与边界条件。计算模型也就是用来求解概念模型的数学方法，按照它来组织概念模型的程序实现。 首先要确定的是模型的维度。水力模型表征在时间上的三维空间，而在本文中我们选择一维模型，主要原因有以下几点。 第一，我们获得的历史数据是一维的，选择二维三维模型需要更多的数据，但这些数据通常较难获得；第二，我们希望计算得到的数据是一维的，指对污染物扩散的预测我们更关注的是它什么时候到达一个特定的断面而不是去研究在此断面上污染物的浓度分布；第三，研究区段长度远远大于河宽，这也使得一维模型客观可行；第四，河床变化剧烈，选择二维或者三维模型对此参数的刻画将十分困难，不仅仅在于实际测量的难度，同样地，模型也要随时更新数据以达到与实际相符的目的。 于是，一维模型成为研究的基础与具体内容。 为了定量研究污染物随水体的迁移过程，首先要先掌握水体演进的过程，将模型分为两部分,分别为水动力学模型和污染物扩散模型。 1.3. 水动力学模型 水动力学模型考虑了动力方程中的各种影响因素，多用连续性方程和动量方程描述。动量方程即纳维斯托克斯方程或雷诺方程，而在此领域应用最多的则是圣维南方程组或与其相似的方程式。水动力学模型在较大程度上克服了水文学方法的缺陷，对水文资料的依赖性较小，计算精度较高，适应性