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

    瞬發(fā)伽馬活化分析與中子層析照相聯(lián)合測(cè)量技術(shù)

    2017-08-16 09:30:27李潤(rùn)東王冠博
    同位素 2017年3期
    關(guān)鍵詞:中子源束流層析

    楊 鑫,李潤(rùn)東,王冠博,唐 彬,孫 勇

    (中國(guó)工程物理研究院 核物理與化學(xué)研究所,四川 綿陽 621900)

    瞬發(fā)伽馬活化分析與中子層析照相聯(lián)合測(cè)量技術(shù)

    楊 鑫,李潤(rùn)東,王冠博,唐 彬,孫 勇

    (中國(guó)工程物理研究院 核物理與化學(xué)研究所,四川 綿陽 621900)

    介紹瞬發(fā)伽馬(中子)活化分析(prompt gamma-ray neutron activation analysis, PGAA/PGNAA)技術(shù)的發(fā)展應(yīng)用情況與國(guó)際上的最新研究成果。將PGAA與中子層析照相(neutron tomography, NT)技術(shù)結(jié)合,利用中子層析照相重建結(jié)果,得到樣品內(nèi)部三維結(jié)構(gòu),研究中子束流在樣品內(nèi)部的衰減過程,確定樣品內(nèi)部特定位置處的中子束流強(qiáng)度,獲得其強(qiáng)度分布。結(jié)合PGAA與NT,可以實(shí)現(xiàn)不對(duì)樣品拆解取樣,對(duì)其內(nèi)部元素成分進(jìn)行三維分布分析。國(guó)際上,德國(guó)FRM Ⅱ研究堆與匈牙利布達(dá)佩斯研究堆的中子束流上已建成相關(guān)裝置,并在天文、考古領(lǐng)域得到應(yīng)用。國(guó)內(nèi)相關(guān)研究也相繼啟動(dòng)。

    瞬發(fā)伽馬活化分析;中子層析;無損檢測(cè);圖像重建

    瞬發(fā)伽馬(中子)活化分析(prompt gamma-ray neutron activation analysis, PGAA/PGNAA)作為一種標(biāo)準(zhǔn)的材料元素成分分析技術(shù),廣泛應(yīng)用于化學(xué)研究、工業(yè)技術(shù)、環(huán)境保護(hù)、生物醫(yī)療等領(lǐng)域[1-2]。受樣品自身對(duì)中子與γ射線的散射、吸收等效應(yīng)的影響,傳統(tǒng)的PGAA多應(yīng)用于小樣品或淺層分析。隨著科學(xué)技術(shù)的發(fā)展,對(duì)分析技術(shù)的要求逐漸提高,在文物保護(hù)、隕石分析、武器分析、生物活體研究等方面需要不對(duì)樣品進(jìn)行拆解取樣既可進(jìn)行分析。近年來,德國(guó)[3]、匈牙利[4]的研究機(jī)構(gòu)在歐盟“Ancient Charm”計(jì)劃的支持下,發(fā)展了將PGAA與中子層析技術(shù)(neutron tomography, NT)相結(jié)合,研究中子束流在樣品內(nèi)部的衰減過程,以確定樣品內(nèi)部特定位置處的中子束流強(qiáng)度,解決PGAA無法深入內(nèi)部測(cè)量中子束流強(qiáng)度分布的問題。建成了PGAA-NT聯(lián)合測(cè)量裝置,并已應(yīng)用于天文、考古領(lǐng)域,取得了良好的效果。國(guó)內(nèi)研究機(jī)構(gòu)也緊跟國(guó)際發(fā)展前沿,啟動(dòng)了相關(guān)研究。

    1 測(cè)量方法

    1.1 瞬發(fā)伽馬活化分析

    瞬發(fā)伽馬活化分析(prompt gamma-ray neutron activation analysis, PGAA/PGNAA)是一種元素含量無損檢測(cè)技術(shù)?;谠雍说乃舶l(fā)中子輻射俘獲(n,γ)反應(yīng)(也可用于半衰期較短的緩發(fā)伽馬射線),通過測(cè)量受熱(冷)中子束流照射后的特征伽馬射線強(qiáng)度,以確定相應(yīng)核素的含量。PGAA檢測(cè)核素含量靈敏度受(n,γ)反應(yīng)的截面影響,對(duì)截面較大的元素更為靈敏,對(duì)H的測(cè)量靈敏度可達(dá)微克量級(jí)。除截面外,其靈敏度與特征峰的本底、計(jì)數(shù)統(tǒng)計(jì)不確定度等有關(guān)[1-2]。對(duì)某一特定元素,理想狀態(tài)下每克產(chǎn)生的特征伽馬射線計(jì)數(shù)率[1]:

    (1)

    式中,NA為阿伏伽德羅常數(shù);Q為俘獲中子的核素的同位素豐度;σ為中子俘獲截面,b;Ф為中子注量率,cm-2·s-1;Γ為伽馬射線產(chǎn)額;ε(E)為能量E的γ射線的探測(cè)效率;M為相對(duì)原子質(zhì)量。

    根據(jù)樣品測(cè)量的計(jì)數(shù)率與S比較,即可得到待測(cè)核素的含量。其中σ、Ф與中子束流的性質(zhì)密切相關(guān)。而ε(E)與樣品幾何、測(cè)量幾何、探測(cè)器自身效率等相關(guān),需要利用刻錄、標(biāo)準(zhǔn)比對(duì)等方法處理。

    典型的PGAA裝置[5](圖 1)主要包含中子源、中子束流導(dǎo)管與整形系統(tǒng)、輻照靶位、樣品支架、γ射線探測(cè)器等元器件,以及各元器件的中子/γ屏蔽結(jié)構(gòu)。其裝置設(shè)計(jì)涉及各部件間的密切、有效配合,以使裝置性能最優(yōu)化。中子源的選擇會(huì)影響PGAA裝置的性能,目前的中子源分為同位素中子源、加速器中子源和反應(yīng)堆中子源。其中,基于同位素或小型加速器中子源的裝置安裝更為便捷、靈活,可設(shè)計(jì)為移動(dòng)式裝置。但由于其中子能量相對(duì)較高,強(qiáng)度偏低,即使經(jīng)過慢化處理,其束流性能也不理想,限制了PGAA裝置的性能。而反應(yīng)堆中子束流能量低、強(qiáng)度高、信噪比好,特別是冷中子源技術(shù)的發(fā)展,使得其束流性能更佳,是進(jìn)行PGAA研究的理想中子源。

    由于PGAA是在線帶束測(cè)量,其γ射線能譜較為復(fù)雜,某些重核素的γ射線多達(dá)成百上千條,且其中包含了中子散射后形成的噪聲信號(hào),需要進(jìn)行大量復(fù)雜的譜分析與解譜研究。

    1.2 PGAA發(fā)展及應(yīng)用

    自1932年中子發(fā)現(xiàn)以來,中子及其引發(fā)的核反應(yīng)被廣泛研究。1934年,Amaldi、Lea等發(fā)現(xiàn)含氫材料與中子發(fā)生俘獲反應(yīng)時(shí)會(huì)放出高穿透力的γ射線,該射線就是1H(n,γ)2H反應(yīng)釋放的能量為2 223.28 keV的γ射線。1936年,Gy?rgy Hevesy提出通過測(cè)量受中子輻照生成的放射性核素進(jìn)行元素分析的方法,這是最早的中子活化分析思想[2]。1966年,Isenhour與Morrison在Cornell大學(xué)建成了首臺(tái)基于反應(yīng)堆[6]、采用NaI伽馬射線探測(cè)器的PGAA裝置。上世紀(jì)60年代末,隨著鍺(Ge)半導(dǎo)體探測(cè)器技術(shù)的發(fā)展,使探測(cè)伽馬射線的能量分辨率相對(duì)NaI探測(cè)器有了數(shù)十倍的提高,特別是反康普頓技術(shù)的應(yīng)用,促進(jìn)了PGAA技術(shù)的發(fā)展。

    上世紀(jì)80年代,大體積Ge探測(cè)器的出現(xiàn)改善了探測(cè)高能伽馬射線的效率,掀起全球范圍內(nèi)的PGAA裝置建設(shè)熱潮,各實(shí)驗(yàn)室都開展了自建PGAA裝置對(duì)各元素分析靈敏度的研究。美國(guó)NIST[7]、密蘇里大學(xué)[8]、日本JAERI建成了基于準(zhǔn)直中子束流的PGAA裝置,由導(dǎo)管引出的低能中子用于PGAA研究的靈敏度高、本底低[9]。90年代,隨著高強(qiáng)度熱(冷)中子束流的出現(xiàn),美國(guó)NIST[5]、德國(guó)KFA[10]、日本JAERI[11]以及匈牙利布達(dá)佩斯中子中心(BNC)[12]都建成了基于高強(qiáng)度中子束流的PGAA裝置。

    部分基于反應(yīng)堆的PGAA裝置信息列于表1[2]。

    PGAA具有無損、適用范圍廣、靈敏度高、精度高等特點(diǎn),與傳統(tǒng)的反應(yīng)堆中子活化分析(reactor neutron activation analysis, ReNAA)相互補(bǔ)充。ReNAA對(duì)元素測(cè)量的靈敏度更高,主要分析緩發(fā)伽馬射線,但其需要在反應(yīng)堆內(nèi)進(jìn)行樣品輻照,因此對(duì)樣品性質(zhì)、形態(tài)等要求更高,制樣過程較為復(fù)雜。且樣品在反應(yīng)堆內(nèi)輻照后,劑量較高,需要一段時(shí)間冷卻后才能測(cè)量,測(cè)量周期較長(zhǎng)。此外,反應(yīng)堆內(nèi)中子能譜復(fù)雜,ReNAA容易受到超熱中子、快中子引起的核反應(yīng)(如(n, p), (n, α)反應(yīng))影響。而PGAA,特別是安裝在冷中子束流上的PGAA裝置,幾乎不含超熱中子和快中子。利用中子束流進(jìn)行在線測(cè)量,測(cè)量時(shí)間控制較易實(shí)現(xiàn),且對(duì)樣品形態(tài)要求較ReNAA更低,達(dá)到無損分析的要求,固態(tài)、液態(tài)或氣態(tài)樣品均可進(jìn)行測(cè)量。對(duì)大部分元素如Ti、Al、Fe、Mn、Mg、Ca、Na、K等,利用PGAA和ReNAA進(jìn)行主成分分析時(shí)可以取得良好效果,在進(jìn)行痕量分析時(shí),由于反應(yīng)堆內(nèi)的中子通量較堆外中子束流通量高5~6個(gè)量級(jí),因此ReNAA的效果更好。但對(duì)一些元素,如H、B、Cl、Si、Nb、Sm、Gd等,由于不存在產(chǎn)生緩發(fā)γ射線的反應(yīng)或反應(yīng)截面較小,而瞬發(fā)γ射線的反應(yīng)截面較大,用PGAA進(jìn)行分析更具優(yōu)勢(shì)。通常,對(duì)低Z元素,PGAA的分析效果更好。因此,PGAA和ReNAA結(jié)合,幾乎可以測(cè)量元素周期表上的所有元素,是標(biāo)準(zhǔn)的材料成分分析方法。

    表1 部分基于反應(yīng)堆中子源的PGAA裝置

    除了對(duì)已知元素進(jìn)行定量分析外,PGAA還可根據(jù)測(cè)量到的γ射線能譜,對(duì)材料中存在的未知雜質(zhì)元素進(jìn)行鑒別和含量測(cè)定,被廣泛應(yīng)用于各類研究。

    在化學(xué)研究中,常用于標(biāo)準(zhǔn)物質(zhì)的認(rèn)證。美國(guó)國(guó)家標(biāo)準(zhǔn)技術(shù)研究院(NIST)利用PGAA對(duì)標(biāo)準(zhǔn)物質(zhì)進(jìn)行認(rèn)證工作已經(jīng)有30多年歷史,認(rèn)證了20多種材料,包括H、B、C、N、Na、Al、Si、S、Cl、K、Ca、Ti、Cd、Sm、Gd等核素的含量信息。標(biāo)準(zhǔn)物質(zhì)也推動(dòng)了PGAA技術(shù)的發(fā)展,如日本原子能研究院(JAERI)利用NIST的相關(guān)標(biāo)準(zhǔn)物質(zhì),對(duì)JAERI的PGAA裝置進(jìn)行了測(cè)量精度驗(yàn)證[13],并在此基礎(chǔ)上發(fā)展、豐富了用于PGAA定量分析的k0因子法,對(duì)各元素的k0因子進(jìn)行了測(cè)量標(biāo)定[14]。

    在工業(yè)生產(chǎn)中,PGAA可用于成品中雜質(zhì)含量測(cè)量、參考材料標(biāo)定等研究[15-17],如石墨烯中的硼雜質(zhì)含量分析。PGAA可以實(shí)現(xiàn)多數(shù)同位素豐度的測(cè)定,對(duì)核材料的鑒定有天然優(yōu)勢(shì)。可用于對(duì)核級(jí)石墨或Be組件中各元素(B、Al、Cl、Sc、Ti、Cr、Mn、Fe、Ni、C)的檢測(cè)[11,18-19],可燃毒物組件、硼涂層及用于反應(yīng)性補(bǔ)償?shù)呐鹚崛芤鸿b定等[20]。此外,結(jié)合束流斬波器,PGAA可對(duì)鈾礦石中的235U及238U進(jìn)行測(cè)定[21]。

    在地理、考古學(xué)研究中,PGAA可用于H、B、C、N、Na、Mg、Al、Si、P、S、Cl、K、Ca、Ti、V、Mn、Fe、Co、Ni、Cu、Cd、Nd、Sm、Gd等元素含量的測(cè)定。為地質(zhì)勘探、斷代、隕石與星際物質(zhì)研究、藝術(shù)品鑒定、人類活動(dòng)區(qū)域劃分與生產(chǎn)力發(fā)展階段等領(lǐng)域的研究提供重要的分析方法。

    在食品安全、環(huán)境科學(xué)研究中,PGAA可用于食品原料跟蹤溯源,農(nóng)作物種子甄別等研究。在大氣污染物排放及顆粒物分析治理研究中,PGAA同樣可以發(fā)揮重要作用。在發(fā)達(dá)國(guó)家,利用PGAA對(duì)大氣污染物進(jìn)行研究已有多年的歷史經(jīng)驗(yàn),日本的Yonezawa等分析了日本多地空氣顆粒物中Gd的含量[22],美國(guó)馬里蘭大學(xué)與NIST合作,對(duì)燃煤排放物中的B、S/SO2、Cl、Si、Cd、Gd含量進(jìn)行了分析,并研究了隨大氣擴(kuò)散的分布流程[23]。PGAA結(jié)合ReNAA,可以給出各主要污染物排放點(diǎn)的元素指紋信息,并對(duì)其進(jìn)行溯源,為環(huán)境監(jiān)督及排放監(jiān)管提供重要的技術(shù)支持。這些技術(shù)同樣適用于水體及土壤污染研究,用于測(cè)定元素成分及主要污染物含量。特別是對(duì)近年來國(guó)內(nèi)大力推動(dòng)的食品安全監(jiān)管、污染源追蹤等有重要促進(jìn)作用。

    在生物與醫(yī)療研究中,PGAA可用于分析生物樣品中H、B、N、K、Ca、S、Cl、Cd、Sm、Gd等元素含量,掌握其在各器官、組織中的分布,特別是重金屬污染物在動(dòng)植物,甚至人體器官內(nèi)的沉積情況,為生命科學(xué)研究提供技術(shù)支持。硼中子俘獲治療(BNCT)是治療癌癥的一種有效方法,世界上建有多套基于反應(yīng)堆中子束流的BNCT裝置,如美國(guó)MIT的TRIGA-Ⅱ反應(yīng)堆[24]、捷克LVR-15反應(yīng)堆[25]、荷蘭Petten高通量堆[26]、韓國(guó)HANARO反應(yīng)堆[27]等。BNCT裝置中患者體內(nèi)的硼分布通過PGAA測(cè)量,以實(shí)現(xiàn)對(duì)輻射劑量精確控制,達(dá)到最佳治療效果。

    PGAA在國(guó)防技術(shù)方面的應(yīng)用鮮有報(bào)道,但結(jié)合其技術(shù)特點(diǎn),可以開展多方面的研究,如航空航天領(lǐng)域重要部件元素含量分析、飛機(jī)機(jī)翼與航空發(fā)動(dòng)機(jī)葉片中氫含量檢測(cè),各類火工品、炸藥中關(guān)鍵元素的含量和比例的無損分析,以及受貧鈾彈或生化武器攻擊地區(qū)地質(zhì)元素的分布和水土污染分布情況等。作為一種標(biāo)準(zhǔn)的元素分析方法,PGAA將為我國(guó)材料科學(xué)、環(huán)境科學(xué)、國(guó)防科技等領(lǐng)域提供先進(jìn)的研究方法,促進(jìn)各學(xué)科的發(fā)展。

    傳統(tǒng)的PGAA技術(shù)多用于分析小尺度樣品(樣品尺度通常為mm量級(jí))或淺層分析,當(dāng)樣品需要分析其內(nèi)部元素成分,且尺寸較大時(shí),超出傳統(tǒng)PGAA方法的適用范圍。與中子層析照相技術(shù)結(jié)合,獲取中子束流在樣品內(nèi)部的分布情況,可以解決這個(gè)問題。

    1.3 中子層析照相

    中子層析照相的基本原理是利用圖像重建算法,對(duì)物體的多幅投影圖像進(jìn)行重建得到物體的內(nèi)部結(jié)構(gòu)圖像,其基礎(chǔ)是傅里葉切片定理:極坐標(biāo)形式下的二維傅里葉變換結(jié)果,其極坐標(biāo)為θ直線上的值,與投影角為θ的投影數(shù)據(jù)的一維傅里葉變換結(jié)果相同。給出投影數(shù)據(jù)時(shí)就可以經(jīng)過反傅里葉變換得到物體的內(nèi)部三維圖像。

    傅里葉切片定理原理圖示于圖2,μ(x,y)為樣品內(nèi)部(x,y)處的中子衰減系數(shù),如果采用平行中子束流穿過樣品,沿著其穿透路徑s的總衰減可以用線積分表示:

    (2)

    投影角為θ的投影直線方程,t=xcosθ+ysinθ,式(2)可重寫為:

    ysinθ-t)dxdy

    (3)

    其傅里葉變換為:

    (4)

    樣品的二維傅里葉變換函數(shù)定義如下:

    圖2 傅里葉切片定理原理圖Fig.2 The Fourier-slice theorem

    (5)

    其中,u=wcosθ,v=wsinθ。因此:

    (6)

    旋轉(zhuǎn)樣品,改變?chǔ)冉牵瑴y(cè)量得到多組pθ(t),即可求解獲得μ(x,y)。分析樣品內(nèi)部圖像重建結(jié)果,得到中子在樣品內(nèi)的衰減過程,進(jìn)一步得到中子強(qiáng)度分布結(jié)果,再將其應(yīng)用于PGAA的測(cè)量。常用的重建算法有直接傅里葉變換重建算法、濾波反投影重建算法、迭代重建算法等。

    1.4 PGAA-NT聯(lián)合測(cè)量方法

    PGAA-NT聯(lián)合測(cè)量示意圖示于圖3[28],PGAA裝置與中子層析裝置安裝在同一條中子束流上。測(cè)量過程先進(jìn)行中子層析照相,再進(jìn)行PGAA測(cè)量。測(cè)量時(shí)首先移除中子準(zhǔn)直器,中子束流以面源形式照射樣品,對(duì)樣品進(jìn)行層析照相,開展樣品結(jié)構(gòu)重建,得到樣品內(nèi)部中子束流衰減參數(shù)。層析照相完成后,在樣品前端安裝束流準(zhǔn)直器,根據(jù)需要分析的樣品內(nèi)部具體位置要求,利用中子準(zhǔn)直器與γ準(zhǔn)直孔實(shí)現(xiàn)對(duì)樣品內(nèi)部位置的定位,開展PGAA測(cè)量。

    圖3 PGAA-NT聯(lián)合測(cè)量示意圖Fig.3 Combined PGAA-NT Setup

    兩套裝置安裝在同一位置進(jìn)行測(cè)量,束流有良好的重復(fù)性,且可以利用層析照相裝置對(duì)樣品的定位進(jìn)行校準(zhǔn),提高定位精度。

    中子層析照相得到的樣品內(nèi)部的三維結(jié)構(gòu)以灰度圖的形式給出,為了直接利用其結(jié)果,將重建圖像的灰度值沿中子束流投影路徑的積分與成像屏測(cè)量到的該束流的衰減結(jié)果進(jìn)行擬合,得到圖像灰度值與衰減系數(shù)的對(duì)應(yīng)關(guān)系(圖 4)[28]。

    圖4 重建灰度值積分gΣ與該路徑上中子總衰減系數(shù)μΣ對(duì)應(yīng)關(guān)系Fig.4 Integration of gray-values in reconstruction gΣ vs total neutron attenuation factors μΣ

    通過與中子層析照相技術(shù)的結(jié)合,可以掌握中子束流在樣品內(nèi)部的衰減過程,并得到樣品內(nèi)部的中子束流強(qiáng)度分布情況,進(jìn)一步通過PGAA分析得到樣品內(nèi)部的元素分布,使得可以在不對(duì)樣品外觀進(jìn)行破壞取樣的情況下開展PGAA測(cè)量,可進(jìn)一步發(fā)展為樣品內(nèi)元素分布的三維掃描分析技術(shù)。

    2 國(guó)內(nèi)外研究及應(yīng)用

    2.1 國(guó)外發(fā)展現(xiàn)狀

    作為一種新型的檢測(cè)方法,近年來,德國(guó)FRM II反應(yīng)堆[3]與匈牙利布達(dá)佩斯反應(yīng)堆[4]上建成了將PGAA與中子層析照相聯(lián)合測(cè)量的裝置,可以對(duì)樣品內(nèi)部元素分布進(jìn)行二維/三維掃描、成像。并開展了初步應(yīng)用,F(xiàn)RM Ⅱ上的裝置成功對(duì)隕石樣品、古代文物進(jìn)行了測(cè)量分析,得到了其主要元素的相對(duì)分布情況[28-29]。

    對(duì)墨西哥Allende地區(qū)的一塊隕石碎片進(jìn)行二維元素分布的掃描分析結(jié)果示于圖 5,重點(diǎn)分析了其中硅、鐵、鎂的分布[29]。

    PGAA-NT對(duì)7世紀(jì)生產(chǎn)的未開封容器內(nèi)皮帶扣的元素構(gòu)成進(jìn)行無損分析的結(jié)果示于圖 6,#2、#3、#5、#6位置檢測(cè)到Cl元素的存在,其生產(chǎn)工藝過程可能引入了外部污染[28]。

    圖7是PGAA-NT對(duì)6世紀(jì)的扣針各區(qū)域內(nèi)的元素構(gòu)成進(jìn)行分析,分析了其中的鐵、硫、金、銅、氫、鋁、銀、鎂、氯等元素的分布情況。鐵、金、銅的分布情況分別示于圖8~圖10,并據(jù)此分析了其工藝方法及達(dá)到的水平[28]。

    此研究發(fā)展了近十年,屬于新型的無損檢測(cè)分析方法,目前仍有一些問題需要解決,如中子衰減系數(shù)分布計(jì)算過程中的中子散射、束流硬化、重建噪聲等影響因素的研究,定量分析方法等。

    2.2 國(guó)內(nèi)發(fā)展現(xiàn)狀

    我國(guó)在ReNAA上的研究已有多年經(jīng)驗(yàn),中國(guó)科學(xué)院高能物理研究所、中國(guó)原子能科學(xué)研究院、中國(guó)工程物理研究院等機(jī)構(gòu)開展了多年的堆內(nèi)中子活化分析研究。在PGAA研究方面,國(guó)內(nèi)多家科研機(jī)構(gòu)與高校開展過相關(guān)研究工作,但大多數(shù)是基于同位素中子源或小型加速器中子源。中國(guó)原子能科學(xué)研究院曾在其重水研究堆(HWRR)上建立PGAA實(shí)驗(yàn)裝置[30],并承擔(dān)國(guó)際原子能機(jī)構(gòu)(IAEA)部分?jǐn)?shù)據(jù)的評(píng)價(jià)工作[31],此外,還在硼中子俘獲治療(BNCT)的小型反應(yīng)堆上建立了專用于測(cè)量硼濃度的PGAA裝置,并計(jì)劃在中國(guó)先進(jìn)研究堆CARR上建造基于冷中子束流的PGAA裝置[32-33]。中國(guó)科學(xué)院高能物理研究所柴之芳院士在2014年北川召開的研究堆技術(shù)應(yīng)用學(xué)術(shù)交流會(huì)上,對(duì)中子活化分析的歷史和發(fā)展進(jìn)行了回顧和評(píng)價(jià),并肯定了PGAA技術(shù)的作用與優(yōu)點(diǎn),提出希望在國(guó)內(nèi)開展先進(jìn)的PGAA技術(shù)研究的期望。

    a——隕石;b——隕石的中子成像結(jié)果; c——硅的二維分布;d——鐵的二維分布;e——鎂的二維分布圖5 PGAA-NT對(duì)隕石碎片元素二維分布的測(cè)量分析結(jié)果a——Allende meteorite;b——Radiography of Allende meteorite;c——Si 2D distribution;d——Fe 2D distribution;e——Mg 2D distributionFig.5 2D distribution of the main elements Fe, Si and Mg of a meteorite fragment by PGAA-NT

    圖6 PGAA-NT對(duì)7世紀(jì)容器內(nèi)皮帶扣主要元素構(gòu)成的分析結(jié)果Fig.6 Result of elements analysis of belt mount from 7th century by PGAA-NT

    隨著中國(guó)綿陽研究堆(CMRR)及其國(guó)內(nèi)首個(gè)冷中子源平臺(tái)的投入使用,中國(guó)工程物理研究院核物理與化學(xué)研究所近年來積極開展冷中子束流的相關(guān)應(yīng)用研究,先后建成了多臺(tái)中子散衍射裝置、中子深度分析(neutron depth profiling, NDP)裝置、熱(冷)中子照相及層析照相裝置,層析照相技術(shù)在國(guó)內(nèi)處于先進(jìn)水平。在PGAA研究方面,已經(jīng)成功在同位素與加速器中子源上開展了PGAA研究,即將建成基于CMRR冷中子束流的PGAA裝置,并進(jìn)一步開展PGAA-NT聯(lián)合測(cè)量研究。PGAA裝置安裝在CMRR的一條冷中子束流末端,與中子深度分析裝置和冷中子照相裝置共用同一束流,其總體布置示于圖 11。在PGAA裝置處安裝一套層析照相裝置(圖12),以實(shí)現(xiàn)與PGAA進(jìn)行聯(lián)合測(cè)量。

    3 小結(jié)

    PGAA-NT的結(jié)合,充分利用了NT的技術(shù)特點(diǎn),彌補(bǔ)了PGAA技術(shù)不能深入樣品內(nèi)部進(jìn)行測(cè)量的不足,可以實(shí)現(xiàn)樣品內(nèi)部元素的三維分布測(cè)量分析。國(guó)際上已有近十年的發(fā)展并建成裝置,在考古、文物保護(hù)、星際物質(zhì)分析等方面獲得應(yīng)用,但技術(shù)仍需完善,特別是從中子層析重建到獲取中子束流衰減參數(shù)的過程中,仍需對(duì)各類影響因素進(jìn)行研究,如中子散射、束流硬化、重建噪聲、中子束流發(fā)散情況等。國(guó)內(nèi)機(jī)構(gòu)也啟動(dòng)了相關(guān)研究,需要在研究其基本原理、實(shí)驗(yàn)技術(shù)、數(shù)據(jù)分析方法的基礎(chǔ)上,深入研究、優(yōu)化中子衰減系數(shù)的計(jì)算方法,提高測(cè)量準(zhǔn)確度。該技術(shù)的發(fā)展將提供一種新型的先進(jìn)中子無損檢測(cè)技術(shù),推動(dòng)相關(guān)學(xué)科發(fā)展。

    圖7 PGAA-NT對(duì)6世紀(jì)的扣針內(nèi)各部位元素構(gòu)成進(jìn)行分層分析Fig.7 Analysis of a disc fibula from 6th century by PGAA-NT

    圖8 PGAA-NT用于鐵元素分析的分布情況Fig.8 Iron distribution by PGAA-NT

    圖9 PGAA-NT用于金元素分析的分布情況Fig.9 Gold distribution by PGAA-NT

    圖10 PGAA-NT用于銅元素分析的分布情況Fig.10 Copper distribution by PGAA-NT

    圖11 CMRR冷中子源上裝置總體布置示意圖Fig.11 Sketch of facilities at one cold neutron beam of CMRR

    圖12 PGAA與NT聯(lián)合測(cè)量布置示意圖Fig.12 Sketch of PGAA-NT measurement

    [1] 丁大釗,葉春堂,趙志祥. 中子物理學(xué)[M]. 北京:原子能出版社,2001.

    [2] Molnar G. Handbook of prompt gamma activation analysis[M]. Dordrecht: KLUWER Academic Publishers, 2004.

    [3] Kudejova P, Meierhofer G, Zeitelhack K, et al. The new PGAA and PGAI facility at the research reactor FRM II in Garching near Munich[J]. Journal of Radioanalytical and Nuclear Chemistry, 2008, 278(3): 691-695.

    [4] Kis Z, Szentmiklosi L, Belgya T. NIPS-NORMA station-A combined facility for neutron-based nondestructive element analysis and imaging at the budapest neutron centre[J]. Nuclear Instruments and Methods in Physics Research, 2015, 779: 116-123.

    [5] Lindstrom R M. Prompt-gamma activation analysis[J]. Journal of Research of the National Institute of Standards and Technology, 1993, 98(1): 127-133.

    [6] Isenhour T L, Morrison G H, Chem A. Modulation technique for neutron capture gamma ray measurements in activation analysis[J]. Analytical Chemistry, 1966, 38(2): 162-167.

    [7] Anderson D L, Failey M P, Zoller W H, et al. Facility for non-destructive nnalysis for major and trace elements using neutron-capture gamma-ray spectrometry[J]. Journal of Radioanalytical Chemistry, 1981, 63(1): 97-119.

    [8] Hanna A G, Brugger R M, Glascock M D, et al. Development of a neutron capture gamma-ray facility at the university of Missouri Research Reactor[J]. Transactions of the American Nuclear Society, 1980, 34: 178-179.

    [9] Lindstrom R M, Zeisler R, Rossbach M. Activation analysis opportunities using cold neutron beams[J]. Journal of Radioanalytical and Nuclear Chemistry, 1987, 112(2): 321-330.

    [10]Rossbach M, Chem A. Multielement prompt gamma cold neutron activation analysis of organic matter[J]. Analytical Chemistry, 1991, 63(19): 2 156-2 162.

    [11]Yonezawa C, Wood A K H, Hoshi M, et al. The Characteristics of the prompt gamma-ray analyzing system at the neutron beam guides of JRR-3M[J]. Nuclear Instruments and Methods in Physics Research, 1993, 329(1-2): 207-216.

    [12]Molnár G L, Révay Z, Paul R L, et al. Prompt-gamma activation analysis using thek0approach[J]. Journal of Radioanalytical and Nuclear Chemistry, 1998, 234(1-2): 21-26.

    [13]Yonezawa C. Multi-Element determination by a cold neutron-induced prompt gamma-ray analysis[J]. Analytical Sciences, 1996, 12(4): 605-613.

    [14]Matsue H, Yonezawa C. k0Standardization approach in neutron-induced prompt gamma-ray analysis at JAERI[J]. Journal of Radioanalytical and Nuclear Chemistry, 2000, 245(1): 189-194.

    [15]Kasztovszky Z, Révay Z, Belgya T, et al. Investigation of impurities in thermoluminescent Al2O3materials by prompt-gamma activation analysis[J]. Journal of Analytical Atomic Spectrometry, 1999, 14(4): 593-596.

    [16]Paul R L, Lindstrom R M, Greenberg R R, et al. Determination of hydrogen in titanium alloys by cold neutron prompt gamma activation analysis[J]. Metallurgical and Materials Transactions A, 1996, 27(11): 3 682-3 687.

    [17]Anderson D L. Ceramic glaze analysis by simultaneous in-beam PGAA and XRFS[J]. Journal of Radioanalytical and Nuclear Chemistry, 1995, 192(2): 281-287.

    [18]Jurney E T. Application of the thermal (n,γ) reaction to elemental analysis[M]. New York: Springer US, 1979.

    [19]Yonezawa C, Wood A K H, Chem A. Prompt gamma ray analysis for boron with cold and thermal neutron guided beams[J]. Analytical Chemistry, 1995, 67(24): 4 466-4 470.

    [20]Pecequilo B R S, Suarez A A. Fuel material analysis using radioactive thermal neutron capture[M]. New York: Springer US, 1979.

    [21]Molnár G L, Révay Z, Belgya T. Non-destructive interrogation of uranium using PGAA[J]. Nuclear Instruments and Methods in Physics Research, 2004, 213: 389-393.

    [22]Yonezawa C, Matsue H, Hoshi M. Multielement analysis of environmental samples by cold and thermal guided neutron induced prompt gamma-ray measurement[J]. Journal of Radioanalytical and Nuclear Chemistry, 1997, 215(1): 81-85.

    [23]Kitto M E, Anderson D L, Zoller W H. Simultaneous collection of particles and gases followed by multielement analysis using nuclear Techniques[J]. Journal of Atmospheric Chemistry, 1988, 7(3): 241-259.

    [24]Matsumoto T, Aizawa O. Prompt gamma-ray neutron activation analysis of boron-10 in biological materials[J]. International Journal of Radiation Applications and Instrumentation Part A Applied Radiation and Isotopes, 1990, 41(9): 897-903.

    [25]Honzátko J, Tomandl I. Boron concentration measurement system for the czech BNCT project[C]∥Capture Gamma-Ray Spectroscopy and Related Topics: 10th International Symposium. American Institute of Physics, 2000: 749-750.

    [26]Raaijmakers C P, Konijnenberg M W, Dewit L, et al. Monitoring of blood-10B concentration for boron neutron capture therapy using prompt gamma-ray analysis[J]. Acta Oncologica, 1995, 34(4): 517-523.

    [27]Byun S H, Choi H D. Design features of A prompt gamma neutron activation analysis system at HANARO[J]. Journal of Radioanalytical and Nuclear Chemistry, 2000, 244(2): 413-416.

    [28]Schulze R. Prompt Gamma-ray 3D-imaging for cultural heritage purposes[D]. K?ln: University of K?ln, 2010.

    [29]Canella L, Kudějová P, Schulze R, et al. PGAA, PGAI and NT with cold neutrons: Test measurement on a meteorite sample[J]. Applied Radiation and Isotopes, 2009, 67(12): 2 070-2 074.

    [30]張?zhí)m芝,倪邦發(fā),田偉之,等. 瞬發(fā)γ射線中子活化分析的現(xiàn)狀與發(fā)展[J]. 原子能科學(xué)技術(shù). 2005,39(3):282-288.

    Zhang Lanzhi, Ni Bangfa, Tian Weizhi, et al. Status and development of prompt γ-ray neutron activation analysis[J]. Atomic Energy Science and Technology, 2005, 39(3): 282-288(in Chinese).

    [31]周春梅. 中子俘獲瞬發(fā)γ活化分析數(shù)據(jù)庫[J]. 原子核物理評(píng)論,2003,20(4):295-298.

    Zhou Chunmei. Database of prompt γ-ray activation analysis of neutron capture[J]. Nuclear Physics Review, 2003, 20(4): 295-298(in Chinese).

    [32]Sun H, Ni B, Xiao C, et al. Design of a prompt-gamma neutron activation analysis system on China Advanced Research Reactor[J]. Nuclear Science and Techniques, 2011, 22(5): 287-292.

    [33]孫洪超,倪邦發(fā),肖才錦,等. CARR堆瞬發(fā)伽馬活化分析系統(tǒng)的初步設(shè)計(jì)[J]. 同位素,2012,25(3):182-188.

    Sun Hongchao, Ni Bangfa, Xiao Caijin, et al. The design of PGNAA facility on CARR[J]. Journal of Isotopes, 2012, 25(3): 182-188(in Chinese).

    Combination of Prompt Gamma-ray Activation Analysis and Neutron Tomography

    YANG Xin, LI Run-dong, WANG Guan-bo, TANG Bin, SUN Yong

    (InstituteofNuclearPhysicsandChemistry,ChinaAcademyofEngineeringPhysics,Mianyang621900,China)

    Applications and progress of prompt gamma-ray neutron activation analysis (PGAA/PGNAA) were presented. PGAA was combined with neutron tomography (NT), which was used to reconstruct the 3D construction inside a sample and the neutron beam attenuation inside the sample could be estimated. Further more, neutron beam intensity at a specific position also could be evaluated by NT. The combination of PGAA and NT was able to analyse 3D distributions of elements inside samples without being disassembled. The PGAA-NT facilities were built at FRM Ⅱ and Budapest research reactor, which were applied to study astronomy and ancient samples. The PGAA-NT combination technique is being developed in China.

    prompt gamma-ray activation analysis; neutron tomography; non-destructive testing; image reconstruction

    2017-03-30;

    2017-05-15

    國(guó)家自然科學(xué)基金資助項(xiàng)目(11475152,11405151)

    楊 鑫(1988—),男,貴州人,博士,主要從事反應(yīng)堆物理及應(yīng)用技術(shù)研究

    TL816+.3

    A

    1000-7512(2017)03-0153-11

    10.7538/tws.2017.youxian.019

    猜你喜歡
    中子源束流層析
    “華龍一號(hào)”取消二次中子源研究
    犬細(xì)小病毒量子點(diǎn)免疫層析試紙條的研制
    “超級(jí)顯微鏡”
    少兒科技(2020年3期)2020-05-13 13:15:25
    高能同步輻射光源低能束流輸運(yùn)線設(shè)計(jì)研究
    中國(guó)散裂中子源加速器注入束流損失調(diào)節(jié)研究
    中國(guó)散裂中子源項(xiàng)目冷源系統(tǒng)設(shè)計(jì)
    中國(guó)散裂中子源首次打靶成功獲得中子束流
    ECR中和器束流引出實(shí)驗(yàn)研究
    A族鏈球菌膠體金免疫層析試紙條的制備及應(yīng)用
    新型B族鏈球菌膠體金免疫層析試紙條的臨床應(yīng)用評(píng)價(jià)
    欧美精品av麻豆av| 日韩,欧美,国产一区二区三区| 久久人人爽av亚洲精品天堂| 欧美人与性动交α欧美软件| 老司机深夜福利视频在线观看 | 日韩免费高清中文字幕av| 国产野战对白在线观看| 麻豆av在线久日| 国产无遮挡羞羞视频在线观看| 欧美中文综合在线视频| 男女边摸边吃奶| 成人免费观看视频高清| 巨乳人妻的诱惑在线观看| 水蜜桃什么品种好| 亚洲欧美日韩高清在线视频 | 一级片免费观看大全| 另类亚洲欧美激情| 久久精品人人爽人人爽视色| 一进一出抽搐动态| 免费在线观看日本一区| 蜜桃国产av成人99| 十八禁网站网址无遮挡| 波多野结衣一区麻豆| 久久久精品94久久精品| 国产精品1区2区在线观看. | 侵犯人妻中文字幕一二三四区| 91成人精品电影| 亚洲国产日韩一区二区| 日韩有码中文字幕| 久久 成人 亚洲| 天堂俺去俺来也www色官网| 精品一区二区三卡| 亚洲欧美一区二区三区黑人| 欧美一级毛片孕妇| 五月天丁香电影| 精品亚洲乱码少妇综合久久| 大陆偷拍与自拍| 亚洲一卡2卡3卡4卡5卡精品中文| av在线播放精品| 亚洲三区欧美一区| h视频一区二区三区| 中文字幕人妻丝袜制服| 国产区一区二久久| 欧美乱码精品一区二区三区| 国产亚洲av片在线观看秒播厂| 欧美激情极品国产一区二区三区| 操出白浆在线播放| 中文字幕高清在线视频| 国产日韩欧美亚洲二区| 国产有黄有色有爽视频| 欧美一级毛片孕妇| 九色亚洲精品在线播放| 国产精品国产av在线观看| 69av精品久久久久久 | 女人爽到高潮嗷嗷叫在线视频| 欧美日韩亚洲高清精品| 国产亚洲欧美在线一区二区| 久久这里只有精品19| 看免费av毛片| 亚洲精品国产av蜜桃| 一二三四社区在线视频社区8| 丝袜人妻中文字幕| 久久人人爽av亚洲精品天堂| 在线亚洲精品国产二区图片欧美| 啦啦啦 在线观看视频| 波多野结衣一区麻豆| 91老司机精品| 三上悠亚av全集在线观看| 99国产精品免费福利视频| 免费高清在线观看视频在线观看| 操美女的视频在线观看| 99国产综合亚洲精品| 老司机午夜十八禁免费视频| 男人舔女人的私密视频| 女人被躁到高潮嗷嗷叫费观| 一个人免费在线观看的高清视频 | 国产成+人综合+亚洲专区| 嫩草影视91久久| 国产亚洲欧美在线一区二区| 一本色道久久久久久精品综合| 99久久99久久久精品蜜桃| 亚洲全国av大片| 国产欧美亚洲国产| 亚洲伊人久久精品综合| 免费日韩欧美在线观看| 欧美变态另类bdsm刘玥| 久久性视频一级片| 免费一级毛片在线播放高清视频 | 成人亚洲精品一区在线观看| 中文字幕色久视频| 老司机在亚洲福利影院| 1024香蕉在线观看| 久久精品成人免费网站| 久久久久久久精品精品| 99久久人妻综合| 五月开心婷婷网| 操出白浆在线播放| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久国产精品麻豆| 亚洲人成电影免费在线| 亚洲专区字幕在线| 午夜免费观看性视频| 国产精品久久久久久人妻精品电影 | 日本vs欧美在线观看视频| 午夜久久久在线观看| 嫩草影视91久久| 女性生殖器流出的白浆| 最新在线观看一区二区三区| 手机成人av网站| av欧美777| 国产成人精品久久二区二区91| 中文字幕人妻丝袜一区二区| 日本wwww免费看| 黄色片一级片一级黄色片| 国产成人免费无遮挡视频| 韩国精品一区二区三区| 肉色欧美久久久久久久蜜桃| 日韩中文字幕视频在线看片| a级毛片黄视频| av视频免费观看在线观看| 国产成人影院久久av| 久久综合国产亚洲精品| 亚洲专区中文字幕在线| 久久99一区二区三区| 夫妻午夜视频| 亚洲中文字幕日韩| 久久中文字幕一级| 亚洲av电影在线观看一区二区三区| 欧美 日韩 精品 国产| 久久久国产精品麻豆| 可以免费在线观看a视频的电影网站| 岛国在线观看网站| 久久久久精品人妻al黑| 久久久久视频综合| 久久久精品区二区三区| 人人妻人人澡人人爽人人夜夜| 亚洲免费av在线视频| 我要看黄色一级片免费的| 亚洲国产看品久久| 最黄视频免费看| 精品国产国语对白av| 男女无遮挡免费网站观看| 嫁个100分男人电影在线观看| xxxhd国产人妻xxx| 欧美在线黄色| 成人黄色视频免费在线看| 午夜91福利影院| 日韩视频在线欧美| 免费观看a级毛片全部| 精品高清国产在线一区| 精品一品国产午夜福利视频| 久久精品久久久久久噜噜老黄| 女人爽到高潮嗷嗷叫在线视频| 在线观看免费视频网站a站| 欧美日韩成人在线一区二区| 欧美人与性动交α欧美软件| www.999成人在线观看| 久久这里只有精品19| 在线 av 中文字幕| 欧美日韩中文字幕国产精品一区二区三区 | 肉色欧美久久久久久久蜜桃| 高潮久久久久久久久久久不卡| 韩国精品一区二区三区| 男女边摸边吃奶| 国产av国产精品国产| 国产在线观看jvid| 国产精品偷伦视频观看了| 91精品国产国语对白视频| 一个人免费看片子| 日本av免费视频播放| 精品高清国产在线一区| 久久毛片免费看一区二区三区| 欧美日韩一级在线毛片| 777久久人妻少妇嫩草av网站| 老汉色av国产亚洲站长工具| 欧美另类亚洲清纯唯美| 久久精品国产a三级三级三级| 亚洲一码二码三码区别大吗| 啦啦啦视频在线资源免费观看| 日本欧美视频一区| 91国产中文字幕| 天天躁日日躁夜夜躁夜夜| 国产精品久久久人人做人人爽| 在线观看免费午夜福利视频| 亚洲欧美精品综合一区二区三区| 免费高清在线观看日韩| 免费久久久久久久精品成人欧美视频| 99国产极品粉嫩在线观看| 交换朋友夫妻互换小说| 中文字幕高清在线视频| 亚洲av美国av| 人妻久久中文字幕网| 精品少妇一区二区三区视频日本电影| 免费观看人在逋| av天堂久久9| 午夜免费成人在线视频| 国产黄频视频在线观看| 美女国产高潮福利片在线看| 午夜福利免费观看在线| 日本猛色少妇xxxxx猛交久久| 婷婷成人精品国产| 纵有疾风起免费观看全集完整版| 91精品国产国语对白视频| 久久精品熟女亚洲av麻豆精品| 青草久久国产| 一边摸一边抽搐一进一出视频| 十八禁网站免费在线| 国产97色在线日韩免费| 国产精品久久久久久精品古装| 亚洲九九香蕉| 黄色毛片三级朝国网站| 久久久精品免费免费高清| 国产成人一区二区三区免费视频网站| 日韩中文字幕欧美一区二区| 性色av乱码一区二区三区2| 国产高清videossex| 亚洲 国产 在线| 美女大奶头黄色视频| 啦啦啦在线免费观看视频4| 国产xxxxx性猛交| 又大又爽又粗| 色婷婷久久久亚洲欧美| 久久国产精品大桥未久av| 亚洲国产中文字幕在线视频| 欧美+亚洲+日韩+国产| 99国产精品一区二区蜜桃av | 国产激情久久老熟女| 成年av动漫网址| 亚洲av欧美aⅴ国产| 日日夜夜操网爽| 人妻 亚洲 视频| 日本av免费视频播放| 欧美性长视频在线观看| 大香蕉久久网| 一级毛片电影观看| 午夜影院在线不卡| 操美女的视频在线观看| 亚洲 国产 在线| av超薄肉色丝袜交足视频| 99国产精品99久久久久| 久久av网站| 少妇的丰满在线观看| 国产精品久久久久久精品古装| 亚洲一码二码三码区别大吗| 亚洲欧美清纯卡通| 女性生殖器流出的白浆| 精品少妇久久久久久888优播| 啦啦啦免费观看视频1| 日本一区二区免费在线视频| 亚洲欧美日韩另类电影网站| 亚洲av日韩在线播放| 少妇猛男粗大的猛烈进出视频| 一区二区日韩欧美中文字幕| 国产区一区二久久| 成人国产av品久久久| tube8黄色片| 亚洲天堂av无毛| 精品福利观看| 国产av精品麻豆| 国产精品一区二区免费欧美 | 国产成人精品久久二区二区91| 视频区图区小说| 两个人看的免费小视频| 国产精品99久久99久久久不卡| 午夜福利一区二区在线看| 少妇裸体淫交视频免费看高清 | 亚洲久久久国产精品| 久9热在线精品视频| 国产免费视频播放在线视频| 久久久久久久国产电影| 成人亚洲精品一区在线观看| 成人国产av品久久久| 国产真人三级小视频在线观看| 精品熟女少妇八av免费久了| 热re99久久精品国产66热6| 日本a在线网址| 免费一级毛片在线播放高清视频 | 一区在线观看完整版| 亚洲午夜精品一区,二区,三区| 成在线人永久免费视频| 国产1区2区3区精品| 亚洲第一青青草原| 18禁黄网站禁片午夜丰满| 丰满迷人的少妇在线观看| 香蕉国产在线看| 国产亚洲欧美精品永久| 欧美激情久久久久久爽电影 | 色综合欧美亚洲国产小说| 精品高清国产在线一区| 欧美精品高潮呻吟av久久| 中文字幕人妻熟女乱码| 不卡一级毛片| 亚洲一区二区三区欧美精品| 下体分泌物呈黄色| 天天操日日干夜夜撸| 久久精品国产综合久久久| 国产精品二区激情视频| 免费av中文字幕在线| 日本欧美视频一区| 老司机午夜福利在线观看视频 | 国产av一区二区精品久久| 两人在一起打扑克的视频| 每晚都被弄得嗷嗷叫到高潮| 波多野结衣一区麻豆| 欧美成人午夜精品| 十八禁人妻一区二区| 欧美日韩精品网址| 国产成人精品久久二区二区91| 一区二区三区乱码不卡18| 啦啦啦 在线观看视频| 欧美久久黑人一区二区| 国产91精品成人一区二区三区 | 久热爱精品视频在线9| 久久99一区二区三区| av电影中文网址| 久久久国产精品麻豆| 中文字幕色久视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲一卡2卡3卡4卡5卡精品中文| 久久精品亚洲熟妇少妇任你| 人人妻人人爽人人添夜夜欢视频| 制服人妻中文乱码| 人妻一区二区av| 亚洲欧洲精品一区二区精品久久久| 国产精品香港三级国产av潘金莲| 中文精品一卡2卡3卡4更新| 亚洲欧美日韩无卡精品| 国产成人欧美在线观看| 亚洲自拍偷在线| 三级毛片av免费| 三级国产精品欧美在线观看 | 中文字幕久久专区| 变态另类丝袜制服| 看片在线看免费视频| 老汉色av国产亚洲站长工具| 日日夜夜操网爽| 天堂动漫精品| 窝窝影院91人妻| 长腿黑丝高跟| 脱女人内裤的视频| 琪琪午夜伦伦电影理论片6080| 夜夜躁狠狠躁天天躁| 久久精品国产清高在天天线| 国产人伦9x9x在线观看| 一区福利在线观看| 国内久久婷婷六月综合欲色啪| 中文字幕人成人乱码亚洲影| 欧美中文综合在线视频| 他把我摸到了高潮在线观看| 级片在线观看| 国产成年人精品一区二区| 精品第一国产精品| 国产欧美日韩精品亚洲av| 淫妇啪啪啪对白视频| 亚洲精品色激情综合| 两个人看的免费小视频| 一二三四社区在线视频社区8| 国产精品一区二区免费欧美| 波多野结衣高清作品| 欧美久久黑人一区二区| 国产av不卡久久| 精品国内亚洲2022精品成人| 亚洲精品粉嫩美女一区| 午夜日韩欧美国产| 国产精品精品国产色婷婷| 国产成人精品久久二区二区免费| 久久中文字幕人妻熟女| 女人爽到高潮嗷嗷叫在线视频| av在线天堂中文字幕| 91九色精品人成在线观看| 一个人免费在线观看的高清视频| 国产午夜精品论理片| 我要搜黄色片| 老熟妇乱子伦视频在线观看| a级毛片在线看网站| 国产乱人伦免费视频| 人成视频在线观看免费观看| 小说图片视频综合网站| 久久午夜亚洲精品久久| 18禁观看日本| 亚洲成a人片在线一区二区| 久久欧美精品欧美久久欧美| 久久久国产成人精品二区| 国产精品久久久av美女十八| 欧美日韩亚洲国产一区二区在线观看| 99久久精品热视频| 99热这里只有是精品50| 哪里可以看免费的av片| 国产成人精品久久二区二区免费| 久久人妻福利社区极品人妻图片| 亚洲欧美日韩东京热| 国产成人精品无人区| 又黄又爽又免费观看的视频| 亚洲精品久久国产高清桃花| 亚洲精品中文字幕在线视频| 欧美日韩福利视频一区二区| 2021天堂中文幕一二区在线观| 亚洲成a人片在线一区二区| 村上凉子中文字幕在线| 精品久久蜜臀av无| 黄色女人牲交| 两个人看的免费小视频| 国产三级黄色录像| 久久精品国产综合久久久| 日韩欧美免费精品| 欧美日韩精品网址| 男插女下体视频免费在线播放| 亚洲国产高清在线一区二区三| 后天国语完整版免费观看| 成人av一区二区三区在线看| 国产视频内射| 日韩精品免费视频一区二区三区| 成人18禁高潮啪啪吃奶动态图| 丰满人妻一区二区三区视频av | 国产真实乱freesex| 又粗又爽又猛毛片免费看| 日本黄大片高清| 99久久无色码亚洲精品果冻| 欧美乱码精品一区二区三区| 少妇裸体淫交视频免费看高清 | 香蕉丝袜av| 国产精品亚洲美女久久久| 淫秽高清视频在线观看| 50天的宝宝边吃奶边哭怎么回事| 精品日产1卡2卡| 国产黄片美女视频| 麻豆国产97在线/欧美 | 午夜影院日韩av| 成熟少妇高潮喷水视频| 黄片小视频在线播放| 人妻丰满熟妇av一区二区三区| 国产一区二区激情短视频| 亚洲av成人一区二区三| а√天堂www在线а√下载| 日韩国内少妇激情av| 国产精品 国内视频| 亚洲电影在线观看av| 99在线视频只有这里精品首页| 97碰自拍视频| 中出人妻视频一区二区| 国产蜜桃级精品一区二区三区| 免费电影在线观看免费观看| 成人一区二区视频在线观看| 香蕉丝袜av| 成年版毛片免费区| 后天国语完整版免费观看| 久久欧美精品欧美久久欧美| 校园春色视频在线观看| 男女床上黄色一级片免费看| 黑人操中国人逼视频| 亚洲av第一区精品v没综合| 午夜激情av网站| 午夜精品在线福利| 男人舔奶头视频| 欧美中文日本在线观看视频| 99久久精品热视频| 亚洲精品av麻豆狂野| 免费看美女性在线毛片视频| 三级男女做爰猛烈吃奶摸视频| 青草久久国产| 男女午夜视频在线观看| 亚洲国产中文字幕在线视频| 久久精品影院6| 亚洲av成人一区二区三| 一级片免费观看大全| 亚洲一区二区三区色噜噜| aaaaa片日本免费| 国产一区二区在线观看日韩 | 麻豆国产av国片精品| 嫩草影视91久久| 精品日产1卡2卡| 岛国视频午夜一区免费看| 中文字幕av在线有码专区| 岛国在线免费视频观看| 国产一区二区三区视频了| 97超级碰碰碰精品色视频在线观看| 午夜a级毛片| 夜夜躁狠狠躁天天躁| 听说在线观看完整版免费高清| 亚洲av五月六月丁香网| 别揉我奶头~嗯~啊~动态视频| 亚洲乱码一区二区免费版| 欧美国产日韩亚洲一区| 久久久久久国产a免费观看| 日韩国内少妇激情av| 黄片大片在线免费观看| 国产av在哪里看| 精品久久久久久成人av| 黄色视频不卡| 免费在线观看视频国产中文字幕亚洲| 这个男人来自地球电影免费观看| 男人舔奶头视频| 亚洲avbb在线观看| 日韩三级视频一区二区三区| 国产免费男女视频| 毛片女人毛片| 精品久久久久久,| 久久午夜综合久久蜜桃| 久久久国产精品麻豆| 国产精品国产高清国产av| 国产av又大| 亚洲成人精品中文字幕电影| 亚洲欧洲精品一区二区精品久久久| 琪琪午夜伦伦电影理论片6080| 欧美激情久久久久久爽电影| 99在线视频只有这里精品首页| 亚洲国产精品sss在线观看| 久久人妻av系列| 国产精品精品国产色婷婷| av福利片在线观看| 在线观看免费午夜福利视频| 日韩精品免费视频一区二区三区| 可以在线观看毛片的网站| 日日干狠狠操夜夜爽| 亚洲中文字幕一区二区三区有码在线看 | 又爽又黄无遮挡网站| 国产精品免费一区二区三区在线| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品久久国产高清桃花| 精品人妻1区二区| 伦理电影免费视频| 成在线人永久免费视频| 午夜两性在线视频| 最新在线观看一区二区三区| 美女高潮喷水抽搐中文字幕| 国产熟女xx| 午夜免费激情av| 亚洲乱码一区二区免费版| 最新美女视频免费是黄的| 午夜福利免费观看在线| 日本成人三级电影网站| 国产视频内射| 亚洲男人天堂网一区| 国产免费男女视频| 真人一进一出gif抽搐免费| 搡老岳熟女国产| 免费在线观看日本一区| 亚洲aⅴ乱码一区二区在线播放 | 成人亚洲精品av一区二区| 五月玫瑰六月丁香| 国内精品一区二区在线观看| 99精品久久久久人妻精品| avwww免费| 99国产精品99久久久久| 久久精品91无色码中文字幕| 精品人妻1区二区| 最新美女视频免费是黄的| 免费看美女性在线毛片视频| 欧美日韩黄片免| 成年女人毛片免费观看观看9| 亚洲av成人av| 国内精品一区二区在线观看| 在线观看日韩欧美| 久久亚洲精品不卡| 一级毛片高清免费大全| 久久久久国内视频| 男女午夜视频在线观看| 精品一区二区三区av网在线观看| 免费观看人在逋| 久久婷婷成人综合色麻豆| 欧美极品一区二区三区四区| 视频区欧美日本亚洲| 国产精品99久久99久久久不卡| x7x7x7水蜜桃| 午夜免费观看网址| 国产免费男女视频| 成人av在线播放网站| 久久这里只有精品中国| 国产成人精品久久二区二区免费| 日韩大码丰满熟妇| 国产精品久久久av美女十八| 丁香六月欧美| 在线观看午夜福利视频| 亚洲欧美日韩无卡精品| 成人欧美大片| 亚洲精品久久国产高清桃花| 日韩欧美精品v在线| 久久久精品大字幕| 日日干狠狠操夜夜爽| 91字幕亚洲| 久久午夜综合久久蜜桃| 亚洲欧美一区二区三区黑人| 亚洲专区中文字幕在线| 一本精品99久久精品77| 亚洲av第一区精品v没综合| 国内少妇人妻偷人精品xxx网站 | 久久精品91蜜桃| 老司机午夜福利在线观看视频| 久久香蕉国产精品| 黄色成人免费大全| 99久久精品国产亚洲精品| 久久天堂一区二区三区四区| 国产伦在线观看视频一区| 一级作爱视频免费观看| 久久这里只有精品中国| 亚洲一区二区三区不卡视频| 别揉我奶头~嗯~啊~动态视频| 国产蜜桃级精品一区二区三区| 一本综合久久免费| 欧美日韩一级在线毛片| 久久精品国产清高在天天线| 少妇熟女aⅴ在线视频| 午夜激情福利司机影院| 99精品欧美一区二区三区四区| 白带黄色成豆腐渣| 国产亚洲欧美98| 国产一区二区在线av高清观看| 国产精品亚洲一级av第二区| 91麻豆av在线| 欧美日韩中文字幕国产精品一区二区三区| 亚洲熟女毛片儿| 91麻豆av在线| 中文字幕熟女人妻在线| 人妻夜夜爽99麻豆av| 国产高清videossex|