關(guān)玉梅，王紫燕，楊錦玲，秦雙龍

（福建省地震局，福州 350003）

0 引言

從2000 年以來，華南各?。ㄗ砸眳^(qū)）陸續(xù)形成以洞體應(yīng)變儀、垂直擺傾斜儀、水管傾斜儀、重力儀、體積式鉆孔應(yīng)變儀為主的定點(diǎn)形變監(jiān)測臺網(wǎng)[1-2]。這些儀器在運(yùn)行期間，記錄了大量觀測數(shù)據(jù)，并積累了豐富的震例總結(jié)資料[3-4]。為了提高華南地區(qū)形變資料在地震預(yù)測上的應(yīng)用效能，在歷史經(jīng)驗(yàn)和多年日常工作經(jīng)驗(yàn)積累的基礎(chǔ)上，對華南定點(diǎn)形變資料進(jìn)行梳理。華南地區(qū)是典型的少震弱震地區(qū)，形變小組收集整理多年的干擾信息、氣象信息和儀器維修記錄，依據(jù)產(chǎn)生的異常形態(tài)分類，歸納總結(jié)出各測項(xiàng)的干擾庫和異常庫，初步建設(shè)華南地區(qū)定點(diǎn)形變預(yù)測指標(biāo)。使用全國形變學(xué)科組推廣的自動(dòng)算法及R值預(yù)測效能檢驗(yàn)再次梳理預(yù)測指標(biāo)，將兩次指標(biāo)結(jié)論進(jìn)行對比印證，利用自動(dòng)算法可有效提高工作效率和準(zhǔn)確性，以漳州洞體應(yīng)變東西分量為例，其自動(dòng)算法掃描結(jié)果與經(jīng)驗(yàn)掃描結(jié)果基本一致[5]。

1 華南形變觀測資料和地震概況

1.1 觀測資料概況

華南形變指標(biāo)建設(shè)采用華南地區(qū)包括福建、廣東、廣西、江西、海南和湖南六省/自冶區(qū)定點(diǎn)形變觀測數(shù)據(jù)，參與定點(diǎn)形變指標(biāo)梳理的臺站和顯著地震分布見圖1。因形變異常判定至少需要累積3年以上平穩(wěn)運(yùn)行數(shù)據(jù)，所以華南形變指標(biāo)建設(shè)基于穩(wěn)定運(yùn)行3年以上的定點(diǎn)形變觀測資料分鐘采樣原始數(shù)據(jù)，共梳理了34套觀測儀器（65個(gè)測項(xiàng)分量）的數(shù)據(jù)，同一套觀測儀器中東西向和北南向在梳理過程過總結(jié)出不同的預(yù)測指標(biāo)，所以以測量分量為形變指標(biāo)數(shù)量單位。其中垂直擺傾斜儀（30個(gè)測項(xiàng)）、水管傾斜儀（20 個(gè)測項(xiàng)分量）、洞體應(yīng)變儀（20 個(gè)測項(xiàng)分量）為主要觀測手段，鉆孔應(yīng)變儀、定點(diǎn)跨斷層、水平擺等儀器觀測數(shù)量較少。此外少量觀測時(shí)長不足或受干擾嚴(yán)重的定點(diǎn)形變資料沒有參與。重點(diǎn)梳理華南地區(qū)ML4.5以上地震相關(guān)的前兆異常指標(biāo)，參與指標(biāo)梳理的測項(xiàng)和地震分布如圖1所示。各測項(xiàng)根據(jù)經(jīng)驗(yàn)和映震效能分別給出映震范圍，基于映震效能最佳的基礎(chǔ)上，統(tǒng)計(jì)得到大部分臺項(xiàng)的映震敏感區(qū)域是測點(diǎn)周邊200～300 km范圍內(nèi)。

圖1 華南地區(qū)顯著地震（2001—2007年）和定點(diǎn)形變測項(xiàng)分布圖Fig.1 Diagram of significant earthquakes from 2001 to 2007 and fixed-point deformation measurement items in South China

1.2 地震概況

參與研究的地震是2001—2017 年，華南形變觀測臺網(wǎng)內(nèi)ML4.5 以上地震（湖南為ML4.0 以上）。華南是少震、弱震地區(qū)，研究時(shí)段區(qū)域17 年內(nèi)僅發(fā)生7次ML5以上中強(qiáng)地震，分別是2004 年9月17日廣東陽江5.2 級地震、2005 年11 月26 日九江6.0級地震、2007 年 3 月 13 日福建順昌 5.1 和 5.0 級雙震、2012 年2 月16 日廣東新豐江5.1 級地震、2013年 2 月 22 日河源 5.0 級地震、2013 年 9 月 4 日福建仙游5.2級地震和2016 年7月31日蒼梧5.4級地震。利用自動(dòng)算法再次梳理福建地區(qū)2018—2021 年數(shù)據(jù)，期間發(fā)生了2018 年11 月的臺灣海峽6.2 級地震，其震中距離海岸線僅109 km，距震中150～210 km 內(nèi)的形變觀測點(diǎn)如福建漳州、福建泉州、廣東潮州等均有明顯同震反應(yīng)。

2 華南形變預(yù)測指標(biāo)建設(shè)思路

2017 年開展的華南地區(qū)預(yù)測指標(biāo)梳理，將華南定點(diǎn)形變異常按曲線形態(tài)分為破年變、趨勢轉(zhuǎn)折、速率改變、潮汐畸變、大幅波動(dòng)等5 類異常。2020 年中國地震局地震預(yù)測研究所印發(fā)的《震情會(huì)商技術(shù)方法業(yè)務(wù)應(yīng)用推薦清單》中，定點(diǎn)形變的異常類型被再次分類，將潮汐畸變異常、大幅波動(dòng)異常和速率改變異常這3 項(xiàng)合并稱為速率改變異常，保留破年變異常和趨勢轉(zhuǎn)折異常。本文開展的基于經(jīng)驗(yàn)掃描的工作早于《震情會(huì)商技術(shù)方法業(yè)務(wù)應(yīng)用推薦清單》印發(fā)時(shí)間，所以仍以5 種異常分類進(jìn)行介紹。

華南地區(qū)所有定點(diǎn)形變觀測均設(shè)在專業(yè)的綜合地震臺內(nèi)，有臺站專職工作人員值守，可以隨時(shí)關(guān)注儀器運(yùn)行狀態(tài)和環(huán)境干擾情況，大部分儀器均可以保證正常觀測。但因城市化進(jìn)程加大，原本地處較偏僻的觀測臺站，現(xiàn)頻繁受到高速公路、住房基建、打井抽水或供電故障等等多種人為干擾影響，沿海臺站受臺風(fēng)、強(qiáng)降雨干擾明顯，且華南地區(qū)氣候潮濕，山洞內(nèi)潮濕滲水明顯，儀器受潮后感應(yīng)器反應(yīng)不靈敏或者過度反應(yīng)時(shí)有發(fā)生，在復(fù)雜多樣的干擾源影響下，觀測過程產(chǎn)生大量異常數(shù)據(jù)，所以華南地區(qū)形變指標(biāo)體系梳理工作第一步需將干擾數(shù)據(jù)與前兆異常進(jìn)行區(qū)分。

基于華南地區(qū)地殼形變定點(diǎn)觀測數(shù)據(jù)的特殊性，形變清理小組在進(jìn)行異常指標(biāo)體系清理工作過程中，收集資料包括：原始數(shù)據(jù)、周邊顯著地震、有記錄以來的氣象數(shù)據(jù)、臺站工作日志、儀器維修記錄等，同時(shí)建立干擾庫與異常庫，達(dá)到定點(diǎn)形變異常指標(biāo)清理的目的，具體包括以下步驟：

步驟一：掃描各測項(xiàng)分量分鐘采樣曲線，發(fā)現(xiàn)潮汐畸變、大幅波動(dòng)、速率改變等中短期異常時(shí)，則查找同一時(shí)期的氣象數(shù)據(jù)、臺站工作日志和儀器維修記錄等觀測資料，以判別該中短期異常是否是由于氣象、環(huán)境變化或儀器故障等因素造成；

步驟二：數(shù)據(jù)異常因干擾而產(chǎn)生，將干擾源起止時(shí)間、受影響起止時(shí)間和受干擾后產(chǎn)生的異常形態(tài)種類記入干擾庫；

步驟三：不因任何干擾而起的數(shù)據(jù)異常，將異常形態(tài)種類和異常起止時(shí)間記入異常庫；

步驟四：逐條掃描日均值采樣曲線，發(fā)現(xiàn)趨勢轉(zhuǎn)折、破年變等中長期異常時(shí)，多次重復(fù)步驟2和步驟3，直至將所有干擾異常與前兆異常分別錄入干擾庫與異常庫；統(tǒng)計(jì)各臺項(xiàng)不同類型異常在地震前出現(xiàn)的頻次，得到不同指標(biāo)的預(yù)報(bào)效能。

圖2 定點(diǎn)形變指標(biāo)體系建設(shè)流程圖Fig.2 The flowchart for fixed-point deformation prediction index system construction

梳理異常高度依賴經(jīng)驗(yàn)性工作，是擺在指標(biāo)清理工作之前的重要預(yù)備階段，目前需要解決的關(guān)鍵性問題有兩個(gè)：

（1）判定數(shù)據(jù)形態(tài)是否異常，需要長期的預(yù)報(bào)工作經(jīng)驗(yàn)，且熟悉異常形態(tài)的分類原則。預(yù)報(bào)工作人員可通過經(jīng)驗(yàn)判斷方式掃描后，再次利用形變學(xué)科組推廣使用的自動(dòng)算法進(jìn)行對照檢驗(yàn)。

（2）判斷曲線形態(tài)異常是否由干擾引起，解決方案是盡可能多地收集當(dāng)?shù)貧庀筚Y料、儀器維護(hù)資料、異常核實(shí)報(bào)告，日常跟蹤數(shù)據(jù)，發(fā)現(xiàn)異常及時(shí)進(jìn)行排查，詳細(xì)記錄人為干擾和周邊環(huán)境干擾。進(jìn)一步規(guī)范臺站日志和省局監(jiān)控日志，將前兆數(shù)據(jù)跟蹤分析與干擾庫的異常判定規(guī)則相結(jié)合，將有助于今后異常識別向標(biāo)準(zhǔn)化發(fā)展。

3 定點(diǎn)形變預(yù)測指標(biāo)梳理成果[6]

3.1 形變干擾信息庫建設(shè)

各測項(xiàng)的干擾庫記錄內(nèi)容包括：臺站代碼、測項(xiàng)及測項(xiàng)分量、干擾開始時(shí)間、干擾結(jié)束時(shí)間、干擾描述、干擾產(chǎn)生的異常類別等信息。此步驟排除了大量有明顯干擾源的數(shù)據(jù)異常，是前兆指標(biāo)梳理的重要預(yù)備步驟。

由大量統(tǒng)計(jì)可見，如臺風(fēng)等氣象因素干擾，通常影響時(shí)長在7 d 以內(nèi)，會(huì)產(chǎn)生潮汐畸變等短期異常形態(tài)；如施工爆破等干擾，通常產(chǎn)生速率改變或大幅波動(dòng)等較顯著的短期異常形態(tài)；如測點(diǎn)周邊抽水干擾或大型基建干擾，通常造成趨勢轉(zhuǎn)折或破年變等長周期異常形態(tài)。定點(diǎn)形變測項(xiàng)的干擾庫，不僅是人工識別異常的基礎(chǔ)，而且也是異常自動(dòng)算法軟件運(yùn)行的前提，排除有明確干擾源的曲線異常時(shí)段，在大量曲線變化中甄別出由地殼運(yùn)動(dòng)產(chǎn)生的形變異常信息，是建設(shè)高效能的預(yù)測指標(biāo)體系的必要條件。

3.2 中期預(yù)測指標(biāo)部分實(shí)例

廣東河源室內(nèi)短水準(zhǔn)（瓷標(biāo)和鋼標(biāo)）為人工觀測，是華南地區(qū)連續(xù)觀測時(shí)長最長的定點(diǎn)形變測項(xiàng)[7]，利用了自1990年以來持續(xù)產(chǎn)出觀測數(shù)據(jù)，自觀測以來河源臺周邊50 km范圍內(nèi)共發(fā)生7次ML4.5以上地震，排除2003—2009 年周邊高速公路修建對觀測的影響后，其中4 次地震前12 月內(nèi)，水準(zhǔn)儀均存在明顯破年變異常現(xiàn)象，河源室內(nèi)短水準(zhǔn)破年變異常對應(yīng)周邊50 km 范圍內(nèi)ML4.5以上地震，地震對應(yīng)率達(dá)80%。

梳理江西南昌垂直擺傾斜儀東西分量破年變異常[8]，得到其對應(yīng)周邊300 km 范圍內(nèi)ML4.5 以上地震對應(yīng)率達(dá)30%。

梧州水管傾斜儀北南分量在梧州ML5.6 地震和北流ML5.2 地震前半年內(nèi)，均發(fā)生趨勢轉(zhuǎn)折異常，除儀器維修該測項(xiàng)分量觀測以來僅發(fā)生3次趨勢轉(zhuǎn)折，梧州水管傾斜儀北南分量趨勢轉(zhuǎn)折類異常對應(yīng)周邊200 km 范圍內(nèi)ML4.5 以上地震對應(yīng)率達(dá)66%。

3.3 短期預(yù)測指標(biāo)實(shí)例

福建漳州洞體應(yīng)變儀2007 年正式入網(wǎng)觀測，自觀測以來在漳州臺周邊200 km 范圍內(nèi)發(fā)生了漳州ML4.8地震、仙游ML4.5地震和仙游ML5.2地震等3 次中強(qiáng)以上地震，排除人為干擾、儀器故障及氣象干擾等因素后，漳州洞體應(yīng)變東西分量非干擾產(chǎn)生的潮汐畸變類異常6 次，其中2 次發(fā)生在地震前一個(gè)月內(nèi)，漳州洞體應(yīng)變東西分量潮汐畸變類異常對應(yīng)周邊200 km 范圍內(nèi)ML4.5以上地震對應(yīng)率達(dá)33%。

廈門臺洞體應(yīng)變儀2005 年正式入網(wǎng)觀測，自觀測以來，廈門臺周邊200 km 范圍內(nèi)發(fā)生了臺灣海峽南部ML4.6地震、永春ML4.9地震、漳州ML4.8地震、仙游ML4.5 地震和仙游ML5.2 地震等5 次中強(qiáng)以上地震，排除人為干擾、氣象和儀器故障等干擾因素后，廈門臺洞體應(yīng)變儀東西分量和北南分量同時(shí)發(fā)生速率改變類異常共發(fā)生20次，其中8次速率改變類異常發(fā)生在上述地震前5個(gè)月內(nèi)，廈門臺洞體應(yīng)變儀速率改變類異常對應(yīng)周邊200 km范圍內(nèi)ML4.5及以上地震對應(yīng)率達(dá)40%。

3.4 指標(biāo)建設(shè)小結(jié)

通過對華南地區(qū)洞體應(yīng)變、體應(yīng)變、垂直擺傾斜儀、水管傾斜儀、連續(xù)重力和定點(diǎn)跨斷層等定點(diǎn)形變觀測數(shù)據(jù)梳理后，梳理得到定點(diǎn)形變的干擾庫和異常庫，表1列出預(yù)測效能較好的華南定點(diǎn)形變預(yù)測指標(biāo)參數(shù)，從定點(diǎn)形變異常指標(biāo)體系清理結(jié)果看，初步可得到以下幾點(diǎn)認(rèn)識：

表1 華南地區(qū)定點(diǎn)形變觀測預(yù)報(bào)指標(biāo)效能統(tǒng)計(jì)表Table 1 Statistical table of the effectiveness of fixed-point deformation prediction indexes

（1）參與定點(diǎn)形變指標(biāo)建設(shè)共65 個(gè)測項(xiàng)，得到3個(gè)月以內(nèi)的短期異常指標(biāo)7項(xiàng)，3個(gè)月至1年的中長期預(yù)測指標(biāo)6項(xiàng)。獲得的指標(biāo)較少，主要原因是數(shù)字化觀測近20 年來，華南形變測網(wǎng)內(nèi)發(fā)生中強(qiáng)地震頻次較低，次要原因是儀器敏感度高且受干擾明顯，導(dǎo)致部分儀器恰好在周邊地震前受到干擾或儀器發(fā)生故障，無法判斷曲線畸變是否前兆異常。

（2）總結(jié)出的13 項(xiàng)指標(biāo)其中8 項(xiàng)指標(biāo)由洞體應(yīng)變分量總結(jié)得出，可以看出洞體應(yīng)變儀運(yùn)行相對穩(wěn)定，映震效果較好，在潮汐畸變、速率變化、趨勢轉(zhuǎn)折、破年變等異常梳理中均可以發(fā)現(xiàn)相應(yīng)預(yù)測指標(biāo)；各觀測臺的洞體應(yīng)變映震指標(biāo)并不一致，且同一臺站洞體應(yīng)變儀不同分量預(yù)測指標(biāo)也不完全相同；

（3）趨勢轉(zhuǎn)折和破年變異常對應(yīng)中長期預(yù)測指標(biāo)較多，破年變類異常的預(yù)測時(shí)間從異常開始至發(fā)震 6 個(gè)月～1 年不等，預(yù)測空間為 200～300 km 范圍內(nèi)，預(yù)測震級以5級地震為主，同時(shí)也需關(guān)注附近200 km 范圍內(nèi)ML4.0以上地震的危險(xiǎn)性。而潮汐畸變、大幅波動(dòng)等短期影響曲線形態(tài)的異常類型較多地對應(yīng)中短期預(yù)測指標(biāo)。

（4）河源的室內(nèi)水準(zhǔn)儀屬于華南僅有的多年持續(xù)觀測定點(diǎn)跨斷層儀器（數(shù)據(jù)），受干擾暫停觀測一段時(shí)間，恢復(fù)觀測后其觀測數(shù)據(jù)質(zhì)量并沒降低，破年變異常對應(yīng)周邊50 km 范圍內(nèi)的ML4.0 以上地震，映震效能較好。

（5）與非沿海地區(qū)相比，沿海地區(qū)觀測儀器受環(huán)境干擾普遍更嚴(yán)重，儀器折舊老化也明顯快于非沿海地區(qū)，因維修和干擾導(dǎo)致多數(shù)儀器長期趨勢規(guī)律性被打亂，所以更適合提取短期異常預(yù)測指標(biāo)，而非沿海地區(qū)觀測數(shù)據(jù)相對穩(wěn)定，可提取到一些中期預(yù)測指標(biāo)。

（6）因華南地震頻次和異常頻次較低，部分測項(xiàng)指標(biāo)信度不高，運(yùn)用于日常會(huì)商時(shí)，應(yīng)該將多臺多測項(xiàng)的觀測資料進(jìn)行綜合跟蹤識別，排除儀器干擾、環(huán)境干擾、人為干擾，將形變、電磁、流體等不同學(xué)科長中短臨異常進(jìn)行綜合研判[9]。

4 定量化預(yù)測指標(biāo)梳理

形變異常曲線形態(tài)的辨識，目前大多基于分析人員主觀經(jīng)驗(yàn)，不同的分析人員對于同一個(gè)異常的判定持不同的意見，異常的判定準(zhǔn)則缺乏定量化的表達(dá)，判定結(jié)果的穩(wěn)定性和可重復(fù)性較低。隨著前兆臺網(wǎng)的密集布設(shè)，采樣率的提高，各測項(xiàng)多年時(shí)間序列數(shù)據(jù)的積累量達(dá)萬億字節(jié)，為實(shí)現(xiàn)海量觀測數(shù)據(jù)的快速高效提取，迫切需要針對不同異常類型、推廣切實(shí)可行的自動(dòng)識別算法。李宏偉、閆偉等研究了定點(diǎn)形變趨勢轉(zhuǎn)折異常的自動(dòng)算法，大大提高了人工判別的工作效率，可快速掃描大量時(shí)序數(shù)據(jù)，從統(tǒng)計(jì)角度定位令指標(biāo)的預(yù)測效能最佳的異常閾值[10]。苑爭一等研究了基于奇異譜分析算法，在去除典型干擾及長周期趨勢變化的基礎(chǔ)上，擬合觀測資料背景年變序列，自動(dòng)提取破年變異常特征時(shí)段[11]。此外地震預(yù)測要求盡可能準(zhǔn)確和精確，即與實(shí)際地震相符且預(yù)測窗口期盡可能短，以利于采取應(yīng)急措施的同時(shí)盡可能地縮小影響范圍，目前自動(dòng)掃描算法大都結(jié)合震例進(jìn)行動(dòng)態(tài)R值檢驗(yàn)，遞歸求解出了最高R值評分對應(yīng)的異常判定閾值[12-13]。華南預(yù)測指標(biāo)采用R值評分作為預(yù)測有效性檢驗(yàn)的具體方法，即扣除隨機(jī)概率的地震預(yù)測效能的R值評分方法：

③ 墻體變形過大原因分析：支撐架設(shè)滯后(且未完全到位)，每道鋼支撐架設(shè)都存在滯后情況；土方開挖過程，基坑內(nèi)積水情況較嚴(yán)重，基坑內(nèi)(含基坑外)的積水沿連續(xù)墻向下滲透，造成土體含水率增加，土體的工程性能降低，特別是造成嵌固端的土體嵌固力減弱，因而墻體向基坑內(nèi)的變形過大。

R＝報(bào)對的地震次數(shù)/應(yīng)預(yù)報(bào)的地震總次數(shù)-預(yù)報(bào)占用時(shí)間/預(yù)報(bào)研究的總時(shí)間

以福建漳州洞體應(yīng)變東西分量為例，采用速率變化異常判定自動(dòng)算法，因漳州臺2019 年前后周邊施工干擾明顯，所以采用2008 年至2019 年漳州洞體應(yīng)變東西分量原始數(shù)據(jù)和漳州臺干擾信息表，通過遍歷掃描0°～90°速率變化范圍，遞歸求解出令R值最高的速率矢量轉(zhuǎn)角，得到漳州臺洞體應(yīng)變東西分量速率變化異常的預(yù)測規(guī)則：當(dāng)速率矢量轉(zhuǎn)角超過閾值2.06°時(shí)，漳州臺周邊200 km 范圍內(nèi)，30 d內(nèi)發(fā)生ML4.5以上的地震，預(yù)測效能R：0.65，R0：0.47，表明漳州洞體應(yīng)變東西分量速率變化的預(yù)測值大于自然概率值，可做為日常跟蹤的預(yù)測指標(biāo)。

2019 年用自動(dòng)算法掃描的福建定點(diǎn)形變測項(xiàng)共有41 個(gè)，每個(gè)測項(xiàng)分量分別利用破年變、趨勢轉(zhuǎn)折、速率變化3個(gè)異常分類自動(dòng)算法掃描，部分常年沒有正常年變的測項(xiàng)不掃描破年變異常，受干擾明顯的測項(xiàng)也不參與異常掃描，最終各測項(xiàng)利分別利用不同算法共掃描110次，得到效能評估R值明顯高于自然概率R0值的指標(biāo)共計(jì)21 條，即通過自動(dòng)算法效能檢驗(yàn)的指標(biāo)。

從福建的自動(dòng)異常掃描結(jié)果可以看出，使用自動(dòng)算法進(jìn)行預(yù)測具有以下特點(diǎn)：

（1）結(jié)論具有可重復(fù)性，在相同的干擾庫基礎(chǔ)上運(yùn)行同一算法，不同的操作者可產(chǎn)生相同的預(yù)報(bào)效果；

（2）預(yù)測三要素?cái)?shù)據(jù)精確度高，為保證預(yù)測效果最優(yōu)化，算法會(huì)根據(jù)不同觀測時(shí)長的數(shù)據(jù)精確調(diào)整預(yù)測要素；

（3）虛報(bào)率較高，目前的算法無法兼顧較低的虛報(bào)率和較高的報(bào)準(zhǔn)率，導(dǎo)致虛報(bào)率過高。日常使用過程中，自動(dòng)算法得到的異常往往只能參考使用，自動(dòng)算法還有待優(yōu)化和改進(jìn)。

（4）對干擾庫依賴性高，自動(dòng)算法預(yù)測干擾庫基礎(chǔ)上才能發(fā)揮優(yōu)勢，完成干擾庫標(biāo)準(zhǔn)化建設(shè)，是自動(dòng)算法得以完善的前提。而目前的干擾庫建設(shè)高度依賴經(jīng)驗(yàn)性總結(jié)，無法實(shí)現(xiàn)標(biāo)準(zhǔn)化和自動(dòng)化，所以現(xiàn)階段自動(dòng)算法掃描可做為經(jīng)驗(yàn)判別方法的有益補(bǔ)充，并不能完全替代經(jīng)驗(yàn)判斷。

5 結(jié)論

2017 年以來由華南地區(qū)各省/自冶區(qū)形變分析預(yù)報(bào)人員共同努力，完成華南定點(diǎn)形變預(yù)測指標(biāo)體系建設(shè)，通過遍歷掃描所有連續(xù)觀測資料較好的定點(diǎn)形變數(shù)據(jù)，統(tǒng)計(jì)各測項(xiàng)自觀測以來的儀器故障、人為干擾和氣象干擾等各種因素，詳細(xì)歸納不同時(shí)段受影響的異常判定類型，總結(jié)所到各測項(xiàng)分量的干擾信息庫。在干擾庫的基礎(chǔ)上，統(tǒng)計(jì)不同類型異常在測點(diǎn)周邊200～300 km 范圍內(nèi)對應(yīng)ML4.5 以上地震的概率，得到預(yù)測效能較好的定點(diǎn)形變預(yù)測指標(biāo)。

通過華南定點(diǎn)形變預(yù)測指標(biāo)的建設(shè)，總結(jié)在少震弱震的華南地區(qū)的形變預(yù)測指標(biāo)特征，筆者分析認(rèn)為洞體應(yīng)變是華南預(yù)測指標(biāo)效能較的觀測儀器，此外模擬觀測的定點(diǎn)跨斷層也可得到較好的中期預(yù)測指標(biāo)。以福建漳州洞體應(yīng)變東西分量為例，再次用《地震分析會(huì)商技術(shù)方法列裝清單》方法，采用全國學(xué)科組推薦使用的自動(dòng)識別算法進(jìn)行檢查，掃描福建各測項(xiàng)的預(yù)測方法效能，能夠通過R值效能檢驗(yàn)的方法與人工判別均認(rèn)為漳州洞體應(yīng)變東西向速率變化異常指標(biāo)較可靠，自動(dòng)軟件掃描的結(jié)果可形成定量化的預(yù)測指標(biāo)。在地震觀測海量數(shù)據(jù)的背景下，使用自動(dòng)算法完成地震預(yù)測工作是提高分析預(yù)報(bào)人員工作效率的必然選擇。

此次華南地區(qū)定點(diǎn)形變預(yù)測指標(biāo)清理中，各測項(xiàng)的預(yù)測空間多以測點(diǎn)為中心，統(tǒng)計(jì)測點(diǎn)周邊近圓形空間區(qū)域發(fā)生的地震，但是大量震例表明地震的孕育和發(fā)生與斷層運(yùn)動(dòng)存在密切關(guān)系，今后應(yīng)考慮測點(diǎn)周邊斷層的分布及走向，作為測點(diǎn)映震區(qū)域的劃定依據(jù)[14]。

致謝：華南各省/自冶區(qū)預(yù)報(bào)工作人員趙愛平、黃元敏、 譚爭光、文翔、原永東、李盛等提供了所在省/自冶區(qū)的定點(diǎn)形變觀測數(shù)據(jù)，共同參與了華南定點(diǎn)形變預(yù)測指標(biāo)體系建設(shè)工作。