进一步依球面三角形余弦定理得，
[0050] cosAΩ=cosN·cosv+sinN·sinv·cosω (10)
[0051]因此，
[0052]
(11)
[0053] 月球在白道上从A点起算的平均经度σM为，
[0054] σΜ=ΑΩ+s-N (12)
[0055]
[0056] 其中eM为月球轨道偏心率即0.054900,m是太阳与月球运行速度的比值，即 0.074804,其余参数见式（1)~（5)所示。至此，依式（6)、（13)可计算得到δM值。
[0057] 步骤2、根据日潮特征海域的Kp两个分潮迟角计算日潮龄值。具体方法如下：
[0058] 日潮龄t可采用下式得到，其中絳、：?分别表示V分潮迟角。
[0059]
[0060] 步骤3、将步骤1得到的月赤炜高度角值作时间尺度因子为步骤2得到的日潮龄值 的平移，得到顾及日潮龄值的月赤炜高度角值。
[0061] 步骤4、将潮汐预报时间间隔为一小时、预报周期为18. 61年的逐时潮汐预报值与 步骤3得到的经时间平移后的逐时月赤炜高度角值作精密配准，以一个回归月为单位，统 计计算多组回归潮高高潮位，最后计算回归潮平均高高潮位。具体方法如下：
[0062] 在本步骤中，潮汐预报时间间隔为一小时、预报周期为18. 61年的逐时潮汐预报 值采用以下方法得到：从当地平均海面起算的任意时刻t潮汐预报值h(t)，采用如下计 算：
[0063]
[0064]式中私和gi表示分潮振幅和迟角，fJPui分别表示分潮交点因数和交点改正角， σi表示分潮角速率，v1表示分潮天文初位相，m为13。
[0065] 在得到逐时潮汐预报值后，将逐时潮汐预报值与逐时月赤炜高度角值作精密配 准，以一个回归月为单位，统计计算多组回归潮高高潮位。假定一个验潮站TideGauge共13 个分潮调和常数如表1所示，由此知该验潮站潮汐类型数为2. 40,所在海域为不正规日潮。
[0066]表1TideGauge的13个分潮的调和常数
[0067]
[0068] 根据日潮龄值t= 44小时，结合18. 61年的月赤炜高度角值，作时间平移配准；进 行以一个回归月为基本单位的回归潮高高潮位统计；对18. 61年周期内所有回归潮高高潮 位计算平均值，其中顾及日潮龄值对月赤炜高度角与回归潮高高潮位作时间配准的结果如 图3所示。图3给出了TideGauge验潮站部分时段（2015-01-01~2015-01-25)的潮汐预 报值与月赤炜值之间的关系，图中水平方向的浅灰色曲线表示月赤炜值，水平方向的深灰 色曲线表示顾及日潮龄作时间平移后的月赤炜值，单位均为度；垂直方向曲线表示潮汐预 报值，单位为cm。
[0069] 经计算，TideGauge验潮站的回归潮平均高高潮位值为77. 2厘米，即该站的海岸 线高度值为77. 2厘米。
[0070] 本发明的回归潮平均高高潮位与传统海岸线定义中的多年平均大潮高潮面在基 本概念、计算原则及方法上均不同，这种差异在不同潮汐类型海域是真实存在且不可替换 的，将海岸线定义修订为平均大潮高潮面或回归潮平均高高潮位与海岸的交接线更科学。 本发明适用于在日潮特征海域开展的海岛礁测绘海岸线高度值计算。
[0071] 需要强调的是，本发明所述的实施例是说明性的，而不是限定性的，因此本发明包 括并不限于【具体实施方式】中所述的实施例，凡是由本领域技术人员根据本发明的技术方案 得出的其他实施方式，同样属于本发明保护的范围。
【主权项】
1. 一种基于月赤缔高度角的回归潮平均高高潮位的计算方法，其特征在于包括W下步 骤： 步骤1、W-个回归月为计算单位，从潮软数据开始时刻起算，逐时精确计算一个完整 潮软运动周期的逐时月赤缔高度角值； 步骤2、根据日潮特征海域的两个分潮迟角计算日潮龄值； 步骤3、将步骤1得到的月赤缔高度角值作时间尺度因子为步骤2得到的日潮龄值的平 移，得到顾及日潮龄值的月赤缔高度角值； 步骤4、将潮软预报时间间隔为一小时、预报周期为一个完整潮软运动周期的逐时潮软 预报值与步骤3得到的经时间平移后的月赤缔高度角值作精密配准，W-个回归月为单位 统计计算多组回归潮高高潮位，最后计算回归潮平均高高潮位。2. 根据权利要求1所述的一种基于月赤缔高度角的回归潮平均高高潮位的计算方法， 其特征在于：所述一个回归月为27. 32158天，所述一个完整潮软运动周期为18. 61年。3. 根据权利要求1或2所述的一种基于月赤缔高度角的回归潮平均高高潮位的计算方 法，其特征在于：所述步骤1计算一个完整潮软运动周期的逐时月赤缔高度角值是通过月 球在白道上的逐时真经度和逐时白赤交角精密计算得到的，具体计算公式为：式中：SM表示月赤缔高度角值，Im为从T'起算的月球在白道上的逐时真经度，！表 示白赤交角；S为从T'起算的月球平均经度，h为从T起算的太阳平均经度，P为从T' 起算的月球近地点平均经度，Pg为从T起算的太阳近地点平均经度，N为从T起算的月球 升交点平均经度；T为春分点、T'为辅助春分点，OM为月球在白道上从A点起算的平均 经度，A为白赤交点；Gm为月球轨道偏屯、率，其取值为0. 054900，m是太阳与月球运行速度的 比值，其取值为0.074804。4. 根据权利要求1或2所述的一种基于月赤缔高度角的回归潮平均高高潮位的计算方 法，其特征在于：所述步骤2计算日潮龄值的方法为： 《=〇.(川(g、,-瑟巧） 式中：斯、扣,分别表示Ki、O汾潮迟角，t为日潮龄值。5. 根据权利要求1或2所述的一种基于月赤缔高度角的回归潮平均高高潮位的计算方 法，其特征在于：所述步骤4逐时潮软预报值的计算方法为：式中：h(t)为从当地平均海面起算的任意时刻t潮软预报值，Hi和gi表示分潮振幅和 迟角，和Ui分别表示分潮交点因数和交点改正角，O1表示分潮角速率，Vi表示分潮天文 初位相，m为13。
【专利摘要】本发明涉及一种基于月赤纬高度角的回归潮平均高高潮位的计算方法，其主要特点是包括：计算一个完整潮汐运动周期的逐时月赤纬高度角值；根据日潮特征海域的两个分潮迟角计算日潮龄值；计算顾及日潮龄值的月赤纬高度角值；将逐时潮汐预报值逐时月赤纬高度角值作精密配准，统计计算多组回归潮高高潮位及回归潮平均高高潮位的计算。本发明实现了日潮特征海域基于月赤纬高度角的回归潮平均高高潮位计算功能，弥补了已有海岸线定义不适用于日潮特征海域的缺陷，拓展了海岸线定义的内涵，适用于日潮特征海域进行的海岸线高度测量。
【IPC分类】G06F17/50
【公开号】CN105373684
【申请号】CN201510955712
【发明人】黄辰虎, 陆秀平, 吴太旗, 欧阳永忠, 边刚, 黄贤源, 邓凯亮, 王耿峰
【申请人】中国人民解放军92859部队
【公开日】2016年3月2日
【申请日】2015年12月16日