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

    利用Lg波譜特征估計(jì)地下核爆炸裝置的埋深①

    2015-06-09 12:36:02何永鋒劉炳燦
    地震工程學(xué)報(bào) 2015年2期
    關(guān)鍵詞:核爆炸低谷波譜

    何永鋒, 李 鍇, 劉炳燦

    (裝甲兵工程學(xué)院,北京 100072)

    利用Lg波譜特征估計(jì)地下核爆炸裝置的埋深①

    何永鋒, 李 鍇, 劉炳燦

    (裝甲兵工程學(xué)院,北京 100072)

    Lg波同時(shí)攜帶地球介質(zhì)結(jié)構(gòu)信息和爆炸源引起的二次源信息,是研究地下爆炸震源機(jī)制及檢測(cè)、識(shí)別手段的重要震相?;趯?duì)某一區(qū)域介質(zhì)結(jié)構(gòu)的了解,將實(shí)測(cè)Lg波譜特征形態(tài)與理論波譜形態(tài)進(jìn)行比對(duì),利用Lg波頻譜曲線的低谷點(diǎn)隨震源深度“左移”的特征,作為輔助識(shí)別手段估計(jì)地下核爆炸裝置的埋深。研究結(jié)果對(duì)深入地理解地下核爆炸震源機(jī)制有重要的參考意義。

    Lg波; 波譜低谷點(diǎn); 譜比

    0 引言

    如何從觀測(cè)資料出發(fā)來區(qū)分地下核爆炸與天然地震是我們所關(guān)心的重要問題, 自從實(shí)施地下核爆炸以來,基于地震學(xué)方法的識(shí)別手段得到了較為深入、廣泛的研究,提出了許多種方法。然而大量的觀測(cè)資料分析和理論研究結(jié)果表明,尋找一個(gè)理想的鑒別判據(jù)存在許多困難[1-2],理論上來講,不存在一個(gè)理想的能夠逃避核查的試驗(yàn)方法。迄今為止,地震學(xué)方法仍然是在遠(yuǎn)區(qū)檢測(cè)、識(shí)別他方隱蔽地下核試驗(yàn)的主要手段。

    天然地震與地下核爆炸的本質(zhì)區(qū)別是震源機(jī)制的差異,兩者的震源矩張量解具有完全不同的表現(xiàn)形式[3-4]。地下核爆炸源輻射出來的地震波形較天然地震源輻射出來的波形而言比較簡(jiǎn)單,可通過臺(tái)陣的聚束技術(shù)來分析;但隨著地下核試驗(yàn)記錄的增多,關(guān)于地下核爆炸的復(fù)雜波形的報(bào)道也逐漸增多[5],使得基于波形復(fù)雜性的判據(jù)的實(shí)用性受到普遍懷疑。

    震源深度是識(shí)別地下核爆炸與天然地震的重要指標(biāo),針對(duì)地震核查問題,通常利用譜分析手段來提高源深度的估計(jì)精度。地震波的頻譜特性不僅與震源有關(guān),與傳播路徑也有相當(dāng)大的關(guān)系。即對(duì)一個(gè)核試驗(yàn)場(chǎng)適用的頻譜判據(jù),對(duì)于另外一個(gè)地區(qū)常常不適用,因此對(duì)于精確的介質(zhì)模型,利用頻譜來估計(jì)震源深度是十分重要的。

    導(dǎo)波以及面波是由地球介質(zhì)的豎向幾何結(jié)構(gòu)特征決定的,對(duì)傳播路徑的橫向非均勻性的敏感程度相對(duì)體波來說要弱得多。地下核爆炸激發(fā)的Lg波在遠(yuǎn)區(qū)記錄中穩(wěn)定且占優(yōu)勢(shì),利用其研究地下核爆炸的震源機(jī)制以及識(shí)別手段具有獨(dú)特優(yōu)勢(shì)。理論上來講,球?qū)ΨQ、短時(shí)的爆炸源在分層地球介質(zhì)中激發(fā)的Lg波相對(duì)天然地震具有“簡(jiǎn)單性”,但是地下核爆炸過程的復(fù)雜性仍然使得Lg波復(fù)雜化,伴隨地下爆炸源區(qū)的構(gòu)造應(yīng)力釋放具有“特殊性”,如經(jīng)過CLVD源調(diào)制的Rg波波譜具有低谷點(diǎn)特征等,這些特征是研究地下核爆炸震源機(jī)制的重要手段。有學(xué)者認(rèn)為地下核爆炸的Lg波源自于近源區(qū)的Rg波[6-8],因?yàn)長(zhǎng)g波波譜也具有低谷點(diǎn)特征,且低谷點(diǎn)頻率隨震源深度變化規(guī)律與Rg波一致,所以Lg波不僅攜帶地球豎向分層結(jié)構(gòu)信息,而且攜帶有伴隨爆炸過程的二次源信息,其波譜中的低谷點(diǎn)普遍被認(rèn)為源自于伴隨CLVD源過程。基于上述研究,本文首先分析Lg波激發(fā)與地球介質(zhì)結(jié)構(gòu)的關(guān)系,然后利用Lg波的特殊性,即其波譜低谷點(diǎn)特征,估計(jì)爆炸裝置的埋深。

    1 不同介質(zhì)模型下Lg波的激發(fā)特征

    形成地下核爆炸地震波Lg波所需的S波主要源于爆炸源引起的二次源,同時(shí)也受到源區(qū)介質(zhì)結(jié)構(gòu)的影響,本文利用理論地震圖方法來說明這一觀點(diǎn)。

    層裂源的表現(xiàn)過程為:通常情況下地下核爆炸會(huì)導(dǎo)致源上方介質(zhì)的破碎,并上拋,該物理過程可用張裂源(TCACK)模型進(jìn)行描述;上拋物的回落會(huì)引起應(yīng)力的重新分布,這種構(gòu)造應(yīng)力釋放可以用逆傾滑形式的斷層來描述,表現(xiàn)形狀為頂點(diǎn)位于爆炸源點(diǎn)的反向圓錐體,錐體低面位于層裂面甚至自由表面[9]。這種形式的構(gòu)造應(yīng)力釋放源可以用線性矢量偶極補(bǔ)償源(CLVD)來描述[10]。因此層裂源可以表現(xiàn)為張裂源和線性矢量偶極補(bǔ)償源的組合。Day用單力和張裂模型來描述層裂過程,結(jié)果表明,在頻段0~5 Hz兩種描述具有等價(jià)性[11]。

    與 CLVD源激發(fā)的近源區(qū)Rg波譜相同,Lg波譜中也具有低谷點(diǎn)特征,且其隨頻率變化的規(guī)律與近源區(qū)Rg波譜低谷點(diǎn)隨頻率變化的規(guī)律相同,有學(xué)者對(duì)此也進(jìn)行了深入研究[8,12]。本文采用適合東哈薩克斯坦區(qū)域的地球介質(zhì)速度模型(模型1)[13],即地表高速模型,并對(duì)該模型進(jìn)行修改,使其分別變?yōu)榈乇淼退倌P?模型 2),及含低速層模型(模型 3)。利用理論地震圖方法,進(jìn)一步研究單力源(SIN)、爆炸源(EXP)、張裂源(TCACK)及線性矢量偶極補(bǔ)償源(CLVD)在不同地球介質(zhì)結(jié)構(gòu)中激發(fā)的區(qū)域震相特征,分析不同源模型激發(fā)的Lg波能量與介質(zhì)模型及源模型之間的關(guān)系,以期更好地理解Lg波的激發(fā)機(jī)理。

    圖1 S波和P波在3種模型下的速度Fig.1 The velocities of S-wave & P-wave in the three models

    理論上講,球?qū)ΨQ的爆炸源只激發(fā)P波,形成Lg波的S波只能源自P波在地表的轉(zhuǎn)換波P-S波或由面波散射而形成的S波,但這樣的激發(fā)機(jī)制對(duì)速度模型較敏感。Lg波的能量依賴于爆炸源處P波速度與上地幔S波速度的比值[14],如果該比值大于1,P-S波則不會(huì)陷俘于地殼中形成導(dǎo)波。而實(shí)際觀測(cè)結(jié)果表明,像東哈薩克斯坦這樣的地表高速地殼結(jié)構(gòu),地下核爆炸激發(fā)的Lg波同樣是區(qū)域震相中的優(yōu)勢(shì)成分,說明激發(fā)Lg波的主要因素并不是球?qū)ΨQ的爆炸源,而是與非爆炸源成分息息相關(guān),或者說與球?qū)ΨQ爆炸源引起的二次源有關(guān)。

    模型2 中地表P波速度低于上地幔S波速度,理論上講,純爆炸源在這樣的速度模型中能激發(fā)能量較大的Lg波;模型3含有一低速層,理論上來講,該低速層利于導(dǎo)波的形成。

    本文采用水平層狀地球模型中理論地震圖的計(jì)算方法[15-16],其本質(zhì)與Yao[17]的經(jīng)典的廣義反射-透射系數(shù)方法是一致的,具有能夠模擬出全波場(chǎng)震相的優(yōu)點(diǎn)。分別計(jì)算SIN、EXP、TCACK及CLVD源在三種模型下的區(qū)域理論地震圖,震中距為300 km,震源時(shí)間函數(shù)均取δ函數(shù),理論計(jì)算時(shí)頻率范圍約束為0~5.0 Hz。數(shù)值模擬結(jié)果如圖2所示。

    圖2分別為SIN、EXP、TCACK及CLVD源在三種模型下的格林函數(shù)??梢钥闯?,CLVD源在三種速度模型下都能激發(fā)能量足夠大的Lg波,而 EXP源在相應(yīng)模型下激發(fā)的Lg波能量要遠(yuǎn)低于CLVD源情況,尤其是在模型1、模型3情況下??梢奅XP源在地表高速模型下,地殼中會(huì)耦合能量很小的S波,而相應(yīng)的CLVD源卻能產(chǎn)生能量很大的用來激發(fā)Lg波的S波。究其原因可以從CLVD源與EXP源空間分布特征得到解釋,確切地說CLVD源對(duì)角分量的偏差直接激發(fā)S波,其速度要低于上地幔S波速,會(huì)在Moho反射形成能量足夠大的導(dǎo)波。EXP源在模型2會(huì)激發(fā)一定能量的Lg波,因?yàn)樵撃P蜑榈乇淼退倌P?,地表轉(zhuǎn)換波P-S波能夠在Moho面反射、疊加形成導(dǎo)波。

    TCACK源也具有對(duì)角分量偏差,相比CLVD源小,因此激發(fā)的Lg波能量較小,但與SIN源激發(fā)的能量相當(dāng),在所考察頻率范圍內(nèi)。支持Day的SIN和TCACK源在描述層裂物理過程中具有一致性觀點(diǎn)[11]。

    從上述分析結(jié)果可以看出,無論是在地表高速、地表低速類型的地殼速度結(jié)構(gòu)的區(qū)域,相對(duì)EXP源而言,CLVD源激發(fā)Lg波對(duì)介質(zhì)模型不敏感,并且與實(shí)際地下核爆炸激發(fā)的區(qū)域震相情況相符。由此可以證明伴隨實(shí)際地下核爆炸的CLVD源擔(dān)當(dāng)了激發(fā)Lg波的主要任務(wù)。

    圖2 三種源在三種模型下的數(shù)值模擬結(jié)果(速度場(chǎng))Fig.2 Numerical simulations of three sources in the three models (velocity field)

    2 波譜形態(tài)比對(duì)

    2.1 譜比方法

    震源機(jī)制和傳播路徑對(duì)地下核爆炸激發(fā)的地震波起著決定性作用。 對(duì)地下核爆炸來說,除了球?qū)ΨQ爆炸源外,需要考慮其他輔助震源模式, 如伴隨層裂過程的TCACK、CLVD源等;傳播路徑影響是造成地下核爆炸地震波波形復(fù)雜化的主要原因,且難以用數(shù)學(xué)形式精確描述,想要從波形資料中提取出感興趣的震源信息,如Lg波中與CLVD源相關(guān)信息,必須消除路徑的影響。利用經(jīng)驗(yàn)格林函數(shù)方法可以實(shí)現(xiàn)這一目的[18-19]:用過比例埋深的地下核爆炸事件作為格林函數(shù),所謂過比例埋深是指爆炸裝置的埋深相對(duì)其當(dāng)量來說大于正常埋深的情況。過比例埋深的地下核爆炸的層裂效應(yīng)很小,在0~3 Hz的頻帶范圍內(nèi)可近似為純爆炸源,它在波譜比中的作用相當(dāng)于經(jīng)驗(yàn)格林函數(shù),用以消除傳播路徑、接受場(chǎng)地及記錄儀器響應(yīng)的影響。因此波譜比給出的結(jié)果主要是震源的影響因素造成的。

    正常埋深地下核爆炸的遠(yuǎn)場(chǎng)位移譜(垂向)可以描述為:

    (1)

    式中:冒號(hào)表示內(nèi)積計(jì)算;Me為爆炸源矩張量;he和hx分別為爆炸源和層裂源的中心深度;G為垂向分量位移場(chǎng)對(duì)應(yīng)的格林函數(shù);Mx為描述層裂源的矩張量。單力模型下,層裂產(chǎn)生的位移場(chǎng)可以描述為:

    (2)

    式中:Fs是σs的面積分結(jié)果。σs為地球表面處的垂向牽引力。在點(diǎn)矩張量模型下,并在波長(zhǎng)遠(yuǎn)大于層裂埋深的條件下,層裂過程中的張裂源的矩張量描述與點(diǎn)力描述在數(shù)學(xué)上是等價(jià)的。

    相應(yīng)地,作為經(jīng)驗(yàn)格林函數(shù)的過比例埋深的地下核爆炸的遠(yuǎn)場(chǎng)位移(垂向)譜可描述為:

    (3)

    (4)

    針對(duì)適合東哈薩克斯坦區(qū)域的地殼速度模型[13],分別計(jì)算EXP、CLVD源的格林函數(shù),所采用參數(shù)同上。CLVD源深度分別取0.2 km、0.3 km、0.6 km、0.9 km,提取Lg波,將EXP的結(jié)果作為格林函數(shù),進(jìn)行式(4)第二項(xiàng)的運(yùn)算。理論譜比結(jié)果如圖3中實(shí)線所示(縱坐標(biāo)已歸一化),顯示出不同震源深度低谷點(diǎn)隨頻率位置的變化規(guī)律。理論上來講,該變化規(guī)律反映了東哈薩克斯坦區(qū)域地下爆炸激發(fā)的Lg波譜低谷點(diǎn)形態(tài)與CLVD源深度之間的關(guān)系,即隨源深度增加低谷點(diǎn)頻率向低頻方向移動(dòng)。這一變化關(guān)系可以用來估計(jì)爆炸源深度。

    2.2 實(shí)際應(yīng)用

    分別對(duì)我國(guó)烏魯木齊臺(tái)站記錄到的兩次哈薩克斯坦地下核爆炸事件(代號(hào)為1988044和1987157)進(jìn)行分析,震中距約950km,震級(jí)分別Mb=6.1、Mb=5.0。記錄波形如圖3所示,為發(fā)育較好的Lg波。將震級(jí)較小事件1987157作為經(jīng)驗(yàn)格林函數(shù),反卷積震級(jí)較大事件1988044的Lg波,將結(jié)果疊加于理論譜比圖4中,與理論值進(jìn)行比對(duì),虛線為1988044事件的譜比結(jié)果。可以看出,譜比值具有明顯的低谷點(diǎn)特征,且其頻率位置與CLVD源深度為0.2km的理論譜比值一致。

    圖3 東哈薩克斯坦地下核爆炸垂向速度場(chǎng)地震波形(單位:counts)Fig.3 Seismic waves in vertical velocity field caused by underground nuclear explosions of East Kazakh (unit:counts)

    圖4 針對(duì)東哈薩克斯坦地殼速度模型的Lg波譜比Fig.4 The Lg-wave spectral ratio calculated by crustal velocity model of East Kazakh

    利用CLVD源深度與低谷點(diǎn)的關(guān)系式HCLVD=V/16fNULL[6],可以得出,低谷點(diǎn)頻率范圍為0.94~1.87 Hz的地下核爆炸的爆炸裝置埋深范圍為0.5~1.0 km,這里取v=5.0 km/s,絕大多數(shù)震級(jí)為Mb=6.0左右的大當(dāng)量地下核爆炸裝置的埋深都在此范圍內(nèi),與圖3理論結(jié)果一致。CLVD源深度約為爆炸源深度的1/3[6],據(jù)此可以判斷出事件1988044地下核爆炸的爆炸裝置的埋深約為0.6 km,與利用HCLVD=V/16fNULL估算出的結(jié)果一致,說明利用適合東哈薩克斯坦地區(qū)的地球介質(zhì)模型得到的地下核爆炸理論Lg波譜比,可以用來估計(jì)實(shí)際地下爆炸裝置的埋深。

    3 結(jié)論

    對(duì)重要震源機(jī)制進(jìn)行相對(duì)深入的研究,是逐步了解地下核爆炸震源物理過程的比較現(xiàn)實(shí)的方法。層裂源輔助源的重要性在于它是區(qū)域震相Lg波的主要激發(fā)因素,對(duì)Lg波具有調(diào)制作用,表現(xiàn)在Lg波譜低谷點(diǎn)的存在。低谷點(diǎn)位置隨頻率變化的客觀事實(shí),是由震源深度及介質(zhì)幾何分層結(jié)構(gòu)特征的影響,所以說Lg波的低谷點(diǎn)同時(shí)攜帶有源和介質(zhì)結(jié)構(gòu)信息。

    本文利用被地下核爆炸二次源“復(fù)雜化了”的Lg波的特殊譜形態(tài),即低谷點(diǎn)特征,估計(jì)爆炸裝置的埋深,將低谷點(diǎn)特征與源深度估計(jì)結(jié)合在一起。從識(shí)別的角度來看,該項(xiàng)工作比單純估計(jì)震源深度更有意義。對(duì)發(fā)生在某一地下核試驗(yàn)區(qū)域的未知事件,如果Lg波譜比具有明顯的低谷點(diǎn),且深度估計(jì)范圍在1 km以內(nèi),則可以作為判斷該事件傾向于地下爆炸的一個(gè)重要的輔助判據(jù)。

    References)

    [1] Blandford R R.Discrimination Between Earthquakes and Underground Explosions[J].Ann Rev Earth Planet Soc,1977,5:111-122.

    [2] Taylor S R,Denny M.Regional Discrimination Between NTS Explosion Western U.S.Earthquakes[J].Bull Seism Soc Am,1989,79:1142-1176 .

    [3] Minson S E,Dreger D S.STableInversions for Complete Moment Tensors[J].Geophysical Journal International,2008,174:585-592.

    [4] 陳繼鋒,林向東,何新社.2013年甘肅岷縣MS6.6地震矩張量反演及發(fā)震構(gòu)造初探[J].地震工程學(xué)報(bào),2013,35(3):425-431.CHEN Ji-feng,LIN Xiang-dong,HE Xin-she.Moment Tensor Inversion and Seismogenic Tectionics of the 2013 MinxianMS6.6 Earthquake in Gansu[J].China Earthquake Engineering Journal,2013,35(3):425-431.(in Chinese)

    [5] 何永鋒,趙克常,張獻(xiàn)兵,等.地下核爆炸地震波二次源特征[J].地球物理學(xué)報(bào),2012,55(5):1742-1748.HE Yong-feng,ZHAO Ke-chang,ZHANG Xian-bing,et al.The Characteristic of the Waveform from the Second Source Induced by Underground Explosion[J].Chinese Journal of Geophysics,2012,55(5):1742-1748.(in Chinese)

    [6] Patton H J,Taylor R S.Analysis of Lg Spectral Ratios From NTS Explosions:Implications for the Source Mechanisms of Spall and Generation of Lg Waves[J].Bull Seism Soc Am,1995,85:220-236.

    [7] Gupta I N,Zhang T R,Wangner R A.Low-Frequency Lg from NTS and Kazakh Nuclear Explosions-observations and Interpretation[J].Bull Seism Soc Am,1997,87:1115-1125.

    [8] 何永鋒,陳曉非,何耀峰,等.地下爆炸Rg波低谷點(diǎn)激發(fā)機(jī)理[J].地球物理學(xué)報(bào),2005,48(3):643-648.HE Yong-feng,CHEN Xiao-fei,HE Yao-feng,et al.Generation of Null in Rg Wave by Underground Explosions[J].Chinese Journal of Geophysics,2005,48(3): 643-648.(in Chinese)

    [9] Masse R P.Review of Seismic Source Models for Underground Nuclear Explosions[J].Bull Seism Soc Am,1981,71:1249-1268.

    [10] Knopoff L,Randall M J.The Compensated Linear Vector Dipole:a Possible Mechanism for Deep Earthquake[J].J Geophys Res,1970,75:4975-4963.

    [11] Day S M,Mclaughlin K L.Seismic Source Representations for Spall[J].Bull Seism Soc Am,1991,81(1):191-201.

    [12] 何永鋒,趙克常,張獻(xiàn)兵.地下核爆炸的主要非爆炸源機(jī)制[J].地球物理學(xué)進(jìn)展,2010,25(3):789-794.HE Yong-feng,ZHAO Ke-chang,ZHANG Xian-bing.Primary Non-explosive Source Mechanisms for Underground Nuclear Explosions[J].Progress in Geophysics,2010,25(3):789-794.(in Chinese)

    [13] McLaughlin K L,Barker T G,Day S M.Effects of Depth of Burial on Explosion and Earthquake Regional Seismograms:Regional Discrimination and Yield Estimation[R].S_CUBED report SSS_R_88_9844,Jolla,California,1988.

    [14] Xie X B,Lay T.The Excitation of Lg Wave by Bxplosions:a Finite-difference Investigation[J].Bull Seism Soc Am,1994,84:324-342.

    [15] Chen X F.A Systematic and Efficient Method of Computing Normal Modes for Multilayered Half-space[J].Geophysical Journal International,1993,115:391-409.

    [16] Chen X F.Seismogram Synthesis in Multi-Layered Half-space Part I[J].Theoretical Formulations,Earthquake Research in China,1999,13:149-174.

    [17] Yao Z X,Harkrider D G.A Generalized Reflection-transmission Coefficient Matrix and Discrete Wave Nnumber Method for Synthetic Seismograms[J].Bull Seismol Soc Am,1983,73:1685-1699.

    [18] Gupta I N,Chan W W,Wanger R A.A Comparison of Regional Phases from Underground Nuclear Explosions at East Kazakh and Nevada Test Sites[J].Bull Seism Soc Am,1992,82:352-382

    [19] He Yong-feng,Chen Xiao-fei,Zhang Hai-ming.The Effect of Spall on Lg Waves[J].Aata Seismological Sinica,2005,18(3):331-336.

    Estimation of Underground Nuclear Explosive Device Depth Using Lg-wave Spectrum

    HE Yong-feng, LI Kai, LIU Bing-can

    (AcademyofArmoredForcesEngineering,Beijing100072,China)

    The efficient frequency-wave number code could synthesize high-frequency (10 Hz) regional seismograms (up to a distance of more than 1 000 km) in a medium consisting of a large number of crustal layers,and investigate the effects of regional wave guides on Lg-waves.Lg-waves contain source information and the effect of the near-source structure.The scattering of explosion-generated Rg into S appears to be the primary contributor to the low-frequency Lg from nuclear explosive devices.For three typical crustal models,the synthetic seismograms calculated for EXP and Spall sources show that the generation mechanism of Lg-waves can be explained by a compensated linear-vector dipole source.If the velocity structure is well known,then the Lg-waveform can be synthesized,so we can compare the theoretical Lg-wave’s spectrum with the observation’s spectrum,and the results can not only be used to estimate the depth of the underground explosion equipment but also a discriminant.A deep understanding of the mechanism of the generation of Lg-waves can help us to create a new discriminant.

    Lg-wave; spectral null; spectral ratios

    2014-07-09

    國(guó)家自然科學(xué)基金項(xiàng)目(40874023, 41374068)

    何永鋒(1966-),男,教授,主要研究方向?yàn)楹吮ǖ卣饘W(xué).E-mail:heyfeng@sina.com

    P315

    A

    1000-0844(2015)02-0601-05

    10.3969/j.issn.1000-0844.2015.02.0601

    猜你喜歡
    核爆炸低谷波譜
    低空核爆炸環(huán)境效應(yīng)模擬研究
    用矩陣分解方法識(shí)別地下核爆炸
    核安全及其技術(shù)的發(fā)展
    越在低谷時(shí),越不要閑著
    基于核爆炸的小行星偏轉(zhuǎn)方案分析
    民族音樂處在“低谷”
    2019平淡不平靜 2020低迷非低谷
    善待低谷,才能迎來翻盤
    海峽姐妹(2018年2期)2018-04-12 07:13:19
    琥珀酸美托洛爾的核磁共振波譜研究
    美國(guó)波譜通訊系統(tǒng)公司
    成人欧美大片| 在线免费观看的www视频| 高清av免费在线| 精品免费久久久久久久清纯| 精品人妻熟女av久视频| 99热这里只有是精品50| 国产精品.久久久| 国产精品福利在线免费观看| 乱人视频在线观看| 国产人妻一区二区三区在| 别揉我奶头 嗯啊视频| 欧美激情国产日韩精品一区| 欧美日韩精品成人综合77777| 国产精品日韩av在线免费观看| av国产久精品久网站免费入址| 特大巨黑吊av在线直播| 看免费成人av毛片| 日日啪夜夜撸| 亚洲av电影在线观看一区二区三区 | 国产午夜精品一二区理论片| 亚洲电影在线观看av| www日本黄色视频网| 乱人视频在线观看| 干丝袜人妻中文字幕| 成人午夜精彩视频在线观看| 麻豆久久精品国产亚洲av| 久久久久久伊人网av| 美女xxoo啪啪120秒动态图| 久久久色成人| 久久99热6这里只有精品| 精品99又大又爽又粗少妇毛片| 国产在线一区二区三区精 | 少妇被粗大猛烈的视频| 99国产精品一区二区蜜桃av| 国产69精品久久久久777片| 久久精品熟女亚洲av麻豆精品 | 熟女电影av网| av免费观看日本| 久久久精品欧美日韩精品| 一本一本综合久久| 啦啦啦观看免费观看视频高清| 亚洲国产高清在线一区二区三| 最新中文字幕久久久久| 国产精品久久久久久av不卡| 又爽又黄a免费视频| 国产黄片美女视频| 美女cb高潮喷水在线观看| 国产精品国产三级国产专区5o | 激情 狠狠 欧美| 在线免费十八禁| 在线播放国产精品三级| 99久久人妻综合| 在线播放国产精品三级| 午夜视频国产福利| 纵有疾风起免费观看全集完整版 | 激情 狠狠 欧美| 国产老妇伦熟女老妇高清| 国产色爽女视频免费观看| 一个人看的www免费观看视频| 亚洲乱码一区二区免费版| 亚洲内射少妇av| 亚洲国产色片| 国产精品野战在线观看| 人人妻人人澡欧美一区二区| 搡老妇女老女人老熟妇| 免费黄网站久久成人精品| 人妻制服诱惑在线中文字幕| 亚洲成人中文字幕在线播放| 免费搜索国产男女视频| 久久久久久久国产电影| 99视频精品全部免费 在线| 国内揄拍国产精品人妻在线| 午夜福利在线观看吧| 久久精品久久久久久久性| 91狼人影院| 国产极品天堂在线| 最近中文字幕高清免费大全6| 一级黄片播放器| 黄色欧美视频在线观看| 国产亚洲5aaaaa淫片| a级毛片免费高清观看在线播放| 国产成人a∨麻豆精品| 久久精品国产亚洲av天美| 免费av不卡在线播放| 草草在线视频免费看| 少妇猛男粗大的猛烈进出视频 | 一级av片app| 亚洲欧美日韩东京热| av国产免费在线观看| 女人十人毛片免费观看3o分钟| 精品99又大又爽又粗少妇毛片| 伦理电影大哥的女人| 成人国产麻豆网| 午夜福利在线观看吧| 国产 一区 欧美 日韩| 熟女电影av网| 男女国产视频网站| 亚洲久久久久久中文字幕| 色综合色国产| 精品午夜福利在线看| 亚洲国产精品成人久久小说| 日日摸夜夜添夜夜添av毛片| 在线观看一区二区三区| 免费搜索国产男女视频| videos熟女内射| 国产一级毛片七仙女欲春2| 国产av不卡久久| 久99久视频精品免费| 秋霞伦理黄片| eeuss影院久久| 久久韩国三级中文字幕| 联通29元200g的流量卡| 久久精品人妻少妇| 免费观看人在逋| 国产精品蜜桃在线观看| 中文亚洲av片在线观看爽| 国产精品一区二区三区四区久久| 青春草国产在线视频| 搞女人的毛片| 美女黄网站色视频| 久久久久性生活片| 国产黄片美女视频| 国产亚洲午夜精品一区二区久久 | 国产精品嫩草影院av在线观看| 国产成人精品一,二区| 成人毛片60女人毛片免费| 久久精品综合一区二区三区| 国产91av在线免费观看| 1000部很黄的大片| 国产高清有码在线观看视频| 又黄又爽又刺激的免费视频.| av在线蜜桃| 97在线视频观看| 色视频www国产| 久久鲁丝午夜福利片| 国产成人91sexporn| 在现免费观看毛片| 国产成人一区二区在线| 国产精品久久久久久精品电影| 日本一二三区视频观看| 色5月婷婷丁香| 婷婷色综合大香蕉| 国产人妻一区二区三区在| 中国美白少妇内射xxxbb| 少妇被粗大猛烈的视频| 乱系列少妇在线播放| 国产成人一区二区在线| 高清日韩中文字幕在线| 国产乱来视频区| 免费看光身美女| 日本wwww免费看| 99久久精品国产国产毛片| 日韩一区二区视频免费看| 免费看a级黄色片| 免费观看性生交大片5| 久久精品国产亚洲av天美| 亚洲中文字幕一区二区三区有码在线看| 欧美一区二区国产精品久久精品| 在线播放国产精品三级| 汤姆久久久久久久影院中文字幕 | 美女cb高潮喷水在线观看| 国产伦精品一区二区三区四那| 国产三级在线视频| 色播亚洲综合网| 国产色婷婷99| 少妇人妻一区二区三区视频| 国产老妇伦熟女老妇高清| 97超碰精品成人国产| 精品国产三级普通话版| 国内精品宾馆在线| 日韩,欧美,国产一区二区三区 | 久久亚洲国产成人精品v| 男女下面进入的视频免费午夜| 国产精华一区二区三区| 18+在线观看网站| 日本色播在线视频| 亚洲国产色片| 看十八女毛片水多多多| 丰满乱子伦码专区| 高清在线视频一区二区三区 | 一二三四中文在线观看免费高清| 久久精品国产亚洲av涩爱| 亚洲精品,欧美精品| 视频中文字幕在线观看| 精品久久久久久久久久久久久| 国产av码专区亚洲av| 97超视频在线观看视频| 国产亚洲一区二区精品| 三级国产精品欧美在线观看| 亚洲欧美日韩卡通动漫| 国产伦精品一区二区三区四那| 免费观看a级毛片全部| 午夜老司机福利剧场| 看黄色毛片网站| 国产一区二区在线观看日韩| 热99re8久久精品国产| 女人久久www免费人成看片 | 一个人看的www免费观看视频| 又粗又硬又长又爽又黄的视频| 美女高潮的动态| 国产高清不卡午夜福利| 97超视频在线观看视频| 一级黄片播放器| av在线亚洲专区| 成人国产麻豆网| 亚洲国产色片| 精品国产一区二区三区久久久樱花 | 少妇猛男粗大的猛烈进出视频 | 中文字幕制服av| 亚洲中文字幕一区二区三区有码在线看| 我要搜黄色片| 中文字幕久久专区| 国产成人91sexporn| 韩国av在线不卡| 欧美最新免费一区二区三区| 综合色丁香网| 国产激情偷乱视频一区二区| 色综合亚洲欧美另类图片| 韩国av在线不卡| 精品不卡国产一区二区三区| 校园人妻丝袜中文字幕| 日韩欧美在线乱码| 七月丁香在线播放| 国产v大片淫在线免费观看| 1000部很黄的大片| 性插视频无遮挡在线免费观看| 国产精品野战在线观看| 国产一区二区在线观看日韩| 99热6这里只有精品| 看十八女毛片水多多多| 日本av手机在线免费观看| 亚洲av成人精品一区久久| 日本与韩国留学比较| 亚洲最大成人手机在线| 色综合色国产| 久久久精品欧美日韩精品| 国产 一区 欧美 日韩| 男人舔奶头视频| 好男人视频免费观看在线| 免费一级毛片在线播放高清视频| 国产私拍福利视频在线观看| 免费黄网站久久成人精品| 91午夜精品亚洲一区二区三区| 18禁在线无遮挡免费观看视频| 亚洲av中文字字幕乱码综合| 中文精品一卡2卡3卡4更新| 婷婷色av中文字幕| 在线观看av片永久免费下载| 在现免费观看毛片| 亚洲在久久综合| 在线免费观看的www视频| 超碰97精品在线观看| 大话2 男鬼变身卡| 久久午夜福利片| 日日摸夜夜添夜夜添av毛片| 成年女人永久免费观看视频| 午夜精品一区二区三区免费看| 色噜噜av男人的天堂激情| 欧美xxxx性猛交bbbb| 中文亚洲av片在线观看爽| 日韩成人av中文字幕在线观看| 午夜精品一区二区三区免费看| 男人舔女人下体高潮全视频| 午夜激情欧美在线| 人妻夜夜爽99麻豆av| 色播亚洲综合网| 国产精品野战在线观看| av又黄又爽大尺度在线免费看 | av在线观看视频网站免费| 美女内射精品一级片tv| 波多野结衣巨乳人妻| 国产成人免费观看mmmm| 午夜精品国产一区二区电影 | 久久久久免费精品人妻一区二区| 久久精品影院6| 免费黄网站久久成人精品| 日产精品乱码卡一卡2卡三| 春色校园在线视频观看| 夫妻性生交免费视频一级片| 色播亚洲综合网| 亚洲无线观看免费| 久久精品影院6| 直男gayav资源| 91av网一区二区| 国产黄色视频一区二区在线观看 | 中文字幕精品亚洲无线码一区| 欧美日韩综合久久久久久| 日本爱情动作片www.在线观看| 看片在线看免费视频| 色综合亚洲欧美另类图片| 免费看日本二区| 亚洲精品日韩在线中文字幕| 久99久视频精品免费| av.在线天堂| 欧美zozozo另类| 亚洲国产日韩欧美精品在线观看| 久久久久久大精品| 国产精品电影一区二区三区| 日韩三级伦理在线观看| 天堂影院成人在线观看| www.色视频.com| 欧美成人免费av一区二区三区| 久久久亚洲精品成人影院| 一本久久精品| 亚洲中文字幕日韩| 最后的刺客免费高清国语| 久久精品国产亚洲av涩爱| 淫秽高清视频在线观看| 波多野结衣巨乳人妻| 亚洲精品色激情综合| 内地一区二区视频在线| 精华霜和精华液先用哪个| 亚洲人与动物交配视频| 久久久亚洲精品成人影院| 日本欧美国产在线视频| 午夜福利在线观看免费完整高清在| 99久久无色码亚洲精品果冻| 看免费成人av毛片| 69av精品久久久久久| 中文资源天堂在线| 亚洲国产欧美人成| 国产欧美日韩精品一区二区| 亚洲国产欧洲综合997久久,| 中文字幕久久专区| 有码 亚洲区| 成人美女网站在线观看视频| 亚洲国产欧洲综合997久久,| 免费观看精品视频网站| 亚洲国产欧洲综合997久久,| 亚洲在线观看片| 丰满人妻一区二区三区视频av| 久久精品国产亚洲av天美| 色噜噜av男人的天堂激情| 女人被狂操c到高潮| 蜜臀久久99精品久久宅男| 精品人妻熟女av久视频| 亚洲精品乱码久久久久久按摩| 2021少妇久久久久久久久久久| 欧美区成人在线视频| 久久韩国三级中文字幕| 亚洲最大成人中文| 卡戴珊不雅视频在线播放| 亚洲欧美日韩东京热| 午夜亚洲福利在线播放| 内射极品少妇av片p| 亚洲国产最新在线播放| 熟女电影av网| 午夜精品在线福利| 三级国产精品欧美在线观看| 国产精品一区www在线观看| 亚洲av日韩在线播放| 插阴视频在线观看视频| 精品国产露脸久久av麻豆 | 人妻少妇偷人精品九色| 国产爱豆传媒在线观看| 午夜爱爱视频在线播放| 男的添女的下面高潮视频| 精品一区二区三区视频在线| ponron亚洲| 国产精品一区www在线观看| 亚洲第一区二区三区不卡| 午夜免费激情av| 欧美日韩综合久久久久久| 亚洲五月天丁香| 亚洲经典国产精华液单| 亚洲国产最新在线播放| 午夜福利在线观看免费完整高清在| 国产在视频线精品| 亚洲精品乱久久久久久| 国产一级毛片在线| 小蜜桃在线观看免费完整版高清| 免费黄色在线免费观看| 国产白丝娇喘喷水9色精品| 舔av片在线| 亚洲丝袜综合中文字幕| 久久久久精品久久久久真实原创| 日韩 亚洲 欧美在线| 国产精品.久久久| 97在线视频观看| 最近最新中文字幕大全电影3| 黄片无遮挡物在线观看| 国产精品久久久久久精品电影小说 | 能在线免费看毛片的网站| 亚洲av中文字字幕乱码综合| 亚洲美女视频黄频| 国产麻豆成人av免费视频| 国产伦理片在线播放av一区| 草草在线视频免费看| 日韩成人av中文字幕在线观看| 国产日韩欧美在线精品| 久久亚洲精品不卡| 99国产精品一区二区蜜桃av| 免费观看精品视频网站| 九九久久精品国产亚洲av麻豆| 成人鲁丝片一二三区免费| 午夜福利高清视频| 好男人视频免费观看在线| 伊人久久精品亚洲午夜| 亚洲av电影在线观看一区二区三区 | 国产精品一区二区性色av| 欧美性猛交╳xxx乱大交人| 国产高清三级在线| 五月伊人婷婷丁香| 狂野欧美白嫩少妇大欣赏| 夜夜看夜夜爽夜夜摸| 欧美性猛交黑人性爽| 欧美97在线视频| 国产精品久久久久久精品电影| 男插女下体视频免费在线播放| 日韩欧美在线乱码| 精品国产三级普通话版| 麻豆av噜噜一区二区三区| 在线天堂最新版资源| 成人午夜高清在线视频| 精品久久久噜噜| 五月伊人婷婷丁香| 亚洲欧美日韩高清专用| 亚洲最大成人中文| 日本猛色少妇xxxxx猛交久久| 欧美成人午夜免费资源| 中文天堂在线官网| videossex国产| 日韩国内少妇激情av| 国产精品日韩av在线免费观看| 国产乱人偷精品视频| 日本午夜av视频| 精品一区二区三区视频在线| 欧美最新免费一区二区三区| 欧美一区二区亚洲| 国产伦精品一区二区三区四那| 国产精品一及| 国产单亲对白刺激| 日韩人妻高清精品专区| 亚洲高清免费不卡视频| 两性午夜刺激爽爽歪歪视频在线观看| 国产高清视频在线观看网站| 天堂影院成人在线观看| 高清av免费在线| 亚洲欧美成人精品一区二区| 91在线精品国自产拍蜜月| 三级经典国产精品| 欧美色视频一区免费| 日韩制服骚丝袜av| 免费黄色在线免费观看| 亚洲成人中文字幕在线播放| 高清在线视频一区二区三区 | 桃色一区二区三区在线观看| 男女下面进入的视频免费午夜| 国产在视频线在精品| 中文精品一卡2卡3卡4更新| 欧美变态另类bdsm刘玥| 嫩草影院新地址| 老司机影院毛片| 午夜精品在线福利| 神马国产精品三级电影在线观看| 中文字幕制服av| 高清日韩中文字幕在线| 小蜜桃在线观看免费完整版高清| 久久99热这里只频精品6学生 | 久久精品久久精品一区二区三区| 小说图片视频综合网站| 亚洲欧美一区二区三区国产| 亚洲成人中文字幕在线播放| 观看免费一级毛片| 日韩强制内射视频| 国产精品一区二区三区四区免费观看| 亚洲最大成人手机在线| 汤姆久久久久久久影院中文字幕 | 成人美女网站在线观看视频| 国产高清国产精品国产三级 | 好男人在线观看高清免费视频| 22中文网久久字幕| 又爽又黄无遮挡网站| 91av网一区二区| 欧美bdsm另类| 超碰97精品在线观看| 国产av在哪里看| 久久精品夜色国产| 亚洲av中文av极速乱| 蜜桃亚洲精品一区二区三区| 七月丁香在线播放| 欧美高清成人免费视频www| 午夜亚洲福利在线播放| 国产精品久久久久久久电影| 午夜激情欧美在线| 亚洲av成人精品一二三区| 亚洲在线观看片| 我要看日韩黄色一级片| 亚洲色图av天堂| 精品欧美国产一区二区三| 波多野结衣高清无吗| 午夜福利成人在线免费观看| 国产亚洲午夜精品一区二区久久 | 午夜视频国产福利| 久久6这里有精品| 免费观看在线日韩| 国产69精品久久久久777片| 99久久精品热视频| 看片在线看免费视频| 成人午夜高清在线视频| 身体一侧抽搐| 桃色一区二区三区在线观看| 国产精品久久久久久精品电影| 国产精品野战在线观看| 久久6这里有精品| 国产一级毛片七仙女欲春2| 日本黄色片子视频| 七月丁香在线播放| 国产亚洲精品av在线| 亚洲五月天丁香| 亚洲欧美精品自产自拍| 女人被狂操c到高潮| 日韩大片免费观看网站 | 日本wwww免费看| 成人综合一区亚洲| 舔av片在线| .国产精品久久| 水蜜桃什么品种好| 精品久久久久久久久亚洲| 国产色婷婷99| 日本爱情动作片www.在线观看| av.在线天堂| 韩国高清视频一区二区三区| 久久久成人免费电影| 高清在线视频一区二区三区 | 国产日韩欧美在线精品| 免费观看人在逋| 日本免费在线观看一区| 久久久国产成人精品二区| 日韩大片免费观看网站 | 久久久久久久久久黄片| 亚洲在线观看片| 亚洲欧洲日产国产| 久久久久久久亚洲中文字幕| 国产成人精品婷婷| 久久久久久久久久久丰满| 日韩中字成人| 国产乱人偷精品视频| 日韩三级伦理在线观看| 亚洲丝袜综合中文字幕| 男女边吃奶边做爰视频| 国产午夜精品一二区理论片| 国产精品精品国产色婷婷| 91精品一卡2卡3卡4卡| 一个人免费在线观看电影| 人妻制服诱惑在线中文字幕| 听说在线观看完整版免费高清| 亚洲性久久影院| 免费观看性生交大片5| 岛国毛片在线播放| 毛片一级片免费看久久久久| 免费无遮挡裸体视频| 午夜福利在线观看吧| 国产极品精品免费视频能看的| 国产精品久久久久久av不卡| 26uuu在线亚洲综合色| 女的被弄到高潮叫床怎么办| 久久精品久久久久久噜噜老黄 | 亚洲国产精品sss在线观看| 国产v大片淫在线免费观看| 日韩制服骚丝袜av| 亚洲三级黄色毛片| 国产高清有码在线观看视频| 1000部很黄的大片| 亚洲久久久久久中文字幕| 有码 亚洲区| 国产精品国产高清国产av| 亚洲四区av| 亚洲成人中文字幕在线播放| 欧美97在线视频| 一二三四中文在线观看免费高清| 欧美一区二区精品小视频在线| 蜜桃久久精品国产亚洲av| 97在线视频观看| 欧美最新免费一区二区三区| 99热6这里只有精品| 高清日韩中文字幕在线| 特级一级黄色大片| 人妻少妇偷人精品九色| 99久久精品热视频| 伦精品一区二区三区| 一本一本综合久久| 久久精品夜夜夜夜夜久久蜜豆| 中文字幕av在线有码专区| 亚洲综合精品二区| 美女大奶头视频| 国产高清三级在线| 久久久精品欧美日韩精品| 成年女人看的毛片在线观看| 女的被弄到高潮叫床怎么办| 亚洲国产欧美人成| 一个人看视频在线观看www免费| 一区二区三区四区激情视频| 大话2 男鬼变身卡| 亚洲人成网站在线播| 久久久精品94久久精品| 亚洲最大成人中文| 男人的好看免费观看在线视频| 国产极品天堂在线| 国产一区有黄有色的免费视频 | 国模一区二区三区四区视频| av播播在线观看一区| 一级爰片在线观看| 国产成人a∨麻豆精品| 久久99精品国语久久久| 日韩视频在线欧美| 亚洲五月天丁香| 深夜a级毛片| 天堂av国产一区二区熟女人妻| 91精品伊人久久大香线蕉| 午夜福利在线在线| 国产色婷婷99| 18禁在线无遮挡免费观看视频| 国产精品无大码| 国产激情偷乱视频一区二区|