余懷忠，程 佳

（中國地震臺網(wǎng)中心，北京 100045）

一種將GPS觀測應(yīng)用于地震中短期預(yù)測的簡單嘗試①

余懷忠，程 佳

（中國地震臺網(wǎng)中心，北京 100045）

本文嘗試采用一種類似于加卸載響應(yīng)比的簡單算法將GPS觀測應(yīng)用于地震中短期預(yù)測。算法根據(jù)孕震區(qū)域在不同階段的加卸載響應(yīng)比變化判斷未來地震的發(fā)生趨勢：當(dāng)?shù)貧ぬ幱诜€(wěn)定狀態(tài)時，得到的時間序列在1.0附近波動；而在臨近地震之前會出現(xiàn)明顯的異常高值。以2004年美國南加州的Parkfield地震為例對方法的有效性進(jìn)行了研究，發(fā)現(xiàn)在地震發(fā)生前半年左右震中附近GPS臺站得到加卸載響應(yīng)比時間序列都發(fā)生了比較明顯的異常變化；而一些距離震中較遠(yuǎn)臺站卻沒有觀察到類似的異常變化。這可能反映了孕震區(qū)地殼介質(zhì)穩(wěn)定性的變化，表明通過處理GPS觀測資料可能獲得新的地震前兆異常。

GPS；加卸載響應(yīng)比；Parkfield地震；地震預(yù)測

Abstract：In this paper we attempt to develop a simple algorithm to apply GPS observations to the short－intermediate－term earthquake prediction practices.By using the Load Unload Response Ratio（LURR）values at different stages，the earthquake potentiality in a seismogenic region is evaluated.The LURR time series fluctuats around 1.0when the source media is in a steady state，and becomes quite high just before a large earthquake.To show the validity of the approach，the 2004Parkfield earthquake in southern California is chosen as the example.We find that about half a year before the earthquake the LURR time series derived from the nearby GPS stations increased obviously；while from some far GPS stations，we can not observe similar anomaly.Therefore，the anomaly may be resulted from the regional crust stability change.The fact suggests that we can get new precursory anomaly to predict large earthquake by properly disposing GPS measurement data.

Key words：GPS；Load unload response ratio；Parkfield earthquake；Earthquake prediction

0 引言

空間測地觀測技術(shù)的出現(xiàn)從根本上突破了傳統(tǒng)大地測量的局限性，為現(xiàn)今地殼運(yùn)動的監(jiān)測開辟了嶄新的途徑［1－10］。其中GPS技術(shù)以其大范圍、高精度、全天候等優(yōu)點(diǎn)為監(jiān)測全球和區(qū)域地殼運(yùn)動帶來了革命性的變化，日益成為三維地殼運(yùn)動觀測的重要技術(shù)手段［11－12］。GPS觀測網(wǎng)絡(luò)實(shí)現(xiàn)了觀測向連續(xù)化、數(shù)據(jù)處理實(shí)時化、技術(shù)手段多元化發(fā)展，使觀測結(jié)果既可把握構(gòu)造運(yùn)動的整體變化，又可勾畫出構(gòu)造運(yùn)動的精細(xì)圖象。與傳統(tǒng)方法相比，不僅觀測效率提高了數(shù)十倍，而且精度也提高了近3個數(shù)量級，使上千公里長的基線觀測精度達(dá)亞厘米量級，這樣的精度足以監(jiān)測大范圍地殼運(yùn)動的微小變化。尤其重要的是，與傳統(tǒng)大地測量方法相比，GPS在不同區(qū)域或不同時期的觀測均可納入全球統(tǒng)一的參考框架，從而使我們能夠定量確定任何相鄰或非相鄰構(gòu)造區(qū)域之間的相互運(yùn)動，從而為各種規(guī)模尺度的地球動力學(xué)的研究提供至關(guān)重要的約束［13－16］。

在我國，自1988年中德合作在滇西試驗(yàn)場建立我國第一個GPS地殼形變監(jiān)測網(wǎng)以來，已經(jīng)建成了以GPS站點(diǎn)為主的“中國地殼運(yùn)動觀測網(wǎng)絡(luò)”，該網(wǎng)絡(luò)包括連續(xù)觀測的25個基準(zhǔn)站，定期復(fù)測的56個基本站和分布于我國多震區(qū)的不定期復(fù)測的1 000個區(qū)域站，初步形成了對我國大陸內(nèi)部較大構(gòu)造塊體運(yùn)動的整體性觀測和對華北、川滇及青藏塊體東北緣等多震區(qū)的區(qū)域性較高空間分辨率的觀測能力［17，4］。建成以來所產(chǎn)出的大量的高精度觀測數(shù)據(jù)，為認(rèn)知中國大陸地殼運(yùn)動特征及其動力學(xué)機(jī)制提供了至關(guān)重要的基礎(chǔ)資料和定量約束［17－20］。目前我國正在實(shí)施中國地殼運(yùn)動觀測網(wǎng)絡(luò)的二期工程——“中國大陸構(gòu)造環(huán)境監(jiān)測網(wǎng)絡(luò)”，將建成大范圍、高密度覆蓋中國大陸的地殼形變觀測網(wǎng)絡(luò)，為研究我國大陸現(xiàn)今地殼運(yùn)動與形變、地球動力學(xué)和地震前兆現(xiàn)象以及探索大陸地震孕育和發(fā)生規(guī)律提供重要基礎(chǔ)。

隨著GPS觀測網(wǎng)絡(luò)的發(fā)展，如何利用地殼形變資料研究地震的孕育過程和規(guī)律，并進(jìn)而進(jìn)行地震預(yù)報研究越來越受到普遍的關(guān)注。地震是地殼介質(zhì)在構(gòu)造力的長期作用下，應(yīng)變能不斷積累，并在一些薄弱部位突然釋放的結(jié)果［21－28］。在地震孕育、發(fā)生過程中，孕震體介質(zhì)必然產(chǎn)生形變，如能采用適當(dāng)?shù)募夹g(shù)手段從GPS地殼形變觀測信息中發(fā)現(xiàn)在不同孕震階段的有特征性的形變信息，無疑可能成為地震預(yù)測的有利判據(jù)。本文嘗試將GPS觀測資料與加卸載響應(yīng)比方法［29－30］相結(jié)合，研究地震發(fā)生的前兆現(xiàn)象。為了檢驗(yàn)方法的有效性，以2004年美國南加州發(fā)生的Parkfield地震為例進(jìn)行了分析。

1 加卸載響應(yīng)比理論

加卸載響應(yīng)比理論是基于非均勻脆性介質(zhì)損傷演化特征提出的一種中短期地震預(yù)測方法［29－30］，研究系統(tǒng)對輸入信號的響應(yīng)，達(dá)到了解該系統(tǒng)的目的。按照定義加卸載響應(yīng)比可以表示為

其中，“＋”表示加載過程，“－”表示卸載過程，X為響應(yīng)率。

假定P和R分別代表一個非線性系統(tǒng)的載荷和響應(yīng)，響應(yīng)率X則可以表示為

ΔR為載荷P發(fā)生微小增量ΔP時響應(yīng)R的變化。

當(dāng)系統(tǒng)處于穩(wěn)定狀態(tài)時，X＋≈X－，Y≈1；當(dāng)系統(tǒng)接近失穩(wěn)破壞時，X＋＞X－，Y＞1（圖1）。圖1顯示當(dāng)系統(tǒng)處于穩(wěn)定狀態(tài)時，加載和卸載響應(yīng)呈線性，當(dāng)系統(tǒng)接近破壞時，響應(yīng)偏離了線性狀態(tài)。因此加卸載響應(yīng)比的演化反映了非線性系統(tǒng)穩(wěn)定狀態(tài)的變化。

圖1 典型的巖石本構(gòu)曲線（P和R分別表示載荷和響應(yīng)）Fig.1 Schematic view of the constitutive law of a brittle mechanic system（Pand Rcorrespond to the load and response of the mechanic system）.

2 GPS觀測中的加卸載響應(yīng)比計算方法

在以往的加卸載響應(yīng)比實(shí)踐時，通常將一定時間窗和空間窗內(nèi)釋放的地震能量被作為響應(yīng)量。加、卸載過程根據(jù)日月潮汐力引起的庫侖破壞應(yīng)力在地震破裂面上的變化ΔCFS進(jìn)行判斷［31－32］，ΔCFS＞0為加載，ΔCFS＜0則為卸載。加卸載響應(yīng)比可以進(jìn)一步定義為加載階段能量釋放和卸載階段能量釋放的之比。為了避免地震數(shù)目太少造成加卸載響應(yīng)比時間序列的強(qiáng)烈波動，計算時間窗內(nèi)通常包含了多個加、卸載循環(huán)過程。

本文所采用的響應(yīng)量不同于傳統(tǒng)的加卸載響應(yīng)比方法，而是直接采GPS觀測的地表形變作為響應(yīng)量。由于GPS觀測時間序列只能按日給出較高精度的形變數(shù)據(jù)，難以滿足加卸載響應(yīng)比方法所需的實(shí)時形變資料（如用以計算潮汐應(yīng)力的加卸載過程），因此本文采用的計算方法延用了加卸載響應(yīng)比的概念，但是實(shí)際操作又有所不同。

假定GPS時間序列觀測的形變總量Di可簡單看作為潮汐力等觸發(fā)因素引起的加卸載形變以及長期構(gòu)造加載產(chǎn)生的形變構(gòu)成，即

其中：Ai和Bi分別表示每天發(fā)生在潮汐應(yīng)力引起的加載和卸載階段的形變總量；ΔC表示長期構(gòu)造應(yīng)力引起的形變。這樣一定時間窗T（T＞1天）內(nèi)的形變總量可以表示為

按照巖石材料的本構(gòu)關(guān)系（圖1），顯然，當(dāng)?shù)貧そ橘|(zhì)處于穩(wěn)定階段時加載階段等同于卸載階段的形變量，即Ai＝Bi；這樣，由式（4）進(jìn)一步得到穩(wěn)定階段不同時間窗內(nèi)的形變總量可表示為

據(jù)此，我們可以定義一個反映地殼介質(zhì)處于穩(wěn)定狀態(tài)的新參量

很明顯當(dāng)系統(tǒng)處于穩(wěn)定階段時，由式（4）和（5）可知Yi＝1，由于長期構(gòu)造加載引起的形變較為穩(wěn)定，這里假定各階段∑TΔC保持相對恒定。而當(dāng)?shù)貧そ橘|(zhì)接近失穩(wěn)破壞時，Ai＞Bi，必然有∑T（Ai－Bi＋ΔC）＞∑TΔC，因而Yi＞1。這樣可以結(jié)合GPS觀測資料，根據(jù)時間序列Yi的變化對地殼介質(zhì)的穩(wěn)定狀況進(jìn)行分析，并對未來地震的發(fā)生趨勢作出判斷。

3 震例分析

為了研究算法的有效性，我們以美國南加州Parkfield地區(qū)2004年9月28日發(fā)生的M6.0地震為例，對震前的加卸載響應(yīng)比時間序列Yi進(jìn)行研究。震中位置如圖2所示。我們收集了震中附近的12個GPS臺站進(jìn)行分析（cand，carh，mnmc，masw，tblp，hunt，pomm，rnch，hogs，mida，lows，land（見圖2）），相關(guān)臺站的GPS時間序列由SOPAC網(wǎng)站下載得到。將各臺站水平位移分量都沿著地震斷層破裂方向投影（走向321°，傾角72°，滑動方向－178°，由Harvard CMT得到），取合成為響應(yīng)量計算加卸載響應(yīng)比。如圖3所示，我們計算得到了這些臺站的加卸載響應(yīng)比時間序列Yi，計算時間窗T＝20天，穩(wěn)定階段形變總量ΔD按照各臺站在穩(wěn)定階段1年內(nèi)的平均形變值計算。

圖2 Parkfield地震震中位置及附近的臺網(wǎng)分布Fig.2 The epicenter location of Parkfield earthquake and nearby GPS stations.

總體來看，采用本算法得到的不同臺站的加卸載響應(yīng)比時間序列Yi在地震發(fā)生之前都發(fā)生了較為明顯地異常變化。加卸載響應(yīng)比值持續(xù)幾年都保持在1.0附近波動（0～2之間），在地震前數(shù)月達(dá)到峰值，而在臨近地震時出現(xiàn)回落。這與加卸載響應(yīng)比方法定義的演化模式相吻合［29－30］。

4 討論與結(jié)論

加卸載響應(yīng)比是依據(jù)潮汐力對地殼的觸發(fā)作用計算的。潮汐力遠(yuǎn)小于構(gòu)造應(yīng)力，因此潮汐力只能觸發(fā)地震，而不能直接導(dǎo)致地震的發(fā)生。當(dāng)構(gòu)造應(yīng)力較低時，潮汐力很難觸發(fā)地震，因此加卸載響應(yīng)比值較小（約為1）；當(dāng)構(gòu)造應(yīng)力較高時，震源區(qū)介質(zhì)對任何微小的擾動都非常敏感，地震就很容易被觸發(fā)，相應(yīng)的加卸載響應(yīng)比值會出現(xiàn)明顯的增大。因此加卸載響應(yīng)比可以看作是反映震源區(qū)地殼介質(zhì)應(yīng)力狀態(tài)的重要標(biāo)志。按照這一理論，震前的應(yīng)力場分布，特別是斷層周圍的應(yīng)力較高區(qū)域?qū)有遁d響應(yīng)時間序列的演化勢必有著更強(qiáng)的調(diào)制作用，其原因在于地震更容易在這些區(qū)域被潮汐力所觸發(fā)，加卸載響應(yīng)比值的異常變化會更為明顯。對一個地震斷層而言，如果其周圍區(qū)域具有較高的應(yīng)力積累，地震就更容易在這些區(qū)域內(nèi)發(fā)生。正因?yàn)榇耍鹎罢鹪磪^(qū)地殼介質(zhì)的應(yīng)力場狀態(tài)及應(yīng)力場分布與加卸載響應(yīng)比的臨界敏感性直接相關(guān)。

圖3 Parkfield地震震中附近臺站得到的加卸載響應(yīng)比時間序列Fig.3 LURR time series derived from the nearby GPS stations before Parkfield earthquake.

從我們的研究結(jié)果來看，很明顯在Parkfield地震發(fā)生前其附近GPS臺站得到的加卸載響應(yīng)比時間序列都發(fā)生了異常明顯的變化。為了進(jìn)行對比，我們對距離震中位置較遠(yuǎn)的一些GPS臺站也進(jìn)行了研究（圖4）。采用同樣的算法和參數(shù)，我們得到了這些臺站的加卸載響應(yīng)比時間序列（圖5）?？梢钥闯鲈诘卣鸢l(fā)生之前這些臺站并沒有出現(xiàn)類似于前述近場臺站的異常變化。這一結(jié)果表明Prakfield周圍臺站觀測到的加卸載響應(yīng)比異常變化帶有明顯的局部化特征，這一局部化特征勢必與隨后的地震發(fā)生密切相關(guān)。這樣震前出現(xiàn)加卸載響應(yīng)比值明顯異常變化的GPS臺站分布可能為地震預(yù)測劃定了一個天然的臨界區(qū)域，地震的震級越高，分布的范圍越大。算法的這一獨(dú)特性質(zhì)使我們有可能使用這一方法針對未來地震進(jìn)行預(yù)測研究。

圖4 部分遠(yuǎn)場的GPS臺站分布Fig.4 Location of some far GPS stations from the epicenter.

本文作為一種將GPS觀測應(yīng)用于中短期地震預(yù)測的簡單嘗試，雖然算法的可行性和有效性還有待各種數(shù)據(jù)的檢驗(yàn)，但是其研究為使用地殼形變資料進(jìn)行地震預(yù)測研究提供了新的參考。算法借鑒了加卸載響應(yīng)比理論的基本概念，延用了輸出信號對輸入信號的響應(yīng)方法，達(dá)到了解地殼介質(zhì)性態(tài)的目的。不同之處就在于響應(yīng)量的選擇：傳統(tǒng)加卸載響應(yīng)比計算的數(shù)據(jù)來源是地震目錄，加載和卸載過程根據(jù)日月潮汐力引起的庫侖破壞應(yīng)力在地震破裂面上的變化進(jìn)行判斷；而我們采用的響應(yīng)量是GPS地殼形變觀測。這一改進(jìn)雖然減少了地下未知因素對加卸載響應(yīng)比計算結(jié)果的影響，但是對GPS的觀測精度以及臺網(wǎng)的分布密度提出了更高要求。我們選擇美國南加州的Parkfield地震進(jìn)行研究，其原因就在于這些區(qū)域具有較高的GPS觀測精度和臺網(wǎng)分布密度。事實(shí)上，算法本身具有良好普適性，隨著國內(nèi)外地殼運(yùn)動觀測網(wǎng)絡(luò)的發(fā)展，將會出現(xiàn)更多的震例。

圖5 遠(yuǎn)場臺站得到的加卸載響應(yīng)比時間序列Fig.5 LURR time series derived from the far GPS stations from the epicenter.

本文研究發(fā)現(xiàn)在美國南加州的Parkfield地震發(fā)生前，其近場GPS臺站的加卸載響應(yīng)比時間序列都發(fā)生了明顯異常變化，這一事實(shí)表明通過適當(dāng)?shù)姆椒軌驅(qū)PS觀測應(yīng)用于地震的中短期預(yù)測研究，極有可能為地震預(yù)報提供一條物理預(yù)報的新途徑。

如果已知發(fā)生異常變化的臺站分布以及加卸載響應(yīng)比的異常程度等，本文提供的算法可以為系統(tǒng)地預(yù)測目標(biāo)地震的震中位置以及發(fā)震的強(qiáng)度提供重要的依據(jù)。即使在GPS網(wǎng)絡(luò)分布較為稀疏的地區(qū)也能通過局部臺站時間序列變化，對各地區(qū)地殼介質(zhì)的穩(wěn)定狀況進(jìn)行分析，為未來地震的發(fā)生趨勢提供參考。

目前算法雖然已經(jīng)有了一些成功的震例研究，然而由于地殼運(yùn)動GPS觀測的密度和精度限制，還難以應(yīng)用到地震預(yù)測研究實(shí)踐，算法本身的發(fā)展也有待未來大量數(shù)據(jù)檢驗(yàn)和證實(shí)。

［1］ Kato T，El－Fiky G S，Oware E N，et al.Crustal Strains in the Japanese Islands as Deduced from Dense GPS Array［J］.Geophys.Res.Lett.，1998，25（18）：3445－3448.

［2］ Holt W E，Chamot－Rooke N，Le Pichon X，et al.Velocity field in Asia inferred from Quaternary fault slip rates and Global Positioning System observations［J］.J.Geophys.Res.，2000，105（B8）：19185－19209.

［3］ Savage J C，W Gan，J L Svarc.Strain accumulation and rotation in the Eastern California Shear Zone［J］.J.Geophys.Res.，（2001，106（B10）：21995－22007.

［4］ 顧國華，申旭輝，王敏，等.中國大陸現(xiàn)今地殼水平運(yùn)動基本特征［J］.地震學(xué)報，2001，23（4）：362－369.

［5］ 黃立人，王敏.中國大陸構(gòu)造塊體的現(xiàn)今活動和變形［J］.地震地質(zhì)，2003，25（1）：23－32.

［6］ 沈正康，王敏，甘衛(wèi)軍，等.中國大陸現(xiàn)今構(gòu)造應(yīng)變率場及其動力學(xué)成因研究［J］.地學(xué)前緣，2003，10（Z1）：93－100.

［7］ 李鐵明，鄧志輝，呂弋培，等.鮮水河斷裂帶北段GPS測量及其運(yùn)動特征［J］.西北地震學(xué)報，2003，25（4）：312－318.

［8］ 陳蜀俊，趙齊樂，曾佐勛，等.利用GPS監(jiān)測和非連續(xù)介質(zhì)力學(xué)研究阿爾金斷裂對中國大陸地殼運(yùn)動變形的影響［J］.西北地震學(xué)報，2005，27（1）：1－7.

［9］ 江在森，楊國華，王敏，等.中國大陸地殼運(yùn)動與強(qiáng)震關(guān)系研究［J］.大地測量與地球動力學(xué)，2006，26（3）：1－9.

［10］ 殷海濤，李杰，張玲，等.基于GPS觀測網(wǎng)的山東地區(qū)地殼運(yùn)動特征分析［J］.西北地震學(xué)報，2008，30（3）：276－281.

［11］ Bennett R A，Wernicke B P，Niemi N A，et al.Contemporary strain rates in the northern basin and range province from GPS data［J］.Tectonics，2003，22（2）：1008，doi：10.1029／2001TC001355.

［12］ Mazzotti S，James T S，Henton J，et al.GPS crustal strain，postglacial rebound，and seismic hazard in eastern North A－merica：The Saint Lawrence valley example［J］.J.Geophys.Res.，2005，110，B11301，doi：10.1029／2004JB003590.

［13］ 江在森，馬宗晉，張希，等.GPS初步結(jié)果揭示的中國大陸水平應(yīng)變場與構(gòu)造變形［J］.地球物理學(xué)報，2003，36（3）：352－ 358.

［14］ 朱文耀，王小亞，符養(yǎng)，等.基于ITRF2000的全球板塊運(yùn)動模型和中國的地殼形變［J］.地球物理學(xué)報，2002，45：197－204.

［15］ 符養(yǎng)，朱文耀，王小亞，等.利用中國地殼運(yùn)動觀測網(wǎng)絡(luò)研究中國大陸相對于ITRF97板塊模型形變［J］.地球物理學(xué)報，2002，45（3）：330－337.

［16］ 李建彪，甘衛(wèi)軍，冉永康，等.青藏高原東部構(gòu)造塊體的運(yùn)動學(xué)及形變特征分析［J］.西北地震學(xué)報，2006，28（2）：97－103.

［17］ 馬宗晉，陳鑫連，葉叔華，等.中國大陸現(xiàn)今地殼運(yùn)動的GPS研究［J］.科學(xué)通報，2001，46（13）：1118－1120.

［18］ Wang Q，Zhang P，F(xiàn)reymueller J，et al.Present－day crustal deformation in China constrained by Global Positioning system measurements［J］.Science，2001，294：574－577.

［19］ Zhang P，Shen Z，Wang M，et al.Continuous deformation of the Tibetan Plateau from Global positioning system data［J］.Geology，2004，32：809－812.

［20］ Shen Z－K，J Lü，M Wang，et al.Contemporary crustal deformation around the southeast borderland of the Tibetan Plateau［J］.J.Geophys.Res.，2005，110，B11409，doi：10.1029／2004JB003421.

［21］ Sornette A，Sornette D.Earthquake rupture as a critical point：Consequences for telluric precursors［J］.Tectonophysics，1990，179：327－334.

［22］ Sornette D，Sammis C G.Complex critical exponents from renormalization group theory of earthquakes：Implications for earthquake predictions［J］.J.Phys.I.France，1995，5：607－619.

［23］ Sammis C G，Smith SW.Seismic cycles and the evolution of stress correlation in cellular automation models of finite fault networks［J］.Pure Appl.Geophys.，1999，155：307－334.

［24］ Bowman D D，Ouillon G，Sammis C G，et al.An observation test of the critical earthquake concept［J］.J.Geophys.Res.，1998，103（24）：359－372.

［25］ Knopoff L，Levshina T，Keilis－Borok V I，et al.Increased long range intermediate magnitude earthquake activity prior to strong earthquakes in California［J］.J.Geophys.Res.，1996，101：5779－5796.

［26］ Rundle J B，Klein W，Turcotte D L，et al.Precursory seismic activity and critical point phenomena［J］.Pure Appl.Geophys.，2000，157：2，165－2，182.

［27］ Jaume S C，Sykes L R.Evolution toward a critical point：A review of accelerating seismic moment／energy release prior to large great earthquakes［J］.Pure Appl.Geophys.，1999，155：279－306.

［28］ 吳忠良.地震震源物理中的臨界現(xiàn)象［M］.北京：地震出版社，2000：1－284.

［29］ Yin Xiang－chu，Chen Xue－zhong，Song Zhi－ping，et al.A New Approach to Earthquake Prediction——The Load／Unload Response Ratio（LURR）Theory［J］.PAGEOPH，1995，145（3／4）：701－715.

［30］ Yin Xiang－chu，Wang Yu－cang，Peng Ke－yin，et al.Development of a New Approach to Earthquake Prediction——Load／unload Response Ratio（LURR）Theory［J］.Pure Appl.Geophys.，2000，157：2365－2383.

［31］ Cochran E S，Vidale J E，Tanaka S.Earthquake tides can trigger shallow thrust fault earthquakes［J］.Science，2004，306：1164－1166.

［32］ Tanaka S，Ohtake M，Sato H.Tidal triggering of earthquakes in Japan related to the regional tectonic stress［J］.Earth Planets Space，2004，56：511－515.

A Simple Attempt to Apply GPS Observation to Short－intermediate－term Earthquake Prediction

YU Huai－zhong，CHENG Jia

（China Earthquake Network Center，Beijing 100045，China）

P315.725

A

1000－0844（2011）01－0009－06

2009－12－14

科技部國際科技合作項目“華北地區(qū)地震構(gòu)造與地震監(jiān)測預(yù)測新方法合作研究”（2010DFB20190）；國家自然科學(xué)基金項目（批準(zhǔn)號：40704007）；中國地震局地質(zhì)研究所基本科研業(yè)務(wù)專項（批準(zhǔn)號：DF－IGCEA－0607－1－16）

余懷忠（1975－），男（漢族），安徽合肥人，主要從事地震預(yù)測理論與前兆方法研究.