張恒璟，趙 亮，王 靜，范文杰

IGS站高程非線性速度場(chǎng)建模初值確定

張恒璟1,2，趙 亮1，王 靜3，范文杰4

（1. 遼寧工程技術(shù)大學(xué)，遼寧 阜新 123000；2. 中國(guó)測(cè)繪科研究院，北京 100830；3. 北京衛(wèi)星制造廠有限公司，北京 100094；4. 北京普達(dá)迪泰科技有限公司 北京 100083）

針對(duì)國(guó)際全球衛(wèi)星導(dǎo)航系統(tǒng)服務(wù)組織（IGS）站高程非線性速度場(chǎng)建模時(shí)周期初值人為假定存在不穩(wěn)定的問題，提出以功率譜分析為手段來探測(cè)高程時(shí)間序列主要周期值的方法，并將探測(cè)得到的主要周期值作為初值建立非線性速度場(chǎng)模型：以國(guó)內(nèi)外8個(gè)IGS站20余年的高程時(shí)間序列數(shù)據(jù)為對(duì)象，利用功率譜分析的方法得到每個(gè)IGS站高程時(shí)間序列數(shù)據(jù)的主周期值；然后作為非線性速度場(chǎng)模型的周期初值進(jìn)行迭代求解；最后與非線性模型線性化人為給定初值的方法進(jìn)行對(duì)比。結(jié)果表明，當(dāng)?shù)蠼獾玫降哪Ｐ蛿M合殘差精度相等時(shí)，功率譜分析給定周期初值的方法比人為給定固定周期初值的方法在迭代次數(shù)、迭代時(shí)間上均可得到降低，且能實(shí)現(xiàn)IGS站高程時(shí)間序列數(shù)據(jù)的有效擬合。

高程時(shí)間序列；擬合模型；功率譜；模型線性化；初值

0 引言

全球衛(wèi)星導(dǎo)航系統(tǒng)（global navigation satellite system, GNSS）已經(jīng)在國(guó)民經(jīng)濟(jì)中發(fā)揮了重要作用，其中美國(guó)的全球定位系統(tǒng)（global positioning system, GPS）是世界上第一個(gè)GNSS。國(guó)際GPS服務(wù)（International GPS Service, IGS）是國(guó)際大地測(cè)量協(xié)會(huì)（The International Association of Geodesy, IAG）為支持大地測(cè)量和地球動(dòng)力學(xué)研究于1993年組建的一個(gè)國(guó)際協(xié)作組織；此后，隨著俄羅斯的全球衛(wèi)星導(dǎo)航系統(tǒng)（global navigation satellite system, GLONASS）等其他全球衛(wèi)星導(dǎo)航定位系統(tǒng)的建成及投入工作，國(guó)際GPS服務(wù)也擴(kuò)大了工作范圍，并改稱為國(guó)際GNSS服務(wù)[1]。IGS站高程數(shù)據(jù)值的變化直接或間接揭示了地球物理現(xiàn)象，如地殼的運(yùn)動(dòng)、潮汐現(xiàn)象等；因此分析和預(yù)測(cè)IGS站高程分量的變化在國(guó)防建設(shè)、社會(huì)發(fā)展中發(fā)揮著重要的作用[2]。文獻(xiàn)[3-4]對(duì)IGS站高程時(shí)間序列數(shù)據(jù)的研究發(fā)現(xiàn)，高程的運(yùn)動(dòng)表現(xiàn)出了明顯的周期性和季節(jié)性等特征且影響因素眾多。文獻(xiàn)[5-6]對(duì)全球IGS站高程時(shí)間序列數(shù)據(jù)的研究發(fā)現(xiàn)，高程具有周期性的震蕩運(yùn)動(dòng)；文獻(xiàn)[7]對(duì)IGS站高程分量進(jìn)行譜分析，發(fā)現(xiàn)高程運(yùn)動(dòng)的周期性特征；文獻(xiàn)[8]對(duì)全球IGS站高程分量的小波分辨率進(jìn)行研究，發(fā)現(xiàn)測(cè)站普遍存在季節(jié)性變化；文獻(xiàn)[9]通過小波分析，發(fā)現(xiàn)高程時(shí)間序列類似正弦函數(shù)的季節(jié)變化，并總結(jié)出了季節(jié)性變化的地球物理因素；文獻(xiàn)[10]分析了近3年川滇地區(qū)的高程時(shí)間序列數(shù)據(jù)和川滇地區(qū)相對(duì)于其他IGS站的運(yùn)動(dòng)趨勢(shì)。

本文依據(jù)多個(gè)IGS站的高程時(shí)間序列，建立具體的非線性速度場(chǎng)數(shù)學(xué)模型。利用功率譜分析獲得高程時(shí)間序列的第一主周期并將其作為非線性模型的周期初值，然后與非線性模型線性化通過人為給定周期初值的方法[11]做對(duì)比發(fā)現(xiàn)，在模型擬合殘差精度相等的情況下，本文提出的功率譜分析探測(cè)主周期的方法與非線性模型線性化人為給定周期的方法在迭代次數(shù)、迭代時(shí)間上面均可得到降低，且所建立的非線性速度場(chǎng)模型也可有效擬合高程時(shí)間序列數(shù)據(jù)。

1 速度場(chǎng)模型初值的獲取

1.1 功率譜估計(jì)

基于周期圖法功率譜估計(jì)時(shí)，為了精確獲取高程序列的第一主周期，需要剔除常數(shù)項(xiàng)和線速度項(xiàng)的影響，只保留周期項(xiàng)。通過功率譜估計(jì)法獲取的頻率-能量圖中，第一主周期所占的能量比例最大；所以只需找到最大能量所對(duì)應(yīng)的頻率即可，即對(duì)應(yīng)著高程時(shí)間序列的主周期的倒數(shù)[12]。

1.2 線性化估計(jì)

1.2.1 線性周期擬合模型

連續(xù)運(yùn)行參考站（continuously operating reference station, CORS）站高程分量的震蕩模式和周期性運(yùn)動(dòng)規(guī)律存在不同，經(jīng)大量研究發(fā)現(xiàn)，各基準(zhǔn)站的運(yùn)動(dòng)周期一年、半年項(xiàng)占主要成分，采用多項(xiàng)式函數(shù)擬合IGS站的高程時(shí)間序列為

其中：為高程時(shí)間序列的已知高程值；

1.2.2 非線性周期擬合模型

非線性速度場(chǎng)模型只是在線性模型式（3）中添加了初相，其模型為

將非線性速度場(chǎng)模型線性化，得

式中：為高程擬合值；

根據(jù)間接平差原理計(jì)算出未知參數(shù)的初值。將非線性擬合模型的初值進(jìn)行迭代以求得擬合系數(shù)。本文通過改變迭代終止條件對(duì)比模型的精度，獲得最佳迭代精度。

2 實(shí)驗(yàn)與結(jié)果分析

本文采用的實(shí)驗(yàn)數(shù)據(jù)來自IGS下設(shè)的SOPAC全球數(shù)據(jù)中心，如表1。采用中誤差計(jì)算公式（8）探測(cè)高程時(shí)間序列數(shù)據(jù)的粗差，將超過3倍中誤差的數(shù)據(jù)進(jìn)行剔除。將剔除后的數(shù)據(jù)分別通過功率譜分析獲取第一主周期，和通過線性化獲取模型初值，即

表1 IGS站高程時(shí)間序列

圖1為8個(gè)IGS站高程時(shí)間序列的頻率-能量譜圖像，橫軸為頻率，以年的倒數(shù)為單位，對(duì)應(yīng)著表2中的第一主周期的頻率；為清晰得到第一主周期，圖中只保留不低于2個(gè)月周期對(duì)應(yīng)的頻率。

表2 功率譜獲取第一主周期對(duì)應(yīng)的頻率

從圖1和表2可知，各IGS站高程時(shí)間序列的第一主周期接近1a，kunm站為368 d，bjfs站是348 d，algo站、bogo站和lhaz站大概為363 d，chat站為364d，vill站為354d，zimm站為359d，表2頻率接近非線性模型年周期頻率擬合系數(shù)。

圖1 功率譜獲取第一主周期項(xiàng)

在非線性擬合模型殘差的均方根誤差相等條件下，表3給出了2種給定初值方法所需要的迭代次數(shù)值。

表3 利用不同初值建立模型的精度對(duì)比

通過實(shí)驗(yàn)發(fā)現(xiàn)，功率譜法的迭代次數(shù)要少于線性迭代法。表明由功率譜分析獲取初值而建立的模型，迭代速度更快，提供的初值更加接近非線性模型迭代后的擬合系數(shù)。2種方法的均方根誤差相同，表明模型的精度相同，且2種方法建立的模型的均方根誤差均在毫米級(jí)，表明模型能夠精確地?cái)M合高程時(shí)間序列。

圖2為由功率譜分析法提供初值建立的非線性周期擬合模型，可以看到真實(shí)值與擬合值互差均未超過2 cm，建立的IGS站高程時(shí)間序列非線性擬合模型達(dá)到了擬合精度要求。

表4為非線性模型的擬合系數(shù)，表中的年周期項(xiàng)擬合值與功率譜獲得的第一主周期值很接近。因此可以得出功率譜分析法提供的周期初值是符合IGS站高程時(shí)間序列數(shù)據(jù)周期運(yùn)動(dòng)情況的，這很好地解決了人為給定周期初值存在的弊端。

圖2 非線性時(shí)間序列擬合模型

表4 非線性周期模型擬合系數(shù)

3 結(jié)束語(yǔ)

1）本文采用線性最小二乘和功率譜估計(jì)2種方法為CORS站高程數(shù)據(jù)非線性速度場(chǎng)模型提供初值，在模型殘差均方根誤差相同的情況下，功率譜估計(jì)提供初值的方法相較于線性最小二乘所需迭代次數(shù)少、迭代效率高，且所提供初值更加接近于非線性速度場(chǎng)模型迭代后的擬合系數(shù)。

2）功率譜分析法給定周期初值得到的結(jié)果與IGS站實(shí)際的周期值非常接近，且建立的非線性速度場(chǎng)模型能很好地?cái)M合高程時(shí)間序列，建立的模型符合IGS站高程分量的運(yùn)動(dòng)特征。

[1] 周亞, 王立峰, 張思慧, 等. IGS連續(xù)運(yùn)行參考站高程時(shí)間序列功率譜分析[J]. 太赫茲科學(xué)與電子信息學(xué)報(bào), 2014, 12(1): 104-106.

[2] 曾波, 張彥芬, 姜衛(wèi)平, 等. 山西CORS網(wǎng)基準(zhǔn)站速度場(chǎng)分析[J]. 武漢大學(xué)學(xué)報(bào)（信息科學(xué)版）, 2012, 37(12): 1401-1404.

[3] 田慧. 我國(guó)CORS站高程時(shí)間序列變化規(guī)律分析[D]. 阜新: 遼寧工程技術(shù)大學(xué), 2012: 17-29.

[4] 程鵬飛. CGCS2000板塊模型構(gòu)建[J]. 測(cè)繪學(xué)報(bào), 2013, 42(2): 159-162.

[5] 符養(yǎng). 中國(guó)大陸現(xiàn)今地殼形變與GPS 坐標(biāo)時(shí)間序列分析[D]. 上海: 中國(guó)科學(xué)院上海天文臺(tái), 2002.

[6] ZHU Guangyao, FU Yang, LI Yan. Global elevation vibration and seasonal changes derived by the analysis of GPS height[J]. Science in China, 2003(8): 765-778.

[7] 韓英, 符養(yǎng). GPS高程數(shù)據(jù)時(shí)間序列分析[J]. 武漢大學(xué)學(xué)報(bào)（信息科學(xué)版）, 2003, 28(4): 425-428.

[8] 楊強(qiáng), 黨亞民, 秘金鐘. 基于IGS連續(xù)跟蹤站的GPS高程時(shí)間序列分析[J]. 測(cè)繪科學(xué), 2007, 32(3): 55-56.

[9] 蔣志浩, 張鵬, 秘金鐘, 等. 基于CGCS2000的中國(guó)地殼水平運(yùn)動(dòng)速度場(chǎng)模型研究[J]. 測(cè)繪學(xué)報(bào), 2009, 38(6): 471-476.

[10] 楊建會(huì), 聶桂根, 薛長(zhǎng)虎. 基于GPS數(shù)據(jù)的川滇地區(qū)地殼水平運(yùn)動(dòng)研究[J]. 測(cè)繪通報(bào), 2015(6): 28-31.

[11] 張恒璟, 改健. CORS站高程非線性速度場(chǎng)估計(jì)方法[J]. 測(cè)繪科學(xué), 2017, 42(12): 146-150.

[12] 張恒璟, 程鵬飛, 郭英. 我國(guó)IGS基準(zhǔn)站高程時(shí)間序列抗差功率譜估計(jì)[J]. 測(cè)繪科學(xué), 2012, 37(3): 52-53, 65.

[13] 丁美玉, 闊永紅, 高新波. 數(shù)字信號(hào)處理-時(shí)域離散隨機(jī)信號(hào)處理[M]. 西安: 西安電子科技大學(xué)出版社, 2002: 78-95.

Initial value determination in modeling nonlinear velocity field of IGS station elevations

ZHANG Hengjing1,2, ZHAO Liang1, WANG Jing3, FAN Wenjie4

（1. Liaoning Technical University, Fuxin, Liaoning 123000, China; 2. The Chinese Academy of Surveying and Mapping, Beijing 100830, China; 3. Beijing Satellite Manufacturing Co., Ltd., Beijing 100094, China; 4. Beijing's Dady Thai Technology Co., Ltd., Beijing 100083, China）

Aiming at the problem of artificial instability in the determination of initial periodic values of nonlinear velocity field modeling for IGS station elevation, the paper proposed a method to detect the main periodic values of elevation time series by means of power spectrum analysis, and established the nonlinear velocity field model by taking the main periodic values obtained from the detection as the initial values: the elevation time series data of more than 20 years from 8 IGS stations at home and abroad were taken as the objects, and the main period values of the elevation time series data of each IGS station were obtained by using the power spectrum analysis method; then the iterative solution was conducted with the initial periodic values of the nonlinear velocity field model; finally the method of linearizing artificial initial value of nonlinear model was compared with the proposed method. Results showed that: the method of power spectrum analysis with given initial value of period could be reduced both in number of iteration and in iteration time compared with the method with fixed initial value when the accuracy of model fitting residuals obtained by iterative solution was equal; moreover, the proposed method could implement the effective fitting of time series data for the IGS station elevations.

elevation time series; fitting model; power spectrum; model linearization; initial value

P228

A

2095-4999(2019)04-0099-06

張恒璟，趙亮，王靜，等.IGS站高程非線性速度場(chǎng)建模初值確定[J].導(dǎo)航定位學(xué)報(bào),2019,7(4): 99-104.（ZHANG Hengjing, ZHAO Liang, WANG Jing, et al.Initial value determination in modeling nonlinear velocity field of IGS station elevations[J].Journal of Navigation and Positioning,2019,7(4): 99-104.）

10.16547/j.cnki.10-1096.20190418.

2018-12-17

高分對(duì)地觀測(cè)系統(tǒng)重大專項(xiàng)（GF-7衛(wèi)星高程基準(zhǔn)轉(zhuǎn)換模型構(gòu)建與應(yīng)用技術(shù)）；2017年遼寧省教育廳青年項(xiàng)目（LJ2017QL008）；中國(guó)測(cè)繪科學(xué)研究院基本科研業(yè)務(wù)費(fèi)資助（AR1918）。

張恒璟（1982—），男，河南鄧州人，中國(guó)測(cè)繪科學(xué)研究院在站博士后，副教授，研究方向?yàn)榭臻g大地測(cè)量數(shù)據(jù)處理與GNSS高程非線性運(yùn)動(dòng)。

趙亮（1993—），男，河北唐山人，碩士生，研究方向?yàn)榭臻g大地測(cè)量數(shù)據(jù)處理與GNSS高程非線性運(yùn)動(dòng)。