尾波分析比較其他分析方法

㈠ 加卸載響應比理論——一種預測地震及其他地質災害的新理論

尹祥礎

（中國地震局分析預報中心，北京100036）

王裕倉

（中國科學院力學研究所非線性連續介質力學開放實驗室（LNM），北京100080）

摘要由於損傷過程的不可逆性，當孕震區介質受到損傷後，其對載入的響應將不同於卸載響應。載入響應與卸載響應的比Y（稱之為加卸載響應比，英文為Load/Unload Response Ratio簡稱LURR）可以度量孕震區介質的損傷程度或接近失穩的程度，因而可以作為一種地震預測的新途徑。對數百個地震震例的檢驗（震級從4級到8.6級）表明：在主震發生前的一段時間里，Y值顯著大於1。而對於7個穩定區（指在較長時期內未發生過強震，而小震資料又較豐富的地區），在長達20年的時間內，Y值始終在1附近作輕微的起伏。近年來，利用本方法對發生在中國大陸的十幾次中強地震及美國北嶺地震（1994.1.17,M_w=6.7）與日本關東地區地震（1996.09.11,M_s=6.6）作出了成功的中期預測。

關鍵詞地震預報加卸載響應比理論

1引言

地震的物理實質是什麼？從力學的觀點看，它正是震源區介質（岩體）的損傷與快速破壞（失穩）過程，並伴隨應力與應變能的快速釋放。讓我們研究孕震區（含斷層或弱化區的岩體）在高溫高壓下的本構關系，如圖1所示。圖1中縱坐標為廣義載荷P，而橫坐標為對於載荷P的響應R。

我們首先定義如下兩個參數：響應率X與加卸載響應比Y。

響應率X定義為：

第30屆國際地質大會論文集第5卷現代岩石圈運動地震地質

式中：△P表示載荷增量，而△R表示相當於△P的響應增量。

加卸載響應比Y定義為：

第30屆國際地質大會論文集第5卷現代岩石圈運動地震地質

式中：X₊、X_-分別表示在載入及卸載條件下的響應率。

眾所周知，材料若處於彈性階段（圖1中的OA段），載入時（△P＞0）的響應率X﹢，將等於卸載時（△P＜0）的響應率，即Y=1。但是若應力超過彈性後，X₊＞X_-，因而Y＞1。當材料逐步趨向破裂時，Y值也隨之逐漸增大。當趨近於圖1中的頂點T時，X₊趨於無限大；而X_-仍保持為有限值，因而Y值也將趨於無窮。所以頂點T可以作為預測失穩的前兆點。

圖1震源區的本構關系

從損傷力學的觀點看，地震孕育過程就是孕震區介質的損傷過程。因此有希望採用損傷力學中的損傷參數D來定量地刻劃地震的孕育進程。損傷D有多種方法定義。最直接的一種是用彈性模量（4階張量）的變化來定義D。為簡單起見，有時只用彈性模量張量的一個分量來定義D。例如Lamaitre^[23]定義D為：

第30屆國際地質大會論文集第5卷現代岩石圈運動地震地質

式中：E_o為未損傷材料的楊氏模量，E為已損傷材料的楊氏模量。如果卸載時的模量為E_o，則（3）式可表示為：

第30屆國際地質大會論文集第5卷現代岩石圈運動地震地質

式（4）意味著D與Y之間存在著密切的內在關系。也就是說，Y也可以作為孕震區損傷程度的度量。

即使損傷D不按式（3）而另行定義，D與Y之間仍然會存在內在關系。這就明參數Y可以定量地刻劃地震孕育進程，因而可以作為地震預報的定量前兆^{[11～18,24～27]}。

要用加卸載響應比理論進行地震預報，必須首先解決幾個主要問題。一是如何對地球進行載入與卸載，以及如何判別載入與卸載。其次是怎樣選擇合適的參數作為響應。以下分別討論這幾個問題。

（1）如何對地殼塊體進行加卸載？

孕震區的線性尺度可以達到幾百甚至上千千米。對其進行加卸載的方法之一是利用潮汐應力。潮汐力不斷地周期性地變化。也就是說它對地球的各部分不停地進行載入與卸載。

（2）用什麼准則判定載入與卸載？

載入與卸載問題，在塑性力學中有詳細的討論。對於不同的材料應選擇不同的准則。對塑性較好的多種金屬（如低中碳鋼、鋁等）Von Mesis准則比較適用；而對地質材料的破壞，則Coulumb准則^[22]更適合。我們在文獻[12,13,24,26]中對此作了詳細的研究與論述。請參閱上述文獻，本文不再贅述。

（3）選擇什麼參數作為計算Y值的響應

地殼形變、井水水位、地震活動性及其他震源參數以及許多其他地球物理參數都可以作為響應，用於計算Y值。從「八五」期間起，我們與國內外許多地球物理學家合作，開展了多學科的研究，國家科委、國家地震局的「八五」、「九五」攻關項目中，均安排了相應的項目，同時還得到了國家自然基金會及地震科學基金會的支持，取得了比較多的成果^{[1,4～6,10,19]}。在本文中主要介紹以地震能量為響應的加卸載響應YE（在本文中很少涉及其他參數的加卸載響應比，所以仍以Y代替Y_E）。Y定義為：

第30屆國際地質大會論文集第5卷現代岩石圈運動地震地質

式中：E為地震能量，符號「＋」「－」分別表示載入與卸載，N₊、N_-分別表示在規定的時空范圍內發生的載入地震及卸載地震（在載入時段內及卸載時段內發生的地震，在以往的文獻中，有時也簡稱為正地震及負地震）的數目，m為常數，取為0、1/3、1/2、2/3及1。m=1時，Y為所有正地震能量之和與負地震能量之比；m=1/2,E^m為貝尼奧夫應變；m=1/3及2/3時，E^m分別表示震源的線尺度及面尺度；當m=0時，Y=N﹢/N﹣，即正、負地震數1之比。在本文中恆取m=1/2。

取一定的時、空范圍（例如2°×2°，幾個月至1年），按式（5），計算出一個Y值。利用該區域Y值隨時間的變化，可能預測出該區域內未來發生地震的危險程度。

2震例檢驗

我們用國內外數百個已發生的地震，對LURR理論進行了廣泛的震例檢驗，震級范圍從M_s=4到M_s=8.6。檢驗的結果是滿意的^{[13,16,24,26]}。以下是1970～1992年期間發生在中國大陸的10個大地震（M_s≥7）的孕震區LURR隨時間的變化情況（圖2）。

在這段時間內，中國大陸共發生M_s>7級以上地震13個，其中有3個地震（青海地震、通海地震及西藏亦基台錯地震）因數據太少無法利用它們計算LURR之外，對其他10個大震全部進行了研究，地震前各孕震區Y隨時間的變化示於圖2，山圖可知，10個震例中，有9個震例，在震前Y明顯大於1,Y>1的時間大約持續1～3年。

除了較系統地研究了中國大陸的震例（4≤M_s≤8.6）外，根據我們所能得到的資料還研究了日本、美國和其他國家的若干震例。都取得了好的結果^[13,18,26]。

此外，我們選擇了中國大陸的7個區域作對比研究，這7個區域在歷史上曾發生過強震，但近20多年來，地震活動性較低，沒有發生過中強地震，處於低地震活動期。其LURR的變化情況（圖3）與圖2形成強烈的對比，在所有這7個區域里，在20多年（1970～1992）的時間里，其Y值均在1附近作輕微起伏。

圖21970～1992年中國大陸所有7級以上大地震前Y隨時間的變化曲線

a—1973.02.06四川爐霍地震（M_s=7.6）;b—1974.05.11雲南永善地震（M_s=7.1）;c—1974.08.11新疆烏恰地震（M_s=7.3）;d—1975.02.04遼寧海城地震（M_s=7.3）;e—1976.05.29雲南龍陵地震（M_s=7.4）;f—1976.07.28河北唐山地震（M_s=7.8）;g—1976.08.16四川松潘地震（M_s=7.2）；h—1985.08.23新疆烏恰地震（M_s=7.1）;i—1988.11.06雲南瀾滄地震（M_s=7.6）;j—1990.04.26青海共和地震（M_s=7.0）

圖31970～1992年間，7個平靜區的Y-t曲線

a—郯廬斷裂帶南段（35.5°N士1°，118°E±1°）;b—陝北（40.5°N±1°，119°E±1°）;c—川東（31.0°N±1°，118E士1°）；d—魯北（37.0°N±1°，119°E±1°）;e—魯西（37.0°N士1°，118°E±1°）;f—豫北（35.0°N±1°，113°E±1°）;g—魯南（35.0°N士1°，117°E±1°）

3地震預報實踐

近年來我們嘗試用本方法進行地震預報，多次成功對發生在中國大陸的6級以上地震成功地作出了中期預報^{[5,7,15～18]}。此外，還成功地預報了1994年1月17日的美國北嶺地震[5.24]及1996年9月11日的日本千葉地震^[18]。部分被預測的地震震前的Y-t曲線示於圖4。

以下對其中幾個典型地震的預測情況略作說明。

1993年夏初，我們得到由USGS所屬NEIC（美國國家地震信息中心）供給的加利福尼亞州地震目錄，利用此目錄研究了該州的加卸載響應比，發現其中有3個地區在較長的時間（長於一年）內Y值顯著大於1。經研究後，於1993年10月28日寫信給提供我們數據的那位科學家（ISOP項目的負責人），在信中提供了分區的加卸載響應比結果，並且預測：在其中3個區域或其附近，在一年內（1993.10.28～1994.10.28）可能發生中強地震（7>M≥6）。在預測後不到3個月，1994年1月17日發生北嶺地震（圖4e），發生在預測的一個地區的邊緣。再後，1994年9月12日，在另兩個地區的中間發生一個M_s=6.0級地震。

圖4用LURR理論成功地預測的某些中外地震的震前Y-t曲線

a—1991.03.26山西大同地震（M_s=6.1）;b—1993.01.27雲南普洱地震（M_s=6.3）;c—1993.10.26青海共和地震（M_s=6.0）;d—1994.01.03青海共和地震（M_s=6.0）;e—1994.01.17美國北嶺地震（M_w=6.7）;f—1996.09.11日本千葉地震（M_s=6.6）

1996年春，應日本氣象廳科學家的要求對日本的關東地區，和歌山地區及兵庫地區的加卸載響應比進行了分析研究（資料由對方提供），研究後得到幾點結果：①關東地區（按對方提供資料的范圍為了35°～36°N,139°～141°E）在一年內發生M_s=6級地震可能性很大；②和歌山地區在近期內不會發生中強以上地震（對方原來預計此地區危險性很高）；③1995.01.17神戶地震前，加卸載響應比異常很顯著。我們於1996年4月1日將有關結果傳真給了對方科學家，同時寫成論文^[18]於1996年5月間投《中國地震》（季刊）。該文於《中國地震》中文版1996年第三期（1996年9月初出版）及英文版（由美國Allerton出版公司出版）第四期刊出。其後，在1996年9月11日發出了M_s=6.6級千葉地震（35.5°N,140.9°E）。

關於國內的地震預測只討論一個震例——1994年12月31日的廣西北部灣地震（M_s=6.1）。

1993年底我們在「1994年中國大陸地震趨勢研究」的報道中將廣西沿海地區列為地震危險區^[15]。直至於1994年11月分析預報中心召開會商會時，該區未發生過任何中強地震，但當我們研究1992.09.01～1994.08.31時段中國大陸LURR的空間分布時，發現該區域的Y值異常仍非常突出^[16]，因此我們認為該區域仍可能在年內發生強烈地震，結果在1994年最後一天發生了北部灣M_s=6.1級地震，並於1995年1月10日再次發生M_s=6.2級地震。

以上震例是成功的例子，但也有些Y值較高的區域，在預測的時段內並未發生強烈地震。這些區域在一定時段內加卸載響應比升高，說明該區域的地震孕育過程正在進行，但隨後卻可能發生了卸載使地震孕育過程推遲甚至中斷，對於這種情況，如何判別是以後要著力研究的課題。

4加卸載響應比的時空演化特徵

大量的震例研究表明，LURR的空間分布圖像是很復雜的。當一個地區未來要發生強烈地震前，將出現一系列高Y值區，這些高Y值區往往連成一個環狀，形成麵包圈圖像^[7]，大部分未來的地震將會在發生麵包圈內或其鄰域。圖5是1979.03.14雲南普洱M_s=7.0級地震震前一年間該區域內Y值的空間分布。由圖可見，高Y值區圍繞著未來的震中形成一個面積約為5°×5°的麵包圈。將LURR的空間分布作成立體圖，每一個高Y值區形成一個高峰，很多個高Y值區就形成群峰突起的圖像。形成鮮明對比的是，在地震活動性低（指未發生強震）的區域內，Y值在1°上下輕微起伏，所以LURR的空間分布立體圖就像平原地區的地形，我們形象地稱這種圖像為「一馬平川」。

圖51979.03.14雲南普洱M_s=7.0級地震震前一年時段內，Y=2.0的等值線圖圖中符號代表未來的震中

對於同一地區，在地震孕育過程中，不同時段的LURR空間分布圖像是不同的，也就是說，空間分布圖像隨時間發生變化。我們發現一個非常有趣的現象：即強震前多個高Y區向未來的震中遷移，稱之為高Y值區的會聚現象^[8]。

研究了1970年後中國大陸的12個M_s≥7.0大震^[8]。12個震例中，有11個發現了會聚現象，且未來的震中處於麵包圈內，只有1992.04.23中緬邊界上的M_s=7.0級地震，未來的震中處於麵包圈外，距圈的外邊界約50km。這可能與該地震發生在兩國邊境地區，緬甸一側的數據不好收集有關。

進一步，我們還發現：強震發生前，高Y值區遷移速度在同一地區，近似不變，大致為100km/a的量級。但不同地區的遷移速度有所差異^[8]。

Scholz曾經撰文說，華北地區的形變鋒（deformation front）傳播速度約為150km/a,Press和Allen則觀測到美國南加州地區的形變波（deformotion wave）速度為100km/a。這三者在物理上是有關系的，而且其數值在數量級上也是彼此相符的。

5展望

前已述及，除了地震能量外，其他許多有關的地球物理參數（如地下水位、地殼體應變、地傾斜、地磁參數、尾波Q…）均可選擇為響應，進行加卸載響應比的研究^[1,2,4,10]。圖6為取尾波Q值的例數作為響應，計算LURR—Y_Q的例子。圖6所示為美國加利福尼亞州南部北嶺地震（1994.01.17,M_w=6.6）前該區Y_Q隨時間的變化情況，將它與圖4e作比較後可以看出，二者在定性上是一致的。

圖6美國加利福尼亞州南部北嶺地震的Y_Q變化圖

對同一時空域用眾多的參數可以計算出眾多的LURR值，然後進行綜合預報，必然會提高用加卸載響應比理論進行地震預報的精度。

簡而言之，加卸載響應比理論可能為地震預報開辟了一條新途徑。現在國內地震界有不少人在研究，應用與改善它^[5]。

近來通過研究，表明北京地區的有感地震（指6＞M≥4）前^[9]甚至北京地區的礦震（M＞2）前，Y值也有較明顯的升高。這說明，除天然地震外，對於誘發地震（礦震、水庫地震^[4]…）以及某些其他地震災害（如岩爆、滑坡、火山噴發…），也可能用LURR理論進行預測。

致謝謹向傅承義、Keiti Aki、秦馨菱、王仁、陳章立、何永年、葛治州、陳鑫連、梅世蓉、羅灼禮、張國民、李宣瑚、張伯民教授及E.A.Bergman、K.Hosono、H.P.Ouyang博士致以誠摯的謝意，感謝他們多方面的支持與幫助。

本項目得到國家自然基金會（批准號19732006）、國家科委及國家地震局「八五」及「九五」攻關項目、地震學基金會以及中國科學院LNM開放實驗室（Lab of Nonlinear Mechanics of Continuous Media,Institute of Mechanics）的資助。

參考文獻

[1]陳建民，張昭棟，楊林章，石榮會，張繼紅.地下水位固體潮響應比的地震異常.地震，1994,10（1）:73～78.

[2]陳學忠，尹祥礎.水庫地震主震前加卸載響應比的變化特徵.中國地震，1994,10（4）:361～367.

[3]陳學忠，尹祥礎，K Ake,H Ouyang，宋志平，王裕倉.以介質參數波Q1作為響應的加卸載響應比研究.中國地震，1996,12（3）:243～249.

[4]國家地震局分析預報中心.關於1994年我國地震趨勢的研究報告.見：國家地震局分析預報中心編.中國地震趨勢預測研究（1994年度）.北京：地震出版社，1994,43.

[5]李宜瑚.「八五」地震預報理論及方法攻關新進展之一：加卸載響應比理論預測洛杉磯地震獲得成功.國際地震動態，1994,10（4）:24～25.

[6]劉桂萍，馬麗，尹祥礎.首都圈地區中等地震前響應比特徵的研究.地震，1994,（6）:34～39.

[7]宋志平，尹祥礎，陳學忠.加卸載響應比的時空演變特徵及其對地震三要素的預測意義.地震學報，1996,18（2）:179～186.

[8]施行覺，許和明，萬永中，盧振剛，陳學忠.模擬引潮力作用下的岩石破裂特徵：加卸載響應比理論的實驗研究之一。地球物理學報，1994,37（5）:633～637.

[9]王丹文.加卸載響應比理論在以磁報震中的應用探索.地震地磁觀測與研究，1995,（16）:26～29.

[10]楊林章，何世海，郗欽文，黎凱武，李松陽.用體應變潮汐加卸載響應比確定岩石彈性的變化.中國地震，1994（增刊）：90～94.

[11]尹祥礎.地震預測新途徑的探索.中國地震，1995,3（1）:1～7.

[12]尹祥礎，尹燦.非線性系統的失穩前兆與地震預報.中國科學，1991,（5）:512～518.

[13]尹祥礎，陳學忠，宋治平，尹燦.加卸載響應比理論（LURR）：一種新的地震預報方法.地球物理學報，1994.37（6）:767～775.

[14]尹祥礎，陳學忠.加卸載響應比理論及其地震預測中的應用研究進展.地球物理學報，1994,37（增刊）：223.

[15]尹祥礎，陳學忠，宋志平.加卸載響應比理論的新進展及其在地震趨勢研究中的應用.見：國家地震局分析預報中心編.中國地震趨勢預測研究（1994年度）.北京：地震出版社.1993.

[16]尹祥礎，陳學忠，宋志平.從加卸載響應比的時空分布，研究中國大陸未來的地震大形勢.國家地震局分析預報中心編：中國地震趨勢預測研究（1995年）.北京：地震出版社，1994.

[17]尹祥礎，陳學忠，宋志平，王裕倉.由加卸載響應比的時空變化預測中國大陸地震趨勢.見：國家地震局分析預報中心編.中國地震趨勢預測研究（1996年）.北京：地震出版社.1995,75～178.

[18]尹祥礎，陳學忠，宋治平，王裕倉.關東地區載入響應比的時間變化及其預測意義.中國地震，1996,12（3）:331～335.

[19]張繼紅.響應比法的地磁異常分析.地震地磁觀測與研究，1995,16:61～63.

[20]B.K.Atkinson.The Theory of Subcritical Crack Growth with Application to Minerals and Rocks.In:Fracture Mechamics of Rock（Ed.B.K.Atkinson）,Academic Press,London,1987.

[21]Ding Zhongyi et al..Seismic triggering effect of tidal stress.Tectonophysics,1994,93:319～335.

[22]J.C.Jaeger,N.G.W.Cook.Fundamentals of Rock Mechanics.Chapman and Hall,London,1976,78～99.

[23]J.Lemaitre.How to Use Damage Mechanics.Nuclear Eng.＆.Design,1984,80:233～245.

[24]Yin Xiangchu,Chen Xuezhong,Song Zhiping and Yin Can.The Load-Unload Response Ratio Theory and its Application to Earthquake Prediction.Journal of Earthquake Prediction Research,1994,3:325～333.

[25]Yin Xiangchu,Yin Can and Chen Xuezhong.The Precursor of Instability for Nonlinear Systen and Its Application to Earthquake Prediction——the load-Unload Response Ratio Theory.Non-linear dynamics and predictability of geophysical phenomena（W.I.Newmana,A.M.Gabrelov and D.L.Turcote eds.）.AGU Geophysical Monograph,1994,83:55～60.

[26]Yin Xiangchu,Chen Xuezhong,Song Zhiping and Yin Can.A New Approach to Earthquake Prediction——The Load/Unload Response Ratio（LURR）Theory.PAGEOPH,1995,145（3/4）:701～715.

[27]Иин Ксянчу（Yin Xiangchu）.Новыйнодходк Ирогногу Землетрясеннй.НРИРОД А，1993,（1）:21～27.

㈡ 五種預報地震的方法

地震學預測方法

一、前兆性地震活動：(1)地震活動增強與平靜；(2)震群活動；(3)前震活動；(4)誘發地震；(5)地震窗；(6)地震震中遷移；(7)相關地震；(8)地震活動重復性。

二、地震空區和條帶：(1)地震空區；(2)地震條帶。

三、強震序列類型和判定：(1)地震序列的類型；(2)序列類型判定。

四、地震活動性參數的應用：(1)N值（頻度）；(2)b值；(3)震情指數（A(b)值）；(4)缺震；(5)應變釋放（能量釋放）；(6)地震集
中度；(7)震群的U、K、ρ值；(8)地震活動熵；(9)地震活動度S；(10)AC（Cn）值；(11)分維數。

五、震源機制、介質參數的應用：(1)震源機制；(2)波速比；(3)振幅比As/Ap；(4)介質品質因子Q值；(5)尾波振動持續時間比和尾波品質因子Qc值；(6)震源參數。

六、小震調製作用/響應比在短臨預測中的作用：(1)小震調制比；(2)加卸載響應比。

七、震級、時間、地點預測的主要地震學方法的歸納：(1)空間圖像方法；(2)時間進程方法；(3)地震序列方法；(4)地震相關方法；(5)震源及介質參數方法；(6)震兆合成方法。

八、震後趨勢判定主要方法：(1)h值；(2)b值截距法；(3)「密集—平靜」現象；(4)等待時間法；(5)應變釋放曲線。

地殼形變學方法

一、定點傾斜、應變預報地震方法：(1)物理參數計算分析方法；(2)濾波分析方法；(3)幾何分析方法；(4)綜合分 析方法。
二、大面積形變測量與跨斷層測量預報地震方法：(1)精密水準復測；(2)水平形變觀測；(3)跨斷層場地觀測方法。 三、重力測量預報方法：強地震一般發生在重力變化梯級帶上或等值線轉折部位，並大多位於低值一邊。
四、
GPS觀測預報地震方法：強地震前可以觀測到GPS的異常變化特徵。一些觀測結果表明，在高剪應變集中及主應變增大區，在其周圍及中心區，往往是潛在強地震的震源區域。大震前地表會大面積地出現向某一方向運動的現象並在震中及鄰近地區產生高剪應變集中區。
地磁地電方法
一、地磁學方法：(1)地磁轉換函數法；(2)地磁加、卸載相應比法；(3)空間相關和加權差分法；(4)低點位移法。
二、地電學方法：(1)常規預報地震方法；(2)地震地點學其他方法(自然電場法、地磁輻射法、大地電場法、大地電磁測探方法、大氣電位方法)。
三、遙感觀測技術的應用：(1)地震前後斷裂帶內外的衛星紅外溫度差異動態過程；(2)地震前後潛熱通量異常； (3)地面溫度的反演；(4)長波輻射信息的短周期地震指標。
地下流體方法

通過分析觀測數據，提取異常，預測地震。

一、地下流體物理類：(1)井水位；(2)井水溫(地熱)；(3)水電導率。

二、地下流體化學類：(1)離子組分；(2)微量組分和放射性組分；(3)氣體組分。

三、地下流體宏觀現象：水色、水氣、氣味、溫度、水位、流量等。

綜合分析方法

一、綜合概率法：(1)概率合成方法；(2)概率增益方法；(3)多元統計組合模型的應用；(4)可拓綜合模型前兆預測 研究。

二、模式識別：(1)基於地震前兆的模糊地震預報；(2)模糊聚類分析；(3)模糊信息檢索方法；(4)模糊綜合評判； (5)CORA—3方法。

三、「故障診斷技術」的應用。

四、專家預測系統。

五、遺傳—神經網路模型：(1)遺傳演算法；(2)神經網路模型。

群測群防常用方法

一、動物宏觀異常觀測方法

二、地下水宏觀異常觀測方法

三、土地電

四、土地磁儀

五、地下水群測網