中国近海潮波运动数值模拟
大连附近海域潮汐潮流的三维数值模拟
大连附近海域潮汐潮流的三维数值模拟
李庆红;赵军方;张永刚
【期刊名称】《海洋测绘》
【年(卷),期】2003(023)002
【摘要】在POM的基础上,建立了一个σ坐标系下的三维斜压预报模式,采用细网格,考虑海底摩擦系数的影响,以实测资料为基础,对大连附近海域的潮汐潮流进行了数值模拟.从单点潮位验证、单点潮流验证、流场的变化规律等多方面的模拟情况看,预报结果具有很高的精度,可以为海洋测量和航海等工作提供三维实时的动态变化规律.
【总页数】5页(P48-52)
【作者】李庆红;赵军方;张永刚
【作者单位】海军大连舰艇学院,海洋测绘系,辽宁,大连,116018;海军大连舰艇学院,海洋测绘系,辽宁,大连,116018;海军大连舰艇学院,海洋测绘系,辽宁,大连,116018【正文语种】中文
【中图分类】P229
【相关文献】
1.流沙湾海域潮汐潮流的三维数值模拟和海域环境容量研究 [J], 张静;刘玉;马志华;孙省利
2.大连湾及附近海域潮流场数值模拟 [J], 韩康;张存智
3.基于FVCOM的泉州湾海域三维潮汐与潮流数值模拟 [J], 林作梁;朱学明;鲍献文;刘钦政
4.基于FVCOM的獐子岛附近海域三维潮汐潮流数值模拟 [J], 齐继峰;曹圣山;郭可采;杨德周;徐振华;尹宝树
5.基于FVCOM的廉州湾及周边海域三维潮汐潮流数值模拟 [J], 宋德海; 鲍献文; 张少峰; 张春华
因版权原因，仅展示原文概要，查看原文内容请购买。
浙江近海潮汐潮流的数值模拟-海洋学报
第25卷第5期海洋学报Vol.25,N o.5 2003年9月ACTA OCEANOLOGICA SINICA September,2003浙江近海潮汐潮流的数值模拟陈倩1,2,3,黄大吉2,章本照1(11浙江大学力学系,浙江杭州310027;21国家海洋局第二海洋研究所浙江杭州310012;31国家海洋局海洋动力过程与卫星海洋学重点实验室,浙江杭州310012)摘要:用三维陆架海模式(HAM SOM)对浙江近海的潮汐、潮流进行了数值模拟,并采用网格嵌套和动边界技术对原模式作了改进,以提高计算的精度,改进后的模式在浙江近海的应用中被证明是成功的.沿岸50个潮位站计算与实测值的比较表明,加入动边界以后的小区域细网格计算较之粗网格以及未加动边界以前精度普遍提高,比较的均方差结果为:M2分潮振幅差416cm,相角差7114b;S2分潮振幅差510cm,相角差514b;K1分潮振幅差2125cm,相角差5176b;O1分潮振幅差1156cm,相角差515b,可见计算与实测符合良好.另外,选取了105个实测潮流点,比较了表层M2和K1分潮流调和常数分量U cos N,U sin N,V cos G,V sin G的实测值与计算值的偏差,结果表明计算与实测的符合程度较好.在此基础上,给出了各主要分潮的潮位同潮图、潮流同潮图、潮汐性质、潮流性质、潮流椭圆和潮流的运动形式等,发现4个主要分潮M2,S2,K1,O1在本区内均未出现无潮点;M2分潮流在29b18c N,122b46c E处有一个圆流点.此外还得到了一些有意义的结论,都与实测情况符合良好,从而对整个浙江沿海区域的潮汐潮流特性有了一个全面认识.关键词:数值模拟;网格嵌套;动边界;潮汐潮流;浙江近海中图分类号:P731123文献标识码:A文章编号:0253-4193(2003)05-0009-121引言浙江近海岸线曲折,地形复杂,港湾众多、岛屿星罗棋布.由北往南的主要港湾有杭州湾、宁波-舟山深水港、象山港、三门湾、乐清湾、温州湾等,均属强潮海区.沿岸和诸岛屿上设有一些验潮站和潮流测点.关于浙江沿岸各海区内潮汐潮流实测资料的研究已有不少成果,本文作者[1]也曾以多年来沿岸各潮位站观测资料以及海岸带和海岛调查的实测海流资料为依收稿日期:2002-09-25;修订日期:2002-12-10.基金项目:浙江省自然科学基金资助项目(499018);国家自然科学基金资助项目(40076010).作者简介:陈倩(1975)),女,浙江省舟山市人,国家海洋局第二研究所和浙江大学力学系联合培养博士,从事海洋动力过程研究.E-mail:dajih2001@10海洋学报25卷据,分析研究了整个浙江沿海区域的潮汐潮流特征,得到了一些规律性结论,但是实际观测既昂贵又费时,而且实测点有限,该海域的流场分布又较为复杂,因此要由点到面准确地掌握大面积潮位和潮流的分布变化规律,数值模拟计算仍不失为一种廉价而有效的研究手段.有关浙江近海各海区的数值模拟,前人已做过不少工作,曹德明等[2,3]用有限差分法对杭州湾的潮汐、潮流进行了二维数值模拟;李身铎等[4]采用R坐标下的三维数值模式模拟了杭州湾三维潮波运动;曹欣中等[5]做了宁波、舟山内海域的潮流场数值模拟计算;董礼先等[6]数值模拟了象山港水域的潮波运动;周大成[7]采用平面二维三角形单元显性有限元潮流模型对椒江河口大潮潮流特性进行了数值模拟;李孟国等[8]建立了时间二次插值的三角形网格显式差分数学模型,对瓯江口海区的潮流场进行了成功的模拟.但是,由于该区地形和岸线条件复杂,为提高数值计算的精度,绝大多数的计算工作是针对某局部海区进行的,极少把浙江近海作为一个整体进行数值计算,而且大多数的研究成果局限于某几点的潮位验证或某时刻的流场分布,未给出整个研究海域的潮分布规律(如各分潮的等振幅线和同潮时线等).也有许多模拟东海或渤海、黄海、东海的数值计算文章[9~11],但由于浙江近海只是其中很小一块区域,故其所提供的该区的潮汐潮流分布信息就不够细致.因此,把浙江近海作为一个整体,在较高精度下进行数值计算,并由此得到一系列规律性成果,是很有意义的,而且为了弄清潮汐潮流对物质输送和扩散以及浙江近海温度、盐度分布的影响,也需要把浙江近海作为一个整体来进行计算.本文采用三维陆架海模式对整个浙江海区的潮汐、潮流进行了数值模拟计算,并引入网格嵌套技术使得小区域内的细网格计算精度提高,引入动边界技术来处理露滩问题,使该区潮模拟更加真实可靠.本文用这两项技术对原有的HAMSOM模式进行了改进,改进后的模式被应用于浙江近海潮汐、潮流的三维数值计算中,经验证效果良好.在验证计算值与实测值符合良好的基础上,给出了各主要分潮的潮位等振幅线和同潮时线、潮汐性质、潮流等振幅线和同潮时线、潮流椭圆、潮流性质及潮流的运动形式.2数值模式及其应用211三维陆架海模式(HAMSOM)简介HAMSOM(H amburg shelf ocean model)是由德国汉堡大学海洋研究所Backhaus和他的同事们发展的三维斜压陆架海模式.自20世纪80年代初发展至今,模式已有了不少改进[12~15].下面对HAM SOM的框架及其特点作简要说明.HAMSOM是一垂向分层模式,控制方程建立在任一垂向层上.这是为简化数值计算,通过对原始三维运动方程和连续方程组进行层积分处理(即把相应的方程对模式的垂向某一层积分),得到层积分的连续方程、运动方程和温度、盐度方程[13],从而将三维问题转化为二维问题,层与层之间通过垂向动量交换联系起来.此外,针对限制时间步长的线性不稳定因子,采取相应的措施,使模式的时间步长不受稳定性限制:(1)科氏力项,在运动方程中通过对科氏力项引入一个稳定的二阶旋转矩阵来克服它在时间迭代过程中产生的线性不稳定[13];(2)外重力波项,对运动方程中的正压梯度力项和连续方程中的水平散度项,采用半隐式的差分格式,以克服由外重力波所引起的对稳定性的限制;(3)垂向黏性(扩散)项,对运动方程中的垂向黏性项和温度、盐度方程中的垂向扩散项,采用半隐式的差分格式,以克服该两项对稳定性的限制.通过以上处理并作进一步推导后,可以得到关于水位的椭圆型方程及输运量的三对角方程.用超松弛法可求解关于水位增量的椭圆型方程,从而得到下一时间步长的水位值,用追赶法可求解三对角方程得到水平输运量和水平流速分量,然后求解层积分方程,确定垂向流速场.由于本文仅探讨浙江近海的主要正压动力过程:潮汐和潮流,故采用HAM SOM 模式的正压模拟部分,对温度、盐度方程暂不考虑.212 模式的应用21211 计算海区和模式安排本计算海域为27b ~31b N,12011b ~12311b E 的整个浙江沿岸水域,区域地形如图1所示.经向和纬向水平网格均取1c ,纬线(x 方向)上网格的大小由纬度来校正,垂向取8层,这是为反映岛屿间水道或海峡通道地形的剧烈变化而设置.从上向下各层的厚度分别为10,20,30,40,50,75,100,150m.时间步长取5m in.陆地边界取法向流速为0;开边界潮位给定,包含5个主要分潮(M 2,S 2,N 2,K 1和O 1),由大区域模式(渤海、黄海、东海,如图2)提供.开边界处流速的边界条件根据Orlanski 辐射条件来确定.图1 区域地形(等值线单位:m)及实测比较点的位置图2 大、小计算海区及网格示意图在模式应用过程中,由于研究海域的特殊性,将引入两项数值技术,即网格嵌套和动边界技术.21212 引入网格嵌套技术为小区提供开边界条件浙江近海是个相对于渤海、黄海、东海等较小的海域,且沿岸地形复杂,岸线曲折.为了得到反映该区特征的数值计算结果,更好地刻划岸线以及尽可能反映影响因素(如江河入海径流),都希望把网格划分得足够细,但是小区域的开边界条件往往难以给定,因为研究区域小,115期 陈 倩等:浙江近海潮汐潮流的数值模拟开边界可能引入误差的影响就比较大,为减小此误差,又希望把计算海域尽量扩大.如果在扩大海域下采用单一的细网格,将会增加计算量,因此我们引入网格嵌套技术来解决这一矛盾,即对外部大区域用粗网格,对内部小区域用细网格分别进行计算,其中大区的计算结果为小区提供所需的开边界条件,这样既保证了小区域计算精度,又可避免开边界上过大的误差影响.本文取24b~41b N,117b~131b E的整个渤海、黄海、东海为大区域,经向和纬向水平网格都取10c,小区域即浙江沿岸海域的大小计算海区如图2所示.开边界上潮位边界条件为F=E5i=1H i cos(R i t-g i),(1)式中,i从1到5分别表示5个主要分潮M2,S2,N2,K1和O1.小区域计算中,开边界上的调和常数由大区的计算结果经线性插值得到[16~18].图3浙江近海的潮滩示意图21213引入动边界模拟潮滩涨落浙江沿岸河口、港湾众多,滩地广泛分布.图3为浙江近海的潮滩示意图,其中黑色区域为潮滩.潮滩沙脊等随着潮位的涨落时而淹没时露出,相应水域的面积也随之增大或减小.为正确模拟这类变化的水域,我们引入动边界技术.本文将采用水位判别法[17~23].在计算过程中,随时对滩地网格点作状态判别,根据该时刻某点的瞬时水深判断其是否淹没,以确定该点是否参与计算.对于所有的潮滩点,由于最大可能的瞬时水深都不超过10m(实测最大潮差为8193m),故垂向仅有一层,我们可以把它当作二维问题来讨论.对于平面上任一潮滩点(i,k),瞬时水深h i,k= D i,k+F i,k,其中D i,k,F i,k分别为(i,k)点的滩地高程及潮位.由于算法的限制,对动边界的判别指标往往不能以单元水深为0来确定,而需要引入一个富裕水深D0(本文中取D0=011m)来保证求解的稳定性.落潮过程中,F i,k<0,当h i,k[D0时,认为该点干出,不参与计算,且令其流速为0;反之,则认为该点淹没,参与计算.涨潮过程中,水边界线随潮位的上升向高滩推进,则计算网格点增多.因新增网格点原无潮位值,故取其周围4点(i+1,k),(i-1,k),(i,k+1),(i,k-1)中已淹水的诸点潮位的平均值:F i,k=(F i+1,k E i+1,k+F i-1,k E i-1,k+F i,k+1E i,k+1+F i,k-1E i,k-1)(E i+1,k+E i-1,k+E i,k+1+E i,k-1),(2)式中,E i,k=1水点,0陆点.涨潮时,F i,k>0,当h i,k\D0时,认为该点淹没,参与计算,反之不参与计算.12海洋学报25卷3 计算结果及与实测的比较根据潮汐调和分析理论,分离主要分潮调和常数所需的最短潮位时间序列为15d,在实际应用中,常用1个月的潮位时间序列以获得较好的结果.计算的初始条件为:当t =0时,u =v =w =F =0,其中u,v 分别为水平流速的东分量和北分量,w 为垂向流速,F 为未扰动海面上的潮高.潮汐和潮流模拟时,模式共运行了32d,在2d 内模式已由零初始状态完全建立起来,模式产生的水位和流速数据每隔1h 进行储存,后30d 的逐时数据用于分析计算潮汐和潮流的调和常数,这里我们用最小二乘法来计算7个主要分潮(M 2,S 2,N 2,K 1,O 1,M 4和MS 3)的调和常数.为了验证计算结果的可靠性,我们从历史观测资料中挑选出50个有代表性的沿岸水位站,其位置如图1中空心五角星所示.此外,选出105个潮流比较点,其位置如图1中实心圆点所示.通过比较模拟计算与实测点的调和常数来验证计算的可靠性.表1为50个潮位站上4个主要分潮调和常数之差值(该值为实测减模拟所得)的均方差比较.为验证网格嵌套和动边界技术应用的效果,我们分别列出了10c 固定边界、1c 固定边界、1c 动边界3种情况下的均方差.由表1可见,计算与实测都符合较好.M 2分潮振幅差值中有82%的站点1c 细网格的计算精度优于10c 粗网格,1c 动边界的计算精度又优于1c 固定边界;10%的站点1c 固定边界效果最好,1c 动边界其次,10c 粗网格最差;另外8%的站点1c 动边界计表1 50个潮位站上4个主要分潮调和常数的模拟与观测之差值的均方差比较M 2S 2K 1O 1$H /cm$g /(b )$H /cm $g /(b )$H /cm $g /(b )$H /cm $g /(b )10c812110160817010167310010100210781701c615081786120719221688100118071151c 动边界41607114510051402125517611565150算结果最好,10c 粗网格次之,1c 固定边界最差.其他分潮的振幅、相位差统计结果与此皆相差不大.因此,综合来看,应用了网格嵌套和动边界技术之后,各点的计算精度普遍得到提高,效果明显.在下面的讨论中,我们将仅以1c 动边界情况下所得的模拟计算结果作为研究对象.图4a,b 分别为50个潮位站上4个主要分潮(M 2,S 2,K 1,O 1)调和常数的振幅和相角的实测值(x 轴)与计算值(y 轴)的对比结果,图中的点均分布在第一象限的角平分线附近,表明计算值与观测值比较一致.浙江近海潮流的数值模拟虽然是在三维的HAMSOM 模式下进行,但除近海底外,潮流的垂向变化都较小,因此仅以表层的潮流特征来讨论.图5a,b 分别给出了105个测流站表层M 2,K 1分潮流调和常数分量U cos N ,U sin N ,V cos G ,V sin G 的实测值(x 轴)和计算值(y 轴)的对比结果.测点均匀分布于计算海区.由图可见,计算值与实测值符合较好.误差统计分析表明,对M 2分潮流来说,两者偏差绝对值小于10cm/s 者占93%,最大偏差为1617cm/s;对K 1分潮流来说,两者偏差绝对值小于3cm/s 者占97%,最大偏差为4159cm/s.135期 陈 倩等:浙江近海潮汐潮流的数值模拟图4 模拟与实测各主要分潮M 2,S 2,K 1,O 1调和常数之比较图5 模拟与实测分潮流分量U cos N ,U sin N ,V cos G ,V sin G 之比较a.M 2分潮流,b.K 1分潮流4 计算结果的讨论411 同潮时线和等振幅图图6a,b,c,d 分别为主要半日分潮M 2,S 2和主要全日分潮K 1,O 1的同潮时线和等振幅线图,由图可见,M 2,S 2分潮同潮时线的走向基本一致.在M 2同潮图中,250b 同潮时线以三门湾为顶点,呈八字形向两旁伸展.图中S 2的300b 同潮时线类似分布.同潮时线的这一分布特征表明:半日潮波进入陆架后,由东南向西北挺进,首先在三门湾口附近达到高潮,然后分南北两支传播.K 1,O 1分潮的等振幅线图大致相似,两者的振幅都是由东向西略有增大,但增幅很小.同潮时线分布表明,K 1分潮由东北向西南传播;O 1分潮以西北-东南向传入本海域,在三门湾附近传向渐变为东北-西南.在本海域内,各分潮都没有出现无潮点.14海洋学报 25卷图6潮位的同潮时线和等振幅线a1M2分潮, b.S2分潮, c.K1分潮, d.O1分潮等振幅线/cm;同潮时线/(b)155期陈倩等:浙江近海潮汐潮流的数值模拟412 潮汐类型潮汐类型是根据潮型数F =(H K 1+H Q 1)/H M 2来划分的,它反映了日分潮与半日分潮的相对重要性.根据值的大小,一般可把潮汐分为4种类型,即规则半日潮(010<F [015)、不规则半日潮(015<F [210)、不规则日潮(210<F [410)和规则日潮(F >410).由图7a 可见,浙江近海以半日分潮为主.大部分海域内F [015,为正规半日潮,如浙北的杭州湾以及浙中、浙南的沿岸海域;小部分地区015<F [210,为非正规半日潮,主要分布在甬江两侧并连及舟山群岛部分海区.其产生原因是由于潮波变形,H M 2减小,而H Q 1和H K 1基本不变,故比值增大.图7 潮型数(F )a.潮位,b.表层潮流413 潮流的同潮时线和等振幅图图8a,b 分别给出了M 2,K 1分潮流表层合成流等振幅线和同潮时线的分布.由图可见,M 2,K 1分潮流等振幅线的分布显示,在本海域内存在几个明显的强流区.浙北的杭州湾和舟山群岛诸水道内潮流最强,M 2分潮流的振幅可达120cm/s 以上,这是由于潮流受地形影响显著.由于杭州湾的喇叭口地形,定海、岱山、嵊泗等海区的一些狭长水道,潮流通道的截面较小,故潮流速很强.另外,在三门湾、温州湾等港湾区,潮流振幅随近岸距离的减小而有明显增强趋势,这种变化趋势主要是由海区的水深条件、岸线和地形等因素造成的.M 2分潮流在29b 18c N,122b 46c E 处有一个圆流点,这与文献[11]中所得结果一致.同潮时线的旋转方向为反时针旋转.16海洋学报 25卷图8 潮流的同潮时线和等振幅线(表层)a 1M 2分潮, b.K 1分潮等振幅线/cm,等潮时线/(b )图9 M 2分潮流旋转率K 的分布(表层)414 潮流类型和潮流的运动形式潮流类型的划分标准与潮汐类型相类似.根据我国5港口工程技术规范6[24]的规定,采用(W O 1+W K 1)/W M 2作为指标,其中W O 1,W K 1,W M 2分别为O 1,K 1,M 2分潮流的椭圆长轴.由图7b 可见,绝大部分海域表层潮流F 值在012左右,均小于015,故浙江近海基本都属于正规半日潮流区.潮流的运动形式由潮流的椭圆旋转率K 来描述.K 值为潮流椭圆的短轴与长轴之比.当K 值大于0125时,潮流表现出较强的旋转性,而当其小于0125时,潮流表现为往复流.K 值前面的正负号表示潮流的旋转方向,正号为左旋,负号为右旋.由于本区域内半日潮流具有支配地位,因此我们给出M 2分潮流旋转率K 的分布来表征潮流的旋转特征,如图9所示.由图可见,在沿岸的港湾、河口水域及潮汐通道等处,潮流运动由于受地形、边界条件的制约,往复流的性质非常明显,K 的绝对值多小于012.外海或离岸较远且较宽敞的海区K 175期 陈 倩等:浙江近海潮汐潮流的数值模拟的绝对值大于014,呈旋转流态,如大目洋、三门湾、温州湾外海等.从旋转方向上看,大多数站点的K 为负号,呈右旋,这是由于地球自转效应产生的结果.另外,杭州湾口、浙闽交界水域各有一明显左旋区,这是潮波干涉区的影响所致.415 潮流椭圆图10a,b 分别给出了5c 为间隔的计算点上主要半日分潮M 2和主要全日分潮K 1的潮流椭圆长短轴分布.潮流的椭圆长轴指示了最大流速和最大流速方向.由图可见M 2,K 1分潮流椭圆长短轴的分布较为相似.从南往北,流速渐增.远岸区域旋转性较强,流速较弱;近岸区域则多为往复流,流速较强.M 2分潮流椭圆长轴的分布比较规则,其变化与地形密切相关.本区东南海域椭圆长轴的走向为西北-东南,这显示了潮波的传入方向.潮波传入后分为两支,往北和往南的传播方向各不相同,这在椭圆长轴的走向上有所反映.往南的一支表现为椭圆长轴在浙南海区多为东-西走向,而往北一支潮波则表现为椭圆长轴由浙中的东北-西南走向渐变为浙北的西北-东南走向,进入杭州湾后,基本为东-西走向.在港湾区或河汊水道处,椭圆长轴的方向一般与岸线或港湾水道走向相一致.K 1分潮流椭圆长轴之走向,在三门湾以北区域多为西北-东南,从三门湾往南渐变为东北-西南走向.在近岸区域,椭圆长轴的方向与岸线平行.图10 潮流椭圆长短轴分布图(表层)a.M 2分潮流,b.K 1分潮流5 小结本文用三维陆架海模式对整个浙江近海的潮汐和潮流进行了三维数值模拟,并针对浙江近海岸线曲折、潮滩广泛分布等地形特点,引入网格嵌套和动边界技术以提高计算的精度.改18海洋学报 25卷进后的HAM SOM 模式被成功地应用于该区的潮流数值模拟中,通过与实测点调和常数的比较,验证了这两项技术引入原模式后的良好效果,计算精度普遍得到提高.这是HAMSOM 模式首次应用于浙江近海的潮汐、潮流计算,该模式在研究陆架海动力学上有着独特的优越性.模式在应用过程中根据实际需要所作的改进也被证明是成功的,这使得该模式得到进一步的完善.改进后的模式可以嵌套用于局部更小区域的精细计算,也可以处理边界变动的问题.除模式的改善和成功应用外,本文将浙江近海作为一个整体进行数值计算并得到了一批反映规律的研究成果.本文给出了全区范围内4个主要分潮M 2,S 2,K 1,Q 1的潮汐同潮图,发现在本海域内,这4个分潮都没有出现无潮点;给出了M 2,K 1分潮潮流同潮图,发现M 2分潮流在29b 18c N,122b 46c E 处有一个圆流点,K 1分潮流在该区内无圆流点;此外还得到了潮型数F ,M 2分潮流椭圆率K 和潮流椭圆的分布规律.以上模拟结果都与实测情况[1]符合良好,而模拟结果由于不受实测点的限制,因此更全面、完整.参考文献:[1] 陈 倩,黄大吉,章本照.浙江近海潮汐特征的研究[J].东海海洋,2003,21(2):1)12.[2] 曹德明,方国洪.杭州湾潮汐潮流的数值计算[J].海洋和湖沼,1986,17(2):93)101.[3] 曹德明,方国洪.杭州湾和钱塘江潮波的联合数值模型[J].海洋学报,1988,10(5):521)530.[4] 李身铎,顾思美.杭州湾潮波三维数值模拟[J].海洋与湖沼,1993,24(1):7)15.[5] 曹欣中,唐龙妹,张月秀.宁波、舟山内海域实测海流分析及潮流场的数值模拟[J].东海海洋,1996,14(2):1)9.[6] 董礼先,苏纪兰.象山港潮波响应和变形研究.Ò.象山港潮波数值研究[J].海洋学报,1999,21(2):1)8.[7] 周大成.椒江河口大潮潮流特性的数值分析[J].浙江水利水电专科学校学报,2001,13(3):13)14,26.[8] 李孟国,王正林.瓯江口潮流数值模拟[J].长江科学院院报,2002,19(2):19)22.[9] 赵保仁,方国洪,曹德明.渤、黄、东海潮汐潮流的数值模拟[J].海洋学报,1994,16(5):1)10.[10] 万振文,乔方利,袁业立.渤、黄、东海三维潮波运动数值模拟[J].海洋与湖沼,1998,29(6):611)616.[11] 王凯,方国洪,冯士.渤海、黄海、东海M 2潮汐潮流的三维数值模拟[J].海洋学报,1999,21(4):1)13.[12] BACKHAUS J O.A sem-i implicit scheme for the shallow water equations for applications to shelf sea modelling[J ].ContS helf Res,1983,(2):243)254.[13] BACKHAUS J O.A three -di m ensional model for the simulation of shelf sea dynamics [J ].Deutsche HydrographischeZeitschrift,1985,38:165)178.[14] BACKHAUS J O,HAINBUCHER D.A finite difference general ci rculation model for shelf seas and i ts application to lowfrequency variability on th e North European shelf[A].NIHOUL J C ,JAM ART B M .Three -Dimensional M odels of M arine and Estuarine Dynamics[M ].Amsterdam:Elsevier Science Publishing B.V.,1987.221)244.[15] 黄大吉,陈宗镛,苏纪兰.三维陆架海模式在渤海中的应用[J].海洋学报,1996,18(5):1)13.[16] 辛文杰.差分模型网格嵌套边界技术在工程潮流计算中的应用[J].水利水运科学研究,1999,12:355)360.[17] 辛文杰.河口、海湾平面潮流数值计算中的几个问题[J].水动力学研究与进展,1993,8(3):348)354.[18] FULT ON S R.An adaptive multigrid barotropic tropical cyclone track model[J].M onthly Weather Review ,2001,129(1):138)151.[19] 曹德明.胶州湾潮汐潮流的数值计算[J].海洋科学集刊(21)[C].北京:科学出版社,1984.157)164.[20] 孙英兰,张越美.胶州湾三维变动边界的潮流数值模拟[J].海洋与湖沼,2001,32(4):355)361.[21] 李燕初,蔡文理.ADI 潮汐模型的活动边界方法及其效应[J].海洋学报,1993,15(2):115)120.[22] 张存智,杨连武,窦振兴.具有潮滩移动边界的浅海环流有限元模型[J].海洋学报,1990,12(1):1)13.[23] 韩 康,吴 冠,张存智.普兰店湾潮流场数值模拟[J].海洋环境科学,2001,20(1):42)46.[24] 中华人民共和国交通部.港口工程技术规范(上册)[M ].北京:人民交通出版社,1987.195期 陈 倩等:浙江近海潮汐潮流的数值模拟20海洋学报25卷Nu merical simulation of tide and tidal currentsin the seas adjacent to ZhejiangCH EN Qian1,2,3,HUANG Da-ji2,ZHANG Ben-zhao1(1.Department of M echa nics,Zhe j iang Univ iersity,Hangz hou310012,China;2.Second I nstitute of Oceanography,State Oceanic A d ministration,Hangz hou310012,China;3.K ey L aboratory of Ocean Dynamic Proc esses and Satellite Oceanogra-p hy of S tate Oce a nic A d ministration,Hangz hou310012,China)Abstract:By means of three-dimensional baroclinic primitive equation model)))Hamburg shelf o cean model (HA M SOM),the tide and tidal cur rents in the seas adjacent to Zhejiang are simulated.Fur thermore,the original model is impro ved by two numerical technologies)))nested g rid and mov ing boundary method,which are intro-duced to increase the computat ional precision.T he impro ved mode is proved to be successful while it is applied to the seas adjacent to paring the computed values w ith t hose of50tidal observ ator ies,it is found that the computat ional precision w ith fine gr ids and moving boundar y is gener ally higher t han that with coarse gr ids or fix ed boundary.T he root-mean-square values of comparativ e results show that t he difference betw een the simulated and t he observed amplitudes of M2constituent is only4.6cm,the differ ence of phase-lags is7.14b;the difference of amplitudes and phase-lags of S2constituent are5.0cm and5.4b;the difference of amplitudes and phase-lags of K1 constituent ar e2.25cm and5.76b;the difference of amplitudes and phase-lags of O1constituent are1.56cm and 5.5b.T hese indicate that the computational r esults agree with the observed ones very well.I n addition,105current stat ions are chosen,and the difference between the calculated and the observed harmonic co nstant,U cos N,U sin N, V cos G,V sin G of M2and K1component curr ents at surface layer is compared,and the results also show a goo d a-gr eement.Based on these results,the co-amplitude and co-phase lag lines of the main co mponent tides,the type of t ide,tidal current ellipse,the type and the mov ing mode of tidal current ar e given.It is found that the four main constituents M2,S2,K1,O1have no tide-fr ee points in this ar ea;and the M2co mponent current has o ne current-amphidromic point at29b18c N,122b46c E.In addition,some meaning ful results are concluded,and ag ree w ith the observed ones w ell.T hereby,a thorough kno wledg e of the character istics of tides and t idal currents is got in the w ho le coastal zone of Z hejiang Prov ince.Key words:numerical simulation;nested gr id;mo ving boundary;tides and tidal cur rents;seas adjacent to Zhejiang。
长江口杭州湾及邻近海区潮汐潮流场三维数值模拟
文章编号:1000-5641(2001)03-0074-11长江口杭州湾及邻近海区潮汐潮流场三维数值模拟杨陇慧, 朱建荣, 朱首贤(华东师范大学河口海岸国家重点实验室,上海 200062)摘要:把长江河口、杭州湾及邻近海区作为整体,应用三维高分辨率非正交曲线网格河口海洋模式,模拟了4个主要分潮M 2,S 2,K 1,O 1。
在长江口外半日分潮M 2、S 2从东南方向传入长江口和杭州湾,全日分潮K 1,O 1从北向南传播。
这4个分潮的振幅在长江口南支向上游逐渐减小,但因杭州湾和长江口北支呈喇叭状,而向上游逐渐增加。
计算潮差变化过程和实测值基本一致,4个分潮潮位振幅和位相的计算值与验潮站观测值相比,误差大部分在10%以内。
结合1996年2—3月长江河口现场观测,考虑了径流的作用,三维数值模拟了计算域内流场。
结果表明,即使在斜压效应不太明显的口门内,流速在垂直方向存在着明显的差异,上层流速明显大于下层流速;潮流具有不对称性,由于径流的作用,落潮时间明显大于涨潮时间,落潮流大于涨潮流,但在象北支涨潮槽中,涨潮流反而比落潮流大。
模拟出的以上结论与观测结论较为一致。
关键词:长江河口; 杭州湾; 潮汐潮流; 三维数值模拟中图分类号:P736.211 文献标识码:A0 引 言长江是我国第一大河,河口地形复杂,三级分汊四口入海(南支和北支,南槽和北槽,南港和北港),口门附近存在水下沙坝(图1)。
长江河口的动力过程也十分复杂,影响因子有径流、潮流、波浪、密度流等,其中潮汐潮流在这个系统中起着重要作用[1～4]。
长江河口是中等强度的潮汐河口,口外是正规半日潮,口内是非正规半日潮。
南支的潮差范围从口门到上游逐渐减小。
口门附近的中浚测站的多年平均潮差为2.66米,最大值是4.62米。
黄浦江附近的吴淞站潮差是2.21米。
由于北支的河槽呈喇叭状,所以其潮差大于南支,并且由口门向上游逐渐增大。
潮流在口门内作往复运动,在口门外为旋转流。
中国近海潮汐数值模型及可视化应用研究的开题报告
中国近海潮汐数值模型及可视化应用研究的开题报告开题报告一、选题背景：潮汐是海洋中的一种规则的周期性运动。
由于潮汐对于海洋与海岸地区的生态系统、天气和气候、海洋交通运输、海洋工程等方面均具有重要的影响，因此研究潮汐运动及其规律对于海洋学和海洋经济的发展至关重要。
中国近海作为全球最大的沿海国家，近海潮汐数值模型及其可视化应用研究具有重要意义。
本研究将通过对中国近海潮汐现状、数值模型及可视化应用的分析，对中国近海潮汐运动及规律进行深入探究，为海洋相关领域的研究和应用提供科学依据。
二、研究目的和意义：本文旨在开展对中国近海潮汐数值模型及可视化应用的研究，具体目的包括：（1）了解中国近海潮汐现状；（2）研究中国近海潮汐数值模型，探究潮汐运动规律；（3）开发潮汐数值模型的可视化应用，方便用户快速准确地了解潮汐情况及变化。
本研究的意义在于：（1）为海洋相关领域的研究提供科学依据；（2）为海洋资源开发、海域利用、海洋环境保护等领域的决策提供科学支持；（3）推进中国近海潮汐数值模型及可视化应用领域的发展。
三、研究内容：本文将主要开展以下研究内容：1. 中国近海潮汐现状分析，包括中国沿海主要港口的潮汐现象、潮周期、潮高、潮差等内容，并绘制潮汐图、潮位表以及潮流图等。
2. 中国近海潮汐数值模型研究，综合采用数值模拟、观测数据等方法，建立中国近海潮汐数值模型，并对模型进行验证和丰富实验。
3. 潮汐数值模型可视化应用开发，通过对潮汐数值模型的可视化处理，开发出一款潮汐数值模型可视化软件，方便用户快速准确地了解潮汐情况及变化。
四、研究方法：1. 文献资料法：通过对近年来与中国近海潮汐数值模型及其可视化应用相关的文献、报告、资料等进行收集、整理和分析，进一步掌握国内外有关研究现状和发展趋势。
2. 数值模拟法：采用数值模拟方法，建立中国近海潮汐数值模型，研究潮汐运动规律。
3. 统计分析法：对中国沿海主要港口的潮汐数据进行统计分析，查明潮汐的周期、潮高、潮差等特征。
南海北部内潮与非线性内波：观测与数值模拟研究
南海北部内潮与非线性内波：观测与数值模拟研究基于水文观测和数值模拟研究了南海北部内潮和非线性内波现象，主要内容有：吕宋海峡的测流观测与线性内波配极关系的比较；南海北部非线性内波背景环境及运动学参数的地理特征及其季节变化；中尺度涡旋对内孤立波传播影响的数值模拟；南海北部陆坡区内孤立波极性转换的数值模拟和吕宋海峡双屋脊上内波产生的数值模拟。
1.吕宋海峡内潮的观测分析鉴于吕宋海峡观测流的稀少，开展了该海域的海流锚定观测。
运用统计和经验正交函数分析方法研究了观测流特征。
谱分析和能量计算均显示全日潮和半日潮为主要的能量频带，近惯性频率峰值仅出现在斜压分量。
在66m以浅，显著分潮和近惯性频率分量都基本符合线性内波能量一致性关系式E+(ω)/E-(ω)=(ω-f)~2/(ω+f)~2；在其它观测深度，并非所有的主要分潮都满足上述关系。
正压潮和内潮均为混合潮，潮日不等现象明显。
全日分潮平行与垂直于陆架方向的分量几乎相等；而半日分潮垂直于陆架方向分量远大于平行于陆架方向分量；M2和S2潮垂向结构主要表现为第一模态，而K1和O1潮则接近于第二模态；在主温跃层附近，它们的短轴与长轴之比与比值f/ω接近。
2.南海北部非线性内波背景环境及运动学参数的地理特征及其季节变化从遥感图像的统计分析表明，南海北部内波呈现显著的季节变化。
基于海洋再分析资料研究了南海北部的层化特征和非线性内波运动学参数的地理分布特征和季节变化。
南海北部季节性密度跃层从2月开始出现，最大浮力频率约在20m；它在6-7月达到最强，自8月开始减弱，在10月消退。
在8-11月出现另一较深的密度跃层，最大浮力频率约在80m，冬季大致在120m。
季节性密度跃层在4-9月十分明显，在8-10月出现双跃层现象，而在冬季仅出现较弱的第二密度跃层。
在1-3月和10-12月深水区最大浮力频率值要大于浅水区，而在5-9月情况则相反。
浮力频率最大值所在深度随季节变化显著，冬季最深；6-7月则最浅。
渤、黄、东海海域9711号风暴潮数值模拟
一
／ ， ／
ｖ８ ｕ
＋ ＋…
—
ＵＶ
＋ ＋一
ａ ，
ｃｓ ａ ｏ妒
ａ 矽
一 喀 ＝ ， 一 ｆ ＝ ｖ一一 ｇ
望
ｕ
ｐａｃ ｓｐ ｏｒ
￣． ２ 一 ｃ ｓ ａ ｏ矽
一一
ｇ ０ （
一
（２）
＋
［０ （券＋一 ＿ ｏ ］ １ｕ ｃ ］ － ｓ ２ ＋
中 图分类 号 ：Ｐ ３ ７２ 文献标 识码 ：Ａ 文章 编号 ： １０—９２２ ０）６０ １—０ ５ ０ １６３ （０ ８０—０５０ ０
风 暴潮 是 由于 台风 周 围的风 推动 产 生 的，风暴 潮 在大 洋或 开 阔海面 上 引起 的水面 异 常升 高一 般不 大 ， 当耦合 天文 潮 向岸行 进 到近 岸 时 由于 地形 的影响 ，这 种水 面异 常升 高往 往会 急剧 增加 ，引起 灾害 发生 。 早期 的风暴 潮数 值 模拟 为纯 风暴 潮模 型 ，往 往采用 辐射 边 界条 件而 不考 虑天 文潮 与风 暴潮 的 非线 性相 互 作用 关 系 。 景庸 等…于 ２ 世 纪 ８ 汪 ０ ０年 代末 ，以 Ｍ２ 潮作 外海 边 界开始 考 虑天 文潮 对风 暴潮 的影 响 。 分 ９ 年代 末 ，周 旭波 和 孙文 心 【采用 Ｍ２ ２ ０ ２ Ｊ 、Ｓ、Ｋ１ Ｏ１ ４ 个 天 文分 潮作 为 外海 开边 界研 究 长江 口风 暴 和 共 潮 增水 ，指 出 台风 过 境 时所 造 成 的风暴 潮增 水过 程 中含有 潮 波成 分 ，只 有采 用考 虑 了风暴 潮与 天文 潮 的非 线 性相 互作 用 的模式 才 能更精 确 地模拟 出风暴 潮增 水过 程 。 风 暴潮 的 发生海 域 一般 是跨 多个 经纬 度 的 ，其 数值模 拟 计算 域 的球面特 征 明显 ，但 国 内采 用球 面坐 标 系研 究 风暴 潮 的模 型并 不多 见 ，较 多 的是在 笛卡 尔坐 标系 基础 上采 用 纬度修 正 。本文 建立 了以整个 东 中 国 海 （ ３７ ２ ．。Ｎ一４ ．。Ｎ； １７ 。Ｅ一１ １１ １ ２ １． ５ ３ ．。Ｅ）为计算 域 ，采用 ６ 个 天文 分潮 为控 制边 界 的球面 坐标 系 ４
南海西北部陆架潮频内波生成过程的数值模拟
南海西北部陆架潮频内波生成过程的数值模拟我国南海海水运动复杂、海底地形多变,是内波的多发区域。
目前关于南海内波的报道、观测和研究主要集中在Luzon海峡到南海东北部东沙群岛附近海域。
2005年4月16日始,中国科学院海洋研究所和中国海洋石油总公司联合在海南
岛东部的文昌站进行了为期近半年的连续定点海洋调查,观测资料清楚地表明,
这一海域陆架上长期存在着比较规则的潮频内波现象。
从航天图片和SAR卫星图片上捕捉到的信息来看,南海西北部陆架上(文昌
站附近海域)有些潮频内波并非从Luzon海峡传播而来,而极有可能是局地正压
潮与陆坡作用生成的。
本文用数值模拟的方法来证明这一观点。
首先,我们从Reynolds方程组出发,用差分方法和有限元方法建立了一个连续层结海水中正
压潮与海底地形相互作用生成潮频内波的二维非线性数值模式。
我们用该模式研究了正压潮流流经下降地形、上升地形、凸起地形、凹陷地形等几种典型的海底地形时生成内波的过程,探讨了该数值模式中各个参数对生成内波的影响。
然后,我们从文昌站海洋科学调查的温盐深数据资料中提取了该海区的密度、跃层厚度和强度等海洋环境参数;用调和分析法从海流流速、流向数据中提取了正压潮潮流成分及各潮流分量的速度、振幅、方向等与数值模式相关的信息。
我们把这些从实测资料中提取出的海洋环境和潮流信息应用到内波生成的非线性数值模式中,模拟了文昌站附近海域正压潮流与陆坡作用生成内波的过程。
模拟结果表明,南海局地潮流与陆坡作用足以生成能够被SAR捕捉到的大振幅内波。
从而我们可以断定,南海西北部陆架上有些潮频内波是天文潮与陆坡作用生成,然后传播到陆架上的。
广西近岸海域潮流数值模拟
广西近岸海域潮流数值模拟
广西近岸海域潮流数值模拟
采用有限元三角形网格的分步杂交方法,建立了广西近岸海域的二维潮流数值模型,计算值与实测资料符合较好.采用主要分潮组合输入,模拟了研究海域的平均潮潮流场.模拟结果表明:涨急时,潮流向为东北方向,最大涨潮流速为74 cm/s左右:落急时,潮流向为西南方向,最大落潮流速约100 cm/s,落潮流速大于涨潮流速.近岸区域潮流为往复流,离岸边越远潮流越接近旋转流.
作者：张燕孙英兰张学庆刘晓丹ZHANG Yan SUN Yinglan ZHANG Xueqing LIU Xiaodan 作者单位：中国海洋大学环境科学与工程学院,山东,青岛,266003 刊名：海洋通报ISTIC PKU英文刊名：MARINE SCIENCE BULLETIN 年，卷(期)：2007 26(5) 分类号：P731.2 关键词：潮流数值模拟广西近海有限元。
基于模式嵌套的近海海浪数值模拟试验
（ 徽 省 巢 湖 水 文水 资源 局 ， 湖 ２ ８ ０ ） 安 巢 ３ ０ ０
摘 要 ： 目前 采 用 的几 种 能 量 平 衡 方 程 的 第 三 代 海 浪 预 报 模 式 中 ， ＡＶＥ ＡＴ Ｈ一 主 要 用 Ｗ Ｗ Ｃ Ⅲ
于分 辨 力低 的计 算 网格 下进 行 大范 围海 区的模 拟 ， ＷＡＮ 模 式则 能 对 目标 海域 进行 较 为 细致 Ｓ
Ｗ ａ ｇ Ｙａ ｗｕ ｎ ｏ
（ ｕｈｕｉＣｈａ ｈｕ Ｈ ｙｄｒ ｌ ｇｙ ｄ ａｔｒ Ｒｅ ｏ ｃ ｓ Ｂｕｒ ａｕ， Ａ ｏ ｏｏ ａｎ Ｗ ｅ ｓ ｕｒ ｅ ｅ Ｃｈａ ｈｕ ２３ ００ ｏ ８ ０）
Ａ ｓｒ ｃ ： ｒｓ ｎ ， ＡＶＥＷ ＡＴＣＨ一 ｉ ｍｏ ｇ ｓ ｖ ｒ ｌｔ ｉｄ ｇ ｎ ｒ ｔｏ ｖ ｒ ｄｃｉｎ ｍｏ ｅｓ ｂ ｔａ ｔ Ａｔ ｐ ｅ ｅ ｔ Ｗ Ⅲ ｓ ａ ｎ ｅ ｅａ ｈ ｒ — ｅ ｅａ ｉｎ ｗａ ｅ ｐ ｅ ｉｔ ｄ ｌ ｏ
第 １ 期
２１ ０ １年 ３月
气 象 水 文 海 洋 仪 器
Ｍ ｅ ｅ ｏｌ ｉａ ， ｙ ｏｌ ｉａ ｎ ａ ｉ ｎ ｔ ｕ ｅ ｓ ｔｏｒ ｏｇｃ ｌ Ｈ ｄｒ ｏｇ ｃ ｌａ ｄ Ｍ ｒｎｅ Ｉ ｓ ｒ ｍ ｎｔ
Ｎｏ １ ． Ｍ ａ ． ０１ ｒ２ １
基 于模 式嵌 套 的近 海海 浪数 值 模 拟试 验
０ 引 言
目前采用 能量 平衡 方程 的第 三 代海 浪 预报 模
式 ，主 要 以 Ｗ ＡＭ 【 、ＷＡ Ｅ Ａ Ｃ 一 ｌ 和 １ ］ Ｖ Ｗ Ｔ Ｈ Ⅲ ＿ ２ ］
Ｓ ＷＡＮ 模 式 采 用 基 于 能 量 守 恒 原 理 的 平 衡
方程， 除了考 虑第 三代 海浪模 式 共有 的特 点 ， 充 还 分 考虑 了模 式在 浅 水 时 的 各种 物 理 过 程 ， 在 进 但 行计 算 时效 率不 高 。本 文 在进 行舟 山海 域高分 辨
乐清湾潮波变形特点及其数值模拟分析
乐清湾潮波变形特点及其数值模拟分析利用乐清湾8个潮位站半个月的观测资料及11条垂线两个航次的实测海流资料，采用调和分析、准调和分析方法，从潮差、涨落潮历时的变化，M<sub>2</sub>、M<sub>4</sub>、M<sub>6</sub>沿湾振幅、迟角的变化等方面分析了乐清湾的潮汐特征；以及从潮流历时、涨落潮平均流速变化，分潮流M<sub>2</sub>、M<sub>4</sub>、M<sub>6</sub>北方向、东方向上的调和常数，潮流运动性质、运动形式、余流的变化等方面分析了乐清湾的潮流特征。
从潮汐的变化来看，潮差沿湾逐渐增加，平均潮差在大门岛为4.13米，到东山头增加到4.60米。
涨潮历时稍大于落潮历时，且向湾顶涨落潮历时差逐渐增大，涨落潮历时差从洞头的4分钟，到东山头增大为15分钟。
M<sub>2</sub>分潮振幅从湾口到湾顶缓慢增加，从湾口的2.米到湾顶增加到2.2米；M<sub>4</sub>振幅增加了近10倍，湾口大门岛M<sub>4</sub>振幅不到1厘米，湾顶的东山头M<sub>4</sub>振幅增加到9厘米，并从乐清湾的地学形态为湾内存在M<sub>4</sub>分潮的共振提供了证据。
对潮位的FFT（傅立叶分析）的结果发现，湾内涨落潮历时不等和非线性而产生的倍潮有关。
倍潮中1／4日分潮的作用较大，湾内的1／4日分潮是造成了涨、落潮历时的变化和涨落流速的不等主要原因，236≤2g<sub>M2</sub>-g<sub>M4</sub>≤328说明湾内一直是涨潮历时大于落潮历时。
宁波舟山近海三维潮汐潮流数值模拟
宁波舟山近海三维潮汐潮流数值模拟熊伟;刘必劲;孙昭晨;梁书秀;张亦飞【摘要】利用FVCOM模型对宁波舟山近海的潮汐潮流进行了三维数值模拟,并对其水动力特性作了相应分析.FVCOM模型采用非结构化三角形网格,很好地解决了精确拟合宁波舟山群岛的复杂岸线的问题.基于计算区域内的海床性质,采用Koutitas公式对FVCOM模型的中海底摩阻系数的计算进行了改进.通过计算域内多个潮位站和海流站的实测资料验证表明,改进的摩阻系数计算公式是合理的,流场的计算结果与实测符合良好,可以用于三维污染物扩散和泥沙输移计算.%By means of FVCOM model,the tide and tidal currents in the seas adjacent to Ningbo and Zhoushan were simulated.Furthermore,the dynamic characteristics of tide were analyzed.The grid system of FVCOM model was unstructured triangle mesh,which solved the accurate fitting of complex of Ningbo and Zhoushan shoreline.Based on the property of seabed in the calculation zone, Koutitas formula was used to instead of the original calculation method of drag coefficient in paring the computed values with those tidal observatories, the improved formula proved to be reasonable and the two values are in good agreement.The computed model can be used for 3D pollutant concentration and sediment transport calculation.【期刊名称】《水道港口》【年(卷),期】2011(032)006【总页数】9页(P399-407)【关键词】摩阻系数;潮汐;潮流;数值模拟;FVCOM模型【作者】熊伟;刘必劲;孙昭晨;梁书秀;张亦飞【作者单位】大连理工大学海岸和近海工程国家重点实验室,大连116023;大连理工大学海岸和近海工程国家重点实验室,大连116023;国家海洋局杭州海洋工程勘测设计研究中心,杭州310012;大连理工大学海岸和近海工程国家重点实验室,大连116023;大连理工大学海岸和近海工程国家重点实验室,大连116023;国家海洋局杭州海洋工程勘测设计研究中心,杭州310012【正文语种】中文【中图分类】P731.23;O242.1Biography：XIONG Wei（1988-），male，master student.宁波舟山及其近海水域岛屿密布，漕滩相间，深水航道众多，地形极其复杂。
基于潮波数值模拟技术的潮位推算及其精度验证
线、 海底摩擦 、 地转偏向力以及浅水效应等因素影 响 ， 动力变化 巨大 ， 潮汐性质复杂 ， 潮差空间分布悬殊 ， 因 此在水深测量 的同时必须进行潮位观测 ， 以提供必要 的水位改正资料 。 传统的水位改正方法是在岸上或海 上设立验潮站 ， 实行双站分带或三（ 站分区水位改正 ， 多） 有效降低 了潮差对测量结果的影响 ， 提高了测深精 度， 但也存在如下问题 ：１海上及孤岛验潮站设立成本较高 ， （） 且无精度及完整性保障 ；２要求各验潮站 的 （） 潮位 数 据为水 深 测量 期 间 的 同步 数 据 ， 能充 分 利用 测 区 内已有 的潮 位 资料 ；３ 原则 上要 求测 量 区域 内潮 不 （） 时、 潮差在 验潮 站之 间或 控制 范 围 内线性 变 化 。
为交点因子与交点订正角 ， 为天文初相角 ， 为分潮角速率 ； ￡ 。 乱（） 为余水位。
２ 基于潮波数值 模拟技术的潮位推算算法基本流程
（ ） 于尺 度不 大 的测 区而言 ， 以通 过调 试 和运行 一个二 维潮 波数 值模 式 ， 算测 区 内网格 点 主要分 １对 可 计
潮 的调 和常数 。
１６ ９
水
道
港 口
第 ３ 卷第 ３ １ 期
＋
＋
：０
ａ ｙ
ａ ａ ａ 一 （ ＋ 一 ｆ ， Ｚ 鲁 ）ａ ｖ ａ ａ ｇ ｘ 。 慧
＋
（ ２ ）
＋ ＝一 ｕ
ｃ＋ 鲁 熹一 ｇ 等
式 中 ：＿＋ ， ｚ ｈ为从平 均海 面起算 的水 深 ， ‘为潮位 ，、 别 为在 方 向和 Ｙ方 向的垂 向平均 流速 分量 为 ｕ 分 科 氏参量 ；和 Ａ分别为拖 曳系数 和水平涡动黏性 系量 ； 为重力加速度 。 ｒ ｇ 数值积分采用 已经成熟 的方法 ，Ｉ Ｄ ￣Ａ Ｉ ］
浙江近海海域潮波数值模拟的开题报告
浙江近海海域潮波数值模拟的开题报告一、课题背景海洋是一个广阔的空间，而浙江近海海域是国内较为繁忙的港口之一，船只和人员的活动密集，海域内的潮汐现象及海浪波动也极容易受到人类活动的干扰及其变化。
因此，通过对浙江近海海域的潮波进行数值模拟，有助于提高对该地区潮汐海浪的认识，加深对自然现象的理解，同时也能更好地规划海洋工程、航行等方面的活动，保障安全。
二、研究目的本研究旨在通过数值模拟技术，研究浙江近海海域的潮波现象，预测该地区潮汐海浪的变化规律，以便为实际工程和航行提供更准确的参考依据。
三、研究内容（1）研究浙江近海海域的潮波数值模拟技术及方法，了解潮汐海浪的形成机制和影响因素。
（2）对浙江近海海域的地理地貌、海底地形进行研究分析，收集相关水文气象数据。
（3）运用计算流体力学(Computational Fluid Dynamics, CFD)数值模拟软件，建立浙江近海海域的潮波数值模型，进行数值模拟和仿真。
（4）分析浙江近海海域潮波的变化规律、影响因素，并对数值模拟结果进行验证和修正。
四、预期成果（1）研究浙江近海海域的潮波现象，深入了解潮汐海浪的形成机制和变化规律。
（2）开发出浙江近海海域的潮波数值模拟系统，预测出潮汐海浪的变化情况，为实际工程和航行提供可靠的参考依据。
（3）通过模拟结果的分析，提出改善该地区潮波海浪的建议和方案，为海上交通和安全提供保障。
五、研究难点（1）模型建立：模型的建立不仅需要考虑到研究区域的地理和气象条件，还需要整合多种专业知识，对模型参数精准处理。
（2）数据采集：数据的准确性和全面性会对模型的准确性产生影响。
如何获取数据、如何筛选数据成为难点。
（3）数值模拟：通过仿真出潮波变化规律也面临诸多困难，如对水流动力学的分析、数值方法选取等。
大海域潮流场数值模拟分析
1 工程概况山心沙岛位于防城港市企沙镇东南海域，东临钦州湾，南濒北部湾。
北纬21°35′18.2″，东经108°30′51.2″，岸线长1.9km，离岸距离为340m。
西北侧陆域为山新村，是一个以海洋捕捞和海水养殖为主的渔村，距钦州港西航道4n mile、企沙镇政府1km、港口区政府25km。
工程海域地处北回归线以南低纬度地区，气候属于我国亚热带海洋性季风气候，冬季温和，夏季多雨，季风明显，受灾害性天气影响较明显。
防城港湾年平均风速为3.1m/s，月平均最大风速出现在12月份，为3.9m/s，其次是1月和2月，为3.7m/s；最小平均风速出现在8月份，为2.3m/s。
防城港的潮汐为正规全日潮，其潮位特征值如下（1985国家高程基准）：最高潮位为3.60m（2008.11.1607：09），最低潮位为-2.03m（2002.12.818：53），平均潮位为0.62m，平均大海域潮流场数值模拟分析刘华峥 杨琴琴（烟台仲伯企业管理咨询有限公司，山东 烟台 264000）摘 要：本文以防城港市山心沙岛生态岛礁建设项目为研究对象，采用数值模拟的手段，建立大海域潮流场平面二维数值模型，对项目建设前后的潮流场状况进行数值模拟。
结果表明，现状潮流流速通常＜50cm/s，工程东南侧海域流速略大；涨急时刻流速通常＜60cm/s；高潮、低潮时刻工程潮流流速均＜10cm/s；工程建成后，西南和北东侧涨急时刻流速主要呈现变小的趋势，最大降幅约17.5cm/s，南、东南和西北侧涨急时刻流速呈现增大的趋势，流速最大降幅约为18.5cm/s；涨急时刻流速变化＞5.0cm/s的区域主要集中在工程区周边0.8km内。
关键词：大海域；潮流场；数值模拟；水动力；防城港中图分类号：U 652 文献标志码：A并同时在智能终端界面显示故障具体信息，从而帮助技术人员及时解决故障问题。
本文设计的煤炭洗选系统工艺流程如图3所示。
渤、黄、东海潮汐、潮流的数值模拟与研究
渤、黄、东海潮汐、潮流的数值模拟与研究朱学明;鲍献文;宋德海;乔璐璐;黄必桂;石新刚【期刊名称】《海洋与湖沼》【年(卷),期】2012(043)006【摘要】基于FVCOM海洋数值模式，采用高分辨率的三角形网格，对渤、黄、东海的潮汐、潮流进行数值模拟，并通过比较120个沿岸验潮站和14个潮流观测站的实测与模拟结果进行模型验证，两者符合较好。
根据模拟结果，给出了四个主要分潮的潮汐同潮图和5m层潮流最大流速及最大潮流同潮时分布。
渤、黄、东海共有5个半日分潮和3个全日分潮的独立旋转潮波系统，且都呈逆时针方向旋转；半日潮流和全日潮流各有12个圆流点：在冲绳岛和奄美岛两侧的4个半日潮流圆流点分别呈对称分布，其中有3个为本文首次给出；在日本九州岛西侧还新给出2个全日潮流圆流点。
有关它们的存在性需要实测资料的进一步检验。
%A successful simulation for the tides and tidal currents has been done based on the Finite Volume Coastal Ocean Model (FVCOM) with unstructured triangle grids in high resolutions in Bohai Sea, Yellow Sea, and East China Sea. The numerical model agreed with the observation well by comparing the observed data with the simulation results from 120 gauge stations and 14 tidal currents stations. The distributions of the co-tidal charts, maximum velocity and co-range charts when maximum velocity occurs of the tidal currents in 5 meter layer are shown in this paper according to the model results, respectively. Model results indicate that five semidiurnal and three diurnal constituents independent rotary tidal systems appear in Bohai Sea,Yellow Sea, and East China Sea, and they are all anti-clockwise; twelve tidal cur- rent-amphidromic points for semidiurnal and diurnal constituents can be found, respectively; the four semidiurnal tidal current-amphidromic points located at the two sides of the Okinawa Island and Amami Island are symmetry, three of them are reported at first time; two diurnal tidal current-amphidromic points located in the west of Kyushu Island of Japan are reported at first time, too. But the existence of them should be validated by observation.【总页数】11页(P1103-1113)【作者】朱学明;鲍献文;宋德海;乔璐璐;黄必桂;石新刚【作者单位】中国海洋大学海洋环境学院,青岛266100 国家海洋环境预报中心国家海洋局海洋灾害预报技术研究重点实验室,北京100081;中国海洋大学海洋环境学院,青岛266100;中国海洋大学海洋环境学院,青岛266100;中国海洋大学海洋地球科学学院,青岛266100;中海油研究总院,北京100027;中海油能源发展股份有限公司北京分公司,北京100027【正文语种】中文【中图分类】P731【相关文献】1.渤黄东海潮流长期物质输运研究 [J], 乔方利;赵伟;袁业立2.渤、黄、东海二维潮汐模式底摩擦系数的反演研究 [J], 张继才;吕咸青3.海平面上升导致渤、黄、东海潮波变化的数值研究Ⅱ——海平面上升后渤、黄、东海潮波的数值模拟 [J], 于宜法;刘兰;郭明克4.海平面上升导致渤、黄、东海潮波变化的数值研究Ⅰ——现有的渤、黄、东海潮波的数值模拟 [J], 于宜法;郭明克;刘兰5.渤,黄,东海潮汐潮流的数值模拟 [J], 赵保仁; 方国洪因版权原因，仅展示原文概要，查看原文内容请购买。
黄河口海域潮汐、潮流、余流、切变锋数值模拟
黄河口海域潮汐、潮流、余流、切变锋数值模拟
答：黄河口海域潮汐、潮流、余流、切变锋数值模拟采用有限体积三维近岸海洋模型,建立无结构三角形网格,对当今地形下的黄河口海域进行高分辨率数值模拟,在成功模拟渤海潮汐潮流的情况下,重点研究了黄河口海域。
当前黄河口附近海域潮汐为不规则半日潮,潮流为往复流,方向近似平行于岸界,潮致欧拉余流在岬角两侧存在成对的涡旋,涡旋的方向为南顺北逆,黄河径流对此涡旋有加强的作用。
由于地形等复杂因素的影响,河口海域附近在涨落潮转换过程中存在内涨外落型和内落外涨型切变锋,其首先出现在浅水区域,然后向深水区域传播,1~2h后消失,它的产生是由于近岸海域潮汐相位领先于外海海域潮汐相位。
温州鳌江近海建设工程环境影响潮汐潮流数值模拟
温州鳌江近海建设工程环境影响潮汐潮流数值模拟吕和娜;汪一航;张钊;单慧洁;岳云飞【摘要】This paper have col ected the measured tides and certain tidal current data in different stages of many projects during past three to five years near the Aojiang River. The harmonic method is used to analyze tide and tidal current data observed at five stations in the sea adjacent to Aojiang River. The results show that the tide is mainly regular and semidiurnal in the sea near Aojiang of Wenzhou. The tidal amplitudes ofM2 constituent are between 170 cm-193 cm and the lags are between 260°-280°. According to the comparison of analytical results of harmonic constants of these stations in 2007, 2010 and 2011, the maximum change of tidal amplitudes and phase-lag range for the main semidiurnal tides (M2, S2, N2), the diurnal tide (K1, O1) and the shal ow water tide (M4, MS4, M6) are 1.8 cm-4.4 cm and 3°-7°, respectively. After analyzing the tide and tidal current characteristics of Aojiang River, this paper uses an unstructured grid and Finite-Volume Coastal Ocean Model (FVCOM) to test the results of simulation. The simulated results agree well with the measured data. The new shoreline and depth which are produced by the construction projects closed in important major years, and the tide and tidal current field for the new shoreline and depth are constructed, which describe the superimposed influences of construction engineering in Aojiang estuary.%收集了近年来鳌江口附近海域多个工程不同阶段5个潮位站的3～5年潮位实测数据和部分海流实测资料,通过对鳌江口附近海域的不同年份的水位资料进行潮汐调和常数分析,鳌江近海海域主要为半日潮区,其中M2分潮的振幅在170 cm ～193 cm；迟角在260°～280°之间,这些站的2007年、2010年、2011年调和常数分析结果相比,主要的半日分潮M2、S2、N2，全日分潮K1、O1及浅水分潮M4、MS4、及M6等分潮振幅、迟角的最大变化分别在1.8cm ～4.4 cm和3°～7°之间。
东中国海潮流三维数值模拟与分析
东中国海潮流三维数值模拟与分析作者：庄桦来源：《中国水运》2014年第01期摘要：本文建立了东中国海潮波三维数值模型，根据模拟结果得到该海域表层、中层、底层的半日分潮流M2和全日分潮流K1分布场以及潮流类型分布图，并对潮流的垂直结构特征进行了分析。
关键词：东中国海潮流类型垂直变化潮波包括潮汐和潮流，指海水在天体（主要是月球和太阳）引潮力作用下所产生的周期性运动，习惯上把海面垂直方向涨落称为潮汐，而水平方向的海水流动称为潮流。
潮流是是海水的基本运动之一，也是海洋动力学研究的重要组成部分。
研究和探讨潮流的基本特征、分布规律有着极其重要的意义。
模型建立模式的计算区域为117°30′（E ～ 131°E，24°30′（N ～ 41°N，包括渤海、黄海、东海、东海陆架坡、琉球以东西北太平洋和部分日本海，开边界选取如图1所示。
研究区域的地形资料，包括来自NOAA的全球地形资料，长江口、江苏北部沿海、渤海的海图资料。
开边界处采用松本国家天文台的海洋潮汐预报模式NAOTIDE提供的水位数据，水位、流速验证采用2008年潮汐表的潮汐和潮流数据。
采用ECOMSED模式，网格水平分辨率为1′×1′，垂向在？滓坐标下分三层。
主要研究天文潮及其传播变形的模拟，因此采用正压模型，温度和盐度为常数，亦不考虑风场和海气交换等因素，故表面热通量、盐通量、风应力分别赋为零值。
经过反复试算调试，90个站点的水位和18个站点的流速、流向验证，所得计算结果较为满意。
潮流垂直结构特征1、分潮流椭圆的垂直变化这里以M2、K1分潮流为例，研究分潮流最大流速及方向、椭圆率随水深的垂直变化，图2—图7为表层、中层、底层M2、K1分潮流椭圆图。
从图2可以看出M2分潮流强潮流区在朝鲜西海岸，最大潮流速超过1.0m/s。
渤海海峡东侧、莱州湾沿岸、东海东南部的比较弱。
M2分潮流的主轴方向与潮波传播的方向一致，潮流主轴在东海为东南向西北，进入黄海后慢慢转为南北方向。
东中国海潮汐能分布数值模拟研究
东中国海潮汐能分布数值模拟研究顾振华;顾晨;姚鹏程【期刊名称】《中国农村水利水电》【年(卷),期】2022()11【摘要】潮汐能是一种可再生、绿色无污染且蕴藏量巨大的可再生新能源,潮汐能资源分布特征的研究对于我国潮汐能的开发利用具有十分重要的现实意义。
当前潮汐能资源的调查较多为对离散点资源的分析统计,缺少对潮汐能资源进行连续全海域的计算研究,为了了解全海域资源的分布,本研究建立了东中国海潮波运动数值模型,数模结果与各沿岸测站以及前人基于实测数据给出的同潮图均吻合较好。
在此基础上结合数值模拟结果和近年发布的国家标准评估了东中国海潮汐能资源技术可开发等级、潮汐能理论平均功率和年蕴藏量等的空间分布规律。
结果表明受月球轨道变化的影响,各验潮站年平均潮差变化周期(约19年)和相位均相等而变幅则各不相同;东中国海大部分海域平均潮差较小,为潮汐能资源贫乏区,但朝鲜半岛西岸、杭州湾和台湾海峡等小部分海域平均潮差达4.5 m以上,为潮汐能资源丰富区或较丰富区,其潮汐能理论平均功率和年蕴藏量可达2×10^(3)kW/km^(2)和2×10^(8)kWh(/a·km^(2))以上;浙江福建和苏北沿海较大部分海域平均潮差达3.0~4.0 m左右,主要为潮汐能资源可利用区,潮汐能理论平均功率和年蕴藏量在1.3×10^(3)kW/km^(2)和1.3×10^(8)kWh(/a·km^(2))左右;与平均潮差相比潮汐能年蕴藏量和理论平均功率在不同海域间空间分布更不均匀。
【总页数】7页(P165-171)【作者】顾振华;顾晨;姚鹏程【作者单位】上海勘测设计研究院有限公司【正文语种】中文【中图分类】TV213.1【相关文献】1.东中国海拉格朗日环流数值模拟Ⅱ.环流模拟2.2001年春季东中国海三维环流的数值模拟3.东中国海潮流结构数值模拟与分析4.东中国海潮能通量与耗散的数值模拟分析因版权原因，仅展示原文概要，查看原文内容请购买。
浙江近海潮汐潮流的数值模拟
浙江近海潮汐潮流的数值模拟陈倩;黄大吉;章本照【期刊名称】《海洋学报（中文版）》【年(卷),期】2003(025)005【摘要】用三维陆架海模式(HAMSOM)对浙江近海的潮汐、潮流进行了数值模拟,并采用网格嵌套和动边界技术对原模式作了改进,以提高计算的精度,改进后的模式在浙江近海的应用中被证明是成功的.沿岸50个潮位站计算与实测值的比较表明,加入动边界以后的小区域细网格计算较之粗网格以及未加动边界以前精度普遍提高,比较的均方差结果为:M2分潮振幅差4.6 cm,相角差7.14°;S2分潮振幅差5.0 cm,相角差5.4°;K1分潮振幅差2.25 cm,相角差5.76°;O1分湖振幅差1.56 cm,相角差5.5°,可见计算与实测符合良好.另外,选取了105个实测潮流点,比较了表层M2和K1分潮流调和常数分量UCosξ,Usinξ,Vcosη,Vsinη的实测值与计算值的偏差,结果表明计算与实测的符合程度较好.在此基础上,给出了各主要分潮的潮位同潮图、潮流同潮图、潮汐性质、潮流性质、潮流椭圆和潮流的运动形式等,发现4个主要分潮M2,S2,K1,O1在本区内均未出现无潮点;M2分潮流在29°18＇N,122°46＇E 处有一个圆流点.此外还得到了一些有意义的结论,都与实测情况符合良好,从而对整个浙江沿海区域的潮汐潮流特性有了一个全面认识.【总页数】12页(P9-20)【作者】陈倩;黄大吉;章本照【作者单位】浙江大学,力学系,浙江,杭州,310027;国家海洋局,第二海洋研究所,浙江,杭州,310012;国家海洋局海洋动力过程与卫星海洋学重点实验室,浙江,杭州,310012;国家海洋局,第二海洋研究所,浙江,杭州,310012;浙江大学,力学系,浙江,杭州,310027【正文语种】中文【中图分类】P731.23【相关文献】1.温州鳌江近海建设工程环境影响潮汐潮流数值模拟 [J], 吕和娜;汪一航;张钊;单慧洁;岳云飞2.温州鳌江近海建设工程环境影响潮汐潮流数值模拟 [J], 单慧洁;张钊;汪一航;吕和娜;申伟杰3.泰国湾及邻近海域潮汐潮流的数值模拟 [J], 吴頔;方国洪;崔欣梅;滕飞4.基于FVCOM的獐子岛附近海域三维潮汐潮流数值模拟 [J], 齐继峰;曹圣山;郭可采;杨德周;徐振华;尹宝树5.基于 FV CO M的温州近海潮汐潮流数值模拟 [J], 岳云飞;王永刚;何善方;汪一航;魏泽勋因版权原因，仅展示原文概要，查看原文内容请购买。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
中国近海潮波运动数值模拟
本文基于球面坐标系下的二维垂线平均潮波运动方程建立中国近海潮波数学模型，模型区域包括渤海、黄海、东海、南海、泰国湾和环台湾岛海域，网格尺寸2′×2′，网格数1175×955。
在考虑引潮力情况下，计算模拟了中国近海的复合潮波运动；并对分布于各个海域的281个潮位站的4个主要分潮
（M<sub>2</sub>、S<sub>2</sub>、K<sub>1</sub>、O<sub>1</sub>）潮位调和常数以及13个海洋预报站的潮流资料进行了验证，验证结果基本合理。
针对计算结果绘制了主要分潮包括浅水分潮（以M<sub>4</sub>为例）的潮汐同潮图和潮流同潮图，对中国近海潮汐和潮流分布即潮波运动进行了分析，并和前人的结果进行比较，结果基本吻合。
整个东中国海的潮波主要是太平洋潮波经台湾和九州之间的水道传入的协振波，南海的潮波主要是太平洋潮波经吕宋海峡传入的协振波。
东海和南海主要通过台湾海峡进行水量和潮能交换。
在东中国海基本以半日潮为主，尤其是
M<sub>2</sub>占优，而在南海基本以全日潮为主。
由于受到地形影响、边界的反射、地转偏向力和陆架浅海的摩阻作用，潮波在各海区或以前进波或以无潮点和圆流点为主要特征的旋转潮波系统组成了复杂的潮波系统。
在模型计算的基础上对台湾海峡的M<sub>2</sub>分潮的潮汐分布特征和传播机制进行了探讨。
认为台湾海峡的M<sub>2</sub>分潮主要是有北部的前进波和南部的前进驻波系统组成，由吕宋海峡进入的太平洋潮波和广东、福建沿海岸线的相互作用形成了南端前进驻波现象和北部地形边界的放大效应产生的潮能幅聚现象是台湾海峡M<sub>2</sub>分潮分布特征的主要原因。
海峡内的潮能消耗以南下的前进潮波为主，其传播到台湾浅滩位置。
另
外，本文在潮波数值模拟技术方面作了以下探讨： 1．利用相关分析讨论周期性函数的一阶差分离散误差和相对步长的关系，得到其中的拟合关系式；2．基于相关分析提出利用计算结果和实测资料的相关分析结果进行模型验证，和其它模型验证方法比较该验证方法具有简洁且可同时考虑振幅和相位误差的优点； 3．通过数学模型和原型的相关理论推导了二维潮波运动方程中曼宁系数的预估校正法。
数值试验和中国近海潮波数学模型的实际应用表明该方法科学有效，既可以节省人工调试的工作量，又能提高计算结果的精度。