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