普兰店电厂温排水数值模拟研究.docVIP

PAGE  PAGE 10 普兰店电厂温排水数值模拟研究 　　摘要：为解普兰店电厂温排水对海水环境的影响，采用水动力及对流-扩散方程的二维潮流及温度扩散数学模型进行了潮流及温度场数值模拟，分析了温排水排放对周边海水的影响，统计了不同温升的包络面积。结果表明：由温排放导致的温度扩散主要集中在排水口附近区域，冬季温升超过0.5℃的面积约为0.0572 km2，温升大于1℃的范围仅0.0001 km2；夏季温升大于1℃的范围均位于水道内，面积为0.0551 km2，不会影响到水道外侧的普兰店湾海域；温升大于0.5℃的海水面积约为0.0118 km2，影响范围为水道内及水道外侧约0.0045 km2的局部海域。 　　关键词：温排水；数学模型；水动力模拟；潮流；普兰店 　　中图分类号：TV143 文献标识码：A 文章编号：1674-9944（2017）2-0083-05 　　1 引言 　　半封闭海湾湾内水体与外海水体交换缓慢，造成海洋生态系统脆弱[1]。近年来，随着我国沿海经济的快速发展，许多半封闭海湾建造了电厂，采用直流供水系统的电厂常年源源不断的将高于环境温度5～12℃的温排水排入受纳海域，会对海水环境产生热影响乃至热污染。 　　目前，半封闭海湾的热影响问题受到学者的广泛关注。周巧菊采用POM模型对大亚湾核电站的温排水进行了模拟研究[2]；林军采用ECOMSED模型对象山港海域温排水的时空输移路径及范围进行了研究[3]；赵懿?B采用Delft3D模型对铁山湾电厂的温度场进行了数值模拟[1]；陈春亮等为了解决温排水排入湛江湾后的随潮输移扩散，应用ECOMSED模型计算温排水的时空输移路径及温升范围[4]；姚姗姗等采用MIKE21-HD模型对舟山海域冷排水的温度扩散进行了数值研究[5]。数学模型作为主要的模拟手段，已得到了广泛的应用。笔者采用MIKE21-HD模型，对普兰店电厂温排水的输移路径、温升范围进行研究，从温排水扩散角度分析项目建设是否可行。 　　2 自然条件及计算工况 　　普兰店海域潮流类型属规则半日潮流性质，日潮不等现象明显，即高、低潮不等较为明显。普兰店湾内地貌主要为海积地貌[6]，海岸地貌都是在海水动力下逐步形成的。海积地貌主要表现为海滩、潮上滩、沿海堤、沙嘴和?_积平原等。大连海湾工业园区主要是坐落在普兰店湾内的海滩上，既是潮间带上堆积的细小泥质物资。工程海区形势图如图1所示。各验潮站潮特征值引自区项目填海可研海洋环评海域使用论证阶段水文测验分析报告》，具体结果列于表1中。海流实测共布设了10个测站，1#测站采用H1潮位，2#、3#、4#测站采用H2潮位，5#～10#测站采用H3潮位，将各个测站的垂线平均流速以落潮为正、涨潮为负绘制潮位及垂线平均流速流向过程线图，具体见图2和图3。根据实测结果，海域潮流类型属规则半日潮流性质。垂线平均潮流的可能最大流速，以6#测站最大，为90.7 cm/s， 10#测站最小，为7.7 cm/s。可能最大流速的流向多为落潮方向。全部8个测站均表现为往复流。余流的变化主要受风场、径流以及地形的支配，从垂线平均的余流结果来看，最大值出现在小潮期间5#测站，达12.6 cm/s，方向为116°，其次是小潮期间7#测站，达11.5 cm/s，方向为71°，其余测站不超过10 cm/s。从各层余流结果看，最大值出现在小潮期间9#测站表层，达15.8 cm/s，方向为305°，其余测站不超过10 cm/s。 　　项目所在地的多年平均气温为9.7 ℃，多年平均最高气温为15.1℃，平均最低气温为4.9℃，极端最高气温为36.1℃（1999年7月29日），极端最低气温-24.8℃（2001年1月15日）。 　　为保证局部流场计算符合潮流场的整体物理特性，便于给出开边界条件，采用网格嵌套方法进行计算，第一模型为普兰店湾，取为56 km×45 km的南北向矩形区域，本计算域基本上将工程区及可能受到影响的区域都包括在内，计算网格尺度采用90 m×90 m的固定网格。第二模型为普兰店内湾海域，计算范围取56 km×45 km，网格尺度为30 m×30 m。第三模型为工程附近的局部区域，取为28 km×18 km的南北向矩形区域，计算网格尺度采用10 m×10 m的固定网格。在第一模型计算后，第二第三模型的开边界均由上一层模型的计算结果提供。 　　模型计算的地形条件分别采用相应的海图中的数据，小范围计算采用工程附近的实测地形。 　　本项目排水口设置于普兰店湾纵深末端海域与附近支流交汇的水道内，水道相对狭窄，排水口距岸边约430 m。根据本项目工艺要求，夏季温排水的排水量为976 m3/h，温升约为6 ℃。冬季温排水的排水量为335 m3/h，温升约为9.6 ℃。排水口位于一