三维水质数学模型修正及其在浙东近海海域水质模拟中的应用
摘要
浙东近海海域大多属于陆架浅海,潮汐、潮流等海洋动力强劲,而污染物的迁移、扩散受水动力的控制。本文在前人研究成果的基础上,开展浙东近海水域实测海洋水文、水质资料的分析,根据海陆水文条件,为数学模型提供边界条件。浙江沿海水域水下地形宽浅,因此底摩阻因子对水动力变化尤为重要,从而影响到污染物的传播和扩散,所以对有限体积海岸海洋模型(FVCOM)模型的源代码中的底摩阻项进行了修改,利用浙东近海实测水下地形、潮位、潮流资料对模型进行率定,并根据典型潮过程对模型潮流验证,选用具有污染物浓度的测验资料进行浙东近海三维水质模型验证。最后,作者结合浙东近海的涉海工程,开展工程对局部水流、水质影响的模拟计算。
本研究根据建立的包含浙东近海各区域的三维潮流、水质数值模型,结合SMS软件的前、后处理模块,用于生成高质量、非结构的三角形网络,并将FVCOM模式结果转移至SMS平台进行流场分析以及结果后处理,结合google earth工具,构建了潮流场和浓度场的三维动态图形显示系统。根据区域的潮流特征,开发了底部阻力参数计算模块,使之更适合目标区域的特点。
本项研究通过对所构建模型的使用,证明该模型可以较好模拟我省沿海的潮汐及三维潮流特点,为浙江省海洋工程建设的前期研究提供了可靠潮流分析模型,提高了模拟的精度及深度。同时应用FVCOM建立的浙东近海海域三维水质模型,对浙江海域三维水质进行了计算,经过进一步深入研究,可以使模型更加完善,也为将来进一步发展打下了坚实的基础。
Offshore areas of Zhejiang Province are mostly shelf seas, so tide and flow played an important role in the controlling of ocean dynamic. Based on the previous achievements, this study analyzed the data of hydrology and water quality of Zhejiang offshore areas, and determined the boundary conditions for the mathematical model based on hydrological conditions. Then, the FVCOM(Finite Volume Coastal Ocean Model) model was developed and rated using the data of underwater topography, tide level and tidal current. And the three-dimensional water quality model was verified according to the monitoring data of water temperature and pollutants concentration in offshore areas of Zhejiang Province. Finally, the effects of some marine projects in Zhejiang offshore areas onto local flow and water quality were simulated.
This study was conducted based on the built three-dimensional tide and water quality numerical models for coastal regions of Zhejiang, and the pre-and post-processing module of the SMS software. The high-quality, non-triangular structure grids were generated, and the results of FVOCM were transferred to the SMS platform for flow analysis and post-processing. Then, combined with google earth tool, three-dimensional graphics display system for tidal current field and concentration field was built. Additionally, according to the characteristics of the regional tide flow, the module for bottom resistance parameter calculation was developed to adapt to the target area.
the established modified model could well simulate the coastal tide and three-dimensional flow of Zhejiang Province. The model also provided reliable tidal flow analysis for preliminary research of the local ocean engineers, and improved the accuracy and depth of the simulation. Besides, the three-dimensional water quality model modefied from FVCOM, was used for the calculation of three-dimensional water quality of Zhejiang Province. After further investigation, the model could be more perfect and would set a solid foundation for the further development.
本研究根据建立的包含浙东近海各区域的三维潮流、水质数值模型,结合SMS软件的前、后处理模块,用于生成高质量、非结构的三角形网络,并将FVCOM模式结果转移至SMS平台进行流场分析以及结果后处理,结合google earth工具,构建了潮流场和浓度场的三维动态图形显示系统。根据区域的潮流特征,开发了底部阻力参数计算模块,使之更适合目标区域的特点。
本项研究通过对所构建模型的使用,证明该模型可以较好模拟我省沿海的潮汐及三维潮流特点,为浙江省海洋工程建设的前期研究提供了可靠潮流分析模型,提高了模拟的精度及深度。同时应用FVCOM建立的浙东近海海域三维水质模型,对浙江海域三维水质进行了计算,经过进一步深入研究,可以使模型更加完善,也为将来进一步发展打下了坚实的基础。
Offshore areas of Zhejiang Province are mostly shelf seas, so tide and flow played an important role in the controlling of ocean dynamic. Based on the previous achievements, this study analyzed the data of hydrology and water quality of Zhejiang offshore areas, and determined the boundary conditions for the mathematical model based on hydrological conditions. Then, the FVCOM(Finite Volume Coastal Ocean Model) model was developed and rated using the data of underwater topography, tide level and tidal current. And the three-dimensional water quality model was verified according to the monitoring data of water temperature and pollutants concentration in offshore areas of Zhejiang Province. Finally, the effects of some marine projects in Zhejiang offshore areas onto local flow and water quality were simulated.
This study was conducted based on the built three-dimensional tide and water quality numerical models for coastal regions of Zhejiang, and the pre-and post-processing module of the SMS software. The high-quality, non-triangular structure grids were generated, and the results of FVOCM were transferred to the SMS platform for flow analysis and post-processing. Then, combined with google earth tool, three-dimensional graphics display system for tidal current field and concentration field was built. Additionally, according to the characteristics of the regional tide flow, the module for bottom resistance parameter calculation was developed to adapt to the target area.
the established modified model could well simulate the coastal tide and three-dimensional flow of Zhejiang Province. The model also provided reliable tidal flow analysis for preliminary research of the local ocean engineers, and improved the accuracy and depth of the simulation. Besides, the three-dimensional water quality model modefied from FVCOM, was used for the calculation of three-dimensional water quality of Zhejiang Province. After further investigation, the model could be more perfect and would set a solid foundation for the further development.
引文
[1]中国海湾志编纂委员会,中国海湾志[M].海洋出版社,1994.
[2]邹志利,海岸动力学[M].人民交通出版社,2009.9.
[3]陈吉作主编,中国海岸侵蚀概要[M].海洋出版社,2010.7.
[4]朱学明,中国近海潮汐潮流的数值模拟与研究[D],中国海洋大学硕士论文,2009.
[5]申霞,基于POM的近海三维水质模型研究及其应用[D],河海大学硕士论文,2006.
[6]章卫胜,中国近海潮波运动数值模拟[D],河海大学硕士论文,2010.
[7]吴伦宇,基于FVCOM的浪、流、泥沙模型耦合及其应用[D],中国海洋大学博士论文,2010.
[8]曹德明,方国洪.杭州湾潮汐潮流的数值计算[J].海洋和湖沼,1986,17(2)：93-101.
[9]曹德明,方国洪.杭州湾和钱塘江潮波的联合数值模型[J].海洋学报,1988,10(5)：521-530.
[10]许卫忆,苏纪兰.杭州湾二维潮波计算及底质分布的动力成因[J].海洋与湖沼,1986(6)：52-56.
[11]曹欣中,唐龙妹张月秀.宁波,舟山内海域实测海流分析及潮流场的数值模拟[J].东海海洋,1996,14(2)：1-9.
[12]董礼先,苏纪兰.象山港潮波响应和变形研究.Ⅱ.象山港潮波数值研究[J].海洋学报,1999,21(2)：1-8.
[13]李孟国,王正林.瓯江口潮流数值模拟[J].长江科学院院报,2002,19(2)：19-22.
[14]李孟国,时钟,范文静.瓯江口外洞头岛北岱后涂围垦工程潮流数值模拟研究[J].海洋通报,2005,24(1)：1-7.
[15]陈晒睿,朱建荣,吴辉等.无结构网格二维河口海岸水动力数值模式的建立及其应用[J].海洋学报,2010,32(2)：31-39
[16]孙文心.超浅海风暴潮的数值模拟.海洋学报[J],1979,1(2)：193-211。
[17]吴培木,许永,李燕初等.台湾海峡台风暴潮非线性数值计算.海洋学报[J],1981,3(1)： 27～43.
[18]陈长胜,秦曾灏.江浙沿海台风暴潮的动力分析.海洋学报[J],1985,7(3)：265-275.
[19]C.S.Chen, Z.H.Qin. Numerical simulation of typhoon surges along the east coast of Zhejiang and Jiangsu Provinces. Advances in Atmosphere Science[J],1985,2(1):8-9.
[20]X.N Wang,F.J Yu, Q.J Yin. Research of Application of Numerical Models of Typhoon Surges in China Seas. The Special Issue of MAUSAM,1997:595-608
[21]赵有皓,张君伦.风暴潮实用数值预报的研发.海洋工程,2001,19(3)：102-107.
[22]Hassan S.Mashriqui,G.Paul Kemp, and Ivor van Heerden. Experimental storm surge simulations for hurricane Katrina.In:Y. J. Xu&V.P.Singh, Coastal Environment and Water Quality. LLC, Highlands Ranch,Co 80163-0026,USA:Water resources Publications,2006:481-490.
[23]Jian Shen,Keqi Zhang.Improved prediction of storm surge inundation with a high-resolution unstructured grid model.Journal of Coastal Research,2006,Vol.22,No.6:1309-1319.
[24]张涤明,陈良圭,乔林.近海水域海洋流体动力学方程的数值模拟[J].力学进展,1988,18(3)：322-331.
[25]李身铎,顾思美.杭州湾潮波三维数值模拟[J].海洋与湖沼,1993,24(1)：7-15.
[26]韩国其,汪德灌,许协庆.潮汐河口三维水流数值模拟[J].水利学报,1989,12(5)：42-60.
[27]万振文,乔方利,袁业立.渤,黄,东海三维潮波运动数值模拟[J].海洋与湖沼,1998,29(6)：611-616.
[28]W.B.Guan,L.A..Wong, D.F. Xu,2001. Modeling nitrogen and phosphorus cycles and dissolved oxygen in the Zhujiang River estuary. Acta Oceanologica Sinica,20(4):493-514.
[29]史峰岩,朱首贤.杭州湾,长江口余流及其物质输运作用的模拟研究[J].海洋学报,2000,22(5)：1一12.
[30]华祖林,顾莉,查玉含,刘晓东.基于COHERENS模型的污染物质输运数值模拟[J].环境科学与技术,2009,32(4)：14-18.
[31]周大成.椒江河口大潮潮流特性的数值分析[J].浙江水利水电专科学校学报,2001,13(3)：13-26.
[32]杨陇慧,朱建荣.长江口杭州湾及邻近海区潮汐潮流场三维数值模拟[J].华东师范大学学报：自然科学版.2001(3)：74-84.
[33]陈倩,黄大吉,章本照.浙东近海潮汐潮流的数值模拟[J].海洋学报,2003,25(5)：9-20.
[34]张威.浙东近海海域潮波数值模拟[D],浙江大学硕士论文,2007.
[35]C.S. Chen,H.D. Liu, Beardsley R C. An unstructured grid, finite-volume. Three-dimensional, primitive equations ocean model:Application to coastal ocean and estuaries[J]. Journal of Atmospheric and Ocean Technology.2003,20(1):159-186.
[36]P.F. Xue, C.S. Chen, P.X. Ding,et al. Saltwater intrusion into the Changjiang River:A model-guided mechanism study[J].Journal of Geophysical Research,2009,114:C02006, doi:10.1029/2008JC004831.
[37]唐建华,赵升伟,刘玮等.基于FVCOM的强潮海湾三维潮流数值模拟[J].水利水运工程学报,2010,4：81-88.
[38]冯兴如,杨德周,尹宝树. FVCOM在龙口海域潮汐潮流模拟中的应用研究[J].海洋科学,2010(6)：
[39]傅国伟,程声通.水污染控制规划[M].北京：清华大学出版社,1985.
[40]谢永明.环境水质模型概论[M].中国科技大学出版社,1996.
[41]申霞.基于POM的近海三维水质模型研究及其应用[D].河海大学硕士论文,2006.
[42]朱建荣,朱首贤.ECOM模式的改进及在长江河口、杭州湾及邻近海区的应用[J].海洋与湖沼.2003,34(4)：364-374.
[43]刘成,李行伟,韦鹤平.长江口水动力及污水扩散稀释模拟[J].海洋与湖 沼.2003,34(5)：474-483.
[44]浙江省水文志编篡委员会.浙江省水文志[M].中华书局,2000.4.
[45]浙江省水利志编篡委员会,浙江省水利志[M].中华书局,1998.10.
[46]钱塘江志编篡委员会,钱塘江志[M].方志出版社,1998.4
[47]韩曾萃,戴泽蘅,李光炳等著,钱塘江河口治理开发[M].中国水利水电出版社,2003.9
[48]瓯江志编篡委员会,瓯江志[M].水利电力出版社,1995.3.
[49]甬江志编委会,甬江志[M].中华书局,2000.1.
[50]鳌江志编篡委员会,鳌江志[M].中华书局,1999.
[51]许雪峰,羊天柱.三门湾多个围垦工程的整体影响数学模型研究[J].海洋学研究,2006,(S1)：49-59.
[52]姚炎明,陈吉余,陈永平等.温州湾水流及其与地形的关系[J].华东师范大学学报(自然科学版),1998(02)：61-68.
[53]杨成浩,廖光洪,罗锋等.乐清湾的潮位、潮流和余流特征[J].海洋学研究,2010(28)：1-13.
[54]张明进,张应华.SMS水动力学软件[J].水道港口,2006,27(1)：57-59.
[55]吴修广,程文龙、赵鑫等.乐清湾进港航道工程海床稳定性及回淤分析专题研究[M].浙江省水利河口研究院报告。
[56]曾剑,谢亚力、赵鑫等.浙江浙能台州第二发电厂工程可行性研究风浪潮、岸滩稳定及循环冷却水温排放数值模拟专题研究[M].浙江省水利河口研究院报告。
[57]潘存鸿,黄世昌、曾剑等,杭州湾大桥工程可行性研究水文分析计算、防洪研究报告[M].浙江省水利河口研究院报告。
[58]郭为民等,钱塘江河口水环境容量及纳污总量控制研究水质取样报告[M].浙江省水利河口研究院报告。
[2]邹志利,海岸动力学[M].人民交通出版社,2009.9.
[3]陈吉作主编,中国海岸侵蚀概要[M].海洋出版社,2010.7.
[4]朱学明,中国近海潮汐潮流的数值模拟与研究[D],中国海洋大学硕士论文,2009.
[5]申霞,基于POM的近海三维水质模型研究及其应用[D],河海大学硕士论文,2006.
[6]章卫胜,中国近海潮波运动数值模拟[D],河海大学硕士论文,2010.
[7]吴伦宇,基于FVCOM的浪、流、泥沙模型耦合及其应用[D],中国海洋大学博士论文,2010.
[8]曹德明,方国洪.杭州湾潮汐潮流的数值计算[J].海洋和湖沼,1986,17(2)：93-101.
[9]曹德明,方国洪.杭州湾和钱塘江潮波的联合数值模型[J].海洋学报,1988,10(5)：521-530.
[10]许卫忆,苏纪兰.杭州湾二维潮波计算及底质分布的动力成因[J].海洋与湖沼,1986(6)：52-56.
[11]曹欣中,唐龙妹张月秀.宁波,舟山内海域实测海流分析及潮流场的数值模拟[J].东海海洋,1996,14(2)：1-9.
[12]董礼先,苏纪兰.象山港潮波响应和变形研究.Ⅱ.象山港潮波数值研究[J].海洋学报,1999,21(2)：1-8.
[13]李孟国,王正林.瓯江口潮流数值模拟[J].长江科学院院报,2002,19(2)：19-22.
[14]李孟国,时钟,范文静.瓯江口外洞头岛北岱后涂围垦工程潮流数值模拟研究[J].海洋通报,2005,24(1)：1-7.
[15]陈晒睿,朱建荣,吴辉等.无结构网格二维河口海岸水动力数值模式的建立及其应用[J].海洋学报,2010,32(2)：31-39
[16]孙文心.超浅海风暴潮的数值模拟.海洋学报[J],1979,1(2)：193-211。
[17]吴培木,许永,李燕初等.台湾海峡台风暴潮非线性数值计算.海洋学报[J],1981,3(1)： 27～43.
[18]陈长胜,秦曾灏.江浙沿海台风暴潮的动力分析.海洋学报[J],1985,7(3)：265-275.
[19]C.S.Chen, Z.H.Qin. Numerical simulation of typhoon surges along the east coast of Zhejiang and Jiangsu Provinces. Advances in Atmosphere Science[J],1985,2(1):8-9.
[20]X.N Wang,F.J Yu, Q.J Yin. Research of Application of Numerical Models of Typhoon Surges in China Seas. The Special Issue of MAUSAM,1997:595-608
[21]赵有皓,张君伦.风暴潮实用数值预报的研发.海洋工程,2001,19(3)：102-107.
[22]Hassan S.Mashriqui,G.Paul Kemp, and Ivor van Heerden. Experimental storm surge simulations for hurricane Katrina.In:Y. J. Xu&V.P.Singh, Coastal Environment and Water Quality. LLC, Highlands Ranch,Co 80163-0026,USA:Water resources Publications,2006:481-490.
[23]Jian Shen,Keqi Zhang.Improved prediction of storm surge inundation with a high-resolution unstructured grid model.Journal of Coastal Research,2006,Vol.22,No.6:1309-1319.
[24]张涤明,陈良圭,乔林.近海水域海洋流体动力学方程的数值模拟[J].力学进展,1988,18(3)：322-331.
[25]李身铎,顾思美.杭州湾潮波三维数值模拟[J].海洋与湖沼,1993,24(1)：7-15.
[26]韩国其,汪德灌,许协庆.潮汐河口三维水流数值模拟[J].水利学报,1989,12(5)：42-60.
[27]万振文,乔方利,袁业立.渤,黄,东海三维潮波运动数值模拟[J].海洋与湖沼,1998,29(6)：611-616.
[28]W.B.Guan,L.A..Wong, D.F. Xu,2001. Modeling nitrogen and phosphorus cycles and dissolved oxygen in the Zhujiang River estuary. Acta Oceanologica Sinica,20(4):493-514.
[29]史峰岩,朱首贤.杭州湾,长江口余流及其物质输运作用的模拟研究[J].海洋学报,2000,22(5)：1一12.
[30]华祖林,顾莉,查玉含,刘晓东.基于COHERENS模型的污染物质输运数值模拟[J].环境科学与技术,2009,32(4)：14-18.
[31]周大成.椒江河口大潮潮流特性的数值分析[J].浙江水利水电专科学校学报,2001,13(3)：13-26.
[32]杨陇慧,朱建荣.长江口杭州湾及邻近海区潮汐潮流场三维数值模拟[J].华东师范大学学报：自然科学版.2001(3)：74-84.
[33]陈倩,黄大吉,章本照.浙东近海潮汐潮流的数值模拟[J].海洋学报,2003,25(5)：9-20.
[34]张威.浙东近海海域潮波数值模拟[D],浙江大学硕士论文,2007.
[35]C.S. Chen,H.D. Liu, Beardsley R C. An unstructured grid, finite-volume. Three-dimensional, primitive equations ocean model:Application to coastal ocean and estuaries[J]. Journal of Atmospheric and Ocean Technology.2003,20(1):159-186.
[36]P.F. Xue, C.S. Chen, P.X. Ding,et al. Saltwater intrusion into the Changjiang River:A model-guided mechanism study[J].Journal of Geophysical Research,2009,114:C02006, doi:10.1029/2008JC004831.
[37]唐建华,赵升伟,刘玮等.基于FVCOM的强潮海湾三维潮流数值模拟[J].水利水运工程学报,2010,4：81-88.
[38]冯兴如,杨德周,尹宝树. FVCOM在龙口海域潮汐潮流模拟中的应用研究[J].海洋科学,2010(6)：
[39]傅国伟,程声通.水污染控制规划[M].北京：清华大学出版社,1985.
[40]谢永明.环境水质模型概论[M].中国科技大学出版社,1996.
[41]申霞.基于POM的近海三维水质模型研究及其应用[D].河海大学硕士论文,2006.
[42]朱建荣,朱首贤.ECOM模式的改进及在长江河口、杭州湾及邻近海区的应用[J].海洋与湖沼.2003,34(4)：364-374.
[43]刘成,李行伟,韦鹤平.长江口水动力及污水扩散稀释模拟[J].海洋与湖 沼.2003,34(5)：474-483.
[44]浙江省水文志编篡委员会.浙江省水文志[M].中华书局,2000.4.
[45]浙江省水利志编篡委员会,浙江省水利志[M].中华书局,1998.10.
[46]钱塘江志编篡委员会,钱塘江志[M].方志出版社,1998.4
[47]韩曾萃,戴泽蘅,李光炳等著,钱塘江河口治理开发[M].中国水利水电出版社,2003.9
[48]瓯江志编篡委员会,瓯江志[M].水利电力出版社,1995.3.
[49]甬江志编委会,甬江志[M].中华书局,2000.1.
[50]鳌江志编篡委员会,鳌江志[M].中华书局,1999.
[51]许雪峰,羊天柱.三门湾多个围垦工程的整体影响数学模型研究[J].海洋学研究,2006,(S1)：49-59.
[52]姚炎明,陈吉余,陈永平等.温州湾水流及其与地形的关系[J].华东师范大学学报(自然科学版),1998(02)：61-68.
[53]杨成浩,廖光洪,罗锋等.乐清湾的潮位、潮流和余流特征[J].海洋学研究,2010(28)：1-13.
[54]张明进,张应华.SMS水动力学软件[J].水道港口,2006,27(1)：57-59.
[55]吴修广,程文龙、赵鑫等.乐清湾进港航道工程海床稳定性及回淤分析专题研究[M].浙江省水利河口研究院报告。
[56]曾剑,谢亚力、赵鑫等.浙江浙能台州第二发电厂工程可行性研究风浪潮、岸滩稳定及循环冷却水温排放数值模拟专题研究[M].浙江省水利河口研究院报告。
[57]潘存鸿,黄世昌、曾剑等,杭州湾大桥工程可行性研究水文分析计算、防洪研究报告[M].浙江省水利河口研究院报告。
[58]郭为民等,钱塘江河口水环境容量及纳污总量控制研究水质取样报告[M].浙江省水利河口研究院报告。