韓曉明 劉 芳 張 帆 胡 博 陳立峰

（中國呼和浩特010010內(nèi)蒙古自治區(qū)地震局）

據(jù)中國地震臺網(wǎng)測定，北京時間2015年4月15日15時39分，在內(nèi)蒙古自治區(qū)阿拉善左旗（39.8°N，106.3°E）發(fā)生MS5.8地震，震源深度為10km，宏觀震中位于巴音木仁蘇木.由于該地震發(fā)生在烏蘭布和沙漠，人口居住密度小，未造成顯著的地表破壞和人員傷亡；其極震區(qū)烈度達Ⅶ度，沿NE向展布，長軸為17km，短軸為9km，面積達126km2；Ⅵ度區(qū)長軸為90km，短軸為62km，面積為4 904km2①內(nèi)蒙古地震局現(xiàn)場工作隊.2015.內(nèi)蒙古自治區(qū)阿拉善左旗5.8地震現(xiàn)場科學考察報告..

阿拉善左旗MS5.8地震發(fā)生在具有右旋剪切破裂特征的巴彥浩特盆地內(nèi)部（鄧起東，尤惠川，1985），屬華北地區(qū)中強地震第四活動期（1815年至今）最新活躍幕發(fā)生的一次中強地震.其余震序列不太發(fā)育，截至2015年5月4日0時，共記錄到余震158次（含單臺記錄），其中ML0.0—0.9余震27次，ML1.0—1.9余震106次，ML2.0—2.9余震20次，ML3.0—3.9余震3次，ML4.0—4.9余震2次，最大余震為ML4.5，為主余型地震.這些余震集中發(fā)生在4月15—17日，共105次，占該序列地震總數(shù)的67%.

本文采用CAP（cut and paste）方法（Zhu，Helmberger，1996）測定了阿拉善左旗MS5.8主震的最佳雙力偶解.結果顯示，該地震為走滑型，最佳矩心深度為13km，如圖1所示.表1給出了部分研究機構所獲得的震源機制解，可以看出，各機構獲得的震源機制類型和震源深度均存在一定差異，本文結果與USGS的結果比較相近.

圖1 采用CAP方法反演的2015年阿拉善左旗MS5.8地震震源機制解及相關結果（a）理論（紅色）與觀測地震圖（黑色）的波形擬合，波形左側為臺站代碼，其上下數(shù)字分別表示震中距（單位：km）和方位角，波形下方第一行數(shù)字為理論地震圖相對觀測地震圖的時移（單位：s），第二行為二者的相關系數(shù)；（b）震源深度為13km時的最佳雙力偶解（下半球投影）；（c）震源機制解誤差隨震源深度的變化Fig.1 Focal mechanism solution of the 2015Alxa Zuoqi MS5.8earthquake by CAP method and related results（a）Waveform fitting of theoretical waveforms（red lines）with observational ones（black lines）.Station codes are given on the left of waveforms with epicentral distances（in unit of km）at the top and azimuths at the bottom，and the numbers of the first line below the waveforms are the time shift（in unit of s），the numbers of the second line are their correlative coefficient.（b）The best double couple solution（lower hemisphere projection，focal depth is 13km）.（c）The misfit error distribution of focal mechanism solution as a function of focal depth

表1 2015年阿拉善左旗MS5.8地震的震源機制解參數(shù)Table 1 Focal mechanism solution parameters of the 2015Alxa Zuoqi MS5.8earthquake

采用自助線性方法（Michael，1991）反演MS5.8地震震中區(qū)（圖2中虛線框）的應力場.結果表明：震中區(qū)應力場的最大主壓應力方位角為52°，傾角為20°，作用方向為NE；其張應力方位角為308°，傾角為33°，作用方向為NW；應力傾角近水平方向，這也是研究區(qū)應力場具有走滑特征的一個表現(xiàn)（圖2）.運用振幅比方法（梁尚鴻等，1984）求解2001年以來鄂爾多斯西北緣226次ML≥2.5地震（不含MS5.8地震序列）的震源機制解，反演計算了鄂爾多斯西北緣的小震滑動角與斷層滑動角的偏差分布（Michael，1991），其與震源機制一致性參數(shù)具有類似含義（陳颙，1978），如圖2所示.可以看出，滑動角偏差值在（105°E—108°E，39.5°N—40.0°N）范圍內(nèi)的磴口—本井斷裂中段存在低值，MS5.8地震則發(fā)生在低值區(qū)西緣.

為了分析余震的空間展布特征，本文基于阿拉善地塊的速度結構模型（楊明芝等，2007）（表2），對阿拉善左旗地震運用共軛梯度算法進行 HypoDD雙差定位（Waldhauser，Ellsworth，2000；楊智嫻等，2003）.經(jīng)篩選，得到能被6個以上臺站記錄到的ML≥1.0地震共78次，重新定位保留64次地震（圖3a）.

圖2 2015年阿拉善左旗MS5.8地震的滑動角偏差分布和震中區(qū)（虛線框）的應力場反演結果Fig.2 Distribution of slip angle deviation of the 2015Alxa Zuoqi MS5.8earthquake and stress field inversion results of its epicentral area delineated by the dashed rectangle

圖3 a給出了阿拉善左旗地震的重定位結果.可以看出：該地震位于（39.78°N（±0.7km），106.34°E（±0.7km）），震源深度為（13±1）km，余震優(yōu)勢分布方向為ENE和NNW，分別延伸約16和18km.由圖3b可以看出，阿拉善左旗地震的震源深度由定位前的3—38km縮小至重定位后的3—19km.

以正態(tài)分布曲線進行擬合考察定位誤差，結果如圖4所示.可以看出：E-W 向、N-S向和 U-D向的誤差分別為0.5—1.6，0.5—1.8和0.5—1.9，均數(shù)分別為0.7，0.8和1.0（圖4a）；E-W 向和 N-S向的定位誤差在0—1.2km區(qū)間的累積分布概率分別為0.97和0.86，U-D向定位誤差在0—1.5km區(qū)間的累積分布概率為0.90（圖4b）；走時殘差的正態(tài)分布均數(shù)為0.17，分布范圍為0.04—0.28s（圖4a），在0—0.25s區(qū)間的累積分布概率為0.93（圖4b）.

表2 阿拉善地塊地殼速度結構（引自楊明芝等，2007）Table 2 Crustal velocity structure of Alxa block（after Yang et al，2007）

圖3 2015年阿拉善左旗MS5.8地震序列重新定位結果（a）及重新定位前（上）、后（下）的震源深度分布（b）Fig.3 Relocation of the 2015Alxa Zuoqi MS5.8earthquake sequence（a）and distribution of focal depth before（upper）and after（lower）relocation（b）

圖4 2015年阿拉善左旗MS5.8地震序列重新定位結果的誤差統(tǒng)計（a）E-W向、N-S向、U-D向誤差和走時殘差的正態(tài)分布統(tǒng)計結果；（b）對應的累積概率密度分布Fig.4 Statistic on error of relocation for the 2015Alxa Zuoqi MS5.8earthquake sequence（a）Normal distribution statistics of relocation error in E-W，N-S，U-D and travel time residual；（b）Their cumulative probability density distribution

阿拉善左旗MS5.8地震發(fā)生在鄂爾多斯西北緣的巴彥浩特盆地，整體處于右旋剪切拉張作用的應力場環(huán)境中.2001年以來，震中區(qū)的應力張量作用比較均勻，且與應力場背景特征協(xié)調，應力場的繼承性作用特征使得NW向的拉張應力占主導作用，其一致性的應力加載作用有利于MS5.8地震發(fā)生右旋剪切破裂，應力場反演結果較為明確地支持N-S向節(jié)面為主破裂面，并與1976年巴音木仁MS6.2地震的發(fā)震機制相似，國家地震局《鄂爾多斯周緣斷裂系》課題組（1988）對巴音木仁MS6.2地震破裂面的主要判據(jù)為，N-S向的破裂面與賀蘭山構造帶走向接近，且震中區(qū)南側E-W走向的紅果子溝長城曾被NE35°展布的斷層錯移產(chǎn)生右旋平移，故N-S向的破裂面較符合該區(qū)構造運動實況.

地震重新定位使用了中國地震局地球物理研究所實時地震學與地震監(jiān)測預報技術研究室提供的“中國地震科學探測臺陣數(shù)據(jù)”，文中圖件由GMT和Zmap軟件包繪制完成，審稿專家提出了科學合理的修改建議，作者在此一并表示感謝.

陳颙.1978.用震源機制一致性作為描述地震活動性的新參數(shù)［J］.地球物理學報，21（2）：146-159.

Chen Y.1978.Consistency of focal mechanism as a new parameter in describing seismic activity［J］.Acta Geophysica Sinica，21（2）：146-159（in Chinese）.

鄧起東，尤惠川.1985.鄂爾多斯周緣斷陷盆地帶的構造活動特征及其形成機制［G］∥現(xiàn)代地殼運動研究（l）.北京：地震出版社：58-78.

Deng Q D，You H C.1985.Mechanism and structure activity of Ordos basin［G］∥Study on the Modern Crustal Movement（1）.Beijing：Seismological Press：58-78（in Chinese）.

國家地震局《鄂爾多斯周緣活動斷裂系》課題組.1988.鄂爾多斯周緣活動斷裂系［M］.北京：地震出版社：39-76.

Research Group on Active Fault System Around Ordos Massif，State Seismological Bureau.1988.Active Fault System Around Ordos Massif［M］.Beijing：Seismological Press：39-76（in Chinese）.

粱尚鴻，李幼銘，束沛鎰，朱碚定.1984.利用區(qū)域地震臺網(wǎng)P、S波振幅比資料測定小震震源參數(shù)［J］.地球物理學報，27（3）：249-257.

Liang S H，Li Y M，Shu P Y，Zhu B D.1984.On the determination of source parameters of small earthquakes by using amplitude ratios of P，S waveform from regional network observation［J］.Acta Geophysica Sinica，27（3）：249-257（in Chinese）.

楊明芝，馬禾青，廖玉華.2007.寧夏地震活動與研究［M］.北京：地震出版社：11-46.

Yang M Z，Ma H Q，Liao Y H.2007.NingxiaEarthquakeActivityandResearch［M］.Beijing：Seismological Press：11-46（in Chinese）.

楊智嫻，陳運泰，鄭月軍，于湘?zhèn)?2003.雙差地震定位法在我國中西部地區(qū)地震精確定位中的應用［J］.中國科學：D輯，33（S1）：129-134.

Yang Z X，Chen Y T，Zheng Y J，Yu X W.2003.Accurate relocation of earthquake in central-western China using the double difference earthquake location algorithm［J］.Science in China：Series D，46（S2）：181-188.

中國地震局地球物理研究所.2015.2015年4月15日內(nèi)蒙古阿拉善左旗5.8級地震［EB／OL］.［2015-04-20］.http：∥www.cea-igp.ac.cn／tpxw／272036.shtml.

Institute of Geophysics，China Earthquake Administration.2015.The MS5.8Alxa，Inner Mongolia，earthquake on April 15，2015［EB／OL］.［2015-04-20］.http：∥www.cea-igp.ac.cn／tpxw／272036.shtml（in Chinese）.

GCMT.2015.M5.4：36km WNW of Wuhai，China［EB／OL］.［2015-04-20］.http：∥www.globalcmt.org／cgi-bin／globalcmt-cgi-bin／CMT4／form？itype＝y(tǒng)md＆yr＝2015＆mo＝4＆day＝14＆otype＝y(tǒng)md＆oyr＝2015＆omo＝4＆oday＝15＆jyr＝1976＆jday＝1＆ojyr＝1976＆ojday＝1＆nday＝1＆lmw＝5＆umw＝10＆lms＝0＆ums＝10＆lmb＝0＆umb＝10＆llat＝35＆ulat＝41＆llon＝105＆ulon＝108＆lhd＝0＆uhd＝1000＆lts＝－9999＆uts＝9999＆lpe1＝0＆upe1＝90＆lpe2＝0＆upe2＝90＆list＝0.

GFZ.2015.Western Nei Mongol，China［EB／OL］.［2015-04-20］.http：∥geofon.gfz-potsdam.de／eqinfo／event.php？id＝gfz2015hipm.

Michael A J.1991.Spatial variations in stress within the 1987Whittier Narrows，California，aftershock sequence：New techniques and results［J］.J Geophys Res，96（B4）：6303-6319.

USGS.2015.M5.4：36km WNW of Wuhai，China［EB／OL］.［2015-04-20］.http：∥earthquake.usgs.gov／earthquakes／eventpage／us200026mr＃scientific＿tensor.

Waldhauser F，Ellsworth W L.2000.A double difference earthquake location algorithm：Method and application to the northern Hayward fault，California［J］.Bull Seismol Soc Am，90（6）：1353-1368.

Zhu L P，Helmberger D V.1996.Advancement in source estimation techniques using broadband regional seismograms［J］.Bull Seismol Soc Am，86（5）：1634-1641.