汶川地震震前龍門山以西地殼垂直隆升機理分析*

蔣鋒云 王慶良 朱良玉 張曉亮 王雙緒

(中國地震局第二監測中心，西安 710054)

汶川地震震前龍門山以西地殼垂直隆升機理分析*

蔣鋒云 王慶良 朱良玉 張曉亮 王雙緒

(中國地震局第二監測中心，西安 710054)

結合流動GPS觀測速度場及層析成像結果，構建跨龍門山斷裂剖面的二維有限元分層模型，分兩種情況討論汶川地震前龍門山前緣地殼垂直隆升的物理機制，以及中、下地殼軟物質垂向和橫向的不均勻性對地殼隆升作用的影響。分析認為：川西高原相對四川盆地的差異抬升和龍門山以西地殼縮短的共同作用是汶川地震震前龍門山前緣地殼垂直隆升的可能原因。

地殼隆升;低速體;黏彈性有限元;GPS速度;川西高原

1 引言

汶川8.0級地震發生以來，許多專家和學者結合地形變資料和其他地球物理資料對其孕震機制展開了大量的研究，獲得了豐碩的研究成果［1-6］。王慶良等［1，2］的研究表明：從阿壩到映秀有一個明顯的地表隆起，在理縣附近隆起量最大，達3.5 mm/a。張竹琪等［3］利用汶川地震近場同震水準變形資料反演了發震斷層隨深度分布的二元結構特征，而有限元模擬研究表明，松潘地塊的地表垂直變形在很大程度上并不是地殼縮短的結果，而是青藏高原東緣相對四川盆地的差異抬升。上述研究提出了松潘地塊隆起的可能原因，但在有限元模型的建立上缺乏和地球物理資料的結合，沒有考慮中下地殼物質的流變特性及橫向不均勻性，模型相對簡單。本文將結合GPS資料、層析成像結果［?］，考慮中下地殼介質垂向和橫向分布不均勻性及流變特征，細化有限元模型，定量分析汶川地震震前跨龍門山水準剖面垂直形變機理。

2 資料分析

從流動GPS觀測得到的地殼水平形變速度場(圖1)可以看出，從震源區開始沿龍門山斷層向北東方向，龍門山斷層呈現明顯的由以逆沖為主轉變為兼具走滑的運動特征，同時可以看出平行于斷層和垂直斷層的GPS速度分量越靠近龍門山量值越小。從垂直跨龍門山斷裂的GPS剖面(圖2)可以看出，垂直于斷裂方向，離斷裂以西200 km范圍內，每年地殼有3～4 mm的縮短，200～400 km地殼縮短不明顯。

滕吉文［4］認為在下地殼的上部和上地殼底部存在約10 km厚的低速層，這些低速體從遠處流向龍門山，并在龍門山前緣堆積。吳建平等人［5］的地震層析成像結果顯示，龍門山斷層前緣附近，下地殼存在低速體的堆積，但在稍遠的地方看不到低速層。

圖1 龍門山地區GPS水平運動速率場(相對歐亞板塊)Fig.1 GPS velocity field of the Longmenshan zone(relative to Eurasia plate)

圖2 垂直龍門山斷裂GPS水平運動分量Fig.2 Horizontal GPS velocity component perpendicular to Longmenshan fault

3 有限元模擬

針對上述情況，結合前人研究結果，參照GPS速度場和長水準剖面測線，選取跨龍門山斷裂的垂直剖面建立二維有限元模型，模型西邊界遠離斷層400 km，東邊界距斷層100 km。模型0 km處為汶川地震發震斷層所在位置。假設研究區巖石圈接近重力均衡，即模型不考慮重力作用。有限元剖面模型參數參照張竹琪［3］和朱守彪［6］等研究成果取值。模型各層參數如下：

模型深120 km，上地殼厚為15 km，低速層厚8 km，下地殼龍門山以西的厚度為35 km，以東厚20 km(圖3);龍門山以西上地殼視為彈性體，彈性模量為72 000 Mpa。龍門山附近及四川盆地上地殼視為彈性體，彈性模量為80 000 Mpa。低速層為黏彈體彈性體，模量為70 000 Mpa，黏滯系數為5.0× 1018pa/s;龍門山以西下地殼視為彈性體，彈性摸量為74 000 Mpa。龍門山以東下地殼視為彈性體，彈性摸量為84 000 Mpa;上地幔為黏彈體，彈性模量為88 000 Mpa，黏滯系數為1.0×1021pa/s;低速體為黏彈體，彈性模量為65 000 Mpa，黏滯系數為5.0× 1017pa/s;低速體的泊松比為0.28外，其余的泊松比為0.25。下地殼和上地幔均勻垂直向上的體力為0.08 N/m3。

通過總結分析前人的研究結果［7，8］，分兩種情況模擬震前的垂直變形：1)只考慮垂直于斷層水平方向的地殼縮短。即考慮遠場效應情況下，以每年6 mm的量值均勻加載于模型西邊界。模型底部垂直方向固定不動，水平方向自由滑動。模型東邊界(四川盆地)垂直方向可以自由滑動，水平方向固定不動。地表為自由表面;2)在考慮東西方向地殼水平縮短的同時，在龍門山以西下地殼和上地幔加入均勻向上的體力(0.08 N/m3)，突出青藏高原東邊緣與四川盆地的差異抬升。邊界條件和第一種情況一致。

圖3 跨龍門山斷裂垂直剖面有限元模型示意圖Fig.3 Sketch of FEM model of the vertical profile across Longmenshan fault

4 模擬結果及分析討論

模擬結果如圖4和5所示，從圖4和5可以看出，只考慮地殼縮短所引起的地殼垂直變形量值(最大0.5～0.6 mm/a)遠遠小于觀測得到的震前地殼垂直形變量值(最大3.5 mm/a)。但模擬所得龍門山斷裂前緣地殼隆起的形態(圖4)和震前水準觀測結果基本相似。這說明以目前GPS觀測得到的地殼水平運動速率量值，僅僅靠垂直于斷層的水平向地殼縮短無法使龍門山以西地殼相對四川盆地的隆起速率達到實際觀測的量值。但水平向地殼縮短對龍門山以西地殼垂直隆升形態的刻畫起到了重要作用。既考慮地殼縮短又考慮青藏高原東邊緣與四川盆地的差異抬升(在龍門山以西下地殼和上地幔垂直方向給定一大小為0.08 N/m3均勻向上的體力)，產生的地殼垂直變形(圖5)在量值上和觀測得到結果較為一致，地殼隆起形態和觀測結果也更為吻合。這說明川西高原相對四川盆地的差異抬升可能是龍門山以西地殼垂直隆升的主要原因。因此我們認為：汶川地震震前龍門山以西的地殼隆起主要是川西高原相對四川盆地抬升的結果。同時在地殼垂直隆起形態上也受到垂直于斷裂方向的地殼縮短的影響。汶川8.0級地震發生時，龍門山斷層兩盤產生的巨大的同震垂直位移，是川西高原相對四川盆地長期垂直抬升運動在龍門山斷層上產生的能量積累釋放的結果。

圖4 只考慮地殼縮短有限元計算垂向運動速率Fig.4 Vertical velocity calculated with the FEM of crustal shortening considered solely

圖5 既考慮地殼縮短又考慮川西高原相對四川盆地的差異抬升有限元計算垂向運動速率Fig.5 Vertical velocity calculated with the FEM both crustal shortening and the vertical uplift between Chuanxi plateau and Sichuan basin considered

為了考察龍門山前緣中下地殼軟物質對垂直形變的影響，再做以下有限元的模擬，不考慮低速層和低速體的黏彈性流變特性，將龍門山以西的低速層、下地殼和低速體統一都看做彈性體，彈性模量為70 000 Mpa，泊松比為0.25，邊界條件和上面第二種情況相同。模擬發現，要使龍門山前緣地殼垂直隆升速率達到3.5 mm/a(圖6)，所加載的龍門山以西下地殼和上地幔的垂直向上的體力大小為0.16 N/ m3，這和張竹琪等［3］的結果一致，是我們前面研究結果的兩倍。另外，從圖7可以看出，模擬的地殼垂直隆升的曲線形態雖然與龍門山地區的地表形態較為接近，但和實際觀測結果在局部地區不太相符，特別是龍門山以西較遠的區域差異較大。這可能說明從一個較長時尺度和一個較大的空間范圍來看，青藏高中下地殼原普遍存在的軟物質層對青藏高原西部的垂直隆升形態的塑造作用不明顯。更重要的影響來自局部地區軟物質的不均勻分布與局部動力環境對地殼細部隆起形態的刻畫。

圖6 不考慮下龍門山以西地殼介質橫向不均勻性有限元計算所得垂向運動速率Fig.6 Vertical velocity calculated by the FEM without considering the latered inhomogereity in the low crust in the west of Longmenshan

5 結語

通過以上的分析，可以看出引起地殼垂直隆升的原因是復雜的，它與地殼深部的動力環境、地殼物質的密度分布，及地殼淺層的構造運動和斷層的幾何形態等方面都有關系。只有把這幾個方面綜合起來考慮，分析才能得到較為可信的結果。本文就是將跨龍門山斷裂的地球物理資料和形變資料，同時考慮斷層的幾何形態及該區域動力構造環境，在已有研究的結果上，通過數值模擬進一步探討汶川8.0級地震震前龍門山前緣地殼垂直隆升的可能原因，認為地殼垂直隆升與龍門山前緣中下地殼的軟物質、垂直于斷裂方向的地殼縮短及川西高原相對四川盆地的差異抬升都有關系。中、下地殼軟物質和地殼縮短刻畫了地殼垂直隆升剖面曲線形態特征。汶川地震震前龍門山以西地殼垂直隆升的主要原因可能是龍門山以西中下地殼相對四川盆地存在垂直向上的力源。地殼縮短對垂向運動也有一定的貢獻，但量值不大。

1 王慶良，等.龍門山及汶川Ms8.0級地震垂直形變場研究［J］.國際地震動態，2010，(6)：11.(Wang Qingliang，et al.Studies on surface vertical deformation of the Longmenshan in the Wenchuan Ms8.0 earthquake［J］.Recent Developments In World Seismology，2010，(6)：11)

2 Wang Q L，Cui D X and Zhang X.Coseismic vertical deformation of theMs8.0 Wenchuan earthquake from repeated levelings and its constraint on listric fault geometry［J］.Earth Sciences，2009，22：595-602)

3 張竹琪，張培震，王慶良.龍門山高傾角逆斷層結構與孕震機制［J］.地球物理學報，2010，53(9)：2 068-2 082.(Zhang Zhuqi，Zhang Peizhen and Wang Qingliang.The structure and seismogenic mechanism of Longmenshan High dip-angle reverse fault［J］.Chinese J Geophys.，2010，53 (9)：2 068-2 082)

4 滕吉文，等.汶川Ms8.0地震發生的深層動力學過程和動力學響應［J］.地球物理學報，2008，51(5)：1 385-1 402.(Teng Jiwen，et al.Deep processes and dynamic responses associated with the Wenchuan Ms8.0 earthquake of 2008［J］.Chinese J Geophys.，2008，51(5)：1 385-1 402)

5 吳建平，等.汶川Ms8.0級地震余震分布及周邊區域P波三維速度結構研究［J］.地球物理學報，2009，52(2)：321-328(Wu Jianping，et al.Aftershock distribution of the Ms8.0 Wenchuan earthquake and three dimensional P-wave velocity structure and around source region［J］.Chinese J Geophys.，2009，52(2)：321-328)

6 朱守彪，張培震.2008年汶川8.0地震發生過程的動力學機制研究［J］.地球物理學報，2009，52(2)：418-427.(Zhu Shoubiao and Zhang Peizhen.A study on the dynamical mechanisms of the Wenchuan Ms8.0 earthquake，2008［J］.Chinese J Geophys.，2009，52(2)：418-427)

7 Hubbard J and Shaw J.Uplift of the Longmenshan and Tibetan plateau，and the 2008 Wenchuan(M=7.9)earthquake［J］.Nature，2009，458：194-197.

8 Burchfiel B C，et al.A geological and geophysical context for the Wenchuan earthquake of 12 May 2008，Sichuan，People Republic of China［J］.GSA Today，18(7)：4-11.

MECHANISM ANALYSIS OF VERTICAL UPLIFT OF CRUST IN THE WEST OF LONGMENSHAN BEFORE WENCHUAN Ms8.0 EARTHQUAKE

Jiang Fengyun，Wang Qingliang，Zhu Liangyu，Zhang Xiaoliang and Wang Shuangxu (Second Crust Monitoring and Application center，CEA，Xi’an 710054)

The mechanism of the vertical uplift of crust in the west of Longmenshan before Wenchuan Ms8.0 earthquake has been debated through constructing the 2-D FEM model for the profile across Longmenshan fault，which is based on the velocity field that had calculated from the GPS observations and the result of seismic tomography obtained by Wu Jiangping et al.In addition，the effect of crustal inhomogeneity on the vertical uplift has been discussed as well.We concluded that the existence of the low velocity anomaly in the low cust below the front of Longmenshan fault play a key role，and vertical uplift is likely caused by both the curst shorten of the westen Longmenshan and the differential uplift between Chuanxi pletau and Sichuan basin.

crust uplift：low velocity body：viscoelastic Finite Element Method(FEM)：GPS velocity;Chuanxi plateau

1671-5942(2011)05-0026-04

2011-05-21

中國地震局監測預報司“南北地震帶北段大震危險性強化跟蹤研究”

蔣鋒云，男，1978年生，工程師，主要研究方向：地震預報與GPS數據處理.E-mail：jfy267862@163.com

P315.71+6

A