引言
年8月4日傍晚，黎巴嫩首都贝鲁特港口区发生剧烈爆炸[1]，巨响震天升起蘑菇云。爆炸释放了巨大的能量，造成数千人伤亡，数十万人无家可归。图1所示的视频是从住宅楼拍摄到的爆炸场面，该视频稍后将用来估计本次爆炸所释放的能量。
图1:从住宅楼拍摄到的贝鲁特港口爆炸事故爆炸发生瞬间一度被认为是地震、核打击。贝鲁特市长更是把爆炸比作年日本广岛爆炸的原子弹。
图2:从海面拍摄到的贝鲁特港口爆炸事故
黎巴嫩总理证实爆炸是6年前被扣押并存放在港口的吨硝酸铵引发的。硝酸铵（
）是一种铵盐，因含氮量高，被常用于制造高氮化肥。硝酸铵遇热会分解，在不同温度下，硝酸铵会发生不同的化学反应。并且温度越高，反应越迅速。在
时：
在-
时：
在
以上，伴有弱光：
在
以上，发生爆炸：
从以上反应方程式可以看出，硝酸铵分解会产生大量的气体。一千克硝酸铵爆炸会产生升气体，这些气体会在短时间内释放，产生巨大的压力和气浪，对周围的建筑物造成破坏。因此，硝酸铵也常被作为采矿、采石和土木建筑中使用的爆炸混合物成分之一。纯硝酸铵在常温下是稳定的，对打击、碰撞或摩擦均不敏感。但如果储存时间过长，并且通风不佳，其自身会发生反应产生热量，进而导致温度升高并且燃烧，引发爆炸。
引发贝鲁特港口爆炸事故的原因正在调查，不过也可以估个大概。若非人为，则是由于硝酸铵发生热分解，存料下层集聚的热量未能及时导出，发生自燃并引起爆炸。港口环境复杂，附近很可能还存放有其它易爆物。加之近期贝鲁特气温高达30多度，更加为爆炸提供了有利条件。
本次爆炸究竟产生了多大的能量？通常我们用释放相同能量的TNT炸药的质量表示爆炸产生的能量，写成TNT当量。一吨TNT完全爆炸释放的能量等于4.
J。而一吨硝酸铵完全爆炸释放的能量等于0.42吨TNT，因此贝鲁特港口爆炸当量的理论上限是0.42
=吨。如果硝酸铵的有效爆炸率为25%，则吨硝酸铵爆炸能量相当于吨TNT。吨TNT是什么概念？可以参考一下原子弹的当量，年8月9日美国在日本长崎投掷的原子弹的TNT当量为1.3万吨。但本次事故中，我们并不清楚硝酸铵的有效爆炸率到底是多少。
历史最有名的爆炸当量估计是对原子弹爆炸的估计。年7月16日美国科学家在新墨西哥州阿拉莫戈多沙漠进行了三位一体核试验（Trinity），但是试爆后原子弹的某些属性却并不明朗[2]。
图3:原子弹爆炸产生的巨大蘑菇云
负责研制核装置的项目组顾问杰弗里·泰勒是一名剑桥大学教授，他亲眼目睹了这场核实验，但是仅靠肉眼并不能估测出原子弹的威力。年之后，三位一体核试验原子弹爆炸的照片逐渐被媒体公开，泰勒根据爆炸后的一系列照片建立数学模型，成功估算出原子弹爆炸时所释放的能量。本文就来展示一下泰勒估计原子弹当量的模型[3]，并应用该模型对本次贝鲁特港口爆炸进行分析。
模型
假设点源爆炸球形冲击波的半径（
）仅仅依赖于爆炸后的时间（
）、爆炸瞬间所释放的能量（
）和空气的密度（
）。即
其中
是无量纲的常数。利用量纲齐次原则（即等式两边的量纲要相同），可以求出
、
、
的值。我们将变量替换为各自的基本量纲，例如将半径
换成它的量纲[
]，本文所涉及到的基本量纲包括质量（
），长度（
）和时间（
），因此变量对应的基本量纲为：
，
，
，
。将变量各自的量纲代入上式，整理后可得到：
由于等式两边的量纲相同，因此
联立以上三式即可得到
=2/5，
=1/5，
=-1/5，所以
其中
是一个无量纲的常数，是热容比（
，或者称作绝热指数）的低阶函数，因此在任何条件下都近似为1[3]。可以推出爆炸所释放的能量
对上式两端同时除以
并取对数可以得到
或
结果三位一体原子弹
根据以上模型得到的公式，我们可以根据原子弹爆炸后某时刻的冲击波半径来估算原子弹爆炸释放的能量。图4是原子弹爆炸后
毫秒后产生的冲击波照片。照片已经根据原比例尺缩放到等比例的坐标系中。通过计算可以得出冲击波的半径大约为m。空气的密度为1.2
，代入模型公式，可以得到原子弹爆炸产生的能量
约为
J，这一能量相当于1.7万吨TNT。这个简单的计算得到的结果同官方公开的数据（两万吨TNT）相差仅15%。
图4:原子弹爆炸25毫秒时刻的冲击波[3]
此外，我们还可以对
和
进行线性拟合。拟合需要原子弹爆炸不同时刻的冲击波半径，这些数据可以在文献[4]中找到（见附录）。拟合出的直线如图5所示。
图5:原子弹爆炸后时间和冲击波半径的拟合
结果显示
与
呈现非常好的线性关系，这也再一次验证了模型的正确性。拟合还得到了直线在纵坐标轴上的截距：
代入空气的密度，可得原子弹爆炸释放的能量为
这相当于1.8万吨TNT。通过拟合得到的结果同官方公开的数据（二万吨TNT）相差仅10%。
贝鲁特港硝酸铵
为了应用前文的模型计算贝鲁特港硝酸铵的爆炸当量，需要知道爆炸后时间与冲击波半径的数据。我们对图1视频中爆炸瞬间的六帧图像进行了提取，并对图像中的爆炸冲击波半径进行了估计，具体如图6，图中的长度单位是像素。
图6:硝酸铵爆炸瞬间的冲击波半径
我们还需要一个能将像素转为真实长度的标尺。最接近爆炸中心的建筑物是贝鲁特港口的谷仓（左侧白色建筑物）。通过Googleearth可以估计出该谷仓的长度
=
30米，并根据该谷仓炸毁前的照片可估计出该谷仓的长高比
=2.6，从而最终估算出该谷仓高约
=/2.6=50米（图7）。
图7:爆炸前的谷仓长度和长宽比估计
该谷仓在图6a中约占据75个像素，因此爆炸点附近的一个像素表示的实际长度约为50/75=2/3米。由此可估算出图6中各帧图像的冲击波半径，例如在第2帧图像（图6b）中，冲击波的半径约为个像素，对应的实际长度约为
2/3=80米。
图6中相邻图像的时间间隔可由视频帧率来计算。该视频是由手机所拍，其帧率为30帧/秒，因此相邻帧的时间间隔为
秒。以第一帧图像（图6a）为0时刻，我们得到了一组不同时刻
的冲击波半径
。假设第一帧图像对应着爆炸发生后
（
）秒，并按照下式对爆炸后时间
和冲击波半径
进行拟合：
拟合结果如图8所示。
图8:硝酸铵爆炸后时间和冲击波半径的拟合
拟合得到
，这表明第一帧图像对应的时间是爆炸后0.s。同时还得到了
=.48，因此可以估算出爆炸释放的能量：
这相当于约吨TNT。荷兰特文特大学Sander博士也做了类似的计算，其给出的计算结果仅为吨左右[5]。而推特上的网友Sina通过计算给出的结果则高达吨[6]。
图9:谷仓高度的估计对爆炸能量计算的影响
以上能量的计算依赖于谷仓高度的估计，该估计是不确定性的最大来源。而由此确定的冲击波半径又极大地影响能量的计算（五次方），因此谷仓高度的估计是能量计算误差的主要来源。考虑谷仓高度的不确定性为
10米，即其高度在40至60米之间。重新应用上述过程进行计算，给出的爆炸当量在至吨（图9）。
结论量纲分析（点源爆炸模型）能够有效地估计出原子弹爆炸所释放的能量。年美国的科学家在新墨西哥州阿拉莫戈多沙漠进行的三位一体核试验的当量约为2万吨。通过量纲分析得到的原子弹爆炸当量与官方公布的数据仅相差10%。本文根据点源爆炸模型估算出贝鲁特港口大爆炸的当量为吨。考虑估算中的误差，当量的不确定性范围为至吨。附录
点击文章左下角“阅读原文”获取本文PDF版和附件。
参考资料[1]
Wikipediacontributors.beirutexplosions—Wikipedia,thefreeencyclopedia,.:
转载请注明：http://www.yinghuochonga.com/lnfz/5905.html
当前位置： 萤火虫 > 萤火虫的天敌 > 黎巴嫩港口大爆炸到底有多厉害
