任 敏 朱萬成 賈瀚文 程關(guān)文 李斕堃 代 風(fēng)

（東北大學(xué)資源與土木工程學(xué)院，遼寧沈陽(yáng)110819）

目前，地下采空區(qū)引起的地表沉降逐漸成為了一種常見地質(zhì)災(zāi)害。以鞍鋼弓長(zhǎng)嶺露天鐵礦為例，早在日偽時(shí)期就進(jìn)行過地下富礦開采，造成地下多處采空區(qū)。隨著露天采場(chǎng)開采范圍拓展與延伸及深部鐵礦資源的開采和利用，地表沉降嚴(yán)重威脅著礦山大型設(shè)備的安全，給正常采礦生產(chǎn)帶來了極大困難。為了能夠?qū)崟r(shí)掌握弓長(zhǎng)嶺露天鐵礦淺埋采空區(qū)頂板巖體的穩(wěn)定性和地表沉降的發(fā)展演化，研究人員在采空區(qū)頂板和露天坑底的地表建立了集錨桿應(yīng)力計(jì)、多點(diǎn)位移計(jì)、液壓水準(zhǔn)儀為一體的聯(lián)合監(jiān)測(cè)系統(tǒng)。因此，基于該聯(lián)合監(jiān)測(cè)系統(tǒng)的監(jiān)測(cè)信息，分析淺埋采空區(qū)頂板巖體內(nèi)應(yīng)力、位移的變化與地表沉降之間的影響關(guān)系與響應(yīng)機(jī)制，對(duì)于多傳感器監(jiān)測(cè)信息有效融合和精準(zhǔn)的預(yù)測(cè)預(yù)警模型構(gòu)建意義重大。

為避免地表沉降等災(zāi)害的發(fā)生，越來越多的監(jiān)測(cè)技術(shù)與監(jiān)測(cè)設(shè)備應(yīng)運(yùn)而生，為地表沉降監(jiān)測(cè)提供了多種多樣的監(jiān)測(cè)手段和監(jiān)測(cè)方法。目前普遍常用的地表沉降監(jiān)測(cè)手段主要有全站儀［1］、水準(zhǔn)儀［2］、GPS站點(diǎn)［3］等定點(diǎn)式監(jiān)測(cè)方法，以及三維激光掃描［4］、雷達(dá)掃描［5］等區(qū)域式監(jiān)測(cè)方法。在獲取現(xiàn)場(chǎng)實(shí)時(shí)監(jiān)測(cè)數(shù)據(jù)的前提下，為避免地表沉降變形引起災(zāi)害及人員、財(cái)產(chǎn)損失，開展地表沉降預(yù)測(cè)研究十分必要。近年來，眾多研究者致力于利用不同的方法手段實(shí)現(xiàn)這一目標(biāo)，主要包括人工智能法、灰色預(yù)測(cè)法、數(shù)學(xué)函數(shù)法等。如李勇發(fā)等［6］基于PS-InSAR 雷達(dá)干涉測(cè)量技術(shù)對(duì)地表沉降進(jìn)行了監(jiān)測(cè)，對(duì)獲取的實(shí)時(shí)沉降數(shù)據(jù)利用遺傳算法優(yōu)化的反饋神經(jīng)網(wǎng)絡(luò)算法進(jìn)行了地表沉降預(yù)測(cè)，獲得了良好的預(yù)測(cè)效果；張濤［7］將DInSAR 技術(shù)與小基線集技術(shù)進(jìn)行有機(jī)結(jié)合對(duì)地表沉降進(jìn)行了監(jiān)測(cè)，并使用粒子群算法優(yōu)化支持向量機(jī)模型對(duì)沉降數(shù)據(jù)進(jìn)行了擬合與預(yù)測(cè)；黨星海等［8］采用卡爾曼濾波方法對(duì)三維激光掃描沉降監(jiān)測(cè)獲取的初始監(jiān)測(cè)數(shù)據(jù)進(jìn)行了去噪處理，而后利用GM（1，1）灰色預(yù)測(cè)模型進(jìn)行了沉降變形預(yù)測(cè)，得到了優(yōu)于傳統(tǒng)灰色預(yù)測(cè)模型的預(yù)測(cè)結(jié)果；李篷［9］采用GPS方法進(jìn)行了煤礦地表非采動(dòng)沉降監(jiān)測(cè)，并利用GM（1，1）灰色預(yù)測(cè)模型對(duì)煤礦主井、副井、煤倉(cāng)附近區(qū)域開展了地表沉降預(yù)測(cè)；Nie 等［10-11］基于采空區(qū)上覆巖體不同破裂帶的“S”型變形特征，以及塌陷坑GPS 監(jiān)測(cè)獲取的“S”型地表沉降曲線，提出了采空區(qū)的反正切函數(shù)預(yù)測(cè)模型；袁興明等［12］利用數(shù)字水準(zhǔn)儀進(jìn)行了沉降觀測(cè)，并分別利用BP 神經(jīng)網(wǎng)絡(luò)算法、改進(jìn)BP 神經(jīng)網(wǎng)絡(luò)算法對(duì)沉降數(shù)據(jù)進(jìn)行了預(yù)測(cè)。就目前研究來看，基于單一監(jiān)測(cè)信息的地表沉降預(yù)測(cè)存在一定的缺陷，是因?yàn)閿?shù)據(jù)源本身既無法充分反映地表沉降的多種誘導(dǎo)因素，也無法全面體現(xiàn)沉降過程中的多種響應(yīng)信息，從而導(dǎo)致預(yù)測(cè)效果不穩(wěn)定且與實(shí)際存在一定的偏差。

為了獲取更全面、精準(zhǔn)的監(jiān)測(cè)數(shù)據(jù)，從多種信息融合角度判斷監(jiān)測(cè)對(duì)象的發(fā)展?fàn)顟B(tài)，采用多種監(jiān)測(cè)手段針對(duì)地表沉降變形建立聯(lián)合監(jiān)測(cè)系統(tǒng)，開展多手段聯(lián)合監(jiān)測(cè)已經(jīng)得到廣泛研究。如何倩等［13］采用三維激光掃描與D-InSAR 聯(lián)合監(jiān)測(cè)礦區(qū)地表動(dòng)態(tài)沉降；臧巍［14］將測(cè)量機(jī)器人與電子水準(zhǔn)儀相結(jié)合應(yīng)用于隧道監(jiān)測(cè)工程中；郭利民等［15］利用GPS 和PS-In-SAR 聯(lián)合監(jiān)測(cè)西秦嶺斷裂北緣地表運(yùn)動(dòng)；李楊楊等［16］將GPS-RTK聯(lián)合水準(zhǔn)儀應(yīng)用于充填工作面地表高精度動(dòng)態(tài)監(jiān)測(cè)中。目前，學(xué)術(shù)界對(duì)于多手段聯(lián)合監(jiān)測(cè)數(shù)據(jù)的處理及應(yīng)用方面存在兩種研究方向：一是將多種監(jiān)測(cè)技術(shù)進(jìn)行有機(jī)結(jié)合，得到精度較高的地表變形監(jiān)測(cè)結(jié)果，如余禮仁［17］在對(duì)淮北袁店二礦煤礦區(qū)地表沉降進(jìn)行監(jiān)測(cè)時(shí)，結(jié)合了D-InSAR 和GIS技術(shù)得到監(jiān)測(cè)數(shù)據(jù)，而后采用三次樣條插值與灰色預(yù)測(cè)組合模型GM（1，1）對(duì)該礦7255工作面中心區(qū)域的下沉趨勢(shì)進(jìn)行了擬合并預(yù)測(cè)；邢學(xué)敏［18］將CR-In-SAR 與PS-InSAR 相結(jié)合，對(duì)礦區(qū)時(shí)序地表變形進(jìn)行了監(jiān)測(cè)，提高了礦區(qū)地表變形監(jiān)測(cè)的精度和可靠性；二是根據(jù)模糊系統(tǒng)理論，將多個(gè)傳感器獲取的信息進(jìn)行融合，進(jìn)而對(duì)地表沉降狀態(tài)進(jìn)行綜合估計(jì)與評(píng)價(jià)，如Wu 等［19］采用野外調(diào)查、地球物理勘探、鉆井勘探和D-InSAR 技術(shù)等聯(lián)合監(jiān)測(cè)手段，應(yīng)用模糊綜合評(píng)價(jià)法對(duì)地表沉降進(jìn)行了風(fēng)險(xiǎn)分區(qū)。由此可見，對(duì)于監(jiān)測(cè)數(shù)據(jù)而言，雖然多傳感器聯(lián)合監(jiān)測(cè)總體上提高了數(shù)據(jù)樣本的多樣性和精度，但是在基于多傳感器聯(lián)合監(jiān)測(cè)結(jié)果的預(yù)測(cè)方法研究與應(yīng)用中，對(duì)多傳感器監(jiān)測(cè)獲得的地表變形時(shí)序數(shù)據(jù)未能進(jìn)行有效融合，對(duì)于提升預(yù)測(cè)預(yù)警效果作用不顯著。

程關(guān)文等［20-21］認(rèn)為深部巖體變形是形成地表變形區(qū)域的原因，并采用數(shù)值模擬方法分析了隨著不同開挖步的增加，深部巖體變形產(chǎn)生地表沉降變形的整個(gè)過程。宋長(zhǎng)根等［22］對(duì)深部巖體變形對(duì)地表沉降變形的影響進(jìn)行了研究，認(rèn)為當(dāng)深部巖體發(fā)生破壞后，地表會(huì)進(jìn)入加速變形階段，地表與構(gòu)筑物不久便出現(xiàn)裂縫，嚴(yán)重影響了礦區(qū)安全生產(chǎn)與構(gòu)筑物正常使用。本研究以弓長(zhǎng)嶺露天鐵礦某空區(qū)頂板為例，根據(jù)錨桿應(yīng)力計(jì)、多點(diǎn)位移計(jì)、液壓水準(zhǔn)儀對(duì)空區(qū)頂板巖體及地表沉降的聯(lián)合監(jiān)測(cè)結(jié)果，采用自回歸分布滯后模型（ARDL）分析淺埋采空區(qū)頂板巖體內(nèi)的應(yīng)力、位移變化對(duì)地表沉降的影響關(guān)系和響應(yīng)機(jī)制，確定該空區(qū)頂板巖體內(nèi)發(fā)生應(yīng)力、位移變化影響到地表沉降的滯后時(shí)間，這種發(fā)展關(guān)系的探究不僅可以為多傳感器的信息融合提供有效途徑，還可以為地表沉降災(zāi)害的預(yù)測(cè)預(yù)警提供支撐。

1 弓長(zhǎng)嶺鐵礦某采空區(qū)概況

近年來，弓長(zhǎng)嶺露采分公司采用地球物理探測(cè)方法研究采空區(qū)的空間分布和大小情況，取得了顯著成效，基本查清了一定深度空區(qū)的分布、大小等信息。但是，僅僅研究采空區(qū)的分布、大小、深度遠(yuǎn)遠(yuǎn)無法滿足礦山安全生產(chǎn)的需要。因此，結(jié)合弓長(zhǎng)嶺鐵礦的工程實(shí)際情況，建立了一套先進(jìn)的監(jiān)測(cè)系統(tǒng)，包括錨桿測(cè)力計(jì)、多點(diǎn)位移計(jì)、液體靜力水準(zhǔn)儀等，實(shí)現(xiàn)對(duì)地表沉降變形和采空區(qū)頂板巖體的實(shí)時(shí)在線監(jiān)測(cè)，掌握空區(qū)頂板巖體內(nèi)應(yīng)力、位移變化及地表沉降的發(fā)展，以確定機(jī)械設(shè)備、車輛能否安全通過，對(duì)于露天礦安全開采至關(guān)重要。系統(tǒng)測(cè)點(diǎn)布置如圖1所示。

本研究選擇2018-06-11—2018-07-01獲取的時(shí)間序列監(jiān)測(cè)數(shù)據(jù)進(jìn)行分析，測(cè)點(diǎn)為研究區(qū)內(nèi)（圖1）地表沉降液壓水準(zhǔn)儀測(cè)點(diǎn)S4、A 孔多點(diǎn)位移計(jì)位于135.2 m 水平的測(cè)點(diǎn)uA4和位于129.2 m 水平的測(cè)點(diǎn)uA1，A 孔錨桿測(cè)力計(jì)位于138.7 m 水平的測(cè)點(diǎn)σA3和位于131.7 m 水平的測(cè)點(diǎn)σA1。選擇該類測(cè)點(diǎn)進(jìn)行分析，原因在于：①3 種傳感器測(cè)點(diǎn)基本處于同一地質(zhì)環(huán)境，上覆巖體中積聚產(chǎn)生的應(yīng)力和變形會(huì)以較單一的形式影響到地表；②各個(gè)測(cè)點(diǎn)空間距離較近，且在垂直方向上多點(diǎn)位移計(jì)測(cè)點(diǎn)與錨桿測(cè)力計(jì)測(cè)點(diǎn)在采空區(qū)頂板巖體中交錯(cuò)布置，基本能夠監(jiān)測(cè)到巖體內(nèi)發(fā)生的應(yīng)力和變形進(jìn)而影響地表產(chǎn)生沉降的全過程。因此，可根據(jù)該類測(cè)點(diǎn)的數(shù)據(jù)，研究上覆巖體內(nèi)的應(yīng)力、位移變化對(duì)地表沉降的影響機(jī)制。

2 模型構(gòu)建及實(shí)證分析

2.1 ARDL模型設(shè)定

本研究采用自回歸分布滯后模型（ARDL）分析由錨桿應(yīng)力計(jì)、多點(diǎn)位移計(jì)監(jiān)測(cè)的淺埋采空區(qū)頂板巖體內(nèi)的應(yīng)力、位移對(duì)液壓水準(zhǔn)儀監(jiān)測(cè)獲得的地表沉降之間的影響關(guān)系、響應(yīng)機(jī)制和滯后時(shí)間。Pesaran等［23］認(rèn)為ADRL 模型的優(yōu)勢(shì)在于不需要變量同階單整，也可以檢驗(yàn)變量之間的長(zhǎng)期關(guān)系，當(dāng)解釋變量為內(nèi)生變量時(shí)，ARDL 模型的估計(jì)也不會(huì)受到影響，并且在小樣本情況下ARDL的估計(jì)結(jié)果也更加穩(wěn)定?；诖耍狙芯繕?gòu)建的地表沉降與頂板巖體應(yīng)力、位移之間的ARDL模型為

式中，S 為地表沉降，為被解釋變量；u 和σ 分別為采空區(qū)頂板巖體的位移和應(yīng)力，為解釋變量；q1、q2、q3分別為各變量對(duì)地表沉降長(zhǎng)期影響的滯后期數(shù)；α1、α2、α3為各變量及其滯后期數(shù)對(duì)被解釋變量的影響程度，若解釋變量的滯后期t 值顯著，則該變量的滯后期對(duì)地表沉降有顯著影響；β0為常數(shù)項(xiàng)；ut為殘差項(xiàng)。

2.2 基于ARDL模型實(shí)證分析

2.2.1 單位根檢驗(yàn)

為了避免在實(shí)證檢驗(yàn)過程中發(fā)生偽回歸，需要首先對(duì)數(shù)據(jù)進(jìn)行單位根檢驗(yàn)。盡管ARDL 模型不要求變量具有同階單整，但應(yīng)用該模型的前提條件需確保變量單整階數(shù)不能超過1。本研究采用最常用的ADF 檢驗(yàn)統(tǒng)計(jì)量對(duì)模型變量進(jìn)行單位根檢驗(yàn)，結(jié)果如表1所示。

注：?jiǎn)挝桓鶛z驗(yàn)形式（C，T，L）中，C為截距項(xiàng)，T為趨勢(shì)項(xiàng)，L為滯后階數(shù)。

由表1可知，在5%的顯著性水平下，包括采空區(qū)頂板巖體的位移、應(yīng)力以及地表沉降在內(nèi)的所有變量均為非平穩(wěn)序列。對(duì)這3 個(gè)變量進(jìn)行一階差分后的ADF 檢驗(yàn)，在5%的顯著性水平上都是平穩(wěn)序列，說明這3 個(gè)變量均是一階單整（即I(1)）序列，符合ARDL模型的應(yīng)用要求。

2.2.2 協(xié)整檢驗(yàn)

ARDL 模型的第一步是檢驗(yàn)變量的協(xié)整關(guān)系。由于協(xié)整關(guān)系對(duì)滯后階數(shù)的選擇非常敏感，因此，需要估計(jì)向量自回歸模型（Vector Autoregression，VAR），確定滯后階數(shù)。通過信息準(zhǔn)則和似然比（LR）檢驗(yàn)法得到VAR 模型滯后階數(shù)的確定結(jié)果如表2 所示。在滯后階數(shù)為1 時(shí)，最終預(yù)測(cè)誤差FPE 和3 個(gè)信息準(zhǔn)則AIC、SC、HQ 均最小，由4 個(gè)參考值確定1 階滯后為最優(yōu)，因此可得滯后階數(shù)為1。

注：LogL 表示對(duì)數(shù)似然函數(shù)；表中用“*”表示從每一列標(biāo)準(zhǔn)中選的滯后階數(shù)；

滯后階數(shù)確定后，為方便起見，本研究采用Johansen 協(xié)整檢驗(yàn)方法進(jìn)行協(xié)整檢驗(yàn)，結(jié)果如表3 所示。

?

由表3可知：地表沉降S4與采空區(qū)頂板巖體內(nèi)的位移uA1、uA4和應(yīng)力σA1、σA3之間跡統(tǒng)計(jì)量和最大特征值統(tǒng)計(jì)量的檢驗(yàn)結(jié)果均存在協(xié)整關(guān)系，因而可得出地表沉降與淺埋采空區(qū)頂板巖層中的位移、應(yīng)力之間存在長(zhǎng)期關(guān)系。

2.2.3 長(zhǎng)期均衡系數(shù)估計(jì)

確定各變量之間存在協(xié)整關(guān)系后，可以用ARDL模型進(jìn)行長(zhǎng)期關(guān)系和短期關(guān)系的估計(jì)。在估計(jì)長(zhǎng)期關(guān)系的變量系數(shù)之前需要確定模型中的滯后階數(shù)，本研究使用SC準(zhǔn)則確定模型中各變量的最優(yōu)滯后階數(shù)。估計(jì)結(jié)果顯示：地表沉降S4與位移uA1、uA4，應(yīng)力σA1、σA3以ARDL( 1,1 ,1,1,1 )模型最為合適，長(zhǎng)期系數(shù)如表4所示。

注：D 表示差分，是因?yàn)閱挝桓鶛z驗(yàn)已得到的各變量均為一階差分平穩(wěn)，故所有變量均采用一階差分形式；表示一階差分后的各變量由ARDL 模型選擇的滯后一階形式表示一階差分的原變量；C 為截距項(xiàng)；@TREND 為時(shí)間趨勢(shì)項(xiàng)；AIC 為赤池AIC信息準(zhǔn)則（Akaike Information Criterion），SC為施瓦茨SIC信息準(zhǔn)側(cè)（Schwarz Criterion），兩者均是衡量回歸模型優(yōu)良性的標(biāo)準(zhǔn)，數(shù)值較小意味著滯后階數(shù)較為合適；R2 統(tǒng)計(jì)量用于衡量模型的估計(jì)值對(duì)實(shí)際值擬合的好壞，數(shù)值越大，表明模型對(duì)因變量的擬合效果越好。

由表4 可知：淺埋采空區(qū)頂板巖體內(nèi)的應(yīng)力、位移滯后量對(duì)地表沉降影響顯著，說明上覆巖體內(nèi)的應(yīng)力、變形影響傳遞到地表需要時(shí)間，以使得滯后量與地表沉降最相關(guān)。變量系數(shù)有正有負(fù)，變量系數(shù)為正，說明采空區(qū)頂板巖體內(nèi)產(chǎn)生的應(yīng)力、位移使得地表沉降加?。蛔兞肯禂?shù)為負(fù)則較為復(fù)雜，存在水平地應(yīng)力增加的情形，使得上覆巖體內(nèi)的應(yīng)力、變形對(duì)地表呈現(xiàn)出與塌陷相反方向變形的影響。估計(jì)得到的地表沉降S4與位移uA1、uA4，應(yīng)力σA1、σA3的ARDL( 1,1 ,1,1,1 )模型中，滯后一期的位移uA4、應(yīng)力σA1和σA3對(duì)地表沉降的影響均顯著，即前1 d 采空區(qū)頂板內(nèi)產(chǎn)生的應(yīng)力、位移與當(dāng)前發(fā)生的地表沉降最相關(guān)，其中應(yīng)力σA1對(duì)地表沉降S4的影響最大，其系數(shù)為0.529 29。

此外，還進(jìn)行了數(shù)值模擬分析，即保證模擬的弓長(zhǎng)嶺露天鐵礦淺埋采空區(qū)頂板巖體所受的豎向應(yīng)力基本控制在監(jiān)測(cè)結(jié)果范圍內(nèi)，測(cè)點(diǎn)所受應(yīng)力為0.38～0.74 MPa 時(shí)，提取地表測(cè)點(diǎn)處的沉降值為0.005 97～0.009 39 m（圖2）。由數(shù)值模擬結(jié)果分析可知：σA1的改變對(duì)地表沉降的影響較大，系數(shù)達(dá)0.572 81 左右，與ARDL 模型的估計(jì)結(jié)果較吻合。故可將淺埋采空區(qū)頂板巖體內(nèi)的應(yīng)力、位移作為地表沉降的長(zhǎng)期影響指標(biāo)，度量地表沉降。

2.2.4 短期動(dòng)態(tài)調(diào)整關(guān)系分析

在估計(jì)得到長(zhǎng)期關(guān)系系數(shù)后，可進(jìn)一步求得基于ARDL模型的誤差修正項(xiàng)（ECM）的系數(shù)，分析各個(gè)解釋變量的短期變動(dòng)對(duì)被解釋變量的動(dòng)態(tài)調(diào)整關(guān)系。ARDL的誤差修正模型具體估計(jì)結(jié)果見表5。

注：EC M( -1 )為滯后一期的誤差修正項(xiàng)。

由表5 可知：在10%的顯著性水平下，滯后一期的地表沉降與滯后一期的位移uA1以及當(dāng)期的應(yīng)力σA1對(duì)地表沉降均存在顯著影響，其中應(yīng)力的影響最為顯著，該變量的正向變動(dòng)會(huì)導(dǎo)致地表沉降加劇。由誤差修正項(xiàng)ECM 估計(jì)結(jié)果表明存在反向調(diào)整機(jī)制，意味著當(dāng)受到短期沖擊后，系統(tǒng)會(huì)較快恢復(fù)至長(zhǎng)期均衡狀態(tài)。

因此，基于上述協(xié)整檢驗(yàn)、ARDL 模型的長(zhǎng)期均衡系數(shù)估計(jì)和短期動(dòng)態(tài)調(diào)整關(guān)系的分析，證實(shí)了地下采空區(qū)頂板巖體內(nèi)產(chǎn)生的應(yīng)力、位移與地表沉降之間存在長(zhǎng)期和短期影響。

3 結(jié)論與建議

在地下開采等多因素的影響下，淺埋采空區(qū)頂板巖體會(huì)產(chǎn)生應(yīng)力、變形，作為直接影響因素導(dǎo)致地表沉降，體現(xiàn)在傳感器上能夠捕捉到的信息即為頂板巖體內(nèi)的應(yīng)力、位移以及地表沉降。為了給地表沉降多傳感器監(jiān)測(cè)信息的有效融合及預(yù)測(cè)預(yù)警模型構(gòu)建提供支撐，本研究利用ARDL模型系統(tǒng)分析了淺埋采空區(qū)頂板巖體內(nèi)應(yīng)力、位移的變化和地表沉降之間的影響關(guān)系。研究表明：地表沉降與淺埋采空區(qū)頂板巖體內(nèi)的應(yīng)力、位移之間存在長(zhǎng)期和短期的影響；利用SC 準(zhǔn)則確定的最優(yōu)模型ARDL( )1,1,1,1,1分析表明，作為地表沉降解釋變量的頂板巖體內(nèi)的應(yīng)力、位移滯后期均為1，即當(dāng)前產(chǎn)生的地表沉降與前1 d 頂板巖體內(nèi)產(chǎn)生的應(yīng)力、位移有關(guān)，且滯后的解釋變量的影響均很顯著。

針對(duì)地表沉降災(zāi)害錨桿應(yīng)力計(jì)、多點(diǎn)位移計(jì)、液壓水準(zhǔn)儀的聯(lián)合監(jiān)測(cè)，本研究提出如下建議：已證實(shí)淺埋采空區(qū)頂板巖體內(nèi)的應(yīng)力、位移對(duì)地表沉降具有長(zhǎng)期影響，且滯后的應(yīng)力、位移影響顯著，故而可以采用各種預(yù)測(cè)模型（如人工智能算法等）將地表沉降時(shí)間序列數(shù)據(jù)作為輸出值，采空區(qū)頂板巖體內(nèi)滯后的應(yīng)力、位移時(shí)間序列數(shù)據(jù)作為輸入值，根據(jù)監(jiān)測(cè)數(shù)據(jù)建立變量之間的線性、非線性關(guān)系進(jìn)行預(yù)測(cè)。基于多傳感器監(jiān)測(cè)數(shù)據(jù)融合構(gòu)建的預(yù)測(cè)模型，可以有效體現(xiàn)出地表沉降過程中淺埋采空區(qū)頂板巖體的響應(yīng)機(jī)理，優(yōu)于具有固定理論模型的數(shù)值模擬方法，因而具有進(jìn)一步深入研究的價(jià)值和應(yīng)用前景。