地震重力觀測技術的新發展方向：絕對重力與重力梯度一體化觀測＊

吳 瓊 許健生

（中國地震局地球物理研究所，北京 100081）

引言

我國大陸地質構造復雜、邊界動力多樣，形成了以活動地塊為主的構造格局。地震是地球內部熱物質運移驅動地塊運動的結果[1-2]，而地球時變重力場包含了豐富的地球內部熱物質分布與運移信息。在強震孕育過程中，震中周邊數百至數千千米范圍內，常伴有殼——幔熱物質運移驅動地塊預滑的現象[3-8]，此現象會引起持續數月至數年的 10 μGal（1 μGal=1×10?8m/s2）級重力場異常變化[9]。地表 μGal級的時變微重力場測量，可為地球內部熱物質運移、地殼形變、火山、地震等地球內部動力學過程提供直接觀測數據。因此，高精度、高時空分辨力的微重力觀測資料能精確描述強震孕育過程中活動地塊的運動特征，為地震預測和地震科學研究提供有效的觀測資料。

目前，我國已建成以流動重力測量為主的重力站網[10-12]，獲得了重力變化異常與地震震級之間的經驗關系，開展了孕震機理探索研究，相關成果已在中強地震的中長期危險區預測和地球科學研究中發揮了重要作用[13-14]。但在對標國際先進水平、自然災害防治需求和技術發展趨勢等方面，目前我國重力站網在絕對重力觀測方面還存在3個方面的明顯不足： ① 絕對重力標定與控制缺失，絕對重力控制能力弱，相對重力儀器格值因子偏離和潮汐基準不統一； ② 缺乏有效的重力場源識別能力； ③ 新型觀測技術特別是重力梯度測量技術尚未被應用于地震監測[15]。主要表現如下：

（1）觀測時間的稀疏和觀測點位的離散。目前中國地震局只有5臺FG5/A10型絕對重力儀，承擔每年150個點位的絕對重力測量任務，每臺儀器每年平均需要完成30個測點的觀測，每個絕對重力觀測點1年只測2次。而絕對重力觀測點之間的距離一般是幾百千米甚至上千千米；

（2）布設空間上的不均勻、不合理。在大地震頻發的西部地區觀測點稀疏，在大地震相對較少的東部地區觀測點相對較密；

（3）目前采用的絕對觀測儀器全部依賴進口，由于價格昂貴、維修周期長、出口限制等因素的制約，不易實現大面積、高時空分辨力的觀測需求。國內多個研究機構自行開發的絕對重力儀還停留在實驗室階段，不能與中國地震臺網其他前兆觀測儀器匹配進行長期連續觀測。

因此，基于已有自主知識產權的高精度絕對重力測量技術，開展絕對重力與重力梯度一體化測量系統研制和陣列式觀測技術研究，針對捕捉大地震前地球內部熱物質運移動態和地下流體遷徙規律這個具體的科學目標，選擇合適的野外臺站，開展示范應用。通過與臺站形變、流體和測震數據的融合分析，探索地球內部熱物質運移的時空演化物理模型，獲取地震孕育、發生過程中地球物理場變化的直接科學證據，具有十分重要的應用意義。

1 高精度絕對重力觀測可為大地震孕育和短臨預測提供直接觀測證據

由于地球的 “不可入性” ，迄今為止，人類對地球內部結構的直接觀察僅達到地球表皮十幾千米的深度。目前，世界上最深的鉆井深達12262 m，而這個超深鉆井深度和地球平均半徑（6370 km）相比，只能算是 “皮毛” 深度。因此，地震學家對地球內部結構、大地震的孕育和發生過程的認識，主要還是間接的。盡管如此，幾十年來地震學家在對地球內部結構和地震成因的研究上還是取得了長足的進步[16]。

對汶川MS8.0地震的研究表明，震前的臨震預滑活動在孕震區周邊的地塊上是普遍存在的。例如，臨震前數天至數小時，在陜西省周至地震臺（Δ=590 km）所處地塊發生了數次預滑（位移），預滑導致周至臺的深井水位發生了持續時間達1小時左右的 “階躍式” 變化（圖1a）；臨震前20多分鐘，在四川省姑咱地震臺（Δ=155.1 km）所處地塊發生的預滑導致高導層中的流體有一次突然釋放，沿巖石裂隙流失。這種質量突然遷徙導致姑咱臺附近的微重力值 “階躍式” 下降約50多μGal（圖1b）[17-18]。

圖1 汶川 MS8.0 地震前周至臺深井水位儀（a） 和姑咱臺相對重力儀記錄（b） [17-18]Fig.1 Recordings from water-level meter in Zhouzhi station （a） and relative gravimeter in Guza station （b） before the Wenchuan earthquake[17-18]

2002年9月至汶川MS8.0地震前，在距震中約50 km的四川郫縣地震臺就觀測到了絕對重力值變化呈增大的趨勢，變化率為（5.01±0.7）μGal/a。特別是2002年9月——2004年1月，顯著增大約15.8 μGal[19]。Tanaka等[20]利用FG5絕對重力儀，在震前1天和震后7天時間段內幸運地觀測到在巖手山西南8 km處MS6.1地震的絕對重力變化，他們結合震區GPS、相對重力的觀測數據，確認了絕對重力觀測值偏大6 μGal是明確的同震重力異常。Chen等[21]利用2002——2008年西藏——青海——新疆——四川地區重力場測量數據，確定了重力場異常數據與2007——2010年13個MS＞6.0的地震之間的統計相關性。Kao等[22]利用臺灣地區24個點位2004——2016年連續絕對重力測量數據，通過扣除地下水等非地球動力學因素的誤差，與GPS觀測結果對比分析確定不同點位的絕對重力測量值的變化與臺風、地震、火山活動、莫霍面沉降等因素的相關關系。Timofeev等[23]基于自主研發的GABL型絕對重力儀，通過長期的絕對重力測量，研究了不同類型地震發生前的絕對重力變化。圖2為Talaya臺站和Irkutsk臺站在2008年8月27日Kultuk地震發生前后的絕對重力測量結果，距離震中地區較近的Talaya臺站絕對重力觀測結果呈現明顯的變化趨勢，最大變化量達到了 25 μGal。

圖2 Talaya 臺站（上部曲線，1992-10——2014-10，最大變化 25 μGal） 與 Irkutsk 臺站（下部曲線，1995-10——2006-10，±2 μGal） 的絕對重力測量結果變化[23]Fig.2 Changes of the absolute gravity measurements at Talaya（Upper curve，1992-10——2014-10，up to maximum 25 μGal changes） and Irkutsk （Lower curve，1995-10——2006-10，±2 μGal changes） observatory[23]

日本神火火山（Shinmoe-dake volcano）的連續絕對重力觀測數據（FG5型絕對重力儀，每小時1組，持續觀測41天）表明（圖3），每次火山噴發前后均有對應絕對重力觀測值的線性減小和恢復過程。其中VE#1——VE#3對應明顯的火山噴發，E#4——E#5對應小規模的火山噴發，T#1——T#9對應地下巖漿的運動但未噴發[24]。

圖3 連續絕對重力測量監測神火火山的孕育和噴發[24]Fig.3 Absolute gravity at the Kirisima Volcano Observatory as a function of Universal Time during the vulcanian eruption stage of Shinmoe-dake volcano in 2011[24]

基于這些高精度、高時空分辨力的絕對重力觀測數據的研究表明，地球內部的熱物質運移可能是地塊運動的動力來源，地球內部的質量遷徙可能是地震形成的主要機制[25-28]。因此，地震成因的基本問題，歸根結底是地球內部熱物質運移在地震孕育過程中的動力學問題，獲得該過程的直接觀測數據是攻克地震成因問題的關鍵所在。而目前重力觀測可能是觀測地球內部十幾公里以下熱物質運移最有效的地球物理學方法之一[29-30]。

2 陣列式高精度重力與重力梯度聯合觀測可有效提高地下場源體重力異常信號的分辨能力

由于直接的高精度重力測量并不能有效區分地下不同埋深、水平向不同間距的場源體。針對特定目標的微重力觀測需要將目標場源體的微重力信號從眾多場源體微重力信號的復雜背景中有效的提取出來。一般的方法是對每種場源體的重力擾動進行建模并從觀測數據中扣除，而更有效的方法是采用共模差分的方式消除長周期共模擾動[31-33]，達到提取目標場源體重力異常的目的。

這種測量模式大多用在航空/衛星重力梯度測量中[34]，基于球諧函數解算的衛星重力數據，雖然具有較大的覆蓋范圍，但精度相對較低[35-37]。衛星重力可觀測到大于MS7.5的張裂型或者大于MS9.0的剪切型地震所產生的同震重力變化，但對于中小地震的監測能力較弱。地面上一般采用一套相對重力儀等間距閉合測量或兩套相對重力儀等間距同步測量的方式來完成，這種方式一般測前需要對儀器進行大量可靠的動態測試和靜態測試，并需要在盡量短的時間內完成多次復測。

如果利用高精度絕對重力觀測儀器，布設陣列式的觀測系統，對于分離不同深度場源體引起的地表微重力時變信號具有明顯的效果。圖4所示的模擬計算表明，對于單位質量場源體的重力異常，水平方向間距80 m的兩臺重力儀同步測量的敏感性和分辨能力明顯提高[32]。圖5所示，垂直方向重力梯度測量的深部場源探測能力明顯提高（圖中的單位可以等比例放大，但不影響仿真計算效果）[38]。

圖4 單位質量重力異常。單臺重力儀測量（a） 、兩臺重力儀水平間距 10 m 和 80 m （b） 測量重力場源敏感性模擬計算[32]Fig.4 Gravity change per unit mass.The simulation of source sensitivity of gravity field by a single gravimeter （a） and two gravimeters separated horizontally by 10 m and 80 m （b）[32]

圖5 單位質量重力異常。垂直向不同高度的梯度測量重力場源敏感性模擬計算[38]Fig.5 Gravity change per unit mass.The simulation of source sensitivity of gravity field by gradient measurement at different heights in vertical direction[38]

由于每套相對重力儀器的格值因子標定誤差和零漂誤差不同，基于相對重力測量儀器的陣列式測量模式將不可避免的引入測量或計算誤差，影響場源體的識別精度。隨著絕對重力測量技術的不斷成熟和完善[39]，以數天為時間間隔的持續高精度絕對重力觀測成為現實。而重力梯度測量的優點在于反應地下密度異常體具有更高的靈敏度，能夠更加直接的突出目標體的邊界[40]。因此，利用陸基臺站絕對意義[41]的重力與重力梯度一體化測量，同步獲得絕對重力加速度與垂直向重力梯度觀測數據，可有效減少相對重力測量數據由于漂移、平差等因素引入的觀測誤差，提高測點區域源自地球內部介質運移引起的地表重力場異常信息的觀測精度[42]。此外，重力與梯度觀測時間的一致性還可以更好地獲取到重力場的高階部分，約束反演過程，提高計算結果的可靠性[43]和重力場源的分辨精度。

3 絕對重力與重力梯度一體化測量技術是未來地震重力監測的新發展方向

高精度的垂直向重力梯度測量數據不僅用于直接描繪地下場源體的邊界，同時也是絕對重力測量結果向不同高度歸算的關鍵數據。在絕對重力測量之前需要完成場地重力梯度的測量，一般的測量方式是借助兩套甚至多套相對重力儀，通過垂直方向不同高度位置多次往復閉合測量的方式完成垂直向重力梯度的測量[44]。由于絕對重力測量的精度一般在5 μGal以內，常用的CG5、貝爾雷斯和拉科斯特相對重力的精度為10 μGal，目前最先進的CG6型相對重力儀的標稱精度為5 μGal，因此，利用相對重力測量得到的梯度值用來歸算絕對重力測量結果時，梯度測量的準確度偏差會直接帶入最終的絕對重力儀測量結果的誤差評估之中[39]。

直接的重力梯度測量需求從1971年美國空軍首次提出 1 E（1 E=0.1 μGal/m）移動級重力梯度測量儀器至今，重力梯度測量儀器的研究得到了世界科學家的高度重視，并取得了迅速發展。測量原理從旋轉差分加速度計、扭力測量、靜電懸浮測量、超導測量等幾個方向發展[45-46]，同時其測量方式也由靜態到動態，從船載到航空，由航空向著星載發展。目前唯一商用的高精度重力梯度儀采用旋轉差分加速度計測量原理，具有高穩定性和測量精度[47]，成功應用于航空和航海重力梯度測量。超導重力測量原理與彈簧式相對重力測量原理相比，具有更小的零漂，是具有重要發展前景的一類靜態重力梯度測量儀器[48]。而基于MEMS加速度計的重力梯度儀以其質量輕和體積小成為了未來航空重力和重力梯度測量的首選[49]。

但上述的梯度測量屬于相對測量范疇，要獲得高精度的重力梯度測量結果，一方面需要加速度測量傳感器技術指標優于梯度測量需求和不同傳感器之間測量性能指標極高的一致性，另一方面還需要對每個加速度傳感器的標度因子進行定期標定以消除漂移誤差累積和不一致帶來的梯度測量誤差。而絕對測量范疇的重力梯度測量儀器，具有無零漂、無需標定、無動態范圍限制的優點，因此，20世紀末美國宇航局的火星探測計劃中，就資助美國從事絕對重力儀生產的Micro-g公司，首次開展絕對測量范疇的激光干涉重力梯度測量的原理性試驗，目的是用于火星的深部結構探測。隨著激光干涉技術和原子干涉技術的不斷發展和完善，激光干涉重力梯度儀和原子干涉重力梯度儀將是未來絕對測量范疇的重力梯度測量儀器的主要發展方向[50]。

4 結論

基于具有自主知識產權的高精度絕對重力測量技術，開展技術成熟度高、結構設計簡單、便于產業化開發和臺站觀測試驗推廣的絕對重力與重力梯度一體化測量儀器的研制，可以打破國際上的行業技術壟斷，有效提升試驗區域重力網的絕對重力控制能力[51]。高精度和高時空分辨力的絕對重力與重力梯度一體化觀測，為我國專業地震臺站的重力觀測技術革新提供了新發展方向。基于地區同類觀測儀器的多源數據融合分析也為開展地球長時間尺度的潮汐響應、局部地區的斷層參數和地殼結構研究，為區域地球動力學背景和孕震機理提供全方位的重力學支撐[52-54]，有助于發現可重復出現，且物理意義明確的前兆信號，對破壞性地震的短、臨預測有 “追本溯源” 和 “臨門一腳” 的作用。