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

    利用常熟地震臺同場地形變儀記錄檢測對比地球自由振蕩信號

    2021-12-07 01:16:18錢文杰劉冬冬陸德明丁建國
    地震地磁觀測與研究 2021年4期
    關(guān)鍵詞:環(huán)型應(yīng)變儀球型

    錢文杰 狄 樑 劉冬冬 陸德明 丁建國

    (中國江蘇 215500 江蘇省地震局常熟地震臺)

    0 引言

    地球自由振蕩是一種重要的地球物理和動力學(xué)現(xiàn)象,是指地球局部受到某種因素激發(fā)時,地球整體產(chǎn)生的連續(xù)振動,如在受到大地震、火山爆發(fā)或地下核爆炸激發(fā)后,會發(fā)生整體振動,并能持續(xù)一段時間。地球自由振蕩可以變換分解為一些特定的頻率,稱為地球的本征頻率,對應(yīng)2 種本征振蕩:①球型地球自由振蕩,其質(zhì)點位移既有徑向運動,也有水平運動;②環(huán)型地球自由振蕩,其質(zhì)點位移只作水平運動。地球自由振蕩的本征頻率可以反映地球內(nèi)部介質(zhì)的固有特征,準(zhǔn)確獲得自由振蕩的頻率特征,可以有效反映地球內(nèi)部的結(jié)構(gòu)信息。19 世紀(jì)末期,Lamb 從理論上推導(dǎo)出地球自由振蕩的存在;1961 年,Benionff 等使用應(yīng)變地震儀和擺式儀,準(zhǔn)確觀測到智利大地震激發(fā)的長周期地球自由振蕩,確認(rèn)了地震自由振蕩的真實存在(邱澤華等,2007;唐磊等,2007)。21 世紀(jì)以來,國內(nèi)外進行地球自由振蕩觀測主要采用寬頻帶地震儀、重力儀、應(yīng)變儀、傾斜儀等靈敏的、可記錄長周期振動的高精度儀器。隨著各類觀測資料的積累,地球自由振蕩的理論研究及觀測資料的分析已成為探查地球內(nèi)部結(jié)構(gòu)的有效手段。

    20 世紀(jì)70 年代,國際上開始利用傾斜儀觀測記錄地球自由振蕩信號。目前,對于大地震激發(fā)的地球自由振蕩,國內(nèi)多針對單臺或多臺單型號的應(yīng)變儀、形變儀記錄進行球型自由振蕩振型的檢測分析,取得不少研究成果,如:邱澤華等(2007)利用泰安地震臺鉆孔差應(yīng)變儀觀測資料、唐磊等(2007)利用中國鉆孔應(yīng)變臺網(wǎng)體應(yīng)變觀測資料,楊躍文等(2010)利用云南水管儀觀測資料,檢測2004 年蘇門答臘地震激發(fā)的地球自由振蕩;任佳等(2009)利用垂直擺儀和水管儀觀測資料,檢測2008 年汶川地震激發(fā)的地球自由振蕩。利用同一場地不同形變測項(傾斜、應(yīng)變)觀測資料進行自由振蕩檢測,特別是環(huán)型自由振蕩振型的檢測,相關(guān)文獻較少。

    常熟地震臺(下文簡稱常熟臺)形變觀測山洞進深220 m,環(huán)境良好,自2007 年起,安裝DSQ 水管傾斜儀(下文簡稱水管儀)、VS 垂直擺傾斜儀(下文簡稱垂直擺)、SS-Y 洞體應(yīng)變儀(下文簡稱洞體應(yīng)變儀)、TJ-2 體應(yīng)變儀(下文簡稱體應(yīng)變儀)等數(shù)字化觀測儀器,記錄的固體潮汐和同震響應(yīng)清晰,各類測項觀測精度均達全國形變Ⅰ類臺標(biāo)準(zhǔn)。

    利用常熟臺水管儀、垂直擺、洞體應(yīng)變儀、體應(yīng)變儀觀測資料,對2011 年3 月11 日日本9.0 級大地震激發(fā)的地球基頻球型振蕩、諧頻球型振蕩、基頻環(huán)型振蕩進行檢測,并與現(xiàn)有PREM 模型理論值進行對比分析,研究同一場地各類形變觀測儀對地球自由振蕩信號的檢測能力及共異特性,檢定該臺站形變觀測儀器的運行狀態(tài)的同時,進一步證明形變儀器對自由振蕩的檢測效能,檢測結(jié)果可為約束和改善地球自由振蕩模型理論值提供參考。

    1 觀測數(shù)據(jù)選取及預(yù)處理

    2011 年3 月11 日13 時46 分在日本本州東海岸附近海域發(fā)生9.0 級大地震,在造成巨大災(zāi)害的同時,為研究地球自由振蕩提供了機會,如姚家俊等(2012)、欒威等(2015)、孟方杰等(2018)均對此次地震激發(fā)的地球自由振蕩進行了相關(guān)分析。常熟臺水管儀、垂直擺、洞體應(yīng)變儀、體應(yīng)變儀清晰記錄到此次大地震,其中垂直擺NS 向受幅值限制,震時出現(xiàn)暫短靠擺現(xiàn)象,各儀器記錄曲線(分鐘值)見圖1,圖中所用數(shù)據(jù)均未經(jīng)干擾消除處理。

    圖1 常熟形變觀測儀對日本大地震的分鐘值記錄曲線(a)水管儀EW、NS 向;(b)垂直擺EW、NS 向;(c)洞體應(yīng)變儀EW、NS 向;(d)體應(yīng)變儀Fig.1 Minute value curves of Japan earthquake recorded by Changshu deformation observatory

    據(jù)Rosat 等(2003)的研究,自由振蕩檢測的最優(yōu)起始時間為震后5 h(Rosat,2005)。因此,選取常熟臺水管儀、垂直擺、洞體應(yīng)變儀、體應(yīng)變儀2011 年3 月11 日17 時46 分—16 日17 時46 分共5 天7 200 個分鐘值采樣數(shù)據(jù),提取日本9 級大地震激發(fā)的地球自由振蕩振型。為去除固體潮影響,選用db4 小波變換方法對觀測數(shù)據(jù)進行6 層濾波預(yù)處理(丁浩等,2013),剔除固體潮汐成分,得到數(shù)據(jù)殘差值。

    大地震激發(fā)的地球自由振蕩信號可持續(xù)數(shù)日至數(shù)月,但振蕩信號是隨著時間衰減的,選取合適的數(shù)據(jù)長度是能否有效檢測更多振蕩信號的前提,同時時間長度越長,觀測資料攜帶的背景噪聲越多,但某些振型之間的頻率間隔太小,若將這些振型清晰分辨出來,要求數(shù)據(jù)長度(丁浩等,2013)為

    式中,N為采樣點數(shù),Δt為采樣間隔,Δω為最小頻率間隔。

    綜合考慮頻譜分辨率與背景噪聲水平因素,通過實測比對,選取一定數(shù)量的殘差數(shù)據(jù)進行功率譜計算,其中:①洞體應(yīng)變儀:4 000 個殘差數(shù)據(jù);②水管儀、垂直擺、體應(yīng)變儀:各5 000 個殘差數(shù)據(jù)。

    功率譜用于描述隨機信號能量特征隨頻率的變化關(guān)系,通常用來描述隨機信號的頻域特征,也是檢測地球自由振蕩的主要方法。采用改進的平均周期圖法,求取隨機信號的功率譜密度估計,并運用信號重疊分段、加窗函數(shù)與FFT 算法等,提高運算效率與效果。設(shè)信號X的自相關(guān)函數(shù)為Rn,定義其Fourier 變換Sk為X的功率譜密度估計(雷湘鄂,2007),即

    其中,自相關(guān)函數(shù)Rn為

    式(1)、(2)中,N為殘差值序列數(shù)值長度,x為殘差值序列數(shù)值,Sk離散值為功率譜密度值。為了消除數(shù)據(jù)不能無限長而必須加窗對功率譜密度估計造成的影響,采用Hanning 窗抑制功率譜旁瓣,突出主瓣。

    對于檢測頻率值是否有效可信,文中參考雷湘鄂等(2007)的甄別方法,利用檢測振型的信噪比(每個檢測振型的譜峰值與該振型附近觀測背景噪聲譜的比值,即該振型的信噪比)予以判別,即若被檢測振型的信噪比大于3,則檢測結(jié)果有效,而且檢測頻率值與理論頻率值不能相差太大(胡小剛,2006)。

    2 功率譜密度計算結(jié)果

    采用上述方法,計算常熟臺4 組形變儀觀測數(shù)據(jù)在2011 年日本大地震后的功率譜密度估計譜值(PSD),結(jié)果見圖2、圖3、圖4,觀測頻段分別為0.4×10-3—1.8×10-3Hz、1.8×10-3—3.5×10-3Hz 和3.5×10-3—5.6×10-3Hz,圖中虛線為地球初步參考模型(PREM)的理論頻率值,實線表示實測結(jié)果。

    圖2 (0.4—1.8)×10-3 Hz 頻段形變觀測數(shù)據(jù)功率譜密度估計曲線(a)水管儀EW 向;(b)水管儀NS 向;(c)垂直擺EW 向;(d)垂直擺NS 向;(e)洞體應(yīng)變儀EW 向;(f)洞體應(yīng)變儀NS 向;(g)體應(yīng)變儀Fig.2 Power spectral density estimation curves of deformation observation data in(0.4—1.8)×10-3 Hz frequency band

    圖3 (1.8—3.5)×10-3 Hz 頻段形變觀測數(shù)據(jù)功率譜密度估計曲線(a)水管儀EW 向;(b)水管儀NS 向;(c)垂直擺EW 向;(d)垂直擺NS 向;(e)洞體應(yīng)變儀EW 向;(f)洞體應(yīng)變儀NS 向;(g)體應(yīng)變儀Fig.3 Power spectral density estimation curves of deformation observation data in(1.8—3.5)×10-3 Hz frequency band

    圖4 (3.5—5.6)×10-3 Hz 頻段形變觀測數(shù)據(jù)功率譜密度估計曲線(a)水管儀EW 向;(b)水管儀NS 向;(c)垂直擺EW 向;(d)垂直擺NS 向;(e)洞體應(yīng)變儀EW 向;(f)洞體應(yīng)變儀NS 向;(g)體應(yīng)變儀Fig.4 Power spectral density estimation curves of deformation observation data in(3.5—5.6)×10-3 Hz frequency band

    圖2 顯示了(0.4—1.8)×10-3Hz 頻段內(nèi)檢測的常熟臺形變儀觀測數(shù)據(jù)功率譜密度(相對值)與PREM 模型給出的自由振蕩頻率的對應(yīng)情況。在此頻段,常熟臺形變觀測資料檢測到基頻球型振型0S3—0S10,基頻環(huán)型振型0T5—0T11以及諧頻球型1S2—1S4、2S4—2S6、3S1。其中水管儀、體應(yīng)變儀檢測到低階振型0S3、0S4,且振型譜峰較清晰,信噪比>3;垂直擺檢測到0S4振型,未檢測到0S3振型;洞體應(yīng)變儀檢測到0S3、0S4,振型觀測效果不明顯,觀測值譜峰相對于周邊噪聲水平未明顯分開。

    圖3 顯示了(1.8—3.5)×10-3Hz 頻段內(nèi)檢測的常熟臺形變儀器觀測數(shù)據(jù)功率譜密度(相對值)與PREM 模型給出的自由振蕩頻率對應(yīng)情況。在此頻段,常熟臺形變觀測資料檢測到基頻球型振型0S11—0S26,基頻環(huán)型振型0T13—0T18、0T24、0T25,諧頻球型2S8、2S10、3S9。該頻段內(nèi)無較大干擾,各個振型附近噪聲水平較低,觀測振型譜峰清晰,與理論值偏離小,信噪比較高。

    圖4 顯示了(3.5—5.6)×10-3Hz 頻段內(nèi)檢測到的常熟臺形變儀器觀測數(shù)據(jù)功率譜密度(相對值)與PREM 模型給出的自由振蕩頻率的對應(yīng)情況。該頻段內(nèi)雖然噪聲水平較大,但仍較清楚地檢測到基頻球型振型0S27—0S49以及諧頻球型振型3S13、3S16—3S18。在此頻段范圍內(nèi),觀測值與理論值的偏離與其他頻段比相對較大。

    3 計算結(jié)果對比分析

    此次日本9 級大地震釋放的能量巨大,常熟臺形變觀測儀不僅檢測到0S3—0S49全部基頻球型振蕩,還檢測到基頻環(huán)型振蕩及諧頻球型振蕩。統(tǒng)計4 組形變儀檢測的地球自由振蕩振型觀測頻率值和PREM 模型理論值的相對偏差,結(jié)果見表1、表2。

    表1 常熟臺4 組形變儀檢測的環(huán)型自由振蕩值與PREM 理論值及其偏差Table 1 Torsional free oscillation values detected by four groups of deformation instruments at Changshu station and PREM theoretical value and their deviations

    表2 常熟臺4 組形變儀檢測的球型自由振蕩值與PREM 理論值及其偏差Table 2 The spheroidal free oscillation values detected by four groups of deformation instruments at Changshu station and the theoretical values of PREM and their deviations

    續(xù)表2

    由表1、表2 可知:①水管儀檢測到45 個基頻球型振型、10 個諧頻球型振蕩、12 個基頻環(huán)型振型;②垂直擺檢測到36 個基頻球型振型、11 個諧頻球型振蕩、14 個基頻環(huán)型振型;③洞體應(yīng)變儀檢測到41 個基頻球型振型、3 個諧頻球型振蕩、9 個基頻環(huán)型振型;④體應(yīng)變儀檢測到40 個基頻球型振型、7 個諧頻球型振蕩、2 個基頻環(huán)型振型。

    4 組形變儀球型振型觀測頻率值與PREM 理論值整體平均偏差分別為:水管儀為0.112%,垂直擺為0.086%,洞體應(yīng)變儀為0.088%,體應(yīng)變儀為0.100%。

    常熟臺形變儀較清晰檢測到與地球內(nèi)部深部結(jié)構(gòu)有密切關(guān)系的低階基頻球型振型0S3、0S4,但觀測值與PREM 理論值相對偏差高于其他振型,可能存在振型譜線分解現(xiàn)象,能量分散到臨近幾個頻率上(姚家俊等,2012)。除0S3、0S4外,常熟臺4 組形變儀在低、中頻段檢測到的基頻球型振型譜峰清晰,與PREM 理論值相對偏差均低于0.22%(表1),低、中頻段觀測振型振幅明顯比高頻段基頻球型振蕩振幅高。但隨著頻率的增加,自由振蕩能量逐步衰減,在4.3×10-3Hz 以上頻段,部分振型譜峰雖能分辨,但并不突出(圖4),檢測的振型觀測頻率值與PREM 理論值相對偏差量明顯變大,信噪比較低(表1),如水管儀0S42振型觀測值與PREM 理論值偏差為0.339%。觀測頻率值與理論值偏差較大的原因可能是,高頻段的地球自由振蕩活動范圍隨著地球基型振蕩階數(shù)的增大趨近于表層,介質(zhì)的不均勻性和各向異性開始變的明顯,即正常觀測下的背景噪聲影響變高。

    與短基線形變儀相比,長基線水管儀、洞體應(yīng)變儀及深井下體應(yīng)變儀更容易受到地表高頻噪聲的干擾。常熟臺垂直擺在高頻段4.0×10-3—5.6×10-3Hz 噪聲水平較高,檢測到的自由振蕩振型較少,水管儀在0.4×10-3—5.6×10-3Hz 頻段背景噪聲干擾較小。由表1、表2 可知,常熟臺形變儀檢測的基頻環(huán)型振型、諧頻球型振型集中分布在1.0×10-3—4.3×10-3Hz 頻段,觀測振型譜峰均較清晰,信噪比>3,振型振幅小于基頻球型振型;基頻環(huán)型振型與PREM 理論值的相對偏差均低于0.27%,部分基頻環(huán)型振型與鄰近球型振型有不同程度的耦合現(xiàn)象;常熟臺形變儀檢測的諧頻球型振型除1S2、1S3相對偏差較大外,其余諧頻球型振型相對偏差小于0.20%。

    如圖5 所示,利用常熟臺日本震后5 天體應(yīng)變觀測數(shù)據(jù),采用小波變換方法,檢測到地震激發(fā)的2 個基頻環(huán)型振蕩:0T9、0T17振型,其觀測頻率值與PREM 理論值相對偏差分別為0.027%、0.041%(表2),信噪比SNR >3,其他3 組形變儀同樣檢測到2 個環(huán)型振蕩振型。

    圖5 TJ-2 體應(yīng)變功率譜密度估計Fig.5 Power spectral density estimation of the observation of TJ-2 volume strainmeter

    檢測結(jié)果顯示:①在傾斜觀測儀中,與垂直擺相比,水管儀檢測振型數(shù)量較多,但整體平均偏差較高。在應(yīng)變觀測儀中,與體應(yīng)變儀相比,洞體應(yīng)變儀檢測振型總數(shù)量較多,整體平均偏差較??;②在4 組形變儀中,水管儀檢測球型振蕩振型(含諧頻)能力最強,垂直擺檢測環(huán)型振蕩振型能力最強,從檢測自由振蕩振型總體數(shù)量來看,水管儀潮汐觀測精度最高、背景噪聲干擾最低,檢測的振型數(shù)量最多,垂直擺次之,體應(yīng)變儀最少,表明常熟臺傾斜觀測儀對自由振蕩信號的檢測能力優(yōu)于應(yīng)變觀測儀。

    4 結(jié)論

    利用常熟臺同一場地觀測的4 組形變觀測儀,檢測提取日本9 級大地震激發(fā)的地球自由振蕩信號,不僅準(zhǔn)確檢測到球型振蕩0S3—0S49全部47 個基頻振型和14 個諧頻振型,還檢測到環(huán)型振蕩0T5—0T25之間15 個基頻振型,與PREM 理論頻率值基本符合,檢測結(jié)果較好。結(jié)果顯示,常熟臺水管儀檢測球型振蕩振型能力最強,垂直擺檢測環(huán)型振蕩振型能力最強,且傾斜儀對自由振蕩信號的檢測能力優(yōu)于應(yīng)變儀。

    利用小波變換進行潮汐信號濾除處理,對形變觀測數(shù)據(jù)提取地球自由振蕩信號效果較好,可有效壓制觀測數(shù)據(jù)中的背景噪聲干擾,在此基礎(chǔ)上,在體應(yīng)變觀測數(shù)據(jù)中檢測到清晰、信噪聲比高的基頻環(huán)型振型。

    4 組形變觀測儀處于同一場地,所檢測振型的振幅隨著頻率的增大出現(xiàn)下降趨勢,證明自由振蕩信號傳播隨著頻率的增加衰減加快。同時,隨著自由振蕩階數(shù)的增加,與PREM 理論值存在偏差是因為地球介質(zhì)的橫向不均勻或各和向異性,對這些數(shù)據(jù)加以研究可以得到地球介質(zhì)的各向異性信息。

    地球自由振蕩頻率可以用來探測地球內(nèi)部結(jié)構(gòu),將觀測結(jié)果與PREM 模型理論值進行對比,研究PREM 模型與真實地球存在的差異,為約束和改善地球模型提供參考,以使其更接近于真實地球。

    猜你喜歡
    環(huán)型應(yīng)變儀球型
    CO2刺激響應(yīng)球型水凝膠的制備及在蛋白質(zhì)分離中的應(yīng)用
    體應(yīng)變儀傳感器更換前后數(shù)據(jù)關(guān)聯(lián)性分析
    后裝治療系統(tǒng)中使用環(huán)型施源器放射源駐留位置偏差值驗證
    TJ-2型體積式應(yīng)變儀故障檢測裝置
    淺談BIM技術(shù)應(yīng)用于球型網(wǎng)架空間坐標(biāo)的智能檢測
    3種鉆孔應(yīng)變儀記錄的瑪多MS7.4地震同震響應(yīng)特征對比分析
    中國地震(2021年3期)2022-01-12 09:42:00
    光纖應(yīng)變儀的研制與應(yīng)用前景
    一種環(huán)型橡膠條防凍測試用拉伸定位治具
    環(huán)型混凝土電桿生產(chǎn)工藝控制及關(guān)鍵技術(shù)
    價值工程(2017年11期)2017-04-18 05:41:23
    基于可操作性指標(biāo)的球型腕優(yōu)化*
    成年美女黄网站色视频大全免费| 1024视频免费在线观看| 亚洲精品自拍成人| 国产欧美亚洲国产| 国产精品偷伦视频观看了| 亚洲天堂av无毛| 桃花免费在线播放| 亚洲精品国产色婷婷电影| 日韩av免费高清视频| 黄色怎么调成土黄色| 美国免费a级毛片| 国产精品久久久久久精品电影小说| 久久精品久久久久久久性| 十八禁高潮呻吟视频| 久久久精品免费免费高清| 亚洲男人天堂网一区| 中文字幕另类日韩欧美亚洲嫩草| 制服诱惑二区| 亚洲精品国产一区二区精华液| 色视频在线一区二区三区| av电影中文网址| 少妇熟女欧美另类| av线在线观看网站| 欧美人与善性xxx| 熟女电影av网| 国产 一区精品| 久久久久视频综合| 美女福利国产在线| 午夜激情av网站| 亚洲精品中文字幕在线视频| 免费看不卡的av| 久久女婷五月综合色啪小说| 黄片无遮挡物在线观看| 色哟哟·www| xxxhd国产人妻xxx| 五月天丁香电影| 美女视频免费永久观看网站| 欧美在线黄色| 最近最新中文字幕大全免费视频 | 久久久久精品久久久久真实原创| 有码 亚洲区| 高清不卡的av网站| 少妇熟女欧美另类| 少妇精品久久久久久久| 日本欧美视频一区| 不卡视频在线观看欧美| 亚洲一码二码三码区别大吗| 精品国产一区二区久久| 精品第一国产精品| 2021少妇久久久久久久久久久| 欧美人与性动交α欧美软件| 乱人伦中国视频| 欧美精品高潮呻吟av久久| 99re6热这里在线精品视频| 女人久久www免费人成看片| 亚洲av日韩在线播放| 只有这里有精品99| 边亲边吃奶的免费视频| 日韩制服骚丝袜av| 国产成人av激情在线播放| 性色avwww在线观看| 麻豆乱淫一区二区| 午夜免费观看性视频| 久久久久久久久久久久大奶| 91aial.com中文字幕在线观看| 91aial.com中文字幕在线观看| 国语对白做爰xxxⅹ性视频网站| 丰满饥渴人妻一区二区三| 欧美+日韩+精品| 寂寞人妻少妇视频99o| 午夜福利视频在线观看免费| 女的被弄到高潮叫床怎么办| 成人影院久久| 中文字幕人妻丝袜制服| 午夜福利视频在线观看免费| 青春草国产在线视频| 丰满乱子伦码专区| 女人精品久久久久毛片| 午夜免费鲁丝| 在现免费观看毛片| 国产成人aa在线观看| 伊人久久国产一区二区| 国产午夜精品一二区理论片| 黄频高清免费视频| 日韩不卡一区二区三区视频在线| 欧美少妇被猛烈插入视频| 校园人妻丝袜中文字幕| 国产日韩欧美视频二区| 欧美人与善性xxx| 91国产中文字幕| 免费av中文字幕在线| 曰老女人黄片| 女人高潮潮喷娇喘18禁视频| 午夜免费观看性视频| 男女午夜视频在线观看| 秋霞伦理黄片| 制服丝袜香蕉在线| av片东京热男人的天堂| 深夜精品福利| 99精国产麻豆久久婷婷| 久久久精品国产亚洲av高清涩受| 两个人看的免费小视频| 久久久久精品性色| 少妇精品久久久久久久| 免费女性裸体啪啪无遮挡网站| 天天躁夜夜躁狠狠久久av| 91国产中文字幕| 熟女少妇亚洲综合色aaa.| 成人影院久久| 免费日韩欧美在线观看| 高清在线视频一区二区三区| 18禁观看日本| 精品一区二区三卡| 丝袜美腿诱惑在线| 亚洲精品自拍成人| 日韩av免费高清视频| av卡一久久| 热re99久久精品国产66热6| 久久久a久久爽久久v久久| 久久久久精品人妻al黑| 街头女战士在线观看网站| 亚洲精品国产一区二区精华液| 啦啦啦在线免费观看视频4| 国产亚洲最大av| 午夜福利影视在线免费观看| 国产成人精品婷婷| 久久精品久久精品一区二区三区| a 毛片基地| 久久精品国产鲁丝片午夜精品| 亚洲三级黄色毛片| 亚洲精品成人av观看孕妇| 一区福利在线观看| 久久精品人人爽人人爽视色| 久久99蜜桃精品久久| 又粗又硬又长又爽又黄的视频| 男男h啪啪无遮挡| 国产精品国产三级国产专区5o| 日韩电影二区| 亚洲人成77777在线视频| 午夜日本视频在线| 欧美中文综合在线视频| 欧美人与性动交α欧美软件| 涩涩av久久男人的天堂| 亚洲中文av在线| 亚洲成色77777| 老汉色av国产亚洲站长工具| 免费大片黄手机在线观看| 精品国产乱码久久久久久男人| 少妇人妻久久综合中文| 欧美日韩精品网址| 久久久久国产网址| 国产精品亚洲av一区麻豆 | 女人被躁到高潮嗷嗷叫费观| 欧美激情极品国产一区二区三区| 在线观看三级黄色| 麻豆精品久久久久久蜜桃| 日韩电影二区| 日韩av不卡免费在线播放| 亚洲av.av天堂| av在线播放精品| 在线观看www视频免费| 精品一区二区三卡| 亚洲国产精品一区二区三区在线| 欧美人与性动交α欧美精品济南到 | 免费观看在线日韩| 亚洲av.av天堂| 美女福利国产在线| 国产成人欧美| 两个人免费观看高清视频| 欧美日韩一区二区视频在线观看视频在线| 一区二区av电影网| 91精品伊人久久大香线蕉| 9191精品国产免费久久| 国产老妇伦熟女老妇高清| 亚洲美女搞黄在线观看| 性高湖久久久久久久久免费观看| 五月伊人婷婷丁香| 日日摸夜夜添夜夜爱| 国产亚洲一区二区精品| 亚洲成av片中文字幕在线观看 | 有码 亚洲区| 中文字幕精品免费在线观看视频| tube8黄色片| 久久久久久久久久人人人人人人| 国产成人91sexporn| 国产1区2区3区精品| 建设人人有责人人尽责人人享有的| 国产不卡av网站在线观看| 一级片免费观看大全| 国产精品人妻久久久影院| 亚洲成人手机| 新久久久久国产一级毛片| 国产男女超爽视频在线观看| 精品国产一区二区久久| 啦啦啦中文免费视频观看日本| 天堂俺去俺来也www色官网| 宅男免费午夜| 这个男人来自地球电影免费观看 | 人人妻人人澡人人爽人人夜夜| 观看av在线不卡| 成人18禁高潮啪啪吃奶动态图| 少妇 在线观看| 亚洲国产看品久久| 母亲3免费完整高清在线观看 | 久久久久久免费高清国产稀缺| 亚洲欧美成人精品一区二区| 亚洲美女黄色视频免费看| 人人妻人人添人人爽欧美一区卜| 丝瓜视频免费看黄片| 久久女婷五月综合色啪小说| 亚洲美女视频黄频| 女性被躁到高潮视频| 欧美人与善性xxx| 精品久久蜜臀av无| 午夜精品国产一区二区电影| 2018国产大陆天天弄谢| 美女视频免费永久观看网站| av有码第一页| 成年女人毛片免费观看观看9 | 欧美成人午夜精品| 欧美国产精品va在线观看不卡| 日日啪夜夜爽| 日本91视频免费播放| 男人添女人高潮全过程视频| 免费少妇av软件| 宅男免费午夜| 精品人妻熟女毛片av久久网站| 欧美日本中文国产一区发布| 深夜精品福利| 精品国产一区二区三区久久久樱花| 亚洲欧美精品自产自拍| 欧美精品一区二区免费开放| 国产精品成人在线| 老鸭窝网址在线观看| 看非洲黑人一级黄片| 日本猛色少妇xxxxx猛交久久| 国产成人免费无遮挡视频| 亚洲综合精品二区| 精品国产一区二区三区四区第35| 日本黄色日本黄色录像| 最黄视频免费看| 久久97久久精品| 久久人人97超碰香蕉20202| 国产男女超爽视频在线观看| 少妇人妻久久综合中文| 在线免费观看不下载黄p国产| 免费看不卡的av| 69精品国产乱码久久久| 美女福利国产在线| 亚洲欧美色中文字幕在线| 边亲边吃奶的免费视频| 久久久久人妻精品一区果冻| 十八禁高潮呻吟视频| 人体艺术视频欧美日本| 欧美中文综合在线视频| 成年av动漫网址| 99国产综合亚洲精品| 国产极品天堂在线| 青春草亚洲视频在线观看| 久热久热在线精品观看| 国产av码专区亚洲av| 久久精品亚洲av国产电影网| 精品国产国语对白av| 国产成人精品一,二区| 日韩av在线免费看完整版不卡| 免费观看性生交大片5| 性色av一级| 99久久中文字幕三级久久日本| 免费观看性生交大片5| 精品视频人人做人人爽| 精品一区二区免费观看| 欧美成人精品欧美一级黄| 成人毛片60女人毛片免费| 日日摸夜夜添夜夜爱| 日韩精品免费视频一区二区三区| 街头女战士在线观看网站| 日韩欧美精品免费久久| 欧美激情极品国产一区二区三区| 男女高潮啪啪啪动态图| 国产色婷婷99| 国产97色在线日韩免费| 久久精品国产a三级三级三级| 久久久亚洲精品成人影院| 王馨瑶露胸无遮挡在线观看| 亚洲成av片中文字幕在线观看 | 最近中文字幕高清免费大全6| 在线观看一区二区三区激情| 中文字幕人妻熟女乱码| 少妇 在线观看| 丝袜美腿诱惑在线| 亚洲国产精品一区三区| 人妻 亚洲 视频| 十分钟在线观看高清视频www| 青春草国产在线视频| 成年美女黄网站色视频大全免费| 性色avwww在线观看| 高清视频免费观看一区二区| 男女高潮啪啪啪动态图| 国产一区有黄有色的免费视频| 婷婷色综合大香蕉| 一区二区三区激情视频| 交换朋友夫妻互换小说| 亚洲精品国产av成人精品| 美女高潮到喷水免费观看| 男女国产视频网站| 成年动漫av网址| 久久久精品国产亚洲av高清涩受| 热re99久久精品国产66热6| 美女大奶头黄色视频| 大片免费播放器 马上看| 亚洲视频免费观看视频| 男女边吃奶边做爰视频| 久久久久久久亚洲中文字幕| 日韩一卡2卡3卡4卡2021年| 飞空精品影院首页| 亚洲熟女精品中文字幕| 国产亚洲一区二区精品| 最新的欧美精品一区二区| 日本欧美国产在线视频| 中文字幕另类日韩欧美亚洲嫩草| 一本久久精品| 国产免费现黄频在线看| 免费观看无遮挡的男女| 亚洲男人天堂网一区| 一级毛片电影观看| 精品人妻一区二区三区麻豆| 国产极品粉嫩免费观看在线| 久热久热在线精品观看| 国产精品偷伦视频观看了| 亚洲欧美一区二区三区国产| 国产免费福利视频在线观看| 99久久精品国产国产毛片| 成人亚洲精品一区在线观看| 午夜免费鲁丝| 高清视频免费观看一区二区| 国产一区二区 视频在线| 午夜福利乱码中文字幕| 在线观看三级黄色| 久久狼人影院| 午夜福利视频精品| 考比视频在线观看| av一本久久久久| 91成人精品电影| 午夜日本视频在线| 亚洲欧美精品综合一区二区三区 | 天美传媒精品一区二区| 欧美黄色片欧美黄色片| 成人18禁高潮啪啪吃奶动态图| 一级爰片在线观看| 99久久中文字幕三级久久日本| 丝瓜视频免费看黄片| 观看av在线不卡| 日韩中文字幕欧美一区二区 | 最近中文字幕2019免费版| 综合色丁香网| 精品少妇黑人巨大在线播放| 美女xxoo啪啪120秒动态图| 两个人免费观看高清视频| 美女国产视频在线观看| 久久精品国产a三级三级三级| 久久精品亚洲av国产电影网| 国产精品久久久久久精品古装| 亚洲激情五月婷婷啪啪| 97人妻天天添夜夜摸| 亚洲欧美精品自产自拍| 高清欧美精品videossex| 国产精品.久久久| 亚洲国产精品一区三区| 亚洲男人天堂网一区| 亚洲av电影在线观看一区二区三区| 青春草视频在线免费观看| 可以免费在线观看a视频的电影网站 | 欧美国产精品一级二级三级| 新久久久久国产一级毛片| 亚洲欧美成人综合另类久久久| 国产无遮挡羞羞视频在线观看| 国产精品国产三级专区第一集| 国产黄色视频一区二区在线观看| 久久免费观看电影| 免费黄色在线免费观看| 性高湖久久久久久久久免费观看| 最近中文字幕高清免费大全6| 久久久久精品人妻al黑| 久久热在线av| 18禁国产床啪视频网站| 久久午夜综合久久蜜桃| 精品少妇黑人巨大在线播放| av一本久久久久| kizo精华| 久久久久人妻精品一区果冻| 国产毛片在线视频| 精品少妇久久久久久888优播| 欧美+日韩+精品| 9色porny在线观看| 中文字幕人妻丝袜一区二区 | 国产精品国产三级国产专区5o| 美女大奶头黄色视频| 日韩 亚洲 欧美在线| 国产成人精品无人区| 2021少妇久久久久久久久久久| 制服丝袜香蕉在线| 亚洲欧美中文字幕日韩二区| 一级,二级,三级黄色视频| 国产精品久久久久久精品电影小说| 国产成人精品久久二区二区91 | 日韩一区二区三区影片| 国产极品天堂在线| 三级国产精品片| 一区二区三区精品91| 韩国高清视频一区二区三区| 黄片无遮挡物在线观看| 成年av动漫网址| 亚洲伊人久久精品综合| 波多野结衣av一区二区av| 亚洲一级一片aⅴ在线观看| av在线app专区| 99久久中文字幕三级久久日本| 免费黄频网站在线观看国产| 自拍欧美九色日韩亚洲蝌蚪91| 欧美在线黄色| 一级毛片电影观看| 欧美+日韩+精品| 久久青草综合色| 欧美激情极品国产一区二区三区| 午夜老司机福利剧场| 黑丝袜美女国产一区| 国语对白做爰xxxⅹ性视频网站| 国产高清不卡午夜福利| 成年人午夜在线观看视频| 欧美精品一区二区大全| 亚洲欧美中文字幕日韩二区| 亚洲精品,欧美精品| 曰老女人黄片| freevideosex欧美| 国产精品蜜桃在线观看| 国产国语露脸激情在线看| 久久久久久人妻| 精品一区二区三区四区五区乱码 | 国产精品久久久久久av不卡| 亚洲中文av在线| 午夜日韩欧美国产| 黄频高清免费视频| 免费黄色在线免费观看| av在线观看视频网站免费| 另类精品久久| 99热国产这里只有精品6| 亚洲图色成人| 日韩大片免费观看网站| 日韩av不卡免费在线播放| 在现免费观看毛片| 一级,二级,三级黄色视频| 人妻系列 视频| 久久精品亚洲av国产电影网| 亚洲精品国产av蜜桃| 日韩制服骚丝袜av| av在线老鸭窝| 欧美另类一区| 深夜精品福利| 毛片一级片免费看久久久久| 丝袜美腿诱惑在线| 狂野欧美激情性bbbbbb| 亚洲精品,欧美精品| 中文字幕亚洲精品专区| 亚洲欧美成人综合另类久久久| 精品亚洲成a人片在线观看| 免费不卡的大黄色大毛片视频在线观看| 精品一品国产午夜福利视频| 国产精品一区二区在线不卡| 国产探花极品一区二区| 9191精品国产免费久久| 午夜福利网站1000一区二区三区| 国产精品国产三级国产专区5o| 一区在线观看完整版| 日韩制服丝袜自拍偷拍| 高清av免费在线| 亚洲精华国产精华液的使用体验| 免费观看性生交大片5| 亚洲av电影在线进入| 国产老妇伦熟女老妇高清| 一级片免费观看大全| 99热全是精品| 最新中文字幕久久久久| 老女人水多毛片| 国产激情久久老熟女| 国产午夜精品一二区理论片| 免费观看性生交大片5| 成人免费观看视频高清| 午夜91福利影院| 国产高清不卡午夜福利| 下体分泌物呈黄色| 亚洲美女黄色视频免费看| 中文字幕制服av| 久久久久视频综合| 18+在线观看网站| 国产av国产精品国产| 久久久国产精品麻豆| 亚洲成人手机| 国产免费福利视频在线观看| 日韩中文字幕欧美一区二区 | 99精国产麻豆久久婷婷| 中文乱码字字幕精品一区二区三区| 丝袜人妻中文字幕| 国产精品蜜桃在线观看| 各种免费的搞黄视频| 亚洲人成网站在线观看播放| 国产1区2区3区精品| 免费高清在线观看日韩| 国产老妇伦熟女老妇高清| 国产一区二区三区av在线| 乱人伦中国视频| √禁漫天堂资源中文www| 在线精品无人区一区二区三| 夜夜骑夜夜射夜夜干| 国产极品天堂在线| 日日啪夜夜爽| 中文字幕最新亚洲高清| 老司机影院毛片| 久久久久精品性色| 国产精品一国产av| 伦精品一区二区三区| 午夜日韩欧美国产| av网站免费在线观看视频| 晚上一个人看的免费电影| 精品午夜福利在线看| 中文字幕人妻丝袜一区二区 | 热99国产精品久久久久久7| 丝袜人妻中文字幕| 九九爱精品视频在线观看| 国产成人精品婷婷| 亚洲在久久综合| 精品国产国语对白av| 制服丝袜香蕉在线| 三级国产精品片| 最新中文字幕久久久久| 精品人妻在线不人妻| 久久 成人 亚洲| 久久99一区二区三区| 久久久久久免费高清国产稀缺| 国产成人91sexporn| 久久精品熟女亚洲av麻豆精品| 久久综合国产亚洲精品| 人体艺术视频欧美日本| 建设人人有责人人尽责人人享有的| 精品国产乱码久久久久久小说| 免费黄频网站在线观看国产| 青春草亚洲视频在线观看| 亚洲久久久国产精品| 国产成人a∨麻豆精品| 亚洲精品aⅴ在线观看| 国产一级毛片在线| 少妇的丰满在线观看| 久久av网站| 午夜老司机福利剧场| 婷婷色麻豆天堂久久| 丝瓜视频免费看黄片| 91久久精品国产一区二区三区| 看十八女毛片水多多多| 9热在线视频观看99| 五月伊人婷婷丁香| 欧美在线黄色| 免费观看性生交大片5| 大话2 男鬼变身卡| 妹子高潮喷水视频| 26uuu在线亚洲综合色| tube8黄色片| 青青草视频在线视频观看| 亚洲欧美精品自产自拍| 久久久久久人人人人人| 亚洲综合色网址| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 男人添女人高潮全过程视频| 免费观看性生交大片5| 精品福利永久在线观看| 黄色 视频免费看| 亚洲综合色网址| 五月开心婷婷网| 在线免费观看不下载黄p国产| 菩萨蛮人人尽说江南好唐韦庄| 久久毛片免费看一区二区三区| 亚洲婷婷狠狠爱综合网| 国产精品偷伦视频观看了| 免费少妇av软件| 免费观看性生交大片5| 母亲3免费完整高清在线观看 | 国产黄频视频在线观看| 亚洲av综合色区一区| 精品久久蜜臀av无| 国产在线视频一区二区| 国产黄色视频一区二区在线观看| 99热网站在线观看| 91精品三级在线观看| 9热在线视频观看99| 性少妇av在线| 99国产精品免费福利视频| 久久久久网色| 久久久久精品久久久久真实原创| 国产又色又爽无遮挡免| 午夜福利影视在线免费观看| 欧美人与性动交α欧美精品济南到 | 一区二区日韩欧美中文字幕| 亚洲成人手机| 街头女战士在线观看网站| 国产精品香港三级国产av潘金莲 | 久久久久久久精品精品| 最近最新中文字幕大全免费视频 | 青春草亚洲视频在线观看| 香蕉丝袜av| 一级a爱视频在线免费观看| 国产黄频视频在线观看| 亚洲av在线观看美女高潮| 人人澡人人妻人| 亚洲欧美一区二区三区黑人 | 亚洲综合色网址| 色94色欧美一区二区| 国产乱来视频区|