深部非充分采動(dòng)巖移參數(shù)反演擬合研究

董杰 尹士獻(xiàn)

摘 要：平頂山礦區(qū)香山礦目前采用的巖移參數(shù)是15年前礦井在回采淺部煤層（采深200～300m）時(shí)通過實(shí)測(cè)獲得的。隨著回采深度的增加，巖移參數(shù)也會(huì)發(fā)生變化。本文通過香山礦己16-17-24061工作面實(shí)測(cè)的下沉值，運(yùn)用反演擬合的方法，擬合出深部煤層非充分采動(dòng)（采深750m左右）的巖移參數(shù)，并與淺部煤層回采巖移參數(shù)相比較，得出兩者之間的差別，為相似地質(zhì)條件下深部非充分采動(dòng)巖移參數(shù)的選取提供借鑒。

關(guān)鍵詞：巖移參數(shù)；反演擬合；誤差分析；開采沉陷

中圖分類號(hào)：TD822 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1003-5168（2018）22-0078-05

Study on Inversion Fitting of Rock Movement Parameters for

Deep Insufficient Mining

DONG Jie1 YIN Shixian 2

（1.Xiangshan Coal Mining of Pingdingshan Tian'an Coal Industry Co.， Ltd.，Pingdingshan Henan 467048；2. Safety Technology Training Institute of Henan Polytechnic University， Jiaozuo Henan 454003）

Abstract： The rock shift parameters used in Xiangshan mine in Pingdingshan mining area are obtained 15 years ago when the mine was measured in the shallow coal seam （200m to 300m）. With the increase of mining depth， the parameters of rock movement will also change. Through the inversion fitting method， the rock shift parameters of the deep coal seam were not fully produced （about 750m in the mining depth） by the inversion fitting method， and the difference between the two in the shallow coal seams was compared with that of the shallow coal seam， which was the deep non sufficient mining rock under the similar geological condition. The selection of moving parameters could be used for reference.

Keywords： rock movement parameters；inversion fitting；error analysis；mining subsidence

1 回采工作面地質(zhì)概況及回采情況

己16-17-24061回采工作面位于己四采區(qū)西翼，東起回風(fēng)下山，西止于59勘探線西287m處。南部為西翼風(fēng)井保護(hù)煤柱，北部為原生煤體。風(fēng)巷標(biāo)高-661.47～-700.66m，平均-697.14m。機(jī)巷標(biāo)高-730.16～-744.05m，平均-737.83m。地面標(biāo)高+167m。可采走向長(zhǎng)796.66m?？刹煞秶鷥?nèi)采面最小傾斜長(zhǎng)度為178.2m，最大傾斜長(zhǎng)度為179.3m。自2016年4月開始回采，2017年3月回采結(jié)束。每月回采煤量、回采進(jìn)度及回采平均厚度見圖1。

根據(jù)己16-17-24061工作面掘進(jìn)揭露資料統(tǒng)計(jì)，己16-17-24061工作面煤層結(jié)構(gòu)簡(jiǎn)單，煤層厚度為4.4～6.7m，平均厚度約為5.77m。其中，風(fēng)巷（有效可采范圍內(nèi)）附近厚度為4.7～6.7m，平均厚度約為5.79m；機(jī)巷（有效可采范圍內(nèi)）附近厚度為4.4～6.3m，平均厚度約為5.75m。根據(jù)機(jī)巷、風(fēng)巷掘進(jìn)期間探煤資料計(jì)算，該工作面回采范圍內(nèi)己16-17煤層可采指數(shù)為1.00（>0.95），變異系數(shù)為10.8%（<25%），屬于穩(wěn)定煤層。

2 觀測(cè)點(diǎn)設(shè)置及觀測(cè)數(shù)據(jù)

2.1 觀測(cè)點(diǎn)設(shè)置

己16-17-24061工作面緊臨鄭石高速和碎石路，為設(shè)置的觀測(cè)站提供了便利條件，同時(shí)減少了設(shè)置觀測(cè)站的費(fèi)用。沿鄭石高速設(shè)置了6個(gè)觀測(cè)點(diǎn)，分別為：2、4、6、8、10和11號(hào)觀測(cè)點(diǎn)；沿碎石路設(shè)置了6個(gè)觀測(cè)點(diǎn)，分別為：1、3、5、7、9和11號(hào)觀測(cè)點(diǎn)。地表觀測(cè)站布置具體見圖2。

沿鄭石高速的觀測(cè)線長(zhǎng)為486.8m，沿碎石路的觀測(cè)線長(zhǎng)為717.7m。各測(cè)點(diǎn)之間的距離見表1和表2。

2.2 兩條觀測(cè)線觀測(cè)數(shù)據(jù)

在沒有回采之前，于2016年1月20日進(jìn)行首次基準(zhǔn)測(cè)量，后又分別于2016年7月28日進(jìn)行第二次測(cè)量，2016年10月30日進(jìn)行第三次測(cè)量，2017年4月30日進(jìn)行第四次測(cè)量，2017年12月18日進(jìn)行第五次測(cè)量。從2016年1月至2017年12月共進(jìn)行了5次地表觀測(cè)，歷時(shí)2年。根據(jù)表3各觀測(cè)點(diǎn)的實(shí)測(cè)值可以擬合出深部回采時(shí)的巖移參數(shù)。

3 影響下沉值的巖移參數(shù)

3.1 時(shí)間、巖性系數(shù)c、下沉系數(shù)[η]、主要影響角[β]、

實(shí)際工程中只測(cè)量了下沉值，那么就要根據(jù)實(shí)測(cè)的下沉來擬合巖移參數(shù)。擬合時(shí)，下沉值按照公式（1）進(jìn)行計(jì)算[1]：

[wx，y，t=wmaxr2Cpe-πx2+y2r2] （1）

[C=1-e-ct] （2）

式中，[wmax]表示最大下沉值，mm；[r]表示地表主要影響半徑，m；p表示回采區(qū)域，m2；t表示開采至計(jì)算的時(shí)間，a；c表示采深、巖性系數(shù)，c的取值如圖3。

考慮時(shí)間因素、采深因素和巖性因素，必須按式（1）計(jì)算。在擬合中按式（1）計(jì)算下沉值，考慮時(shí)間因素：最后一次測(cè)量的時(shí)間為2017年12月18日，此時(shí)間節(jié)點(diǎn)與每一塊段回采時(shí)間之差，即為時(shí)間t的值；在初次擬合時(shí)，時(shí)間、巖性系數(shù)c取1.0～3.0，把這兩個(gè)參數(shù)代入式（2）中，計(jì)算綜合巖性系數(shù)c，再代入式（1）中計(jì)算下沉值，按照這樣方法，消除了時(shí)間、巖性因素引起的預(yù)計(jì)誤差[2]。

除此之外，還有4種因素影響下沉預(yù)計(jì)精度。分別為：最大下沉值[wmax]、回采區(qū)域p、最大下沉偏移角和主要影響半徑r。其中最大下沉值[wmax]的計(jì)算公式為：

[wmax=ηmcosα] （3）

式中，[η]表示下沉系數(shù)；[m]表示回采厚度，[m]；[α]表示煤層傾角，°；

主要影響半徑r的計(jì)算公式：

[r=Htanβ] （4）

式中，[H]表示采深，m；[β]表示主要影響角，°；

從式（3）可知，在回采厚度與煤層傾角一定的情況下，最大下沉值的大小只與下沉系數(shù)有關(guān)；從式（4）可知，在埋深一定的情況下，半徑大小只與主要影響角有關(guān)。

3.2 拐點(diǎn)偏移距s

拐點(diǎn)偏移距s對(duì)預(yù)計(jì)下沉值的影響可用圖4說明。當(dāng)拐點(diǎn)偏移距為正值時(shí)，有效開采區(qū)域減小，按照拐點(diǎn)偏移距s的大小，有效開采區(qū)域從實(shí)際開采邊界向采空區(qū)偏移，預(yù)計(jì)時(shí)按一區(qū)域面積進(jìn)行計(jì)算；拐點(diǎn)偏移距為零時(shí)，實(shí)際開采區(qū)域與預(yù)計(jì)所采用的面積相同，預(yù)計(jì)時(shí)按二區(qū)域面積計(jì)算；拐點(diǎn)偏移距為負(fù)時(shí)，按照拐點(diǎn)偏移距s的大小，有效開采區(qū)域從實(shí)際開采邊界向?qū)嶓w煤一側(cè)偏移，有效開采區(qū)域增大，預(yù)計(jì)時(shí)按三區(qū)域面積進(jìn)行計(jì)算[2，3]。

3.3 最大下沉偏角[θ0]

最大下沉偏角是指采空區(qū)中心與下沉盆地中心在地表投影的連線與水平線在下山方向的夾角。最大下沉偏角對(duì)地表下沉盆地有較大影響。在采深相同的情況下，下沉偏移越小，向下山方向偏移的距離越遠(yuǎn)，反之越近，如圖5所示。因此，需要對(duì)最大下沉偏移角進(jìn)行反演擬合[2，4]。

3.4 開采影響傳播角[θ]

開采影響傳播角的定義為：在傾向主剖面上，最大下沉值一半的點(diǎn)在地表的投影與有效開采邊界點(diǎn)的連線與水平線在下山方向的夾角，如圖6所示。由此可知，只要求得了拐點(diǎn)偏移距和傾斜主剖面最大下沉值一半的點(diǎn)位，就可以求得開采影響傳播角。在工程應(yīng)用中，只要知道這三者中的任意2個(gè)，即可求得另一個(gè)巖移參數(shù)[2，5]。

4 擬合方法

通過以上分析可知，有4個(gè)巖移參數(shù)影響下沉預(yù)計(jì)精度，這4個(gè)巖移參數(shù)依次為：拐點(diǎn)偏移距s、下沉系數(shù)[η]、主要影響角[β]、最大下沉偏移角[θ0]。開采影響傳播角是通過傾斜主剖面上的下沉值和拐點(diǎn)偏移距求得。

第一步，擬合拐點(diǎn)偏移距s，其他巖移參數(shù)固定不變，按照香山礦給出的巖移參數(shù)作為擬合的初始值進(jìn)行計(jì)算。拐點(diǎn)偏移距s逐漸變化，當(dāng)擬合達(dá)到最小誤差時(shí)，此時(shí)的拐點(diǎn)偏移距s即為所求，程序框圖如圖7所示；第二步，擬合下沉系數(shù)[η]，拐點(diǎn)偏移距s采用擬合出的值，其他另外2個(gè)巖移參數(shù)，采用香山礦給出巖移參數(shù)作為擬合的初始值進(jìn)行計(jì)算，逐漸變化下沉系數(shù)[η]，當(dāng)擬合達(dá)到最小誤差時(shí)，此時(shí)的下沉系數(shù)[η]即為所求；第三步，擬合主要影響角[β]，下沉系數(shù)[η]、拐點(diǎn)偏移距s采用擬合出的數(shù)值，其他巖移參數(shù)采用香山礦給出的數(shù)值作為擬合的初始值進(jìn)行計(jì)算，逐漸變化主要影響角[β]，當(dāng)擬合達(dá)到最小誤差時(shí)，此時(shí)的主要影響角[β]即為所求；第四步，擬合最大下沉偏移角[θ0]，其他3個(gè)巖移參數(shù)采用擬合出的數(shù)值，逐漸變化最大下沉偏移角[θ0]，當(dāng)擬合達(dá)到最小誤差時(shí)，此時(shí)的最大下沉偏移角[θ0]即為所求[6]。

以上介紹的為一次完整循環(huán)。第一次循環(huán)結(jié)束后，若精度達(dá)不到工程要求，再按照上述方法進(jìn)行第二次循環(huán)，一直循環(huán)到達(dá)到工程精度要求為止。利用C++語(yǔ)言自主研發(fā)的擬合程序進(jìn)行深部巖移參數(shù)擬合。

5.1 巖移參數(shù)反演擬合初次賦值

擬合中要用到巖移參數(shù)的初次賦值，為了減少循環(huán)次數(shù)，節(jié)省計(jì)算時(shí)間，初次賦值在香山礦提供的巖移參數(shù)基礎(chǔ)上，適當(dāng)上下浮動(dòng)，作為初次賦值的最大值與最小值。如表4所示，其中主要影響角和拐點(diǎn)偏移距分為上山、下山和走向，三個(gè)方向分別賦值擬合。在擬合出了拐點(diǎn)偏移距后，開采影響傳播角根據(jù)傾向主剖面的下沉形態(tài)作圖求得[7]。

5.2 反演擬合結(jié)果

運(yùn)用上述循環(huán)擬合的方法，擬合出香山礦己16-17-24061工作面回采非充分采動(dòng)時(shí)的巖移參數(shù)如表5所示。表5中開采影響傳播角是通過上下山拐點(diǎn)偏移距及傾斜主剖面的形態(tài)，用作圖法求得，限于篇幅在這里不再給出求解過程，只給出求得的結(jié)果[8]。

沿碎石路擬合下沉曲線與實(shí)測(cè)下沉曲線如圖8所示。從圖8中可以直觀地看出，9號(hào)觀測(cè)點(diǎn)和3號(hào)觀測(cè)點(diǎn)誤差較??；沿鄭石高速擬合下沉曲線與實(shí)測(cè)下沉曲線如圖9所示。從圖9中可以直觀地看出，11號(hào)觀測(cè)點(diǎn)和10號(hào)觀測(cè)點(diǎn)誤差較小。

為了更好地驗(yàn)證擬合的精確程度，需要計(jì)算擬合的相對(duì)誤差和絕對(duì)誤差。誤差計(jì)算見表6。從表6中可以看出，9號(hào)測(cè)點(diǎn)絕對(duì)誤差及相對(duì)誤差最小，8號(hào)測(cè)點(diǎn)絕對(duì)誤差最大，11號(hào)測(cè)點(diǎn)相對(duì)誤差最大。最大絕對(duì)誤差不超過4mm，最大相對(duì)誤差不超過26%，擬合精度達(dá)到工程要求[9]。

6 結(jié)論

由于香山礦在測(cè)量中沒有測(cè)量水平移動(dòng)值，所以，其他巖移參數(shù)無法求得，如水平移動(dòng)系數(shù)等。通過反演擬合求得的深部非充分采動(dòng)巖移參數(shù)與淺部非充分采動(dòng)相比，一些巖移參數(shù)有較大差別，如下沉系數(shù)[η]、主要影響角[β]和拐點(diǎn)偏移距s，而另一些巖移參數(shù)變化不大，如時(shí)間、巖性系數(shù)c、最大下沉偏角[θ0]、開采影響傳播角[θ]。筆者總結(jié)的深部回采與淺部回采巖移參數(shù)差別如下，以期為相似地質(zhì)條件下的工程提供借鑒。

①下沉系數(shù)[η]，淺部回采為0.3，深部回采為0.12，那么，若知道淺部非充分采動(dòng)的下沉系數(shù)[η淺]和淺部回采時(shí)的采深[H淺]，就可以用式（5）估算采深為[H深]時(shí)非充分采動(dòng)的下沉系數(shù)[η深]：

[η深=η淺-H深-H淺×2.93×10-4] （5）

②香山礦只給出了走向主要影響角，未給出上山和下山方向，且擬合出的走向主要影響角正切與淺部非充分采動(dòng)的主要影響角正切值相差較大，待以后的工程實(shí)驗(yàn)證。

③根據(jù)擬合的拐點(diǎn)偏移距及采深，可以算出深部非充分采動(dòng)時(shí)拐點(diǎn)偏移距與采深的關(guān)系：

[S=0.024～0.049H] （6）

式中，s為深部非充分采動(dòng)拐點(diǎn)偏移距，m；[H]為回采深度，m。

參考文獻(xiàn)：

[1]余學(xué)義，張恩強(qiáng)，開采損害學(xué)[M].北京：煤炭工業(yè)出版社，2004.

[2]尹士獻(xiàn).構(gòu)造應(yīng)力場(chǎng)與采動(dòng)應(yīng)力場(chǎng)協(xié)同作用下覆巖破壞區(qū)域研究[D].焦作：河南理工大，2015.

[3]國(guó)家煤炭工業(yè)局.建筑物、水體、鐵路及主要井巷煤柱留設(shè)與壓煤開采規(guī)范[M].北京：煤炭工業(yè)出版社，2017.

[4]郝從娜，張子平，田富超.國(guó)家地籍?dāng)?shù)據(jù)庫(kù)的設(shè)計(jì)及關(guān)鍵技術(shù)分析[J].河南理工大學(xué)學(xué)報(bào)（自然科學(xué)版），2009（4）：455-458.

[5]寧永香.程莊煤礦2305工作面開采沉陷觀測(cè)站設(shè)計(jì)[J].河南科技，2017（2）：94-98.

[6]朱沖朋.煤礦開采沉陷預(yù)計(jì)方法研究[J].河南科技，2015（9）：71-72.

[7]孟磊，王建坤.高瓦斯特厚煤層采動(dòng)影響下巷道快速掘進(jìn)研究[J].陜西煤炭，2012（6）：35-39.

[8]冉小軍.試論礦山生產(chǎn)過程中的幾種常見環(huán)境難題[J].河南科技，2014（7）：165.

[9]尹士獻(xiàn)，程志勇，李德海.裂隙帶布置高抽巷實(shí)現(xiàn)巷道快速掘進(jìn)[J].河南理工大學(xué)學(xué)報(bào)（自然科學(xué)版），2012（5）：512-517.