低空核爆炸環境效應模擬研究

李曉菲 李帆 尹祿高 王翔 錢玙 周晨

（1 北京宇航系統工程研究所，北京 100076; 2 湖北大學, 武漢 430062; 3 武漢大學, 武漢 430072）

0 引言

爆炸是指在極短時間內由高能物質（如炸藥等）釋放出大量能量，并產生高溫高壓，對周圍介質的物理狀態產生急劇的變化[1]。核爆炸時，彈體內釋放出巨大的能量，因此核武器是一種大規模的殺傷武器。從核爆炸的破壞作用而言，核爆炸具有極大的直接和間接破壞作用，直接破壞作用主要有沖擊波和光輻射等，又稱為瞬時殺傷因素，而間接破壞作用主要有核輻射和放射性沾染[2-5]。直接破壞作用的時間和空間作用范圍相對于間接破壞作用小，間接破壞作用中的核輻射可能導致大范圍的空氣附加電離，這就可能導致局部的無線電信號干擾[3-6]。

核武器是利用能自持進行的核裂變或聚變反應，瞬間釋放的巨大能量并產生爆炸作用，是具有大規模殺傷破壞效應的武器的總稱[7]。核武器由于其巨大的威力和破壞性，已成為國際軍事戰略平衡和戰略反制的核心籌碼。因此，越來越多的國家都在努力發展先進的核武器技術，以維護其國家的戰略利益。

核爆炸產生的沖擊波和光輻射都是在爆后幾秒到幾十秒時間內起殺傷破壞作用，又稱之為瞬時殺傷破壞因素[2,5]，其作用范圍與當量有關，但最大不會超過二十公里范圍[5]。低空核爆炸中，空氣沖擊波是主要殺傷破壞因素，能在大范圍內殺傷人員、破壞物體；光沖量是重要殺傷破壞因素，占爆炸能量的30%～40%。低空核爆炸所釋放出來的瞬發核輻射和緩發核輻射均能引起大氣電離；形成距火球幾百米以及達到電離層D層的電離，具體的影響范圍與爆炸高度和爆炸當量有關。核爆炸引起的附加電離持續時間也是很長的，對于低空核爆炸，僅考慮對電波傳播具有的影響的時間，也可以持續數小時[4-6,8-10]。由于開展核試驗受到限制，因此利用理論分析和數值模擬、結合已公開的試驗數據研究核爆炸效應，進而分析核爆炸發生時產生的破壞作用，可以為核爆炸效應研究提新的研究方法。本文采用理論分析和數值模擬方法研究了低空核爆炸的沖擊波、光輻射和附加電離等效應的時空分布。

1 低空核爆炸沖擊波數值模擬

核爆炸屬于高溫、高壓、瞬時的物理過程[11]，核爆炸瞬間會形成高溫火球，其猛烈向外膨脹、壓縮周圍空氣形成的高壓氣浪就是沖擊波。高壓氣浪壓縮區的前界面是沖擊波陣面，其壓力最大，沖擊波陣面后的壓力逐漸變小。沖擊波以超音速向四周傳播，隨著傳播距離的增大，傳播速度變慢，壓力逐漸減小最后變為聲波。

沖擊波是導彈、運載火箭、衛星等航天器在總裝、運輸、發射與飛行等全壽命周期階段中需要經受的嚴酷的力學環境之一[12]。沖擊波約占核爆炸能量的一半，是核爆炸的主要殺傷破壞因素[2-5]。沖擊波陣面與它前面的未擾動空氣相比，是一個具有很陡的壓強和很高的密度、質點速度以及溫度變化的突變峰面[10]。核爆炸沖擊波與一般爆炸沖擊波在大氣中的傳播過程是類似的，是大振幅擾動的非線性波在氣體中的傳播過程。它的基本特征是：波的傳播速度是超聲速的，并且波內氣體受到顯著的壓縮。因此在理論上可用可壓縮流體力學來描述。實踐也證明流體力學能夠定量地反映沖擊波的傳播規律。由于爆炸力學所考慮的運動變化快，以致流線之間不會發生明顯的動量和能量的遷移，所以可以忽略粘性和熱傳導。

沖擊波遇到地面會發生反射，從而形成反射沖擊波，又稱為反射波。距爆心投影點一倍爆高距離范圍內的區域稱為規則反射區，其形成過程如下：當球形沖擊波傳到地面時，與地面相撞，便產生了反射沖擊波，反射沖擊波也是以球形向空間傳播。由于反射波是在入射波壓縮、加熱后的空氣中傳播，傳播的速度比入射波快，它能很快追上入射波，并形成合成波，這個位置距爆心投影點的距離大約等于一倍爆高，如圖1所示。圖中t1、t2、t3、t4是沖擊波傳播的四個不同時刻，實線表示入射波，虛線表示反射波。α0為入射角，α2為反射角。在規則反射區內，入射波與反射波的交點始終在地面上。

圖1 沖擊波在地面上的規則反射[10]Fig.1 The regular reflection of shock waves from the ground [10]

對于核爆炸沖擊波的數值研究可以采用以點爆炸的自模擬解為初始條件，在大氣環境中求解二維歐拉方程組得到沖擊波流場的演化過程[13-15]。本文設置模擬參數如下：爆炸當量為20kt TNT，爆高3 km，風速0 m/s。模擬結果如圖2～圖7所示，分別是沖擊波波陣面爆后500 ms、1500 ms、5000 ms、7000 ms、10000 ms、11000 ms的壓力分布圖，各圖中橫軸是地面水平距離，單位是km；縱軸是距離地面的高度，單位是km；圖中顏色代表了壓力值，單位是×105Pa。

圖2 沖擊波波陣面壓強在爆后500 ms的空間分布Fig.2 Spatial distribution of the pressure on the front of the shock wave at 500ms

從圖中可以看出，低空核爆炸產生的沖擊波隨著時間不斷向外擴張，沖擊波的壓強在爆炸瞬間達到極值，然后隨時間而逐漸減小，之后快速衰減接近正常水平。在圖5中，此次模擬中沖擊波波陣面約在7000ms之前接觸地面，之后出現了反射現象，形成了反射波，反射波也是以球形向空間傳播，呈現出饅頭狀構型，這與文獻[5]中描述的空爆外觀現象是一致的。如圖6-圖7所示，沖擊波發生了規則反射，入射波與反射波的交點始終在地面上。

圖3 沖擊波波陣面壓強在爆后1500ms的空間分布Fig.3 Spatial distribution of the pressure on the front of the shock wave at 1500ms

圖4 沖擊波波陣面壓強在爆后5000ms的空間分布Fig.4 Spatial distribution of the pressure on the front of the shock wave at 5000ms

圖5 沖擊波波陣面壓強在爆后7000ms的空間分布 Fig.5 Spatial distribution of the pressure on the front of the shock wave at 7000ms

圖6 沖擊波波陣面壓強在爆后10000 ms的空間分布Fig.6 Spatial distribution of the pressure on the front of the shock wave at 10000ms

圖7 沖擊波波陣面壓強在爆后11000 ms的空間分布Fig.7 Spatial distribution of the pressure on the front of the shock wave at 11000ms

2 低空核爆炸光輻射數值模擬

大氣層核爆炸火球迅速向外膨脹，壓縮周圍的空氣形成沖擊波，同時向外發出光輻射。光輻射是核爆炸的主要殺傷因素之一，主要由核爆炸火球產生。核爆炸火球的發展過程可粗略的分為三個階段，第一階段是輻射擴展階段，第二階段是沖擊波擴張階段，第三階段是復燃冷卻階段[10]?；鹎虬l展過程中將出現兩個極大亮度和一個極小亮度，在此過程中輻射擴張階段和沖擊波擴張階段形成了核爆炸的第一亮度極大，極小亮度時刻沖擊波脫離火球，第三階段出現了第二亮度極大，各個階段時間和爆炸當量有關。 火球亮度第一極大時間為

火球最小亮度時間為

火球亮度第二極大時間為

式中，時間單位為毫秒。

對于爆高在30 km以下的低空核爆炸，在核爆炸發生瞬間，出現極強烈的閃光，隨即出現一個圓而明亮，高溫高壓的火球，這個火球逐漸膨脹并向外擴展，但亮度逐漸下降，直至最小值。之后，隨著沖擊波離開火球越遠，波后空氣對輻射的吸收越來越弱，逐漸變為透明，里面的火球一方面逐漸冷卻，另一方面卻越來越看的清楚，從外部觀察來看好像火球在此重新燃燒，這稱之為“復燃”；在復燃過程中，火球亮度回升，逐漸達到亮度極值。

本文通過對核爆炸光沖量的計算來模擬核爆炸光輻射的發展過程，核爆炸光沖量與環境因素和爆炸當量[14]有很強的相關性。光沖量定義為：在火球整個發光過程中，投射到垂直于光輻射傳播方向的物體表面單位面積的總能量，也稱為曝輻射量。以4.19J/cm2為單位的光沖量公式如下

式中Q為核爆炸當量，單位為kt，R為距離，單位為km，J為大氣透射比。討論這個公式關鍵在于給出J。大氣透射比J的公式為

其中JR(λ)是氣體分子透射比，Jas(λ)是氣溶膠粒子透射比，JN(λ)是氣體污染物吸收透射比，Ja(λ)是水蒸氣吸收透射比，Jo(λ)是臭氧吸收透射比。 氣體分子透射比公式

其中h1是爆高，h2是觀察點的高度，α是觀察點至爆高的天頂角，H0是大氣消光系數從爆高起向上衰減的標高

其中取a=1.057，b=0.437，V為水平能見度。μR(λ,0)是波長為λ時的氣體分子消光系數（單位km-1）

ρ0是地面空氣密度（單位mg/cm3）。

氣溶膠粒子透射比公式

其中Ha是大氣消光系數從爆高起向上衰減的標高，取為1.25 km，μas(λ,0)Ha是波長為λ時的氣溶膠粒子消光系數（單位km-1）。

氣體污染物吸收一般只考慮二氧化氮對光沖量輻射傳輸的影響，氣體污染物吸收透射比公式

其中，μN是氣體污染物吸收系數（單位km-1）

[NO2]是大氣中二氧化氮濃度（單位10-6g/m3）。水蒸氣吸收透射比公式

其中A、B為與吸收帶有關的常數，取A= 0.0279，B= 0.871；Sω為觀察點高度的比濕，定義為

是同體積下水蒸氣質量與空氣質量的比值，Sω0是地面的比濕；p1是觀察點高度的大氣壓，p0是地面的大氣壓。

臭氧吸收透射比公式

設置模擬參數如下：爆炸當量為20kt TNT，波長為0.57 μm，爆高為3 km，觀察點高度為10.5 km，二氧化氮濃度為36×10-3μg/m3，地面比濕為0.2，觀察點比濕為0.15，水平能見度為15 km時，光沖量如圖8所示，各階段光輻射計算結果如圖9-12所示，其中圖9是起爆時刻光輻射，圖10是第一最大亮度時光輻射，圖11是第一最小亮度時刻光輻射，圖12是第二最大亮度時刻光輻射。

圖8 光沖量結果Fig.8 Light impulse results

圖9 起爆時刻光輻射Fig.9 Light radiation at initiation time

圖10 第一最大亮度時刻光輻射Fig.10 Light radiation at the first maximum brightness moment

圖11 第一最小亮度時刻光輻射Fig.11 Light radiation at the first minimum brightness moment

圖12 第二最大亮度時刻光輻射Fig.12 Light radiation at the second maximum brightness moment

圖8 與文獻[5]中的計算結果趨勢、量級一致。通過圖9-圖12可以看出，模擬結果顯示的光沖量變化過程與試驗觀測結果一致，即爆炸瞬間出現了一個圓而明亮的火球，隨后火球逐漸膨脹且亮度增加，達到第一亮度極大；之后火球亮度下降，達到極小值；在此只有火球復燃，且達到第二亮度極大。

3 低空核爆炸附加電離效應模擬

核爆炸釋放的大部分能量都能在空間大氣層中不同高度上產生不同程度的電離，其中一部分能量是以微粒物質（中子、β粒子、核彈頭以及其附近的高速碎片）形式釋放出來的；另一部分能量是以光（核武器的熱輻射和γ射線）形式釋放出來。其中90%～95%的能量是瞬間（微秒級）釋放的，約5%～10%的能量是以裂變的放射性衰變形式延遲輻射的。

核爆炸所釋放出來的輻射射線按照事件發生的先后可分為兩類[4-5]：一類是瞬發核輻射（如X射線、中子流、γ射線等），是在核爆炸發生的瞬間由核武器所釋放出來的，占據了核爆炸的大部分能量，因此可以在大氣層中形成很高的電離并且釋放這些射線的源點位于核爆炸爆心的位置。在這些瞬發核輻射中X射線占50%以上，而γ射線、中子僅占核武器當量的千分之一到百分之一之間。但因為X射線的自由程很小，當核爆炸發生在60km以下時，只在爆心附近形成高溫等離子體，對D區影響不大。對于瞬發γ射線和中子，雖然它們所占核武器總當量的比例很小，但它們在大氣層中的自由程卻遠遠大于X射線。有研究表明，只要核爆炸發生在20km以上，γ射線和中子就可以作用到電離層D區。因此，相對于X射線，當核爆炸發生在60km以下時，瞬發核輻射在電離層D區形成的附加電離區主要由瞬發γ射線、瞬發中子流來決定[9]。另一類是緩發核輻射（如γ射線、β粒子等），是核輻射的裂變產物在很長一段時間內逐漸釋放出來的，緩發核輻射在大氣中形成附加電離的機制，基本上與瞬發核輻射相同；雖然這類輻射僅占核爆炸釋放能量的很少一部分，但是由于它的持續時間長，因此會對大氣層造成長時間的電離效應[16]。

核爆炸所釋放出來的各種輻射射線均含有很高的能量。當這些射線在媒質中傳播時，會與媒質中的分子、原子發生相互作用，把它們所攜帶的能量傳遞給媒質中的分子、原子。當傳遞的能量足夠大時，會使得原子中的某個電子脫離原子力影響而成為自由電子，從而造成媒質的電離。對于空氣中的核爆炸，核輻射射線是在大氣中傳播，因而會導致大氣分子的電離，在大氣中形成電離區。根據輻射傳輸理論，做如下假設[8]：（1）射線在傳播過程中無能量損失；（2）核輻射由某一點源釋放；（3）大氣中所沉淀的射線能量全部用來造成大氣電離。采用參考文獻[9]中的數值方法，模擬參數選取如下：爆高3km，爆炸當量為20kt，爆炸地點地理坐標（30°N，114°E），爆炸時間為2000年3月21日10:00LT，模擬結果如圖13-15所示，其中圖13是瞬發核輻射產生的附加電離的時空分布，圖14是緩發γ射線產生的附加電離的時空分布，圖15是緩發β粒子產生的附加電離的時空分布。各圖中橫軸為距離爆心的水平距離，單位是km；縱軸是距離地面的高度，單位是km；色條代表了核輻射產生的電子密度的常用對數值，單位是：10n個/立方厘米。

圖13 瞬發核輻射產生的附加電離時空分布Fig.13 Spatial and temporal distribution of additional ionization from transient nuclear radiation

圖14 緩發γ射線產生的附加電離時空分布Fig.14 Spatiotemporal distribution of additional ionization produced by delayed γ-rays

圖15 緩發β粒子產生的附加電離時空分布Fig.15 Spatiotemporal distribution of additional ionization produced by delayed β particles

從圖15中可以看出，在低空核爆炸中，由于釋放源的位置位于稠密大氣區，在爆炸初期的輻射受到稠密大氣的阻礙，大部分輻射射線和粒子不能到達電子、正離子有效復合率較低的上層大氣而形成附加電離。它們只能在爆心周圍形成附加電離區。從模擬結果來看，球內電離區的電子密度分布最初具有球狀形態，幾秒鐘后隨著密度的擴散發展成為圓環形，并逐漸形成均勻分布。

對于爆高3km的低空核爆炸而言，由于自由電子非常迅速地附著到較稠密大氣中的中性分子上，因而火球內電離區存在的時間是非常短暫的，其持續時間約為1分鐘左右。大部分能量沉積在鄰近爆心的大氣中，形成火球和空氣沖擊波?；鹎蛑械碾娮用芏冗_到了1012cm-3以上，這比天然電離層F2層峰的電子密度（106cm-3）要高得多。緩發核輻射產生的附加電離主要出現在火球外，發生的高度集中在30-60km高度上，但是產生的電子密度值較低，與自然電離層D層密度相當，但是持續時間不長，約10分鐘后附加電離的電子密度值與自然電離層D層電子密度值相比更低。

4 結論

空中核爆炸是自然空間中的一個巨大的能量源，對空間環境造成了嚴重的影響。從模擬結果來看，低空核爆炸沖擊波的影響范圍比較大，傳播速度很快，超壓變化很快；光輻射變化與以往試驗現象一致；低空核爆炸的附加電離在爆炸瞬間量級很高，超過了自然電離層的強度，但受到周圍稠密大氣的影響，附加電離的持續時間和影響范圍都有限，但仍然能在爆后極短時間內嚴重影響通信信道。

本文完成了低空核爆炸產生的沖擊波、光輻射和附加電離效應的時空演化的數值模擬，實現了低空核爆炸多物理場的數值模擬。