盧繼鷹　田亞坤　劉　永　周薛淼

(1.南華大學核資源工程學院；2.南華大學研究生處)

?

基于概率積分法的某煤礦地表塌陷區(qū)范圍預(yù)測*

盧繼鷹1田亞坤1劉永2周薛淼1

(1.南華大學核資源工程學院；2.南華大學研究生處)

為了科學合理地預(yù)測地表沉陷范圍，結(jié)合礦山開采沉陷理論，應(yīng)用概率積分法，以某煤礦周圍的實測數(shù)據(jù)為例，建立基本數(shù)學預(yù)測模型，對該礦山塌陷區(qū)范圍進行預(yù)測。將預(yù)測結(jié)果與實測結(jié)果進行對比分析，預(yù)測結(jié)果基本覆蓋全部實測結(jié)果，并有一定的安全冗余，故概率預(yù)測法能在一定程度上滿足沉陷區(qū)預(yù)測要求，為塌陷區(qū)預(yù)測提供科學依據(jù)，對塌陷區(qū)災(zāi)害防治具有重要意義。

煤礦采空區(qū)塌陷影響范圍傳統(tǒng)經(jīng)驗法概率積分法數(shù)值模擬

礦產(chǎn)資源開采后遺留形成的空洞即為采空區(qū)[1]，也是開采作業(yè)引發(fā)的圍巖變形、失穩(wěn)，進而產(chǎn)生位移、開裂、破碎、垮落，甚至上覆巖土層整體彎曲、下沉，引起地表變形和破壞的區(qū)域[2-3]，是地面塌陷和地裂縫形成的主要原因[4-5]。采空區(qū)引發(fā)的地表變形分為連續(xù)變形和不連續(xù)變形兩種，對地表形態(tài)、土地資源及農(nóng)田、地面建筑、鐵路、公路、地下水、地表水、滑坡及井田內(nèi)的植被等造成較大的影響和破壞，嚴重影響采空區(qū)上腹及其周圍環(huán)境內(nèi)人們正常的生產(chǎn)和生活。采空區(qū)塌陷還導(dǎo)致礦震頻發(fā)：據(jù)中國地震臺網(wǎng)中心監(jiān)測發(fā)布的消息顯示，2012年2月5日，發(fā)生由煤礦采空區(qū)塌陷引發(fā)的3.0級地震；2015年1月4日上午，榆陽區(qū)麻黃梁鎮(zhèn)二墩煤礦采空區(qū)發(fā)生2.5級塌陷地震，塌陷面積約8.7hm2；2014年4月21日，阿拉善盟阿左旗敖倫布拉格鎮(zhèn)境內(nèi)萬泉礦業(yè)廢棄采空區(qū)發(fā)生一起塌陷事故，塌方區(qū)域面積約400m2，最深處約47m，致居住在此的祖孫3人失蹤，一處房屋被埋壓。

為此，確定采空區(qū)地面塌陷的影響范圍，對合理規(guī)劃周圍居民的正常生產(chǎn)生活尤為重要。本文以湖南某煤礦為實例，采用概率積分法對該礦采空區(qū)地面塌陷的影響范圍進行預(yù)測[6]，期望能準確地描述采空區(qū)地面塌陷范圍[7]，以保證人民群眾的生命財產(chǎn)安全。

1　工程概況

該礦井建于1959年，當時分-20，-55，-90m3個水平，至1989年底，-90m水平以上煤層基本采完，僅剩難以利用的儲量和保安煤柱。1990年進行技術(shù)改造，新增-150m和-200m2個水平。1993年設(shè)計新開主斜井，1997年新主斜井落底于-200m，同年投入生產(chǎn)，主采28～31線-200m以上4、6、7，3個煤層。至1999年底，-150m標高以上煤層基本采完，2004年底，-200m以上煤層已基本采完，開拓工程已達-270m。2004年以后礦山主要開采-200～-300m水平，開采方案為-300m水平暗主斜井延深開拓，暗主斜井口落平標高-200m。

區(qū)內(nèi)地形為丘陵，起伏不大，地勢總體西高東低，丘陵間發(fā)育沖溝。褶斷剝蝕中丘地貌及侵蝕沖溝地貌，最高點海拔178m，最低點海拔85m，自然坡角15°～28°，植被較發(fā)育，溝谷一般為40～150m，海拔80～150m，稻田耕作區(qū)相對標高4～12m，地形由南向北傾斜。根據(jù)調(diào)研結(jié)果，在平面圖中確定采空區(qū)上覆區(qū)域開采沉陷危害的分布位置，繪制了基于實測的采空區(qū)地表塌陷范圍圖，如圖1所示。

圖1　實測的采空區(qū)地表塌陷范圍

2　概率積分法圈定影響范圍

概率積分法是把巖層看作由大量松散的顆粒體介質(zhì)組成，通過隨機介質(zhì)理論，把巖層移動看作一種服從統(tǒng)計規(guī)律的隨機過程來研究巖層與地表移動[8-9]。按隨機介質(zhì)理論，單元開采引起的地表單元盆地呈正態(tài)分布，且與概率密度的分布一致。整個開采引起的下沉剖面方程可以表示為概率密度函數(shù)的積分公式。

2.1概率積分法的基本原理

單元開采引起地表任意點(x，y)的最終下沉值為：

(1)

(2)

式中，m為采高；α為煤層傾角；r為主要影響半徑；H為地面上待計算點(x,y)與煤層上微元點(s,t)的標高差；q為地表下沉系數(shù)；d為拐點偏移距離。

沿x及y方向的水平位移：

(3)

(4)

地表任意一點沿?方向的傾斜值為：

(5)

地表任意一點水平位移為：

(6)

地表任意一點水平變形：

(7)

(8)

地表任意一點沿?方向的曲率為：

(9)

(10)

2.2概率積分法參數(shù)選定

在使用概率積分法對地表移動變形預(yù)測過程中，需要確定的參數(shù)有[10]：下沉系數(shù)q；主要影響角正切tanβ；拐點偏移距離d；水平偏動系數(shù)b。

2.2.1下沉系數(shù)(q)

在充分采動條件下，下沉系數(shù)q由地表出現(xiàn)的最大下沉值[11]與平均開采厚度之比來確定：

(11)

2.2.2主要影響角正切值(tanβ)

主要影響角正切tanβ是沿走向主斷面上開采深度與其主要影響半徑的比值。影響角正切值的大小與煤層上方覆巖的軟硬有關(guān)，巖體越軟、tanβ值越大；巖體巖性越硬、tanβ值越小。我國礦山礦區(qū)的影響正切tanβ值為1.0～3.8，常見的為1.4～2.6。

2.2.3拐點偏移距(d)

2.2.4水平移動系數(shù)(b)

當采動完成后形成移動盆地，沿走向的主斷面上地表水平移動最大值與地表下沉最大值的比值為水平移動系數(shù)，即

(12)

本次預(yù)計參數(shù)采用2006年由地質(zhì)隊完成的地質(zhì)環(huán)境報告里的沉陷值的應(yīng)用特征點法進行反演得到的，見表1。

表1　各參數(shù)的計算結(jié)果

2.3影響范圍確定(r)

將表1中各參數(shù)代入式(2)，得到采空區(qū)塌陷區(qū)影響半徑值，見表2。

將表2所得的采空區(qū)塌陷區(qū)影響半徑值連成一個閉合圈，最終得到用概率積分法圈定的地表塌陷區(qū)影響范圍，如圖2所示。

圖2　概率積分法圈定塌陷區(qū)影響范圍

由圖2可知，概率積分法的預(yù)測范圍基本覆蓋了實測范圍，對于預(yù)測沉陷范圍的準確性較高，雖然預(yù)測范圍比實測范圍大，但為沉陷區(qū)安全管理提供了安全冗余。

表2　各邊界點的影響半徑計算結(jié)果

3　結(jié)　論

(1)以某煤礦為例，采用概率積分法對該礦周圍采空區(qū)沉陷區(qū)進行預(yù)測，將預(yù)測結(jié)果與實測結(jié)果對比發(fā)現(xiàn)，預(yù)測結(jié)果基本覆蓋全部實測結(jié)果，并有一定的安全冗余，故概率預(yù)測法能在一定程度上滿足沉陷區(qū)預(yù)測要求。

(2)以概率積分法進行沉陷范圍預(yù)測形式簡單，計算較快，便于節(jié)省預(yù)測時間和費用。

(3)概率積分法參數(shù)明確，采用實測數(shù)據(jù)作為參數(shù)依據(jù),既能得到精度較高且較可靠的預(yù)測參數(shù)及預(yù)測結(jié)果,又便于研究認識開采沉陷規(guī)律，給沉陷區(qū)范圍預(yù)測提供了科學依據(jù)。

[1]彭林軍,趙曉東,李術(shù)才.深部開采地表沉陷規(guī)律模擬研究[J].巖土力學,2011(6):1910-1914.

[2]楊梅忠,任秀芳,于遠祥.概率積分法在煤礦采空區(qū)地表變形動態(tài)評價中的應(yīng)用[J].西安科技大學學報,2007(1):39-42.

[3]吳韶艷.煤礦采空沉陷原理及地表移動變形預(yù)測研究[D].西安：長安大學,2011.

[4]顧葉,宋振柏,張勝偉.基于概率積分法的開采沉陷預(yù)計研究[J].山東理工大學學報：自然科學版,2011(1):33-36.

[5]張應(yīng)昌,孫亞峰,何堯.基于概率積分法的開采沉陷預(yù)計分析[J].山西建筑,2009(20):119-120.

[6]趙啟峰,孟祥瑞,趙光明,等.厚松散層下開采沉陷預(yù)測模型的研究及應(yīng)用[J].中國煤炭,2008(5):36-38,43.

[7]查劍鋒,郭廣禮,趙海濤,等.概率積分法修正體系現(xiàn)狀及發(fā)展展望[J].金屬礦山,2008(1):15-18.

[8]卞和方,楊化超,張書畢.概率積分法預(yù)計參數(shù)的智能優(yōu)化選擇方法研究[J].采礦與安全工程學報,2013(3):385-389.

[9]張正興,趙愛軍,李威,等.基于空間分析的概率積分法在預(yù)測多煤層采空區(qū)塌陷中的應(yīng)用——以青海塔妥煤礦為例[J].中國地質(zhì)災(zāi)害與防治學報,2012(4):94-98.

[10]郭文兵,鄧喀中,鄒友峰.概率積分法預(yù)計參數(shù)選取的神經(jīng)網(wǎng)絡(luò)模型[J].中國礦業(yè)大學學報,2004(3):88-92.

[11]姚琦,馮濤,李石林,等.基于概率積分法的煤礦“三下”開采沉陷預(yù)計[J].煤礦安全,2012(7):188-190,193.

PredictionoftheEarthSurfaceSubsidenceScopeofaCoalMineBasedontheProbabilityIntegralMethod

LuJiying1TianYakun1LiuYong2ZhouXuemiao1

(1.SchoolofNuclearResourcesEngineering,UniversityofSouthChina;2.DepartmentofGraduateStudent,UniversityofSouthChina)

Inordertocarryoutscientificandreasonablepredictionoftheearthsurfacesubsidencescope,combingwiththeminingsubsidencetheoryandprobabilityintegralmethod,thebasicmathematicalpredictionmodelofminesubsidenceareaisestablishedtocarryoutthepredictionofthesurfcesubsidencerange.Combingwiththemeasureddataaroundacoalmine,theprobabilityintegralmethodisusedtopredicttheminegoaf.Thepredictionresultsiscomparedwiththemeasuredresults,thecomparisionanalysisresultsshowthatthemeasuredresultsisalmostcoveredbythepredictionresults,andacertaindegreeofsafetyredundancyisexisted.So,theprobabilityintegralmethodcanmeettherequirementsofsubsidencepredictioninacertaindegree,besidesthat,italsocanprovideascientificbasisforpredictingtheminesubsidenceareaandthedisasterpreventionofminesubsidencearea.

Goafofcoalmine,Traditionalempiricalmethod,Probabilityintegralmethod,Numericalsimulation,Subsidencescope

劉永(1971—)，男，教授，處長，博士。

2016-03-09)

*湖南省研究生科研創(chuàng)新項目(編號：2015SCX13)。

盧繼鷹(1990—),男，碩士研究生,421001 湖南省衡陽市常勝西路28號。