• <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ù)中尺寸公差的研究
    亚洲专区字幕在线| 亚洲色图av天堂| 亚洲精品在线美女| 国产伦理片在线播放av一区| 日本a在线网址| 日韩一卡2卡3卡4卡2021年| 国产亚洲午夜精品一区二区久久| 91精品国产国语对白视频| 超碰成人久久| 国产xxxxx性猛交| 极品教师在线免费播放| 最黄视频免费看| 黄色丝袜av网址大全| 久久精品亚洲精品国产色婷小说| 免费黄频网站在线观看国产| 成在线人永久免费视频| 日本vs欧美在线观看视频| 国产精品电影一区二区三区 | 亚洲视频免费观看视频| 精品久久久久久久毛片微露脸| 日韩大码丰满熟妇| 色尼玛亚洲综合影院| 怎么达到女性高潮| 国产成+人综合+亚洲专区| 一级a爱视频在线免费观看| 一级a爱视频在线免费观看| 99re在线观看精品视频| 黑丝袜美女国产一区| 男女下面插进去视频免费观看| 精品国产乱码久久久久久小说| 人人澡人人妻人| 九色亚洲精品在线播放| 成年动漫av网址| 性色av乱码一区二区三区2| 久久久精品国产亚洲av高清涩受| 精品少妇久久久久久888优播| 一区二区三区激情视频| 国产一区二区三区综合在线观看| 国产精品久久久久久人妻精品电影 | 一边摸一边抽搐一进一小说 | 777久久人妻少妇嫩草av网站| 免费少妇av软件| 狠狠婷婷综合久久久久久88av| 久久这里只有精品19| 精品国内亚洲2022精品成人 | 国产野战对白在线观看| 亚洲精品中文字幕一二三四区 | 在线永久观看黄色视频| 母亲3免费完整高清在线观看| 精品午夜福利视频在线观看一区 | 亚洲精品粉嫩美女一区| 如日韩欧美国产精品一区二区三区| av网站在线播放免费| 丁香六月欧美| 亚洲精品在线观看二区| 久久久久久久精品吃奶| 夫妻午夜视频| 一级黄色大片毛片| 亚洲成av片中文字幕在线观看| 国产成人影院久久av| av一本久久久久| 变态另类成人亚洲欧美熟女 | 啦啦啦中文免费视频观看日本| 国产精品久久久久久精品古装| 高清视频免费观看一区二区| 精品国产乱码久久久久久小说| 国产1区2区3区精品| 超碰97精品在线观看| 欧美黑人欧美精品刺激| 免费在线观看黄色视频的| 午夜福利视频精品| 午夜福利视频在线观看免费| 日韩视频一区二区在线观看| 成人国语在线视频| 国产精品99久久99久久久不卡| 久久久精品免费免费高清| 人人妻人人澡人人爽人人夜夜| 午夜老司机福利片| 在线观看舔阴道视频| 国产精品影院久久| 亚洲七黄色美女视频| 一区二区av电影网| 色尼玛亚洲综合影院| 国产激情久久老熟女| 精品国产亚洲在线| 一级a爱视频在线免费观看| 黄色毛片三级朝国网站| 高清欧美精品videossex| 久久久久久久国产电影| 大型黄色视频在线免费观看| 国产欧美日韩一区二区三| 日本黄色日本黄色录像| 精品少妇黑人巨大在线播放| 91成年电影在线观看| 一级,二级,三级黄色视频| 欧美成狂野欧美在线观看| 大陆偷拍与自拍| 成年人免费黄色播放视频| 亚洲国产看品久久| avwww免费| 另类亚洲欧美激情| 国产主播在线观看一区二区| 在线观看一区二区三区激情| 久久久国产欧美日韩av| 80岁老熟妇乱子伦牲交| 中文字幕制服av| 成人特级黄色片久久久久久久 | 亚洲欧美日韩另类电影网站| 老汉色∧v一级毛片| 岛国毛片在线播放| 国产精品秋霞免费鲁丝片| 亚洲第一青青草原| 一边摸一边抽搐一进一小说 | 一级黄色大片毛片| av有码第一页| 国产男女内射视频| 国产精品久久久久久精品古装| 19禁男女啪啪无遮挡网站| 老司机午夜十八禁免费视频| 国产在视频线精品| 一区二区日韩欧美中文字幕| 中文欧美无线码| 国精品久久久久久国模美| 婷婷成人精品国产| 国产精品av久久久久免费| 色综合欧美亚洲国产小说| 成年版毛片免费区| 99国产精品免费福利视频| 别揉我奶头~嗯~啊~动态视频| 热re99久久国产66热| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲一码二码三码区别大吗| 精品一区二区三区av网在线观看 | 老鸭窝网址在线观看| 久久午夜综合久久蜜桃| 纯流量卡能插随身wifi吗| 国产免费现黄频在线看| 成人18禁在线播放| 成人av一区二区三区在线看| 国产三级黄色录像| 色综合欧美亚洲国产小说| 女同久久另类99精品国产91| 精品国产一区二区三区久久久樱花| 女人精品久久久久毛片| 一区二区三区激情视频| 少妇裸体淫交视频免费看高清 | 亚洲av国产av综合av卡| 国产av国产精品国产| 国产精品影院久久| 久久精品成人免费网站| 王馨瑶露胸无遮挡在线观看| 色在线成人网| 成人黄色视频免费在线看| 国产成人免费无遮挡视频| 亚洲少妇的诱惑av| 18在线观看网站| 国产在线精品亚洲第一网站| 国产野战对白在线观看| 视频区欧美日本亚洲| 成人国产av品久久久| 久久久国产精品麻豆| 女人精品久久久久毛片| 777久久人妻少妇嫩草av网站| 99国产精品免费福利视频| aaaaa片日本免费| 欧美人与性动交α欧美软件| 一二三四社区在线视频社区8| 99热网站在线观看| 一边摸一边做爽爽视频免费| 性高湖久久久久久久久免费观看| 熟女少妇亚洲综合色aaa.| 欧美成人午夜精品| 精品免费久久久久久久清纯 | 亚洲欧美精品综合一区二区三区| 亚洲av日韩精品久久久久久密| 两人在一起打扑克的视频| 国产精品国产高清国产av | 如日韩欧美国产精品一区二区三区| 午夜91福利影院| 在线观看免费视频网站a站| 中国美女看黄片| 窝窝影院91人妻| 美女视频免费永久观看网站| 亚洲精品中文字幕一二三四区 | 亚洲九九香蕉| 9色porny在线观看| 国产亚洲精品第一综合不卡| 欧美亚洲 丝袜 人妻 在线| 丝袜美足系列| 他把我摸到了高潮在线观看 | 美女午夜性视频免费| 国产精品偷伦视频观看了| 精品高清国产在线一区| 人妻一区二区av| 国产麻豆69| 国产精品一区二区精品视频观看| 国产精品自产拍在线观看55亚洲 | 亚洲精品一卡2卡三卡4卡5卡| 亚洲,欧美精品.| 在线观看免费午夜福利视频| 一级a爱视频在线免费观看| 性高湖久久久久久久久免费观看| 久久国产精品人妻蜜桃| 2018国产大陆天天弄谢| 99国产精品99久久久久| 国产91精品成人一区二区三区 | 国产一区有黄有色的免费视频| 国产黄频视频在线观看| 不卡av一区二区三区| 欧美激情 高清一区二区三区| 操美女的视频在线观看| 天堂俺去俺来也www色官网| 色94色欧美一区二区| 亚洲欧洲日产国产| 韩国精品一区二区三区| 大香蕉久久成人网| 99九九在线精品视频| 久久久欧美国产精品| 三级毛片av免费| 侵犯人妻中文字幕一二三四区| 啦啦啦 在线观看视频| 一级片'在线观看视频| 黑人巨大精品欧美一区二区mp4| 亚洲精品美女久久av网站| av福利片在线| 成年版毛片免费区| 免费观看人在逋| 久久久久精品国产欧美久久久| 91字幕亚洲| 欧美精品人与动牲交sv欧美| 日日摸夜夜添夜夜添小说| 日本欧美视频一区| 亚洲天堂av无毛| 亚洲午夜精品一区,二区,三区| 精品人妻熟女毛片av久久网站| 人妻久久中文字幕网| 下体分泌物呈黄色| 成在线人永久免费视频| 久久 成人 亚洲| 日韩精品免费视频一区二区三区| 欧美日韩视频精品一区| 在线av久久热| xxxhd国产人妻xxx| 另类精品久久| 欧美日韩中文字幕国产精品一区二区三区 | 国产免费福利视频在线观看| 久久久久国产一级毛片高清牌| 国产精品1区2区在线观看. | 亚洲少妇的诱惑av| a级毛片在线看网站| 午夜免费鲁丝| 亚洲av成人不卡在线观看播放网| 香蕉丝袜av| 亚洲欧洲日产国产| 成人av一区二区三区在线看| 免费观看人在逋| 久久久久精品国产欧美久久久| 国产日韩欧美亚洲二区| 久久香蕉激情| 天堂俺去俺来也www色官网| 亚洲欧美一区二区三区黑人| 亚洲色图综合在线观看| 亚洲男人天堂网一区| 亚洲综合色网址| 亚洲av成人不卡在线观看播放网| 中文字幕人妻熟女乱码| 另类精品久久| 高清毛片免费观看视频网站 | 18禁裸乳无遮挡动漫免费视频| 亚洲国产欧美日韩在线播放| 十八禁人妻一区二区| 久久精品亚洲精品国产色婷小说| 国产精品九九99| 正在播放国产对白刺激| 精品亚洲乱码少妇综合久久| 欧美人与性动交α欧美精品济南到| 亚洲精品一卡2卡三卡4卡5卡| 国内毛片毛片毛片毛片毛片| 国产人伦9x9x在线观看| 成人国产一区最新在线观看| 亚洲熟女精品中文字幕| 国产日韩欧美视频二区| 亚洲一区二区三区欧美精品| 黄色视频在线播放观看不卡| 久久精品国产亚洲av高清一级| 中文字幕高清在线视频| 亚洲国产毛片av蜜桃av| 黄色视频,在线免费观看| 亚洲专区国产一区二区| 91麻豆精品激情在线观看国产 | 搡老熟女国产l中国老女人| 热99久久久久精品小说推荐| 美女主播在线视频| 欧美精品亚洲一区二区| videosex国产| 夫妻午夜视频| 97在线人人人人妻| 欧美亚洲 丝袜 人妻 在线| 老司机深夜福利视频在线观看| 日本黄色日本黄色录像| √禁漫天堂资源中文www| 国产深夜福利视频在线观看| 久久久久国内视频| 香蕉久久夜色| 亚洲精品一二三| 桃花免费在线播放| 欧美成人午夜精品| 久久人妻福利社区极品人妻图片| 欧美久久黑人一区二区| 午夜福利欧美成人| 国产欧美日韩综合在线一区二区| 侵犯人妻中文字幕一二三四区| 亚洲av日韩精品久久久久久密| 日韩大片免费观看网站| 一本综合久久免费| 久久亚洲精品不卡| 超色免费av| 亚洲国产欧美网| 国产精品香港三级国产av潘金莲| 男女下面插进去视频免费观看| 欧美另类亚洲清纯唯美| aaaaa片日本免费| 中文字幕精品免费在线观看视频| 免费少妇av软件| 又大又爽又粗| 欧美另类亚洲清纯唯美| 亚洲精品中文字幕在线视频| 人人妻人人澡人人爽人人夜夜| 夜夜夜夜夜久久久久| 不卡av一区二区三区| 亚洲精品久久成人aⅴ小说| 欧美亚洲 丝袜 人妻 在线| 9色porny在线观看| 久久国产精品影院| 国产精品九九99| 成人国语在线视频| 国产精品熟女久久久久浪| 精品乱码久久久久久99久播| 午夜老司机福利片| 女人爽到高潮嗷嗷叫在线视频| 精品国产一区二区三区四区第35| 亚洲精品自拍成人| 国产精品电影一区二区三区 | 欧美在线黄色| 国产黄色免费在线视频| 免费一级毛片在线播放高清视频 | 午夜久久久在线观看| 欧美av亚洲av综合av国产av| 18禁裸乳无遮挡动漫免费视频| 国产男靠女视频免费网站| 国产熟女午夜一区二区三区| 狠狠精品人妻久久久久久综合| 69精品国产乱码久久久| 亚洲午夜精品一区,二区,三区| 精品少妇黑人巨大在线播放| 午夜免费鲁丝| 十八禁网站免费在线| www.熟女人妻精品国产| 在线永久观看黄色视频| 精品久久久久久电影网| 精品一区二区三区视频在线观看免费 | 一级毛片电影观看| 国产一卡二卡三卡精品| 亚洲精品国产色婷婷电影| 免费看a级黄色片| 欧美精品一区二区大全| 亚洲欧洲精品一区二区精品久久久| 热re99久久精品国产66热6| 在线天堂中文资源库| 亚洲国产av新网站| 激情在线观看视频在线高清 | 91精品三级在线观看| 十分钟在线观看高清视频www| 一边摸一边抽搐一进一出视频| 国产一区二区激情短视频| 一本—道久久a久久精品蜜桃钙片| 在线观看66精品国产| 成人免费观看视频高清| 两人在一起打扑克的视频| 丝袜美足系列| 精品少妇黑人巨大在线播放| 别揉我奶头~嗯~啊~动态视频| 日日夜夜操网爽| 国产不卡av网站在线观看| 欧美黄色淫秽网站| 欧美精品亚洲一区二区| 精品免费久久久久久久清纯 | 午夜福利欧美成人| 丰满人妻熟妇乱又伦精品不卡| 午夜福利一区二区在线看| av网站在线播放免费| 成年动漫av网址| 欧美国产精品一级二级三级| 宅男免费午夜| xxxhd国产人妻xxx| 精品久久久久久久毛片微露脸| 午夜91福利影院| 久久人妻福利社区极品人妻图片| www日本在线高清视频| 在线看a的网站| 在线 av 中文字幕| 一级片免费观看大全| 精品少妇一区二区三区视频日本电影| 十八禁人妻一区二区| 热re99久久精品国产66热6| 麻豆国产av国片精品| 丝袜在线中文字幕| 午夜老司机福利片| 丁香欧美五月| 国产精品久久电影中文字幕 | videosex国产| 亚洲自偷自拍图片 自拍| 满18在线观看网站| 国产精品亚洲一级av第二区| 亚洲天堂av无毛| 91字幕亚洲| e午夜精品久久久久久久| 国产男靠女视频免费网站| 成年女人毛片免费观看观看9 | 精品人妻熟女毛片av久久网站| 日日爽夜夜爽网站| 久久人妻福利社区极品人妻图片| 手机成人av网站| 男女下面插进去视频免费观看| 97人妻天天添夜夜摸| 欧美在线一区亚洲| 一区二区三区激情视频| 国产欧美日韩综合在线一区二区| 丁香六月天网| √禁漫天堂资源中文www| 久久中文看片网| 脱女人内裤的视频| 丝袜美腿诱惑在线| 热99国产精品久久久久久7| 亚洲国产欧美一区二区综合| 亚洲五月婷婷丁香| 日本黄色视频三级网站网址 | 国产99久久九九免费精品| 成人亚洲精品一区在线观看| 久久精品亚洲熟妇少妇任你| 久久亚洲真实| 欧美成人午夜精品| 国产精品美女特级片免费视频播放器 | a在线观看视频网站| 动漫黄色视频在线观看| 国产欧美日韩综合在线一区二区| 一个人免费看片子| 高清在线国产一区| 精品少妇黑人巨大在线播放| 黑人操中国人逼视频| 国产免费福利视频在线观看| 麻豆乱淫一区二区| 99香蕉大伊视频| 老熟女久久久| 少妇的丰满在线观看| 久久精品国产亚洲av香蕉五月 | 少妇粗大呻吟视频| av片东京热男人的天堂| 国产又色又爽无遮挡免费看| 亚洲情色 制服丝袜| 视频区欧美日本亚洲| 老司机深夜福利视频在线观看| 亚洲精华国产精华精| 成人特级黄色片久久久久久久 | 动漫黄色视频在线观看| 精品一区二区三区四区五区乱码| 男女无遮挡免费网站观看| 国产日韩一区二区三区精品不卡| 色94色欧美一区二区| 色婷婷久久久亚洲欧美| 国产色视频综合| 日韩免费av在线播放| 亚洲成国产人片在线观看| av视频免费观看在线观看| 人人澡人人妻人| 人人妻人人澡人人爽人人夜夜| 免费日韩欧美在线观看| 国产成人影院久久av| 999久久久国产精品视频| 少妇的丰满在线观看| 一进一出好大好爽视频| 国产精品电影一区二区三区 | 欧美成人免费av一区二区三区 | 欧美老熟妇乱子伦牲交| 免费一级毛片在线播放高清视频 | 亚洲av第一区精品v没综合| av电影中文网址| 久久中文字幕人妻熟女| 丁香六月欧美| 国产一区二区激情短视频| 啦啦啦视频在线资源免费观看| 成在线人永久免费视频| 一级黄色大片毛片| 午夜老司机福利片| 久热这里只有精品99| 亚洲 欧美一区二区三区| 少妇裸体淫交视频免费看高清 | 精品国产乱子伦一区二区三区| 极品教师在线免费播放| 国产精品亚洲一级av第二区| 国产一区二区在线观看av| 久久精品aⅴ一区二区三区四区| 精品少妇久久久久久888优播| 国产一区二区三区在线臀色熟女 | av电影中文网址| 麻豆av在线久日| 国产高清视频在线播放一区| 成人三级做爰电影| 国产不卡av网站在线观看| 精品高清国产在线一区| 中亚洲国语对白在线视频| 妹子高潮喷水视频| 夜夜夜夜夜久久久久| 热99久久久久精品小说推荐| 大片免费播放器 马上看| 母亲3免费完整高清在线观看| 999精品在线视频| 久久久久视频综合| 99国产精品免费福利视频| 又黄又粗又硬又大视频| 十八禁人妻一区二区| 人人妻,人人澡人人爽秒播| 国产激情久久老熟女| 亚洲精品粉嫩美女一区| 国产精品久久久久久精品电影小说| 久久久国产成人免费| 十八禁高潮呻吟视频| 久久ye,这里只有精品| 美女高潮到喷水免费观看| 曰老女人黄片| bbb黄色大片| 久热这里只有精品99| 在线观看免费高清a一片| 亚洲欧美一区二区三区黑人| 免费在线观看影片大全网站| 久久久国产一区二区| 啦啦啦 在线观看视频| 色老头精品视频在线观看| 高清毛片免费观看视频网站 | 国产精品熟女久久久久浪| 制服诱惑二区| 精品国产一区二区三区四区第35| 日日爽夜夜爽网站| 亚洲国产欧美一区二区综合| 久久精品国产a三级三级三级| 久久久国产一区二区| 国产一区二区激情短视频| 黄色视频在线播放观看不卡| 精品国产一区二区久久| 亚洲精品一二三| 成人永久免费在线观看视频 | 精品亚洲成a人片在线观看| 人妻久久中文字幕网| 国产午夜精品久久久久久| 亚洲色图av天堂| 一区二区av电影网| 亚洲精品一二三| 又黄又粗又硬又大视频| 久9热在线精品视频| 国产又爽黄色视频| 日日摸夜夜添夜夜添小说| 女人高潮潮喷娇喘18禁视频| 国产高清videossex| 亚洲人成77777在线视频| 亚洲三区欧美一区| 中文字幕色久视频| 黄色视频不卡| 激情在线观看视频在线高清 | 亚洲欧美一区二区三区久久| 久久精品亚洲av国产电影网| 久久国产精品男人的天堂亚洲| 久久青草综合色| 国产成人精品久久二区二区91| 国产一区二区三区视频了| 大片电影免费在线观看免费| 91成年电影在线观看| 亚洲三区欧美一区| 国产单亲对白刺激| www.999成人在线观看| 欧美 亚洲 国产 日韩一| 久久九九热精品免费| 国产不卡一卡二| 国产成人精品久久二区二区免费| 亚洲人成电影观看| 色在线成人网| 国精品久久久久久国模美| 麻豆av在线久日| 波多野结衣av一区二区av| 亚洲 欧美一区二区三区| 久久天躁狠狠躁夜夜2o2o| 极品少妇高潮喷水抽搐| 久久久久久久久免费视频了| 人成视频在线观看免费观看| 亚洲专区中文字幕在线| 久久午夜综合久久蜜桃| 狠狠狠狠99中文字幕| 丝瓜视频免费看黄片| 黑人欧美特级aaaaaa片| 人成视频在线观看免费观看| 97在线人人人人妻| 一级毛片女人18水好多| av网站免费在线观看视频| 免费看a级黄色片| 搡老乐熟女国产| 一本一本久久a久久精品综合妖精| 国产成人精品久久二区二区免费| 国产精品.久久久| 香蕉国产在线看| 欧美+亚洲+日韩+国产| 日本vs欧美在线观看视频| 亚洲美女黄片视频| 免费av中文字幕在线| 91麻豆精品激情在线观看国产 | 日韩三级视频一区二区三区| 脱女人内裤的视频|