車 平 刁志翔 王奇石 商翔宇

(1.江蘇省有色金屬華東地質勘查局,江蘇 南京 210007;2.中國礦業(yè)大學 力學與土木工程學院,江蘇 徐州 221116;3.中國礦業(yè)大學深部巖土與地下工程國家重點實驗室,江蘇 徐州 221116)

華東礦區(qū)是我國非常重要的煤炭能源基地。自1987年以來,該礦區(qū)相繼發(fā)生井筒破裂100余處,嚴重影響了礦區(qū)的安全生產(chǎn)。為此,學者們圍繞立井井筒的破壞機理開展了大量研究。目前,附加應力破壞理論已經(jīng)得到廣泛認可[1]:即底部含水層疏水引發(fā)周圍土層固結和沉降,沉降地層在立井井筒上施加豎直向下的附加應力,附加應力與井筒原有的應力疊加,導致井筒破裂。在此基礎上,學者們開發(fā)出多種有效技術措施來防治井壁破裂[2,3]。然而,經(jīng)修復加固后的破裂井壁僅能在不確定的一段時間內安全運行,并且?guī)缀醪豢杀苊獾貢俅伟l(fā)生破裂[4]。歸根結底,這是由于井筒周圍地層長期持續(xù)沉降造成的。

已有大量研究表明,底部含水層疏水是造成立井井筒周圍深厚地層沉降的主要原因[1,5,6]。也有關于立井井筒周圍不同深度地層變形對地表長期沉降貢獻的研究。李文平和于雙忠分析認為,深埋黏土天然含水量小于吸附結合水含量,因而可能不具有壓縮性[7];但李定龍分析認為,底部含水層上、下的微弱隔水黏土層中部分弱結合水會隨底部含水層和下覆裂隙基巖水位的下降而逐漸釋放[6]。王秀艷等試驗研究指出,深層黏土中以結合水為主,但仍能釋水壓縮[8];許延春等[5]的實測結果指出,底部含水層深部砂土及相鄰深埋黏土的壓縮變形是地表沉降的主要來源。柴敬和楊金宏等[9,10]基于光纖技術的實測數(shù)據(jù)也表明,厚度較大的深部砂土和黏土層均壓縮明顯。由以上文獻分析可知,深部黏土地層的壓縮變形很可能對井筒周圍地表長期沉降有貢獻,但具體貢獻占比如何? 伴隨著井壁的局部破裂,其在井筒整個服役周期內將如何變化? 這些問題的正確解答對于井筒治理工作具有重要意義[4],但目前尚未有相關研究報道。

筆者擬利用模擬方法,對華東礦區(qū)孔莊煤礦主井井筒周圍地表長期沉降進行模擬對比研究,以期初步回答前述問題,所獲研究結論對破裂立井井筒的長期安全管理有借鑒意義。

1 立井井筒周圍地表長期沉降數(shù)值模型

1.1 孔莊主井基本地質條件

孔莊煤礦主井于1975年4月投入生產(chǎn),井筒深度465.8 m,井筒凈直徑5.00 m。根據(jù)實際統(tǒng)計資料[14],貫穿主井井筒的第四系表土層厚度約為156 m,主要由黏土、砂質黏土及砂層組成。根據(jù)實際情況及研究需要,為簡化計算,將相似相近地層進行合并,從而將實際地層簡化為7個含水層及6 個弱透水層,其組成、厚度及埋深如表1所示。

表1 簡化后的地層分布

1.2 數(shù)值模型采用的簡化假設

由于該井筒建設年代較遠,缺少必要的土層力學和水力試驗結果,本文假設:①各弱透水層和含水層分別假設為同類均質黏土和砂土,除了先期固結壓力和初始孔隙比不同外,其他力學參數(shù)均相同,滲透系數(shù)由孔隙比決定;②假定含水層及弱透水層均為正常固結狀態(tài)的飽和土體;③假設模擬土層的力學性質和滲流性質均為各向同性,土中水的運動符合達西定律;④忽略地表建筑對地面沉降的影響;⑤根據(jù)該煤礦水文地質資料,僅考慮地層靜水壓力;⑥主要關注井筒周圍地層的長期沉降變形,不考慮大氣降水對淺部含水層孔隙水壓的波動影響[1]。

1.3 數(shù)值模型的建立

計算基于商用有限元軟件ABAQUS 進行,模型按常用的比奧固結理論,考慮滲流固結耦合作用;采用空間軸對稱模型,選用二階精度的減縮積分位移-孔壓耦合單元;模型高度取為156 m(與實際表土層厚度相同),模型寬度取為50 m(根據(jù)最遠沉降觀測點距井筒的距離),因模型網(wǎng)格劃分相對規(guī)則,故不再給出圖示。參考該煤礦設計生產(chǎn)年限,模型總模擬時間為100 a。

需要確定表1所示土層的力學參數(shù)以及滲流參數(shù)。土體采用修正劍橋模型,其中黏土力學參數(shù)見表2,系由前期同地區(qū)深部黏土的實驗結果確定[15]。其中,λ是臨界狀態(tài)線在v-lnp′平面(v為比容,p′為有效平均應力)上的斜率,κ是卸載曲線在v-lnp′平面上的斜率,M是臨界狀態(tài)線在p′-q平面(q為等效偏應力)上的斜率,e1是正常固結曲線在p′=1 kPa時的孔隙比,β是控制p′-q平面上屈服面形狀的參數(shù),K是控制偏平面上屈服面形狀的參數(shù),ν是泊松比。

表2 黏土力學參數(shù)

模型的砂土力學參數(shù)則參考文獻中典型標準砂的臨界土力學參數(shù)[16],具體見表3。

表3 砂土力學參數(shù)

由于不同深度土層的滲透系數(shù)不同,因此,考慮滲透系數(shù)隨孔隙比的非線性變化,具體變化規(guī)律采用作者前期的實驗結果[17]:

SPME獲得的樣品于安捷倫GC*GC-Q-TOF(7890B/7200型)進行分析。調制周期：3 s；熱噴時間：350 ms。一維色譜柱為HP-5ms，30 m×0.25 mm×0.25 μm；二維色譜柱為DB-17HT，1.8 m×0.1 mm×0.1 μm。溶劑延遲時間：3 min，不分流進樣，載氣為He，流速：1 mL/min。進樣口溫度為250 ℃。升溫程序為40 ℃保持2 min，5 ℃/min升至270 ℃，保持10 min，共58 min。質譜為EI源，接口溫度：280 ℃，掃描范圍：40～700 m/z。

黏土的滲透系數(shù)由公式(1)計算,而砂土則參照經(jīng)驗,取有相同孔隙比黏土滲透系數(shù)的3倍。

模型的定界條件包括初始條件和邊界條件。其中模型初始條件包括各地層的初始豎向有效應力、初始孔隙比、初始孔壓和初始干密度等,根據(jù)工程實際地質情況以及前期深部土壓縮試驗結果[17]確定,具體見表4。

表4 各地層的初始條件

模型邊界條件包括位移邊界條件以及水力邊界條件,其中,位移邊界條件較為常規(guī),模型上邊界豎向可自由移動,而水平方向固定;模型底部設定為剛性基座,水平與垂直方向均不可移動;模型左、右邊界垂直方向可移動,水平方向固定。水力邊界條件對地表長期沉降模擬相當關鍵,因此,下文著重介紹模型水力邊界的考慮與確定。

模型上部水力邊界較易確定,設為排水邊界,而底部水力邊界則需要根據(jù)工程實際區(qū)分不同情況。首先,工程實測資料表明,該煤礦底部含水層水位多年來一直持續(xù)下降,很可能是造成地層沉降的主要因素,因此,僅將底部設為孔壓降低邊界條件,記為邊界條件1。其次,該煤礦主井井筒內實際淋水觀測表明,靠近底部含水層的上臨含水層很可能沿破裂井筒外壁與底部含水層形成水力聯(lián)系,李定龍也曾分析指出,地層水可能在底部含水層附近產(chǎn)生向井筒微裂隙的滲流運動[6]。因此,本模型將最靠近底部的含水層⑥也設置為疏水地層,其水位下降量為底部含水層的1/4,將底部含水層和含水層⑥孔隙水壓力同時下降,記為邊界條件2?？梢?上述底部水力邊界條件1對應于“單純底含失水”疏水條件,而上述水力邊界條件2則對應于“底含失水+井筒破裂失水”的疏水條件。

目前,孔莊煤礦主井僅有近10 a的附近水文觀測數(shù)據(jù)(見圖1)。數(shù)據(jù)表明,2008年到2018年底含水層水位下降8 m,即其水頭降低80 k Pa;此外,實際運營結果顯示,隨著井筒服役時間的不斷增長,井壁損壞程度日益加大,井筒的涌水量不斷增加,水位下降速度越來越快,有理由推斷近10 a是水位下降最多的10 a。據(jù)此分析,水力邊界條件1可由圖2所示的底含上邊界孔壓變化曲線表示,即前40 a的水位每10 a分別下降1 m、2 m、4 m、8 m,相應底部含水層水頭降低分別為:10 k Pa、20 kPa、40 k Pa和80 k Pa。水力邊界條件2的具體數(shù)值可以根據(jù)前述假設由邊界條件1確定。

圖1 底含水位變化

圖2 底含上邊界孔壓變化

2 模擬結果與實測數(shù)據(jù)的對比與分析

通過模型前40 a的模擬結果與已有實測數(shù)據(jù)的對比,以檢驗模型的合理性。

2.1 模擬結果與實測數(shù)據(jù)的對比

自2005年以來,礦方在孔莊主井井塔四周布置有4個沉降測點(即主1～主4),并在距井筒約30 m 處布置沉降測點N6,由于觀測期間部分測點被破壞,目前,只有主2和N6測點近10 a的數(shù)據(jù)最為完整。由于主井井塔的基礎和井筒本身相互分離,即井塔的沉降和井筒沉降無關,故監(jiān)測數(shù)據(jù)可反映地層沉降。

模擬含水層疏水的第30～40 a可以大致對應為2008～2018年(目前僅有該時間段的數(shù)據(jù)),為了與監(jiān)測數(shù)據(jù)進行對比,將井筒周圍主2、N6測點的實際監(jiān)測值和地層沉降模擬值分別繪于圖3中。

圖3 模擬和實測地表沉降對比

由圖3可以看到,主2和N6兩點的實際監(jiān)測值盡管不完全重合,但其趨勢基本相同,并且過去十年的總沉降量也基本相同。這意味著,實測地表沉降有相當比例源自于深部土層壓縮變形的貢獻,進而使得井筒-地層相互作用對上述兩個測點的實際沉降變形影響差異較小,因此,兩者沉降規(guī)律大致相當。

對比模擬結果與實際監(jiān)測值可以看到,前述兩類水力邊界條件所對應的模擬沉降規(guī)律和實際沉降規(guī)律均較為接近,模擬結果與實際沉降的最大相對誤差為-16.7%～+28.3%?？紤]到前述所采用的模型參數(shù)均由類比法選取,而且水力邊界條件是由實測水位下降推出,不可避免地會造成模擬結果和實際值存在偏差,但兩者總體規(guī)律較為接近,因此,可以認為模型基本合理。由圖3可以看到,近10 a內,實際監(jiān)測沉降量約為60 mm,邊界條件1的模擬沉降量約為50 mm,而邊界條件2模擬沉降量約為87 mm,由此可見,考慮井筒破裂新增失水通道會明顯增大總地表沉降量。

含水層疏水不僅會引起其自身的壓縮,而且還導致相鄰[5,9,13]甚至更遠含水層和弱透水層的固結壓縮。為了分析各土層特別是黏土層對地表沉降的貢獻量,對比了這兩種水力邊界條件下前40 a各層的具體貢獻值(如圖4所示),圖4中橫軸h代表“含水層”,r代表“弱透水層”?？梢?無論何種水力邊界條件情況下,深部弱透水層即深部黏土地層對地表沉降的貢獻均較大,具體而言,在邊界條件1下,深部黏土地層對地表沉降的貢獻值約為80%,而在邊界條件2下,深部黏土地層對地表沉降的貢獻值約為70%,而含水層的貢獻值較低,其原因在于,含水層是由砂土組成,其埋深較深,砂土較為密實,砂土的孔隙總體較小,因此,沉降量相對較小?？偠灾?砂土與黏土交替出現(xiàn)的深厚地層中,總的地表沉降中深部黏土層的貢獻相當可觀。

圖4 過去40 a各地層沉降貢獻值

此外,可以看到,如果只有底部含水層疏水,緊臨底部含水層的弱透水層⑥會成為對地層沉降貢獻最大的地層,并且與其他地層有著很大的差異,導致這種差異的主要原因是:由于弱透水層⑥的埋深相對較大,有效應力大、孔隙比小,相應的滲透系數(shù)非常低,孔隙水滲流速度極其緩慢,短期內不會影響上覆含水層和弱透水層中孔隙水的運動,更不會在短時間內進一步造成更遠的上覆地層的大量沉降;而如果含水層⑥也疏水,即為邊界條件2,則緊臨含水層⑥的弱透水層⑤和⑥都將成為對地表沉降貢獻最大的兩個地層,并且與其他地層之間的差異也有逐漸縮小的趨勢;可以預測,當井筒中存在多處漏水或下部多個含水層之間存在水力聯(lián)系時,地層沉降速率將更快,并且每個地層的貢獻將相對更加均勻。

2.2 服役周期內井筒未來地表沉降預測分析

上述模型前40 a的模擬結果與已有實測數(shù)據(jù)的對比表明了模型的合理性。以此為基礎,模型后60 a的計算結果可以作為預測該井筒周圍地層長期沉降規(guī)律的合理依據(jù)。

考慮一種理想的極限情況:第40 a之后,地層經(jīng)過加固、井筒經(jīng)過治理,底部含水層原有的以及井筒破裂新增的疏水通道均被關閉,不再疏水,即上述相應邊界孔隙水壓不再變化。圖2 中40～100 a的水平線段給出了相應的底部含水層孔壓。

模擬獲得的后60 a及前40 a地層的平均沉降速率見表5,各個地層沉降貢獻值見圖5,總體沉降曲線見圖6?？梢钥吹?盡管含水層不再疏水,但其內部孔隙水壓在重新平衡的過程中,其上覆地層的壓縮潛力逐漸釋放,沉降量持續(xù)增加。由于砂土滲透系數(shù)相對較大、變形較易穩(wěn)定,因此,每個含水層的后期貢獻值均較小,甚至之前疏水的含水層已經(jīng)停止下沉,而黏土層由于滲透系數(shù)小、壓縮性大,隨著時間的增加,其貢獻值逐漸增加。具體而言,在未來60 a內,弱透水層對地層沉降的總貢獻率將達到90%。由表5 和圖6可知,在水力邊界條件1情況下,井筒周圍地表的沉降速率在整個井筒工程服役期限內變化幅度很小;在水力邊界條件2情況下,后60 a的沉降速率較前40 a的降低也不足30%,這歸因于圖5中所示的黏土地層對井筒周圍地表沉降的貢獻具有長期性。由此可見,含深厚黏土地層的長期沉降性質決定了井筒破裂治理也將是個長期、多次的過程[4]。

圖5 未來60 a各地層沉降貢獻值

圖6 0～100 a模擬沉降曲線

表5 地層平均沉降速率

除上述計算模型外,筆者采用類似方法對徐莊煤礦主井周圍地層的長期沉降變形也進行了模擬計算,所獲結論與本文結果相近[18]。

3 結 論

筆者對孔莊煤礦主井周圍地表長期沉降進行了數(shù)值模擬研究,圍繞深立井周圍地表沉降中深厚黏土層的貢獻、井筒破裂新增失水通道對地表沉降的貢獻等問題,開展了數(shù)值模擬,并與實測結果對比分析,獲得了如下結論:

對孔莊煤礦主井的模擬結果表明,2008～2018年的井筒周圍地層沉降模擬結果和實測結果趨勢一致,量值相當,表明模型較為合理;模擬結果顯示:深厚黏土地層在立井井筒周圍地表沉降中貢獻顯著,前40 a內,與疏水含水層相鄰的深部黏土層壓縮變形顯著,其壓縮量占地層總沉降量的70～80%;即使是在考慮含水層不再疏水的理想極限情況下,后60 a內,地層沉降速率依然可觀,且與前40 a內疏水含水層相距較遠的弱透水層也產(chǎn)生了壓縮變形,其中,黏性弱透水層的壓縮變形對立井周圍地表沉降的總貢獻率達到了90%。此外,由于井筒破裂造成的新失水通道會加速總的地表沉降,對未來長期沉降有顯著影響。