鄧世廣 王月 馬亞偉

中國(guó)地震臺(tái)網(wǎng)中心,北京 100045

0 引言

地震目錄是地震危險(xiǎn)性分析中的關(guān)鍵資料,也是研究巖石圈動(dòng)力學(xué)過程的重要數(shù)據(jù),在國(guó)民經(jīng)濟(jì)建設(shè)、地震分析預(yù)報(bào)、工程抗震、地震減災(zāi)等領(lǐng)域發(fā)揮著重要的作用(李圣強(qiáng)等,2004；王海濤等,2006；張衛(wèi)華等,2007；何琳等,2017)。地震目錄又是一種時(shí)效性要求較高的數(shù)據(jù),特別是在強(qiáng)震發(fā)生后的應(yīng)急工作中,需要及時(shí)更新地震目錄,跟蹤地震序列的發(fā)展演化。在日常的地震活動(dòng)性分析工作中,為了跟蹤分析最新的地震活動(dòng)情況,也需要頻繁地更新地震目錄。準(zhǔn)確、實(shí)時(shí)和完整地獲取地震目錄對(duì)地震科研人員以及防震減災(zāi)、城建規(guī)劃等工作都至關(guān)重要(張衛(wèi)華等,2007)。

由于歷史原因,地震目錄格式多種多樣,如CapSeis系統(tǒng)使用的EQT或EQ2格式、SSEPC所用的WKF或SHK格式、地震速報(bào)目錄、地震月報(bào)目錄、Apnet網(wǎng)絡(luò)報(bào)送目錄及其他特征格式等(李圣強(qiáng)等,2004)。為了支持更高精度的地震定位結(jié)果,同時(shí)包含地震參數(shù)數(shù)據(jù)和震源機(jī)制解數(shù)據(jù),李圣強(qiáng)等(2004)設(shè)計(jì)了EQ3格式的地震目錄,這種格式由2個(gè)相互關(guān)聯(lián)的二進(jìn)制結(jié)構(gòu)來表示記錄,與ASCII碼格式相比,這種二進(jìn)制結(jié)構(gòu)的地震目錄格式具有數(shù)據(jù)存儲(chǔ)量少、存取訪問速度快的優(yōu)勢(shì)。由于地震目錄格式類型多樣,不同機(jī)構(gòu)提供的地震目錄格式不盡相同(馬寶君等,2005),而不同地震分析系統(tǒng)軟件支持的地震目錄格式也各有不同,因此，不同地震目錄格式之間的轉(zhuǎn)換是地震活動(dòng)性分析中經(jīng)常遇到的問題。

為提高地震目錄格式轉(zhuǎn)換及更新的工作效率,前人做了較為豐富的工作。考慮到數(shù)據(jù)交換或不同地震系統(tǒng)軟件分析處理的需要,李圣強(qiáng)等(2004)分析了應(yīng)用較多的目錄格式,并提供了格式轉(zhuǎn)換工具，可以將多種格式的地震目錄轉(zhuǎn)換為EQ3格式,同時(shí)也支持EQ3格式地震目錄轉(zhuǎn)換為其他格式地震目錄。馬寶君等(2005)使用VB編制了地震目錄快速轉(zhuǎn)換軟件,將網(wǎng)上檢索的地震目錄轉(zhuǎn)換成TXT、EQT、WKF等格式目錄。為了進(jìn)一步提高地震目錄更新的時(shí)效性,張衛(wèi)華等(2007)開發(fā)了一種能從FTP、HTTP服務(wù)器上直接下載地震目錄資料的軟件,實(shí)現(xiàn)了地震目錄的自動(dòng)下載并轉(zhuǎn)換為EQT、TXT、WKF、CSV等地震目錄格式。由于SuperSeis提供的地震目錄格式轉(zhuǎn)換工具支持的地震目錄格式有限,于洪波(2009)采用Delphi語言實(shí)現(xiàn)了Excel和TXT格式地震目錄向EQ3格式的轉(zhuǎn)換,有效提高了震后趨勢(shì)會(huì)商的速度和準(zhǔn)確性,但該方法沒有考慮EQ3格式地震目錄的輔助記錄部分,因此未將發(fā)震地點(diǎn)寫入EQ3地震目錄內(nèi)。梁阿全等(2016)應(yīng)用Excel VBA實(shí)現(xiàn)了從Excel、Q01、EQT等格式至WKF格式之間的轉(zhuǎn)換。何琳等(2017)通過對(duì)安徽地震編目數(shù)據(jù)庫(kù)的訪問,實(shí)現(xiàn)了地震目錄查詢以及數(shù)據(jù)提取、轉(zhuǎn)換和導(dǎo)出等功能,減輕了數(shù)據(jù)維護(hù)人員的工作量,提高了數(shù)據(jù)的即時(shí)性和準(zhǔn)確性,能夠有效滿足地震應(yīng)急的迫切要求。由于EQT、TXT、WKF、CSV等地震目錄格式屬于ASCII碼格式,讀寫相對(duì)簡(jiǎn)單,因此針對(duì)這些地震目錄的格式轉(zhuǎn)換方法研究相對(duì)較多。EQ3地震目錄格式屬于二進(jìn)制結(jié)構(gòu),具有數(shù)據(jù)存儲(chǔ)量少、存取訪問速度快的優(yōu)勢(shì),但這種結(jié)構(gòu)讀寫方式相對(duì)復(fù)雜,目前還沒有針對(duì)EQ3格式地震目錄的自動(dòng)更新方法。

隨著地震監(jiān)測(cè)定位水平的不斷提高,中國(guó)地震編目網(wǎng)站提供的地震目錄中經(jīng)緯度的精度可達(dá)到小數(shù)點(diǎn)后3位,更高精度的地震經(jīng)緯度能夠更加清晰地刻畫地震分布、揭示斷裂構(gòu)造形態(tài),有助于地震活動(dòng)性的深入分析。目前常用的EQT、WKF、CSV等ASCII碼格式的地震目錄,一般僅支持經(jīng)緯度精確到小數(shù)點(diǎn)后2位,難以滿足更加精細(xì)的地震活動(dòng)性分析工作。而EQ3格式的地震目錄屬于二進(jìn)制結(jié)構(gòu),震中經(jīng)緯度參數(shù)采用Float類型,其參數(shù)值精度保持在7位有效數(shù)值的精度范圍內(nèi),經(jīng)緯度的有效精度可以表示到小數(shù)點(diǎn)后4位,即相當(dāng)于10m級(jí)(李閩峰等,2006)。因此,EQ3格式地震目錄是一種能夠適應(yīng)更加精細(xì)地震活動(dòng)性分析工作的地震目錄格式,具有更好的應(yīng)用前景。然而,EQ3格式地震目錄雖然能夠支持較高的地震定位精度,但目前尚沒有針對(duì)高精度定位結(jié)果的地震目錄的格式轉(zhuǎn)換方法,高精度的地震目錄還未能在地震活動(dòng)性分析工作中充分發(fā)揮作用。

為實(shí)現(xiàn)高定位精度EQ3格式地震目錄的自動(dòng)更新,設(shè)計(jì)并實(shí)現(xiàn)了地震目錄自動(dòng)更新程序。該程序首先基于Python爬蟲技術(shù)實(shí)現(xiàn)地震目錄的自動(dòng)下載；然后針對(duì)地震編目系統(tǒng)中多個(gè)省局的快報(bào)地震目錄中存在重復(fù)定位地震的情況,提出了一種篩除重復(fù)定位地震的方法；最后研究分析EQ3格式地震目錄的二進(jìn)制結(jié)構(gòu),在此基礎(chǔ)上實(shí)現(xiàn)了高精度定位結(jié)果的地震目錄的格式轉(zhuǎn)換、目錄拼接等功能。實(shí)際應(yīng)用表明，該程序能夠完成高精度定位結(jié)果的地震目錄的自動(dòng)更新,提高了地震目錄更新的工作效率,能夠滿足地震應(yīng)急對(duì)地震目錄的高實(shí)效性要求。

1 技術(shù)方法

為實(shí)現(xiàn)高精度定位結(jié)果的地震目錄的自動(dòng)下載,首先調(diào)用預(yù)設(shè)的地震目錄下載相關(guān)參數(shù),利用Python爬蟲技術(shù)從地震編目網(wǎng)站自動(dòng)下載指定的地震目錄,針對(duì)鄰省快報(bào)目錄存在重復(fù)定位地震的問題,采用本文提出的重復(fù)定位地震篩除方法對(duì)地震目錄內(nèi)的重復(fù)定位地震進(jìn)行篩除,保證地震目錄的唯一性,最后分析EQ3格式地震目錄的數(shù)據(jù)結(jié)構(gòu),實(shí)現(xiàn)地震目錄向EQ3、EQT等格式的轉(zhuǎn)換,整體技術(shù)流程如圖1 所示。需要說明的是,地震編目網(wǎng)站提供的正式地震目錄的更新時(shí)間是滯后的,而快報(bào)目錄的實(shí)時(shí)性較好。為了提高地震目錄的完整性,部分地震目錄使用最近一段時(shí)間的快報(bào)目錄來補(bǔ)充正式目錄的不足。為了在目錄更新時(shí)將快報(bào)目錄替換為正式目錄,下載地震目錄的起始時(shí)間應(yīng)早于待更新目錄內(nèi)省快報(bào)地震目錄內(nèi)最早的發(fā)震時(shí)刻。另外,為了保證地震目錄更新時(shí)不遺漏地震事件,一般也將下載地震目錄的起始時(shí)間設(shè)置早于地震目錄的最近1次地震發(fā)震時(shí)刻。在地震目錄更新時(shí),程序?qū)⒆詣?dòng)刪除待更新地震目錄內(nèi)下載起始時(shí)間之后的地震目錄,再將其與下載的最新地震目錄進(jìn)行拼接。

圖1 地震目錄自動(dòng)更新技術(shù)流程

1.1 基于Python爬蟲技術(shù)的高精度定位結(jié)果的地震目錄自動(dòng)下載

隨著互聯(lián)網(wǎng)的快速發(fā)展,爬蟲技術(shù)已經(jīng)成為高效準(zhǔn)確地獲取數(shù)據(jù)的重要方式之一,在各個(gè)領(lǐng)域都有著廣泛的應(yīng)用(周中華等,2014；熊暢,2017；余豪士等,2018)。Python中的第三方HTTP庫(kù)——Requests庫(kù)被爬蟲開發(fā)者廣泛應(yīng)用(余豪士等,2018)。Requests集成了定制請(qǐng)求頭部、發(fā)送請(qǐng)求、傳遞URL參數(shù)、獲取相應(yīng)內(nèi)容等多種函數(shù)(謝克武,2017)。通過Requests庫(kù),能夠方便實(shí)現(xiàn)網(wǎng)絡(luò)數(shù)據(jù)的爬取。

目前,我國(guó)地震活動(dòng)性分析使用的地震目錄主要通過中國(guó)地震臺(tái)網(wǎng)中心的地震編目系統(tǒng)網(wǎng)站下載獲取。地震編目系統(tǒng)網(wǎng)站主要通過GET方式發(fā)送請(qǐng)求、傳遞參數(shù)。Requests集成了多種HTTP協(xié)議中的請(qǐng)求方式,包括最普遍的GET和POST請(qǐng)求。因此，利用Requests庫(kù)可以方便地訪問地震編目網(wǎng)站并獲取指定的地震目錄數(shù)據(jù)。

地震編目系統(tǒng)網(wǎng)站的用戶登陸通過GET方式完成,實(shí)現(xiàn)過程如下：

s=requests. Session()

login_params={′tb000.t1′:username,′tb000.t3′:password}

s.get(“http://10.5.160.18/console/login.action”,params=login_params)

登陸之后即可繼續(xù)通過GET方式傳遞下載地震目錄的參數(shù),從而獲取地震目錄。下載地震目錄的參數(shù)涉及地震目錄類型、地震目錄格式、震級(jí)類型、起止時(shí)間、震級(jí)范圍、經(jīng)緯度范圍等,可以通過字典的形式構(gòu)建1個(gè)參數(shù)字典,隨后通過GET方式獲取下載參數(shù)規(guī)定范圍內(nèi)的地震目錄,實(shí)現(xiàn)過程如下：

DownloadParams={′key1′:′value1′,′key2′:′value2′,′key3′:′value3′,…}

response=s.get(“http://10.5.160.18/export/download.action”,params=DownloadParams)

值得注意的是,在地震編目系統(tǒng)網(wǎng)站提供的地震目錄格式中,通用格式文本記錄的地震經(jīng)緯度精度能夠達(dá)到小數(shù)點(diǎn)后3位,定位精度高于其他格式。如需下載高精度定位結(jié)果的地震目錄,應(yīng)將地震目錄格式參數(shù)設(shè)置為通用格式文本。

1.2 重復(fù)定位地震篩除

根據(jù)《測(cè)震臺(tái)網(wǎng)運(yùn)行管理細(xì)則》(1)中國(guó)地震局監(jiān)測(cè)預(yù)報(bào)司,2015,中震測(cè)函〔2015〕115號(hào)測(cè)震臺(tái)網(wǎng)運(yùn)行管理細(xì)則(2015修訂版)。,地震編目工作產(chǎn)出的數(shù)據(jù)包括快報(bào)和正式報(bào)目錄。省級(jí)臺(tái)網(wǎng)中心應(yīng)在每日12時(shí)前報(bào)送前1日本省及周邊地區(qū)的省局快報(bào)目錄,在每周日上午12時(shí)以前報(bào)送上周的省局正式目錄。國(guó)家測(cè)震臺(tái)網(wǎng)中心收集有人值守的國(guó)家測(cè)震臺(tái)站和省級(jí)臺(tái)網(wǎng)中心處理的地震目錄和觀測(cè)報(bào)告,并進(jìn)行全國(guó)統(tǒng)一編目,保證目錄和觀測(cè)報(bào)告的唯一性。每日24時(shí)前產(chǎn)出前1日的全國(guó)快報(bào)目錄,每周編輯國(guó)內(nèi)及全球地震目錄和觀測(cè)報(bào)告，于第3周的周日24時(shí)前發(fā)布全國(guó)正式目錄和觀測(cè)報(bào)告。由此可見,省局快報(bào)地震目錄的實(shí)時(shí)性最高,且在實(shí)際工作中,省級(jí)臺(tái)網(wǎng)中心一般會(huì)及時(shí)發(fā)布最近的地震事件至省局快報(bào)目錄內(nèi),因此省局快報(bào)的滯后時(shí)間一般僅在數(shù)小時(shí)內(nèi)。而其他類型的地震目錄則至少會(huì)滯后1天的時(shí)間,難以滿足地震活動(dòng)性分析的實(shí)時(shí)性要求。

為了跟蹤最新的地震活動(dòng),在地震目錄更新時(shí),需要將省局快報(bào)目錄進(jìn)行匯總來補(bǔ)充正式目錄尚未更新的部分。然而，由于各省級(jí)臺(tái)網(wǎng)中心會(huì)記錄本省周邊地區(qū)的地震事件,因此，鄰省之間的地震目錄會(huì)存在重復(fù)定位地震。由于直接獲得各臺(tái)站的實(shí)時(shí)波形數(shù)據(jù)存在困難,我們僅能根據(jù)地震目錄內(nèi)的信息來識(shí)別重復(fù)定位地震并予以剔除,從而保證地震目錄的唯一性。傳統(tǒng)的工作方式需要通過人工判斷,如果2個(gè)地震的發(fā)震時(shí)刻、震中位置以及震級(jí)均相近,且屬于不同省局上傳的地震事件,則認(rèn)為2個(gè)地震為重復(fù)定位地震并刪除。這種方式需要人工對(duì)每個(gè)地震事件逐一判斷,操作繁瑣且易出錯(cuò)。本文基于傳統(tǒng)人工識(shí)別重復(fù)定位地震的經(jīng)驗(yàn),建立重復(fù)定位地震的識(shí)別準(zhǔn)則,通過計(jì)算機(jī)高效地自動(dòng)識(shí)別并篩除重復(fù)定位地震,技術(shù)流程如圖2 所示?；镜闹貜?fù)定位地震識(shí)別準(zhǔn)則為：

圖2 重復(fù)定位地震識(shí)別方法技術(shù)流程

(1)發(fā)震時(shí)刻小于給定閾值；

(2)震中距小于給定閾值；

(3)震級(jí)差小于給定閾值。

如果2次地震同時(shí)滿足以上3點(diǎn)標(biāo)準(zhǔn),即認(rèn)為是重復(fù)定位地震。閾值需要根據(jù)實(shí)際地震監(jiān)測(cè)水平設(shè)定。例如,當(dāng)中國(guó)臺(tái)灣發(fā)生較強(qiáng)地震時(shí),多個(gè)省局會(huì)監(jiān)測(cè)并記錄該地震,但目前我國(guó)大陸對(duì)中國(guó)臺(tái)灣地區(qū)的監(jiān)測(cè)能力相對(duì)較弱,不同省局確定的地震震級(jí)、發(fā)震時(shí)刻以及發(fā)震位置偏差較大,因此,對(duì)于這種情況應(yīng)給予較寬松的閾值。需注意的是,同一省局的快報(bào)目錄內(nèi)不存在重復(fù)定位地震,因此,對(duì)于同一省局的快報(bào)目錄內(nèi)的2次地震,即使?jié)M足以上標(biāo)準(zhǔn),也不認(rèn)為是重復(fù)定位地震。另外,針對(duì)以下特殊情況給出了相應(yīng)的解決方法：

(1)重復(fù)出現(xiàn)的地震未被本省記錄而被其他多個(gè)省局記錄，則保留發(fā)震時(shí)刻最早的記錄。

(2)重復(fù)出現(xiàn)的地震發(fā)生區(qū)域(如海域)沒有明確省局歸屬，則保留記錄到該區(qū)域地震數(shù)最多的省快報(bào)目錄,如果各省快報(bào)記錄的地震數(shù)相同，則保留發(fā)震時(shí)刻最早的記錄。

1.3 EQ3地震目錄格式轉(zhuǎn)換

為將下載得到的最新地震目錄轉(zhuǎn)換至EQ3格式,我們對(duì)EQ3格式地震目錄的二進(jìn)制結(jié)構(gòu)進(jìn)行了分析研究。李閩峰等(2006)指出,EQ3格式地震目錄由2個(gè)相互關(guān)聯(lián)的二進(jìn)制結(jié)構(gòu)來表示記錄,這能夠有效地控制數(shù)據(jù)冗余,并為應(yīng)用程序節(jié)省空間資源提供了有效的機(jī)制。因此,完整的EQ3格式地震目錄包含EQ3和EQB 2個(gè)文件,分別記錄了主要信息和輔助參數(shù)(參考地名和地震震源機(jī)制解參數(shù))。

EQ3主記錄結(jié)構(gòu)定義如下：

typedef struct{

long Date; //地震日期(年*10000+月*100+日)

long Time; //地震時(shí)間(時(shí)*1000000+分*10000+秒*100)

float Lat; //震中緯度(單位:°)

float Lon; //震中經(jīng)度(單位:°)

charMS; //震級(jí)(MS*10)

charML; //震級(jí)(ML*10)

charMb; //震級(jí)(Mb*10)

charMW; //震級(jí)(MW*10)

float Depth; //震源深度(單位:km)

long Sequence; //地震序列編號(hào)

long Index; //對(duì)應(yīng)輔助記錄位置的偏離記錄數(shù)

}EQ3AUX;//EQB格式

根據(jù)EQ3格式的結(jié)構(gòu)定義可知,震中經(jīng)緯度采用Float類型,其參數(shù)值精度保持在7位有效數(shù)值的精度范圍內(nèi),這樣經(jīng)緯度的有效精度就可以表示到小數(shù)點(diǎn)后4位,即相當(dāng)于10m級(jí)(李閩鋒等,2006)。因此，EQ3格式地震目錄是一種能夠有效支持高精度定位結(jié)果的地震目錄格式。EQ3格式中每個(gè)地震記錄主記錄結(jié)構(gòu)的最后部分記錄的信息為輔助記錄索引Index,它表示對(duì)應(yīng)于該地震記錄的輔助參數(shù)在輔助記錄文件中的順序號(hào)(李閩鋒等,2006)。目前EQ3格式地震目錄采用的輔助記錄定義為EQB格式,參考李閩鋒等(2006)給出的EQA格式的二進(jìn)制結(jié)構(gòu),EQB格式的輔助記錄結(jié)構(gòu)定義如下：

typedef struct{

char Location[32]; //地震參考地名(16個(gè)漢字)

float StrikeI; //節(jié)面I走向

float DipI; //節(jié)面I傾角

float StrikeII; //節(jié)面Ⅱ走向

float DipII; //節(jié)面Ⅱ傾角

float Pazim; //P軸方位角

float Pdip; //P軸傾角

float Tazim; //T軸方位角

float Tdip; //T軸傾角

float RakeI; //節(jié)面I滑動(dòng)角

float RakeII; //節(jié)面Ⅱ滑動(dòng)角

}EQ3AUX;//EQA格式

這種相互關(guān)聯(lián)的結(jié)構(gòu),使得EQ3格式既有效地描述了較多的地震參數(shù),同時(shí)也將數(shù)據(jù)冗余降到最低。

下載使用的地震目錄資料主要包括地震的發(fā)震時(shí)刻、震中經(jīng)緯度、震源深度、震級(jí)及參考地名信息。地震編目系統(tǒng)提供的高精度定位結(jié)果的地震目錄的經(jīng)緯度信息精確到了小數(shù)點(diǎn)后3位,EQ3格式地震目錄能夠有效支持高精度的定位結(jié)果。在格式轉(zhuǎn)換過程中,應(yīng)將參考地名信息存儲(chǔ)于EQB格式的輔助記錄內(nèi),其他信息寫入EQ3格式的主記錄內(nèi),同時(shí)根據(jù)輔助信息所在文件的位置,將相應(yīng)的Index順序號(hào)寫入主記錄。如果沒有對(duì)應(yīng)的輔助參數(shù),應(yīng)將Index賦值為-1。

2 應(yīng)用效果

使用Python的圖形工具包PyQt編寫了程序的可視化界面(圖3),并且預(yù)設(shè)了常用地震目錄的下載參數(shù)。設(shè)置好待更新的地震目錄路徑后,更新地震目錄時(shí)只需選擇目錄類型以及起止時(shí)間即可自動(dòng)完成地震目錄的更新工作,無需其他人工操作。為驗(yàn)證所提方法的有效性,以全國(guó)快報(bào)地震目錄更新一周時(shí)間內(nèi)的地震為例,分別測(cè)試人工手動(dòng)和計(jì)算機(jī)自動(dòng)2種更新方式的用時(shí)。表1 對(duì)比了人工手動(dòng)和計(jì)算機(jī)自動(dòng)更新全國(guó)快報(bào)地震目錄所涉及的操作,可以看出人工手動(dòng)更新方式的操作步驟較為繁瑣。表2 列出了同一測(cè)試人員采用不同更新方式在多次測(cè)試中使用的時(shí)間。從對(duì)比結(jié)果可以看出,本文所提的地震目錄自動(dòng)更新方法能夠有效提高地震目錄更新的工作效率,且耗時(shí)穩(wěn)定,受人為因素影響較小,而人工手動(dòng)更新方式的耗時(shí)受人為因素影響較大。另外,由于操作人員的熟練程度不同,更新地震目錄的耗時(shí)也會(huì)有較大差異。

圖3 地震目錄更新程序界面

表1 人工手動(dòng)更新和計(jì)算機(jī)自動(dòng)更新全國(guó)快報(bào)目錄的操作對(duì)比

表2 全國(guó)快報(bào)地震目錄更新對(duì)比測(cè)試結(jié)果

本文所提方法支持高精度EQ3格式地震目錄的自動(dòng)更新,以全國(guó)統(tǒng)一正式報(bào)地震目錄為例,對(duì)比傳統(tǒng)地震目錄與高精度定位結(jié)果的地震目錄的定位結(jié)果。分別利用2種地震目錄選取了九寨溝余震區(qū)2018年1～12月的ML>0.0的地震,繪制了相應(yīng)的震中分布圖(圖4)。從對(duì)比結(jié)果可以看出,通過高精度定位結(jié)果的地震目錄繪制的震中分布(圖4(b))更符合地質(zhì)構(gòu)造背景。特別是在研究小區(qū)域尺度地震活動(dòng)性時(shí),如果地震定位精度過低,震中分布圖上會(huì)呈現(xiàn)較明顯的網(wǎng)格狀分布,更高精度的地震目錄能夠刻畫更加細(xì)致的地震活動(dòng)分布,揭示斷裂構(gòu)造形態(tài),有利于分析認(rèn)識(shí)區(qū)域地震活動(dòng)特征。

圖4 九寨溝7.0級(jí)地震余震區(qū)2018年ML>0.0地震分布對(duì)比

3 結(jié)語

本文設(shè)計(jì)了一種針對(duì)高定位精度EQ3格式地震目錄自動(dòng)更新的程序,該程序利用Python爬蟲技術(shù),從地震編目網(wǎng)站自動(dòng)下載高精度定位結(jié)果的地震目錄,并轉(zhuǎn)換至EQ3格式。針對(duì)相鄰省局快報(bào)地震目錄內(nèi)存在重復(fù)定位地震的問題,提出了一種自動(dòng)篩除重復(fù)定位地震的方法，并應(yīng)用于地震目錄的自動(dòng)更新中。采用PyQt編寫了程序的可視化界面,同時(shí)預(yù)設(shè)地震目錄下載參數(shù),使該程序具有良好的實(shí)用性。實(shí)際應(yīng)用表明,該程序能夠有效提高地震目錄的更新效率,相比傳統(tǒng)的人工手動(dòng)更新方式,能夠提高近30倍的工作效率。

利用高精度的地震目錄能夠更細(xì)致地刻畫地震分布、揭示斷裂構(gòu)造形態(tài),有助于對(duì)地震活動(dòng)性進(jìn)行深入分析。目前,高精度的地震目錄尚未在地震活動(dòng)性分析工作中充分發(fā)揮作用,本文提出的高精度定位結(jié)果的地震目錄自動(dòng)更新程序?qū)⒂兄趯⒏呔榷ㄎ唤Y(jié)果的地震目錄進(jìn)一步應(yīng)用于地震活動(dòng)性分析中。

致謝：審稿專家對(duì)本文提出了寶貴的建議,中國(guó)地震臺(tái)網(wǎng)中心孟令媛副研究員為本項(xiàng)工作及論文撰寫給予了指導(dǎo)并進(jìn)行了有益的討論,薛艷副研究員對(duì)本項(xiàng)工作提供了有益的建議,本文部分圖件采用GMT繪制,在此一并表示感謝。