管明雷，黃辰虎，李清泉，汪馳升,4，丁 凱

(1. 深圳大學海岸帶地理環(huán)境監(jiān)測國家測繪地理信息局重點實驗室，廣東 深圳 518060； 2. 深圳大學空間信息智能感知與服務(wù)深圳市重點實驗室，廣東 深圳 518060； 3. 海軍海洋測繪研究所，天津 300061； 4. 武漢大學地理國情監(jiān)測國家測繪地理信息局重點實驗室，湖北 武漢430079)

?

利用水位修正模型精化潮位數(shù)值模擬

管明雷1,2，黃辰虎3，李清泉1,2，汪馳升1,2,4，丁 凱1,2

(1. 深圳大學海岸帶地理環(huán)境監(jiān)測國家測繪地理信息局重點實驗室，廣東 深圳 518060； 2. 深圳大學空間信息智能感知與服務(wù)深圳市重點實驗室，廣東 深圳 518060； 3. 海軍海洋測繪研究所，天津 300061； 4. 武漢大學地理國情監(jiān)測國家測繪地理信息局重點實驗室，湖北 武漢430079)

利用海洋動力模式FVCOM(finite volume coast and ocean model)模式建立了高分辨率的海州灣及其外海域潮汐數(shù)值模型，提取了整個場內(nèi)各個網(wǎng)格點上天文潮位值的時間序列。利用該區(qū)域內(nèi)分布的5個驗潮站資料提取水位修正模型，對場內(nèi)的網(wǎng)格點進行水位修正，得到高精度的格點水位時間序列，然后對修正前后的水位值進行比較。結(jié)果表明，基于水位修正后的潮位數(shù)值模擬結(jié)果水位精度提高了10 cm左右。

FVCOM；潮汐數(shù)值模擬；余水位

目前，在高精度海道測量中，驗潮站提供潮位觀測值通常受有限距離的限定，驗潮站的有效范圍取決于深度測量的精度、驗潮站與深度點上瞬時(同一時間)的最大潮高差[2]。本文利用FVCOM模式模擬計算得到模擬海域整個場的天文潮潮位值。然后根據(jù)模擬區(qū)域內(nèi)分布的長短期驗潮站資料提取水位修正模型，將修正水位利用反距離加權(quán)算法配置到模擬計算的網(wǎng)格中，得到整個網(wǎng)格每個格點的潮位時間序列，對于區(qū)域內(nèi)的任意點的潮位都可以通過插值得到高精度的潮位。

1 FVCOM水動力模型

1.1 模型方程

模型控制方程包括連續(xù)性方程、動量方程、溫鹽擴散方程和密度方程。σ坐標下的控制方程如下

(1)

(2)

(3)

(4)

(5)

ρ=ρ(T,S)

(6)

海表面和底部速度邊界條件為

(7)

(8)

式中，(τsx,τsy)和(τbx,τby)分別為海表面風切應(yīng)力和海底摩擦切應(yīng)力在x和y方向上的分量。

海表面和底部溫鹽邊界條件為

(9)

(10)

(11)

(12)

1.2 潮位數(shù)值模擬

研究區(qū)域位于南黃海，包括海州灣及其附近海域，整體范圍為119°42′—120°，35°78′—34°41′，水深數(shù)據(jù)采用電子海圖數(shù)據(jù)，調(diào)和常數(shù)選用潮汐數(shù)據(jù)反演模式OTIS。計算網(wǎng)格對岸邊和島嶼的網(wǎng)格進行了局部加密，開邊界平均分辨率約0.01°，岸線及島嶼附近平均分辨率約0.002°，共包含49 333個網(wǎng)格節(jié)點、96 252個水平網(wǎng)格單元。開邊界水位采用M2、S2、N2、K2、K1、O1、P1和Q18個分潮調(diào)和常數(shù)給出的預(yù)報水位。模式模擬了2007年1月16至2007年2月15日期間海州灣及其外海海域內(nèi)的潮汐、潮流運動狀況，計算過程中每小時輸出一次結(jié)果。如圖1所示。

圖1 區(qū)域不規(guī)則三角網(wǎng)結(jié)構(gòu)

2 水位修正模型

2.1 驗潮站資料

本文采用了日照、嵐山頭、連云港、開山島和車牛山島5個驗潮站的資料，其中日照和連云港為長期驗潮站，嵐山頭、開山島和車牛山島為短期驗潮站，短期驗潮站進行了一個月的連續(xù)驗潮。

2.2 水位修正模型設(shè)計

本文設(shè)計的水位修正模型主要包括短周期分潮的水位修正、長周期分潮的水位修正和余水位修正。

2.2.1 長周期分潮水位修正模型

長周期分潮必須利用6個月以上的潮位資料才能調(diào)和分析出長周期分潮。對日照、連云港兩個長期站一年水位資料進行調(diào)和分析，得到了Sa、SSa兩個主要的長期分潮調(diào)和常數(shù)，然后對長周期分潮調(diào)和常數(shù)進行水位預(yù)報，從而獲取兩個長期驗潮站的長周期分潮引起的水位起伏時間序列，截取2007年1月27日7時至2007年2月14日16時的水位序列作為長周期分潮水位訂正。

2.2.2 短周期分潮及余水位修正模型

本文研究區(qū)域分布2個長期驗潮站和3個短期驗潮站，截取5個驗潮站2007年1月16日0時至2007年2月15日23時的水位資料進行調(diào)和分析，得到了32個常數(shù)調(diào)和常數(shù)。在32個常數(shù)中選取了M2、S2、N2、K2、K1、O1、P1和Q18個主要分潮的調(diào)和常數(shù)進行了水位預(yù)報。修正水位算法是利用實測水位減去8個分潮的預(yù)報水位的方法得到，利用該算法提取了5個驗潮站的修正水位時間序列，其數(shù)學表達式為

R(t)=h(t)-T8(t)

(13)

式中，h(t)為實測水位；T8(t)為8個分潮的預(yù)報水位；R(t)為修正水位，修正水位包括了剩余短周期天文潮和余水位。

2.3 修正水位配置算法

通過修正水位算法得到了兩組水位修正模型，利用反距離加權(quán)算法將兩組修正模型的修正水位值插入網(wǎng)格的格點水位值中。在確定的空間尺度L和時間尺度T范圍內(nèi)，單個網(wǎng)格點的水位值的數(shù)學表達式為

(14)

(15)

3 修正模型驗證

選用了車牛山島站、連云港站2個具有代表性的驗潮站作為實測資料進行對比試驗。其中，車牛山島驗潮站為島礁短期驗潮站，距陸地43 km；連云港站為岸基長期驗潮站。提取驗潮站站點位置的模擬潮位資料，潮高基準采用理論深度基準面，將模擬潮位資料和驗潮站潮位資料的水位和調(diào)和常數(shù)分別作對比。

3.1 潮位驗證

圖2和圖3分別為截取兩個驗潮站部分模擬水位和實測水位的對比，時間為2007年1月20日0時到2007年2月5日23時?？梢郧宄吹侥M結(jié)果和實測結(jié)果吻合較好，水位誤差不超過50 cm，相位基本一致，模擬的高潮和低潮的發(fā)生時間與實測的基本一致。

圖2 車牛山站模擬水位與實測水位對比

圖3 連云港站模擬水位與實測水位對比

3.2 調(diào)和常數(shù)驗證

模型計算與驗潮站實測調(diào)和常數(shù)的差值比較見表1。

表1 修正前后調(diào)和常數(shù)誤差對比

3.3 中誤差比較

對于修正后的水位數(shù)據(jù)精度評估，采用測量數(shù)據(jù)處理理論中的中誤差指標表示，即

(16)

表2 兩個站位修正前后中誤差對比 cm

通過圖1和圖2、圖3比較可以看出，經(jīng)過水位修正后的潮位相對于模擬結(jié)果得到了明顯改進，與實測結(jié)果吻合度提高，水位誤差精度提高了10 cm左右，相位和實測結(jié)果也無差別。從表2中的中誤差對比也可以看出，模擬結(jié)果的修正前后，潮位得到了明顯改善。潮位剩余誤差主要來自于兩方面：一方面主要來自數(shù)值模擬精度不夠；另一方面驗潮站分布和數(shù)量還比較少，使得改正水位模型的精度還不夠完善。

4 結(jié) 語

本文利用了FVCOM對海州灣及其外海進行了潮汐數(shù)值模擬，利用區(qū)域內(nèi)驗潮站資料對數(shù)值進行了模擬，得到了各個網(wǎng)格點的潮位，并進行了水位修正，使得修正后的潮位精度相對于數(shù)值模擬結(jié)果有了明顯的改善。通過對修正后的網(wǎng)格水位進行插值可以得到區(qū)域內(nèi)任一點的潮位值，為近海離岸工程和海上作業(yè)提供了精準的潮位信息，解決了臨時驗潮站布設(shè)的難題。

[1] 黃祖珂.潮汐原理與計算[M].青島：中國海洋大學出版社,2005.

[2] 劉軍.余水位的主要特性分析及其在海洋測繪中的應(yīng)用[J].測繪通報,2013(9):105-107.

[3] 李家彪.多波束勘測原理技術(shù)與方法[M].青島：海洋出版社,1999.

[4] 劉雷,李寶森,李冬,等.基于余水位修正的海洋潮位推算關(guān)鍵技術(shù)研究[J].海洋測繪,2012,32(2):11-14.

[5] 張雙成, 南陽, 李振宇,等. GNSS-MR技術(shù)用于潮位變化監(jiān)測分析[J]. 測繪學報, 2016, 45(9):1042-1049.

[6] 陳春,黃辰虎,王智明,等.海道測量中幾種水位改正方法的比較與分析[J].海洋測繪,2014,34(1):16-20.

[7] 許軍,暴景陽,劉雁春,等.基于區(qū)域潮汐場模型的水位控制可行性研究[J].海洋測繪,2011,31(4):8-12.

[8] 侯世喜,黃辰虎,陸秀平,等.基于余水位配置的海洋潮汐推算研究[J].海洋測繪,2005,25(6):29-33.

[9] 倪曉東, 鄒德華, 鄧德標. 薄板樣條函數(shù)支持下的等深線追蹤算法研究[J]. 測繪通報, 2016(1):132-135.

[10] 黃辰虎,陸秀平,申家雙,等.海道測量水位改正通用模式研究[J].海洋測繪,2011,31(4):13-16.

[11] 劉洪博,潘國富,應(yīng)元康,等.基于海洋潮汐動力模型的水位改正方法研究[J].海洋學研究,2014,32(2):35-39.

[12] 吳頔. 泰國灣潮汐潮流的數(shù)值模擬與理論模型研究[D]. 青島：國家海洋局第一海洋研究所, 2013.

[13] 桂新, 秦海波, 王勝平,等. 基于相關(guān)系數(shù)迭代法的水下地形測量時間延遲探測方法研究[J]. 測繪通報, 2015(5):57-59,62.

[14] 許軍, 暴景陽, 章傳銀,等. 最小二乘水位擬合法的應(yīng)用優(yōu)化設(shè)置[J]. 海洋測繪, 2016, 36(1):21-24.

[15] 劉雷, 董玉磊, 曲萌,等. 基于潮汐模型與余水位監(jiān)控法的實例分析[J]. 海洋測繪, 2015, 35(4):36-39.

Improvement of Tidal Level Numerical Simulation Based on Water Level Correction Model

GUAN Minglei1,2,HUANG Chenhu3,LI Qingquan1,2,WANG Chisheng1,2,4,DING Kai1,2

(1. Key Laboratory for Geo-Environmental Monitoring of Coastal Zone of the National Administration of Surveying, Mappingand GeoInformation, Shenzhen University, Shenzhen 518060, China; 2. Shenzhen Key Laboratory of Spatial Smart Sensing and Services, Shenzhen University, Shenzhen 518060, China; 3. Naval Institute of Hydrographic Surveying and Charting, Tianjin 300061, China; 4. Key Laboratory for National Geographic Census and Monitoring，National Administration of Surveying, Mapping and Geoinformation, Wuhan University, Wuhan 430079, China)

The finite volume coastal ocean model (FVCOM) has been applied to establish a high-resolution numerical model of tide in Haizhou Bay. The time series of astronomical tidal values at each grid point are extracted. The water level correction model is used to extract the water level correction data from the five tide station data in this area, and the water level of the grid site is corrected. The water level time series of the grid is obtained with high precision, then the modified water level value is interpolated. The instantaneous sea surface height is obtained at a certain time. The results show that the water level accuracy is improved by about 10 cm based on the numerical simulation of water level correction and the sea surface height model has been improved obviously.

FVCOM; tidal numerical simulation; residual water level

管明雷，黃辰虎，李清泉，等.利用水位修正模型精化潮位數(shù)值模擬[J].測繪通報，2017(6)：68-71.

10.13474/j.cnki.11-2246.2017.0192.

2016-12-13

深圳市科創(chuàng)委深港創(chuàng)新圈項目(SGLH20150206152559032)；深圳市未來產(chǎn)業(yè)“創(chuàng)新鏈+產(chǎn)業(yè)鏈”融合專項扶持計劃(201507211219247860)；地理國情監(jiān)測國家測繪地理信息局重點實驗室開放基金 (2015NGCM)

管明雷(1990—)，男，博士生，主要從事潮汐模型、海域垂直基準的理論及其應(yīng)用研究。E-mail:lyggml@126.com

李清泉

P229

A

0494-0911(2017)06-0068-04