動水條件下孔距對多孔注漿擴(kuò)散的影響

摘要：以城門山銅礦礦區(qū)水文地質(zhì)條件和注漿材料及漿液性能實(shí)驗(yàn)研究成果為基礎(chǔ)，模擬不同配比的漿液在溶洞內(nèi)存在動水條件下的注漿效果，深入探究孔距對于多孔注漿搭接的影響，最后根據(jù)實(shí)際的地層條件和動水條件，為合理的漿液配比、孔距設(shè)置以及針對性的注漿方案提供參考依據(jù)。

關(guān)鍵詞：注漿； 動水； 孔距； 數(shù)值模擬

中國分類號：TU475+.9A

［定稿日期］2022-04-11

［作者簡介］王穎卓（1997—），男，碩士，研究為巖土工程。

1 動水溶洞注漿漿液擴(kuò)散數(shù)值模擬分析

我國與巖溶相關(guān)的礦區(qū)、礦坑突水問題時有發(fā)生，且后果都較為嚴(yán)重。我國北方的石炭系、二疊系煤層發(fā)生突水量在600 m3以上的有200多次。究其原因，都離不開地下水的影響，或是受到奧陶系灰?guī)r中的巖溶水的影響，或是奧陶系中的巖溶水沿?cái)嗔褞А⑾萋渲嚷┤搿Ｎ覈戏剑▉啛釒А駶檸r溶區(qū)）的一些礦區(qū)，突水情況也時有發(fā)生，礦井水主要來自煤層頂板長興灰?guī)r和底板茅口組灰?guī)r的巖溶水。有些礦區(qū)的水由于河水補(bǔ)給，涌水量大，經(jīng)常造成災(zāi)難性的事件。礦坑突水中，還時有突泥現(xiàn)象發(fā)生，是填充型溶洞（填充物為泥砂）中的泥砂隨水涌出，其數(shù)量可達(dá)數(shù)千立方米，堵塞巷道，影響生產(chǎn)，造成嚴(yán)重的生命財(cái)產(chǎn)損失，這種情況多在南方發(fā)生。礦區(qū)的巖溶發(fā)育特征為，在平面上接近礦體的地方巖溶發(fā)育程度高，在豎直方向上0～200 m處發(fā)育強(qiáng)烈，礦區(qū)內(nèi)主要分布著小溶洞，在中帶與外帶分界線的西南部分布著較多的大溶洞，充填率高，同時礦區(qū)地質(zhì)構(gòu)造發(fā)育水網(wǎng)密布，且?guī)r漿巖多次入侵礦區(qū)巖體，導(dǎo)致礦區(qū)水文地質(zhì)條件十分復(fù)雜，對于一般的裂隙型和管道型巖溶可以采用滲透注漿以及壓密注漿的方法進(jìn)行處理，對于溶洞一般采用全充填注漿，但是對于動水條件下大溶洞注漿的注漿材料和注漿工藝還需深入探索。巖溶水資源的枯竭問題礦山開采或者地下工程排水改變了巖溶水文系統(tǒng)內(nèi)的平衡狀態(tài)，過量開采和礦坑疏排水造成巖溶水資源的枯竭，區(qū)域地下水大幅度下降，這是我國很多地區(qū)泉井?dāng)嗔骺萁叩闹饕颉?/p>

巖溶地下水在礦山開采過程中的危害不言而喻，礦坑突水、巖溶塌陷是礦山開采過程中最常見、最嚴(yán)重的巖溶地下水災(zāi)害。我國作為礦產(chǎn)資源大國，在長期生產(chǎn)實(shí)踐中積累豐富的經(jīng)驗(yàn)和大量的資料，為解決礦坑突水和巖溶塌陷等問題打下堅(jiān)實(shí)的基礎(chǔ)。相關(guān)專家提出帷幕注漿堵水為主、排疏水相結(jié)合的理論方法，也初見成效。

本文結(jié)合城門山銅礦地層條件和溶洞的特點(diǎn)，針對不同的漿液配比開展動水條件下的擴(kuò)散數(shù)值模擬，根據(jù)室內(nèi)注漿材料試驗(yàn)和模型試驗(yàn)選取合適的參數(shù)，采用多物理場數(shù)值模擬分析軟件Comsol Multiphysics進(jìn)行漿液擴(kuò)散數(shù)值分析、注漿帷幕整體穩(wěn)定性數(shù)值分析以及溶洞內(nèi)成墻后墻體穩(wěn)定性數(shù)值分析。將不同注漿孔距的漿液擴(kuò)散形態(tài)特征通過數(shù)值模擬軟件進(jìn)行可視化，更加方便為實(shí)際注漿施工提供針對性的漿液配比方案。

2 數(shù)值模擬控制方程

2.1 本構(gòu)方程

湖泥-水泥膏漿在溶洞中的擴(kuò)散機(jī)制十分復(fù)雜，且湖泥-水泥膏漿運(yùn)移受漿液自身理化性質(zhì)影響顯著。對于水灰比較低的普通水泥漿液以及許多化學(xué)類漿液而言，其黏度很高，因此不能將其視為簡單的牛頓流體。并在膏漿流動過程中會形成流核，故而其流動規(guī)律與水在溶洞中的流動規(guī)律相比存在很大的差別。

注漿過程中所用到的湖泥-水泥膏漿，可將其視為假塑性流體或賓漢姆流體。對于假塑性流體，隨剪切速率的增大其黏度逐漸下降，在流動過程中表現(xiàn)出剪切稀釋的特性。若按照某一固定的初始黏度值進(jìn)行計(jì)算，那么所得到的擴(kuò)散半徑將遠(yuǎn)大于實(shí)際值，且流體黏度與溶洞粗糙性之間的耦合效應(yīng)也會產(chǎn)生很大差別。根據(jù)李術(shù)才等的研究成果，膏漿是剪切變稀的，可視為Carreau模型，其本構(gòu)方程表達(dá)式為式（1）。

漿液在動水作用下的擴(kuò)散運(yùn)動我們可以簡化為漿液在注漿壓力作用下的垂向運(yùn)動和漿液在動水作用下的直線運(yùn)動的合運(yùn)動。故將三維空間網(wǎng)絡(luò)結(jié)構(gòu)簡化為由水平面和豎直面相交的結(jié)構(gòu)面，則水平面漿液的擴(kuò)散模型參照單一裂隙中漿液的擴(kuò)散模型；本節(jié)針對豎直結(jié)構(gòu)面漿液的擴(kuò)散進(jìn)行分析。

假設(shè)漿液在注漿壓力下注入空洞中，同時受到注漿壓力、重力和水流的沖刷力，在空洞中擴(kuò)散，現(xiàn)從流域中取任意寬度為b、長度為100 m的矩形流段作為研究對象，任意取出任意微元體，其受力如圖1所示。

研究動水條件下漿液擴(kuò)散規(guī)律就是研究漿液在動水條件下漿液顆粒運(yùn)動軌跡規(guī)律。在物質(zhì)傳遞過程中，由擴(kuò)散引起的通量符合菲克第一定律，這一物理量只取決于擴(kuò)散系數(shù)。可以簡單地將擴(kuò)散系數(shù)理解為，某一平面的單位面外濃度梯度通過該平面的摩爾通量的大小。在多孔介質(zhì)中，有效擴(kuò)散系數(shù)與實(shí)際擴(kuò)散系數(shù)有所不同。這是因?yàn)椋行U(kuò)散截面小于自由流體的截面，而多孔材料中一個點(diǎn)到另一個點(diǎn)的距離也小于分子必須在這2點(diǎn)之間移動的距離（因?yàn)榉肿颖仨氃谶@種材料的固體部分之間移動）。這樣一來，實(shí)際濃度梯度就小于表觀濃度梯度。在菲克第一定律中，通過將擴(kuò)散系數(shù)乘以迂曲度（τgt;1）引入了這一效應(yīng)。在此，通量和濃度梯度是指將這一多孔介質(zhì)作為整體進(jìn)行測量的面積和長度，而不僅僅針對可以發(fā)生擴(kuò)散的流體成分。因此，迂曲度是分子必須沿流體通道在2個點(diǎn)之間移動的實(shí)際距離與這兩點(diǎn)的直線距離的比值。迂曲度與孔隙率ε之間存在標(biāo)準(zhǔn)相關(guān)性：Ψ=ε-13。

3 注漿模擬建立

3.1 模型尺寸與參數(shù)取值

（1） 模型尺寸。

模型為100 m×30 m的各向同性不可壓縮且弱透水介質(zhì)模擬溶洞室，模擬一序孔多孔注漿時分別設(shè)置注漿間距為5、10 m，如圖所示為注漿孔距為10 m。

（2） 參數(shù)選取見表1。

3.2 網(wǎng)格劃分

在網(wǎng)格剖分的選擇中用極密來進(jìn)行處理，共產(chǎn)生234 765個域單元、103 572個邊界單元和8 456個邊單元，在對于其他細(xì)化方案中，其計(jì)算的結(jié)果更加圓潤，與實(shí)際形態(tài)相似，且各網(wǎng)格的數(shù)據(jù)處理較為精細(xì)，更加便于后期數(shù)據(jù)分析記錄，以及與實(shí)驗(yàn)進(jìn)行對比。

4 模擬結(jié)果及分析

多孔動水條件下孔距和漿液配比對擴(kuò)散的影響。

4.1 孔距為5 m，水流速度為0.5 m/s

圖2（a）～圖2（d）是坍落度分別為130、170、210、250 mm的湖泥-水泥膏漿在動水流速為0.5 m/s，孔距為5 m條件下的多孔注漿擴(kuò)散情況，相比于靜水條件下多孔注漿，坍落度為130 mm的漿液中間堆積體開始出現(xiàn)拖尾現(xiàn)象；坍落度170、210和250 mm的漿液堆積體堆積緊密，受水流影響較小。

4.2 孔距為5 m，水流速度為1.0 m/s

圖3（a）～圖3（d）是坍落度分別為130、170、210、250 mm的湖泥-水泥膏漿在動水流速為1.0 m/s，孔距為5 m條件下的多孔注漿擴(kuò)散情況，相比于靜水條件下多孔注漿，坍落度為130 mm的漿液中間堆積體拖尾現(xiàn)象較明顯；坍落度170、210、250 mm的漿液堆積體堆積緊密，但受水流影響，有背水面后移的跡象，整體穩(wěn)定。

4.3 孔距為10 m，流速為0.5 m/s

圖4（a）～圖4（d）是坍落度分別為130、170、210、250 mm的湖泥-水泥膏漿在動水流速為0.5 m/s，和孔距為10 m條件下的多孔注漿擴(kuò)散情況。坍落度為130 mm，相鄰堆積體無搭接，堆積體擴(kuò)散情況和單孔的一致；坍落度為170 mm、210 mm的漿液堆積體受水流影響，搭接處和整體背水面有向水流方向移動的趨勢；坍落度為250 mm時，相鄰堆積體接觸緊密，受水流影響較小。

4.4 孔距為10 m，流速1.0 m/s

圖5（a）～圖5（d）是坍落度分別為130、170、210、250 mm的湖泥-水泥膏漿在動水流速為0.5 m/s，和孔距為10 m條件下的多孔注漿擴(kuò)散情況。坍落度為130 mm，相鄰堆積體無搭接，堆積體擴(kuò)散情況和單孔的一致；坍落度為170 mm、210 mm的漿液堆積體受水流影響，搭接處和整體背水面有向水流方向移動的趨勢加大；坍落度為250 mm時，相鄰堆積體接觸緊密，受水流影響整體背水面有向水流方向移動的趨勢。

5 結(jié)論

對于模擬多孔注漿在動水條件下，注漿孔距小于漿液擴(kuò)散距離時，漿液之間接觸緊密，在水流沖刷下?lián)p失的漿液相對于單孔漿液損失更少；注漿孔距與漿液擴(kuò)散距離時相近時，在水流沖刷下，堆積體能較好保留下來。注漿孔距大于漿液擴(kuò)散距離時，相鄰堆積體之間沒有搭接，和單孔注漿時一樣被水流沖刷。實(shí)際注漿中，動水條件下，注漿成墻第一序孔應(yīng)結(jié)合動水流速下漿液的留存率選擇比漿液擴(kuò)散范圍稍大的注漿孔距。

參考文獻(xiàn)

［1］ 張丁陽. 裂隙巖體動水注漿擴(kuò)散多場耦合機(jī)理研究［D］.中國工業(yè)大學(xué)， 2018.

［2］ 張改玲， 湛鎧瑜， 隋旺華. 水流速度對單裂隙化學(xué)注漿漿液擴(kuò)散影響的試驗(yàn)研究［J］. 煤炭學(xué)報(bào)， 2011， 36（3）：403-406.

［3］ 張更生. 海水入侵機(jī)理及防治措施的三維數(shù)值模擬［D］. 大連： 大連海事大學(xué)， 2007.

［4］ 張貴金， 胡榮宗， 鐘平，等. 新型可控性粘土水泥膏漿試驗(yàn)研究［J］. 水利水電技術(shù)， 2013， 44（2）： 60-70.

［5］ 張貴金， 劉杰， 胡大可，等. 黏土水泥膏漿流變性能及其對灌漿的影響［J］. 長江科學(xué)院院報(bào)， 2017（3）：119-125.

［6］ 張海濤， 許光泉， 朱代雙.目標(biāo)規(guī)劃在礦山地下水疏放管理中應(yīng)用——以淮南潘北煤礦西翼采區(qū)為例［J］.地質(zhì)論評，2018，64（3）：753-763.