【引言】

在頁巖油開采的競技場上，水力壓裂技術曾被譽為“解鎖地下黑金的鑰匙”，但隨之而來的水資源消耗、化學污染和低效裂縫預測等問題，正讓這把鑰匙逐漸生銹。當全球能源行業將目光投向更清潔的超臨界CO?壓裂技術時，一個更棘手的難題浮出水面：如何馴服這種介于氣液之間的“暴躁流體”，精準預判它在千米地層下的裂巖軌跡？

傳統實驗手段如同“盲人摸象”——物理模擬成本高昂，現場試錯風險巨大，而經典數值模型又難以刻畫CO?與頁巖間復雜的相變交互。直到相場法（Phase-Field Method）的出現，這場博弈迎來了轉機。這項起源于材料科學的數學工具，正將超臨界CO?轉化為可被方程描述的“數字流體”，在虛擬空間中重構裂縫生長的每一個細節：從CO?分子穿透巖石孔隙的微觀動力學，到宏觀裂縫網絡的混沌分形演化，原本不可見的流體暴力被解構成萬億次計算的優雅舞蹈。

當算法的精度突破物理實驗的邊界，頁巖壓裂正在從“經驗驅動”邁向“預測驅動”的新紀元。這場由相場法引領的仿真革命，或將重新定義非常規油氣開采的底層邏輯——用數字孿生代替盲目試錯，用計算預見性取代經驗不確定性。而我們，正站在這場技術范式轉移的臨界點上。

【超臨界CO2壓裂模擬存在難點】

難點一：CO?的“流體變身記”——超臨界態下的物性捉迷藏

超臨界CO?既非氣體也非液體，密度接近液體，黏度類似氣體，擴散性極強。這種“雙重人格”導致它在頁巖孔隙中的流動規律難以捉摸。傳統本構方程（如Forchheimer方程）在水力壓裂中表現穩定，但面對超臨界CO?時，其滲透率-壓力關系會因微小的溫壓波動（±5℃或±0.5MPa）而發生劇變，就像在仿真中埋下無數個“數值地雷”。更棘手的是，CO?與頁巖有機質接觸時可能引發萃取效應（如瀝青質溶解），進一步改變巖石力學參數——流體在改造巖石，巖石也在反向馴化流體。

難點二：裂縫演化的“混沌預言”——多物理場耦合陷阱

CO?壓裂并非單純的流體驅動裂縫問題，而是熱-流-固-化（THMC）四場耦合的混沌系統：

*熱場：CO?注入時的焦耳-湯姆遜效應會引發近井地帶溫度驟降（可能達30℃），導致巖石脆性增強；

*化學場：CO?-地層水-頁巖的化學反應（如碳酸鹽化）會動態改變孔隙結構和裂縫表面能；

*流固耦合：超臨界CO?的低黏度特性使其更易滲入納米級孔隙，誘發“裂縫分支”（多級次生裂縫無序增殖）。

難點三：相場法的“精度圍城”——從方程到現實的最后一公里 盡管相場法通過序參量（Phase-Field Parameter）優雅地統一了裂縫形核、擴展與交匯的數學描述，但其在CO?壓裂中的落地仍面臨三重圍剿：

*參數：CO?-頁巖界面能、裂縫表面耗散率等關鍵參數在高溫高壓下極難通過實驗獲取，許多模型被迫使用“借來的數據”（如借用甲烷或水的數據）；

*尺度：納米孔隙（~10nm）中的CO?分子吸附效應與宏觀裂縫（~10m）的擴展需跨10個數量級建模，常規網格剖分技術如同“用漁網撈病毒”；

*驗證困局：地下壓裂過程無法直接觀測，即便仿真結果與微地震監測數據吻合，也可能只是“錯誤的勝利”（不同機理導致相似表象）。

【采用COMSOL模擬實現】

采用COMSOL實現相場法模擬裂縫擴展過程，考慮以下幾個因素：

（1）超臨界CO2物性隨壓力變化，包括流體粘度和密度；

（2）超臨界CO2流動能力強，很容易進入孔縫之中，改變巖體有效應力，減弱巖體的剪切強度，因此，需要考慮拉伸和剪切強度變化，建立相場驅動方程，從而實現多分支裂縫擴展形態；

（3）進一步，可以分析各類因素的綜合影響，包括：天然裂縫、地應力、紋層以及多簇裂縫應力干擾作用等。

采用滑溜水壓裂，求解獲得簡單裂縫形態

采用超臨界CO2壓裂，求解獲得多分支裂縫形態

https://mp.weixin.qq.com/s/i42LhN_9amAAOlMHBc7lTA?