米麗倩，查劍鋒，劉丙方

（1.中國礦業(yè)大學(xué)，江蘇徐州 221116;2.江蘇省資源環(huán)境信息工程重點實驗室，江蘇徐州 221116; 3.國土環(huán)境與災(zāi)害監(jiān)測國家測繪局重點實驗室，江蘇徐州 221116;4.兗州煤業(yè)股份公司南屯煤礦，山東鄒城 273515）

概率積分法參數(shù)對預(yù)計地表下沉的影響度分析

米麗倩1，2，3，查劍鋒2，劉丙方4

（1.中國礦業(yè)大學(xué)，江蘇徐州 221116;2.江蘇省資源環(huán)境信息工程重點實驗室，江蘇徐州 221116; 3.國土環(huán)境與災(zāi)害監(jiān)測國家測繪局重點實驗室，江蘇徐州 221116;4.兗州煤業(yè)股份公司南屯煤礦，山東鄒城 273515）

地表沉陷預(yù)計的常用算法——概率積分法的預(yù)計參數(shù)及參數(shù)間交互作用對變形值不同的影響程度，造成參數(shù)選取、反演的不確定性。針對這一問題，采用五元二次正交旋轉(zhuǎn)組合設(shè)計，研究了概率積分法中預(yù)計參數(shù)對地表點最大下沉的影響關(guān)系，得出了單因素與兩兩因素交互作用對地表點最大下沉值的影響規(guī)律。研究結(jié)果表明:沉陷預(yù)計參數(shù)中的下沉系數(shù)、主要影響角正切、拐點偏移距對預(yù)計最大下沉具有顯著影響，影響程度排序為:下沉系數(shù)＞主要影響角正切＞拐點偏移距;且q與tanβ、tanβ與S/H的交互作用影響明顯。

概率積分法參數(shù);開采沉陷;最大下沉值;二次正交旋轉(zhuǎn)組合設(shè)計

下沉系數(shù)q、水平移動系數(shù)b、主要影響角正切tanβ、拐點偏移距S、影響傳播角θ，是概率積分法模型的主要預(yù)計參數(shù)。根據(jù)地表移動監(jiān)測站實測數(shù)據(jù)反演是獲取這些參數(shù)的主要方法，由于參數(shù)之間的相關(guān)性，使得反演結(jié)果過度依賴初值，且反演結(jié)果不惟一，造成使用過程中的困難。針對這一問題，采用試驗設(shè)計方法研究了概率積分法參數(shù)對預(yù)計結(jié)果的影響度，并分析了各因素之間相關(guān)性。

開采沉陷穩(wěn)定后，地表最大下沉值及其他移動變形值能從整體上反映地表變形的劇烈程度，而最大下沉值又是所有其他變量的決定參量之一［1］。因此，將最大下沉值作為影響度指標(biāo)來評價各參數(shù)對變形值的影響程度。本文采用二次正交旋轉(zhuǎn)組合設(shè)計建立了各因素與最大下沉值的回歸模型，并分析了各參數(shù)及其交互作用對最大下沉值的影響規(guī)律，對概率積分法參數(shù)選取和反演具有指導(dǎo)意義。

1 概率積分法模型基本原理

概率積分法是因其所用的移動和變形預(yù)計公式中含有概率積分 （或其導(dǎo)數(shù)）而得名，由于這種方法的基礎(chǔ)是隨機介質(zhì)理論，所以又叫隨機介質(zhì)理論法［2］。

1.1 最大下沉值計算基本方程

充分采動條件下最大下沉值為

式中，q為下沉系數(shù);m為煤層的開采厚度，m;α為煤層傾角，（°）。

實際上，大多數(shù)情況下地表沉陷預(yù)計區(qū)域為非充分采動區(qū)域，此時按照概率積分法理論，可將其視為兩個半無限開采的疊加，不失一般性。

傾斜煤層開采中半無限開采時，單元 ［坐標(biāo)為 （s，t）］的開采引起地表任意點［坐標(biāo)為 （x，y）］的下沉值為:

式中，r為主要影響半徑，r=H/tanβ;H為采深。

有限開采時煤層的開采等效于兩個半無限開采的幾何疊加，地表下沉值計算公式為:

式中，D1，D2分別為工作面傾向長和走向長;S上，S下，S左，S右為上、下、左、右邊界的拐點偏距。

1.2 參數(shù)影響的理論分析

由上述公式可知，W與q始終呈線性關(guān)系; tanβ與S/H參與運算較多，且變化范圍較大，故其變化可能會引起下沉值的較大改變，所以3個參數(shù)對下沉值的影響應(yīng)該較為顯著。而q，tanβ，S/ H的變化都會引起地表移動盆地邊界和面積的改變，故3個參數(shù)之間可能具有交互作用，但從繁雜的計算公式中難以準(zhǔn)確觀察出3個參數(shù)及其交互作用對下沉值的具體影響程度。計算公式中參數(shù)b不參與運算，可知b的變化對最大下沉值無影響。θ是傾向主斷面預(yù)計特有的參數(shù)，θ的變化使下沉盆地和傾斜方向上的拐點位置都產(chǎn)生一定的偏移，因此，當(dāng)θ的變化不改變移動盆地采動的充分性或非充分性時，其只改變傾向主斷面上的最大下沉值位置，對其數(shù)值無影響;但當(dāng)θ的變化改變移動盆地的采動程度時，將會使最大下沉值產(chǎn)生變化，由于θ的變化范圍較小，其對下沉值的影響甚微。

2 試驗方案

二次正交旋轉(zhuǎn)組合設(shè)計是正交試驗設(shè)計的一種，其既能分析處理單因素及因子間交互作用對指標(biāo)的影響，又能建立定量的數(shù)學(xué)模型，屬于更高級的試驗設(shè)計技術(shù)［3－4］。采用五因素五水平的試驗，根據(jù)二次正交旋轉(zhuǎn)組合設(shè)計的原則，僅需36次試驗就能全面分析各單因子及因子間的交互作用對最大下沉值的影響。

模擬試驗區(qū)的地質(zhì)采礦條件如圖1所示，若采厚m=3m，工作面傾向長D1=400m，走向長D2= 400m，下邊界采深H1=450m，上邊界采深H2= 350m，已知地面標(biāo)高為30m。

圖1 工作面采礦條件

試驗因素選取為q—X1，b—X2，tanβ—X3，S/H—X4，θ—X5，試驗指標(biāo)為最大下沉值W（mm）。采用二次正交旋轉(zhuǎn)組合設(shè)計進行試驗，因素水平編碼如表1所示。

表1 因素水平編碼

表2 試驗方案及結(jié)果

表2 試驗方案及結(jié)果（續(xù)）

3 試驗結(jié)果分析

3.1 數(shù)學(xué)模型的建立

利用DPS［5］軟件處理試驗數(shù)據(jù)，得到試驗因素對最大下沉值的影響分析表，如表3所示。

表3 試驗結(jié)果方差分析

由表3可得，因素X1，X3，X4，1X3，X3X4極顯著，各因素與最大下沉值之間的回歸模型為:

3.2 主因子效應(yīng)分析

由表3中F值可知，影響最大下沉值的各因素的主次順序依次為q，tanβ，S/H，θ，b。分別將回歸模型 （5）中的4個因素固定在零水平，另一因素固定在－2，－1，0，1，2五個水平，得到單因素對最大下沉值影響的變化曲線 （見圖2）。

圖2 主因子效應(yīng)

由圖2可知，X1與W的關(guān)系曲線為直線，W隨著X1的增加而增大;X3與W的關(guān)系曲線為拋物線，在［－2.0，1.3］水平區(qū)域內(nèi)，隨著X3水平的增加，W隨之增大，X3在1.3水平，W達到最大值;在［1.3，2.0］水平區(qū)域內(nèi)，隨著X3水平的增加，W隨之減小;X4與W的關(guān)系曲線為拋物線，隨著X4水平的增加，W隨之減小;X5與W的關(guān)系曲線為近似與X軸平行的直線，說明θ的變化對W基本無影響。X2與W的關(guān)系為與X軸平行的直線，說明b對W無影響。

3.3 兩因素交互作用分析

從回歸系數(shù)的絕對值可以看出，各因子間交互作用對最大下沉值的影響大小依次為:X3X4＞X1X3＞X1X4＞X2X5＞X1X5＞X1X2＞X4X5＞X3X5＞X2X3＞X2X4。從回歸系數(shù)顯著性檢驗可看出，q與tanβ，tanβ與S/H的交互作用對最大下沉值的影響明顯，且tanβ與S/H的交互作用大于q與tanβ的交互作用。

將X2，X4，X5固定在零水平，得到X1，X3兩因素交互作用對最大下沉值的影響模型為:

通過MATLAB的圖形繪制工具［6］，得到X1，X3兩因素交互作用對最大下沉值影響的響應(yīng)曲面圖，如圖3。

圖3 q和tanβ對最大下沉值的影響

由圖3可知:當(dāng)X3（tanβ）一定時，隨著X1（q）的增加，W逐漸增大;X1（q）一定時，隨著X3（tanβ）的增加，W有一個先增大后減小的過程。但q對結(jié)果的影響更為顯著，且當(dāng)X1，X3都為－2，即q=0.40，tanβ=1.2時，W取得最小值。

同理將X1，X2，X5固定在零水平，得到X3，X4兩因素交互作用對最大下沉值的影響模型為:

圖4為X3，X4兩因素交互作用對最大下沉值影響的響應(yīng)曲面圖。

由圖4可知，當(dāng)X4（S/H）為定值時，隨著X3（tanβ）的增加，W有一個先增大后減小的過程;當(dāng)tanβ為定值時，隨著S/H的增加，W逐漸減小。但tanβ對結(jié)果的影響更為顯著。當(dāng)X3在－2～－1，X4在1.5～2.0，即tanβ在1.2～1.65，S/H在0.17～0.20時，W值較小。

圖4 tanβ和S/H對最大下沉值的影響

4 結(jié)論

（1）通過五元二次正交旋轉(zhuǎn)組合設(shè)計試驗，建立了q，b，tanβ，S/H，θ五因素對最大下沉值的影響模型。

（2）根據(jù)單因子及兩因子交互作用分析表明，在試驗范圍內(nèi)，得出5個因素對最大下沉值影響的顯著性大小依次為q＞tanβ＞S/H＞θ＞b。q與tanβ，tanβ與S/H的交互作用對最大下沉值的影響明顯，且tanβ與S/H的交互作用大于q與tanβ的交互作用。此結(jié)論與對概率積分法模型的理論分析結(jié)果相吻合。

（3）參數(shù)選取和參數(shù)反演的優(yōu)化。應(yīng)用概率積分法預(yù)計最大下沉值時，根據(jù)參數(shù)的影響度依次選擇參數(shù)，抓住主要因素，忽略次要因素，提高預(yù)計的精度。由于q與tanβ，tanβ與S/H的交互作用對最大下沉值的影響明顯，應(yīng)用最大下沉值進行概率積分法參數(shù)反演時應(yīng)將tanβ作為首要反演值，以減少反演時參數(shù)之間的相關(guān)性帶來的反演結(jié)果的不確定性。

［1］鄒友峰，鄧喀中，馬偉民，等.礦山開采沉陷工程［M］.徐州:中國礦業(yè)大學(xué)出社，2003.

［2］何國清，等.礦山開采沉陷學(xué)［M］.徐州:中國礦業(yè)大學(xué)出版社，1991.

［3］洪 偉，吳承禎.試驗設(shè)計與分析—原理·操作·案例［M］.北京:中國林業(yè)出版社，2004.

［4］茆詩松，周紀(jì)薌，陳 穎，等.試驗設(shè)計［M］.北京:中國統(tǒng)計出版社，2004.

［5］唐起義，馮光明，等.DPS數(shù)據(jù)處理系統(tǒng):實驗設(shè)計、統(tǒng)計分析及數(shù)據(jù)挖掘［M］.北京:科學(xué)出版社，2007.

［6］周建興，豈興明，矯津益，等.MATLAB從入門到精通［M］.北京:人民郵電出社，2008.

Analysis of Probability Integral Parameters'Influence on Predicting Surface Subsidence

MILi-qian1，2，3，ZHA Jian-feng2，LIU Bing-fang4

（1.China University of Mining＆Technology，Xuzhou 221116，China;2.Jiangsu Provincial Key Laboratory of Resources Environment Information Engineering，Xuzhou 221116，China; 3.Land Environment and Disaster Monitoring Key Laboratory of State Survey Bureau of China，Xuzhou 221116，China; 4.Nantun Colliery，Yanzhou Coal Corporation，Zoucheng 273515，China）

Different influence degree of parameters of probability integralmethod for surface subsidence prediction and mutual action on deformation value would cause the uncertain problem of parameter selection and back-analysis.In order to solve the problem，this paper applied secondary orthogonal rotary combined design with 5 factors to researching influence of prediction parameters on surface maximum subsidence.Influence rule of single factor and random double factors on maximum subsidence was obtained.Results showed that subsidence ratio，tangent ofmain influence angle and inflexion offsetweremain influence factors ofmaximum subsidence.The influence degree was listed as:subsidence ratio＞tangent of main influence angle＞inflexion offset.Mutual influence of q and tanβ，tanβand S/H was evident.

mining subsidence;parameters of probability integralmethod;maximum subsidence value;secondary orthogonal rotary combined design

TD327

A

1006-6225（2011）04-0013-04

2011－03－31

國家自然科學(xué)基金重點項目 （50834004）;地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護國家重點實驗室開放基金資助項目（SKLGP2010K002）;國土環(huán)境與災(zāi)害監(jiān)測國家測繪局重點實驗室開放基金 （LEDM2009B01）

米麗倩 （1989－），女，山東聊城人，研究生，主要研究方向為開采沉陷預(yù)計，變形監(jiān)測數(shù)據(jù)處理。

［責(zé)任編輯:徐乃忠］