• <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ù)新進展
    日韩 亚洲 欧美在线| 国产免费现黄频在线看| 亚洲精品乱久久久久久| 国产精品免费大片| www.自偷自拍.com| 国产在线免费精品| 美国免费a级毛片| 国产精品二区激情视频| 91精品国产国语对白视频| 男女免费视频国产| 国产欧美日韩一区二区三区在线| 丝袜人妻中文字幕| 成人国产一区最新在线观看 | 99国产精品一区二区三区| 在线看a的网站| 少妇 在线观看| 欧美日韩成人在线一区二区| 午夜精品国产一区二区电影| 久久久精品国产亚洲av高清涩受| 国产主播在线观看一区二区 | 久久久久网色| 超色免费av| 国产欧美日韩一区二区三区在线| av天堂久久9| 国产成人精品无人区| 欧美日韩黄片免| 日韩一卡2卡3卡4卡2021年| 久久久久久久大尺度免费视频| 中文欧美无线码| 国产免费又黄又爽又色| 一本色道久久久久久精品综合| 青春草视频在线免费观看| 免费在线观看日本一区| 国产福利在线免费观看视频| 精品国产超薄肉色丝袜足j| 可以免费在线观看a视频的电影网站| 狠狠精品人妻久久久久久综合| 一级片'在线观看视频| 亚洲中文av在线| 欧美精品亚洲一区二区| 成人国语在线视频| 国产不卡av网站在线观看| 欧美 亚洲 国产 日韩一| 大香蕉久久成人网| 久久人人97超碰香蕉20202| 日本猛色少妇xxxxx猛交久久| 亚洲av电影在线进入| 婷婷色av中文字幕| 亚洲情色 制服丝袜| 久久免费观看电影| 国产野战对白在线观看| 黑人欧美特级aaaaaa片| 人成视频在线观看免费观看| 乱人伦中国视频| 欧美日韩视频高清一区二区三区二| 狠狠婷婷综合久久久久久88av| 亚洲av成人不卡在线观看播放网 | 交换朋友夫妻互换小说| 精品国产一区二区三区久久久樱花| 女人爽到高潮嗷嗷叫在线视频| 91国产中文字幕| 成人亚洲精品一区在线观看| 午夜福利乱码中文字幕| 国产无遮挡羞羞视频在线观看| 无限看片的www在线观看| 国产日韩一区二区三区精品不卡| 久久久久久亚洲精品国产蜜桃av| 一本一本久久a久久精品综合妖精| 中文乱码字字幕精品一区二区三区| 国产精品成人在线| 午夜福利乱码中文字幕| 国产精品一区二区精品视频观看| 大话2 男鬼变身卡| 丝袜在线中文字幕| 国产女主播在线喷水免费视频网站| 青草久久国产| 婷婷色av中文字幕| 黑人欧美特级aaaaaa片| 亚洲人成电影观看| 国产精品国产av在线观看| av在线播放精品| 在线观看免费午夜福利视频| 免费看不卡的av| 日韩免费高清中文字幕av| 在线精品无人区一区二区三| 国产av国产精品国产| 天天躁夜夜躁狠狠久久av| av天堂在线播放| 亚洲av电影在线观看一区二区三区| 人体艺术视频欧美日本| 男人爽女人下面视频在线观看| 日本wwww免费看| 国产视频一区二区在线看| 精品国产一区二区久久| 日韩,欧美,国产一区二区三区| 脱女人内裤的视频| 国产成人精品久久二区二区91| 国产xxxxx性猛交| 精品一区二区三区av网在线观看 | 波野结衣二区三区在线| 超碰97精品在线观看| 欧美国产精品一级二级三级| 欧美激情极品国产一区二区三区| 国精品久久久久久国模美| 亚洲色图 男人天堂 中文字幕| a级片在线免费高清观看视频| 精品久久久精品久久久| 亚洲伊人久久精品综合| kizo精华| 色播在线永久视频| 久久性视频一级片| 国产欧美日韩精品亚洲av| svipshipincom国产片| 日韩熟女老妇一区二区性免费视频| 另类精品久久| 午夜久久久在线观看| 老司机影院毛片| 亚洲 国产 在线| 国产日韩一区二区三区精品不卡| 中文欧美无线码| 亚洲第一青青草原| 亚洲熟女毛片儿| 久久青草综合色| 99热网站在线观看| 日韩熟女老妇一区二区性免费视频| 精品亚洲乱码少妇综合久久| 婷婷色综合www| 亚洲国产中文字幕在线视频| 日韩av不卡免费在线播放| 尾随美女入室| 一本大道久久a久久精品| 国产av国产精品国产| 这个男人来自地球电影免费观看| 女人被躁到高潮嗷嗷叫费观| 天堂8中文在线网| 视频区图区小说| 色94色欧美一区二区| 久久久久网色| av在线播放精品| 操出白浆在线播放| 人人妻人人澡人人爽人人夜夜| 国产一区有黄有色的免费视频| 欧美日韩国产mv在线观看视频| 在线av久久热| 老鸭窝网址在线观看| www.999成人在线观看| 午夜免费成人在线视频| a 毛片基地| 久久精品久久精品一区二区三区| 亚洲 欧美一区二区三区| 又粗又硬又长又爽又黄的视频| 纵有疾风起免费观看全集完整版| 91成人精品电影| 欧美精品亚洲一区二区| 50天的宝宝边吃奶边哭怎么回事| 少妇被粗大的猛进出69影院| 亚洲,一卡二卡三卡| 亚洲av男天堂| 日本一区二区免费在线视频| 日韩伦理黄色片| bbb黄色大片| 中文字幕人妻熟女乱码| 18在线观看网站| 伊人亚洲综合成人网| 成在线人永久免费视频| 精品国产一区二区三区久久久樱花| 欧美日本中文国产一区发布| 在线精品无人区一区二区三| 中文字幕人妻丝袜一区二区| 国产伦人伦偷精品视频| 国产成人av教育| 国产男人的电影天堂91| 国产女主播在线喷水免费视频网站| 桃花免费在线播放| 人妻人人澡人人爽人人| 精品亚洲乱码少妇综合久久| 午夜福利免费观看在线| 欧美亚洲 丝袜 人妻 在线| 国产人伦9x9x在线观看| 丰满饥渴人妻一区二区三| 久久天堂一区二区三区四区| 免费女性裸体啪啪无遮挡网站| 成人国语在线视频| av电影中文网址| 午夜福利视频在线观看免费| 国产女主播在线喷水免费视频网站| 国产精品欧美亚洲77777| 一边摸一边抽搐一进一出视频| 久久久精品免费免费高清| 18在线观看网站| 亚洲欧美日韩高清在线视频 | 性高湖久久久久久久久免费观看| 国产成人一区二区三区免费视频网站 | 亚洲av日韩在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲自偷自拍图片 自拍| 丝袜在线中文字幕| 黄片小视频在线播放| 亚洲人成77777在线视频| 下体分泌物呈黄色| 亚洲三区欧美一区| 在线亚洲精品国产二区图片欧美| 亚洲 国产 在线| 国产深夜福利视频在线观看| 我的亚洲天堂| 男女边摸边吃奶| 高清视频免费观看一区二区| 欧美亚洲 丝袜 人妻 在线| 欧美精品一区二区大全| 国产av精品麻豆| 黄频高清免费视频| 久久精品成人免费网站| 别揉我奶头~嗯~啊~动态视频 | 成年人免费黄色播放视频| 成人免费观看视频高清| 欧美日韩综合久久久久久| 成人手机av| 日韩av免费高清视频| 一区二区三区乱码不卡18| 99九九在线精品视频| 热99国产精品久久久久久7| 亚洲欧美日韩另类电影网站| 久久精品久久久久久久性| 亚洲欧美精品综合一区二区三区| 欧美 亚洲 国产 日韩一| 国产成人a∨麻豆精品| av电影中文网址| 91九色精品人成在线观看| www日本在线高清视频| 好男人视频免费观看在线| 亚洲国产最新在线播放| 精品高清国产在线一区| 99九九在线精品视频| 亚洲国产最新在线播放| 男人爽女人下面视频在线观看| 午夜激情av网站| 日韩中文字幕欧美一区二区 | 免费在线观看日本一区| 免费观看a级毛片全部| 亚洲一卡2卡3卡4卡5卡精品中文| 精品久久久久久电影网| xxx大片免费视频| 久久午夜综合久久蜜桃| 91精品伊人久久大香线蕉| 2021少妇久久久久久久久久久| 亚洲专区国产一区二区| 另类精品久久| 亚洲伊人色综图| 青春草视频在线免费观看| 男女床上黄色一级片免费看| 五月开心婷婷网| 国产亚洲欧美在线一区二区| 国产精品一区二区在线观看99| 亚洲一码二码三码区别大吗| 后天国语完整版免费观看| 首页视频小说图片口味搜索 | 精品人妻一区二区三区麻豆| 人人妻人人澡人人爽人人夜夜| 欧美精品一区二区大全| 十八禁人妻一区二区| 少妇人妻 视频| 人人妻人人澡人人看| 久久久国产一区二区| 在线观看国产h片| 亚洲中文字幕日韩| 日本五十路高清| 两个人免费观看高清视频| 国产亚洲欧美在线一区二区| 交换朋友夫妻互换小说| 如日韩欧美国产精品一区二区三区| e午夜精品久久久久久久| 亚洲免费av在线视频| 欧美日韩亚洲国产一区二区在线观看 | 色网站视频免费| 日韩人妻精品一区2区三区| 国产精品成人在线| 亚洲免费av在线视频| 久久久久国产一级毛片高清牌| 搡老乐熟女国产| 欧美乱码精品一区二区三区| 国产无遮挡羞羞视频在线观看| av在线老鸭窝| www.av在线官网国产| 亚洲一区二区三区欧美精品| 免费不卡黄色视频| 国产精品一区二区免费欧美 | 精品第一国产精品| 高清欧美精品videossex| 悠悠久久av| 精品人妻在线不人妻| 香蕉丝袜av| 中文字幕人妻熟女乱码| 黑人猛操日本美女一级片| 国产免费一区二区三区四区乱码| 又紧又爽又黄一区二区| 十八禁高潮呻吟视频| 91麻豆精品激情在线观看国产 | 国产片内射在线| 王馨瑶露胸无遮挡在线观看| 一区二区三区激情视频| av在线老鸭窝| 啦啦啦 在线观看视频| 日本a在线网址| a级毛片在线看网站| 日韩大片免费观看网站| 青青草视频在线视频观看| 国产在视频线精品| 精品福利观看| 一个人免费看片子| 亚洲久久久国产精品| 国产91精品成人一区二区三区 | 男女高潮啪啪啪动态图| 黄色怎么调成土黄色| 久久青草综合色| 老司机亚洲免费影院| 一个人免费看片子| 老司机影院毛片| 久久免费观看电影| 桃花免费在线播放| 亚洲成人手机| 亚洲午夜精品一区,二区,三区| 97人妻天天添夜夜摸| www.自偷自拍.com| 日韩电影二区| 亚洲国产精品一区二区三区在线| 男女免费视频国产| 大话2 男鬼变身卡| 少妇人妻 视频| 丰满迷人的少妇在线观看| 麻豆乱淫一区二区| 青春草视频在线免费观看| 日韩,欧美,国产一区二区三区| 欧美精品人与动牲交sv欧美| 91麻豆av在线| 久久精品国产综合久久久| 国产成人欧美| 99re6热这里在线精品视频| 久久精品成人免费网站| 91麻豆av在线| 久久久久视频综合| 久久国产精品人妻蜜桃| 国产精品99久久99久久久不卡| 免费少妇av软件| 90打野战视频偷拍视频| 香蕉国产在线看| 成人影院久久| 久久国产精品影院| 少妇精品久久久久久久| 成人亚洲精品一区在线观看| 国产精品成人在线| 成人18禁高潮啪啪吃奶动态图| 亚洲一区中文字幕在线| 亚洲第一av免费看| 视频区欧美日本亚洲| 美女主播在线视频| 亚洲色图综合在线观看| 国产一级毛片在线| 蜜桃在线观看..| 男女边摸边吃奶| 欧美大码av| 一区二区三区激情视频| 欧美黑人欧美精品刺激| 大香蕉久久成人网| 色视频在线一区二区三区| 久久精品久久久久久久性| 尾随美女入室| 免费看十八禁软件| 免费av中文字幕在线| 国产精品一区二区精品视频观看| av欧美777| 亚洲一码二码三码区别大吗| av又黄又爽大尺度在线免费看| 老司机深夜福利视频在线观看 | 一边摸一边做爽爽视频免费| 久久精品熟女亚洲av麻豆精品| 满18在线观看网站| 2018国产大陆天天弄谢| 欧美日韩av久久| av在线老鸭窝| 国产伦人伦偷精品视频| 精品欧美一区二区三区在线| 欧美人与性动交α欧美软件| 亚洲欧美色中文字幕在线| 天堂中文最新版在线下载| 国产成人欧美在线观看 | 精品一区在线观看国产| 国产日韩欧美视频二区| 欧美激情 高清一区二区三区| 国产精品成人在线| 91老司机精品| 麻豆乱淫一区二区| 高清视频免费观看一区二区| 麻豆国产av国片精品| 久久99热这里只频精品6学生| 热re99久久精品国产66热6| 国产91精品成人一区二区三区 | 女性被躁到高潮视频| 久热这里只有精品99| 国产黄色免费在线视频| 国产高清国产精品国产三级| 国产精品一二三区在线看| 国产91精品成人一区二区三区 | 不卡av一区二区三区| 在线观看人妻少妇| 黄频高清免费视频| 大片电影免费在线观看免费| 777久久人妻少妇嫩草av网站| 丝袜喷水一区| 国产一区二区激情短视频 | 看十八女毛片水多多多| 超色免费av| 久久影院123| 老司机靠b影院| 女人久久www免费人成看片| 日韩伦理黄色片| 99久久综合免费| 国产欧美日韩一区二区三区在线| 精品少妇黑人巨大在线播放| 日本wwww免费看| 女性被躁到高潮视频| 精品亚洲乱码少妇综合久久| 久久国产精品大桥未久av| 日韩大片免费观看网站| 成年人午夜在线观看视频| 在线观看免费视频网站a站| 亚洲人成电影观看| 9色porny在线观看| 悠悠久久av| 精品国产一区二区三区久久久樱花| 欧美国产精品一级二级三级| 亚洲熟女毛片儿| 老汉色av国产亚洲站长工具| 精品久久久久久电影网| 女人爽到高潮嗷嗷叫在线视频| 在线观看一区二区三区激情| 久久99精品国语久久久| 国产av国产精品国产| 婷婷色av中文字幕| 国产免费现黄频在线看| 爱豆传媒免费全集在线观看| 亚洲国产精品成人久久小说| 亚洲精品一二三| 欧美人与性动交α欧美精品济南到| 国产成人av激情在线播放| 亚洲精品国产av蜜桃| 成人亚洲欧美一区二区av| 日韩一卡2卡3卡4卡2021年| 性色av一级| 侵犯人妻中文字幕一二三四区| 丝袜喷水一区| 亚洲精品美女久久av网站| 国产精品成人在线| 国产视频首页在线观看| 一区二区三区乱码不卡18| 一本久久精品| 2021少妇久久久久久久久久久| 精品一区二区三区四区五区乱码 | 精品熟女少妇八av免费久了| 精品第一国产精品| 国产精品国产三级国产专区5o| 深夜精品福利| 久9热在线精品视频| 99国产精品一区二区三区| 亚洲,欧美,日韩| 精品国产一区二区久久| a级毛片黄视频| 在线天堂中文资源库| 啦啦啦在线观看免费高清www| 秋霞在线观看毛片| 亚洲av电影在线进入| av一本久久久久| 成人国产一区最新在线观看 | 真人做人爱边吃奶动态| 欧美av亚洲av综合av国产av| 亚洲五月婷婷丁香| 免费在线观看视频国产中文字幕亚洲 | 一级毛片我不卡| av视频免费观看在线观看| 交换朋友夫妻互换小说| 啦啦啦在线观看免费高清www| 亚洲欧美日韩另类电影网站| 国产精品一区二区在线观看99| 成人三级做爰电影| 悠悠久久av| 母亲3免费完整高清在线观看| 啦啦啦中文免费视频观看日本| 国产xxxxx性猛交| 婷婷色综合大香蕉| 欧美人与性动交α欧美软件| 男女床上黄色一级片免费看| 国产一区二区在线观看av| 丝袜美足系列| 免费看不卡的av| 操出白浆在线播放| 一边摸一边抽搐一进一出视频| 欧美黑人欧美精品刺激| 国产精品亚洲av一区麻豆| 蜜桃国产av成人99| 在线看a的网站| 亚洲精品久久午夜乱码| 亚洲国产毛片av蜜桃av| 99香蕉大伊视频| 国产黄色视频一区二区在线观看| 久久女婷五月综合色啪小说| 欧美大码av| 国产又爽黄色视频| 中文欧美无线码| 国产精品久久久久久精品电影小说| av有码第一页| 亚洲欧洲精品一区二区精品久久久| 91成人精品电影| 国产黄色免费在线视频| 国语对白做爰xxxⅹ性视频网站| 日本黄色日本黄色录像| 欧美精品av麻豆av| 久久精品熟女亚洲av麻豆精品| 99热网站在线观看| 久久久精品区二区三区| 一级毛片黄色毛片免费观看视频| 9色porny在线观看| 日韩电影二区| 国产成人免费观看mmmm| 水蜜桃什么品种好| 人人妻人人爽人人添夜夜欢视频| 午夜老司机福利片| 一级黄色大片毛片| 久久精品国产亚洲av高清一级| 中文字幕高清在线视频| 精品少妇黑人巨大在线播放| 国产精品欧美亚洲77777| 精品亚洲乱码少妇综合久久| 免费观看a级毛片全部| 国产精品99久久99久久久不卡| 久久人妻福利社区极品人妻图片 | 91字幕亚洲| 一区在线观看完整版| 午夜两性在线视频| 亚洲av日韩在线播放| 国产免费一区二区三区四区乱码| 亚洲精品美女久久久久99蜜臀 | 欧美在线黄色| 日韩大片免费观看网站| 天天躁夜夜躁狠狠躁躁| 制服诱惑二区| 婷婷成人精品国产| 日韩大码丰满熟妇| 777久久人妻少妇嫩草av网站| 欧美xxⅹ黑人| 国产在视频线精品| 国产在线免费精品| av国产精品久久久久影院| 亚洲国产精品999| 搡老乐熟女国产| 国产主播在线观看一区二区 | 国产高清视频在线播放一区 | 男女免费视频国产| 最近手机中文字幕大全| 久久精品人人爽人人爽视色| 亚洲 欧美一区二区三区| 18禁国产床啪视频网站| svipshipincom国产片| 赤兔流量卡办理| 欧美大码av| 久久av网站| 男女国产视频网站| 久久精品aⅴ一区二区三区四区| 在线av久久热| 国产亚洲精品第一综合不卡| 视频区图区小说| 亚洲 国产 在线| 中文字幕人妻熟女乱码| 免费看av在线观看网站| 国产精品国产三级专区第一集| 亚洲欧美一区二区三区国产| 久久99一区二区三区| 久久精品亚洲熟妇少妇任你| 国产免费福利视频在线观看| 高清欧美精品videossex| 精品国产乱码久久久久久男人| 好男人视频免费观看在线| 久久女婷五月综合色啪小说| 国产免费一区二区三区四区乱码| 久久久国产精品麻豆| 国产成人a∨麻豆精品| 看十八女毛片水多多多| 在线观看免费高清a一片| 亚洲精品美女久久久久99蜜臀 | 国产真人三级小视频在线观看| 国产成人欧美在线观看 | 亚洲 国产 在线| 视频区图区小说| 亚洲国产精品国产精品| 国产精品秋霞免费鲁丝片| 精品人妻1区二区| 国产成人av教育| 啦啦啦啦在线视频资源| 高清黄色对白视频在线免费看| 国产爽快片一区二区三区| 岛国毛片在线播放| 美女福利国产在线| 久久国产精品大桥未久av| 亚洲精品国产色婷婷电影| av视频免费观看在线观看| 久久久欧美国产精品| svipshipincom国产片| 精品人妻在线不人妻| 视频区欧美日本亚洲| 国产男女内射视频| 老司机靠b影院| av网站免费在线观看视频| 777久久人妻少妇嫩草av网站| 美女大奶头黄色视频| 国产免费视频播放在线视频| 视频区欧美日本亚洲| 精品国产国语对白av|