围海工程堵口水力计算与研究
|
摘要
我国人多地少,土地资源非常紧张。当前又正处于工业化、城镇化快速发展时期,建设用地供需矛盾十分突出,人口增加及城市化进程的加快,对土地的需求越来越大。围海造地对沿海港口及城市的发展有非常重要的作用,是我国尤其是东部沿海地区解决土地紧张问题和拓展生存空间较好的途径之一,世界上许多国家和地区都开发海岸带资源来解决人口与土地供给之间的矛盾。
围海工程一般包括海堤、水闸、堵口三部分。海堤是围海工程主体,直接承受风浪、风暴潮、急流作用,工程量较大,以软基筑堤和堵口闭气为其重要技术问题；围海工程中的水闸主要起到分流做用,有双向过流的特点,海堤和水闸都是围海工程先期施工的部分；堵口是海堤修筑的最后阶段,堵口合拢是围海工程成败的关键,堵口施工中所遇到的技术问题,是整个围海工程中最复杂和最困难的问题。山于此阶段所遇到的水力条件十分恶劣,口门落差可达2-3米,流速最大可达到6-8米/秒,且水流极其复杂,根据我国《围海工程技术规范》规定围海堵口必须进行龙口水力计算。
本文首先介绍围海工程的发展进程与国内外围海的开发现状,介绍了我国目前围海工程的特点及在围海工程实践中积累的经验和理论,利用CFX软件对堵口工程龙口水力要素进行数值模拟,重点研究了堵口快速合拢期间,流速值在二维平面的分布规律以及流速值随潮流的变化规律,提出合理合拢方案。并运用该数值模拟方法对某海上围堰工程进行堵口水力计算,模拟堵口过程中的水动力情况,为制定堵口合拢方案提供科学依据。
China's inadequate per-capita share of land resources is becoming more and more serious. Particularly, when china's economic development and social progress accelerated in the past decades. The gap between land supply and demand is still expanding. Land reclamation from sea is a better method to solve the problem especially in china's eastern coastal areas. Many countries and regions in the world have exploited their coastal zones and tide land resources to increase provision of land.
Land reclamation project encircles the sea to make land, It has three main sections:dike、sluice gate and dyke enclosure.Sea dyke is the main part of the reclamation project,which bears the load of waves、storms and streams. Soft clay ground treatment and infiltration closure are the two key technique problems. Sluice gate's main function is to separate the stream of entrance, usually constructed before enclosure begins, together with dyke's main section. Dyke enclosure is the most crucial and the latest part of sea reclamation engineering. At this stage,dyke enclosure confronts very severe hydraulic conditions, head drop of the flow can reach 2-3 m, the maximum flow velocity probably can reach up to 6-8 m/s. Appropriate enclosure steps and precise hydraulic calculation are needed before enclosure project begins.
This dissertation firstly introduces the present status of reclamation project. Secondly, analyzes the character of China's reclamation project, summarizes the hydraulic calculation theory and practical experience;Thirdly, uses commercial software CFX to simulate the hydraulic parameter of fast enclosure stage, especially simulates velocity distribution patterns on horizontal two-dimensional plane. Preliminarily studies the law or velocity variation pattern in pace with tidal action. Lastly, using the CFX model, calculates the hydraulic parameter of a maritime park project.
围海工程一般包括海堤、水闸、堵口三部分。海堤是围海工程主体,直接承受风浪、风暴潮、急流作用,工程量较大,以软基筑堤和堵口闭气为其重要技术问题；围海工程中的水闸主要起到分流做用,有双向过流的特点,海堤和水闸都是围海工程先期施工的部分；堵口是海堤修筑的最后阶段,堵口合拢是围海工程成败的关键,堵口施工中所遇到的技术问题,是整个围海工程中最复杂和最困难的问题。山于此阶段所遇到的水力条件十分恶劣,口门落差可达2-3米,流速最大可达到6-8米/秒,且水流极其复杂,根据我国《围海工程技术规范》规定围海堵口必须进行龙口水力计算。
本文首先介绍围海工程的发展进程与国内外围海的开发现状,介绍了我国目前围海工程的特点及在围海工程实践中积累的经验和理论,利用CFX软件对堵口工程龙口水力要素进行数值模拟,重点研究了堵口快速合拢期间,流速值在二维平面的分布规律以及流速值随潮流的变化规律,提出合理合拢方案。并运用该数值模拟方法对某海上围堰工程进行堵口水力计算,模拟堵口过程中的水动力情况,为制定堵口合拢方案提供科学依据。
China's inadequate per-capita share of land resources is becoming more and more serious. Particularly, when china's economic development and social progress accelerated in the past decades. The gap between land supply and demand is still expanding. Land reclamation from sea is a better method to solve the problem especially in china's eastern coastal areas. Many countries and regions in the world have exploited their coastal zones and tide land resources to increase provision of land.
Land reclamation project encircles the sea to make land, It has three main sections:dike、sluice gate and dyke enclosure.Sea dyke is the main part of the reclamation project,which bears the load of waves、storms and streams. Soft clay ground treatment and infiltration closure are the two key technique problems. Sluice gate's main function is to separate the stream of entrance, usually constructed before enclosure begins, together with dyke's main section. Dyke enclosure is the most crucial and the latest part of sea reclamation engineering. At this stage,dyke enclosure confronts very severe hydraulic conditions, head drop of the flow can reach 2-3 m, the maximum flow velocity probably can reach up to 6-8 m/s. Appropriate enclosure steps and precise hydraulic calculation are needed before enclosure project begins.
This dissertation firstly introduces the present status of reclamation project. Secondly, analyzes the character of China's reclamation project, summarizes the hydraulic calculation theory and practical experience;Thirdly, uses commercial software CFX to simulate the hydraulic parameter of fast enclosure stage, especially simulates velocity distribution patterns on horizontal two-dimensional plane. Preliminarily studies the law or velocity variation pattern in pace with tidal action. Lastly, using the CFX model, calculates the hydraulic parameter of a maritime park project.
引文
[1]裘江海.开发滩涂资源拓展生存空间[J].杭州：浙江水利水电专科学校学报,2005(3)：5-6.
[2]严恺主编.海岸工程[M].北京：海洋出版社,2002(2)674-683.
[3]陈吉余主编.中国围海工程[M].北京：中国水利水电出版社,2001.11：77-223.
[4]国家行动计划编委会.保护海洋环境免受陆基活动影响全球行动方案.国家环境保护总局.中国环境科学研究院,2006年4月9日.
[5]陈吉余.开发浅海滩涂资源拓宽我国的生存空间[J].北京：中国工程科学,2000,2(3)：27-31.
[6]冯利华,鲍毅新.滩涂围垦的负面影响与可持续发展策略[J].青岛：海洋科学,2004,28(4)：76-77.1-20.
[7]严恺,周家苞.中国海岸带开发利用[J].北京：中国水利,水利电力出版社,1991.
[8]潘丽红,朱建荣,等.河口大型围垦工程中龙口水动力特点[J].南京：海洋工程,2010,28(1)：110-116.
[9]严恺,梁其荀.海港工程[M].北京：海洋出版社,1996.6：689-696,674-676.
[10]陈德春,吴继伟,李宇,等.围海工程堵口合龙技术研究[J].南京：河海大学学报(自然科学版),2002,30(5)：67-70.
[11]福建省水力水电厅.福建省围垦工程设计技术规程[M].附录L,1992：121-125.
[12]王光谦,杨文俊,夏军强,等.三峡工程明渠截流水流数学模型研究及其应用(Ⅰ)：模型的建立与率定[J].武汉：长江科学院院报,2005,22(2)：1-4.
[13]夏军强,王光谦,杨文俊,等.三峡工程明渠截流水流数学模型研究及其应用(Ⅱ)：方案计算与反演计算[J].武汉：长江科学院院报,2005,22(3)：1-5.
[14]迟耀瑜,朱旦生,王秀民,等.立堵截流规律再探讨[J].泥沙研究,1998,(4)：56-64.
[15]李家星,陈立德.水力学(上册)[M].南京：河海大学出版社,1996：335-336；398-401.
[16]李炜,等.水力计算手册[M].北京：中国水利水电出版社,2006：73-76.
[17]汪龙腾,等.围海堵口工程龙口水力条件及堵口程序[J].南京：华东水力学院学报,1979,(4).
[18]李开运.围海工程三十五年[J].南京：华水科技情报,1984,(3).
[19]谭浩强,田淑清Fortran语言程序设计[M].北京：高等教育出版社,1996.10：1-2.
[20]刘维东.一、二维围海堵口水力计算的应用研究[硕士学位论文].南京：河海大学交通与海洋工程学院,2007.
[21]吴继伟.围海工程堵口水力计算软件的开发与应用[硕士学位论文].南京：河海大学交通与海洋工程学院,2003.
[22]李家星、赵振兴.水力学(上、下)[M].南京：河海大学出版社,2001.
[23]王绍成、虞和莹.河流动力学[M].人民交通出版社,1991.
[24]Waltes, R. A. and Cheng, R. T. A Two-Dimensional Hydrodynamics model of a Tidal Estuary, Finete Elements in Water Resources, PENTECH PRESS, London: Plymoutuh,1978.
[25]史宏达,李水鱼.青岛浮山湾数值化海区建立的研究[J].海岸工程,2006,25(2)：1-12.
[26]戴会超,槐文信,等.水利水电工程水流精细模拟理论与应用[M].北京：科学出版社2006.4：17-20
[27]季永兴,卢永金,姚华生.浦东国际机场围海大堤龙口水力数学模型研究[J].水利水电科技进展,2000,20(6)：36-38.
[28]王斌.曹妃甸围海造地二维潮流数值计算及滩槽稳定性研究[硕士学位论文].南京：河海大学交通与海洋工程学院,2007.
[29]李人宪.有限体积法基础[M].北京：国防工业出版社2005.7：23-27.
[30]刘臻.岸式振荡水柱波能发电装置的实验及数值模拟研究[硕士学位论文].青岛：中国海洋大学工程学院,2008.
[31]约翰D.安德森.计算流体力学基础及其应用[M].北京：机械工业出版社2007.6：116-121.
[32]刘超群.多重网格法及其在计算流体力学中的应用.北京：清华大学出版社1995.6：115-124.
[33]阎超.计算流体力学方法及应用.北京：北京航空航天大学出版社2006.6：180-185.
[34]Peter Monk. A Mixed Finite Element Method for the Biharmonic Equation[J]. SIAM J NUMER ANAL,1987,24(4):737-749.
[35]Jim Douglas Jr, Thomas F Russell. Numerical Method for Convection-dominated Diffusion Problems Based on Combining the Method of Characteristics with Finite Element or Finite Difference Procedures[J]. SIAM J Numer Anal,1982,19(5):871-885.
[36]Russell T F. Time-stepping along Characteristics with Incomplete Iteration for a Galerkin Approximation of Miscible Displacement in Porous Media[J].SIAM J Numer Anal,1985,22:970-1013.
[37]V Sundar,P M Koola. Dynamic pressure on inclined cylinders due to freak waves [J]. Ocean Engineering,1999,26:841-863.
[38]Hirt C W, Nichols B D. Volume of Fluid (VOF) Mehtod for the Dynamics of Free Boundaries. Journal of Computational Physics,1981,39(1):201-225.
[39]Rosenblum L J.1994. Research Issues in Scientific Visualization. IEEE Computer Graphics and Applications. March,61-85.
[40]李孟国,蒋德才.关于波浪缓坡方程的研究.海洋通报[J].1999,18：70-92.
[41]文圣常,余宙文.海浪理论与计算原理.北京：科学出版社.1984.1-20.
[42]谭维炎.计算水动力学—有限体积法的应用.北京：清华大学出版社.1998.1-32.
[43]钦文婷Boussinesq方程数学模型的改进及其工程应用：[硕士学位论文].天津：天津大学.2003.
[44]赵艳艳.金山湖二维水流水质数学模型研究[硕士学位论文].南京：河海大学,2005.
[45]Gabbai, R. D. and H. Benaroya, An overview of modeling and experiments of vortex-induced. Journal of Sound and Vibration,2005,282:575-616.
[46]Zdravkovich, M. M., Flow Around Circular Cylinders, Vol.2:Applications.2003, Oxford University Press:Oxford, England.
[47]Sinmons T J. Development of three-dimensional numerical model of the Great Lakes[J]. Calnlan Water Branch, Sci Ser,1973(12):269.
[48]王同科,时红霞.二维不可压Navier-stokes方程的特征混合有限元算法[J].山东大学学报(理学版),2002(6)：221-224
[49]Leederts J J, liu S K. A Three-Dimensional Turbulent Energy Model for Nonhomogeneous Estuaries and Coastal Sea System, Hydraulics of Estuaries and Fjouls,1978.
[50]Warzee G,Sterling W A, Adam Y. Application of the Finite Element Method to Periodic Tidal computation, Finite Elements in Water Resources,1980.
[51]黄胜,卢启苗等.河口动力学.北京：水利水电出版社,1995.
[52]周发毅,等.大规模多岛屿海域潮流场的数值模拟[J].海洋工程,1995(4)
[53]王船海,等.河道二维非恒定流计算[J].河海大学学报,1987(4).
[54]王金瑞,卢玫君,伊华伟.平面二维河道水流泥沙数学模型[J].华东水利水电学院院报,1999(6)：13-17.
[2]严恺主编.海岸工程[M].北京：海洋出版社,2002(2)674-683.
[3]陈吉余主编.中国围海工程[M].北京：中国水利水电出版社,2001.11：77-223.
[4]国家行动计划编委会.保护海洋环境免受陆基活动影响全球行动方案.国家环境保护总局.中国环境科学研究院,2006年4月9日.
[5]陈吉余.开发浅海滩涂资源拓宽我国的生存空间[J].北京：中国工程科学,2000,2(3)：27-31.
[6]冯利华,鲍毅新.滩涂围垦的负面影响与可持续发展策略[J].青岛：海洋科学,2004,28(4)：76-77.1-20.
[7]严恺,周家苞.中国海岸带开发利用[J].北京：中国水利,水利电力出版社,1991.
[8]潘丽红,朱建荣,等.河口大型围垦工程中龙口水动力特点[J].南京：海洋工程,2010,28(1)：110-116.
[9]严恺,梁其荀.海港工程[M].北京：海洋出版社,1996.6：689-696,674-676.
[10]陈德春,吴继伟,李宇,等.围海工程堵口合龙技术研究[J].南京：河海大学学报(自然科学版),2002,30(5)：67-70.
[11]福建省水力水电厅.福建省围垦工程设计技术规程[M].附录L,1992：121-125.
[12]王光谦,杨文俊,夏军强,等.三峡工程明渠截流水流数学模型研究及其应用(Ⅰ)：模型的建立与率定[J].武汉：长江科学院院报,2005,22(2)：1-4.
[13]夏军强,王光谦,杨文俊,等.三峡工程明渠截流水流数学模型研究及其应用(Ⅱ)：方案计算与反演计算[J].武汉：长江科学院院报,2005,22(3)：1-5.
[14]迟耀瑜,朱旦生,王秀民,等.立堵截流规律再探讨[J].泥沙研究,1998,(4)：56-64.
[15]李家星,陈立德.水力学(上册)[M].南京：河海大学出版社,1996：335-336；398-401.
[16]李炜,等.水力计算手册[M].北京：中国水利水电出版社,2006：73-76.
[17]汪龙腾,等.围海堵口工程龙口水力条件及堵口程序[J].南京：华东水力学院学报,1979,(4).
[18]李开运.围海工程三十五年[J].南京：华水科技情报,1984,(3).
[19]谭浩强,田淑清Fortran语言程序设计[M].北京：高等教育出版社,1996.10：1-2.
[20]刘维东.一、二维围海堵口水力计算的应用研究[硕士学位论文].南京：河海大学交通与海洋工程学院,2007.
[21]吴继伟.围海工程堵口水力计算软件的开发与应用[硕士学位论文].南京：河海大学交通与海洋工程学院,2003.
[22]李家星、赵振兴.水力学(上、下)[M].南京：河海大学出版社,2001.
[23]王绍成、虞和莹.河流动力学[M].人民交通出版社,1991.
[24]Waltes, R. A. and Cheng, R. T. A Two-Dimensional Hydrodynamics model of a Tidal Estuary, Finete Elements in Water Resources, PENTECH PRESS, London: Plymoutuh,1978.
[25]史宏达,李水鱼.青岛浮山湾数值化海区建立的研究[J].海岸工程,2006,25(2)：1-12.
[26]戴会超,槐文信,等.水利水电工程水流精细模拟理论与应用[M].北京：科学出版社2006.4：17-20
[27]季永兴,卢永金,姚华生.浦东国际机场围海大堤龙口水力数学模型研究[J].水利水电科技进展,2000,20(6)：36-38.
[28]王斌.曹妃甸围海造地二维潮流数值计算及滩槽稳定性研究[硕士学位论文].南京：河海大学交通与海洋工程学院,2007.
[29]李人宪.有限体积法基础[M].北京：国防工业出版社2005.7：23-27.
[30]刘臻.岸式振荡水柱波能发电装置的实验及数值模拟研究[硕士学位论文].青岛：中国海洋大学工程学院,2008.
[31]约翰D.安德森.计算流体力学基础及其应用[M].北京：机械工业出版社2007.6：116-121.
[32]刘超群.多重网格法及其在计算流体力学中的应用.北京：清华大学出版社1995.6：115-124.
[33]阎超.计算流体力学方法及应用.北京：北京航空航天大学出版社2006.6：180-185.
[34]Peter Monk. A Mixed Finite Element Method for the Biharmonic Equation[J]. SIAM J NUMER ANAL,1987,24(4):737-749.
[35]Jim Douglas Jr, Thomas F Russell. Numerical Method for Convection-dominated Diffusion Problems Based on Combining the Method of Characteristics with Finite Element or Finite Difference Procedures[J]. SIAM J Numer Anal,1982,19(5):871-885.
[36]Russell T F. Time-stepping along Characteristics with Incomplete Iteration for a Galerkin Approximation of Miscible Displacement in Porous Media[J].SIAM J Numer Anal,1985,22:970-1013.
[37]V Sundar,P M Koola. Dynamic pressure on inclined cylinders due to freak waves [J]. Ocean Engineering,1999,26:841-863.
[38]Hirt C W, Nichols B D. Volume of Fluid (VOF) Mehtod for the Dynamics of Free Boundaries. Journal of Computational Physics,1981,39(1):201-225.
[39]Rosenblum L J.1994. Research Issues in Scientific Visualization. IEEE Computer Graphics and Applications. March,61-85.
[40]李孟国,蒋德才.关于波浪缓坡方程的研究.海洋通报[J].1999,18：70-92.
[41]文圣常,余宙文.海浪理论与计算原理.北京：科学出版社.1984.1-20.
[42]谭维炎.计算水动力学—有限体积法的应用.北京：清华大学出版社.1998.1-32.
[43]钦文婷Boussinesq方程数学模型的改进及其工程应用：[硕士学位论文].天津：天津大学.2003.
[44]赵艳艳.金山湖二维水流水质数学模型研究[硕士学位论文].南京：河海大学,2005.
[45]Gabbai, R. D. and H. Benaroya, An overview of modeling and experiments of vortex-induced. Journal of Sound and Vibration,2005,282:575-616.
[46]Zdravkovich, M. M., Flow Around Circular Cylinders, Vol.2:Applications.2003, Oxford University Press:Oxford, England.
[47]Sinmons T J. Development of three-dimensional numerical model of the Great Lakes[J]. Calnlan Water Branch, Sci Ser,1973(12):269.
[48]王同科,时红霞.二维不可压Navier-stokes方程的特征混合有限元算法[J].山东大学学报(理学版),2002(6)：221-224
[49]Leederts J J, liu S K. A Three-Dimensional Turbulent Energy Model for Nonhomogeneous Estuaries and Coastal Sea System, Hydraulics of Estuaries and Fjouls,1978.
[50]Warzee G,Sterling W A, Adam Y. Application of the Finite Element Method to Periodic Tidal computation, Finite Elements in Water Resources,1980.
[51]黄胜,卢启苗等.河口动力学.北京：水利水电出版社,1995.
[52]周发毅,等.大规模多岛屿海域潮流场的数值模拟[J].海洋工程,1995(4)
[53]王船海,等.河道二维非恒定流计算[J].河海大学学报,1987(4).
[54]王金瑞,卢玫君,伊华伟.平面二维河道水流泥沙数学模型[J].华东水利水电学院院报,1999(6)：13-17.