数值模拟在嘉定蕴南中部地区引调水试运行准备工作中的应用

2012年9月第9期城市道桥与防洪相关专业317数值模拟在嘉定蕴南中部地区引调水试运行准备工作中的应用包小燕，李念斌(1．上海市嘉定区水文站，上海201800；2．上海嘉定水务工程设计有限公司，上海201：800)摘要：该文通过结合上海市嘉定区蕴南中部地区引调水试运行准备工作，主要阐述了数值模拟的发展现状，介绍了一维河网非恒定流的数值解法、主要步骤，并结合工程实际进行了模拟运算。关键词：数值模拟；水质分析；引调水中图分类号：TV67、TV68文献标识码：B文章编号：1009—7716(2012)09-0317-03模型、水环境容量模型等)的基础和前提。水动力模0引言型的建立涉及控制方程组的简化、方程组的离散上海市嘉定区蕴南中部地区涉及南翔、江桥、和求解、初边条件的确定、模型的率定和验证等一黄渡三镇，是嘉定区经济最发达的地区之一，但也系列问题。平原河网模型不同于单一河流的特点，是嘉定区水环境质量最差的地区。为此嘉定区水在于河网错综复杂性，以及由此带来的方程组离务局于2007年组织编制了《嘉定蕴南中部地区水散和求解上的困难，这是多年来人们研究河网问环境治理——引调水专项规划》，经过2008年至题的一大难点。平原河网水动力数学模型目前大2012年的建设，规划中的水利控制工程已基本建体可以分为节点一河道模型、单元划分模型、混合设完成，嘉定区水务局拟定于2012年10月进行模型、人工神经网络模型，以及蒙特卡罗随机游动调水试运行。针对规划中设置的调水试运行方案，模型，其中最常用的是节点一河道模型。使用数值模拟的方法进行运算对比，选取合理的2河网非恒定流数值解法水质分析比对点，为调水试运行提供参考依据。2．1水动力方程1河网水动力数值模拟发展现状2．1．1基本方程组数值模拟也叫计算机模拟。它以电子计算机水流在平底、棱柱形明渠中一维非恒定流动为手段，通过数值计算和图像显示的方法，达到对的基本方程组一圣维南方程组，包括反映动量守工程问题和物理问题乃至自然界各类问题研究的恒定律的运动方程和反映质量守恒定律的连续性目的。20世纪50年代以前，数值模拟的基本理论方程：已经建立，运用这些理论也曾解决过一些简单的警+=g工程问题。但是，这些基本理论真正应用于解决工程问题却是在电子计算机发明以后。1952年一+(仅)一=o1954年，E．1saacson、J．J．stoker和B．A．Troesch首次式中，A为河道内任意断面的面积，ITI；Q为断建立了俄亥俄河和密西西比河部分河段的数学模面水流量，1TI；t、X分别为时间与距离的坐标；Z为型，并进行了实际洪水过程的模拟。到上世纪60水位，m；g为重力加速度，m·S-2；q为旁侧人流，m；年代中期，为了解决各种各样的设计和规划问题，Rh为水力半径，in；Ck为谢才系数，ITI沈·S；a为动水动力数学模型再次得到重视。随着计算机技术量校正系数。的发展，水动力数学模型的模拟功能也大大增强，2．1．2方程组求解方法可以对整个流域、洪泛区、已建或规划中的水利工如图1所示，利用Abbott六点隐式差分格式求程进行系统模拟。水动力模型是其它模型f如水质解圣维南方程，在每个网格节点按顺序交替计算水位和流量。该格式为无条件稳定，可以在相当大的收稿日期：2012—07—06作者简介：包小燕(1975一)，女，上海人，工程师，从事水文水Courant数下保持计算稳定，因此，可以取较长的时质研究工作。间步长以节省模型运算时间。 318相关专业城市道桥与防洪2012年9月第9期3．3方案设置根据《嘉定蕴南中部地区水环境治理——引调水专项规划》的方案规划，将通过沿蕴藻浜南侧沿线设站，提抽内部河道水进入蕴藻浜，待内部河水图1Abott格式水位点、流量点布置示意图位降低后，强行引入苏州河清水实行水体交换与2．1．3定解条件流动。并利用东侧现有的连浦闸站向内引水，走马定解条件为水流的初值与边界值，具体如下：塘水闸开启。这样形成“北排、南引”的有序水流大水流初始条件：t=0，Z(x，t)=z(x，0)，Q(x，t)=Q(x，0)。边格局。界条件：当x=0时，z(x，t)：z(O，t)；当x=L时，Z(x，t)=该方案沿蕴藻浜在浅江、封浜、横沥设泵闸，Z(L，t)。浅江泵闸流量为5m3／s、封浜泵闸流量为10m3／s、河网的内边界条件比较复杂，其中河网支汊横沥泵闸流量为10m3／s。为形成良好的换水条件，节点的处理最为关键，支汊点应满足z1=z2=⋯⋯与该方案配套的工程有：西线沿盐铁塘在狗颈骨河、北线沿蕴藻浜在花园江处设闸控制。沿蕴藻浜=Zn，Qi=O。控制工程如图2所示。2．2水质扩散方程皿根据物质在水体中对流扩散的过程，考虑污染物的对流扩散与线性降解，其基本方程为ll-2l：+一鲁lAD】=-AKC+式中：C为污染物浓度，mg／L；D为污染物弥散系数；A为断面过水面积，133；Q为流量，m3／s；K为降解系数，s；C为污染物的点源浓度，mg／L；q为污染物的点源流量，m3／s；X为空间步长，m；t为时间步长，s。3区域概况及方案设置图2蕴南中部地区引调水控制工程分布图3．1区域情况嘉定区位于上海市的西北部，市级河道蕴藻4数值模拟方法浜横贯其南部，在穿越沪宁铁路后折向西南与苏4．1河网概化州河交汇。由此，嘉定区水系形成了界线分明的南在充分掌握详细的天然河网、湖泊的水动力、北两片。其中，蕴藻浜以南地区的面积约为106．4水文资料的基础上，以主干河道为基础，按照河网km，其范围是：东以嘉定行政区域为界；西以蕴藻概化的基本原则，对蕴南中部地区河网进行合理浜为界；南北两侧各以苏州河和蕴藻浜为界，形成地概化。所涉及的主要河流有：吴淞江、蕴藻浜、盐了一个相对独立的区域。这一区域俗称嘉定的蕴铁塘、浅江、封浜、横沥、西虬江、中槎浦等。南地区。4．2边界条件3．2水环境现状4．2．1外边界条件该区属平原感潮河网地区，内部水系绝大部平原河网地区在洪水期往往同时遭遇平原区分天然形成，且弯曲流长；河道水流流向不定，主产流和潮汐等多重影响，其边界条件大致分为3要受南部感潮河流苏州河的水位影响，水流方向类：水位(潮位)边界、流量边界(水位流量关系边由北向南、由西向东居多。但由于区内河道向苏州界可以转化为流量边界)，以及计算区域的净雨输河排水受到限制，阻断了区域内外水体自然交换人。该项研究采用吴淞江黄渡站水位及蕴藻浜闸的道路，改变了地区原先“由北向南流”的水力特内站水位为外边界。性。此外，由于铁路、公路的建设，交叉建筑物也阻4_2．2内边界条件碍了水流的畅通，使该地区基本成为死水区，水质在内部边界条件处，圣维南方程组不再适用，情况常年处于劣V类水。必须依据其水力特性作特殊处理，以保证河道追 2012年9月第9期城市道桥与防洪相关专业319赶系数的连续传递。4．3时间空间步长的选取考虑到模型计算稳定性要求及模型运行时间予·‘；，／，要求，水动力模型的时间步长可以取5～15min，，，l／j／空间步长取50ITI～1500nl不等。／{4．4水质分析点选取——CSlI一一一一(：S2I该项计算主要选取封浜、横沥、西虬江的汊点⋯⋯CS3l⋯-·CS4l断面处，由于封浜、西虬江被铁路阻隔，在封浜被i铁路阻隔的南北、西虬江被铁路阻隔的东西分别；0122436柏6072选取水质分析点，共计20个，详见表1所列。时间／h图3横沥选取断面72hCOO表1水质分析点一览表相对减少量随时间变化过程曲线图5建议在蕴南中部地区引调水试运行前，针对一维河网非恒定流复杂的水流特性，建立了嘉定蕴南口／0＼姆靛靶o8∞∞∞∞伸0中部地区的水量水质模型，选取了水质分析点，并进行了数值计算，为试运行提供了参考依据。现针对将进行的试运行，提出以下建议：(1)做好水质取样分析工作，掌握蕴南中部地区的水质在引调水前、后的情况，尤其是数值计算涉及的这些分析点，为模型的校准提供数据。(2)及时采集引调水的数据，分析试运行中存在的问题，为蕴南中部地区引调水常态化做好准4．5模拟结果备。由于溶解氧、氨氮等水质指标会沿程和随时间变化，所以采用综合的CODer作为计算指标，用(3)在已建立的数值模型：基础上，通过试运行f换水效果的比较。限于篇幅，以横沥河为例，72的数据，对此区域的水利控制工程调度问题开展hCOD相对减少量的平均值(％)见表2所列，随时探索。间变化过程曲线如图3所示。表272hCOD相对减少量的平均值一览表(单位：％)参考文献【1]欧剑，马进荣，张行南，张晓艳．大通至长江13整体水动力模型IJ1．河海大学学报：自然科学版，2009，(5)：258—262．[2】匡翠萍，李行伟，刘曙光．大规模围垦对香港维多利亚港水动力环流的影响【J1．同济大学学报(自然科学版1，2009，(2)：176-181．[3]郜会彩，李义天，何用，陈建．改善汉阳湖群水环境的调水方案研究【J】．水资源保护，2006，(9)：41—44．【4]黄伟，何用，李荣，陈钢锋．感潮地区引清调水方案研究IJJ．水电能源科学，2004，(6)：19—22．