亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        福建仙游地震序列的震源機(jī)制解

        2014-09-04 07:43:32邱毅李軍康蘭池袁麗文
        中國(guó)地震 2014年2期
        關(guān)鍵詞:仙游震區(qū)張量

        邱毅 李軍 康蘭池 袁麗文

        福建省地震局,福州市華鴻路7號(hào) 350003

        0 引言

        福建仙游自2010年8月以來(lái)陸續(xù)出現(xiàn)小震活動(dòng),截至2013年10月4日共記錄到地震1209次,其中5級(jí)地震1次,4級(jí)地震3次,3級(jí)地震11次,最大震級(jí)為ML5.0。該區(qū)域是2012、2013年福建地區(qū)地震最活躍的區(qū)域,多次發(fā)生有感地震。2013年9月4日福建仙游ML5.0地震也是福建地區(qū)近16年來(lái)發(fā)生的最大地震,地震發(fā)生時(shí)福建全省多處有明顯震感,震中區(qū)域震感強(qiáng)烈,震中區(qū)烈度達(dá)到Ⅵ度,出現(xiàn)房屋倒塌、墻壁開裂等災(zāi)情。因而開展對(duì)仙游地震序列的震源特征研究,進(jìn)一步了解發(fā)震構(gòu)造及震區(qū)區(qū)域應(yīng)力場(chǎng)具有一定的科學(xué)意義。

        震源機(jī)制解對(duì)于震源過(guò)程研究和地震預(yù)測(cè)研究具有重要意義,是研究現(xiàn)代構(gòu)造應(yīng)力場(chǎng)的基礎(chǔ)資料之一。通過(guò)震源機(jī)制研究,可以分析地震發(fā)生的力學(xué)過(guò)程,推斷地震的主要成因(趙翠萍等,2008;王勤彩等,2009)。P波初動(dòng)法和地震矩張量反演方法是目前常用的震源機(jī)制反演手段。美國(guó)加州大學(xué) Dreger等(Dreger et al,1993、1999;Dreger,2002)從1993年就開始開發(fā)利用區(qū)域臺(tái)網(wǎng)寬頻帶地震資料實(shí)現(xiàn)矩張量反演的程序,提出了利用區(qū)域Pnl波進(jìn)行地震矩張量反演的方法。本文基于此方法,利用福建數(shù)字地震臺(tái)網(wǎng)寬頻帶地震波形記錄反演得到福建仙游地震序列ML>3.5的6次地震的震源機(jī)制解。

        1 計(jì)算原理

        Derger等(1993、1999)提出了利用區(qū)域臺(tái)網(wǎng)寬頻帶臺(tái)站長(zhǎng)周期地震記錄來(lái)反演中等強(qiáng)度地震矩張量的方法。此后,Derger(2002)對(duì)該方法作了改進(jìn),使用Saikai(1994)改進(jìn)的離散波數(shù)積分方法計(jì)算理論地震波形,并用該方法反演了加州北部M3.5以上地震的矩張量。

        Saikai(1994)改進(jìn)的離散波數(shù)積分方法基于Wang等(1980)和Herrmann等(1985)給出的層狀介質(zhì)中任意一個(gè)位錯(cuò)點(diǎn)源和爆炸源所產(chǎn)生的地震波場(chǎng)表達(dá)式及所需要的10個(gè)格林函數(shù)。該表達(dá)式在圓柱坐標(biāo)系下,傅立葉變換后在自由表面距震源r處的垂直、徑向和切向位移可分別表示為

        式中,r是臺(tái)站到震源的距離;α為臺(tái)站相對(duì)震中位置的方位角;θ、λ、δ分別是斷層的走向、滑動(dòng)角和傾角;ZSS、ZDS、ZDD分別相當(dāng)于 1 個(gè)直立的純走滑斷層(θ=0°,λ =0°,δ=90°)、1 個(gè)直立的純逆斷層(θ=0°,λ =90°,δ=90°) 以及一個(gè) 45°的斜滑斷層(θ=0°,λ =90°,δ=45°)在坐標(biāo)為(r,α,0)點(diǎn)產(chǎn)生的垂直分量位移;RSS、RDS、RDD相當(dāng)于上述3種斷層產(chǎn)生的徑向位移;TSS、TDS為直立的走滑斷層和逆斷層產(chǎn)生的切向分量位移;ZEP和REP為一個(gè)爆炸源產(chǎn)生的垂直、徑向位移。

        若僅考慮一個(gè)純偏矩張量震源,有ZEP=0,REP=0,Mzz=-(Mxx+Myy),上式中只包含5個(gè)矩張量元素。矩張量的反演問(wèn)題可以用矩陣表示為

        式中,d為n行觀測(cè)位移矢量,G 為n×5階的格林函數(shù)矩陣,m 為包含Mxx、Myy、Mzz、Mxy、Myz這5個(gè)矩張量元素的矢量??梢姰?dāng)參加反演的觀測(cè)數(shù)據(jù)大于5個(gè)時(shí),反演問(wèn)題的解是超定問(wèn)題。采用最小二乘法進(jìn)行反演,通過(guò)考察目標(biāo)函數(shù)RES/Pdc(偏量與最佳雙力偶的比值)和VR來(lái)確定最佳解。VR為觀測(cè)波形和理論波形的擬合程度的參數(shù),VR值越大,表示殘差越小,波形擬合效果越好,求解出來(lái)的地震矩張量越可靠。RES/Pdc及VR的計(jì)算公式如下

        式中,data表示觀測(cè)波形數(shù)據(jù),synth表示理論波形數(shù)據(jù)。

        反演得到矩張量后,經(jīng)過(guò)坐標(biāo)變換,把計(jì)算得出的地震矩張量旋轉(zhuǎn)到主軸坐標(biāo)系下

        式中,M1、M2、M3是矩張量的本征值,M1對(duì)應(yīng)于張應(yīng)力軸(T),M2對(duì)應(yīng)于零軸(N),M3對(duì)應(yīng)于壓應(yīng)力軸(P)。它們的傾角和方位角以及斷層的2個(gè)節(jié)面的走向、傾角和滑動(dòng)角都可由P軸、T軸和N軸計(jì)算得到。

        2 數(shù)據(jù)處理和模型選取

        2.1 地震事件

        福建仙游地震序列共發(fā)生ML>3.5的地震6次,具體地震事件見表1。

        表1 福建仙游地震序列ML>3.5地震

        2.2 速度結(jié)構(gòu)模型

        速度模型的選取影響理論波形的生成,進(jìn)而影響震源機(jī)制解的反演結(jié)果和震源深度的確定,因此選取一個(gè)較為符合福建地區(qū)實(shí)際速度結(jié)構(gòu)的模型顯得較為重要。張路(2008)根據(jù)近20年來(lái)福建地區(qū)進(jìn)行的人工地震測(cè)深、大地電磁測(cè)深、重磁測(cè)量、地?zé)釡y(cè)量、地學(xué)大斷面的編制、天然地震的層析成像以及地球物理資料綜合分析等大量深部構(gòu)造研究工作,提出福建東南沿海地殼結(jié)構(gòu)。本文使用的福建地區(qū)速度結(jié)構(gòu)模型即在此基礎(chǔ)上給出的(表2)。

        表2 福建地區(qū)速度結(jié)構(gòu)模型

        2.3 數(shù)據(jù)選取與處理

        Derger(2002)方法基于區(qū)域臺(tái)網(wǎng)寬頻帶臺(tái)站長(zhǎng)周期地震記錄來(lái)反演中等強(qiáng)度地震矩張量。對(duì)于不同的震級(jí),要選取不同震中距、不同頻帶范圍、不同長(zhǎng)度的波形記錄(Kubo et al,2002),如表3所示。

        表3 震中距、頻帶范圍、數(shù)據(jù)長(zhǎng)度與震級(jí)關(guān)系(據(jù)Kubo等(2002)文獻(xiàn))

        福建數(shù)字地震臺(tái)網(wǎng)現(xiàn)有寬頻帶臺(tái)站70個(gè)(圖1)。其中23個(gè)臺(tái)的帶寬為50Hz/60s,如圖1中空心三角形所示,47個(gè)臺(tái)帶寬為50Hz/120s,如圖1實(shí)心三角形所示,震中四周臺(tái)站分布較好,這些臺(tái)為仙游地震震源機(jī)制解反演提供了豐富可靠的波形資料。福建仙游地震序列的6次ML>3.5地震的震級(jí)都小于等于5.0級(jí),考慮到若震中距太遠(yuǎn),波形數(shù)據(jù)信噪比會(huì)降低。因此,本文選取福建臺(tái)網(wǎng)震中距在50~400km內(nèi)的寬頻帶臺(tái)站;對(duì)這些臺(tái)三分量波形記錄進(jìn)行數(shù)據(jù)質(zhì)量檢測(cè),剔除信噪比差的臺(tái)站數(shù)據(jù);對(duì)剩下的臺(tái)站記錄進(jìn)行儀器響應(yīng)校正,將記錄變?yōu)閷?shí)際的速度值,再將儀器響應(yīng)校正后的實(shí)際速度記錄積分為位移,并旋轉(zhuǎn)到Z、R、T分量,使用Butterworth 4階濾波器濾波,濾波頻帶為0.02~0.05 Hz;采用前述2.2的速度結(jié)構(gòu)模型,利用Saikai(1994)改進(jìn)的離散波數(shù)積分方法計(jì)算不同深度、不同震中距的理論地震波形,采用與觀測(cè)地震波形記錄同樣的濾波范圍0.02~0.05 Hz對(duì)理論地震波形進(jìn)行濾波;由于Derger方法使用3臺(tái)三分量的波形記錄進(jìn)行反演,可以得到較為穩(wěn)定的結(jié)果(Dreger,2002;Kubo et al,2002),對(duì)剩下的臺(tái)站進(jìn)行不同深度、不同臺(tái)站組合的3臺(tái)組合反演,利用VR值作為第1目標(biāo)函數(shù),選取VR值最大的幾個(gè)反演結(jié)果,然后用RES/Pdc的值作為第2目標(biāo)函數(shù),對(duì)已選取VR值最大的幾個(gè)反演結(jié)果,選取其中RES/Pdc值最小的最佳解為地震的震源機(jī)制解。

        圖1 福建寬頻帶臺(tái)站分布

        3 計(jì)算結(jié)果與分析

        3.1 反演結(jié)果

        福建仙游2012年4月15日ML4.1、11月25日ML3.8、2013年8月3日ML4.2、8月19日ML3.8、8月23日ML4.5、9月4日ML5.0地震的反演結(jié)果如圖2、3和表4所示。

        3.2 結(jié)果分析

        福建仙游地震序列發(fā)生于福建省莆田市仙游縣,歷史上該區(qū)域?qū)儆谏僬稹⑷跽饏^(qū)。仙游震群自2010年8月開始活動(dòng),在2012、2013年發(fā)生多次有感地震,其中2013年9月4日的ML5.0地震為福建地區(qū)近16年來(lái)發(fā)生的最大地震。整個(gè)序列地震活動(dòng)呈現(xiàn)逐漸增強(qiáng)的趨勢(shì),地震的頻次和強(qiáng)度逐漸提高。

        表4 最佳雙力偶解(單位:°)

        本文利用區(qū)域臺(tái)網(wǎng)寬頻帶地震記錄,采用地震矩張量反演方法得到福建仙游2012年4月15日ML4.1、11月25日ML3.8、2013年8月3日ML4.2、8月19日ML3.8、8月23日ML4.5、9月4日ML5.0地震的震源機(jī)制解(圖3)。仙游震群發(fā)生在NW向的沙縣-南日島斷裂帶上,結(jié)合地震序列的震中分布圖(圖3),判定斷層走向?yàn)镹W向,斷層節(jié)面為節(jié)面B。通過(guò)反演結(jié)果可知,這6次地震的震源機(jī)制解具有較高的相似性,都為走滑型,斷層走向都為NW向,具有高傾角的性質(zhì),主壓應(yīng)力軸都為近SN向。

        通過(guò)震源機(jī)制解可以得知震區(qū)應(yīng)力狀態(tài)。福建地區(qū)地處歐亞板塊的東南邊緣地帶,區(qū)域應(yīng)力場(chǎng)主要受菲律賓板塊對(duì)歐亞板塊的擠壓碰撞,福建地區(qū)的現(xiàn)代區(qū)域應(yīng)力場(chǎng)為NW向(汪素云等,1985;謝富仁等,2004;丁學(xué)仁等,1999;林松建等,2009;袁麗文等,2009),近些年福建地區(qū)的中等強(qiáng)度地震震源機(jī)制解的主壓應(yīng)力軸也主要為NW向(袁麗文等,2009;邱毅等,2013)。本文反演得到的福建仙游地震序列的6次ML>3.5地震的主壓應(yīng)力軸都為近SN向,與福建地區(qū)現(xiàn)代區(qū)域應(yīng)力場(chǎng)并不一致。金鐘水庫(kù)位于仙游震區(qū),于2012年5月開始下閘蓄水,庫(kù)容1.06億m3,壩高97.5m,水壩離這6次地震最大的震中距為4km,最小的震中距不到1km。仙游地區(qū)歷史上屬于少震、弱震區(qū),但在金鐘水庫(kù)開始蓄水3個(gè)月后(2012年8月),該區(qū)出現(xiàn)震群活動(dòng),地震的頻次和強(qiáng)度逐漸提高,從2010年8月4日至2013年10月4日共記錄到地震1209次,其中ML5.0以上地震1次,ML4.0~4.9地震3次。本文認(rèn)為仙游地震是金鐘水庫(kù)蓄水引起仙游震區(qū)應(yīng)力調(diào)整而誘發(fā)的地震,所以仙游地震震源機(jī)制解主壓應(yīng)力軸方向與福建的現(xiàn)代構(gòu)造應(yīng)力場(chǎng)的方向不同。

        圖2 地震波形擬合及矩張量反演結(jié)果

        圖3 仙游地震序列的震中分布與ML>3.5地震的震源機(jī)制解

        4 結(jié)論

        (1)本文利用矩張量反演方法求解了福建仙游地震序列ML>3.5的6次地震的震源機(jī)制解。這6次地震發(fā)生在NW向的沙縣-南日島斷裂帶上,結(jié)合地震序列的震中分布圖,確定這6次地震的斷層走向?yàn)镹W向,發(fā)震斷層面為節(jié)面B。這6次地震的震源機(jī)制解具有較高的相似性,都為走滑斷層,高傾角,斷層走向都為 NW向,主壓應(yīng)力軸方向都為近 SN向,震區(qū)應(yīng)力場(chǎng)并未出現(xiàn)明顯調(diào)整變化的情況。

        (2)福建區(qū)域應(yīng)力場(chǎng)的方向?yàn)镹W向,而本文得到的仙游地震震源機(jī)制解主壓應(yīng)力軸方向?yàn)榻黃N向,兩者存在差異。本文認(rèn)為仙游震區(qū)的金鐘水庫(kù)蓄水引起震區(qū)應(yīng)力調(diào)整,進(jìn)而誘發(fā)了仙游地震。仙游地震是由震區(qū)的局部應(yīng)力場(chǎng)控制,故仙游震區(qū)的應(yīng)力場(chǎng)與福建地區(qū)的區(qū)域應(yīng)力場(chǎng)存在差異。

        致謝:本文使用的地震觀測(cè)波形記錄由福建省地震局監(jiān)測(cè)中心提供;使用的矩張量反演程序TDMT_INVC由美國(guó)BSL的Douglas S Dreger等開發(fā),其中格林函數(shù)計(jì)算使用的是URS的 Chandan Saikia開發(fā)的FKRPROG程序;兩位審稿專家對(duì)本文提出寶貴建議,在此表示謝意。

        猜你喜歡
        仙游震區(qū)張量
        流浪衛(wèi)星
        偶數(shù)階張量core逆的性質(zhì)和應(yīng)用
        四元數(shù)張量方程A*NX=B 的通解
        仙游百松村勁搖“清風(fēng)”
        紅土地(2018年12期)2018-04-29 09:16:52
        提高高中數(shù)學(xué)課堂效率的教學(xué)方法研究
        蘆山震區(qū)大田壩崩塌發(fā)育特征及其防治措施
        法律服務(wù)新窗口群眾訴求好幫手——仙游度尾鎮(zhèn)成立全市首家村級(jí)司法室
        紅土地(2016年11期)2017-01-15 13:46:25
        擴(kuò)散張量成像MRI 在CO中毒后遲發(fā)腦病中的應(yīng)用
        貢茶第一寨仙游伍家臺(tái)
        民族大家庭(2016年5期)2016-03-20 14:52:42
        工程中張量概念的思考
        河南科技(2014年19期)2014-02-27 14:15:33
        欧美a级在线现免费观看| 亚洲综合一区二区三区蜜臀av| 女同舌吻互慰一区二区| 成午夜福利人试看120秒| 亚洲精品美女久久久久99| 最新国产午夜福利| 手机在线免费av网址| 午夜福利一区在线观看中文字幕| 老司机亚洲精品影院| 欧美在线播放一区二区| 亚洲区一区二区三区四| 女优一区二区三区在线观看 | 国产亚洲2021成人乱码| 日日人人爽人人爽人人片av| 五月天亚洲av优女天堂| 亚洲精品国产av成人精品| 日本乱偷人妻中文字幕在线| 97人妻视频妓女网| 蜜桃成人精品一区二区三区| 少妇真人直播免费视频| 国产av无码专区亚洲av琪琪| 91精品综合久久久久m3u8 | 欧美熟妇另类久久久久久不卡| 日韩人妻精品无码一区二区三区 | 国产黑色丝袜在线观看下| 丰满人妻中文字幕乱码| 国产主播一区二区三区蜜桃| 无码av无码天堂资源网| 亚洲国产精品悠悠久久琪琪| 亚洲av男人的天堂在线| a级毛片免费观看在线播放| 乌克兰少妇xxxx做受6| 国产高跟丝袜在线诱惑| 国产亚洲自拍日本亚洲| 一本加勒比hezyo无码人妻| 春色成人在线一区av| 亚洲免费在线视频播放| 99爱在线精品免费观看| 香蕉视频毛片| 国产精品成年人毛片毛片| 亚洲av成人无码一区二区三区在线观看|