目前使用comsol實現水力壓裂的方法主要是相場法與連續介質損傷方法，相場法的實現比較復雜，不過一些學者已經把模型代碼部分開源，幫助我們學習。連續介質損傷方法發展的比較早，其中國產軟件RFPA在這方面做的比較好。目前線彈性損脆性或者軟化模型使用的比較多，對于頁巖、花崗巖水力壓裂一般使用脆性損傷模型。而對于煤這種軟巖，脆性模型有時候并不適用。基于煤破壞峰后軟化行為，軟化損傷模型比脆性模型更適用煤。相場法模型應用在彈性模量與強度比較低的巖石壓裂過程中，很容易出現模型不收斂現象。相場法主要用在彈模比較大的且以張拉破壞為主的巖石壓裂過程中，對于軟煤可能存在失效的問題。

RFPA比較適用于脆性巖石的壓裂或者破壞，模擬出來的效果也比較好，但是應用在煤的壓裂時，形成的裂縫很寬，并不能很好的反映壓裂效果。我目前借助使用比較多的COMSOL with Matlab平臺，初步實現了實驗室和現場中裂隙煤體中復雜裂縫擴展的模擬。模型中很大的問題，也是收斂問題，主要的參數與方程來自與公開發表的文獻。該模型使用的主要方程是線彈性軟化損傷方程與裂隙本構方程。水力裂縫與天然裂縫之間的相互作用，是模型的難點。comsol中的裂隙流模塊，可以實現裂隙中水流動。在5.6之前的版本中，固體力學模塊中有彈性薄層接口，這個接口可以自定義裂隙的本構方程。基于裂隙的本構模型，可以獲得裂隙表面的法向應力與剪切應力，從而實現裂隙的閉合與張開，具體方程可以參考Qinghua Lei在IJRMMS上發表的論文。使用零厚度的線段或者平面來代替裂隙，煤巖的損傷主要發生在基質中，天然裂隙或其他節理不會出現損傷。使用矩形或者很薄的長方體表征裂隙，可以設置裂隙的強度參數和根據破壞準則判斷破壞類型。不過，使用成百上千的矩形或者長方體的話，網格單元數量比較多，對計算機配置有較高的要求。COMSOL中比較容易生成離散裂隙網絡（DFN)，模型計算量會小一些。

下面幾幅圖是實驗室、現場水力壓裂裂縫擴展的效果展示圖。

Qinghua Lei的文獻題目：

Modelling fluid injection-induced fracture activation, damage growth, seismicity occurrence and connectivity change in naturally fractured rocks