• <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諧波分析
    保持信噪比的相位分解反褶積方法研究
    “雙重打擊”致恐龍滅絕
    亚洲精品国产av蜜桃| 久久青草综合色| 一级爰片在线观看| 春色校园在线视频观看| 日本猛色少妇xxxxx猛交久久| 亚洲国产精品专区欧美| 成人毛片a级毛片在线播放| 亚洲国产毛片av蜜桃av| 亚洲图色成人| 最近中文字幕高清免费大全6| 国产国语露脸激情在线看| 日韩人妻精品一区2区三区| 熟女av电影| 久久久a久久爽久久v久久| 欧美丝袜亚洲另类| 五月天丁香电影| 国产综合精华液| 各种免费的搞黄视频| 这个男人来自地球电影免费观看 | 69精品国产乱码久久久| 午夜日本视频在线| 日韩一本色道免费dvd| 国产伦理片在线播放av一区| 最后的刺客免费高清国语| 国产白丝娇喘喷水9色精品| 久久国产亚洲av麻豆专区| 草草在线视频免费看| 丝袜人妻中文字幕| 午夜福利网站1000一区二区三区| 亚洲久久久国产精品| 亚洲国产精品国产精品| 国产高清国产精品国产三级| 国产精品一国产av| 成人亚洲精品一区在线观看| 国产精品麻豆人妻色哟哟久久| 人妻少妇偷人精品九色| 9色porny在线观看| 亚洲av男天堂| 一级,二级,三级黄色视频| 国产成人精品在线电影| 久久鲁丝午夜福利片| 一个人免费看片子| 国产无遮挡羞羞视频在线观看| 国产成人av激情在线播放| 国产免费福利视频在线观看| 在线亚洲精品国产二区图片欧美| 国产熟女欧美一区二区| 久久精品久久精品一区二区三区| 亚洲国产精品999| 亚洲一级一片aⅴ在线观看| 国产有黄有色有爽视频| 欧美少妇被猛烈插入视频| 赤兔流量卡办理| 亚洲综合色网址| 尾随美女入室| 校园人妻丝袜中文字幕| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲天堂av无毛| 欧美97在线视频| 亚洲国产精品成人久久小说| 久久久久人妻精品一区果冻| 一区在线观看完整版| 中文字幕人妻丝袜制服| 欧美激情极品国产一区二区三区 | 最近最新中文字幕大全免费视频 | 91成人精品电影| 国产xxxxx性猛交| 少妇人妻精品综合一区二区| 欧美日韩视频精品一区| av福利片在线| 久久99蜜桃精品久久| 最近最新中文字幕大全免费视频 | 久久久精品区二区三区| 久久久久精品性色| 欧美性感艳星| 又粗又硬又长又爽又黄的视频| 午夜福利在线观看免费完整高清在| 日韩一区二区视频免费看| 亚洲精品成人av观看孕妇| 最近中文字幕高清免费大全6| 中文字幕av电影在线播放| 黄片播放在线免费| 色视频在线一区二区三区| 天美传媒精品一区二区| 欧美精品国产亚洲| 日本午夜av视频| 亚洲一码二码三码区别大吗| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 成人亚洲欧美一区二区av| 18禁国产床啪视频网站| 新久久久久国产一级毛片| 国产日韩欧美视频二区| 国产1区2区3区精品| 一二三四中文在线观看免费高清| 天堂8中文在线网| 欧美精品av麻豆av| av播播在线观看一区| 精品国产一区二区三区久久久樱花| 欧美老熟妇乱子伦牲交| 国产av一区二区精品久久| 另类亚洲欧美激情| 久久婷婷青草| 久久精品aⅴ一区二区三区四区 | xxx大片免费视频| 最近最新中文字幕免费大全7| av又黄又爽大尺度在线免费看| 日韩制服骚丝袜av| 国产精品.久久久| 少妇的逼水好多| 国产又爽黄色视频| 狠狠精品人妻久久久久久综合| 黄色视频在线播放观看不卡| av国产久精品久网站免费入址| 国产精品一区二区在线不卡| 国产精品国产三级国产av玫瑰| 亚洲激情五月婷婷啪啪| a级毛片在线看网站| 国产精品久久久久久久电影| 亚洲国产日韩一区二区| 中国美白少妇内射xxxbb| 日韩精品免费视频一区二区三区 | 99热6这里只有精品| 欧美人与性动交α欧美软件 | 精品视频人人做人人爽| 国产在线一区二区三区精| 久久久精品区二区三区| 天堂8中文在线网| 两个人看的免费小视频| 男女高潮啪啪啪动态图| 在线观看美女被高潮喷水网站| 狂野欧美激情性xxxx在线观看| 欧美日韩视频高清一区二区三区二| www.熟女人妻精品国产 | 黄色 视频免费看| 晚上一个人看的免费电影| 飞空精品影院首页| 黄色毛片三级朝国网站| 成人漫画全彩无遮挡| 久久久精品94久久精品| 桃花免费在线播放| 国产精品一国产av| freevideosex欧美| 满18在线观看网站| 欧美日韩综合久久久久久| 欧美激情国产日韩精品一区| 国产淫语在线视频| 又黄又爽又刺激的免费视频.| 男女午夜视频在线观看 | 秋霞伦理黄片| 熟妇人妻不卡中文字幕| 97精品久久久久久久久久精品| 黑人猛操日本美女一级片| 看十八女毛片水多多多| 国产精品久久久久久精品古装| 一级a做视频免费观看| 满18在线观看网站| 美女中出高潮动态图| 日韩免费高清中文字幕av| 人人妻人人澡人人看| 女人久久www免费人成看片| 九色亚洲精品在线播放| 亚洲成国产人片在线观看| 免费av不卡在线播放| 婷婷成人精品国产| av国产精品久久久久影院| 日本色播在线视频| 亚洲激情五月婷婷啪啪| 亚洲激情五月婷婷啪啪| 一区二区三区精品91| 国产精品一国产av| 欧美日韩成人在线一区二区| 日本av免费视频播放| 26uuu在线亚洲综合色| 欧美精品一区二区免费开放| 美女主播在线视频| 最黄视频免费看| 99香蕉大伊视频| 观看av在线不卡| 免费播放大片免费观看视频在线观看| 又大又黄又爽视频免费| 日韩不卡一区二区三区视频在线| 成人国产麻豆网| 国产一区亚洲一区在线观看| 国产精品免费大片| 两个人看的免费小视频| 视频中文字幕在线观看| 2021少妇久久久久久久久久久| 性色av一级| 夜夜爽夜夜爽视频| 日韩在线高清观看一区二区三区| 国产欧美亚洲国产| 亚洲av日韩在线播放| 久久99蜜桃精品久久| 黑人高潮一二区| 亚洲精品国产色婷婷电影| 精品国产国语对白av| a级毛片黄视频| 久久久精品区二区三区| 欧美成人午夜免费资源| 精品人妻一区二区三区麻豆| videos熟女内射| 精品人妻在线不人妻| 伦精品一区二区三区| 国产av精品麻豆| 国产永久视频网站| 亚洲色图 男人天堂 中文字幕 | 日本wwww免费看| 我的女老师完整版在线观看| 精品午夜福利在线看| www.熟女人妻精品国产 | 在线观看免费视频网站a站| 亚洲欧美成人精品一区二区| 亚洲中文av在线| 午夜福利在线观看免费完整高清在| 精品福利永久在线观看| 日产精品乱码卡一卡2卡三| 亚洲色图 男人天堂 中文字幕 | 美女xxoo啪啪120秒动态图| 大片免费播放器 马上看| 超色免费av| 国产精品一国产av| av片东京热男人的天堂| 男女高潮啪啪啪动态图| 人人妻人人添人人爽欧美一区卜| 毛片一级片免费看久久久久| 咕卡用的链子| 国产成人a∨麻豆精品| 男男h啪啪无遮挡| 国产精品99久久99久久久不卡 | 99精国产麻豆久久婷婷| 边亲边吃奶的免费视频| 精品一区二区免费观看| 精品一区在线观看国产| 免费观看性生交大片5| 美国免费a级毛片| 男人舔女人的私密视频| 国产不卡av网站在线观看| 大香蕉久久网| 欧美日韩成人在线一区二区| 国产精品人妻久久久影院| 国产免费一区二区三区四区乱码| 免费看av在线观看网站| 欧美变态另类bdsm刘玥| 99国产综合亚洲精品| 一级毛片黄色毛片免费观看视频| 蜜臀久久99精品久久宅男| 中国三级夫妇交换| 久久青草综合色| 99热6这里只有精品| 久久久久久久亚洲中文字幕| 伦精品一区二区三区| 中文字幕人妻丝袜制服| 最近最新中文字幕大全免费视频 | 亚洲伊人色综图| 久久久a久久爽久久v久久| 亚洲精品一区蜜桃| 天天影视国产精品| 日本免费在线观看一区| 国产熟女午夜一区二区三区| 欧美精品国产亚洲| 国产男女内射视频| 一区二区三区四区激情视频| 我要看黄色一级片免费的| 亚洲精品成人av观看孕妇| 久久精品国产a三级三级三级| 欧美精品高潮呻吟av久久| 亚洲性久久影院| 免费观看av网站的网址| 女人精品久久久久毛片| xxxhd国产人妻xxx| 久久久久精品性色| 久久人人97超碰香蕉20202| 国产伦理片在线播放av一区| 日本-黄色视频高清免费观看| 99香蕉大伊视频| 精品一区在线观看国产| 亚洲情色 制服丝袜| 日本爱情动作片www.在线观看| xxxhd国产人妻xxx| 久久国产亚洲av麻豆专区| 午夜91福利影院| 久久精品夜色国产| 日本欧美视频一区| 99精国产麻豆久久婷婷| 国产精品一国产av| 最近中文字幕高清免费大全6| 香蕉精品网在线| 日本黄色日本黄色录像| 熟女电影av网| 免费人妻精品一区二区三区视频| 久久这里只有精品19| 黑人欧美特级aaaaaa片| 亚洲婷婷狠狠爱综合网| 亚洲精品第二区| 十八禁高潮呻吟视频| a级片在线免费高清观看视频| 岛国毛片在线播放| 最新的欧美精品一区二区| 日韩av免费高清视频| 亚洲 欧美一区二区三区| 亚洲精品456在线播放app| 午夜福利,免费看| 亚洲天堂av无毛| 色视频在线一区二区三区| 熟女人妻精品中文字幕| 国产日韩欧美亚洲二区| 亚洲精品美女久久久久99蜜臀 | 看免费成人av毛片| 精品国产乱码久久久久久小说| 亚洲av电影在线观看一区二区三区| 国产老妇伦熟女老妇高清| 国产淫语在线视频| 亚洲一区二区三区欧美精品| 欧美少妇被猛烈插入视频| 亚洲欧美清纯卡通| 老熟女久久久| 国产欧美亚洲国产| 国产精品秋霞免费鲁丝片| 精品国产国语对白av| 国产无遮挡羞羞视频在线观看| 日韩电影二区| 日本-黄色视频高清免费观看| 爱豆传媒免费全集在线观看| 美国免费a级毛片| 在线观看一区二区三区激情| 美女大奶头黄色视频| 国国产精品蜜臀av免费| 18禁动态无遮挡网站| 久久女婷五月综合色啪小说| 国产极品粉嫩免费观看在线| 美女大奶头黄色视频| 成人18禁高潮啪啪吃奶动态图| 亚洲精品中文字幕在线视频| 成人亚洲欧美一区二区av| 18禁裸乳无遮挡动漫免费视频| www日本在线高清视频| 久久这里有精品视频免费| 成人毛片60女人毛片免费| 国产亚洲欧美精品永久| 亚洲国产精品成人久久小说| 午夜视频国产福利| 最新的欧美精品一区二区| 免费大片18禁| 亚洲图色成人| 少妇被粗大的猛进出69影院 | 精品熟女少妇av免费看| 国产精品秋霞免费鲁丝片| 最近中文字幕高清免费大全6| 国产精品久久久久久久久免| 欧美老熟妇乱子伦牲交| 1024视频免费在线观看| 国产精品免费大片| 一本久久精品| 亚洲精品国产av蜜桃| 免费少妇av软件| 交换朋友夫妻互换小说| 又粗又硬又长又爽又黄的视频| 色视频在线一区二区三区| 99久国产av精品国产电影| 国产精品不卡视频一区二区| 免费女性裸体啪啪无遮挡网站| 哪个播放器可以免费观看大片| 亚洲国产精品一区三区| 曰老女人黄片| 欧美激情 高清一区二区三区| 久久久久久久久久久久大奶| 爱豆传媒免费全集在线观看| 国产精品一区www在线观看| 精品国产一区二区三区四区第35| 久久av网站| 亚洲精品自拍成人| 免费人妻精品一区二区三区视频| 色哟哟·www| 欧美人与性动交α欧美软件 | 国产av精品麻豆| 22中文网久久字幕| 各种免费的搞黄视频| 国产亚洲最大av| 欧美国产精品va在线观看不卡| tube8黄色片| 全区人妻精品视频| 久久久久久伊人网av| 午夜老司机福利剧场| 熟妇人妻不卡中文字幕| 成人免费观看视频高清| 波多野结衣一区麻豆| 黄色怎么调成土黄色| 欧美亚洲日本最大视频资源| 99国产精品免费福利视频| 香蕉国产在线看| 丰满迷人的少妇在线观看| 国产黄色视频一区二区在线观看| 亚洲国产看品久久| 纵有疾风起免费观看全集完整版| 亚洲国产欧美在线一区| 久久久久网色| 国产一区有黄有色的免费视频| 老司机影院毛片| 成人国产av品久久久| 边亲边吃奶的免费视频| 夫妻午夜视频| 另类精品久久| 18禁国产床啪视频网站| 亚洲av中文av极速乱| 青春草亚洲视频在线观看| 亚洲精品第二区| 亚洲人成网站在线观看播放| 国产精品 国内视频| 亚洲情色 制服丝袜| 免费黄网站久久成人精品| 高清黄色对白视频在线免费看| av天堂久久9| 韩国精品一区二区三区 | 亚洲成人手机| 亚洲精品久久成人aⅴ小说| 交换朋友夫妻互换小说| 欧美少妇被猛烈插入视频| 欧美3d第一页| 人人妻人人澡人人爽人人夜夜| 一本—道久久a久久精品蜜桃钙片| 日韩免费高清中文字幕av| 国产一区有黄有色的免费视频| 满18在线观看网站| 在线观看美女被高潮喷水网站| 极品少妇高潮喷水抽搐| 久久99一区二区三区| 九色成人免费人妻av| 国产片特级美女逼逼视频| 欧美变态另类bdsm刘玥| 国产精品 国内视频| 久久午夜综合久久蜜桃| 丝袜在线中文字幕| 国产免费现黄频在线看| 日韩三级伦理在线观看| 18禁在线无遮挡免费观看视频| 日韩视频在线欧美| 亚洲色图 男人天堂 中文字幕 | 少妇熟女欧美另类| 国产不卡av网站在线观看| 日韩精品有码人妻一区| 大香蕉97超碰在线| 97超碰精品成人国产| 日本午夜av视频| 嫩草影院入口| 国产精品麻豆人妻色哟哟久久| 视频在线观看一区二区三区| 啦啦啦中文免费视频观看日本| 久久精品国产亚洲av天美| 高清黄色对白视频在线免费看| 在线观看美女被高潮喷水网站| 精品一区二区三区四区五区乱码 | 国产av一区二区精品久久| 母亲3免费完整高清在线观看 | 国产精品一区二区在线观看99| 精品午夜福利在线看| 天堂俺去俺来也www色官网| 亚洲成色77777| 精品久久国产蜜桃| 精品99又大又爽又粗少妇毛片| 99国产综合亚洲精品| 午夜久久久在线观看| 在线亚洲精品国产二区图片欧美| 美女脱内裤让男人舔精品视频| 色吧在线观看| a 毛片基地| 亚洲成av片中文字幕在线观看 | 免费少妇av软件| 超色免费av| 99久久综合免费| 80岁老熟妇乱子伦牲交| 亚洲国产精品专区欧美| 久热久热在线精品观看| 日韩av在线免费看完整版不卡| 一本久久精品| 宅男免费午夜| 九色成人免费人妻av| 18禁观看日本| 国产精品国产三级专区第一集| 国产精品一区www在线观看| 在线天堂最新版资源| 中国国产av一级| 国产一区有黄有色的免费视频| 波多野结衣一区麻豆| 熟女人妻精品中文字幕| 99精国产麻豆久久婷婷| 校园人妻丝袜中文字幕| 精品人妻在线不人妻| 成人国语在线视频| 三上悠亚av全集在线观看| 午夜日本视频在线| 免费日韩欧美在线观看| 日韩熟女老妇一区二区性免费视频| 久久午夜综合久久蜜桃| 欧美97在线视频| 搡女人真爽免费视频火全软件| 丰满饥渴人妻一区二区三| 成人毛片a级毛片在线播放| 丰满少妇做爰视频| 亚洲久久久国产精品| 97在线视频观看| 在线观看一区二区三区激情| 久久久久精品久久久久真实原创| 国产亚洲一区二区精品| 22中文网久久字幕| 99久国产av精品国产电影| 亚洲人成77777在线视频| 日本爱情动作片www.在线观看| 国产午夜精品一二区理论片| 爱豆传媒免费全集在线观看| 欧美成人午夜免费资源| 日本与韩国留学比较| 狂野欧美激情性bbbbbb| 亚洲av日韩在线播放| 人妻人人澡人人爽人人| 五月开心婷婷网| 色视频在线一区二区三区| 国产午夜精品一二区理论片| 大片电影免费在线观看免费| 涩涩av久久男人的天堂| 国产国语露脸激情在线看| 黄色 视频免费看| 大香蕉久久成人网| 亚洲五月色婷婷综合| 美女大奶头黄色视频| 我要看黄色一级片免费的| 国产成人精品在线电影| 一区二区三区四区激情视频| 亚洲性久久影院| 搡老乐熟女国产| av视频免费观看在线观看| 日韩中文字幕视频在线看片| 国产深夜福利视频在线观看| 欧美日韩精品成人综合77777| a级毛色黄片| 男的添女的下面高潮视频| 欧美国产精品va在线观看不卡| 插逼视频在线观看| 精品卡一卡二卡四卡免费| 日韩熟女老妇一区二区性免费视频| 啦啦啦中文免费视频观看日本| 欧美+日韩+精品| 久久ye,这里只有精品| 日韩一区二区视频免费看| 少妇的逼好多水| 成人综合一区亚洲| 五月玫瑰六月丁香| 91成人精品电影| 美女国产视频在线观看| 亚洲美女搞黄在线观看| 免费在线观看黄色视频的| 晚上一个人看的免费电影| 国产 一区精品| 国产一级毛片在线| 波多野结衣一区麻豆| 蜜桃国产av成人99| 久久99蜜桃精品久久| 欧美成人午夜精品| av在线观看视频网站免费| 免费黄色在线免费观看| 22中文网久久字幕| 水蜜桃什么品种好| 一级,二级,三级黄色视频| 日韩 亚洲 欧美在线| 精品人妻偷拍中文字幕| 国产成人aa在线观看| 22中文网久久字幕| 九色亚洲精品在线播放| 一级,二级,三级黄色视频| 最近的中文字幕免费完整| 91午夜精品亚洲一区二区三区| 亚洲国产色片| 母亲3免费完整高清在线观看 | 久久久国产一区二区| 精品国产一区二区三区四区第35| 国产成人精品婷婷| 亚洲精品久久成人aⅴ小说| 欧美另类一区| 赤兔流量卡办理| 久久久久人妻精品一区果冻| 国产伦理片在线播放av一区| 亚洲情色 制服丝袜| 亚洲,欧美,日韩| 亚洲av电影在线观看一区二区三区| 尾随美女入室| 免费av中文字幕在线| av线在线观看网站| 99热国产这里只有精品6| 乱人伦中国视频| 亚洲一级一片aⅴ在线观看| 亚洲第一av免费看| 有码 亚洲区| 宅男免费午夜| 久久久久久久精品精品| 女性生殖器流出的白浆| 最后的刺客免费高清国语| 天堂中文最新版在线下载| 亚洲国产色片| 亚洲av国产av综合av卡| 亚洲成人一二三区av| 最黄视频免费看| 亚洲伊人色综图| 五月玫瑰六月丁香| 成人午夜精彩视频在线观看| 中国国产av一级| 免费大片黄手机在线观看| 18禁在线无遮挡免费观看视频| 插逼视频在线观看| 国产精品久久久久久精品电影小说| 男女啪啪激烈高潮av片| 国产69精品久久久久777片| 少妇被粗大猛烈的视频| 高清在线视频一区二区三区| 精品久久蜜臀av无| 亚洲丝袜综合中文字幕| 亚洲精品国产色婷婷电影|