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

    海拉爾臺(tái)陣數(shù)據(jù)功率譜分析

    2020-07-01 09:02:38白永福奚沖霄
    關(guān)鍵詞:海拉爾譜估計(jì)背景噪聲

    白永福,高 翔,張 偉,奚沖霄,王 凌

    (北京市地震局,北京 100080)

    0 引言

    國(guó)際監(jiān)測(cè)系統(tǒng)(IMS)常用4 種監(jiān)測(cè)方法(地震監(jiān)測(cè)、水聲監(jiān)測(cè)、次聲監(jiān)測(cè)、放射性核素監(jiān)測(cè))對(duì)核試驗(yàn)進(jìn)行監(jiān)測(cè)[1],其中地震監(jiān)測(cè)方法采用全球50 個(gè)地震臺(tái)陣(站),這其中就包括我國(guó)的海拉爾臺(tái)陣(PS12)和蘭州臺(tái)陣(PS13)。臺(tái)陣研究最早的目的是為了地下核試驗(yàn)地震信號(hào)的捕獲,檢測(cè)定位低震級(jí)事件,提高信號(hào)檢測(cè)能力。在國(guó)內(nèi),郝春月、鄭重、王娟等在地震臺(tái)陣方面做了很多研究[1,8-9],對(duì)MS型小孔徑臺(tái)陣的臺(tái)陣設(shè)計(jì)、勘址數(shù)據(jù)分析處理、臺(tái)陣技術(shù)等進(jìn)行了一系列研究[13]。隨著臺(tái)陣技術(shù)的逐漸成熟,地震臺(tái)陣不僅要是完成最初目的,而且現(xiàn)在也已經(jīng)為地震學(xué)和地球物理學(xué)的發(fā)展做出了重要的貢獻(xiàn),地震臺(tái)陣可通過(guò)壓低干擾背景,疊加增強(qiáng)地震信號(hào),對(duì)于弱信號(hào)的識(shí)別有重要作用,有助于獲取地下速度結(jié)構(gòu)、提高定位精度、了解有關(guān)震源及地球內(nèi)部結(jié)構(gòu)信息[2]。從上地幔到地球內(nèi)核,從地震震源到接近地表接收臺(tái)站的結(jié)構(gòu),地震臺(tái)陣都有方法進(jìn)行探究。地震臺(tái)陣的主要技術(shù)有聚束、速度譜分析、f-k分析等[3、5]。

    數(shù)據(jù)質(zhì)量控制是所有地震學(xué)、地球物理學(xué)基礎(chǔ)研究;所有地震預(yù)警、海嘯預(yù)警、地震預(yù)測(cè)、地震監(jiān)測(cè)、核爆炸監(jiān)測(cè)、臺(tái)站、臺(tái)網(wǎng)、臺(tái)陣建設(shè)等工作,都依靠地震數(shù)據(jù),數(shù)據(jù)質(zhì)量不好直接影響工作質(zhì)量,所以數(shù)據(jù)質(zhì)量控制是國(guó)際各大臺(tái)網(wǎng)都做的必要工作。數(shù)據(jù)質(zhì)量參數(shù)主要包括,背景噪聲水平、數(shù)據(jù)接收率、數(shù)據(jù)可用率、數(shù)據(jù)延時(shí)平均數(shù)、Gap 間隔、Gap 長(zhǎng)度,Spike 間隔、Spike 振幅、均方根等,其中背景噪聲水平是最基本的參數(shù),噪聲水平直接影響臺(tái)站的監(jiān)測(cè)能力,所以對(duì)噪聲的持續(xù)性分析對(duì)評(píng)價(jià)整個(gè)臺(tái)站具有重要的意義,如果噪聲水平達(dá)不到監(jiān)測(cè)要求,則存在遷址的可能。在本文中,主要對(duì)海拉爾臺(tái)陣的背景噪聲進(jìn)行了分析和研究。

    1 海拉爾臺(tái)陣概況

    海拉爾臺(tái)陣位于內(nèi)蒙古呼倫貝爾陳巴爾虎旗,距離海拉爾市38 公里。海拉爾區(qū)位于海拉爾構(gòu)造盆地的東段,在東北震區(qū)的兩個(gè)邊緣,其特點(diǎn)從屬東北震區(qū)[4]。地質(zhì)上主要是晚古生代的褶皺帶——大興安嶺褶皺帶,松遼平原主要是中生代的構(gòu)造盆地。海拉爾臺(tái)陣投入了九套地震儀,部署在中心臺(tái)附近(中心臺(tái)為三分項(xiàng)),為小孔徑臺(tái)陣(3 ~5km)。臺(tái)陣采用CMG-ESPV 型地震儀,Nanometrics 公司生產(chǎn)的Europa 型數(shù)采,在0.02~16Hz 頻帶范圍內(nèi)呈加速度平坦型幅頻特性。9 個(gè)子臺(tái)位于40~50m 的鉆孔中,其分布如圖1 所示。海拉爾日均溫差10 度左右,附近廠礦較多。

    圖1 臺(tái)站分布Fig.1 Distribution of the seismic stations

    海拉爾臺(tái)陣自2001 年建設(shè)完成后,由于各種原因一直無(wú)法實(shí)時(shí)向IDC(國(guó)際數(shù)據(jù)中心)傳送數(shù)據(jù),2004 年海拉爾臺(tái)陣通過(guò)我國(guó)驗(yàn)收,開(kāi)始試運(yùn)行,并于2007 年7 月向禁核試北京國(guó)家數(shù)據(jù)中心傳送實(shí)時(shí)數(shù)據(jù)[12]。2013 年12 月地球物理研究所與CTBTO (全面禁止核試驗(yàn)條約組織)簽訂合同,開(kāi)始執(zhí)行檢測(cè)與評(píng)估(Test&Estimation)計(jì)劃,實(shí)時(shí)向IDC 傳送數(shù)據(jù)[1]。

    2 資料處理

    2.1 數(shù)據(jù)收集與整理

    本文使用數(shù)據(jù)來(lái)源于中國(guó)地震局地球物理研究所,由于臺(tái)陣建立較早,臺(tái)陣設(shè)備數(shù)據(jù)采用wfdisk 格式進(jìn)行存儲(chǔ),為了便于文件后期使用和進(jìn)行更多研究,通過(guò)Python 語(yǔ)言將wfdisk文件轉(zhuǎn)換成sac 格式文件。由于缺乏管理與運(yùn)行維護(hù),自2001 年到2013 年間的數(shù)據(jù)質(zhì)量存在很多問(wèn)題。作為我國(guó)僅有的IMS 小孔徑臺(tái)陣歷史數(shù)據(jù),對(duì)它們進(jìn)行數(shù)據(jù)質(zhì)量評(píng)估是非常必要的。

    2.2 功率譜基本原理

    經(jīng)典譜估計(jì)方法在應(yīng)用中都是以離散傅里葉變換為近處的,它有著分辨率不高、能量泄露的缺點(diǎn),為了克服這些缺點(diǎn),經(jīng)過(guò)研究提出了各種算法。目前經(jīng)典譜估計(jì)計(jì)算方法有周期圖法、Barlett 算法、Welch 算法、Nattall 算法等。Welch 算法是有由Welch 提出的對(duì)周期圖的修正算法,是經(jīng)典譜估計(jì)中過(guò)的有效應(yīng)用的一種算法。功率譜密度的計(jì)算采用經(jīng)典譜估計(jì)中的Welch 算法,Welch 算法是在周期圖法的基礎(chǔ)上改進(jìn)得到的一種修正周期圖功率譜密度估計(jì)方法。Welch 算法譜估計(jì)采用數(shù)據(jù)分段加窗再求平均的方法,先分別求出每段的譜估計(jì),然后進(jìn)行總平均。即:把某一長(zhǎng)度為N 的數(shù)據(jù)x(n)分成L 段,每一段為M(在分段時(shí)可允許每一段的數(shù)據(jù)有部分的重疊),分別求每段的功率譜,然后加以平均,每一段的數(shù)據(jù)窗口可使用哈寧窗或者哈明窗。這樣可以改善由于矩形窗邊瓣較大所產(chǎn)生的譜失真。用P(ω)表示用Welch 平均周期圖法算法估計(jì)的功率譜。即:

    公式(1)中:N 為數(shù)據(jù)噪聲長(zhǎng)度

    L 為數(shù)據(jù)N 的分段數(shù)

    U 為歸一化因子

    d(n)是數(shù)據(jù)窗

    e-jam為傅里葉變換

    2.3 計(jì)算方法

    對(duì)海拉爾臺(tái)陣的中心臺(tái)進(jìn)行計(jì)算,每天選出4 個(gè)時(shí)間點(diǎn),00/06/12/18 點(diǎn)進(jìn)行計(jì)算,每個(gè)時(shí)間點(diǎn)計(jì)算一個(gè)小時(shí)的功率譜,如果存在地震,則避開(kāi)地震,利用前后1 小時(shí)代替。應(yīng)用Matlab 軟件進(jìn)行處理,本文對(duì)每小時(shí)40 采樣率共144000 點(diǎn)進(jìn)行計(jì)算繪圖,并把計(jì)算結(jié)果和圖形進(jìn)行存儲(chǔ)。數(shù)據(jù)處理過(guò)程如下:

    首先對(duì)轉(zhuǎn)換后的數(shù)據(jù)進(jìn)行波形讀取,無(wú)地震、連續(xù)性高且干擾不嚴(yán)重的時(shí)段。

    對(duì)數(shù)據(jù)進(jìn)行去均值處理,得到速度記錄(m/s),根據(jù)海拉爾臺(tái)陣的數(shù)采采用Nanometrics公司生產(chǎn)的Europa 型數(shù)采,在0.02~16Hz 頻帶范圍內(nèi)呈加速度平坦型幅頻特性,確定轉(zhuǎn)換因子。

    對(duì)每段數(shù)據(jù)利用Welch 平均周期圖法進(jìn)行功率譜計(jì)算。

    將功率譜轉(zhuǎn)換成加速度譜,根據(jù)地震計(jì)的傳遞函數(shù)(海拉爾臺(tái)陣采用CMG-ESPV 型地震儀),去除儀器影響,得到加速度功率譜。

    為了能夠讓結(jié)果客觀的反映海拉爾臺(tái)陣的背景噪聲水平,本文選取了海拉爾臺(tái)陣2002 年1—4 月的數(shù)據(jù)進(jìn)行了計(jì)算,作為數(shù)據(jù)樣本,是因?yàn)橄鄬?duì)其他時(shí)間此時(shí)間段比較平靜且中斷率低。

    圖2 2002 年某日功率譜圖Fig.2 The power spectrum of one day in 2002

    3 噪聲功率譜分析

    發(fā)生地震時(shí),功率譜密度與平靜狀態(tài)下的噪聲有顯著不同,長(zhǎng)周期的噪聲不容忽視,受溫度與氣壓的影響較大。除了地震活動(dòng)和周期影響外,人為因素、儀器故障、維修等非噪聲事件也會(huì)引起功率譜異常。根據(jù)功率譜圖(圖3)可以看出1 到10Hz 區(qū)域比較穩(wěn)定,噪聲水平低,在1Hz、5Hz 和10Hz 頻段噪聲水平平均達(dá)到了-150dB、-144dB、-155dB。

    3.1 噪聲隨時(shí)間的變化

    為了更好地研究周期性變化和異常,對(duì)得到的每日功率譜值進(jìn)行繪制0.02Hz、0.2Hz、1Hz、2.5Hz、6.5Hz 頻段的單頻曲線。對(duì)比白天和夜間各兩個(gè)時(shí)間段的功率譜進(jìn)行分析(圖2),發(fā)現(xiàn)高頻曲線2.5Hz 頻段中白天時(shí)段的功率譜密度值高于夜間的功率譜密度值(低值相差5dB,高值相差7dB),而且形態(tài)相似。由于樣本數(shù)據(jù)時(shí)段是冬季,海拉爾的氣溫屬于冬季嚴(yán)寒、漫長(zhǎng),地面積雪時(shí)間長(zhǎng),1 月(最冷月)平均低溫為-30.83℃,極寒天氣主要出現(xiàn)在冬季,高發(fā)期在12 月下旬至1 月下旬,2 月份出現(xiàn)極寒天氣也較多,所以極寒天氣影響人們的出行和各類(lèi)工作活動(dòng)。海拉爾極寒天氣多,在全市出現(xiàn)的極寒天氣中,21—23 時(shí)是極寒天氣開(kāi)始增多的時(shí)段,11—12 時(shí)和16—20 時(shí)氣溫相對(duì)溫暖。由于氣溫原因,11—12 點(diǎn)16—20點(diǎn)也成為冬季出行、工作、工廠、活動(dòng)等相對(duì)頻繁時(shí)段。這也符合人的作息規(guī)律,人為噪聲是來(lái)自各個(gè)方面的波疊加而形成相對(duì)穩(wěn)定的。

    圖3 白天和夜間功率譜Fig.3 The power spectrum for daytime and nighttime

    通過(guò)單頻曲線發(fā)現(xiàn)各個(gè)時(shí)間段在曲線形態(tài)變化上一致性很高,但發(fā)現(xiàn)在2 月1 日至2 月27 日間,各頻率噪聲相對(duì)于別的時(shí)段變化較大(圖4),具體如下:2.5Hz 背景噪聲分別在0點(diǎn)、6 點(diǎn)、12 點(diǎn)、18 點(diǎn)時(shí)出現(xiàn)噪聲降低的現(xiàn)象,0 時(shí)存在10dB 差值;6 時(shí)存在12dB 差值;12 時(shí)存在10dB 差值;18 時(shí)存在18dB 差值。因?yàn)?002 年2 月正是春節(jié)期間,工廠和單位停工、人文活動(dòng)減少,表示高頻能量減少,導(dǎo)致在高頻的變化曲線上,2 月份出現(xiàn)低值。就季節(jié)變化而言,具體還需要增大數(shù)據(jù)量進(jìn)行分析對(duì)比,比如用同年其他時(shí)間段與1-4 月對(duì)比,再用往年同時(shí)間段進(jìn)行對(duì)比。

    3.2 噪聲源分析

    通過(guò)功率譜圖發(fā)現(xiàn)在2002 年1—4 月期間每日12 時(shí)和18 時(shí)都在4Hz 附近處出現(xiàn)大約10dB 的單頻率脈沖(圖5),如圖5(a)中橢圓部分是2002 年1 月21 日12 時(shí)在4Hz 處發(fā)現(xiàn)的單頻率脈沖。

    圖5 2002 年1 月至4 月單頻率脈沖Fig.5 The single frequency pulse from January to April in 2002

    陣址噪聲源會(huì)對(duì)地震臺(tái)陣的監(jiān)測(cè)造成影響,但是統(tǒng)計(jì)好噪聲源的頻率、時(shí)間、大小,可作為地震監(jiān)測(cè)背景,對(duì)監(jiān)測(cè)工作提供重要的參考依據(jù)[11]。海拉爾附近有很多大的礦區(qū)(圖6),寶日希勒煤礦距臺(tái)陣約20km;伊敏煤礦距臺(tái)陣約75km;大雁煤礦距臺(tái)陣約70km;免渡河煤礦、五九煤礦距臺(tái)陣約120km;扎拉諾爾煤礦距里臺(tái)陣大約178km。

    經(jīng)考察,脈沖主要與寶日希勒煤礦有關(guān),此煤礦為海拉爾較大的煤礦,距離海拉爾臺(tái)陣近,經(jīng)調(diào)查,該煤礦每天中午12 點(diǎn)和晚上18點(diǎn)左右進(jìn)行人工爆炸,主要用于開(kāi)采煤礦。

    圖6 海拉爾附近礦區(qū)Fig.6 The mining area near Hailar

    4 結(jié)論與討論

    海拉爾臺(tái)陣位于陳巴爾虎旗,距離海拉爾市38km。海拉爾區(qū)位于海拉爾構(gòu)造盆地的東段,在東北震區(qū)的兩個(gè)邊緣,其特點(diǎn)從屬東北震區(qū)[12]。地質(zhì)上主要是晚古生代的褶皺帶——大興安嶺褶皺帶,松遼平原主要是中生代的構(gòu)造盆地。通過(guò)海拉爾臺(tái)陣噪聲背景分析,從時(shí)間上看,對(duì)夜間和白天的噪聲功率譜計(jì)算結(jié)果表明,背景噪聲水平?jīng)]有明顯的變化,從一致性上看,人為噪聲基本不影響背景噪聲水平的變化。海拉爾臺(tái)陣噪聲譜均處于全球噪聲的標(biāo)準(zhǔn)曲線之間,且明顯偏向于低噪聲一側(cè)曲線,平均低于噪聲曲線高值30dB,說(shuō)明觀測(cè)點(diǎn)噪聲水平低,臺(tái)陣背景噪聲符合并優(yōu)于地震觀測(cè)對(duì)背景噪聲的要求。

    下一步有待繼續(xù)積累資料進(jìn)行研究,通過(guò)增大臺(tái)陣內(nèi)其他觀測(cè)點(diǎn)的數(shù)據(jù)量,對(duì)臺(tái)陣整體和各子臺(tái)間對(duì)比等方面進(jìn)行分析,從而達(dá)到對(duì)該臺(tái)陣背景噪聲的進(jìn)一步跟蹤和研究。

    致謝:在本文的撰寫(xiě)過(guò)程中,獲得中國(guó)地震局地球物理研究所郝春月副研究員的支持,在此表示感謝。

    猜你喜歡
    海拉爾譜估計(jì)背景噪聲
    窄帶電力線通信信道背景噪聲抑制方法
    海拉爾油田低產(chǎn)低效區(qū)塊引效壓裂治理技術(shù)
    應(yīng)用背景噪聲成像研究祁連山地區(qū)地殼S波速度結(jié)構(gòu)
    地震研究(2017年3期)2017-11-06 23:38:05
    基于MATLAB實(shí)現(xiàn)的AR模型功率譜估計(jì)
    “海拉爾杯”草原星第五屆內(nèi)蒙古青年歌手電視大獎(jiǎng)賽完美落幕
    草原歌聲(2017年4期)2017-04-28 08:20:48
    “海拉爾杯”草原星第五屆內(nèi)蒙古青年歌手電視大賽圓滿落下帷幕
    草原歌聲(2017年4期)2017-04-28 08:20:39
    海上單道地震勘探中船舶等背景噪聲的影響分析及壓制
    經(jīng)典功率譜估計(jì)方法的研究
    Welch譜估計(jì)的隨機(jī)誤差與置信度
    青藏高原東北緣背景噪聲特征分析
    地震研究(2014年4期)2014-02-27 09:31:15
    成人av一区二区三区在线看| 亚洲五月天丁香| 99精品久久久久人妻精品| 自拍欧美九色日韩亚洲蝌蚪91| 欧美成人性av电影在线观看| 国产成人欧美在线观看| av电影中文网址| 国产高清视频在线播放一区| 亚洲精品国产精品久久久不卡| а√天堂www在线а√下载| 极品人妻少妇av视频| 午夜免费鲁丝| 老汉色∧v一级毛片| 午夜免费观看网址| 禁无遮挡网站| 黄片播放在线免费| 黄色 视频免费看| 黄色女人牲交| 电影成人av| 成熟少妇高潮喷水视频| 9热在线视频观看99| 亚洲国产毛片av蜜桃av| 国产成人精品久久二区二区91| 精品熟女少妇八av免费久了| 精品久久久久久久人妻蜜臀av | 久久精品国产亚洲av香蕉五月| 视频区欧美日本亚洲| 国产熟女午夜一区二区三区| 日日夜夜操网爽| 在线观看www视频免费| 久久天堂一区二区三区四区| 亚洲五月婷婷丁香| 人妻丰满熟妇av一区二区三区| 国产精品99久久99久久久不卡| 人人妻人人澡欧美一区二区 | 成人永久免费在线观看视频| 亚洲国产精品999在线| 久久人妻福利社区极品人妻图片| 97超级碰碰碰精品色视频在线观看| 精品久久久久久久久久免费视频| 国产精品亚洲一级av第二区| www.999成人在线观看| 成人手机av| 国产精品一区二区在线不卡| 黄网站色视频无遮挡免费观看| 日韩精品免费视频一区二区三区| 大型av网站在线播放| 可以免费在线观看a视频的电影网站| 看黄色毛片网站| 亚洲国产欧美日韩在线播放| 成人免费观看视频高清| 最新在线观看一区二区三区| 国产熟女xx| av天堂在线播放| 一个人免费在线观看的高清视频| 热99re8久久精品国产| 一级毛片精品| 最好的美女福利视频网| 免费女性裸体啪啪无遮挡网站| 精品一区二区三区四区五区乱码| 国产1区2区3区精品| 99久久久亚洲精品蜜臀av| www.999成人在线观看| 国产在线精品亚洲第一网站| 久热爱精品视频在线9| 啦啦啦免费观看视频1| 久久久国产精品麻豆| 女生性感内裤真人,穿戴方法视频| 久久人人爽av亚洲精品天堂| 9色porny在线观看| 天天一区二区日本电影三级 | 免费av毛片视频| 成人免费观看视频高清| 久久亚洲真实| 欧美一级毛片孕妇| 我的亚洲天堂| 国产精品98久久久久久宅男小说| 一级a爱视频在线免费观看| 亚洲午夜精品一区,二区,三区| 欧美大码av| 精品高清国产在线一区| 中文字幕精品免费在线观看视频| 丝袜美足系列| 欧美激情久久久久久爽电影 | 国产高清有码在线观看视频 | 757午夜福利合集在线观看| 国产亚洲精品久久久久久毛片| 国产精品免费视频内射| 国产人伦9x9x在线观看| 久久人妻av系列| 欧美绝顶高潮抽搐喷水| 亚洲欧美日韩另类电影网站| 国产高清视频在线播放一区| 老司机靠b影院| 国产免费av片在线观看野外av| 国产欧美日韩一区二区三| 极品教师在线免费播放| 日韩高清综合在线| 欧美最黄视频在线播放免费| 色在线成人网| 国产成人精品无人区| 亚洲国产中文字幕在线视频| 欧美亚洲日本最大视频资源| 久久久国产精品麻豆| 国产亚洲欧美98| 少妇熟女aⅴ在线视频| 亚洲片人在线观看| 中文字幕人妻熟女乱码| 老司机深夜福利视频在线观看| 午夜福利18| 亚洲aⅴ乱码一区二区在线播放 | 夜夜爽天天搞| 啦啦啦观看免费观看视频高清 | 最新在线观看一区二区三区| 亚洲欧美日韩另类电影网站| 成人免费观看视频高清| 亚洲精品国产一区二区精华液| 国产精品乱码一区二三区的特点 | 91大片在线观看| 丝袜在线中文字幕| 制服丝袜大香蕉在线| 校园春色视频在线观看| 成人av一区二区三区在线看| 国产av精品麻豆| 老司机福利观看| 精品少妇一区二区三区视频日本电影| 天堂动漫精品| 非洲黑人性xxxx精品又粗又长| 中文字幕av电影在线播放| 国产精品一区二区精品视频观看| 香蕉国产在线看| 90打野战视频偷拍视频| 啦啦啦 在线观看视频| 高清毛片免费观看视频网站| 国产成人一区二区三区免费视频网站| 国内精品久久久久久久电影| 国产精品亚洲美女久久久| 美女大奶头视频| 麻豆av在线久日| 我的亚洲天堂| 国产真人三级小视频在线观看| 亚洲aⅴ乱码一区二区在线播放 | 午夜老司机福利片| 国产精品一区二区精品视频观看| 99在线人妻在线中文字幕| 国产亚洲欧美98| 免费无遮挡裸体视频| 国产精品二区激情视频| 国产成人欧美在线观看| 亚洲欧美激情在线| 欧美激情高清一区二区三区| 国产日韩一区二区三区精品不卡| www.自偷自拍.com| 1024视频免费在线观看| 日日夜夜操网爽| 无遮挡黄片免费观看| 午夜福利在线观看吧| 看免费av毛片| 午夜a级毛片| 好男人在线观看高清免费视频 | 99国产极品粉嫩在线观看| 男女下面插进去视频免费观看| 国内久久婷婷六月综合欲色啪| 一边摸一边做爽爽视频免费| 中文字幕另类日韩欧美亚洲嫩草| 天堂影院成人在线观看| 国产av在哪里看| 人成视频在线观看免费观看| 久久久久久久久中文| 国产一区在线观看成人免费| av福利片在线| 老司机深夜福利视频在线观看| 精品一品国产午夜福利视频| 青草久久国产| 亚洲午夜精品一区,二区,三区| 国产一卡二卡三卡精品| 狂野欧美激情性xxxx| 又紧又爽又黄一区二区| 欧美日韩亚洲综合一区二区三区_| 色综合欧美亚洲国产小说| 在线观看www视频免费| 一个人免费在线观看的高清视频| 亚洲欧美精品综合久久99| 国产视频一区二区在线看| 一夜夜www| 热99re8久久精品国产| 久久精品亚洲熟妇少妇任你| 亚洲中文av在线| 亚洲欧美精品综合久久99| 在线观看www视频免费| 亚洲性夜色夜夜综合| 亚洲国产中文字幕在线视频| 女人高潮潮喷娇喘18禁视频| 视频在线观看一区二区三区| 中文字幕精品免费在线观看视频| 午夜免费观看网址| 精品午夜福利视频在线观看一区| 亚洲精品中文字幕在线视频| cao死你这个sao货| 欧美日本视频| 99在线人妻在线中文字幕| 91麻豆av在线| 人人妻,人人澡人人爽秒播| 国产av精品麻豆| 日本撒尿小便嘘嘘汇集6| 日韩大尺度精品在线看网址 | 亚洲av日韩精品久久久久久密| 国产欧美日韩一区二区精品| 怎么达到女性高潮| 亚洲电影在线观看av| 两个人免费观看高清视频| 精品熟女少妇八av免费久了| 亚洲欧美激情在线| cao死你这个sao货| 国产亚洲精品一区二区www| 免费av毛片视频| aaaaa片日本免费| 在线观看免费视频日本深夜| 51午夜福利影视在线观看| 91老司机精品| 黄色女人牲交| 精品午夜福利视频在线观看一区| 老汉色∧v一级毛片| 高清黄色对白视频在线免费看| 制服人妻中文乱码| 欧美在线黄色| 在线观看免费午夜福利视频| 动漫黄色视频在线观看| 日本 欧美在线| 午夜日韩欧美国产| 在线天堂中文资源库| 国产激情久久老熟女| 大型黄色视频在线免费观看| 国产av又大| 岛国在线观看网站| 琪琪午夜伦伦电影理论片6080| 欧美另类亚洲清纯唯美| 精品人妻1区二区| 99riav亚洲国产免费| 亚洲专区字幕在线| 亚洲av电影在线进入| 丰满人妻熟妇乱又伦精品不卡| 大型黄色视频在线免费观看| 后天国语完整版免费观看| 亚洲av五月六月丁香网| 黄频高清免费视频| 精品卡一卡二卡四卡免费| 国产私拍福利视频在线观看| 久久久久国产一级毛片高清牌| 国产一区二区三区视频了| 两个人视频免费观看高清| 窝窝影院91人妻| 他把我摸到了高潮在线观看| 亚洲男人的天堂狠狠| 人人妻,人人澡人人爽秒播| 大香蕉久久成人网| 女性被躁到高潮视频| 久久久国产精品麻豆| 国产aⅴ精品一区二区三区波| av天堂在线播放| 如日韩欧美国产精品一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 久久精品国产清高在天天线| 亚洲国产中文字幕在线视频| 国产真人三级小视频在线观看| 好男人在线观看高清免费视频 | 一a级毛片在线观看| 热99re8久久精品国产| 一区二区三区激情视频| 欧美成人一区二区免费高清观看 | 中亚洲国语对白在线视频| 国产熟女xx| 在线天堂中文资源库| 色尼玛亚洲综合影院| 97碰自拍视频| 亚洲aⅴ乱码一区二区在线播放 | 男人舔女人的私密视频| 免费在线观看完整版高清| 国产亚洲av嫩草精品影院| 亚洲在线自拍视频| 日韩免费av在线播放| or卡值多少钱| 老司机午夜福利在线观看视频| 精品国内亚洲2022精品成人| 女人高潮潮喷娇喘18禁视频| 麻豆国产av国片精品| 国产精品1区2区在线观看.| 99国产综合亚洲精品| 欧美一级a爱片免费观看看 | 最好的美女福利视频网| 日韩一卡2卡3卡4卡2021年| 久久婷婷成人综合色麻豆| a级毛片在线看网站| 一本综合久久免费| 亚洲精品粉嫩美女一区| 女人被狂操c到高潮| 中文字幕精品免费在线观看视频| 麻豆成人av在线观看| 99精品欧美一区二区三区四区| 99国产精品一区二区蜜桃av| 丰满人妻熟妇乱又伦精品不卡| 国产视频一区二区在线看| 久久久水蜜桃国产精品网| 99久久综合精品五月天人人| 黑人巨大精品欧美一区二区mp4| 88av欧美| 黄色片一级片一级黄色片| 两性午夜刺激爽爽歪歪视频在线观看 | 国产黄a三级三级三级人| 性欧美人与动物交配| 欧美成人午夜精品| 亚洲无线在线观看| 91麻豆av在线| 欧美成人午夜精品| 国产黄a三级三级三级人| 亚洲免费av在线视频| 日韩国内少妇激情av| 在线观看一区二区三区| 日韩av在线大香蕉| 中文字幕最新亚洲高清| 亚洲一区二区三区色噜噜| 一二三四社区在线视频社区8| 久久久国产成人免费| 香蕉久久夜色| 少妇被粗大的猛进出69影院| 国产精品美女特级片免费视频播放器 | 日韩欧美国产在线观看| 色播在线永久视频| 侵犯人妻中文字幕一二三四区| 国产精品免费一区二区三区在线| 欧美日韩一级在线毛片| 91在线观看av| 国产免费av片在线观看野外av| 亚洲精华国产精华精| 亚洲va日本ⅴa欧美va伊人久久| 国产一区二区在线av高清观看| 中文字幕最新亚洲高清| 亚洲成人国产一区在线观看| 欧美大码av| 精品国产美女av久久久久小说| 久久草成人影院| 久久国产乱子伦精品免费另类| 在线观看免费午夜福利视频| 97碰自拍视频| e午夜精品久久久久久久| 亚洲精品在线观看二区| 操美女的视频在线观看| 亚洲国产精品sss在线观看| 亚洲成a人片在线一区二区| 人人妻,人人澡人人爽秒播| bbb黄色大片| 成人av一区二区三区在线看| 国产麻豆成人av免费视频| 亚洲成人久久性| 久久午夜亚洲精品久久| 老司机靠b影院| 后天国语完整版免费观看| 国产99白浆流出| 两性夫妻黄色片| 在线观看www视频免费| 男人的好看免费观看在线视频 | 9191精品国产免费久久| 99国产综合亚洲精品| 久久精品国产综合久久久| 亚洲久久久国产精品| 热99re8久久精品国产| 可以在线观看的亚洲视频| 中文字幕av电影在线播放| 97人妻精品一区二区三区麻豆 | 欧美不卡视频在线免费观看 | 成人精品一区二区免费| 国产主播在线观看一区二区| 在线观看免费午夜福利视频| 欧美精品亚洲一区二区| 女生性感内裤真人,穿戴方法视频| 女人被狂操c到高潮| 欧美日韩黄片免| 国产高清有码在线观看视频 | 男女做爰动态图高潮gif福利片 | 又紧又爽又黄一区二区| 在线观看免费日韩欧美大片| 国产一级毛片七仙女欲春2 | 看片在线看免费视频| 国产精品日韩av在线免费观看 | 久久人妻福利社区极品人妻图片| 亚洲av成人一区二区三| 久久人妻熟女aⅴ| 天天一区二区日本电影三级 | 亚洲熟女毛片儿| 禁无遮挡网站| 天天躁夜夜躁狠狠躁躁| 亚洲免费av在线视频| 亚洲国产精品999在线| 色综合欧美亚洲国产小说| 国产又色又爽无遮挡免费看| 天堂影院成人在线观看| 长腿黑丝高跟| 在线国产一区二区在线| 亚洲av第一区精品v没综合| 女性被躁到高潮视频| 久久亚洲真实| 精品少妇一区二区三区视频日本电影| 一本久久中文字幕| 欧美亚洲日本最大视频资源| 人人妻人人澡人人看| www.自偷自拍.com| 老司机午夜福利在线观看视频| 亚洲国产欧美日韩在线播放| 视频区欧美日本亚洲| 国产91精品成人一区二区三区| 这个男人来自地球电影免费观看| 亚洲国产中文字幕在线视频| 身体一侧抽搐| 久久天躁狠狠躁夜夜2o2o| 国产伦一二天堂av在线观看| 欧美日本中文国产一区发布| 欧美中文综合在线视频| 国产男靠女视频免费网站| 性欧美人与动物交配| 在线十欧美十亚洲十日本专区| 91麻豆精品激情在线观看国产| 精品国产一区二区久久| 亚洲中文字幕日韩| 69av精品久久久久久| 精品国产国语对白av| 久久精品影院6| 国产aⅴ精品一区二区三区波| 久久国产乱子伦精品免费另类| 国产色视频综合| 成年版毛片免费区| a在线观看视频网站| 亚洲精品在线美女| bbb黄色大片| 亚洲美女黄片视频| 国内久久婷婷六月综合欲色啪| 国产免费av片在线观看野外av| 亚洲第一电影网av| 天堂√8在线中文| 窝窝影院91人妻| 黄色视频,在线免费观看| 757午夜福利合集在线观看| 高清黄色对白视频在线免费看| 亚洲一区二区三区不卡视频| 久久久久久国产a免费观看| 亚洲九九香蕉| 欧美中文综合在线视频| 在线观看免费视频日本深夜| 日韩精品青青久久久久久| 中文字幕人妻熟女乱码| 亚洲精品在线美女| 国产av又大| 国产成人欧美| 一区在线观看完整版| 精品国产亚洲在线| 免费看十八禁软件| 一级a爱视频在线免费观看| 两个人视频免费观看高清| 国产熟女xx| 淫秽高清视频在线观看| 脱女人内裤的视频| 日本精品一区二区三区蜜桃| 亚洲天堂国产精品一区在线| 男女下面进入的视频免费午夜 | 久久久久久国产a免费观看| 国产在线精品亚洲第一网站| 亚洲视频免费观看视频| 午夜成年电影在线免费观看| 巨乳人妻的诱惑在线观看| 夜夜夜夜夜久久久久| 久久欧美精品欧美久久欧美| 免费在线观看黄色视频的| 亚洲全国av大片| 国产亚洲av嫩草精品影院| 日韩中文字幕欧美一区二区| 母亲3免费完整高清在线观看| 高潮久久久久久久久久久不卡| 亚洲精品久久成人aⅴ小说| 一个人观看的视频www高清免费观看 | 两个人视频免费观看高清| 一本综合久久免费| 黑丝袜美女国产一区| 一级毛片高清免费大全| 黄色视频不卡| 久久青草综合色| 天堂动漫精品| 夜夜躁狠狠躁天天躁| 免费在线观看完整版高清| 国产精品一区二区在线不卡| 丰满的人妻完整版| 国内久久婷婷六月综合欲色啪| 国产精品久久久av美女十八| 在线av久久热| 女人高潮潮喷娇喘18禁视频| 日本vs欧美在线观看视频| 国产精品久久久久久亚洲av鲁大| 69精品国产乱码久久久| 国产三级在线视频| 国产精品一区二区在线不卡| 91成年电影在线观看| 欧美激情久久久久久爽电影 | 神马国产精品三级电影在线观看 | 非洲黑人性xxxx精品又粗又长| 日韩成人在线观看一区二区三区| √禁漫天堂资源中文www| 国产欧美日韩精品亚洲av| 一本综合久久免费| 男人的好看免费观看在线视频 | 妹子高潮喷水视频| 国产av又大| 午夜视频精品福利| 欧美色视频一区免费| 国语自产精品视频在线第100页| 90打野战视频偷拍视频| 可以在线观看毛片的网站| 一二三四社区在线视频社区8| 亚洲成国产人片在线观看| 在线视频色国产色| 少妇的丰满在线观看| 免费看a级黄色片| 国产成人一区二区三区免费视频网站| 欧美亚洲日本最大视频资源| 亚洲色图 男人天堂 中文字幕| 国产一区二区三区视频了| 首页视频小说图片口味搜索| 久久中文看片网| 国产亚洲精品综合一区在线观看 | 久久精品人人爽人人爽视色| 国产精品98久久久久久宅男小说| 99香蕉大伊视频| 啦啦啦免费观看视频1| 日韩大尺度精品在线看网址 | 在线观看免费视频日本深夜| 免费不卡黄色视频| 免费在线观看完整版高清| 好看av亚洲va欧美ⅴa在| 嫁个100分男人电影在线观看| 日韩精品青青久久久久久| 丰满人妻熟妇乱又伦精品不卡| 啦啦啦 在线观看视频| 国产精品国产高清国产av| 99香蕉大伊视频| 男女下面进入的视频免费午夜 | 不卡av一区二区三区| 欧美乱妇无乱码| 亚洲第一av免费看| 久久青草综合色| 亚洲国产精品成人综合色| 久久久久久人人人人人| 国产精品久久电影中文字幕| 久久久久久大精品| 侵犯人妻中文字幕一二三四区| 日本一区二区免费在线视频| 亚洲片人在线观看| 在线国产一区二区在线| 一本久久中文字幕| 国产在线观看jvid| 成人欧美大片| 久久精品91无色码中文字幕| av中文乱码字幕在线| 久久人人爽av亚洲精品天堂| 51午夜福利影视在线观看| 午夜视频精品福利| 色尼玛亚洲综合影院| 嫩草影院精品99| 精品国产亚洲在线| 成人三级做爰电影| 欧美日韩中文字幕国产精品一区二区三区 | 大型黄色视频在线免费观看| 亚洲第一电影网av| 免费看美女性在线毛片视频| 麻豆久久精品国产亚洲av| 美女扒开内裤让男人捅视频| 久久久久久大精品| 夜夜看夜夜爽夜夜摸| 久久久久国产一级毛片高清牌| 亚洲男人天堂网一区| 欧美绝顶高潮抽搐喷水| 欧美日韩一级在线毛片| 他把我摸到了高潮在线观看| 成人永久免费在线观看视频| 69精品国产乱码久久久| 日韩欧美一区视频在线观看| 国产成人精品久久二区二区免费| 色精品久久人妻99蜜桃| 国产野战对白在线观看| 啦啦啦韩国在线观看视频| 色精品久久人妻99蜜桃| 亚洲美女黄片视频| 99精品在免费线老司机午夜| 午夜福利在线观看吧| 久久久国产欧美日韩av| 国产99久久九九免费精品| 久久天躁狠狠躁夜夜2o2o| 亚洲美女黄片视频| 可以在线观看毛片的网站| 在线观看午夜福利视频| 久久人人精品亚洲av| 亚洲 欧美一区二区三区| 色精品久久人妻99蜜桃| 最好的美女福利视频网| 免费观看精品视频网站| 色精品久久人妻99蜜桃| 波多野结衣av一区二区av| 日韩大尺度精品在线看网址 | 亚洲熟女毛片儿| 欧美成人免费av一区二区三区| 老汉色av国产亚洲站长工具| 亚洲精品一卡2卡三卡4卡5卡| 美女 人体艺术 gogo| 亚洲成av片中文字幕在线观看| 又黄又粗又硬又大视频| 69精品国产乱码久久久| 亚洲激情在线av| 在线十欧美十亚洲十日本专区| 波多野结衣一区麻豆|