假如那團氣雲真的點燃了，破壞半徑會到哪裡？

假如那團氣雲真的點燃了，破壞半徑會到哪裡？

入門篇我們學會了一句心法：風險 = 頻率 × 後果。當時我們把「後果」輕輕帶過——說它「用輸送方程算出來」就算交代了。但真正的製程安全工程師，每天面對的最棘手、也最迷人的問題，恰恰是這個被輕輕帶過的「後果」：

一座儲槽破了一個 50 mm 的洞，裡面是 20 bar 的液態丙烷。它會洩漏多快？洩出來的液體有多少瞬間閃蒸（flash）成蒸氣？這團比空氣重的蒸氣雲會貼著地面飄多遠、在多遠的距離仍維持在爆炸下限（LEL）以上？如果它在 200 公尺外遇到一輛汽車的火花而引爆，產生的爆震超壓（overpressure）足以震碎廠界外住宅的玻璃嗎？還是會壓垮整面承重牆？

這一連串問題，沒有一個能靠「貼標語」或「多裝一個閥」回答。它們需要把孔口流、閃蒸熱力學、重氣擴散、燃燒爆轟、熱輻射一路串起來，算出一個個帶單位的數字，再用這些數字去回答最現實的工程決策：廠房之間該隔多遠？控制室的牆要做多厚？社區該退到哪條線以外？ 這就是進階製程安全的核心——後果模型（consequence modeling）。本文假設你已熟悉危害與風險的區別、HAZOP 與保護層的概念，我們直接潛入「後果」這口深井。

後果分析的骨架：從源項到傷害的因果鏈

任何一場洩漏事故的後果，都可以拆解成一條清晰的因果鏈，每一環都對應一套物理模型。理解這條鏈，是進階後果分析的地圖：

$$\text{源項（source term）} \rightarrow \text{擴散（dispersion）} \rightarrow \text{點燃／毒性} \rightarrow \text{物理效應} \rightarrow \text{傷害（harm）}$$

• 源項：破口多大、洩漏速率多少、洩出的是液體還是氣體、有多少閃蒸。這是質能平衡與流體力學的問題。
• 擴散：洩出的物質在大氣中如何被風與紊流稀釋、濃度場如何隨距離與時間演變。這是質量輸送的問題。
• 點燃或毒性：若是可燃物，遇火源會閃火（flash fire）、爆炸（explosion）或噴射火（jet fire）；若是毒物，則直接以濃度造成傷害。
• 物理效應：爆炸產生超壓、火災產生熱輻射、毒氣產生致死濃度暴露。這是能量釋放與傳遞的問題。
• 傷害：超壓、熱通量、毒性劑量如何轉換成人員傷亡與設備損壞的機率。這由機率單位模型（probit model）量化。

入門篇把整條鏈濃縮成一個孔口流公式加一句「擴散模型」。進階的功夫，在於把每一環都建立可計算的方程，並理解它們各自的適用範圍與陷阱。我們逐環來看。

源項：洩漏的不只是流量，還有閃蒸

入門篇給過孔口流公式 $\dot{m} = C_d A \sqrt{2\rho \Delta P}$，那是液體穩態洩漏的近似。但真實源項遠比這複雜，這裡點出兩個進階關鍵。

第一，氣體洩漏會「壅塞（choke）」。 當容器內壓力高到一定程度，氣體通過破口的速度會達到該處的音速而無法再增加——這叫壅塞流（choked flow）。判據是壓力比超過臨界值：

$$\frac{P_{\text{atm}}}{P_{\text{in}}} \le \left(\frac{2}{\gamma+1}\right)^{\gamma/(\gamma-1)}$$

其中 $\gamma = C_p/C_v$ 是比熱比（空氣約 1.4）。代入可得臨界壓力比約 0.53。也就是說，只要容器內絕對壓力大於外界約兩倍（如多數加壓儲槽），洩漏就處於壅塞狀態，質量流率只取決於上游壓力與溫度，與下游大氣壓無關：

$$\dot{m}_{\text{choked}} = C_d A P_{\text{in}} \sqrt{\frac{\gamma M}{R T_{\text{in}}}\left(\frac{2}{\gamma+1}\right)^{\frac{\gamma+1}{\gamma-1}}}$$

這個「與下游無關」的特性很反直覺，卻是高壓氣體洩漏估算的基礎。

第二，過熱液體會「閃蒸」。 許多危險物（丙烷、LPG、液氨）在儲槽中是加壓液化的。一旦破口讓壓力驟降到常壓，這些液體瞬間處於遠高於其常壓沸點的「過熱」狀態，會有一部分立刻汽化來吸收多餘的顯熱。這個閃蒸分率（flash fraction）可由能量平衡估算：

$$f_v = \frac{C_p (T_{\text{store}} - T_{\text{boil}})}{\Delta H_{\text{vap}}}$$

直觀地說：分子是液體從儲存溫度降到常壓沸點所釋放的顯熱，分母是汽化潛熱，兩者之比就是必須汽化掉以帶走顯熱的那部分。這一步至關重要，因為閃蒸產生的蒸氣往往還會夾帶等量甚至更多的液滴霧（aerosol），使實際進入氣雲的物質量遠大於單純的閃蒸分率——這是很多事故後果被低估的根源。源項算錯，後面整條鏈全錯。

擴散：當氣雲比空氣重，高斯模型就失效了

洩漏物進入大氣後，被風與紊流帶著走、同時被稀釋。入門篇提過「擴散模型本質是對流–擴散方程」，進階的關鍵在於分辨兩種截然不同的擴散行為。

中性或輕氣體（與空氣密度相近）遵循經典的高斯煙羽模型（Gaussian plume model）。對一個連續點源，下風處 $(x,y,z)$ 的濃度為：

$$C(x,y,z) = \frac{\dot{m}}{2\pi u\, \sigma_y \sigma_z} \exp\!\left(-\frac{y^2}{2\sigma_y^2}\right)\left[\exp\!\left(-\frac{(z-H)^2}{2\sigma_z^2}\right) + \exp\!\left(-\frac{(z+H)^2}{2\sigma_z^2}\right)\right]$$

其中 $u$ 是風速、$H$ 是有效釋放高度，$\sigma_y, \sigma_z$ 是橫向與垂直擴散參數，隨下風距離與大氣穩定度（Pasquill stability class A–F）增長。第二個指數項是地面反射的鏡像源，代表氣體碰到地面後反彈。這個模型乾淨、解析、好用——但它有一個致命前提：被釋放物的密度與空氣相近、會乖乖隨大氣紊流飄散。

問題是，化工最危險的洩漏物——氯、氨、丙烷、大多數碳氫化合物蒸氣——都比空氣重（尤其閃蒸後溫度極低、密度更大）。重氣（dense gas）的行為完全不同：它不會升空稀釋，而是像看不見的洪水一樣貼著地面流動、向四周塌陷蔓延，在低窪處積聚，被稀釋的速度遠比高斯模型預測的慢。博帕爾的 MIC 氣雲、許多氯氣洩漏的高致死範圍，都是重氣行為的後果。對這類洩漏，工程上改用重氣擴散模型，最經典的是 Britter–McQuaid 關聯式，它以無因次群整理風洞與實場數據，預測重氣雲沿地面的濃度衰減。若用高斯模型去算一團重氣，會嚴重低估近場的致死範圍——這是後果分析最常見、也最危險的誤用。

進階工程師因此會先做一個判斷：算一個理查森數（Richardson number）之類的密度判據，決定該用高斯還是重氣模型。選錯模型，比不算還危險，因為它給了一個看似精確、實則樂觀的假安全感。

物理效應之一：氣雲爆炸的超壓有多大

假設那團可燃氣雲飄到一個有火源的地方點燃了。如果它在開闊處快速燒掉，是一場熱輻射為主的閃火；但若氣雲處在有管架、設備、建築物的壅塞空間，火焰加速、產生顯著的壓力波——這就是蒸氣雲爆炸（vapour cloud explosion, VCE），工業界最具破壞力的事故型態之一。

後果分析要回答的是：爆炸在不同距離產生多大的超壓（overpressure）$\Delta P$？最簡單的工程模型是 TNT 當量法（TNT equivalence）。它把氣雲中參與爆炸的可燃物能量，折算成等效的 TNT 質量：

$$W_{\text{TNT}} = \frac{\eta \, m_f \, \Delta H_c}{E_{\text{TNT}}}$$

其中 $m_f$ 是參與爆炸的可燃物質量、$\Delta H_c$ 是其燃燒熱、$E_{\text{TNT}} \approx 4.6\ \text{MJ/kg}$ 是 TNT 的爆炸能、$\eta$ 是爆炸效率（explosion yield）——一個經驗折減因子，通常只取 1%–10%，因為氣雲爆炸遠不如固體炸藥有效率。接著用一個比例距離（scaled distance）把不同當量的爆炸整理到同一條曲線上：

$$z_{\text{scaled}} = \frac{R}{W_{\text{TNT}}^{1/3}}$$

這個 $W^{1/3}$ 的縮放（Hopkinson–Cranz scaling）是爆震物理的核心：爆炸的破壞範圍隨能量的立方根增長。也就是說，能量放大 8 倍，相同超壓的破壞半徑只放大 2 倍。查一條標準的超壓 vs 比例距離曲線，就能得到任意距離的 $\Delta P$。

有了超壓，再對照破壞門檻就能評估後果：

TNT 當量法的局限與更進階的模型：TNT 法把氣雲爆炸當成「一個點的固體炸藥」，但氣雲爆炸的壓力波形狀其實平緩得多、且高度依賴空間壅塞程度。因此現代後果分析更常用 TNO 多能量法（multi-energy method） 或 Baker–Strehlow–Tang 法，它們明確考慮「壅塞區體積」與「火焰加速程度」，對近場超壓的預測比 TNT 法合理得多。但 TNT 法因其簡單，仍是手算與初步篩選的常用工具。

物理效應之二：火災的熱輻射與致命的 BLEVE

不是所有後果都來自爆震。火災透過熱輻射（thermal radiation）傷人毀物。後果分析的關鍵量是某距離處接收到的入射熱通量（incident heat flux）$q''$（單位 $\text{kW/m}^2$）。對一個池火（pool fire）或火球，常用點源或固體火焰模型估算：

$$q'' = \frac{\tau \, \eta_r \, \dot{m}_f \, \Delta H_c}{4\pi R^2}$$

分子是火源輻射出的總熱功率（燃燒速率 × 燃燒熱 × 輻射分率 $\eta_r$，再乘大氣穿透率 $\tau$），分母的 $4\pi R^2$ 代表能量隨距離平方稀釋地散布到球面上——這就是輻射的平方反比律。對照人體與設備的熱通量門檻（如 $1.6\ \text{kW/m}^2$ 長時間暴露才會疼痛、$12.5\ \text{kW/m}^2$ 足以點燃木材並造成重度燒傷、$37.5\ \text{kW/m}^2$ 使設備失效），就能畫出火災的危害等高線。

火災後果中最戲劇性的是 BLEVE（Boiling Liquid Expanding Vapour Explosion，沸騰液體膨脹蒸氣爆炸）。當一個盛裝過熱液體（如 LPG）的儲槽長時間被外部火焰炙烤，槽壁強度因高溫下降、同時內壓因液體受熱上升，最終槽體結構性破裂。瞬間，整槽過熱液體暴露於常壓、爆發性閃蒸，若是可燃物便立刻點燃成一個巨大火球（fireball），同時拋出致命的金屬碎片。BLEVE 之所以可怕，在於它同時產生超壓、巨大熱輻射火球、與飛射物三種效應，且常發生在消防人員已抵達現場、自以為控制住火勢的時刻。理解 BLEVE 的機制，正是「源項中的閃蒸熱力學」與「火災熱輻射」兩條線索的交會——這也是為什麼進階後果分析要把整條因果鏈打通來看。

從物理效應到傷害：機率單位（probit）模型

算出了超壓、熱通量或毒物濃度，最後一步是把這些物理量轉換成傷害機率。人不是「超過某值就死、低於就活」的開關，而是有統計分布的脆弱性。後果分析用機率單位模型（probit model）描述這種劑量–反應關係：

$$Y = a + b \ln(V)$$

其中 $V$ 是致害負荷（causative variable）——對毒物是「濃度的 $n$ 次方乘暴露時間」$\int C^n\,dt$，對熱輻射是「熱通量的 4/3 次方乘暴露時間」，對超壓則是超壓本身。$a, b$ 是各種傷害（死亡、肺出血、玻璃破裂）的實驗回歸常數。算出的 $Y$（probit 值）再透過標準常態分布轉成傷害機率 $P$。

probit 模型的價值，在於它把「物理效應等高線」轉成「傷害機率等高線」，讓後果分析能與頻率相乘，最終回到入門篇那條 $\text{風險} = \text{頻率} \times \text{後果}$，只是現在「後果」是一個帶機率、帶空間分布的嚴謹結果，而不只是一個籠統的形容詞。

看一個例子：估算一場丙烷氣雲爆炸的破壞半徑

讓我們把上面幾條鏈打通，做一個量級估算。某管線洩漏形成可燃氣雲，估計參與爆炸的丙烷質量 $m_f = 1000\ \text{kg}$，丙烷燃燒熱 $\Delta H_c = 46\ \text{MJ/kg}$，取保守的爆炸效率 $\eta = 0.05$（5%）。問：產生「玻璃破裂」（$\Delta P \approx 1\ \text{kPa}$，對應比例距離約 $z \approx 15\ \text{m/kg}^{1/3}$）的距離大約在哪裡？

第一步，算 TNT 當量。 取 $E_{\text{TNT}} = 4.6\ \text{MJ/kg}$：

$$W_{\text{TNT}} = \frac{\eta\, m_f\, \Delta H_c}{E_{\text{TNT}}} = \frac{0.05 \times 1000 \times 46}{4.6} = \frac{2300}{4.6} = 500\ \text{kg TNT}$$

第二步，由比例距離反推實際距離。 已知 $z_{\text{scaled}} = R / W_{\text{TNT}}^{1/3}$，先算立方根：

$$W_{\text{TNT}}^{1/3} = 500^{1/3} \approx 7.94$$

於是玻璃破裂距離：

$$R = z_{\text{scaled}} \times W_{\text{TNT}}^{1/3} \approx 15 \times 7.94 \approx 119\ \text{m}$$

結論：這場爆炸大約在 120 公尺處仍足以震碎玻璃。若廠界距此洩漏點不到 120 公尺，廠外民宅就會受影響——這個數字直接決定了「廠界該退多遠」「控制室該設在哪」的工程決策。

注意這個估算的靈敏度：爆炸效率 $\eta$ 從 5% 改成 10%，當量加倍到 1000 kg，但破壞半徑只乘上 $2^{1/3} \approx 1.26$，變成約 150 m。立方根縮放讓後果對能量不那麼敏感，卻也意味著「想用減少存量來縮小破壞圈」的效益是遞減的——存量砍一半，破壞半徑只縮小 20%。這正是後果分析能給設計者的深刻洞見，也再次印證入門篇「減量」策略的價值與其物理上限。

後果模型如何餵養設施配置與風險等高線

把後果分析做完，工程師手上會有一張張危害等高線：超壓 21 kPa 的圈、熱通量 $12.5\ \text{kW/m}^2$ 的圈、致死濃度 LC50 的圈。這些圈疊上廠區與社區地圖，就成了製程安全最具體的產出——以後果為基礎的設施配置（consequence-based facility siting）。

更進一步，把每一種洩漏情境的後果等高線與其發生頻率結合，對廠區周圍每一個點累加所有情境的致死機率，就得到一張個人風險等高線圖（individual risk contour）——地圖上每一點標注「一個長期站在這裡的人每年的死亡機率」。再把人口分布納入，就能畫出 F–N 曲線（社會風險曲線）：橫軸是死亡人數 $N$、縱軸是「造成至少 $N$ 人死亡的事故年頻率 $F$」。法規會在 F–N 圖上劃出「不可接受」與「可忽略」兩條線，中間是 ALARP（As Low As Reasonably Practicable，合理可行下盡可能低）區。

這就把入門篇抽象的「個人風險 $\le 10^{-6}/\text{yr}$」變成了一張可以直接拿來做決策的地圖——而支撐這張地圖的每一條等高線，都是前面源項、擴散、爆炸、輻射、probit 一路算出來的。後果模型不是學術練習，它是把「廠界該畫在哪、社區該退到哪、控制室牆該多厚」這些攸關人命的問題，從直覺猜測變成可辯護的工程計算。

後果優化：把破壞半徑寫進設計目標

進階後果分析最終要回到優化學。入門篇說「安全是優化的硬約束」，這裡我們把它精確化:後果模型讓「安全約束」變成設計變數的顯函數，從而真正可被優化器處理。

考慮一個典型的設施配置優化問題:在廠區內安排 $n$ 個設備單元的位置 $\{x_i, y_i\}$，目標是最小化「總管線成本 + 期望災損」，同時滿足安全間距約束：

$$\min_{\{x_i,y_i\}} \;\; \sum_{i,j} c_{ij}\,d_{ij} \;+\; \sum_{k} f_k \cdot \text{Loss}_k(\{x_i,y_i\})$$

$$\text{s.t.}\quad d_{ij} \ge D_{ij}^{\min}(\Delta P_{\text{threshold}}),\qquad \text{IR}(p) \le 10^{-6}/\text{yr}\ \ \forall p \in \text{廠界外}$$

其中 $d_{ij}$ 是單元間距、$D_{ij}^{\min}$ 是由爆炸超壓後果模型反推出的最小安全間距（單元 $i$ 爆炸時不致引發單元 $j$ 二次破壞的距離），而 $\text{Loss}_k$ 中的災損正比於後果模型算出的破壞圈面積。這裡的關鍵是：後果模型把「安全」翻譯成了距離與面積的可微（或至少可計算）函數，優化器於是能在「設備靠近以省管線」與「設備分散以降連鎖爆炸風險」之間找到數學上的最佳平衡。

更深一層，本質安全的「減量」在此有了量化的優化解讀。前面例子顯示破壞半徑隨存量的立方根增長，因此「為了縮小破壞圈而減量」的邊際效益遞減；優化器會自動在「減量的製程成本」與「破壞圈縮小帶來的災損與用地節省」之間權衡，找出最優存量——這比入門篇「減量總是好的」更精確：減量有最優點，而後果模型告訴你那個點在哪。製程模擬軟體（如 Aspen、PHAST、Safeti）的後果模組，本質上就是在每組候選設計上跑完整條「源項→擴散→效應→傷害」鏈，把結果回饋給優化器當約束與目標——這正是現代化工把安全、經濟、配置三者整合求解的縮影。

重點回顧

• 後果分析是一條因果鏈：源項 → 擴散 → 點燃/毒性 → 物理效應 → 傷害。每一環各有物理模型，進階的功夫在於把整條鏈打通並理解各環的適用範圍。
• 源項不只看流量：高壓氣體洩漏會壅塞（質量流率與下游壓力無關）；過熱液體洩漏會閃蒸並夾帶液滴霧，實際進入氣雲的量常被嚴重低估。
• 擴散要先分密度：中性氣體用高斯煙羽模型，但化工最危險的重氣會貼地蔓延、稀釋極慢，須用 Britter–McQuaid 等重氣模型；用錯模型會給出樂觀的假安全感。
• 爆炸與火災各有度量：VCE 用 TNT 當量法（破壞半徑隨能量立方根增長）或更精確的多能量法算超壓；火災用平方反比的熱輻射通量；BLEVE 同時帶來超壓、火球與飛射物。
• probit 把物理效應轉成傷害機率：再結合頻率與人口分布，產出個人風險等高線與 F–N 社會風險曲線，支撐設施配置與 ALARP 決策。
• 後果模型讓安全成為可優化的顯函數：設施配置優化在「省管線/省地」與「降連鎖風險」間求最佳解；減量的邊際效益隨立方根遞減，存在數學最優點。

深入探討（研究所視角）

研究所層級的後果模型，會從「經驗關聯式與查曲線」走向「第一原理的高解析數值模擬」與「不確定性下的嚴謹決策」。幾個值得深入的前沿：

1. CFD 取代積分模型（high-fidelity dispersion & explosion CFD）。 高斯與 Britter–McQuaid 是積分模型（integral model）——用簡化關聯式快速給答案，但無法處理真實廠區的複雜地形與設備壅塞。前沿研究改用計算流體力學（CFD）直接求解雷諾平均（RANS）或大渦模擬（LES）下的 Navier–Stokes 方程，配合燃燒模型（如 FLACS 軟體），逐格算出氣雲在真實管架間如何流動、火焰如何在壅塞中加速、超壓場的真實空間分布。這把後果分析從「平坦開闊地的近似」提升到「這座具體工廠的數值風洞」，但代價是巨大的計算量與對亂流/燃燒次網格模型的依賴——準確度與計算成本的取捨，本身就是活躍的研究課題。

2. 重氣與紊流燃燒的次網格物理（sub-grid physics of dense gas & turbulent combustion）。 重氣擴散的核心難點是密度分層抑制紊流混合，而 VCE 的核心是火焰–紊流交互作用（flame–turbulence interaction）導致火焰加速乃至爆轟轉變（deflagration-to-detonation transition, DDT）。這兩者都涉及跨尺度的物理：大尺度的流場與小尺度的混合/反應強烈耦合。研究者用 LES 配合火焰面密度（flame surface density）或增厚火焰（thickened flame）模型來捕捉這些次網格現象，是燃燒學與安全工程的交界前沿。

3. 不確定性量化與貝氏更新（uncertainty quantification & Bayesian inference）。 後果模型的每個輸入——破口尺寸、風速、穩定度、爆炸效率 $\eta$——都帶有巨大不確定性，$\eta$ 從 1% 到 10% 的範圍就能讓災損估計差一個量級。前沿研究用不確定性量化（UQ）——蒙地卡羅、多項式混沌展開（polynomial chaos）、代理模型（surrogate model）——把這些不確定性傳遞到最終風險，給出的不是一個點估計而是一個風險分布；並用事故後的實測數據透過貝氏推論反過來更新模型參數，讓後果模型隨經驗持續校準。

4. 連鎖效應與系統層級的「骨牌風險」（domino effects & spatial risk propagation）。 真實災難常是一個單元的爆炸/火災引燃鄰近單元，形成連鎖效應（domino effect）——這正是前面設施配置約束 $D_{ij}^{\min}$ 要防的事。研究者用圖論與機率傳播模型，把廠區建模成一個「節點為設備、邊為事故傳播機率」的網路，分析整廠的連鎖風險，並反過來優化配置以阻斷傳播路徑。這把後果分析從「單一情境」提升到「系統層級的空間風險網路」，與入門篇提到的 STAMP 系統觀遙相呼應。

把這四個方向串起來，你會看到後果模型早已不是一本查表手冊，而是融合計算流體力學、燃燒物理、不確定性量化與網路科學的活躍研究領域。它的核心命題始終如一：在災難尚未發生時，盡可能精確地預測「假如它發生，傷害會到哪裡」——因為唯有把這個數字算清楚，我們才談得上理性地把廠界、社區與生命，安放在安全的那一側。