布爾莎七參數在深水港疏浚檢測中的應用研究

王 海

(浙江省紹興市上虞區水利局，浙江 紹興 312300)

1 概述

近年來我國的經濟發展突飛猛進，港口作為物資集散及周轉的最基礎設施的重要性不言而喻，集裝箱吞吐量已成為現代港口地位的主要標志，港口航道的水深條件決定了港口自身的發展規模，當水深條件不能滿足通航能力時需要進行疏浚，疏浚后須進行水下地形檢測，以了解建成后的航道水深是否滿足設計及規范要求。

深水港碼頭通常位于近海區域，一側為陸地，海域進行疏浚施工，傳統海域水下地形測量平面采用信標機進行定位，定位精度在亞米級；水位改正通常采用潮位站改正，測量船位置潮位通常采用驗潮站所觀測潮位結合潮位模型計算，在狹長航道測量時由于由于測量水域距岸邊潮位站的距離較遠，一般通過在測區附近拋設驗潮儀，再采取分帶內插的方式獲得[1- 4]，潮位改正存在潮位站水尺零點測定誤差、潮位觀測誤差、潮位改正誤差等，測深包含聲速改正誤差、時間誤差等，環境效應包含船舶姿態引起的誤差、船舶動吃水引起的誤差。JTS 257—2008《水運工程質量檢驗標準》[5]規定港口港池及航道水深大于14m時疏浚允許淺值不得大于0.3m，小于10m水深允許淺值僅為0.1m，可見疏浚質量檢測的精度要求之高，采用傳統的信標機+測深儀測深模式難以高效了解港口疏浚效果。近年來，部分科研院所在港口碼頭近海區域采用無驗潮測深作業模式實施疏浚質量檢測，控制點的選擇及解算出高精度的布爾莎七參數是獲取理想水深成果的關鍵。

2 測量坐標系的分類

地球坐標系根據原點的不同，即參考橢球的不同分為地心坐標系和參心坐標系。經過局部定位和定向，與某一地區大地水準面最佳密合的地球橢球為參考橢球，該橢球上的坐標系為參心坐標系，適用于地球局部應用，坐標原點不位于地球質心；除滿足地心定位和雙平行定向條件外，在全球范圍內與大地體最密合的地球橢球，該橢球上的坐標系為地心坐標系，適用于全球應用，坐標原點位于位于地球質心，目前主要的地心坐標系有WGS- 84、ITRF- 2000、PZ- 90及中國2000。

地球坐標系有空間大地坐標系、空間直角坐標系、平面直角坐標系三種表達形式。空間大地坐標系采用大地緯度B、大地經度L和大地高H來描述空間位置，如圖1所示。空間直角坐標系坐標原點位于參考橢球的中心，X軸指向起始子午面與赤道交點，Z軸指向參考橢球的北極，Y軸位于赤道面上，且按右手系與X軸成90°夾角，如圖2所示。平面直角坐標系是利用投影變換將空間坐標通過某種數學變換映射到平面上，如UTM、Lambuda等投影，我國常用的是高斯克呂格投影，如圖3所示。

圖1 空間大地坐標系

圖2 空間直角坐標系

圖3 高斯平面直角坐標系

3 布爾莎七參數轉換模型

不同基準下的兩個空間直角坐標系如圖4所示。坐標變換存在三個平移參數和三個旋轉參數，顧及兩個坐標系尺度不一致，還存在一個尺度變化參數。相應的坐標變換公式為：

(1)

式中：ΔX0、ΔY0、ΔZ0—平移參數；εX、εY、εZ—旋轉參數；m—尺度變化參數。

圖4 布爾莎七參數轉換模型

由空間直角坐標變換公式可以看出，用于求解布爾莎七參數的重合點至少需要3個[6- 9]，且重合點的分布最好要覆蓋整個測區。參數的轉換精度與重合點的內符合精度密切相關，為求得準確的坐標轉換參數，應選擇測區內分布均勻、精度較高的重合點，當重合點大于3個時，采用最小二乘法求最或然值[10- 12]。

一般求解七參數前應在測區采用GPS靜態相對定位測量獲取控制點的平面坐標，采用高精度的幾何水準測量方法獲取控制點的正常高，根據靜態控制自由網平差成果中的大地坐標與當地坐標求解測區布爾莎七參數。部分測區僅收集到測區的地方坐標及高程成果，可在測區中部穩定區域架設基準站或在控制點架設基準站，利用流動站采集高精度的差分大地坐標，根據控制點的兩套坐標求解七參數；城市CORS定位技術已逐漸普及，獲取的差分大地坐標能滿足對精度不高的一般工程項目；近年來千尋定位的精度已達厘米級，筆者在幾個工程項目中應用其定位的大地坐標成果求解七參數殘差均能滿足要求。

轉換模型采用的是兩套空間直角坐標求解七參，通常GPS直接獲取的是大地坐標，已有控制成果為地方坐標及幾何水準高程成果。對大地坐標先將其轉化為源橢球空間直角坐標，對地方坐標通過高斯反算將其轉化為目標橢球大地坐標，鑒于工程中使用的是正常高，建議使用正常高代替目標橢球中的大地高[13- 15]，將目標橢球中的大地坐標及大地高轉化為空間直角坐標，這樣就形成了兩對空間直角坐標求解七參數。

大地坐標轉空間直角坐標公式：

(2)

a—地球橢球的長半軸；b—地球橢球的短半軸。

高斯坐標反算大地坐標公式：

(3)

(4)

其中：

(5)

(6)

tf=tanBf

(7)

(8)

式中，L0—中央子午線經度；Bf—底點緯度，即當x=X時的子午線弧長所對應的緯度。

4 工程案例

4.1 測區概況

某深水港碼頭總長2350m。碼頭前沿停泊水域寬度為102.4m，停泊水域東、西邊線與碼頭岸線成30°角布置，設計水深17.5m，設計底標高-17.1m，滿足萬噸級集裝箱船滿載靠泊要求；船舶回旋及連接水域總面積約3.46km2，水域寬度250～1553m，設計水深16.0m，設計底高程-15.5m，滿足萬噸級集裝箱船營運吃水14.0m全潮通航；在現有港內航道基礎上向西延伸2150m，航道寬度250m，設計底高程-16.0m。

4.2 控制點分布與選擇

本項目已有控制點均位于深水港測區北側，采用D級靜態控制網獲取網點平面坐標，高程由二等幾何水準法獲取。由于建網時間相對較早，部分控制點已被破壞，完好的控制點均不在測區內。根據測區圖形選擇6個控制點用于布爾莎七參數的解算，如圖5所示。控制點在東西方向大于測區長度，南北方向解算區域部分與測區重疊，解算區域邊緣距回旋及連接水域邊緣約1175m，距港口航道邊緣約1630m，解算區域最大寬度約1359m，測區偏離解算區域的最大距離較解算區域寬度達1倍之多。鑒于深水港特殊的地形條件，重合點的布設無法包圍測區，采用經典的布爾莎七參數模型進行轉標轉換對解算區域外坐標點外擴能力需進行驗證，這是本項目的難點。

圖5 控制點與測區相對關系圖

4.3 重合點大地坐標的采集方法

本項目深水港距海岸線約32km，難以采用城市CORS系統實施控制點大地坐標的采集，因此單基站RTK作業模式相對合理。基站可架設在測區中部任意點或已知控制點上，基站設置于測區中部主要考慮采集各控制點的精度相對均勻，架設于控制點使得可參與解算的控制點相應減少，兩種作業模式下流動站采集的重合點大地坐標精度相當，單基站RTK作業模式下流動站的精度與兩者的距離密切相關，隨著流動站與基準站的距離越來越遠，兩者的電離層及對流層等誤差各項異性，流動站實時差分成果精度逐漸衰減。將基站架設于解算區域中部近測區中部控制點SK15附近，除考慮重合點大地坐標采集的精度高及均勻外，同時還考慮了測區水深測量時水深點定位的準確性。

4.4 布爾莎七參數解算

采用6個重合點的WGS- 84大地坐標與地方坐標計算布爾莎七參數，計算成果顯示，轉換平面殘差均在5mm以內，具有較高的精度，高程殘差在18mm以內，滿足相應規范要求，轉換殘差信息見表1，兩對坐標見表2。

表1 重合點空間坐標轉換殘差信息表

注：dx為北方向平面殘差；dy為東方向轉換殘差；dh為高程轉換殘差。

4.5 布爾莎七參數外業檢核

通常在參數解算完成后要進行外業檢核，作用主要有兩點：一是在輸入重合點兩對坐標時會產生偶然誤差甚至粗差，求出的參數有系統偏差；二是對測區部分區域在解算區域外圍的項目需驗證參數的外擴能力，本項目受地形條件的限制，控制點分布在測區外圍，且呈帶狀分布，必須進行參數外擴能力的驗證。根據解算的參數在測區選擇了3個控制點，在測區外選擇東南側選擇控制點ZK11，其中ZK11距解算區域邊緣控制點SK01距離2858m，約占解算區域長度的38%，ZK11距基準站距離約6116m。

表2 重合點WGS-84坐標與當地坐標

采用布爾莎七參數現場進行采集了未參與解算的控制點的地方坐標，并將實測坐標與控制點已知坐標進行了比對。檢核成果見表3。

表3 檢核坐標與已有坐標較差成果表

注：ΔX為北方向較差；ΔY為東方向較差；Δh為高程較差。

從表3可以看出，參數解算區域內檢核坐標與已有坐標平面較差最大僅為12mm，高程較差最大僅17mm，且較差大小與流動站和基站間的距離正相關，高程具有明顯的代表性；原控制網點ZK11在距解算區域邊緣2858m的條件下其平面較差僅為15mm，表現在東西方向，高程較差22mm，表明本項目的布爾莎七參數成果準確可靠，空間轉換模型對平面及高程均具有一定的外擴能力。受地形條件的限制，參數外擴能力僅能在東西方向作了驗證，南北方向無條件進行驗證，但東西方向驗證的外擴長度已達南北方向非解算區域最大寬度的1.75倍，說明采用本項目的布爾莎七參數對現有測區的大地坐標進行轉換獲取的地方坐標及高程是準確可靠的。

4.6 海底高程點精度計算

本項目采用無驗潮RTK+測深儀作業模式獲取海底高程，測量點高程誤差包含RTK-GPS測量誤差、控制點自身誤差、聲速改正誤差、時間誤差、船舶姿態引起的誤差、船舶動吃水引起的誤差等，誤差源相對較多，因此需通過主測深線與檢查線交叉點測量成果較差評價海底高程的測量精度，共檢查263個水深點，按同精度原則計算水深點高程測量中誤差，計算公式如下：

(9)

式中，n—測點數；Δi—各點較差。

M高程=±0.07m

可見，海底水深高程點具有相對的高精度，疏浚檢測成果準確，可靠地反應了浚后深水港碼頭前沿各功能區域的水下泥面高程現狀。

5 結語

(1)布爾莎七參數解算殘差建議以2cm為控制指標，如殘差過大，應剔除殘差較大的重合點，增加新的重合點。解算完成后，應實地采集未參與解算控制點的地方坐標，檢核參數的正確性，當兩者較差偏大時，應擴大驗證范圍，判斷是參數輸入誤差或是控制點自身誤差問題，直至驗證成果可靠。

(2)深水港碼頭由于地形條件的限制，解算區域通常不能覆蓋測區，參與解算七參數的控制點應盡可能擴大范圍。解算完成后應在解算范圍外長度和寬度方向進行參數外擴能力的檢核，檢核長度和寬度以未解算測區的長度和寬度為基本要求，根據已有的工程經驗，建議參數的外擴在解算區域長度和寬度1倍范圍內。