薛建華，劉 超

(1.山西大同大學建筑與測繪學院，山西大同037003；2.國核電力規(guī)劃設計研究院，北京100095)

重力是影響地球表層空間一切相關物理運動的基本場，重力異常作為表征區(qū)域物質分布勻質性的最優(yōu)量，已經成為地球重力場研究的重要對象。在針對全球或大范圍變化的衛(wèi)星探測技術，特別是GOCE衛(wèi)星因其大尺度，不接觸的方法，能夠采用高精度的重力梯度儀精密測定地球外部重力分布情況，更有效的為相關地球表面物理活動研究提供數據基礎[1]。

1 重力場的原理

衛(wèi)星重力梯度測量是利用低軌衛(wèi)星裝載的重力梯度儀直接測量引力梯度張量的各分量，已知地球重力場模型求解地球上某點（r,θ,λ,）（球坐標表示）的重力異常。綜合地球物理大地測量中的理論結果，重力異常帶諧系數球函數級數展開表達式為：

式中：ae為地球長半徑；N為地球重力場模型的最高階數；n和m為球諧展開的階和次；Cnm,Snm為完全規(guī)格化球諧位系數；GM為引力常數與地球質量的乘積；Pnm( )sinθ為完全規(guī)格化的締合勒讓德多項式函數[2]。

利用地球各種物理現象發(fā)生時間內重力場與背景重力場的疊積求差，可推測地震地質活動引起的質量變化。

2 算例與結果

針對2011年03月11日的里氏9.0級東日本大地震，選取北緯35～40°、東經140～145°作為研究區(qū)域，網格間隔0.1×0.1°。具體采用由ICGEM組織提供的GO CONS GCF 2 TIM R5重力場模型（完全基于GOCE觀測數據，最大280階）；對比數據選用Tongji GRACE02月模型（基于GRACE數據，最大60階）。嘗試性的計算了日本及鄰近區(qū)域重力異常同震變化結果，單位為uGal。每組計算結果分別有2 601個重力異常數據，其中GOCE計算的該區(qū)域最大值為169.870 20 mGal，最小值為-200.487 8 mGal，均值為9.392 150 3 mGal，標準偏差為87.420 20 mGal；震中位置的重力異常值為54.193 747 mGal；GRACE重力異常經統(tǒng)計顯示：最大值為74.388 265 mGal，最小值為-37.202 04 mGal，均值為7.473 010 5 mGal，標準偏差為29.625 211 mGal；震中重力異常值為12.199 234 mGal。

同時為了更好體現兩種探測技術在解算重力場結果的直觀性，通過surfer8.0軟件生成了表征方向和梯度的矢量等值線，見圖1和圖2。

圖中明顯發(fā)現，基于GOCE的重力場比基于GRACE的計算結果的變化起伏更突出，而GRACE較平緩。圖中箭頭表示重力梯度的方向，觀察震源中心緯度的點，該區(qū)域梯度值大或者方向陡峭，表明這些區(qū)域物質分布非勻質性大，表示巨大的地震往往發(fā)生在地球物理場變異帶和地殼陡變帶，等值線由密集變?yōu)槭婢徎蛐逼聨?，圖中表明此地區(qū)附近存在斷層，正是斷層出現大的位移導致地震的發(fā)生，這與一般大地震需要產生的地質條件比較吻合，說明衛(wèi)星重力異常等值線圖能夠反映地球深部構造變化特征[3]。

圖1 GOCE衛(wèi)星重力異常結果矢量等值線圖

圖2 GRACE衛(wèi)星重力異常結果矢量等值線圖

同時，為了更好的分析兩種重力異常結果的差距，將兩組結果作差并生成新的等值線，見圖3。圖中可得：同位置上兩組重力異常數值上相差較大，甚至在差異大的區(qū)域可達±100 mGal，而其他部分區(qū)域起伏不大，差異不明顯。地震震源中心差異值為42 mGal,說明本身數值不大，但卻位于差異值密集分布帶上，梯度變化大，而其他區(qū)域相對平坦。

圖3 GOCE和GRACE重力異常結果差值圖

3 結論

GRACE數據計算的重力異常結果整體呈現比較平緩趨勢，震中區(qū)域處于正負變化的邊緣帶上，說明其對地震等地球表面活動重力場變化不敏感，與實際地震變化場分布情況相關性不好，一般更好地對應于地?；蚋钐幍睾速|量分布的不均勻及不對稱性；GOCE數據計算結果中，整體變化起伏較大，區(qū)域中各塊體劃分細致明顯，震源基本處于一個變化越發(fā)突出的隆起和凹陷塊體中間，說明其間重力變化相對劇烈，與板塊的分布特點更加對稱，能更清楚地表現上地幔和地殼質量體的局部密度變化的不均勻性，這與地震主要是由斷層活動引起這一物理信息的解釋基本一致，與強烈地震的孕育、發(fā)生和發(fā)展是由于地殼中深部的物質運移所造成的這一動力學原因相對應[4-5]，驗證了GOCE衛(wèi)星更具有分析探測大地震引起重力異常變化的可行性。