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

    紫坪鋪水庫(kù)水位變化對(duì)剪切波分裂參數(shù)的影響

    2015-06-06 12:40:33劉莎吳朋
    地球物理學(xué)報(bào) 2015年11期
    關(guān)鍵詞:方向

    劉莎, 吳朋

    1 中國(guó)地震局地球物理研究所, 北京 100081 2 四川省地震局, 成都 610041

    ?

    紫坪鋪水庫(kù)水位變化對(duì)剪切波分裂參數(shù)的影響

    劉莎1, 吳朋2

    1 中國(guó)地震局地球物理研究所, 北京 100081 2 四川省地震局, 成都 610041

    本文通過(guò)對(duì)2006—2009年四川紫坪鋪水庫(kù)庫(kù)區(qū)8個(gè)地震臺(tái)站記錄的地震事件,采用剪切波分裂方法獲得了水庫(kù)庫(kù)區(qū)剪切波分裂參數(shù),并結(jié)合地震活動(dòng)性與水庫(kù)水位之間的變化關(guān)系,分析了紫坪鋪水庫(kù)庫(kù)區(qū)地殼應(yīng)力的變化特征.剪切波分裂結(jié)果顯示該研究區(qū)域快波偏振方向有兩個(gè),分別為北東向和北西向,充分體現(xiàn)了紫坪鋪水庫(kù)地區(qū)地殼應(yīng)力是由北西向的區(qū)域主壓應(yīng)力與南東走向的龍門(mén)山斷裂帶綜合作用的結(jié)果.慢波延遲時(shí)間平均值為5.8 ms·km-1,慢波延遲時(shí)間較大的地區(qū)位于庫(kù)壩和庫(kù)尾,分別是水庫(kù)蓄水排水引起地殼應(yīng)力變化最大的區(qū)域.對(duì)比慢波延遲時(shí)間的變化和水庫(kù)水位的變化顯示了慢波延遲時(shí)間與水庫(kù)水位之間的一致變化關(guān)系,揭示了水庫(kù)的蓄水排水對(duì)地殼應(yīng)力的影響.

    快波偏振方向; 慢波延遲時(shí)間; 水位; 地殼應(yīng)力

    1 引言

    近些年來(lái),隨著國(guó)內(nèi)外大型高壩水庫(kù)的不斷興建,水庫(kù)誘發(fā)地震事件也隨之增多.水庫(kù)蓄水或排水通常會(huì)引起地殼巖體中應(yīng)力集中和能量釋放而產(chǎn)生地震(秦嘉政等,2009).這些地震位置一般鄰近重要的水庫(kù)工程設(shè)施,而且震源較淺,震中烈度相對(duì)較高,往往具有很大的破壞性,造成大壩及附近建筑物的破壞和人員傷亡,甚至引起滑坡、坍塌等嚴(yán)重的次生災(zāi)害.因此,水庫(kù)地區(qū)地震的研究不僅是水利水電工程研究的內(nèi)容,也是地震學(xué)研究的重要內(nèi)容之一.

    剪切波分裂作為一個(gè)有效的觀(guān)測(cè)手段,可用于對(duì)地殼應(yīng)力、流體壓力等方面的探索(張永久等,2010; Gao et al., 2011; Tang et al., 2005).剪切波分裂參數(shù)快波偏振方向和慢波延遲時(shí)間是用于分析地殼應(yīng)力場(chǎng)的狀態(tài)和各向異性參數(shù)的靜態(tài)和動(dòng)態(tài)探討的重要依據(jù)(Crampin, 1984; 高原等,1999).剪切波分裂對(duì)地殼中各向異性介質(zhì)較為敏感,與地殼應(yīng)力場(chǎng)和斷裂線(xiàn)的分布密切相關(guān).快波偏振方向與地殼微裂隙排列方向一致,與原地主壓應(yīng)力的方向一致(Crampin, 1978; Zhang et al., 2000; Musumeci et al.,2005).慢波延遲時(shí)間通常是評(píng)價(jià)各向異性程度的物理量,對(duì)地殼中微裂隙的幾何形態(tài)和密度分布非常敏感,可以反映地震發(fā)生前后應(yīng)力的積累與釋放過(guò)程,而且也可以很好地反映地殼中水壓變化所引起的微裂隙變化.

    紫坪鋪水庫(kù)位于四川省都江堰市北西9 km處的岷江上游,是一座壩高156 m庫(kù)容11.12億立方米的大型水利工程.紫坪鋪水庫(kù)已經(jīng)建立了專(zhuān)門(mén)用來(lái)監(jiān)測(cè)水庫(kù)地區(qū)中小地震的數(shù)字遙測(cè)地震臺(tái)網(wǎng),該地震臺(tái)網(wǎng)由7個(gè)臺(tái)站組成.在紫坪鋪水庫(kù)庫(kù)區(qū)還有油榨坪臺(tái)站,屬于四川地震臺(tái)網(wǎng).自2005年9月紫坪鋪水庫(kù)下閘蓄水之后,庫(kù)區(qū)弱震活動(dòng)頻繁,為紫坪鋪水庫(kù)庫(kù)區(qū)地震相關(guān)研究提供了豐富的數(shù)據(jù)資料(胡先明,2007;雷興林等,2008).

    之前有學(xué)者(張永久等,2010)已經(jīng)對(duì)紫坪鋪水庫(kù)庫(kù)區(qū)剪切波分裂進(jìn)行探討,獲得了水庫(kù)庫(kù)區(qū)各向異性特征.該結(jié)果采用的是紫坪鋪水庫(kù)庫(kù)區(qū)2004年8月—2008年5月汶川地震之前的地震數(shù)據(jù),本文選擇2006—2009年紫坪鋪水庫(kù)庫(kù)區(qū)的地震波形資料,獲得水庫(kù)庫(kù)區(qū)的各向異性特征,進(jìn)而分析水庫(kù)水位的變化與地殼應(yīng)力的變化關(guān)系,并對(duì)比張永久等(2010)的結(jié)果,探討兩次剪切波分裂獲得的地震各向異性特征有何異同.

    2 地震和水位數(shù)據(jù)

    紫坪鋪水庫(kù)位于龍門(mén)山斷裂帶中段,龍門(mén)山斷裂帶處于中國(guó)大陸南北地震帶的中部,東南側(cè)為揚(yáng)子準(zhǔn)地臺(tái),西側(cè)為松潘—甘孜地塊.地震震源機(jī)制解結(jié)果表明龍門(mén)山斷裂帶所處地區(qū)的現(xiàn)代應(yīng)力場(chǎng)主壓應(yīng)力軸方位為北西西向(鄧起東等,1994).紫坪鋪水庫(kù)自2005年9月下閘蓄水之后弱震活動(dòng)較為活躍,為剪切波分裂研究提供了豐富的數(shù)據(jù)資料.本文選取了2006—2009年紫坪鋪水庫(kù)庫(kù)區(qū)8個(gè)地震臺(tái)站記錄的地震數(shù)據(jù),包括紫坪鋪水庫(kù)數(shù)字地震臺(tái)網(wǎng)的7個(gè)臺(tái)站(八角BAJ,白巖BAY,桂花樹(shù)GHS,靈巖寺LYS,廟子坪MZP,桃子坪TZP,鉆洞子ZDZ)和四川數(shù)字地震臺(tái)網(wǎng)的油榨坪YZP臺(tái)站,8個(gè)臺(tái)站平均臺(tái)間距10 km左右,采樣率為100 sps (samples per second)(圖1).

    在2006—2009年整個(gè)觀(guān)測(cè)時(shí)間段內(nèi),紫坪鋪水庫(kù)水位起伏較大,水庫(kù)庫(kù)區(qū)弱震活動(dòng)較強(qiáng)(圖2).M-T圖顯示水庫(kù)庫(kù)區(qū)沒(méi)有大于4.0級(jí)的地震,地震震級(jí)主要集中在3.0級(jí)以下.而且水庫(kù)庫(kù)區(qū)地震活動(dòng)與水庫(kù)水位的升降變化有一定的相關(guān)性.2006年上半年水庫(kù)水位處于830 m左右,8月13日開(kāi)始,水庫(kù)水位開(kāi)始出現(xiàn)大幅度上升,到10月7日累計(jì)上升54 m,此后水位開(kāi)始下降,地震活動(dòng)從2006年8月開(kāi)始逐漸增強(qiáng).2007年5月20日水庫(kù)水位再次上升,12月14日水位達(dá)到高值后開(kāi)始下降,直到2008年6月降至820 m.2008年2月庫(kù)區(qū)地震活動(dòng)性較強(qiáng).2008年11月,水庫(kù)水位再次達(dá)到峰值隨后下降.2009年6月水庫(kù)水位重新開(kāi)始上升,地震活動(dòng)性在2009年下半年也表現(xiàn)出逐漸增強(qiáng)的趨勢(shì).

    3 剪切波分裂方法

    本文采用偏振分析法獲得剪切波分裂參數(shù),偏振分析法雖然計(jì)算速度較慢,但是準(zhǔn)確率較高,直觀(guān)性強(qiáng)(Crampin et al., 2006).

    剪切波分裂方法研究地殼各向異性要求選取的地震波記錄必須在剪切波窗口內(nèi).利用剪切波分裂探討地殼介質(zhì)的各向異性,還要求地震數(shù)據(jù)有較高的信噪比.根據(jù)淺層勘探給出的速度模型(表1)(趙珠等,1997),該地區(qū)深度12 km,震中距20 km的地震入射角大約58°時(shí),實(shí)際地震射線(xiàn)的入射角僅35°.因此本文選取58°的入射角作為剪切波窗口.對(duì)數(shù)據(jù)進(jìn)行進(jìn)一步的篩選,選取地震事件震源深度在5~12 km范圍內(nèi),而且剪切波信噪比大于6.0.最后符合剪切波分裂計(jì)算條件的地震數(shù)據(jù)共512條(表2).

    圖1 紫坪鋪水庫(kù)庫(kù)區(qū)地震(2006—2009年)和臺(tái)站分布Fig.1 Distribution of earthquakes (2006—2009) and stations in Zipingpu reservoir

    圖2 紫坪鋪水庫(kù)庫(kù)區(qū)水庫(kù)水位變化(a),地震M-T(b)及N-T(c)圖Fig.2 The water level (a), M-T (b) and N-T (c) of Zipingpu reservoir

    圖3 BAJ(八角臺(tái))剪切波分裂分析示意圖

    圖4 剪切波偏振方向等面積極射投影與等面積玫瑰圖

    表1 剪切波速度模型Table 1 The velocity model of shear wave

    在確定剪切波分裂參數(shù)時(shí),文中選取剪切波到時(shí)前20個(gè)采樣點(diǎn)和到時(shí)后40個(gè)采樣點(diǎn)的波形數(shù)據(jù)窗口(時(shí)長(zhǎng)0.6 s),將這60個(gè)采樣點(diǎn)數(shù)據(jù)分為三段分別繪制質(zhì)點(diǎn)偏振圖進(jìn)行分析(圖3).畫(huà)出所有滿(mǎn)足數(shù)據(jù)要求的質(zhì)點(diǎn)偏振圖,然后從這些圖中挑選出快剪切波質(zhì)點(diǎn)振動(dòng)為線(xiàn)性或近于線(xiàn)性的質(zhì)點(diǎn)偏振圖.由于各向異性介質(zhì)的存在,剪切波偏振方向第一個(gè)突然轉(zhuǎn)變的位置被認(rèn)為是快剪切波的到達(dá),慢剪切波的到達(dá)同樣也是極化偏振方向的突然轉(zhuǎn)變,快、慢剪切波到時(shí)差即為慢波延遲時(shí)間(圖3).

    4 快波偏振方向結(jié)果

    快波偏振方向通常反映了原地主壓應(yīng)力的方向,與地殼中應(yīng)力誘發(fā)的定向排列的微裂隙面的走向一致(賴(lài)院根等,2006;高原等,1995;石玉濤等,2009).從表2可以看出,8個(gè)臺(tái)站均得到了30條以上的可靠結(jié)果,特別是BAJ臺(tái)有172個(gè)結(jié)果.剪切波分裂結(jié)果統(tǒng)計(jì)8個(gè)臺(tái)站的快波偏振方向均有兩個(gè),本文分別計(jì)算剪切波分裂參數(shù)平均值以及對(duì)應(yīng)的均方根誤差.快波偏振方向的誤差結(jié)果都小于30°,表明各個(gè)臺(tái)站的快波偏振方向優(yōu)勢(shì)取向都較為明顯.結(jié)果顯示研究區(qū)域內(nèi)快波偏振方向主要有兩個(gè),分別為北西向和北東向.

    表2 四川紫坪鋪水庫(kù)庫(kù)區(qū)地震臺(tái)站剪切波分裂參數(shù)Table 2 The parameters of shear wave splitting in Zipingpu reservoir

    快波偏振方向等面積投影圖(圖4)顯示,研究區(qū)域內(nèi)8個(gè)臺(tái)站快波偏振方向有兩個(gè),但是其優(yōu)勢(shì)取向都較為明顯.BAJ臺(tái),BAY臺(tái),GHS臺(tái),LYS臺(tái),MZP臺(tái),TZP臺(tái)和ZDZ臺(tái)7個(gè)臺(tái)站快波偏振方向顯示為北西向和北東向,而YZP臺(tái)的結(jié)果為北東向和近似的南北向.北東向的快波偏振方向與龍門(mén)山斷裂帶的走向一致,北西向的快波偏振方向與研究區(qū)域主壓應(yīng)力的方向一致.

    5 慢波延遲時(shí)間與水庫(kù)水位

    慢波延遲時(shí)間對(duì)地殼中微裂隙的幾何形態(tài)較為敏感,因此它在一定程度上反映了地殼中應(yīng)力的分布特點(diǎn)(Gao and Crampin, 2008; Liu et al., 2014).研究區(qū)域內(nèi)512條地震事件的慢波延遲時(shí)間平均值為5.8 ms·km-1.慢波延遲時(shí)間的最大值位于八角臺(tái)(BAJ),兩個(gè)方向上的慢波延遲時(shí)間分別是(7.7±3.8) ms·km-1,(8.1±3.1) ms·km-1.其次是靈巖寺臺(tái)(LYS),兩個(gè)方向上的慢波延遲時(shí)間分別是(6.5±2.9) ms·km-1,(5.4±2.0) ms·km-1.最小值位于油榨坪臺(tái)(YZP),兩個(gè)方向上的慢波延遲時(shí)間分別是(3.8±1.6) ms·km-1,(4.2±1.7) ms·km-1.其余5個(gè)臺(tái)站慢波延遲時(shí)間多集中在5 ms·km-1左右(表2).慢波延遲時(shí)間較大的靈巖寺臺(tái)和八角臺(tái)分別位于紫坪鋪水庫(kù)的庫(kù)壩區(qū)和庫(kù)尾區(qū)(圖1),正好是處于水庫(kù)蓄水引起地殼應(yīng)力變化較大的兩個(gè)區(qū)域.慢波延遲時(shí)間最小的油榨坪臺(tái)距離水庫(kù)最遠(yuǎn),地殼應(yīng)力受水庫(kù)水位變化影響最小.可見(jiàn),水庫(kù)的庫(kù)壩和庫(kù)尾地區(qū)地殼應(yīng)力受水庫(kù)水位的影響較大.

    由于慢波延遲時(shí)間的離散性,本文采用 7點(diǎn)滑動(dòng)平均來(lái)表示慢波延遲時(shí)間在時(shí)間上的變化.慢波延遲時(shí)間的7點(diǎn)滑動(dòng)平均曲線(xiàn)顯示慢波延遲時(shí)間變化較為明顯,特別是BAJ,BAY, GHS和LYS臺(tái)站(圖5).BAJ臺(tái)站慢波延遲時(shí)間從2006年初開(kāi)始升高直到10月底升至12 ms·km-1,隨后慢波延遲時(shí)間開(kāi)始下降直到2007年4月.之后,慢波延遲時(shí)間再次升高直到2007年10月.2007年年底,慢波延遲時(shí)間一直維持在11 ms·km-1,隨后降低至7 ms·km-1.2008年2月,慢波延遲時(shí)間再次升高直到2008年10月,12月慢波延遲時(shí)間又一次降低,2009年7月降至4 ms·km-1,隨后開(kāi)始持續(xù)升高.BAY臺(tái)站在2006年和2007年,慢波延遲時(shí)間變化不明顯,2008年年初開(kāi)始稍有降低.2008年7月慢波延遲時(shí)間上升直到2009年初,隨后開(kāi)始降低,7月降至最低3.5 ms·km-1.GHS臺(tái)站慢波延遲時(shí)間也可以觀(guān)測(cè)到變化的趨勢(shì),在2006年底和2008年11月分別可以觀(guān)測(cè)到慢波延遲時(shí)間的兩個(gè)峰值,達(dá)到7 ms·km-1.在LYS臺(tái)站,慢波延遲時(shí)間的變化較小,在2007年底慢波延遲時(shí)間出現(xiàn)低谷,低至4 ms·km-1.

    圖5 慢波延遲時(shí)間與水庫(kù)水位對(duì)比圖

    慢波延遲時(shí)間對(duì)地殼應(yīng)力場(chǎng)所引起的微裂隙變化非常敏感.紫坪鋪水庫(kù)庫(kù)區(qū)慢波延遲時(shí)間在觀(guān)測(cè)時(shí)間段內(nèi)表現(xiàn)出了很好的一致性變化.為了分析水庫(kù)水位對(duì)地殼應(yīng)力變化的影響,文中將慢波延遲時(shí)間變化與水庫(kù)水位變化進(jìn)行對(duì)比分析(圖5).結(jié)果顯示慢波延遲時(shí)間的變化和水庫(kù)水位之間存在很好的一致性.每當(dāng)水庫(kù)水位出現(xiàn)大幅度上升,慢波延遲時(shí)間也同步出現(xiàn)相應(yīng)的增大;水位下降之后,慢波延遲時(shí)間隨之減小.這一結(jié)果與之前的結(jié)果(張永久等,2010)較為吻合.水庫(kù)水位與慢波延遲時(shí)間的一致變化關(guān)系在Gao和Crampin(2003)關(guān)于慢剪切波時(shí)間延遲變化與壓力變化的試驗(yàn)中也得到了證實(shí).水庫(kù)的蓄水與排水通常會(huì)引起地殼流體壓力的變化,剪切波分裂慢波延遲時(shí)間可能會(huì)由于地殼應(yīng)力的微小變化而產(chǎn)生顯著的變化,該結(jié)論也與Crampin和Zatsepin(1997)所提出的APE(Anisotropic Poro-elasticity, 各向異性孔隙彈性)理論相吻合.

    6 地震各向異性的討論

    紫坪鋪水庫(kù)位于龍門(mén)山斷裂帶上,龍門(mén)山斷裂帶作為2008年汶川地震的發(fā)震斷裂,在構(gòu)造上極其復(fù)雜.快波偏振方向通常與區(qū)域主壓應(yīng)力方向或者活動(dòng)斷裂的走向表現(xiàn)一致.研究區(qū)域中,8個(gè)臺(tái)站均顯示有兩個(gè)快波偏振方向,BAJ臺(tái),BAY臺(tái),GHS臺(tái),LYS臺(tái),MZP臺(tái),TZP臺(tái)和 ZDZ臺(tái)7個(gè)臺(tái)站剪切波分裂快波偏振方向顯示北東和北西兩個(gè)方向,而YZP臺(tái)快波偏振方向?yàn)楸睎|向和近似南北向.YZP臺(tái)站北東向的快波偏振方向與斷裂帶的走向一致,南北向的快波偏振方向與其他臺(tái)站結(jié)果存在較大差異,與斷裂帶走向和最大主壓應(yīng)力方向也不一致.南北向的快波偏振方向結(jié)果與張永久等(2010)的結(jié)果相一致,證實(shí)了本文結(jié)果的可靠性.YZP臺(tái)站位于紫坪鋪水庫(kù)東南方向,慢波延遲時(shí)間的結(jié)果顯示該臺(tái)站受到水庫(kù)影響較小,該臺(tái)站南北向快波偏振方向反映了YZP臺(tái)站附近復(fù)雜的應(yīng)力場(chǎng),該區(qū)域應(yīng)力方向并不僅僅是受到龍門(mén)山活動(dòng)斷裂帶單一的影響,也暗示了附近隱伏斷層的情況.以往的研究結(jié)果(高原等,1995;賴(lài)院根等,2006;石玉濤等2009)也表明快波偏振方向也揭示了新的構(gòu)造.

    對(duì)比張永久等(2010)的結(jié)果(圖6),BAJ臺(tái),LYS臺(tái),ZDZ臺(tái),YZP臺(tái)和BAY臺(tái)快波偏振方向均有兩個(gè)方向.本文的結(jié)果中MZP臺(tái),TZP臺(tái)和GHS臺(tái)仍然有兩個(gè)快波偏振方向,而之前的結(jié)果(張永久等,2010)顯示只有一個(gè)快波偏振方向.本文的地震數(shù)據(jù)采用的是2006—2009年,跨越汶川地震.通常大地震的發(fā)生會(huì)使得震區(qū)地殼應(yīng)力進(jìn)行調(diào)整,進(jìn)而快波偏振方向表現(xiàn)的較為零散,所以本文觀(guān)測(cè)到這三個(gè)臺(tái)站快波偏振方向有兩個(gè).較之前的結(jié)果,MZP臺(tái)和TZP臺(tái)快波偏振方向多了北東向,與龍門(mén)山活動(dòng)斷裂的走向一致.而GHS臺(tái)多了北西向的快波偏振方向.

    紫坪鋪水庫(kù)8個(gè)地震臺(tái)站慢波延遲時(shí)間平均值為5.8 ms·km-1,與張永久等(2010)的結(jié)果5.38 ms·km-1近似.位于庫(kù)尾的白巖臺(tái)(BAY)和庫(kù)壩的靈巖寺臺(tái)(LYS)慢波延遲時(shí)間明顯大于其他幾個(gè)臺(tái)站.水庫(kù)的庫(kù)壩和庫(kù)尾相比于水庫(kù)其他地區(qū)地殼水壓較大,使得微裂隙孔隙流體壓力增強(qiáng),慢波延遲時(shí)間將表現(xiàn)出較大的值,這一結(jié)果表明紫坪鋪水庫(kù)庫(kù)壩和庫(kù)尾附近是地殼應(yīng)力相對(duì)較大的區(qū)域.對(duì)比慢波延遲時(shí)間與水庫(kù)水位之間的關(guān)系表明兩者之間的變化有一定的相關(guān)性,水庫(kù)水位的增加往往伴隨著慢波延遲時(shí)間的增大,這是由于水庫(kù)蓄水導(dǎo)致了地殼中孔隙壓力的擴(kuò)散(Talwani,1997;陶瑋等,2014),使得地殼裂隙孔隙壓強(qiáng)增大,進(jìn)而觀(guān)測(cè)到慢波延遲時(shí)間的增大.

    圖6 紫坪鋪水庫(kù)庫(kù)區(qū)快剪切波偏振方向

    本文使用的地震數(shù)據(jù)資料跨越汶川地震,對(duì)比汶川地震之前的剪切波分裂結(jié)果(張永久等,2010)發(fā)現(xiàn),紫坪鋪庫(kù)區(qū)8個(gè)地震臺(tái)站快波偏振方向結(jié)果基本一致,慢波延遲時(shí)間反映的各向異性程度也相當(dāng).水庫(kù)周?chē)穆ㄑ舆t時(shí)間主要受到水庫(kù)水位的影響,與地殼水壓的變化相關(guān).本文得到的剪切波分裂結(jié)果初步反映了水庫(kù)庫(kù)區(qū)的區(qū)域應(yīng)力情況以及水庫(kù)的蓄水排水對(duì)地殼應(yīng)力的影響.

    7 結(jié)論

    紫坪鋪水庫(kù)自2005年9月下閘蓄水之后,水庫(kù)庫(kù)區(qū)弱震活動(dòng)頻繁. 而且紫坪鋪水庫(kù)已經(jīng)建立了專(zhuān)門(mén)用于監(jiān)測(cè)水庫(kù)地區(qū)中小地震的數(shù)字遙測(cè)地震臺(tái)網(wǎng),積累了大量的地震數(shù)據(jù)資料. 本文利用2006—2009年紫坪鋪水庫(kù)庫(kù)區(qū)8個(gè)地震臺(tái)站記錄的地震事件,采用剪切波分裂方法分析了水庫(kù)庫(kù)區(qū)地殼應(yīng)力的分布特征,同時(shí)也探討了水庫(kù)水位的變化對(duì)剪切波分裂參數(shù)的影響.

    紫坪鋪水庫(kù)庫(kù)區(qū)剪切波分裂快波偏振方向主要有兩個(gè),北東向和北西向,反映了研究區(qū)域內(nèi)地殼各向異性的分布特征. 本文的結(jié)果與之前的剪切波分裂結(jié)果表現(xiàn)一致,紫坪鋪水庫(kù)地區(qū)快波偏振方向和區(qū)域主壓應(yīng)力的方向、龍門(mén)山斷裂的走向一致. 紫坪鋪水庫(kù)8個(gè)臺(tái)站慢波延遲時(shí)間平均值為5.8 ms·km-1. 水庫(kù)的庫(kù)壩和庫(kù)尾相比于水庫(kù)其他地區(qū)慢波延遲時(shí)間要大,表明庫(kù)壩和庫(kù)尾是地殼應(yīng)力相對(duì)較大的地區(qū).水庫(kù)庫(kù)區(qū)慢波延遲時(shí)間與水庫(kù)水位之間表現(xiàn)出了很好的一致性變化. 水庫(kù)的蓄水和排水通常會(huì)引起地殼水壓的變化,進(jìn)而使得地殼裂隙孔隙壓強(qiáng)發(fā)生相應(yīng)的變化.

    致謝 感謝四川省地震局為本文的研究提供數(shù)據(jù)的支持和幫助;感謝兩位匿名評(píng)審專(zhuān)家提出的寶貴意見(jiàn).

    Crampin S. 1978. Seismic-wave propagation through a cracked solid: polarization as a possible dilatancy diagnostic.Geophys.J.Int., 53(3): 467-496. Crampin S. 1984. An introduction to wave propagation in anisotropic media.Geophys.J.Int., 76(1): 17-28.

    Crampin S, Gao Y. 2006. A review of techniques for measuring shear-wave splitting above small earthquakes.Phys.Earth.Planet.Inter., 159: 1-14.

    Crampin S, Zatsepin S V. 1997. Modelling the compliance of crustal rock—II. Response to temporal changes before earthquakes.Geophys.J.Int., 129(3): 495-506.

    Deng Q D, Chen S F, Zhao X L. 1994. Tectonics, seismicity and dynamics of Longmenshan Mountains and its adjacent regions.SeismologicalandGeology(in Chinese), 16(4): 389-403.

    Gao Y, Zheng S H, Sun Y. 1995. Crustal crack anisotropy in Tangshan.ActaSeismologicaSinica(in Chinese), 17(3): 283-293.

    Gao Y, Zheng S H, Zhou H L. 1999. Polarization patterns of fast shear wave in Tangshan region and their variations.ChineseJ.Geophys. (in Chinese), 42(2): 228-232.

    Gao Y, Crampin S. 2003. Temporal variations of shear-wave splitting in field and laboratory studies in China.J.Appl.Geophys., 54(3-4): 279-287.

    Gao Y, Crampin S. 2008. Shear-wave splitting and earthquake forecasting.TerraNova, 20(6): 440-448.

    Gao Y, Wu J, Fukao Y, et al. 2011. Shear wave splitting in the crust in North China: stress, faults and tectonic implications.Geophys.J.Int., 187(2): 642-654.

    Hu X M. 2007. Seismicity background around Zipingpu Dam before impounding.EarthquakeResearchinSichuan(in Chinese), (2): 16-21.

    Lai Y G, Liu Q Y, Chen J H, et al. 2006. Shear wave splitting and the features of the crustal stress field in the Capital circle.ChineseJ.Geophys. (in Chinese), 49(1): 189-196.

    Lei X L, Ma S L, Wen X Z, et al. 2008. Integrated analysis of stress and regional seismicity by surface loading—A case study of Zipingpu reservoir.SeismologyandGeology(in Chinese), 30(4): 1046-1064.

    Liu S, Crampin S, Luckett R, et al. 2014. Changes in shear wave splitting before the 2010 Eyjafjallaj?kull eruption in Iceland.Geophys.J.Int., 199(1): 102-112.

    Musumeci C, Patanè D, Scarfi L, et al. 2005. Stress directions and shear-wave anisotropy: observations from local earthquakes in southeastern Sicily, Italy.Bull.Seismol.Soc.Am., 95(4): 1395-1374.

    Qin J Z, Liu L F, Qian X D. 2009. Research on characteristics and prediction of reservoir induced seismicity.JournalofSeismologicalResearch(in Chinese), 32(2): 105-113.

    Shi Y T, Gao Y, Zhao C P, et al. 2009. A study of seismic anisotropy of Wenchuan earthquake sequence.ChineseJ.Geophys. (in Chinese), 52(2): 398-407.

    Talwani P. 1997. On the nature of reservoir-induced seismicity.PureAppl.Geophys., 150(3-4): 473-492, doi: 10.1007/s000240050089.

    Tang C H, Rial J A, Lees J M. 2005. Shear-wave splitting: a diagnostic tool to monitor fluid pressure in geothermal fields.Geophys.Res.Lett., 32(21): L21317, doi: 10.1029/2005GL023551.

    Tao W, Masterlark T, Shen Z K, et al. 2014. Triggering effect of the Zipingpu Reservoir on the 2008MW7.9 Wenchuan, China, Earthquake due to poroelastic coupling.ChineseJ.Geophys. (in Chinese), 57(10): 3318-3331, doi: 10.6038/cjg20141019.

    Zhang Y J, Gao Y, Shi Y T, et al. 2010. The shear-wave splitting study of Sichuan Zipingpu reservoir region.ChineseJ.Geophys. (in Chinese), 53(9): 2091-2101, doi: 10.3969/j.issn.0001-5733.2010.09.009.

    Zhang Z J, Li Y K, Lu D Y, et al. 2000. Velocity and anisotropy structure of the crust in the Dabieshan orogenic belt from wide-angle seismic data.Phys.EarthPlanet.Inter., 122(1-2): 115-131.

    Zhao Z, Fan J, Zheng S H, et al. 1997. Precision determination of the crustal structure and hypocentral locations in the Longmenshan thrust nappe belt.ActaSeismologicaSinica(in Chinese), 19(6): 615-622.

    附中文參考文獻(xiàn)

    鄧起東, 陳社發(fā), 趙小麟. 1994. 龍門(mén)山及其鄰區(qū)的構(gòu)造和地震活

    動(dòng)及動(dòng)力學(xué). 地震地質(zhì), 16(4): 389-403.

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

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

    胡先明. 2007. 紫坪鋪水庫(kù)蓄水前天然地震活動(dòng). 四川地震, (2): 16-21.

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

    雷興林, 馬勝利, 聞學(xué)澤等. 2008. 地表水體對(duì)斷層應(yīng)力與地震時(shí)空分布影響的綜合分析——以紫坪鋪水庫(kù)為例. 地震地質(zhì), 30(4): 1046-1064.

    秦嘉政, 劉麗芳, 錢(qián)曉東. 2009. 水庫(kù)誘發(fā)地震活動(dòng)特征及其預(yù)測(cè)方法研究. 地震研究, 32(2): 105-113.

    石玉濤, 高原, 趙翠萍等. 2009. 汶川地震余震序列的地震各向異性. 地球物理學(xué)報(bào), 52(2): 398-407.

    陶瑋, Masterlark T, 沈正康等. 2014. 紫坪鋪水庫(kù)造成孔隙彈性耦合變化及其對(duì)2008年汶川地震觸發(fā)作用. 地球物理學(xué)報(bào), 57(10): 3318-3331, doi: 10.6038/cjg20141019.

    張永久, 高原, 石玉濤等. 2010. 四川紫坪鋪水庫(kù)庫(kù)區(qū)地震剪切波分裂研究. 地球物理學(xué)報(bào), 53(9): 2091-2101, doi: 10.3969/j.issn.0001-5733.2010.09.009.

    趙珠, 范軍, 鄭斯華等. 1997. 龍門(mén)山斷裂帶地殼速度結(jié)構(gòu)和震源位置的精確修定. 地震學(xué)報(bào), 19(6): 615-622.

    (本文編輯 何燕)

    The effect of water level changes in Zipingpu reservoir on the parameters of shear wave splitting

    LIU Sha1, WU Peng2

    1InstituteofGeophysics,ChinaEarthquakeAdministration,Beijing100081,China2EarthquakeAdministrationofSichuanProvince,Chengdu610041,China

    The water storage and drainage of the reservoir usually cause the stress change in the crust, which can induce earthquakes. The Zipingpu reservoir is located in the middle section of Longmenshan fault with complex structure. The digital seismic network has been established in the reservoir to monitor moderate and small earthquakes, and accumulated a wealth of data, facilitating the study of earthquakes in the reservoir.

    Shear wave splitting can be used to detect the stress change in the crust, and obtain the anisotropic static and dynamic gauge. The polarization direction of fast shear waves is consistent with the direction of cracks in the crust, indicating the regional stress direction in the crust. The time delay of the slow shear waves is related to the distribution of the micro-cracks in the crust which implies the change of water pressure in the crust. And the delay time can be used to study the effect of reservoir water level on the crustal anisotropy. This paper obtained the shear wave splitting parameters by the polarization analysis with high accuracy and strong visuality, using seismic data recorded by 8 seismic stations in the Zipingpu reservoir area from 2006 to 2009. According to the velocity model of the shallow exploration, the shear wave window chosen in this paper is 58 degree. The focal depths of the earthquakes are in the range of 5~12 km, and the signal to noise ratio of shear waves is greater than 6.0. 512 seismic events are justified by the data selection criteria. On the basis of the particle moving figures, we choose the position where the particle trajectory deflects vertically to ascertain the arrival time of fast and slow shear waves.

    Every seismic station recorded more than 30 events in the study area, and the polarization direction at each station is obvious. Combining with the relationship between the seismic activity and water level of the Zipingpu reservoir, this paper obtains the variation characteristics of the stress in the crust. The results of shear wave splitting show two polarization directions of fast shear waves, which are northeast and northwest. The northeast polarization direction of fast shear waves is consistent with the strike of the Longmenshan fault. And the northwest direction agrees with the direction of the regional principal compressive stress. The average delay time of the slow shear waves in the reservoir is 5.8 ms·km-1. The large delay time of the slow shear waves is from the station BAJ and LYS at the dam and reservoir tail, respectively, where the stress change is large because of the water storage and drainage in the reservoir. The small delay time is from the station YZP far from the reservoir, where the stress in the crust is affected least by the water level of the reservoir.

    Comparison with the previous results of shear wave splitting in the Zipingpu reservoir shows agreement with our results. The polarization direction of fast shear waves is consistent with the direction of the regional principal compressive stress and the strike of active faults. The relationship between the delay time of the slow shear waves and the water level of the reservoir suggests the increase of water level in the reservoir is often accompanied with the increase of the delay time. The reservoir storage causes the diffusion of the pore pressure in the crust, which makes the aspect ratio of the microcracks in the crust increase.

    Polarization direction of fast shear waves; Delay time of slow shear waves; Water level; Stress in the crust

    10.6038/cjg20151118.Liu S, Wu P. 2015. The effect of water level changes in Zipingpu reservoir on the parameters of shear wave splitting.ChineseJ.Geophys. (in Chinese),58(11):4106-4114,doi:10.6038/cjg20151118.

    中國(guó)地震局基本科研業(yè)務(wù)費(fèi)(DQJB13B07),“2013年臺(tái)站三結(jié)合——四川地區(qū)S波分裂研究”(2013-48)共同資助.

    劉莎,女,1985年生,助理研究員,主要從事地球介質(zhì)地震各向異性研究.E-mail:liusha@cea-igp.ac.cn

    10.6038/cjg20151118

    P315

    2015-04-28,2015-09-28收修定稿

    劉莎,吳朋. 2015. 紫坪鋪水庫(kù)水位變化對(duì)剪切波分裂參數(shù)的影響.地球物理學(xué)報(bào),58(11):4106-4114,

    猜你喜歡
    方向
    2023年組稿方向
    方向
    青年運(yùn)動(dòng)的方向(節(jié)選)
    2022年組稿方向
    2022年組稿方向
    2021年組稿方向
    如何確定位置與方向
    2021年組稿方向
    2021年組稿方向
    大自然中的方向
    人妻一区二区av| 亚洲欧美日韩另类电影网站 | 国产男女超爽视频在线观看| 精品国产乱码久久久久久小说| 伦理电影免费视频| 精品人妻熟女av久视频| 国产在线男女| 亚洲精品乱码久久久久久按摩| 日韩大片免费观看网站| 精品久久久久久久末码| videos熟女内射| 欧美日韩视频高清一区二区三区二| 一级毛片我不卡| 久热久热在线精品观看| 内射极品少妇av片p| 亚洲一区二区三区欧美精品| 一级毛片电影观看| 麻豆精品久久久久久蜜桃| 精品久久久噜噜| 美女脱内裤让男人舔精品视频| 亚洲欧美日韩无卡精品| 噜噜噜噜噜久久久久久91| 亚洲国产精品成人久久小说| 亚洲国产毛片av蜜桃av| 黑丝袜美女国产一区| 久久精品国产自在天天线| 国产综合精华液| 99热6这里只有精品| 国产淫语在线视频| 男女边摸边吃奶| 国产亚洲欧美精品永久| 在现免费观看毛片| 又黄又爽又刺激的免费视频.| 波野结衣二区三区在线| 亚洲av男天堂| 国产精品久久久久成人av| 国产乱人偷精品视频| 女性生殖器流出的白浆| 日本vs欧美在线观看视频 | 成人二区视频| 少妇高潮的动态图| 亚洲欧美成人综合另类久久久| 久久久久久久精品精品| 久久久久久人妻| 色综合色国产| 婷婷色综合大香蕉| 亚洲,一卡二卡三卡| 亚洲av综合色区一区| 日日摸夜夜添夜夜添av毛片| 亚洲精品久久久久久婷婷小说| 丰满乱子伦码专区| 性色av一级| 日韩在线高清观看一区二区三区| 日本色播在线视频| 91在线精品国自产拍蜜月| 亚洲精品中文字幕在线视频 | 人人妻人人澡人人爽人人夜夜| 精品久久久久久久久亚洲| 天堂8中文在线网| 青春草国产在线视频| 久久毛片免费看一区二区三区| 精品人妻偷拍中文字幕| 美女xxoo啪啪120秒动态图| 99热6这里只有精品| 亚洲av中文字字幕乱码综合| 亚洲精品,欧美精品| 激情 狠狠 欧美| 99热6这里只有精品| 蜜桃在线观看..| 亚洲精品久久久久久婷婷小说| 久久影院123| 日韩一区二区三区影片| 色综合色国产| 国产男人的电影天堂91| 秋霞在线观看毛片| 亚洲欧美成人精品一区二区| 久久精品久久久久久噜噜老黄| 亚洲,一卡二卡三卡| 成年av动漫网址| 看十八女毛片水多多多| 在线观看三级黄色| 秋霞在线观看毛片| 伦精品一区二区三区| 国产真实伦视频高清在线观看| 日本爱情动作片www.在线观看| 91精品国产九色| 亚洲国产日韩一区二区| www.色视频.com| 性色avwww在线观看| 国产精品麻豆人妻色哟哟久久| 精品熟女少妇av免费看| 国产欧美另类精品又又久久亚洲欧美| 久久久久国产精品人妻一区二区| 人妻制服诱惑在线中文字幕| 十分钟在线观看高清视频www | 国产乱来视频区| 久久亚洲国产成人精品v| 亚洲欧美精品专区久久| 亚洲aⅴ乱码一区二区在线播放| 中文字幕免费在线视频6| 久久女婷五月综合色啪小说| 日日撸夜夜添| 欧美xxxx性猛交bbbb| 久久毛片免费看一区二区三区| 自拍偷自拍亚洲精品老妇| 超碰av人人做人人爽久久| 干丝袜人妻中文字幕| 亚洲最大成人中文| 一本色道久久久久久精品综合| 老司机影院成人| 成年女人在线观看亚洲视频| 久久精品人妻少妇| 午夜福利视频精品| 国产免费又黄又爽又色| 中文欧美无线码| 久久99热这里只频精品6学生| 中文字幕免费在线视频6| 精品国产一区二区三区久久久樱花 | 午夜免费鲁丝| 久久久久久人妻| 插逼视频在线观看| 青春草亚洲视频在线观看| 香蕉精品网在线| 少妇高潮的动态图| 国产精品免费大片| 中文字幕免费在线视频6| 精品国产一区二区三区久久久樱花 | 日韩免费高清中文字幕av| 亚洲第一av免费看| 精品人妻偷拍中文字幕| 免费观看无遮挡的男女| 中文字幕久久专区| 日本午夜av视频| 亚洲婷婷狠狠爱综合网| 国产又色又爽无遮挡免| 搡女人真爽免费视频火全软件| 国模一区二区三区四区视频| 一级毛片aaaaaa免费看小| 久久精品人妻少妇| 亚洲高清免费不卡视频| 久久毛片免费看一区二区三区| 精品人妻视频免费看| 美女内射精品一级片tv| 欧美高清成人免费视频www| 哪个播放器可以免费观看大片| 美女内射精品一级片tv| kizo精华| 十分钟在线观看高清视频www | 极品教师在线视频| 日本黄大片高清| 亚洲精品国产成人久久av| 亚洲成人一二三区av| 黄色日韩在线| a级毛片免费高清观看在线播放| 亚洲精品456在线播放app| 久久久久久久久久人人人人人人| 男女国产视频网站| 特大巨黑吊av在线直播| 我的老师免费观看完整版| 亚洲精品国产成人久久av| 欧美 日韩 精品 国产| av福利片在线观看| 18禁在线播放成人免费| 日本欧美视频一区| 黄片wwwwww| 亚洲精品国产av成人精品| h视频一区二区三区| 亚洲熟女精品中文字幕| a 毛片基地| 在线精品无人区一区二区三 | h视频一区二区三区| 国产一区二区三区av在线| 午夜老司机福利剧场| 成人午夜精彩视频在线观看| 亚洲av男天堂| 国产精品伦人一区二区| 三级国产精品欧美在线观看| 亚洲欧美日韩卡通动漫| 久热这里只有精品99| 欧美日韩在线观看h| 欧美xxxx黑人xx丫x性爽| 男女无遮挡免费网站观看| 中文欧美无线码| 一级毛片我不卡| 另类亚洲欧美激情| 黑人高潮一二区| 韩国高清视频一区二区三区| 国产精品一区二区性色av| 国产综合精华液| 精品亚洲成国产av| 亚洲精品视频女| 男男h啪啪无遮挡| 成人综合一区亚洲| 九九久久精品国产亚洲av麻豆| av女优亚洲男人天堂| 亚洲国产日韩一区二区| 精品一区在线观看国产| 丝瓜视频免费看黄片| av网站免费在线观看视频| av国产精品久久久久影院| 久久国产精品男人的天堂亚洲 | 在线观看人妻少妇| 美女主播在线视频| 欧美xxxx性猛交bbbb| 女人久久www免费人成看片| 九色成人免费人妻av| 超碰97精品在线观看| 亚洲国产高清在线一区二区三| 成人漫画全彩无遮挡| 久久久久久久久久久丰满| 午夜日本视频在线| 97热精品久久久久久| 亚洲精华国产精华液的使用体验| 777米奇影视久久| 国产一级毛片在线| 午夜福利高清视频| 高清欧美精品videossex| 欧美+日韩+精品| 十分钟在线观看高清视频www | 亚洲国产精品国产精品| 国产一区二区三区av在线| 99久国产av精品国产电影| 色5月婷婷丁香| 色吧在线观看| 欧美极品一区二区三区四区| 日本与韩国留学比较| 一区二区三区乱码不卡18| 亚洲无线观看免费| 国产精品欧美亚洲77777| 国产爽快片一区二区三区| 男的添女的下面高潮视频| 亚洲aⅴ乱码一区二区在线播放| 国产欧美日韩精品一区二区| 九九久久精品国产亚洲av麻豆| 国产中年淑女户外野战色| 亚洲av二区三区四区| 国产精品欧美亚洲77777| 久久精品国产鲁丝片午夜精品| 免费高清在线观看视频在线观看| 在线 av 中文字幕| 一级毛片黄色毛片免费观看视频| 色视频在线一区二区三区| 联通29元200g的流量卡| 最新中文字幕久久久久| 性色av一级| 国产精品一区二区性色av| 国产成人免费无遮挡视频| 久久久欧美国产精品| 国产成人91sexporn| 国产精品国产三级国产专区5o| 日本免费在线观看一区| 久久青草综合色| a级一级毛片免费在线观看| 国内揄拍国产精品人妻在线| 一级av片app| 1000部很黄的大片| 我要看黄色一级片免费的| 日韩,欧美,国产一区二区三区| 男女边吃奶边做爰视频| 欧美97在线视频| 99久久精品热视频| 中文字幕免费在线视频6| 亚洲精品视频女| 亚洲综合精品二区| 老司机影院毛片| 色视频在线一区二区三区| 国产欧美日韩精品一区二区| 日韩三级伦理在线观看| 新久久久久国产一级毛片| 国产午夜精品久久久久久一区二区三区| 一本一本综合久久| 亚洲电影在线观看av| 日日摸夜夜添夜夜爱| kizo精华| 激情五月婷婷亚洲| 午夜免费观看性视频| 美女国产视频在线观看| 老熟女久久久| 男女下面进入的视频免费午夜| 麻豆乱淫一区二区| 高清黄色对白视频在线免费看 | 涩涩av久久男人的天堂| 黄片wwwwww| 久久精品国产亚洲网站| 亚洲精品亚洲一区二区| 国产精品99久久99久久久不卡 | 国产在线免费精品| 丰满迷人的少妇在线观看| 黄色欧美视频在线观看| 国产欧美另类精品又又久久亚洲欧美| 中文在线观看免费www的网站| 观看av在线不卡| 少妇 在线观看| .国产精品久久| 中国国产av一级| 久久国产精品男人的天堂亚洲 | av专区在线播放| 亚洲一级一片aⅴ在线观看| 99热这里只有精品一区| 黄片wwwwww| 久久99热6这里只有精品| 中文天堂在线官网| 国产淫片久久久久久久久| 久久久久久久久久人人人人人人| 99久久综合免费| 欧美人与善性xxx| 精品一区在线观看国产| 夜夜爽夜夜爽视频| 黑人高潮一二区| 亚洲国产精品国产精品| 99久久精品热视频| 久久久久久久久久久丰满| 日韩中字成人| 特大巨黑吊av在线直播| 看非洲黑人一级黄片| 欧美丝袜亚洲另类| 2022亚洲国产成人精品| 国产乱来视频区| 日韩人妻高清精品专区| 国产精品久久久久久精品电影小说 | 视频中文字幕在线观看| 伦精品一区二区三区| 国产视频首页在线观看| av视频免费观看在线观看| 亚洲国产欧美在线一区| 黄色欧美视频在线观看| 国产一区亚洲一区在线观看| 日本欧美视频一区| 国产亚洲欧美精品永久| 一级毛片aaaaaa免费看小| 亚洲自偷自拍三级| 欧美一级a爱片免费观看看| 一级毛片 在线播放| 99热这里只有精品一区| 国产精品秋霞免费鲁丝片| 蜜桃亚洲精品一区二区三区| 精华霜和精华液先用哪个| 色视频www国产| 汤姆久久久久久久影院中文字幕| 麻豆精品久久久久久蜜桃| 久久精品人妻少妇| 久久人人爽人人片av| 久久久亚洲精品成人影院| 精品一区二区三区视频在线| av网站免费在线观看视频| 久久久久网色| 中文天堂在线官网| 97精品久久久久久久久久精品| 91久久精品电影网| 国产精品嫩草影院av在线观看| 亚洲国产高清在线一区二区三| 97热精品久久久久久| 人人妻人人澡人人爽人人夜夜| 久久国产亚洲av麻豆专区| 精品少妇黑人巨大在线播放| 黄色怎么调成土黄色| 日韩在线高清观看一区二区三区| 美女xxoo啪啪120秒动态图| 熟女av电影| 久久影院123| 亚洲欧美成人精品一区二区| 免费大片18禁| 精品少妇黑人巨大在线播放| 啦啦啦视频在线资源免费观看| 黄色配什么色好看| 大片免费播放器 马上看| 三级经典国产精品| 卡戴珊不雅视频在线播放| 少妇高潮的动态图| 毛片一级片免费看久久久久| 夫妻性生交免费视频一级片| 王馨瑶露胸无遮挡在线观看| 插逼视频在线观看| 五月天丁香电影| 欧美xxⅹ黑人| 乱系列少妇在线播放| 国产爱豆传媒在线观看| 久久精品国产亚洲网站| 久久国产乱子免费精品| 午夜精品国产一区二区电影| 亚洲欧美成人综合另类久久久| 蜜桃久久精品国产亚洲av| 视频区图区小说| 老熟女久久久| 日韩三级伦理在线观看| 久久这里有精品视频免费| 欧美另类一区| 国产伦在线观看视频一区| 综合色丁香网| 成人免费观看视频高清| 亚洲va在线va天堂va国产| 精品国产露脸久久av麻豆| 美女cb高潮喷水在线观看| 亚洲欧美精品自产自拍| 日韩 亚洲 欧美在线| 亚洲精品自拍成人| 婷婷色av中文字幕| 久久久久精品久久久久真实原创| 中文字幕人妻熟人妻熟丝袜美| 中国三级夫妇交换| 联通29元200g的流量卡| 久久ye,这里只有精品| 国产伦理片在线播放av一区| 国产av国产精品国产| 成人高潮视频无遮挡免费网站| 久久精品久久精品一区二区三区| 黄色日韩在线| 亚洲四区av| 亚洲精品aⅴ在线观看| av女优亚洲男人天堂| 日韩av不卡免费在线播放| 久久精品国产亚洲网站| 婷婷色综合大香蕉| 新久久久久国产一级毛片| 日韩三级伦理在线观看| 午夜视频国产福利| 青青草视频在线视频观看| 亚洲国产日韩一区二区| 亚洲国产最新在线播放| 亚洲精品aⅴ在线观看| 51国产日韩欧美| 免费人妻精品一区二区三区视频| 在线亚洲精品国产二区图片欧美 | 国产精品福利在线免费观看| 香蕉精品网在线| 亚洲av成人精品一二三区| 老师上课跳d突然被开到最大视频| 我的女老师完整版在线观看| 免费人成在线观看视频色| 男女国产视频网站| 亚洲最大成人中文| 最后的刺客免费高清国语| 天堂中文最新版在线下载| 国产精品一区www在线观看| 汤姆久久久久久久影院中文字幕| 18禁动态无遮挡网站| 菩萨蛮人人尽说江南好唐韦庄| 一个人免费看片子| 丰满乱子伦码专区| 久久久久国产网址| 夜夜骑夜夜射夜夜干| 91久久精品电影网| 国产69精品久久久久777片| 黄片wwwwww| av黄色大香蕉| 色哟哟·www| 91狼人影院| 一区二区三区四区激情视频| 你懂的网址亚洲精品在线观看| 亚洲国产精品成人久久小说| 中文欧美无线码| 熟妇人妻不卡中文字幕| 丰满人妻一区二区三区视频av| 中文字幕久久专区| 超碰av人人做人人爽久久| 国产精品.久久久| 一区二区三区乱码不卡18| 午夜福利网站1000一区二区三区| 久久久久人妻精品一区果冻| 亚洲在久久综合| 99久久综合免费| 国产免费福利视频在线观看| 免费黄色在线免费观看| 欧美精品一区二区免费开放| 亚洲精品一区蜜桃| 麻豆国产97在线/欧美| 亚洲性久久影院| 天天躁日日操中文字幕| 国产女主播在线喷水免费视频网站| 麻豆乱淫一区二区| 国产精品av视频在线免费观看| 成人美女网站在线观看视频| 丰满迷人的少妇在线观看| 性高湖久久久久久久久免费观看| av不卡在线播放| 男女啪啪激烈高潮av片| 国产欧美日韩一区二区三区在线 | 日韩中字成人| 亚洲av成人精品一区久久| av一本久久久久| 国产欧美亚洲国产| 国内精品宾馆在线| av在线观看视频网站免费| 欧美高清性xxxxhd video| 国产男人的电影天堂91| 欧美人与善性xxx| 97热精品久久久久久| 男的添女的下面高潮视频| 婷婷色av中文字幕| 一本色道久久久久久精品综合| 国产成人91sexporn| 青春草视频在线免费观看| 偷拍熟女少妇极品色| av福利片在线观看| 国产黄色视频一区二区在线观看| 3wmmmm亚洲av在线观看| 欧美97在线视频| 免费大片18禁| 全区人妻精品视频| a级毛色黄片| 菩萨蛮人人尽说江南好唐韦庄| 日韩电影二区| 久久97久久精品| 国产免费福利视频在线观看| 欧美3d第一页| 一本—道久久a久久精品蜜桃钙片| 另类亚洲欧美激情| 嫩草影院入口| 男的添女的下面高潮视频| 国产精品久久久久久久久免| 亚洲自偷自拍三级| 久久久久久久久久成人| 国产精品久久久久成人av| 国产爱豆传媒在线观看| 欧美极品一区二区三区四区| 日韩av免费高清视频| 日韩,欧美,国产一区二区三区| 久久国内精品自在自线图片| 国产午夜精品久久久久久一区二区三区| 97在线人人人人妻| 噜噜噜噜噜久久久久久91| 久久女婷五月综合色啪小说| 久久久久国产精品人妻一区二区| 亚洲av成人精品一二三区| 秋霞伦理黄片| 欧美zozozo另类| 寂寞人妻少妇视频99o| 亚洲av综合色区一区| 各种免费的搞黄视频| 国产毛片在线视频| 免费观看的影片在线观看| 黄片无遮挡物在线观看| 在线观看人妻少妇| 免费大片18禁| 日本-黄色视频高清免费观看| 黄色视频在线播放观看不卡| 日韩精品有码人妻一区| 成人午夜精彩视频在线观看| 国产精品久久久久久av不卡| 成人毛片a级毛片在线播放| 国产一级毛片在线| 国产国拍精品亚洲av在线观看| 人妻制服诱惑在线中文字幕| 在线免费十八禁| 免费观看a级毛片全部| 久久国产精品男人的天堂亚洲 | 久久人妻熟女aⅴ| 欧美3d第一页| 爱豆传媒免费全集在线观看| 亚洲久久久国产精品| 亚洲四区av| 哪个播放器可以免费观看大片| 夜夜爽夜夜爽视频| 亚洲aⅴ乱码一区二区在线播放| av网站免费在线观看视频| 免费黄网站久久成人精品| 久久久精品94久久精品| 色网站视频免费| 色视频www国产| 亚洲精品色激情综合| 伊人久久精品亚洲午夜| 国产精品久久久久久精品电影小说 | 亚洲人与动物交配视频| 九九久久精品国产亚洲av麻豆| 国产精品人妻久久久久久| 亚洲va在线va天堂va国产| 亚洲美女黄色视频免费看| 亚洲av成人精品一区久久| 免费观看性生交大片5| 男的添女的下面高潮视频| 国产深夜福利视频在线观看| 久久影院123| 亚洲欧美精品自产自拍| 国产欧美日韩一区二区三区在线 | 久久久久视频综合| 成人高潮视频无遮挡免费网站| 人妻少妇偷人精品九色| 亚洲伊人久久精品综合| 尾随美女入室| 国产老妇伦熟女老妇高清| 国内揄拍国产精品人妻在线| 欧美高清成人免费视频www| 亚洲精品国产av成人精品| 哪个播放器可以免费观看大片| 在线观看免费视频网站a站| 好男人视频免费观看在线| 自拍偷自拍亚洲精品老妇| 日本欧美视频一区| 午夜福利视频精品| 国产亚洲5aaaaa淫片| 观看免费一级毛片| 免费观看性生交大片5| 日韩一本色道免费dvd| 乱码一卡2卡4卡精品| 免费观看性生交大片5| 久久精品人妻少妇| 免费黄频网站在线观看国产| 久久青草综合色| 男女边摸边吃奶| 久久久a久久爽久久v久久| 国产av国产精品国产| 国产黄色免费在线视频| 久久久久久久亚洲中文字幕| 91久久精品国产一区二区成人| 国产极品天堂在线| 黄色欧美视频在线观看| 99久久综合免费| 国产精品无大码| 99久久精品一区二区三区| 亚洲综合精品二区| 少妇被粗大猛烈的视频| 亚洲精品一二三| 中文乱码字字幕精品一区二区三区| 日本wwww免费看| 91狼人影院|