QGIS鐵路地震影響場(chǎng)自動(dòng)計(jì)算和發(fā)布

孫文韜，王友彪，周學(xué)影

（中國(guó)鐵道科學(xué)研究院集團(tuán)有限公司 鐵道科學(xué)技術(shù)研究發(fā)展中心，北京 100081）

0 引言

地震是影響鐵路行車(chē)安全的自然災(zāi)害之一，為增強(qiáng)我國(guó)高速鐵路防震減災(zāi)能力，自2012年起原鐵道部和中國(guó)地震局聯(lián)合攻關(guān)，開(kāi)始共同研發(fā)高速鐵路地震預(yù)警系統(tǒng)［1-3］，2015年底高速鐵路地震預(yù)警監(jiān)測(cè)系統(tǒng)研發(fā)基本完成，并在福廈、成灌、大西等高速鐵路開(kāi)展了地震預(yù)警監(jiān)測(cè)系統(tǒng)試驗(yàn)，目前已在大西［4］、太焦［5］、京張、京雄等高速鐵路運(yùn)用，并逐步向全國(guó)推廣。

地震強(qiáng)弱程度和范圍是描述地震影響的2個(gè)指標(biāo)，可以通過(guò)地震峰值加速度或烈度等地震動(dòng)參數(shù)的分布反映出來(lái)。根據(jù)《高速鐵路地震預(yù)警監(jiān)測(cè)系統(tǒng)技術(shù)條件》［6］等相關(guān)標(biāo)準(zhǔn)要求，將峰值加速度（Peak Ground Acceleration，PGA）40、80、120 gal作為高速鐵路地震預(yù)警緊急處置Ⅰ、Ⅱ、Ⅲ級(jí)的加速度閾值，故選取鐵路行業(yè)最關(guān)心的PGA作為地震影響場(chǎng)指標(biāo)。

地理信息系統(tǒng)（GIS）在防震減災(zāi)中作用的研究起始于20世紀(jì)80年代末。經(jīng)過(guò)30年的發(fā)展，GIS已在災(zāi)害評(píng)估、地震危險(xiǎn)性分析、地震預(yù)警、地震烈度速報(bào)和緊急救援等領(lǐng)域發(fā)揮了至關(guān)重要的作用［7］。隨著我國(guó)鐵路建設(shè)蓬勃發(fā)展，高速鐵路已四通八達(dá)，鐵路線在空間上呈交錯(cuò)分布，地震發(fā)生的隨機(jī)性導(dǎo)致無(wú)法提前預(yù)估受地震影響的鐵路位置。利用GIS強(qiáng)大的空間分析能力能更好、更快地展現(xiàn)地震的震中和周邊鐵路的分布，再結(jié)合地震動(dòng)傳播衰減規(guī)律可快速估算地震對(duì)周邊鐵路的影響。通過(guò)第一時(shí)間計(jì)算出地震對(duì)鐵路的影響區(qū)段并將影響信息自動(dòng)發(fā)送，可以使鐵路工作者和管理單位最快、最精準(zhǔn)地了解可能受到地震威脅的鐵路，為震后處置和防災(zāi)、減災(zāi)提供數(shù)據(jù)支撐。根據(jù)前期調(diào)研情況，目前鐵路行業(yè)尚未應(yīng)用地震后快速計(jì)算鐵路影響場(chǎng)并實(shí)現(xiàn)信息自動(dòng)發(fā)布的系統(tǒng)。在QGIS（Quantum GIS）軟件中設(shè)計(jì)并開(kāi)發(fā)一個(gè)插件，根據(jù)國(guó)家地震臺(tái)網(wǎng)中心地震目錄，利用鐵路地理信息數(shù)據(jù)和地震動(dòng)衰減關(guān)系實(shí)現(xiàn)地震后快速估算地震對(duì)最近鐵路的影響并自動(dòng)發(fā)布信息，是鐵路地震分析需要進(jìn)行的工作。該工作可為鐵路地震應(yīng)急處置和災(zāi)后救援提供數(shù)據(jù)支撐。

1 QGIS發(fā)展現(xiàn)狀

QGIS是一個(gè)方便使用的開(kāi)源地理信息系統(tǒng)，是地理空間開(kāi)源基金會(huì)（OSGeo）的一個(gè)正式專案軟件，可在Linux、Unix、Mac OS X和Android等系統(tǒng)上應(yīng)用，并支持眾多矢量、柵格、數(shù)據(jù)庫(kù)格式及功能。QGIS目前的測(cè)試版本已經(jīng)更新至3.22，穩(wěn)定版更新至3.16。與成熟但昂貴的商業(yè)軟件ArcGIS不同，QGIS最友好的特點(diǎn)之一就是開(kāi)源，這使其在短短20年中發(fā)展迅速。QGIS是開(kāi)源GIS的集大成者，整合了GRASS、SAGA GIS等多個(gè)開(kāi)源桌面軟件工具，開(kāi)發(fā)者可以自行檢閱與調(diào)整程序代碼，并保障讓所有使用者可以免費(fèi)且自由地修改程序。QGIS擁有自己的開(kāi)源插件社區(qū)，供世界各地的學(xué)者相互應(yīng)用和分享插件［8］。QGIS中各算法之間可以相互調(diào)用，基于QGIS開(kāi)發(fā)鐵路地震加速度影響計(jì)算和發(fā)布程序簡(jiǎn)便易行。

2 國(guó)內(nèi)外研究現(xiàn)狀

地震防災(zāi)、減災(zāi)領(lǐng)域中，國(guó)內(nèi)外相關(guān)學(xué)者利用GIS開(kāi)展了不同方向的研究應(yīng)用。孫哲等［9］基于陜西省不同分區(qū)的橢圓和線源兩種烈度衰減模型，采用Python自動(dòng)化技術(shù)和Arcpy空間繪制技術(shù)自動(dòng)生成地震影響場(chǎng)方法，通過(guò)ArcGIS Server實(shí)現(xiàn)影響場(chǎng)地圖服務(wù)的自動(dòng)發(fā)布；馮駿等［10］基于Arc Engine、C#編程語(yǔ)言開(kāi)發(fā)設(shè)計(jì)了畢節(jié)市地震應(yīng)急分類響應(yīng)研究中地震危險(xiǎn)性模型、建筑物易損性模型及系統(tǒng)，研發(fā)出適合研究區(qū)域震情的地震應(yīng)急分類響應(yīng)系統(tǒng)，該系統(tǒng)實(shí)現(xiàn)了四大功能，包括基礎(chǔ)數(shù)據(jù)管理、地震信息管理、地震知識(shí)庫(kù)管理和震時(shí)評(píng)估；美國(guó)灣區(qū)政府協(xié)會(huì)（ABAG）利用GIS對(duì)舊金山灣區(qū)交通系統(tǒng)開(kāi)展了易損性評(píng)估，并將成果用作該區(qū)域交通規(guī)劃的依據(jù)；日本建筑署建筑研究所在橫濱大地震后，同樣建立基于GIS的防震減災(zāi)系統(tǒng)，開(kāi)展了橫濱大地震的震后損失評(píng)估，利用GIS數(shù)據(jù)庫(kù)統(tǒng)計(jì)了地質(zhì)、斷層、建筑物震害等資料［11］。

3 地震動(dòng)衰減關(guān)系應(yīng)用

地震動(dòng)衰減關(guān)系是建立震級(jí)與加速度關(guān)系的經(jīng)驗(yàn)公式。GB 18306—2015《中國(guó)地震動(dòng)參數(shù)區(qū)劃圖》基于國(guó)內(nèi)外大量強(qiáng)震記錄和美國(guó)NGA（Next Generation Attenuation）的研究成果，同時(shí)利用我國(guó)區(qū)域觀測(cè)數(shù)據(jù)和實(shí)際情況換算統(tǒng)計(jì)，得到適用于我國(guó)的地震動(dòng)衰減關(guān)系模型（見(jiàn)式（1））［12］。

式中：Y為地震動(dòng)參數(shù)PGA，A、B、C、D、E為回歸系數(shù)，M為震級(jí)，R為震中距。

標(biāo)準(zhǔn)中按照地震帶將我國(guó)分為4個(gè)大區(qū)：青藏區(qū)、新疆區(qū)、東部強(qiáng)震區(qū)和中強(qiáng)地震區(qū)，每個(gè)區(qū)域?qū)?yīng)不同的回歸系數(shù)。

根據(jù)地震發(fā)生的震源三要素信息，可以在QGIS中定位震源所在區(qū)域。將震源位置與最近鐵路的距離帶入式（1）R中，可以得到鐵路最近位置的PGA；將PGA=40、80、120 gal分別帶入Y中，也可以得到相應(yīng)的震中距R。以R為半徑圓的范圍，即對(duì)應(yīng)了PGA=40、80、120 gal的影響范圍。

根據(jù)鐵路行業(yè)應(yīng)用的安全余量，衰減關(guān)系通常運(yùn)用長(zhǎng)軸結(jié)果。對(duì)應(yīng)公式中，將表1中不同區(qū)域的長(zhǎng)軸回歸系數(shù)帶入式（1）計(jì)算，可對(duì)應(yīng)獲得該地震的地震影響場(chǎng)。

地震動(dòng)衰減關(guān)系是基于一定樣本量的經(jīng)驗(yàn)推導(dǎo)，估算結(jié)果可能與實(shí)際測(cè)量值有一定偏差。由于我國(guó)絕大多數(shù)地區(qū)并非密集布設(shè)地震臺(tái)站，所以根據(jù)地震動(dòng)衰減關(guān)系預(yù)估我國(guó)任一地區(qū)的地震動(dòng)PGA是目前行業(yè)內(nèi)的通用做法。特定區(qū)域中有符合局部特征更好的衰減關(guān)系，在此不做分析。分區(qū)PGA衰減關(guān)系系數(shù)見(jiàn)表1。

表1 分區(qū)PGA衰減關(guān)系系數(shù)

4 鐵路地震影響場(chǎng)自動(dòng)計(jì)算和發(fā)布

實(shí)現(xiàn)鐵路地震影響場(chǎng)自動(dòng)計(jì)算和發(fā)布需要2個(gè)模塊。一是信息處理模塊，可對(duì)地震目錄進(jìn)行抓取、圖層生成、圖文撰寫(xiě)、整合排版，并推送至短信服務(wù)器；二是地震計(jì)算模塊，可基于信息處理模塊抓取的信息，進(jìn)行數(shù)據(jù)分析處理和數(shù)據(jù)輸出。2個(gè)模塊配合工作，QGIS鐵路地震影響場(chǎng)自動(dòng)計(jì)算和發(fā)布流程見(jiàn)圖1。

圖1 QGIS鐵路地震影響場(chǎng)自動(dòng)計(jì)算和發(fā)布流程

4.1 數(shù)據(jù)處理

利用QGIS實(shí)現(xiàn)鐵路地震PGA影響的計(jì)算主要依靠以下基礎(chǔ)數(shù)據(jù)：

（1）我國(guó)鐵路線分布矢量數(shù)據(jù)。我國(guó)已建成鐵路線的數(shù)據(jù)為開(kāi)源信息，應(yīng)用類地圖軟件（例如高德地圖、百度地圖等）對(duì)于鐵路線的位置均有明確標(biāo)記。研究選取OpenStreet Map上中國(guó)地區(qū)的鐵路線數(shù)據(jù)，進(jìn)行合并矢量、幾何校正、名稱校對(duì)等人工處理工作，作為基礎(chǔ)鐵路線數(shù)據(jù)，開(kāi)源鐵路線數(shù)據(jù)見(jiàn)圖2。

圖2 開(kāi)源鐵路線數(shù)據(jù)

（2）衰減關(guān)系分區(qū)數(shù)據(jù)。如利用第五代區(qū)劃圖中的衰減關(guān)系進(jìn)行運(yùn)算，需按照表1系數(shù)判斷震中所處位置的分區(qū)屬性。相關(guān)信息可查閱GB18306—2015和中國(guó)標(biāo)準(zhǔn)出版社出版的《中國(guó)地震動(dòng)參數(shù)區(qū)劃圖》。衰減關(guān)系分區(qū)數(shù)據(jù)見(jiàn)圖3，其中，深藍(lán)色為東部強(qiáng)震區(qū)、淺藍(lán)色為中強(qiáng)地震區(qū)、粉色為新疆區(qū)、橙色為青藏區(qū)。

圖3 衰減關(guān)系分區(qū)數(shù)據(jù)

如利用其他衰減關(guān)系，可替換表1中的數(shù)據(jù)。

（3）底圖數(shù)據(jù)。底圖數(shù)據(jù)主要為展現(xiàn)地震震源、鐵路線等要素的具體參考位置和周邊興趣點(diǎn)（POI）信息，不參與地震分析和計(jì)算。底圖需要在線或離線加載，研究建議選取天地圖的瓦片服務(wù)形式在線加載，具體可查詢天地圖官網(wǎng)（https：//www.tianditu.gov.cn/）。

以上數(shù)據(jù)為開(kāi)源數(shù)據(jù)，程序的功能實(shí)現(xiàn)不受數(shù)據(jù)準(zhǔn)確性的影響。如需在路網(wǎng)環(huán)境應(yīng)用，則應(yīng)將相應(yīng)各矢量數(shù)據(jù)替換為專業(yè)工程數(shù)據(jù)。

4.2 設(shè)計(jì)流程

（1）從中國(guó)地震臺(tái)網(wǎng)中心界面抓取地震信息目錄。抓取的時(shí)間間隔可人為設(shè)定，也可人工點(diǎn)擊進(jìn)行抓取，并將結(jié)果單獨(dú)保存至本地。

（2）在抓取到新地震數(shù)據(jù)后進(jìn)行自動(dòng)分析，可對(duì)震中位置進(jìn)行區(qū)域判斷，并根據(jù)地震位置運(yùn)用表1中不同的PGA系數(shù)，計(jì)算結(jié)果為PGA影響范圍的半徑。根據(jù)高速鐵路地震預(yù)警監(jiān)測(cè)系統(tǒng)的規(guī)則習(xí)慣，需計(jì)算當(dāng)PGA=40、80、120 gal或人工輸入值時(shí)對(duì)應(yīng)的范圍半徑，將所有結(jié)果保存至本地。此處預(yù)留衰減關(guān)系替換功能，可根據(jù)運(yùn)用的不同算法替換衰減關(guān)系系數(shù)，以達(dá)到精準(zhǔn)預(yù)測(cè)的目的。

（3）將所得影響分布結(jié)果在QGIS中以新建圖層的方式進(jìn)行可視化展示，震中以“☆”標(biāo)注，大小和顏色隨震級(jí)而定，并配以描述性標(biāo)簽文字；影響范圍是以震中為圓心、步驟（2）中計(jì)算結(jié)果為半徑的圓，在圓外側(cè)的適當(dāng)位置標(biāo)注半徑大小。將以上結(jié)果保存在本地。

（4）人工輸入距離參數(shù)，如100，200 km等。尋找震中到這個(gè)參數(shù)范圍內(nèi)鐵路線最近的位置，并計(jì)算最近的距離；根據(jù)地震衰減公式計(jì)算出鐵路線最近處的地震動(dòng)參數(shù)PGA。同時(shí)，支持判斷該位置處于哪條鐵路線上，將最近位置的經(jīng)緯度坐標(biāo)輸出，并反解析地址。將以上結(jié)果保存在本地。

（5）將步驟（4）的結(jié)果以新建圖層的方式可視化展示，鐵路線最近處默認(rèn)以“×”顯示，與震中連線并標(biāo)注距離。將PGA計(jì)算結(jié)果標(biāo)注在“×”旁邊。

（6）步驟（5）圖層創(chuàng)建后，選取適當(dāng)?shù)慕缑姹壤哌M(jìn)行截圖保存，要求截圖中包含所有計(jì)算結(jié)果生成的圖形。將計(jì)算結(jié)果串連文字輸出。每個(gè)地震輸出一段文字配一個(gè)截圖。將此內(nèi)容保存至本地。

（7）將步驟（6）的文字和截圖鏈接至網(wǎng)站接口，自動(dòng)以網(wǎng)頁(yè)形式進(jìn)行排版編輯，完成后將網(wǎng)址鏈接發(fā)送給短信服務(wù)器。

（8）由短信服務(wù)器根據(jù)信息模板將簡(jiǎn)短文字和網(wǎng)頁(yè)鏈接發(fā)送給訂閱的用戶。

經(jīng)以上步驟，最終實(shí)現(xiàn)地震信息自動(dòng)抓取、分析、自動(dòng)可視化展示、輸出文字，自動(dòng)發(fā)布等一體化全部功能。

4.3 實(shí)現(xiàn)流程

4.3.1 創(chuàng)建基礎(chǔ)化文件

工作計(jì)劃基于QGIS利用Python語(yǔ)言進(jìn)行開(kāi)發(fā)，根據(jù)4.2的設(shè)計(jì)流程，創(chuàng)建插件需要的基礎(chǔ)化文件，包括：

（1）__init__.py，主要描述插件的基本信息、開(kāi)發(fā)者、版本號(hào)、程序主體類、載入方式等。

（2）resources.qrc和resources.py，描述引用的相關(guān)資源文件路徑。

（3）earthquake_analyzer_dialog_base.ui和earth?quake_analyzer_dialog.py，是地震插件界面的ui文件。

（4）earthquake Data.py（地震信息下載和格式化數(shù)據(jù)）、earthquake_analyzer.py（界面交互邏輯）、geo.py（經(jīng)緯度反解析）、mapTool.py（坐標(biāo)轉(zhuǎn)換、測(cè)量、創(chuàng)建圖層等）、worker.py（主要計(jì)算功能）、scheduler.py（間隔定時(shí)自動(dòng)采集）等各類運(yùn)用的算法及函數(shù)。

4.3.2 操作插件

4.3.2.1 插件安裝

在QGIS界面菜單欄中點(diǎn)擊“插件（P）”選項(xiàng)中的“管理并安裝插件...”菜單欄（見(jiàn)圖4）。在插件安裝界面中選擇左側(cè)“從ZIP文件安裝”欄目，選擇安裝包文件后點(diǎn)擊“安裝插件”按鈕（見(jiàn)圖5）。安裝成功后提示“插件已成功安裝”。插件中會(huì)顯示“”圖標(biāo)。

圖4 插件安裝菜單欄

圖5 選擇安裝插件

4.3.2.2 基礎(chǔ)數(shù)據(jù)載入

載入源數(shù)據(jù)圖層，包括鐵路線信息、地震區(qū)劃圖第五代分區(qū)圖和地圖底圖（見(jiàn)圖6），建議數(shù)據(jù)的坐標(biāo)系設(shè)定為WGS 84（EPSG：4326）。

圖6 載入數(shù)據(jù)源圖層

4.3.2.3 插件使用

點(diǎn)擊數(shù)據(jù)采集配置，插件配置參數(shù)見(jiàn)圖7，數(shù)據(jù)采集配置參數(shù)有采集間隔、鐵路線圖層、地震帶圖層、判別距離、考慮深度、存放路徑，短信服務(wù)URL。

圖7 插件配置參數(shù)

（1）輸入?yún)?shù)?！安杉g隔”代表自動(dòng)采集時(shí)從網(wǎng)頁(yè)抓取地震目錄的時(shí)間間隔；“鐵路線圖層”和“地震帶圖層”選擇第4.3.2.2節(jié)中載入的數(shù)據(jù)源；“判別距離”定義為：若震中在此距離范圍內(nèi)存在鐵路線則進(jìn)行分析和計(jì)算，如果此距離范圍內(nèi)沒(méi)有鐵路線，則不執(zhí)行后續(xù)流程；“考慮深度”勾選時(shí)輸入地震目錄中的震源深度參與計(jì)算，此時(shí)式（1）的R由震中距替換為根據(jù)深度幾何關(guān)系計(jì)算得到的震源距；“存放路徑”代表用戶存放截圖和數(shù)據(jù)表格的本地地址；如勾選“自定義半徑”，可根據(jù)用戶設(shè)定的半徑距離，以震中為圓心畫(huà)圓，如不勾選則跳過(guò)此步驟；“短信服務(wù)URL”為短信推送的接口鏈接地址。

（2）采集數(shù)據(jù)。系統(tǒng)軟件采集數(shù)據(jù)支持自動(dòng)采集和手動(dòng)采集2種模式。

手動(dòng)采集數(shù)據(jù)模式是指用戶手動(dòng)點(diǎn)擊系統(tǒng)軟件“手動(dòng)采集”按鈕，會(huì)有窗口彈出，此時(shí)根據(jù)用戶需求輸入數(shù)字n，代表主動(dòng)獲取國(guó)家地震臺(tái)網(wǎng)中國(guó)境內(nèi)地震最新n條目錄，并執(zhí)行后續(xù)步驟。

自動(dòng)采集數(shù)據(jù)模式是指系統(tǒng)軟件根據(jù)第4.3.2.3（1）節(jié)中設(shè)定的采集間隔，自動(dòng)采集獲取國(guó)家地震臺(tái)網(wǎng)中國(guó)境內(nèi)地震最新目錄，并執(zhí)行后續(xù)步驟。長(zhǎng)期監(jiān)測(cè)時(shí)使用自動(dòng)采集模式，實(shí)現(xiàn)鐵路地震影響場(chǎng)的自動(dòng)計(jì)算和發(fā)布。

4.3.2.4 數(shù)據(jù)記錄

執(zhí)行插件后，存放目錄將生成若干個(gè)本地文件夾（見(jiàn)圖8）。文件夾名稱為采集時(shí)間，文件夾中包括當(dāng)時(shí)的地震截圖和數(shù)據(jù)（見(jiàn)圖9），記錄了地震事件基礎(chǔ)信息，包含震級(jí)、發(fā)震時(shí)刻、經(jīng)度、緯度、參考位置、最近的鐵路線名稱、最近鐵路線經(jīng)緯度信息、最近鐵路距離，最近位置預(yù)測(cè)PGA大小及40、80、120 gal半徑和文字描述等數(shù)據(jù)。

圖8 采集數(shù)據(jù)本地文件目錄

圖9 采集數(shù)據(jù)本地文件

以2021年09月21日14時(shí)45分29秒發(fā)生的M3.0級(jí)河南省鄭州市新密市地震為例，輸出文件為地震信息和影響情況，具體為：“據(jù)地震臺(tái)網(wǎng)中心測(cè)定，2021年09月21日14時(shí)45分29秒河南省鄭州市新密市（113.50，34.53）發(fā)生3.0級(jí)地震（疑似塌陷），震源深度0 km。距離附近鐵路線最近位置為新密線河南省鄭州市新密市岳村鎮(zhèn)，最近位置約2.5 km，預(yù)估PGA為54.75 gal。預(yù)估地震40 gal半徑范圍為3.86 km；預(yù)估地震80 gal半徑范圍為1.07 km?！币貓D層生成后地震信息整體可視化展示見(jiàn)圖10。

圖10 地震信息整體可視化展示

4.3.2.5 信息發(fā)布

設(shè)計(jì)短信管理系統(tǒng)，該系統(tǒng)提供短信號(hào)碼訂閱和管理功能，訂閱用戶在地震后可收到鐵路影響信息網(wǎng)頁(yè)鏈接。手機(jī)號(hào)碼狀態(tài)“是”為推送狀態(tài)，手機(jī)號(hào)碼狀態(tài)“否”則為不推送狀態(tài)，短信用戶管理頁(yè)面見(jiàn)圖11。點(diǎn)擊“短信日志信息”可以看到已推送的用戶和手機(jī)短信內(nèi)容。

圖11 短信用戶管理頁(yè)面

信息發(fā)布及推送截圖見(jiàn)圖12，用戶接收的信息見(jiàn)圖12（a），點(diǎn)擊鏈接（圖中已隱去）即可查看最近一次鐵路影響信息（見(jiàn)圖12（b））。

圖12 信息發(fā)布及推送

4.4 測(cè)試結(jié)果

基于QGIS的“鐵路地震影響加速度計(jì)算及發(fā)布插件”開(kāi)發(fā)完成后，開(kāi)展了3個(gè)月的內(nèi)部測(cè)試。因開(kāi)發(fā)所用數(shù)據(jù)和軟件皆為開(kāi)源性質(zhì)，為穩(wěn)定運(yùn)行，環(huán)境部署在互聯(lián)網(wǎng)云服務(wù)器上。

（1）穩(wěn)定性測(cè)試中，插件經(jīng)過(guò)2 000余h測(cè)試無(wú)執(zhí)行錯(cuò)誤、崩潰、內(nèi)存溢出等運(yùn)行故障，能及時(shí)獲取地震信息，并通過(guò)鏈接短信正確發(fā)送文字、圖片等信息至用戶手機(jī)。

（2）準(zhǔn)確性測(cè)試中，插件在周期內(nèi)共計(jì)分析、發(fā)送地震信息1 000余次，其中153次為自動(dòng)獲取的臺(tái)網(wǎng)中心天然地震數(shù)據(jù)，其余為人工添加的歷史地震測(cè)試數(shù)據(jù)，地震加速度影響計(jì)算無(wú)誤。

5 結(jié)論

在詳細(xì)分析鐵路地震需求和合理設(shè)計(jì)實(shí)現(xiàn)流程的基礎(chǔ)上，利用鐵路基礎(chǔ)地理信息數(shù)據(jù)，基于地震動(dòng)衰減關(guān)系開(kāi)發(fā)了一個(gè)開(kāi)源QGIS插件，主要實(shí)現(xiàn)了以下功能：

（1）自動(dòng)抓取國(guó)家地震臺(tái)網(wǎng)中心地震目錄頁(yè)面的地震三要素，并應(yīng)用地震動(dòng)衰減關(guān)系對(duì)地震影響場(chǎng)進(jìn)行分析。

（2）自動(dòng)計(jì)算震源與最近鐵路線的距離，并預(yù)估鐵路線最近位置的地震動(dòng)峰值加速度（PGA）。

（3）對(duì)計(jì)算分析的數(shù)據(jù)和地圖可視化展示，并實(shí)時(shí)發(fā)送給訂閱用戶。

此插件可在地震后及時(shí)通知相關(guān)人員地震影響范圍和最近鐵路的預(yù)估峰值加速度，可為鐵路地震分析、應(yīng)急處置和災(zāi)后救援提供數(shù)據(jù)支撐。