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

    使用P波快速測定國家臺網(wǎng)大震標準震級1

    2015-10-14 05:18:07翟璐媛
    震災防御技術(shù) 2015年2期
    關(guān)鍵詞:平均偏差速報寬頻

    翟璐媛 楊 陳

    ?

    使用P波快速測定國家臺網(wǎng)大震標準震級1

    翟璐媛 楊 陳

    (中國地震臺網(wǎng)中心,北京100045)

    本文針對國家臺網(wǎng)速報面波震級測定時間偏長和中深源地震震級速報有一定偏差的問題,采用IASPEI推薦的寬頻帶體波震級B及寬頻帶P波矩震級WP對2009—2013年國家臺網(wǎng)地震速報的大震進行了對比分析。對于經(jīng)過轉(zhuǎn)換成W后的B和WP震級來說,其結(jié)果均與我國速報地震發(fā)布的震級有一定的偏差,一般表現(xiàn)為偏小。其中,對于6.0—6.9級地震,B偏差相對較小,但離散度相對較大(整體偏差要比平均偏差大不少);對于7.0—7.9級地震,WP偏差相對較小;而對于8.0級以上地震,由于震級飽和等原因,B偏差較大,但WP偏差相對較小,一般主要表現(xiàn)為偏小??傮w來說,W(WP)的穩(wěn)定性要比W(B)更好一些(線性回歸的相關(guān)系數(shù)更大,標準誤差更?。?。對于綜合B和WP震級來說,由于采取分段平均的方法,結(jié)果的穩(wěn)定性有了一定的提高,但較大地震仍以偏小為主,如果在綜合震級P上加0.2,則可以得出與震級較為接近的結(jié)果。通過W(B)、W(WP)、P()、與W(GCMT)的對比,可以驗證綜合標準震級P()和國家臺網(wǎng)速報震級具備一定的可信度,而P()可作為P波快速測定的震級,所以用P()作為大震速報初報震級,在某種程度上是可行的。

    P波震級 體波震級B寬頻帶P波矩震級WP標準震級

    引言

    目前,根據(jù)“地震速報技術(shù)管理規(guī)定”(中國地震局監(jiān)測預報司,2013),國家臺網(wǎng)地震速報使用的震級是中華人民共和國國家標準即《國家標準震級(GB 17740-1999)》(中華人民共和國國家標準,1999)。對于淺源地震,通常以測定面波震級S作為國家標準震級;而對于中深源地震,通常以測定體波震級B表示,然后通過公式轉(zhuǎn)換成國家標準震級。這在實際運用過程中往往會存在以下兩方面的問題:第一種是體波震級B在轉(zhuǎn)換后與最終速報震級存在一定的偏差,這種偏差與震源深度及震級大小有關(guān),要么表現(xiàn)為中小地震的最終速報震級偏小,要么表現(xiàn)為大地震的最終速報震級偏大;第二種是面波震級S等待時間較長,尤其是對于發(fā)生在遠距離的地震,其等待時間更長。

    造成第一種現(xiàn)象的主要原因是,現(xiàn)有的地震速報系統(tǒng)所使用的體波震級B的計算公式為模擬公式,并未針對數(shù)字地震記錄進行更新,僅在數(shù)學上做了計算速度記錄的處理,其使用的值資料并不符合目前所使用的寬頻帶數(shù)字地震儀的需求。造成第二種現(xiàn)象的原因與面波震級S的性質(zhì)有關(guān),盡管使用面波測量震級較為穩(wěn)定,但由于其傳播速度的限制,已經(jīng)滿足不了對時間要求越來越高的地震速報需求,尤其是對于遠震表現(xiàn)的更為突出。

    針對上述兩種現(xiàn)象,本文試圖利用IASPEI(International Association of Seismology and Physics of the Earth's Interior,國際地震學與地球內(nèi)部物理學協(xié)會)推薦的寬頻帶體波震級B及寬頻帶P波矩震級WP,對滿足國家臺網(wǎng)地震速報的大震案例進行分析,并通過一定的方式進行轉(zhuǎn)換,以求得到滿足地震速報的P波快速測定震級的方法。

    1 原理及方法

    1.1 寬頻帶體波震級B測定方法

    IASPEI新標準推薦的寬頻帶體波震級B的計算公式為(Bormann等,2008):

    B=lg()max+(,) 0.2s<<30s,21°<<100° (1)

    式中,為體波質(zhì)點運動最大速度所對應的地動位移的最大振幅,單位為μm;為相應的周期,單位為s;(,)為震中距和震源深度的函數(shù)。

    當?shù)孛孢\動速度=2π/時,一般在實際計算中使用:

    B=lg(max/2π)+(,) (2)

    在減去可能的偏移值后,直接用與速度成正比的地震圖計算出最大速度的振幅max。根據(jù)IASPEI(2005)的規(guī)定:測算寬頻帶體波震級B的時間窗必須包含P波波列的最大振幅,包括P波、pP波、sP波,另外還有PcP波以及它們的尾波,但最好在PP波前結(jié)束。

    1.2 P波矩震級WP的測定方法

    寬頻帶P波矩震級是從遠場P波或pP波位移的垂直分量計算出的(Tsuboi,2000)。而矩震級WP的引入是基于這樣的假設(shè),即可以從寬頻帶位移波形u的P波段獲取地震矩0,其計算公式為:

    式中,和為P波沿傳播路徑的密度和波速;為震中距;P為輻射圖形。

    對于通常采用不加儀器響應校正的STS-1型寬頻帶地震儀獲取的地震圖,在積分速度之前需要十分小心地消除偏差。為了包含較遲的波至(如pP波),還需讀取合成位移地震圖中的第二峰,然后取第一峰和第二峰的差值,并將其與第一峰比較,取兩者中較大的來確定矩震級WP,這樣式(3)應寫為:

    其中,P1和P2分別為積分位移地震圖中P波的第一峰和第二峰值。

    由于每個地震臺站計算的矩震級WP都沒有加P校正,所以采用的標準矩震級計算公式為:

    (5)

    其中,0的單位為N·m,然后將W平均后加0.2,即得到WP。這是因為。為此需要對平均的0乘以2以補償,這就等同于對W加0.2,即WP=W+0.2。

    2 數(shù)據(jù)選取及預處理

    由于本文是針對國家臺網(wǎng)大震速報進行的P波快速測定震級的研究,所以筆者選取了2009年1月1日至2013年12月31日期間,國家地震臺網(wǎng)中心速報的全球6.0級以上地震共計698條作為基礎(chǔ)數(shù)據(jù)(全國地震速報目錄,2014),其震中分布圖1所示。

    為保證對比結(jié)果的一致,對于這698條地震目錄,筆者使用國家地震臺站的波形數(shù)據(jù),采用地震速報時所用的震中位置(只計算震級,不進行定位),并利用第一節(jié)中的方法計算了對應的B及WP,然后將其轉(zhuǎn)換成各自對應的W震級,最后使用轉(zhuǎn)換后的W震級與地震速報時采用的震級進行對比。同時為了適合國家臺網(wǎng)震級的快速測定,筆者對B的計算并沒有限制震中距;而對于WP的計算,也并未對儀器種類進行限制,其原因是目前中國國家臺網(wǎng)所使用的地震計以CTS-1系列和BBVS-120為主,并沒有使用STS-1型地震計。

    筆者使用Saul等(2007)和Whitmore等(2002)給出的W震級轉(zhuǎn)換公式:

    W(B)=1.33B-2.34 (6)

    W(WP)=1.31WP-1.91 (7)

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

    3.1 寬頻帶體波震級B計算結(jié)果分析

    在完成測定體波震級B后,可根據(jù)Saul等(2007)給出的計算公式再轉(zhuǎn)換成W,這樣即可進行W與的對比,圖2為具體的對比結(jié)果。

    從圖2可以看出,由B轉(zhuǎn)換成的W震級與相比,平均偏小了0.18,整體偏差(絕對值)為0.25,其計算結(jié)果的離散度較大。當在6.0—6.5之間時,兩者的偏差較為接近;而當>7.5級尤其是8級以上時,W(B)偏差就小得多,但由于B沒有8級以上的,所以這與B的飽和有一定的關(guān)系。表1為W(B)與分區(qū)間震級差。

    表1 不同震級區(qū)間的震級差

    從表1可以看出,隨著震級的增加,兩者的偏差隨之變大,平均偏差和整體偏差均是一樣;而對于8.0級以上地震,由于樣本較少和B震級飽和的原因,偏差較大達到1.06。

    利用本文的數(shù)據(jù)并結(jié)合線性回歸的方法,筆者得到的兩者之間的轉(zhuǎn)換公式為:

    =0.89W(B)+0.87 (8)

    其中,相關(guān)系數(shù)為0.8337;截距和斜率的標準誤差分別為0.1407和0.0224。

    由于以上僅使用了國家臺網(wǎng)的數(shù)據(jù)測定B,為了檢測B震級的方向性,筆者以西安(108.94°E,34.26°N)為中心,分四個象限對W(B)與進行了對比,結(jié)果如圖3所示。

    從圖3可以看出,對于測定的698條地震,W(B)與四個象限的對比并無明顯區(qū)別,與總體的對比結(jié)果基本一致。

    3.2 P波矩震級WP計算結(jié)果分析

    在完成測定寬頻帶P波矩震級WP后,可根據(jù)Whitmore等(2002)給出的計算公式轉(zhuǎn)換成W,再利用W與進行對比,圖4為具體的對比結(jié)果。

    從圖4可以看出,由WP轉(zhuǎn)換成的W震級與震級相比,平均偏小0.25,整體偏差(絕對值)為0.28,其計算結(jié)果以W(WP)偏小為主,尤其當為6—7級時,W(B)偏小較多,但當>7.5級時,與W(B)相比偏小的程度要少一些。表2為W(WP)與分區(qū)間震級偏差。

    表2 不同震級區(qū)間的震級差

    從表2可以看出,在6.0—6.9級區(qū)間的震級偏差較大,平均偏差和整體偏差均是一樣的;7.0—7.9級區(qū)間的震級偏差相對較小,而對于8.0級以上地震,平均偏差和整體偏差均為0.43,這與WP僅使用P波階段測量有一定的關(guān)系,同時與樣本數(shù)較少也有關(guān)系。

    利用本文的數(shù)據(jù)并結(jié)合線性回歸的方法,筆者得到的兩者之間的轉(zhuǎn)換公式為:

    =0.81W(WP)+1.45 (9)

    其中,相關(guān)系數(shù)為0.9275;截距和斜率的標準誤差分別為0.0767和0.0123。

    3.3 結(jié)合B和WP的綜合分析

    對于同一地震,取W(B)和W(WP)平均值,再與進行對比,對比結(jié)果如圖5所示。

    從圖5可以看出,相對于W(B)和W(WP)單獨的對比,如果簡單地取兩者的平均,對于7.5級以下地震,其震級的偏差有一定程度的提高,但仍以偏小為主;而對于7.5級以上地震,由于B本身飽和等因素的影響,這部分震級的偏差要比單獨的WP大。針對這種情況,筆者對震級進行了分段處理,具體原則為:對于W(WP)震級小于7.5級的地震,取W(B)和W(WP)的平均值;而對于W(WP)震級大于等于7.5級的地震,直接取W(WP)的值。筆者將此平均值命名為P(P波震級),再與進行對比,其結(jié)果如圖6所示。

    從圖6可以看出,與單獨的W(B)和W(WP)以及它們簡單的平均值相比,其結(jié)果的離散程度要好一些。對于大地震的震級偏差也相應的較小,其平均偏差為0.21,整體偏差為0.23。如果要以此作為速報震級的快速產(chǎn)出來說,直接在P上加0.2,則可以得出較為可靠的結(jié)果,這對于較大震級的地震也不會有太大的偏差,即可用P()=P+0.2表示。

    3.4W(B)、W(WP)、P()、與Global CMT震級的對比分析

    對于全球6級以上的地震,Global CMT的目錄中已經(jīng)提供了W震級,它是通過擬合全球臺網(wǎng)波形數(shù)據(jù)得到的結(jié)果,在國際上具有一定的認知度和權(quán)威性,Bormann等(2009)使用中國臺網(wǎng)48個臺站的資料進行了多種震級的對比。為了進一步驗證使用P波快速測定的W(B)、W(WP)、綜合震級P()以及國家臺網(wǎng)大震速報震級的結(jié)果,筆者將這四種震級與W(GCMT)進行了對比。在全部698條國家臺網(wǎng)速報目錄中找到了與Global CMT震級對應的687條,其對比結(jié)果如圖7所示。

    從圖7可以看出,相對于同國家臺網(wǎng)速報震級的對比,W(B)和W(MP)與W(GCMT)的對比結(jié)果要更好一些,但總體偏差仍不小;而轉(zhuǎn)換后的P波綜合震級P()和國家臺網(wǎng)速報震級與W(GCMT)對比的結(jié)果偏差較小,但P()仍有個別點的偏差較大,如日本“3·11”大地震與P波階段震級測定偏小有一定關(guān)系。

    使用線性回歸可得出四種震級與W(GCMT)的轉(zhuǎn)換公式,如表3所示。

    表3 四種轉(zhuǎn)換公式

    4 討論與結(jié)論

    對于經(jīng)過轉(zhuǎn)換成W后的B和WP震級來說,其結(jié)果均與我國速報地震發(fā)布的震級有一定的偏差,總體表現(xiàn)為偏小。其中,對于6.0—6.9級地震,B偏差相對較小,但離散度相對較大,整體偏差要比平均偏差大不少;對于7.0—7.9級地震,WP偏差相對較?。欢鴮τ?.0級以上地震,由于震級飽和等原因,B偏差較大,WP偏差相對較小,但仍以偏小為主??傮w來說,W(WP)的穩(wěn)定性要比W(B)更好一些(線性回歸的相關(guān)系數(shù)更大,標準誤差更?。?。

    對于綜合B和WP震級來說,由于采取分段平均的方法,結(jié)果的穩(wěn)定性有了一定的提高,但對于較大地震仍以偏小為主,但如果在綜合震級P上加0.2,則可以得出與震級較為接近的結(jié)果。筆者設(shè)想如果對兩者合成的權(quán)重進行調(diào)整,有可能會得出更加合理的結(jié)果。

    通過W(B)、W(WP)、P()、與W(GCMT)的對比后發(fā)現(xiàn),轉(zhuǎn)換后的P波綜合震級P()和國家臺網(wǎng)速報震級與W(GCMT)的對比結(jié)果偏差較小,說明這兩種結(jié)果具有一定程度的可信度,而P()作為P波快速測定的震級,用其作為大震速報震級的初步結(jié)果,這在某種層度上來說是可行的。

    由于受統(tǒng)計樣本數(shù)量的限制,尤其是較大地震的樣本數(shù)較少,所以本文對較大地震的統(tǒng)計還存在一定的偶然性,如果在較大地震的樣本數(shù)足夠多的情況下再進行分析,則可能得出更加合理的結(jié)論。

    致謝:國家地震臺網(wǎng)中心為本文提供了數(shù)據(jù)資料,審稿專家對本文提出了有益的建議,在此一并表示感謝!

    全國地震速報目錄(OL),2014.http://10.5.160.37:8080/index.jsp(內(nèi)網(wǎng)地址).

    中華人民共和國國家標準,1999.地震震級的規(guī)定(GB 17740-1999).北京:中國標準出版社.

    中國地震局監(jiān)測預報司,2013.地震速報技術(shù)管理規(guī)定(2013修訂版)[2013]26.北京:中國地震局.

    Bormann P. and Saul J., 2008. The New IASPEI Standard Broadband MagnitudeB.Seismol. Res. Lett., 79 (5): 698—705.

    Bormann P., Liu R.F., Xu Z.G. et al., 2009. First Application of the New IASPEI Teleseismic Magnitude Standards to Data of the China National Seismographic Network. Bulletin of the Seismological Society of America, 99 (3): 1868—1891. doi: 10.1785/0120080010.

    Saul J. and Bormann P., 2007. Rapid estimation of earthquake size using the broadband P-wave magnitudeB. 2007AGUFM. S53A1035S.

    Tsuboi S., 2000. Application ofWPto Tsunami Earthquake. Geophysical Research Letters, 27 (19): 3105—3108. doi:10.1029/2000GL011735.

    Whitmore P.M., Tsuboi S., Hirshorn B. and Sokolowski T.J., 2002. Magnitude-dependent Correction forWP. Science of Tsunami Hazards, 20 (4): 187—192.

    Rapid Determination of Standard Magnitude in National Network by Using P Wave

    Zhai Luyuan and Yang Chen

    (China Earthquake Networks Center,Beijing 100045,China)

    Usually the determinate time of surface wave is long and the magnitude of deep source earthquake is instable when National network earthquake information releases. In order to solve this problem we used body wave magnitudeBwhich IASPEI recommends and broadband P wave moment magnitudeWPcontrast the standard magnitudewhich National network determinate during 2009—2013. For the converted magnitudeWbyBandWP, the results have a certain deviation to National network earthquake information releases standard magnitude, are characterized by small. For magnitude from 6.0 to 6.9,Bdeviation is relatively smaller, but the discrete degree is relatively larger (the overall deviation is greater than the average deviation); for magnitude from 7.0 to 7.9,WPdeviation is relatively smaller; for the magnitude more than 8.0, due to the reason such as magnitude saturation, the deviation ofBis relatively larger than that ofWP, but is still smaller than. In general, the stability ofW(WP) is better thanW(B) (the linear regression correlation coefficient is bigger, and the standard error is smaller). For the magnitude comprehensiveBandWP, due to the piecewise average method, the stability of the results is improved in certain degree, but is still smaller for large earthquakes. If add 0.2 to the comprehensive magnitudeP, we can get results close to the standard magnitude. By comparison withW(GCMT),we can verifyP() and(CENC) have a certain degree of credibility.P() as a rapid determination of P-wave, serve as a preliminary result of quick report of great earthquake is also possible.

    P wave magnitude; Body wave magnitude mB; Broad band P wave moment magnitudeWP; Standard magnitude

    中國地震局測震臺網(wǎng)青年骨干培養(yǎng)專項(20130217)

    2014-12-05

    翟璐媛,女,生于1982年。2006年畢業(yè)于中國科技大學,工程師。主要從事地震監(jiān)測及臺網(wǎng)管理等工作。 E-mail:zly@seis.ac.cn

    楊陳,男,生于1979年。高級工程師。主要從事地震監(jiān)測與地震預警方面的研究。E-mail:yangchen@seis.ac.cn

    猜你喜歡
    平均偏差速報寬頻
    河北地方性震級量規(guī)函數(shù)與方位角校正值研究1
    寬頻高磁導率R10k軟磁材料的開發(fā)
    山東冶金(2022年2期)2022-08-08 01:50:52
    FY-3C/VIRR西北太平洋區(qū)域海表溫度精度評估?
    每月速報
    空中之家(2018年1期)2018-01-31 01:50:30
    每月速報
    空中之家(2017年11期)2017-11-28 05:27:41
    每月速報
    空中之家(2016年1期)2016-05-17 04:47:41
    基于矢量匹配法的扼流變壓器的寬頻建模
    電氣化鐵道(2016年4期)2016-04-16 05:59:40
    寬頻鎖相的一種實現(xiàn)方法
    電測與儀表(2015年8期)2015-04-09 11:50:10
    脛前動脈穿刺可行性及心肺流轉(zhuǎn)下脛前動脈與橈動脈壓力監(jiān)測的一致性研究
    數(shù)控機床自動編程技術(shù)中尺寸公差的研究
    午夜精品久久久久久毛片777| 亚洲欧美日韩东京热| 国产精品乱码一区二三区的特点| 国产午夜精品论理片| 乱码一卡2卡4卡精品| 国产主播在线观看一区二区| .国产精品久久| 久久久久国产精品人妻aⅴ院| 国产免费男女视频| 亚洲欧美日韩东京热| 亚洲,欧美,日韩| 国产探花极品一区二区| 国产精品亚洲美女久久久| 在线国产一区二区在线| 啦啦啦观看免费观看视频高清| 精品午夜福利视频在线观看一区| 亚洲人成伊人成综合网2020| 在线看三级毛片| 国产亚洲欧美98| 国产69精品久久久久777片| 亚洲久久久久久中文字幕| 在线播放无遮挡| 免费在线观看影片大全网站| 97人妻精品一区二区三区麻豆| 久久久久久久午夜电影| 国产精品一区www在线观看 | 国产私拍福利视频在线观看| 在线观看午夜福利视频| 免费黄网站久久成人精品| 两性午夜刺激爽爽歪歪视频在线观看| 国产三级中文精品| 国产精品精品国产色婷婷| 久99久视频精品免费| 十八禁国产超污无遮挡网站| 999久久久精品免费观看国产| 日本在线视频免费播放| 亚洲图色成人| 亚洲色图av天堂| 99视频精品全部免费 在线| 午夜免费成人在线视频| 18禁在线播放成人免费| 亚洲欧美精品综合久久99| 国产成人一区二区在线| 国产 一区精品| 舔av片在线| 亚洲va在线va天堂va国产| 亚洲国产精品久久男人天堂| 国产精品国产高清国产av| 少妇人妻一区二区三区视频| 日韩av在线大香蕉| 成人亚洲精品av一区二区| 久久精品国产亚洲av涩爱 | 日本黄色视频三级网站网址| 亚洲熟妇中文字幕五十中出| 麻豆成人午夜福利视频| 亚洲国产精品久久男人天堂| 中文字幕免费在线视频6| 91麻豆精品激情在线观看国产| 国产精品无大码| 国产男靠女视频免费网站| 午夜福利在线在线| 日本黄大片高清| 成人国产一区最新在线观看| 亚洲国产精品sss在线观看| 亚洲综合色惰| 精品久久国产蜜桃| 一个人免费在线观看电影| 久久亚洲真实| 麻豆成人午夜福利视频| 校园春色视频在线观看| 亚洲av中文字字幕乱码综合| 亚洲成人久久性| 精品午夜福利在线看| 一本久久中文字幕| 国产探花在线观看一区二区| 国产高清视频在线播放一区| 俄罗斯特黄特色一大片| 欧美另类亚洲清纯唯美| 夜夜看夜夜爽夜夜摸| 99久久中文字幕三级久久日本| 亚洲五月天丁香| 日韩 亚洲 欧美在线| 亚洲欧美激情综合另类| 亚洲成人精品中文字幕电影| 日本五十路高清| 日本与韩国留学比较| 18禁黄网站禁片免费观看直播| 黄色女人牲交| 亚洲熟妇中文字幕五十中出| 性插视频无遮挡在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 成人高潮视频无遮挡免费网站| 国产午夜精品久久久久久一区二区三区 | 中文亚洲av片在线观看爽| 久久99热这里只有精品18| 一区二区三区免费毛片| 久久久久久久久久久丰满 | 亚洲av熟女| 欧美成人免费av一区二区三区| 人妻制服诱惑在线中文字幕| 一进一出好大好爽视频| 国产成人影院久久av| 中文在线观看免费www的网站| 亚洲国产高清在线一区二区三| 十八禁网站免费在线| 在线免费十八禁| 成人欧美大片| 无遮挡黄片免费观看| 嫩草影院新地址| 欧美日韩亚洲国产一区二区在线观看| 亚洲人成网站在线播放欧美日韩| 国产一区二区在线av高清观看| bbb黄色大片| 国产精品不卡视频一区二区| 一区福利在线观看| 麻豆国产av国片精品| 日韩大尺度精品在线看网址| 自拍偷自拍亚洲精品老妇| 亚洲av一区综合| 如何舔出高潮| 亚洲中文字幕日韩| 我要看日韩黄色一级片| 一级av片app| 欧美日韩瑟瑟在线播放| 一区二区三区高清视频在线| 男人和女人高潮做爰伦理| 俺也久久电影网| 成人鲁丝片一二三区免费| 啪啪无遮挡十八禁网站| 欧美性猛交╳xxx乱大交人| 99久久无色码亚洲精品果冻| 亚洲av免费高清在线观看| 一级av片app| 99精品在免费线老司机午夜| 欧美精品啪啪一区二区三区| ponron亚洲| 一进一出好大好爽视频| 我要搜黄色片| 亚洲四区av| 老司机福利观看| 国产精品野战在线观看| 黄色日韩在线| 国产成年人精品一区二区| 能在线免费观看的黄片| 18+在线观看网站| 亚洲乱码一区二区免费版| 亚洲七黄色美女视频| 欧洲精品卡2卡3卡4卡5卡区| 国产爱豆传媒在线观看| 69av精品久久久久久| 国产高清三级在线| 2021天堂中文幕一二区在线观| 超碰av人人做人人爽久久| 精品午夜福利视频在线观看一区| 欧美黑人巨大hd| 又黄又爽又刺激的免费视频.| 日韩人妻高清精品专区| 一级毛片久久久久久久久女| 国产成人a区在线观看| ponron亚洲| 一级毛片久久久久久久久女| 亚洲久久久久久中文字幕| 全区人妻精品视频| 在线观看美女被高潮喷水网站| 国产久久久一区二区三区| 夜夜看夜夜爽夜夜摸| 久久久久久国产a免费观看| 免费在线观看成人毛片| av.在线天堂| 久久久久久久午夜电影| 久久久久久久久久久丰满 | 少妇的逼好多水| 成年免费大片在线观看| 网址你懂的国产日韩在线| 国产私拍福利视频在线观看| videossex国产| 三级毛片av免费| 国产男人的电影天堂91| a级毛片a级免费在线| 亚洲精品影视一区二区三区av| 最近在线观看免费完整版| 性欧美人与动物交配| 欧美在线一区亚洲| 亚洲va在线va天堂va国产| 亚洲av一区综合| 91久久精品国产一区二区成人| 99久久中文字幕三级久久日本| 女生性感内裤真人,穿戴方法视频| 国产精品自产拍在线观看55亚洲| 精品欧美国产一区二区三| 日本撒尿小便嘘嘘汇集6| av在线老鸭窝| 欧美精品啪啪一区二区三区| 日本五十路高清| 色哟哟·www| 春色校园在线视频观看| 国产亚洲欧美98| 亚洲av中文字字幕乱码综合| 亚洲电影在线观看av| 久久国产乱子免费精品| 一个人免费在线观看电影| 中出人妻视频一区二区| 精品国内亚洲2022精品成人| 精品人妻偷拍中文字幕| 国产精品久久视频播放| 亚洲av中文字字幕乱码综合| 少妇丰满av| 国内精品宾馆在线| 亚洲欧美日韩东京热| 五月玫瑰六月丁香| 国产精品久久久久久久久免| 日韩欧美 国产精品| 欧美日韩乱码在线| 又黄又爽又刺激的免费视频.| 人人妻,人人澡人人爽秒播| 岛国在线免费视频观看| 国产主播在线观看一区二区| 能在线免费观看的黄片| 波野结衣二区三区在线| 我要看日韩黄色一级片| 联通29元200g的流量卡| 热99在线观看视频| 久久久国产成人免费| 一个人看视频在线观看www免费| 免费电影在线观看免费观看| 一级av片app| 日本-黄色视频高清免费观看| 亚洲最大成人av| 免费黄网站久久成人精品| 两个人的视频大全免费| 欧美在线一区亚洲| 色5月婷婷丁香| 午夜福利欧美成人| 少妇人妻精品综合一区二区 | 99精品在免费线老司机午夜| 成人国产麻豆网| av女优亚洲男人天堂| 欧美绝顶高潮抽搐喷水| 露出奶头的视频| 亚洲美女搞黄在线观看 | 国产成人福利小说| 久久精品国产鲁丝片午夜精品 | 老司机福利观看| 麻豆久久精品国产亚洲av| 亚洲av一区综合| 久久久精品大字幕| 久久精品国产亚洲av香蕉五月| 国产精品美女特级片免费视频播放器| 国产高清不卡午夜福利| 一区二区三区四区激情视频 | 久久精品国产鲁丝片午夜精品 | 久久精品国产鲁丝片午夜精品 | 国产av一区在线观看免费| 悠悠久久av| 日日干狠狠操夜夜爽| 国产色婷婷99| 精品国产三级普通话版| 国产亚洲欧美98| 国产一级毛片七仙女欲春2| 免费看美女性在线毛片视频| 久久国产精品人妻蜜桃| 男女下面进入的视频免费午夜| 国产美女午夜福利| 不卡一级毛片| 欧美在线一区亚洲| 国产免费男女视频| 国产一区二区在线av高清观看| 精品久久久久久成人av| 亚洲专区中文字幕在线| 美女cb高潮喷水在线观看| 一个人看视频在线观看www免费| 国产亚洲91精品色在线| 国产在视频线在精品| 日本与韩国留学比较| 日本一二三区视频观看| 亚洲av五月六月丁香网| 深爱激情五月婷婷| 麻豆一二三区av精品| 国产精品,欧美在线| 欧美3d第一页| 国产乱人视频| 亚洲av免费高清在线观看| 欧美日韩瑟瑟在线播放| 天堂√8在线中文| 亚洲欧美精品综合久久99| 成人av一区二区三区在线看| a在线观看视频网站| a级一级毛片免费在线观看| 我的女老师完整版在线观看| 一级av片app| 国产乱人视频| 99精品久久久久人妻精品| 国产精品1区2区在线观看.| 亚洲av日韩精品久久久久久密| 欧洲精品卡2卡3卡4卡5卡区| 精华霜和精华液先用哪个| 很黄的视频免费| 999久久久精品免费观看国产| 我要搜黄色片| 免费一级毛片在线播放高清视频| 美女高潮喷水抽搐中文字幕| 精品久久久久久成人av| 日韩中字成人| 午夜激情福利司机影院| 国产高清视频在线播放一区| 久久久国产成人免费| 又紧又爽又黄一区二区| 日本 欧美在线| 色5月婷婷丁香| 久久久久久久午夜电影| 免费观看精品视频网站| 国产av不卡久久| 精品久久久久久久久久久久久| 久久人人精品亚洲av| 99精品久久久久人妻精品| 99视频精品全部免费 在线| 日本a在线网址| 免费av观看视频| 哪里可以看免费的av片| 精品一区二区三区av网在线观看| 久久精品国产清高在天天线| 最近最新中文字幕大全电影3| 少妇被粗大猛烈的视频| 久久草成人影院| 人人妻,人人澡人人爽秒播| 精品人妻熟女av久视频| 久久精品91蜜桃| 91久久精品国产一区二区成人| 成人永久免费在线观看视频| 日韩强制内射视频| 日韩欧美在线乱码| 国产一区二区在线av高清观看| 国产色爽女视频免费观看| 国内揄拍国产精品人妻在线| 老司机午夜福利在线观看视频| 久久精品人妻少妇| 51国产日韩欧美| 国产午夜福利久久久久久| 久久午夜福利片| 国内精品一区二区在线观看| 国产伦人伦偷精品视频| 国产成人一区二区在线| 美女xxoo啪啪120秒动态图| 麻豆一二三区av精品| 国产欧美日韩精品亚洲av| 老女人水多毛片| 国产精品av视频在线免费观看| 欧美中文日本在线观看视频| 亚洲国产精品sss在线观看| aaaaa片日本免费| 国内少妇人妻偷人精品xxx网站| 成年人黄色毛片网站| 在线观看66精品国产| 一区二区三区四区激情视频 | 免费人成在线观看视频色| 国产蜜桃级精品一区二区三区| 色av中文字幕| 国产不卡一卡二| 夜夜夜夜夜久久久久| 欧美又色又爽又黄视频| 成人av在线播放网站| 永久网站在线| 日韩亚洲欧美综合| 欧美一区二区精品小视频在线| 给我免费播放毛片高清在线观看| 亚洲第一区二区三区不卡| 久久久久九九精品影院| 国产白丝娇喘喷水9色精品| 搡老妇女老女人老熟妇| 欧美日韩瑟瑟在线播放| 亚洲黑人精品在线| 麻豆成人av在线观看| 日本a在线网址| 一个人看视频在线观看www免费| 国内毛片毛片毛片毛片毛片| 成人永久免费在线观看视频| 两个人视频免费观看高清| 国产成人一区二区在线| 两个人视频免费观看高清| 日韩亚洲欧美综合| 永久网站在线| 能在线免费观看的黄片| 日韩精品青青久久久久久| 熟女人妻精品中文字幕| 999久久久精品免费观看国产| 国产精华一区二区三区| 99视频精品全部免费 在线| 桃红色精品国产亚洲av| 国产一级毛片七仙女欲春2| 免费观看的影片在线观看| 在线a可以看的网站| 国内毛片毛片毛片毛片毛片| 两性午夜刺激爽爽歪歪视频在线观看| 日韩欧美免费精品| 欧美绝顶高潮抽搐喷水| 日韩av在线大香蕉| 黄色日韩在线| 国模一区二区三区四区视频| 国产精品亚洲一级av第二区| 男人的好看免费观看在线视频| 在线天堂最新版资源| 亚洲五月天丁香| 国产乱人伦免费视频| 99国产精品一区二区蜜桃av| 精品福利观看| 久久久久九九精品影院| 三级国产精品欧美在线观看| 日韩 亚洲 欧美在线| 99国产极品粉嫩在线观看| 夜夜夜夜夜久久久久| 亚洲乱码一区二区免费版| 亚洲成a人片在线一区二区| 美女大奶头视频| 如何舔出高潮| 少妇丰满av| 别揉我奶头 嗯啊视频| 香蕉av资源在线| 国产aⅴ精品一区二区三区波| 久久久久久伊人网av| 一本精品99久久精品77| 美女高潮喷水抽搐中文字幕| 成人欧美大片| 极品教师在线免费播放| 看免费成人av毛片| 欧美在线一区亚洲| 欧美区成人在线视频| 亚洲欧美日韩高清专用| 男女啪啪激烈高潮av片| 免费看av在线观看网站| 婷婷六月久久综合丁香| 十八禁国产超污无遮挡网站| 久久人人爽人人爽人人片va| 在线国产一区二区在线| 日本爱情动作片www.在线观看 | 精品久久久久久久久亚洲 | 哪里可以看免费的av片| 亚洲午夜理论影院| 亚洲aⅴ乱码一区二区在线播放| 国产精品一区二区三区四区免费观看 | 欧美日本亚洲视频在线播放| 久久这里只有精品中国| 国产av一区在线观看免费| 国产精品免费一区二区三区在线| 啦啦啦观看免费观看视频高清| 久久久久久国产a免费观看| 小蜜桃在线观看免费完整版高清| 一级毛片久久久久久久久女| 国内精品一区二区在线观看| 99久久精品国产国产毛片| 国产精品一区二区性色av| 日韩 亚洲 欧美在线| 色在线成人网| 内地一区二区视频在线| 美女高潮喷水抽搐中文字幕| 赤兔流量卡办理| av在线天堂中文字幕| 老司机深夜福利视频在线观看| 欧美成人a在线观看| 永久网站在线| 免费在线观看成人毛片| 一a级毛片在线观看| 九色国产91popny在线| 亚洲av成人av| 在线播放无遮挡| 亚洲国产精品成人综合色| 国产高清不卡午夜福利| 熟女人妻精品中文字幕| 亚洲国产精品sss在线观看| 亚洲图色成人| 亚洲18禁久久av| 小蜜桃在线观看免费完整版高清| 一区二区三区四区激情视频 | 免费看光身美女| 69av精品久久久久久| 国产熟女欧美一区二区| 亚洲人与动物交配视频| 婷婷色综合大香蕉| 精品人妻一区二区三区麻豆 | 18禁在线播放成人免费| 少妇高潮的动态图| 久久久国产成人免费| 亚洲成人精品中文字幕电影| 尤物成人国产欧美一区二区三区| 九色成人免费人妻av| 亚洲国产日韩欧美精品在线观看| 毛片一级片免费看久久久久 | 99久久无色码亚洲精品果冻| 国产亚洲精品久久久久久毛片| 国产欧美日韩一区二区精品| 午夜福利在线在线| 欧美高清性xxxxhd video| 国产高清视频在线观看网站| 韩国av一区二区三区四区| 天天躁日日操中文字幕| 欧美一区二区精品小视频在线| 成年女人毛片免费观看观看9| 国产精品久久久久久精品电影| 啦啦啦韩国在线观看视频| 免费看美女性在线毛片视频| 久久香蕉精品热| 狂野欧美白嫩少妇大欣赏| 国产不卡一卡二| 日韩大尺度精品在线看网址| 一级毛片久久久久久久久女| 国产亚洲精品久久久久久毛片| 非洲黑人性xxxx精品又粗又长| 一a级毛片在线观看| 中文字幕久久专区| 嫩草影院新地址| 夜夜看夜夜爽夜夜摸| 99热6这里只有精品| 精品久久久久久久久亚洲 | 欧美一区二区亚洲| 欧美日韩综合久久久久久 | 麻豆久久精品国产亚洲av| 最新在线观看一区二区三区| 一进一出好大好爽视频| 亚洲av中文字字幕乱码综合| 国产精品久久久久久av不卡| 国产精品综合久久久久久久免费| 亚洲av不卡在线观看| 国产成年人精品一区二区| 国产伦一二天堂av在线观看| 一边摸一边抽搐一进一小说| 国产爱豆传媒在线观看| 校园春色视频在线观看| 国产av不卡久久| 免费看光身美女| 1024手机看黄色片| 国产真实伦视频高清在线观看 | 国产精品精品国产色婷婷| 国产精品98久久久久久宅男小说| 18禁黄网站禁片午夜丰满| 日本 av在线| 啦啦啦韩国在线观看视频| 成人午夜高清在线视频| 身体一侧抽搐| 蜜桃亚洲精品一区二区三区| 88av欧美| 亚洲va日本ⅴa欧美va伊人久久| 给我免费播放毛片高清在线观看| 精品久久久久久,| 国产精品一区二区免费欧美| 一级毛片久久久久久久久女| 日本黄大片高清| 欧美3d第一页| 国产男人的电影天堂91| 日本熟妇午夜| 午夜爱爱视频在线播放| 国内揄拍国产精品人妻在线| 看黄色毛片网站| 最近最新免费中文字幕在线| 精品人妻1区二区| 亚洲美女搞黄在线观看 | xxxwww97欧美| 日日撸夜夜添| 搡老熟女国产l中国老女人| 欧美成人一区二区免费高清观看| 亚洲精品亚洲一区二区| 成人亚洲精品av一区二区| 亚洲自偷自拍三级| 日本熟妇午夜| 简卡轻食公司| 黄色女人牲交| 91在线观看av| 在现免费观看毛片| 午夜福利欧美成人| 午夜视频国产福利| 97热精品久久久久久| 亚洲精品在线观看二区| 久久久久精品国产欧美久久久| 俄罗斯特黄特色一大片| 亚洲电影在线观看av| 婷婷丁香在线五月| 欧美一区二区国产精品久久精品| 欧美最新免费一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 国内少妇人妻偷人精品xxx网站| 国内久久婷婷六月综合欲色啪| 狂野欧美激情性xxxx在线观看| 女同久久另类99精品国产91| 亚洲一区高清亚洲精品| 久久中文看片网| 久久久久久久久久成人| 精品一区二区三区视频在线| 99久久精品热视频| 久久久久久久久久成人| 午夜精品久久久久久毛片777| 黄色女人牲交| 伦精品一区二区三区| 精品欧美国产一区二区三| 国产主播在线观看一区二区| 国产亚洲精品综合一区在线观看| 给我免费播放毛片高清在线观看| 少妇丰满av| 亚洲午夜理论影院| 男人和女人高潮做爰伦理| 日本 av在线| 在线观看舔阴道视频| 国产欧美日韩一区二区精品| eeuss影院久久| 日日摸夜夜添夜夜添小说| 国产精品久久久久久av不卡| 人人妻人人澡欧美一区二区| 精品一区二区三区av网在线观看| 啪啪无遮挡十八禁网站| 91午夜精品亚洲一区二区三区 | 国产av不卡久久| 亚洲精品一卡2卡三卡4卡5卡| 亚洲成a人片在线一区二区| 国产国拍精品亚洲av在线观看| 亚洲av二区三区四区| 午夜福利视频1000在线观看| 日韩一本色道免费dvd| 九色国产91popny在线|