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

    雙平方根方程偏移反假頻研究

    2023-02-09 12:42:50劉守偉周靜毅
    石油物探 2023年1期
    關鍵詞:波場平方根間距

    劉守偉,楊 浩,李 列,隋 波,周靜毅

    (1.南方海洋科學與工程廣東省實驗室(湛江),廣東湛江524088;2.上海青鳳致遠地球物理地質勘探科技有限公司,上海200093;3.中海石油(中國)有限公司湛江分公司,廣東湛江524088;4.中石化上海海洋油氣分公司勘探開發(fā)研究院,上海200120)

    地震數據采集觀測系統(tǒng)設計的重要目標之一是平衡成本和控制數據假頻。即使是近年來發(fā)展起來的“兩寬一高”地震數據采集,對于中高頻數據空間采樣仍舊不夠,假頻現象仍然不可忽視[1]。假頻的存在給地震數據去噪、偏移和反演都帶來了很大困難。在偏移成像中,地球物理學家分別對Kirchhoff積分偏移、單程波偏移、逆時偏移和高斯束偏移等成像技術的反假頻進行了系列研究,改善了相關算法的成像效果[2-7]。假頻通常劃分為數據空間假頻、算子假頻和成像條件假頻[4]。不同的成像方法面對的假頻也不相同,如Kirchhoff積分偏移通常要處理算子假頻[5-6]、炮域單程波方程偏移和逆時偏移通常要處理數據假頻和成像條件假頻等[2,4]。當前,幾乎所有的工業(yè)偏移成像方法都會對假頻進行處理[3]。

    雙平方根方程偏移由于其計算效率和精度較高、偏移孔徑無需人為選擇、偏移角度道集容易輸出等優(yōu)點,成為一種重要的成像方法[8-10]。近年來,業(yè)界研究人員持續(xù)研究該成像方法的各向異性介質和吸收衰減介質的推廣、角度道集生成、傳播精度改進以及實際數據的應用[11-16]。對于雙平方根方程偏移反假頻研究,ZHANG等[3]指出需要以半道間距成像來避免假頻。在偏移實踐中,通常我們在中點-半偏移距域進行外推及成像[8-12],中點方向采樣間隔為道距的一半,看上去這種雙平方根偏移不受假頻影響。但是,在實際偏移應用中,需要同時滿足中點-半偏移距兩個維度(三維地震數據為四個維度,分別為Inline方向和Crossline方向中點、Inline和Crossline方向半偏移距)均沒有假頻。ZHANG等[3]僅指出以半道距成像還遠遠不夠。對于雙平方根方程偏移,滿足中點方向以半道間距采樣情況下沒有成像條件假頻,但對半偏移距方向的數據假頻并沒有展開論述,這才是產生雙平方根方程偏移假頻的關鍵因素,只保證波場以半道間距間隔傳播和成像無法完全消除該算法的假頻。對于當前三維實際數據,炮線間距和檢波線間距較大,面元間距通常選擇炮點和檢波點間距的一半,偏移距間距為炮線和檢波線間距的兩倍,就造成一個面元中Inline方向和Crossline方向偏移距間距很大,導致在中點-半偏移距域數據中假頻現象非常嚴重。如果半偏移距方向的假頻處理不到位,雙平方根方程波場延拓會引入大量噪聲,嚴重影響成像質量。這也是該方法多年來沒有被工業(yè)界大規(guī)模應用的重要原因這一。因此,開展雙平方根方程偏移反假頻的研究及消除十分必要。當前,針對雙平方根方程偏移的假頻研究除了在文獻[3]中簡單提及,并沒有引起地球物理業(yè)界的廣泛關注,針對性的假頻研究文獻還較為少見。

    本文首先簡單回顧了雙平方根方程偏移成像算法的基本原理,然后從采樣定理角度分析了成像過程中假頻的來源,給出了解決中點-半偏移距域數據空間假頻問題的三角平滑濾波方法,指出了該方法相較數據插值加密的優(yōu)點。最后,利用Marmousi模型數據驗證了本文方法的有效性,并對比了與常規(guī)方法計算需求的差異,進一步利用實際數據對本文方法進行了驗證。

    1 雙平方根偏移

    為了論述需要,首先,簡單回顧雙平方根方程偏移方法。在橫向變速情況下,雙平方根波場外推方程為[8]:

    (1)

    式中:Xs=(xs,ys),Xg=(xg,yg)分別為炮點和檢波點坐標;P(t,Xs,Xg,z)為z深度t時刻波場;v(Xs,z)和v(Xg,z)分別表示炮、檢點層速度。ZHANG等[3]指出該方程無成像假頻的條件是需要滿足成像間隔為檢波點間距的一半。對于雙平方根方程偏移,直接在中點-半偏移距域實現更為自然。轉化方程(1)到中點-半偏移距域,并將波場表示為頻率-波數域:

    (2)

    式中:ω表示圓頻率;M=(Xs+Xg)/2和H=(Xg-Xs)/2分別表示中點和半偏移距;FT表示傅里葉變換;Δs為炮點方向坐標增量;k0=ω/v(z),kz,Km=(kmx,kmy),Kh=(khx,khy)分別表示總波數、深度方向波數、中心點方向波數和半偏移距方向波數,且有:

    (3)

    顯然炮點和檢波點方向波數滿足:

    (4)

    式中:Δg為檢波點方向坐標增量。

    對(2)式進行分解,得:

    (5)

    (6)

    (5)式是背景介質中背景波場的相移外推方程。(6)式描述了慢度攝動引起的散射場。背景場和散射場的和構成總場。對(6)式進行不同的近似,得到不同精度的外推方程。

    (5)式和(6)式的雙平方根算子將激發(fā)點和接收點同時向下外推,并根據零時間零偏移距成像條件:

    (7)

    實現成像值提取,完成疊前深度偏移[8]。

    2 數據空間假頻的分析及消除

    (8)

    式中:fmax為數據最高頻率。即:

    (9)

    式中:ph=?t/?h為半偏移距方向的射線參數。同理,中點方向采樣間隔Δm需要滿足:

    (10)

    式中:pm為中點方向的射線參數。

    (9)式和(10)式限定了滿足中點-半偏移距域雙平方根延拓的空間采樣間隔。對于橫向緩變速介質,在CMP道集中,中點-半偏移距域時距關系近似滿足:

    (11)

    式中:vrms為均方根速度;t0為自激自收時間;t為地震波旅行時間。將(11)式對半偏移距h求偏導數并代入(9)式,得:

    (12)

    (12)式為半偏移距方向無假頻必須滿足的條件,即時間越小、均方根速度越小、頻率越高、偏移距越大,要求采樣間距越小。即,淺層、低速區(qū)域和大偏移距數據容易產生假頻,偏移距間距越大越容易產生假頻。

    在二維共偏移距道集或者三維OVT道集中,反射數據中點方向射線參數滿足:

    (13)

    式中:α為反射界面的時間域傾角。通常,傾角α較小,因此,(13)式較容易滿足。在中點-半偏移距域波場延拓中,中點間距為道距一半,因此,(13)式的采樣條件幾乎總是成立的。分析當前主流的觀測數據CMP道集,在中點間距為道距一半的情況下,對于二維規(guī)則觀測數據,半偏移距間距恰好等于炮間距,它通常是道間距的兩倍或者三倍;對于三維正交規(guī)則觀測數據,Inline方向和Crossline方向的半偏移距間距分別為炮線間距和檢波線間距,通常為數百米的量級。即現行主流觀測系統(tǒng),二維數據的炮間距和三維數據的炮線間距/檢波線間距都較大,對于中高頻數據非常容易違反公式(12)的采樣條件,從而產生假頻。

    消除數據空間假頻通常有兩種方法,方法一為壓制數據的高頻信息,從而使數據滿足無假頻采樣條件;方法二為進行數據反假頻插值,加密數據的空間采樣,從而壓制假頻。方法一需要采取局部濾波措施,避免全局濾波器濾波造成波場的過度低頻化,本文采用方法一進行高頻壓制。方法二會加大數據量,從而增加算法的計算量和內存需求量,影響算法的效率。

    LUMLEY等[5]根據三角平滑濾波器的局部低通性質通過選擇恰當的濾波器長度進行低通平滑濾波,壓制高頻數據,實現無假頻的Kirchhoff偏移。結合三角平滑濾波器性質和(12)式CMP域半偏移距采樣條件,設計合適的濾波器,實現雙平方根方程偏移的反假頻處理。根據三角平滑濾波器性質,濾波器長度N需要滿足:

    (14)

    該濾波器在對地震道數據進行疊加的基礎上,通過三點運算完成,計算效率非常高[5]。在雙平方根方程偏移波場延拓之前根據公式(14)設計的濾波器對CMP輸入數據進行三角平滑濾波,壓制中、淺層低速區(qū)域數據假頻,從而實現無假頻波場外推。

    圖1為Marmousi模型某CMP道集反假頻前后對比結果,圖1a的半偏移距間距為50m,可以看到淺層大偏移距位置處明顯存在假頻。圖1b為根據(14)式設計的三角平滑濾波算子長度濾波后的數據。淺層大偏移距數據高頻得到壓制,小偏移距和中深層沒有假頻的同相軸都得到了很好的保持。

    圖1 CMP道集數據反假頻前(a)、后(b)的對比結果

    3 模型數據及實際資料測試

    3.1 Marmousi模型數據

    Marmousi模型數據共計240炮,記錄長度為3s,時間采樣間隔為4ms。單邊接收,每炮96道接收,炮間距為50m,道間距為25m,最大偏移距為2575m。成像中心點間距為12.5m,半偏移距間距為50m。對于該模型數據,最高頻率為40Hz,淺層0.5s均方根速度為1750m/s。對于該時間,根據(12)式,半偏移距需要滿足Δh≤4785/h,即對于半偏移距為250.0m的地震道,需要滿足半偏移距間距小于19.1m。以半偏移距50m為間距進行成像對淺層顯然會出現明顯假頻。對于中深層2.1s均方根速度約為2700m/s,需要滿足Δh≤47840/h,即對于最大半偏移距1287.5m的地震道,半偏移距需要小于37.0m,對于最大半偏移距1000m的地震道,半偏移距間距小于47.8m即可。

    圖2a為平滑后的Marmousi速度模型,圖2b為利用成像點間距50m且沒有進行反假頻處理雙平方根方程偏移剖面,圖2c為利用本文方法進行反假頻處理后的雙平方根方程偏移剖面。可以看出,圖2b存在明顯的假頻噪聲干擾,尤其是中淺層,假頻非常明顯,隨著深度增加,假頻逐步減少。圖2c中已經基本看不到假頻噪聲的影響。圖2d為將偏移距間距加密到12.5m后的成像結果,可以看到,加密采樣后也能夠消除假頻的影響,改善成像質量。圖2e至圖2g分別為圖2b至圖2d右上角放大后的顯示,可以更明顯地看到假頻噪聲對成像效果的影響和反假頻后成像效果的改善。

    圖2 Marmousi模型雙平方根方程偏移反假頻效果對比a 速度模型; b 進行反假頻前成像剖面; c 進行反假頻后剖面; d 加密偏移距方向采樣后剖面; e b圖右上角放大顯示; f c圖右上角放大顯示; g d圖右上角放大顯示

    根據(8)式和圖2成像結果的對比可知,本文方法和加密采樣均能夠很好地實現假頻噪聲的壓制,但是加密空間采樣會加大計算內存需求和計算量,減弱了雙平方根成像方法計算高效的優(yōu)勢。表1為圖2c和圖2d計算內存需求和計算時間統(tǒng)計。內存需求方面,常規(guī)加密采樣方法將偏移距加密為原來的4倍,整個數據規(guī)模也增加為原來的4倍,需求的內存隨之膨脹為原來的4倍。計算量方面,隨著內存需求的增加,計算量也增大為原來的4倍。

    表1 計算需求對比

    結合表1計算需求數據以及圖2成像結果對比可知,加密采樣方法和本文方法在成像上均可以有效壓制假頻,但在存儲和計算時間上加密采樣方法均隨著加密采樣而數據規(guī)模顯著增加,本文方法是兼顧效率和效果的方法。

    3.2 實際資料

    某三維地震數據Inline方向以及Crossline方向半偏移距間隔250m,道間距為50m,中心點間距為25m,數據最高頻率60Hz。分別利用本文反假頻方法處理前、后的雙平方根方程偏移成像的結果如圖3和圖4 所示,其中,圖3a和圖4a為未進行反假頻處理的雙平方根方程偏移Crossline方向和Inline方向剖面,圖3b和圖4b為采用本文方法進行反假頻處理后的Crossline方向和Inline方向成像剖面。對比可以看出,反假頻處理之前,成像剖面中淺層假頻噪聲嚴重,有效反射波同相軸被假頻覆蓋。利用本文方法進行反假頻壓制,在剖面1km上下范圍內的假頻噪聲得到了有效壓制,有效信號得到凸顯。對于該數據,常規(guī)加密采樣方法需要加密半偏移距間距到25m,數據規(guī)模是原來的100倍,計算內存和計算量也膨脹100倍,已經遠遠超過了測試計算機的內存,因此加密采樣對于三維數據不具備可行性。

    圖3 實際三維地震數據Crossline方向雙平方根方程偏移反假頻前、后效果對比a 進行反假頻前成像剖面; b 進行反假頻后剖面

    圖4 實際三維地震數據Inline方向雙平方根方程偏移反假頻前、后效果對比a 進行反假頻前成像剖面; b 進行反假頻后剖面

    4 討論

    本文的壓制假頻方法以地下緩變介質的時距曲線為基礎,對于絕大多數介質是適用的。對于復雜介質,理論上只有準確把握同相軸的曲率才能設計準確的反假頻公式,在實際處理中既不可能做到,也沒有實際意義。在實際應用過程中,可以在假頻公式中加入調節(jié)系數,根據地下構造傾角反假頻濾波器長度進行調節(jié),實現不同傾角構造的假頻壓制,提高方法對復雜構造的適應性。

    本文方法的優(yōu)勢在于能夠很好地平衡假頻的壓制和計算需求,保障偏移的效率。眾所周知,對于假頻的壓制,最直接的方式是通過反假頻插值加密空間采樣密度。但是加密空間采樣帶來的問題是計算需求的急劇提高。對于三維雙平方根方程偏移,需要同時在中點-半偏移距域進行波場延拓,這里的中點方向和半偏移距方向分別為一個兩維空間,兩個方向組成一個四維數組進行波場延拓。這個四維數組需要較大的計算機內存,如果通過反假頻插值加密半偏移距方向的采樣,數據量會急劇膨脹,造成偏移算法無法在當前主流的計算機上進行運算或者計算效率低下。結合觀測系統(tǒng)實際情況,適度地加密采樣并利用本文反假頻方法,減小對淺層大偏移距高頻數據的壓制,得到計算量和計算效果均衡的結果。

    5 結論

    偏移假頻通常分為數據假頻、算子假頻和成像條件假頻,對于中點-半偏移距域雙平方根方程偏移本身不會產生傳播算子及成像條件假頻,數據空間假頻是該成像方法假頻的唯一來源。對于地震數據,偏移距方向的射線參數與時距曲線相關,中點方向的射線參數與界面時間域傾角相關,當前主流的地震觀測系統(tǒng),半偏移距采樣間距通常是中點間距的數倍,中點方向不容易產生假頻,主要假頻來源為半偏移距方向。

    本文提出了針對中點-半偏移距域雙平方根方程偏移數據空間假頻消除方法。該方法在偏移成像之前對地震數據進行三角平滑濾波器濾波,壓制偏移距方向的數據假頻,從而避免波場延拓過程中偏移距方向的假頻輸入,實現無假頻傳播及成像,提高成像品質。該方法計算效率非常高,與偏移成像相比,計算量幾乎可以忽略不計。模型數據和實際資料測試表明,利用本文方法進行反假頻處理后的成像剖面假頻噪聲得到了有效壓制,明顯改善了地震數據的成像質量。

    猜你喜歡
    波場平方根間距
    宋代君臣殿上間距考論
    高速公路指揮中心小間距LED應用探討
    彈性波波場分離方法對比及其在逆時偏移成像中的應用
    平方根易錯點警示
    “平方根”學習法升級版
    如何學好平方根
    幫你學習平方根
    交錯網格與旋轉交錯網格對VTI介質波場分離的影響分析
    地震學報(2016年1期)2016-11-28 05:38:36
    基于Hilbert變換的全波場分離逆時偏移成像
    算距離
    亚洲第一欧美日韩一区二区三区| 日本黄大片高清| 草草在线视频免费看| 亚洲欧美日韩东京热| 精品人妻一区二区三区麻豆 | 一本一本综合久久| 欧美性猛交黑人性爽| 丰满的人妻完整版| 国产野战对白在线观看| 欧美精品啪啪一区二区三区| www.熟女人妻精品国产| 99精品久久久久人妻精品| 淫秽高清视频在线观看| 日韩成人在线观看一区二区三区| 国产免费男女视频| 久久久久久久亚洲中文字幕 | 精品日产1卡2卡| 搡老妇女老女人老熟妇| 亚洲美女视频黄频| 久久国产精品人妻蜜桃| 日韩有码中文字幕| 日韩免费av在线播放| 真人做人爱边吃奶动态| 非洲黑人性xxxx精品又粗又长| 亚洲欧美日韩无卡精品| 国产精品人妻久久久久久| 久久午夜福利片| 老熟妇仑乱视频hdxx| 成人亚洲精品av一区二区| 丰满人妻一区二区三区视频av| 国产熟女xx| 久久久久国内视频| 日本一本二区三区精品| 无遮挡黄片免费观看| 国产欧美日韩精品一区二区| 青草久久国产| 少妇丰满av| 免费在线观看日本一区| 国产精品久久久久久久久免 | 亚洲男人的天堂狠狠| 黄色视频,在线免费观看| 人人妻,人人澡人人爽秒播| 国产亚洲精品av在线| 热99re8久久精品国产| 国产三级黄色录像| 97碰自拍视频| 一本综合久久免费| 精品午夜福利视频在线观看一区| 国产精品女同一区二区软件 | 国产三级在线视频| 亚洲av五月六月丁香网| 啦啦啦观看免费观看视频高清| 亚洲国产色片| 久久久久久大精品| 村上凉子中文字幕在线| 麻豆成人av在线观看| 欧美潮喷喷水| 麻豆av噜噜一区二区三区| 两个人视频免费观看高清| 亚洲美女黄片视频| ponron亚洲| 国产精品一区二区性色av| 日本黄大片高清| 成人毛片a级毛片在线播放| 欧美区成人在线视频| 亚洲 国产 在线| 国产乱人视频| 91字幕亚洲| 一区二区三区激情视频| 亚洲欧美日韩高清在线视频| 如何舔出高潮| 热99re8久久精品国产| 日本撒尿小便嘘嘘汇集6| 久久久国产成人免费| 亚洲五月天丁香| 十八禁人妻一区二区| 国产精品女同一区二区软件 | 午夜福利欧美成人| 此物有八面人人有两片| 欧美+亚洲+日韩+国产| 黄色日韩在线| 看十八女毛片水多多多| 女同久久另类99精品国产91| 欧美激情久久久久久爽电影| 熟妇人妻久久中文字幕3abv| 看片在线看免费视频| 乱码一卡2卡4卡精品| 99热这里只有精品一区| 日本撒尿小便嘘嘘汇集6| 十八禁人妻一区二区| av在线天堂中文字幕| 国产精品久久久久久精品电影| 午夜福利高清视频| 看片在线看免费视频| 精品久久久久久久久久久久久| 亚洲人与动物交配视频| 性插视频无遮挡在线免费观看| 天堂√8在线中文| 国产精品久久久久久精品电影| 精华霜和精华液先用哪个| 中文亚洲av片在线观看爽| 伊人久久精品亚洲午夜| 天堂网av新在线| 黄色丝袜av网址大全| 99精品在免费线老司机午夜| 在线播放国产精品三级| 欧美3d第一页| 夜夜爽天天搞| 极品教师在线免费播放| 久久人人爽人人爽人人片va | 亚洲经典国产精华液单 | 午夜激情福利司机影院| 亚洲欧美日韩卡通动漫| 黄色女人牲交| 久久久久性生活片| 狠狠狠狠99中文字幕| 亚洲无线在线观看| 久久精品影院6| 亚洲精品亚洲一区二区| 日本黄大片高清| 精品国产亚洲在线| 国产高清视频在线播放一区| 全区人妻精品视频| 看十八女毛片水多多多| 亚洲三级黄色毛片| 蜜桃久久精品国产亚洲av| 黄色丝袜av网址大全| 午夜福利成人在线免费观看| 老司机午夜福利在线观看视频| 极品教师在线视频| 啦啦啦韩国在线观看视频| 亚洲一区二区三区不卡视频| 免费在线观看日本一区| 91狼人影院| 欧美激情久久久久久爽电影| 看十八女毛片水多多多| 在线a可以看的网站| 国语自产精品视频在线第100页| 中文资源天堂在线| 午夜免费男女啪啪视频观看 | 成人av在线播放网站| 美女xxoo啪啪120秒动态图 | 如何舔出高潮| 欧美3d第一页| 欧美黑人巨大hd| 国产精品自产拍在线观看55亚洲| 精品久久久久久久久亚洲 | 亚洲自偷自拍三级| 精品不卡国产一区二区三区| 伦理电影大哥的女人| 国产精品久久久久久久电影| 成人欧美大片| 国产在线男女| 欧美精品国产亚洲| 亚洲va日本ⅴa欧美va伊人久久| 亚洲七黄色美女视频| 听说在线观看完整版免费高清| 日本五十路高清| 亚洲午夜理论影院| 亚洲在线观看片| 成年人黄色毛片网站| 嫁个100分男人电影在线观看| 色综合站精品国产| 国产精品久久久久久精品电影| 国产一区二区三区视频了| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品一卡2卡三卡4卡5卡| 亚洲成人免费电影在线观看| 国产精品久久久久久久电影| 九九在线视频观看精品| 精品国产亚洲在线| 国产精品久久久久久亚洲av鲁大| 男人舔奶头视频| 亚洲国产精品sss在线观看| 亚洲精品日韩av片在线观看| 久久精品国产清高在天天线| 免费电影在线观看免费观看| 日韩欧美免费精品| 嫩草影院精品99| or卡值多少钱| 18禁在线播放成人免费| 极品教师在线免费播放| 午夜福利免费观看在线| 美女xxoo啪啪120秒动态图 | 成年人黄色毛片网站| 午夜精品在线福利| 精品欧美国产一区二区三| 女人被狂操c到高潮| 性欧美人与动物交配| 色视频www国产| 欧美另类亚洲清纯唯美| 国产久久久一区二区三区| 老司机午夜福利在线观看视频| 日本三级黄在线观看| 精品一区二区三区av网在线观看| 免费观看精品视频网站| 精品国内亚洲2022精品成人| 午夜精品一区二区三区免费看| 亚洲无线观看免费| 国产精品98久久久久久宅男小说| 亚洲精品一卡2卡三卡4卡5卡| 亚洲电影在线观看av| 国产亚洲精品综合一区在线观看| 午夜日韩欧美国产| 欧美日本亚洲视频在线播放| av在线蜜桃| 成年女人看的毛片在线观看| 日本与韩国留学比较| 免费av毛片视频| 午夜免费成人在线视频| 午夜福利18| 99久久精品热视频| 男女床上黄色一级片免费看| 亚洲欧美日韩高清专用| 一级作爱视频免费观看| 日韩高清综合在线| 国产人妻一区二区三区在| 久久久国产成人免费| 色播亚洲综合网| 亚洲欧美日韩高清专用| 蜜桃亚洲精品一区二区三区| 久久伊人香网站| 啦啦啦观看免费观看视频高清| 亚洲午夜理论影院| 亚洲经典国产精华液单 | 国产一区二区在线观看日韩| 色吧在线观看| 国产精品日韩av在线免费观看| 亚洲美女视频黄频| 久久国产精品人妻蜜桃| 直男gayav资源| 久久久久精品国产欧美久久久| 丝袜美腿在线中文| 日本一二三区视频观看| 午夜免费激情av| 色精品久久人妻99蜜桃| 一级毛片久久久久久久久女| 我要搜黄色片| 欧美一区二区国产精品久久精品| 亚洲18禁久久av| 国产av不卡久久| 最近视频中文字幕2019在线8| 69av精品久久久久久| 中亚洲国语对白在线视频| 给我免费播放毛片高清在线观看| 99热6这里只有精品| 一区二区三区免费毛片| 免费观看人在逋| 性色av乱码一区二区三区2| 91在线精品国自产拍蜜月| 成人午夜高清在线视频| 日日干狠狠操夜夜爽| 亚洲黑人精品在线| 夜夜爽天天搞| 欧美zozozo另类| 制服丝袜大香蕉在线| 乱人视频在线观看| 性色av乱码一区二区三区2| 91麻豆精品激情在线观看国产| 国产一区二区三区在线臀色熟女| 国产三级黄色录像| 我要搜黄色片| 一边摸一边抽搐一进一小说| 青草久久国产| а√天堂www在线а√下载| 嫩草影院入口| 99在线视频只有这里精品首页| 热99re8久久精品国产| 亚洲欧美精品综合久久99| 脱女人内裤的视频| 久久久久免费精品人妻一区二区| 久久国产精品影院| 91午夜精品亚洲一区二区三区 | 日本a在线网址| 免费搜索国产男女视频| 国产亚洲精品久久久com| 成年人黄色毛片网站| 熟女人妻精品中文字幕| 首页视频小说图片口味搜索| 国产v大片淫在线免费观看| 欧美极品一区二区三区四区| 99热这里只有是精品50| 日韩国内少妇激情av| 久久精品国产亚洲av涩爱 | 少妇人妻精品综合一区二区 | 动漫黄色视频在线观看| 久久热精品热| 日韩人妻高清精品专区| 婷婷色综合大香蕉| 亚洲真实伦在线观看| 午夜视频国产福利| 我的女老师完整版在线观看| 国产精品日韩av在线免费观看| 亚洲五月婷婷丁香| 国产精品亚洲av一区麻豆| 99热这里只有精品一区| 成人一区二区视频在线观看| 99国产精品一区二区三区| 成人国产综合亚洲| 老熟妇仑乱视频hdxx| 最新中文字幕久久久久| 精品国产三级普通话版| 色吧在线观看| 综合色av麻豆| 俄罗斯特黄特色一大片| 51午夜福利影视在线观看| 一夜夜www| 久久久久国内视频| 精品人妻熟女av久视频| 色av中文字幕| 亚洲不卡免费看| 欧美一级a爱片免费观看看| 99热精品在线国产| 18+在线观看网站| 夜夜看夜夜爽夜夜摸| 中文字幕av成人在线电影| av国产免费在线观看| 亚洲av电影在线进入| а√天堂www在线а√下载| 最后的刺客免费高清国语| 无人区码免费观看不卡| 国产单亲对白刺激| 亚洲av成人精品一区久久| 在线国产一区二区在线| 午夜免费成人在线视频| 深夜精品福利| 国产精品久久视频播放| 免费在线观看成人毛片| 老女人水多毛片| 国产伦在线观看视频一区| 欧美三级亚洲精品| 精品久久久久久,| 久久中文看片网| 搡老妇女老女人老熟妇| 久久久久国内视频| x7x7x7水蜜桃| 亚洲美女搞黄在线观看 | 国产熟女xx| 久久久久久九九精品二区国产| 亚洲国产精品999在线| 欧美黑人欧美精品刺激| 女生性感内裤真人,穿戴方法视频| 国产精品爽爽va在线观看网站| 少妇裸体淫交视频免费看高清| 亚洲av一区综合| 日本一二三区视频观看| 午夜日韩欧美国产| 国产精品久久久久久亚洲av鲁大| 国产黄片美女视频| 午夜福利18| 精品久久国产蜜桃| 久久久久九九精品影院| 色视频www国产| 中文字幕高清在线视频| 特级一级黄色大片| 中文字幕av成人在线电影| 免费高清视频大片| 国产亚洲精品久久久com| 亚洲人成网站在线播| 国产aⅴ精品一区二区三区波| 亚洲成av人片在线播放无| 午夜福利高清视频| 国产91精品成人一区二区三区| 日韩av在线大香蕉| 久久99热6这里只有精品| 淫妇啪啪啪对白视频| 亚洲第一区二区三区不卡| 国产精品av视频在线免费观看| 天美传媒精品一区二区| 中文字幕精品亚洲无线码一区| 亚洲人成网站高清观看| 黄色一级大片看看| 18禁黄网站禁片免费观看直播| 最近中文字幕高清免费大全6 | 精品一区二区三区av网在线观看| 麻豆成人av在线观看| 免费电影在线观看免费观看| 色播亚洲综合网| 国产av一区在线观看免费| 国产精品1区2区在线观看.| 亚洲 国产 在线| 亚洲无线在线观看| 男人和女人高潮做爰伦理| 国产亚洲欧美98| 国产精品一区二区免费欧美| 中文字幕av在线有码专区| 一区二区三区激情视频| 一进一出好大好爽视频| 精品免费久久久久久久清纯| 国产黄a三级三级三级人| 深爱激情五月婷婷| 九九久久精品国产亚洲av麻豆| 亚洲性夜色夜夜综合| 美女cb高潮喷水在线观看| 久久午夜亚洲精品久久| 高清毛片免费观看视频网站| 好男人在线观看高清免费视频| 91久久精品国产一区二区成人| 欧美xxxx性猛交bbbb| av黄色大香蕉| 亚洲,欧美,日韩| 日韩免费av在线播放| 日本黄色片子视频| 别揉我奶头 嗯啊视频| 亚洲无线观看免费| 亚洲真实伦在线观看| 欧美性猛交╳xxx乱大交人| 丰满人妻熟妇乱又伦精品不卡| 久久香蕉精品热| 久久久久久久久久黄片| 俺也久久电影网| 午夜久久久久精精品| 免费无遮挡裸体视频| 亚洲国产精品999在线| 日本免费a在线| 少妇人妻精品综合一区二区 | 久久久久久久久久成人| 国产av在哪里看| 亚洲综合色惰| 欧美成狂野欧美在线观看| 亚洲天堂国产精品一区在线| 亚洲经典国产精华液单 | 免费在线观看亚洲国产| 亚洲av五月六月丁香网| 丰满乱子伦码专区| 色综合欧美亚洲国产小说| 久久人人精品亚洲av| 观看美女的网站| 老女人水多毛片| 日韩免费av在线播放| 亚洲精品粉嫩美女一区| 熟女人妻精品中文字幕| av专区在线播放| 久久精品国产亚洲av涩爱 | 人妻制服诱惑在线中文字幕| 亚洲综合色惰| 一边摸一边抽搐一进一小说| 天堂网av新在线| 中文字幕av成人在线电影| 舔av片在线| 麻豆av噜噜一区二区三区| 成人国产综合亚洲| 婷婷六月久久综合丁香| 夜夜爽天天搞| 午夜久久久久精精品| 永久网站在线| 免费高清视频大片| 少妇人妻精品综合一区二区 | 国产av在哪里看| 欧美zozozo另类| 中文字幕久久专区| 最新中文字幕久久久久| bbb黄色大片| 成年版毛片免费区| 国产一区二区三区在线臀色熟女| 91在线精品国自产拍蜜月| 国产久久久一区二区三区| 白带黄色成豆腐渣| 好男人电影高清在线观看| 国产精品一区二区三区四区久久| 在线观看午夜福利视频| 精品久久久久久成人av| 午夜福利在线在线| 性色av乱码一区二区三区2| av黄色大香蕉| 好男人电影高清在线观看| 亚洲精华国产精华精| 淫妇啪啪啪对白视频| 丁香欧美五月| 中国美女看黄片| 国产伦一二天堂av在线观看| 欧美最新免费一区二区三区 | 国产精品久久久久久久电影| 欧美午夜高清在线| 偷拍熟女少妇极品色| 亚洲久久久久久中文字幕| 国产av一区在线观看免费| 国产视频一区二区在线看| 少妇的逼好多水| 国产单亲对白刺激| 亚洲专区国产一区二区| 国产精品免费一区二区三区在线| 色播亚洲综合网| 免费在线观看日本一区| 亚洲精品一卡2卡三卡4卡5卡| 美女黄网站色视频| 精品无人区乱码1区二区| 夜夜躁狠狠躁天天躁| 午夜免费成人在线视频| 亚洲七黄色美女视频| 国产精品一及| 亚洲七黄色美女视频| x7x7x7水蜜桃| 免费在线观看日本一区| 精品久久久久久久久亚洲 | 国产精品不卡视频一区二区 | 床上黄色一级片| 一级av片app| 亚洲 欧美 日韩 在线 免费| 桃红色精品国产亚洲av| 男插女下体视频免费在线播放| 最近最新中文字幕大全电影3| 国产极品精品免费视频能看的| 欧美性猛交黑人性爽| 亚洲久久久久久中文字幕| 亚洲 欧美 日韩 在线 免费| 免费观看人在逋| 日韩欧美国产在线观看| 久久精品人妻少妇| 99国产综合亚洲精品| 中出人妻视频一区二区| 亚洲精品亚洲一区二区| 亚洲欧美清纯卡通| 精品福利观看| 日本免费a在线| 中文字幕高清在线视频| 天堂√8在线中文| 九色国产91popny在线| 男人舔女人下体高潮全视频| 桃色一区二区三区在线观看| 神马国产精品三级电影在线观看| 好男人在线观看高清免费视频| 欧美最新免费一区二区三区 | 禁无遮挡网站| 在线a可以看的网站| 69av精品久久久久久| 好看av亚洲va欧美ⅴa在| 国产一区二区激情短视频| 亚洲在线自拍视频| 久久久久亚洲av毛片大全| 亚洲五月天丁香| 国产伦精品一区二区三区四那| 丁香欧美五月| 国产欧美日韩一区二区精品| 日本熟妇午夜| 国产精品三级大全| 校园春色视频在线观看| 精品不卡国产一区二区三区| 黄色视频,在线免费观看| 国产aⅴ精品一区二区三区波| 人人妻,人人澡人人爽秒播| 狠狠狠狠99中文字幕| 国产在线男女| 九九久久精品国产亚洲av麻豆| 免费在线观看亚洲国产| 国产熟女xx| 麻豆国产97在线/欧美| 久久精品国产亚洲av天美| 神马国产精品三级电影在线观看| 少妇裸体淫交视频免费看高清| 天堂影院成人在线观看| 丁香欧美五月| 欧美中文日本在线观看视频| 此物有八面人人有两片| 亚洲欧美清纯卡通| 亚洲人成网站在线播放欧美日韩| а√天堂www在线а√下载| 赤兔流量卡办理| 国产蜜桃级精品一区二区三区| 可以在线观看的亚洲视频| 国产精品亚洲一级av第二区| 日日干狠狠操夜夜爽| 亚洲不卡免费看| 午夜福利在线观看免费完整高清在 | 淫秽高清视频在线观看| 成人一区二区视频在线观看| 久久久久久久久久黄片| 欧美高清性xxxxhd video| 又粗又爽又猛毛片免费看| 欧美成人一区二区免费高清观看| 一本综合久久免费| 久久久久久大精品| 国产成人啪精品午夜网站| 五月伊人婷婷丁香| 国产成人欧美在线观看| 欧美日韩黄片免| 90打野战视频偷拍视频| 国产在视频线在精品| avwww免费| 国产单亲对白刺激| 亚洲自拍偷在线| 久久99热6这里只有精品| 好男人电影高清在线观看| 国产亚洲精品久久久com| 一级黄色大片毛片| 亚洲av电影在线进入| 最新在线观看一区二区三区| 成人特级黄色片久久久久久久| 噜噜噜噜噜久久久久久91| 亚洲av不卡在线观看| 日本免费一区二区三区高清不卡| 亚洲av免费高清在线观看| 小说图片视频综合网站| 色哟哟·www| 国产大屁股一区二区在线视频| 男女做爰动态图高潮gif福利片| 蜜桃亚洲精品一区二区三区| 日韩大尺度精品在线看网址| www.999成人在线观看| 波多野结衣巨乳人妻| 国内毛片毛片毛片毛片毛片| 日韩欧美国产在线观看| 国产成人福利小说| 特级一级黄色大片| 欧美xxxx黑人xx丫x性爽| 两个人的视频大全免费| 国产亚洲精品久久久com| 国产乱人视频| 熟女人妻精品中文字幕| ponron亚洲| 久久午夜亚洲精品久久| 精品福利观看| 国产一区二区激情短视频| 99热精品在线国产| 天堂√8在线中文|