郭海英，王長(zhǎng)建

(1.山西陽泉固莊煤礦，山西 陽泉 045060;2.河南理工大學(xué)能源學(xué)院，河南 焦作 454000)

在大多情況下，礦山壓力顯現(xiàn)會(huì)對(duì)采礦工程造成不同程度的危害。為了使礦山壓力顯現(xiàn)對(duì)正常的采礦工作不產(chǎn)生影響，必須采取各種措施把礦山壓力顯現(xiàn)控制在一定范圍內(nèi)。以避免造成巨大的人員傷亡和財(cái)產(chǎn)損失。因此，事先對(duì)各個(gè)煤層處的礦山壓力進(jìn)行估測(cè)顯得越來越重要。

礦山壓力現(xiàn)場(chǎng)觀測(cè)研究法是在生產(chǎn)現(xiàn)場(chǎng)對(duì)礦山壓力大小進(jìn)行觀測(cè)和記錄，利用各種儀表、工具觀測(cè)采場(chǎng)圍巖變形、位移、頂?shù)装迤茐?、支架受載荷及壓縮等情況，然后從動(dòng)態(tài)分析中得出對(duì)采場(chǎng)礦山壓力顯現(xiàn)有顯著影響的巖層運(yùn)動(dòng)的條件、時(shí)間和方向以及采場(chǎng)周圍支承壓力的變化狀況，為解決現(xiàn)場(chǎng)礦山壓力控制問題提供依據(jù)。由于現(xiàn)場(chǎng)觀測(cè)獲得的資料反映多種因素的相互作用，用已經(jīng)測(cè)得淺部礦山壓力數(shù)據(jù)預(yù)測(cè)深處礦山壓力的大小，具有可靠方便的特點(diǎn)。故利用這種資料預(yù)測(cè)分析深部礦山壓力大小較為可靠，效果顯著，越來越為生產(chǎn)現(xiàn)場(chǎng)認(rèn)可采用。

為此，本文運(yùn)用常微分方程的數(shù)值解為方法，根據(jù)實(shí)驗(yàn)數(shù)據(jù)，研究得出煤炭埋深與其所受礦山壓力大小的變化規(guī)律，進(jìn)而對(duì)煤炭埋深處的礦山壓力值進(jìn)行預(yù)測(cè)，為煤炭開采時(shí)開挖巷道所需支護(hù)強(qiáng)度進(jìn)行預(yù)測(cè)并提供理論依據(jù)。

1 問題分析

某煤礦煤炭賦存深度平均約為500 m，現(xiàn)在需要初步預(yù)測(cè)煤層處所受的礦山壓力大小。根據(jù)煤層上方0～50 m測(cè)量得到一組數(shù)據(jù)，見表1，根據(jù)測(cè)得的數(shù)據(jù)得到此處礦山壓力大小，大體可用微分方程x2+x－y來表示，其中y是井下某煤層處的礦山壓力，x是對(duì)應(yīng)的煤層埋深，且假定地表平面的壓力為零。

2 常微分方程的數(shù)值解法

常微分方程(組)的解是一個(gè)函數(shù)或一族函數(shù)。所謂數(shù)值解法就是尋求解函數(shù)自變量的一系列離散點(diǎn)上的近似值。常微分方程的數(shù)值解法很多，常用的是歐拉公式?？紤]微分方程的初值問題:

把微分方程中的導(dǎo)數(shù)用折線的斜率代替，得到:

假設(shè)對(duì)于不同的n有xn+1－xn相同，記為h=xn+1－xn得到y(tǒng)n+1=yn+hf(xn，yn)基本的歐拉公式，其中，h稱為步長(zhǎng)。

3 算法的Matlab實(shí)現(xiàn)

3.1 實(shí)驗(yàn)數(shù)據(jù)

在本實(shí)驗(yàn)中，礦山工作人員已測(cè)得u在50 m內(nèi)的礦山壓力大小，見表1。

表1 50 m內(nèi)礦山壓力大小測(cè)量值

3.2 Matlab程序代碼

常微分方程數(shù)值解法Matlab程序:

4 計(jì)算結(jié)果及分析

4.1 計(jì)算結(jié)果

為了便于對(duì)實(shí)驗(yàn)數(shù)據(jù)和數(shù)值解進(jìn)行比較，在此把實(shí)驗(yàn)中得到的兩種數(shù)據(jù)放在一起，具體數(shù)據(jù)見表2。

表2 不同u值得到的兩種解析值和數(shù)值

4.2 Matlab圖形

歐拉公式計(jì)算的數(shù)值解與解析解示意圖見圖1。

圖1 歐拉公式計(jì)算的數(shù)值解與解析解示意圖

4.3 結(jié)果分析

從表2和圖1可知，用數(shù)值法得到的數(shù)據(jù)與用解析方法得到的數(shù)據(jù)基本是一致的，當(dāng)u=5時(shí)，誤差為3，當(dāng)u=50時(shí)，誤差為25。隨著u的增大，數(shù)值解與解析解的誤差也逐漸增大，但是變化量比較小。從圖1可以看出，隨著埋深的增大，礦山壓力的值越來越大，而且，開始時(shí)礦山壓力增加比較緩慢，深度越大礦山壓力增長(zhǎng)越快。由此得出的礦山壓力變化情況，進(jìn)而得出不同埋深處所受礦上壓力的大體范圍。在此可以預(yù)測(cè)u為500 m時(shí)，礦山壓力值為21 MPa。

5 結(jié)論

本文運(yùn)用常微分方程數(shù)值解法，分析了埋深變化對(duì)礦山壓力的影響，根據(jù)預(yù)測(cè)方程得出了不同埋深與其所受礦山壓力大小的關(guān)系，結(jié)果表明:

1)不同埋深的煤層所受的礦山壓力大小不同，埋深越深，所受到的礦山壓力越大。

2)在一定范圍內(nèi)，可以根據(jù)較淺深度的礦山壓力值預(yù)測(cè)較深處的礦山壓力大小，預(yù)測(cè)結(jié)果與真實(shí)值比較接近。

［1］錢鳴高，石平五，許家林.礦山壓力與巖層控制［M］.徐州:中國(guó)礦業(yè)大學(xué)出版社，2010:353－371.

［2］王高雄，周之銘.常微分方程［M］.北京:高等教育出版社，2006:19－21.

［3］孫 祥，徐流美，吳 清.matlab7.0基礎(chǔ)教程［M］.北京:清華大學(xué)出版社，2005:382－385.

［4］王 輝.大埋深采煤工作面礦壓特征的數(shù)值模擬［J］.煤炭科學(xué)技術(shù)，2010(6):73－75.