• <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男天堂| 夫妻午夜视频| 最近中文字幕2019免费版| 午夜91福利影院| 最近最新免费中文字幕在线| avwww免费| 久久久精品94久久精品| 51午夜福利影视在线观看| 丝袜在线中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 99国产精品免费福利视频| 黑人操中国人逼视频| 成人国语在线视频| 国产精品 国内视频| 精品福利观看| 法律面前人人平等表现在哪些方面 | 可以免费在线观看a视频的电影网站| 国产福利在线免费观看视频| 热99久久久久精品小说推荐| 久久久久久久大尺度免费视频| 欧美日韩国产mv在线观看视频| 亚洲综合色网址| 老司机影院毛片| 久久人人爽人人片av| 国产男人的电影天堂91| 欧美在线黄色| 国产日韩一区二区三区精品不卡| 伊人亚洲综合成人网| 91成人精品电影| 国产亚洲av高清不卡| 国产1区2区3区精品| 国产精品亚洲av一区麻豆| 中亚洲国语对白在线视频| 久久精品国产a三级三级三级| 欧美精品av麻豆av| 国产一区二区三区在线臀色熟女 | 免费日韩欧美在线观看| www.精华液| 久久av网站| 一本—道久久a久久精品蜜桃钙片| 国产一区二区三区综合在线观看| 午夜免费成人在线视频| 国产激情久久老熟女| 操出白浆在线播放| 午夜视频精品福利| 在线av久久热| 女人久久www免费人成看片| 久久久久久人人人人人| 日本wwww免费看| 男女边摸边吃奶| 男人添女人高潮全过程视频| 亚洲三区欧美一区| 丰满少妇做爰视频| 欧美另类一区| 菩萨蛮人人尽说江南好唐韦庄| 51午夜福利影视在线观看| 在线天堂中文资源库| 丝袜喷水一区| 欧美性长视频在线观看| 久久99热这里只频精品6学生| 桃花免费在线播放| 国产有黄有色有爽视频| 日本欧美视频一区| 午夜福利影视在线免费观看| 久久这里只有精品19| 日韩,欧美,国产一区二区三区| 亚洲精品国产区一区二| 久久ye,这里只有精品| 亚洲黑人精品在线| 国产日韩欧美在线精品| 色老头精品视频在线观看| 亚洲国产欧美在线一区| 人成视频在线观看免费观看| 亚洲国产看品久久| 91精品国产国语对白视频| 两人在一起打扑克的视频| 亚洲国产av新网站| 欧美日韩一级在线毛片| 99精品欧美一区二区三区四区| 国产一区有黄有色的免费视频| 少妇粗大呻吟视频| 五月天丁香电影| 91精品伊人久久大香线蕉| 精品第一国产精品| 国产欧美日韩一区二区精品| 一本一本久久a久久精品综合妖精| 建设人人有责人人尽责人人享有的| 国产成人精品久久二区二区免费| 国产国语露脸激情在线看| 免费观看人在逋| 亚洲精品一区蜜桃| 99久久99久久久精品蜜桃| 欧美日韩亚洲国产一区二区在线观看 | 精品久久蜜臀av无| 国产精品一区二区在线观看99| 老熟妇乱子伦视频在线观看 | 别揉我奶头~嗯~啊~动态视频 | 久久久久国内视频| 欧美日本中文国产一区发布| 亚洲少妇的诱惑av| 亚洲欧美激情在线| 黄色毛片三级朝国网站| 亚洲精品国产一区二区精华液| 精品一区在线观看国产| 免费女性裸体啪啪无遮挡网站| 国产淫语在线视频| 亚洲第一av免费看| 一区福利在线观看| 亚洲成人免费电影在线观看| 两个人看的免费小视频| 久久久精品免费免费高清| e午夜精品久久久久久久| 亚洲性夜色夜夜综合| 18禁国产床啪视频网站| 99久久人妻综合| 五月天丁香电影| 99久久99久久久精品蜜桃| xxxhd国产人妻xxx| 欧美另类一区| 欧美中文综合在线视频| 99热网站在线观看| 中文欧美无线码| 两性夫妻黄色片| 69av精品久久久久久 | 久久久久久久大尺度免费视频| 嫩草影视91久久| av天堂在线播放| 国产xxxxx性猛交| 男人添女人高潮全过程视频| av视频免费观看在线观看| 久久精品aⅴ一区二区三区四区| 国产一级毛片在线| 黄色a级毛片大全视频| tocl精华| 成人亚洲精品一区在线观看| 亚洲 国产 在线| 大片免费播放器 马上看| 一进一出抽搐动态| 国产伦人伦偷精品视频| 久热这里只有精品99| 热re99久久精品国产66热6| xxxhd国产人妻xxx| 在线 av 中文字幕| 超色免费av| 中文字幕人妻丝袜一区二区| 婷婷成人精品国产| 色播在线永久视频| 国产区一区二久久| 国产在线视频一区二区| 亚洲国产精品一区三区| 亚洲精品一二三| 午夜福利乱码中文字幕| 一级黄色大片毛片| 91精品国产国语对白视频| 热re99久久国产66热| 18禁裸乳无遮挡动漫免费视频| 久久久久视频综合| 成人三级做爰电影| 亚洲精品国产一区二区精华液| 久久香蕉激情| 亚洲国产欧美一区二区综合| 成人手机av| 免费在线观看黄色视频的| 女人精品久久久久毛片| 午夜久久久在线观看| 亚洲av男天堂| 国产一区二区在线观看av| tocl精华| 国产黄频视频在线观看| 亚洲国产欧美一区二区综合| 男女之事视频高清在线观看| 香蕉丝袜av| 欧美日韩亚洲综合一区二区三区_| 亚洲人成77777在线视频| 日韩熟女老妇一区二区性免费视频| 啦啦啦 在线观看视频| 国产在线免费精品| 日韩,欧美,国产一区二区三区| 亚洲av电影在线进入| 国产男女超爽视频在线观看| 高清视频免费观看一区二区| 国产三级黄色录像| 亚洲第一青青草原| 久久av网站| 蜜桃在线观看..| 一本一本久久a久久精品综合妖精| www.av在线官网国产| h视频一区二区三区| 十八禁高潮呻吟视频| 欧美在线黄色| 久久久水蜜桃国产精品网| 51午夜福利影视在线观看| 超色免费av| 成人三级做爰电影| 高清av免费在线| 国产高清国产精品国产三级| 人人澡人人妻人| 亚洲精品av麻豆狂野| 99久久99久久久精品蜜桃| 亚洲专区字幕在线| 另类精品久久| 丝瓜视频免费看黄片| 涩涩av久久男人的天堂| 国产伦理片在线播放av一区| 亚洲欧美精品自产自拍| 热99re8久久精品国产| 国产亚洲午夜精品一区二区久久| 国产成人av激情在线播放| 久久 成人 亚洲| 色老头精品视频在线观看| 国产高清视频在线播放一区 | 免费观看av网站的网址| 亚洲精品久久久久久婷婷小说| 亚洲国产av新网站| 制服人妻中文乱码| 欧美日韩国产mv在线观看视频| 国产精品一区二区在线不卡| 亚洲中文日韩欧美视频| 欧美另类一区| 男女高潮啪啪啪动态图| 女人精品久久久久毛片| 99久久精品国产亚洲精品| 三级毛片av免费| 伦理电影免费视频| 久久国产精品人妻蜜桃| 9191精品国产免费久久| 亚洲欧美激情在线| 久久久国产欧美日韩av| 18禁国产床啪视频网站| 欧美激情高清一区二区三区| 久久性视频一级片| 桃花免费在线播放| 色视频在线一区二区三区| 夜夜骑夜夜射夜夜干| 精品国产超薄肉色丝袜足j| 日本五十路高清| 久久久久精品人妻al黑| 日本av免费视频播放| 桃花免费在线播放| 欧美另类亚洲清纯唯美| 亚洲精品国产色婷婷电影| 亚洲国产精品一区三区| 精品国产乱子伦一区二区三区 | 十八禁高潮呻吟视频| 精品国产一区二区三区四区第35| 青春草视频在线免费观看| 国产91精品成人一区二区三区 | 考比视频在线观看| 人人妻人人爽人人添夜夜欢视频| 男女免费视频国产| 成人黄色视频免费在线看| 免费观看a级毛片全部| 中文精品一卡2卡3卡4更新| 亚洲欧洲日产国产| 久久 成人 亚洲| 老司机影院成人| 亚洲熟女毛片儿| 日韩电影二区| 丝袜美腿诱惑在线| 欧美亚洲 丝袜 人妻 在线| 午夜久久久在线观看| 欧美日韩国产mv在线观看视频| 亚洲精品中文字幕在线视频| 宅男免费午夜| 一级片免费观看大全| 成年女人毛片免费观看观看9 | 成人影院久久| 精品国产乱子伦一区二区三区 | 国产男女超爽视频在线观看| 国产欧美日韩一区二区精品| 高清av免费在线| 久久99热这里只频精品6学生| 国产精品熟女久久久久浪| 不卡一级毛片| 国产免费av片在线观看野外av| 天天添夜夜摸| 国产极品粉嫩免费观看在线| 久久综合国产亚洲精品| 欧美激情久久久久久爽电影 | 国产成人av教育| 正在播放国产对白刺激| 不卡av一区二区三区| 一个人免费看片子| 99精国产麻豆久久婷婷| 亚洲精品中文字幕一二三四区 | 十八禁网站网址无遮挡| 久久亚洲国产成人精品v| 国产又色又爽无遮挡免| 久热这里只有精品99| 午夜福利一区二区在线看| 免费在线观看影片大全网站| 国产精品麻豆人妻色哟哟久久| 91成人精品电影| 日本av手机在线免费观看| 国产精品一区二区精品视频观看| 欧美av亚洲av综合av国产av| 国产三级黄色录像| 免费人妻精品一区二区三区视频| 欧美日韩亚洲国产一区二区在线观看 | 日韩一卡2卡3卡4卡2021年| av在线播放精品| 午夜免费成人在线视频| 亚洲九九香蕉| 亚洲av男天堂| 精品国产乱子伦一区二区三区 | 黄色视频不卡| 伊人亚洲综合成人网| 老司机午夜十八禁免费视频| 9色porny在线观看| 久久久久网色| 精品人妻1区二区| 国产91精品成人一区二区三区 | 国产精品九九99| 丝袜人妻中文字幕| 国产男人的电影天堂91| 亚洲中文日韩欧美视频| 动漫黄色视频在线观看| 国产精品香港三级国产av潘金莲| 啦啦啦 在线观看视频| 国产精品 国内视频| 麻豆国产av国片精品| 欧美大码av| 丝袜脚勾引网站| 欧美+亚洲+日韩+国产| 国产不卡av网站在线观看| 永久免费av网站大全| av又黄又爽大尺度在线免费看| 一边摸一边抽搐一进一出视频| 精品亚洲成国产av| 久久国产精品大桥未久av| 99久久精品国产亚洲精品| 真人做人爱边吃奶动态| 午夜免费成人在线视频| 日韩欧美一区二区三区在线观看 | 女人精品久久久久毛片| 亚洲精品第二区| 久久精品国产亚洲av香蕉五月 | 午夜福利影视在线免费观看| 韩国高清视频一区二区三区| 久热这里只有精品99| 久久精品国产亚洲av高清一级| 亚洲国产av影院在线观看| 日日爽夜夜爽网站| 好男人电影高清在线观看| 成人免费观看视频高清| 日本一区二区免费在线视频| 高潮久久久久久久久久久不卡| 久热爱精品视频在线9| 国产精品国产av在线观看| 欧美日韩视频精品一区| 热re99久久国产66热| av欧美777| 亚洲色图 男人天堂 中文字幕| 免费在线观看日本一区| 18禁观看日本| 国产免费现黄频在线看| 亚洲av男天堂| 中文字幕精品免费在线观看视频| 精品国产一区二区三区四区第35| 国产精品久久久久久人妻精品电影 | kizo精华| 欧美日韩视频精品一区| 日日爽夜夜爽网站| 丁香六月天网| 欧美激情高清一区二区三区| 老司机午夜福利在线观看视频 | 久久久久久久久久久久大奶| 免费观看a级毛片全部| 两个人看的免费小视频| 男女边摸边吃奶| 91老司机精品| 午夜精品国产一区二区电影| 国产欧美日韩精品亚洲av| 18禁黄网站禁片午夜丰满| 欧美乱码精品一区二区三区| 无遮挡黄片免费观看| 久久精品国产亚洲av高清一级| 可以免费在线观看a视频的电影网站| 99久久99久久久精品蜜桃| 男女午夜视频在线观看| 狠狠狠狠99中文字幕| 老司机福利观看| 一边摸一边做爽爽视频免费| www.av在线官网国产| 国产成人免费观看mmmm| 亚洲欧美清纯卡通| 波多野结衣av一区二区av| 久久久久久久久久久久大奶| 可以免费在线观看a视频的电影网站| 丝袜脚勾引网站| 亚洲男人天堂网一区| 99精品久久久久人妻精品| 中文字幕精品免费在线观看视频| avwww免费| 中文欧美无线码| 午夜激情久久久久久久| 一本久久精品| 日本欧美视频一区| a在线观看视频网站| 欧美国产精品一级二级三级| 丰满少妇做爰视频| 精品国产一区二区三区四区第35| 五月开心婷婷网| 亚洲人成电影免费在线| 伦理电影免费视频| 国产精品一二三区在线看| 精品一品国产午夜福利视频| 成人亚洲精品一区在线观看| 可以免费在线观看a视频的电影网站| xxxhd国产人妻xxx| 精品亚洲成国产av| 精品久久久久久久毛片微露脸 | 国产精品一区二区在线观看99| 欧美精品一区二区免费开放| 午夜免费鲁丝| 亚洲精品国产av蜜桃| 国产精品自产拍在线观看55亚洲 | 伊人久久大香线蕉亚洲五| 成年美女黄网站色视频大全免费| 欧美少妇被猛烈插入视频| 亚洲精品第二区| 国产精品.久久久| 亚洲国产欧美在线一区| 日韩中文字幕欧美一区二区| 色94色欧美一区二区| 亚洲avbb在线观看| 不卡一级毛片| 91成年电影在线观看| 国产男人的电影天堂91| 狠狠精品人妻久久久久久综合| 18在线观看网站| 一本久久精品| 国产欧美日韩精品亚洲av| 高清av免费在线| 韩国精品一区二区三区| 亚洲黑人精品在线| 久久青草综合色| 午夜精品国产一区二区电影| 欧美激情高清一区二区三区| 国产精品免费大片| 超碰97精品在线观看| 日韩三级视频一区二区三区| 亚洲伊人久久精品综合| 在线 av 中文字幕| 久久中文看片网| 99香蕉大伊视频| 99久久99久久久精品蜜桃| 男女无遮挡免费网站观看| 视频在线观看一区二区三区| 精品少妇久久久久久888优播| 韩国精品一区二区三区| 国产成人精品无人区| 亚洲人成电影观看| 可以免费在线观看a视频的电影网站| 99久久综合免费| 久久久国产欧美日韩av| 一区二区三区精品91| 黑人操中国人逼视频| h视频一区二区三区| 国产欧美日韩综合在线一区二区| 国产精品久久久久久精品电影小说| 女性生殖器流出的白浆| 国产成人一区二区三区免费视频网站| 精品久久久精品久久久| 国产不卡av网站在线观看| 日本撒尿小便嘘嘘汇集6| 久热爱精品视频在线9| 日韩熟女老妇一区二区性免费视频| 亚洲精品国产色婷婷电影| 18禁观看日本| 电影成人av| 精品卡一卡二卡四卡免费| 国产又色又爽无遮挡免| 操出白浆在线播放| 欧美国产精品va在线观看不卡| 久久久久精品人妻al黑| 亚洲精品国产一区二区精华液| 中文字幕制服av| 日韩人妻精品一区2区三区| 国产亚洲精品一区二区www | 亚洲精品日韩在线中文字幕| 国产日韩一区二区三区精品不卡| 另类精品久久| 一级,二级,三级黄色视频| 青春草亚洲视频在线观看| 国产一级毛片在线| 久久性视频一级片| 精品国内亚洲2022精品成人 | tube8黄色片| 成年美女黄网站色视频大全免费| 女人精品久久久久毛片| 黑人巨大精品欧美一区二区mp4| 国产真人三级小视频在线观看| 国产黄色免费在线视频| 午夜福利影视在线免费观看| 久久久水蜜桃国产精品网| 人妻久久中文字幕网| 久久人人97超碰香蕉20202| 久久九九热精品免费| a级毛片黄视频| 中文字幕另类日韩欧美亚洲嫩草| 久久久久久久久免费视频了| 免费少妇av软件| 亚洲国产av影院在线观看| 激情视频va一区二区三区| 亚洲国产精品一区二区三区在线| 国产1区2区3区精品| 久久中文看片网| 免费久久久久久久精品成人欧美视频| 999精品在线视频| 免费高清在线观看视频在线观看| 下体分泌物呈黄色| 国内毛片毛片毛片毛片毛片| 久久精品国产综合久久久| 亚洲精品美女久久av网站| 国产精品免费大片| 新久久久久国产一级毛片| 一区二区三区精品91| 人妻人人澡人人爽人人| 午夜福利,免费看| 精品久久蜜臀av无| 老熟妇仑乱视频hdxx| 最近最新免费中文字幕在线| 精品国产国语对白av| 久久亚洲国产成人精品v| 久久久久国产一级毛片高清牌| 在线亚洲精品国产二区图片欧美| 久久精品国产亚洲av高清一级| av国产精品久久久久影院| 80岁老熟妇乱子伦牲交| 欧美黄色淫秽网站| 久久久久久久大尺度免费视频| 中国美女看黄片| 免费日韩欧美在线观看| 色94色欧美一区二区| 久久久久久久久免费视频了| 美女午夜性视频免费| 一区二区av电影网| 国产无遮挡羞羞视频在线观看| 亚洲七黄色美女视频| 亚洲中文日韩欧美视频| 国产91精品成人一区二区三区 | 香蕉国产在线看| 丝袜脚勾引网站| 少妇裸体淫交视频免费看高清 | 亚洲自偷自拍图片 自拍| 国产在线视频一区二区| 亚洲精品国产色婷婷电影| 欧美黑人精品巨大| 免费不卡黄色视频| 日韩制服丝袜自拍偷拍| 亚洲人成电影观看| av免费在线观看网站| 香蕉丝袜av| 美国免费a级毛片| 色94色欧美一区二区| 亚洲 欧美一区二区三区| 另类亚洲欧美激情| 最新在线观看一区二区三区| 精品国产超薄肉色丝袜足j| 欧美精品人与动牲交sv欧美| 在线十欧美十亚洲十日本专区| 丝袜美腿诱惑在线| 99久久综合免费| 国产精品九九99| 高潮久久久久久久久久久不卡| 国产真人三级小视频在线观看| 亚洲免费av在线视频| 亚洲专区中文字幕在线| 欧美亚洲 丝袜 人妻 在线| 777久久人妻少妇嫩草av网站| 男女下面插进去视频免费观看| 女人爽到高潮嗷嗷叫在线视频| 女人久久www免费人成看片| 丁香六月欧美| www.av在线官网国产| 久久久久久久精品精品| 国产精品二区激情视频| 9色porny在线观看| 少妇被粗大的猛进出69影院| 国产精品欧美亚洲77777| 午夜免费观看性视频| 777米奇影视久久| 精品亚洲成a人片在线观看| 日本vs欧美在线观看视频| 精品人妻1区二区| 日韩,欧美,国产一区二区三区| 久久热在线av| 男女床上黄色一级片免费看| 青青草视频在线视频观看| 精品国产一区二区三区四区第35| 嫩草影视91久久| tocl精华| 午夜精品国产一区二区电影| 69精品国产乱码久久久| 欧美精品亚洲一区二区| 成年美女黄网站色视频大全免费| 成年人午夜在线观看视频| 99精国产麻豆久久婷婷| 91精品国产国语对白视频| 国产在线免费精品| 少妇猛男粗大的猛烈进出视频| 国产精品一区二区在线观看99| 丰满少妇做爰视频| 中文字幕人妻丝袜制服| 两性午夜刺激爽爽歪歪视频在线观看 | 视频区欧美日本亚洲| 国产av国产精品国产| 在线天堂中文资源库| 啦啦啦在线免费观看视频4| 精品少妇黑人巨大在线播放| 精品国内亚洲2022精品成人 | 青草久久国产| 婷婷色av中文字幕| 欧美日韩福利视频一区二区|