电厂温排水数值模拟
|
摘要
发电厂的废热通过温排水的形式排入自然水体所造成的影响是近年来人们十分关注的环境问题，同时也是科技工作者研究的一个重要方向。为了较好的评估电厂温排水的影响，本文建立了一个用于预测评估电厂温排水对受纳水体水温分布影响情况的数学模型—平面二维有限单元水动力热输运数学模型，对温排水进入水体后的水流运动及温度扩散规律进行了研究。
主要内容有：
(1) 对温排水的研究方法和前人主要研究成果进行了较为全面的总结，指出运用数学模型研究温排水运动扩散规律的必要性及可行性。
(2) 采用一种简单便捷的三角形网格生成方法生成有限单元三角形网格，运用距离最小二乘法插值补足节点高程值，良好的拟合计算区域边界并快速获得水下地形。
(3) 基于无结构三角形网格，建立了沿水深平均二维有限单元水动力数学模型，对河段的水流流态进行了研究。
(4) 在水动力计算的基础上，建立了六节点三角形单元的有限节点法温度扩散数学模型，研究水流运动对温度场的影响。
(5) 作为案例，采用本文所建的数学模型分别在大潮、小潮情况下对华能电厂温排水的水流和温度进行了计算，对大、小潮的涨潮和落潮的流场及温度场进行了分析。
(6) 采用已建成的电厂一期工程运行中水温同步实测数据，将计算所得到的流速场、潮位和温度场与其实测值进行了对比经验证，对拟建的二期工程的温排水温度影响范围及取水口温升进行了预测计算，验证结果比较理想，说明本模型可以应用于相似的工程和研究，预测结果可以为电厂设计和环境评价提供科学依据。
The problem of waste heat discharge into natural bodies of water is one of the environmental problems that has received considerable attention in recent years, which also is an major research content for the environment scholar. A 2-D numerical model, which using the finite element methods, is established in order to evaluate how the thermal discharge affects the water body. The main conclusions include as follows:
1. The advances in the research of thermal discharge are reviewed overall. The necessary and feasibility of the application of mathematical models in the research of thermal discharge into water body are pointed out.
2. A simple grid generation method is applied to gain the finite element network, and least square method was adopted for the elevation interpolation.
3. A 2-Dimension mathematical model is established on unstructured triangular grids, and used to simulate hydrodynamic field.
4. A mathematical model is estiblished to simulate the temperature transport and diffusion after the hydrodynamic computation. A finite nodes method based on quadratic triangular elements is adopted.
5. As an application case, the model is applied to the cooling water project of Hua'neng Power Plant. The computation is conducted at situations of spring tide and neap tide. An analysis is made according to the computation results.
6. Compared the computation results to the real mesured data of the first-stage project of plant, the verification of flow velocities, water depths and temperature rises is complished. As refer to the second-stage project, numerical prediction about temperature influence area and intake temperature rise of thermal effluent was made. The results of verification are merged well and the conclusions of the prediction provide the scientifical references for project design of power station and environmental assessment.
主要内容有：
(1) 对温排水的研究方法和前人主要研究成果进行了较为全面的总结，指出运用数学模型研究温排水运动扩散规律的必要性及可行性。
(2) 采用一种简单便捷的三角形网格生成方法生成有限单元三角形网格，运用距离最小二乘法插值补足节点高程值，良好的拟合计算区域边界并快速获得水下地形。
(3) 基于无结构三角形网格，建立了沿水深平均二维有限单元水动力数学模型，对河段的水流流态进行了研究。
(4) 在水动力计算的基础上，建立了六节点三角形单元的有限节点法温度扩散数学模型，研究水流运动对温度场的影响。
(5) 作为案例，采用本文所建的数学模型分别在大潮、小潮情况下对华能电厂温排水的水流和温度进行了计算，对大、小潮的涨潮和落潮的流场及温度场进行了分析。
(6) 采用已建成的电厂一期工程运行中水温同步实测数据，将计算所得到的流速场、潮位和温度场与其实测值进行了对比经验证，对拟建的二期工程的温排水温度影响范围及取水口温升进行了预测计算，验证结果比较理想，说明本模型可以应用于相似的工程和研究，预测结果可以为电厂设计和环境评价提供科学依据。
The problem of waste heat discharge into natural bodies of water is one of the environmental problems that has received considerable attention in recent years, which also is an major research content for the environment scholar. A 2-D numerical model, which using the finite element methods, is established in order to evaluate how the thermal discharge affects the water body. The main conclusions include as follows:
1. The advances in the research of thermal discharge are reviewed overall. The necessary and feasibility of the application of mathematical models in the research of thermal discharge into water body are pointed out.
2. A simple grid generation method is applied to gain the finite element network, and least square method was adopted for the elevation interpolation.
3. A 2-Dimension mathematical model is established on unstructured triangular grids, and used to simulate hydrodynamic field.
4. A mathematical model is estiblished to simulate the temperature transport and diffusion after the hydrodynamic computation. A finite nodes method based on quadratic triangular elements is adopted.
5. As an application case, the model is applied to the cooling water project of Hua'neng Power Plant. The computation is conducted at situations of spring tide and neap tide. An analysis is made according to the computation results.
6. Compared the computation results to the real mesured data of the first-stage project of plant, the verification of flow velocities, water depths and temperature rises is complished. As refer to the second-stage project, numerical prediction about temperature influence area and intake temperature rise of thermal effluent was made. The results of verification are merged well and the conclusions of the prediction provide the scientifical references for project design of power station and environmental assessment.
引文
[1] 周鑫玉．．电厂废水的环境效应及其处理技术[J]．环境科学与技术，1997．
[2] 张育兰．火力发电厂与水境[J]．湖南电力，2000．
[3] WRIGHT S. J. Mean behavior of buoyant jets in a cross flow[J]. Proc. ASCE, J.Hydraul. Div., 1977, 103(HYS): 499～513.
[4] LEE J. H., JIRKA G. H. Multiport diffiser as line source of momentum in shallow water[J].Water Res Res, 1980, 16(4): 695～705.
[5] Engelman, M. S., Strang" G., and Bathe, K. J. (1981). "The Application of Quasi-Newton Methodsin Fluid Mechanics."Intemational Journal for Numerical Methods in Engineering, 17(5), 707～718.
[6] 黄明游，刘播等．数值计算方法[M]：科学出版社，2005．
[7] Klaus-Jurgen Bathe. Finite Element Procedures[M]:Prentice Hall; 2r.e. edition (1995)
[8] FEPG有限元自动生成程序系统参考手册[M]：北京飞箭软件公司，2005．
[9] 金忠青．N-S方程的数值解和紊流模型[M]：河海大学出版社，1989．
[10] 周学漪．计算水力学[M]：清华大学出版社，1995．
[11] 谭维炎．计算浅水动力学[M]：清华大学出版社，1998．
[12] 汪德爟．计算水力学[M]：河海大学出版社，1989．
[13] Anderson J. D.Computationai fluid dynamics: the basics with applications/John[M]: McGraw-Hill, Inc., 1995.
[14] Kjjellgren P. A semi-implicit fractional step finite element method for viscous incompressible flows[J].computional mechanics, 1997, 20: 541～550.
[15] 施祖蓉，施麟宝．萧山电厂温排水二维热污染数值模型[J]．浙江水利科技，1991(01)．
[16] 帕坦卡．传热与流动的数值计算[M]．北京：科学出版社，1984．
[17] 杨芳丽，谢作涛，张小峰，等．非正交曲线坐标系平面二维电厂温排水模拟[J]．水利水运工程学报，2005(02)．
[18] 陶文铨．数值传热学[M]．西安：西安交通大学出版社，2001
[19] Mesh Generation Version 1.0[M]: Edinburgh Parallel Computing Centre (The University of Edinburgh), 1996.
[20] 关振群，宋超，顾元宪等．有限元网格生成方法研究的新进展[J]．计算机辅助设计与图形学学报，2003，15(1)．
[21] 刘智，林秉南，何少苓．三角形网格在二维不恒定流计算中的应用[J]．水利学报，1987(09)．
[22] 罗小峰，长江口水流盐度数值模拟，南京水利科学研究院博士学位论文，2003，39-39
[23] 潘存鸿．三角形网格在潮汐河流水质计算中的应用[J]．热带海洋，1996(04)．
[24] 唐洪武，徐久荣，等．粒子图象测速技及其在垂直进水口漩涡流场中的应用[J]．水动力学研究进展，1999．14(1)：128～134．
[25] 何生存，郭三刚，等．离散点数据插值方法及等值线绘制技术的应用[J]．青海科技，2006(04)．
[26] Buell, W. R., and Bush, B. A.. "Mesh generation—a survey." Transactions of the American Society of Mechanical Engineers, Journal of Engineering for Industry, set. B.,95(1), 332-338.
[27] 华祖林，褚克坚．基于三角形网格的潮汐水域水流水质的一种计算模式[J]．河海大学学报 (自然科学版)，2001(04)．
[28] 张细兵，吴腾，卢金友，等．潮汐河口平面二维非恒定流数值模拟[J]．水动力学研究与进展A辑，2003(04)．
[29] 李光炽，饶光辉．分叉型海湾温水排放数学模型[J]．水利水运工程学报，2005(03)：20～25．
[30] Harleman D. R. F., Hall L. C. Thermal Diffusion of Condenser Water in A River During Steady and Unsteady Flows with Application to the T.V.A. Browns Ferry Nuclear Power Plant[J].Hydrodynamics Laboratory Report, 1968(3): 98～115.
[31] McGuirk J. J., Rodi W.'A Depth averaged Mathematical Model for Near Field of Side Discharges into Open channel Flow'[J].Fluid Mech., 1978( 86): 761～781.
[32] Casulli V. Numerical simulation of shallow water flow[J].Computational Methods in Surface Hydrology, 1990: 13-22.
[33] Hamriek JV.Analysis of water temperatures in Conowrngo Pond as influenced by the Peach Bottom atomic power plant thermal discharge[J].Environmental Science&Policy, 2000, 3.
[34] J. Jiang, D.B. Fissel and A. Taylor. Burrard Generating Station Cooling Water Recirculation Study[J].ASL Technical Report, 2001.
[35] J. Jiang, D. B. Fissel and D. D.Lemon. Modeling Cooling Water Discharges From the Burrard Generating Station[M]: Presented at oceans 2002, 2002.
[36] 吴江航，韩庆书，张继春，等．潮汐河道中热电厂冷却水系统水力热力数值模拟[J]计算物理，1987(01)．
[37] 吴江航，陈凯麒，韩庆书．核电站冷却水远区热、核污染数值计算的一种新方法[J]．水利学报，1986(10)．
[38] 李燕初，蔡文理．沿海港口电厂温捧水、废水远区影响数值模拟[J]．台湾海峡，1988(9)：49～61．
[39] 吴时强．剖开算子法解具有自由表面的温排水平面紊流流场．：南京水利科学研究所；1989．
[40] 黄平．哑铃湾电厂温排水扩散预测[J]．海洋环境科学，1992(04)．
[41] 黄平．汕头港水域温排水热扩散的三维数值模拟[J]．海洋环境科学，1996(01)．
[42] 华祖林．电厂温排放对感潮河段环境水体影响预测研究[J]．电力环境保护，1995(04)．
[43] 华祖林．贴体边界下潮汐河口热(核)电厂温排放数值计算[J]．电力环境保护，1997(02)．
[44] 董耀华．河道温度及浓度场的平面二维数值模拟[J]．长江科学院院报，1997．
[45] 王丽霞，孙英兰，郑连远．三维热扩散预测模型[J]．中国海洋大学学报(自然科学版)，1998(01)．
[46] 王丽霞，孙英兰，田晖．热扩散预测方法研究概况Ⅰ．影响海洋水温的因素[J]．海洋科学，1997(05)．
[47] 王丽霞，孙英兰，田晖．热扩散预测方法研究概况Ⅱ．热扩散的研究现状[J]．海洋科学，1997(06)．
[48] 韩康，张存智，张砚峰，等．三亚电厂温排水数值模拟[J]．海洋环境科学，1998(02)．
[49] 徐啸，匡翠萍，顾杰．漳州后石电厂温排水数学模型[J]．台湾海峡，1998(02)．
[50] 罗斌，朱京兴，唐辉．有限节点法在电厂温排水计算中的应用[J]．电力勘测，1999(01)．
[51] 江洧．惠州LNG电厂冷却水工程数值模拟研究[J]．广东水利水电，2001(04)．
[52] 江洧，林佑金，陆耀辉．惠州LNG电厂循环冷却水工程模型试验研究[J]．人民珠江，2001(04)．
[53] 程杭平、韩曾萃．热污染的一、二维耦合模型及其应有[J]．水动力研究与进展 2002(06)： 648～655．
[54] 李振海，祝秋梅．填海工程影响下的电厂冷却水工程布置数值模拟研究[J]．电力环境保护，2003(03)．
[55] 郝瑞霞，韩新生．潮汐水域电厂温排水的水流和热传输准三维数值模拟[J]．水利学报，2004(08)．
[56] 郝瑞霞，齐伟，李海香，等．潮汐水域流速场和温度场的数值模拟研究[J]．太原理工大学学报，2005(03)．
[2] 张育兰．火力发电厂与水境[J]．湖南电力，2000．
[3] WRIGHT S. J. Mean behavior of buoyant jets in a cross flow[J]. Proc. ASCE, J.Hydraul. Div., 1977, 103(HYS): 499～513.
[4] LEE J. H., JIRKA G. H. Multiport diffiser as line source of momentum in shallow water[J].Water Res Res, 1980, 16(4): 695～705.
[5] Engelman, M. S., Strang" G., and Bathe, K. J. (1981). "The Application of Quasi-Newton Methodsin Fluid Mechanics."Intemational Journal for Numerical Methods in Engineering, 17(5), 707～718.
[6] 黄明游，刘播等．数值计算方法[M]：科学出版社，2005．
[7] Klaus-Jurgen Bathe. Finite Element Procedures[M]:Prentice Hall; 2r.e. edition (1995)
[8] FEPG有限元自动生成程序系统参考手册[M]：北京飞箭软件公司，2005．
[9] 金忠青．N-S方程的数值解和紊流模型[M]：河海大学出版社，1989．
[10] 周学漪．计算水力学[M]：清华大学出版社，1995．
[11] 谭维炎．计算浅水动力学[M]：清华大学出版社，1998．
[12] 汪德爟．计算水力学[M]：河海大学出版社，1989．
[13] Anderson J. D.Computationai fluid dynamics: the basics with applications/John[M]: McGraw-Hill, Inc., 1995.
[14] Kjjellgren P. A semi-implicit fractional step finite element method for viscous incompressible flows[J].computional mechanics, 1997, 20: 541～550.
[15] 施祖蓉，施麟宝．萧山电厂温排水二维热污染数值模型[J]．浙江水利科技，1991(01)．
[16] 帕坦卡．传热与流动的数值计算[M]．北京：科学出版社，1984．
[17] 杨芳丽，谢作涛，张小峰，等．非正交曲线坐标系平面二维电厂温排水模拟[J]．水利水运工程学报，2005(02)．
[18] 陶文铨．数值传热学[M]．西安：西安交通大学出版社，2001
[19] Mesh Generation Version 1.0[M]: Edinburgh Parallel Computing Centre (The University of Edinburgh), 1996.
[20] 关振群，宋超，顾元宪等．有限元网格生成方法研究的新进展[J]．计算机辅助设计与图形学学报，2003，15(1)．
[21] 刘智，林秉南，何少苓．三角形网格在二维不恒定流计算中的应用[J]．水利学报，1987(09)．
[22] 罗小峰，长江口水流盐度数值模拟，南京水利科学研究院博士学位论文，2003，39-39
[23] 潘存鸿．三角形网格在潮汐河流水质计算中的应用[J]．热带海洋，1996(04)．
[24] 唐洪武，徐久荣，等．粒子图象测速技及其在垂直进水口漩涡流场中的应用[J]．水动力学研究进展，1999．14(1)：128～134．
[25] 何生存，郭三刚，等．离散点数据插值方法及等值线绘制技术的应用[J]．青海科技，2006(04)．
[26] Buell, W. R., and Bush, B. A.. "Mesh generation—a survey." Transactions of the American Society of Mechanical Engineers, Journal of Engineering for Industry, set. B.,95(1), 332-338.
[27] 华祖林，褚克坚．基于三角形网格的潮汐水域水流水质的一种计算模式[J]．河海大学学报 (自然科学版)，2001(04)．
[28] 张细兵，吴腾，卢金友，等．潮汐河口平面二维非恒定流数值模拟[J]．水动力学研究与进展A辑，2003(04)．
[29] 李光炽，饶光辉．分叉型海湾温水排放数学模型[J]．水利水运工程学报，2005(03)：20～25．
[30] Harleman D. R. F., Hall L. C. Thermal Diffusion of Condenser Water in A River During Steady and Unsteady Flows with Application to the T.V.A. Browns Ferry Nuclear Power Plant[J].Hydrodynamics Laboratory Report, 1968(3): 98～115.
[31] McGuirk J. J., Rodi W.'A Depth averaged Mathematical Model for Near Field of Side Discharges into Open channel Flow'[J].Fluid Mech., 1978( 86): 761～781.
[32] Casulli V. Numerical simulation of shallow water flow[J].Computational Methods in Surface Hydrology, 1990: 13-22.
[33] Hamriek JV.Analysis of water temperatures in Conowrngo Pond as influenced by the Peach Bottom atomic power plant thermal discharge[J].Environmental Science&Policy, 2000, 3.
[34] J. Jiang, D.B. Fissel and A. Taylor. Burrard Generating Station Cooling Water Recirculation Study[J].ASL Technical Report, 2001.
[35] J. Jiang, D. B. Fissel and D. D.Lemon. Modeling Cooling Water Discharges From the Burrard Generating Station[M]: Presented at oceans 2002, 2002.
[36] 吴江航，韩庆书，张继春，等．潮汐河道中热电厂冷却水系统水力热力数值模拟[J]计算物理，1987(01)．
[37] 吴江航，陈凯麒，韩庆书．核电站冷却水远区热、核污染数值计算的一种新方法[J]．水利学报，1986(10)．
[38] 李燕初，蔡文理．沿海港口电厂温捧水、废水远区影响数值模拟[J]．台湾海峡，1988(9)：49～61．
[39] 吴时强．剖开算子法解具有自由表面的温排水平面紊流流场．：南京水利科学研究所；1989．
[40] 黄平．哑铃湾电厂温排水扩散预测[J]．海洋环境科学，1992(04)．
[41] 黄平．汕头港水域温排水热扩散的三维数值模拟[J]．海洋环境科学，1996(01)．
[42] 华祖林．电厂温排放对感潮河段环境水体影响预测研究[J]．电力环境保护，1995(04)．
[43] 华祖林．贴体边界下潮汐河口热(核)电厂温排放数值计算[J]．电力环境保护，1997(02)．
[44] 董耀华．河道温度及浓度场的平面二维数值模拟[J]．长江科学院院报，1997．
[45] 王丽霞，孙英兰，郑连远．三维热扩散预测模型[J]．中国海洋大学学报(自然科学版)，1998(01)．
[46] 王丽霞，孙英兰，田晖．热扩散预测方法研究概况Ⅰ．影响海洋水温的因素[J]．海洋科学，1997(05)．
[47] 王丽霞，孙英兰，田晖．热扩散预测方法研究概况Ⅱ．热扩散的研究现状[J]．海洋科学，1997(06)．
[48] 韩康，张存智，张砚峰，等．三亚电厂温排水数值模拟[J]．海洋环境科学，1998(02)．
[49] 徐啸，匡翠萍，顾杰．漳州后石电厂温排水数学模型[J]．台湾海峡，1998(02)．
[50] 罗斌，朱京兴，唐辉．有限节点法在电厂温排水计算中的应用[J]．电力勘测，1999(01)．
[51] 江洧．惠州LNG电厂冷却水工程数值模拟研究[J]．广东水利水电，2001(04)．
[52] 江洧，林佑金，陆耀辉．惠州LNG电厂循环冷却水工程模型试验研究[J]．人民珠江，2001(04)．
[53] 程杭平、韩曾萃．热污染的一、二维耦合模型及其应有[J]．水动力研究与进展 2002(06)： 648～655．
[54] 李振海，祝秋梅．填海工程影响下的电厂冷却水工程布置数值模拟研究[J]．电力环境保护，2003(03)．
[55] 郝瑞霞，韩新生．潮汐水域电厂温排水的水流和热传输准三维数值模拟[J]．水利学报，2004(08)．
[56] 郝瑞霞，齐伟，李海香，等．潮汐水域流速场和温度场的数值模拟研究[J]．太原理工大学学报，2005(03)．