劉 軍 孫甲寧 宋立軍畢雪梅 譚 明 何金剛

1) 中國烏魯木齊830011新疆維吾爾自治區(qū)地震局2) 中國南京210014江蘇省地震局

?

2014年新疆于田MS7.3地震強地面運動模擬及震動圖預測*

劉 軍1),*孫甲寧1)宋立軍1)畢雪梅2)譚 明1)何金剛1)

1) 中國烏魯木齊830011新疆維吾爾自治區(qū)地震局2) 中國南京210014江蘇省地震局

基于地震動的時空衰減規(guī)律和傳播特征， 采用鄰近地震監(jiān)測臺站地震動時程對1 km×1 km尺度的網(wǎng)格點進行近實時插值計算， 同時結合場地效應對震區(qū)地震動參數(shù)進行修正， 并以2014年2月12日新疆于田MS7.3地震為例， 計算震區(qū)格網(wǎng)內(nèi)各點的地震時程， 同時以8 s為時間間隔繪制出地震動峰值等值線圖并將其連續(xù)播放， 得到了于田MS7.3地震峰值地震動(PGV， PGA)的空間分布． 結果表明， 于田縣東部至民豐縣北部地區(qū)受場地條件影響， 震區(qū)震害在軟弱地基土層及淺地下水位等因素作用下對震區(qū)地震動具有明顯的放大效應， 預測的地震動特征與現(xiàn)場宏觀調查結果是一致的． 在當前強震臺網(wǎng)分布不均勻的情形下， 本文方法能較好地描述震區(qū)地震動特征， 較客觀地反映災區(qū)的強地面運動特征．

于田MS7.3地震 地震動 監(jiān)測臺站 震動圖

引言

2014年2月12日新疆于田縣發(fā)生MS7.3地震， 該地震影響范圍涉及和田地區(qū)于田縣、 民豐縣、 策勒縣、 洛浦縣等區(qū)域， 受災戶達到12.4萬戶， 共45萬余人， 直接經(jīng)濟損失為10.8億元， 屬較大破壞性地震， 社會影響較大．

于田MS7.3地震發(fā)生在阿爾金斷裂帶的西段末端， 極震區(qū)位于昆侖山內(nèi)部的高山區(qū)． 阿爾金斷裂帶是青藏高原的北部邊界， 是歐亞大陸內(nèi)部的巨型活動斷裂帶(姚遠等， 2014； 張勇等， 2014)， 為全新世活動斷裂． 李海兵等(2014)的研究結果顯示該斷裂帶以左旋走滑水平運動為主． 利用中國地震烈度表(GB/T 17742—2008)(中華人民共和國國家質量監(jiān)督檢驗檢疫總局， 中國國家標準化管理委員會， 2008)中烈度與強地面運動參數(shù)之間的關系， 在震后黑箱期內(nèi)模擬震區(qū)地震動分布和烈度分布， 可為震后應急救援和災情規(guī)模判定提供參考依據(jù)(陳鯤等， 2010； 孟令媛等， 2014)． 由于新疆“十五”期間布設的強震臺站主要集中在北天山中段和南天山西段， 呈帶狀分布， 臺站分布極不均勻， 孔徑較大(張振斌， 唐麗華， 2010)． 在缺少歷史強震記錄和震后強震臺站觀測數(shù)據(jù)的背景下， 本研究采用數(shù)字地震監(jiān)測臺網(wǎng)觀測數(shù)據(jù)， 利用震區(qū)鄰近臺站數(shù)據(jù)進行插值進而計算格網(wǎng)內(nèi)各點的地震動時程， 便于客觀地描述地震動衰減的特征， 并彌補因震區(qū)強震臺網(wǎng)建設稀疏所造成的信息匱乏． 在此基礎上， 使用鄰近測震臺站記錄計算地震動時程， 以8 s為時間間隔得到各間隔內(nèi)的峰值地震動(PGV， PGA)分布， 進而得到于田MS7.3地震的震動圖． 本研究提供的震動圖不僅應用于新疆地震局于田現(xiàn)場工作隊， 還為地方政府和應急主管部門提供災區(qū)地表震動特征和烈度估計， 為研究地震預期破壞情況及災區(qū)開展應急救援提供了幫助．

1 理論方法

1.1 地震動預測模型

地震波在地球介質中傳播的過程中所產(chǎn)生地震動的強度會隨傳播距離的增加而衰減， 為了客觀地反映地震動傳播的時空衰減特征， 采用地震動峰值的衰減規(guī)律對其進行描述(孟令媛等， 2014； 張振國等， 2014)． 根據(jù)地震發(fā)震構造和加速度記錄資料的分布情況， 本研究選取霍俊榮和胡聿賢(1992)提出的模型， 其形式為

(1)

(2)

式中，Y為地震波加速度、 速度、 位移等地震動峰值，M為震級，R為震中距，R0為與震級相關的近場距離飽和因子，C1，C2，C3，C4，C5為回歸系數(shù)．

基于上述衰減關系模型， 利用新疆天山地區(qū)1985—2014年MS≥4.0歷史地震動數(shù)據(jù)記錄， 運用兩步回歸法進行回歸， 得到新疆天山地區(qū)的地震動衰減關系的系數(shù) (張振斌， 唐麗華， 2010) ．

在同一地震作用下， 震中距分別為R和RA的兩個臺站的地震動峰值Y(R)與Y(RA)之間的關系為(金星等， 2010)

(3)

Y(R)用Y(RA)表示為

(4)

由上式可以看出， 震中距為R的地震動峰值可采用RA處的地震動峰值Y(RA)與補償函數(shù)f(R，RA)相乘進行計算．

根據(jù)張紅才(2008)提出的方法， 任一點D的地震動時程可以采用鄰近的A，B，C等3個臺站的地震動時程進行插值計算得到， 可表示為

(5)

式中，ωA，ωB，ωC為插值計算的加權系數(shù)， 可通過震中距倒數(shù)平方加權的方法得到較好的展現(xiàn)， 且具有較高的計算效率， 表達式如下:

(6)

A，B，C這3個臺站相對于插值點D的時間延遲ΔtAD， ΔtBD， ΔtCD可表示為

(7)

式中：ΔAD，ΔBD，ΔCD分別為臺站A，B，C與插值點D之間的距離；RA，RB，RC，RD分別為A，B，C這3個臺站和插值點D的震中距． 為客觀反映地震動的衰減特性， 將式(5)中震中距換為震源距， 式(5)可改寫為

(8)

如圖1b(張紅才， 2008)所示， 若臺站B所在位置為需要插值計算的某處， 其P波到時、 最大值到時分別為tP和tmax． 計算該處地震動時程時，P波到達至最大值到達之間的時程采用上述方法進行計算； 而對于P波未到之前的地震動時程則為周圍臺站脈動記錄的加權平均； 最大值到達之后時刻的時程則認為地震動以S波速度vS進行傳播， 如式(9)所表示．

圖1 插值計算方法(a)及插值原理(b)示意圖(引自張紅才， 2008)

(9)

1.2 不同場地條件下地震動參數(shù)校正

為了客觀地反映震區(qū)地震動特征分布， 獲取震區(qū)土層的放大系數(shù)， 本文采用第五代地震區(qū)劃圖(簡稱為五代圖)(GB/T 18306—2015)(中華人民共和國國家質量監(jiān)督檢驗檢疫總局， 中國國家標準化管理委員會， 2015)中不同場地條件下的放大系數(shù)． 五代圖計算出了不同場地條件下短周期分檔加速度的地震動參數(shù)的放大系數(shù)， 如表1所示． 本文計算中， 首先根據(jù)1.1節(jié)中所述方法進行初步插值計算得到PGA， 然后采用五代圖中的場地放大系數(shù)對地震動參數(shù)進行校正．

表1 場地放大系數(shù)

2 于田MS7.3地震震動圖

2.1 計算區(qū)域

圖2 本文所用監(jiān)測臺站分布圖Fig.2 Distribution of monitoring stations used in this paper

本文計算區(qū)域如圖2所示， 主要涉及新疆和田地區(qū)及周邊區(qū)域． 根據(jù)前期準備的千米網(wǎng)格(1 km×1 km)數(shù)據(jù)判斷， 于田MS7.3地震共有18萬6067個千米網(wǎng)格地震動數(shù)據(jù)待計算． 由圖2可以看到該地震周邊監(jiān)測臺站分布較為均勻．

本文選取14個臺站進行插值計算， 其相關參數(shù)列于表2．

為了彌補臺站稀疏地區(qū)測震數(shù)據(jù)質量較低的缺陷， 首先根據(jù)新疆地區(qū)地震動衰減關系估計方法， 計算了震區(qū)基巖面的PGA和PGV地震動參數(shù); 然后利用五代圖在不同場地條件下的放大系數(shù)將地震動參數(shù)校正到相應土層地表上， 計算得到震區(qū)周邊地區(qū)地表的PGA和PGV； 最后根據(jù)測震臺站對于本次地震的觀測記錄， 采用距離加權插值法對以震區(qū)臺站為中心、 30 km×30 km范圍內(nèi)的千米格網(wǎng)點的PGA和PGV值進行重新計算， 得到網(wǎng)格節(jié)點的地震動參數(shù)并繪制出地震烈度圖(Borcherdt， 1994)． 在計算地震動參數(shù)過程中， 加大了實際觀測數(shù)據(jù)參數(shù)的權重， 便于發(fā)揮震區(qū)強震臺站的作用．

表2 本文所用臺站及其主要參數(shù)

2.2 震動圖計算

本文以1 km×1 km為尺度對震區(qū)周邊30 km×30 km區(qū)域進行網(wǎng)格化離散后， 結合新疆及青海周邊地區(qū)數(shù)字地震臺網(wǎng)記錄到的于田地震資料， 采用上文所述方法計算區(qū)域內(nèi)所有網(wǎng)格點的地震動時程， 并結合震區(qū)場地條件對地震動參數(shù)進行校正， 最終提取 以8 s 為時間間隔所計算的速度、 加速度產(chǎn)生的地震動參數(shù)快照． 臺站稀疏地區(qū)則結合地震動衰減關系進行補充， 發(fā)震斷層主要用于判別長軸的方向，MS≥7.0地震采用線源衰減關系模型(李鉑等， 2012)． 于田地震加速度和速度震動圖結果分別如圖3和圖4所示． 每幅子圖為快照提取圖， 時間間隔為8 s， 從圖中可以看出于田地震的加速度和速度衰減特性.

圖3 2014年于田MS7.3地震加速度震動圖

2.3 峰值地震動分布

通過上述方法插值計算得到于田MS7.3地震的PGA和PGV地震動分布， 如圖5所示． 震動圖描繪了該地震引起的峰值地面運動分布， 震后初期可以根據(jù)地震動預測的地震動強度判定地震災害的影響和規(guī)模． 利用中國地震烈度表(GB/T 17742—2008)(中華人民共和國國家質量監(jiān)督檢驗檢疫總局， 中國國家標準化管理委員會， 2008)中地震動參數(shù)與烈度之間的關系將地震動值轉換為烈度值， Ⅶ度及以上區(qū)域采用PGV進行計算， Ⅶ度以下區(qū)域則采用PGA進行衡量．

根據(jù)和田地區(qū)水文資料顯示， 于田縣英巴格鄉(xiāng)以東至奧依牙依拉克鄉(xiāng)屬于塔里木盆地， 地下水埋深普遍為1.8—2.5 m， 其中尼雅鄉(xiāng)地下水埋深平均約1.3 m， 震區(qū)位于昆侖山北麓平原區(qū)， 總體地勢南高北低， 沉積物顆粒由南向北逐漸變細， 靠近山前屬于礫質平原區(qū)(劉軍等， 2014)． 利用場地覆蓋層厚度與土層的平均等效剪切波速的雙參數(shù)法， 確定于田縣英巴格鄉(xiāng)至奧依牙依拉克鄉(xiāng)以東區(qū)域為Ⅲ類場地， 如圖6所示． 該區(qū)域理論PGA在65—86 cm/s2之間， 具有優(yōu)勢分布． 由于場地條件影響， Ⅲ類場地對該地區(qū)PGA具有明顯放大效應. 從圖7中可以看出， 本次計算所得到的地震烈度分布在于田縣、 民豐縣以北大部分地區(qū)， 這些地區(qū)為軟弱地基土層， 地下水位淺， 在地震動的作用下， 震害顯著加?。?Ⅶ度區(qū)努爾鄉(xiāng)、 博斯坦鄉(xiāng)及葉亦克鄉(xiāng)的房屋破壞較為嚴重， 于田縣城、 民豐縣城以北的Ⅵ度區(qū)影響范圍較大． 圖7給出了本文所計算的烈度分布與現(xiàn)場調查烈度的對比， 可以看出， 利用震動圖預測的地震動特征與現(xiàn)場宏觀調查結果是一致的．

圖4 2014年于田MS7.3地震速度震動圖

圖5 插值網(wǎng)格點峰值加速度(a)和峰值速度(b)等值線圖

圖6 研究區(qū)場地條件分布圖

圖7 預測與現(xiàn)場調查烈度對比圖

3 討論與結論

本文依據(jù)構造地震動時程的插值方法， 應用鄰近測震記錄臺站地震動時程插值地震動時程的模型， 并采用千米格網(wǎng)離散化進行插值計算， 得到2014年2月12日于田MS7.3地震影響區(qū)內(nèi)格網(wǎng)點的地震動時程， 繪制了該地震的震動圖， 并對震動圖計算結果予以分析． 本文設計的震動圖計算方法具有以下特征：

1) 嘗試利用新疆及鄰近省份的測震臺站觀測數(shù)據(jù)對缺少臺站的阿爾金地震帶的地震動時程進行插值計算， 不僅較好地發(fā)揮了現(xiàn)有測震臺網(wǎng)資源的作用， 同時又較客觀地模擬了震區(qū)強地面運動情況．

2) 通過計算得到新疆于田MS7.3地震的震動圖， 該圖可以直觀地描述震區(qū)地面運動的整個過程， 而其地震動(PGA， PGV)等值線圖能較全面地反映該地震對地表的影響和地震動時空衰減特征．

3) 利用新疆地區(qū)不同場地條件下的地震動放大效應對臺站地震動時程插值進行校正， 考慮了發(fā)震斷層、 局部場地地形條件以及震源復雜性影響的震動圖更能反映出災區(qū)的地震動空間分布特征． 由于當前新疆地區(qū)強震觀測臺站主要分布在天山地區(qū)， 局部地區(qū)的結果因當前臺站分布不均等因素， 其可靠性在一定程度上受到影響， 但其影響必將隨著預警工程臺站的建設得到逐漸改善. 此外， 地震動衰減關系的引入在一定程度上能夠彌補局部地區(qū)臺站稀疏等對計算結果的影響．

另一方面， 本研究仍存在一些問題亟待解決， 也是進一步深入研究的方向：

1) 采用最近距離原則選用臺站， 即單純選用與待插值點距離最近的觀測臺站， 在插值過程中并未考慮臺站與插值點之間的空間分布和相對位置， 所以其結果的準確性有待于進一步研究與驗證， 可嘗試引入臺站相對位置加權項等方法來減弱由此造成的影響.

2) 在新疆北天山中段和南天山西段強震臺站分布較為密集的情況下， 可以嘗試采用實時觀測強震數(shù)據(jù)與測震數(shù)據(jù)相結合的方式進行插值計算， 并可利用地震動衰減關系和地震臺站記錄的數(shù)據(jù)研究地震動強度衰減中地震事件內(nèi)部的不確定性． 通過臺站實際觀測值與理論計算值之間的殘差分析， 應用測震臺站觀測數(shù)據(jù)和偏差校正方法， 對本文設計的插值方法進行校正， 以提高觀測臺站空缺地區(qū)的震動圖精度． 同時， 在臺站稀疏的地區(qū)也可以結合區(qū)域地震動衰減關系、 經(jīng)驗模型和斷層方向等因素進行計算和補充．

陳鯤， 俞言祥， 高孟潭. 2010. 考慮場地效應的ShakeMap系統(tǒng)研究[J]. 中國地震， 26(1): 92--102.

Chen K， Yu Y X， Gao M T. 2010. Research on ShakeMap system in terms of the site effect[J].EarthquakeResearchinChina， 26(1): 92--102 (in Chinese).

霍俊榮， 胡聿賢. 1992. 地震動峰值參數(shù)衰減規(guī)律的研究[J]. 地震工程與工程振動， 12(2)： 1--11.

Huo J R， Hu Y X. 1992. Study on attenuation laws of ground motion parameters[J].EarthquakeEngineeringandEngineeringVibration， 12(2): 1--11 (in Chinese).

李鉑， 董翔， 陳亞紅， 何琳， 于澄. 2012. 線源地震影響場計算模型[J]. 地震地磁觀測與研究， 33(3): 14--19.

Li B， Dong X， Chen Y H， He L， Yu C. 2012. Realization of calculation model of line-source effect fields[J].SeismicGeomagneticObservationandResearch, 33(3): 14--19 (in Chinese).

李海兵， 孫知明， 潘家偉， 劉棟梁， 張佳佳， 李成龍， 劉亢， 云錕， 龔正. 2014. 2014年于田MS7.3地震野外調查： 特殊的地表破裂帶[J]. 地球學報， 35(3): 391--394.

Li H B， Sun Z M， Pan J W， Liu D L， Zhang J J， Li C L， Liu K， Yun K， Gong Z. 2014. Field study of the 12 February 2014 YutianMS7.3 earthquake: A special surface rupture zone[J].ActaGeoscienticaSinica， 35(3): 391--394 (in Chinese).

劉軍， 宋立軍， 胡偉華， 李帥. 2014． 新疆地區(qū)場地條件分區(qū)及其矢量化[J]． 內(nèi)陸地震， 28(1)： 70--74．

Liu J， Song L J， Hu W H， Li S. 2014． Zoning and vectorize of site condition in Xinjiang[J].InlandEarthquake， 28(1)： 70--74 (in Chinese).

金星， 張紅才， 韋永祥， 李軍. 2010. 基于地震監(jiān)測臺網(wǎng)資料近實時插值計算震動圖的初步研究[J]. 防災減災學報， 26(1): 1--11.

Jin X， Zhang H C， Wei Y X， Li J. 2010. Preliminary study of near real-time interpolation calculation of ShakeMap based on seismic monitoring network’s information[J].JournalofDisasterPreventionandReduction， 26(1): 1--11 (in Chinese).

孟令媛， 周龍泉， 劉杰. 2014. 2014年新疆于田MS7.3地震近斷層強地面運動模擬及烈度分布估計[J]. 地震學報， 36(3)： 362--371.

Meng L Y， Zhou L Q， Liu J. 2014. Estimation of near-fault strong ground motion and intensity distribution of the 2014 Yutian， XinjiangMS7.3 earthquake[J].ActaSeismologicaSinica, 36(3): 362--371 (in Chinese).

姚遠， 伊力亞爾·阿不力孜， 吳傳勇， 吳國棟， 陳建波. 2014. 2014年2月12日新疆于田MS7.3地震發(fā)震構造初步研究[J]. 內(nèi)陸地震， 28(2)： 121--126．

Yao Y， Yiliyar·Abulizi， Wu C Y, Wu G D, Chen J B. 2014. Preliminary study on seismogenic structure of YutianMS7.3 earthquake on Feb.12， 2014[J].InlandEarthquake, 28(2): 121--126 (in Chinese).

張紅才. 2008. 基于地震監(jiān)測臺網(wǎng)資料的震動圖及震動烈度研究[D]. 哈爾濱： 中國地震局工程力學研究所： 17--23.

Zhang H C. 2008.ResearchonShakeMapsandShakingIntensityBasedonRecordsofMonitoringNetwork[D]. Harbin: Institute of Engineering Mechanics, China Earthquake Administration: 17--23 (in Chinese).

張勇， 許力生， 陳運泰， 汪榮江. 2014. 2014年2月12日于田MW6.9地震破裂過程初步反演: 兼論震源機制對地震破裂過程反演的影響[J]. 地震學報， 36(2): 159--164.

Zhang Y， Xu L S， Chen Y T， Wang R J. 2014. Fast inversion for the rupture process of the 12 February 2014 YutianMW6.9 earthquake: Discussion on the impacts of focal mechanisms on rupture process inversions[J].ActaSeismologicaSinica， 36(2): 159--164 (in Chinese).

張振斌， 唐麗華. 2010. 新疆南天山地區(qū)土層場地地震動峰值加速度衰減關系[J]. 中國地震， 26(4): 401--406.

Zhang Z B， Tang L H. 2010. The attenuation relationship of peak ground motion acceleration on soil site at the south Tianshan region in Xinjiang[J].EarthquakeResearchinChina， 26(4): 401--406 (in Chinese).

張振國， 張偉， 孫耀充, 朱耿尚， 文健， 陳曉非. 2014. 2014年2月12日新疆于田地震強地面運動初步模擬及烈度預測[J]. 地球物理學報， 57(2): 685--689．

Zhang Z G， Zhang W， Sun Y C， Zhu G S, Wen J, Chen X F. 2014. Preliminary simulation pf strong ground motion for Yutian， Xinjiang earthquake of 12 February 2014， and hazard implication[J].ChineseJournalofGeophysics， 57(2): 685--689 (in Chinese).

中華人民共和國國家質量監(jiān)督檢驗檢疫總局， 中國國家標準化管理委員會. 2008. 中國地震烈度表(GB/T17742—2008)[S]. 北京: 中國標準出版社： 5--7.

General Administration of Quality Supervision， Inspection and Quarantine of the People’s Republic of China， Standardization Administration of the People’s Republic of China. 2008.TheChineseSeismicIntensityScale(GB/T17742--2008)[S]. Beijing: Standards Press of China: 5--7.

中華人民共和國國家質量監(jiān)督檢驗檢疫總局， 中國國家標準化管理委員會. 2015. 中國地震動參數(shù)區(qū)劃圖(GB/T18306—2015)[S]. 北京： 中國標準出版社: 6--10.

General Administration of Quality Supervision， Inspection and Quarantine of the People’s Republic of China， Standardization Administration of the People’s Republic of China. 2015.ChinaSeismicParameterZoningMap(GB/T18306--2015)[S]. Beijing: Standards Press of China: 6--10.

Borcherdt R D. 1994. Estimates of site-dependent response spectra for design (methodology and justification)[J].EarthquakeSpectra， 10(4): 617--653．

Research on strong ground motion simulation and ShakeMap of 2014 YutianMS7.3 earthquake

Liu Jun1),*Sun Jianing1)Song Lijun1)Bi Xuemei2)Tan Ming1)He Jingang1)

1)EarthquakeAdministrationofXinjiangUygurAutonomousRegion， ürümqi830011，China2)EarthquakeAdministrationofJiangsuProvince,Nanjing210014,China

Based on attenuation law and space-time propagation characteristics of seismic ground motion, we use recorded time histories from adjacent monitoring stations to make near real-time interpolation calculation in a 1 km×1 km grid, and correct ground motion parameters according to the site effects. Taking the 12th February, 2014 YutianMS7.3 earthquake as an example, we calculated seismic time history in each grid near the epicenter, and then plot peak ground motion contour maps with 8 s intervals continuously. Finally we showed spatial distribution of peak ground motions (PGV and PGA) of this earthquake. The results show that, influenced by site conditions as well as soft ground soil and shallow groundwater level, ground motions in the region from eastern Yutian County to northern Minfeng County were amplified obviously, which were well consistent with field survey results. On the conditions of uneven distribution of strong ground motion observation stations, the method can be used to describe the characteristics of earthquake ground motions, and reflect the characteristics of strong ground motion in disaster areas objectively.

YutianMS7.3 earthquake; seismic ground motion; monitoring station; ShakeMap

中國地震局 “三結合”項目(163101)與地震科技星火攻關項目(XH15051)聯(lián)合資助．

2015-09-28收到初稿， 2016-07-22決定采用修改稿.

10.11939/jass.2016.06.011

P315.9

A

劉軍, 孫甲寧, 宋立軍, 畢雪梅， 譚明, 何金剛. 2016. 2014年新疆于田MS7.3地震強地面運動模擬及震動圖預測. 地震學報, 38(6): 914--923. doi:10.11939/jass.2016.06.011.

Liu J, Sun J N, Song L J, Bi X M, Tan M, He J G. 2016. Research on strong ground motion simulation and ShakeMap of 2014 YutianMS7.3 earthquake.ActaSeismologicaSinica, 38(6): 914--923. doi:10.11939/jass.2016.06.011.

*通訊作者 e-mail: liujun_eq@sina.com