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

    基于分頻編碼的彈性波全波形反演

    2020-12-09 01:51:02邵祥奇何兵壽史才旺
    石油地球物理勘探 2020年6期
    關(guān)鍵詞:編碼方法炮點(diǎn)子波

    邵祥奇 何兵壽* 史才旺

    (①中國海洋大學(xué)海底科學(xué)與探測技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,山東青島 266100; ②青島海洋科學(xué)與技術(shù)國家實(shí)驗(yàn)室海洋礦產(chǎn)資源評價(jià)與探測技術(shù)功能實(shí)驗(yàn)室,山東青島 266100; ③南方科技大學(xué)地球與空間科學(xué)系,廣東深圳 518000)

    0 引言

    近年來,全波形反演(FWI)方法在理論和實(shí)際應(yīng)用方面均取得了較大進(jìn)展[1-5]。但是,F(xiàn)WI的龐大計(jì)算量和過低的計(jì)算效率是限制該方法進(jìn)入三維地震勘探領(lǐng)域的瓶頸之一。目前,業(yè)界提高反演效率的思路主要有兩種: 一是提高收斂速度; 二是提升正演(反演需要不斷進(jìn)行多次正演)計(jì)算速度。前者通過改進(jìn)反演算法減少迭代次數(shù),降低計(jì)算量[6-10]; 后者采用高性能的計(jì)算機(jī)硬件提高計(jì)算效率或采用降維思路降低計(jì)算量[11-26]。

    本文主要討論第二種思路。FWI大部分運(yùn)算量集中在正演計(jì)算過程中(FWI的每次迭代至少需要三次正演),因此業(yè)界普遍認(rèn)為高效的正演計(jì)算是解決FWI反演效率問題的有效途徑,采用降維思路減少運(yùn)算量是提高正演效率的有效方法。降維技術(shù)的本質(zhì)是采用某種壓縮方式減少參與反演的炮數(shù),主要方法包括震源編碼法[16-20]、炮采樣法[21-23]和主成分分析法[24-28]等。

    Romero等[16]在地震資料的逆時(shí)偏移處理中提出了相位編碼技術(shù),將多炮數(shù)據(jù)壓縮為幾個(gè)超級炮進(jìn)行計(jì)算,但編碼算法會在偏移結(jié)果中產(chǎn)生串?dāng)_噪聲。Krebs等[17]將相位編碼技術(shù)引入FWI領(lǐng)域,利用隨機(jī)的極性編碼函數(shù)壓制反演中的串?dāng)_噪聲,將所有炮轉(zhuǎn)換成一個(gè)超級炮集記錄并進(jìn)行反演,提升了運(yùn)算效率。此外,學(xué)者們還提出了平面波編碼[18]、頻域隨機(jī)相位編碼[19]、頻率組編碼[20]、頻域分割編碼[27]等方法,完善并豐富了編碼方案。

    炮采樣法也是一種高效的加速策略。Díaz等[21]提出了一種隨機(jī)抽取炮集進(jìn)行反演的方法,每次迭代時(shí)僅使用部分炮,在保證效果的同時(shí)能夠減少運(yùn)算。Ha等[22]提出了在Laplace域循環(huán)炮采樣,不同于隨機(jī)炮采樣,該方法循環(huán)使用炮集中的所有炮,使反演中每一炮的使用次數(shù)較均衡。Shi等[23]將炮采樣技術(shù)應(yīng)用到頻率多尺度FWI中,在每個(gè)頻率段均隨機(jī)抽取炮集參與反演,該方法不僅能加速運(yùn)算,而且不會引入額外的串?dāng)_噪聲。

    主成分分析法也是一種經(jīng)典的降維方法,Liu等[24]將主成分分析法應(yīng)用于FWI,減少了炮集的數(shù)量,提高了計(jì)算效率。段超然等[25]在高頻部分使用主成分分析法,使FWI在缺少低頻成分時(shí)也能高效收斂。

    震源編碼技術(shù)實(shí)現(xiàn)簡單、效果明顯,是目前應(yīng)用最廣泛的加速方法,但同時(shí)也存在一些缺陷,主要表現(xiàn)在傳統(tǒng)震源編碼技術(shù)嚴(yán)重依賴于觀測系統(tǒng)。震源編碼技術(shù)應(yīng)用的前提是采用檢波點(diǎn)固定的方式觀測,因而在模擬同時(shí)激發(fā)震源時(shí),無法對單炮獲取特定炮檢距的地震記錄。對于滾動排列觀測系統(tǒng),震源編碼會造成模擬與觀測資料的不匹配。鑒于此,Choi 等[29]提出了基于歸一化互相關(guān)目標(biāo)函數(shù)的波形反演,提升了震源編碼對海上拖纜資料的反演精度。夏冬明等[30]在此基礎(chǔ)上發(fā)展了基于歸一化局部互相關(guān)目標(biāo)函數(shù)的反演方法,進(jìn)一步改善了對拖纜數(shù)據(jù)的反演效果。Huang等[31]提出了基于分頻編碼的最小二乘偏移方法,能夠真正實(shí)現(xiàn)拖纜資料中多炮數(shù)據(jù)的完全分離。Zhang等[32]將分頻編碼應(yīng)用于聲波FWI,并且討論了其在計(jì)算效率、內(nèi)存消耗、抗噪性方面的優(yōu)點(diǎn)。

    在前人研究的基礎(chǔ)上,本文主要解決目前編碼方法只適用于特殊觀測系統(tǒng)的缺陷,提出了一種基于震源分頻編碼的混合域彈性波FWI方法。首先為同時(shí)進(jìn)行正向延拓的每一炮賦以不同的頻率,并利用該頻率對應(yīng)的諧波代替原時(shí)間域子波進(jìn)行正演,通過相位靈敏度檢測(Phase Sensitivity Detection,PSD)提取單頻波場,從而實(shí)現(xiàn)各炮頻率域波場的完全分離;然后實(shí)現(xiàn)多炮編碼反演。這樣既能壓制炮間串?dāng)_噪聲,又能對每一炮波場和記錄實(shí)現(xiàn)靈活處理,可以適應(yīng)排列滾動的觀測系統(tǒng)。

    1 震源編碼技術(shù)

    混合域彈性波FWI的運(yùn)算量和炮集數(shù)成正比,運(yùn)用震源編碼技術(shù)將多炮數(shù)據(jù)進(jìn)行組合、正演,可有效減少計(jì)算量。震源編碼技術(shù)借助于編碼函數(shù),將炮集疊加、組合得到超級炮記錄,利用常規(guī)方法對該超級炮記錄進(jìn)行反演。

    在頻率域,超級炮記錄可表示為

    (1)

    式中: 上標(biāo)“~”表示超級炮;dn是第n炮得到的觀測地震記錄;γn表示第n炮的編碼函數(shù)序列;NS為總炮數(shù)。

    相對應(yīng)的超級炮震源可表示為

    (2)

    (3)

    (4)

    式中:Pn是第n炮的正傳波場;Bn是第n炮的殘差波場。

    推導(dǎo)得出超級炮對應(yīng)的梯度

    (5)

    式中:E(m)為目標(biāo)函數(shù),m為待反演參數(shù),本文指速度; R(·)表示取實(shí)部;上角標(biāo)“*”表示取共軛;A是波動方程的阻抗矩陣,滿足AP=s,s為震源子波。

    2 分頻編碼技術(shù)

    2.1 分頻編碼策略

    本文把分頻編碼技術(shù)引入彈性波波形反演中。其核心思想是:在正演過程中,同時(shí)對多個(gè)炮記錄進(jìn)行正向延拓,給每一炮賦以不同頻率的諧波震源,使同時(shí)正向延拓的各炮波場在頻譜上互不重合。同時(shí),在炮點(diǎn)波場正傳和殘差反傳過程中,采用相位靈敏度檢測技術(shù)提取每炮的單頻波場,實(shí)現(xiàn)各炮波場的完全分離,避免炮間串?dāng)_。

    2.2 基于相位靈敏度檢測的分頻多炮正演

    分頻編碼技術(shù)的核心在于給每炮賦以不同頻帶的震源函數(shù),避免多炮同時(shí)延拓中的炮間串?dāng)_。震源的頻帶一般由震源函數(shù)的中心頻率控制,不同炮可以選擇不同主頻的震源函數(shù)。但由于不同主頻的震源函數(shù)往往具有一定的頻帶寬度,當(dāng)炮數(shù)較多時(shí),難以完全分離多炮頻譜,利用諧波作為震源可以避免這一缺陷。

    利用諧波震源進(jìn)行正演,首先需要提取正確的波場信息。本文引入PSD正演方法實(shí)現(xiàn)頻率域波場的準(zhǔn)確提取。Nihei等[26]提出的基于PSD的正演方法用諧波作為震源,在正演過程中可以很方便地提取出不同頻率的單頻波場,實(shí)現(xiàn)多震源的同步模擬且互不干擾。

    PSD方法是利用一個(gè)參考信號和參考信號的90°相移信號計(jì)算原諧波信號的頻譜。對于含有兩個(gè)頻率的信號,原信號εsig用振幅Esig和相位θsig表示為

    εsig=Esig1cos(ω1t+θsig1)+Esig2cos(ω2t+θsig2)

    (6)

    式中下標(biāo)“1”和“2”分別對應(yīng)這兩個(gè)頻率。

    選取兩個(gè)不同的頻率ω1和ω2, 令它們滿足以下關(guān)系

    (7)

    式中k為正整數(shù)。

    以ω1為例,設(shè)計(jì)如下參考信號和相移信號

    (8)

    式中 ref1表示參考信號。

    ω1對應(yīng)的單頻信號的振幅和相位[26]可以表示為

    (9)

    其中

    (10)

    Nihei等[26]證明PSD算法與傅里葉變換的功能是相同的,但是前者更適用于對諧波信號進(jìn)行處理。對只包含幾個(gè)等間隔頻率的信號,PSD算法僅需要選擇很小長度的信號進(jìn)行運(yùn)算就能夠得到其頻率響應(yīng),效率更高。

    在時(shí)間域正演過程中,通常使用Ricker子波作為震源,但PSD算法需要以諧波震源作為輸入。為了得到與常規(guī)正演方法等價(jià)的結(jié)果,需要建立Ri-cker子波和諧波震源之間的數(shù)學(xué)聯(lián)系。實(shí)際上,諧波震源能夠視為Ricker子波的某一個(gè)或幾個(gè)頻率成分,所以僅需要將Ricker子波做離散傅里葉變換后,提取出某幾個(gè)頻率的響應(yīng),再做傅里葉反變換,就能得到諧波震源。處理實(shí)際資料時(shí),對真實(shí)的子波做類似操作即可。當(dāng)然,這樣做的前提是能夠獲取較為準(zhǔn)確的震源子波,否則使用不正確的子波將會導(dǎo)致錯(cuò)誤的反演結(jié)果。子波的提取是波形反演需要解決的一個(gè)重要問題,但是本文的研究目標(biāo)是加快反演速度,暫不考慮如何提取子波的問題。另外,本文方法可以與當(dāng)前的不依賴于子波的反演方法相結(jié)合[33-35],即使利用錯(cuò)誤的震源子波進(jìn)行反演也能得到正確的反演結(jié)果。

    本文基于PSD的彈性波正演方法流程如圖1所示。

    利用PSD算法實(shí)現(xiàn)分頻多炮正演實(shí)驗(yàn)可以驗(yàn)證PSD算法的有效性。以簡單的均勻介質(zhì)為例,網(wǎng)格數(shù)設(shè)為150×200,網(wǎng)格間距為10m,縱、橫波速度分別為2500、1500m/s。設(shè)有4個(gè)炮點(diǎn),炮點(diǎn)深度為20m,炮點(diǎn)等間隔分布,橫向坐標(biāo)分別為400、800、1200和1600m,原始子波是主頻8Hz的Ricker子波,時(shí)間長度設(shè)為4s,四個(gè)炮點(diǎn)分別同時(shí)正演2、3、4、5Hz的單頻波場。

    圖1 基于PSD的彈性波分頻編碼正演流程

    利用PSD方法進(jìn)行正演所得到的各炮點(diǎn)頻率響應(yīng)如圖2所示,其中幾乎見不到多炮的串?dāng)_。

    為精確驗(yàn)證PSD算法能夠有效避免炮間串?dāng)_,做了如下對照實(shí)驗(yàn):每次只設(shè)置一個(gè)震源一個(gè)頻率,其他參數(shù)與圖2保持一致,得到不含炮間串?dāng)_的波場。將兩次實(shí)驗(yàn)結(jié)果中橫坐標(biāo)1km處不同深度的波場值繪制成圖3。對比可知,兩次實(shí)驗(yàn)得到的波場幾乎一致,這也說明PSD算法具有足夠的精度,可以在分頻率的多炮正演中得到準(zhǔn)確的頻率域波場值。

    2.3 分頻編碼的彈性波混合域FWI

    從理論上講,分頻編碼技術(shù)不會引入額外噪聲[31]。本文應(yīng)用分頻編碼方法以實(shí)現(xiàn)混合域的彈性波反演。與Huang等[31]實(shí)現(xiàn)的最小二乘偏移方法不同,混合域FWI一般采用多尺度反演策略,在每一個(gè)尺度下,可供使用的頻段是受限的,所以頻率的選擇在分頻編碼混合域反演中非常重要。

    圖2 分頻編碼單次正演的各炮頻率響應(yīng)(實(shí)部)

    圖3 分頻編碼與單諧波震源波場值對比

    通常而言,隨多尺度反演的遞進(jìn),可供選擇的頻段逐漸拓寬。如果對每一尺度反演都設(shè)置相同的頻率數(shù)目,結(jié)合上述兩個(gè)限制條件,就不難得出如下結(jié)論:小尺度反演時(shí),頻率間隔大,Δω較大而PSD計(jì)算時(shí)長??;大尺度反演時(shí),頻率選擇較密集,Δω較小而PSD計(jì)算時(shí)長大。反之,若Δω固定不變,則在大尺度反演時(shí),可供選擇的頻率更少,進(jìn)一步造成能同時(shí)正演的炮數(shù)變少,正演次數(shù)就會明顯增加。本文選擇用較多的頻率個(gè)數(shù)完成多炮同時(shí)正演。具體到反演過程中,首先確定反演使用的頻率個(gè)數(shù)M,把每M炮數(shù)據(jù)組合成為一個(gè)超級炮(組合方法可以有多種,本文選擇把空間上鄰近的M炮數(shù)據(jù)進(jìn)行組合);然后每次迭代前把M個(gè)頻率的數(shù)據(jù)隨機(jī)分配給組合中的M炮,讓迭代過程中各炮使用的頻率不斷變化,這樣可以使原始數(shù)據(jù)得到充分利用。

    3 模型算例

    3.1 分頻編碼對FWI的加速效果

    為證明本文方法的加速效果,用部分Overthrust模型的正演記錄進(jìn)行測試。圖4為部分Overthrust縱、橫波速度模型,該模型網(wǎng)格數(shù)為500×160,網(wǎng)格間隔為25m×25m。正演記錄的炮點(diǎn)均勻分布于地表,第一炮位于0m處,炮間距為125m,共得到100炮合成地震記錄,全排列接收,道間距為25m。反演所用的初始模型(圖5)是將圖4進(jìn)行高斯平滑得到的。

    圖4 部分Overthrust速度模型

    圖5 反演使用的初始速度模型

    采用多尺度反演策略,每個(gè)頻率組包含10個(gè)頻率(表1),反演分5個(gè)頻段,每個(gè)頻段迭代15次。分別采用常規(guī)全部炮反演、傳統(tǒng)震源編碼(極性編碼)反演和分頻編碼反演,并將三者反演結(jié)果(圖6)進(jìn)行對比。

    表1 分頻編碼反演的頻率組信息

    圖7為三種方法的誤差對比, 反演誤差定義為[19]

    (11)

    式中:mFWI和mTRUE分別為反演結(jié)果和真實(shí)模型的彈性參數(shù); ‖·‖2表示L2范數(shù)。

    從圖6、圖7可以看出,分頻編碼方法和極性編碼方法的反演精度與常規(guī)方法相近。但就計(jì)算量而言,由于每次迭代需要4次正演(梯度計(jì)算需要兩次,步長估計(jì)需要兩次),因此常規(guī)全部炮反演需要計(jì)算30000(100×5×15×4)次正演,而本文方法與極性編碼一樣,將10炮數(shù)據(jù)組合為1炮進(jìn)行反演,因此正演次數(shù)僅為常規(guī)方法的十分之一。

    當(dāng)然,計(jì)算效率的提升與組合的炮數(shù)有關(guān)。通常反演中能夠進(jìn)行組合的炮數(shù)受頻帶寬度的影響,不能把任意多的炮集進(jìn)行組合。此算例中,每個(gè)超級炮都由固定的10炮數(shù)據(jù)組成。而在一般的多尺度反演中,大尺度反演時(shí)給定的頻帶寬度往往較小,能夠選擇的頻率點(diǎn)(亦即能夠進(jìn)行組合的炮數(shù))也比較少;隨著尺度變小,頻帶可以更寬,從而可以用更多的炮進(jìn)行組合。但是由于混合域算法中需要存儲每個(gè)頻率的波場,更多的頻率點(diǎn)意味著存儲量的增加,同時(shí)也意味著需更多的計(jì)算完成PSD的頻率分離,因此最終選擇多少個(gè)頻率需視情況而定。

    圖6 全排列接收時(shí)不同方法縱波速度(左)、橫波速度(右)反演結(jié)果

    圖7 誤差分析對比

    3.2 滾動排列觀測系統(tǒng)反演測試

    應(yīng)用傳統(tǒng)震源編碼方法反演時(shí),每炮的接收排列必須相同且覆蓋整個(gè)工區(qū),但是排列滾動的數(shù)據(jù)只在炮點(diǎn)附近一定范圍內(nèi)接收,遠(yuǎn)炮點(diǎn)并沒有數(shù)據(jù),所以面對滾動排列觀測系統(tǒng)傳統(tǒng)編碼方法反演效果很差。

    仍以圖4所示的部分Overthrust模型為例,采用中間放炮的方式進(jìn)行觀測,炮點(diǎn)均勻分布于地表,第一炮位于0m處,炮間距為125m,每炮300道接收,道間距為25m,在模型兩邊排列固定不動,炮點(diǎn)滾動,共得到100炮兩分量記錄。首先采用傳統(tǒng)震源編碼反演方法(極性編碼)對合成記錄進(jìn)行反演,反演結(jié)果如圖8a所示。顯然,常規(guī)震源編碼方法反演得到的結(jié)果與真實(shí)模型相差很大。由圖8a可以看出,橫波的反演結(jié)果優(yōu)于縱波,這與本文的實(shí)驗(yàn)條件有關(guān)。傳統(tǒng)震源編碼方法應(yīng)用于排列滾動的地震數(shù)據(jù)時(shí),誤差主要源于排列范圍之外(遠(yuǎn)炮檢距)的數(shù)據(jù)缺失。具體到本文的實(shí)驗(yàn),采用爆炸震源激發(fā),橫波分量均來自于波的類型轉(zhuǎn)換,因此占比較少。遠(yuǎn)炮檢距存在很強(qiáng)的縱波分量,尤其是直達(dá)縱波,這一部分?jǐn)?shù)據(jù)的缺失會對縱波反演產(chǎn)生極大的不利影響。結(jié)合圖6給出的反演結(jié)果,可以看出,傳統(tǒng)編碼方法在全排列接收時(shí)可以給出正確的反演結(jié)果,但是在滾動排列情況下會受到嚴(yán)重的干擾。

    利用本文的分頻編碼方法對排列滾動的數(shù)據(jù)進(jìn)行反演,結(jié)果如圖8b所示。由圖可見本文方法的反演結(jié)果結(jié)構(gòu)清晰,與真實(shí)模型基本一致,遠(yuǎn)遠(yuǎn)優(yōu)于圖8a的反演結(jié)果。計(jì)算效率方面,本文基于GPU集群環(huán)境利用MPI+CUDA實(shí)現(xiàn)多炮并行反演,采用16張Tesla K80運(yùn)算。反演共進(jìn)行75次迭代,分頻編碼方法完成反演用時(shí)70min,極性編碼方法用時(shí)66min,分頻編碼方法耗時(shí)略高。分頻編碼方法首先需要對記錄做分頻處理;其次,本文在實(shí)施相位靈敏度檢測時(shí)需要適當(dāng)延長正演時(shí)間,因此產(chǎn)生了一些額外計(jì)算量。但是這些額外的計(jì)算量相比于波場正演和反傳而言影響較小,所以整體上兩種編碼方法計(jì)算效率相近。

    圖8 流動排列接收時(shí)不同方法縱波(左)、橫波(右)速度反演結(jié)果

    4 結(jié)束語

    分頻編碼方法在不明顯降低反演精度的基礎(chǔ)上,能夠?qū)⑦\(yùn)算效率提升數(shù)倍,有效減少正演次數(shù),減小全波形反演的計(jì)算量。本文的分頻編碼方法首先利用諧波震源實(shí)現(xiàn)分頻的多炮同步正演,再利用基于相位靈敏度檢測的波場分離方法,實(shí)現(xiàn)不同炮波場的完全分離,從而能夠在梯度計(jì)算中避免引入多炮串?dāng)_噪聲。另外,本文使用的這種分離方法僅依靠頻率信息,與震源、接收點(diǎn)的排列布置無關(guān),能夠適用于滾動排列觀測系統(tǒng),彌補(bǔ)了傳統(tǒng)震源編碼方法在觀測系統(tǒng)適應(yīng)性上存在的不足。

    猜你喜歡
    編碼方法炮點(diǎn)子波
    一類非線性動力系統(tǒng)的孤立子波解
    可變摩擦力觸感移動終端的漢語盲文編碼設(shè)計(jì)
    基于最小炮檢距道快速檢測炮點(diǎn)偏移方法
    對二維地震覆蓋次數(shù)和炮點(diǎn)距的思考
    山西冶金(2018年6期)2018-03-04 12:53:26
    無樁號施工中炮點(diǎn)COG現(xiàn)場快速偏移技術(shù)
    一二八團(tuán)開展“夏送清涼”慰問
    毫米波大規(guī)模MIMO系統(tǒng)中低復(fù)雜度混合預(yù)編碼方法
    地震反演子波選擇策略研究
    基于倒雙譜的地震子波估計(jì)方法
    一種新的星載InSAR直接地理編碼方法
    日本av手机在线免费观看| 国产男女超爽视频在线观看| 国产精品一区二区免费欧美| 亚洲熟女精品中文字幕| 蜜桃在线观看..| 曰老女人黄片| 亚洲国产欧美在线一区| 国产成人精品久久二区二区91| 欧美 日韩 精品 国产| 久久人妻av系列| 最新的欧美精品一区二区| 欧美变态另类bdsm刘玥| 亚洲欧美日韩另类电影网站| 亚洲精品乱久久久久久| 欧美精品人与动牲交sv欧美| 久久久精品免费免费高清| 黄色成人免费大全| av片东京热男人的天堂| 欧美精品一区二区大全| 亚洲熟女毛片儿| 两个人看的免费小视频| 麻豆国产av国片精品| 国内毛片毛片毛片毛片毛片| 午夜福利欧美成人| 黄色 视频免费看| 国产成人av激情在线播放| 18禁黄网站禁片午夜丰满| 女人久久www免费人成看片| 天堂俺去俺来也www色官网| 午夜免费鲁丝| 女人精品久久久久毛片| 国产免费福利视频在线观看| 欧美在线黄色| 亚洲精品成人av观看孕妇| 亚洲 欧美一区二区三区| 精品福利永久在线观看| 国产一卡二卡三卡精品| 久久久久久人人人人人| 窝窝影院91人妻| av片东京热男人的天堂| 肉色欧美久久久久久久蜜桃| 丝袜喷水一区| 老司机深夜福利视频在线观看| 久久久久视频综合| 国产成人精品久久二区二区91| 日本精品一区二区三区蜜桃| 999久久久国产精品视频| 午夜91福利影院| 久久久久国产一级毛片高清牌| 日韩人妻精品一区2区三区| 捣出白浆h1v1| 人人妻,人人澡人人爽秒播| 欧美国产精品va在线观看不卡| 国产精品 欧美亚洲| 久9热在线精品视频| 香蕉丝袜av| 午夜福利乱码中文字幕| 日韩一区二区三区影片| 黄色怎么调成土黄色| 国产成人免费观看mmmm| 纵有疾风起免费观看全集完整版| 欧美黑人精品巨大| 欧美 日韩 精品 国产| 午夜日韩欧美国产| 黑人猛操日本美女一级片| 2018国产大陆天天弄谢| 日本a在线网址| 日韩视频在线欧美| 男人舔女人的私密视频| 国产视频一区二区在线看| 丝袜美足系列| 两个人看的免费小视频| 少妇的丰满在线观看| tocl精华| 国产有黄有色有爽视频| 精品亚洲乱码少妇综合久久| 一本综合久久免费| 国产精品 欧美亚洲| 99热网站在线观看| 99riav亚洲国产免费| 亚洲九九香蕉| 免费在线观看完整版高清| av欧美777| 国产不卡av网站在线观看| 国产淫语在线视频| 中文亚洲av片在线观看爽 | 18禁美女被吸乳视频| 久久毛片免费看一区二区三区| 成人av一区二区三区在线看| 国产精品98久久久久久宅男小说| 两个人看的免费小视频| 18禁黄网站禁片午夜丰满| 91av网站免费观看| 欧美黑人精品巨大| 精品亚洲成a人片在线观看| 啦啦啦视频在线资源免费观看| 欧美国产精品一级二级三级| 国产精品免费大片| 青草久久国产| 性色av乱码一区二区三区2| 高清黄色对白视频在线免费看| 免费观看人在逋| 国产亚洲欧美在线一区二区| 1024香蕉在线观看| 80岁老熟妇乱子伦牲交| 久久精品国产亚洲av香蕉五月 | 午夜福利在线免费观看网站| 成人av一区二区三区在线看| 国产成人一区二区三区免费视频网站| 麻豆av在线久日| 大型av网站在线播放| 日韩熟女老妇一区二区性免费视频| 亚洲视频免费观看视频| 我的亚洲天堂| 人人澡人人妻人| 国产在线精品亚洲第一网站| 黑人猛操日本美女一级片| 亚洲精品国产色婷婷电影| 最新美女视频免费是黄的| 在线av久久热| 欧美人与性动交α欧美精品济南到| 女性生殖器流出的白浆| 日韩欧美一区二区三区在线观看 | 午夜视频精品福利| 欧美 亚洲 国产 日韩一| 高清欧美精品videossex| 91九色精品人成在线观看| 亚洲午夜理论影院| 亚洲 欧美一区二区三区| 国产激情久久老熟女| 亚洲国产欧美在线一区| 美女扒开内裤让男人捅视频| 成人18禁在线播放| 久久天躁狠狠躁夜夜2o2o| 中文字幕人妻丝袜制服| 欧美黑人精品巨大| 国产视频一区二区在线看| 久久99热这里只频精品6学生| 熟女少妇亚洲综合色aaa.| 精品福利永久在线观看| 久久久久久久久久久久大奶| 露出奶头的视频| avwww免费| 少妇被粗大的猛进出69影院| 国产免费av片在线观看野外av| 亚洲精品久久成人aⅴ小说| 国产淫语在线视频| av视频免费观看在线观看| 日韩免费av在线播放| 老熟妇乱子伦视频在线观看| 亚洲一区中文字幕在线| 99精国产麻豆久久婷婷| 亚洲一码二码三码区别大吗| 人成视频在线观看免费观看| 成人免费观看视频高清| 热99久久久久精品小说推荐| 1024香蕉在线观看| 啦啦啦免费观看视频1| 亚洲精品中文字幕在线视频| 久久国产精品影院| 美女国产高潮福利片在线看| 一本久久精品| 国产亚洲精品一区二区www | 亚洲av日韩在线播放| 国产伦人伦偷精品视频| 免费看十八禁软件| 欧美成狂野欧美在线观看| 欧美日韩一级在线毛片| 欧美另类亚洲清纯唯美| 9色porny在线观看| 免费一级毛片在线播放高清视频 | 午夜免费成人在线视频| 国产精品影院久久| 国产在线观看jvid| 午夜福利一区二区在线看| 91麻豆精品激情在线观看国产 | 少妇精品久久久久久久| 欧美亚洲日本最大视频资源| 亚洲欧美日韩另类电影网站| 黄片播放在线免费| 2018国产大陆天天弄谢| 菩萨蛮人人尽说江南好唐韦庄| 成人国语在线视频| 十分钟在线观看高清视频www| 欧美成人免费av一区二区三区 | 黄色丝袜av网址大全| 久久午夜综合久久蜜桃| 亚洲av成人一区二区三| 国产单亲对白刺激| 久久午夜亚洲精品久久| 午夜福利在线观看吧| 啦啦啦 在线观看视频| 国产在线观看jvid| 妹子高潮喷水视频| 又黄又粗又硬又大视频| 久久人人爽av亚洲精品天堂| e午夜精品久久久久久久| 久久免费观看电影| 国产精品偷伦视频观看了| 十分钟在线观看高清视频www| 黑人欧美特级aaaaaa片| 无限看片的www在线观看| 999久久久精品免费观看国产| 亚洲少妇的诱惑av| 日韩欧美一区视频在线观看| 久久久水蜜桃国产精品网| 亚洲av日韩精品久久久久久密| 欧美在线一区亚洲| 麻豆国产av国片精品| 亚洲人成伊人成综合网2020| 成人国语在线视频| 成人18禁在线播放| 淫妇啪啪啪对白视频| 国产免费视频播放在线视频| 色播在线永久视频| 超碰成人久久| 热99国产精品久久久久久7| 不卡av一区二区三区| 色老头精品视频在线观看| 亚洲专区字幕在线| 欧美日韩国产mv在线观看视频| 一区二区三区乱码不卡18| 咕卡用的链子| 亚洲一区中文字幕在线| 国产片内射在线| 亚洲精品国产色婷婷电影| av一本久久久久| 久久久久久久久免费视频了| 久久国产亚洲av麻豆专区| 天堂8中文在线网| 丁香六月欧美| 精品国内亚洲2022精品成人 | 成人永久免费在线观看视频 | 日韩制服丝袜自拍偷拍| 十八禁高潮呻吟视频| 新久久久久国产一级毛片| 国产亚洲精品第一综合不卡| 久久亚洲精品不卡| 久热这里只有精品99| 18禁黄网站禁片午夜丰满| 99久久精品国产亚洲精品| 法律面前人人平等表现在哪些方面| 久久午夜综合久久蜜桃| 精品卡一卡二卡四卡免费| 亚洲国产欧美在线一区| www.精华液| 99久久精品国产亚洲精品| 美女主播在线视频| 亚洲国产av新网站| 中文欧美无线码| av超薄肉色丝袜交足视频| 日本av免费视频播放| 精品福利永久在线观看| 69精品国产乱码久久久| 777米奇影视久久| 人人澡人人妻人| 两性午夜刺激爽爽歪歪视频在线观看 | 91字幕亚洲| 亚洲美女黄片视频| 久久毛片免费看一区二区三区| 丝瓜视频免费看黄片| 国产老妇伦熟女老妇高清| 色视频在线一区二区三区| 人人妻人人添人人爽欧美一区卜| 国产99久久九九免费精品| 丁香六月欧美| 成人18禁高潮啪啪吃奶动态图| 免费av中文字幕在线| 国产精品一区二区免费欧美| 精品国内亚洲2022精品成人 | 在线观看免费视频网站a站| 老汉色av国产亚洲站长工具| 99精国产麻豆久久婷婷| 在线天堂中文资源库| 精品久久久久久久毛片微露脸| 1024视频免费在线观看| 桃红色精品国产亚洲av| 日韩有码中文字幕| 大码成人一级视频| 淫妇啪啪啪对白视频| 一本久久精品| 女同久久另类99精品国产91| 欧美日韩黄片免| 国产真人三级小视频在线观看| 91成人精品电影| 亚洲精品一二三| 国产男女超爽视频在线观看| 桃红色精品国产亚洲av| 中文字幕色久视频| 一区在线观看完整版| e午夜精品久久久久久久| 欧美精品av麻豆av| 法律面前人人平等表现在哪些方面| 亚洲色图综合在线观看| a级毛片黄视频| 亚洲av成人不卡在线观看播放网| 亚洲国产av新网站| 精品熟女少妇八av免费久了| av一本久久久久| 美国免费a级毛片| 欧美乱码精品一区二区三区| 午夜91福利影院| 国产日韩欧美视频二区| 亚洲精品在线美女| 国产日韩欧美亚洲二区| 天天操日日干夜夜撸| 男人操女人黄网站| 久久这里只有精品19| 成人18禁在线播放| 久久青草综合色| 亚洲va日本ⅴa欧美va伊人久久| 精品国产亚洲在线| 黄色视频不卡| 无人区码免费观看不卡 | 免费一级毛片在线播放高清视频 | 久久 成人 亚洲| 中文亚洲av片在线观看爽 | 丝袜在线中文字幕| 欧美黑人欧美精品刺激| av超薄肉色丝袜交足视频| 欧美久久黑人一区二区| 亚洲熟妇熟女久久| 欧美日韩亚洲综合一区二区三区_| av超薄肉色丝袜交足视频| 亚洲情色 制服丝袜| 超碰97精品在线观看| 成年动漫av网址| 嫁个100分男人电影在线观看| 国产精品秋霞免费鲁丝片| 欧美久久黑人一区二区| 欧美成人午夜精品| 大片电影免费在线观看免费| 久久人妻av系列| 国产精品99久久99久久久不卡| 搡老岳熟女国产| 丝袜美腿诱惑在线| 国产欧美日韩综合在线一区二区| 久久午夜综合久久蜜桃| 一级毛片精品| 99国产精品一区二区三区| 999精品在线视频| 99久久人妻综合| 欧美午夜高清在线| 黄色视频不卡| av福利片在线| 一本色道久久久久久精品综合| 精品第一国产精品| 9色porny在线观看| 一级黄色大片毛片| 新久久久久国产一级毛片| 女同久久另类99精品国产91| 九色亚洲精品在线播放| 97人妻天天添夜夜摸| 国产精品影院久久| 国产在线视频一区二区| 久久精品亚洲av国产电影网| av网站在线播放免费| av电影中文网址| 丁香六月天网| 一级片免费观看大全| 熟女少妇亚洲综合色aaa.| 757午夜福利合集在线观看| 国产又爽黄色视频| 日韩视频一区二区在线观看| 午夜久久久在线观看| 久久婷婷成人综合色麻豆| 99在线人妻在线中文字幕 | 亚洲天堂av无毛| 真人做人爱边吃奶动态| 妹子高潮喷水视频| 老汉色av国产亚洲站长工具| 波多野结衣一区麻豆| 看免费av毛片| 18禁观看日本| 欧美激情极品国产一区二区三区| 别揉我奶头~嗯~啊~动态视频| 精品第一国产精品| 美女视频免费永久观看网站| www.999成人在线观看| 视频在线观看一区二区三区| 正在播放国产对白刺激| 手机成人av网站| 久久精品熟女亚洲av麻豆精品| 欧美精品人与动牲交sv欧美| 2018国产大陆天天弄谢| 日韩免费高清中文字幕av| 色在线成人网| 国内毛片毛片毛片毛片毛片| 亚洲成a人片在线一区二区| 大片免费播放器 马上看| 丝瓜视频免费看黄片| 国产欧美亚洲国产| 亚洲国产av新网站| 精品国产国语对白av| 波多野结衣av一区二区av| 香蕉丝袜av| 最新在线观看一区二区三区| 伦理电影免费视频| 亚洲中文字幕日韩| 老汉色av国产亚洲站长工具| 亚洲一卡2卡3卡4卡5卡精品中文| 人人妻人人添人人爽欧美一区卜| 久久久久久久国产电影| 欧美激情高清一区二区三区| 99久久精品国产亚洲精品| 在线观看舔阴道视频| 人人妻,人人澡人人爽秒播| 90打野战视频偷拍视频| 亚洲免费av在线视频| 亚洲一区中文字幕在线| 一区二区日韩欧美中文字幕| 亚洲一区二区三区欧美精品| 9色porny在线观看| 18禁观看日本| 亚洲男人天堂网一区| 国产精品亚洲av一区麻豆| 久久精品亚洲熟妇少妇任你| 国产真人三级小视频在线观看| 一夜夜www| 国产aⅴ精品一区二区三区波| 精品国产超薄肉色丝袜足j| 欧美日本中文国产一区发布| 99精品久久久久人妻精品| 国产1区2区3区精品| 在线观看人妻少妇| 69av精品久久久久久 | 一区二区日韩欧美中文字幕| a级毛片在线看网站| 欧美精品av麻豆av| 国产片内射在线| 飞空精品影院首页| 大陆偷拍与自拍| 91麻豆av在线| 在线观看免费视频日本深夜| 国产午夜精品久久久久久| 三级毛片av免费| 超碰成人久久| 国产一卡二卡三卡精品| 黄色片一级片一级黄色片| netflix在线观看网站| 中文字幕制服av| 国产精品av久久久久免费| 日本vs欧美在线观看视频| 9热在线视频观看99| 十八禁网站免费在线| 欧美日韩视频精品一区| 亚洲专区字幕在线| 亚洲第一青青草原| 久久天躁狠狠躁夜夜2o2o| 精品国产乱子伦一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 一夜夜www| 人人澡人人妻人| 日韩免费av在线播放| 久久99一区二区三区| 黑人猛操日本美女一级片| 麻豆乱淫一区二区| 中文字幕另类日韩欧美亚洲嫩草| av网站在线播放免费| 久久久精品区二区三区| 久久人妻福利社区极品人妻图片| 免费av中文字幕在线| 一区二区三区激情视频| 亚洲人成77777在线视频| 黑人猛操日本美女一级片| www.自偷自拍.com| 亚洲黑人精品在线| 高潮久久久久久久久久久不卡| 丁香六月天网| a在线观看视频网站| 这个男人来自地球电影免费观看| 亚洲色图 男人天堂 中文字幕| 日韩免费高清中文字幕av| 亚洲精品美女久久av网站| 欧美中文综合在线视频| 色在线成人网| 99热网站在线观看| 天天躁夜夜躁狠狠躁躁| 欧美乱妇无乱码| 色老头精品视频在线观看| 51午夜福利影视在线观看| 亚洲欧美一区二区三区久久| 黄色视频不卡| 国精品久久久久久国模美| 精品乱码久久久久久99久播| 人妻久久中文字幕网| 人人妻,人人澡人人爽秒播| 国产国语露脸激情在线看| 另类精品久久| 亚洲精品自拍成人| www日本在线高清视频| 亚洲精品粉嫩美女一区| 久久影院123| 欧美 日韩 精品 国产| 成人手机av| 动漫黄色视频在线观看| 免费看十八禁软件| 激情在线观看视频在线高清 | 肉色欧美久久久久久久蜜桃| 啦啦啦中文免费视频观看日本| 日韩视频在线欧美| 老司机午夜十八禁免费视频| 欧美 日韩 精品 国产| 91国产中文字幕| 午夜两性在线视频| 1024视频免费在线观看| 在线观看66精品国产| 91av网站免费观看| 午夜91福利影院| 色94色欧美一区二区| 啦啦啦在线免费观看视频4| 亚洲视频免费观看视频| 多毛熟女@视频| av视频免费观看在线观看| 大码成人一级视频| 国产野战对白在线观看| 91国产中文字幕| 变态另类成人亚洲欧美熟女 | 国产区一区二久久| 嫩草影视91久久| 国产av又大| 悠悠久久av| 美女主播在线视频| 亚洲精品中文字幕在线视频| 大香蕉久久成人网| 亚洲熟女毛片儿| 亚洲国产欧美一区二区综合| netflix在线观看网站| 国产av精品麻豆| 99久久99久久久精品蜜桃| 精品久久久久久久毛片微露脸| 91精品国产国语对白视频| 黑丝袜美女国产一区| 久久精品国产综合久久久| 免费女性裸体啪啪无遮挡网站| 亚洲自偷自拍图片 自拍| 俄罗斯特黄特色一大片| 99九九在线精品视频| 精品国产一区二区三区久久久樱花| 男男h啪啪无遮挡| 大香蕉久久网| 日日摸夜夜添夜夜添小说| 色94色欧美一区二区| 色综合婷婷激情| 免费高清在线观看日韩| 精品高清国产在线一区| 久久国产精品大桥未久av| 91大片在线观看| 久久精品国产99精品国产亚洲性色 | tube8黄色片| 中文字幕色久视频| 天天躁夜夜躁狠狠躁躁| 国产免费现黄频在线看| 淫妇啪啪啪对白视频| 午夜福利乱码中文字幕| 91大片在线观看| 国产一区有黄有色的免费视频| 久久久久国产一级毛片高清牌| avwww免费| √禁漫天堂资源中文www| 精品国产一区二区三区四区第35| 青草久久国产| 99九九在线精品视频| 精品免费久久久久久久清纯 | 亚洲中文字幕日韩| 一区在线观看完整版| 欧美 日韩 精品 国产| bbb黄色大片| 精品国产超薄肉色丝袜足j| 一区在线观看完整版| 大型黄色视频在线免费观看| 国产麻豆69| 精品亚洲成国产av| 亚洲人成伊人成综合网2020| 极品少妇高潮喷水抽搐| 婷婷丁香在线五月| 精品视频人人做人人爽| 在线播放国产精品三级| 99久久人妻综合| 可以免费在线观看a视频的电影网站| 久久国产精品男人的天堂亚洲| svipshipincom国产片| 热99久久久久精品小说推荐| 国产成人av激情在线播放| 五月开心婷婷网| 咕卡用的链子| 极品教师在线免费播放| 不卡av一区二区三区| 成在线人永久免费视频| 我要看黄色一级片免费的| 国产精品av久久久久免费| 久久久久网色| 高清毛片免费观看视频网站 | 母亲3免费完整高清在线观看| 亚洲五月色婷婷综合| 欧美人与性动交α欧美精品济南到| 午夜视频精品福利| 日本撒尿小便嘘嘘汇集6| 免费少妇av软件| 午夜激情久久久久久久| 久久毛片免费看一区二区三区| 久久99热这里只频精品6学生| 岛国毛片在线播放| 9191精品国产免费久久| 欧美日韩亚洲综合一区二区三区_| 美女国产高潮福利片在线看| 十八禁人妻一区二区| 在线观看免费视频日本深夜| 久久香蕉激情| 免费黄频网站在线观看国产| 国产精品 欧美亚洲| 欧美黑人欧美精品刺激| 美女福利国产在线| 亚洲一码二码三码区别大吗| 亚洲成人免费电影在线观看|