基于GRACE衛星數據和位錯模型的2007年印尼明古魯MW8.4地震同震重力梯度信號對比分析

李靜　萬曉云　徐銘

摘要：為評估基于GRACE重力衛星提供的時變重力場數據計算地震同震重力變化梯度的可行性，利用GRACE衛星數據計算得到2007年印尼明古魯MW8.4地震前后12個月的重力梯度變化值，使用黏彈性半空間層狀位錯模型正演計算得到同震重力梯度變化的理論值，并與GRACE重力衛星數據處理結果進行對比分析。結果表明：GRACE重力衛星數據處理結果與位錯模型正演結果量級一致，但具體數值有差異；相關性分析表明隨著距震中越近，兩種方法得到的梯度變量相關性越高，在震中區域的相關性可高達90%以上。

關鍵詞：GRACE；位錯模型；重力梯度；印尼明古魯地震

中圖分類號：P315.726文獻標識碼：A文章編號：1000-0666（2023）04-0521-08

doi：10.20015/j.cnki.ISSN1000-0666.2023.0031

0引言

地震對世界所造成的災難性影響，時刻提醒著人們要加強地震的監測和預測。然而，地面監測數據因實時性差、儀器設備容易受損等缺點無法完全滿足要求，GRACE衛星可以獲得高精度的地區時變重力場。因此，可使用衛星重力測量彌補地面地震資料的不足，提高大地震的信號檢測與監測能力。采用衛星重力數據研究地震引起的重力場變化，對于研究地震的孕育機制和預測探測都有重要意義。

基于GRACE衛星數據的研究中，大多數學者主要采用重力異常研究大地震的同震變化，如張赤軍（1999）通過重力異常對區域性擾動物質的分布進行研究；Chen等（2007）利用GRACE時變重力場模型，將2004年蘇門答臘9.3級地震前后兩年的重力變化分別作為震前震后的平均重力場變化，利用該地震的差值提取同震重力變化信號，研究表明采用傳統的差分方法可以消除背景中的平均引力，減小周期性的水文信號影響，從而提取同震重力變化。由于存在南北條帶噪聲，通常使用Wahr等（1998，2000）提出的高斯平滑濾波方法來消除GRACE衛星高頻誤差帶來的影響；鄒正波（2016）通過檢測證實了重力衛星具有探測地震前重力異常的能力；付廣裕等（2018）通過重力異常對川西地區的地殼密度結構進行研究。

事實上，重力梯度數據更能反映重力場的細部結構，對于研究地震孕震過程具有較大潛力。Wolfgang（1993）指出重力梯度以重力場的曲率描述重力場的區域結構，包含了地球物理學和大地測量學中十分重要的局部重力場信息；Wang等（2006）發現重力場變化中的信號可以通過推導重力梯度得到放大；姜永濤等（2015）利用GFZ Release 05重力衛星對2010年智利8.8級特大地震同震重力和重力梯度變化進行分析，發現利用GRACE衛星觀測數據與基于SNREI地球模型的位錯模型理論計算的同震重力變化一致；Rahimi等（2018）通過提取2006—2007年千島（Kuril Island）地震同震重力梯度變化，發現相比于觀測重力變化，只要考慮適當的信號恢復方法，重力梯度變化更適用于提取同震信號；鄭增記等（2019）利用GRACE衛星觀測數據提取到2012年蘇門答臘8.6級特大地震的同震重力梯度變化，發現其在空間形態上與斷層模型模擬結果符合較好，振幅差異較大，且重力和重力梯度變化的北向分量可以很好地抑制條帶誤差。

為探究GRACE重力衛星時變模型提取地震重力梯度變化信息的可行性，本文通過GRACE重力衛星時變模型計算得到2007年印尼明古魯8.4級地震前后12個月的重力梯度變化，再利用位錯模型進行正演計算，并將兩種方法的計算結果作對比分析，從而評估利用GRACE衛星數據探測地震同震重力梯度變化的可行性。

1地震與數據基本情況

2007年9月12日，印尼蘇門答臘東南部明古魯市附近150 km海溝水域發生MW8.4大地震，震中位置（4.517°S，101.382°E），震源深度為30 km（圖1）。9月13日又發生MW7.9余震。

郭飛霄等（2020）通過對美國德克薩斯大學空間研究中心（CSR）、德國地學研究中心（GFZ）和美國宇航局噴氣推進實驗室（JPL）三大官方機構發布的GRACE RL06版本月重力場模型數據進行比對，發現CSR發布的模型數據階方差最小。故本文選取CSR RL06版本的GRACE衛星觀測數據進行同震信號提取和分析，數據截斷至60階次，采用300 km高斯空間平滑濾波去除南北條帶噪聲影響。

本文選用美國地質調查局USGS發布的2007年印尼明古魯MW8.4地震的滑動斷層模型，由于M_W7.9余震對主震的影響不可忽略（Dai et al，2016），所以也選用M_W7.9余震的滑動斷層模型進行計算分析。在位錯模型計算過程中，根據CRUST2.0（Bassin et al，2000）數據選取主震震中地區的地球分層模型進行地震同震重力梯度變化正演計算，見表1。

2原理與方法

2.1GRACE重力衛星數據處理

2.2位錯模型計算方法

Steketee（1958）最早在地震學中引入位錯模型，并迅速建立起基于斷層研究的地球位錯理論。Okada（1992）推導得出基于半無限空間均勻介質地球模型的位錯公式。Wang等（2006）對該公式進行改進，提出考慮到層狀黏彈性介質的矩形位錯理論，獨立研發了PSGRN/PSCMP計算源程序。Sun等（2010）提出球體位錯模型理論，能夠高精度解譯200 km外的遠場同震重力變化。

本文采用Wang等（2006）提出的黏彈性半空間層狀位錯模型及其相應的配套計算程序PSGRN/PSCMP，對2007年明古魯M_W8.4大地震的同震重力變化作了正演計算。對通過上述程序計算得到重力異常后進行球諧展開，將得到的系數帶入式（2）即可計算得到重力梯度。為使位錯模型計算結果與GRACE衛星數據計算結果具有一致的空間分辨率，將位錯正演結果進行球諧截斷至GRACE衛星數據相同階次，并采用相同濾波半徑的高斯濾波進行處理。

3結果與分析

3.1GRACE計算結果

圖2給出了使用GRACE衛星數據計算得到明古魯M_W8.4地震前后12個月不同重力梯度變量的重力梯度均值變化。從圖中可以看到，在震中兩側重力梯度變量均具有非常顯著的正、負的變化。各梯度分量重力變化的極大值、極小值也均出現在震中兩側附近。從圖2d～f中，還能觀測到較明顯的高斯濾波未消除的南北異常條帶噪聲。若消除異常條帶噪聲，這將極大地減弱地球的物理信號，從而使空間分辨率下降，干擾濾波分析，因此本文后續未作進一步的濾波處理。

3.2位錯模型正演結果

為了與 GRACE衛星數據計算結果進行比較分析，本文使用黏彈性半空間層狀位錯模型進行正演計算，得到了明古魯MW8.4地震的同震重力梯度變化分布（圖3）。由圖3可知，在斷層的上盤隆起帶中，各重力梯度變量存在正向變化，觀測到的重力梯度有顯著升高的跡象，而在斷層的上盤俯沖帶存在重力負極值區。

3.3比較與分析

圖4給出了基于GRACE衛星數據和位錯模型的重力梯度變化差值分布。根據圖4可知，6個梯度變量的重力梯度差值都在同一量級且數值都較小?？傮w看，正、負重力梯度差異的空間分布范圍類似，這表明 GRACE衛星能夠探測到地震重力梯度信號，但在具體數值上，用GRACE衛星探測得到的重力梯度變化與用位錯模型的正演計算結果有一定差別。 需要說明的是，高斯濾波半徑的選擇，會影響梯度變量南北條帶噪聲的削弱情況，如果適當增大濾波半徑，可以進一步消除異常條帶，但會嚴重削弱地球真實物理信號，從而導致空間分辨率降低，不利于濾波分析。本文選用的濾波半徑為300 km，導致Vyy、Vzy、Vzz梯度變量差值存在南北條帶噪聲。

表2為使用GRACE衛星數據和位錯模型所得的各梯度變量的重力梯度最大值、最小值、標準差及差值。由表可知，各梯度變量最大值和最小值的數量級基本一致，差值在同一數量級或更小。進一步對這兩種方法的處理結果進行相關性分析，得到GRACE計算結果與位錯模型正演結果的相關系數矩陣及相關性見表3，相關性隨區域的變化如圖5所示。

根據皮爾遜設計的統計指標（謝?；ǖ?，2021；陳天異，2021）可知，相關系數0.00～±0.30是微相關，±0.30～±0.50是實相關，±0.50～±0.80是顯著相關，±0.80～±1.00是高度相關。因此，由于GRACE計算結果和位錯模型正演結果的6個梯度變量相關系數為0.38～0.51（表3），可以確定兩種方法處理地震數據的結果具有相關性，且相關關系的密切程度較高。為衡量相同區域使用這兩種方法得出的重力梯度之間的相關性，計算各梯度變量之間的皮爾遜相關系數，圖5給出了兩種方法各梯度變量相關性的區域分布，在（95°～105°E，0°～10°S）范圍內相關性較高，震中西南方向的海水域相關性效果比較好。為進一步探究距震中遠近與相關性的關系，選?。?00°～102°E，2.4°～4.4°S）為初始范圍，以1°為步長向四周擴散，對得到的13個形狀相同大小不一的矩形進行相關性分析，依次得到關于GRACE衛星數據和位錯模型正演結果的各梯度變量的相關系數，各梯度變量的相關系數隨震中距的變化情況如圖6所示。

以V_xx分量為例，在包括了震中的初始區域范圍中，兩種方法得到的各梯度變量相關性最大，為0.98。隨著區域逐漸擴大，相關系數值越來越小，最終整個區域的相關系數約為0.50。6個梯度變量相關系數變化幅度較為接近，除V_zx外，其余分量均表現出距離震中越近的區域相關性越高的特點，震中區域的相關性均在0.8之上；V_zx分量在震中區域相關性稍低，但也高于0.6。上述結果表明：GRACE衛星數據探測結果與位錯模型正演結果在震中區域具有強相關性，能夠用于地震的同震重力梯度變化研究。

4結論

本文以2007年印尼明古魯MW8.4地震為例，對 GRACE重力衛星在地震同震重力梯度變化中的探測效果進行了分析，得出以下結論：

（1）選取地震前后12個月數據為例，基于GRACE衛星數據計算得到的重力梯度變化與位錯模型正演所得同震重力梯度變化雖然在數值上存在一定差異，但在整體空間分布上較為類似；GRACE衛星與位錯模型所得數據基本保持在同一數量級，兩組數據為顯著相關或實相關。

（2）各梯度變量的區域相關系數，在震中位置附近達到最大值，隨著震中距增大，相關系數逐漸降低。

需要說明的是，雖然本文中GRACE衛星對于震中位置的重力梯度提取效果較好，而對于較大范圍的重力梯度變化探測，衛星重力產品的精度有待進一步提高。另一方面，受制于震級大小、發震類型以及發震地點的不同，GRACE衛星對于同震或震后變化的探測也會有差異，因此對于其它地震重力梯度變化以及震中位置的探測均有待進一步研究。

參考文獻：

陳天異.2021.試驗工作中統計學知識的應用［J］.市政技術，39（12）：118-123.

付廣裕，金紅林，王灼華，等.2018.汶川MW7.9地震周邊地區布格重力異常與巖石圈垂向構造應力場［J］.地震，38（2）：28-36.

郭飛霄，孫中苗，任飛龍，等.2020.GRACE RL06與RL05時變重力場模型數據初步比較分析［J］.大地測量與地球動力學，40（5）：546-550.

姜永濤，張永志，王帥，等.2015.智利MW8.8地震同震重力梯度變化［J］.地震地磁觀測與研究，36（1）：60-64.

謝?；?，馬曉飛，趙月平，等.2021.基于回歸預測的公路施工揚塵擴散預測［J］.公路，66（5）：89-92.

張赤軍.1999.用地形數據確定重力異常垂直梯度［J］.科學通報，44（6）：656-661.

鄭增記，金雙根，范麗紅.2019.利用GRACE重力和重力梯度變化估計2012年蘇門答臘地震斷層參數［J］.地球物理學報，62（11）：4129-4141.

鄒正波.2016.利用GRACE衛星重力場研究地震重力變化［D］.武漢：武漢大學.

Wolfgang T.重力測量學［M］.徐菊生譯.1993.北京：地震出版社.

Bassin C，Laske G，Masters G.2000.The current limits of resolution for surface wave tomography in North America，Eos Trans［J］.AGU，81：F897.

Chen J L，Wilson C R，Tapley B D，et al.2007.GRACE detects coseismic and postseismic deformation from the Sumatra-Andaman earthquake［J］.Geophysical Research Letters，34（13）：173-180.

Dai C，Shum C K，Guo J，et al.2016.Improved source parameter constraints for five undersea earthquakes from north component of GRACE gravity and gravity gradient change measurements［J］.Earth and Planetary Science Letters，443：118-128.

Okada Y.1992.Internal deformation due to shear and tensile faults in a half-space［J］.Bulletin of the Seismological Society of America，82（2）：1018-1040.

Rahimi A，Li J，Raoofian N M，et al.2018.On the extraction of co-seismic signal for the Kuril Island earthquakes using GRACE observations［J］.Geophysical Journal International，215（1）：346-362.

Steketee J A.1958.On volterras dislocations in a semi-infinite elastic medium［J］. Canadian Journal of Physics，36（2）：192-205.

Sun W，Okubo S，Fu G，et al.2010.General formulations of global co-seismic deformations caused by an arbitrary dislocation in a spherically symmetric earth model-applicable to deformed earth surface and space-fixed point［J］.Geophysical Journal International，177（3）：817-833.

Wahr J，Molenaar M，Bryan F.1998.Time variability of the earths gravity field：Hydrological and oceanic effects and their possible detection using GRACE［J］.Journal of Geophysical Research：Solid Earth，103（B12）：30205-30229.

Wahr J，Wingham D，Bentley C.2000.A method of combining ICESat and GRACE Satellite data to cinstrain Antarctic mass balance［J］.Journal of Geophysical Research：Solid Earth，105（B7）：16279-16294.

Wang L，Shum C K，Christopher J.2012.Gravitational gradient changes following the 2004 December 26 Sumatra-Andaman Earthquake inferred from GRACE［J］.Geophysical Journal International，191（3）：1109-1118.

Wang R，Lorenzo-Martín F，Roth F.2006.PSGRN/PSCMP—a new code for calculating co-and post-seismic deformation，geoid and gravity changes based on the viscoelastic-gravitational dislocation theory［J］.Computers & Geosciences，32（4）：527-541.

Comparative Analysis of the Coseismic Gravity Gradient Signals

of the 2007 Bengkulu，Indonesia MW8.4 Earthquake Based

on GRACE Satellite Data and the Dislocation Model

LI Jing WAN Xiaoyun XU Ming

（1.School of Land Science and Technology，China University of Geosciences（Beijing），Beijing 100083，China）

（2.National Engineering Research Center for Gas Hydrate Exploration and Development，

Guangzhou Marine Geological Survey，Guangzhou 511466，Guangdong，China）

Abstract

To assess the feasibility of the time-variable gravity field data provided by GRACE satellite in detection of the coseismic gravity gradient changes，the gravity gradient changes in 12 months before and after the 2007 Bengkulu，Indonesia MW8.4 earthquake in the epicenter and its vicinity are calculated by GRACE satellite data.The theoretical values of coseismic gravity gradient changes are calculated by the forward of the dislocation model，and then compared with the GRACE results.At the level of the order of magnitude，the GRACE results are consistent with the results from the forward of the dislocation model，though the specific values are a little different.Correlation analysis shows that the closer to the epicenter，the higher the correlation of the gradient values obtained by the two methods are，and the correlation in the epicentral area are over 90%.

Keywords：GRACE；the dislocation model；gravity gradient；the Bengkulu，Indonesia MW8.4 earthquake