房子蓋好那天才是開始：為什麼有些大樓會慢慢「沉」十年

房子蓋好那天才是開始：為什麼有些大樓會慢慢「沉」十年

入門篇談的是「土壤會不會破壞」——基礎下的土被剪壞、地基瞬間崩塌。但在台灣的軟弱沖積平原上，工程師更常半夜睡不著覺的，其實是另一個問題：土壤明明沒有破壞，建築物卻在落成後的十年、二十年間，一公分一公分地往下沉，而且不同位置沉得不一樣。雲林、彰化的地層下陷讓高鐵橋墩需要長期監測，台北盆地的軟弱黏土讓不少老建築門框歪斜、磁磚龜裂——這些都不是「強度不足」造成的破壞，而是時間（time）與水（water）共同主導的緩慢變形。

這篇進階文章要處理的，是入門篇刻意略過的另一半土壤力學：變形與時間。我們不再問「土能撐多重」，而是問「在撐得住的前提下，它會沉多少、沉多久、沉得均不均勻」。這需要把有效應力原理從靜態的力平衡，推進到動態的、隨時間演化的擴散方程式。同時，我們也會把入門篇只點到名的液化定量評估，真正算一遍。

沉陷的三個成分：彈性、主壓密、次壓密

當基礎把載重加到飽和黏土上，總沉陷量 $S_t$ 由三個物理機制疊加而成：

$$ S_t = S_e + S_c + S_s $$

第一項 $S_e$ 是彈性沉陷（immediate / elastic settlement），載重一加上去土體立刻發生的瞬時壓縮，主要來自顆粒骨架的彈性變形，與水無關。

第二項 $S_c$ 是主壓密沉陷（primary consolidation settlement），這是飽和黏土的主角。載重加上去的瞬間，因為水來不及排出，幾乎全部由孔隙水承受，產生超額孔隙水壓（excess pore water pressure）$\Delta u$。隨著時間推移，水慢慢從低滲透性的黏土中被擠出去，$\Delta u$ 逐漸消散，原本由水承受的應力一點一點轉移到顆粒骨架上——有效應力 $\sigma'$ 上升，土體隨之壓縮。這個過程在砂土中幾秒內完成，在黏土中卻可能拖上數年甚至數十年。

第三項 $S_s$ 是次壓密沉陷（secondary consolidation / creep），發生在超額孔隙水壓已經完全消散之後，土骨架在固定有效應力下仍持續的潛變變形。對有機質土、軟黏土特別明顯。

入門篇談的破壞是「瞬間、與強度有關」；這裡談的壓密是「漫長、與滲透性和時間有關」。同一塊飽和黏土，兩種行為由不同方程式主宰。

壓密參數：壓縮曲線怎麼讀

要量化主壓密沉陷，核心工具是壓密試驗（oedometer test）得到的 $e$–$\log \sigma'$ 曲線：把試體在側向受限的鋼環裡逐級加載，記錄每級有效應力下穩定後的孔隙比 $e$。畫在半對數座標上，正常壓密黏土會呈現一條近似直線，其斜率定義為壓縮指數（compression index）：

$$ C_c = \frac{\Delta e}{\Delta (\log \sigma')} = \frac{e_0 - e_1}{\log(\sigma'_1 / \sigma'_0)} $$

$C_c$ 越大，土壤被壓縮的潛力越大、越「軟」。台灣西部沖積黏土的 $C_c$ 常落在 0.3 到 0.6 之間，屬於高壓縮性。

這裡有一個關鍵概念：前期最大壓密應力（preconsolidation pressure）$\sigma'_p$。它是這塊土在歷史上曾經承受過的最大有效應力。若目前有效應力 $\sigma'_0$ 等於 $\sigma'_p$，稱為正常壓密（normally consolidated, NC）；若 $\sigma'_0 < \sigma'_p$（土曾被更大的應力壓過，例如冰河時期的覆冰或後來被侵蝕掉的上覆土層），稱為過壓密（overconsolidated, OC）。過壓密土在 $\sigma'_p$ 以內的加載只沿較平緩的回脹線（re-compression line, 斜率 $C_r$）走，沉陷量小得多；一旦載重把有效應力推過 $\sigma'_p$，才跌落到陡峭的處女壓縮線（virgin compression line, 斜率 $C_c$），沉陷量急遽放大。判斷一塊地是 NC 還是 OC，往往決定了沉陷預測差幾倍。

正常壓密黏土的主壓密沉陷量為：

$$ S_c = \frac{C_c \, H}{1 + e_0} \log \frac{\sigma'_0 + \Delta \sigma'}{\sigma'_0} $$

其中 $H$ 是可壓縮黏土層厚度，$e_0$ 是初始孔隙比，$\sigma'_0$ 是該層中點的初始有效應力，$\Delta \sigma'$ 是基礎載重在該層中點引起的應力增量。

Terzaghi 一維壓密理論：沉陷需要多久

知道「最終會沉多少」還不夠，工程上更需要知道「任一時刻沉了多少」——這牽涉到水從黏土中擴散排出的速率。Terzaghi 在 1925 年把這個過程寫成一條擴散型偏微分方程，這是古典土壤力學最漂亮的成果之一：

$$ \frac{\partial u}{\partial t} = c_v \, \frac{\partial^2 u}{\partial z^2} $$

其中 $u$ 是超額孔隙水壓，$z$ 是深度，$t$ 是時間。控制整個速率的是壓密係數（coefficient of consolidation）：

$$ c_v = \frac{k}{m_v \, \gamma_w} $$

$k$ 是滲透係數（permeability，黏土極小，約 $10^{-9}\ \mathrm{m/s}$ 量級），$m_v$ 是體積壓縮係數（coefficient of volume compressibility）。這條式子告訴我們：滲透性越低、土越可壓縮，水越難排出，壓密越慢。

這個擴散方程的解，用兩個無因次量整理最為簡潔。定義時間因數（time factor）：

$$ T_v = \frac{c_v \, t}{H_{dr}^2} $$

$H_{dr}$ 是最長排水路徑（drainage path）——若黏土層上下都能排水（雙面排水），$H_{dr}$ 等於層厚的一半；若只有單面排水，$H_{dr}$ 等於整個層厚。注意 $H_{dr}$ 是平方出現的：排水路徑加倍，所需時間變成四倍。這正是工程上插設「排水砂樁／塑膠排水帶（PVD）」加速壓密的理論依據——人為製造大量短排水路徑，把原本要十年的沉陷壓縮到幾個月內完成，讓建築物在交屋前就「沉到位」。

平均壓密度（degree of consolidation）$U$ 與 $T_v$ 有近似關係，常用的工程公式為：

$$ T_v \approx \frac{\pi}{4} U^2 \quad (U \le 60\%), \qquad T_v \approx 1.781 - 0.933 \log(100 - U\%) \quad (U > 60\%) $$

看一個例子：軟黏土上的填土要等多久

某海埔新生地工程，在 6 公尺厚的飽和正常壓密黏土上填築路堤。黏土層上方為透水砂層、下方為不透水基岩，屬單面排水。已知 $C_c = 0.45$、初始孔隙比 $e_0 = 1.1$、黏土層中點初始有效應力 $\sigma'_0 = 50\ \mathrm{kN/m^2}$、填土在中點引起的應力增量 $\Delta\sigma' = 80\ \mathrm{kN/m^2}$、壓密係數 $c_v = 2.0\ \mathrm{m^2/yr}$。

第一步，算最終主壓密沉陷量。

$$ S_c = \frac{0.45 \times 6}{1 + 1.1} \log \frac{50 + 80}{50} = \frac{2.7}{2.1} \log(2.6) $$

$$ S_c = 1.286 \times 0.415 \approx 0.534\ \mathrm{m} $$

也就是說，這塊地最終會沉下約 53 公分——非常可觀，若不預先處理，路堤完工後將持續下陷、開裂。

第二步，算沉到 90% 需要多久。

單面排水，所以 $H_{dr} = 6\ \mathrm{m}$。$U = 90\%$ 屬 $U > 60\%$ 範圍：

$$ T_v = 1.781 - 0.933 \log(100 - 90) = 1.781 - 0.933 \times 1 = 0.848 $$

由 $T_v = c_v t / H_{dr}^2$ 反解時間：

$$ t = \frac{T_v \, H_{dr}^2}{c_v} = \frac{0.848 \times 6^2}{2.0} = \frac{0.848 \times 36}{2.0} \approx 15.3\ \mathrm{yr} $$

等 15 年才完成 90% 沉陷，顯然無法接受。這正是預壓（preloading）加堆載排水帶（PVD）工法登場的時機：若插設 PVD 把排水改成徑向、有效排水路徑縮短到 1 公尺等級，$H_{dr}^2$ 從 36 降到接近 1，同樣的 90% 沉陷可壓縮到數個月內完成。台灣高鐵沿線、桃園與台中的許多軟弱地盤大型填方，都採用這套「先讓它沉夠、再蓋上去」的策略。

把液化算到底：CSR 對 CRR

入門篇提到飽和砂土在地震中會因孔隙水壓飆升而液化，並點到用安全係數 $FS_L = CRR / CSR$ 評估。進階一點，我們把這個比值真正算一遍，因為它是台灣每一份大型建築地質調查報告都要做的工作。

循環應力比（cyclic stress ratio, CSR） 代表地震施加給該深度土壤的「需求」，由 Seed–Idriss 簡化式給出：

$$ CSR = 0.65 \, \frac{a_{max}}{g} \, \frac{\sigma_{v0}}{\sigma'_{v0}} \, r_d $$

其中 $a_{max}$ 是地表最大加速度（peak ground acceleration, PGA），$\sigma_{v0}$、$\sigma'_{v0}$ 分別是該深度的總垂直應力與有效垂直應力，$r_d$ 是隨深度遞減的應力折減係數。注意 $\sigma_{v0}/\sigma'_{v0}$ 這一項：地下水位越高、有效應力越小，這個比值越大，CSR 越高——再次呼應入門篇的核心，水位是液化風險的關鍵推手。

循環阻抗比（cyclic resistance ratio, CRR） 代表土壤抵抗液化的「能力」，由標準貫入試驗修正後的 $(N_1)_{60}$ 經驗推估。當 $FS_L = CRR/CSR < 1$，判定該深度可能液化。

動手試試：判斷一個深度會不會液化

某基地地下水位在地表下 1 公尺，欲評估地表下 5 公尺處飽和砂層。已知該深度總應力 $\sigma_{v0} = 95\ \mathrm{kN/m^2}$、有效應力 $\sigma'_{v0} = 56\ \mathrm{kN/m^2}$，設計地震 $a_{max} = 0.32g$（約對應台灣多數都會區的設計地震），深度折減 $r_d = 0.97$，現地修正後 $(N_1)_{60} = 12$ 對應的 $CRR_{7.5} = 0.13$。

先算 CSR：

$$ CSR = 0.65 \times 0.32 \times \frac{95}{56} \times 0.97 = 0.65 \times 0.32 \times 1.696 \times 0.97 \approx 0.342 $$

安全係數：

$$ FS_L = \frac{CRR}{CSR} = \frac{0.13}{0.342} \approx 0.38 $$

$FS_L = 0.38$ 遠小於 1，判定此深度在設計地震下高度可能液化。實務上 $(N_1)_{60} = 12$ 的鬆砂在 0.32g 地震下確實非常危險，需採地盤改良（如夯實、灌漿、置換）或改用深基礎穿越液化層。台灣經濟部地質調查所公布的全台液化潛勢圖，背後正是把這套逐深度計算疊加成的成果。

重點回顧

• 入門篇處理「強度與破壞」，進階篇處理「變形與時間」：飽和黏土的沉陷由彈性、主壓密、次壓密三個機制疊加，主壓密是軟弱地盤的主角。
• 主壓密沉陷量由壓縮指數 $C_c$ 與應力比的對數決定；判斷土壤是正常壓密（NC）還是過壓密（OC），常使沉陷預測相差數倍。
• Terzaghi 一維壓密把孔隙水壓消散寫成擴散方程 $\partial u/\partial t = c_v\,\partial^2 u/\partial z^2$；時間因數 $T_v = c_v t / H_{dr}^2$ 中排水路徑以平方出現，是 PVD 排水帶能大幅加速沉陷的理論根據。
• 液化評估把地震需求 CSR 與土壤阻抗 CRR 相比，$FS_L = CRR/CSR < 1$ 即判定可能液化；地下水位越高、$(N_1)_{60}$ 越低，風險越大。
• 台灣軟弱沖積平原與高地下水位，使壓密沉陷與液化成為大地工程必須同時面對的兩大本土課題。

深入探討（研究所視角）

把上面的古典框架推到研究前沿，有幾條值得追的線。

第一條是壓密理論本身的修正。Terzaghi 一維理論假設小應變、線性 $e$–$\log\sigma'$ 關係、$c_v$ 為常數，但軟黏土的大變形壓密下，孔隙比改變會讓滲透係數 $k$ 隨之劇烈下降，$c_v$ 不再是常數。Gibson 等人發展的大應變壓密理論（large-strain consolidation）用 Lagrangian 座標重寫方程，對淤泥、疏濬土、尾礦壩這類超軟材料才能準確預測。此外，Terzaghi 理論把主壓密與次壓密截然分開，但 Bjerrum 的「等時間線（isotache）」模型主張兩者其實是同一潛變過程的連續展現，這在預測長期沉陷時差異顯著。

第二條是把壓密與強度統一在臨界狀態框架下。入門篇提到的劍橋模型（Cam-Clay）正是這個統一的橋梁：它在 $p'$–$q$–$v$（平均有效應力、偏差應力、比容）空間中同時描述壓密的體積變化與剪切的強度變化，用同一組參數（$\lambda$、$\kappa$、$M$）刻畫。現代有限元素軟體（如 PLAXIS）內建的硬化土模型（Hardening Soil）、軟土潛變模型（Soft Soil Creep）都是 Cam-Clay 思想的工程化延伸，能在同一次分析裡算出開挖、加載、消散的完整時程。

第三條與台灣最切身的是液化評估方法的演進與地層下陷的耦合。Seed–Idriss 簡化法是 1970 年代的經驗框架，後續 Idriss–Boulanger（2008、2014）大幅修訂了 $r_d$、震矩規模修正係數 $MSF$、與上覆應力修正 $K_\sigma$，並把 CPT 的 $q_{c1N}$ 納為更可靠的判據。更前沿的方向，是用有效應力動態分析（effective-stress dynamic analysis）直接在時域中模擬地震過程裡孔隙水壓的累積與消散，而非用單一安全係數判斷有無——這能進一步預測液化後的「側向流動（lateral spreading）」與「沉陷量」，這正是 1999 年集集地震沿濁水溪、2016 年美濃地震在台南安南區造成嚴重災損的機制。再加上雲彰地區因超抽地下水導致的長期壓密下陷，使有效應力、壓密、與液化這三條看似獨立的線，在台灣的地下水脈絡裡其實緊緊纏在一起——這也是為什麼大地工程在這座島上，永遠不只是力學，更是一門關於水與時間的學問。