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

    虛擬波場變換方法在電磁法中的進展

    2020-04-08 13:26:54嚴家斌
    物探化探計算技術 2020年1期
    關鍵詞:方法

    林 超, 嚴家斌

    (1.中南大學 地球科學與信息物理學院,長沙 410083;2.有色金屬成礦預測與地質環(huán)境監(jiān)測教育部重點實驗室,長沙 410083)

    0 引言

    電磁法勘探是基于分析電磁波在大地介質中的傳播特性,達到研究地下地質體的賦存或構造特性的目的。通過使用天然或人工場源在大地中激發(fā)的交變電磁場,研究電磁場的時間和空間分布,分析觀測到的電磁響應信號,來獲得地下目標體電性分布的一種勘探方法[1]。在工程勘查[2-3]、資源勘探[4-5]、地下水探測[6-7]、地質填圖[8-9]和礦井災害預測[10-11]等方面都有廣泛地應用。電磁法方法的研究與應用是基于麥克斯韋方程組,為了避免電磁場特性的復雜分析以及降低各種干擾,提高電磁法方法的實用性,大部分電磁法方法都是采用準靜態(tài)近似的麥克斯韋方程組(如瞬變電磁法、海洋可控源電磁法等),但擴散場方程是不利于成像的,為了滿足高分辨率和精細成像以及快速解決復雜問題的要求,必須尋找新的數據處理方法,由此開展了虛擬波場變換研究。

    虛擬波場變換方法是指通過數學上的積分變換,將麥克斯韋擴散方程轉換為與電磁波類似的虛擬波動方程,在虛擬波場中進行計算和分析,該方法的初始研究是借鑒波動方程的性質,通過虛擬波場變換將地震中一些成熟的方法引入到電磁法的數據解釋中,提高電磁法數據解釋的分辨率,后來研究發(fā)現它在提高數值模擬的計算效率以及其他方面也可發(fā)揮重要作用。

    電磁法虛擬波場變換的方法研究開始于上世紀70年代。Kunetz[12]首先發(fā)現了擴散的電磁場和波動方程之間存在著聯系;Lavrent'ev[13]用數學方法證明并提出了擴散方程和波動方程之間的轉換關系;Lee等[14]通過在準靜態(tài)擴散的麥克斯韋方程組中引入“q函數”實現了擴散電磁場到虛擬波場的轉換,將虛擬波場變換方法運用到電磁法中;de Hoop[15]通過借助格林函數,在麥克斯韋方程組中應用拉普拉斯變換,推導出了擴散電磁場和虛擬波場的變換關系;Maa?[16]將最初的麥克斯韋方程通過數學變換,轉換到較低頻率依賴性的復頻率域中,提出一種基于復頻率域的虛擬波場方法;Mitte[17]在對麥克斯韋方程組進行變換時利用傅里葉變換,使得虛擬波場的變換過程更簡單。伴隨著虛擬波場變換理論方法的逐漸完善,虛擬波場的應用也越來越受到重視。為了充分認識波場變換技術和它的應用價值,這里對虛擬波場變換的方法原理以及應用做了簡要的分析與綜述。

    1 虛擬波場變換理論

    虛擬波場變換方法一般用在低頻電磁法中,從麥克斯韋基本的擴散方程出發(fā),對于準靜態(tài)情況下的時域擴散場,麥克斯韋方程組可以表示為:

    ▽×H(r,t)-σ(r)E(r,t)=-J(r,t)

    (1)

    ▽×E(r,t)+μ?tH(r,t)=-K(r,t)

    (2)

    式中:H(r,t)表示時間域磁場分量;E(r,t)表示時間域電場分量;J(r,t)表示外加電性源;K(r,t)表示外加磁性源;r表示空間變量;σ(r)表示電導率;μ表示磁導率。式(1)和式(2)經過數學物理變換最終可變?yōu)槿缦滦问剑?/p>

    (3)

    (4)

    (5)

    (6)

    式中:FE(r,g)為虛擬波場中的電場;FH(r,g)為虛擬波場中的磁場;SE(r,g)和SH(r,g)分別為虛擬波場中電場和磁場的源時間函數;c為虛擬波波速,與電導率、磁導率有關;g為虛擬波場中的時間變量與擴散場的時間對應;r為虛擬波場中的空間變量和擴散場中的一致。需要注意的是,從擴散場到虛擬波場的變換方法并不是唯一的,主要體現在引入的函數上,式(1)~式(6)是麥克斯韋擴散方程到虛擬波動方程的變換形式,通過積分變換可得虛擬波場恢復到時域擴散場的對應關系,其對應關系由如下方程表示:

    (7)

    (8)

    (9)

    (10)

    2 虛擬波場變換的方法

    虛擬波場變換中引入不同的函數就會有不同的變換形式以及波和擴散場之間的對應關系。目前,已知的波場變換方法有三種:Lee[14]提出的“q域法”方法、Maa?[16]提出的基于復頻率的方法及Mittet[17]提出的方法。

    2.1 q域法

    q域法是通過在時域擴散的麥克斯韋方程組中引入定義的“q函數”,實現了虛擬波場的變換,然后通過積分變換得到波場到擴散場的對應關系,大致過程如下:

    定義“q函數”:

    (11)

    式中:ω0為比例伸縮系數;t′為虛擬波場的時間。將“q函數”引入到時域的擴散麥克斯韋方程組(1)、方程組(2)中,即用q替代方程中的時間t,得到一組虛擬波場方程:

    (12)

    (13)

    虛擬波場與擴散場的對應關系為:

    (14)

    H″(r,q)=H′(r,q)

    (15)

    J″(r,q)=J′(r,q)

    (16)

    (17)

    將式(14)~式(17)代入到式(12)~式(13)之中可得重構的虛擬波動方程:

    -▽×H″(r,q)+σ(r)?qE″(r,q)=-J″(r,q)

    (18)

    ▽×E″(r,q)+μ?qH″(r,q)=-K″(r,q)

    (19)

    最后通過積分變換可從虛擬波場恢復到擴散場:

    (20)

    (21)

    (22)

    2.2 基于復頻率的方法

    復頻率域方法是通過在頻率域的麥克斯韋擴散方程組中引入定義的復頻率函數實現虛擬波場的變換,再通過積分變換從虛擬波場恢復到擴散場,其過程如下:

    復頻率函數定義為:

    iω=iω′(1-iω′a)

    (23)

    式中:a為常數,它的取值影響著時間步長的大小,一般根據實際情況選取合適的值,但不能取小于零的數;ω′為虛擬的角頻率;ω為擴散場的角頻率。將復頻率函數iω帶入到頻率域麥克斯韋擴散方程組中,其他量保持不變:

    -▽×H(r,ω)+σ(r)E(r,ω)=-J(r,ω)

    (24)

    ▽×E(r,ω)-iω′(1-iω′a)H(r,ω)=

    -K(r,ω)

    (25)

    在虛、實波場中對應關系為:

    E′(r,ω′)=E(r,ω)

    (26)

    H′(r,ω′)=(1-iω′a)H(r,ω)

    (27)

    J′(r,ω′)=(1-iω′a)J(r,ω)

    (28)

    K′(r,ω′)=K(r,ω)

    (29)

    利用對應關系可得重構的虛擬波動方程:

    -▽×H′(r,ω′)+(1-iω′a)σ(r)E′(r,ω)

    =-J′(r,ω′)

    (30)

    ▽×E′(r,ω′)-iω′μH′(r,ω) =

    -K′(r,ω′)

    (31)

    通過積分變換得到虛擬波場恢復到擴散場的對應關系為:

    (32)

    (33)

    (34)

    2.3 Mittet方法

    通過定義一個虛擬的介電常數,將虛擬介電常數代入到頻率域擴散麥克斯韋方程組中,然后直接在方程中乘以一個適當的函數,最后通過定義一個虛擬的角頻率以及對應關系實現虛擬波場的變換,該方法從虛擬波場恢復到擴散場需要利用格林函數,其過程如下:

    定義一個虛擬的介電常數:

    (35)

    式中:ω0為比例伸縮系數,它的取值也影響時間步長的大小,需要根據實際情況選取合適的值,且不能小于零。將虛擬的介電常數代入到頻率的擴散麥克斯韋方程組中可得:

    -▽×H(r,ω)+2ω0ε′(r)E(r,ω)=

    -J(r,ω)

    (36)

    ▽×E(r,ω)-iωμH(r,ω)=-K(r,ω)

    (37)

    (38)

    =-K(r,ω)

    (39)

    定義虛擬波動場中的角頻率以及對應關系:

    (40)

    E′(r,ω′)=E(r,ω)

    (41)

    (42)

    (43)

    K′(r,ω′)=K(r,ω)

    (44)

    式中:ω′為虛擬波場的角頻率,將式(40)~式(44)代入到式(38)或式(39)之中可得頻率域的虛擬波動方程,經過傅里葉變換可得時域的虛擬波動方程:

    -▽×H′(r,t′)+ε′(r)?t′E′(r,t′)=

    -J′(r,t′)

    (45)

    ▽×E′(r,t′)+?t′H′(r,t′)=-K′(r,t′)

    (46)

    從時域波場恢復到時域擴散場的過程如下:

    首先,根據類傅里葉變換積分將虛擬波場中的時域電磁場響應變換到虛擬波場中的頻率域電磁場響應:

    (47)

    (48)

    (49)

    其次,構造格林函數:

    (50)

    (51)

    最后,通過對所得格林函數進行傅里葉變換可得到真實擴散場的時域電磁響應:

    (52)

    (53)

    3 虛擬波場的應用

    虛擬波場方法可以有效地提高數值模擬的計算效率,目前主要針對三維的數值模擬和擬地震反演成像等數據處理與解釋方面。

    3.1 在提高數值模擬計算效率方面的應用

    虛擬波場變換方法能提高數值模擬的計算效率主要是基于它滿足的波動方程(式(5)~式(6))的性質,在擴散電磁場中,由于電磁波在地下介質中的傳播的速度很大,利用有限差分進行數值計算時,需要很小的時間步長才能滿足穩(wěn)定性條件。在虛擬波場中滿足的方程是波動方程,虛擬波的波速不再是現實中電磁波波速,可以選擇較小的虛擬波波速來獲得較大的時間步長且滿足它的穩(wěn)定性條件,這使得在虛擬波場中有限差分法的計算效率很高,最后再恢復到擴散場。積分方程和有限元等數值模擬方法也可以求解,只是顯示有限差分法具有占用內存小,建模簡單等優(yōu)勢,所以在該方法中比較適用。

    Maa?[16]提出了一種基于復頻率的波場變換方法之后,在虛擬波場中利用有限差分方法求解的巨大時間效益才體現出來。St?ren 等[18]在利用局部梯度優(yōu)化方法對海洋可控源數據進行三維反演時,在三維正演建模過程中采用Maa?[16]提出的虛擬波場變換方法,加快了失配梯度的計算速度。Mittet[17]提出了一種相對更簡單的虛擬波場變換方法,其主要是在Maa?提出的基于復頻率的波場變換方法上進行改進,并采用傅里葉變換去代替Lee等[14]的拉普拉斯變換。Mittet[22]將該方法引入到海洋大地三維數據反演的正演建模之中,同時,在正演計算過程中提出了一種非均勻網格垂直節(jié)點間距設計的標準以及采用了卷積完全匹配層進行邊界條件的處理,減少了計算區(qū)域的尺寸,提高了計算效率。Imamura等[19-20]在對海洋可控源電磁法三維全波形反演研究時,在正演建模時采用了Mittet[17]提出的虛擬波場變換方法,提高了正演的計算效率。Liu等[21]將Mittet[17]提出的虛擬波場變換方法用于瞬變電磁實測數據的二維正演,取得了良好效果。Hu等[23]將該方法引用到三維瞬變電磁的數值模擬,在正演計算中利用卷積完全匹配層處理邊界條件,并且在計算時將空氣層作為有限高阻引入到計算區(qū)域內,實現了對帶地形模型的模擬,最后將所得結果和SLDMA方法的對比,兩者具有良好的一致性,證明了該方法的有效性和計算的高效性。

    3.2 虛擬波場擬地震反演成像的應用

    虛擬波場擬地震反演成像主要是利用虛擬波場中與電磁波波傳播的有關特性,它的實質是通過數學上的積分變換提取電磁波中與波傳播有關的信息,使得虛擬波場類似于地震子波具有諸如反射、折射、衍射和透射等特征。虛擬波場擬地震成像可簡單的概述為:根據擴散場和虛擬波場的對應關系,利用奇異值分解等優(yōu)化算法將擴散電磁場信號轉換成虛擬波信號,然后借助于偏移成像、層析成像等地震方法以及一些其他的方法來對獲得的虛擬波信號進行成像等數據處理。

    虛擬波場擬地震成像的應用最早開始于20世紀90年代,Lee等[24]提出利用奇異值分解來進行波場逆變換,從獲得的穩(wěn)定波形中估計源到接收器的旅行時間,再用獲得的旅行時間進行層析成像來反映地下電阻率結構,并通過費馬原理的兩點光線追蹤算法來進行高效迭代。Lee等[25-26]提出一種直接從頻率擴散場垂直磁場中提取旅行時間的方法,在虛擬波場中利用提取的旅行時間進行非線性層析成像,成功地估算了合成電磁數據的電阻率,并在實測數據中應用,證明了該方法的有效性。Kusuda等[27]在研究如何提高海洋可控源電磁法對甲烷水合物檢測時,利用奇異值分解的方法,將海洋可控源數據轉換成虛擬波場數據,分離出甲烷水合物異常區(qū)的波形。利用分離出的波形進行成像分析,提高了海洋可控源電磁法對甲烷水合物的分辨率。Amani等[28]提出利用地震中的逆時偏移成像和克希霍夫偏移成像的方法對虛擬波場變換得到的穩(wěn)定波形成像,并將這兩種方法用于海洋可控源電磁法數據的二維反演成像,通過綜合比較得出逆時偏移成像的精度較高。

    國內的瞬變電磁擬地震成像開始于20世紀90年代,研究方向可分為:波場優(yōu)化算法、反褶積壓縮子波寬度、三維曲面延拓成像、合成孔徑成像。波場優(yōu)化算法主要是用于在波場逆變換中,減少不適定問題帶來的高度欠定的影響,從而使波場逆變換獲得的波形穩(wěn)定,在這方面的研究可參考文獻[29-32];反褶積壓縮子波寬度是削弱波場逆變換中核函數隨時間的增加分布范圍增大帶來的波形較寬的問題,其實是對波場逆變換中獲得的波形進行濾波處理,在這方面的研究可參考文獻[33-34];三維曲面延拓是為了消除地形起伏帶來的曲面以及測線不規(guī)則對反演成像帶來的影響,提高反演成像的分辨率。在這方面的研究可參考文獻[35-37];合成孔徑處理是借鑒合成孔徑雷達成像的方法,即通過利用機載真實孔徑發(fā)射線圈與目標的相對運動,把尺寸較小的真實天線孔徑用數據處理的方法合成一較大孔徑的發(fā)射線圈,提高它的分辨能力和穿透力,由于雷達滿足的偏微分方程和虛擬波場滿足的偏微分方程都是波動方程,所以該方法可用于虛擬波場逆變換反演成像中,在這方面的研究可參考文獻[38-40]。

    3.3 其他方面的應用

    在利用有限差分求解方程時,對函數的一階偏導數的離散一般用含算子系數的中心差分近似來表示, Pasalic等[41]研究發(fā)現,可以通過對算子系數進行最優(yōu)化然后取高階的差分算子來減少離散網格的尺寸,并將該方法應用到虛擬波動方程的離散中,進一步提高了虛擬波場-有限差分法的計算效率。Imamura等[42]在虛擬波動方程離散上采用了粒子法,并提出了交錯粒子法和非交錯粒子法,通過模擬對比,交錯粒子在數值上比非交錯粒子有效率更高,主要是在同等的空間間隔中,交錯粒子只需要一半的數目就能達到和非交錯一樣的精度。Mittet[43]在對虛擬波域的波場和地震波波場關系研究過程中,為了使兩者條件一致,利用聲波來近似代替地震波,通過對虛擬波場中的波形和聲波波形的分析,證明了地震波場和海洋可控源電磁波場之間的相似性,以此來了解從虛擬波域到擴散域的轉換效果。

    4 虛擬波場存在的問題

    虛擬波場變換是通過數學上的變換實現的,從式(7)~式(10)可知,它有兩個過程:①波場到擴散場的正變換;②擴散場到虛擬波場的逆變換。虛擬波場方法的主要應用思路基本上都是圍繞這兩個不同的過程。但是這兩個不同的過程各有不足,在逆變換中,獲得的虛擬波類似于地震子波一樣具有傳播、反射、透射的特征,雖然在理論模型和一些實際應用中驗證確實可行,實際上它是存在著一些問題的,例如:逆變換是第一類Fredholm型算子方程的反問題,它本身是非常復雜的,現在一般的解法都是求線性逼近系統(tǒng)的最小二乘解,但是這種求解方法會帶來一定的誤差,該誤差很難估計;兩種波場的物理背景不同,一個是客觀存在的波場,一個是由數學變換得到的虛擬波場,將地震方法用于虛擬的波場中是否能反映出地下真實情況;擴散場具有頻散和吸收的特性,十分復雜,用簡單的虛擬波場來代替是否符合實際情況,以及擴散場的波速跟頻率有關,而虛擬波場中的波速跟頻率無關,這種不一致帶來的影響等。在正變換中,雖然最終恢復到擴散場,但是在恢復過程中,積分變換帶來的誤差也很難估計,以及可能會引入高頻成分、噪聲等干擾因素。

    5 結論

    筆者詳細介紹了虛擬波場變換的基本原理和方法,評價了它的優(yōu)缺點,并對其在地球物理中的應用做了綜述。虛擬波場在現實中是不存在的,在應用中,從擴散場到虛擬波場的變換相當于一個提取電磁波中有用信息的過程,虛擬波場相當于一個包含必要信息的容器。該方法為基于麥克斯韋擴散方程的電磁法方法的數據解釋提供了高效的數值模擬和新的思路,有助于提高數據解釋的效率和分辨率,豐富了電磁法數據解釋的方法。

    致謝:

    感謝評審專家們提出的中肯的有建設性的修改意見,對本文的改進有很大的幫助。

    猜你喜歡
    方法
    中醫(yī)特有的急救方法
    中老年保健(2021年9期)2021-08-24 03:52:04
    高中數學教學改革的方法
    河北畫報(2021年2期)2021-05-25 02:07:46
    化學反應多變幻 “虛擬”方法幫大忙
    變快的方法
    兒童繪本(2020年5期)2020-04-07 17:46:30
    學習方法
    可能是方法不對
    用對方法才能瘦
    Coco薇(2016年2期)2016-03-22 02:42:52
    最有效的簡單方法
    山東青年(2016年1期)2016-02-28 14:25:23
    四大方法 教你不再“坐以待病”!
    Coco薇(2015年1期)2015-08-13 02:47:34
    賺錢方法
    综合色av麻豆| 一边摸一边抽搐一进一小说| 国产v大片淫在线免费观看| 精品国产乱子伦一区二区三区| 五月伊人婷婷丁香| 夜夜夜夜夜久久久久| 亚洲真实伦在线观看| 99久久国产精品久久久| 欧美+亚洲+日韩+国产| 亚洲人与动物交配视频| 午夜福利高清视频| 久9热在线精品视频| 国产爱豆传媒在线观看| 香蕉国产在线看| 国内精品美女久久久久久| 麻豆国产av国片精品| 午夜视频精品福利| 十八禁人妻一区二区| 精品国产超薄肉色丝袜足j| 日韩欧美国产在线观看| 91麻豆av在线| 无遮挡黄片免费观看| 亚洲av电影不卡..在线观看| 日本 av在线| 国产高清videossex| 亚洲熟妇中文字幕五十中出| 免费在线观看日本一区| 免费高清视频大片| 国产伦精品一区二区三区四那| 色综合站精品国产| 国产亚洲精品综合一区在线观看| 制服丝袜大香蕉在线| 手机成人av网站| 久久精品91无色码中文字幕| 欧美日韩黄片免| 亚洲人成伊人成综合网2020| 日本五十路高清| 老司机午夜福利在线观看视频| 国产伦精品一区二区三区视频9 | 午夜两性在线视频| 91久久精品国产一区二区成人 | 女生性感内裤真人,穿戴方法视频| 波多野结衣巨乳人妻| 欧美日韩乱码在线| 一个人看的www免费观看视频| 嫁个100分男人电影在线观看| 亚洲中文字幕一区二区三区有码在线看 | 此物有八面人人有两片| 18禁观看日本| 午夜激情福利司机影院| 国产av在哪里看| 999精品在线视频| 亚洲国产高清在线一区二区三| 久久国产精品影院| 国产视频一区二区在线看| 免费av不卡在线播放| 好看av亚洲va欧美ⅴa在| 香蕉久久夜色| 高潮久久久久久久久久久不卡| 啪啪无遮挡十八禁网站| 欧美黄色片欧美黄色片| ponron亚洲| 久久久久久人人人人人| 午夜福利高清视频| 人人妻人人澡欧美一区二区| 国产精品久久久久久久电影 | 日韩欧美三级三区| 最新中文字幕久久久久 | 国产三级在线视频| 午夜福利在线在线| 日韩欧美 国产精品| 欧美+亚洲+日韩+国产| 一级a爱片免费观看的视频| 久久久久久久午夜电影| 可以在线观看的亚洲视频| 中文字幕最新亚洲高清| 亚洲电影在线观看av| 国产高清激情床上av| 男女之事视频高清在线观看| 中出人妻视频一区二区| 老司机深夜福利视频在线观看| 亚洲 国产 在线| 岛国视频午夜一区免费看| 男女视频在线观看网站免费| 亚洲国产中文字幕在线视频| 在线观看舔阴道视频| 2021天堂中文幕一二区在线观| 国产黄a三级三级三级人| 一本久久中文字幕| 久久人妻av系列| 黄色日韩在线| 18禁国产床啪视频网站| 最近在线观看免费完整版| 久久热在线av| 身体一侧抽搐| 成年女人永久免费观看视频| 午夜精品久久久久久毛片777| 啦啦啦观看免费观看视频高清| 黑人巨大精品欧美一区二区mp4| 婷婷六月久久综合丁香| 久久精品91无色码中文字幕| www.精华液| 亚洲国产精品合色在线| 巨乳人妻的诱惑在线观看| cao死你这个sao货| bbb黄色大片| 欧美中文综合在线视频| 怎么达到女性高潮| 性欧美人与动物交配| 又大又爽又粗| 非洲黑人性xxxx精品又粗又长| 丝袜人妻中文字幕| 美女 人体艺术 gogo| 一个人观看的视频www高清免费观看 | 极品教师在线免费播放| 我要搜黄色片| 国产精品女同一区二区软件 | 亚洲 国产 在线| 午夜影院日韩av| 久久精品国产清高在天天线| 亚洲精品国产精品久久久不卡| 啪啪无遮挡十八禁网站| 精品欧美国产一区二区三| 精品日产1卡2卡| 91九色精品人成在线观看| 精品一区二区三区视频在线观看免费| 精品人妻1区二区| 1024香蕉在线观看| 狠狠狠狠99中文字幕| 麻豆成人午夜福利视频| 国产日本99.免费观看| 国产久久久一区二区三区| 窝窝影院91人妻| 特级一级黄色大片| 欧美日韩黄片免| 亚洲国产精品sss在线观看| 国产成人福利小说| av视频在线观看入口| 真人一进一出gif抽搐免费| 人妻久久中文字幕网| 人妻久久中文字幕网| 欧美色欧美亚洲另类二区| 国产97色在线日韩免费| 亚洲人成电影免费在线| 久久九九热精品免费| 最近最新免费中文字幕在线| 亚洲av熟女| av视频在线观看入口| 亚洲欧美日韩卡通动漫| 色综合站精品国产| 国产成人福利小说| 无遮挡黄片免费观看| 国产精品久久久久久亚洲av鲁大| 他把我摸到了高潮在线观看| 91av网站免费观看| 啪啪无遮挡十八禁网站| 黄色片一级片一级黄色片| 国产成人啪精品午夜网站| 人妻丰满熟妇av一区二区三区| 丁香六月欧美| 日日夜夜操网爽| 亚洲av片天天在线观看| 91在线观看av| 国产成+人综合+亚洲专区| 好看av亚洲va欧美ⅴa在| 无限看片的www在线观看| 1000部很黄的大片| 国产又色又爽无遮挡免费看| 午夜福利免费观看在线| 十八禁人妻一区二区| 欧美丝袜亚洲另类 | 亚洲天堂国产精品一区在线| 国产精品女同一区二区软件 | 丰满人妻一区二区三区视频av | 国产一区二区在线观看日韩 | 欧美成人一区二区免费高清观看 | 在线十欧美十亚洲十日本专区| 免费无遮挡裸体视频| 搡老熟女国产l中国老女人| 悠悠久久av| 香蕉av资源在线| 亚洲无线在线观看| 99精品久久久久人妻精品| 国产高清videossex| 欧美不卡视频在线免费观看| 国产av在哪里看| 色综合站精品国产| 国产亚洲精品久久久com| 黄色片一级片一级黄色片| 亚洲av电影不卡..在线观看| 亚洲av电影不卡..在线观看| 国产在线精品亚洲第一网站| 波多野结衣巨乳人妻| 久久久久免费精品人妻一区二区| 一个人看视频在线观看www免费 | 亚洲在线观看片| 欧美日本视频| 亚洲欧美日韩东京热| 一级黄色大片毛片| 成人特级黄色片久久久久久久| 欧美+亚洲+日韩+国产| 日韩成人在线观看一区二区三区| 精品免费久久久久久久清纯| 成人18禁在线播放| 亚洲精品在线美女| 国产精品自产拍在线观看55亚洲| 国产亚洲欧美98| 精品久久久久久久人妻蜜臀av| 日韩免费av在线播放| 综合色av麻豆| 日本在线视频免费播放| 变态另类成人亚洲欧美熟女| 伊人久久大香线蕉亚洲五| 国产精品九九99| 中文字幕精品亚洲无线码一区| 欧美乱码精品一区二区三区| 中文字幕人成人乱码亚洲影| 禁无遮挡网站| 亚洲一区二区三区不卡视频| 少妇人妻一区二区三区视频| 十八禁人妻一区二区| 日日夜夜操网爽| 日本与韩国留学比较| cao死你这个sao货| 亚洲人成网站高清观看| 在线视频色国产色| tocl精华| 国产熟女xx| 一a级毛片在线观看| 日本一二三区视频观看| 可以在线观看毛片的网站| 悠悠久久av| 香蕉av资源在线| 男女之事视频高清在线观看| 精品久久久久久久末码| 给我免费播放毛片高清在线观看| 国内精品久久久久久久电影| 99久久成人亚洲精品观看| 一本综合久久免费| 老司机在亚洲福利影院| 999久久久国产精品视频| 身体一侧抽搐| 午夜影院日韩av| 国产一区二区在线av高清观看| 午夜亚洲福利在线播放| 欧美色视频一区免费| 亚洲avbb在线观看| 黑人操中国人逼视频| 午夜影院日韩av| 天天躁日日操中文字幕| 免费在线观看日本一区| 少妇人妻一区二区三区视频| 18美女黄网站色大片免费观看| 黑人欧美特级aaaaaa片| 亚洲精品在线观看二区| 亚洲无线在线观看| 桃色一区二区三区在线观看| 亚洲国产欧美网| 久久久国产欧美日韩av| 日韩欧美免费精品| 美女高潮的动态| 亚洲专区国产一区二区| 国产一区二区激情短视频| 国产一区在线观看成人免费| 免费高清视频大片| 精品熟女少妇八av免费久了| 国产精品一区二区精品视频观看| 午夜福利18| 亚洲精品久久国产高清桃花| 高潮久久久久久久久久久不卡| 国内少妇人妻偷人精品xxx网站 | 一个人看视频在线观看www免费 | 亚洲av成人av| 18禁黄网站禁片午夜丰满| 色老头精品视频在线观看| 国产免费av片在线观看野外av| 久久精品91无色码中文字幕| 看片在线看免费视频| 久久人人精品亚洲av| 最近最新免费中文字幕在线| 国产亚洲精品av在线| 国产免费av片在线观看野外av| 99精品久久久久人妻精品| 99国产精品一区二区三区| 久久久久久久久久黄片| 久久这里只有精品19| 身体一侧抽搐| 在线观看午夜福利视频| av福利片在线观看| 日日干狠狠操夜夜爽| 国产亚洲精品综合一区在线观看| 久久久国产精品麻豆| 久久久久久久久免费视频了| 久久99热这里只有精品18| 国产单亲对白刺激| 亚洲国产精品999在线| 熟女人妻精品中文字幕| 12—13女人毛片做爰片一| 国产成人影院久久av| 在线十欧美十亚洲十日本专区| 嫩草影院精品99| 窝窝影院91人妻| 99热只有精品国产| 最近最新中文字幕大全电影3| 亚洲国产看品久久| 后天国语完整版免费观看| 又紧又爽又黄一区二区| 亚洲av成人av| 女人被狂操c到高潮| 日韩中文字幕欧美一区二区| 中文字幕人妻丝袜一区二区| 禁无遮挡网站| 夜夜躁狠狠躁天天躁| 亚洲黑人精品在线| 两个人视频免费观看高清| 国产伦人伦偷精品视频| 精品99又大又爽又粗少妇毛片 | 国产黄a三级三级三级人| 91av网一区二区| 99在线人妻在线中文字幕| 精品欧美国产一区二区三| 9191精品国产免费久久| 午夜福利在线在线| 叶爱在线成人免费视频播放| 性欧美人与动物交配| 国产亚洲精品一区二区www| av中文乱码字幕在线| 国产精品一区二区精品视频观看| 高清在线国产一区| 日本黄大片高清| 国产av不卡久久| 久9热在线精品视频| 啦啦啦韩国在线观看视频| 在线观看免费视频日本深夜| 免费看a级黄色片| 久久精品aⅴ一区二区三区四区| 亚洲性夜色夜夜综合| 亚洲色图 男人天堂 中文字幕| 精品久久久久久久久久免费视频| 亚洲自偷自拍图片 自拍| 小说图片视频综合网站| 国产一区在线观看成人免费| 一夜夜www| 无人区码免费观看不卡| 两个人看的免费小视频| 午夜两性在线视频| 亚洲,欧美精品.| 国产欧美日韩一区二区精品| 免费观看精品视频网站| 99国产极品粉嫩在线观看| 午夜两性在线视频| 18禁黄网站禁片免费观看直播| av欧美777| 国产成人精品久久二区二区91| 成人鲁丝片一二三区免费| 国内精品一区二区在线观看| 国内揄拍国产精品人妻在线| 亚洲欧美日韩卡通动漫| 一个人观看的视频www高清免费观看 | 日韩欧美 国产精品| 国产午夜精品久久久久久| 欧美日韩一级在线毛片| 高潮久久久久久久久久久不卡| 亚洲专区国产一区二区| 国产欧美日韩一区二区三| 一边摸一边抽搐一进一小说| 欧美中文日本在线观看视频| 亚洲午夜理论影院| 欧美日本亚洲视频在线播放| 免费大片18禁| 中文字幕久久专区| 国产精品,欧美在线| 午夜福利成人在线免费观看| 午夜福利18| 午夜免费激情av| 国内精品一区二区在线观看| 叶爱在线成人免费视频播放| 精品国产乱子伦一区二区三区| 搡老岳熟女国产| 嫩草影院入口| 精品久久久久久久人妻蜜臀av| 男人舔女人下体高潮全视频| 亚洲欧美日韩高清专用| 国产毛片a区久久久久| 90打野战视频偷拍视频| 啦啦啦韩国在线观看视频| 国产91精品成人一区二区三区| 99久国产av精品| av天堂在线播放| 国产欧美日韩一区二区精品| 一二三四在线观看免费中文在| 制服丝袜大香蕉在线| 亚洲 欧美一区二区三区| 久久国产乱子伦精品免费另类| 亚洲成人久久爱视频| 亚洲人成网站在线播放欧美日韩| 亚洲欧美日韩卡通动漫| 成人永久免费在线观看视频| АⅤ资源中文在线天堂| 精品久久久久久久久久免费视频| 日韩欧美一区二区三区在线观看| 成人三级做爰电影| av国产免费在线观看| xxxwww97欧美| 伊人久久大香线蕉亚洲五| 后天国语完整版免费观看| 日韩中文字幕欧美一区二区| 中国美女看黄片| 成人国产综合亚洲| 91久久精品国产一区二区成人 | 亚洲第一电影网av| 欧美黑人欧美精品刺激| 色综合婷婷激情| 欧美日韩乱码在线| 国产黄片美女视频| 一本一本综合久久| 好男人电影高清在线观看| 桃色一区二区三区在线观看| avwww免费| 18禁裸乳无遮挡免费网站照片| 国语自产精品视频在线第100页| 日韩高清综合在线| 午夜视频精品福利| 国产精品亚洲av一区麻豆| avwww免费| 国产又黄又爽又无遮挡在线| 亚洲中文av在线| 亚洲av成人一区二区三| 亚洲av熟女| 国产亚洲av高清不卡| 欧美成狂野欧美在线观看| 91字幕亚洲| 国产97色在线日韩免费| 男女做爰动态图高潮gif福利片| 最近最新中文字幕大全电影3| 亚洲av成人av| 亚洲av成人一区二区三| 九色成人免费人妻av| 久久人妻av系列| 国内少妇人妻偷人精品xxx网站 | 麻豆av在线久日| 久久精品国产综合久久久| 亚洲激情在线av| 欧美黄色淫秽网站| 精品一区二区三区视频在线观看免费| 亚洲欧美日韩卡通动漫| 国产激情欧美一区二区| 亚洲激情在线av| 校园春色视频在线观看| 伊人久久大香线蕉亚洲五| 最近最新免费中文字幕在线| 国产伦在线观看视频一区| www.999成人在线观看| 欧美日韩亚洲国产一区二区在线观看| 成人国产一区最新在线观看| 欧美乱妇无乱码| 麻豆成人av在线观看| 天天添夜夜摸| 国产精品野战在线观看| 色吧在线观看| 欧洲精品卡2卡3卡4卡5卡区| 在线视频色国产色| 黄色片一级片一级黄色片| 97超级碰碰碰精品色视频在线观看| 国产毛片a区久久久久| 欧美性猛交╳xxx乱大交人| 精品久久久久久久人妻蜜臀av| 哪里可以看免费的av片| 九色国产91popny在线| 香蕉av资源在线| 老司机午夜福利在线观看视频| 波多野结衣高清无吗| 最近在线观看免费完整版| 亚洲成人免费电影在线观看| 亚洲五月天丁香| 又黄又爽又免费观看的视频| 精品国产乱码久久久久久男人| 波多野结衣巨乳人妻| 天堂√8在线中文| 黑人欧美特级aaaaaa片| 日韩欧美三级三区| 男女那种视频在线观看| 91麻豆精品激情在线观看国产| 国产av一区在线观看免费| 国产精品一及| 国产人伦9x9x在线观看| 热99在线观看视频| 法律面前人人平等表现在哪些方面| 搡老岳熟女国产| 亚洲成a人片在线一区二区| 国产爱豆传媒在线观看| 久久午夜亚洲精品久久| 亚洲 欧美 日韩 在线 免费| 俺也久久电影网| 女同久久另类99精品国产91| 老鸭窝网址在线观看| 色综合婷婷激情| 久久久国产精品麻豆| 国产精品野战在线观看| 哪里可以看免费的av片| 看免费av毛片| 日韩成人在线观看一区二区三区| 国产淫片久久久久久久久 | 日本 av在线| 狂野欧美白嫩少妇大欣赏| 在线观看免费视频日本深夜| 中国美女看黄片| 网址你懂的国产日韩在线| 观看美女的网站| 欧美绝顶高潮抽搐喷水| 观看美女的网站| 亚洲自偷自拍图片 自拍| 午夜亚洲福利在线播放| 亚洲欧美日韩东京热| 色播亚洲综合网| 亚洲欧美精品综合一区二区三区| 色播亚洲综合网| 国产黄片美女视频| 成人无遮挡网站| 全区人妻精品视频| 禁无遮挡网站| av在线蜜桃| 国产美女午夜福利| 欧美又色又爽又黄视频| 首页视频小说图片口味搜索| 一级a爱片免费观看的视频| 色综合欧美亚洲国产小说| 日韩 欧美 亚洲 中文字幕| 蜜桃久久精品国产亚洲av| 他把我摸到了高潮在线观看| 九色国产91popny在线| 国产av不卡久久| 9191精品国产免费久久| 久久久国产精品麻豆| 欧美乱码精品一区二区三区| 叶爱在线成人免费视频播放| 一本综合久久免费| 18禁美女被吸乳视频| 欧美激情久久久久久爽电影| 免费搜索国产男女视频| 国产精品国产高清国产av| 亚洲性夜色夜夜综合| 国产伦精品一区二区三区四那| 亚洲激情在线av| 日韩欧美在线二视频| 丁香欧美五月| 色视频www国产| 欧美又色又爽又黄视频| 国产日本99.免费观看| 欧美又色又爽又黄视频| 美女 人体艺术 gogo| 91字幕亚洲| 99精品久久久久人妻精品| 精品国产亚洲在线| 又大又爽又粗| 欧美日韩综合久久久久久 | 亚洲精品色激情综合| 国产精品 欧美亚洲| 午夜免费观看网址| 欧美国产日韩亚洲一区| 一级a爱片免费观看的视频| 久久久国产成人精品二区| 国产三级在线视频| 很黄的视频免费| 国产高清视频在线播放一区| 18禁裸乳无遮挡免费网站照片| 小说图片视频综合网站| 亚洲精品久久国产高清桃花| 黄色成人免费大全| 日韩av在线大香蕉| 日韩国内少妇激情av| 性色av乱码一区二区三区2| 亚洲一区二区三区不卡视频| 成人午夜高清在线视频| 老司机在亚洲福利影院| 久久亚洲真实| 动漫黄色视频在线观看| 国产真实乱freesex| 最好的美女福利视频网| 校园春色视频在线观看| 国产精品,欧美在线| 国产激情欧美一区二区| 偷拍熟女少妇极品色| 99国产精品一区二区三区| 最近最新免费中文字幕在线| 老汉色∧v一级毛片| 在线免费观看不下载黄p国产 | 国内毛片毛片毛片毛片毛片| 老熟妇乱子伦视频在线观看| 国产亚洲av嫩草精品影院| www.自偷自拍.com| 成年免费大片在线观看| 国产亚洲精品久久久com| 午夜福利在线在线| 亚洲国产精品久久男人天堂| 国产av不卡久久| 国产亚洲精品av在线| 亚洲aⅴ乱码一区二区在线播放| 精品国产露脸久久av麻豆 | 免费在线观看成人毛片| 日本色播在线视频| 久久人妻av系列| 少妇的逼好多水| 91狼人影院| 欧美zozozo另类| 老司机影院毛片| 国产亚洲精品久久久com| 91久久精品电影网| 人妻系列 视频| 男人舔女人下体高潮全视频| 国产免费又黄又爽又色| 99九九线精品视频在线观看视频| 国产精品.久久久| 久久久久网色| 少妇人妻一区二区三区视频|