李進(jìn)璽，邱孟通，程引會(huì)，吳 偉，來(lái)定國(guó)，馬 良，趙 墨，郭景海（西北核技術(shù)研究所強(qiáng)脈沖輻射環(huán)境模擬與效應(yīng)國(guó)家重點(diǎn)實(shí)驗(yàn)室，陜西西安 710024）

四路并聯(lián)二極管輻射X射線場(chǎng)參數(shù)計(jì)算

李進(jìn)璽，邱孟通，程引會(huì)，吳 偉，來(lái)定國(guó)，馬 良，趙 墨，郭景海
（西北核技術(shù)研究所強(qiáng)脈沖輻射環(huán)境模擬與效應(yīng)國(guó)家重點(diǎn)實(shí)驗(yàn)室，陜西西安 710024）

摘要：利用二極管的電壓、電流計(jì)算了發(fā)射電子束能譜參數(shù)，建立了四路并聯(lián)二極管陽(yáng)極靶蒙特卡羅粒子輸運(yùn)計(jì)算模型，給出了輻射X射線場(chǎng)參數(shù)；將四路并聯(lián)二極管的每個(gè)二極管劃分為若干小單元，將其作為點(diǎn)源，采用數(shù)值積分的方法計(jì)算了輻射X射線劑量分布，并分析了空間不同位置處每路二極管對(duì)劑量的貢獻(xiàn)。結(jié)果表明：真空中，距離四路并聯(lián)二極管陽(yáng)極靶5cm位置處，X射線注量為3．55mJ／cm2，光子平均能量為62．18keV，120keV以下的光子占輻射X射線譜總能量的81．84%，電子束轉(zhuǎn)換效率為0．30%；在2 700cm2范圍內(nèi)，中軸線和對(duì)角線上的劑量均勻性分別為3．20和6．31；在2 000cm2范圍內(nèi)，中軸線和對(duì)角線上的劑量均勻性均小于2。

關(guān)鍵詞：并聯(lián)二極管；X射線場(chǎng)；蒙特卡羅方法；劑量均勻性

脈沖X射線環(huán)境，特別是大面積、高能注量脈沖X射線，在電離輻照效應(yīng)測(cè)試方面有重要的應(yīng)用價(jià)值［1－4］。如何獲得滿足要求的X射線輻射場(chǎng)環(huán)境一直是脈沖功率技術(shù)領(lǐng)域研究的問(wèn)題［5－8］。大面積多環(huán)二極管技術(shù)［6－7］是產(chǎn)生強(qiáng)脈沖硬X射線的技術(shù)途徑之一，其主要用于產(chǎn)生輸出窗面積大且輻射均勻的硬X射線。大面積多環(huán)二極管的陰極由多個(gè)直徑不等的圓環(huán)構(gòu)成，陽(yáng)極為平面輻射靶。多路脈沖功率源各自驅(qū)動(dòng)輻射二極管，各路二極管輻射場(chǎng)在空間疊加形成大面積的輻射X射線場(chǎng)，同時(shí)，通過(guò)降低脈沖功率源的電壓可降低輻射X射線的能量。

本文利用二極管電子束能譜參數(shù)，采用蒙特卡羅方法對(duì)四路并聯(lián)二極管輻射產(chǎn)生的X射線輻射場(chǎng)參數(shù)進(jìn)行數(shù)值模擬，并采用點(diǎn)源數(shù)值積分的方法計(jì)算輻射X射線劑量分布，分析每路二極管對(duì)劑量的貢獻(xiàn)，給出劑量隨空間距離的變化曲線。

1　計(jì)算模型

四路并聯(lián)二極管計(jì)算模型結(jié)構(gòu)如圖1所示。每個(gè)陰極環(huán)外徑200mm，環(huán)內(nèi)、外徑相差10mm，相鄰兩個(gè)環(huán)外徑之間的距離為30mm。四路二極管電子束能譜參數(shù)相同，單路二極管脈沖功率源輸出電壓、電流如圖2所示。

圖1　 四路并聯(lián)二極管計(jì)算模型結(jié)構(gòu)Fig．1 Structure of calculation model for four parallel diodes

利用二極管的電壓、電流在時(shí)間上的對(duì)應(yīng)關(guān)系可換算得到對(duì)應(yīng)的電子能譜［9］，圖3示出了圖2中的二極管電壓、電流所對(duì)應(yīng)的電子能譜，四路并聯(lián)電子束總能量為3．2kJ。

圖2　 單路脈沖功率源輸出電壓、電流Fig．2 Current and voltage of single pulse power supply

圖3　 四路并聯(lián)二極管電子能譜Fig．3 Electron energy spectrum of four parallel diodes

2　計(jì)算結(jié)果

2．1 X射線能譜

輻射靶結(jié)構(gòu)如圖4所示，陽(yáng)極轉(zhuǎn)換靶由20μm的鉭和3mm的聚乙烯構(gòu)成。

圖4　 輻射靶結(jié)構(gòu)Fig．4 Structure of radiation targets

鉭、聚乙烯以及真空中距陽(yáng)極靶5cm處截面的輻射X射線能譜和前向電子能譜如圖5所示。表1列出了不同截面處的輻射X射線和前向電子參數(shù)的比較。

圖5　 不同截面處輻射X射線（a）和前向電子（b）能譜Fig．5 Energy spectra of X－ray（a）and forward electron（b）at different sections

表1　 不同截面處輻射X射線和前向電子參數(shù)的比較Table 1 Comparison of parameters of X－ray and forward electron at different sections

由圖5及表1可看出：

1）輻射X射線能譜由軔致輻射和特征X射線兩部分組成；

2）只有鉭靶時(shí)，輻射場(chǎng)中的電子數(shù)和電子總能量的份額均比X射線的大得多；

3）聚乙烯靶對(duì)鉭靶產(chǎn)生的輻射X射線的能譜分布的影響很小，只會(huì)衰減部分低能光子，但可明顯衰減鉭靶的前向電子；

4）真空中距陽(yáng)極靶5cm截面處，輻射X射線轉(zhuǎn)換效率為0．30%，平均光子能量為62．18keV，在邊長(zhǎng)為52cm的正方形平面上，X射線注量為3．55mJ／cm2，光子與電子數(shù)之比為1．96×104，120keV以下光子占X射線譜總能量的81．84%。

2．2 劑量分布

設(shè)S0為點(diǎn)源強(qiáng)度，a為點(diǎn)源到測(cè)量點(diǎn)的距離，φ為測(cè)量點(diǎn)的強(qiáng)度，不考慮衰減時(shí)，有［10］：

將圖1中的每個(gè)二極管劃分為若干小單元，將其作為點(diǎn)源，利用上式采用數(shù)值積分的方法計(jì)算了輻射X射線劑量分布，并分析了每路二極管對(duì)空間不同距離處劑量的貢獻(xiàn)。計(jì)算中，考慮了點(diǎn)源在不同距離平面上投影角度不同。

4個(gè)二極管順序以及對(duì)角線和中軸線的定義如圖1所示，測(cè)量點(diǎn)在邊長(zhǎng)為52cm的正方形平面上。對(duì)于中軸線上的劑量，環(huán)1和環(huán)2、環(huán)3和環(huán)4是對(duì)稱的，其貢獻(xiàn)也相同，因此，只給出環(huán)1和環(huán)3的貢獻(xiàn)；對(duì)于對(duì)角線上的劑量，環(huán)1和環(huán)4是對(duì)稱的，其貢獻(xiàn)也相同，因此，只給出環(huán)1、環(huán)2和環(huán)3的貢獻(xiàn)。假設(shè)源強(qiáng)度為1，圖6～8分別示出了距陽(yáng)極靶1、7、30cm位置處的劑量分布。

由圖6～8可看出：

1）環(huán)1、環(huán)3對(duì)中軸線上劑量的貢獻(xiàn)沿中軸線對(duì)稱，4個(gè)環(huán)對(duì)中軸線上劑量的貢獻(xiàn)相同。

圖6　 距離陽(yáng)極靶1cm處劑量分布Fig．6 Dose distribution from anode target 1cm

圖7　 距離陽(yáng)極靶7cm處劑量分布Fig．7 Dose distribution from anode target 7cm

2）環(huán)1對(duì)對(duì)角線上劑量的貢獻(xiàn)沿對(duì)角線方向?qū)ΨQ，環(huán)2、環(huán)3對(duì)對(duì)角線上劑量的貢獻(xiàn)沿對(duì)角線方向?qū)ΨQ；距環(huán)較近時(shí)，環(huán)2、環(huán)3對(duì)對(duì)角線上劑量的貢獻(xiàn)遠(yuǎn)大于環(huán)1、環(huán)4的貢獻(xiàn)，隨著距離的增大，4個(gè)環(huán)對(duì)對(duì)角線上劑量的貢獻(xiàn)趨于一致。

3）距環(huán)較近時(shí)，中軸線上劑量呈馬鞍形分布，而對(duì)角線上劑量呈多個(gè)馬鞍形分布，對(duì)角線上劑量大于中軸線上劑量，隨著距離的增大，兩個(gè)方向均逐漸過(guò)渡為余弦分布，劑量大小也趨于一致。

圖9示出了不同位置中軸線以及對(duì)角線上歸一化劑量的比較，圖10示出了總的歸一化劑量與空間距離D的關(guān)系曲線。

圖8　 距離陽(yáng)極靶30cm處劑量分布Fig．8 Dose distribution from anode target 30cm

圖9　 不同位置處歸一化劑量Fig．9 Normalized dose for different locations

圖10　 總的歸一化劑量與距離的關(guān)系曲線Fig．10 Total normalized dosefor different distances

由圖9可看出：

1）中軸線上劑量均勻性好于對(duì)角線。

2）5cm位置處，在2 700cm2范圍內(nèi)，中軸線和對(duì)角線上的劑量均勻性分別為3．20、6．31；在2 000cm2范圍內(nèi)，中軸線和對(duì)角線上的劑量均勻性均小于2。

3）20cm位置處，在2 700cm2范圍內(nèi)，中軸線和對(duì)角線上的劑量均勻性分別為1．73、2．60；在2 000cm2范圍內(nèi)，中軸線和對(duì)角線上的劑量均勻性分別為1．41、1．90。

由圖10可看出，歸一化劑量隨著空間距離的增大呈指數(shù)衰減，可用以下擬合公式計(jì)算不同距離處的歸一化劑量：

其中：A＝0．85；t＝12．79cm；y0＝0．17；D為空間距離，cm。

表2列出了強(qiáng)光一號(hào)［5］、閃光二號(hào)［11］以及四路并聯(lián)二極管產(chǎn)生的硬X射線參數(shù)的比較。表中，閃光二號(hào)和四路并聯(lián)二極管的光子能量為平均光子能量，強(qiáng)光一號(hào)的為光子能譜范圍；光子份額指120keV以下光子占X射線總能量的份額；強(qiáng)光一號(hào)和閃光二號(hào)的數(shù)據(jù)為測(cè)量結(jié)果，四路并聯(lián)二極管的數(shù)據(jù)為計(jì)算結(jié)果。

表2　 不同加速裝置X射線參數(shù)的比較Table 2 Comparison of X－ray parameters of different accelerators

3　結(jié)論

利用二極管的電壓、電流計(jì)算得到了電子的能譜分布；建立了四路并聯(lián)二極管蒙特卡羅粒子輸運(yùn)計(jì)算模型，計(jì)算得到了X射線輻射場(chǎng)參數(shù)；采用點(diǎn)源數(shù)值積分的方法計(jì)算了輻射X射線劑量分布特性。給出的二極管電壓電流和陽(yáng)極靶厚度條件下，距離陽(yáng)極靶5cm截面處，輻射X射線轉(zhuǎn)換效率為0．30%，平均光子能量為62．18keV，在邊長(zhǎng)為52cm的正方形平面上，X射線注量為3．55mJ／cm2，光子與電子數(shù)之比為1．96×104，120keV以下光子占X射線譜總能量的81．84%，在2 000cm2范圍內(nèi)，中軸線和對(duì)角線上的劑量均勻性均小于2。

模擬結(jié)果表明：四路并聯(lián)二極管產(chǎn)生的脈沖硬X射線環(huán)境滿足電離輻照效應(yīng)實(shí)驗(yàn)要求［1－4，11－12］，能應(yīng)用于系統(tǒng)電磁脈沖等電離損傷效應(yīng)研究中。

參考文獻(xiàn)：

［1］ 李進(jìn)璽，程引會(huì)，吳偉，等．系統(tǒng)電磁脈沖綜合環(huán)境中導(dǎo)線響應(yīng)的計(jì)算［J］．強(qiáng)激光與粒子束，2010，22（7）：1 615－1 618．LI Jinxi，CHENG Yinhui，WU Wei，et al．Calculation of wires responses caused by SGEMP comprehensive environments［J］．High Power Laser and Particle Beams，2010，22（7）：1 615－1 618（in Chinese）．

［2］ 李進(jìn)璽，程引會(huì)，李寶忠，等．不同能量光子輻照電纜響應(yīng)規(guī)律研究［J］．原子能科學(xué)技術(shù)，2011，45（7）：893－896．LI Jinxi，CHENG Yinhui，LI Baozhong，et al．Responses laws of coaxial cable to different X－ray energy［J］．Atomic Energy Science and Technology，2011，45（7）：893－896（in Chinese）．

［3］ 林鵬，王肖鈞．強(qiáng)脈沖X射線輻射熱力學(xué)效應(yīng)的數(shù)值模擬［J］．中國(guó)科學(xué)技術(shù)大學(xué)學(xué)報(bào)，2007，37（7）：732－737．LIN Peng，WANG Xiaojun．Numerical simulation of thermal－mechanical effects by intense pulsed X－ray irradiation［J］．Journal of University of Science and Technology of China，2007，37（7）：732－737（in Chinese）．

［4］ 郭紅霞，陳雨生，張義門(mén)，等．穩(wěn)態(tài)、瞬態(tài)X射線

輻照引起的互補(bǔ)性金屬－氧化物－半導(dǎo)體器件劑量增強(qiáng)效應(yīng)研究［J］．物理學(xué)報(bào)，2001，50（12）：2 279－2 283．GUO Hongxia，CHEN Yusheng，ZHANG Yimen，et al．Study of relative dose－enhancement effects on CMOS device irradiated by steady－state and transient pulsed X－rays［J］．Acta Physica Sinica，2001，50（12）：2 279－2 283（in Chinese）．

［5］ 邱愛(ài)慈．脈沖X射線模擬源技術(shù)的發(fā)展［J］．中國(guó)工程科學(xué)，2000，2（9）：24－28．QIU Aici．The development of technology for pulsed X－ray simulators［J］．Engineering Science，2000，2（9）：24－28（in Chinese）．

［6］ 蒯斌，邱愛(ài)慈，王亮平，等．強(qiáng)脈沖超硬X射線產(chǎn)生技術(shù)研究［J］．強(qiáng)激光與粒子束，2005，17（11）：1 739－1 743．KUAI Bin，QIU Aici，WANG Liangping，et al．Generation of intense pulsed super hard X－ray ［J］．High Power Laser and Particle Beams，2005，17（11）：1 739－1 743（in Chinese）．

［7］ BLOOMQUIST D D，STINNETT R W，Mc－DANIEL D H，et al．Saturn，a large area X－ray simulation accelerator［C］∥Proc of 6th IEEE Power Pulsed Conference．US：IEEE，1987：310－317．

［8］ 樊亞軍，石磊，邱愛(ài)慈．相對(duì)論電子束在角向磁場(chǎng)中產(chǎn)生硬X射線的實(shí)驗(yàn)研究［J］．強(qiáng)激光與粒子束，2002，14（6）：929－932．FAN Yajun，SHI Lei，QIU Aici．Experimental study on hard X－ray generation of relativistic electron beams in azimuthal magnetic field［J］．High Power Laser and Particle Beams，2002，14（6）：929－932（in Chinese）．

［9］ 全林，張永民，屠荊，等．強(qiáng)流二極管產(chǎn)生脈沖X

射線能譜場(chǎng)穩(wěn)定性研究［J］．強(qiáng)激光與粒子束，2006，18（2）：317－320．QUAN Lin，ZHANG Yongmin，TU Jing，et al．Stability of pulse X－ray spectrum field generated by intense diod［J］．High Power Laser and Particle Beams，2006，18（2）：317－320（in Chinese）．

［10］馮江平，陳羽，孫慧斌，等．屏蔽材料對(duì)γ射線屏蔽情況的蒙特卡羅模擬［J］．核電子學(xué)與探測(cè)技術(shù)，2011，31（1）：106－112．FENG Jiangping，CHENG Yu，SUN Huibin，et al．Study on theγ－ray shielding effect of shielding material with MCNP－4C［J］．Nuclear Electronics ＆Detection Technology，2011，31（1）：106－112（in Chinese）．

［11］來(lái)定國(guó)，張永民，李進(jìn)璽，等．強(qiáng)流電子束軔致輻射復(fù)合薄靶設(shè)計(jì)［J］．強(qiáng)激光與粒子束，2013，25（6）：1 396－1 400．LAI Dingguo，ZHANG Yongmin，LI Jinxi，et al．Design of bremstrahlung composite thin converter for high current electron beams［J］．High Power Laser and Particle Beams，2013，25（6）：1 396－1 400（in Chinese）．

［12］周輝，李寶忠，王立君，等．不同注量X射線系統(tǒng)電磁脈沖響應(yīng)的數(shù)值計(jì)算［J］．計(jì)算物理，1999，16（2）：157－161．ZHOU Hui，LI Baozhong，WANG Lijun，et al．The calculation of SGEMP response in various ranges of X－ray fluence［J］．Chinese Journal of Computational Physic，1999，16（2）：157－161（in Chinese）．

Calculation of X－ray Field Parameters Generated by Four Parallel Diodes

LI Jin－xi，QIU Meng－tong，CHENG Yin－h(huán)ui，WU Wei，LAI Ding－guo，MA Liang，ZHAO Mo，GUO Jing－h(huán)ai
（National Key Laboratory of Intense Pulsed Radiation Simulation and Effect，Northwest Institute of Nuclear Technology，Xi’an710024，China）

Abstract：The electron energy distribution of diode was calculated using current and voltage．Moreover，the characteristics of X－ray spectrum in parallel diodes environment were calculated using Monte Carlo simulation model．The dose distribution of X－ray of four parallel diodes environment was calculated using numerical integral method．In the model，the small cells，which were divided by diodes，were used as radiation source．The results indicate that the X－ray fluence in the position from the target 5cm is 3．55mJ／cm2．The energy of photons with less than 120keV in the X－ray spectrum is less than 81．84%of the total energy．The average photon energy of the X－ray spectrum is 62．18keV，and the diode electron beams conversion efficiency is 0．30%．In the area of 2 700cm2，the dose uniformities of axial and diagonal orientation are 3．20and 6．31 separately，and in the area of 2 000cm2，the values are less than 2．

Key words：parallel diode；X－ray field；Monte Carlo method；dose uniformity

作者簡(jiǎn)介：李進(jìn)璽（1978—），男，甘肅白銀人，副研究員，碩士，從事輻射效應(yīng)研究

收稿日期：2014－04－10；修回日期：2014－09－25

doi：10．7538／yzk．2015．49．08．1460

文章編號(hào)：1000－6931（2015）08－1460－07

文獻(xiàn)標(biāo)志碼：A

中圖分類號(hào)：TL501