黃建平 胡 彬

1)中國地震局地震預測研究所,北京 100036

2)中國地震臺網(wǎng)中心,北京 100045

臨界加速地震矩釋放模型的網(wǎng)格計算＊

黃建平1)胡 彬2)

1)中國地震局地震預測研究所,北京 100036

2)中國地震臺網(wǎng)中心,北京 100045

根據(jù)臨界加速地震矩釋放模型的基本原理,根據(jù)研究目的的不同,采取兩種不同的方式進行基于網(wǎng)格的并行處理,即按時間分解和按空間分解。如果要根據(jù)歷史地震目錄資料預測某地區(qū)未來的地震活動,可以采用空間分解的方式;如果要根據(jù)某地區(qū)歷史上特定時間段的地震活動情況來檢驗地震加速破裂模型,則可以采用時間分解的方式。最終,能夠在同樣計算精確度下縮短計算時間或者是同樣計算時間下取得更大的計算精確度。

臨界加速地震矩釋放模型;網(wǎng)格計算;Globus

引言

地震科學作為一門以觀測和實驗為基礎的學科,從19世紀開始,不斷累積了海量觀測數(shù)據(jù)。與此同時,地震學研究中也引入了越來越復雜的算法和模型[1-6]。這樣,人們不得不面臨這樣一個問題:如果以現(xiàn)有的計算方式進行關(guān)于地震的科學計算和地震過程的模擬,就不得不花越來越長的時間和投入更多的精力。但是,網(wǎng)格技術(shù)的出現(xiàn),使人們通過網(wǎng)格技術(shù),可以充分利用各種計算資源,并將這些資源轉(zhuǎn)化為隨處可得的、可靠的、標準的數(shù)據(jù),同時還得到經(jīng)濟的計算和處理能力。這里所謂的計算資源包括網(wǎng)絡通信能力、數(shù)據(jù)資料、儀器設備等各種資源。

目前,國內(nèi)外已有許多大中型的應用網(wǎng)格計算的項目[7]。應用網(wǎng)格技術(shù)使海量的地震觀測數(shù)據(jù)處理和復雜模型運算得以實現(xiàn),從而將使得地震研究方式和方法發(fā)生革命性變革。

2003年6月,在國家發(fā)展和改革委員會的資助下,中國地震局和中國航天二院聯(lián)合構(gòu)建“中國地震減災仿真網(wǎng)格試驗系統(tǒng)CEDASGrid(ChinaEarthquakeDisaster AlleviationandSimulationGrid)”,系統(tǒng)研究網(wǎng)格技術(shù)及其在地震減災領域內(nèi)的應用[8]。本項目的目標是利用現(xiàn)代網(wǎng)格技術(shù)和現(xiàn)代數(shù)字地震觀測技術(shù),通過對網(wǎng)格關(guān)鍵技術(shù)的研究和開發(fā),檢驗網(wǎng)格技術(shù)解決地震問題的技術(shù)可行性,驗證地震仿真網(wǎng)格技術(shù)對提高首都圈地區(qū)觀測資料的分析處理能力、綜合信息顯示能力和分析預報能力,實現(xiàn)對地震震情準實時監(jiān)測和災情仿真的重要作用,為全國的地震仿真網(wǎng)格系統(tǒng)項目的提出和論證提供依據(jù)。

這里介紹的是該項目下的一個地震預報模型,即臨界地震矩加速釋放模型的網(wǎng)格計算的實現(xiàn)。

1 模型的分解

建立網(wǎng)格計算的目標是把分散在各地的大型計算機資源利用網(wǎng)絡和網(wǎng)格中間件連接起來,構(gòu)成一個虛擬的超大型計算機,獲得前所未有的處理能力,從而解決更為復雜的科學計算問題。為了實現(xiàn)上述目標,必須將計算模型按照某種方式分解。常用的分解方式包括:

(1)按步驟分解。某些科學計算由一系列處理序列組成,可以將處理序列分解到網(wǎng)格上不同的計算節(jié)點上,各個計算節(jié)點處理完成后的結(jié)果匯集到計算管理節(jié)點進行最終的合成處理。例如,地震觀測臺站獲得地震相關(guān)數(shù)據(jù),并實時地將數(shù)據(jù)不斷地送到本地分析預報中心進行預處理,然后再送到國家地震局分析預報中心進行處理。

(2)按功能分解。某些科學計算任務能夠分解成相互獨立的不同功能單元。例如,全球氣候模型可以分解為大氣模型和海洋模型,不同的功能單元可以分別在不同的大型計算機上處理,計算節(jié)點之間的數(shù)據(jù)傳輸不是很頻繁。

(3)按數(shù)據(jù)分解。有很多科學計算(仿真)任務的數(shù)據(jù)量很大,可以把不同數(shù)據(jù)塊按其數(shù)據(jù)相關(guān)性的強弱分解到不同的計算機上處理。例如,模擬一個艦隊的戰(zhàn)爭態(tài)勢,雖然理論上艦隊所有軍艦戰(zhàn)斗單位(動力系統(tǒng)、通信系統(tǒng)和武器系統(tǒng)等)之間的行為是相關(guān)的,但是在絕大多數(shù)情況下,相同軍艦戰(zhàn)斗單位之間的行為相關(guān)性比不同軍艦戰(zhàn)斗單位之間的行為相關(guān)性要強。因此,可按其數(shù)據(jù)相關(guān)性的強弱進行分解,每條軍艦在一個網(wǎng)格節(jié)點上仿真。

但并不是所有的科學計算問題都能按照上述方式進行分解。另外,即使能夠被分解,在下述情況下,也不太適合用網(wǎng)格計算:

(1)科學計算問題可以分解成多個步驟,但是每個步驟只處理一次;

(2)科學計算問題需要多個功能組件,但功能組件之間密切關(guān)聯(lián);

(3)科學計算問題需要處理大量的數(shù)據(jù),但數(shù)據(jù)之間的全局相關(guān)性很強。

本文介紹“松散型”科學計算問題,它是指其串行計算中至少存在一個模塊可以分解為多個彼此相對獨立的子模塊,如圖1所示。

圖1 “松散型”科學計算問題特征示意圖

在圖1所示的“松散型”科學計算問題實例中,一個串行模塊C被進一步分解為C1、C2和C3等多個彼此獨立的子模塊。

2 流程設計

下面是地震加速釋放模型的計算流程:

(1)選定一個較大的區(qū)域,并將它網(wǎng)格化。對于網(wǎng)格中每一個節(jié)點,我們都可以認為是未來強震發(fā)生的位置。

(2)根據(jù)地震目錄提供的資料,用timeto-failture公式對累積地震矩進行擬合;

(3)根據(jù)冪函數(shù)與直線擬合誤差的比值來確定震源區(qū)的尺度;

(4)根據(jù)(3),對每一個節(jié)點尋找其最佳震源區(qū)尺度(一般使用圓的半徑表示);

(5)根據(jù)LogR∝0.5M,從而確定未來地震的震級;

(6)多數(shù)震源區(qū)(用圓表示)相重疊的地方往往是未來該區(qū)域發(fā)生強震的地方。

根據(jù)上面的地震加速釋放模型計算流程,按照研究目的的不同,可以采取兩種不同的方式進行基于網(wǎng)格計算的并行處理,即按時間分解和按空間分解。如果要根據(jù)歷史地震目錄資料預測某地區(qū)未來的地震活動,可以采用空間分解的方式,能夠在同樣計算精確度下縮短計算時間或者是在同樣計算時間下取得更大的計算精確度,在用地震加速模型分析全中國的地震活動情況時就使用了空間分解的方式,把全中國分成了6個大區(qū)進行處理。如果要根據(jù)某地區(qū)歷史上特定時間段的地震活動情況來檢驗地震加速破裂模型,則可以采用時間分解的方式,把特定時間段劃分成幾個地震加速期進行處理,同樣可以達到縮短計算時間、提高計算精確度的效果。圖2是基于時間的地震加速釋放計算模型并發(fā)處理的流程圖。

圖2 加速釋放模型并發(fā)流程圖

在上圖中,第5步可以分解為多個沒有數(shù)據(jù)交互的并行計算單元,所以該問題屬于“松散型”科學計算問題,適合在上面設計和開發(fā)的網(wǎng)格計算平臺上運行。

3 南天山地區(qū)的運行實例

我國的南天山地區(qū)自有記錄以來,本地區(qū)地震活動的水平一直處于頻度與震級都很高的狀態(tài)。這為我們研究和檢驗地震加速破裂模型提供了一個寶貴的實例?；诰W(wǎng)格技術(shù)的南天山地震加速釋放模型計算頁面流程如圖3。

4 串行計算與網(wǎng)格計算的效率比較

圖3 計算頁面流程

在1臺PentiumIII1GPC單機上,上面實例采用串行計算的方式進行處理,運行時間大約為6min;在3臺PentiumIII1GPC組成的網(wǎng)格計算系統(tǒng)上,采用Globus網(wǎng)格中間件,對實例進行處理,所有網(wǎng)格節(jié)點都參與了模型的處理,按照模型的分解方法和程序的算法,中間件根據(jù)網(wǎng)格實際運行狀況選擇其中1臺網(wǎng)格計算節(jié)點(資源)作為地震加速釋放計算的預處理節(jié)點,將另兩臺網(wǎng)格計算節(jié)點作為并行計算節(jié)點,Globus中間件管理所有加入網(wǎng)格的計算就資源,整個處理過程對用戶是透明的。上面實例采用網(wǎng)格計算方式處理運行時間大約為3min。加速釋放計算本身越復雜,網(wǎng)格的計算資源就越豐富,網(wǎng)格計算的執(zhí)行效率會體現(xiàn)得越明顯。

(作者電子信箱,黃建平:huangjianping@seis.ac.cn)

[1]Yang W Z,Ma L.Aproposed method for locating the critical region of afuture earthquake using the critical earthquake concept.J.Geophys.Res.,2001,106(B3):4121-4128

[2]楊文政,馬麗.地震活動加速模型及其在中國的應用.地震學報,1999,21(1):32-41

[3]楊文政,邵輝成.用臨界地震概念確定未來地震臨界區(qū)的方法.世界地震譯叢,2002(4):14-23

[4]劉文兵,馬麗.臺灣集集和華北大同地震序列的定量模型分析.地震,2004,24(1):155-162

[5]Bowman DD,Ouillon G,Sammis C G,et al.An observational test of the critical earthquake concept.J.Geophys.Res.,1998,103:24359-24372

[6]Saleur H,Sammis C G,Sornette D.Discrete scale invariance,complex fractal dimensions and log-periodic fluctuations in seismicity.J.Geophys.Res.,1996,101:17661-17677

[7]劉文兵,馬麗,陳會忠,等.地學和地震科學研究中的應用網(wǎng)格計算.國際地震動態(tài),2005(6):9-13

[8]侯建民,劉瑞豐,任裊,等.網(wǎng)格技術(shù)研究進展.地震地磁觀測與研究,2006,27(1):84-89

Grid computation of the CAMR model

Huang Jianping1),Hu Bin2)

1)Institute of Earthquake Science,CEA,Beijing100036,China

2)China Earthquake Networks Center,Beijing100045,China

Based on the CAMR model,two kinds of method are used in grid computation with different purposes,by time decomposition and by space decomposition.The space decomposition would be used if the historic catalogue is used to predict the regional seismicity.And the time decomposition method would betaken to test the CAMR model based on the selected seismicity in certain time interval.Finally,the time in computation could be shortened with same precision or the more precision could be achieved within the same computation time.

CAMR model;grid computation;globus

P315;

A;

10.3969/j.issn.0235-4975.2011.07.008

2010-07-13;

2010-08-21。

中國地震局地震預測研究所基本科研業(yè)務項目0207690236和自然基金項目40574020。