呼義翔，孫鳳舉，曾江濤，叢培天

（西北核技術(shù)研究所，西安710024；強脈沖輻射環(huán)境模擬與效應(yīng)國家重點實驗室，西安710024）

多級串聯(lián)磁絕緣感應(yīng)電壓疊加器次級阻抗優(yōu)化

呼義翔，孫鳳舉，曾江濤，叢培天

（西北核技術(shù)研究所，西安710024；強脈沖輻射環(huán)境模擬與效應(yīng)國家重點實驗室，西安710024）

基于Thevenin原理，建立了一種多級串聯(lián)磁絕緣感應(yīng)電壓疊加器（magneticallyinsulated induction voltage adder，MIVA）等效電路模型。依據(jù)磁絕緣限定條件，給出了MIVA次級磁絕緣傳輸線（magnetically－insulated transmission line，MITL）最小磁絕緣電流估計方法。仿真分析了12級串聯(lián)MIVA電壓和電流輸出特性，結(jié)合X射線劑量率計算式，給出了X射線劑量率隨次級MITL運行阻抗變化的分布規(guī)律。結(jié)果表明，多級串聯(lián)MIVA輸出X射線劑量率隨次級運行阻抗的增大呈先增大后減小的趨勢，MITL存在一個最佳運行阻抗，可使輸出的X射線劑量率最大?；贛endel和Creedon磁絕緣模型，對比給出了次級MITL幾何阻抗的兩種估計模型，結(jié)果表明，除第1級外，其余各級用兩種模型估計的幾何阻抗偏差系數(shù)均小于2%。

感應(yīng)電壓疊加器；磁絕緣傳輸線；運行阻抗；幾何阻抗；X射線劑量率

磁絕緣感應(yīng)電壓疊加器（magneticallyinsulated induction voltage adder，MIVA）作為一種高電壓、大電流產(chǎn)生系統(tǒng)，是高能閃光X射線照相裝置的核心組成部分［1-3］。為獲得透視對象高品質(zhì)圖像，一般要求閃光X射線源具有高輻射劑量和小焦斑尺寸［4-5］，因此，提高X射線輻射劑量對于發(fā)展MIVA型閃光照相技術(shù)具有重要的應(yīng)用價值。

一般而言，在給定二極管條件下，X射線劑量隨二極管工作電壓和電流的增大而增大［6］。但對于某一確定的MIVA裝置，二極管電壓和電流受限于前級驅(qū)動源輸出功率，且電流與電壓間的分配比例主要取決于MIVA次級磁絕緣傳輸線（magneticallyinsulated transmission line，MITL）的運行阻抗。因此，合理調(diào)整次級MITL運行阻抗，實現(xiàn)X射線劑量最大化，成為值得探索的技術(shù)途徑。

本文以12級感應(yīng)腔串聯(lián)MIVA裝置為研究對象，理論分析了影響該裝置電壓及電流脈沖輸出機理的因素，給出了以最大化X射線輻射劑量為目標(biāo)的次級MITL運行阻抗優(yōu)化方法。

1 多級串聯(lián)MIVA

1.1 多級串聯(lián)MIVA的基本結(jié)構(gòu)

通常，多級串聯(lián)MIVA由幾級至幾十級感應(yīng)腔通過次級磁絕緣傳輸線串聯(lián)組成，以滿足二極管不同驅(qū)動參數(shù)需求，其結(jié)構(gòu)如圖1所示。前級驅(qū)動源包括Marx、主開關(guān)、形成線等，通常采用典型的幾級脈沖壓縮技術(shù)產(chǎn)生感應(yīng)腔饋入電壓。次級MITL由各級感應(yīng)腔內(nèi)壁和一根沿感應(yīng)腔軸線延伸的次級內(nèi)筒構(gòu)成，每一級驅(qū)動源饋入脈沖通過感應(yīng)腔耦合疊加至次級內(nèi)筒，使內(nèi)筒電壓逐級遞增并輸出至負(fù)載區(qū)域。

多級串聯(lián)MIVA裝置的次級MITL電長度同傳輸于其中的電脈沖傳輸時間相當(dāng)，甚至更長，因此設(shè)計MITL時，需考慮其阻抗從源至負(fù)載的過渡方式，以獲得MIVA多級感應(yīng)脈沖的高效疊加與傳輸。

目前，國際上已建的多級串聯(lián)MIVA裝置中［1-3，7-9］，次級MITL阻抗多采用“階梯”變化，即從第1級至負(fù)載端，次級MITL運行阻抗等比例遞增（次級內(nèi)筒直徑逐級遞減），以確保各級感應(yīng)腔饋入次級MITL系統(tǒng)的電壓和電流波形（即功率脈沖）完全一致。同時，負(fù)載與次級內(nèi)筒之間通過匹配阻抗的MITL延伸段連接，在方便照相診斷的同時，可有效緩解二極管工作產(chǎn)物對感應(yīng)腔絕緣堆棧的污染。

1.2 多級串聯(lián)MIVA電路模型

通常要求驅(qū)動MIVA各級感應(yīng)腔的電脈沖以理想IVA時序饋入［10-11］，即下游感應(yīng)腔脈沖饋入時間較緊鄰上游感應(yīng)腔脈沖饋入時間滯后τc，τc為次級電脈沖沿單級感應(yīng)腔軸向的傳輸時間。在此時序下，次級MITL采用“階梯”阻抗變化時，各級感應(yīng)腔饋入次級MITL的電壓與電流波形完全相同，負(fù)載電壓脈沖為各級饋入脈沖的完全線性疊加。針對上述工作特點，從負(fù)載向驅(qū)動源方向，多級串聯(lián)MIVA可視為電壓源電路，基于Thevenin原理，等效電路模型如圖2所示。

圖2中各參數(shù)滿足如下關(guān)系：

式中，Vs為驅(qū)動源等效電壓，MV；n為多級感應(yīng)腔串聯(lián)級數(shù)；VPFL為形成線電壓脈沖峰值，MV；Rs為驅(qū)動源等效內(nèi)阻，Ω；ZPFL為形成線阻抗，Ω；m為驅(qū)動單級感應(yīng)腔的形成線并聯(lián)數(shù)；Zload為負(fù)載等效阻抗，Ω；k為阻尼系數(shù)。

基于該電路模型，在給定負(fù)載阻抗條件下，可估計負(fù)載電壓和電流，通過調(diào)整負(fù)載阻抗大小，實現(xiàn)對X射線輻射劑量的優(yōu)化。優(yōu)化過程中必須考慮給定負(fù)載阻抗參數(shù)條件下，次級MITL最小磁絕緣限制條件，即實際傳輸線電流大于最小磁絕緣電流。

1.3 次級MITL磁絕緣限定條件

基于單純電路分析，圖2中Zload可以是任意的，且當(dāng)Zload與Rs相等時，負(fù)載上可以獲得最大的電功率。但由于MIVA次級為磁絕緣傳輸線，其傳輸電流存在一個最小磁絕緣電流，以滿足磁絕緣條件。針對該電流，Ottinger等人開展的相關(guān)理論研究［12］結(jié)果表明，雖然最小磁絕緣電流與自限制電流并非完全一致，但二者之間的差別非常小，可認(rèn)為二者近似一致。

對于自限制電流，Ottinger等人通過對Mendel模型的粒子仿真修正，給出了滿足工程實踐的計算模型（見式（2）至式（5）），并用于BERTHA程序［1315］。美國海軍實驗室利用該模型順利將Mercury裝置輸出電壓峰值從6MV提升至8 MV［16］，較好地檢驗了模型的有效性。

式中，η和g均為修正系數(shù)；V為MITL傳輸線電壓，MV；IASL為自限制電流，kA；Zop為MITL運行阻抗，Ω。

工程實際中，在給定次級MITL運行阻抗和傳輸線電壓時，利用式（2）至式（5）估計MITL最小磁絕緣電流，通過將該電流與次級MITL實際傳輸電流進行比對，實現(xiàn)對次級MITL工作狀態(tài)的有效判斷。

2 12級串聯(lián)MIVA脈沖輸出特性

針對極端能量密度條件下的材料流體動力學(xué)診斷需求，以12級感應(yīng)腔串聯(lián)MIVA為例，分析輸出電壓和電流脈沖隨次級阻抗的變化。已有研究，如RITS－12計劃［17］認(rèn)為：12級感應(yīng)腔串聯(lián)產(chǎn)生的電壓脈沖，基本上能夠滿足高質(zhì)量面密度材料照相需求。假設(shè)饋入脈沖形成線阻抗ZPFL為6.0Ω，傳輸于其中的電壓脈沖峰值VPFL為1.0MV；每級感應(yīng)腔饋入單條PFL，即m＝1，共12級感應(yīng)腔串聯(lián)，即n＝12；次級MITL運行阻抗采用“階梯”阻抗變化形式，且與負(fù)載Zload匹配，即Zop，i＝i·Zload／n。其中，i為介于1～n之間的整數(shù)；Zop，i為第i級感應(yīng)腔對應(yīng)MITL段的運行阻抗，假設(shè)遠離負(fù)載端的感應(yīng)腔為第1級。

目前，國際上已有的二極管阻抗一般為25～150Ω［18］，本文二極管阻抗在該范圍內(nèi)。基于圖2電路模型，計算獲得如圖3所示的電流和電壓關(guān)系。

可以看出，在次級MITL運行阻抗確定時（以150Ω為例），即運行阻抗同負(fù)載匹配，第1級為12.5Ω，第2級為25.0Ω，第3級為37.5Ω，依次類推，第12級為150.0Ω，傳輸線電流隨負(fù)載電壓的上升而線性減?。▓D中實線所示，定義該曲線為MIVA工作曲線），但其最小磁絕緣電流隨負(fù)載電壓增大呈單調(diào)增長趨勢。因此，依據(jù)磁絕緣限定條件，MIVA運行存在一個最值點，即，傳輸線電流與最小磁絕緣電流相等位置（圖中的A點）。當(dāng)負(fù)載電壓繼續(xù)升高時，即負(fù)載阻抗大于216Ω，傳輸線電流小于最小磁絕緣電流，次級MITL將磁絕緣失效。

當(dāng)次級MITL運行阻抗改變，且始終保持與負(fù)載匹配時，最值點A將沿圖3中實線移動。因此，給定MIVA裝置饋入脈沖參數(shù)下，通過調(diào)整次級MITL運行阻抗，理論上可以實現(xiàn)整個MIVA工作曲線上各電壓、電流參數(shù)的輸出。此外，通過改變驅(qū)動源饋入感應(yīng)腔電壓，也可以使MIVA工作曲線平行移動，實現(xiàn)對負(fù)載電脈沖參數(shù)調(diào)整，如圖4所示。需要指出，改變驅(qū)動源電壓，通常需要調(diào)整裝置多個部件的工作狀態(tài)，包括Marx、主開關(guān)、多針開關(guān)、PFL等。

3 次級MITL運行阻抗優(yōu)化

對于給定MIVA裝置，調(diào)整次級MITL運行阻抗，二極管電壓、電流峰值會發(fā)生改變，但其脈沖前沿和脈寬基本保持不變，因此，裝置輸出X射線脈沖波形基本一致。針對上述特點，以X射線劑量率表征X射線劑量指標(biāo)（扣除了共同的時間積分項），評判MIVA的輻射輸出特性。

已有研究表明，電子束二極管X射線劑量率是其工作電壓、電流的非線性函數(shù)［6］，即

實際上，通過進一步的理論推導(dǎo)和數(shù)據(jù)統(tǒng)計發(fā)現(xiàn)，對于負(fù)載匹配型MIVA裝置（次級MITL與負(fù)載匹配），次級MITL最佳運行阻抗ZOPT是驅(qū)動源電參數(shù)和α的顯性函數(shù)，即

利用式（7）可直接獲得最佳次級MITL運行阻抗，從而為工程實踐提供指導(dǎo)。

4 次級MITL幾何阻抗估計

依據(jù)第3節(jié)分析，已知α條件下，對于負(fù)載匹配型次級MITL，可以獲得一個理論上的最佳運行阻抗值，以實現(xiàn)最大X輻射劑量率。但通過最佳運行阻抗，如何準(zhǔn)確估計MITL幾何阻抗，是工程實踐中必須解決的問題。

國際上關(guān)于運行阻抗與幾何阻抗之間關(guān)系的描述，主要以Mendel和Creedon穩(wěn)態(tài)流模型為主。近年，Ottinger等人基于粒子仿真手段對Mendel模型進行了修正［15］，獲得了更為接近實際的修正模型。本節(jié)將分別介紹修正Mendel理論模型和Creedon理論模型的具體應(yīng)用方法，并對二者計算結(jié)果進行對比分析。

修正Mendel理論模型［15，19］為

式中，ZVAC為次級MITL幾何阻抗，Ω。

Creedon理論模型［20］為

式中，γ0為MITL工作電壓對應(yīng)的相對論因子；VMITL為磁絕緣傳輸線電壓，kV；γl為空間電子層邊緣對應(yīng)的相對論因子；IMITL為磁絕緣傳輸線總電流，kA。

依據(jù)上述理論模型，在已知磁絕緣傳輸線電壓和運行阻抗條件下，可直接獲得幾何阻抗。本文Zload＝120Ω，定義偏差系數(shù)為σ，即

式中，ZVAC－CR為基于Creedon模型計算的MITL幾何阻抗，Ω；ZVAC－ME為基于修正Mendel模型計算的MITL幾何阻抗，Ω。計算的12級串聯(lián)次級MITL幾何阻抗，如表1所列。

由表1可以看出，除第1級外（偏差系數(shù)約7.5%），其余各級偏差系數(shù)均小于2%，兩模型獲得了較好的相互印證?？傮w而言，修正Mendel模型經(jīng)過粒子模型校正，考慮了非穩(wěn)態(tài)因素造成的影響（體現(xiàn)在修正系數(shù)η和g上），計算結(jié)果相對更為準(zhǔn)確，特別是在傳輸線電壓小于7MV的條件下，計算結(jié)果更為合理。

5 小結(jié)

通過調(diào)整多級串聯(lián)MIVA裝置次級MITL運行阻抗，可以在負(fù)載上獲得一系列峰值連續(xù)變化的電壓、電流脈沖。在輻射轉(zhuǎn)換系數(shù)α確定的條件下，存在一個最佳次級MITL運行阻抗，使得負(fù)載匹配型MIVA裝置具有最大X射線劑量率輸出。最佳運行阻抗可表示為輻射轉(zhuǎn)換系數(shù)α、感應(yīng)腔串聯(lián)級數(shù)n、單級感應(yīng)腔PFL并聯(lián)數(shù)m和PFL等效阻抗ZPFL的顯性函數(shù)。工程實踐中，基于修正Mendel模型，可估計各級感應(yīng)腔最佳運行阻抗對應(yīng)的MITL幾何阻抗。

需要指出，限于目前對二極管物理機制認(rèn)識不足和二極管工作的復(fù)雜性，輻射轉(zhuǎn)換系數(shù)α尚難以用解析函數(shù)描述，只能通過對輻射場數(shù)值分析或?qū)嶒灲Y(jié)果數(shù)據(jù)擬合給出。因此，工程設(shè)計中，需要將次級MITL最優(yōu)阻抗同二極管輻射場數(shù)值計算結(jié)果進行相互迭代才能得到。

［1］RAMIREZ J J，PRESTWICH K R，JOHNSON D L，et al.Performance of the Hermes－III gamma ray simulator［C］／／Proc of the 7th IEEE International Pulsed Power Conference，Monterey，1989：26 31.

［2］BRUNER N，MOSTROM C，ROSE D V，et al.Modeling the RITS－6transmission line［C］／／Proc of the 16th IEEE International Pulsed Power Conference，Albuquerque，2007：807 810.

［3］THOMAS K，BEECH P，BROWN S，et al.Status of the AWE hydrus IVA fabrication［C］／／Proc of the 18th IEEE International Pulsed Power Conference，Chicago，2011：1 042 1 047.

［4］LISOWSKI P W，PAISNER J A.An advanced hydrotest facility for the stockpile stewardship program［C］／／Proc of the 4th International Meeting on Nuclear Apply of Accelerator Technology，2000：13 25.

［5］EKDAHL C.Modern electron accelerators for radiography［C］／／Proc of the 13th IEEE International Pulsed Power Conference，Las Vegas，2001：21 29.

［6］OLIVER B V，HAHN K，JOHNSTON M D，et al.Recent advances in radiographic X－ray source development at Sandia［C］／／Proc of the 17th International Conference on High Power Particle Beams，Xi'an，2008：1 5.

［7］CORCORAN P A，ALTES R，DOUGLAS J，et al.Design of a driver for the Cygnus X－ray source［C］／／Proc of the 13th IEEE International Pulsed Power Conference，Las Vegas，2001：591 595.

［8］OLIVER B V，GENONI T C.Two and three－dimensional MITL power－flow studies on RITS［C］／／Proc of the 14th IEEE International Pulsed Power Conference，Dallas，2003：395 398.

［9］COMMISSO R J，ALLEN R J，BOLLER J R，et al.Status of the Mercury pulsed power generator，a 6－MV，360kA，magnetically－insulated inductive voltage adder［C］／／Proc of the 14th IEEE International Pulsed Power Conference，Dallas，2003：383 386.

［10］PATE R C，PATTERSON J C，DOWDICAN M C，et al.Self－magnetically insulated transmission line（MITL）system design for the 20－stage Hermes－III accelerator［C］／／Proc of 6th IEEE International Pulsed Power Conference，Arlington，1987：478 481.

［11］LONG F W，MCKEE G R，STOLTZFUS B S.Shaping the output pulse of a linear－transformer driver module［J］.Phys Rev ST Accel Beams，2009，12（3）：030402.

［12］OTTINGER P，SCHUMER J W.Self－consistent modeling of power flow in a recyclable transmission line for Z pinch driven IFE system［J］.IEEE Trans Plasma Sci，2007，35（2）：154 164.

［13］OTTINGER P F，SCHUMER J W，ALLEN R J，et al.Modeling magnetically insulated power flow in Mercury［C］／／Proc of 14th International Pulsed Power Conference，Dallas，2003：849 852.

［14］ALLEN R J，OTTINGER P F，SCHUMER J W，et al.Circuit simulation of the MITL in an IVA with a non－ideal center conductor［C］／／Proc of 16th International Pulsed Power Conference，Albuquerque，2007：1 095 1 101.

［15］OTTINGER P F，SCHUMER J W.Rescaling of equilibrium magnetically insulated flow theory based on results from particle－in－cell simulations［J］.Phys Plasma，2006，13（6）：063109.

［16］ALLEN R J，COMMISSO R J，COOPERSTEIN G，et al.Extension of the operating point of the Mercury IVA from 6to 8MV［C］／／Proc of the 18th IEEE International Pulsed Power Conference，Chicago，2011：1 036 1 041.

［17］SMITH I D，BAILEY V L，F(xiàn)OCKLER J，et al.Design of a radiographic integrated test stand（RITS）based on a voltage adder，to drive a diode immersed in a high magnetic field［J］.IEEE Trans Plasma Sci，2000，28（5）：1 653 1 659.

［18］SMITH I D.Induction voltage adders and the induction accelerator family［J］.Phys Rev ST Accel Beams，2004，7（6）：064801.

［19］MENDEL C W.Planar one－dimensional magnetically insulated electron flow for arbitrary canonical－momentum distribution［J］.J Appl Phys，1979，50（6）：3 830 3 837.

［20］CREEDON J.Magnetic cutoff in high－current diodes［J］.J Appl Phys，1977，48（3）：1 070 1 077.

Optimization of the Secondary Impedance of a Multi－Stage Magnetically－Insulated Induction Voltage Adder

HU Yi－xiang，SUN Feng－ju，ZENG Jiang－tao，CONG Pei－tian
（Northwest Institute of Nuclear Technology，Xi'an 710024，China；State Key Laboratory of Intense Pulsed Radiation Simulation and Effect，Xi'an 710024，China）

Based on the Thevenin theorem，an equivalent circuit model was developed for a multi－stage magnetically－insulated induction voltage adder（MIVA）.According to the limitation of the magnetic insulation theory，a calculation method of the minimum magneticallyinsulated current was described for the magnetically－insulated transmission line（MITL）.The characteristics of the output voltages and currents were analyzed for a 12－stage MIVA.Also，by using the dose estimation equation，the output X－ray dose rate was quantified as a function of the operation impedance of the MITL.Simulation results show that the output X－ray dose rate increases firstly and then decreases with the increase of the MITL operation impedance.There would be an optimized operation impedance which will lead to the maximum of the X－ray dose.Finally，based on the Mendel and Creedon model，two methods were presented for estimations of the MITL geometrical impedance.The calculation resultsindicate that，except for the first stage，deviations in the MITL geometrical impedances obtained by these two models for the other stages are less than 2.0%.

magnetically－insulated induction voltage adders；magnetically－insulated transmission line；operation impedance；geometrical impedance；X－ray radiation dose rate

TM501.2；TM564.1

A

2095 6223（2015）03 191 06

2015 04 04；

2015 07 15

呼義翔（1984－），男，陜西榆林人，工程師，博士，主要從事高功率脈沖能量傳輸與匯聚技術(shù)研究。

E－mail：huyixiang＠nint.ac.cn