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

    利用湖州臺(tái)形變觀測(cè)資料檢測(cè)地震激發(fā)的地球球形自由振蕩

    2015-09-03 09:33:56方燕勛卞根發(fā)惠若愚孫宗強(qiáng)
    中國(guó)地震 2015年3期
    關(guān)鍵詞:應(yīng)變儀傾斜儀球型

    方燕勛 卞根發(fā) 惠若愚 孫宗強(qiáng)

    1)浙江省地震局湖州地震臺(tái),浙江省湖州市峴山路15號(hào) 313000

    2)浙江省地震局,杭州 310013

    3)云南省水富縣防震減災(zāi)局,云南水富 657800

    0 引言

    大地震能激發(fā)產(chǎn)生涉及地球局部運(yùn)動(dòng)的體波和面波,也能激發(fā)全球規(guī)模的地球自由振蕩。巨大遠(yuǎn)震的震時(shí)、震后變化在形變觀測(cè)資料中反映為階變或振蕩型波動(dòng)現(xiàn)象的疊加。最早采用應(yīng)變觀測(cè)方法對(duì)地球自由振蕩進(jìn)行準(zhǔn)確測(cè)量是在1960年智利MS8.3地震后,由Benioff等(1961)用Isabella應(yīng)變儀實(shí)現(xiàn)的,因其觀測(cè)結(jié)果與理論值一致性良好,故此確定了長(zhǎng)周期自由振蕩的存在。近年來(lái),一些學(xué)者根據(jù)傾斜或應(yīng)變觀測(cè)資料深入研究了地球自由振蕩,邱澤華等(2007)利用泰安地震臺(tái)鉆孔差應(yīng)變儀觀測(cè)資料、唐磊等(2007)利用中國(guó)鉆孔體應(yīng)變臺(tái)網(wǎng)觀測(cè)資料研究了2004年12月26日蘇門(mén)答臘巨大地震激發(fā)的地球自由振蕩,Bolt等(1975)、Zadro等(1999)根據(jù)應(yīng)變儀和傾斜儀的觀測(cè)結(jié)果對(duì)自由振蕩進(jìn)行了研究,任佳等(2009)利用寬城臺(tái)垂直擺傾斜儀和涉縣臺(tái)水管傾斜儀的觀測(cè)資料提取了汶川地震激發(fā)的地球球型振蕩信息。目前已有的研究結(jié)果表明,能檢測(cè)到自由振蕩的形變觀測(cè)儀主要有體應(yīng)變儀、水管傾斜儀、伸縮儀、垂直擺傾斜儀等。通過(guò)檢索以往相關(guān)文獻(xiàn)發(fā)現(xiàn),大多數(shù)以一個(gè)臺(tái)站的一種觀測(cè)資料進(jìn)行自由振蕩分析,任佳等(2009)雖然同時(shí)使用了垂直擺和水管儀的觀測(cè)資料提取汶川地震激發(fā)的球型振蕩,但使用的是不同臺(tái)站的資料,并沒(méi)有對(duì)同一臺(tái)站的多種傾斜或應(yīng)變儀的資料進(jìn)行對(duì)比分析。

    2011年3月11日,日本發(fā)生9.0級(jí)地震。該地震激發(fā)了全球規(guī)模的地球自由振蕩,這也為利用同一臺(tái)站的不同地傾斜或應(yīng)變儀觀測(cè)資料對(duì)同一地球自由振蕩進(jìn)行檢測(cè)與對(duì)比分析提供了難得的機(jī)遇。本文重點(diǎn)對(duì)浙江省湖州地震臺(tái)(后簡(jiǎn)稱(chēng)湖州臺(tái))工作狀態(tài)良好的垂直擺傾斜儀、水管傾斜儀、伸縮儀和體應(yīng)變儀針對(duì)此次地震的觀測(cè)資料進(jìn)行分析,目的是通過(guò)同一臺(tái)站不同形變儀器觀測(cè)的資料,進(jìn)一步證明利用傾斜和形變的觀測(cè)資料分析地震激發(fā)地球自由振蕩的可行性,并分析不同儀器對(duì)地球自由振蕩的檢測(cè)能力和差異特征。

    1 資料

    1.1 資料選取

    地殼形變觀測(cè)具有明確的幾何意義和地球物理意義(國(guó)家地震局科技監(jiān)測(cè)司,1995)。中國(guó)地震局武漢地震科學(xué)儀器研究院研制的VS型垂直擺傾斜儀、DSQ型水管傾斜儀、SSY型伸縮儀和中國(guó)地震局地殼應(yīng)力研究所研制的TJ-2型體應(yīng)變儀是目前國(guó)內(nèi)地傾斜、應(yīng)變潮汐觀測(cè)臺(tái)網(wǎng)中使用的主要觀測(cè)儀器。

    湖州臺(tái)水管傾斜儀、垂直擺傾斜儀和伸縮儀于2002年投入正常觀測(cè),體應(yīng)變儀于2008年投入正常觀測(cè)。這些觀測(cè)儀器均安裝在完整性極好的花崗巖上,能清晰記錄到潮汐波和同震響應(yīng)(方燕勛等,2013)。筆者從4套形變觀測(cè)儀所產(chǎn)出的數(shù)據(jù)中提取了日本9.0級(jí)地震激發(fā)的地球球型自由振蕩的振型,采用地震發(fā)生后受地震影響明顯的5天(共7200min)觀測(cè)數(shù)據(jù)進(jìn)行對(duì)比分析。圖1給出了湖州臺(tái)垂直擺NS、EW向,水管儀NS、EW向,伸縮儀NS、EW向和應(yīng)變儀記錄的日本9.0級(jí)大地震后5天的分鐘值曲線(xiàn),數(shù)據(jù)未經(jīng)任何消除干擾處理。

    圖1 湖州臺(tái)地傾斜、應(yīng)變儀記錄到日本大地震的分鐘值曲線(xiàn)

    1.2 觀測(cè)資料預(yù)處理

    形變觀測(cè)儀運(yùn)行時(shí)的觀測(cè)背景噪聲非常復(fù)雜,主要包括兩部分:①儀器運(yùn)行時(shí)儀器本身的觀測(cè)噪聲;②臺(tái)站周?chē)皡^(qū)域性的背景噪聲等。為了消除觀測(cè)背景噪聲的影響,我們作了簡(jiǎn)單的背景噪聲去除處理。即采用大地震前后的一段地震平靜且儀器正常運(yùn)行時(shí)的觀測(cè)殘差資料,來(lái)模擬資料的正常觀測(cè)背景噪聲,再對(duì)其進(jìn)行傅氏分析最終得到觀測(cè)背景噪聲譜。因這次日本地震震級(jí)巨大且大于6.0級(jí)的余震較多,震后較長(zhǎng)時(shí)間仍不平靜,所以我們采用震前平靜時(shí)的觀測(cè)殘差資料模擬觀測(cè)儀的背景噪聲,對(duì)其進(jìn)行傅氏分析,獲得觀測(cè)背景噪聲譜。

    數(shù)據(jù)預(yù)處理后,我們采用通常的功率譜密度估計(jì)方法來(lái)區(qū)分球型自由振蕩振型的頻率(萬(wàn)永革,2007),其定義為記錄到波形資料的自相關(guān)函數(shù)Rn的傅里葉變換,即

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

    式中,N=7200為所采用數(shù)據(jù)的個(gè)數(shù),x為湖州臺(tái)去背景噪聲后的地傾斜或應(yīng)變數(shù)據(jù),Sk離散值為功率譜密度值。其所對(duì)應(yīng)的頻率值運(yùn)用觀測(cè)數(shù)據(jù)的采樣率60/s乘以序列號(hào)k確定。最后可得到功率譜密度估計(jì)。為了消除數(shù)據(jù)不能無(wú)限長(zhǎng)而必須加窗造成的對(duì)功率譜密度估計(jì)的影響,本文采用漢寧窗來(lái)抑制旁瓣,以突出主瓣。

    2 地球自由振蕩計(jì)算

    按照此前的計(jì)算方法,筆者得到了湖州臺(tái)4套形變觀測(cè)儀觀測(cè)資料的功率譜密度估計(jì)(圖 2~5)。

    圖2 0.2~1.3MHz頻段形變觀測(cè)數(shù)據(jù)功率譜密度估計(jì)曲線(xiàn)

    圖3 1.3~4.5MHz頻段形變觀測(cè)數(shù)據(jù)功率譜密度估計(jì)曲線(xiàn)

    圖4 4.5~6.0MHz頻段形變觀測(cè)數(shù)據(jù)功率譜密度估計(jì)曲線(xiàn)

    圖2為0.2×10-3~1.3×10-3Hz頻段的功率譜密度估計(jì),由圖2可見(jiàn),利用湖州臺(tái)形變觀測(cè)儀所記錄的數(shù)據(jù),除了能檢測(cè)到球型振蕩的基型振蕩0S5、0S6、0S7,外,還能清楚地檢測(cè)到1S3、2S3。但1S3、2S3振型在SSY、TJ-2兩套應(yīng)變儀中并不能明顯地從噪聲背景中分離出來(lái)。SSY、TJ-2兩套應(yīng)變儀可分辨 0.8MHz以下的頻段內(nèi)的0S2、0S3、0S4振型。由此看來(lái),運(yùn)用湖州臺(tái)形變觀測(cè)儀所記錄的數(shù)據(jù),能準(zhǔn)確檢測(cè)出大多數(shù)地球球型振蕩的基本振型,且能檢測(cè)到信號(hào)較強(qiáng)的諧頻成分。為此,作為初步研究結(jié)果,我們只給出PREM模型的基型球型振蕩的頻率與2011年日本9.0級(jí)地震后湖州臺(tái)形變觀測(cè)數(shù)據(jù)的功率譜密度進(jìn)行對(duì)比。

    圖5 6.0~8.0MHz頻段形變觀測(cè)數(shù)據(jù)功率譜密度估計(jì)曲線(xiàn)

    圖3為1.3×10-3~4.5×10-3Hz頻段范圍內(nèi)湖州臺(tái)形變觀測(cè)數(shù)據(jù)的功率譜密度(相對(duì)值)與PREM模型的球型基型振蕩頻率的對(duì)應(yīng)情況。在這一頻段范圍內(nèi)可看到,由湖州臺(tái)形變觀測(cè)數(shù)據(jù)可清楚地檢測(cè)到球型振蕩的基型振蕩0S8~0S37。但0S29和0S36的頻率值呈少許偏離。

    圖4為4.5×10-3~6.0×10-3Hz頻段范圍內(nèi)湖州臺(tái)形變觀測(cè)數(shù)據(jù)的功率譜密度(相對(duì)值)與PREM模型的球型基型振蕩0S38~0S54頻率值的對(duì)應(yīng)情況??煽吹胶菖_(tái)形變數(shù)據(jù)均能準(zhǔn)確反映地球自由振蕩0S38~0S50的信息,但0S51~0S54頻段范圍內(nèi)的譜峰并不明顯,較難分辨。

    圖5為6.0×10-3~8.0×10-3Hz頻段范圍的 PREM模型給出的地球球型自由振蕩振型頻率與2011年日本大地震形變觀測(cè)數(shù)據(jù)功率譜密度的對(duì)應(yīng)關(guān)系。湖州臺(tái)形變觀測(cè)數(shù)據(jù)有較豐富的諧頻成分,但大部分譜峰尚不能與理論計(jì)算值相對(duì)應(yīng),只有0S60振型與PREM模型基型球型振蕩頻率成分對(duì)應(yīng)較好。

    3 計(jì)算結(jié)果及不同觀測(cè)項(xiàng)的對(duì)比分析

    此次日本大地震中,湖州臺(tái)4套形變觀測(cè)儀所記錄的數(shù)據(jù)均能檢測(cè)到0S5~0S50地球自由振蕩振型(表1)。通常淺源地震主要激發(fā)基頻地球自由振蕩,但此次日本地震釋放出的能量巨大,不僅激發(fā)了基頻地球自由振蕩,還激發(fā)了諧頻地球自由振蕩。我們檢測(cè)到了2個(gè)諧頻振型,并與PREM地球模型的理論預(yù)測(cè)值作比較,發(fā)現(xiàn)實(shí)測(cè)值和理論值基本一致,表明對(duì)這2個(gè)振型的檢測(cè)是可靠的(表1)。此次日本地震激發(fā)的振型只有在1級(jí)、2級(jí)諧頻較好地檢測(cè)到,表明日本地震激發(fā)的振型并不是各級(jí)振型都能檢測(cè)到的,僅有某些種類(lèi)的諧頻振型在激發(fā)上占優(yōu)勢(shì),可被觀測(cè)儀器檢測(cè)到。

    為了便于對(duì)比分析湖州臺(tái)4套形變觀測(cè)儀記錄資料檢測(cè)到的地球自由振蕩振型與PREM模型理論值偏差,我們計(jì)算了其每一振型與理論值的相對(duì)偏差η,圖6中繪出4組觀測(cè)值在0S2~0S50振型與理論值的相對(duì)偏差(其中 VS、DSQ未檢測(cè)到0S2、0S3和0S4振型,該振型的相對(duì)偏差數(shù)據(jù)為空)。VS型垂直擺傾斜儀觀測(cè)數(shù)據(jù)檢測(cè)的振型與理論值相對(duì)偏差均小于0.83‰,其相對(duì)偏差較大的振型是0S31、0S35;DSQ型水管傾斜儀觀測(cè)數(shù)據(jù)檢測(cè)的振型與理論值相對(duì)偏差均小于1.2‰,其相對(duì)偏差較大的振型是0S29、0S36;SSY型伸縮儀觀測(cè)數(shù)據(jù)檢測(cè)的振型與理論值相對(duì)偏差均小于0.89‰,其相對(duì)偏差較大的振型是0S37;TJ-2型體應(yīng)變儀觀測(cè)數(shù)據(jù)檢測(cè)的振型與理論值相對(duì)偏差均小于0.54‰。湖州臺(tái)4套形變觀測(cè)儀記錄資料檢測(cè)到的地球自由振蕩振型與PREM模型理論值偏差較大的振型集中在0S29~0S37。

    表1 湖州臺(tái)4套形變儀所記錄到的地球基頻球型振蕩值、PREM模型理論值

    圖6 湖州臺(tái)形變觀測(cè)數(shù)據(jù)檢測(cè)結(jié)果和PREM模型理論值相對(duì)偏差

    將地傾斜觀測(cè)儀檢測(cè)到的地球自由振蕩振型與理論值相對(duì)偏差對(duì)比分析后發(fā)現(xiàn),DSQ型水管傾斜儀高于VS型垂直擺傾斜儀。從地應(yīng)變觀測(cè)儀檢測(cè)到的地球自由振蕩振型與理論值相對(duì)偏差對(duì)比分析發(fā)現(xiàn),SSY型伸縮儀高于TJ-2型體應(yīng)變儀。這是由于兩者的工作環(huán)境、工作條件不同所決定的,一個(gè)是數(shù)十米的大跨度,一個(gè)是數(shù)十毫米的小跨度。它們獲得的信號(hào)變幅相差幾個(gè)量級(jí),且地殼應(yīng)變儀檢測(cè)到的地球自由振蕩振型與理論值相對(duì)偏差小于地殼傾斜觀測(cè)儀檢測(cè)到的。其中SSY、TJ-2這兩套應(yīng)變儀能清晰地檢測(cè)到與地球內(nèi)部的深部結(jié)構(gòu)有密切關(guān)系的低階振型:0S2、0S3和0S4。而VS、DSQ這兩套傾斜儀未檢測(cè)到這幾個(gè)振型,是由于正常運(yùn)行時(shí)的觀測(cè)背景噪聲在該頻段內(nèi)功率譜密度很高,地球自由振蕩波被淹沒(méi)了。而0S29和0S36的頻率值呈少許偏離,一方面是由于臺(tái)站可能在振型的節(jié)線(xiàn)附近,信噪比較??;另一方面也可能是由于這個(gè)高階基頻振型對(duì)上地幔的剪切波速較敏感。

    4 結(jié)論和討論

    (1)本研究利用同一臺(tái)站的垂直擺傾斜儀、水管傾斜儀、伸縮儀和體應(yīng)變儀的觀測(cè)資料成功地對(duì)日本9.0級(jí)大地震激發(fā)的球型自由振蕩信息進(jìn)行了提取,準(zhǔn)確檢測(cè)到了0S5~0S50地球自由振蕩振型,其中SSY、TJ-2這兩套應(yīng)變儀能清晰檢測(cè)到與地球內(nèi)部的深部結(jié)構(gòu)有密切關(guān)系的低階振型(0S2、0S3和0S4)。VS、DSQ檢測(cè)到 2個(gè)諧頻振型(1S3、2S3),且與 PREM預(yù)測(cè)的振蕩頻率基本符合。這既反映了形變觀測(cè)儀器運(yùn)行狀態(tài)良好,也表明本文的研究路徑可行。

    (2)在125~176s頻段內(nèi)(即大于0S51的頻域),我們觀測(cè)到較豐富的基頻振型,同時(shí)也注意到一些譜峰尚不能與理論計(jì)算值相對(duì)應(yīng)。這一方面與地幔的非均勻性有關(guān),另一方面可能是受傾斜或應(yīng)變觀測(cè)儀幅頻特征和采樣率的限制,不能記錄更高頻率的緣故。

    (3)湖州臺(tái)TJ-2型體應(yīng)變儀觀測(cè)數(shù)據(jù)對(duì)日本9.0級(jí)大地震激發(fā)的球型自由振蕩信息檢測(cè)結(jié)果與理論值的相對(duì)偏差最小,這可能與該型號(hào)儀器安裝在井下,受地表影響小有關(guān)(邱澤華等,2004)。

    (4)隨著對(duì)地球自由振蕩問(wèn)題探討的逐漸深入,除了能為地球結(jié)構(gòu)問(wèn)題的研究提供約束條件外,還可以利用自由振蕩觀測(cè)資料研究震源過(guò)程、鑒別慢地震(Beroza et al,1990)。近來(lái)的研究結(jié)果表明,地球的自由振蕩是一直存在的,其激發(fā)機(jī)制可能與大氣和海洋的運(yùn)動(dòng)有關(guān)(Suda et al,1998;Tanimoto et al,1998)。另外,如何用這些觀測(cè)資料來(lái)檢驗(yàn)傾斜或應(yīng)變儀的頻響特性,也是值得進(jìn)一步研究的問(wèn)題。

    致謝:本項(xiàng)研究得到萬(wàn)永革、任佳等老師的指點(diǎn),PREM模型數(shù)據(jù)由萬(wàn)永革老師提供,匿名審稿人對(duì)論文的修改提出了詳細(xì)而中肯的建議,在此一并表示衷心的感謝。

    猜你喜歡
    應(yīng)變儀傾斜儀球型
    CO2刺激響應(yīng)球型水凝膠的制備及在蛋白質(zhì)分離中的應(yīng)用
    體應(yīng)變儀傳感器更換前后數(shù)據(jù)關(guān)聯(lián)性分析
    TJ-2型體積式應(yīng)變儀故障檢測(cè)裝置
    淺談BIM技術(shù)應(yīng)用于球型網(wǎng)架空間坐標(biāo)的智能檢測(cè)
    3種鉆孔應(yīng)變儀記錄的瑪多MS7.4地震同震響應(yīng)特征對(duì)比分析
    安西地震臺(tái)洞體長(zhǎng)短基線(xiàn)傾斜儀同震響應(yīng)對(duì)比分析
    高原地震(2021年2期)2021-10-08 01:30:00
    光纖應(yīng)變儀的研制與應(yīng)用前景
    宜昌地震臺(tái)VP 型寬頻帶傾斜儀和DSQ 型水管傾斜儀同震響應(yīng)對(duì)比分析
    融合傾斜儀數(shù)據(jù)的盾構(gòu)姿態(tài)嚴(yán)密解算模型
    基于可操作性指標(biāo)的球型腕優(yōu)化*
    欧美日韩黄片免| 国产极品粉嫩免费观看在线| 日韩大尺度精品在线看网址 | 国产成人精品在线电影| 老司机午夜十八禁免费视频| 欧美大码av| 欧美国产精品va在线观看不卡| 波多野结衣高清无吗| 亚洲av成人av| 久久天躁狠狠躁夜夜2o2o| 可以免费在线观看a视频的电影网站| 中文字幕高清在线视频| 久久欧美精品欧美久久欧美| 国产精品自产拍在线观看55亚洲| 国内精品久久久久久久电影| 丁香六月欧美| 97碰自拍视频| 亚洲av日韩精品久久久久久密| 欧美成狂野欧美在线观看| 日本在线视频免费播放| 亚洲中文字幕日韩| 岛国视频午夜一区免费看| 日韩高清综合在线| 欧美一级a爱片免费观看看 | 欧美不卡视频在线免费观看 | 亚洲专区中文字幕在线| 国产精品二区激情视频| 欧美人与性动交α欧美精品济南到| 欧美大码av| 又紧又爽又黄一区二区| 午夜福利一区二区在线看| 中文字幕另类日韩欧美亚洲嫩草| 桃色一区二区三区在线观看| 欧美国产日韩亚洲一区| 国产成+人综合+亚洲专区| 搞女人的毛片| 免费人成视频x8x8入口观看| 免费久久久久久久精品成人欧美视频| 国语自产精品视频在线第100页| 免费少妇av软件| 精品免费久久久久久久清纯| 一卡2卡三卡四卡精品乱码亚洲| 国产精品亚洲美女久久久| 国产精品影院久久| 国产精品av久久久久免费| 如日韩欧美国产精品一区二区三区| 日韩欧美国产在线观看| 亚洲一区二区三区色噜噜| 日韩欧美三级三区| 国产极品粉嫩免费观看在线| 欧美乱码精品一区二区三区| 国产三级在线视频| 国产成人精品久久二区二区91| 可以免费在线观看a视频的电影网站| 无限看片的www在线观看| 村上凉子中文字幕在线| 国产成人精品久久二区二区91| 大码成人一级视频| 侵犯人妻中文字幕一二三四区| 久久精品aⅴ一区二区三区四区| 母亲3免费完整高清在线观看| 亚洲国产中文字幕在线视频| 精品国产亚洲在线| 亚洲av成人av| 97超级碰碰碰精品色视频在线观看| 99久久99久久久精品蜜桃| www.自偷自拍.com| 免费观看人在逋| 国产成人欧美在线观看| 色老头精品视频在线观看| 久久国产精品男人的天堂亚洲| 久久精品国产亚洲av高清一级| 成年女人毛片免费观看观看9| 老司机深夜福利视频在线观看| 久久久国产欧美日韩av| 日韩欧美免费精品| 操出白浆在线播放| 国产成+人综合+亚洲专区| 国产成人系列免费观看| 一级,二级,三级黄色视频| 久久久久精品国产欧美久久久| 国产精品99久久99久久久不卡| 两个人免费观看高清视频| 亚洲国产欧美一区二区综合| 最近最新中文字幕大全电影3 | 黑人巨大精品欧美一区二区蜜桃| 免费看美女性在线毛片视频| 国产精品久久电影中文字幕| 亚洲色图av天堂| 国产蜜桃级精品一区二区三区| 十八禁人妻一区二区| 级片在线观看| 日韩免费av在线播放| 午夜福利免费观看在线| 一卡2卡三卡四卡精品乱码亚洲| 国产精品一区二区精品视频观看| 亚洲成av片中文字幕在线观看| 99在线视频只有这里精品首页| 国产精品野战在线观看| 搡老岳熟女国产| 国产精品久久久久久人妻精品电影| 精品久久久久久成人av| 国语自产精品视频在线第100页| 他把我摸到了高潮在线观看| 真人一进一出gif抽搐免费| 中文字幕久久专区| 99国产极品粉嫩在线观看| 亚洲性夜色夜夜综合| 久久久久久亚洲精品国产蜜桃av| 国产精品99久久99久久久不卡| 宅男免费午夜| bbb黄色大片| 搡老熟女国产l中国老女人| 无遮挡黄片免费观看| 两个人视频免费观看高清| 亚洲国产精品久久男人天堂| 国产精品综合久久久久久久免费 | 91老司机精品| 老司机靠b影院| 国产精品爽爽va在线观看网站 | 高清在线国产一区| 韩国精品一区二区三区| 亚洲五月天丁香| 国产aⅴ精品一区二区三区波| 亚洲中文字幕一区二区三区有码在线看 | 在线观看午夜福利视频| 热re99久久国产66热| av中文乱码字幕在线| 一夜夜www| 国内久久婷婷六月综合欲色啪| 999久久久精品免费观看国产| 日韩三级视频一区二区三区| 欧美一级毛片孕妇| 亚洲欧美激情在线| 露出奶头的视频| 久久精品国产99精品国产亚洲性色 | 丁香欧美五月| 午夜精品在线福利| 欧美日韩瑟瑟在线播放| 欧美黑人欧美精品刺激| 亚洲第一欧美日韩一区二区三区| 精品不卡国产一区二区三区| 黄色视频,在线免费观看| 亚洲欧美精品综合久久99| 精品少妇一区二区三区视频日本电影| 久久久精品欧美日韩精品| 久久久久国产一级毛片高清牌| 黄色成人免费大全| 国产精品爽爽va在线观看网站 | 黑人欧美特级aaaaaa片| 国产三级在线视频| 亚洲精品久久国产高清桃花| 丁香六月欧美| av视频免费观看在线观看| 国产熟女午夜一区二区三区| 色婷婷久久久亚洲欧美| 99在线人妻在线中文字幕| 亚洲人成伊人成综合网2020| 在线十欧美十亚洲十日本专区| 国产私拍福利视频在线观看| 欧美色欧美亚洲另类二区 | 一区在线观看完整版| 欧美乱码精品一区二区三区| 久久中文看片网| 一进一出抽搐动态| 欧美成人午夜精品| 日本免费a在线| 国产精品 欧美亚洲| 在线观看www视频免费| 久久精品人人爽人人爽视色| 久久午夜综合久久蜜桃| 国产精品爽爽va在线观看网站 | 淫妇啪啪啪对白视频| 婷婷精品国产亚洲av在线| 69av精品久久久久久| 我的亚洲天堂| 啦啦啦免费观看视频1| 国产色视频综合| 黑人操中国人逼视频| 日本 av在线| 一区在线观看完整版| 中文字幕精品免费在线观看视频| 天堂√8在线中文| 狂野欧美激情性xxxx| 亚洲中文字幕日韩| 两个人视频免费观看高清| 人人澡人人妻人| 精品人妻1区二区| 一级,二级,三级黄色视频| 中文字幕av电影在线播放| 搡老岳熟女国产| 国产色视频综合| 咕卡用的链子| 亚洲av成人av| 精品久久久久久成人av| 免费看a级黄色片| 中国美女看黄片| 午夜日韩欧美国产| www.999成人在线观看| e午夜精品久久久久久久| 亚洲欧美日韩另类电影网站| 亚洲av第一区精品v没综合| 50天的宝宝边吃奶边哭怎么回事| 日本 av在线| 亚洲人成电影免费在线| 一二三四在线观看免费中文在| 亚洲熟妇熟女久久| 国内久久婷婷六月综合欲色啪| 精品电影一区二区在线| 亚洲五月婷婷丁香| 9热在线视频观看99| 久久香蕉国产精品| 黑人欧美特级aaaaaa片| 精品一品国产午夜福利视频| 桃色一区二区三区在线观看| 国产高清视频在线播放一区| 高清毛片免费观看视频网站| 欧美精品啪啪一区二区三区| 成在线人永久免费视频| 亚洲欧美日韩高清在线视频| 在线永久观看黄色视频| 三级毛片av免费| 一二三四在线观看免费中文在| 久久久国产精品麻豆| 亚洲男人的天堂狠狠| 午夜福利高清视频| 久久久久九九精品影院| 欧美不卡视频在线免费观看 | 在线免费观看的www视频| 十分钟在线观看高清视频www| 国产亚洲av高清不卡| 俄罗斯特黄特色一大片| 国产精品一区二区免费欧美| 久热这里只有精品99| 18禁黄网站禁片午夜丰满| 午夜激情av网站| 女性生殖器流出的白浆| 黄片播放在线免费| 日韩欧美一区视频在线观看| 亚洲成a人片在线一区二区| 日韩中文字幕欧美一区二区| 国产高清有码在线观看视频 | 麻豆国产av国片精品| 成人国产综合亚洲| 日本免费a在线| 啦啦啦观看免费观看视频高清 | 精品国内亚洲2022精品成人| 国产精品,欧美在线| 国产区一区二久久| 一进一出抽搐gif免费好疼| 免费人成视频x8x8入口观看| 久久久久久久久中文| 丰满的人妻完整版| 亚洲欧美精品综合久久99| 18禁观看日本| 在线观看午夜福利视频| 超碰成人久久| 91成年电影在线观看| 日韩欧美国产在线观看| 一区福利在线观看| 精品久久蜜臀av无| 午夜免费成人在线视频| 黄色片一级片一级黄色片| 精品少妇一区二区三区视频日本电影| 69av精品久久久久久| 亚洲av片天天在线观看| 国产精品一区二区三区四区久久 | 人人澡人人妻人| 国产真人三级小视频在线观看| 99精品久久久久人妻精品| 黄色成人免费大全| 色播亚洲综合网| 中文字幕高清在线视频| 长腿黑丝高跟| 久久国产精品人妻蜜桃| 午夜福利一区二区在线看| 国产精品亚洲av一区麻豆| 老熟妇仑乱视频hdxx| 亚洲精品久久国产高清桃花| 九色亚洲精品在线播放| 亚洲七黄色美女视频| 丝袜美足系列| 久久精品人人爽人人爽视色| 中文字幕av电影在线播放| 1024香蕉在线观看| 国产高清videossex| 天天躁狠狠躁夜夜躁狠狠躁| 极品人妻少妇av视频| 动漫黄色视频在线观看| 欧美成人一区二区免费高清观看 | 一本久久中文字幕| 精品少妇一区二区三区视频日本电影| 电影成人av| 欧美乱码精品一区二区三区| 嫩草影院精品99| 欧美最黄视频在线播放免费| 中文字幕人妻丝袜一区二区| 久久午夜综合久久蜜桃| 国产一卡二卡三卡精品| 亚洲天堂国产精品一区在线| 一区福利在线观看| 亚洲男人的天堂狠狠| 欧美av亚洲av综合av国产av| 两人在一起打扑克的视频| 午夜精品国产一区二区电影| av福利片在线| 午夜精品国产一区二区电影| 欧美日韩福利视频一区二区| 午夜两性在线视频| 1024视频免费在线观看| 一二三四在线观看免费中文在| 18禁美女被吸乳视频| 黑人操中国人逼视频| 黄色视频不卡| 脱女人内裤的视频| 99精品欧美一区二区三区四区| 欧美国产日韩亚洲一区| av视频免费观看在线观看| 成人亚洲精品av一区二区| 91成人精品电影| 精品欧美国产一区二区三| 亚洲午夜精品一区,二区,三区| 91av网站免费观看| 91精品国产国语对白视频| 一二三四社区在线视频社区8| 高清在线国产一区| 国产精品久久久久久精品电影 | 国产精品乱码一区二三区的特点 | 亚洲av成人一区二区三| 91麻豆精品激情在线观看国产| 在线观看日韩欧美| 777久久人妻少妇嫩草av网站| 可以在线观看毛片的网站| 午夜福利18| av网站免费在线观看视频| 精品一区二区三区av网在线观看| 国产成人系列免费观看| 国产熟女午夜一区二区三区| 亚洲一区中文字幕在线| 在线天堂中文资源库| 夜夜夜夜夜久久久久| 老熟妇乱子伦视频在线观看| 国产三级黄色录像| 久久伊人香网站| 欧美亚洲日本最大视频资源| 天天添夜夜摸| 级片在线观看| 精品久久久久久,| 中亚洲国语对白在线视频| 99精品在免费线老司机午夜| 国产不卡一卡二| 久久久久九九精品影院| 欧美绝顶高潮抽搐喷水| 99在线人妻在线中文字幕| 成年女人毛片免费观看观看9| 国产精品亚洲一级av第二区| 少妇裸体淫交视频免费看高清 | 99精品久久久久人妻精品| 欧美日韩黄片免| 青草久久国产| 亚洲精品国产区一区二| 99香蕉大伊视频| 亚洲精品在线美女| 午夜免费激情av| 国产亚洲av高清不卡| 视频在线观看一区二区三区| 夜夜看夜夜爽夜夜摸| 国产成人影院久久av| 在线观看66精品国产| 最近最新免费中文字幕在线| 一级a爱视频在线免费观看| 国产精品久久久久久人妻精品电影| 757午夜福利合集在线观看| 国产亚洲精品久久久久5区| 成人永久免费在线观看视频| 黄频高清免费视频| 成人国产综合亚洲| 18美女黄网站色大片免费观看| 国产精品98久久久久久宅男小说| 又大又爽又粗| 激情视频va一区二区三区| 高清毛片免费观看视频网站| 亚洲一区二区三区色噜噜| 精品欧美国产一区二区三| 色哟哟哟哟哟哟| 精品欧美国产一区二区三| 身体一侧抽搐| 真人一进一出gif抽搐免费| 久久婷婷成人综合色麻豆| 好男人在线观看高清免费视频 | 国产精品二区激情视频| 日日爽夜夜爽网站| 免费高清在线观看日韩| 韩国av一区二区三区四区| 国产精品精品国产色婷婷| netflix在线观看网站| 亚洲性夜色夜夜综合| 美国免费a级毛片| 免费观看人在逋| 两性午夜刺激爽爽歪歪视频在线观看 | 一区二区三区激情视频| 99在线视频只有这里精品首页| 12—13女人毛片做爰片一| 一级作爱视频免费观看| 欧美+亚洲+日韩+国产| 无人区码免费观看不卡| 91成人精品电影| 久久草成人影院| 亚洲第一av免费看| 国产精品永久免费网站| 国产成年人精品一区二区| 看免费av毛片| 亚洲av日韩精品久久久久久密| 亚洲午夜理论影院| 欧美成人性av电影在线观看| 久久国产精品影院| 国产色视频综合| 每晚都被弄得嗷嗷叫到高潮| 日韩有码中文字幕| 嫩草影视91久久| 夜夜夜夜夜久久久久| 97人妻精品一区二区三区麻豆 | 国产精品一区二区免费欧美| 亚洲精品国产精品久久久不卡| 男男h啪啪无遮挡| 午夜福利一区二区在线看| 在线播放国产精品三级| 免费在线观看完整版高清| 国产成人av激情在线播放| 在线天堂中文资源库| 丝袜美足系列| 狠狠狠狠99中文字幕| 亚洲欧美精品综合一区二区三区| 午夜精品在线福利| 亚洲av成人不卡在线观看播放网| 精品欧美一区二区三区在线| 久久天堂一区二区三区四区| 欧美最黄视频在线播放免费| 俄罗斯特黄特色一大片| 色婷婷久久久亚洲欧美| www.www免费av| 丝袜在线中文字幕| 操美女的视频在线观看| 国产又爽黄色视频| 高清在线国产一区| 国产黄a三级三级三级人| 9191精品国产免费久久| 久久久久久久精品吃奶| 国产亚洲精品久久久久久毛片| 亚洲成av人片免费观看| 天天躁夜夜躁狠狠躁躁| 首页视频小说图片口味搜索| 午夜福利一区二区在线看| 黄片播放在线免费| 国产一区二区三区在线臀色熟女| 99国产精品一区二区三区| 日韩av在线大香蕉| 黄网站色视频无遮挡免费观看| 99香蕉大伊视频| 女人高潮潮喷娇喘18禁视频| 欧洲精品卡2卡3卡4卡5卡区| 亚洲熟妇中文字幕五十中出| 亚洲中文日韩欧美视频| 久久久久久久精品吃奶| 国产99白浆流出| 美女 人体艺术 gogo| 成人三级黄色视频| 久久午夜综合久久蜜桃| 亚洲av日韩精品久久久久久密| 中文字幕精品免费在线观看视频| 免费看十八禁软件| 每晚都被弄得嗷嗷叫到高潮| 在线av久久热| 日韩欧美一区二区三区在线观看| 高潮久久久久久久久久久不卡| 欧美日韩中文字幕国产精品一区二区三区 | 欧美+亚洲+日韩+国产| 亚洲成国产人片在线观看| 男女做爰动态图高潮gif福利片 | 两个人免费观看高清视频| 欧美一级a爱片免费观看看 | 天堂√8在线中文| 国产熟女xx| 大型av网站在线播放| 欧美成人午夜精品| 久久久精品国产亚洲av高清涩受| 精品国产一区二区三区四区第35| 九色国产91popny在线| 女生性感内裤真人,穿戴方法视频| 久久精品亚洲精品国产色婷小说| 久久久久国内视频| 日韩精品青青久久久久久| 欧美成人免费av一区二区三区| 在线永久观看黄色视频| 精品久久久久久久人妻蜜臀av | 热re99久久国产66热| av视频在线观看入口| 亚洲三区欧美一区| 三级毛片av免费| 国产不卡一卡二| 淫妇啪啪啪对白视频| 黄色视频,在线免费观看| 国产精品99久久99久久久不卡| 99在线视频只有这里精品首页| 色在线成人网| av天堂久久9| 成人免费观看视频高清| 激情在线观看视频在线高清| 麻豆国产av国片精品| 精品久久久久久,| 性少妇av在线| 色综合亚洲欧美另类图片| 久久 成人 亚洲| 满18在线观看网站| 亚洲国产欧美日韩在线播放| 性欧美人与动物交配| 国产私拍福利视频在线观看| 9191精品国产免费久久| 视频在线观看一区二区三区| 伊人久久大香线蕉亚洲五| 亚洲第一欧美日韩一区二区三区| 深夜精品福利| 又黄又粗又硬又大视频| 一区二区三区精品91| 好男人电影高清在线观看| 两人在一起打扑克的视频| 精品国产美女av久久久久小说| 波多野结衣巨乳人妻| 国产成人精品无人区| 麻豆久久精品国产亚洲av| 女人被狂操c到高潮| 日韩欧美一区二区三区在线观看| 热re99久久国产66热| a级毛片在线看网站| 伊人久久大香线蕉亚洲五| 国产激情久久老熟女| 国产一区二区三区综合在线观看| 色婷婷久久久亚洲欧美| 别揉我奶头~嗯~啊~动态视频| 国产又爽黄色视频| 99精品在免费线老司机午夜| 欧美成人性av电影在线观看| 少妇粗大呻吟视频| 国产精品久久电影中文字幕| 国产精品日韩av在线免费观看 | 国产熟女午夜一区二区三区| 久久国产精品男人的天堂亚洲| 精品国产乱子伦一区二区三区| 精品国产乱码久久久久久男人| 亚洲午夜精品一区,二区,三区| av福利片在线| 男女做爰动态图高潮gif福利片 | 国产成年人精品一区二区| 搞女人的毛片| 久9热在线精品视频| 国产不卡一卡二| 国产成人免费无遮挡视频| 免费人成视频x8x8入口观看| 成人国产一区最新在线观看| 国产av一区在线观看免费| 天天添夜夜摸| 亚洲电影在线观看av| 在线av久久热| 一二三四社区在线视频社区8| www日本在线高清视频| 一级作爱视频免费观看| 亚洲片人在线观看| 色播在线永久视频| 久久久久九九精品影院| 亚洲国产毛片av蜜桃av| 午夜免费鲁丝| 69av精品久久久久久| 自线自在国产av| 欧美国产日韩亚洲一区| 国产精品精品国产色婷婷| 亚洲精品在线美女| 成熟少妇高潮喷水视频| 久久久久久国产a免费观看| 在线观看66精品国产| 两人在一起打扑克的视频| 午夜福利,免费看| 日韩国内少妇激情av| 日韩大码丰满熟妇| 12—13女人毛片做爰片一| 亚洲成a人片在线一区二区| 精品国产一区二区三区四区第35| 亚洲激情在线av| 高清毛片免费观看视频网站| 大型黄色视频在线免费观看| 女警被强在线播放| 国产一区在线观看成人免费| 国产高清videossex| 美女高潮到喷水免费观看| 超碰成人久久| 国产精品久久视频播放| 老司机靠b影院| 三级毛片av免费| 巨乳人妻的诱惑在线观看| 国产单亲对白刺激| 国产亚洲精品av在线| а√天堂www在线а√下载| 成人特级黄色片久久久久久久| 99久久精品国产亚洲精品| 精品国产美女av久久久久小说| 看黄色毛片网站| 亚洲国产精品成人综合色| 又紧又爽又黄一区二区| 在线观看免费午夜福利视频| 国产精品久久久久久亚洲av鲁大| 夜夜爽天天搞| 欧美国产精品va在线观看不卡| 国内毛片毛片毛片毛片毛片| 真人做人爱边吃奶动态| 美女高潮喷水抽搐中文字幕| 丰满人妻熟妇乱又伦精品不卡|