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

    低信噪比地震數(shù)據圖像結構引導去噪方法與應用

    2024-02-03 13:18:46鄭啟明都小芳吳高奎
    石油物探 2024年1期
    關鍵詞:信噪比剖面濾波器

    鄭啟明,李 琦,都小芳,吳高奎

    (1.北京派特杰奧科技有限公司,北京100044;2.中國地質大學(北京)海洋學院,北京100083;3.中石油東方地球物理公司研究院,河北涿州072751;4.中石化石油勘探開發(fā)研究院,北京102249)

    基于地震資料刻畫斷層、裂縫等結構特征,對于發(fā)現(xiàn)和描述油藏具有重要意義,可支撐決策最佳開發(fā)方案、減少探井風險,提高石油勘探成功率。然而,在地震勘探區(qū)域,由于地質結構比較復雜,斷裂帶的地層形態(tài)多變,因而常導致采集的地震數(shù)據信噪比較低,使得地震數(shù)據處理、解釋十分困難。高品質的地震資料是高質量解釋、油氣儲層精準預測的重要依據和保證。因此,在地震數(shù)據處理時需要采用有效的去噪技術來提高地震資料的品質和可靠性。地震資料處理和解釋人員探索了一系列的地震數(shù)據濾波技術以提高地震資料信噪比,如頻率-波數(shù)域濾波、時域中值濾波等[1],這些技術可以有效去除地震數(shù)據中的部分噪聲,從而提高地震資料解釋和油氣儲層預測的準確性。相比通用的圖像處理算法,地震數(shù)據處理更關注利用反射同相軸的橫向一致性,處理過程中可能會壓制傾斜斷層陰影部分的幅值、破壞小尺度斷層和裂縫等非連續(xù)性的結構信息,導致錯誤匹配斷層兩側的同相軸。因此,進一步探索研究可以保護斷裂、裂縫、尖滅等地質信息的疊后地震保邊去噪技術對地震解釋工作尤為重要。

    為更好地利用地震反射的局部信息,LUO等[2-3]首次將保邊光滑濾波的思想引入地震資料解釋性處理領域。采用Kuwahara多窗分析技術得出,此類方法可以考慮局部結構信息設計自適應光滑濾波器[4-5],在提高信噪比的同時,保護地震圖像中的有效地質信息。此后,國內外學者引入基于各向異性擴散濾波的構造導向濾波技術[6-10],通過結構張量評估局部結構走向、計算各向異性算子,設計增強地震資料沿同相軸方向的信噪比、抑制斷點位置光滑強度的自適應濾波器。在此基礎上,國內外學者優(yōu)化了相關算法,提高了結構張量評估局部層位走向的有效性[11-16],通過偏微分方程的近似解法提高了計算效率[17],并將相關技術擴展應用到三維地震數(shù)據的處理[18-20]。考慮到擴散濾波的有效性依賴于結構張量的準確估計,HALE[21]提出了引入高度自適應的雙邊濾波技術來實現(xiàn)保邊光滑處理的方法,但該方法的計算代價相對較高。國內學者也探索了在變換域[22-24],如曲波域、變分模態(tài)域或多道奇異譜上的數(shù)據濾波去噪處理技術。TRINH等[25]提出了利用Bessel函數(shù)逼近各向異性Laplace濾波的結構導向光滑算法;在圖像處理領域,HE等[26]提出了利用引導圖像濾波進行自適應保邊光滑處理的方法,計算效率顯著高于雙邊濾波方法,該方法尚未在地震數(shù)據處理中廣泛應用。此外,結構導向光滑濾波技術在全波形反演中也具有較明顯的應用效果[27-28]。近年來人工智能技術飛速發(fā)展,高性能濾波算法[29]、神經網絡地震數(shù)據去噪方法[30-31]和智能化斷層表征技術[32-33]也得到了深入的探索。

    在前人研究的基礎上,將“二維圖像濾波”首次擴展應用到二維或三維地震數(shù)據處理中,該計算過程有別于傳統(tǒng)的二維圖像濾波。為抑制隨機與相干噪聲、增強地震反射同相軸的一致性,建立了雙邊顯式加權濾波器的框架,引入了分時窗的引導濾波核函數(shù)與正規(guī)化算子,建立輸入與輸出的目標方程,利用線性回歸顯式求解方程系數(shù),迭代計算直至收斂。該算法采用高效率的保邊光滑處理,使得處理后的地震資料能更加準確地刻畫地質結構,從而更好地服務于地質勘探開發(fā)。

    1 圖像結構引導濾波原理

    針對地震數(shù)據去噪開發(fā)的濾波器可分為非平均濾波器和加權平均濾波器,前者以中值濾波為代表,其濾波過程一般為非線性過程,后者可進一步分為顯式濾波器和隱式濾波器。顯式濾波器通過設計濾波核函數(shù),然后根據核函數(shù)對待處理的數(shù)據進行加權平均實現(xiàn)特定濾波的目的,其中較典型的有均值濾波、高斯濾波。在地震數(shù)據處理中,以隱式加權平均的擴散濾波方法較為常見,其將包含在矩陣中的濾波器作用于待處理數(shù)據,可以較為直接地在矩陣中引入各向異性參數(shù)并求解隱式矩陣來實現(xiàn)結構導向濾波。需要注意的是,各向同性擴散濾波與高斯濾波是等價濾波器。

    定義問題Iout=Iin+n,其中,Iin,Iout,n分別為含噪聲數(shù)據、去噪后的數(shù)據和噪聲。參考顯式加權平均濾波器的設計思路,可利用線性加權平均定義濾波核函數(shù):

    (1)

    (2)

    (3)

    其中,k為窗wk的中心坐標,(ak,bk)是表征該線性關系的參數(shù),假定其在窗wk中為常量。若輸入數(shù)據Iin和輸出數(shù)據Iout的差異為隨機或相干噪聲,則確定參數(shù)(ak,bk)的最小化目標函數(shù)為:

    (4)

    其中,ξ是正則化算子系數(shù),用以確保ak的解在可接受的數(shù)值范圍內。利用線性回歸算法,可顯式地求解(ak,bk)的值,即:

    (5)

    (6)

    分時窗處理過程中,由于在時窗交疊區(qū)的輸出數(shù)據Iout具有非唯一性,因此,在同一點處多個交疊窗的平均輸出表達式為:

    (7)

    (8)

    傳統(tǒng)的雙邊濾波器可以利用引導數(shù)據通過顯式加權平均進行保邊結構光滑,其核心在于使用線性加權平均濾波器,可以引入濾波核函數(shù)。引導濾波器采用同樣思想,但通過構建最優(yōu)化問題引入了隱式的加權平均,通過線性最優(yōu)化解獲得了引導濾波器算子,同時在優(yōu)化問題中引入正則化算子形成超參數(shù),可控制對于邊界結構的敏感程度。

    以上推導過程適用于二維或三維地震數(shù)據體。在具體應用過程中,需要準備引導數(shù)據G(原始輸入數(shù)據或濾波處理后的窄帶高信噪比數(shù)據),試驗確定半窗寬度r和正則化因子ξ兩個參數(shù)。當r和ξ變大時,處理后的數(shù)據具有更強的光滑效果,可能會損傷地震反射的斷裂結構信息,實際計算中需要根據保邊處理的目的與需求,綜合考慮確定合適的處理參數(shù)。由公式(5)至公式(7)可知,該算法可多次迭代應用方框濾波器實現(xiàn)地震數(shù)據的去噪處理,自動具備O(n)的計算效率。

    具體應用于實際地震數(shù)據處理時,采用以下技術流程:

    1) 準備待處理的地震數(shù)據Iin(二維或三維地震數(shù)據)、結構引導數(shù)據G(可以與Iin相同)、設置算法參數(shù)半窗寬度r和正則化因子ξ;

    2) 計算f(Iin),f(G),f(Iin·Iin)與f(Iin·G),其中,f為均值濾波器;

    3) 計算參數(shù)a=[f(Iin·G)]/[f(Iin·Iin)-f(Iin)f(Iin)+ξ];

    4) 計算參數(shù)b=f(G)-af(Iin);

    5) 計算濾波后的數(shù)據Iout=f(a)Iin+f(b)。

    此流程主要計算6次均值濾波,由于均值濾波器可以使用方框濾波器的方法構建,因此本流程具有O(n)復雜度的計算效率。

    該流程中結構信息的引入依賴于結構引導數(shù)據G的選擇。引導濾波應用過程中利用均值濾波器提取G的結構信息。假如G由多個均質的區(qū)塊構成,那么引導濾波會在各個區(qū)塊內分別進行操作,而不會在光滑區(qū)塊的邊界進行均值濾波,從而有效保留邊緣信息,實現(xiàn)保邊濾波的效果。若G是輸入數(shù)據本身,則引導濾波器會強化數(shù)據本身的走向或者傾斜結構,同時具有保留邊界、突出主要結構的功能。導向濾波器中的引導數(shù)據G的選取可根據具體應用調整,本次研究選擇輸入數(shù)據作為引導數(shù)據,是較為直接和有效的方法。

    2 模型數(shù)據測試

    為驗證圖像結構引導濾波算法的穩(wěn)定性與結果的可靠性,對比該算法與傳統(tǒng)濾波器及主流濾波器的優(yōu)劣,基于Marmousi模型通過卷積15Hz主頻的Ricker子波和反射系數(shù),在設計的二維合成地震剖面(圖1a)中加入一定噪聲得到信噪比(signal to noise ratio,SNR)約8dB的含噪數(shù)據(圖1b),以此評估后續(xù)各類濾波算法的去噪效果。其中,SNR定義如下:

    (9)

    圖1 基于Marmousi模型的合成地震剖面去噪測試對比a 無噪剖面; b 含噪剖面; c 圖像結構引導濾波后的結果; d 構造導向濾波后的結果; e 中值濾波后的結果; f 均值濾波后的結果; g 高斯濾波后的結果

    其中,Isignal表示不含噪聲的有效信號,Isignal+noise表示含噪聲的信號,SNR表示計算得到的信噪比。

    為方便對比不同方法的去噪效果,試驗中選用了相同的半窗寬度,其中圖像結構引導濾波與構造導向濾波皆以原始數(shù)據本身作為引導數(shù)據提取結構張量進行處理,含噪地震數(shù)據的自適應引導圖像濾波、構造導向濾波、中值濾波、均值濾波和高斯濾波的處理結果分別如圖1c至圖1g所示。

    由模型去噪結果可知,圖像結構引導濾波(圖1c)、構造導向濾波(圖1d)、中值濾波(圖1e)、均值濾波(圖1f)、高斯濾波(圖1g)均有效消除了部分噪聲,提高了模型數(shù)據的信噪比。前兩種方法的信噪比由8.00dB分別提高至18.31,14.77dB,與構造導向濾波相比,圖像結構引導濾波的信噪比高約4dB,對地層邊界具有更好的保護作用;而后3種方法的信噪比由8.00dB分別提高至9.19,9.58,10.09dB,信噪比雖有所提高,但在一定程度上模糊了地層與斷裂的反射邊界。

    為便于清晰展示上述去噪效果,圖2更詳細地對比了不同方法分離的噪聲剖面。圖1與圖2的截斷幅值不同。在層位相對平緩的區(qū)域,如距離為8km,時間為1s處,圖像結構引導濾波和構造導向濾波均較好地保持了地震數(shù)據的結構信息(圖2b、圖2c),而中值濾波、均值濾波和高斯濾波因過度光滑破壞了有效地震反射的真實結構信息(圖2d、圖2e、圖2f)。針對本合成數(shù)據中其它斜率較大的結構,從噪聲剖面可知,構造導向濾波(圖2c)保留邊界的能力相對中值濾波(圖2d)、均值濾波(圖2e)和高斯濾波(圖2f)更強,但是圖像結構引導濾波(圖2b)更好地保留了地震反射細微的結構。總之,相比其它處理方法,圖像結構引導濾波不但能有效抑制噪聲,而且在處理方法與計算效率上也更加簡單高效。

    圖2 不同方法分離的噪聲剖面對比a 原始噪聲; b 圖像結構引導濾波; c 構造導向濾波; d 中值濾波; e 均值濾波; f 高斯濾波

    3 實際應用效果

    以信噪比較低、含有斷層的三維地震數(shù)據為例,利用含噪數(shù)據為引導,基于核濾波函數(shù)自適應地對其進行濾波處理。圖3和圖4分別為自適應濾波前、后不同頻帶三維地震剖面和沿層振幅切片。由圖3和圖4可見,自適應濾波處理前,低頻地震數(shù)據的地震剖面和沿層振幅切片中存在明顯低頻噪聲(圖3a、圖4a);受噪聲影響,中頻信號與小尺度斷層反射信息相混疊,導致屬性邊界模糊不清(圖3c、圖4c);高頻地震數(shù)據的地震剖面和沿層振幅切片中的高頻有效信號受噪聲影響更加嚴重,信噪比極低,同相軸連續(xù)性極差(圖3e、圖4e)。

    圖3 自適應濾波前、后不同頻帶三維地震剖面a 低頻含噪信號; b 低頻圖像結構引導濾波結果; c 中頻含噪信號; d 中頻圖像結構引導濾波結果; e 高頻含噪信號; f 高頻圖像結構引導濾波結果

    圖4 濾波前、后不同頻帶沿層振幅切片屬性a 低頻含噪信號; b 低頻圖像結構引導濾波結果; c 中頻含噪信號; d 中頻圖像結構引導濾波結果; e 高頻含噪信號; f 高頻圖像結構引導濾波結果

    與圖3a、圖4a低頻信息相比,圖像結構引導濾波很好地壓制了低頻數(shù)據中的散射噪聲,同時很好地保持大尺度斷層的錯斷邊界(圖3b、圖4b);其次,對比圖3c、圖3d和圖4c、圖4d可見,本文濾波處理方法可有效保邊光滑去除噪聲、提高信噪比,同時保持小尺度斷層、沉積體邊界的有效反射信息(圖3d、圖4d 標識區(qū)域);對比圖3e、圖3f與圖4e、圖4f可見,本文方法可自適應利用高信噪比優(yōu)勢頻帶信號作為引導圖像數(shù)據,有效壓制高頻噪聲、提高地震同相軸連續(xù)性與信噪比(圖3f與圖4f方框區(qū)域),有助于解釋人員提高小斷層和地震屬性的解釋與描述精度。

    構造導向與圖像結構引導濾波的全頻帶處理結果如圖5所示。在含噪地震數(shù)據(圖5a)上提取結構張量的傾斜構造導向濾波結果如圖5b與圖5c所示;以含噪數(shù)據自身為引導的圖像結構引導濾波處理結果如圖5d 所示。由圖5可見,采用構造導向濾波與自適應圖像結構引導濾波處理后的地震剖面信噪比顯著提高,除大斷層特征更加明顯之外,小斷裂的細節(jié)更加清晰。對比圖5b、圖5c可見,當構造導向濾波的相關半徑較大時,會導致計算效率降低、橫向保邊結果略差;當相關半徑較小時,局部產生相干噪聲異常(圖5c中紅圈),散射噪聲壓制不徹底(圖5b、圖5c黃框)。

    圖5 原始地震數(shù)據、構造導向濾波與自適應圖像結構引導濾波結果a 原始地震數(shù)據; b 小半徑構造導向濾波結果; c 大半徑構造導向濾波結果; d 自適應圖像結構引導濾波結果

    自適應圖像結構引導濾波可較為徹底地壓制相干噪聲、保護有效反射信號(圖5d)。在算法上,其相比于構造導向濾波,可自適應地從原始含噪地震數(shù)據中優(yōu)選引導圖像數(shù)據與處理參數(shù),無需象構造導向濾波那樣為獲取傾角信息而多次調整參數(shù)來計算結構張量并重復進行引導濾波;在計算效率上,其利用的線性函數(shù)可有效降低計算復雜度、提高運算效率,計算代價僅相當于數(shù)次均值濾波。

    對原始地震數(shù)據、構造導向濾波和自適應圖像結構引導濾波結果分別利用相同的特征值參數(shù)提取相干體剖面與沿層相干屬性信息進一步評價處理效果(圖6、圖7)。

    圖6 原始地震數(shù)據、構造導向濾波與自適應引導圖像濾波后相干體剖面a 原始地震數(shù)據相干體剖面; b 小半徑構造導向濾波后相干體剖面; c 大半徑構造導向濾波后相干體剖面; d 自適應圖像結構引導濾波后相干體剖面

    圖7 原始地震數(shù)據(a)與自適應圖像結構引導濾波后(b)沿層相干屬性

    在斷裂位置附近,對比原始數(shù)據與不同構造導向半徑獲取的相干結果(圖6a、圖6b與圖6c),因受散射噪聲影響,導致其對斷層的分辨能力大大降低,給后續(xù)的解釋工作造成困難。此類污染主要由于斷裂帶附近的地震反射資料信噪比低、同相軸的連續(xù)性不夠明顯,導致計算的相干體屬性難以區(qū)分地震層位與斷裂結構特征;與其它方法相比,本文方法可有效壓制地震數(shù)據中的噪聲,因此可獲得信噪比更高的相干體,從而能夠更加清晰地描述復雜的斷裂結構(圖6d);同時,沿層相干屬性也能夠很好地描述細微斷裂結構,多級斷裂在橫向上的層次關系也更加清晰(圖7)。該結果有助于后續(xù)的層位與斷裂解釋。

    4 討論

    圖像結構引導濾波的引導數(shù)據G的設計具有關鍵作用。傳統(tǒng)的構造導向濾波需要利用結構張量等方法,顯式地提取結構體在某個區(qū)域內的傾角和方位角,而本文方法的引導圖像數(shù)據不需要顯式表征結構傾角和方位角,只需要其梯度具備與期望輸出圖像Iout一致的邊界信息。同時由于其具有滑動窗實現(xiàn)的特點,以含噪聲的圖像作為引導圖像,可在有效保留結構邊界的前提下壓制隨機噪聲。然而,在實際地震數(shù)據處理中,相干噪聲、散射噪聲是影響地震圖像質量的關鍵,僅以含噪圖像作為引導可能會增強相關噪聲的強度。但由于引導圖像本身設計的靈活性,可采用以下步驟壓制相干噪聲和散射噪聲:①先采用自引導的方式壓制隨機噪聲;②再利用均值、中值、高斯或構造導向濾波處理剖面,獲得邊界受到一定破壞但相干、散射噪聲得到壓制的剖面;③再以步驟②的輸出結果作為引導圖像處理步驟①的輸出結果。但是,此方法仍具有一定局限性,其有效范圍僅限于相干或散射噪聲相對較弱的狀態(tài),否則其它濾波器在壓制相干或散射噪聲時會過度傷害主要結構信息。需要在未來工作中進一步改進圖像結構引導濾波,形成自適應壓制相干或散射噪聲的方法。

    5 結論

    在進行地震資料的斷層、裂縫解釋時,噪聲的存在會明顯影響解釋精度。通過引入圖像結構引導濾波方法,能自適應地提取有效地震結構引導信息,快速、有效地對地震數(shù)據進行保邊、保幅去噪處理。圖像結構引導濾波基于理論推導可通過方窗濾波器實現(xiàn),對于大尺度地震數(shù)據不存在顯著的計算瓶頸,具有O(n)的計算效率,相比工業(yè)界更常用的構造導向濾波器也有一定優(yōu)勢。含噪合成模型數(shù)據的處理分析結果證實了圖像結構引導濾波可保留反射結構邊界、有效壓制噪聲,結果顯著優(yōu)于常規(guī)非結構性光滑濾波器。針對實際地震數(shù)據的處理,利用引導圖像可靈活設計引導特征,在一定程度上壓制相干噪聲與隨機噪聲,提高地震數(shù)據信噪比,增強地震反射同相軸的一致性、保留并突出斷裂與地層反射的非連續(xù)特征。計算的相干屬性較其它方法斷裂結構特征更加清晰,地質規(guī)律更加明顯,進一步證明了方法的有效性與實用性。

    猜你喜歡
    信噪比剖面濾波器
    三點法定交叉剖面方法
    ——工程地質勘察中,一種做交叉剖面的新方法
    基于深度學習的無人機數(shù)據鏈信噪比估計算法
    從濾波器理解卷積
    電子制作(2019年11期)2019-07-04 00:34:38
    開關電源EMI濾波器的應用方法探討
    電子制作(2018年16期)2018-09-26 03:26:50
    低信噪比下LFMCW信號調頻參數(shù)估計
    電子測試(2018年11期)2018-06-26 05:56:02
    低信噪比下基于Hough變換的前視陣列SAR稀疏三維成像
    雷達學報(2017年3期)2018-01-19 02:01:27
    基于曲線擬合的投棄式剖面儀電感量算法
    電子測試(2017年12期)2017-12-18 06:35:46
    復雜多約束條件通航飛行垂直剖面規(guī)劃方法
    基于Canny振蕩抑制準則的改進匹配濾波器
    基于TMS320C6678的SAR方位向預濾波器的并行實現(xiàn)
    欧美极品一区二区三区四区| 欧美日本亚洲视频在线播放| 老熟妇乱子伦视频在线观看| 亚洲av中文字字幕乱码综合| 人妻丰满熟妇av一区二区三区| 欧美丝袜亚洲另类 | 男女那种视频在线观看| 欧美成人性av电影在线观看| 女人高潮潮喷娇喘18禁视频| 国语自产精品视频在线第100页| 精品久久久久久久人妻蜜臀av| 夜夜躁狠狠躁天天躁| 成在线人永久免费视频| 亚洲精品色激情综合| 日本黄色视频三级网站网址| 12—13女人毛片做爰片一| 嫩草影视91久久| 两个人的视频大全免费| 久久精品人妻少妇| 麻豆久久精品国产亚洲av| 久热爱精品视频在线9| 国产一级毛片七仙女欲春2| 国产亚洲精品一区二区www| 午夜激情av网站| 老司机在亚洲福利影院| 欧美三级亚洲精品| 人妻夜夜爽99麻豆av| 日本黄色视频三级网站网址| 麻豆久久精品国产亚洲av| 男人舔奶头视频| 国产av在哪里看| 欧美日本亚洲视频在线播放| 国产精品香港三级国产av潘金莲| 男人舔女人下体高潮全视频| 日韩欧美一区二区三区在线观看| 亚洲精品美女久久久久99蜜臀| 免费电影在线观看免费观看| 黄色视频,在线免费观看| 免费搜索国产男女视频| a在线观看视频网站| 免费av毛片视频| 在线观看免费视频日本深夜| 成人av一区二区三区在线看| 国产视频一区二区在线看| 身体一侧抽搐| 日本免费一区二区三区高清不卡| x7x7x7水蜜桃| 男人舔奶头视频| 99久久精品热视频| 久久这里只有精品中国| 他把我摸到了高潮在线观看| 日韩欧美三级三区| 精品久久久久久久毛片微露脸| 大型av网站在线播放| 亚洲色图 男人天堂 中文字幕| 成人三级做爰电影| svipshipincom国产片| 久久精品国产综合久久久| 欧美绝顶高潮抽搐喷水| 精品国产超薄肉色丝袜足j| 人人妻人人澡欧美一区二区| 免费av毛片视频| 18禁美女被吸乳视频| 久9热在线精品视频| 成人国语在线视频| 蜜桃久久精品国产亚洲av| 一区二区三区国产精品乱码| 国产伦一二天堂av在线观看| 亚洲av成人不卡在线观看播放网| 欧美黑人欧美精品刺激| 一个人免费在线观看电影 | 男女视频在线观看网站免费 | 18禁观看日本| 两个人免费观看高清视频| 亚洲成人国产一区在线观看| 两个人视频免费观看高清| 午夜激情福利司机影院| 日本一本二区三区精品| 亚洲午夜精品一区,二区,三区| 亚洲中文日韩欧美视频| 精品久久久久久久久久免费视频| 亚洲av五月六月丁香网| 18禁观看日本| 丰满人妻一区二区三区视频av | 18禁裸乳无遮挡免费网站照片| 亚洲成a人片在线一区二区| 一级毛片女人18水好多| 在线a可以看的网站| 婷婷精品国产亚洲av在线| 亚洲,欧美精品.| 一级片免费观看大全| 亚洲av电影不卡..在线观看| 91在线观看av| 久久欧美精品欧美久久欧美| 男女床上黄色一级片免费看| 日本撒尿小便嘘嘘汇集6| 国产精品野战在线观看| 欧美不卡视频在线免费观看 | 国产高清有码在线观看视频 | 国产亚洲精品综合一区在线观看 | 国产熟女午夜一区二区三区| 蜜桃久久精品国产亚洲av| 亚洲av中文字字幕乱码综合| 国产一区二区三区视频了| 久久天堂一区二区三区四区| 听说在线观看完整版免费高清| 91国产中文字幕| 亚洲男人的天堂狠狠| 亚洲在线自拍视频| 成人高潮视频无遮挡免费网站| 欧美日韩黄片免| 国产在线观看jvid| 久久久国产精品麻豆| 亚洲欧美一区二区三区黑人| 两性夫妻黄色片| xxxwww97欧美| 色av中文字幕| 啦啦啦观看免费观看视频高清| 亚洲欧美日韩东京热| 久久国产精品影院| 男人舔女人的私密视频| 欧美黑人精品巨大| 亚洲欧美一区二区三区黑人| 精品熟女少妇八av免费久了| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲人成电影免费在线| 午夜久久久久精精品| 国产午夜精品论理片| 超碰成人久久| 757午夜福利合集在线观看| 亚洲精品在线美女| 国产av在哪里看| 色噜噜av男人的天堂激情| 一区福利在线观看| 国产精品一区二区三区四区久久| 久久99热这里只有精品18| 国产欧美日韩精品亚洲av| 中文字幕熟女人妻在线| 免费观看人在逋| 99久久无色码亚洲精品果冻| 一级作爱视频免费观看| 国产真实乱freesex| 国产69精品久久久久777片 | 国产黄色小视频在线观看| 国产成人一区二区三区免费视频网站| 一本精品99久久精品77| 熟妇人妻久久中文字幕3abv| 午夜福利在线在线| 国产精品九九99| 俺也久久电影网| 法律面前人人平等表现在哪些方面| 国产激情欧美一区二区| 又黄又粗又硬又大视频| 999久久久国产精品视频| 亚洲专区中文字幕在线| 色综合亚洲欧美另类图片| 久久中文字幕人妻熟女| 国产成人av教育| 欧美大码av| 国产乱人伦免费视频| 亚洲成人久久性| 精品乱码久久久久久99久播| 俄罗斯特黄特色一大片| 免费在线观看黄色视频的| 免费观看精品视频网站| 不卡av一区二区三区| 国内少妇人妻偷人精品xxx网站 | ponron亚洲| 国产高清视频在线观看网站| 亚洲国产精品成人综合色| 国产成人影院久久av| 日本三级黄在线观看| 欧美绝顶高潮抽搐喷水| 国产精品久久久av美女十八| 最近视频中文字幕2019在线8| 久久精品成人免费网站| 欧洲精品卡2卡3卡4卡5卡区| 超碰成人久久| 天天躁狠狠躁夜夜躁狠狠躁| 露出奶头的视频| 一级片免费观看大全| 欧美人与性动交α欧美精品济南到| 熟女少妇亚洲综合色aaa.| 精品国产乱子伦一区二区三区| 国产精品免费视频内射| 日韩欧美精品v在线| 色综合亚洲欧美另类图片| 国产成人精品无人区| 一进一出抽搐动态| 日本三级黄在线观看| 99精品在免费线老司机午夜| 麻豆成人av在线观看| av天堂在线播放| 欧美黑人欧美精品刺激| 亚洲一码二码三码区别大吗| 久久精品成人免费网站| 久9热在线精品视频| 国产精品永久免费网站| 高潮久久久久久久久久久不卡| 精品福利观看| 999久久久国产精品视频| 床上黄色一级片| 久久精品成人免费网站| 伦理电影免费视频| 最近最新中文字幕大全免费视频| 午夜福利视频1000在线观看| 欧美黑人巨大hd| 99久久综合精品五月天人人| 亚洲av成人一区二区三| 成人av一区二区三区在线看| 精品乱码久久久久久99久播| 99久久无色码亚洲精品果冻| 老司机福利观看| 国产激情偷乱视频一区二区| 全区人妻精品视频| 啦啦啦免费观看视频1| 精品久久久久久久毛片微露脸| 国产欧美日韩一区二区精品| 亚洲国产精品久久男人天堂| 婷婷亚洲欧美| 美女大奶头视频| 中文字幕熟女人妻在线| tocl精华| 色尼玛亚洲综合影院| 国产午夜精品论理片| 亚洲美女黄片视频| 久久久水蜜桃国产精品网| tocl精华| 国产人伦9x9x在线观看| 亚洲精品一卡2卡三卡4卡5卡| 在线视频色国产色| 欧美日韩乱码在线| 精品无人区乱码1区二区| 免费一级毛片在线播放高清视频| 91在线观看av| av天堂在线播放| 巨乳人妻的诱惑在线观看| 少妇裸体淫交视频免费看高清 | 夜夜夜夜夜久久久久| 日韩欧美在线乱码| 久久人人精品亚洲av| 亚洲国产欧美人成| 热99re8久久精品国产| 国内精品久久久久精免费| 久久婷婷成人综合色麻豆| 19禁男女啪啪无遮挡网站| 亚洲人与动物交配视频| 男男h啪啪无遮挡| 波多野结衣高清无吗| 国产精品爽爽va在线观看网站| 亚洲欧美一区二区三区黑人| 久99久视频精品免费| xxxwww97欧美| 熟妇人妻久久中文字幕3abv| 成人三级黄色视频| 亚洲狠狠婷婷综合久久图片| av视频在线观看入口| 精品久久久久久成人av| 欧美日韩瑟瑟在线播放| 啦啦啦免费观看视频1| 欧美日韩亚洲综合一区二区三区_| 欧美色视频一区免费| xxx96com| 欧美日韩黄片免| 免费电影在线观看免费观看| 男女视频在线观看网站免费 | 日本熟妇午夜| 午夜影院日韩av| 一个人免费在线观看电影 | 免费在线观看成人毛片| 香蕉久久夜色| 精华霜和精华液先用哪个| a级毛片在线看网站| 欧美国产日韩亚洲一区| 国产乱人伦免费视频| 亚洲,欧美精品.| 香蕉国产在线看| 一本精品99久久精品77| 在线观看免费视频日本深夜| 亚洲熟妇熟女久久| 美女 人体艺术 gogo| 熟妇人妻久久中文字幕3abv| 久久精品夜夜夜夜夜久久蜜豆 | 丰满人妻熟妇乱又伦精品不卡| 黑人巨大精品欧美一区二区mp4| 丰满人妻一区二区三区视频av | 国语自产精品视频在线第100页| www日本黄色视频网| 天堂动漫精品| 久久精品国产亚洲av高清一级| 亚洲欧美日韩高清专用| 脱女人内裤的视频| 日韩成人在线观看一区二区三区| 一个人免费在线观看电影 | 麻豆av在线久日| 久久人妻av系列| 亚洲一区二区三区色噜噜| 床上黄色一级片| 99riav亚洲国产免费| xxxwww97欧美| 国产精品一区二区免费欧美| 中文资源天堂在线| 一边摸一边做爽爽视频免费| 亚洲乱码一区二区免费版| a级毛片在线看网站| 啦啦啦观看免费观看视频高清| 亚洲人成网站在线播放欧美日韩| 久久中文看片网| 少妇裸体淫交视频免费看高清 | 一级片免费观看大全| 亚洲av熟女| 亚洲七黄色美女视频| 极品教师在线免费播放| 午夜福利免费观看在线| 国产欧美日韩精品亚洲av| 精品乱码久久久久久99久播| 国产免费av片在线观看野外av| 丁香欧美五月| 久久人人精品亚洲av| 精品免费久久久久久久清纯| 久久中文字幕人妻熟女| 露出奶头的视频| 亚洲无线在线观看| 亚洲精品粉嫩美女一区| 精品乱码久久久久久99久播| 丰满的人妻完整版| av超薄肉色丝袜交足视频| 免费观看精品视频网站| 午夜福利18| 国产精品精品国产色婷婷| 在线免费观看的www视频| 国产av又大| 欧美不卡视频在线免费观看 | 欧美3d第一页| 国产精品久久久久久久电影 | 国产精品久久久久久亚洲av鲁大| 91字幕亚洲| 日本黄色视频三级网站网址| 日本五十路高清| 成人午夜高清在线视频| 在线播放国产精品三级| 国产精品久久久久久久电影 | 国产久久久一区二区三区| av福利片在线观看| 美女黄网站色视频| 亚洲黑人精品在线| 丰满人妻熟妇乱又伦精品不卡| 老司机靠b影院| 久久中文字幕一级| 国内精品久久久久久久电影| 精品一区二区三区av网在线观看| 少妇熟女aⅴ在线视频| 黄片小视频在线播放| 午夜激情福利司机影院| 日韩精品中文字幕看吧| 国产精品一及| 国产亚洲欧美98| 亚洲中文字幕日韩| 舔av片在线| 91大片在线观看| 精品久久久久久,| 在线观看午夜福利视频| netflix在线观看网站| 宅男免费午夜| 在线看三级毛片| 亚洲成人免费电影在线观看| 俄罗斯特黄特色一大片| 91大片在线观看| 国产一区二区在线观看日韩 | 国产精品一区二区免费欧美| www.自偷自拍.com| cao死你这个sao货| 欧美最黄视频在线播放免费| 国产三级在线视频| 午夜精品在线福利| 搞女人的毛片| 深夜精品福利| 热99re8久久精品国产| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品野战在线观看| 亚洲自拍偷在线| 免费看十八禁软件| 一本综合久久免费| 国产一区二区在线观看日韩 | 国产精品99久久99久久久不卡| 欧美性猛交黑人性爽| 99热这里只有精品一区 | 亚洲一区中文字幕在线| 亚洲av美国av| 色综合欧美亚洲国产小说| 999精品在线视频| 成人午夜高清在线视频| 成人精品一区二区免费| 欧美黑人精品巨大| 亚洲一卡2卡3卡4卡5卡精品中文| 久久精品国产综合久久久| 免费av毛片视频| 午夜激情福利司机影院| 免费无遮挡裸体视频| av有码第一页| 日本在线视频免费播放| 日本免费一区二区三区高清不卡| 免费人成视频x8x8入口观看| 亚洲国产中文字幕在线视频| 国产爱豆传媒在线观看 | 国产成人av激情在线播放| 男女午夜视频在线观看| 这个男人来自地球电影免费观看| 亚洲av中文字字幕乱码综合| 两个人免费观看高清视频| 国产97色在线日韩免费| 亚洲成人久久性| 国产成人精品无人区| 日日摸夜夜添夜夜添小说| 国产在线观看jvid| 国产欧美日韩一区二区精品| 婷婷亚洲欧美| 国产精品爽爽va在线观看网站| svipshipincom国产片| 国产精品影院久久| 岛国在线免费视频观看| 看免费av毛片| 国产欧美日韩一区二区三| 桃色一区二区三区在线观看| 三级国产精品欧美在线观看 | 搡老妇女老女人老熟妇| 亚洲av成人精品一区久久| 国产男靠女视频免费网站| av视频在线观看入口| 看免费av毛片| 国内揄拍国产精品人妻在线| 精品熟女少妇八av免费久了| 老熟妇乱子伦视频在线观看| 久久精品综合一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 97碰自拍视频| 一本久久中文字幕| 高潮久久久久久久久久久不卡| 亚洲精品一卡2卡三卡4卡5卡| 天天添夜夜摸| 九色国产91popny在线| 国产男靠女视频免费网站| 毛片女人毛片| 婷婷亚洲欧美| 淫秽高清视频在线观看| 国产又黄又爽又无遮挡在线| 国产高清激情床上av| 色av中文字幕| 亚洲美女黄片视频| 99国产精品一区二区蜜桃av| 岛国在线免费视频观看| 黄色片一级片一级黄色片| 亚洲电影在线观看av| 国产私拍福利视频在线观看| 亚洲第一欧美日韩一区二区三区| 精品一区二区三区视频在线观看免费| 免费人成视频x8x8入口观看| 制服人妻中文乱码| 亚洲专区中文字幕在线| 俺也久久电影网| 成年版毛片免费区| 天天添夜夜摸| 人人妻人人看人人澡| 欧美日韩中文字幕国产精品一区二区三区| 制服人妻中文乱码| 久久婷婷人人爽人人干人人爱| www国产在线视频色| 婷婷精品国产亚洲av在线| 午夜精品久久久久久毛片777| 亚洲精品在线美女| 国产伦人伦偷精品视频| 亚洲一区二区三区不卡视频| 久久久精品欧美日韩精品| 国产成人av激情在线播放| 午夜福利欧美成人| e午夜精品久久久久久久| 国产高清视频在线播放一区| 午夜精品久久久久久毛片777| 亚洲精品中文字幕在线视频| 丝袜美腿诱惑在线| 人妻丰满熟妇av一区二区三区| 色在线成人网| 黄片大片在线免费观看| 亚洲欧美精品综合一区二区三区| 成人午夜高清在线视频| 国内精品一区二区在线观看| 亚洲精品中文字幕一二三四区| 国产伦一二天堂av在线观看| 精品久久久久久久毛片微露脸| 亚洲av中文字字幕乱码综合| 床上黄色一级片| 成人一区二区视频在线观看| 亚洲精品美女久久久久99蜜臀| 国产在线精品亚洲第一网站| 婷婷精品国产亚洲av在线| 97超级碰碰碰精品色视频在线观看| 国产一级毛片七仙女欲春2| 一个人观看的视频www高清免费观看 | 国产精品亚洲一级av第二区| 日韩中文字幕欧美一区二区| 国产三级中文精品| 嫩草影院精品99| 欧美乱码精品一区二区三区| 免费在线观看完整版高清| 日韩欧美国产在线观看| 色综合亚洲欧美另类图片| 国产三级黄色录像| 看免费av毛片| 国产精品,欧美在线| 国产aⅴ精品一区二区三区波| 中文字幕人成人乱码亚洲影| 成年免费大片在线观看| 黄色成人免费大全| а√天堂www在线а√下载| 成人国语在线视频| 日韩欧美精品v在线| 欧美中文日本在线观看视频| 成人av一区二区三区在线看| 欧美成人一区二区免费高清观看 | 亚洲人成77777在线视频| 男插女下体视频免费在线播放| 无限看片的www在线观看| 午夜免费观看网址| av片东京热男人的天堂| 久久欧美精品欧美久久欧美| 午夜精品久久久久久毛片777| 亚洲中文字幕一区二区三区有码在线看 | 91字幕亚洲| 中文字幕精品亚洲无线码一区| 成人av一区二区三区在线看| 亚洲欧美日韩高清专用| 亚洲专区国产一区二区| 99精品欧美一区二区三区四区| 亚洲人成网站高清观看| 国产高清激情床上av| 久久香蕉激情| 午夜视频精品福利| 男插女下体视频免费在线播放| 成人特级黄色片久久久久久久| 18禁裸乳无遮挡免费网站照片| 国产精品1区2区在线观看.| 免费人成视频x8x8入口观看| 性欧美人与动物交配| 免费观看精品视频网站| 99在线人妻在线中文字幕| 亚洲精品国产一区二区精华液| 又粗又爽又猛毛片免费看| 亚洲精品久久国产高清桃花| 中文资源天堂在线| 亚洲av成人av| 女生性感内裤真人,穿戴方法视频| 丝袜美腿诱惑在线| 在线观看舔阴道视频| 亚洲精品美女久久久久99蜜臀| 亚洲免费av在线视频| 亚洲成人中文字幕在线播放| 欧美黄色淫秽网站| 亚洲va日本ⅴa欧美va伊人久久| 亚洲人成77777在线视频| 国产免费男女视频| 久久久久久久精品吃奶| av福利片在线观看| 一本久久中文字幕| 亚洲精品中文字幕在线视频| 男插女下体视频免费在线播放| 精品午夜福利视频在线观看一区| 国产一区在线观看成人免费| 法律面前人人平等表现在哪些方面| av在线播放免费不卡| 亚洲欧洲精品一区二区精品久久久| 黄片小视频在线播放| 日日摸夜夜添夜夜添小说| 国产精品自产拍在线观看55亚洲| netflix在线观看网站| 久久中文字幕一级| 精品第一国产精品| 国产一区二区激情短视频| 香蕉久久夜色| 99国产精品一区二区蜜桃av| 亚洲人成网站在线播放欧美日韩| 18禁裸乳无遮挡免费网站照片| 亚洲精品中文字幕在线视频| 热99re8久久精品国产| 久久久国产欧美日韩av| 两性夫妻黄色片| 三级男女做爰猛烈吃奶摸视频| 青草久久国产| 亚洲精品一卡2卡三卡4卡5卡| 禁无遮挡网站| 麻豆久久精品国产亚洲av| 日本一区二区免费在线视频| 成人三级黄色视频| 级片在线观看| 国产高清视频在线观看网站| 亚洲国产日韩欧美精品在线观看 | 国产精品乱码一区二三区的特点| 啪啪无遮挡十八禁网站| 一本大道久久a久久精品| 夜夜爽天天搞| 亚洲成人免费电影在线观看| 亚洲美女黄片视频| 国产三级黄色录像| 久久久久久久久中文| 亚洲成人国产一区在线观看| 精品欧美国产一区二区三| 女人被狂操c到高潮| 久久亚洲精品不卡| xxx96com| 久久久久久久午夜电影| 中文资源天堂在线| 中国美女看黄片| 欧美+亚洲+日韩+国产| 欧美另类亚洲清纯唯美| 国产精品一区二区三区四区久久| 欧美乱码精品一区二区三区|