基于相似模擬的淺埋深煤層覆巖采動裂隙分析

李曉璐，歐陽振華，李少剛

(1.煤炭科學技術研究院有限公司安全分院，北京 100013；2.煤炭資源高效開采與潔凈利用國家重點實驗室(煤炭科學研究總院)，北京 100013 )

采選技術

基于相似模擬的淺埋深煤層覆巖采動裂隙分析

李曉璐1，2，歐陽振華1，2，李少剛1，2

(1.煤炭科學技術研究院有限公司安全分院，北京 100013；2.煤炭資源高效開采與潔凈利用國家重點實驗室(煤炭科學研究總院)，北京 100013 )

摘要：根據(jù)相似模型試驗，探討了運用視電阻率法的物探技術來反演淺埋深覆巖采動裂隙發(fā)育演化規(guī)律的可行性，分析了采動過程中淺埋深煤層覆巖變形破壞前后電性參數(shù)的變化及其與覆巖裂隙發(fā)育范圍、程度之間的關系。研究發(fā)現(xiàn)：視電阻率法可以探測淺埋深煤層覆巖采動裂隙的演化過程，電阻率變化率的等值線圖可以體現(xiàn)裂隙發(fā)育形態(tài)，提出以電阻率變化率作為劃分覆巖破壞“三帶”的標準。在淺埋深煤層的采動過程中，上覆巖層導水裂縫帶的范圍持續(xù)向上方、前方擴展，但裂隙在水平與垂直方向的發(fā)育具有不平衡性，薄基巖層明顯阻礙了裂隙向上的發(fā)育，水平方向的裂隙發(fā)育范圍明顯大于垂直方向。而一旦裂隙發(fā)育突破了基巖層，就會瞬間貫通該基巖層上方的附著巖層或地表。

關鍵詞：淺埋深；視電阻率法；采動裂隙；電阻率變化率；相似模型試驗

在我國西部地區(qū)，相當一部分礦區(qū)具有煤層埋藏淺、基巖薄、地表有松散覆蓋層的特殊地質和采礦條件，在開采的過程中，頂板基巖會在工作面上方形成裂縫，裂縫的貫通會造成地表砂土的涌入，導致井下潰砂等，威脅礦井生產安全，同時這些裂縫易貫通地表，使得地表形成塌陷坑，加劇對地表植被的毀壞，并嚴重影響地面建筑物的安全。

對于覆巖裂隙的探測，我國學者已經開展了大量的研究工作，其中許宏發(fā)等[1]為研究深部巷道圍巖破裂情況，在深部巷道布置鉆孔，對鉆孔內巖體電阻率進行測試，測量圍巖電阻率沿鉆孔深度的變化；康永華等[2]通過彩色鉆孔電視和鉆孔沖洗液消耗量觀測，取得了高水壓含水層作用于原生縱向裂隙發(fā)育巖體條件下的垮落帶、導水斷裂帶高度及發(fā)育特征的現(xiàn)場實測結果；張玉軍等[3]采用鉆孔沖洗液漏失量觀測和鉆孔彩色電視系統(tǒng)，探測了高強度綜放開采覆巖破壞高度；孫亞軍等[4]通過并行網絡電法CT觀測小浪底水庫下采煤導水裂隙發(fā)育情況。劉貴等[5]針對河下開采，通過在工作面不同位置布孔進行鉆孔沖洗液漏失量和彩色電視觀測，并通過物理模擬研究了在各工作面間留設一定寬度隔離煤柱的開采方式的覆巖破壞過程。而對于淺埋深薄基巖的覆巖的采動裂隙的相關研究還是很缺乏的。

煤層覆巖的裂隙發(fā)育可以通過鉆孔沖洗液消耗量法、鉆孔電視或鉆孔成像技術[6]、物探等方法進行探測。任奮華等[7]據(jù)現(xiàn)場鉆孔沖洗液漏失量、鉆孔水位變化以及井下彩色鉆孔電視影像，探明了煤層開采后上覆巖層內垮落帶和導水裂縫帶的發(fā)育高度及其分布形態(tài)；張玉軍等[8]采用鉆孔彩色電視系統(tǒng)對鉆孔內部原生裂隙及采動裂隙的分布及其發(fā)育特點進行了研究；文學寬[9]應用電磁波層析技術(即CT)探測了覆巖的破壞高度；張平松、劉盛東等[10-12]采用電阻率法、立體直流電法、孔巷電法和并行電法觀測，獲得覆巖變形與破壞的發(fā)育規(guī)律。

鉆孔沖洗液消耗量法可以探測覆巖破壞的高度，但當觀測鉆孔內觀測段的巖層中有原生裂隙時，沖洗液消耗量法觀測就不能順利進行，甚至由于對原生裂隙的堵漏不成功而造成觀測失敗。鉆孔電視或鉆孔成像技術雖具有形象、直觀的優(yōu)點，但當巖層完整性較好時，常導致進行觀測的鉆孔內仍然有水，影響觀測效果。因此，筆者根據(jù)相似模型試驗，利用高密度電法，分析淺埋薄基巖的覆巖在采動影響下的裂隙演化規(guī)律，為指導西部淺埋深薄基巖礦區(qū)高效、經濟、安全生產和保護地表生態(tài)環(huán)境提供依據(jù)。

1相似模擬試驗方案

本文以內蒙古雙欣礦業(yè)有限公司楊家村煤礦222203工作面為原型，該工作面是22采區(qū)第二個工作面，埋深約140m，傾角小于5°，煤層結構較簡單，平均煤厚7.3m。工作面四周及上下煤層均未開采，為實體煤區(qū)。工作面采用綜合機械化走向長壁后退式一次采全高采煤法回采，全部垮落法管理頂板。

1.1相似模型試驗設計

相似模型尺寸3200mm×250mm×1500mm，為平面應力模型，如圖1所示。依據(jù)Hoke-Brown經驗強度準則及小試件室內試驗資料，對本次相似模擬實驗所涉及到的煤巖層進行估算，具體參照根據(jù)現(xiàn)場地質資料、巖土體物理力學性質試驗以及上述轉化關系公式，并考慮本次相似模擬材料選取的參考指標主要是抗壓強度和模擬，故得到研究區(qū)域內相似材料模型的物理力學性質參數(shù)如表1所示。

表1　各層物理力學性質

圖1　相似模型圖

1.2試驗內容

模擬煤層回采過程，其中煤層切眼寬度取10m，工作面每次推進10m，共模擬推進237m。通過開挖前在模型表面布置位移變形監(jiān)測點，在煤層上覆巖層中布置視電阻率測試點，來監(jiān)測煤層上方頂板及地表的位移變化、覆巖破斷過程中的視電阻率的變化。

布置6條水平測線進行位移觀測，如圖2所示，其中回采煤層上方5條，煤層下方1條，相鄰測線間距150mm，每條測線13個測點，相鄰測點間距100mm，合計共78個位移測點。

在模型頂表面，沿中心線布設11個視電阻率測試點(D1～D11)，如圖3所示。測量D1點的視電阻率時，以D1點為中心，將電阻率測試儀的兩個接收探頭M、N對稱布置在D1點兩側并不變(MN=50mm)，兩個發(fā)射探頭A、B對稱于測試中心D1點向兩旁按一定間距增加(A、B兩點的起始距離為100mm，向兩邊移動的步距各為50mm，直到移動到模型邊界為止)，并記錄下A、B每次移動后，M、N兩點處測得的電壓和電流，可計算出該點不同深度處的視電阻率[13-15]。D2～D11采用相同方式測試相應點不同深度電阻率。

圖2　位移測點布置示意圖(單位：mm)

圖3　視電阻率測點布置示意圖

2視電阻率測試方法

視電阻率測試方法，是一種陣列勘探方法，源于電剖面法和電測深法，主要是通過高密度布點，進行二維地電斷面測量[16]。由于其具有探測數(shù)據(jù)量大、速度快和觀測精度高的優(yōu)點，現(xiàn)已在巷道工作面、頂板、底板及工作面富水構造探測等方面發(fā)

揮了重要作用，并嘗試用于覆巖“三帶”的探測。

視電阻率測試方法的物理前提是地下介質間的導電性差異。一般來說，在開采前，不同煤巖類型的電阻率也不同。開采過程中，當巖體受拉伸或破碎時，顆粒間接觸松散，產生間隙，孔隙度增加，孔隙充填電阻率為無窮大的空氣，因而巖體的視電阻率升高。覆巖的這一電性特征是覆巖在煤層開采過程中力學性質變化的結果。因而，通過比較煤巖采動前后電阻率的變化，就可以探測到覆巖的裂隙發(fā)育情況。

2.1視電阻率的模塊試驗

在配制模型材料時，適當加入Nacl溶液可模擬不同電性特征的實際地層，而NaCl的濃度和模型的干燥時間都會影響模型的電阻率，為了確保物理模型能夠反映實際的地層情況，首先進行了模塊實驗，以確定合適的材料配比方案。

模塊大小為10cm×10cm，針對實際地層情況，共設置表土層到細粒砂巖共8組模塊，如圖4所示，每組模塊具有不同濃度的NaCl溶液，各層組分詳細配比見表2。

為了降低測量誤差，每組制作三塊相同模塊，測量求平均值。在不同時間段測量電阻率值，記錄不同時間段的電阻率變化，測量結果見表3。

由表3可見，NaCl溶液很好的控制了模塊的電阻率，經過139小時的自然風干后，模塊電阻率與預期電阻率相符，據(jù)此，表2配比能夠滿足實驗要求。

圖4　模塊實物圖

巖性沙子重量/kg碳酸鈣重量/kg石膏重量/kg水重量/kgNaCl溶液/%NaCl重量/kg期望電阻率/Ω.m表土層0.5460.0530.0110.0632.50.01710礫巖0.5360.0420.0320.0630.10.0007120粉砂巖0.5250.0420.0420.06320.01315細粒砂巖0.5250.0530.0210.0630.50.00340粉砂巖0.5250.0420.0320.06320.01315砂質泥巖0.53550.05250.01050.06330.01984510煤0.5460.05250.01050.0630.150.00100890細粒砂巖0.5250.05250.0210.0630.50.00330840

表3　模塊電阻率測量結果

2.2相似模型視電阻率測試準備

為了模擬覆巖破壞過程中電阻率的變化，相似模型不僅需要滿足與研究對象的幾何、容重、強度和彈模等力學性質的相似條件，還要符合電性參數(shù)的物理相似準則，即滿足式(1)。

(1)

式中：ρ1為某一固定層電阻率；ρs為實測視電阻率；M、N分別代表實際地層和模型，才能保證實驗室條件下的模擬結果與現(xiàn)場相同。

在配制模型材料時，根據(jù)模塊實驗驗證的配比方案，添加NaCl溶液。同時，在模型上表面中線位置設置30個測量電極，極距0.1m，采用高密度觀測系統(tǒng)進行測量，根據(jù)模型覆巖破壞程度，共測量22組高密度視電阻率數(shù)據(jù)，從中選取典型時刻的視電阻率數(shù)據(jù)進行反演分析。

3相似模型試驗結果分析

圖5　回采過程中電阻率變化率等值線圖

圖6　回采過程中位移圖

根據(jù)表2可以看出，該工作面上方存在兩個薄基巖層，一個是距工作面頂板10m、厚46m的細粒砂巖層，另一個是距工作面頂板64m、厚58m的礫巖層。由于基巖層限制了裂隙的發(fā)育，在兩個薄基巖層的下方分別形成了拱形的導水裂縫帶和馬鞍形的導水裂縫帶。

當工作面回采100m時(圖5(a))，有小裂隙發(fā)育至地表，裂隙第一次導通地表，地表發(fā)生彎曲下沉。導水裂縫帶上部發(fā)育至頂板上方的90m左右，進入礫巖層；而導水裂縫帶下部初步形成兩邊高、中間底的“馬鞍形”的裂隙形態(tài)，發(fā)育范圍僅限于頂板上方的50m內，位于細粒砂巖層內。

工作面回采115m時(圖5(b))，導水裂縫帶上部發(fā)育至頂板上方的100m左右，向上擴展了10m，但仍處于礫巖層內；導水裂縫帶下部的“馬鞍形”也小幅度的向上發(fā)展，且工作面頂板上方的電阻率變化率也在增大，新的裂隙在產生。

工作面回采140m時(圖5(c))，由于礫巖層限制了導水裂縫帶上部的向上發(fā)育，導水裂縫帶的高度變化不大，依然處于礫巖層內，導水裂縫帶繼續(xù)在水平方向擴展；導水裂縫帶下部的“馬鞍形”的頂部由距頂板的40m向上發(fā)育至55m左右(圖6)。

工作面回采150m和160m時(圖5(d)與圖5(e))，導水裂縫帶的高度沒有變化，還位于礫巖層內；導水裂縫帶下部的“馬鞍形”的頂部向上發(fā)展了20m，貫穿細粒砂巖層，薄基巖的貫通導致其上的粉砂巖層也隨之切落，導水裂縫帶直接擴展過粉砂巖層，“馬鞍形”的頂部距頂板達的90m左右，已接近礫巖層。

當工作面回采170m時(圖5(f))，礫巖層的薄基巖層發(fā)生破斷，導水裂縫帶瞬間發(fā)育至地表，上部松散層也隨基巖層下沉，地表形成臺階下沉(圖6)。同時由于淺埋深的特點，裂隙多次貫通地表，形成多個彼此獨立的裂隙發(fā)育區(qū)。

可以看出，在工作面的回采過程中，隨著覆巖裂隙的發(fā)育，視電阻率發(fā)生規(guī)律性的變化，通過視電阻率的變化來劃分覆巖破壞高度是可行的。在彎曲下沉帶內，巖體的電阻率變化不大，在導水裂縫帶中，上部裂隙發(fā)育弱，電阻率值一般是正常值的1.5倍，下部裂隙發(fā)育，電阻率值是正常值的2.5倍左右；在垮落帶中，電阻率遠比正常值大得多，會達到 4～5倍以上。針對本相似模型，可提出以電阻率變化率作為劃分“三帶”的標準，見式(2)。

(2)

因此，只要選擇合適的電阻率變化率值，通過電阻率變化率的分布情況，可以得到回采過程中覆巖裂隙的發(fā)育規(guī)律。

4結論

1)通過相似模型試驗可以看出，運用視電阻率法的物探技術手段，可以有效克服原生裂隙和覆巖完整性的影響，監(jiān)測到采動裂隙的發(fā)育演化過程。

2)提出以電阻率變化率作為劃分覆巖變形破壞“三帶”的標準，可以有效避免模型風干控制的不確定性和反演中的等值現(xiàn)象(多解性)等對反演電阻率結果的影響。

3)在淺埋深煤層的采動過程中，上覆巖層導水裂縫帶的范圍持續(xù)向上方、前方擴展，但裂隙在水平與垂直方向的發(fā)育具有不平衡性，薄基巖層明顯阻礙了裂隙向上的發(fā)育，水平方向的裂隙發(fā)育范圍明顯大于垂直方向。而一旦裂隙發(fā)育突破了基巖層，就會瞬間貫通該基巖層上方的附著巖層或地表。

參考文獻

[1]許宏發(fā)，錢七虎，王發(fā)軍，等.電阻率法在深部巷道分區(qū)破裂探測中的應用[J].巖石力學與工程學報，2009，28(1):111-119.

[2]康永華，趙開全，劉治國，等.高水壓裂隙巖體綜采覆巖破壞規(guī)律[J].煤炭學報，2009，34(6):721-725.

[3]張玉軍，李鳳明.高強度綜放開采采動覆巖破壞高度及裂隙發(fā)育演化監(jiān)測分析[J].巖石力學與工程學報，2011，30(S1):2994-3001.

[4]孫亞軍，徐智敏，董青紅.小浪底水庫下采煤導水裂隙發(fā)育監(jiān)測與模擬研究[J].巖石力學與工程學報，2009，28(2):238-245.

[5]劉貴，張華興，劉治國，等.河下綜放開采覆巖破壞發(fā)育特征實測及模擬研究[J].煤炭學報，2013，38(6):987-993.

[6]王傳嬰，LAWKTim.鉆孔攝像技術的發(fā)展與現(xiàn)狀[J].巖石力學與工程學報，2005，24(19):3440-3448.

[7]任奮華，蔡美峰，來興平，等.采空區(qū)覆巖破壞高度監(jiān)測分析[J].北京科技大學學報，2004，26(2):115-117.

[8]張玉軍，張華興，陳佩佩.覆巖及采動巖體裂隙場分布特征的可視化探測[J].煤炭學報，2008，33(11):1216-1218.

[9]文學寬.CT探測覆巖破壞高度的試驗研究[J].煤炭學報，1998，23(3):300-304.

[10]劉盛東，吳榮新，張平松.高密度電阻率法觀測煤層上覆巖層破壞[J].煤炭科學技術，2001，29(4)：18-20.

[11]張平松，劉盛東，吳榮新.采煤面覆巖變形與破壞立體電法動態(tài)測試[J].巖石力學與工程學報，2009，28(9):1870-1875.

[12]吳榮新，張衛(wèi)，張平松.并行電法監(jiān)測工作面“垮落帶”巖層動態(tài)變化[J].煤炭學報，2012，37(4):571-577.

[13]于師建，程久龍，王玉和． 覆巖破壞視電阻率變化特征研究[J]．煤炭學報，1999，24( 5) ：457-460．

[14]程久龍，于師建.覆巖變形破壞電阻率響應特征的模擬實驗研究[J].地球物理學報，2000，43(5)：699-706.

[15]于師建，程久龍.覆巖導水裂縫帶電阻率響應特征[J].巖土工程學報，2000，22(3):336-339.]

[16]劉盛東，劉靜，岳建華.中國礦井物探技術發(fā)展現(xiàn)狀和關鍵問題[J].煤炭學報，2014，39(1):19-25.

Analysis of mining-induced fractures evolution in overburden about shallow seam based on simulation model test

LIXiao-lu1，2，OUYANGZhen-hua1，2，LIShao-gang1，2

(1.MineSafetyTechnologyBranch，ChinaCoalResearchInstitute，Beijing100013，China；2.StateKeyLaboratoryofCoalResourceHigh-efficiency&CleanUtilization，

ChinaCoalResearchInstitute，Beijing100013，China)

Abstract:On the basis of similar model test，the feasibility which mining-induced fractures evolution in overburden about shallow seam could be inversed by apparent resistivity method was discussed.The electric parameter changes in overlying strata during the course of deformation and failure and the relation between it and mining-induced fractures are canalized.The results show that apparent resistivity method can detect mining-induced fractures evolution in overlying strata about shallow seam.The isoline map of change rate in electrical resistivity could give expression to fracture development，so it can be the new division standard of the three failured strata zones.During mining active process in shallow seam，the range of fractured-zone in overlying strata is extended upward and forward，but it is imbalance in horizontal and vertical direction.The horizontal range is obviously larger than vertical，because the propagation of fractures is impeded by the thin bedrock.But once the thin bedrock is broken through by fractures，the upper rock or surface is instantly transfixion.

Key words：shallow seam；apparent resistivity method；mining-induced fractures；change rate in electrical resistivity；similar model test

收稿日期：2015-06-17

基金項目：國家重點基礎研究發(fā)展計劃(973) 項目資助 (編號：2012CB724208)；國家自然科學基金項目資助(編號：51174272；51404140)

作者簡介：李曉璐(1977-)，女，河北邯鄲人，副研究員。E-mail: lxl_bj@126.com。

中圖分類號：TD821

文獻標識碼：A

文章編號：1004-4051(2016)04-0072-05