突水過程中顆粒運移對破碎砂巖水力特性影響試驗研究

2021-06-18 07:04:56陳學文呂閏生

中國礦業 2021年6期

陳學文，呂閏生

(1.河南建筑職業技術學院，河南鄭州 450064；2.河南理工大學資源環境學院，河南焦作 454003)

近年來，我國中東部礦井淺部煤炭資源已開采殆盡，因此紛紛轉向深部煤炭資源開采[1-3]，但受“三高一擾動”的影響，煤礦突水事故日趨頻繁，嚴重影響了煤礦的安全生產和經濟效益[4-7]。針對礦山突水災變問題，許多專家學者采用理論分析、數值模擬和室內試驗等方法開展了研究并取得了豐碩的成果。在理論方面，錢鳴高等[8]從不同斷層破碎帶類型角度分別研究了相應的突水機制；羅雄文等[9]分析研究了斷層破碎帶內部不良地質含水體特性與突水災害的關系，將觸發突水事故的地質構造類型及災變模式分為五種；WANG等[10]基于有效多孔介質理論研究分析了采動條件下裂隙巖體滲透率變化規律及其對地下水突水影響。數值模擬方面，李廷春等[11]基于連續介質力學和變質量動力學理論，推導了飽和斷層破碎巖體變質量滲流-變形耦合模型并將其編譯入FLAC3D數值軟件中，模擬揭示了富水斷層破碎帶圍巖突水機理；陳國慶等[12]利用COMSL軟件將Darcy流、Brinkman和N-S滲流方程耦合，實現破碎帶隧道突水全過程模擬。室內試驗方面，MA等[13]用采用自發設計加工的試驗系統研究了采動過程中小顆粒運移對誘發斷層突水機理；郭信山等[14]采用斷層影響因子等值線和斷層分維等值線相互融合技術，得到新汶煤田協莊煤礦奧灰巖溶裂隙分布規律。

以上研究大多依托不同的工程背景下的破碎帶突水事故而得出適用于特定條件結論和防治方案。事實上，滲流情況下顆粒運移也是導致破碎帶水害事故發生的重要原因，本文利用自發設計和研制的試驗裝置，定量研究了不同初始孔隙率以及水壓對破碎巖體水力特性的影響，進而揭示致使煤礦突水事故發生的重要內因，由此加深對破碎帶煤礦水害問題的認識和理解，完善和豐富突水機理。

1 室內試驗

1.1 試驗系統及方案

破碎巖體顆粒運移試驗系統由4個部分組成：應力控制系統、顆粒收集系統、顆粒滲流系統和數據采集系統(圖1)。通過控制水壓大小即可進行不用滲透壓下顆粒運移試驗，通過控制應力大小從而控制試樣高度進而可進行不同孔隙率下顆粒運移試驗。一般而言，破碎巖體可視為巖體框架和細顆粒成分兩部分組成，本文選用粒徑為4～8 mm的砂巖顆粒作為試樣框架，選用1～2 mm和2～4 mm粒徑的砂巖顆粒作為細顆粒成分(充填顆粒)，兩部分按照質量比1∶1∶1混合加工成1 500 g的圓柱形試樣。具體試驗方案見表1。

表1 試樣級配及試驗設計Table 1 Test scheme and particle sizes

圖1 破碎巖體顆粒運移試驗系統Fig.1 Particle migration test system for broken rock mass

1.2 試驗過程

試驗首先采用電子稱稱取相應質量的破碎巖石顆粒；然后裝入模具并根據試驗方案的設計采用鐵錘將其搗實，壓實且靜置24 h；隨后將模具拆除，把試樣整體裝入試驗系統中，啟動應力控制器，液壓油推動活塞運動使得試樣高度完全吻合試驗設計時的試樣高度；最后啟動滲流系統，通過水箱及其穩壓器施加穩定的滲透壓(水壓)，電腦實時記錄顆粒收集系統中遷移排出的顆粒質量。

1.3 水力參數計算

1) 滲透率k。根據一維達西定律可得式(1)。

-Δp/Δh=μvk-1

(1)

式中：Δp/Δh為壓力梯度，Pa/m；μ為水在常溫下的黏度，本文取值為8.94×10-4Pa·s；v為平均水流速度，m/s；根據v=Q/(πr2)，Q為水流量，m3，r為試樣半徑，m。

由于本試驗中試樣的頂部與大氣相連通且試樣的高度保持不變，因此式(1)可改寫為式(2)。

ki=μvih0/p

(2)

式中：i為第i時刻相應的參數值；h0為試驗過程中試樣高度，m；p為滲透壓，Pa。

2) 孔隙率n0。為進行不同孔隙率下的顆粒運移試驗，需要施加一定的應力荷載，進而改變試樣度從而達到預期的孔隙率，破碎砂巖試樣試驗前的孔隙率n0計算公式見式(3)。

n0=1-ms/(πr2ρsh0)

(3)

式中：ms為試驗前試樣質量，本文取值為1.5 kg；ρs和h0分別為試驗前試樣的密度和高度，ρs取值2.32×103kg/m3。

試驗過程中，每隔10 s記錄一次排出的顆粒質量，每個試樣的試驗時間為200 s，因此記錄的最大次數為j=20。試驗過程中收集的顆粒總質量計算公式見式(4)。

mj=Δm1+Δm2+…+Δmj(j=1，2，3，…，20)

(4)

試驗過程中顆粒的運移會導致試樣孔隙率的增加，因此孔隙率在試驗過程中的計算公式見式(5)。

nj=1-(ms-mj)/(πr2ρsh0)

(5)

2 試驗結果分析與討論

2.1 滲透率演化分析

破碎巖體滲透率隨時間的變化規律如圖2所示。在整個試驗過程中滲透率的演化特征可分為4個階段，分別是振蕩增長階段、快速增長階段、緩慢增長階段、穩定階段。試驗初期即振蕩增長階段，由于試樣中大顆粒物質的重新排列，引起試樣結構的變化，大顆粒物質在重新排列過程中可能造成其顆粒菱角以及邊緣的破壞進而萌生新的、隨機的導水通道，導致試樣滲透率呈振蕩型增長；隨著試驗的進行，萌生出來的導水通道會搭接、連通形成完整的突水通道，該階段水流速度明顯增大(突水)，使得大量細顆粒物質遷移排出甚至沖出試樣，造成快速增長階段中試樣滲透率急劇增加；在緩慢增長階段中，導水通道明顯擴寬，顆粒在其中遷移的阻力也明顯減少，但是由于大量細顆粒物質已經被沖出試樣，使得殘余的少量細顆粒物質沖出試樣，導致試樣的滲透率進一步增大，但是由于沖出的細顆粒物質的量不多，具有較為緩慢、均勻的速度，因此該階段的滲透率增長速率近似成緩慢的線性增長；最后由于導水通道明顯趨于穩定且幾乎沒有顆粒物質被沖出試樣，滲透率保持穩定。

根據本文所設計的試驗裝置，不宜將破碎砂巖試樣視為雙重結構(骨架-充填物)去研究試樣顆粒和裂縫網絡單因素對滲透率的影響。雖然滲透率的變化可能會影響砂巖樣本內部水力特性的整體變化，但如果想要研究單個因素的變化對滲透率的影響就需要更復雜的測試設備，因為在試驗過程中，一些大顆粒也會分裂成較小的顆粒，所以根據本試驗裝置的功能以及本文的研究目的，認為從整體上研究破碎砂巖試樣水力特性的整體變化較為合理。

2.2 影響因素討論

1) 初始孔隙率。試樣1、試樣2和試樣3分別代表在相同水壓不同初始孔隙率條件下，試樣滲透率變化的結果。在滲透壓相同的情況下，試樣的初始孔隙率越高，試驗結束后其滲透率就越高(圖2)。在振蕩增長階段中，大顆粒物質(試樣框架)重新排列所需時間會隨著初始孔隙率的增加而減少，這種情況與緩慢增長階段完全類似。此外，孔隙率越高，試驗結束后試樣的最終孔隙率也會增大(圖3)。

圖2 試樣滲透率隨時間的變化曲線Fig.2 Change laws of permeability with time

圖3 試樣孔隙率隨時間的變化曲線Fig.3 Change laws of porosity with time

2) 滲透壓。試樣2、試樣4和試樣5分別代表在相同初始孔隙率不同水壓條件下，試樣滲透率的變化結果。滲透壓越高，試驗結束后滲透率和孔隙率均會相應增大。在振蕩增長階段中，由于大顆粒物質的重新排列，滲透壓越高，孔隙率的變化速率會越快，這是因為較大水壓具有較強的攜帶細顆粒物質的能力。

3 數值模擬分析

3.1 數值模型和方案

為了進一步研究突水過程中顆粒運移規律及其對破碎砂巖水力特性的影響，利用離散元數值軟件PFC3D開展不同水壓下破碎砂巖流-固耦合數值分析。根據室內試驗材料和方案設計特點，數值分析中存在的基本特點和假定如下所述。

1) 當模型運行時，允許顆粒相互重疊，重疊量的大小和接觸力的大小有關，且重疊量相對顆粒的尺寸是極其微小的。

2) 薄墻與薄墻之間的間隙為裂縫，薄墻是相互平行、光滑且平直的理想墻面。

3) 流體的滲流速度緩慢，符合達西定律，流體采用網格代替。

4) 滲流過程中僅考慮顆粒在不同水壓作用下的運移效應，忽略滲流過程中顆粒和骨架、顆粒和顆粒之間的碰撞。圖4是建立的離散元數值模型圖，數值試驗設計方案見表2，模型參數見表3。

圖4 離散元數值模型圖Fig.4 DEM model

表2 數值試驗設計Table 2 Designer of numerical simulation

表3 數值試驗材料參數Table 3 Parameters of model

3.2 數值結果分析

根據模擬結果可知，在剛開始施加水頭時有極少部分顆粒立即進入裂隙，絕大部分顆粒往模型中間移動，同時試樣頂部拐角位置的顆粒運動類似于渦流。此外，裂縫中的顆粒垂直運移，裂縫外的顆粒往裂縫處集中(圖5)。隨著模型運行時間的增加，顆粒持續進入裂隙并在水壓作用下遷移出模型，同時顆粒的流動方向會出現一定的變化，渦流消失，拐角處顆粒流動方向指向裂隙，其余顆粒沿垂直方向流動，運移直至從模型中排出(圖6)。

圖5 初期渦流現象和顆粒運動方向分布Fig.5 Initial vortex phenomenon and distribution ofparticle movement direction

圖6 顆粒流動方向的變化Fig.6 Changes of particle movement direction

為進一步研究不同水壓下滲流過程中顆粒的運移總體積與時間的關系，利用FISH語言程序記錄不同運行時間下的裂隙顆粒滲流量，利用WRITE命令將記錄的各物理量以ASCll格式寫入dat文件，再將該dat文件讀入后處理軟件ORIGIN中，所得計算結果如圖7所示。由圖7可知，不同水壓下顆粒體積遷移隨時間的變化曲線形狀基本一致，并且可細分為四個階段：緩慢增加階段、快速增加階段、緩慢增長階段和穩定階段，這四個階段剛好與室內試驗中滲透率與時間的變化關系吻合。因此，數值試驗說明了顆粒的遷移和排出是引起破碎帶滲透率變化的本質因素。除此之外，隨著滲透壓的增大，達到穩定階段的耗時越少，如當壓力梯度由0.3 MPa增加至1.0 MPa時，達到穩定階段的耗時從0.28 s縮減至0.13 s，這就說明在遷移體積相同的情況下，水壓越大，顆粒的體積遷移速率越快。