基于數理分析法進行的云南地區地震分析與預測

馬江濤，周維云，晁智華

（大理學院天文與科技史研究所，云南大理 671003）

基于數理分析法進行的云南地區地震分析與預測

馬江濤，周維云，晁智華

（大理學院天文與科技史研究所，云南大理 671003）

結合基于數理分析與可公度法得到的云南地區地震分析結果來預測云南地區地震，結果顯示用兩種方法分析得出的結果大致相同，這種將不同方法結合起來研究地震的方法將為地震的研究提供途徑和方法，使預測更加精準。

云南；地震；分析預測

太陽作為地球能量的源泉，它的變化直接或間接地影響和調制著地球各部分和整體的運動與變化〔1〕。隨著對太陽活動與地震的深入研究，一些結果表明：太陽黑子活動與地震有很大的聯系〔2-4〕。太陽黑子活動在地殼內產生的感應電流、地磁異常，以及通過對大氣環流等的作用影響地球自轉，都可能成為觸發和影響大地震的因素。雖然這些結論還不夠明確，并且不同方法研究得出的結果存在一定的差異，但它們都表明太陽活動對地震是有影響的〔5〕。以下將使用數理分析法分析云南地區的地震，并用可公度法對分析得出的結果〔6〕進行對比，對云南地區的地震進行綜合分析與預測。

1 用數理分析法對數據進行分析

1.1 數據的來源與選取 人類地震科學技術發展近百年的歷史和我國地震監測預報科學技術發展近百年的歷史均表明：準確掌握歷史地震資料是對地震統計研究的關鍵所在〔7〕。本文收集了2001年7月中國地震局監制預測司編制的《中國地震活動分布統計圖籍》（圖紙版）的1900年-1999年21世紀云南省及其周邊地區5級以上的地震數據和張元東著中國華僑出版社的《太陽黑子》（網絡版）的太陽黑子年相對數進行分析研究。選取MS≥5.0的地震，是為了避免余震帶來的影響；樣本所占時間要長，所取樣本的時間段上要盡量多覆蓋幾個太陽活動周期；另外選擇地震多發地區，這樣的地區從構造上說容易發生地震，以便將太陽活動作為一個觸發因素來探討其間的關系；最后所選擇地域應相同或相似，以使發震構造相同或相似，盡可能排除不同構造類型地震的干擾。研究數據的選取范圍如表1所示。

表1 數據選取范圍

1.2 頻數統計分析

圖1 1900年-1999年云南省及其周邊地區5、6、7級地震和太陽黑子活動分布統計對比圖

從圖1的對比分析中發現5級、6級地震在1930年和1976年附近發震頻數比較高，同時這兩年也是發震的高峰年，從太陽黑子活動來看，這兩年是太陽黑子活動的下降年（1930峰年后3年）和谷年（1976年）；而7級以上的地震云南省共發生12次，主要集中在谷年和峰年的前后3年。

1.3 位相統計分析 取太陽活動峰年為零年（即M年），以此年為準，往前各年分別為-1、-2、-3、-4、-5年，往后各年分別為+l、+2、+3、+4、+5年，統計在這些年內發生M≥5.0級地震的次數。

圖2 云南省及其周邊地區5、6、7級地震頻數位相曲線圖

從圖2中可以看出：5、6、7級地震在太陽活動峰年前后3年達最大值，且有較好的對稱關系，在峰年和前后1年地震較少，主要集中在峰年的前后3年。所以，云南省及其周邊地區在太陽黑子活動的上升和下降年都會發生地震，且在太陽活動峰年的前后3年發生較多，峰年不發生較大地震。

1.4 地震波能量統計分析 對地震數據和太陽黑子相對數進行了頻數統計分析，下面就地震活動性中另一個重要指標──地震波能量進行分析。震級的測定：地震強度的估計，最初以極震區最大烈度為依據，但烈度受外界的影響很大，不是一個很穩定的物理量。微觀地震學發展后，人們已知道地震釋放大量能量，并認為以地震能量規定地震大小較為合適。若以E為地震能量，M為地震大小的等級，則可用下式表示：

1956年古登堡和李克特根據大量觀測資料確定了底值和系數為

現在通用的公式〔8〕：

圖3 1900年-1999年地震波能量和太陽黑子活動統計分布對比圖

從圖3中可以直觀看出20世紀后50年里太陽黑子活動比較強，地震能量釋放也比較多。在太陽黑子活動最強的第19個太陽黑子活動周和太陽黑子活動最弱的第14個太陽黑子活動周內地震能量釋放較小（幾乎沒有釋放能量）。進一步發現地震釋放能量大多出現在太陽黑子活動的下降年，有兩次在谷年，兩次在上升年。在波峰年沒有發生，多出現在波峰的前后3年。

1.5 基于數理分析法得出地震結論與預測 通過以上分析，可以得出以下結論并作出預測：①太陽黑子活動對云南省及其周邊地區的地震有影響，云南省及其周邊地區的地震與太陽黑子活動有很好的相關性。太陽黑子活動峰年前后3、4、5年內地震發生頻繁。7級以上地震多發生在谷年附近（上升年和下降年都有），從未發生在峰年。地震所釋放的能量在太陽活動的下降年較多，20世紀后50年釋放能量相對較多；②太陽黑子活動的變化均有明顯的11年周期，但11年活動周期僅僅是其規律性表現的一個平均周期，而地震總能量和頻數的周期性變化比較復雜，并非簡單的l1年周期；根據太陽黑子活動的周期性〔9〕，2013年會是一個太陽黑子活動的峰年，經過上面的三個分析再綜合預測，峰年的前后3年是地震發生的頻繁年，所以2010年、2016年可能是地震發生的頻繁年。

2 可公度法分析得到的結果及與本文預測結果比較

這里引用《基于可公度方法的云南地區地震數據分析與預測》（譚鍇）中的結果。文中選用的數據是云南地區（22°～29°N，98°～104°E）1800年到2007年發生且震級在6.7級以上的地震的地震大小和發生的時間。數據處理方法如下所示，見表2～表5。

表2 三元可公度法得出的結果

從表2中取m值較大的年份，則發生地震年份可能性2013＞2009＞2010=2016。

表3 四元可公度法得出的結果

從表3取m值較大的年份，則發生地震年份可能性2016＞2010＞2013＞2015。

表4 五元可公度法得出的結果

從表4取m值較大的年份，則發生地震年份可能性2010＞2011＞2016＞2012。

由于有些年份出現次數較少，所以為虛報，如2009年、2011年、2012年、2015年。所以預測結果應該是2010年、2013年、2016年。

表5 兩種方法得出的預測結果

從表5中發現兩種方法預測到的只有一年2013年不一致，其余兩年2010年和2016都一致，所以2010年和2016年可能是云南地區地震發生的頻繁年。

3 討論

地震的發生是多因素的〔10〕，掌握地震的規律和確定地震的發生地點仍然有許多科學問題需要進行長期深入的研究。地球處于一個時刻變化的環境中，單從某一個方面進行分析都有局限性，應該將多種方法研究的結果進行綜合的分析和比較，從而對地震的研究提供新的途徑和方法，使地震的預測更加準確。

〔1〕王家龍，孫靜蘭.太陽活動及其對地球環境的影響〔J〕.第四紀研究，2002，6（4）：81-84.

〔2〕樊忠玉.中國地震與太陽活動相關的初步分析〔J〕.紫金山天文臺臺刊，1997，16（3）：151-154.

〔3〕張兆鑫，杜升云.云南省地震與太陽活動的統計關系研究〔J〕.云南地理環境研究，2008，20（4）：118-120.

〔4〕方煒，劉春，張春生，等.太陽黑子與全球強震活動〔J〕.高于地震，2003，15（4）：27-31.

〔5〕蔣伯琴.太陽黑子、磁暴與地震活動的關系〔J〕.地震學報，1985，7（4）：452-460.

〔6〕譚鍇.基于可公度方法的云南地區地震數據分析與預測〔J〕.地震地磁觀測與研究，2009，30（6）：12-16.

〔7〕中國地震局監測預報司.中國地震活動分布統計圖集〔M〕.長沙：湖南地圖出版社，2001：257-258.

〔8〕郭履燦，趙風竹.震級與震源參數測定〔M〕.北京：中國科學技術出版社，1986：177-179.

〔9〕占臘生，葉藝林，袁文亮，等.太陽活動周期及其數學描述〔J〕.天文研究與技術，2008，2（15）：106-121.

〔10〕章海峰.揭秘大地震〔J〕.飛碟探索，2010（5）：52-53.

Analysis and Prediction of Earthquake in Yunnan Based on Mathematical Analysis

MA Jiangtao,ZHOU Weiyun,CHAO Zhihua
（Astronomical and Technology History Institute,Dali University,Dali,Yunnan 671003,China）

Mathematical analysis and the method of commensurability are adopted to predict earthquake in Yunnan,which shows that the results of the earthquake analyzed with two different methods are much the same.Thus combining the two different methods to study earthquake will provide new approach to the earthquake research,which makes the prediction more accurate.

Yunnan;earthquake;analysis and prediction

P65［文獻標志碼］A［文章編號］1672-2345（2012）04-0039-03

2011-05-17

2011-08-24

馬江濤，講師，主要從事控制理論與控制工程研究.

（責任編輯 袁 霞）