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

    4π平坦效率3He中子探測(cè)器陣列的設(shè)計(jì)與模擬分析

    2020-11-17 08:40:38郝子銳范功濤劉龍祥王宏偉胡新榮李鑫祥王俊文戈松雨
    核技術(shù) 2020年11期
    關(guān)鍵詞:束流中子探測(cè)器

    郝子銳 范功濤 劉龍祥 王宏偉,3 張 岳 胡新榮 李鑫祥王俊文 匡 攀 戈松雨

    1(中國(guó)科學(xué)院上海應(yīng)用物理研究所 上海 201800)

    2(中國(guó)科學(xué)院大學(xué) 北京 100049)

    3(中國(guó)科學(xué)院上海高等研究院 上海 201210)

    4(上海科技大學(xué)物質(zhì)科學(xué)與技術(shù)學(xué)院 上海 201210)

    5(上海潤(rùn)理真空有限公司 上海 201806)

    上海激光電子 γ 源(Shanghai Laser Electron Gamma Source,SLEGS)是在建的國(guó)家重大科技基礎(chǔ)設(shè)施-上海光源二期建設(shè)的16 條線(xiàn)站之一,SLEGS 采用外部CO2激光和上海光源儲(chǔ)存環(huán)的3.5 GeV 電子束以一定的夾角發(fā)生激光康普頓散射(Laser Compton Scattering,LCS)產(chǎn)生γ光,并通過(guò)連續(xù)改變激光和電子束夾角的辦法,提供光通量在105~107s-1,能量在 0.4~20 MeV 范圍內(nèi)的 γ 束[1-2]。這個(gè)能區(qū)覆蓋了核結(jié)構(gòu)研究、核數(shù)據(jù)測(cè)量、核天體物理等熱點(diǎn)研究區(qū)域,可以開(kāi)展光核反應(yīng)截面測(cè)量,核結(jié)構(gòu)相關(guān)的核共振熒光(Nuclear Resonance Fluorescence,NRF)測(cè)量,也為核廢料的γ 嬗變、γ 成像、醫(yī)學(xué)同位素產(chǎn)生等應(yīng)用研究提供了實(shí)驗(yàn)平臺(tái)。目前國(guó)際上主要運(yùn)行的LCS 光源有基于日本SPring-8 儲(chǔ)存環(huán)的 LEPS/LEPS2 線(xiàn)站[3-4]和位于同一園區(qū)的 NewSUBARU B01 線(xiàn)站[5-8],美國(guó)高通量γ 源HIGS裝置[9],以及歐盟正在羅馬尼亞建設(shè)的極端光物理設(shè)施(ELI-NP)上基于電子直線(xiàn)加速器(LINAC)的 LCS γ 源 VEGA 系統(tǒng)(Variable-Energy Gamma Ray System)[10-12]。

    γ 引起的光核反應(yīng)有核共振熒光(γ,γ/γ′)、光致中子反應(yīng)(γ,n)、光致帶電粒子反應(yīng)(γ,p/d/t/α)、光致裂變(γ,fission)等,其中光致中子反應(yīng)(簡(jiǎn)稱(chēng)光中子反應(yīng))是光核反應(yīng)最重要的一類(lèi)反應(yīng),截面在10~30 MeV能區(qū)的寬共振峰(即偶極巨共振GDR峰)反映出原子核集體運(yùn)動(dòng)的性質(zhì),在中子發(fā)射閾值附近的可分辨共振峰(即矮共振PDR 峰)反映出原子核能級(jí)結(jié)構(gòu)和中子-質(zhì)子相對(duì)分布(即中子皮)等信息。恒星演化中重元素主要通過(guò)快中子俘獲過(guò)程(rprocess)和慢中子俘獲過(guò)程(s-process)產(chǎn)生,核天體物理中重元素合成的網(wǎng)絡(luò)計(jì)算對(duì)中子俘獲截面和光中子截面精度提出較高的要求,同時(shí)這些截面數(shù)據(jù)也是核反應(yīng)研究、核技術(shù)應(yīng)用,核能發(fā)展的重要參數(shù)之一[13-15]。中子俘獲截面可以在白光中子源裝置上直接測(cè)量,但是對(duì)于一些短壽命的核素來(lái)說(shuō)中子俘獲反應(yīng)無(wú)法直接測(cè)量,而光中子截面作為中子俘獲反應(yīng)的逆反應(yīng),可以利用細(xì)致平衡原理通過(guò)測(cè)量(γ,n)反應(yīng)獲得γ強(qiáng)度函數(shù)等信息來(lái)反推,例如反應(yīng)道166Er(γ,n)165Er(γ,n)164Er(γ,n)163Er(γ,n)162Er,Γ(165Er)=10.36 h,Γ(163Er)=75 min,其中162Er 和164Er為自然界穩(wěn)定存在但是豐度很低的p 核(p-nuclei),其來(lái)源和產(chǎn)生機(jī)制還存在很多不確定性。4π 平坦效率3He中子探測(cè)器已廣泛應(yīng)用于LCS光源上光中子(γ,n)截面測(cè)量[8,16-21]。4π 平坦效率3He 中子探測(cè)器設(shè)計(jì)為在1 keV~6 MeV 能量范圍內(nèi),探測(cè)器效率較為穩(wěn)定,這一特點(diǎn)使得利用環(huán)比方法(Ring-ratio method)和直接中子多重?cái)?shù)排序法(Direct Neutronmultiplicity Sorting Method)[16]測(cè)量光中子截面成為可能。

    國(guó)際上的4π平坦效率3He中子探測(cè)器除了用于β緩發(fā)中子實(shí)驗(yàn)測(cè)量外,還用于光中子實(shí)驗(yàn)測(cè)量,日本 NewSUBARU B01 γ 源裝置上的早期 4π 平坦效率3He 中子探測(cè)器裝置[8]采用 20 支3He 正比計(jì)數(shù)管(Eurisys Mesures 公司生產(chǎn)的 94NH45 型3He 正比計(jì)數(shù)管),每個(gè)計(jì)數(shù)管直徑25 mm,長(zhǎng)度500 mm,內(nèi)部充有1 MPa的3He氣體,這20支3He管在聚乙烯慢化體(360 mm×360 mm×500 mm)中分成三環(huán)且在每個(gè)環(huán)內(nèi)均勻分布。內(nèi)、中、外環(huán)分別含有4、8、8只計(jì)數(shù)管,分別距離束流中心38 mm、70 mm、100 mm。中間的通孔為束流通道,直徑30 mm。整個(gè)裝置對(duì)單能中子的效率為75%(1 keV)~43%(3 MeV)。ELINP γ 源 裝 置 上 的 平 坦 效 率 中 子 探 測(cè) 器[17]與NewSUBARU 的裝置類(lèi)似,但在探測(cè)器結(jié)構(gòu)上有所區(qū)別。它整體由聚乙烯慢化體與20支94NH45型和11 支來(lái)自美國(guó)LND 公司的充有10 個(gè)大氣壓3He 氣體的正比計(jì)數(shù)管組成,三個(gè)同心圓環(huán)分別放置4、9、18 只計(jì)數(shù)管,分別距離束流管道中心55 mm、130 mm、160 mm,探測(cè)器效率為35%~40%。

    本文利用Geant4 程序模擬設(shè)計(jì)了SLEGS 線(xiàn)站的4π 平坦效率3He 中子探測(cè)器陣列,計(jì)算了平均效率曲線(xiàn),同時(shí)模擬了γ 作用于209Bi 靶后產(chǎn)生的中子能譜、角分布、到達(dá)時(shí)間分布(Arrival Time Distribution)等信息,為裝置建成后光中子反應(yīng)截面的實(shí)驗(yàn)測(cè)量和數(shù)據(jù)分析提供了方法和參考。

    1 環(huán)比方法與直接中子多重?cái)?shù)排序法

    環(huán)比方法(Ring-ratio method)最早是由Berman提出用于確定光中子(γ,n)反應(yīng)中平均中子能量的方法[18],Ring-ratio指的是利用慢化體中具有多同心圓布局結(jié)構(gòu)的中子計(jì)數(shù)管,計(jì)數(shù)管最外環(huán)的計(jì)數(shù)率與最內(nèi)環(huán)計(jì)數(shù)率之比。中子到達(dá)最外環(huán)需要穿過(guò)的慢化材料最厚,故最外環(huán)對(duì)高能中子敏感,相比而言,最內(nèi)環(huán)收集的更多是低能中子,二者計(jì)數(shù)率之比與平均中子能量相關(guān)。光中子截面測(cè)量需要準(zhǔn)確的中子探測(cè)器效率,而通常中子探測(cè)器的效率不能視為常數(shù)的情況下,通過(guò)Ring-ratio方法可以降低中子探測(cè)器效率的影響。首先通過(guò)Ring-ratio 方法得到平均中子能量,一般而言,探測(cè)器效率與中子能量的關(guān)系是已知的,利用該關(guān)系得到探測(cè)器對(duì)于當(dāng)前中子的探測(cè)效率。這是準(zhǔn)確測(cè)量光中子(γ,n)截面的要求[19]。當(dāng)γ能量大于靶原子的雙中子分離能(S2n)時(shí),Ring-ratio 方法只能給出平均中子能量,實(shí)際上(γ,xn)反應(yīng)(x> 1)產(chǎn)生中子的能量是不同的,這樣一來(lái)增加了多光中子截面的測(cè)量誤差。

    直接中子多重?cái)?shù)排序法是由日本甲南大學(xué)(Konan University)Ustunomiya教授[20-21]領(lǐng)導(dǎo)的合作組提出的多光中子截面測(cè)量方法,它改進(jìn)了Ringratio方法,降低了探測(cè)器效率對(duì)截面的影響。式(1)表示單能近似下的光中子截面:

    式中:Nx為在(γ,xn)反應(yīng)的數(shù)目(x=1,2,3,…);Nγ為入射 γ 的數(shù)目;NT為單位面積的靶核數(shù);ε為探測(cè)器效率。當(dāng)使用的靶較厚時(shí),需要進(jìn)行厚靶修正ξ=(1-e-μt)/(μt)。g因子表示γ 束流中能量高于中子激發(fā)閾值的比例,因?yàn)槭饔幸欢▽挾?,?dāng)束流中存在能量小于中子激發(fā)閾值的γ時(shí),這部分γ無(wú)法引發(fā)光中子反應(yīng),所以需要g因子去除束流寬度對(duì)截面測(cè)量的影響。

    當(dāng)γ 能量大于雙中子分離能時(shí),測(cè)量時(shí)存在多中子符合情況,用Nj表示j重中子符合事件數(shù),為了區(qū)別反應(yīng)產(chǎn)生的中子多重?cái)?shù)Nx與探測(cè)器收集的中子多重?cái)?shù)Nj,j的取值記為single(s),double(d),triple(t),…,m。

    假設(shè)γ能量高于雙中子分離能(S2n)且低于三中子分離能(S3n),即同時(shí)存在(γ,n)和(γ,2n)反應(yīng),此時(shí)探測(cè)器最多測(cè)量到雙中子符合事件j=d(double),式(3)可以展開(kāi)成:

    對(duì)單中子事件,Ns可能來(lái)自于(γ,n)反應(yīng),或者(γ,2n)反應(yīng)的一個(gè)中子,且另一個(gè)中子沒(méi)有在探測(cè)器中產(chǎn)生信號(hào)。此時(shí)需要測(cè)量的量不僅有Ns、Nd,還需分別測(cè)量(γ,n)、(γ,2n)反應(yīng)的中子平均能量,如果使用平坦效率探測(cè)器,將探測(cè)器效率視為常數(shù)ε(E1)=ε(E2)=ε,簡(jiǎn)化了實(shí)驗(yàn)步驟,式(4)和式(5)簡(jiǎn)寫(xiě)為:

    這時(shí),只需要測(cè)量單中子事件數(shù)Ns,雙中子符合事件數(shù)Nd,即可計(jì)算(γ,n)、(γ,2n)反應(yīng)的事件數(shù),進(jìn)而利用式(1)計(jì)算光中子截面。

    2 4π 平坦效率3He 中子探測(cè)器結(jié)構(gòu)設(shè)計(jì)及效率模擬

    除了NewSUBARU 和ELI-NP 用于光中子截面測(cè)量的4π平坦效率探測(cè)器以外,還有一些類(lèi)似的裝置,如加拿大粒子與核物理國(guó)家實(shí)驗(yàn)室(TRIUMF)的 BELEN(BEta-deLayEd Neutron detector)[22],美國(guó)橡樹(shù)嶺國(guó)家實(shí)驗(yàn)室(Oak Ridge National Laboratory)的 HRIBF(Holifield Radioactive Ion Beam Facility)平坦效率中子探測(cè)器裝置3Hen[23]和日本理化學(xué)研究所(RIkagaku KENkyusho/Institute of Physical and Chemical Research,RIKEN)的 BRIKEN 裝置[24]等。這些裝置主要用于β 緩發(fā)中子測(cè)量,探測(cè)器的平坦效率區(qū)間大多低于1 MeV,但在平坦區(qū)間內(nèi)效率高達(dá)60%~80%。表1 總結(jié)了幾種相似結(jié)構(gòu)的4π 平坦效率中子探測(cè)器,SLEGS 的探測(cè)器達(dá)不到BRIKEN的探測(cè)器效率,但SLEGS探測(cè)器的平坦效率區(qū)間較高,對(duì)快中子的測(cè)量有明顯優(yōu)勢(shì)。相比ELI-NP的平坦效率中子探測(cè)器,我們采用26 只正比計(jì)數(shù)管,內(nèi)部充有0.2 MPa 的3He 氣體,即可實(shí)現(xiàn)40%的平坦效率。

    本文的模擬使用Geant4 4.10.04版本,由于模擬使用的中子低于20 MeV,為得到更精確結(jié)果,采用QGSP_BERT_HP 物理過(guò)程,并利用 G4NDL4.5 作為該物理過(guò)程的數(shù)據(jù)庫(kù),該物理過(guò)程和數(shù)據(jù)庫(kù)通過(guò)重復(fù)NewSUBARU 和ELI-NP 的設(shè)計(jì)方案和效率曲線(xiàn)獲得驗(yàn)證。

    圖1 是為SLEGS 線(xiàn)站經(jīng)過(guò)Geant4 模擬計(jì)算優(yōu)化設(shè)計(jì)的4π 平坦效率3He 中子探測(cè)器布局,優(yōu)化探測(cè)器在慢化體中的位置布局,合理設(shè)置中子計(jì)數(shù)管的數(shù)量和降低管內(nèi)氣壓,從而滿(mǎn)足探測(cè)器陣列的平坦效率和降低成本的設(shè)計(jì)要求,通過(guò)表1 的比較可以看到,用較少的探測(cè)器數(shù)量和較低的氣壓達(dá)到高探測(cè)器效率,我們的優(yōu)化結(jié)果滿(mǎn)足了線(xiàn)站的設(shè)計(jì)需求。探測(cè)器采用聚乙烯(450 mm×450 mm×550 mm)作為中子慢化體。慢化后的中子進(jìn)入3He正比計(jì)數(shù)管后通過(guò)n+3He→p+3H+764 keV反應(yīng)獲得中子的計(jì)數(shù)。中心通道是γ 束通道,同時(shí)用于放置實(shí)驗(yàn)靶。探測(cè)器共使用了26 只3He 正比計(jì)數(shù)管,每個(gè)管內(nèi)充有0.2 MPa 的3He 氣體。以束流通道為中心形成三組同心圓,內(nèi)、中、外環(huán)分別由長(zhǎng)度均為500 mm 的6、8、12 只3He 正比計(jì)數(shù)管組成,放置在距離軸心65 mm、110 mm、175 mm 的位置,其中內(nèi)環(huán)3He 正比計(jì)數(shù)管的直徑為25 mm,中、外環(huán)的3He 正比計(jì)數(shù)管的直徑均為50 mm,由0.5 mm 厚的不銹鋼管封裝。模擬中使用單能、各向同性、4π發(fā)射的中子源,產(chǎn)生于γ束通道中心位置。探測(cè)器效率的模擬結(jié)果如圖2 所示,空心圓點(diǎn)表示不同中子能量下的探測(cè)器總效率,其余三條線(xiàn)分別表示內(nèi)環(huán)(Ring-1)、中環(huán)(Ring-2)、外環(huán)(Ring-3)的效率。在1 keV~6 MeV中子能量范圍內(nèi),模擬結(jié)果給出探測(cè)器總效率介于34.4%~45.9%之間,達(dá)到了平坦效率的設(shè)計(jì)要求。

    圖1 SLEGS設(shè)計(jì)的4π平坦效率3He中子探測(cè)器陣列布局Fig.1 The layout of 4π flat-efficiency 3He neutron detector array for SLEGS

    圖2 4π平坦效率3He中子探測(cè)器陣列效率曲線(xiàn)Fig.2 Curves of the total detection efficiency and individual ring's contribution in 4π flat-efficiency 3He neutron array detector

    從圖2 的總探測(cè)器效率曲線(xiàn)可以看出,在1 keV~3 MeV區(qū)間,探測(cè)器效率高于ELI-NP的探測(cè)器效率。而 ELI-NP 采用的是 1.0 MPa,25 mm 直徑的3He管31只,這大大增加了陣列探測(cè)器的造價(jià),而SLEGS的3He陣列探測(cè)器采用了20只50 mm和6只25 mm(均為0.2 MPa)的正比計(jì)數(shù)管,通過(guò)Geant4模擬計(jì)算調(diào)整了布局位置和數(shù)量,增加了中子接收面積,同時(shí)增加了中子穿過(guò)正比計(jì)數(shù)管的厚度,彌補(bǔ)了3He 氣體壓強(qiáng)的降低對(duì)探測(cè)器效率減弱的影響,通過(guò)合理布局,使得效率有小幅度提升,達(dá)到與多只高氣壓3He管類(lèi)似的探測(cè)器效率。

    表1 同類(lèi)型的平坦效率中子探測(cè)器陣列比較Table 1 Comparison of similar flat efficiency detectors

    圖3展示了探測(cè)器聚乙烯慢化體及支架部分設(shè)計(jì)方案,中心通道及實(shí)驗(yàn)靶直徑27 mm,25 mm3He正比計(jì)數(shù)管孔道直徑27 mm,50 mm3He正比計(jì)數(shù)管孔道直徑52.4 mm,探測(cè)器及慢化體放在可移動(dòng)支架上,中心高度距離地面1 300 mm,支架可以在±30 mm范圍內(nèi)調(diào)整高度。調(diào)節(jié)好探測(cè)器位置后,固定支架,避免位置移動(dòng)。探測(cè)器前面板處,以束流通道為中心,標(biāo)記以準(zhǔn)直的十字線(xiàn),實(shí)驗(yàn)前用激光校正相對(duì)位置。

    圖3 4π平坦效率3He中子探測(cè)器支架(a)及內(nèi)部結(jié)構(gòu)(b)示意圖Fig.3 Diagram of 4π flat-efficiency 3He neutron detector:the support(a)and internal structure(b)

    3 209Bi(γ,n)反應(yīng)產(chǎn)生的中子分布

    利用Geant4 模擬光中子反應(yīng)和探測(cè)器布局,從而研究出射中子的能譜和角分布,檢驗(yàn)平坦效率3He探測(cè)器測(cè)量中子的能力,首先對(duì)中子能譜、角分布做了模擬。模擬中以209Bi靶為例,靶尺寸為?20 mm×4 mm。

    利用 8~37 MeV 單能 γ 束轟擊209Bi 靶,γ 能量步長(zhǎng)為1 MeV,模擬分析得到30組中子能譜數(shù)據(jù)(圖4所示)。隨著γ能量的增加,出射中子的能量也在增加,峰值對(duì)應(yīng)的中子能量從200 keV 增加至峰值的550 keV,在γ 入射能量超過(guò)雙中子分離能(S2n=14.3 MeV)時(shí),部分中子通過(guò)(γ,2n)反應(yīng)生成,這些中子能量普遍較低,在中子能譜的低能部分形成一個(gè)峰。隨著γ 能量的繼續(xù)增加,雙中子出射的中子產(chǎn)額逐步提高,能譜的峰值對(duì)應(yīng)的中子能量也逐漸增加,重復(fù)上述分布規(guī)律直到三中子分離能。在線(xiàn)性坐標(biāo)系下對(duì)圖4 中的每組中子能譜進(jìn)行分析,計(jì)算峰值對(duì)應(yīng)的能量,得到如圖5(a)所示的出射中子峰值能量分布(圖4顯示為對(duì)數(shù)橫坐標(biāo),是為了能清晰得看到低能中子的分布,其能譜峰值的分布和線(xiàn)性坐標(biāo)下的分布有很大的不同,圖5(a)為線(xiàn)性橫坐標(biāo)下的峰值擬合數(shù)據(jù)),橫坐標(biāo)為入射γ 能量,縱坐標(biāo)為產(chǎn)生的中子峰值能量。圖5(a)中可以看到中子能量有3 個(gè)極值點(diǎn),即在約15 MeV、24 MeV、34 MeV 處,分別對(duì)應(yīng)209Bi 的雙中子分離能(S2n=14.3 MeV)、三中子分離能(S3n=22.4 MeV)和四中子分離能(S4n=29.5 MeV),這主要是隨著入射γ能量的增加,雙中子和三中子反應(yīng)道相繼開(kāi)放導(dǎo)致的。當(dāng)γ能量大于單中子分離能時(shí),僅通過(guò)(γ,n)反應(yīng)產(chǎn)生中子,由于激發(fā)一個(gè)中子需要的能量是固定的,隨入射γ能量增加,傳遞給中子的能量也不斷增加。當(dāng)γ能量略高于雙中子分離能時(shí),中子可以通過(guò)(γ,2n)反應(yīng)產(chǎn)生,此時(shí)γ傳遞給中子的能量最少,中子能量下降至最低。隨著γ 能量繼續(xù)增加,(γ,2n)的中子能量也增加,但當(dāng)γ 能量增加至剛好大于三中子分離能時(shí),中子能量沒(méi)有立即下降,反而在25 MeV 時(shí)才至最低點(diǎn),這種現(xiàn)像在γ 能量大于四中子分離能時(shí)更為明顯。這是因?yàn)椋é茫?n)截面本身較小,在γ 能量略大于三中子分離能時(shí),即使γ 能量足以打開(kāi)該反應(yīng)道,相比(γ,2n)的中子產(chǎn)額,(γ,3n)的中子產(chǎn)額增加不明顯,大于四中子分離能時(shí)亦是如此。對(duì)比圖5(a)中出現(xiàn)的三個(gè)極小值點(diǎn),分別在γ 能量為15 MeV、25 MeV 和35 MeV,對(duì)應(yīng)的中子能量分別為0.11 MeV、0.27 MeV 和0.37 MeV,極小值點(diǎn)的中子能量是不斷增加的,這是因?yàn)棣?能量增加至一定程度后,中子不僅來(lái)源于(γ,xn),有一部分來(lái)源于(γ,np)、(γ,nα)和光致裂變(γ,fission)過(guò)程。當(dāng)γ能量處于8~37 MeV 范圍內(nèi),產(chǎn)生的中子能量大多集中于4 MeV 以下,處于3He 平坦效率探測(cè)器的平坦效率區(qū)間,且在這個(gè)能量范圍內(nèi),探測(cè)器總效率40%~46%更穩(wěn)定。圖5(b)是出射中子的角分布,實(shí)驗(yàn)室系中出射角(θ)定義為靶點(diǎn)處中子動(dòng)量方向與束流正向的夾角。從圖5(b)中可以看出中子是90°對(duì)稱(chēng)出射,只有極少數(shù)中子以0°或180°出射。這樣的分布使得絕大多數(shù)中子都能經(jīng)過(guò)慢化體慢化后進(jìn)入到探測(cè)器中,使探測(cè)器對(duì)中子有4π立體角收集。

    圖4 不同能量的γ產(chǎn)生的中子能譜Fig.4 Energy spectra of neutrons produced by a series of gamma of different energies

    圖5 中子峰值能量與入射γ能量的關(guān)系(209Bi單中子分離能S1n=7.5 MeV,雙中子分離能S2n=14.3 MeV,三中子分離能S3n=22.4 MeV,四中子分離能S4n=29.5 MeV)(a),靶點(diǎn)處出射中子角分布(b)Fig.5 The relationship of neutron peak energy with photon energy(209Bi neutron emission threshold for 1n(Sn=7.46 MeV),2n(S2n=14.3 MeV),3n(S3n=22.4 MeV),and 4n(S4n=29.5 MeV))(a),the number of emitted neutron vs.θ scattering angle in target region(b)

    到達(dá)時(shí)間(Arrival time distribution)是指從中子產(chǎn)生時(shí)刻起至被探測(cè)器收集的時(shí)間,可近似視為中子在慢化體中被慢化的時(shí)間。實(shí)驗(yàn)中γ能量可能高于靶核的多中子分離能,故探測(cè)器陣列可能同時(shí)產(chǎn)生多個(gè)中子符合事件。由于中子是經(jīng)過(guò)慢化以后被3He正比計(jì)數(shù)管記錄到,因此3He正比計(jì)數(shù)管信號(hào)不能直接反映中子的能量和位置信息。此時(shí)需要對(duì)多個(gè)中子的情況分別統(tǒng)計(jì),計(jì)算每個(gè)中子的到達(dá)時(shí)間。由于SLEGS 線(xiàn)站目前提供的γ 能量上限為20 MeV,高于209Bi的雙中子分離能且低于三中子分離能,故最多只能產(chǎn)生雙中子符合事件。例如圖4中入射γ 能量為15 MeV,在雙中子剛出現(xiàn)的時(shí)候,雙中子的能量較低,約為100 keV,而單中子能量較高,約為1 MeV,因此不同能量的中子在聚乙烯慢化體中的慢化時(shí)間不一樣,即不同能量中子的到達(dá)時(shí)間有很大的差異,利用這一時(shí)間差異從數(shù)據(jù)獲取電子學(xué)中來(lái)輔助區(qū)分不同的中子來(lái)源。圖6為單中子與雙中子的到達(dá)時(shí)間分布,在雙中子事件中以到達(dá)探測(cè)器次序區(qū)分兩個(gè)中子,給出每個(gè)中子到達(dá)時(shí)間。通過(guò)對(duì)中子到達(dá)時(shí)間譜的分析可以得到單中子事件數(shù)Ns、雙中子符合事件數(shù)Nd。為了區(qū)分來(lái)自于一次反應(yīng)的雙中子事件和來(lái)自于兩次反應(yīng)的單中子事件,相鄰兩束γ 的時(shí)間間隔不能過(guò)短。同時(shí)為了提高測(cè)量效率,縮短測(cè)量時(shí)間,可以將兩束γ時(shí)間間隔定為400~500 μs。模擬計(jì)算結(jié)果和New SUBARUB B01 的209Bi 實(shí)驗(yàn)測(cè)量[16]有明顯的不同,這種差異來(lái)自于探測(cè)器布局結(jié)構(gòu)的不同,New SUBARUB B01實(shí)驗(yàn)的探測(cè)器三層3He 管分別距離束流管道中心55 mm、130 mm、160 mm,中環(huán)和外環(huán)更近,且內(nèi)環(huán)距離束流中心較近且只有4 只計(jì)數(shù)管,這就使得其探測(cè)器陣列內(nèi)環(huán)對(duì)keV 中子的探測(cè)效率較高,而當(dāng)中子能量大于2 MeV時(shí)探測(cè)器陣列外環(huán)的效率相對(duì)升高。而SELGS的探測(cè)器三環(huán)3He管分別距離束流管道中心65 mm、110 mm、175 mm,內(nèi)、中、外環(huán)分布均勻,內(nèi)環(huán)由于距離靶點(diǎn)較近,采用25 mm小尺寸計(jì)數(shù)管,從而保證每環(huán)的探測(cè)器效率均接近,因此使得中子的慢化時(shí)間均勻,到達(dá)時(shí)間分布較為一致。

    圖6 到達(dá)時(shí)間分布 (a)單中子事件,(b)雙中子事件Fig.6 Arrival time distribution (a)Single neutron event,(b)Double neutron events

    4 結(jié)語(yǔ)

    SLEGS線(xiàn)站用于光中子(γ,n)截面測(cè)量的4π平坦效率3He 探測(cè)器陣列需要滿(mǎn)足低成本、高探測(cè)效率的設(shè)計(jì)需求,通過(guò)Geant4 蒙特卡羅程序模擬,優(yōu)化3He 正比計(jì)數(shù)器在聚乙烯慢化體中的位置分布,減少計(jì)數(shù)管數(shù)量和降低計(jì)數(shù)管內(nèi)的充氣氣壓達(dá)到設(shè)計(jì)要求,同時(shí)通過(guò)模擬光中子反應(yīng)的中子能譜和角分布,到達(dá)時(shí)間等檢驗(yàn)和確認(rèn)了模擬設(shè)計(jì)的合理性和可行性,并獲得了以下的研究結(jié)論;

    1)只使用了20只50 mm和6只25 mm3He正比計(jì)數(shù)管,管內(nèi)充0.2 MPa3He氣體。在降低了制造成本的同時(shí)實(shí)現(xiàn)了與已有光中子3He中子探測(cè)器陣列相近(40%~46%)的平坦效率。

    2)出射中子能譜峰值能量隨中子分離能的變化和中子出射角分布的模擬計(jì)算結(jié)果,對(duì)于SLEGS線(xiàn)站光中子反應(yīng)道出射中子的特點(diǎn)有了更深入的認(rèn)識(shí),研究得到的1~4 MeV 的中子能量分布范圍,90°的對(duì)稱(chēng)角分布特征對(duì)于后續(xù)設(shè)計(jì)飛行時(shí)間中子探測(cè)器[25]具有非常好的參考作用。

    3)3 環(huán)計(jì)數(shù)管位置均勻分布與近距離處采用25 mm 計(jì)數(shù)管的設(shè)計(jì)使得中子的平均慢化時(shí)間接近,探測(cè)器效率接近,保證數(shù)據(jù)獲取的三環(huán)具有一樣的觸發(fā)效率,有利于多中子符合事件的判選和分析。

    4)為4π平坦效率3He探測(cè)器陣列建立了通用的Geant4模擬程序,結(jié)合逆康普頓散射的γ產(chǎn)生程序,就可以完整的模擬和分析光中子截面測(cè)量實(shí)驗(yàn),為未來(lái)的線(xiàn)站用戶(hù)提供標(biāo)準(zhǔn)的光中子實(shí)驗(yàn)?zāi)M工具包。

    平坦效率對(duì)于實(shí)現(xiàn)反應(yīng)單中子和雙中子,三中子事件的區(qū)分和截面的計(jì)算具有非常關(guān)鍵的作用,探測(cè)器陣列對(duì)1 keV~4 MeV的中子具有較高的探測(cè)效率且穩(wěn)定的特點(diǎn),適用于直接中子多重?cái)?shù)排序法測(cè)量,這也是逆康普頓散射γ 源的光中子截面測(cè)量的主要方法。通過(guò)模擬209Bi 靶,尺寸為?20 mm×4 mm,發(fā)現(xiàn)出射中子90°對(duì)稱(chēng)角分布,能譜隨中子分離能規(guī)律變化的特點(diǎn),使得絕大多數(shù)中子能夠進(jìn)入到慢化體中并被探測(cè)器俘獲,且中子能量處于探測(cè)器的平坦效率區(qū)間內(nèi)。

    猜你喜歡
    束流中子探測(cè)器
    第二章 探測(cè)器有反應(yīng)
    EN菌的引力波探測(cè)器
    3D打印抗中子輻照鋼研究取得新進(jìn)展
    高能同步輻射光源低能束流輸運(yùn)線(xiàn)設(shè)計(jì)研究
    中國(guó)散裂中子源加速器注入束流損失調(diào)節(jié)研究
    第二章 探測(cè)器有反應(yīng)
    基于PLC控制的中子束窗更換維護(hù)系統(tǒng)開(kāi)發(fā)與研究
    DORT 程序進(jìn)行RPV 中子注量率計(jì)算的可靠性驗(yàn)證
    ECR中和器束流引出實(shí)驗(yàn)研究
    有7顆彗星已經(jīng)被探測(cè)器造訪(fǎng)過(guò)
    太空探索(2014年9期)2014-07-10 13:06:30
    巨乳人妻的诱惑在线观看| 91av网站免费观看| 亚洲精品久久成人aⅴ小说| 日本免费一区二区三区高清不卡 | 欧美老熟妇乱子伦牲交| 免费av毛片视频| 亚洲,欧美精品.| 国产精品久久久人人做人人爽| 久久狼人影院| 亚洲在线自拍视频| 91精品国产国语对白视频| 精品高清国产在线一区| 国产精品国产高清国产av| 一级黄色大片毛片| 12—13女人毛片做爰片一| 欧美成人午夜精品| 久久欧美精品欧美久久欧美| 露出奶头的视频| 后天国语完整版免费观看| 日韩欧美国产在线观看| 国产99白浆流出| 宅男免费午夜| 如日韩欧美国产精品一区二区三区| av天堂久久9| 一二三四在线观看免费中文在| 巨乳人妻的诱惑在线观看| 免费观看人在逋| 亚洲人成电影免费在线| 欧美日本视频| 精品久久久久久,| 露出奶头的视频| 一进一出抽搐动态| 脱女人内裤的视频| 99国产极品粉嫩在线观看| 成人免费观看视频高清| av片东京热男人的天堂| 亚洲一区二区三区不卡视频| 精品国产一区二区三区四区第35| 亚洲五月天丁香| 法律面前人人平等表现在哪些方面| 天天一区二区日本电影三级 | √禁漫天堂资源中文www| 久久久久国产精品人妻aⅴ院| 在线视频色国产色| 69av精品久久久久久| 婷婷丁香在线五月| 97人妻天天添夜夜摸| 国产精品香港三级国产av潘金莲| 满18在线观看网站| 一卡2卡三卡四卡精品乱码亚洲| 日韩国内少妇激情av| 国产三级在线视频| cao死你这个sao货| 91大片在线观看| 精品少妇一区二区三区视频日本电影| 亚洲第一青青草原| 久久久久九九精品影院| 国产在线精品亚洲第一网站| 每晚都被弄得嗷嗷叫到高潮| 国产黄a三级三级三级人| 久久香蕉激情| netflix在线观看网站| 久久久久亚洲av毛片大全| 一个人免费在线观看的高清视频| 十分钟在线观看高清视频www| 香蕉久久夜色| 精品一品国产午夜福利视频| 黄片大片在线免费观看| 丁香六月欧美| 午夜影院日韩av| 精品国产乱码久久久久久男人| 女人精品久久久久毛片| 老熟妇乱子伦视频在线观看| 黄片小视频在线播放| 国产色视频综合| 黑人欧美特级aaaaaa片| 亚洲 国产 在线| 亚洲av电影在线进入| 日本vs欧美在线观看视频| 国产亚洲精品av在线| 亚洲美女黄片视频| 大型av网站在线播放| 国产激情久久老熟女| av免费在线观看网站| 少妇被粗大的猛进出69影院| 男人操女人黄网站| 丁香六月欧美| 欧美在线黄色| 青草久久国产| 国产亚洲精品综合一区在线观看 | 女人被躁到高潮嗷嗷叫费观| 亚洲狠狠婷婷综合久久图片| 亚洲精品国产一区二区精华液| 国语自产精品视频在线第100页| av视频免费观看在线观看| 麻豆av在线久日| 成人18禁高潮啪啪吃奶动态图| 亚洲成人免费电影在线观看| 国产精品爽爽va在线观看网站 | www.熟女人妻精品国产| 亚洲 欧美 日韩 在线 免费| 亚洲精品美女久久久久99蜜臀| 久久久久九九精品影院| 亚洲国产欧美网| 色在线成人网| 久久久久九九精品影院| 可以在线观看毛片的网站| 婷婷精品国产亚洲av在线| 免费观看人在逋| 亚洲精品国产一区二区精华液| 又紧又爽又黄一区二区| 99热只有精品国产| 欧美日韩福利视频一区二区| 欧美性长视频在线观看| 美女高潮喷水抽搐中文字幕| 性少妇av在线| 亚洲va日本ⅴa欧美va伊人久久| 一区福利在线观看| 午夜精品久久久久久毛片777| 91精品国产国语对白视频| 最好的美女福利视频网| 又大又爽又粗| 日日爽夜夜爽网站| 男女下面进入的视频免费午夜 | 国产麻豆69| 日韩国内少妇激情av| 国产av一区二区精品久久| 一区二区三区精品91| 欧美最黄视频在线播放免费| 日本欧美视频一区| 男人操女人黄网站| 久久久久国产精品人妻aⅴ院| 免费av毛片视频| 亚洲一卡2卡3卡4卡5卡精品中文| 天天一区二区日本电影三级 | 91精品国产国语对白视频| 无遮挡黄片免费观看| 一个人观看的视频www高清免费观看 | 国产精品乱码一区二三区的特点 | 久久婷婷人人爽人人干人人爱 | 日韩人妻高清精品专区| 欧美zozozo另类| 国产综合懂色| 日韩欧美精品v在线| 最近中文字幕高清免费大全6 | 村上凉子中文字幕在线| 国产成人影院久久av| 欧美潮喷喷水| 亚洲国产精品sss在线观看| 男女视频在线观看网站免费| 波多野结衣巨乳人妻| 美女免费视频网站| 99精品久久久久人妻精品| 成人永久免费在线观看视频| 亚洲自偷自拍三级| 成人鲁丝片一二三区免费| 欧美最新免费一区二区三区| 国产精品爽爽va在线观看网站| 人妻久久中文字幕网| 亚洲国产高清在线一区二区三| 国国产精品蜜臀av免费| 在线播放国产精品三级| 最好的美女福利视频网| 在现免费观看毛片| 欧美不卡视频在线免费观看| 亚洲精品国产成人久久av| 99久国产av精品| 欧美区成人在线视频| 国产成人福利小说| 久久精品国产99精品国产亚洲性色| 三级男女做爰猛烈吃奶摸视频| 少妇的逼好多水| 桃色一区二区三区在线观看| 亚洲电影在线观看av| or卡值多少钱| 校园春色视频在线观看| 成人av一区二区三区在线看| 性插视频无遮挡在线免费观看| 精品一区二区三区av网在线观看| eeuss影院久久| 给我免费播放毛片高清在线观看| 国产精品不卡视频一区二区| 国产蜜桃级精品一区二区三区| 成人性生交大片免费视频hd| 变态另类成人亚洲欧美熟女| 我要搜黄色片| 午夜精品一区二区三区免费看| 欧美日本亚洲视频在线播放| 麻豆成人午夜福利视频| .国产精品久久| 成人特级av手机在线观看| 欧美极品一区二区三区四区| 此物有八面人人有两片| 久久午夜亚洲精品久久| 国产亚洲精品综合一区在线观看| 老熟妇仑乱视频hdxx| 日韩中文字幕欧美一区二区| 真实男女啪啪啪动态图| 国产在线精品亚洲第一网站| 免费无遮挡裸体视频| 深夜a级毛片| 国产乱人伦免费视频| 久久精品国产自在天天线| 男人的好看免费观看在线视频| 亚洲va在线va天堂va国产| 啦啦啦观看免费观看视频高清| 国产爱豆传媒在线观看| 久久久久久大精品| 免费搜索国产男女视频| 国产亚洲欧美98| 欧美极品一区二区三区四区| 久久久久免费精品人妻一区二区| 人妻夜夜爽99麻豆av| 久久久久久国产a免费观看| 国产黄a三级三级三级人| 观看美女的网站| 国产精品爽爽va在线观看网站| 日韩亚洲欧美综合| 久久久久九九精品影院| 久久久久精品国产欧美久久久| 亚洲经典国产精华液单| 97超视频在线观看视频| 国产亚洲精品久久久com| 国产精品电影一区二区三区| 麻豆久久精品国产亚洲av| 亚洲内射少妇av| 可以在线观看的亚洲视频| 午夜视频国产福利| 亚洲专区中文字幕在线| 亚洲18禁久久av| 久久天躁狠狠躁夜夜2o2o| 人妻久久中文字幕网| 中文字幕av在线有码专区| ponron亚洲| 亚洲无线在线观看| 桃红色精品国产亚洲av| 中文字幕久久专区| 欧美又色又爽又黄视频| 国产黄色小视频在线观看| 干丝袜人妻中文字幕| 国产黄片美女视频| 97人妻精品一区二区三区麻豆| 亚洲国产欧美人成| 男人和女人高潮做爰伦理| 日韩大尺度精品在线看网址| 网址你懂的国产日韩在线| 村上凉子中文字幕在线| 亚洲无线在线观看| 别揉我奶头~嗯~啊~动态视频| 直男gayav资源| 亚洲国产精品合色在线| 国产精品福利在线免费观看| 舔av片在线| 亚洲avbb在线观看| 精品福利观看| 国内揄拍国产精品人妻在线| 久久久久国内视频| 精品99又大又爽又粗少妇毛片 | 日本成人三级电影网站| 久久精品国产自在天天线| 女人十人毛片免费观看3o分钟| 欧洲精品卡2卡3卡4卡5卡区| h日本视频在线播放| 搡老岳熟女国产| 精品免费久久久久久久清纯| 97热精品久久久久久| 99热6这里只有精品| a级一级毛片免费在线观看| 亚洲一区二区三区色噜噜| 国产精品日韩av在线免费观看| 桃红色精品国产亚洲av| 国产极品精品免费视频能看的| 男女做爰动态图高潮gif福利片| 国产精品自产拍在线观看55亚洲| 亚洲黑人精品在线| 午夜激情欧美在线| 国产精品一区二区免费欧美| 久久精品国产亚洲av涩爱 | 午夜精品在线福利| 精品久久久久久久人妻蜜臀av| 日本免费一区二区三区高清不卡| 老熟妇乱子伦视频在线观看| 啪啪无遮挡十八禁网站| 干丝袜人妻中文字幕| 国产精品99久久久久久久久| 欧美一区二区精品小视频在线| 长腿黑丝高跟| 免费观看的影片在线观看| 尤物成人国产欧美一区二区三区| 精品福利观看| 色综合亚洲欧美另类图片| 夜夜爽天天搞| 欧美高清性xxxxhd video| 成人高潮视频无遮挡免费网站| 国产 一区精品| 日本爱情动作片www.在线观看 | 午夜免费激情av| 欧美色欧美亚洲另类二区| 国产色婷婷99| 白带黄色成豆腐渣| 午夜影院日韩av| 日韩强制内射视频| 男人狂女人下面高潮的视频| 在线播放无遮挡| 国产精品久久久久久亚洲av鲁大| 亚洲精品国产成人久久av| 午夜精品久久久久久毛片777| 免费在线观看日本一区| 日本黄色视频三级网站网址| 日韩欧美一区二区三区在线观看| 久久午夜福利片| 久久久久精品国产欧美久久久| 在线观看免费视频日本深夜| 天堂动漫精品| 三级毛片av免费| 国产精品野战在线观看| 99国产极品粉嫩在线观看| 偷拍熟女少妇极品色| 精品久久久久久成人av| 欧美日韩综合久久久久久 | 成人无遮挡网站| 国产在线精品亚洲第一网站| 男女啪啪激烈高潮av片| 极品教师在线免费播放| 成人二区视频| 免费av不卡在线播放| 嫩草影视91久久| 久久这里只有精品中国| 亚洲熟妇中文字幕五十中出| 深爱激情五月婷婷| 午夜免费激情av| 麻豆av噜噜一区二区三区| 九色国产91popny在线| 天天一区二区日本电影三级| 免费在线观看影片大全网站| 婷婷六月久久综合丁香| 久久人人爽人人爽人人片va| 色精品久久人妻99蜜桃| 欧美中文日本在线观看视频| 在线国产一区二区在线| 国内少妇人妻偷人精品xxx网站| 欧美性猛交黑人性爽| 欧美3d第一页| 欧美激情国产日韩精品一区| 18禁黄网站禁片免费观看直播| 色精品久久人妻99蜜桃| av在线老鸭窝| 精品人妻偷拍中文字幕| 乱人视频在线观看| 热99re8久久精品国产| 欧美日韩瑟瑟在线播放| 在线a可以看的网站| 午夜福利18| 91av网一区二区| 亚洲熟妇中文字幕五十中出| 久久中文看片网| 嫩草影院入口| 欧美国产日韩亚洲一区| av在线亚洲专区| 日韩精品青青久久久久久| 欧美zozozo另类| 国产色婷婷99| 成人av一区二区三区在线看| 精品久久国产蜜桃| 可以在线观看的亚洲视频| 联通29元200g的流量卡| 国内精品美女久久久久久| 日韩欧美三级三区| 亚洲欧美清纯卡通| 天堂网av新在线| 国产精品综合久久久久久久免费| 天堂√8在线中文| 亚洲欧美激情综合另类| 国产精品不卡视频一区二区| 日韩欧美国产一区二区入口| 天堂动漫精品| 99热6这里只有精品| 欧美+日韩+精品| 亚洲美女视频黄频| 美女cb高潮喷水在线观看| 亚洲在线观看片| 在线观看av片永久免费下载| 国产精品女同一区二区软件 | 亚洲va日本ⅴa欧美va伊人久久| 精品人妻1区二区| 国产精品嫩草影院av在线观看 | 我的老师免费观看完整版| 极品教师在线免费播放| 国产精品永久免费网站| 日韩欧美 国产精品| 校园人妻丝袜中文字幕| 中文字幕人妻熟人妻熟丝袜美| 日本免费一区二区三区高清不卡| 一个人免费在线观看电影| 中文字幕人妻熟人妻熟丝袜美| 可以在线观看毛片的网站| 亚洲色图av天堂| 国产精品久久久久久精品电影| 欧美xxxx黑人xx丫x性爽| 午夜爱爱视频在线播放| 赤兔流量卡办理| 欧美+亚洲+日韩+国产| 麻豆国产97在线/欧美| 干丝袜人妻中文字幕| 久久香蕉精品热| 琪琪午夜伦伦电影理论片6080| 午夜激情福利司机影院| 亚洲色图av天堂| 国产精品三级大全| 欧美成人免费av一区二区三区| 亚洲成人免费电影在线观看| 乱人视频在线观看| 久久久久精品国产欧美久久久| 97碰自拍视频| 久久亚洲真实| 99久久精品一区二区三区| 少妇被粗大猛烈的视频| 精品一区二区三区视频在线| 欧美绝顶高潮抽搐喷水| 国产精品国产三级国产av玫瑰| 两人在一起打扑克的视频| 悠悠久久av| 久久久久国内视频| 亚洲成av人片在线播放无| 国产精品久久久久久精品电影| 久99久视频精品免费| 看免费成人av毛片| 亚洲熟妇中文字幕五十中出| 久久草成人影院| 免费在线观看影片大全网站| 久久精品国产清高在天天线| 国模一区二区三区四区视频| 国产大屁股一区二区在线视频| 免费av不卡在线播放| 亚洲国产欧美人成| 国产成人aa在线观看| 亚洲va日本ⅴa欧美va伊人久久| 搡老岳熟女国产| 国产高清有码在线观看视频| 亚洲av不卡在线观看| 日韩,欧美,国产一区二区三区 | av.在线天堂| 综合色av麻豆| 国产av一区在线观看免费| 国产高清视频在线播放一区| videossex国产| 久久久久久久久久成人| 99热只有精品国产| 国产精品久久久久久久久免| 狠狠狠狠99中文字幕| 精品久久久噜噜| 日韩精品有码人妻一区| 午夜激情欧美在线| .国产精品久久| 18禁在线播放成人免费| 精品一区二区三区av网在线观看| eeuss影院久久| 少妇人妻精品综合一区二区 | 国产av在哪里看| 观看美女的网站| 国产aⅴ精品一区二区三区波| 国产女主播在线喷水免费视频网站 | 99九九线精品视频在线观看视频| 亚洲av不卡在线观看| 一区二区三区高清视频在线| 他把我摸到了高潮在线观看| 非洲黑人性xxxx精品又粗又长| 变态另类成人亚洲欧美熟女| 嫩草影院新地址| 午夜日韩欧美国产| 色吧在线观看| 日本黄色片子视频| 变态另类成人亚洲欧美熟女| 欧美一区二区国产精品久久精品| 美女高潮的动态| 日韩中文字幕欧美一区二区| 91久久精品国产一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 欧美一区二区精品小视频在线| 成人午夜高清在线视频| 99久久精品一区二区三区| 村上凉子中文字幕在线| 给我免费播放毛片高清在线观看| 一个人看视频在线观看www免费| 国国产精品蜜臀av免费| 亚洲天堂国产精品一区在线| 无遮挡黄片免费观看| 亚洲欧美日韩高清专用| .国产精品久久| 欧美在线一区亚洲| 两性午夜刺激爽爽歪歪视频在线观看| av在线天堂中文字幕| 亚洲性久久影院| 狂野欧美激情性xxxx在线观看| 免费在线观看日本一区| 香蕉av资源在线| 亚洲久久久久久中文字幕| 国产av在哪里看| 97热精品久久久久久| 成人国产综合亚洲| 成人午夜高清在线视频| 欧美精品啪啪一区二区三区| 搡老妇女老女人老熟妇| 18+在线观看网站| 久久久久久国产a免费观看| 最近最新中文字幕大全电影3| 精华霜和精华液先用哪个| xxxwww97欧美| 久久精品人妻少妇| 能在线免费观看的黄片| 两个人的视频大全免费| 99在线视频只有这里精品首页| 免费观看在线日韩| a在线观看视频网站| 国产熟女欧美一区二区| 春色校园在线视频观看| 全区人妻精品视频| 久久午夜福利片| 三级男女做爰猛烈吃奶摸视频| 一个人观看的视频www高清免费观看| 不卡视频在线观看欧美| 人人妻人人看人人澡| 精品久久久久久久久av| 免费一级毛片在线播放高清视频| 亚洲熟妇中文字幕五十中出| 大又大粗又爽又黄少妇毛片口| 十八禁网站免费在线| 亚洲在线自拍视频| 99在线人妻在线中文字幕| 成人欧美大片| 九九热线精品视视频播放| 亚洲无线观看免费| 欧美激情在线99| 校园人妻丝袜中文字幕| 亚洲av五月六月丁香网| 露出奶头的视频| 熟妇人妻久久中文字幕3abv| 少妇猛男粗大的猛烈进出视频 | 国产探花在线观看一区二区| 少妇丰满av| aaaaa片日本免费| 国产免费一级a男人的天堂| 亚洲最大成人手机在线| 国产高清不卡午夜福利| 91在线精品国自产拍蜜月| 18+在线观看网站| 长腿黑丝高跟| 国产综合懂色| 婷婷六月久久综合丁香| 成人鲁丝片一二三区免费| 大型黄色视频在线免费观看| 人妻少妇偷人精品九色| 国产精品久久久久久久电影| 一个人观看的视频www高清免费观看| 色5月婷婷丁香| 日韩强制内射视频| 日韩一区二区视频免费看| 亚洲av免费在线观看| 日韩欧美在线二视频| 五月伊人婷婷丁香| 日韩欧美在线二视频| 久久中文看片网| 99热精品在线国产| 九色成人免费人妻av| 久久久久久久亚洲中文字幕| 他把我摸到了高潮在线观看| 精品久久久久久成人av| 亚洲精品成人久久久久久| 麻豆av噜噜一区二区三区| 亚洲av中文av极速乱 | 亚洲中文字幕日韩| 亚洲aⅴ乱码一区二区在线播放| 国产伦在线观看视频一区| 露出奶头的视频| 美女被艹到高潮喷水动态| 国产男人的电影天堂91| 国产日本99.免费观看| 日韩一本色道免费dvd| 久9热在线精品视频| 三级国产精品欧美在线观看| 免费观看的影片在线观看| 999久久久精品免费观看国产| 乱系列少妇在线播放| 91午夜精品亚洲一区二区三区 | 在线观看免费视频日本深夜| 精品久久久久久久末码| 国产精品98久久久久久宅男小说| 日日干狠狠操夜夜爽| 欧美日韩瑟瑟在线播放| 精品免费久久久久久久清纯| 国产精品亚洲美女久久久| 久久6这里有精品| 少妇熟女aⅴ在线视频| 我要搜黄色片| 久久精品国产亚洲av涩爱 | 午夜免费男女啪啪视频观看 | 亚洲内射少妇av| 久久久久久久亚洲中文字幕| 国产私拍福利视频在线观看| 乱码一卡2卡4卡精品| 99热只有精品国产| 中文字幕av在线有码专区| а√天堂www在线а√下载| 亚洲国产高清在线一区二区三| 亚洲va在线va天堂va国产| 51国产日韩欧美| 亚洲一区高清亚洲精品| 亚洲av免费在线观看| 欧美人与善性xxx| 午夜福利视频1000在线观看| 午夜影院日韩av| 男女那种视频在线观看| 免费看美女性在线毛片视频| 在线a可以看的网站| 一进一出抽搐动态|