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

    P波衰減規(guī)律應用于地震預警研究*

    2016-04-11 05:39:02ItzhakLiorAlonZivRaulMadariaga
    地震科學進展 2016年2期
    關鍵詞:震級方根震源

    Itzhak Lior, Alon Ziv, Raul Madariaga

    1) Department of Geosciences, Tel-Aviv University,P.O. Box 39040,

    Tel-Aviv 6997801,Israel

    2) Laboratoire de Géologie,Ecole Normale Supérieure, 24 rue Lhomond,

    75231 Paris cedex 05,F(xiàn)rance

    ?

    P波衰減規(guī)律應用于地震預警研究*

    Itzhak Lior1), Alon Ziv1), Raul Madariaga2)

    1) Department of Geosciences, Tel-Aviv University,P.O. Box 39040,

    Tel-Aviv 6997801,Israel

    2) Laboratoire de Géologie,Ecole Normale Supérieure, 24 rue Lhomond,

    75231 Paris cedex 05,F(xiàn)rance

    摘要一些廣泛應用和正在測試的地震預警算法是采用與地面運動峰值振幅和震源距離相關的震級經(jīng)驗公式。 這些方法的有效性在一定程度上依賴于可用于設定那些方程的擬合系數(shù)的離線數(shù)據(jù)集的數(shù)量和質量是否充分。 然而,獲取控制P波衰減主要因素的物理涵義對解決數(shù)據(jù)集有限問題具有指導意義。 本研究中,理論表達式來自P波位移drms和波速vrms的均方根與地震矩、 應力降和震源距離的關系。

    利用南加州和日本的地震數(shù)據(jù),將觀測衰減用于驗證理論衰減規(guī)律。 地面運動的觀測值和預測值顯示出較好的一致性。 在加州和日本相似的地面運動衰減表明,這兩個區(qū)域的衰減規(guī)律相似,推斷它們也可以用于其他沒有長期觀測數(shù)據(jù)校準的區(qū)域。

    由于drms比vrms更易受地震矩的影響,將衰減規(guī)律用于drms比用于vrms能得到更好的震級預測結果。 結果表明,drms/vrms比值與破裂長度特征成比例,應力降是地震矩和drms/vrms立方的函數(shù)。 此結果可以用來確定應力降,該方法與之前所用過的方法無關。 最后,結果表明地面運動的均方根與其峰值成比例。

    引言

    地震預警系統(tǒng)(earthquake early warning systems,EEWS)利用前幾秒的地震記錄通過實時算法預測地震動烈度或評估警報水平。 具有代表性的算法是利用離線地震記錄通過回歸分析獲得經(jīng)驗關系[1-8]。 P波衰減規(guī)律就是這類經(jīng)驗關系中的一種,它描述P波峰值速度或位移與震級和震源距離的關系。 由于場地和特有震源特性影響地面運動,一種構造背景下所得到的衰減規(guī)律并不適用于別的區(qū)域。 此外,除非數(shù)據(jù)集非常大,這些相關性需要在不同參數(shù)間進行取舍,導致擬合系數(shù)不穩(wěn)定。 為了解決數(shù)據(jù)集有限問題,且便于參數(shù)調整,在校準不同地質構造時,獲取控制P波衰減主要因素的物理涵義具有指導意義。

    本研究中,理論表達式來自與地震矩、 應力降和震源距離相關的P波地面運動(位移和速度)的均方根。 利用南加州和日本的地震數(shù)據(jù),針對觀測衰減驗證理論衰減規(guī)律。 本文還討論了將位移和波速衰減規(guī)律用于實時震級的確定。 最后,結果表明,地面運動均方根和峰值呈線性關系。

    1理論衰減規(guī)律

    本研究的基本原理是將地震破裂高度簡化地看成在完全彈性介質中發(fā)生的動態(tài)破裂過程[9]。 依據(jù)此觀點,此問題涉及固有長度尺度和受應力降與破裂特征長度尺度控制的輻射波場。 因此,衰減規(guī)律應該綜合以下參數(shù)。

    1.1位移和速度譜

    P波遠場位移譜可以利用歐米伽平方(omega-squared)模型[10-11]表示:

    (1)

    其中,ω是角頻率,ω0是拐角頻率,Ω0是低頻譜。 類似地,P波速度譜可以表示為:

    (2)

    因此,當位移譜振幅平緩,低于拐角頻率,且隨著ω-2衰減時,速度譜振幅呈現(xiàn)峰值,在拐角頻率輻射大部分能量。 拐角頻率與破裂尺度的關系為[11]:

    (3)

    其中,f0為拐角頻率,r是破裂半徑,CS是剪切波速,k是常數(shù)。 低頻譜與地震矩M0的關系表示式為:

    (4)

    其中,Uφθ為輻射圖案,F(xiàn)S為自由表面校正系數(shù),CP為P波波速,R為震源距離,ρ為密度。 對于圓形破裂,地震矩與應力降Δτ和破裂半徑的關系為[12]:

    (5)

    滯彈性和近地表衰減無法反映在方程(1)和(2)中。 的確,對于在地震預警中應用的震源距離(R<60 km),P波損耗很小,可以忽略[13-14]。 這將在下一部分進行核實。 另一方面,近地表效應不能被校正,除非能收集到地表響應的詳細數(shù)據(jù),這里將這些因素認為是觀測誤差的一部分。

    1.2衰減規(guī)律

    從P波地面運動均方根和地震矩、 應力降和震源距離的關系獲得表達式。 給定時間序列y(t)的均方根可以利用Parseval定理計算其頻譜Y(f):

    (6)

    其中,TS-P是數(shù)據(jù)區(qū)間,它等于S波和P波到時之差:

    (7)

    利用方程(1)和(2)中的位移和波速頻譜分別求解方程(6),算出位移均方根:

    (8a)

    和波速均方根:

    (8b)

    通過將方程(3)~(5)和(7)代入方程(8a)和(8b),drms和vrms的衰減關系可以根據(jù)應力降和地震矩表示:

    (9a)

    (9b)

    或可以根據(jù)應力降和破裂長度表示:

    (10a)

    (10b)

    (11)

    統(tǒng)計理論和模型預測可以得到給定時間序列內的均方根和絕對值的直接關系[15-17]。 因此,原則上,峰值和地面震動均方根的衰減相似。 許多地震學家曾在地面運動預測上采取此類推測[18-21]。 此后,經(jīng)驗證實,vrms和drms確實與峰值速度Pv和峰值位移Pd成比例,并將此應用于幾個地震預警算法中[4-5,7-8,22-25]。

    2理論與觀測

    在本部分中,根據(jù)觀測衰減,利用南加州和日本的地震數(shù)據(jù)驗證了理論衰減規(guī)律。 文中參數(shù)設定為:Uφθ=0.52,F(xiàn)S=2,ρ=2600 kg/m3,CS=3200 m/s,CP=5333 m/s,k=0.32[9],地震矩從Hanks和Kanamori[26]矩震級目錄中獲取。 利用以上參數(shù)得到,η和分別等于s/km和7.5×10-131/Pa。 最后一個參數(shù)是應力降,其值設定為應力降中值。 應力降將在應力降估算部分進行闡述。

    2.1數(shù)據(jù)

    動圖產生了震級偏差,例如日本中強地震(圖3中的圓形)的震級比加州小地震(圖3中的三角形)降低了許多。 因此,注意,數(shù)據(jù)區(qū)間在S波到達前終止(傳感器傾斜是不可能的),決定不采用高通濾波位移震動圖的方法,包括那些通過雙綜合加速度時間序列的數(shù)據(jù)。 采用信噪比(signal-to-noise ratios,SNRs)評估效果:

    圖1日本(a)和南加州(b)的地震(正方形)和臺站(三角形)的位置圖

    圖2根據(jù)震級和震源距離所形成的波形分布。 三角形表示加州綜合地震臺網(wǎng)記錄的寬頻帶速度,圓形表示K-NET和KiK-net地表加速度計記錄的加速度圖。 當應力降小于7.9 MPa(實線)時,數(shù)據(jù)區(qū)間大于破裂周期

    圖3drms的原始數(shù)據(jù)和高通濾波結果對比圖。 雙對數(shù)圖表示,在0.02 Hz時,drms高通濾波與原始drms的函數(shù)關系。 三角形表示加州綜合地震臺網(wǎng)記錄的寬頻帶速度,圓形表示K-NET和KiK-net地表加速度計記錄的加速度圖。 虛線和綠色實線分別是1:1線和最佳線性擬合線。 右下角為擬合系數(shù)和線性擬合標準差。 彩色表示震級大小

    (12)

    其中,AZi是垂直地面運動的原始數(shù)據(jù)(速度或加速度),l和m分別是數(shù)據(jù)和前置信號區(qū)間的樣本數(shù)量,后者可延長達20 s。 數(shù)據(jù)集排除掉信噪比小于20的記錄。

    2.2應力降評估

    近來,運用drms和vrms的衰減規(guī)律來表示應力降。 將衰減關系式(9b)除以(9a),求解應力降表達式[27-28]:

    (13)

    此表達式中的重要特征為括號內的關系式是破裂半徑特征長度r的倒數(shù)。 此外,這個方程式與廣泛應用于地震預警系統(tǒng)中的基于頻率的震級替代品τc密切相關[29]。 應力降的另一個表達式可通過將方程(13)中的M0代入方程(9a)中得到:

    (14)

    在方程(13)和(14)兩個表達式中,應力降是vrms,drms和CS的函數(shù); 然而,在方程(13)中,它也是地震矩的函數(shù); 在方程(14)中,它是和震源距離的函數(shù)。 因此,方程(13)的顯著優(yōu)勢是的獨立性。 利用方程(14)可以很好地將約束在較好的區(qū)域,因為距離不確定性遠小于地震矩的影響。

    利用方程(13)和(14)可以確定應力降,有別于之前利用遠場地面運動頻譜確定應力降。 然而,由于數(shù)據(jù)范圍遠小于破裂周期使得震源參數(shù)不可能通過數(shù)據(jù)區(qū)間進行可靠估算,因此,以上表達式不能任意應用于所有數(shù)據(jù)集。 為了篩選出與破裂周期相一致的數(shù)據(jù)點,單次地震的破裂周期可以通過計算破裂直徑和破裂速度的比率估算,后者設為0.9CS:

    (15)

    在估算破裂直徑時,可以尋找一種以過濾掉一些質量好的數(shù)據(jù)點為代價,高估(overestimating)破裂持續(xù)時間。 因此,采用相對小的應力降1 MPa,通過方程(5)評估破裂直徑。 本研究中67%震動圖分析滿足條件TS-P>Tr。 對于震動圖的子集,利用方程(13)和(14)計算應力降。 圖4顯示單個數(shù)據(jù)點(三角形和圓形分別表示加州和日本地震)和平均事件(用╳標識)的震級目錄和應力降函數(shù)關系。 圖4顯示了詳細的應力降中值。 這些值與之前利用完全不同的方法估算的應力降一致[30-35]。 下一部分,利用與單獨數(shù)據(jù)點中值(7.9 MPa)相一致的應力降和將應力降評估值代入方程(5)得到的破裂直徑特征計算理論衰減曲線。

    圖4應力降評估。 圖(a)和圖(b)分別表示方程(13)與方程(14)計算結果與震級的關系。 三角形和圓形分別表示加州和日本單臺評估。 ╳表示平均事件評估,其表示具有4個或更多地震計記錄事件的一個子集。 圖中顯示了平均和單一數(shù)據(jù)點的中位值。 虛線表示單個臺站評估的中位值。 圖中沒有描述數(shù)據(jù)點為TS-P

    2.3觀測衰減值

    圖5和圖6分別是利用比TS-P更長(空心符號)或更短(實心符號)的破裂周期相一致的數(shù)據(jù)點構成的vrms和drms的觀測值和預測值對比圖。 破裂周期通過采用方程(13)中應力降7.9 MPa估算,破裂周期的數(shù)據(jù)點比TS-P圖短,位于圖2中震級與距離對數(shù)關系圖實線7.9 MPa之下。 在圖5中,vrms和drms與震源距離呈函數(shù)關系,理論曲線為方程(9a)和(9b)中vrms和drms的表達式。 圖6分別為vrms和drms與R/r和R/r5/3的函數(shù)關系,理論曲線為方程(10a)和(10b)中vrms和drms的表達式。 利用震源半徑劃分距離軸能降低圖5中一大群數(shù)據(jù)點高度離散狀態(tài),使之成為圖6中高度集中于其理論曲線周圍(綠虛線)的狹窄范圍內。 毫無意外,呈現(xiàn)在理論曲線之下的數(shù)據(jù)點子集正好位于那些破裂周期遠大于TS-P區(qū)間(圖5和圖6中的空心符號)。 為了進一步評估理論值與觀測值的一致性,圖6a中的數(shù)據(jù)利用下式擬合:

    圖5(a) 地面運動位移與震源距離的均方根; (b) 波速與震源距離的均方根。 比TS-P長的破裂周期用空心符號表示,比TS-P短的破裂周期用實心符號表示。 三角形和圓形分別代表加州和日本的數(shù)據(jù)。 彩色代表震級大小。 實線是與方程(9a)和(9b)相一致的理論曲線

    log10(drms)=ad+bdlog10(R/r5/3)

    (16a)

    圖6b中的數(shù)據(jù)利用下式擬合:

    log10(vrms)=av+bvlog10(R/r)

    (16b)

    其中,ad,bd,av,bv為擬合系數(shù),圖6左下角顯示具體數(shù)值。 以上衰減規(guī)律僅能用于當滯彈性衰減較小且數(shù)據(jù)區(qū)間并非遠小于破裂周期時有效。 為確保符合后者條件,位于圖2中7.9 MPa之上的數(shù)據(jù)點(圖5和圖6中空心符號)需從圖6的線性擬合中排除。 圖6中紅線和綠線分別表示R<30 km和R<60 km的最佳擬合曲線。 類似地,不同距離范圍的最佳擬合曲線和理論曲線具有較好的一致性。 可以推斷方程(1)和(2)中滯彈性衰減的影響是可校正的,假設TS-P大于破裂周期,理論衰減規(guī)律對于距離達到至少60 km是有效的。 觀察圖6中許多位于理論線之下的空心符號,結果表明震源大小不能利用遠小于破裂周期的數(shù)據(jù)區(qū)間進行估算。

    圖6(a) 地面運動位移與R/r5/3的均方根; (b) 波速與R/r的均方根。 比TS-P長的破裂周期用空心符號表示,比TS-P短的破裂周期用實心符號表示。 三角形和圓形分別代表加州和日本的數(shù)據(jù)。 彩色表示震級大小。 藍色虛線是與方程(10a)和(10b)相一致的理論曲線。 紅線和綠線分別表示R≤30 km和R≤60 km的線性擬合結果。 左下角顯示方程(16a)和(16b)的擬合系數(shù)和標準差

    圖7顯示日本和加州的衰減情況,并與理論衰減進行對比。 加州和日本具有相似的地面運動衰減,表明衰減規(guī)律在這兩個區(qū)域近似可用,推斷可將其推廣到相似的地質構造區(qū)域。 將最新獲取的衰減規(guī)律整合入地震預警算法,對于那些利用離線數(shù)據(jù)進行參數(shù)調整,觀測數(shù)據(jù)的數(shù)量和質量都受到局限的地方是最有利的。

    圖7對比加州(點線)和日本(虛線)地面運動衰減。 灰線為log(drms)與log(R/r5/3)的最佳線性擬合結果(左下角),黑線為log(vrms)與log(R/r)的最佳線性擬合結果(右上角)。 實灰線和實黑線分別是與方程(10a)和(10b)相一致的理論曲線。 利用數(shù)據(jù)滿足TS-P>Tr要求繼續(xù)計算擬合。 圖中顯示了最佳擬合系數(shù)

    3衰減規(guī)律應用于地震預警

    3.1實時震級測定

    對3個獨立震級估算結果進行評價。 前兩個表達式通過求解理論衰減規(guī)律[方程(9)]獲取地震矩:

    (17a)

    (17b)

    第3個表達式通過將方程(13)代入方程(9a)獲取并求解地震矩:

    (18)

    圖8震級預測值和觀測值差異分布圖。 利用應力降獨立方程組(17)和(18)計算預測震級,在這些方程中設置應力降等于圖4中7.9 MPa中值。 (a) 用drms和vrms利用方程(18)得到震級。 (b) 用drms利用方程(17a)得到震級。 (c) 用vrms利用方程(17b)得到震級。 淡色柱形表示整個數(shù)據(jù)庫,深色柱形表示數(shù)據(jù)子集,其數(shù)據(jù)區(qū)間大于破裂周期,即為TS-P>Tr。 圖中表示震級差別的平均值和標準差與大于破裂周期的數(shù)據(jù)區(qū)間相一致

    利用P波衰減規(guī)律進行實時震級測定,其準確性受實時震源距離準確性影響。 檢測方程(17)和(18)表明,方程(17b)比方程(17a)和(18)對距離敏感度更高。 然而,由于矩震級關系的對數(shù)特征和實時距離測定的準確性[36-38],可以推斷求解這些方程估算震級并不會很大程度地受距離誤差的影響。

    3.2從vrms和drms到Pv和Pd

    (19a)

    (19b)

    此經(jīng)驗關系與統(tǒng)計理論和模型相一致[15-17],這些方法曾用于地面運動預測方程[18-21]。 從加州和日本的數(shù)據(jù)擬合結果到回歸曲線相似性,可以推斷以上峰值與均方根關系幾乎不受場地條件和地質背景的影響。 將以上經(jīng)驗關系式代入方程(9a)和(9b)或者(10a)和(10b)獲取Pd和Pv的半經(jīng)驗衰減規(guī)律。

    3.3理論τc

    信號前幾秒的特征周期τc是基于頻率實時震級的替代品,其與Pd一起,構成地震預警系統(tǒng)Pd-τc的基礎[4,29,40]。 本項研究一個有趣的副產品是揭示了τc與drms/vrms成比例:

    圖9雙對數(shù)圖顯示地面運動峰值與地面運動均方根的函數(shù)關系。 (a)log(Pd)和log(drms)。 (b)log(Pv)與log(vrms)。 綠線表示log(drms)=ad+bd×log(Pd)和log(vrms)=av+bv×log(Pv)的最小平方擬合,每幅圖右下角為擬合系數(shù)和標準方差。 三角形和圓形分別表示加州和日本的數(shù)據(jù)。 彩色表示震級大小

    (20)

    其中,u和v分別是地面位移和速度,t0是P波初至,Δt是預先設定的區(qū)間。 把衰減規(guī)律公式(8a)除以(8b),表明τc與拐角頻率倒數(shù)相等,因此解釋了其可以表征地震大小的原因。 將上式代入方程(13)得到τc理論表達式:

    (21)

    圖10是τc理論值和觀測值的對比。 盡管海量數(shù)據(jù)離散以及τc理論值和觀測值間一致性較差,但τc的對數(shù)和震級目錄仍具有統(tǒng)計意義。 最近,Ziv[41]提出了τlog,一個基于時頻的實時震級替代品,通過求波速譜對數(shù)平均值獲取,并表明其與震級的相關性優(yōu)于τc。

    4結論

    利用日本和加州當?shù)氐卣痱炞C理論衰減規(guī)律。 震源距離達30 km和60 km的地面運動觀測值和預測值具有很好的一致性。 因此,推斷在這些距離范圍內可校正滯彈性效應的影響。 結果表明,vrms/drms與特征破裂長度成比例,應力降與地震矩和drms/vrms的立方成函數(shù)關系[方程(13)]。 利用這個結果,可以獲取應力降中值7.9 MPa。

    圖10特征周期τc與震級的函數(shù)關系。 三角形和圓形分別代表加州和日本的數(shù)據(jù)。 利用方程(21)計算的理論直線與不同應力降相一致。 大或小于數(shù)據(jù)區(qū)間TS-P的破裂周期分別用空心和實心符號表示。 粗線表示實心符號的最小平方擬合log(τc)=a+b×M,其中M是編目震級,右下角為擬合系數(shù)和標準偏差

    為了評價這些結果在地震預警應用中的效果,利用理論衰減規(guī)律(9a)和(9b)進行震級估算,如設置這些方程中應力降等于中值,或是將方程(13)替換成(9a)。 當數(shù)據(jù)區(qū)間大于破裂周期時,震級預測值和觀測值呈現(xiàn)很好的一致性。 反之,預測震級將小于觀測值。 由于drms比vrms更容易受地震矩影響,利用方程(17a)和(18)比方程(17b)能獲得更好的震級預測結果。

    5數(shù)據(jù)和來源

    本研究所用數(shù)據(jù)來自南加州地震數(shù)據(jù)中心(http:∥scedc.caltech. edu/,截止到2015年1月),以及來自K-NET和KiK-net強震動觀測網(wǎng)(http:∥www.kyoshin.bosai.go.jp/,截止到2015年1月)。

    文獻來源: Itzhak Lior,Alon Ziv,Raul Madariaga. P-wave attenuation with implications for earthquake early warning. Bull. Seismol. Soc. Amer.,2016,106(1): 13-22. doi:10.1785/0120150087

    (福建省地震局王林譯; 黃宏生校)

    (譯者電子信箱,王林: wl_0117@163.com)

    參 考 文 獻

    [1] Nakamura Y. Development of earthquake early-warning system for the Shinkansen,some recent earthquake engineering research and practical in Japan. The Japanese National Committee of the International Association for Earthquake Engineering.1984: 224-238

    [2] Nakamura Y. On the Urgent Earthquake Detection and Alarm System (UrEDAS). Proc. 9thWorld Conference on Earthquake Engineering VII. Tokyo,Japan,2-9 August 1988: 673-678

    [3] Allen R M,Kanamori H. The potential for earthquake early warning in southern California. Science, 2003, 300(5620): 786-789

    [4] Wu Y M,Kanamori H. Rapid assessment of damaging potential of earthquakes in Taiwan from the beginning of P waves. Bull. Seismol. Soc. Amer., 2005, 95(3): 1181-1185

    [5] Cua G,Heaton T. The virtual seismologist (VS) method: A Bayesian approach to earthquake early warning. in Seismic Early Warning. Gasparini P,Manfredi G,Zschau J (Editors). New York: Springer,2007

    [6] Allen R M, Brown H,Hellweg M,et al. Real-time earthquake detection and hazard assessment by ElarmS across California. Geophys. Res. Lett., 2009, 36(5): 277-291. L00B08. doi:10.1029/2008GL0367

    66

    [7] B?se M,Hauksson E,Solanki K,et al. A new trigger criterion for improved real-time performance of on-site earthquake early warning in southern California. Bull. Seismol. Soc. Amer., 2009, 99(2A): 897-905

    [8] Zollo A,Iannaccone G,Lancieri M,et al. Earthquake early warning system in southern Italy: Methodologies and performance evaluation. Geophys. Res. Lett., 2009, 36(4): 121-136. L00B07. doi:10.1029/2008GL036689

    [9] Madariaga R. Dynamics of an expanding circular fault. Bull. Seismol. Soc. Amer., 1976, 66(3): 639-666

    [10] Aki K. Scaling law of seismic spectrum. J. Geophys. Res., 1967, 72(4): 1217-1231

    [11] Brune J N. Tectonic stress and the spectra of seismic shear waves from earthquakes. J. Geophys. Res., 1970, 75(75): 4997-5009

    [12] Eshelby J D. The determination of the elastic field of an ellipsoidal inclusion,and related problems. Proc. Math. Phys. Sci., 1957, 241(1226): 376-396

    [13] Wu Y M,Allen R M,Wu C F. RevisedMLdetermination for crustal earthquakes in Taiwan. Bull. Seismol. Soc. Amer., 2005, 95(6): 2517-2524

    [14] Wu Y M,Zhao L. Magnitude estimation using the first three seconds P-wave amplitude in earthquake early warning. Geophys. Res. Lett., 2006, 33(16): 271-284. L16312. doi:10.1029/2006GL026871

    [15] Cartwright D E,Longuet-Higgins M S. The statistical distribution of the maxima of a random function. Proc. Math. Phys. Sci., 1956, 237(1209): 212-232

    [16] Davenport A G. Note on the distribution of the largest value of a random function with application to gust loading. Proc. Inst. Civ. Eng., 1964, 28(2): 187-196

    [17] Clough R W,Penzien J. Dynam. Struct. McGraw-Hill. New York,1975: 634

    [18] Hanks T C,McGuire R K. The character of high frequency strong ground motion. Bull. Seismol. Soc. Amer., 1981, 71(6): 2071-2095

    [19] Boore D M. Stochastic simulation of high-frequency ground motions based on seismological models of the radiated spectra. Bull. Seismol. Soc. Amer., 1983, 73(6): 1865-1894

    [20] Boore D M. Simulation of ground motion using the stochastic method. Pure Appl. Geophys., 2003, 160(3-4): 635-676

    [21] Boore D M,Joyner W B. A note on the use of random vibration theory to predict peak amplitudes of transient signals. Bull. Seismol. Soc. Amer., 1984, 74(5): 2035-2039

    [22] Lancieri M,Zollo A. A Bayesian approach to the real-time estimation of magnitude from the early P and S wave displacement peaks. J. Geophys. Res., 2008, 113(113): 178-196

    [23] Brown H M,Allen R M,Hellweg M,et al. Development of the ElarmS methodology for earthquake early warning: Real time application in California and offline testing in Japan. Soil Dynam. Earthq. Eng., 2011, 31(2): 188-200

    [24] Lancieri M,F(xiàn)uenzalida A,Ruiz S,et al. Magnitude scaling of early-warning parameters for theMW7.8 Tocopilla,Chile,earthquake and its aftershocks. Bull. Seismol. Soc. Amer., 2011, 101(2): 447-463

    [25] Sadeh M,Ziv A,Wust-Bloch H. Real-time magnitude proxies for earthquake early warning in Israel. Geophys. J. Int., 2014, 196(2): 939-950

    [26] Hanks T C,Kanamori H. A moment magnitude scale. J. Geophys. Res., 1979, 84(B5): 2348-2350

    [27] Andrews D J. Objective determination of source parameters and similarity of earthquakes of different size. in Earthquake Source Mechanics. Das S,Boatwright J,Scholz C H(Editors). Washington D.C: American Geophysical Union,1986: 259-268

    [28] Snoke J A. Stable determination of (Brune) stress drops. Bull. Seismol. Soc. Amer., 1987, 77(2): 530-538

    [29] Wu Y M,Kanamori H. Experiment on an onsite early warning method for the Taiwan early warning system. Bull. Seismol. Soc. Amer., 2005, 95(1): 347-353

    [30] Hanks T C,Thatcher W. A graphical representation of seismic source parameters. J. Geophys. Res., 1972, 77(23): 4393-4405

    [31] Mayeda K,Walter W R. Moment,energy,stress drop,and source spectra of western United States earthquakes from regional coda envelopes. J. Geophys. Res., 1996, 101(B5): 11195-11208

    [32] Oth A,Bindi D,Parolai S,et al. Earthquake scaling characteristics and the scale-(in)dependence of seismic energy-to-moment ratio: Insights from KiK-net data in Japan. Geophys. Res. Lett., 2010, 37(19): 470-479. doi:10.1029/2010GL044572

    [33] Baltay A,Ide S,Prieto G,et al. Variability in earthquake stress drop and apparent stress. Geophys. Res. Lett., 2011, 38. L06303. doi:10.1029/2011GL046698

    [34] Baltay A S,Hanks T C,Beroza G C. Stable stress-drop measurements and their variability: Implications for ground-motion prediction. Bull. Seismol. Soc. Amer., 2013,103(1): 211-222

    [35] Abercrombie R E. Comparison of direct and coda wave stress drop measurements for the Wells,Nevada,earthquake sequence. J. Geophys. Res., 2013, 118(4): 1458-1470

    [36] Satriano C,Lomax A,Zollo A. Real-time evolutionary earthquake location for seismic early warning. Bull. Seismol. Soc. Amer., 2008, 98(3): 1482-1494

    [37] Kurzon I,Vernon F L,Rosenberger A,et al. Real-time automatic detectors of P and S waves using singular value decomposition. Bull. Seismol. Soc. Amer., 2014, 104(4): 1696-1708

    [38] Eisermann A S,Ziv A,Wust-Bloch G H. Real-time back azimuth for earthquake early warning. Bull. Seismol. Soc. Amer., 2015, 105(4): 2274-2285. doi:10.1785/0120140298

    [39] Doi K. The operation and performance of earthquake early warnings by the Japan Meteorological Agency. Soil Dynam. Earthq. Eng., 2011, 31(2): 119-126

    [40] Kanamori H. Real-time seismology and earthquake damage mitigation. Ann. Rev. Earth Planet. Sci., 2005, 33(5): 195-214

    [41] Ziv A. New frequency-based real-time magnitude proxy for earthquake early warning. Geophys. Res. Lett., 2014, 41(20): 7035-7040

    中圖分類號:P315.9;

    文獻標識碼:A;

    doi:10.3969/j.issn.0235-4975.2016.02.009

    * 收稿日期:2016-02-03; 采用日期: 2016-02-15。

    猜你喜歡
    震級方根震源
    方根拓展探究
    基于累積絕對位移值的震級估算方法
    地震后各國發(fā)布的震級可能不一樣?
    新震級國家標準在大同臺的應用與評估
    山西地震(2020年1期)2020-04-08 07:34:26
    震源的高返利起步
    均方根嵌入式容積粒子PHD 多目標跟蹤方法
    自動化學報(2017年2期)2017-04-04 05:14:28
    揭開心算方根之謎
    可控震源地震在張掖盆地南緣逆沖斷裂構造勘探中的應用
    華北地質(2015年3期)2015-12-04 06:13:25
    中國地震臺網(wǎng)面波震級與矩震級的統(tǒng)計關系
    同步可控震源地震采集技術新進展
    淫秽高清视频在线观看| bbb黄色大片| 人人妻,人人澡人人爽秒播| xxxhd国产人妻xxx| 99精品久久久久人妻精品| 亚洲av日韩精品久久久久久密| 99久久人妻综合| svipshipincom国产片| 长腿黑丝高跟| 一进一出抽搐gif免费好疼 | 午夜免费鲁丝| av福利片在线| 高清黄色对白视频在线免费看| 国内毛片毛片毛片毛片毛片| 乱人伦中国视频| 亚洲国产欧美网| 成人18禁在线播放| 性少妇av在线| 1024视频免费在线观看| 国产精品九九99| 啦啦啦免费观看视频1| 国产高清视频在线播放一区| 欧美不卡视频在线免费观看 | 国产亚洲欧美98| 国产熟女午夜一区二区三区| 免费看a级黄色片| 黄色怎么调成土黄色| 99国产精品一区二区蜜桃av| 亚洲精品久久成人aⅴ小说| 777久久人妻少妇嫩草av网站| 亚洲成国产人片在线观看| 亚洲国产欧美一区二区综合| 91精品三级在线观看| 国产xxxxx性猛交| 成人国产一区最新在线观看| 777久久人妻少妇嫩草av网站| 韩国av一区二区三区四区| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品爽爽va在线观看网站 | 亚洲av美国av| 女性被躁到高潮视频| 精品久久蜜臀av无| 国产精品久久视频播放| 亚洲av美国av| 亚洲少妇的诱惑av| 老司机在亚洲福利影院| 人人澡人人妻人| 免费在线观看黄色视频的| 男女之事视频高清在线观看| 午夜免费观看网址| 伊人久久大香线蕉亚洲五| 午夜两性在线视频| 国产伦一二天堂av在线观看| www.熟女人妻精品国产| 国产精品成人在线| 国产亚洲欧美在线一区二区| 天天影视国产精品| 欧美日韩黄片免| 亚洲av片天天在线观看| 脱女人内裤的视频| 国产精品永久免费网站| 成人三级做爰电影| 国产极品粉嫩免费观看在线| 国产一区二区三区在线臀色熟女 | 后天国语完整版免费观看| 国产精品 国内视频| 精品福利永久在线观看| 天天添夜夜摸| 日韩三级视频一区二区三区| 18禁国产床啪视频网站| 国产av一区二区精品久久| 免费少妇av软件| 午夜影院日韩av| 国产三级在线视频| 久久人人97超碰香蕉20202| 免费观看精品视频网站| 色综合欧美亚洲国产小说| 欧美成人性av电影在线观看| 在线视频色国产色| 女性被躁到高潮视频| av国产精品久久久久影院| 日本撒尿小便嘘嘘汇集6| 一区二区三区精品91| 俄罗斯特黄特色一大片| 欧美激情极品国产一区二区三区| 国产欧美日韩一区二区三| 90打野战视频偷拍视频| 亚洲自拍偷在线| 大型黄色视频在线免费观看| 亚洲av成人不卡在线观看播放网| 丰满饥渴人妻一区二区三| 亚洲色图av天堂| 亚洲人成网站在线播放欧美日韩| 丁香六月欧美| 中文字幕av电影在线播放| 老司机福利观看| 69精品国产乱码久久久| 人人妻人人爽人人添夜夜欢视频| 高潮久久久久久久久久久不卡| 国产欧美日韩综合在线一区二区| 国产免费现黄频在线看| 国产午夜精品久久久久久| 成人三级黄色视频| 一边摸一边做爽爽视频免费| 亚洲一区二区三区不卡视频| 亚洲欧洲精品一区二区精品久久久| 国产亚洲精品久久久久5区| 国产亚洲精品第一综合不卡| 激情视频va一区二区三区| 91九色精品人成在线观看| 中文字幕最新亚洲高清| 99精品久久久久人妻精品| 免费高清在线观看日韩| 日本五十路高清| 老熟妇乱子伦视频在线观看| www.熟女人妻精品国产| a在线观看视频网站| 五月开心婷婷网| 欧美黑人精品巨大| av中文乱码字幕在线| 99国产极品粉嫩在线观看| 成人特级黄色片久久久久久久| 精品国产国语对白av| 一本综合久久免费| 国产成人精品无人区| 欧美一区二区精品小视频在线| 国产精品 欧美亚洲| 国产av精品麻豆| 日韩欧美一区二区三区在线观看| 岛国在线观看网站| 成人18禁在线播放| 一区二区三区精品91| 欧美精品一区二区免费开放| 国产精品 国内视频| 国产一区二区激情短视频| 亚洲熟妇熟女久久| 成人三级黄色视频| 欧美性长视频在线观看| 在线观看日韩欧美| 露出奶头的视频| 久久人人精品亚洲av| 久久精品亚洲熟妇少妇任你| 国内久久婷婷六月综合欲色啪| 十分钟在线观看高清视频www| 99riav亚洲国产免费| 丝袜美腿诱惑在线| 嫩草影院精品99| 国产蜜桃级精品一区二区三区| 男人舔女人的私密视频| 好男人电影高清在线观看| 欧美精品亚洲一区二区| 在线国产一区二区在线| 一区二区三区国产精品乱码| www国产在线视频色| 波多野结衣一区麻豆| 亚洲九九香蕉| 免费搜索国产男女视频| 亚洲性夜色夜夜综合| 日本黄色日本黄色录像| 一夜夜www| 宅男免费午夜| 国产精品综合久久久久久久免费 | 身体一侧抽搐| 欧美人与性动交α欧美软件| 午夜免费鲁丝| 久久热在线av| 黄色毛片三级朝国网站| 午夜成年电影在线免费观看| 成人亚洲精品av一区二区 | 免费日韩欧美在线观看| 不卡av一区二区三区| 国产精品永久免费网站| 日韩欧美三级三区| 久久人人爽av亚洲精品天堂| 真人一进一出gif抽搐免费| 亚洲黑人精品在线| 男女做爰动态图高潮gif福利片 | 日本a在线网址| 十八禁网站免费在线| 琪琪午夜伦伦电影理论片6080| 免费在线观看完整版高清| 国产av一区二区精品久久| 亚洲久久久国产精品| 少妇被粗大的猛进出69影院| 亚洲国产欧美日韩在线播放| 欧美精品亚洲一区二区| 99精品在免费线老司机午夜| 精品久久蜜臀av无| 日韩大码丰满熟妇| 成年人免费黄色播放视频| 男人舔女人下体高潮全视频| 嫩草影视91久久| 嫩草影视91久久| 久久久久久人人人人人| 亚洲色图av天堂| 侵犯人妻中文字幕一二三四区| 色婷婷久久久亚洲欧美| 岛国在线观看网站| 欧美成人午夜精品| av电影中文网址| 免费在线观看日本一区| 国产精品 欧美亚洲| 成人三级做爰电影| 欧美成人免费av一区二区三区| 精品久久久久久成人av| 在线观看午夜福利视频| 国产成人精品无人区| 女生性感内裤真人,穿戴方法视频| 亚洲成人国产一区在线观看| 国产深夜福利视频在线观看| 午夜视频精品福利| 亚洲国产毛片av蜜桃av| 久久久久国产精品人妻aⅴ院| 777久久人妻少妇嫩草av网站| 亚洲欧美日韩高清在线视频| 色综合欧美亚洲国产小说| 后天国语完整版免费观看| 黄频高清免费视频| 久久精品影院6| 欧美黄色片欧美黄色片| 久久久精品欧美日韩精品| 91麻豆精品激情在线观看国产 | 国产精品免费一区二区三区在线| 最好的美女福利视频网| 成人国产一区最新在线观看| 麻豆成人av在线观看| 亚洲精品成人av观看孕妇| 国产乱人伦免费视频| 又黄又爽又免费观看的视频| 午夜视频精品福利| 好男人电影高清在线观看| 少妇粗大呻吟视频| 亚洲 欧美 日韩 在线 免费| 又紧又爽又黄一区二区| 50天的宝宝边吃奶边哭怎么回事| 免费久久久久久久精品成人欧美视频| 又大又爽又粗| 亚洲中文av在线| 日韩欧美一区视频在线观看| 91成人精品电影| www.自偷自拍.com| 999久久久精品免费观看国产| 成人永久免费在线观看视频| 亚洲国产欧美网| 在线观看66精品国产| 久久国产精品影院| 一边摸一边抽搐一进一出视频| 国产成人欧美| 国产高清视频在线播放一区| 99国产精品一区二区三区| 咕卡用的链子| 色综合婷婷激情| 自拍欧美九色日韩亚洲蝌蚪91| 日韩有码中文字幕| 婷婷丁香在线五月| 母亲3免费完整高清在线观看| 国产有黄有色有爽视频| 午夜福利在线观看吧| 人成视频在线观看免费观看| 天堂中文最新版在线下载| 中国美女看黄片| 日韩大码丰满熟妇| videosex国产| 午夜福利在线观看吧| xxx96com| 国产三级在线视频| 身体一侧抽搐| 九九在线视频观看精品| 午夜免费成人在线视频| 日本与韩国留学比较| eeuss影院久久| 最近视频中文字幕2019在线8| 亚洲国产欧美人成| 99热只有精品国产| 国产精品自产拍在线观看55亚洲| 国产成年人精品一区二区| 噜噜噜噜噜久久久久久91| 精品一区二区三区视频在线| 国产精品嫩草影院av在线观看 | 久久久久国内视频| 国产黄色小视频在线观看| 成人性生交大片免费视频hd| 日韩 亚洲 欧美在线| 好男人在线观看高清免费视频| 免费在线观看亚洲国产| 精品熟女少妇八av免费久了| 热99re8久久精品国产| 午夜精品一区二区三区免费看| 小蜜桃在线观看免费完整版高清| 午夜两性在线视频| 精品一区二区三区人妻视频| 日本免费一区二区三区高清不卡| 五月玫瑰六月丁香| 黄色女人牲交| 人人妻,人人澡人人爽秒播| 一级黄色大片毛片| 99国产极品粉嫩在线观看| 日本黄大片高清| 自拍偷自拍亚洲精品老妇| 亚洲精品亚洲一区二区| av在线蜜桃| 国产欧美日韩一区二区三| 国内揄拍国产精品人妻在线| 久久精品国产清高在天天线| 69av精品久久久久久| 别揉我奶头~嗯~啊~动态视频| 日本撒尿小便嘘嘘汇集6| 成熟少妇高潮喷水视频| 韩国av一区二区三区四区| 成人永久免费在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看| 久久热精品热| 国内少妇人妻偷人精品xxx网站| 亚洲熟妇中文字幕五十中出| 热99在线观看视频| av国产免费在线观看| 在线观看一区二区三区| 淫秽高清视频在线观看| 国产高清有码在线观看视频| bbb黄色大片| 午夜老司机福利剧场| 观看美女的网站| 少妇的逼好多水| 男女之事视频高清在线观看| 丝袜美腿在线中文| 俺也久久电影网| 亚洲精品456在线播放app | www.色视频.com| 精品乱码久久久久久99久播| 亚洲av第一区精品v没综合| 高清毛片免费观看视频网站| 国产在视频线在精品| 亚洲精品在线美女| 变态另类成人亚洲欧美熟女| 久久精品人妻少妇| 最新中文字幕久久久久| 欧美一区二区国产精品久久精品| 免费观看人在逋| 午夜免费激情av| 国产精品伦人一区二区| 精品99又大又爽又粗少妇毛片 | 欧美一区二区国产精品久久精品| 可以在线观看的亚洲视频| 欧美另类亚洲清纯唯美| 中文字幕高清在线视频| 又黄又爽又免费观看的视频| 久久久久久九九精品二区国产| 在线观看免费视频日本深夜| 两个人的视频大全免费| 精品一区二区三区人妻视频| 久久99热这里只有精品18| 国产不卡一卡二| 日韩欧美三级三区| 国产成人福利小说| 精品久久国产蜜桃| 午夜免费激情av| 欧美黄色淫秽网站| 一本精品99久久精品77| 91在线观看av| 国产三级在线视频| 精品乱码久久久久久99久播| 久久久成人免费电影| 亚洲国产精品成人综合色| 国产成人av教育| 99精品在免费线老司机午夜| 赤兔流量卡办理| 全区人妻精品视频| 欧美日本亚洲视频在线播放| 国产欧美日韩精品亚洲av| 又爽又黄a免费视频| 久久国产精品人妻蜜桃| 99视频精品全部免费 在线| 成人毛片a级毛片在线播放| 3wmmmm亚洲av在线观看| 天堂动漫精品| 亚洲成av人片免费观看| 性欧美人与动物交配| 国产av一区在线观看免费| 国语自产精品视频在线第100页| 久久精品人妻少妇| 国产一区二区亚洲精品在线观看| 欧美成人性av电影在线观看| 久99久视频精品免费| 日日摸夜夜添夜夜添小说| 欧美最新免费一区二区三区 | 99久久精品国产亚洲精品| av黄色大香蕉| 国产乱人视频| 日本 av在线| 精品久久国产蜜桃| 欧美色欧美亚洲另类二区| 高潮久久久久久久久久久不卡| 一区二区三区高清视频在线| 搡老岳熟女国产| 免费高清视频大片| 麻豆成人午夜福利视频| 高清在线国产一区| 亚洲 欧美 日韩 在线 免费| 亚洲美女黄片视频| 国产av不卡久久| 久99久视频精品免费| 99热这里只有是精品50| 国产爱豆传媒在线观看| 国产高清视频在线播放一区| 久久精品久久久久久噜噜老黄 | 国产探花在线观看一区二区| 精品福利观看| 精品人妻1区二区| 国产av不卡久久| 免费在线观看亚洲国产| 小蜜桃在线观看免费完整版高清| 精品人妻偷拍中文字幕| 少妇丰满av| 床上黄色一级片| 午夜激情福利司机影院| 啦啦啦观看免费观看视频高清| 在线十欧美十亚洲十日本专区| 国产一区二区在线观看日韩| 久久人人爽人人爽人人片va | 中文字幕av成人在线电影| 神马国产精品三级电影在线观看| 一级毛片久久久久久久久女| 国产中年淑女户外野战色| 日日摸夜夜添夜夜添av毛片 | 日本一本二区三区精品| 国产成人a区在线观看| 精品福利观看| 国产精品一区二区性色av| 色在线成人网| 最近视频中文字幕2019在线8| 国产成人影院久久av| 18禁黄网站禁片午夜丰满| 亚洲无线观看免费| 日韩中文字幕欧美一区二区| 九九热线精品视视频播放| 十八禁人妻一区二区| 欧美成人免费av一区二区三区| 小说图片视频综合网站| 男女做爰动态图高潮gif福利片| 精品欧美国产一区二区三| www.熟女人妻精品国产| 在线十欧美十亚洲十日本专区| 免费av毛片视频| 在线观看午夜福利视频| 成人精品一区二区免费| 国产高清视频在线播放一区| 日日干狠狠操夜夜爽| 亚洲人成伊人成综合网2020| 真人做人爱边吃奶动态| 97人妻精品一区二区三区麻豆| 色综合站精品国产| 男人狂女人下面高潮的视频| 国产午夜精品论理片| 草草在线视频免费看| 成人午夜高清在线视频| 国产精品伦人一区二区| 亚洲精品日韩av片在线观看| 精品久久久久久久久亚洲 | 在线观看舔阴道视频| 国内揄拍国产精品人妻在线| 亚洲性夜色夜夜综合| 亚洲精华国产精华精| 狂野欧美白嫩少妇大欣赏| 尤物成人国产欧美一区二区三区| 又黄又爽又刺激的免费视频.| 99视频精品全部免费 在线| 亚洲美女视频黄频| 亚洲av第一区精品v没综合| 天堂网av新在线| 日本五十路高清| 最近在线观看免费完整版| 啦啦啦观看免费观看视频高清| 亚洲自拍偷在线| 18美女黄网站色大片免费观看| 美女高潮的动态| 在线观看66精品国产| 男女做爰动态图高潮gif福利片| 午夜激情福利司机影院| 日日摸夜夜添夜夜添av毛片 | 亚洲午夜理论影院| ponron亚洲| 99国产精品一区二区蜜桃av| 精品日产1卡2卡| 一本综合久久免费| 久久久久久久久久黄片| 毛片女人毛片| 一个人看视频在线观看www免费| eeuss影院久久| 很黄的视频免费| 美女被艹到高潮喷水动态| 99久久久亚洲精品蜜臀av| 色综合亚洲欧美另类图片| 热99re8久久精品国产| 精品日产1卡2卡| 伦理电影大哥的女人| 久久精品夜夜夜夜夜久久蜜豆| 免费看美女性在线毛片视频| 99热这里只有是精品在线观看 | 国产av在哪里看| 国产精品日韩av在线免费观看| 国产av一区在线观看免费| 亚洲第一区二区三区不卡| 日日干狠狠操夜夜爽| 国产成人啪精品午夜网站| 97碰自拍视频| 免费大片18禁| 欧美3d第一页| 欧美不卡视频在线免费观看| 亚洲va日本ⅴa欧美va伊人久久| 久久婷婷人人爽人人干人人爱| 久久久久久久精品吃奶| 国产精品一区二区免费欧美| 九九在线视频观看精品| 久久这里只有精品中国| 岛国在线免费视频观看| 亚洲人成网站高清观看| 每晚都被弄得嗷嗷叫到高潮| 色吧在线观看| 日韩亚洲欧美综合| 中文字幕av在线有码专区| 国产精品98久久久久久宅男小说| 自拍偷自拍亚洲精品老妇| 免费看美女性在线毛片视频| 脱女人内裤的视频| 亚洲乱码一区二区免费版| 好男人在线观看高清免费视频| 免费一级毛片在线播放高清视频| 国产欧美日韩一区二区三| 一本久久中文字幕| 免费人成在线观看视频色| 欧美+亚洲+日韩+国产| 搡老妇女老女人老熟妇| 小蜜桃在线观看免费完整版高清| 美女大奶头视频| 一区二区三区免费毛片| 色av中文字幕| 成人永久免费在线观看视频| 亚洲精品成人久久久久久| 国产一区二区三区在线臀色熟女| 久久国产精品人妻蜜桃| 久久热精品热| 嫁个100分男人电影在线观看| 欧美潮喷喷水| 丰满人妻一区二区三区视频av| 久久久国产成人免费| 又紧又爽又黄一区二区| 麻豆国产97在线/欧美| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 91麻豆精品激情在线观看国产| 日本一二三区视频观看| 级片在线观看| 日韩中文字幕欧美一区二区| 有码 亚洲区| 性欧美人与动物交配| 在线观看66精品国产| 长腿黑丝高跟| av黄色大香蕉| 嫩草影视91久久| 成人性生交大片免费视频hd| 日韩国内少妇激情av| 亚洲性夜色夜夜综合| 又爽又黄a免费视频| avwww免费| 一进一出抽搐动态| 免费在线观看影片大全网站| 99国产综合亚洲精品| 久久久久久久久中文| 少妇人妻一区二区三区视频| 久久6这里有精品| 欧美乱色亚洲激情| 欧美午夜高清在线| 欧美极品一区二区三区四区| 国产极品精品免费视频能看的| 亚洲精华国产精华精| 国产在线男女| 黄色视频,在线免费观看| 成年免费大片在线观看| 欧美一区二区精品小视频在线| 嫩草影院精品99| 波多野结衣巨乳人妻| 草草在线视频免费看| 午夜福利成人在线免费观看| av在线蜜桃| 在现免费观看毛片| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲性夜色夜夜综合| 99热这里只有是精品50| 高清日韩中文字幕在线| 国产亚洲欧美98| 美女大奶头视频| 国产亚洲av嫩草精品影院| 波野结衣二区三区在线| 国产精品爽爽va在线观看网站| 国内精品久久久久精免费| 日本五十路高清| 国产一区二区三区在线臀色熟女| 99热这里只有是精品50| 成人三级黄色视频| 一级黄片播放器| 国产乱人伦免费视频| 免费搜索国产男女视频| 男女那种视频在线观看| 亚洲综合色惰| 亚洲一区高清亚洲精品| 99国产极品粉嫩在线观看| a在线观看视频网站| 成人特级黄色片久久久久久久| 最后的刺客免费高清国语| 国产视频一区二区在线看| 亚洲成av人片免费观看| 亚洲第一区二区三区不卡| 搡老妇女老女人老熟妇| 久久久久国产精品人妻aⅴ院| 欧美激情在线99| 国产精品不卡视频一区二区 |