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

    紫坪鋪水庫區(qū)域地殼QS動態(tài)變化及其與水庫蓄水關系的研究*

    2012-09-15 08:14:52王惠琳張曉東周龍泉盧顯徐曉楓
    地震學報 2012年5期
    關鍵詞:紫坪小震蓄水

    王惠琳張曉東周龍泉盧 顯徐曉楓

    1)中國???70203海南省地震局

    2)中國北京100045中國地震臺網(wǎng)中心

    紫坪鋪水庫區(qū)域地殼QS動態(tài)變化及其與水庫蓄水關系的研究*

    王惠琳1),張曉東2)周龍泉2)盧 顯2)徐曉楓1)

    1)中國???70203海南省地震局

    2)中國北京100045中國地震臺網(wǎng)中心

    利用2004年8月16日—2008年5月11日雙差精定位后的地震S波資料,采用衰減層析成像方法對紫坪鋪水庫區(qū)域地殼QS進行動態(tài)衰減成像,反演得到了該區(qū)的QS動態(tài)演化特征.以1個月為步長、1年為窗長,計算了紫坪鋪水庫地區(qū)不同大小區(qū)域按月滑動的平均QS值,統(tǒng)計了不同震級范圍按月滑動的地震頻次,并分別與水位進行了對比分析.研究結果表明,水位首次快速加載后,紫坪鋪水庫低QS值區(qū)域范圍擴大,QS值下降,水庫西南區(qū)、東北區(qū)、庫區(qū)中段、東北岸QS值的變化較為明顯.紫坪鋪水庫不同大小區(qū)域按月滑動的平均QS值隨蓄水時間的增加、水位的升高而減小,庫區(qū)中段QS值下降最為劇烈.水庫蓄水增強了該區(qū)的地震活動性,對ML0.0—2.0的小震活動影響較為顯著.

    S波衰減因子二維成像 QS值動態(tài)變化 地震活動性 紫坪鋪水庫

    Abstract:Using a two-dimensional seismic attenuation tomography technique,the crustal QSdynamic variation in Zipingpu reservoir area,Sichuan Province,was obtained by inverting S-wave data from relocated earthquakes from 16 Aug.2004 to 11 May 2008 with the double-difference algorithm.Monthly moving average QSvalue in different-sized regions and earthquake frequency for different magnitude range were compared with the water level of the Zipingpu reservoir,respectively(the window length is 1 year and the step length is 1 month).Theresult shows that low-QSregion expanded and QSvalue dropped significantly around the reservoir area during the first fast water-loading period.In particular,QSvalue varied remarkably in the southwestern,northeastern and middle parts of the reservoir area.Monthly moving average QSvalues for different-sized regions in the Zipingpu reservoir area decrease with elapsed impoundment time and increased water level,and the QSvalue in the middle reservoir area dropped sharply.Furthermore,the impoundment enhanced seismic activity in the reservoir area,and mainly increased activity of small earthquakes of ML0.0—2.0.

    Key words:two-dimensional tomography of S-wave attenuation;QSdynamic variation;seismic activity;Zipingpu reservoir

    引言

    在地殼介質(zhì)中傳播的地震波振幅的衰減,與構造斷裂帶的破碎程度、活動強度,以及巖漿、水、石油等流體活動密切相關(汪素云等,2008).實驗研究表明,巖石中流體的存在將會影響巖石介質(zhì)的地震波傳播特性.陳颙等(2009)研究發(fā)現(xiàn),地震波衰減主要取決于巖石的微觀性質(zhì)諸如巖石內(nèi)部裂紋的密度、分布、構造以及孔隙流體的相互作用.王大興等(2006)認為對于橫波,巖石孔隙流體黏滯運動引起的能量耗散是造成能量衰減的主要因素.劉建華等(2004)通過對地震波衰減的物理機制研究認為,造成地殼介質(zhì)中地震波衰減的主要原因是地殼內(nèi)存在大量裂隙,裂紋中飽含水或部分含水,地震波傳播時引起裂隙中的流體運動,從而造成地震波的衰減.在特定的地質(zhì)環(huán)境下,例如水庫蓄水,水體增加庫區(qū)范圍地殼的荷載,庫水沿著巖石破碎帶、節(jié)理和裂隙向深部滲透,增加了地下巖體的孔隙水壓力,降低了發(fā)震斷層的抗剪強度,使得天然構造地震已積蓄的能量提前釋放(郭永剛等,2008).由于巖石中的裂紋、流體含量等均與地震波的衰減密切相關,因此,通過研究反映地震波衰減的Q值(王惠琳等,2012)及其動態(tài)變化,可進一步探索水庫蓄水區(qū)域庫水的滲透作用及其對淺層地殼介質(zhì)的影響.

    本研究利用2004年8月16日—2008年5月11日的雙差精定位后的ML0.0—4.0地震S波形資料,對紫坪鋪水庫區(qū)域淺層地殼(1—15km)QS進行動態(tài)衰減成像,反演得到了該區(qū)動態(tài)的QS分布特征.以1個月為步長、1年為窗長,計算了紫坪鋪水庫地區(qū)不同大小區(qū)域按月滑動的平均QS值,統(tǒng)計了不同震級范圍按月滑動的地震頻次,并將二者分別與水位進行了對比分析,進而討論了紫坪鋪水庫地區(qū)S波衰減及地震活動性與水庫蓄水的關系.

    1 紫坪鋪水庫及臺網(wǎng)概況

    紫坪鋪水庫是岷江上的一座大型水利樞紐,位于四川省都江堰市區(qū)西北方向9km處的麻溪鄉(xiāng),距離成都市60km.紫坪鋪水庫為大(一)型工程,最大壩高156 m,總庫容11.12×108m3,是一座以灌溉和供水為主,兼有發(fā)電、防洪、環(huán)境保護、旅游等綜合效益的大型水利樞紐工程.2001年3月29日正式動工興建,2005年9月30日下閘蓄水,2006年12月竣工.紫坪鋪水庫位于青藏高原東緣的龍門山造山帶的中南段,橫跨龍門山造山帶自西向東發(fā)育的汶川—茂縣(WMF)、北川—映秀(BYF)、安縣—灌縣(AGF)和廣元—大邑(GDF)4條主干斷裂控制了紫坪鋪水庫地區(qū)的基本構造格架(圖1).其中,位于北川—映秀與安縣—灌縣斷裂之間的兩條次級斷裂直接穿越庫區(qū),它們是通濟場斷裂在平面上的兩條分支,在深部最終匯聚到通濟場斷裂的主滑脫面上(周斌等,2010).

    圖1 紫坪鋪水庫臺網(wǎng)及斷裂分布(王惠琳等,2012)LYS:靈隱寺臺;ZDZ:鉆洞子臺;BAY:白巖臺;BAJ:八角臺;TZP:桃子坪臺;GHS:桂花樹臺;MZP:廟子坪臺;YZP:油榨坪臺;WMF:汶川—茂縣斷裂;BYF:北川—映秀斷裂;TJCF:通濟場斷裂;AGF:安縣—灌縣斷裂;GDF:廣元—大邑斷裂Fig.1 Faults and seismic stations in Zipingpu reservoir area LYS:Lingyinsi station;ZDZ:Zuandongzi station;BAY:Baiyan station;BAJ:Bajiao station;TZP:Taoziping station;GHS:Guihuashu station;MZP:Miaoziping station;YZP:Youzhaping station;WMF:Wenchuan--Maoxian fault;BYF:Beichuan--Yingxiu fault;TJCF:Tongjichang fault;AGF:Anxian--Guanxian fault;GDF:Guangyuan--Dayi fault(Wang et al,2012)

    紫坪鋪水庫數(shù)字遙測地震臺網(wǎng)由靈隱寺(LYS)、鉆洞子(ZDZ)、白巖(BAY)、八角(BAJ)、桃子坪(TZP)、桂花樹(GHS)、廟子坪(MZP)和成都區(qū)域數(shù)字臺網(wǎng)中的油榨坪(YZP)8個數(shù)字遙測地震臺組成(圖1).該臺網(wǎng)于2004年8月試運行,2005年6月通過驗收.臺站大致均勻地分布在庫區(qū)周邊,平均臺距為10km左右.臺站使用JC-V104型短周期地震計,頻帶寬度為1—40Hz,配備EDAS-24L型24位數(shù)據(jù)采集器,數(shù)據(jù)采用超短波段的單工數(shù)字無線傳輸與數(shù)字數(shù)據(jù)網(wǎng)(Digital Data Network)接力傳輸、中心統(tǒng)一記錄的方式.觀測系統(tǒng)動態(tài)范圍優(yōu)于100 dB,采樣率每秒100點.采用了GPS統(tǒng)一授時和定位系統(tǒng),數(shù)據(jù)采集的時間服務精度優(yōu)于1 ms(張永久,趙翠萍,2009;王惠琳等,2012).

    2 方法原理

    地震儀記錄到的數(shù)字地震波形是一種綜合信息,包含了地震震源、地震波的傳播路徑及臺站的場地響應特征等信息.地震波傳播路徑的衰減包括幾何擴散和介質(zhì)的非彈性衰減(用Q值度量).在第j個臺站上觀測到的第i個事件的振幅譜可以寫成如下形式(Scherbaum,1990;Haberland,Rietbrock,2001):

    式中,f為頻率,Si(f)為震源譜,Ij(f)為儀器響應,Rj(f)為場地響應,Gij(r)為沿路徑r的幾何擴散因子,Bij(f)為衰減譜.

    震源譜采用Brune類型的ω2震源模式(Brune,1970;Hanks,Wyss,1972)表示:

    式中,Ω0為與頻率無關的振幅,包括地震矩、幾何擴散等靜態(tài)因素;f為給定頻率;fc為拐角頻率.

    由于本文采用的是直達S波形數(shù)據(jù),最大震中距為49km,因此采用Atkinson和Mereu(1992)提出的直達波幾何衰減函數(shù),即

    式中,r為震源距,系數(shù)b1與頻率無關.對于沿整個射線路徑的衰減譜表示為(Scherbaum,1990;周龍泉等,2009)

    式中,i為第i次地震,j為第j個臺站,tij為沿射線路徑的射線走時,Qij為品質(zhì)因子,t*>ij為沿整個射線路徑的衰減算子.

    由于僅模擬有限頻段內(nèi)的觀測振幅譜,不考慮絕對振幅大小,儀器響應可設為1(Rietbrock,2001).此外,由于大多數(shù)區(qū)域臺站都是建在基巖上,其場地響應可以假定為接近1的常數(shù).因此,方程(1)可以寫成如下形式:

    另外,t*>ij算子可以寫成1/(Q(s)·v(s))沿震源i到臺站j的射線路徑的投影(Cormier,1982;Wittlinger et al,1983),即

    式中,v為S波速度,d s為射線路徑單元.

    如式(5)所示,方程中存在Ω0,fc以及t*>ij三個未知變量.利用式(5)對單條觀測譜進行擬合,3個未知變量取不同的組合時都可能得到最佳擬合解,即解的非唯一性.由于Ω0和fc只與震源有關,因此對同一震源采用多臺觀測譜進行擬合,可以減少未知變量Ω數(shù)量.為限制解的非唯一性,采用Levenberg-Marqurdt或迭代阻尼最小二乘(iterative damped least-squares)法進行譜反演(Hansen et al,2004).通過對多臺記錄的波形聯(lián)合進行譜反演,可以確定同一震源到各個臺站的路徑衰減t*>.假定Q值不變,則t*>與震源至臺站的走時成正比,因此震源深度對t*>數(shù)據(jù)具有較大的影響.利用t*>反演Q值的過程中,必須消除震源深度的影響.為此在計算理論走時時,震中距用震源至臺站的距離來代替(周龍泉等,2009);然后,根據(jù)式(6)采用衰減層析成像方法確定Q結構.

    3 數(shù)據(jù)預處理

    地震波衰減層析成像需要截取S波時間窗進行譜反演,獲得震源到臺站的t*>算子.“S窗”定義為從S波開始到包括S波總能量的90%的時間段.本文“S窗”只包含了直達S波.對于截取后的“S窗”,采用平移窗譜法(劉杰等,2003;蘇有錦等,2006)獲得穩(wěn)定的觀測信號的傅里葉譜.按“每次地震至少有3個以上臺記錄到,每個臺至少有3條記錄”的原則,挑選波形較好,能經(jīng)過信噪比檢驗(SNR>2)的地震記錄.采用多臺記錄譜聯(lián)合反演方法獲得多個臺站的t*>衰減算子.假定一個地震被N個臺站記錄,則進行多臺觀測譜聯(lián)合反演時共需反演N+2個未知變量,即Ω0,fc和t*>1,t*>2,…,t*>N,這比單臺觀測譜反演少反演2 N-2個變量,從而可以降低解的非唯一性.假定共有M個頻率點,則第i個臺站的觀測譜在第j個頻率點的值可以表示為(fj),理論振幅譜為(fj).多臺觀測譜反演就是尋找使殘差最小時的未知變量值.為了尋找全局最優(yōu)解,采用遺傳算法反演N+2個未知量,使多臺觀測譜與理論譜之間在不同頻率點的殘差之和為最小.消除噪聲和幾何擴散影響后,不同臺站位移譜的高頻衰減主要受不同路徑的Q值影響(圖2).本研究分析了1—15Hz的S波譜,且有10Hz以上連續(xù)頻帶寬.圖2為2007年2月9日記錄到的ML2.8地震經(jīng)噪聲校正和幾何校正后在GHS,LYS,MZP,ZDZ,TZP和BAJ臺站的觀測位移譜和擬合譜.

    圖2 2007年2月9日ML2.8地震經(jīng)噪聲校正和幾何校正后在不同臺站的觀測位移譜和擬合譜Fig.2 Observed and fitted displacement spectra of the 9 Feb.2007 ML2.8 earthquake recorded at different stations after noise and geometric spreading were corrected

    4 紫坪鋪水庫地區(qū)動態(tài)QS成像

    4.1 動態(tài)時段劃分及分辨率測試

    為了分析紫坪鋪水庫區(qū)域S波衰減特征隨水庫蓄、放水過程的演化規(guī)律,依據(jù)2004年1月1日—2008年5月11日的壩前水位資料,將庫水加、卸載過程劃分為4個不同的時段進行QS成像研究.如圖3所示,4個時段分別為:a時段(2004-08-16—2005-09-30)為水庫下閘蓄水前,水位在755—760 m之間;b時段(2005-10-01—2006-10-14)為水庫首次快速加載期,水位從776 m上升至第一次蓄水高峰875 m(正常蓄水位為877 m);c時段(2006-04-19—2007-05-21)為水庫第一次高幅加載、卸載期,水位從820 m升至第一次蓄水高峰875 m后快速卸載,水位下降至819 m的低水位;d時段(2007-05-22—2008-05-11)為第二次高幅加載、卸載期,水位從819 m升至第二次蓄水高峰873 m后快速卸載,水位下降至826 m的較低水位.

    動態(tài)QS反演初始輸入模型采用平均vS和平均Q0,即ˉvS=3.2km/s,ˉQ0=90,4個時段均進行了初始模型的歸算統(tǒng)一.表1為各時段經(jīng)過10次迭代反演后的參數(shù),t*>均方根殘差迭代后比迭代前均有所下降.本研究將紫坪鋪水庫地區(qū)地殼在平面上劃分為0.05°× 0.05°大小的均勻網(wǎng)格進行QS成像.圖4、圖5分別為上述4個時段的地震射線分布和檢測板分辨率測試結果.本文僅對反演效果較好的區(qū)域(圖5中虛線框內(nèi)區(qū)域)進行QS成像.其中較大的矩形框為水庫區(qū)域(30.78°—31.16°N,103.30°—103.72°E,以下簡稱研究區(qū)),較小矩形框為庫區(qū)(30.972°—31.062°N,103.444°—103.583°E).從圖5中各時段的檢測板分辨率試驗可見,a,b和d時段研究區(qū)分辨效果較好;c時段地震射線覆蓋不均勻(圖4c),東南區(qū)覆蓋較少,分辨率相對要差;而庫區(qū)各時段地震射線覆蓋密集,分辨效果均較好.

    表1 迭代反演后的參數(shù)Table 1 Parameters after iterative inversion

    圖5 4個時段檢測板分辨率試驗結果(圖(a)對應a時段,以此類推)BYF:北川—映秀斷裂;TJCF:通濟場斷裂;AGF:安縣—灌縣斷裂;GDF:廣元—大邑斷裂Fig.5 Checkerboard resolution test results for the 4 time intervals a,b,c and d,as shown in Fig.3 BYF:Beichuan--Yingxiu fault;TJCF:Tongjichang fault;AGF:Anxian--Guanxian fault;GDF:Guangyuan--Dayi fault

    4.2 動態(tài)QS成像結果及分析

    圖6a--d分別對應研究區(qū)a,b,c和d四個蓄水時段的QS成像結果.水庫下閘蓄水前,即a時段(圖6a),研究區(qū)QS總體較高,庫區(qū)周圍存在東北、東南、西南、西北4個低QS區(qū)域(圖中虛線所示).其中,東北低QS區(qū)范圍較大,數(shù)值較低,沿斷裂帶附近區(qū)域小震活動較為活躍,成叢分布在離庫岸1—10km范圍內(nèi).水庫下閘蓄水后,水位首次快速加載期間(圖6b),東北低QS區(qū)范圍擴大,數(shù)值降低,離庫岸1—7km范圍內(nèi)小震活動有所抑制,僅在7—10km范圍出現(xiàn)小震叢集現(xiàn)象;西南低QS區(qū)范圍有所擴大,且出現(xiàn)了明顯的小震叢集現(xiàn)象;東南、西北低QS區(qū)數(shù)值有所升高.第一次水位高幅加載、卸載期間(圖6c),東北、西南低QS區(qū)范圍繼續(xù)擴大,數(shù)值大幅下降,水庫周邊形成大范圍的低QS值區(qū)域.水庫西南區(qū)小震叢集現(xiàn)象顯著,東北區(qū)小震活動減弱,東南區(qū)出現(xiàn)少量地震活動.第二次水位高幅加載、卸載期間(圖6d),東北低QS區(qū)數(shù)值升高,而西南低QS區(qū)范圍繼續(xù)擴大,數(shù)值大幅下降,區(qū)域地震活動整體減弱,只有東南區(qū)出現(xiàn)小震密集現(xiàn)象.

    圖6 研究區(qū)動態(tài)QS成像結果(圖(a)對應a時段,以此類推)三角為臺站,五角星為汶川主震,空心圓為地震.BYF:北川—映秀斷裂;TJCF:通濟場斷裂;AGF:安縣—灌縣斷裂;GDF:廣元—大邑斷裂Fig.6 Imaged dynamic QSvariation in study area Triangles represent stations,pentagram denotes epicenter of Wenchuan earthquake,open circles indicate earthquakes.BYF:Beichuan--Yingxiu fault;TJCF:Tongjichang fault;AGF:Anxian--Guanxian fault;GDF:Guangyuan--Dayi fault

    圖7a--d分別對應庫區(qū)a,b,c,d四個蓄水時段放大后的QS成像結果.水庫下閘蓄水前,即a時段(圖7a),庫區(qū)整體QS值較高.水庫下閘蓄水后,水庫水位首次快速加載期間(圖7b),庫區(qū)中段區(qū)域QS值開始下降.第一次高幅加載、卸載期間(圖7c),庫區(qū)中段區(qū)域QS保持低值狀態(tài),且低QS區(qū)沿著庫體向兩端延伸.另外,水庫東北岸QS值大幅下降,出現(xiàn)明顯的低QS區(qū)域.第二次高幅加載、卸載期間(圖7d),庫區(qū)中段區(qū)域QS值升高,水庫東北岸低QS區(qū)域范圍縮小,數(shù)值有所回升.庫區(qū)4個時段地震活動均較為平靜,幾乎為空白.

    5 平均QS、地震頻次與水位的關系

    5.1 QS值與水位的關系

    根據(jù)地震精定位后震中分布的特點,將紫坪鋪水庫地區(qū)從大到小依次細分為A區(qū)、B區(qū)、C區(qū)、庫區(qū)中段4個區(qū)域(圖8),以1個月為步長、1年為窗長,空間域上插值求平均分別計算了4個區(qū)域按月滑動的QS值,并與按月滑動的水位進行了對比(圖9).4個區(qū)域經(jīng)緯度范圍分別為A區(qū):30.7°—31.2°N,103.2°—103.8°E;B區(qū):30.91°—31.15°N,103.28°—103.69°E;C區(qū):30.972°—31.062°N,103.444°—103.583°E;庫區(qū)中段(用矩形經(jīng)緯度表示):30.98°—31.02°N,103.48°—103.52°E.

    如圖9所示,橫坐標為蓄水時間,-2表示蓄水前2個月(時間窗為2004-08-01—2005-07-31),0表示蓄水0個月(時間窗為2004-10-01—2005-09-30),2表示蓄水后2個月(時間窗為2004-12-01—2005-11-30),以此類推.從圖9可見,隨著蓄水時間的增加,水位逐漸升高,而A區(qū)、B區(qū)、C區(qū)和庫區(qū)中段按月滑動的平均QS值則逐漸降低,表明地震S波的衰減逐漸增大.當蓄水15個月時,水位達到850 m左右后基本保持不變,各分區(qū)的平均QS值依然繼續(xù)下降.其中,庫區(qū)中段區(qū)域平均QS值從蓄水2個月后隨水位的增加而大幅降低,其受蓄水影響最為顯著.

    圖9 滑動平均QS值、水位隨蓄水時間的變化曲線Fig.9 Variation of moving average QSfor different area scopes(color curves)and water level change(black line)

    根據(jù)李亞林等(2001)、施行覺等(1995)及王大興等(2006)實驗研究發(fā)現(xiàn),S波Q值隨著滲透率的增加而減小,隨含水飽和度的增加剛開始迅速下降,當飽和度達到某個值時,基本不變或略有上升.因此,可能是隨著紫坪鋪水庫蓄水時間的增加,水位逐漸升高.當庫水逐漸沿著巖石破碎帶、節(jié)理、裂隙發(fā)育地區(qū)向地下滲透,使巖石孔隙中充滿流體,內(nèi)摩擦增大,從而使地震波大大衰減,導致QS值大幅下降.庫區(qū)中段區(qū)域被通濟場次級斷裂穿過,且位于庫區(qū)中部,是庫水滲透較為有利的場所,故庫區(qū)中段區(qū)域的地震波衰減對庫水的滲透作用反映較快也較為明顯.

    5.2 地震頻次與水位的關系

    為了解水庫蓄水對紫坪鋪水庫地區(qū)地震活動性的影響,研究了水庫區(qū)域(圖8中B區(qū))不同震級范圍的地震頻次與水位隨蓄水時間的變化關系(圖10).

    如圖10所示,橫坐標為蓄水時間,時間窗選取方法同圖9.隨著蓄水水位的逐漸升高,蓄水4個月時,ML0.0—4.0的小震開始迅速增多;蓄水13—15個月、水位接近840 m時小震最為活躍,地震頻次達到峰值;蓄水15個月后、水位基本保持在850 m不變時,地震頻次則開始減少.其中,ML0.0—1.0的地震數(shù)量占總數(shù)的83%,因此其頻次變化趨勢與ML0.0—4.0地震的變化趨勢基本一致.ML1.0—2.0的地震頻次蓄水5個月時開始增加,蓄水15個月后呈緩慢下降的趨勢;ML2.0—3.0,ML3.0—4.0的地震頻次隨著水位的升高無明顯變化.這表明,紫坪鋪水庫蓄水增強了該區(qū)的地震活動性,且主要對ML0.0—2.0的小震活動影響較為明顯.

    圖10 滑動地震頻次和水位隨蓄水時間的變化曲線Fig.10 Moving averaged earthquake numbers(color curves)and water level change(black curve)

    6 討論與結論

    1)QS動態(tài)成像結果顯示,水庫下閘蓄水前(a時段),研究區(qū)QS總體較高,庫區(qū)周圍存在東北、東南、西南、西北4個低QS區(qū)域.水庫下閘蓄水后,水位首次快速加載期間(b時段),東北低QS區(qū)范圍擴大,數(shù)值降低;西南低QS區(qū)范圍有所擴大.第一次和第二次高幅加載、卸載期間(c,d時段),東北、西南低QS區(qū)范圍繼續(xù)擴大,數(shù)值變化幅度增大,水庫周邊形成大范圍的低QS值區(qū)域.由于c和d時段研究區(qū)地震射線分布相對較少,根據(jù)檢測板分辨試驗結果,c和d時段研究區(qū)的QS變化誤差可能相對要大.庫區(qū)位于研究區(qū)中部,4個時段均有大量地震射線覆蓋,檢測板分辨試驗結果較好.其QS動態(tài)成像結果顯示,蓄水前(a時段),庫區(qū)基本為高QS區(qū)域;下閘蓄水后,水位首次快速加載期間(b時段),庫區(qū)中段區(qū)域QS開始下降;第一次高幅加載、卸載期間(c時段),水庫東北岸出現(xiàn)明顯的低QS區(qū);第二次高幅加、卸載期間(d時段),庫區(qū)中段區(qū)域和東北岸QS值有所回升.此外,按月滑動的平均QS值與水位的關系顯示,紫坪鋪水庫地區(qū)的QS值隨著蓄水時間的增加、水位的升高而減小.其中,庫區(qū)中段區(qū)域的平均QS值下降最為劇烈.上述QS動態(tài)研究結果表明,紫坪鋪水庫區(qū)域S波衰減總體上隨著蓄水而逐漸增大.

    根據(jù)李亞林等(2001)、施行覺等(1995)、王大興等(2006)實驗研究發(fā)現(xiàn),S波Q值隨著滲透率的增加而減??;隨含水飽和度的增加一開始迅速下降,當飽和度達到一定值時,變化趨緩或略有上升.Walsh(1966)和Johnston等(1979)在巖石衰減實驗研究基礎上提出,地震波衰減的重要原因之一是由于巖石中微裂縫的存在及裂隙中存在的流體.在一定的應力波作用下,裂紋表面和晶粒邊界產(chǎn)生滑動摩擦使一部分應力波能量由于克服摩擦產(chǎn)生滑動而轉(zhuǎn)化為熱能耗散掉.施行覺等(1995)認為,孔隙介質(zhì)中的流體引起波的衰減,其機理之一是流體的黏滯運動引起能量耗散(Mavko,Nur,1979);衰減的另一機理,可能是波傳播時引起液體與固體的接觸表面能的變化.

    本研究結果與上述實驗研究成果相符.其中,紫坪鋪水庫西南區(qū)、東北區(qū)、東北岸位于北川—映秀斷裂、通濟場斷裂、安縣—灌縣斷裂之間,節(jié)理、裂隙發(fā)育,巖石較為破碎.隨著紫坪鋪水庫蓄水時間的增加,水位逐漸升高.當庫水逐漸沿著巖石破碎帶、節(jié)理、裂隙發(fā)育地區(qū)向地下滲透,使巖石孔隙中充滿流體,內(nèi)摩擦增大,使地震波大大衰減,導致QS值下降.庫區(qū)中段區(qū)域被通濟場次級斷裂穿過,且位于庫區(qū)中部,是庫水滲透較為有利的場所(王惠琳等,2012),故庫區(qū)中段區(qū)域的地震波衰減對庫水的滲透作用反映較快也較為明顯.

    2)地震活動性方面,蓄水前(a時段)水庫東北區(qū)離庫岸1—10km范圍有小震叢集,且小震群震源深度大多在8km左右,可能與原地斷層活動有關;蓄水后,水位首次快速加載期間(b時段),東北區(qū)離庫岸1—3km區(qū)域的小震活動有所抑制,離庫岸7—10km有小震活動,且震源深度在5km左右,可能與水庫蓄水有關,西南區(qū)小震活動則明顯增強,形成顯著的小震叢集區(qū),這可能與水庫蓄水關系密切,周斌等(2010)、盧顯等(2010)研究認為,水庫西南小震群屬于“快速響應型”水庫地震;第一次和第二次高幅加載、卸載期間(c,d時段),東北區(qū)小震活動減弱,西南區(qū)小震活動繼續(xù)增強,東南區(qū)也出現(xiàn)了部分小震叢集現(xiàn)象,周斌等(2010)認為,東南區(qū)小震群為“滯后型”水庫地震.

    按月滑動的不同震級的地震頻次與水位的關系表明,紫坪鋪水庫蓄水增強了該區(qū)的地震活動性.隨著蓄水水位的逐漸升高,蓄水4個月時ML0.0—2.0的小震迅速增多;蓄水13—15個月、水位接近840 m時,地震頻次達到峰值;蓄水15個月后、水位保持在850 m左右時,地震頻次則開始減少.可能是由于蓄水后隨著庫水的滲透作用,增加了巖石空隙水壓力.當巖石裂紋飽含水或部分含水后,水對裂紋邊界起潮濕和潤滑作用,降低了發(fā)震斷層的抗剪強度,使滑動容易產(chǎn)生,從而誘發(fā)地震.

    3)本文反演的QS與頻率無關,初始平均QS值比頻率相關的反演結果略高(張永久,趙翠萍,2009;周連慶等,2009).假設與頻率無關雖會影響QS值的數(shù)值大小,但其分布不會改變.因此,依然可以用來分析介質(zhì)的物性分布特征(Eberhart-Phillips,Chadwick,2002;汪素云等,2008;周龍泉等,2009;王惠琳等,2012).

    四川省地震局為本研究提供數(shù)據(jù)資料,劉杰研究員、蔣海昆研究員對論文給出具體的指導,張永久研究員、張致偉碩士對本文給予幫助,審稿專家對本文提出很好的修改建議和意見.作者在此一并表示感謝.

    陳颙,黃庭芳,劉恩儒.2009.巖石物理學[M].合肥:中國科學技術大學出版社:1--496.

    郭永剛,常廷改,蘇克忠.2008.汶川8.0級特大地震與紫坪鋪水庫蓄水關系的討論[J].震害防御技術,3(3):259--265.

    李亞林,賀振華,黃德濟,單鈺銘.2001.巖石孔滲特性與地震波衰減、傳播速度的相互關系[J].天然氣工業(yè),21(4):7--11.

    劉建華,胥頤,郝天珧.2004.地震波衰減的物理機制研究[J].地球物理學進展,19(1):1--7.

    劉杰,鄭斯華,黃玉龍.2003.利用遺傳算法反演非彈性衰減系數(shù)、震源參數(shù)和場地響應[J].地震學報,25(2):211--218.

    盧顯,張曉東,周龍泉,張致偉.2010.紫坪鋪水庫庫區(qū)地震精定位研究及分析[J].地震,30(2):10--19.

    施行覺,徐果明,靳平,盧振剛,劉文忠.1995.巖石的含水飽和度對縱、橫波速及衰減影響的實驗研究[J].地球物理學報,38(1):281--287.

    蘇有錦,劉杰,鄭斯華,劉麗芳,付虹,徐彥.2006.云南地區(qū)S波非彈性衰減Q值研究[J].地震學報,28(2):206--212.

    汪素云,裴順平,Thomas M Heam,許忠淮,James F Ni,俞言祥.2008.利用M L振幅研究地殼橫波Q值Ⅱ:Q橫向變化特征[J].地球物理學報,51(1):133--139.

    王大興,辛可鋒,李幼銘,高靜懷,伍向陽.2006.地層條件下砂巖含水飽和度對波速及衰減影響的實驗研究[J].地球物理學報,49(3):908--914.

    王惠琳,張曉東,周龍泉,徐曉楓,楊志高,盧顯,李智蓉.2012.紫坪鋪水庫區(qū)域地殼QS成像及其與滲透關系研究[J].地球物理學報,55(2):526--537.

    張永久,趙翠萍.2009.紫坪鋪水庫庫區(qū)介質(zhì)衰減、臺站響應和震源參數(shù)研究[J].地震地質(zhì),31(4):664--675.

    周斌,薛世峰,鄧志輝,孫峰,蔣海昆,張曉東,盧顯.2010.水庫誘發(fā)地震時空演化與庫水加卸載及滲透過程的關系:以紫坪鋪水庫為例[J].地球物理學報,53(11):2651--2670.

    周連慶,趙翠萍,陳章立.2009.紫坪鋪水庫地區(qū)尾波Qc值研究[J].地震,29(4):44--51.

    周龍泉,劉杰,蘇有錦,馬宏生,周俊杰.2009.利用S波高頻衰減參數(shù)對云南地區(qū)地殼Q值成像[J].地球物理學報,52(2):1500--1507.

    Atkinson G M,Mereu R F.1992.The shape of ground motion attenuation curves in southeastern Canada[J].Bull Seism Soc Amer,82(5):2014--2031.

    Brune J N.1970.Tectonic stress and the spectra of seismic shear waves from earthquakes[J].J Geophys Res,75(26):4997--5009.

    Cormier V F.1982.The effect of attenuation on seismic body waves[J].Bull Seism Soc Amer,72(B6):169--200.

    Eberhart-Phillips D,Chadwick M.2002.Three-dimensional attenuation model of the shallow Hikuangi subduction zone in the Raukumara Peninsula,New Zealand[J].J Geophys Res,107:2033,doi:10.1029/2000JB000046.

    Hanks T C,Wyss M.1972.The use of body wave spectra in the determination of seismic source parameters[J].Bull Seism Soc Amer,62(2):561--589.

    Hansen S,Thurber C,Mandernach M,Haslinger F,Doran C.2004.Seismic velocity and attenuation structure of the east rift zone and south Flank of Kilauea volcano,Hawaii[J].Bull Seism Soc Amer,94(4):1430--1440.

    Haberland C,Rietbrock A.2001.Attenuation tomography in the western central Andes:A detailed insight into the structure of a magmatic arc[J].J Geophys Res,106(B6):11151--11167.

    Johnston D H,Toksoz M N,Timur A.1979.Attenuation of seismic wave in dry and saturated rocks:Ⅱ.Mechanisms[J].Geophysics,44(4):691--711.

    Mavko G M,Nur A.1979.Wave attenuation in partially saturated rocks[J].Geophysics,44(2):161--178.

    Rietbrock A.2001.P wave attenuation structure in the fault area of the 1995 Kobe earthquake[J].J Geophys Res,106(B3):4141--4154.

    Scherbaum F.1990.Combined inversion for the three-dimensional Q structure and source parameters using microear thquake spectra[J].J Geophys Res,95(B8):12423--12438.

    Walsh J B.1966.Seismic waves attenuation in rock due to friction[J].J Geophys Res,71(10):2591--2599.

    Wittlinger G H,Haessler H,Granet M.1983.Three dimensional inversion of Q P from low magnitude earthquakes analysis[J].Ann Geophys,1(6):427--438.

    Relationship between the reservoir storage and crustal QSdynamic variation in Zipingpu reservoir area

    Wang Huilin1),Zhang Xiaodong2)Zhou Longquan2)Lu Xian2)Xu Xiaofeng1)
    1)Earthquake Administration of Hainan Province,Haikou 570203,China
    2)China Earthquake Networks Center,Beijing 100045,China

    10.3969/j.issn.0253-3782.2012.05.009

    P315.5

    A

    地震行業(yè)科研專項(200808036)、國家科技支撐計劃(2008BAC38B03)和地震科技星火計劃(XH12037Y)共同資助.

    2011-09-22收到初稿,2012-04-10決定采用修改稿.

    http:∥www.cnki.net/kcms/detail/11.2021.P.20120830.1420.007.html

    王惠琳,張曉東,周龍泉,盧顯,徐曉楓.2012.紫坪鋪水庫區(qū)域地殼QS動態(tài)變化及其與水庫蓄水關系的研究.地震學報,34(5):676--688.

    Wang Huilin,Zhang Xiaodong,Zhou Longquan,Lu Xian,Xu Xiaofeng.2012.Relationship between the reservoir storage and crustal QSdynamic variation in Zipingpu reservoir area.Acta Seismologica Sinica,34(5):676--688.

    猜你喜歡
    紫坪小震蓄水
    2022年巴西水電站水庫蓄水率達50%
    云南小江地區(qū)小震震源機制及構造應力場研究
    2017年9月2日固原4.6級地震前后小震調(diào)制比分析
    西藏在建大型水電站工程通過蓄水驗收
    水泵技術(2021年3期)2021-01-22 21:55:04
    跟蹤導練(五)3
    金馬河安全泄流管理體系的構建及建議
    四川水利(2017年1期)2017-03-31 09:25:13
    強震下紫坪鋪壩前大型古滑坡體變形破壞效應
    小震調(diào)制比在寧夏及鄰區(qū)映震能力的研究
    紫坪鋪電廠3號水輪發(fā)電機組下導振擺超標分析與處理
    北京20座下凹式立交橋汛前可蓄水
    麻豆国产av国片精品| 国产淫片久久久久久久久| 简卡轻食公司| 十八禁国产超污无遮挡网站| 国产毛片a区久久久久| 亚洲精品自拍成人| 长腿黑丝高跟| 亚洲精品久久国产高清桃花| 性欧美人与动物交配| 免费观看精品视频网站| 久久精品91蜜桃| 午夜a级毛片| 听说在线观看完整版免费高清| 国内精品宾馆在线| 中文在线观看免费www的网站| 伊人久久精品亚洲午夜| 亚洲精品日韩av片在线观看| 麻豆国产97在线/欧美| 亚洲欧美清纯卡通| 精品久久久久久久久久久久久| 亚洲欧美日韩东京热| 国产亚洲5aaaaa淫片| 国产成人福利小说| 国产老妇伦熟女老妇高清| 国产综合懂色| 日本-黄色视频高清免费观看| 草草在线视频免费看| 日本爱情动作片www.在线观看| 51国产日韩欧美| 精品国内亚洲2022精品成人| 你懂的网址亚洲精品在线观看 | 欧美一区二区亚洲| 国产黄色视频一区二区在线观看 | 国产女主播在线喷水免费视频网站 | 国产精品一二三区在线看| 亚洲丝袜综合中文字幕| 男的添女的下面高潮视频| 免费不卡的大黄色大毛片视频在线观看 | 只有这里有精品99| 亚洲最大成人av| 日韩 亚洲 欧美在线| 大又大粗又爽又黄少妇毛片口| 熟女电影av网| 波多野结衣高清作品| 亚洲欧美清纯卡通| 国产在线精品亚洲第一网站| 国产亚洲精品久久久久久毛片| 久久亚洲精品不卡| av卡一久久| 99热这里只有精品一区| 三级经典国产精品| 久久这里只有精品中国| 久久午夜福利片| 午夜精品国产一区二区电影 | a级毛片a级免费在线| 国产精品美女特级片免费视频播放器| 18+在线观看网站| 天天躁夜夜躁狠狠久久av| 长腿黑丝高跟| 麻豆乱淫一区二区| 欧美日韩国产亚洲二区| 最好的美女福利视频网| 欧美日韩乱码在线| 久久久久性生活片| 国产伦在线观看视频一区| 免费看a级黄色片| 黄色视频,在线免费观看| 国产麻豆成人av免费视频| 日韩一区二区视频免费看| 麻豆久久精品国产亚洲av| 国产精品国产三级国产av玫瑰| 日韩av在线大香蕉| 蜜桃久久精品国产亚洲av| 国产精品国产高清国产av| 亚洲成人中文字幕在线播放| 超碰av人人做人人爽久久| 精品无人区乱码1区二区| 联通29元200g的流量卡| 韩国av在线不卡| 又黄又爽又刺激的免费视频.| 女同久久另类99精品国产91| 国产成人91sexporn| 在线观看午夜福利视频| av在线老鸭窝| 12—13女人毛片做爰片一| 又粗又硬又长又爽又黄的视频 | 大香蕉久久网| 国产私拍福利视频在线观看| 99久久精品国产国产毛片| 又黄又爽又刺激的免费视频.| 欧美三级亚洲精品| 国产成人福利小说| 国产视频内射| 久久久久久久久久久免费av| 国产精品av视频在线免费观看| 亚洲婷婷狠狠爱综合网| 久久久久久久亚洲中文字幕| 国产免费一级a男人的天堂| 亚洲精品粉嫩美女一区| 日本撒尿小便嘘嘘汇集6| 欧美激情国产日韩精品一区| 日本色播在线视频| 能在线免费看毛片的网站| 国产大屁股一区二区在线视频| 亚洲国产高清在线一区二区三| 99热只有精品国产| 能在线免费观看的黄片| 亚洲,欧美,日韩| 欧美日韩精品成人综合77777| 观看美女的网站| 蜜桃亚洲精品一区二区三区| 国语自产精品视频在线第100页| 久久久久久伊人网av| 男人和女人高潮做爰伦理| 可以在线观看毛片的网站| 亚洲精品日韩在线中文字幕 | 日本成人三级电影网站| 深夜精品福利| 美女高潮的动态| 九草在线视频观看| 久久久久久久久大av| 久久久色成人| 久久久久久久久久成人| 丰满的人妻完整版| 国产片特级美女逼逼视频| 天堂av国产一区二区熟女人妻| 欧美又色又爽又黄视频| 国产高清有码在线观看视频| 精品久久久久久久久久久久久| 国内精品美女久久久久久| 亚洲国产欧美人成| 亚洲欧美精品综合久久99| 黄色日韩在线| 国产视频首页在线观看| 精品99又大又爽又粗少妇毛片| 亚洲成人av在线免费| 五月玫瑰六月丁香| 夜夜夜夜夜久久久久| 亚洲欧美日韩东京热| 美女国产视频在线观看| 国产精品美女特级片免费视频播放器| 伦理电影大哥的女人| 99久久中文字幕三级久久日本| 欧美xxxx黑人xx丫x性爽| 免费看日本二区| 亚洲欧美精品专区久久| 有码 亚洲区| 精品久久久久久久末码| 中文字幕精品亚洲无线码一区| 久久久a久久爽久久v久久| 全区人妻精品视频| 亚洲精品亚洲一区二区| av在线天堂中文字幕| 国产黄色视频一区二区在线观看 | 日韩一区二区视频免费看| 你懂的网址亚洲精品在线观看 | 男女边吃奶边做爰视频| 亚洲四区av| 亚洲第一电影网av| 黄色配什么色好看| 哪个播放器可以免费观看大片| 亚洲精品国产成人久久av| 麻豆成人午夜福利视频| 精品久久久久久久久久免费视频| 亚洲av二区三区四区| 男女下面进入的视频免费午夜| 亚洲av成人av| 亚洲精品久久国产高清桃花| 高清毛片免费看| h日本视频在线播放| 国产极品精品免费视频能看的| 男人舔女人下体高潮全视频| 九九热线精品视视频播放| 高清毛片免费观看视频网站| 精品99又大又爽又粗少妇毛片| 国产探花极品一区二区| av天堂在线播放| 国产一区二区亚洲精品在线观看| 国产精品嫩草影院av在线观看| 国产乱人偷精品视频| 网址你懂的国产日韩在线| 五月伊人婷婷丁香| 午夜亚洲福利在线播放| 女人被狂操c到高潮| 日日摸夜夜添夜夜添av毛片| 国产精品av视频在线免费观看| 国产真实乱freesex| 日本熟妇午夜| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 午夜精品一区二区三区免费看| 国产白丝娇喘喷水9色精品| 国产成人精品一,二区 | 精品午夜福利在线看| 成人无遮挡网站| 久久久国产成人精品二区| 成人欧美大片| 1000部很黄的大片| 亚洲精品国产成人久久av| 一级av片app| h日本视频在线播放| 日韩欧美三级三区| 亚洲四区av| 国内精品久久久久精免费| 国产一级毛片七仙女欲春2| 能在线免费看毛片的网站| 婷婷色av中文字幕| 久久精品国产亚洲av涩爱 | 免费看av在线观看网站| 国产免费男女视频| 国产一区二区三区av在线 | 九草在线视频观看| 亚洲欧洲国产日韩| videossex国产| 99热这里只有精品一区| 少妇熟女aⅴ在线视频| 男人舔女人下体高潮全视频| 男插女下体视频免费在线播放| 久久精品国产亚洲av天美| 欧美最新免费一区二区三区| 人人妻人人澡欧美一区二区| 久久精品综合一区二区三区| 亚洲国产精品成人久久小说 | 国产黄a三级三级三级人| 久久久久网色| 国产av一区在线观看免费| 女同久久另类99精品国产91| www.色视频.com| 日本av手机在线免费观看| 精品无人区乱码1区二区| 大又大粗又爽又黄少妇毛片口| 精品一区二区三区人妻视频| АⅤ资源中文在线天堂| 成人av在线播放网站| 国产欧美日韩精品一区二区| 国产精品综合久久久久久久免费| 国产成人精品婷婷| 亚洲国产色片| 欧美潮喷喷水| 久久欧美精品欧美久久欧美| 中国国产av一级| 桃色一区二区三区在线观看| 欧美+日韩+精品| 精品人妻视频免费看| 1024手机看黄色片| 黄色日韩在线| 最近中文字幕高清免费大全6| 黄色欧美视频在线观看| 亚洲综合色惰| 国产91av在线免费观看| 只有这里有精品99| 熟女人妻精品中文字幕| 亚洲电影在线观看av| 天美传媒精品一区二区| 日韩欧美在线乱码| 天堂√8在线中文| 午夜福利成人在线免费观看| 男女啪啪激烈高潮av片| 男插女下体视频免费在线播放| 日本熟妇午夜| 我要搜黄色片| 尾随美女入室| 久久久精品欧美日韩精品| 亚洲av成人精品一区久久| 少妇丰满av| 麻豆乱淫一区二区| 在线a可以看的网站| 在线免费观看的www视频| 国产亚洲精品久久久com| 午夜免费男女啪啪视频观看| avwww免费| 日韩欧美精品免费久久| 亚洲三级黄色毛片| 国产三级在线视频| 狠狠狠狠99中文字幕| 国产探花在线观看一区二区| 亚洲自拍偷在线| 久久九九热精品免费| 国产免费男女视频| 菩萨蛮人人尽说江南好唐韦庄 | 国产精品av视频在线免费观看| 狠狠狠狠99中文字幕| 亚洲人成网站在线播| 在线a可以看的网站| 深爱激情五月婷婷| 可以在线观看毛片的网站| 高清日韩中文字幕在线| 99热这里只有精品一区| 五月伊人婷婷丁香| 干丝袜人妻中文字幕| 国产女主播在线喷水免费视频网站 | 精华霜和精华液先用哪个| 51国产日韩欧美| 成人特级黄色片久久久久久久| 一区二区三区四区激情视频 | 天天躁日日操中文字幕| 免费看光身美女| 一级毛片电影观看 | 日韩在线高清观看一区二区三区| 哪个播放器可以免费观看大片| 亚洲欧美日韩无卡精品| av天堂在线播放| 国产高清有码在线观看视频| 亚洲国产精品久久男人天堂| 国产极品精品免费视频能看的| 91在线精品国自产拍蜜月| 黄色欧美视频在线观看| 91久久精品电影网| 久久久久久国产a免费观看| 欧美色欧美亚洲另类二区| 国产高清不卡午夜福利| 色综合站精品国产| 如何舔出高潮| 99视频精品全部免费 在线| 成年女人永久免费观看视频| 婷婷色av中文字幕| 久久久久国产网址| 淫秽高清视频在线观看| 直男gayav资源| 天天躁夜夜躁狠狠久久av| 91久久精品国产一区二区成人| 一级毛片aaaaaa免费看小| 少妇丰满av| 精品一区二区三区视频在线| 国产激情偷乱视频一区二区| 国产精品福利在线免费观看| 青青草视频在线视频观看| a级毛色黄片| 尾随美女入室| 国产成人91sexporn| 能在线免费观看的黄片| 欧美日韩国产亚洲二区| 日韩欧美一区二区三区在线观看| 亚洲国产精品成人久久小说 | 91麻豆精品激情在线观看国产| 成人性生交大片免费视频hd| 欧美+日韩+精品| 色吧在线观看| 99热全是精品| 亚洲色图av天堂| 色哟哟·www| 亚洲欧美精品综合久久99| 欧美精品一区二区大全| 简卡轻食公司| 此物有八面人人有两片| 韩国av在线不卡| 国产三级中文精品| 免费在线观看成人毛片| 亚洲精品色激情综合| 高清毛片免费看| 最新中文字幕久久久久| 蜜桃久久精品国产亚洲av| 欧美bdsm另类| 免费黄网站久久成人精品| 黄片wwwwww| 成人亚洲精品av一区二区| 亚洲国产精品sss在线观看| 亚洲成a人片在线一区二区| 蜜臀久久99精品久久宅男| 三级毛片av免费| 国内少妇人妻偷人精品xxx网站| av女优亚洲男人天堂| 寂寞人妻少妇视频99o| 亚洲av免费高清在线观看| 亚洲欧美精品自产自拍| 成人性生交大片免费视频hd| 国产乱人视频| 99热全是精品| 成人三级黄色视频| 欧美另类亚洲清纯唯美| 久久婷婷人人爽人人干人人爱| 成人亚洲精品av一区二区| 看非洲黑人一级黄片| 校园春色视频在线观看| 精品一区二区免费观看| 久久久a久久爽久久v久久| 丰满人妻一区二区三区视频av| av黄色大香蕉| 国产精品人妻久久久影院| 天堂网av新在线| 久久人人爽人人爽人人片va| 青春草视频在线免费观看| 床上黄色一级片| 成人午夜精彩视频在线观看| a级毛片a级免费在线| 日韩在线高清观看一区二区三区| 最近视频中文字幕2019在线8| 91精品一卡2卡3卡4卡| 欧美一区二区亚洲| 欧美不卡视频在线免费观看| 六月丁香七月| 亚洲中文字幕日韩| 国产人妻一区二区三区在| 91久久精品电影网| a级毛色黄片| 亚洲自拍偷在线| 欧美成人精品欧美一级黄| 婷婷色综合大香蕉| 亚洲第一电影网av| 国产高清视频在线观看网站| 欧美最黄视频在线播放免费| 亚洲无线在线观看| 一个人看视频在线观看www免费| 精品免费久久久久久久清纯| 看十八女毛片水多多多| 国产精华一区二区三区| 国产精品综合久久久久久久免费| 99热6这里只有精品| 国产亚洲5aaaaa淫片| 国产单亲对白刺激| 久久精品国产亚洲av香蕉五月| 级片在线观看| 一边亲一边摸免费视频| 欧美日本视频| 欧美激情在线99| 人人妻人人澡人人爽人人夜夜 | 99热只有精品国产| 99久久精品一区二区三区| 一边亲一边摸免费视频| 久久人妻av系列| 成人国产麻豆网| 夜夜夜夜夜久久久久| 人人妻人人澡人人爽人人夜夜 | 国产v大片淫在线免费观看| 2021天堂中文幕一二区在线观| 又爽又黄无遮挡网站| 在线观看一区二区三区| 中文在线观看免费www的网站| av又黄又爽大尺度在线免费看 | 久久久久国产网址| 午夜老司机福利剧场| 不卡一级毛片| 日韩一区二区三区影片| 欧美日韩精品成人综合77777| 亚洲图色成人| 国产一区亚洲一区在线观看| 精品久久久久久久人妻蜜臀av| 狂野欧美激情性xxxx在线观看| 欧美日韩国产亚洲二区| 六月丁香七月| 国产一级毛片七仙女欲春2| 51国产日韩欧美| 国产精品无大码| 午夜精品国产一区二区电影 | 免费在线观看成人毛片| 天堂av国产一区二区熟女人妻| 男女边吃奶边做爰视频| 一个人看视频在线观看www免费| 毛片一级片免费看久久久久| 三级男女做爰猛烈吃奶摸视频| 一区二区三区免费毛片| av在线蜜桃| 国产中年淑女户外野战色| 深夜a级毛片| 亚洲不卡免费看| 五月玫瑰六月丁香| 日日摸夜夜添夜夜添av毛片| 国产精品久久久久久av不卡| 麻豆久久精品国产亚洲av| 久久人人爽人人片av| 五月伊人婷婷丁香| 国产在线精品亚洲第一网站| 综合色av麻豆| 国产精品一区二区三区四区免费观看| 久久国内精品自在自线图片| 亚洲真实伦在线观看| 午夜福利成人在线免费观看| 五月玫瑰六月丁香| 91狼人影院| 亚洲高清免费不卡视频| 国产女主播在线喷水免费视频网站 | 精品久久国产蜜桃| 国产伦一二天堂av在线观看| 97热精品久久久久久| 两个人的视频大全免费| 干丝袜人妻中文字幕| 久久99热这里只有精品18| 精品久久久噜噜| 黄色视频,在线免费观看| 白带黄色成豆腐渣| 亚洲国产精品成人久久小说 | 国产亚洲av片在线观看秒播厂 | 亚洲真实伦在线观看| 免费黄网站久久成人精品| 亚洲国产欧洲综合997久久,| 国产精品99久久久久久久久| 欧美成人一区二区免费高清观看| 久久精品人妻少妇| 亚洲久久久久久中文字幕| 亚洲熟妇中文字幕五十中出| 特级一级黄色大片| 一本久久精品| 亚洲真实伦在线观看| 色尼玛亚洲综合影院| 中文精品一卡2卡3卡4更新| 久久亚洲国产成人精品v| 美女内射精品一级片tv| 天天躁夜夜躁狠狠久久av| 老司机福利观看| 精品午夜福利在线看| 国产成人一区二区在线| 久久久久九九精品影院| 亚洲精品久久久久久婷婷小说 | 免费人成在线观看视频色| 午夜亚洲福利在线播放| 免费观看在线日韩| 国产视频内射| 美女黄网站色视频| 亚洲精华国产精华液的使用体验 | 人妻夜夜爽99麻豆av| 国产69精品久久久久777片| 国产成人91sexporn| 少妇的逼水好多| 国产成人a∨麻豆精品| 一进一出抽搐gif免费好疼| 晚上一个人看的免费电影| 亚洲自偷自拍三级| 一个人看视频在线观看www免费| 最近手机中文字幕大全| 日本与韩国留学比较| 又爽又黄a免费视频| or卡值多少钱| 人体艺术视频欧美日本| 韩国av在线不卡| 久久久成人免费电影| 少妇的逼水好多| 嫩草影院入口| 国产精品久久久久久亚洲av鲁大| 高清毛片免费看| 丰满乱子伦码专区| 22中文网久久字幕| 亚洲av不卡在线观看| 国产男人的电影天堂91| 亚洲国产精品合色在线| 成人av在线播放网站| 少妇丰满av| 99久久精品一区二区三区| 久久韩国三级中文字幕| 男的添女的下面高潮视频| 国产91av在线免费观看| 色综合站精品国产| 变态另类成人亚洲欧美熟女| 久久久久性生活片| 美女黄网站色视频| 国产精品人妻久久久久久| 97热精品久久久久久| 蜜臀久久99精品久久宅男| 国产综合懂色| 天美传媒精品一区二区| 国产女主播在线喷水免费视频网站 | 国产高清有码在线观看视频| 边亲边吃奶的免费视频| 久久精品久久久久久久性| 黄片wwwwww| 99久久成人亚洲精品观看| 99在线人妻在线中文字幕| 黄色一级大片看看| 久久久久九九精品影院| 国产爱豆传媒在线观看| 女人被狂操c到高潮| 国产精品一区二区三区四区久久| 九草在线视频观看| 国产麻豆成人av免费视频| 国产精品久久久久久精品电影| 97超碰精品成人国产| 国模一区二区三区四区视频| 人人妻人人澡欧美一区二区| 亚洲在线观看片| 69人妻影院| 国产高潮美女av| 亚洲图色成人| 美女xxoo啪啪120秒动态图| 永久网站在线| 午夜精品国产一区二区电影 | 免费观看精品视频网站| 夜夜爽天天搞| 波多野结衣高清无吗| 级片在线观看| 搡老妇女老女人老熟妇| 悠悠久久av| 内地一区二区视频在线| 亚洲精品国产成人久久av| 高清午夜精品一区二区三区 | 久久久久久久久久黄片| 看免费成人av毛片| 一进一出抽搐动态| 我的女老师完整版在线观看| 特级一级黄色大片| 精品人妻熟女av久视频| 日韩制服骚丝袜av| 久久久久久久久久成人| 人妻少妇偷人精品九色| 久久久久久伊人网av| 日韩国内少妇激情av| 国产精品美女特级片免费视频播放器| 久久久久久伊人网av| 不卡视频在线观看欧美| 国产麻豆成人av免费视频| 99久久中文字幕三级久久日本| 亚洲最大成人中文| 国产中年淑女户外野战色| 深夜精品福利| 网址你懂的国产日韩在线| 桃色一区二区三区在线观看| 99在线人妻在线中文字幕| 一级黄色大片毛片| 国产精品久久久久久久久免| 久久久久久久午夜电影| 九九在线视频观看精品| 国产v大片淫在线免费观看| 成人特级av手机在线观看| 免费观看在线日韩| 欧美高清性xxxxhd video| 国产一区二区三区在线臀色熟女| 日本-黄色视频高清免费观看| 欧美高清性xxxxhd video|