深圳河减污的水力条件计算和分析.pdfVIP

·182· 中国水利学会2003学术年会论文集 深圳河减污的水力条件计算与分析 郑金海1昊良冰2童朝峰1 (1，河海大学海岸及海洋工程研究所2．深圳市治理探圳河办公室) 摘要采用拉格朗日法研究深圳河上游的罗芳污水处理厂和滨河污水处理厂所捧出污水 沿深圳河的运动过程。建立了深圳河非恒定流数值模型，计算分析了河口建闸和上游调水对 深圳河减污的水力条件。结果表明，河口建闸对于减轻深圳河的污染情况作用不大，增加上 游来水量对减轻深圳河的污染状况是有利的。 关■词深圳河减污河口建闸数值模型 深圳河是深圳经济特区与香港特别行政区之间的界河，位于珠江口东侧，在东经1140一114。12’50”， 北纬22027’一22。397之间。属海湾水系。深圳河发源于梧桐山牛尾岭，三岔河口以上河道走向为东南向， 以下河道折向西南，流经深圳市区后注入深圳湾，河道全长约37km。其主流为沙湾河，在三岔河口与莲塘 河汇合后为深圳河干流，支流梧桐河在罗湖桥下汇入，布吉河在渔民村汇人，从三岔河口到深圳河口的干 流河段为感潮河段，长约16kin。近年来，随着深圳市的城市建设不断发展，城市人口数量急剧增加，向深 圳河排人的污染物不断增加，对周围环境影响极大，与深圳市高速发展的经挤和美丽的市容极不协调。 受深圳市水务局的委托，河海大学海岸及海洋工程研究所在深圳河治理三期工程完工的前提条件下， 建立数值模型从水动力条件上对如何利用潮汐动力来减轻深圳河的污染问题进行分析研究，通过计算深 圳河上游的罗芳污水处理厂和滨河污水处理厂所排放的污水在深圳河运移的变化情况，对减轻深圳河污 染的水力条件进行研究和评价，提出减污的可行方案。 一、数值模型 程组进行有限差分运算——采用Preissann格式，得单一河道方程——以各断面水位及流量为自变量的 差分方程组，再利用河道上下游的边界条件求解差分方程组，计算出各断面的水位、流速和流量等水力 参数…。 1．1基本方程 描述水流在明渠中运动的一维非恒定流基本方程为： 等+风警=口 (1) 警+2u等+(卵一舰2)警U2罢+g与垆=o (2) 式中：I为时间坐标，z为空间坐标，Q为流量，Z为水位，u为断厩平均流速，n为糙率，A为过流断面面积，B 为主流断面宽度，见为水面宽度(包括主流断面宽度B及仅起调蓄作用的附加宽度)，R为水力半径。口为 深圳河减污的水力条件计算与分析 ·183· 旁侧人流流量。此方程组属于二元一阶双曲型拟线性方程组，采用有限差分法求数值解。为了加大时间 步长，节省计算时间，在求解上述方程组中，采用四点线性隐式差分格式…。 1．2数值离散 465m。数值模型的断面距离 数值模型范围包括河口桩号0+000至上游13+465的河道，总长为13 的位置布置在断面上(图1)。 三岔河口●————T————1_—————————T_—————．r—————]—————●深捌I河口 断面号 1 2 3 …… i i+I 270 271 图1河道节点及变量布设 最终得到的差分万程组卿r， 1+C,ZT：?一口?“+Q?l=D， (3) C,Z7 E。Q■+G；QT：?一只刃“+I,z7：?=皿 (4) 2哦．。+学(噜掣)h ”急协‰。+学‘可n2lul