張官進(jìn)，閆 威，江克貴

( 1.安徽科技學(xué)院 資源與環(huán)境學(xué)院，安徽 滁州 233100；2.安徽理工大學(xué) 測(cè)繪學(xué)院，安徽 淮南 232001 )

開采沉陷預(yù)測(cè)、防治理論及相關(guān)技術(shù)是我國(guó)礦山領(lǐng)域研究的熱點(diǎn)，求解精確、可靠的開采沉陷預(yù)計(jì)參數(shù)也是預(yù)計(jì)模型應(yīng)用的難點(diǎn)。傳統(tǒng)的開采沉陷預(yù)計(jì)參數(shù)求解方法主要有線性近似法[1-2]、特征點(diǎn)法和正交試驗(yàn)法[3]。其中，線性近似法對(duì)地表監(jiān)測(cè)站的布設(shè)形式以及求參模型初值準(zhǔn)確性要求較高，在工程應(yīng)用中具有較大的困難；對(duì)于特征點(diǎn)法，由于移動(dòng)與變形曲線復(fù)雜，曲線特征點(diǎn)難以精確確定，進(jìn)而導(dǎo)致求參誤差較大；正交試驗(yàn)法存在試驗(yàn)次數(shù)多、求參速度慢、計(jì)算機(jī)不易實(shí)施等缺點(diǎn)。針對(duì)常規(guī)求參方法存在的問題，大量學(xué)者把智能優(yōu)化算法引入到開采沉陷預(yù)計(jì)參數(shù)求取中來(lái)，葛家新[4]等將模矢法( 步長(zhǎng)加速法 )應(yīng)用到求取開采沉陷預(yù)計(jì)參數(shù)中，解決了任意形狀工作面和動(dòng)態(tài)實(shí)測(cè)數(shù)據(jù)求參問題；査劍鋒[5]等提出采用遺傳算法反演概率積分參數(shù)，通過工程應(yīng)用證實(shí)了遺傳算法求參的準(zhǔn)確性與可靠性；后來(lái)，眾多學(xué)者又引入粒子群算法[6-7]、改進(jìn)人工蜂群算法[8]、果蠅算法[9]、BFGS算法[10]等方法。相較于常規(guī)求參方法，基于智能優(yōu)化算法的求解參數(shù)精度得到了大大提高。但是，由于該類算法尋優(yōu)過程大同小異，自身求參能力已經(jīng)達(dá)到瓶頸，同時(shí)受概率積分模型擬合度的限制，求參精度很難得到進(jìn)一步提高。

綜上，為了進(jìn)一步提高求參的準(zhǔn)確性，本文針對(duì)開采沉陷預(yù)計(jì)模型重點(diǎn)開展研究。文獻(xiàn)[11-13]表明，組合預(yù)計(jì)模型綜合利用各種預(yù)計(jì)模型的信息，以適當(dāng)?shù)募訖?quán)形式得出組合預(yù)計(jì)模型，相對(duì)于單項(xiàng)預(yù)計(jì)模型，組合預(yù)計(jì)模型能夠有效提高擬合精度和預(yù)測(cè)能力。另外，一些學(xué)者[14]提出了基于Boltzmann函數(shù)的開采沉陷預(yù)計(jì)模型，通過分析曲線形態(tài)和參數(shù)特點(diǎn)發(fā)現(xiàn)，相對(duì)于概率積分模型，Boltzmann函數(shù)預(yù)計(jì)模型預(yù)計(jì)曲線在邊緣部分收斂緩慢，且預(yù)計(jì)參數(shù)與概率積分參數(shù)可以相互轉(zhuǎn)化，2種預(yù)計(jì)模型能夠共用1套參數(shù)系統(tǒng)。因此，基于變權(quán)組合預(yù)計(jì)模型原理，融合概率積分模型和Boltzmann函數(shù)預(yù)計(jì)模型，構(gòu)建了1種變權(quán)PB組合預(yù)計(jì)模型，鑒于該模型屬于高度非線性模型，耦合煙花算法，提出了1種基于變權(quán)PB組合預(yù)計(jì)模型的開采沉陷預(yù)計(jì)參數(shù)反演方法。

1 單項(xiàng)預(yù)計(jì)模型原理

1.1 概率積分模型

概率積分模型由我國(guó)學(xué)者提出，目前已成為我國(guó)較為成熟、應(yīng)用最為廣泛的開采沉陷預(yù)計(jì)方法之一[15]。根據(jù)概率積分模型開采沉陷預(yù)計(jì)原理，半無(wú)限和有限開采引起地表下沉的曲線公式為

式中，

其中，W0為最大下沉值；m為工作面采厚；q為下沉系數(shù)；α為工作面傾角；為概率積分函數(shù)；r為主要影響半徑，r＝H/tan β；H為工作面采深；tan β為主要影響角正切；L為工作面計(jì)算開采長(zhǎng)度。

走向計(jì)算開采長(zhǎng)度為

傾向計(jì)算開采長(zhǎng)度為

式中，D3，D1分別為走向和傾向開采長(zhǎng)度；S1，S2，S3，S4分別為傾向下山、傾向上山、開切眼和終采線處拐點(diǎn)偏移距；θ為開采影響傳播角。

水平移動(dòng)預(yù)計(jì)公式是由單元盆地下沉影響函數(shù)推導(dǎo)而來(lái)，具體推導(dǎo)過程可參考文獻(xiàn)[16]，半無(wú)限和有限開采引起地表水平移動(dòng)的曲線公式為

式中，b為水平移動(dòng)系數(shù)。

根據(jù)有限開采預(yù)計(jì)公式，整個(gè)工作面開采引起地表任意點(diǎn)( x，y )沿φ方向( x軸的正向逆時(shí)針旋轉(zhuǎn)到指定方向的角值 )的下沉和水平移動(dòng)為

1.2 Boltzmann函數(shù)預(yù)計(jì)模型

近年來(lái)，一些學(xué)者基于Boltzmann函數(shù)提出了一種新的開采沉陷預(yù)計(jì)模型。對(duì)非充分開采，充分開采以及超充分開采的實(shí)測(cè)結(jié)果分析表明，Boltzmann函數(shù)可以較好地預(yù)計(jì)不同開采條件下的地表開采沉陷形變。根據(jù)文獻(xiàn)[14]，半無(wú)限開采地表形變盆地下沉預(yù)計(jì)公式為

式中，R為新的主要影響半徑；S為拐點(diǎn)偏移距。

對(duì)式( 4 )求微分，可得Boltzmann函數(shù)預(yù)計(jì)模型的單元影響函數(shù)為

仿照概率積分法推導(dǎo)經(jīng)驗(yàn)和彈性力學(xué)知識(shí)[17]，單元水平移動(dòng)函數(shù)和半無(wú)限開采地表水平移動(dòng)預(yù)計(jì)公式為

式中，b'為新的水平移動(dòng)系數(shù)。

有限開采情況下，基于Boltzmann函數(shù)的地表下沉和水平移動(dòng)預(yù)計(jì)公式Wb( x )，Wb( y )，Ub( x )和Ub( y )可按照1.1節(jié)公式推導(dǎo)過程推理得到。因此，工作面開采引起地表任意點(diǎn)( x，y )沿φ方向的下沉和水平移動(dòng)預(yù)計(jì)公式為

1.3 模型參數(shù)及形變擬合分析

綜上所述，結(jié)合文獻(xiàn)[14]的研究，開采沉陷預(yù)計(jì)參數(shù)體系為P＝[q，tan β，b，θ，S1，S2，S3，S4]。概率積分模型和Boltzmann函數(shù)預(yù)計(jì)模型在tan β和b上參數(shù)意義相同，數(shù)值不同，但參數(shù)數(shù)值可以通過經(jīng)驗(yàn)公式轉(zhuǎn)化；新的主要影響半徑R≈r/4.13，由于tan β＝H/r，則tan β'＝4.13tan β，又 由 于b ＝B/r，則b'＝4.13b。另外，2種預(yù)計(jì)模型在其余預(yù)計(jì)參數(shù)上意義相同，數(shù)值相等。因此，通過一定的數(shù)值轉(zhuǎn)化，2種預(yù)計(jì)模型可共用1套預(yù)計(jì)參數(shù)系統(tǒng)P，這給2種模型變權(quán)組合進(jìn)而反演預(yù)計(jì)參數(shù)提供了理論前提。

研究表明，概率積分模型在大變形區(qū)域預(yù)計(jì)效果較優(yōu)，在變形盆地邊緣區(qū)域收斂過快，與實(shí)際變形有較大差別[18]；而Boltzmann函數(shù)預(yù)計(jì)模型在變形盆地邊緣區(qū)域收斂緩慢。理論上，融合2種單項(xiàng)預(yù)計(jì)模型的組合預(yù)計(jì)模型能夠改進(jìn)形變擬合程度，進(jìn)而提高求解參數(shù)的精度。在1套預(yù)計(jì)參數(shù)體系下，2種預(yù)計(jì)模型預(yù)計(jì)下沉( 半無(wú)限開采 )比較如圖1所示。

2 變權(quán)PB組合預(yù)計(jì)模型及反演方法構(gòu)建

2.1 變權(quán)組合預(yù)計(jì)模型

概率積分模型和Boltzmann函數(shù)預(yù)計(jì)模型組合優(yōu)勢(shì)在于，2種模型可共用1套概率積分參數(shù)；同時(shí)Boltzmann函數(shù)模型預(yù)計(jì)曲線在盆地邊緣區(qū)域收斂緩慢，以及概率積分模型在大變形區(qū)域的擬合能力，理論上可使得組合模型在全盆地?cái)M合能力顯著提升。

基于變權(quán)組合理論，融合概率積分模型和Boltzmann函數(shù)預(yù)計(jì)模型，構(gòu)造的變權(quán)PB組合預(yù)計(jì)模型的下沉和水平移動(dòng)預(yù)計(jì)公式分別為

圖1 半無(wú)限開采下沉曲線比較Fig.1 Comparison of subsidence curves in semi-infinite mining

式中，Wp( x，y )，Wb( x，y )分別為概率積分模型和Boltzmann函數(shù)預(yù)計(jì)模型預(yù)計(jì)下沉；Up( x，y )，Ub( x，y )分別為概率積分模型和Boltzmann函數(shù)預(yù)計(jì)模型預(yù)計(jì)水平移動(dòng)；ρ為權(quán)重。

綜上，概率積分參數(shù)體系為P，假設(shè)地質(zhì)采礦概況參數(shù)為G＝[m，α，θ，H，D3，D1]，變權(quán)PB組合預(yù)計(jì)模型函數(shù)可以概括為

2.2 參數(shù)反演方法構(gòu)建

在工作面地質(zhì)采礦概況和地表形變量已知的情況下，如何利用構(gòu)建的變權(quán)PB組合預(yù)計(jì)模型( 高度非線性模型 )精確地反演開采沉陷預(yù)計(jì)參數(shù)，本文引入了煙花算法。

近年來(lái)，TAN和ZHU[19]根據(jù)煙花爆炸產(chǎn)生火花這一自然現(xiàn)象提出了煙花算法。煙花算法是1種具有局部和全局搜索能力的自調(diào)節(jié)機(jī)制群體智能優(yōu)化算法[20]。煙花種群中各個(gè)煙花根據(jù)其相對(duì)于其他煙花的適應(yīng)度進(jìn)行資源分配和信息交互，使得種群能夠在全局和局部搜索能力之間達(dá)到平衡，其強(qiáng)大的尋優(yōu)搜索能力使得煙花算法成為1種新型智能優(yōu)化算法。煙花算法主要由爆炸算子、變異算子、映射規(guī)則和選擇策略組成，它們決定了煙花算法尋優(yōu)性能的優(yōu)劣。具體的煙花算法尋優(yōu)過程和相關(guān)公式可參考文獻(xiàn)[21-22]。

綜合變權(quán)PB組合預(yù)計(jì)模型特點(diǎn)、地表三維形變以及煙花算法原理，提出了1種基于變權(quán)PB組合預(yù)計(jì)模型的開采沉陷預(yù)計(jì)參數(shù)反演方法研究，其具體步驟如下：

( 1 ) 數(shù)據(jù)準(zhǔn)備。獲取目標(biāo)礦區(qū)工作面地表實(shí)測(cè)形變以及工作面地質(zhì)采礦概況參數(shù)；選定適當(dāng)?shù)臒熁ㄋ惴▍?shù)；以[P，ρ]＝[q，tan β，b，θ，S1，S2，S3，S4，ρ]作為種群的個(gè)體( 煙花或火花 )，同時(shí)確定種群活動(dòng)的可行域范圍，進(jìn)而生成煙花初始種群。

( 2 ) 適應(yīng)度評(píng)價(jià)。根據(jù)煙花初始種群和地質(zhì)概況參數(shù)，基于組合預(yù)計(jì)模型，煙花種群中第i組參數(shù)預(yù)計(jì)的下沉和水平移動(dòng)分別為Wi和Ui；假設(shè)實(shí)測(cè)下沉和水平移動(dòng)分別為Wm和Um，以預(yù)計(jì)值和實(shí)測(cè)值之差的平方和為準(zhǔn)則進(jìn)行適應(yīng)度評(píng)價(jià)，適應(yīng)度函數(shù)F的構(gòu)造如式( 10 )所示。判斷適應(yīng)度是否達(dá)到要求，否則，執(zhí)行步驟( 3 )。

( 3 ) 生成火花。根據(jù)煙花算法原理，結(jié)合開采沉陷預(yù)計(jì)參數(shù)特點(diǎn)，分別生成爆炸火花和高斯變異火花，其中高斯變異火花增加了火花種群的多樣性，使得算法尋優(yōu)能力進(jìn)一步提高。

( 4 ) 火花映射。步驟( 3 )生成的火花可能會(huì)超出可行域邊界，通過映射規(guī)則把超出邊界的火花映射到參數(shù)范圍內(nèi)。

( 5 ) 后代選擇。種群中適應(yīng)度最小的候選者會(huì)被確定性地選擇到下一代中，剩下的后代采用輪盤賭法在候選者中進(jìn)行選擇，候選者xi被選中的概率計(jì)算公式為

式中，R( xi)為當(dāng)前個(gè)體xi與其他個(gè)體的距離之和；d( xi，yi)為任意2個(gè)個(gè)體xi和yi之間的歐氏距離；xj∈K表示第j個(gè)位置屬于集合K。

( 6 ) 迭代計(jì)算。循環(huán)執(zhí)行步驟( 2 )～( 5 )，直至達(dá)到適應(yīng)度要求，輸出參數(shù)最優(yōu)解。

基于上述思想，提出了如圖2所示的基于變權(quán)PB組合預(yù)計(jì)模型的開采沉陷預(yù)計(jì)參數(shù)反演方法技術(shù)路線，并編制了求參程序。

3 模擬試驗(yàn)

3.1 數(shù)據(jù)模擬

圖2 構(gòu)建方法技術(shù)路線Fig.2 Technology roadmap of constructed method

以淮南礦區(qū)采礦背景為例，模擬工作面上覆巖層巖性為中硬，煤層采厚m＝3 m，采高H＝400 m，煤層傾角α＝6°，工作面傾向方位角θ＝0°；模擬工作面尺寸D3×D1＝800 m×300 m，工作面走向?yàn)槌浞植蓜?dòng)，傾向?yàn)榉浅浞植蓜?dòng)，整體為非充分采動(dòng)。此外，分別沿工作面走向和傾向模擬布設(shè)了2條監(jiān)測(cè)線，走向共布設(shè)35個(gè)監(jiān)測(cè)點(diǎn)( E1～E35 )，傾向共布設(shè)23個(gè)監(jiān)測(cè)點(diǎn)( N1～N23 )，監(jiān)測(cè)點(diǎn)間距均為40 m；模擬工作面及監(jiān)測(cè)線如圖3所示。模擬的開采沉陷預(yù)計(jì)參數(shù)分別為：q＝0.8，tan β＝2.0，b＝0.3，θ＝87.0°，S1＝S2＝S3＝S4＝20 m；模擬權(quán)值ρ＝0.6。根據(jù)模擬的地質(zhì)采礦概況參數(shù)和預(yù)計(jì)參數(shù)，利用組合預(yù)計(jì)模型預(yù)計(jì)走向和傾向模擬監(jiān)測(cè)點(diǎn)的下沉和水平移動(dòng)。

圖3 模擬工作面和觀測(cè)線示意Fig.3 Schematic diagram of simulated working face and observation lines

3.2 反演及結(jié)果分析

以3.1節(jié)中預(yù)計(jì)的地表形變作為參數(shù)反演基礎(chǔ)數(shù)據(jù)，結(jié)合模擬的工作面地質(zhì)概況參數(shù)，利用本文構(gòu)建的方法對(duì)工作面進(jìn)行開采沉陷預(yù)計(jì)參數(shù)反演試驗(yàn)，為了排除偶然誤差的干擾，分別獨(dú)立進(jìn)行10次試驗(yàn)，10次試驗(yàn)求參取平均值作為最終試驗(yàn)結(jié)果，并計(jì)算求參相對(duì)誤差及中誤差，試驗(yàn)結(jié)果見表1。

表1 模擬試驗(yàn)求參結(jié)果Table 1 Results of the solved parameters in simulation experiment

從表1中的數(shù)據(jù)可以看出，① 構(gòu)建方法反演q，tan β，b，θ的相對(duì)誤差控制在1.7%以內(nèi)，以及最大參數(shù)中誤差不超過1.6；反演拐點(diǎn)偏移距S的最大相對(duì)誤差約為19.6%( 在傾向下山處取得 )，最大參數(shù)中誤差不超過8.6 m。② 構(gòu)建方法求取ρ的值為0.507 3，與模擬值比較，相對(duì)誤差不超過15.5%。③構(gòu)建方法擬合下沉和水平移動(dòng)誤差在-22.3～30.2 mm，擬合中誤差約為12.3 mm，工作面走向和傾向下沉擬合效果如圖4，5所示；工作面走向和傾向水平移動(dòng)擬合效果如圖6，7所示。綜上，本文構(gòu)建方法能夠精確可靠地反演開采沉陷預(yù)計(jì)參數(shù)。

圖4 走向方向?qū)崪y(cè)和擬合下沉曲線比較Fig.4 Comparison of measured and fitted subsidence curves along the strike

4 工程應(yīng)用

4.1 試驗(yàn)區(qū)概況

圖5 傾向方向?qū)崪y(cè)和擬合下沉曲線比較Fig.5 Comparison of measured and fitted subsidence curves along the inclination

圖6 走向方向?qū)崪y(cè)和擬合水平移動(dòng)曲線比較Fig.6 Comparison of measured and fitted horizontal movement curves along the strike

圖7 傾向方向?qū)崪y(cè)和擬合水平移動(dòng)曲線比較Fig.7 Comparison of measured and fitted horizontal movement curves along the inclination

以淮南顧橋煤礦某工作面地表為目標(biāo)研究區(qū)域，該工作面采用綜合機(jī)械化采煤，一次采全高，全部垮落法管理頂板。工作面煤層平均采高為3 m，平均埋深為735 m，傾角平均為5°，工作面開采尺寸為2 120 m×241 m，走向?yàn)槌浞植蓜?dòng)，傾向?yàn)榉浅浞植蓜?dòng)，整體為非充分采動(dòng)。為了動(dòng)態(tài)監(jiān)測(cè)地表移動(dòng)與變形情況，沿走向觀測(cè)線布設(shè)3個(gè)控制點(diǎn)和95個(gè)監(jiān)測(cè)點(diǎn)，點(diǎn)間距為30/60 m，長(zhǎng)度為3 480 m；沿傾向布設(shè)3個(gè)控制點(diǎn)和50個(gè)監(jiān)測(cè)點(diǎn)，長(zhǎng)度為1 500 m；地表移動(dòng)變形穩(wěn)定后，由于后期變形盆地積水以及點(diǎn)位破壞等情況，部分點(diǎn)位不能達(dá)到使用要求，在本試驗(yàn)中予以刪除。

4.2 試驗(yàn)過程及結(jié)果分析

為了比較組合預(yù)計(jì)模型求參方法和單個(gè)預(yù)計(jì)模型求參方法擬合及求參效果，基于試驗(yàn)區(qū)地質(zhì)采礦概況和地表觀測(cè)站監(jiān)測(cè)形變，利用本文構(gòu)建方法、基于概率積分模型的求參方法( 方法Ⅰ )和基于Boltzmann函數(shù)預(yù)計(jì)模型的求參方法( 方法Ⅱ )分別對(duì)該工作面開采沉陷預(yù)計(jì)參數(shù)進(jìn)行求取，求參試驗(yàn)結(jié)果見表2。

表2 工程應(yīng)用求參結(jié)果比較Table 2 Comparison of the solved parameters in engineering applications

由表2可知，① 就求參差異來(lái)看，本文構(gòu)建方法與方法Ⅰ相比，求參差異主要集中在下沉系數(shù)q和下山拐點(diǎn)偏移距S1處，參數(shù)q的差異百分比超過9%，參數(shù)S1的求參差異超過14 m；本文構(gòu)建方法與方法Ⅱ相比，求參差異主要集中在主要影響傳播角θ以及下山和上山拐點(diǎn)偏移距S1，S2處，求參差異分別約為2.1°，9 m和13 m。② 工作面走向和傾向下沉和水平移動(dòng)擬合效果分別如圖8，9所示，本文構(gòu)建方法擬合優(yōu)于方法Ⅰ和方法Ⅱ( 相較于藍(lán)點(diǎn)和綠點(diǎn)，紅點(diǎn)接近于黑線的數(shù)量較多 )，3種方法擬合下沉和水平移動(dòng)中誤差分別約為82，107，100 mm。③ 求取權(quán)值ρ約為0.7，表明概率積分模型參與擬合的形變比例大于Boltzmann函數(shù)預(yù)計(jì)模型。試驗(yàn)結(jié)果表明，本文構(gòu)建方法綜合了方法Ⅰ和方法Ⅱ的擬合形變優(yōu)點(diǎn)，能夠精確可靠地反演開采沉陷預(yù)計(jì)參數(shù)。

5 結(jié) 論

圖8 工程應(yīng)用走向方向擬合下沉和水平移動(dòng)曲線及誤差比較Fig.8 Comparison of curves and errors of fitting sinking and horizontal movement along the strike in engineering applications

圖9 工程應(yīng)用傾向方向擬合下沉和水平移動(dòng)曲線及誤差比較Fig.9 Comparison of curves and errors of fitting sinking and horizontal movement along the inclination in engineering applications

( 1 ) 基于變權(quán)組合理論，融合概率積分模型和Boltzmann模型，構(gòu)建了變權(quán)PB組合預(yù)計(jì)模型，考慮 到該模型屬于高度非線性模型，耦合煙花算法，提出了基于變權(quán)PB組合預(yù)計(jì)模型的開采沉陷預(yù)計(jì)參數(shù)反演方法。

( 2 ) 模擬試驗(yàn)結(jié)果表明，本文構(gòu)建方法反演全部參數(shù)的平均相對(duì)誤差約為4.0%，參數(shù)中誤差約為3.9；構(gòu)建方法擬合下沉和水平移動(dòng)誤差為-22.3～30.2 mm，擬合中誤差約為12.3 mm；構(gòu)建方法能夠精確可靠地反演概率積分模型參數(shù)。

( 3 ) 將本文構(gòu)建方法、基于概率積分模型的求參方法和基于Boltzmann函數(shù)模型的求參方法應(yīng)用到淮南顧橋礦某工作面的開采沉陷預(yù)計(jì)模型參數(shù)反演中，試驗(yàn)結(jié)果表明，本文構(gòu)建方法擬合程度遠(yuǎn)優(yōu)于其他2種方法，3種方法擬合下沉和水平移動(dòng)中誤差分別約為82，107，100 mm，求取權(quán)值ρ約為0.7，表明概率積分模型參與擬合的形變比例大于Boltzmann函數(shù)模型，構(gòu)建方法求參結(jié)果為：q＝0.99，tan β＝1.90，b＝0.42，θ＝87.03°，S1＝6.29 m，S2＝-29.52 m，S3＝49.78 m，S4＝32.56 m。

參考文獻(xiàn)( References )：

[1] 郭廣禮，汪云甲.概率積分法參數(shù)的穩(wěn)健估計(jì)模型及其應(yīng)用研究[J].測(cè)繪學(xué)報(bào)，2000，29( 2 )：162-165.GUO Guangli，WANG Yunjia.Study of robust determining parameters model for probability-integral method and its application[J].Acta Geodaetica et Cartographica Sinica，2000，29( 2 )：162-165.

[2] 沈震，徐良驥，劉哲，等.基于Matlab的概率積分法開采沉陷預(yù)計(jì)參數(shù)解算[J].金屬礦山，2015，44( 9 )：170-174.SHEN Zhen，XU Liangji，LIU Zhe，et al.Calculating on the prediction parameters of mining subsidence with probability integral method based on Matlab[J].Metal Mine，2015，44( 9 )：170-174.

[3] 吳侃，周鳴.礦區(qū)沉陷預(yù)測(cè)預(yù)報(bào)系統(tǒng)[M].徐州：中國(guó)礦業(yè)大學(xué)出版社，1999.WU Kan，ZHOU Ming.Subsidence prediction system of mining area[M].Xuzhou：China University of Mining and Technology Press，1999.

[4] 葛家新.地表沉陷預(yù)計(jì)參數(shù)求取及其分析[J].礦山壓力與頂板管理，2004，21( 1 )：78-79.GE Jiaxin.Surface subsidence prediction parameters evaluation and analysis[J].Ground Pressure and Strata Control，2004，21( 1 )：78-79.

[5] 查劍鋒，馮文凱，朱曉峻.基于遺傳算法的概率積分法預(yù)計(jì)參數(shù)反演[J].采礦與安全工程學(xué)報(bào)，2011，28( 4 )：655-659.ZHA Jianfeng，F(xiàn)ENG Wenkai，ZHU Xiaojun.Research on parameters inversion in probability integral method by genetic algorithm[J].Journal of Mining & Safety Engineering，2011，28( 4 )：655-659.

[6] 徐孟強(qiáng)，查劍鋒，李懷展.基于PSO算法的概率積分法預(yù)計(jì)參數(shù)反演[J].煤炭工程，2015，47( 7 )：117-119，123.XU Mengqiang，ZHA Jianfeng，LI Huaizhan.Parameters inversion in probability integral method by particle swarm optimization[J].Coal Engineering，2015，47( 7 )：117-119，123.

[7] 魏宗海.粒子群算法在概率積分法沉陷預(yù)計(jì)模型參數(shù)反演中的應(yīng)用[J].測(cè)繪工程，2017，26( 10 )：36-39.WEI Zonghai.Application of particle swarm optimization in parameter inversion of probabilistic integral subsidence prediction model[J].Engineering of Surveying and Mapping，2017，26( 10 )：36-39.

[8] 劉 奇，朱建軍，蘇軍明，等.基于改進(jìn)人工蜂群算法的概率積分法參數(shù)反演[J].測(cè)繪工程，2017，26( 10 )：18-21，29.LIU Qi，ZHU Jianjun，SU Junming，et al.Parameters inversion in probability integral method based on self-adaptive artificial bee colony[J].Engineering of Surveying and Mapping，2017，26( 10 )：18-21，29.

[9] 陳濤，郭廣禮，朱曉峻，等.利用果蠅算法反演概率積分法開采沉陷預(yù)計(jì)參數(shù)[J].金屬礦山，2016，46( 6 )：185-188.CHEN Tao，GUO Guangli，ZHU Xiaojun，et al.Mining subsidence prediction parameters inversion of the probability integral method based on fruit flies algorithm[J].Metal Mine，2016，46( 6 )：185-188.

[10] 楊靖宇，劉超，王彬.基于BFGS算法的概率積分模型的參數(shù)反演[J].煤炭學(xué)報(bào)，2019，44( 10 )：3058-3068.YANG Jingyu，LIU Chao，WANG Bin.BFGS method based inversion of parameters in probability integral model[J].Journal of China Coal Society，2019，44( 10 )：3058-3068.

[11] 王金濤，郭廣禮，郭慶彪，等.相關(guān)性的組合預(yù)測(cè)模型在開采沉陷中的應(yīng)用[J].測(cè)繪科學(xué)，2017，42( 10 )：181-185.WANG Jintao，GUO Guangli，GUO Qingbiao，et al.A combined model for mining subsidence forecasting based on correlativity[J].Science of Surveying and Mapping，2017，42( 10 )：181-185.

[12] 羅亦泳，姚宜斌，趙慶志，等.利用優(yōu)化的組合核相關(guān)向量機(jī)算法構(gòu)建地表下沉預(yù)測(cè)模型[J].武漢大學(xué)學(xué)報(bào)( 信息科學(xué)版 )，2018，43( 9 )：1295-1301.LUO Yiyong，YAO Yibin，ZHAO Qingzhi，et al.Prediction of surface subsidence of underground mining based on HIOA and MK-RVM[J].Geomatics and Information Science of Wuhan University，2018，43( 9 )：1295-1301.

[13] 魯小紅.最優(yōu)權(quán)組合預(yù)測(cè)法在采煤沉陷變形預(yù)測(cè)中的應(yīng)用[J].測(cè)繪通報(bào)，2020( 4 )：111-115.LU Xiaohong.Application of optimal weight combination method in predicting coal mining subsidence deformation[J].Bulletin of Surveying and Mapping，2020( 4 )：111-115.

[14] 王寧，吳侃，劉錦，等.基于Boltzmann函數(shù)的開采沉陷預(yù)測(cè)模型[J].煤炭學(xué)報(bào)，2013，38( 8 )：1352-1356.WANG Ning，WU Kan，LIU Jin，et al.Model for mining subsidence prediction based on Boltzmann function[J].Journal of China Coal Society，2013，38( 8 )：1352-1356.

[15] YANG Z F，LI Z W，ZHU Jianjun，et al.An extension of the insarbased probability integral method and its application for predicting 3-d mining-induced displacements under different extraction conditions[J].IEEE Transactions on Geoscience and Remote Sensing，2017，99：1-11.

[16] 何國(guó)清，楊倫，凌賡娣，等.礦山開采沉陷學(xué)[M].徐州：中國(guó)礦業(yè)大學(xué)出版社，1991.HE Guoqing，YANG Lun，LING Gengdi，et al.Mining subsidence in mine[M].Xuzhou：China University of Mining and Technology Press，1991.

[17] 劉玉成，曹樹剛.基于關(guān)鍵層理論的地表下沉盆地模型初探[J].巖土力學(xué)，2012，33( 3 )：719-724.LIU Yucheng，CAO Shugang.Preliminary study of surface subsidence model based on theory of key rock stratum[J].Rock and Soil Mechanics，2012，33( 3 )：719-724.

[18] FAN H D，CHENG D，DENG K Z，et al.Subsidence monitoring using D-InSAR and probability integral prediction modelling in deep mining areas[J].Survey Review，2015，47( 345 )：1752270614Y.000.

[19] TAN Y，ZHU Y.Fireworks algorithm for optimization[J].In Lecture Notes in Computer Science ( including subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics )，2010，6145 LNCS：355-364.

[20] 譚營(yíng)，鄭少秋.煙花算法研究進(jìn)展[J].智能系統(tǒng)學(xué)報(bào)，2014，9( 5 )：515-528.TAN Ying ，ZHENG Shaoqiu.Recent advances in fireworks algorithm[J].CAAI Transactions on Intelligent Systems，2014，9( 5 )：515-528.

[21] SUDHAKAR Babu T，PRASANTH Ram J，SANGEETHA K，et al.Parameter extraction of two diode solar PV model using Fireworks algorithm[J].Solar Energy，2016( 140 )：265-276.

[22] EI MAJDOULI M A，RBOUH I，BOUGRINE S，et al.Fireworks algorithm framework for big data optimization[J].Memetic Computing，2016，8( 4 )：333-347.