《Cleaner Engineering and Technology》：The Effect of Stress History on the Hydraulic Properties and Geomechanical Behavior of Caprock During CO
2 Breakthrough
編輯推薦：
本文聚焦地質碳封存（GCS）中蓋層完整性評估的關鍵挑戰。為解決應力歷史對蓋層水力-力學耦合行為影響不清的問題，研究人員系統研究了不同超固結比（OCR）對高嶺石泥巖在CO2注入下的水力特性和變形響應。結果表明，OCR增加會降低孔隙度和滲透率，但提高CO2突破壓力，且不同流體（鹽水/CO2）注入引發截然不同的變形模式。該研究揭示了應力歷史與多相流耦合對蓋層密封能力的影響，為GCS的安全設計提供了重要實驗依據。
隨著全球氣候變化問題日益嚴峻，如何有效減少大氣中的二氧化碳（CO2）成為科學家們亟待破解的難題。地質碳封存（Geologic Carbon Storage, GCS）作為一種頗具前景的規模化解決方案，其核心是將捕獲的CO2注入深層地下巖層，例如咸水層或廢棄油氣藏，以期實現長達地質時間尺度的安全封存。然而，這個看似完美的方案背后隱藏著一個關鍵的“守門人”問題：上覆的蓋層（caprock）是否足夠堅固，能永久阻擋具有浮力的CO2向上逃逸？
蓋層，通常是低滲透性的泥巖，如同一個巨大的地下密封蓋。當CO2在儲層中聚集并向上運移時，最終會抵達蓋層底部，并對其施加壓力。一旦這個壓力超過了蓋層巖石孔隙的毛細管進入壓力（capillary entry pressure），CO2就可能找到“突破口”，開始滲透，從而威脅整個封存項目的安全。因此，準確評估蓋層的“密封”能力至關重要。
以往的研究多側重于測量蓋層的突破壓力，評估其毛細管密封性能。但現實情況更為復雜，因為CO2的注入不僅是一個流體流動問題，還會引發一系列地質力學（geomechanical）響應。注入的流體會改變巖石內部的孔隙壓力，進而改變其有效應力狀態，可能導致巖石變形甚至破壞。尤其對于富含粘土的泥巖蓋層，其力學行為強烈依賴于其所經歷的應力歷史——就像黏土被反復揉捏后，其軟硬和孔隙結構會發生變化一樣。這種應力歷史通常用超固結比（Overconsolidation Ratio, OCR）來量化。OCR等于1表示巖石當前承受的應力就是其歷史上承受過的最大應力（正常固結狀態），而OCR大于1則表示巖石曾經承受過比現在更大的應力（超固結狀態）。顯然，一個歷經滄桑、被深埋后又經歷抬升剝蝕的“老”蓋層（高OCR），與一個沉積不久、相對“年輕”的蓋層（低OCR），其內部結構和力學特性必然不同。然而，關于不同的應力歷史如何具體影響蓋層在CO2注入和突破過程中的水力特性與變形行為，人們所知甚少。這種知識空白限制了我們對蓋層長期完整性的準確預測，也增加了GCS工程的不確定性。
為了填補這一空白，一篇發表在《Cleaner Engineering and Technology》上的研究，由Sangyeong Park、Hyeontae Park、Hangseok Choi和Kiseok Kim（通訊作者）共同完成，進行了一項系統而深入的實驗研究。他們的目標很明確：揭示應力歷史（以OCR表征）對泥巖蓋層在CO2突破過程中的水力特性和地質力學行為的具體影響。研究人員想知道，經歷過不同“命運”的巖石，在面對CO2的“進攻”時，會如何應對？它們的“防線”（滲透率和突破壓力）有多堅固？內部結構又會發生怎樣的“形變”？
為了回答這些問題，研究人員開展了一項精密的實驗。他們采用定制的高壓固結系統，以高嶺石和鹽水制備了成分均一的人工合成泥巖試樣，以排除天然巖石礦物不均一性的干擾。通過控制加載和卸載路徑，他們成功制備了具有三種不同應力歷史（OCR分別為1、2和4）的試樣，并讓它們最終都在14.50兆帕的垂直有效應力下達到平衡，以模擬約1.45公里埋深的地下條件。在固結完成后，研究人員對這些試樣先后進行了鹽水滲透率測試和CO2突破壓力測試。在整個流體注入過程中，他們利用激光位移傳感器實時、高精度地監測試樣的軸向變形（反映為孔隙比e的變化），從而實現了水力過程與力學變形的耦合監測。CO2突破壓力采用逐步升壓法（Step-by-Step, SBS）進行測定，該方法通過逐步增加CO2注入壓力并監測出口端體積變化，來精確確定突破發生的壓力區間。
主要研究結果如下：
1. 應力相關的壓縮性和時間依賴性響應
通過固結實驗，研究人員量化了材料的壓縮指數（Cc）、回彈/再壓縮指數（Cr）和次壓縮指數（Cα）。結果清晰顯示了泥巖的彈塑性行為：在初始加載（ virgin loading）階段壓縮性大（Cc= 0.21），而在卸載/再加載階段壓縮性顯著變小（Cr= 0.01）。OCR越高，在相同最終應力下的孔隙比越低，表明材料更致密。次壓縮指數Cα隨著垂直應力的增加而減小，表明在高應力下，時間依賴性變形（蠕變）的潛力降低。
2. 鹽水滲透率和CO2突破壓力
實驗結果清晰地展示了應力歷史對水力特性的強大影響：
- •
孔隙度與滲透率：隨著OCR從1增加到4，試樣的孔隙度從0.34顯著降低至0.20。相應的，其固有滲透率（intrinsic permeability）也從5.0×10-18m2下降到1.5×10-18m2。值得注意的是，滲透率的下降在OCR從1到2時最為劇烈，之后趨于平緩。
- •
CO2突破壓力：與滲透率趨勢相反，CO2突破壓力隨OCR增加而單調上升。對于OCR=1、2和4的試樣，其突破壓力分別為4.5、5.5和6.9兆帕。這表明，經歷過更大應力歷史的蓋層，其抵抗CO2突破的毛細管密封能力更強。
研究人員將這種現象歸因于孔隙結構的演化。在加載（固結）過程中，較大的孔隙（宏觀孔隙）優先被壓縮，這些孔隙正是流體（尤其是非潤濕相CO2）運移的主要通道。因此，OCR的增加（意味著經歷了更高的歷史應力）會導致大孔隙減少，從而顯著降低滲透率。同時，更小、更不連通的孔隙網絡增加了CO2突破所需的毛細管力，從而提高了突破壓力。汞侵入孔隙度測定法（MIP）的微觀分析支持了這一解釋，顯示OCR=4的試樣相比OCR=1的試樣，其孔隙喉道半徑分布更偏向小尺寸端，大喉道的貢獻減少。
3. 鹽水與CO2注入過程中的地質力學響應
這是本研究揭示耦合行為的關鍵部分：
- •
鹽水注入：鹽水注入提升了試樣內的平均孔隙壓力，降低了作用在巖石骨架上的有效應力，導致試樣發生彈性膨脹，孔隙比增加。這種膨脹沿著固結曲線中的卸載路徑（與Cr相關）進行。
- •
CO2注入：CO2注入引發了截然不同的力學響應，且突破前后行為不同。
- •
突破前：在CO2突破發生之前，試樣表現出輕微的壓縮（孔隙比減小），其變形速率遠低于單純的次壓縮（蠕變）預期值。這表明，CO2的注入壓力并未完全轉化為孔隙壓力作用于骨架，可能涉及CO2在鹽水中的部分溶解、氣-液界面的形成以及相關的化學-力學耦合過程，這些耦合效應抑制了時間依賴性變形。
- •
突破后：一旦CO2突破發生，試樣開始膨脹（孔隙比增加），但膨脹的幅度小于在更低注入壓力下鹽水注入所引起的膨脹。此外，OCR越高的試樣，在突破前和突破后的變形速率都越低，表現出更小的壓縮性。
研究結論與討論
本研究通過精密的耦合實驗，系統闡明了應力歷史對泥巖蓋層在CO2注入場景下水力-力學行為的關鍵控制作用。主要結論可歸納為三點：
- 1.
應力歷史重塑孔隙結構，對水力特性產生非對稱影響：OCR增加通過漸進式孔隙坍塌，顯著降低蓋層的孔隙度和固有滲透率，尤其是滲透率在OCR初步增加時下降最為明顯。與之協同的是，CO2突破壓力隨OCR增加而持續升高，表明超固結蓋層具有更強的毛細管密封能力。這源于高應力歷史優先壓縮了控制流體運移的大孔隙和連通喉道。
- 2.
流體類型與應力歷史共同控制地質力學響應：鹽水和CO2注入會引發截然不同的變形機制。鹽水注入通過提高孔隙壓力、降低有效應力導致蓋層發生彈性膨脹。而CO2作為非潤濕相，其注入過程伴隨復雜的多相流與界面效應，在突破前觀測到的變形速率異常低，揭示了毛細管壓力、部分溶解和化學相互作用等耦合過程抑制了傳統意義上的固結變形。突破后雖也引起膨脹，但機制和幅度不同于鹽水。
- 3.
高OCR蓋層表現出更強的綜合密封性能：高OCR（超固結）蓋層不僅具有更低的滲透率和更高的突破壓力，其力學變形（壓縮和膨脹）的幅度和速率也更小，表現出更高的剛度和更弱的時間依賴性。這意味著，在經歷過大應力歷史的蓋層下方進行CO2封存，可能同時享有更優的水力密封性和力學穩定性。
這項研究的意義在于，它將應力歷史這一地質背景因素定量化地納入了蓋層完整性評估框架，強調了一個完整的評估必須同時考慮毛細管突破準則和水力-力學耦合變形。研究結果挑戰了僅依靠突破壓力測試來評估密封能力的傳統做法，指出對于低OCR（正常固結）的年輕蓋層，即使其突破壓力尚可，也可能因注入引起的有效應力降低和變形而面臨更高的完整性風險。該工作為地質碳封存場址的篩選、蓋層密封能力的精細化評價以及長期安全性的數值模擬提供了關鍵的實驗數據和機理認識。未來研究需要拓展至礦物組成更復雜的天然泥巖、考慮更長的時間尺度和化學反應（如CO2-水-巖相互作用）的影響，并將實驗室尺度的發現通過建模升級至場地尺度，以更好地服務于安全、高效的碳中和地下實踐。
