叢 森，程建遠(yuǎn)，王云宏，段建華

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

導(dǎo)水裂隙帶發(fā)育高度的微震監(jiān)測(cè)研究

叢 森1，程建遠(yuǎn)2，王云宏2，段建華2

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

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

微震監(jiān)測(cè)；微震事件；導(dǎo)水裂隙帶；數(shù)值模擬

0 引言

煤礦開(kāi)采過(guò)程中，采掘活動(dòng)引起采掘工作面周圍巖體應(yīng)力的重新分布，圍巖發(fā)生變形和破壞。微震監(jiān)測(cè)技術(shù)是利用煤巖受力變形和破壞時(shí)產(chǎn)生的地震波來(lái)監(jiān)測(cè)煤巖穩(wěn)定性的技術(shù)方法[1]，通過(guò)在監(jiān)測(cè)區(qū)域不同方位設(shè)置傳感器，對(duì)微震事件進(jìn)行實(shí)時(shí)記錄，然后計(jì)算微震發(fā)生的時(shí)間和位置，進(jìn)而推斷煤巖體的破壞情況。微震事件記錄的是煤巖體變形和破壞逐漸發(fā)生、發(fā)展的過(guò)程，利用微震監(jiān)測(cè)就可以記錄和描述導(dǎo)水裂隙帶形成和發(fā)展的過(guò)程，為導(dǎo)水裂隙帶高度的確定提供依據(jù)。

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

微震監(jiān)測(cè)技術(shù)自21世紀(jì)初開(kāi)始在我國(guó)礦山水害防治、沖擊地壓防治等方面開(kāi)展了廣泛的應(yīng)用[9-13]，具有連續(xù)、實(shí)時(shí)、三維和動(dòng)態(tài)觀測(cè)的特點(diǎn)，可實(shí)現(xiàn)對(duì)整個(gè)工作面導(dǎo)水裂隙帶發(fā)育的長(zhǎng)期動(dòng)態(tài)監(jiān)測(cè)，因此微震監(jiān)測(cè)技術(shù)是導(dǎo)水裂隙帶發(fā)育高度探測(cè)的有力手段。本文采用中煤科工集團(tuán)西安研究院自主研發(fā)的防爆型YTZ3微震監(jiān)測(cè)系統(tǒng)，在某礦首采面開(kāi)展了導(dǎo)水裂隙帶發(fā)育高度的現(xiàn)場(chǎng)監(jiān)測(cè)，并結(jié)合數(shù)值模擬對(duì)導(dǎo)水裂隙帶的發(fā)育高度進(jìn)行了綜合研究。

1 工作面概況及監(jiān)測(cè)系統(tǒng)布置

1.1 工作面概況

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

1.2 YTZ3微震監(jiān)測(cè)系統(tǒng)布置

以往微震監(jiān)測(cè)系統(tǒng)在煤礦上應(yīng)用均以井下部署為主，但由于井下空間有限，觀測(cè)系統(tǒng)幾乎在一個(gè)平面內(nèi)，導(dǎo)致震源定位在深度方向存在一定誤差，這給微震監(jiān)測(cè)技術(shù)在頂板裂隙帶發(fā)育高度監(jiān)測(cè)上的應(yīng)用帶來(lái)了一定的困難。為提高微震監(jiān)測(cè)垂直方向的定位精度，監(jiān)測(cè)臺(tái)網(wǎng)應(yīng)該在已有平面布置的基礎(chǔ)上向立體監(jiān)測(cè)方向發(fā)展，即形成立體的微震監(jiān)測(cè)臺(tái)網(wǎng)。本次監(jiān)測(cè)采用中煤科工集團(tuán)西安研究院自主研發(fā)的YTZ3微震監(jiān)測(cè)系統(tǒng)，在地面和井下聯(lián)合布置監(jiān)測(cè)設(shè)備，形成立體的井-地聯(lián)合微震監(jiān)測(cè)系統(tǒng)，提高微震監(jiān)測(cè)系統(tǒng)在垂直方向的定位精度。

本次監(jiān)測(cè)共設(shè)置測(cè)點(diǎn)14個(gè)，在工作面輔助運(yùn)輸平巷和回風(fēng)平巷頂板各布置3個(gè)微震監(jiān)測(cè)點(diǎn)，地面布置8個(gè)微震監(jiān)測(cè)點(diǎn)。監(jiān)測(cè)系統(tǒng)布置如圖1所示，圖1中紅色點(diǎn)代表井下監(jiān)測(cè)點(diǎn)，綠色點(diǎn)代表地面監(jiān)測(cè)點(diǎn)。

圖1 監(jiān)測(cè)系統(tǒng)布置圖

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

圖2 標(biāo)定炮波形圖

圖3 微震事件波形圖

2 微震監(jiān)測(cè)結(jié)果及分析

2.1 垂直方向微震事件數(shù)量特征

在2016年3月16日～2016年4月26日期間，共接收到有效微震事件8 547個(gè)，微震事件波形如圖3所示，前八道為地面記錄，后六道為井下記錄。對(duì)監(jiān)測(cè)到的微震事件按深度每隔10 m統(tǒng)計(jì)微震事件個(gè)數(shù)，圖4為在各深度區(qū)間內(nèi)對(duì)應(yīng)的微震事件個(gè)數(shù)柱狀圖，煤層底板標(biāo)為+100 m。

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

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

2.2 微震事件空間分布特征

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

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

圖5 微震事件定位結(jié)果3D顯示圖

圖6 微震事件定位結(jié)果XY平面圖

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

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

3 導(dǎo)水裂隙帶發(fā)育高度的數(shù)值模擬

3.1 數(shù)值計(jì)算模型及參數(shù)

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

圖7 微震事件定位結(jié)果走向XZ剖面圖

圖8 三維地質(zhì)模型圖

3.2 模擬結(jié)果分析

模擬工作面回采長(zhǎng)度為360 m，采高4 m，工作面每次推進(jìn)10 m，直到推進(jìn)至360 m。計(jì)算結(jié)果處理過(guò)程中，將采空區(qū)上方相互貫通的拉應(yīng)力塑性區(qū)作為導(dǎo)水裂隙帶，覆巖塑性區(qū)發(fā)育的最高點(diǎn)作為導(dǎo)水裂隙帶發(fā)育的最大高度[14]?；夭山Y(jié)束后，沿工作面走向按一定距離截取工作面傾斜剖面，觀察其塑性區(qū)發(fā)育的最高點(diǎn)。根據(jù)模擬結(jié)果可知，隨著工作面的推進(jìn)導(dǎo)水裂隙帶發(fā)育高度逐漸增加，當(dāng)工作面推進(jìn)至一定位置時(shí)發(fā)育到最高點(diǎn)，之后保持不變，但在停采線后方一定距離內(nèi)呈下降趨勢(shì)，導(dǎo)水裂隙帶發(fā)育高度走勢(shì)如圖9所示，數(shù)值模擬導(dǎo)水裂隙帶發(fā)育的最大高度為108 m。

導(dǎo)水裂隙帶的發(fā)育具有一定的時(shí)間滯后性和距離滯后性(時(shí)間滯后：某一位置導(dǎo)水裂隙帶的發(fā)育要經(jīng)過(guò)一段時(shí)間后才會(huì)趨于穩(wěn)定；距離滯后：某一位置導(dǎo)水裂隙帶發(fā)育到最大高度時(shí)工作面已向前推進(jìn)了一段距離)。因此，在工作面回采結(jié)束后，導(dǎo)致工作面后方一定距離內(nèi)的頂板裂隙沒(méi)有發(fā)育至最高點(diǎn)，即整個(gè)工作面導(dǎo)水裂隙帶的發(fā)育呈現(xiàn)出“發(fā)生→發(fā)展→最高→穩(wěn)定→回降”的發(fā)育規(guī)律。

表1 煤巖體力學(xué)參數(shù)表

圖9 導(dǎo)水裂隙帶發(fā)育高度統(tǒng)計(jì)

4 結(jié)論

根據(jù)上述對(duì)微震事件的特征分析，結(jié)合數(shù)值模擬，可有如下結(jié)論。

1)井-地聯(lián)合微震監(jiān)測(cè)系統(tǒng)形成的立體監(jiān)測(cè)網(wǎng)絡(luò)具有較高的定位精度，能夠滿足對(duì)導(dǎo)水裂隙帶發(fā)育高度的監(jiān)測(cè)要求。

2)基于礦井開(kāi)采的地質(zhì)條件，系統(tǒng)分析微震監(jiān)測(cè)結(jié)果，綜合研究研究得到某礦首采面導(dǎo)水裂隙帶頂界面發(fā)育的最大高度為標(biāo)高+220 m，即導(dǎo)水裂隙帶的發(fā)育高度為110 m。

3)數(shù)值模擬得到某礦首采面導(dǎo)水裂隙帶發(fā)育的最大高度為108 m，對(duì)數(shù)值模擬結(jié)果與微震監(jiān)測(cè)結(jié)果進(jìn)行對(duì)比驗(yàn)證，二者得出的導(dǎo)水裂隙帶發(fā)育最大高度基本一致，微震監(jiān)測(cè)結(jié)果具有較高的可靠性。

4)研究成果為導(dǎo)水裂隙帶的探查提供了一種新的、有效的監(jiān)測(cè)方法。

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

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

[3] 國(guó)家煤炭工業(yè)局. 建筑物、水體、鐵路及主要井巷煤柱留設(shè)與壓煤開(kāi)采規(guī)程[M]. 北京: 煤炭工業(yè)出版社,2000．

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

[5] 張玉軍, 康永華. 覆巖破壞規(guī)律探測(cè)技術(shù)的發(fā)展及評(píng)價(jià)[J]. 煤礦開(kāi)采, 2005, 10(2): 10-12．

[6] 任奮華, 蔡美峰, 來(lái)興平, 等. 采空區(qū)覆巖破壞高度監(jiān)測(cè)分析[J]. 北京科技大學(xué)學(xué)報(bào), 2004, 26(2):40-45.

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

[8] 范志勝, 程敏. 應(yīng)用井下仰孔分段注水觀測(cè)導(dǎo)水?dāng)嗔褞Ц叨萚J]. 礦業(yè)安全與環(huán)保, 2010: 37(6): 75-77.

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

[10] 張興民, 于克君, 席京德, 等. 微地震技術(shù)在煤礦“兩帶”監(jiān)測(cè)領(lǐng)域的研究與應(yīng)用[J]. 煤炭學(xué)報(bào), 2000,25(6): 566-570.

[11] 姜福興, XUN Luo．微震監(jiān)測(cè)技術(shù)在礦井巖層破裂監(jiān)測(cè)中的應(yīng)用[J]．巖土工程學(xué)報(bào), 2002, 24(2):147-149.

[12] 陸菜平, 竇林名, 曹安業(yè), 等. 深部高應(yīng)力集中區(qū)域礦震活動(dòng)規(guī)律研究[J]. 巖石力學(xué)與工程學(xué)報(bào),2008, 27(11): 2302-2308.

[13] 夏永學(xué), 康立軍, 齊慶新, 等. 基于微震監(jiān)測(cè)的5個(gè)指標(biāo)及其在沖擊地壓預(yù)測(cè)中的應(yīng)用[J]. 煤炭學(xué)報(bào), 2010, 35(12):2012-2016.

[14] 韓軍, 張宏偉, 高照宇, 等. 巨厚煤層軟弱覆巖分層綜放開(kāi)采覆巖破壞高度研究[J]. 采礦與安全工程學(xué)報(bào), 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

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

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

TD76

A

1004-4051(2017)03-0126-06