抽一口看不見的水：為什麼台灣高鐵會「沉下去」？

抽一口看不見的水：為什麼台灣高鐵會「沉下去」？

如果你搭過台灣高鐵經過雲林、彰化一帶，可能不會注意到一件事：這段軌道底下的地面，正以每年數公分的速度緩慢下沉。最嚴重的時期，雲林土庫一帶的累積沉陷曾逼近高鐵營運安全的紅線，逼得政府投入巨資監測、限制地下水抽取。地面為什麼會「消失」？答案不在地表，而在你看不見的地底——那裡有一層層飽含水的砂石，當人們把水抽走，土壤骨架就被壓密、地表隨之塌陷。

入門篇我們談的是「水在地表怎麼流」——河川、管路、排水溝。但水資源工程真正棘手的一半，發生在地表之下：地下水（groundwater）佔了台灣總用水量約兩成、在中南部枯水期更是農業命脈。要理解地層下陷、海水入侵、井群干擾這些問題，我們需要一套全新的物理框架：多孔介質中的流動（flow through porous media）。這篇文章，我們就鑽進地底，看看水如何在砂石的縫隙間滲流，以及工程師如何在「用水」與「護地」之間走鋼索。

達西定律：滲流世界的歐姆定律

地下水不像河水那樣奔流，而是以極緩慢的速度在土壤顆粒間的孔隙裡滲透。描述這種流動的根本定律，是 1856 年法國工程師達西（Henry Darcy）在研究迪戎市供水砂濾池時提出的達西定律（Darcy's law）：

$$Q = -K A \frac{dh}{dl}$$

其中 $Q$ 是滲流量、$A$ 是垂直於流向的斷面積、$\dfrac{dh}{dl}$ 是水力梯度（hydraulic gradient，即水頭沿流向的變化率），而 $K$ 是水力傳導度（hydraulic conductivity，又稱滲透係數）。負號代表水從高水頭流向低水頭。

這條式子的結構，和電學裡的歐姆定律 $I = \dfrac{V}{R}$ 驚人地相似：水頭差就像電壓、滲流量就像電流、而 $1/K$ 就像電阻。$K$ 反映介質「讓水通過的難易度」——乾淨粗砂的 $K$ 可達 $10^{-3}\text{ m/s}$，而緊密的黏土卻可能低到 $10^{-9}\text{ m/s}$，相差百萬倍。這也解釋了為什麼黏土層能當作天然的「不透水底」，把地下水封在特定地層裡。

這裡要小心一個觀念陷阱。達西定律裡的 $\dfrac{Q}{A}$ 稱為達西流速（Darcy velocity）或比流量 $q$，但它不是水分子真正的移動速度。因為水只能走孔隙、走不了固體顆粒，真實的孔隙流速（seepage velocity）$v_s$ 要把孔隙率 $\phi$（porosity，孔隙體積佔比）考慮進來：

$$v_s = \frac{q}{\phi} = \frac{Q}{\phi A}$$

由於 $\phi$ 通常只有 $0.2$ 到 $0.4$，真實流速比達西流速快上好幾倍。當我們追蹤一個污染物在地下水中要多久抵達下游水井時，該用的是孔隙流速，不是達西流速——這個區別在地下水污染整治裡攸關成敗。

含水層、水井與洩降漏斗

地下水儲存在含水層（aquifer）裡。依上方是否被不透水層封住，分成兩類：

• 自由含水層（unconfined aquifer）：頂部是自由地下水面，補注直接來自上方入滲，水位就是地下水面本身。
• 受壓含水層（confined aquifer）：上下都被不透水層夾住，水承受著超過大氣壓的壓力。若鑿井打穿頂部封層，水會自動上升到一個假想高度，這個高度連成的面叫測壓水位面（piezometric surface）。當測壓水位高過地面時，水甚至會自己噴出來——這就是「自流井」（artesian well）的由來。

當我們開始抽水，井周圍的水位（或測壓水位）會下降，形成一個以井為中心、向四周漸淺的洩降漏斗（cone of depression）。穩定狀態下，受壓含水層的徑向滲流可由蒂姆方程式（Thiem equation）描述：

$$Q = \frac{2\pi T (h_2 - h_1)}{\ln(r_2 / r_1)}$$

其中 $T = Kb$ 是導水度（transmissivity，$b$ 為含水層厚度），$h_1$、$h_2$ 分別是離井距離 $r_1$、$r_2$ 處的水頭。這條式子的威力在於：只要在抽水井附近設兩口觀測井、量測它們的水位與抽水量，就能反推出含水層的導水度 $T$——這正是「抽水試驗」（pumping test）的理論基礎，是水文地質調查最常用的手段。

至於非穩態、剛開始抽水那段時間的水位變化，則由更複雜的泰斯方程式（Theis equation）描述，它引入了貯水係數（storativity）$S$ 來刻畫含水層釋水的快慢。這兩個參數 $T$ 與 $S$，就是描述任何含水層的兩張身分證。

看一個例子

讓我們用蒂姆方程式做一次完整的抽水試驗分析。

問題：某受壓含水層厚度 $b = 25\text{ m}$。在一口井以定流量 $Q = 0.03\text{ m}^3/\text{s}$ 抽水並達穩定後，量測到距井 $r_1 = 10\text{ m}$ 的觀測井水頭為 $h_1 = 48.0\text{ m}$，距井 $r_2 = 100\text{ m}$ 的觀測井水頭為 $h_2 = 50.5\text{ m}$。試求含水層的導水度 $T$ 與水力傳導度 $K$。

第一步：整理蒂姆方程式求 $T$。

$$T = \frac{Q \, \ln(r_2/r_1)}{2\pi (h_2 - h_1)}$$

代入數值：

$$T = \frac{0.03 \times \ln(100/10)}{2\pi \times (50.5 - 48.0)}$$

先算各項：$\ln(10) = 2.303$、$h_2 - h_1 = 2.5\text{ m}$、$2\pi \times 2.5 = 15.71$。

$$T = \frac{0.03 \times 2.303}{15.71} = \frac{0.0691}{15.71} \approx 4.40 \times 10^{-3}\text{ m}^2/\text{s}$$

第二步：由 $T = Kb$ 求水力傳導度。

$$K = \frac{T}{b} = \frac{4.40 \times 10^{-3}}{25} \approx 1.76 \times 10^{-4}\text{ m/s}$$

第三步：判讀結果。 $K \approx 1.76 \times 10^{-4}\text{ m/s}$ 落在「中細砂」的典型範圍內，代表這是一層出水條件不錯的含水層。

延伸思考：注意到觀測井離井越近（$r_1 = 10\text{ m}$），水位降得越多（$h_1$ 越低）——這正是洩降漏斗「越靠近井越陡」的特徵。如果今天在這口井旁邊 $50\text{ m}$ 處又鑿一口井同時抽水，兩個洩降漏斗會疊加，使中間區域的水位降幅大於任一口井單獨抽水之和。這種「井群干擾」（well interference）正是濱海地區密集抽水導致海水入侵的隱形推手。

抽水的代價：地層下陷與有效應力

現在回到開頭的高鐵問題。為什麼抽地下水會讓地面沉下去？關鍵在於土壤力學裡最重要的一條原理——有效應力原理（principle of effective stress），由現代土壤力學之父太沙基（Karl Terzaghi）提出：

$$\sigma' = \sigma - u$$

其中 $\sigma$ 是土壤承受的總應力（total stress，來自上覆土層與水的重量）、$u$ 是孔隙水壓（pore water pressure，孔隙中的水所撐起的壓力），而 $\sigma'$ 是有效應力——真正由土壤顆粒骨架承擔、決定土壤強度與變形的那部分應力。

把這三者連起來想：當你抽走地下水，孔隙水壓 $u$ 下降；但上覆土層重量造成的總應力 $\sigma$ 不變。由 $\sigma' = \sigma - u$ 可知，$u$ 一降，有效應力 $\sigma'$ 就上升。原本由水「分擔」的那部分重量，現在全壓到土壤顆粒骨架上，迫使顆粒重新排列、孔隙被擠扁——這就是壓密（consolidation）。對於含黏土的細顆粒地層，這種壓密往往是不可逆的：就算日後地下水位回升，被壓扁的黏土層也彈不回原狀。雲林、彰化沿海持續抽水灌溉與養殖，正是把這層「水墊」一點一滴抽掉，讓地表年復一年地永久下沉。

這也說明了為什麼地層下陷如此難解。它不是把抽出的水補回去就能復原的工程問題，而是一場土壤骨架的不可逆形變。因此防治策略只能是「源頭管制」——減少抽取、改用地面水或海水淡化、推動水資源循環，把孔隙水壓守住，別讓有效應力繼續攀升。

把水庫也算進來：水資源系統的調度

水利工程的尺度，可以從一口井放大到一整座水庫、甚至整個流域的水資源調配。這裡的核心工具，是質量平衡（mass balance）——其實就是入門篇連續方程式在「儲蓄體」尺度上的版本。對一座水庫，任一時段的蓄水量變化等於入流減出流：

$$\frac{dV}{dt} = I(t) - O(t)$$

寫成離散的時段差分式（水庫操作最實用的形式）：

$$V_{t+1} = V_t + I_t - O_t - E_t - L_t$$

其中 $V$ 是蓄水量、$I$ 是入流、$O$ 是放流（含供水、發電、生態基流）、$E$ 是水面蒸發、$L$ 是滲漏。水庫操作員每天都在解這道帳：上游來了多少水、下游要放多少、還要為下一場旱或洪預留多少空間。

這套質量平衡延伸出兩個關鍵的水資源規劃概念：

• 可靠度（reliability）：在一段長期紀錄中，水庫能滿足需求的時段比例。台灣水庫常以「供水可靠度 95%」為目標，意思是每 20 個時段最多容許 1 次缺水。
• 歷線法／序列模擬（sequent peak / simulation）：給定一串歷史或合成的入流序列與需水序列，反覆套用上面的質量平衡式，找出「要多大的水庫容量才能撐過最嚴重的乾旱段」。這是決定水庫設計庫容的標準方法。

台灣的水資源調度尤其艱難：降雨在時間上極度集中（豐枯比可達 9:1，雨幾乎全下在五到十月），在空間上又被中央山脈切割。河川坡陡流急、水「留不住」，使得水庫蓄豐濟枯的角色格外吃重，卻又面臨上游沖蝕導致的水庫淤積（reservoir sedimentation）逐年侵蝕有效庫容。如何在防洪、供水、發電、生態四個常常彼此衝突的目標間求取平衡，正是水資源系統分析（systems analysis）要用最佳化與模擬來回答的問題。

重點回顧

1. 達西定律（$Q = -KA\frac{dh}{dl}$） 是地下水滲流的根本，結構與歐姆定律相似；要區分「達西流速」與真正的「孔隙流速 $v_s = q/\phi$」，後者才是污染物的真實傳輸速度。
2. 含水層分自由與受壓兩類；抽水形成洩降漏斗，穩態徑向流由蒂姆方程式描述，是抽水試驗反推導水度 $T$ 的依據。
3. 有效應力原理（$\sigma' = \sigma - u$） 是理解地層下陷的鑰匙：抽水使孔隙水壓下降、有效應力上升，導致黏土層不可逆壓密——這是台灣高鐵沿線下陷的成因。
4. 水庫質量平衡（$V_{t+1} = V_t + I_t - O_t - E_t - L_t$） 把單井問題放大到流域調度，衍生出可靠度與設計庫容的規劃方法。
5. 台灣降雨時空分布極不均、河川留不住水、又有地層下陷與水庫淤積夾擊，使得「地表水—地下水聯合運用」成為無可迴避的課題。

深入探討（研究所視角）

進入研究所層級，地下水與水資源工程會從「套用解析公式」走向「偏微分方程的數值模擬」與「跨系統的最佳化決策」。

從達西定律到地下水流控制方程式。 蒂姆與泰斯方程式都只適用於理想化的徑向流，真實含水層是非均質、各向異性、邊界複雜的三維場域。完整描述地下水流動的，是把達西定律代入質量守恆後得到的地下水流方程式（groundwater flow equation），一條拋物型偏微分方程：

$$\frac{\partial}{\partial x}\!\left(K_x \frac{\partial h}{\partial x}\right) + \frac{\partial}{\partial y}\!\left(K_y \frac{\partial h}{\partial y}\right) + \frac{\partial}{\partial z}\!\left(K_z \frac{\partial h}{\partial z}\right) = S_s \frac{\partial h}{\partial t} - W$$

其中 $S_s$ 是比貯水量、$W$ 是源匯項（抽水或補注）。這條方程式沒有解析解，業界標準工具 MODFLOW（美國地質調查所開發）用有限差分法數值求解它，能模擬整個盆地的水位變化、井群干擾與海水入侵。

溶質傳輸與地下水污染。 當地下水帶著污染物移動，問題就從「水往哪流」升級為「污染物往哪擴散」。這需要在流場之上再疊一條對流—延散方程式（advection–dispersion equation），同時刻畫隨水流移動的對流項、因孔隙路徑分歧造成的延散項，以及吸附、降解等化學反應。台灣的農業硝酸鹽污染、加油站底下的油品洩漏整治，背後都是這套理論在支撐。

海水入侵與聯合運用最佳化。 濱海地區過度抽取淡水會破壞淡水—海水的密度平衡，使鹹水楔形向內陸推進，污染水井——這正是雲嘉南沿海面臨的威脅。研究前沿一方面用變密度流模式（如 SEAWAT）模擬鹹淡水界面的移動，另一方面把「抽多少、何時抽、哪裡抽」設計成一個最佳化問題：在「滿足供水需求」與「界面不可越過某警戒線」的限制下，求總成本最低的抽水排程。再結合地面水庫的調度，就形成地表水—地下水聯合運用（conjunctive use）——豐水期多用河水、把地下水當「天然水庫」存著，枯水期再有限度地動用。

氣候變遷下的補注不確定性。 地下水的天然補注高度依賴降雨入滲，而氣候變遷正讓台灣的降雨更集中於少數強降雨事件——雨來得又急又猛，反而更容易化為地表逕流直接入海，能慢慢入滲補注含水層的比例下降。這意味著未來「可永續抽取量」（safe yield）本身就是個會縮水的移動標靶，如何在不確定的補注下訂定動態的抽水管制，是水資源永續經營最前沿、也最貼近台灣生存處境的活躍課題。

值得延伸的思考：當你下次擰開水龍頭，或看到一片稻田正用抽水機灌溉，不妨想像地底那層看不見的含水層——它的水位、它的有效應力、它承載的整片土地，正隨著我們每一次的抽取悄悄改變。理解地表之下的水如何流動、又會付出什麼代價，或許是這座海島在缺水與下陷雙重夾擊下，最該補上的一堂課。