地殼形變GPS監測數據的統計分析研究

引文格式：酒正綱，張林廣，郁雯. 地殼形變GPS監測數據的統計分析研究[J].測繪通報，2015(2)：67-69.DOI:10.13474/j.cnki.11-2246.2015.0044

地殼形變GPS監測數據的統計分析研究

酒正綱，張林廣，郁雯

(河北建筑工程學院，河北 張家口 075024)

TheStudyofGPSCrustalDeformationMonitoringStatisticalAnalysis

JIUZhenggang，ZHANGLinguang，YUWen

摘要：地殼形變與地震孕育有非常密切的關系，地殼形變監測可為地震預報提供重要的數據依據。本文論述了地殼形變和地震孕育的相互關系，對中國地殼運動觀測網絡的建立進行說明，并采用統計軟件SPSS對區域網的GPS點位變化數據進行了統計分析，采用箱圖的形式表示數據分析的結果，得到了統計數據中一些奇異點數據，這些點位數據需要進一步研究和探討。

關鍵詞：GPS；地殼形變監測；地震預報；SPSS；統計分析

中圖分類號：P228.4文獻標識碼：B

收稿日期：2013-12-04

基金項目：華北地區地殼形變及地震預報研究(11276927)

作者簡介：酒正綱(1965—)，男，碩士，副教授，主要研究方向為GPS理論與應用。E-mail：zjkzlg001@sohu.com

一、引言

全球定位系統(GPS) 觀測可提供高精度、大范圍和準實時的地殼運動定量數據, 使得在短時間內獲取大范圍地殼運動速度場成為可能, 已廣泛地用于監測現今地殼運動[1]。本文將通過研究中國地殼運動觀測網絡各GPS觀測站點的坐標變化來找到點位變化奇異的點，這些奇異點將是地殼形變監測中重點關注的點。

前期的大量研究表明，地殼形變與地震孕育有著非常密切的關系，而地殼形變的具體數值在某種程度上可以通過GPS站點三維坐標的變化表現出來。對于地殼形變和地震孕育的相互關系，國內專家和學者論述較多，汪在森(2007)認為從地殼運動速度場的分布來看，強烈地震通常發生在地殼運動差異顯著的區域，也就是運動速度的矢量大小、方向存在顯著差異的地區[2]。王雙緒(2009)有以下論述：通過研究2001年昆侖山口西8.1級地震和汶川8.0級地震發現，通常人們認為強烈地震或巨大地震應該發生在構造運動非常強烈的區域，而這兩個強震的發生區域卻并不如此，2001年昆侖山口西8.1級強震發生在青藏高原塊體內部第二弧形帶上，該區域斷裂錯動速率不高，但兩側的昆侖山地塊和柴達木地塊運動速率南強、北弱差異十分顯著；汶川8.0級地震發生在青藏高原塊體東部邊界的龍門山斷裂中段水平區域，該區域相對運動量不大，但其兩側的羌塘地塊與華南地塊運動速率西強、東弱差異十分顯著[3]。冒愛泉(2008)利用GPS監測地殼形變得出了以下結論：地殼水平運動的方向和大小顯著變化的區域是需要密切關注的區域，強烈地震可能就發生在這樣的區域[4]。

中國地殼運動觀測網絡是1997年國家批準實施的“九五”國家重大科學工程，該工程于2000年12月全面建成并正式投入運行。中國地殼運動觀測網絡建成了由27個連續觀測基準站組成的基準網、55個定期觀測的基本站組成的基本網和1000個不定期觀測的區域站組成的區域網3部分構成的觀測網絡，以及由一個數據中心與3個數據共享子系統構成的數據傳輸與分析處理系統。本文對中國地殼運動觀測網絡GPS區域網1999年、2001年、2004年、2007年觀測資料的數據處理結果分析點位異常的點，即采用統計分析的方法，對中國地殼運動觀測網絡的地殼形變數據進行統計分析，通過統計分析尋找奇異點。

下面通過箱圖表示中國地殼運動觀測網絡區域網點坐標變化的統計特性。

二、區域網點坐標變化統計分析界面

1. 區域網點X坐標變化箱圖

從中國地殼運動觀測網絡GPS數據共享服務(http:∥www.igs.org.cn:8080/)獲取846個區域網點坐標在ITRF2000框架下X、Y和Z方向的變化數據，并輸入到SPSS系統中，編輯成文件區域網點*.sav。圖1為區域網點位變化分析界面。圖2是經過運行統計分析程序后得到的區域網點X坐標變化箱圖。

從圖2區域網點X坐標變化箱圖中可以明顯看到區域網點X坐標變化數據的最大值、最小值、中位數和兩個四分位數5個特征點，以及所有異常值的點，這些異常值點分別是H057、J019、J010、H062、H067、H056、H075、H053、H042、C044、A109、A100、A101、H218、H21、H213、C011、A103、C029、A096、A105、A146、B012、C003、C022，A137和C002。

圖1　區域網點位變化分析界面

圖2　區域網點X坐標變化箱圖

2. 區域網點Y坐標變化箱圖

圖3是經過運行統計分析程序后得到的區域網點Y坐標變化箱圖。

圖3　區域網點Y坐標變化箱圖

圖3顯示了區域網點Y坐標變化數據的最大值、最小值、中位數和兩個四分位數5個特征值，并顯示出了所有異常值的點號，這些點號分別是C002、A137、C022、、A105、C029、C007、C003、B012、A096、A103、C011、A112、C058、A138、A146、C010、A100、B021、C021、D029、C034、I027、H078、H186、H107和H137。

3. 區域網點Z坐標變化箱圖

圖4是經過運行統計分析程序后得到的區域網點Z坐標變化箱圖，圖4區域網Z坐標變化箱圖顯示了區域網點位Z坐標變化數據的最大值、最小值、中位數和兩個四分位數5個特征值，并顯示了所有異常值的點號，這些點號分別是C002、A137、C022、A105、C030、C003、B012、A103、H205、A096、C044、A135、C010、A138、A136、I009、I068、G156、I017、I075、J043、I090、I084和I088。

圖4　區域網點Z坐標變化箱圖

以上圖2—圖4分別顯示了區域網X、Y和Z坐標變化異常的點號，在以上點號中， X、Y和Z坐標變化均為異常的點號為A096、A103、A105、A137、C002、C003、C022；將這些X、Y和Z坐標變化均為異常的點確定為奇異點，這些奇異點的位置及坐標變化量見表1。

表1　奇異點的位置及坐標變化量

在數據計算之前曾認為這些奇異點會分布在全國各地，但是計算結果卻出人意外地顯示這些點幾乎全部位于河北省中南部地區(A137寧河除外)，這是非常令人關注的，希望有關研究人員對此進一步深入研究和探討。

三、結束語

本文論述了地震的孕育和地殼形變存在著密切的關系，并通過相關震例論證了部分強震震前地殼形變的特點，即強烈地震常常發生在地殼形變差異較大的區域；同時本文又論述了中國地殼運動觀測網絡的建立過程，采用了其中846個區域網點的X、Y、和Z方向變化的數據，并采用統計軟件SPSS對數據進行統計計算，對統計計算的結果進行分析，得到了一些奇異點數據，這些數據具有明顯的統計異常，本文給出了這些奇異點數據的編號和點位，以便進行深入的研究。

致謝：本文計算原始數據來自中國地殼運動觀測網絡GPS數據共享服務(http:∥www.igs.org.cn:8080/)，在此表示感謝。

參考文獻：

[1]姜衛平,周曉慧,劉經南，等.青藏高原地殼運動與應變的GPS監測研究[J]. 測繪學報，2008,37(3):285-292.

[2]汪在森,楊國華,方穎，等. 利用GPS觀測結果研究地殼運動分布動態及其與強震關系[J].國際地震動態,2007(7):32-41.

[3]王雙緒,張希,張四新，等.基于地殼形變的強地震中長期預測的針對性探討[J].國際地震動態,2009,24(7):234-236.

[4]冒愛泉, 張金華,胡伍生，等.GPS地殼形變和前兆信息與地震預測[J].現代測繪,2008,31(5):3-6.

[5]張曉亮,張希,蔣峰云，等. 汶川Ms8．0地震孕震機制及同震影響研究[J]. 國際地震動態,2010,25(8):112-116.

[6]錢進,韓志云,林學峰，等.GPS在地殼形變和板塊運動監測中的應用[J]. 測繪與空間地理信息,2013,6(36):21-24.

[7]田云鋒.GPS位置時間序列中的中長期誤差研究[J]. 國際地震動態,2012(5):46-48.

[8]劉峽.華北地區現今地殼運動及形變動力學數值模擬[D].北京：中國科學技術大學，2007.

[9]符養.中國大陸現今地殼形變與GPS坐標時間序列分析[D].上海：中國科學院上海天文臺，2002.

[10]中國地震局監測預報司組織編著. 汶川8.0級地震科學研究報告[M].北京：地震出版社,2009.