張正興，趙愛(ài)軍，李 威，李小林，孫 瑩，彭 亮

(1.青海省水文地質(zhì)工程地質(zhì)環(huán)境地質(zhì)調(diào)查院，青海西寧 810008;2.青海省環(huán)境地質(zhì)勘查局，青海西寧 810007;3.青海省地質(zhì)環(huán)境監(jiān)測(cè)總站，青海西寧 810003)

0 引言

煤炭是我國(guó)的主要能源資源，對(duì)國(guó)民經(jīng)濟(jì)的發(fā)展起著舉足輕重的作用，但煤礦開(kāi)采引起的地表移動(dòng)和變形將影響位于下沉盆地范圍內(nèi)的建筑物、河流、鐵路、管道及其他構(gòu)筑物。煤層開(kāi)采前對(duì)開(kāi)采沉陷做出預(yù)計(jì)，對(duì)保護(hù)地面建筑物和地面環(huán)境具有十分重要的意義。關(guān)于礦區(qū)開(kāi)采沉陷國(guó)內(nèi)外已進(jìn)行了大量的研究工作，原蘇聯(lián)學(xué)者早在20世紀(jì)30年代就開(kāi)始了相關(guān)研究，我國(guó)則是從20世紀(jì)50年代開(kāi)始的［1-2］。近年來(lái)，由于數(shù)學(xué)、力學(xué)、計(jì)算機(jī)等學(xué)科的新成果不斷地引入到開(kāi)采沉陷預(yù)測(cè)中，形成了許多開(kāi)采沉陷變形預(yù)計(jì)的新方法和新理論，并取得了一定的研究成果。目前我國(guó)用于地面變形預(yù)測(cè)的方法中，概率積分法是最為成熟、應(yīng)用最廣泛的方法之一。但該方法也存在一些缺陷，尤其是對(duì)于地層結(jié)構(gòu)比較復(fù)雜的地區(qū)［3-5］，若以概率積分法為基礎(chǔ)，運(yùn)用一些空間分析軟件，將礦區(qū)不同煤組的破壞作用進(jìn)行疊加，便得到了更加符合實(shí)際的預(yù)測(cè)結(jié)果［6-11］。

1 概率積分法

概率積分法是以正態(tài)分布函數(shù)為影響函數(shù)，用積分式表示地表盆地下沉的方法。這種方法的基礎(chǔ)是隨機(jī)介質(zhì)理論。其數(shù)學(xué)模型為:

下沉(W，單位mm):

傾斜(i，單位mm/m)

曲率(K，單位10-3/m)

水平位移(U，單位mm):

水平變形(ξ，單位mm/m):

式中:r——主要影響半徑(m);

r=H/tanβ

H——采深(m);

tanβ-主要影響角正切;

Wcm、Ucm——地表充分采動(dòng)下的最大下沉值和最大水平移動(dòng)值(mm)。

該方法的預(yù)測(cè)模型相對(duì)簡(jiǎn)單，參數(shù)比較穩(wěn)定，計(jì)算精度也較高。在《建筑物、水體、鐵路及主要井巷煤柱留設(shè)與壓煤開(kāi)采規(guī)程》、《礦山環(huán)境保護(hù)與綜合治理方案技術(shù)規(guī)范》(DZ/T223-2007)等標(biāo)準(zhǔn)中，有明確的規(guī)定［2］。

從其數(shù)學(xué)模型可以看出，該方法適用于連續(xù)的，分布規(guī)整的單一煤層采空區(qū)。但在實(shí)際情況下，煤層分布的規(guī)則性一般都較弱，其采空區(qū)分布凌亂，這樣就導(dǎo)致了地表移動(dòng)預(yù)測(cè)的困難，單靠概率積分法難以解決問(wèn)題。為了使預(yù)測(cè)計(jì)算結(jié)果更加接近實(shí)際，將復(fù)雜煤層開(kāi)采形成的采空區(qū)概化成一系列比較連續(xù)、規(guī)整的采空區(qū)組合體，先運(yùn)用概率積分法預(yù)測(cè)計(jì)算出每一個(gè)采空區(qū)產(chǎn)生的地表移動(dòng)，然后利用一些空間分析軟件，將各個(gè)采空區(qū)產(chǎn)生的地面破壞進(jìn)行疊加，從而得到最終的預(yù)測(cè)結(jié)果。

2 塔妥煤礦

塔妥煤礦位于青海省都蘭縣溝里鄉(xiāng)，自建國(guó)初期就開(kāi)始開(kāi)采，目前年生產(chǎn)能力6×104t。由于長(zhǎng)年的開(kāi)采，淺層地下水已完全疏干。礦區(qū)地層缺失較多，可分為三個(gè)部分:煤系地層;基底地層以及第四系地層。基底地層為石炭系和二疊系，巖性主要為砂巖、板巖、千枚巖;煤系地層巖性主要為侏羅系砂巖、碳質(zhì)頁(yè)巖、頁(yè)巖為主;第四系地層以砂礫石、砂土、角礫為主。

研究區(qū)含煤層分布復(fù)雜，自下而上沉積有一、二、三、四4個(gè)煤組，其中三煤組最為發(fā)育，厚度大，分布廣;二煤組次之;一煤組因受古地理影響，沉積面積較少，僅礦區(qū)東部有沉積;四煤組因成煤后期剝蝕作用的影響，保存的范圍較小，僅在礦區(qū)中部可見(jiàn)。依據(jù)礦山的探礦鉆孔以及塌陷區(qū)勘探鉆孔資料，所揭示的各煤組產(chǎn)狀特征如圖1、表1所示。

圖1 塔妥煤礦地質(zhì)剖面圖Fig.1 The cross-section of Tatuo coal mine

表1 塔妥煤礦各煤組產(chǎn)狀特征表Table 1 The attitude characteristics of Tatuo coal mine's all coal groups

3 地面破壞預(yù)測(cè)

3.1 計(jì)算參數(shù)

在參考前人研究成果的基礎(chǔ)上，根據(jù)以往的鉆孔資料、相關(guān)圖件以及本次勘查的鉆探、物探成果，結(jié)合勘查區(qū)內(nèi)四個(gè)主要含煤組段的形態(tài)特征，得出塔妥煤礦的相關(guān)計(jì)算參數(shù)。其中下沉系數(shù)、主要影響角正切以及開(kāi)采影響傳播角與上覆圍巖的特性關(guān)系密切，根據(jù)當(dāng)?shù)貒鷰r的相關(guān)力學(xué)實(shí)驗(yàn)和經(jīng)驗(yàn)參數(shù)給定，水平移動(dòng)系數(shù)是最大水平移動(dòng)值與最大下沉值的比值，根據(jù)經(jīng)驗(yàn)確定(表2)。

表2 各煤組計(jì)算參數(shù)表Table 2 Counting parameters of every coal group

3.2 預(yù)測(cè)計(jì)算過(guò)程

塔妥煤礦并非單煤層開(kāi)采，其含煤層主要分為四個(gè)組段，地表沉陷也非一個(gè)煤組的開(kāi)采所造成的，而是四個(gè)煤組共同作用的結(jié)果。傳統(tǒng)的概率積分法通常會(huì)將多個(gè)煤層概化為單一含煤層，然后再進(jìn)行地面沉陷預(yù)測(cè)，這樣難以突出主采煤層的影響，同時(shí)計(jì)算結(jié)果的誤差較大。因此在本次預(yù)測(cè)計(jì)算中，為了充分考慮各煤組采空區(qū)對(duì)地表的破壞作用，先運(yùn)用概率積分法計(jì)算出各采空區(qū)對(duì)地表的作用，然后利用GIS軟件將各煤組的預(yù)測(cè)結(jié)果進(jìn)行疊加。

3.3 預(yù)測(cè)結(jié)果

圖2～圖4為運(yùn)用概率積分法計(jì)算出的二煤組采空區(qū)產(chǎn)生的地面破壞結(jié)果。圖2反映了二煤組采空區(qū)產(chǎn)生的地面沉陷值，最大沉陷量為671mm;圖3體現(xiàn)了采空區(qū)產(chǎn)生的地面傾斜程度，其最大曲率值為0.095×10-3/m，出現(xiàn)在離計(jì)算邊界44.24m的位置;圖4反映了采空區(qū)導(dǎo)致地面水平變形值的大小，在距離西部計(jì)算邊界0～16.49m的范圍內(nèi)，地面呈拉伸態(tài)，其余地區(qū)均為壓縮態(tài)。最大拉伸變形為1.97mm/m，出現(xiàn)在西側(cè)計(jì)算邊界處，最大壓縮變形為4.32mm/m，出現(xiàn)在距離東側(cè)計(jì)算邊界53.00m的位置。同理運(yùn)用概率積分法分別計(jì)算預(yù)測(cè)出一煤組、三煤組和四煤組采空區(qū)產(chǎn)生的地面破壞。

圖2 二煤組采空區(qū)地面沉陷預(yù)測(cè)圖Fig.2 The forecasting chart for surface depression of the second coal group

圖3 二煤組采空區(qū)地面曲率預(yù)測(cè)圖Fig.3 The forecasting chart for surface curvature of the second coal group

圖4 二煤組采空區(qū)地面水平變形預(yù)測(cè)圖Fig.4 The forecasting chart for surface horizontal distortion of the second coal group

由二煤組的預(yù)測(cè)結(jié)果可以推出若按照傳統(tǒng)的概率積分法，將塔妥煤礦的所有含煤層概化為單一煤層而進(jìn)行地面沉陷預(yù)測(cè)的結(jié)果。屆時(shí)礦區(qū)中部的大面積區(qū)域都將呈均勻沉降趨勢(shì)，都將達(dá)到地面塌陷的最大值，地面彎曲變形和水平變形也僅發(fā)生在礦區(qū)的邊緣地帶。煤礦開(kāi)采形成的采空區(qū)對(duì)礦區(qū)中部的地表建筑物和道路等影響不大，然而通過(guò)現(xiàn)場(chǎng)的實(shí)際調(diào)查發(fā)現(xiàn)，塔妥煤礦中部地帶的房屋等設(shè)施目前也發(fā)生嚴(yán)重的變形，由此可見(jiàn)按照傳統(tǒng)概率積分法進(jìn)行地面沉陷預(yù)測(cè)的結(jié)果與現(xiàn)實(shí)情況存在較大差異。

運(yùn)用GIS軟件的矢量疊加功能，將各煤組采空區(qū)產(chǎn)生的破壞結(jié)果進(jìn)行疊加，便可得到假定塔妥煤礦中的含煤層完全開(kāi)采時(shí)，所有采空區(qū)造成的地面破壞情況。其地面沉陷預(yù)測(cè)結(jié)果如圖5，其最大沉陷為4002mm，出現(xiàn)在四個(gè)煤組重疊地帶。圖6反映了采空區(qū)致使的地面曲率值，其預(yù)測(cè)最大值為0.88×10-3/m，分布在三煤組邊界附近，呈條帶狀，礦區(qū)地面曲率的大小主要受到三煤組和四煤組的影響。采空區(qū)產(chǎn)生的地面水平變形預(yù)測(cè)值反映在圖7中，由圖可見(jiàn)，采空區(qū)地面水平變形的大小也主要是受到三煤組和四煤組的影響，其最大值為28.03mm/m，出現(xiàn)在距離東側(cè)計(jì)算邊界約124.2m的位置。將此預(yù)測(cè)結(jié)果與目前的實(shí)際房屋損壞調(diào)查相比較，此預(yù)測(cè)結(jié)果更具有實(shí)際意義。

由此可見(jiàn)，塔妥煤礦一、二、三、四煤組的開(kāi)采對(duì)地表都會(huì)有一定的影響，只是影響程度的大小不一。其中三煤組的開(kāi)采對(duì)地表的破壞程度最大，它產(chǎn)生的地面沉陷值、曲率值和水平變形值較其他煤組大，這是由于三煤組的分布范圍廣，煤層厚度較大，它的開(kāi)采會(huì)產(chǎn)生較大的地下采空區(qū)，而且上覆巖層較薄;四煤開(kāi)采組對(duì)地表的破壞程度居于次席，這主要是由于該煤組的采深較小，它的影響之所以比三煤組小，主要是因?yàn)樗姆植挤秶?，煤層的厚度也比三煤組的薄，因此其產(chǎn)生的地表破壞較三煤組小;二煤組和一煤組隨著地下埋深的增加，對(duì)地表的影響逐漸減弱。

圖5 塔妥煤礦地面沉陷預(yù)測(cè)圖Fig.5 The forecasting chart for surface depression of the Tatuo Coal Mine

圖6 塔妥煤礦地面曲率預(yù)測(cè)圖Fig.6 The forecasting chart for surface curvature of the Tatuo Coal Mine

圖7 塔妥煤礦地面水平變形預(yù)測(cè)圖Fig.7 The forecasting chart for surface horizontal distortion of the Tatuo Coal Mine

4 結(jié)論

概率積分法計(jì)算模型簡(jiǎn)單，參數(shù)穩(wěn)定，對(duì)于平原區(qū)的單一煤層，計(jì)算精度較高。但由于其采用的參數(shù)中大多數(shù)以平均值代替，難以客觀完全的反映實(shí)際煤層的特征，從而致使誤差較大，尤其是對(duì)于地層復(fù)雜的地區(qū)，煤層參數(shù)更加難以準(zhǔn)確把握，計(jì)算精度較低。在本文中，對(duì)于地層結(jié)構(gòu)非常復(fù)雜的塔妥煤礦，改變傳統(tǒng)的平均值法，運(yùn)用概率積分法分別計(jì)算出每一個(gè)煤組的地面破壞，然后利用GIS軟件，將其破壞結(jié)果進(jìn)行合理疊加，而不是先將煤層概化成單層煤層之后，再運(yùn)用概率積分法進(jìn)行預(yù)測(cè)計(jì)算。根據(jù)目前的監(jiān)測(cè)資料和煤礦的開(kāi)采進(jìn)程來(lái)看，這樣預(yù)測(cè)方法的計(jì)算結(jié)果離實(shí)際情況又靠近一步，對(duì)該區(qū)地面沉陷的防護(hù)和治理有很好的指導(dǎo)意義。

［1］鄒友峰，鄧喀中，馬偉民.礦山開(kāi)采沉陷工程［M］.北京:中國(guó)礦業(yè)大學(xué)出版社，2003.ZOU Youfeng，DENG Kazhong，MA Weimin.Mining Subsidence Project［M］.Beijing:Mines University China Publishing Company，2003.

［2］建筑物、水體、鐵路及主要井巷煤柱留設(shè)與壓煤開(kāi)采規(guī)程［S］.北京:煤炭工業(yè)出版社，2000.The regulations of supposing coal wedge for building，water，railway together with main tunnel and mining pressing coal［S］.Beijing:Coal Industry Publishing Company，2000.

［3］何國(guó)清，楊倫，等.礦山開(kāi)采沉陷學(xué)［M］.北京:中國(guó)礦業(yè)大學(xué)出版社，1991.HE Guoqing，YANG Lun，et al.Mining Subsidence［M］. Beijing:Mns University China Publishing Company，1991.

［4］李豐亮，郝兵元，楊雙鎖，等.大寧礦區(qū)開(kāi)采沉陷規(guī)律研究［J］.科技情報(bào)開(kāi)發(fā)與經(jīng)濟(jì)，2009，19(6):227-228.LI Fengliang，HAO Bingyun，YANG Shuangsuo，et al.Research on the mining subsidence laws of Daning mining area［J］. Sci-Tech Information Developmentand Economy，2009，19(6):227-228.

［5］張應(yīng)昌，孫亞峰，何堯.基于概率積分法的開(kāi)采沉陷預(yù)計(jì)分析［J］.山西建筑，2009，35(20):119-120.ZHANG Yingchang，SUN Yafeng，HE Rao.On mining subsidence anticipate analysis on the base of probabilityintegral method［J］.Shanxi Architecture，2009，35(20):119-120.

［6］李文秀，趙勝濤，梁旭黎，等.魯中礦區(qū)地下開(kāi)采對(duì)豎井井塔樓的影響分析［J］.巖石力學(xué)與工程學(xué)報(bào)，2006，25(1):74-78.LI Wenxiu，ZHAO Shengtao，LIANG Xuli，et al.Influence of underground mining on shaft towers in Luzhong Mining Areas［J］.Chinese Journal of Rock Mechanics and Engineering，2006，25(1):74-78.

［7］李新強(qiáng)，高延法，張慶松.開(kāi)采沉陷動(dòng)態(tài)數(shù)值仿真研究［J］. 巖石力學(xué)與工程學(xué)報(bào)，2004，23(1):86-90.LI Xinqiang，GAO Yanfa，ZHANG Qingsong.Study on dynamic numerical simulation of mining-induced subsidence［J］.Chinese Journal of Rock Mechanics and Engineering，2004，23(1):86-90.

［8］劉立民，劉漢龍，李建剛，等.基于 GIS的地下開(kāi)挖沉陷計(jì)算的有限層法及數(shù)據(jù)場(chǎng)表達(dá)［J］.巖石力學(xué)與工程學(xué)報(bào)，2005，24(12):2170-2175.LIU Limin，LIU Hanlong，LI Jiangang，et al.FLEM for ground subsidence calculation and data field expression model based on GIS［J］.Chinese Journal of Rock Mechanics and Engineering，2005，24(12):2170-2175.

［9］MARBLE D F.Some thoughts on the integration of spatial analysis and geographic information systems［J］.Journal of Geographical System，2000(2):31-35.

［10］柴華彬，鄒友峰，劉景艷.DTM在開(kāi)采沉陷可視化預(yù)計(jì)中的應(yīng)用［J］.遼寧工程技術(shù)大學(xué)學(xué)報(bào)，2004，23(2):171-174.CHAI Huabin，ZOU Youfeng ，LIU Jingyan.Application of DTM in visualization prediction of mining subsidence［J］.Journal of Liaoning Technical University，2004，23(2):171-174.

［11］李春雷，謝謨文，李曉璐.基于GIS和概率積分法的北洺河鐵礦開(kāi)采沉陷預(yù)測(cè)及應(yīng)用［J］.巖石力學(xué)與工程學(xué)報(bào)，2007，26(6):1243-1250.LI Chunlei，XIE Mowen，LI Xiaolu.Gis and probability integral based approach for subsidence prediction and application to beiminghe iron mine［J］.Chinese Journal of Rock Mechanics and Engineering，2007，26(6):1243-1250.