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

    低強(qiáng)度脈沖中子束的數(shù)字式n/γ分辨測(cè)量

    2015-12-02 05:37:03歐陽曉平渠紅光張顯鵬劉金良李海濤
    核技術(shù) 2015年6期
    關(guān)鍵詞:測(cè)量方法

    田 耕 歐陽曉平 渠紅光 張顯鵬,2 劉金良 李海濤

    低強(qiáng)度脈沖中子束的數(shù)字式n/γ分辨測(cè)量

    田 耕1歐陽曉平1渠紅光1張顯鵬1,2劉金良1李海濤1

    1(西北核技術(shù)研究所 西安 710024)
    2(西安交通大學(xué) 核科學(xué)與技術(shù)學(xué)院 西安 710049)

    針對(duì)低強(qiáng)度脈沖中子束測(cè)量,使用高速數(shù)字示波器作為數(shù)據(jù)采集設(shè)備,配合BC501A液體閃爍體探測(cè)器組建了數(shù)字式脈沖形狀甄別(Digital Pulse Shape Discrimination, DPSD)測(cè)量系統(tǒng),實(shí)現(xiàn)了中子的n/γ分辨測(cè)量。系統(tǒng)工作時(shí)采集并存儲(chǔ)探測(cè)器輸出的中子與γ射線的脈沖波形及其記錄時(shí)刻,利用DPSD方法甄別中子實(shí)現(xiàn)了中子脈沖高度譜統(tǒng)計(jì);系統(tǒng)具有連續(xù)記錄和具備時(shí)間戳的采集窗記錄兩種工作方式以適應(yīng)不同的脈沖中子束強(qiáng)度,并通過分析數(shù)據(jù)記錄中脈沖波形的位置或時(shí)間戳,實(shí)現(xiàn)了中子事件的時(shí)間信息統(tǒng)計(jì)。使用該系統(tǒng)在Am-Be中子源上使用采集窗工作模式開展了實(shí)驗(yàn),成功獲得中子脈沖幅度譜、中子時(shí)間譜以及n/γ甄別譜。

    脈沖中子束,中子能譜,時(shí)間譜,n/γ分辨,數(shù)字式脈沖形狀甄別

    在圍繞加速器和托卡馬克裝置等大型物理設(shè)施所開展的各類實(shí)驗(yàn)研究中,脈沖中子束測(cè)量是重要的研究方向,對(duì)于物理過程診斷具有重要意義。中子束流中通常伴隨有γ射線,對(duì)中子測(cè)量形成干擾。在脈沖中子束的強(qiáng)度較高時(shí),各個(gè)粒子在探測(cè)器中形成的電流脈沖大量堆積,無法區(qū)分單個(gè)粒子的脈沖波形。這種情況下通常采用電流型探測(cè)系統(tǒng)進(jìn)行測(cè)量,通過結(jié)構(gòu)設(shè)計(jì)或材料選擇等途徑,使探測(cè)裝置在物理原理上只對(duì)中子靈敏,以此排除γ射線對(duì)測(cè)量的干擾[1]。在脈沖中子束的強(qiáng)度較低、單個(gè)粒子在探測(cè)器中形成的電流波形多數(shù)不堆積時(shí),上述電流型探測(cè)系統(tǒng)輸出信號(hào)的幅度往往過低而不能有效測(cè)量。此時(shí)可以采用具備脈沖形狀甄別(Pulse Shape Discrimination, PSD)能力的探測(cè)器,組建計(jì)數(shù)型PSD系統(tǒng)進(jìn)行n/γ分辨測(cè)量。但傳統(tǒng)模擬式PSD系統(tǒng)的計(jì)數(shù)率較低,在測(cè)量脈沖中子束時(shí)難以獲得足夠數(shù)量的事件。同時(shí),測(cè)量中有時(shí)需要獲得中子束流強(qiáng)度隨時(shí)間變化的信息,在束流強(qiáng)度較低時(shí)可通過測(cè)量每個(gè)粒子的入射時(shí)刻來實(shí)現(xiàn)。模擬式PSD系統(tǒng)實(shí)現(xiàn)該功能需要額外的電路,增加了系統(tǒng)的復(fù)雜程度。目前PSD技術(shù)向數(shù)字式測(cè)量的方向發(fā)展,即采用高速數(shù)字化設(shè)備記錄探測(cè)器輸出的電流脈沖,并利用數(shù)字式PSD (Digital PSD, DPSD)等數(shù)字信號(hào)處理方法進(jìn)行n/γ分辨測(cè)量,以期獲得比模擬系統(tǒng)更高的計(jì)數(shù)率和更優(yōu)的甄別效果,并實(shí)現(xiàn)粒子時(shí)間信息統(tǒng)計(jì)等功能[2?4]。因此,采用數(shù)字式n/γ分辨測(cè)量技術(shù)有望較好實(shí)現(xiàn)低強(qiáng)度脈沖中子束測(cè)量。

    本文利用BC501A液體閃爍體探測(cè)器和高帶寬數(shù)字示波器建立了數(shù)字式低強(qiáng)度脈沖中子束測(cè)量系統(tǒng),實(shí)現(xiàn)了n/γ甄別DPSD算法以及粒子入射時(shí)刻統(tǒng)計(jì)算法。該系統(tǒng)具備連續(xù)記錄和采集窗記錄兩種工作方式,以適應(yīng)不同強(qiáng)度的脈沖中子束以及穩(wěn)態(tài)中子束的測(cè)量。使用該系統(tǒng)在Am-Be中子源上開展實(shí)驗(yàn),采用具備時(shí)間戳的采集窗記錄方式獲取了BC501A探測(cè)器在中子和γ射線入射時(shí)產(chǎn)生的電流脈沖波形;利用離線數(shù)據(jù)處理的方式用DPSD等算法得到n/γ甄別譜和中子脈沖幅度譜,并通過處理時(shí)間戳得到中子的時(shí)間譜。

    1 測(cè)量原理與實(shí)驗(yàn)裝置

    1.1 測(cè)量原理

    某些探測(cè)器在不同種類粒子入射時(shí)輸出脈沖的形狀不同,PSD技術(shù)即利用該性質(zhì)實(shí)現(xiàn)粒子分辨測(cè)量。對(duì)于n/γ分辨,目前液體閃爍體探測(cè)器效果較好,常用閃爍體牌號(hào)有BC501A、EJ301和較早的NE213等,其性能相近。這些液體閃爍體在粒子入射時(shí)所產(chǎn)生的熒光脈沖包含持續(xù)時(shí)間不同的快、慢成分,其中慢成分的衰減時(shí)間常數(shù)與粒子種類有關(guān)。反映在探測(cè)器輸出的電流波形上,即脈沖下降沿的速度與入射粒子的種類相關(guān)。具體到n/γ分辨測(cè)量,表現(xiàn)為中子產(chǎn)生的脈沖信號(hào)的下降沿比γ射線產(chǎn)生的要慢。在DPSD測(cè)量中,利用高速模擬數(shù)字轉(zhuǎn)換設(shè)備采集探測(cè)器輸出的電流波形,并用數(shù)字信號(hào)處理的方法提取脈沖波形的形狀信息即可實(shí)現(xiàn)粒子分辨測(cè)量。

    對(duì)于低強(qiáng)度脈沖中子束的測(cè)量,按照具體的束流強(qiáng)度,數(shù)據(jù)采集設(shè)備的工作方式可分為連續(xù)記錄和采集窗記錄兩種。連續(xù)記錄方式即在觸發(fā)信號(hào)的控制下開始連續(xù)采集數(shù)據(jù),直至存儲(chǔ)器記滿或在指令信號(hào)控制下停止采集。此方法適合測(cè)量強(qiáng)度相對(duì)較高的中子束,這種情況下有效的波形數(shù)據(jù)在整體測(cè)量數(shù)據(jù)中占較大成分、基線數(shù)據(jù)相對(duì)較少,例如歐洲聯(lián)合環(huán)(Joint European Torus, JET)和國內(nèi)的HL-2A等裝置的放電過程中D-D及D-T聚變中子測(cè)量[2?3]。而對(duì)于強(qiáng)度較低的中子束,有效事件占總時(shí)間的比例較低,用連續(xù)記錄方法測(cè)量到的數(shù)據(jù)中會(huì)包含大量無用的基線數(shù)據(jù),浪費(fèi)了存儲(chǔ)空間,使得能夠獲取的有效事件數(shù)目減小。在這種情況下,更為有效的采集方法是采集窗記錄方式。該方法在有符合記錄條件的信號(hào)到來時(shí)進(jìn)行采集并存儲(chǔ)一定長(zhǎng)度的數(shù)據(jù),而在無有效信號(hào)的時(shí)間段中不進(jìn)行采集記錄。用這種方法能夠有效實(shí)現(xiàn)較低強(qiáng)度中子束的測(cè)量。且該方法可以通過給每一段記錄數(shù)據(jù)標(biāo)記時(shí)間戳的方法記錄每一個(gè)事件的采集時(shí)刻信息,在測(cè)量中可利用該時(shí)間戳獲取束流強(qiáng)度隨時(shí)間變化的情況。例如Riva等[4]利用該方法在JET上開展粒子束測(cè)量。以上方法均使用數(shù)字化儀作為數(shù)據(jù)采集設(shè)備,其特點(diǎn)是采樣率較低(數(shù)百M(fèi)S·s?1)而幅度分辨率(也稱垂直分辨率,以下簡(jiǎn)稱分辨率)較高(12bits及以上)。

    1.2 實(shí)驗(yàn)裝置

    實(shí)驗(yàn)裝置由Am-Be中子源和測(cè)量系統(tǒng)組成,其中測(cè)量系統(tǒng)主要包括液體閃爍體探測(cè)器、高速數(shù)字示波器、實(shí)驗(yàn)控制PC計(jì)算機(jī)以及其它輔助設(shè)備等,如圖1所示。探測(cè)器使用的閃爍體為BC501A,封裝的尺寸為?50.8 mm×50.8 mm;光電倍增管為濱松R329-02,與閃爍體一起封裝在鋁外殼中,工作時(shí)加高壓?1700 V,在探測(cè)Am-Be源時(shí)陽極輸出脈沖信號(hào)的上升沿tr≈5 ns,底寬tw≈50 ns。

    圖1 實(shí)驗(yàn)裝置結(jié)構(gòu)Fig.1 Structure of experimental setup.

    數(shù)據(jù)采集設(shè)備的指標(biāo)對(duì)測(cè)量有直接影響,最重要的指標(biāo)是采樣率、分辨率和存儲(chǔ)長(zhǎng)度。S?derstr?m等[5]指出,影響DPSD算法的n/γ甄別效果的參數(shù)是分辨率和采樣率。存儲(chǔ)長(zhǎng)度直接決定能夠測(cè)量的脈沖中子束的持續(xù)時(shí)間,無論是在連續(xù)記錄方式或是采集窗記錄方式下,增加存儲(chǔ)深度都將延長(zhǎng)測(cè)量時(shí)間。目前常用于數(shù)字式粒子測(cè)量的數(shù)據(jù)采集設(shè)備主要為數(shù)字化儀和數(shù)字示波器,前者具有量化位數(shù)高、記錄長(zhǎng)度長(zhǎng)的優(yōu)點(diǎn),但通常觸發(fā)控制、基線調(diào)節(jié)、增益設(shè)置等能力較弱。而這些能力對(duì)于實(shí)現(xiàn)某些測(cè)量功能較為重要。數(shù)字示波器通常設(shè)置靈活,在這些方面具備優(yōu)勢(shì),但多數(shù)數(shù)字示波器的分辨率較低,限制了DPSD算法的效果,同時(shí)大存儲(chǔ)深度的示波器價(jià)格較高。因此在脈沖中子束的測(cè)量中,為獲得最優(yōu)的測(cè)量效果,理想的數(shù)據(jù)采集設(shè)備應(yīng)兼具數(shù)字示波器的設(shè)置靈活,以及數(shù)字化儀的量化精度高、存儲(chǔ)長(zhǎng)度大的特點(diǎn)。

    本文選用數(shù)字示波器作為數(shù)據(jù)采集設(shè)備是因?yàn)槠湓诓蓸勇省⒂|發(fā)功能、偏置設(shè)置和增益調(diào)節(jié)等方面具有優(yōu)勢(shì),適用于n/γ分辨測(cè)量和DPSD算法研究[6?8]。具體型號(hào)選用美國泰克公司Tektronix的DPO7104,因其具備較強(qiáng)的采集窗功能,能夠以自觸發(fā)的方式進(jìn)行采集并生成精確到1 ps的時(shí)間戳;該示波器的帶寬(?3 dB)為1 GHz,分辨率為8 bits,實(shí)驗(yàn)中采樣率為10 GS·s?1。裝置中BC501A探測(cè)器輸出脈沖信號(hào)的上升沿tr≈5 ns,可以估算出其頻譜主要處在70 MHz以下,故實(shí)際采樣率遠(yuǎn)高于量化該脈沖信號(hào)所需的最低采樣率,因此數(shù)據(jù)采集設(shè)備指標(biāo)對(duì)n/γ甄別效果的影響主要來源于分辨率。實(shí)驗(yàn)中采用運(yùn)行Window XP的PC計(jì)算機(jī)通過局域網(wǎng)(Local Area Network, LAN)控制示波器,在程序控制下完成數(shù)據(jù)自動(dòng)采集。該實(shí)驗(yàn)控制程序利用Visual C++和虛擬儀器軟件架構(gòu)(Virtual Instrument Software Architecture, VISA)開發(fā),可兼容多種型號(hào)示波器[9]。數(shù)據(jù)處理采用離線方式進(jìn)行,這是由于數(shù)據(jù)采集設(shè)備的采樣率較高、數(shù)據(jù)的吞吐率較大,實(shí)時(shí)數(shù)據(jù)處理較困難。利用Visual BASIC.net開發(fā)接口程序?qū)崿F(xiàn)數(shù)據(jù)讀取、格式轉(zhuǎn)換、計(jì)算結(jié)果存盤和部分圖形繪制等功能,使用MATLAB實(shí)現(xiàn)DPSD算法,在數(shù)據(jù)處理工作中由接口程序調(diào)用MATLAB計(jì)算引擎完成甄別數(shù)據(jù)計(jì)算。

    測(cè)量系統(tǒng)經(jīng)測(cè)試在采集窗方式下最大計(jì)數(shù)率可達(dá)3×105s?1,該指標(biāo)主要受到記錄設(shè)備死時(shí)間的限制。在連續(xù)記錄模式下能夠達(dá)到更高的計(jì)數(shù)率,具體指標(biāo)與探測(cè)器的響應(yīng)速度有關(guān),探測(cè)器在粒子入射下輸出脈沖的寬度越小計(jì)數(shù)率越高。利用LaBr3:Ce無機(jī)閃爍體探測(cè)器在穩(wěn)態(tài)γ射線源上對(duì)系統(tǒng)的數(shù)據(jù)采集部分進(jìn)行考核,得到當(dāng)輸入信號(hào)脈沖底寬為100 ns、堆疊事件率小于20%時(shí),計(jì)數(shù)率可達(dá)到2.2×106s?1。測(cè)量系統(tǒng)可滿足上述探測(cè)器事件率范圍內(nèi)的低強(qiáng)度脈沖中子束的測(cè)量。

    2 數(shù)據(jù)處理算法

    2.1 DPSD算法

    目前國內(nèi)外的研究人員[5?8,10?13]研究了各種DPSD算法,其效果各異。DPSD算法的甄別效果與數(shù)據(jù)采集設(shè)備的采樣率和分辨率有關(guān)[5],而本文組建的數(shù)據(jù)采集系統(tǒng)的特點(diǎn)是分辨率低但采樣率較高,在這種情況下,不同DPSD算法的n/γ甄別效果差異較明顯。本文實(shí)現(xiàn)了三種常用的DPSD算法。

    2.1.1 電荷比較法(Charge Comparison Method, CC)

    閃爍體探測(cè)器在粒子入射時(shí)產(chǎn)生的脈沖信號(hào)的積分表征了該脈沖包含的電荷量。γ射線和中子在液體閃爍體探測(cè)器中形成的脈沖,其下降沿的速度不同,則在給定的兩個(gè)不同的時(shí)間區(qū)間上分別對(duì)脈沖波形積分,積分值的比與粒子種類相關(guān)。設(shè)脈沖信號(hào)經(jīng)數(shù)字化后得到長(zhǎng)度為N的離散序列w(n),甄別因子QCC可按式(1)定義:

    式中,n1、n2、n3、n4為采樣點(diǎn)序號(hào),其中n1,n3∈[1,N),n2,n4∈(1,N],且[n1,n2]≠[n3,n4];kCC為比例系數(shù)。本文利用脈沖下降沿積分與總積分之比作為甄別參數(shù)。

    2.1.2 脈沖梯度分析法(Pulse Gradient Analysis Method, PGA)

    由于γ射線和中子產(chǎn)生的脈沖信號(hào)的下降速度不同,故在脈沖上給定兩個(gè)時(shí)間點(diǎn),所對(duì)應(yīng)波形上的兩點(diǎn)可確定一條直線,該直線的斜率與粒子種類相關(guān)。設(shè)脈沖信號(hào)經(jīng)數(shù)字化后得到長(zhǎng)度為N的離散序列w(n),甄別因子QPGA可按式(2)定義:

    式中,n1、n2為采樣點(diǎn)序號(hào),其中n1,n2∈[1,N],且n1≠n2;kPGA為比例系數(shù)。本文以距脈沖波形峰值某一固定距離的波形與峰值幅度的比作為甄別參數(shù)。

    2.1.3 積分上升時(shí)間法(Integrated Rise Time Method, IRT)

    將脈沖按時(shí)間積分,獲得的積分曲線的上升時(shí)間可反映脈沖下降沿的速度,即積分曲線上升時(shí)間與粒子種類相關(guān)。設(shè)脈沖信號(hào)經(jīng)數(shù)字化后得到長(zhǎng)度為N的離散序列w(n)。按式(3)定義wI(n):

    式中,n∈[1,N]。記wI(n)的最大值為Amax,則甄別因子QIRT可按式(4)定義:

    式中,kIRT、k1、k2為比例系數(shù),其中k1,k2∈(0,1),且k1≠k2。本文選取k1=0.10,k2=0.75。

    2.2 時(shí)間信息統(tǒng)計(jì)

    在連續(xù)記錄方式下,通過確定每個(gè)脈沖事件在整個(gè)數(shù)據(jù)記錄長(zhǎng)度中的位置即可完成事件的時(shí)間信息統(tǒng)計(jì)。在具備時(shí)間戳的采集窗記錄方式下,該功能需要通過分析每個(gè)事件的時(shí)間戳來完成。本文采用的DPO7104數(shù)字示波器所產(chǎn)生的時(shí)間戳不但包含年、月、日信息,且具有1 ps的精度,無法直接利用軟件開發(fā)工具提供的函數(shù)庫進(jìn)行處理。因此在事件時(shí)間信息統(tǒng)計(jì)程序中,采用面向?qū)ο蟮腃++程序設(shè)計(jì)技術(shù)實(shí)現(xiàn)了此類時(shí)間戳的判讀、比較、時(shí)間間隔計(jì)算及統(tǒng)計(jì)等功能[9]。

    3 實(shí)驗(yàn)及數(shù)據(jù)處理

    在Am-Be中子源上開展測(cè)量,實(shí)驗(yàn)裝置在PC計(jì)算機(jī)控制下以自動(dòng)測(cè)量的方式獲取了50組數(shù)據(jù),每組數(shù)據(jù)記錄5000個(gè)事件。使用上述三種DPSD算法對(duì)數(shù)據(jù)進(jìn)行處理,其結(jié)果顯示,IRT方法的甄別效果較好,獲得了n/γ甄別譜;并且可根據(jù)甄別參數(shù)與脈沖幅度的分布情況對(duì)甄別參數(shù)進(jìn)行修正,以提高甄別的效果。而CC方法和PGA方法的甄別效果均不如IRT方法,無法有效區(qū)分中子和γ射線。通過對(duì)判讀為中子的事件利用直接比較法進(jìn)行脈沖幅度統(tǒng)計(jì)可以得到中子脈沖幅度譜[14];對(duì)其時(shí)間戳進(jìn)行統(tǒng)計(jì)可得中子的時(shí)間譜。各組數(shù)據(jù)處理結(jié)果的情況一致。

    以其中的某一組數(shù)據(jù)的處理為例,分別利用CC、PGA和IRT方法處理每個(gè)事件的波形數(shù)據(jù),得到各個(gè)粒子的甄別參數(shù),并按照甄別參數(shù)-計(jì)數(shù)繪制甄別譜。n/γ分辨測(cè)量的效果可以用優(yōu)質(zhì)因子(Figure of Merit, FOM)來衡量。設(shè)甄別譜中γ射線形成的峰的半高寬為Wγ,中子形成的峰的半高寬為Wn,兩個(gè)峰的距離為ΔD,則FOM由式(5)定義[15]:

    利用雙高斯函數(shù)擬合甄別譜,由式(5)得到FOM。各算法的FOM及運(yùn)算耗時(shí)如表1所示。

    表1 CC、PGA和IRT算法的FOM與運(yùn)算耗時(shí)Table 1 FOMs and time consumptions of CC, PGA and IRT.

    按甄別參數(shù)與脈沖幅度進(jìn)行統(tǒng)計(jì)得到DPSD散點(diǎn)圖、按甄別參數(shù)與計(jì)數(shù)統(tǒng)計(jì)得到n/γ甄別譜,如圖2?4所示。

    根據(jù)表1中不同算法的耗時(shí)數(shù)據(jù)可知,三種算法的計(jì)算量相當(dāng)。從FOM值可見IRT方法的甄別效果優(yōu)于CC方法和PGA方法。在圖2?4給出的散點(diǎn)圖和甄別譜中,IRT方法得到的散點(diǎn)圖分成較清晰的兩部分,在甄別譜曲線上表現(xiàn)為中子峰與γ峰基本沒有重疊;而CC方法和PGA方法得到的散點(diǎn)圖中表示中子與γ射線的散點(diǎn)不能清晰分開,在甄別譜曲線上表現(xiàn)為中子峰與γ峰有部分重疊。

    圖2 CC方法的n/γ甄別結(jié)果(a) QCC—脈沖高度散點(diǎn)圖,(b) n/γ甄別譜Fig.2 n/γ discrimination with CC method. (a) QCC vs. pulse height, (b) Counts vs. QCC

    圖3 PGA方法的n/γ甄別結(jié)果(a) QPGA—脈沖高度散點(diǎn)圖,(b) n/γ甄別譜Fig.3 n/γ discrimination with PGA method. (a) QPGA vs. pulse height, (b) Counts vs. QPGA

    圖4 IRT方法的n/γ甄別結(jié)果(a) QIRT—脈沖高度散點(diǎn)圖,(b) n/γ甄別譜Fig.4 n/γ discrimination with IRT method. (a) QIRT vs. pulse height, (b) Counts vs. QIRT

    按照甄別效果最佳的IRT方法的數(shù)據(jù)處理結(jié)果對(duì)粒子種類進(jìn)行判讀。選取甄別譜的中子峰與γ峰之間計(jì)數(shù)率最小值處的甄別參數(shù)值作為界限區(qū)分中子和γ射線,得到該組數(shù)據(jù)的中子計(jì)數(shù)共2330個(gè)。將判讀為中子的事件利用直接比較法進(jìn)行脈沖高度統(tǒng)計(jì),得到中子脈沖高度譜。根據(jù)時(shí)間戳數(shù)據(jù)計(jì)算得到此次數(shù)據(jù)采集的時(shí)間間隔約為45.92 s。將此時(shí)間間隔等分為50個(gè)時(shí)間段,按每個(gè)時(shí)間段內(nèi)的中子事件計(jì)數(shù)統(tǒng)計(jì)得到中子時(shí)間譜。由于所測(cè)的源在測(cè)量時(shí)間內(nèi)可認(rèn)為是穩(wěn)態(tài)的,因此中子時(shí)間譜以一均值為中心呈隨機(jī)漲落,其統(tǒng)計(jì)量為均值Na=46.6,標(biāo)準(zhǔn)差σ=6.9,符合輻射測(cè)量計(jì)數(shù)的1/aN規(guī)律。中子脈沖高度譜和時(shí)間譜如圖5所示。當(dāng)束流強(qiáng)度在測(cè)量時(shí)間內(nèi)發(fā)生變化時(shí),利用該方法可測(cè)得強(qiáng)度隨時(shí)間變化的情況。

    圖5 中子脈沖高度譜(a)和時(shí)間譜(b)Fig.5 Pulse height spectrum (a) and time spectrum (b) of neutrons.

    4 結(jié)語

    采用BC501A液體閃爍體探測(cè)器、高速數(shù)字示波器、PC計(jì)算機(jī)以及基于VISA技術(shù)開發(fā)的實(shí)驗(yàn)控制程序,組建了適用于低強(qiáng)度脈沖中子束測(cè)量的數(shù)字式n/γ分辨測(cè)量系統(tǒng)。該系統(tǒng)具備連續(xù)記錄和具有時(shí)間戳的采集窗記錄等兩種工作方式,能夠適用于不同事件率的測(cè)量場(chǎng)合,最大計(jì)數(shù)率達(dá)到2.2×106s?1。研究中開發(fā)了CC、PGA和IRT三種DPSD算法以及時(shí)間信息統(tǒng)計(jì)算法,實(shí)現(xiàn)了對(duì)脈沖中子束的n/γ分辨測(cè)量和中子時(shí)間信息測(cè)量。在Am-Be中子源上開展了驗(yàn)證實(shí)驗(yàn),獲得了甄別散點(diǎn)圖、n/γ甄別譜、中子脈沖幅度譜以及中子時(shí)間譜。實(shí)驗(yàn)結(jié)果表明,當(dāng)數(shù)據(jù)采集設(shè)備的分辨率較低而采樣率較高時(shí),IRT方法的n/γ甄別效果較好,獲得甄別譜的優(yōu)質(zhì)因子FOM=1.38。

    基于FPGA的硬件并行算法等新興數(shù)字信號(hào)處理技術(shù)在提高數(shù)據(jù)處理速度方面具有良好前景。通過在低強(qiáng)度脈沖中子束數(shù)字式n/γ分辨測(cè)量研究中開發(fā)此類技術(shù),并結(jié)合使用具有高分辨率、高采樣率和大存儲(chǔ)深度的數(shù)據(jù)采集設(shè)備,將可進(jìn)一步擴(kuò)展測(cè)量功能并提高測(cè)量的質(zhì)量和效率,而且有望實(shí)現(xiàn)測(cè)量數(shù)據(jù)的在線實(shí)時(shí)處理。

    1 劉慶兆. 脈沖輻射場(chǎng)診斷技術(shù)[M]. 北京: 科學(xué)出版社, 1994: 492?517

    LIU Qingzhao. Pulsed radiation field diagnosis technology[M]. Beijing: Science Press, 1994: 492?517

    2 Esposito B, Bertalot L, Marocco D, et al. Neutron measurements on Joint European Torus using NE213 scintillator with digital pulse shape discrimination[J]. Review of Scientific Instruments, 2004, 75(10): 3550?3552

    3 Yuan X, Zhang X, Xie X F, et al. First neutron spectrometry measurement at the HL-2A Tokamak[J]. Chinese Physics C, 2013, 37(12): 126001

    4 Riva M, Esposito B, Marocco D. FPGA-based pulse-oriented digital acquisition system for nuclear detectors[C]. IEEE 22ndSymposium on Fusion Engineering, Albuquerque, NM, USA, 2007

    5 S?derstr?m P A, Nyberg J, Wolters R. Digital pulse-shape discrimination of fast neutrons and γ rays[J]. Nuclear Instruments and Methods in Physics Research, 2008, A594: 79?89

    6 Yousefi S, Lucchese L, Aspinall M D. Digital discrimination of neutrons and gamma-rays in liquid scintillators using wavelets[J]. Nuclear Instruments and Methods in Physics Research, 2009, A598: 551?555

    7 Liu G, Aspinall M D, Ma X, et al. An investigation of the digital discrimination of neutrons and γ rays with organic scintillation detectors using an artificial neural network[J]. Nuclear Instruments and Methods in Physics Research, 2009, A607: 620?628

    8 Liu G F, Joyce M J, Ma X D, et al. A digital method for the discrimination of neutrons and γ rays with organic scintillation detectors using frequency gradient analysis[J]. IEEE Transactions on Nuclear Science, 2010, 57(3): 1682?1691

    9 田耕, 渠紅光, 李海濤, 等. 一種基于數(shù)字示波器的能譜和時(shí)間聯(lián)合測(cè)量系統(tǒng)[J]. 現(xiàn)代應(yīng)用物理, 2014, 5(3): 245?249

    TIAN Geng, QU Hongguang, LI Haitao, et al. A spectrum and time measurement system based on digital oscilloscope[J]. Modern Applied Physics, 2014, 5(3): 245?249

    10 Hawkes N P, Gamage K A A, Taylor G C. Digital approaches to field neutron spectrometry[J]. Radiation Measurements, 2010, 45: 1305?1308

    11 左廣霞, 何彬, 許鵬, 等. n-γ脈沖信號(hào)的仿真及其波形甄別技術(shù)研究[J]. 核技術(shù), 2012, 35(5): 365?368

    ZUO Guangxia, HE Bin, XU Peng, et al. Simulation of neutrons and gamma pulse signal and research on the pulse shape discrimination technology[J]. Nuclear Techniques, 2012, 35(5): 365?368

    12 劉大治, 李元景, 李玉蘭, 等. 粒子脈沖形狀甄別的向量投影法[J]. 高能物理與核物理, 2003, 27(11): 943?948

    LIU Dazhi, LI Yuanjing, LI Yulan, et al. Vector projection method in particles pulse shape discrimination[J]. High Energy Physics and Nuclear Physics, 2003, 27(11): 943?948

    13 羅曉亮, 劉國福, 楊俊. 基于模糊c均值聚類的液體閃爍體探測(cè)器n-γ射線甄別方法[J]. 原子能科學(xué)技術(shù), 2011, 45(6): 736?740

    LUO Xiaoliang, LIU Guofu, YANG Jun. Discrimination of neutrons and γ-rays in liquid scintillators based on fuzzy c-means clustering[J]. Atomic Energy Science and Technology, 2004, 45(6): 736?740

    14 張懷強(qiáng), 湯彬, 吳和喜. 數(shù)字核譜儀系統(tǒng)中脈沖幅度提取方法研究[J]. 核技術(shù), 2013, 36(5): 050401

    ZHANG Huaiqiang, TANG Bin, WU Hexi. Pulse amplitude extraction in digital nuclear spectrometer system[J]. Nuclear Techniques, 2013, 36(5): 050401

    15 吳治華, 趙國慶, 陸福全, 等. 原子核物理實(shí)驗(yàn)方法[M].北京: 原子能出版社, 1997: 392?393

    WU Zhihua, ZHAO Guoqing, LU Fuquan, et al. Experimental method of nuclear physics[M]. Beijing: Atomic Energy Press, 1997: 392?393

    CLC O571.53, TN152

    Digital n/γ discrimination measurement of low intensity pulsed neutron

    TIAN Geng1OUYANG Xiaoping1QU Hongguang1ZHANG Xianpeng1,2LIU Jinliang1LI Haitao1
    1(Northwest Institute of Nuclear Technology, Xi'an 710024, China)
    2(School of Nuclear Science and Technology, Xi'an Jiaotong University, Xi'an 710049, China)

    Background: The traditional measurement methods in which the detectors are working in counting mode or current mode all have limitations in the measurement of low intensity pulsed neutron. Purpose: We aim to establish a method for low intensity pulsed neutron measurement to acquire the spectra of energy and time by digitalizing and analyzing the fast current pulse generated by detector as each single neutron induced. Methods: A digital pulse shape discrimination (DPSD) system for low intensity pulsed neutron measurement has been developed, which employs wideband digital oscilloscope as data acquisition device. With BC501A liquid scintillator detector, the system can acquire and store the waveforms of neutrons and γ-rays, and discriminate neutrons from all waveforms by DPSD algorithms. The system has two operation modes as “continuous acquisition” and “acquisition window with time stamp” for different event rates according to the intensity of pulsed neutron. Results: The function of pulse height analysis of neutrons is achieved, and time information of neutron's arriving can be acquired by the analysis of the position of the waveform in the record or the time stamps. Experiment has been carried out with Am-Be neutron source with the operation mode of acquisition window, and the neutron pulse height spectrum, time spectrum and n/γ discrimination spectrum have been acquired. Conclusion: The spectra of energy and time of low intensity pulsed neutron can be measured by the digital method which employees wideband digital oscilloscope and digital signal processing algorithms, and has the advantage that all original waveforms of neutrons and γ-rays can be stored for further analysis.

    Pulsed neutron, Neutron spectrum, Time spectrum, n/γ discrimination, DPSD

    O571.53,TN152

    10.11889/j.0253-3219.2015.hjs.38.060204

    No.11305129)資助

    田耕,男,1978年出生,2007年于西北核技術(shù)研究所獲碩士學(xué)位,研究領(lǐng)域?yàn)楹穗娮訉W(xué)和脈沖射線束測(cè)量

    2015-02-10,

    2015-04-16

    猜你喜歡
    測(cè)量方法
    把握四個(gè)“三” 測(cè)量變簡(jiǎn)單
    學(xué)習(xí)方法
    滑動(dòng)摩擦力的測(cè)量和計(jì)算
    滑動(dòng)摩擦力的測(cè)量與計(jì)算
    測(cè)量的樂趣
    可能是方法不對(duì)
    用對(duì)方法才能瘦
    Coco薇(2016年2期)2016-03-22 02:42:52
    測(cè)量
    四大方法 教你不再“坐以待病”!
    Coco薇(2015年1期)2015-08-13 02:47:34
    賺錢方法
    国产伦精品一区二区三区四那| 欧美成人性av电影在线观看| 一区二区三区激情视频| xxx96com| 欧美大码av| 成人亚洲精品av一区二区| 欧美丝袜亚洲另类 | 别揉我奶头~嗯~啊~动态视频| 成人无遮挡网站| 日韩欧美免费精品| 91av网一区二区| 久久久久亚洲av毛片大全| 小蜜桃在线观看免费完整版高清| 亚洲久久久久久中文字幕| 99热6这里只有精品| www.www免费av| 香蕉丝袜av| 国产不卡一卡二| 日本成人三级电影网站| 欧美乱色亚洲激情| 成人欧美大片| 日韩大尺度精品在线看网址| 老司机午夜福利在线观看视频| 欧美中文日本在线观看视频| 桃色一区二区三区在线观看| 久久伊人香网站| 成人性生交大片免费视频hd| 免费在线观看亚洲国产| 丰满人妻熟妇乱又伦精品不卡| 色综合欧美亚洲国产小说| 一边摸一边抽搐一进一小说| 欧美激情在线99| www.色视频.com| 国产主播在线观看一区二区| 亚洲av免费在线观看| 亚洲熟妇中文字幕五十中出| 91麻豆精品激情在线观看国产| 亚洲美女视频黄频| 日韩成人在线观看一区二区三区| 美女免费视频网站| 搡女人真爽免费视频火全软件 | 一区二区三区高清视频在线| 欧美三级亚洲精品| 每晚都被弄得嗷嗷叫到高潮| 一进一出抽搐gif免费好疼| 欧美在线黄色| 午夜福利在线在线| 国产av在哪里看| 精品乱码久久久久久99久播| 男女床上黄色一级片免费看| 日韩欧美一区二区三区在线观看| 波多野结衣高清作品| 国产精品一及| 欧美xxxx黑人xx丫x性爽| 亚洲国产精品成人综合色| 亚洲欧美日韩高清专用| 嫁个100分男人电影在线观看| 最近最新免费中文字幕在线| 久久久精品欧美日韩精品| 两个人视频免费观看高清| 免费大片18禁| 最近视频中文字幕2019在线8| 日韩av在线大香蕉| 欧美zozozo另类| 露出奶头的视频| 搡女人真爽免费视频火全软件 | 韩国av一区二区三区四区| 熟妇人妻久久中文字幕3abv| 又紧又爽又黄一区二区| 人人妻,人人澡人人爽秒播| 久久婷婷人人爽人人干人人爱| 国产伦精品一区二区三区视频9 | 少妇裸体淫交视频免费看高清| 亚洲精品影视一区二区三区av| 亚洲无线在线观看| 国产精品久久久久久精品电影| 亚洲av美国av| 岛国在线免费视频观看| 欧美乱妇无乱码| 亚洲成人久久爱视频| 日本黄色片子视频| 怎么达到女性高潮| 9191精品国产免费久久| 9191精品国产免费久久| 老司机在亚洲福利影院| 国产伦精品一区二区三区四那| 韩国av一区二区三区四区| 国产精品99久久久久久久久| 日韩精品中文字幕看吧| 搡老熟女国产l中国老女人| 可以在线观看毛片的网站| 亚洲国产日韩欧美精品在线观看 | 长腿黑丝高跟| 亚洲精品美女久久久久99蜜臀| 久久久国产精品麻豆| 午夜亚洲福利在线播放| 一边摸一边抽搐一进一小说| 国产欧美日韩一区二区精品| 中文字幕精品亚洲无线码一区| 精品无人区乱码1区二区| 老司机午夜福利在线观看视频| 午夜久久久久精精品| 真实男女啪啪啪动态图| 成人18禁在线播放| av黄色大香蕉| 免费一级毛片在线播放高清视频| 在线天堂最新版资源| 村上凉子中文字幕在线| 男女做爰动态图高潮gif福利片| 噜噜噜噜噜久久久久久91| or卡值多少钱| 日本一二三区视频观看| 蜜桃亚洲精品一区二区三区| 黄色丝袜av网址大全| 免费在线观看成人毛片| 欧美中文综合在线视频| 波多野结衣巨乳人妻| 亚洲成人免费电影在线观看| 国产不卡一卡二| 97超视频在线观看视频| 丝袜美腿在线中文| 国产成人av教育| 欧美日韩亚洲国产一区二区在线观看| 最新在线观看一区二区三区| 日韩欧美国产在线观看| 免费观看精品视频网站| 麻豆一二三区av精品| 亚洲成人免费电影在线观看| 久久久久久人人人人人| 男人舔女人下体高潮全视频| 九九热线精品视视频播放| 亚洲av电影在线进入| 欧美黑人巨大hd| 精品熟女少妇八av免费久了| 欧美色欧美亚洲另类二区| 成人亚洲精品av一区二区| 国产 一区 欧美 日韩| 99久久精品一区二区三区| 成人一区二区视频在线观看| 99热这里只有是精品50| 亚洲无线观看免费| 欧美区成人在线视频| 999久久久精品免费观看国产| 成人三级黄色视频| 高清日韩中文字幕在线| 免费观看人在逋| 日韩高清综合在线| 欧美日韩一级在线毛片| 国产色爽女视频免费观看| 一级毛片高清免费大全| 亚洲中文字幕日韩| 日本黄大片高清| 日本与韩国留学比较| 久久精品夜夜夜夜夜久久蜜豆| 国产av不卡久久| 国产欧美日韩一区二区精品| av福利片在线观看| 天天一区二区日本电影三级| 欧美+亚洲+日韩+国产| 成年女人永久免费观看视频| 18禁裸乳无遮挡免费网站照片| 国产成人影院久久av| 欧美激情在线99| 国产真实乱freesex| 一进一出抽搐动态| 亚洲在线观看片| 国内精品美女久久久久久| 亚洲成人中文字幕在线播放| 日韩欧美免费精品| 97超级碰碰碰精品色视频在线观看| 国产高清videossex| 手机成人av网站| 免费av毛片视频| 久9热在线精品视频| 人人妻人人看人人澡| 国产免费一级a男人的天堂| 3wmmmm亚洲av在线观看| 国产精品久久久久久精品电影| 亚洲国产欧美人成| 99在线视频只有这里精品首页| 99热这里只有是精品50| 99国产精品一区二区三区| 国产午夜精品久久久久久一区二区三区 | 丰满人妻一区二区三区视频av | 成年女人毛片免费观看观看9| 成人欧美大片| 亚洲人成网站在线播放欧美日韩| 日韩欧美国产一区二区入口| 国产精品av视频在线免费观看| av视频在线观看入口| 一级作爱视频免费观看| 蜜桃亚洲精品一区二区三区| 国产成人av激情在线播放| 激情在线观看视频在线高清| 亚洲电影在线观看av| 日本与韩国留学比较| 超碰av人人做人人爽久久 | av黄色大香蕉| h日本视频在线播放| 91av网一区二区| 成年女人永久免费观看视频| 国产精品女同一区二区软件 | 久久精品人妻少妇| 高清日韩中文字幕在线| 国产高清有码在线观看视频| 一本久久中文字幕| 国产久久久一区二区三区| 亚洲中文字幕日韩| 日韩 欧美 亚洲 中文字幕| 变态另类丝袜制服| 亚洲成av人片免费观看| 日日摸夜夜添夜夜添小说| 色在线成人网| 亚洲狠狠婷婷综合久久图片| 欧美乱妇无乱码| 首页视频小说图片口味搜索| 亚洲av中文字字幕乱码综合| 亚洲专区中文字幕在线| 久久久国产精品麻豆| 禁无遮挡网站| 久久久精品大字幕| 小蜜桃在线观看免费完整版高清| 人妻夜夜爽99麻豆av| 99久久精品热视频| 亚洲精华国产精华精| 99久久成人亚洲精品观看| 欧美+日韩+精品| 2021天堂中文幕一二区在线观| 长腿黑丝高跟| 国产亚洲精品久久久com| 国产又黄又爽又无遮挡在线| 岛国在线观看网站| 搡女人真爽免费视频火全软件 | 欧美午夜高清在线| 亚洲性夜色夜夜综合| 免费观看的影片在线观看| 99久久精品国产亚洲精品| 欧美最新免费一区二区三区 | 淫秽高清视频在线观看| 精品久久久久久久久久免费视频| 欧美bdsm另类| 女警被强在线播放| 精品国产亚洲在线| 天天躁日日操中文字幕| 国产精品99久久久久久久久| 久久久国产精品麻豆| 中文在线观看免费www的网站| 亚洲成人精品中文字幕电影| 国产美女午夜福利| 亚洲精品在线美女| 成人无遮挡网站| 国产免费男女视频| 女生性感内裤真人,穿戴方法视频| 校园春色视频在线观看| 亚洲成人精品中文字幕电影| av女优亚洲男人天堂| 免费无遮挡裸体视频| 九九热线精品视视频播放| 国产aⅴ精品一区二区三区波| 在线视频色国产色| 99久久精品国产亚洲精品| 午夜福利在线观看免费完整高清在 | 日日干狠狠操夜夜爽| 国产伦精品一区二区三区四那| 国产精品久久久久久精品电影| 久久久国产成人精品二区| av在线蜜桃| 久久天躁狠狠躁夜夜2o2o| netflix在线观看网站| 99热只有精品国产| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 他把我摸到了高潮在线观看| 成人亚洲精品av一区二区| 国产精品影院久久| 国产成人av教育| 成年女人看的毛片在线观看| 亚洲片人在线观看| 无人区码免费观看不卡| 人妻夜夜爽99麻豆av| 国产精品久久视频播放| 亚洲色图av天堂| 国产精品女同一区二区软件 | 蜜桃久久精品国产亚洲av| av中文乱码字幕在线| 精品欧美国产一区二区三| 国语自产精品视频在线第100页| 深爱激情五月婷婷| 亚洲中文日韩欧美视频| 色综合站精品国产| 亚洲av中文字字幕乱码综合| 国产黄a三级三级三级人| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | www日本黄色视频网| 非洲黑人性xxxx精品又粗又长| 欧美日韩黄片免| 香蕉丝袜av| 免费一级毛片在线播放高清视频| 少妇的逼好多水| 日韩欧美国产在线观看| www.色视频.com| 成熟少妇高潮喷水视频| 中文字幕人成人乱码亚洲影| 国产单亲对白刺激| 午夜影院日韩av| 51国产日韩欧美| 麻豆久久精品国产亚洲av| 亚洲精华国产精华精| 亚洲最大成人中文| 51午夜福利影视在线观看| 午夜免费男女啪啪视频观看 | 搡老熟女国产l中国老女人| 亚洲精品粉嫩美女一区| 国产精品电影一区二区三区| 变态另类丝袜制服| 日日夜夜操网爽| 国产成年人精品一区二区| 一卡2卡三卡四卡精品乱码亚洲| 午夜福利在线在线| 三级毛片av免费| 全区人妻精品视频| 性色avwww在线观看| 欧美日韩精品网址| 亚洲国产精品成人综合色| 啦啦啦韩国在线观看视频| 久久久成人免费电影| 两人在一起打扑克的视频| 国产爱豆传媒在线观看| 99精品久久久久人妻精品| 日本黄色视频三级网站网址| 桃红色精品国产亚洲av| 国产视频一区二区在线看| 国产伦人伦偷精品视频| 国产熟女xx| 露出奶头的视频| 天堂av国产一区二区熟女人妻| АⅤ资源中文在线天堂| 免费在线观看日本一区| 久久香蕉国产精品| 天堂网av新在线| 高清日韩中文字幕在线| xxx96com| 岛国视频午夜一区免费看| 欧美大码av| 免费av毛片视频| 日本在线视频免费播放| 亚洲激情在线av| 国产精品 欧美亚洲| avwww免费| 欧美日韩亚洲国产一区二区在线观看| 日本 av在线| 校园春色视频在线观看| 黄色女人牲交| 精品人妻一区二区三区麻豆 | 国产亚洲欧美98| 婷婷精品国产亚洲av在线| 亚洲成人久久爱视频| 国产高清视频在线播放一区| 国产成人av教育| 欧美+亚洲+日韩+国产| 狂野欧美白嫩少妇大欣赏| 精品人妻偷拍中文字幕| 中文字幕精品亚洲无线码一区| 亚洲专区中文字幕在线| 欧美色欧美亚洲另类二区| av天堂在线播放| 无人区码免费观看不卡| 最好的美女福利视频网| 女人被狂操c到高潮| 亚洲自拍偷在线| 桃色一区二区三区在线观看| 搡老岳熟女国产| 又爽又黄无遮挡网站| 亚洲国产精品999在线| 亚洲一区高清亚洲精品| 国产一区在线观看成人免费| 亚洲人成网站在线播| 亚洲片人在线观看| 制服人妻中文乱码| 免费高清视频大片| 午夜视频国产福利| 亚洲狠狠婷婷综合久久图片| 内射极品少妇av片p| 亚洲av第一区精品v没综合| 天天一区二区日本电影三级| 国产精品 国内视频| 婷婷精品国产亚洲av| 欧美黑人巨大hd| 免费av不卡在线播放| 18禁裸乳无遮挡免费网站照片| 熟女电影av网| 久9热在线精品视频| av在线蜜桃| 成人av在线播放网站| 日韩人妻高清精品专区| 日本黄色片子视频| 一a级毛片在线观看| 一卡2卡三卡四卡精品乱码亚洲| 国内揄拍国产精品人妻在线| 国产综合懂色| 久久精品影院6| 久久精品国产清高在天天线| av片东京热男人的天堂| 精品国产超薄肉色丝袜足j| 亚洲最大成人手机在线| 在线观看av片永久免费下载| 桃色一区二区三区在线观看| 欧美bdsm另类| 久久伊人香网站| 一a级毛片在线观看| 1000部很黄的大片| 久久久久亚洲av毛片大全| 无遮挡黄片免费观看| 欧美大码av| 天堂av国产一区二区熟女人妻| 熟女电影av网| 搡女人真爽免费视频火全软件 | 国产男靠女视频免费网站| 国产精品一区二区三区四区久久| 国产精华一区二区三区| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 免费高清视频大片| 一级黄片播放器| 成人欧美大片| 99视频精品全部免费 在线| www日本黄色视频网| 免费高清视频大片| 久9热在线精品视频| 欧美最黄视频在线播放免费| 男人舔奶头视频| 99久久无色码亚洲精品果冻| 国产成人av激情在线播放| 成人国产综合亚洲| 亚洲成人精品中文字幕电影| 一区二区三区免费毛片| 国产高清三级在线| 国产三级中文精品| 亚洲性夜色夜夜综合| 中国美女看黄片| 久久国产精品人妻蜜桃| 精品国产超薄肉色丝袜足j| 国产一区二区激情短视频| 亚洲人与动物交配视频| 国产毛片a区久久久久| 天堂√8在线中文| 老司机福利观看| 成人av一区二区三区在线看| 亚洲精品美女久久久久99蜜臀| 午夜精品在线福利| 欧美乱码精品一区二区三区| 免费av不卡在线播放| 日韩免费av在线播放| 性色av乱码一区二区三区2| 一本精品99久久精品77| 亚洲乱码一区二区免费版| 亚洲狠狠婷婷综合久久图片| xxx96com| 久久国产乱子伦精品免费另类| 欧美黑人巨大hd| 日日夜夜操网爽| 久久久久亚洲av毛片大全| 一区二区三区免费毛片| 日日夜夜操网爽| 一a级毛片在线观看| 在线观看日韩欧美| 欧美绝顶高潮抽搐喷水| 丰满人妻一区二区三区视频av | 狂野欧美激情性xxxx| 免费av不卡在线播放| 男女做爰动态图高潮gif福利片| 日韩欧美精品免费久久 | 99国产精品一区二区蜜桃av| 免费在线观看亚洲国产| 欧美不卡视频在线免费观看| 五月玫瑰六月丁香| 最新在线观看一区二区三区| 黄色日韩在线| 国产一区二区在线av高清观看| 国产亚洲欧美98| 小说图片视频综合网站| 99在线人妻在线中文字幕| 国产精品亚洲av一区麻豆| 色老头精品视频在线观看| 草草在线视频免费看| 亚洲欧美日韩高清专用| 亚洲一区二区三区不卡视频| 宅男免费午夜| 狂野欧美激情性xxxx| 免费在线观看影片大全网站| 老鸭窝网址在线观看| 久久精品亚洲精品国产色婷小说| 深夜精品福利| 狂野欧美激情性xxxx| 精品国产亚洲在线| 久久久久久大精品| 欧美另类亚洲清纯唯美| 亚洲中文字幕日韩| 亚洲av一区综合| 男女那种视频在线观看| 亚洲国产精品999在线| 禁无遮挡网站| 99国产精品一区二区三区| avwww免费| 很黄的视频免费| 在线播放无遮挡| 日韩欧美三级三区| 国产成人欧美在线观看| 亚洲成人免费电影在线观看| 国内毛片毛片毛片毛片毛片| 怎么达到女性高潮| 国产一区二区在线观看日韩 | 亚洲欧美日韩高清专用| 国产亚洲精品综合一区在线观看| 欧美黑人欧美精品刺激| 欧美另类亚洲清纯唯美| 久久国产乱子伦精品免费另类| 国产麻豆成人av免费视频| 国产极品精品免费视频能看的| 69av精品久久久久久| 成人一区二区视频在线观看| av欧美777| 免费在线观看成人毛片| 亚洲成人免费电影在线观看| 十八禁网站免费在线| e午夜精品久久久久久久| 欧美日韩瑟瑟在线播放| 天堂√8在线中文| 观看免费一级毛片| 久久精品夜夜夜夜夜久久蜜豆| 国产单亲对白刺激| 不卡一级毛片| 亚洲精华国产精华精| 亚洲一区二区三区不卡视频| 欧美一区二区国产精品久久精品| 国产精品一区二区免费欧美| 九九久久精品国产亚洲av麻豆| xxxwww97欧美| 两个人的视频大全免费| 久久久久国内视频| 国产乱人伦免费视频| 欧美日韩乱码在线| 中文字幕av在线有码专区| 亚洲精品美女久久久久99蜜臀| 国产又黄又爽又无遮挡在线| 久久精品国产自在天天线| 国内精品久久久久久久电影| 九九在线视频观看精品| 欧美一区二区亚洲| 亚洲成av人片在线播放无| 天堂√8在线中文| 国内精品久久久久精免费| 97超视频在线观看视频| 国产又黄又爽又无遮挡在线| 日韩欧美三级三区| 国产成人啪精品午夜网站| 久久欧美精品欧美久久欧美| 99在线视频只有这里精品首页| 少妇熟女aⅴ在线视频| 18禁在线播放成人免费| av专区在线播放| 亚洲av免费在线观看| 国产高清videossex| 久久99热这里只有精品18| 午夜福利成人在线免费观看| 国产成人系列免费观看| 首页视频小说图片口味搜索| 97碰自拍视频| 日韩亚洲欧美综合| 精品一区二区三区视频在线观看免费| 嫩草影院入口| 欧美不卡视频在线免费观看| 国产精品 国内视频| 一级黄色大片毛片| 岛国在线免费视频观看| 91在线观看av| 有码 亚洲区| 久久久精品大字幕| 国产精品,欧美在线| 亚洲不卡免费看| 欧美在线黄色| 99久国产av精品| 午夜福利18| 亚洲一区高清亚洲精品| 波多野结衣高清无吗| 久99久视频精品免费| 日本黄大片高清| 国产亚洲精品久久久com| 全区人妻精品视频| 最后的刺客免费高清国语| h日本视频在线播放| 免费人成视频x8x8入口观看| 亚洲国产精品合色在线| 51午夜福利影视在线观看| 搞女人的毛片| 国产精品自产拍在线观看55亚洲| 18+在线观看网站| a级毛片a级免费在线| 免费av观看视频| 一个人免费在线观看的高清视频| 国产探花在线观看一区二区| 久久中文看片网| 熟女少妇亚洲综合色aaa.| 美女免费视频网站| 亚洲av电影不卡..在线观看| 国产一区二区在线观看日韩 | 一级毛片女人18水好多| 99在线人妻在线中文字幕| av专区在线播放| 国产精品爽爽va在线观看网站| 波多野结衣巨乳人妻| 天天添夜夜摸| 熟妇人妻久久中文字幕3abv| 观看美女的网站| 久久精品亚洲精品国产色婷小说| 两个人视频免费观看高清| 黄色日韩在线|