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

    基于雙重Bregman迭代的地震數(shù)據(jù)重構(gòu)與去噪

    2020-10-11 07:25:40張會(huì)星劉明珠
    石油物探 2020年5期
    關(guān)鍵詞:雙重插值信噪比

    郭 萌,張會(huì)星,劉明珠

    (1.中國(guó)海洋大學(xué)海底科學(xué)與探測(cè)技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,山東青島266100;2.青島海洋科學(xué)與技術(shù)國(guó)家實(shí)驗(yàn)室海洋礦產(chǎn)資源評(píng)價(jià)與探測(cè)技術(shù)功能實(shí)驗(yàn)室,山東青島266071)

    陸上地震勘探受地形、障礙物、禁采區(qū)等采集環(huán)境以及采集成本、采集設(shè)備等的影響,采集到的地震數(shù)據(jù)常常出現(xiàn)道間距過(guò)大、壞道、部分偏移距數(shù)據(jù)缺失等問題,使得數(shù)據(jù)連續(xù)性變差,用這種不連續(xù)的地震數(shù)據(jù)進(jìn)行成像處理時(shí)往往會(huì)產(chǎn)生空間假頻或虛假繞射,從而降低了地震剖面的分辨率[1]。海上拖纜和海底地震勘探受成本、設(shè)備與洋流等影響,采集到的數(shù)據(jù)也常常出現(xiàn)空間方向的不規(guī)則、道距稀疏等現(xiàn)象,不規(guī)則的稀疏數(shù)據(jù)同樣會(huì)在偏移剖面中產(chǎn)生空間假頻、繞射噪聲或虛假同相軸等問題[2],造成這些問題的原因在于現(xiàn)有的許多處理技術(shù)都假設(shè)采集到的數(shù)據(jù)是完整連續(xù)的,當(dāng)實(shí)際輸入數(shù)據(jù)不滿足這一假設(shè)條件時(shí)必然會(huì)產(chǎn)生誤差,因此需要在成像處理之前采用合適的技術(shù)對(duì)這種缺道地震數(shù)據(jù)進(jìn)行重構(gòu),使其滿足后續(xù)成像處理所需的假設(shè)條件。

    地震數(shù)據(jù)插值重構(gòu)方法的研究始于20世紀(jì)80年代,LARNER等[3]基于同相軸傾角一致的假設(shè),給出了一種地震數(shù)據(jù)分時(shí)窗傾角掃描并內(nèi)插的方法。之后業(yè)界先后發(fā)展出一系列地震數(shù)據(jù)插值技術(shù),主要分為以下4類,第一類為基于預(yù)測(cè)濾波的方法,如SPITZ[4]提出的f-x域缺失地震數(shù)據(jù)插值恢復(fù)方法,該方法能夠一定程度上克服空間假頻,但只對(duì)規(guī)則采樣且當(dāng)?shù)卣鹩涗浀耐噍S為線性時(shí)的數(shù)據(jù)有用。國(guó)九英等[5]在f-x域預(yù)測(cè)濾波的基礎(chǔ)上將其拓展到f-k域,提高了插值計(jì)算效率。隨后GULUNAY等[6]提出抗假頻f-k域缺失地震道插值方法,并對(duì)該方法進(jìn)行擴(kuò)展,提出適用性更廣泛的廣義f-k域缺失地震道插值方法[7]。基于預(yù)測(cè)濾波的重構(gòu)方法可以加密空間采樣率,防止空間假頻,但要求輸入數(shù)據(jù)必須規(guī)則采樣。第二類為基于稀疏變換的方法,如傅里葉變換[2]、Radon變換[8]、曲波變換[9-10]等,這類方法對(duì)于非規(guī)則缺失數(shù)據(jù)插值效果較好,但該類方法大都需要人為設(shè)定閾值,插值效果依賴于處理人員的經(jīng)驗(yàn)。第三類為矩陣降秩的方法[11],該類方法將高維數(shù)據(jù)通過(guò)低秩矩陣進(jìn)行表示然后重建,計(jì)算簡(jiǎn)單快速,并發(fā)展到了五維地震數(shù)據(jù)重構(gòu)[12],但該方法的計(jì)算量會(huì)隨著地震數(shù)據(jù)維度的增加而成倍地增長(zhǎng),對(duì)硬件要求很高。第四類為基于波動(dòng)方程的方法,FOMEL[13]基于有限差分法進(jìn)行波動(dòng)方程缺失數(shù)據(jù)重構(gòu),LIU等[14]基于快速傅里葉變換的共軛梯度并結(jié)合快速矢量矩陣乘法的數(shù)值算法進(jìn)行了空間帶限信號(hào)的波場(chǎng)重建。這類方法需要速度信息,并且計(jì)算量大,因而沒有得到廣泛應(yīng)用,而且該類方法主要是對(duì)滿足傳統(tǒng)Nyquist采樣定律的地震數(shù)據(jù)進(jìn)行插值加密處理,對(duì)于不滿足傳統(tǒng)Nyquist采樣定理的地震數(shù)據(jù)插值效果不甚理想。

    壓縮感知理論的提出,給信號(hào)恢復(fù)和重構(gòu)帶來(lái)新的思路,地球物理學(xué)者將其引入地震數(shù)據(jù)恢復(fù)領(lǐng)域[15-17]。該理論以信號(hào)的稀疏性為前提,可以對(duì)遠(yuǎn)低于Nyquist標(biāo)準(zhǔn)的欠采樣數(shù)據(jù)進(jìn)行準(zhǔn)確恢復(fù),打破了Nyquist采樣定理的約束。FOMEL等[18]和LIU等[19]基于該理論提出了適用于地震數(shù)據(jù)的Seislet變換,實(shí)現(xiàn)了地震數(shù)據(jù)插值重構(gòu)。此外基于壓縮感知的Bregman迭代[20]、迭代閾值法[21]等地震數(shù)據(jù)重構(gòu)方法得到不斷研究和應(yīng)用。

    由于采集環(huán)境等因素的影響,野外地震數(shù)據(jù)往往存在隨機(jī)噪聲干擾,這種噪聲使得缺失數(shù)據(jù)得不到充分重建,或重構(gòu)精度較低,因此在重構(gòu)過(guò)程中需要對(duì)隨機(jī)噪聲進(jìn)行適當(dāng)壓制[22-23]。ZHANG等[24]利用二維曲波變換對(duì)地震數(shù)據(jù)同時(shí)進(jìn)行重構(gòu)和去噪,取得了較好的效果。本文基于壓縮感知理論,提出雙重Bregman迭代算法,即將分裂Bregman迭代作為內(nèi)部迭代用于壓制隨機(jī)噪聲,將線性Bregman迭代作為外部迭代用于數(shù)據(jù)重構(gòu),更高效地實(shí)現(xiàn)含噪缺失地震數(shù)據(jù)的同步重構(gòu)與去噪。

    1 地震數(shù)據(jù)重構(gòu)與去噪的數(shù)學(xué)反問題

    地震數(shù)據(jù)重建可以歸納為一類線性估計(jì)問題:假設(shè)需重建的地震數(shù)據(jù)為m,實(shí)際觀測(cè)到的數(shù)據(jù)為d,則d和m之間的關(guān)系可表示為:

    d=Lm

    (1)

    式中:L為空間采樣算子。由觀測(cè)到的已知信號(hào)d求解需重建數(shù)據(jù)m的過(guò)程就是地震數(shù)據(jù)重構(gòu)過(guò)程。

    FOMEL[25]指出:對(duì)缺失地震數(shù)據(jù)進(jìn)行重構(gòu),可設(shè)置一個(gè)掩膜矩陣K代替空間采樣算子L,則公式(1)可表示為:

    d=Km

    (2)

    掩膜矩陣K為對(duì)角矩陣,K中待重構(gòu)位置處的元素為0,其它位置處的元素為1。

    本文目標(biāo)是利用壓縮感知理論由已知的觀測(cè)數(shù)據(jù)d求取未知的完整數(shù)據(jù)m。由于壓縮感知理論假設(shè)信號(hào)m具有稀疏特性,因此設(shè)計(jì)一個(gè)稀疏基Φ對(duì)m進(jìn)行稀疏表示:

    m=Φ-1x

    (3)

    式中:x為m的變換域系數(shù)。將公式(3)代入公式(2)可得:

    d=KΦ-1x=Ax

    (4)

    式中:A=KΦ-1為感知矩陣,掩膜矩陣K與稀疏矩陣Φ滿足等距約束性準(zhǔn)則[26]。

    由此,求解完整數(shù)據(jù)m的過(guò)程轉(zhuǎn)換為求解稀疏信號(hào)x的過(guò)程,根據(jù)壓縮感知理論,將公式(4)轉(zhuǎn)化成L0范數(shù)問題進(jìn)行求解:

    (5)

    式中:‖·‖0為L(zhǎng)0范數(shù)。

    L0范數(shù)的求解較困難,利用L1范數(shù)代替L0范數(shù),將L0范數(shù)最優(yōu)化問題轉(zhuǎn)化為L(zhǎng)1范數(shù)的凸優(yōu)化問題:

    (6)

    構(gòu)造懲罰函數(shù)將公式(6)統(tǒng)一轉(zhuǎn)化成下列形式:

    (7)

    (8)

    式中:μ和λ為罰參數(shù),用于平衡正則項(xiàng)和保真項(xiàng)。

    對(duì)(7)式或(8)式進(jìn)行求解可獲得稀疏信號(hào)x,進(jìn)而通過(guò)公式(3)得到完整地震數(shù)據(jù)m。如果利用單位矩陣E代替掩膜矩陣K,則地震數(shù)據(jù)重構(gòu)問題便轉(zhuǎn)換成了地震數(shù)據(jù)去噪問題。問題的關(guān)鍵是如何準(zhǔn)確求解公式(7)或公式(8)。

    2 基于雙重Bregman迭代的地震數(shù)據(jù)重構(gòu)與去噪

    2.1 Bregman迭代重構(gòu)算法

    對(duì)于公式(7)或公式(8)的求解,YIN等[27]以Bregman迭代正則化算法為基礎(chǔ),結(jié)合不動(dòng)點(diǎn)延續(xù)(fixed-point continuation,FPC)算法形成了一種新的Bregman迭代框架,即:

    (9)

    (10)

    (11)

    i=1,…,n

    其中,閾值算子shrink的定義為:設(shè)實(shí)數(shù)空間R中,y∈R,α∈R+,則

    (12)

    為提高公式(10)的求解效率,使用FPC方法求解。根據(jù)文獻(xiàn)[28]可知,FPC算法具體迭代流程如表1所示。

    表1 FPC算法實(shí)現(xiàn)步驟

    YIN等[27]的研究結(jié)果表明,由公式(9)迭代算法重構(gòu)出的數(shù)據(jù)更精確,圖像對(duì)比度更高,且由于使用了FPC算法,計(jì)算速度更快,因此本文采用公式(9)進(jìn)行缺失地震數(shù)據(jù)插值處理。

    表1迭代過(guò)程中,迭代次數(shù)k是決定能否得到最優(yōu)值的關(guān)鍵參數(shù),k過(guò)小無(wú)法收斂到最優(yōu)解,k過(guò)大則會(huì)降低計(jì)算效率,因此需要通過(guò)設(shè)定迭代停止準(zhǔn)則來(lái)控制迭代次數(shù),本文采用的迭代停止準(zhǔn)則為:

    (13)

    2.2 分裂Bregman迭代去噪算法

    分裂Bregman迭代算法由GOLDSTEIN等[29]提出,本文基于Rudin-Osher-Fatemi(ROF)模型[29]推導(dǎo)分裂Bregman迭代去噪框架。ROF模型利用數(shù)據(jù)的全變分將圖像降噪問題轉(zhuǎn)化為求解最優(yōu)化問題,對(duì)于高斯噪聲,具體的ROF模型表示形式為:

    (14)

    (15)

    參照公式(8)加入罰函數(shù)λ1,可得:

    (16)

    (17)

    可以看出公式(17)與公式(8)形式相同。參照上述Bregman的迭代步驟,利用公式(9)即可將公式(16)轉(zhuǎn)化為以下多變量的Bregman迭代形式:

    (18)

    上述公式含有兩個(gè)變量x,D,根據(jù)文獻(xiàn)[29],為獲得公式(18)的最優(yōu)解,解耦公式(18)中的L1部分和L2部分,分別對(duì)D和x進(jìn)行迭代最小化,迭代過(guò)程如下。

    第一步:固定Dk,更新xk。公式為:

    (19)

    第二步:固定xk,更新Dk。公式為:

    (20)

    第三步:更新bk。公式為:

    (21)

    (22)

    分裂Bregman迭代的迭代停止準(zhǔn)則為:

    (23)

    本文采用傅里葉變換對(duì)地震數(shù)據(jù)進(jìn)行稀疏表示,采用掩膜算子作為觀測(cè)矩陣,以YIN等[27]提出的Bregman迭代作為外部迭代,GOLDSTEIN等[29]提出的分裂Bregman迭代作為內(nèi)部迭代,形成雙重Bregman迭代,在滿足迭代停止準(zhǔn)則的條件下完成對(duì)含噪聲的缺失地震記錄進(jìn)行重構(gòu)和去噪。具體實(shí)現(xiàn)流程見圖1。

    圖1 雙重Bregman迭代流程

    3 模型數(shù)據(jù)測(cè)試

    為驗(yàn)證雙重Bregman算法的迭代效果,選擇Marmousi部分模型(圖2)正演數(shù)據(jù)的縱波分量進(jìn)行算法測(cè)試。正演記錄由彈性波方程有限差分模擬得到,模擬所用的參數(shù)如下:地表激發(fā)和接收,激發(fā)震源為主頻30Hz的Ricker子波,接收道為601道,道間距為5m,采樣間隔為0.35ms,記錄長(zhǎng)度為2100ms。圖3a為模擬單炮記錄,圖3b為加入信噪比為3.46dB的高斯隨機(jī)噪聲后的單炮記錄,隨機(jī)抽取30%的地震道充零(圖3c),對(duì)充零后的地震記錄進(jìn)行雙重Bregman迭代。根據(jù)迭代準(zhǔn)則(13)式和(23)式,設(shè)置內(nèi)部迭代允許最大迭代誤差為10-3,外部迭代允許最大誤差為10-5;為避免出現(xiàn)程序無(wú)限迭代不能正常退出的情況,經(jīng)多次試驗(yàn),同時(shí)限定了內(nèi)部和外部迭代的最大迭代次數(shù)為20次。計(jì)算得到的重構(gòu)結(jié)果如圖3d所示,圖3e所示為算法誤差剖面,是圖3a和圖3d相減后的結(jié)果。結(jié)合圖3d和圖3e可以看出,對(duì)含隨機(jī)噪聲的缺失地震數(shù)據(jù),采用本文算法實(shí)現(xiàn)了地震數(shù)據(jù)的重構(gòu)與去噪同步進(jìn)行,可以有效壓制隨機(jī)噪聲,并同時(shí)獲得完整的重構(gòu)數(shù)據(jù),提高了信噪比。

    圖2 Marmousi部分模型

    利用本文方法對(duì)不同缺失程度下地震數(shù)據(jù)的重構(gòu)與去噪效果以及該方法計(jì)算誤差剖面如圖4所示。其中,圖4a和圖4d分別為對(duì)信噪比為3.46dB的Marmousi地震數(shù)據(jù)進(jìn)行50%地震道充零和70%地震道充零,然后進(jìn)行與30%地震道充零相同的雙重Bregman迭代重構(gòu)與去噪處理,分別得到如圖4b和圖4e所示的重構(gòu)與去噪結(jié)果,以及圖4c和圖4f的誤差剖面。圖5為3種情況下的誤差曲線。另外,為對(duì)比本文方法與線性Bregman迭代的重構(gòu)效果與計(jì)算效率,對(duì)上述3種地震數(shù)據(jù)同時(shí)進(jìn)行了線性Bregman迭代計(jì)算。計(jì)算環(huán)境參數(shù)為:CPU 2.60GHz,RAM 4.0GB,64位Windows 7系統(tǒng)。重構(gòu)和去噪時(shí)間以及處理后的信噪比量化評(píng)價(jià)結(jié)果見表2。綜合圖3、圖4、圖5以及表2,可得如下結(jié)論:當(dāng)欠采樣數(shù)據(jù)為完整數(shù)據(jù)的70%時(shí),利用本文方法能夠得到相對(duì)完整的重構(gòu)數(shù)據(jù),處理后的信噪比為8.73,得到了大幅度提升;當(dāng)欠采樣數(shù)據(jù)為完整數(shù)據(jù)的50%時(shí),本文方法基本上實(shí)現(xiàn)了對(duì)缺失數(shù)據(jù)的重構(gòu),并且重構(gòu)與處理后的信噪比由原來(lái)的3.46提升到4.46;當(dāng)欠采樣數(shù)據(jù)為完整數(shù)據(jù)的30%時(shí),本文方法能夠較完整地重建出缺失數(shù)據(jù),但仍有部分?jǐn)?shù)據(jù)未能得到準(zhǔn)確重構(gòu),而且可能由于缺失程度較大、迭代過(guò)程中引入重構(gòu)噪聲等因素,迭代處理后的數(shù)據(jù)信噪比有所降低。經(jīng)測(cè)試發(fā)現(xiàn),當(dāng)欠采樣數(shù)據(jù)為完整數(shù)據(jù)的一半及以上時(shí),能夠利用本文方法得到相對(duì)完整的重構(gòu)與去噪成果,并且數(shù)據(jù)缺失越少,重構(gòu)和去噪所需計(jì)算時(shí)間越少,誤差曲線下降越快,信噪比提升也越高;在不同缺失程度下,雙重Bregman迭代與線性Bregman迭代的運(yùn)行時(shí)間相差不大。這是因?yàn)榉至袯regman迭代的速度相當(dāng)快。GOLDSTEIN等[29]曾指出,分裂Bregman迭代速度很快的主要原因是由其計(jì)算方法決定的,分裂Bregman迭代是將L1和L2部分進(jìn)行解耦,然后分別用快速收斂的Gauss-Seidel方法和閾值算子方法求解,這種算法加快了分裂Bregman迭代的速度;另外,分裂Bregman的迭代速度與參數(shù)選擇有很大關(guān)系,參數(shù)選擇得當(dāng),該算法非常快速。綜上,雙重Bregman迭代不僅有效去除了隨機(jī)噪聲,其結(jié)果的信噪比明顯高于線性Bregman迭代結(jié)果,而且運(yùn)算時(shí)間并沒有大量增加,因此,對(duì)于含有隨機(jī)噪聲的地震數(shù)據(jù),本文方法優(yōu)于線性Bregman迭代方法,具有一定的研究意義及應(yīng)用價(jià)值。

    圖3 Marmousi模型含噪聲的缺失地震記錄重構(gòu)與去噪

    圖4 不同缺失程度下雙重Bregman迭代對(duì)含噪聲地震數(shù)據(jù)的重構(gòu)與去噪結(jié)果和誤差剖面

    表2 雙重Bregman迭代和線性Bregman迭代在不同缺失程度下的重構(gòu)和去噪時(shí)間及信噪比

    圖5 不同缺失程度下含噪聲地震數(shù)據(jù)的重構(gòu)與去噪誤差曲線

    4 實(shí)際數(shù)據(jù)測(cè)試

    利用海上某工區(qū)實(shí)際地震數(shù)據(jù)進(jìn)行測(cè)試,原始地震數(shù)據(jù)共301道,每道10001個(gè)采樣點(diǎn),重構(gòu)工作在CDP道集(圖6a)上進(jìn)行,對(duì)該數(shù)據(jù)隨機(jī)抽取30%的地震道充零后的結(jié)果如圖6b,對(duì)抽稀后的地震數(shù)據(jù)進(jìn)行雙重Bregman迭代,迭代停止準(zhǔn)則采用公式(13)和公式(23),經(jīng)多次試驗(yàn),同時(shí)設(shè)置內(nèi)部迭代允許最大迭代次數(shù)為50次,外部迭代允許最大迭代次數(shù)為100次,得到的結(jié)果如圖6c。為說(shuō)明本文算法的重構(gòu)與去噪效果,對(duì)圖6b進(jìn)行線性Bregman迭代計(jì)算,迭代停止準(zhǔn)則為公式(13),允許最大迭代次數(shù)為100次,迭代結(jié)果如圖6d。結(jié)合圖6中藍(lán)色方框?qū)?yīng)的局部放大波形顯示(圖7)可以看出,兩種迭代算法都能夠?qū)θ笔?shù)據(jù)進(jìn)行較好恢復(fù),但本文迭代算法與原始數(shù)據(jù)的重合度更高。將兩種算法迭代后的地震數(shù)據(jù)與原始地震數(shù)據(jù)求差,其中1.5~2.5s時(shí)間段第100~150道(即藍(lán)色方框?qū)?yīng)范圍)的誤差如圖8所示,整體地震數(shù)據(jù)的重構(gòu)誤差如圖9所示,對(duì)比可以看出,雙重Bregman算法迭代誤差遠(yuǎn)遠(yuǎn)小于線性Bregman迭代,進(jìn)一步顯現(xiàn)出本文迭代算法的優(yōu)越性。

    圖6 實(shí)測(cè)地震記錄的重構(gòu)與去噪

    圖7 圖6藍(lán)色方框的放大顯示

    圖8 兩種Bregman迭代局部結(jié)果的精度比較

    圖9 兩種Bregman迭代整體結(jié)果的精度比較

    5 結(jié)論

    本文結(jié)合Bregman迭代重構(gòu)算法與分裂Bregman迭代去噪算法,提出了一種基于雙重Bregman迭代的地震數(shù)據(jù)重構(gòu)與去噪方法。Marmousi模型數(shù)值模擬實(shí)驗(yàn)和實(shí)際數(shù)據(jù)應(yīng)用結(jié)果表明:本文算法沒有因兩種獨(dú)立算法的結(jié)合而減弱對(duì)數(shù)據(jù)的處理效果,且只需較少次數(shù)的迭代就能高效實(shí)現(xiàn)含噪聲缺道地震數(shù)據(jù)的同時(shí)重構(gòu)與去噪,為地震數(shù)據(jù)恢復(fù)提供了一種新的缺失地震數(shù)據(jù)處理方法。

    基于壓縮感知的地震數(shù)據(jù)處理方法需要先將地震數(shù)據(jù)通過(guò)數(shù)學(xué)變換進(jìn)行稀疏表示,本文采用了最常用的傅里葉稀疏表示方法,而該方法只適用于同相軸近似線性或者平穩(wěn)變化的地震數(shù)據(jù),因此本文進(jìn)一步的工作是尋找復(fù)雜地質(zhì)構(gòu)造下地震數(shù)據(jù)的最佳稀疏表示方法,結(jié)合雙重Bregman迭代,進(jìn)一步研究復(fù)雜地質(zhì)構(gòu)造下含噪聲缺失地震數(shù)據(jù)的高精度重構(gòu)與去噪方法。

    猜你喜歡
    雙重插值信噪比
    自然與成長(zhǎng)的雙重變奏
    基于深度學(xué)習(xí)的無(wú)人機(jī)數(shù)據(jù)鏈信噪比估計(jì)算法
    化解“雙重目標(biāo)”之困
    基于Sinc插值與相關(guān)譜的縱橫波速度比掃描方法
    低信噪比下LFMCW信號(hào)調(diào)頻參數(shù)估計(jì)
    低信噪比下基于Hough變換的前視陣列SAR稀疏三維成像
    一種改進(jìn)FFT多譜線插值諧波分析方法
    基于四項(xiàng)最低旁瓣Nuttall窗的插值FFT諧波分析
    保持信噪比的相位分解反褶積方法研究
    “雙重打擊”致恐龍滅絕
    日本与韩国留学比较| 国产成人精品婷婷| 又大又黄又爽视频免费| 伊人亚洲综合成人网| 亚洲精品aⅴ在线观看| 中文天堂在线官网| 日韩人妻高清精品专区| 日本黄大片高清| 日本-黄色视频高清免费观看| 亚洲av二区三区四区| 国产精品.久久久| 久久久国产精品麻豆| 你懂的网址亚洲精品在线观看| 国产欧美日韩综合在线一区二区| 久久韩国三级中文字幕| 久久精品国产亚洲av天美| 不卡视频在线观看欧美| 26uuu在线亚洲综合色| 日韩强制内射视频| 国产国拍精品亚洲av在线观看| 精品亚洲成国产av| 国产在线一区二区三区精| 老熟女久久久| 久久久国产欧美日韩av| 嘟嘟电影网在线观看| 免费看av在线观看网站| 精品人妻熟女av久视频| 亚洲欧美一区二区三区国产| 狂野欧美白嫩少妇大欣赏| a级毛色黄片| 午夜激情福利司机影院| 国产精品99久久99久久久不卡 | 特大巨黑吊av在线直播| 亚洲精品乱码久久久v下载方式| 国产成人freesex在线| 国产精品偷伦视频观看了| a级毛片免费高清观看在线播放| 国产黄色视频一区二区在线观看| 国产不卡av网站在线观看| 亚洲第一区二区三区不卡| 女人精品久久久久毛片| 大香蕉97超碰在线| 少妇被粗大的猛进出69影院 | 久久人人爽av亚洲精品天堂| 亚洲人成网站在线播| 女人久久www免费人成看片| 婷婷色综合大香蕉| videossex国产| 亚洲国产精品国产精品| 久久久久久久久久久免费av| 母亲3免费完整高清在线观看 | 成人综合一区亚洲| 你懂的网址亚洲精品在线观看| 亚洲一级一片aⅴ在线观看| 精品久久久精品久久久| 午夜影院在线不卡| 日韩人妻高清精品专区| 亚洲精品久久午夜乱码| 亚洲不卡免费看| 看十八女毛片水多多多| 各种免费的搞黄视频| 最近手机中文字幕大全| 成年女人在线观看亚洲视频| 欧美日韩成人在线一区二区| 中文字幕av电影在线播放| 国产精品久久久久久久久免| av有码第一页| 国产欧美另类精品又又久久亚洲欧美| 人妻少妇偷人精品九色| 大又大粗又爽又黄少妇毛片口| 亚洲av日韩在线播放| 一级片'在线观看视频| av女优亚洲男人天堂| h视频一区二区三区| av网站免费在线观看视频| 亚洲内射少妇av| 97精品久久久久久久久久精品| 赤兔流量卡办理| 2022亚洲国产成人精品| 91在线精品国自产拍蜜月| 欧美成人精品欧美一级黄| 伊人久久精品亚洲午夜| 各种免费的搞黄视频| 在线 av 中文字幕| 中国美白少妇内射xxxbb| 日本av免费视频播放| 精品少妇久久久久久888优播| 99热全是精品| 国产高清三级在线| 丝袜喷水一区| 亚洲av福利一区| 日本猛色少妇xxxxx猛交久久| 黑人猛操日本美女一级片| 午夜久久久在线观看| 91久久精品国产一区二区三区| 最近中文字幕高清免费大全6| 国产高清有码在线观看视频| 免费日韩欧美在线观看| 久久精品久久久久久久性| 高清在线视频一区二区三区| 中文字幕久久专区| 日韩在线高清观看一区二区三区| 国产免费现黄频在线看| 伊人亚洲综合成人网| 最新中文字幕久久久久| 日本爱情动作片www.在线观看| av女优亚洲男人天堂| 日韩一区二区视频免费看| 日韩一本色道免费dvd| 亚洲情色 制服丝袜| 国产精品蜜桃在线观看| 制服丝袜香蕉在线| 18禁在线播放成人免费| 高清在线视频一区二区三区| 亚洲精品自拍成人| 啦啦啦视频在线资源免费观看| 熟妇人妻不卡中文字幕| 人妻夜夜爽99麻豆av| 成人影院久久| 亚洲精品乱码久久久v下载方式| 免费观看av网站的网址| 国产片内射在线| 久久精品久久久久久久性| 精品酒店卫生间| 99久久中文字幕三级久久日本| av女优亚洲男人天堂| 久久亚洲国产成人精品v| 亚洲精品日韩av片在线观看| 久久国内精品自在自线图片| 91成人精品电影| 日韩一区二区视频免费看| 黑人高潮一二区| 日韩中文字幕视频在线看片| 欧美日韩视频高清一区二区三区二| 男男h啪啪无遮挡| 免费观看的影片在线观看| 久热久热在线精品观看| 国产一级毛片在线| 国产色婷婷99| 成人亚洲精品一区在线观看| 国产精品偷伦视频观看了| 一本—道久久a久久精品蜜桃钙片| 国产精品久久久久久久久免| 日韩电影二区| 国产成人精品福利久久| 久久久久久久久久久久大奶| 一级毛片 在线播放| 国产成人91sexporn| 啦啦啦视频在线资源免费观看| 日韩 亚洲 欧美在线| 欧美精品高潮呻吟av久久| 一级毛片黄色毛片免费观看视频| 69精品国产乱码久久久| 丰满少妇做爰视频| 18+在线观看网站| 亚洲一区二区三区欧美精品| 激情五月婷婷亚洲| 亚洲不卡免费看| 伦理电影大哥的女人| 久久99精品国语久久久| 国产成人a∨麻豆精品| 国产成人精品一,二区| 国产黄色免费在线视频| 插阴视频在线观看视频| 国产精品不卡视频一区二区| 亚洲国产av新网站| 99热这里只有是精品在线观看| 人人妻人人澡人人爽人人夜夜| av又黄又爽大尺度在线免费看| 青春草国产在线视频| 人人妻人人澡人人爽人人夜夜| 日本猛色少妇xxxxx猛交久久| 蜜桃国产av成人99| 亚洲欧美精品自产自拍| av国产精品久久久久影院| tube8黄色片| 国产极品天堂在线| 夜夜骑夜夜射夜夜干| 日韩一区二区视频免费看| 69精品国产乱码久久久| 97精品久久久久久久久久精品| 国产在线免费精品| 精品久久蜜臀av无| 国产一级毛片在线| 亚洲欧美一区二区三区黑人 | 亚洲精品乱码久久久久久按摩| 日韩成人av中文字幕在线观看| 新久久久久国产一级毛片| 成年美女黄网站色视频大全免费 | 久久精品人人爽人人爽视色| 丰满迷人的少妇在线观看| 亚洲国产精品国产精品| 亚洲美女黄色视频免费看| 丰满少妇做爰视频| 人妻制服诱惑在线中文字幕| 免费久久久久久久精品成人欧美视频 | 亚洲美女黄色视频免费看| 99九九在线精品视频| 成年女人在线观看亚洲视频| 国产成人免费观看mmmm| 国产69精品久久久久777片| 看十八女毛片水多多多| 日韩强制内射视频| 亚洲三级黄色毛片| 国产伦精品一区二区三区视频9| 国产精品99久久99久久久不卡 | 美女国产高潮福利片在线看| 老司机影院毛片| 国产亚洲午夜精品一区二区久久| 91成人精品电影| 国产精品偷伦视频观看了| videossex国产| 熟女av电影| 亚洲伊人久久精品综合| 狂野欧美白嫩少妇大欣赏| 亚洲经典国产精华液单| 中文字幕最新亚洲高清| 性色av一级| 免费黄网站久久成人精品| 欧美日韩视频高清一区二区三区二| 亚洲四区av| 久久久国产欧美日韩av| 你懂的网址亚洲精品在线观看| 99国产综合亚洲精品| 精品一区二区免费观看| 久久97久久精品| 蜜桃国产av成人99| 91精品三级在线观看| 这个男人来自地球电影免费观看 | 啦啦啦中文免费视频观看日本| 妹子高潮喷水视频| 日韩在线高清观看一区二区三区| 亚洲av国产av综合av卡| 成人18禁高潮啪啪吃奶动态图 | 少妇丰满av| 少妇人妻精品综合一区二区| 亚洲成色77777| 制服人妻中文乱码| 久久这里有精品视频免费| 五月伊人婷婷丁香| 一边摸一边做爽爽视频免费| 欧美精品人与动牲交sv欧美| 一本一本综合久久| 美女脱内裤让男人舔精品视频| 久久久久人妻精品一区果冻| 久久狼人影院| 99久久中文字幕三级久久日本| 在线观看www视频免费| 亚洲国产精品一区二区三区在线| 王馨瑶露胸无遮挡在线观看| 99国产精品免费福利视频| 精品午夜福利在线看| 亚洲人成网站在线播| 免费高清在线观看视频在线观看| 午夜激情av网站| 大陆偷拍与自拍| 综合色丁香网| 免费看不卡的av| 国产精品人妻久久久影院| 特大巨黑吊av在线直播| 制服诱惑二区| 日韩一区二区三区影片| 天天躁夜夜躁狠狠久久av| 性色avwww在线观看| 亚洲,欧美,日韩| 久久人人爽人人爽人人片va| 亚洲少妇的诱惑av| 最新中文字幕久久久久| 18禁动态无遮挡网站| 男女边摸边吃奶| 精品少妇内射三级| 一个人看视频在线观看www免费| 色哟哟·www| av黄色大香蕉| 国产一区二区在线观看日韩| 国产一级毛片在线| 午夜av观看不卡| 国产精品99久久99久久久不卡 | 性色av一级| 伊人亚洲综合成人网| 一本久久精品| 国产白丝娇喘喷水9色精品| 九九久久精品国产亚洲av麻豆| 99热这里只有是精品在线观看| 天天操日日干夜夜撸| av电影中文网址| 成年人免费黄色播放视频| 伦理电影免费视频| 少妇的逼水好多| 777米奇影视久久| 高清在线视频一区二区三区| 日韩制服骚丝袜av| 九九爱精品视频在线观看| 免费人妻精品一区二区三区视频| 狂野欧美白嫩少妇大欣赏| 寂寞人妻少妇视频99o| 国模一区二区三区四区视频| 久久久久久久久久久免费av| 国产精品久久久久久av不卡| 国产精品.久久久| 校园人妻丝袜中文字幕| 精品人妻熟女毛片av久久网站| 97在线视频观看| 国产免费又黄又爽又色| 国产男人的电影天堂91| 少妇精品久久久久久久| 精品久久久久久电影网| 日日啪夜夜爽| 极品人妻少妇av视频| 中文字幕制服av| 国产午夜精品久久久久久一区二区三区| 午夜激情av网站| 国产成人精品久久久久久| 男男h啪啪无遮挡| 亚洲精品美女久久av网站| 久久影院123| 秋霞伦理黄片| 九色成人免费人妻av| 一级爰片在线观看| 七月丁香在线播放| 丝袜喷水一区| 五月玫瑰六月丁香| 中文字幕久久专区| 有码 亚洲区| 亚洲三级黄色毛片| 免费久久久久久久精品成人欧美视频 | 最近2019中文字幕mv第一页| 一级,二级,三级黄色视频| a级毛片黄视频| 日本免费在线观看一区| 国产精品熟女久久久久浪| 夜夜看夜夜爽夜夜摸| 成人无遮挡网站| 国产精品一区二区三区四区免费观看| 久久韩国三级中文字幕| 夫妻性生交免费视频一级片| 王馨瑶露胸无遮挡在线观看| 国产成人精品在线电影| 午夜老司机福利剧场| videosex国产| 嫩草影院入口| 国产av码专区亚洲av| videossex国产| 国产男女超爽视频在线观看| 一区二区三区免费毛片| 国产精品久久久久久久电影| 国产精品嫩草影院av在线观看| 一级毛片黄色毛片免费观看视频| 久久国内精品自在自线图片| 国产免费一级a男人的天堂| 亚洲欧美清纯卡通| 欧美xxⅹ黑人| 国产精品一区www在线观看| 天天操日日干夜夜撸| 精品一区二区三区视频在线| 一级,二级,三级黄色视频| 九色亚洲精品在线播放| 免费高清在线观看视频在线观看| 妹子高潮喷水视频| 看非洲黑人一级黄片| 特大巨黑吊av在线直播| 亚洲怡红院男人天堂| 国产精品蜜桃在线观看| 少妇人妻 视频| 亚洲精品国产av蜜桃| 日韩中字成人| 亚洲国产精品国产精品| 黑人高潮一二区| 天天影视国产精品| 日韩av免费高清视频| 美女中出高潮动态图| xxx大片免费视频| av在线老鸭窝| 免费久久久久久久精品成人欧美视频 | 在线免费观看不下载黄p国产| 大陆偷拍与自拍| 国产淫语在线视频| 久久久国产精品麻豆| 国产男女内射视频| 久久久久久久大尺度免费视频| 成人亚洲精品一区在线观看| 国产国拍精品亚洲av在线观看| 亚洲精品国产av蜜桃| 男女边摸边吃奶| 能在线免费看毛片的网站| 熟女av电影| 草草在线视频免费看| 热99国产精品久久久久久7| 男女无遮挡免费网站观看| 制服诱惑二区| 久久久久网色| 人妻 亚洲 视频| 亚洲精品一区蜜桃| 久久热精品热| 亚洲高清免费不卡视频| av视频免费观看在线观看| www.色视频.com| 水蜜桃什么品种好| 精品少妇黑人巨大在线播放| 日韩三级伦理在线观看| 亚洲成人手机| 亚洲精品乱码久久久v下载方式| av网站免费在线观看视频| 毛片一级片免费看久久久久| 午夜免费男女啪啪视频观看| 亚洲av.av天堂| 另类精品久久| 久久人妻熟女aⅴ| 最近中文字幕高清免费大全6| 中文字幕人妻熟人妻熟丝袜美| 最新的欧美精品一区二区| 精品卡一卡二卡四卡免费| 久久亚洲国产成人精品v| 日韩av免费高清视频| 80岁老熟妇乱子伦牲交| 综合色丁香网| 亚洲人成77777在线视频| 我要看黄色一级片免费的| 国产免费福利视频在线观看| 一级黄片播放器| 欧美精品一区二区免费开放| 久久精品国产亚洲av天美| 一本色道久久久久久精品综合| 日本欧美视频一区| 国产片内射在线| 伦精品一区二区三区| 亚洲精品av麻豆狂野| 久久国产精品大桥未久av| 亚洲四区av| 一级,二级,三级黄色视频| 久久精品久久久久久久性| 精品国产乱码久久久久久小说| 伊人久久国产一区二区| 精品久久蜜臀av无| 18+在线观看网站| 色视频在线一区二区三区| 亚洲中文av在线| 性色avwww在线观看| 三级国产精品欧美在线观看| 午夜激情福利司机影院| 最黄视频免费看| 99久久人妻综合| 一本色道久久久久久精品综合| 亚洲精品久久成人aⅴ小说 | 亚洲国产成人一精品久久久| 我的女老师完整版在线观看| 人妻 亚洲 视频| 黄片播放在线免费| 国产白丝娇喘喷水9色精品| 婷婷色综合www| 久久免费观看电影| 亚洲精品456在线播放app| 99热国产这里只有精品6| 老女人水多毛片| av免费在线看不卡| 久久久久久人妻| 久久免费观看电影| 亚洲欧美一区二区三区黑人 | 亚洲久久久国产精品| 日韩av免费高清视频| 亚洲人成网站在线播| 亚洲精品国产av成人精品| 欧美一级a爱片免费观看看| www.色视频.com| 亚洲av不卡在线观看| 国产成人精品无人区| 最近手机中文字幕大全| 久久久久国产网址| 久久影院123| 亚洲av成人精品一区久久| 欧美人与性动交α欧美精品济南到 | 最近最新中文字幕免费大全7| 国产精品久久久久久精品古装| 亚洲美女黄色视频免费看| 精品99又大又爽又粗少妇毛片| 日韩三级伦理在线观看| 中国美白少妇内射xxxbb| 热re99久久国产66热| 人妻夜夜爽99麻豆av| 国产亚洲一区二区精品| 久久久久久久久大av| 精品人妻在线不人妻| 免费观看性生交大片5| 全区人妻精品视频| 久久精品夜色国产| 国产一区二区在线观看日韩| 黄色欧美视频在线观看| 天堂8中文在线网| 亚洲成色77777| 丝袜脚勾引网站| 91aial.com中文字幕在线观看| 97在线人人人人妻| 97在线视频观看| 国产在线免费精品| 亚洲国产精品成人久久小说| 精品久久久久久久久亚洲| 又黄又爽又刺激的免费视频.| 国产精品一区二区在线不卡| 久久国产精品大桥未久av| 国产成人aa在线观看| 久久精品国产自在天天线| 夜夜爽夜夜爽视频| 久久久久久久精品精品| 成年美女黄网站色视频大全免费 | 日本vs欧美在线观看视频| 在线观看免费视频网站a站| 国产综合精华液| 久久国产精品大桥未久av| 亚洲欧美一区二区三区国产| 亚洲欧洲国产日韩| 亚洲丝袜综合中文字幕| 91在线精品国自产拍蜜月| 亚洲欧美精品自产自拍| 亚洲精品中文字幕在线视频| 最近最新中文字幕免费大全7| 在线播放无遮挡| 成人无遮挡网站| 婷婷色麻豆天堂久久| 老熟女久久久| 日韩强制内射视频| 99精国产麻豆久久婷婷| a级毛色黄片| 日日撸夜夜添| 色5月婷婷丁香| 日韩欧美精品免费久久| 免费高清在线观看视频在线观看| 国产亚洲精品第一综合不卡 | 亚洲av综合色区一区| 一本色道久久久久久精品综合| 尾随美女入室| 99久久精品国产国产毛片| 91aial.com中文字幕在线观看| av在线老鸭窝| 日韩成人av中文字幕在线观看| 久久人妻熟女aⅴ| 婷婷色av中文字幕| 国产成人精品婷婷| 岛国毛片在线播放| 亚洲国产精品成人久久小说| 最新中文字幕久久久久| 考比视频在线观看| 欧美日韩国产mv在线观看视频| 汤姆久久久久久久影院中文字幕| 精品亚洲成国产av| 观看美女的网站| 亚洲经典国产精华液单| 日韩 亚洲 欧美在线| 国产伦精品一区二区三区视频9| 国产午夜精品一二区理论片| 午夜91福利影院| freevideosex欧美| 免费黄频网站在线观看国产| 观看美女的网站| 亚洲av日韩在线播放| 中文字幕制服av| 国国产精品蜜臀av免费| 99精国产麻豆久久婷婷| 另类精品久久| 一边亲一边摸免费视频| 性高湖久久久久久久久免费观看| 亚洲精品第二区| 最近最新中文字幕免费大全7| 久久久精品免费免费高清| 青春草亚洲视频在线观看| 日韩一区二区视频免费看| 欧美最新免费一区二区三区| 黑人欧美特级aaaaaa片| 不卡视频在线观看欧美| 91精品国产国语对白视频| 老司机影院毛片| 美女cb高潮喷水在线观看| 午夜福利在线观看免费完整高清在| 久久久久网色| 亚洲精品自拍成人| 美女视频免费永久观看网站| 少妇丰满av| 精品人妻一区二区三区麻豆| 午夜激情久久久久久久| 国产不卡av网站在线观看| 国产高清有码在线观看视频| 国产精品一区二区在线观看99| 99久久精品国产国产毛片| 亚洲激情五月婷婷啪啪| 丝袜在线中文字幕| www.色视频.com| 久久精品熟女亚洲av麻豆精品| 国产亚洲午夜精品一区二区久久| 国产成人一区二区在线| 丝袜美足系列| 国产女主播在线喷水免费视频网站| 国产乱来视频区| 久久精品久久久久久噜噜老黄| 国产精品麻豆人妻色哟哟久久| 青春草亚洲视频在线观看| 性高湖久久久久久久久免费观看| 美女cb高潮喷水在线观看| 久久精品国产自在天天线| 国产不卡av网站在线观看| 国产亚洲精品久久久com| 久久精品久久久久久噜噜老黄| 97超碰精品成人国产| 中文字幕精品免费在线观看视频 | 欧美少妇被猛烈插入视频| 色网站视频免费| 国产精品久久久久久久久免| 男女国产视频网站| 天堂中文最新版在线下载| 免费黄网站久久成人精品| 美女xxoo啪啪120秒动态图| 久久99一区二区三区| 九九爱精品视频在线观看| 久久人人爽av亚洲精品天堂| 欧美精品人与动牲交sv欧美| 国产亚洲精品第一综合不卡 | 亚洲不卡免费看|