王煒棟,王明孝

(1.蘭州交通大學 測繪與地理信息學院,甘肅 蘭州 730070; 2.甘肅省地理國情監(jiān)測工程實驗室,甘肅 蘭州 730070; 3.陸軍研究院科技創(chuàng)新研究中心，江蘇 無錫 214100)

0 引 言

在經濟與科學技術快速發(fā)展的今天,我國的超高層建筑呈現蓬勃發(fā)展的趨勢．超高層建筑相比于其他高層建筑可以節(jié)省空間,但是相比于其他建筑具有更明顯的動態(tài)變形,受風荷載、日照作用、溫度、地震等因素的影響更加明顯．隨著測繪新技術的發(fā)展,全球衛(wèi)星導航系統(tǒng)(GNSS) 技術可以精確地獲取監(jiān)測點的三維坐標,并且能對建筑物實現連續(xù)變形監(jiān)測．它高精度、高效、全天候和不受天氣影響的優(yōu)點更適合工程實際的應用．相比于實時動態(tài)載波相位差分(RTK)技術,精密單點定位(PPP)技術只需要一臺接收機,不需要參考站,只需要根據IGS提供的精密星歷以及鐘差等產品,即可獲取監(jiān)測點的高精度三維坐標更具有不受距離限制,操作簡單,成本較低的優(yōu)點．我國北斗衛(wèi)星導航系統(tǒng)(BDS)能提供高精度導航與定位服務,因此將BDS PPP定位技術應用于建筑物的變形監(jiān)測將是今后發(fā)展的必然趨勢．目前PPP技術主要被應用于滑坡變形監(jiān)測,畢航權[1]改進了一種PPP的數學模型,并驗證了PPP技術在滑坡變形監(jiān)測中應用的可行性;王利等[2]證明PPP技術在滑坡變形監(jiān)測中精度可以達到厘米級;李黎等[3]證明動態(tài)PPP不能滿足毫米級變形監(jiān)測,2 h以上的靜態(tài)PPP可以滿足毫米級變形監(jiān)測．本文采用BDS/GPS雙模接收機,利用PPP技術對國內某超高層建筑進行變形監(jiān)測分析為今后超高層建筑變形監(jiān)測提供一種新的監(jiān)測方法．

1 PPP技術

PPP技術只需要利用一臺GNSS接收機,根據精密星歷、高精度衛(wèi)星軌道和鐘差產品,通過參數估計,可實現高精度定位[4]．

1.1 觀測方程

在GNSS測量中,主要有兩種觀測量,偽距觀測量和載波相位觀測量．其中偽距觀測量精度較低,載波相位觀測量精度較高．二者觀測方程表示如下:

(1)

(2)

(3)

1.2 函數模型

UofC模型是PPP定位中常用的函數模型,表達式如下:

(4)

ΦIF=ρ+c(δtr-δts)+T+λIFbIF+δm+

(5)

式中的符號表示量與式(1)和(2)相同．

1.3 參數隨機模型

在PPP定位中,卡爾曼濾波參數估計模型是常用的參數估計方法,其過程如下[5]:

Xk+1=Φk+1Xk+wk,wk～N(0,Qk),

(6)

Lk+1=Hk+1Xk+1+vk+1,vk+1～N(0,Rk+1).

(7)

式中:k表示觀測歷元時刻;Xk表示狀態(tài)向量;Φk+1為轉移矩陣;wk為動態(tài)噪聲;Qk為wk的協(xié)方差陣;Lk+1為觀測向量;Hk+1為系數陣;vk+1為觀測噪聲;Rk+1為vk+1的協(xié)方差陣．

將式(6)和(7)進行線性化,進一步得到:

(8)

(9)

利用式(8)和(9)進行濾波估計,過程如下:

第一步定權:

(10)

第二部進行改正:

(11)

(12)

第三步進行預測:

(13)

(14)

式中,-表示預測,濾波實際是不斷的“預測-修正的過程”． 卡爾曼濾波在實際應用中要根據實際環(huán)境選擇合適的參數才能發(fā)揮其最佳性能．

2 數據質量分析

GNSS定位性能、精度以及可靠性很大程度上取決于GNSS觀測數據質量,因此在進行數據處理之前需要進行數據質量評估．GNSS數據的質量常規(guī)的評估指標主要有多路徑效應、信噪比(SNR)、衛(wèi)星可見數和位置精度因子(PDOP)值等幾個方面[6-7]．

本文對國內某超高層建筑進行變形監(jiān)測采樣的是BDS/GPS雙模接收機,可以同時接收到BDS和GPS觀測數據,首先對BDS與GPS數據質量進行評估．

2.1 SNR分析

SNR是信號強度噪聲的比值,是評估GNSS數據質量的一項重要指標．主要用來評估數據質量好壞,SNR越高,證明數據質量越好．

如圖1、2所示,BDS和GPS的SNR都大于30 dB-Hz,能達到GNSS數據質量評估SNR的最小值．由圖可以看出,BDS的SNR最大值可以達到49 dB-Hz,而GPS的SNR最大值可以達到52 dB-Hz．

圖1 GPS信噪比

圖2 BDS信噪比

2.2 多路徑分析

在GNSS測量中,信號在傳播過程中會受到傳播路徑上各種物質的折射,因此導致接收機收到了不同折射信號的疊加,從而導致觀測值偏離真值,觀測值與真實值之間存在誤差,這種誤差即為多路徑誤差,由于多路徑的信號傳播所引起的干涉時延效應稱為多路徑效應．

如圖3、4所示,BDS和GPS多路徑誤差都在2 m之內,但是對比BDS和GPS多路徑誤差可以發(fā)現,在BDS觀測數據中存在明顯的系統(tǒng)偏差[8-10],而在GPS觀測數據中未發(fā)現這種系統(tǒng)偏差．

圖3 GPS多路徑

圖4 BDS多路徑

2.3 衛(wèi)星數以及精度因子分析

衛(wèi)星可見數是指在監(jiān)測點接收機在監(jiān)測時間段內可以接收到的衛(wèi)星數．PDOP是衛(wèi)星分布的空間幾何強度因子,一般衛(wèi)星分布越好時,PDOP值越小,一般小于3為比較理想狀態(tài)．

如圖5、6所示,BDS和GPS衛(wèi)星數都大于6顆,而在監(jiān)測時間到12 h時BDS衛(wèi)星數增多,多于GPS衛(wèi)星數,因為監(jiān)測區(qū)域位于我國,因此BDS衛(wèi)星數多于GPS衛(wèi)星數．BDS和GPS的PDOP都要小于3,衛(wèi)星分布狀態(tài)理想．

圖5 衛(wèi)星可見數

圖6 幾何精度因子PDOP

3 實例分析

本文所進行的數據分析是基于國內某超高層建筑變形監(jiān)測工程實例,監(jiān)測時間是2017年,監(jiān)測時長為24 h,采樣隔為10 s,截止高度角為15°,接收機為華測1758接收機．而監(jiān)測點布設則是在超高層建筑樓頂四個方位固定四臺接收機,以便對超高層建筑的整體變形進行監(jiān)測分析．本文選取開始監(jiān)測之后第13-15 h共三個小時時長的數據進行分析,而監(jiān)測點則選取西側點位(D1)和北側點位(D2)兩個方向的監(jiān)測定位進行分析,只對BDS數據進行處理．

在對BDS PPP數據進行處理分析之前,先對此次PPP定位精度進行分析．由于此次采用的接收機為華測接收機,利用華測公司自帶的軟件對四個監(jiān)測點進行靜態(tài)基線解算,基線長度為26 km,基準點則采用北京房山的CORS站,進行網平差解算,得到四個點的靜態(tài)精密坐標．PPP計算則利用目前最常用的軟件RTKLIB,PPP計算坐標則選取最后一個歷元坐標作為PPP的最終結果．二者的差值則作為點位中誤差,如表1所示;統(tǒng)計兩個監(jiān)測點三個方向的RMS值如表2所示．

表1 監(jiān)測點三個方向中誤差

表2 RMS統(tǒng)計表m

監(jiān)測點ENU D10.011 70.015 50.020 7 D20.012 80.013 00.021 4

如表1、2所示,在收斂之后水平精度可以達到1 cm,豎直精度可以達到2 cm．

接下來對收斂之后連續(xù)3 h的監(jiān)測數據進行分析,主要分析監(jiān)測點坐標三個方向的形變信息,具體如圖7、8所示.

圖7 D1點變形情況

圖8 D2點變形情況

圖7、8示出了經過卡爾曼濾波后2個監(jiān)測點的位移變化情況,從圖上可以看出,超高層建筑的最大位移為6 cm,N方向和U方向的形變趨勢一致,而E方向的變形比另外兩個方向的變形小,超高層建筑變形呈周期性變化．

4 結 論

本文利用BDS PPP技術,以國內某超高層建筑實測數據為依托,先對BDS數據質量進行分析,然后分析超高層建筑復雜環(huán)境下BDS PPP定位精度及變形趨勢,得出以下結論:

1)BDS觀測數據質量良好,但是在北斗二號(BDS-2)偽距觀測數據中存在系統(tǒng)偏差;

2)在超高層建筑復雜環(huán)境下,BDS PPP水平定位能達到1 cm,垂直精度可達到2 cm,能滿足超高層建筑變形監(jiān)測的精度要求,能完整測量出超高層建筑的整體變形情況,為今后的超高層建筑變形監(jiān)測提供了一種可靠的技術手段．