留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

利用时变潮时差进行水位改正

刘聚 暴景阳 许军

刘聚, 暴景阳, 许军. 利用时变潮时差进行水位改正[J]. 武汉大学学报 ● 信息科学版, 2019, 44(5): 675-681, 689. doi: 10.13203/j.whugis20170302
引用本文: 刘聚, 暴景阳, 许军. 利用时变潮时差进行水位改正[J]. 武汉大学学报 ● 信息科学版, 2019, 44(5): 675-681, 689. doi: 10.13203/j.whugis20170302
LIU Ju, BAO Jingyang, XU Jun. Water Correction Method Based on Temporal Variation of Tidal Time Difference[J]. Geomatics and Information Science of Wuhan University, 2019, 44(5): 675-681, 689. doi: 10.13203/j.whugis20170302
Citation: LIU Ju, BAO Jingyang, XU Jun. Water Correction Method Based on Temporal Variation of Tidal Time Difference[J]. Geomatics and Information Science of Wuhan University, 2019, 44(5): 675-681, 689. doi: 10.13203/j.whugis20170302

利用时变潮时差进行水位改正

doi: 10.13203/j.whugis20170302
基金项目: 

国家自然科学基金 41074002

国家自然科学基金 41501500

详细信息
    作者简介:

    刘聚, 博士生, 主要从事海道测量水位控制理论与应用方面的研究。liuju9752@163.com

    通讯作者: 暴景阳, 博士, 教授。jingyangb@sina.com
  • 中图分类号: P229

Water Correction Method Based on Temporal Variation of Tidal Time Difference

Funds: 

The National Natural Science Foundation of China 41074002

The National Natural Science Foundation of China 41501500

More Information
    Author Bio:

    LIU Ju, PhD candidate, specializes in water level control in hydrographic survey. E-mail:liuju9752@163.com

    Corresponding author: BAO Jingyang, PhD, professor. E-mail:jingyangb@sina.com
  • 摘要: 针对目前实施时差法水位改正时使用的潮时差未考虑时变现象的影响以及利用潮时差确定水位观测时段选择混乱等问题,提出利用逐时滑动时段确定潮时差以进行水位改正的方法,分析了潮时差时变现象的原因及潮时差逐时滑动时段与逐日独立时段确定结果的差异,推导了时差法水位改正的误差公式。以某海域验潮站为例,进行了水位改正的数据实验,比较了逐时滑动时段与逐日独立时段确定的潮时差的水位改正效果。实验结果表明,利用逐时滑动确定的潮时差进行水位改正,在理论与精度方面均对目前使用的时差法有所完善与提高,尤其在潮差差异较大的海域,改进效果更加显著。
  • 图  1  逐日独立时段与逐时滑动时段划分示意图

    Figure  1.  Daily and Hourly Time-Period-Selecting Methods

    图  2  某海域各验潮站分布图

    Figure  2.  Distribution of Tidal Stations

    图  3  随从分潮对主分潮调制作用

    Figure  3.  Minor Constituent's Modulation Effect on Main Constituent

    图  4  3种时段潮时差确定结果比较图

    Figure  4.  Comparison of Time Difference Results of Various Time-Period-Selecting Methods

    图  5  水位曲线导数差异

    Figure  5.  Derivative Variations of Water Level Curve

    表  1  验潮站主要分潮调和常数及其他实验数据

    Table  1.   Experiment Data of Tidal Stations

    分潮 A B 单分潮时差/min
    H/cm g/(°) H/cm g/(°)
    O1 29.5 237.9 30.7 243.5 -24.1
    K1 36.8 313.9 38.6 319.9 -23.9
    M2 44.0 274.0 50.1 290.3 -33.7
    S2 17.3 329.3 19.2 347.0 -35.4
    下载: 导出CSV

    表  2  逐日独立时段与逐时滑动时段潮时差差异统计/min

    Table  2.   Variations of Tidal Time Differences Determined by Daily and Hourly/min

    差异参数 基准站-待定站
    B-I I-H E-I D-I A-I F-I C-B B-A B-D C-D F-H E-F F-G
    中误差 2.54 1.67 0.95 0.93 2.34 0.86 2.32 1.90 1.58 3.98 0.79 1.03 1.32
    变幅 25.9 14.1 9.4 9.9 7.2 8.1 42.4 18.6 16.8 60.0 7.0 6.9 8.4
    最大 10.9 6.0 3.4 3.1 4.8 2.6 14.0 6.4 5.8 18.1 2.8 2.9 3.1
    平均 5.4 4.3 1.3 2.3 3.6 1.8 10.4 4.0 4.1 14.3 1.5 1.7 1.8
    下载: 导出CSV

    表  3  各组实验水位改正精度统计/cm

    Table  3.   Precise of Water Correction Experiments/cm

    精度指标 (1)AIH (2)ABC (3)DEH (4)EFG
    max 11.8 14.7 22.3 15.7 21.2 26.4 9.8 11.3 12.3 5.6 5.4 9.7
    $\hat \sigma $ 1.1 1.8 10.8 2.7 5.8 12.3 0.9 1.8 2.6 1.1 1.0 2.1
    下载: 导出CSV
  • [1] 暴景阳, 刘雁春.海道测量水位控制方法研究[J].测绘科学, 2006, 31(6): 49-51 doi:  10.3771/j.issn.1009-2307.2006.06.013

    Bao Jingyang, Liu Yanchun. A Study of Water Level Control for Hydrographic Survey [J]. Science of Surveying and Mapping, 2006, 31(6): 49-51 doi:  10.3771/j.issn.1009-2307.2006.06.013
    [2] 刘雁春.海洋测深空间结构及其数据处理[M].北京:测绘出版社, 2002

    Liu Yanchun. Space Structure and Data Processing in Marine Sounding [M]. Beijing: Surveying and Mapping Press, 2002
    [3] 许军.水下地形测量的水位改正效应研究[D].大连: 海军大连舰艇学院, 2009

    Xu Jun. Research on the Effects of Water Level Correction on Underwater Topographic Survey [D]. Dalian: Dalian Naval Academy, 2009
    [4] Lea S H. Hydrographic Surveying Methods, Tables and Forms of Notes [M]. Canada: Andesite Press, 2015
    [5] Lachapelle G. Hydrographic Surveying [D]. Canada: The University of Calgary, 1996
    [6] 暴景阳, 许军.海道测量水位控制的技术体系及标准更新[J].海洋测绘, 2016, 36(6):1-6 doi:  10.3969/j.issn.1671-3044.2016.06.001

    Bao Jingyang, Xu Jun. Technical Framework and Specification Revision for Water Level Control in Hydrographic Survey [J]. Hydrographic Surveying and Charting, 2016, 36(6):1-6 doi:  10.3969/j.issn.1671-3044.2016.06.001
    [7] 谢锡君, 翟国君, 黄谟涛.时差法水位改正[J].海洋测绘, 1988, 8(3):22-26 http://d.old.wanfangdata.com.cn/Conference/5903555

    Xie Xijun, Zhai Guojun, Huang Motao. Water Correction in Time-Difference Method [J]. Hydrographic Surveying and Charting, 1988, 8(3):22-26 http://d.old.wanfangdata.com.cn/Conference/5903555
    [8] 刘雁春, 暴景阳.海道测量水位改正数学模型[J].海洋测绘, 1992, 12(4): 17-22 http://www.cnki.com.cn/Article/CJFDTOTAL-HYCH201104003.htm

    Liu Yanchun, Bao Jingyang. Mathematical Model of Water Correction in Hydrographic Surveying [J]. Hydrographic Surveying and Charting, 1992, 12(4):17-22 http://www.cnki.com.cn/Article/CJFDTOTAL-HYCH201104003.htm
    [9] 许军, 刘雁春, 暴景阳, 等.最小二乘水位拟合模型的参数时变规律研究[J].武汉大学学报·信息科学版, 2010, 35(12): 1418-1422 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=whchkjdxxb201012008

    Xu Jun, Liu Yanchun, Bao Jingyang, et al. Temporal Variation Law of Parameters of Least Square Fitting Model [J]. Geomatics and Information Science of Wuhan University, 2010, 35(12):1418-1422 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=whchkjdxxb201012008
    [10] 许军, 暴景阳, 章传银, 等.最小二乘水位拟合法的应用优化设置[J].海洋测绘, 2016, 36(1): 21-24 doi:  10.3969/j.issn.1671-3044.2016.01.005

    Xu Jun, Bao Jingyang, Zhang Chuanyin, et al. Optimal Application of the Least Square Fitting Method [J]. Hydrographic Surveying and Charting, 2016, 36(1):21-24 doi:  10.3969/j.issn.1671-3044.2016.01.005
    [11] 方国洪, 郑文振, 陈宗镛, 等.潮汐和潮流的分析和预报[M].北京:海洋出版社, 1986

    Fang Guohong, Zheng Wenzhen, Chen Zongyong, et al. Analysis and Prediction of Tide and Tidal Current [M]. Beijing: Ocean Press, 1986
    [12] 黄祖珂, 黄磊.潮汐原理与计算[M].青岛:中国海洋大学出版社, 2005

    Huang Zuke, Huang Lei. Tidal Theory and Calculation [M]. Qingdao: Publishing House of Ocean University of China, 2005
    [13] Boon J D. Tide and Tidal Current Analysis and Applications, Storm Surges and Sea Level Trends [M]. USA: Horwood Publishing Ltd, 2004
    [14] 刘聚, 暴景阳, 许军.基于验潮站组网的深度基准面确定[J].海洋测绘, 2015, 35(6):24-28 doi:  10.3969/j.issn.1671-3044.2015.06.006

    Liu Ju, Bao Jingyang, Xu Jun. Determination of Chart Datum Based on the Network of Tidal Station [J]. Hydrographic Surveying and Charting, 2015, 35(6):24-28 doi:  10.3969/j.issn.1671-3044.2015.06.006
    [15] 柯灏.海洋无缝垂直基准构建理论和方法研究[D].武汉: 武汉大学, 2012 http://cdmd.cnki.com.cn/article/cdmd-10486-1012397375.htm

    Ke Hao. Research on the Theory and Implementation Method of Marine Seamless Vertical Reference Surface [D]. Wuhan: Wuhan University, 2012 http://cdmd.cnki.com.cn/article/cdmd-10486-1012397375.htm
    [16] 刘雁春.海道测量基准面传递的数学模型[J].测绘学报, 2000, 30(4): 310-316 doi:  10.3321/j.issn:1001-1595.2000.04.006

    Liu Yanchun. Mathematical Models for Hydrographic Datum Transfer [J]. Acta Geodaetica et Cartographica Sinica, 2000, 30(4): 310-316 doi:  10.3321/j.issn:1001-1595.2000.04.006
    [17] 刘聚, 暴景阳, 许军.基于半参数模型的水位数据预处理方法研究[C].第二十八届海洋综合性学术研讨会, 沈阳, 2016 http://www.cnki.com.cn/Article/CJFDTotal-HYCH201802005.htm

    Liu Ju, Bao Jingyang, Xu Jun. Application of Semi-Parametric Model in Water Level Data Pre-Processing [C]. Hydrographic Surveying and Charting Symposium of the 28th General Meeting, Shenyang, China, 2016 http://www.cnki.com.cn/Article/CJFDTotal-HYCH201802005.htm
    [18] 暴景阳, 许军.卫星测高数据的潮汐提取与建模应用[M].北京:测绘出版社, 2013

    Bao Jingyang, Xu Jun. Tide Analysis from Altimeter Data and the Establishment and Application of Tide Model[M]. Beijing: Surveying and Mapping Press, 2013
    [19] 孙湘平.中国近海区域海洋[M].北京:海洋出版社, 2006

    Sun Xiangping. China's Offshore Waters [M]. Beijing: Ocean Press, 2006
    [20] 薛鸿超, 谢金赞.中国海岸带水文[M].北京:海洋出版社, 1995

    Xue Hongchao, Xie Jinzan. Chinese Coastal Zone Hydrology [M]. Beijing: Ocean Press, 1995
    [21] 武汉大学测绘学院测量平差科学组.误差理论与测量平差基础[M].武汉:武汉大学出版社, 2009

    Discipline Group of Surveying Adjustment of School of Geodesy and Geomatics of Wuhan University. Error Theory and Foundation of Surveying Adjustment [M]. Wuhan: Wuhan University Press, 2009
    [22] 海军司令部.海道测量规范: GB 12327-1988 [S].北京: 中国标准出版社, 1999

    Naval Command. Specifications for Hydrographic Survey: GB 12327-1998 [S]. Beijing: Standards Press of China, 1999
  • [1] 刘聚, 暴景阳, 许军, 周唯.  中国香港验潮站1962?2017年水位相对变化分析 . 武汉大学学报 ● 信息科学版, 2020, 45(7): 1065-1072. doi: 10.13203/j.whugis20180325
    [2] 陈欣, 翟国君, 暴景阳, 欧阳永忠, 陆秀平, 邓凯亮.  航空重力向下延拓的最小二乘配置Tikhonov正则化法 . 武汉大学学报 ● 信息科学版, 2018, 43(4): 578-585. doi: 10.13203/j.whugis20150728
    [3] 谢建, 龙四春, 李黎, 李博超.  不等式约束加权整体最小二乘的凝聚函数法 . 武汉大学学报 ● 信息科学版, 2018, 43(10): 1526-1530. doi: 10.13203/j.whugis20160507
    [4] 王乐洋, 吴飞, 吴良才.  GPS高程转换的总体最小二乘拟合推估模型 . 武汉大学学报 ● 信息科学版, 2016, 41(9): 1259-1264. doi: 10.13203/j.whugis20140421
    [5] 马洪磊, 刘成龙, 宋韬, 邹浜, 杨雪峰.  一种铁路既有线纵断面圆曲线拟合新方法 . 武汉大学学报 ● 信息科学版, 2016, 41(12): 1696-1700. doi: 10.13203/j.whugis20140612
    [6] 李晶, 李冬海, 赵拥军.  利用角度和时差的单站外辐射源定位方法 . 武汉大学学报 ● 信息科学版, 2015, 40(2): 227-232.
    [7] 姚宜斌, 黄书华, 孔建, 何军泉.  空间直线拟合的整体最小二乘算法 . 武汉大学学报 ● 信息科学版, 2014, 39(5): 571-574. doi: 10.13203/j.whugis20120104
    [8] 张清华, 隋立芬, 贾小林, 田亮, 朱永兴.  利用高精度PPS测量进行GPS-GLONASS时差监测 . 武汉大学学报 ● 信息科学版, 2014, 39(11): 1347-1351.
    [9] 陈传法, 蔡乾广.  DEM快速构建的最小二乘配置法 . 武汉大学学报 ● 信息科学版, 2013, 38(1): 86-89.
    [10] 向巍, 郭际明, 傅露.  基于垂直距离最小二乘拟合的双曲线沉降模型 . 武汉大学学报 ● 信息科学版, 2013, 38(5): 571-574.
    [11] 薛树强, 党亚民, 陈武.  最小二乘估值均方差计算的矩阵体积法 . 武汉大学学报 ● 信息科学版, 2009, 34(9): 1106-1109.
    [12] 张立华, 刘雁春, 朱庆, 韩元利.  基于动态水位的航线可行性分析方法 . 武汉大学学报 ● 信息科学版, 2008, 33(9): 892-895.
    [13] 陆秀平, 黄辰虎, 黄谟涛, 欧阳永忠.  浅水多波束测深潮汐改正技术研究 . 武汉大学学报 ● 信息科学版, 2008, 33(9): 922-925.
    [14] 李建成, 褚永海, 姜卫平, 徐新禹.  利用卫星测高资料监测长江中下游湖泊水位变化 . 武汉大学学报 ● 信息科学版, 2007, 32(2): 144-147.
    [15] 罗长林, 张正禄, 梅文胜, 邓勇.  三维直角坐标转换的一种阻尼最小二乘稳健估计法 . 武汉大学学报 ● 信息科学版, 2007, 32(8): 707-710.
    [16] 雷伟刚, 童小华, 刘大杰.  基于曲线拟合的线要素综合数据整体处理方法 . 武汉大学学报 ● 信息科学版, 2006, 31(10): 896-899.
    [17] 张勤, 陶本藻.  基于同伦法的非线性最小二乘平差统一模型 . 武汉大学学报 ● 信息科学版, 2004, 29(8): 708-710.
    [18] 李桂苓, 万剑华, 陶华学.  用差商代替导数的非线性最小二乘估计 . 武汉大学学报 ● 信息科学版, 2001, 26(2): 118-121.
    [19] 赵建虎, 刘经南, 周丰年.  GPS测定船体姿态方法研究 . 武汉大学学报 ● 信息科学版, 2000, 25(4): 353-357.
    [20] 陶本藻.  论最小二乘拟合 . 武汉大学学报 ● 信息科学版, 1980, 5(1): 27-34.
  • 加载中
图(5) / 表(3)
计量
  • 文章访问数:  703
  • HTML全文浏览量:  52
  • PDF下载量:  100
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-02-21
  • 刊出日期:  2019-05-05

利用时变潮时差进行水位改正

doi: 10.13203/j.whugis20170302
    基金项目:

    国家自然科学基金 41074002

    国家自然科学基金 41501500

    作者简介:

    刘聚, 博士生, 主要从事海道测量水位控制理论与应用方面的研究。liuju9752@163.com

    通讯作者: 暴景阳, 博士, 教授。jingyangb@sina.com
  • 中图分类号: P229

摘要: 针对目前实施时差法水位改正时使用的潮时差未考虑时变现象的影响以及利用潮时差确定水位观测时段选择混乱等问题,提出利用逐时滑动时段确定潮时差以进行水位改正的方法,分析了潮时差时变现象的原因及潮时差逐时滑动时段与逐日独立时段确定结果的差异,推导了时差法水位改正的误差公式。以某海域验潮站为例,进行了水位改正的数据实验,比较了逐时滑动时段与逐日独立时段确定的潮时差的水位改正效果。实验结果表明,利用逐时滑动确定的潮时差进行水位改正,在理论与精度方面均对目前使用的时差法有所完善与提高,尤其在潮差差异较大的海域,改进效果更加显著。

English Abstract

刘聚, 暴景阳, 许军. 利用时变潮时差进行水位改正[J]. 武汉大学学报 ● 信息科学版, 2019, 44(5): 675-681, 689. doi: 10.13203/j.whugis20170302
引用本文: 刘聚, 暴景阳, 许军. 利用时变潮时差进行水位改正[J]. 武汉大学学报 ● 信息科学版, 2019, 44(5): 675-681, 689. doi: 10.13203/j.whugis20170302
LIU Ju, BAO Jingyang, XU Jun. Water Correction Method Based on Temporal Variation of Tidal Time Difference[J]. Geomatics and Information Science of Wuhan University, 2019, 44(5): 675-681, 689. doi: 10.13203/j.whugis20170302
Citation: LIU Ju, BAO Jingyang, XU Jun. Water Correction Method Based on Temporal Variation of Tidal Time Difference[J]. Geomatics and Information Science of Wuhan University, 2019, 44(5): 675-681, 689. doi: 10.13203/j.whugis20170302
  • 水位改正是海道测量中获取稳态水深数据的重要技术环节,其目标是尽量获得待测点测深时刻的真实海面位置,并将瞬时水深值归算至深度基准面[1-2]。随着卫星定位技术的广泛应用和测深仪器的智能化发展,测深时刻平面位置的精确获取及瞬时水深的全覆盖探测已取得根本突破,而水位改正涉及对潮汐的认识、对海面的观测等若干非仪器所能处理的问题。这些技术难题需通过水位改正方法解决,而目前水位改正方法仍停留在历史阶段[3-5],因此水位改正已成为约束水深测量结果精度的主要因素,而且精密海底地形测量对水位改正技术提出了更高的精度要求[6]

    20世纪80年代,经过对水深测量中水位改正的研究,文献[7]提出了基于相关系数的时差法,实现了水位改正由传统的手工作业方式向计算机处理的转变,并且将离散的潮汐分带法改进为水位的连续分带[3]。文献[8]提出了潮汐曲线最小二乘比较法,利用潮差比、潮时差、基准面偏差3个潮汐比较参数对验潮站间的水位进行配准,为水位改正提供了一种新方法。潮汐曲线最小二乘比较法的优势在于可对水位改正结果进行精度评价。在此基础上,文献[9]研究了潮汐比较参数的时变规律的成因与影响,指出潮时差的变化周期与分潮的会合周期有关,文献[10]根据参数时变规律研究了模型应用的优化设置,对利用最小二乘曲线比较法确定潮汐比较参数时观测时长的选择提出了一些具体方案。

    目前,在实施时差法水位改正作业时确定的潮时差通常为稳定值或一天内为稳定值,传统时差法尚不能解决时变潮时差下的潮位内插问题。为降低时变潮时差对水位改正的影响,本文在现有海道测量水位改正研究的基础上,顾及潮时差的时变特性,对如何合理地划分时段确定潮时差以提高水位改正精度进行了分析和实例验证。

    • 不考虑非潮汐水位的影响,观测水位可看作若干分潮的叠加,假设共有m个分潮,则水位h(t)可表示为:

      $$h\left( t \right) = \mathop {\mathop \sum\limits_{i = 1} }\limits^m {H_i}{\rm{cos}}\left( {{\theta _i} - {g_i}} \right)$$ (1)

      式中,Hg分别表示分潮的振幅与迟角;θi=μ1iτ+μ2is+μ3ih+μ4ip+μ5iN′+μ6ip′+μ0iπ/2,θi表示各分潮的相角,其中μji为各分潮的Doodson数,τshpN′、p′为6个基本天文参数[11-13]

      在构成水位的众多分潮中,随从分潮对主分潮变化过程有调制作用,假设m个分潮中的主要分潮振幅与迟角分别为Hkgk,设其余分潮与主分潮具有差比关系κi=Hi/Hkφi=gk-gi,则式(1)可表示为:

      $$h\left( t \right) = \mathop {\mathop \sum\limits_{i = 1, i \ne k} }\limits^{m - 1} {\kappa _i}{H_k}{\rm{cos}}\left( {{\theta _k} - {g_k} + {\rm{\Delta }}{\theta _i} + {\varphi _i}} \right) + {H_k}{\rm{cos}}\left( {{\theta _k} - {g_k}} \right)$$ (2)

      式中,Δθi表示相角差,为各随从分潮与主分潮Doodson数之差与对应天文变量乘积的和。对于主分潮κk=Hk /Hk=1,φk=gk-gk=0,Δθk=0,故式(2)可以描述为:

      $$h\left( t \right) = \mathop {\mathop \sum\limits_{i = 1} }\limits^m {\kappa _i}{H_k}{\rm{cos}}\left( {{\theta _k} - {g_k} + {\rm{\Delta }}{\theta _i} + {\varphi _i}} \right)$$ (3)

      $$\left\{ {\begin{array}{*{20}{c}} {D{\rm{cos}}d = \mathop {\mathop \sum\limits_{i = 1} }\limits^m {\kappa _i}{\rm{cos}}\left( {{\rm{\Delta }}{\theta _i} + {\varphi _i}} \right) = a}\\ {D{\rm{sin}}d = \mathop {\mathop \sum\limits_{i = 1} }\limits^m {\kappa _i}{\rm{sin}}\left( {{\rm{\Delta }}{\theta _i} + {\varphi _i}} \right) = b} \end{array}} \right.$$ (4)

      则式(3)中的各分潮可合并为:

      $$h\left( t \right) = D{H_k}{\rm{cos}}\left( {{\theta _k} + d - {g_k}} \right)$$ (5)

      式中,Dd体现了随从分潮对主分潮的调制影响,分别为随从分潮对主分潮振幅的调幅因子及对迟角的调相因子,确定公式为:

      $$\left\{ {\begin{array}{*{20}{c}} {D = \sqrt {{a^2} + {b^2}} }\\ {d = {\rm{ta}}{{\rm{n}}^{ - 1}}\left( {b/a} \right)} \end{array}} \right.$$ (6)
    • 假设基准验潮站AB的水位观测分别为wA(t)、wB(t),基准站AB之间有一待定站P,水位值为wP(t),则时差法水位改正的模型[7]为:

      $${w_P}\left( t \right) = \frac{{{S_{PB}}}}{{{S_{AB}}}}{w_A}\left( {t + {\tau _{AP}}} \right) + \frac{{{S_{AP}}}}{{{S_{AB}}}}{w_B}\left( {t - {\tau _{PB}}} \right)$$ (7)

      式中,τAPτPB分别为验潮站APPB之间的潮时差; SAPSPB分别表示待定站P与验潮站AB之间的距离,本文提到的待定站与基准站之间的距离均表示待定站P在两基准站连线AB上的投影与基准站之间的距离; SAB为两验潮站间距离,且SAB =SAP + SPB

      潮位最小二乘曲线比较法的数学模型[8]为:

      $${w_P}\left( t \right){\rm{}} = {\rm{}}{\gamma _{AP}}{w_A}\left( {t{\rm{}} + {\rm{}}{\tau _{AP}}} \right){\rm{}} + {\rm{}}{\varepsilon _{AP}}$$ (8)

      式中,γAPτAPεAP分别为验潮站AP之间的潮差比、潮时差以及基准偏差。

      时差法与最小二乘曲线比较法都是对传统潮汐分带法的改进,均实现了验潮站间水位的连续时空内插,两者具有一定的等效性[1],但基本思想不同。时差法水位改正是对两个基准验潮站水位的反距离加权,自变量是待定站与两基准站间的潮时差。最小二乘曲线比较法是在考虑到相位变化的基础上,通过将基准站水位进行伸缩、平移变换至待定站水位,潮时差为基准站水位的水平移动量,并且最小二乘曲线比较法可以给出验潮站之间水位的代表误差用于水位改正的精度评定。最小二乘曲线比较法的时差参数是最小二乘原理下的最优解,受水位观测误差影响较小,相较于时差法确定潮时差的方法更稳定。因此,本文采用最小二乘比较法[3, 14-16]确定潮时差。

      确定潮时差时,水位时段划分方法包括利用一个月水位的整时段法、利用每天水位的逐日独立时段法以及以每个水位观测时刻为中间时刻的逐时滑动时段法等。一个月水位的整时段法的潮时差确定结果数量较少且为稳定量,因此本文主要对逐日独立时段与逐时滑动时段进行分析。图 1为逐日独立时段与逐时滑动时段确定潮时差的示意图。图 1wAwB为两个验潮站的水位观测曲线,τττN为独立时段法确定的潮时差,Ⅰ,Ⅱ…N表示天数,τ1τ2τn为滑动时段法确定的潮时差,1,2…n表示小时数。

      图  1  逐日独立时段与逐时滑动时段划分示意图

      Figure 1.  Daily and Hourly Time-Period-Selecting Methods

      图 1所示,逐日独立时段法是利用验潮站整天的同步水位观测数据每日计算一个潮时差,此时确定的潮时差是每日中间时刻的潮时差,即以每日中间时刻的潮时差代替同天内其他时刻的潮时差进行水位改正。逐时滑动时段法是以每个水位观测时刻为中间时刻,利用同步观测数据逐时向后滑动确定潮时差,即每个水位观测值对应一个潮时差值。通常验潮站的观测间隔为1 h,利用逐时滑动时段法可以每天算得24个潮时差值。图 1中,逐日独立时段与逐时滑动时段时差确定结果每隔24 h相等一次,即τ1= ττ25= ττ24(n-1) +1= τN,逐日独立时段法确定的潮时差是逐时滑动时段法确定潮时差的子集,逐时滑动时段法是对逐日独立时段法观测间隔内的补充,反映了逐日独立时段观测间隔内潮时差的变化。

    • 对于基准站AB及待定站P,设逐时滑动确定AB间的潮时差τAB,令r = SAP /SAB,则逐时滑动时段确定的APPB间潮时差可根据空间内插原理[2]表示为:

      $$ \left\{ {\begin{array}{*{20}{c}} {{\tau _{AP}} = r{\tau _{AB}}\begin{array}{*{20}{c}} {}&{} \end{array}}\\ {{\tau _{PB}} = \left( {1 - r} \right){\tau _{AB}}} \end{array}} \right. $$ (9)

      设逐日独立时段确定的AB间潮时差为τAB0,则两种方法引起的潮时差确定误差为ΔτAB=τAB-τAB0,由潮时差确定误差引起的水位改正误差为ΔwP(t)=wP(t)-wP0(t),将式(9)代入式(7)可得:

      $$ {w_P}\left( t \right) = \left( {1 - r} \right){\rm{}}{w_A}\left( {t + r{\tau _{AB}}} \right) + r{w_B}\left[ {t - \left( {1 - r} \right){\tau _{AB}}} \right]{\rm{}} $$ (10)

      式(10)可看作τAB的函数,设为F(τAB),对式(10)在τAB0处按泰勒级数展开,由于二次以上项很微小,可以忽略,则式(10)可表示为:

      $$ {w_P}\left( t \right) = F\left( {\tau _{AB}^0} \right) + {\left( {\frac{{\partial F}}{{\partial {\tau _{AB}}}}} \right)_0}\left( {{\tau _{AB}} - \tau _{AB}^0} \right) $$ (11)

      对式(11)整理,得:

      $$ {\rm{\Delta }}{w_P}\left( t \right) = r\left( {1 - r} \right){\rm{\Delta }}{\tau _{AB}}{\rm{\Delta }}{{w'}_{AB}} $$ (12)

      式中,w′A(t+AB0)、w′B(t+AB0-τAB0)分别为AB验潮站水位曲线的导数,由式(12)可知,两种时段划分方式下潮时差的水位改正差异与待定站P的距离系数r(1-r)、基准站水位曲线导数的差Δw′AB以及潮时差确定偏差ΔτAB有关。

    • 本文潮位数据来自某不规则半日潮海域,共有9个长期验潮站,分布如图 2所示,该海域潮波传播方向如箭头所示。选择各验潮站连续30 d的同步观测数据作为本文的实验数据,各验潮站水位观测间隔为10 min。各验潮站无缺测数据,但观测水位存在不规则抖动,为保证数据质量,使用前,本文利用半参数模型对水位数据进行了预处理[17]

      图  2  某海域各验潮站分布图

      Figure 2.  Distribution of Tidal Stations

      O1K1M2S2为实验海域振幅较大的4个分潮,验潮站AB的调和常数如表 1所示。选择M2分潮作为研究分潮间调制作用的主分潮,其余分潮视为随从分潮,根据式(4)-(6),计算一个月内随从分潮对主分潮的调幅、调相因子,分析潮时差变化与水位相互作用的关系,计算结果及对应量的变化周期如图 3所示。

      表 1  验潮站主要分潮调和常数及其他实验数据

      Table 1.  Experiment Data of Tidal Stations

      分潮 A B 单分潮时差/min
      H/cm g/(°) H/cm g/(°)
      O1 29.5 237.9 30.7 243.5 -24.1
      K1 36.8 313.9 38.6 319.9 -23.9
      M2 44.0 274.0 50.1 290.3 -33.7
      S2 17.3 329.3 19.2 347.0 -35.4

      图  3  随从分潮对主分潮调制作用

      Figure 3.  Minor Constituent's Modulation Effect on Main Constituent

      由于验潮站的K1分潮略小于M2分潮,因此图 3(a)中调幅与调相变化的过程比单纯只有一个主要分潮的情况复杂,但图 3(a)调幅与调相因子的变化仍呈现周期性。图 3(b)3(c)分别为利用谱分析法确定的调幅、调相因子变化周期。调幅、调相因子均具有约1 d、15 d的周期,因为O1K1S2M2分潮之间的会合周期分别为1.0 d、1.1 d、14.8 d。图 3(d)为实测水位数据确定的潮时差值,以验潮站A为基准站、验潮站B为待定站,确定间隔为1 h,同步时长取24 h,利用最小二乘曲线拟合模型计算两验潮站之间的潮时差。图 3(e)为实测水位确定的潮时差的谱分析结果,体现的是潮时差变化的周期。由图 3(d)3(e)可知,验潮站之间的潮时差也具有周期性变化,变化周期约为1 d、15 d,与调相因子的变化周期相符,若不存在分潮的调制作用,即当实验中仅存在M2分潮时,验潮站间的潮时差为一个稳定值,不存在周期性变化,因此可认为潮时差的变化与水位之间各分潮的相互作用有关,潮时差的变化周期由分潮间的会合周期[18]决定。另外,对AB两站间不同时长、分辨率的潮时差确定结果进行了周期项的分析,发现1 a的潮时差数据除具有与图 3(e)中相同周期项外,还具有更长的182.5 d周期项,原因是由S2K2分潮间的会合周期182.6 d所致;10 min与1 h分辨率的潮时差数据在周期项的确定上基本一致。

    • 图 4为实验海域验潮站I相对于验潮站B整时段、逐日独立时段及逐时滑动时段潮时差确定结果的比较图。如图 4(a)所示,利用整时段确定的潮时差值接近于主分潮的潮时差[9],体现不出潮时差的时变现象,显然不能用较长的数据求定潮时差用于水位改正。因此,目前通常采用逐日独立时段的水位数据确定潮时差,图 4(b)4(c)中,逐日独立时段确定的潮时差呈现阶梯形变化,而逐时滑动时段更细致地反映了潮时差的变化过程。水位观测间隔越短,滑动时段法描述的潮时差变化过程越精细,目前自动记录水位数据的压力式验潮仪水位观测间隔可达到5、6、10 min等,实例中,验潮站BI之间潮时差10 min间隔的滑动时段确定结果与逐时滑动的结果最大相差2.35 min,平均相差0.25 min,此时利用滑动时段法确定的潮时差接近连续变化。

      图  4  3种时段潮时差确定结果比较图

      Figure 4.  Comparison of Time Difference Results of Various Time-Period-Selecting Methods

      利用最小二乘曲线比较法分别计算逐日独立时段与逐时滑动时段的潮时差,将计算结果统计于表 2表 2中,“中误差”表示逐时滑动时段确定的潮时差最大的中误差值;“变幅”表示逐时滑动潮时差的变化幅度,取逐时滑动潮时差最大值与最小值的差;“最大”表示逐日与逐时滑动潮时差结果差异的最大值;“平均”表示在如图 4(c)中72~216 h、432~576 h等潮时差变化迅速的时段内逐日与逐时滑动潮时差的平均差异,取两种潮时差确定结果差的绝对值的平均值。

      表 2  逐日独立时段与逐时滑动时段潮时差差异统计/min

      Table 2.  Variations of Tidal Time Differences Determined by Daily and Hourly/min

      差异参数 基准站-待定站
      B-I I-H E-I D-I A-I F-I C-B B-A B-D C-D F-H E-F F-G
      中误差 2.54 1.67 0.95 0.93 2.34 0.86 2.32 1.90 1.58 3.98 0.79 1.03 1.32
      变幅 25.9 14.1 9.4 9.9 7.2 8.1 42.4 18.6 16.8 60.0 7.0 6.9 8.4
      最大 10.9 6.0 3.4 3.1 4.8 2.6 14.0 6.4 5.8 18.1 2.8 2.9 3.1
      平均 5.4 4.3 1.3 2.3 3.6 1.8 10.4 4.0 4.1 14.3 1.5 1.7 1.8

      表 2所示,C站至D站的潮时差存在接近60 min的变化幅度,另外,CB站的潮时差变化幅度为42 min左右,B站至I站、B站至A站、B站至D站等3组潮时差的变化幅度超过了15 min,其余各组潮时差均存在10 min左右的变化幅度,理论上顺着潮波传播方向,潮时差变化幅度较大,垂直于潮波传播方向,潮时差变化幅度较小。本文中的实验海域潮波变化复杂,在受岸线地形影响较小的开阔海域,验潮站间潮时差的变化幅度与潮波传播方向的规律还需进一步研究。表 2中,B站至I站、C站至B站、C站至D站逐日独立时段与逐时滑动时段潮时差的最大差异超过10 min,C站至D站潮时差的最大差异接近20 min,在潮时差变化剧烈的时段内,C站至B站、C站至D站潮时差差异的平均值超过了10 min,此时利用逐日独立时段计算潮时差已不合理。因此,逐日确定的潮时差只能大致反映海区验潮站之间的变化趋势,当两验潮站之间的潮时差变幅较小时尚为合理的确定值,当两验潮站之间潮时差变幅过大时,逐日潮时差的代表性较差。

    • 两基准站水位曲线及导数之差的关系如图 5所示。图 5中, P_A表示验潮站A经平移后的水位曲线wA(t+AB0),P_B表示验潮站B经平移后的水位曲线wB(t+AB0AB0),数据时长为1 d。如图 5所示,验潮站AB的水位曲线之间存在时差,而经过平移之后曲线P_AP_B间无时差,基准站水位曲线导数的差与待定站P的位置无关,只与基准站潮差差异有关。根据式(12),可得到以下结论:

      图  5  水位曲线导数差异

      Figure 5.  Derivative Variations of Water Level Curve

      1) 在Δw′AB、ΔτAB均相同的情况下,当待定站P位于基准站AB中点时,水位改正误差最大。

      在利用两个相同基准站对其连线上不同待定站进行水位改正时,对于各待定站,同一时刻的Δw′AB、ΔτAB均相同,水位改正效果仅与待定站与基准站之间的距离系数r(1-r)有关。已知P位于AB之间,则0 < r < 1,因此当r=0.5时,系数r(1-r)达到最大值,即基准站AB中点处较其他位置水位改正效果差,因为当待定站P位于中间位置时,两个基准站对水位的决定作用均较弱。

      2) 当待定站距离系数及Δw′AB相同时,潮时差确定偏差越大,水位改正误差越大。

      待定站距离系数、水位曲线导数差主要由验潮站位置及海区潮汐性质等自然属性决定,而潮时差确定偏差取决于计算潮时差时水位观测时段的划分方法。

      假设待定站位于两基准站中点,当Δw′AB=3cm/min时,ΔτAB每增加1min,水位改正误差增加约0.8cm;当Δw′AB=4cm/min时,ΔτAB每增加1min,水位改正误差增加1.0cm;当Δw′AB=5cm/min时,ΔτAB每增加1min,水位改正误差增加约1.3cm;当Δw′AB=6cm/min时,ΔτAB每增加1 min,水位改正误差增加1.5 cm。目前,中国水位改正的精度标准通常取为10 cm[6],当Δw′AB分别为3、4、5、6cm/min时,潮时差误差的限差分别为13.3、10.0、8.0、6.7 min。

      3) 当待定站距离系数及ΔτAB相同时,基准站水位曲线导数的差相差越大,水位改正误差越大。

      水位曲线的导数表示水位变化的速率,在高低潮时,水位变化较慢,导数趋于零;在涨落潮期间,导数出现极大或极小值,一般验潮站水位导数极值的大小与潮差有关。如图 5所示,在某观测时刻t处,对应相同的ΔτAB,当两基准站潮差相差越大时,导数之差越大,水位改正偏差越大。

      在本文实验海域,各验潮站之间平均潮差最大相差1.83 m,各验潮站之间水位导数最大相差约为1.86 cm/min。实验海域验潮站的潮差相对较小,但在中国一些潮汐变化较剧烈的海区,潮差可达4~7 m,甚至个别验潮站的潮差超过8 m[12, 19-20]。以杭州湾的王盘山、渔山验潮站为例,在极值情况下,两站潮差之差最大可达3.37 m,涨潮期间,王盘山的水位5 min内迅速抬升了近3 m,两站水位导数之差最大可达14.7 cm/min,此时即使逐日独立时段与逐时滑动时段确定的潮时差仅相差3 min,也可导致两者水位改正结果的偏差超过10 cm。

      本文采用图 1海域中的4组验潮站进行对时差法改进的实例计算,分别是:(1)AHI,(2)ACB,(3)DHE,(4)EGF,各组中前两个验潮站是基准站,最后一个验潮站是待定站,利用一个月的水位数据,分别使用逐日独立时段与逐时滑动时段确定的潮时差实施水位改正,比较两种时间采样间隔的水位改正精度,结果统计于表 3。其中,“时”表示逐时滑动时段法,“日”表示逐日独立时段法,“固”表示一个月整时段确定的潮时差固定值,max表示两种潮时差确定方法下水位改正值与实测值最大互差,$\hat \sigma $表示水位改正值与实测水位的拟合中误差[21],即$\hat \sigma = \sqrt {\left[ {VV} \right]/\left( {n - 1} \right)} $其中[VV]表示待定站水位改正值与实测水位残差的平方和。

      表 3  各组实验水位改正精度统计/cm

      Table 3.  Precise of Water Correction Experiments/cm

      精度指标 (1)AIH (2)ABC (3)DEH (4)EFG
      max 11.8 14.7 22.3 15.7 21.2 26.4 9.8 11.3 12.3 5.6 5.4 9.7
      $\hat \sigma $ 1.1 1.8 10.8 2.7 5.8 12.3 0.9 1.8 2.6 1.1 1.0 2.1

      表 3的4组实验中,逐时滑动潮时差的水位改正结果与实测水位的最大互差均小于逐日独立时段潮时差的水位改正与实测水位的最大互差。在第3组实验中,逐时滑动的潮时差水位改正结果与实测水位的差值均在10 cm以内,说明使用逐时滑动确定的潮时差得到的水位改正值与实测水位足够接近,利用逐时滑动潮时差进行水位改正改善了水位差异的极值情况。第4组实验中由于基准站与待定站的位置相对于其他组实验都较近,因此两种时段划分方式下的潮时差水位内插结果精度接近,但均比一个月整时段的潮时差水位改正结果的精度高。在中误差方面,4组实验中逐时滑动方法的中误差均小于逐日方法的中误差,而且逐时滑动潮时差的水位改正中误差比逐日独立时段确定的潮时差缩小近一倍。《海道测量规范》[22]中关于用多项式拟合实测水位的精度要求为水位拟合的中误差不大于3 cm,可认为表 3中逐时滑动潮时差的水位改正结果已达到验潮站水位记录成果的精度要求,逐时滑动法确定的潮时差在水位改正中起到了提高精度的作用。

      在利用时差法实施水位改正之前应对待定站位置、验潮站潮差、水位变化速率、潮时差变化幅度等情况进行具体分析。当基准站潮差相差较大或者潮时差变化剧烈或者对水深测量成果具有较高精度要求时,应该顾及潮时差的时变特性,利用逐时滑动时段法确定潮时差实施水位改正;当测深海区潮差变化较小且对测深成果精度要求足够宽松时,可采用传统的时差法实施水位改正,但为了水位改正结果的科学性,仍建议考虑潮时差的变化影响。

    • 本文利用分潮合并原理分析了潮时差时变现象与分潮会合周期的关系,比较了逐日独立时段与逐时滑动时段确定潮时差的差异,逐时滑动时段法确定的潮时差能够更细致地反映潮时差变化过程,对时差法水位改正进行了误差分析,时差法水位改正误差与待定站位置、基准站水位曲线导数差、潮时差确定偏差有关。实验分析结果表明,利用逐时滑动时段确定的潮时差实施水位改正时,理论及精度方面较传统方法均有所完善和提高。关于潮时差变化幅度与潮波传播的关系及如何对时差法水位改正进行合理的精度评估还需进一步研究。

参考文献 (22)

目录

    /

    返回文章
    返回