自然災害與防災（進階）：把天災當成可量化的隨機過程

為什麼「百年一遇」的洪水，可能連續兩年都來？

你大概聽過氣象主播說某場暴雨是「兩百年一遇」，結果隔沒幾年同一條溪又淹了一次。許多人因此覺得專家在亂講。但問題其實不在專家，而在我們對「機率」的直覺出了錯。「兩百年一遇」並不是「每兩百年發生一次」的時刻表，而是「在任一年發生的機率約為 $1/200$」的統計陳述。一個機率為 $0.5\%$ 的事件，今年發生、明年再發生，並不違反任何規律——就像連續兩次擲到六點一樣完全可能。

入門篇談過地震、颱風、土石流的成因與基本防減災觀念。這一篇要往上走一層：把災害當作一個可量化的隨機過程（stochastic process）來分析。我們會處理三個核心問題——災害「多久來一次」（重現期與機率）、「來了會多嚴重」（規模—頻率律）、以及「我們能不能在它抵達前搶到幾秒」（地震早期預警的物理）。最後把這些拼成現代防災科學的骨架：風險 = 危害度 × 暴露度 × 脆弱度。

重現期與年超越機率：把直覺翻譯成數學

防災工程的第一個語言是重現期（return period, $T$）與年超越機率（annual exceedance probability, $P$）。兩者互為倒數：

$$P = \frac{1}{T}$$

所謂「$T = 200$ 年的洪水」，指的是流量達到或超過某個門檻值的事件，在任一年發生的機率是 $P = 1/200 = 0.5\%$。注意「達到或超過」這個詞很關鍵——重現期定義的是一個超越事件，不是恰好等於某值。

關鍵的盲點在於：人們把 $T$ 當成週期，但真實災害更接近無記憶（memoryless）的過程。如果每年發生與否彼此獨立，那麼「在 $n$ 年內至少發生一次」的機率不是 $n/T$，而要用補集來算：

$$P(\text{至少一次} \mid n\text{ 年}) = 1 - (1 - P)^n$$

動手算一下：你家的房子在使用年限內會遇到幾次大洪水？

假設你買了一棟設計壽命 $50$ 年的房子，位於「百年一遇」洪水（$P = 0.01$）的氾濫平原上。許多人直覺認為「我住不到一百年，應該安全」。我們算算看：

$$P(\text{50 年內至少一次}) = 1 - (1 - 0.01)^{50} = 1 - 0.99^{50}$$

$$= 1 - 0.605 = 0.395$$

也就是說，在房子壽命內遭遇至少一次百年洪水的機率高達 約 $40\%$——這跟「住不到一百年就安全」的直覺天差地遠。如果是設計壽命 $50$ 年的橋梁，工程師通常不會用百年重現期設計，而要求更高的標準（例如 $475$ 年或更長），原因正是這條公式。

我們也可以反過來問：若要把「使用年限內遭遇」的機率壓到 $10\%$ 以下，需要多大的重現期？由 $1 - (1 - 1/T)^{50} = 0.10$ 解得 $T \approx 475$ 年。這正是地震工程裡「$475$ 年重現期」這個常見數字的由來——它對應「$50$ 年內超越機率 $10\%$」這個國際通用的設計地震動標準。

規模—頻率律：大地震為什麼「稀有但致命」

第二個語言處理「災害有多大」。對地震而言，最著名的經驗律是古登堡—芮克特定律（Gutenberg–Richter law）：

$$\log_{10} N = a - b M$$

其中 $N$ 是規模大於等於 $M$ 的地震數目（在某時段、某區域內），$a$ 反映該區域的整體地震活躍度，$b$（稱為 $b$ 值）描述大小地震的比例關係，全球平均約為 $b \approx 1$。

$b \approx 1$ 意味著什麼？每當規模增加 $1$，地震次數就減少約十倍。規模 $5$ 的地震若一年發生 $100$ 次，規模 $6$ 約只有 $10$ 次，規模 $7$ 約 $1$ 次。但別忘了，規模每增加 $1$，釋放的地震矩能量增加約 $32$ 倍（因為矩規模 $M_w$ 與能量的關係是 $\log_{10} E \propto 1.5 M_w$）。所以大地震「次數十倍稀有、能量數十倍巨大」——這正是它們在統計上罕見卻在災害上主宰一切的原因。

看一個例子：用 b 值讀台灣的地震潛勢

$b$ 值不只是個常數，它本身就帶有物理意義。當一個區域的應力高度集中、或岩體趨於均質、接近破裂臨界時，大地震相對變多，$b$ 值會偏低（小於 $1$）；反之，破碎、應力分散的區域 $b$ 值偏高。火山活動區因岩漿活動產生大量小地震，$b$ 值常明顯大於 $1$。

對台灣這種位於歐亞板塊與菲律賓海板塊聚合帶的島嶼，地震學家會持續監測各斷層帶的 $b$ 值變化。若某段斷層的 $b$ 值在一段時間內系統性下降，可能暗示應力正在累積、進入孕震後期。當然，這只是潛勢評估的一環，不能直接拿來「預測」某日某時的地震——但它讓我們從「會不會有大地震」進化到「哪裡的大地震相對機率較高」。台灣的地震危害潛勢圖（如國震中心 NCREE 與中央氣象署發布的機率式地震危害分析）背後，正大量使用這類規模—頻率統計。

從「次數」到「搖晃」：機率式地震危害分析

把重現期與規模—頻率律結合，就得到現代防災的核心工具——機率式地震危害分析（Probabilistic Seismic Hazard Analysis, PSHA）。它要回答的不是「哪裡會地震」，而是更工程化的問題：

在未來 $50$ 年內，某地點的地表加速度（PGA）超過某個值的機率是多少？

PSHA 的邏輯可以拆成三步：

1. 震源：列出周邊所有可能的斷層與震源區，用古登堡—芮克特律給出各規模地震的年發生率。
2. 地動衰減：用地動預測方程式（Ground Motion Prediction Equation, GMPE）估計「規模 $M$ 的地震、距離 $R$」會在目標地點造成多大的地表震動。地震波能量隨距離擴散與吸收而衰減，這是一個帶不確定性的統計關係。
3. 積分：把所有可能的（規模、距離）組合，依其機率加權整合，得到「地表加速度超越某值」的年機率，畫成危害曲線（hazard curve）。

工程師再從危害曲線上讀取「$50$ 年超越機率 $10\%$」對應的設計地震加速度，據此訂定建築規範的耐震係數。這就是為什麼前面那條 $475$ 年重現期的公式會一路長進你頭頂的鋼筋混凝土裡。台灣的《建築物耐震設計規範》正是建立在 PSHA 的成果之上，並依不同地區的地震潛勢劃分震區係數。

搶那關鍵的幾秒：地震早期預警的物理

如果地震無法準確預測，防災還剩什麼能做？答案之一是早期預警（Earthquake Early Warning, EEW）——它不預測地震何時發生，而是在地震已經發生之後，搶在破壞性震波抵達前發出警報。這聽起來矛盾，物理上卻完全成立，關鍵在於兩種地震波的速度差。

地震發生時會輻射出兩種體波：

• P 波（縱波）：傳播快（地殼中約 $6$–$7$ km/s），但振動小、破壞力低，最先抵達。
• S 波（橫波）：傳播較慢（約 $3.5$–$4$ km/s），振動大、是造成搖晃與破壞的主力，較晚抵達。

預警系統靠近震央的測站先偵測到 P 波，立刻估算地震規模與位置，再趕在 S 波抵達遠處城市之前，透過電磁波（幾乎光速）把警報傳出去。警報跑得比地震波快——這就是 EEW 的物理本質。

動手算一下：你能搶到幾秒？

設震央距離某城市 $D = 100$ km。S 波抵達該城市所需時間：

$$t_S = \frac{D}{v_S} = \frac{100}{4} = 25 \text{ 秒}$$

P 波抵達時間：

$$t_P = \frac{D}{v_P} = \frac{100}{6.5} \approx 15.4 \text{ 秒}$$

假設預警系統需要約 $3$–$5$ 秒接收 P 波、判定參數並發出警報，而電磁波傳輸時間可忽略，則城市可獲得的預警時間（lead time）約為：

$$t_{\text{warn}} = t_S - t_P - t_{\text{處理}} \approx 25 - 15.4 - 4 \approx 5.6 \text{ 秒}$$

短短幾秒看似微不足道，卻足以讓高鐵自動減速、電梯停靠最近樓層開門、手術中止、瓦斯阻斷、學生就地掩護。距離越遠，預警時間越長——但這也帶出 EEW 無可迴避的弱點：

盲區（blind zone）。震央正上方與鄰近區域，P 波與 S 波幾乎同時抵達，警報還來不及發出搖晃就到了。偏偏震央附近往往正是搖得最猛、災情最重的地方。台灣中央氣象署的強震即時警報，在近震源區同樣受此物理限制。因此 EEW 是「遠距防護工具」，不是萬靈丹——這也是為什麼結構耐震設計（讓建築物本身扛得住）永遠是防災的第一道防線。

把全部拼起來：風險的三要素

到這裡，我們已經有足夠的工具回答防災最根本的問題：該把有限資源投在哪裡？現代災害科學用一條看似簡單、卻極具威力的關係式來組織思考：

$$\text{風險（Risk）} = \text{危害度（Hazard）} \times \text{暴露度（Exposure）} \times \text{脆弱度（Vulnerability）}$$

• 危害度（Hazard）：災害事件本身的強度與機率——也就是前面 PSHA、重現期所量化的部分。這是「天」的部分，人類幾乎無法改變。
• 暴露度（Exposure）：暴露在危害下的人、財產、基礎設施的數量與價值。同樣規模的地震，發生在無人荒野與發生在千萬人都會區，風險天差地別。
• 脆弱度（Vulnerability）：暴露對象受損的容易程度——老舊建築 vs 耐震新建築、有無預警系統、社區應變能力。

這條乘法式有個深刻的政策含義：因為是相乘，任一項趨近於零，總風險就趨近於零。 我們改變不了斷層（危害度幾乎固定），但可以透過都市計畫引導開發遠離高潛勢區（降低暴露度），更可以透過建築補強、土地利用管制、保險與演練大幅降低脆弱度。換句話說，防災的主戰場不在預測天災，而在降低我們自己的脆弱度。這也是為什麼聯合國《仙台減災綱領（Sendai Framework）》把重心從「災後應變」轉向「災前的脆弱度治理」。

重點回顧

• 重現期 $T$ 是機率的倒數（$P = 1/T$），不是時刻表；「$n$ 年內至少一次」要用 $1 - (1-P)^n$ 計算，結果常遠高於直覺。
• 古登堡—芮克特定律 $\log_{10} N = a - bM$ 描述規模—頻率關係；$b \approx 1$ 表示規模每增 $1$、次數減十倍，但能量增約 $32$ 倍，故大地震稀有卻主宰災情。
• PSHA 結合震源發生率與地動衰減，產出「危害曲線」，是耐震規範（如 $475$ 年重現期設計）的科學基礎。
• 地震早期預警 利用 P 波快、S 波慢的速度差搶得數秒至數十秒，但震央附近有無法避免的盲區。
• 風險 = 危害度 × 暴露度 × 脆弱度；因為相乘，降低人為可控的暴露度與脆弱度，是最有效的減災策略。

深入探討（研究所視角）

若想再往前推一層，有幾個方向值得探索。

地震是否真的「無記憶」？ 前面為簡化把地震當成泊松過程（Poisson process）——事件彼此獨立、發生率恆定，這也是 PSHA 的標準假設。但物理上，斷層在大地震後會釋放應力、需時間重新累積，這暗示了時間相依（time-dependent）的特性。彈性回跳理論（elastic rebound theory）與地震週期（seismic cycle）催生了 BPT（Brownian Passage Time）等更新過程模型（renewal model），其危害率（hazard rate）在剛發生大震後偏低、隨應力累積而上升，呈現「鐘形」的條件機率，與泊松過程的常數率截然不同。台灣車籠埔等主要斷層的長期潛勢評估，部分即採用此類時間相依模型。然而餘震序列又遵循大森定律（Omori's law）——餘震頻率隨時間 $t$ 衰減為 $n(t) \propto 1/(t+c)^p$，呈現強烈的時間群聚（clustering）。如何把「主震的長期更新過程」與「餘震的短期叢集」整合進統一框架，是當代地震統計（如 ETAS 模型，Epidemic-Type Aftershock Sequence）的核心課題。

重尾分布與極值統計。 災害損失往往不服從常態分布，而是重尾（heavy-tailed）——少數巨災主宰了長期總損失。處理「百年一遇」「千年一遇」這類稀有極端事件，需要極值理論（Extreme Value Theory），特別是廣義極值分布（GEV）與廣義柏拉圖分布（GPD）。它讓我們能在僅有數十年觀測資料的情況下，外推估計遠超過觀測期的重現水準（return level），但也伴隨巨大的外推不確定性——這正是「兩百年一遇」估計值彼此差異甚大的統計根源。

從危害到韌性（resilience）。 風險三要素是靜態快照，而真實社會是動態系統。近年研究將焦點從「降低脆弱度」延伸到「韌性」——系統在受擾後維持功能、快速回復的能力，常以韌性三角形（resilience triangle，功能損失隨時間的恢復曲線下面積）量化。把地球科學的危害量化、工程的脆弱度評估、與社會科學的恢復力治理整合成跨尺度的耦合模型，是減災科學從「抵抗災害」邁向「與災害共存」的前沿方向。