走滑斷層地震地表變形與破裂寬度數值分析

李 紅鄧志輝陳連旺周 慶冉洪流邢成起中國地震局地質研究所，北京 100029北京市地震局，北京 100080廣東省地震局，廣州 510070中國地震局地殼應力研究所，北京 100085

走滑斷層地震地表變形與破裂寬度數值分析

李 紅1,2）鄧志輝3）陳連旺4）周 慶1）冉洪流1）邢成起2）
1）中國地震局地質研究所，北京 100029
2）北京市地震局，北京 100080
3）廣東省地震局，廣州 510070
4）中國地震局地殼應力研究所，北京 100085

震后野外考察表明，走滑斷層同震引起的地表變形主要集中在斷層兩側很窄的范圍，為定量研究同震地表強變形分布特征與破裂帶寬度，以鮮水河斷裂北西段為原型設計了一個簡化的三維黏彈有限元模型，斷層遵從庫侖滑動定律，斷層長度20 km，走向320°，厚度11 km，兩側各寬3 km。由已知地震斷層面的突然錯動在地表引起的位移分布，求出任意相鄰兩點之間的位移量差值，即變形量，它是依賴于兩點間距離的變量。對于空間相鄰兩點di和di+1，其水平位移差hi為hi+1－hi；垂向位移差為Zi+1－Zi。總位移差di由求得，di空間位置位于相鄰兩點di和di+1的中間，可以表示為：

以斷層相對左旋位錯模擬強震（如1973年爐霍MS7.6地震的最大地表水平位錯為3.6 m），模型南西與北東邊界在垂直斷層走向方向約束，平行斷層走向方向自由，底面垂向約束。取地表上自南西盤向北東盤經過斷層長度中點（假定為強震震中投影）O，垂直于斷層走向的AB測線（AO、BO均為100 m，沿AB單元尺度為2 m，網格為六面體）的地表變形數據，以每米0.02 m的地表變形量為破裂閾值，即2 m跨度變形量超過0.04 m時認為可導致地表破裂，換算出寬度。結果如下。

（1）不同位錯量（1.5～8.5 m，增量約為0.5 m，斷層傾向北東，傾角80°）：①AB測線的變形量曲線為約以斷層為中心的對稱單峰型，強變形主要集中在斷層兩側各50 m范圍內，達到破裂閾值的破裂寬度更窄。②不同位錯量的同震地表強變形曲線形態相似，隨位錯量增大，強變形寬度有所增加，但增加速率逐漸減小，斷層兩側50 m范圍內變形量隨位錯量增大而增大。③地表破裂寬度隨位錯量增大而增大，但斜率逐漸較小，并趨于一個漸進值，位錯量為2 m與8.5 m時，破裂寬度分別為12 m與66 m。

（2）不同斷層傾角（傾向北東，傾角為45°、50°～90°,位錯量為3.6 m，增量為10°）：傾角對變形量曲線形態影響微小，對破裂寬度影響不大，最大變形量峰值位置隨傾角減小而向北東方向略有偏移。

（3）位錯量為3.6 m時隱伏垂直走滑斷層：①上覆沉積層為松散土層時，沿AB測線的地表變形曲線的單峰隨沉積層厚度（5～40 m，厚度增量為5 m；42～45 m，增量為1 m）增加逐漸變平緩，最大變形量隨厚度增加呈現二次下凹曲線，并趨于一漸近值，破裂寬度隨厚度增加約44 m時，已不能在地表形成破裂（圖略）。②不同土質類型不同沉積厚度的情況見表1。

表1　不同土質類型在不同厚度的地震地表最大變形量與破裂帶寬度 （單位：m）