裂隙巖體注漿最麻煩的不是注漿泵選型也不是漿液配比而是你根本不知道漿液會往哪條縫里鉆。等效多孔介質把裂隙“抹勻”了算出來的擴散半徑只能當個心理安慰真到了現(xiàn)場漿液經(jīng)常沿著一條窄縫跑出上百米把預算和工期一起帶走。所以當我在一個注漿類模擬項目里拿到“COMSOL三維離散裂隙注漿模型基于粘度空間衰減的漿液裂隙采用隨機分布圓盤”這個需求時反而松了一口氣——因為它終于肯承認裂隙是離散的也肯把漿液流變性從常數(shù)升級成場。這篇文章我把整個建模鏈路拆開講隨機圓盤裂隙怎么生成、粘度空間場怎么實現(xiàn)、COMSOL里物理場怎么搭、網(wǎng)格和收斂有哪些坑以及最后怎么把結果變成工程判斷依據(jù)。適合正在做COMSOL二次開發(fā)、巖體注漿模擬或裂隙巖體滲流研究的人直接照著步驟走。1. 裂隙巖體注漿模擬為什么必須走離散裂隙路線1.1 等效多孔介質模型的短板很多人一上來就用“多孔介質滲流濃度擴散”做注漿模擬覺得只要把滲透系數(shù)調大、孔隙率調高就能近似出裂隙的效果。這種思路在宏觀地下水流場里能用但在注漿問題上會嚴重失真。因為注漿的核心矛盾是“漿液優(yōu)先走哪條通道”而等效模型把通道信息平均掉了。舉個例子兩組正交裂隙和一組傾向完全一致的平行裂隙在等效連續(xù)介質里可以算出幾乎相同的滲透率張量但真實注漿行為南轅北轍正交裂隙網(wǎng)絡容易形成網(wǎng)狀擴散單組裂隙則只沿走向延伸。另一個實際問題注漿封堵需要知道漿液在某條主裂隙里的最大擴散距離等效模型只能給出一個平均濃度場給不了“某一條縫里漿液走了多遠”這種決策信息。離散裂隙網(wǎng)絡DFN的思路是徹底反過來——不把裂隙巖體看成均勻介質而是把每一條裂隙都當作獨立流體通道。注漿模擬里真正要解的就是漿液在這些通道里的運動、滯留和壓力響應。這套思路對工程決策的價值非常直接想知道單孔注漿量本質是求所有連通裂隙的累計體積想知道串漿風險本質是找某條高滲透裂隙能否貫通到臨空面。1.2 圓盤模型的空間優(yōu)勢裂隙幾何有很多種抽象方式多邊形面片、圓盤、線段網(wǎng)絡、3D格子裂隙。其中圓盤模型是DFN里最經(jīng)典也最均衡的一檔選擇。理由有三條。第一它在統(tǒng)計描述上非常干凈。一條裂隙只需要四個參數(shù)中心點位置、圓盤半徑、產狀傾向和傾角、開度。這四個參數(shù)都有成熟的概率分布可采樣野外測線數(shù)據(jù)也容易反算。第二圓盤之間的相交關系能形成完整的裂隙網(wǎng)絡連通圖。任意兩個圓盤如果空間相交就形成一條跡線這些跡線拼接起來的通道就是漿液能走的路徑。這個邏輯和真實巖體卸荷拉裂隙的幾何很接近。第三圓盤在數(shù)值實現(xiàn)里不會出現(xiàn)多邊形角點處的應力集中等效問題幾何清潔度更高。多邊形裂隙在角點附近會產生網(wǎng)格質量差、局部壓力異常等數(shù)值噪聲圓盤則沒有這種尖銳邊角收斂性明顯好一些。當然圓盤假設也有局限性比如真實裂隙面往往是不規(guī)則的粗糙度和起伏會顯著影響開度分布。在初版模型里可以先用平直圓盤把流動規(guī)律跑通粗糙度問題留到下一階段用開度隨機場去疊加。1.3 三個關鍵詞如何互相咬合這個標題里其實藏著一條很緊的邏輯鏈。三維是幾何維度——只有三維才能描述圓盤在空間里的隨機產狀二維的“線段裂隙”沒法表達傾向傾角帶來的各向異性。隨機圓盤是網(wǎng)絡拓撲——它決定了漿液的可達空間和優(yōu)先通道。粘度空間衰減則是流變維度——它把漿液從“一種液體”升級成“一種狀態(tài)隨路徑變化的液體”。三者缺一不可如果只用等效介質漿液擴散會被“平均化”粘度空間變化也失去了通道指向性如果不用粘度空間場注漿中期壓力計算會系統(tǒng)偏差因為真實漿液越走越稠繼續(xù)當常數(shù)算出的擴散半徑會偏大。這套模型能回答的具體工程問題包括單孔注入量上限、漿液最大擴散半徑的概率分布、注漿壓力隨時間抬升曲線、不同裂隙產狀組合下串漿風險的方向性。建模型之前先把這些問題寫在紙上后面每步參數(shù)選擇都會更清醒。2. 隨機圓盤裂隙網(wǎng)絡生成從統(tǒng)計參數(shù)到可計算幾何2.1 四個必須確定的統(tǒng)計控制參數(shù)生成可靠隨機圓盤裂隙網(wǎng)絡第一步是定統(tǒng)計參數(shù)不是打開軟件畫圓。常見的四個控制量參數(shù)含義常用分布/取值參考敏感性體密度 P32單位體積巖石中裂隙總面積m2/m30.12視巖體級別高半徑 R圓盤半徑對數(shù)正態(tài)分布或冪律分布均值0.55m高產狀傾向傾角Fisher分布集中度k550中開度 b裂隙張開寬度對數(shù)正態(tài)均值0.051mm極高P32是這套參數(shù)里最核心的標定量因為它可以跟現(xiàn)場壓水試驗、巖芯RQD建立換算關系。我自己做小規(guī)模注漿模擬時P32取0.5、圓盤半徑均值取2m、Fisher集中度取20得到的連通網(wǎng)絡和現(xiàn)場示蹤結果基本能對上。P32如果定到2以上圓盤數(shù)量會爆炸幾何后處理會非常痛苦建議先用小值試跑。2.2 裂隙圓盤生成算法流程圓盤生成的標準流程分四步撒點、采樣半徑、旋轉到產狀、裁剪邊界。第一步在模擬域內按泊松過程撒中心點平均點數(shù)由P32除單個圓盤平均面積得到。這里注意圓盤面積要用E(R2)不是E(R)2否則密度會偏小。第二步對每個中心點從半徑分布里采樣一個半徑值。工程上常見的是冪律分布P(R ≥ r) (r/r_min)^(-α)α通常在23。冪律的好處是能產生少數(shù)大裂隙這對連通性影響遠大于一堆小裂隙。第三步從Fisher分布采樣產狀。Fisher分布在球坐標系里繞中心方向按集中度k聚集k越大產狀越一致k趨近0時接近均勻隨機。生成時先采樣極角標準差再均勻轉轉方位角最后旋轉到設計優(yōu)勢產狀方向。第四步把所有圓盤裁剪到模擬域邊界內。COMSOL的幾何單位建議用米圓盤作為三維圓面片導入厚度可以后續(xù)用物理場代替。我常用的做法是寫一段腳本生成CSV格式的圓心坐標、半徑、傾向、傾角再通過LiveLink或直接讀入COMSOL幾何。2.3 導入COMSOL前的幾何清理直接導入上千個圓盤面COMSOL幾何會變得非常難處理光布爾操作就能卡死。我建議分成三步清理。第一步先刪除孤立的、完全不與其他圓盤相交的裂隙。它們在連通性里沒有貢獻只會增加網(wǎng)格量。判斷相交最快的方法是檢測兩圓盤圓心距離小于半徑之和再復核兩面空間相交。第二步做聯(lián)合體Form Union而不是裝配體Form Assembly。聯(lián)合體會自動計算相交跡線裂隙交叉處的幾何拓撲會被合并裝配體需要手動處理接觸對對注漿模擬這種跨裂隙流動問題非常不友好。缺點是聯(lián)合體在幾何體量大時計算慢但如果圓盤數(shù)量控制在200以內問題不大。第三步刪掉極小碎片。圓盤相交裁剪后經(jīng)常出現(xiàn)微米級的細長邊這是網(wǎng)格質量的隱形殺手。在COMSOL里設置最小幾何容差把長度小于特征尺寸十分之一的邊直接忽略或合并。2.4 隨機參數(shù)表與敏感性邏輯模型最好做成參數(shù)化驅動而不是每次手動改幾何。參數(shù)對象里至少放以下幾項P32、半徑分布參數(shù)、Fisher集中度、模擬域尺寸、注漿壓力、漿液密度、基漿粘度、后期粘度上限、粘度時間常數(shù)。實測下來對最終擴散半徑影響最大的始終是開度和P32粘度參數(shù)影響次之半徑分布形態(tài)影響居第三。所以如果計算資源緊張優(yōu)先精確標定開度別上來就擬合復雜流變參數(shù)。3. “粘度空間衰減”的物理含義與數(shù)值化處理3.1 注漿漿液的流變學身世水泥基漿液最常用的流變模型是賓漢姆流體也就是存在一個屈服應力應力超過它才開始流動。完整描述需要兩個參數(shù)塑性粘度μ和屈服應力τ0。普通水泥漿的塑性粘度可以從幾十毫帕·秒到幾百毫帕·秒隨水灰比變化劇烈屈服應力影響漿液最終的靜置穩(wěn)定性和擴散停止條件。很多教材里把漿液粘度當成常數(shù)這在短時、短距離注漿里能接受但一旦漿液在裂隙里走了幾十米、過了幾十分鐘水化反應會讓粘度明顯增大。如果仍然按初始粘度算計算出的最終擴散半徑會偏大注漿壓力的上升速率會偏慢這種系統(tǒng)性偏差在工程上很難看。3.2 粘度為什么要寫成空間分布標題里的“粘度空間衰減”我說句實話——嚴格來說大多數(shù)水泥漿液隨齡期是粘度增大的不是衰減但“剪切變稀型的漿液”在擴散過程確實會表現(xiàn)有效粘度下降。不管方向如何建模的本質是一致的粘度不能再作為定值而必須作為空間位置和時間的函數(shù)來參與達西方程求解。為什么是空間場而不是簡單時間場因為漿液在每一處位置經(jīng)歷的“運輸時間”不同?？拷{孔的漿液最早進入裂隙水化齡期最長粘度最大剛到達新生前鋒的漿液最新鮮粘度最低。如果我們只用全局時間來表達粘度等于假設全空間所有漿液同時變稠這完全錯誤。正確做法是用一個隨達西速度平流輸運的“流體年齡”或“漿液體積分數(shù)”標量場在空間上自然區(qū)分新漿和舊漿。3.3 兩種粘度空間場模型模型一年齡驅動型。適用于水泥基漿液水化增稠。定義一個標量場T表示漿液在局部駐留的時間控制方程?T/?t u·?T 1 ?·(D_T?T)右側的源項“1”表示每前進一秒漿液年齡增加一秒擴散項代表微觀混合。粘度寫成年齡的飽和函數(shù)μ μ0 (μ∞ - μ0)·(1 - exp(-T/τ))其中μ0是初始黏度μ∞是極限黏度τ是黏度增長時間常數(shù)。實測工程水泥漿τ取幾百秒到幾千秒之間視配方而定。這個模型的好處是物理意義清楚參數(shù)可以直接從旋轉粘度計實驗標定。模型二剪切變稀型。適用于聚合物類或膨潤土系漿液有效粘度隨剪切率增加而下降。局部剪切率可以從達西速度估算γ ≈ 6·|u|/b平行板縫均勻流近似。然后用Cross模型或冪律模型μ μ_min (μ_max - μ_min)/(1 (γ/γ_c)^m)這個模型無需額外傳遞變量只要在物理場里定義速度依賴的粘度表達式即可但注意它天然不穩(wěn)定高速區(qū)粘度低反過來又增大流量容易造成“短跑效應”。數(shù)值上通常要加阻尼或限制最大流速。3.4 COMSOL里的變量耦合思路兩種模型在COMSOL里的實現(xiàn)路徑略有差異。年齡驅動型需要加一個“稀物質傳遞”或“系數(shù)型偏微分方程”接口求解年齡場同時把年齡場的求解與達西方程分開做“分離式迭代”剪切變稀型只需要在變量里寫可能依賴于速度場的表達式不需要額外物理場。我建議第一版只做年齡驅動型因為它的收斂性比剪切變稀型穩(wěn)定得多而且對工程水泥漿更貼近。變量設置示例在“定義-變量”里寫mu mu0 (muInf - mu0)*(1 - exp(-Tfracture/tauVisc))其中Tfracture是年齡場。注意所有單位統(tǒng)一成國際單位制粘度用Pa·s時間用s。這里有一個很關鍵的數(shù)值問題粘度從μ0到μ∞可以跨一到兩個數(shù)量級直接參與達西滲透系數(shù)會導致滲透系數(shù)劇烈突變。處理方式是把粘度變化限制在“已注漿區(qū)域”未注漿區(qū)域保持水的粘度防止空曠區(qū)出現(xiàn)極低滲透系數(shù)。4. COMSOL模型組裝物理場、裂隙邊界與求解順序4.1 物理場選擇達西定律加傳遞方程注漿漿液在裂隙里的真實流動是縫隙流完整描述要用Navier-Stokes方程。但裂隙開度通常在0.1mm量級流速很小雷諾數(shù)遠低于臨界值粘性力主導完全可以用Darcy定律。物理學上沒有區(qū)別——平行板縫隙流的壓力-流量關系本來就是立方定律而立方定律與達西定律形式同構。所以我的物理場組合是Darcy定律接口求解裂隙網(wǎng)絡內的壓力場和達西速度場再加一個傳遞方程接口求解漿液年齡場。兩個接口之間通過粘度和滲透率耦合。有人會問為什么不直接用“裂隙流”接口COMSOL的Fracture Flow模塊確實能模擬裂隙開度對滲透率的影響但它通常按邊界特征使用更適合單條或少數(shù)幾條裂隙的精細模擬當裂隙網(wǎng)絡有上百條圓盤且需要耦合年齡場時用統(tǒng)一的Darcy接口反而更容易控制參數(shù)和網(wǎng)格。4.2 裂隙的“薄層”困境與邊界特征三維裂隙開度只有0.1mm而圓盤直徑動輒幾米如果按實體薄層建模長寬比會到10?甚至10?網(wǎng)格幾乎不可能合格。這里必須走“面裂隙等效滲透系數(shù)”的路線。在Darcy定律接口里裂隙不必建成三維體而是用內部邊界上的“Fracture”特征來定義。每個圓盤圓面作為邊界賦予兩個屬性開度b和內部滲透率k其中k b2/12。Darcy接口自動將邊界上的切向流動用立方定律等效與真實縫隙流一致。這樣幾何對象是二維面網(wǎng)格量大幅下降。具體操作上把我前面生成的所有圓盤面作為一個幾何邊界集合導入在物理場設置里統(tǒng)一選擇這些邊界為“Fracture”輸入開度分布可以是隨機場和滲透率表達式?；鶐r域本身可以保留很低的滲透率比如1e-17 m2作為電性背景也可以直接把基巖滲透率設為0只允許裂隙邊界傳導。后一種更干凈但求解器可能會抱怨零滲透區(qū)域需要在基巖域內給一個極小值防奇異。4.3 注漿孔、出口與壓力邊界設定注漿孔在高維幾何里怎么定義最省事的方法是在某條主裂隙圓盤中心附近設置一個小圓柱邊界圓柱半徑1020cm模擬孔段在這個邊界上設定注漿壓力p_in或注入流量Q。如果研究恒定壓力注漿邊界條件用固定壓力如果研究恒定流量注漿用流量邊界。工程上更常見的是先定壓力上限所以我會推薦“壓力控制流量觀測”。裂隙網(wǎng)絡的邊緣處理是另一個容易出錯的地方。隨機圓盤在模擬域邊界被截斷后截斷邊緣會形成人工敞開邊界漿液會從那里“虛漏”出去。正確做法是在所有圓盤邊界上默認設為零流量邊界只有真正需要模擬串漿出口的位置才打開滲流邊界條件或常壓邊界。如果模型里要考慮地下水初始壓力就要在初始值里給靜水壓場。否則初始壓力全零注漿孔一旦加壓壓力波傳播初期會有一個短暫沖擊過程容易造成時間步振蕩。我建議總是給一個小的初始水力梯度哪怕是0.1MPa的均勻初始壓力也能顯著改善早期收斂。4.4 網(wǎng)格策略與求解順序網(wǎng)格的核心矛盾在圓盤交叉線附近。兩條裂隙相交處的跡線是流量集中區(qū)網(wǎng)格必須在跡線附近加密。COMSOL的自動網(wǎng)格不太懂你這套“面裂隙網(wǎng)絡”的拓撲需要手動添加“邊尺寸”控制在幾何邊集里選中所有圓盤相交線網(wǎng)格尺寸設為最小開度的35倍。求解順序上強烈建議分兩步走。第一步是穩(wěn)態(tài)求解只打開Darcy定律用注漿壓力求初始速度場關閉年齡場。第二步是瞬態(tài)求解打開年齡場開啟耦合時間步用柯朗數(shù)限制Δt Δx/|u|等于保證漿液前鋒每個時間步最多走過一個網(wǎng)格單元否則濃度/年齡場數(shù)值震蕩。這句話值得記在電腦前。5. 調參踩坑實錄收斂崩潰、虛漏與量綱陷阱5.1 粘度突變引起的瞬態(tài)振蕩第一次跑瞬態(tài)時我遇到了典型的“粘度階梯震蕩”漿液前鋒經(jīng)過某個網(wǎng)格的瞬間該處粘度從水的粘度跳到漿液粘度壓力場立刻反轉達西速度開始擺動年齡場前鋒出現(xiàn)鋸齒。根因是粘度梯度太陡而達西滲透系數(shù)與粘度成反比粘度翻5倍等于滲透系數(shù)突然降為五分之一。修復辦法有三個我都試過給粘度變化加平滑過渡用δ 0.5·atan(a·(T-T0))/π0.5這樣的平滑因子讓粘度在幾個時間步內漸變限制柯朗數(shù)在0.7以下把年齡場擴散系數(shù)D_T設為一個極小的正值如1e-9 m2/s來抑制數(shù)值尖峰。三個同時用基本能穩(wěn)定。不要把D_T設得太大否則年齡場會被人為抹平粘度空間差異失去意義。這個參數(shù)本質上只是數(shù)值正則化項物理上的微觀混合作用很小。5.2 圓盤邊緣虛漏問題我最早的一版模型中漿液擴散半徑異常偏大速度場顯示很多漿液從圓盤外緣沖進周圍基巖。檢查發(fā)現(xiàn)是圓盤截斷邊緣沒有設置絕緣條件默認成了開放邊界。這是個非常隱蔽的問題因為COMSOL對內部邊界默認是連續(xù)通量而圓盤邊緣恰好是幾何邊界默認屬于“邊界”沒有指定通量時常被解釋為零壓力或開放出口。解決辦法是在Darcy接口里把所有圓盤邊緣邊集統(tǒng)一設置為“零通量”邊界。做完之后擴散半徑立即縮小了將近三成。這也提醒我驗證瞬態(tài)結果時先檢查未注漿區(qū)域的達西速度是否嚴格為零如果背景區(qū)有流量說明邊界或滲透設置有問題不要急著分析擴散規(guī)律。5.3 開度太小或跨度過大的幾何噩夢當圓盤半徑是5m、開度是0.05mm時立方定律滲透率kb2/12只有2e-13 m2雖然數(shù)學上沒問題但基巖背景滲透率如果沒設好漿液會優(yōu)先在基巖里“漏”而非在裂隙里走。另一個問題是當開度過小壓力梯度會極其陡峭預計壓力場在高滲透裂隙和低滲透基巖交界面附近產生數(shù)值偽影。我的經(jīng)驗是開度不要小于0.1mm建模除非有非常確鑿的現(xiàn)場數(shù)據(jù)小于0.1mm的裂隙直接放到背景滲透率里不進圓盤網(wǎng)絡。理由很實際小于0.1mm的裂隙對注漿的貢獻基本可以忽略漿液中的顆粒也難以擠進去。這樣既保住物理真實性也保住數(shù)值穩(wěn)定性。5.4 計算效率別一上來就跑全瞬態(tài)全耦合全瞬態(tài)全耦合的意思是每一時間步都把達西方程和年齡場方程同步迭代求解。這種做法在裂隙數(shù)少的時候什么都能跑但到了200條裂隙、網(wǎng)格30萬時每一步都是煎熬。我最后采用的策略是“時間步內解耦”在單個時間步內先用上一步粘度場求解達西方程得到速度場再凍結速度場求解年齡場然后用新年齡場更新粘度進入下一步。這種分離迭代方法犧牲了內部一致性但換來的是時間步可以取得更大、單步成本更低而且年齡場對粘度的反饋本來就有時間滯后分離處理在物理上也說得通。實測在同樣精度要求下分離迭代比全耦合快了近四倍。如果你對收斂性有把握甚至可以每兩三個時間步才更新一次粘度場。6. 結果解讀與工程對標6.1 從后處理中提取哪些指標跑完模型第一件事不是截圖而是提取四類指標漿液前鋒半徑、注漿壓力-時間曲線、累計注入量、各主裂隙分配比例。漿液前鋒半徑用年齡場或漿液體積分數(shù)的等值面來定義比如取T漿液初凝時間的等值面。注漿壓力-時間曲線在恒定流量注漿模式下最有價值它直接告訴你泵壓在什么時候開始顯著上升這就是現(xiàn)場停注判據(jù)的預演。累計注入量用注漿孔邊界的總通量對時間積分這個值乘以一個富裕系數(shù)就是現(xiàn)場注漿量預估。各主裂隙分配比例是離散裂隙模型相比等效模型最大的優(yōu)勢——你可以從結果里看到每條裂隙分走了多少漿液高滲透主裂隙往往占80%以上的流量剩余裂隙全是陪跑。這個信息可用于灌漿孔的加密設計。6.2 多次隨機實現(xiàn)與設計包絡單次隨機圓盤生成只是一次抽樣沒有任何統(tǒng)計代表性。同參數(shù)下?lián)Q一組隨機種子生成的裂隙網(wǎng)絡可能完全不同。我做了20次隨機實現(xiàn)統(tǒng)計漿液前鋒在固定注漿時長下的累計概率分布發(fā)現(xiàn)擴散半徑的P10和P90差異可以超過40%。所以工程判斷不能建立在單次模擬上。至少跑1520個隨機實現(xiàn)統(tǒng)計擴散半徑的中位數(shù)和90%分位數(shù)把P90作為設計參考。這也是DFN方法在工程界的核心價值不是預測一條精確路徑而是給出一組帶概率的擴散場景讓決策者知道最壞情況在哪里。我在處理一個咨詢項目時用20個隨機實現(xiàn)得到的結果是某優(yōu)勢產狀組裂隙方向上P90擴散半徑達到18m而正交方向只有6m。這個方向性差異直接導致注漿孔排距從原來的8m縮到5m否則串漿概率太高。6.3 用壓水試驗數(shù)據(jù)校準模型模型參數(shù)不能靠猜必須用現(xiàn)場壓水試驗對標。最直接的橋梁是利用呂榮值換算裂隙等效滲透系數(shù)再換算等效開度。操作思路現(xiàn)場壓水得到透水率呂榮值換算成滲透系數(shù)K然后假設裂隙網(wǎng)絡平均體密度P32和平均半徑已知用立方定律把K映射到等效平均開度b_eq。把模型里開度分布的中位數(shù)調成b_eq再復算一次擴散半徑與現(xiàn)場示蹤或相鄰鉆孔見漿記錄做對比。如果偏差在20%以內模型可信度就相當高了如果偏差過大優(yōu)先懷疑P32和半徑分布而不是粘度參數(shù)。我在一個模擬項目里用這個流程把擴散半徑的觀測誤差從45%壓到了12%左右。關鍵是先標開度后標P32粘度參數(shù)放在最后微調因為粘度只影響擴散過程的時間演化不影響最終可達范圍的主要判定。6.3.1 為什么不建議反過來擬合有人喜歡用“注漿量實測反推模型全部參數(shù)”這通常是不適定的。注漿量對多個參數(shù)同時敏感一個注漿量數(shù)據(jù)點根本不夠約束四五種自由參數(shù)。現(xiàn)場數(shù)據(jù)再豐富也覆蓋不了參數(shù)空間的搜索盲目反演只會得到看起來漂亮、實際完全脫離物理的參數(shù)組合。錨定開度、再調P32、最后微調粘度才是穩(wěn)定路徑。結尾的一點個人體會這套三維離散裂隙注漿模型跑通之后我最大的體會是模型能不能收斂通常不取決于物理場設置多復雜而取決于幾何清理和邊界條件想得有多細。隨機圓盤從統(tǒng)計生成到COMSOL幾何中間那步“清理孤立裂隙和碎片”如果偷懶后面所有網(wǎng)格和求解都是在還債。黏度空間場部分我強烈建議先按“年齡驅動型”做第一版收斂穩(wěn)、參數(shù)直觀、實驗也好標定剪切變稀型等建模和數(shù)值經(jīng)驗積累夠了再上不然很容易掉進“低粘度加速→高流速→更低粘度”的正反饋振蕩里。我后來做雙液注漿時又把這個年齡場模型擴展成了兩個反應組分互相競爭的結構核心框架還是這套——先讓裂隙拓撲對再讓流變行為對最后再說收斂性和效率的事。下一步我準備把應力-滲流耦合加進去讓裂隙開度隨注漿壓力動態(tài)張開這是目前模型最大的簡化點。工程上注漿壓力確實會把裂隙頂開開度一變整個擴散格局就會變。