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

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

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

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

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

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

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

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

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

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

0 引言

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

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

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

根據(jù)中國活動地塊劃分及活動斷裂分布，結合地殼上地幔的波速結構、地震活動和GPS資料等，確定模型的幾何邊界范圍是 99.8°～121.4°E，27.9°～42.3°N，區(qū)域總面積約為3,840,000km2(圖1)。研究區(qū)域包含的斷裂分布非常密集，為了更加準確地模擬活動斷裂的運動特征，模型中加入了已知的所有活動斷裂。由于斷層空間分布復雜，在建模過程中充分考慮了斷層的走向及延伸的幾何形態(tài)，對其局部的細節(jié)進行一定的簡化。另外，華北地區(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)計算得到的1°×1°中國大陸地區(qū)地殼上地幔波速結構，對數(shù)據(jù)進行篩選，求出平均值，計算得到上述8個區(qū)域分別在上地殼、中地殼、下地殼及上地幔4層上的楊氏模量 (表1)。在模型中將斷層處理為軟弱帶，取地塊介質(zhì)的1/10。地塊的泊松比取0.25，斷層帶的泊松比為0.33。

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

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

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

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

2 庫侖應力變化分析

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

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

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

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

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

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

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

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

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

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

(4)Ⅳ區(qū)位于廣元-大邑斷裂和深部原地斷層交會的區(qū)域。由于該區(qū)距離庫底較遠，在庫體荷載作用下產(chǎn)生的彈性附加應力場本身很微弱，有效附加應力場的變化主要受控于孔隙壓擴散過程，在2008年2—4月期間，ΔCFS上升到0.1MPa左右，該區(qū)相繼發(fā)生了都江堰震群活動，依此判斷該震群活動可能與水庫蓄水有一定的關系。

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

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

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

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

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

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

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

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

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

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

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

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

4 結論與討論

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

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

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

曹建玲，石耀霖，張懷，等.2009.青藏高原GPS位移繞喜馬拉雅東構造結順時針旋轉成因的數(shù)值模擬［J］.科學通報，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).

陳連旺，張培震，陸遠忠，等.2008.川滇地區(qū)強震序列庫侖破裂應力加卸載效應的數(shù)值模擬［J］.地球物理學報，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.從地震模擬看勻阻段與大震的關系［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.華北不同構造塊體地殼結構及其對比研究［J］.地球物理學報，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)_斷帶構造特征研究主要進展及存在問題探討［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年爐霍大地震(MS7.6)最大余震(MS6.3)的庫侖破裂應力觸發(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.地震活動性的細胞自動機模型［J］.西北地震學報，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ū)形變場和構造應力場［J］.大地測量與地球動力學，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］.震災防御技術，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.庫侖應力計算及應用過程中若干問題的討滄:以汶川地震為例［J］.地球物理學報，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ù)值模擬及地震學類比［J］.地球物理學進展，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.巖石力學參數(shù)手冊［Z］.北京:水利電力出版社.

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

王凱英.2009.斷-塊模型中走滑型地震應力觸發(fā)研究:以青藏高原北部幾次強震為例［J］.地球物理學報，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).

王瓊，陳學忠，王海濤.2006.2003年2月24日新疆巴楚-伽師MS6.8地震的靜態(tài)庫侖應力觸發(fā)作用［J］.中國科學院研究生院學報，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)應力觸發(fā)模型的全球檢驗［J］.地震學報，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ū)域應力場確定大震斷層面參數(shù)方法及其在唐山地震序列中的應用［J］.地球物理學報，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).

聞學澤.1990.鮮水河斷裂帶未來三十年內(nèi)地震復發(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).

聞學澤.1998.時間相依的活動斷裂分段地震危險性評估及其問題［J］.科學通報，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.構造塊體的成組孕震模型和前兆場某些特征的分析［J］.地震學報，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ù)和實時概率模型評估川滇南部中長期地震危險性［J］.大地測量與地球動力學，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］.地球物理學報，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.利用三維細胞自動機模型模擬地震活動性［J］.北京大學學報(自然科學版)，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］.地球物理學報，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ā)生過程的動力學機制研究［J］.地球物理學報，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è)務專項(IGCEA1001)和國家自然科學基金(40841016、40372131、40702056)共同資助。

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