• <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
    精品亚洲乱码少妇综合久久| 看十八女毛片水多多多| 亚洲成人av在线免费| 晚上一个人看的免费电影| 女的被弄到高潮叫床怎么办| 中文字幕人妻熟人妻熟丝袜美| 国产成人精品一,二区| 午夜免费观看性视频| 特级一级黄色大片| av.在线天堂| 亚洲欧美精品自产自拍| 亚洲精品自拍成人| 亚洲精品视频女| 久久99蜜桃精品久久| 国产精品av视频在线免费观看| 黑人高潮一二区| 少妇的逼好多水| 国产一区二区亚洲精品在线观看| 精品久久久久久久末码| 狂野欧美白嫩少妇大欣赏| 天堂俺去俺来也www色官网 | 日本爱情动作片www.在线观看| 啦啦啦中文免费视频观看日本| 欧美日韩国产mv在线观看视频 | 美女黄网站色视频| 国产成人午夜福利电影在线观看| 欧美激情在线99| 色哟哟·www| 欧美bdsm另类| 毛片女人毛片| 免费大片18禁| 丰满人妻一区二区三区视频av| 亚洲自拍偷在线| 男女边摸边吃奶| 欧美日韩一区二区视频在线观看视频在线 | 成人午夜精彩视频在线观看| 中文资源天堂在线| 精品一区二区三区人妻视频| 亚洲成人中文字幕在线播放| 哪个播放器可以免费观看大片| 黄片wwwwww| 国产 亚洲一区二区三区 | 蜜臀久久99精品久久宅男| 日韩不卡一区二区三区视频在线| 久久久精品欧美日韩精品| 欧美zozozo另类| 嘟嘟电影网在线观看| 国产色爽女视频免费观看| 久久久久免费精品人妻一区二区| 夜夜看夜夜爽夜夜摸| 自拍偷自拍亚洲精品老妇| 亚洲国产精品专区欧美| 搡老乐熟女国产| 在线观看免费高清a一片| 在线免费十八禁| 免费看av在线观看网站| 黄色欧美视频在线观看| 亚洲av二区三区四区| 日韩人妻高清精品专区| 国产乱人视频| 亚洲av男天堂| 日日摸夜夜添夜夜爱| 亚洲国产色片| 亚洲精品乱码久久久v下载方式| 青春草视频在线免费观看| 久久久久国产网址| 亚洲精品,欧美精品| 亚洲av男天堂| 嫩草影院入口| 少妇熟女欧美另类| 成年女人在线观看亚洲视频 | 国产麻豆成人av免费视频| 亚洲精品国产av蜜桃| 成人一区二区视频在线观看| 网址你懂的国产日韩在线| 国产高潮美女av| 最新中文字幕久久久久| 一夜夜www| 日韩亚洲欧美综合| 99热这里只有精品一区| 久久久久久久久久黄片| 一级毛片我不卡| 精品亚洲乱码少妇综合久久| 欧美一级a爱片免费观看看| 国产一级毛片七仙女欲春2| 男女啪啪激烈高潮av片| 欧美性猛交╳xxx乱大交人| 最近手机中文字幕大全| 亚洲丝袜综合中文字幕| 亚洲经典国产精华液单| 自拍偷自拍亚洲精品老妇| 黄片wwwwww| 七月丁香在线播放| 国产综合精华液| 成人无遮挡网站| 国产单亲对白刺激| 亚洲精华国产精华液的使用体验| 一区二区三区高清视频在线| 亚洲自偷自拍三级| 一本久久精品| 亚洲人成网站高清观看| 亚洲电影在线观看av| 美女大奶头视频| 久久久国产一区二区| 少妇猛男粗大的猛烈进出视频 | 免费看av在线观看网站| 免费少妇av软件| 久久精品熟女亚洲av麻豆精品 | 男人舔女人下体高潮全视频| 中国国产av一级| 国产精品熟女久久久久浪| 成人亚洲精品av一区二区| 小蜜桃在线观看免费完整版高清| 亚洲精品aⅴ在线观看| 国产伦精品一区二区三区视频9| 国产亚洲最大av| 成人无遮挡网站| 18禁动态无遮挡网站| 天堂av国产一区二区熟女人妻| 欧美最新免费一区二区三区| 亚洲欧美成人综合另类久久久| 深爱激情五月婷婷| 男人爽女人下面视频在线观看| 内地一区二区视频在线| 亚洲精品一区蜜桃| 成年版毛片免费区| 婷婷色综合大香蕉| 国产欧美另类精品又又久久亚洲欧美| 搡老妇女老女人老熟妇| 极品教师在线视频| 日韩,欧美,国产一区二区三区| 日韩av在线免费看完整版不卡| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产爱豆传媒在线观看| 免费av毛片视频| 久久综合国产亚洲精品| 日日啪夜夜爽| 精品国产露脸久久av麻豆 | 中文天堂在线官网| 亚洲综合色惰| 国产免费视频播放在线视频 | 99久久人妻综合| 中文字幕人妻熟人妻熟丝袜美| 欧美+日韩+精品| 国产精品一区二区在线观看99 | 久久精品久久精品一区二区三区| 乱系列少妇在线播放| 亚洲在线自拍视频| 全区人妻精品视频| av国产免费在线观看| 欧美另类一区| 久久久久九九精品影院| 国产91av在线免费观看| 国产一区二区三区综合在线观看 | 亚洲aⅴ乱码一区二区在线播放| 男女视频在线观看网站免费| 国产在线一区二区三区精| 真实男女啪啪啪动态图| 国内精品美女久久久久久| 亚洲国产精品专区欧美| 一级毛片我不卡| 亚洲不卡免费看| 欧美+日韩+精品| 伦理电影大哥的女人| 一个人看的www免费观看视频| 久久久久久久久久久免费av| 免费看光身美女| 亚洲最大成人中文| 国产成人精品福利久久| 最近中文字幕2019免费版| 国产老妇女一区| 久久久亚洲精品成人影院| 街头女战士在线观看网站| 最新中文字幕久久久久| av又黄又爽大尺度在线免费看| 亚洲熟妇中文字幕五十中出| 最近视频中文字幕2019在线8| 91在线精品国自产拍蜜月| 一本久久精品| 天堂俺去俺来也www色官网 | 精品久久久久久电影网| freevideosex欧美| 国产一区二区在线观看日韩| 亚洲18禁久久av| 又粗又硬又长又爽又黄的视频| 精品久久国产蜜桃| 在线观看美女被高潮喷水网站| 91精品伊人久久大香线蕉| 一区二区三区乱码不卡18| 欧美zozozo另类| 亚洲av二区三区四区| 91在线精品国自产拍蜜月| 国产视频首页在线观看| 美女黄网站色视频| 久久综合国产亚洲精品| 99久久精品国产国产毛片| 午夜激情欧美在线| 秋霞在线观看毛片| 午夜精品国产一区二区电影 | 大香蕉久久网| 国产亚洲精品av在线| 国产白丝娇喘喷水9色精品| 午夜福利在线观看免费完整高清在| 男的添女的下面高潮视频| 久久久午夜欧美精品| 亚洲怡红院男人天堂| 有码 亚洲区| 亚洲av一区综合| 久久97久久精品| 22中文网久久字幕| 一区二区三区高清视频在线| 日韩一本色道免费dvd| 国产一级毛片七仙女欲春2| eeuss影院久久| 有码 亚洲区| 乱人视频在线观看| 国语对白做爰xxxⅹ性视频网站| 好男人在线观看高清免费视频| 亚洲国产成人一精品久久久| 80岁老熟妇乱子伦牲交| 日韩,欧美,国产一区二区三区| 精品酒店卫生间| 精品少妇黑人巨大在线播放| 在线天堂最新版资源| 亚洲av国产av综合av卡| 日韩一本色道免费dvd| 一级毛片黄色毛片免费观看视频| 日韩av不卡免费在线播放| 国产乱人视频| 精品久久久久久久久亚洲| 免费不卡的大黄色大毛片视频在线观看 | 真实男女啪啪啪动态图| 小蜜桃在线观看免费完整版高清| 亚洲第一区二区三区不卡| 天天躁夜夜躁狠狠久久av| 日日啪夜夜撸| 有码 亚洲区| 丝袜美腿在线中文| 日韩成人av中文字幕在线观看| 99热全是精品| 亚洲精品国产av蜜桃| 日韩,欧美,国产一区二区三区| 亚洲精品色激情综合| 禁无遮挡网站| 日韩电影二区| 日日啪夜夜撸| 女人十人毛片免费观看3o分钟| 亚洲美女搞黄在线观看| 美女脱内裤让男人舔精品视频| 日韩av免费高清视频| 亚洲一区高清亚洲精品| 欧美日韩综合久久久久久| 国产精品三级大全| 久久久久久久国产电影| 国产免费福利视频在线观看| 一级二级三级毛片免费看| 九九爱精品视频在线观看| 久久久久免费精品人妻一区二区| 国产成人一区二区在线| 最近最新中文字幕大全电影3| 国产色婷婷99| 欧美最新免费一区二区三区| 哪个播放器可以免费观看大片| 国产国拍精品亚洲av在线观看| 欧美zozozo另类| 国产v大片淫在线免费观看| 熟妇人妻久久中文字幕3abv| 一区二区三区高清视频在线| 亚洲精品色激情综合| 亚洲欧洲日产国产| 有码 亚洲区| 国产人妻一区二区三区在| 国产成人精品一,二区| 边亲边吃奶的免费视频| 国产黄色小视频在线观看| 91精品国产九色| 一个人看的www免费观看视频| 日韩制服骚丝袜av| 一区二区三区乱码不卡18| 中文字幕人妻熟人妻熟丝袜美| 高清午夜精品一区二区三区| 777米奇影视久久| 中文乱码字字幕精品一区二区三区 | 欧美潮喷喷水| 亚洲成人精品中文字幕电影| 成人无遮挡网站| 国产免费又黄又爽又色| 亚洲电影在线观看av| 91av网一区二区| av卡一久久| 精品人妻一区二区三区麻豆| 美女脱内裤让男人舔精品视频| 免费电影在线观看免费观看| 在线免费十八禁| 欧美+日韩+精品| 中文天堂在线官网| 两个人视频免费观看高清| 亚洲精品国产av蜜桃| 成人亚洲欧美一区二区av| xxx大片免费视频| 国产乱来视频区| 亚洲色图av天堂| 日日撸夜夜添| 国产精品美女特级片免费视频播放器| 白带黄色成豆腐渣| 久久久久久久久久黄片| 在线观看免费高清a一片| 非洲黑人性xxxx精品又粗又长| 日本猛色少妇xxxxx猛交久久| 国产综合精华液| 国产精品美女特级片免费视频播放器| 日本黄色片子视频| 午夜免费激情av| 久久草成人影院| 中文字幕久久专区| 肉色欧美久久久久久久蜜桃 | 成年人午夜在线观看视频 | 久久久久性生活片| 色尼玛亚洲综合影院| 不卡视频在线观看欧美| 亚洲成人精品中文字幕电影| 少妇的逼好多水| 国语对白做爰xxxⅹ性视频网站| 内地一区二区视频在线| 一个人看视频在线观看www免费| 男女啪啪激烈高潮av片| 亚洲av成人av| 久久人人爽人人片av| 91精品一卡2卡3卡4卡| 婷婷色麻豆天堂久久| 最近中文字幕高清免费大全6| 久久精品国产鲁丝片午夜精品| 国产免费一级a男人的天堂| 国产精品爽爽va在线观看网站| 国产单亲对白刺激| 欧美日韩综合久久久久久| 久久久精品94久久精品| 亚洲精品乱码久久久v下载方式| 成年女人看的毛片在线观看| 久久人人爽人人爽人人片va| 午夜福利成人在线免费观看| 中文天堂在线官网| 床上黄色一级片| 秋霞在线观看毛片| 国产精品1区2区在线观看.| 国产亚洲午夜精品一区二区久久 | 亚洲av免费在线观看| 麻豆精品久久久久久蜜桃| 国产又色又爽无遮挡免| 搡女人真爽免费视频火全软件| 男人和女人高潮做爰伦理| 特级一级黄色大片| 亚洲精品乱久久久久久| 国产成人aa在线观看| 在线观看人妻少妇| 亚洲国产精品成人久久小说| 一二三四中文在线观看免费高清| 韩国av在线不卡| 熟妇人妻不卡中文字幕| 波多野结衣巨乳人妻| 欧美精品一区二区大全| 国产成人91sexporn| 18禁在线播放成人免费| 好男人在线观看高清免费视频| 成人亚洲精品av一区二区| 国产亚洲av片在线观看秒播厂 | 日本午夜av视频| 午夜老司机福利剧场| 国产成人午夜福利电影在线观看| 国产在视频线精品| 国产探花在线观看一区二区| 免费观看a级毛片全部| 国产美女午夜福利| 熟女电影av网| 少妇裸体淫交视频免费看高清| 中文字幕久久专区| 色综合亚洲欧美另类图片| 又爽又黄无遮挡网站| av网站免费在线观看视频 | 欧美日韩亚洲高清精品| 亚洲精品aⅴ在线观看| 亚洲精品影视一区二区三区av| 久久久久久久久久成人| 亚洲av成人av| 又爽又黄a免费视频| 97人妻精品一区二区三区麻豆| 能在线免费观看的黄片| 熟女电影av网| 日韩av在线大香蕉| 午夜视频国产福利| 日韩欧美三级三区| 久久精品国产亚洲网站| 午夜福利在线观看免费完整高清在| 69人妻影院| 精品国内亚洲2022精品成人| 久久99蜜桃精品久久| 成人av在线播放网站| 国产成人午夜福利电影在线观看| 亚洲伊人久久精品综合| 麻豆精品久久久久久蜜桃| av在线蜜桃| 日日摸夜夜添夜夜爱| 国产成人a区在线观看| 亚洲乱码一区二区免费版| 国产精品久久视频播放| 色5月婷婷丁香| 国产精品一区二区三区四区久久| 国产一级毛片在线| 午夜日本视频在线| 2021少妇久久久久久久久久久| 欧美不卡视频在线免费观看| 91久久精品国产一区二区三区| 国产大屁股一区二区在线视频| 精华霜和精华液先用哪个| 国产人妻一区二区三区在| 九九在线视频观看精品| 日日干狠狠操夜夜爽| 免费av观看视频| 国产成年人精品一区二区| 亚洲一区高清亚洲精品| 看非洲黑人一级黄片| 直男gayav资源| 日韩欧美三级三区| 80岁老熟妇乱子伦牲交| 天美传媒精品一区二区| 看十八女毛片水多多多| 99久国产av精品国产电影| 国产爱豆传媒在线观看| 女人被狂操c到高潮| 26uuu在线亚洲综合色| 日日啪夜夜爽| 国产成人福利小说| 亚洲,欧美,日韩| 一区二区三区四区激情视频| .国产精品久久| 毛片女人毛片| 午夜福利在线在线| 中文字幕制服av| 久久久久久久久大av| 夜夜爽夜夜爽视频| 国产精品精品国产色婷婷| 久久鲁丝午夜福利片| 午夜福利在线观看吧| 日本免费在线观看一区| 午夜激情欧美在线| 看非洲黑人一级黄片| 午夜久久久久精精品| 亚洲自偷自拍三级| 久久亚洲国产成人精品v| 成人午夜高清在线视频| 欧美日韩综合久久久久久| 五月天丁香电影| 久久久午夜欧美精品| 欧美 日韩 精品 国产| 天天躁日日操中文字幕| 亚洲成人av在线免费| 国产精品嫩草影院av在线观看| 久久久久久九九精品二区国产| 久久99热6这里只有精品| 久久精品久久久久久噜噜老黄| 欧美成人午夜免费资源| 色综合亚洲欧美另类图片| 亚洲三级黄色毛片| 日本猛色少妇xxxxx猛交久久| av在线亚洲专区| 国产精品久久久久久av不卡| 两个人的视频大全免费| 女人被狂操c到高潮| 精品久久久久久电影网| 国产欧美日韩精品一区二区| 国产精品女同一区二区软件| h日本视频在线播放| 春色校园在线视频观看| 国产精品国产三级专区第一集| 嘟嘟电影网在线观看| 九九久久精品国产亚洲av麻豆| 午夜精品一区二区三区免费看| 色综合站精品国产| 色综合色国产| 国内精品宾馆在线| 亚洲欧美精品专区久久| videos熟女内射| 成人特级av手机在线观看| 少妇熟女aⅴ在线视频| 国产精品久久视频播放| 久久久亚洲精品成人影院| 国产v大片淫在线免费观看| 99久久中文字幕三级久久日本| 80岁老熟妇乱子伦牲交| 日韩av在线大香蕉| 边亲边吃奶的免费视频| 极品教师在线视频| 精品午夜福利在线看| 久久久欧美国产精品| 少妇猛男粗大的猛烈进出视频 | 亚洲av成人av| 日韩国内少妇激情av| 国产精品久久久久久av不卡| 国产精品综合久久久久久久免费| 欧美 日韩 精品 国产| 日韩三级伦理在线观看| 亚洲性久久影院| 日韩欧美精品免费久久| 日韩av不卡免费在线播放| 97精品久久久久久久久久精品| 美女大奶头视频| 男女那种视频在线观看| 超碰av人人做人人爽久久| 欧美区成人在线视频| 在线免费观看不下载黄p国产| 久久久久精品久久久久真实原创| 久久精品熟女亚洲av麻豆精品 | 精品国产三级普通话版| 成人欧美大片| 两个人视频免费观看高清| 国产成人精品久久久久久| 欧美3d第一页| 亚洲av免费在线观看| 亚洲性久久影院| 国产在视频线精品| 久久99蜜桃精品久久| 国产亚洲av片在线观看秒播厂 | 成人欧美大片| 欧美zozozo另类| 午夜福利网站1000一区二区三区| 午夜激情福利司机影院| 久久99蜜桃精品久久| 老女人水多毛片| 夫妻午夜视频| 少妇高潮的动态图| 亚洲国产色片| 久久精品夜夜夜夜夜久久蜜豆| 欧美97在线视频| 最近中文字幕高清免费大全6| 国产免费福利视频在线观看| 国产v大片淫在线免费观看| 亚洲国产欧美在线一区| 搞女人的毛片| 亚洲精品日本国产第一区| 哪个播放器可以免费观看大片| 黄色欧美视频在线观看| av在线天堂中文字幕| 欧美最新免费一区二区三区| 秋霞在线观看毛片| 99久久中文字幕三级久久日本| 国产伦理片在线播放av一区| 99久久精品国产国产毛片| 少妇人妻一区二区三区视频| 一级爰片在线观看| 一区二区三区四区激情视频| 久久久久久九九精品二区国产| 熟妇人妻不卡中文字幕| av一本久久久久| 黄色日韩在线| 国产视频内射| 日韩大片免费观看网站| 国产探花极品一区二区| 精品欧美国产一区二区三| 久久国产乱子免费精品| 青春草视频在线免费观看| 欧美zozozo另类| 91午夜精品亚洲一区二区三区| 精品亚洲乱码少妇综合久久| 国产亚洲最大av| 国产熟女欧美一区二区| 一级av片app| 777米奇影视久久| 精品酒店卫生间| 免费黄色在线免费观看| 久久久成人免费电影| 国产激情偷乱视频一区二区| 毛片一级片免费看久久久久| 汤姆久久久久久久影院中文字幕 | 精品久久久久久久久亚洲| 联通29元200g的流量卡| 国产综合精华液| 亚洲第一区二区三区不卡| 国产精品蜜桃在线观看| 久久精品综合一区二区三区| 色综合色国产| 三级国产精品片| 精品午夜福利在线看| 精品酒店卫生间| 亚洲精品成人av观看孕妇| 亚洲精品乱久久久久久| 十八禁网站网址无遮挡 | 精品不卡国产一区二区三区| 简卡轻食公司| 99热网站在线观看| 十八禁网站网址无遮挡 | 欧美bdsm另类| 欧美高清成人免费视频www| 99久久精品国产国产毛片| 国产一区二区在线观看日韩| 人人妻人人看人人澡| 免费观看的影片在线观看| 亚洲精品日韩av片在线观看| 国产激情偷乱视频一区二区| 18禁裸乳无遮挡免费网站照片| 日本色播在线视频| 赤兔流量卡办理| 国产亚洲av嫩草精品影院| 哪个播放器可以免费观看大片| 久久久精品免费免费高清| 欧美成人一区二区免费高清观看| 欧美日韩在线观看h| 午夜免费观看性视频| 一区二区三区乱码不卡18| 婷婷色综合大香蕉| 国产精品人妻久久久久久| 国产精品嫩草影院av在线观看| 亚洲美女视频黄频| 欧美日本视频| 在线免费十八禁| 一级毛片 在线播放|