• <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香蕉五月| 成人国产一区最新在线观看| 在线十欧美十亚洲十日本专区| 亚洲中文日韩欧美视频| 国产探花在线观看一区二区| 精品人妻一区二区三区麻豆 | 亚洲国产欧洲综合997久久,| 亚洲 国产 在线| 婷婷亚洲欧美| 成年版毛片免费区| 黄色配什么色好看| 啦啦啦韩国在线观看视频| 51国产日韩欧美| 简卡轻食公司| 久久草成人影院| 久久精品国产99精品国产亚洲性色| 一级作爱视频免费观看| 看黄色毛片网站| 在线播放国产精品三级| 国产单亲对白刺激| 麻豆成人午夜福利视频| 69av精品久久久久久| 精品免费久久久久久久清纯| 成人特级av手机在线观看| 国产精品影院久久| 夜夜看夜夜爽夜夜摸| 黄色日韩在线| 伊人久久精品亚洲午夜| 国产伦一二天堂av在线观看| 好看av亚洲va欧美ⅴa在| 脱女人内裤的视频| 国产大屁股一区二区在线视频| 成熟少妇高潮喷水视频| 最新在线观看一区二区三区| 日本撒尿小便嘘嘘汇集6| 99精品久久久久人妻精品| 一区二区三区免费毛片| 亚洲七黄色美女视频| 欧美日韩乱码在线| 欧美性猛交╳xxx乱大交人| 99国产综合亚洲精品| 舔av片在线| 国产一区二区激情短视频| 亚洲激情在线av| 国产欧美日韩一区二区三| 黄色女人牲交| 少妇人妻一区二区三区视频| 亚洲成人中文字幕在线播放| 国产高清视频在线播放一区| 麻豆久久精品国产亚洲av| 99久久久亚洲精品蜜臀av| 国产乱人视频| 国产乱人视频| 99国产极品粉嫩在线观看| 欧美又色又爽又黄视频| 亚洲精品一区av在线观看| 精品久久久久久,| 18禁黄网站禁片免费观看直播| 成人国产一区最新在线观看| 午夜免费成人在线视频| 日韩欧美国产一区二区入口| 日本五十路高清| 婷婷精品国产亚洲av在线| 男人舔奶头视频| 在线观看av片永久免费下载| 日韩人妻高清精品专区| 亚洲自拍偷在线| 免费看a级黄色片| 精品人妻视频免费看| 美女大奶头视频| 2021天堂中文幕一二区在线观| 亚洲美女黄片视频| 18禁黄网站禁片免费观看直播| 精品久久久久久,| avwww免费| 久久精品夜夜夜夜夜久久蜜豆| avwww免费| 久久精品影院6| 欧美性猛交黑人性爽| 深爱激情五月婷婷| 久久6这里有精品| 久久6这里有精品| 亚洲精品影视一区二区三区av| 久久久久国内视频| 国产精品嫩草影院av在线观看 | 在线观看一区二区三区| 亚洲精品日韩av片在线观看| 午夜a级毛片| 亚洲,欧美精品.| 九色成人免费人妻av| 中文字幕熟女人妻在线| 亚洲狠狠婷婷综合久久图片| 亚洲av日韩精品久久久久久密| 中文字幕熟女人妻在线| 日本a在线网址| 一本久久中文字幕| 国产成人av教育| 午夜激情福利司机影院| 久久久久性生活片| 亚洲avbb在线观看| 亚洲天堂国产精品一区在线| 看免费av毛片| 搡老熟女国产l中国老女人| 内射极品少妇av片p| 亚洲av二区三区四区| 久久国产精品影院| 91久久精品电影网| 深夜a级毛片| 亚洲综合色惰| 国产国拍精品亚洲av在线观看| 国产黄a三级三级三级人| 国产精品自产拍在线观看55亚洲| 午夜福利免费观看在线| 国产伦人伦偷精品视频| 欧美xxxx黑人xx丫x性爽| 国模一区二区三区四区视频| 久久久国产成人免费| av天堂中文字幕网| 99久久九九国产精品国产免费| 成人特级黄色片久久久久久久| 三级国产精品欧美在线观看| 一级a爱片免费观看的视频| 久久精品影院6| 中亚洲国语对白在线视频| 一级黄片播放器| 夜夜躁狠狠躁天天躁| 99久久九九国产精品国产免费| 99国产极品粉嫩在线观看| 天堂动漫精品| 成人无遮挡网站| 欧美3d第一页| 天堂av国产一区二区熟女人妻| 一区二区三区激情视频| 国产主播在线观看一区二区| 国产麻豆成人av免费视频| 两性午夜刺激爽爽歪歪视频在线观看| 日韩欧美三级三区| 蜜桃亚洲精品一区二区三区| 亚洲第一区二区三区不卡| 麻豆国产97在线/欧美| x7x7x7水蜜桃| 99久久九九国产精品国产免费| 日韩欧美三级三区| 国产精品伦人一区二区| 亚洲国产高清在线一区二区三| 有码 亚洲区| 99热这里只有是精品在线观看 | 亚洲五月天丁香| 九色国产91popny在线| 久久久久久九九精品二区国产| 国产精品自产拍在线观看55亚洲| 精品久久久久久久久久久久久| 身体一侧抽搐| 女人十人毛片免费观看3o分钟| av黄色大香蕉| 亚洲成人精品中文字幕电影| 国产精品嫩草影院av在线观看 | 又紧又爽又黄一区二区| 乱人视频在线观看| 欧美极品一区二区三区四区| 国产精品亚洲一级av第二区| a级毛片免费高清观看在线播放| 如何舔出高潮| 大型黄色视频在线免费观看| 性色avwww在线观看| 国产欧美日韩精品一区二区| 国产野战对白在线观看| 男人和女人高潮做爰伦理| 午夜日韩欧美国产| 亚洲av免费在线观看| 国产中年淑女户外野战色| 天天一区二区日本电影三级| 哪里可以看免费的av片| 欧美日韩亚洲国产一区二区在线观看| 国产高清视频在线播放一区| 久久午夜福利片| av专区在线播放| 日韩中字成人| 18禁黄网站禁片午夜丰满| 欧美日韩瑟瑟在线播放| 国产私拍福利视频在线观看| 欧美3d第一页| 日本免费一区二区三区高清不卡| 久久久久国内视频| 欧美激情国产日韩精品一区| 少妇的逼好多水| 两人在一起打扑克的视频| 激情在线观看视频在线高清| 国产亚洲欧美在线一区二区| 欧美国产日韩亚洲一区| 自拍偷自拍亚洲精品老妇| 97热精品久久久久久| 美女被艹到高潮喷水动态| 精华霜和精华液先用哪个| 一边摸一边抽搐一进一小说| 成人永久免费在线观看视频| 亚洲精品色激情综合| 女人被狂操c到高潮| 麻豆成人av在线观看| 亚洲国产欧美人成| 午夜久久久久精精品| 欧美又色又爽又黄视频| 精品欧美国产一区二区三| 好看av亚洲va欧美ⅴa在| 我的女老师完整版在线观看| 亚洲第一欧美日韩一区二区三区| 嫩草影视91久久| 免费人成视频x8x8入口观看| 国产av麻豆久久久久久久| 少妇被粗大猛烈的视频| 午夜激情福利司机影院| 国产精品一区二区性色av| 91在线精品国自产拍蜜月| 夜夜躁狠狠躁天天躁| 搡老岳熟女国产| 麻豆成人午夜福利视频| 99视频精品全部免费 在线| 一个人免费在线观看的高清视频| 男女那种视频在线观看| 国产主播在线观看一区二区| 国产成人欧美在线观看| 我的老师免费观看完整版| 一夜夜www| 国产在线精品亚洲第一网站| 在线观看一区二区三区| 性欧美人与动物交配| 亚洲 欧美 日韩 在线 免费| 国产精品精品国产色婷婷| 真实男女啪啪啪动态图| 少妇人妻一区二区三区视频| 精品无人区乱码1区二区| 免费av不卡在线播放| 欧美一区二区国产精品久久精品| 乱码一卡2卡4卡精品| 久久国产乱子伦精品免费另类| 国产精品嫩草影院av在线观看 | 国产精品一区二区性色av| 国产高清有码在线观看视频| netflix在线观看网站| 国产私拍福利视频在线观看| 亚洲男人的天堂狠狠| 日韩精品青青久久久久久| 国产精品久久视频播放| 好看av亚洲va欧美ⅴa在| 亚洲国产欧洲综合997久久,| 亚洲美女搞黄在线观看 | 亚洲中文字幕一区二区三区有码在线看| 欧美日韩乱码在线| 欧美日本视频| 色哟哟·www| 高清在线国产一区| 亚洲欧美日韩无卡精品| av女优亚洲男人天堂| 99久久九九国产精品国产免费| 午夜福利免费观看在线| 窝窝影院91人妻| 亚洲美女视频黄频| 亚洲av不卡在线观看| 国内少妇人妻偷人精品xxx网站| 欧美成人一区二区免费高清观看| 亚洲va日本ⅴa欧美va伊人久久| 色哟哟·www| 日韩欧美国产在线观看| 成人三级黄色视频| 最新在线观看一区二区三区| 久久久久亚洲av毛片大全| 桃红色精品国产亚洲av| 亚洲av二区三区四区| 日本一本二区三区精品| 久久久成人免费电影| 欧美色欧美亚洲另类二区| 夜夜看夜夜爽夜夜摸| 99国产精品一区二区蜜桃av| 日本熟妇午夜| 美女xxoo啪啪120秒动态图 | 日本 欧美在线| 看黄色毛片网站| 美女黄网站色视频| 国产av一区在线观看免费| 亚洲国产精品999在线| 久久久国产成人免费| 一级黄色大片毛片| 色播亚洲综合网| 国产精品亚洲美女久久久| 99在线人妻在线中文字幕| 欧美高清性xxxxhd video| av天堂在线播放| 国产淫片久久久久久久久 | 一夜夜www| 人人妻人人看人人澡| 日韩欧美一区二区三区在线观看| 夜夜夜夜夜久久久久| 国内少妇人妻偷人精品xxx网站| 高清毛片免费观看视频网站| 激情在线观看视频在线高清| 亚洲av熟女| 9191精品国产免费久久| 亚洲av电影在线进入| 脱女人内裤的视频| 97超视频在线观看视频| av在线老鸭窝| 在线免费观看的www视频| 欧洲精品卡2卡3卡4卡5卡区| 国产一区二区在线观看日韩| 男人狂女人下面高潮的视频| 亚洲国产精品999在线| 国内精品久久久久久久电影| 精品人妻偷拍中文字幕| 亚洲人与动物交配视频| 国产精品综合久久久久久久免费| 免费观看的影片在线观看| 欧美日韩瑟瑟在线播放| 日日干狠狠操夜夜爽| 韩国av一区二区三区四区| 国产精品日韩av在线免费观看| 丰满人妻熟妇乱又伦精品不卡| 深夜a级毛片| 99久久精品国产亚洲精品| 日本免费一区二区三区高清不卡| 久久久国产成人免费| 欧美最新免费一区二区三区 | 精品一区二区三区av网在线观看| 亚州av有码| 日韩欧美国产在线观看| 免费电影在线观看免费观看| 国产高清三级在线| 国产精品一及| 亚洲人成网站高清观看| 波多野结衣高清作品| 国产综合懂色| 国产成年人精品一区二区| 国产综合懂色| 日本精品一区二区三区蜜桃| 国产在线精品亚洲第一网站| 国产高清视频在线观看网站| 脱女人内裤的视频| 国产精品日韩av在线免费观看| 免费观看精品视频网站| 亚洲国产欧洲综合997久久,| 看免费av毛片| 色综合站精品国产| 亚洲中文日韩欧美视频| 成人特级av手机在线观看| 国产av麻豆久久久久久久| 哪里可以看免费的av片| 日日摸夜夜添夜夜添小说| АⅤ资源中文在线天堂| 日本一本二区三区精品| 长腿黑丝高跟| 天堂av国产一区二区熟女人妻| 亚洲内射少妇av| 亚洲在线观看片| 精品一区二区免费观看| 级片在线观看| 亚洲自偷自拍三级| 国产精品电影一区二区三区| a级毛片免费高清观看在线播放| 国产欧美日韩精品一区二区| 美女黄网站色视频| 怎么达到女性高潮| 国产黄a三级三级三级人| 又爽又黄a免费视频| 草草在线视频免费看| 久久精品国产亚洲av天美| 久久这里只有精品中国| 精品人妻一区二区三区麻豆 | 女同久久另类99精品国产91| 一边摸一边抽搐一进一小说| 精品国产亚洲在线| 日韩亚洲欧美综合| 成人国产一区最新在线观看| 日本黄色片子视频| 亚洲精品日韩av片在线观看| 亚洲久久久久久中文字幕| 国产不卡一卡二| 国产精品久久电影中文字幕| 久久精品国产亚洲av香蕉五月| 一卡2卡三卡四卡精品乱码亚洲| 一a级毛片在线观看| 精品久久久久久成人av| 久久婷婷人人爽人人干人人爱| 色尼玛亚洲综合影院| 中文字幕人成人乱码亚洲影| 亚洲第一电影网av| 日本精品一区二区三区蜜桃| 天堂av国产一区二区熟女人妻| 九九热线精品视视频播放| 亚洲午夜理论影院| 亚洲七黄色美女视频| 黄色一级大片看看| 亚洲精品影视一区二区三区av| 久久国产精品影院| 欧美乱妇无乱码| 国产免费一级a男人的天堂| 成人性生交大片免费视频hd| 亚洲真实伦在线观看| 性欧美人与动物交配| av国产免费在线观看| 一卡2卡三卡四卡精品乱码亚洲| 热99在线观看视频| 日本一本二区三区精品| 久久久久久国产a免费观看| 免费av观看视频| 免费看美女性在线毛片视频| 老熟妇乱子伦视频在线观看| 久久久久九九精品影院| 每晚都被弄得嗷嗷叫到高潮| 亚洲人成伊人成综合网2020| 亚洲成人中文字幕在线播放| 亚洲美女搞黄在线观看 | 少妇被粗大猛烈的视频| 99热这里只有精品一区| 欧美黄色片欧美黄色片| 精品久久久久久久人妻蜜臀av| 听说在线观看完整版免费高清| 熟妇人妻久久中文字幕3abv| 亚洲五月天丁香| 日韩中文字幕欧美一区二区| 亚洲国产色片| 午夜福利在线在线| 午夜福利成人在线免费观看| 久久中文看片网| 99热这里只有是精品50| 90打野战视频偷拍视频| 在线十欧美十亚洲十日本专区| 日本 欧美在线| 亚洲电影在线观看av| 别揉我奶头~嗯~啊~动态视频| 日本免费一区二区三区高清不卡| 成人鲁丝片一二三区免费| 18禁黄网站禁片免费观看直播| 亚洲不卡免费看| 97超视频在线观看视频| 亚洲中文日韩欧美视频| 日日干狠狠操夜夜爽| 亚洲专区国产一区二区| 麻豆成人午夜福利视频| 色在线成人网| 精品一区二区免费观看| 国产aⅴ精品一区二区三区波| 天堂√8在线中文| 看黄色毛片网站| 国产亚洲精品久久久久久毛片| 我要搜黄色片| 欧美最黄视频在线播放免费| 香蕉av资源在线| 成熟少妇高潮喷水视频| 国产男靠女视频免费网站| 成人av在线播放网站| 亚洲 欧美 日韩 在线 免费| 最后的刺客免费高清国语| 婷婷精品国产亚洲av| 嫁个100分男人电影在线观看| 国产高清视频在线播放一区| 久久久国产成人精品二区| 国产伦精品一区二区三区视频9| 国产成人a区在线观看| 九九热线精品视视频播放| 欧美绝顶高潮抽搐喷水| 18美女黄网站色大片免费观看| 国产免费一级a男人的天堂| 在线看三级毛片| 在线观看一区二区三区| 国内精品久久久久精免费| 久久久精品欧美日韩精品| 成人高潮视频无遮挡免费网站| www.www免费av| 日韩免费av在线播放| 久久久国产成人精品二区| 精品一区二区免费观看| 免费一级毛片在线播放高清视频| 亚洲七黄色美女视频| 中文亚洲av片在线观看爽| 亚洲成人久久爱视频| 亚洲国产精品999在线| 精品人妻一区二区三区麻豆 | 国产三级中文精品| 欧美成人免费av一区二区三区| 美女黄网站色视频| 在线播放无遮挡| 嫩草影院入口| 悠悠久久av| 性欧美人与动物交配| avwww免费| 日韩大尺度精品在线看网址| 国产精品自产拍在线观看55亚洲| 亚洲精品一卡2卡三卡4卡5卡| 久久久色成人| 成人三级黄色视频| 99热这里只有是精品50| 免费大片18禁| 久久国产精品影院| 中文字幕精品亚洲无线码一区| 欧美bdsm另类| 天堂网av新在线| 日韩欧美免费精品| 在线观看免费视频日本深夜| x7x7x7水蜜桃| www.色视频.com| 真人做人爱边吃奶动态| 91久久精品国产一区二区成人| 国内精品美女久久久久久| 美女高潮的动态| 宅男免费午夜| 91狼人影院| 午夜福利在线观看吧| 大型黄色视频在线免费观看| 免费黄网站久久成人精品 | 欧美精品啪啪一区二区三区| 一本一本综合久久| 美女免费视频网站| 亚洲人成网站高清观看| 少妇裸体淫交视频免费看高清| 国产亚洲精品久久久久久毛片| 99热6这里只有精品| 国产高清激情床上av| 免费在线观看日本一区| 成人午夜高清在线视频| 一进一出抽搐gif免费好疼| 成人鲁丝片一二三区免费| 黄片小视频在线播放| 日韩欧美在线乱码| 亚洲精华国产精华精| 91字幕亚洲| 亚洲不卡免费看| 欧美日韩中文字幕国产精品一区二区三区| 欧美日韩瑟瑟在线播放| 国产真实乱freesex| 3wmmmm亚洲av在线观看| 窝窝影院91人妻| 色播亚洲综合网| 国产乱人视频| 直男gayav资源| or卡值多少钱| 熟女电影av网| 男人狂女人下面高潮的视频| 十八禁网站免费在线| 在线十欧美十亚洲十日本专区| 欧美色视频一区免费| 3wmmmm亚洲av在线观看| 日日摸夜夜添夜夜添小说| 免费人成视频x8x8入口观看| 又粗又爽又猛毛片免费看| 色尼玛亚洲综合影院| 综合色av麻豆| 熟女电影av网| 久久精品国产99精品国产亚洲性色| 日本精品一区二区三区蜜桃| 成人高潮视频无遮挡免费网站| 久9热在线精品视频| 久久精品久久久久久噜噜老黄 | 国产精品影院久久| 乱码一卡2卡4卡精品| 亚洲av一区综合| 色哟哟哟哟哟哟| 特大巨黑吊av在线直播| 色哟哟·www| 免费av不卡在线播放| 日本撒尿小便嘘嘘汇集6| 美女黄网站色视频| 日本熟妇午夜| 国产美女午夜福利| 国产精品久久久久久精品电影| 色综合婷婷激情| 亚洲欧美清纯卡通| 一级a爱片免费观看的视频| 国产精品电影一区二区三区| 亚洲精品影视一区二区三区av| av天堂在线播放| 精品久久久久久久久久免费视频| 一区二区三区四区激情视频 | 亚洲avbb在线观看| 淫秽高清视频在线观看| 亚洲 欧美 日韩 在线 免费| 特大巨黑吊av在线直播| 欧美性猛交黑人性爽| 国产精品久久久久久久电影| 亚洲精品456在线播放app | 国内精品久久久久久久电影| 日韩欧美三级三区| 欧美日韩瑟瑟在线播放| 在线免费观看不下载黄p国产 | 久久久色成人| 热99在线观看视频| 国产熟女xx| 最近在线观看免费完整版| 国产高清激情床上av| 日日摸夜夜添夜夜添av毛片 | 中文字幕久久专区| 欧美黑人巨大hd| 亚洲精品一区av在线观看| 午夜精品一区二区三区免费看| 一个人观看的视频www高清免费观看| 亚洲精品在线观看二区| 欧美一级a爱片免费观看看| 欧美激情在线99| 精品人妻1区二区| 十八禁国产超污无遮挡网站| 久久久久久久精品吃奶| 午夜免费成人在线视频| 啦啦啦观看免费观看视频高清| 日本 欧美在线| 精品人妻1区二区| 中文字幕久久专区| 久久午夜福利片| 搡老岳熟女国产| 在线国产一区二区在线| 一本一本综合久久| 黄色女人牲交| 国产探花极品一区二区| bbb黄色大片| 欧美黑人欧美精品刺激|