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

    低緯度磁異常波數(shù)域化極方法機理分析與功能定位

    2020-08-18 08:01:24柴玉璞萬海珍
    石油地球物理勘探 2020年4期
    關(guān)鍵詞:低緯度波數(shù)算子

    柴玉璞 萬海珍

    (中國石油東方地球物理公司綜合物化探處,河北涿州072751)

    0 引言

    非垂直向下的磁化強度(矢量)和非垂直向下的磁異常場強度(矢量),致使磁異常形態(tài)復雜化,從而使磁異常與其場源的關(guān)系變得很不直觀。隨著磁傾角變小,此問題越來越嚴重,給磁異常解釋帶來許多困難。Baranov[1]提出的磁異?;瘶O概念(非垂直向下的磁化強度和磁異常場強度均轉(zhuǎn)換為垂直向下,就像將磁源體搬到北極所產(chǎn)生的磁異常場一樣,故稱化極),奠定了解決這一問題的理論基礎(chǔ);Bhattacharyya[2]提出的位場波數(shù)域轉(zhuǎn)換理論和方法,使磁異常化極得以實現(xiàn)。然而,隨著磁傾角的減小,化極過程不穩(wěn)定現(xiàn)象越發(fā)嚴重,使低緯度磁異?;瘶O效果很不理想。于是各種低緯度波數(shù)域化極方法相繼出現(xiàn)。這些方法可分為三類:①波數(shù)域壓制型化極法,如偽傾角化極法[3]、平方余弦壓制化極法[4]、維納濾波化極法[5]及雙曲正弦壓制化極法[6]等;②波數(shù)域反演化極法[7];③波數(shù)域偏移抽樣化極法[8-10]。

    本文主要討論三個問題:①從理論上證明壓制型化極法肯定不能用于低緯度磁異常定量解釋,但可用于定性解釋;②從理論上揭示,采用等效源法實現(xiàn)的波數(shù)域反演化極和采用理論化極算子實現(xiàn)的波數(shù)域反演化極,前者效果總是優(yōu)于后者的根本原因,給出該方法難以用于低緯度磁測資料定量解釋的功能定位;③闡明波數(shù)域偏移抽樣化極法的基本思想和方法原理,指出該方法相當苛刻的應(yīng)用條件。

    1 可用于定性解釋的低緯度化極方法——壓制型化極法

    從位場波數(shù)域轉(zhuǎn)換算法誤差方程[8-10]出發(fā),導出了實用含噪磁異常波數(shù)域化極算法誤差方程[9](附錄A)

    式中:Re(·)表示取實部;DFT(·)、IDFT(·)分別表示正、反離散傅里葉變換;R(·)、R′(·)分別為理論化極算子和實用化極算子;g′1(nΔ)、g2(nΔ)分別為含噪待化極磁異常場和理論化極磁異常場;n(n1,n2)、m(m1,m2)分別為空間域和波數(shù)域坐標矢量,n1=1,2,…,N1表示空間域x方向的離散坐標點,n2=1,2,…,N2表示空間域y方向的離散坐標點,m1=1,2,…,N1表示波數(shù)域u方向的離散坐標點,m2=1,2,…,N2表示波數(shù)域v方向的離散坐標點;是空間域采樣間隔對角矩陣,Δ1、Δ2分別為空間域x、y方向的采樣間隔;是波數(shù)域采樣間隔對角矩陣,d1、d2分別為波數(shù)域u、v方向的采樣間隔;表示波數(shù)域和空間域數(shù)據(jù)規(guī)模對角矩陣;e(nΔ)為噪聲,也稱觀測誤差。

    利用式(1),可分析各種壓制型化極法的基本功能。式中左端是以R′(md)為化極算子的實用化極算法表達式;右端第一項是理論化極異常場,第二項是化極損失項,第三項是觀測誤差。在低緯度區(qū),理論化極算子R(md)實際是一放大作用極強的扇通濾波器,其中心線方位角θ0=90°+D(D為磁偏角);實用化極算子R′(md)是R(md)經(jīng)壓制的結(jié)果,也是一放大作用很強的扇通濾波器,中心線與R(md)的中心線重合。數(shù)學推導和計算可證明

    也是一種濾波器。不過,式(2)對偽傾角化極法[3]和平方余弦壓制化極法[4]而言是扇通濾波器(圖1),對維納濾波化極法[5]和雙曲正弦壓制化極法[6]而言是準扇通濾波器。這四個濾波器的中心線都與R(md)、R′(md)的中心線重合,但濾波強度弱很多。

    筆者在研究磁異?;瘶O中的誤差規(guī)律時發(fā)現(xiàn),像理論化極算子R(md)這類扇通濾波器,可分解為一個方向低通濾波器(中心線方位角為θ0=90°+D)和兩個高通濾波器(其中濾波作用較強的濾波器中心方位角為θ0=90°+D,濾波作用較弱的濾波器中心方位角為θ0=D)[10]。其實實用化極算子R′(md)也可做同樣的分解。扇通濾波器可做同樣的分解,只是中心線方位角為θ0=D的濾波器消失,因為該扇通濾波器的極小值為零,而R(md)、R′(md)的極小值都不為零。

    圖1 扇通濾波器濾波等值線圖

    下面對分解出的各個濾波器的功能及其對化極結(jié)果的影響作分析。

    (3)濾波器R′(md)中較強的方向低通濾波器(濾波器3),會導致化極結(jié)果圖中出現(xiàn)沿磁偏角方向、正負相間的平行線狀等值線。這是濾波器3對隨機誤差所形成的、遍布全圖的正負等值線小圈沿磁偏角方向極度拉長的結(jié)果。因為隨機誤差所形成的小圈,高頻成分較弱,所以R′(md)中的兩個高通濾波器表現(xiàn)不明顯。

    上述三個濾波器的濾波強度都與壓制因子的強度相關(guān)。壓制因子增強,則濾波器3的作用減弱,濾波器1和濾波器2的作用增強。也就是說,壓制因子增強,會削弱隨機誤差對化極結(jié)果的干擾,但同時會使化極損失效應(yīng)增強。據(jù)此可得如下結(jié)論。

    (1)當增強壓制強度,使濾波器3所產(chǎn)生的平行線異常減弱到可忽略不計的情況下,濾波器2將大幅度降低化極結(jié)果中正異常幅度,而小幅度地增強磁偏角方向上的邊部負異常。

    (2)此種情況下,濾波器1對化極結(jié)果在垂直于磁偏角方向上有一定程度的圓滑作用。

    (3)在上述條件下,濾波器1、濾波器2將一定程度地影響化極異常的形態(tài),但不影響化極異常的中心,也基本不影響化極結(jié)果中正異常的分布范圍。

    為了驗證上述理論的正確性,設(shè)計一方柱體組合模型進行試驗,模型參數(shù)見文獻[11]。模型的磁傾角和磁偏角均為0°,觀測范圍為50km×50km,點、線距均為1km。圖2是該模型的總磁異常(ΔT),圖3是磁傾角為90°的理論化極異常。圖4是采用偽傾角法[3]計算的化極結(jié)果,其余三種壓制型化極方法的化極結(jié)果與偽傾角法化極結(jié)果大同小異,故本文未展示。粗看上去,圖3與圖4差異非常大:圖4中化極結(jié)果的正異常幅值比圖3小得多,而磁偏角方向上的邊部負異常則比圖3強;此外,垂直于磁偏角方向上的異常等值線出現(xiàn)一定程度的圓滑,而且變得稀疏。然而仔細觀察可以發(fā)現(xiàn),圖4所示的化極結(jié)果中正異常的中心位置和分布范圍與圖3基本一致。

    圖2 方柱體組合模型磁源的總磁異常(ΔT)分布圖

    理論分析和模型計算結(jié)果表明:各種波數(shù)域壓制型化極法不能用于低緯度磁異常數(shù)據(jù)的定量解釋,因為這些方法的化極結(jié)果與理論化極結(jié)果,無論異常形態(tài)還是量值,差異都相當大;但是這些方法完全可用于定性解釋,因為化極結(jié)果中正異常的中心位置和分布范圍與理論化極異常的基本一致。

    圖5是秘魯A地區(qū)實測磁異常(ΔT)。該地區(qū)的地磁傾角和磁偏角分別為-2°和0°,觀測點、線距均為2km,資料覆蓋范圍為100km×288km。圖6是圖5采用偽傾角法的化極結(jié)果,該圖為秘魯A地區(qū)的地質(zhì)勘探提供了重要信息。

    圖3 圖2的理論化極異常(磁傾角為90°)

    圖4 采用偽傾角法對圖2磁異常的化極結(jié)果

    圖5 秘魯A地區(qū)實測磁異常圖

    圖6 利用偽傾角法對圖5的化極異常圖

    2 有可能用于定量解釋的低緯度化極方法:波數(shù)域反演化極法和波數(shù)域偏移抽樣化極法

    2.1 波數(shù)域反演化極法

    從原理上看,波數(shù)域反演化極法[7]可用于低緯度磁異常數(shù)據(jù)的定量解釋。該方法的核心是求解依據(jù)基本關(guān)系

    建立的方程組

    式中:G是以1/R(md)為元素的對角矩陣;m′和d′分別為模型向量和觀測值向量;W為加權(quán)函數(shù)矩陣;S為加權(quán)差分矩陣;上標“H”和“T”分別表示復共軛轉(zhuǎn)置和轉(zhuǎn)置;μ表示正則化參數(shù)。由理論化極算法誤差方程[10]可知,式(3)嚴格成立的條件是G以T(md)為元素,而不是1/R(md),因為式(3)中的d′和m′分別是待化極異常的有限離散譜和化極目標異常的有限離散譜。T(md)與1/R(md)的關(guān)系(推導過程詳見附錄A)為

    Li等[7]用1/R(md)作為G的元素所獲得的化極結(jié)果不理想,改用等效源法計算G的元素,結(jié)果大為改善。式(5)從理論上解釋了這一現(xiàn)象產(chǎn)生的原因:采用等效源法得到的G中,至少包含了計算窗外數(shù)據(jù)的部分影響,而采用1/R(md)計算的G,則完全忽略了計算窗外數(shù)據(jù)的影響。文獻[7]中的模型算例是對一直立方柱體的斜磁化磁異常場進行化極,所以采用一小直立方柱體為等效源計算G的元素,化極效果相當好。因為小方柱體磁異常的窗外數(shù)據(jù)與大方柱體磁異常的窗數(shù)據(jù)有很強的相似性。

    由理論化極算法誤差方程導出的式(5),為波數(shù)域反演化極過程中采用等效源法計算G的元素,提供了理論依據(jù)。然而對實際資料而言,窗外數(shù)據(jù)難以較精確地模擬。這可能就是該方法難以應(yīng)用于實際資料處理的根本原因。

    2.2 波數(shù)域偏移抽樣化極法

    由含噪實用化極算法方程的推導過程(附錄A)可知:式(1)中,若將DFT替換為移樣離散傅里葉正變換算子DFT0η,將IDFT替換為移樣離散傅里葉反變換算子IDFT0η(下標“0η”表示僅波數(shù)域偏移抽樣,偏移參數(shù)η為二維行向量),將m替換為m+η,并令R′[(m+η)d]=R[(m+η)d],即得波數(shù)域偏移抽樣化極算法方程

    式中,左端為化極算法表達式,右端第一項為理論化極結(jié)果項,第二項為隨機誤差。波數(shù)域偏移抽樣化極方法的基本思想是:測線沿磁偏角方向布設(shè),使磁偏角為0°,然后通過沿u軸方向偏移v軸半個采樣間隔(0.5d1)抽樣,大大減小化極算子的量值,以控制隨機誤差對化極結(jié)果的影響。

    由上述基本思想和式(6)的化極算法表達式可知,波數(shù)域偏移抽樣化極法的具體做法包括如下步驟:

    (1)對磁偏角為0°的磁測資料做DFT0η變換,偏移抽樣參數(shù)η=(0.5,0);

    (2)將步驟(1)的結(jié)果乘以理論化極算子R[(m+η)d];

    (3)對步驟(2)的結(jié)果做IDFT0η變換,偏移參數(shù)η=(0.5,0);

    (4)對步驟(3)的結(jié)果取實部,即得偏移抽樣化極結(jié)果。

    波數(shù)域偏移抽樣化極算法方程右端只有兩項,一項是理論化極結(jié)果項,一項是觀測誤差項。因此把觀測誤差項的量值減至可忽略的程度,是該方法成功的關(guān)鍵。而誤差項量值的減小,是通過沿u軸方向偏移v軸半個采樣間隔(0.5d1)抽樣實現(xiàn)的(d1越大,減小量越大)。這種抽樣方式可使化極算子的量值大大減小,而且隨著N1(d1=1/N1Δ1)的減小,量值還會進一步減小。然而當沿磁偏角方向上的網(wǎng)格數(shù)N1降至64時,化極算子的量值仍然相當大,以至于要求資料相對精度高于0.2%,才能保證觀測誤差項量值小到可忽略不計。

    由上述分析可以看出,波數(shù)域偏移抽樣化極法可用于低緯度磁異常的定量解釋,但條件相當苛刻:磁偏角方向上網(wǎng)格數(shù)受限,且資料精度要求很高。因此在目前資料精度條件下,波數(shù)域偏移抽樣化極法不可能用于低緯度區(qū)大面積磁測資料的定量解釋。隨著資料精度的不斷提高,將來該方法有可能用于低緯度區(qū)小面積高精度磁測資料的定量解釋。

    3 結(jié)論

    (1)理論研究表明,各種波數(shù)域壓制型化極法不能應(yīng)用于低緯度磁測資料的定量解釋,但可以成為低緯度磁測資料定性解釋的有效工具。這一理論成果的獲得,實用化極算法誤差方程的導出是基礎(chǔ),扇通濾波器定性分解(可分解為一個方向低通濾波器和兩個或一個高通濾波器)觀點的提出和運用是關(guān)鍵。這一源于化極問題的原創(chuàng)性概念和觀點,為其他涉及扇通濾波器問題的解決提供了新的思路和方法。

    (2)文中的式(5)為波數(shù)域反演化極中采用等效源法計算G的元素,提供了嚴格的理論依據(jù)。理論上,該方法可用于低緯度磁測資料的定量解釋,然而實際資料的窗外數(shù)據(jù)難以較精確地模擬。因此,該方法難以應(yīng)用于實際數(shù)據(jù)的化極。該方法需求解超大型方程組,可能是難以實用化的另一原因。

    (3)波數(shù)域偏移抽樣化極法,理論上可用于低緯度磁測資料的定量解釋,但由于應(yīng)用條件過于苛刻,在目前數(shù)據(jù)精度條件下不可能用于低緯度、大面積磁測資料的定量解釋。但是隨著未來資料精度的提高,該方法有可能用于低緯度區(qū)、小面積、高精度磁測資料的定量解釋。

    感謝中國石油集團東方地球物理公司賈繼軍高級工程師在研究初期在編程方面給予的幫助,楊戰(zhàn)軍高級工程師在圖件繪制方面的幫助,以及孫喜明高級工程師在資料準備方面提供的幫助!

    附錄A 實用化極算法方程及T(md)表達式的推導

    假設(shè)g1(x,y)和G1(u,v)分別表示無噪待化極異常及其譜,g2(x,y)和G2(u,v)分別表示化極目標異常及其譜,那么g2(x,y)和G2(u,v)的二維IDFT算法誤差方程[9]為

    g1(x,y)和G1(u,v)的二維DFT算法誤差方程[9]可寫作

    式(A-2)兩邊乘以R(md),并利用關(guān)系

    將式(A-4)代入式(A-1),并移項,得理論化極算法誤差方程

    在采用實用化極算子的情況下,式(A-6)右端前四項量值很小,可忽略不計。所以實用無噪磁異?;瘶O算法誤差方程簡化為

    并將g2(nΔ)由方程左端移到右端,即得正文中的式(1),即含噪實用化極算法方程。

    對式(A-5)做如下變換:①將左端的g2(nΔ)移到右端;②對方程兩端作DFT;③方程兩端同除以R(md)DFT[g2(nΔ)],可得正文中的式(5),即T(md)的表達式。

    猜你喜歡
    低緯度波數(shù)算子
    聲場波數(shù)積分截斷波數(shù)自適應(yīng)選取方法
    一種基于SOM神經(jīng)網(wǎng)絡(luò)中藥材分類識別系統(tǒng)
    電子測試(2022年16期)2022-10-17 09:32:26
    擬微分算子在Hp(ω)上的有界性
    基于FY-3C/MWHTS觀測資料反演中低緯度海面氣壓
    各向異性次Laplace算子和擬p-次Laplace算子的Picone恒等式及其應(yīng)用
    高海拔低緯度條件下的8個家蠶品種比較試驗
    一類Markov模算子半群與相應(yīng)的算子值Dirichlet型刻畫
    低緯度海上平臺天然氣冷卻系統(tǒng)的研究
    Roper-Suffridge延拓算子與Loewner鏈
    低緯度、短基線大落差、高原山區(qū)網(wǎng)絡(luò)RTK技術(shù)的實現(xiàn)與應(yīng)用*
    地礦測繪(2015年2期)2015-04-19 07:25:56
    国产成人免费无遮挡视频| 晚上一个人看的免费电影| 免费观看a级毛片全部| 亚洲精华国产精华液的使用体验| 黄色视频在线播放观看不卡| 亚洲一卡2卡3卡4卡5卡精品中文| 纯流量卡能插随身wifi吗| 青青草视频在线视频观看| 日韩大码丰满熟妇| 妹子高潮喷水视频| 嫩草影院入口| 女人精品久久久久毛片| 日本午夜av视频| 汤姆久久久久久久影院中文字幕| 亚洲成人av在线免费| 男女午夜视频在线观看| 男人爽女人下面视频在线观看| 亚洲 欧美一区二区三区| 观看av在线不卡| 欧美精品一区二区免费开放| 国产伦理片在线播放av一区| av在线老鸭窝| 婷婷成人精品国产| 黄色视频在线播放观看不卡| 亚洲情色 制服丝袜| 国产精品偷伦视频观看了| 考比视频在线观看| 精品亚洲乱码少妇综合久久| 高清欧美精品videossex| 精品午夜福利在线看| 中文字幕精品免费在线观看视频| 69精品国产乱码久久久| 黄片小视频在线播放| 少妇人妻 视频| 国产欧美亚洲国产| www.自偷自拍.com| 亚洲一区中文字幕在线| 亚洲国产欧美在线一区| 色精品久久人妻99蜜桃| av在线老鸭窝| 国产亚洲av高清不卡| 亚洲欧美一区二区三区久久| 少妇猛男粗大的猛烈进出视频| 国产精品 国内视频| 天天躁夜夜躁狠狠躁躁| 少妇人妻 视频| 两个人免费观看高清视频| 极品人妻少妇av视频| 一个人免费看片子| 男男h啪啪无遮挡| 成人毛片60女人毛片免费| 97人妻天天添夜夜摸| 成人黄色视频免费在线看| 黑人巨大精品欧美一区二区蜜桃| 国产精品亚洲av一区麻豆 | 国产成人系列免费观看| 五月开心婷婷网| 久久久久精品国产欧美久久久 | 免费日韩欧美在线观看| 久久久久久久久久久久大奶| 1024香蕉在线观看| 国产精品麻豆人妻色哟哟久久| 亚洲国产成人一精品久久久| 国产伦人伦偷精品视频| 久久久精品免费免费高清| 亚洲av成人不卡在线观看播放网 | 亚洲成人一二三区av| 一本色道久久久久久精品综合| 人人妻人人爽人人添夜夜欢视频| 搡老乐熟女国产| 人人妻人人澡人人看| 亚洲,欧美,日韩| 国产一区二区激情短视频 | 人人妻人人爽人人添夜夜欢视频| 亚洲天堂av无毛| 国产99久久九九免费精品| 亚洲成人免费av在线播放| 欧美精品人与动牲交sv欧美| 好男人视频免费观看在线| 黄片小视频在线播放| 男女边摸边吃奶| 咕卡用的链子| 涩涩av久久男人的天堂| e午夜精品久久久久久久| 最近手机中文字幕大全| 国产精品二区激情视频| 伦理电影免费视频| 免费看不卡的av| 一区二区三区四区激情视频| 中文字幕人妻丝袜一区二区 | 校园人妻丝袜中文字幕| 日韩欧美一区视频在线观看| 亚洲自偷自拍图片 自拍| 美女主播在线视频| 免费黄色在线免费观看| 黄频高清免费视频| 日韩一本色道免费dvd| 美女主播在线视频| 免费黄频网站在线观看国产| 日韩一卡2卡3卡4卡2021年| 免费日韩欧美在线观看| 中文字幕人妻丝袜制服| 久久综合国产亚洲精品| 亚洲图色成人| 日韩不卡一区二区三区视频在线| 久久久久国产精品人妻一区二区| 99国产精品免费福利视频| 国产一区亚洲一区在线观看| 国产精品国产三级国产专区5o| 国产精品.久久久| 多毛熟女@视频| 搡老岳熟女国产| 欧美精品人与动牲交sv欧美| 精品一区二区三区av网在线观看 | 精品少妇黑人巨大在线播放| netflix在线观看网站| 欧美在线一区亚洲| 黄色一级大片看看| 少妇精品久久久久久久| 老司机在亚洲福利影院| 久久热在线av| 国产成人免费观看mmmm| 色网站视频免费| avwww免费| 悠悠久久av| 好男人视频免费观看在线| 国产成人91sexporn| 又大又爽又粗| 免费高清在线观看日韩| 下体分泌物呈黄色| 成人漫画全彩无遮挡| 国产精品免费大片| 国产 一区精品| 大话2 男鬼变身卡| 国产一区二区三区综合在线观看| 午夜激情av网站| 卡戴珊不雅视频在线播放| 黄片播放在线免费| 亚洲国产欧美在线一区| 欧美日韩成人在线一区二区| 新久久久久国产一级毛片| 2021少妇久久久久久久久久久| 中文精品一卡2卡3卡4更新| 国产精品一区二区在线观看99| 欧美成人精品欧美一级黄| 久久精品久久精品一区二区三区| 少妇被粗大猛烈的视频| 国产精品.久久久| 欧美av亚洲av综合av国产av | 亚洲五月色婷婷综合| 啦啦啦在线免费观看视频4| 日韩av在线免费看完整版不卡| 国产伦理片在线播放av一区| 韩国精品一区二区三区| 亚洲国产精品成人久久小说| 久久99热这里只频精品6学生| 一本—道久久a久久精品蜜桃钙片| 成人黄色视频免费在线看| 国产又色又爽无遮挡免| 精品亚洲乱码少妇综合久久| 丝袜脚勾引网站| 亚洲四区av| 国语对白做爰xxxⅹ性视频网站| 成年美女黄网站色视频大全免费| 亚洲在久久综合| 一区二区日韩欧美中文字幕| 国产乱来视频区| 一级黄片播放器| 久久久久精品性色| 一级毛片我不卡| 亚洲,一卡二卡三卡| 两个人看的免费小视频| 欧美国产精品一级二级三级| 国产精品成人在线| 亚洲激情五月婷婷啪啪| 亚洲国产精品一区二区三区在线| 人人妻人人爽人人添夜夜欢视频| 精品视频人人做人人爽| 男人舔女人的私密视频| 美女视频免费永久观看网站| 亚洲精品国产av成人精品| 搡老乐熟女国产| svipshipincom国产片| 搡老乐熟女国产| 国产亚洲午夜精品一区二区久久| 日本一区二区免费在线视频| 久久精品亚洲av国产电影网| 国产熟女欧美一区二区| 黄色视频不卡| 一本一本久久a久久精品综合妖精| 欧美变态另类bdsm刘玥| 亚洲图色成人| 国产日韩欧美亚洲二区| 亚洲五月色婷婷综合| 久久精品国产综合久久久| 黄色毛片三级朝国网站| 校园人妻丝袜中文字幕| 欧美久久黑人一区二区| 免费不卡黄色视频| 爱豆传媒免费全集在线观看| 国产男女内射视频| www.精华液| 成人漫画全彩无遮挡| 男女国产视频网站| 国产精品欧美亚洲77777| 天美传媒精品一区二区| av天堂久久9| 精品一区在线观看国产| 国产成人免费观看mmmm| 熟女少妇亚洲综合色aaa.| 2018国产大陆天天弄谢| 街头女战士在线观看网站| 中文字幕人妻丝袜制服| 日本一区二区免费在线视频| 日本午夜av视频| 最近最新中文字幕免费大全7| 国产成人91sexporn| 日韩精品有码人妻一区| 波多野结衣一区麻豆| 人体艺术视频欧美日本| 亚洲,欧美精品.| 亚洲av日韩精品久久久久久密 | 老司机影院毛片| 中文字幕另类日韩欧美亚洲嫩草| 欧美日韩亚洲综合一区二区三区_| 亚洲av日韩精品久久久久久密 | 又黄又粗又硬又大视频| 青青草视频在线视频观看| 精品一区在线观看国产| av福利片在线| 亚洲伊人久久精品综合| 欧美日韩视频高清一区二区三区二| 精品午夜福利在线看| 国产成人欧美在线观看 | h视频一区二区三区| tube8黄色片| 久久久精品免费免费高清| 国产乱人偷精品视频| 亚洲,欧美精品.| 国产一区二区三区综合在线观看| 中文字幕人妻丝袜一区二区 | 波野结衣二区三区在线| 国产毛片在线视频| 精品人妻熟女毛片av久久网站| 少妇被粗大的猛进出69影院| 欧美 日韩 精品 国产| 成人漫画全彩无遮挡| 中国国产av一级| 亚洲欧美一区二区三区国产| netflix在线观看网站| 中文字幕精品免费在线观看视频| 精品一区在线观看国产| 中文字幕人妻丝袜制服| 亚洲国产日韩一区二区| 女人被躁到高潮嗷嗷叫费观| 久久99热这里只频精品6学生| 亚洲国产欧美在线一区| 免费久久久久久久精品成人欧美视频| 多毛熟女@视频| 又粗又硬又长又爽又黄的视频| 少妇被粗大的猛进出69影院| 老司机影院毛片| 国产精品久久久久久人妻精品电影 | 在线观看一区二区三区激情| 久久国产精品男人的天堂亚洲| 熟女av电影| 亚洲av日韩精品久久久久久密 | 中国三级夫妇交换| 91aial.com中文字幕在线观看| 女性生殖器流出的白浆| 人体艺术视频欧美日本| 深夜精品福利| 视频在线观看一区二区三区| 久久99一区二区三区| 满18在线观看网站| 国产黄频视频在线观看| 两性夫妻黄色片| 成人国产麻豆网| 亚洲国产欧美在线一区| 精品卡一卡二卡四卡免费| 日日摸夜夜添夜夜爱| 9热在线视频观看99| 在线免费观看不下载黄p国产| 欧美日韩av久久| 免费看av在线观看网站| 秋霞在线观看毛片| 自线自在国产av| 国产成人精品福利久久| 欧美日韩综合久久久久久| 国产精品.久久久| 九草在线视频观看| 十八禁网站网址无遮挡| 一级片免费观看大全| 亚洲国产欧美一区二区综合| 男女边摸边吃奶| 91aial.com中文字幕在线观看| 欧美日韩一级在线毛片| 少妇被粗大猛烈的视频| 操出白浆在线播放| 久久久久人妻精品一区果冻| 久久久久久久精品精品| 久久韩国三级中文字幕| 又黄又粗又硬又大视频| 久久精品久久精品一区二区三区| 免费观看性生交大片5| 免费av中文字幕在线| 久久久久精品国产欧美久久久 | 国产高清不卡午夜福利| 十八禁网站网址无遮挡| 中文字幕另类日韩欧美亚洲嫩草| 精品国产乱码久久久久久男人| 美女脱内裤让男人舔精品视频| 1024视频免费在线观看| 18禁国产床啪视频网站| 久久精品人人爽人人爽视色| 日韩电影二区| 熟妇人妻不卡中文字幕| 亚洲图色成人| 成年动漫av网址| 久久精品人人爽人人爽视色| 国产av精品麻豆| 欧美日韩一区二区视频在线观看视频在线| 亚洲成国产人片在线观看| xxx大片免费视频| 久久99热这里只频精品6学生| 国产免费一区二区三区四区乱码| 在线 av 中文字幕| 另类亚洲欧美激情| 女性生殖器流出的白浆| 免费高清在线观看日韩| 国产又色又爽无遮挡免| 电影成人av| 精品一区二区三区四区五区乱码 | 高清黄色对白视频在线免费看| 亚洲精品aⅴ在线观看| 一区二区三区精品91| 精品久久久久久电影网| www.熟女人妻精品国产| 精品第一国产精品| 最近的中文字幕免费完整| 18禁观看日本| 精品一品国产午夜福利视频| 捣出白浆h1v1| 日韩一本色道免费dvd| 日本色播在线视频| 日本猛色少妇xxxxx猛交久久| 国产国语露脸激情在线看| 美女福利国产在线| 飞空精品影院首页| 精品一区二区三卡| 国产成人精品福利久久| 欧美精品人与动牲交sv欧美| 亚洲美女视频黄频| 丰满乱子伦码专区| 婷婷成人精品国产| 亚洲成人免费av在线播放| 青春草视频在线免费观看| 欧美av亚洲av综合av国产av | 国产在线免费精品| 叶爱在线成人免费视频播放| 日韩av不卡免费在线播放| 看非洲黑人一级黄片| 色婷婷av一区二区三区视频| 十分钟在线观看高清视频www| 国产成人啪精品午夜网站| 久久久精品94久久精品| 国产片内射在线| 老司机在亚洲福利影院| 亚洲三区欧美一区| 十八禁高潮呻吟视频| 日韩视频在线欧美| 成人亚洲精品一区在线观看| 黄频高清免费视频| 满18在线观看网站| 满18在线观看网站| 美女主播在线视频| 天美传媒精品一区二区| 国产人伦9x9x在线观看| 午夜久久久在线观看| 欧美日韩一级在线毛片| 精品久久久久久电影网| 中文字幕最新亚洲高清| 日日摸夜夜添夜夜爱| 午夜激情av网站| 少妇精品久久久久久久| 日本黄色日本黄色录像| 在线看a的网站| 日本色播在线视频| 夫妻午夜视频| 一级毛片我不卡| 亚洲,欧美精品.| 欧美在线一区亚洲| 七月丁香在线播放| 黄色怎么调成土黄色| 国产成人a∨麻豆精品| 久久国产亚洲av麻豆专区| 婷婷色综合www| 亚洲av成人不卡在线观看播放网 | 亚洲美女搞黄在线观看| 美国免费a级毛片| av女优亚洲男人天堂| 亚洲国产欧美日韩在线播放| 成年人午夜在线观看视频| 亚洲三区欧美一区| 亚洲国产精品国产精品| 老司机亚洲免费影院| 亚洲欧美清纯卡通| 伊人久久大香线蕉亚洲五| 无遮挡黄片免费观看| av视频免费观看在线观看| 亚洲国产欧美日韩在线播放| 精品人妻一区二区三区麻豆| 19禁男女啪啪无遮挡网站| 天天躁夜夜躁狠狠久久av| 视频区图区小说| 日日爽夜夜爽网站| a 毛片基地| 伊人久久国产一区二区| 精品国产一区二区三区久久久樱花| 国产精品成人在线| 国产日韩一区二区三区精品不卡| 最近手机中文字幕大全| 晚上一个人看的免费电影| 亚洲久久久国产精品| 欧美日韩福利视频一区二区| 成人黄色视频免费在线看| 免费高清在线观看视频在线观看| 9色porny在线观看| 一边摸一边做爽爽视频免费| 9191精品国产免费久久| 国产黄色视频一区二区在线观看| 乱人伦中国视频| 亚洲男人天堂网一区| 日本av免费视频播放| 激情五月婷婷亚洲| 免费黄色在线免费观看| 18禁观看日本| 老汉色av国产亚洲站长工具| 久久女婷五月综合色啪小说| 中文天堂在线官网| 一二三四中文在线观看免费高清| 亚洲精品在线美女| 国产不卡av网站在线观看| 天堂中文最新版在线下载| 香蕉国产在线看| 天美传媒精品一区二区| 欧美国产精品va在线观看不卡| 少妇 在线观看| 国产免费现黄频在线看| 欧美在线一区亚洲| 69精品国产乱码久久久| 国产亚洲最大av| 欧美日韩精品网址| 精品国产一区二区三区四区第35| 黑丝袜美女国产一区| 在线天堂中文资源库| 久热这里只有精品99| av免费观看日本| 日韩成人av中文字幕在线观看| 欧美亚洲 丝袜 人妻 在线| 免费观看a级毛片全部| 9色porny在线观看| 天天躁夜夜躁狠狠躁躁| 人人妻人人爽人人添夜夜欢视频| 黄色一级大片看看| 亚洲免费av在线视频| 久久99热这里只频精品6学生| 国产黄色视频一区二区在线观看| 亚洲成av片中文字幕在线观看| 久久天躁狠狠躁夜夜2o2o | 美女视频免费永久观看网站| 欧美黑人精品巨大| 婷婷色麻豆天堂久久| 亚洲精品乱久久久久久| 毛片一级片免费看久久久久| 男人舔女人的私密视频| 精品人妻熟女毛片av久久网站| 国产黄频视频在线观看| 大片电影免费在线观看免费| 亚洲精品久久午夜乱码| 女人被躁到高潮嗷嗷叫费观| 男女床上黄色一级片免费看| 成人黄色视频免费在线看| 日韩av免费高清视频| 国产精品熟女久久久久浪| 伊人久久国产一区二区| 叶爱在线成人免费视频播放| 国产精品久久久久久人妻精品电影 | 18禁观看日本| 精品亚洲成国产av| 男女之事视频高清在线观看 | 免费观看性生交大片5| 高清欧美精品videossex| 一边摸一边做爽爽视频免费| 亚洲精品成人av观看孕妇| 日韩一区二区三区影片| 国产乱来视频区| 久久久久久人人人人人| 亚洲伊人色综图| 欧美日韩亚洲高清精品| a 毛片基地| 日韩制服丝袜自拍偷拍| 久久久久久久国产电影| 日本黄色日本黄色录像| 人人妻人人澡人人爽人人夜夜| 最新的欧美精品一区二区| 啦啦啦在线观看免费高清www| 免费在线观看视频国产中文字幕亚洲 | 成年人免费黄色播放视频| 青草久久国产| 精品午夜福利在线看| 亚洲精品国产av成人精品| tube8黄色片| 午夜日本视频在线| 亚洲精品av麻豆狂野| 国产精品嫩草影院av在线观看| 国产日韩欧美在线精品| 亚洲精品美女久久av网站| 亚洲av国产av综合av卡| 日韩欧美精品免费久久| 久久精品国产亚洲av涩爱| 在线观看免费午夜福利视频| 欧美另类一区| 日韩视频在线欧美| 国产男人的电影天堂91| 在线观看人妻少妇| 亚洲免费av在线视频| 少妇猛男粗大的猛烈进出视频| 久久毛片免费看一区二区三区| 亚洲色图 男人天堂 中文字幕| 丁香六月欧美| 看免费成人av毛片| 免费观看性生交大片5| 亚洲欧美成人精品一区二区| 高清av免费在线| 欧美97在线视频| 日本猛色少妇xxxxx猛交久久| 搡老岳熟女国产| 午夜激情久久久久久久| 999久久久国产精品视频| 国产 一区精品| 亚洲精品中文字幕在线视频| 国产成人精品久久久久久| 在线观看国产h片| 亚洲久久久国产精品| 国产片内射在线| kizo精华| 亚洲人成电影观看| 国产成人免费无遮挡视频| 妹子高潮喷水视频| 亚洲欧美精品自产自拍| av网站在线播放免费| 天美传媒精品一区二区| 欧美成人精品欧美一级黄| 亚洲人成网站在线观看播放| 亚洲精品,欧美精品| 久久久久人妻精品一区果冻| 亚洲欧美中文字幕日韩二区| 久久精品aⅴ一区二区三区四区| 久久天躁狠狠躁夜夜2o2o | 另类亚洲欧美激情| 精品国产乱码久久久久久小说| 黑人欧美特级aaaaaa片| a 毛片基地| 亚洲,欧美,日韩| 在线看a的网站| 色视频在线一区二区三区| 欧美黑人精品巨大| 这个男人来自地球电影免费观看 | 蜜桃在线观看..| 免费人妻精品一区二区三区视频| 欧美中文综合在线视频| 亚洲欧美中文字幕日韩二区| 免费黄频网站在线观看国产| 天堂中文最新版在线下载| av又黄又爽大尺度在线免费看| 久久 成人 亚洲| tube8黄色片| 十八禁人妻一区二区| 男女边摸边吃奶| 少妇人妻 视频| 中文字幕高清在线视频| av卡一久久| 精品亚洲乱码少妇综合久久| 亚洲成人国产一区在线观看 | 美女脱内裤让男人舔精品视频| 国产免费福利视频在线观看| 国产精品偷伦视频观看了| 日韩中文字幕视频在线看片| 91精品伊人久久大香线蕉| 中文字幕精品免费在线观看视频| 丝袜脚勾引网站| 日韩精品免费视频一区二区三区| 色网站视频免费| 亚洲av成人精品一二三区| 男女边吃奶边做爰视频| 777米奇影视久久| 国产精品一区二区精品视频观看| 在线观看免费高清a一片| 热99国产精品久久久久久7| 女人高潮潮喷娇喘18禁视频| www.av在线官网国产| 久久久亚洲精品成人影院| 少妇人妻久久综合中文| 亚洲精品国产一区二区精华液| 美女福利国产在线| 午夜免费男女啪啪视频观看| 亚洲欧洲国产日韩| 热re99久久国产66热| 最近的中文字幕免费完整| videosex国产| 欧美少妇被猛烈插入视频| 国产精品久久久av美女十八| 男女免费视频国产| 大片免费播放器 马上看|