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

    考慮活動(dòng)地塊邊界帶的分區(qū)震級(jí)經(jīng)驗(yàn)關(guān)系

    2016-12-02 05:23:50李小軍梁建宏
    地震學(xué)報(bào) 2016年2期
    關(guān)鍵詞:震級(jí)關(guān)系式分區(qū)

    李 敏 李小軍 梁建宏

    1) 中國(guó)北京100045中國(guó)地震臺(tái)網(wǎng)中心 2) 中國(guó)北京100081中國(guó)地震局地球物理研究所

    ?

    考慮活動(dòng)地塊邊界帶的分區(qū)震級(jí)經(jīng)驗(yàn)關(guān)系

    李 敏1),李小軍2)梁建宏1)

    1) 中國(guó)北京100045中國(guó)地震臺(tái)網(wǎng)中心 2) 中國(guó)北京100081中國(guó)地震局地球物理研究所

    地方性震級(jí)ML和面波震級(jí)MS是兩種常用的震級(jí)標(biāo)度, 二者經(jīng)常需要依據(jù)經(jīng)驗(yàn)公式進(jìn)行相互轉(zhuǎn)換. 基于歷史地震整體統(tǒng)計(jì)的經(jīng)驗(yàn)公式由于資料有限, 且沒有考慮區(qū)域地質(zhì)條件的差異, 在實(shí)際使用中所產(chǎn)生的誤差較大. 本文整理了40多年來同時(shí)包含ML和MS的1萬2621條歷史地震數(shù)據(jù), 通過回歸分析得到了震級(jí)ML與MS的經(jīng)驗(yàn)關(guān)系; 并進(jìn)一步按照震源所處的Ⅱ級(jí)活動(dòng)地塊邊界帶進(jìn)行了分區(qū)回歸分析, 得到了27個(gè)分區(qū)內(nèi)地震震級(jí)ML與MS的經(jīng)驗(yàn)關(guān)系. 結(jié)果顯示, 經(jīng)過分區(qū)轉(zhuǎn)換后的ML與MS具有更好的相關(guān)性, 且考慮活動(dòng)地塊邊界帶所獲取的地震震級(jí)能有效地降低ML與MS之間的轉(zhuǎn)換誤差.

    地方性震級(jí)ML面波震級(jí)MS震級(jí)轉(zhuǎn)換 活動(dòng)地塊邊界帶

    引言

    地方性震級(jí)ML和面波震級(jí)MS均為表示地震規(guī)模大小的常規(guī)標(biāo)度.ML通常由短周期地震儀采集的1000 km區(qū)域內(nèi)的近震震相來測(cè)定, 計(jì)算公式為

    式中:AN和AE分別為NS向、 EW向S波或Lg波的最大振幅, 單位為μm;Δ為震中距, 單位為km;R(Δ)為地方性震級(jí)的量規(guī)函數(shù).MS通常由中、 長(zhǎng)周期地震儀采集的面波質(zhì)點(diǎn)最大位移確定, 是震級(jí)國(guó)家標(biāo)準(zhǔn)中規(guī)定的地震震級(jí), 其計(jì)算公式(陳運(yùn)泰, 劉瑞豐, 2004)為

    式中:A為兩水平分向地動(dòng)位移的矢量和, 單位為μm;Δ為震中距, 單位為度;T為A所對(duì)應(yīng)的周期, 單位為s.

    不同震級(jí)標(biāo)度反映了地震波在不同周期范圍內(nèi)所輻射的地震能量大小, 我國(guó)實(shí)際地震速報(bào)采用MS作為標(biāo)準(zhǔn)震級(jí). 在地震臺(tái)網(wǎng)的日常工作中, 雖然所測(cè)定的不同類型震級(jí)之間一律不進(jìn)行轉(zhuǎn)換, 但是在實(shí)際工作中, 又需要各類震級(jí)之間的關(guān)系, 通常使用經(jīng)驗(yàn)公式將一類震級(jí)換算成另一類震級(jí). 郭履燦(1971)通過整理邢臺(tái)地震的余震資料, 提出經(jīng)驗(yàn)關(guān)系MS=1.13ML-1.08, 這一公式多年來一直被普遍采用; 甚至在目前的地震速報(bào)與編目中, 區(qū)域臺(tái)網(wǎng)使用的也是這一換算關(guān)系. 近年來張宏志等(2008)和汪素云等(2010b)利用1990—2008年的地震資料對(duì)這一經(jīng)驗(yàn)關(guān)系進(jìn)行改進(jìn), 但均未涉及不同區(qū)域地殼結(jié)構(gòu)特性對(duì)地震波傳播衰減特性的影響. 由于區(qū)域地質(zhì)環(huán)境對(duì)ML測(cè)定的影響很大(汪素云等, 2007), 適用于某個(gè)地區(qū)的ML與MS經(jīng)驗(yàn)關(guān)系并不一定適用于其它地區(qū), 也就是說, 區(qū)域地質(zhì)環(huán)境差異對(duì)建立震級(jí)轉(zhuǎn)換關(guān)系的影響不容忽視. 近年來, 隨著數(shù)字化地震觀測(cè)技術(shù)的發(fā)展, 大量地震資料中同時(shí)給出了ML和MS的地震記錄, 使得更系統(tǒng)地研究ML與MS之間的關(guān)系成為可能.

    本文首先根據(jù)活動(dòng)地塊將40多年來數(shù)十萬條地震目錄按區(qū)域進(jìn)行劃分, 分析現(xiàn)有方法所產(chǎn)生的ML與MS之間的轉(zhuǎn)換偏差, 并進(jìn)一步討論該誤差與區(qū)域地質(zhì)特征的相關(guān)性; 然后參照已有的關(guān)于ML與MS關(guān)系的研究成果(汪素云等, 2010a), 嘗試按照不同區(qū)域特征建立新的經(jīng)驗(yàn)關(guān)系, 為區(qū)域臺(tái)網(wǎng)和國(guó)家臺(tái)網(wǎng)的速報(bào)、 編目以及地震活動(dòng)性研究提供依據(jù).

    1 資料

    本文所選用資料引自中國(guó)地震臺(tái)網(wǎng)中心(2015)的地震目錄, 共包含1970年1月—2012年2月約14.4萬條地震數(shù)據(jù), 其中7.2萬條有MS記錄, 3.3萬條有ML記錄; 而同時(shí)具有MS和ML震級(jí)記錄的數(shù)據(jù)僅有1萬2621條, 占全部記錄的8.77%; 超過一半的歷史地震僅有MS或ML記錄. 圖1給出了各類震級(jí)數(shù)據(jù)所占的比例. 圖2給出了本文所選取地震目錄中同時(shí)具有MS和ML記錄的震級(jí)分布, 可以看出低震級(jí)(M<4.5)地震的MS記錄較多.

    圖1 震級(jí)數(shù)據(jù)組成

    圖2 同時(shí)具有MS和ML記錄的地震震級(jí)分布

    2 MS與ML經(jīng)驗(yàn)關(guān)系

    分析萬余條針對(duì)同一地震獨(dú)立獲得的MS和ML震級(jí)數(shù)據(jù)可知, 二者之間存在良好的線性關(guān)系, 可用

    (1)

    表達(dá). 通過線性回歸方法歸納獲取二者的經(jīng)驗(yàn)關(guān)系式是合理的做法, 郭履燦(1971)很早就總結(jié)出了MS與ML的經(jīng)驗(yàn)關(guān)系式:

    (2)

    該關(guān)系式在很多年內(nèi)均得到廣泛使用. 然而, 由于所使用的地震記錄有限, 該關(guān)系式與近40年來的觀測(cè)結(jié)果之間已存在明顯偏差(圖3). 近年來, 汪素云等(2010b)基于1990—2008年長(zhǎng)達(dá)19年的地震數(shù)據(jù)采用雙隨機(jī)變量法(Madansky, 1959; Wyss, Habermann, 1982)將該關(guān)系式修正為

    圖3 MS與ML震級(jí)的經(jīng)驗(yàn)關(guān)系 Fig.3 Empirical relationships between MS and ML The orange and the green lines represent the empirical equations from Eqs.(2) and (3), respectively. The blue and the black lines represents the linear regres-sion relations from Eqs.(4) and (5), respectively

    (3)

    本文則基于整理所得的1萬2621條地震記錄, 采用最小二乘法進(jìn)行一元線性回歸得到ML與MS的關(guān)系式為

    (4)

    基于同樣的數(shù)據(jù)集, 采用與汪素云等(2010b)相同的方法得到的ML與MS的關(guān)系式為

    (5)

    將本文選取的MS和ML數(shù)據(jù)繪于圖3, 并將由式(2)—(5)得到的回歸直線繪于圖中, 以觀察各直線與觀測(cè)點(diǎn)之間的擬合情況. 可以看到, 式(3)和式(5)

    圖4 基于MS和ML共有記錄得到的4組線性經(jīng)驗(yàn)關(guān)系的誤差分布Fig.4 The error distribution of the four linear empirical relationships based on the records with both MS and ML

    由于采用雙隨機(jī)變量法來估計(jì)MS與ML的經(jīng)驗(yàn)關(guān)系式, 因此能較好地反映中小地震MS與ML的經(jīng)驗(yàn)關(guān)系. 二者的細(xì)微差別源自所用數(shù)據(jù)集的更新, 這說明利用雙隨機(jī)變量法來統(tǒng)計(jì)MS與ML的線性轉(zhuǎn)換關(guān)系較為合理; 然而, 當(dāng)ML>6.0時(shí), 二者均傾向于低估對(duì)應(yīng)的MS. 從圖3中數(shù)據(jù)本身的離散性也可以看出, 無論采用哪個(gè)關(guān)系式, 均存在大量偏離真實(shí)的點(diǎn), 偏離程度最大可達(dá)1.0以上(圖4). 從圖4中基于式(2)—(5)的4組MS與ML線性關(guān)系的誤差分布可以看出, 式(3)—(5)的計(jì)算誤差分布相差不大.

    地方性震級(jí)ML的測(cè)量取決于區(qū)域臺(tái)網(wǎng)所獲取的地震最大振幅,ML的大小會(huì)受到局部地震傳播介質(zhì)的影響(張國(guó)民等, 2004; 汪素云等, 2007, 2008), 因此對(duì)數(shù)據(jù)進(jìn)行整體線性回歸而不考慮震源地的實(shí)際情況, 必然會(huì)導(dǎo)致MS與ML的相關(guān)性減弱. 為了更精確地確定MS與ML的經(jīng)驗(yàn)關(guān)系, 有必要將不同地質(zhì)條件對(duì)ML測(cè)定的影響考慮在內(nèi). 由于地震多集中于活動(dòng)地塊之間的邊界帶區(qū)域, 兩側(cè)活動(dòng)地塊相同的地震震源所處地可能具有更相似的地質(zhì)環(huán)境, 因此本文考慮按活動(dòng)地塊劃分區(qū)域來研究MS與ML的相關(guān)性.

    3 活動(dòng)地塊及其邊界帶

    活動(dòng)地塊是具有相對(duì)統(tǒng)一運(yùn)動(dòng)方式的地質(zhì)單元(張培震等, 2003). 我國(guó)大陸包含青藏、 西域、 南華、 滇緬、 華北和東北亞等6個(gè)Ⅰ級(jí)活動(dòng)地塊區(qū), Ⅰ級(jí)活動(dòng)地塊又可分成22個(gè)Ⅱ級(jí)活動(dòng)地塊區(qū), 詳見表1(張國(guó)民等, 2005).

    本文依照張國(guó)民等(2005)劃定的我國(guó)大陸26條活動(dòng)地塊邊界帶和臺(tái)灣地區(qū)所在的板塊邊界, 來分區(qū)研究MS與ML之間的經(jīng)驗(yàn)關(guān)系. 如圖5所示, 將本文所收集的1萬2621條地震記錄按其經(jīng)緯度繪于圖中, 可見我國(guó)歷史上大多數(shù)地震均發(fā)生在活動(dòng)地塊之間的邊界帶上. 根據(jù)活動(dòng)地塊邊界帶的分布特點(diǎn), 并考慮本文所提出的MS與ML轉(zhuǎn)換關(guān)系式的實(shí)用性, 將邊界帶近似處理成線段; 然后根據(jù)震中位置, 將地震劃分至不同的邊界帶分區(qū)上. 由于活動(dòng)地塊邊界帶已近似為線段, 為便于應(yīng)用, 按照張國(guó)民等(2005)提供的邊界帶經(jīng)緯度信息, 本文將邊界帶處理為由起始經(jīng)緯度確定的、 寬約3°的狹窄矩形區(qū)域, 具體的劃分標(biāo)準(zhǔn)詳見表1. 對(duì)于待轉(zhuǎn)換地震, 只需已知震中的經(jīng)緯度信息, 即可得到其所屬邊界帶的相關(guān)信息.

    按照表1所示的邊界帶和分區(qū)標(biāo)準(zhǔn), 將同時(shí)具有MS和ML記錄的地震按其所處邊界帶分區(qū), 落入邊界帶兩側(cè)1.5°范圍內(nèi)的地震視為兩側(cè)活動(dòng)地塊相互作用所產(chǎn)生, 并根據(jù)該范圍內(nèi)地震記錄的MS和ML數(shù)據(jù)分別統(tǒng)計(jì)適用于該邊界帶分區(qū)的關(guān)系式. 在總體經(jīng)驗(yàn)關(guān)系式的推導(dǎo)中, 本文依然采用誤差較小的雙隨機(jī)變量法分別計(jì)算不同邊界帶分區(qū)的經(jīng)驗(yàn)關(guān)系式, 結(jié)果列于表2. 可以看到, 所得分區(qū)關(guān)系式的斜率分布較為分散, 從北天山的0.68到岷山—龍門山的0.99不等, 說明MS與ML的經(jīng)驗(yàn)關(guān)系與所處地質(zhì)條件的相關(guān)性較大.

    圖5 活動(dòng)地塊及其邊界帶

    圖6 基于鮮水河邊界帶和北天山邊界帶地震數(shù)據(jù) 分別擬合得到的MS與ML經(jīng)驗(yàn)關(guān)系曲線 Fig.6 Empirical relationships fitting between MSand ML based on the earthquake data within Xian-shuihe and northern Tianshan boundary areas The red dots represent earthquake data from Xian-shuihe boundary area, the blue dots represent data from northern Tianshan boundary area, and the grey dots represent data from other areas

    此外, 局部地質(zhì)條件對(duì)MS與ML經(jīng)驗(yàn)關(guān)系的影響還體現(xiàn)在考慮分區(qū)前后MS與ML的相關(guān)性變化上. 未考慮分區(qū)時(shí),MS與ML的總體皮爾森(Pearson)相關(guān)系數(shù)為0.834; 將分區(qū)考慮在內(nèi)分別計(jì)算27個(gè)邊界帶分區(qū)數(shù)據(jù)的MS與ML相關(guān)系數(shù), 其中20個(gè)分區(qū)的相關(guān)系數(shù)均高于總體相關(guān)系數(shù). 對(duì)各邊界帶分區(qū)的相關(guān)系數(shù)與總體相關(guān)系數(shù)的差值(提升值)求和, 得到相關(guān)系數(shù)總的提升值高達(dá)0.828. 這充分說明分組以后的MS和ML數(shù)據(jù)更加收斂, 按活動(dòng)地塊邊界帶分區(qū)分析MS與ML經(jīng)驗(yàn)關(guān)系能在很大程度上排除地質(zhì)條件差異所帶來的擾動(dòng). 圖6給出了基于北天山和鮮水河兩條邊界帶的MS和ML數(shù)據(jù)擬合得到的經(jīng)驗(yàn)關(guān)系. 可以看到, 由兩條邊界帶分區(qū)的地震數(shù)據(jù)得到的相應(yīng)關(guān)系曲線有明顯差異; 若利用基于其它邊界帶的背景數(shù)據(jù)(灰色圓點(diǎn))所得關(guān)系式分析北天山邊界帶的地震數(shù)據(jù), 則會(huì)顯示多數(shù)MS數(shù)值被高估.

    表2 不同活動(dòng)地塊邊界帶分區(qū)地震MS與ML的經(jīng)驗(yàn)關(guān)系

    注:r表示各邊界帶分區(qū)內(nèi)地震震級(jí)MS與ML的皮爾森相關(guān)系數(shù); Δr表示由邊界帶分區(qū)內(nèi)地震統(tǒng)計(jì)得到的相關(guān)系數(shù)與總體相關(guān)系數(shù)的差值.

    4 震級(jí)分區(qū)轉(zhuǎn)換與誤差分析

    4.1 震級(jí)分區(qū)

    根據(jù)表1和表2, 可對(duì)國(guó)內(nèi)任一地震進(jìn)行分區(qū)轉(zhuǎn)換, 具體步驟為: ① 由地震震中位置, 根據(jù)表1所示經(jīng)緯度范圍查找到相應(yīng)的邊界帶, 并通過計(jì)算與其經(jīng)度對(duì)應(yīng)的緯度范圍確定其是否符合相應(yīng)的分區(qū)標(biāo)準(zhǔn), 最多允許1.5°差距; ② 若該地震符合多個(gè)邊界帶分區(qū)標(biāo)準(zhǔn), 則選擇實(shí)際緯度最接近緯度范圍中值的邊界帶分區(qū)所對(duì)應(yīng)的轉(zhuǎn)換關(guān)系式; ③ 若該地震無法劃分至任一邊界帶分區(qū)(約占總地震35%), 則按不分區(qū)的關(guān)系式即式(5)進(jìn)行轉(zhuǎn)換.

    為便于對(duì)歷史地震進(jìn)行大規(guī)模分析研究, 本文研究開發(fā)了震級(jí)轉(zhuǎn)換程序(詳見https:∥git.oschina.net/littlemin/ML-MS_caculator), 可快速實(shí)現(xiàn)MS與ML的相互轉(zhuǎn)換.

    4.2 誤差分析

    圖7 按邊界帶分區(qū)計(jì)算結(jié)果與式(3)計(jì)算結(jié)果的誤差分布對(duì)比Fig.7 Comparison of the conversion error of equations taking active tectonic boun-daries into consideration (red) with that by using Eq.(3)

    由于式(3)—(5)的計(jì)算誤差分布接近, 本文選取汪素云等(2010b)提出的式(3)作為MS與ML整體分布關(guān)系的代表, 將其與考慮活動(dòng)地塊邊界帶的計(jì)算結(jié)果進(jìn)行比較, 結(jié)果如圖7所示. 可以看出, 考慮邊界帶分區(qū)所得計(jì)算結(jié)果的估計(jì)誤差更?。?在參與比較的1萬1712條地震記錄中, 邊界帶分區(qū)計(jì)算的估計(jì)誤差在-0.1-0.1范圍內(nèi)的記錄共計(jì)4653條, 占總記錄的40%; 基于總體數(shù)據(jù)通過式(3)計(jì)算震級(jí), 其誤差為-0.1-0.1范圍內(nèi)的記錄有4315條, 占總記錄的37%. 另一方面, 式(3)估計(jì)結(jié)果中誤差>0.5的地震記錄為1147條, 而邊界帶分區(qū)估計(jì)的結(jié)果中誤差>0.5的地震記錄為1048條. 由此可見, 邊界帶分區(qū)估計(jì)結(jié)果中誤差≤0.1的地震次數(shù)較總體估計(jì)所得約多10%, 同時(shí)誤差>0.5的錯(cuò)誤估計(jì)數(shù)比總體估計(jì)所得約少10%, 故考慮邊界帶分區(qū)的MS與ML經(jīng)驗(yàn)關(guān)系式的精度更高.

    5 討論與結(jié)論

    1935年, 美國(guó)加州理工學(xué)院的里克特最早提出使用ML來表征美國(guó)南加州地震的大小(Richter, 1935), 此后ML在全世界范圍內(nèi)被廣泛使用. 然而, 由于局部地質(zhì)環(huán)境差異較大, 不同地區(qū)的ML震級(jí)難以直接比較, 因此需要將ML轉(zhuǎn)換成其它震級(jí)標(biāo)度, 例如, 體波震級(jí)Mb, 矩震級(jí)MW(Ristau, 2003, 2009; Castelloetal, 2007; Yenieretal, 2008)或面波震級(jí)MS(汪素云等, 2010b).

    對(duì)同時(shí)具有多個(gè)震級(jí)標(biāo)度的地震記錄進(jìn)行回歸分析是推導(dǎo)不同震級(jí)標(biāo)度經(jīng)驗(yàn)關(guān)系的常規(guī)手段(Bindi, 2005; Deniz, Yucemen, 2010; Scordilis, 2006). 以往的研究通常僅對(duì)全國(guó)地震記錄進(jìn)行總體統(tǒng)計(jì), 而忽視了諸如地殼速度、 是否存在俯沖板塊等局部地質(zhì)環(huán)境差異對(duì)ML計(jì)算的影響(Miao, Langston, 2007). 近年來已有研究注意到該問題, 例如, Ristau(2003)和Ristau等(2005)將地震波傳播介質(zhì)的差異考慮在內(nèi), 分別統(tǒng)計(jì)了加拿大西海岸線和西部陸地ML與MW的經(jīng)驗(yàn)關(guān)系.

    本文按照地震震源所處的活動(dòng)地塊邊界帶, 研究了近40余年來我國(guó)歷史地震記錄中MS與ML震級(jí)的對(duì)應(yīng)關(guān)系, 給出了我國(guó)總體的MS與ML震級(jí)經(jīng)驗(yàn)關(guān)系式和27個(gè)活動(dòng)地塊邊界帶分區(qū)的MS與ML震級(jí)經(jīng)驗(yàn)關(guān)系式. 不同MS與ML震級(jí)轉(zhuǎn)換經(jīng)驗(yàn)關(guān)系式的相關(guān)性表明, 按活動(dòng)地塊邊界帶將地震分區(qū), 能夠顯著提高M(jìn)S與ML的相關(guān)性, 有助于實(shí)現(xiàn)MS與ML之間更準(zhǔn)確的轉(zhuǎn)換. 本文關(guān)于歷史地震的MS與ML經(jīng)驗(yàn)關(guān)系的比較分析結(jié)果表明, 分區(qū)震級(jí)轉(zhuǎn)換在降低轉(zhuǎn)換誤差的同時(shí), 也提高了轉(zhuǎn)換精度.

    由于缺乏精確的活動(dòng)地塊邊界帶數(shù)據(jù), 本文劃定震源地所處Ⅱ級(jí)活動(dòng)地塊邊界帶的方法還較為粗放. 此外, 由于不同邊界帶分區(qū)中可統(tǒng)計(jì)的地震次數(shù)差異較大, 其中6個(gè)分區(qū)內(nèi)的可統(tǒng)計(jì)地震次數(shù)小于100, 可能會(huì)存在統(tǒng)計(jì)偏差; 對(duì)于這些區(qū)域, 可使用總體轉(zhuǎn)換公式進(jìn)行轉(zhuǎn)換.

    無論是總體轉(zhuǎn)換關(guān)系式還是分區(qū)域轉(zhuǎn)換關(guān)系式, 在ML>6.0時(shí)均傾向于低估相應(yīng)的ML震級(jí)(圖3, 6). 這是由于低震級(jí)地震記錄較多, 使得經(jīng)驗(yàn)關(guān)系的統(tǒng)計(jì)受低震級(jí)記錄影響較大所致; 同時(shí)也說明MS與ML的關(guān)系可能并非呈完全線性.ML在強(qiáng)震發(fā)生時(shí)會(huì)趨于飽和, 其最大值不會(huì)超過7.0(Kanamori, 1983), 因此不宜使用線性回歸關(guān)系式進(jìn)行強(qiáng)震的震級(jí)轉(zhuǎn)換, 可考慮采用非線性方法.

    本研究給出了提高M(jìn)S與ML震級(jí)轉(zhuǎn)換經(jīng)驗(yàn)關(guān)系精度的一個(gè)技術(shù)途徑. 隨著地質(zhì)資料的不斷完善和地震數(shù)據(jù)的繼續(xù)擴(kuò)充, 可綜合考慮震中位置、 震源深度、 活動(dòng)地塊、 地殼厚度等信息, 建立更精細(xì)的轉(zhuǎn)換模型, 這不僅有利于震級(jí)轉(zhuǎn)換的實(shí)際應(yīng)用, 也為更深入地認(rèn)識(shí)地震傳播規(guī)律提供參考.

    陳運(yùn)泰, 劉瑞豐. 2004. 地震的震級(jí)[J]. 地震地磁觀測(cè)與研究, 25(6): 1--12.

    Chen Y T, Liu R F. 2004. Earthquake magnitude[J].SeismologicalandGeomagneticObservationandResearch, 25(6): 1--12 (in Chinese).

    郭履燦. 1971. 華北地區(qū)的地方性震級(jí)ML和面波震級(jí)MS經(jīng)驗(yàn)關(guān)系[C]∥全國(guó)地震工作會(huì)議: 1--10.

    Guo L C. 1971. Empirical relationship of local magnitudeMLand surface magnitudeMS[C]∥ConferenceRecordofNationalSeismologyWork: 1--10 (in Chinese).

    汪素云, 裴順平, 許忠淮, 張國(guó)民, 俞言祥, 潘華. 2007. 利用ML振幅研究地殼橫波Q值Ⅰ: 不同構(gòu)造區(qū)的衰減特征[J]. 地球物理學(xué)報(bào), 50(6): 1740--1747.

    Wang S Y, Pei S P, Xu Z H, Zhang G M, Yu Y X, Pan H. 2007. Crustal S-waveQestimated fromMLamplitude Ⅰ: Attenuation in different tectonic regions of China[J].ChineseJournalofGeophysics, 50(6): 1740--1747 (in Chinese).

    汪素云, 裴順平, Hearn T M, 許忠淮, Ni J F, 俞言祥. 2008. 利用ML振幅研究地殼橫波Q值Ⅱ:Q橫向變化特征[J]. 地球物理學(xué)報(bào), 51(1): 133--139.

    Wang S Y, Pei S P, Hearn T M, Xu Z H, Ni J F, Yu Y X. 2008. Crustal S-waveQestimated fromMLamplitude Ⅱ:Qlateral variation in China[J].ChineseJournalofGeophysics, 51(1): 133--139 (in Chinese).

    汪素云, 高阿甲, 馮義鈞, 和銳. 2010a. 中國(guó)地震目錄間的對(duì)比及標(biāo)準(zhǔn)化[J]. 地震, 30(2): 38--45.

    Wang S Y, Gao A J, Feng Y J, He R. 2010a. Comparison and standardization of the Chinese earthquake catalogs[J].Earthquake, 30(2): 38--45 (in Chinese).

    汪素云, 王健, 俞言祥, 吳清, 高阿甲, 高孟潭. 2010b. 基于中國(guó)地震臺(tái)網(wǎng)觀測(cè)報(bào)告的ML與MS經(jīng)驗(yàn)關(guān)系[J]. 中國(guó)地震, 26(1): 14--22.

    Wang S Y, Wang J, Yu Y X, Wu Q, Gao A J, Gao M T. 2010b. The empirical relation betweenMLandMSbased on bulletin of seismological observations of Chinese stations[J].EarthquakeResearchinChina, 26(1): 14--22 (in Chinese).

    張國(guó)民, 馬宏生, 王輝, 李麗. 2004. 中國(guó)大陸活動(dòng)地塊與強(qiáng)震活動(dòng)關(guān)系[J]. 中國(guó)科學(xué): D輯, 34(7): 591--599.

    Zhang G M, Ma H S, Wang H, Li L. 2004. Relation between active-tectonic blocks and strong earthquakes in the China mainland[J].ScienceinChina:SeriesD, 34(7): 591--599 (in Chinese).

    張國(guó)民, 馬宏生, 王輝, 王新嶺. 2005. 中國(guó)大陸活動(dòng)地塊邊界帶與強(qiáng)震活動(dòng)[J]. 地球物理學(xué)報(bào), 48(3): 602--610.

    Zhang G M, Ma H S, Wang H, Wang X L. 2005. Boundaries between active-tectonic blocks and strong earthquakes in the China mainland[J].ChineseJournalofGeophysics, 48(3): 602--610 (in Chinese).

    張宏志, 刁桂苓, 趙明淳, 王勤彩, 張驍, 黃媛. 2008. 不同標(biāo)度震級(jí)關(guān)系和臺(tái)基影響問題探討[J]. 中國(guó)地震, 23(2): 141--146.

    Zhang H Z, Diao G L, Zhao M C, Wang Q C, Zhang X, Huang Y. 2008. Discussion on relation between different earthquake magnitude scales and effect of seismic station site on magnitude estimation[J].EarthquakeResearchinChina, 23(2): 141--146 (in Chinese).

    張培震, 鄧起東, 張國(guó)民, 馬瑾, 甘衛(wèi)軍, 閔偉, 毛鳳英, 王琪. 2003. 中國(guó)大陸的強(qiáng)震活動(dòng)與活動(dòng)地塊[J]. 中國(guó)科學(xué): D輯, 33(增刊): 12--20.

    Zhang P Z, Deng Q D, Zhang G M, Ma J, Gan W J, Min W, Mao F Y, Wang Q. 2003. Active tectonic blocks and strong earthquakes in the continent of China[J].ScienceinChina:SeriesD, 46(S2): 13--24.

    中國(guó)地震臺(tái)網(wǎng).2015. 國(guó)家臺(tái)網(wǎng)大震速報(bào)目錄[EB/OL]. [2015-01-14]. http:∥www.csndmc.ac.cn/newweb/data.htm.

    China Seismic Networks Center. 2015. China seismic networks earthquake catalogue[EB/OL]. [2015-01-14]. http:∥www.csndmc.ac.cn/newweb/data.htm.

    Bindi D. 2005. Local and duration magnitudes in northwestern Italy, and seismic momentversusmagnitude relationships[J].BullSeismolSocAm, 95(2): 592--604.

    Castello B, Olivieri M, Selvaggi G. 2007. Local and duration magnitude determination for the Italian earthquake catalog, 1981—2002[J].BullSeismolSocAm, 97(1B): 128--139.

    Deniz A, Yucemen M S. 2010. Magnitude conversion problem for the Turkish earthquake data[J].NatHazards, 55(2): 333--352.

    Kanamori H. 1983. Magnitude scale and quantification of earthquakes[J].Tectonophysics, 93(3/4): 185--199.

    Madansky A. 1959. The fitting of straight lines when both variables are subject to error[J].JAmStatistAssoc, 54(285): 173--205.

    Miao Q, Langston C A. 2007. Empirical distance attenuation and the local-magnitude scale for the central United States[J].BullSeismolSocAm, 97(6): 2137--2151.

    Richter C F. 1935. An instrumental earthquake magnitude scale[J].BullSeismolSocAm, 25(1): 1--32.

    Ristau J. 2003. Moment magnitude-local magnitude calibration for earthquakes off Canada’s west coast[J].BullSeismolSocAm, 93(5): 2296--2300.

    Ristau J, Rogers G C, Cassidy J F. 2005. Moment magnitude-local magnitude calibration for earthquakes in western Canada[J].BullSeismolSocAm, 95(5): 1994--2000.

    Ristau J. 2009. Comparison of magnitude estimates for New Zealand earthquakes: Moment magnitude, local magnitude, and teleseismic body-wave magnitude[J].BullSeismolSocAm, 99(3): 1841--1852.

    Scordilis E M. 2006. Empirical global relations convertingMSandmbto moment magnitude[J].JSeismol, 10(2): 225--236.

    Wyss M, Habermann R E. 1982. Conversion ofmbtoMSfor estimating the recurrence time of large earthquakes[J].BullSeismolSocAm, 72(5): 1651--1662.

    姚道平, 張藝峰, 閆培, 李普春, 王筍, 郭曉然. 2016. 臺(tái)灣海峽大容量氣槍震源海陸聯(lián)測(cè)初探. 地震學(xué)報(bào), 38(2): 167--178. doi:10.11939/jass.2016.02.002.

    Empirical relationships of earthquake magnitudes for local regions considering the active tectonic boundary areas

    Li Min1),Li Xiaojun2)Liang Jianhong1)

    1)ChinaEarthquakeNetworksCenter,Beijing100045,China2)InstituteofGeophysics,ChinaEarthquakeAdministration,Beijing100081,China

    Local magnitudeMLand surface wave magnitudeMSare two widely used scales of earthquake magnitude, which usually need to be converted to each other using empirical equations. Due to a limit of data, and ignorance the difference in geological conditions, empirical equations based on historical earthquakes suffer from conversion errors in practice. Based on the historical earthquake data of 12 621 records with bothMLandMSrecorded in the past four decades, this paper derived the empirical conversion equations betweenMLandMSby regression analysis. Furthermore, empirical conversion equations were derived from the partitioned data for the 27 regions classified according to the secondary active tectonic boundary areas. The results show that the correlation betweenMSandMLis more significant after the magnitude data are partitioned, and the magnitudes that were converted with consideration of active tectonic boundaries showed lower conversion error.

    local magnitudeML; surface wave magnitudeMS; magnitude conversion; active tectonic boundary area

    李敏, 李小軍, 梁建宏. 2016. 考慮活動(dòng)地塊邊界帶的分區(qū)震級(jí)經(jīng)驗(yàn)關(guān)系. 地震學(xué)報(bào), 38(2): 157--166. doi:10.11939/jass.2016.02.001.

    Li M, Li X J, Liang J H. 2016. Empirical relationships of earthquake magnitudes for local regions considering the active tectonic boundary areas.ActaSeismologicaSinica, 38(2): 157--166. doi:10.11939/jass.2016.02.001.

    2015-09-07收到初稿, 2016-01-18決定采用修改稿.

    e-mail: limin@seis.ac.cn

    國(guó)家自然科學(xué)基金項(xiàng)目(U1434210, 51421005)、 國(guó)家國(guó)際科技合作項(xiàng)目(2012DFG20510)和北京市屬高等學(xué)校創(chuàng)新團(tuán)隊(duì)建設(shè)提升計(jì)劃共同資助.

    10.11939/jass.2016.02.001

    P315.3

    A

    猜你喜歡
    震級(jí)關(guān)系式分區(qū)
    基于累積絕對(duì)位移值的震級(jí)估算方法
    上海實(shí)施“分區(qū)封控”
    例談同角三角函數(shù)基本關(guān)系式的應(yīng)用
    地震后各國(guó)發(fā)布的震級(jí)可能不一樣?
    新震級(jí)國(guó)家標(biāo)準(zhǔn)在大同臺(tái)的應(yīng)用與評(píng)估
    山西地震(2020年1期)2020-04-08 07:34:26
    浪莎 分區(qū)而治
    速尋關(guān)系式巧解計(jì)算題
    明確關(guān)系式
    中國(guó)地震臺(tái)網(wǎng)面波震級(jí)與矩震級(jí)的統(tǒng)計(jì)關(guān)系
    基于SAGA聚類分析的無功電壓控制分區(qū)
    黄片大片在线免费观看| 此物有八面人人有两片| 亚洲国产欧洲综合997久久, | 亚洲国产毛片av蜜桃av| 久热这里只有精品99| 日韩免费av在线播放| 国产精品永久免费网站| 色综合欧美亚洲国产小说| 欧美久久黑人一区二区| 怎么达到女性高潮| 国产aⅴ精品一区二区三区波| 欧美一级a爱片免费观看看 | 人人妻人人澡人人看| 精品熟女少妇八av免费久了| 精品熟女少妇八av免费久了| 禁无遮挡网站| 这个男人来自地球电影免费观看| 久久久久国内视频| 国产成人影院久久av| 黄频高清免费视频| 成人午夜高清在线视频 | 色婷婷久久久亚洲欧美| 精品国产超薄肉色丝袜足j| 国产精品久久久久久精品电影 | 999久久久精品免费观看国产| 国产av一区二区精品久久| 国产伦人伦偷精品视频| 精品久久久久久成人av| 免费无遮挡裸体视频| 国产又爽黄色视频| 色精品久久人妻99蜜桃| 丰满人妻熟妇乱又伦精品不卡| 白带黄色成豆腐渣| 午夜免费观看网址| 99riav亚洲国产免费| svipshipincom国产片| 手机成人av网站| 精品国产一区二区三区四区第35| 久久精品91无色码中文字幕| 男女那种视频在线观看| 12—13女人毛片做爰片一| 日本精品一区二区三区蜜桃| 欧美日韩一级在线毛片| 亚洲一区二区三区不卡视频| 一a级毛片在线观看| 最近最新中文字幕大全免费视频| 国内久久婷婷六月综合欲色啪| 欧美不卡视频在线免费观看 | 国产成人精品无人区| 在线av久久热| 亚洲精品一卡2卡三卡4卡5卡| 久久狼人影院| 高清毛片免费观看视频网站| 日韩欧美国产在线观看| 国产成人欧美| 变态另类丝袜制服| 久久久精品国产亚洲av高清涩受| 又黄又粗又硬又大视频| 国产免费av片在线观看野外av| 欧美黑人巨大hd| 亚洲精品国产一区二区精华液| 女性被躁到高潮视频| 18禁黄网站禁片免费观看直播| 香蕉丝袜av| 国产成人欧美在线观看| 日韩精品中文字幕看吧| 国产精品免费一区二区三区在线| 久99久视频精品免费| 免费一级毛片在线播放高清视频| 老熟妇仑乱视频hdxx| 国产一区二区激情短视频| 黄色片一级片一级黄色片| 国产激情偷乱视频一区二区| 神马国产精品三级电影在线观看 | xxxwww97欧美| 亚洲午夜精品一区,二区,三区| 自线自在国产av| 91成人精品电影| 天堂动漫精品| 两人在一起打扑克的视频| 久久国产亚洲av麻豆专区| 亚洲av成人不卡在线观看播放网| 免费观看人在逋| 色综合欧美亚洲国产小说| 亚洲色图av天堂| 丰满的人妻完整版| 麻豆成人午夜福利视频| 人人妻人人看人人澡| 午夜福利在线在线| 日韩av在线大香蕉| 黄色丝袜av网址大全| 亚洲熟女毛片儿| 亚洲av美国av| 操出白浆在线播放| 99精品在免费线老司机午夜| 精品久久久久久久人妻蜜臀av| 国产精华一区二区三区| 国产色视频综合| 性欧美人与动物交配| 成人18禁高潮啪啪吃奶动态图| 级片在线观看| 国产激情欧美一区二区| 亚洲午夜理论影院| 无限看片的www在线观看| 在线十欧美十亚洲十日本专区| 日本熟妇午夜| 欧美精品亚洲一区二区| 亚洲一区二区三区色噜噜| 国产一级毛片七仙女欲春2 | 日韩av在线大香蕉| 好男人在线观看高清免费视频 | 日本a在线网址| 午夜成年电影在线免费观看| 亚洲中文字幕一区二区三区有码在线看 | 免费在线观看视频国产中文字幕亚洲| 久久亚洲精品不卡| 一级作爱视频免费观看| 757午夜福利合集在线观看| 男女做爰动态图高潮gif福利片| 免费看十八禁软件| 精品熟女少妇八av免费久了| 在线国产一区二区在线| av欧美777| 黄片大片在线免费观看| 女人爽到高潮嗷嗷叫在线视频| 搡老岳熟女国产| 一级片免费观看大全| 国产精品久久电影中文字幕| 久久久久久九九精品二区国产 | 久久久久久亚洲精品国产蜜桃av| 一级a爱视频在线免费观看| 免费在线观看完整版高清| 国产精品美女特级片免费视频播放器 | 给我免费播放毛片高清在线观看| 9191精品国产免费久久| 国产不卡一卡二| a在线观看视频网站| 免费观看精品视频网站| 国产成人系列免费观看| 一区二区三区精品91| 亚洲 欧美一区二区三区| 亚洲美女黄片视频| 禁无遮挡网站| www日本在线高清视频| 男女做爰动态图高潮gif福利片| 成年人黄色毛片网站| 国产成人精品久久二区二区免费| 国产主播在线观看一区二区| 久久午夜亚洲精品久久| 久久精品国产综合久久久| 精品国产一区二区三区四区第35| 亚洲成人久久爱视频| 老司机深夜福利视频在线观看| 亚洲国产毛片av蜜桃av| 俄罗斯特黄特色一大片| 午夜福利高清视频| 亚洲第一青青草原| 又大又爽又粗| 国产97色在线日韩免费| 亚洲va日本ⅴa欧美va伊人久久| 欧美日本亚洲视频在线播放| 亚洲av第一区精品v没综合| 午夜成年电影在线免费观看| 中文亚洲av片在线观看爽| 丰满的人妻完整版| 国产成人精品久久二区二区91| 99热这里只有精品一区 | 亚洲精品美女久久av网站| 欧美人与性动交α欧美精品济南到| 精品国产乱子伦一区二区三区| 欧美日韩精品网址| 国产91精品成人一区二区三区| 日韩高清综合在线| 亚洲第一电影网av| 国产精品一区二区精品视频观看| 亚洲第一av免费看| 99国产综合亚洲精品| 精品无人区乱码1区二区| 中文字幕高清在线视频| 亚洲 国产 在线| 天天添夜夜摸| 欧美中文日本在线观看视频| 757午夜福利合集在线观看| 国内揄拍国产精品人妻在线 | 国产在线精品亚洲第一网站| 伊人久久大香线蕉亚洲五| 精品乱码久久久久久99久播| 亚洲五月婷婷丁香| 丁香六月欧美| 日本五十路高清| 99国产精品一区二区蜜桃av| 亚洲国产中文字幕在线视频| 超碰成人久久| 一区二区日韩欧美中文字幕| 久久九九热精品免费| 黑人欧美特级aaaaaa片| 成年版毛片免费区| 国产91精品成人一区二区三区| 国产欧美日韩一区二区精品| 在线国产一区二区在线| 亚洲国产精品合色在线| 欧美午夜高清在线| 欧美人与性动交α欧美精品济南到| 麻豆一二三区av精品| 神马国产精品三级电影在线观看 | 国产欧美日韩一区二区精品| 波多野结衣巨乳人妻| 国产精品久久久久久亚洲av鲁大| 欧美日韩中文字幕国产精品一区二区三区| 亚洲成人久久爱视频| 国产亚洲精品一区二区www| 亚洲欧美日韩无卡精品| 亚洲aⅴ乱码一区二区在线播放 | 9191精品国产免费久久| 亚洲国产精品sss在线观看| 亚洲性夜色夜夜综合| 桃色一区二区三区在线观看| 国产精品乱码一区二三区的特点| 午夜福利欧美成人| 男人的好看免费观看在线视频 | 午夜日韩欧美国产| 精品国产超薄肉色丝袜足j| 国产野战对白在线观看| 人成视频在线观看免费观看| 一本精品99久久精品77| 久久婷婷成人综合色麻豆| 视频区欧美日本亚洲| 桃色一区二区三区在线观看| 欧美 亚洲 国产 日韩一| 日韩欧美一区视频在线观看| 午夜福利在线观看吧| 中文字幕久久专区| 波多野结衣高清无吗| 日韩av在线大香蕉| 成人18禁高潮啪啪吃奶动态图| 国产高清视频在线播放一区| 叶爱在线成人免费视频播放| a级毛片在线看网站| 最新美女视频免费是黄的| 久久久久久人人人人人| 亚洲狠狠婷婷综合久久图片| 亚洲av电影不卡..在线观看| 亚洲专区国产一区二区| 国产一卡二卡三卡精品| 黄色视频,在线免费观看| 国产精华一区二区三区| 精品欧美一区二区三区在线| 淫妇啪啪啪对白视频| 日韩国内少妇激情av| 国产区一区二久久| 久久国产亚洲av麻豆专区| 韩国精品一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 最新在线观看一区二区三区| 香蕉久久夜色| 91字幕亚洲| av在线天堂中文字幕| 99热6这里只有精品| 好看av亚洲va欧美ⅴa在| 国产精品久久久av美女十八| 欧美激情 高清一区二区三区| 欧美性猛交黑人性爽| 午夜福利在线观看吧| 女人被狂操c到高潮| 精品一区二区三区av网在线观看| 欧美乱色亚洲激情| 亚洲精品国产区一区二| 久久狼人影院| 色综合站精品国产| 免费看十八禁软件| 三级毛片av免费| 久久午夜综合久久蜜桃| 亚洲精品在线美女| 久久精品影院6| 日本一区二区免费在线视频| 看黄色毛片网站| 欧美日韩黄片免| 美女免费视频网站| 老汉色av国产亚洲站长工具| 日本a在线网址| 精品午夜福利视频在线观看一区| 成年版毛片免费区| 1024香蕉在线观看| av福利片在线| 一级毛片精品| 搡老岳熟女国产| 亚洲av成人av| 最近最新免费中文字幕在线| 狂野欧美激情性xxxx| av电影中文网址| 在线免费观看的www视频| 国产麻豆成人av免费视频| 国产精品乱码一区二三区的特点| 亚洲男人天堂网一区| 夜夜夜夜夜久久久久| 亚洲熟女毛片儿| 国产黄a三级三级三级人| av福利片在线| 一本大道久久a久久精品| av在线播放免费不卡| 黄色视频,在线免费观看| 免费在线观看亚洲国产| 亚洲色图av天堂| 日韩欧美国产一区二区入口| 久久精品国产亚洲av香蕉五月| 亚洲av电影在线进入| 中国美女看黄片| 99久久综合精品五月天人人| 国产午夜福利久久久久久| 欧美日韩中文字幕国产精品一区二区三区| 午夜激情福利司机影院| 又黄又爽又免费观看的视频| 亚洲国产欧洲综合997久久, | 国产视频一区二区在线看| av在线天堂中文字幕| 成人18禁在线播放| 男人舔女人的私密视频| 日本 欧美在线| 热re99久久国产66热| 男女做爰动态图高潮gif福利片| 99热6这里只有精品| 在线观看66精品国产| 国产欧美日韩一区二区三| 制服诱惑二区| 精品卡一卡二卡四卡免费| 国产欧美日韩一区二区三| 中文字幕精品亚洲无线码一区 | 国产黄片美女视频| 最好的美女福利视频网| 免费高清视频大片| 桃红色精品国产亚洲av| 校园春色视频在线观看| 无限看片的www在线观看| 精品久久久久久,| 啪啪无遮挡十八禁网站| 久久久久久亚洲精品国产蜜桃av| 黄片播放在线免费| 国产亚洲精品久久久久久毛片| 国产精品香港三级国产av潘金莲| 精品国产乱码久久久久久男人| 国产91精品成人一区二区三区| 国产高清有码在线观看视频 | 久久精品国产清高在天天线| 国产高清有码在线观看视频 | 国产精品香港三级国产av潘金莲| 色综合站精品国产| 欧美人与性动交α欧美精品济南到| 国产精品 欧美亚洲| 精品国产亚洲在线| 久久中文字幕人妻熟女| 精品久久久久久久久久久久久 | 免费看a级黄色片| 啪啪无遮挡十八禁网站| 亚洲欧美精品综合久久99| 丝袜在线中文字幕| 亚洲国产精品成人综合色| 一级毛片精品| 婷婷精品国产亚洲av在线| 搡老妇女老女人老熟妇| 怎么达到女性高潮| 欧美在线黄色| 好男人在线观看高清免费视频 | 国产精品精品国产色婷婷| 丰满人妻熟妇乱又伦精品不卡| 大型黄色视频在线免费观看| 免费人成视频x8x8入口观看| 男女那种视频在线观看| 波多野结衣高清无吗| 国产精品 国内视频| 精品第一国产精品| 亚洲 国产 在线| 亚洲九九香蕉| 精品久久久久久久久久免费视频| 精品久久久久久,| 51午夜福利影视在线观看| 日本三级黄在线观看| 一进一出抽搐gif免费好疼| 国产av不卡久久| 亚洲男人的天堂狠狠| 此物有八面人人有两片| 香蕉国产在线看| 黄色成人免费大全| 欧美性猛交╳xxx乱大交人| 欧美一级a爱片免费观看看 | 亚洲第一av免费看| 最近在线观看免费完整版| 国语自产精品视频在线第100页| 久久人妻福利社区极品人妻图片| 亚洲人成网站高清观看| 男人舔女人下体高潮全视频| 欧美乱妇无乱码| 激情在线观看视频在线高清| 男女床上黄色一级片免费看| 国产单亲对白刺激| 妹子高潮喷水视频| 久久中文字幕一级| 日日夜夜操网爽| 男女做爰动态图高潮gif福利片| 女警被强在线播放| 一区二区三区激情视频| 亚洲狠狠婷婷综合久久图片| 最新美女视频免费是黄的| 国产在线观看jvid| 黄色片一级片一级黄色片| 女性被躁到高潮视频| 一二三四社区在线视频社区8| 国产成人av教育| 波多野结衣巨乳人妻| 久久精品国产亚洲av香蕉五月| 久久久久国内视频| 精品久久久久久久久久免费视频| 亚洲国产精品999在线| 香蕉国产在线看| 国产精品99久久99久久久不卡| 国产精品精品国产色婷婷| 高潮久久久久久久久久久不卡| 亚洲国产看品久久| 夜夜夜夜夜久久久久| 午夜福利高清视频| 国产av又大| 国产爱豆传媒在线观看 | 亚洲国产欧美网| av在线天堂中文字幕| 757午夜福利合集在线观看| 99久久无色码亚洲精品果冻| 亚洲成人免费电影在线观看| 亚洲成人久久爱视频| 最好的美女福利视频网| 国产精品久久电影中文字幕| 欧美黑人欧美精品刺激| 国产激情偷乱视频一区二区| 国产伦人伦偷精品视频| 女同久久另类99精品国产91| 亚洲一卡2卡3卡4卡5卡精品中文| 国产av在哪里看| 中出人妻视频一区二区| 国产蜜桃级精品一区二区三区| 国产又黄又爽又无遮挡在线| 女性生殖器流出的白浆| 免费在线观看视频国产中文字幕亚洲| 国产精品久久久久久人妻精品电影| 成人三级黄色视频| 国产97色在线日韩免费| 国产三级在线视频| 两人在一起打扑克的视频| 久久久久久九九精品二区国产 | 麻豆av在线久日| 成年人黄色毛片网站| 精品少妇一区二区三区视频日本电影| 一边摸一边抽搐一进一小说| 亚洲五月色婷婷综合| 叶爱在线成人免费视频播放| 99久久99久久久精品蜜桃| 嫩草影院精品99| 中文资源天堂在线| 亚洲男人天堂网一区| 成人18禁在线播放| 丁香欧美五月| 可以在线观看的亚洲视频| 亚洲中文字幕一区二区三区有码在线看 | 亚洲第一欧美日韩一区二区三区| 久久欧美精品欧美久久欧美| 亚洲色图 男人天堂 中文字幕| 日本黄色视频三级网站网址| 欧美黑人欧美精品刺激| 成人手机av| 成人特级黄色片久久久久久久| 一个人观看的视频www高清免费观看 | www.熟女人妻精品国产| 人人妻人人澡欧美一区二区| 99热6这里只有精品| 日本免费一区二区三区高清不卡| 亚洲欧美一区二区三区黑人| 亚洲在线自拍视频| 男人舔女人下体高潮全视频| 无人区码免费观看不卡| 不卡一级毛片| 嫁个100分男人电影在线观看| 国产精品免费视频内射| 一个人观看的视频www高清免费观看 | 伊人久久大香线蕉亚洲五| av福利片在线| 丝袜在线中文字幕| 制服诱惑二区| 麻豆久久精品国产亚洲av| 午夜免费激情av| 国产精品98久久久久久宅男小说| 亚洲一码二码三码区别大吗| 亚洲va日本ⅴa欧美va伊人久久| 黄频高清免费视频| 亚洲国产日韩欧美精品在线观看 | 一本综合久久免费| 最近最新免费中文字幕在线| 亚洲 国产 在线| 午夜福利成人在线免费观看| 日本精品一区二区三区蜜桃| 无限看片的www在线观看| 亚洲精品久久国产高清桃花| 久久欧美精品欧美久久欧美| 国内久久婷婷六月综合欲色啪| 黄色视频不卡| 亚洲午夜精品一区,二区,三区| 真人一进一出gif抽搐免费| 亚洲精品一区av在线观看| 国产精品永久免费网站| 久久久久久大精品| 一二三四在线观看免费中文在| 免费在线观看黄色视频的| 18禁黄网站禁片免费观看直播| 日本一本二区三区精品| 欧美激情高清一区二区三区| 成人特级黄色片久久久久久久| 日韩一卡2卡3卡4卡2021年| 男女做爰动态图高潮gif福利片| 亚洲精品av麻豆狂野| 成人午夜高清在线视频 | 欧美中文日本在线观看视频| 亚洲第一青青草原| 国产亚洲精品一区二区www| 99riav亚洲国产免费| 午夜免费成人在线视频| 中文在线观看免费www的网站 | 成年免费大片在线观看| 久久天堂一区二区三区四区| 热re99久久国产66热| a级毛片在线看网站| 成人国产一区最新在线观看| 国产激情久久老熟女| 国产一区二区激情短视频| 白带黄色成豆腐渣| 三级毛片av免费| 国产免费av片在线观看野外av| 欧美大码av| 国产精品精品国产色婷婷| 亚洲人成网站在线播放欧美日韩| 99热6这里只有精品| 免费av毛片视频| 男男h啪啪无遮挡| 色精品久久人妻99蜜桃| 12—13女人毛片做爰片一| 狂野欧美激情性xxxx| 波多野结衣高清作品| 色婷婷久久久亚洲欧美| 欧美日韩一级在线毛片| 久热这里只有精品99| 热re99久久国产66热| 视频在线观看一区二区三区| 99国产极品粉嫩在线观看| 国产av一区在线观看免费| 嫩草影院精品99| 亚洲中文字幕日韩| 少妇的丰满在线观看| 免费高清视频大片| 黄色毛片三级朝国网站| 黑丝袜美女国产一区| 无限看片的www在线观看| 国产高清有码在线观看视频 | 亚洲精品国产一区二区精华液| 国产私拍福利视频在线观看| 怎么达到女性高潮| or卡值多少钱| 国产伦在线观看视频一区| 国产成人精品无人区| 黑人巨大精品欧美一区二区mp4| 俺也久久电影网| 亚洲电影在线观看av| 国产真实乱freesex| 精品人妻1区二区| 午夜福利成人在线免费观看| 久久99热这里只有精品18| 日本一本二区三区精品| 亚洲av日韩精品久久久久久密| av电影中文网址| 99久久无色码亚洲精品果冻| 少妇 在线观看| 1024手机看黄色片| 人人妻,人人澡人人爽秒播| 91麻豆av在线| 国产精品,欧美在线| 精品国产乱子伦一区二区三区| 精品久久久久久久毛片微露脸| 99国产精品一区二区三区| 99国产精品一区二区蜜桃av| 欧美日韩精品网址| 国产色视频综合| 久久久久亚洲av毛片大全| 精品一区二区三区四区五区乱码| 美女大奶头视频| 国产高清videossex| 天堂动漫精品| netflix在线观看网站| 久久久久久久精品吃奶| 精品国产亚洲在线| 成人三级做爰电影| 露出奶头的视频| 国产午夜福利久久久久久| 亚洲精品国产一区二区精华液| 久久久精品欧美日韩精品| 不卡一级毛片| 久久伊人香网站| 免费观看人在逋| 日本免费一区二区三区高清不卡| 中文字幕最新亚洲高清| 狠狠狠狠99中文字幕| 此物有八面人人有两片| 黄片播放在线免费| 母亲3免费完整高清在线观看| 国产亚洲精品第一综合不卡| 99精品在免费线老司机午夜| 久久精品成人免费网站| 美女免费视频网站| 人人妻人人看人人澡| 桃红色精品国产亚洲av| 一a级毛片在线观看| netflix在线观看网站|