• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    遼寧蓋州海城地區(qū)地殼剪切波分裂特征

    2016-02-13 06:03:19
    地震地質(zhì) 2016年4期
    關(guān)鍵詞:海城營(yíng)口偏振

    李 金 高 原 張 博 王 亮 高 榮

    1)新疆維吾爾自治區(qū)地震局,烏魯木齊 830011 2)中國(guó)地震局地震預(yù)測(cè)研究所,北京 100036 3)遼寧省地震局,沈陽(yáng) 110034

    遼寧蓋州海城地區(qū)地殼剪切波分裂特征

    李 金1)高 原2)張 博3)王 亮3)高 榮1)

    1)新疆維吾爾自治區(qū)地震局,烏魯木齊 830011 2)中國(guó)地震局地震預(yù)測(cè)研究所,北京 100036 3)遼寧省地震局,沈陽(yáng) 110034

    2012年2月2日,在遼寧蓋州(40.56°N,122.36°E)發(fā)生ML4.7地震,此后該地區(qū)小震不斷,地震活動(dòng)持續(xù)至今仍未停止,發(fā)生多次4級(jí)以上地震。與此同時(shí),海城地震老震區(qū)小震持續(xù)活動(dòng),自蓋州震群活動(dòng)以來(lái),亦發(fā)生1,100多次ML≥1.0地震。基于這2叢地震活動(dòng),研究了蓋州-海城地區(qū)的剪切波分裂特征。初步研究結(jié)果表明,海城老震區(qū)快剪切波優(yōu)勢(shì)偏振方向較為穩(wěn)定,與該地區(qū)區(qū)域應(yīng)力場(chǎng)方向基本一致。由于蓋縣臺(tái)下方存在活動(dòng)斷裂,其快剪切波的優(yōu)勢(shì)偏振方向較為復(fù)雜,存在2個(gè)優(yōu)勢(shì)偏振方向,分別與臺(tái)站下方的金州斷裂走向和該區(qū)最大主應(yīng)力方向一致。此外,2013年12月22日之后,蓋州震群活動(dòng)突然增強(qiáng),與此同時(shí),自該時(shí)間之后,蓋縣臺(tái)(GAX)所反映出的快剪切波優(yōu)勢(shì)偏振方向?yàn)镾EE向,與營(yíng)口臺(tái)(YKO)的快剪切波優(yōu)勢(shì)方向基本接近,同時(shí)與該區(qū)最大主應(yīng)力方位較為一致。由此可以推斷,蓋州震群自2013年12月22日的這組增強(qiáng)活動(dòng)或與區(qū)域應(yīng)力的局部增強(qiáng)有關(guān)。但這2個(gè)震群的慢波時(shí)間延遲變化相對(duì)較小,與以往相關(guān)研究中強(qiáng)震前慢波時(shí)間延遲的變化差異較大。

    蓋州震群 剪切波分裂 快波偏振方向

    0 引言

    2012年2月2日,在遼寧蓋州(40.56°N,122.36°E)發(fā)生ML4.7地震,截至2014年10月30日,蓋州地區(qū)共發(fā)生ML≥1.0地震1,233次,其中ML1.0~1.9地震934次,ML2.0~2.9地震247次,ML3.0~3.9地震45次,ML≥4.0地震7次。7次≥4級(jí)地震分別為2012年2月2日ML4.7、ML4.2,7月12日ML4.2,11月1日ML4.0,2013年3月30日ML4.0,2014年6月26日ML4.1以及8月22日ML4.3,其地震活動(dòng)的表現(xiàn)形式為震群。與此同時(shí),1975年發(fā)生7.3級(jí)地震的海城地震老震區(qū)小震持續(xù)活動(dòng),自蓋州震群活動(dòng)以來(lái),亦發(fā)生1,100多次ML≥1.0地震。該區(qū)距蓋州震群不遠(yuǎn)。由于這2組地震活動(dòng)展示了遼寧地區(qū)近期地震活動(dòng)的新特點(diǎn),因而本研究利用海城地震老震區(qū)的營(yíng)口臺(tái)(YKO)記錄到的波形資料和蓋縣臺(tái)(GAX)記錄到的蓋州震群波形資料,開(kāi)展2個(gè)臺(tái)站剪切波分裂參數(shù)變化特征研究。

    大量觀測(cè)研究表明,地殼介質(zhì)中廣泛存在各向異性特征,其主要是由大量充滿(mǎn)流體的微裂隙的定向排列引起的,剪切波穿過(guò)這種含有微裂隙的介質(zhì)時(shí)會(huì)產(chǎn)生分裂現(xiàn)象,其中速度較快的波列稱(chēng)為快剪切波,速度較慢的波列稱(chēng)為慢剪切波(Crampin,1978,1981)??旒羟胁ǖ钠穹较蚺c定向排列的直立裂隙面平行,與原地最大主應(yīng)力方向一致,慢剪切波的偏振方向則近似垂直于裂隙面(Zatsepinetal.,1997; Gaoetal.,1998,2004; Crampinetal.,2002)。由于快剪切波的偏振方向?qū)^(qū)域主應(yīng)力環(huán)境的變化非常敏感,因此,剪切波分裂參數(shù)能夠較好地反映出區(qū)域應(yīng)力狀態(tài)及當(dāng)?shù)氐木植繕?gòu)造特征,尤其是對(duì)分析小尺度應(yīng)力環(huán)境特征非常有效(高原等,1995,1999;賴(lài)院根等,2006;鄭秀芬等,2008)。

    圖1 2012年以來(lái)遼寧蓋州海城地區(qū)ML1.5以上地震震中分布圖及歷史強(qiáng)震(據(jù)遼寧區(qū)域地震臺(tái)網(wǎng)正式目錄)Fig. 1 The distribution of Gaizhou earthquake swarm with magnitude larger than ML1.5 from 2012 and historical strong earthquake(after Liaoning regional seismic network catalog).圓圈表示地震: 紅色 ML1.5~1.9; 粉色 ML2.0~2.9; 黃色 ML3.0~3.9; 藍(lán)色 ML4.0~4.9

    1 構(gòu)造背景與資料

    蓋州震群和海城地震老震區(qū)附近有3條主要的活動(dòng)斷裂,分別為NE向的海城-營(yíng)口斷裂、 NNE向的金州斷裂以及NW向的海城河斷裂(圖1)。海城-營(yíng)口斷裂以及金州斷裂都屬于中國(guó)東部最大的斷裂帶,即郯廬斷裂。金州斷裂為正斷層,在第四紀(jì)時(shí)期表現(xiàn)出明顯的活動(dòng)性,走向 N15°~30°E,傾向NW,傾角40°~80°,為遼南地區(qū)主要的控震斷裂,但實(shí)際的地震破裂多發(fā)生在NW向斷裂上(雷清清等,2008)。大洋河斷裂則是1975年海城地震和1999年岫巖地震的發(fā)震斷裂。NW向的大洋河斷裂是1條正在形成和發(fā)育的斷裂,該斷裂從營(yíng)口高坎起,往東南經(jīng)岫巖偏嶺鎮(zhèn)延伸到東港市境內(nèi),走向?yàn)?20°(呂培苓等,2004)。金州斷裂主要展布在普蘭店西側(cè),在斷裂附近有一系列相關(guān)聯(lián)、 近平行的小型分支斷裂發(fā)育(鐘以章等,1999)。依據(jù)震源機(jī)制解、 小震綜合斷層面解、 原地應(yīng)力測(cè)量等資料,已經(jīng)確定遼東半島處在NEE或近EW向水平主應(yīng)力場(chǎng)作用下(高常波,1998; 徐杰等,1999,2012; 葛榮峰,2009)。2012年蓋州ML4.7地震距離1975年2月4日海城7.3級(jí)地震約37km,距離1999年11月29日岫巖5.9級(jí)地震約60km(圖1)。

    2 研究方法

    本研究采用 “剪切波分裂系統(tǒng)分析法”——即SAM方法(高原等,2004a),以相關(guān)函數(shù)法為基礎(chǔ),包括3個(gè)部分: 相關(guān)函數(shù)計(jì)算、 時(shí)間延遲校正和偏振分析檢驗(yàn); 該方法可以自我檢驗(yàn),提高了分析結(jié)果的準(zhǔn)確性。

    使用地表臺(tái)站記錄波形分析剪切波時(shí),要求地震射線(xiàn)滿(mǎn)足剪切波窗口的約束條件。剪切波入射到地表時(shí),若入射角大于臨界角會(huì)發(fā)生全反射現(xiàn)象。為了避免發(fā)生全反射現(xiàn)象,在挑選數(shù)據(jù)時(shí)應(yīng)選擇剪切波窗口內(nèi)的數(shù)據(jù)。本研究針對(duì)遼寧地區(qū)的一維速度模型查閱了多份資料,其中盧造勛等(2002)給出了遼寧地區(qū)的速度模型(表1)。由于蓋州震群大多發(fā)生于6~9km深度范圍內(nèi)(曹鳳娟等,2013),本研究在盧造勛(2002)速度模型的基礎(chǔ)上對(duì)該區(qū)速度模型進(jìn)行了適當(dāng)?shù)暮?jiǎn)化,將該區(qū)地殼分為3層,層厚分別為4km、 9km、 11km,每層P波速度分別為 4.5km/s、 5.8km/s、 6.1km/s,平均地幔速度設(shè)為 7.8km/s。在理論上,對(duì)于泊松介質(zhì)(泊松比為0.25),臨界角約為35°,由于低速地表沉積層的存在,實(shí)際地震射線(xiàn)可放寬至45°~50°(Crampin,2005)。本研究中,蓋州—海城地區(qū)沉積層較厚,而蓋州震群和海城老震區(qū)地震的深度多為6~9km,其沉積層的影響仍較為顯著,因此本研究刪除了入射角>50°的地震。由此確定了蓋州震群附近營(yíng)口臺(tái)(YKO)、 蓋縣臺(tái)(GAX)剪切波窗口內(nèi)的地震數(shù)據(jù)。

    表1 盧造勛(2002)給出的遼寧地區(qū)一維速度模型

    Table1 One dimensional velocity model in Liaoning Province (after Lu Zaoxun,2002)

    深度/km-1041324333571P波速度/km·s-14.24.55.86.16.46.87.88.0

    圖2 分別顯示了對(duì)營(yíng)口臺(tái)(YKO)(a)和蓋縣臺(tái)(GAX)(b)記錄到的1次地震事件使用SAM方法進(jìn)行剪切波分裂的實(shí)例。 圖2a地震為2013年1月16日發(fā)生在營(yíng)口臺(tái)附近的1.6級(jí)地震,震中距為6.23km,震源深度為6.0km; 圖2b地震為2014年5月11日發(fā)生在蓋縣臺(tái)附近的2.9級(jí)地震,震中距為7.91km,震源深度為10km。使用SAM軟件對(duì)原始波形數(shù)據(jù)進(jìn)行濾波處理并選取含完整S波列的SN分量(a3,b3)和EW分量(a4,b4),產(chǎn)生剪切波質(zhì)點(diǎn)運(yùn)動(dòng)軌跡(a1,b1)。SAM方法最大的優(yōu)點(diǎn)是可以自我檢驗(yàn),將慢剪切波(a6,b6)按得到的延遲值前移,本例的地震計(jì)采樣率為100采樣點(diǎn)/s,則1個(gè)采樣點(diǎn)為0.01s。慢剪切波分別提前0.03s(a)、 0.04s(b),剪切波質(zhì)點(diǎn)運(yùn)動(dòng)軌跡(a2,b2)顯示出線(xiàn)性關(guān)系。

    圖2 YKO(營(yíng)口臺(tái))記錄的數(shù)據(jù)(201301161708)(a)和GAX(蓋縣臺(tái))記錄的數(shù)據(jù)(201405111057)(b)進(jìn)行的剪切波分裂分析過(guò)程Fig. 2 Shear-wave splitting analysis of seismic wave(201301161708,201405111057)recorded by station YKO(a) and GAX(b).a1 為(YKO)營(yíng)口臺(tái)SN和EW方向分量的剪切波的質(zhì)點(diǎn)運(yùn)動(dòng)軌跡以及剪切波波形(a3、 a4),a1中的S1和S2分別表示快剪切波和慢剪切波的到時(shí); a2為已消除了時(shí)間延遲影響后快(F)和慢(S)剪切波的質(zhì)點(diǎn)運(yùn)動(dòng)軌跡以及快、 慢剪切波的波形(a5、 a6),其中橫坐標(biāo)表示采樣點(diǎn)的個(gè)數(shù),縱坐標(biāo)表示振幅count值,豎線(xiàn)表示用來(lái)做質(zhì)點(diǎn)運(yùn)動(dòng)軌跡的剪切波部分。b圖GAX臺(tái)的數(shù)據(jù)處理過(guò)程與YKO相同; a1、 a2、 b1、 b2的縱、 橫坐標(biāo)都表示count值

    表2 蓋縣臺(tái)地震計(jì)方位角計(jì)算結(jié)果

    Table 2 Azimuth calculation results of seismometer of GAX

    時(shí) 段地震計(jì)方位角/(°)2008-05-08前02008-05-08—2012-12-31-792013-01-01—2013-03-141012013-03-14—2013-10-141802013-10-15至今0

    3 研究結(jié)果

    通過(guò)對(duì)營(yíng)口臺(tái)(YKO)和蓋縣臺(tái)(GAX)波形資料分析,得到了營(yíng)口臺(tái)(YKO)、 蓋縣臺(tái)(GAX)的快波偏振方向,其有效數(shù)據(jù)個(gè)數(shù)分別為154和27個(gè)。在此尤其需要提及的是,在計(jì)算蓋縣臺(tái)(GAX)剪切波分裂時(shí),發(fā)現(xiàn)該臺(tái)站極性存在偏轉(zhuǎn)(地震計(jì)的N方向與真N方向存在夾角),參考鄭秀芬等*鄭秀芬等,2015,遼寧蓋縣臺(tái)地震計(jì)方位角校正過(guò)程(內(nèi)部報(bào)告)。的計(jì)算結(jié)果,其具體偏轉(zhuǎn)角度見(jiàn)表2,因此在計(jì)算蓋縣臺(tái)(GAX)剪切波分裂之前,首先對(duì)該臺(tái)2012—2015年底的數(shù)據(jù)進(jìn)行了方位角校正。

    圖3 YKO臺(tái)(a)、 GAX臺(tái)(b)的快剪切波偏振方向等面積極射投影與等面積投影玫瑰圖Fig. 3 The homolographic projection and homolographic projection rose graphs of the polarizations of the fast shear-wave splittings at station YKO(a)and GAX(b).a YKO(154)為營(yíng)口臺(tái)全部154個(gè)有效數(shù)據(jù)結(jié)果,b GAX(27)為蓋縣臺(tái)全部27個(gè)有效數(shù)據(jù)的結(jié)果;短線(xiàn)段的方向是地震臺(tái)記錄的每個(gè)地震事件的快剪切波偏振方向,外圓的半徑是剪切波窗口

    營(yíng)口臺(tái)(YKO)位于NE向的海城-營(yíng)口活動(dòng)斷裂東側(cè),附近還有NNE向的金州斷裂以及NW向的大洋河斷裂。根據(jù)營(yíng)口臺(tái)(YKO)的等面積投影玫瑰圖(圖3a)顯示,營(yíng)口臺(tái)(YKO)快波偏振的優(yōu)勢(shì)方向?yàn)镹EE-SWW向。太齡雪等(2009)使用遼寧區(qū)域地震臺(tái)網(wǎng)1999年6月至2004年12月的地震記錄,研究了遼寧地區(qū)的剪切波分裂特征,得到的營(yíng)口臺(tái)(YKO)的快剪切波平均偏振方向?yàn)镹77.6°E±26°,而1999年11月29日遼寧岫巖MS5.9地震前后營(yíng)口臺(tái)(YKO)的快剪切波優(yōu)勢(shì)偏振方向從NEE-SWW向變?yōu)榻麰W方向(太齡雪等,2008),該結(jié)果與本研究結(jié)果有很好的一致性。張萍等(2001a,b)利用岫巖地震序列得到該區(qū)的最大主應(yīng)力方位與本研究得到的營(yíng)口臺(tái)(YKO)快剪切波優(yōu)勢(shì)偏振方向比較一致,也與華北北部近EW 向的構(gòu)造應(yīng)力場(chǎng)方向較為一致。而曹鳳娟等(2005)對(duì)岫巖地震序列中營(yíng)口臺(tái)(YKO)記錄的1999年11—12月的26個(gè)數(shù)據(jù)進(jìn)行剪切波分裂進(jìn)行分析,得出了震前快剪切波平均偏振方向?yàn)镹58.7°E,震后快剪切波優(yōu)勢(shì)方向有2個(gè),其中1個(gè)平均偏振方向?yàn)镹68°E。最近根據(jù)剪切波分裂結(jié)果得到華北地區(qū)(至少在華北北部)的最大主應(yīng)力方向約為ENE(Gaoetal.,2011),而在考慮了盡量去除斷裂等因素影響后的區(qū)域背景壓應(yīng)力方向?yàn)镹95°E ±15°,近似為EW向(高原等,2010)。

    此外,由圖1 可以看出,營(yíng)口臺(tái)(YKO)附近的小震活動(dòng)空間分布近EW向。此前關(guān)于1999年岫巖MS5.9地震的相關(guān)研究也表明,岫巖地震序列的空間分布呈EW向(太齡雪,2008)。有研究表明,位于活動(dòng)斷層上的地震臺(tái)站得到的快剪切波優(yōu)勢(shì)偏振方向與斷層走向一致(石玉濤等,2006; 吳晶等,2007)。進(jìn)一步的研究揭示快剪切波優(yōu)勢(shì)偏振方向可用來(lái)判斷地下斷層是否為活動(dòng)斷層,并可能尋找隱伏活動(dòng)斷層(Gaoetal.,2011)。營(yíng)口臺(tái)(YKO)位于NE走向的海城-營(yíng)口斷裂東側(cè),在NE向金州斷裂的NE方向的前方,似乎并未直接處于大斷裂帶上,由該臺(tái)站附近小震空間分布可以顯示出臺(tái)站下方可能存在近EW向的活動(dòng)斷裂,營(yíng)口臺(tái)(YKO)的快波優(yōu)勢(shì)偏振方向可能與此條斷裂有關(guān)。

    由圖1 可以看出,蓋縣臺(tái)(GAX)位于NE向的金州斷裂的上方。蓋縣臺(tái)(GAX)的等面積投影玫瑰圖(圖3b)顯示,蓋縣臺(tái)(GAX)快波偏振優(yōu)勢(shì)方向?yàn)镹NE向和SEE向。SEE向的優(yōu)勢(shì)方向與當(dāng)?shù)刈畲笾鲬?yīng)力方向較為一致。曹鳳娟等(2013)對(duì)蓋縣臺(tái)(GAX)附近蓋州震群定位結(jié)果顯示,蓋州震群震中分布近似NWW-SEE向,與本研究蓋縣臺(tái)(GAX)的快剪切波NWW-SEE向優(yōu)勢(shì)偏振方向較為一致,由此或可推斷,該臺(tái)站下方也可能存在1條近EW向的隱伏斷裂。此外NNE向的優(yōu)勢(shì)偏振方向和臺(tái)站下方的金州斷裂走向較為接近,根據(jù)對(duì)較大規(guī)模的活動(dòng)斷裂的研究結(jié)果,活動(dòng)走滑斷裂上的臺(tái)站快剪切波偏振方向與周?chē)_(tái)站的結(jié)果不一致,往往與斷裂走向一致(Crampinetal.,2003; 吳晶等,2007; Gaoetal.,2011)。蓋縣臺(tái)(GAX)的研究結(jié)果支持了這個(gè)結(jié)論。

    圖4 海城老震區(qū)(a)和蓋州震群(b)m-t圖Fig. 4 The m-t of the previous Haicheng earthquake area(a)and Gaizhou earthquake swarm(b).a,b中的黑色箭頭為圖5 中的時(shí)間分段點(diǎn)

    分析海城老震區(qū)小震和蓋州震群的m-t圖可以看出(圖4),海城老震區(qū)ML≥1.0地震自2012年以來(lái)活動(dòng)頻度變化并不明顯,且該震群活動(dòng)強(qiáng)度不大,最大震級(jí)僅為ML3.4。而蓋州震群2012—2013年12月21日地震頻度相對(duì)較低,但發(fā)生了5次ML≥4.0以上地震,活動(dòng)強(qiáng)度相對(duì)較高; 2013年12月22日以來(lái),蓋州震群的地震頻度明顯升高,后期又發(fā)生2次ML≥4.0地震。本研究以2013年12月22日為時(shí)間節(jié)點(diǎn),分別研究此時(shí)間節(jié)點(diǎn)前(時(shí)段I)、 后(時(shí)段Ⅱ)營(yíng)口臺(tái)(YKO)和蓋縣臺(tái)(GAX)快剪切波優(yōu)勢(shì)偏振方向的變化(圖5)。營(yíng)口臺(tái)(YKO)在時(shí)段I(圖5a)、 時(shí)段Ⅱ(5b)快波的優(yōu)勢(shì)偏振方向均為NEE-SWW向,與此前營(yíng)口臺(tái)(YKO)全時(shí)段數(shù)據(jù)結(jié)果(圖3a)完全一致。由此或可推斷,海城老震區(qū)的區(qū)域應(yīng)力場(chǎng)較為均勻,或處于應(yīng)力均勻加載的過(guò)程中,該區(qū)小地震活動(dòng)或?yàn)楸尘靶缘卣鸹顒?dòng)。蓋縣臺(tái)(GAX)在時(shí)段I(圖5c)、 時(shí)段Ⅱ(5d)快剪切波的優(yōu)勢(shì)偏振方向差異較大,時(shí)段I蓋縣臺(tái)(GAX)的快剪切波優(yōu)勢(shì)偏振方向?yàn)镹NE向,與臺(tái)站下方的金州斷裂的走向大體相近。而時(shí)段Ⅱ蓋縣臺(tái)(GAX)快剪切波的優(yōu)勢(shì)偏振方向?yàn)镹WW-SEE,與營(yíng)口臺(tái)(YKO)的快剪切波優(yōu)勢(shì)方向基本接近,同時(shí)與該區(qū)最大主應(yīng)力方位較為一致,蓋縣臺(tái)(GAX)快剪切波的優(yōu)勢(shì)偏振方向在時(shí)段Ⅰ和時(shí)段Ⅱ變化近90°。關(guān)于臺(tái)站極性發(fā)生近90°偏轉(zhuǎn)的現(xiàn)象此前也有相關(guān)報(bào)道,Crampin (2002)根據(jù)冰島的小震資料進(jìn)行剪切波分裂研究時(shí)發(fā)現(xiàn),大部分臺(tái)站的快剪切波偏振方向?yàn)镹E-SW方向,與區(qū)域應(yīng)力場(chǎng)方向一致,但位于冰島北部1條活動(dòng)斷層附近的3個(gè)臺(tái)站的快剪切波偏振方向?yàn)镹W-SE,垂直于其他地區(qū)的快剪切波偏振方向,認(rèn)為這種現(xiàn)象可能是地震射線(xiàn)在裂隙中傳播時(shí)受到高孔隙流體壓力的影響而引起的,稱(chēng)為90°翻轉(zhuǎn)現(xiàn)象。由此推斷,蓋州震群自2013年12月22日的這組增強(qiáng)活動(dòng)或與區(qū)域應(yīng)力的局部增強(qiáng)(調(diào)整)所導(dǎo)致的高孔隙壓變化有關(guān)。

    此外,我們研究了營(yíng)口臺(tái)(YKO)和蓋縣臺(tái)(GAX)的慢波時(shí)間延遲(Δt)的變化。營(yíng)口臺(tái)(YKO)的Δt分布在0.52~7.45ms/km范圍內(nèi),均值為3.21ms/km,其在時(shí)段I內(nèi)的均值為3.23ms/km,在時(shí)段Ⅱ內(nèi)的均值為3.12ms/km; 蓋縣臺(tái)(GAX)的Δt分布在0.76~5.06ms/km范圍內(nèi),均值為2.36ms/km,其在時(shí)段I內(nèi)的均值為2.46ms/km,在時(shí)段Ⅱ內(nèi)的均值為2.22ms/km。由此可以看出無(wú)論是營(yíng)口臺(tái)(YKO)還是蓋縣臺(tái)(GAX),其剪切波分裂參數(shù)時(shí)間延遲并未出現(xiàn)明顯的波動(dòng)變化。而根據(jù)震前應(yīng)力釋放的結(jié)論(Gaoetal.,2008)在強(qiáng)震前慢波時(shí)間延遲(Δt)值會(huì)有個(gè)下降的過(guò)程。以往的相關(guān)研究結(jié)果也表明,強(qiáng)震前幾十分鐘至幾個(gè)月應(yīng)力開(kāi)始釋放,微裂紋的合并與聚合形成觸發(fā)巖石破壞水平的破裂(Gaoetal.,2004)。如2014年5月30日盈江6.1級(jí)地震的前震序列剪切波分裂慢波時(shí)間延遲總體上呈現(xiàn)較低的水平,而6.1級(jí)地震后的慢波時(shí)間延遲值相對(duì)較高(李金等,2015)。此外,在2001年6月8日云南施甸5.3級(jí)強(qiáng)余震前(高原等,2004b)、 2000年6月17日冰島6.6級(jí)地震前(Volitietal.,2003)也都出現(xiàn)了類(lèi)似的現(xiàn)象。由蓋縣臺(tái)(GAX)所得到的剪切波分裂參數(shù)慢波時(shí)間延遲或可推斷,蓋州地區(qū)的地震集中活動(dòng),其表現(xiàn)形式為震群形式,盡管其小震頻度相對(duì)較高,但從目前的結(jié)果尚不能看出該震群有演變成前震序列的可能性。

    圖5 2013年12月22日前、 后營(yíng)口臺(tái)(YKO)(a,b)、 蓋縣臺(tái)(GAX)(c,d)的快剪切波偏振方向等面積極射投影與等面積投影玫瑰圖Fig. 5 The homolographic projection and homolographic projection rose graphs of the polarizations of the fast shear-wave splittings at stations YKO and GAX before December 22,2013(a,b) and after December 22,2013(c,d).

    4 結(jié)論與討論

    本研究采用SAM分析方法,針對(duì)2012年2月12日發(fā)生的蓋州震群,研究其剪切波分裂特征,同時(shí)對(duì)比分析距其不遠(yuǎn)的海城老震區(qū)小震活動(dòng)的剪切波分裂特征。得到了蓋州—海城地區(qū)地殼各向異性的初步結(jié)果。

    研究結(jié)果表明,營(yíng)口臺(tái)(YKO)反映的海城老震區(qū)小震活動(dòng)的快剪切波偏振的優(yōu)勢(shì)方向?yàn)镹EE向,該優(yōu)勢(shì)偏振方向與太齡雪(2008)利用1999年11月29日遼寧岫巖MS5.9地震前后營(yíng)口臺(tái)(YKO)資料得到的優(yōu)勢(shì)偏振方向完全一致。與該地區(qū)的最大主應(yīng)力的方位基本一致。即快剪切波的偏振方向可以反映該區(qū)的構(gòu)造應(yīng)力方向。由海城老震區(qū)小震活動(dòng)的空間分布近EW向的優(yōu)勢(shì)分布可以推測(cè)營(yíng)口臺(tái)(YKO)下方可能存在EW向的活動(dòng)斷裂。

    蓋縣臺(tái)(GAX)所反映的蓋州震群的快波偏振優(yōu)勢(shì)方向?yàn)镹NE向和NWW向,NWW向的優(yōu)勢(shì)方向與當(dāng)?shù)刈畲笾鲬?yīng)力方向較為一致。根據(jù)蓋州震群精定位后的震中分布所表現(xiàn)出的EW向的優(yōu)勢(shì)分布和快剪切波NWW向的優(yōu)勢(shì)方向,推斷該臺(tái)站下方可能存在1條近EW向的隱伏斷裂。最新對(duì)該震群的精定位及沿不同剖面的分析結(jié)果亦支持該震群活動(dòng)發(fā)生在金州斷裂北段與海城-營(yíng)口斷裂南段的1個(gè)新的活動(dòng)斷層上,認(rèn)為該斷層破裂尺度較小,長(zhǎng)度不超過(guò)6km,深度較淺,大約在3~8km左右(王亮等,2015)。但該斷層是否真的存在,仍需更多野外地質(zhì)資料進(jìn)行證實(shí)。

    針對(duì)蓋州震群2013年底以來(lái)的增強(qiáng)活動(dòng),研究了2013年12月22日前(時(shí)段I)、 后(時(shí)段Ⅱ)營(yíng)口臺(tái)(YKO)和蓋縣臺(tái)(GAX)快剪切波優(yōu)勢(shì)偏振方向的變化。結(jié)果顯示,營(yíng)口臺(tái)(YKO)快剪切波優(yōu)勢(shì)偏振方向在該時(shí)間節(jié)點(diǎn)前、 后沒(méi)有變化。由此或可推斷,海城老震區(qū)區(qū)域應(yīng)力場(chǎng)較為均勻,或處于應(yīng)力均勻加載的過(guò)程中,該區(qū)小地震活動(dòng)或?yàn)楸尘靶缘卣鸹顒?dòng)。蓋縣臺(tái)(GAX)在時(shí)段I(圖5c)、 時(shí)段Ⅱ(圖5d)快剪切波的優(yōu)勢(shì)偏振方向差異較大,時(shí)段I蓋縣臺(tái)(GAX)的快剪切波優(yōu)勢(shì)偏振方向?yàn)镹NE向,與臺(tái)站下方的金州斷裂的走向較為接近。時(shí)段Ⅱ蓋縣臺(tái)(GAX)快剪切波的優(yōu)勢(shì)偏振方向?yàn)镹WW,與營(yíng)口臺(tái)(YKO)的快剪切波優(yōu)勢(shì)方向基本接近,同時(shí)與該區(qū)最大主應(yīng)力方位較為一致。由此可以推斷,蓋州震群自2013年12月22日以來(lái)的這組增強(qiáng)活動(dòng)或與區(qū)域應(yīng)力的局部增強(qiáng)有關(guān)。

    此外,針對(duì)營(yíng)口臺(tái)(YKO)和蓋縣臺(tái)(GAX)的慢波時(shí)間延遲的研究結(jié)果顯示,這2個(gè)震群的慢波時(shí)間延遲就平均值而言在2013年12月22日前后變化并不明顯,目前的初步結(jié)果尚不能確定其有演變成前震序列的可能,對(duì)剪切波分裂參數(shù)時(shí)間變化的細(xì)節(jié)探討尚待今后的進(jìn)一步分析。

    致謝 感謝蔣海昆研究員多年悉心的指導(dǎo)。遼寧省地震局監(jiān)測(cè)中心為本研究提供了波形數(shù)據(jù); 本項(xiàng)工作得到了中國(guó)地震局 “地震分析預(yù)測(cè)研究青年工作組”的支持; 研究過(guò)程中,中國(guó)地震局地震預(yù)測(cè)研究所石玉濤副研究員、 中國(guó)地震局第二監(jiān)測(cè)中心劉庚助理工程師、 中國(guó)地震局地球物理研究所肖卓博士、 安徽省地震局鮑子文助理工程師提供了諸多幫助: 在此一并表示衷心的感謝。文中部分圖件采用GMT繪制。

    曹鳳娟,王亮,焦明若,等. 2013. 2012年2月2日蓋州ML4.7震群序列特征及發(fā)震構(gòu)造探討 [J]. 地震地質(zhì),35(4): 842—852. doi: 10.3969/j.issn.0253-4967.2013.04.03.

    CAO Feng-juan,WANG Liang,JIAO Ming-ruo,etal. 2013. Discussion on the characteristics and seismo-tectonics of the GaizhouML4.7 earthquake swarm sequence on February 2,2012 [J]. Seismology and Geology,35(4): 842—852(in Chinese).

    曹鳳娟,王連權(quán),張萍,等. 2005. 岫巖MS5.4地震前后S波分裂研究 [J]. 東北地震研究,21(4): 39— 44.

    CAO Feng-juan,WANG Lian-quan,ZHANG Ping,etal. 2005. Study on S-wave fission around XiuyanMS5.4 earthquake [J]. Seismological Research of Northeast China,21(4): 39— 44(in Chinese).

    高常波,鐘以章. 1998. 東北書(shū)友管道場(chǎng)地?cái)鄬踊顒?dòng)性與地震危險(xiǎn)性研究 [M]. 北京: 地震出版社. 103—159.

    GAO Chang-bo,ZHONG Yi-zhang. 1998. Fault Activity and Seismic Risk Research for Petroleum Transmission Pipeline in NE China [M]. Seismological Press,Beijing. 103—159(in Chinese).

    高原,梁維,丁香,等. 2004b. 云南2001年施甸地震的剪切波分裂參數(shù)變化特征 [J]. 地震學(xué)報(bào),26(6): 576—582.

    GAO Yuan,LIANG Wei,DING Xiang,etal. 2004b. Variational characteristics of shear-wave splitting on the 2001 Shidian earthquakes in Yunnan,China[J]. Acta Seismologica Sinica,26(6): 576—582(in Chinese).

    高原,劉希強(qiáng),梁維. 2004a. 剪切波分裂系統(tǒng)分析方法(SAM)軟件系統(tǒng) [J]. 中國(guó)地震,20(1): 101—107.

    GAO Yuan,LIU Xi-qiang,LIANG Wei. 2004a. Systematic analysis method of shear-wave splitting: SAM software system [J]. Earthquake Research in China,20(1): 101—107(in Chinese).

    高原,吳晶,易桂喜,等. 2010. 從殼幔地震各向異性初探華北地區(qū)殼幔耦合關(guān)系 [J]. 科學(xué)通報(bào),55(29): 2837—2843.

    GAO Yuan,WU Jing,Yi Gui-xi,etal. 2010. Crust-mantle coupling in North China: Preliminary analysis from seismic anisotropy [J]. Chinese Sci Bull,55(29): 2837—2843(in Chinese).

    高原,鄭斯華,孫勇. 1995. 唐山地區(qū)地殼裂隙各向異性 [J]. 地震學(xué)報(bào),17(3): 283—293.

    GAO Yuan,ZHENG Si-hua,SUN Yong. 1995. Crustal crack anisotropy in Tangshan [J]. Acta Seismologica Sinica,17(3): 283—293(in Chinese).

    高原,鄭斯華,周慧蘭. 1999. 唐山地區(qū)快剪切波偏振圖像及其變化 [J]. 地球物理學(xué)報(bào),42(2): 228—232.

    GAO Yuan,ZHENG Si-hua,ZHOU Hui-lan. 1999. Polarization patterns of fast shear wave in Tangshan region and their variations [J]. Chinese Journal of Geophysics,42(2): 228—232(in Chinese).

    葛榮峰,張慶龍,解國(guó)愛(ài),等. 2009. 郯廬斷裂帶北段及鄰區(qū)現(xiàn)代地震活動(dòng)性與應(yīng)力狀態(tài) [J]. 地震地質(zhì),31(1): 141—154. doi: 10.3969/j.issn.0253-4967.2009.01.013.

    GE Rong-feng,ZHANG Qing-long,XIE Guo-ai,etal. 2009. Seismic activity and stress state of the northern Tanlu fault zone and its adjacent areas [J]. Seismology and Geology,31(1): 141—154(in Chinese).

    賴(lài)院根,劉啟元,陳九輝,等. 2006. 首都圈地區(qū)橫波分裂與地殼應(yīng)力場(chǎng)特征 [J]. 地球物理學(xué)報(bào),49(1): 189—196.

    LAI Yuan-gen,LIU QI-Yuan,CHEN Jiu-hui,etal. 2006. Shear wave splitting and the features of crustal stress field in the Capital circle [J]. Chinese Journal of Geophysics,49(1): 189—196(in Chinese).

    雷清清,廖旭,董曉燕. 2008. 遼寧省地震構(gòu)造研究 [J]. 東北地震研究,24(4): 1—10.

    LEI Qing-qing,LIAO Xu,DONG Xiao-yan. 2008. Seismotectonic study on Liaoning Province [J]. Seismological Research of Northeast China,24(4): 1—10(in Chinese).

    李金,高原,徐甫坤,等. 2015. 2014年5月30日盈江6.1級(jí)地震序列剪切波分裂研究 [J]. 中國(guó)地震,31(2): 245—252.

    LI Jin,GAO Yuan,XU Fu-kun,etal. 2015. The shear-wave splitting study on the YingjiangMS6.1 earthquake sequence on May 30,2014 [J]. Earthquake Research in China,31(2): 245—252(in Chinese).

    呂培苓,楊智嫻,王曉青,等. 2004. 1999年岫巖地震序列研究 [J]. 地震,24(4): 45—50.

    Lü Pei-ling,YANG Zhi-xian,WANG Xiao-qing,etal. 2004. Investigation of the 1999 Xiuyan earthquake sequence,Liaoning Province,China [J]. Earthquake,24(4): 45—50(in Chinese).

    盧造勛,蔣秀琴,潘科,等. 2002. 中朝地臺(tái)東北緣地區(qū)的地震層析成像 [J]. 地球物理學(xué)報(bào),45(3): 338—351.

    LU Zao-xun,JIANG Xiu-qin,PAN Ke,etal. 2002. Seismic tomography in the northeast margin area of Sino-Korean platform [J]. Chinese Journal of Geophysics,45(3): 338—351(in Chinese).

    石玉濤,高原,吳晶,等. 2006. 云南地區(qū)地殼介質(zhì)各向異性: 快剪切波偏振特性 [J]. 地震學(xué)報(bào),28(6): 574—585.

    SHI Yu-tao,GAO Yuan,WU Jing,etal. 2006. Seismic anisotropy of the crust in Yunnan,China: Polarizations of fast shear-waves [J]. Acta Seismologica Sinica,28(6): 574—585(in Chinese).

    太齡雪,高原,曹鳳娟,等. 2008. 遼寧1999年MS5.9岫巖地震的剪切波分裂特征 [J]. 地震學(xué)報(bào),30(4): 340—354.

    TAI Ling-xue,GAO Yuan,CAO Feng-juan,etal. 2008. Shear-wave splitting before and after the 1999 Xiuyan earthquake in Liaoning,China [J]. Acta Seismologica Sinica,30(4): 340—354(in Chinese).

    太齡雪,高原,石玉濤,等. 2009. 遼寧區(qū)域地震臺(tái)網(wǎng)的地殼剪切波分裂研究 [J]. 地震地質(zhì),31(3): 401— 414.

    TAI Ling-xue,GAO Yuan,SHI Yu-tao,etal. 2009. A study of shear-wave splitting in the crust by Liaoning telemetry seismic network of China [J]. Seismology and Geology,31(3): 401— 414(in Chinese).

    王亮,黃金水,焦明若,等. 2015. 蓋州震群的重新定位 [J]. 中國(guó)地震,31(2): 439— 446.

    WANG Liang,HUANG Jin-shui,JIAO Ming-ruo,etal. 2015. Accurate relocation of earthquake swarm in Gaizhou,Liaoning,China [J]. Earthquake Research in China,31(2): 439— 446(in Chinese).

    吳晶,高原,蔡晉安,等. 2007. 華夏地塊東南部地殼地震各向異性特征初步研究 [J]. 地球物理學(xué)報(bào),60(6): 1748—1756.

    WU Jing,GAO Yuan,CAI Jin ̄an,etal. 2007. Preliminary study on seismic anisotropy in the crust of Cathaysia Block [J].Chinese Journal of Geophysics,60(6): 1748—1756(in Chinese).

    徐杰,宋長(zhǎng)青,高戰(zhàn)武,等. 1999. 營(yíng)口-濰坊斷裂帶新生代活動(dòng)特征 [J]. 地震地質(zhì),21(4): 289—300.

    XU Jie,SONG Chang-qing,GAO Zhan-wu,etal. 1999. Active characters of the Yingkou-Weifang fault zone in the Cenozoic era [J]. Seismology and Geology,21(4): 289—300(in Chinese).

    徐杰,周本剛,計(jì)鳳桔,等. 2012. 華北渤海灣盆地區(qū)大震發(fā)震構(gòu)造的基本特征 [J]. 地震地質(zhì),34(4): 618— 636. doi: 10.3969/j.issn.0253-4967.2012.04.008.

    XU Jie,ZHOU Ben-gang,JI Feng-ju,etal. 2012. Features of seismogenic structures of great earthquakes in the Bohai Bay Basin area,North China [J]. Seismology and Geology,34(4): 618— 636(in Chinese).

    張萍,谷光裕,高艷玲. 2001b. 岫巖-海城MS5.4地震序列震源機(jī)制解[J]. 地震,21(1): 98—101.

    ZHANG Ping,GU Guang-yu,GAO Yan-ling. 2001b. Focal mechanism of Xiuyan-Haicheng(MS5.4)earthquake sequence [J]. Earthquake,21(1): 98—101(in Chinese).

    張萍,蔣秀琴. 2001a. 岫巖-海城MS5.4地震序列的震源機(jī)制解及應(yīng)力場(chǎng)特征[J]. 地震地磁觀測(cè)與研究,22(2): 76—82.

    ZHANG Ping,JIANG Xiu-qin. 2001a. The focal mechanism solutions and the crust stress field characteristics in Xiuyan-Haicheng(MS5.4)earthquake sequence [J]. Seismological and Geomagnetic Observation and Research,22(2): 76—82(in Chinese).

    鄭秀芬,陳朝輝,張春賀. 2008. 1999年臺(tái)灣集集地震余震區(qū)-嘉義地區(qū)地震的剪切波分裂參數(shù)隨時(shí)間的變化的研究[J]. 地球物理學(xué)報(bào),51(1): 149—157.

    ZHENG Xiu-fen,CHEN Chao-hui,ZHANG Chun-he. 2008. Study on temporal variations of shear wave splitting in the Chiayi area,aftershock zone of 1999 Chichi earthquake,Taiwan [J]. Chinese J Geophys,51(1): 149—157(in Chinese).

    鐘以章,白云. 1999. 遼東半島及兩側(cè)海域深部構(gòu)造、 活動(dòng)構(gòu)造與地震關(guān)系研究 [J]. 華北地震科學(xué),17(2): 59— 65.

    ZHONG YI-zhang,BAI Yun. 1999. Study on deep tectonics,active tectonics and their relation with earthquakes in Liaodong Peninsula and its surrounding sea area [J]. North China Earthquake Science,17(2): 59— 65(in Chinses).

    Crampin S. 1978. Seismic-wave propagation through a cracked solid: Polarization as a possible dilatancy diagnostic [J]. Geophys J R Astr Soc,53: 467— 496.

    Crampin S. 1981. A review of wave motion in anisotropic and cracked elastic-media [J]. Wave Motion,3: 343—391.

    Crampin S,Chastin S. 2003. A review of shear-wave splitting in the crack-critical crust [J]. Geophys J Int,155: 221—240.

    Crampin S,Peacock S. 2005. A review of shear-wave splitting in the compliant crack-critical anisotropic earth [J]. Wave Motion,41: 59—77.

    Crampin S,Volti T,Chastin,etal. 2002. Indication of high pore-fluid pressures in a seismically-active fault zone [J]. Geophys J Int,151: F1—F5.

    Gao Y,Crampin S. 2004. Observation of stress relaxation before earthquake [J]. Geophys J Int,157(2): 578—582.

    Gao Y,Crampin S. 2008. Shear-wave splitting and earthquake forecasting [J]. Terra Nova,20(6): 440— 448.

    Gao Y,Wang P D,Zheng S H,etal. 1998. Temporal changes in shear wave splitting at an isolated swarm of small earthquakes in 1992 near Dongfang,Hainan Island,southern China [J]. Geophys J Int,135(1): 102—112.

    Gao Y,Wu J,F(xiàn)ukao Y,etal. 2011. Shear-wave splitting in the crust in North China: Stress,faults and tectonic implications [J]. Geophys J Int,187(2): 642— 654.

    Voliti T,Crampin S. 2003. A four-year study of shear-wave splitting in Iceland: 1. Background and preliminary anylysis [M]∥Nieuwland D A(ed). New Insights into Structural Interpretation and Modelling. Geol Soc Lond Spec Publ,212: 117—133.

    Zatsepin S V,Crampin S. 1997. Modeling the compliance of crustal rock: I-Response of shear-wave splitting to differential stress [J]. Geophys J Int,129: 477— 494.

    THE SHEAR-WAVE SPLITTING OF GAIZHOU EARTHQUAKE SWARM IN LIAONING

    LI Jin1)GAO Yuan2)ZHANG Bo3)WANG Liang3)GAO Rong1)

    1)EarthquakeAdministrationofXinjiangUygurAutonomousRegion,Urumqi830011,China2)InstituteofEarthquakeScience,ChinaEarthquakeAdiministration,Beijing100036,China3)EarthquakeAdministrationofLiaoningProvince,Shenyang110034,China

    AnML4.7 earthquake occurred on February 2,2012 in Liaoning Gaizhou(40.56°N,122.36°E),since then,small earthquakes are frequent in this area,and until now the seismic activity does not stop,several earthquakes with magnitude larger than 4.0 have occurred. As of October 30,2014,1223 earthquakes have happened in the Gaizhou area,including 934 earthquakes with the magnitudeML1.0~1.9,247 with the magnitudeML2.0~2.9 and 45 with the magnitudeML3.0~3.9. Meanwhile,earthquakes are continuously active in Haicheng area where theMS7.3 earthquake happened in 1975,and there are over 1100 earthquakes(ML≥1.0)having occurred since the Gaizhou earthquake swarm activity. Because the polarization direction of the fast shear wave is very sensitive to the variation of the principal stress environment,the shear wave splitting parameter can reflect the regional stress state and the local structural features,especially effective for the analysis of small-scale stress environment characteristics. So based on the seismic activities of the two earthquake clusters,this study analyzes the characteristics of shear-wave splitting in Gaizhou-Haicheng area. Preliminary results show that predominant polarization direction of fast shear-waves in the old earthquake region of Haicheng is stable,consistent with the direction of regional stress field. Due to presence of active fault below the Gaixian station(GAX),the predominant polarization direction of fast shear-waves is more complicated. There are two predominant polarizations,consistent respectively with Jinzhou Fault strike which is below the station and the maximum principal stress direction in this area. In addition,Gaizhou earthquake swarm activity increased after December 22,2013,and after the time node,the predominant polarization direction of fast shear-waves in Gaixian station is SEE,which is close to the predominant polarization direction of fast shear-waves in Yingkou station,at the same time consistent with the maximum principal stress direction of this region. Thus it can be inferred,the enhanced activity of Gaizhou earthquake swarm since December 22,2014 may be related to local enhancement of regional stress. In addition,the average time-delays of slow waves in station YKO and GAX show that there are no obvious changes before and after the time point of December 22,2013,which is different greatly with the previous related researches on the variation of slow wave time-delays,and there is no possibility that the Gaizhou earthquake swarm evolved into foreshock sequences from current preliminary results. We should do more work to study the details of the time delay variation of shear wave splitting parameter.

    shear-wave splitting,Gaizhou earthquake swarm,polarization of fast shear-waves

    10.3969/j.issn.0253- 4967.2016.04.020

    2015-09-24收稿,2016-05-04改回。

    地震科技星火計(jì)劃(XH17041Y)、 2015年度震情跟蹤合同制定向工作任務(wù)(2015010113)與新疆維吾爾自治區(qū)青年自然科學(xué)基金(2014211B39)共同資助。

    P315.6

    A

    0253-4967(2016)04-1058-12

    李金,男,1986年出生,2012年畢業(yè)于中國(guó)地震局地震預(yù)測(cè)研究所固體地球物理專(zhuān)業(yè),獲理學(xué)碩士學(xué)位,助理研究員,主要從事數(shù)字地震學(xué)及地震預(yù)測(cè)研究工作,電話(huà): 15299475246,E-mail: lijin6205@163.com。

    猜你喜歡
    海城營(yíng)口偏振
    嫁給愛(ài)情
    偏振糾纏雙光子態(tài)的糾纏特性分析
    電子制作(2019年12期)2019-07-16 08:45:20
    基于LabVIEW的偏振調(diào)制激光測(cè)距測(cè)量控制系統(tǒng)
    SUPPLEMENTAL DESCRIPTION OF MYXOBOLUS HAICHENGENSIS CHEN, 1958(MYXOZOA: MYXOSPOREA) INFECTING THE GILLS OF ABBOTTINA RIVULARIS BASILEWAKY: MORPHOLOGICAL ANDMOLECULAR DATA
    海城河非均勻流水面線(xiàn)的推算研究
    偏振旋轉(zhuǎn)效應(yīng)在全光緩存器中的應(yīng)用與實(shí)現(xiàn)
    《營(yíng)口古代文明展》的內(nèi)容設(shè)計(jì)與特色
    海城的容顏和氣質(zhì)
    今日遼寧(2015年4期)2015-04-18 07:52:11
    高分-2衛(wèi)星營(yíng)口融合影像
    營(yíng)口機(jī)場(chǎng)建設(shè)對(duì)生態(tài)環(huán)境的影響研究
    精品乱码久久久久久99久播| 国产精品爽爽va在线观看网站 | 成在线人永久免费视频| 成人av一区二区三区在线看| 一本大道久久a久久精品| 亚洲欧洲精品一区二区精品久久久| 非洲黑人性xxxx精品又粗又长| 欧美日本亚洲视频在线播放| 日韩 欧美 亚洲 中文字幕| 精品一品国产午夜福利视频| 亚洲av成人av| 亚洲七黄色美女视频| 在线观看午夜福利视频| 国产欧美日韩综合在线一区二区| 亚洲成国产人片在线观看| 老司机福利观看| 国内久久婷婷六月综合欲色啪| 人妻丰满熟妇av一区二区三区| 中文字幕久久专区| 一级片免费观看大全| 久久热在线av| 中文字幕人妻熟女乱码| 无限看片的www在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 色综合婷婷激情| 亚洲国产精品久久男人天堂| 日韩三级视频一区二区三区| 亚洲色图av天堂| 国产激情欧美一区二区| 午夜精品久久久久久毛片777| 男女下面进入的视频免费午夜 | 一区二区三区精品91| 亚洲国产看品久久| 国产伦一二天堂av在线观看| 久久影院123| 亚洲人成电影观看| 国产真人三级小视频在线观看| 丁香欧美五月| 啦啦啦观看免费观看视频高清 | 国产成+人综合+亚洲专区| 亚洲人成电影观看| 久久婷婷人人爽人人干人人爱 | 国产真人三级小视频在线观看| 12—13女人毛片做爰片一| 日日摸夜夜添夜夜添小说| 麻豆av在线久日| 欧美一级毛片孕妇| 天天添夜夜摸| 亚洲精品一区av在线观看| 久久精品国产99精品国产亚洲性色 | 国产成人精品久久二区二区91| 亚洲精品中文字幕在线视频| 夜夜夜夜夜久久久久| 黄网站色视频无遮挡免费观看| 精品久久久久久久久久免费视频| 嫩草影院精品99| 侵犯人妻中文字幕一二三四区| 黄片大片在线免费观看| 色在线成人网| 看黄色毛片网站| 亚洲精品国产一区二区精华液| 男女之事视频高清在线观看| 每晚都被弄得嗷嗷叫到高潮| 久久伊人香网站| 成在线人永久免费视频| 日韩精品中文字幕看吧| 两个人看的免费小视频| 中国美女看黄片| 国产亚洲精品第一综合不卡| 日日摸夜夜添夜夜添小说| 亚洲成国产人片在线观看| videosex国产| 久久精品亚洲熟妇少妇任你| 人人澡人人妻人| 18禁裸乳无遮挡免费网站照片 | 亚洲性夜色夜夜综合| 99国产精品99久久久久| 国产区一区二久久| 精品欧美国产一区二区三| 天堂√8在线中文| 日韩欧美国产一区二区入口| 日日摸夜夜添夜夜添小说| 欧美黄色片欧美黄色片| 亚洲精品av麻豆狂野| 精品一区二区三区四区五区乱码| 丝袜人妻中文字幕| 精品人妻在线不人妻| 9191精品国产免费久久| 国产成人av激情在线播放| 两个人看的免费小视频| 国产99久久九九免费精品| 国产1区2区3区精品| 女人爽到高潮嗷嗷叫在线视频| 性少妇av在线| 成年人黄色毛片网站| 黄色视频,在线免费观看| 一区二区日韩欧美中文字幕| 大型黄色视频在线免费观看| 不卡av一区二区三区| 欧美一级a爱片免费观看看 | 动漫黄色视频在线观看| 咕卡用的链子| 最近最新中文字幕大全电影3 | 国产欧美日韩精品亚洲av| 久久久国产欧美日韩av| 日本黄色视频三级网站网址| 搞女人的毛片| 夜夜躁狠狠躁天天躁| 自拍欧美九色日韩亚洲蝌蚪91| 女人被狂操c到高潮| 精品熟女少妇八av免费久了| av视频在线观看入口| videosex国产| 国产又爽黄色视频| 精品人妻1区二区| 少妇熟女aⅴ在线视频| 成人欧美大片| 黄色丝袜av网址大全| 日本在线视频免费播放| 精品日产1卡2卡| 可以免费在线观看a视频的电影网站| 国内精品久久久久久久电影| 一边摸一边抽搐一进一小说| 欧美成人一区二区免费高清观看 | 亚洲,欧美精品.| 天堂√8在线中文| 天堂√8在线中文| 少妇被粗大的猛进出69影院| 国产欧美日韩一区二区三| 日韩 欧美 亚洲 中文字幕| 久久国产精品男人的天堂亚洲| 成人av一区二区三区在线看| av免费在线观看网站| 757午夜福利合集在线观看| 国产精品一区二区精品视频观看| 亚洲欧美日韩另类电影网站| 亚洲三区欧美一区| 男人舔女人的私密视频| 亚洲国产看品久久| 免费少妇av软件| 日韩免费av在线播放| 一级a爱视频在线免费观看| 超碰成人久久| 久久国产亚洲av麻豆专区| 自拍欧美九色日韩亚洲蝌蚪91| 日本vs欧美在线观看视频| 亚洲人成电影免费在线| 欧美日韩乱码在线| 日本a在线网址| 午夜亚洲福利在线播放| 51午夜福利影视在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 天天一区二区日本电影三级 | 久久狼人影院| 中文字幕av电影在线播放| 亚洲成a人片在线一区二区| 一个人观看的视频www高清免费观看 | 国产成人啪精品午夜网站| 国产极品粉嫩免费观看在线| 老司机福利观看| 成年人黄色毛片网站| 99精品久久久久人妻精品| 最近最新中文字幕大全免费视频| 19禁男女啪啪无遮挡网站| 日韩三级视频一区二区三区| av视频免费观看在线观看| 亚洲,欧美精品.| 好男人电影高清在线观看| 女同久久另类99精品国产91| 亚洲国产日韩欧美精品在线观看 | 国产91精品成人一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 亚洲国产精品合色在线| 人人澡人人妻人| x7x7x7水蜜桃| 国产人伦9x9x在线观看| 欧美日韩福利视频一区二区| 精品免费久久久久久久清纯| 一边摸一边抽搐一进一小说| 琪琪午夜伦伦电影理论片6080| 一二三四社区在线视频社区8| 99re在线观看精品视频| 亚洲精品中文字幕在线视频| 精品国产亚洲在线| 咕卡用的链子| 一级毛片女人18水好多| 国产成+人综合+亚洲专区| 高清在线国产一区| 成年女人毛片免费观看观看9| 热99re8久久精品国产| 嫁个100分男人电影在线观看| 中文字幕另类日韩欧美亚洲嫩草| 亚洲av日韩精品久久久久久密| 久久久国产成人免费| 亚洲成a人片在线一区二区| 女生性感内裤真人,穿戴方法视频| 女生性感内裤真人,穿戴方法视频| 99久久国产精品久久久| 中文字幕久久专区| 老司机靠b影院| 露出奶头的视频| 亚洲一区高清亚洲精品| 又黄又粗又硬又大视频| 午夜影院日韩av| 一个人免费在线观看的高清视频| 国产精品美女特级片免费视频播放器 | 精品久久久久久成人av| 9热在线视频观看99| 少妇粗大呻吟视频| 久久久久久大精品| 午夜老司机福利片| 欧洲精品卡2卡3卡4卡5卡区| 色播在线永久视频| 妹子高潮喷水视频| av欧美777| 久久人妻福利社区极品人妻图片| 香蕉国产在线看| 丰满人妻熟妇乱又伦精品不卡| 亚洲欧美日韩另类电影网站| 国产精品香港三级国产av潘金莲| 国产不卡一卡二| aaaaa片日本免费| 国产麻豆成人av免费视频| 亚洲黑人精品在线| av在线天堂中文字幕| 91成年电影在线观看| 变态另类丝袜制服| 高清黄色对白视频在线免费看| 欧美激情高清一区二区三区| 最近最新免费中文字幕在线| 亚洲性夜色夜夜综合| 久久性视频一级片| 两性夫妻黄色片| 中文字幕最新亚洲高清| 亚洲色图综合在线观看| 啦啦啦观看免费观看视频高清 | 国产精品影院久久| 国产成人欧美| 国产真人三级小视频在线观看| 国产男靠女视频免费网站| 免费在线观看影片大全网站| 女同久久另类99精品国产91| 老熟妇乱子伦视频在线观看| 精品国产超薄肉色丝袜足j| 国产激情欧美一区二区| 午夜激情av网站| 久久精品亚洲精品国产色婷小说| 美国免费a级毛片| 好看av亚洲va欧美ⅴa在| 一区二区三区国产精品乱码| 女性被躁到高潮视频| netflix在线观看网站| 亚洲色图av天堂| 91九色精品人成在线观看| 国产99久久九九免费精品| 搞女人的毛片| 窝窝影院91人妻| 久久国产精品影院| 成人国产综合亚洲| 久热爱精品视频在线9| 国产精品免费视频内射| 亚洲精品一卡2卡三卡4卡5卡| 日韩欧美国产一区二区入口| 黑人巨大精品欧美一区二区蜜桃| 国产精品久久久人人做人人爽| 欧美在线黄色| 美女免费视频网站| 日韩精品免费视频一区二区三区| 午夜福利欧美成人| 国产亚洲欧美在线一区二区| 88av欧美| 亚洲精华国产精华精| 一级,二级,三级黄色视频| 最近最新免费中文字幕在线| 精品一区二区三区av网在线观看| 国产av一区在线观看免费| 91大片在线观看| 精品国内亚洲2022精品成人| 亚洲一区中文字幕在线| 欧美精品亚洲一区二区| 啦啦啦免费观看视频1| 日韩精品免费视频一区二区三区| 大型av网站在线播放| 电影成人av| 老司机午夜福利在线观看视频| 国产亚洲欧美98| 制服诱惑二区| 亚洲精品中文字幕一二三四区| 波多野结衣一区麻豆| 久久香蕉激情| 高潮久久久久久久久久久不卡| 高潮久久久久久久久久久不卡| 亚洲av片天天在线观看| 色播在线永久视频| 女人爽到高潮嗷嗷叫在线视频| 精品不卡国产一区二区三区| 久热爱精品视频在线9| 亚洲精品久久成人aⅴ小说| 日韩欧美一区二区三区在线观看| 久久精品成人免费网站| 国产亚洲av嫩草精品影院| 欧美国产日韩亚洲一区| 黄片播放在线免费| 国产99白浆流出| 不卡一级毛片| 日韩欧美一区二区三区在线观看| 18禁国产床啪视频网站| 国产精华一区二区三区| 国产精品美女特级片免费视频播放器 | 老熟妇仑乱视频hdxx| 久久久久久国产a免费观看| 久久久精品欧美日韩精品| 好男人在线观看高清免费视频 | 99国产精品一区二区蜜桃av| 日韩高清综合在线| 好男人在线观看高清免费视频 | 老司机午夜福利在线观看视频| 在线十欧美十亚洲十日本专区| 无遮挡黄片免费观看| 成人免费观看视频高清| 亚洲欧美激情综合另类| 波多野结衣巨乳人妻| 亚洲精品中文字幕在线视频| 他把我摸到了高潮在线观看| 大型av网站在线播放| 久久久国产成人精品二区| 日日爽夜夜爽网站| 麻豆成人av在线观看| 1024香蕉在线观看| 男人的好看免费观看在线视频 | 国产高清videossex| 精品不卡国产一区二区三区| АⅤ资源中文在线天堂| 午夜福利免费观看在线| 亚洲午夜精品一区,二区,三区| 每晚都被弄得嗷嗷叫到高潮| 99国产综合亚洲精品| 欧美色欧美亚洲另类二区 | 18禁观看日本| 91成年电影在线观看| 色尼玛亚洲综合影院| 国产欧美日韩一区二区三区在线| 亚洲av电影在线进入| 狂野欧美激情性xxxx| 女同久久另类99精品国产91| 淫秽高清视频在线观看| 香蕉国产在线看| 久久九九热精品免费| 亚洲在线自拍视频| 中文亚洲av片在线观看爽| 亚洲精华国产精华精| 国产三级黄色录像| 国产极品粉嫩免费观看在线| 无限看片的www在线观看| 18禁国产床啪视频网站| 精品无人区乱码1区二区| 亚洲av成人一区二区三| 国产欧美日韩一区二区三区在线| 亚洲中文av在线| 97超级碰碰碰精品色视频在线观看| 色婷婷久久久亚洲欧美| 在线永久观看黄色视频| 一个人观看的视频www高清免费观看 | 成人精品一区二区免费| 久热爱精品视频在线9| 丁香欧美五月| 十八禁网站免费在线| 亚洲中文av在线| 欧美人与性动交α欧美精品济南到| 老司机福利观看| 一本大道久久a久久精品| 久久人人97超碰香蕉20202| 久久香蕉国产精品| 色精品久久人妻99蜜桃| 亚洲人成伊人成综合网2020| 欧美丝袜亚洲另类 | 国产av在哪里看| av有码第一页| 俄罗斯特黄特色一大片| 久久久精品国产亚洲av高清涩受| 欧美一级毛片孕妇| 黑人欧美特级aaaaaa片| 在线十欧美十亚洲十日本专区| 日韩欧美一区视频在线观看| 亚洲精品中文字幕在线视频| 久久精品91蜜桃| www.精华液| 波多野结衣巨乳人妻| 91成人精品电影| 免费女性裸体啪啪无遮挡网站| 自拍欧美九色日韩亚洲蝌蚪91| 老司机深夜福利视频在线观看| 97人妻精品一区二区三区麻豆 | 黄色视频不卡| 免费看十八禁软件| 亚洲第一欧美日韩一区二区三区| 午夜视频精品福利| 一个人观看的视频www高清免费观看 | 午夜老司机福利片| 亚洲欧美精品综合久久99| 亚洲av成人一区二区三| 涩涩av久久男人的天堂| 美女高潮到喷水免费观看| 两人在一起打扑克的视频| 十分钟在线观看高清视频www| 久久久久国内视频| 久久久久久亚洲精品国产蜜桃av| www国产在线视频色| 在线十欧美十亚洲十日本专区| 男女下面进入的视频免费午夜 | 丁香六月欧美| 亚洲 国产 在线| 亚洲成av片中文字幕在线观看| 如日韩欧美国产精品一区二区三区| 美女 人体艺术 gogo| 中文字幕精品免费在线观看视频| 亚洲中文字幕日韩| www.999成人在线观看| 亚洲精品久久国产高清桃花| 50天的宝宝边吃奶边哭怎么回事| 一边摸一边做爽爽视频免费| 19禁男女啪啪无遮挡网站| 91国产中文字幕| 国产区一区二久久| 国语自产精品视频在线第100页| 亚洲中文字幕日韩| 一个人免费在线观看的高清视频| 亚洲精品美女久久久久99蜜臀| 不卡av一区二区三区| 欧美av亚洲av综合av国产av| 久久久精品欧美日韩精品| 久9热在线精品视频| 天堂影院成人在线观看| 色老头精品视频在线观看| 男女做爰动态图高潮gif福利片 | 国产免费男女视频| 999精品在线视频| 婷婷六月久久综合丁香| 久久亚洲真实| 精品一区二区三区四区五区乱码| 一级毛片高清免费大全| 日韩精品中文字幕看吧| 国产一区在线观看成人免费| 精品人妻在线不人妻| 一个人免费在线观看的高清视频| 亚洲人成电影免费在线| 在线播放国产精品三级| 九色亚洲精品在线播放| 波多野结衣高清无吗| 亚洲片人在线观看| 中亚洲国语对白在线视频| 女人被狂操c到高潮| 91在线观看av| 亚洲熟妇中文字幕五十中出| 狠狠狠狠99中文字幕| 亚洲情色 制服丝袜| 黄色毛片三级朝国网站| 天天躁狠狠躁夜夜躁狠狠躁| 免费在线观看完整版高清| 日韩有码中文字幕| 亚洲国产精品久久男人天堂| 亚洲第一电影网av| 香蕉久久夜色| 亚洲欧美日韩无卡精品| 亚洲精品美女久久av网站| 日韩av在线大香蕉| 免费不卡黄色视频| 亚洲色图综合在线观看| 亚洲精品一卡2卡三卡4卡5卡| 99热只有精品国产| 天天添夜夜摸| 精品日产1卡2卡| 久久 成人 亚洲| www.自偷自拍.com| 99国产精品免费福利视频| 国产欧美日韩综合在线一区二区| 久久久久精品国产欧美久久久| 欧美乱码精品一区二区三区| 在线视频色国产色| 国产一区二区三区在线臀色熟女| 99在线视频只有这里精品首页| 我的亚洲天堂| 黄网站色视频无遮挡免费观看| 中文字幕人妻丝袜一区二区| 午夜亚洲福利在线播放| √禁漫天堂资源中文www| 国产一级毛片七仙女欲春2 | 国产av又大| www日本在线高清视频| 老汉色∧v一级毛片| 天堂√8在线中文| 国产成人精品久久二区二区91| 日本在线视频免费播放| 国产亚洲精品第一综合不卡| 极品教师在线免费播放| 国产日韩一区二区三区精品不卡| 成人av一区二区三区在线看| cao死你这个sao货| 在线国产一区二区在线| 久久久久久免费高清国产稀缺| 黄频高清免费视频| av天堂久久9| 嫩草影视91久久| 91成年电影在线观看| 色av中文字幕| 一级a爱视频在线免费观看| 日本免费一区二区三区高清不卡 | 国产主播在线观看一区二区| 久久久久久亚洲精品国产蜜桃av| 午夜精品国产一区二区电影| 欧美老熟妇乱子伦牲交| 大码成人一级视频| 国产精品影院久久| 欧美在线一区亚洲| 黄色a级毛片大全视频| 一二三四在线观看免费中文在| av视频在线观看入口| 国产亚洲精品一区二区www| 黄色女人牲交| 99久久国产精品久久久| 久久 成人 亚洲| 国产不卡一卡二| 欧美成人一区二区免费高清观看 | 身体一侧抽搐| or卡值多少钱| 久久青草综合色| 色综合亚洲欧美另类图片| 老鸭窝网址在线观看| 又黄又粗又硬又大视频| 国产一区二区在线av高清观看| 国产成人一区二区三区免费视频网站| 女同久久另类99精品国产91| АⅤ资源中文在线天堂| 国产精品永久免费网站| 欧美激情久久久久久爽电影 | 亚洲片人在线观看| 精品熟女少妇八av免费久了| 啦啦啦韩国在线观看视频| 亚洲欧美日韩无卡精品| 色精品久久人妻99蜜桃| 少妇熟女aⅴ在线视频| 国产精品爽爽va在线观看网站 | 欧美不卡视频在线免费观看 | 中亚洲国语对白在线视频| 两个人看的免费小视频| 91成人精品电影| 久久亚洲真实| 精品卡一卡二卡四卡免费| 在线视频色国产色| 亚洲成国产人片在线观看| 自线自在国产av| 两个人视频免费观看高清| 中文字幕av电影在线播放| 亚洲av成人不卡在线观看播放网| 两个人视频免费观看高清| www.熟女人妻精品国产| 成人三级黄色视频| 亚洲五月婷婷丁香| 黄色成人免费大全| 日韩三级视频一区二区三区| 欧美大码av| 好男人电影高清在线观看| www.www免费av| 久久伊人香网站| 亚洲精品一卡2卡三卡4卡5卡| 精品人妻在线不人妻| www国产在线视频色| 人人妻,人人澡人人爽秒播| 亚洲成人国产一区在线观看| 日韩高清综合在线| 真人做人爱边吃奶动态| 精品国产美女av久久久久小说| 满18在线观看网站| 国产伦一二天堂av在线观看| 亚洲av成人av| 夜夜躁狠狠躁天天躁| 精品久久久久久,| 久久久国产成人精品二区| 久久精品91蜜桃| 久久久久国内视频| 精品卡一卡二卡四卡免费| 琪琪午夜伦伦电影理论片6080| 非洲黑人性xxxx精品又粗又长| 久久婷婷成人综合色麻豆| 99在线人妻在线中文字幕| 久久影院123| 久久天堂一区二区三区四区| 日韩一卡2卡3卡4卡2021年| 久久欧美精品欧美久久欧美| 99国产精品一区二区蜜桃av| 成人永久免费在线观看视频| 久久天堂一区二区三区四区| 国产精品 欧美亚洲| 满18在线观看网站| 一本大道久久a久久精品| 欧美日韩乱码在线| 精品日产1卡2卡| 国产成人系列免费观看| 禁无遮挡网站| 每晚都被弄得嗷嗷叫到高潮| 欧美色视频一区免费| 美国免费a级毛片| 国内毛片毛片毛片毛片毛片| 国产精品亚洲美女久久久| 亚洲中文字幕日韩| 亚洲视频免费观看视频| 国产亚洲欧美98| 色在线成人网| 国产av又大| 久久婷婷成人综合色麻豆| 嫩草影院精品99| 精品国产国语对白av| 国产免费av片在线观看野外av| 欧美日韩一级在线毛片| 亚洲精品国产区一区二|