陳向軍，劉杰超，滕海濤

(新疆維吾爾自治區(qū)地震局，新疆 烏魯木齊 830011）

0 引言

震情會(huì)商是當(dāng)前預(yù)報(bào)水平下形成地震預(yù)報(bào)意見的一個(gè)重要環(huán)節(jié)[1]。做好震情會(huì)商工作勢(shì)必需要分析處理大量的前兆和地震目錄數(shù)據(jù)。目前，基于GIS的地震分析預(yù)報(bào)系統(tǒng)MAPSIS，是一個(gè)全新的、在Windows平臺(tái)上，整合Mapinfo6.X的軟件系統(tǒng)。它技術(shù)新，功能齊全，方便實(shí)用，是分析預(yù)報(bào)人員的好幫手[2]，MAPSIS軟件用到的數(shù)據(jù)格式主要有TXT格式、EQT格式、WKF格式，而對(duì)來自地震目錄周報(bào)、月報(bào)、年報(bào)的DBF文件不能直接應(yīng)用。地震尾波QC值是描述地球介質(zhì)非均勻程度的確定性指標(biāo)，其大小及對(duì)頻率的依賴程度反映了介質(zhì)的非均勻性和非彈性特征[3]。QC值研究作為地震學(xué)的一項(xiàng)研究方法，其前期篩選符合條件地震事件的工作量較大。從事地震分析預(yù)報(bào)，用程序求解震源機(jī)制是經(jīng)常性的工作，原始資料的數(shù)據(jù)格式需要進(jìn)行相應(yīng)轉(zhuǎn)換才能使用。為了方便地震分析預(yù)報(bào)人員開展研究，研制了一套地震日常分析會(huì)商輔助軟件，實(shí)現(xiàn)地震事件篩選和相應(yīng)的數(shù)據(jù)格式轉(zhuǎn)換，其中DBF表數(shù)據(jù)轉(zhuǎn)換為EQT文件格式、QC值篩選地震、單臺(tái)小震綜合機(jī)制解數(shù)據(jù)格式轉(zhuǎn)換已在新源地震臺(tái)用于會(huì)商和地震研究工作中。

1 軟件的功能特點(diǎn)

Delphi是一套集成開發(fā)環(huán)境的程序語言開發(fā)軟件，它提供了包括設(shè)計(jì)、開發(fā)、測(cè)試、調(diào)試及部署應(yīng)用程序所需的全部程序工具，讓設(shè)計(jì)器很容易根據(jù)自己的需求，開發(fā)出實(shí)用的應(yīng)用程序（Delphi7完美經(jīng)典，2003）。

地震日常分析會(huì)商工作中，需要將地震目錄DBF月報(bào)數(shù)據(jù)文件格式轉(zhuǎn)換成EQT文件，日常工作多使用CAPSeis2.0或升級(jí)版的Super Seis軟件，為了減少用戶跳轉(zhuǎn)使用多個(gè)軟件，本軟件在功能上增加了DBF文件格式轉(zhuǎn)換，便于用戶操作使用。QC值篩選地震，以往主要依靠人工瀏覽篩選，效率低且容易遺漏地震事件。單臺(tái)小震綜合機(jī)制解數(shù)據(jù)格式轉(zhuǎn)換，以往采用人工排版編輯，隨著地震事件數(shù)量的增大，人工編輯操作更加困難。地震日常分析會(huì)商輔助軟件針對(duì)以上3個(gè)問題，優(yōu)化了數(shù)據(jù)處理流程，實(shí)現(xiàn)窗口化操作界面，提高了工作效率。

軟件主要包括3個(gè)模塊，分別是：DBF數(shù)據(jù)表文件轉(zhuǎn)換為EQT文件格式、QC值篩選地震事件、單臺(tái)小震綜合機(jī)制解數(shù)據(jù)格式轉(zhuǎn)換（圖1）。

圖1 地震日常分析會(huì)商輔助軟件構(gòu)成Fig.1 Structure of the assistant software of earthquake routine analysis and consultation

2 軟件的設(shè)計(jì)與實(shí)現(xiàn)

軟件采用模塊化設(shè)計(jì)，3個(gè)模塊之間相互獨(dú)立。在軟件設(shè)計(jì)過程中，為了便于用戶操作，將常用參數(shù)設(shè)置放在軟件主界面上。最終實(shí)現(xiàn)DBF數(shù)據(jù)表文件轉(zhuǎn)換為EQT格式；QC值篩選地震事件,主要依據(jù)臺(tái)站的月報(bào)文件從地震事件保存目錄中篩選出符合條件的地震事件EVT文件；單臺(tái)小震綜合機(jī)制解數(shù)據(jù)格式轉(zhuǎn)換，按震源機(jī)制求解程序的要求，將單臺(tái)記錄的小震P波初動(dòng)按一定格式排列，并將實(shí)測(cè)資料轉(zhuǎn)換為相應(yīng)格式[4～6]。

地震目錄DBF數(shù)據(jù)文件是從地震目錄數(shù)據(jù)庫服務(wù)器中每月定時(shí)導(dǎo)出并在地震系統(tǒng)內(nèi)部網(wǎng)站上發(fā)布（見表1），為專業(yè)技術(shù)人員開展地震研究工作提供數(shù)據(jù)資料。本軟件利用Delphi的數(shù)據(jù)庫技術(shù)，實(shí)現(xiàn)讀取地震目錄EXCEL數(shù)據(jù)文件（來自地震目錄DBF文件的轉(zhuǎn)換）的字段內(nèi)容，提取地震事件主要參數(shù)（發(fā)震時(shí)刻、震中經(jīng)緯度、深度、震級(jí)），以EQT文件格式輸出。

在Delphi下讀取EXCEL文件使用的主要控件有：ExcelApplication，ExcelWorksheet，ExcelOLEObject，ExcelWorkbook。具體代碼如下：

EQT文件格式為：年 月 日 時(shí)分 秒 緯度 經(jīng)度震級(jí) 深度。例如，1956 3 8 44950 34.35 104.115.00 110

地震波傳播品質(zhì)因子QC值是量度介質(zhì)衰減的基本物理參數(shù)之一，品質(zhì)因子與區(qū)域構(gòu)造活動(dòng)、地震活動(dòng)密切相關(guān)，可以作為評(píng)價(jià)區(qū)域構(gòu)造活動(dòng)的基本指標(biāo)，用以分析區(qū)域地震活動(dòng)水平[7]。采用朱新運(yùn)等研制的基于Sato模型和Matlab語言的近震S波尾波QC值求解及分析軟件[8]，在進(jìn)行地震事件資料選取時(shí)，根據(jù)震級(jí)、震中距等條件挑選地震波形信噪比高，事件波形文件截取時(shí)間足夠長，震相清晰的地震事件。在成百上千的地震事件中挑選符合條件的地震事件，其工作量非常大。根據(jù)臺(tái)站月報(bào)文件的處理結(jié)果，結(jié)合地震事件EVT文件的命名規(guī)則，通過軟件快速將符合條件的地震事件EVT文件挑選到指定目錄下，提高了研究工作的效率。在Delphi下以地震目錄月報(bào)文件（y.dat，TXT格式）內(nèi)容為依據(jù)，將所有符合條件的目錄事件編輯生成一個(gè)DOS批處理文件，并最后自動(dòng)執(zhí)行該批處理文件，完成地震事件的篩選工作。

表1 地震目錄DBF數(shù)據(jù)文件格式Table 1 Data Format of earthquake catalog

y.dat地震目錄月報(bào)文件內(nèi)容如下：

DOS批處理文件內(nèi)容如下：

以上內(nèi)容，必須提前對(duì)EVT事件目錄和挑選結(jié)果的保存路徑進(jìn)行指定，本軟件已配有相應(yīng)功能，可靈活修改。

小震綜合震源機(jī)制解可以反映局部地區(qū)應(yīng)力狀態(tài)，求解震源機(jī)制的程序使用固定格式的數(shù)據(jù)文件，將單臺(tái)記錄的小震P波初動(dòng)按要求進(jìn)行人工編輯排列費(fèi)時(shí)費(fèi)力，且容易出錯(cuò)，隨著研究區(qū)域的擴(kuò)大，時(shí)間范圍的加長，大量的基礎(chǔ)資料已是人工編輯無法實(shí)現(xiàn)的。為了提高多個(gè)單臺(tái)小震綜合震源機(jī)制解研究的工作效率，在Delphi下以讀取文本文件和截取字符串技術(shù)為主，對(duì)單臺(tái)小震初動(dòng)符號(hào)原始文件進(jìn)行循環(huán)遍歷，截取字符串并重新排序結(jié)合，最終生成sta和eth文件，完成相應(yīng)格式的轉(zhuǎn)換。

單臺(tái)小震初動(dòng)符號(hào)原始文件數(shù)據(jù)格式為：年月日時(shí)分秒初動(dòng)符號(hào)（+表示向上，-表示向下）震中緯度經(jīng)度震級(jí)震源深度。

例如，2010122345+022.220055.570.80018000

求解震源機(jī)制的程序需要兩個(gè)基礎(chǔ)數(shù)據(jù)文件sta和eth。

sta文件的數(shù)據(jù)格式為：地震事件編號(hào)，度 分 秒 （地震事件經(jīng)度）度 分 秒 （地震事件緯度）。

例如，AAA,83 34 12 44 25 12

eth文件的數(shù)據(jù)格式為：

首 行 年月日-時(shí)分 例如，20000603-0041

第二行 年 月 日 時(shí) 分 經(jīng)度 (臺(tái)站)緯度 (臺(tái)站)震源深度 震級(jí)

例如，2000 06 03 00 41 83 33 43 33 20 5.0

第三行地震事件編號(hào)，P波初動(dòng)方向

例如，AAA,-1

3 軟件的使用及安裝

3.1 軟件的使用

（1）地震目錄DBF月報(bào)數(shù)據(jù)文件格式轉(zhuǎn)換模塊。先使用Microsoft Visual FoxPro 6.0或更高版本打開地震目錄DBF月報(bào)數(shù)據(jù)文件，使用“文件”菜單中的導(dǎo)出功能將DBF數(shù)據(jù)表文件轉(zhuǎn)換為EXCEL文件并保存。打開本軟件地震目錄DBF月報(bào)數(shù)據(jù)文件格式轉(zhuǎn)換模塊，直接選擇導(dǎo)出的EXCEL表文件，點(diǎn)擊“轉(zhuǎn)換”按鈕，將地震目錄DBF月報(bào)數(shù)據(jù)文件格式轉(zhuǎn)換為EQT 文件（圖 2）。

圖2 震中參數(shù)EQT文件格式轉(zhuǎn)換界面Fig.2 Format conversion interface of epicenter parameters EQT file

（2）QC值篩選地震事件模塊。打開本軟件QC值篩選地震事件模塊，在當(dāng)前位置中選擇或輸入臺(tái)站月報(bào)文件（如，blt1005y.dat）；在源目錄項(xiàng)中設(shè)置地震事件EVT文件的保存目錄；在目標(biāo)目錄項(xiàng)中設(shè)置篩選出的地震事件EVT文件保存目錄路徑；在震級(jí)ml和震中距項(xiàng)中設(shè)置所要挑選地震事件的條件。點(diǎn)擊“篩選地震”按鈕，挑選出符合條件的地震事件EVT 文件（圖 3）。

（3）小震綜合震源機(jī)制解數(shù)據(jù)格式轉(zhuǎn)換模塊。打開本軟件小震綜合震源機(jī)制解數(shù)據(jù)格式轉(zhuǎn)換模塊，在待轉(zhuǎn)換數(shù)據(jù)文件項(xiàng)中選擇單臺(tái)小震初動(dòng)符號(hào)原始文件；在時(shí)間、地殼厚度、震級(jí)MS、單臺(tái)臺(tái)站經(jīng)緯度項(xiàng)中輸入相應(yīng)參數(shù)；在保存文件名中輸入轉(zhuǎn)換后sta和eth兩個(gè)文件名的前綴；在結(jié)果保存地址項(xiàng)中輸入轉(zhuǎn)換后的sta和eth兩個(gè)文件的保存目錄。點(diǎn)擊“轉(zhuǎn)換”按鈕，生成求解震源機(jī)制軟件可以直接使用的兩個(gè)基礎(chǔ)數(shù)據(jù)文件sta和eth（圖4）。

圖3 Qc值篩選地震事件界面Fig.3 Interface of selecting earthquake-event by counting Qc Value

圖4 單臺(tái)小震綜合機(jī)制解數(shù)據(jù)格式轉(zhuǎn)換Fig.4 Interface of data format conversion of the composite mechanism solution of small earthquake of single station

3.2 軟件的安裝

地震日常分析會(huì)商輔助軟件屬于綠色軟件，在Windows操作系統(tǒng)（如Win2000、Windows XP）下，直接將本軟件拷貝至任意文件夾下即可使用。

4 結(jié)束語

地震日常分析會(huì)商輔助軟件利用Delphi集成軟件開發(fā)環(huán)境，解決地震日常分析會(huì)商工作中要將地震目錄DBF月報(bào)數(shù)據(jù)文件格式轉(zhuǎn)換成EQT文件、QC值篩選地震事件、小震綜合震源機(jī)制解數(shù)據(jù)格式轉(zhuǎn)換的3個(gè)問題，意在提高地震日常分析會(huì)商工作效率，簡化數(shù)據(jù)處理流程和步驟，實(shí)現(xiàn)相關(guān)功能并使操作簡單易行。本軟件的研制有助于專業(yè)技術(shù)人員進(jìn)行地震日常分析會(huì)商工作，對(duì)于促進(jìn)專業(yè)技術(shù)人員開展地震研究發(fā)揮了輔助作用。

[1]王建國，聶永安，董洪軍，等.地震日常分析會(huì)商數(shù)據(jù)庫管理軟件[J].西北地震學(xué)報(bào)，2008，30(1)：94-95.

[2]馬寶君，馬莉，崔東源，等.地震目錄快速轉(zhuǎn)換軟件[J].地震地磁觀測(cè)與研究，2005，26(3)：87-88.

[3]朱新運(yùn)，張帆.浙江及鄰區(qū)尾波Q值研究[J].西北地震學(xué)報(bào)，2006，28(2)：108-113.

[4]楊志國，王小瓊，李世姣，編著.Delphi應(yīng)用程序開發(fā)精彩實(shí)例[M].北京:機(jī)械工業(yè)出版社，2003

[5]梁水，李方超，賽奎春.Delphi開發(fā)技術(shù)大全[M].北京：人民郵電出版社，2007.

[6]張?jiān)鰪?qiáng)編著.突破Delphi7.0編程實(shí)例五十講[M].北京：中國水利水電出版社，2002

[7]王日風(fēng)，張秀萍，韓和平，等.基于Aki模型的張北地震序列QC值特征研究[J].防災(zāi)科技學(xué)院學(xué)報(bào)，2010， 12(2)： 37-38.

[8]魏紅梅，賀曼秋，黃世源，等.重慶榮昌地區(qū)尾波QC值特征[J].西北地震學(xué)報(bào)，2009，31(1)：99-100.