• <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
    在线观看舔阴道视频| 免费一级毛片在线播放高清视频| bbb黄色大片| 天堂影院成人在线观看| 久久国产精品影院| 久久精品国产清高在天天线| 一级毛片高清免费大全| 成人18禁在线播放| 这个男人来自地球电影免费观看| 久久久久久久精品吃奶| 丁香欧美五月| 欧美日本视频| 国产精品爽爽va在线观看网站| 日韩免费av在线播放| 成人高潮视频无遮挡免费网站| 久久久精品大字幕| 熟女少妇亚洲综合色aaa.| 三级毛片av免费| www.www免费av| 搡老妇女老女人老熟妇| 俄罗斯特黄特色一大片| 国产黄a三级三级三级人| 久久天躁狠狠躁夜夜2o2o| 少妇的丰满在线观看| 久久久久久久久免费视频了| 麻豆成人av在线观看| 亚洲精品一区av在线观看| 免费高清视频大片| 色播亚洲综合网| 国产成人aa在线观看| 久久热在线av| 身体一侧抽搐| 中文字幕人妻丝袜一区二区| 国产成人影院久久av| 免费看日本二区| 99国产精品99久久久久| 国产三级中文精品| 国产黄a三级三级三级人| 毛片女人毛片| 国产成人av激情在线播放| 一级黄色大片毛片| 香蕉av资源在线| 欧美zozozo另类| 国产精品亚洲av一区麻豆| 真人做人爱边吃奶动态| 91字幕亚洲| 国产精品一区二区三区四区久久| 精品欧美国产一区二区三| 国产精品 欧美亚洲| 91九色精品人成在线观看| 国产爱豆传媒在线观看| 亚洲av日韩精品久久久久久密| 午夜福利在线观看吧| 国产v大片淫在线免费观看| 真人做人爱边吃奶动态| 人人妻,人人澡人人爽秒播| www.999成人在线观看| 又大又爽又粗| 淫妇啪啪啪对白视频| 麻豆成人av在线观看| 精品国产三级普通话版| 18美女黄网站色大片免费观看| 成熟少妇高潮喷水视频| 黄色日韩在线| 久久天堂一区二区三区四区| 男人的好看免费观看在线视频| 丰满人妻熟妇乱又伦精品不卡| 岛国在线观看网站| 精品国产三级普通话版| www.www免费av| 高清在线国产一区| 日韩中文字幕欧美一区二区| 波多野结衣高清作品| 免费观看的影片在线观看| 中文字幕最新亚洲高清| 色哟哟哟哟哟哟| 99热只有精品国产| 精品一区二区三区视频在线 | 国产成人精品无人区| e午夜精品久久久久久久| 国产成人精品久久二区二区免费| 最好的美女福利视频网| 欧美成人性av电影在线观看| 老司机福利观看| 丁香欧美五月| 国产精品自产拍在线观看55亚洲| x7x7x7水蜜桃| 亚洲最大成人中文| 久久99热这里只有精品18| 国产精品 国内视频| 精品久久久久久成人av| 在线观看免费视频日本深夜| 天堂动漫精品| 久久这里只有精品19| 99国产精品99久久久久| 欧美一区二区精品小视频在线| 亚洲乱码一区二区免费版| 九色成人免费人妻av| 五月玫瑰六月丁香| 美女 人体艺术 gogo| 国产成人精品久久二区二区免费| 亚洲欧美日韩无卡精品| 亚洲国产欧美一区二区综合| 国产在线精品亚洲第一网站| 日韩人妻高清精品专区| 少妇裸体淫交视频免费看高清| 免费看日本二区| 欧美激情久久久久久爽电影| 免费电影在线观看免费观看| 亚洲人成电影免费在线| 亚洲无线观看免费| 久久精品aⅴ一区二区三区四区| 男女视频在线观看网站免费| 久久久久久久久免费视频了| 日本五十路高清| 亚洲狠狠婷婷综合久久图片| 精品99又大又爽又粗少妇毛片 | 桃红色精品国产亚洲av| 脱女人内裤的视频| 午夜福利在线在线| 91av网站免费观看| 黄色丝袜av网址大全| 久久国产精品影院| 国产精品日韩av在线免费观看| 美女 人体艺术 gogo| 欧美日本视频| 国内揄拍国产精品人妻在线| 香蕉国产在线看| 日本一本二区三区精品| 神马国产精品三级电影在线观看| 亚洲av熟女| 99热这里只有精品一区 | 成在线人永久免费视频| 亚洲国产中文字幕在线视频| 悠悠久久av| 国产一区二区激情短视频| 久99久视频精品免费| 成人精品一区二区免费| 日韩 欧美 亚洲 中文字幕| 啪啪无遮挡十八禁网站| 欧美一级a爱片免费观看看| av在线蜜桃| 国产精品日韩av在线免费观看| 精品一区二区三区四区五区乱码| 婷婷亚洲欧美| 日韩欧美三级三区| 99久久精品国产亚洲精品| 亚洲av成人一区二区三| 久久这里只有精品中国| 欧美色视频一区免费| 免费在线观看日本一区| 国产三级在线视频| 亚洲精华国产精华精| 国产高清videossex| 亚洲电影在线观看av| www.999成人在线观看| 久99久视频精品免费| 人人妻人人看人人澡| 女人高潮潮喷娇喘18禁视频| 色精品久久人妻99蜜桃| 国产伦精品一区二区三区四那| 99国产极品粉嫩在线观看| 我要搜黄色片| 制服人妻中文乱码| av天堂在线播放| 免费一级毛片在线播放高清视频| 精品久久久久久,| 久久久久久久精品吃奶| 99精品在免费线老司机午夜| 91av网一区二区| 久久精品91无色码中文字幕| 真人做人爱边吃奶动态| 热99在线观看视频| 熟女人妻精品中文字幕| 亚洲欧美日韩东京热| 男女下面进入的视频免费午夜| 日韩高清综合在线| 小蜜桃在线观看免费完整版高清| 一区福利在线观看| 久99久视频精品免费| 91在线观看av| 国产伦精品一区二区三区视频9 | 欧美成人免费av一区二区三区| 免费无遮挡裸体视频| 国产91精品成人一区二区三区| 禁无遮挡网站| 亚洲国产日韩欧美精品在线观看 | 国产成人精品久久二区二区91| 日韩欧美在线二视频| 国产极品精品免费视频能看的| 国产精品九九99| 亚洲一区高清亚洲精品| 国产99白浆流出| 国产精品日韩av在线免费观看| 日本黄色片子视频| 久久精品91蜜桃| 亚洲精品国产精品久久久不卡| 免费在线观看成人毛片| 免费av不卡在线播放| 五月玫瑰六月丁香| 欧美一级毛片孕妇| 一进一出抽搐动态| 欧美乱色亚洲激情| 97超视频在线观看视频| av片东京热男人的天堂| 午夜两性在线视频| 国产精品国产高清国产av| 欧美午夜高清在线| aaaaa片日本免费| 精品福利观看| 哪里可以看免费的av片| 国产69精品久久久久777片 | 亚洲成av人片免费观看| 中文资源天堂在线| 久久中文看片网| 亚洲成人久久爱视频| 中文资源天堂在线| 在线免费观看的www视频| av女优亚洲男人天堂 | 国产蜜桃级精品一区二区三区| 亚洲五月婷婷丁香| 精品一区二区三区av网在线观看| 色哟哟哟哟哟哟| 亚洲九九香蕉| xxx96com| av国产免费在线观看| 亚洲精华国产精华精| 欧美最黄视频在线播放免费| 少妇丰满av| 欧美色视频一区免费| 国产免费男女视频| a级毛片a级免费在线| 亚洲男人的天堂狠狠| 国产午夜福利久久久久久| 叶爱在线成人免费视频播放| 国产在线精品亚洲第一网站| 国产真人三级小视频在线观看| 男插女下体视频免费在线播放| 99久久99久久久精品蜜桃| 久久婷婷人人爽人人干人人爱| 成人av在线播放网站| 亚洲va日本ⅴa欧美va伊人久久| 国产91精品成人一区二区三区| 99riav亚洲国产免费| 日日摸夜夜添夜夜添小说| 此物有八面人人有两片| 中国美女看黄片| 午夜日韩欧美国产| 日本精品一区二区三区蜜桃| 精品99又大又爽又粗少妇毛片 | 中出人妻视频一区二区| 亚洲中文字幕一区二区三区有码在线看 | 亚洲国产精品成人综合色| 国产精品一区二区免费欧美| av国产免费在线观看| 麻豆国产av国片精品| 国内精品一区二区在线观看| 综合色av麻豆| 一卡2卡三卡四卡精品乱码亚洲| 亚洲中文字幕日韩| 91在线精品国自产拍蜜月 | 1000部很黄的大片| 精品久久久久久久末码| 成人永久免费在线观看视频| 国产三级黄色录像| 一级毛片女人18水好多| 男人舔女人下体高潮全视频| 国产97色在线日韩免费| 我的老师免费观看完整版| 51午夜福利影视在线观看| 日韩大尺度精品在线看网址| 九九久久精品国产亚洲av麻豆 | 悠悠久久av| 欧美zozozo另类| 日日摸夜夜添夜夜添小说| 精品国产三级普通话版| 国产单亲对白刺激| av天堂在线播放| 午夜福利视频1000在线观看| 国产精品98久久久久久宅男小说| 美女午夜性视频免费| 亚洲成人久久爱视频| 黑人巨大精品欧美一区二区mp4| 国内毛片毛片毛片毛片毛片| 免费一级毛片在线播放高清视频| 亚洲av美国av| av黄色大香蕉| 午夜a级毛片| 亚洲欧美激情综合另类| 悠悠久久av| 午夜福利在线观看吧| 亚洲一区二区三区不卡视频| 亚洲精品美女久久av网站| 精品久久久久久成人av| 91av网一区二区| 女警被强在线播放| 99久久久亚洲精品蜜臀av| 欧美成人一区二区免费高清观看 | 亚洲aⅴ乱码一区二区在线播放| 国产av麻豆久久久久久久| 日韩欧美一区二区三区在线观看| 一卡2卡三卡四卡精品乱码亚洲| 成人国产一区最新在线观看| 性色avwww在线观看| 色老头精品视频在线观看| or卡值多少钱| 精品久久久久久久人妻蜜臀av| 久久精品国产清高在天天线| 在线播放国产精品三级| 国产又色又爽无遮挡免费看| 久久精品国产99精品国产亚洲性色| 人妻夜夜爽99麻豆av| 99riav亚洲国产免费| 亚洲专区字幕在线| 国内毛片毛片毛片毛片毛片| 亚洲午夜理论影院| 国产精品日韩av在线免费观看| 亚洲人成网站高清观看| 色综合亚洲欧美另类图片| 亚洲第一欧美日韩一区二区三区| 免费看a级黄色片| 国产美女午夜福利| 人妻夜夜爽99麻豆av| 一边摸一边抽搐一进一小说| 手机成人av网站| 国产精品女同一区二区软件 | 日韩欧美免费精品| 搡老熟女国产l中国老女人| АⅤ资源中文在线天堂| 两个人看的免费小视频| 他把我摸到了高潮在线观看| svipshipincom国产片| 搡老岳熟女国产| 99热6这里只有精品| 国产99白浆流出| 国模一区二区三区四区视频 | 俺也久久电影网| 手机成人av网站| 国产黄色小视频在线观看| 免费一级毛片在线播放高清视频| 国产精品日韩av在线免费观看| 亚洲av成人精品一区久久| 国产又色又爽无遮挡免费看| 老司机福利观看| 午夜福利在线在线| 两人在一起打扑克的视频| 国产爱豆传媒在线观看| 亚洲色图 男人天堂 中文字幕| 国产精品一区二区精品视频观看| 久久久久久久精品吃奶| 欧美激情久久久久久爽电影| 亚洲av成人av| 99精品在免费线老司机午夜| 国产黄a三级三级三级人| 怎么达到女性高潮| 美女黄网站色视频| 制服人妻中文乱码| 黄片小视频在线播放| www.熟女人妻精品国产| 免费av毛片视频| www.精华液| 美女免费视频网站| 成人亚洲精品av一区二区| 午夜影院日韩av| 亚洲人成伊人成综合网2020| www.自偷自拍.com| 在线播放国产精品三级| 一进一出好大好爽视频| 国内久久婷婷六月综合欲色啪| 亚洲美女视频黄频| 欧美大码av| 久久天堂一区二区三区四区| 一本精品99久久精品77| 岛国视频午夜一区免费看| 欧美大码av| 国产精品久久久久久精品电影| 成熟少妇高潮喷水视频| 日韩 欧美 亚洲 中文字幕| 中文字幕精品亚洲无线码一区| 老汉色∧v一级毛片| a在线观看视频网站| 一个人看视频在线观看www免费 | 成人av一区二区三区在线看| ponron亚洲| 亚洲狠狠婷婷综合久久图片| 人妻丰满熟妇av一区二区三区| 美女午夜性视频免费| 1024手机看黄色片| 亚洲专区字幕在线| 亚洲国产欧洲综合997久久,| 美女 人体艺术 gogo| 免费看光身美女| 国产亚洲欧美在线一区二区| 中文字幕高清在线视频| 观看美女的网站| 欧美日本视频| 一夜夜www| 亚洲欧美日韩东京热| 91久久精品国产一区二区成人 | 黑人欧美特级aaaaaa片| 免费av毛片视频| 久久国产精品人妻蜜桃| 亚洲精品美女久久久久99蜜臀| 中文字幕久久专区| 色综合站精品国产| netflix在线观看网站| 国产av不卡久久| 国产高清三级在线| 波多野结衣高清作品| 精品国内亚洲2022精品成人| 亚洲国产欧美人成| 真人做人爱边吃奶动态| 亚洲九九香蕉| 老司机在亚洲福利影院| 国产精品久久视频播放| 一区二区三区激情视频| tocl精华| 日韩欧美国产一区二区入口| 欧美中文日本在线观看视频| 国产不卡一卡二| 久久草成人影院| 巨乳人妻的诱惑在线观看| 亚洲av电影不卡..在线观看| 欧美午夜高清在线| 国产激情欧美一区二区| 午夜福利免费观看在线| 韩国av一区二区三区四区| 黄色视频,在线免费观看| 国产熟女xx| 黄色视频,在线免费观看| 国产精品乱码一区二三区的特点| 日本 欧美在线| 国产精品亚洲一级av第二区| 久久久色成人| 国内少妇人妻偷人精品xxx网站 | 国产不卡一卡二| 99热这里只有精品一区 | 成年女人毛片免费观看观看9| 免费看a级黄色片| 色视频www国产| 午夜视频精品福利| 三级毛片av免费| 免费看a级黄色片| 少妇的丰满在线观看| 日本 av在线| 国内久久婷婷六月综合欲色啪| 丰满人妻一区二区三区视频av | 国产精品乱码一区二三区的特点| a级毛片在线看网站| 全区人妻精品视频| 可以在线观看的亚洲视频| 国产精品影院久久| 国产v大片淫在线免费观看| avwww免费| 日韩欧美在线二视频| 香蕉av资源在线| 精品国产乱子伦一区二区三区| 成年女人永久免费观看视频| 黑人操中国人逼视频| а√天堂www在线а√下载| 亚洲国产精品合色在线| 十八禁网站免费在线| 国产在线精品亚洲第一网站| 午夜福利免费观看在线| 国产精品野战在线观看| 欧美日韩乱码在线| 亚洲专区国产一区二区| 男人和女人高潮做爰伦理| 18禁观看日本| 国产亚洲欧美98| 18禁黄网站禁片免费观看直播| 国产精品永久免费网站| 国产av不卡久久| 看黄色毛片网站| 亚洲欧美一区二区三区黑人| 在线观看66精品国产| 精品一区二区三区视频在线 | 久久中文字幕人妻熟女| 国产乱人伦免费视频| 久久久色成人| 国产三级在线视频| 日本在线视频免费播放| 老汉色av国产亚洲站长工具| 亚洲无线观看免费| 又大又爽又粗| 伊人久久大香线蕉亚洲五| 狂野欧美激情性xxxx| 丰满的人妻完整版| 国产激情偷乱视频一区二区| 制服丝袜大香蕉在线| 日本成人三级电影网站| 久久久国产成人免费| 免费看a级黄色片| 久久伊人香网站| 老汉色av国产亚洲站长工具| 一本综合久久免费| 久久久久国内视频| 久久精品人妻少妇| 欧洲精品卡2卡3卡4卡5卡区| 亚洲熟妇中文字幕五十中出| 91在线观看av| 最近最新免费中文字幕在线| ponron亚洲| 成熟少妇高潮喷水视频| 中文字幕人妻丝袜一区二区| 国产午夜精品论理片| 精品久久久久久久末码| 18禁观看日本| 国产成人啪精品午夜网站| 亚洲精品美女久久久久99蜜臀| 亚洲午夜理论影院| 夜夜爽天天搞| 亚洲成人久久性| 又爽又黄无遮挡网站| 亚洲,欧美精品.| 久久久久久久久中文| 成人午夜高清在线视频| 亚洲 欧美 日韩 在线 免费| 免费在线观看成人毛片| 久久久成人免费电影| 亚洲中文日韩欧美视频| 久久久久久久精品吃奶| 999精品在线视频| 免费在线观看成人毛片| 亚洲美女视频黄频| 天堂av国产一区二区熟女人妻| 91久久精品国产一区二区成人 | 一二三四社区在线视频社区8| 精品久久久久久久人妻蜜臀av| 亚洲精品中文字幕一二三四区| 三级男女做爰猛烈吃奶摸视频| 亚洲国产色片| 国产成人系列免费观看| 国产精品一区二区免费欧美| 黄频高清免费视频| 欧美日韩综合久久久久久 | 亚洲一区二区三区色噜噜| 桃色一区二区三区在线观看| 一进一出好大好爽视频| 中文字幕高清在线视频| 色播亚洲综合网| ponron亚洲| 日韩欧美在线二视频| 亚洲精品国产精品久久久不卡| 亚洲最大成人中文| tocl精华| 国产成人系列免费观看| 午夜精品一区二区三区免费看| 少妇裸体淫交视频免费看高清| 99视频精品全部免费 在线 | 欧美黑人欧美精品刺激| 久久久久免费精品人妻一区二区| 男人的好看免费观看在线视频| 亚洲成人精品中文字幕电影| 制服人妻中文乱码| 国产精品精品国产色婷婷| 麻豆成人av在线观看| 桃色一区二区三区在线观看| 国产一区二区在线观看日韩 | 亚洲无线观看免费| 男插女下体视频免费在线播放| 男女午夜视频在线观看| 欧美大码av| 久久中文字幕一级| 色哟哟哟哟哟哟| 精品午夜福利视频在线观看一区| www.999成人在线观看| 欧美绝顶高潮抽搐喷水| 欧美乱妇无乱码| 一级a爱片免费观看的视频| 天堂网av新在线| 一夜夜www| 身体一侧抽搐| aaaaa片日本免费| 精品国产美女av久久久久小说| 9191精品国产免费久久| 日本一二三区视频观看| 久久国产乱子伦精品免费另类| 最近最新中文字幕大全电影3| 亚洲av电影在线进入| 久久香蕉精品热| 女人被狂操c到高潮| 亚洲av免费在线观看| 国产免费男女视频| 亚洲欧美日韩高清专用| 18禁黄网站禁片免费观看直播| 午夜精品一区二区三区免费看| 欧美激情在线99| 日韩免费av在线播放| 欧洲精品卡2卡3卡4卡5卡区| 亚洲av日韩精品久久久久久密| 国产精品一区二区精品视频观看| 亚洲人成伊人成综合网2020| 欧美黑人欧美精品刺激| 日韩精品中文字幕看吧| 国产淫片久久久久久久久 | 变态另类成人亚洲欧美熟女| 亚洲欧美激情综合另类| 国产精品一区二区三区四区久久| 少妇熟女aⅴ在线视频| 一本一本综合久久| 亚洲国产看品久久| 亚洲 欧美一区二区三区| 亚洲欧美日韩卡通动漫| 久久久久亚洲av毛片大全| 中文在线观看免费www的网站| 日韩欧美免费精品| av在线蜜桃| 日本三级黄在线观看| 我要搜黄色片| 12—13女人毛片做爰片一| 欧美中文综合在线视频| 噜噜噜噜噜久久久久久91| 全区人妻精品视频| 精品99又大又爽又粗少妇毛片 |