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

    榮昌及周邊三維速度結(jié)構(gòu)初步研究

    2022-04-16 03:14:08李翠平唐茂云郭衛(wèi)英王小龍
    地震地質(zhì) 2022年1期
    關(guān)鍵詞:雙差榮昌震源

    李翠平 唐茂云 郭衛(wèi)英 王小龍 董 蕾

    (重慶市地震局, 重慶 401147)

    0 引言

    經(jīng)中國地震臺網(wǎng)中心測定, 北京時間2016年12月27日08時17分, 重慶榮昌區(qū)發(fā)生MS4.9 地震, 震中位于四川盆地東緣的華鎣山基底斷裂系中北段。該區(qū)域地震活動性原本很弱, 自20世紀(jì)80年代末起, 由于天然氣田廢水回注的影響, 出現(xiàn)4.0級地震頻發(fā)的現(xiàn)象, 目前已發(fā)生了14次MS≥4.0地震, 最大為1997年榮昌MS5.2 地震(圖 1)。從地質(zhì)構(gòu)造背景來看, 華鎣山基底斷裂以西的川中穹窿地區(qū)變形微弱、 斷層發(fā)育不明顯, 以鼻狀或短穹隆狀背斜構(gòu)造為特點。而其東側(cè)為川東陷褶束, 表現(xiàn)為沉積蓋層在結(jié)晶基底之上的滑脫變形, 產(chǎn)生一系列狹長背斜和寬大平緩向斜相間排列的侏羅山式隔擋式褶皺以及與背斜構(gòu)造共生的中等規(guī)模斷續(xù)分布的早—中更新世地表和隱伏斷層, 長度多為幾千米至幾十千米(丁仁杰等, 2004; 李峰等, 2013)。區(qū)域內(nèi)中等地震頻發(fā)被認(rèn)為是人為注水使這些中等規(guī)模的斷層發(fā)生再活動所致(Leietal., 2008; 王小龍等, 2015; 王志偉等, 2018)。多名學(xué)者先后分析了該區(qū)域的地震活動與注水的關(guān)系, Lei等(2008)基于ETAS模型對1988—2006年榮昌地區(qū)的地震活動進行了詳細(xì)研究, 認(rèn)為榮昌誘發(fā)地震序列的物理機制是伴隨深井注水的孔隙壓力擴散作用和先前地震本身引起的庫侖破裂應(yīng)力變化聯(lián)合作用的結(jié)果。此外, 有研究通過傳統(tǒng)統(tǒng)計學(xué)方法分析榮昌地區(qū)注水和地震活動(1993—2006)的關(guān)系, 發(fā)現(xiàn)兩者存在很好的相關(guān)性(黃世源等, 2007; 朱麗霞等, 2007)。2010年9月10日MS4.5 和2014年2月23日ML3.8地震的震源參數(shù)結(jié)果表明其震源深度與注水井深度一致, 為注水活動誘發(fā)的構(gòu)造地震活動(王小龍等, 2012, 2015)。王志偉等(2018)采用匹配定位方法對榮昌地區(qū)2008—2011年的微小地震進行檢測和定位, 結(jié)合ML≥3.5地震的震源機制解, 表明地震主要分布在幾條隱伏斷層附近, 與注水活動導(dǎo)致已有斷層的再活動有關(guān)。

    圖 1 研究區(qū)臺站、 歷史震中分布圖Fig. 1 Seismic stations and historical earthquakes in the study area.

    前人多對研究區(qū)2012年以前的地震活動進行分析, 之后該區(qū)域又相繼發(fā)生了5次4.0級以上地震, 特別是2016年榮昌MS4.9 地震, 對當(dāng)?shù)鼐用裨斐闪艘欢ǖ挠绊懀?針對此次地震的孕震環(huán)境以及人為注水活動對震區(qū)地殼介質(zhì)結(jié)構(gòu)的影響都亟待開展進一步研究。此外, 現(xiàn)有的研究結(jié)果在該區(qū)域震源分布與三維速度結(jié)構(gòu)之間的關(guān)系方面涉及較少, 這為認(rèn)識榮昌地區(qū)深部介質(zhì)環(huán)境以及破裂過程的深部構(gòu)造成因機制帶來困難。本文基于2008年1月—2020年6月重慶區(qū)域地震臺網(wǎng)、 流動地震臺網(wǎng)和自貢地方臺網(wǎng)記錄的地震觀測報告, 開展榮昌及周邊精細(xì)的三維速度結(jié)構(gòu)和地震重定位研究, 為更深入地了解該區(qū)域中等地震的發(fā)震構(gòu)造和注水活動影響范圍提供基礎(chǔ)信息。

    1 方法和數(shù)據(jù)

    1.1 地震數(shù)據(jù)

    表 1 榮昌及周邊斷裂性質(zhì)Table1 The fault properties in the Rongchang and its adjacent areas

    1.2 雙差層析成像方法

    地震層析成像是一種利用大量地震觀測數(shù)據(jù)反演地下介質(zhì)的速度以及震源參數(shù)等重要信息的方法。本文基于獲得的絕對到時和相對到時數(shù)據(jù), 采用Zhang等(2003, 2006)提出的雙差層析成像方法進行速度結(jié)構(gòu)反演。該方法在雙差定位算法(Waldhauseretal., 2000)的基礎(chǔ)上加入速度參數(shù)項, 減少了雙差定位中假設(shè)臺站到事件對之間路徑為恒定速度所引起的誤差, 通過改變速度模型重新計算震源參數(shù), 使得定位精度進一步提高。同時相對精確的地震定位可以降低由于射線路徑的不同造成的速度結(jié)構(gòu)反演的誤差, 從而提高速度結(jié)構(gòu)反演的精度, 其獲得的分辨尺度可達3~5km, 能夠刻畫出較小尺度范圍內(nèi)比較精細(xì)的速度結(jié)構(gòu)(Zhangetal., 2003; 王小娜等, 2015), 可用于研究地震活動與速度結(jié)構(gòu)的關(guān)系(肖卓等, 2017; 韓曉明等, 2018; 劉白云等, 2018; 王長在等, 2018; Longetal., 2020)和庫水滲透影響作用(Dixitetal., 2015; 吳海波等, 2018), 所得結(jié)果將為分析地震孕育發(fā)生環(huán)境等提供重要信息。

    圖 2 a 研究區(qū)斷層、 注水井、 網(wǎng)格節(jié)點; b 初始一維速度模型Fig. 2 The distribution map of faults, injection wells and grid nodes(a)and initial 1D velocity model(b) in the study area.綠色六邊形代表注水井, 六邊形內(nèi)的數(shù)字代表注水井編號, 其中注水井1螺4井, 2螺2井, 3包11井, 4包24井, 5包18井, 6 包23井(圖7和圖8相同)

    同時, 雙差層析成像方法在反演解算中采用阻尼最小二乘分解算法(LSQR)求解, 該算法中的光滑權(quán)重(smooth)和阻尼(damp)參數(shù)可約束地震位置和慢度的變化量, 影響反演結(jié)果的穩(wěn)定性。因此, 在實際數(shù)據(jù)反演前, 要對不同的光滑權(quán)重和阻尼參數(shù)進行權(quán)衡分析(Eberhart, 1986; 王小娜等, 2015), 利用L曲線法進行不同光滑權(quán)重和阻尼值的測試, 通過歸一化模型與走時殘差關(guān)系曲線分析獲得最佳的光滑權(quán)重和阻尼參數(shù)。分析顯示, 當(dāng)光滑權(quán)重和阻尼參數(shù)分別為20和300時(圖 3)模型較平滑, 且走時殘差也相對較小, 故在反演成像中采用這2個值。

    圖 3 利用L曲線選擇最佳的光滑權(quán)重和阻尼參數(shù)Fig. 3 The best smoothing weight and damping parameters selected by trade-off curves.

    在實際數(shù)據(jù)反演過程中共進行7輪迭代, 不斷地調(diào)整迭代參數(shù), 使WRCT和WDCT值不斷降低, 并根據(jù)迭代結(jié)果調(diào)整阻尼值, 使得表征雙差方程條件數(shù)的CND為40~80。被定位到空氣中的個數(shù)(AQ)從初始的12個減少到結(jié)束時的0個, 地震走時殘差均方根(RMS)也由初始的299ms降為36ms, 說明迭代過程趨于穩(wěn)定合理。重定位前地震走時殘差均方根的主要分布范圍為200~800ms, 重新定位后地震走時殘差均方根主要集中在0~150ms(圖 4), 定位精度明顯提高。

    圖 4 定位前(a)、 后(b)0~1200ms范圍內(nèi)地震走時殘差均方根直方圖Fig. 4 Histograms of travel time residual RMS between 0~1200ms for seismic events before(a)and after(b)the relocation.

    2 結(jié)果分析

    2.1 分辨率及可靠性分析

    圖 5 研究區(qū)內(nèi)不同深度的檢測板結(jié)果Fig. 5 Checkerboard test results at different depths in the study area.

    圖 6 沿X、 Y方向速度剖面的檢測板結(jié)果Fig. 6 Checkerboard test results along the direction of X and Y.

    圖 7 深度為0~16km處的P波速度結(jié)構(gòu)及地震分布Fig. 7 Distribution of P wave velocity structure and earthquake in 0~16km deep layer.紅色五角星表示層面上下各1.5km范圍MS>4.0的地震(圖8相同)

    2.2 P波速度結(jié)構(gòu)

    采用雙差層析成像方法得到了榮昌及周邊不同深度層和不同剖面上P波速度結(jié)構(gòu)分布圖??紤]該區(qū)域地震的震源深度和檢測板結(jié)果, 本文主要針對上地殼20km以內(nèi)的速度結(jié)構(gòu)特征進行分析。圖 7 為研究區(qū)不同深度層的P波速度分布圖, 將每條速度剖面上下1.5km范圍內(nèi)的地震事件投影到剖面上。其中1km深度層的P波高速區(qū)分布在螺觀山背斜的北端, 該區(qū)域發(fā)育燕子巖斷層、 螺觀山南翼隱伏斷層和廣順橫斷層等多條斷層, 螺4井也分布在該區(qū)域, 該注水井自1984年開始向約2km的深度注入廢水, 一直持續(xù)到2001年。但由于1km深度處的分辨率結(jié)果有限, 該高速區(qū)的分布僅作參考。4km深度層P波高速區(qū)分布在仁義—榮昌一帶, 該區(qū)域分布著4口注水井, 是自2008年以來榮昌持續(xù)注水的主要集中區(qū)域(魏紅梅等, 2014; 王小龍等, 2015)。7km深度層的速度結(jié)構(gòu)與淺表層相比有明顯變化, 仁義—榮昌高速區(qū)的范圍明顯縮小, 華鎣山基底斷裂的西側(cè)盤龍和隆昌地區(qū)表現(xiàn)為高速區(qū), 這與基于接收函數(shù)反演得到的盤龍測震臺在6~8km出現(xiàn)高速層的結(jié)果一致(王小龍等, 2012)。而華鎣山斷裂東側(cè)榮昌—雙河地區(qū)則表現(xiàn)為大范圍的低速區(qū), 這種高速區(qū)和低速區(qū)的差異分布一直持續(xù)到13km深度。華鎣山斷裂以東的川東陷褶束基底主要由一套巨厚的沉積變質(zhì)碎屑巖夾碳酸鹽巖與火山碎屑巖的復(fù)理石建造所組成, 屬低密度和弱至無磁性的塑性基底結(jié)構(gòu), 而以西的川中臺拱基底是由密度高、 磁性強的一套基性中性及較強磁性的火成巖組成(趙從俊, 1984)。本文反演的速度結(jié)構(gòu)結(jié)果顯示以華鎣山斷裂帶為界, 7~13km深度處兩側(cè)呈現(xiàn)截然不同的速度分布, 這與前人的研究結(jié)果基本吻合。

    圖 8 研究區(qū)不同剖面的P波速度結(jié)構(gòu)Fig. 8 P wave velocity structure along the different profiles.紅色倒三角表示斷層在地表的位置

    圖 8 展示了不同剖面上的P波速度結(jié)構(gòu)分布, 顯示出3個明顯特征: 1)12km深度處為明顯的速度界面, 其上速度起伏變化較大, 其下速度變化較平穩(wěn)。根據(jù)榮昌地震反射剖面結(jié)果可知, 榮昌地區(qū)的上地殼為典型的薄皮構(gòu)造, 表現(xiàn)為沉積蓋層在基底面之上滑脫、 變形, 并發(fā)育多條逆沖斷層和褶皺, 基底埋深10~11km, 最深可達12km(王贊軍等, 2018)。本文結(jié)果中12km深處的速度界面對應(yīng)了上地殼沉積蓋層和結(jié)晶基底的分界深度層。2)在螺2井位置下3~7km深度有明顯的高速體。螺2井注水時間為1991—2012年, 周邊發(fā)育埋深1.7km的螺觀山北翼隱伏斷層, 高速體沿著該斷層的傾向分布, 可能表明螺觀山北翼隱伏斷層是廢水下滲的主要斷裂通道。3)Y=-10km和Y=0km剖面的速度結(jié)構(gòu)能夠較好地反映螺觀山褶皺的形態(tài), 也進一步論證了結(jié)果的可靠性。

    2.3 地震活動與速度結(jié)構(gòu)的關(guān)系

    基于雙差層析成像反演得到了榮昌及周邊2008年以來相對精確的地震時空展布特征, 重定位后的地震叢集性更加明顯, 除了2020年分布在榮昌與瀘縣交界附近的地震外, 其余地震主要分布在變形強烈的螺觀山背斜軸部, 呈明顯的NE向帶狀分布, 與背斜和斷層走向一致(圖7a, b), 優(yōu)勢深度為0~6km(圖 9)。該區(qū)域受廢水回注的影響, 周邊的地表和隱伏斷層易活化, 而使地震活動與斷層表現(xiàn)較強的空間相關(guān)性。2009年8月8日MS4.3 、 2010年2月22日MS4.3 和2010年9月10日MS4.5 地震位于地表燕子巖斷層上, 震源機制解較為一致, 均顯示為NE走向的逆沖型地震(表2)。垂直于斷層走向的P波速度剖面(圖8c)顯示地震呈SE傾向的帶狀分布, 位于高、 低速分界線上, 發(fā)生在其底端的2010年9月10日MS4.5 地震則更偏向于北西側(cè)的高速區(qū)。利用根據(jù)小地震分布確定斷層面參數(shù)的方法對該地震條帶進行斷層面走向和傾角參數(shù)的求解(萬永革等, 2008), 得到地震斷層面的走向為42°, 傾角為80°, 傾向SE, 斷層深度為1~6km, 對應(yīng)MS4.5 地震震源機制解的節(jié)面Ⅱ。該地震條帶的發(fā)震斷層與螺觀山南翼隱伏斷層和燕子巖斷層的產(chǎn)狀相吻合, 由于地震主要位于斷層的南東盤且1km深度以淺地震分布較少, 與螺觀山南翼隱伏斷層的位置和埋深一致, 故推測螺觀山南翼隱伏斷層應(yīng)為其發(fā)震斷層。2009年8月8日MS4.5 地震也有類似的規(guī)律(圖8e)。2016年12月27日MS4.9 和2016年12月28日MS4.0 地震則分布在以螺觀山北翼隱伏斷層為界限的高、 低速交界地區(qū)(圖8b, f), 該處地表未見斷層出露, 僅存在螺觀山北翼隱伏斷層, 該斷層走向225°, 傾角45°, 傾向NW(丁仁杰等, 2004), 與2016年12月27日MS4.9 的震源機制解的節(jié)面Ⅰ相吻合。

    圖 9 初始定位(a)和重定位(b)震源深度統(tǒng)計直方圖Fig. 9 The statistical histogram of focal depth before(a) and after(b) relocation.

    表 2 研究區(qū)MS≥4.0地震的震源機制解Table2 Focal mechanism solutions of MS≥4.0 earthquake in the study area

    2.4 注水對榮昌地區(qū)速度結(jié)構(gòu)的影響

    3 結(jié)論

    本文基于2008年1月—2020年6月重慶區(qū)域臺網(wǎng)固定臺站、 流動臺網(wǎng)及自貢地方臺網(wǎng)記錄到的大量地震事件的P波、 S波到時數(shù)據(jù), 采用雙差層析成像方法反演獲得了榮昌及周邊精細(xì)的三維P波速度結(jié)構(gòu)和地震事件的重定位結(jié)果, 得到的主要認(rèn)識如下:

    (1)榮昌及周邊上地殼為典型的薄皮構(gòu)造, 為變形強烈的沉積蓋層與結(jié)晶基底的雙層結(jié)構(gòu), 本文的P波速度結(jié)構(gòu)結(jié)果顯示該地區(qū)12km深度為2層的分界深度, 與榮昌地震反射剖面結(jié)果一致。4km深度以內(nèi)的高、 低速區(qū)分布與7~13km以下深度層存在明顯差別, 7~13km深度層的速度結(jié)構(gòu)顯示華鎣山基底斷裂東、 西兩側(cè)的川中臺拱和川東陷褶束分別表現(xiàn)為低速區(qū)和高速區(qū), 與前人的研究結(jié)果基本一致。

    (2)榮昌淺表層P波速度結(jié)構(gòu)橫向差異變化較大, 在螺觀山背斜軸部的螺2注水井附近為高速區(qū), 分布深度為3~7km, 該井周邊發(fā)育埋深1.7km的螺觀山北翼隱伏斷層, P波高速體主要沿著螺觀山北翼隱伏斷層的傾向分布, 可能表明該斷層是廢水下滲的主要斷裂通道, 下滲深度達7km, 使該斷層兩側(cè)存在較大的速度差異。2016年12月27日MS4.9 和2016年12月28日MS4.0 地震就分布在該速度過渡帶上。

    (3)地震重定位結(jié)果顯示地震主要分布在變形強烈的螺觀山背斜軸部, 呈明顯的NE向帶狀分布, 震源優(yōu)勢深度為0~6km。以廣順橫斷層為界, 南側(cè)的地震表現(xiàn)為傾向SE的地震條帶, 結(jié)合震源機制解和區(qū)域地震構(gòu)造環(huán)境分析, 廣順橫斷層南側(cè)的MS>4.0地震的發(fā)震斷層應(yīng)為螺觀山南翼隱伏斷層, 而廣順橫斷層北側(cè)的地震活動可能與螺觀山北翼隱伏斷層有關(guān)。

    致謝中國科學(xué)技術(shù)大學(xué)張海江教授提供了雙差層析成像程序; 重慶地震臺網(wǎng)中心提供了地震觀測報告; 審稿專家為本文提供了修改意見和建議; 本文圖件使用GMT及Origin軟件繪制。在此一并表示感謝!

    猜你喜歡
    雙差榮昌震源
    虛擬地震臺陣雙差測深法及應(yīng)用
    呂繼成:給榮昌陶注入新活力
    工會博覽(2022年3期)2022-06-30 03:59:12
    BDS中長基線三頻RTK算法研究
    BDS參考站間低高度角衛(wèi)星三頻整周模糊度解算方法
    全國首個榮昌豬數(shù)據(jù)標(biāo)準(zhǔn)體系發(fā)布
    確保十九大精神全面落實在榮昌大地上開創(chuàng)“成渝城市群新興戰(zhàn)略支點”新篇章——專訪中共榮昌區(qū)委書記曹清堯
    基于雙差的波動方程反射波旅行時反演方法
    震源的高返利起步
    人類學(xué)視角中的榮昌陶藝——關(guān)于榮昌陶藝的田野調(diào)查
    可控震源地震在張掖盆地南緣逆沖斷裂構(gòu)造勘探中的應(yīng)用
    精品久久久久久久久久久久久| 成年av动漫网址| 国产探花在线观看一区二区| 国语自产精品视频在线第100页| 亚洲欧美精品自产自拍| 18+在线观看网站| 99久久人妻综合| 亚洲久久久久久中文字幕| 国产色爽女视频免费观看| 国产av一区在线观看免费| 亚洲在线观看片| a级一级毛片免费在线观看| 伦精品一区二区三区| 不卡视频在线观看欧美| 伦理电影大哥的女人| 国产白丝娇喘喷水9色精品| 岛国毛片在线播放| 亚洲av免费高清在线观看| 男人的好看免费观看在线视频| 成人午夜高清在线视频| 久久韩国三级中文字幕| 一区二区三区免费毛片| av在线播放精品| av卡一久久| 丝袜美腿在线中文| 91久久精品国产一区二区成人| 亚洲欧美日韩高清专用| 免费观看精品视频网站| ponron亚洲| 又粗又硬又长又爽又黄的视频| 最后的刺客免费高清国语| 亚洲国产精品久久男人天堂| av福利片在线观看| 九九久久精品国产亚洲av麻豆| 男女视频在线观看网站免费| 国产乱人偷精品视频| 白带黄色成豆腐渣| 亚洲av.av天堂| 国产精品久久视频播放| 成人美女网站在线观看视频| 欧美一区二区精品小视频在线| 国产男人的电影天堂91| 久99久视频精品免费| 精品国内亚洲2022精品成人| 亚洲av中文av极速乱| 色哟哟·www| 亚洲国产精品成人久久小说| 少妇的逼水好多| 久久热精品热| 中国国产av一级| 久久人人爽人人片av| 免费观看精品视频网站| 亚洲精品自拍成人| 日韩强制内射视频| 18禁在线无遮挡免费观看视频| 成人鲁丝片一二三区免费| 免费黄网站久久成人精品| 91精品一卡2卡3卡4卡| 波多野结衣高清无吗| 综合色av麻豆| 国产精品久久久久久久久免| av福利片在线观看| 国产黄a三级三级三级人| 国产亚洲av嫩草精品影院| 观看免费一级毛片| 亚洲第一区二区三区不卡| 在线免费观看的www视频| 午夜福利视频1000在线观看| 看片在线看免费视频| 天天躁夜夜躁狠狠久久av| 国产精品三级大全| 最近的中文字幕免费完整| 国产午夜精品一二区理论片| 亚洲精品成人久久久久久| 91在线精品国自产拍蜜月| 女的被弄到高潮叫床怎么办| 欧美日韩综合久久久久久| 亚洲av.av天堂| 久久国内精品自在自线图片| 麻豆一二三区av精品| 男人舔奶头视频| 国产亚洲最大av| 人妻系列 视频| 尾随美女入室| 成人国产麻豆网| 亚洲五月天丁香| 日韩欧美国产在线观看| 日本免费一区二区三区高清不卡| 欧美性感艳星| 精品久久久久久久久久久久久| 日本免费在线观看一区| 国产精品一区二区三区四区久久| 少妇人妻精品综合一区二区| 九草在线视频观看| 亚洲av免费在线观看| 久久久久网色| 亚洲无线观看免费| 搡女人真爽免费视频火全软件| videossex国产| 久久久久久久亚洲中文字幕| ponron亚洲| 亚洲av电影在线观看一区二区三区 | 久久精品国产99精品国产亚洲性色| 18禁在线播放成人免费| 男插女下体视频免费在线播放| 亚洲美女搞黄在线观看| 少妇高潮的动态图| 亚洲精品日韩av片在线观看| 国产一区亚洲一区在线观看| 日韩强制内射视频| 久久午夜福利片| 99热全是精品| 蜜桃亚洲精品一区二区三区| 久久精品国产鲁丝片午夜精品| 国产伦在线观看视频一区| 人体艺术视频欧美日本| 最新中文字幕久久久久| av在线播放精品| 国产乱来视频区| 国产淫语在线视频| 色网站视频免费| 91aial.com中文字幕在线观看| 国产精品三级大全| 国产亚洲精品久久久com| 久久精品国产亚洲av天美| 亚洲在久久综合| 日日啪夜夜撸| 高清av免费在线| 午夜福利视频1000在线观看| 99久久成人亚洲精品观看| 一个人看视频在线观看www免费| 超碰av人人做人人爽久久| 人妻少妇偷人精品九色| 熟妇人妻久久中文字幕3abv| 在线免费观看不下载黄p国产| 丰满乱子伦码专区| 黄片wwwwww| 国产午夜精品一二区理论片| 麻豆成人午夜福利视频| 免费人成在线观看视频色| 乱系列少妇在线播放| 欧美日韩国产亚洲二区| 成人鲁丝片一二三区免费| 国产精品久久视频播放| 精品久久久久久久久av| 男人舔奶头视频| 亚洲精品乱码久久久久久按摩| 日韩强制内射视频| 久久久国产成人精品二区| 亚洲欧美日韩东京热| 色噜噜av男人的天堂激情| 日韩高清综合在线| 在线播放无遮挡| 亚洲精品日韩av片在线观看| 能在线免费看毛片的网站| 老女人水多毛片| 边亲边吃奶的免费视频| 久久精品久久精品一区二区三区| 日韩欧美国产在线观看| 汤姆久久久久久久影院中文字幕 | 淫秽高清视频在线观看| 内地一区二区视频在线| 男的添女的下面高潮视频| 一级毛片电影观看 | 免费观看在线日韩| av在线天堂中文字幕| 午夜精品国产一区二区电影 | 一级爰片在线观看| 精品人妻视频免费看| 亚洲aⅴ乱码一区二区在线播放| 久久精品国产鲁丝片午夜精品| h日本视频在线播放| 亚洲一级一片aⅴ在线观看| 简卡轻食公司| 久久久久久久久久久免费av| 久久久久久大精品| 在线播放国产精品三级| 蜜臀久久99精品久久宅男| 91在线精品国自产拍蜜月| 久久精品国产鲁丝片午夜精品| av在线蜜桃| av在线天堂中文字幕| 日本三级黄在线观看| 免费观看在线日韩| 一区二区三区四区激情视频| 联通29元200g的流量卡| 中文欧美无线码| 久久午夜福利片| 一级毛片aaaaaa免费看小| 欧美日本亚洲视频在线播放| 国产精华一区二区三区| 久久鲁丝午夜福利片| 99热网站在线观看| 欧美3d第一页| 日韩欧美 国产精品| 久久久精品欧美日韩精品| 成人亚洲精品av一区二区| 插逼视频在线观看| 精品免费久久久久久久清纯| 美女被艹到高潮喷水动态| av播播在线观看一区| 国产高清有码在线观看视频| 久久精品熟女亚洲av麻豆精品 | 高清日韩中文字幕在线| 国产69精品久久久久777片| 亚洲色图av天堂| 国产成人精品婷婷| 极品教师在线视频| 男女啪啪激烈高潮av片| 欧美成人一区二区免费高清观看| 国产伦理片在线播放av一区| 日韩三级伦理在线观看| 中文字幕亚洲精品专区| 又粗又爽又猛毛片免费看| 亚洲欧美日韩卡通动漫| 麻豆av噜噜一区二区三区| 欧美激情在线99| 国产三级在线视频| 久久这里只有精品中国| 蜜臀久久99精品久久宅男| 最近中文字幕2019免费版| 国产黄片美女视频| 永久网站在线| 人人妻人人澡欧美一区二区| 国产一区二区在线观看日韩| 淫秽高清视频在线观看| 美女xxoo啪啪120秒动态图| 如何舔出高潮| 亚洲精品乱码久久久久久按摩| av视频在线观看入口| 亚洲三级黄色毛片| 亚洲伊人久久精品综合 | 国产在线男女| 在线免费观看的www视频| 三级国产精品片| .国产精品久久| 亚洲国产欧美人成| 亚洲天堂国产精品一区在线| 国产一区二区在线观看日韩| 青春草国产在线视频| 久久久久久久久久久免费av| 亚洲第一区二区三区不卡| 免费黄网站久久成人精品| 观看美女的网站| 99九九线精品视频在线观看视频| 国产视频首页在线观看| 亚洲综合色惰| 天堂√8在线中文| 搡老妇女老女人老熟妇| 日本一二三区视频观看| 欧美成人精品欧美一级黄| 亚洲欧美成人精品一区二区| 在线观看66精品国产| 久久久久久大精品| 亚洲三级黄色毛片| 亚洲欧洲国产日韩| 色视频www国产| 国产极品精品免费视频能看的| 一夜夜www| 国产免费又黄又爽又色| 久久婷婷人人爽人人干人人爱| 中文字幕精品亚洲无线码一区| 九草在线视频观看| 偷拍熟女少妇极品色| 国产淫语在线视频| 爱豆传媒免费全集在线观看| 亚洲图色成人| 亚洲乱码一区二区免费版| 2021少妇久久久久久久久久久| 亚洲在久久综合| 高清av免费在线| 久久精品国产亚洲av天美| 蜜臀久久99精品久久宅男| 国产精品熟女久久久久浪| 大香蕉久久网| 久久精品久久精品一区二区三区| 麻豆成人av视频| 男女边吃奶边做爰视频| 国产淫语在线视频| 国产成人精品婷婷| 国产乱人偷精品视频| 天天躁夜夜躁狠狠久久av| 亚洲熟妇中文字幕五十中出| 日韩av在线免费看完整版不卡| 亚洲av.av天堂| 国语对白做爰xxxⅹ性视频网站| 一级爰片在线观看| 日韩一本色道免费dvd| 免费观看性生交大片5| 日韩欧美在线乱码| 国产精品日韩av在线免费观看| 亚洲精品一区蜜桃| 亚洲av成人精品一二三区| 久久精品久久久久久久性| 三级国产精品欧美在线观看| 成人无遮挡网站| 亚洲在久久综合| 日韩强制内射视频| 免费av观看视频| 天堂√8在线中文| 国产 一区 欧美 日韩| 日韩中字成人| 欧美成人免费av一区二区三区| 日韩av在线免费看完整版不卡| 色播亚洲综合网| 99热这里只有是精品50| 波多野结衣高清无吗| 69人妻影院| 久久精品国产亚洲av天美| 日本黄色视频三级网站网址| 51国产日韩欧美| 久久久久久久亚洲中文字幕| 成人毛片a级毛片在线播放| 国产一区二区在线av高清观看| 女人被狂操c到高潮| 国产精品.久久久| 长腿黑丝高跟| 麻豆精品久久久久久蜜桃| 国产午夜精品论理片| 久久精品国产自在天天线| 97热精品久久久久久| 国产日韩欧美在线精品| 婷婷色麻豆天堂久久 | 国产伦精品一区二区三区四那| 2021天堂中文幕一二区在线观| 国产免费视频播放在线视频 | 麻豆一二三区av精品| 免费观看a级毛片全部| 美女内射精品一级片tv| 一级毛片aaaaaa免费看小| 2021天堂中文幕一二区在线观| 国产精品久久久久久av不卡| 免费人成在线观看视频色| 美女大奶头视频| 最近手机中文字幕大全| 尾随美女入室| 精品一区二区免费观看| 人妻少妇偷人精品九色| 久久鲁丝午夜福利片| 久久久久久伊人网av| 国产乱来视频区| 女人久久www免费人成看片 | 久久精品久久久久久久性| 伊人久久精品亚洲午夜| 国产私拍福利视频在线观看| 国产黄色视频一区二区在线观看 | 两个人视频免费观看高清| 国产 一区 欧美 日韩| 亚洲真实伦在线观看| 99久国产av精品国产电影| 国产老妇女一区| 欧美成人a在线观看| 成人高潮视频无遮挡免费网站| 欧美性感艳星| 内射极品少妇av片p| 欧美性感艳星| 日韩国内少妇激情av| 日本免费a在线| 日日摸夜夜添夜夜添av毛片| 免费大片18禁| 亚洲国产精品合色在线| 一二三四中文在线观看免费高清| 一级毛片aaaaaa免费看小| 精品国产三级普通话版| 中文乱码字字幕精品一区二区三区 | 久久精品久久久久久噜噜老黄 | 两性午夜刺激爽爽歪歪视频在线观看| 中文欧美无线码| 一二三四中文在线观看免费高清| 秋霞伦理黄片| 3wmmmm亚洲av在线观看| 美女cb高潮喷水在线观看| av在线亚洲专区| 久久久久久久久久黄片| 免费在线观看成人毛片| 成人高潮视频无遮挡免费网站| 色吧在线观看| 精品人妻熟女av久视频| 久久久久久久久久成人| 亚洲精品乱久久久久久| 精品欧美国产一区二区三| 欧美高清成人免费视频www| 亚洲av成人精品一二三区| 99热6这里只有精品| 国产精品99久久久久久久久| 中文精品一卡2卡3卡4更新| 女人被狂操c到高潮| 久久久久久久久中文| 亚洲精品成人久久久久久| 久久久久九九精品影院| 国产亚洲精品av在线| 国产精品人妻久久久影院| 韩国高清视频一区二区三区| 色网站视频免费| 中文字幕熟女人妻在线| 男女视频在线观看网站免费| 少妇人妻精品综合一区二区| 国产毛片a区久久久久| 中文字幕熟女人妻在线| 国产激情偷乱视频一区二区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产淫语在线视频| 中文欧美无线码| 在现免费观看毛片| 秋霞伦理黄片| 久久久午夜欧美精品| 欧美xxxx黑人xx丫x性爽| 欧美+日韩+精品| 一级毛片电影观看 | 波多野结衣高清无吗| 国产亚洲精品久久久com| 国产精品一区二区在线观看99 | 男女边吃奶边做爰视频| 日韩欧美国产在线观看| 国产三级在线视频| 亚洲va在线va天堂va国产| 九草在线视频观看| 免费观看人在逋| 看黄色毛片网站| 最近最新中文字幕免费大全7| 99久久人妻综合| 深夜a级毛片| 建设人人有责人人尽责人人享有的 | 直男gayav资源| 国产精品福利在线免费观看| 人妻少妇偷人精品九色| 午夜福利网站1000一区二区三区| 能在线免费看毛片的网站| 免费无遮挡裸体视频| 日韩成人av中文字幕在线观看| 国产午夜精品一二区理论片| 国产亚洲91精品色在线| 九九爱精品视频在线观看| 97人妻精品一区二区三区麻豆| 狂野欧美白嫩少妇大欣赏| 亚洲中文字幕日韩| 免费观看人在逋| 欧美高清性xxxxhd video| 亚洲国产精品久久男人天堂| 欧美一区二区精品小视频在线| 91在线精品国自产拍蜜月| av专区在线播放| 国产精品人妻久久久影院| 亚洲欧美一区二区三区国产| 国产白丝娇喘喷水9色精品| 亚洲av电影在线观看一区二区三区 | 联通29元200g的流量卡| 国产三级在线视频| 中文乱码字字幕精品一区二区三区 | 欧美激情在线99| av在线亚洲专区| 国产精品精品国产色婷婷| 好男人视频免费观看在线| 插逼视频在线观看| 51国产日韩欧美| 国产真实伦视频高清在线观看| 中文字幕亚洲精品专区| 亚洲av免费在线观看| 国产欧美日韩精品一区二区| 搡老妇女老女人老熟妇| 欧美区成人在线视频| 老司机影院成人| 午夜日本视频在线| 99热这里只有精品一区| av又黄又爽大尺度在线免费看 | 国产一区亚洲一区在线观看| 看黄色毛片网站| 一区二区三区四区激情视频| 亚洲国产高清在线一区二区三| 成人性生交大片免费视频hd| 日韩欧美精品v在线| 亚洲国产精品久久男人天堂| 亚洲内射少妇av| 国产午夜精品论理片| 亚洲高清免费不卡视频| 国内少妇人妻偷人精品xxx网站| 免费在线观看成人毛片| 99久久精品热视频| 麻豆一二三区av精品| 两个人视频免费观看高清| 欧美日韩一区二区视频在线观看视频在线 | 波多野结衣高清无吗| 老司机福利观看| 国产精品一及| 久久久欧美国产精品| 日日干狠狠操夜夜爽| 日本黄色视频三级网站网址| a级毛色黄片| 日韩制服骚丝袜av| 秋霞伦理黄片| 欧美zozozo另类| 大话2 男鬼变身卡| 国产69精品久久久久777片| 菩萨蛮人人尽说江南好唐韦庄 | 99久久中文字幕三级久久日本| 精品久久久久久久人妻蜜臀av| 搡老妇女老女人老熟妇| 天堂中文最新版在线下载 | 日本av手机在线免费观看| 中文亚洲av片在线观看爽| 亚洲人与动物交配视频| 在线观看一区二区三区| 可以在线观看毛片的网站| 日韩成人av中文字幕在线观看| 国产av码专区亚洲av| 亚洲人与动物交配视频| 草草在线视频免费看| 可以在线观看毛片的网站| av在线播放精品| 成年女人看的毛片在线观看| 三级国产精品片| 日韩制服骚丝袜av| 日韩欧美在线乱码| 亚洲国产最新在线播放| 亚洲精品日韩在线中文字幕| 在线播放国产精品三级| 噜噜噜噜噜久久久久久91| av天堂中文字幕网| 热99re8久久精品国产| 亚洲av成人精品一二三区| 日韩一区二区视频免费看| 免费观看性生交大片5| 久久精品91蜜桃| 亚洲欧美日韩高清专用| 国产精品一区二区性色av| 岛国在线免费视频观看| 热99re8久久精品国产| 国内少妇人妻偷人精品xxx网站| 熟女人妻精品中文字幕| 久久久亚洲精品成人影院| 亚洲成人久久爱视频| 一区二区三区高清视频在线| 禁无遮挡网站| 全区人妻精品视频| 精品人妻熟女av久视频| 国产精品无大码| 欧美极品一区二区三区四区| 亚洲国产精品专区欧美| av在线老鸭窝| 欧美性猛交黑人性爽| 亚洲精华国产精华液的使用体验| 日本爱情动作片www.在线观看| 三级男女做爰猛烈吃奶摸视频| 大香蕉97超碰在线| 精品熟女少妇av免费看| 亚洲国产色片| 日韩一本色道免费dvd| 日韩 亚洲 欧美在线| 亚洲欧美精品专区久久| 国产一级毛片在线| 国产不卡一卡二| 亚洲国产成人一精品久久久| 日本wwww免费看| 国产日韩欧美在线精品| 久久久久久国产a免费观看| 国产午夜精品论理片| 天美传媒精品一区二区| 精品少妇黑人巨大在线播放 | 丝袜喷水一区| 丰满人妻一区二区三区视频av| 国产乱人偷精品视频| 免费电影在线观看免费观看| 午夜精品一区二区三区免费看| 国产成人免费观看mmmm| av在线老鸭窝| 男女边吃奶边做爰视频| 伦精品一区二区三区| 女人十人毛片免费观看3o分钟| av福利片在线观看| 99久久精品国产国产毛片| 国产爱豆传媒在线观看| 国产精品熟女久久久久浪| 男人舔奶头视频| 又爽又黄无遮挡网站| 精品久久久久久电影网 | 亚洲av免费在线观看| 看黄色毛片网站| 女人十人毛片免费观看3o分钟| 伦理电影大哥的女人| 丰满少妇做爰视频| 大话2 男鬼变身卡| 国产老妇伦熟女老妇高清| 日韩精品青青久久久久久| 国产精品伦人一区二区| 免费播放大片免费观看视频在线观看 | 麻豆乱淫一区二区| 久久99精品国语久久久| 国产精品久久久久久久久免| 精品久久久久久电影网 | 男人舔女人下体高潮全视频| 亚洲精品乱码久久久v下载方式| 亚洲怡红院男人天堂| 国产69精品久久久久777片| 两个人的视频大全免费| 国产精华一区二区三区| 身体一侧抽搐| 日本黄色片子视频| 亚洲国产精品sss在线观看| 男人狂女人下面高潮的视频| 26uuu在线亚洲综合色| 国产久久久一区二区三区| 美女国产视频在线观看| 真实男女啪啪啪动态图| 97在线视频观看| 免费看光身美女| 成人无遮挡网站| 国产精品综合久久久久久久免费| 国产中年淑女户外野战色| 欧美日本亚洲视频在线播放| 午夜福利成人在线免费观看| 国产欧美另类精品又又久久亚洲欧美| 免费一级毛片在线播放高清视频| 欧美性猛交黑人性爽| 内射极品少妇av片p| 六月丁香七月|