颱風夜的一場暴雨，雨水究竟流去哪裡？

颱風夜的一場暴雨，雨水究竟流去哪裡？

想像一個颱風夜，台北盆地三小時內降下超過 200 公釐的暴雨。屋頂上的雨水順著落水管流進街道側溝，側溝匯入雨水下水道，下水道再排進河川；河川水位逐漸抬升，逼近堤防頂端；如果上游某座水庫此刻正在洩洪，下游的流量還會再疊加。與此同時，自來水廠裡的清水正透過密布全市的壓力管網，從淨水場一路送到你家的水龍頭。

這整套「水從天上落下、在地表與管路中流動、最後排放或被利用」的過程，正是水利工程（hydraulic engineering）與水資源工程（water resources engineering）的核心。土木工程師在這個領域要回答的問題非常具體：這條排水溝該設計多大斷面，才能在 25 年一遇的暴雨下不溢淹？這根供水主幹管的管徑要多粗，水壓才足以送到頂樓？這座堤防要築多高，才能防住百年洪水？

要回答這些問題，工程師必須掌握水在不同情境下的兩種根本流態：在管路中承受壓力的封閉流動，以及在河川、渠道中有自由水面的明渠流動。這篇文章，我們就從支撐一切的兩條守恆定律談起，一路走到台灣最切身的防洪課題。

一切的起點：質量守恆與能量守恆

水流分析的地基，是兩條物理守恆定律。

第一條是連續方程式（continuity equation），它其實就是質量守恆。對不可壓縮流體（水可視為不可壓縮），通過任一封閉斷面的體積流率（discharge）$Q$ 必須相等：

$$Q = A_1 V_1 = A_2 V_2$$

其中 $A$ 是斷面積、$V$ 是平均流速。這條式子告訴我們一件直覺的事：管路變細時，水流必然加速；河道束窄處，水流變急。流量 $Q$ 一旦進入一段沒有分支的管路，從頭到尾都守恆。

第二條是白努利方程式（Bernoulli equation），它是能量守恆在流體中的展現。對一條流線上的兩點，在理想（無摩擦）情況下：

$$\frac{p_1}{\gamma} + \frac{V_1^2}{2g} + z_1 = \frac{p_2}{\gamma} + \frac{V_2^2}{2g} + z_2$$

式中 $p$ 是壓力、$\gamma$ 是水的單位重（約 $9.81\text{ kN/m}^3$）、$V$ 是流速、$g$ 是重力加速度、$z$ 是高程。三個項分別代表三種「水頭」（head，單位是長度公尺）：

• 壓力水頭（pressure head）$\dfrac{p}{\gamma}$
• 速度水頭（velocity head）$\dfrac{V^2}{2g}$
• 位置水頭（elevation head）$z$

白努利方程式的核心精神是：三種能量的總和守恆，彼此之間可以互相轉換。水流加速時速度水頭上升、壓力水頭就下降（這正是文氏管 Venturi 測流量的原理）；水從高處流下，位置水頭轉成速度水頭，水流就變急。

但真實的水流有摩擦。工程上必須在等式右邊補上一個水頭損失（head loss）$h_L$，這就把理想的白努利方程式修正為實用的能量方程式：

$$\frac{p_1}{\gamma} + \frac{V_1^2}{2g} + z_1 = \frac{p_2}{\gamma} + \frac{V_2^2}{2g} + z_2 + h_L$$

水利設計的大半工作，就是在估算這個 $h_L$ 到底有多大。

壓力管流：水頭損失怎麼算

當水充滿整根管子、靠壓力推動（例如自來水管），這叫壓力管流（pipe flow）。它的主要能量損失來自管壁摩擦，由 Darcy–Weisbach 方程式描述：

$$h_f = f \cdot \frac{L}{D} \cdot \frac{V^2}{2g}$$

其中 $h_f$ 是摩擦水頭損失、$L$ 是管長、$D$ 是管徑、$V$ 是流速、$f$ 是摩擦係數（friction factor）。這個 $f$ 不是常數，它取決於流體是層流還是紊流，而判斷依據是雷諾數（Reynolds number）：

$$Re = \frac{VD}{\nu}$$

$\nu$ 是水的運動黏滯度（約 $1.0 \times 10^{-6}\text{ m}^2/\text{s}$）。一般以 $Re < 2300$ 為層流、$Re > 4000$ 為紊流。對層流，摩擦係數有簡潔的理論解 $f = 64/Re$；對紊流，則需查 Moody 圖（Moody diagram）或用 Colebrook 公式迭代求解，並會牽涉到管壁的相對粗糙度 $\varepsilon/D$。

除了沿管的摩擦損失，管路的彎頭、閥門、縮放等配件還會造成次要損失（minor loss）：

$$h_m = K \cdot \frac{V^2}{2g}$$

$K$ 是各配件的損失係數。供水管網設計時，工程師要把所有 $h_f$ 與 $h_m$ 加總，確認最不利端（通常是最高樓層、最遠端）的剩餘壓力仍足以正常供水。

要提醒一個常見迷思：水管裡的壓力並不等於水的「流速」。一根末端被關閉的水管，內部可以有很高的靜水壓，但流速為零、水頭損失也為零。壓力與流動是兩回事——能量方程式同時追蹤壓力水頭與速度水頭，正是為了把這兩者分清楚。

明渠流：有自由水面的世界

河川、排水溝、灌溉渠道、雨水下水道——這些水流的上方是開放的大氣，水面可以自由升降，稱為明渠流（open channel flow）。它和壓力管流最大的差別是：明渠流的驅動力不是壓力，而是重力沿渠床坡度的分量，而且水深本身是未知數，會隨流量自我調整。

工程實務上，最常用的明渠流公式是 Manning 公式：

$$V = \frac{1}{n} R^{2/3} S^{1/2}$$

或以流量表示：

$$Q = \frac{1}{n} A R^{2/3} S^{1/2}$$

其中各符號意義為：

• $n$：曼寧粗糙係數（Manning's roughness coefficient），混凝土渠約 $0.013$、天然土渠約 $0.025$、雜草叢生的河道可達 $0.05$ 以上。渠面越粗糙，$n$ 越大、水流越慢。
• $R$：水力半徑（hydraulic radius），定義為通水斷面積 $A$ 除以濕周 $P$，即 $R = A/P$。它代表斷面「過水效率」的高低。
• $S$：渠床坡度（無因次）。

明渠流還有一個壓力管流沒有的核心概念：流態。我們用福祿數（Froude number）來區分：

$$Fr = \frac{V}{\sqrt{gy}}$$

$y$ 是水深。$Fr < 1$ 為緩流（subcritical flow，水深大、流速慢，下游擾動能往上游傳遞）；$Fr > 1$ 為急流（supercritical flow，水淺流急，擾動無法逆流而上）；$Fr = 1$ 為臨界流。當水流從急流突然轉為緩流時，會發生水躍（hydraulic jump）——水面驟然抬升、劇烈翻騰並消散大量能量。工程師反而常常「故意製造」水躍，把溢洪道高速下洩的水流能量在消能池中耗掉，避免下游河床被沖刷掏空。

防洪：把水文與水力縫在一起

前面講的是「水怎麼流」（水力學 hydraulics），但防洪設計還需要先回答「會來多少水」（水文學 hydrology）。

最基礎的橋樑是合理化公式（Rational Method），用來推估小集水區的洪峰流量：

$$Q_p = \frac{1}{360} \, C \, i \, A$$

（採台灣慣用單位：$Q_p$ 為 $\text{m}^3/\text{s}$、$C$ 為逕流係數、$i$ 為降雨強度 $\text{mm/hr}$、$A$ 為集水面積 $\text{ha}$。）

逕流係數 $C$ 反映地表的「不透水程度」：森林約 $0.1$、農田約 $0.3$，而都市的柏油與屋頂可高達 $0.9$。這條公式背後藏著台灣防洪最沉重的一課——都市化讓 $C$ 大幅上升。同樣一場暴雨，過去農田能吸收大半雨水、慢慢入滲，如今水泥地讓九成雨水瞬間化為地表逕流，又快又猛地灌入下水道，這正是都市內澇日益嚴重的根本原因。

防洪設計還必須選定一個設計重現期（return period）。台灣的區域排水通常以 25 年重現期設計、主要河川堤防多以 100 年重現期設計，而像曾文水庫這類大型水庫的溢洪道，則要能安全宣洩萬年重現期等級的可能最大洪水（PMF）。重現期越長代表防護標準越高、工程造價也越貴，這之間的取捨是公共政策與工程經濟的交會點。

值得強調的是，現代防洪觀念已從「把水盡快排掉」轉向「與水共存」。海綿城市（sponge city）、低衝擊開發（low impact development, LID）、滯洪池、透水鋪面、雨水花園等做法，都是刻意降低逕流係數、把暴雨「留下來慢慢放」，減輕下游排水系統的瞬間壓力。在台灣這個颱風頻繁、坡陡流急、都市又高度集中於少數平原的環境下，這種思維轉變格外重要。

看一個例子

讓我們把水文與水力串起來，設計一條都市排水箱涵。

問題：某都市集水區面積 $A = 50\text{ ha}$，逕流係數 $C = 0.7$，設計降雨強度 $i = 100\text{ mm/hr}$。試求設計洪峰流量；接著設計一條混凝土矩形箱涵（$n = 0.013$，渠床坡度 $S = 0.001$）來宣洩此流量，假設斷面採用「水力最佳斷面」（寬 $b$ 為水深 $y$ 的兩倍，即 $b = 2y$）。

第一步：用合理化公式求洪峰流量。

$$Q_p = \frac{1}{360} \, C \, i \, A = \frac{1}{360} \times 0.7 \times 100 \times 50 \approx 9.7\text{ m}^3/\text{s}$$

第二步：建立箱涵的幾何關係。 矩形斷面寬 $b = 2y$，則通水斷面積與濕周為：

$$A = b \cdot y = 2y \cdot y = 2y^2$$ $$P = b + 2y = 2y + 2y = 4y$$

水力半徑：

$$R = \frac{A}{P} = \frac{2y^2}{4y} = \frac{y}{2}$$

第三步：代入 Manning 公式求水深。

$$Q = \frac{1}{n} A R^{2/3} S^{1/2}$$

$$9.7 = \frac{1}{0.013} \times 2y^2 \times \left(\frac{y}{2}\right)^{2/3} \times (0.001)^{1/2}$$

先整理常數項：$\dfrac{1}{0.013} = 76.9$、$(0.001)^{1/2} = 0.0316$、$\left(\dfrac{1}{2}\right)^{2/3} = 0.63$。代入：

$$9.7 = 76.9 \times 2 \times 0.63 \times 0.0316 \times y^{2 + 2/3}$$

$$9.7 = 3.06 \, y^{8/3}$$

$$y^{8/3} = \frac{9.7}{3.06} = 3.17 \implies y = 3.17^{3/8} \approx 1.56\text{ m}$$

於是水深約 $1.56\text{ m}$、渠寬 $b = 2y \approx 3.12\text{ m}$。

第四步：檢核流速與流態。

$$V = \frac{Q}{A} = \frac{9.7}{2 \times 1.56^2} = \frac{9.7}{4.87} \approx 1.99\text{ m/s}$$

$$Fr = \frac{V}{\sqrt{gy}} = \frac{1.99}{\sqrt{9.81 \times 1.56}} = \frac{1.99}{3.91} \approx 0.51$$

$Fr < 1$，為緩流，水流穩定不會產生急流沖刷的問題；流速約 $2\text{ m/s}$ 也在混凝土渠道可接受範圍內（既不會慢到淤積，也不會快到磨損）。設計合理。

延伸思考：如果這片集水區未來持續開發，$C$ 從 $0.7$ 上升到 $0.85$，洪峰流量會增加約兩成，原本剛好夠用的箱涵可能就會溢淹。這正說明了為什麼防洪設施不能只看「現在」，還必須預留未來都市化與氣候變遷的成長空間。

重點回顧

1. 連續方程式（$Q = AV$）與能量方程式是水流分析的兩大基石；能量方程式比理想的白努利方程式多了一個水頭損失 $h_L$，而估算 $h_L$ 正是水利設計的核心工作。
2. 壓力管流靠壓力驅動，摩擦損失用 Darcy–Weisbach 公式（$h_f = f \frac{L}{D}\frac{V^2}{2g}$）計算，摩擦係數 $f$ 隨雷諾數與管壁粗糙度而變。
3. 明渠流靠重力沿坡度驅動，常用 Manning 公式（$V = \frac{1}{n}R^{2/3}S^{1/2}$）；福祿數 $Fr$ 區分緩流與急流，急流轉緩流會發生消能的水躍。
4. 防洪設計要先以水文方法（如合理化公式 $Q_p = \frac{1}{360}CiA$）推估會來多少水，再以水力方法設計排水斷面；都市化推升逕流係數 $C$ 是內澇加劇的主因。
5. 現代防洪從「快排」轉向「與水共存」——海綿城市與低衝擊開發刻意降低逕流、滯洪緩排，對颱風頻繁、地形陡峭的台灣尤其重要。

深入探討（研究所視角）

進入研究所層級，水利工程的視野會從「會用公式設計斷面」擴展到「用偏微分方程描述非穩態流動」，並走向數值模擬與不確定性分析。

從穩態到非穩態：聖維南方程組。 前文的 Manning 公式假設流況是穩定均勻流，但真實洪水是一波隨時間移動的水體。完整描述明渠非穩態流動的是聖維南方程組（Saint-Venant equations），由連續方程式與動量方程式構成的一階雙曲型偏微分方程組：

$$\frac{\partial A}{\partial t} + \frac{\partial Q}{\partial x} = 0$$

$$\frac{\partial Q}{\partial t} + \frac{\partial}{\partial x}\left(\frac{Q^2}{A}\right) + gA\frac{\partial h}{\partial x} = gA(S_0 - S_f)$$

其中 $S_0$ 是渠床坡度、$S_f$ 是摩擦坡降。這組方程式沒有解析解，現代洪水演算與河川模擬軟體（如 HEC-RAS、SOBEK、SRH-2D）都是用有限差分或有限體積法數值求解它。依保留項目的多寡，又衍生出運動波、擴散波、動力波等簡化模型，工程師需依問題特性（是緩坡長河還是陡峭山溪）選擇適當的近似。

水文模式的參數不確定性。 合理化公式只是入門，真正的集水區降雨—逕流模擬會用到單位歷線（unit hydrograph）、SCS 曲線數法（curve number）、或概念式水筒模式。這些模式的參數（入滲率、滯時、儲蓄係數）難以直接量測，必須透過率定（calibration）與驗證（validation）反推。台灣山高水急、集水區反應時間極短，加上雨量站空間分布不均，使得參數不確定性與率定難度都偏高，這也是為何防洪設計要保守地納入安全係數。

氣候變遷下的非定態水文。 傳統重現期分析建立在「水文序列是定態（stationary）」的假設上——也就是假設過去的統計分布能代表未來。但氣候變遷正在打破這個前提：極端降雨的頻率與強度都在改變，過去的「百年洪水」可能在未來變成「五十年洪水」。研究前沿因此轉向非定態頻率分析（non-stationary frequency analysis），讓水文分布的參數隨時間演變，並結合區域氣候模式（RCM）的降尺度輸出來推估未來情境。對於整個本島都暴露在颱風與短延時強降雨下的台灣，如何在設計規範中納入氣候變遷調適，是當前水利工程界最迫切的活躍課題。

值得延伸的思考：當你下次經過一條河川或一座滯洪池，可以試著想像那條看不見的能量線（energy grade line）如何沿著水流逐漸下降，以及設計者在「防護標準」與「工程造價」之間做了怎樣的權衡。在氣候越來越極端的時代，理解水如何流動、又該如何被引導，已不只是工程專業，更是每個生活在這片土地上的人都該具備的素養。