桩和桩群的水流阻力及其在潮流数值模拟中的应用
摘要
桩基码头是港口工程中最常用的码头结构型式之一，应用极其广泛。由于桩基所在区域水流结构的特殊性，目前理论上研究尚少，实践上则通常采用物理模型试验或数值模拟来了解桩基码头前的水流情况。由于模型试验范围一般很大，物理模型往往都采用变态模型，而数值模拟则由于尚无合适的数学模型，因此，目前在这两种模型中，通常都将桩基码头简化为直立岸壁，通过模拟来预演工程实施后流态改变的情况。此种概化方法不能反映桩基对水流的影响，无法较好地描述桩基码头前沿真实的水流情况。
桩群的水流阻力比单桩的水流阻力要复杂得多，理论上解决桩群水流阻力目前还较困难。因此，本研究采用局部物理模型试验，通过分析研究，提出在桩式码头数学模型计算中，考虑阻力系数的计算方法。并提出在数值模拟控制方程中，如何考虑加入桩(群)影响的阻力项。
本文研究包括物理模型、数值模拟和工程应用三大部分。
1．物理模型试验包括三类模型：
(1)单桩和群桩阻力的水槽试验；试验中采用不同桩截面形状和不同排列方式。
(2)水流中墩群阻力的物理模型试验(在水池内进行)；
(3)码头群桩阻力的宽水槽试验；
通过上述三个物理模型试验，得出不同截面形状和不同排列时的阻力系数，并应用于数值模拟中。
大连理工大学博士学位论文
2．数值模拟：
通过理论分析，建立了桩体影响下的平面二维潮流数值模拟基本控制方程，
并建立了数值模拟水流槽。在数值模拟槽中，桩基区域用本研究中物模所得的
附加糙率系数代替桩基，水深也按过水断面相等原理修正。由此计算得出的码
头前沿及港池调头区水域的涨、落潮流速变化与对应的物模水槽试验结果一致，
证明用附加阻力系数和修正水深的方法可以很好地模拟桩群对水流的影响。
3．工程应用：
本文以上述研究成果为基础，对上海国际航运中心洋山港区一期工程5个
桩基泊位进行二维潮流数值模拟计算，在该模拟系统中，引入本项研究成果用
以模拟码头桩群对水流的影响，使研究结果更加符合实际，并已提供给工程设
计单位应用。
Piled wharf (wharf with pile group underwater structure) is a kind of structure, which is very commonly used in port engineering. The flow pattern in the region of pile group is complicate, and is still rarely to be studied at present. Usually, in a real project, the flow pattern in front of the piled wharf is studied by means of physical model test or numerical simulation. Because the water area to be studied in such project is generally so large that an abnormal model should be used in the physical model test and also there is no suitable numerical model for such kind of structure in numerical simulation up to now. In order to preview the variations of flow pattern in the water region after construction, the piled wharf is usually simplified as a vertical quay in these two kinds of models at present. Making such simplification can not reflect the effects of pile group on water flow completely and so can not describe the real flow pattern in front of the piled wharf either.
The flow resistance of pile group is more complicated than that of single pile. It is difficulty to solve flow resistance of pile group theoretically at present. A set of physical model tests were made here and a calculation method considering resistance coefficient of pile group was applied in the numerical simulation of piled wharf through analysis and research. In the governing equation of numerical simulation, factors of resistance considering effects of pile or pile group were applied.
Three parts of research work are included in this paper. They are a set of physical model tests, numerical simulation and project application.
1. Three physical model experiments:
(1) The tests of the resistance of single pile and group pile in water flume: different shape of pile section with different array style were adopted in the experiment;
(2) Model tests of pile group resistance in water flow in a water basin ;
(3) Model tests for a piled wharf in a wide water flume.
The resistance coefficients with different shape section and different array style have been obtained through above-mentioned three physical model experiments and have been applied in numerical simulation as well.
2. Numerical simulation:
The fundamental governing equation of two-dimensional tidal flow under the influence of pile body has been induced through theoretical analysis, and a numerical flow flume for simulation has been built as well. Within the numerical simulation flume, the pile group has been replaced by the additional roughness coefficient obtained from the physical model tests in pile foundation region, and the water depth has been revised according to the principle for equality of flow section as well. Then the variation of flow velocity for flood and ebb tide flow in the water region of the wharf frontage and vessel turning basin has been obtained. The obtained results are quite coincident with the results of the equivalent physical model tests. It proves that using the method of additional resistance coefficients and revised water depth can simulate the effects on water flow by pile group quite well.
3. Project application:
Depending on the above research results, the numerical simulation method for two dimensional water flow has been applied in calculating the five piled wharf berths of the first phase project of Yangshan port of Shanghai international shipping center. Among the simulation system, the above mentioned research results have been used to simulate the effects on the water flow by piled wharf, which makes the result more tallies with the actual situation. The results have been applied in the project design and research.
桩群的水流阻力比单桩的水流阻力要复杂得多，理论上解决桩群水流阻力目前还较困难。因此，本研究采用局部物理模型试验，通过分析研究，提出在桩式码头数学模型计算中，考虑阻力系数的计算方法。并提出在数值模拟控制方程中，如何考虑加入桩(群)影响的阻力项。
本文研究包括物理模型、数值模拟和工程应用三大部分。
1．物理模型试验包括三类模型：
(1)单桩和群桩阻力的水槽试验；试验中采用不同桩截面形状和不同排列方式。
(2)水流中墩群阻力的物理模型试验(在水池内进行)；
(3)码头群桩阻力的宽水槽试验；
通过上述三个物理模型试验，得出不同截面形状和不同排列时的阻力系数，并应用于数值模拟中。
大连理工大学博士学位论文
2．数值模拟：
通过理论分析，建立了桩体影响下的平面二维潮流数值模拟基本控制方程，
并建立了数值模拟水流槽。在数值模拟槽中，桩基区域用本研究中物模所得的
附加糙率系数代替桩基，水深也按过水断面相等原理修正。由此计算得出的码
头前沿及港池调头区水域的涨、落潮流速变化与对应的物模水槽试验结果一致，
证明用附加阻力系数和修正水深的方法可以很好地模拟桩群对水流的影响。
3．工程应用：
本文以上述研究成果为基础，对上海国际航运中心洋山港区一期工程5个
桩基泊位进行二维潮流数值模拟计算，在该模拟系统中，引入本项研究成果用
以模拟码头桩群对水流的影响，使研究结果更加符合实际，并已提供给工程设
计单位应用。
Piled wharf (wharf with pile group underwater structure) is a kind of structure, which is very commonly used in port engineering. The flow pattern in the region of pile group is complicate, and is still rarely to be studied at present. Usually, in a real project, the flow pattern in front of the piled wharf is studied by means of physical model test or numerical simulation. Because the water area to be studied in such project is generally so large that an abnormal model should be used in the physical model test and also there is no suitable numerical model for such kind of structure in numerical simulation up to now. In order to preview the variations of flow pattern in the water region after construction, the piled wharf is usually simplified as a vertical quay in these two kinds of models at present. Making such simplification can not reflect the effects of pile group on water flow completely and so can not describe the real flow pattern in front of the piled wharf either.
The flow resistance of pile group is more complicated than that of single pile. It is difficulty to solve flow resistance of pile group theoretically at present. A set of physical model tests were made here and a calculation method considering resistance coefficient of pile group was applied in the numerical simulation of piled wharf through analysis and research. In the governing equation of numerical simulation, factors of resistance considering effects of pile or pile group were applied.
Three parts of research work are included in this paper. They are a set of physical model tests, numerical simulation and project application.
1. Three physical model experiments:
(1) The tests of the resistance of single pile and group pile in water flume: different shape of pile section with different array style were adopted in the experiment;
(2) Model tests of pile group resistance in water flow in a water basin ;
(3) Model tests for a piled wharf in a wide water flume.
The resistance coefficients with different shape section and different array style have been obtained through above-mentioned three physical model experiments and have been applied in numerical simulation as well.
2. Numerical simulation:
The fundamental governing equation of two-dimensional tidal flow under the influence of pile body has been induced through theoretical analysis, and a numerical flow flume for simulation has been built as well. Within the numerical simulation flume, the pile group has been replaced by the additional roughness coefficient obtained from the physical model tests in pile foundation region, and the water depth has been revised according to the principle for equality of flow section as well. Then the variation of flow velocity for flood and ebb tide flow in the water region of the wharf frontage and vessel turning basin has been obtained. The obtained results are quite coincident with the results of the equivalent physical model tests. It proves that using the method of additional resistance coefficients and revised water depth can simulate the effects on water flow by pile group quite well.
3. Project application:
Depending on the above research results, the numerical simulation method for two dimensional water flow has been applied in calculating the five piled wharf berths of the first phase project of Yangshan port of Shanghai international shipping center. Among the simulation system, the above mentioned research results have been used to simulate the effects on the water flow by piled wharf, which makes the result more tallies with the actual situation. The results have been applied in the project design and research.
引文
[1] 方国洪，杨景飞，渤海潮运动的一个二维数值模型，海洋与潮沼，1985，16(5)337～346．
[2] 杭州湾潮汐和潮流的数值计算，海洋与潮沼，1986，17(2)：93～101．
[3] 曹祖德。胶州湾潮汐潮流的数值计算，海洋科学集刊，第三产业1集，北京科学出版社，1984．157～164．
[4] 方国洪，杨景飞，赵绪才，台湾海峡潮汐和潮流的一个数值模型，海洋学报，1985，7(1)：12～20．
[5] 赵士清，张镜潮，连云港潮流的数值模拟，海洋学报，1981，3(3)：500～515．
[6] 孙献清，尤玉明等厦门港东渡二期工程水流泥沙的初步研究，海洋工程，1991，9(1)：78～86．
[7] 辛文杰，高栏港10万吨级煤炭中转港区方案比选工程潮流数值模拟，第七届全国海岸工程学术讨论会论文集，海洋出版社，1993，93～102．
[8] 徐敏福，黄晋鹏，汕尾港鱼船避风港二维潮流数值计算及泥沙淤积分析，第七届全国海岸工程学术讨论会论文集，海洋出版社，1993，518～524．
[9] 窦振兴，罗远诠等，渤海潮流污染扩散的数值模拟，海洋学报，1982，4(6)：667～678．
[10] 姚炎明，耿兆铨，北仑港三期工程废水排放对北仓海域水质影响的数值分析，海洋通报，1995，14(3)：99～105．
[11] 顾杰，匡翠萍，潮汐流动中污染物排放的一种远区计算模型，海洋工程4(1)：88～96．
[12] 耿兆铨，非恒定流数值解空间失稳问题分析，水利学报，1985，(7)：1～5．
[13] 何少苓，王连祥等，平面二维不恒定流数值计算的若干经验，水利学报，1986，(9)：44～48．
[14] 何少苓，王连祥，窄缝法在二维边界变动水域计算中的应用，水利学报，1986，(12)：11～19．
[15] 辛文杰，河口、海湾平面潮流数值计算中的几个问题，水动力学研究与进展，A辑，1993，8(3)．
[16] Leendertse JJ, RC Alexader, S K Lin A Three-dimensional Model for Estuaries and Coastal Seas. Vol. L, Principles of Computations, R-1417-Owrr. Rand Corp,Santa Monica Califomin, 1973.
[15] J.Lou and P.v. Ridd, Modelling of suspended sediment transport in coastal under waves and currents, Estuarine Coastal and Sheif Science Vol. 45, 1997.
[16] BinLangLin, Roger A. Falconer, Numerical modeling of three-dimensional suspended sediment for estuarine and coastal waters, J. Hydraulic Research, Vol. 34(4)1996.
[17] 孙文心，三维浅海流体力学的一种数值方法——流速分解法，物理海洋数值
计算，河南科学技术出版社，1992．
[18] 方国洪，朱耀华，海洋流体动力学的一种三维数值模式，物理海洋数值计算，河南科学技术出版社，1992．
[19] 李孟国，海岸工程中实用三维流场数学模型研究，海洋通报，1995，14(5)：1～12．
[20] 卢启苗，Onyx，Wai W H．海岸河口三维潮流数学模型，1995，13(4)，47～60．
[21] 孙文心，环境流体力学模型，山东海洋学院学报，1988，18(2)，12～24．
[22] 杨景飞，方国洪，朝鲜海峡潮汐和潮流的数值计算，海洋与湖沼，1987，18(2)：197～207．
[23] 赵士清，张镜潮，一种简便的二维潮流数学模型，水利水运科学研究，1983，(2)：15～23．
[24] 李蓓，二维不规则三角形网格的潮流数学模型，水道港口，1988，(2)：10～15．
[25] 刘智，林秉南等，三角形网格在二维不恒定流计算中的应用，水利学报，1987，(9)：25～33．
[26] 赵士清，长江口潮流的一种数值模式，海洋与湖沼，1985，16(1)：18～27．
[27] 周发毅，郭勇等，大规模多岛屿海域潮流场的数值模拟，海洋工程，1995，13(4)：61～69．
[28] 林建国，时间二次插值求解浅水方程的显式模型，水动力学研究与进展，1989，492)：95～102．
[29] 窦希萍，李提来，三角形网格生成法在海岸工程潮流数学模型中的应用，水利水运科学研究，1995，(3)：65～69．
[30] 窦希萍，李提来，二维潮流数学模型的四边形等参量单元法，海洋工程，1995，13(1)：47～53．
[31] 张廷芳，段小宁，二维潮流数值计算的一种显式模式，大连理工大学学报1992，32(6)：702～706．
[32] 陈伟，温生辉等，厦门港潮动力的数值模拟，台湾海峡，1992，11(4)：328～332．
[33] Flater R A, Heaps N S Tidal computation for the MorecambeBey. Geophys .J.R. Astr. Soc., 1975, (42) :489～517.
[34] Leendertse J J.A Water Quality Simulation Model for Well-mixed Estuaries and Coastal Seas,Vol.1 Principle of Computation,The Rand Corporation RM-6230,1970.
[35] 潘海，方国洪，海洋流体动力学的一种交替方向隐式二维数值模式，海洋学报，1995，17(5)：21～23．
[36] 严以新，诸裕良等，黄茅海域及上游河网地区一、二维联网整体数学模型，第七届全国海岸工程学术讨论会论文集，海洋出版社，1993．
[37] 曹祖德，王桂芬。波浪掀沙、潮流输沙的数值模拟，海洋学报，1993，15(1)：107～118．
[38] 曹祖德，王运洪，水动力泥沙数值模拟，天津大学出版社，1994．
[39] 张存智，杨连武等，极浅海域潮流数值模型，海洋与湖沼，1994，25(5)：671～676．
[40] 韩康，张存智等，大连湾及附近海域潮流场数值模拟，海洋通报，1994，13(4)：20～25．
[41] 华秀菁，吕玉麟，二维浅水域潮流数值模拟的ADI-QUICK格式，水动力学研究与进展，A辑，1996，11(3)：77～92．
[42] Vincenzo Casulli. Semi-implicit finite difference methods for the two-dimensional shallow water equations. J. Comput. Phys.,1990,86:56～74.
[43] Vincenzo Casulli, Filippo Notamicola. An Eulerian-Lagrangian method for tidal current computation. In:Computer Modeling in Ocean Engineering (Sehrefler &.Zienkiewicz,eds),Balkema,Rotterdam, 1988.237～244.
[44] 温生辉，陈季良，厦门附近海域潮流场的数值模拟，海洋学报，1996，18(2)：15～25．
[45] 耿兆铨，二维非恒定流的有限元计算模式，华东水利学院学报，1981，(1)：1～12．
[46] Grotkop G.Finite element analysis of long period water waves.Comp. Meth.in Applied. Mech.and Eng., 1973.(2): 147～160.
[47] 吕玉麟，赖国璋，近海浅水环流问题的数值模拟，大连工学院学报，1981，20(1)：39～52．
[48] 谭维炎，赵枥华，二维浅水渐变不恒定明流的有限元算法及程序包，水利学报，1984，(10)：1～13．
[49] Gray W G. An efficient finite element scheme for two-dimensional surface computation.In: Proceedings of the 1st Intn. Conf. on Finite Elements in Water Resources,Pentech Press, 1977,433～449.
[50] 于天常，二维浅水环流问题的一个有限元模型，海洋与湖沼，1984，15(1)：37～45．
[51] 史金松，计算二维非恒定流的守恒型有限元法，水利学报，1987，(10)：41～48．
[52] 谭维炎，赵枥华，二维非恒定浅水明流的三种算法，水利学报，1985，(9)：1～9．
[53] 李浩麟，易家豪，河口浅水方程的隐式和显式有限元解法，水利水运科学研究，1983，(10)：15～26．
[54] 耿兆铨，二维非恒定流的显式迎流有限元模型，全国第一届水动力学学术会议论文集，北京：海洋出版社，1984．
[55] 张存智，杨连武等，具有潮滩移动边界的浅海环流有限元模型，海洋学报，1990，12(1)：1～13．
[56] 李浩麟等，河口二维潮流单元积分的数值计算，水利学报，1984，(11)：66～72．
[57] 郭庆超，何明民等，控制体积法在二维潮流计算中的应用，水动力学研究与进展，A辑，1995，10(6)：602～609．
[58] 张华庆，乐培九，钦州港港区流场的数值模拟，水利水运科学研究，1994，(1)：81～88．
[59] 董耀华，黄煜龄，河口潮流河段二维非恒定流数模研究及应用，长江科学院院报，1995，12(1)：31～39．
[60] 张东生，蒋勤，江苏北部灌河口悬沙输送数学模型，海洋学报，1991，13(1)：125～136．
[61] 金忠青，N-S方程的数值解及紊流模型，南京：河海大学出版社，1989．
[62] 林秉南，赵雪华等，河口建坝对毗邻海潮波影响的计算，水利学报，1980，(20)：16～26．
[63] 汪德灌，陈新建等，长江口南支流场的平面二维有限差分模拟，河海大学学报，1987，15(5)：34～46．
[64] 何少苓，林秉南，破开算子法在二维潮流计算中的应用，海洋学报，6(2)：260～271．
[65] 何少苓，龚振瀛等，式隐破开算子法在二维潮流计算中的应用，海洋学报，1985，7(2)：225～232．
[66] 吴江航，陈凯麒等，核电站冷却水远区热、核污染数值计算的一种新方法，水利学报，1986，(10)：16～25．
[67] 宋志尧，严以新等，二维潮流方程的一种分裂算子格式及算例，第七届全国海岸工程学术讨论会论文集，海洋出版社，1993，384～392．
[68] 周朦，方国洪，二维长波方程的一个无条件稳定有限差分格式，海洋与湖沼，1988，164～172．
[69] 吴江航，韩庆书，计算流体力学的理论、方法及应用，北京：科学出版社，1988．
[70] 张二骏，张东生等，一种新的数值方法——准分析法在二维非恒定流问题中的应用，海洋学报，1986，8(5)：636～643．
[71] 高飞，张二骏，局部一维准分析法的二维非恒定流的数学模型，水利水运科学研究，1986，(1)：65～75．
[72] 李家星，张镜潮，用潮波能谱法计算二维非恒定流，海洋工程，1988，6(2)：45～55．
[73] 张镜潮，李家星，连云港港区二维潮流计算，海洋工程，1991，10(1)：49～57．
[74] 华祖林，边界拟合坐标下流场计算模式，海洋工程，1994，12(3)：65～72．
[75] 汪德灌，杨艳艳，边界拟合坐标法的应用，河海大学学报，1989，17(1)：50～57．
[76] 王船海，程文辉，河道二维非恒定流计算，河海大学学报，1987，15(3)：39～53．
[77] Wei Wenli,Jin Zhongqing. Numerical solrtion for unsteady 2-D flow using the transformed shallow water quations. Journal of Hydrodynamics, Ser. B, 1995,7(3):65～71.
[78] 赵士清，长江口三维潮流的数值模拟，水利水运科学研究，1985，(1)：23～31．
[79] 张涤明，乔林等，近海环流三维样条分层有限元模式，水动力学研究与进展，1988，3(2)：83～93．
[80] Koutitas C.Three-dimensional models of coastal circulation:an engineering viewpoint. Coastal and Estuarine Sciences,4.American Geophysical Union,Washington,D.C., 1987.
[81] Koutitas C,BO"connor. Modeling 3-D wind-induced flows.J.Hydraulics Div.,ASCE, 1980, 106(11):472～485.
[82] 易家豪，叶雪祥，长江口南港航道三维水流数值模拟，第二次河流泥沙国际学术讨论会论文集，水利电力出版社，1983．
[83] Tee K T. The structure of three-dimensional tide generating currents. Estuarine, Coastal and ShelfScience, 1982,14(1):70～82.
[84] Heaps N S. On the numerical solution of the three-dimensional equations for tides and storm
surges.Mem. Soc. R.Sci.Liege., 1972,6(2): 143～180.
[85] Davies A M.Spectral models in continemtal shelf sea Three-dimensional Coastal oceanography. Ocean Models Norman S. Heaps,ed.,American Geophysical Union,Washington,D. C., 1987,71～106.
[86] 石磊，奚盘根，三维浅海域流体动力学方程的分步杂交解法，青岛海洋大学学报，1995，20(2)：162～172．
[87] 窦振兴，杨连武等，渤海三维潮流数值模拟，海洋学报，1993，15(5)：1～15．
[88] 奚盘根，冯士作等，非线性三维潮波边值问题模型，山东海洋学院学报，1988，18[2(1)]：37～47．
[89] 王晓建，张廷芳，非正交网格有限体积法在三维潮流场计算中的应用，水动力学研究与进展，A辑，1997，12(1)：56～61．
[90] 李孟国，伶仃洋三维流场数值模拟，水动力学研究与进展，A辑，1996，11(3)：342～351．
[91] 刘子龙，王船海等，长江口三维水流模拟，河海大学学报，1996，24(5)：108～110．
[92] Noye J,Stevean M.A three dimensional model of tidal propagatian using transformations and variable grids.Three-dirnen-sional Coastal Ocean Models,AGU, Washington,D.C., 1987.
[93] 李蓓，珠海电厂煤港海区海域二维波浪、潮流、泥沙淤积数学模型，第七届全国海岸工程学术讨论会论文集，海洋出版社，1993，149～158．
[94] 中华人民共和国交通部，海岸与河口潮流泥沙模拟技术规程(JTJ／T233-98)，人民交通出版社，1999．
[95] 李孟国，珠江口伶仃洋三维潮流数学模型，1993．
[96] 李孟国，泥灰礁三维潮流数学模型，2000．
[2] 杭州湾潮汐和潮流的数值计算，海洋与潮沼，1986，17(2)：93～101．
[3] 曹祖德。胶州湾潮汐潮流的数值计算，海洋科学集刊，第三产业1集，北京科学出版社，1984．157～164．
[4] 方国洪，杨景飞，赵绪才，台湾海峡潮汐和潮流的一个数值模型，海洋学报，1985，7(1)：12～20．
[5] 赵士清，张镜潮，连云港潮流的数值模拟，海洋学报，1981，3(3)：500～515．
[6] 孙献清，尤玉明等厦门港东渡二期工程水流泥沙的初步研究，海洋工程，1991，9(1)：78～86．
[7] 辛文杰，高栏港10万吨级煤炭中转港区方案比选工程潮流数值模拟，第七届全国海岸工程学术讨论会论文集，海洋出版社，1993，93～102．
[8] 徐敏福，黄晋鹏，汕尾港鱼船避风港二维潮流数值计算及泥沙淤积分析，第七届全国海岸工程学术讨论会论文集，海洋出版社，1993，518～524．
[9] 窦振兴，罗远诠等，渤海潮流污染扩散的数值模拟，海洋学报，1982，4(6)：667～678．
[10] 姚炎明，耿兆铨，北仑港三期工程废水排放对北仓海域水质影响的数值分析，海洋通报，1995，14(3)：99～105．
[11] 顾杰，匡翠萍，潮汐流动中污染物排放的一种远区计算模型，海洋工程4(1)：88～96．
[12] 耿兆铨，非恒定流数值解空间失稳问题分析，水利学报，1985，(7)：1～5．
[13] 何少苓，王连祥等，平面二维不恒定流数值计算的若干经验，水利学报，1986，(9)：44～48．
[14] 何少苓，王连祥，窄缝法在二维边界变动水域计算中的应用，水利学报，1986，(12)：11～19．
[15] 辛文杰，河口、海湾平面潮流数值计算中的几个问题，水动力学研究与进展，A辑，1993，8(3)．
[16] Leendertse JJ, RC Alexader, S K Lin A Three-dimensional Model for Estuaries and Coastal Seas. Vol. L, Principles of Computations, R-1417-Owrr. Rand Corp,Santa Monica Califomin, 1973.
[15] J.Lou and P.v. Ridd, Modelling of suspended sediment transport in coastal under waves and currents, Estuarine Coastal and Sheif Science Vol. 45, 1997.
[16] BinLangLin, Roger A. Falconer, Numerical modeling of three-dimensional suspended sediment for estuarine and coastal waters, J. Hydraulic Research, Vol. 34(4)1996.
[17] 孙文心，三维浅海流体力学的一种数值方法——流速分解法，物理海洋数值
计算，河南科学技术出版社，1992．
[18] 方国洪，朱耀华，海洋流体动力学的一种三维数值模式，物理海洋数值计算，河南科学技术出版社，1992．
[19] 李孟国，海岸工程中实用三维流场数学模型研究，海洋通报，1995，14(5)：1～12．
[20] 卢启苗，Onyx，Wai W H．海岸河口三维潮流数学模型，1995，13(4)，47～60．
[21] 孙文心，环境流体力学模型，山东海洋学院学报，1988，18(2)，12～24．
[22] 杨景飞，方国洪，朝鲜海峡潮汐和潮流的数值计算，海洋与湖沼，1987，18(2)：197～207．
[23] 赵士清，张镜潮，一种简便的二维潮流数学模型，水利水运科学研究，1983，(2)：15～23．
[24] 李蓓，二维不规则三角形网格的潮流数学模型，水道港口，1988，(2)：10～15．
[25] 刘智，林秉南等，三角形网格在二维不恒定流计算中的应用，水利学报，1987，(9)：25～33．
[26] 赵士清，长江口潮流的一种数值模式，海洋与湖沼，1985，16(1)：18～27．
[27] 周发毅，郭勇等，大规模多岛屿海域潮流场的数值模拟，海洋工程，1995，13(4)：61～69．
[28] 林建国，时间二次插值求解浅水方程的显式模型，水动力学研究与进展，1989，492)：95～102．
[29] 窦希萍，李提来，三角形网格生成法在海岸工程潮流数学模型中的应用，水利水运科学研究，1995，(3)：65～69．
[30] 窦希萍，李提来，二维潮流数学模型的四边形等参量单元法，海洋工程，1995，13(1)：47～53．
[31] 张廷芳，段小宁，二维潮流数值计算的一种显式模式，大连理工大学学报1992，32(6)：702～706．
[32] 陈伟，温生辉等，厦门港潮动力的数值模拟，台湾海峡，1992，11(4)：328～332．
[33] Flater R A, Heaps N S Tidal computation for the MorecambeBey. Geophys .J.R. Astr. Soc., 1975, (42) :489～517.
[34] Leendertse J J.A Water Quality Simulation Model for Well-mixed Estuaries and Coastal Seas,Vol.1 Principle of Computation,The Rand Corporation RM-6230,1970.
[35] 潘海，方国洪，海洋流体动力学的一种交替方向隐式二维数值模式，海洋学报，1995，17(5)：21～23．
[36] 严以新，诸裕良等，黄茅海域及上游河网地区一、二维联网整体数学模型，第七届全国海岸工程学术讨论会论文集，海洋出版社，1993．
[37] 曹祖德，王桂芬。波浪掀沙、潮流输沙的数值模拟，海洋学报，1993，15(1)：107～118．
[38] 曹祖德，王运洪，水动力泥沙数值模拟，天津大学出版社，1994．
[39] 张存智，杨连武等，极浅海域潮流数值模型，海洋与湖沼，1994，25(5)：671～676．
[40] 韩康，张存智等，大连湾及附近海域潮流场数值模拟，海洋通报，1994，13(4)：20～25．
[41] 华秀菁，吕玉麟，二维浅水域潮流数值模拟的ADI-QUICK格式，水动力学研究与进展，A辑，1996，11(3)：77～92．
[42] Vincenzo Casulli. Semi-implicit finite difference methods for the two-dimensional shallow water equations. J. Comput. Phys.,1990,86:56～74.
[43] Vincenzo Casulli, Filippo Notamicola. An Eulerian-Lagrangian method for tidal current computation. In:Computer Modeling in Ocean Engineering (Sehrefler &.Zienkiewicz,eds),Balkema,Rotterdam, 1988.237～244.
[44] 温生辉，陈季良，厦门附近海域潮流场的数值模拟，海洋学报，1996，18(2)：15～25．
[45] 耿兆铨，二维非恒定流的有限元计算模式，华东水利学院学报，1981，(1)：1～12．
[46] Grotkop G.Finite element analysis of long period water waves.Comp. Meth.in Applied. Mech.and Eng., 1973.(2): 147～160.
[47] 吕玉麟，赖国璋，近海浅水环流问题的数值模拟，大连工学院学报，1981，20(1)：39～52．
[48] 谭维炎，赵枥华，二维浅水渐变不恒定明流的有限元算法及程序包，水利学报，1984，(10)：1～13．
[49] Gray W G. An efficient finite element scheme for two-dimensional surface computation.In: Proceedings of the 1st Intn. Conf. on Finite Elements in Water Resources,Pentech Press, 1977,433～449.
[50] 于天常，二维浅水环流问题的一个有限元模型，海洋与湖沼，1984，15(1)：37～45．
[51] 史金松，计算二维非恒定流的守恒型有限元法，水利学报，1987，(10)：41～48．
[52] 谭维炎，赵枥华，二维非恒定浅水明流的三种算法，水利学报，1985，(9)：1～9．
[53] 李浩麟，易家豪，河口浅水方程的隐式和显式有限元解法，水利水运科学研究，1983，(10)：15～26．
[54] 耿兆铨，二维非恒定流的显式迎流有限元模型，全国第一届水动力学学术会议论文集，北京：海洋出版社，1984．
[55] 张存智，杨连武等，具有潮滩移动边界的浅海环流有限元模型，海洋学报，1990，12(1)：1～13．
[56] 李浩麟等，河口二维潮流单元积分的数值计算，水利学报，1984，(11)：66～72．
[57] 郭庆超，何明民等，控制体积法在二维潮流计算中的应用，水动力学研究与进展，A辑，1995，10(6)：602～609．
[58] 张华庆，乐培九，钦州港港区流场的数值模拟，水利水运科学研究，1994，(1)：81～88．
[59] 董耀华，黄煜龄，河口潮流河段二维非恒定流数模研究及应用，长江科学院院报，1995，12(1)：31～39．
[60] 张东生，蒋勤，江苏北部灌河口悬沙输送数学模型，海洋学报，1991，13(1)：125～136．
[61] 金忠青，N-S方程的数值解及紊流模型，南京：河海大学出版社，1989．
[62] 林秉南，赵雪华等，河口建坝对毗邻海潮波影响的计算，水利学报，1980，(20)：16～26．
[63] 汪德灌，陈新建等，长江口南支流场的平面二维有限差分模拟，河海大学学报，1987，15(5)：34～46．
[64] 何少苓，林秉南，破开算子法在二维潮流计算中的应用，海洋学报，6(2)：260～271．
[65] 何少苓，龚振瀛等，式隐破开算子法在二维潮流计算中的应用，海洋学报，1985，7(2)：225～232．
[66] 吴江航，陈凯麒等，核电站冷却水远区热、核污染数值计算的一种新方法，水利学报，1986，(10)：16～25．
[67] 宋志尧，严以新等，二维潮流方程的一种分裂算子格式及算例，第七届全国海岸工程学术讨论会论文集，海洋出版社，1993，384～392．
[68] 周朦，方国洪，二维长波方程的一个无条件稳定有限差分格式，海洋与湖沼，1988，164～172．
[69] 吴江航，韩庆书，计算流体力学的理论、方法及应用，北京：科学出版社，1988．
[70] 张二骏，张东生等，一种新的数值方法——准分析法在二维非恒定流问题中的应用，海洋学报，1986，8(5)：636～643．
[71] 高飞，张二骏，局部一维准分析法的二维非恒定流的数学模型，水利水运科学研究，1986，(1)：65～75．
[72] 李家星，张镜潮，用潮波能谱法计算二维非恒定流，海洋工程，1988，6(2)：45～55．
[73] 张镜潮，李家星，连云港港区二维潮流计算，海洋工程，1991，10(1)：49～57．
[74] 华祖林，边界拟合坐标下流场计算模式，海洋工程，1994，12(3)：65～72．
[75] 汪德灌，杨艳艳，边界拟合坐标法的应用，河海大学学报，1989，17(1)：50～57．
[76] 王船海，程文辉，河道二维非恒定流计算，河海大学学报，1987，15(3)：39～53．
[77] Wei Wenli,Jin Zhongqing. Numerical solrtion for unsteady 2-D flow using the transformed shallow water quations. Journal of Hydrodynamics, Ser. B, 1995,7(3):65～71.
[78] 赵士清，长江口三维潮流的数值模拟，水利水运科学研究，1985，(1)：23～31．
[79] 张涤明，乔林等，近海环流三维样条分层有限元模式，水动力学研究与进展，1988，3(2)：83～93．
[80] Koutitas C.Three-dimensional models of coastal circulation:an engineering viewpoint. Coastal and Estuarine Sciences,4.American Geophysical Union,Washington,D.C., 1987.
[81] Koutitas C,BO"connor. Modeling 3-D wind-induced flows.J.Hydraulics Div.,ASCE, 1980, 106(11):472～485.
[82] 易家豪，叶雪祥，长江口南港航道三维水流数值模拟，第二次河流泥沙国际学术讨论会论文集，水利电力出版社，1983．
[83] Tee K T. The structure of three-dimensional tide generating currents. Estuarine, Coastal and ShelfScience, 1982,14(1):70～82.
[84] Heaps N S. On the numerical solution of the three-dimensional equations for tides and storm
surges.Mem. Soc. R.Sci.Liege., 1972,6(2): 143～180.
[85] Davies A M.Spectral models in continemtal shelf sea Three-dimensional Coastal oceanography. Ocean Models Norman S. Heaps,ed.,American Geophysical Union,Washington,D. C., 1987,71～106.
[86] 石磊，奚盘根，三维浅海域流体动力学方程的分步杂交解法，青岛海洋大学学报，1995，20(2)：162～172．
[87] 窦振兴，杨连武等，渤海三维潮流数值模拟，海洋学报，1993，15(5)：1～15．
[88] 奚盘根，冯士作等，非线性三维潮波边值问题模型，山东海洋学院学报，1988，18[2(1)]：37～47．
[89] 王晓建，张廷芳，非正交网格有限体积法在三维潮流场计算中的应用，水动力学研究与进展，A辑，1997，12(1)：56～61．
[90] 李孟国，伶仃洋三维流场数值模拟，水动力学研究与进展，A辑，1996，11(3)：342～351．
[91] 刘子龙，王船海等，长江口三维水流模拟，河海大学学报，1996，24(5)：108～110．
[92] Noye J,Stevean M.A three dimensional model of tidal propagatian using transformations and variable grids.Three-dirnen-sional Coastal Ocean Models,AGU, Washington,D.C., 1987.
[93] 李蓓，珠海电厂煤港海区海域二维波浪、潮流、泥沙淤积数学模型，第七届全国海岸工程学术讨论会论文集，海洋出版社，1993，149～158．
[94] 中华人民共和国交通部，海岸与河口潮流泥沙模拟技术规程(JTJ／T233-98)，人民交通出版社，1999．
[95] 李孟国，珠江口伶仃洋三维潮流数学模型，1993．
[96] 李孟国，泥灰礁三维潮流数学模型，2000．