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

    REFTEK連續(xù)波形截取程序

    2015-12-19 02:09:36查小惠張廣偉楊雪超
    地震地磁觀測(cè)與研究 2015年1期
    關(guān)鍵詞:走時(shí)文件夾臺(tái)站

    查小惠 張廣偉 楊雪超

    1)中國南昌 330039 江西省地震局

    2)中國北京100085 中國地震局地殼應(yīng)力研究所

    0 引言

    近年來,流動(dòng)地震觀測(cè)越來越多,流動(dòng)地震觀測(cè)數(shù)據(jù)的存儲(chǔ)管理和使用問題,受到越來越多的關(guān)注。2010年4月到2011年7月,中國地震局地殼應(yīng)力研究所在云南地區(qū)布設(shè)21個(gè)寬頻帶野外流動(dòng)觀測(cè)臺(tái)站,所用地震計(jì)為CMG-3ESPC,數(shù)據(jù)采集器為REFTEK-130B,最終獲得233.7 G連續(xù)波形資料(Lei et al, 2012)。在資料處理階段,由于研究目的和手段的差異,有時(shí)需要采用不同的地震目錄,選用不同的地震參數(shù),反復(fù)從連續(xù)波形數(shù)據(jù)中截取所需事件波形資料。基于研究需要,嘗試開發(fā)一套簡(jiǎn)易的REFTEK連續(xù)波形截取程序,以便從REFTEK格式的連續(xù)波形數(shù)據(jù)中截取任意長(zhǎng)度的SAC格式的地震事件波形。

    1 程序運(yùn)行平臺(tái)和功能

    REFTEK連續(xù)波形截取程序可以在個(gè)人計(jì)算機(jī)上安裝運(yùn)行,基于Linux系統(tǒng)下的Matlab平臺(tái)開發(fā),因?yàn)椋孩賚inux系統(tǒng)下有現(xiàn)成的可以截取REFTEK連續(xù)波形數(shù)據(jù)的arcfetch程序和將REFTEK格式數(shù)據(jù)轉(zhuǎn)換為SAC格式數(shù)據(jù)的ref2sac程序;②Matlab平臺(tái)有很多高級(jí)函數(shù),比如:datenum和datestr函數(shù)可以很好地解決波形截取過程中遇到的時(shí)間運(yùn)算問題。在處理SAC文件上,也有現(xiàn)成的saclab工具包可以使用,而且在Matlab中可以利用Unix命令方便調(diào)用Linux系統(tǒng)的程序?;贚inux系統(tǒng)下的Matlab平臺(tái)編寫程序,可以同時(shí)利用Matlab高級(jí)函數(shù)和Linux系統(tǒng)程序兩者的優(yōu)勢(shì),高效實(shí)現(xiàn)程序功能。

    本程序的功能是從REFTEK連續(xù)波形數(shù)據(jù)中截取事件波形,原始數(shù)據(jù)為流動(dòng)觀測(cè)記錄的REFTEK格式的連續(xù)波形數(shù)據(jù),截取得到SAC格式的三分量事件地震波形數(shù)據(jù)。目前程序的測(cè)試平臺(tái)為Ubuntu12.04和Matlab2009b,只要確保Matlab相關(guān)程序包的安裝和Linux系統(tǒng)子程序的成功調(diào)用,程序即可順利運(yùn)行。程序的運(yùn)行時(shí)間和地震目錄的個(gè)數(shù)及截取波形的長(zhǎng)度有關(guān),雖然調(diào)用了較多外部程序,但是程序整體運(yùn)行效率可以接受的。

    2 程序運(yùn)行流程

    在調(diào)用Matlab程序截取波形前,需要進(jìn)行相關(guān)數(shù)據(jù)準(zhǔn)備。首先是原始連續(xù)波形數(shù)據(jù),程序要求每個(gè)臺(tái)站記錄的REFTEK原始波形數(shù)據(jù)單獨(dú)存放在一個(gè)文件夾下,而且為了方便訪問該文件夾,使用阿拉伯?dāng)?shù)字給文件夾編號(hào)命名。然后是地震目錄的整理,地震目錄不限行數(shù),每行8列,1到3列分別為年、月、日,第4列為地震發(fā)震時(shí)刻,格式為時(shí)分秒,5到8列分別為震源的緯度、經(jīng)度、深度和震級(jí)大小,各列之間以空格分開。最后是準(zhǔn)備臺(tái)站信息,包括4列,分別為臺(tái)站編號(hào)、臺(tái)站經(jīng)度、臺(tái)站緯度和臺(tái)站高程。臺(tái)站信息文件直接控制截取哪些臺(tái)站的事件波形,如果不需要截取相關(guān)臺(tái)站的波形,則將該臺(tái)站的信息從臺(tái)站信息文件中刪除。臺(tái)站信息文件中的臺(tái)站編號(hào)應(yīng)該和存放該臺(tái)原始連續(xù)波形的文件夾命名編號(hào)一致,方便在波形截取中尋找對(duì)應(yīng)臺(tái)站的數(shù)據(jù)文件。

    程序?qū)嶋H為一個(gè)Matlab函數(shù),調(diào)用該函數(shù)需要輸入4個(gè)參數(shù),分別為P波(實(shí)際為初至波)到達(dá)前時(shí)間、截取波形的總長(zhǎng)度、相對(duì)于臺(tái)站的最小震中距和最大震中距(用來篩選地震目錄)。截取事件波形需要一個(gè)參考時(shí)間,本程序中參考時(shí)間是P波(初至波)到達(dá)臺(tái)站的理論到時(shí),函數(shù)輸入的第1個(gè)參數(shù)是相對(duì)于該參考時(shí)間的提前時(shí)間,結(jié)合截取波形的總時(shí)間長(zhǎng)度,計(jì)算截取波形的時(shí)間終點(diǎn)。程序運(yùn)行流程見圖1。

    圖1 程序流程Fig.1 Program fl ow diagram

    (1)從存放臺(tái)站信息的文件中,讀入臺(tái)站信息,包括臺(tái)站編號(hào)、臺(tái)站經(jīng)緯度和高程。然后從第1行開始循環(huán),即波形數(shù)據(jù)的截取是逐臺(tái)進(jìn)行的。

    (2)讀入地震目錄,結(jié)合臺(tái)站經(jīng)緯度信息和地震目錄中地震事件的經(jīng)緯度信息,可以對(duì)地震目錄進(jìn)行挑選,選擇標(biāo)準(zhǔn)是調(diào)用主程序時(shí)輸入的最小和最大震中距。如果截取的事件波形用于接收函數(shù)研究,一般將震中距限制在30°—95°。

    (3)選取地震目錄后,計(jì)算到時(shí)目錄。主要計(jì)算需要截取波形的前后兩個(gè)時(shí)間點(diǎn)的時(shí)間。利用TAUP軟件計(jì)算P波理論走時(shí),發(fā)震時(shí)刻加上P波理論走時(shí)就得到P波理論到時(shí),即波形截取的參考時(shí)間。然后根據(jù)調(diào)用函數(shù)時(shí)輸入的截取波形提前時(shí)間和截取波形總時(shí)間長(zhǎng)度,計(jì)算截取波形前后兩個(gè)端點(diǎn)的時(shí)間。

    (4)從第1個(gè)地震目錄開始,使用arcfetch程序,從REFTEK連續(xù)波形中截取已計(jì)算兩時(shí)間點(diǎn)間的事件波形,使用ref2sac程序,將REFTEK格式數(shù)據(jù)轉(zhuǎn)換為SAC格式數(shù)據(jù)。

    (5)使用saclab程序包處理SAC數(shù)據(jù),主要重寫SAC數(shù)據(jù)頭文件,需要注意SAC頭文件中的時(shí)間處理。

    (6)利用相關(guān)文件處理命令,對(duì)SAC文件進(jìn)行存儲(chǔ),并返回流程(4)截取下一個(gè)地震事件,地震目錄截取完畢,返回流程(1)截取下一個(gè)臺(tái)站。最后獲得的SAC數(shù)據(jù)仍然以臺(tái)站為單位進(jìn)行存放,也就是,同一個(gè)臺(tái)站截取的不同地震事件的SAC文件存放在同一個(gè)文件夾。如果想改為一個(gè)地震事件文件夾下存放多個(gè)臺(tái)站文件的數(shù)據(jù)格式,需要編寫相應(yīng)的腳本。

    3 關(guān)鍵子程序

    3.1 arcfetch

    arcfetch程序用來從原始的連續(xù)REFTEK波形數(shù)據(jù)中截取一段波形數(shù)據(jù),波形截取的起始時(shí)間和波形長(zhǎng)度可以在程序輸入?yún)?shù)中指定,且截取數(shù)據(jù)格式仍然為REFTEK數(shù)據(jù)格式。使用該程序時(shí),需要指定原始連續(xù)波形數(shù)據(jù)的存放路徑。為保證程序成功執(zhí)行,在存放連續(xù)波形數(shù)據(jù)的文件夾中需要一個(gè)archive.sta文件,該文件可以通過arccreate程序生成。以上所述兩個(gè)程序和程序使用的說明文檔可以在REFTEK官網(wǎng)下載。arcfetch程序截取連續(xù)波形會(huì)有3種結(jié)果:①截取到完整正確的REFTEK數(shù)據(jù)文件;②沒有截取到REFTEK數(shù)據(jù)文件;③截取到REFTEK數(shù)據(jù)文件,但數(shù)據(jù)為空。對(duì)于后兩種情況需編寫相關(guān)測(cè)試程序,以避免后續(xù)的ref2sac程序出錯(cuò)。

    3.2 ref2sac

    使用arcfetch截取得到REFTEK事件波形,需要使用ref2sac程序轉(zhuǎn)換為SAC格式。該程序?yàn)镻ASSCAL(大陸巖石圈地震臺(tái)陣研究計(jì)劃)提供的rpm軟件包中的一個(gè)子程序,可以將該子程序包在Linux系統(tǒng)下單獨(dú)編譯使用。該程序有兩個(gè)參數(shù),分別是輸入文件夾和輸出文件夾。它將輸入文件夾下REFTEK格式數(shù)據(jù)轉(zhuǎn)換為SAC格式數(shù)據(jù),并存放在輸出文件夾下。當(dāng)該程序運(yùn)行正確時(shí),返回值為0,如果返回值為其他值,表明數(shù)據(jù)格式轉(zhuǎn)換不成功。

    3.3 taup_time

    波形截取時(shí),以理論P(yáng)波到達(dá)臺(tái)站的時(shí)間點(diǎn)為截取波形的參考時(shí)間,選擇截取該參考時(shí)間點(diǎn)前后各多少秒的數(shù)據(jù)。TAUP軟件是目前地震學(xué)中計(jì)算地震波走時(shí)的常用軟件(Crotwell et al,1999)。本程序利用taup_time計(jì)算地震目錄中地震事件傳播到臺(tái)站的理論走時(shí),將理論走時(shí)加上發(fā)震時(shí)刻得到理論P(yáng)波到時(shí),作為波形截取的參考時(shí)間。具體實(shí)現(xiàn)過程中發(fā)現(xiàn),并不是所有震中距都是P波為初至波,實(shí)際上參考時(shí)間應(yīng)該是理論初至波到時(shí)。在使用taup_time程序時(shí),只給定參考模型和震中距參數(shù),比時(shí)TAUP會(huì)計(jì)算出所有震相走時(shí),而且初至波走時(shí)顯示在第1行,本程序提取第1行的震相走時(shí)作為實(shí)際使用的走時(shí)。所以,文中所說P波和初至波是等價(jià)的。

    3.4 saclab程序包

    ref2sac程序只是將REFTEK格式的地震數(shù)據(jù)轉(zhuǎn)換為SAC格式,但是SAC數(shù)據(jù)的頭文件并不正確,所以在轉(zhuǎn)換得到SAC格式地震數(shù)據(jù)后,需要進(jìn)一步修改SAC數(shù)據(jù)的頭文件。saclab是由Andreas Wustefeld使用Matlab語言編寫的程序包,便于SAC文件讀寫和頭文件修改。需要注意的是,O、B、A、E這4個(gè)SAC頭文件,因?yàn)樯婕暗讲ㄐ螖?shù)據(jù)的時(shí)間信息,一旦錯(cuò)誤,后續(xù)利用該SAC波形數(shù)據(jù)進(jìn)行的相關(guān)走時(shí)研究都將會(huì)出現(xiàn)系統(tǒng)性錯(cuò)誤。

    3.5 datestr 和datenum

    Matlab中表示日期的時(shí)間有3種格式:日期字符串、連續(xù)的日期數(shù)值和日期向量(王正林等,2007)。日期字符串格式使用較多,也較為直觀,如2013-10-1 10∶10∶10,但是不方便進(jìn)行時(shí)間的加減運(yùn)算。連續(xù)的日期數(shù)值格式以公元元年1月1日為起點(diǎn)的,用一個(gè)數(shù)值表示當(dāng)前時(shí)間到該起點(diǎn)以天為單位的時(shí)間距離,該數(shù)值為浮點(diǎn)數(shù),所以方便進(jìn)行時(shí)間的加減運(yùn)算。Matlab中的datenum函數(shù)可以實(shí)現(xiàn)將日期字符串和日期向量格式轉(zhuǎn)化為連續(xù)的日期數(shù)字格式,而datestr函數(shù)可以實(shí)現(xiàn)把日期數(shù)字和日期向量格式轉(zhuǎn)換成日期字符串格式,波形截取程序中,需要進(jìn)行儒略日的計(jì)算,利用datenum函數(shù)使用一條語句就可以實(shí)現(xiàn), julianday=datenum(year,month,day)-datenum(year,1,1)+1。波形截取程序需要進(jìn)行時(shí)間的加減運(yùn)算,以確定截取波形的時(shí)間起始點(diǎn),但arcfetch在截取波形時(shí)需要輸入的時(shí)間格式是日期字符串格式,因?yàn)槿掌谧址硎镜臅r(shí)間進(jìn)行加減運(yùn)算比較困難。本程序采取的方法是,在Matlab中,先將日期字符串格式轉(zhuǎn)化為連續(xù)的日期數(shù)值格式,然后進(jìn)行相關(guān)的時(shí)間運(yùn)算,再轉(zhuǎn)換回日期字符串的格式,便于時(shí)間運(yùn)算。

    4 結(jié)束語

    REFTEK連續(xù)波形截取程序,是一個(gè)Matlab函數(shù),功能就是從REFTEK連續(xù)波形數(shù)據(jù)中截取事件波形。該程序連接多個(gè)重要的子程序,并使用相應(yīng)的Matlab工具包,分別完成特定功能,利用Matlab的高端平臺(tái)特性,以簡(jiǎn)單編程高效完成波形截取任務(wù)。在Linux系統(tǒng)下,利用Matlab平臺(tái)連接各種處理程序的編程方式,減少編程時(shí)間,快速完成特定功能,但程序的移植性和易用性較差,不過作為個(gè)人的科研使用仍然具有一定使用價(jià)值。利用該程序截取得到的事件波形,可以進(jìn)行震相到時(shí)讀取,進(jìn)而開展走時(shí)相關(guān)研究工作(黎源等,2012),也可以利用三分量波形數(shù)據(jù)開展包括接收函數(shù)(查小惠等,2013)在內(nèi)的一系列波形研究工作。目前使用該程序截取的波形所開展的相關(guān)研究工作,得到了較好的結(jié)果(Lei et al,2012;黎源等,2012; 查小惠等,2013),進(jìn)一步證明了該程序的可靠性。

    程序編寫得到中國科學(xué)院青藏高原研究所劉啟民博士和中國地震局地殼應(yīng)力研究所孫長(zhǎng)青博士的幫助,在此表示感謝。

    黎源, 雷建設(shè). 青藏高原東緣上地幔頂部Pn波速度結(jié)構(gòu)及各向異性研究[J]. 地球物理學(xué)報(bào), 2012, 55(11):3 615-3 624.

    王正林, 劉明. 精通MATLAB7[M].北京:電子工業(yè)出版社, 2007.

    查小惠, 雷建設(shè). 云南地區(qū)地殼厚度和泊松比研究[J]. 中國科學(xué):地球科學(xué), 2013, 43(3): 446-456.

    Crotwell H P, Owens T J and Ritsema J. The Taup ToolKit: Flexible Seismic Travel-Time and Raypath Utilities[J]. Seismological Research Letters, 1999,70(2):154-160.

    Lei J S, Zhang G W, Xie F R, et al. .Relocation of the 10 March 2011 Yingjiang, China, earthquake sequence and its tectonic implications [J]. Earthquake Science, 2012, 25: 103-110.

    猜你喜歡
    走時(shí)文件夾臺(tái)站
    磁力文件夾
    中國科學(xué)院野外臺(tái)站檔案工作回顧
    氣象基層臺(tái)站建設(shè)
    西藏科技(2021年12期)2022-01-17 08:46:38
    來了晃一圈,走時(shí)已鍍金 有些掛職干部“假裝在基層”
    調(diào)動(dòng)右鍵 解決文件夾管理三大難題
    基層臺(tái)站綜合觀測(cè)業(yè)務(wù)管理之我見
    西藏科技(2015年6期)2015-09-26 12:12:13
    TC一鍵直達(dá)常用文件夾
    電腦迷(2015年1期)2015-04-29 21:24:13
    MDOS平臺(tái)臺(tái)站級(jí)使用方法及技巧
    仰望云天
    意林(2007年20期)2007-05-14 08:14:55
    少妇的逼水好多| 国产白丝娇喘喷水9色精品| 麻豆精品久久久久久蜜桃| 国产在视频线精品| 综合色av麻豆| 床上黄色一级片| 久久久久久大精品| 麻豆国产97在线/欧美| 国产伦精品一区二区三区视频9| 国产女主播在线喷水免费视频网站 | 波野结衣二区三区在线| 岛国在线免费视频观看| 欧美日韩精品成人综合77777| 国产精品久久久久久久久免| 国产 一区精品| 最近最新中文字幕大全电影3| 国产一区二区在线av高清观看| av天堂中文字幕网| 午夜久久久久精精品| 国产成人aa在线观看| 国产高清视频在线观看网站| 天天躁夜夜躁狠狠久久av| 亚洲不卡免费看| 综合色av麻豆| 天堂中文最新版在线下载 | 搡老妇女老女人老熟妇| 成人毛片60女人毛片免费| 啦啦啦观看免费观看视频高清| 国产精品久久久久久久久免| 日韩av在线大香蕉| 天堂√8在线中文| 久久久午夜欧美精品| 欧美三级亚洲精品| 亚洲成色77777| 国产精品1区2区在线观看.| 欧美+日韩+精品| 日韩欧美在线乱码| 国产欧美日韩精品一区二区| 免费黄网站久久成人精品| 久久精品影院6| 男女啪啪激烈高潮av片| 久久久久久久亚洲中文字幕| 久久99热这里只频精品6学生 | 免费观看精品视频网站| 精品久久国产蜜桃| 国产成人福利小说| 国产精品一区二区三区四区免费观看| 成人性生交大片免费视频hd| 变态另类丝袜制服| 全区人妻精品视频| 色视频www国产| 亚洲综合色惰| 成人高潮视频无遮挡免费网站| 国产淫片久久久久久久久| 国产一区亚洲一区在线观看| av天堂中文字幕网| 国产在线一区二区三区精 | 国产片特级美女逼逼视频| 最后的刺客免费高清国语| 国语对白做爰xxxⅹ性视频网站| 午夜福利在线观看吧| 国产精品久久久久久av不卡| 免费搜索国产男女视频| 国产黄片视频在线免费观看| 中文字幕熟女人妻在线| 欧美一区二区亚洲| 中国国产av一级| 最近视频中文字幕2019在线8| 少妇裸体淫交视频免费看高清| av卡一久久| 亚洲av熟女| 国产精品爽爽va在线观看网站| 亚洲精品亚洲一区二区| 在线观看66精品国产| 在线观看66精品国产| 欧美日本视频| 国产片特级美女逼逼视频| 国产一区二区在线观看日韩| 欧美三级亚洲精品| 看十八女毛片水多多多| 天堂av国产一区二区熟女人妻| a级毛色黄片| 国产黄片视频在线免费观看| 久久婷婷人人爽人人干人人爱| 国产亚洲精品av在线| 久久精品夜色国产| 岛国在线免费视频观看| 免费不卡的大黄色大毛片视频在线观看 | 中文天堂在线官网| 岛国毛片在线播放| 亚洲av二区三区四区| h日本视频在线播放| 欧美最新免费一区二区三区| 在线天堂最新版资源| 我的老师免费观看完整版| 男的添女的下面高潮视频| 人人妻人人澡欧美一区二区| 亚洲精品,欧美精品| 村上凉子中文字幕在线| 色噜噜av男人的天堂激情| 热99在线观看视频| 有码 亚洲区| 观看美女的网站| 国产高清有码在线观看视频| 国产男人的电影天堂91| 国产三级中文精品| 大香蕉久久网| 国产在线男女| 老司机福利观看| 国产麻豆成人av免费视频| 变态另类丝袜制服| 神马国产精品三级电影在线观看| 1024手机看黄色片| 欧美日韩精品成人综合77777| 国产高潮美女av| 一个人看视频在线观看www免费| 亚洲国产精品久久男人天堂| 成年免费大片在线观看| 最近手机中文字幕大全| 少妇的逼好多水| 亚洲欧美成人综合另类久久久 | 少妇丰满av| 2022亚洲国产成人精品| 国产毛片a区久久久久| 黄色欧美视频在线观看| 最后的刺客免费高清国语| 69人妻影院| 丝袜美腿在线中文| 国模一区二区三区四区视频| 亚洲四区av| 亚洲一区高清亚洲精品| 国产av一区在线观看免费| 日韩视频在线欧美| 亚洲在久久综合| 精品不卡国产一区二区三区| 人妻夜夜爽99麻豆av| 欧美又色又爽又黄视频| www日本黄色视频网| 变态另类丝袜制服| 一边摸一边抽搐一进一小说| 又爽又黄无遮挡网站| 欧美潮喷喷水| 色网站视频免费| 在线天堂最新版资源| 国产老妇女一区| 午夜视频国产福利| 成人鲁丝片一二三区免费| 亚洲av免费高清在线观看| 国产免费视频播放在线视频 | 男的添女的下面高潮视频| 亚洲丝袜综合中文字幕| 欧美成人午夜免费资源| 69人妻影院| 少妇的逼好多水| 国产精品一及| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 成人高潮视频无遮挡免费网站| 如何舔出高潮| 亚洲欧美一区二区三区国产| 日本免费a在线| 日本猛色少妇xxxxx猛交久久| 日韩制服骚丝袜av| 22中文网久久字幕| 国产精品一区二区三区四区免费观看| 国产高清不卡午夜福利| 成人毛片60女人毛片免费| 欧美日韩在线观看h| 七月丁香在线播放| 亚洲av电影在线观看一区二区三区 | 晚上一个人看的免费电影| videos熟女内射| 欧美+日韩+精品| 天堂中文最新版在线下载 | 国产乱来视频区| 久久精品国产鲁丝片午夜精品| 国产在线男女| 天堂中文最新版在线下载 | 中文欧美无线码| 99热精品在线国产| av在线老鸭窝| 亚洲国产精品国产精品| 桃色一区二区三区在线观看| 中文字幕精品亚洲无线码一区| 建设人人有责人人尽责人人享有的 | 精品99又大又爽又粗少妇毛片| 国产老妇伦熟女老妇高清| 18禁在线播放成人免费| 青青草视频在线视频观看| 床上黄色一级片| 在线观看一区二区三区| 日本wwww免费看| 久久精品久久久久久噜噜老黄 | 日本爱情动作片www.在线观看| 色综合站精品国产| 青青草视频在线视频观看| 日韩av在线大香蕉| videossex国产| 18禁在线无遮挡免费观看视频| 亚洲国产精品成人综合色| 欧美色视频一区免费| 嘟嘟电影网在线观看| 国产淫语在线视频| 欧美日韩一区二区视频在线观看视频在线 | 亚洲av中文av极速乱| 男女下面进入的视频免费午夜| 成年女人永久免费观看视频| 免费看美女性在线毛片视频| 精品一区二区三区视频在线| 午夜久久久久精精品| 麻豆成人av视频| 在线观看av片永久免费下载| 少妇的逼水好多| 欧美成人午夜免费资源| 天堂影院成人在线观看| 午夜久久久久精精品| 国产一区有黄有色的免费视频 | 精品久久国产蜜桃| 有码 亚洲区| 日本爱情动作片www.在线观看| 一级毛片我不卡| videossex国产| 欧美xxxx黑人xx丫x性爽| 丝袜美腿在线中文| 丝袜美腿在线中文| 18禁在线无遮挡免费观看视频| 91久久精品国产一区二区成人| 一级毛片我不卡| 精品一区二区三区视频在线| 性插视频无遮挡在线免费观看| 日韩av在线大香蕉| 午夜免费激情av| 免费观看的影片在线观看| 亚洲综合精品二区| 亚洲人与动物交配视频| 国产乱人偷精品视频| 一级毛片aaaaaa免费看小| 国产免费男女视频| 91狼人影院| 天堂av国产一区二区熟女人妻| 在线观看66精品国产| 欧美又色又爽又黄视频| 久久精品影院6| 日本一二三区视频观看| 精品久久久久久成人av| 91狼人影院| 身体一侧抽搐| 又粗又硬又长又爽又黄的视频| 最近中文字幕高清免费大全6| 狂野欧美白嫩少妇大欣赏| 精品久久久久久久久久久久久| 岛国毛片在线播放| 久久精品人妻少妇| 国产av不卡久久| videossex国产| 成年av动漫网址| 特级一级黄色大片| 亚洲精品,欧美精品| 欧美性感艳星| 狂野欧美白嫩少妇大欣赏| 国产亚洲av片在线观看秒播厂 | av在线天堂中文字幕| 天美传媒精品一区二区| 亚洲欧美日韩卡通动漫| 一级毛片我不卡| 亚洲国产色片| 女人被狂操c到高潮| 亚洲人成网站高清观看| 国产伦一二天堂av在线观看| 国产成人一区二区在线| 黄色一级大片看看| 丰满人妻一区二区三区视频av| 看片在线看免费视频| 久久久久久久久中文| 最近中文字幕2019免费版| 日本色播在线视频| 国产精品国产三级国产av玫瑰| 亚洲国产欧美人成| 一二三四中文在线观看免费高清| 精品一区二区三区视频在线| 国产在视频线在精品| 在线播放无遮挡| 22中文网久久字幕| 日韩大片免费观看网站 | 建设人人有责人人尽责人人享有的 | 最后的刺客免费高清国语| 婷婷色av中文字幕| 国产伦一二天堂av在线观看| 亚洲精品自拍成人| 色哟哟·www| 熟女人妻精品中文字幕| 久久精品熟女亚洲av麻豆精品 | 国产精品不卡视频一区二区| 欧美一区二区精品小视频在线| 日韩三级伦理在线观看| 国产精品久久久久久av不卡| 亚洲国产精品成人久久小说| 亚洲av电影不卡..在线观看| 看片在线看免费视频| 国产av码专区亚洲av| 成人三级黄色视频| 精品国产三级普通话版| 亚洲自偷自拍三级| 在线观看美女被高潮喷水网站| 国产中年淑女户外野战色| 十八禁国产超污无遮挡网站| 久久6这里有精品| 国产麻豆成人av免费视频| 可以在线观看毛片的网站| 国产免费一级a男人的天堂| 美女黄网站色视频| 18禁裸乳无遮挡免费网站照片| 观看免费一级毛片| 日本免费在线观看一区| av免费在线看不卡| 久久久色成人| 久久久精品大字幕| 天堂中文最新版在线下载 | 国产人妻一区二区三区在| 晚上一个人看的免费电影| 色5月婷婷丁香| 欧美激情久久久久久爽电影| 男人舔奶头视频| 日本wwww免费看| 日日摸夜夜添夜夜爱| 国产69精品久久久久777片| 亚洲欧美成人精品一区二区| 日本与韩国留学比较| 久久99热这里只频精品6学生 | 成人欧美大片| 天天躁夜夜躁狠狠久久av| 国产一区亚洲一区在线观看| 老司机影院成人| 99热精品在线国产| 插逼视频在线观看| 精品午夜福利在线看| 精品久久久久久久末码| 亚洲精品乱码久久久v下载方式| 欧美激情国产日韩精品一区| 久久精品人妻少妇| 最近最新中文字幕免费大全7| 男女视频在线观看网站免费| 蜜桃亚洲精品一区二区三区| 一个人看的www免费观看视频| 最近最新中文字幕免费大全7| 亚洲va在线va天堂va国产| 成人美女网站在线观看视频| 午夜福利高清视频| 精品久久久久久久久av| 一个人观看的视频www高清免费观看| 亚洲人与动物交配视频| 日韩成人伦理影院| 国产三级中文精品| 舔av片在线| 日本欧美国产在线视频| 国产精品国产三级国产专区5o | 男女下面进入的视频免费午夜| 国产高清有码在线观看视频| 久久精品国产亚洲网站| 大又大粗又爽又黄少妇毛片口| 白带黄色成豆腐渣| 国产精品人妻久久久影院| 搡老妇女老女人老熟妇| 嘟嘟电影网在线观看| 午夜视频国产福利| 2022亚洲国产成人精品| 中国国产av一级| 亚洲欧美日韩高清专用| 男人狂女人下面高潮的视频| 亚洲av成人精品一区久久| 中文亚洲av片在线观看爽| 亚洲精品国产成人久久av| 麻豆av噜噜一区二区三区| 人妻制服诱惑在线中文字幕| 亚洲av.av天堂| 午夜激情欧美在线| 国产探花极品一区二区| 在线观看一区二区三区| 久久午夜福利片| 日本一本二区三区精品| 日韩一本色道免费dvd| 国产精品久久久久久av不卡| 免费无遮挡裸体视频| 日韩成人av中文字幕在线观看| 久久久欧美国产精品| 久久久久久国产a免费观看| 夜夜爽夜夜爽视频| 国产av不卡久久| 成人美女网站在线观看视频| 中文精品一卡2卡3卡4更新| 国产高清三级在线| 免费不卡的大黄色大毛片视频在线观看 | 国产老妇伦熟女老妇高清| 九九在线视频观看精品| 日韩欧美精品免费久久| 最近中文字幕2019免费版| 精品久久久久久久末码| 黄色日韩在线| 国产免费视频播放在线视频 | 色噜噜av男人的天堂激情| 久热久热在线精品观看| 五月伊人婷婷丁香| 欧美另类亚洲清纯唯美| 国产精品av视频在线免费观看| 欧美另类亚洲清纯唯美| 日本wwww免费看| 亚洲国产精品久久男人天堂| 精品一区二区免费观看| 国产亚洲av片在线观看秒播厂 | 欧美一区二区国产精品久久精品| 国产午夜精品一二区理论片| 青春草亚洲视频在线观看| 97超碰精品成人国产| 噜噜噜噜噜久久久久久91| 小蜜桃在线观看免费完整版高清| 舔av片在线| 国产亚洲精品av在线| 亚洲av电影在线观看一区二区三区 | 老司机影院成人| 国语自产精品视频在线第100页| 久久久久久久午夜电影| 欧美丝袜亚洲另类| 亚洲性久久影院| 成人三级黄色视频| 中国国产av一级| 国产精品久久久久久精品电影小说 | 精品国产三级普通话版| 18禁裸乳无遮挡免费网站照片| 国产精品永久免费网站| 亚洲国产高清在线一区二区三| 在线免费观看的www视频| av国产久精品久网站免费入址| 久久久久久久久久久丰满| 欧美日韩综合久久久久久| 美女国产视频在线观看| 国产老妇伦熟女老妇高清| 久久久欧美国产精品| 国产一区二区三区av在线| 午夜久久久久精精品| 网址你懂的国产日韩在线| 热99在线观看视频| 毛片女人毛片| 国产精品麻豆人妻色哟哟久久 | 欧美日本视频| 亚洲激情五月婷婷啪啪| 淫秽高清视频在线观看| 欧美3d第一页| 国产精品一区www在线观看| 亚洲人与动物交配视频| .国产精品久久| 久久综合国产亚洲精品| 成人综合一区亚洲| 最新中文字幕久久久久| 三级毛片av免费| 中文亚洲av片在线观看爽| 嫩草影院新地址| 天堂av国产一区二区熟女人妻| 国产成人免费观看mmmm| 成人鲁丝片一二三区免费| av在线蜜桃| 欧美激情在线99| 午夜精品一区二区三区免费看| 久久99蜜桃精品久久| 毛片一级片免费看久久久久| 中文字幕人妻熟人妻熟丝袜美| 小说图片视频综合网站| 国产欧美另类精品又又久久亚洲欧美| 成人性生交大片免费视频hd| 国产伦精品一区二区三区四那| 婷婷色麻豆天堂久久 | 国产在线一区二区三区精 | eeuss影院久久| 大香蕉97超碰在线| 99热精品在线国产| 国产爱豆传媒在线观看| av在线天堂中文字幕| 别揉我奶头 嗯啊视频| 成年女人看的毛片在线观看| 欧美性猛交黑人性爽| 亚洲精品456在线播放app| 亚洲色图av天堂| 久久人人爽人人爽人人片va| 建设人人有责人人尽责人人享有的 | 麻豆国产97在线/欧美| 免费看日本二区| 人体艺术视频欧美日本| av在线天堂中文字幕| 亚洲国产精品国产精品| 日韩人妻高清精品专区| 日本爱情动作片www.在线观看| 久久久久久久亚洲中文字幕| 日日啪夜夜撸| 视频中文字幕在线观看| 国产在线一区二区三区精 | 热99在线观看视频| 老司机影院成人| 国产精品1区2区在线观看.| 99九九线精品视频在线观看视频| 国产一区有黄有色的免费视频 | av天堂中文字幕网| 嫩草影院新地址| 亚洲精华国产精华液的使用体验| АⅤ资源中文在线天堂| 国产综合懂色| 亚洲怡红院男人天堂| 亚洲欧美精品综合久久99| 搡女人真爽免费视频火全软件| 又黄又爽又刺激的免费视频.| 国产精品日韩av在线免费观看| 亚洲无线观看免费| 亚洲aⅴ乱码一区二区在线播放| 熟女电影av网| 午夜激情欧美在线| 在线观看一区二区三区| 欧美一级a爱片免费观看看| 狂野欧美激情性xxxx在线观看| 尾随美女入室| 性插视频无遮挡在线免费观看| 亚洲中文字幕一区二区三区有码在线看| 乱系列少妇在线播放| 深夜a级毛片| 精品久久久久久久久久久久久| 一边亲一边摸免费视频| 久久久久久久久久成人| 岛国在线免费视频观看| 18+在线观看网站| 精品久久久久久久久av| 真实男女啪啪啪动态图| 日本午夜av视频| 中文字幕精品亚洲无线码一区| 人体艺术视频欧美日本| 欧美性感艳星| 国产成人a区在线观看| 欧美日韩精品成人综合77777| 久久久色成人| 搞女人的毛片| a级毛片免费高清观看在线播放| 中文在线观看免费www的网站| 深夜a级毛片| 能在线免费观看的黄片| av.在线天堂| 久久精品久久精品一区二区三区| 秋霞伦理黄片| 久久久久久久久久黄片| 国产在视频线在精品| 在线观看一区二区三区| 成人鲁丝片一二三区免费| 亚洲精品乱码久久久久久按摩| 国产精品美女特级片免费视频播放器| 麻豆久久精品国产亚洲av| 毛片女人毛片| 永久网站在线| 日韩亚洲欧美综合| 日本-黄色视频高清免费观看| 亚洲精品国产av成人精品| 色网站视频免费| 亚洲欧美中文字幕日韩二区| 亚洲精品成人久久久久久| 日本一本二区三区精品| 欧美一区二区亚洲| 欧美不卡视频在线免费观看| 欧美成人a在线观看| 欧美最新免费一区二区三区| 免费搜索国产男女视频| 日日摸夜夜添夜夜添av毛片| 九色成人免费人妻av| 老师上课跳d突然被开到最大视频| 91aial.com中文字幕在线观看| 亚洲欧美日韩高清专用| 国产白丝娇喘喷水9色精品| 综合色av麻豆| 久久精品熟女亚洲av麻豆精品 | kizo精华| 国产成人福利小说| 久久久久久久亚洲中文字幕| 最近的中文字幕免费完整| 91精品一卡2卡3卡4卡| 97超视频在线观看视频| 欧美成人午夜免费资源| 精品人妻视频免费看| 秋霞在线观看毛片| 纵有疾风起免费观看全集完整版 | 免费看a级黄色片| 麻豆成人午夜福利视频| 免费观看在线日韩| 一区二区三区免费毛片| 国产成人freesex在线| 欧美性猛交黑人性爽| 国产色婷婷99| 国产精品av视频在线免费观看| 麻豆av噜噜一区二区三区| 久久久亚洲精品成人影院| 久久国内精品自在自线图片| 波野结衣二区三区在线| 免费看美女性在线毛片视频| 99热精品在线国产| 亚洲一级一片aⅴ在线观看| 99热这里只有精品一区| 亚洲成人中文字幕在线播放| 亚洲一区高清亚洲精品| 91午夜精品亚洲一区二区三区| 午夜爱爱视频在线播放| 日韩欧美 国产精品| 男人和女人高潮做爰伦理| 网址你懂的国产日韩在线| 久久6这里有精品| 天天一区二区日本电影三级| 日韩成人av中文字幕在线观看| 国内精品一区二区在线观看| 日韩av在线大香蕉| 在线a可以看的网站| 亚洲va在线va天堂va国产| 亚洲欧美精品自产自拍| 精品人妻熟女av久视频| 又粗又爽又猛毛片免费看| 国产精品一二三区在线看|