新疆中強(qiáng)地震前ECRS方法的異常研究①

唐蘭蘭，王海濤，2，王 瓊

（1.新疆維吾爾自治區(qū)地震局，新疆 烏魯木齊 830011；2.中國(guó)地震局地震預(yù)測(cè)研究所蘭州科技創(chuàng)新基地，甘肅 蘭州 730000）

0 引言

地震學(xué)參數(shù)是地震預(yù)報(bào)研究中普遍使用的分析方法。各種地震學(xué)參數(shù)從不同側(cè)面揭示了地震孕育過(guò)程的某些規(guī)律，并在較大地震前可能出現(xiàn)前兆異常。以往工作中大多是將這些單項(xiàng)異常進(jìn)行經(jīng)驗(yàn)性的、定性的預(yù)報(bào)分析，很多研究人員也嘗試對(duì)單項(xiàng)方法異常特征進(jìn)行定量化描述。如王海濤等［1］利用分布概率法對(duì)地下水氡觀測(cè)曲線(xiàn)進(jìn)行概率化預(yù)處理，進(jìn)而進(jìn)行綜合信息熵的計(jì)算；薄萬(wàn)舉等［2－3］提出多點(diǎn)組斜率法和合成信息流方法，將形變觀測(cè)資料轉(zhuǎn)換成“標(biāo)準(zhǔn)化信息曲線(xiàn)”；鄭熙銘等［4］利用模糊數(shù)學(xué)方法將前兆觀測(cè)資料時(shí)間進(jìn)程變化轉(zhuǎn)化為異常動(dòng)態(tài)從屬函數(shù)；杜學(xué)彬等［5］利用月速率（斜率）法對(duì)地電阻率資料進(jìn)行歸一化處理，生成新的無(wú)量綱數(shù)據(jù)曲線(xiàn)。

通過(guò)將可能與孕震不同階段物理過(guò)程相關(guān)聯(lián)的各種單項(xiàng)方法在不同數(shù)學(xué)或物理模型上組合起來(lái)，有可能可以較為客觀的提取出地震孕育的綜合異常特征，更有效地進(jìn)行地震預(yù)測(cè)。王海濤等［6］提出了地震對(duì)應(yīng)概率譜和累計(jì)滑動(dòng)平均概率方法，并結(jié)合實(shí)際算例，分析了北天山地區(qū)地震學(xué)參數(shù)η值不同考察時(shí)長(zhǎng)的前兆異常時(shí)效屬性，并提出利用多種地震學(xué)參數(shù)進(jìn)行基于地震對(duì)應(yīng)相關(guān)譜的綜合異常研究的方法，即是一種新的探索嘗試。王瓊等［7－8］利用該方法進(jìn)一步研究了新疆天山地區(qū)中強(qiáng)震前多種地震學(xué)參數(shù)綜合異常時(shí)間特征。在該方法深入研究過(guò)程中，王海濤等于2009年將地震對(duì)應(yīng)概率譜重新命名為地震對(duì)應(yīng)相關(guān)譜（Earthquake Corresponding Relevancy Spectrum，（ECRS））［9］方法。

目前常用的地震學(xué)參數(shù)預(yù)報(bào)方法有20多種［10－11］，主要是利用地震目錄上的信息來(lái)描述區(qū)域地震活動(dòng)性的統(tǒng)計(jì)特征，包括時(shí)段區(qū)間、空間范圍、震級(jí)區(qū)間三大要素整體的分布特征。不同區(qū)域的震前活動(dòng)性參數(shù)表現(xiàn)也各不相同。本文采用不同的統(tǒng)計(jì)窗長(zhǎng)和滑動(dòng)步長(zhǎng)進(jìn)行時(shí)空掃描計(jì)算，并參考文獻(xiàn)［12］中高國(guó)英、曲延軍等對(duì)新疆天山地區(qū)地震學(xué)參數(shù)短期前兆異常特征的統(tǒng)計(jì)分析結(jié)果，最終選取了適于短期預(yù)報(bào)且物理意義明確并相對(duì)獨(dú)立的7類(lèi)活動(dòng)性參數(shù)進(jìn)行研究，包括描述地震活動(dòng)的震級(jí)和頻度關(guān)系的b值、η值、缺震；描述地震時(shí)空狀態(tài)的空區(qū)參數(shù)σ值；描述地震活動(dòng)本身強(qiáng)度變化的參數(shù)頻度N值、強(qiáng)度因子Mf值；環(huán)境因子調(diào)制參數(shù)Rm值。在王海濤等［6，9］工作和理論的基礎(chǔ)上，本研究利用上述7種參數(shù)進(jìn)一步分析識(shí)別新疆中強(qiáng)地震前的地震學(xué)參數(shù)綜合異?？臻g演化特征。

1 ECRS方法

（1）值域譜計(jì)算與分析

首先由MAPSIS軟件全時(shí)空掃描程序得到地震學(xué)參數(shù)的逐月空間掃描格點(diǎn)文件，使用本文作者在MATLAB平臺(tái)上研制開(kāi)發(fā)的ECRS1.0軟件，將各參數(shù)的逐月空間格點(diǎn)數(shù)據(jù)文件轉(zhuǎn)換為掃描空間內(nèi)所有格點(diǎn)相對(duì)應(yīng)的參數(shù)值時(shí)間序列文件，時(shí)間序列文件的個(gè)數(shù)即是任一空間格點(diǎn)文件中行與列的乘積。計(jì)算每個(gè)參數(shù)值時(shí)間序列文件的最大值xmax和最小值xmin，從而得到該序列的均值ˉx（i＝1，2，…，n）和標(biāo)準(zhǔn)差σ。

在均值和標(biāo)準(zhǔn)差計(jì)算的基礎(chǔ)上，對(duì)值域分布區(qū)間按一定規(guī)則進(jìn)行定義，統(tǒng)計(jì)觀測(cè)值序列分布在相應(yīng)區(qū)間的頻次，構(gòu)建值域譜曲線(xiàn)。

（2）地震對(duì)應(yīng)相關(guān)譜

通過(guò)上述值域譜的分析與計(jì)算，進(jìn)一步計(jì)算出不同研究考察時(shí)段的分值域地震對(duì)應(yīng)相關(guān)譜。采用地震對(duì)應(yīng)相關(guān)譜可以分析原始參數(shù)值序列基本前兆異常特征，判定不同值域數(shù)據(jù)的異常信度屬性，得到各參數(shù)值序列不同值域的地震對(duì)應(yīng)相關(guān)譜pi。

（3）滑動(dòng)平均相關(guān)度

在地震對(duì)應(yīng)相關(guān)譜計(jì)算的基礎(chǔ)上，將各值域的地震對(duì)應(yīng)相關(guān)譜pi替換相應(yīng)值域內(nèi)的參數(shù)值，每一個(gè)地震對(duì)應(yīng)相關(guān)譜值仍然保留被替換參數(shù)的時(shí)間屬性，再將值域內(nèi)的地震對(duì)應(yīng)譜值按照時(shí)間的先后順序重新組合成一個(gè)新的序列，即相關(guān)譜時(shí)間序列。再根據(jù)不同考察時(shí)段t的長(zhǎng)度（單位：月），選擇多點(diǎn)累加平均、逐點(diǎn)滑動(dòng)計(jì)算方法，最后得到不同地震學(xué)參數(shù)的滑動(dòng)平均相關(guān)度序列ˉPij。

（4）多參數(shù)滑動(dòng)極值平均相關(guān)度

不同參數(shù)在中強(qiáng)地震前表現(xiàn)出不同的前兆異常形態(tài)，使得異常的起始時(shí)間、峰值、結(jié)束時(shí)間存在差異。為了突出中強(qiáng)地震前中短期前兆異常過(guò)程，在單項(xiàng)參數(shù)異常識(shí)別的基礎(chǔ)上，逐點(diǎn)求解不同參數(shù)滑動(dòng)平均相關(guān)度ˉPij序列的多點(diǎn)（本文選取6點(diǎn)作為滑動(dòng)窗長(zhǎng)）滑動(dòng)極值相關(guān)度序列Mij，將不同參數(shù)的滑動(dòng)極值相關(guān)度序列Mj求算術(shù)平均，得到多參數(shù)滑動(dòng)極值平均相關(guān)度Pc。再依據(jù)多參數(shù)滑動(dòng)極值平均相關(guān)度的時(shí)間和坐標(biāo)屬性，將其投影到相應(yīng)的空間位置上，就得到了多參數(shù)滑動(dòng)極值平均相關(guān)度空間格點(diǎn)數(shù)據(jù)。

以上各步驟的具體計(jì)算公式及方法參見(jiàn)文獻(xiàn)［6］～［9］。

2 ECRS方法分析軟件

基于王海濤等［6］提出的地震對(duì)應(yīng)相關(guān)譜方法，在Windows平臺(tái)上利用MATLAB2007語(yǔ)言開(kāi)發(fā)了地震對(duì)應(yīng)相關(guān)譜分析分析軟件——ECRS1.0（圖1）。該軟件具有以下特點(diǎn)：

（1）格式識(shí)別與文件批量讀取功能。ECRS1.0軟件能夠識(shí)別并批量讀取由MAPSIS軟件地震學(xué)參數(shù)空間掃描功能產(chǎn)出的格點(diǎn)數(shù)據(jù)文件，軟件根據(jù)數(shù)據(jù)文件頭信息中的經(jīng)緯度和行列信息及每個(gè)文件名中包含的時(shí)間信息，將數(shù)據(jù)讀入后賦值給相同維數(shù)的三維矩陣變量。

（2）參數(shù)值空間格點(diǎn)文件轉(zhuǎn)換為時(shí)間序列文件。為了構(gòu)建值域譜曲線(xiàn)，需要將空間格點(diǎn)文件轉(zhuǎn)換為時(shí)間序列文件。每個(gè)空間格點(diǎn)文件中的所有格點(diǎn)數(shù)據(jù)都具有一個(gè)相同的時(shí)間屬性和不同的坐標(biāo)屬性，利用其頭文件信息中的經(jīng)緯度信息將每一個(gè)格點(diǎn)數(shù)據(jù)賦予對(duì)應(yīng)的坐標(biāo)屬性，再?gòu)囊幌盗锌臻g格點(diǎn)文件中將具有相同坐標(biāo)屬性的參數(shù)值抽取出來(lái)，然后根據(jù)時(shí)間屬性的先后順序重新組成新的時(shí)間序列文件。時(shí)間序列文件的個(gè)數(shù)是任一空間格點(diǎn)文件中行與列的乘積。

圖1 ECRS程序主界面及參數(shù)配置界面Fig.1 The interface of software and parameter configuration of ECRS.

（3）值域譜曲線(xiàn)的正態(tài)分布檢驗(yàn)。我們?cè)诘玫街涤蜃V曲線(xiàn)后發(fā)現(xiàn)有少數(shù)曲線(xiàn)形態(tài)與標(biāo)準(zhǔn)正態(tài)分布曲線(xiàn)完全不符，說(shuō)明構(gòu)建這條曲線(xiàn)的時(shí)間序列值不符合正態(tài)分布，為了保證后續(xù)的計(jì)算具有統(tǒng)計(jì)學(xué)意義，需將這樣的數(shù)據(jù)在計(jì)算過(guò)程中剔除。結(jié)合MATLAB標(biāo)準(zhǔn)的正態(tài)分布檢驗(yàn)函數(shù)，在程序中設(shè)計(jì)了一種對(duì)稱(chēng)度檢驗(yàn)方法：根據(jù)正態(tài)分布曲線(xiàn)關(guān)于中心對(duì)稱(chēng)的特征首先計(jì)算值域譜曲線(xiàn)的中間點(diǎn)與曲線(xiàn)第一個(gè)點(diǎn)（最小值區(qū)間頻次）連線(xiàn)的斜率及中心點(diǎn)與倒數(shù)第二個(gè)點(diǎn)連線(xiàn)的斜率（由于m值總是偶數(shù)）；之后將兩個(gè)斜率通過(guò)反正切計(jì)算得到兩條直線(xiàn)的方位角，若是符合正態(tài)分布，則兩個(gè)方位角應(yīng)該一個(gè)正值一個(gè)負(fù)值，并且兩個(gè)方位角的絕對(duì)值要非常接近，程序中規(guī)定絕對(duì)值之差小于4°且能通過(guò)matlab的Lilliefors檢驗(yàn)方法（該檢驗(yàn)方法適用于小樣本的正態(tài)分布檢驗(yàn)）lillietest函數(shù)的檢驗(yàn)，即被認(rèn)為是符合正態(tài)分布，能夠繼續(xù)進(jìn)行下一步驟的計(jì)算，不能通過(guò)檢驗(yàn)的序列賦值為0。

（4）地震對(duì)應(yīng)相關(guān)譜時(shí)間序列的空間屬性。若要得到地震對(duì)應(yīng)相關(guān)譜的空間格點(diǎn)數(shù)據(jù)文件，就需要利用每一個(gè)相關(guān)譜值的時(shí)間和空間屬性將地震對(duì)應(yīng)相關(guān)譜重新投影到空間坐標(biāo)上去，構(gòu)建新的地震對(duì)應(yīng)相關(guān)譜空間格點(diǎn)文件。每個(gè)地震對(duì)應(yīng)相關(guān)譜時(shí)間序列對(duì)應(yīng)空間的一個(gè)格點(diǎn)，程序通過(guò)將這些時(shí)間序列文件的尾行寫(xiě)入對(duì)應(yīng)格點(diǎn)的經(jīng)緯度數(shù)據(jù)的方式，給每個(gè)數(shù)據(jù)文件賦予了空間坐標(biāo)屬性，再利用該序列文件的時(shí)間屬性，即可實(shí)現(xiàn)地震對(duì)應(yīng)相關(guān)譜的空間映射。

（5）地震對(duì)應(yīng)相關(guān)譜不同區(qū)域本底值的扣除。從新疆歷史地震活動(dòng)情況來(lái)看，南、北天山地震活動(dòng)性存在較大差異，比如烏?。κ步粎R區(qū)中強(qiáng)震的發(fā)震強(qiáng)度和頻度明顯高于北天山地區(qū)，所以地震對(duì)應(yīng)相關(guān)譜的空間分布值對(duì)于不同區(qū)域存在較大差異，只有扣除不同區(qū)域的本底值才能將其放在同一個(gè)平臺(tái)上分析整個(gè)新疆地區(qū)有無(wú)異常區(qū)存在。如圖2所示，依據(jù)地質(zhì)構(gòu)造及地震活動(dòng)性特征將新疆中強(qiáng)地震活動(dòng)較為頻繁的地區(qū)大致劃分了幾個(gè)區(qū)域，根據(jù)經(jīng)驗(yàn)值將各個(gè)區(qū)域分別扣除本底值（圖上箭頭指示數(shù)值）。一些區(qū)域（如阿勒泰地震帶、阿爾金、東天山）臺(tái)站分布稀疏，使得歷史地震（特別是3級(jí)以下小震）目錄不全，這些區(qū)域空間譜值相對(duì)較低，故未扣除本底值，其結(jié)果僅做為參考。程序通過(guò)預(yù)先得到各個(gè)分區(qū)的邊界信息（經(jīng)度和緯度），將邊界內(nèi)的空間譜值減去相應(yīng)的本底值。

圖2 研究區(qū)目標(biāo)地震震中分布圖（箭頭指示數(shù)字為各個(gè)區(qū)域扣除的背景值）Fig.2 The distribution of target earthquake in the research region（it is background value deducted from various region that the arrow point at number）.

（6）空間色板圖批量繪制。為了更直觀的驗(yàn)證該方法的映震效能，需要將研究時(shí)段內(nèi)逐月的地震對(duì)應(yīng)相關(guān)譜空間格點(diǎn)文件繪制空間色板圖。本軟件利用MATLAB自帶的等高線(xiàn)繪制函數(shù)，同時(shí)建立了新疆的省界及新疆地名數(shù)據(jù)文件，使用循環(huán)的嵌套實(shí)現(xiàn)了圖形的批量繪制及輸出與存儲(chǔ)。文件名以每張圖對(duì)應(yīng)的年月來(lái)命名；使用統(tǒng)一的圖例；標(biāo)注研究區(qū)內(nèi)一個(gè)月后發(fā)生的若干歷史中強(qiáng)地震，同時(shí)標(biāo)注該圖對(duì)應(yīng)的年月。

3 ECRS方法的震例檢驗(yàn)

3.1 資料和參數(shù)選取

選取35°～49°N，73°～96°E范圍作為研究區(qū)域；考慮到地震目錄中地震數(shù)量的完備性，選取研究區(qū)內(nèi)1978年1月－2009年11月的研究時(shí)段內(nèi)MS2以上地震目錄，并對(duì)所選取的地震目錄剔除了余震、震群。對(duì)選定的7項(xiàng)指標(biāo)以1年為時(shí)間窗長(zhǎng)，1個(gè)月為時(shí)間步長(zhǎng)，進(jìn)行2°×2°的全時(shí)空掃描，從而得到研究區(qū)研究時(shí)段內(nèi)選取的7種地震學(xué)參數(shù)的逐月空間格點(diǎn)數(shù)據(jù)文件。

研究中為分析不同時(shí)效、不同窗長(zhǎng)的前兆異常時(shí)空演化特征，分別選取了3個(gè)月、6個(gè)月和12個(gè)月作為研究考察時(shí)長(zhǎng)，以及1°×1°、2°×2°作為研究考察窗長(zhǎng)，通過(guò)調(diào)整不同時(shí)長(zhǎng)、窗長(zhǎng)選定合理時(shí)間和空間尺度。限于篇幅，下文以2°×2°空間窗長(zhǎng)，0.2°經(jīng)緯度滑動(dòng)步長(zhǎng)，12個(gè)月考察時(shí)長(zhǎng)，1個(gè)月時(shí)間步長(zhǎng)的計(jì)算結(jié)果進(jìn)行分析討論。除喀什—烏恰交匯區(qū)預(yù)測(cè)目標(biāo)地震定為該區(qū)域內(nèi)5.5級(jí)以上地震，其他地區(qū)預(yù)測(cè)目標(biāo)地震均定為5級(jí)以上地震。

利用ECRS1.0軟件，基于1978年12月－2007年12月新疆地區(qū)的7個(gè)地震學(xué)參數(shù)的逐月空間格點(diǎn)數(shù)據(jù)文件，計(jì)算得到7個(gè)參數(shù)的地震對(duì)應(yīng)相關(guān)譜數(shù)據(jù)，將此作為單項(xiàng)參數(shù)外推預(yù)測(cè)學(xué)習(xí)庫(kù)數(shù)據(jù)。將2008年1月－2009年11月的7個(gè)地震學(xué)參數(shù)的逐月空間格點(diǎn)數(shù)據(jù)文件進(jìn)行外推預(yù)測(cè)計(jì)算，利用學(xué)習(xí)庫(kù)的地震對(duì)應(yīng)相關(guān)譜進(jìn)行對(duì)比分析，將參數(shù)值替換為相應(yīng)值域區(qū)間的地震對(duì)應(yīng)譜值。為了檢驗(yàn)該方法對(duì)2008年－2009年發(fā)生的中強(qiáng)震的真實(shí)預(yù)測(cè)效果，避免預(yù)測(cè)實(shí)踐中使用目標(biāo)地震發(fā)生以后的資料而逆向影響震前預(yù)測(cè)，因此將2008年1月－2009年11月的實(shí)際資料與地震對(duì)應(yīng)相關(guān)譜學(xué)習(xí)庫(kù)進(jìn)行對(duì)比分析后賦值。在上述基礎(chǔ)上，將1978年12月－2007年12月的空間逐月地震相關(guān)譜數(shù)據(jù)進(jìn)行滑動(dòng)平均概率計(jì)算和多點(diǎn)滑動(dòng)極值平均概率計(jì)算。利用得到的最終數(shù)據(jù)可進(jìn)行歷史中強(qiáng)震的回溯性震例研究及地震危險(xiǎn)性外推式預(yù)測(cè)。考慮文章篇幅限制，本文選取1982－1986年這一時(shí)段做地震對(duì)應(yīng)相關(guān)譜空間特征回溯性震例分析。

3.2 回溯性地震對(duì)應(yīng)相關(guān)譜空間異常特征分析

本文定義有震報(bào)準(zhǔn)率為震前有異常的目標(biāo)地震次數(shù)與目標(biāo)地震總次數(shù)之比。

選取研究區(qū)內(nèi)1982－2007年72次5級(jí)以上地震進(jìn)行回顧性目標(biāo)地震震例研究，通過(guò)對(duì)比分析逐月地震對(duì)應(yīng)相關(guān)譜空間演化圖像，對(duì)地震前地震對(duì)應(yīng)相關(guān)譜空間異常特征進(jìn)行了統(tǒng)計(jì)。統(tǒng)計(jì)結(jié)果表明，72次地震中有62次在震前半年到1年半的時(shí)間開(kāi)始出現(xiàn)較為明顯的高值異常集中過(guò)程，其中39次地震在震后半年到2年的時(shí)間高值區(qū)域逐漸消散或弱化，其余23次在高值異常集中區(qū)或其附近發(fā)震后高值區(qū)域沒(méi)有消散或者弱化，是由于在其后1個(gè)月～1年的時(shí)間里震源區(qū)附近再次有中強(qiáng)震發(fā)生；其余10次地震前沒(méi)有出現(xiàn)高值異常集中過(guò)程，這10次地震主要位于塔里木盆地、西昆侖、阿勒泰等地區(qū)，其原因可能是因?yàn)檫@些地區(qū)臺(tái)站分布稀疏導(dǎo)致歷史地震目錄記錄不全造成。統(tǒng)計(jì)結(jié)果顯示，在1982－2007年這個(gè)統(tǒng)計(jì)時(shí)段內(nèi)的有震報(bào)準(zhǔn)率為83%，地震對(duì)應(yīng)情況較好，漏報(bào)率較低。

圖3 1982－1986年新疆地區(qū)目標(biāo)地震前多參數(shù)滑動(dòng)極值相關(guān)譜空間演化圖像Fig.3 The space image evolvement of multi－parameter sliding extreme－value relevancy spectrum before target earthquake in Xinjiang region from 1982to 1986.

選取1982年－1986年新疆地區(qū)目標(biāo)地震前地震對(duì)應(yīng)相關(guān)譜空間演化圖像進(jìn)行分析。圖3顯示了研究區(qū)內(nèi)7次目標(biāo)地震前的空間圖像，可以明顯的看到震前異常增強(qiáng)、震后異常弱化的演化趨勢(shì)。圖3（a）顯示，1982年8月在天山中段及喀什—烏恰交匯區(qū)分別出現(xiàn)小片異常集中區(qū)域；圖3（b）顯示，1983年1月上述的兩個(gè)異常區(qū)域有集中擴(kuò)展的趨勢(shì)，并在喀什—烏恰交匯區(qū)的高值集中區(qū)的邊緣發(fā)生了2月13日的烏恰6.7級(jí)地震；圖3（c）顯示，1983年2月天山中段的高值集中區(qū)繼續(xù)增強(qiáng)擴(kuò)展，并在其邊緣附近發(fā)生了3月3日的呼圖壁5.0級(jí)地震；圖3（d）和圖3（e）顯示，1983年3月的呼圖壁5.0級(jí)地震后，高值異常集中區(qū)域并未逐漸消散而是仍在擴(kuò)展，9個(gè)月內(nèi)在高值區(qū)的邊緣相繼發(fā)生了6月1日的阜康5.0級(jí)地震和12月15日的烏魯木齊5.0級(jí)地震。從這兩張圖同時(shí)還可以看出1983年2月13日的烏恰6.7級(jí)地震的高值集中區(qū)東部邊緣開(kāi)始弱化而西部邊緣又有所增強(qiáng)，烏恰6.7級(jí)地震之后10個(gè)月又在西部邊境地區(qū)發(fā)生了1983年12月16日的塔吉克斯坦6.1級(jí)地震；圖3（f）和圖3（g）顯示，天山中段的高值區(qū)域在逐漸縮小弱化，而在1985年2－7月喀什—烏恰交匯區(qū)又出現(xiàn)新的異常高值區(qū)在向東逐漸增強(qiáng)擴(kuò)展，并在其東部邊緣發(fā)生了8月23日的疏附7.1級(jí)地震；圖3（h）顯示，1986年3月喀什—烏恰交匯區(qū)和天山中段的高值集中區(qū)域已縮小弱化，這張圖上能夠看出在1986年4月23日富蘊(yùn)5.4級(jí)地震前沒(méi)有出現(xiàn)異常高值集中區(qū)域。

圖3中的7次目標(biāo)地震前有5次在震前震中附近出現(xiàn)較為明顯的高值異常集中區(qū)，它們較為一致的空間異常演化特征是：震前高值區(qū)域在震中附近逐漸增強(qiáng)擴(kuò)展，震后高值區(qū)域逐漸縮小消散或轉(zhuǎn)為低值。

4 認(rèn)識(shí)與討論

本文應(yīng)用基于地震對(duì)應(yīng)相關(guān)譜方法編制的ECRS1.0軟件，回溯了新疆地區(qū)中強(qiáng)震前的多種地震學(xué)參數(shù)的綜合空間異常演化特征，并對(duì)這些綜合異常特征進(jìn)行了識(shí)別和分析，通過(guò)上述分析得到以下幾點(diǎn)認(rèn)識(shí)：

（1）中強(qiáng)地震前1年，震區(qū)及其周?chē)貐^(qū)滑動(dòng)極值平均概率普遍出現(xiàn)高值異常區(qū)，且高值異常區(qū)呈現(xiàn)出由外圍相對(duì)分散演變?yōu)橄蛘鹬袇^(qū)集中擴(kuò)大的時(shí)空變化過(guò)程。

（2）中強(qiáng)地震多數(shù)發(fā)生在地震對(duì)應(yīng)相關(guān)譜高值異常區(qū)內(nèi)部或其邊緣。

（3）中強(qiáng)地震發(fā)生后1年左右，地震對(duì)應(yīng)相關(guān)譜高值異常集中區(qū)趨于分散、弱化或消失。

（4）盡管基于地震對(duì)應(yīng)相關(guān)譜的多參數(shù)滑動(dòng)極值平均相關(guān)度空間異常特征在新疆地區(qū)進(jìn)行的回顧性震例研究的預(yù)測(cè)效能較高，但由于是震后總結(jié)，故異常閾值的判定一定程度受人為調(diào)整的影響，而且不同區(qū)域發(fā)震概率的本底值存在較大差異，所以幾個(gè)區(qū)域選取不同的閾值作為經(jīng)驗(yàn)閾值其合理性仍待檢驗(yàn)。

（5）研究區(qū)多參數(shù)滑動(dòng)極值地震對(duì)應(yīng)相關(guān)譜方法的異常對(duì)應(yīng)的信度較高，但仍存在漏報(bào)現(xiàn)象?？紤]到不同參數(shù)預(yù)測(cè)效能的差異，為提高異常提取的效能，下一步工作將嘗試在對(duì)不同單項(xiàng)參數(shù)的預(yù)測(cè)效能評(píng)估的基礎(chǔ)上，分別賦予各參數(shù)不同的權(quán)，依此計(jì)算多參數(shù)滑動(dòng)極值平均概率。

（6）盡管ECRS方法的研究目前還存在著一定問(wèn)題，但此項(xiàng)研究對(duì)于充分提取地震學(xué)參數(shù)所包含綜合前兆異常信息具有積極的意義。因此，期望通過(guò)地震相關(guān)譜分析ECRS軟件的及時(shí)研制能夠有力推動(dòng)更深入地開(kāi)展這方面的研究。目前ECRS分析軟件在已開(kāi)展的地震預(yù)測(cè)研究中已取得了一些初步的研究成果，但還需要處理更多的不同地區(qū)的地震學(xué)參數(shù)資料，獲得更廣泛的處理應(yīng)用結(jié)果，實(shí)現(xiàn)地震對(duì)應(yīng)相關(guān)譜分析軟件的進(jìn)一步檢驗(yàn)和完善。

［1］王海濤，朱令人，楊馬陵，等.地震前兆觀測(cè)數(shù)據(jù)的信息熵分析［J］.地震，1991，11（5）：13－18.

［2］薄萬(wàn)舉，王彥.單項(xiàng)資料異常變化信息提取的一種新方法［J］.大地形變測(cè)量，1993，9（3）：17－20.

［3］薄萬(wàn)舉，吳翼麟.異常信息流的標(biāo)準(zhǔn)化方法及其應(yīng)用［J］.地殼形變與地震，1993，13（2）：9－15.

［4］鄭熙銘，馮德益.前兆值異常的時(shí)空演化與地震活動(dòng)［J］.地殼形變與地震，1993，13（2）：80－87.

［5］杜學(xué)彬，阮?lèi)?ài)國(guó)，范世宏，等.強(qiáng)震近震巾區(qū)地電阻率變化速率的各向異性［J］.地震學(xué)報(bào)，2001，23（3）：289－297.

［6］王海濤，王瓊.基于地震對(duì)應(yīng)概率譜分析的前兆異常識(shí)別研究［J］.地震研究，2008，31（4）：330－334.

［7］王瓊，王海濤，唐蘭蘭.基于地震對(duì)應(yīng)概率譜的多參數(shù)綜合異常研究［J］.中國(guó)地震，2009，25（3）：294－302.

［8］王瓊，王海濤，唐蘭蘭，等.基于地震對(duì)應(yīng)概率譜的多參數(shù)綜合異?？臻g特征研究［J］.地震，2010，30（1）：98－107.

［9］Haitao Wang，Qiong Wang，Lanlan Tang.Comprehensive precursor anomaly research based on earthquake corresponding relevancy spectrum［J］.Earthq Sci，2009，22（6）：675－681.

［10］高立新，曹井泉.強(qiáng)震前震中區(qū)附近地震學(xué)參數(shù)中期震兆標(biāo)志研究［J］.西北地震學(xué)報(bào)，2004，26（4）：347－352.

［11］中國(guó)地震局監(jiān)測(cè)預(yù)報(bào)司.強(qiáng)地震中短期預(yù)報(bào)方法及其效能研究［M］.北京：地震出版社，2002.

［12］中國(guó)地震局監(jiān)測(cè)預(yù)報(bào)司.西北地區(qū)強(qiáng)地震短期前兆特征和預(yù)測(cè)方法研究［M］.北京：地震出版社，2006.