六盤山地區地殼重力均衡與隆升機制研究

佘雅文付廣裕蘇小寧孟國杰祝意青郭樹松徐云馬中國地震局地震預測研究所地震預測重點實驗室，北京 100036中國地震局第二監測中心，西安 710054

六盤山地區地殼重力均衡與隆升機制研究

佘雅文1）付廣裕1）蘇小寧1）孟國杰1）祝意青2）郭樹松2）徐云馬2）
1）中國地震局地震預測研究所地震預測重點實驗室，北京 100036
2）中國地震局第二監測中心，西安 710054

近年來很多學者對六盤山地區的隆起機制給出了不同的見解，為了探明這一問題，利用重力數據對六盤山地區的隆起機制進行研究。通過對該地區實測重力數據與模型重力數據的分析，發現六盤山地區的隆起主要是由于青藏高原東北緣的向外延伸而造成的。

2014年8月跨越六盤山東麓斷裂帶，在長度約為200 km的剖面上展開了流動重力/GPS聯合剖面觀測。觀測結果表明，測線東端的布格重力異常約為-1.9×10-3ms-2，西端則為-2.5×10-3ms-2左右。在假設地殼均一的前提下，基于Airy均衡模型，利用布格重力異常和GPS觀測數據，分別計算了測線所在剖面的莫霍面深度與均衡面深度。發現六盤山地區均衡深度大于莫霍面深度，表明該區域處于正均衡異常狀態。考慮地殼分層的前提下，以研究區域地殼P波波速結構為約束，使用布格重力異常數據反演獲取六盤山地區的地殼密度結構（圖1）。利用上述地殼密度結構，依舊基于Airy均衡模型，計算六盤山地區均衡面與莫霍面深度，對比顯示該區域亦處于均衡正異常狀態。

假設地殼為彈性板的前提下，利用EGM2008自由空氣重力異常和SRTM V18.1 DEM數據，在頻域下研究了以六盤山為中心的青藏高原東北緣地區的巖石圈有效彈性厚度（Te）和加載機制。使用的方法為自由空氣重力異常導納方法。此方法利用MTM（多窗譜分析方法）對重力數據和地形數據進行頻譜估計，進而獲取自由空氣重力異常導納。使用此重力導納與彈性板模型導納在波數域（k=2π/λ，k為波數，λ為波長）中擬合，從而求解巖石圈撓曲的初始加載位置和有效彈性厚度。計算結果表明六盤山地區的Te為5 km，巖石圈加載主要來自于地表，占總加載的95%。依據獲取的Te和加載狀態，計算研究區域的布格模型導納。使用此導納和頻域地形數據可計算研究區域的頻域布格重力異常，再進行傅里葉逆變換，可得到區域的彈性板布格重力異常，進而獲取彈性板均衡異常。此異常表明六盤山地區在假設地殼為彈性板的前提下基本處于均衡狀態。

Airy均衡異常與彈性板均衡異常的差異即是地殼分別為塑性體和彈性體的差異，此差異代表彈性地殼對其自身的承載。對比六盤山地區的Airy均衡異常和彈性板均衡異常的差異，發現六盤山東麓斷層處地殼承載梯度值最大（圖2），表明該斷裂帶吸收了較多的應變能。

圖1　跨六盤山地區布格重力異常與地殼密度結構反演結果

圖2　彈性板模型與Airy模型均衡異常的對比