北京土層GNSS連續觀測站沉降影響的識別和改正

吳培稚　梁　芳　胡樂銀　徐　平　邢成起　孫建寶　黃雨蕊　施玉芳　季　紅

1　北京市地震局，北京市蘇州街28號，100080 2　中國地震局地質研究所地震動力學國家重點實驗室，北京市華嚴里甲1號，100029

?

北京土層GNSS連續觀測站沉降影響的識別和改正

吳培稚1梁芳1胡樂銀1徐平1邢成起1孫建寶2黃雨蕊1施玉芳1季紅1

1北京市地震局，北京市蘇州街28號，100080 2中國地震局地質研究所地震動力學國家重點實驗室，北京市華嚴里甲1號，100029

以基巖GNSS連續觀測站為背景，考察各土層GNSS連續觀測站垂直、水平分量的變化特征。觀測結果表明，排除沉降干擾，嚴重沉降站一般需要改正，一般沉降站視研究需要決定是否改正。最后給出受沉降影響的測項的改正方法。

GNSS；區域基準；沉降干擾；改正

北京GNSS連續觀測網于2005年建成[1]，大部分測站建在土層上，部分測站還設在沉降區。它們為研究沉降提供了采樣率最高、連續性最好的觀測手段，也為GNSS監測地殼運動，研究排除沉降干擾提供了有利條件。本文以基巖站為背景，考察各土層站3個分量的變化特征。

1　背景選擇

在全球框架下，北京地區GNSS站的水平分量存在南移、東移的趨勢，基巖站垂直分量呈現小幅度上升的趨勢。仿照構造區域基準的方法，選擇沉降區及周邊的基巖站構成穩定點組。北京GNSS連續觀測網共有4個基巖站，其中房山、十三陵、薊縣3個站缺數很少，由此3個站組成穩定點組；牛口峪站作為檢驗站，用以檢驗基巖區域基準的穩定性。圖1是區域基準下4個基巖站的日值變化曲線。可以看出，牛口峪站變化與其他3個站整體相同，4個測站所有測項的年平均變化量均小于1.5 mm，低于各分量日值的平均誤差，基巖區域基準是穩定的(圖中各測站曲線的標識符見表1)。

圖1　基巖站U、N、E分量變化曲線Fig.1　Time series of U, N and E components in bedrock stations

2　垂直分量沉降影響的識別和改正

基巖區域基準下，全部12個土層站的垂直分量日值變化曲線見圖2(a)，圖中缺數原因為缺測及凍土、大樹干擾等。可以看出，這12個站的下降變化各不相同。由水準、InSAR、地下水和宏觀調查得知，北京市郊存在大面積的土層沉降[2-3]，上述下降變化是由地面沉降引起的。

綜合日值觀測誤差、年平均下降量等因素可以初步設定，年平均下降量小于5 mm的測站為不沉降站，5～10 mm的為一般沉降站，超過10 mm的為嚴重沉降站。據此，12個土層站可分成3組，不沉降站有昌平、延慶、湯河口、齋堂，它們不在沉降區內；一般沉降站有大興、西集、石景山、測繪院，它們大部分在沉降區邊緣；嚴重沉降站有牛欄山、平谷、朝陽、東三旗，其中東三旗站年平均沉降量80 mm左右，為連續GNSS觀測到的沉降最嚴重的臺站。

圖2　土層測站U分量變化曲線Fig.2　Time series of U components in soil stations

在時間序列曲線圖上，土層沉降使原來的水平趨勢變化變為緩慢下降，沉降改正、期望改正后的趨勢變化接近原來的水平狀態。圖2(d)中嚴重沉降的4個測站，經多種方法試算、多項指標比較發現，東三旗、牛欄山、朝陽3個站適于用二次多項式擬合，平谷站適于用線性模型擬合。擬合改正后的曲線見圖3，各曲線年平均變化量接近0，基本排除了沉降干擾的影響。

圖3　土層嚴重沉降組各站垂直分量改正后的變化曲線Fig.3　Corrected time series of U components in severely subsidence soil stations

3　水平分量沉降影響的識別和改正

水平分量實測結果見圖4。圖4(a)、圖4(b)為不沉降組的N、E曲線，圖4(c)、圖4(d)為一般沉降組的N、E曲線，圖4(e)、4(f)為嚴重沉降組的N、E曲線。仿照上文垂直分量設定沉降指標的方法，以年平均變化量1.5 mm為界，不沉降組和一般沉降組各站的年平均變化量均小于1.5 mm(測繪院站因人為干擾，U分量曲線未繪)，可不作改正。

嚴重沉降組中平谷站年平均變化量小于1.5 mm，可不作改正。東三旗、朝陽、牛欄山站的平均年變化量一般大于1.5 mm，變化特點也不同。南北分量中，東三旗站和朝陽站北移，牛欄山站南移，東三旗站偏移量最大；東西分量中，朝陽站東移，牛欄山站和東三旗站西移，朝陽站偏移量最大。3個測站偏移的方向和大小受該站在沉降區中位置的制約，這種制約符合沉降漏斗的假設，屬沉降影響所致。排除沉降對水平分量的干擾，改正后的趨勢變化應在時序曲線圖中接近水平。改正方法可用本測站垂直分量的改正數學模型，即二次多項式。圖5是東三旗、朝陽、牛欄山站改正后的時序曲線，它們的趨勢變化接近水平。

年平均變化量的計算方法如下：在選定時段內繪日值曲線，時間以a為單位，對所有測站各分量曲線進行直線回歸，直線的斜率即是年平均變化量。表1給出了所有測站在基巖區域基準下日值曲線的斜率值、沉降改正的數學模型及所屬組別，改正后曲線斜率都接近0，改正效果較好。

表1　測站日值曲線的斜率及沉降改正模型

圖4　土層站水平分量各組變化曲線Fig.4　Time series of horizontal components in soil stations

圖5　土層嚴重沉降組N分量、E分量改正后的變化曲線Fig.5　Corrected time series of N and E components in severely subsidence soil stations

4　結　語

在以日值為基礎的GNSS測值變化曲線上，地面沉降和地殼運動均會呈現緩慢的趨勢變化，消除沉降干擾必須慎重[4-6]。測站沉降的識別和改正可分為以下幾步：1)在排除凍土、樹木、超差、遷址、設備和人為等已知干擾后，土層測站在垂直方向上表現出下降趨勢，并積累了一定程度，可確定為疑似沉降站；2)收集疑似沉降站能反映沉降的資料，找出肯定、否定沉降的證據；3)由基巖等不沉降站組成穩定點組，構造區域基準；4)適當擴大時空范圍，在區域基準下計算疑似沉降站及周邊其他測站的3個分量值；5)繪制區域基準下垂直分量曲線，土層站年平均下降量大于設定值(如10 mm)時視為嚴重沉降站；6)考慮到沉降的局部、緩慢、多變等特點，分時段選擇合適的數學方法消除垂直分量的沉降干擾，改正后的趨勢變化大體與基巖站接近；7)沉降區內測站的水平分量也會受沉降影響，這種影響的數量、方向、大小與沉降不是一一對應的，需逐站逐分量逐時段地識別和改正；8)對于一般沉降站是否需要改正，應視研究目標而定，一般應盡可能少改正。

本文給出的沉降識別指標還需要在更大范圍、更長時段內得到證實，沉降改正算法也需隨時間、地點的改變而變化。在使用GNSS數據監測地殼運動時，考慮測站有無沉降干擾是必要的。

致謝：感謝顧國華研究員、王敏研究員的悉心幫助和指導！

[1]吳培稚，孫建寶，徐平，等.利用北京GNSS連續站監測地殼運動遇到的問題及解決辦法[J].地震，2014，34(1):144-151(Wu Peizhi, Sun Jianbao, Xu Ping, et al. Problems and Solutions Encountered in Continuous GNSS Monitoring of Crustal Movement [J]. Earthquake, 2014, 34(1):144-151)[2]楊艷，賈三滿，王海剛. 北京平原區地面沉降現狀及發展趨勢分析[J].上海地質，2010，31(4):23-28(Yang Yan, Jia Sanman, Wang Haigang, et al. The Status and Development of Land Subsidence in Beijing Plain [J]. Shanghai Geology, 2010, 31(4):23-28)

[3]梁芳，孫建寶，沈正康，等.L波段InSAR數據觀測到的北京及其周邊2007-2010年間累計地殼形變基本特征[J].地震，2013，33(4):43-54(Liang Fang, Sun Jianbao, Shen Zhengkang, et al. Accumulated Crustal Deformation and Its Characteristics in Beijing and Surrounding Regions in 2007-2010 from L-band InSAR [J].Earthquake, 2013, 31(4): 43-54)

[4]黃立人. GPS基準站坐標分量時間序列的噪聲特性分析[J].大地測量與地球動力學，2006，26(2)：31-33(Huang Liren. Noise Properties in Time Series of Coordinate Component at GPS Fiducial Stations [J].Journal of Geodesy and Geodynamics, 2006, 26(2):31-33)

[5]王敏，沈正康，董大南. 非構造形變對GPS連續站位置時間序列的影響和修正[J].地球物理學報，2005，48(5)：1 045-1 052(Wang Min, Shen Zhengkang, Dong Danan, et al. Effects of Non-Tectonic Crustal Deformation on Continuous GPS Position Time Series and Correction to Them [J]. Chinese Journal of Geophysics, 2005, 48(5):1 045-1 052)

[6]顧國華. 基準站GPS連續觀測得到的垂直位移時間序列[J].地震地質，2005，27(2)：332-340(Gu Guohua. Time Series of Vertical Displacements Obtained from Continuous GPS Measurement at Fiducial Stations in the Crustal Movement of Observation Network of China [J]. Seismology and Geology, 2005, 27(2):332-340)

Foundation support:The Old Experts Scientific Research Fund of CEA, No.201502；Natural Science Foundation of Beijing, No.8041001,8092012.

About the first author:WU Peizhi, researcher, majors in earthquake monitoring and prediction, E-mail:quakewu@163.com.

Identification and Correction of Subsidence Signals for Soil-Based GNSS Continuous Stations in Beijing Area

WUPeizhi1LIANGFang1HULeyin1XUPing1XINGChengqi1SUNJianbao2HUANGYurui1SHIYufang1JIHong1

1Earthquake Administration of Beijing Municipality, 28 Suzhou Street, Beijing 100080，China 2State Key Laboratory of Earthquake Dynamics, Institute of Geology,CEA, A1 Huayanli,Beijing 100029, China

In this paper, we take rock-based GNSS continuous stations as referents, and inspect the observed deformation variations in both horizontal and vertical components of every soil-based GNSS continuous station. Observations show that severe subsidence stations must be corrected, slight subsidence stations should be corrected depending on the research goals, and no correction is needed on non-subsidence stations for continuous GPS observations. We also present some methods for subsidence correction on GNSS station observations.

GNSS; regional base station; subsidence interference; correction

2015-07-15

吳培稚，研究員，主要從事地震監測預報研究，E-mail:quakewu@163.com。

10.14075/j.jgg.2016.09.008

1671-5942(2016)09-0789-04

P315

A

項目來源：中國地震局老專家科研基金(201502)；北京市自然科學基金(8041001,8092012)。