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

董杰 尹士獻(xiàn)

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

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

中圖分類號：TD822 文獻(xiàn)標(biāo)識碼：A 文章編號：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。可采走向長796.66m。可采范圍內(nèi)采面最小傾斜長度為178.2m，最大傾斜長度為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)簡單，煤層厚度為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 觀測點(diǎn)設(shè)置及觀測數(shù)據(jù)

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

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

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

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

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

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

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

實(shí)際工程中只測量了下沉值，那么就要根據(jù)實(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ì)算的時間，a；c表示采深、巖性系數(shù)，c的取值如圖3。

考慮時間因素、采深因素和巖性因素，必須按式（1）計(jì)算。在擬合中按式（1）計(jì)算下沉值，考慮時間因素：最后一次測量的時間為2017年12月18日，此時間節(jié)點(diǎn)與每一塊段回采時間之差，即為時間t的值；在初次擬合時，時間、巖性系數(shù)c取1.0～3.0，把這兩個參數(shù)代入式（2）中，計(jì)算綜合巖性系數(shù)c，再代入式（1）中計(jì)算下沉值，按照這樣方法，消除了時間、巖性因素引起的預(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對預(yù)計(jì)下沉值的影響可用圖4說明。當(dāng)拐點(diǎn)偏移距為正值時，有效開采區(qū)域減小，按照拐點(diǎn)偏移距s的大小，有效開采區(qū)域從實(shí)際開采邊界向采空區(qū)偏移，預(yù)計(jì)時按一區(qū)域面積進(jìn)行計(jì)算；拐點(diǎn)偏移距為零時，實(shí)際開采區(qū)域與預(yù)計(jì)所采用的面積相同，預(yù)計(jì)時按二區(qū)域面積計(jì)算；拐點(diǎn)偏移距為負(fù)時，按照拐點(diǎn)偏移距s的大小，有效開采區(qū)域從實(shí)際開采邊界向?qū)嶓w煤一側(cè)偏移，有效開采區(qū)域增大，預(yù)計(jì)時按三區(qū)域面積進(jìn)行計(jì)算[2，3]。

3.3 最大下沉偏角[θ0]

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

3.4 開采影響傳播角[θ]

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

4 擬合方法

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

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

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

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

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

5.2 反演擬合結(jié)果

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

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

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

6 結(jié)論

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

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

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

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

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

[S=0.024～0.049H] （6）

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

參考文獻(xiàn)：

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

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

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

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

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

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

[7]孟磊，王建坤.高瓦斯特厚煤層采動影響下巷道快速掘進(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é)報（自然科學(xué)版），2012（5）：512-517.