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

    利用區(qū)域地震波形振幅包絡(luò)約束朝鮮地下核試驗(yàn)的埋深和當(dāng)量

    2016-07-28 07:04:27林鑫姚振興
    地球物理學(xué)報 2016年6期

    林鑫, 姚振興

    1 華中科技大學(xué)地球物理研究所,華中科技大學(xué)物理學(xué)院, 武漢 430074 2 中國科學(xué)院地質(zhì)與地球物理研究所, 地球與行星物理重點(diǎn)實(shí)驗(yàn)室, 北京 100029

    ?

    利用區(qū)域地震波形振幅包絡(luò)約束朝鮮地下核試驗(yàn)的埋深和當(dāng)量

    林鑫1, 2, 姚振興2

    1 華中科技大學(xué)地球物理研究所,華中科技大學(xué)物理學(xué)院, 武漢4300742 中國科學(xué)院地質(zhì)與地球物理研究所, 地球與行星物理重點(diǎn)實(shí)驗(yàn)室, 北京100029

    摘要核爆當(dāng)量和埋藏深度是地下核試驗(yàn)的兩個重要參數(shù).根據(jù)中國東北地區(qū)區(qū)域范圍內(nèi)地震臺站觀測記錄,利用Pn,Pg,Sn和Lg波波形,計算水平分量的尾波振幅包絡(luò),調(diào)查2006年10月9日,2009年5月25日,2013年2月13日和2016年1月6日四次朝鮮地下核試驗(yàn)的爆炸當(dāng)量和埋藏深度.以牡丹江(MDJ)臺站記錄為例,對兩個水平分量波形進(jìn)行帶通濾波,計算平均波形的振幅包絡(luò).最終得到區(qū)域地震臺站水平分量振幅包絡(luò),振幅穩(wěn)定,包絡(luò)振幅的變化清晰地顯示區(qū)域震相的位置.區(qū)域震相的時間域振幅包絡(luò)由震源譜函數(shù)、傳播效應(yīng)、臺基響應(yīng)和傳遞函數(shù)及尾波形狀函數(shù)構(gòu)成.通過網(wǎng)格搜索的方法,擬合水平分量記錄的波形包絡(luò),可以獲得核爆當(dāng)量和埋藏深度的估計.結(jié)果表明,朝鮮四次地下核試驗(yàn)爆炸當(dāng)量以時間為序從0.6±0.2 kt到3.0±1.5 kt,增加至10.0±2.0 kt,再降到8.0±2.0 kt.2006年核爆的埋藏深度較淺,為150±100 m,2009年朝鮮核試驗(yàn)的埋藏深度為350±100 m,2013年和2016年朝鮮核試驗(yàn)的深度非常一致,均為500±200 m.這些結(jié)果與前人的調(diào)查結(jié)果一致性較高,說明使用單一地震臺站時間域水平分量尾波振幅包絡(luò)是可能同時約束地下核試驗(yàn)爆炸當(dāng)量和埋藏深度的.

    關(guān)鍵詞區(qū)域震相; 振幅包絡(luò); 朝鮮核試驗(yàn); 爆炸當(dāng)量; 埋藏深度

    1引言

    從2006年10月9日朝鮮第一次地下核試驗(yàn)以來,分別于2009年5月25日,2013年2月13日和2016年1月6日相繼進(jìn)行了三次規(guī)模較大的核爆試驗(yàn)(以下簡稱 NKT1, 2,3和4).四次朝鮮核爆的試驗(yàn)地點(diǎn)非常接近,在中國東北的地震記錄中能夠清楚地識別出區(qū)域震相中的Pn,Pg,Lg和Rayleigh波(例如, Zhao et al., 2008,2012,2014).由于NKT2和NKT3的震級較大,地震記錄的信噪比相對較高,甚至在跨過日本海的日本F-NET臺網(wǎng)的臺站也記錄到清楚的P波震相.朝鮮核試驗(yàn)所產(chǎn)生的地震數(shù)據(jù),因?yàn)檎鹪春唵?,記錄質(zhì)量高,能夠?yàn)楣浪惚ó?dāng)量、分析朝鮮半島核爆事件與地震事件之間的震源差異提供重要的信息,還可以用于研究核爆記錄覆蓋區(qū)域的地殼和上地幔的地震速度和衰減結(jié)構(gòu).因此,在全世界范圍內(nèi)激發(fā)了許多研究者的極大興趣(Kim and Richards,2007;Kvrna et al.,2007;Salzberg and Marshall,2007;Bonner et al.,2008;Koper et al.,2008;Hong et al.,2008;Patton and Taylor,2008;Salzberg,2008;Tibuleac et al.,2008;Zhao et al.,2008;Hong and Rhie,2009;Kim et al.,2009;Chi et al.,2010;Murphy et al.,2010, 2011;Ni et al.,2010;Schlittenhardt et al.,2010;Shin et al.,2010;Rougier等,2011;Zhang and Wen,2013).

    朝鮮核爆發(fā)生以后,很多研究者給出了地震學(xué)當(dāng)量的估計.但是,因?yàn)槭褂觅Y料和方法不同,當(dāng)量估計結(jié)果差異較大.例如,對于2006年10月9日朝鮮地下核試驗(yàn),Salzberg和Marshall(2007)利用同一臺站不同事件地震記錄的能量包絡(luò)比值法,參照中國東北地區(qū)用于地震測深的已知當(dāng)量的化爆事件,獲得NKT1的當(dāng)量估計為0.45 kt.基于最小埋藏深度假定,Zhao等(2008)利用中國東北及鄰近地區(qū)的地震Lg波資料得到的當(dāng)量估計是0.48 kt.但是,Chun等(2011)使用中朝邊界附近的流動臺網(wǎng)記錄的區(qū)域Lg波資料得到的結(jié)果是1.22 kt.Murphy等(2013)利用遠(yuǎn)震P波振幅譜得到的核爆當(dāng)量與埋深分別為0.9 kt和200 m.

    早期核試驗(yàn)的核爆當(dāng)量較高,通常是計算遠(yuǎn)震P波震級,利用震級與當(dāng)量之間的經(jīng)驗(yàn)關(guān)系獲得爆炸當(dāng)量的估計(例如, Ringdal et al.,1992; Douglas and Marshall,1996; Murphy and Barker, 2001).也可以通過遠(yuǎn)震P波波形反演求解核爆等效震源參數(shù),利用相對靜力學(xué)強(qiáng)度給出核試驗(yàn)的當(dāng)量估計(趙連鋒等,2005).隨著核武器小型化,地下核試驗(yàn)的當(dāng)量逐漸減小,區(qū)域地震記錄成為研究核爆當(dāng)量的首選資料.例如,Vergino和 Mensing(1990)發(fā)展了使用區(qū)域震相Pn波計算體波震級的方法;Nuttli (1973,1986a,1986b)以Lg波列中第三峰值振幅為基礎(chǔ)發(fā)展了一種區(qū)域地震級系統(tǒng),并建立了Lg波震級與爆炸當(dāng)量之間的關(guān)系.Patton(1988)將Nuttli的方法推廣應(yīng)用到Nevada測試點(diǎn)(NTS)的更大的數(shù)據(jù)體.Hansen等(1990)和Ringdal等(1992)發(fā)現(xiàn)基于Lg波均方根(rms)振幅的震級估計更加穩(wěn)定,即使臺站很少,也能給出令人滿意的當(dāng)量估計.

    最近,一些研究者嘗試使用區(qū)域地震波形的高頻波列,建立了計算震級和估計爆炸當(dāng)量的尾波法(Mayeda and Walter,1996;Mayeda et al.,2003;Mayeda,1993;Murphy et al.,2009).尾波是地方震記錄圖中體波或面波之后的某些震蕩,通常被認(rèn)為是由于分布于地殼和上地幔的橫向不均勻性引起的散射波,尾波的激發(fā)依賴于臺站的地質(zhì)情況(Aki,1969;Aki and Chouet,1975).在地方性地震范圍內(nèi),體波震相Pn,Pg,Sn和短周期面波Lg較容易識別,但波形較為復(fù)雜且容易受到源和路徑不均勻性的影響,但其尾波對輻射花樣、傳播路徑等不敏感,波形特征穩(wěn)定.因此,地方性地震范圍內(nèi)單個地震臺站體波震相的尾波有可能用于測量地下核試驗(yàn)的當(dāng)量和埋深.Pasyanos等(2012)將這種使用單一震相尾波方法擴(kuò)展到多個震相,以等效震源函數(shù)、傳播效應(yīng)、臺基響應(yīng)和傳遞函數(shù)的物理模型為基礎(chǔ),建立了一種通過單一地震臺站擬合區(qū)域波形全波列的包絡(luò),同時獲得地下核試驗(yàn)當(dāng)量和埋藏深度的方法.

    2數(shù)據(jù)

    選取中國東北地區(qū)距離朝鮮核試驗(yàn)場(表1)震中距范圍在200 km到500 km之間的11個地震臺站以及2002年到2013年間該區(qū)域發(fā)生的27個天然地震記錄(ML≥2.5)來進(jìn)行數(shù)據(jù)分析(表2).

    表1 朝鮮地下核試驗(yàn)事件目錄

    注:數(shù)據(jù)源1. USGS (Wen and Long, 2010; Zhao et al., 2014).

    2. 衛(wèi)星圖像(Wen and Long, 2010; Zhao et al., 2014).

    3方法

    基于尾波單一散射模型(Aki,1969; Aki and Chouet,1975),Mayeda等(2003)給出在中心頻率f和震中距r的單個區(qū)域震相尾波包絡(luò)(圖2)計算公式:Ac(t,f,r)=S(f)·G(r)·B(f,r)·P(f)

    圖1 研究區(qū)域背景圖. 朝鮮地下核試驗(yàn)點(diǎn)(五角星),天然地震記錄(黑色十字形),區(qū)域臺站(黑色三角形)Fig.1 Map showing the background of study area and North Korean underground nuclear tests (star), earthquakes (black cross), and stations (black triangles)

    ·T(f)·E(t,f,r),

    (1)

    其中S(f)為震源譜,G(r)為幾何擴(kuò)展項(xiàng),B(f,r)為衰減項(xiàng),P(f)為臺基響應(yīng),T(f)為傳遞函數(shù),E(t,f,r)為描述尾波形狀的函數(shù),

    E(t,f,r)=H(t-ta)·(t-ta)-γ(f,r)

    ·exp[b(f,r)·(t-ta)],

    (2)

    H為單位階躍函數(shù);ta為到時;γ(f,r)和b(f,r)為尾波形狀參數(shù).

    圖2 單個震相理論包絡(luò)形狀Fig.2 The shape of single phase theory envelope

    地方性或區(qū)域地震范圍地震全記錄的包絡(luò)(圖3a)為該區(qū)域噪聲水平Anoise與所有區(qū)域震相包絡(luò)之和(Pasyanos et al.,2012):

    鬼子軍官又說一聲有戲。然后做了個讓他們不要怕的手勢,說了一大通生硬的中國話。大意是,你們幾個都是支那的精英,我們十分歡迎各位和大日本皇軍合作,如果同意,金票大大的有,花姑娘大大的有。

    (3)

    應(yīng)用于實(shí)際觀測資料,首先將兩個水平分量記錄經(jīng)過反褶積去掉儀器響應(yīng),獲得地面速度記錄,然后,利用八階零相(四個極點(diǎn),雙通)巴特沃茲濾波器得到窄頻帶波形,應(yīng)用窄帶包絡(luò)公式(f為中心頻率,v(t)為一個改正后的帶通濾波地震圖地面速度,h(t)為v(t)對應(yīng)的希爾伯特變換):

    (4)

    求出兩水平分量的窄帶包絡(luò),再求取兩個水平包絡(luò)的平均值,經(jīng)過平滑處理,最后獲得該頻帶的尾波包絡(luò)(圖3b).這個過程也能用單分量(如垂直分量)來完成,但比起單分量,兩個水平分量能提供一個更平穩(wěn)光滑的包絡(luò).對于天然地震事件和地下核試驗(yàn)事件,其尾波包絡(luò)(1)式中的各參數(shù)項(xiàng)只有震源項(xiàng)不同,其他參數(shù)項(xiàng)都可以由天然地震記錄資料獲得.應(yīng)用于地下核試驗(yàn)事件中,通過天然地震記錄構(gòu)建地下核爆的理論包絡(luò),與實(shí)際結(jié)果比較,可以給出地下核試驗(yàn)的爆炸當(dāng)量和埋藏深度(Pasyanos et al.,2012).

    我們用事件前噪聲水平來經(jīng)驗(yàn)地估計Anoise,取Pn到達(dá)之前的噪聲窗,進(jìn)行統(tǒng)計處理,得到均值.對于殼幔模型,地殼用單層均勻介質(zhì)層表示,殼幔邊界水平,忽略坡度變化,地幔水平方向分塊,每個方塊內(nèi)速度均勻,射線在地幔頂部走直線,忽略速度橫向變化引起的實(shí)際射線彎曲,忽略地球曲率和地幔頂部垂向速度梯度的影響(裴順平,2002).由于所用數(shù)據(jù)中天然地震和四次核試驗(yàn)事件的震源非常淺,震中距離較近,我們將Pg和Lg波的傳播路徑近似為震源到臺站的直線(圖4).Pn和Sn波的走時方程可描述為:

    圖3 區(qū)域事件理想包絡(luò)形狀(a)和實(shí)際波形包絡(luò)形狀(b).(b)為MDJ臺站記錄的2002年4月16日天然地震的不同頻帶的包絡(luò),橫軸為時間(取發(fā)震時刻后250 s的記錄),縱軸為速度振幅的對數(shù).Fig.3 (a) Expected envelope shapes for regional events; (b) Shapes of real waveforms envelopes in different frequency bands of the earthquake event on 16 April 2002 recorded at MDJ station. X-axes are times (250 s after the origin time). Y-axes are the log amplitude of velocity.

    (5)

    (6)

    其中Hs為震源處地殼厚度,hs為震源深度,Hr為臺站處地殼厚度,hr為臺站高程,r為震中距.

    圖4 區(qū)域震相射線路徑圖及殼幔模型Fig.4 Ray paths of regional phases and the crust-mantle model

    中國東北地區(qū)Pn波平均波速度約為8.0~8.05 km·s-1,單層水平地殼P波平均速度約為6.3 km·s-1(賈軍等,2011;裴順平等,2004a;張國民等,2012),Sn波平均速度約為4.55 km·s-1,地殼S波速度約為3.6 km·s-1(裴順平等,2004b), Pg波平均速度約為6.0 km·s-1,Lg波群速度為3.5 km·s-1(趙連鋒,2008).

    (7)

    (8)

    (9)

    采用Mueller-Murphy模型(Mueller and Murphy,1971)來描述爆炸震源的震源譜,爆炸源P波部分震源譜的表達(dá)式為:

    (10)

    尾波形狀參數(shù)b(f,r)和γ(f,r)代表尾波振幅的衰減,分別依賴于介質(zhì)內(nèi)在吸收引起的尾波幾何擴(kuò)展和衰減速率,b(f,r)控制尾波包絡(luò)后半部分的斜率,γ(f,r)控制著緊跟直達(dá)波之后的尾波初期的衰減,γ(f,r)值越大表明尾波初期更陡峭或是有更大的直達(dá)峰值.根據(jù)Mayeda等(2003)給出的在地方性及區(qū)域地震范圍內(nèi)b(f,r)和γ(f,r)的函數(shù)關(guān)系:

    (11)

    (12)

    其中γ0,γ1,γ2和b0,b1,b2為常量,可用網(wǎng)格搜索法確定.我們利用區(qū)域內(nèi)27個天然地震的地震記錄資料,在選定的頻帶內(nèi)為每個地震事件-臺站對用網(wǎng)格搜索找出γ的最佳參數(shù)值,γ的變化范圍從0到1.6每次增加0.1.當(dāng)γ固定之后用最小二乘法來反演求解b(Yoo et al.,2011).通過匹配E(t,f,r)與實(shí)際波形的衰減速率,由S波尾波和Lg波尾波校準(zhǔn)得到我們所用區(qū)域臺站的尾波衰減參數(shù)b(f,r)和γ(f,r)的數(shù)值(圖5).

    表3 MDJ臺站各震相的Q0值和η值

    尾波包絡(luò)(1)式中幾何擴(kuò)展、衰減和尾波形狀函數(shù)項(xiàng)確定后,將S(f)歸一化,兩邊取對數(shù)log10,應(yīng)用L1范數(shù)改正并通過網(wǎng)格搜索比較MDJ臺站記錄的2002年4月16日和2004年12月16日兩次天然地震事件來得到MDJ臺站臺基響應(yīng)和傳遞函數(shù),即對理論包絡(luò)做出改正,再應(yīng)用于地下核試驗(yàn)事件中(圖6).

    4結(jié)果

    圖5 MDJ臺站尾波衰減參數(shù)γ(f,r)和b(f,r)在頻帶2~4 Hz內(nèi)所得結(jié)果Fig.5 Coda-decay parameters γ(f,r) and b(f,r) in the 2~4 Hz frequency band at MDJ station

    圖6 MDJ臺站理論包絡(luò)改正2002年(a)和2004年(b)地震事件MDJ臺站的擬合結(jié)果(2~4 Hz頻帶內(nèi)),實(shí)際波形包絡(luò)(黑色波形),未改正的地震理論包絡(luò)(藍(lán)色虛線),改正后的地震理論包絡(luò)(綠色實(shí)線).橫軸為時間(取發(fā)震時刻后250s的記錄),縱軸為速度振幅的對數(shù).Fig.6 Corrected theoretical envelopes at MDJ station(a) Seismic event of 2002 (2~4 Hz); (b) Seismic event of 2004 (2~4 Hz) with envelopes of real data (black waveforms), uncorrected theoretical seismic envelopes (blue dash lines), corrected theoretical envelopes (green lines). X-axes are times (250 s after the origin time). Y-axes are log amplitudes of velocity.

    圖7 MDJ臺站四次朝鮮地下核試驗(yàn)事件垂直分量速度記錄波形上標(biāo)出了視群速度值.Fig.7 Vertical velocities of four North Korean underground nuclear tests recorded at MDJ station Apparent group velocity values are marked on waveforms.

    MDJ臺站記錄到朝鮮四次地下核試驗(yàn)的垂直分量速度記錄如圖7所示.對于朝鮮的四次地下核試驗(yàn)事件,根據(jù)中國東北地區(qū)地震臺網(wǎng)的區(qū)域地震臺站全波列尾波的理論包絡(luò),采用爆炸源模型,介質(zhì)類型為花崗巖,通過網(wǎng)格搜索比較不同當(dāng)量、不同埋深事件的理論包絡(luò)與實(shí)際結(jié)果,確定四次地下核試驗(yàn)事件的最佳擬合理論包絡(luò)(圖8,9).從圖8中可以看出,尾波包絡(luò)的振幅對于當(dāng)量和埋深的變化是比較敏感的.當(dāng)量1 kt、埋深100 m的理論包絡(luò)(藍(lán)色實(shí)線)比2006年核試驗(yàn)事件實(shí)際包絡(luò)的當(dāng)量和埋深要大,但比2009年核試驗(yàn)事件實(shí)際包絡(luò)的當(dāng)量和埋深要小得多.對于2013年和2016年兩次核試驗(yàn)事件,當(dāng)量15 kt、埋深500 m的理論包絡(luò)(藍(lán)色實(shí)線)大于實(shí)際包絡(luò)的當(dāng)量和埋深,而當(dāng)量5 kt、埋深500 m的理論包絡(luò)(黑色虛線)小于實(shí)際包絡(luò)的當(dāng)量和埋深,對比兩次核試驗(yàn)事件實(shí)際波形包絡(luò)可以看出兩次核試驗(yàn)事件的實(shí)際當(dāng)量和埋深差別較小.我們將當(dāng)量和深度分別從0.01 kt到100 kt、10 m到1000 m的范圍內(nèi)進(jìn)行搜索,找出最佳擬合范圍(圖9)獲得2006年核試驗(yàn)事件的當(dāng)量為0.6 kt(0.4~0.8 kt),埋藏深度為150 m(50~250 m),2009年核試驗(yàn)事件的當(dāng)量為3 kt(1.5~4.5 kt),埋藏深度為350 m(250~450 m),2013年核試驗(yàn)事件的當(dāng)量為10 kt(8~12 kt),埋藏深度為500 m(300~700 m),而2016年核試驗(yàn)事件的當(dāng)量為8 kt(6~10 kt),埋藏深度為500 m(300~700 m).

    我們選用除MDJ臺站外其他10個區(qū)域地震臺站采用爆炸源模型來研究2013年地下核試驗(yàn)事件,結(jié)果如圖10所示.從圖10中可以看出,當(dāng)量10 kt、埋藏深度500 m時,10個區(qū)域臺站尾波的理論包絡(luò)與實(shí)際資料都擬合較好,表明整個區(qū)域地震臺網(wǎng)內(nèi)的臺站結(jié)果具有高度的相似性,說明地方性及區(qū)域地震范圍內(nèi)單個地震臺站體波震相的尾波有可能用于測量地下核試驗(yàn)的當(dāng)量和埋深.

    Chi等(2010),Murphy等(2010),Salzberg(2008),Zhang和Wen(2013),Pasyanos 等(2012)及Zhao等(2008,2012,2014,2015)都曾給出朝鮮地下核試驗(yàn)核爆當(dāng)量和深度的研究結(jié)果.表4給出了本文與前人調(diào)查結(jié)果的比較,可以看出本文結(jié)果與表中幾位國內(nèi)外專家的研究一致性較高,說明地方性及區(qū)域地震范圍內(nèi)單個地震臺站體波震相的尾波是可能同時約束地下核試驗(yàn)的爆炸當(dāng)量和埋藏深度的.

    5結(jié)論

    尾波對輻射花樣、傳播路徑等不敏感,且波形特征穩(wěn)定,本文通過構(gòu)建一個理論的區(qū)域震相尾波振幅包絡(luò),來擬合地方性及區(qū)域地震范圍內(nèi)單一地震臺站實(shí)際波形記錄,調(diào)查朝鮮四次地下核試驗(yàn)的爆炸當(dāng)量和埋藏深度.所有區(qū)域震相的時間域振幅包絡(luò)由其震源譜函數(shù)、傳播特征、臺基響應(yīng)、傳遞函數(shù)及尾波形狀函數(shù)組成,我們通過研究區(qū)域內(nèi)天然地震記錄資料對各參數(shù)項(xiàng)進(jìn)行實(shí)際分析,得到計算尾波理論包絡(luò)所需各參數(shù)值,例如27個天然地震記錄資料來獲得MDJ臺站所選用頻帶內(nèi)的尾波形狀參數(shù)b(f,r)和γ(f,r)值.我們對當(dāng)量和埋藏深度進(jìn)行網(wǎng)格搜索來擬合臺站水平分量記錄的尾波包絡(luò),得到2006年核試驗(yàn)事件的當(dāng)量為0.6 kt±0.2 kt,埋藏深度為150 m±100 m,2009年核試驗(yàn)事件的當(dāng)量為3 kt±1.5 kt,埋藏深度為350 m±100 m,2013年核試驗(yàn)事件的當(dāng)量為10 kt±2 kt,埋藏深度范圍為500 m±200 m,而2016年核試驗(yàn)事件的當(dāng)量為8 kt±2 kt,埋藏深度范圍為500 m±200 m.

    圖8 MDJ臺站四次朝鮮地下核試驗(yàn)事件擬合結(jié)果(2~4 Hz頻帶內(nèi))實(shí)際波形包絡(luò)(黑色波形),最佳擬合理論包絡(luò)(紅色實(shí)線). 2006年(a)和2009年(b)圖中藍(lán)色實(shí)線為當(dāng)量1 kt、埋深100 m的理論包絡(luò),2013年(c)和2016年(d)圖中藍(lán)色實(shí)線和黑色虛線分別為當(dāng)量15 kt、埋深500 m及當(dāng)量5 kt、埋深500 m的理論包絡(luò).橫軸為時間(取發(fā)震時刻后250 s的記錄),縱軸為速度振幅的對數(shù).Fig.8 Results of four North Korean underground nuclear tests in 2006 (a), 2009 (b), 2013 (c) and 2016 (d) at MDJ station (2~4 Hz)The envelopes of real data (black waveforms), best fitting theoretical envelopes (red lines) are plotted. Blue lines in 2006 (a) and 2009 (b) represent the theoretical envelopes of 1 kt with 100 m depth, while for 2013 (c) and 2016 (d) the blue lines and black dash lines represent the theory envelopes of 15 kt with 500 m and 5 kt with 500 m depths, respectively. X-axes show times (250 s after the origin time). Y-axes show the log amplitudes of velocity.

    2006-05-252009-10-092013-02-122016-01-06當(dāng)量(kt)深度(m)當(dāng)量(kt)深度(m)當(dāng)量(kt)深度(m)當(dāng)量(kt)深度(m)Chi0.4(0.2~0.8)-3.2(1.6~6.4)-----Murphy0.92004.6550----Salzberg0.85±0.1-------Zhang--7.0±1.961012.2±3.843011.13±4.0*-Pasyanos0.5±0.3100(20~300)2(1~5)200(70~600)----Zhao0.48-2.35-7.47-3.58**-本文0.6±0.2150±1003.0±1.5350±10010.0±2.0500±2008.0±2.0500±200

    注:*中國科學(xué)院大學(xué)精確測定朝鮮核爆(http:∥news.cri.cn/201617/14d96c7a-fe28-4307-fcee-d4b4e2ed8a5a.html).**研究快報:2016年1月6日朝鮮核試驗(yàn)定位、識別和當(dāng)量估計的初步調(diào)查結(jié)果 ( http:∥www.igg.cas.cn/xwzx/kyjz/201601/t20160107_4515052.html).

    圖9 2006年(a)、2009年(b)、2013年(c)和2016(d)年四次朝鮮地下核試驗(yàn)事件的MDJ臺站網(wǎng)格搜索結(jié)果黑色五角星為最佳擬合,色標(biāo)范圍表示不擬合度由最佳擬合處向四周延伸越來越大.Fig.9 Results of grid search for 2006 (a), 2009 (b), 2013 (c) and 2016 (d) North Korean underground nuclear tests at MDJ stationThe best fit is indicated by the star. The color scale from the best fit extending all around means misfit getting greater.

    由于核爆當(dāng)量與埋藏深度依賴于爆炸源模型和傳遞函數(shù),而傳遞函數(shù)如前所述是由該區(qū)域內(nèi)天然地震記錄資料中經(jīng)驗(yàn)地推導(dǎo)出的,需要在未來能更加透徹地理解爆炸源模型和傳遞函數(shù),并希望能給出更好、更合理的解釋.另外,區(qū)域范圍內(nèi)四個區(qū)域震相Pn,Pg,Sn,Lg在實(shí)際波形中能較清楚地觀測到,而在地方震距離內(nèi)(Δ<100 km)的地震臺站上只有Pg和Lg兩個震相,所以我們下一步的工作是研究從近場到區(qū)域地震范圍內(nèi),該方法在更多更密集的臺站上的具體表現(xiàn).

    致謝感謝美國地震學(xué)研究機(jī)構(gòu)聯(lián)合會(IRIS)提供的四次朝鮮地下核試驗(yàn)的數(shù)據(jù)資料,感謝趙連鋒副研究員、郝金來副研究員和趙旭副研究員提供的寶貴建議和幫助.

    圖10 2013年朝鮮地下核試驗(yàn)事件10個區(qū)域臺站擬合結(jié)果(2~4 Hz頻帶內(nèi))實(shí)際波形包絡(luò)(黑色波形),最佳擬合理論包絡(luò)(紅色實(shí)線).橫軸為時間(取發(fā)震時刻后250 s的記錄),縱軸為速度振幅的對數(shù).Fig.10 Fitting results of North Korean underground nuclear test in 2013 at 10 regional stations (2~4 Hz)Envelopes of real data (black waveforms), best fitting theoretical envelopes (red lines) are displayed. X-axes show the time (250 s after the origin time). Y-axes show log amplitudes of velocity.

    References

    Aki K. 1969. Analysis of the seismic coda of local earthquakes as scattered waves.J.Geophys.Res., 74(2): 615-631.

    Aki K, Chouet B. 1975. Origin of coda waves: source, attenuation, and scattering effects.J.Geophys.Res., 80(23): 3322-3342.

    Bonner J, Herrmann R B, Harkrider D, et al. 2008. The surface wave magnitude for the 9 October 2006 North Korean nuclear explosion.Bull.Seismol.Soc.Am., 98(5): 2498-2506.

    Brune J N. 1970. Tectonic stress and the spectra of seismic shear waves from earthquakes.J.Geophys.Res., 75(26): 4997-5009.Chi H C, Park J H, Kim G Y, et al. 2010. Monitoring North Korea nuclear tests: Comparison of 1st and 2nd tests.GeophysicsandGeophysicalExploration, 13(3): 243-248.

    Chun K Y, Wu Y, Henderson G A. 2011. Magnitude estimation and source discrimination: A close look at the 2006 and 2009 North Korean underground nuclear explosions.Bull.Seismol.Soc.Am., 101(3): 1315-1329.

    Douglas A, Marshall P D. 1996. Seismic source size and yield for nuclear explosions. ∥ Monitoring a Comprehensive Test Ban Treaty. Netherlands: Springer, 309-353. Fisk M D. 2006. Source spectral modeling of regional P/S discriminants at nuclear test sites in China and the former Soviet Union.Bull.Seismol.Soc.Am., 96(6): 2348-2367.Hansen R A, Ringdal F, Richards P G. 1990. The stability of RMS Lg measurements and their potential for accurate estimation of the yields of Soviet underground nuclear explosions.Bull.Seismol.Soc.Am., 80(6B): 2106-2126.

    Hong T K, Baag C E, Choi H, et al. 2008. Regional seismic observations of the 9 October 2006 underground nuclear explosion in North Korea and the influence of crustal structure on regional phases.J.Geophys.Res.:SolidEarth(1978—2012), 113(B3),doi: 10.1029/2007JB004950.Hong T K, Rhie J. 2009. Regional source scaling of the 9 October 2006 underground nuclear explosion in North Korea.Bull.Seismol.Soc.Am., 99(4): 2523-2540.

    Jia J, Zhou S Y, Meng X S, et al. 2011. Study on Pn-wave velocity structure in the northeast China.JournalofDisasterPreventionandReduction(in Chinese), 27(2): 53-59.Kim T S, Kang I B, Kim G Y. 2009. Yield ratio estimates using regional Pnand Pgfrom North Korea′s underground nuclear explosions.Geophys.Res.Lett., 36(22), doi:10.1029/2009GL040495.Kim W Y, Richards P G. 2007. North Korean nuclear test: Seismic discrimination low yield.Eos, 88(14): 158-161.

    Koper K D, Herrmann R B, Benz H M. 2008. Overview of open seismic data from the North Korean event of 9 October 2006.SeismologicalResearchLetters, 79(2): 178-185.

    Mayeda K. 1993. mb(LgCoda): A stable single station estimator of magnitude.Bull.Seismol.Soc.Am., 83(3): 851-861.

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

    Mayeda K, Hofstetter A, O′Boyle J L, et al. 2003. Stable and transportable regional magnitudes based on coda-derived moment rate spectra.Bull.Seismol.Soc.Am., 93(1): 224-239. Mueller R A, Murphy J R. 1971. Seismic characteristics of underground nuclear detonations, Part I: Seismic-spectrum scaling.Bull.Seismol.Soc.Am., 61(6): 1675-1692.

    Murphy J R, Barker B W. 2001. Application of network-averaged teleseismic P-wave spectra to seismic yield estimation of underground nuclear explosions.PureandAppliedGeophysics, 158(11): 2123-2171. Murphy J R, Kohl B C, Stevens J L, et al. 2010. Exploitation of the IMS and other data for a comprehensive, advanced analysis of the North Korean nuclear tests. ∥ Proceedings of the 2010 Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies 1. 456-465.

    Murphy J R, Kohl B C, Stevens J L, et al. 2011. Supplemental analysis of the seismic characteristics of the 2006 and 2009 North Korean nuclear tests. ∥ Proceedings of the 2010 Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies 3. 513-523. Murphy K R, Mayeda K, Walter W R. 2009. Lg-coda methods applied to Nevada Test Site events: Spectral peaking and yield estimation.Bull.Seismol.Soc.Am., 99(1): 441-448.

    Murphy J R, Stevens J L, Kohl B C, et al. 2013. Advanced seismic analyses of the source characteristics of the 2006 and 2009 North Korean nuclear tests.Bull.Seismol.Soc.Am., 103(3): 1640-1661.

    Ni S D, Helmberger D, Pitarka A. 2010. Rapid source estimation from global calibrated paths.SeismologicalResearchLetters, 81(3): 498-504.

    Nuttli O W. 1973. Seismic wave attenuation and magnitude relations for eastern North America.J.Geophys.Res., 78(5): 876-885. Nuttli O W. 1986a. Yield estimates of Nevada test site explosions obtained from seismic Lg waves.J.Geophys.Res.:SolidEarth(1978—2012), 91(B2): 2137-2151.Nuttli O W. 1986b. Lg magnitudes of selected East Kazakhstan underground explosions.Bull.Seismol.Soc.Am., 76(5): 1241-1251.Pasyanos M E, Walter W R, Mayeda K M. 2012. Exploiting regional amplitude envelopes: a case study for earthquakes and explosions in the Korean Peninsula.Bull.Seismol.Soc.Am., 102(5): 1938-1948.

    Patton H J. 1988. Application of Nuttli′s method to estimate yield of Nevada Test Site explosions recorded on Lawrence Livermore National Laboratory′s digital seismic system.Bull.Seismol.Soc.Am., 78(5): 1759-1772.

    Patton H J, Taylor S R. 2008. Effects of shock-induced tensile failure on mb-Ms discrimination: Contrasts between historic nuclear explosions and the North Korean test of 9 October 2006.Geophys.Res.Lett., 35(14),doi: 10.1029/2008GL034211.Pei S P. 2002. Pn & Sn tomography of uppermost mantle in China [Ph. D. thesis] (in Chinese). Beijing: Institute of Geophysics, China Earthquake Administration.

    Pei S P, Xu Z H, Wang S Y. 2004a. Discussion on origin of Pn velocity variation in China and adjacent region.ActaSeismologicaSinica(in Chinese), 6(1): 1-10.Pei S P, Xu Z H, Wang S Y. 2004b. Sn wave tomography in the uppermost mantle beneath the China continent and adjacent regions.ChineseJ.Geophys. (in Chinese), 47(2): 250-256.

    Ringdal F, Marshall P D, Alewine R W. 1992. Seismic yield determination of Soviet underground nuclear explosions at the Shagan River test site.Geophys.J.Int., 109(1): 65-77.

    Rougier E, Patton H J, Knight E E, et al. 2011. Constraints on burial depth and yield of the 25 May 2009 North Korean test from hydrodynamic simulations in a granite medium.Geophys.Res.Lett., 38(16),doi: 10.1029/2011GL048269.

    Salzberg D H. 2008. Semi-empirical yield estimates for the 2006 North Korean Explosion. Science Applications International Corp San Diego CA.

    Salzberg D H, Marshall M E. 2007. Seismic source locations and parameters for sparse networks by matching observed seismograms to semi-empirical synthetic seismograms. Science Applications International Corp Mclean VA. Schlittenhardt J, Canty M, Grünberg I. 2010. Satellite Earth observations support CTBT monitoring: A case study of the nuclear test in North Korea of Oct. 9, 2006 and comparison with seismic results.PureandAppliedGeophysics, 167(4-5): 601-618.

    Shin J S, Sheen D H, Kim G. 2010. Regional observations of the second North Korean nuclear test on 2009 May 25.Geophys.J.Int., 180(1): 243-250.

    Tibuleac I M, von Seggern D H, Anderson J G, et al. 2008. Location and magnitude estimation of the 9 October 2006 Korean nuclear explosion using the southern Great Basin digital seismic network as a large-aperture array.Bull.Seismol.Soc.Am., 98(2): 756-767.

    Vergino E S, Mensing R W. 1990. Yield estimation using regional mb(Pn).Bull.Seismol.Soc.Am., 80(3): 656-674.

    Walter W R, Taylor S R. 2001. A revised magnitude and distance amplitude correction (MDAC2) procedure for regional seismic discriminants: Theory and testing at NTS. Lawrence Livermore National Laboratory (LLNL), Livermore, CA.

    Wen L X, Long H. 2010. High-precision location of North Korea′s 2009 nuclear test.SeismologicalResearchLetters, 81(1): 26-29. Xie J, Wu Z, Liu R, et al. 2006. Tomographic regionalization of crustal LgQin eastern Eurasia.Geophys.Res.Lett., 33(3),doi: 10.1029/2005GL024410.

    Yoo S H, Rhie J, Choi H, et al. 2011. Coda-derived source parameters of earthquakes and their scaling relationships in the Korean Peninsula.Bull.Seismol.Soc.Am., 101(5): 2388-2398.Zhang G M, Wang S Y, Wang H. 2012. Tectonic implication of crust and uppermost mantle Pnvelocity structure of eastern Asia.SeismologyandGeology(in Chinese), 34(4): 586-596.

    Zhang M, Wen L X. 2013. High-precision location and yield of North Korea′s 2013 nuclear test.Geophys.Res.Lett., 40(12): 2941-2946. Zhao L F, Wang W M, Li J, et al. 2005. Yield estimation of the 11 May 1998 Indian nuclear test using relative static source strength.ChineseJ.Geophys. (in Chinese), 48(5): 1092-1097. Zhao L F. 2008. Regional seismic characteristics of the North Korean nuclear explosion and LgQtomography [Postdoctoral research report] (in Chinese). Beijing: Institute of Geology and Geophysics, Chinese Academy of Sciences.

    Zhao L F, Xie X B, Wang W M, et al. 2008. Regional seismic characteristics of the 9 October 2006 North Korean nuclear test.Bull.Seismol.Soc.Am., 98(6): 2571-2589.

    Zhao L F, Xie X B, Wang W M, et al. 2012. Yield estimation of the 25 May 2009 North Korean nuclear explosion.Bull.Seismol.Soc.Am., 102(2): 467-478.

    Zhao L F, Xie X B, Wang W M, et al. 2014. The 12 February 2013 North Korean underground nuclear test.SeismologicalResearchLetters, 85(1): 130-134.

    Zhao L F, Xie X B, Tian B F, et al. 2015. Pn wave geometrical spreading and attenuation in Northeast China and the Korean Peninsula constrained by observations from North Korean nuclear explosions.J.Geophys.Res.:SolidEarth, 120(11): 7558-7571.

    Zhu T F, Chun K Y, West G F. 1991. Geometrical spreading andQof Pn waves: An investigative study in eastern Canada.Bull.Seismol.Soc.Am., 81(3): 882-896.

    附中文參考文獻(xiàn)

    賈軍, 周仕勇, 孟憲森等. 2011. 東北地區(qū)Pn波速度結(jié)構(gòu)研究. 防災(zāi)減災(zāi)學(xué)報, 27(2): 53-59.

    裴順平. 2002. 中國大陸上地幔頂部體波速度層析成像[博士論文]. 北京: 中國地震局地球物理研究所.

    裴順平, 許忠淮, 汪素云. 2004a. 中國及鄰區(qū)Pn波速度結(jié)構(gòu)成因探討. 地震學(xué)報, 6(1): 1-10.

    裴順平, 許忠淮, 汪素云. 2004b. 中國大陸及鄰近地區(qū)上地幔頂部Sn波速度層析成像. 地球物理學(xué)報, 47(2): 250-256.

    張國民, 汪素云, 王輝. 2012. 東亞地區(qū)上地幔頂部Pn速度結(jié)構(gòu)的構(gòu)造動力學(xué)意義. 地震地質(zhì), 34(4): 586-596.

    趙連鋒, 王衛(wèi)民, 李娟等. 2005. 利用相對靜力學(xué)強(qiáng)度測定1998年5月11日印度核爆的當(dāng)量. 地球物理學(xué)報, 48(5): 1092-1097.

    趙連鋒. 2008. 北朝鮮核爆的區(qū)域地震特征與Lg波Q值成像[博士后出站報告]. 北京: 中國科學(xué)院地質(zhì)與地球物理研究所.

    (本文編輯何燕)

    基金項(xiàng)目國家自然科學(xué)基金重點(diǎn)項(xiàng)目(41390455)資助.

    作者簡介林鑫,男,1988年生,博士,主要從事朝鮮地下核試驗(yàn)的監(jiān)測和估計等相關(guān)研究. E-mail:alyssalinxin@163.com

    doi:10.6038/cjg20160613 中圖分類號P315

    收稿日期2016-02-01,2016-04-27收修定稿

    Yield and burial depth of the North Korean underground nuclear tests constrained by amplitude envelopes of regional seismic waveforms

    LIN Xin1, 2, YAO Zhen-Xing2

    1InstituteofGeophysics,HuazhongUniversityofScienceandTechnology,SchoolofPhysics,HuazhongUniversityofScienceandTechnology,Wuhan430074,China2KeyLaboratoryofEarthandPlanetaryPhysics,InstituteofGeologyandGeophysics,ChineseAcademyofSciences,Beijing100029,China

    AbstractBoth the yield and burial depth of a nuclear explosion are important parameters of an underground nuclear explosion test. According to the regional seismic stations in northeast China, we investigate these two parameters of the North Korean nuclear tests on 9 October 2006, 25 May 2009, 13 February 2013 and 6 January 2016 by calculating the horizontal component waveforms of coda amplitude envelopes using Pn,Pg,Sn and Lg waves. As an example, we use the seismograms recorded at the Mudanjiang (MDJ) station to show the procedure of calculating amplitude envelopes by applying a bandpass filter to horizontal component waveforms. Thus, we obtain stable estimations of amplitude envelopes at regional seismic stations, from which the regional phases are apparently displayed. The time-domain amplitude envelopes of regional phases can be modeled by using source spectral function, geometric spreading, site response, transfer function and coda shape function. Based on the grid searching method, we simultaneously estimate the yields and burial depths of the North Korean nuclear tests. The results show that the yields increase from 0.6±0.2 kt to 3.0±1.5 kt and then 10.0±2.0 kt, decrease to 8.0±2.0 kt according to their origin time. The first test has a shallow burial depth of 150±100 m, for the second one, the depth is estimated to be 350±100 m, while both the 2013 and 2016 North Korean nuclear tests have very similar depths of 500±200 m. Our results are consistent with the previous findings, which suggest that the yield and burial depth of the underground nuclear explosion tests can be simultaneously constrained by using the time-domain horizontal component of amplitude envelopes at a single station or a few stations.

    KeywordsRegional phase; Amplitude envelope; North Korean nuclear test; Yield of explosion; Burial depth

    林鑫, 姚振興. 2016. 利用區(qū)域地震波形振幅包絡(luò)約束朝鮮地下核試驗(yàn)的埋深和當(dāng)量.地球物理學(xué)報,59(6):2066-2079,doi:10.6038/cjg20160613.

    Lin X, Yao Z X. 2016. Yield and burial depth of the North Korean underground nuclear tests constrained by amplitude envelopes of regional seismic waveforms.ChineseJ.Geophys. (in Chinese),59(6):2066-2079,doi:10.6038/cjg20160613.

    男插女下体视频免费在线播放| 亚洲国产高清在线一区二区三| 成人亚洲精品av一区二区| 国产伦人伦偷精品视频| 精华霜和精华液先用哪个| 精品不卡国产一区二区三区| 久久精品亚洲精品国产色婷小说| 亚洲欧美日韩高清专用| 18禁国产床啪视频网站| 1000部很黄的大片| 亚洲国产日韩欧美精品在线观看 | 波多野结衣高清无吗| 一级毛片高清免费大全| 夜夜看夜夜爽夜夜摸| 啦啦啦观看免费观看视频高清| 麻豆av在线久日| 伦理电影免费视频| 亚洲国产日韩欧美精品在线观看 | 亚洲av电影在线进入| 伊人久久大香线蕉亚洲五| 亚洲专区字幕在线| 国产激情久久老熟女| 欧美3d第一页| 给我免费播放毛片高清在线观看| 亚洲精品美女久久av网站| 午夜福利在线在线| 欧美激情久久久久久爽电影| 免费看光身美女| 黄频高清免费视频| 久久人人精品亚洲av| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲国产中文字幕在线视频| 青草久久国产| 亚洲av电影在线进入| 宅男免费午夜| 精华霜和精华液先用哪个| 亚洲电影在线观看av| 好看av亚洲va欧美ⅴa在| 亚洲在线自拍视频| 成人18禁在线播放| 免费在线观看亚洲国产| 一边摸一边抽搐一进一小说| 亚洲国产欧美人成| 在线免费观看的www视频| 午夜a级毛片| 国产乱人视频| 激情在线观看视频在线高清| 九色国产91popny在线| 精品一区二区三区四区五区乱码| 搡老妇女老女人老熟妇| 亚洲无线在线观看| 少妇丰满av| 2021天堂中文幕一二区在线观| 国产亚洲精品综合一区在线观看| 亚洲 欧美一区二区三区| 久久久久性生活片| 岛国在线免费视频观看| 成年女人永久免费观看视频| 亚洲成人中文字幕在线播放| 一级黄色大片毛片| 一级毛片高清免费大全| 国产成人av激情在线播放| 麻豆成人午夜福利视频| 欧美日韩乱码在线| 亚洲av免费在线观看| 2021天堂中文幕一二区在线观| 欧美日韩综合久久久久久 | 男女之事视频高清在线观看| 亚洲片人在线观看| 国产男靠女视频免费网站| 成年女人看的毛片在线观看| 亚洲国产高清在线一区二区三| 日韩欧美在线乱码| 亚洲国产精品久久男人天堂| 欧美日韩一级在线毛片| 色播亚洲综合网| 国产精品日韩av在线免费观看| 国产蜜桃级精品一区二区三区| 97碰自拍视频| 动漫黄色视频在线观看| 色综合婷婷激情| 制服人妻中文乱码| 亚洲人成伊人成综合网2020| 午夜福利在线观看吧| 久久久精品欧美日韩精品| 99久久精品热视频| 99精品在免费线老司机午夜| 日本黄大片高清| 国产黄片美女视频| 久久久久免费精品人妻一区二区| 一本综合久久免费| 韩国av一区二区三区四区| 午夜精品一区二区三区免费看| 亚洲男人的天堂狠狠| 欧美中文综合在线视频| 黄色丝袜av网址大全| 国产高清videossex| e午夜精品久久久久久久| 很黄的视频免费| 国产精品美女特级片免费视频播放器 | 国产精品一区二区三区四区免费观看 | 中文字幕高清在线视频| 99热这里只有精品一区 | 成人18禁在线播放| 性欧美人与动物交配| 日本一二三区视频观看| 亚洲av熟女| 精品久久久久久久毛片微露脸| 国产私拍福利视频在线观看| 美女大奶头视频| 成人国产一区最新在线观看| av福利片在线观看| 亚洲va日本ⅴa欧美va伊人久久| 看黄色毛片网站| 精品久久久久久久久久免费视频| 少妇的丰满在线观看| 午夜激情欧美在线| 最近在线观看免费完整版| 色播亚洲综合网| 亚洲狠狠婷婷综合久久图片| 国产精品一区二区三区四区久久| 1024香蕉在线观看| 床上黄色一级片| 一区福利在线观看| 国产一级毛片七仙女欲春2| 成人特级黄色片久久久久久久| 欧美激情久久久久久爽电影| 欧美高清成人免费视频www| 国产探花在线观看一区二区| 级片在线观看| 法律面前人人平等表现在哪些方面| 性欧美人与动物交配| 亚洲人成伊人成综合网2020| tocl精华| 午夜福利视频1000在线观看| 亚洲七黄色美女视频| 亚洲国产欧美一区二区综合| av视频在线观看入口| 日韩有码中文字幕| 久久精品91无色码中文字幕| 嫩草影院精品99| 国产一区二区激情短视频| 亚洲,欧美精品.| 他把我摸到了高潮在线观看| 久久精品国产亚洲av香蕉五月| 久久久成人免费电影| 草草在线视频免费看| 18禁国产床啪视频网站| 香蕉丝袜av| 日韩欧美三级三区| 1024香蕉在线观看| 国产伦精品一区二区三区视频9 | 国产人伦9x9x在线观看| 亚洲第一电影网av| 丰满人妻熟妇乱又伦精品不卡| 伦理电影免费视频| 久久午夜综合久久蜜桃| 欧美国产日韩亚洲一区| 99热这里只有精品一区 | 亚洲av中文字字幕乱码综合| 国产精品一区二区免费欧美| 国产精品永久免费网站| 国产成人欧美在线观看| 桃色一区二区三区在线观看| 亚洲成a人片在线一区二区| 女警被强在线播放| 精品一区二区三区视频在线 | 欧美绝顶高潮抽搐喷水| 一区二区三区国产精品乱码| 国产激情欧美一区二区| 免费观看精品视频网站| 啦啦啦韩国在线观看视频| 99久久99久久久精品蜜桃| 亚洲av片天天在线观看| 欧美日韩瑟瑟在线播放| 久久香蕉精品热| 国内少妇人妻偷人精品xxx网站 | 一个人观看的视频www高清免费观看 | 美女黄网站色视频| 国产av不卡久久| 亚洲性夜色夜夜综合| 少妇裸体淫交视频免费看高清| 狠狠狠狠99中文字幕| 黄频高清免费视频| a级毛片在线看网站| 香蕉国产在线看| 午夜激情福利司机影院| 亚洲国产中文字幕在线视频| 级片在线观看| 午夜成年电影在线免费观看| 精品国内亚洲2022精品成人| 精品久久久久久久人妻蜜臀av| 国产伦一二天堂av在线观看| 女同久久另类99精品国产91| 97碰自拍视频| 九九久久精品国产亚洲av麻豆 | 欧美日韩中文字幕国产精品一区二区三区| 国产av不卡久久| 欧美性猛交黑人性爽| 香蕉丝袜av| cao死你这个sao货| 男女视频在线观看网站免费| 国产成人影院久久av| 久久九九热精品免费| 午夜精品久久久久久毛片777| 亚洲人成电影免费在线| xxxwww97欧美| 免费看日本二区| 亚洲第一电影网av| 91老司机精品| 成人精品一区二区免费| 亚洲专区字幕在线| 看黄色毛片网站| 欧美日韩国产亚洲二区| 黄色成人免费大全| 久久欧美精品欧美久久欧美| 日韩 欧美 亚洲 中文字幕| 男女午夜视频在线观看| 一a级毛片在线观看| 老汉色av国产亚洲站长工具| av在线天堂中文字幕| 制服丝袜大香蕉在线| 最近最新免费中文字幕在线| 男人的好看免费观看在线视频| 欧美日韩福利视频一区二区| 桃色一区二区三区在线观看| 亚洲avbb在线观看| 黑人欧美特级aaaaaa片| 精品国产乱子伦一区二区三区| 国产精品野战在线观看| 老熟妇乱子伦视频在线观看| 99国产极品粉嫩在线观看| 精品久久久久久,| 国产又色又爽无遮挡免费看| 国产91精品成人一区二区三区| 日本精品一区二区三区蜜桃| 国产一区二区激情短视频| 搡老妇女老女人老熟妇| 一本久久中文字幕| 免费在线观看成人毛片| 亚洲av电影不卡..在线观看| 亚洲av成人精品一区久久| 少妇熟女aⅴ在线视频| 老汉色∧v一级毛片| 五月伊人婷婷丁香| 亚洲国产欧美人成| 制服丝袜大香蕉在线| 成人午夜高清在线视频| 久久午夜亚洲精品久久| 免费在线观看亚洲国产| 欧美中文综合在线视频| 看免费av毛片| h日本视频在线播放| 国产午夜福利久久久久久| 99re在线观看精品视频| 淫秽高清视频在线观看| 男人和女人高潮做爰伦理| 欧美大码av| 国产亚洲精品一区二区www| h日本视频在线播放| 嫩草影视91久久| 国产精品av视频在线免费观看| 成年女人永久免费观看视频| 啪啪无遮挡十八禁网站| 国产成人aa在线观看| 国产成人啪精品午夜网站| 国产伦精品一区二区三区四那| 99国产精品一区二区三区| 国产极品精品免费视频能看的| 亚洲中文字幕一区二区三区有码在线看 | 亚洲人成网站高清观看| 国产亚洲精品av在线| 亚洲国产高清在线一区二区三| 国产精品久久久av美女十八| 国产精品一区二区三区四区久久| 中亚洲国语对白在线视频| 色综合欧美亚洲国产小说| 看免费av毛片| av在线蜜桃| 亚洲片人在线观看| 久久久久久国产a免费观看| 欧美黑人巨大hd| 国产成人av教育| 亚洲天堂国产精品一区在线| 亚洲国产日韩欧美精品在线观看 | 熟妇人妻久久中文字幕3abv| 国产精品久久久久久久电影 | 国产高清videossex| 国产欧美日韩精品一区二区| 国产成人精品久久二区二区免费| 欧美不卡视频在线免费观看| 999精品在线视频| 亚洲av成人不卡在线观看播放网| 国产成人系列免费观看| 欧美黄色片欧美黄色片| 久久精品综合一区二区三区| 亚洲av熟女| 在线观看日韩欧美| 欧美一级毛片孕妇| 欧美大码av| 久久久久免费精品人妻一区二区| 欧美日韩精品网址| 日韩欧美 国产精品| 99re在线观看精品视频| 黄色女人牲交| 亚洲av日韩精品久久久久久密| 18禁美女被吸乳视频| 男人和女人高潮做爰伦理| 欧洲精品卡2卡3卡4卡5卡区| 国产欧美日韩精品一区二区| 欧美不卡视频在线免费观看| 日韩免费av在线播放| 日韩欧美三级三区| 国产亚洲欧美在线一区二区| 亚洲黑人精品在线| 精品一区二区三区视频在线 | 丰满的人妻完整版| 免费电影在线观看免费观看| 国产一区二区三区在线臀色熟女| 国产淫片久久久久久久久 | 91av网站免费观看| 好男人在线观看高清免费视频| 日本三级黄在线观看| 亚洲欧美日韩东京热| 99久久久亚洲精品蜜臀av| 国产91精品成人一区二区三区| 男女视频在线观看网站免费| 99久久综合精品五月天人人| 久久久久免费精品人妻一区二区| 两个人的视频大全免费| 美女黄网站色视频| 久久九九热精品免费| 亚洲第一电影网av| 国产免费男女视频| 久久久久亚洲av毛片大全| 国产欧美日韩精品亚洲av| 亚洲精品久久国产高清桃花| 成人三级做爰电影| 国产亚洲精品av在线| 悠悠久久av| 九九热线精品视视频播放| 免费搜索国产男女视频| 九色国产91popny在线| 在线十欧美十亚洲十日本专区| 国产一区在线观看成人免费| 听说在线观看完整版免费高清| x7x7x7水蜜桃| 亚洲一区二区三区不卡视频| 午夜福利18| 免费人成视频x8x8入口观看| 午夜久久久久精精品| 午夜两性在线视频| 久久精品夜夜夜夜夜久久蜜豆| 少妇裸体淫交视频免费看高清| 美女高潮的动态| 亚洲成av人片在线播放无| 国产精品综合久久久久久久免费| 亚洲精品美女久久久久99蜜臀| 香蕉国产在线看| 亚洲欧美日韩高清在线视频| 麻豆久久精品国产亚洲av| 久久精品影院6| 熟女人妻精品中文字幕| 看片在线看免费视频| 好男人电影高清在线观看| 色综合站精品国产| 无人区码免费观看不卡| 午夜日韩欧美国产| 色在线成人网| 国产成人精品久久二区二区91| 国产精品 欧美亚洲| 日本黄色视频三级网站网址| 国产伦精品一区二区三区视频9 | 欧美日韩黄片免| 日韩欧美三级三区| 久久午夜亚洲精品久久| 亚洲 欧美一区二区三区| 午夜免费成人在线视频| 日本免费a在线| 成年免费大片在线观看| 欧美高清成人免费视频www| 神马国产精品三级电影在线观看| 国产精品一区二区三区四区久久| 法律面前人人平等表现在哪些方面| 欧美日韩国产亚洲二区| 亚洲成人精品中文字幕电影| 制服丝袜大香蕉在线| 少妇的丰满在线观看| 老司机深夜福利视频在线观看| 97人妻精品一区二区三区麻豆| 午夜亚洲福利在线播放| 国产v大片淫在线免费观看| www日本在线高清视频| 人人妻,人人澡人人爽秒播| 午夜福利免费观看在线| 老司机午夜十八禁免费视频| 好看av亚洲va欧美ⅴa在| 国产成人系列免费观看| 69av精品久久久久久| 韩国av一区二区三区四区| 亚洲国产精品999在线| 91av网站免费观看| 色av中文字幕| 网址你懂的国产日韩在线| 久久精品国产亚洲av香蕉五月| 欧美成狂野欧美在线观看| 黑人操中国人逼视频| 搡老岳熟女国产| 嫩草影院入口| 91老司机精品| 男人舔女人的私密视频| 国产成人精品久久二区二区91| 午夜久久久久精精品| 欧美日韩黄片免| 国产精品九九99| 桃色一区二区三区在线观看| 免费在线观看视频国产中文字幕亚洲| 国产成人系列免费观看| 99久久无色码亚洲精品果冻| 制服人妻中文乱码| 久久久久久人人人人人| or卡值多少钱| 欧美一级毛片孕妇| 国产精品1区2区在线观看.| 91麻豆精品激情在线观看国产| 给我免费播放毛片高清在线观看| 日日摸夜夜添夜夜添小说| 免费看a级黄色片| 亚洲黑人精品在线| 国产综合懂色| 午夜福利在线在线| 色播亚洲综合网| 国产精品电影一区二区三区| 色在线成人网| 青草久久国产| 99久久99久久久精品蜜桃| 在线视频色国产色| 男女视频在线观看网站免费| 亚洲在线自拍视频| 国产精品影院久久| 99国产精品一区二区蜜桃av| 男女那种视频在线观看| 国产精品影院久久| 久久久久国产精品人妻aⅴ院| 在线永久观看黄色视频| 在线观看日韩欧美| 欧美不卡视频在线免费观看| 床上黄色一级片| 老司机午夜十八禁免费视频| 99国产精品一区二区三区| 美女免费视频网站| 村上凉子中文字幕在线| 日韩欧美三级三区| 国产欧美日韩一区二区精品| 日韩国内少妇激情av| 国产91精品成人一区二区三区| 亚洲成人久久爱视频| 久久久久国产精品人妻aⅴ院| 国产精品九九99| 午夜成年电影在线免费观看| 日韩人妻高清精品专区| 夜夜爽天天搞| 国产亚洲精品av在线| 日韩欧美三级三区| 欧美中文日本在线观看视频| 麻豆久久精品国产亚洲av| 色播亚洲综合网| 日日摸夜夜添夜夜添小说| 午夜a级毛片| 亚洲人与动物交配视频| 手机成人av网站| 免费高清视频大片| 99久久综合精品五月天人人| 12—13女人毛片做爰片一| 亚洲国产精品成人综合色| 小蜜桃在线观看免费完整版高清| 99久久精品热视频| e午夜精品久久久久久久| 亚洲熟女毛片儿| 在线观看美女被高潮喷水网站 | 99精品在免费线老司机午夜| 午夜a级毛片| av片东京热男人的天堂| 欧美日韩黄片免| 99久久无色码亚洲精品果冻| 丁香六月欧美| 国产伦精品一区二区三区视频9 | 变态另类丝袜制服| 婷婷丁香在线五月| 亚洲 国产 在线| 亚洲五月婷婷丁香| 欧美激情久久久久久爽电影| 露出奶头的视频| 亚洲在线观看片| 久久婷婷人人爽人人干人人爱| 亚洲成人久久爱视频| 亚洲欧美精品综合一区二区三区| 九九热线精品视视频播放| 级片在线观看| 亚洲国产精品sss在线观看| 亚洲色图 男人天堂 中文字幕| a级毛片a级免费在线| 亚洲五月婷婷丁香| 亚洲欧美日韩东京热| 国产亚洲欧美在线一区二区| 丰满人妻熟妇乱又伦精品不卡| 午夜福利在线在线| 美女cb高潮喷水在线观看 | 国产午夜精品久久久久久| 日本在线视频免费播放| 色在线成人网| 日韩有码中文字幕| 欧美三级亚洲精品| 99国产精品一区二区蜜桃av| 黑人操中国人逼视频| 久久久久久国产a免费观看| 久久人妻av系列| 国产高清激情床上av| 久久草成人影院| 男人舔奶头视频| 欧美乱妇无乱码| 欧美日韩综合久久久久久 | 不卡一级毛片| 欧美zozozo另类| 中文资源天堂在线| 99在线人妻在线中文字幕| 丰满人妻熟妇乱又伦精品不卡| 亚洲成人精品中文字幕电影| 亚洲人成电影免费在线| 小说图片视频综合网站| 亚洲七黄色美女视频| 99久久99久久久精品蜜桃| 免费高清视频大片| 精品久久久久久久久久久久久| 91av网一区二区| 日本黄色片子视频| 欧美+亚洲+日韩+国产| 亚洲五月婷婷丁香| 国产视频一区二区在线看| 国产精品久久久久久人妻精品电影| 亚洲一区二区三区不卡视频| 九九热线精品视视频播放| 国产av麻豆久久久久久久| 日韩三级视频一区二区三区| 国产人伦9x9x在线观看| 一a级毛片在线观看| 99久国产av精品| 午夜福利18| 免费在线观看亚洲国产| 日韩免费av在线播放| 欧美色视频一区免费| 国产一区二区三区视频了| 午夜福利视频1000在线观看| 小说图片视频综合网站| 中亚洲国语对白在线视频| 成年版毛片免费区| 亚洲av熟女| 国内精品一区二区在线观看| 一夜夜www| 一本综合久久免费| 最近视频中文字幕2019在线8| 一区二区三区国产精品乱码| 国产精品日韩av在线免费观看| 美女被艹到高潮喷水动态| 国产精品一区二区三区四区久久| 国产精品 欧美亚洲| www日本黄色视频网| 国产亚洲欧美98| 婷婷精品国产亚洲av在线| 熟女人妻精品中文字幕| av在线蜜桃| 欧美色视频一区免费| 亚洲精品在线美女| 欧美成人一区二区免费高清观看 | 国产午夜精品久久久久久| 高潮久久久久久久久久久不卡| 97超视频在线观看视频| 国产又黄又爽又无遮挡在线| 丰满的人妻完整版| 99久久成人亚洲精品观看| 亚洲片人在线观看| 亚洲avbb在线观看| 国产一区二区在线观看日韩 | 国产精品自产拍在线观看55亚洲| 99热精品在线国产| 国产成人福利小说| 脱女人内裤的视频| 日韩三级视频一区二区三区| 啦啦啦韩国在线观看视频| 免费高清视频大片| 免费观看精品视频网站| 国产欧美日韩精品亚洲av| 免费在线观看日本一区| 九色成人免费人妻av| 亚洲 欧美一区二区三区| 床上黄色一级片| 看片在线看免费视频| 久久久国产精品麻豆| 国产成人精品久久二区二区免费| 噜噜噜噜噜久久久久久91| 亚洲av片天天在线观看| 嫩草影院精品99| 国产熟女xx| 性色av乱码一区二区三区2| 免费在线观看视频国产中文字幕亚洲| 18禁美女被吸乳视频| 色老头精品视频在线观看| 亚洲五月天丁香| 成人特级黄色片久久久久久久| 黄片大片在线免费观看| 夜夜看夜夜爽夜夜摸| 最近在线观看免费完整版| 97碰自拍视频| 国产精品亚洲av一区麻豆| 99re在线观看精品视频|