天角潭水利枢纽工程近、远期供水方案对北门江河口水环境影响研究.docxVIP

? ? 天角潭水利枢纽工程近、远期供水方案对北门江河口水环境影响研究 ? ? 张志广，李 准，金 弈，王 浪，邓景军，赵高磊 (1.中国电建北京勘测设计研究院有限公司，北京 100024；2.海南省水利电力集团有限公司，海南 海口 570203；3.四川大学水力学与山区河口开发保护国家重点实验室，四川 成都 610065) 0 引 言 入海河口作为一种特殊的水体，往往成为人类活动密集的地区或经济中心[1-2]。如，我国的长江口、珠江口和黄河口等。近年来，人类活动对河口带来的一系列生态环境问题已成为河口研究的热点。Andrew等[3]采用Mike21模型得到旧金山河口水质恶化的原因为上游水资源的过度开发利用；赵高磊等[4]研究了北门江径流变化对河口段水环境的影响；周荣香等[5]采用一、三维耦合水流盐度模型研究了珠江河口咸潮上溯机制；Burchard等[6]提出河口环流强度随盐度梯度的增大而明显增大；潘明婕等[7]通过采用三维水流盐度模型研究了台风过境期间磨刀门水道咸潮上溯的动力机制；姚姗姗等[8]采用二维数学模型研究了海棠湾规划水系水体交换能力。天角潭水利枢纽建成后，改变了坝下河道的水文过程，丰水期水量有所减少，枯水期水量有所增加，这些变化可能引起北门江河口咸潮上溯、盐度、河口海域冲淡水区变化等生态环境问题。为此，需要定量分析天角潭水利枢纽运行后对河口水环境的影响程度和范围。 本文以天角潭水利枢纽工程近期、远期供水方案为依据，采用平面二维数学模型定量分析了天角潭水利枢纽工程两种供水方案对河口水环境的影响。 1 工程概况 1.1 流域概况 北门江流域位于海南岛西北部，东经109°18′～109°39′和北纬19°24′～19°48′之间，流域集雨面积为621 km2，干流河长67 km，干流平均坡降0.21%。北门江从东南向西北流，地形总的趋势东南高，西北低，由东南内陆向西北沿海逐渐降低。流域上游属丘陵区，下游属阶地和海滨平原区，地势平坦。北门江有东西两源，其中西支沙河为干流，发源于儋州市的纱帽岭，向北流经西培农场至沙河村从左岸汇入牙拉河，流域面积84.7 km2；东支为南茶河，发源于莲花岭，向北流经南洋农场、新村至牙拉河苗圃从右岸汇入牙拉河，流域面积110 km2，河长19.4 km，干流平均坡降0.43%。东西两支于儋州市沙河村汇合后称“牙拉河”，再北流与巴黎河汇合后始称“北门江”，经长坡、东成、中和等乡镇，于新英湾出海。沿途小支流较多，较大的分叉支流主要位于西庆二十五队、西联光明队以及天角潭水陂附近，北门江流域水系图见图1。 图1 北门江流域水系 1.2 工程概况 天角潭水利枢纽工程位于北门江干流天角潭水陂上游500 m处，是国家172项重大水利工程项目，也是《海南省北门江流域规划(修编)》《海南水网规划》中确定的近期建设的大(2)型水利枢纽工程。工程任务以工业供水、农业灌溉为主，兼顾发电等综合利用。水库正常蓄水位58.0 m，兴利库容为1.54亿m3，总库容1.94亿m3，具有多年调节性能。工程近期(2030年)多年平均供水量为7 974万m3，全部用于农业灌溉；远期(2040年)多年平均供水量为13 557万m3，其中工业供水量为6 138万m3，农业灌溉供水量为7 419万m3。工程主要建设内容包括1座主坝、4座副坝、引水隧洞、渠首电站、坝后电站、天角潭总干渠、过鱼设施、增殖放流站等。其中，主坝最大坝高52 m，副坝最大坝高32 m，渠首和坝后电站总装机容量5 000 kW。 2 数学模型 2.1 水动力方程 水动力控制方程为经Navier-Stokes方程沿深积分的二维浅水方程组 (1) (2) (3) 式中，h=η+d为总水深；η为自由面高程；d为静水深；x和y分别表示横轴和纵轴坐标；t为时间；g为重力加速度；U和V分别为沿x和y方向的垂向平均流速；ρ=ρ(T，S)为与温度、盐度相关的水体密度；ρ0为参考密度；Tij为应力项，包括粘性应力、紊流应力和对流等，与流速梯度相关。 (4) 式中，cf为拖曳力系数，cf=g/C2。 水平涡粘性系数采用Samagorinsky亚网格尺度模型求解，涡粘系数 (5) 式中，Cs为可调系数；Sij为变形速率，与速度梯度相关。 2.2 盐度扩散方程 盐度扩散模拟采用对流扩散方程 (6) 式中，S为盐度；Dh为水平扩散系数，Dh=A/σT，σT为Prandtl数。 表1 计算工况 m3/s 3 计算结果分析 3.1 研究范围 根据北门江河口洋浦港潮位观测资料和北门江实测地形数据，并结合历史水质观测资料，萍塘村断面是感潮河段的主要分界断面。本文以萍塘村断面上游约2 km处为盐度扩散模拟计算的上游边界，下游至新英湾至开源大道外7.5 km的部分海域(洋浦鼻外海域约3.5 km)，图