導水裂隙帶發育高度的微震監測研究

叢 森，程建遠，王云宏，段建華

(1.西安科技大學地質與環境學院，陜西 西安 710054；2.中煤科工集團西安研究院有限公司，陜西 西安 710077)

導水裂隙帶發育高度的微震監測研究

叢 森1，程建遠2，王云宏2，段建華2

(1.西安科技大學地質與環境學院，陜西 西安 710054；2.中煤科工集團西安研究院有限公司，陜西 西安 710077)

基于井-地聯合微震監測工作面導水裂隙帶發育高度的監測成果，分析了微震事件的空間分布特征和在垂直方向的數量特征，并對微震事件的數量和空間分布特征與導水裂隙帶發育高度的關系進行了研究；采用FLAC3D數值模擬技術，對比研究了導水裂隙帶發育高度。研究結果表明：井-地聯合微震監測具有較高的定位精度；微震事件的平面位置主要集中在采空區和回采巷道附近，垂直方向上主要集中在煤層頂板以上110 m范圍內，微震事件數量在頂板以上110 m后急劇減少；數值模擬結果與微震監測結果基本吻合；研究成果為導水裂隙帶發育高度的監測提供了新的方法。

微震監測；微震事件；導水裂隙帶；數值模擬

0 引言

煤礦開采過程中，采掘活動引起采掘工作面周圍巖體應力的重新分布，圍巖發生變形和破壞。微震監測技術是利用煤巖受力變形和破壞時產生的地震波來監測煤巖穩定性的技術方法[1]，通過在監測區域不同方位設置傳感器，對微震事件進行實時記錄，然后計算微震發生的時間和位置，進而推斷煤巖體的破壞情況。微震事件記錄的是煤巖體變形和破壞逐漸發生、發展的過程，利用微震監測就可以記錄和描述導水裂隙帶形成和發展的過程，為導水裂隙帶高度的確定提供依據。

確定導水裂隙帶發育高度的常用方法有理論計算法、數值模擬法、相似材料模擬法和現場實測法。理論計算常用經驗公式法，原國家煤炭工業局2000年頒發的《建筑物、水體、鐵路及主要井巷煤柱留設與壓煤開釆規范》中給出了相關的經驗公式，但該公式具有一定的局限性，已經不能適應特厚煤層的綜放開采[2-4]。由于數值模擬與相似模擬受模擬參數等條件的限制，通常模擬結果只作為一種定性的趨勢分析，而不能作為定量的精確分析，大都結合其他方法使用。近年來發展了多種以鉆孔為基礎的現場實測法，如井下仰孔雙端封堵測漏法、鉆孔漏失量觀測法、鉆孔聲速法和鉆孔電視法等。然而這些方法受鉆孔位置的限制，在空間上為“一孔之見”的瞬時值，無法滿足對整個工作面導水裂隙帶全面、動態和高精度的探查要求[5-8]。

微震監測技術自21世紀初開始在我國礦山水害防治、沖擊地壓防治等方面開展了廣泛的應用[9-13]，具有連續、實時、三維和動態觀測的特點，可實現對整個工作面導水裂隙帶發育的長期動態監測，因此微震監測技術是導水裂隙帶發育高度探測的有力手段。本文采用中煤科工集團西安研究院自主研發的防爆型YTZ3微震監測系統，在某礦首采面開展了導水裂隙帶發育高度的現場監測，并結合數值模擬對導水裂隙帶的發育高度進行了綜合研究。

1 工作面概況及監測系統布置

1.1 工作面概況

某礦位于彬長礦區西北部，主采煤層為4#煤層，首采面自2015年12月開始回采，工作面埋深1 000 m左右，其走向長360 m，寬130 m，煤層為近水平煤層，平均厚度10 m，采用綜放開采技術。主要威脅含水層為白堊系洛河組砂巖含水層，洛河組砂巖下段與煤層頂面距離約60 m，工作面范圍內洛河組厚度約430 m，巖性多為砂巖、砂礫巖，以中-粗粒砂巖為主。

1.2 YTZ3微震監測系統布置

以往微震監測系統在煤礦上應用均以井下部署為主，但由于井下空間有限，觀測系統幾乎在一個平面內，導致震源定位在深度方向存在一定誤差，這給微震監測技術在頂板裂隙帶發育高度監測上的應用帶來了一定的困難。為提高微震監測垂直方向的定位精度，監測臺網應該在已有平面布置的基礎上向立體監測方向發展，即形成立體的微震監測臺網。本次監測采用中煤科工集團西安研究院自主研發的YTZ3微震監測系統，在地面和井下聯合布置監測設備，形成立體的井-地聯合微震監測系統，提高微震監測系統在垂直方向的定位精度。

本次監測共設置測點14個，在工作面輔助運輸平巷和回風平巷頂板各布置3個微震監測點，地面布置8個微震監測點。監測系統布置如圖1所示，圖1中紅色點代表井下監測點，綠色點代表地面監測點。

圖1 監測系統布置圖

為檢驗微震監測系統的定位精度，根據工作面開采地質條件，確定工作面內深孔爆破的標定方案，共實施四次爆破，圖2所示為標定炮波形圖，前八道為地面記錄，后五道為井下記錄。爆破結束后對接收到的微震數據計算處理，震源實際坐標與計算求得坐標之間最大誤差為4.12 m，最小為1.89 m，四炮的平均誤差為3.04 m，能夠滿足5 m以內的定位精度，因此該結果滿足工程的精度要求。

圖2 標定炮波形圖

圖3 微震事件波形圖

2 微震監測結果及分析

2.1 垂直方向微震事件數量特征

在2016年3月16日～2016年4月26日期間，共接收到有效微震事件8 547個，微震事件波形如圖3所示，前八道為地面記錄，后六道為井下記錄。對監測到的微震事件按深度每隔10 m統計微震事件個數，圖4為在各深度區間內對應的微震事件個數柱狀圖，煤層底板標為+100 m。

由圖4可知，在標高0～+300 m范圍內微震事件數目分布大致呈“單峰”狀突起形態。煤層底板的微震事件數量較少，自標高+100 m向上微震事件數目急劇增加，開采煤層所在標高+100～+110 m

區段內事件數目達722個，由煤層開采而引起的微震事件較多；在+160～+170 m區段內事件數目達到最大值之后開始下降，大于頂板標高+210 m后微震事件數目急劇減少，+220 m以上只有零星微震事件發生；在+100～+210 m區段內微震事件數目均大于450個，該區段內微震事件十分集中，說明該區段內巖層變形與破壞較為劇烈；標高+210～+300 m區段微震事件數目隨標高的增加而急劇減少，該區段內巖層變形與破壞逐漸趨于穩定。

2.2 微震事件空間分布特征

圖5、圖6分別為微震事件定位結果3D分布圖和XY平面分布圖，圖5、圖6中每一點代表一個微震事件。由3D分布圖和XY平面圖可知，本次監測到的微震事件主要集中在101工作面采空區及回采巷道上方，回采巷道名稱見圖2(b)。回采巷道附近微震事件較為集中，并且在超前停采線一定距離內也密集分布著微震事件，經分析主要由以下原因造成：①受工作面超前支撐應力的影響，停采線前方頂板巖層發生破裂；②由圖2(b)可知，停采線前方聯絡巷道較多，受工作面回采影響易產生應力集中造成覆巖的變形與破壞；③巷道內機械設備運行及移動、材料搬運和人員活動較為頻繁。除接收到101工作面內部及其附近的微震事件外，圖7中綠色、紅色和黃色橢圓內也接收到較為集中的微震事件，結合煤礦掘進進度與微震事件空間位置綜合分析可知上述區域微震事件分別由103工作面巷道掘進、盤區水倉掘進和201工作面巷道掘進所引起。

圖4 微震事件垂向分布柱狀圖

圖5 微震事件定位結果3D顯示圖

圖6 微震事件定位結果XY平面圖

圖7為沿工作面走向繪制的微震事件XZ方向剖面圖，圖7中背景是以鉆孔地質資料為基礎的地層信息，每一個球形點代表一個微震事件。由圖8可知，微震事件在垂直方向主要分布在標高+100～+210 m的范圍內，該范圍內微震事件十分密集。對應的地層范圍是煤層向上至洛河組下段砂巖，由密集分布的微震事件可知該范圍內的覆巖發生了充分的破壞。+210 m以上僅存在少量的微震事件，空間內分布較為稀疏。

微震事件在平面上主要分布在工作面采空區及其回采巷道附近，垂向上主要分布在煤層底板至標高+210 m的范圍內，該空間內覆巖變形與破壞較為充分，是導水裂隙帶發育較充分區段；標高+210～+300 m范圍內不規則分布著少量微震事件，事件數目隨標高的增加而急劇減少，分布范圍相對分散，覆巖變形與破壞程度較為微弱，因此該范圍內頂板巖層僅在局部區域發生了小規模的破裂，這些破裂尚不能構成貫通性的裂隙，應處于導水裂隙帶與彎曲下沉帶的過渡區段。為安全起見，導水裂隙帶的頂界面應選在微震事件數目隨標高而急劇減小的區段內，綜合以上分析判定某礦首采面導水裂隙帶頂界面的最大發育高度為標高+220 m(煤層底板標高為+100 m，煤層厚度10 m)，即導水裂隙帶的發育高度為110 m。

3 導水裂隙帶發育高度的數值模擬

3.1 數值計算模型及參數

本次數值模擬采用FLAC3D有限差分軟件進行，由于煤層傾角較小，計算模型設為水平模型，模型由Generate命令生成，尺寸為長(X)×寬(Y)×高(Z)=560 m×330 m×300 m，模型沿X軸正向開挖360 m，切眼與停采線距模型左右邊界各100 m；工作面傾斜長130 m，上、下順槽距模型前后邊界各100 m；按照煤層地質綜合柱狀圖選定煤層底板向下取21 m，煤層厚度取10 m，模型總高300 m。三維數值計算模型如圖8所示。整個模型由164 472個單元、174 420個節點組成，網格規格為10 m×10 m。基于現場勘探工程獲取的頂板覆巖組合結構及地層性質的相似性概括為23層模型材料，依據已有的煤巖物理力學性質資料，模型中煤巖物理力學參數如表1所示。采用Mohr-Coulomb本構模型，模型頂部以等效應力的方式設置頂界面應力邊界條件，模型的底部邊界固定，前后和兩側限制水平運動。

圖7 微震事件定位結果走向XZ剖面圖

圖8 三維地質模型圖

3.2 模擬結果分析

模擬工作面回采長度為360 m，采高4 m，工作面每次推進10 m，直到推進至360 m。計算結果處理過程中，將采空區上方相互貫通的拉應力塑性區作為導水裂隙帶，覆巖塑性區發育的最高點作為導水裂隙帶發育的最大高度[14]。回采結束后，沿工作面走向按一定距離截取工作面傾斜剖面，觀察其塑性區發育的最高點。根據模擬結果可知，隨著工作面的推進導水裂隙帶發育高度逐漸增加，當工作面推進至一定位置時發育到最高點，之后保持不變，但在停采線后方一定距離內呈下降趨勢，導水裂隙帶發育高度走勢如圖9所示，數值模擬導水裂隙帶發育的最大高度為108 m。

導水裂隙帶的發育具有一定的時間滯后性和距離滯后性(時間滯后：某一位置導水裂隙帶的發育要經過一段時間后才會趨于穩定；距離滯后：某一位置導水裂隙帶發育到最大高度時工作面已向前推進了一段距離)。因此，在工作面回采結束后，導致工作面后方一定距離內的頂板裂隙沒有發育至最高點，即整個工作面導水裂隙帶的發育呈現出“發生→發展→最高→穩定→回降”的發育規律。

表1 煤巖體力學參數表

圖9 導水裂隙帶發育高度統計

4 結論

根據上述對微震事件的特征分析，結合數值模擬，可有如下結論。

1)井-地聯合微震監測系統形成的立體監測網絡具有較高的定位精度，能夠滿足對導水裂隙帶發育高度的監測要求。

2)基于礦井開采的地質條件，系統分析微震監測結果，綜合研究研究得到某礦首采面導水裂隙帶頂界面發育的最大高度為標高+220 m，即導水裂隙帶的發育高度為110 m。

3)數值模擬得到某礦首采面導水裂隙帶發育的最大高度為108 m，對數值模擬結果與微震監測結果進行對比驗證，二者得出的導水裂隙帶發育最大高度基本一致，微震監測結果具有較高的可靠性。

4)研究成果為導水裂隙帶的探查提供了一種新的、有效的監測方法。

[1] 王桂峰, 竇林名, 李振雷, 等. 沖擊礦壓空間孕育機制及其微震特征分析[J]. 采礦與安全工程學報, 2014, 31(1): 41-48.

[2] 許延春, 李俊成, 劉世奇, 等. 綜放開采覆巖“兩帶”高度的計算公式及適用性分析[J]. 煤礦開采, 2011,16(2): 4-7．

[3] 國家煤炭工業局. 建筑物、水體、鐵路及主要井巷煤柱留設與壓煤開采規程[M]. 北京: 煤炭工業出版社,2000．

[4] 白利民, 尹尚先, 李文. 綜采一次采全高頂板導水裂縫帶發育高度的計算公式及適用性分析[J]. 煤田地質與勘探, 2013, 41(5): 36-39.

[5] 張玉軍, 康永華. 覆巖破壞規律探測技術的發展及評價[J]. 煤礦開采, 2005, 10(2): 10-12．

[6] 任奮華, 蔡美峰, 來興平, 等. 采空區覆巖破壞高度監測分析[J]. 北京科技大學學報, 2004, 26(2):40-45.

[7] 劉英峰. 深埋特厚煤層綜放開采覆巖破壞高度綜合探查技術[J]. 中國煤炭地質, 2013, 25(7): 26-30.

[8] 范志勝, 程敏. 應用井下仰孔分段注水觀測導水斷裂帶高度[J]. 礦業安全與環保, 2010: 37(6): 75-77.

[9] 于克君, 駱循, 張興民. 煤層頂板“兩帶”高度的微地震監測技術[J]. 煤田地質與勘探, 2002, 30(1):47-51.

[10] 張興民, 于克君, 席京德, 等. 微地震技術在煤礦“兩帶”監測領域的研究與應用[J]. 煤炭學報, 2000,25(6): 566-570.

[11] 姜福興, XUN Luo．微震監測技術在礦井巖層破裂監測中的應用[J]．巖土工程學報, 2002, 24(2):147-149.

[12] 陸菜平, 竇林名, 曹安業, 等. 深部高應力集中區域礦震活動規律研究[J]. 巖石力學與工程學報,2008, 27(11): 2302-2308.

[13] 夏永學, 康立軍, 齊慶新, 等. 基于微震監測的5個指標及其在沖擊地壓預測中的應用[J]. 煤炭學報, 2010, 35(12):2012-2016.

[14] 韓軍, 張宏偉, 高照宇, 等. 巨厚煤層軟弱覆巖分層綜放開采覆巖破壞高度研究[J]. 采礦與安全工程學報, 2016, 33(2): 226-237.

Study on microseismic monitoring of height of water flowing fracture zone

CONG Sen1， CHENG Jianyuan2， WANG Yunhong2， DUAN Jianhua2

(1.College of Geology & Environment，Xi’an University of Science and Technology，Xi’an 710054, China; 2.Xi’an Research Institute of China Coal Technology & Engineering Group Corp，Xi’an 710077, China)

Based on microseismic monitoring results of the working face water flowing fractured zone development, the authors analyzed the spatial distribution characteristic of microseismic events and the quantitative characteristics in vertical direction. The relationship between the spatial distribution characteristics and the number of seismic events and height of water flowing fracture zone has been analyzed. A numerical simulation was performed to study height of water flowing fracture zone with FLAC3D by contrast. The research results have shown that microseismic monitoring has been proved to be of high positioning accuracy. The plane position of microseismic events mainly concentrated in the goaf and near the extraction opening, vertical position of microseismic events is mainly concentrated in the 110m above the coal roof, above of which, the number of microseismic events drastically reduced; the numerical simulation results agreed with the results of microseismic monitoring; the analysis result provides a new method for microseismic monitoring of height of water flowing fracture zone.

microseismic monitoring；microseismic event；water flowing fractured zone； numerical simulation

2016-12-21

“十三五”科技重大專項資助項目資助(編號：2016ZX05045-003-05)

叢森(1987-)，男，內蒙古通遼人，博士研究生，主要從事微震監測技術的研究，E-mail：congsenck@163.com。

TD76

A

1004-4051(2017)03-0126-06