許 冬 王臨清 吳 侃

(1.中國礦業(yè)大學(xué)環(huán)境與測繪學(xué)院，江蘇 徐州 221116；2.環(huán)境保護部環(huán)境工程評估中心，北京 100012)

任意形狀工作面沉陷預(yù)計計算方法

許 冬1王臨清2吳 侃1

(1.中國礦業(yè)大學(xué)環(huán)境與測繪學(xué)院，江蘇 徐州 221116；2.環(huán)境保護部環(huán)境工程評估中心，北京 100012)

為彌補現(xiàn)有任意形狀工作面沉陷預(yù)計計算方法存在的不足，更加準確地預(yù)計煤炭開采引起的地表移動和變形，基于概率積分法理論，將工作面按照角點劃分為若干梯形，再利用變步長Simpson二重積分法在各梯形區(qū)域按預(yù)計公式進行積分計算，最后將各梯形單元產(chǎn)生的影響進行疊加處理，實現(xiàn)了任意形狀工作面的沉陷預(yù)計計算。另外，通過研究拐點偏移距對工作面位置的影響，探討了任意形狀工作面實際計算邊界角點坐標的解算方法。在此基礎(chǔ)上設(shè)計、開發(fā)了開采沉陷預(yù)計程序，實現(xiàn)了地表移動變形預(yù)計、預(yù)計結(jié)果分析及可視化輸出。以某礦區(qū)實際開采工作面為例，對計算方法及預(yù)計程序進行了驗證。結(jié)果表明：程序設(shè)計合理，預(yù)計結(jié)果準確可靠，符合實際情況，能夠滿足大范圍沉陷預(yù)計的需要，可為“三下”采煤、土地復(fù)墾規(guī)劃提供必要的技術(shù)支持。

任意形狀工作面 開采沉陷 概率積分法 移動變形 變步長Simpson方法 拐點偏移距

煤炭資源的開采在我國高速發(fā)展的國民經(jīng)濟中發(fā)揮了重要的作用，同時也給礦區(qū)生態(tài)環(huán)境造成了比較嚴重的影響[1-2]。礦區(qū)開采沉陷預(yù)計是礦山開采沉陷學(xué)科的核心內(nèi)容之一[3]，目前針對任意形狀工作面沉陷預(yù)計計算的問題，國內(nèi)學(xué)者從不同角度出發(fā)，提出了各具特色的解決辦法[4-12]。文獻[4-6]利用數(shù)學(xué)手段將面積分轉(zhuǎn)化為線積分進行非矩形工作面的沉陷預(yù)計，雖然實現(xiàn)了任意煤層傾角任意形狀工作面的地表移動與變形計算，但存在數(shù)據(jù)準備繁瑣，不能快速實現(xiàn)多工作面預(yù)計計算等問題。文獻[7-8]提出將不規(guī)則工作面按照角點劃分成若干三角形,用多重積分的高斯方法求解地表移動和變形的方法，雖然預(yù)計結(jié)果精度較高，但存在積分上下限確定困難，計算速度慢等問題。文獻[9-10]提出采用將非矩形工作面劃分為若干個近似小矩形工作面，最后將預(yù)計結(jié)果進行疊加處理的方法進行預(yù)計計算。雖然方法操作簡單、便于理解，但存在預(yù)計結(jié)果帶有近似性、不能處理凹多邊形工作面、下沉等值線在開采邊界處呈現(xiàn)“鋸齒狀”等問題。

筆者借鑒前人研究成果，應(yīng)用概率積分法的基本原理，提出了任意形狀工作面沉陷預(yù)計計算的具體算法，并開發(fā)了一個能夠?qū)θ我庑螤疃喙ぷ髅孢M行沉陷預(yù)計與分析的程序，以期為評價“三下”采煤的可行性提供依據(jù)。

1 理論模型及關(guān)鍵算法

1.1 理論模型

概率積分法因其所用的移動和變形預(yù)計公式中含有概率積分而得名[2]，是以隨機介質(zhì)理論為基礎(chǔ)，以正態(tài)分布函數(shù)為影響函數(shù)，用積分式準確合理描述下沉盆地的方法?；诟怕史e分法的全盆地移動與變形預(yù)計計算公式如下[11]。

下沉：

(1)

水平移動：

(2)

(3)

式中，W(x,y)為地表任意點(坐標為x、y)的下沉值，mm；Ux(x,y)為地表任意點走向方向的水平移動值，mm；Uy(x,y)為地表任意點傾斜方向的水平移動值，mm；Wcm、Ucm為地表充分采動的最大下沉值和最大水平移動值，mm；r為走向主斷面上采空區(qū)邊界，m；x,y為計算點相對坐標(考慮拐點偏移距)；D為開采煤層區(qū)域；η，ξ分別為x,y方向的變形值，mm。

1.2 梯形劃分與積分計算

由式(1)～式(3)可以看出，按照全盆地移動與變形公式進行預(yù)計計算主要是在區(qū)域D上進行概率積分運算。對于這類復(fù)雜積分問題的求解，只有給出被積函數(shù)和積分上、下限才可進行數(shù)值積分計算。積分限與開采區(qū)域D有關(guān)，當開采區(qū)域是非矩形工作面時很難直接確定其積分限[8]。本研究嘗試采用將工作面按照角點劃分為若干梯形的方法來解決概率積分計算困難的問題，具體方法如下。

設(shè)開采區(qū)域D為任意凸多邊形工作面(如圖1(a)所示)，從A～G各頂點向X軸做垂線，按照A—B—C—E—F—G—A的角點順序(逆時針方向)將開采區(qū)域劃分為AA′B′B、BB′C′C、CC′EE′、EE′FF′、FF′GG′、GG′AA′共6個梯形單元，對于每個梯形單元可以唯一確定其積分的上下限，如AA′B′B單元的積分限如式(4)所示：

(4)

圖1 積分區(qū)域劃分Fig.1 Integral area division

每個梯形單元的積分上、下限確定后，分別采用高效的變步長Simpson二重積分法進行數(shù)值積分計算，計算完成后將各梯形單元產(chǎn)生的影響進行疊加處理來計算整個采空區(qū)造成的地表移動變形。由于設(shè)定工作面角點按照逆時針順序(A—B—C—E—F—G—A)排列，根據(jù)正負判斷法則(如果x1x2，則S11′2′2前符號為正)，開采區(qū)域D的移動變形值可用下式表示：

一處狹窄的巷子入口，巷子里面黑黢黢的，黑的像是要把我吞沒。入口處一盞昏暗的路燈，仿佛隨時要被夜風(fēng)吹滅。踩踏板的腳頓了一下。但我終究還是向這吃人般的黑暗騎去。

(5)

當工作面為凹多邊形時，仍可按凸多邊形的處理方法進行處理，凹多邊形工作面的梯形劃分如圖1(b)所示。因此，無論工作面是凸多邊形還是凹多邊形，均可使用梯形劃分方法，確定積分限進行概率積分計算，求解地表移動變形，使預(yù)計計算不再受開采工作面布置形式的影響。

1.3 工作面計算邊界求取

拐點偏移距是開采沉陷預(yù)計過程中的一個重要參數(shù)，它的大小直接影響到地表下沉盆地的形狀、范圍以及移動變形值的大小[13-15]。考慮拐點偏移距后，工作面的開采邊界將會發(fā)生變化，求取任意形狀工作面實際開采邊界的主要步驟是：①根據(jù)工作面角點坐標計算工作面的重心；②以重心為交點分別做X、Y軸的平行線，將工作面劃分為4個子區(qū)域(Ⅰ～Ⅳ)；③先判斷角點屬于哪一個子區(qū)域，再根據(jù)其所屬區(qū)域的不同選擇相應(yīng)的計算方法，如圖2所示。圖2中S1～S4分別為走向左、右邊界拐點偏移距和上、下山邊界拐點偏移距。若在采空區(qū)一側(cè)，取正值，在煤柱一側(cè)，取負值[2]。

圖2 分區(qū)域計算方法Fig.2 Calculation method of partition

由于煤層存在一定的傾角，使計算邊界與考慮拐點偏移距后的實際開采邊界不完全一致，兩者在煤層傾斜方向存在一個偏移量，該偏移量d可按式(6)計算：

(6)

式中，h為下邊界采深，m;S4為下山邊界拐點偏移距，m;α為煤層傾角，(°);θ為開采影響傳播角,(°)。

當偏移量d為正值時，從實際下山邊界向下山方

2 程序?qū)崿F(xiàn)流程

為了能夠更加方便、準確地進行開采沉陷預(yù)計計算，應(yīng)用軟件工程思想并結(jié)合工程實際需求，以“數(shù)據(jù)”為核心，按照數(shù)據(jù)準備、處理和分析的原則對任意形狀工作面沉陷預(yù)計程序進行設(shè)計與開發(fā)，主要實現(xiàn)地表移動變形預(yù)計、預(yù)計結(jié)果分析及可視化輸出等功能。程序?qū)崿F(xiàn)的基本流程是：①構(gòu)建預(yù)計文件。從采掘工程平面圖上獲得工作面尺寸信息、開采參數(shù)、沉陷預(yù)計參數(shù)，并圈定預(yù)計范圍、設(shè)定格網(wǎng)間距。②沉陷預(yù)計計算。根據(jù)預(yù)計文件中參數(shù)進行坐標轉(zhuǎn)換、求取工作面計算邊界，將計算邊界按角點劃分為若干梯形，分別對每一個梯形進行二重積分將積分結(jié)果進行疊加處理，得到下沉、水平移動等預(yù)計結(jié)果，進而繪制移動變形等值線、云圖等圖件；其中，對于多個工作面的沉陷預(yù)計采用將每個工作面產(chǎn)生的影響進行疊加的方法得到最終的地表移動變形值[16-17]。③沉陷結(jié)果分析。根據(jù)得到的預(yù)計結(jié)果可計算煤炭開采導(dǎo)致的沉陷體積和沉陷面積、制作開采影響后的礦區(qū)地形圖、評價區(qū)域內(nèi)建筑物的損壞等級、分析沉陷范圍內(nèi)的積水情況[18-19]。程序?qū)崿F(xiàn)的流程如圖3所示。

圖3 程序?qū)崿F(xiàn)流程Fig.3 The program implementation process

3 實例應(yīng)用

3.1 實驗區(qū)概況

為測試預(yù)計結(jié)果的準確性，選取某礦區(qū)實際開采工作面進行任意形狀、多工作面的沉陷預(yù)計計算，工作面如圖4所示。該測試算例既包含凹多邊形工作面(如工作面1、2)又有凸多邊形工作面(如工作面3)，能夠滿足程序測試的要求。

圖4 工作面布置Fig.4 Layout of working face

結(jié)合礦區(qū)實測資料，通過綜合分析，確定了沉陷預(yù)計時所使用的預(yù)計參數(shù)，見表1。

3.2 沉陷預(yù)計結(jié)果

按照構(gòu)建預(yù)計文件、沉陷預(yù)計計算、預(yù)計結(jié)果分析的流程對實驗區(qū)的3個工作面進行處理，生成下沉、傾斜、曲率、水平移動、水平變形等值線及下沉云圖。由于篇幅有限，只給出部分結(jié)果圖件，其中圖5為下沉等值線圖，圖6為下沉云圖，預(yù)計數(shù)值結(jié)果見表2。實例應(yīng)用表明，通過采用基于梯形劃分的任意形狀工作面沉陷預(yù)計計算方法可使預(yù)計結(jié)果與實際情況更加符合。

4 結(jié) 論

(1)本研究采用將工作面按照角點劃分為若干梯形單元，用預(yù)計計算公式計算每一個梯形單元造成的地表移動變形，最后再將各小單元的影響進行疊加來計算整個采空區(qū)造成的地表移動變形的方法實現(xiàn)了任意形狀工作面的沉陷預(yù)計計算，解決了概率積分求解困難的問題，使預(yù)計計算不再受開采工作面布置形式的影響，為任意形狀工作面開采沉陷預(yù)計計算提供了新的解決思路。

(2)通過研究拐點偏移距對工作面實際位置的影響，探討了求取任意形狀工作面實際計算邊界角點坐標的方法，使預(yù)計結(jié)果更符合地表移動與變形的實際情況。實現(xiàn)了多工作面沉陷預(yù)計計算，使預(yù)計結(jié)果能夠滿足礦區(qū)大范圍開采沉陷預(yù)計的需要。

表1 開采沉陷預(yù)計參數(shù)Table 1 Prediction parameters of mining subsidence

圖5 下沉等值線 (單位：m)Fig.5 Subsidence equivalent curve

圖6 下沉云圖(單位：mm)Fig.6 Subsidence nephogram表2 部分預(yù)計數(shù)值結(jié)果

Table 2 Part of prediction values mm

(3)通過研究相關(guān)沉陷預(yù)計理論、算法，設(shè)計、開發(fā)了沉陷預(yù)計程序，實現(xiàn)了地表移動變形預(yù)計、預(yù)計結(jié)果分析及可視化輸出。通過實例應(yīng)用表明：沉陷預(yù)計結(jié)果準確、可靠，程序具有良好的實用性與可擴展性，能夠簡單方便地進行開采沉陷預(yù)計，為“三下”采煤方案優(yōu)化設(shè)計、開展礦區(qū)土地復(fù)墾、保護礦區(qū)生態(tài)環(huán)境提供科學(xué)依據(jù)。

[1] 呂偉才,楊世清.開采沉陷對礦區(qū)生態(tài)環(huán)境的影響評價[J].江蘇環(huán)境科技,2006,19(S1):57-59. Lu Weicai，Yang Shiqing.The effect evaluation on the mining area ecological environment by the mining subsidence[J].Jiangsu Environmental Science and Techonlogy,2006,19(S1):57-59.

[2] 何國清,楊 倫,凌賡娣，等.礦山開采沉陷學(xué)[M].徐州:中國礦業(yè)大學(xué)出版社,1994. He Guoqing，Yang Lun，Ling Gengdi，et al.Mining Subsidence Theory[M].Xuzhou：China University of Mining and Technology Press,1994.

[3] 吳 侃,黃珍珍,韓奎峰.礦區(qū)土地破壞預(yù)測[J].山東煤炭科技,2006(3):48-49. Wu Kan，Huang Zhenzhen，Han Kuifeng.Prediction of the destruction on mining land[J].Shandong Coal Science and Technology,2006(3):48-49.

[4] 張華興.地表移動計算程序的可視化[J].煤炭科學(xué)技術(shù),1998(11):29-30. Zhang Huaxing.The visualization of earth′s surface mobile computing program[J].Coal Science and Technology，1998(11):29-30.

[5] 張華興.地表移動與變形的等價轉(zhuǎn)換線積分計算[J].煤礦開采,1991(1):14-19. Zhang Huaxing.The movement and deformation of ground equivalent conversion line integral calculation[J].Coal Mining Technology,1991(1):14-19.

[6] 張華興,焦傳武.地表移動精確預(yù)計法——線積分運算[J].礦山測量,1988(12):21-29. Zhang Huaxing，Jiao Chuanwu.Precise method of surface movement expected：line integral operation[J].Mine Surveying,1988(12):21-29.

[7] 康建榮.采動覆巖動態(tài)移動破壞規(guī)律及開采沉陷預(yù)計系統(tǒng)(MSPS)研究[D].北京:中國礦業(yè)大學(xué),1999. Kang Jianrong.The Rule of Overburden Dynamic Movement and Damage and Mining Subsidence Prediction System[D].Beijing：China University of Mining and Technology,1999.

[8] 康建榮,王金莊,溫澤民.任意形多工作面多線段開采沉陷預(yù)計系統(tǒng)(MSPS)[J].礦山測量,2000,1(3):24-27. Kang Jianrong，Wang Jinzhuang，Wen Zemin.Arbitrarily shaped multi-face multi-line of mining subsidence prediction system (MSPS)[J].Mine Surveying,2000,1(3):24-27.

[9] 吳 侃,周 鳴.開采沉陷預(yù)測預(yù)報系統(tǒng)[M].徐州:中國礦業(yè)大學(xué)出版社,1999. Wu Kan，Zhou Ming.Mining Subsidence Prediction System[M].Xuzhou：China University of Mining and Technology Press,1999.

[10] 吳 侃,葛家新,王玲丁,等.開采沉陷預(yù)計一體化方法[M].徐州:中國礦業(yè)大學(xué)出版社,1998. Wu Kan，Ge Jiaxin，Wang Lingding，et al.Integration Method of Mining Subsidence Prediction[M].Xuzhou：China University of Mining and Technology Press,1998.

[11] 國家煤炭工業(yè)局.建筑物、水體、鐵路及主要井巷煤柱留設(shè)與壓煤開采規(guī)程[S].北京:煤炭工業(yè)出版社,2004. State Bureau of Coal Industry.Mining Regulations of Buildings，Water，Railway and Main Roadway Pillar Design and Pressure Coal[S].Beijing：Coal Industry Press,2004.

[12] 顧 偉,譚志祥,鄧喀中.基于Active X Automation 的三維可視開采沉陷預(yù)計[J].金屬礦山,2013(3):119-122. Gu Wei，Tan Zhixiang，Deng Kazhong.3D visualization prediction of mining subsidence system based on the Active X Automation[J].Metal Mine,2013(3):119-122.

[13] 陳 巖,奚硯濤.基于Surfer的金屬礦山開采沉陷預(yù)測和三維可視化[J].金屬礦山,2012(5):149-151. Chen Yan，Xi Yantao.Metal mine mining subsidence prediction and 3D visualization based on surfer[J].Metal Mine,2012(5):149-151.

[14] 周萬茂,張華興,何瑞華.任意形狀工作面拐點移動距求取方法[J].煤礦開采,2000(12):13-16. Zhou Wanmao，Zhang Huaxing，He Ruihua.A method for determining the displacement distance of the inflection point of the mining faces with arbitrary shapes[J].Coal Mining Technology,2000(12):13-16.

[15] 張合兵,王世東,李春意.煤礦區(qū)土地復(fù)墾中土地損毀預(yù)測方法研究[J].金屬礦山,2013(9):142-147. Zhang Hebing，Wang Shidong，Li Chunyi.Research on land destruction prediction methods of land reclamation in coal mining area[J].Metal Mine,2013(9):142-147.

[16] 朱劉娟,陳俊杰,鄒友峰.任意形狀工作面開采地表移動變形預(yù)計的算法實現(xiàn)[J].遼寧工程技術(shù)大學(xué)學(xué)報:自然科學(xué)版,2005,24(3):337-339. Zhu Liujuan，Chen Junjie，Zou Youfeng.Algorithm of ground movement and deformation prediction in mining of random shape face[J].Journal of Liaoning Technical University：Natural Science Edition,2005，24(3):337-339.

[17] 李 明,李琰慶,王紅梅.任意形狀工作面開采沉陷預(yù)計系統(tǒng)開發(fā)[J].礦業(yè)安全與環(huán)保,2008,35(5):18-21. Li Ming，Li Yanqing，Wang Hongmei.Development of predicting system for mining subsidence of arbitrarily shaped face[J].Mining Safety & Environmental Protection,2008,35(5):18-21.

[18] 劉 虎,吳 侃,陳冉麗.基于AutoCAD的開采沉陷預(yù)計分析軟件研究[J].礦業(yè)安全與環(huán)保,2009,36(2):22-23. Liu Hu，Wu Kan，Chen Ranli.Research on mining subsidence prediction and analysis system based on AutoCAD[J].Mining Safety & Environmental Protection,2009,36(2):22-23.

[19] 石秀偉,李 晶,趙換新,等.基于ArcGIS的煤礦開采沉陷預(yù)測與可視化分析[J].金屬礦山,2012(9):103-106. Shi Xiuwei，Li Jing，Zhao Huanxin，et al.The prediction and visualization of coal mining subsidence based on ArcGIS[J].Metal Mine，2012(9):103-106.

(責(zé)任編輯 徐志宏)

Mining Subsidence Prediction Calculation Methods of Random Shape Working Face

Xu Dong1Wang Linqing2Wu Kan1

(1.ChinaUniversityofMiningandTechnology,SchoolofEnvironmentalScienceandSpatialInformaticsXuzhou221116，China；2.AppraisalCenterforEnvironment&EngineeringMinistryofEnvironmentalProtection，Beijing100012，China)

In order to make up for the shortage in mining subsidence prediction calculation methods of existing random shape working face，and more accurately predict the value of surface movement and deformation caused by coal mining，the working face is divided into several trapezoids according to corners based on the theory of probability integral method.The variable step-size Simpson double integral method is adopted to conduct the integration in every trapezoidal region according to the prediction formula，and deal with the effects which were produced by every trapezoidal unit to realize the subsidence prediction calculation of random shape working face through superposition method.In addition，through the study of inflection point displacement effect on the working face position，the solution method for actually calculating angular coordinate at boundary of the random shape working face is discussed.On this basis，the mining subsidence prediction program is designed and developed，realizing prediction of surface movement and distortion，analysis of prediction results and visual output.The actual mine working faces in a mining area was taken as a case to verify the calculation method and prediction program.The results showed that the program design is reasonable，and the prediction result is accurate and reliable，and accords with the fact.It can meet the needs for a wide range of subsidence prediction and provide the necessary technical support for “under three” coal mining and land reclamation planning.

Random shape working face，Mining subsidence，Probability integral method，Movement and deformation，Variable step-size Simpson-method，Inflection point displacement

2014-01-02

“十二五”國家科技支撐計劃項目(編號：2012BAC04B03)。

許 冬(1990—)，男，碩士研究生。

TD325+.2

A

1001-1250(2014)-05-055-05