刘文杰,马庆鹏,王鹏来

(中国工程物理研究院流体物理研究所，四川 绵阳 621999)

爆轰加载下银气溶胶源项实验研究*

刘文杰,马庆鹏,王鹏来

(中国工程物理研究院流体物理研究所，四川 绵阳 621999)

为研究钚气溶胶生成规律及源项分布，提出了一种密封环境下的金属气溶胶源项实验方法，利用爆炸容器及气溶胶采样器等实验设备，开展了钚替代材料银的气溶胶源项生成实验，分析了空气动力学直径小于10 μm的银气溶胶源项分布特性，并与外场扩散试验中钚气溶胶数据进行了对比。进一步采用不同炸药及加载装置，对相同银片样品进行加载实验，研究了加载峰值压力对银气溶胶源项分布的影响。研究结果表明，金属银可作为替代材料研究爆轰加载条件下钚气溶胶的源项参数，在一定实验条件下，二者归一化积累质量分布具有较高的一致性；爆轰加载产生的银气溶胶中，较小粒径气溶胶质量含量较高。气溶胶总量与加载峰值压力呈二次函数分布，当峰值压力超过某个临界值时，气溶胶总量将达到极值；此外，通过对同发次实验的3次不同采样数据进行分析，发现了较明显的气溶胶凝并和沉降现象。

爆炸力学；粒径分布；归一化积累质量；银气溶胶；源项；钚气溶胶

钚-239(下称钚)材料一般通过反应堆中产生的慢中子轰击铀-238人工产生，是武器级核装置的重要裂变材料和核电站的初始加浓原料。作为高活性和高毒性的放射性核素，钚材料具有强烈地诱发骨肿瘤的生物毒性，且半衰期极长，约为24 400年。钚材料的化学性质极为活泼，常见的非合金态α位相是一种室温下坚硬的金属，但仅在小于122 ℃时稳定。当含钚材料的核装置处于高能炸药化学爆炸和燃料燃烧等事故环境中时，钚材料将发生氧化反应并气溶胶化[1-2]。而钚气溶胶扩散至环境后极易因呼吸进入生物肺部和血液循环系统，使生物机体在存活期间遭受持续地内照射[3-4]。R.H.Condit[5]将可能对人体产生吸入危害的钚气溶胶颗粒的空气动力学直径上限定为10 μm，急性致死剂量为4.3 mg。因此，掌握爆轰加载下钚气溶胶的源项参数(产生气溶胶的总质量及粒径分布)，对意外核事故的预防及应急响应措施的制定具有重大的意义。

对于核装置化爆条件下钚气溶胶的研究，仅有在内华达核试验场进行的大型外场安全试验ORC[4,6]。迄今为止，仅有一次系统完备地测量炸药化爆条件下钚气溶胶的扩散试验，ORC为几十年间针对化爆条件下钚气溶胶扩散的数值模拟提供了重要且唯一的源项及扩散分布数据[7-8]。但是，考虑对生物和环境的危害，之后对于钚气溶胶的源项及扩散研究转入核事故、核试验场地的钚气溶胶再悬浮[9-10]、沉积评估[11-12]以及实验室实验、数值模拟。M.J.Sagartz[13]通过PBX-9404炸药加载银球壳实验，模拟核装置在事故条件下钚气溶胶源项生成规律，研究了不同炸药构型条件下银样品气溶胶化的转化率，分析了炸药爆能对银气溶胶转化率的影响规律，但未说明使用金属银作为钚替代材料的依据。伍怀龙等[14]对钚气溶胶环境中放射性气体的迁移及扩散过程开展了实验研究，建立了描述放射性气体迁移过程的物理模型，但对钚气溶胶的源项特征未做明确描述。在数值模拟方面，相关研究机构相继开发了多种针对钚气溶胶扩散研究的数值计算模型，并对核事故扩散场景进行了数值计算。主要有Hotspot模型[15]、DIFOUT模型[16]、MACCS2模型[17]和ERAD模型[18]，前3种是高斯烟云模型，后1种是采用蒙特卡洛方法的三维数值模型。

在模拟核事故条件下钚气溶胶的扩散机制时，所有模型的输入参数均依赖于事故发生的源项数据。然而出于安全性等方面的考虑，直接开展钚材料化爆实验研究缺乏现实可行性。鉴于银是一种低熔点 (961.78 ℃)、高沸点 (2 163 ℃)的金属，除密度相对较低外，其基本物理化学性质与钚具有一定的相似性，且不存在放射毒性等问题。因此，本文中，拟采用银作为钚的替代材料，开展炸药化爆条件下银气溶胶的源项实验，研究粒径10 μm以下 (PM10) 的银气溶胶源项分布特性，探索使用银模拟钚材料气溶胶源项分布的可行性，并进一步利用多种爆轰模型及加载方式，研究加载峰值压力对银气溶胶源项参数的影响特性。

1 实验方法

图1 爆轰加载金属银气溶胶实验装置示意图Fig.1 Illustration of experimental setup

图1为爆轰加载金属银的实验装置示意图，装置由高密封爆炸容器、气溶胶采集系统(多级冲击采样器、压力传感器、热电偶及质量流量计)和抽气系统(真空泵、角阀及波纹管道)组成。其中，爆炸容器为卧式柱形，两端为球头，设计TNT当量为2.5 kg，容器采用高密封设计，600 h总泄漏份额不大于10-8。气溶胶采集系统的核心装置是安德森八级撞击式气溶胶采样器TE-20-800，可对0 ～ 10 μm粒径范围内的气溶胶进行分段采集，图2(a)为该采样器的实物图。该采样器具有一个10 μm粒径的前置切割级，用来切割、分离所有粒径大于等于10 μm的颗粒，如图2(b)所示；中间八级切割粒径分别为9.0、5.8、4.7、3.3、2.1、1.1、0.7和0.4 μm，对应冲击板为0～7，分别采集9～10 μm至0.4～0.7 μm粒径的气溶胶颗粒；后置级用来收集所有粒径小于0.4 μm的气溶胶颗粒。采样器利用惯性原理对气溶胶进行粒径切割，其切割原理如图2(c)所示。当含气溶胶气流由各级加速喷嘴喷出时，受到冲击板的阻挡，粒径小、质量轻的气溶胶将随气流改变流动方向流入下一级，粒径大、质量重的气溶胶将打在冲击板上进而被滤膜收集。利用该原理，每级滤膜收集到的气溶胶粒径分布在该级和上一级截断直径范围内，从而实现了气溶胶多级采样。选用对金属颗粒物吸附性能最优的高密度石英纤维滤膜作为气溶胶样品采集基底。通过压力传感器和热电偶监测容器环境压力及温度；通过真空抽气系统的可调角阀控制抽速，并通过质量流量计监测采样流量，将它限制为(28.3±1.5) L/min，以保证采样器的粒径切割精度。

图2 安德森八级撞击式气溶胶采样器Fig.2 Anderson 8 stage non-viable cascade impactor

图3 平面对碰爆轰加载装置示意图Fig.3 Sketch of detonative device

为了获得较高的加载压力，主要采用平面对碰加载模型(见图3)对银片进行爆轰加载，该构型与ORC-DT钚气溶胶外场扩散试验加载构型[18]类似。单片炸药尺寸为∅100 mm×20 mm，银片放置于两片炸药之间；起爆装置使用BL-21雷管，传爆药为泰安∅32 mm×11 mm；装置外用厚度5 mm的有机玻璃外壳进行约束；实验银样品均为高纯银片(含银99.99%)，质量为45 g。

实验状态参数见表1，共开展7发次爆轰加载条件下的银气溶胶实验。其中前3发次为标定实验，用于评估采样器的气溶胶采样误差；后4发次为变加载压力的银气溶胶源相实验，是为了获得多种加载压力条件下银气溶胶生成总量及粒径分布。加载分为单面和对碰两种方式，单面加载实验装置仍如图3所示，但拆除了底部起爆雷管2和传爆药2。p为各加载条件下作用在银样品上的峰值压力的理论计算结果。其中，平面对碰及单面加载两种加载方式下作用在银片上的峰值压力以炸药的CJ爆轰参数为输入参数，根据Rankine-Hugoniot曲线得到。

表1 实验状态参数Table 1 Experimental details

实验前，先清洗容器内壁，并尽可能减少容器内可能影响气溶胶样品分析的实验辅助材料，如加载装置采用棉绳悬挂而取消木支架，雷管套采用有机玻璃而不用铝材料。

2 结果及分析

2.1 采样数据汇总及误差

炸药爆轰零时后，立即打开容器密封阀和真空抽气系统，利用采样器开始对气相爆轰产物进行分级采集。每次采样时间均为10 min，单次累积抽气量为283 L，约为待测气溶胶体系总量(容器总体积)的3.877 %。采样结束后，回收采样滤膜对沉积在各级滤膜上的气溶胶颗粒物进行物化分析：利用粒度分析仪或电化学消解后能谱分析，给出容器内气相银气溶胶粒径分布范围；利用原子吸收光谱方法，给出气相银气溶胶采集总量。7发次实验的银气溶胶各粒径的采样分析结果如表2所示，其中实验3进行了3次间隔采样。

对气溶胶源项生成的分析主要从两个方面开展：一是各粒径气溶胶质量分布，二是气溶胶的归一化积累质量分布(小于某一粒径的气溶胶质量占采样总质量的百分比)。

图4为实验1～3标定得到的银气溶胶归一化积累质量分布。对比采样1、2可知，在3.3～10 μm粒径范围内，银气溶胶归一化积累质量分布基本吻合；在小于3.3 μm粒径范围内，因小粒径颗粒悬浮运动具有一定随机性且易受气流扰动的影响，3次采样结果稍有差异，但变化趋势一致。综上，3次采样银气溶胶归一化积累质量分布趋势一致，仅在数值上稍有差异。因此，可以认为，采用密封容器开展爆轰加载条件下银气溶胶源项的研究方法是可行的，具有较高的可重复性，其数据是合理、有效的。

表2 银气溶胶采样数据Table 2 Data of collected samples of silver aerosol

图4 爆轰加载下银气溶胶源项生成数据标定Fig.4 Calibration of silver aerosol source-term after explosion

表3 电化学消解银标样试液浓度Table 3 Calibration error of samples

2.2 银气溶胶源项分布及沉降、凝并

图5 银气溶胶源项分布Fig.5 Distribution of silver aerosol source-term

图5为7发次实验爆轰零时后采样银气溶胶粒径质量分布，这7发次实验银气溶胶的粒径质量分布具有一定的相似性。各次采样结果表明，爆轰加载产生的银气溶胶含量峰值均出现在1.1～3.3 μm粒径范围内，在不同的加载方式下，该峰值在1.1～2.1 μm和2.1～3.3 μm两个粒径范围内变化。爆轰零时后，采样1、2、3a、4、5、6、7中，0.7～3.3 μm银气溶胶百分含量分别为66.8%、76.9%、91.0%、93.4%、79.8%、82.2%、60.1%。在整个采样区间内，粒径在0.7～3.3 μm范围内的气溶胶含量占总量比例均高于60%。这表明，在对金属银进行爆轰加载时，产生较小粒径 (0.7～3.3 μm) 气溶胶含量较高，在RHT-901高猛炸药加载下，这个粒径范围内的气溶胶含量甚至可达90%以上。

图6为实验3(RHT-901平面对碰加载)中3次间隔采样银气溶胶归一化积累质量分布及粒径质量分布。由图6(a)可知，爆轰加载零后2 h采样3c的银气溶胶归一化积累质量分布，与ORC-DT外场扩散试验钚气溶胶数据[19-20]具有较高的一致性。这表明，采用银作为替代材料对爆轰加载条件下钚气溶胶源项的模拟是完全可行的。此外，对比3次采样数据(见图6(a))可知，三者在归一化积累质量分布方面有较大的差别，采样时间越早，小粒径气溶胶含量越高。为此，对3次采样得到的粒径质量分布(见图6(b))做进一步地分析。在0～2.1 μm粒径范围内，随着采样时间的推移，该范围内气溶胶质量呈现明显的下降趋势；而在2.1～10 μm粒径范围内，第2次采样得到的气溶胶质量比第1次有明显下降，第3次采样得到的气溶胶质量比第2次有明显上升；在5.8～10 μm粒径范围内，第3次的结果甚至高于第1次采样结果。这表明，随着时间的推移，容器内的银气溶胶发生了较显著的凝并及沉降现象。

图6 实验3气溶胶源项分布Fig.6 Source-term of silver aerosol of experiment 3

综上所述，在采用的实验密封环境实验条件下，采用RHT-901炸药对银金属进行平面对碰爆轰加载，产生的银气溶胶经过2 h稳定后，较接近开放环境化爆条件下钚气溶胶源项分布。同时，爆轰零时后产生的0.7～3.3 μm粒径范围内银气溶胶占PM10银气溶胶比例较高，且随着时间推移，容器内小粒径气溶胶发生较明显的凝并和沉降现象。

2.3 加载峰值压力的影响

图7 不同峰值压力下采样中银气溶胶总量Fig.7 Amount of silver aerosol under different pressures

参照表1～2，得到加载峰值压力与采样中银气溶胶总量的对应关系，并对数据点进行二阶多项式拟合，如图7所示，其中RHT-901平面对碰加载实验1～3的气溶胶总量取平均值。

由图7可知，随着峰值压力增大，银气溶胶总量呈先增大后减小的趋势，表明气溶胶的浓度不是随峰值压力而渐增的。通过对比数值拟合曲线与实验采样结果可知，二者一致性很高，这表明银气溶胶总量与加载峰值压力存在典型的二次函数关系。依据拟合曲线可知，在45.85～83.11 GPa的峰值压力区间内，存在一个银气溶胶产物总量极值，通过计算，得出该极值对应的压力约为68.69 GPa。对于这个结果，文献[21]中提到，在对多个核爆试验场和次临界试验场的扰动再悬浮钚气溶胶浓度测量时得出，核爆试验场的气溶胶浓度明显低于次临界试验场。经对核爆试验场爆心土壤进行筛选、电镜分析后得到大量的微小的琉璃化钚化合物。这表明，当加载峰值压力超过一定值时，被加载金属材料化合物可能会发生琉璃化现象，抑制小粒径气溶胶颗粒的形成，进而导致气溶胶浓度降低。

3 结 论

设计了爆炸容器内金属材料在爆轰加载条件下的气溶胶源项实验方案，通过对银材料爆轰加载获得了PM10银气溶胶的源项参数，实现了对ORC-DT试验中钚气溶胶源项分布的模拟。具体获得了以下结论。

(1) 实验结果表明，在爆炸容器约束空间内，爆轰加载下产生的银气溶胶与钚气溶胶的归一化积累质量分布具有较高的一致性。在一定条件下，可以使用银作为替代材料研究钚气溶胶的源项分布，为后续气溶胶扩散的数值计算提供输入参数。

(2) 在对金属银进行爆轰加载时，较小粒径 (0.7～3.3 μm) 气溶胶含量较高，在RHT-901高猛炸药加载下，这个粒径范围内的气溶胶含量甚至可达总量的90%以上。

(3) 结果表明，爆轰零时后，银气溶胶发生显著的凝并和沉降现象：爆轰零时后1 h内，各粒径范围气溶胶含量均有下降，较小粒径气溶胶含量下降更明显；零后2 h内，较小粒径气溶胶含量进一步减少，但较大粒径气溶胶含量不降反升。

(4) 得到加载峰值压力和生成PM10银气溶胶总量的变化规律，银气溶胶总量与加载峰值压力呈二次函数分布，该现象发生的原因可能在于高加载压力时银材料发生了琉璃化。进一步研究，在该实验装置上，加载峰值压力约为68.69 GPa时，可能获得的PM10银气溶胶质量最大。对该极值，拟在后续的实验工作中予以确认。

[1] Martz J C, Haschke J M. A mechanism for combustive heating and explosive dispersal of plutonium[J]. Journal of Alloys and Compounds, 1998,266(1):90-103.

[2] Haschke J M, Martz J C. Oxidation kinetics of plutonium in air from 500 to 3 500 ℃: Application to source terms for dispersal[J]. Journal of Alloys and Compounds, 1998,266(1):81-89.

[3] 刘文杰,胡八一,李庆忠.核事故条件下钚气溶胶源项研究综述[J].安全与环境学报,2011,11(5):259-263. Liu Wenjie, Hu Bayi, Li Qingzhong. Review of plutonium aerosol source-term in nuclear accident[J]. Journal of Safety and Environment, 2011,11(5):259-263.

[4] Stewart K, Thomas D C, Terry J L, et al. Atomic weapons research establishment[R]. AWRE REPORT NO.0-29/65, 1965.

[5] Condit R H. Plutonium dispersal in fires: Summary of what is known[R]. UCRL-ID-114164, 1993.

[6] Titus R W. Operation roller coaster: Project officers report-project 2.4[R]. POR-2504, 1965.

[7] Boughton B A, DeLaurentis J M. An integral model of plume rise from high explosive detonations[C]∥Davis L R, Catton I. Buoyant plumes: Proceedings of the 24th National Heat Transfer Conference. Pittsburgh, PA, USA: ASME, 1987:27-32.

[8] Boughton B A, DeLaurentis J M, Dunn W E. A note on a stochastic model of turbulent dispersion in the atmosphere[J]. Boundary-Layer Meteorology, 1989,48(4):443-444.

[9] Shinn J H, Homan D N, Hofmann C B. A summary of plutonium aerosol studies: Resuspension at the Nevada Test Site[R]. UCRL-90746, 1986.

[10] Shinn J H. Complementary Pu resuspension study at Palomares, Spain[R]. UCRL-ID-150980, 2002.

[11] Johansen M P, Kamboj S, Kuhne W W. Whole-organism concentration ratios for plutonium in wildlife from past US nuclear research data[J]. Journal of Environmental Radioactivity, 2013,126(6):412-419.

[12] Johansen M P, Child D P, Davis E, et al. Plutonium in wildlife and soils at the Maralinga legacy site: Persistence over decadal time scales[J]. Journal of Environmental Radioactivity, 2014,131(2):72-80.

[13] Sagartz M J. Violent reaction source term study[R]. SAND-943252, 1995.

[14] 伍怀龙,田东风,郝樊华,等.钚气溶胶环境中惰性气体氪迁移过程研究[J].物理学报,2011,60(3):162-165. Wu Huailong, Tian Dongfeng, Hao Fanhua, et al. Investigation on inert gases transport process in plutonium aerosol[J]. Acta Physica Sinica, 2011,60(3):162-165.

[15] Homann S G, Wilson D V. Hotspot health physics code: Version8.0[R]. UCRL-MA-118617, 1995.

[16] Shreve J D. Operation roller coaster: Scientific directors report[R]. DASA-1644, 1965.

[17] Steele C M, Wald T L, Chanin D I. Plutonium explosive dispersal modeling using the MACCS2 computer vode[R]. LA-UR-98-1901, 1998.

[18] Boughton B A, Delaurentis J M. Description and validation of ERAD: An atmospheric dispersion model for high explosive detonations[R]. SAND92-2069-UC-702, 1992.

[19] Beasely R R. Operation roller coaster: Technical photography[R]. POR-2519, 1965.

[20] Church H W. Cloud rise from high-explosives detonation[R]. SC-RR-68-903, 1969.

[21] Shinn J H. Technical basis for air pathway assessment of resuspended radioactive aerosols: LLNL experiences at seven sites around the world[R]. UCRL-JC-115045, 1993.

(责任编辑 丁 峰)

Experimental research of silver aerosol source-term under explosive detonation

Liu Wenjie, Ma Qingpeng, Wang Penglai

(InstituteofFluidPhysics,ChinaAcademyofEngineeringPhysics,Mianyang621999,Sichuan,China)

An experimental method in the confinement vessel for silver aerosol source-term investigation of metal under high explosive detonation was designed to simulate the plutonium aerosol source-term of nuclear devices under the circumstance of high explosive detonation. Seven silver source-term experiments under variable explosive pressure were conducted in the aerosol facility and the size distribution characteristics of the silver aerosol with an aerodynamic diameter less than 10 μm were examined and compared with the data from the field test. The results show that the silver aerosol source-term under specific conditions has good agreement with the plutonium aerosol source-term obtained by the U. S. government under test conditions simulating an open-air detonation accident of nuclear device. The amount of the silver aerosolized in the experiment is determined by the peak value of high explosive pressure that acted on the silver plate and exhibited a quadratic polynomial. The silver aerosol with an aerodynamic diameter ranging from 0.7 to 3.3 μm account for the major part of the silver aerosol induced by high explosive detonation. When the peak value exceeds a certain critical value, the total amount of aerosol reaches the maximum. Settlement and combination of aerosol were observed in the comparative analysis of the size-mass distribution of the samples collected at different times after the explosion.

mechanics of explosion; size distribution; cumulative mass; silver aerosol; source-term; plutonium aerosol

10.11883/1001-1455(2016)06-0774-07

2015-04-09； < class="emphasis_bold">修回日期：2015-05-25

2015-05-25

国家自然科学基金项目(21307119)

刘文杰(1982— )，男，博士，副研究员;

马庆鹏,maqingpeng@caep.cn。

O381 <国标学科代码：1303510 class="emphasis_bold"> 国标学科代码：1303510 文献标志码：A国标学科代码：1303510

A