• <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
    免费人成视频x8x8入口观看| 免费电影在线观看免费观看| 国产男靠女视频免费网站| 女同久久另类99精品国产91| 窝窝影院91人妻| 精品久久久久久成人av| 中文在线观看免费www的网站| 老司机午夜福利在线观看视频| 亚洲欧美精品综合久久99| 中文字幕人妻丝袜一区二区| 99在线人妻在线中文字幕| 亚洲18禁久久av| 国产精品自产拍在线观看55亚洲| 欧美日韩中文字幕国产精品一区二区三区| 好男人电影高清在线观看| 久久精品亚洲精品国产色婷小说| 久久久久久久久中文| 看免费av毛片| 日本一本二区三区精品| 99久久成人亚洲精品观看| 三级毛片av免费| 麻豆国产av国片精品| 午夜福利视频1000在线观看| 午夜免费激情av| 亚洲五月天丁香| 欧美成人性av电影在线观看| 日韩av在线大香蕉| 欧美色视频一区免费| 欧美国产日韩亚洲一区| 亚洲无线观看免费| 国产午夜精品久久久久久一区二区三区 | 亚洲欧美日韩高清在线视频| 日本黄色视频三级网站网址| 亚洲精品色激情综合| 夜夜爽天天搞| 丰满人妻熟妇乱又伦精品不卡| 老司机深夜福利视频在线观看| 免费av不卡在线播放| 日韩成人在线观看一区二区三区| 亚洲av一区综合| 少妇人妻精品综合一区二区 | 国产成人a区在线观看| 欧美日韩精品网址| 国内揄拍国产精品人妻在线| 亚洲精品影视一区二区三区av| 又爽又黄无遮挡网站| 中文字幕高清在线视频| 男人舔奶头视频| 亚洲av成人av| 亚洲黑人精品在线| 国产69精品久久久久777片| 国产伦精品一区二区三区四那| 一区二区三区激情视频| 五月玫瑰六月丁香| 日韩欧美精品免费久久 | 老司机深夜福利视频在线观看| 此物有八面人人有两片| 天天添夜夜摸| 成人性生交大片免费视频hd| 脱女人内裤的视频| 91久久精品电影网| 女生性感内裤真人,穿戴方法视频| 国产亚洲精品久久久com| 亚洲人成网站在线播| 熟妇人妻久久中文字幕3abv| 3wmmmm亚洲av在线观看| 一个人免费在线观看电影| 日本撒尿小便嘘嘘汇集6| 丁香六月欧美| 精品电影一区二区在线| 90打野战视频偷拍视频| 成年免费大片在线观看| 亚洲欧美日韩东京热| 亚洲精品一区av在线观看| 国产精品日韩av在线免费观看| 最近在线观看免费完整版| 国产 一区 欧美 日韩| 婷婷精品国产亚洲av| 久久精品国产自在天天线| 精品免费久久久久久久清纯| 性欧美人与动物交配| 精品久久久久久久人妻蜜臀av| 欧美日韩精品网址| 国产黄片美女视频| 久久久久亚洲av毛片大全| 亚洲成人免费电影在线观看| 精品电影一区二区在线| 亚洲精品日韩av片在线观看 | 久久精品影院6| 国产色婷婷99| 亚洲欧美日韩高清专用| 亚洲欧美日韩卡通动漫| 亚洲精品粉嫩美女一区| 国产av一区在线观看免费| 国产一区二区三区视频了| 国内精品久久久久久久电影| 久久6这里有精品| 波多野结衣巨乳人妻| 亚洲av熟女| 少妇人妻一区二区三区视频| 香蕉久久夜色| 午夜激情福利司机影院| 18禁在线播放成人免费| 久久这里只有精品中国| 全区人妻精品视频| 非洲黑人性xxxx精品又粗又长| 在线免费观看的www视频| 老司机深夜福利视频在线观看| 一进一出抽搐动态| 别揉我奶头~嗯~啊~动态视频| 国产69精品久久久久777片| 免费人成在线观看视频色| 国产久久久一区二区三区| 成人精品一区二区免费| 国产老妇女一区| 亚洲国产精品sss在线观看| 一区福利在线观看| 天堂动漫精品| 欧美性猛交黑人性爽| 老鸭窝网址在线观看| 午夜精品在线福利| 国产中年淑女户外野战色| 亚洲欧美精品综合久久99| 中文在线观看免费www的网站| 国产伦一二天堂av在线观看| 国产爱豆传媒在线观看| 精品国产美女av久久久久小说| 欧美日韩乱码在线| 欧美zozozo另类| 国产中年淑女户外野战色| 久久久精品大字幕| 九九在线视频观看精品| 久久精品91蜜桃| 精品国产亚洲在线| 12—13女人毛片做爰片一| 免费看日本二区| 国产精品日韩av在线免费观看| 男人舔女人下体高潮全视频| 12—13女人毛片做爰片一| netflix在线观看网站| av天堂中文字幕网| 在线观看66精品国产| 色综合婷婷激情| 久久久国产成人精品二区| 窝窝影院91人妻| 国产精品久久久久久久电影 | 日韩高清综合在线| 老司机福利观看| 青草久久国产| 亚洲精品456在线播放app | 国产激情偷乱视频一区二区| 成人亚洲精品av一区二区| 久久久久性生活片| 看黄色毛片网站| 大型黄色视频在线免费观看| 深夜精品福利| 高清毛片免费观看视频网站| 久99久视频精品免费| www.熟女人妻精品国产| 少妇熟女aⅴ在线视频| 悠悠久久av| 久久这里只有精品中国| 国产黄a三级三级三级人| 日韩欧美免费精品| 国产精品美女特级片免费视频播放器| 日日夜夜操网爽| 色哟哟哟哟哟哟| 午夜免费成人在线视频| 亚洲18禁久久av| 欧美色视频一区免费| 亚洲在线自拍视频| 久久精品91蜜桃| 日韩成人在线观看一区二区三区| 女生性感内裤真人,穿戴方法视频| 午夜老司机福利剧场| 特大巨黑吊av在线直播| xxx96com| 国产又黄又爽又无遮挡在线| 国产黄色小视频在线观看| 别揉我奶头~嗯~啊~动态视频| 欧美中文综合在线视频| 波多野结衣巨乳人妻| 亚洲成人久久性| 国产精品久久久久久亚洲av鲁大| 一二三四社区在线视频社区8| 深爱激情五月婷婷| 久久久精品欧美日韩精品| а√天堂www在线а√下载| 搡女人真爽免费视频火全软件 | 亚洲精品亚洲一区二区| 中文字幕高清在线视频| 中出人妻视频一区二区| svipshipincom国产片| 成人18禁在线播放| 亚洲成av人片在线播放无| 悠悠久久av| 成年女人永久免费观看视频| 精品欧美国产一区二区三| 黑人欧美特级aaaaaa片| 白带黄色成豆腐渣| 久久精品91蜜桃| 波多野结衣高清无吗| 国产精品av视频在线免费观看| 精品国产美女av久久久久小说| 九九久久精品国产亚洲av麻豆| 久久久久久国产a免费观看| 伊人久久大香线蕉亚洲五| 日韩精品中文字幕看吧| 国产精品亚洲美女久久久| 成人欧美大片| 99久久精品一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 亚洲成人久久爱视频| 最近最新免费中文字幕在线| 精品免费久久久久久久清纯| 嫩草影院精品99| 99国产极品粉嫩在线观看| 在线观看舔阴道视频| 亚洲第一电影网av| 在线观看免费午夜福利视频| x7x7x7水蜜桃| 精品久久久久久久久久免费视频| 尤物成人国产欧美一区二区三区| 国产精品一区二区三区四区免费观看 | 色综合婷婷激情| 免费av不卡在线播放| 亚洲 国产 在线| 国产高清激情床上av| 亚洲,欧美精品.| 不卡一级毛片| 成人特级黄色片久久久久久久| 婷婷亚洲欧美| 国产伦一二天堂av在线观看| 国产亚洲精品av在线| 好看av亚洲va欧美ⅴa在| 又黄又粗又硬又大视频| 91九色精品人成在线观看| 亚洲av成人不卡在线观看播放网| 久久久久久久午夜电影| 精品无人区乱码1区二区| 久久久久久久久中文| 婷婷精品国产亚洲av在线| 嫩草影院精品99| 淫秽高清视频在线观看| 亚洲熟妇熟女久久| 欧美av亚洲av综合av国产av| 午夜福利免费观看在线| 久久国产乱子伦精品免费另类| 亚洲国产精品合色在线| 国产伦精品一区二区三区四那| 动漫黄色视频在线观看| 亚洲真实伦在线观看| 在线观看av片永久免费下载| 白带黄色成豆腐渣| 男女做爰动态图高潮gif福利片| 亚洲aⅴ乱码一区二区在线播放| 精品人妻一区二区三区麻豆 | 韩国av一区二区三区四区| 亚洲精品亚洲一区二区| 色综合站精品国产| 精品免费久久久久久久清纯| 久久久久性生活片| 淫妇啪啪啪对白视频| 国产午夜精品论理片| 亚洲国产日韩欧美精品在线观看 | 久久久久久久久大av| 18禁国产床啪视频网站| 九九久久精品国产亚洲av麻豆| 国产精品一及| 国产黄a三级三级三级人| 99在线人妻在线中文字幕| 少妇的逼好多水| 精品欧美国产一区二区三| 看片在线看免费视频| 好男人电影高清在线观看| 深爱激情五月婷婷| 亚洲avbb在线观看| 全区人妻精品视频| 亚洲内射少妇av| 国内精品久久久久久久电影| 日本免费a在线| 91av网一区二区| 真实男女啪啪啪动态图| 亚洲精品影视一区二区三区av| 在线观看一区二区三区| 国产探花极品一区二区| 一级作爱视频免费观看| 黄片小视频在线播放| 久久香蕉国产精品| 亚洲av熟女| 18禁黄网站禁片午夜丰满| 日本在线视频免费播放| 99久久99久久久精品蜜桃| 99久久九九国产精品国产免费| 久久久久久久亚洲中文字幕 | 18禁国产床啪视频网站| 精品免费久久久久久久清纯| 69人妻影院| 日韩欧美精品免费久久 | 免费观看人在逋| 日本 欧美在线| 成人性生交大片免费视频hd| 日日夜夜操网爽| 老司机深夜福利视频在线观看| 中文亚洲av片在线观看爽| 1000部很黄的大片| 人人妻,人人澡人人爽秒播| 久久久久久人人人人人| 97碰自拍视频| 亚洲精品日韩av片在线观看 | 国产精品影院久久| 波多野结衣高清无吗| 国产精品女同一区二区软件 | 久久欧美精品欧美久久欧美| 99国产综合亚洲精品| 精品99又大又爽又粗少妇毛片 | 最近在线观看免费完整版| 丰满人妻熟妇乱又伦精品不卡| 国产精品一区二区三区四区免费观看 | 村上凉子中文字幕在线| 最好的美女福利视频网| www.色视频.com| 欧美av亚洲av综合av国产av| 一个人免费在线观看电影| 日本一二三区视频观看| 麻豆成人午夜福利视频| 一级黄片播放器| 人妻久久中文字幕网| 久久久久国内视频| 欧美黑人巨大hd| 欧美xxxx黑人xx丫x性爽| 久久久久久久久久黄片| 欧美一区二区精品小视频在线| 伊人久久精品亚洲午夜| 悠悠久久av| 男女床上黄色一级片免费看| 欧美xxxx黑人xx丫x性爽| 亚洲avbb在线观看| 搞女人的毛片| 亚洲真实伦在线观看| 国产精品av视频在线免费观看| 亚洲专区中文字幕在线| 久久精品国产自在天天线| 搡女人真爽免费视频火全软件 | 久久精品影院6| 岛国视频午夜一区免费看| 天天添夜夜摸| bbb黄色大片| 国产亚洲欧美98| 久久亚洲精品不卡| 高潮久久久久久久久久久不卡| ponron亚洲| 久久精品综合一区二区三区| 欧美日韩福利视频一区二区| 国产伦精品一区二区三区四那| 搡老岳熟女国产| 精品乱码久久久久久99久播| 国产精品久久久久久精品电影| 免费电影在线观看免费观看| 香蕉丝袜av| 午夜精品久久久久久毛片777| 啦啦啦观看免费观看视频高清| 12—13女人毛片做爰片一| 91麻豆精品激情在线观看国产| 一二三四社区在线视频社区8| 老熟妇乱子伦视频在线观看| 成人国产一区最新在线观看| 国产av不卡久久| 成年免费大片在线观看| 国产又黄又爽又无遮挡在线| 看片在线看免费视频| 在线a可以看的网站| 亚洲美女黄片视频| 国产精品三级大全| 国产精品av视频在线免费观看| 搡女人真爽免费视频火全软件 | 国产精品影院久久| 亚洲国产精品999在线| 母亲3免费完整高清在线观看| 美女被艹到高潮喷水动态| 热99在线观看视频| 精品欧美国产一区二区三| 激情在线观看视频在线高清| 露出奶头的视频| 婷婷丁香在线五月| 亚洲电影在线观看av| 观看美女的网站| 美女被艹到高潮喷水动态| 亚洲 国产 在线| 午夜亚洲福利在线播放| 欧美黄色淫秽网站| 国产一级毛片七仙女欲春2| 久久久国产精品麻豆| 1024手机看黄色片| 亚洲国产精品久久男人天堂| 日韩人妻高清精品专区| 午夜福利在线在线| 色综合亚洲欧美另类图片| 欧美日韩乱码在线| 在线看三级毛片| 欧美高清成人免费视频www| 中文亚洲av片在线观看爽| 亚洲美女视频黄频| 免费看十八禁软件| 亚洲va日本ⅴa欧美va伊人久久| 久久久久久大精品| 欧美大码av| 久久精品夜夜夜夜夜久久蜜豆| 91在线精品国自产拍蜜月 | 日韩国内少妇激情av| 欧美+亚洲+日韩+国产| 久久香蕉国产精品| 日本 欧美在线| 在线观看免费午夜福利视频| 久久精品国产99精品国产亚洲性色| 国产高清三级在线| 三级男女做爰猛烈吃奶摸视频| 高清日韩中文字幕在线| 91在线观看av| 手机成人av网站| 免费在线观看日本一区| 久久国产乱子伦精品免费另类| 1024手机看黄色片| 一级黄色大片毛片| 黄色视频,在线免费观看| 男女那种视频在线观看| 精品国产三级普通话版| 日本黄色片子视频| 亚洲专区中文字幕在线| 午夜福利视频1000在线观看| 手机成人av网站| 香蕉丝袜av| 色播亚洲综合网| 国语自产精品视频在线第100页| 美女高潮喷水抽搐中文字幕| 日韩欧美在线二视频| 国内少妇人妻偷人精品xxx网站| 成人午夜高清在线视频| 日韩大尺度精品在线看网址| 国产私拍福利视频在线观看| 日韩欧美在线二视频| 在线视频色国产色| 久久午夜亚洲精品久久| 色老头精品视频在线观看| 日本三级黄在线观看| 日韩欧美一区二区三区在线观看| 丰满的人妻完整版| 亚洲无线在线观看| 亚洲黑人精品在线| 国产成人av教育| 亚洲精品456在线播放app | 韩国av一区二区三区四区| 欧美一区二区精品小视频在线| 男女床上黄色一级片免费看| 给我免费播放毛片高清在线观看| 午夜久久久久精精品| 亚洲精品在线美女| 麻豆国产97在线/欧美| aaaaa片日本免费| 又黄又爽又免费观看的视频| 久久久成人免费电影| 久久性视频一级片| av专区在线播放| 亚洲美女黄片视频| 欧美日韩一级在线毛片| 精品国产亚洲在线| 国产精品久久视频播放| 99久久九九国产精品国产免费| 我要搜黄色片| 看黄色毛片网站| 国产成人福利小说| 国产av一区在线观看免费| 3wmmmm亚洲av在线观看| 国产69精品久久久久777片| 欧美又色又爽又黄视频| 在线观看免费午夜福利视频| 老汉色∧v一级毛片| 日本一二三区视频观看| av福利片在线观看| 欧美成人a在线观看| 手机成人av网站| 久久久久久久精品吃奶| 久久精品国产亚洲av涩爱 | 亚洲真实伦在线观看| 成人鲁丝片一二三区免费| 黄色女人牲交| 看片在线看免费视频| 亚洲欧美一区二区三区黑人| 夜夜看夜夜爽夜夜摸| 黄色片一级片一级黄色片| 国产成人系列免费观看| 男女那种视频在线观看| av在线天堂中文字幕| 黄色日韩在线| 欧美日韩一级在线毛片| 亚洲 国产 在线| 日本撒尿小便嘘嘘汇集6| 99国产极品粉嫩在线观看| 偷拍熟女少妇极品色| 一个人免费在线观看电影| 国产伦精品一区二区三区四那| 99久久精品热视频| 亚洲欧美精品综合久久99| 精品国产三级普通话版| 国产精品永久免费网站| 国产伦精品一区二区三区视频9 | 两性午夜刺激爽爽歪歪视频在线观看| 精品久久久久久久毛片微露脸| 亚洲成av人片在线播放无| 变态另类丝袜制服| 欧美一级a爱片免费观看看| 亚洲av免费在线观看| 99热这里只有是精品50| 欧美黑人欧美精品刺激| 国产探花在线观看一区二区| 国产三级黄色录像| 精品一区二区三区av网在线观看| 十八禁网站免费在线| 精品人妻一区二区三区麻豆 | 高清日韩中文字幕在线| 女人被狂操c到高潮| 日韩精品青青久久久久久| 真人一进一出gif抽搐免费| 久久香蕉精品热| 最后的刺客免费高清国语| 久久香蕉精品热| 久久九九热精品免费| 三级男女做爰猛烈吃奶摸视频| 最近最新免费中文字幕在线| 成人高潮视频无遮挡免费网站| av视频在线观看入口| 国产精品一区二区三区四区免费观看 | 精品乱码久久久久久99久播| 亚洲第一电影网av| 免费看a级黄色片| 国产亚洲精品久久久com| 国产免费男女视频| 九色成人免费人妻av| 欧美极品一区二区三区四区| 一进一出抽搐gif免费好疼| 午夜福利在线观看吧| 熟妇人妻久久中文字幕3abv| 狠狠狠狠99中文字幕| 白带黄色成豆腐渣| 性色avwww在线观看| 国产99白浆流出| 别揉我奶头~嗯~啊~动态视频| 好男人电影高清在线观看| xxxwww97欧美| 日韩 欧美 亚洲 中文字幕| 国产久久久一区二区三区| 国产私拍福利视频在线观看| 高潮久久久久久久久久久不卡| 19禁男女啪啪无遮挡网站| 看黄色毛片网站| 香蕉丝袜av| 欧美激情久久久久久爽电影| 午夜福利在线观看吧| 精品国产亚洲在线| 亚洲人成电影免费在线| 精品人妻1区二区| 成人鲁丝片一二三区免费| eeuss影院久久| 在线播放国产精品三级| 亚洲精品色激情综合| 免费一级毛片在线播放高清视频| 天天添夜夜摸| 亚洲国产精品合色在线| 在线观看午夜福利视频| 我的老师免费观看完整版| 美女 人体艺术 gogo| 日韩欧美精品v在线| 精品久久久久久成人av| 国产私拍福利视频在线观看| 99精品在免费线老司机午夜| 亚洲精品国产精品久久久不卡| 国产亚洲精品av在线| 午夜老司机福利剧场| 免费看a级黄色片| 国产在线精品亚洲第一网站| 狠狠狠狠99中文字幕| 久久久久久九九精品二区国产| 亚洲人成网站在线播| 久久久国产精品麻豆| 国产视频一区二区在线看| 性欧美人与动物交配| 老汉色av国产亚洲站长工具| 激情在线观看视频在线高清| 非洲黑人性xxxx精品又粗又长| 免费电影在线观看免费观看| 老司机深夜福利视频在线观看| 国产淫片久久久久久久久 | 亚洲精品美女久久久久99蜜臀| 一个人看视频在线观看www免费 | 在线视频色国产色| 啦啦啦韩国在线观看视频| 欧美黑人巨大hd| 2021天堂中文幕一二区在线观| 夜夜躁狠狠躁天天躁| 99久久无色码亚洲精品果冻| 精品人妻1区二区| 夜夜躁狠狠躁天天躁| 啪啪无遮挡十八禁网站| 国产精品影院久久| 亚洲在线观看片| 男女视频在线观看网站免费| 久久欧美精品欧美久久欧美| 男女那种视频在线观看| 久久久久久久精品吃奶| 人人妻,人人澡人人爽秒播| 亚洲在线观看片| 日本 欧美在线| 日韩欧美精品v在线| 国产欧美日韩精品一区二区| 99riav亚洲国产免费|