• <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)計關系
    同步可控震源地震采集技術新進展
    国产视频一区二区在线看| 中出人妻视频一区二区| 久久九九热精品免费| 亚洲七黄色美女视频| 亚洲第一青青草原| 国产精品久久久av美女十八| 国产成人欧美| 天堂俺去俺来也www色官网| 中文字幕精品免费在线观看视频| 一级毛片精品| 黄色怎么调成土黄色| √禁漫天堂资源中文www| 精品一区二区三区四区五区乱码| 亚洲熟女精品中文字幕| 国产精品 国内视频| 久久人妻福利社区极品人妻图片| 黑人操中国人逼视频| 久久这里只有精品19| 日本一区二区免费在线视频| 欧美精品人与动牲交sv欧美| 俄罗斯特黄特色一大片| 国产av精品麻豆| 久久精品人人爽人人爽视色| 天堂√8在线中文| 啦啦啦 在线观看视频| 国产免费现黄频在线看| 欧美日本中文国产一区发布| 男女床上黄色一级片免费看| 99国产精品免费福利视频| 国产不卡av网站在线观看| 村上凉子中文字幕在线| 国产高清视频在线播放一区| 老汉色∧v一级毛片| 亚洲精品在线美女| 9热在线视频观看99| 黑人操中国人逼视频| 色播在线永久视频| 免费在线观看日本一区| 久久青草综合色| 夜夜躁狠狠躁天天躁| 国产高清videossex| 成人18禁在线播放| 国产精品欧美亚洲77777| 久久久国产成人精品二区 | 黄频高清免费视频| 黄色视频,在线免费观看| 女人久久www免费人成看片| 黄色视频不卡| 国产免费现黄频在线看| 成人免费观看视频高清| 国产一区二区激情短视频| 天堂√8在线中文| 亚洲av美国av| 国产淫语在线视频| 人妻 亚洲 视频| 国产一卡二卡三卡精品| 久久精品国产清高在天天线| 91麻豆精品激情在线观看国产 | 欧美+亚洲+日韩+国产| 嫩草影视91久久| 伊人久久大香线蕉亚洲五| 12—13女人毛片做爰片一| 欧美日韩亚洲高清精品| 亚洲精品国产区一区二| 日韩免费av在线播放| 好男人电影高清在线观看| 国产无遮挡羞羞视频在线观看| 一二三四在线观看免费中文在| 热re99久久精品国产66热6| 亚洲成人免费电影在线观看| 少妇 在线观看| 国产欧美日韩一区二区三区在线| 12—13女人毛片做爰片一| 丝袜在线中文字幕| 成人av一区二区三区在线看| 三上悠亚av全集在线观看| 18禁观看日本| av线在线观看网站| tube8黄色片| 欧美精品高潮呻吟av久久| 国产真人三级小视频在线观看| 精品电影一区二区在线| 新久久久久国产一级毛片| 村上凉子中文字幕在线| 自线自在国产av| 悠悠久久av| 飞空精品影院首页| av天堂在线播放| 日本欧美视频一区| 色在线成人网| 精品亚洲成国产av| 99国产极品粉嫩在线观看| 老熟妇仑乱视频hdxx| 12—13女人毛片做爰片一| 99国产极品粉嫩在线观看| 村上凉子中文字幕在线| 免费少妇av软件| 免费看a级黄色片| 欧美在线黄色| 啦啦啦免费观看视频1| 久久精品熟女亚洲av麻豆精品| 欧美黑人精品巨大| 久久中文看片网| 午夜福利,免费看| 黄色视频不卡| 露出奶头的视频| 在线看a的网站| a级毛片在线看网站| 免费在线观看黄色视频的| 国产亚洲欧美98| 亚洲精品乱久久久久久| 男女高潮啪啪啪动态图| 国产成人av激情在线播放| 人人澡人人妻人| av一本久久久久| xxx96com| 丁香六月欧美| 韩国av一区二区三区四区| 成年人午夜在线观看视频| 久久国产精品影院| 好男人电影高清在线观看| 一a级毛片在线观看| 亚洲男人天堂网一区| 午夜亚洲福利在线播放| 亚洲精华国产精华精| 欧美色视频一区免费| 成年女人毛片免费观看观看9 | 午夜福利影视在线免费观看| 午夜激情av网站| 久久久久久人人人人人| www.自偷自拍.com| 国产在线观看jvid| a级片在线免费高清观看视频| 国产一区二区三区综合在线观看| 精品国产乱码久久久久久男人| 欧美日韩成人在线一区二区| 性少妇av在线| 国内久久婷婷六月综合欲色啪| 美女福利国产在线| 超碰成人久久| 一二三四在线观看免费中文在| 午夜精品国产一区二区电影| 村上凉子中文字幕在线| 丰满的人妻完整版| 老汉色av国产亚洲站长工具| 99国产综合亚洲精品| 国产精品 国内视频| 国产欧美日韩一区二区精品| av片东京热男人的天堂| 亚洲avbb在线观看| xxxhd国产人妻xxx| 99国产精品一区二区三区| 欧美成人午夜精品| 男人舔女人的私密视频| 亚洲美女黄片视频| 亚洲精品中文字幕在线视频| 国产男女内射视频| 久久亚洲精品不卡| 一级片免费观看大全| 欧美久久黑人一区二区| а√天堂www在线а√下载 | 国产日韩一区二区三区精品不卡| 欧美最黄视频在线播放免费 | 国产成人免费观看mmmm| 欧美中文综合在线视频| 精品久久久精品久久久| 水蜜桃什么品种好| 丝瓜视频免费看黄片| 巨乳人妻的诱惑在线观看| 天堂√8在线中文| 51午夜福利影视在线观看| 大香蕉久久成人网| 欧美一级毛片孕妇| 日日夜夜操网爽| 97人妻天天添夜夜摸| 三上悠亚av全集在线观看| 搡老岳熟女国产| 亚洲av美国av| 亚洲av日韩在线播放| 无人区码免费观看不卡| www.精华液| 精品国内亚洲2022精品成人 | 超碰97精品在线观看| 亚洲国产看品久久| 人人妻人人澡人人爽人人夜夜| 久久人妻熟女aⅴ| 一边摸一边抽搐一进一小说 | 国产精品久久久久久精品古装| 黑人巨大精品欧美一区二区蜜桃| 无遮挡黄片免费观看| 99re6热这里在线精品视频| 国产一区二区三区在线臀色熟女 | 黄片小视频在线播放| 亚洲,欧美精品.| 久久亚洲精品不卡| 99在线人妻在线中文字幕 | av不卡在线播放| 精品一区二区三区视频在线观看免费 | 欧美午夜高清在线| 视频区欧美日本亚洲| 欧美亚洲 丝袜 人妻 在线| 国产成人欧美在线观看 | 久久久精品免费免费高清| 91成年电影在线观看| 日韩大码丰满熟妇| 欧美黄色淫秽网站| 日韩欧美一区二区三区在线观看 | 欧美日韩亚洲综合一区二区三区_| 久久久国产精品麻豆| 啦啦啦 在线观看视频| 国产av又大| 99热网站在线观看| 免费观看人在逋| 18禁黄网站禁片午夜丰满| 人妻丰满熟妇av一区二区三区 | 国产精品98久久久久久宅男小说| 波多野结衣av一区二区av| 成人特级黄色片久久久久久久| 一区二区三区国产精品乱码| 国产无遮挡羞羞视频在线观看| 在线永久观看黄色视频| 国产精品av久久久久免费| 久久九九热精品免费| 国产成人欧美在线观看 | 在线视频色国产色| 日韩熟女老妇一区二区性免费视频| 又大又爽又粗| 日韩欧美在线二视频 | 成人亚洲精品一区在线观看| 99香蕉大伊视频| 黄片播放在线免费| 下体分泌物呈黄色| 国产成人系列免费观看| 成在线人永久免费视频| 一进一出抽搐动态| 黑人欧美特级aaaaaa片| 老司机在亚洲福利影院| 久久国产精品男人的天堂亚洲| 久久人妻熟女aⅴ| 久久久久国产一级毛片高清牌| 美国免费a级毛片| 久久精品国产a三级三级三级| 国产精品.久久久| 免费观看a级毛片全部| 久久国产精品人妻蜜桃| 国产精品1区2区在线观看. | 精品福利永久在线观看| 大码成人一级视频| 国产免费男女视频| 嫩草影视91久久| 亚洲熟妇中文字幕五十中出 | 午夜精品久久久久久毛片777| 中文字幕最新亚洲高清| 一级黄色大片毛片| 久久精品熟女亚洲av麻豆精品| 一级片免费观看大全| 亚洲精品国产色婷婷电影| 午夜福利影视在线免费观看| 成熟少妇高潮喷水视频| 亚洲av片天天在线观看| 又黄又爽又免费观看的视频| 黑人巨大精品欧美一区二区蜜桃| 久久国产精品大桥未久av| 国产成人一区二区三区免费视频网站| 一级毛片精品| 色在线成人网| 99re6热这里在线精品视频| 老熟妇仑乱视频hdxx| 色94色欧美一区二区| 十分钟在线观看高清视频www| 黄色怎么调成土黄色| 亚洲熟女精品中文字幕| 色尼玛亚洲综合影院| 精品福利观看| 久久久国产精品麻豆| 久久久久久亚洲精品国产蜜桃av| 国产乱人伦免费视频| 午夜福利,免费看| 中文字幕人妻熟女乱码| 久久 成人 亚洲| 欧美亚洲日本最大视频资源| 老汉色av国产亚洲站长工具| 欧美色视频一区免费| 一级毛片女人18水好多| 欧美 日韩 精品 国产| 啦啦啦 在线观看视频| 久久久久久亚洲精品国产蜜桃av| 久久青草综合色| 亚洲专区中文字幕在线| 狠狠婷婷综合久久久久久88av| 国产精品免费大片| 香蕉久久夜色| 王馨瑶露胸无遮挡在线观看| 国产av一区二区精品久久| 日本一区二区免费在线视频| 国产精品99久久99久久久不卡| 亚洲欧美色中文字幕在线| 久久精品国产a三级三级三级| 国产精品二区激情视频| 精品一区二区三区四区五区乱码| 久久国产精品影院| 亚洲熟妇中文字幕五十中出 | 免费日韩欧美在线观看| 亚洲精品一卡2卡三卡4卡5卡| 欧美久久黑人一区二区| 侵犯人妻中文字幕一二三四区| 日本黄色日本黄色录像| svipshipincom国产片| 十八禁高潮呻吟视频| 久久午夜亚洲精品久久| 亚洲午夜精品一区,二区,三区| 国产午夜精品久久久久久| 亚洲人成伊人成综合网2020| 一区二区三区国产精品乱码| 国产蜜桃级精品一区二区三区 | 真人做人爱边吃奶动态| 日本五十路高清| 中亚洲国语对白在线视频| 黄频高清免费视频| 亚洲美女黄片视频| 99国产精品免费福利视频| 自线自在国产av| 美女高潮喷水抽搐中文字幕| 亚洲一卡2卡3卡4卡5卡精品中文| 精品高清国产在线一区| 久久国产乱子伦精品免费另类| 欧美最黄视频在线播放免费 | 人妻久久中文字幕网| 水蜜桃什么品种好| 亚洲成人国产一区在线观看| 成人黄色视频免费在线看| 国产1区2区3区精品| 伊人久久大香线蕉亚洲五| 亚洲精品中文字幕一二三四区| 国产成人啪精品午夜网站| 欧美亚洲日本最大视频资源| 欧美激情极品国产一区二区三区| x7x7x7水蜜桃| av天堂久久9| 亚洲成人免费电影在线观看| 久久精品人人爽人人爽视色| 99热网站在线观看| 欧美日韩瑟瑟在线播放| 免费在线观看完整版高清| 色婷婷av一区二区三区视频| 国产精品一区二区精品视频观看| 成年版毛片免费区| 老熟女久久久| 狠狠婷婷综合久久久久久88av| 日韩欧美在线二视频 | 欧美 亚洲 国产 日韩一| 无限看片的www在线观看| 日日摸夜夜添夜夜添小说| 最近最新中文字幕大全电影3 | 法律面前人人平等表现在哪些方面| 在线观看免费视频日本深夜| 国产欧美日韩精品亚洲av| 国产成人一区二区三区免费视频网站| 久久精品国产综合久久久| 欧美黑人精品巨大| 麻豆国产av国片精品| 人人妻人人添人人爽欧美一区卜| 老熟妇乱子伦视频在线观看| 亚洲av美国av| 99热网站在线观看| 成人永久免费在线观看视频| 国产主播在线观看一区二区| 超色免费av| 国产99久久九九免费精品| 伊人久久大香线蕉亚洲五| 女性生殖器流出的白浆| 日本wwww免费看| 看片在线看免费视频| 极品少妇高潮喷水抽搐| 天天操日日干夜夜撸| 久久午夜综合久久蜜桃| 黄网站色视频无遮挡免费观看| 十八禁网站免费在线| 久久九九热精品免费| 黑人操中国人逼视频| 亚洲国产精品sss在线观看 | 亚洲国产欧美网| 亚洲美女黄片视频| 在线看a的网站| 婷婷丁香在线五月| 9热在线视频观看99| 欧美激情 高清一区二区三区| 午夜精品在线福利| 看片在线看免费视频| 人人妻,人人澡人人爽秒播| 久久中文字幕一级| 少妇的丰满在线观看| 亚洲一区中文字幕在线| 一进一出好大好爽视频| 制服诱惑二区| 精品高清国产在线一区| 国产成人系列免费观看| 男人的好看免费观看在线视频 | 老熟妇乱子伦视频在线观看| 精品国产美女av久久久久小说| 中文字幕人妻熟女乱码| 一级毛片精品| 欧美精品高潮呻吟av久久| 一级,二级,三级黄色视频| 免费在线观看日本一区| av中文乱码字幕在线| 亚洲精品中文字幕一二三四区| 欧洲精品卡2卡3卡4卡5卡区| 女人被躁到高潮嗷嗷叫费观| 欧美日韩亚洲综合一区二区三区_| 国产高清videossex| 女人被躁到高潮嗷嗷叫费观| 国产精品自产拍在线观看55亚洲 | 国产精品久久久久久精品古装| 中文字幕另类日韩欧美亚洲嫩草| av天堂久久9| 无遮挡黄片免费观看| 最新的欧美精品一区二区| 久久精品国产综合久久久| 国产精品国产av在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 又黄又爽又免费观看的视频| 欧美激情 高清一区二区三区| 久久草成人影院| 麻豆乱淫一区二区| 亚洲中文av在线| 法律面前人人平等表现在哪些方面| 久久婷婷成人综合色麻豆| 国产成人精品久久二区二区91| 免费观看精品视频网站| 国产成人免费无遮挡视频| 午夜福利欧美成人| 嫁个100分男人电影在线观看| 欧美乱码精品一区二区三区| 国产不卡一卡二| 色婷婷av一区二区三区视频| 建设人人有责人人尽责人人享有的| 一夜夜www| 99国产精品99久久久久| 国产欧美日韩一区二区三| 国产成人欧美在线观看 | 久久精品亚洲av国产电影网| av不卡在线播放| 欧美日韩中文字幕国产精品一区二区三区 | 美女扒开内裤让男人捅视频| 精品一区二区三区四区五区乱码| 一区二区三区国产精品乱码| 国产蜜桃级精品一区二区三区 | 美女视频免费永久观看网站| 亚洲精品中文字幕一二三四区| av国产精品久久久久影院| 日本欧美视频一区| 国产成人一区二区三区免费视频网站| 不卡av一区二区三区| 一个人免费在线观看的高清视频| 国产日韩一区二区三区精品不卡| 国产成人系列免费观看| 老司机午夜十八禁免费视频| 妹子高潮喷水视频| 久久亚洲真实| 国产精品乱码一区二三区的特点 | 色精品久久人妻99蜜桃| 在线国产一区二区在线| 一本一本久久a久久精品综合妖精| 久久草成人影院| 亚洲精品在线观看二区| 欧美日韩视频精品一区| 成年版毛片免费区| 国产有黄有色有爽视频| 精品少妇久久久久久888优播| 在线永久观看黄色视频| 精品欧美一区二区三区在线| av线在线观看网站| 香蕉国产在线看| 国产精品一区二区精品视频观看| 51午夜福利影视在线观看| 久久精品国产综合久久久| 啦啦啦 在线观看视频| 久久人人爽av亚洲精品天堂| 久热这里只有精品99| 大码成人一级视频| 人人妻,人人澡人人爽秒播| 免费在线观看黄色视频的| 999精品在线视频| 午夜精品在线福利| 一边摸一边做爽爽视频免费| 日本a在线网址| 女人爽到高潮嗷嗷叫在线视频| 色94色欧美一区二区| 丁香六月欧美| 天堂√8在线中文| 精品国产乱子伦一区二区三区| 一进一出好大好爽视频| 亚洲国产中文字幕在线视频| 丁香欧美五月| 99久久精品国产亚洲精品| 亚洲 国产 在线| 久9热在线精品视频| 又紧又爽又黄一区二区| 免费人成视频x8x8入口观看| 亚洲国产欧美网| videosex国产| 亚洲 欧美一区二区三区| 精品乱码久久久久久99久播| 国产精品免费大片| 国产一区有黄有色的免费视频| 精品无人区乱码1区二区| 国产精品综合久久久久久久免费 | 天天躁日日躁夜夜躁夜夜| 久久久精品免费免费高清| 国产亚洲欧美精品永久| 亚洲片人在线观看| 中文字幕色久视频| 国产日韩欧美亚洲二区| 精品无人区乱码1区二区| 18禁观看日本| 桃红色精品国产亚洲av| 黑人巨大精品欧美一区二区mp4| 欧美精品人与动牲交sv欧美| 老司机在亚洲福利影院| 男女下面插进去视频免费观看| av不卡在线播放| 俄罗斯特黄特色一大片| 韩国av一区二区三区四区| 91成年电影在线观看| 国产精品电影一区二区三区 | 久久国产精品人妻蜜桃| 午夜福利乱码中文字幕| 大片电影免费在线观看免费| 啪啪无遮挡十八禁网站| 法律面前人人平等表现在哪些方面| 巨乳人妻的诱惑在线观看| 中出人妻视频一区二区| tube8黄色片| 精品欧美一区二区三区在线| 国产成人系列免费观看| 最新在线观看一区二区三区| 精品久久久久久,| 久久久久国产一级毛片高清牌| 天天影视国产精品| 亚洲精品乱久久久久久| 国产成人系列免费观看| 精品午夜福利视频在线观看一区| 天天躁狠狠躁夜夜躁狠狠躁| 女人被狂操c到高潮| 午夜福利欧美成人| 午夜成年电影在线免费观看| 国产不卡一卡二| 国产精品免费大片| 久久精品aⅴ一区二区三区四区| 国产精品 欧美亚洲| 国产高清激情床上av| 法律面前人人平等表现在哪些方面| 精品国产一区二区三区四区第35| 欧美一级毛片孕妇| 久久久精品免费免费高清| 久久ye,这里只有精品| 欧美老熟妇乱子伦牲交| 黑人欧美特级aaaaaa片| 91麻豆精品激情在线观看国产 | 国产成人啪精品午夜网站| 久久天堂一区二区三区四区| xxx96com| 丰满饥渴人妻一区二区三| 最近最新免费中文字幕在线| 国产精品一区二区精品视频观看| 亚洲精品自拍成人| 久热爱精品视频在线9| 这个男人来自地球电影免费观看| 麻豆乱淫一区二区| 大香蕉久久成人网| 久久婷婷成人综合色麻豆| 视频区图区小说| 极品少妇高潮喷水抽搐| 中文字幕色久视频| 一a级毛片在线观看| 多毛熟女@视频| 国产免费av片在线观看野外av| 国产高清国产精品国产三级| 亚洲久久久国产精品| 啦啦啦在线免费观看视频4| 色94色欧美一区二区| 欧美乱色亚洲激情| 亚洲av日韩在线播放| 黄色成人免费大全| 国产精品成人在线| 欧美日韩亚洲高清精品| bbb黄色大片| 亚洲少妇的诱惑av| 欧美一级毛片孕妇| 久热这里只有精品99| 亚洲在线自拍视频| 日韩熟女老妇一区二区性免费视频| 王馨瑶露胸无遮挡在线观看| 极品教师在线免费播放| 国产精品乱码一区二三区的特点 | 建设人人有责人人尽责人人享有的| 亚洲精品自拍成人| 一区二区三区激情视频| 成人手机av| 嫩草影视91久久| 欧美黄色淫秽网站| 欧美在线黄色| 男女高潮啪啪啪动态图| 日本精品一区二区三区蜜桃| 啦啦啦视频在线资源免费观看| 俄罗斯特黄特色一大片| 成年人午夜在线观看视频| av超薄肉色丝袜交足视频| 激情视频va一区二区三区| 久9热在线精品视频| 日本五十路高清| 精品午夜福利视频在线观看一区| 另类亚洲欧美激情| 黄色怎么调成土黄色|