船用潮汐、潮流图表编制方法

船用潮汐、潮流图表编制方法

国家标准 GB/T13474一92 船用潮汐、潮流图表编制方法 Themethodofeditingtidaltablesand tidalstreamchartsforship 1992-06-04发布 1993-03-01实施 发布 国家技术监督局国家标准
国家标准 船用潮汐、潮流图表编制方法 GB/T13474-92 Themethodofeditingtidaltablesand tidalstreammchartsforship 主题内容与适用范围 本标准规定了主要港口和水域的船用潮汐、潮流图表的计算及编制方去〔所用符号见附录E(参考 件)) 本标准适用于我国沿海及感潮江河各主要港口和水域 术语 2.1潮汐Tides 由于月球和太阳等天体引潮力的作用产生的海面周期性升降运动 2.2潮流Tidalstreamandtidalcurrent 伴随着潮汐运动,海水在水平方向的周期性流动 2.3潮高基准面Tidaldatum 潮位高度的起算面 2.4潮高(潮位Theheightofthetide 任意时刻水面到潮高基准面的垂直距离 2.5潮差Tdaltange 两相邻高潮和低潮的(潮高)差 2.6潮振幅(半潮差》Tdalamplitude 相邻高潮和低潮间潮高差的二分之- 27海图深度基准面Chart.datum 海图中水深的起算面 2.8水尺零点(验潮仪零点Zeroofgaugeorzeroofpole 进行水位观测的水位高度起算点 2.9平均海面Meansealevel 某一时期(如一日、一月、一年和数年)海面的平均高度 2.10高低潮时差Timediffterences 主港与副港出现高(低)潮的时间差 2.11潮差比 Ratioofranges 半日潮港的副港与主港的平均潮差之比,日潮港的副港与主港的回归潮最大潮差之比 2.12改正值Correction 由于主,副港的潮高基准面不一致,所引入的修正量 2.13潮汐间除Tidalinterval 国家技术监督局1992-06-04批准 1993-03-01实施
GB/T13474-92 月球经过某地子午线至该地发生第一个高(低)潮的时间间隔 2.14平均大潮升Meanspringrise 大潮期间的平均高潮潮高 2.15平均小潮升Meanneaprise 小潮期间的平均高潮潮高 2.16潮信Tidalinformation 反映港口潮汐特性的有关信息 2.17差比数Ratio 主、副港间潮时差、潮差比和改正值的总称 2.18正规半日潮Regularsemi-diurnaltide 在一个太阴日里(约24h50min)海水出现两次涨落,且相邻高低)潮的时间间隔约为12h25min, 两相邻的高低)潮潮高接近,涨落潮时间约相等 2.19正规日潮Regulardiurnaltide 在半个太阴月中,至少有一半的天数每天出现一个高渐和低潮 2.20不正规半日潮混合潮 guarmixed Irre semidiurnaltide 在一个太阴日中发生两次高潮和两次低潮,但相邻的高低)潮的潮高差显著不等 2.21不正规日潮混合潮 Irregularmixeddiurnaltide 在半个太阴月中每天发生一次高潮和一次低潮的天数少于七天 222回归潮Tropictide 当月球赤纬最大时,在混合潮区两相邻高<低)潮的潮高不等达到最大的现象 223分点潮Equinoetialtide 当月球位于赤道附近时,在日潮或混合潮区,日潮变得很小,潮汐呈半日潮特征 2.24大潮和小潮Springandneap 半日潮港在朔或望及其后的-至三日,由于月球和太阳引力所引起的潮位相合,使潮差很大,称为 大潮,在上弦或下弦及其后一至三日两者位相相反,潮差很小,称为小潮 2.25潮汐表Tidetables 刊载有若干潮汐港口和站点的潮时和潮高预报值的表册 2.26潮流表Tidalstreamtables 刊载有若干潮流站点的流况预报表册 2.27潮流图Tidalstreamcharts 以矢量围的形式提供某海区的潮流预报 编制潮汐表时采用验潮资料的技术要束 潮位观测资料的质量直接影响到调和分析的结果,经整编的验潮资料必须包括下述内容(按月报表 整编资料). 时间(年、月、日、时); 台站名,经、纬度C以(),'),(")表示; 基准面位置<以m为单位); 使用时间标准[我国沿海台站采用北京标准时(东8时); 每天正点潮高,高低潮潮时,潮高(潮高精确到cm,潮时精确到min) 以月为单位提供下列统计特征值: 最高潮高及对应的日、时、分， 最低潮高及对应的日,时、分;
GB/T13474一92 平均高潮潮高， 平均低潮潮高， 月平均海平面 备注;用以记录潮位观测过程中异常情况及技术处理过醒,对该潮位观测时段进行质量品评 并记录潮高基面的变动情况 每一份潮汐观测月报表必须标明台站负责人、抄录人、校对者、复核人的姓名 潮汐调和分析 根据所采用潮位观测资料的时间系列不同,可分为; 短期资料;有一天、二天、三天、七天,十五天、三十天等 b 长期资料;有三百五十五天、三百六十九天或多年资料 短期潮位观测资料的调和分析 对于一次,二次,三次或三次以上周日观测的潮位观测资料的调和分析方法与短期观测资料的潮流 调和分析方法类同(参照7.1条). 4.2长期潮位观测资料的调和分析 长期潮位观测资料按最小二乘法进行调和分析(见4.3条) 4.3溯汐最小二乘法调和分析 4.3.1潮高方程的建立 某时间计算潮高可用式(1)表达 十 :Mco+w. H,()=a,- u)j一g (1 式中;H,(e) -时间的计算潮高,cm; 水尺零点起算的平均海面,cmr ao” 分潮序号; 分潮总数; j分溯的交点因数， f, 分潮的平均振幅,cm; H j分潮的角速度,rad/s; o 平太阳时,h; w,十n)j分潮的天文初位相角,), -j分潮的专用迟角,". 8- 43.2潮汐最小二乘法分析主要步骤 潮汐现象是一随时间变化的连续过程,故又可用式(2)表达: H,C(m十n)]=a,十 晶叫(-曾+盏叫(-凯叫) 2 式中 “一HcaG,-w十)-(n+尝p 3
GB/T13474-92 ("+“们 占=手,H,sin[s,-(V,十n),- 实测资料第一个观测值距该年元旦零时的小时数,h; n -1h; m-自然序数,m=0,1,2,,M M实测潮位资料的总个数 对参与分析的M+1个实测资料作平方通近,使实测潮高与计算潮高间的偏差平方和最小,即 (5 H,(m十n))'= min 0a(a十)n- 式中;H[(m十n)某时间的实测潮高 根据最小二乘法原理,取 [2D ao 6 p 0 [品 4.3.3关于a,a,b的正规方程 ,N H.(m十n)mce("-小) i=0,1,2, (7 M ,N 一 i=1,2,* C-习t(m十n)力ai(( 2 式中;>" 一表示求和时首尾两项的权系数为 ,N -====- F i夭i,j=1,2， * ,N F，= [u+性雨 i=l,2，" E =M 9十色 G= i夭j;i,j=1,2,,N lH半一半闪 9 ,N G M一性 i=l,2,"
GB/T13474-92 sin c0s (10 U o g" 对正规方程的右端项,拟采用watt迭代法,以节省运算时间,按常规算法求解方程 4.3.4各分潮调和常数计算公式 A,=a,十AA, H=(;十时)tJ (11 -es[(/e)+?,十幻+(n+啊] lg= 式中;A,由潮高基准面起算的平均海面计算值,cm A水尺零点与潮高基准面之差,em 4.3.5对调和常数检验的要求 基本符合测站主要分潮调和常数变化的历史趋势; b. 基本符合所在海区的潮汐变化规律， 分析所得调和常数用于实际预报时,在正常情况下,预报误差潮时在20一30min以内,潮高在 2030cm以内 4.4潮汐计算 潮汐计算公式 4.4.1 Mmco H,()=A 十+ 十(V 十u);一gs (12 4.4.2正点潮高的计算 结合具体计算站点,选定分潮数N及其各分潮的调和常数,确定计算的初始时间(e,)及终止时间 e),增量为1h,应用式(12)计算逐时正点潮高 4.4.3高低潮潮时、潮高的计算 在非单调变化的三个正点潮高所夹的两个连续区间中,采用压缩区间二分法逐步逼近求出极值潮 时;或采用拉格朗日插值公式拟合求解极值潮时,再将极值潮时代入式(12)求得相应的极值潮高 5 差比数和潮信 5.1潮汐分类(潮汐性质 采用主要全日分潮(K,O.)与太阴主要半日分潮(M.)的平均振幅之比值进行分类 Hk，+H Hr (13 HM 式中;Hm -潮汐类型特征值 s -K,分潮平均振幅;
GB/T13474-92 -O分潮平均振幅; Ho HM -M,分潮平均振幅 H<0.5为半日潮， Q.54.0为全日潮 5.2 浅水影响指标 以太阴浅水四分之一日分潮(M)和太阴主要半日分潮(M.)平均振幅之比(H/H)值确定浅水 影响的程度 5.3高低)潮间隙 高低)潮间隙按式(14)计算 半日潮、分点潮平均高潮间隙w值计算公式; (14 wI=(s的,一v)/28.98 -M,分潮的专用迟角,(); 式中，发叫 v半日潮平均高潮间除浅海改正辅助角,(),CG计算公式见附录A(补充件)). 半日潮,分点潮平均低潮间隙wL值计算公式 15 WI=gM，+180一w/28.98 式中:W " -半日潮,分点潮平均低潮间隙浅海改正辅助角,('),(计算公式见附录A>. 回归潮平均高高潮间隙wImm计算公式; c. (16 wlm=cGs十)/2一w/14.6 式中h,一K分潮专用迟角,o, -O分潮专用迟角,() go -回归潮平均高潮间隙浅海改正辅助角,(),(G计算公式见附录A). mmin 迷，当(s,十e)2与瓜,或狐之差超过9时,应减去180当wlmw大于24hs0mn时应去2h50 d. 回归潮平均低低潮间隙wIwmL计算公式 17 wIT (G,十,)2十w-w>h4a 式中;w回归潮平均低潮间隙浅海改正辅助角,("),(计算公式见附录A. 注，当求出值大于24h50min时应减去24h50min. 5.4平均高(低)潮高(在潮高基准面上 541回归潮平均高高精高Hmm值计算公式
GB/T13474-92 -cea必十&e试,+战 十H.cosg HTPH w.oasH v,sin s,十&)一 G-(,十鸟,)十4 5.4.2回归潮平均低低潮高Hm值计算公式: Hmu (ow十0.021)(H十H,)+H,oeCs叫 s,十,]一Q.035H,w,sin 19 s十s)十A sx 5.43分点潮平均潮差4H,值计算公式: n-u,(+H'(最'+'(能" H没'朵-皮震7 Q.17(H十/H4十2H, 20 5.44分点潮平均高潮潮高Hm值计算公式 Hn=A 十AH/2 (21 5.45分点潮平均低潮高H值计算公式: H礼=A,一AH,/2 (22 5.5平均大(小潮升 5.5.1平均潮差AH值计算公式 [是 H我=2Ha+习 -N,A,Aa月a, + H(oov,十cew,一)+0.05H(?,-w) aincs的 23 g十2HM,cos(3g 8M. 5.5.2平均大(小)潮差AH,(AH,)值计算公式 十H AH, Ho -0.oe )战 最-儿"- -AH一 24 H (AH， 5.5.3平均大潮升(在潮高基准面上)H,值计算公式
13474-92 GB/T 战-山十aH" 25 5.5.4平均小潮升(在潮高基准面上)H,值计算公式 26 H =A,AH,/2 5.5.5回归潮大潮差AH,值计算公式 27 AH,一Hrwmw一星mt 5.6差比数 5.6.1主、副港潮时差值按下列公式计算 半日潮,不正规半日潮混合潮 T =WIp一 -wie- (S，一S2十n，12h50nin 28 (S,一S2十n (12h50min [4T一ww-wua- 式中;Ar -主、副港高潮时差,h; ATa -主、副港低潮时差,h; Saxu) 主(副)港标准子午线,以h表示， -日潮不等指标(取值法见附录A); n0 wIn -主港平均高潮间隙， -副港平均高潮间隙; WIe -主港平均低潮间隙， WI wu -副港平均低潮间隙 注，w彻、,wm,wa和w心都需完成浅水订正 不正规日潮混合潮、日潮 b -WwI 4T ，=wI. (S,一S) 上TH 1TP .29 lT=wImu,一wIm GB/厘13474-92 Pa (AH);/(AHx" 31 -主(副)港回归潮大潮差 式中(AHm) 6 编制潮流图表时对潮流观测资料的技术要求 潮流观测资料的排序按下列要求; 按年、月、日先后排序 6.2以i表示潮流周日观测的序数(i=1,2,,S,. S,表示潮流周日观测的次数 以j表示每次周日观测中每组记录的编号C-1.,2, 6.3 ,n(i) n()表示第i次潮流周日观测资料的组数 潮流调和分析的方法 7.1短期潮流观测资料的调和分析 当潮流周日观测的天数少于30d时,按潮流准调和分析法或潮流分析的调和算法求得调和常数 7.2长期潮流观测资料的调和分析 当潮流周日观测的天数大于30d,可将观测资料按潮汐最小二乘法求得各分潮的调和常数(参照 4.3条) 潮流准调和分析的方法 本方法分两种,即引入差比数法(简称“方法I")和不引人差比数法(简称“方法【”). 注，1)以一个分潮为主,对其余分潮的调和常数给出它们与这个分潮调和常数的一定关系,并将它们合并至这主要 分潮为一项,称这一项为准调和分潮 8.1良好天文日期的条件 8.1.1“方法1”良好天文日期的条件,对于能计算出观测海区附近四个主要分潮(Ms,K,O,)潮汐 调和常数之间的差比数,且周日观测天数小于7d的测站采用此法选择良好天文条件对应的测流日期 12° 32 AC，=do一dk Ag - 12 e== (33 AC，=dM,一d，一Ag 式中dM,d,dsdo -准调和分潮的迟角订正,('); Ag-K和O分潮的迟角差,(°); -M和S分潮的迟角差,() Ag， 我国沿海各海区,和4g,的取值见附录B(补充件). 270"GB/T13474-92 (D,入,(D,,cosC(d入,-(d)) 习(Do,(D,sin(d),一(d),))}'=0 34 半日分潮最优天文条件的指标量o计算公式 b. (D,(D,),cos[(d,-(d))) D,(D,),sinc(d,)-(d) =0 (35 式中;Du,D，、Dk、Do1 准调和分潮的振幅系数; 表示第i次周日观测 8.2潮流观测资料的矢量分解 将实测流向、流速资料按直角坐标系分解为北、东分量 yU(i,)=wi,)eos((,7 36 v(,)-w(,)sia(,) 第i次周日观测第j个记录的流速,m/s, 式中，w(i,) 第次周日观测第j个记录的流向,()1 (i,) 第次周日观测第j个记录的北分量 U(i,) 第次周日观测第j个记录的东分量 v(i,) 8.3天文变量的计算 天文变量计算见附录A. 8.4计算潮流调和常数方法I 8.41计算观测日零时的天文海区参数按下列公式 全日分潮 瓦=可,十a,4Hcs山s,一a,AHsimAs 互一可+可,AHsinAx+a,AH,coaAe (38 半日分潮 o 万=西,+a,AHcosAg，一a,AHsins" 39 瓦，=a+a,AH,sing,+a,AH,eosg 1/4日分潮 c. 10
GB/T13474一92 万=,十a,4Hea么g,-a,aHsinAk" 4 瓦=a十a,AHsinAg，十nAHlcosAs. 42 式中,友天文海区参数(K-1,2,,>. 8.4.2计算观测时刻对应的天文海区参数按下列公式 醒24-硬o7 秘》=秘(o)+ td,)" ? e48)+eo)一2P24 .sncoa 24们 43 1152 一当D=0时第i次观测对应的天文海区参数 式中，硬o一 (24) -当D=1时第i次观测对应的天文海区参数 硬(48) -当D=2时第i次观测对应的天文海区参数 8.4.3建立矛盾方程组 r(i)N'X+r(i)Xt=r(iu(i,j 44 r(6)Y十r()Yr=r()vG,) 45 r()的取值见附录D(参考件) 8.4.4求解正规方程 Ck.X，=UK=1,2,,6十S 46 s Cx.Y=Vx(K=1,2,,6+S 8.4.5观测误差的方差及系数按下列公式计算 用正规方程的解计算各观测时刻的潮流分敢 人 [U(i,)=Xsti 硬少Xx.(i=1,2,,S 47 no [V(i,j)=Y十习:Yx(j=1,2,,n . 计算方差 11
GB/T13474-92 d (o6)-G)+vG)-'(G,)m) 48 -l-1 n( 计算误差系数 f=2.56十1.738la[习n(i)]十0.0096ln n() 49 8.4.6当某组观测数据满足下列不等式时， U(i,)一乙.(i,)]'十v(i,)-(i,))>eo叫 50 则应将对应这组的观测数据舍去,并以,(G,),Pd,)代替它们 &47调和常数和余流的计算 将经过处理后的数据重复8.4.38.4.4的计算过程 北分量调和常数U.值计算公式 将求解正规方程的结果进行代码转换[见附录c(参考件刀,而后求解北分量调和常数; U.=Cx(x",)力(c=O,,K,M,,s,,M.,Ms 51 [x/2 若X，=0,X' >0 3x/2 若X =0,X'.<0 0 若X.>0,X'.> g (52 tg 若X <0 若x.>0,x',GB/T13474-92 灭/2 0,Y4>0 若X+ 3x/2 若Xst，=0,Y1 tg 若X， 0,Yt> 6+ a(i 54 若X tg 十T Y 十2x 若X>0,Y,4<0 tg-' 着4i 8.5计算潮流调和常数-方法I 8.5.1计算观测时刻对应的天文参数 创》=e0)+壁24-eon 24 a(48)十a(0)一2a(24 2G.) 1152 Ctd,) 24] (55 式中;aP0)当D=0时第i次观测对应的天文参数; a(24)--当D=1时第i次观测对应的天文参数， aP(48) -当D=2时第i次观测对应的天文参数 8.5.2建立矛盾方程组 党 r(G)ax".十r()X，=r(i)U(i, 56 " Br(6ae"y".十(oY. (57 'ti=r(i)v(i, 8.5.3求解正规方程 C'xX",=U'x(K=1,2,,12+S, 12+s 58 C"xaY"=v'x(K=1,2,,12十S, 8.5.4观测误差的方差按下列公式计算 用正规方程的解计算各观测时刻的潮流分量 ,G,) =X"这ti+ afX". 59 [,G,) =Y" a爱y" "+ 13
GB/T13474一92 计算方差 b 计算公式与&.4.b项的式(8)类同,将式中变量.(G)和vG,)相应置换成o,Gi,)和 .(i,). a.5.5当满是下列不等式时,则将对应这组的观测数据含去,在以后的计算中不再使用、 (60 wG,)-o,G.)+vG.)-y,G,)y>小ei 含去不合理的观测数据,并以U,(G,),,G,)代替它们,然后再按8.5.2求解矛盾方程组 8.5.6调和常数和余流的计算参照&.4.7的方法,公式中变量作相应的置换 g 潮流分析的调和算法 采用本算法可直接推算出七个分潮的潮流调和常数,其中包括两个半日分潮,两个全日分潮和三个 浅水分潮 在投入计算前要求将s,个周日非整点观测资料按8.2提供的方程进行矢量分解,然后按三点拉格 朗日插值公式;求北分量相东分盖的整点值 9.1北分量调和常数的计算 某时间湘流北分量的计算值可用下式表达 U,(o)一对+习x'es(G-12十r))a)+习?"sin(G一12+Ar)a. n=1,2,,7 (61 式中 62 X”=R;"os{喊一Cn,(1)+12S十d.Ja, 63 Y:”一e'sin{喊一Cn,(Q)十125s,十d.Ja. R"=.U 64 65 究 =e 一(V，十u). Ar()=An,()一12(S,一1)一d 66 (67 d. 68 d=0 69 d，=n,(Gi)一n,(i一1)一24， 70 An,(i)=n,()一n,(1),i=1,2,,S n()为第i次周日观测,第一个整点观测值由元且零时起算的小时数 将测站邻近海域潮港对应六个分潮(M,,s.,K,,o,,M,Ms,)的调和常数,组成三个分潮对(M与 s为一分潮对,K,与O为一分潮对,M与Ms,为一分潮对),速立大个但等式作为约束条件 14
GB/T13474-92 X一C.(Xcosa+YH'sina, r-Ca(xyosea Ysina X”一C(XYcosas+Y?'sin ina4 71 =Yy一C(YYcosds一X"sina A=x?一Ca(xfwea十Ysine, B=Y一C,(Ycosdss X"'sinan 式中;aa,a,as -对应分潮的位相差; 对应分潮的振幅比. C,C,C一 引人Lagrang乘子、,K取12,34,56(下标12表示M与s分潮对,下标34表示K,与o分 潮对,下标56表示M与Ms,分潮对). 按最小二乘法原理 w.G.D )一x:-习x'cosG一12十Ar(G))a =min(72 习xamG一2十aro叫'+习aA十所A 式中，U,(,)实测潮流北分量的整点值 对x“和Y?求导,通过变换,使原先的21阶线性方程组简化为15阶线性方程组,按常规算法即可 求得北分量的调和常数 U =xy U,=c(Xx;)'十(Y)幻/. 73 'cG:)/(x:)十(Cn(1)+13s,十d4+?十n). =tg g.2东分量调和常数的计算 采用计算北分量调和常数同样的工作流程,即可求得东分量的调和常数 v=x; v.=[(x:+(Y5)l/f. 74 .=g-'cY)/x:)十(Cn(a)十12s,十d)e.+(w,十). [刀. 9.3余流w,,a的计算公式与(53)式相同 10 潮流裤圆要素的计算 0.1c分潮流最大流速wL()值(潮流椭圆长半轴)的计算 w.(c)=(M十N. (75 15
GB/T13474一92 76 M.=U+V十2U.Vsin(G 一%刀t N =U.+v.一2U.v.sin( 一死)) 10.2c分潮流最小流速ws,(c)值(潮流椭圆短半轴)的计算 w,(e)= IM.一N .78 一= 当 0.3c分潮流椭率K.值计算 79 K =(M 一N/(MN. 当,K.>0,表示c分潮流按逆时针方向旋转 K.<0,表示c分潮流按顺时针方向旋转 10.4c分潮流藏大常夜出现时间值计算 80 A =U cos十ecos*仇 81 战=Ugsin装十sin'？ }g-'(B/4当4>0,B>0 82 tg-'(B/A,)们当A 0,B0,B'>0 g-(B'/A.'十x当A!'<0 85 u, ltg-'(B'/A,'十2x当A,>0,B'<0 16
GB/T13474-92 11 潮流计算 11.1分潮的选择 采用短期观测资料的分析结果,结合计算要求,选择不同的分潮组合,进行计算 常用分潮组合表 数 分 潮 代 分 M,S,K,0 M,S.,K,O,M,MS M,S.,K,O,N,K,P,Q1 M,s,,K,O,,N,,K.,P,Q,M,MS. 10 M.,s,,N,,K.,K,,O,P,Q,M,Ms,M 1 Q,O,M,x,P,K,p,J,0O,2N 22 N,哟,M,L,T,,s,,K;,MN,M.,MS,s 1.2正点溯流计算 结合具体计算站点,选定分潮数及其调和常数,确定计算的初始时间(t)及终止时间),增量为 h,应用正点潮流计算公式计算 正点潮流计算公式如下 U,=U,十 一, eac十w十以 86 v=v,+习v,cosC4+(V+u),一们 w,=(+vw) 87 x/2 若=0,v>0 |3x/2 若U，=0,V,<0 =tg-(V/U,若U,>0,V> 0 88 tg-1(v,/U,十其若U,<0 若u>o,W,<0 -1V/U.十2x [tg 11.3旋转流极值的计算 11.3.1旋转流极值的判别公式 89 Cw(e一w(,一1)Cw(e,)一w(,1))>0 一时的流速 式中w) w(-1)(t-1)时的流速 17
GB/T13474一92 w(十1) (e+1)时的流速 上式成立时说明在(,-1,6,+1)这个时间区间内有流速极值 1.3.2旋转流极值时间(M值的计算公式 wG二土w土 t(M=,十+ 90 wI一w半w 11.3.3计算旋转流极值,将已算得的()值代入正点潮流计算公式,即求得极值时刻t(M)对应的流 向、流速 11.4计算往复流的极值 11.4.1计算往复流主轴方向各分潮的调和常数,将各分潮的调和常数投影至往复流的主轴方向,得出 -组新的调和常数 U.sin.cosp十V.sinn.sing a,(c)=tg 十n其 91 UrocopVCsim C,(c)=[(w.sine.cosp十v.cosn.sinp)° 十(U.sine.cosp+v.sinn,sinp) (92 式中a,(e)-e分潮迟角在罗方向上的投影 c,e)c分潮的幅度在伊方向上的投影 11.4.2余流在往复流主轴方向上的投影 C,(c)=w,cos(Q,一p 93 当对余流情况不了解时,取C,(O)=0 1.43往复流每小时流速值we)值计算公式 c,e w(=C,()十 cos[o+(V +u),一a,(i) 94 144往复流极值时间判别公式与1l.3.1相同 145往复流极值流速w)值计算公式 w(dM一w(- (ex一,)'[w(一1一2w(e,)+w(十1)刀95 云 11.4.6往复流转流时刻的判断 第一判断式 a. w(e,)w(,+1)<0 96 18
GB/T13474-92 当相邻二小时流速满足该式时,表示在(o,t+1)间有流速零值 第二判断式 b w(e,w(en<0 97 表示在(o,tM)间有流速零值 11.4.7往复流流速零值时刻(0)值的计算 -W wG1- t(0=to十 TwC"一T-w竿wF功 w)土w十1)2] ICw二)一w土)]'一8w)CwG 98 士 IwwC半w觉 应用上式可算出两个t(0)值,应取满足不等式 0一t<1的值为上式的解 当相邻两整点的流速同号且它们之间有一极值,此极值的前后各有一转流时刻,则流速零值时刻的 两个解都有意义 12 半日潮流围的编制 2.1潮流类型 以主要全日分湖(K,O,)与太阴主要半日分潮(M)的潮流椭圆长半轴之比值确定潮流性质(潮流 椭圆要素的计算方法见第10章 W 99 w= L十w)/w w<0.5为半日潮流; 0.5GB/T13474一92 2.2.2.1平均流速流向w.值计算公式 w=(U'十v't (101 若U=0,V>0 x/2 3x/2 若U=0,v<0 tg-(v/)若U>0,V>0 102 - tg'(V/U)十x若u0,vGB/T13474一92 附 录 A 天文变量及非调和常数的有关计算公式 (补充件 天文变量的计算公式 A1 S =277".02+129".3848(y一1900)13".1764 D,+D十D十一A A1 个 A' =280°.19一0°.2387y一1900)十0°.9857 A2 D,十D十D十一A P=334.39+40".6625(y一1900)0'.1114 A3 [B,+马+n+式一A ”=100".84十19".3282(y一1900+0".0530 (A4 (2+马十马+式-A P =281".22十0".0172(y一1900)十0".00005 D,DD十一A A5 月球平均经度， 式中;S -太阳平均经度 h 月球近地点平均经度 P ' 白道升交点经度 -太阳近地点平均经度; P -年分; 从y年1月1日以后所经过的日数， D D -1900年至y年的国年数; -天文变量日期指标(潮汐分析及潮汐潮流推算时,取值为零); Dr 当;D,=0,表示计算观测日期零点的值 D一1,表示计算观测日期第二天零点的值 D=2,表示计算观测日期第三天零点的值 第一个潮位资料所对应的平太阳时(在潮流分析与潮汐潮流推算时取值为零) 21
GB/T13474一92 -潮流准调和分析潮令指标量(在潮汐分析、潮流分析的调和算法和潮汐潮流推算时取值为 零),它的取值法有三种 测流点附近有长期验潮资料的半日潮海区 A=1.837(g (A6 gN， 测流点附近有长期验潮资料的全日潮海区 b A=0.911(s Bo 测流点附近缺乏长期验潮资料的我国近海区 A=2.0 (A8 A2关于交点因数<)及天文相角(w)的计算公式 A2.1L分潮了及让的计算公式 fcos以=1一0.2505cos2P一0.1102eos(2P一"' 一0.0156cos(2P一2'一0.0370cosN" A9 fsinu=-0.2505sin2P一0.1102sin(2P一N' a.0156sin(2P一2w)-0.0370sin" (A10 A2.2M,分潮f及u的计算公式 fcosu=2cosP十0.4cos(P一N' (A11 Jfsinu=sinP+0.2sin(P一N' A12 表A1主要分潮的f及u计算公式 f的乘子系列 x的乘子系列 cosN os2N" cos3” sinN" sin2N" sin3N" M 1.000o -0.1300 0.0013 0.0o 0.00 0.00 Mr 1.0429 0.4135 -0.0040 -23.74 2.68 -0.38 O -0.0147 -1.34 1.0089 0.1871 0.001 10.80 0.19 22
GB/T13474-92 续表A1 f的乘子系列 的乘子系列 分 潮 cosN" cos2N” cos3” sinN sin2N sin3N 1.0060 0.115o 一0.0088 0.0006 8.86 0.68 -0.07 K 1.0129 0.1676 -0.0170 0.0016 -12.94 -0.19 1.34 oo 1.1027 0.6504 0.0317 -0.0014 -36.68 4.02 -0.57 M 1.0004 -0.0373 0.0002 一2.14 0.0082 -0.001 5 -0.04 K 0.2863 -17.74 1 .0241 0.68 A3分潮参数表 f,vu的参数表达见表A2 表A2f,v,丝参数表 号 角速度系数 分潮名称 A' “11”“63”|"16”"128" S. 0.0410686 s 0.0821373 S 0.1232095 S. 0.1642746 M. M 0 MS 4622374 0.5443747 M. M. M. SM 0.5854433 M. M. 一M M., MS 1.0158958 1.0980331 M M M 10 S.M 1.139101？ M M M.十M 11 Mm 1.6424078 MmXMn 90 0QK 2.3006278 0,一K O×K -5 12 2 O 12.8542862 270 O Q. 2OK 12.8450025 90 2O一K O×K， o. 270 2.9271398 1 o1 Q 0 10 13.3986609 270 0. 14 15 13.4715145 270 23
GB/T13474一92 续表A2 序 号 分潮名称 角速度系数 A 11"“63”|"116""128" O. O,s. 270 12 13.8608983 0 0 0S 13.9019669 270 13 一2 2 10 14 16 0 13.9430356 2 270 0 0 s.o. 270 O O 15 13.9841042 14.0251728 90 M; 16 MP M 1 11 17 12 18 M 14.4920521 90 M M x 90 13 18 19 14.5695476 -2 14.9178647 270 19 20 x 14 P 20 14.9589314 270 15 21 180 16 21 22 S 15.0000000 90 K 22 23 KK 5.040686 K 17 9g 15.0821353 18 23 24 g0 19 24 25 15.1232059 P1 20 5 26 15.5125897 90 J 6 27 J 15.585443 90 21 26 27 0 28 2PO. 15.9748272 270 16.056964s O 2 90 -O. SO 8 2 22 23 29 30 0O. 16.139101？7 90 0O. 0O. 2K-O Kxo. 30 2KQ. 16.6834764 270 2 2Ns 6.879 2M 4590 6 M 31 32 2NK.S. 26.9615963 -6 2 2M+K; M×K Of 331 33 24 0Q: 27.3416964 180 2O. MNs. 2M M 34 25 32 27,4238337 35 27.5059710 2M,+K; M×K MNK.S. 36 2MS,K27.8039338 2M;一2K M×K N M. h 37 26 33 27.8953548 24
GB/T13474一92 续表A2 号 序 分潮名称 角速度系数 pf “11” “63”“116”“128" 2M 34 27.8953548 M M 0. 2 35 27.886071l 180 2O. OH 27 36 38 27.9682084 M M NS 37 28.3986609 M M 28 38 39 28.4397295 N M. M 39 S,N 28.4807981 M M M 29 40 40 28.5125831 M: 4 -3 NKS 28.5218668 M+K; MXK 41 2KNS 28.6040041 -3 2K;+M K×M 30 42 42 OP 28.901967 0 180 O, O. 43 28.9430356 Ms. M M 2 43 2 31 44 M. 28.9841042 M M. -2 45 S.M: 29.0251728 M M 29.0662415 -2 32 6 4 MKs M+K; MXK 45 M.(KS).29.1483788 --2 M+2K M.×K 46 12SN(MK.29.3734880 -K; M×K2 47 33 47 29.4556253 18o M M 入 34 48 48 L 29.5284789 180 L L K 2 49 NKM 29.5377626 M×K2 5o 29.917862" -2 s,s 35 51 49 T 29.9589333 36 52 50 S 30.0000000 30.o41o66" 18o R 37 53 51 54 38 52 K 30.0821373 2 K K 39 55 53 MSN 30.5443747 M 18o KJ 30.6265119 K;+小 KxJ 40 56 54 25
GB/T13474--92 续表A2 日 V 序 分潮名称角速度系数 h “63”“116”“128” 552KM(SND.30.7086493 2K K×M 41 57 56 2SM 31.0158958 -M M" 31.0980331 57 sKM, K;一M. K,XM 58 59 O. 41.8291068 90 3O. O 58 NO 42.3827651 270 M+0 MXO 270 M+O M.xO. MO 42.9271398 60 42 59 2MK 42.9271398 270 2M-K MXK 180o -M. M!? 43 61 60 M, 43.4761563 . 270 O. so, 43.9430356 44 62 61 63 62 44.0251728 3 90 M+K M,XK, 5 MK K K 64 45.0410686 46 63 SK 90 K 27o 3K K 65 45.1232059 2O+M OXMs 66 0M56.8701754 67 64 N 56.8794590 2M M 3M. M 65 3MS. 56.9523127 57.4238337 2M MN Mi 66 447 68 67 MNKS 57.5059710 2M+K M×K 3 Ms. 2M. M 69 57.9271398 70 68 57.9682084 2M M室 M 48 71 S.M 58.0092771 2M M号 sN M M 49 72 69 58.4397295 58.5218668 K,+M K.XM 70 KN 73 MS 58.9430356 M M M Ms 10 50 74 71 MS 58.9841042 59.0251728 3 M S,.M M -2 75 76 72 MK 59.0662415 M+K M×K 51 26
GB/T13474一92 续表A2 号 序 分潮名称 角速度系数 “11”l“63”l“116”“128" A 18o .528 sL 4789 73 59. L 77 74 s 60.0000000 52 K K 53 78 75 SK 60.0821373 MNO. 8693 270 2M+O. Mxo. 76 71.366 79 77 71.9112440 270 2M+O. MxO, 2MO. 3M. M 9C 78 3MP 71.9933813 72.4649023 79 MNK 90 M.+K MxK 8o MSO, 72.9271398 270 M+O. M×0. 2M. Mi 80 2MP 72.9271398 270 2MK 3.0092771 90 M+K MxK 81 81 82 74.0251729 -2 90 M+K M×K MSK, 82 83 3KM 74.1073101 90K+K+MaK×K×M. M,+K MxK 84 3NKS 85,4013258 8 85.8645632 85 2NM. 3M M 86 2NMKS. 85.9457005 3M;+K; M×K 3M. M 54 83 87 2MN 86.4079379 2MNKs. 86.4900752 8 M+K MXK M 55 11 84 89 M 86.9523127 3M 2M. M 56 85 90 MSN 87.4238337 MKN 87.5059710 91 2M,+K MXK 53 86 92 2MS 87.9682084 2M M 2M+K 58 87 93 2MK 88.0503457 M×K; 88.5218668 NSK M+K M,×K. 94 88 59 95 2SM 88.9841042 M M; 60 89 96 MSK 89.0662415 M+K M×K; 97 90.0000000 S 27
cB/T13474一92 续表A2 序 号 V 分潮名称 角速度系数 P “63”|“116”“128" 1” 98 2MNO. 100.3509735 270 3M+O MxO. 3MO, 3M+O. 90 100.8953483 270 M×O. 99 2NMK， 100.9046318 90 3M+K MxK 91 10G 2MSO 101.9112440 270 2M+O MxO. 3MK， 92 101.9933813 90 3M+K M×K 103.0092771 90 93 2MSK， MM+ M×K 101 MSKO 103.0092771 270M+K+OM×KxO. 2MN), 102 114.8476674 -10 4M: M金 94 103 115.3920422 MN M., Mt 104 3MNKS115,4741795 10 4M+K， M×K 61 95 105 M 115.9364169 4M: M 106 2MSN 116.407938o 96 M., M 2MNK 116.4900752 3M十K 107 M×K 97 108 3MSn 116.9523127 3M M 109 3MK117.0344500 3M,十K; M×K 110 MSNK，117.505971o 一5 2M+K M×K; 98 111 2MS). 17.9682084 2M M3 122MsK 2M+K 99 118.0503457 M×K 113 2M,NK， 129.8887360 -10 90 4M;+K Mt×K 130.4331108 M×K 114 3MNK， 90 4M+K 3MsO,130.8953483 2703M+O. 100 M×O 115 130.9774855 Mt×K 4MK 90 4M;+K 101 2MSO，131.9112441 270 2M+O. MxO 3MsK 02 116 131.9933813 90 3M;+K M×K 03 90 2M,SK，133.0092771 2M+K Mi×K 117 4MNo 144.3761464 1 10 5M， M 28
GB/T13474一92 续表A2 号 序 角速度系数 分潮名称 “11" “63”|“116”“128" 4.8208a11 -1o 10 5M M M 62 118 119g 3MSN 145.3920422 M M 104 4M Mt 105 120 4MSo0 145.9364169 .4079380 2M,sN146. 3M. M 106 3M+K， M×K 1212MSNK146.4900752 M 107 3M,S0 146.9523127 122 3M 270 MxO. Mso M;+O, -10 108 159.8794525 109 3M.SO160.8953483 270 3M+O. M×O, M,十K MtxK 110 123 4MSK160.9774855 90 90 .983381 3M,sK161. 3M,+K MxK 111 173.9046253 -12 12 6M M 124 Ms 63 M 112 125 4MSN174.3761464一11 10 5M 1o 5M. M 174.9205211 -10 113 126 5MSnz M M 114 3M.,SN175.3920423 M+K 127 3MSNKa175.4741795 -9 10 MXK M. Mt 115 128 4M.S175.9364169 175.9523128 -8 6 3M M 16 3M,Su A4潮流调和分析天文变量的计算公式 A4.1计算入,C/r.,A,C/r,I,,s 入=A'十0.0335sin(h'一P')十0.0004sin2(h'一P') A13 A14 C,/r，=1+0.0168cos(h'P'+0.0003cos2(A'一P') 29