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

    EJ339A型液體閃爍探測(cè)系統(tǒng)光輸出響應(yīng)函數(shù)的理論模擬與實(shí)驗(yàn)研究

    2016-10-13 03:53:44王亭亭莫釗洪趙德山熊忠華
    核技術(shù) 2016年5期
    關(guān)鍵詞:響應(yīng)函數(shù)中子質(zhì)子

    王亭亭 王 亮 莫釗洪 趙德山 熊忠華

    ?

    EJ339A型液體閃爍探測(cè)系統(tǒng)光輸出響應(yīng)函數(shù)的理論模擬與實(shí)驗(yàn)研究

    王亭亭1王 亮2莫釗洪2趙德山2熊忠華2

    1(中國(guó)工程物理研究院 材料研究所表面物理與化學(xué)重點(diǎn)實(shí)驗(yàn)室 江油 621907)2(中國(guó)工程物理研究院材料研究所 綿陽(yáng) 621700)

    對(duì)液體閃爍探測(cè)器EJ339A中子與伽馬輻射測(cè)量問(wèn)題,采用理論模擬與實(shí)驗(yàn)分析方法,結(jié)合自主設(shè)計(jì)搭建的基于LabVIEW的數(shù)字化信號(hào)處理系統(tǒng),分別完成22Na、133Ba、137Cs和60Co等4種不同能量的γ源等效電子能量測(cè)量與刻度。在此基礎(chǔ)上,利用中子飛行時(shí)間測(cè)量原理,在不同時(shí)間窗下將锎(252Cf)源近似分化成若干個(gè)單能中子源,獲得1.9?7.8MeV范圍內(nèi)中子在探測(cè)器中的光輸出響應(yīng)函數(shù)。結(jié)果表明,理論模擬與實(shí)驗(yàn)測(cè)量值在低能段比較吻合,而高能段存在約7.3%的誤差。

    液體閃爍體探測(cè)器,光輸出響應(yīng)函數(shù),能量刻度

    近年來(lái),全球核安全態(tài)勢(shì)日趨復(fù)雜,國(guó)家發(fā)展先進(jìn)、新型核安保技術(shù)以遏制核材料的非法使用。核材料的裂變過(guò)程通常伴隨有一定量的中子和伽馬射線出射,并且出射射線具有一定的時(shí)間關(guān)聯(lián)特性,通過(guò)分析中子伽馬之間的關(guān)聯(lián)特性,可以反推核材料本身某些物理性質(zhì),比如質(zhì)量、同位素豐度等。EJ339A液體閃爍探測(cè)器不僅對(duì)快中子有響應(yīng),而且還能同時(shí)兼顧伽馬射線的探測(cè),在對(duì)混合輻射場(chǎng)射線的探測(cè)方面具有其特有的優(yōu)勢(shì)。探測(cè)器n/γ甄別性能、光響應(yīng)函數(shù)等技術(shù)指標(biāo)決定了該探測(cè)器的探測(cè)精度,一定程度上決定著探測(cè)器的應(yīng)用前景。研究發(fā)現(xiàn)閃爍體的光輸出受到很多因素的影響,如閃爍體的形狀、自吸收特性[1]等。Enqvist等[2]模擬了非摻雜閃爍體EJ309的光輸出響應(yīng)函數(shù)和探測(cè)效率,它比EJ339A具有更高的光輸出和更好的n/γ甄別能力。Naeem等[3]利用Geant4程序模擬了EJ309的光輸出函數(shù),并與Pino的結(jié)果進(jìn)行比較,發(fā)現(xiàn)兩種模擬方法均與實(shí)驗(yàn)結(jié)果符合較好。美國(guó)俄勒岡州立大學(xué)的Liao等[4]對(duì)比了EJ301與EJ339A的脈沖甄別性能,并認(rèn)為含硼液體閃爍體的脈沖形狀甄別效果比普通液體閃爍體差。國(guó)內(nèi)也有研究探討EJ339A的雙脈沖時(shí)間譜分布和中子俘獲時(shí)間測(cè)量的相關(guān)問(wèn)題[5],但未見(jiàn)EJ339A光輸出響應(yīng)函數(shù)的模擬與測(cè)量的相關(guān)報(bào)導(dǎo)。

    本文就EJ339A作為液體閃爍體的基本特性光輸出函數(shù)進(jìn)行了實(shí)驗(yàn)和理論模擬研究。首先通過(guò)探測(cè)器對(duì)不同能量γ源進(jìn)行測(cè)量,從而實(shí)現(xiàn)能量刻度。然后利用中子飛行時(shí)間的方法對(duì)252Cf源進(jìn)行測(cè)量,得到準(zhǔn)單能中子的脈沖幅度譜,結(jié)合γ源能量刻度結(jié)果,得到中子在EJ339A中的光輸出函數(shù),并將實(shí)驗(yàn)測(cè)量數(shù)據(jù)與MCNPX模擬結(jié)果進(jìn)行對(duì)比驗(yàn)證。

    1 中子飛行時(shí)間實(shí)驗(yàn)和模擬方法

    1.1 實(shí)驗(yàn)系統(tǒng)

    本實(shí)驗(yàn)系統(tǒng)主要包括液體閃爍探測(cè)器、信號(hào)處理單元及分析軟件三部分。探測(cè)器使用EJ339A、BC501兩種不同的液體閃爍體,尺寸分別為?12.7cm×7.62 cm、?5.08 cm×5.08 cm,其中,EJ339A中含有5%的硼,其中10B的濃縮度為90%。它的最大發(fā)光波長(zhǎng)為425 nm,光輸出為蒽的65%,碳、氫、氧的原子比例約為3.58:6.21:1[6];信號(hào)處理單元采用意大利CAEN公司生產(chǎn)的DT5751桌面型數(shù)字脈沖波形分析器,其采樣頻率、采樣精度分別為1GS·s?1、10 bit,快脈沖輸入脈沖信號(hào)最低值到最高值的電壓差為2 V[7];分析軟件主要功能包括數(shù)據(jù)采集、數(shù)據(jù)預(yù)處理、數(shù)據(jù)庫(kù)管理、數(shù)據(jù)選擇、數(shù)據(jù)分析、界面管理等模塊,實(shí)現(xiàn)波形數(shù)字化儀的控制、入射粒子信號(hào)的時(shí)幅分析以及n/γ信號(hào)甄別等功能。圖1為本實(shí)驗(yàn)中子飛行時(shí)間測(cè)量的基本布局。圖1中252Cf源的中子產(chǎn)額為106n·s?1,中子源到BC501和EJ339A的距離分別為10 cm和90 cm。

    圖1 中子飛行時(shí)間測(cè)量實(shí)驗(yàn)布局 Fig.1 Detectors and instruments arrangement used in the TOF experiment.

    1.2 實(shí)驗(yàn)方法

    1.2.1 γ射線刻度

    由于液閃探測(cè)器主要由低原子序數(shù)的碳、氫等元素組成,其光電吸收截面很小,因而γ射線與探測(cè)器相互作用后不存在明顯的光電峰,這為此類探測(cè)器的能量刻度帶來(lái)困難。但γ射線入射到液體閃爍探測(cè)器中發(fā)生康普頓散射的概率卻最大,從而形成康普頓邊緣[8],這為能量刻度提供了另一種方法。γ能量刻度數(shù)據(jù)處理中最關(guān)鍵的步驟就是確定康普頓邊緣所在的確切道數(shù),本文通過(guò)康普頓邊緣半高點(diǎn)來(lái)讀出峰位。康普頓電子的最大能量與入射γ射線的能量存在如下關(guān)系[9]:

    式中:02表示電子的靜止能量,即0.511MeV。本次實(shí)驗(yàn)中所用的γ放射源包括22Na、133Ba、137Cs和60Co,對(duì)應(yīng)的γ能量由低到高依次為0.356MeV、0.511MeV、0.661MeV、1.25MeV和1.274MeV。

    1.2.2 中子飛行時(shí)間方法

    中子飛行時(shí)間方法是以BC501探測(cè)到的γ信號(hào)作為事件飛行時(shí)間的起點(diǎn),EJ339A探測(cè)到的中子信號(hào)作為事件飛行時(shí)間的終點(diǎn),獲得252Cf源自發(fā)裂變時(shí)間關(guān)聯(lián)符合實(shí)驗(yàn)測(cè)量譜,其能量由中子的飛行距離和與γ光子的飛行時(shí)間差值確定。一般情況下,被探測(cè)的裂變事件均發(fā)生在50ns內(nèi),因此,依據(jù)實(shí)驗(yàn)的空間布局、探測(cè)器時(shí)間特性、數(shù)字化儀的采樣速率、源能譜特征等將短符合時(shí)間窗控制在100ns內(nèi)。兩種探測(cè)器的信號(hào)是在同一個(gè)數(shù)字化儀中進(jìn)行脈沖信號(hào)處理,探測(cè)到的每一個(gè)脈沖在進(jìn)入數(shù)字化儀后就會(huì)有一個(gè)唯一的時(shí)間戳,兩個(gè)信號(hào)的時(shí)間戳的差值小于100ns時(shí)就會(huì)被當(dāng)作是同一次裂變事件記錄下來(lái),從而由EJ339A探測(cè)到中子飛行時(shí)間分布關(guān)聯(lián)252Cf 源自發(fā)裂變的中子能譜。

    在中子飛行時(shí)間譜中不同的時(shí)間窗對(duì)應(yīng)著不同能量的中子,中子能量由飛行距離和γ峰位確定,由于γ射線速度近似光速,因此以γ射線到達(dá)BC501作為時(shí)間起點(diǎn),以中子到達(dá)EJ339A作為時(shí)間終點(diǎn),以此近似獲得中子飛行時(shí)間。252Cf是自發(fā)裂變中子源,由于不同能量的中子對(duì)應(yīng)著不同的飛行時(shí)間,在90 cm飛行距離內(nèi),14.2 MeV和1 MeV能量的中子飛行時(shí)間分別為17.3 ns和65 ns。低能中子的飛行時(shí)間與能量有如下關(guān)系[10]:

    式中:為中子飛行距離;為中子靜止質(zhì)量;為中子的動(dòng)能。通過(guò)中子飛行時(shí)間即可獲知該時(shí)刻中子能量,極短時(shí)間窗內(nèi)中子脈沖可近似為單能中子脈沖。根據(jù)準(zhǔn)單能中子的脈沖幅度分布圖,找到其最大反沖質(zhì)子光輸出的位置,進(jìn)而獲得中子響應(yīng)函數(shù)。

    1.3 MCNPX模擬方法

    MCNPX是一種被廣泛應(yīng)用的粒子輸運(yùn)計(jì)算模擬軟件[11]。它可以準(zhǔn)確模擬中子在液體閃爍體中發(fā)生的每次碰撞,統(tǒng)計(jì)碰撞過(guò)程中中子能量損失,進(jìn)而給出相對(duì)應(yīng)脈沖幅度譜。本文針對(duì)EJ339A探測(cè)器,在MCNPX中準(zhǔn)確建立其幾何模型及材料組成,模擬真實(shí)情況下單能中子在探測(cè)器中的響應(yīng)。通過(guò)f8卡統(tǒng)計(jì)反沖質(zhì)子能量沉積,并利用MATLAB軟件對(duì)模擬PTRAC文件進(jìn)行后處理,提取入射粒子光輸出分布圖譜,光輸出的模擬結(jié)果以等效電子能量為坐標(biāo)數(shù)值。由于中子在有機(jī)閃爍體中與氫原子發(fā)生碰撞產(chǎn)生的反沖質(zhì)子的能量為:p=cos2n,其中:為反沖質(zhì)子與入射中子之間的夾角,當(dāng)這個(gè)夾角為0時(shí),最大反沖質(zhì)子的能量p與入射中子能量n相同,達(dá)到最大值,反沖質(zhì)子的能量越高,激發(fā)的熒光分子的發(fā)光越強(qiáng),產(chǎn)生的光輸出就高,因此下降沿位置就是最大反沖質(zhì)子的光輸出。通過(guò)模擬的數(shù)據(jù)判斷出圖譜中下降沿位置,即為最大反沖質(zhì)子光輸出位置。因此擬合一系列不同單能中子所對(duì)應(yīng)的最大質(zhì)子能量與光輸出,進(jìn)而獲得質(zhì)子響應(yīng)函數(shù)的數(shù)學(xué)表達(dá)式。

    2 結(jié)果分析與討論

    2.1 γ射線能量刻度

    圖2為不同能量γ源的能譜實(shí)驗(yàn)測(cè)量結(jié)果。由于60Co產(chǎn)生的兩種能量的伽馬射線能量相差較小,因此利用二者的平均能量1.25MeV參與擬合。根據(jù)式(1)計(jì)算得出5種不同能量的γ源在液閃中的光輸出分別為:0.207 MeV、0.341 MeV、0.478 MeV、1.038 MeV、1.061 MeV。利用Origin軟件找到能譜圖中康普頓邊下降沿位置,以下降沿10%和90%位置的中值處作為康普頓邊沿的道址數(shù)。

    圖2 不同能量γ源的能譜測(cè)量結(jié)果 Fig.2 Energy spectrum of different gamma sources.

    式中:為道數(shù)值;e為電子能量,MeV。γ射線能量刻度的實(shí)驗(yàn)結(jié)果表明,其線性相關(guān)系數(shù)為0.998。

    圖3 γ射線能量刻度擬合曲線 Fig.3 Fitting of energy of Compton edge calibration spectrum.

    2.2 中子響應(yīng)函數(shù)測(cè)量

    圖4為實(shí)驗(yàn)測(cè)量獲得的每一次符合事件的中子飛行時(shí)間-計(jì)數(shù)分布圖。記錄下來(lái)的符合事件是以脈沖甄別結(jié)果為前提的,圖4中所示結(jié)果包括4種不同的符合類型n-γ、γ-γ、γ-n和n-n,前者為EJ339A探測(cè)到的粒子,后者則為BC501探測(cè)到的粒子。因此,圖4中位于負(fù)半軸的峰為n-γ符合峰,位于零點(diǎn)位置代表的是γ-γ和n-n符合事件,而位于正半軸的峰就是γ-n符合峰。由于BC501的尺寸小于EJ339A,探測(cè)效率受制于探測(cè)器尺寸,所以γ-n符合峰小于n-γ符合峰。圖5為飛行時(shí)間與等效電子能量的二維分布圖。由圖5也可以看到不同類型符合事件峰位置的差異。

    圖4252Cf飛行時(shí)間脈沖高度譜 Fig.4 Pulse height distribution of TOF for252Cf.

    為研究不同能量中子在探測(cè)器EJ339A內(nèi)產(chǎn)生的光輸出響應(yīng)函數(shù),分別選取中子的飛行時(shí)間和能量沉積位置作圖,考慮到252Cf的實(shí)際能譜以及測(cè)量時(shí)間內(nèi)記錄下的中子數(shù)量有限,我們僅研究了飛行時(shí)間為21?42ns時(shí)間窗內(nèi)的中子,其能量為1.94?7.78 MeV。不同飛行時(shí)間的中子脈沖幅度譜如圖6所示。

    從不同能量單能中子的脈沖高度譜提取出最大能量反沖質(zhì)子所對(duì)應(yīng)的道址位置,此處采取下降沿終點(diǎn)作為該區(qū)間內(nèi)最大能量反沖質(zhì)子產(chǎn)生的能量沉積位置。通過(guò)式(3)轉(zhuǎn)換為等效電子能量,得到不同

    能量單能中子產(chǎn)生的最大反沖質(zhì)子的光輸出,如圖7所示。利用多項(xiàng)式擬合得到光輸出函數(shù)為:

    圖7 EJ339A的光輸出函數(shù) Fig.7 Light output function on different neutron energies of EJ339A.

    2.3 MCNPX模擬

    圖8是利用MCNPX的光輸出模擬結(jié)果。單能中子入射到EJ339A探測(cè)器中,產(chǎn)生的反沖質(zhì)子并不是單能的,不同能量的反沖質(zhì)子將產(chǎn)生不同的光輸出。因此選取最大能量反沖質(zhì)子及對(duì)應(yīng)的光輸出作圖。利用三種方式對(duì)其進(jìn)行擬合[2, 13?14]:

    式中:L(p)是等效電子能量,MeV;、、、、、、均為擬合參數(shù)。圖9列出了這些參數(shù)的擬合值。這三種擬合方式的相關(guān)系數(shù)都大于0.9999。

    圖9 EJ339A的光輸出模擬曲線 Fig.9 Simulated light output function of EJ339A.

    將實(shí)驗(yàn)所得光輸出函數(shù)曲線與模擬結(jié)果進(jìn)行比較,發(fā)現(xiàn)在低能段實(shí)驗(yàn)值與模擬值符合得較好。但是隨著中子能量不斷上升,二者的差異越來(lái)越大。由于EJ339A的n/γ甄別效果不夠理想,某些低能中子被誤判為γ光子,進(jìn)而影響中子的能量分布,造成模擬值與實(shí)驗(yàn)值的差異。此外,由于在實(shí)驗(yàn)時(shí)以1 ns內(nèi)到達(dá)探測(cè)器的中子近似為準(zhǔn)單能中子,但隨著中子能量增大,這種近似與模擬值越不相符。因此,在高能段采用中子飛行時(shí)間的方法進(jìn)行中子光響應(yīng)函數(shù)刻度并不是很理想。

    3 結(jié)語(yǔ)

    1) EJ339A型液體閃爍體探測(cè)器對(duì)電子能量響應(yīng)的線性關(guān)系。通過(guò)對(duì)不同能量γ射線源能譜的實(shí)驗(yàn)測(cè)量,確定其康普頓邊緣位置,實(shí)現(xiàn)了探測(cè)器γ射線能量刻度,擬合結(jié)果線性較好。

    2) EJ339A型探測(cè)器對(duì)質(zhì)子的非線性響應(yīng)。采用蒙特卡羅方法模擬得到的質(zhì)子光響應(yīng)函數(shù)與中子飛行時(shí)間實(shí)驗(yàn)方法得到的函數(shù)在高能端有約7.3%的誤差存在。

    3) 探測(cè)器的電子和質(zhì)子光輸出函數(shù),為全吸收中子能譜測(cè)量的開(kāi)展提供了技術(shù)參考和理論基礎(chǔ)。

    1 Moszynski M, Kapusta M, Mayhugh M,. Absolute light output of scintillators[J]. IEEE Transactions on Nuclear Science, 1997, 44(3): 1052?1061. DOI: 10.1109/23.603803

    2 Enqvist A, Lawrence C C, Wieger B M,. Neutron light output response and resolution functions in EJ-309 liquid scintillation detectors[J]. Nuclear Instruments and Methods in Physics Research, 2013, A715: 79?86. DOI: 10.1016/j.nima.2013.03.032

    3 Naeem S, Clarke S, Pozzi S. Validation of Geant4 and MCNPX-PoliMi simulations of fast neutron detection with the EJ-309 liquid scintillator[J]. Nuclear Instruments and Methods in Physics Research, 2013, A714: 98?104. DOI: 10.1016/j.nima.2013.02.017

    4 Liao C, Yang H. n/γ pulse shape discrimination comparison of EJ301 and EJ339A liquid scintillation detectors[J]. Annals of Nuclear Energy, 2014, 69: 57?61. DOI: 10.1016/j.anucene.2014.01.039

    5 常樂(lè), 劉龍祥, 王宏偉, 等. EJ339A 俘獲門(mén)控中子探測(cè)器雙脈沖特性[J]. 核技術(shù), 2015, 38(5): 050403. DOI: 10.11889/j.0253-3219.2015.hjs.38.050403
    CHANG Le, LIU Longxiang, WANG Hongwei,. Double pulse waveform spectrum of EJ339A capture-gated neutrondetector[J]. Nuclear Techniques, 2015, 38(5): 050403. DOI: 10.11889/j.0253-3219.2015. hjs.38.050403

    6 EJ-339A boron loaded liquid scintillator[OL]. 2015. http://www.eljentechnology.com/index.php/products/loaded-scintillators/82-ej-339a

    7 CAEN - costruzioni apparecchiature elettroniche nucleari s.p.a. DT5751 2/4 channel 10 bit 2/1 GS/s digitizer webpage[OL]. 2015. http://www.caen.it/csite/CaenProd. jsp?parent=14&idmod=632

    8 Krane K S. Introductory nuclear physics[M]. 3rd ed. New York: Wiley & Sons, 1987

    9 Knoll G F. Radiation detection and measurement[M]. John Wiley & Sons, 2010

    10 沈冠仁. 中子飛行時(shí)間方法及其應(yīng)用[M]. 北京: 原子能出版社, 2007
    SHEN Guanren. The time of flight method and applications[M]. Beijing: Atomic Energy Press, 2007

    11 Pelowitz D. MCNPX user’s manual 2.7.0[R]. Los Alamos National Laboratory, 2011

    12 Novotny T, Büermann L, Guldbakke S,. Response of NE213 liquid scintillation detectors to high-energy photons (7 MeV<γ< 20 MeV)[J]. Nuclear Instruments and Methods in Physics Research, 1997, A400(2): 356?366. DOI: 10.1016/S0168-9002(97)00904-2

    13 Katz R, Sharma S, Homayoonfar M. Detection of energetic heavy ions[J]. Nuclear Instruments and Methods, 1972, 100(1): 13?32. DOI: 10.1016/0029-554X(72) 90264-9

    14 Kornilov N, Fabry I, Oberstedt S,. Total characterization of neutron detectors with a252Cf source and a new light output determination[J]. Nuclear Instruments and Methods in Physics Research, 2009, A599(2): 226?233. DOI: 10.1016/j.nima.2008.10.032

    Theoretical simulation and experimental research of the light output function of liquid scintillator EJ339A

    WANG Tingting1WANG Liang2MO Zhaohong2ZHAO Deshan2XIONG Zhonghua2

    1(Science and Technology on Surface Physics and Chemistry Laboratory, China Academy of Engineering Physics, Jiangyou 621907, China)2(Institute of Material Research, China Academy of Engineering Physics, Mianyang 621700, China)

    Background: The light output function is an elementary aspect for neutron spectrometry measured by an organic scintillator to unfold the neutron spectrum from the pulse height distribution. Purpose: This paper aims to investigate the properties of a boron-loaded liquid scintillator EJ339A, using photon sources and neutrons from a252Cf source. Methods: The light output was defined using Compton electron spectra. Pulse shape discrimination (PSD) and time of flight (TOF) were used to distinguish neutron andg-rays. Combining theoretical modeling and experimental research, four different gamma sources are applied to achieve equivalent calibration of electronic energy, namely22Na,133Ba,137Cs, and60Co, respectively. Accordingly, based on TOF method,252Cf, a spontaneous fission neutron source, is treated as several quasi-monoenergetic neutron sources under different time-windows. Results: The electron light output scale was calibrated by the measured detector response. Neutrons deposit energy in the detector and give a light output function with the neutron energy in the range of 1.9?7.8MeV. Conclusion: It shows that experimental data of light output function are in good agreement with simulation result in the low energy section, while discrepancies arise in the high energy section.

    Liquid scintillator, Light output function, Energy calibration

    WANG Tingting, female, born in 1990, graduated from University of Science and Technology of China in 2013, master student, major in radiation protection and environmental protection

    XIONG Zhonghua, E-mail:xiongzhonghua@caep.cn

    TL812+.2

    10.11889/j.0253-3219.2016.hjs.39.050403

    王亭亭,女,1990年出生,2013年畢業(yè)于中國(guó)科學(xué)技術(shù)大學(xué),現(xiàn)為碩士研究生,輻射防護(hù)與環(huán)境保護(hù)專業(yè)

    熊忠華,E-mail: xiongzhonghua@caep.cn

    2015-10-15,

    2016-03-15

    猜你喜歡
    響應(yīng)函數(shù)中子質(zhì)子
    不同探源距下241Am的α譜模擬與解析
    一類具有Beddington-DeAngelis響應(yīng)函數(shù)的階段結(jié)構(gòu)捕食模型的穩(wěn)定性
    質(zhì)子束放療在腫瘤中的研究新進(jìn)展
    3D打印抗中子輻照鋼研究取得新進(jìn)展
    相機(jī)響應(yīng)函數(shù)定標(biāo)的正則化方法
    淺談質(zhì)子守恒
    克服動(dòng)態(tài)問(wèn)題影響的相機(jī)響應(yīng)函數(shù)標(biāo)定
    基于PLC控制的中子束窗更換維護(hù)系統(tǒng)開(kāi)發(fā)與研究
    DORT 程序進(jìn)行RPV 中子注量率計(jì)算的可靠性驗(yàn)證
    “質(zhì)子”號(hào)一箭發(fā)雙星
    太空探索(2014年6期)2014-07-10 13:06:11
    后天国语完整版免费观看| 亚洲av电影在线进入| 中文字幕另类日韩欧美亚洲嫩草| 大香蕉久久成人网| av一本久久久久| 午夜91福利影院| 热re99久久精品国产66热6| 国产一区有黄有色的免费视频| 成人亚洲精品一区在线观看| 色婷婷av一区二区三区视频| 美女福利国产在线| 欧美日韩亚洲高清精品| 午夜福利,免费看| 麻豆av在线久日| 久久久久国产一级毛片高清牌| 99久久99久久久精品蜜桃| 亚洲精品国产区一区二| 人人妻人人爽人人添夜夜欢视频| 日本av免费视频播放| 在线播放国产精品三级| 侵犯人妻中文字幕一二三四区| 搡老岳熟女国产| 国产精品自产拍在线观看55亚洲 | 后天国语完整版免费观看| 国产片内射在线| 在线 av 中文字幕| av欧美777| 亚洲精品中文字幕一二三四区 | 欧美黑人精品巨大| 久久免费观看电影| 丰满人妻熟妇乱又伦精品不卡| 精品国内亚洲2022精品成人 | a级毛片在线看网站| 交换朋友夫妻互换小说| 高清黄色对白视频在线免费看| 精品福利永久在线观看| 欧美亚洲 丝袜 人妻 在线| 美女主播在线视频| 精品一区二区三区四区五区乱码| 日本vs欧美在线观看视频| av网站在线播放免费| 成人国语在线视频| 免费久久久久久久精品成人欧美视频| 免费久久久久久久精品成人欧美视频| 人人妻人人添人人爽欧美一区卜| 女人爽到高潮嗷嗷叫在线视频| 老熟妇乱子伦视频在线观看| 国产男靠女视频免费网站| 大香蕉久久网| 国产精品美女特级片免费视频播放器 | 欧美大码av| 97人妻天天添夜夜摸| 精品一区二区三卡| 亚洲av欧美aⅴ国产| 亚洲综合色网址| 国产成+人综合+亚洲专区| 国产在线一区二区三区精| 性高湖久久久久久久久免费观看| 女人精品久久久久毛片| 欧美精品啪啪一区二区三区| 久久亚洲精品不卡| 十八禁人妻一区二区| 正在播放国产对白刺激| 老司机靠b影院| 高清视频免费观看一区二区| 在线观看免费午夜福利视频| 91九色精品人成在线观看| av线在线观看网站| 香蕉丝袜av| 三级毛片av免费| 久久午夜综合久久蜜桃| 色综合婷婷激情| 一本大道久久a久久精品| 欧美人与性动交α欧美软件| 亚洲精品一卡2卡三卡4卡5卡| 黄色 视频免费看| 久久久久久人人人人人| 亚洲精品在线观看二区| 日韩视频在线欧美| 久久毛片免费看一区二区三区| 午夜久久久在线观看| 亚洲国产中文字幕在线视频| 欧美久久黑人一区二区| 一级毛片女人18水好多| 一本久久精品| 国产精品99久久99久久久不卡| 乱人伦中国视频| 高清视频免费观看一区二区| 老鸭窝网址在线观看| 国产精品偷伦视频观看了| 在线av久久热| 色婷婷久久久亚洲欧美| 18禁国产床啪视频网站| 国产av又大| 精品少妇久久久久久888优播| 国产高清激情床上av| 亚洲黑人精品在线| 亚洲精品久久午夜乱码| 少妇精品久久久久久久| 午夜福利在线观看吧| 亚洲精品国产色婷婷电影| 国产精品亚洲一级av第二区| 久久久水蜜桃国产精品网| 午夜精品久久久久久毛片777| 精品少妇久久久久久888优播| 欧美 日韩 精品 国产| 19禁男女啪啪无遮挡网站| 操出白浆在线播放| 国产老妇伦熟女老妇高清| h视频一区二区三区| 日韩人妻精品一区2区三区| 在线播放国产精品三级| 大型黄色视频在线免费观看| 欧美日韩亚洲综合一区二区三区_| 在线av久久热| 亚洲三区欧美一区| 亚洲av电影在线进入| 久久久久久亚洲精品国产蜜桃av| 一二三四在线观看免费中文在| 宅男免费午夜| 国产一区二区在线观看av| 国产在视频线精品| 欧美午夜高清在线| 69精品国产乱码久久久| 香蕉久久夜色| 国产91精品成人一区二区三区 | 悠悠久久av| 女性生殖器流出的白浆| 色综合欧美亚洲国产小说| 水蜜桃什么品种好| 在线天堂中文资源库| 久久这里只有精品19| 天天躁日日躁夜夜躁夜夜| 亚洲欧美色中文字幕在线| 久久亚洲真实| 成人18禁高潮啪啪吃奶动态图| 中文亚洲av片在线观看爽 | 久久午夜亚洲精品久久| 久久精品亚洲精品国产色婷小说| 亚洲第一青青草原| 日韩大片免费观看网站| 精品国产亚洲在线| 一个人免费在线观看的高清视频| 老司机福利观看| 黄色成人免费大全| 91成年电影在线观看| 亚洲精品乱久久久久久| 欧美日韩黄片免| 性高湖久久久久久久久免费观看| 国产精品 欧美亚洲| 曰老女人黄片| 悠悠久久av| 欧美变态另类bdsm刘玥| 老司机影院毛片| 日韩人妻精品一区2区三区| 亚洲国产欧美日韩在线播放| 亚洲第一青青草原| 亚洲一区二区三区欧美精品| 国产精品 欧美亚洲| 黄网站色视频无遮挡免费观看| 一区二区日韩欧美中文字幕| 亚洲av日韩在线播放| 首页视频小说图片口味搜索| 欧美日韩国产mv在线观看视频| 精品人妻1区二区| 久久中文字幕一级| av福利片在线| 啦啦啦中文免费视频观看日本| 欧美 日韩 精品 国产| 日本精品一区二区三区蜜桃| 母亲3免费完整高清在线观看| 久久久久国产一级毛片高清牌| 99国产精品99久久久久| 国产精品98久久久久久宅男小说| 免费看a级黄色片| 十八禁网站免费在线| 最新的欧美精品一区二区| 操美女的视频在线观看| 黄色毛片三级朝国网站| 亚洲欧美精品综合一区二区三区| 亚洲第一av免费看| 女性生殖器流出的白浆| 超色免费av| 老司机福利观看| 色播在线永久视频| 免费在线观看影片大全网站| 欧美黄色片欧美黄色片| 首页视频小说图片口味搜索| 美女视频免费永久观看网站| 99国产精品一区二区蜜桃av | 麻豆成人av在线观看| 捣出白浆h1v1| 欧美一级毛片孕妇| 三上悠亚av全集在线观看| 少妇裸体淫交视频免费看高清 | av超薄肉色丝袜交足视频| 50天的宝宝边吃奶边哭怎么回事| 免费少妇av软件| 在线观看免费日韩欧美大片| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美精品啪啪一区二区三区| 国产亚洲欧美在线一区二区| 欧美中文综合在线视频| 亚洲国产成人一精品久久久| 69av精品久久久久久 | 天天添夜夜摸| 成人18禁在线播放| 99久久人妻综合| a级毛片在线看网站| 国产免费视频播放在线视频| 久久天堂一区二区三区四区| 成年女人毛片免费观看观看9 | 精品人妻1区二区| 1024视频免费在线观看| 成人永久免费在线观看视频 | 精品国产国语对白av| 人妻 亚洲 视频| 一区二区av电影网| 建设人人有责人人尽责人人享有的| 99国产综合亚洲精品| 欧美大码av| 国产男靠女视频免费网站| 丝袜人妻中文字幕| 免费日韩欧美在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 大码成人一级视频| 美女国产高潮福利片在线看| 亚洲性夜色夜夜综合| 老司机午夜十八禁免费视频| 亚洲自偷自拍图片 自拍| 在线永久观看黄色视频| 国产97色在线日韩免费| 亚洲专区中文字幕在线| 黄片大片在线免费观看| 手机成人av网站| 女性生殖器流出的白浆| 国产精品一区二区免费欧美| 午夜福利乱码中文字幕| 一区二区av电影网| 不卡一级毛片| 久久天躁狠狠躁夜夜2o2o| 91精品国产国语对白视频| 国产在线观看jvid| 精品久久久久久电影网| 老司机靠b影院| 纵有疾风起免费观看全集完整版| 王馨瑶露胸无遮挡在线观看| 亚洲av美国av| 免费av中文字幕在线| 老熟女久久久| 69精品国产乱码久久久| 日韩一区二区三区影片| 夜夜爽天天搞| 亚洲va日本ⅴa欧美va伊人久久| 捣出白浆h1v1| 两个人看的免费小视频| 日韩精品免费视频一区二区三区| 两人在一起打扑克的视频| 中文亚洲av片在线观看爽 | 99热网站在线观看| 国产欧美日韩一区二区三区在线| 青草久久国产| 国产精品1区2区在线观看. | 黄片大片在线免费观看| av免费在线观看网站| 天天添夜夜摸| 美女福利国产在线| 国产一区二区三区综合在线观看| 国产aⅴ精品一区二区三区波| 国产亚洲精品第一综合不卡| 免费在线观看黄色视频的| 少妇 在线观看| 极品教师在线免费播放| 一区二区三区国产精品乱码| 国产在线视频一区二区| 国产在线精品亚洲第一网站| 国产不卡av网站在线观看| 十分钟在线观看高清视频www| 亚洲三区欧美一区| 人成视频在线观看免费观看| 亚洲色图综合在线观看| 亚洲男人天堂网一区| 国产欧美日韩一区二区三| 午夜福利在线免费观看网站| 一区二区av电影网| 亚洲欧美一区二区三区久久| 天天躁夜夜躁狠狠躁躁| 叶爱在线成人免费视频播放| 欧美中文综合在线视频| 午夜两性在线视频| 成年人免费黄色播放视频| 国产麻豆69| 免费日韩欧美在线观看| 国产伦理片在线播放av一区| 亚洲一卡2卡3卡4卡5卡精品中文| 咕卡用的链子| 国产真人三级小视频在线观看| 亚洲国产欧美一区二区综合| 国产精品一区二区在线不卡| avwww免费| 亚洲成a人片在线一区二区| 飞空精品影院首页| 亚洲熟女精品中文字幕| 激情视频va一区二区三区| cao死你这个sao货| 超碰成人久久| 久久中文字幕一级| 亚洲午夜理论影院| 国产欧美亚洲国产| 99精品在免费线老司机午夜| 亚洲久久久国产精品| www.熟女人妻精品国产| 欧美黑人精品巨大| 午夜久久久在线观看| 一个人免费看片子| 18禁裸乳无遮挡动漫免费视频| 999久久久国产精品视频| 中文欧美无线码| 亚洲专区国产一区二区| 性色av乱码一区二区三区2| 少妇猛男粗大的猛烈进出视频| 人妻一区二区av| 丰满人妻熟妇乱又伦精品不卡| 国产又爽黄色视频| 精品欧美一区二区三区在线| 日本欧美视频一区| 亚洲综合色网址| 色在线成人网| 中文字幕色久视频| 亚洲第一av免费看| 亚洲国产欧美在线一区| 国产在线观看jvid| 婷婷丁香在线五月| 国产精品香港三级国产av潘金莲| 亚洲精品国产一区二区精华液| 日韩成人在线观看一区二区三区| 欧美黄色淫秽网站| 国产野战对白在线观看| 精品久久久久久电影网| 久久精品国产a三级三级三级| 亚洲熟女精品中文字幕| 建设人人有责人人尽责人人享有的| 亚洲国产欧美一区二区综合| 亚洲专区国产一区二区| 久久久久国内视频| 亚洲中文字幕日韩| 欧美成狂野欧美在线观看| 日韩一区二区三区影片| 青草久久国产| 老汉色∧v一级毛片| 亚洲人成电影观看| 动漫黄色视频在线观看| 欧美成狂野欧美在线观看| 一边摸一边做爽爽视频免费| 亚洲五月婷婷丁香| 超碰97精品在线观看| 免费人妻精品一区二区三区视频| 激情视频va一区二区三区| 十八禁人妻一区二区| 一级黄色大片毛片| 免费人妻精品一区二区三区视频| 亚洲中文字幕日韩| tocl精华| 每晚都被弄得嗷嗷叫到高潮| 女警被强在线播放| 老熟女久久久| 成人国产av品久久久| 久久性视频一级片| 欧美亚洲日本最大视频资源| 欧美日韩一级在线毛片| 中文亚洲av片在线观看爽 | 久久午夜综合久久蜜桃| 亚洲av美国av| av国产精品久久久久影院| www.999成人在线观看| 久久影院123| 国产精品久久久人人做人人爽| 国产精品 欧美亚洲| 精品人妻在线不人妻| 国产亚洲午夜精品一区二区久久| 欧美午夜高清在线| 91老司机精品| 国产在线免费精品| 成人精品一区二区免费| 叶爱在线成人免费视频播放| 色婷婷久久久亚洲欧美| av超薄肉色丝袜交足视频| 老司机亚洲免费影院| 丝袜美腿诱惑在线| 757午夜福利合集在线观看| 满18在线观看网站| 九色亚洲精品在线播放| 亚洲国产欧美一区二区综合| 中文亚洲av片在线观看爽 | 高清av免费在线| 日本欧美视频一区| 欧美国产精品一级二级三级| 欧美午夜高清在线| 50天的宝宝边吃奶边哭怎么回事| 国产一区二区激情短视频| 成人精品一区二区免费| 久久人妻福利社区极品人妻图片| 欧美激情久久久久久爽电影 | 天堂俺去俺来也www色官网| 一级片免费观看大全| 亚洲色图av天堂| 一级,二级,三级黄色视频| 又紧又爽又黄一区二区| 日本av免费视频播放| 悠悠久久av| 国产欧美日韩综合在线一区二区| 精品亚洲乱码少妇综合久久| 精品欧美一区二区三区在线| 另类亚洲欧美激情| 久久人人爽av亚洲精品天堂| 国产欧美亚洲国产| 91九色精品人成在线观看| 91av网站免费观看| 亚洲午夜精品一区,二区,三区| 中亚洲国语对白在线视频| 亚洲精品成人av观看孕妇| 嫁个100分男人电影在线观看| 国产亚洲一区二区精品| 国产欧美日韩综合在线一区二区| 色在线成人网| √禁漫天堂资源中文www| 国产一卡二卡三卡精品| 亚洲精品久久午夜乱码| 天天添夜夜摸| 一边摸一边抽搐一进一出视频| av超薄肉色丝袜交足视频| 欧美成人午夜精品| 中文字幕人妻熟女乱码| 人人妻,人人澡人人爽秒播| 成人国产av品久久久| 女人精品久久久久毛片| 男女无遮挡免费网站观看| www.熟女人妻精品国产| 精品人妻熟女毛片av久久网站| 在线观看免费高清a一片| 久久中文字幕一级| 午夜激情av网站| 国产成人影院久久av| 超碰97精品在线观看| 操出白浆在线播放| 亚洲国产毛片av蜜桃av| 国产片内射在线| 99国产精品免费福利视频| 亚洲 欧美一区二区三区| 国产精品一区二区在线观看99| 日韩免费高清中文字幕av| 久久久久久免费高清国产稀缺| 日韩视频一区二区在线观看| 欧美精品av麻豆av| 国产深夜福利视频在线观看| 久久久久久久精品吃奶| 国产激情久久老熟女| 久久久久久免费高清国产稀缺| 中文字幕另类日韩欧美亚洲嫩草| 国精品久久久久久国模美| 宅男免费午夜| bbb黄色大片| 亚洲精品中文字幕在线视频| 日韩欧美三级三区| 国产精品久久电影中文字幕 | 在线 av 中文字幕| 一个人免费在线观看的高清视频| 日本av手机在线免费观看| 国产精品国产av在线观看| 首页视频小说图片口味搜索| 三上悠亚av全集在线观看| 国产熟女午夜一区二区三区| 免费看a级黄色片| 亚洲精品国产一区二区精华液| 18禁国产床啪视频网站| 国产1区2区3区精品| 肉色欧美久久久久久久蜜桃| 水蜜桃什么品种好| 国产av又大| 欧美+亚洲+日韩+国产| 国产精品亚洲av一区麻豆| 欧美黑人精品巨大| 亚洲精品久久午夜乱码| 精品国产乱码久久久久久男人| 水蜜桃什么品种好| 国产熟女午夜一区二区三区| 国产野战对白在线观看| www.熟女人妻精品国产| 国产午夜精品久久久久久| 欧美日韩国产mv在线观看视频| 国产免费福利视频在线观看| 一级a爱视频在线免费观看| 欧美成狂野欧美在线观看| 99国产极品粉嫩在线观看| 国产亚洲欧美精品永久| 99国产精品免费福利视频| 久久精品国产综合久久久| 亚洲黑人精品在线| 老司机靠b影院| 成人国产av品久久久| 亚洲综合色网址| 亚洲男人天堂网一区| 中文字幕高清在线视频| 日本a在线网址| 欧美成人午夜精品| 亚洲国产av影院在线观看| avwww免费| xxxhd国产人妻xxx| 久久中文字幕一级| 国产视频一区二区在线看| 无人区码免费观看不卡 | 黄色a级毛片大全视频| 黑丝袜美女国产一区| 下体分泌物呈黄色| 精品熟女少妇八av免费久了| 国产成人系列免费观看| 国产在线精品亚洲第一网站| 丝袜美腿诱惑在线| 女同久久另类99精品国产91| 精品久久久久久电影网| 日韩欧美一区二区三区在线观看 | 夜夜骑夜夜射夜夜干| 国产黄色免费在线视频| 99久久人妻综合| 国产精品欧美亚洲77777| 亚洲精品国产精品久久久不卡| 在线十欧美十亚洲十日本专区| 亚洲人成电影观看| 97在线人人人人妻| av天堂久久9| 老司机影院毛片| 变态另类成人亚洲欧美熟女 | 久久国产精品男人的天堂亚洲| 久久热在线av| 叶爱在线成人免费视频播放| 亚洲中文字幕日韩| 日韩大片免费观看网站| 宅男免费午夜| 午夜福利影视在线免费观看| 免费人妻精品一区二区三区视频| 国产野战对白在线观看| 国产一区二区 视频在线| 久久中文字幕一级| 久久久精品区二区三区| 久久久久视频综合| 捣出白浆h1v1| 欧美激情久久久久久爽电影 | 在线观看人妻少妇| 免费av中文字幕在线| 90打野战视频偷拍视频| 亚洲专区国产一区二区| 日韩一卡2卡3卡4卡2021年| 欧美变态另类bdsm刘玥| 9191精品国产免费久久| 久久精品成人免费网站| 午夜免费成人在线视频| 国产高清国产精品国产三级| 免费高清在线观看日韩| 亚洲熟女毛片儿| 十分钟在线观看高清视频www| 欧美激情 高清一区二区三区| av电影中文网址| 欧美乱妇无乱码| 亚洲欧美色中文字幕在线| 精品国产乱子伦一区二区三区| 久久精品国产亚洲av高清一级| 婷婷成人精品国产| 高清av免费在线| 女性生殖器流出的白浆| 侵犯人妻中文字幕一二三四区| 性少妇av在线| 亚洲成人国产一区在线观看| 老司机午夜福利在线观看视频 | 亚洲第一欧美日韩一区二区三区 | 国产高清国产精品国产三级| 丝袜美腿诱惑在线| 亚洲精品国产区一区二| av片东京热男人的天堂| 99精品欧美一区二区三区四区| 在线十欧美十亚洲十日本专区| 中文字幕人妻熟女乱码| 80岁老熟妇乱子伦牲交| 色播在线永久视频| 18在线观看网站| videos熟女内射| 人人妻人人添人人爽欧美一区卜| 久久99一区二区三区| 在线亚洲精品国产二区图片欧美| 色尼玛亚洲综合影院| 亚洲人成77777在线视频| 一边摸一边做爽爽视频免费| 国产成人欧美在线观看 | 汤姆久久久久久久影院中文字幕| 91成人精品电影| 黄色a级毛片大全视频| 天天操日日干夜夜撸| 大型黄色视频在线免费观看| 亚洲精品国产色婷婷电影| 91av网站免费观看| 欧美日韩亚洲高清精品| 在线播放国产精品三级| 少妇裸体淫交视频免费看高清 | 熟女少妇亚洲综合色aaa.| 久久久国产精品麻豆| 国产日韩欧美视频二区| 婷婷丁香在线五月| 国产免费现黄频在线看| 动漫黄色视频在线观看| 美女扒开内裤让男人捅视频| 日韩视频一区二区在线观看| 亚洲精品国产区一区二| svipshipincom国产片| 香蕉丝袜av| 欧美 日韩 精品 国产| 亚洲精品国产色婷婷电影|