王宏偉,王 晴,石瑞明,姜耀東,2,田 政
(1.中國礦業大學(北京) 力學與建筑工程學院,北京 100083;2.中國礦業大學(北京) 煤炭資源與安全開采國家重點實驗室,北京 100083;3.中國礦業大學(北京) 深部巖土力學與地下工程國家重點實驗室,北京 100083)
沖擊地壓作為典型的動力災害,由于其頻繁發生的危害性、極短前兆的突發性、機理研究和預測防治的復雜性等特點,給煤炭資源安全高效開采帶來巨大的挑戰,同時也是國內外學者高度關注的熱點課題。
在眾多誘發沖擊地壓的因素中,斷層等復雜地質構造賦存環境是沖擊地壓等動力災害的重要誘因。斷層失穩誘發沖擊地壓是由于斷層構造等不連續面的存在,在開采擾動下出現突然的相對滑動和能量的瞬間釋放,從而導致礦井劇烈振動和破壞的動力現象的發生。與其他典型的沖擊地壓類型相比,斷層失穩誘發的沖擊地壓所影響的范圍更加廣泛,災害發生時所釋放的能量強度更大、震級更高。
近年來,沖擊地壓等典型動力災害事故在煤礦開采中頻繁發生。例如,2018-10-20山東能源集團龍鄆煤業有限公司發生嚴重的沖擊地壓事故,由于受采掘、施工擾動和斷層滑動等影響導致地應力更加集中從而誘發動力災害,造成數名人員傷亡和巨大經濟損失;2019-06-09吉林煤業集團龍家堡礦業有限責任公司的沖擊地壓事故,該事故的直接原因是高水平構造應力、采掘活動影響和斷層活化等多因素共同導致,造成巷道嚴重破壞和人員傷亡;2019-08-02開灤(集團)有限責任公司唐山礦業分公司井田內斷層和褶皺構造附近發生了井下沖擊地壓事故,給煤礦的安全生產帶來了巨大損失;此外,2020-02-22山東新巨龍能源有限責任公司在開采工作面與斷層構造形成的強沖擊危險區域發生了較大沖擊地壓事故等。事故分析報告顯示,上述沖擊地壓事故的發生多與斷層等地質構造賦存狀態影響相關,事故高發區域地應力集中明顯。因此,研究斷層結構失穩誘發沖擊地壓的特征和機理對沖擊地壓的預測和防治至關重要。
關于斷層構造和沖擊地壓的相關性,國內外學者從各個角度研究了沖擊地壓機理,并提出一系列經典理論,為斷層構造失穩誘沖機理研究奠定了理論基礎。MICHALSHI最早提出了在臨近斷層開采時更容易誘發沖擊地壓的觀點,OHNAKA給出了斷層滑移動態摩擦定律假說;RICE等則認為沖擊地壓的發生主要是因為斷層面的摩擦和剪切作用,并分析了斷層滑移過程中溫度場變化特征;于廣明,謝和平等將分形理論應用于斷層活化研究,分析了斷層面分型界面效應與巖體破壞的關系;考慮斷層構造的影響,王學濱、林遠東等應用應變梯度塑性理論得到了斷層帶與圍巖失穩的判據及解釋,并推導了斷層滑移錯動位移理論公式;宋振騏在“實用礦山壓力控制”理論的基礎上,建立了過斷層開采的控制力學模型以及相關控制準則;紀洪廣等將庫侖應力的概念引入開采擾動下斷層滑移失穩分析中;HOFMANN等應用摩爾-庫倫準則,研究了開采擾動下斷層滑移失穩時滑移帶的黏聚力變化情況等。
此外,部分學者從斷層面滑移趨勢入手研究了斷層活化機理。SAINOKI和MITRI等通過推導出新的滑移減弱距離準則對斷層滑移的摩擦阻力進行了標定,考慮斷層傾角、落差、斷層面粗糙度等因素的影響,揭示了逆沖斷層擠壓型地質構造誘發沖擊地壓的機理;TAGHIPOUR等采用數值模擬手段研究了斷層的滑移趨勢,分析了斷層活化時不同部位的應力積累和塑性應變的發展情況;WEI等發現斷層面上不同位置的滑移趨勢大不相同,分析了動摩擦和臨界滑移距離對斷層結構線性滑動弱化的影響。
為進一步揭示開采擾動影響下斷層失穩的誘沖機理,國內外學者對斷層的賦存狀態,以及斷層影響區域位移場、應力場和能量場等多物理場的演化特征也進行了大量研究。目前,對于斷層面位移場的研究手段主要依靠相似模擬試驗和數字散斑監測技術。張寧博、趙善坤等分析了逆沖斷層在單側卸載條件下失穩瞬間位移場的變化規律,探討了斷層滑移失穩機理;張科學、何滿潮等研究了在采動影響下斷層活化過程中斷層附近圍巖的變形情況,并總結出斷層滑移失穩誘發巷道沖擊地壓的機理;宋義敏等通過直剪摩擦滑動試驗分析了斷層滑移失穩過程中位移場的變化規律,發現位移場具有時間間隔滑動特征和空間差異性特征。王宏偉等通過相似模擬試驗給出了斷層面上相對位移的時空演化特征,認為上覆巖層的存在對斷層結構位移具有重要的影響。
開采擾動影響下,應力場的分布與演化特征以及能量場積聚與釋放在斷層結構滑移失穩過程中起著決定性作用。呂進國等認為高水平構造應力和采動應力等多應力場的耦合疊加作用為斷層滑移失穩誘發沖擊地壓提供了動力力源條件;蔡武、竇林名等提出開采擾動引起的水平采空側卸載和豎直方向加載是導致斷層滑移的力學本質;羅浩等采用數值模擬方法與物理相似試驗,得出了圍巖應力場在開采擾動和斷層構造影響下的變化規律;王宏偉等通過建立理論力學模型研究了斷層面上應力場和能量場的時空演化特征,給出斷層滑移判別力學依據;朱廣安等研究了斷層面上應力場的變化情況,并給出了斷層活化可能危險性判別的應力條件;蔣金泉等分析了在巨厚堅硬巖層與逆斷層耦合作用下采動應力的演化規律與斷層活化特征;王濤等通過對開采條件下斷層應力變化特征的研究,發現開采擾動是導致斷層滑移失穩的重要因素之一;譚云亮等發現深部開采活動引起斷層發生滑移失穩時會釋放大量能量,同時煤巖應力瞬間增大,并且認為深部動力災害的監測應當以能量判據為主;姜福興等采用微地震監測手段分析了斷層構造在滑移失穩過程中能量變化特征,總結出構造控制型沖擊地壓的分類依據。
綜上所述,煤礦沖擊地壓的發生與斷層構造賦存條件密切相關,開采擾動下斷層構造的失穩過程可以從圍巖位移場、應力場和能量場的動態演化過程中獲得前兆規律,是研究沖擊地壓精準預測和高效防控的前提。筆者總結了在承擔國家重點研發計劃(2016YFC0801401)、國家自然科學基金(41872205)和北京市自然科學基金(8202041)等課題期間的研究成果,從斷層構造成因、地質賦存狀態和高水平原巖應力環境等角度介紹了沖擊地壓發生時斷層構造的地質賦存環境等特征,建立了考慮工作面開采時的斷層結構力學模型和厚層頂板賦存的物理模型,推導了斷層面上下盤相對滑動位移的逐點積分計算方法,提出了可表征斷層滑動摩擦與能量釋放相關性的斷層結構勢能的概念,分析了開采擾動下斷層構造位移場、應力場和能量場的動態演化特征,系統闡述了斷層構造失穩的發生特點和誘發沖擊地壓的力學機理。筆者的研究結果期望能為我國煤礦沖擊地壓的預測與防治提供理論支撐。
由于巖體受到地殼運動的強烈擠壓作用,使得在斷層等地質構造成型期間積聚了大量的彈性應變能,是其誘發煤礦沖擊地壓的主要動力環境。在現場監測中,盡管可以利用鉆孔應力、地質雷達探測、微震監測和電磁輻射監測等技術手段獲得工作面開采過程中斷層結構附近及上覆巖層的位移場、應力場和能量場的宏觀變化特征。然而,目前能夠表征斷層面多物理場時空演化特征的力學機理解釋尚未清晰,無法準確獲知斷層結構失穩的多物理場動態演化特征和前兆信息。
筆者建立了考慮工作面開采時的斷層結構力學模型和厚層頂板賦存的物理模型,推導了斷層面上下盤相對滑動位移計算的逐點積分方法,提出了表征斷層滑動摩擦與能量釋放相關性的斷層結構勢能的概念,分析了開采擾動下斷層構造位移場、應力場和能量場的動態演化特征,總結出斷層滑移失穩動態過程中摩擦狀態發生轉變的臨界位置和應力條件,系統闡述了斷層構造失穩的發生特點和誘發沖擊地壓的力學機理。
采動條件下斷層構造及其影響區域的應力重分布動態演化過程十分復雜,斷層面的靜態滑移、整體劇烈滑動以及斷層結構位移場的動態演化過程,是研究斷層失穩誘發沖擊地壓的重要前兆信息,而該過程勢必伴隨斷層面上正應力和剪應力的動態響應。同時,由于受到來自側向壓力和上覆巖層下沉壓力的共同作用,斷層滑動是斷層面摩擦效應的體現。在斷層未滑動之前,斷層面由于未達到最大滑動靜摩擦,其斷層區域圍巖彈性應變能處于長期積聚而未釋放的狀態,斷層失穩后的能量釋放也是當前微震事件數、聲發射事件數和電磁輻射強度值監測的重要目標,也可為沖擊地壓事故的發生做出預測。
筆者在研究開采擾動下斷層滑移失穩的力學機理基礎上,研究斷層附近上覆巖層運移與斷層滑動破壞的相互作用,分析斷層側壓系數對斷層面應力場的影響,給出斷層滑動危險性判別條件,著重分析斷層面位移場緩急突變、應力場正剪分化以及能量場激增滯后的動態演化特征,分析斷層結構及影響區域多物理場異常演化現象與斷層滑移失穩誘沖機理之間的聯系,從而獲得能夠預測斷層滑移失穩的前兆特征。
斷層構造的成型過程、地質賦存狀態,是斷層影響區域高水平原巖應力環境的主要原因,了解斷層構造的地質特征有助于探究其失穩誘發沖擊地壓的機理。
經過長期的地殼運動和強烈的擠壓作用,巖層中較為薄弱的位置發生斷裂并與相鄰巖體產生相對運動,從而形成斷層構造。根據斷層構造形成時所受到的應力狀態以及巖石物理性質的差異,可將其分為正斷層、逆斷層和平移斷層。斷層構造成型后沿其走向上不同位置的巖體變形能力和能量積聚程度大不相同,在水平構造應力作用下主斷層兩側巖體上衍生出大量的節理和次生斷裂面,并與主斷層面形成具有一定寬度的斷層帶影響區域。然而,斷層等地質構造由于其深部賦存的特點,在地層中并不容易被發現。目前,國內外主要應用地球物理探測技術對煤層賦存條件和復雜地質構造進行勘探。比較常見的探測方法有地震法和電磁法,例如,彭蘇萍等將三維地震勘探技術應用于淮南煤礦,并證明可以有效地探測小斷層、陷落柱、地質異常體和活動斷層等地質構造,成為煤礦高效安全生產的重要技術手段之一。王揚州等將瞬變電磁法應用到礦井中的全空間范圍,實現井下采空區多角度立體探測技術。因此,預先探測煤層的賦存環境和斷層等復雜地質構造的幾何形態,有助于探究采動影響下斷層等復雜地質構造附近多物理場的時空演化特征,為深入探索斷層滑移的力學機制提供了技術前提條件。
斷層等復雜地質構造在我國各地區煤田中普遍存在。筆者以河南義馬煤業集團股份有限公司義馬礦區、北京昊華能源股份有限公司京西礦區、黑龍江龍煤礦業集團股份有限公司雞西礦區以及開灤(集團)唐山礦業分公司唐山煤礦作為工程案例,分析斷層等復雜地質構造的成因及其賦存地質環境特征。
圖1為河南義馬礦區地質構造。義馬礦區地處褶皺-逆沖斷層的構造區域范圍內,形成具有大量走向、傾向、斜交的斷層和褶曲等復雜地質構造,煤層出現嚴重的分叉與合并現象。義馬礦區從西向東依次分布著楊村煤礦、耿村煤礦、千秋煤礦、躍進煤礦和常村煤礦,向南以F陜石—義馬逆斷層為界線(圖2)。據現場資料統計,千秋煤礦和躍進煤礦中多個工作面的多次沖擊地壓事故均與F逆沖斷層有直接關系。

圖1 義馬礦區地質構造

圖2 義馬礦區各煤礦地理位置及2006—2015年各煤礦沖擊地壓發生頻率
由于地殼運動過程中比較軟弱的煤層受擠壓應力作用易發生斷裂和相對滑動,且在巨厚礫巖層的共同作用下形成壓扭性逆沖斷層F。斷層走向近東西,傾向南,略偏東,淺部傾角75°,深部傾角15°~35°,逆沖面上陡下緩,落差50~500 m,水平錯距120~1 080 m,北接千秋煤礦,向東延入躍進煤礦。根據圖3可知,F斷層緊鄰多個采煤工作面,其存在不僅切斷了巖層的連續性,還伴有應力集中和能量積聚現象,增大了開采的難度。因此,斷層構造的存在被認為是該礦區沖擊地壓發生的主要誘因。

圖3 義馬礦區千秋煤礦21221工作面采掘平面示意
位于北京西部地區的北京昊華能源股份有限公司的京西礦區是我國另一典型的復雜地質構造賦存區域。由于受燕山和喜馬拉雅山多次造山運動、構造運動以及侵蝕作用的劇烈影響,在強擠壓力作用下形成地層倒轉型褶皺和大中型逆沖斷層相當發育,如圖4所示。水平構造力使得礦區內形成多個向斜和背斜構造從而造成地層發生倒轉現象,導致褶皺軸部及兩翼形成了多種次級褶曲構造,同時伴有大量逆沖斷層的產生。以大安山煤礦的大寒嶺背斜和向斜為例,其中背斜傾角范圍為10°~90°,向斜傾角范圍為50°~90°,斷層的最大落差為150 m,最大傾角為77°。圖5給出了京西礦區大安山和木城澗煤礦的構造剖面示意。逆沖斷層的存在導致煤層和覆巖厚度以及煤層傾角急劇變化,在長期的地質演變過程中斷層和褶皺等構造中蘊藏著巨大的變形能,開采擾動條件下極易誘發沖擊地壓等惡性動力災害事故。

圖4 京西礦區地質構造示意

圖5 京西礦區大安山煤礦和木城澗煤礦構造剖面示意
黑龍江龍煤礦業集團股份有限公司的雞西礦區和開灤(集團)唐山礦業分公司的唐山煤礦的煤炭資源同樣賦存著復雜的地質構造區域。圖6為雞西礦區中部礦井開拓示意。由于長期受南北向壓力和左行扭力的作用,雞西礦區內出現南北2條單斜構造拗陷和中部恒山基地隆起的現象,從而構成了兩拗夾一隆的構造格局。井田整體上為一單斜構造,煤層傾角為8°~45°,呈現淺部陡深部緩的總體趨勢。礦區內地質構造復雜程度為中等,以北東向和北西向的2組斜交正斷層為主,東西方向和南北方向的斷層系統為次,主、次級斷層互相交接或切割,使井田呈現類似棋盤格式構造,導致煤層走向和傾角變化極大,增加了巷道的布置難度。

圖6 雞西礦區中部礦井開拓示意
開灤(集團)唐山礦業分公司的唐山煤礦位于河北省唐山市,主要構造為平行于地層走向的斷層,如圖7所示。由于受近似東西向的擠壓作用,井田西部賦存著復雜的褶皺構造和多條縱向斷層,并伴隨著大量的中心斷層的衍生,數條壓扭性斷裂在深部合為陡河大斷裂。由于斷裂區域內構造應力相對較大,從而成為唐山地區的主要的發震構造。統計分析發現,井田向斜軸部是發生沖擊地壓的高危險區,且煤巖沖擊失穩與斷層構造密切相關。

圖7 唐山煤礦地質構造
綜上所述,地殼運動和水平構造擠壓力有助于斷層構造等復雜地質構造的形成和發育,并且褶皺構造通常會造成地層倒轉進而會衍生出大量逆沖斷層,導致地層原巖應力場變得更加復雜,在開采擾動下容易發生斷層滑移失穩從而誘發嚴重的沖擊破壞事故。因此,斷層等復雜地質構造對多物理場的影響是不容忽視的。
原巖應力場是地殼巖石圈應力的原始分布狀態。隨礦井開采深度的增加,原巖應力也隨之增大,導致礦山壓力顯現強烈,為能量積聚提供了理想的內部環境,在增加巷道維護和煤層開采難度的同時,沖擊地壓發生的可能性也隨之增大。因此,考察原巖應力場的分布特征是研究復雜地質構造特征與誘發沖擊動力失穩之間相關性的前提條件。
高水平原巖應力環境是我國煤礦復雜地質賦存環境下誘發沖擊地壓的主要原因之一。以河南義馬礦區和京西礦區為例,根據原巖應力現場測試結果可知,河南義馬礦區內千秋煤礦是以水平應力為主導的構造型原巖應力區,而躍進煤礦原巖應力場則是由垂直應力占主導。考慮斷層構造的存在,水平構造應力依舊是義馬礦區高應力環境中的主要應力分量。與華北地區和中國大陸平均應力水平相比而言,義馬煤田原巖應力水平較高,如圖8所示。尤其是垂直應力分別高于華北地區和全國中國大陸的42.7%和42.1%。以上數據進一步說明了義馬礦區復雜地質構造條件下的高水平原巖應力特征,也是頻繁發生沖擊地壓事故的重要原因。

圖8 義馬礦區原巖應力與華北地區和中國大陸原巖應力對比
京西礦區各礦井的開采水平均屬于深部開采范疇,原巖應力測試結果可知該礦區原巖應力場同樣以水平應力為主導。如圖9所示,京西礦區中主應力隨深度增加呈現不同程度的線性增長趨勢。由于礦區內存在大量的倒轉型褶皺構造和大中型逆斷層等復雜地質構造,使得原巖應力水平要高于北京和華北等地區的平均水平,以深度800 m為例,其中三大主應力分別高于中國大陸地區平均水平的34.1%,75.3%和57.2%。因此,在現有開采水平條件下,高水平原巖應力釋放勢必會引起斷層等構造的瞬間失穩,從而誘發沖擊地壓等動力災害。

圖9 京西礦區原巖應力和北京地區、華北地區和中國大陸原巖應力平均水平對比
所以,義馬礦區和京西礦區的高水平應力環境,以及礦區內頻繁發生的沖擊地壓事故,證實了高水平原巖應力場與沖擊地壓之間的相關性,為沖擊地壓機理的研究提供了支撐。
斷層沖擊地壓作為一種與地質構造賦存狀態相關的沖擊類型,其發生機理是目前沖擊地壓研究領域的熱點。為了更加清晰地了解斷層構造失穩誘發沖擊地壓的力學機理,學者們提出了與斷層相關的誘沖模式來揭示沖擊災害發生的原因。
齊慶新等根據力源條件將斷層沖擊地壓歸納為構造型沖擊地壓,認為斷層等地質構造的存在使得地層中構造應力大幅提高,造成煤巖突然失穩和沖擊破壞。潘一山等根據煤巖失穩的物理特征對沖擊地壓機理進行了分類,總結出斷層錯動型沖擊地壓一般發生在斷層構造附近,并且只有在深部開采和上覆巖層壓力足夠大時才會引起斷層圍巖剪切失穩從而誘發沖擊地壓。姜耀東等從應力狀態導致煤巖突然破壞的本質出發,根據現有的開采工藝和常見地質構造賦存狀態,將由斷層構造或其他結構面的滑移錯動引起的劇烈破壞的動力現象稱為滑移錯動型沖擊地壓。譚云亮等提出深部開采斷層滑移型沖擊地壓,認為深部斷層影響區域受采掘活動影響更為顯著且構造應力明顯大于淺部開采。姜福興等依據微震數據監測結果和構造附近能量釋放情況,將復雜地質構造引起的沖擊地壓分為增壓型和減壓型,當工作面臨近斷層構造開采期間,微震事件在斷層附近出現積聚現象并且不隨工作面遷移時,很可能引起構造增壓型沖擊地壓在斷層附近發生,造成煤巖破壞。
斷層構造的滑移失穩與采掘活動和高水平應力是密切相關的。竇林名等根據煤層開采活動引起的斷層活化類型針對斷層型沖擊地壓提出了4種概念模型,即遠場和近場礦震動載觸發的斷層活化、逼近斷層和遠離斷層開采的采動應力主導型斷層活化。呂進國等根據采掘工作面與斷層構造的相對空間位置不同引起的應力場變化的差異性,將斷層沖擊地壓分為構造應力型、斷層活化型和斷層失穩滑動型等。
筆者在研究斷層構造位移場、應力場和能量場動態演化特征的基礎上,也將重點探討斷層構造和上覆厚層頂板耦合賦存時的工作面逼近斷層構造的采動應力主導型斷層沖擊地壓的力學機理,研究煤巖受采動應力和上覆巖層垂直壓力共同作用下引起的斷層失穩滑動型沖擊地壓。
探究斷層構造滑移失穩區域在開采擾動過程中位移場、應力場和能量場的動態演化特征,建立斷層影響區域沖擊地壓發生前、后多物理場的互饋關系,得到能夠預測斷層失穩誘發沖擊地壓的多物理場前兆特征,研究斷層結構失穩誘發沖擊地壓力學機理,從而達到準確監測以及有效預防沖擊地壓的目的。
斷層影響區上覆巖層運移特征
為了能夠多參量監測開采擾動下斷層影響區域的位移場演化特征,以河南義馬礦區千秋煤礦21221工作面為工程地質背景,在中國礦業大學(北京)深部巖土力學與地下工程國家重點實驗室的深部工程破壞過程模型試驗系統上開展了考慮工作面開采工況的斷層滑動物理相似模擬試驗,如圖10所示。該物理模型的幾何相似比為120,容重相似比為0.67,模型尺寸為1 600 mm×1 600 mm×400 mm。

圖10 斷層滑動物理試驗模型和監測方案
研究顯示,除復雜地質構造賦存條件和高水平原巖應力環境外,上覆巖層失穩也是沖擊地壓發生的重要原因之一。隨著工作面推進,上覆巖層的空間結構發生變化,造成大面積的覆巖運動甚至出現斷裂現象。因此,考察巖層運移與斷層結構滑移失穩相關性對研究斷層沖擊地壓發生機理具有重要意義。
圖11給出了上覆巖層位移隨工作面向斷層推進過程中的變化規律,圖11中表示試驗模型中工作面到斷層的距離。從物理性質上,堅硬巖層的厚度和強度遠大于煤層,上覆巖層通常被視為不易發生斷裂破壞的堅硬頂板。隨著煤層開采過程中,上覆巖層在很長一段時間內處于懸空狀態。上覆巖層發生破裂前,采空區圍巖變形影響范圍持續擴大,巖層位移量也隨之增加,靠近斷層的巖層運動相對于遠處更為活躍,并且在開采擾動條件下發生多次不規律的周期性垮落,導致斷層頻繁受到擾動。

圖11 工作面開采過程中上覆巖層位移變化
此外,研究還發現在開采條件下,斷層并非一直處于滑動狀態,只有當上覆巖層發生垮落破壞時斷層才出現嚴重的滑動。究其原因是由于上覆巖層未破壞前斷層結構處于相對穩定的狀態,而在采動應力和巖層運移共同作用下,斷層結構附近出現了部分離層空間,使得施加到斷層面上的水平荷載減小,造成在斷層構造一側呈現出卸載效應,加劇了斷層的活化程度。同時,斷層滑動又會反過來以非穩態的沖擊和加卸載反作用于煤巖,使得上覆巖層發生更嚴重的垮落破壞,最終誘發沖擊地壓事故。以上結果與分析初步給出了巖層運動和斷層滑移的相互作用的關系。
斷層結構的相對滑動位移的計算
斷層滑移一般表現為上盤與下盤的相對錯動,不僅會引起周圍圍巖運動的改變,還會誘發嚴重沖擊地壓事故。因此,進一步研究在開采擾動影響下斷層面上下盤之間的相對滑移過程,分析發生最危險滑動的位置及其演化特征,能夠更加全面地了解斷層的滑移機制。
在煤礦開采過程中,斷層周圍煤巖受上覆巖層影響呈現整體向下移動的趨勢,隨著開采工作面逐漸靠近斷層,下沉位移量逐漸增大,最終引起斷層構造上、下盤發生錯動。圖12給出了30 cm寬斷層帶區域的滑移情況。研究發現,斷層構造的滑動位移最先出現在斷層頂部位置,并逐漸向下擴展到煤層附近,而斷層最大滑移量最終在斷層與厚層頂板交界處的中部區域集中發展。從圖13可以發現,斷層面上的位移隨著開采過程中由緩慢增加轉為突然急劇增加的現象,斷層面最大位移量也由斷層頂部突然轉移至斷層中部區域,造成斷層結構滑移失穩破壞。

圖12 煤礦開采過程中斷層面上位移變化情況

圖13 工作面推進過程中斷面上斷層位移變化
為進一步解釋斷層面滑動位移突變的現象,筆者提出了基于逐點積分的相對滑移算法來求解斷層面上盤和下盤之間的相對位移。斷層附近區域上盤與下盤在水平和垂直方向上的相對位移矢量如圖14所示。

圖14 斷層面上下盤相對位移計算示意


(1)

(2)

計算表明,開采擾動下斷層面的相對滑動是由局部靜態滑移向整體劇烈滑動突然轉變的非線性動態演化過程。在這個動態演變過程中,斷層附近的圍巖運動經歷了沿垂直向下、近似垂直于斷層和平行于斷層3個運動階段,如圖15所示。圖16給出了工作面在開采過程中斷層上下盤相對滑動的動態演化過程。在開采前期階段,斷層上下盤由于覆巖壓力作用整體向下移動。隨著工作面的推進,斷層面上開始出現局部滑動,且斷層下盤相對于上盤沿斷層傾向向下運動,此時斷層面上不同位置的相對位移呈現基本同步的緩慢增長趨勢,如圖17所示。臨近斷層開采時,斷層上下盤的擠壓效應代表了此時斷層結構能量處于積聚狀態,斷層中部區域的相對位移量大幅度增加,預示斷層發生嚴重的整體滑移破壞并釋放大量的應變能。

圖15 斷層面上下盤相對滑移矢量

圖16 工作面開采過程中斷層相對滑動動態演化

圖17 工作面開采過程中斷層下盤相對上盤的位移曲線
綜上可知,斷層面最大相對位移量出現在斷層中部靠近煤層的區域,而相對位移緩急突變發生時伴隨著上覆巖層嚴重垮落現象。上覆巖層運移失穩會瞬時釋放大量應變能,從而促進了厚層頂板與斷層構造交界區域的相對滑動,特別是未破壞巖層的失穩破壞可能會進一步引起嚴重的斷層滑移失穩。筆者從斷層上下盤相對位移角度入手揭示了開采擾動下斷層滑移由局部滑移轉向整體失穩的動態過程,證實了斷層滑移失穩對上覆巖層運移的影響作用。
斷層面應力分布的監測結果
采動應力的存在打破了采場原位應力的平衡。工作面支承壓力與斷層處高構造應力疊加形成應力集中,導致斷層滑移失穩。
因此,進一步探究采動影響下斷層構造應力場的演化特征,能夠更加深入地了解研究斷層滑移失穩的力學機制。
圖18為部分監測點應力隨工作面開采變化曲線(監測點位置如圖10所示)。在開采初期階段,工作面距離斷層較遠時(距斷層處100 cm),斷層面上的應力開始發生緩慢變化,斷層開始受到擾動,與上述斷層滑移情況基本一致,可將該位置視為斷層影響區和非影響區的分界線。

圖18 工作面推進過程中斷層面上正應力和剪應力變化
臨近斷層開采時(距斷層處60 cm),斷層因開采活動影響頻繁受到擾動,特別是斷層中部區域的應力變化尤其劇烈。當工作面在斷層附近開采時(距斷層處20 cm),斷層發生嚴重的滑移失穩,此時應力水平明顯降低。
開采擾動影響下,斷層面正應力和剪應力從同步變化轉為突然分化的急劇變化。斷層面上正應力變化相比剪應力而言更為強烈,且近煤層處的應力值也要大于遠離煤層處。除外,在臨近斷層開采時,斷層面中部區域正應力和剪應力呈現截然相反的變化趨勢,即正應力增加而剪應力降低或剪應力增加而正應力降低的分化現象。
考慮水平構造應力和采動應力的影響,斷層構造在側向壓力和上覆巖層下沉壓力的共同作用下而處于壓剪受力狀態。斷層面上的摩擦力隨著正應力增加而增大,剪應力和正應力持續增大導致斷層結構由靜態摩擦轉為動態摩擦,此時斷層面上剪應力會相應減小。加之斷層中部區域巖層被壓實且不易松,進而出現正應力增加而剪應力降低的分化現象。為進一步證實這一監測結果,筆者提出了考慮采動影響的斷層滑動的力學模型,解釋了正應力和剪應力的分化特征,并做詳細闡述。
考慮采動影響的斷層結構的力學模型
由于斷層面上應力場分布和變化較為復雜,僅依靠現場監測和數值仿真技術還不能完全揭示開采擾動下斷層面上的應力場演化特征。為進一步驗證和解釋斷層面上正應力和剪應力的分化現象,筆者建立了考慮采動影響的斷層表面應力分布和斷層滑動特征的力學模型,如圖19所示。綜合考慮上覆巖層壓力、重力附加應力和采動應力的影響,得到式(3)給出的斷層面上任意一點處的應力分量,進而計算得到如式(4)所示的斷層面上正應力和剪應力的表達式。

(3)

(4)
其中,,和分別為沿方向的正應力分量、沿方向的正應力分量和剪切應力分量;和為斷層面任意點的正應力和剪應力;和分別為任一點的水平和鉛直距離;和分別為模型的長度和高度;為采動應力峰值點到斷層的水平距離;為垂直分布應力;為巖石密度;為重力加速度;為斷層傾角;()和()分別為重力附加應力和采動應力分布的水平應力函數;()d為距坐標原點處微小長度d上所受的微小集中力。
圖19中,為斷層面上任意點到模型頂部的鉛直距離;′為力學模型中工作面與斷層之間的距離;為側向系數。

圖19 考慮采動影響的斷層面上任意點應力計算的力學模型
斷層面應力場的動態演化特征
在采動應力影響下,斷層面上不同位置應力分布各異,圖20給出了隨工作面向斷層推進過程中斷層面應力值的分布和變化規律。圖21為斷層面應力計算的試驗校對。對比試驗與理論結果可知,關于斷層面上正應力和剪應力的分布特征,理論模型計算結果和試驗結果基本吻合。結果顯示,斷層面下段區域應力水平高于斷層面上段區域,隨著工作面的推進,應力變化幅度更加明顯。沿斷層面向下,斷層面的正應力值和剪應力值變化呈相反趨勢,出現正應力增加而剪應力降低或剪應力增加而正應力降低的分化現象。應力場變化的轉折點代表著斷層表面的最小正應力值和最大剪應力值,說明了應力極值與斷層滑移失穩具有一定的聯系。

圖20 斷層面應力分布的理論模型計算結果

圖21 斷層面應力分布的試驗結果校對
根據正應力和剪應力的極值變化,可將斷層表面劃分為R,R和R三個區域,如圖22所示。工作面在到達斷層擾動影響區繼續向前推進時,斷層面的應力極值逐漸向底部移動。沿斷層面向下,根據正應力和剪應力的變化狀態,得出R為斷層上、下盤的靜摩擦穩定狀態區域,R區域中剪應力極值位置處于即將動態滑移的臨界平衡狀態,在R區域斷層為已經發生了動態相對滑動。因此,斷層滑動面的剪應力極值出現的位置可視為斷層面從靜態轉向動態滑動的臨界區域,同時也是斷層面上最先發生滑動的位置。


圖22 開采過程中斷層面上的應力演化
綜上可知,當工作面臨近斷層開采時,由于采動應力和上覆巖層的影響,斷層結構面和厚層頂板的交界區域的應力水平較高,此時該區域積累了大量的應變能,為斷層滑移失穩提供了應力條件。此外,上述結果還給出了應力極值分布演化特征,并從摩擦狀態轉變的角度分析了斷層滑動的力學機制。
斷層結構勢能概念的提出
開采擾動下斷層應力的變化會引起斷層表面應變能的積累和釋放。斷層周圍煤巖能量變化實質上就是從穩定積聚狀態逐漸向非穩定釋放狀態轉化的過程。由于較高應力水平、斷層構造以及開采條件等諸多影響因素,造成開采過程中煤巖的能量變化具有明顯的時空特征。在現場監測中,微震事件數、聲發射事件數和電磁輻射強度值等數據在沖擊地壓發生前常出現監測數據消失的現象。為進一步解釋這一工程現象,在斷層面應力分布的理論模型基礎上提出了斷層結構勢能的概念,并對斷層結構勢能的時空演化特征進行了分析。
根據彈性理論中的應力變分方程,通過斷層面上的應力來計算斷層結構勢能為

(5)

斷層結構勢能演化特征
圖23為斷層面上結構勢能的計算結果。工作面過斷層前,斷層結構勢能分布特征與剪應力分布特征相似,對照圖22可知,結構勢能最大值與剪應力最大值類似,均出現在R區域的底部附近,并且在R區域逐漸減小。當采煤工作面穿過斷層后,斷層結構勢能變化明顯,斷層頂部結構勢能突然增大,R區域的結構勢能急劇下降。從圖24可以看出,隨著工作面向斷層推進,斷層面上部區域的結構勢能逐漸減小,而中下部區域結構勢能則呈增大的趨勢。

圖23 煤層開采過程中斷層面上結構勢能的分布

圖24 工作面開采過程中斷層面上結構勢能演化
綜上可知,開采工作面過斷層前結構勢能最大的位置出現在斷層結構面和上覆厚層頂板的交界區域,同時也是最大剪應力出現的位置,即斷層上下盤運動由靜摩擦轉變為動摩擦的臨界區域,可將該位置認為是斷層滑動的起始點。由于斷層滑移錯動導致斷層表面應力對開采活動產生不同的響應,從應變能角度解釋了正應力和剪應力變化過程中的分化特征。
由于斷層區域的上覆巖層較為堅硬不易破壞,能量激增前斷層處于相對穩定的狀態,導致監測數據缺失現象,但此時能量則一直處于積聚而未能釋放的狀態。在開采擾動和斷層滑移作用下,未破壞上覆巖層發生失穩導致斷層頂部的能量積聚和瞬時釋放,從而引起斷層的嚴重滑動破壞誘發沖擊地壓。筆者從結構勢能角度解釋了斷層滑移失穩過程中能量的穩態積聚和非穩態釋放的變化,同時證實了上覆巖層運移和破裂對斷層滑移失穩具有重要的貢獻。
開采擾動下斷層結構發生滑移錯動引起斷層圍巖位移場、應力場和能量場發生改變,繼而誘發嚴重的沖擊地壓。筆者在研究斷層結構滑移失穩多物理場演化特征的基礎上,進一步分析采動影響下斷層沖擊地壓的力學機理。
斷層失穩誘沖的位移場突變特征
由前文可知,在整個煤層開采過程中,斷層的滑動位移通常會經歷靜止、緩慢滑移和快速滑動3個階段。斷層最大位移出現位置也發生了由斷層頂部轉向中部的空間變化。在臨近斷層開采時,斷層位移量突然明顯增加,代表斷層發生嚴重滑移失穩,可將位移場這種由穩定到急劇增加的突變現象作為斷層滑移失穩的前兆特征。
斷層面上位移場的突變特征還表現在采動應變和應變率的監測結果中,圖25和圖26分別為斷層面應變和應變率的變化曲線。當斷層開始擾動后,斷層附近圍巖應變出現急劇上升和穩定增加的相互交替變化階段,隨后斷層發生嚴重滑移,采動應變再次急劇減小。同時,斷層附近應變率也由低水平穩定狀態轉變為急劇波動階段,預示上覆巖層中裂縫急劇發展,最終引起上覆巖層垮落破壞。因此,可將斷層面上采動應變急劇增加和維持穩定交替出現的階段或應變率急劇變化的現象視為斷層滑移的另一前兆特征。

圖25 煤層開采過程中斷層面應變變化

圖26 斷層面上應變率變化曲線以及試驗對比
斷層失穩誘沖的應力場分化特征
斷層面上中部區域的應力場正剪分化特征表現為正應力增加而剪應力降低或剪應力增加而正應力降低的現象。為進一步分析斷層面應力場這一特征,圖27給出了斷層中部應力場變化以及對應的位移場變化。正應力和剪應力的共同作用使得斷層面發生滑移失穩,此時由于水平應力的擠壓作用,正應力持續增大,代表能量逐漸積聚,而剪應力隨斷層滑移錯動開始呈下降趨勢。另外,當斷層出現滑動時伴隨著正應力的急劇下降和剪應力的突然增加,此時能量開始釋放,正應力的減小會導致斷層表面剪切強度的降低,從而加速了嚴重的斷層滑移。隨著工作面向斷層方向推進,斷層相對滑移量逐漸增加,最大位移量出現的位置也開始向下部轉移,直到斷層滑移失穩造成斷層區域圍巖的破壞,此時正應力和剪應力均處于較小值。因此,可以將應力場的正剪分化現象視為斷層滑移的前兆特征。

圖27 工作面回采期間斷層中部區域正應力和剪應力變化及相對應的斷層滑移情況
為進一步解釋斷層中部發生劇烈相對滑動的現象,根據探地雷達監測到的物理模型前、后兩側的電磁波幅值差值,得出開采過程中斷層中部區域振幅衰減變化,如圖28所示。煤層開采過程中,正應力和剪應力開始發生分化現象時,振幅衰減量從極小值開始增長并隨之保持穩定。在這個過程中,斷層附近圍巖中有大量裂隙產生并擴展,正應力和剪應力的急劇變化同時也說明了斷層由局部滑動演化為整體滑動。以上從巖石損傷破壞程度的角度驗證了應力場正剪分化現象作為斷層滑移前兆信息的正確性。

圖28 開采過程中斷層面正應力、剪應力曲線及相應位置的振幅衰減變化
斷層失穩誘沖的能量場激增特征
聲發射技術能夠監測到能量場在煤層開采過程中的釋放情況,從而判斷斷層滑移破壞程度。
如圖29所示,聲發射事件數在工作面推進過程中出現多次激增,分別是斷層加載初期、臨近距斷層開采和到達斷層處開采時。斷層結構釋放的能量主要聚集在斷層與上覆厚層頂板交界處的中部區域,這與斷層最大滑動位移和剪應力最大值出現位置近似,同斷層結構勢能計算結果基本一致。

圖29 工作面開采事件數分布
此外,聲發射事件在激增前會出現事件缺失的現象,與工程監測結果基本一致。以河南義馬礦區的重大沖擊地壓案為例。圖30給出了2011年義馬礦區千秋煤礦和2010年躍進煤礦的電磁輻射強度和支架阻力的監測結果。電磁輻射強度和支架阻力值呈現基本一致的增大與減小循環變化特征,表明沖擊地壓發生前斷層和上覆巖層積聚了大量變形能,而電磁輻射強度和支架阻力突降說明大量應變能開始釋放,證實了斷層結構滑移失穩時能量場的激增現象。圖31為2011年義馬礦區千秋煤礦監測到的微震事件頻次與能量釋放隨時間變化的現場監測結果。監測結果中存在微震事件的缺失現象,并且對應時間段的微震事件頻次也會相對較低,此時能量處于積聚而未釋放的狀態,與本文斷層結構勢能計算結果基本一致。

圖30 千秋煤礦和躍進煤礦電磁輻射強度和支架阻力變化曲線

圖31 千秋煤礦微震能量及頻次累計情況
斷層滑移失穩過程中能量場變化與應力場、位移場緊密相關,圖32和33分別給出了斷層面聲發射事件與位移場、應力場之間的變化關系。研究表明,斷層面上的滑移并非時時刻刻在發生的,只有在滑移較為嚴重并釋放大量能量時才可捕捉到聲發射事件。而上覆巖層突然垮落或斷層發生嚴重滑移之前,巖層中裂縫較少,聲發射事件也相對穩定,此時能量場處于一直積聚的狀態。因此,根據位移場緩急突變和能量場激增現象先后出現,可將聲發射事件急劇增加前的恒定不變階段視為斷層滑移的前兆信息。

圖32 煤層開采過程中聲發射事件和最大斷層位移的綜合分析

圖33 斷層面上剪應力和聲發射事件數之間的聯系
根據斷層面聲發射事件與應力場的變化關系可知,當工作面開采至斷層附近時,斷層聲發射事件數激增,此時斷層面剪應力出現逐步降低后陡然增加的現象,斷層面即將發生滑動。當工作面開采至斷層處,應力降至最小值,聲發射事件數再次激增,斷層滑移失穩并釋放出大量的能量。因此,斷層面應力在逐步降低的過程中陡然增加和聲發射事件數恒定不變時突然激增的特征可作為斷層滑移失穩的前兆信息。此外,從上述結果中可以發現,能量場的激增現象一般滯后于位移場的突變現象和應力場的分化現象。
綜上可知,斷層的滑動失穩過程中位移場的緩急突變,應力場的正剪分化以及能量場的激增滯后出現的位置基本一致,斷層面最大位移量代表斷層出現嚴重滑移破壞,應力場正剪分化現象是斷層的滑移錯動造成的異常響應,并且伴隨著大量積聚的應變能瞬間釋放,以上結果均表明了斷層結構在開采擾動影響下多物理場的互饋關系。
斷層滑移失穩滑動危險性判別
斷層滑移失穩滑動危險性主要與斷層構造有關,尤其是斷層構造的幾何形態對滑移失穩具有一定的影響作用。比如,斷層切向剛度主要控制著斷層剪切變形滑移,隨著斷層切向剛度的增加,剪應力值變化越明顯,斷層滑移危險性越高。落差與斷層正應力呈正相關關系,落差越大,正應力越大,斷層更易發生滑移失穩。另外,斷層傾角不僅代表了斷層的傾斜程度,描述了斷層的大致形態和走向,還是影響斷層滑移失穩的一個重要指標。
為考察斷層構造的傾角與斷層滑移失穩的關系,圖34給出了不同斷層傾角作用下位移場的變化。研究發現,斷層相對位移在開采后期突變隨傾角增大而越發明顯,尤其是傾角為75°時。開采初期釋放能量較少代表了應變能在持續積聚,此時斷層結構面相對位移幾乎為0,并且隨著工作面推進呈現緩慢增長趨勢。當工作面臨近斷層開采時,大量應變能的突然釋放導致斷層面相對位移的急劇增加,最終造成嚴重的斷層滑移失穩。

圖34 不同斷層傾角模型的斷層附近上覆巖層相對位移變化
為進一步解釋采動過程中傾角變化對斷層結構相對位移的影響,提出了斷層滑移危險性的判別條件。圖35給出了斷層面不同位置相對位移的變化情況。斷層面底部區域相對位移值較大,且斷層的相對位移隨傾角變化逐漸減小甚至趨近于0,表明斷層上、下盤由于強擠壓作用下彈性應變能處于長期積聚而未釋放的狀態。斷層的傾角越大,出現擠壓應變積聚的區域范圍越大,斷層滑移失穩時瞬時釋放的能量越多,誘發沖擊地壓的危險性越高。因此,當斷層傾角較大時,來自水平側向壓力和上覆巖層下沉壓力作用越顯著,斷層結構處應力集中和能量積聚越明顯,斷層滑動更加劇烈,發生沖擊地壓的危險性越高。

圖35 煤層開采過程中不同斷層傾角斷層附近上覆巖層的相對位移
高應力環境下斷層構造失穩誘沖機理
復雜多變的地質因素對斷層面應力場演化和滑移量突變具有重要的影響作用,尤其是水平構造應力對斷層的擠壓作用。在水平構造應力作用下,隨著采空區范圍增大,斷層上下盤的相對滑動在一定程度上穩步釋放了水平構造應力擠壓積聚的應變能。圖36為斷層滑移量和滑移速率(每計算100步作為單位時間)隨工作面向斷層方向推進的變化情況,發現在臨近斷層面開采時,滑移速率絕對值突增和滑移量陡增,此時斷層滑移危險極高。

圖36 斷層滑移量及滑移速率
為進一步研究水平構造應力對斷層滑移失穩的影響,圖37給出了不同水平側壓系數作用下應力場的演化特征。當側壓系數>1.0時,斷層下段的正應力值較大,剪應力極大值也處于較高水平,此時上盤與下盤之間擠壓力明顯,為逆斷層的形成創造了壓力條件。對于側壓系數<1.0的情況,正應力在斷層上段隨側壓系數的減小而增大,下段應力水平相對較低,有利于正斷層的形成。因此,筆者初步給出了正、逆斷層形成的力學基本條件,但對于斷層形成的力學機制還需進一步的研究和推導。

圖37 不同側壓系數下斷層面正應力和剪應力的變化

斷層構造和上覆厚層頂板耦合誘沖機理
斷層構造和上覆厚層頂板賦存的復雜地質環境中,在采動應力影響作用下,巖層運移和斷層滑移相互作用增加了沖擊地壓發生的可能性。由于堅硬頂板的厚度和強度相較于煤層很大,當巖層發生下沉或垮落現象,采場發生沖擊破壞的危險性更高。
根據前文可知,上覆巖層運移與斷層滑移失穩兩者存在相互作用的關系。上覆巖層擠壓斷層產生滑動,而斷層的滑移又作用于上覆巖層致使其大范圍運動甚至垮落,二者的聯合作用很可能會誘發沖擊地壓。圖38給出了義馬礦區沖擊地壓成因示意。逆沖斷層作為地質環境中的最不穩定的因素,不僅反轉了煤層方向,并且在其成型期間積累了巨大的應變能,在開采擾動下斷層積聚的應變能會在短時間內突然釋放。由于上覆厚層頂板對煤層產生高強度的動壓力,導致能量場在上覆巖層與斷層構造的交界處出現激增滯后的異常現象。因此,煤田高水平原巖應力為上覆巖層運移和斷層的活化創造了有利條件,綜合以上多因素耦合作用誘發沖擊動力災害。

圖38 義馬礦區沖擊地壓成因示意
此外,根據探地雷達監測到未破壞礫巖的低波速帶可知,如圖39所示。由于未受破壞巖層可以保護斷層不受開采擾動和頂板垮落破壞的影響,斷層在開采前期存在穩定階段。當巖層發生沉降和分離時釋放了大量的應變能,從而加速了未破壞巖層與斷層構造相互影響區域的相對滑動。當上覆巖層坍塌或失穩時,反而會進一步引起更加劇烈的斷層相對滑動導致沖擊地壓的發生。上述研究結果說明,在復雜斷層構造賦存、高水平原巖應力環境和開采擾動條件下,上覆巖層運動與斷層滑移的耦合作用是斷層發生滑移失穩的重要原因。

圖39 煤層開采過程中雷達探測結果及巖層位移矢量
(1)在開采擾動下,斷層面相對滑動位移最早出現在斷層頂部并逐漸向下發展,且由緩慢滑動變化為突然急劇滑移的非線性動態失穩過程,最后在斷層結構面和上覆厚層頂板交界的中部區域發生最大相對滑動位移。斷層面上不同位置應力分布各異,且正應力變化相比剪應力更為劇烈,尤其是正應力和剪應力在同步演化發展過程中發生了分化的現象,并表現出正應力增加而剪應力降低或剪應力增加而正應力降低的特征。斷層失穩前,結構勢能處于不斷積聚而未釋放的狀態,斷層失穩后,結構勢能出現由幾乎為零轉為突然激增的現象。
(2)斷層結構面和上覆厚層頂板的交界區域出現的最大滑動靜摩擦向滑動失穩摩擦的轉換,是斷層面位移場突變和應力場分化的時刻,也是斷層結構勢能從零到激增的時刻。因此,斷層面位移場緩急突變、應力場正剪分化以及能量場激增滯后的現象均可視為斷層滑移失穩的前兆特征。
(3)在復雜斷層構造賦存、高水平原巖應力環境以及開采擾動條件下,由于斷層附近厚層頂板下沉運動和垮落破壞的影響,加劇了斷層的活化和滑移失穩程度,導致斷層結構面和厚層頂板的交界區域出現了由滑動靜摩擦轉為劇烈動態摩擦的現象。斷層最大位移量的出現預示著斷層的嚴重失穩,釋放該區域積聚的大量應變能,從而誘發沖擊地壓災害。
(4)該研究是對煤礦復雜地質構造賦存環境下沖擊地壓發生機理和斷層構造失穩的多物理場前兆特征的研究探索。斷層失穩過程中,斷層滑動面及周邊圍巖的幾何結構和溫度場也應有相應的變化,本文的成果為這一方向研究打下理論基礎,同時為沖擊地壓的預測和防治提供理論支撐。