Matlab在天津市強(qiáng)震動(dòng)臺(tái)網(wǎng)烈度速報(bào)中的應(yīng)用①

劉雙慶，聶永安，高武平，李雅靜

（天津市地震局，天津 300201）

0 引言

烈度速報(bào)是一項(xiàng)非常重要的地震災(zāi)害快速評(píng)估工作，快速產(chǎn)出準(zhǔn)確的烈度評(píng)估圖可以有效地指導(dǎo)地震救援，從而極大地減輕地震造成的二次災(zāi)害。隨著我國(guó)強(qiáng)震動(dòng)臺(tái)網(wǎng)的建設(shè)與完善，烈度速報(bào)工作已獲得了長(zhǎng)足的進(jìn)展，并在近幾年向地震預(yù)警工作方面進(jìn)行開(kāi)拓［1］。“九五”期間，中國(guó)地震局建立了我國(guó)第一個(gè)由72個(gè)實(shí)時(shí)數(shù)據(jù)傳輸強(qiáng)震臺(tái)和80個(gè)電話撥號(hào)數(shù)據(jù)傳輸強(qiáng)震臺(tái)構(gòu)成的地震動(dòng)強(qiáng)度（烈度）速報(bào)臺(tái)網(wǎng)。在“十五”期間，中國(guó)地震局除了繼續(xù)對(duì)北京市、天津市的地震烈度速報(bào)臺(tái)網(wǎng)進(jìn)行加密建設(shè)外，還在昆明、蘭州、烏魯木齊等城市新建地震烈度速報(bào)臺(tái)網(wǎng)，部分省市地震局也正在計(jì)劃建設(shè)具有地震烈度速報(bào)功能的強(qiáng)震臺(tái)網(wǎng)［2］。烈度速報(bào)的形式也從過(guò)去的宏觀調(diào)查逐步向儀器烈度側(cè)重［3－4］，發(fā)揮儀器烈度在震害快速評(píng)估中的作用。烈度等震線也加入了場(chǎng)地因素并對(duì)不同烈度段的PGA、PGV、PGD進(jìn)行了融合，這些技術(shù)方法在美國(guó)的Shakemap和日本的UrEDAS都有體現(xiàn)［5－7］。澤仁志瑪還就等震線的顯示問(wèn)題進(jìn)行了討論［8］，劉雙慶［9］則利用蒙特卡羅法分析了強(qiáng)震臺(tái)網(wǎng)布局對(duì)等震線的繪制精度的影響。

隨著強(qiáng)震動(dòng)臺(tái)網(wǎng)的高密度建設(shè)，不同廠家不同型號(hào)的強(qiáng)震儀器得到了使用，從而產(chǎn)生了多種數(shù)據(jù)回收方式及數(shù)據(jù)保存格式。如何快速地處理這些不同型號(hào)不同存儲(chǔ)格式的大批量強(qiáng)震數(shù)據(jù)成為了快速烈度速報(bào)的阻礙，就目前天津市強(qiáng)震動(dòng)臺(tái)網(wǎng)110個(gè)強(qiáng)震臺(tái)而言，就涵蓋了四種類(lèi)型的儀器。其部分回傳的數(shù)據(jù)還存在著壞死而無(wú)法用廠家提供的軟件打開(kāi)的問(wèn)題，且直接利用廠家提供的軟件無(wú)法批量處理強(qiáng)震數(shù)據(jù)并快速歸算烈度值。因此本文擬利用已成為國(guó)際控制界的標(biāo)準(zhǔn)計(jì)算軟件Matlab來(lái)解決上述的主要難題，并對(duì)其中的核心處理技巧進(jìn)行說(shuō)明。

1 Matlab及天津市強(qiáng)震動(dòng)臺(tái)網(wǎng)簡(jiǎn)介

1．1 Matlab簡(jiǎn)介

20世紀(jì)70年代，美國(guó)新墨西哥大學(xué)Cleve Moler為了減輕學(xué)生編程的負(fù)擔(dān)，用FORTRAN編寫(xiě)了最早的Matlab。1984年由Little、Moler、Steve Bangert合作成立的MathWorks公司正式把Matlab推向市場(chǎng)。到20世紀(jì)90年代，Matlab已成為國(guó)際控制界的標(biāo)準(zhǔn)計(jì)算軟件，并在2011年推出了Matlab2011b（7．13版本）。從 Matlab6．5開(kāi)始，其主要代碼已改由Java語(yǔ)言編寫(xiě)，并與Fortran，C，C＋＋，Java，Vb等語(yǔ)言兼容。由于 MathWorks公司允許世界各地的專(zhuān)家將自己開(kāi)發(fā)的Matlab子程序上傳到該公司網(wǎng)站，并由MathWorks公司的工程師對(duì)程序代碼進(jìn)行復(fù)核，通過(guò)全面測(cè)試后在新版的Matlab中追加上去，因此目前的Matlab涵蓋了多個(gè)模塊的應(yīng)用程序，主要包括：Web服務(wù)、信息通信、財(cái)政金融、控制系統(tǒng)設(shè)計(jì)、曲線擬合、數(shù)據(jù)庫(kù)使用、濾波器設(shè)計(jì)、模糊邏輯分析、仿生算法、圖像處理、地圖工具、神經(jīng)網(wǎng)絡(luò)、Mu－分析、優(yōu)化分析、Robust控制、信號(hào)處理、樣條工具箱、統(tǒng)計(jì)工具箱、符號(hào)數(shù)學(xué)、小波分析、偏微分方程分析、線性矩陣不等式控制工具箱等多個(gè)模塊，并有良好的界面操作功能和非常優(yōu)秀的使用指導(dǎo)說(shuō)明。

1．2 天津市強(qiáng)震動(dòng)臺(tái)網(wǎng)簡(jiǎn)介

天津市強(qiáng)震臺(tái)網(wǎng)是在“九五”實(shí)時(shí)數(shù)據(jù)傳輸強(qiáng)震臺(tái)的基礎(chǔ)上，在“十五”期間重新建設(shè)起來(lái)的，總共110個(gè)強(qiáng)震臺(tái)（工程力學(xué)研究所的SLJ加速度儀暫停用），涵蓋了美國(guó)凱尼、瑞士GeoSig、SysCom、港震GSMA－2400IP四種類(lèi)型的儀器，平均間距14．5 km。其中有31個(gè)臺(tái)（光纖通訊）具有準(zhǔn)實(shí)時(shí)波形回傳的條件，另外79臺(tái)當(dāng)觸發(fā)條件滿足后能在15分鐘內(nèi)自動(dòng)回傳數(shù)據(jù)到天津市地震局?jǐn)?shù)據(jù)處理中心數(shù)據(jù)庫(kù)。基于Oracle數(shù)據(jù)庫(kù)、Java運(yùn)行環(huán)境及 Web管理的強(qiáng)震動(dòng)臺(tái)網(wǎng)實(shí)現(xiàn)了強(qiáng)震動(dòng)數(shù)據(jù)自動(dòng)匯集、存儲(chǔ)，準(zhǔn)實(shí)時(shí)或定時(shí)監(jiān)控儀器運(yùn)行狀態(tài)和臺(tái)站遠(yuǎn)程控制與管理等功能。強(qiáng)震臺(tái)網(wǎng)運(yùn)行以來(lái)，獲得了2008年5月汶川M8．0大地震較好的加速度記錄和2010年4月9日豐南M4．1地震加速度記錄。但此前烈度速報(bào)效率低下的情況比較嚴(yán)重。

2 天津市強(qiáng)震動(dòng)臺(tái)網(wǎng)烈度速報(bào)的幾個(gè)技術(shù)難點(diǎn)及Matlab的應(yīng)用

通過(guò)分析研究，歸納天津市強(qiáng)震動(dòng)臺(tái)網(wǎng)速報(bào)的技術(shù)難點(diǎn)如下：

（1）數(shù)據(jù)格式不統(tǒng)一，需要廠家的界面化軟件對(duì)事件文件進(jìn)行一個(gè)一個(gè)人工交互式處理（除個(gè)別廠家儀器數(shù)據(jù)能批處理外）。

（2）數(shù)據(jù)文件內(nèi)容存在壞死情況，軟件無(wú)法打開(kāi)，用解碼程序運(yùn)行則顯示錯(cuò)誤。

（3）數(shù)據(jù)文件命名不規(guī)范，數(shù)據(jù)處理后存放路徑不統(tǒng)一。

（4）用 Web界面調(diào)用Oracle數(shù)據(jù)，操作不便利，影響下載速度。

（5）數(shù)據(jù)預(yù)處理操作及峰值提取無(wú)法批處理。

（6）等值線繪制方法單一且沒(méi)有融合理論烈度。

以上過(guò)程嚴(yán)重影響了烈度速報(bào)的速度和效率。針對(duì)以上難點(diǎn)，本文利用Matlab處理了下面幾個(gè)問(wèn)題：

（1）Matlab對(duì)Oracle數(shù)據(jù)庫(kù)的直接調(diào)用

這個(gè)部分主要解決Web調(diào)用的數(shù)據(jù)不能快速下載的缺陷。由于對(duì)Oracle數(shù)據(jù)文件命名作了新的設(shè)定（圖1），將觸發(fā)事件和儀器標(biāo)定事件的命名做了區(qū)分，因此可以通過(guò)時(shí)間段和事件命名的差異直接提取出Oracle某個(gè)時(shí)間段的地震觸發(fā)事件文件或標(biāo)定文件。其涉及的Matlab代碼有：

當(dāng)在本地機(jī)器安裝好Oracle客戶服務(wù)端后，配置Matlab對(duì)Oracle的驅(qū)動(dòng)，通過(guò)以上代碼就可以直接從遠(yuǎn)程O(píng)racle數(shù)據(jù)庫(kù)批量讀取數(shù)據(jù)文件。

（2）Matlab與exe文件的融合

該部分主要是充分利用廠家提供的exe格式的解碼程序，對(duì)從Oracle數(shù)據(jù)庫(kù)讀出的數(shù)據(jù)文件進(jìn)行批量轉(zhuǎn)換成ascii碼。由于臺(tái)站代碼和命名固定，所以直接通過(guò)文件名就可以確定采用哪個(gè)廠家的exe文件進(jìn)行數(shù)據(jù)轉(zhuǎn)換。命名方法見(jiàn)圖1，如果是標(biāo)定文件則在上述文件名的基礎(chǔ)上加前綴FT＿，即function test的縮寫(xiě)。

Matlab與exe的調(diào)用可以通過(guò)Matlab提供的system或dos兩個(gè)命令來(lái)操作．由于存在數(shù)據(jù)文件內(nèi)容壞死而無(wú)法解碼的情況，因此可以進(jìn)一步使用Matlab提供的try．．．catch．．．end的命令來(lái)嘗試解碼，當(dāng)解碼不成功時(shí)直接跳過(guò)，從而保證后續(xù)程序的連續(xù)運(yùn)行。

圖1 強(qiáng)震事件文件命名規(guī)則Fig．1 Naming rule of strong motion record file．

（3）Matlab的批量顯示與通道選擇

這部分可通過(guò)Matlab界面化函數(shù)uicontrol（gcf，‘style’，‘slider’，‘unit’，‘normalized’，‘position’），uicontrol（gcf，‘style’，‘listbox’，‘Position’）組合的方式，提供一個(gè)左側(cè)的滑動(dòng)選擇框及右側(cè)的圖形曲線顯示框（圖5）來(lái)確認(rèn)哪道數(shù)據(jù)可以使用。

（4）儀器烈度的換算

該部分先對(duì)第（3）部分選擇出來(lái)的數(shù)據(jù)道保存，取每道觸發(fā)點(diǎn)前20s的數(shù)據(jù)均值作為基線，以去除各道基線漂移。然后利用ButterWorth帶通濾波器壓制低頻和高頻的部分波形能量，考慮強(qiáng)震儀器的寬頻與高通特性以及建筑物反應(yīng)譜設(shè)計(jì)曲線［10］，本文通頻帶設(shè)為20Hz～2s。之后程序?qū)⑻崛∶總€(gè)臺(tái)站兩個(gè)水平道上的第20個(gè)最大點(diǎn)作為有效峰值，按每套儀器滿量程時(shí)對(duì)應(yīng)的gal值，計(jì)算出每個(gè)儀器兩水平道第20個(gè)最大點(diǎn)對(duì)應(yīng)的gal值的均值。提取絕對(duì)值排序后的第20個(gè)最大值作為峰值，可以有效避免儀器尖脈沖的影響（單道記錄有時(shí)存在多個(gè)尖脈沖干擾），同時(shí)因?yàn)椴蓸勇蕿?00sps，不會(huì)嚴(yán)重漏用記錄峰值及其附近區(qū)數(shù)據(jù)。程序經(jīng)多次測(cè)試，該選點(diǎn)一般可取第5～25點(diǎn)，換算出的gal值較穩(wěn)定。然后按中國(guó)地震烈度表 GB／T 17742－2008［11］中烈度與gal的數(shù)值表，采用線性內(nèi)插的方式（Matlab命令：interp1）獲取gal對(duì)應(yīng)的儀器烈度值。當(dāng)然interp1命令還提供了最近鄰、三次多項(xiàng)式、樣條法等多種插值方法。由于烈度表不提供V以下的PGA參考值，因此這些外推的低烈度區(qū)儀器等值線僅作參考，且其本身因烈度低，其災(zāi)評(píng)意義亦不明顯（Ⅵ以上才有明顯地震災(zāi)情）。另外，本文亦嘗試?yán)妹绹?guó)ATC－3規(guī)范中的加速度有效峰值反應(yīng)譜EPA（線性化方法計(jì)算程序見(jiàn)附錄）轉(zhuǎn)換儀器烈度，但由于地震樣本量太少，尚未能深入分析本地區(qū)PGA，EPA的具體差異。

（5）Matlab的儀器烈度顯示

對(duì)于烈度的顯示，除了利用Delaunay三角剖分三角形顯示之外［9］，本部分還提供了兩種繪制方案．一是利用Matlab的scatter命令來(lái)對(duì)不同空間點(diǎn)上的烈度值按數(shù)據(jù)大小決定顯示點(diǎn)大小及顏色的離散方式顯示（圖7（a）），這種方式比等值線法更能體現(xiàn)單點(diǎn)的儀器烈度。二是利用華北地區(qū)理論烈度衰減模型（需要人工輸入地震震中的經(jīng)緯度、震級(jí)和長(zhǎng)軸方向這些參數(shù)，圖4（c））構(gòu)建大網(wǎng)格點(diǎn)（本地區(qū)取0．2°）的理論烈度場(chǎng)，然后將儀器記錄烈度值加密并取代附近的理論烈度，之后利用griddata命令對(duì)非規(guī)則的加密后的數(shù)據(jù)插值，生成烈度等值線。這種方法既利用了儀器烈度也使用了理論烈度，是一種綜合效果（圖7（b））。而地理底圖的顯示則直接利用Matlab與arcgis地圖格式兼容的命令shaperead和mapshow進(jìn)行地理圖讀取與顯示。當(dāng)然，低烈度區(qū)（IV以下）的理論計(jì)算值亦僅作為畫(huà)圖參考，其本身災(zāi)評(píng)意義較弱。

（6）可移植性問(wèn)題

考慮程序的可移植性，程序使用過(guò)程中所涉及的程序路徑都是相對(duì)路徑，不存在放在不同路徑無(wú)法運(yùn)行的問(wèn)題。其次，對(duì)于臺(tái)站名命名、地理底圖、廠家提供的解碼程序等靜態(tài)資料，程序統(tǒng)一從一個(gè)cfg文件夾獲取，cfg隨整個(gè)程序綁定，只要修改cfg的內(nèi)容并按圖1的命名方法對(duì)數(shù)據(jù)文件命名，且替換相應(yīng)的地理底圖，則本文編寫(xiě)的程序可以移植到其他省臺(tái)網(wǎng)使用。同時(shí)在最后烈度成圖的操作上本文采用了調(diào)用計(jì)算過(guò)程中生成的文件內(nèi)容方式，因此當(dāng)有野外實(shí)地烈度調(diào)查結(jié)果時(shí)可以補(bǔ)充到這個(gè)烈度結(jié)果文件中，從而進(jìn)一步突出烈度等值線對(duì)災(zāi)情現(xiàn)場(chǎng)烈度的反映。

利用Matlab人工處理天津市強(qiáng)震臺(tái)網(wǎng)數(shù)據(jù)與烈度顯示流程及使用的Matlab主要命令如圖2所示。

（7）“事件偵聽(tīng)”

由于地震是突發(fā)事件，為了保證烈度速報(bào)的產(chǎn)出速度，本文的Matlab程序除了提供人機(jī)交互處理的方式外，還提供了自動(dòng)觸發(fā)處理方式——“事件偵聽(tīng)”（圖4）。這是建立在天津市強(qiáng)震動(dòng)臺(tái)網(wǎng)數(shù)據(jù)上傳的具體特征的基礎(chǔ)上。一般地，地震到來(lái)后在5～13分鐘內(nèi)觸發(fā)的臺(tái)站基本把數(shù)據(jù)自動(dòng)上傳到了中心數(shù)據(jù)庫(kù)。而無(wú)地震的情況下，上傳到數(shù)據(jù)庫(kù)的文件在5～13分鐘內(nèi)一般小于2個(gè)，因此可以每5分鐘對(duì)數(shù)據(jù)庫(kù)上傳事件的個(gè)數(shù)進(jìn)行統(tǒng)計(jì)，如果超過(guò)一定的限度，則認(rèn)為有叢集事件記錄上傳，并按如下4個(gè)步驟進(jìn)行自動(dòng)處理（圖3）。

圖2 Matlab處理天津市強(qiáng)震動(dòng)臺(tái)網(wǎng)數(shù)據(jù)流程圖Fig．2 Flow chart of dealing with the data of Tianjin strong motion network using Matlab．

圖3 自動(dòng)掃描處理示意圖Fig．3 Sketch of auto－scaning and auto－processing．

①如果有地震事件發(fā)生，下次循環(huán)時(shí)間由5分鐘改為4．5分鐘并等待8分鐘，預(yù)留0．5分鐘提供程序自動(dòng)處理前14分鐘的數(shù)據(jù)，并且提取滿足掃描條件時(shí)刻的前1分鐘的數(shù)據(jù)，因此總共處理8＋1＋5（其中5＝4．5＋0．5）分鐘內(nèi)上傳到數(shù)據(jù)庫(kù)里的數(shù)據(jù)。當(dāng)沒(méi)有重復(fù)事件在13（8＋5）分鐘內(nèi)發(fā)生時(shí)，數(shù)據(jù)將提取14分鐘的內(nèi)容，與上一個(gè)5分鐘的循環(huán)有1分鐘的重疊。目前我們?cè)O(shè)定5分鐘內(nèi)上傳數(shù)大于5個(gè)則觸發(fā)，這個(gè)約束條件將可能丟失4個(gè)文件，因?yàn)?分鐘內(nèi)最多可以上傳4個(gè)文件。提前一分鐘的重疊，至少可以挽回2個(gè)文件。

②如果事件發(fā)生后，在13分鐘之后又發(fā)生一次地震事件（即第13～18分鐘再次出現(xiàn)叢集事件上傳，則提取這個(gè)后續(xù)事件的前14分鐘的事件文件附加到上一次事件中，此時(shí)有1．5分鐘的重疊）。總共處理14＋14－1．5＝26．5分鐘內(nèi)上傳的文件數(shù)據(jù)。

③如果在18分鐘后才再次發(fā)生地震事件叢集上傳，則作為第二個(gè)地震事件處理。

④如果不發(fā)生其他事件，則正常進(jìn)入下一個(gè)5分鐘的循環(huán)。

經(jīng)過(guò)一年的實(shí)踐運(yùn)行結(jié)果顯示，未出現(xiàn)誤自動(dòng)處理現(xiàn)象。當(dāng)“設(shè)定5分鐘內(nèi)上傳數(shù)大于2個(gè)”時(shí)，偶爾出現(xiàn)自動(dòng)處理結(jié)果，這些結(jié)果常是由儀器工作不正常，或職工當(dāng)時(shí)進(jìn)行儀器維修造成。

3 應(yīng)用實(shí)例

本實(shí)例采用2010年4月9日豐南M4．1地方震加速度記錄，觸發(fā)臺(tái)站36個(gè)，可利用的數(shù)據(jù)臺(tái)27個(gè)（由于當(dāng)時(shí)外圍網(wǎng)絡(luò)數(shù)據(jù)傳輸程序有缺陷，SYSCOM公司的儀器記錄文件壞死情況較嚴(yán)重）。地震發(fā)生后，當(dāng)時(shí)數(shù)據(jù)處理過(guò)程耗用了多個(gè)小時(shí)。本程序編寫(xiě)完成后，在3分鐘內(nèi)就可以通過(guò)界面化（圖4～6）的鼠標(biāo)操作完成了所有的操作，其中點(diǎn)擊“人工數(shù)據(jù)提取”可以彈出圖4（b）的參數(shù)輸入界面。數(shù)據(jù)繪制的烈度顯示結(jié)果（圖7）也與當(dāng)時(shí)各區(qū)縣群眾打電話詢問(wèn)的頻度很相近，其中寧河縣，塘沽區(qū)和市中心高層住戶有明顯震感。需說(shuō)明的是，圖5、圖6的縱坐標(biāo)做了如下處理：先對(duì)各臺(tái)站記錄數(shù)據(jù)歸一化以壓制異常記錄，然后對(duì)每臺(tái)記錄的基線遞增進(jìn)行顯示以避免曲線堆疊。數(shù)據(jù)提取按內(nèi)存變量值為準(zhǔn)。

圖4 天津強(qiáng)震動(dòng)臺(tái)網(wǎng)速報(bào)軟件界面（a），獲取事件選擇框（b）和理論烈度計(jì)算參數(shù)輸入框（c）Fig．4 GUI of intensity rapid report of Tianjin strong motion（a），GUI of selecting events（b）and GUI of input parameters to calculate theoretical intensity（c）．

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

本文將Matlab軟件提供的調(diào)用數(shù)據(jù)庫(kù)命令，調(diào)用儀器廠家的exe命令，數(shù)據(jù)分析處理命令和界面化操作等命令進(jìn)行有機(jī)組合，從而比較成功地解決了天津市強(qiáng)震動(dòng)臺(tái)網(wǎng)在烈度速報(bào)中遇到的速度與效率問(wèn)題，大大縮減了儀器烈度發(fā)布的時(shí)間，由此得到以下認(rèn)識(shí)：

圖5 Uicontrol函數(shù)實(shí)現(xiàn)的數(shù)據(jù)顯示及通道選擇操作Fig．5 Data imaging and selecting with uicontrol function．

圖6 Uicontrol函數(shù)實(shí)現(xiàn)的數(shù)據(jù)顯示局部放大（由右滑鍵控制）Fig．6Local zooming－in with uicontrol function（Operated by right－side slider）．

（1）Matlab簡(jiǎn)明的程序代碼書(shū)寫(xiě)格式和豐富的命令式程序庫(kù)，可以有效降低應(yīng)用程序開(kāi)發(fā)難度，大大縮減程序代碼量。但Matlab軟件再?gòu)?qiáng)大，也需要和研究者所研究對(duì)象的物理模型、數(shù)學(xué)模型充分融合才能發(fā)揮其作用。本文中涉及的理論烈度計(jì)算的非線性迭代方程，數(shù)據(jù)ButterWorth濾波頻段設(shè)計(jì)，“事件偵聽(tīng)”的循環(huán)判斷方案等等都是具體物理模型的數(shù)學(xué)化，是Matlab本身不具備的。

圖7 豐南4．1級(jí)地震的儀器烈度散點(diǎn)顯示（a）和等值線顯示（b）Fig．7 Equipment intensity imaging by scattering（a）and contouring method（b）for M4．1earthquake at Fennan county．

（2）儀器烈度與傳統(tǒng)烈度的“多因子綜合判定”、“宏觀分級(jí)”、“以結(jié)果表述原因”三個(gè)特征尚有一定的銜接困難，本文的PGA直接轉(zhuǎn)化儀器烈度亦有較多值得商榷之處，利用EPA方法尚還在嘗試，由于缺乏足夠的地震樣本未能深入研究其與PGA的差別。日本目前已不再過(guò)分強(qiáng)調(diào)儀器烈度與傳統(tǒng)烈度的災(zāi)害評(píng)估差異，并減少現(xiàn)場(chǎng)宏觀調(diào)查的工作量，突顯儀器烈度在快速災(zāi)評(píng)的指導(dǎo)，他們這種處理方式不失為一種明智之舉［3，12］。

由于在本文中并未對(duì)天津地區(qū)場(chǎng)地進(jìn)行烈度校正，且沒(méi)有進(jìn)一步討論利用反應(yīng)譜的結(jié)果來(lái)刻畫(huà)儀器烈度的可行性，因此儀器烈度速報(bào)的結(jié)果還是比較粗糙。另外，Matlab的通訊工具箱直接對(duì)串口實(shí)時(shí)數(shù)據(jù)，對(duì)數(shù)字采集器直接訪問(wèn)的研究也未深入展開(kāi)，在人機(jī)交互處理界面上也不太美觀，這些工作都是后續(xù)需要努力的。

致謝：本文程序編寫(xiě)過(guò)程中得到了長(zhǎng)安大學(xué)鄧亞虹博士，北京航空航天大學(xué)吳鵬博士的幫助，在此非常感謝！

［1］ 李山有，金星，馬強(qiáng)，等．地震預(yù)警系統(tǒng)與智能應(yīng)急控制系統(tǒng)研究［J］．世界地震工程，2004，20（4）：21－26．

［2］ 金星，廖詩(shī)榮，陳緋雯．區(qū)域數(shù)字地震臺(tái)網(wǎng)實(shí)時(shí)速報(bào)系統(tǒng)研究［J］．地震地磁觀測(cè)與研究，2007，28（1）：64－72．

［3］ 張敏政．地震烈度及其評(píng)定［J］．防災(zāi)科技學(xué)院學(xué)報(bào)，2010，12（1）：1－6．

［4］ 金星，張紅才，韋永祥．基于地震臺(tái)網(wǎng)資料快速發(fā)布的震動(dòng)烈度標(biāo)準(zhǔn)及其應(yīng)用研究［J］．國(guó)際地震動(dòng)態(tài)，2008，（10）：20－27．

［5］ 李俊，蘇楓，米宏亮，等．ShakeMap及其在地震動(dòng)快速預(yù)估中的應(yīng)用［J］．中國(guó)地震，2010，26（1）：103－111．

［6］ D J Wald，V Quitoriano，T H Heaton，et al．TriNet“ShakeMaps”：Rapid Generation of Peak Ground Motion and Intensity Maps for Earthquakes in Southern California［J］．Earthquake Spectra，1999，15（3）：537－555．

［7］ K T Shabestari，F(xiàn) Yamazaki．A proposal of instrumental seismic intensity scale compatible with MMI evaluated from three－component acceleration records［J］．Earthquake Spectra，2001，17（4）：711－723．

［8］ 澤仁志瑪，陳會(huì)忠，何加勇，等．震動(dòng)圖快速生成系統(tǒng)研究［J］．地球物理學(xué)進(jìn)展，2006，21（3）：809－813．

［9］ 劉雙慶，謝靜，聶永安，等．Delaunay三角剖分算法在地震烈度速報(bào)中的應(yīng)用－－以天津市強(qiáng)震臺(tái)網(wǎng)為例［J］．地震地磁觀測(cè)與研究，2011，32（5）：115－122．

［10］ 國(guó)家標(biāo)準(zhǔn)建筑抗震設(shè)計(jì)規(guī)范管理組 編．建筑抗震設(shè)計(jì)規(guī)范（GB50011－2010）統(tǒng)一培訓(xùn)教材［M］．北京：地震出版社，2010．

［11］ 中國(guó)國(guó)家標(biāo)準(zhǔn)化管理委員會(huì)．中國(guó)地震烈度表GB／T 17742－2008［S］．北京：中國(guó)標(biāo)準(zhǔn)出版社，2009．

［12］ 顏心儀．利用臺(tái)灣寬頻地震網(wǎng)從事強(qiáng)震預(yù)警研究［D］．臺(tái)灣省：臺(tái)灣大學(xué)地質(zhì)科學(xué)研究所，2007．

附錄：Newmark方法計(jì)算地震動(dòng)反應(yīng)譜子程序（Matlab版）：