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

    基于Python語言的瑞利面波頻散反演實行方案

    2022-06-16 10:24:18吳為治婁利王鵬王賓
    物探與化探 2022年3期
    關(guān)鍵詞:面波瑞利橫波

    吳為治,婁利,王鵬,王賓

    (安徽省勘查技術(shù)院,安徽 合肥 230041)

    0 引言

    Python語言作為一種解釋型的高級語言,其語法簡單、代碼可讀性較高且易擴展移植,近年來迅速發(fā)展成為科學、工程研究領(lǐng)域內(nèi)流行的編程語言。Python語言還具有功能強大且開源標準函數(shù)庫以及第三方函數(shù)庫,眾多科學計算軟件都提供有Python的安裝調(diào)用接口,Github網(wǎng)站也提供Python軟件包的開源以方便開發(fā)者搭建展開研究工作的平臺。正是這樣,調(diào)用高效完善的函數(shù)庫,既能滿足專業(yè)化研究的需求,又能進一步提高開發(fā)效率。彭一波等[1]提出利用Python 語言地震學數(shù)據(jù)處理軟件包obspy進行背景噪聲研究的工作方案,實現(xiàn)了兩個地震臺站之間噪聲互相關(guān)格林函數(shù)波頻散曲線的提取,并進行背景噪聲源的時空性研究。這是利用Python語言建立地球物理研究平臺一次成功的案例,然而其并未進行進一步頻散反演。我們知道面波頻散曲線含有豐富的地下層狀結(jié)構(gòu)信息,尤其對橫波速度結(jié)構(gòu)特別敏感,所以通過對頻散曲線的反演可以得到地下層狀介質(zhì)的物理性質(zhì)[2]。因此基于Python語言的頻散反演方法的研究很有必要,可以為今后繼續(xù)研究背景噪聲成像、面波勘探提供開源的支持。

    1 研究基礎(chǔ)

    1.1 頻散正演算法

    頻散曲線的反演基于正演算法,面波頻散的正演計算方法眾多,如Knopoff算法[3],廣義反射—透射系數(shù)法[4]等,均能快速準確地求取面波頻散。Python語言第三方函數(shù)庫提供有disba軟件包(https://pypi.org/project/disba)用于計算頻散,該軟件包使用Dunkin矩陣或者快速增量矩陣算法計算瑞利波相和群速度頻散曲線;使用Thomson-Haskell方法計算勒夫波相和群速度頻散曲線。GitHub上則有提供pysurf96軟件包(https://github.com/miili/pysurf96)用來計算頻散,其算法基于圣路易斯大學的CPS(computer programs in seismology)中surf96程序[5],雖然需要Fortran編譯環(huán)境,但穩(wěn)定性更好。這里我們選擇用pysurf96程序來測試(附錄示例1),選取一個簡單的兩層模型(表1)來正演瑞利面波的頻散曲線(圖1)。

    表1 兩層層狀介質(zhì)模型參數(shù)

    圖1 基于pysurf96計算的兩層層狀介質(zhì)模型的瑞利面波頻散曲線Fig.1 Rayleigh wave dispersion curve of two layered media model based on pysurf96

    從圖1中可以看出基階瑞利面波相速度在低頻部分接近模型半空間橫波速度的0.9倍左右,高頻部分接近模型上層橫波速度的0.9倍左右;二階面波相速度在高頻部分趨近模型半空間橫波速度;三階面波相速度則在低頻接近模型上層的橫波速度。這與夏江海在《高頻面波方法》中提到的現(xiàn)象一致[2]。同時,我們在測試不同速度模型中發(fā)現(xiàn)當?shù)貙幽P椭邪肟臻g的波速過低時,會無法得到某些頻率的瑞利面波相速度,如使用軟件包中surf96函數(shù)計算頻散會報錯,使用surfdisp96函數(shù)則需要先定義返回值為復數(shù)。這是因為頻散方程在這些頻率的根是復數(shù),這里可以采用在原半空間的下方加上一層速度更高的層作為新半空間來計算頻散。但這樣會帶來一定的誤差[6],最好是取頻散方程復數(shù)根解的實部作為相速度來解決這一問題[7]。

    1.2 頻散反演理論和算法

    頻散曲線的反演是通過某種方式調(diào)節(jié)模型參數(shù)達到對實際頻散曲線的最好擬合。這里我們構(gòu)造目標函數(shù):

    Φ(m)=‖d-d(m)‖2,

    (1)

    即實際頻散相速度d與速度結(jié)構(gòu)模型m正演得到頻散相速度d(m)的L2范數(shù),來定義理論計算數(shù)據(jù)對實際數(shù)據(jù)的擬合程度。反演的過程實際上就是采用算法搜索調(diào)整m,使得目標函數(shù)Φ(m)為極小的過程,若m有多個,則取其目標函數(shù)值最小的那一個作為反演的最終結(jié)果。最初頻散曲線反演方法主要利用半波長法、拐點法、漸進線法等。隨著計算能力的進步,算法的開發(fā),越來越多反演方法應(yīng)用到頻散反演中,尤其是啟發(fā)式算法的出現(xiàn),讓反演效率和速度提升到一個新高度:凡友華等[8]使用的BFGS擬牛頓法引入到頻散曲線的反演當中,提出了基于基階模的擬牛頓法反演,在理論上進行了反演模擬;石耀霖等[9]以勒夫波為例討論從面波頻散反演地球內(nèi)部構(gòu)造的遺傳算法;翟佳羽等[10]在蟻群算法引入集中度的概念并對層狀介質(zhì)模型以及實測數(shù)據(jù)的頻散曲線進行反演,驗證了方法的有效性。這些頻散反演的方法大致可以分為以微分為基礎(chǔ)的方法、枚舉算法和隨機算法。前面提到的BFGS擬牛頓法即利用了目標函數(shù)的一階導數(shù)信息,應(yīng)用于頻散反演這種非線性化反演時,依賴初始模型的選擇,容易出現(xiàn)不收斂的現(xiàn)象。遺傳算法屬于全局搜索方法,容易陷入局部極小值中,且反演用時較長、效率較低,但可進行優(yōu)化。模擬退火法則是一種偽隨機算法,通過用溫度代表概率進行演算來搜尋最優(yōu)解,更依賴初始模型和參數(shù)選擇。

    郭飛在Github上提供的scikit-opt(https://github.com/guofei9987/scikit-opt)軟件包封裝了多種啟發(fā)式算法函數(shù),并且在軟件包的網(wǎng)絡(luò)在線說明文檔(https://scikit-opt.github.io/scikit-opt/#/zh/README)中給出了算法特性、詳細的參數(shù)說明和使用案例。通過對提供的算法函數(shù)分析,其中可用于頻散反演的有:遺傳算法、差分進化算法、粒子群算法和模擬退火法。

    2 算法設(shè)計和模型測試

    首先需要實現(xiàn)頻散正演的子程序,輸入速度模型和需要計算的周期數(shù)數(shù)組,調(diào)用pysurf96程序包函數(shù)進行正演計算:如使用的是surf96函數(shù),則要在模型下方增加一個更高速的半空間進行計算;如使用的是surfdisp96函數(shù),則要取返回值的實數(shù)部分。然后實現(xiàn)輸入不同模型時計算目標函數(shù)的子程序,最后在主程序?qū)崿F(xiàn)文件數(shù)據(jù)的讀寫、參數(shù)的設(shè)定和反演算法的運算。其中,Python的numpy擴展庫提供豐富的數(shù)組計算函數(shù),為我們計算提供了方便,代碼也顯得更加簡潔。例如在計算目標函數(shù)時需要進行L2范數(shù)運算,python語言可以用numpy.linalg.norm()方便快捷算出,省去了復雜數(shù)組計算的代碼(附錄示例2)。

    a—差分進化算法;b—遺傳算法;c—粒子群算法;d—模擬退火算法a—differential evolution algorithm;b—genetic algorithm;c—particle swarm optimization;d—simulated annealing algorithm圖2 不同啟發(fā)式算法反演的結(jié)果Fig.2 Results of different heuristic algorithms

    可以看出深度越淺,橫波速度結(jié)構(gòu)反演的精度越高,例如圖2a差分進化算法左圖第一、二、三層反演的效果都很好,且無論是高速層還是低速層,都很好地反演出來。這是因為瑞利面波高頻對淺層敏感;最后一層即半空間層的橫波速度也都反演的較好。圖2a右圖則是反演得到模型的頻散與實際頻散曲線的差別,可以看出擬合很好。至于中間兩層誤差較大,我們可以利用disba軟件計算瑞利面波的相速度敏感核進行探討,圖3是頻率為2 Hz的瑞利面波相速度敏感核。從圖中紅線和黃線可以看出,縱波速度和密度的改變總體都對相速度影響很?。痪G線的變化說明橫波速度在0~20 m淺層變化對相速度影響大,30~50 m范圍影響則較小,半空間層橫波速度變化則又對相速度起一定影響。這從另一個方向解釋了上述反演的結(jié)果,也可以看出各層反演得到速度結(jié)果的可靠性。

    圖3 頻率為2 Hz瑞利面波相速度敏感核Fig.3 Phase sensitivity kernel of Rayleigh wave on 2 Hz

    3 算法穩(wěn)定性和效率

    為了測試算法穩(wěn)定性和運行效率,這里所有運算過程均使用單CPU核心并且沒有其他計算任務(wù)在系統(tǒng)平臺上運行,系統(tǒng)平臺選擇Ubuntu以避免其他后臺程序干擾。測試使用的處理器為英特爾二代酷睿i5(i5-6200U)主頻2.3 GHz,且并沒有讓GPU參與運算。算法種群規(guī)模、迭代次數(shù)等參數(shù)均使用默認值,且均不加入外部判斷是否繼續(xù)或跳出反演代碼。測試結(jié)果如圖4所示。可以看出差分進化算法擬合度最高,但算法耗時最長;遺傳算法和粒子群算法反演耗時大幅度縮短,但擬合度較差;模擬退火法擬合度較好,耗時由于初始模型的選擇,這里接近差分進化算法??紤]到默認的種群規(guī)模較大,差分進化算法實際使用可以考慮減少數(shù)值來提升計算效率。遺傳算法和粒子群算法則可以加入合適判定增加反演迭代次數(shù)來提高擬合度的穩(wěn)定性。模擬退火法效率十分依靠初始模型的選擇,很適合已有初始模型的反演,但要注意模擬退火其他參數(shù)選擇也很重要,否則容易出現(xiàn)迭代次數(shù)或者溫度達到后跳出反演。

    a—差分進化算法;b—遺傳算法;c—粒子群算法;d—模擬退火算法a—differential evolution algorithm;b—genetic algorithm;c—particle swarm optimization;d—simulated annealing algorithm圖4 不同啟發(fā)式算法的效率和穩(wěn)定性Fig.4 Efficiency and stability of different heuristic algorithms

    scikit-opt軟件包還提供多CPU 目標函數(shù)計算和GPU加速功能,能成倍提升反演計算速度,并且可以進行斷點續(xù)行,能有效避免運算過程中部分函數(shù)出錯導致程序崩潰問題。隨著算力的提升,計算機科學的發(fā)展,能更好地進行地球物理學反演問題的研究。Zhang等[11]2021年應(yīng)用大數(shù)據(jù)先行計算好區(qū)域模型的格林函數(shù)并儲存,然后利于人工智能快速識別并反演震源破裂機制參數(shù),在智能預警地震上取得了突破的進展。如今Python語言已有眾多大數(shù)據(jù)和人工智能算法庫支持,需要我們提供封包好的地球物理學正演和反演方法。引入大數(shù)據(jù)和人工智能,這將是未來地球物理研究重要的方向。

    4 實例

    前文中提到彭一波[1]在研究中利用obspy軟件包提取了海拉爾盆地的兩個地震臺站間背景噪聲互相關(guān)格林函數(shù)的頻散曲線,這里我們對其頻散曲線進行一個多層厚度5 km的速度模型反演(圖5)??梢钥闯鰞蓚€臺站之間的地殼以及上地幔平均速度結(jié)構(gòu):區(qū)域的淺層橫波速度較小,波速在3 km/s附近,這可能是因為沉積層較厚;在15~25 km處存在低速速層,莫霍面深度在40 km附近。這些與李英康[12]、李皎皎[13]等人研究結(jié)論一致??紤]到背景噪聲提取的長周期面波頻散可靠性較低,上地幔以及深部結(jié)構(gòu)還需結(jié)合地震面波進行反演確定[14]。

    圖5 遺傳算法反演海拉爾盆地面瑞利波頻散結(jié)果Fig.5 Inversion of Rayleigh wave dispersion curve in Hailar Basin by Genetic Algorithm

    5 結(jié)論及討論

    本文基于Python語言scikit-opt軟件包進行瑞利面波頻散曲線反演實行方案的研究,給出了頻散正演算法在實際反演過程中所遇到頻散函數(shù)解虛數(shù)問題的解決方法,搭建了瑞利面波頻散反演地下層狀速度結(jié)構(gòu)的平臺,為今后研究背景噪聲成像、面波勘探提供開源基礎(chǔ)。同時對比了scikit-opt軟件包中不同啟發(fā)式算法的效率和穩(wěn)定性,并利用相速度敏感核探討了反演結(jié)果的可靠性,為使用者在反演算法的選擇上提供一些依據(jù)。最后對彭一波[1]在海拉爾盆地背景噪聲研究提取的頻散曲線進行了反演,得到了該區(qū)域的地殼和上地幔層狀速度模型。

    可以看出地球物理正演和反演算法的封包,能為廣大研究者提供一個可靠且高效的平臺,也為今后接入大數(shù)據(jù)和人工智能打下基礎(chǔ)。本文只進行了基階瑞利面波頻散的一維速度模型反演,現(xiàn)在研究趨勢是加入高階頻散進行二維、三維模型的反演,這是將今后繼續(xù)研究的方向。

    致謝:感謝R.Hermann、郭飛等人算法的開源共享。

    附錄

    示例1:基于pysurf96程序包計算瑞利面波基階頻散腳本示例,通過修改mode數(shù)值可以計算高階瑞利面波頻散。

    示例2:基于scikit-opt程序包中遺傳算法反演一維層狀速度模型腳本示例。腳本目錄需包含三個文件:模型上下邊界模型l.model和u.model;反演的目標頻散文件m.disp。反演的速度模型結(jié)果寫入本地文件fit.model。

    猜你喜歡
    面波瑞利橫波
    亞瑞利散斑場的二階累積量鬼成像
    橫波技術(shù)在工程物探中的應(yīng)用分析
    gPhone重力儀的面波頻段響應(yīng)實測研究
    地震研究(2021年1期)2021-04-13 01:04:56
    自適應(yīng)相減和Curvelet變換組合壓制面波
    馬瑞利推出多項汽車零部件技術(shù)
    汽車零部件(2015年4期)2015-12-22 05:32:56
    瑞利波頻散成像方法的實現(xiàn)及成像效果對比研究
    揚眉一顧,妖嬈橫波處
    橫波一顧,傲殺人間萬戶侯
    火花(2015年1期)2015-02-27 07:40:24
    橫波淺層地震在城市勘探中的應(yīng)用
    淺析工程勘探的面波勘探方法
    河南科技(2014年8期)2014-02-27 14:07:40
    成人亚洲精品一区在线观看| 宅男免费午夜| 国产熟女欧美一区二区| 两性夫妻黄色片| xxx大片免费视频| 男男h啪啪无遮挡| 婷婷色av中文字幕| 日韩电影二区| 精品国产乱码久久久久久男人| 欧美日韩一级在线毛片| www.熟女人妻精品国产| 国产精品熟女久久久久浪| 亚洲精品一区蜜桃| 成年人免费黄色播放视频| 日韩av在线免费看完整版不卡| 国产成人欧美在线观看 | 久久 成人 亚洲| 97精品久久久久久久久久精品| 2021少妇久久久久久久久久久| 亚洲精品日本国产第一区| 嫩草影视91久久| 黄片小视频在线播放| a级毛片在线看网站| 中文精品一卡2卡3卡4更新| 久久精品国产综合久久久| 三上悠亚av全集在线观看| 国产精品二区激情视频| 精品免费久久久久久久清纯 | 2018国产大陆天天弄谢| 丁香六月欧美| 黄色 视频免费看| 高清视频免费观看一区二区| 国产精品国产三级专区第一集| 免费av中文字幕在线| 久久久久网色| av在线老鸭窝| 国产免费现黄频在线看| 99香蕉大伊视频| 久久久久久久精品精品| 一区二区三区激情视频| 亚洲视频免费观看视频| 母亲3免费完整高清在线观看| xxx大片免费视频| 超碰成人久久| 在线观看人妻少妇| 免费高清在线观看日韩| 午夜老司机福利片| 美女福利国产在线| 成年av动漫网址| 黄色一级大片看看| 日本午夜av视频| 观看美女的网站| 激情视频va一区二区三区| 亚洲第一青青草原| 又大又爽又粗| www.精华液| 国产精品久久久久久精品电影小说| 国产色婷婷99| 亚洲国产精品一区二区三区在线| 狂野欧美激情性xxxx| 欧美人与性动交α欧美软件| 国产亚洲av片在线观看秒播厂| 欧美日韩综合久久久久久| 一区二区av电影网| 毛片一级片免费看久久久久| 国产在线一区二区三区精| 色94色欧美一区二区| 久久久久精品久久久久真实原创| 久久国产亚洲av麻豆专区| 女的被弄到高潮叫床怎么办| 十八禁网站网址无遮挡| av在线观看视频网站免费| 精品久久蜜臀av无| 男女高潮啪啪啪动态图| 午夜福利一区二区在线看| 国产成人精品在线电影| av在线观看视频网站免费| 亚洲av福利一区| av卡一久久| 男人爽女人下面视频在线观看| 老汉色∧v一级毛片| 男人添女人高潮全过程视频| 人妻一区二区av| 国产av精品麻豆| 亚洲精品久久久久久婷婷小说| 极品人妻少妇av视频| 久久狼人影院| 亚洲婷婷狠狠爱综合网| 国产片特级美女逼逼视频| 亚洲精品中文字幕在线视频| 97在线人人人人妻| 欧美变态另类bdsm刘玥| 午夜福利免费观看在线| 日韩av不卡免费在线播放| 亚洲一卡2卡3卡4卡5卡精品中文| 一本色道久久久久久精品综合| 伦理电影免费视频| 国产99久久九九免费精品| 国产在视频线精品| 国产日韩欧美视频二区| 纯流量卡能插随身wifi吗| 免费黄频网站在线观看国产| 18禁国产床啪视频网站| 久久女婷五月综合色啪小说| 国产女主播在线喷水免费视频网站| av在线观看视频网站免费| 各种免费的搞黄视频| 午夜免费鲁丝| 高清欧美精品videossex| 777久久人妻少妇嫩草av网站| 亚洲精品国产av蜜桃| 一区福利在线观看| 欧美日本中文国产一区发布| 狠狠婷婷综合久久久久久88av| 一级片免费观看大全| 波多野结衣一区麻豆| 久久久久视频综合| 黄色怎么调成土黄色| 国产精品嫩草影院av在线观看| 国产精品99久久99久久久不卡 | 欧美日本中文国产一区发布| 免费观看人在逋| 亚洲精品视频女| 极品人妻少妇av视频| 中文精品一卡2卡3卡4更新| 久久人人爽av亚洲精品天堂| 久久精品国产综合久久久| 男人爽女人下面视频在线观看| 国产精品国产三级国产专区5o| 肉色欧美久久久久久久蜜桃| 亚洲一区二区三区欧美精品| 亚洲国产av影院在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久人妻精品一区果冻| 中文字幕最新亚洲高清| 热re99久久国产66热| 午夜激情av网站| 国产成人欧美| 老熟女久久久| 美女中出高潮动态图| 18禁国产床啪视频网站| 久久久国产精品麻豆| 国产精品 国内视频| 国产一区二区三区av在线| 一本色道久久久久久精品综合| 校园人妻丝袜中文字幕| 精品卡一卡二卡四卡免费| 欧美 亚洲 国产 日韩一| 久久久精品区二区三区| 亚洲熟女毛片儿| av线在线观看网站| 精品亚洲成国产av| 亚洲男人天堂网一区| 亚洲成人av在线免费| 国产熟女欧美一区二区| 纵有疾风起免费观看全集完整版| 成人黄色视频免费在线看| 国产精品免费大片| 99久久99久久久精品蜜桃| 少妇 在线观看| 日韩 欧美 亚洲 中文字幕| 成年人免费黄色播放视频| 国产黄色免费在线视频| 自拍欧美九色日韩亚洲蝌蚪91| 热re99久久国产66热| 欧美亚洲日本最大视频资源| 欧美变态另类bdsm刘玥| 肉色欧美久久久久久久蜜桃| 日韩不卡一区二区三区视频在线| 久久久久人妻精品一区果冻| 啦啦啦在线观看免费高清www| 亚洲情色 制服丝袜| 亚洲国产欧美网| 国产精品一二三区在线看| 日韩 欧美 亚洲 中文字幕| 黑丝袜美女国产一区| 一区二区三区精品91| 一边亲一边摸免费视频| 少妇猛男粗大的猛烈进出视频| 在线免费观看不下载黄p国产| 亚洲国产毛片av蜜桃av| 久久女婷五月综合色啪小说| 亚洲av国产av综合av卡| 成人国语在线视频| 丰满迷人的少妇在线观看| a 毛片基地| 久久精品国产a三级三级三级| 欧美黑人欧美精品刺激| 汤姆久久久久久久影院中文字幕| 爱豆传媒免费全集在线观看| 国产成人a∨麻豆精品| 999精品在线视频| 男女床上黄色一级片免费看| 精品福利永久在线观看| 免费不卡黄色视频| 美女福利国产在线| 亚洲国产成人一精品久久久| 宅男免费午夜| 男女边摸边吃奶| 国产亚洲欧美精品永久| 亚洲综合精品二区| 久久久久久人妻| 免费在线观看完整版高清| 亚洲精品aⅴ在线观看| 久久这里只有精品19| 观看美女的网站| 麻豆乱淫一区二区| 欧美97在线视频| 国产高清不卡午夜福利| 观看av在线不卡| 精品第一国产精品| www.av在线官网国产| 亚洲七黄色美女视频| 99久久精品国产亚洲精品| 久久韩国三级中文字幕| 亚洲精品国产区一区二| 免费在线观看完整版高清| 人成视频在线观看免费观看| 亚洲精品乱久久久久久| 国产1区2区3区精品| 国产一级毛片在线| 亚洲,欧美精品.| 欧美精品人与动牲交sv欧美| 国产男人的电影天堂91| 一二三四中文在线观看免费高清| 黄网站色视频无遮挡免费观看| 国产精品一二三区在线看| 国产乱来视频区| 久久天堂一区二区三区四区| 亚洲av日韩在线播放| 亚洲欧美成人综合另类久久久| 9色porny在线观看| 一级a爱视频在线免费观看| 国产免费又黄又爽又色| 性高湖久久久久久久久免费观看| 久久99精品国语久久久| 黑人猛操日本美女一级片| 一本大道久久a久久精品| 悠悠久久av| 18禁国产床啪视频网站| 日韩av在线免费看完整版不卡| 成人三级做爰电影| 一级毛片我不卡| 18禁动态无遮挡网站| 日韩 欧美 亚洲 中文字幕| 成人亚洲欧美一区二区av| 最黄视频免费看| 色吧在线观看| 一级毛片电影观看| 一本久久精品| 波多野结衣av一区二区av| 久久久国产精品麻豆| 久久久久精品久久久久真实原创| 高清不卡的av网站| 纯流量卡能插随身wifi吗| 久久精品国产亚洲av涩爱| 少妇被粗大猛烈的视频| 美女福利国产在线| 国产爽快片一区二区三区| 老汉色av国产亚洲站长工具| 国产精品偷伦视频观看了| 久久午夜综合久久蜜桃| 一级黄片播放器| 午夜福利一区二区在线看| 在线天堂最新版资源| 女人精品久久久久毛片| 一级毛片我不卡| 亚洲国产精品国产精品| 国产在线视频一区二区| 91国产中文字幕| 国产成人av激情在线播放| 国产精品一二三区在线看| 人人妻人人爽人人添夜夜欢视频| 亚洲精品日韩在线中文字幕| 欧美中文综合在线视频| www日本在线高清视频| 搡老岳熟女国产| 亚洲精品aⅴ在线观看| 在线观看三级黄色| 成年人午夜在线观看视频| 精品亚洲乱码少妇综合久久| 老司机影院毛片| 国产一区二区在线观看av| 好男人视频免费观看在线| 在线免费观看不下载黄p国产| 少妇被粗大的猛进出69影院| 纵有疾风起免费观看全集完整版| 熟女少妇亚洲综合色aaa.| 亚洲精品日本国产第一区| 国产色婷婷99| 亚洲国产欧美一区二区综合| 久久久国产一区二区| 免费黄色在线免费观看| 亚洲欧美成人精品一区二区| 婷婷色麻豆天堂久久| 国产精品三级大全| 老司机在亚洲福利影院| 另类亚洲欧美激情| 大香蕉久久成人网| 国产成人av激情在线播放| 丝袜在线中文字幕| 性高湖久久久久久久久免费观看| 日韩伦理黄色片| 中文字幕高清在线视频| 亚洲图色成人| 国产成人精品福利久久| 亚洲美女视频黄频| 国语对白做爰xxxⅹ性视频网站| 亚洲伊人色综图| 五月开心婷婷网| 亚洲欧洲国产日韩| netflix在线观看网站| 精品亚洲成国产av| 肉色欧美久久久久久久蜜桃| 少妇的丰满在线观看| 免费日韩欧美在线观看| 国产免费又黄又爽又色| 久久久久久久久久久久大奶| 国产片内射在线| 国产无遮挡羞羞视频在线观看| 久久av网站| 国产 一区精品| 日本爱情动作片www.在线观看| 三上悠亚av全集在线观看| 日韩熟女老妇一区二区性免费视频| 少妇人妻久久综合中文| 亚洲四区av| 日韩电影二区| 免费高清在线观看视频在线观看| 我的亚洲天堂| 亚洲欧美激情在线| 久久久久精品人妻al黑| 哪个播放器可以免费观看大片| 日日爽夜夜爽网站| 久久久久久久大尺度免费视频| 91精品三级在线观看| 人人妻人人澡人人看| 一本久久精品| 久久精品熟女亚洲av麻豆精品| 欧美精品一区二区大全| 99热网站在线观看| 2021少妇久久久久久久久久久| 精品视频人人做人人爽| 纵有疾风起免费观看全集完整版| 国产人伦9x9x在线观看| 欧美日韩亚洲高清精品| 国产av精品麻豆| 一边亲一边摸免费视频| 日韩一区二区三区影片| 日本黄色日本黄色录像| 国产97色在线日韩免费| 水蜜桃什么品种好| 国产成人一区二区在线| 精品一区二区三区四区五区乱码 | 丝袜美腿诱惑在线| 午夜久久久在线观看| 另类精品久久| 国产精品.久久久| 最近最新中文字幕免费大全7| 午夜老司机福利片| 日韩中文字幕视频在线看片| 日本欧美国产在线视频| 人妻 亚洲 视频| 成人黄色视频免费在线看| 精品卡一卡二卡四卡免费| 一区二区日韩欧美中文字幕| 亚洲美女黄色视频免费看| 一区在线观看完整版| 国产成人a∨麻豆精品| 一个人免费看片子| 久久久久久久久久久久大奶| 最近中文字幕高清免费大全6| 天天躁夜夜躁狠狠躁躁| 免费观看人在逋| 精品一区在线观看国产| 如何舔出高潮| 99国产综合亚洲精品| √禁漫天堂资源中文www| 女人被躁到高潮嗷嗷叫费观| 亚洲国产日韩一区二区| 十八禁人妻一区二区| 国产极品粉嫩免费观看在线| 交换朋友夫妻互换小说| 亚洲成人av在线免费| 国产精品av久久久久免费| 日韩熟女老妇一区二区性免费视频| 亚洲精品美女久久av网站| 日韩精品有码人妻一区| 精品国产乱码久久久久久小说| 夜夜骑夜夜射夜夜干| 欧美变态另类bdsm刘玥| 欧美日韩成人在线一区二区| 卡戴珊不雅视频在线播放| 亚洲五月色婷婷综合| 久久精品aⅴ一区二区三区四区| 久久精品国产亚洲av高清一级| 美女脱内裤让男人舔精品视频| 精品一区在线观看国产| 99久久人妻综合| 亚洲精品视频女| 久久精品国产综合久久久| kizo精华| 日本一区二区免费在线视频| 中文乱码字字幕精品一区二区三区| 亚洲欧美成人精品一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| netflix在线观看网站| 日韩欧美一区视频在线观看| 国产极品粉嫩免费观看在线| 久热爱精品视频在线9| 国产高清不卡午夜福利| av天堂久久9| 亚洲少妇的诱惑av| 又大又爽又粗| 午夜福利影视在线免费观看| 深夜精品福利| 亚洲av福利一区| 国产成人欧美| 大话2 男鬼变身卡| 天堂8中文在线网| 欧美av亚洲av综合av国产av | 国产 精品1| 九草在线视频观看| www.精华液| 亚洲图色成人| 最新的欧美精品一区二区| 丝袜美足系列| 亚洲人成网站在线观看播放| 日韩制服骚丝袜av| 我要看黄色一级片免费的| 99国产精品免费福利视频| 777米奇影视久久| 一区二区三区激情视频| 天堂中文最新版在线下载| 色视频在线一区二区三区| 91成人精品电影| 大香蕉久久成人网| 99久久99久久久精品蜜桃| 黄片小视频在线播放| 少妇的丰满在线观看| 成年av动漫网址| 成年美女黄网站色视频大全免费| 亚洲专区中文字幕在线 | 国产熟女午夜一区二区三区| 国产片特级美女逼逼视频| 久久久国产精品麻豆| 精品第一国产精品| 秋霞伦理黄片| 亚洲国产av新网站| 亚洲精品在线美女| 亚洲av欧美aⅴ国产| 一级毛片我不卡| 国产在线免费精品| 久久ye,这里只有精品| 丝袜喷水一区| 免费在线观看黄色视频的| 麻豆乱淫一区二区| 久久精品国产亚洲av高清一级| 中文字幕精品免费在线观看视频| 国产精品人妻久久久影院| 纯流量卡能插随身wifi吗| 天天躁夜夜躁狠狠久久av| 精品国产国语对白av| 国产探花极品一区二区| 欧美少妇被猛烈插入视频| 久久人人97超碰香蕉20202| 只有这里有精品99| 亚洲精品日本国产第一区| 日本黄色日本黄色录像| 精品一区在线观看国产| 欧美av亚洲av综合av国产av | 丝袜脚勾引网站| 性高湖久久久久久久久免费观看| 青春草亚洲视频在线观看| 亚洲av在线观看美女高潮| 91国产中文字幕| 国产精品99久久99久久久不卡 | 精品少妇内射三级| 1024香蕉在线观看| 国产一区亚洲一区在线观看| 丝瓜视频免费看黄片| 亚洲精品,欧美精品| 黄色视频在线播放观看不卡| 老司机影院毛片| 久久99热这里只频精品6学生| 国产av一区二区精品久久| 大陆偷拍与自拍| 久久影院123| 在线观看免费视频网站a站| 成人国语在线视频| 啦啦啦啦在线视频资源| 一级,二级,三级黄色视频| 在线 av 中文字幕| 午夜福利免费观看在线| 国产日韩一区二区三区精品不卡| 97在线人人人人妻| 男女下面插进去视频免费观看| 日日摸夜夜添夜夜爱| 亚洲精华国产精华液的使用体验| 操美女的视频在线观看| 麻豆精品久久久久久蜜桃| 国产一区二区激情短视频 | 午夜福利一区二区在线看| 一级,二级,三级黄色视频| 亚洲婷婷狠狠爱综合网| 99re6热这里在线精品视频| 精品一区在线观看国产| 日韩电影二区| 各种免费的搞黄视频| 在线观看人妻少妇| 七月丁香在线播放| 中文字幕亚洲精品专区| 国产精品偷伦视频观看了| 男人爽女人下面视频在线观看| 大陆偷拍与自拍| 亚洲色图综合在线观看| 91精品三级在线观看| 欧美精品av麻豆av| 国产色婷婷99| 中文字幕色久视频| 涩涩av久久男人的天堂| 操出白浆在线播放| 亚洲国产成人一精品久久久| 久久性视频一级片| 久久久久网色| 黑丝袜美女国产一区| 日日爽夜夜爽网站| 老汉色av国产亚洲站长工具| 观看美女的网站| 亚洲第一青青草原| 中国国产av一级| 午夜免费男女啪啪视频观看| 日韩av免费高清视频| 日韩电影二区| 2018国产大陆天天弄谢| 日韩 欧美 亚洲 中文字幕| 成年女人毛片免费观看观看9 | 91aial.com中文字幕在线观看| 国产精品一二三区在线看| 欧美日韩成人在线一区二区| 麻豆乱淫一区二区| 亚洲人成电影观看| 麻豆乱淫一区二区| 汤姆久久久久久久影院中文字幕| 狂野欧美激情性bbbbbb| 高清黄色对白视频在线免费看| 99久国产av精品国产电影| 欧美xxⅹ黑人| av女优亚洲男人天堂| 亚洲国产欧美一区二区综合| 老司机在亚洲福利影院| 日韩一本色道免费dvd| 日韩 亚洲 欧美在线| 日韩人妻精品一区2区三区| 欧美日韩av久久| 久久免费观看电影| 亚洲五月色婷婷综合| 免费人妻精品一区二区三区视频| 少妇人妻精品综合一区二区| 亚洲av欧美aⅴ国产| 亚洲欧美精品综合一区二区三区| 男女边摸边吃奶| 麻豆乱淫一区二区| 亚洲精品视频女| 亚洲欧美激情在线| 大香蕉久久成人网| 精品少妇一区二区三区视频日本电影 | 新久久久久国产一级毛片| 一本色道久久久久久精品综合| av在线老鸭窝| 亚洲综合色网址| 久久精品熟女亚洲av麻豆精品| 国产男女内射视频| 成人午夜精彩视频在线观看| 在线观看国产h片| 在线精品无人区一区二区三| 亚洲 欧美一区二区三区| 90打野战视频偷拍视频| 王馨瑶露胸无遮挡在线观看| 欧美日韩亚洲综合一区二区三区_| 美女福利国产在线| 美女脱内裤让男人舔精品视频| 久久精品亚洲熟妇少妇任你| 美女大奶头黄色视频| 亚洲欧美日韩另类电影网站| 成人手机av| 国产一级毛片在线| 国产一区二区在线观看av| 在线观看一区二区三区激情| 日本猛色少妇xxxxx猛交久久| 另类亚洲欧美激情| 成人漫画全彩无遮挡| 亚洲国产日韩一区二区| 王馨瑶露胸无遮挡在线观看| 亚洲一码二码三码区别大吗| 免费黄色在线免费观看| 欧美日韩视频精品一区| 妹子高潮喷水视频| 久久久久精品久久久久真实原创| 亚洲少妇的诱惑av| 2018国产大陆天天弄谢| 国产又色又爽无遮挡免| 日本爱情动作片www.在线观看| 午夜福利免费观看在线| 国产精品嫩草影院av在线观看| 日韩一本色道免费dvd| 亚洲精品在线美女| 又大又黄又爽视频免费| 国产精品麻豆人妻色哟哟久久| 一级毛片电影观看| 十八禁网站网址无遮挡| 亚洲成人手机| 永久免费av网站大全| 青春草国产在线视频| 欧美日韩福利视频一区二区|