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

    內(nèi)蒙古東部地區(qū)面波震級與近震震級轉(zhuǎn)換關(guān)系

    2016-11-26 04:50:41芳劉瑞豐帆張暉趙艷紅李彬娜
    地震地磁觀測與研究 2016年5期

    劉 芳劉瑞豐,2)張 帆張 暉趙艷紅李 彬娜 熱

    1)中國呼和浩特010010內(nèi)蒙古自治區(qū)地震局

    2)中國北京100081中國地震局地球物理研究所

    內(nèi)蒙古東部地區(qū)面波震級與近震震級轉(zhuǎn)換關(guān)系

    劉 芳1)劉瑞豐1),2)張 帆1)張 暉1)趙艷紅1)李 彬1)娜 熱1)

    1)中國呼和浩特010010內(nèi)蒙古自治區(qū)地震局

    2)中國北京100081中國地震局地球物理研究所

    基于震級殘差統(tǒng)計、線性回歸和正交回歸方法,選取2008年至2016 年3月內(nèi)蒙古測震臺網(wǎng)記錄的內(nèi)蒙古東部地區(qū)地震觀測報告和波形事件,確立該區(qū)量規(guī)函數(shù),建立近震震級與面波震級間的轉(zhuǎn)換關(guān)系。研究表明,HLH、HLR和AGL等5個臺站震級偏差較大,可能由地震臺站場地響應(yīng)衰減或放大、臺基風(fēng)化等造成;全國量規(guī)函數(shù)僅在140 km<⊿≤280 km時適用于內(nèi)蒙古東部地區(qū);采用正交回歸方法(OR)得到古東部地區(qū)MS東部=0.99ML-0.12。

    震級;量規(guī)函數(shù);近震震級ML;面波震級MS;震級轉(zhuǎn)換

    0 引言

    震級是表征地震大小的量,是地震基本參數(shù)之一,是地震預(yù)報和其他有關(guān)地震研究的一個重要參數(shù)(陳運(yùn)泰等,2000;陳運(yùn)泰等,2004)。表征地震波通過和介質(zhì)有地區(qū)性差異或特殊構(gòu)造的量規(guī)函數(shù)是影響震級測定的重要因素之一。多年實測結(jié)果表明,中國使用1971年建立的R1(⊿)和R2(⊿)尚不能充分反映中國地殼的地區(qū)性差異,不同區(qū)域的量規(guī)函數(shù)急需修定。陳培善等(1983)、嚴(yán)尊國(1982,1983,1986,1992,1995)獲得中國不同地質(zhì)、地震單元的差異和中國東部地區(qū)和西部量規(guī)函數(shù);郭履燦(1984)、秦嘉政等(1986)、彭美鳳(1982)、張誠(1981)、黃緯瓊(1982)、薛志照(1992)、高國英(1987)、趙明淳等(2005)、孟曉琴等(2008)等開展測定震級和量規(guī)函數(shù)關(guān)系研究,取得一些成果。

    中國目前使用的面波震級MS和近震震級ML兩種標(biāo)度(張宏志等,2007;汪素云,2010)間的轉(zhuǎn)換關(guān)系,是1971年郭履燦得到的適用于中國華北地區(qū)經(jīng)驗關(guān)系式MS經(jīng)驗=1.13ML-1.08(震中距⊿≤1 000 km時)①郭履燦.華北地區(qū)的地方性震級ML和面波震級MS經(jīng)驗關(guān)系.全國地震工作會議資料,1971:1-10.。考慮到不同震級測量的方法不同,因此在地震臺網(wǎng)的震級測定中,不同震級之間一律不進(jìn)行換算(劉瑞豐等,2006)。然而,目前區(qū)域臺網(wǎng)大震速報、地震編目、地震預(yù)測預(yù)報及相關(guān)地球科學(xué)研究中使用的震級M,均是由經(jīng)驗公式轉(zhuǎn)換得到的。多年實測表明,ML和MS間轉(zhuǎn)換關(guān)系已不符合經(jīng)驗公式,存在系

    統(tǒng)偏差,需要矯正。理論上,采用正交回歸方法研究各種震級標(biāo)度間的關(guān)系接近實際震級測定 (Madansky,1959; Fuller,1987; Carroll R I et al,1996;劉瑞豐等,2007)。劉瑞豐等(1996,2006,2007)、楊晶瓊等(2013)和張誠(1981)等均開展了相關(guān)研究,取得較好成果。

    內(nèi)蒙古東部地區(qū)(40°—53°N,114°—126°E)位于華北板塊與東北板塊交匯處,構(gòu)造復(fù)雜、構(gòu)造運(yùn)動強(qiáng)烈。據(jù)歷史地震記載及近代中強(qiáng)地震記錄,曾發(fā)生M≥5.0地震20次,M≥6.0地震3次,且M≥5.0地震主要集中于該區(qū)域東部和北部的阿倫河斷裂、雅魯河斷裂及東南部的八里罕斷裂、赤峰—開原斷裂和北票—朝陽斷裂一帶(張暉等,2013)。這是該區(qū)域歷史地震發(fā)生的關(guān)鍵因素,同時表明該區(qū)域內(nèi)具備發(fā)生中強(qiáng)地震的地質(zhì)構(gòu)造條件(曹鳳娟,2013)。

    遵循繼承性、區(qū)域性、實用性震級確定三原則,基于震級殘差統(tǒng)計分析方法、線性回歸和正交回歸方法,采用2008年至2016年3月內(nèi)蒙古東部及鄰省地震觀測報告和ML≥3.8數(shù)字地震波形資料,建立內(nèi)蒙古東部地區(qū)區(qū)域量規(guī)函數(shù)、近震震級與面波震級間的轉(zhuǎn)換關(guān)系。同時,開展新量規(guī)函數(shù)的檢驗、近震震級與面波震級間轉(zhuǎn)換關(guān)系高斯隨機(jī)擾動檢驗,驗證其適用性、可靠性,為內(nèi)蒙古地震預(yù)報及地震學(xué)相關(guān)研究提供科學(xué)參考。

    1 資料選取

    選用2008—2016年3月內(nèi)蒙古東部地區(qū)及邊境線50 km以內(nèi)2 149次ML≥0.0地震事件[圖1(a)]作為建立新量規(guī)函數(shù)的基礎(chǔ)數(shù)據(jù)。震級測定中,采用MSDP定位軟件,進(jìn)行WA仿真,量取兩水平向最大位移值,使得周期值得到較好約束,基本消除周期測定對震級計算造成的偏差,根據(jù)近震震級計算公式求取各子臺震級值,由多臺震級平均值得到臺網(wǎng)震級值?;谡鸺墯埐罱y(tǒng)計方法,篩選參與定位的至少4個的臺站資料,得到臺網(wǎng)多臺震級平均值,重新修定震級的量規(guī)函數(shù)。選取68次ML≥3.8數(shù)字地震波形資料[圖1(b)],作為建立近震震級與面波震級間轉(zhuǎn)換關(guān)系的基礎(chǔ)數(shù)據(jù),近震震級ML、面波震級MS均按照新國標(biāo)《地震震級的規(guī)定》(GB17740)進(jìn)行測定。

    圖1 地震和臺站分布(a) 2 149次ML≥0.0地震; (b) 68次ML≥3.8地震Fig.1 Distribution of stations and earthquakes

    2 計算方法

    2.1 量規(guī)函數(shù)

    根據(jù)《測震臺網(wǎng)運(yùn)行管理細(xì)則》(中震測函[2015]115號),采用(1)式來計算近震震級。所有的數(shù)字記錄均仿真成短周期地震儀WA后,在位移記錄上量取兩水平向最大位移值。

    式中: Aμ為水平分向記錄的最大地動位移的算術(shù)平均值(單位μm); R(⊿)為短周期地震儀測定ML的量規(guī)函數(shù)。

    對第i個地震,求出對所有臺站的震級平均值

    式中,Ns為臺站個數(shù),Ne為地震個數(shù) 。

    將單個臺站對應(yīng)的單個地震震級減去單個地震震級平均值,可得到單臺震級的殘差值。第j個臺站記錄地震的震級標(biāo)準(zhǔn)偏差δj和平均偏差為

    單臺震級偏差值隨震中距的變化

    式中: k是震中距間隔序數(shù),N為間隔范圍內(nèi)震級偏差樣本數(shù)。

    2.2 回歸方法

    2.2.1 線性回歸方法。對于兩個或多個存在統(tǒng)計相關(guān)的隨機(jī)變量,可依據(jù)大量觀測數(shù)據(jù)來確定它們之間統(tǒng)計的定量關(guān)系,建立回歸方程。將 N 個數(shù)據(jù)點 (xi,yi),其中i=1,2,3,…,N,擬合的直線模型為

    隨著自變量和因變量的不同選擇,對于相同樣本點會得到不同的回歸方程。選擇X方向或Y方向擬合,所得回歸直線是不同的(孫彥清,2002;劉瑞豐,2007)。根據(jù)測量誤差判斷誤差方向,確定式(6)中的A和B,用誤差方向確定擬合方向(黃杰等,2000;孫彥清,2002)。

    通常利用線性最小二乘回歸(SR)方法確定系數(shù)A和B。古登堡和里克特(1956)采用SR方法得到MS和ML之間的關(guān)系式。確定系數(shù) A和B擬合線性公式(6)有以下兩種可能 (Draper N R et al,1998)。

    適用條件為:σ2xx→0,σ2yy>0(σ2xx、σ2yy分別是x、y的方差)。

    SR2稱為反標(biāo)準(zhǔn)回歸(Carroll R I et al,1996),適用條件為:

    線性回歸方法只考慮一個變量比另一個變量產(chǎn)生的偏差大,線性回歸關(guān)系不能使用等號,因此在公式(7)和公式(8)中不使用等號而使用箭頭符號。

    2.2.2 正交回歸方法。當(dāng)X和Y方向上測量誤差均不可忽略,不能在單一方向上進(jìn)行擬合時,應(yīng)滿足各測量點到擬合直線垂直距離的平方和最小,即正交擬合(吳俊林等,1992;劉魏等,1994;李雄軍,2005a,2005b;姜慧等,2006)。正交回歸方法能夠克服固定單方向最優(yōu)帶來的擬合穩(wěn)定性差的弊端。

    采用正交回歸的方法擬合公式(6),以O(shè)R表示正交回歸方法 ,則

    通常采用 Hesse方法來表示正交回歸 (Carroll R I et al,1996) ,即把兩個變量放在等號右邊,表示兩個變量均在變化,即

    把上述關(guān)系應(yīng)用于震級標(biāo)度Mx和My,則SR1適用于在Mx測量偏差較大情況下對My產(chǎn)生的影響,而 SR2適用于在My測量偏差較大的情況下對Mx產(chǎn)生的影響。然而,震級測量時均可能存在一定誤差,SR1和SR2均會與實際測定震級有一定差別。理論上,使用正交回歸方法研究各種震級標(biāo)度之間的關(guān)系更能接近實際的震級測定 (Madansky A,1959; Carroll R I et al,1996;劉瑞豐等,2007)。

    3 計算結(jié)果及分析

    3.1 量規(guī)函數(shù)

    采用內(nèi)蒙古東部地區(qū)36個子臺記錄的11 101個地震數(shù)據(jù),基于震級殘差統(tǒng)計方法,計算得到各子臺震級相對于臺網(wǎng)震級Mi的偏差、平均偏差和標(biāo)準(zhǔn)偏差δi。對臺站引起的震級偏差進(jìn)行臺基校正,并對比校正前后單臺震級偏差值隨震中距的變化、震級與臺網(wǎng)平均震級偏差值分布,依據(jù)所得結(jié)果建立東部地區(qū)量規(guī)函數(shù)。

    3.1.1 臺基影響。表1給出東部地區(qū)36個子臺震級相對于臺網(wǎng)震級Mi的偏差、平均偏差和標(biāo)準(zhǔn)偏差δi。從表1可知,其中18個臺站不存在震級偏差;有13個臺站震級偏差范圍為0.20—0.26,占總數(shù)的36.1%;有5個臺站(內(nèi)蒙古HLH、HLR和AGL;河北LOH和WEC)的震級偏差較大,偏差范圍在-0.321—0.376,占總數(shù)的13.9%??梢?,臺基是引起近震級偏差的重要因素之一。

    3.1.2 單臺震級與臺網(wǎng)平均震級偏差值分布。分析臺基校正前后單臺震級與平均震級偏差值(簡稱偏差值σ)變化可知:①σ分布形態(tài)。臺基校正前[圖2(a)],偏差值基本呈正態(tài)分布,離散度比較大,σ分布范圍是-1.8—1.4。臺基校正后[圖2(b)],偏差值基本呈正態(tài)分布,離散度較小,σ分布范圍是-1.5—1.2,分布更集中;②震級偏差的標(biāo)準(zhǔn)差。臺基校正前為0.295 0,臺基校正校正后為0.271 0,降低0.024。

    表1 各子臺計算得到的平均偏差和標(biāo)準(zhǔn)偏差δiTable 1 Average deviation and standard deviation calculated from each station

    表1 各子臺計算得到的平均偏差和標(biāo)準(zhǔn)偏差δiTable 1 Average deviation and standard deviation calculated from each station

    序號 臺網(wǎng)代碼臺站代碼δi樣本數(shù) 序號 臺網(wǎng)代碼臺站代碼δi樣本數(shù)1 NM BAC 0.08 0.23 311 19 JL BCT 0.16 0.22 263 2 NM XLT 0.02 0.29 337 20 NM AGL-0.36 0.32 155 3 NM XIQ 0.08 0.22 244 21 NM ZLT-0.03 0.19 748 4 NM MZL-0.09 0.26 207 22 NM IDR-0.20 0.21 483 5 NM JIP-0.08 0.24 393 23 NM NJT 0.18 0.20 615 6 NM LIX 0.00 0.32 351 24 HE LOH-0.18 0.31 228 7 NM NIC-0.07 0.25 365 25 HE WEC-0.13 0.31 324 8 NM CHF-0.12 0.27 240 26 LN LYA 0.40 0.24 193 9 NM HLH-0.30 0.31 161 27 LN BEP 0.14 0.29 229 10 NM HLR-0.32 0.26 313 28 LN FXI 0.27 0.21 233 11 NM XIH-0.08 0.32 273 29 LN HXQ 0.11 0.26 153 12 NM ARS 0.08 0.22 570 30 LN FKU 0.38 0.20 135 13 NM TIS-0.16 0.23 379 31 LN XFN 0.17 0.25 113 14 NM MDG 0.12 0.20 282 32 HL NZN 0.16 0.20 674 15 NM LUB-0.07 0.23 292 33 HL JGD 0.08 0.23 289 16 NM CHR-0.10 0.22 758 34 HL ZHY-0.21 0.29 113 17 NM GNH-0.06 0.23 346 35 HL NEH-0.05 0.21 268 18 NM WLT-0.02 0.22 280 36 HL MOH-0.07 0.24 56

    圖2 單臺震級與平均震級偏差標(biāo)準(zhǔn)差值統(tǒng)計(a)臺基校正前; (b)臺基校正后Fig.2 Statistics of magnitude deviation between single station and average one

    3.2 東部地區(qū)量規(guī)函數(shù)校正

    扣除臺基影響,修定震中距0—500 km的全國量規(guī)函數(shù),得到內(nèi)蒙古東部地區(qū)量規(guī)函數(shù),見表2和圖4。從表2和圖4可知,當(dāng)0≤⊿≤140 km時,全國量規(guī)函數(shù)偏高;當(dāng)140 km<⊿≤280 km時,全國量規(guī)函數(shù)基本適合;當(dāng)⊿>280 km時,東部量規(guī)函數(shù)高于全國量規(guī)函數(shù)。

    圖4 內(nèi)蒙古東部地區(qū)新量規(guī)函數(shù)曲線Fig.4 New calibration function in eastern Inner Mongolia region

    表2 全國及內(nèi)蒙古東部地區(qū)新量規(guī)函數(shù)R(⊿)Table 2 The new eastern Inner Mongolia calibration function R(⊿) and the China calibration function

    3.3 東部地區(qū)量規(guī)函數(shù)檢驗

    以震級偏差的標(biāo)準(zhǔn)差δi和震級偏差分布集中度做為檢驗指標(biāo),分別對臺基校正前、臺基校正后以及新量規(guī)函數(shù)3個階段進(jìn)行對比,檢驗東部量規(guī)函數(shù)的適用性、可靠性。結(jié)果表明,臺基校正前δi=0.295 0,臺基校正后δi=0.271 0,使用新量規(guī)函數(shù)測定的δi=0.237 4。由圖5可見,盡管在臺基校正前、臺基校正后、以及使用新量規(guī)函數(shù)3種情況下,震級偏差分布集中度均呈正態(tài)分布,但是采用新量規(guī)函數(shù)后,震級偏差分布集中度較前兩種情況有了明顯提高。分析認(rèn)為,使用東部量規(guī)函數(shù)測定的δi有一定降低,震級偏差分布集中度明顯提高。

    圖5 單臺震級與平均震級偏差標(biāo)準(zhǔn)差值統(tǒng)計(a)臺基校正前; (b)臺基校正后; (c)新量規(guī)函數(shù)Fig.5 Statistics of magnitude deviation between single station and average one

    4 ML與MS轉(zhuǎn)換關(guān)系

    4.1 東部地區(qū)ML和MS關(guān)系的建立

    依據(jù)上述線性回歸方法(SR1和SR2)和正交回歸方法(OR),重新測定內(nèi)蒙古東部地區(qū)2008—2016年3月68次MS≥3.8地震震級ML和MS(共132個震級對)。表3列出經(jīng)過線性回歸和正交回歸得到的ML和MS震級關(guān)系式,圖6(a)展示線性回歸、正交回歸以及經(jīng)驗關(guān)系式得到4條ML與MS關(guān)系曲線。由表3和圖6(a)可見:正交回歸方法(OR)得到的紅色曲線位于數(shù)據(jù)點中間,均方根誤差最小,為最佳擬合曲線; OR曲線位于經(jīng)驗曲線上方。通過擬合,得到內(nèi)蒙古東部地區(qū)近震震級ML與面波震級MS的轉(zhuǎn)換關(guān)系

    分析MS經(jīng)驗與MS東部差值數(shù)量分布[圖6(b)],發(fā)現(xiàn)絕大多數(shù)地震震級差值分布于0.15—0.45,0.30—0.45范圍內(nèi)的地震數(shù)量最多。

    表3 ML和MS關(guān)系式Table 3 Relation of MLand MS

    圖6 震級關(guān)系及震級差值數(shù)量分布(a) ML和MS震級關(guān)系; (b) MS東部與MS經(jīng)驗差值數(shù)量分布Fig.6 Magnitude relation and magnitude difference numerical distribution

    4.2 兩種震級轉(zhuǎn)換關(guān)系對比分析

    通過對比MS東部和MS經(jīng)驗及差值數(shù)量分布,獲得以下認(rèn)識:①從兩種震級轉(zhuǎn)換關(guān)系式看,差別較大。關(guān)系式中系數(shù)A由1.13降為0.96,降低0.17,變化不大;系數(shù)B由1.08降到0.10,降低0.98,變化較大,兩者明顯存在系統(tǒng)偏差,可能是由于區(qū)域構(gòu)造特征差異性所致;②對比MS經(jīng)驗與MS東部差值數(shù)量分布[圖6(b)]可知,大多數(shù)地震震級差值分布在0.30—0.45。

    5 結(jié)論

    采用內(nèi)蒙古測震臺網(wǎng)地震觀測報告和數(shù)字地震波形,建立內(nèi)蒙古東部地區(qū)量規(guī)函數(shù)、近震震級與面波震級間轉(zhuǎn)換關(guān)系,得到以下結(jié)論。

    (1)由36個子臺震級相對于臺網(wǎng)震級Mi的偏差、平均偏差和標(biāo)準(zhǔn)偏差δi可知,18個臺站不存在震級偏差;13個臺站震級偏差范圍在0.20—0.26,占總數(shù)的36.1%;5個臺站(內(nèi)蒙古HLH、HLR和AGL,河北LOH和WEC)震級偏差較大,偏差范圍在-0.321—0.376,占總數(shù)的13.9%。分析造成震級偏差的原因認(rèn)為,HLH臺站的場地響應(yīng)低頻放大,這可能是受臺基位于軟的沉積層場地影響,臺基風(fēng)化較為嚴(yán)重;AGL臺站場地響應(yīng)呈現(xiàn)高頻衰減效應(yīng);HLR臺場在響應(yīng)呈現(xiàn)低頻明顯衰減(王鑫等,2015);河北境內(nèi)的LOH和WEC臺站,場地響應(yīng)值較大。

    (2)對比可知,當(dāng)0≤⊿≤140 km時,全國量規(guī)函數(shù)偏高;當(dāng)140<⊿≤280 km時,全國量規(guī)函數(shù)基本適合;當(dāng)⊿>280 km時,東部量規(guī)函數(shù)高于全國量規(guī)函數(shù)。通過對東部量規(guī)函數(shù)檢驗可知,使用東部地區(qū)量規(guī)函數(shù)測定的震級標(biāo)準(zhǔn)差有了一定降低,且震級偏差分布集中度較臺基校正前后有了明顯提高。可見,全國量規(guī)函數(shù)不符合內(nèi)蒙古東部地區(qū)區(qū)域特征。

    (3)對比MS經(jīng)驗=1.13ML-1.08和MS東部=0.99ML-0.12 可知,兩者明顯存在系統(tǒng)偏差,可能是由于區(qū)域構(gòu)造特征差異性所致。結(jié)果表明,經(jīng)驗轉(zhuǎn)換關(guān)系已不符合內(nèi)蒙古東部地區(qū)區(qū)域特征,建議使用本研究建立的ML與MS轉(zhuǎn)換關(guān)系式。

    (4)2008年以來,東部地區(qū)ML≥3.8地震數(shù)量有限,且地震分布、臺站布局不均勻,因此,新建立的東部地區(qū)量規(guī)函數(shù)、近震震級與面波震級間的轉(zhuǎn)換關(guān)系有待完善,以便符合內(nèi)蒙古東部地區(qū)區(qū)域特征。

    中國地震局地球物理研究所劉瑞豐研究員提供計算程序,并對本研究給予支持和幫助,在此表示感謝。

    曹鳳娟,焦明若,王亮,等.東北地區(qū)M≥5.0地震前中等地震集中活動的震兆研究[J].地震工程學(xué)報,2013,35(2):385-393.

    陳培善,秦嘉政.量規(guī)函數(shù)、臺站方位、臺基及不同測量方法對近震震級ML的影響[J].地震學(xué)報,1983,5(1):87-98.

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

    陳運(yùn)泰,吳忠良,王培德,等.數(shù)字地震學(xué)[M].北京:地震出版社,2000:1-30

    高國英,艾尼瓦爾.新疆地區(qū)量規(guī)函數(shù)、臺站校正值的分析計算[J].西北地震學(xué)報,1987,9(1):103-107.

    本工程采用的混凝土是商用混凝土,要求混凝土坍落度控制在180mm,上下波動不大于30mm。由于坡屋面施工時用水養(yǎng)護(hù)不易,需要在混凝土中適量混入緩凝劑,并使用薄膜包裹。養(yǎng)護(hù)不少于14d,避免混凝土出現(xiàn)開裂問題。

    郭履燦,趙鳳竹,趙共玲,等.震級與震源參數(shù)測定[M].北京:中國科學(xué)技術(shù)出版社,1984.

    黃杰,蔡希潔,林尊琪.一種對稱的線性擬合數(shù)據(jù)處理方法[J].計量技術(shù),2000,(5):47-49.

    黃緯瓊.近震衰減特征的地區(qū)差異及近震震級的估算[J].地震學(xué)報,1982,4(4):337-349.

    姜慧,唐麗華,梁海華.回歸分析方法在地震科學(xué)應(yīng)用中的問題與探討[J].華南地震,2006,26(3):1-7.

    李雄軍.對X和Y方向最小二乘線性回歸的討論[J].計量技術(shù),2005a,(1):50-52.

    李雄軍.幾種線性回歸方法的比較[J].計量技術(shù),2005b,(8):52-54.

    劉瑞豐,陳運(yùn)泰,Peter Bormann,等.中國地震臺網(wǎng)與美國地震臺網(wǎng)測定震級的對比Ⅱ—面波震級[J].地震學(xué)報,2006,28(1):1-7.

    劉瑞豐,陳運(yùn)泰,任梟,等.中國地震臺網(wǎng)震級的對比[J].地震學(xué)報,2007,29(5):467-476.

    劉瑞豐,黨京平,陳培善.利用速度數(shù)字地震儀記錄測定面波震級[J].地震地磁觀測與研究,1996,17(2):1-4.

    孟曉琴,梁向軍,張玲,等.山西地震臺網(wǎng)速報震級與量規(guī)函數(shù)的研究[J].山西地震,2008,2:19-23.

    彭美鳳,閩臺地區(qū)量規(guī)函數(shù)R(⊿)的探討[J].地震地磁觀測與研究,1982,3(3):71-75.

    秦嘉政,等.不同儀器以及地震波的區(qū)域性衰減特征對測定ML的影響[J].地震研究,1986,9(1):13-25.

    孫彥清.最小二乘線性擬合應(yīng)注意的兩個問題[J].漢中師范學(xué)院學(xué)報(自然科學(xué)),2002,20(1):58-61.

    汪素云,王健,俞言祥,等.基于中國地震臺網(wǎng)觀測報告的ML與MS經(jīng)驗關(guān)系[J].中國地震,2010,26(1):14-22.

    王鑫,尹戰(zhàn)軍,郝美仙,等.內(nèi)蒙古地區(qū)場地響應(yīng)區(qū)域特征分析[J].震害防御技術(shù),2015,10(1):163-172.

    吳俊林,王社柱,李樹華.直線擬合方向選擇的應(yīng)用分析[J].陜西師大學(xué)報,1992,20(3):72-75.

    薛志照.近震震級ML量規(guī)函數(shù)和臺站校正值[J].華北地震科學(xué),1992,108 (1):60-651.

    嚴(yán)尊國,等,地方性均勻震級系統(tǒng)量規(guī)函數(shù)的區(qū)域特征[J].地震地磁觀測與研究,1982,4(4):4-10.

    嚴(yán)尊國,等,湖北區(qū)域臺網(wǎng)測定ML震級中量規(guī)函數(shù)的修正[J].地震學(xué)刊,1983,2:9-17.

    嚴(yán)尊國,建立華北南部地方性均勻震級系統(tǒng)的嘗試[J].中國地震,1986,2(2):57-63.

    嚴(yán)尊國,李普麗,薛軍蓉.中國近震震級量規(guī)函數(shù)檢驗研究[J].華南地震,1995,15(1):1-11.

    嚴(yán)尊國,李普麗,薛軍蓉.中國近震震級量規(guī)函數(shù)研究[J].中國地震,1992,8(4):76-91.

    楊晶瓊,楊周勝 蔡明軍,等.云南地區(qū)近震震級與面波震級轉(zhuǎn)換關(guān)系研究[J].中國地震,2013,29(4):513-521

    張誠.西北地區(qū)測定震級ML中的某些問題[J].地震地磁觀測與研究,1981,3(2):28-31.

    張宏志,刁桂苓 ,趙明淳,等.不同標(biāo)度震級關(guān)系和臺基影響問題探討[J].中國地震,2007,23(2):141-146.

    張暉,高立新,韓曉明.內(nèi)蒙古東部地區(qū)波速比變化特征研究[J].地震工程學(xué)報,2013,35(4):921-927.

    趙明淳,毛國良,張叢珍,等.河北遙測臺網(wǎng)速度震級偏差與量規(guī)函數(shù)關(guān)系的研究[J].地震地磁觀測與研究,2005,26(5):87-94.

    Carroll R I,Ruppert D.The use and misuse of orthogonal regression in linear errors-in-variables models[J].The American Statistician,1996,50(1):1-6.

    Draper N R,Smith H.Applied Regression Analysis[M].New York:John Wiley & Sons,1998: 100-125

    Gutenberg B,Richter C F.Earthquake magnitude,intensity,energy and acceleration[J].Bull Seism Soc Amer,1956a,46: 105-145.

    Gutenberg B,Richter C F.Magnitude and energy of earthquakes[J].Annali di Geofisica,1956b,91: 1-15.

    Madansky A.The fitting of straight lines when both variables are subject to error[J].J Amer Statist Assoc,1959,54: 173-205.

    Contrast of MLand MSin eastern Inner Mongolia Region

    Liu Fang1),Liu Ruifeng1),2),Zhang Fan1),Zhang Hui1),Zhao Yanhong1),Li Bin1)and Na Re1)
    1) Earthquake Administration of Inner Mongolia Autonomous,Huhhot 010010,China
    2) Institute of Geophysics,China Earthquake Administration,Beijing 100081,China

    Using report and seismic events recorded on eastern Inner Mongolia digital seismic network since 2008 to March 2016 and residual statistical,linear regression and orthogonal regression methods,MLcalibration function was modified and the relationship between MLand MSis corrected.The higher residual in fi ve stations,such as HLH,HLR and AGL,may be caused by the amplifi cation or attenuation of station site response,or stylobate weathering.The china calibration function is suitable for Inner Mongolia region with epicentral distance between 140 km and 60 km,yet unsuitable in other case.Using orthogonal regression method (OR) in the eastern regions of Inner Mongolia,MS=0.9 ML-0.12 relationship is found.

    magnitude,calibration function,near earthquake magnitude ML,surface wave magnitude MS,magnitude conversion

    10.3969/j.issn.1003-3246.2016.05.001

    劉芳(1963—),女,遼寧大連人,副研究員,主要從事地震監(jiān)測、地震預(yù)報和地震學(xué)研究工作。

    E-mail: lfnm88@163.com

    中國地震局地震科技星火項目(項目編號:XH15008)

    本文收到日期:2015-07-28

    欧美+日韩+精品| 搡老熟女国产l中国老女人| 亚洲在线观看片| 精品人妻一区二区三区麻豆 | 欧美激情在线99| 亚洲国产精品成人综合色| 亚洲国产欧洲综合997久久,| 哪里可以看免费的av片| 有码 亚洲区| 乱人视频在线观看| а√天堂www在线а√下载| 女人十人毛片免费观看3o分钟| 老鸭窝网址在线观看| 色精品久久人妻99蜜桃| 性色avwww在线观看| tocl精华| 国产美女午夜福利| 国产美女午夜福利| 在线免费观看不下载黄p国产 | 五月玫瑰六月丁香| 国产亚洲欧美在线一区二区| 欧美一区二区国产精品久久精品| 亚洲av日韩精品久久久久久密| 听说在线观看完整版免费高清| 91在线精品国自产拍蜜月 | 亚洲专区国产一区二区| 不卡一级毛片| av天堂在线播放| 国产一区在线观看成人免费| 一进一出好大好爽视频| 国产精品一区二区三区四区久久| 一个人免费在线观看的高清视频| 亚洲av熟女| 禁无遮挡网站| 在线看三级毛片| 亚洲人成网站在线播| 精品久久久久久久久久免费视频| 人人妻人人看人人澡| 少妇熟女aⅴ在线视频| 在线观看免费视频日本深夜| 国内精品美女久久久久久| 国产精品影院久久| 午夜a级毛片| 天堂动漫精品| 黄色视频,在线免费观看| 91久久精品国产一区二区成人 | 国产激情欧美一区二区| 国产免费男女视频| 一区二区三区激情视频| 好男人电影高清在线观看| 搡老岳熟女国产| 18禁国产床啪视频网站| 一个人免费在线观看的高清视频| 天天躁日日操中文字幕| 久9热在线精品视频| 欧美黄色淫秽网站| 黄色丝袜av网址大全| 可以在线观看毛片的网站| 伊人久久精品亚洲午夜| 日本黄色视频三级网站网址| 欧美又色又爽又黄视频| 国产主播在线观看一区二区| 亚洲熟妇熟女久久| 成人午夜高清在线视频| 1024手机看黄色片| 日韩欧美一区二区三区在线观看| 91麻豆精品激情在线观看国产| 国产一区二区亚洲精品在线观看| 18禁美女被吸乳视频| 午夜免费激情av| av片东京热男人的天堂| 级片在线观看| 国产熟女xx| 欧美日韩综合久久久久久 | 日韩有码中文字幕| 精品午夜福利视频在线观看一区| 九色国产91popny在线| 欧美乱码精品一区二区三区| 老司机在亚洲福利影院| 国产精品综合久久久久久久免费| 18禁黄网站禁片免费观看直播| 美女高潮喷水抽搐中文字幕| 欧美性感艳星| 男人的好看免费观看在线视频| e午夜精品久久久久久久| 最近最新中文字幕大全电影3| 在线国产一区二区在线| 欧美区成人在线视频| 国产男靠女视频免费网站| 亚洲激情在线av| 国产精品精品国产色婷婷| 十八禁人妻一区二区| 国产日本99.免费观看| 黄色丝袜av网址大全| 香蕉久久夜色| 精品一区二区三区人妻视频| 91av网一区二区| 九色成人免费人妻av| 日韩av在线大香蕉| 成人av在线播放网站| 熟女少妇亚洲综合色aaa.| 美女被艹到高潮喷水动态| 一个人看的www免费观看视频| 97碰自拍视频| 国产精品一及| 在线观看美女被高潮喷水网站 | 性色av乱码一区二区三区2| 久久午夜亚洲精品久久| 色av中文字幕| 久久久久久久久中文| 三级男女做爰猛烈吃奶摸视频| 国产av在哪里看| 日韩欧美三级三区| 男女之事视频高清在线观看| 欧美色视频一区免费| 国产精品久久电影中文字幕| 日韩欧美在线乱码| 成人特级黄色片久久久久久久| 国产三级在线视频| 国产欧美日韩精品亚洲av| 色播亚洲综合网| 一级黄片播放器| 日本一本二区三区精品| 搞女人的毛片| 日韩人妻高清精品专区| 宅男免费午夜| 两个人视频免费观看高清| 热99re8久久精品国产| 久久人妻av系列| 熟妇人妻久久中文字幕3abv| 悠悠久久av| 他把我摸到了高潮在线观看| 免费看美女性在线毛片视频| 99国产综合亚洲精品| 亚洲熟妇熟女久久| www.999成人在线观看| 色综合亚洲欧美另类图片| 欧美bdsm另类| 色在线成人网| 伊人久久大香线蕉亚洲五| 亚洲精品在线美女| 九九在线视频观看精品| 国语自产精品视频在线第100页| 99国产精品一区二区三区| 欧美+日韩+精品| 国产伦在线观看视频一区| 18+在线观看网站| 精品久久久久久久人妻蜜臀av| 高潮久久久久久久久久久不卡| 免费看光身美女| 国产精品影院久久| 亚洲精品国产精品久久久不卡| 午夜免费激情av| 国产成人啪精品午夜网站| 香蕉久久夜色| 此物有八面人人有两片| 亚洲精品美女久久久久99蜜臀| 欧美黑人巨大hd| 又爽又黄无遮挡网站| 国产亚洲精品综合一区在线观看| 色哟哟哟哟哟哟| 一进一出好大好爽视频| 色播亚洲综合网| 蜜桃亚洲精品一区二区三区| 欧美高清成人免费视频www| 日韩欧美 国产精品| 欧美黄色片欧美黄色片| 午夜福利在线观看吧| 亚洲欧美精品综合久久99| 欧美日本视频| 91字幕亚洲| 床上黄色一级片| 毛片女人毛片| 精品人妻一区二区三区麻豆 | 国产一级毛片七仙女欲春2| 日韩精品青青久久久久久| 村上凉子中文字幕在线| 色综合站精品国产| 真人一进一出gif抽搐免费| 日韩欧美国产在线观看| 成年版毛片免费区| 成人三级黄色视频| 日韩欧美免费精品| 免费在线观看成人毛片| 国产伦精品一区二区三区四那| 国内精品美女久久久久久| 久久久久久久久久黄片| 久久精品91无色码中文字幕| 国产单亲对白刺激| 又黄又粗又硬又大视频| 午夜福利在线在线| 国产亚洲欧美在线一区二区| 国产淫片久久久久久久久 | 琪琪午夜伦伦电影理论片6080| 欧美日韩中文字幕国产精品一区二区三区| 一个人免费在线观看的高清视频| 亚洲av一区综合| 欧美在线一区亚洲| 18美女黄网站色大片免费观看| 久久亚洲真实| 69av精品久久久久久| 天天躁日日操中文字幕| 日韩高清综合在线| 禁无遮挡网站| 欧美日韩一级在线毛片| bbb黄色大片| av片东京热男人的天堂| 中出人妻视频一区二区| 亚洲人成网站在线播| 亚洲天堂国产精品一区在线| 怎么达到女性高潮| 真人做人爱边吃奶动态| 久久精品国产综合久久久| 又黄又粗又硬又大视频| 亚洲av五月六月丁香网| 亚洲成人免费电影在线观看| 国产一区二区三区视频了| 91久久精品电影网| 99国产极品粉嫩在线观看| 香蕉av资源在线| 中亚洲国语对白在线视频| 国产久久久一区二区三区| 亚洲av五月六月丁香网| 国产精品av视频在线免费观看| 在线观看美女被高潮喷水网站 | 一本一本综合久久| 亚洲av成人精品一区久久| 午夜福利在线观看免费完整高清在 | 精品一区二区三区视频在线 | 国产av不卡久久| 精品一区二区三区视频在线 | 国产一区二区亚洲精品在线观看| 操出白浆在线播放| 真人做人爱边吃奶动态| 国产精品久久久久久亚洲av鲁大| 成人特级av手机在线观看| bbb黄色大片| 日韩欧美在线乱码| 久久人人精品亚洲av| 俄罗斯特黄特色一大片| 成人鲁丝片一二三区免费| 国产亚洲精品一区二区www| 国产男靠女视频免费网站| 我要搜黄色片| 午夜福利成人在线免费观看| 白带黄色成豆腐渣| 草草在线视频免费看| 久久精品91无色码中文字幕| 最新美女视频免费是黄的| 国产精品久久久久久久电影 | 精品久久久久久久久久免费视频| 免费在线观看影片大全网站| 女生性感内裤真人,穿戴方法视频| 熟女人妻精品中文字幕| 少妇的逼水好多| 人人妻,人人澡人人爽秒播| a级毛片a级免费在线| 狂野欧美白嫩少妇大欣赏| 亚洲成av人片在线播放无| 色在线成人网| 亚洲国产精品成人综合色| 久久久久精品国产欧美久久久| 小说图片视频综合网站| 国产精品一区二区三区四区久久| 精品熟女少妇八av免费久了| 午夜日韩欧美国产| 精品久久久久久,| 久久久久久久精品吃奶| 丰满人妻熟妇乱又伦精品不卡| 亚洲电影在线观看av| 免费高清视频大片| 国产美女午夜福利| 国产伦在线观看视频一区| 国产一区二区在线观看日韩 | 国产成人av激情在线播放| 国产精品综合久久久久久久免费| 亚洲狠狠婷婷综合久久图片| 久久精品国产99精品国产亚洲性色| 99久久九九国产精品国产免费| 午夜免费观看网址| 麻豆国产av国片精品| 成人精品一区二区免费| 午夜福利18| 丝袜美腿在线中文| 两个人看的免费小视频| 一区二区三区高清视频在线| a级毛片a级免费在线| 99久久综合精品五月天人人| 又紧又爽又黄一区二区| 天堂动漫精品| 一级黄片播放器| 久久久久国产精品人妻aⅴ院| 久久精品国产清高在天天线| 成人性生交大片免费视频hd| 国产亚洲精品综合一区在线观看| www日本黄色视频网| 久久久久久久精品吃奶| 一本精品99久久精品77| 国产爱豆传媒在线观看| 日韩欧美免费精品| www.熟女人妻精品国产| 国产一区二区三区视频了| 舔av片在线| 欧美另类亚洲清纯唯美| 国产午夜福利久久久久久| 国产一区二区激情短视频| 成人鲁丝片一二三区免费| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品美女久久久久99蜜臀| aaaaa片日本免费| 岛国在线观看网站| 亚洲,欧美精品.| 精品无人区乱码1区二区| 欧美激情久久久久久爽电影| 久久精品国产清高在天天线| 国产美女午夜福利| 欧美日韩瑟瑟在线播放| 在线免费观看的www视频| 女警被强在线播放| 久久精品国产亚洲av涩爱 | 成人av一区二区三区在线看| 偷拍熟女少妇极品色| 亚洲精品乱码久久久v下载方式 | 亚洲欧美日韩高清专用| 国产亚洲精品一区二区www| 国产精品久久久人人做人人爽| 成人性生交大片免费视频hd| 伊人久久精品亚洲午夜| 国产精品免费一区二区三区在线| 一进一出抽搐动态| 婷婷亚洲欧美| 长腿黑丝高跟| 欧美3d第一页| 欧美日韩综合久久久久久 | 国产伦精品一区二区三区四那| 成人三级黄色视频| 美女cb高潮喷水在线观看| 久久精品国产综合久久久| 国产免费av片在线观看野外av| 一级作爱视频免费观看| 欧美精品啪啪一区二区三区| 黄片小视频在线播放| 99热6这里只有精品| 国产精品1区2区在线观看.| АⅤ资源中文在线天堂| 久久草成人影院| 国内精品美女久久久久久| 性色av乱码一区二区三区2| 黑人欧美特级aaaaaa片| 国产野战对白在线观看| 18禁裸乳无遮挡免费网站照片| 国产亚洲精品综合一区在线观看| 黄色女人牲交| 亚洲人与动物交配视频| 国产黄片美女视频| 麻豆国产97在线/欧美| 色综合欧美亚洲国产小说| 亚洲精品在线美女| 国产黄色小视频在线观看| 黑人欧美特级aaaaaa片| 无遮挡黄片免费观看| 午夜福利在线观看吧| 午夜日韩欧美国产| 精品国产美女av久久久久小说| 少妇人妻精品综合一区二区 | 午夜激情福利司机影院| 欧美最新免费一区二区三区 | 亚洲激情在线av| e午夜精品久久久久久久| 在线天堂最新版资源| 听说在线观看完整版免费高清| 亚洲天堂国产精品一区在线| 国产男靠女视频免费网站| 一区二区三区免费毛片| 小蜜桃在线观看免费完整版高清| 日韩成人在线观看一区二区三区| 亚洲性夜色夜夜综合| 国产麻豆成人av免费视频| 99热精品在线国产| 亚洲成人久久性| 亚洲精品美女久久久久99蜜臀| 青草久久国产| 一本精品99久久精品77| 麻豆成人av在线观看| 久久这里只有精品中国| 岛国在线免费视频观看| 久久精品国产99精品国产亚洲性色| 在线观看舔阴道视频| av中文乱码字幕在线| 淫秽高清视频在线观看| 熟女人妻精品中文字幕| a级一级毛片免费在线观看| 又粗又爽又猛毛片免费看| 综合色av麻豆| 美女 人体艺术 gogo| 久久久久久久久中文| 久久久久亚洲av毛片大全| 亚洲精品粉嫩美女一区| 狂野欧美白嫩少妇大欣赏| 免费av观看视频| 亚洲成人精品中文字幕电影| 一夜夜www| 亚洲人成网站在线播| av天堂在线播放| 好男人在线观看高清免费视频| 婷婷亚洲欧美| 亚洲性夜色夜夜综合| 脱女人内裤的视频| 他把我摸到了高潮在线观看| 国产成人欧美在线观看| 大型黄色视频在线免费观看| 男人的好看免费观看在线视频| 又黄又爽又免费观看的视频| 国产精品自产拍在线观看55亚洲| 国产aⅴ精品一区二区三区波| 亚洲精品日韩av片在线观看 | 国产免费男女视频| 欧美一区二区国产精品久久精品| 九色国产91popny在线| avwww免费| 欧美3d第一页| 91字幕亚洲| 成年人黄色毛片网站| 少妇高潮的动态图| 我要搜黄色片| 小蜜桃在线观看免费完整版高清| 国产成年人精品一区二区| 亚洲精品乱码久久久v下载方式 | 特大巨黑吊av在线直播| 久久精品人妻少妇| 日韩亚洲欧美综合| 欧美最新免费一区二区三区 | 久久久久久国产a免费观看| xxx96com| 黄色女人牲交| 亚洲最大成人手机在线| 国产精品99久久久久久久久| 国产一区二区亚洲精品在线观看| 国产在线精品亚洲第一网站| 亚洲片人在线观看| 国产亚洲精品一区二区www| 久久久久久久久久黄片| 国产精品久久久久久久久免 | 动漫黄色视频在线观看| 久久亚洲真实| 亚洲最大成人中文| 免费大片18禁| 国产亚洲精品综合一区在线观看| 99在线人妻在线中文字幕| 天堂影院成人在线观看| 国产真实伦视频高清在线观看 | 欧美最黄视频在线播放免费| 97碰自拍视频| 欧美激情在线99| 国产成人a区在线观看| 波多野结衣高清无吗| 久久6这里有精品| 观看免费一级毛片| 欧美+亚洲+日韩+国产| 精品免费久久久久久久清纯| 国产av一区在线观看免费| 悠悠久久av| 在线观看美女被高潮喷水网站 | 国产精品久久久久久精品电影| 久久久色成人| netflix在线观看网站| 男女下面进入的视频免费午夜| 国产精品 国内视频| 日本在线视频免费播放| 熟女电影av网| 久久香蕉国产精品| 国产精品亚洲av一区麻豆| h日本视频在线播放| 97超视频在线观看视频| 国产麻豆成人av免费视频| 无限看片的www在线观看| 精品久久久久久,| 国产成人影院久久av| 久久久久久久久中文| 最新中文字幕久久久久| 好男人在线观看高清免费视频| 香蕉久久夜色| 一进一出好大好爽视频| 国产成人av教育| www日本在线高清视频| 日本a在线网址| 欧美日韩国产亚洲二区| 国产精品爽爽va在线观看网站| 波多野结衣巨乳人妻| 欧美一区二区国产精品久久精品| 免费一级毛片在线播放高清视频| 床上黄色一级片| av专区在线播放| 亚洲内射少妇av| 国产av麻豆久久久久久久| 国产成人福利小说| 精品国内亚洲2022精品成人| 日本精品一区二区三区蜜桃| 国产在线精品亚洲第一网站| 深夜精品福利| 久久久久久久亚洲中文字幕 | 亚洲av免费在线观看| 亚洲国产精品久久男人天堂| 搡女人真爽免费视频火全软件 | 久久人妻av系列| 在线免费观看的www视频| 国产高清三级在线| 日韩欧美三级三区| 精品久久久久久久久久免费视频| 欧美黑人欧美精品刺激| 国产色婷婷99| 哪里可以看免费的av片| 美女被艹到高潮喷水动态| 亚洲国产欧美人成| 少妇的逼水好多| 嫩草影院精品99| 成人永久免费在线观看视频| 天美传媒精品一区二区| 757午夜福利合集在线观看| 成人午夜高清在线视频| 亚洲色图av天堂| 国模一区二区三区四区视频| 国产亚洲精品久久久久久毛片| 午夜老司机福利剧场| 久久精品91蜜桃| 日韩亚洲欧美综合| 国产爱豆传媒在线观看| 午夜激情欧美在线| 人人妻人人澡欧美一区二区| 久久精品国产99精品国产亚洲性色| 丁香欧美五月| 亚洲国产中文字幕在线视频| 久久亚洲精品不卡| 中文在线观看免费www的网站| 欧美黄色片欧美黄色片| 精品国内亚洲2022精品成人| 久久久国产精品麻豆| 好看av亚洲va欧美ⅴa在| 欧美一区二区亚洲| 成人高潮视频无遮挡免费网站| 最近最新中文字幕大全免费视频| 白带黄色成豆腐渣| 日韩高清综合在线| 日本黄大片高清| 黄色片一级片一级黄色片| 国产男靠女视频免费网站| 亚洲精品456在线播放app | 国产成+人综合+亚洲专区| 狠狠狠狠99中文字幕| 激情在线观看视频在线高清| 亚洲精品影视一区二区三区av| 亚洲人成网站在线播放欧美日韩| 免费观看的影片在线观看| 成人av一区二区三区在线看| x7x7x7水蜜桃| 欧美av亚洲av综合av国产av| 精品99又大又爽又粗少妇毛片 | www日本在线高清视频| 亚洲精华国产精华精| 精品一区二区三区视频在线观看免费| 1024手机看黄色片| 国产精品影院久久| 亚洲第一欧美日韩一区二区三区| 亚洲美女视频黄频| 亚洲av成人精品一区久久| 成年女人永久免费观看视频| 一本综合久久免费| 免费搜索国产男女视频| 一区二区三区激情视频| 操出白浆在线播放| 无遮挡黄片免费观看| 日本三级黄在线观看| 欧美黑人欧美精品刺激| 中文字幕av在线有码专区| 国产黄色小视频在线观看| 婷婷精品国产亚洲av| 午夜精品在线福利| 久久久久亚洲av毛片大全| 特大巨黑吊av在线直播| 高清日韩中文字幕在线| 夜夜躁狠狠躁天天躁| 久久久国产精品麻豆| 天天一区二区日本电影三级| 成人无遮挡网站| 日韩欧美 国产精品| 国产69精品久久久久777片| 午夜a级毛片| 天堂网av新在线| 亚洲性夜色夜夜综合| 婷婷亚洲欧美| 亚洲av免费高清在线观看| 人人妻,人人澡人人爽秒播| 午夜福利成人在线免费观看| 国产精品免费一区二区三区在线| 麻豆一二三区av精品| 亚洲七黄色美女视频| 亚洲中文字幕日韩| 性色av乱码一区二区三区2| 免费无遮挡裸体视频| 久久久色成人| 男女午夜视频在线观看| 国产精品 欧美亚洲| 国产av不卡久久| 国产精品亚洲美女久久久| 色精品久久人妻99蜜桃| 一区二区三区高清视频在线| 国产精品亚洲av一区麻豆| 香蕉av资源在线| 麻豆一二三区av精品| 国产高清有码在线观看视频| 国产精品亚洲一级av第二区| 日韩有码中文字幕| 五月伊人婷婷丁香| 精品久久久久久久久久免费视频| 最近视频中文字幕2019在线8|