鄧志輝 胡勐乾 周 斌 陸遠(yuǎn)忠陶京玲 馬曉靜 姜 輝 李 紅

1)中國地震局地質(zhì)研究所，北京 100029

2)中國地震局地殼應(yīng)力研究所，北京 100085

數(shù)值模擬方法在地震預(yù)測研究中應(yīng)用的初步探討(Ⅱ)

鄧志輝1)胡勐乾1)周 斌1)陸遠(yuǎn)忠2)陶京玲1)馬曉靜1)姜 輝1)李 紅1)

1)中國地震局地質(zhì)研究所，北京 100029

2)中國地震局地殼應(yīng)力研究所，北京 100085

在地震危險性分析和預(yù)測研究中，隨著數(shù)值模擬技術(shù)的進(jìn)步和對地球動力學(xué)過程認(rèn)識的不斷深入，有可能發(fā)揮越來越重要的作用。數(shù)值模擬一個區(qū)域內(nèi)部的構(gòu)造變形，在正常時期與觀測結(jié)果應(yīng)該是比較接近的，而當(dāng)局部地方處于孕震臨界失穩(wěn)狀態(tài)時，其實(shí)測結(jié)果就有可能與模擬結(jié)果不協(xié)調(diào)，據(jù)此有可能為地震的中期至短期預(yù)測提供有價值的信息;應(yīng)用數(shù)值模擬方法研究地震庫侖應(yīng)力變化，可模擬分析各種因素對庫侖應(yīng)力變化的大小和分布的影響，獲得更接近實(shí)際的結(jié)果;計(jì)算機(jī)模擬合成地震目錄是解決建立區(qū)域地震預(yù)測概率模型數(shù)據(jù)不足的有效方法，未來的發(fā)展方向應(yīng)該是建立更接近實(shí)際的三維動力學(xué)模型，同時考慮多場耦合作用，模擬斷層活動的規(guī)律。

地震預(yù)測 數(shù)值模擬 庫侖應(yīng)力 發(fā)震概率模型 變形不協(xié)調(diào)

0 引言

在地震危險性分析和預(yù)測研究中，隨著數(shù)值模擬技術(shù)的進(jìn)步和對地球動力學(xué)過程認(rèn)識的深入，研究對象不斷增加，研究思路不斷創(chuàng)新，研究問題不斷深化。本文在數(shù)值模擬方法在地震預(yù)測研究中應(yīng)用的初步探討(Ⅰ)的基礎(chǔ)上，進(jìn)一步討論數(shù)值模擬方法在地震預(yù)測研究的變形不協(xié)調(diào)分析、庫侖應(yīng)力變化分析和發(fā)震概率模型建立方面的可能應(yīng)用。

1 變形不協(xié)調(diào)分析

地震孕育是一個相當(dāng)長的時間過程，在相對短的觀測時間段內(nèi)，很難觀測到能量積累的顯著變化。實(shí)驗(yàn)結(jié)果表明，巖石變形進(jìn)入彈性屈服階段以后，系統(tǒng)處于臨界失穩(wěn)狀態(tài)，各質(zhì)點(diǎn)的運(yùn)動將可能偏離正常的軌跡，變得復(fù)雜化。對一個區(qū)域構(gòu)造變形的動力學(xué)模擬，其結(jié)果主要反映正常的應(yīng)力、應(yīng)變和位移等的時空變化。如果以區(qū)域運(yùn)動背景為邊界條件，模擬區(qū)域內(nèi)部的構(gòu)造變形，在正常時期與觀測結(jié)果應(yīng)該是比較接近的，而當(dāng)局部地方處于孕震臨界失穩(wěn)狀態(tài)時，其實(shí)測結(jié)果就有可能與模擬結(jié)果不協(xié)調(diào)，據(jù)此有可能為地震的中期至短期預(yù)測提供有價值的信息。本文以2006年7月4日發(fā)生在華北地區(qū)的文安5.1級地震為例說明分析過程。

根據(jù)中國活動地塊劃分及活動斷裂分布，結(jié)合地殼上地幔的波速結(jié)構(gòu)、地震活動和GPS資料等，確定模型的幾何邊界范圍是 99.8°～121.4°E，27.9°～42.3°N，區(qū)域總面積約為3,840,000km2(圖1)。研究區(qū)域包含的斷裂分布非常密集，為了更加準(zhǔn)確地模擬活動斷裂的運(yùn)動特征，模型中加入了已知的所有活動斷裂。由于斷層空間分布復(fù)雜，在建模過程中充分考慮了斷層的走向及延伸的幾何形態(tài)，對其局部的細(xì)節(jié)進(jìn)行一定的簡化。另外，華北地區(qū)的斷裂大多傾角較大，為了簡化模型及工作量，將所有斷層按垂直統(tǒng)一處理。

圖1 研究區(qū)活動斷裂分布與介質(zhì)分區(qū)圖Fig.1 Distribution of active faults and crustalmedia in study area.

根據(jù)Zhang等(2003)對中國大陸活動地塊的劃分，模型的介質(zhì)參數(shù)在每一層上劃分為8個區(qū)域:鄂爾多斯地塊、華北平原地塊、魯東-黃海地塊、阿拉善地塊、青藏地塊、燕山地塊、華南地塊以及地塊邊界帶部分(圖1)。模型在沿縱向上分為水平的4層，包括上地殼、中地殼、下地殼以及部分上地幔(嘉世旭等，2005)，總厚度為100km。利用Huang等(2003)計(jì)算得到的1°×1°中國大陸地區(qū)地殼上地幔波速結(jié)構(gòu)，對數(shù)據(jù)進(jìn)行篩選，求出平均值，計(jì)算得到上述8個區(qū)域分別在上地殼、中地殼、下地殼及上地幔4層上的楊氏模量 (表1)。在模型中將斷層處理為軟弱帶，取地塊介質(zhì)的1/10。地塊的泊松比取0.25，斷層帶的泊松比為0.33。

有限元模型的網(wǎng)格劃分采用ANSYS軟件中的solid 187三維實(shí)體單元，為保證模型網(wǎng)格密度和計(jì)算質(zhì)量，單元邊長＜25km，模型共劃分為單元416,582個，節(jié)點(diǎn)582,392個。

利用GPS測量數(shù)據(jù)作為數(shù)值模型的邊界條件已有眾多研究者進(jìn)行了嘗試，并取得了一定的成果(劉峽，2006;曹建玲等，2009)。本文以1999—2004年 GPS速度場數(shù)據(jù)(地殼運(yùn)動第一觀測中心提供)為邊界條件，它是由1999、2001和2004年3期的觀測結(jié)果得到的運(yùn)動速度。為了得到較好的GPS插值結(jié)果，選取了研究區(qū)及其鄰區(qū)共979個GPS觀測臺站的數(shù)據(jù)進(jìn)行插值計(jì)算，得到有限元模型邊界節(jié)點(diǎn)的運(yùn)動速度約束條件。模型的底面在垂向上固定約束。

圖2為以1999—2004年GPS插值結(jié)果為邊界運(yùn)動速度約束條件下，模擬得到的地表運(yùn)動速度場和GPS站點(diǎn)的觀測結(jié)果之間的對比。從圖2中可以看出，數(shù)值模擬運(yùn)動速度與GPS觀測運(yùn)動速度大部分比較吻合，反映了比較好的一致性。模擬結(jié)果顯示，地殼運(yùn)動速率的整體分布呈現(xiàn)由東向西逐漸減小、自北向南逐漸增大的特點(diǎn)。華北地區(qū)北部、燕山地塊的地殼運(yùn)動速率總體較小，＜2mm/a，而在華北地區(qū)中部、南部以及華南地塊，運(yùn)動速率逐漸增大。在燕山地塊，地殼運(yùn)動速度方向總體為SW向，而在華北地區(qū)及華南地塊地殼運(yùn)動速度方向逐漸轉(zhuǎn)為SE向。

模擬結(jié)果與GPS觀測結(jié)果在個別區(qū)域也存在明顯的不協(xié)調(diào)現(xiàn)象，其中，差異較大的點(diǎn)主要分布在華北平原塊體的東北部(圖2)，其原因可能與2006年7月4日文安5.1級地震的孕育和發(fā)生有關(guān)。

2 庫侖應(yīng)力變化分析

Chinnery(1963)對走滑斷層相應(yīng)位移場的應(yīng)力分布的研究，拉開了地震應(yīng)力觸發(fā)研究的序幕。應(yīng)用靜態(tài)庫侖破裂應(yīng)力模型研究地震間應(yīng)力觸發(fā)作用是當(dāng)前地震研究的熱點(diǎn)問題。國內(nèi)外許多學(xué)者在強(qiáng)震對后續(xù)地震的觸發(fā)作用、靜態(tài)庫侖破裂應(yīng)力與余震發(fā)生、靜態(tài)庫侖破裂應(yīng)力與區(qū)域地震活動的關(guān)系等方面做了大量工作，取得了很多研究成果，大地震產(chǎn)生的庫侖破裂應(yīng)力變化對后續(xù)地震活動影響的研究受到廣泛關(guān)注(Harris，1998)。許多震例表明，大地震后庫侖應(yīng)力的變化與后續(xù)地震活動有促進(jìn)作用，很小的靜態(tài)庫侖應(yīng)力變化就可能對地震活動性、地震的時空分布產(chǎn)生較大的影響(Deng et al.，1997)。大多數(shù)庫侖破裂應(yīng)力增加的地區(qū)地震活動增強(qiáng);而庫侖破裂應(yīng)力減小的地區(qū)，地震活動減弱(Reasenberg et al.，1992;Hardebeck et al.1998;Toda，1998;劉桂萍等，2002)。Stein等(1997)計(jì)算了在1939—1992年發(fā)生在土耳其North Anatolian斷裂帶10個6.7級以上地震的庫侖破裂應(yīng)力演化過程，發(fā)現(xiàn)其中90%的地震是被先前地震觸發(fā)的。庫侖破裂應(yīng)力變化起觸發(fā)作用的典型值是0.1～1MPa，也有一些地區(qū)庫侖破裂應(yīng)力增加量＞0.01～0.03MPa就有可能觸發(fā)地震活動(Parsons，2002)。陳連旺等(2008)通過建立川滇地區(qū)的三維有限元模型，數(shù)值模擬了川滇地區(qū)1981—2000年MS≥6.5的地震序列，結(jié)果顯示，后續(xù)地震全部位于已發(fā)生地震所引起的庫侖破裂應(yīng)力正值區(qū)。庫侖應(yīng)力對地震的觸發(fā)作用有一個延遲時間，延遲時間可為幾秒到幾十年。

圖2 數(shù)值模擬地表速度場與GPS觀測速度場對比圖①胡勐乾，2010，并行計(jì)算三維數(shù)值模擬在華北地區(qū)現(xiàn)今構(gòu)造變形分析中的應(yīng)用研究，碩士學(xué)位論文。Fig.2 Comparison of the numerical simulation speed and GPS velocity.

水庫誘發(fā)地震也是庫侖應(yīng)力變化觸發(fā)地震的很好例證。這里以四川紫坪鋪水庫為例，說明庫侖應(yīng)力變化對地震活動的影響。在垂直于庫區(qū)主構(gòu)造線的方向，選用二維地質(zhì)剖面作為有限元數(shù)值模擬實(shí)驗(yàn)的地質(zhì)模型，共劃分出12,513個三角形單元，6,434個節(jié)點(diǎn)。12種地質(zhì)體力學(xué)參數(shù)與物性參數(shù)的選取主要參考了中石化南方勘探開發(fā)分公司在龍門山地區(qū)及四川盆地開展的巖石力學(xué)參數(shù)與物性參數(shù)測試成果，以及巖石力學(xué)手冊中研究區(qū)附近的映秀灣水電站和四川“紅層”地層的實(shí)測數(shù)據(jù)(葉金漢等，1991)(表2)。地質(zhì)構(gòu)造主要參考中石化CX-NW-03-167、CX-NW-03-161和CX-NW-05-26地震解釋剖面(金文正等，2007)。

依據(jù)2005年3月31日至2008年4月14日期間壩前水位變化資料，將相對于2005年3月31日基準(zhǔn)水位(海拔高度約754.5m)的水位變化，用分段直線加載函數(shù)在庫體處施加動水頭邊界和動荷載邊界 (圖3)。

圖3 水位變化分段直線加載方式①周斌，2010，水庫誘發(fā)地震時空演化及其機(jī)理研究，博士學(xué)位論文。Fig.3 Load-unloaded process in different period of reservoir water change.

圖4 是數(shù)值模擬2005年3月31日水庫開始蓄水，60d、570d、990d和1,140d后有效附加應(yīng)力(ΔCFS)水平分量和垂直分量的剖面圖。從圖4中可以看出，隨著水庫蓄水時間的延長，孔隙壓擴(kuò)散對有效附加應(yīng)力場的影響可逐漸擴(kuò)大到地下10km以上的范圍。圖5為紫坪鋪水庫水位變化與各分區(qū)地震頻次統(tǒng)計(jì)。在4個典型的區(qū)域(Ⅰ，Ⅱ，Ⅲ，Ⅳ)有效附加應(yīng)力狀態(tài)表現(xiàn)出不同的特征:

圖4 紫坪鋪水庫(位于Ⅰ區(qū)上方)蓄水不同階段有效附加應(yīng)力(ΔCFS)水平分量(左圖)與垂直分量(右圖)剖面圖①同674頁①。Fig.4 Distribution of additional effective stress field at different stages of impounding.a 60d后;b 570d后;c 990d后;d 1140d后

(1)Ⅰ區(qū)位于庫底正下方、通濟(jì)場滑脫斷裂面以上的區(qū)域，巖體破碎程度高，次級斷裂發(fā)育且與水庫有直接的水力聯(lián)系。該區(qū)地震活動對庫水加卸載過程的時間響應(yīng)表現(xiàn)出伴隨著庫水位的高低變化，小震發(fā)生頻度呈同步變化的特點(diǎn)。經(jīng)過長時間以后，該區(qū)的水平向和垂向有效附加應(yīng)力分量逐漸由壓應(yīng)力狀態(tài)轉(zhuǎn)變?yōu)閺垜?yīng)力狀態(tài)，該區(qū)可能會誘發(fā)新的微震活動。

(2)Ⅱ區(qū)位于安縣－灌縣斷裂的上盤，發(fā)育三疊系須家河組陸相碎屑巖和雷口坡組—嘉陵江組膏巖，巖性軟弱，滲透條件差。ΔCFS(特別是其垂直分量)總體上在負(fù)值區(qū)間變化，庫體荷載作用對巖體和斷層的穩(wěn)定起到了決定性的作用，因而該區(qū)小震活動非常之少。

圖5 紫坪鋪水庫水位變化與各分區(qū)地震頻次統(tǒng)計(jì)①同674頁① 。Fig.5 Relationship between the water level change and earthquake frequency in different areas.

(3)Ⅲ區(qū)位于安縣-灌縣斷裂與通濟(jì)場斷裂在地腹的交會區(qū)域，主要發(fā)育泥盆系—下三疊統(tǒng)飛仙關(guān)組海相碳酸鹽巖，巖性堅(jiān)硬，裂隙相對發(fā)育，滲透條件好。該區(qū)遠(yuǎn)離庫體覆蓋區(qū)，庫體荷載作用對斷層穩(wěn)定性的影響已不占主導(dǎo)地位。由于NW向張扭性斷裂的溝通使該區(qū)深部的斷裂與地表水體之間建立了良好的流體運(yùn)移通道。在2005年9月30日至2006年11月30日壩前水位高幅、快速上升時段，無論是有效附加應(yīng)力場的水平向分量還是垂向分量都快速地從壓應(yīng)力狀態(tài)轉(zhuǎn)變?yōu)閺垜?yīng)力狀態(tài)，表明了附加水頭壓力擴(kuò)散對該區(qū)有較大的影響，該區(qū)小震活動增強(qiáng)現(xiàn)象說明了孔隙壓力擴(kuò)散作用對該區(qū)誘發(fā)地震活動的促進(jìn)作用。

(4)Ⅳ區(qū)位于廣元-大邑?cái)嗔押蜕畈吭財(cái)鄬咏粫膮^(qū)域。由于該區(qū)距離庫底較遠(yuǎn)，在庫體荷載作用下產(chǎn)生的彈性附加應(yīng)力場本身很微弱，有效附加應(yīng)力場的變化主要受控于孔隙壓擴(kuò)散過程，在2008年2—4月期間，ΔCFS上升到0.1MPa左右，該區(qū)相繼發(fā)生了都江堰震群活動，依此判斷該震群活動可能與水庫蓄水有一定的關(guān)系。

隨著研究的深入，不同的作者也發(fā)現(xiàn)了地震應(yīng)力觸發(fā)模型存在許多問題，1989年Loma Prieta地震發(fā)生后，Beroza等(1993)的研究表明，余震大部分出現(xiàn)在主震產(chǎn)生的庫侖應(yīng)力變化的影區(qū)中，即庫侖應(yīng)力變化為負(fù)的區(qū)域。萬永革等(2002)從全球角度出發(fā)，研究了不同類型斷層的靜態(tài)應(yīng)力觸發(fā)效應(yīng)，發(fā)現(xiàn)非走滑斷層地震具有明顯的觸發(fā)效應(yīng)，而走滑斷層地震觸發(fā)效應(yīng)不明顯。王凱英(2009)也得到相似的結(jié)果。Heaton(1982)的研究發(fā)現(xiàn)，潮汐應(yīng)力對淺源(＜30km)的傾滑型地震和斜滑型地震具有明顯的觸發(fā)作用，而對走滑型斷層沒有觸發(fā)作用。

有的作者還研究了影響計(jì)算結(jié)果的因素。王瓊等(2006)分析了不同震中定位精度對靜態(tài)庫侖應(yīng)力觸發(fā)的影響，得到了定位精度直接影響判定結(jié)果的結(jié)論，認(rèn)為選擇定位精度高的資料才能提高判定的客觀性和正確性。近年來，不少學(xué)者致力于更精確地確定地震斷層面參數(shù)的研究，如萬永革等(2008)研究了利用小震密集程度求解主震斷層面走向、傾角和位置，研發(fā)了根據(jù)小震位置確定斷層面走向、傾角和位置，以及根據(jù)局部應(yīng)力場參數(shù)進(jìn)一步確定斷層上滑動角的程序。并將所研究方法應(yīng)用到具有一定線性分布的叢集地震實(shí)際震例中，取得了較好的效果。石耀霖等(2010)也發(fā)現(xiàn)震后主應(yīng)力方向可能改變，對剪應(yīng)力變化量的計(jì)算結(jié)果有影響。

數(shù)值模擬可在地震預(yù)測研究的地震庫侖應(yīng)力變化分析中發(fā)揮重要作用。在建立三維動力學(xué)模型的基礎(chǔ)上，可模擬分析各種因素對庫侖應(yīng)力變化的大小和分布的影響，獲得更接近實(shí)際的結(jié)果。

3 地震發(fā)生概率模型建立

美國南加州地震中心發(fā)展了一個叫“主模型”的新概念，也可稱為地球巖石層的模型，希望能像氣象學(xué)家通過觀測氣壓、溫度、風(fēng)速等，用計(jì)算機(jī)進(jìn)行天氣預(yù)報(bào)一樣，使用“主模型”能夠依據(jù)各種地球物理觀測數(shù)據(jù)計(jì)算構(gòu)造應(yīng)力的時空分布，并以此來預(yù)測地震的發(fā)生。就目前的工作而言，南加州地震中心把建立區(qū)域地震預(yù)測概率模型(RELM-Regional Earthquake Likelihood Models)作為研究目標(biāo)?！爸髂Ｐ汀备拍畹淖饔檬亲鳛橐粋€框架來整合南加州地震研究的多學(xué)科資料，盡其所能給予科學(xué)家在創(chuàng)造性科學(xué)中的領(lǐng)導(dǎo)權(quán)，同時會給公眾帶來最好的結(jié)果(Aki，2002)。這一模型和工作方式值得中國地震研究工作借鑒。

地震預(yù)測概率模型是基于“彈性回跳理論”的有關(guān)物理模型(特征地震模型、準(zhǔn)周期模型、時間可預(yù)報(bào)模型等)，采用隨離逝時間增長的發(fā)震概率和預(yù)測的震級共同描述斷裂段的地震危險性，與現(xiàn)階段人們對地震復(fù)發(fā)機(jī)理的普遍理解一致，因此在國內(nèi)外都得到廣泛應(yīng)用(WGCEP，1995，2003;聞學(xué)澤，1990;張培震等，1996;聞學(xué)澤，1998;冉洪流等，2006;張效亮等，2009)。當(dāng)然，目前的地震預(yù)測概率模型仍然存在不少問題，需要不斷檢驗(yàn)和完善，而數(shù)據(jù)的不完備是進(jìn)行這些工作的最大難題。

應(yīng)用地震預(yù)測概率模型計(jì)算某一研究區(qū)的復(fù)發(fā)條件概率時，關(guān)鍵問題是平均復(fù)發(fā)間隔的確定精度(聞學(xué)澤，1990)。而解決這一問題的關(guān)鍵是資料的完備性，以往地震復(fù)發(fā)間隔的計(jì)算大都是基于活斷層地質(zhì)學(xué)定量資料，包括歷史地震資料、古地震、滑動速率、GPS觀測結(jié)果等，然而，對于歷史地震資料缺乏和活斷層定量研究程度較低的地區(qū)，歷史地震資料大震記錄太少，一般不能完全反映出大震的復(fù)發(fā)規(guī)律，古地震事件的揭露往往具有不完整性，滑動速率受野外地質(zhì)條件和測年精度的限制，GPS觀測資料時間太短或站點(diǎn)密度不夠。

為了解決大震統(tǒng)計(jì)樣本不足的問題，Ward(1996，1997，2000)利用計(jì)算機(jī)模擬合成加州地區(qū)5 000a的地震目錄。Ward(2000)在編制南加州斷層系統(tǒng)地震概率分布圖時，針對資料完備性的差異問題，分4類資料計(jì)算了4種結(jié)果，它們分別是完全基于形變測量、地質(zhì)調(diào)查、歷史地震和計(jì)算機(jī)模擬數(shù)據(jù)計(jì)算的。計(jì)算機(jī)模擬數(shù)據(jù)在地震概率分析中發(fā)揮了重要作用。周仕勇(2008)改進(jìn)了分析模型，也模擬了中國川西地區(qū)長達(dá)10 000a的理論地震目錄，計(jì)算出了強(qiáng)震在各斷層間的轉(zhuǎn)移概率。

絕大多數(shù)構(gòu)造地震是沿著板塊邊界或者已經(jīng)存在的斷層發(fā)生的，一些地震可能是斷層上黏滑摩擦的不穩(wěn)定性的結(jié)果(Brace et al.，1966)。為了模擬斷層黏滑運(yùn)動，Burridge等(1967)提出了模擬斷層滑動的彈簧滑塊模型，隨后許多學(xué)者進(jìn)行了斷層黏滑現(xiàn)象的模擬研究，模擬結(jié)果可以用來討論和解釋黏滑過程和一些實(shí)驗(yàn)現(xiàn)象。洪漢凈(1994)認(rèn)為斷層摩擦阻力的大小不是決定因素，而具有較長的勻阻段才是發(fā)生大震的首要條件。但是滑塊不是連續(xù)巖石體，滑塊本身不能變形等問題，難于與實(shí)際地質(zhì)構(gòu)造類比。

近年來，細(xì)胞自動機(jī)被許多學(xué)者用來計(jì)算人工合成地震目錄，取得了一些有意義的結(jié)果(Bak etal.，1989;張國民等，1995;劉桂萍等，1995;Liu etal.，1998;Shi et al.，1998;朱守彪等，2006)。但是到目前為止，細(xì)胞自動機(jī)在模擬復(fù)雜的斷層系統(tǒng)，比如刻畫復(fù)雜的斷層結(jié)構(gòu)、復(fù)雜加載等方面，還需要進(jìn)一步的改進(jìn)。

更加接近地質(zhì)過程的有限元分析可能是未來數(shù)值模擬地震活動性的發(fā)展方向，特別是隨著超算系統(tǒng)、并行計(jì)算、多場耦合模擬技術(shù)的完善，復(fù)雜的地質(zhì)、地球物理過程的有限元模擬分析已經(jīng)成為可能。Xing等(2002，2004)利用有限元程序模擬了走滑型斷層系統(tǒng)活動和走滑型地震的孕育、發(fā)生過程，取得了一些有意義的結(jié)果。朱守彪等(2008)模擬了2004年蘇門答臘地震所在俯沖帶上俯沖板片與上伏板塊之間的閉鎖、解鎖、滑動到再閉鎖這一準(zhǔn)周期性過程。朱守彪等(2009)通過模擬2008年汶川MS8.0地震發(fā)生的動力學(xué)過程，計(jì)算了龍門山斷裂帶地震的平均復(fù)發(fā)周期約為3 163a，與其他手段的研究結(jié)果比較一致。薛霆鳩等(2009)利用Ansys有限元軟件模擬了在斷層摩擦系數(shù)非均勻分布時，斷層黏滑運(yùn)動的時空分布特征，分析了地震觸發(fā)、震中遷移、大震重復(fù)和斷層位錯分段特征。

未來的發(fā)展方向應(yīng)該是建立更接近實(shí)際的三維動力學(xué)模型，同時考慮熱、流體等的多場耦合作用，改進(jìn)硬軟件條件，縮短計(jì)算時間步長，獲取更完整的斷層運(yùn)動信息，為建立地震預(yù)測概率模型服務(wù)。

4 結(jié)論與討論

巖石變形進(jìn)入彈性屈服階段以后，系統(tǒng)處于臨界失穩(wěn)狀態(tài)，各質(zhì)點(diǎn)的運(yùn)動將可能偏離正常的軌跡，變得復(fù)雜化。如果以區(qū)域運(yùn)動背景為邊界條件，模擬區(qū)域內(nèi)部的構(gòu)造變形，在正常時期與觀測結(jié)果應(yīng)該是比較接近的，而當(dāng)局部地方處于孕震臨界失穩(wěn)狀態(tài)時，其實(shí)測結(jié)果就有可能與模擬結(jié)果不協(xié)調(diào)，據(jù)此有可能為地震的中期至短期預(yù)測提供有價值的信息。研究結(jié)果顯示，在2006年7月4日文安5.1級地震前，華北地區(qū)地殼運(yùn)動數(shù)值模擬速度與GPS觀測運(yùn)動速度大部分比較吻合，而在個別區(qū)域也存在一定差異，特別是在華北平原塊體的東北部，兩者表現(xiàn)出非常不協(xié)調(diào)的運(yùn)動特點(diǎn)，其原因可能與文安5.1級地震的孕育和發(fā)生有關(guān)。

庫侖應(yīng)力變化觸發(fā)地震是當(dāng)前國內(nèi)外研究的熱點(diǎn)。數(shù)值模擬可在地震預(yù)測研究的地震庫侖應(yīng)力變化分析中發(fā)揮重要作用。在建立三維動力學(xué)模型的基礎(chǔ)上，可模擬分析各種因素對庫侖應(yīng)力變化的大小和分布的影響，獲得更接近實(shí)際的結(jié)果。通過數(shù)值模擬四川紫坪鋪水庫蓄水庫侖應(yīng)力變化對地震活動的影響，發(fā)現(xiàn)隨著水庫蓄水時間的延長，孔隙壓擴(kuò)散對有效附加應(yīng)力場的影響可逐漸擴(kuò)大到地下10km以上的范圍。在幾個典型的區(qū)域有效附加應(yīng)力狀態(tài)和地震活動表現(xiàn)出不同的特征。

美國南加州地震中心把建立區(qū)域地震預(yù)測概率模型作為研究目標(biāo)的工作方式值得中國地震研究工作借鑒。計(jì)算機(jī)模擬合成地震目錄是解決數(shù)據(jù)不足的有效方法。未來的發(fā)展方向應(yīng)該是建立更接近實(shí)際的三維動力學(xué)模型，同時考慮熱、流體等的多場耦合作用，改進(jìn)硬軟件條件，縮短計(jì)算時間步長，獲取更完整的斷層運(yùn)動信息，模擬斷層活動的規(guī)律。

曹建玲，石耀霖，張懷，等.2009.青藏高原GPS位移繞喜馬拉雅東構(gòu)造結(jié)順時針旋轉(zhuǎn)成因的數(shù)值模擬［J］.科學(xué)通報(bào)，54(2):224—234.

CAO Jian-ling，SHIYao-lin，ZHANG Huai，etal.2009.Numerical simulation of GPSobserved clockwise rotation around the eastern Himalayan syntax in the Tibetan Plateau ［J］.Chinese Science Bulletin，54(2):224—234(in Chinese).

陳連旺，張培震，陸遠(yuǎn)忠，等.2008.川滇地區(qū)強(qiáng)震序列庫侖破裂應(yīng)力加卸載效應(yīng)的數(shù)值模擬［J］.地球物理學(xué)報(bào)，51(5):1411—1421.

CHEN Lian-wang，ZHANG Pei-zhen，LU Yuan-zhong，et al.2008.Numerical simulation of loading/unloading effect on Coulomb failure stress among strong earthquakes in Sichuan-Yunnan area ［J］.Chinese J Geophys，51(5):1411—1421(in Chinese).

洪漢凈.1994.從地震模擬看勻阻段與大震的關(guān)系［J］.地震地質(zhì)，16(2):109—114.

HONG Han-jing.1994.Seismic simulation:Correlation of resistance-homogenous fault segment with large earthquakes［J］.Seismology and Geology，16(2):109—114(in Chinese).

嘉世旭，張先康.2005.華北不同構(gòu)造塊體地殼結(jié)構(gòu)及其對比研究［J］.地球物理學(xué)報(bào)，48(3):611—620.

JIA Shi-xu，ZHANG Xian-kang.2005.Crustal structure and comparison of different tectonic blocks in North China［J］.Chinese JGeophys，48(3):611—620(in Chinese).

金文正，湯良杰，楊克明，等.2007.龍門山?jīng)_斷帶構(gòu)造特征研究主要進(jìn)展及存在問題探討［J］.地質(zhì)論評，54(1):37—46.

JINWen-zheng，TANG Liang-jie，YANG Ke-ming，et al.2007.Deformation and zonation of the Longmenshan fold and thrust zone in the western Sichuan Basin［J］.Acta Geologica Sinica，81(8):1072—1080.

劉桂萍，傅征祥.2002.1973年?duì)t霍大地震(MS7.6)最大余震(MS6.3)的庫侖破裂應(yīng)力觸發(fā)［J］.中國地震，18(2):175—182.

LIU Gui-ping，F(xiàn)U Zheng-xiang.2002.The triggering mechanism of the largest aftershock(MS6.3)of the 1973 Luhuo great earthquake(MS7.6)［J］.Earthquake Research in China，18(2):175—182(in Chinese).

劉桂萍，石耀霖，馬麗.1995.地震活動性的細(xì)胞自動機(jī)模型［J］.西北地震學(xué)報(bào)，17(2):20—25.

LIU Gui-ping，SHIYao-lin，MA Li.1995.A cellular automation model of seismicity ［J］.Northwestern Seismological Journal，17(2):20—25(in Chinese).

劉峽，傅容珊，楊國華，等.2006.用GPS資料研究華北地區(qū)形變場和構(gòu)造應(yīng)力場［J］.大地測量與地球動力學(xué)，26(3):33—39.

LIU Xia，F(xiàn)U Rong-shan，YANGGuo-hua，etal.2006.Deformation field and tectonic stress field constrained by GPSobservations in North China［J］.Journal of Geodesy and Geodynamics，26(3):33—39(in Chinese).

冉洪流.2006.時間相依模型與泊松模型計(jì)算結(jié)果的對比分析:以道孚及爐霍潛源為例［J］.震災(zāi)防御技術(shù)，1(3):245—250.

RAN Hong-liu.2006.Comparison of computational results derived from time-dependentand Possionmodels respectively［J］.Technology for Earthquake Disaster Prevention，1(3):245—250(in Chinese).

石耀霖，曹建玲.2010.庫侖應(yīng)力計(jì)算及應(yīng)用過程中若干問題的討滄:以汶川地震為例［J］.地球物理學(xué)報(bào)，53(1):102—110.

SHIYao-lin，CAO Jian-ling.2010.Some aspects in static stress change calculation-Case study on Wenchuan earthquake［J］.Chinese JGeophys，53(1):102—110(in Chinese).

薛霆鳩，傅容珊，陳宇衛(wèi)，等.2009.大尺度斷層活動性數(shù)值模擬及地震學(xué)類比［J］.地球物理學(xué)進(jìn)展，24(5):1616—1626.

XUE Ting-jiu，F(xiàn)U Rong-shan，CHEN Yu-wei.etal.2009.Numerical simulation of large scale fault activity and its seismological analogy［J］.Progress in Geophysics，24(5):1616—1626(in Chinese).

葉金漢，郗綺霞，夏萬人，等.1991.巖石力學(xué)參數(shù)手冊［Z］.北京:水利電力出版社.

YE Jin-han，XIQi-xia，XIAWan-ren，etal.1991.Handbook onmechanical properties of rocks［Z］.ChinaWater Power Press，Beijing(in Chinese).

王凱英.2009.斷-塊模型中走滑型地震應(yīng)力觸發(fā)研究:以青藏高原北部幾次強(qiáng)震為例［J］.地球物理學(xué)報(bào)，52(7):1776—1781.

WANG Kai-ying.2009.Stress triggering of strike-slip earthquakes in fault-blockmodels:A case of several strong earthquakes in Northern Qinghai-Xizang Plateau ［J］.Chinese Journal of Geophysics，52(7):1776—1781(in Chinese).

王瓊，陳學(xué)忠，王海濤.2006.2003年2月24日新疆巴楚-伽師MS6.8地震的靜態(tài)庫侖應(yīng)力觸發(fā)作用［J］.中國科學(xué)院研究生院學(xué)報(bào)，23(2):213—220.

WANG Qing，CHEN Xue-zhong，WANG Hai-tao.2006.Stress triggering of static coulomb stress of Xinjiang Bachu-Jiashi earthquake of MS6.8，F(xiàn)eb.24，2003 ［J］.Journal of the Graduate School of the Chinese Academy of Sciences，23(2):213—220(in Chinese).

萬永革，吳忠良，周公威，等.2002.地震靜態(tài)應(yīng)力觸發(fā)模型的全球檢驗(yàn)［J］.地震學(xué)報(bào)，24(3):302—316.

WAN Yong-ge，WU Zhong-liang，ZHOU Gong-wei，et al.2002.Global testof seismic static stress triggeringmodel［J］.Acta Seismologica Sinica，24(3):302—316(in Chinese).

萬永革，沈正康，刁桂苓，等.2008.利用小震分布和區(qū)域應(yīng)力場確定大震斷層面參數(shù)方法及其在唐山地震序列中的應(yīng)用［J］.地球物理學(xué)報(bào)，51(3):793—804.

WAN Yong-ge，SHEN Zheng-kang，DIAO Gui-ling，et al.2008.An algorithm of fault parameter determination using distribution of small earthquakes and parameters of regional stress field and its application to Tangshan earthquake sequence［J］.Chinese JGeophys，51(3):793—804(in Chinese).

聞學(xué)澤.1990.鮮水河斷裂帶未來三十年內(nèi)地震復(fù)發(fā)的條件概率［J］.中國地震，6(1):8—16.

WEN Xue-ze.1990.Conditional probabilities for the recurrence of earthquakes on the Xianshuihe Fault zone within the coming three decades［J］.Earthquake Research in China，6(1):8—16(in Chinese).

聞學(xué)澤.1998.時間相依的活動斷裂分段地震危險性評估及其問題［J］.科學(xué)通報(bào)，43(14):1457—1466.

WEN Xue-ze.1998.Assessmentof time-dependent seismic hazards on segments of active fault，and its problems［J］.Chinese Science Bulletin，43(23):1937—1950.

張國民，耿魯明，張永仙，等.1995.構(gòu)造塊體的成組孕震模型和前兆場某些特征的分析［J］.地震學(xué)報(bào)，17(1):1—10.

ZHANG Guo-min，GENG Lu-ming，ZHANG Yong-xian，et al.1995.Seismogenic model of earthquakes in groups in tectonic block and analysis for some features ofearthquake precursory field［J］.Acta Seismologica Sinica，17(1):1—10(in Chinese).

張培震，毛鳳英.1996.活動斷裂定量研究與中長期地震危險性概率評價［A］.見:中國地震局地質(zhì)研究所編.活動斷裂研究(5).北京:地震出版社.

ZHANG Pei-zhen，MAO Feng-ying.1996.A quantitative study of active faultand the long-andmid-term seismic risk evaluation［A］.In:Institue of Geology，CSB(ed).Research on Active Faults(5).Seismological Press，Beijing(in Chinese).

張效亮，謝富仁.2009.利用GPS數(shù)據(jù)和實(shí)時概率模型評估川滇南部中長期地震危險性［J］.大地測量與地球動力學(xué)，29(3):38—41.

ZHANG Xiao-liang，XIE Fu-ren.2009.Middle and long-term seismic risk evaluation of southern Sichuan-Yunnan region based on GPS data and real-time probabilistic model［J］.Journal of Geodesy and Geodynamics，29(3):38—41(in Chinese).

周仕勇.2008.川西及鄰近地區(qū)地震活動性模擬和斷層間相互作用研究［J］.地球物理學(xué)報(bào)，51(1):165—174.

ZHOU Shi-Yong.2008.Seismicity simulation in Western Sichuan of China based on the fault interactions and its implication on the estimation of the regional earthquake risk［J］.Chinese Journal of Geophysics，51(1):165—174(in Chinese).

朱守彪，蔡永恩，劉杰，等.2006.利用三維細(xì)胞自動機(jī)模型模擬地震活動性［J］.北京大學(xué)學(xué)報(bào)(自然科學(xué)版)，42(2):206—211.

ZHU Shou-biao，CAIYong-en，LIU Jie，et al.2006.Modeling seismicity by 3-D cellular automata［J］.Acta Scientiarum Naturalium Universitatis Pekinensis，42(2):206—211(in Chinese).

朱守彪，邢會林，謝富仁，等.2008.地震發(fā)生過程的有限單元法模擬:以蘇門答臘俯沖帶上的大地震為例［J］.地球物理學(xué)報(bào)，51(2):460—468.

ZHU Shou-biao，XING Hui-lin，XIE Fu-ren，et al.2008.Simulation of earthquake processes by finite elementmethod:The case ofmegathrust earthquakes on the Sumatra subduction zone［J］.Chinese Journal of Geophysics，51(2):460—468(in Chinese).

朱守彪，張培震.2009.2008年汶川MS8.0地震發(fā)生過程的動力學(xué)機(jī)制研究［J］.地球物理學(xué)報(bào)，52(2):418—427.

ZHU Shou-biao，ZHANG Pei-zhen.2009.A study on the dynamicalmechanisms of the Wenchuan MS8.0 earthquake，2008 ［J］.Chinese JGeophys，52(2):418—427(in Chinese).

Aki K.2002.Synthesis of earthquake science information and its public transfer:A history of the Southern California Earthquake Center［J］.International Handbook of Earthquake and Engineering Seismology，81A:39—49.

Bak P，Tang C.1989.Earthquakes as a self-organized critical phenomenon［J］.JGeophys Res，(94):15635—15637.

Beroza G C，Zoback M D.1993.Mechanism diversity of the Loma Prieta aftershocks and themechanics ofmainshock-aftershock interaction［J］.Science，259:210—213.

Brace WF，Byerlee JD.1966.Stick-slip as amechanism for earthquakes［J］.Science，153(3739):990—992.

Burridge R，Knopoff L.1967.Model and theoretical seismicity［J］.Bull Seismo Soc Am，(57):341—371.

Chinnery M A.1963.The stress changes that accompany strike slip faulting［J］.Bull Seism Soc Amer，53(1):921—932.

Deng J，Sykes L R.1997.Evolution of the stress field in southern California and triggering of moderate-size earthquakes:A 2002 year perspective［J］.JGeophys Res，102:9859—9886.

Hardebeck JL，Nazareth JJ，Hauksson E.1998.The static stress change trigeringmodel:Constraints from two southern California aftershock sequences［J］.JGeophys Res，98:4461—4472.

Harris R A.1998.Introduction to special section:Stress triggers，stress shadows，and implications for seismic hazard［J］.Geophys Res，103:24347—24358.

Heaton T H.1982.Tidal triggering of earthquakes ［J］.Bulletin of the Seismological Society of America，72(6):2181—2200.

Huang Zhongxian，Su Wei，Peng Yanju，e tal.2003.Rayleigh wave tomography of China and adjacent regions［J］.Journal of Geophysical Research，108(B2):2073.

LIU Jie，Vere-Jones D，MA Li，et al.1998.The principle of coupled stress releasemodel and its application ［J］.Acta Seismologica Sinica，11(3):273—281.

Parsons T.2002.Global Omori law decay of triggered earthquakes:Large aftershocks outside the classical aftershock zone［J］.Geophys Res，107.doi:10.1029/200lJB000646.

Reasenberg PA，Simpson RW.1992.Response of regional seismicity to the static stress change produced by the Lama Prieta earthquake［J］.Science，255:1687—1690.

SHIYao-lin，LIU Jie，Vere-Jones D，etal.1998.Application ofmechanicaland statisticalmodels to the study of seismicity and of synthetic earthquakes and the prediction of natural ones［J］.Acta Seismologica Sinica，11(4):421—430.

Stein R S，Barka A A，Dieterich JH.1997.Progressive failure on the North Anatolian Fault since l 939 by earthquake stress triggering［J］.Geophys J Int，128(3):594—604.

Toda S.1998.Stress transferred by the 1995 MW=6.9 Kobe.Japan.shock:Effect on aftershocks and future earthquake probabilities［J］.JGeophys Res，103:24543—24565.

Ward SN.1996.A synthetic seismicity model for southern California:Cycles，probabilities，hazards［J］.JGeophys Res，101:22393—22418.

Ward SN.1997.Dogtails versus rainbows:Synthetic earthquake rupturemodels as an aid in interpreting geological data［J］.Bull Seism Soc Am，87:1422—1441.

Ward SN.2000.San Francisco Bay area earthquake simulations:A step toward a standard physical earthquakemodel［J］.Bull Seism Soc Am，90:370—386.

Xing H，Makinouchi A.2002.Finite element analysis of a sandwich friction experimentmodel of rocks［J］.Pure Appl Geophys，159:1985—2009.

Xing H，Mora P，Makinouchi A.2004.Finite element analysis of fault bend influence on stick-slip instability along an intra-plate fault［J］.Pure Appl Geophys，161:2091—2102.

Wen Xueze.1998.Assessment of time-dependent seismic hazards on segments of active fault，and its problems［J］.Chinese Science Bulletin，43(23):1937—1950.

Working Group on California Earthquake Probabilities.1995.Seismic hazards in Southern California:Probable earthquakes:1994-2024［J］.Bull Seism Soc Am，85:379—439.

Working Group on California Earthquake Probabilities(WGCEP).2003.Earthquake probabilities in the San Francisco Bay region:2002-2031 ［R］.USGeol Surv Open-File Rept，03—214.

Zhang Peizhen，Deng Qidong，Zhang Guomin，etal.2003.Active tectonic blocks and strong earthquakes in the continent of China［J］.Science in China(Ser D)，46(suppl):13—24.

PRELIM INARY STUDY ON APPLICATION OF NUMERICAL SIMULATION METHODS TO EARTHQUAKE PREDICTION RESEARCH(Ⅱ)

DENG Zhi-hui1)HU Meng-qian1)ZHOU Bin1)LU Yuan-zhong2)TAO Jing-ling1)MA Xiao-jing1)JIANG Hui1)LIHong1)
1)Institute of Geology，China Earthquake Administration，Beijing 100029，China
2)Institute of Crustal Dynamics，China Earthquake Administration，Beijing 100085，China

With the advances in simulation techniques and understanding of geodynamic processes，numerical simulation is likely to play an increasingly important role in the research of seismic hazard analysis and earthquake prediction.In this paper，on the basis of the paper“A preliminary study on the application of numerical simulation methods to earthquake prediction research(Ⅰ)”，the possible application of uncoordinated deformation analysis，Coulomb stress changes and earthquake probability modeling to the study of earthquake prediction is further discussed.

When rock deforms from the elastic into the yield stage，the system is in a critical unstable state，the rockmovementmay deviate from the normal track and become complicated.The study results show that，before Wenan earthquake(MS5.1)on July 4，2006，GPS velocity was well consistent with the numerical simulation speed in most areas of North China，while there were some differences in some regions，especially in the northeast of the North China Plain block，where big inconsistency in movement characteristics occurred，resulting perhaps from the preparation ofWenan earthquake.

Research on earthquakes triggered by Coulomb stress change is a focus problem now.Numerical simulationmay play an important role in the analysis of Coulomb stress changes.By constructing threedimensional dynamic model，the effect of various factors on the value and distribution of Coulomb stress change can be simulated，and more realistic results can be obtained.By numerical simulation of Coulomb stress changes to seismic activities beneath Sichuan Zipingpu reservoirs，it is found thatwith the increase of reservoirwater storage time，the pore pressure diffusion in the effective additional stress field will be gradually expanded to the range ofmore than 10km underground.The regional effective additional stress field and seismic activities show different characteristics in several typical regions.

The United States Southern California Earthquake Center has tried to study the earthquake probability as research objectives.It isworthy of referencing in China's earthquake research.Computer simulation of synthetic earthquake catalog is an effective way to solve the lack of data.The future direction of development should be a more realistic three-dimensional dynamic model，taking into account the multi-field coupling between heat，fluid and etc.，improving hardware and software conditions and shortening the calculation time step，obtaining more complete information on fault movement，and simulating the fault activities.

earthquake prediction，numerical simulation，Coulomb stress，earthquake probability modeling，uncoordinated deformation

P315.61

A

0253-4967(2011)03-0670-14

10.3969/j.issn.0253-4967.2011.03.016

2011-06-01收稿，2011-08-22改回。

中國地震局地質(zhì)研究所基本科研業(yè)務(wù)專項(xiàng)(IGCEA1001)和國家自然科學(xué)基金(40841016、40372131、40702056)共同資助。

鄧志輝，男，1962年生，1992年在國家地震局地質(zhì)研究所獲得大地構(gòu)造物理專業(yè)博士學(xué)位，研究員，現(xiàn)主要研究方向?yàn)闃?gòu)造物理、地震預(yù)測方法研究，電話0106-2009089，E-mail:deng 6789@163.com。