自2008年汶川大地震發(fā)生以來,我所陳祖安副研究員等人與中國地震局地球物理研究所林邦慧研究員合作,采用三維流變非連續(xù)變形(塊體邊界)與有限元(塊體內(nèi))相結(jié)合的方法(DDA+FEM),在青藏高原及其東側(cè)的四川盆地、鄂爾多斯塊體地區(qū)三維構(gòu)造塊體相互制約的大環(huán)境中,考慮了龍門山斷裂帶東西兩側(cè)地勢、地殼厚度和分層的明顯變化,以及斷裂帶東側(cè)四川盆地與鄂爾多斯塊體堅(jiān)硬地殼阻擋的影響,通過用GPS資料做位移速率邊界約束和震源機(jī)制約束,計(jì)算得到研究區(qū)的速度場和應(yīng)力場與該地區(qū)GPS測量結(jié)果和震源機(jī)制分布結(jié)果基本相似的結(jié)論。
在此基礎(chǔ)之上,他們計(jì)算了現(xiàn)今構(gòu)造塊體邊界斷層上表征剪應(yīng)力及法向應(yīng)力等綜合影響的危險(xiǎn)度分布,揭示出汶川與蘆山震源區(qū)正處在危險(xiǎn)度較大地區(qū)(見圖a,))。隨后,他們數(shù)值模擬了2008年汶川大地震的破裂過程,研究此次大震引起青藏高原及鄰區(qū)各構(gòu)造塊體邊界斷層應(yīng)力狀態(tài)變化特征,結(jié)果表明:在上地殼層汶川發(fā)震斷層南北兩個(gè)端部地段外側(cè)(即鮮水河斷裂帶南段、東昆侖斷裂帶東南段)庫侖破裂應(yīng)力增加幅度最大,達(dá)到約0.2MPa(見圖b,))。而蘆山正處于汶川發(fā)震斷層南端外側(cè)。由此說明2013年4月20日蘆山7.0級地震所處區(qū)域是滿足危險(xiǎn)度大和庫侖破裂應(yīng)力增加幅度大兩個(gè)條件的地區(qū),因而是大震發(fā)生的危險(xiǎn)區(qū)。