• <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歐陽(yáng)曉平渠紅光張顯鵬劉金良李海濤
    核技術(shù) 2015年6期
    關(guān)鍵詞:測(cè)量方法

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

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

    田 耕1歐陽(yáng)曉平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è)器中形成的電流脈沖大量堆積,無(wú)法區(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)和國(guó)內(nèi)的HL-2A等裝置的放電過程中D-D及D-T聚變中子測(cè)量[2?3]。而對(duì)于強(qiáng)度較低的中子束,有效事件占總時(shí)間的比例較低,用連續(xù)記錄方法測(cè)量到的數(shù)據(jù)中會(huì)包含大量無(wú)用的基線數(shù)據(jù),浪費(fèi)了存儲(chǔ)空間,使得能夠獲取的有效事件數(shù)目減小。在這種情況下,更為有效的采集方法是采集窗記錄方式。該方法在有符合記錄條件的信號(hào)到來時(shí)進(jìn)行采集并存儲(chǔ)一定長(zhǎng)度的數(shù)據(jù),而在無(wú)有效信號(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í)陽(yáng)極輸出脈沖信號(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í)間,無(wú)論是在連續(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)選用美國(guó)泰克公司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無(wú)機(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算法

    目前國(guó)內(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的精度,無(wú)法直接利用軟件開發(fā)工具提供的函數(shù)庫(kù)進(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方法,無(wú)法有效區(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é)語(yǔ)

    采用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 羅曉亮, 劉國(guó)福, 楊俊. 基于模糊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 吳治華, 趙國(guó)慶, 陸福全, 等. 原子核物理實(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新网站| 国内揄拍国产精品人妻在线| 又粗又硬又长又爽又黄的视频| 久久国产亚洲av麻豆专区| 欧美精品亚洲一区二区| 老师上课跳d突然被开到最大视频| 一级毛片aaaaaa免费看小| 亚洲色图av天堂| 日韩人妻高清精品专区| www.色视频.com| 99九九线精品视频在线观看视频| 国产精品一区二区三区四区免费观看| 亚洲精品乱码久久久久久按摩| 国产成人免费观看mmmm| 18禁在线无遮挡免费观看视频| 久久99热6这里只有精品| 国产免费又黄又爽又色| 日本-黄色视频高清免费观看| 久久精品国产自在天天线| 亚洲精品国产色婷婷电影| 国产精品久久久久成人av| 久久午夜福利片| 妹子高潮喷水视频| 国产av一区二区精品久久 | 国产国拍精品亚洲av在线观看| 国产精品久久久久久av不卡| 如何舔出高潮| 高清视频免费观看一区二区| 欧美精品一区二区免费开放| 黑人猛操日本美女一级片| 丰满少妇做爰视频| 日本黄色片子视频| 午夜精品国产一区二区电影| 日韩国内少妇激情av| 成人高潮视频无遮挡免费网站| 欧美亚洲 丝袜 人妻 在线| 夫妻午夜视频| 精品一区二区三区视频在线| 国产亚洲av片在线观看秒播厂| 国产成人精品福利久久| 五月伊人婷婷丁香| 久久ye,这里只有精品| a 毛片基地| 免费少妇av软件| 夫妻性生交免费视频一级片| 国内揄拍国产精品人妻在线| 丰满少妇做爰视频| videos熟女内射| 性高湖久久久久久久久免费观看| 国产伦精品一区二区三区视频9| 偷拍熟女少妇极品色| 国产美女午夜福利| 欧美高清成人免费视频www| 久久久久久九九精品二区国产| 久久久久视频综合| 亚洲四区av| 国产精品麻豆人妻色哟哟久久| 最近中文字幕2019免费版| 国内揄拍国产精品人妻在线| 一区二区三区免费毛片| 精品亚洲成国产av| 午夜福利在线观看免费完整高清在| 亚洲av男天堂| 美女cb高潮喷水在线观看| av黄色大香蕉| 国产有黄有色有爽视频| 亚洲国产av新网站| 最近手机中文字幕大全| 最近最新中文字幕大全电影3| 亚洲成人av在线免费| 亚洲欧洲国产日韩| 一区二区三区乱码不卡18| 亚洲人成网站在线播| 亚洲欧美日韩东京热| 国产一区二区三区av在线| 噜噜噜噜噜久久久久久91| 少妇熟女欧美另类| 在线精品无人区一区二区三 | 亚洲美女搞黄在线观看| 成人高潮视频无遮挡免费网站| 成人黄色视频免费在线看| 美女xxoo啪啪120秒动态图| 最近的中文字幕免费完整| 久久这里有精品视频免费| 男女啪啪激烈高潮av片| 亚洲国产精品一区三区| 久久久久久久国产电影| 大陆偷拍与自拍| 高清日韩中文字幕在线| 国内精品宾馆在线| av专区在线播放| 国产永久视频网站| 搡老乐熟女国产| 成年免费大片在线观看| 不卡视频在线观看欧美| 亚洲国产av新网站| 久久久久网色| 亚洲精品一二三| 久久99精品国语久久久| 最近中文字幕高清免费大全6| 王馨瑶露胸无遮挡在线观看| 纯流量卡能插随身wifi吗| 久久久精品免费免费高清| 亚洲国产日韩一区二区| 成人国产麻豆网| 男男h啪啪无遮挡| 久久久久久久国产电影| 熟女电影av网| 日本猛色少妇xxxxx猛交久久| 乱系列少妇在线播放| 亚洲人成网站高清观看| 中文天堂在线官网| 久久久久人妻精品一区果冻| 亚洲成人一二三区av| 久久精品国产亚洲av涩爱| 99热国产这里只有精品6| 91精品国产国语对白视频| 国产久久久一区二区三区| 国产在线免费精品| 各种免费的搞黄视频| av不卡在线播放| 日韩电影二区| 亚洲精品久久久久久婷婷小说| av一本久久久久| 少妇熟女欧美另类| 岛国毛片在线播放| 亚洲综合色惰| 国产 精品1| 欧美变态另类bdsm刘玥| 哪个播放器可以免费观看大片| 人人妻人人添人人爽欧美一区卜 | 91久久精品电影网| 99久久中文字幕三级久久日本| 欧美性感艳星| 一级a做视频免费观看| 国产成人免费观看mmmm| 免费大片18禁| 久久久久国产精品人妻一区二区| 大又大粗又爽又黄少妇毛片口| 成人免费观看视频高清| 色5月婷婷丁香| 夜夜看夜夜爽夜夜摸| 观看av在线不卡| 日本爱情动作片www.在线观看| 国产视频内射| 久久亚洲国产成人精品v| 身体一侧抽搐| 人人妻人人澡人人爽人人夜夜| 亚洲美女黄色视频免费看| 超碰av人人做人人爽久久| 最近2019中文字幕mv第一页| 在线亚洲精品国产二区图片欧美 | 一区二区三区四区激情视频| 在线亚洲精品国产二区图片欧美 | 在现免费观看毛片| 亚洲第一av免费看| 免费观看av网站的网址| 一本色道久久久久久精品综合| 欧美一区二区亚洲| 久久精品夜色国产| 男人添女人高潮全过程视频| 一级毛片黄色毛片免费观看视频| 欧美bdsm另类| 中国国产av一级| 黄色怎么调成土黄色| 高清日韩中文字幕在线| 一区二区三区四区激情视频| 在线免费观看不下载黄p国产| 亚洲不卡免费看| 国产精品一区二区三区四区免费观看| 一区二区三区精品91| 久久久国产一区二区| 亚洲综合精品二区| 内地一区二区视频在线| 91狼人影院| av在线观看视频网站免费| 777米奇影视久久| 国产免费视频播放在线视频| 妹子高潮喷水视频| 日韩欧美一区视频在线观看 | 国产一区二区三区综合在线观看 | 国国产精品蜜臀av免费| 国内精品宾馆在线| 精品人妻一区二区三区麻豆| 亚洲中文av在线| 国产精品偷伦视频观看了| 一边亲一边摸免费视频| 毛片女人毛片| 国产日韩欧美亚洲二区| 看免费成人av毛片| 欧美日韩视频精品一区| 青青草视频在线视频观看| 国产精品一二三区在线看| 联通29元200g的流量卡| 精品熟女少妇av免费看| 一级毛片 在线播放| 18禁裸乳无遮挡免费网站照片| 国产乱人偷精品视频| 免费av不卡在线播放| 亚洲美女视频黄频| 毛片一级片免费看久久久久| 久久久久久久国产电影| 国产在线免费精品| av国产免费在线观看| 嫩草影院新地址| 成人二区视频| 亚洲无线观看免费| 日产精品乱码卡一卡2卡三| 亚洲精品国产色婷婷电影| 久久精品夜色国产| 人人妻人人看人人澡| 午夜日本视频在线| 伦理电影免费视频| 国产精品国产三级国产专区5o| 中国美白少妇内射xxxbb| 久久精品人妻少妇| 欧美国产精品一级二级三级 | 久久久久久久国产电影| 伦理电影大哥的女人| 男女下面进入的视频免费午夜| 一二三四中文在线观看免费高清| 特大巨黑吊av在线直播| 不卡视频在线观看欧美| 亚洲成人中文字幕在线播放| 国产av国产精品国产| 亚洲精品乱码久久久v下载方式| 少妇高潮的动态图| 免费播放大片免费观看视频在线观看| 91久久精品国产一区二区成人| 特大巨黑吊av在线直播| 女性生殖器流出的白浆| 又大又黄又爽视频免费| 亚洲高清免费不卡视频| 国产91av在线免费观看| 伦理电影大哥的女人| 又大又黄又爽视频免费| 一区二区三区免费毛片| 国产色爽女视频免费观看| 制服丝袜香蕉在线| 国产亚洲91精品色在线| a 毛片基地| 久久久欧美国产精品| 在线观看免费日韩欧美大片 | 午夜福利在线在线| 成人影院久久| 免费看av在线观看网站| 国产在线男女| 亚洲一区二区三区欧美精品| 舔av片在线| 婷婷色av中文字幕| 女人久久www免费人成看片| 色婷婷av一区二区三区视频| 中文字幕精品免费在线观看视频 | 日产精品乱码卡一卡2卡三| 成人18禁高潮啪啪吃奶动态图 | kizo精华| 欧美精品一区二区大全| 国国产精品蜜臀av免费| 三级国产精品片| 亚洲在久久综合| 日韩不卡一区二区三区视频在线| 一区二区三区免费毛片| tube8黄色片| 在线看a的网站| 麻豆国产97在线/欧美| 丝袜脚勾引网站| 亚洲三级黄色毛片| 久久热精品热| 成人漫画全彩无遮挡| 如何舔出高潮| 老司机影院毛片| av又黄又爽大尺度在线免费看| 久久久久国产精品人妻一区二区| 国产伦在线观看视频一区| 在线观看一区二区三区| videossex国产| 亚洲第一av免费看| 亚洲精品亚洲一区二区| 青青草视频在线视频观看| 蜜桃久久精品国产亚洲av| 夜夜看夜夜爽夜夜摸| 久久久国产一区二区| 久久人人爽人人片av| 国产成人一区二区在线| 亚洲精品日韩av片在线观看| 女人十人毛片免费观看3o分钟| 亚洲国产精品国产精品| 99视频精品全部免费 在线| 国产亚洲91精品色在线| 乱系列少妇在线播放| 国产伦理片在线播放av一区| 自拍偷自拍亚洲精品老妇| 一级毛片久久久久久久久女| 极品教师在线视频| 三级经典国产精品| 色吧在线观看| 在线 av 中文字幕| 秋霞伦理黄片| 国产精品偷伦视频观看了| 久久97久久精品| 99久久人妻综合| 日本黄大片高清| 色哟哟·www| 久久久久国产精品人妻一区二区| 亚洲性久久影院| 又大又黄又爽视频免费| 春色校园在线视频观看| av专区在线播放| 久久99蜜桃精品久久| 久久人人爽av亚洲精品天堂 | 亚洲美女黄色视频免费看| 卡戴珊不雅视频在线播放| 新久久久久国产一级毛片| 高清黄色对白视频在线免费看 | 少妇猛男粗大的猛烈进出视频| 日日啪夜夜爽| 在线观看av片永久免费下载| 偷拍熟女少妇极品色| 国产在线一区二区三区精| 国产成人精品久久久久久| 青春草亚洲视频在线观看| 高清毛片免费看| 你懂的网址亚洲精品在线观看| 精品久久久久久久久av| 麻豆乱淫一区二区| 色综合色国产| 黑人高潮一二区| 少妇熟女欧美另类| 大香蕉久久网| 18禁裸乳无遮挡动漫免费视频| 午夜激情福利司机影院| 欧美激情国产日韩精品一区| 久久99热这里只有精品18| 亚洲av免费高清在线观看| 国产精品国产av在线观看| 国产 一区精品| 国精品久久久久久国模美| 国产一区有黄有色的免费视频| 国产精品偷伦视频观看了| 嫩草影院入口| 久久精品熟女亚洲av麻豆精品| 蜜臀久久99精品久久宅男| 最近中文字幕高清免费大全6| 国产日韩欧美亚洲二区| 国产高潮美女av| av在线播放精品| 少妇精品久久久久久久| 久久久亚洲精品成人影院| av在线播放精品| av在线老鸭窝| 国产精品国产三级专区第一集| 女人十人毛片免费观看3o分钟| 亚洲人成网站高清观看| 亚洲精品视频女| 国产高清有码在线观看视频| tube8黄色片| 亚洲熟女精品中文字幕| 五月天丁香电影| 黄色一级大片看看| 国产精品人妻久久久影院| 少妇 在线观看| 99视频精品全部免费 在线| 免费大片18禁| 亚洲国产毛片av蜜桃av| 国产亚洲5aaaaa淫片| 精品久久久精品久久久| 欧美高清成人免费视频www| 欧美bdsm另类| 亚洲国产色片| 免费大片黄手机在线观看| 久久久国产一区二区| 亚洲国产成人一精品久久久| 交换朋友夫妻互换小说| 亚洲aⅴ乱码一区二区在线播放| 亚洲三级黄色毛片| 五月伊人婷婷丁香| 国产片特级美女逼逼视频| 一个人看视频在线观看www免费| 国产精品人妻久久久影院| 91精品国产九色| 交换朋友夫妻互换小说| 国产免费视频播放在线视频| 久热这里只有精品99| 亚洲伊人久久精品综合| 午夜激情福利司机影院| 大陆偷拍与自拍| 中国国产av一级| 99热这里只有是精品50| 亚洲天堂av无毛| 男人爽女人下面视频在线观看| 十八禁网站网址无遮挡 | 人人妻人人添人人爽欧美一区卜 | 国产精品国产三级专区第一集| 日韩免费高清中文字幕av| 国产成人一区二区在线| 日韩三级伦理在线观看| 亚洲av不卡在线观看| 国产女主播在线喷水免费视频网站| 日日摸夜夜添夜夜爱| 人妻 亚洲 视频| 日韩制服骚丝袜av| 亚洲国产精品一区三区| 国产乱人偷精品视频| 免费高清在线观看视频在线观看| 91久久精品国产一区二区成人| 美女中出高潮动态图| 亚洲av成人精品一二三区| 久久精品夜色国产| 熟女电影av网| 久久青草综合色| 国产午夜精品一二区理论片| 日本与韩国留学比较| 最近手机中文字幕大全| 欧美+日韩+精品| xxx大片免费视频| 精品99又大又爽又粗少妇毛片| 午夜日本视频在线| 能在线免费看毛片的网站| 91精品一卡2卡3卡4卡| 黑人高潮一二区| 老司机影院成人| 街头女战士在线观看网站| 国产成人免费无遮挡视频| 国产精品久久久久久精品电影小说 | 久久久亚洲精品成人影院| 青青草视频在线视频观看| 亚洲av综合色区一区| videos熟女内射| 久久国产精品大桥未久av | 国产精品一区二区三区四区免费观看| 自拍欧美九色日韩亚洲蝌蚪91 | 午夜免费鲁丝| 亚洲精品一区蜜桃| 伦精品一区二区三区| 国产爱豆传媒在线观看| 亚洲欧洲日产国产| 亚洲高清免费不卡视频| 精品亚洲成a人片在线观看 | 日韩伦理黄色片| 午夜老司机福利剧场| 老师上课跳d突然被开到最大视频| 日韩国内少妇激情av| 亚洲精品亚洲一区二区| 爱豆传媒免费全集在线观看| 久久久久人妻精品一区果冻| 欧美xxxx黑人xx丫x性爽| 久久精品国产亚洲av天美| 精品少妇黑人巨大在线播放| 九九久久精品国产亚洲av麻豆| 午夜老司机福利剧场| 一个人免费看片子| 大陆偷拍与自拍| 最新中文字幕久久久久| 亚洲色图综合在线观看| 亚洲av在线观看美女高潮| 人体艺术视频欧美日本| 爱豆传媒免费全集在线观看| 欧美丝袜亚洲另类| kizo精华| 亚洲欧美精品专区久久| 噜噜噜噜噜久久久久久91| 乱系列少妇在线播放| 精品人妻一区二区三区麻豆| 亚洲色图av天堂| 街头女战士在线观看网站| 视频中文字幕在线观看| 香蕉精品网在线| 九九在线视频观看精品| 哪个播放器可以免费观看大片| 精品午夜福利在线看| 视频中文字幕在线观看| 国产在视频线精品| 美女高潮的动态| 热re99久久精品国产66热6| 韩国av在线不卡| tube8黄色片| www.av在线官网国产| 亚洲真实伦在线观看| 免费av中文字幕在线| 国产精品久久久久成人av| 欧美精品亚洲一区二区| 干丝袜人妻中文字幕| 国产白丝娇喘喷水9色精品| 色网站视频免费| 久久久久视频综合| 成年av动漫网址| 欧美三级亚洲精品| 青春草亚洲视频在线观看| 狠狠精品人妻久久久久久综合| 亚洲电影在线观看av| 久久人人爽人人爽人人片va| 免费看日本二区| 夜夜爽夜夜爽视频| 亚洲天堂av无毛| 又大又黄又爽视频免费| 免费不卡的大黄色大毛片视频在线观看| 天堂8中文在线网| 亚洲国产精品国产精品| 久久久亚洲精品成人影院| 免费观看无遮挡的男女| 极品教师在线视频| 成人免费观看视频高清| 日日撸夜夜添| 久久人人爽人人片av| 青春草视频在线免费观看| 人妻夜夜爽99麻豆av| 五月伊人婷婷丁香| 亚洲av电影在线观看一区二区三区| av女优亚洲男人天堂| 久热这里只有精品99| 亚洲欧美一区二区三区黑人 | 国产男女超爽视频在线观看| 九九久久精品国产亚洲av麻豆| 国产精品人妻久久久影院| 在线观看美女被高潮喷水网站| 51国产日韩欧美| 在线观看免费视频网站a站| 亚洲经典国产精华液单| 97在线视频观看| 九九久久精品国产亚洲av麻豆| 国产精品嫩草影院av在线观看| 少妇熟女欧美另类| 国产午夜精品一二区理论片| 2022亚洲国产成人精品| 色视频在线一区二区三区| 亚洲精品一区蜜桃| 寂寞人妻少妇视频99o| 波野结衣二区三区在线| 日韩av在线免费看完整版不卡| 国产黄片视频在线免费观看| 精品亚洲成国产av| 在线天堂最新版资源| 插逼视频在线观看| 免费人妻精品一区二区三区视频| 日本vs欧美在线观看视频 | 国产精品不卡视频一区二区| 欧美变态另类bdsm刘玥| 黄色一级大片看看| 亚洲精品视频女| 亚洲国产欧美在线一区| 久久99热这里只频精品6学生| 国产av国产精品国产| 久久 成人 亚洲| av在线app专区| 欧美日韩综合久久久久久| 国产成人精品婷婷| 国产成人a∨麻豆精品| 美女高潮的动态| 极品少妇高潮喷水抽搐| 纵有疾风起免费观看全集完整版| 亚洲av成人精品一二三区| 91精品国产国语对白视频| 国产爽快片一区二区三区| 一级毛片aaaaaa免费看小| 国产黄频视频在线观看| 激情 狠狠 欧美| 国产色爽女视频免费观看| 国产免费福利视频在线观看| 国产av国产精品国产| 九九爱精品视频在线观看| 欧美三级亚洲精品| 日本猛色少妇xxxxx猛交久久| 精品亚洲成国产av| 2022亚洲国产成人精品| 久久午夜福利片| 日韩,欧美,国产一区二区三区| 亚洲自偷自拍三级| 欧美日韩亚洲高清精品| 观看免费一级毛片| 纵有疾风起免费观看全集完整版| 成人影院久久| 国产 精品1| 黄色日韩在线| 黑人猛操日本美女一级片| 国产精品99久久99久久久不卡 | 校园人妻丝袜中文字幕| 老女人水多毛片| 少妇人妻一区二区三区视频| 少妇猛男粗大的猛烈进出视频| 最近的中文字幕免费完整| 国产黄片美女视频| 国产精品一区二区在线不卡| av黄色大香蕉| 最新中文字幕久久久久| 中文精品一卡2卡3卡4更新| 日本猛色少妇xxxxx猛交久久| 日本av手机在线免费观看| 丝瓜视频免费看黄片| 日韩av不卡免费在线播放| 精品人妻一区二区三区麻豆| 三级国产精品欧美在线观看| 久久精品久久精品一区二区三区| 久久久久久久亚洲中文字幕| 高清黄色对白视频在线免费看 | 一级黄片播放器| 欧美日韩综合久久久久久| 一级二级三级毛片免费看| 人人妻人人爽人人添夜夜欢视频 | 精品一区在线观看国产| 80岁老熟妇乱子伦牲交| 爱豆传媒免费全集在线观看| 乱码一卡2卡4卡精品| 人人妻人人澡人人爽人人夜夜| 黑人猛操日本美女一级片| 又爽又黄a免费视频| 哪个播放器可以免费观看大片| 久久久久精品性色| 国产黄片美女视频| 亚洲精品久久午夜乱码| 精品一区二区三卡| 深爱激情五月婷婷| 人人妻人人看人人澡| 日本av免费视频播放| 熟女人妻精品中文字幕| 中文资源天堂在线|