一种海域理论最低潮面计算方法、系统、设备及介质与流程

1.本发明涉及海洋测绘技术领域，特别是涉及一种海域理论最低潮面计算方法、系统、设备及介质。


背景技术：

2.理论最低潮面是我国法定的海图深度基准面，通常通过建设验潮站并进行一年以上连续潮位观测后，通过调和分析方法得到13个分潮的调和常数，再利用弗拉基米尔法计算理论最低潮面。
3.然而，现有的理论最低潮面计算方法往往直接将长周期分潮的贡献直接折入短周期潮汐变化中，使得长周期分潮对最低潮的贡献改正存在较大的随机性和不确定性，导致理论最低潮面的计算结果准确度有限。


技术实现要素：

4.针对上述技术问题，本发明提供一种海域理论最低潮面计算方法，通过分别计算短周期分潮和长周期分潮的最低潮绝对值，综合得到二者对理论最低潮面的贡献，从而提高理论最低潮面计算结果的准确度。
5.第一方面，本发明提供一种海域理论最低潮面计算方法，包括：
6.获取预设时间段内的潮位观测数据；
7.通过潮汐调和分析法对所述潮位观测数据进行调和分析，得到十三个分潮的调和常数；所述十三个分潮包括天文分潮、浅水分潮和长周期分潮；
8.分别计算所述天文分潮与所述浅水分潮的高频分潮潮高最低潮面，以及所述长周期分潮的低频分潮潮高最低潮面；
9.根据所述高频分潮潮高最低潮面与所述低频分潮潮高最低潮面的和，确定理论最低潮面。
10.可选的，所述计算所述长周期分潮的低频分潮潮高最低潮面，具体为：
11.基于所述长周期分潮的调和常数以及所述长周期分潮中的sa分潮相角的预设步长，计算所述sa分潮周期内所述长周期分潮的最低潮面绝对值，得到低频分潮潮高最低潮面。
12.可选的，所述计算所述sa分潮周期内所述长周期分潮的最低潮面绝对值的过程中，所述长周期分潮中的s
sa
分潮相角通过下式得到：
[0013][0014]
其中，表示所述sa分潮相角，表示所述sa分潮的迟角，表示所述s
sa
分潮的迟角。
[0015]
可选的，所述计算所述天文分潮与所述浅水分潮的高频分潮潮高最低潮面，具体为：
[0016]
基于所述天文分潮与所述浅水分潮的调和常数，以及所述天文分潮中的k1分潮相角的预设步长，计算所述k1分潮周期内所述天文分潮与所述浅水分潮组合的最低潮面的绝对值，得到高频分潮潮高最低潮面。
[0017]
第二方面，本发明还提供一种海域理论最低潮面计算系统，包括：
[0018]
数据获取单元，用于获取预设时间段内的潮位观测数据；
[0019]
第一计算单元，用于通过潮汐调和分析法对所述潮位观测数据进行调和分析，得到十三个分潮的调和常数；所述十三个分潮包括天文分潮、浅水分潮和长周期分潮；
[0020]
第二计算单元，用于分别计算所述天文分潮与所述浅水分潮的高频分潮潮高最低潮面，以及所述长周期分潮的低频分潮潮高最低潮面；
[0021]
第三计算单元，用于根据所述高频分潮潮高最低潮面与所述低频分潮潮高最低潮面的和，确定理论最低潮面。
[0022]
可选的，所述第二计算单元具体用于：
[0023]
基于所述长周期分潮的调和常数以及所述长周期分潮中的sa分潮相角的预设步长，计算所述sa分潮周期内所述长周期分潮的最低潮面绝对值，得到低频分潮潮高最低潮面。
[0024]
可选的，所述计算所述sa分潮周期内所述长周期分潮的最低潮面绝对值的过程中，所述长周期分潮中的s
sa
分潮相角通过下式得到：
[0025][0026]
其中，表示所述sa分潮相角，表示所述sa分潮的迟角，表示所述s
sa
分潮的迟角。
[0027]
可选的，所述第二计算单元具体还用于：
[0028]
基于所述天文分潮与所述浅水分潮的调和常数，以及所述天文分潮中的k1分潮相角的预设步长，计算所述k1分潮周期内所述天文分潮与所述浅水分潮组合的最低潮面的绝对值，得到高频分潮潮高最低潮面。
[0029]
第三方面，本发明提供一种数据处理设备，包括处理器，所述处理器和存储器耦合，所述存储器存储有程序，所述程序由所述处理器执行，使得所述数据处理设备执行第一方面所述的海域理论最低潮面计算方法。
[0030]
第四方面，本发明还提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时，实现如上述第一方面所述的海域理论最低潮面计算方法。
[0031]
相比现有技术，本发明的有益效果在于：
[0032]
本发明提供的海域理论最低潮面计算方法基于两步极小值叠加法，分别计算天文分潮与浅水分潮的高频分潮潮高最低潮面，以及长周期分潮的低频分潮潮高最低潮面，再对二者进行求和，以将短于日周期的海面变化折叠于长周期的长期海面变化，从而消除了潮位观测数据计算得到的长周期分潮调和常数的误差，有效提高理论最低潮面计算精度，确保得到准确度更高的理论最低潮面。
附图说明
[0033]
为了更清楚地说明本发明的技术方案，下面将对实施方式中所需要使用的附图作
简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0034]
图1是本发明实施例提供的海域理论最低潮面计算方法的流程示意图；
[0035]
图2是本发明实施例提供的海域理论最低潮面计算系统的结构示意图。
具体实施方式
[0036]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0037]
应当理解，在本发明说明书中所使用的术语仅仅是出于描述特定实施例的目的而并不意在限制本发明。如在本发明说明书和所附权利要求书中所使用的那样，除非上下文清楚地指明其它情况，否则单数形式的“一”、“一个”及“该”意在包括复数形式。
[0038]
术语“包括”和“包含”指示所描述特征、整体、步骤、操作、元素和/或组件的存在，但并不排除一个或多个其它特征、整体、步骤、操作、元素、组件和/或其集合的存在或添加。
[0039]
术语“和/或”是指相关联列出的项中的一个或多个的任何组合以及所有可能组合，并且包括这些组合。
[0040]
理论最低潮面在中国被规定作为海域的深度基准面，现有的方法通常采用弗拉基米尔法计算理论最低潮面。然而，现有方法在利用海面的年周期变化和半年周期变化修正一天内的潮高短期变化过程中，由于在短期变化的极小值状态时，长周期变化可能处于最大值附近，这使得长周期变化的影响在不同地点存在随机性，从而影响了理论最低潮面确定的准确性。
[0041]
对此，本发明一个实施例提供了一种海域理论最低潮面计算方法，能够有效提高理论最低潮面值的计算准确度。
[0042]
请参照图1，第一方面，本发明实施例所提供的海域理论最低潮面计算方法具体包括下述步骤。
[0043]
s1：获取预设时间段内的潮位观测数据。
[0044]
需要说明的是，对理论最低潮面的计算需采集至少一年以上的潮位观测数据，故预设时间段应大于等于一年。
[0045]
s2：通过潮汐调和分析法对所述潮位观测数据进行调和分析，得到十三个分潮的调和常数；所述十三个分潮包括天文分潮、浅水分潮和长周期分潮。
[0046]
具体地，十三个分潮包括八个天文分潮、三个浅水分潮m4、ms4、m6和两个长周期分潮sa、s
sa
；其中，所述八个天文分潮包括四个半日分潮q1、o1、p1、k1和四个全日分潮n2、m2、s2、k2。
[0047]
可以理解的是，所述调和常数又称为分潮调和常数，具体为各分潮的平均振幅和迟角。
[0048]
s3：分别计算所述天文分潮与所述浅水分潮的高频分潮潮高最低潮面，以及所述长周期分潮的低频分潮潮高最低潮面。
[0049]
在本实施例中，可根据长周期分潮的调和常数，计算在一年周期分潮sa的一个周
期内，一年周期分潮sa与半年周期分潮s
sa
的最低潮面的绝对值，得到低频分潮潮高最低潮面。
[0050]
同时，基于天文分潮与浅水分潮的调和常数以及预设步长，计算k1分潮周期内的天文分潮与浅水分潮组合的最低潮面的绝对值，得到高频分潮潮高最低潮面。
[0051]
s4：根据所述高频分潮潮高最低潮面与所述低频分潮潮高最低潮面的和，确定理论最低潮面。
[0052]
需要说明的是，每一个分潮实质上为余弦形式简谐振动，通常以英文字母或英文字母组合为分潮名，各主要分潮名在海洋潮汐学中均已统一约定，对应引潮力位调和展开式的具体一项，其下标表示振动频次或振动周期。
[0053]
在主要分潮中，下标为1的各分潮q1、o1、p1、k1，其相角变化一个周期用时接近1天，归于全日分潮族。
[0054]
下标为2的各分潮n2、m2、s2、k2，其相角变化一个周期接近半天，或分潮一天内约完成两个周期振动，归于半日分潮族。
[0055]
下标为4和6的三项主要振动m4、ms4、m6分别表示一天内变化4个周期和6个周期的分潮，其主要由浅水摩擦效应引起，分别归于1/4日和1/6日分潮族，并根据其诱因称为浅水分潮。
[0056]
下标为a和sa的两个分潮sa、s
sa
对应周期分别为一年和半年，归类为长周期分潮。
[0057]
需要说明的是，作为简谐振动，每一分潮的潮高变化过程可表示为：
[0058][0059]
式中，h(t)表示随时间变化的潮高，表示随时间以对应角速率σ线性变化的相角，其中为分潮周期；v(t)表示天体引潮作用相应分量的相角，v0则为v(t)选定参考时刻的分潮相位，即分潮初相位。h、g分别为分潮振动的幅度和迟角，合称为分潮调和常数。
[0060]
可以理解的是，利用多个分潮求任一时刻潮高或求特征潮高的过程，实际上是各分潮变化过程的线性叠加。
[0061]
在每一分潮族中，各分潮的相角变化周期接近，分潮周期决定于角速率或频率，各分潮周期是确定的，可根据太阳和月球轨道运动的有关天文变量进行计算。需要说明的是，对应特定时刻，也可由一组天文变量计算分潮相角。
[0062]
在本实施例中，全日潮族和半日潮族中分别选定的4个分潮，不同分潮之间的角速率差大于月球公转角速率或地球绕日公转角速率，每一分潮在同族内均存在分群现象，但根据一年的观测数据可以分析计算出对应的分潮调和常数。
[0063]
可以理解的是，在一年内两个长周期分潮分别变化一个周期和半个周期，其调和常数同样可由一年的观测数据计算得到。
[0064]
需要说明的是，根据分潮角速率排列，在每一主要分潮附近还存在有量值不可忽略的振动项：小分潮。在海洋潮汐研究和分析中，将分潮角速率差异小于地球绕日公转运动角速率的所有分潮称为分潮亚群。
[0065]
受限于信号分解rayleigh准则的限制，用一年或更短时段的潮汐观测数据并不能有效分辨同一亚群的分潮，故通常将同一亚群分潮的贡献合并于其中的主要分潮，合并过
程中所产生的振幅修正项被称为交点(订正)因子，记为f，f的主要周期为18.613年，即月球轨道面进动周期或月球升交点周期。
[0066]
因此，本实施例考虑亚群内分潮合并后，可通过下式得到分潮亚群分潮潮高：
[0067][0068]
可以理解的是，多分潮潮高的组合实际为所在分潮亚群潮高的组合。
[0069]
可以理解的是，全日分潮和半日分潮之间存在引潮力相角的相互联系，即对于任一时刻t，存在：
[0070][0071][0072][0073][0074]
考虑海洋对相应频率的响应延迟(即迟角)，对应分潮的实际相角之间满足下述关系：
[0075][0076][0077][0078][0079]
即：
[0080][0081][0082][0083][0084]
为了简洁表达，可将有关迟角和常数项组合为常数符号。
[0085]
四个半日分潮的相角表达式均与分潮k1的相角有关，可实现m2与o1、s2与p1以及n2与q1三对分潮的配对，以及k2与k1的单独配对。即：将6个分潮组合潮高配对为3个变振幅和变相角的振动项，并用变量r表示这类组合振动的振幅。
[0086]
对应的，对于前三个配对分潮，其组合潮高具有相同的形式：
[0087][0088]
进一步地，可将上式变换为如下形式的变振幅、变迟角振动：
[0089]
其中：依赖于k1分潮的时变振幅和时变相角分别表示为：
[0090][0091][0092]
在本实施例中，每一对合成潮高的振幅和迟角均为k1分潮相角的函数，为了简化计算，可将三对分潮的组合潮高分别标号为1、2和3，并略去自变量标识，则对根号内关系进行整合后得到：
[0093][0094][0095][0096][0097][0098][0099]
因此，在考虑全日分潮和半日分潮共8个分潮情况下，组合潮高随时间的变化可表示为:
[0100][0101]
进一步地，可对合成潮高的3个分潮取最小值形式，即令：
[0102][0103]
则面向求潮高最小值的变化过程可简化为：
[0104][0105]
上式右端的各部分均退化为k1分潮相角的函数，因此，可在该分潮一个周期内(23.93小时，约1天)，通过预设的角度步长计算低潮潮高，以求得最小值。
[0106]
由于得到的最小值结果为负值，需取绝对值作为8个分潮的深度基准值。
[0107]
进一步地，本实施例还考虑了以更快角速率变化的浅水分潮的潮高贡献。由于浅水分潮主要由半日分潮波与海底的摩擦效应产生，这些主要分潮即为浅水分潮的源分潮，根据浅水分潮的来源，浅水分潮的理论相角和实际相角间分别存在下述关联：
[0108][0109][0110][0111]
[0112][0113][0114]
设定源分潮的相角分别为极小值形式的量值ε1 180
°
和ε2 180
°
，则三个浅水分潮的相角可通过源分潮间接与k1分潮的相角变化建立联系：
[0115][0116][0117][0118]
因此，本实施例考虑全日分潮、半日分潮和浅水分潮共11个分潮贡献的求解高频分潮潮高最低潮面的计算方式具体为：
[0119][0120]
可以理解的是，浅水分潮相对其源分潮量值较小且周期短，因此对最低潮高的求解有合理的调节作用。
[0121]
进一步地，对于长周期分潮的贡献，由于半日分潮和浅水分潮的周期为全日分潮k1的2倍至6倍，因此利用上述11个分潮计算高频分潮潮高最低潮面时，将短周期高频率振动的贡献折入长周期低频率的振动是合理的。
[0122]
另一方面，将潮高根据时间变化归算为依据单一分潮的相角变化，不仅可避免时段的选择，且根据k1分潮的相角在一个周期内变化时，按设定的步长取得的最小值进行计算，还可有效保证计算效率。
[0123]
然而，将两个长周期分潮的贡献折入周期约一天的潮汐变化则具有较大的不合理性。尽管两个长周期分潮的相角与s2分潮相角存在有形式上的关联，但考虑s2分潮已与p1分潮进行合成，并采用了其综合潮高极小值状态的相角进行取值，进一步考虑各有关分潮的实际迟角时，将使得长周期分潮对最低潮的贡献改正在不同的地点时产生较大的随机性和不确定性。
[0124]
因此，本发明实施例将两个长周期分潮sa、s
sa
的潮高进行单独组合，在sa分潮变化的一个周期内计算长周期分潮的综合最低潮的绝对值，以得到低频分潮潮高最低潮面。
[0125]
具体计算方式如下式：
[0126][0127]
进一步地，结合高频和低频的潮汐成分对理论最低潮面的贡献，本实施例将上述低频分潮潮高最低潮面与高频分潮潮高最低潮面进行求和计算，得到最终的理论最低潮面数：
[0128]
l＝l
hf
l
lf
[0129]
可以理解的是，周期为一日以内的短周期海面变化基本每日发生的，最多顾及了
每月的大小潮变化，而周期为一年或半年的长周期变化本不受控于短周期振动。因此，本发明实施例考虑了海面的季节性周期变化规律，将日周期振动现象定位于季节性周期变化之中，从而消除因每一时段潮汐观测数据计算求得的长周期分潮调和常数误差，特别是其迟角误差与s2分潮迟角耦合过程中的随机影响，从而有效提高理论最低潮面值计算的精度。
[0130]
本发明上述实施例提供的海域理论最低潮面计算方法基于两步极小值叠加法，分别计算了天文分潮与浅水分潮的高频分潮潮高最低潮面，以及长周期分潮的低频分潮潮高最低潮面，再对二者进行求和，以将短于日周期的海面变化折叠于长周期的长期海面变化，从而消除了潮位观测数据计算得到的长周期分潮调和常数的误差，有效提高理论最低潮面计算精度，确保得到准确度更高的理论最低潮面结果。
[0131]
请参照图2，第二方面，本发明另一实施例提供一种海域理论最低潮面计算系统，包括数据获取单元101、第一计算单元102、第二计算单元103和第三计算单元104。
[0132]
数据获取单元101用于获取预设时间段内的潮位观测数据。
[0133]
第一计算单元102用于通过潮汐调和分析法对所述潮位观测数据进行调和分析，得到十三个分潮的调和常数；所述十三个分潮包括天文分潮、浅水分潮和长周期分潮。
[0134]
第二计算单元103用于分别计算所述天文分潮与所述浅水分潮的高频分潮潮高最低潮面，以及所述长周期分潮的低频分潮潮高最低潮面。
[0135]
第三计算单元104用于根据所述高频分潮潮高最低潮面与所述低频分潮潮高最低潮面的和，确定理论最低潮面。
[0136]
上述系统内的各单元之间信息交互、执行过程等内容，由于与本发明的海域理论最低潮面计算方法实施例基于同一构思，具体内容可参见本发明方法实施例中的叙述，此处不再赘述。
[0137]
第三方面，本发明提供一种数据处理设备，包括处理器，所述处理器和存储器耦合，所述存储器存储有程序，所述程序由所述处理器执行，使得所述数据处理设备执行第一方面所述的海域理论最低潮面计算方法。
[0138]
第四方面，本发明还提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时，实现如上述第一方面所述的海域理论最低潮面计算方法。
[0139]
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于计算机可监听存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体(read-only memory，rom)或随机存储记忆体(random access memory，ram)等。
[0140]
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也视为本发明的保护范围。