高頻GNSS地震監(jiān)測中頻率混疊現(xiàn)象的仿真實驗研究

劉剛 聶兆生 賈治革 張銳 喬學(xué)軍

摘要：利用振動臺實驗，通過設(shè)置多個采樣頻率監(jiān)測仿真的正弦波形及天然地震波形，研究高頻GNSS出現(xiàn)的頻率混疊現(xiàn)象及其時頻特征，并討論合理的規(guī)避混疊現(xiàn)象的采樣率設(shè)置方法。結(jié)果表明：混疊效應(yīng)在時域和頻域中均有所體現(xiàn)，混疊頻率的幅度若高于GNSS的誤差水平，其導(dǎo)致的失真將不能忽略；震級與震中距是導(dǎo)致混疊現(xiàn)象的主要因素，活動斷裂的孕震能力與地震危險性及站點與斷層面距離是決定高頻GNSS監(jiān)測站采樣頻率的重要依據(jù)。

關(guān)鍵詞：高頻GNSS；采樣率；頻率混疊；振動臺仿真實驗

中圖分類號：P3157文獻標識碼：A文章編號：1000-0666（2017）02-0186-07

0引言

逐歷元解算使GNSS能夠捕獲瞬時動態(tài)變形信號，也使得空間大地測量形變監(jiān)測從傳統(tǒng)的“零頻”（如靜態(tài)GPS、InSAR）擴展至地震學(xué)觀測窗口。自Larson等（2003）成功利用1 Hz采樣的GPS數(shù)據(jù)恢復(fù)了2002年Denali地震的遠場地震波形后，高頻GNSS數(shù)據(jù)被大量用于大地震研究中，如快速確定大震矩震級及地震預(yù)警（Blewitt et al，2009；Yang et al，2011；Wright et al，2012；Melgar et al，2012），為大震的破裂過程提供近場的強地面運動約束（Delouis et al，2010；Yue，Lay，2011；Galetzka et al，2015），開展遠場地震波面波震相的捕獲及其頻散特性研究（Davis，Smalley，2009；Hung，Rau，2013）。

高頻GNSS被譽為一種新的“地震儀”（Bilich et al，2008；方榮新等，2011），與傳統(tǒng)地震儀在監(jiān)測方式上具有顯著的不同：高頻GNSS的直接觀測對象為衛(wèi)星發(fā)射的載波信號，具體定位方法是以動態(tài)解算站-星間的位置關(guān)系獲取不同歷元的點位坐標，并最終通過歷元間坐標差分獲取包含地震動信號的點位運動軌跡，但其本質(zhì)是以地震動信號為實際觀測對象、等步長采樣的數(shù)字化傳感器，其應(yīng)用必須遵守采樣定理（Nyquists Law）。依據(jù)采樣定理，數(shù)字信號中所包含的被采集模擬信號的最高頻率小于采樣率一半時，可完整重構(gòu)被采集模擬信號，而當(dāng)被采集模擬信號的最高頻率超過儀器采樣率一半時，高出部分的信號將折疊至正常的低頻部分，造成采集信號不可逆轉(zhuǎn)的失真，即頻率混疊現(xiàn)象（萬永革，2007）。頻率混疊發(fā)生的必要條件為等步長、離散采樣的數(shù)字記錄，并且考慮到具體地震信號頻率的未知性，高頻GNSS既滿足混疊現(xiàn)象發(fā)生的條件，也存在發(fā)生的風(fēng)險。Smalley（2009）曾利用已有的強震記錄重采樣至部分低采樣率，借以說明高頻GPS若以類似低頻采樣觀測也可能產(chǎn)生混疊現(xiàn)象，但并未對真實的高頻GPS數(shù)據(jù)進行分析，因此高頻GPS存在的頻率混疊現(xiàn)象僅停留在理論分析層面，尚無實例佐證。強震加速度記錄與高頻GPS位移記錄兩者直接觀測的物理量不同，兩種儀器相互補充而無法取代彼此，且轉(zhuǎn)化的積分過程中亦存在不可忽略的誤差，利用前者作為后者的近似也有其局限性。

通過觀測衛(wèi)星信號解算出位置的瞬時坐標以得到地震信號的“間接”[HJ2mm]觀測方式?jīng)Q定了高頻GNSS必須存儲中間觀測量（載波相位），而目前的存儲、傳輸手段尚難以支持高頻GNSS以絕對高于地震信號頻率范圍的采樣率（如強震儀設(shè)計的200 Hz）進行連續(xù)觀測。目前大量的GNSS觀測網(wǎng)絡(luò)的高頻采樣率設(shè)置在1～5 Hz之間，遠低于地震儀器的采樣率，從而導(dǎo)致未能覆蓋地震信號頻率的可能性大大增加，若產(chǎn)生混疊現(xiàn)象將對以近場高頻GNSS為數(shù)據(jù)基礎(chǔ)的研究，如破裂模型反演、矩震級的快速確定等，產(chǎn)生難以忽略的影響。本文主要以高頻GNSS中采樣率設(shè)計為主要研究內(nèi)容，以振動臺仿真實驗為主要方法，分析高頻GNSS中頻率混疊現(xiàn)象的時頻特征。

1振動臺仿真實驗及數(shù)據(jù)處理

劉剛等（2014）以振動臺實驗研究高頻GNSS捕捉的震相，本文實驗設(shè)置與之相似，即將高頻GNSS天線設(shè)置于GSK-166單向線性振動臺上以驅(qū)動GNSS扼流圈天線發(fā)生位移，GNSS接收機分別以1、2、5、10、20和50 Hz的等間隔采樣率同步記錄。由于室外條件難以滿足垂直向振動，我們僅在水平軸向進行實驗（圖1）。此次實驗設(shè)計為2個部分：（1）將振動信號頻率固定為2 Hz的正弦波作為振動臺輸入；[KG-*3]（2）將1999年Izmit MW74地震震中距20 km的YPT強震儀的加速度記錄作為振動臺輸入。

我們采用動態(tài)精密單點定位策略解算高頻GNSS數(shù)據(jù)，其原理為利用衛(wèi)星軌道和鐘差對單臺接收機相位數(shù)據(jù)逐歷元進行非差定位處理，獲取測站相對于全球參考框架、以歷元間隔的高精度位移（Larson et al，2003）。具體實現(xiàn)采用武漢大學(xué)研發(fā)的高精度數(shù)據(jù)處理軟件PANDA，能夠處理采樣率高達50 Hz的GPS動態(tài)數(shù)據(jù)（Fang et al，2013）。衛(wèi)星軌道和鐘差選用IGS發(fā)布的精密產(chǎn)品，處理得到的瞬時坐標精度水平向優(yōu)于1 cm，垂直向優(yōu)于3 cm。

2實驗結(jié)果及分析

采樣定理是數(shù)字信號處理學(xué)的基本理論，即被采集信號的頻率范圍不能高于采樣率一半，否則采集信號將因頻率混疊而失真，這個頻率的限制被稱為Nyquist頻率。頻率混疊現(xiàn)象具體表現(xiàn)為高出Nyquist頻率部分的信號將偽裝成低頻部分，并折疊至低于Nyquist頻率的信號部分，從而導(dǎo)致采集數(shù)據(jù)記錄為兩部分信號能量的疊加而無法恢復(fù)出真實的信號?；殳B部分的能量造成的影響還依賴于記錄儀器所能達到的精度范圍，若高于Nyquist頻率部分的能量所造成的影響小于儀器的誤差水平，則可因記錄儀器無法分辨而忽略，若高于儀器的誤差水平，則將對信號產(chǎn)生不可忽略的影響。本文處理的高頻GNSS位移時序水平向精度為毫米級，因此我們將對儀器記錄產(chǎn)生影響的位移閾值設(shè)計為1 mm，高于Nyquist頻率部分產(chǎn)生的位移大于此閾值時，則認為混疊效應(yīng)產(chǎn)生的影響不可忽略。

圖2給出了信號卓越頻率為2 Hz的正弦波作為振動臺輸入時高頻GNSS不同采樣率下的位移時序：10 Hz及以上采樣率時序完整地恢復(fù)了正弦波形，5 Hz時序則基本上刻畫出了正弦波形的輪廓，而2 Hz及1 Hz時序近似顯示為直線。依據(jù)采樣定理，[HJ2mm]2 Hz及1 Hz時序因出現(xiàn)頻率混疊而導(dǎo)致采集信號失真。從時域上看，造成2 Hz及1 Hz采樣率失真的直接原因是相鄰采樣點間存在高于觀測精度的、可分辨的細節(jié)變化，而粗糙的采樣率不能記錄到這些變化。從頻率上看（圖3a），輸入信號卓越頻率為2 Hz時，5 Hz及以上采樣率時序覆蓋了此頻率，而2 Hz及1 Hz采樣率的Nyquist頻率分別為1 Hz及05 Hz，明顯低于2 Hz，并且高于Nyquist頻率的部分可造成位移大于1 mm并接近1 cm的幅度。高頻部分偽裝成低頻成分并折疊至低于Nyquist頻率部分的信號中，且產(chǎn)生的幅度值高頻GNSS完全可分辨時，可使得2 Hz及1 Hz采樣率的時序在低頻部分（低于2 Hz的頻帶內(nèi)）與5 Hz及以上采樣率時序有所差別（幅度值偏高），因此2 Hz及1 Hz高頻GNSS時序中包含了不可忽略的混疊效應(yīng)。為進一步研究混疊產(chǎn)生的影響，我們對50 Hz采樣時序進行了05、1、25、5、10 Hz（對應(yīng)1、2、5、10、20 Hz采樣率）以上的高通濾波（圖3b）。對05 Hz及1 Hz（對應(yīng)1、2 Hz采樣率）進行高通濾波后重構(gòu)的時序基本上呈現(xiàn)了整個正弦波形；對25、5及10 Hz（對應(yīng)5、10、20 Hz采樣率）進行高通濾波后重構(gòu)時序在GPS觀測噪聲以下，且無正弦波特征。

上述的正弦波振動臺實驗驗證了高頻GNSS出現(xiàn)頻率混疊的可能性。為進一步研究高頻GNSS在實際地震監(jiān)測中的混疊現(xiàn)象，以1999年Izmit MW74地震震中距20 km的YPT強震儀的加速度記錄為輸入數(shù)據(jù)進行振動臺實驗。考慮到振動臺對復(fù)現(xiàn)YPT加速度記錄存在難以避免的畸變，將振動臺自帶的以200 Hz采樣的拾振器記錄作為此次試驗的真值（與YPT采樣率相同），并對其進行校正（Li et al，2012），積分至位移后與高頻GNSS時序進行比較（圖4中藍色曲線）。積分后的位移時序在70 s后基本回落至零值，與實驗結(jié)尾振動平臺歸位于初始位置相對應(yīng)，表明此次試驗僅包含了地震產(chǎn)生的波動效應(yīng)，振幅為GNSS完全可分辨的厘米級。[JP2]

實驗結(jié)果表明：在時域中，10 Hz及其以上的采樣率的GNSS時序與拾振器時序擬合的很好，二者的互相關(guān)系數(shù)達90%（圖4a），完整地恢復(fù)了振動臺產(chǎn)生的位移變化。5 Hz及2 Hz采樣率時序較好地擬合了拾振器記錄，但1 Hz采樣率時序則缺失了振動中某些劇烈的變化：在10～30 s的P波震相和S波初動區(qū)間中，部分脈沖出現(xiàn)了遺漏（圖4b）。

譜分析表明相同的地震波信號，其加速度譜、速度譜及位移譜的特征不盡相同（圖5），加速度譜的高頻部分的能量豐富，卓越頻率集中在02～5 Hz區(qū)間，且在10 Hz以上的高頻部分依然存在較大的幅度；速度譜相對于加速譜，其卓越頻率有所降低，集中在02～1 Hz左右，高于1 Hz的高頻部分幅度相對較低；位移譜的卓越頻率為03 Hz，高于03 Hz頻率的幅度相對于速度譜及加速度譜為最低，表明位移時序以低頻為主，包含的高頻能量較少，因此直接以位移為觀測量的高頻GNSS與加速度計或速度計相較，其采樣率可設(shè)置為相對較低值。在位移時序上，輸入信號的卓越頻率在03 Hz左右，幅度大于1 mm閾值的頻帶為002～1 Hz（圖6a）。采樣率為1 Hz的GNSS時序無法捕捉05～1 Hz頻帶范圍的信號，而此范圍內(nèi)的信號可產(chǎn)生GNSS能分辨的頻率混疊效應(yīng)，因此采樣率1 Hz時序受到頻率混疊效應(yīng)影響而出現(xiàn)失真現(xiàn)象。大于1 Hz采樣率的其他采樣率數(shù)據(jù)雖然也受到高頻混疊效應(yīng)的影響，但混疊產(chǎn)生的幅度值低于GNSS的觀測誤差，因此產(chǎn)生的影響可忽略。此外，我們對拾振器記錄分別進行了05、1、25 Hz的高通濾波（分別對應(yīng)1、2及5 Hz采樣率）（圖6b），重構(gòu)高頻信號結(jié)果表明高于25 Hz頻率的信號產(chǎn)生的幅度值小于閾值1 mm，高于1 Hz頻率信號產(chǎn)生的幅度值為小于5 mm，而高于05 Hz頻率信號產(chǎn)生的幅度則大于1 cm，基本顯示了地震波的輪廓。

3討論

高頻GNSS出現(xiàn)頻率混疊現(xiàn)象的根本原因是由采樣率決定的Nyquist頻率低于地震信號頻率，因此調(diào)高GNSS采樣率成為避免混疊現(xiàn)象有效辦法之一。然而在哪些條件下需要調(diào)高采樣率是需要討論的問題，可從儀器接收到的地震信號位移譜頻帶特征方面考慮（相對于速度譜和加速度譜，位移譜的特征頻率更低）。Smalley（2009）提出震級是影響GNSS混疊現(xiàn)象的重要因素，6～75級地震頻譜特征包含較多可被高頻GNSS觀測到的高頻分量，這些高頻分量將導(dǎo)致距離震中較近且采樣率不高（如1 Hz）的站點產(chǎn)生混疊現(xiàn)象；此外，震中距亦是另一個重要因素，震中距越小，地下介質(zhì)濾除的高頻信息越少，儀器接收到的信號頻率相對越高，則要求采樣率設(shè)置越高。我國大陸強震頻發(fā)，6級以上地震分布具有明顯的規(guī)律性，7級以上強震基本分布于大型活動斷裂上（張培震等，2013）。因此活動斷裂的孕震能力與地震危險性及站點與斷層面距離是決定高頻GNSS監(jiān)測站采樣頻率的主要因素，而依據(jù)活動斷裂分布來布局具有不同采樣率的GNSS監(jiān)測網(wǎng)絡(luò)亦具有重要意義。

本文對Izmit地震的實驗結(jié)果支持Smalley（2009）給出的10～20 km內(nèi)站點需要采樣5 Hz的結(jié)論。同時依據(jù)蘆山地震SCTQ（震中距29 km）站點1 Hz與50 Hz GPS位移序列的異同（Lou et al，2013），認為距離斷層50 km以內(nèi)的站點的采樣率亦應(yīng)高于1 Hz。

提高采樣率將會對數(shù)據(jù)儲存和傳輸產(chǎn)生壓力。最為合理的解決方法為實時處理，用數(shù)據(jù)流的方式將原始數(shù)據(jù)存儲至數(shù)據(jù)中心以降低傳輸壓力。考慮到實時處理方法暫未成熟，且震時可能導(dǎo)致數(shù)據(jù)傳輸中斷，因此有必要增加接收機本身的存儲能力。

4結(jié)論

通過振動臺實驗，本文展示了高頻GNSS出現(xiàn)的頻率混疊現(xiàn)象并分析了其特征。與定位方式無關(guān)，高頻GNSS離散化等間隔的采樣方式滿足頻率混疊產(chǎn)生的條件，一旦信號頻率高于Nyuist頻率將導(dǎo)致混疊現(xiàn)象產(chǎn)生?；殳B效應(yīng)在時域和頻域中均有所體現(xiàn)，混疊頻率的幅度若高于GNSS的誤差水平，其導(dǎo)致的失真將不能忽略，這對震源物理過程、地震預(yù)警等基于近場數(shù)據(jù)的研究領(lǐng)域具有重大的影響。

[HTK]感謝武漢地震科學(xué)儀器研究院對仿真實驗的協(xié)助，感謝兩位匿名審稿專家的建設(shè)性修改意見和編輯的辛勤工作。

參考文獻：

方榮新，施闖，陳克杰，等2011GPS地震儀：PANDA軟件測試結(jié)果與驗證[J].武漢大學(xué)學(xué)報（信息科學(xué)版），36（4）：453-456

劉剛，聶兆生，方榮新，等2014高頻GNSS形變波的震相識別：模擬實驗與實例分析[J].地球物理學(xué)報，57（9）：2813-2825，doi：106038/cjg20140908

萬永革2007數(shù)字信號處理的MATLAB實現(xiàn)[M].北京：科學(xué)出版社

張培震，鄧起東，張竹琪，等2013中國大陸的活動斷裂、地震災(zāi)害及其動力過程[J].中國科學(xué)：地球科學(xué)，43（10）：1607-1620

BILICH A，CASSIDY J，LARSON K M2008GPS Seismology：Application to the 2002 MW79 Denali Fault Earthquake[J].Bull Seismol Soc Am，98（2）：593-606

BLEWITT G，HAMMOND W C，KREEMER C，et al2009GPS for real-time earthquake source determination and tsunami warning systems[J].Geod，83（3）：335-343

DAVIS J P，SMALLEY J R2009Love wave dispersion in central North America determined using absolute displacement seismograms from high-rate GPS[J].Geophys Res Lett，114（B11）：292-310，doi：101029/2009JB006288

DELOUIS B，NOCQUET J-M，VALLE M2010Slip distribution of the February 27，2010 MW88 Maule Earthquake，central Chile，from static and high-rate GPS，InSAR and broadband teleseismic data[J].Geophys Res Lett，37（17），L17305，doi：101029/2010GL043899

FANG R，SHI C，SONG W W，et al2013Determination of earthquake magnitude using GPS displacement waveforms from real-time precise point positioning[J].Geophys J Int，196（1）：461-472，doi：101093/gji/ggt378

GALETZKA J，MELGAR D，GENRICH J F，et al2015Slip pulse and resonance of the Kathmandu basin during the 2015 Gorkha earthquake，Nepal[J].Science，349（6252）：1091-1095

HUNG H K，RAU R J2013Surface waves of the 2011 Tohoku earthquake：Observations of Taiwans dense high-rate GPS network[J].J Geophys Res，118（1）：332-345，doi：101029/2012JB009689

LARSON KM，BODIN P，GOMBERG J2003Using 1-Hz GPS Data to Measure Deformations Caused by the Denali Fault Earthquake[J].Science，300（5624）：1421，doi：101126/science1084531

LI H，WANG D，CAI Y J，et al2012Baseline Correction of Digital Strong-Motion Data-Examples from the 2008 Wenchuan，China，Earthquake[J].Adv Mater Res，378-379，247-250

LOU Y，ZHANG W，SHI C，et al2013High-rate（1-Hz and 50-Hz）GPS Seismology：Application to the 2013 MW66 Lushan earthquake[J].Journal of Asian Earth Sciences，79：426-431，doi：http：//dxdoirog/101016/jjseases201310016

MELGAR D，BOCK Y，CROWELL W B2012Real-time centroid moment tensor determination for large earthquakes from local and regional displacement records[J].Geophys Res Lett，188（2）：703-718

SMALLEY R2009High-rate GPS：How High do We Need to go？[J].Seismol Res Lett，80（6）：1054-1061

WRIGHT T J，HOULI N，HILDYARD M，et al2012Real-time，reliable magnitudes for large earthquakes from 1Hz GPS precise point positioning：The 2011 Tohoku-Oki（Japan）earthquake[J].Geophys Res Lett，39（1），L12302，doi：101029/2012GL051894

YANG S M，NIE Z S，JIAZ G，et al2011Co-seismic displacements of 2011 Japan MW90 earthquake recorded by far-field GPS stations[J].Geodesy and Geodynamics，2（3），12-15，doi：103724/SPJ1246201100012

YUE H，LAY T2011Inversion of high-rate（1 sps）GPS data for rupture process of the 11 March 2011 Tohoku earthquake（MW91）[J].GeophysRes Lett，38（7）：752-767，doi：101029/2011GL048700[ZK）][HJ][FL）]

[STHZ][WT4HZ][JZ]Simulation Experiment on the Aliasing of High Rate GNSS[JZ]Used for Earthquake Monitoring

[WT5B1][STBZ][JZ]LIU Gang1，NIE Zhaosheng1，JIA Zhige1，ZHANG Rui2，QIAO Xuejun1

[WT5"B1X][JZ]（1 Key Laboratory of Earthquake Geodesy，Institute of Seismology，CEA，Wuhan 430071，Hubei，China）[JZ]

（2 National Earthquake Infrastructure Service，Beijing 100036，China）

[WT5HZ][JZ]Abstract[WTB1]

High rate GNSS measuring the kinematic displacements caused by earthquake with discrete and equal interval sampling rate According the Nyquist Law，if analog signals have resolve energy at the frequencies higher than the half of sampling rate，the higher frequency signals masquerading as the low frequency signal in the time and frequency domains will contaminate the sampled data This effect is called aliasing In this paper，we investigate the aliasing of high rate GNSS by using shake table tests First we used the sinusoids as input signals to confirm the aliasing effects in the high rate GNSS Then we used an acceleration record of real earthquake as the input signal to study the mechanism of aliasing effects in the high rate GNSSThe results show that the aliasing occurs both in time domain and frequency domain，and shows non-negligible effect if the masqueraded amplitude is greater than the uncertainty of GNSS positioning The earthquake magnitude and epicenter distance of GNSS site are two major factors which determine the amplitude of the aliasing effect Therefore， sampling rate of GNSS should be set according to the seismogenic ability of monitoring fault and the distances between the fault plane and GNSS sites