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

    水力壓裂誘發(fā)的剪張型微地震震源機制矩張量反演方法

    2019-08-06 08:54:38戚瑞軒
    石油地球物理勘探 2019年4期
    關(guān)鍵詞:張量震源振幅

    唐 杰 溫 雷 李 聰 戚瑞軒

    (①中國石油大學(華東)地球科學與技術(shù)學院,山東青島 266580;②海洋國家實驗室海洋礦產(chǎn)資源評價與探測技術(shù)功能實驗室,山東青島 266071)

    0 引言

    本文首先介紹利用P波初動振幅信息反演水力壓裂裂縫震源機制解的方法原理,然后通過理論模型記錄測試、分析方法的應用效果。

    1 矩張量理論分析

    1.1 源張量與矩張量

    采用斷面法向矢量n=(n1,n2,n3)T及位錯方向矢量v=(v1,v2,v3)T表示源張量[15]

    (1)

    式中P為源強度。

    由剪張源模型(圖1)可知

    (2)

    (3)

    走向角φ、 傾角η、 滑動角θ和張裂角α(圖1)通常是空間位置和時間的函數(shù),為了簡化起見,在研究剪張錯動的點源近似時可視為常量。走向角、傾角和滑動角描述了斷面法向和斷面位錯矢量的切向,張裂角描述了位錯矢量偏離斷面的程度。張裂角α∈[-90°,90°],即:當α>0°時為剪張源,對于純張裂源,α=90°;當α<0°時,為壓裂源,對于純壓裂源,α=-90°;對于純剪切源,α=0°。

    圖1 剪張源模型

    微地震震源機制可用矩張量的形式表示[16]

    (4)

    式中:M為對稱矩陣,滿足Mij=Mji,其中獨立元素有6個,Mij表示沿著±i方向大小相等并且在j方向相距無限小的力偶(i、j=1,2,3);M0為地震矩幅度。

    對于各向同性介質(zhì)而言,矩張量分量Mij和源張量分量Dij滿足

    Mij=λDkkδij+2μDij

    (5)

    式中:λ和μ為震源區(qū)介質(zhì)的拉梅常數(shù);Dkk表示下標相同的源張量分量;δ為克羅內(nèi)克函數(shù),當i=j時,δij=1,當i≠j時,δij=0。輻射花樣可觀測破裂面的幾何分布和破裂機制類型,本質(zhì)上是對沙灘球顯示的拓展。圖2為不同張裂角震源的縱、橫波輻射花樣圖。

    1.2 矩張量分解

    對微地震震源矩張量M采用特征值分解法分解為雙力偶部分MDC、補償線性矢量偶極成分MCLVD和各向同性部分MISO[17]

    M=MISO+MCLVD+MDC

    (6)

    其中

    圖2 不同張裂角震源的縱、橫波輻射花樣圖

    各分量所占的比例滿足

    (7)

    式中Mmax為M的最大特征值。在各向同性介質(zhì)中,非DC分量決定源的類型,即:當非DC分量都為正值時,源是張裂源(裂紋張開);當非DC分量都為負值時,源是壓縮源(裂紋閉合);DC分量總是正的;當為純剪切源時DC分量占比為100%。

    1.3 矩張量圖示

    Hudson等[19]定義了一種震源類型的圖解表示法,可實現(xiàn)各種破裂模式的概念化和可視化,可以研究矩張量中剪切組分和張性組分的分布情況。首先對矩張量進行正則化分解,然后計算其特征值Mj,并按大小排序

    M1≥M2≥M3

    (8)

    (9)

    (10)

    通過下列參數(shù)化運算,可以計算張性組分k和剪切組分T

    (11)

    圖3 不同張裂角(α=90°、45°、10°、0°、-10°、-45°、-90°)

    橫坐標表示剪切組分T值,縱坐標表示張性組分k值,兩者的取值范圍均為-1~1。各向同性爆炸型震源(ISO)中的膨脹震源(E)和壓縮震源(I)的T-k坐標分別為(0,1)和(0,-1);純剪切源的坐標為(0,0),位于圖的中心位置;兩類錯動方向相反的補償線性偶極子震源(CLVD+與CLVD-)的T-k坐標分別為(-1,0)和(1,0)

    (12)

    從而將矩張量的6個獨立元素簡化為Hudson圖中的點(T,k)[20](圖3)。

    2 基于振幅信息的矩張量反演

    微地震震源矩張量反演方法很多,包括采用極性、振幅以及全波形或集中考慮特定的體波相位(如P波和S波)的反演方法等。根據(jù)輸入數(shù)據(jù)的質(zhì)量和數(shù)量,可獲得全張量或者部分張量,如果限制反演為偏量部分能夠減少未知參數(shù)的數(shù)目,從而增強反演問題的穩(wěn)定性[21-22]。

    矩張量反演利用點源在均勻彈性介質(zhì)中P波和S波的質(zhì)點運動方程計算格林函數(shù)。P波和S波在x點、t時刻的振幅分別滿足

    (13)

    (14)

    式中:i為分量號;r為源到接收器的距離;ρ為密度;VP為P波速度;VS為S波速度;γi、γj和γk為從源到接收器的方位余弦;Mjk為矩張量分量;δij為克羅內(nèi)克函數(shù)。

    忽略震源時間函數(shù)時求解式(13)、式(14)的問題簡化為

    d=Gm

    式中:d為包含N個振幅的矢量;G為格林矩陣,表征震源和接收器之間的振幅變化;m為地震矩矢量,包含矩張量的6個獨立參數(shù)。上式的矩陣形式為

    (15)

    通過

    m=G-1d

    計算矩張量。上式為超定方程,求解需要更多的觀測數(shù)據(jù)。為此,采用G-g代替G-1進行反演

    G-g=[GTG]-1GT

    (16)

    由于實際地震數(shù)據(jù)中含有噪聲,因此需要6個以上的振幅值準確求解m。在獲得m后,需要穩(wěn)定性分析以證明解的可靠性。反演的穩(wěn)定性依賴于很多方面,如接收器的分布以及信號的信噪比等。如果要計算剪張源的角度參數(shù),可以通過非線性反演方法獲得,或者先求出Mij,再計算各個角度參數(shù)。針對震源機制的反演問題,Li等[23]和Tan等[24]研究了高頻波形匹配和非線性反演方法 。

    實際矩張量和反演矩張量的誤差為

    (17)

    3 微地震矩張量反演測試分析

    矩張量反演的置信度取決于多種因素。信噪比的局限性造成的輸入信號振幅的不確定性、信號初動估計不準與定位不準導致的不確定性以及對穿過界面的反射(透射)損失考慮不周全、射線路徑聚焦效應考慮不當、地震衰減與各向異性解釋不當?shù)榷紩绊懢貜埩糠囱萁Y(jié)果。輻射花樣的震源球采樣對矩張量反演結(jié)果影響很大,測井資料由于井的數(shù)據(jù)觀測范圍有限和不規(guī)則方向孔徑使矩張量反演結(jié)果不穩(wěn)定。在沒有先驗約束條件時,利用單觀測井數(shù)據(jù)通常不足以唯一地確定矩張量;當有3口觀測井數(shù)據(jù)時單獨使用P波或P波與S波聯(lián)合反演可以獲得全部矩張量;當有2口觀測井數(shù)據(jù)時則須采用P波與S波聯(lián)合反演才可獲得全部矩張量[25]。地面監(jiān)測能夠在各個方向提供良好的震源球覆蓋,單分量地面接收器可獲得縱波振幅和初動,若沒有橫波信息約束,也會影響矩張量反演的穩(wěn)定性。

    3.1 數(shù)據(jù)誤差對反演結(jié)果的影響

    設(shè)定速度模型(圖4a)的觀測系統(tǒng)臺站為星形排列(圖4b),震源位于地下2600m處,選取相應震源機制(圖4c),獲得的MDC、MISO和MCLVD的百分比分別為61.34%、21.48%和17.18%(圖4d)。圖5為理論模型矩張量反演結(jié)果。由圖可見:基于振幅信息的偏量解在源中沒有體積變化,因此不包含MISO,僅包含MCLVD和MDC(圖5c);基于振幅信息的MDC解在源中沒有體積變化和MCLVD成分,因此僅包含MDC成分(圖5d)。

    圖4 理論模型的Hudson圖

    圖5 理論模型矩張量反演結(jié)果

    本文的反演結(jié)果由振幅信息獲得,可根據(jù)反演結(jié)果繪制沙灘球,沙灘球使用縱波初動表征的震源機制圖。震源機制反演可得到一組節(jié)面,其中一個為斷面,另一個為輔助面。對于純剪切源(DC)來說,這兩個面是等價的,還需要借助其他輔助信息(如應力或地質(zhì)信息等)確定具體節(jié)面類型。對于地面微地震監(jiān)測數(shù)據(jù)處理,可以采用P波振幅信息反演震源機制。地面微地震監(jiān)測具有較高的覆蓋次數(shù)、較寬的方位角以及較大炮檢距等特點,因此利用P波振幅信息可以得到穩(wěn)定的震源機制解,文中采用全張量反演、偏量解和DC解反演方法。為了測試不同參數(shù)誤差的影響,采用Monte-Carlo測試方法,對于不同誤差均進行100次隨機試驗,采用MISO分量和MCLVD分量百分比分布圖以及矩張量Hudson圖與鉆石圖展示反演結(jié)果。根據(jù)輸入數(shù)據(jù)的質(zhì)量和數(shù)量,可獲得全張量或者部分張量,如果限制反演為偏量部分可減少未知參數(shù)的數(shù)目、增強反演的穩(wěn)定性。

    受環(huán)境干擾影響,地面微地震監(jiān)測資料信噪比低,考慮隨機噪聲對振幅的影響,首先分析極性誤差和振幅誤差的影響。 圖6是振幅誤差為20%時的全矩張量反演結(jié)果,圖7是極性誤差為5%時的全矩張量反演結(jié)果。對比兩圖可見,振幅誤差導致反演結(jié)果偏差較小(圖6),極性誤差導致反演結(jié)果偏差較大(圖7)。采用P波振幅反演結(jié)果的MISO反演誤差相對較小,而全矩張量反演的MCLVD誤差相對較大,觀測誤差主要影響MCLVD的結(jié)果。原因為:矩張量分解為MISO和MCLVD時,MCLVD由計算最小特征值/最大特征值的比值獲得,MCLVD對反演矩張量的精度更敏感;MISO和矩張量特征值之間滿足線性關(guān)系,因此相對于MCLVD,MISO對振幅誤差不敏感。

    3.2 模型、定位精度及觀測儀器對反演結(jié)果的影響

    研究表明,模型參數(shù)的正確性、震源—接收器的幾何分布以及定位精度等影響反演結(jié)果的可靠性,接收器的分布和選擇影響反演結(jié)果的穩(wěn)定性。圖8為接收器20%重采樣時的全矩張量反演結(jié)果。由圖可見,將數(shù)據(jù)隨機舍去20%,實現(xiàn)20%重采樣,觀測系統(tǒng)對非DC分量的反演結(jié)果非常重要,合適的觀測系統(tǒng)可增加地震矩張量DC和非DC分量分析的可信性。

    圖6 振幅誤差為20%時的全矩張量反演結(jié)果

    圖9為定位誤差小于30m時的全矩張量反演結(jié)果。由圖可見,當存在定位誤差時,理論計算結(jié)果和反演結(jié)果的均方根誤差較小,非DC分量也受影響,采用P波時誤差較大,非DC分量易被虛假的DC分量所掩蓋。

    圖10為速度參數(shù)誤差小于10%時的全矩張量反演結(jié)果。由圖可見,當速度參數(shù)誤差較大時,采用P波振幅時非DC分量的誤差較大。

    圖7 極性誤差為5%時的全矩張量反演結(jié)果

    圖8 接收器20%重采樣時的全矩張量反演結(jié)果

    圖9 定位誤差小于30m時的全矩張量反演結(jié)果

    圖10 速度參數(shù)誤差小于10%時的全矩張量反演結(jié)果

    3.3 不同震源類型的反演結(jié)果分析

    不同震源類型的輻射花樣間存在差異,因此有限的震源球覆蓋導致反演多解性。圖11為不同震源類型在振幅誤差為20%及極性誤差為5%時的全矩張量反演結(jié)果。由圖可知,張裂角α=90°時受到極性和振幅誤差的影響較大,節(jié)面分布較雜亂(圖11b、圖11e),與前文的結(jié)果類似,相對于MCLVD的結(jié)果而言,MISO的誤差更小。

    圖11 不同震源類型在振幅誤差為20%及極性誤差為5%時的全矩張量反演結(jié)果

    4 結(jié)論

    本文研究了利用P波初動振幅信息反演水力壓裂裂縫震源機制解的方法原理,并通過理論模型記錄測試、分析了方法的應用效果。基于P波輻射花樣進行微地震震源機制反演,避免了傳統(tǒng)P波初動極性反演方法的不確定性,可以得到更可靠的壓裂裂縫解釋結(jié)果。獲得以下認識:

    (1)矩張量反演的置信度取決于多種因素,振幅誤差導致反演結(jié)果偏差較小,極性誤差導致反演結(jié)果偏差較大。相對于CLVD分量,采用全矩張量反演時ISO分量對誤差不敏感,這是由于在矩張量分解時,ISO分量是矩張量特征值的線性函數(shù),而CLVD分量是其非線性函數(shù)。

    (2)觀測系統(tǒng)對于非DC分量的反演結(jié)果非常重要,合適的觀測系統(tǒng)可增加地震矩張量DC和非DC分量分析的可信度。當模型誤差較大時,采用P波振幅反演時非DC分量的誤差較大。

    (3)輻射花樣能夠觀測破裂面的幾何分布和破裂機制類型,不同震源類型的輻射花樣間存在差異,因此有限的震源球覆蓋導致反演多解性。

    本文利用P波初至振幅信息反演震源機制全矩張量,由于地面資料信噪比低、數(shù)據(jù)量大、觀測范圍大,只利用P波資料反演可以得到相對穩(wěn)定的震源機制解,若利用P波、S波振幅信息聯(lián)合反演可進一步提高反演的穩(wěn)定性,同時在一定程度上能消除傳播路徑、近地表效應和儀器響應等因素對反演結(jié)果的影響。

    猜你喜歡
    張量震源振幅
    偶數(shù)階張量core逆的性質(zhì)和應用
    四元數(shù)張量方程A*NX=B 的通解
    震源的高返利起步
    十大漲跌幅、換手、振幅、資金流向
    十大漲跌幅、換手、振幅、資金流向
    擴散張量成像MRI 在CO中毒后遲發(fā)腦病中的應用
    十大漲跌幅、換手、振幅、資金流向
    滬市十大振幅
    可控震源地震在張掖盆地南緣逆沖斷裂構(gòu)造勘探中的應用
    同步可控震源地震采集技術(shù)新進展
    欧美日韩亚洲国产一区二区在线观看 | 亚洲中文av在线| 亚洲黑人精品在线| 两人在一起打扑克的视频| 成人影院久久| 不卡av一区二区三区| 黄色 视频免费看| h视频一区二区三区| 99re6热这里在线精品视频| 午夜福利免费观看在线| 欧美精品一区二区免费开放| 亚洲精品国产色婷婷电影| 又黄又粗又硬又大视频| 精品人妻一区二区三区麻豆| 亚洲精品国产精品久久久不卡| a级毛片在线看网站| 18禁黄网站禁片午夜丰满| 青春草亚洲视频在线观看| 久久精品亚洲熟妇少妇任你| 少妇精品久久久久久久| 国产三级黄色录像| 麻豆国产av国片精品| 精品一区在线观看国产| 在线 av 中文字幕| 免费少妇av软件| 色94色欧美一区二区| 中文欧美无线码| 可以免费在线观看a视频的电影网站| 我的亚洲天堂| 啦啦啦啦在线视频资源| 国产欧美日韩一区二区三 | 久久国产亚洲av麻豆专区| 久久久精品区二区三区| 亚洲精品美女久久av网站| 国产精品成人在线| 亚洲精品久久午夜乱码| 国产高清视频在线播放一区 | 69精品国产乱码久久久| 在线av久久热| 免费黄频网站在线观看国产| 欧美老熟妇乱子伦牲交| 女人久久www免费人成看片| 久久精品国产亚洲av高清一级| 亚洲精品自拍成人| 丝袜人妻中文字幕| 久久av网站| 国产av一区二区精品久久| 日韩三级视频一区二区三区| 亚洲av电影在线进入| 国产福利在线免费观看视频| 亚洲一码二码三码区别大吗| 最新在线观看一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 纯流量卡能插随身wifi吗| 十八禁网站网址无遮挡| 最黄视频免费看| 免费在线观看日本一区| 国产欧美日韩一区二区三区在线| 性高湖久久久久久久久免费观看| av片东京热男人的天堂| a级毛片黄视频| 国产精品久久久av美女十八| 久久免费观看电影| 亚洲欧洲日产国产| 五月天丁香电影| 久久国产精品男人的天堂亚洲| 美女脱内裤让男人舔精品视频| 国产精品秋霞免费鲁丝片| 丁香六月欧美| 新久久久久国产一级毛片| 91成年电影在线观看| 中文字幕人妻熟女乱码| 亚洲国产欧美网| 乱人伦中国视频| 国产精品九九99| 99香蕉大伊视频| 精品国产乱子伦一区二区三区 | 啪啪无遮挡十八禁网站| 亚洲欧洲精品一区二区精品久久久| 亚洲精品成人av观看孕妇| 人妻 亚洲 视频| 男女边摸边吃奶| 久久久久精品国产欧美久久久 | 亚洲欧美精品自产自拍| 国产精品欧美亚洲77777| 手机成人av网站| 欧美黑人精品巨大| 最近最新免费中文字幕在线| 久久ye,这里只有精品| 亚洲成人免费电影在线观看| www.av在线官网国产| 99香蕉大伊视频| 国产免费视频播放在线视频| av天堂在线播放| 免费高清在线观看日韩| 日本vs欧美在线观看视频| 精品国产国语对白av| videos熟女内射| 中文字幕人妻丝袜制服| 老司机福利观看| 下体分泌物呈黄色| 亚洲成av片中文字幕在线观看| 亚洲美女黄色视频免费看| 高潮久久久久久久久久久不卡| 超碰97精品在线观看| 十八禁网站网址无遮挡| 色婷婷av一区二区三区视频| 99久久99久久久精品蜜桃| 极品少妇高潮喷水抽搐| 精品少妇黑人巨大在线播放| 亚洲美女黄色视频免费看| av天堂久久9| a级毛片黄视频| 国产av一区二区精品久久| 久久天躁狠狠躁夜夜2o2o| 久久久精品94久久精品| 亚洲成人免费电影在线观看| 久久国产精品大桥未久av| 欧美黄色淫秽网站| 又紧又爽又黄一区二区| 久久久国产精品麻豆| 最黄视频免费看| 亚洲国产欧美网| 美女福利国产在线| 99国产极品粉嫩在线观看| 久久国产精品大桥未久av| 国产成人欧美在线观看 | 美女福利国产在线| 精品视频人人做人人爽| 波多野结衣一区麻豆| 真人做人爱边吃奶动态| 精品视频人人做人人爽| 精品熟女少妇八av免费久了| 9191精品国产免费久久| 午夜影院在线不卡| 美女大奶头黄色视频| av一本久久久久| 亚洲欧美精品综合一区二区三区| 亚洲av成人不卡在线观看播放网 | 欧美日韩精品网址| av免费在线观看网站| 成人av一区二区三区在线看 | 国产男女内射视频| 亚洲中文日韩欧美视频| 午夜91福利影院| 国产男女内射视频| 岛国在线观看网站| 涩涩av久久男人的天堂| 19禁男女啪啪无遮挡网站| 狠狠狠狠99中文字幕| 成人手机av| 亚洲第一欧美日韩一区二区三区 | 亚洲精品中文字幕一二三四区 | 成人手机av| 操出白浆在线播放| 国产97色在线日韩免费| 一级毛片电影观看| 啦啦啦在线免费观看视频4| 欧美乱码精品一区二区三区| 免费在线观看视频国产中文字幕亚洲 | 婷婷色av中文字幕| 悠悠久久av| 人妻 亚洲 视频| 国产免费福利视频在线观看| 亚洲精品日韩在线中文字幕| 日本wwww免费看| 亚洲综合色网址| 欧美精品av麻豆av| 欧美人与性动交α欧美软件| 亚洲精品成人av观看孕妇| 亚洲第一av免费看| 老司机影院毛片| 狠狠狠狠99中文字幕| 国产成人精品久久二区二区91| a在线观看视频网站| 视频区图区小说| 国产亚洲午夜精品一区二区久久| 人妻人人澡人人爽人人| 午夜视频精品福利| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品自产拍在线观看55亚洲 | 亚洲午夜精品一区,二区,三区| 曰老女人黄片| 亚洲欧美激情在线| 国产在线一区二区三区精| 日本91视频免费播放| 国产一卡二卡三卡精品| 岛国毛片在线播放| 肉色欧美久久久久久久蜜桃| 男女边摸边吃奶| 美女主播在线视频| 一本色道久久久久久精品综合| 成人影院久久| 热99国产精品久久久久久7| 国产高清视频在线播放一区 | 免费观看a级毛片全部| 国产区一区二久久| 日本猛色少妇xxxxx猛交久久| 窝窝影院91人妻| 日韩 亚洲 欧美在线| 国产精品国产三级国产专区5o| 十八禁高潮呻吟视频| 1024香蕉在线观看| 精品国产一区二区三区久久久樱花| 国产黄色免费在线视频| 麻豆av在线久日| 欧美国产精品一级二级三级| 99久久综合免费| 国产精品1区2区在线观看. | 亚洲成人手机| 精品国内亚洲2022精品成人 | 一级毛片电影观看| 我的亚洲天堂| 老司机在亚洲福利影院| 欧美 日韩 精品 国产| 波多野结衣av一区二区av| 久久精品国产亚洲av香蕉五月 | 精品国产乱码久久久久久男人| 精品人妻熟女毛片av久久网站| 亚洲av电影在线进入| 超色免费av| 丝袜人妻中文字幕| 亚洲熟女精品中文字幕| 少妇人妻久久综合中文| 亚洲国产欧美在线一区| 可以免费在线观看a视频的电影网站| 多毛熟女@视频| 久久久久久久精品精品| 天天躁日日躁夜夜躁夜夜| av在线播放精品| 美国免费a级毛片| 久久ye,这里只有精品| 人人妻人人添人人爽欧美一区卜| 一级毛片女人18水好多| 久久热在线av| 黄色毛片三级朝国网站| 色视频在线一区二区三区| 久久久久久久精品精品| 少妇猛男粗大的猛烈进出视频| 午夜激情av网站| 精品亚洲成a人片在线观看| 亚洲成人免费av在线播放| 亚洲成人国产一区在线观看| 男女之事视频高清在线观看| 香蕉国产在线看| 亚洲第一欧美日韩一区二区三区 | 久久久久久久国产电影| 黑人巨大精品欧美一区二区mp4| 动漫黄色视频在线观看| 日韩 欧美 亚洲 中文字幕| 侵犯人妻中文字幕一二三四区| 久久久久久人人人人人| 亚洲欧美色中文字幕在线| 日韩三级视频一区二区三区| 天堂中文最新版在线下载| 欧美日韩一级在线毛片| 人妻久久中文字幕网| 搡老乐熟女国产| 一本综合久久免费| 亚洲中文字幕日韩| 午夜福利视频在线观看免费| 91av网站免费观看| 热99re8久久精品国产| 日日爽夜夜爽网站| 老司机在亚洲福利影院| 国产深夜福利视频在线观看| 少妇被粗大的猛进出69影院| 亚洲精品国产精品久久久不卡| 精品免费久久久久久久清纯 | 美女扒开内裤让男人捅视频| 午夜久久久在线观看| √禁漫天堂资源中文www| 国产免费一区二区三区四区乱码| 黑人操中国人逼视频| 午夜福利视频精品| 黄色毛片三级朝国网站| 亚洲专区字幕在线| 在线亚洲精品国产二区图片欧美| 男女无遮挡免费网站观看| 久久久欧美国产精品| 久久中文字幕一级| 一区二区av电影网| 久久人人爽av亚洲精品天堂| a 毛片基地| 老熟妇仑乱视频hdxx| 国产主播在线观看一区二区| 久久久国产精品麻豆| 99久久国产精品久久久| netflix在线观看网站| 男女午夜视频在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 99国产精品一区二区蜜桃av | 亚洲国产毛片av蜜桃av| 日本欧美视频一区| 黄色视频不卡| 中国美女看黄片| 亚洲精品乱久久久久久| 国产亚洲精品一区二区www | 国产精品麻豆人妻色哟哟久久| 男女床上黄色一级片免费看| 国产欧美日韩综合在线一区二区| 啦啦啦在线免费观看视频4| 国产精品免费视频内射| 97在线人人人人妻| 美女高潮到喷水免费观看| 久久国产精品大桥未久av| 丝袜美腿诱惑在线| 99国产精品一区二区蜜桃av | 国产免费一区二区三区四区乱码| 精品福利观看| 一区二区三区激情视频| 免费黄频网站在线观看国产| 久久亚洲国产成人精品v| 国产亚洲精品久久久久5区| 日本av免费视频播放| 男女国产视频网站| 侵犯人妻中文字幕一二三四区| 久久人人97超碰香蕉20202| 欧美午夜高清在线| 永久免费av网站大全| 亚洲免费av在线视频| 久久久久久久精品精品| 曰老女人黄片| tube8黄色片| 国产三级黄色录像| 亚洲激情五月婷婷啪啪| www.精华液| 中文字幕精品免费在线观看视频| 老熟妇仑乱视频hdxx| 在线精品无人区一区二区三| 久久久久久久久久久久大奶| 成人国产av品久久久| 黑人巨大精品欧美一区二区蜜桃| 美女中出高潮动态图| 欧美激情高清一区二区三区| 最黄视频免费看| 久久久久精品国产欧美久久久 | 一区二区三区四区激情视频| 精品少妇一区二区三区视频日本电影| 国产野战对白在线观看| 大香蕉久久网| 亚洲人成电影观看| 丰满人妻熟妇乱又伦精品不卡| 美女国产高潮福利片在线看| 日本欧美视频一区| 成年美女黄网站色视频大全免费| 一本—道久久a久久精品蜜桃钙片| 日日爽夜夜爽网站| 国产精品久久久久久精品古装| 精品乱码久久久久久99久播| 极品人妻少妇av视频| 狠狠狠狠99中文字幕| 国产成人一区二区三区免费视频网站| 久久久久久久精品精品| 麻豆乱淫一区二区| www.av在线官网国产| 老鸭窝网址在线观看| 成年女人毛片免费观看观看9 | 在线十欧美十亚洲十日本专区| 成人国产一区最新在线观看| 午夜福利,免费看| 亚洲七黄色美女视频| 免费不卡黄色视频| 丁香六月欧美| 亚洲全国av大片| 精品久久久久久久毛片微露脸 | 亚洲色图 男人天堂 中文字幕| 正在播放国产对白刺激| 一级黄色大片毛片| 涩涩av久久男人的天堂| 日韩,欧美,国产一区二区三区| 两性夫妻黄色片| 久久久久精品国产欧美久久久 | 亚洲五月色婷婷综合| 巨乳人妻的诱惑在线观看| 亚洲精品自拍成人| 在线av久久热| 亚洲国产精品成人久久小说| 久久狼人影院| 母亲3免费完整高清在线观看| 搡老乐熟女国产| 久久国产精品影院| 国产欧美日韩一区二区三 | 国产色视频综合| 国产成人av激情在线播放| 曰老女人黄片| 国产av一区二区精品久久| 在线观看免费日韩欧美大片| 国产精品免费大片| 手机成人av网站| av线在线观看网站| 国产精品久久久久久人妻精品电影 | 成人av一区二区三区在线看 | 成人国语在线视频| 大片免费播放器 马上看| 天天影视国产精品| 最近最新中文字幕大全免费视频| 桃花免费在线播放| 国产成+人综合+亚洲专区| 欧美av亚洲av综合av国产av| 国产有黄有色有爽视频| 两人在一起打扑克的视频| 精品高清国产在线一区| 丝袜喷水一区| 久久女婷五月综合色啪小说| 国产精品秋霞免费鲁丝片| 午夜福利影视在线免费观看| 午夜福利一区二区在线看| 大香蕉久久成人网| 正在播放国产对白刺激| 久热爱精品视频在线9| 99国产精品免费福利视频| 菩萨蛮人人尽说江南好唐韦庄| 美国免费a级毛片| 久久人人97超碰香蕉20202| 亚洲性夜色夜夜综合| 亚洲人成电影免费在线| 夜夜骑夜夜射夜夜干| 新久久久久国产一级毛片| 国产精品偷伦视频观看了| 欧美久久黑人一区二区| 亚洲精品国产av蜜桃| 老司机影院毛片| 国产亚洲欧美精品永久| 亚洲精品成人av观看孕妇| 一二三四社区在线视频社区8| 热99re8久久精品国产| 不卡av一区二区三区| 午夜福利一区二区在线看| 久久久国产成人免费| 亚洲精品av麻豆狂野| 亚洲国产毛片av蜜桃av| 国产福利在线免费观看视频| 伦理电影免费视频| 99热全是精品| 成人亚洲精品一区在线观看| 十分钟在线观看高清视频www| 欧美精品亚洲一区二区| 亚洲第一青青草原| 久久久精品国产亚洲av高清涩受| 国产熟女午夜一区二区三区| 热re99久久精品国产66热6| 亚洲精品美女久久久久99蜜臀| 桃花免费在线播放| 国产老妇伦熟女老妇高清| 久久久久久久国产电影| 岛国毛片在线播放| 一进一出抽搐动态| 美女福利国产在线| 91成人精品电影| 免费女性裸体啪啪无遮挡网站| 搡老乐熟女国产| 国产欧美日韩一区二区三区在线| 最近中文字幕2019免费版| 久久九九热精品免费| 人人澡人人妻人| 精品国产超薄肉色丝袜足j| 美女高潮喷水抽搐中文字幕| 免费不卡黄色视频| 热re99久久精品国产66热6| 精品免费久久久久久久清纯 | 下体分泌物呈黄色| 久久天堂一区二区三区四区| 久久久精品国产亚洲av高清涩受| 久久久久国产精品人妻一区二区| 女人久久www免费人成看片| 国产男人的电影天堂91| 亚洲精品中文字幕在线视频| 91av网站免费观看| 欧美少妇被猛烈插入视频| 好男人电影高清在线观看| 久久人妻福利社区极品人妻图片| 看免费av毛片| 国产精品久久久av美女十八| 波多野结衣一区麻豆| 一本大道久久a久久精品| 日韩 欧美 亚洲 中文字幕| 一区二区三区精品91| 我的亚洲天堂| 精品亚洲成a人片在线观看| a在线观看视频网站| 老司机亚洲免费影院| 美女高潮到喷水免费观看| 亚洲情色 制服丝袜| 巨乳人妻的诱惑在线观看| 69精品国产乱码久久久| 老熟妇乱子伦视频在线观看 | 51午夜福利影视在线观看| 中文字幕av电影在线播放| 久久精品熟女亚洲av麻豆精品| 麻豆国产av国片精品| 亚洲成人手机| 另类亚洲欧美激情| 老司机午夜十八禁免费视频| 欧美性长视频在线观看| 成人亚洲精品一区在线观看| 亚洲情色 制服丝袜| 久久精品国产亚洲av香蕉五月 | 亚洲午夜精品一区,二区,三区| 日韩三级视频一区二区三区| 狂野欧美激情性xxxx| 999久久久精品免费观看国产| 纵有疾风起免费观看全集完整版| 国产成人a∨麻豆精品| 考比视频在线观看| 高清av免费在线| 亚洲成人免费av在线播放| 国产成人精品久久二区二区91| 一边摸一边抽搐一进一出视频| 亚洲精品自拍成人| 日韩中文字幕欧美一区二区| 黄色视频,在线免费观看| 亚洲色图综合在线观看| 日韩欧美一区视频在线观看| 国产极品粉嫩免费观看在线| 97精品久久久久久久久久精品| 国产三级黄色录像| 亚洲精品国产一区二区精华液| 亚洲av成人不卡在线观看播放网 | 99国产精品一区二区三区| 国产欧美亚洲国产| 两个人看的免费小视频| videos熟女内射| 成人黄色视频免费在线看| 欧美激情久久久久久爽电影 | 亚洲第一欧美日韩一区二区三区 | 久久午夜综合久久蜜桃| 在线观看一区二区三区激情| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久精品区二区三区| 高清黄色对白视频在线免费看| av电影中文网址| 亚洲精品国产精品久久久不卡| 成人影院久久| 日韩精品免费视频一区二区三区| 亚洲精品乱久久久久久| 久久久久久免费高清国产稀缺| 国产黄色免费在线视频| 欧美日韩精品网址| 免费观看人在逋| 亚洲精品国产av成人精品| 捣出白浆h1v1| 亚洲av电影在线观看一区二区三区| 十八禁高潮呻吟视频| 精品国产乱码久久久久久小说| 水蜜桃什么品种好| 狠狠婷婷综合久久久久久88av| 免费一级毛片在线播放高清视频 | 欧美日韩亚洲国产一区二区在线观看 | 9191精品国产免费久久| 五月开心婷婷网| av在线播放精品| 十八禁网站网址无遮挡| 国产高清国产精品国产三级| kizo精华| 免费av中文字幕在线| 久久久精品免费免费高清| 婷婷丁香在线五月| 久久人人97超碰香蕉20202| 午夜精品久久久久久毛片777| 午夜两性在线视频| 色精品久久人妻99蜜桃| 老熟妇乱子伦视频在线观看 | 男人操女人黄网站| 999精品在线视频| 别揉我奶头~嗯~啊~动态视频 | 一级片免费观看大全| 下体分泌物呈黄色| 久久久久久久久久久久大奶| 美女国产高潮福利片在线看| 天天操日日干夜夜撸| 欧美日韩精品网址| 亚洲精品中文字幕一二三四区 | 丰满迷人的少妇在线观看| 欧美黑人欧美精品刺激| 久久国产精品人妻蜜桃| 热99久久久久精品小说推荐| 免费av中文字幕在线| 午夜免费观看性视频| 免费久久久久久久精品成人欧美视频| 欧美在线黄色| 国产成人啪精品午夜网站| 精品人妻在线不人妻| 国产高清国产精品国产三级| 国产欧美日韩一区二区三 | 欧美乱码精品一区二区三区| 男女之事视频高清在线观看| 99九九在线精品视频| 国产精品自产拍在线观看55亚洲 | www.av在线官网国产| 大片免费播放器 马上看| 视频区欧美日本亚洲| 一本综合久久免费| 国产福利在线免费观看视频| 一边摸一边抽搐一进一出视频| 国产福利在线免费观看视频| 不卡av一区二区三区| 大型av网站在线播放| 国产av一区二区精品久久| 精品一区在线观看国产| 国产免费av片在线观看野外av| 国产人伦9x9x在线观看| 日本精品一区二区三区蜜桃| 一区二区三区四区激情视频| 国产av精品麻豆| 手机成人av网站| 一区二区三区精品91| 亚洲精品国产av蜜桃| 桃花免费在线播放| 巨乳人妻的诱惑在线观看| 免费av中文字幕在线| 欧美激情 高清一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 久久久久精品人妻al黑| 热re99久久国产66热|