基于Bertalanffy時間函數(shù)的地表動態(tài)沉陷預(yù)測模型

高 超，徐乃忠，孫萬明，鄧偉男，韓科明

(1.煤炭科學(xué)研究總院 開采研究分院，北京 100013; 2.天地科技股份有限公司 開采設(shè)計事業(yè)部，北京 100013)

地表沉陷是一個非常復(fù)雜的時間-空間運動問題，地表任意時間段內(nèi)的變形量值若達到地表建(構(gòu))筑物的極限承受值，將可能誘發(fā)地質(zhì)災(zāi)害，威脅到地面人員的生命與財產(chǎn)安全。為準(zhǔn)確科學(xué)地評估工作面開采對地面沉陷的影響程度，需要事先建立科學(xué)的地表沉陷動態(tài)沉陷預(yù)計模型，獲得該地區(qū)準(zhǔn)確的模型參數(shù)，進而分析地表任意點、任意時刻的地表動態(tài)沉陷變形量。國內(nèi)外相關(guān)學(xué)者針對不同的動態(tài)沉陷預(yù)計模型進行了探索。

1952年，波蘭學(xué)者KNOTHE[1]基于“某一時刻的地表下沉速度與地表最終下沉值之差成比例”原理推導(dǎo)了Knothe動態(tài)下沉?xí)r間函數(shù)。崔希民等[2]通過實測資料對比研究得出Knothe時間函數(shù)地表移動下沉速度、下沉加速度與地表移動變形的實際過程的差異性。常占強等[3]基于Knothe時間函數(shù)改進并提出了分段Knothe時間函數(shù)模型;后續(xù)張兵等[4]對該函數(shù)進行了優(yōu)化。劉玉成等[5]針對Knothe時間函數(shù)的不足，提出冪指數(shù)-Knothe時間函數(shù)模型。王軍保等[6]將S型生長曲線MMF模型用于地表沉降時間預(yù)測中。王正帥與鄧喀中[7]將Richards生長模型引入地表沉降預(yù)計中，分析了模型參數(shù)與地質(zhì)采礦條件之間的定性、定量關(guān)系。高延法等[8]基于Knothe時間函數(shù)，導(dǎo)出了適用于地表移動衰退期及其以后的地表下沉衰減函數(shù)并對沉陷區(qū)的穩(wěn)定性進行了分類。胡青峰與崔希民等[9]對Knothe時間函數(shù)的各參數(shù)規(guī)律進行了分析。彭小沾與崔希民等[10]對Knothe時間函數(shù)、Sroka-Schober時間函數(shù)的優(yōu)缺點及其相互關(guān)系進行了詳細分析。劉玉成[11]將概率積分法與Weibull時間函數(shù)相結(jié)合，構(gòu)建了沉陷主剖面線的地表動態(tài)下沉函數(shù)模型。朱廣軼等[12]運用概率積分法提出地表動態(tài)移動變形的坐標(biāo)-時間函數(shù)。張文志與鄒友峰等[13]將Logistic模型引入礦區(qū)單點動態(tài)沉陷預(yù)計中，并通過實測數(shù)據(jù)分析了模型的精度與適用性。李德海[14]將覆巖巖性與Knothe時間函數(shù)中的參數(shù)相結(jié)合，探討了覆巖巖性與時間影響參數(shù)的關(guān)系。余闖等[15]將Gompertz函數(shù)引入路堤沉降預(yù)測中，并結(jié)合相關(guān)實踐證明了該函數(shù)的適用性。

沉陷動態(tài)預(yù)計一般由開采引起的最終地表移動變形乘以反映巖層移動時間效應(yīng)的時間函數(shù)模型來進行，該時間函數(shù)往往由實測數(shù)據(jù)擬合求得。受測量精度的限制，開采沉陷學(xué)中，規(guī)定的地表初始移動時間t=0時的地表下沉值為10 mm，以往學(xué)者大多以通過原點的函數(shù)作為動態(tài)沉陷預(yù)計函數(shù)，并基于實測成果反算求得各模型中的參數(shù)值，沉陷實測-理論模型構(gòu)建-實踐應(yīng)用的過程，往往由不對應(yīng)的時間函數(shù)帶來一定誤差。

大量的實測資料表明，Knothe,Logistic等時間函數(shù)與實際的地表沉陷的動態(tài)過程不完全符合，不能科學(xué)地描述地表動態(tài)下沉過程。筆者基于朔州地區(qū)特厚煤層高強度綜放開采的地表移動變形的特殊性，在分析常用時間函數(shù)地表點沉陷變形的特征基礎(chǔ)上，引入Von Bertalanffy時間函數(shù)，基于此函數(shù)改進該函數(shù)，使之成為三參量的時間函數(shù)，并應(yīng)用Matlab擬合求取各參數(shù)，最后基于該模型并結(jié)合概率密度影響函數(shù)對井工一礦高壓線塔下回采進行地表動態(tài)沉陷計算分析，并成功指導(dǎo)地面線塔的基礎(chǔ)加固與糾偏治理。

1 動態(tài)沉陷時間函數(shù)

在許多時間序列預(yù)測中常用“S型”曲線描述事物的增長或衰減過程，由于“S型”曲線形式多樣、結(jié)構(gòu)不同，對同樣觀測值預(yù)測結(jié)果略有差別。在時間序列里，有些變量的增長量最初比較小，隨時間的增加逐漸增長而達到一個快速增長時期，而后增長速度趨緩，最終達到穩(wěn)定的總增長量，這一過程若用曲線來表示，則是一種沿時間軸方向拉長的“S型”曲線。

1.1 常見函數(shù)模型

若用wm表示地表任意點沉陷穩(wěn)定后的最大下沉量(單位為mm)，w(t)表示任意時刻的動態(tài)下沉，v(t)表示任意時刻的地表下沉速度(單位為mm/d)，a(t)表示任意時刻的地表下沉加速度(單位為mm/d2)。

(1)Knothe函數(shù)。該模型假定地表下沉速率dW(t)/dt與地表最終下沉值wm和某一時刻t的動態(tài)下沉值w(t)之差成比例，推導(dǎo)并提出:

w(t)=wm(1-e-ct)

(1)

式中，c為時間影響參數(shù)，與采場覆巖巖性有關(guān)。

為進一步研究該模型在煤礦開采沉陷預(yù)計中的適用性，對該函數(shù)進行一階和二階函數(shù)求導(dǎo)，分析該下沉速度v(t)=wmce-ct和下沉加速度a(t)=-wmc2e-ct的函數(shù)，可知Knothe時間函數(shù)具有以下特征:① 地表下沉曲線形狀非“S型”;② 地表剛開始受到沉陷擾動，即t=0時下沉速度最大，然后逐漸減小至0;③ 下沉加速度恒為負值，下沉加速度一開始以最大加速度逐漸衰減;④ 僅有1個參數(shù)c，實際的動態(tài)沉陷表達性不足。

(2)Logistic模型。Logistic增長模型[13]在生態(tài)學(xué)、人口學(xué)等領(lǐng)域得到廣泛地應(yīng)用;它反映了事物發(fā)生、發(fā)展、成熟并趨于飽和(極限)的過程。其函數(shù)式為

(2)

式中，g,f為與采場上覆巖層巖性有關(guān)的參數(shù)。

(3)Gomepertz模型。由英國數(shù)學(xué)家B.Gomepertz提出并應(yīng)用于動植物的消亡規(guī)律，現(xiàn)多用于非采動影響下的路堤或地基的沉降，其函數(shù)式[15]為

w(t)=wme-e-ct+d

(3)

式中，c,d為與采場上覆巖層巖性有關(guān)的參數(shù)。

對該函數(shù)進行一階和二階函數(shù)求導(dǎo)，分析該下沉速度v(t)=wmce-ed-cted-ct和下沉加速度a(t)=-wmc2e-ct的函數(shù)，可知Gomepertz時間函數(shù)具有以下特征:① 下沉曲線是“S型”，與煤礦開采地表沉陷特征較相符;② 地表移動初始期，t=0時，地表下沉速度v≠0、下沉加速度a≠0，這與實際情況不符。

(4)Weibull模型。Weibull曲線函數(shù)在軟土地基的沉降預(yù)測中適用性較好，其函數(shù)式[11]為

w(t)=wm(1-e-ntk)

(4)

式中，n,k為與采場上覆巖層巖性有關(guān)的參數(shù)。

對該函數(shù)進行一階和二階函數(shù)求導(dǎo)，分析該下沉速度v(t)=wmnktk-1e-ntk和下沉加速度a(t)=wmnke-ntk[(k-1)tk-2-nkt2k-2]的函數(shù)，可知Weibull時間函數(shù)具有以下特征:① 該函數(shù)是對Knothe時間函數(shù)的改進;② 收斂速度較快，衰退期的下沉“拖尾”性質(zhì)表現(xiàn)不足;③ 時間函數(shù)中的參數(shù)物理意義尚不明確，難推廣。

(5)Richards時間函數(shù)。Richard生長方程[7]為一個3參數(shù)的非線性方程，主要用于生物生長的預(yù)測，其函數(shù)方程為

w(t)=wm(1-b1e-b2t)1/(1-b3)

(5)

式中，b1為初始沉降值參數(shù);b2為初始沉降速度參數(shù);b3為曲線形狀參數(shù)。

對該函數(shù)進行一階和二階函數(shù)求導(dǎo)，分析該下沉速度和下沉加速度的函數(shù)，可知Richard時間函數(shù)具有以下特征:① Richards函數(shù)可較好地表達煤礦開采后地表點任意時刻的下沉、下沉速度、下沉加速度等;② 參數(shù)對曲線形狀影響敏感性差，不同采礦地質(zhì)條件下曲線形狀變化規(guī)律性不明顯。

(6)雙曲線模型。雙曲線函數(shù)模型主要用于預(yù)測地基和路基的沉降，其函數(shù)模型在一定程度上可反映次固結(jié)的影響。

(6)

式中，b為與采場上覆巖層巖性有關(guān)的參數(shù)。

對該函數(shù)進行一階和二階函數(shù)求導(dǎo)，分析該下沉速度v(t)=b/(b+t/wm)2和下沉加速度a(t)=-2b/[(b+t/wm)3wm]的函數(shù)，可知雙曲線時間函數(shù)具有以下特征:① 地表動態(tài)下沉曲線形狀非S型;② 與Knothe等時間函數(shù)相似，t=0時，下沉速度最大，然后逐漸減小為0;③ 下沉加速度恒為負值，然后逐漸衰減為0。

(7)分段Knothe時間函數(shù)。

(7)

式中，c為與采場上覆巖層巖性有關(guān)的參數(shù);T,τ分別為地表總移動周期與地表下沉速度達到最大的時刻。

與Knothe函數(shù)相似，但其又有新的特點:① 人為對時間函數(shù)進行分段劃分，當(dāng)最大下沉速度時的τ時刻時，地表下沉值并非為1/(2wm)，與實際情況的相差較大;② 下沉的加速過程和減速過程是對稱的，沒有反映出沉陷中活躍期的快速性、衰退期的長期性和拖尾性。

(8)冪函數(shù)-Knothe時間函數(shù)。考慮到巖層與地表移動的非線性過程，在Knothe函數(shù)式中增加一個冪級系數(shù)k。

w(t)=wm(1-e-ct)k

(8)

式中，c,k為與采場上覆巖層巖性有關(guān)的參數(shù)。

分析下沉速度v(t)=wmcke-ct(1-e-ct)k-1和下沉加速度a(t)=wmc2k(2e-2ct-e-ct)(1-e-ct)k-2的函數(shù)，可知雙曲線時間函數(shù)具有以下特征:① 在Knothe函數(shù)的基礎(chǔ)之上，增加了冪級系數(shù)k，可以適當(dāng)?shù)母淖兿鲁燎€的形狀;② 在一定程度上克服了Knothe時間函數(shù)某些方面缺點，但缺少地表下沉速度對函數(shù)的控制影響研究。

(9)MMF時間函數(shù)模型。MMF模型是于1975年提出是一種“S型”生長曲線模型，王軍保等[1]將該模型用于地表沉降時間預(yù)測中。

(9)

式中，c,b為時間影響參數(shù)，與采場覆巖巖性有關(guān)。

對該函數(shù)進行一階和二階函數(shù)求導(dǎo)，分析該下沉速度v(t)=wmbctb-1/(c+tb)2和下沉加速度a(t)=wmbctb-2(ba-c-btb-tb)/(c+tb)3的函數(shù)，可知MMF函數(shù)具有以下特征:① 該函數(shù)可以較好地表示不同地質(zhì)采礦條件下的地表下沉;②b與c的物理意義有待研究;③ 時間t=0時下沉速度和加速度都無法表示。

(10)Sroka-Schober時間函數(shù)。Sroka以Knothe函數(shù)為基礎(chǔ)建立的時間函數(shù)為

(10)

式中，f為巖層的相對收斂率;ξ為上覆巖層的時間系數(shù)，用于描述巖層的滯后效應(yīng)。

該函數(shù)具有以下特征:① 該函數(shù)曲線非S型;② 下沉速度更接近于理想的時間函數(shù)，但下沉加速度曲線形態(tài)不正確;③ 實際應(yīng)用中，f和ξ不是已有參數(shù)，參數(shù)難以獲取，并且f和ξ難以同時求出。

大量的實測資料表明，Knothe,Logistic,Gomepertz與Weibull等常見時間函數(shù)與實際的地表沉陷的動態(tài)過程不完全符合，這些函數(shù)雖能在一定程度上對地表動態(tài)移動變形進行預(yù)計，但這些模型或不能反映地表下沉速度和加速度的變化規(guī)律，或函數(shù)中的參數(shù)物理意義不明確，不能反映各參數(shù)隨地質(zhì)采礦條件的變化。

1.2 理想函數(shù)模型特征

理想的時間函數(shù)模型應(yīng)該能夠反映:工作面回采尺寸逐步增大，巖層移動變形傳播至地表，地表在可監(jiān)測范圍內(nèi)隨時間變化應(yīng)為規(guī)律性地移動變形。依據(jù)煤礦開采后地表移動變形的動態(tài)變形實測情況，結(jié)合前述學(xué)者對地表動態(tài)沉陷時間函數(shù)的研究成果，理想的時間函數(shù)模型應(yīng)具有以下特點:①t=0時，地表某點的下沉量w(t)=0或為可監(jiān)測值(受監(jiān)測精度限制，一般取值為10 mm)、下沉速度v(t)=0、下沉加速度a(t)=0;② 在初始期的下沉量、下沉速度與下沉加速度較小、為正值且逐漸增大;③ 活躍期下沉速度與下沉加速度較大，下沉加速度變化過程為先增大至峰值后、逐漸減小至0、然后負值、負值達到峰值后逐漸減小;④ 衰退期下沉速度為正、逐漸變小且逐漸趨于0，下沉加速度由負值逐漸增大且逐漸趨于0;⑤ 整個下沉-時間曲線近似“S”型曲線，由于巖塊與裂隙的壓實、蠕變歷時較長，衰退期長而平緩;⑥ 針對不同地質(zhì)采礦條件，可通過形狀參數(shù)調(diào)整動態(tài)下沉曲線形狀的陡峭或平緩、可通過調(diào)整沉降速度參數(shù)調(diào)整地表動態(tài)下沉速度。

2 基于Bertalanffy的時間函數(shù)模型

2.1 模型及其改進

Von Bertalanffy模型[17-18]為常見的生物生長模型，主要用于預(yù)測動植物的生長規(guī)律，也常應(yīng)用于經(jīng)濟與統(tǒng)計學(xué)中。該模型是Logistic的推廣，且具有較好的靈活性，其函數(shù)方程為

w(t)=wm(1-be-ct)3

(11)

式中，b與c為時間影響參數(shù)，其中b可反映t=0時的初始沉降值，c反映t=0時的初始沉降加速度值。

模型改進(三參模型):考慮到巖層與地表移動過程的非線性，尤其對于近淺埋深特厚煤層綜放開采，地表沉陷變形集中、盆地邊緣陡峭等特征，將Bertalanffy函數(shù)式中的3次冪變?yōu)闀r間影響參數(shù)d(反映曲線的形狀，可反映不同地質(zhì)采礦條件下的曲線陡峭、收斂速度、拖尾性等)，用于煤礦開采地表動態(tài)沉陷研究中，可將Bertalanffy函數(shù)改進為三參數(shù)模型:

w(t)=wm(1-be-ct)d

(12)

根據(jù)模型中參數(shù)對曲線形態(tài)或量值的控制性作用，式中的參數(shù)b為初始沉降值參數(shù)(可反映地表初始沉降值及地表移動初始期的長短)，c為初始沉降速度參數(shù)，d為曲線形狀參數(shù)。

2.2 模型特點分析

2.2.1不通過原點性

煤礦開采地表沉陷動態(tài)預(yù)計的實際應(yīng)用中，在選擇好模型情況后，將實測參數(shù)代入函數(shù)模型，反算求得各參數(shù)。在實際“三下”采煤或沉陷保護應(yīng)用中，受測量精度的限制，開采沉陷預(yù)計參數(shù)求取及沉陷預(yù)計所涉及的邊界角、超前影響角等都是基于地表下沉量值為10 mm點作為計算依據(jù)的。即多個角量參數(shù)通常以地表下沉量達到10 mm為地表移動起算點，t=0時，地表下沉取值w(0)=10 mm。

Knothe時間函數(shù)的收斂速度隨地質(zhì)采礦條件的調(diào)整性較差，僅通過一個參數(shù)c很難進行調(diào)整;在e-ct項之前缺少一個沉降初始值的控制參數(shù)，改進Bertalanffy三參數(shù)函數(shù)將這個缺陷進行了補充。

2.2.2“S型”曲線特征

為進一步研究該模型在煤礦開采沉陷預(yù)計中的適用性，對該函數(shù)進行一階和二階函數(shù)求導(dǎo)(各函數(shù)的時間函數(shù)曲線如圖1所示)。

v(t)=wmbcd(1-be-ct)d-1e-ct

(13)

a(t)=wmbc2de-ct(1-be-ct)d-2(bde-ct-1)

(14)

圖1 改進Bertalanffy三參數(shù)模型動態(tài)下沉、下沉速度及下沉加速度時間函數(shù)Fig.1 Dynamic subsidence,velocity and acceleration time function of improved three-parameters model

結(jié)合圖1所示的形態(tài)曲線，分析改進Bertalanffy三參數(shù)時間函數(shù)的動態(tài)下沉、下沉速度和下沉加速度時間函數(shù)，具有以下特征:

(1)單調(diào)性。因參數(shù)b,c與d均為正值，(1-be-ct)d>0，所以w(t)是t的增函數(shù)。

(2)凹凸性。令w″(t)=a(t)=0，求得拐點的時間坐標(biāo)t=ln(bd)/c，此時w(t0)=wm(1-1/d)d，并且當(dāng)tt0時，下沉曲線是“凸”的;即w(t)為時間t的“S型”曲線。

(3)漸進性。為有界函數(shù)，漸近線為w(t)=wm。

2.3 模型參數(shù)討論

依據(jù)多個煤礦地表移動變形實測數(shù)據(jù)表明，改進Bertalanffy三參模型參數(shù)與地質(zhì)采礦條件關(guān)系密切，不同的采深/采厚比值、采動充分程度、控制沉陷的關(guān)鍵層破斷位置高度、工作面推進速度等因素，地表移動變形越劇烈程度、地表移動的初始期與活躍期的時間跨度、地表點的動態(tài)下沉速度、動態(tài)下沉加速度有一定的差異。Bertalanffy改進三參模型中的參數(shù)取值亦可通過相應(yīng)調(diào)整適應(yīng)不同地質(zhì)采礦條件下的地表動態(tài)變形計算。

2.3.1初始沉降參數(shù)b

改進Bertalanffy三參模型中參數(shù)b反映t=0時的初始沉降值與地表移動初始期的長短。b的經(jīng)驗取值范圍(0,1]，當(dāng)b=1時，函數(shù)可通過原點;相同地質(zhì)采礦條件下，工作面推進速度越快，地表移動越劇烈，b越大。該參數(shù)的變化規(guī)律如圖2所示(圖中各曲線的其他變量c=0.032 5，d=36)。

圖2 初始沉降參數(shù)b對模型的影響規(guī)律Fig.2 Influence law of initial settlement parameter b

2.3.2動態(tài)沉降加速度參數(shù)

c反映t=0時的初始沉降加速度值。c的經(jīng)驗取值范圍[0.005,0.500];地表及覆巖移動變形越劇烈(工作面推進速度越快等因素)，地表動態(tài)沉降加速度越大，曲線越陡，c值越大;反之，地表移動劇烈程度越低，c值越小。該參數(shù)的變化規(guī)律如圖3所示(其他變量b=0.992 0，d=36)。

圖3 加速度參數(shù)c對模型的影響規(guī)律Fig.3 Influence law of acceleration parameter c

2.3.3形狀系數(shù)d

d反映S型曲線的形狀，可反映不同地質(zhì)采礦條件下的地表動態(tài)下沉初始期的長短、沉陷收斂速度與拖尾性等。d的經(jīng)驗取值范圍為[1,100];曲線形狀控制參數(shù)d越大，地表初始期越長，受采動影響反映至地表的時間越長;同時d越大，地表收斂速度越小，地表移動延續(xù)時間越長;當(dāng)d=1時(同時b=1)，即為常見的Knothe時間函數(shù)。該參數(shù)的變化規(guī)律如圖4所示(其他變量b=0.992 0，c=0.032 5)。

圖4 形狀系數(shù)d對模型的影響規(guī)律Fig.4 Influence law of curve shape factor d

3 參數(shù)擬合求取與動態(tài)沉陷應(yīng)用驗證

3.1 研究區(qū)域地質(zhì)采礦及觀測站概況

研究區(qū)域東坡煤礦位于山西朔州市，井田內(nèi)覆蓋有40～60 m厚的黃土，914走向長壁工作面開采9號煤，工作面傾向長約240 m，連續(xù)推進長度約1 250 m，開切眼處平均采深265 m，平均采出煤層厚度為14.4 m;采煤方法為特厚煤層綜放一次采全高，平均煤層傾角2.3°;工作面推進速度約為2.8 m/d，開采強度較大[19-20]。

914工作面A測線為走向半盆地觀測線，總長約為480 m，布置測點24個;自2011年11月至2014年5觀測15次，歷時約30個月。

3.2 時間函數(shù)求參及擬合優(yōu)度分析

擬合優(yōu)度是指回歸曲線對觀測值的擬合程度;度量擬合優(yōu)度的統(tǒng)計量是可決系數(shù)(亦稱確定系數(shù))R2;R2取值范圍為(0,1)，R2越接近1，說明回歸曲線對觀測值的擬合程度越好，反之R2值越小，擬合程度越差。對于元素個數(shù)為n的序列Yi，擬合值為yi，則擬合優(yōu)度R2可表示為

(15)

根據(jù)東坡煤礦A測線25號測點地表移動觀測站的實測數(shù)據(jù)，應(yīng)用Matlab對改進Bertalanffy三參模型時間函數(shù)進行曲線擬合求參，同時與Knothe時間函數(shù)作對比。

3.2.1Knothe時間函數(shù)

Knothe時間函數(shù)是我國煤礦動態(tài)沉陷預(yù)計中應(yīng)用最早也是最為常見的時間函數(shù)，東坡煤礦A測線25號測點實測數(shù)據(jù)擬合求取的Knothe時間函數(shù)(圖5)為

w(t)=wm(1-e-0.005 7t)

(16)

圖5 25號測點不同時刻地表下沉曲線Fig.5 Subsidence curves at different times of No.25 points

Knothe時間函數(shù)中僅有一個參數(shù)c，不能有效地對“S型”動態(tài)下沉曲線形態(tài)、下沉加速度的曲線陡峭程度及初始沉降數(shù)值進行調(diào)整。25號測點應(yīng)用Knothe時間函數(shù)的擬合優(yōu)度R2為0.817。

3.2.2改進Bertalanffy三參模型時間函數(shù)

現(xiàn)實實測中時間t=0對應(yīng)下沉等于0 mm的時刻受監(jiān)測精度的限制，難以捕捉，實際應(yīng)用中若仍基于過原點的時間函數(shù)進行擬合求參，再將該函數(shù)應(yīng)用于動態(tài)沉陷預(yù)計工作中，這一過程將對時間函數(shù)有一定的誤差影響;改進Bertalanffy三參模型時間函數(shù)可不經(jīng)過原點(初始沉降值參數(shù)b=1時，可經(jīng)過原點)，并將(t0=0,w0=0)=(0,10)作為時間函數(shù)的初始值進行擬合求參(圖5)，25號測點應(yīng)用改進Bertalanffy三參模型時間函數(shù)擬合優(yōu)度R2為0.984。擬合求取的改進Bertalanffy三參模型時間函數(shù)為

w(t)=wm(1-0.992e-0.032 5t)36

(17)

由圖5可直觀地看出改進Bertalanffy三參模型時間函數(shù)模型與實測值很接近;應(yīng)用Knothe時間函數(shù)的擬合優(yōu)度R2為0.817，應(yīng)用改進Bertalanffy三參模型時間函數(shù)擬合優(yōu)度R2為0.982;后者時間函數(shù)相對誤差較小，同時該模型中3個參數(shù)均被賦予了不同的物理意義，各個參數(shù)變化均與地質(zhì)采礦條件有一定程度的聯(lián)系。該時間函數(shù)模型方便簡潔、符合理想時間函數(shù)模型的要求，可較好地模擬采動區(qū)地表動態(tài)沉降過程，能夠計算出地表移動持續(xù)時間和地表點在某一時刻的下沉值、下沉速度、加速度等動態(tài)參數(shù)，且模型的可塑性較強，具有較廣泛的適用性。

3.3 高壓線塔受動態(tài)沉陷影響現(xiàn)場驗證

隨著高壓輸電網(wǎng)的規(guī)?；ㄔO(shè)，高壓線塔下伏煤炭開采及對地面線塔的動態(tài)變形影響與抗采動治理工作越來越受到重視。中煤平朔井工一礦試驗工作面采煤工藝為特厚煤層綜放開采，開采4號煤，機采高度為3.2～3.5 m，平均采出煤厚為8.7 m;工作面中部對應(yīng)地表有1個基塔，線塔附近工作面采深為200～260 m，開切眼側(cè)工作面平均推進速度為4.2 m/d。井工一礦與東坡煤礦工作面最近平面距離約2 km。

對于高壓線下采煤地面控制變形一般可分為井下特殊采煤控制與地面線塔加固與糾偏措施。為保證高壓線壓覆煤炭的安全回采，首先基于地表動態(tài)沉陷變形預(yù)計，對線塔周邊一定范圍(沿工作面推進方向的開切眼側(cè)60 m、沿工作面推進方向的終采線側(cè)30 m、工作面回采巷道兩側(cè)各50 m)采取只采不放的措施。地面線塔為4個獨立體基礎(chǔ)，平行工作面推進方向(垂直于線路延伸方向)的兩個基礎(chǔ)間距為10 m，沿工作面推進方向的兩個基礎(chǔ)間距為8 m;為提高線塔的抗采動能力，將線塔原獨立基礎(chǔ)加固改造為鋼筋混凝土連體的“井字形”基礎(chǔ)，同時實現(xiàn)工作面開采過程中的糾斜。

參考巖性對地表沉陷的影響規(guī)律[21]，結(jié)合2個礦井差異，依據(jù)改進Bertalanffy三參模型時間函數(shù)各參數(shù)隨地質(zhì)采礦條件的差異規(guī)律(4號煤層位于9號煤層的上部，比9號煤層工作面地表移動初始期短、盆地邊緣相對平緩;同時4號煤試驗工作面采厚約8.7 m，比東坡煤礦914工作面平均采厚為14.4 m的數(shù)值小)，對各參數(shù)進行適當(dāng)調(diào)整，井工一礦試驗工作面選用的動態(tài)下沉?xí)r間函數(shù)為

w(t)=wm(1-0.994e-0.152t)32

(18)

考慮到部分區(qū)域的只采不放范圍，結(jié)合工作面推進速度，對試驗工作面沿水平方向?qū)⒉沙雒簩雍穸?.7 m，劃分為機采高度部分3.2 m和放煤高度部分5.5 m，地表下沉系數(shù)取值q=0.83;根據(jù)實測資料當(dāng)工作面超前線塔約96 m時，線塔地表首次擾動，t=0;當(dāng)工作面超前線塔為60 m時，t=8.6 d;當(dāng)工作面推過線塔30 m時，t=21.5 d，因此井工一礦試驗工作面的動態(tài)下沉?xí)r間函數(shù)分段表述為

(19)

其中，f1(x)與f2(x)為線塔開切眼側(cè)和終采線側(cè)距離線塔為x時設(shè)計規(guī)劃的放煤部分對線塔處的影響函數(shù)，如圖6所示。

圖6 線塔兩側(cè)放煤區(qū)對線塔的影響函數(shù)Fig.6 Influence function of top coal caving area on tower

線塔兩側(cè)放煤部分距離線塔越近，對線塔的采動影響越大，反之越小;其中f1(x)與f2(x)可由含主要影響半徑r的函數(shù)表示為

(20)

借鑒東坡煤礦實測資料及巖移數(shù)據(jù)成果取tanβ=2.6，則r=96 m，代入式(20)并經(jīng)積分運算可求得f1(x)=0.059，f2(x)=0.217。應(yīng)用式(19)所示的分段函數(shù)繪制的地表動態(tài)下沉曲線與線塔基礎(chǔ)的實測值對比情況如圖7所示。

圖7 時間函數(shù)+影響函數(shù)法的變采高地表動態(tài)沉陷曲線Fig.7 Dynamic surface subsidence curve of variable mining height based on time function and influence function method

由于式(20)僅考慮了沿工作面推進方向的各影響函數(shù)疊加，無法對只采不放段x=(-60,30 m)沿工作面布置方向放煤段對影響函數(shù)的影響數(shù)值;同時特厚煤層綜采放頂煤開采工藝復(fù)雜性，采放煤層厚度不容易控制等;圖7繪制的基于影響函數(shù)的變采高地表動態(tài)沉陷曲線與實測數(shù)據(jù)相比雖有一定誤差，但在實踐生產(chǎn)中，應(yīng)用該思路安全有效地采出了線塔壓覆煤炭資源，并應(yīng)用基于影響函數(shù)結(jié)合改進Bertalanffy三參模型時間函數(shù)指導(dǎo)線塔的動態(tài)傾斜糾偏，為該類工作提供借鑒意義。

4 結(jié) 論

(1)分析了常見Knothe等時間函數(shù)各時間段的地表下沉曲線形態(tài)、下沉速度與下沉加速度，對各函數(shù)優(yōu)缺點及適用性進行了分析。

(2)引入Bertalanffy時間函數(shù)并對其進行改進，分析了改進Bertalanffy三參時間函數(shù)中各參數(shù)的物理意義及隨地質(zhì)采礦條件的聯(lián)系。該函數(shù)與實測值符合性較好，擬合優(yōu)度較高;該函數(shù)模型的可塑性較強，具有較廣泛的適用性。

(3)針對特厚煤層綜放開采的地面線塔等建(構(gòu))筑物的動態(tài)沉陷預(yù)計與動態(tài)糾偏，提出基于改進Bertalanffy三參時間函數(shù)結(jié)合影響函數(shù)法，以解決只采不放等區(qū)域性變采高條件下地表動態(tài)沉陷預(yù)計問題。