韓曉明劉芳張帆魏建民徐巖侯迪李娟梁沙沙
1 內(nèi)蒙古自治區(qū)地震局,呼和浩特市哲里木路80號,010010
?
河套地震帶中小地震的震源機制測定及應力場特征
韓曉明1劉芳1張帆1魏建民1徐巖1侯迪1李娟1梁沙沙1
1內(nèi)蒙古自治區(qū)地震局,呼和浩特市哲里木路80號,010010
摘要:運用振幅比和CAP方法測定224次地震的震源機制解,進行歸一頻度統(tǒng)計和震源機制類型時空分析,以考察河套地震帶的應力場在模式轉(zhuǎn)換和應力作用方式兩方面的變化形態(tài)。河套地震帶的壓應力方位在北東向存在優(yōu)勢分布,張應力方位在北西向存在優(yōu)勢分布,滑動角主要分布在近水平方向,應力場整體表現(xiàn)出繼承性運動特征,表明水平剪切拉張運動總體控制了帶內(nèi)斷裂的作用模式而多產(chǎn)出走滑型地震,并貫穿整個研究時段,屬穩(wěn)態(tài)應力作用模式。在太平洋板塊向西俯沖和青藏塊體向北運移的聯(lián)合作用下,河套地震帶兩側(cè)的陰山隆起和鄂爾多斯塊體作相對旋性扭動,最終驅(qū)使帶內(nèi)眾多北東-近東西向展布的斷裂產(chǎn)生水平旋性剪切,構(gòu)成相對穩(wěn)定的應力場演化模式。
關鍵詞:河套地震帶;震源機制解;振幅比;CAP;應力場特征
河套地震帶屬于中國東部華北至東北應力區(qū)的次級應力區(qū)域,最大主壓應力方向為北東(臨河地區(qū))-北東東(呼和浩特至包頭地區(qū)),構(gòu)造應力張量結(jié)構(gòu)以走滑型為主(圖1)。本文基于河套地震帶的構(gòu)造應力場背景特征,根據(jù)2001年以來發(fā)生在河套地震帶的中小地震事件,運用振幅比和CAP方法測定震源機制,通過歸一頻度統(tǒng)計和斷層滑動類型的時空分析,考察河套地震帶的應力場在模式轉(zhuǎn)換和應力作用方式兩方面的變化形態(tài),并結(jié)合外圍動力學環(huán)境和區(qū)域斷裂構(gòu)造的運動性質(zhì)闡釋其變化原因。
1研究區(qū)地質(zhì)構(gòu)造背景及地震活動
河套地震帶位于陰山隆起與鄂爾多斯臺地之間,是鄂爾多斯周緣斷陷帶中規(guī)模最大、構(gòu)造活動和垂直差異運動十分顯著的地區(qū)。第四紀以來,陰山隆起與鄂爾多斯塊體之間作相對左旋的剪切拉張運動[1],這種運動特征通過1970年以來河套地震帶發(fā)生的4次6級以上地震的震源機制解能得到體現(xiàn),震源機制解所表達的張應力方位為NW-NNW,壓應力方位為NE-NEE。其中發(fā)生在臨河盆地周緣的兩次地震的張應力為NW,壓應力為NE,發(fā)生在呼包盆地周緣的兩次地震的張應力為NNW,壓應力為NEE,表明河套地震帶次級構(gòu)造單元的應力場方向特征并非完全一致[2]。
相對于歷史地震的高強度,研究區(qū)近代以來的中強地震表現(xiàn)出高頻度的特點[3-4]。1920年以來,華北地區(qū)的中強地震出現(xiàn) “填空性”活動特征,而位于鄂爾多斯塊體北緣的河套地震帶也在長期平靜后發(fā)生了5次6級以上地震。其中,1976~1979年連續(xù)發(fā)生3次6級以上地震,形成小規(guī)模叢集活躍,使其成為鄂爾多斯塊體周緣中強地震最為活躍的區(qū)域(圖1)。由此可見,河套地震帶的斷裂運動模式和應力場積累規(guī)模具備產(chǎn)出強烈地震的條件。
2震源機制求解
選擇振幅比[5-6]和CAP(cut and paste)[7]方法進行震源機制測定,這兩種方法均可削弱臺站布局、介質(zhì)非均勻性和速度結(jié)構(gòu)模型的影響而得到相對精確的震源機制解。ML≥4.0地震采用基于波形擬合反演的CAP方法求解,ML<4.0地震采用基于P波初動的振幅比方法求解。CAP反演震源機制的過程中,東段的地震采用魏文博等[8]的速度模型,西段的地震采用李祥等[9]的速度模型(圖2)。
2.1振幅比方法和CAP方法
振幅比計算過程中,運用層狀介質(zhì)中一點源位錯模型計算綜合地震圖得到P、S最大振幅,并對P、S最大振幅的比值與觀測數(shù)據(jù)進行擬合來確定震源機制。為保證震源機制解的穩(wěn)定性和可靠性,實際操作中約定參與計算的臺站數(shù)目為5個以上,振幅讀取精度控制在0.2 mm以內(nèi),記錄振幅大于0.5 mm且不限幅。同時,采用限定震中距(150 km以內(nèi))的辦法排除Pn、P11、Sn、S11等震相。
根據(jù)CAP方法和相關約束原則[7],挑選震中距在150~500 km范圍內(nèi)地震事件的EVT格式波形記錄,首先將臺站記錄的儀器響應去除,把地震波形旋轉(zhuǎn)至徑向、切向和垂向,分解為Pnl和Snl兩部分,用4階帶通濾波分別在0.05~0.2 Hz和0.05~0.1 Hz頻段內(nèi)對Pnl和Snl的噪聲部分進行壓制并賦予不同的權重,然后采用頻率-波數(shù)方法(F-K)計算震中距的格林函數(shù)并得到合成地震圖,與實際觀測波形進行互相關,最終通過格點搜索方法進行震源機制反演求解。
2.2計算結(jié)果
2001年以來河套地震帶224次ML≥2.8地震事件震源機制結(jié)果見圖3。圖4給出了CAP方法求解震源機制的示例。對于波形擬合效果較好的地震(圖4(a)),可以通過反演不同矩心深度下的震源機制解(圖4(c))獲取最佳震源深度和最佳雙力偶解(圖4(b))。
圖4(a)為理論地震圖(粗灰線)和觀測地震圖(細黑線)的波形擬合效果,每行波形最左側(cè)為臺站代碼和震中距,擬合曲線下方第1行數(shù)字為理論地震圖相對觀測地震圖的移動時間,第2行數(shù)字為理論地震圖和觀測地震圖的相關系數(shù),PnlV、PnlR分別為Pnl垂直分量和徑向分量。圖4(b)為震源深度為4 km時的最佳雙力偶解(下半球投影),黑色為拉張區(qū),白色為壓縮區(qū);圖4(c)為不同震源深度下的擬合殘差分布,其中深度為4 km時對應的殘差值最小。
為整體把握震源機制解所表達的應力場方向特征,分別以10°、20°和30°為基本統(tǒng)計單位,對震源機制參數(shù)進行歸一頻數(shù)統(tǒng)計和累計比率分析(圖5)。結(jié)果顯示,在研究時段內(nèi),斷層的兩個節(jié)面走向分別在北東至北東東向和北西向存在優(yōu)勢分布,根據(jù)河套地震帶所轄構(gòu)造的分布特點,北東向節(jié)面更能代表應力場方向特征;節(jié)面的傾角分布在30°~80°,覆蓋范圍較大;滑動角在±180°和0°附近存在優(yōu)勢分布,契合研究區(qū)以水平走滑為主的斷層作用機制特征;壓應力方位整體分布呈北東向,張應力T軸方位整體分布呈北西向,同樣符合應力場背景特征。
3震源機制類型的時空分布
在震源機制參數(shù)中,斷層節(jié)面的滑動角最易表現(xiàn)其運動情形。為考察河套地震帶的震源機制類型構(gòu)成情況,依據(jù)滑動角和震源機制類型的對應關系[9],對震源機制類型進行歸類統(tǒng)計。其中,斜滑型占46.2%、走滑型占33.3%、逆沖型占12.5%、正斷層占8.0%,斜滑型和純走滑型共占79.5%。可以看出,研究時段內(nèi)河套地震帶的構(gòu)造應力張量結(jié)構(gòu)以走滑為主。
不同的震源機制類型對應不同的應力作用模式,其中正斷層對應拉張應力,逆斷層對應擠壓應力,走滑斷層對應水平剪切應力[10],因此,震源機制類型的時序分布也可反映應力作用方式的時序變化形態(tài)。從圖6看出,2001年以來,河套地震帶主要產(chǎn)出走滑型和斜滑型地震,正斷層和逆沖型震源機制類型分布較為稀疏;從空間分布看(圖7),除在盆地向隆起過渡地段出現(xiàn)少量的正斷層和逆沖型震源機制外,河套地震帶內(nèi)地震主要受到水平剪切作用而多表現(xiàn)為走滑型震源機制。
4結(jié)語
震源機制類型構(gòu)成結(jié)果顯示,2001年以來,河套地震帶的斜滑型地震占46.2%、走滑型地震占33.3%、逆沖型地震占12.5%、正斷層地震占8.0%,走滑型和斜滑型地震總比例達到79.5%。歸一頻度統(tǒng)計結(jié)果顯示,2001年以來,河套地震帶的壓應力方位在北東向存在優(yōu)勢分布,張應力方位在北西向存在優(yōu)勢分布,滑動角主要分布在-30°~30°,即水平剪切拉張運動總體控制了斷裂的作用模式而多產(chǎn)出走滑型地震。結(jié)合外圍動力學環(huán)境和區(qū)域構(gòu)造運動性質(zhì)分析認為,在太平洋板塊向西俯沖和青藏塊體向北運移的聯(lián)合作用下,河套地震帶兩側(cè)的陰山隆起和鄂爾多斯塊體相對作旋性扭動,驅(qū)使北東向-近東西向展布的斷裂產(chǎn)生水平旋性剪切而主要產(chǎn)出正走滑型震源機制。另外,河套地震帶西段應力場的壓應力方向為北北東-北東向,東段應力場的壓應力方位為北東東-近東西向,應力場方向的區(qū)域性偏轉(zhuǎn)現(xiàn)象應該是河套地震帶應力場長期以來的一個典型特征[2]。除了受外圍動力學作用形式以及構(gòu)造塊體的運動模式約束外,其還可能與西段的斷裂構(gòu)造多呈北東向展布、東段的斷裂構(gòu)造多為近東西向展布有關(圖1)。斷裂的展布形態(tài)也影響到了震源機制斷層節(jié)面的方位分布,使其產(chǎn)生構(gòu)造相依的分布特點。
致謝:感謝美國圣路易斯大學朱露培教授提供CAP程序。
參考文獻
[1]國家地震局《鄂爾多斯周緣活動斷裂系》課題組.鄂爾多斯周緣活動斷裂系[M].北京:地震出版社,1988(The Research Group on “Active Fault System around Ordos Massif”. Active Fault System around Ordos Massif[M].Beijing:Seismological Press,1988)
[2]謝富仁,陳群策,崔效鋒,等.中國大陸地殼應力環(huán)境研究[M].北京:地質(zhì)出版社,2003(Xie Furen,Chen Qunce,Cui Xiaofeng,et al.Research on Crustal Stress State in China and Adjacent Area[M].Beijing: Geological Publish House,2003)
[3]聶宗笙,任云,劉志明,等.內(nèi)蒙古包頭市區(qū)大青山山前斷裂地震活動斷層初步研究[J].現(xiàn)代地質(zhì),2011,25(5):938-957(Nie Zongsheng,Ren Yun,Liu Zhiming,et al.Preliminary Study on Seismic Active Faults of the Daqingshan Frontal Fault in Baotou City,Inner Mongolia[J].Geoscience,2011,25(5):938-957)
[4]聶宗笙.公元前7年內(nèi)蒙古包頭地區(qū)8級地震的初步研究[J].地震學報,2013,35(4):584-603(Nie Zongsheng.Preliminary Investigation on the Historical M8 Earthquake Occurred in 7 BC at Baotou,Inner Mongolia[J].Acta Seismologica Sinica,2013,35(4):584-603)
[5]粱尚鴻,李幼銘,束沛鎰,等.利用區(qū)域地震臺網(wǎng)P、S波振幅比資料測定小震震源參數(shù)[J].地球物理學報,1984,27(3):249-257(Liang Shanghong,Li Youming,Shu Peiyi,et al.On the Determining of Source Parameters of Small Earthquakes by Using Amplitude Ratios of P and S from Regional Network Observations[J].Chinese J Geophys,1984,27(3):249-257)
[6]刁桂苓,于利明,寧杰遠,等.1989年大同震群的破裂特征[J].地球物理學報,1992,36(3):360-368(Diao Guiling,Yu Liming,Ning Jieyuan,et al.Fracture Characteristic of the 1989 Datong Earthquake Swarm[J].Chinese J Geophys,1992,36(3):360-368)
[7]Zhu L P,Helmberger D V.Advancement in Source Estimation Techniques Using Broadband Regional Seismograms[J].Bull Seism Soc Amer,1996,86(5):1 634-1 641
[8]李祥,包東健,陳元平,等.內(nèi)蒙古中部地區(qū)的殼幔介質(zhì)分層結(jié)構(gòu)[J].華北地震科學,1987,5(4):58-65(Li Xiang,Bao Dongjian,Chen Yuanping,et al.The Layered Structure of the Curst and Mantle Medium in Central Inner Mongolia[J].North China Earthquake Sciences,1987,5(4):58-65)
[9]陳運泰,顧浩鼎.震源理論基礎[M].北京:地震出版社,2008(Chen Yuntai,Gu Haoding.Basis of Focal Theory[M].Beijing: Seismological Press, 2008)
[10]刁桂苓,徐錫偉,陳于高,等.汶川MW7.9和集集MW7.6地震前應力場轉(zhuǎn)換現(xiàn)象及其可能的前兆意義[J].地球物理學報,2011,54(1):128-136(Diao Guiling,Xu Xiwei,Chen Yugao,et al.The Precursory Significance of Tectonic Stress Field Transformation before the Wenchuan MW7.9 Earthquake and the Chi-Chi MW7.6 Earthquake[J].Chinese J Geophys,2011, 54(1):128-136)
[11]張致偉,周龍泉,龍峰,等.汶川8.0和蘆山7.0級地震序列應力場時空特征[J].地震地質(zhì),2015,37(3):804-817(ZhangZhiwei,ZhouLongquan,LongFeng,etal.SpatialandTemporalCharacateristicofStressFieldforWenchuanMS8.0andLushanMS7.0EarthquakeSequence[J].SeismologyandGeology, 2015, 37(3):804-817)
FocalMechanismDeterminationofMedium-SmallEarthquakesand
Foundationsupport:TheSparkProgramofEarthquakeTechnologyofCEA,No.XH15010Y.
Aboutthefirstauthor:HANXiaoming,associateresearcher,majorsindigitalseismology,E-mail:hxmpower@126.com.
CharacteristicofStressFieldinHetaoSeismicBelt
HAN Xiaoming1LIU Fang1ZHANG Fan1WEI Jianmin1XU Yan1HOU Di1LI Juan1LIANG Shasha1
1EarthquakeAdministrationofInnerMongolia, 80ZhelimuRoad,Hohhot010010,China
Abstract:InordertoinspectthechangeformofschematransformationandthestressformofHetaoseismicstressfield,weuseamplituderatioandCAPmethodstocalculate224focalmechanismsolutions,basedonnormalizedfrequencystatisticsandanalysisoffaultsliptype.TheresultsshowthattheadvantagedistributionofcompressivestressorientationisNE,theadvantagedistributionoftensilestressorientationisNW,andtherakeangleisdistributednearlyhorizontally.Further,thewholestressfieldshowsinheritancefeatures,horizontalshearandtensioncontrolthestressformoffaultsinHetaoseismicbelt,andthemainoutputsarestrike-sliptypeearthquakes.Throughoutthewholestudyperiod,asteadystressmodelispresented.Combiningperipheraldynamicpropertieswithregionaltectonicmovementdeterminesthat,underthejointactionofthewestwardsubductionofthePacificplateandtheQinghai-Tibetblocknorthwardmigration,YinshanupliftandOrdosblockrelativetodospintwist.Finally,theNE-EWdistributionfaultspresenthorizontalandrotaryshear,andconstitutearelativelystableevolutionmodelofstressfield.
Keywords:Hetaoseismicbelt;focalmechanismsolution;amplituderatio;CAP;characteristicsofstressfield
收稿日期:2015-06-26
第一作者簡介:韓曉明,副研究員,主要從事數(shù)字地震學研究,E-mail:hxmpower@126.com。
DOI:10.14075/j.jgg.2016.06.006
文章編號:1671-5942(2016)06-0495-04
中圖分類號:P315
文獻標識碼:A
項目來源:中國地震局地震科技星火計劃(XH15010Y)。