• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    煤礦沖擊地壓與斷層構(gòu)造失穩(wěn)的多物理場(chǎng)互饋機(jī)制研究進(jìn)展

    2022-04-06 08:53:36王宏偉石瑞明姜耀東
    煤炭學(xué)報(bào) 2022年2期
    關(guān)鍵詞:剪應(yīng)力應(yīng)力場(chǎng)巖層

    王宏偉,王 晴,石瑞明,姜耀東,2,田 政

    (1.中國礦業(yè)大學(xué)(北京) 力學(xué)與建筑工程學(xué)院,北京 100083;2.中國礦業(yè)大學(xué)(北京) 煤炭資源與安全開采國家重點(diǎn)實(shí)驗(yàn)室,北京 100083;3.中國礦業(yè)大學(xué)(北京) 深部巖土力學(xué)與地下工程國家重點(diǎn)實(shí)驗(yàn)室,北京 100083)

    沖擊地壓作為典型的動(dòng)力災(zāi)害,由于其頻繁發(fā)生的危害性、極短前兆的突發(fā)性、機(jī)理研究和預(yù)測(cè)防治的復(fù)雜性等特點(diǎn),給煤炭資源安全高效開采帶來巨大的挑戰(zhàn),同時(shí)也是國內(nèi)外學(xué)者高度關(guān)注的熱點(diǎn)課題。

    在眾多誘發(fā)沖擊地壓的因素中,斷層等復(fù)雜地質(zhì)構(gòu)造賦存環(huán)境是沖擊地壓等動(dòng)力災(zāi)害的重要誘因。斷層失穩(wěn)誘發(fā)沖擊地壓是由于斷層構(gòu)造等不連續(xù)面的存在,在開采擾動(dòng)下出現(xiàn)突然的相對(duì)滑動(dòng)和能量的瞬間釋放,從而導(dǎo)致礦井劇烈振動(dòng)和破壞的動(dòng)力現(xiàn)象的發(fā)生。與其他典型的沖擊地壓類型相比,斷層失穩(wěn)誘發(fā)的沖擊地壓所影響的范圍更加廣泛,災(zāi)害發(fā)生時(shí)所釋放的能量強(qiáng)度更大、震級(jí)更高。

    近年來,沖擊地壓等典型動(dòng)力災(zāi)害事故在煤礦開采中頻繁發(fā)生。例如,2018-10-20山東能源集團(tuán)龍鄆煤業(yè)有限公司發(fā)生嚴(yán)重的沖擊地壓事故,由于受采掘、施工擾動(dòng)和斷層滑動(dòng)等影響導(dǎo)致地應(yīng)力更加集中從而誘發(fā)動(dòng)力災(zāi)害,造成數(shù)名人員傷亡和巨大經(jīng)濟(jì)損失;2019-06-09吉林煤業(yè)集團(tuán)龍家堡礦業(yè)有限責(zé)任公司的沖擊地壓事故,該事故的直接原因是高水平構(gòu)造應(yīng)力、采掘活動(dòng)影響和斷層活化等多因素共同導(dǎo)致,造成巷道嚴(yán)重破壞和人員傷亡;2019-08-02開灤(集團(tuán))有限責(zé)任公司唐山礦業(yè)分公司井田內(nèi)斷層和褶皺構(gòu)造附近發(fā)生了井下沖擊地壓事故,給煤礦的安全生產(chǎn)帶來了巨大損失;此外,2020-02-22山東新巨龍能源有限責(zé)任公司在開采工作面與斷層構(gòu)造形成的強(qiáng)沖擊危險(xiǎn)區(qū)域發(fā)生了較大沖擊地壓事故等。事故分析報(bào)告顯示,上述沖擊地壓事故的發(fā)生多與斷層等地質(zhì)構(gòu)造賦存狀態(tài)影響相關(guān),事故高發(fā)區(qū)域地應(yīng)力集中明顯。因此,研究斷層結(jié)構(gòu)失穩(wěn)誘發(fā)沖擊地壓的特征和機(jī)理對(duì)沖擊地壓的預(yù)測(cè)和防治至關(guān)重要。

    關(guān)于斷層構(gòu)造和沖擊地壓的相關(guān)性,國內(nèi)外學(xué)者從各個(gè)角度研究了沖擊地壓機(jī)理,并提出一系列經(jīng)典理論,為斷層構(gòu)造失穩(wěn)誘沖機(jī)理研究奠定了理論基礎(chǔ)。MICHALSHI最早提出了在臨近斷層開采時(shí)更容易誘發(fā)沖擊地壓的觀點(diǎn),OHNAKA給出了斷層滑移動(dòng)態(tài)摩擦定律假說;RICE等則認(rèn)為沖擊地壓的發(fā)生主要是因?yàn)閿鄬用娴哪Σ梁图羟凶饔?,并分析了斷層滑移過程中溫度場(chǎng)變化特征;于廣明,謝和平等將分形理論應(yīng)用于斷層活化研究,分析了斷層面分型界面效應(yīng)與巖體破壞的關(guān)系;考慮斷層構(gòu)造的影響,王學(xué)濱、林遠(yuǎn)東等應(yīng)用應(yīng)變梯度塑性理論得到了斷層帶與圍巖失穩(wěn)的判據(jù)及解釋,并推導(dǎo)了斷層滑移錯(cuò)動(dòng)位移理論公式;宋振騏在“實(shí)用礦山壓力控制”理論的基礎(chǔ)上,建立了過斷層開采的控制力學(xué)模型以及相關(guān)控制準(zhǔn)則;紀(jì)洪廣等將庫侖應(yīng)力的概念引入開采擾動(dòng)下斷層滑移失穩(wěn)分析中;HOFMANN等應(yīng)用摩爾-庫倫準(zhǔn)則,研究了開采擾動(dòng)下斷層滑移失穩(wěn)時(shí)滑移帶的黏聚力變化情況等。

    此外,部分學(xué)者從斷層面滑移趨勢(shì)入手研究了斷層活化機(jī)理。SAINOKI和MITRI等通過推導(dǎo)出新的滑移減弱距離準(zhǔn)則對(duì)斷層滑移的摩擦阻力進(jìn)行了標(biāo)定,考慮斷層傾角、落差、斷層面粗糙度等因素的影響,揭示了逆沖斷層擠壓型地質(zhì)構(gòu)造誘發(fā)沖擊地壓的機(jī)理;TAGHIPOUR等采用數(shù)值模擬手段研究了斷層的滑移趨勢(shì),分析了斷層活化時(shí)不同部位的應(yīng)力積累和塑性應(yīng)變的發(fā)展情況;WEI等發(fā)現(xiàn)斷層面上不同位置的滑移趨勢(shì)大不相同,分析了動(dòng)摩擦和臨界滑移距離對(duì)斷層結(jié)構(gòu)線性滑動(dòng)弱化的影響。

    為進(jìn)一步揭示開采擾動(dòng)影響下斷層失穩(wěn)的誘沖機(jī)理,國內(nèi)外學(xué)者對(duì)斷層的賦存狀態(tài),以及斷層影響區(qū)域位移場(chǎng)、應(yīng)力場(chǎng)和能量場(chǎng)等多物理場(chǎng)的演化特征也進(jìn)行了大量研究。目前,對(duì)于斷層面位移場(chǎng)的研究手段主要依靠相似模擬試驗(yàn)和數(shù)字散斑監(jiān)測(cè)技術(shù)。張寧博、趙善坤等分析了逆沖斷層在單側(cè)卸載條件下失穩(wěn)瞬間位移場(chǎng)的變化規(guī)律,探討了斷層滑移失穩(wěn)機(jī)理;張科學(xué)、何滿潮等研究了在采動(dòng)影響下斷層活化過程中斷層附近圍巖的變形情況,并總結(jié)出斷層滑移失穩(wěn)誘發(fā)巷道沖擊地壓的機(jī)理;宋義敏等通過直剪摩擦滑動(dòng)試驗(yàn)分析了斷層滑移失穩(wěn)過程中位移場(chǎng)的變化規(guī)律,發(fā)現(xiàn)位移場(chǎng)具有時(shí)間間隔滑動(dòng)特征和空間差異性特征。王宏偉等通過相似模擬試驗(yàn)給出了斷層面上相對(duì)位移的時(shí)空演化特征,認(rèn)為上覆巖層的存在對(duì)斷層結(jié)構(gòu)位移具有重要的影響。

    開采擾動(dòng)影響下,應(yīng)力場(chǎng)的分布與演化特征以及能量場(chǎng)積聚與釋放在斷層結(jié)構(gòu)滑移失穩(wěn)過程中起著決定性作用。呂進(jìn)國等認(rèn)為高水平構(gòu)造應(yīng)力和采動(dòng)應(yīng)力等多應(yīng)力場(chǎng)的耦合疊加作用為斷層滑移失穩(wěn)誘發(fā)沖擊地壓提供了動(dòng)力力源條件;蔡武、竇林名等提出開采擾動(dòng)引起的水平采空側(cè)卸載和豎直方向加載是導(dǎo)致斷層滑移的力學(xué)本質(zhì);羅浩等采用數(shù)值模擬方法與物理相似試驗(yàn),得出了圍巖應(yīng)力場(chǎng)在開采擾動(dòng)和斷層構(gòu)造影響下的變化規(guī)律;王宏偉等通過建立理論力學(xué)模型研究了斷層面上應(yīng)力場(chǎng)和能量場(chǎng)的時(shí)空演化特征,給出斷層滑移判別力學(xué)依據(jù);朱廣安等研究了斷層面上應(yīng)力場(chǎng)的變化情況,并給出了斷層活化可能危險(xiǎn)性判別的應(yīng)力條件;蔣金泉等分析了在巨厚堅(jiān)硬巖層與逆斷層耦合作用下采動(dòng)應(yīng)力的演化規(guī)律與斷層活化特征;王濤等通過對(duì)開采條件下斷層應(yīng)力變化特征的研究,發(fā)現(xiàn)開采擾動(dòng)是導(dǎo)致斷層滑移失穩(wěn)的重要因素之一;譚云亮等發(fā)現(xiàn)深部開采活動(dòng)引起斷層發(fā)生滑移失穩(wěn)時(shí)會(huì)釋放大量能量,同時(shí)煤巖應(yīng)力瞬間增大,并且認(rèn)為深部動(dòng)力災(zāi)害的監(jiān)測(cè)應(yīng)當(dāng)以能量判據(jù)為主;姜福興等采用微地震監(jiān)測(cè)手段分析了斷層構(gòu)造在滑移失穩(wěn)過程中能量變化特征,總結(jié)出構(gòu)造控制型沖擊地壓的分類依據(jù)。

    綜上所述,煤礦沖擊地壓的發(fā)生與斷層構(gòu)造賦存條件密切相關(guān),開采擾動(dòng)下斷層構(gòu)造的失穩(wěn)過程可以從圍巖位移場(chǎng)、應(yīng)力場(chǎng)和能量場(chǎng)的動(dòng)態(tài)演化過程中獲得前兆規(guī)律,是研究沖擊地壓精準(zhǔn)預(yù)測(cè)和高效防控的前提。筆者總結(jié)了在承擔(dān)國家重點(diǎn)研發(fā)計(jì)劃(2016YFC0801401)、國家自然科學(xué)基金(41872205)和北京市自然科學(xué)基金(8202041)等課題期間的研究成果,從斷層構(gòu)造成因、地質(zhì)賦存狀態(tài)和高水平原巖應(yīng)力環(huán)境等角度介紹了沖擊地壓發(fā)生時(shí)斷層構(gòu)造的地質(zhì)賦存環(huán)境等特征,建立了考慮工作面開采時(shí)的斷層結(jié)構(gòu)力學(xué)模型和厚層頂板賦存的物理模型,推導(dǎo)了斷層面上下盤相對(duì)滑動(dòng)位移的逐點(diǎn)積分計(jì)算方法,提出了可表征斷層滑動(dòng)摩擦與能量釋放相關(guān)性的斷層結(jié)構(gòu)勢(shì)能的概念,分析了開采擾動(dòng)下斷層構(gòu)造位移場(chǎng)、應(yīng)力場(chǎng)和能量場(chǎng)的動(dòng)態(tài)演化特征,系統(tǒng)闡述了斷層構(gòu)造失穩(wěn)的發(fā)生特點(diǎn)和誘發(fā)沖擊地壓的力學(xué)機(jī)理。筆者的研究結(jié)果期望能為我國煤礦沖擊地壓的預(yù)測(cè)與防治提供理論支撐。

    1 斷層構(gòu)造失穩(wěn)誘沖的關(guān)鍵科學(xué)問題

    1.1 開采擾動(dòng)下斷層滑移失穩(wěn)的力學(xué)機(jī)理

    由于巖體受到地殼運(yùn)動(dòng)的強(qiáng)烈擠壓作用,使得在斷層等地質(zhì)構(gòu)造成型期間積聚了大量的彈性應(yīng)變能,是其誘發(fā)煤礦沖擊地壓的主要?jiǎng)恿Νh(huán)境。在現(xiàn)場(chǎng)監(jiān)測(cè)中,盡管可以利用鉆孔應(yīng)力、地質(zhì)雷達(dá)探測(cè)、微震監(jiān)測(cè)和電磁輻射監(jiān)測(cè)等技術(shù)手段獲得工作面開采過程中斷層結(jié)構(gòu)附近及上覆巖層的位移場(chǎng)、應(yīng)力場(chǎng)和能量場(chǎng)的宏觀變化特征。然而,目前能夠表征斷層面多物理場(chǎng)時(shí)空演化特征的力學(xué)機(jī)理解釋尚未清晰,無法準(zhǔn)確獲知斷層結(jié)構(gòu)失穩(wěn)的多物理場(chǎng)動(dòng)態(tài)演化特征和前兆信息。

    筆者建立了考慮工作面開采時(shí)的斷層結(jié)構(gòu)力學(xué)模型和厚層頂板賦存的物理模型,推導(dǎo)了斷層面上下盤相對(duì)滑動(dòng)位移計(jì)算的逐點(diǎn)積分方法,提出了表征斷層滑動(dòng)摩擦與能量釋放相關(guān)性的斷層結(jié)構(gòu)勢(shì)能的概念,分析了開采擾動(dòng)下斷層構(gòu)造位移場(chǎng)、應(yīng)力場(chǎng)和能量場(chǎng)的動(dòng)態(tài)演化特征,總結(jié)出斷層滑移失穩(wěn)動(dòng)態(tài)過程中摩擦狀態(tài)發(fā)生轉(zhuǎn)變的臨界位置和應(yīng)力條件,系統(tǒng)闡述了斷層構(gòu)造失穩(wěn)的發(fā)生特點(diǎn)和誘發(fā)沖擊地壓的力學(xué)機(jī)理。

    1.2 斷層構(gòu)造失穩(wěn)誘沖的多物理場(chǎng)前兆特征

    采動(dòng)條件下斷層構(gòu)造及其影響區(qū)域的應(yīng)力重分布動(dòng)態(tài)演化過程十分復(fù)雜,斷層面的靜態(tài)滑移、整體劇烈滑動(dòng)以及斷層結(jié)構(gòu)位移場(chǎng)的動(dòng)態(tài)演化過程,是研究斷層失穩(wěn)誘發(fā)沖擊地壓的重要前兆信息,而該過程勢(shì)必伴隨斷層面上正應(yīng)力和剪應(yīng)力的動(dòng)態(tài)響應(yīng)。同時(shí),由于受到來自側(cè)向壓力和上覆巖層下沉壓力的共同作用,斷層滑動(dòng)是斷層面摩擦效應(yīng)的體現(xiàn)。在斷層未滑動(dòng)之前,斷層面由于未達(dá)到最大滑動(dòng)靜摩擦,其斷層區(qū)域圍巖彈性應(yīng)變能處于長期積聚而未釋放的狀態(tài),斷層失穩(wěn)后的能量釋放也是當(dāng)前微震事件數(shù)、聲發(fā)射事件數(shù)和電磁輻射強(qiáng)度值監(jiān)測(cè)的重要目標(biāo),也可為沖擊地壓事故的發(fā)生做出預(yù)測(cè)。

    筆者在研究開采擾動(dòng)下斷層滑移失穩(wěn)的力學(xué)機(jī)理基礎(chǔ)上,研究斷層附近上覆巖層運(yùn)移與斷層滑動(dòng)破壞的相互作用,分析斷層側(cè)壓系數(shù)對(duì)斷層面應(yīng)力場(chǎng)的影響,給出斷層滑動(dòng)危險(xiǎn)性判別條件,著重分析斷層面位移場(chǎng)緩急突變、應(yīng)力場(chǎng)正剪分化以及能量場(chǎng)激增滯后的動(dòng)態(tài)演化特征,分析斷層結(jié)構(gòu)及影響區(qū)域多物理場(chǎng)異常演化現(xiàn)象與斷層滑移失穩(wěn)誘沖機(jī)理之間的聯(lián)系,從而獲得能夠預(yù)測(cè)斷層滑移失穩(wěn)的前兆特征。

    2 斷層構(gòu)造的地質(zhì)特征及誘沖模式

    斷層構(gòu)造的成型過程、地質(zhì)賦存狀態(tài),是斷層影響區(qū)域高水平原巖應(yīng)力環(huán)境的主要原因,了解斷層構(gòu)造的地質(zhì)特征有助于探究其失穩(wěn)誘發(fā)沖擊地壓的機(jī)理。

    2.1 斷層構(gòu)造的成因及地質(zhì)賦存狀態(tài)

    經(jīng)過長期的地殼運(yùn)動(dòng)和強(qiáng)烈的擠壓作用,巖層中較為薄弱的位置發(fā)生斷裂并與相鄰巖體產(chǎn)生相對(duì)運(yùn)動(dòng),從而形成斷層構(gòu)造。根據(jù)斷層構(gòu)造形成時(shí)所受到的應(yīng)力狀態(tài)以及巖石物理性質(zhì)的差異,可將其分為正斷層、逆斷層和平移斷層。斷層構(gòu)造成型后沿其走向上不同位置的巖體變形能力和能量積聚程度大不相同,在水平構(gòu)造應(yīng)力作用下主斷層兩側(cè)巖體上衍生出大量的節(jié)理和次生斷裂面,并與主斷層面形成具有一定寬度的斷層帶影響區(qū)域。然而,斷層等地質(zhì)構(gòu)造由于其深部賦存的特點(diǎn),在地層中并不容易被發(fā)現(xiàn)。目前,國內(nèi)外主要應(yīng)用地球物理探測(cè)技術(shù)對(duì)煤層賦存條件和復(fù)雜地質(zhì)構(gòu)造進(jìn)行勘探。比較常見的探測(cè)方法有地震法和電磁法,例如,彭蘇萍等將三維地震勘探技術(shù)應(yīng)用于淮南煤礦,并證明可以有效地探測(cè)小斷層、陷落柱、地質(zhì)異常體和活動(dòng)斷層等地質(zhì)構(gòu)造,成為煤礦高效安全生產(chǎn)的重要技術(shù)手段之一。王揚(yáng)州等將瞬變電磁法應(yīng)用到礦井中的全空間范圍,實(shí)現(xiàn)井下采空區(qū)多角度立體探測(cè)技術(shù)。因此,預(yù)先探測(cè)煤層的賦存環(huán)境和斷層等復(fù)雜地質(zhì)構(gòu)造的幾何形態(tài),有助于探究采動(dòng)影響下斷層等復(fù)雜地質(zhì)構(gòu)造附近多物理場(chǎng)的時(shí)空演化特征,為深入探索斷層滑移的力學(xué)機(jī)制提供了技術(shù)前提條件。

    斷層等復(fù)雜地質(zhì)構(gòu)造在我國各地區(qū)煤田中普遍存在。筆者以河南義馬煤業(yè)集團(tuán)股份有限公司義馬礦區(qū)、北京昊華能源股份有限公司京西礦區(qū)、黑龍江龍煤礦業(yè)集團(tuán)股份有限公司雞西礦區(qū)以及開灤(集團(tuán))唐山礦業(yè)分公司唐山煤礦作為工程案例,分析斷層等復(fù)雜地質(zhì)構(gòu)造的成因及其賦存地質(zhì)環(huán)境特征。

    圖1為河南義馬礦區(qū)地質(zhì)構(gòu)造。義馬礦區(qū)地處褶皺-逆沖斷層的構(gòu)造區(qū)域范圍內(nèi),形成具有大量走向、傾向、斜交的斷層和褶曲等復(fù)雜地質(zhì)構(gòu)造,煤層出現(xiàn)嚴(yán)重的分叉與合并現(xiàn)象。義馬礦區(qū)從西向東依次分布著楊村煤礦、耿村煤礦、千秋煤礦、躍進(jìn)煤礦和常村煤礦,向南以F陜石—義馬逆斷層為界線(圖2)。據(jù)現(xiàn)場(chǎng)資料統(tǒng)計(jì),千秋煤礦和躍進(jìn)煤礦中多個(gè)工作面的多次沖擊地壓事故均與F逆沖斷層有直接關(guān)系。

    圖1 義馬礦區(qū)地質(zhì)構(gòu)造

    圖2 義馬礦區(qū)各煤礦地理位置及2006—2015年各煤礦沖擊地壓發(fā)生頻率

    由于地殼運(yùn)動(dòng)過程中比較軟弱的煤層受擠壓應(yīng)力作用易發(fā)生斷裂和相對(duì)滑動(dòng),且在巨厚礫巖層的共同作用下形成壓扭性逆沖斷層F。斷層走向近東西,傾向南,略偏東,淺部傾角75°,深部傾角15°~35°,逆沖面上陡下緩,落差50~500 m,水平錯(cuò)距120~1 080 m,北接千秋煤礦,向東延入躍進(jìn)煤礦。根據(jù)圖3可知,F(xiàn)斷層緊鄰多個(gè)采煤工作面,其存在不僅切斷了巖層的連續(xù)性,還伴有應(yīng)力集中和能量積聚現(xiàn)象,增大了開采的難度。因此,斷層構(gòu)造的存在被認(rèn)為是該礦區(qū)沖擊地壓發(fā)生的主要誘因。

    圖3 義馬礦區(qū)千秋煤礦21221工作面采掘平面示意

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

    圖4 京西礦區(qū)地質(zhì)構(gòu)造示意

    圖5 京西礦區(qū)大安山煤礦和木城澗煤礦構(gòu)造剖面示意

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

    圖6 雞西礦區(qū)中部礦井開拓示意

    開灤(集團(tuán))唐山礦業(yè)分公司的唐山煤礦位于河北省唐山市,主要構(gòu)造為平行于地層走向的斷層,如圖7所示。由于受近似東西向的擠壓作用,井田西部賦存著復(fù)雜的褶皺構(gòu)造和多條縱向斷層,并伴隨著大量的中心斷層的衍生,數(shù)條壓扭性斷裂在深部合為陡河大斷裂。由于斷裂區(qū)域內(nèi)構(gòu)造應(yīng)力相對(duì)較大,從而成為唐山地區(qū)的主要的發(fā)震構(gòu)造。統(tǒng)計(jì)分析發(fā)現(xiàn),井田向斜軸部是發(fā)生沖擊地壓的高危險(xiǎn)區(qū),且煤巖沖擊失穩(wěn)與斷層構(gòu)造密切相關(guān)。

    圖7 唐山煤礦地質(zhì)構(gòu)造

    綜上所述,地殼運(yùn)動(dòng)和水平構(gòu)造擠壓力有助于斷層構(gòu)造等復(fù)雜地質(zhì)構(gòu)造的形成和發(fā)育,并且褶皺構(gòu)造通常會(huì)造成地層倒轉(zhuǎn)進(jìn)而會(huì)衍生出大量逆沖斷層,導(dǎo)致地層原巖應(yīng)力場(chǎng)變得更加復(fù)雜,在開采擾動(dòng)下容易發(fā)生斷層滑移失穩(wěn)從而誘發(fā)嚴(yán)重的沖擊破壞事故。因此,斷層等復(fù)雜地質(zhì)構(gòu)造對(duì)多物理場(chǎng)的影響是不容忽視的。

    2.2 高水平原巖應(yīng)力環(huán)境

    原巖應(yīng)力場(chǎng)是地殼巖石圈應(yīng)力的原始分布狀態(tài)。隨礦井開采深度的增加,原巖應(yīng)力也隨之增大,導(dǎo)致礦山壓力顯現(xiàn)強(qiáng)烈,為能量積聚提供了理想的內(nèi)部環(huán)境,在增加巷道維護(hù)和煤層開采難度的同時(shí),沖擊地壓發(fā)生的可能性也隨之增大。因此,考察原巖應(yīng)力場(chǎng)的分布特征是研究復(fù)雜地質(zhì)構(gòu)造特征與誘發(fā)沖擊動(dòng)力失穩(wěn)之間相關(guān)性的前提條件。

    高水平原巖應(yīng)力環(huán)境是我國煤礦復(fù)雜地質(zhì)賦存環(huán)境下誘發(fā)沖擊地壓的主要原因之一。以河南義馬礦區(qū)和京西礦區(qū)為例,根據(jù)原巖應(yīng)力現(xiàn)場(chǎng)測(cè)試結(jié)果可知,河南義馬礦區(qū)內(nèi)千秋煤礦是以水平應(yīng)力為主導(dǎo)的構(gòu)造型原巖應(yīng)力區(qū),而躍進(jìn)煤礦原巖應(yīng)力場(chǎng)則是由垂直應(yīng)力占主導(dǎo)??紤]斷層構(gòu)造的存在,水平構(gòu)造應(yīng)力依舊是義馬礦區(qū)高應(yīng)力環(huán)境中的主要應(yīng)力分量。與華北地區(qū)和中國大陸平均應(yīng)力水平相比而言,義馬煤田原巖應(yīng)力水平較高,如圖8所示。尤其是垂直應(yīng)力分別高于華北地區(qū)和全國中國大陸的42.7%和42.1%。以上數(shù)據(jù)進(jìn)一步說明了義馬礦區(qū)復(fù)雜地質(zhì)構(gòu)造條件下的高水平原巖應(yīng)力特征,也是頻繁發(fā)生沖擊地壓事故的重要原因。

    圖8 義馬礦區(qū)原巖應(yīng)力與華北地區(qū)和中國大陸原巖應(yīng)力對(duì)比

    京西礦區(qū)各礦井的開采水平均屬于深部開采范疇,原巖應(yīng)力測(cè)試結(jié)果可知該礦區(qū)原巖應(yīng)力場(chǎng)同樣以水平應(yīng)力為主導(dǎo)。如圖9所示,京西礦區(qū)中主應(yīng)力隨深度增加呈現(xiàn)不同程度的線性增長趨勢(shì)。由于礦區(qū)內(nèi)存在大量的倒轉(zhuǎn)型褶皺構(gòu)造和大中型逆斷層等復(fù)雜地質(zhì)構(gòu)造,使得原巖應(yīng)力水平要高于北京和華北等地區(qū)的平均水平,以深度800 m為例,其中三大主應(yīng)力分別高于中國大陸地區(qū)平均水平的34.1%,75.3%和57.2%。因此,在現(xiàn)有開采水平條件下,高水平原巖應(yīng)力釋放勢(shì)必會(huì)引起斷層等構(gòu)造的瞬間失穩(wěn),從而誘發(fā)沖擊地壓等動(dòng)力災(zāi)害。

    圖9 京西礦區(qū)原巖應(yīng)力和北京地區(qū)、華北地區(qū)和中國大陸原巖應(yīng)力平均水平對(duì)比

    所以,義馬礦區(qū)和京西礦區(qū)的高水平應(yīng)力環(huán)境,以及礦區(qū)內(nèi)頻繁發(fā)生的沖擊地壓事故,證實(shí)了高水平原巖應(yīng)力場(chǎng)與沖擊地壓之間的相關(guān)性,為沖擊地壓機(jī)理的研究提供了支撐。

    2.3 斷層沖擊地壓誘發(fā)模式

    斷層沖擊地壓作為一種與地質(zhì)構(gòu)造賦存狀態(tài)相關(guān)的沖擊類型,其發(fā)生機(jī)理是目前沖擊地壓研究領(lǐng)域的熱點(diǎn)。為了更加清晰地了解斷層構(gòu)造失穩(wěn)誘發(fā)沖擊地壓的力學(xué)機(jī)理,學(xué)者們提出了與斷層相關(guān)的誘沖模式來揭示沖擊災(zāi)害發(fā)生的原因。

    齊慶新等根據(jù)力源條件將斷層沖擊地壓歸納為構(gòu)造型沖擊地壓,認(rèn)為斷層等地質(zhì)構(gòu)造的存在使得地層中構(gòu)造應(yīng)力大幅提高,造成煤巖突然失穩(wěn)和沖擊破壞。潘一山等根據(jù)煤巖失穩(wěn)的物理特征對(duì)沖擊地壓機(jī)理進(jìn)行了分類,總結(jié)出斷層錯(cuò)動(dòng)型沖擊地壓一般發(fā)生在斷層構(gòu)造附近,并且只有在深部開采和上覆巖層壓力足夠大時(shí)才會(huì)引起斷層圍巖剪切失穩(wěn)從而誘發(fā)沖擊地壓。姜耀東等從應(yīng)力狀態(tài)導(dǎo)致煤巖突然破壞的本質(zhì)出發(fā),根據(jù)現(xiàn)有的開采工藝和常見地質(zhì)構(gòu)造賦存狀態(tài),將由斷層構(gòu)造或其他結(jié)構(gòu)面的滑移錯(cuò)動(dòng)引起的劇烈破壞的動(dòng)力現(xiàn)象稱為滑移錯(cuò)動(dòng)型沖擊地壓。譚云亮等提出深部開采斷層滑移型沖擊地壓,認(rèn)為深部斷層影響區(qū)域受采掘活動(dòng)影響更為顯著且構(gòu)造應(yīng)力明顯大于淺部開采。姜福興等依據(jù)微震數(shù)據(jù)監(jiān)測(cè)結(jié)果和構(gòu)造附近能量釋放情況,將復(fù)雜地質(zhì)構(gòu)造引起的沖擊地壓分為增壓型和減壓型,當(dāng)工作面臨近斷層構(gòu)造開采期間,微震事件在斷層附近出現(xiàn)積聚現(xiàn)象并且不隨工作面遷移時(shí),很可能引起構(gòu)造增壓型沖擊地壓在斷層附近發(fā)生,造成煤巖破壞。

    斷層構(gòu)造的滑移失穩(wěn)與采掘活動(dòng)和高水平應(yīng)力是密切相關(guān)的。竇林名等根據(jù)煤層開采活動(dòng)引起的斷層活化類型針對(duì)斷層型沖擊地壓提出了4種概念模型,即遠(yuǎn)場(chǎng)和近場(chǎng)礦震動(dòng)載觸發(fā)的斷層活化、逼近斷層和遠(yuǎn)離斷層開采的采動(dòng)應(yīng)力主導(dǎo)型斷層活化。呂進(jìn)國等根據(jù)采掘工作面與斷層構(gòu)造的相對(duì)空間位置不同引起的應(yīng)力場(chǎng)變化的差異性,將斷層沖擊地壓分為構(gòu)造應(yīng)力型、斷層活化型和斷層失穩(wěn)滑動(dòng)型等。

    筆者在研究斷層構(gòu)造位移場(chǎng)、應(yīng)力場(chǎng)和能量場(chǎng)動(dòng)態(tài)演化特征的基礎(chǔ)上,也將重點(diǎn)探討斷層構(gòu)造和上覆厚層頂板耦合賦存時(shí)的工作面逼近斷層構(gòu)造的采動(dòng)應(yīng)力主導(dǎo)型斷層沖擊地壓的力學(xué)機(jī)理,研究煤巖受采動(dòng)應(yīng)力和上覆巖層垂直壓力共同作用下引起的斷層失穩(wěn)滑動(dòng)型沖擊地壓。

    3 采動(dòng)下斷層構(gòu)造的多物理場(chǎng)特征

    探究斷層構(gòu)造滑移失穩(wěn)區(qū)域在開采擾動(dòng)過程中位移場(chǎng)、應(yīng)力場(chǎng)和能量場(chǎng)的動(dòng)態(tài)演化特征,建立斷層影響區(qū)域沖擊地壓發(fā)生前、后多物理場(chǎng)的互饋關(guān)系,得到能夠預(yù)測(cè)斷層失穩(wěn)誘發(fā)沖擊地壓的多物理場(chǎng)前兆特征,研究斷層結(jié)構(gòu)失穩(wěn)誘發(fā)沖擊地壓力學(xué)機(jī)理,從而達(dá)到準(zhǔn)確監(jiān)測(cè)以及有效預(yù)防沖擊地壓的目的。

    3.1 開采擾動(dòng)下斷層影響區(qū)位移場(chǎng)演化特征

    斷層影響區(qū)上覆巖層運(yùn)移特征

    為了能夠多參量監(jiān)測(cè)開采擾動(dòng)下斷層影響區(qū)域的位移場(chǎng)演化特征,以河南義馬礦區(qū)千秋煤礦21221工作面為工程地質(zhì)背景,在中國礦業(yè)大學(xué)(北京)深部巖土力學(xué)與地下工程國家重點(diǎn)實(shí)驗(yàn)室的深部工程破壞過程模型試驗(yàn)系統(tǒng)上開展了考慮工作面開采工況的斷層滑動(dòng)物理相似模擬試驗(yàn),如圖10所示。該物理模型的幾何相似比為120,容重相似比為0.67,模型尺寸為1 600 mm×1 600 mm×400 mm。

    圖10 斷層滑動(dòng)物理試驗(yàn)?zāi)P秃捅O(jiān)測(cè)方案

    研究顯示,除復(fù)雜地質(zhì)構(gòu)造賦存條件和高水平原巖應(yīng)力環(huán)境外,上覆巖層失穩(wěn)也是沖擊地壓發(fā)生的重要原因之一。隨著工作面推進(jìn),上覆巖層的空間結(jié)構(gòu)發(fā)生變化,造成大面積的覆巖運(yùn)動(dòng)甚至出現(xiàn)斷裂現(xiàn)象。因此,考察巖層運(yùn)移與斷層結(jié)構(gòu)滑移失穩(wěn)相關(guān)性對(duì)研究斷層沖擊地壓發(fā)生機(jī)理具有重要意義。

    圖11給出了上覆巖層位移隨工作面向斷層推進(jìn)過程中的變化規(guī)律,圖11中表示試驗(yàn)?zāi)P椭泄ぷ髅娴綌鄬拥木嚯x。從物理性質(zhì)上,堅(jiān)硬巖層的厚度和強(qiáng)度遠(yuǎn)大于煤層,上覆巖層通常被視為不易發(fā)生斷裂破壞的堅(jiān)硬頂板。隨著煤層開采過程中,上覆巖層在很長一段時(shí)間內(nèi)處于懸空狀態(tài)。上覆巖層發(fā)生破裂前,采空區(qū)圍巖變形影響范圍持續(xù)擴(kuò)大,巖層位移量也隨之增加,靠近斷層的巖層運(yùn)動(dòng)相對(duì)于遠(yuǎn)處更為活躍,并且在開采擾動(dòng)條件下發(fā)生多次不規(guī)律的周期性垮落,導(dǎo)致斷層頻繁受到擾動(dòng)。

    圖11 工作面開采過程中上覆巖層位移變化

    此外,研究還發(fā)現(xiàn)在開采條件下,斷層并非一直處于滑動(dòng)狀態(tài),只有當(dāng)上覆巖層發(fā)生垮落破壞時(shí)斷層才出現(xiàn)嚴(yán)重的滑動(dòng)。究其原因是由于上覆巖層未破壞前斷層結(jié)構(gòu)處于相對(duì)穩(wěn)定的狀態(tài),而在采動(dòng)應(yīng)力和巖層運(yùn)移共同作用下,斷層結(jié)構(gòu)附近出現(xiàn)了部分離層空間,使得施加到斷層面上的水平荷載減小,造成在斷層構(gòu)造一側(cè)呈現(xiàn)出卸載效應(yīng),加劇了斷層的活化程度。同時(shí),斷層滑動(dòng)又會(huì)反過來以非穩(wěn)態(tài)的沖擊和加卸載反作用于煤巖,使得上覆巖層發(fā)生更嚴(yán)重的垮落破壞,最終誘發(fā)沖擊地壓事故。以上結(jié)果與分析初步給出了巖層運(yùn)動(dòng)和斷層滑移的相互作用的關(guān)系。

    斷層結(jié)構(gòu)的相對(duì)滑動(dòng)位移的計(jì)算

    斷層滑移一般表現(xiàn)為上盤與下盤的相對(duì)錯(cuò)動(dòng),不僅會(huì)引起周圍圍巖運(yùn)動(dòng)的改變,還會(huì)誘發(fā)嚴(yán)重沖擊地壓事故。因此,進(jìn)一步研究在開采擾動(dòng)影響下斷層面上下盤之間的相對(duì)滑移過程,分析發(fā)生最危險(xiǎn)滑動(dòng)的位置及其演化特征,能夠更加全面地了解斷層的滑移機(jī)制。

    在煤礦開采過程中,斷層周圍煤巖受上覆巖層影響呈現(xiàn)整體向下移動(dòng)的趨勢(shì),隨著開采工作面逐漸靠近斷層,下沉位移量逐漸增大,最終引起斷層構(gòu)造上、下盤發(fā)生錯(cuò)動(dòng)。圖12給出了30 cm寬斷層帶區(qū)域的滑移情況。研究發(fā)現(xiàn),斷層構(gòu)造的滑動(dòng)位移最先出現(xiàn)在斷層頂部位置,并逐漸向下擴(kuò)展到煤層附近,而斷層最大滑移量最終在斷層與厚層頂板交界處的中部區(qū)域集中發(fā)展。從圖13可以發(fā)現(xiàn),斷層面上的位移隨著開采過程中由緩慢增加轉(zhuǎn)為突然急劇增加的現(xiàn)象,斷層面最大位移量也由斷層頂部突然轉(zhuǎn)移至斷層中部區(qū)域,造成斷層結(jié)構(gòu)滑移失穩(wěn)破壞。

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

    圖13 工作面推進(jìn)過程中斷面上斷層位移變化

    為進(jìn)一步解釋斷層面滑動(dòng)位移突變的現(xiàn)象,筆者提出了基于逐點(diǎn)積分的相對(duì)滑移算法來求解斷層面上盤和下盤之間的相對(duì)位移。斷層附近區(qū)域上盤與下盤在水平和垂直方向上的相對(duì)位移矢量如圖14所示。

    圖14 斷層面上下盤相對(duì)位移計(jì)算示意

    (1)

    (2)

    計(jì)算表明,開采擾動(dòng)下斷層面的相對(duì)滑動(dòng)是由局部靜態(tài)滑移向整體劇烈滑動(dòng)突然轉(zhuǎn)變的非線性動(dòng)態(tài)演化過程。在這個(gè)動(dòng)態(tài)演變過程中,斷層附近的圍巖運(yùn)動(dòng)經(jīng)歷了沿垂直向下、近似垂直于斷層和平行于斷層3個(gè)運(yùn)動(dòng)階段,如圖15所示。圖16給出了工作面在開采過程中斷層上下盤相對(duì)滑動(dòng)的動(dòng)態(tài)演化過程。在開采前期階段,斷層上下盤由于覆巖壓力作用整體向下移動(dòng)。隨著工作面的推進(jìn),斷層面上開始出現(xiàn)局部滑動(dòng),且斷層下盤相對(duì)于上盤沿?cái)鄬觾A向向下運(yùn)動(dòng),此時(shí)斷層面上不同位置的相對(duì)位移呈現(xiàn)基本同步的緩慢增長趨勢(shì),如圖17所示。臨近斷層開采時(shí),斷層上下盤的擠壓效應(yīng)代表了此時(shí)斷層結(jié)構(gòu)能量處于積聚狀態(tài),斷層中部區(qū)域的相對(duì)位移量大幅度增加,預(yù)示斷層發(fā)生嚴(yán)重的整體滑移破壞并釋放大量的應(yīng)變能。

    圖15 斷層面上下盤相對(duì)滑移矢量

    圖16 工作面開采過程中斷層相對(duì)滑動(dòng)動(dòng)態(tài)演化

    圖17 工作面開采過程中斷層下盤相對(duì)上盤的位移曲線

    綜上可知,斷層面最大相對(duì)位移量出現(xiàn)在斷層中部靠近煤層的區(qū)域,而相對(duì)位移緩急突變發(fā)生時(shí)伴隨著上覆巖層嚴(yán)重垮落現(xiàn)象。上覆巖層運(yùn)移失穩(wěn)會(huì)瞬時(shí)釋放大量應(yīng)變能,從而促進(jìn)了厚層頂板與斷層構(gòu)造交界區(qū)域的相對(duì)滑動(dòng),特別是未破壞巖層的失穩(wěn)破壞可能會(huì)進(jìn)一步引起嚴(yán)重的斷層滑移失穩(wěn)。筆者從斷層上下盤相對(duì)位移角度入手揭示了開采擾動(dòng)下斷層滑移由局部滑移轉(zhuǎn)向整體失穩(wěn)的動(dòng)態(tài)過程,證實(shí)了斷層滑移失穩(wěn)對(duì)上覆巖層運(yùn)移的影響作用。

    3.2 斷層結(jié)構(gòu)的力學(xué)模型及應(yīng)力場(chǎng)演化特征

    斷層面應(yīng)力分布的監(jiān)測(cè)結(jié)果

    采動(dòng)應(yīng)力的存在打破了采場(chǎng)原位應(yīng)力的平衡。工作面支承壓力與斷層處高構(gòu)造應(yīng)力疊加形成應(yīng)力集中,導(dǎo)致斷層滑移失穩(wěn)。

    因此,進(jìn)一步探究采動(dòng)影響下斷層構(gòu)造應(yīng)力場(chǎng)的演化特征,能夠更加深入地了解研究斷層滑移失穩(wěn)的力學(xué)機(jī)制。

    圖18為部分監(jiān)測(cè)點(diǎn)應(yīng)力隨工作面開采變化曲線(監(jiān)測(cè)點(diǎn)位置如圖10所示)。在開采初期階段,工作面距離斷層較遠(yuǎn)時(shí)(距斷層處100 cm),斷層面上的應(yīng)力開始發(fā)生緩慢變化,斷層開始受到擾動(dòng),與上述斷層滑移情況基本一致,可將該位置視為斷層影響區(qū)和非影響區(qū)的分界線。

    圖18 工作面推進(jìn)過程中斷層面上正應(yīng)力和剪應(yīng)力變化

    臨近斷層開采時(shí)(距斷層處60 cm),斷層因開采活動(dòng)影響頻繁受到擾動(dòng),特別是斷層中部區(qū)域的應(yīng)力變化尤其劇烈。當(dāng)工作面在斷層附近開采時(shí)(距斷層處20 cm),斷層發(fā)生嚴(yán)重的滑移失穩(wěn),此時(shí)應(yīng)力水平明顯降低。

    開采擾動(dòng)影響下,斷層面正應(yīng)力和剪應(yīng)力從同步變化轉(zhuǎn)為突然分化的急劇變化。斷層面上正應(yīng)力變化相比剪應(yīng)力而言更為強(qiáng)烈,且近煤層處的應(yīng)力值也要大于遠(yuǎn)離煤層處。除外,在臨近斷層開采時(shí),斷層面中部區(qū)域正應(yīng)力和剪應(yīng)力呈現(xiàn)截然相反的變化趨勢(shì),即正應(yīng)力增加而剪應(yīng)力降低或剪應(yīng)力增加而正應(yīng)力降低的分化現(xiàn)象。

    考慮水平構(gòu)造應(yīng)力和采動(dòng)應(yīng)力的影響,斷層構(gòu)造在側(cè)向壓力和上覆巖層下沉壓力的共同作用下而處于壓剪受力狀態(tài)。斷層面上的摩擦力隨著正應(yīng)力增加而增大,剪應(yīng)力和正應(yīng)力持續(xù)增大導(dǎo)致斷層結(jié)構(gòu)由靜態(tài)摩擦轉(zhuǎn)為動(dòng)態(tài)摩擦,此時(shí)斷層面上剪應(yīng)力會(huì)相應(yīng)減小。加之?dāng)鄬又胁繀^(qū)域巖層被壓實(shí)且不易松,進(jìn)而出現(xiàn)正應(yīng)力增加而剪應(yīng)力降低的分化現(xiàn)象。為進(jìn)一步證實(shí)這一監(jiān)測(cè)結(jié)果,筆者提出了考慮采動(dòng)影響的斷層滑動(dòng)的力學(xué)模型,解釋了正應(yīng)力和剪應(yīng)力的分化特征,并做詳細(xì)闡述。

    考慮采動(dòng)影響的斷層結(jié)構(gòu)的力學(xué)模型

    由于斷層面上應(yīng)力場(chǎng)分布和變化較為復(fù)雜,僅依靠現(xiàn)場(chǎng)監(jiān)測(cè)和數(shù)值仿真技術(shù)還不能完全揭示開采擾動(dòng)下斷層面上的應(yīng)力場(chǎng)演化特征。為進(jìn)一步驗(yàn)證和解釋斷層面上正應(yīng)力和剪應(yīng)力的分化現(xiàn)象,筆者建立了考慮采動(dòng)影響的斷層表面應(yīng)力分布和斷層滑動(dòng)特征的力學(xué)模型,如圖19所示。綜合考慮上覆巖層壓力、重力附加應(yīng)力和采動(dòng)應(yīng)力的影響,得到式(3)給出的斷層面上任意一點(diǎn)處的應(yīng)力分量,進(jìn)而計(jì)算得到如式(4)所示的斷層面上正應(yīng)力和剪應(yīng)力的表達(dá)式。

    (3)

    (4)

    其中,,分別為沿方向的正應(yīng)力分量、沿方向的正應(yīng)力分量和剪切應(yīng)力分量;為斷層面任意點(diǎn)的正應(yīng)力和剪應(yīng)力;和分別為任一點(diǎn)的水平和鉛直距離;和分別為模型的長度和高度;為采動(dòng)應(yīng)力峰值點(diǎn)到斷層的水平距離;為垂直分布應(yīng)力;為巖石密度;為重力加速度;為斷層傾角;()和()分別為重力附加應(yīng)力和采動(dòng)應(yīng)力分布的水平應(yīng)力函數(shù);()d為距坐標(biāo)原點(diǎn)處微小長度d上所受的微小集中力。

    圖19中,為斷層面上任意點(diǎn)到模型頂部的鉛直距離;′為力學(xué)模型中工作面與斷層之間的距離;為側(cè)向系數(shù)。

    圖19 考慮采動(dòng)影響的斷層面上任意點(diǎn)應(yīng)力計(jì)算的力學(xué)模型

    斷層面應(yīng)力場(chǎng)的動(dòng)態(tài)演化特征

    在采動(dòng)應(yīng)力影響下,斷層面上不同位置應(yīng)力分布各異,圖20給出了隨工作面向斷層推進(jìn)過程中斷層面應(yīng)力值的分布和變化規(guī)律。圖21為斷層面應(yīng)力計(jì)算的試驗(yàn)校對(duì)。對(duì)比試驗(yàn)與理論結(jié)果可知,關(guān)于斷層面上正應(yīng)力和剪應(yīng)力的分布特征,理論模型計(jì)算結(jié)果和試驗(yàn)結(jié)果基本吻合。結(jié)果顯示,斷層面下段區(qū)域應(yīng)力水平高于斷層面上段區(qū)域,隨著工作面的推進(jìn),應(yīng)力變化幅度更加明顯。沿?cái)鄬用嫦蛳?,斷層面的正?yīng)力值和剪應(yīng)力值變化呈相反趨勢(shì),出現(xiàn)正應(yīng)力增加而剪應(yīng)力降低或剪應(yīng)力增加而正應(yīng)力降低的分化現(xiàn)象。應(yīng)力場(chǎng)變化的轉(zhuǎn)折點(diǎn)代表著斷層表面的最小正應(yīng)力值和最大剪應(yīng)力值,說明了應(yīng)力極值與斷層滑移失穩(wěn)具有一定的聯(lián)系。

    圖20 斷層面應(yīng)力分布的理論模型計(jì)算結(jié)果

    圖21 斷層面應(yīng)力分布的試驗(yàn)結(jié)果校對(duì)

    根據(jù)正應(yīng)力和剪應(yīng)力的極值變化,可將斷層表面劃分為R,R和R三個(gè)區(qū)域,如圖22所示。工作面在到達(dá)斷層擾動(dòng)影響區(qū)繼續(xù)向前推進(jìn)時(shí),斷層面的應(yīng)力極值逐漸向底部移動(dòng)。沿?cái)鄬用嫦蛳?,根?jù)正應(yīng)力和剪應(yīng)力的變化狀態(tài),得出R為斷層上、下盤的靜摩擦穩(wěn)定狀態(tài)區(qū)域,R區(qū)域中剪應(yīng)力極值位置處于即將動(dòng)態(tài)滑移的臨界平衡狀態(tài),在R區(qū)域斷層為已經(jīng)發(fā)生了動(dòng)態(tài)相對(duì)滑動(dòng)。因此,斷層滑動(dòng)面的剪應(yīng)力極值出現(xiàn)的位置可視為斷層面從靜態(tài)轉(zhuǎn)向動(dòng)態(tài)滑動(dòng)的臨界區(qū)域,同時(shí)也是斷層面上最先發(fā)生滑動(dòng)的位置。

    圖22 開采過程中斷層面上的應(yīng)力演化

    綜上可知,當(dāng)工作面臨近斷層開采時(shí),由于采動(dòng)應(yīng)力和上覆巖層的影響,斷層結(jié)構(gòu)面和厚層頂板的交界區(qū)域的應(yīng)力水平較高,此時(shí)該區(qū)域積累了大量的應(yīng)變能,為斷層滑移失穩(wěn)提供了應(yīng)力條件。此外,上述結(jié)果還給出了應(yīng)力極值分布演化特征,并從摩擦狀態(tài)轉(zhuǎn)變的角度分析了斷層滑動(dòng)的力學(xué)機(jī)制。

    3.3 斷層失穩(wěn)的能量場(chǎng)時(shí)空演化特征

    斷層結(jié)構(gòu)勢(shì)能概念的提出

    開采擾動(dòng)下斷層應(yīng)力的變化會(huì)引起斷層表面應(yīng)變能的積累和釋放。斷層周圍煤巖能量變化實(shí)質(zhì)上就是從穩(wěn)定積聚狀態(tài)逐漸向非穩(wěn)定釋放狀態(tài)轉(zhuǎn)化的過程。由于較高應(yīng)力水平、斷層構(gòu)造以及開采條件等諸多影響因素,造成開采過程中煤巖的能量變化具有明顯的時(shí)空特征。在現(xiàn)場(chǎng)監(jiān)測(cè)中,微震事件數(shù)、聲發(fā)射事件數(shù)和電磁輻射強(qiáng)度值等數(shù)據(jù)在沖擊地壓發(fā)生前常出現(xiàn)監(jiān)測(cè)數(shù)據(jù)消失的現(xiàn)象。為進(jìn)一步解釋這一工程現(xiàn)象,在斷層面應(yīng)力分布的理論模型基礎(chǔ)上提出了斷層結(jié)構(gòu)勢(shì)能的概念,并對(duì)斷層結(jié)構(gòu)勢(shì)能的時(shí)空演化特征進(jìn)行了分析。

    根據(jù)彈性理論中的應(yīng)力變分方程,通過斷層面上的應(yīng)力來計(jì)算斷層結(jié)構(gòu)勢(shì)能為

    (5)

    斷層結(jié)構(gòu)勢(shì)能演化特征

    圖23為斷層面上結(jié)構(gòu)勢(shì)能的計(jì)算結(jié)果。工作面過斷層前,斷層結(jié)構(gòu)勢(shì)能分布特征與剪應(yīng)力分布特征相似,對(duì)照?qǐng)D22可知,結(jié)構(gòu)勢(shì)能最大值與剪應(yīng)力最大值類似,均出現(xiàn)在R區(qū)域的底部附近,并且在R區(qū)域逐漸減小。當(dāng)采煤工作面穿過斷層后,斷層結(jié)構(gòu)勢(shì)能變化明顯,斷層頂部結(jié)構(gòu)勢(shì)能突然增大,R區(qū)域的結(jié)構(gòu)勢(shì)能急劇下降。從圖24可以看出,隨著工作面向斷層推進(jìn),斷層面上部區(qū)域的結(jié)構(gòu)勢(shì)能逐漸減小,而中下部區(qū)域結(jié)構(gòu)勢(shì)能則呈增大的趨勢(shì)。

    圖23 煤層開采過程中斷層面上結(jié)構(gòu)勢(shì)能的分布

    圖24 工作面開采過程中斷層面上結(jié)構(gòu)勢(shì)能演化

    綜上可知,開采工作面過斷層前結(jié)構(gòu)勢(shì)能最大的位置出現(xiàn)在斷層結(jié)構(gòu)面和上覆厚層頂板的交界區(qū)域,同時(shí)也是最大剪應(yīng)力出現(xiàn)的位置,即斷層上下盤運(yùn)動(dòng)由靜摩擦轉(zhuǎn)變?yōu)閯?dòng)摩擦的臨界區(qū)域,可將該位置認(rèn)為是斷層滑動(dòng)的起始點(diǎn)。由于斷層滑移錯(cuò)動(dòng)導(dǎo)致斷層表面應(yīng)力對(duì)開采活動(dòng)產(chǎn)生不同的響應(yīng),從應(yīng)變能角度解釋了正應(yīng)力和剪應(yīng)力變化過程中的分化特征。

    由于斷層區(qū)域的上覆巖層較為堅(jiān)硬不易破壞,能量激增前斷層處于相對(duì)穩(wěn)定的狀態(tài),導(dǎo)致監(jiān)測(cè)數(shù)據(jù)缺失現(xiàn)象,但此時(shí)能量則一直處于積聚而未能釋放的狀態(tài)。在開采擾動(dòng)和斷層滑移作用下,未破壞上覆巖層發(fā)生失穩(wěn)導(dǎo)致斷層頂部的能量積聚和瞬時(shí)釋放,從而引起斷層的嚴(yán)重滑動(dòng)破壞誘發(fā)沖擊地壓。筆者從結(jié)構(gòu)勢(shì)能角度解釋了斷層滑移失穩(wěn)過程中能量的穩(wěn)態(tài)積聚和非穩(wěn)態(tài)釋放的變化,同時(shí)證實(shí)了上覆巖層運(yùn)移和破裂對(duì)斷層滑移失穩(wěn)具有重要的貢獻(xiàn)。

    4 斷層結(jié)構(gòu)失穩(wěn)誘發(fā)沖擊地壓機(jī)理

    開采擾動(dòng)下斷層結(jié)構(gòu)發(fā)生滑移錯(cuò)動(dòng)引起斷層圍巖位移場(chǎng)、應(yīng)力場(chǎng)和能量場(chǎng)發(fā)生改變,繼而誘發(fā)嚴(yán)重的沖擊地壓。筆者在研究斷層結(jié)構(gòu)滑移失穩(wěn)多物理場(chǎng)演化特征的基礎(chǔ)上,進(jìn)一步分析采動(dòng)影響下斷層沖擊地壓的力學(xué)機(jī)理。

    4.1 斷層結(jié)構(gòu)失穩(wěn)誘發(fā)沖擊地壓的前兆特征

    斷層失穩(wěn)誘沖的位移場(chǎng)突變特征

    由前文可知,在整個(gè)煤層開采過程中,斷層的滑動(dòng)位移通常會(huì)經(jīng)歷靜止、緩慢滑移和快速滑動(dòng)3個(gè)階段。斷層最大位移出現(xiàn)位置也發(fā)生了由斷層頂部轉(zhuǎn)向中部的空間變化。在臨近斷層開采時(shí),斷層位移量突然明顯增加,代表斷層發(fā)生嚴(yán)重滑移失穩(wěn),可將位移場(chǎng)這種由穩(wěn)定到急劇增加的突變現(xiàn)象作為斷層滑移失穩(wěn)的前兆特征。

    斷層面上位移場(chǎng)的突變特征還表現(xiàn)在采動(dòng)應(yīng)變和應(yīng)變率的監(jiān)測(cè)結(jié)果中,圖25和圖26分別為斷層面應(yīng)變和應(yīng)變率的變化曲線。當(dāng)斷層開始擾動(dòng)后,斷層附近圍巖應(yīng)變出現(xiàn)急劇上升和穩(wěn)定增加的相互交替變化階段,隨后斷層發(fā)生嚴(yán)重滑移,采動(dòng)應(yīng)變?cè)俅渭眲p小。同時(shí),斷層附近應(yīng)變率也由低水平穩(wěn)定狀態(tài)轉(zhuǎn)變?yōu)榧眲〔▌?dòng)階段,預(yù)示上覆巖層中裂縫急劇發(fā)展,最終引起上覆巖層垮落破壞。因此,可將斷層面上采動(dòng)應(yīng)變急劇增加和維持穩(wěn)定交替出現(xiàn)的階段或應(yīng)變率急劇變化的現(xiàn)象視為斷層滑移的另一前兆特征。

    圖25 煤層開采過程中斷層面應(yīng)變變化

    圖26 斷層面上應(yīng)變率變化曲線以及試驗(yàn)對(duì)比

    斷層失穩(wěn)誘沖的應(yīng)力場(chǎng)分化特征

    斷層面上中部區(qū)域的應(yīng)力場(chǎng)正剪分化特征表現(xiàn)為正應(yīng)力增加而剪應(yīng)力降低或剪應(yīng)力增加而正應(yīng)力降低的現(xiàn)象。為進(jìn)一步分析斷層面應(yīng)力場(chǎng)這一特征,圖27給出了斷層中部應(yīng)力場(chǎng)變化以及對(duì)應(yīng)的位移場(chǎng)變化。正應(yīng)力和剪應(yīng)力的共同作用使得斷層面發(fā)生滑移失穩(wěn),此時(shí)由于水平應(yīng)力的擠壓作用,正應(yīng)力持續(xù)增大,代表能量逐漸積聚,而剪應(yīng)力隨斷層滑移錯(cuò)動(dòng)開始呈下降趨勢(shì)。另外,當(dāng)斷層出現(xiàn)滑動(dòng)時(shí)伴隨著正應(yīng)力的急劇下降和剪應(yīng)力的突然增加,此時(shí)能量開始釋放,正應(yīng)力的減小會(huì)導(dǎo)致斷層表面剪切強(qiáng)度的降低,從而加速了嚴(yán)重的斷層滑移。隨著工作面向斷層方向推進(jìn),斷層相對(duì)滑移量逐漸增加,最大位移量出現(xiàn)的位置也開始向下部轉(zhuǎn)移,直到斷層滑移失穩(wěn)造成斷層區(qū)域圍巖的破壞,此時(shí)正應(yīng)力和剪應(yīng)力均處于較小值。因此,可以將應(yīng)力場(chǎng)的正剪分化現(xiàn)象視為斷層滑移的前兆特征。

    圖27 工作面回采期間斷層中部區(qū)域正應(yīng)力和剪應(yīng)力變化及相對(duì)應(yīng)的斷層滑移情況

    為進(jìn)一步解釋斷層中部發(fā)生劇烈相對(duì)滑動(dòng)的現(xiàn)象,根據(jù)探地雷達(dá)監(jiān)測(cè)到的物理模型前、后兩側(cè)的電磁波幅值差值,得出開采過程中斷層中部區(qū)域振幅衰減變化,如圖28所示。煤層開采過程中,正應(yīng)力和剪應(yīng)力開始發(fā)生分化現(xiàn)象時(shí),振幅衰減量從極小值開始增長并隨之保持穩(wěn)定。在這個(gè)過程中,斷層附近圍巖中有大量裂隙產(chǎn)生并擴(kuò)展,正應(yīng)力和剪應(yīng)力的急劇變化同時(shí)也說明了斷層由局部滑動(dòng)演化為整體滑動(dòng)。以上從巖石損傷破壞程度的角度驗(yàn)證了應(yīng)力場(chǎng)正剪分化現(xiàn)象作為斷層滑移前兆信息的正確性。

    圖28 開采過程中斷層面正應(yīng)力、剪應(yīng)力曲線及相應(yīng)位置的振幅衰減變化

    斷層失穩(wěn)誘沖的能量場(chǎng)激增特征

    聲發(fā)射技術(shù)能夠監(jiān)測(cè)到能量場(chǎng)在煤層開采過程中的釋放情況,從而判斷斷層滑移破壞程度。

    如圖29所示,聲發(fā)射事件數(shù)在工作面推進(jìn)過程中出現(xiàn)多次激增,分別是斷層加載初期、臨近距斷層開采和到達(dá)斷層處開采時(shí)。斷層結(jié)構(gòu)釋放的能量主要聚集在斷層與上覆厚層頂板交界處的中部區(qū)域,這與斷層最大滑動(dòng)位移和剪應(yīng)力最大值出現(xiàn)位置近似,同斷層結(jié)構(gòu)勢(shì)能計(jì)算結(jié)果基本一致。

    圖29 工作面開采事件數(shù)分布

    此外,聲發(fā)射事件在激增前會(huì)出現(xiàn)事件缺失的現(xiàn)象,與工程監(jiān)測(cè)結(jié)果基本一致。以河南義馬礦區(qū)的重大沖擊地壓案為例。圖30給出了2011年義馬礦區(qū)千秋煤礦和2010年躍進(jìn)煤礦的電磁輻射強(qiáng)度和支架阻力的監(jiān)測(cè)結(jié)果。電磁輻射強(qiáng)度和支架阻力值呈現(xiàn)基本一致的增大與減小循環(huán)變化特征,表明沖擊地壓發(fā)生前斷層和上覆巖層積聚了大量變形能,而電磁輻射強(qiáng)度和支架阻力突降說明大量應(yīng)變能開始釋放,證實(shí)了斷層結(jié)構(gòu)滑移失穩(wěn)時(shí)能量場(chǎng)的激增現(xiàn)象。圖31為2011年義馬礦區(qū)千秋煤礦監(jiān)測(cè)到的微震事件頻次與能量釋放隨時(shí)間變化的現(xiàn)場(chǎng)監(jiān)測(cè)結(jié)果。監(jiān)測(cè)結(jié)果中存在微震事件的缺失現(xiàn)象,并且對(duì)應(yīng)時(shí)間段的微震事件頻次也會(huì)相對(duì)較低,此時(shí)能量處于積聚而未釋放的狀態(tài),與本文斷層結(jié)構(gòu)勢(shì)能計(jì)算結(jié)果基本一致。

    圖30 千秋煤礦和躍進(jìn)煤礦電磁輻射強(qiáng)度和支架阻力變化曲線

    圖31 千秋煤礦微震能量及頻次累計(jì)情況

    斷層滑移失穩(wěn)過程中能量場(chǎng)變化與應(yīng)力場(chǎng)、位移場(chǎng)緊密相關(guān),圖32和33分別給出了斷層面聲發(fā)射事件與位移場(chǎng)、應(yīng)力場(chǎng)之間的變化關(guān)系。研究表明,斷層面上的滑移并非時(shí)時(shí)刻刻在發(fā)生的,只有在滑移較為嚴(yán)重并釋放大量能量時(shí)才可捕捉到聲發(fā)射事件。而上覆巖層突然垮落或斷層發(fā)生嚴(yán)重滑移之前,巖層中裂縫較少,聲發(fā)射事件也相對(duì)穩(wěn)定,此時(shí)能量場(chǎng)處于一直積聚的狀態(tài)。因此,根據(jù)位移場(chǎng)緩急突變和能量場(chǎng)激增現(xiàn)象先后出現(xiàn),可將聲發(fā)射事件急劇增加前的恒定不變階段視為斷層滑移的前兆信息。

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

    圖33 斷層面上剪應(yīng)力和聲發(fā)射事件數(shù)之間的聯(lián)系

    根據(jù)斷層面聲發(fā)射事件與應(yīng)力場(chǎng)的變化關(guān)系可知,當(dāng)工作面開采至斷層附近時(shí),斷層聲發(fā)射事件數(shù)激增,此時(shí)斷層面剪應(yīng)力出現(xiàn)逐步降低后陡然增加的現(xiàn)象,斷層面即將發(fā)生滑動(dòng)。當(dāng)工作面開采至斷層處,應(yīng)力降至最小值,聲發(fā)射事件數(shù)再次激增,斷層滑移失穩(wěn)并釋放出大量的能量。因此,斷層面應(yīng)力在逐步降低的過程中陡然增加和聲發(fā)射事件數(shù)恒定不變時(shí)突然激增的特征可作為斷層滑移失穩(wěn)的前兆信息。此外,從上述結(jié)果中可以發(fā)現(xiàn),能量場(chǎng)的激增現(xiàn)象一般滯后于位移場(chǎng)的突變現(xiàn)象和應(yīng)力場(chǎng)的分化現(xiàn)象。

    綜上可知,斷層的滑動(dòng)失穩(wěn)過程中位移場(chǎng)的緩急突變,應(yīng)力場(chǎng)的正剪分化以及能量場(chǎng)的激增滯后出現(xiàn)的位置基本一致,斷層面最大位移量代表斷層出現(xiàn)嚴(yán)重滑移破壞,應(yīng)力場(chǎng)正剪分化現(xiàn)象是斷層的滑移錯(cuò)動(dòng)造成的異常響應(yīng),并且伴隨著大量積聚的應(yīng)變能瞬間釋放,以上結(jié)果均表明了斷層結(jié)構(gòu)在開采擾動(dòng)影響下多物理場(chǎng)的互饋關(guān)系。

    4.2 斷層結(jié)構(gòu)失穩(wěn)誘發(fā)沖擊地壓的力學(xué)機(jī)理

    斷層滑移失穩(wěn)滑動(dòng)危險(xiǎn)性判別

    斷層滑移失穩(wěn)滑動(dòng)危險(xiǎn)性主要與斷層構(gòu)造有關(guān),尤其是斷層構(gòu)造的幾何形態(tài)對(duì)滑移失穩(wěn)具有一定的影響作用。比如,斷層切向剛度主要控制著斷層剪切變形滑移,隨著斷層切向剛度的增加,剪應(yīng)力值變化越明顯,斷層滑移危險(xiǎn)性越高。落差與斷層正應(yīng)力呈正相關(guān)關(guān)系,落差越大,正應(yīng)力越大,斷層更易發(fā)生滑移失穩(wěn)。另外,斷層傾角不僅代表了斷層的傾斜程度,描述了斷層的大致形態(tài)和走向,還是影響斷層滑移失穩(wěn)的一個(gè)重要指標(biāo)。

    為考察斷層構(gòu)造的傾角與斷層滑移失穩(wěn)的關(guān)系,圖34給出了不同斷層傾角作用下位移場(chǎng)的變化。研究發(fā)現(xiàn),斷層相對(duì)位移在開采后期突變隨傾角增大而越發(fā)明顯,尤其是傾角為75°時(shí)。開采初期釋放能量較少代表了應(yīng)變能在持續(xù)積聚,此時(shí)斷層結(jié)構(gòu)面相對(duì)位移幾乎為0,并且隨著工作面推進(jìn)呈現(xiàn)緩慢增長趨勢(shì)。當(dāng)工作面臨近斷層開采時(shí),大量應(yīng)變能的突然釋放導(dǎo)致斷層面相對(duì)位移的急劇增加,最終造成嚴(yán)重的斷層滑移失穩(wěn)。

    圖34 不同斷層傾角模型的斷層附近上覆巖層相對(duì)位移變化

    為進(jìn)一步解釋采動(dòng)過程中傾角變化對(duì)斷層結(jié)構(gòu)相對(duì)位移的影響,提出了斷層滑移危險(xiǎn)性的判別條件。圖35給出了斷層面不同位置相對(duì)位移的變化情況。斷層面底部區(qū)域相對(duì)位移值較大,且斷層的相對(duì)位移隨傾角變化逐漸減小甚至趨近于0,表明斷層上、下盤由于強(qiáng)擠壓作用下彈性應(yīng)變能處于長期積聚而未釋放的狀態(tài)。斷層的傾角越大,出現(xiàn)擠壓應(yīng)變積聚的區(qū)域范圍越大,斷層滑移失穩(wěn)時(shí)瞬時(shí)釋放的能量越多,誘發(fā)沖擊地壓的危險(xiǎn)性越高。因此,當(dāng)斷層傾角較大時(shí),來自水平側(cè)向壓力和上覆巖層下沉壓力作用越顯著,斷層結(jié)構(gòu)處應(yīng)力集中和能量積聚越明顯,斷層滑動(dòng)更加劇烈,發(fā)生沖擊地壓的危險(xiǎn)性越高。

    圖35 煤層開采過程中不同斷層傾角斷層附近上覆巖層的相對(duì)位移

    高應(yīng)力環(huán)境下斷層構(gòu)造失穩(wěn)誘沖機(jī)理

    復(fù)雜多變的地質(zhì)因素對(duì)斷層面應(yīng)力場(chǎng)演化和滑移量突變具有重要的影響作用,尤其是水平構(gòu)造應(yīng)力對(duì)斷層的擠壓作用。在水平構(gòu)造應(yīng)力作用下,隨著采空區(qū)范圍增大,斷層上下盤的相對(duì)滑動(dòng)在一定程度上穩(wěn)步釋放了水平構(gòu)造應(yīng)力擠壓積聚的應(yīng)變能。圖36為斷層滑移量和滑移速率(每計(jì)算100步作為單位時(shí)間)隨工作面向斷層方向推進(jìn)的變化情況,發(fā)現(xiàn)在臨近斷層面開采時(shí),滑移速率絕對(duì)值突增和滑移量陡增,此時(shí)斷層滑移危險(xiǎn)極高。

    圖36 斷層滑移量及滑移速率

    為進(jìn)一步研究水平構(gòu)造應(yīng)力對(duì)斷層滑移失穩(wěn)的影響,圖37給出了不同水平側(cè)壓系數(shù)作用下應(yīng)力場(chǎng)的演化特征。當(dāng)側(cè)壓系數(shù)>1.0時(shí),斷層下段的正應(yīng)力值較大,剪應(yīng)力極大值也處于較高水平,此時(shí)上盤與下盤之間擠壓力明顯,為逆斷層的形成創(chuàng)造了壓力條件。對(duì)于側(cè)壓系數(shù)<1.0的情況,正應(yīng)力在斷層上段隨側(cè)壓系數(shù)的減小而增大,下段應(yīng)力水平相對(duì)較低,有利于正斷層的形成。因此,筆者初步給出了正、逆斷層形成的力學(xué)基本條件,但對(duì)于斷層形成的力學(xué)機(jī)制還需進(jìn)一步的研究和推導(dǎo)。

    圖37 不同側(cè)壓系數(shù)下斷層面正應(yīng)力和剪應(yīng)力的變化

    斷層構(gòu)造和上覆厚層頂板耦合誘沖機(jī)理

    斷層構(gòu)造和上覆厚層頂板賦存的復(fù)雜地質(zhì)環(huán)境中,在采動(dòng)應(yīng)力影響作用下,巖層運(yùn)移和斷層滑移相互作用增加了沖擊地壓發(fā)生的可能性。由于堅(jiān)硬頂板的厚度和強(qiáng)度相較于煤層很大,當(dāng)巖層發(fā)生下沉或垮落現(xiàn)象,采場(chǎng)發(fā)生沖擊破壞的危險(xiǎn)性更高。

    根據(jù)前文可知,上覆巖層運(yùn)移與斷層滑移失穩(wěn)兩者存在相互作用的關(guān)系。上覆巖層擠壓斷層產(chǎn)生滑動(dòng),而斷層的滑移又作用于上覆巖層致使其大范圍運(yùn)動(dòng)甚至垮落,二者的聯(lián)合作用很可能會(huì)誘發(fā)沖擊地壓。圖38給出了義馬礦區(qū)沖擊地壓成因示意。逆沖斷層作為地質(zhì)環(huán)境中的最不穩(wěn)定的因素,不僅反轉(zhuǎn)了煤層方向,并且在其成型期間積累了巨大的應(yīng)變能,在開采擾動(dòng)下斷層積聚的應(yīng)變能會(huì)在短時(shí)間內(nèi)突然釋放。由于上覆厚層頂板對(duì)煤層產(chǎn)生高強(qiáng)度的動(dòng)壓力,導(dǎo)致能量場(chǎng)在上覆巖層與斷層構(gòu)造的交界處出現(xiàn)激增滯后的異常現(xiàn)象。因此,煤田高水平原巖應(yīng)力為上覆巖層運(yùn)移和斷層的活化創(chuàng)造了有利條件,綜合以上多因素耦合作用誘發(fā)沖擊動(dòng)力災(zāi)害。

    圖38 義馬礦區(qū)沖擊地壓成因示意

    此外,根據(jù)探地雷達(dá)監(jiān)測(cè)到未破壞礫巖的低波速帶可知,如圖39所示。由于未受破壞巖層可以保護(hù)斷層不受開采擾動(dòng)和頂板垮落破壞的影響,斷層在開采前期存在穩(wěn)定階段。當(dāng)巖層發(fā)生沉降和分離時(shí)釋放了大量的應(yīng)變能,從而加速了未破壞巖層與斷層構(gòu)造相互影響區(qū)域的相對(duì)滑動(dòng)。當(dāng)上覆巖層坍塌或失穩(wěn)時(shí),反而會(huì)進(jìn)一步引起更加劇烈的斷層相對(duì)滑動(dòng)導(dǎo)致沖擊地壓的發(fā)生。上述研究結(jié)果說明,在復(fù)雜斷層構(gòu)造賦存、高水平原巖應(yīng)力環(huán)境和開采擾動(dòng)條件下,上覆巖層運(yùn)動(dòng)與斷層滑移的耦合作用是斷層發(fā)生滑移失穩(wěn)的重要原因。

    圖39 煤層開采過程中雷達(dá)探測(cè)結(jié)果及巖層位移矢量

    5 結(jié)論與展望

    (1)在開采擾動(dòng)下,斷層面相對(duì)滑動(dòng)位移最早出現(xiàn)在斷層頂部并逐漸向下發(fā)展,且由緩慢滑動(dòng)變化為突然急劇滑移的非線性動(dòng)態(tài)失穩(wěn)過程,最后在斷層結(jié)構(gòu)面和上覆厚層頂板交界的中部區(qū)域發(fā)生最大相對(duì)滑動(dòng)位移。斷層面上不同位置應(yīng)力分布各異,且正應(yīng)力變化相比剪應(yīng)力更為劇烈,尤其是正應(yīng)力和剪應(yīng)力在同步演化發(fā)展過程中發(fā)生了分化的現(xiàn)象,并表現(xiàn)出正應(yīng)力增加而剪應(yīng)力降低或剪應(yīng)力增加而正應(yīng)力降低的特征。斷層失穩(wěn)前,結(jié)構(gòu)勢(shì)能處于不斷積聚而未釋放的狀態(tài),斷層失穩(wěn)后,結(jié)構(gòu)勢(shì)能出現(xiàn)由幾乎為零轉(zhuǎn)為突然激增的現(xiàn)象。

    (2)斷層結(jié)構(gòu)面和上覆厚層頂板的交界區(qū)域出現(xiàn)的最大滑動(dòng)靜摩擦向滑動(dòng)失穩(wěn)摩擦的轉(zhuǎn)換,是斷層面位移場(chǎng)突變和應(yīng)力場(chǎng)分化的時(shí)刻,也是斷層結(jié)構(gòu)勢(shì)能從零到激增的時(shí)刻。因此,斷層面位移場(chǎng)緩急突變、應(yīng)力場(chǎng)正剪分化以及能量場(chǎng)激增滯后的現(xiàn)象均可視為斷層滑移失穩(wěn)的前兆特征。

    (3)在復(fù)雜斷層構(gòu)造賦存、高水平原巖應(yīng)力環(huán)境以及開采擾動(dòng)條件下,由于斷層附近厚層頂板下沉運(yùn)動(dòng)和垮落破壞的影響,加劇了斷層的活化和滑移失穩(wěn)程度,導(dǎo)致斷層結(jié)構(gòu)面和厚層頂板的交界區(qū)域出現(xiàn)了由滑動(dòng)靜摩擦轉(zhuǎn)為劇烈動(dòng)態(tài)摩擦的現(xiàn)象。斷層最大位移量的出現(xiàn)預(yù)示著斷層的嚴(yán)重失穩(wěn),釋放該區(qū)域積聚的大量應(yīng)變能,從而誘發(fā)沖擊地壓災(zāi)害。

    (4)該研究是對(duì)煤礦復(fù)雜地質(zhì)構(gòu)造賦存環(huán)境下沖擊地壓發(fā)生機(jī)理和斷層構(gòu)造失穩(wěn)的多物理場(chǎng)前兆特征的研究探索。斷層失穩(wěn)過程中,斷層滑動(dòng)面及周邊圍巖的幾何結(jié)構(gòu)和溫度場(chǎng)也應(yīng)有相應(yīng)的變化,本文的成果為這一方向研究打下理論基礎(chǔ),同時(shí)為沖擊地壓的預(yù)測(cè)和防治提供理論支撐。

    猜你喜歡
    剪應(yīng)力應(yīng)力場(chǎng)巖層
    高應(yīng)力巖層巷道鉆孔爆破卸壓技術(shù)
    變截面波形鋼腹板組合箱梁的剪應(yīng)力計(jì)算分析
    地球故事之復(fù)理石巖層
    鋁合金多層多道窄間隙TIG焊接頭應(yīng)力場(chǎng)研究
    焊接(2016年9期)2016-02-27 13:05:22
    考慮斷裂破碎帶的丹江口庫區(qū)地應(yīng)力場(chǎng)與水壓應(yīng)力場(chǎng)耦合反演及地震預(yù)測(cè)
    基于位移相關(guān)法的重復(fù)壓裂裂縫尖端應(yīng)力場(chǎng)研究
    斷塊油氣田(2014年5期)2014-03-11 15:33:49
    三噴兩錨一注漿+U型鋼聯(lián)合支護(hù)在松軟巖層中的應(yīng)用
    河南科技(2014年22期)2014-02-27 14:18:11
    瀝青路面最大剪應(yīng)力分析
    河南科技(2014年13期)2014-02-27 14:11:25
    岸坡應(yīng)力場(chǎng)及卸荷帶劃分量化指標(biāo)研究
    復(fù)合式路面層間最大剪應(yīng)力影響因素研究
    色视频在线一区二区三区| 国产成人一区二区在线| 国产男女超爽视频在线观看| 少妇丰满av| 搡老乐熟女国产| 精品一区二区三区视频在线| 国产成人a区在线观看| 不卡视频在线观看欧美| 国产亚洲5aaaaa淫片| 国国产精品蜜臀av免费| 少妇 在线观看| 男人添女人高潮全过程视频| 国产高潮美女av| 久久久久久伊人网av| 免费观看无遮挡的男女| 成年av动漫网址| 日本熟妇午夜| 国内揄拍国产精品人妻在线| 69av精品久久久久久| 婷婷色av中文字幕| 一个人看视频在线观看www免费| 国产亚洲5aaaaa淫片| 亚洲精品国产成人久久av| 免费少妇av软件| 伦理电影大哥的女人| 亚洲天堂av无毛| 欧美性猛交╳xxx乱大交人| 亚洲国产色片| 小蜜桃在线观看免费完整版高清| 日韩人妻高清精品专区| 亚洲最大成人av| 国产精品一区www在线观看| 欧美一级a爱片免费观看看| 欧美日韩一区二区视频在线观看视频在线 | 少妇丰满av| 黄色怎么调成土黄色| a级毛色黄片| 免费看不卡的av| 婷婷色麻豆天堂久久| 色吧在线观看| av女优亚洲男人天堂| 国产日韩欧美在线精品| 免费观看在线日韩| 日本猛色少妇xxxxx猛交久久| 免费观看的影片在线观看| 极品少妇高潮喷水抽搐| 精品熟女少妇av免费看| 有码 亚洲区| 亚洲欧美中文字幕日韩二区| 欧美97在线视频| 寂寞人妻少妇视频99o| 亚洲av免费在线观看| 一区二区三区乱码不卡18| 亚洲精品一区蜜桃| 国产黄片美女视频| 爱豆传媒免费全集在线观看| 日韩 亚洲 欧美在线| 人妻 亚洲 视频| 久久久亚洲精品成人影院| 国产黄a三级三级三级人| 男人和女人高潮做爰伦理| 亚洲人成网站高清观看| 日韩人妻高清精品专区| 亚洲精品久久午夜乱码| 99久久人妻综合| 国产精品一区二区在线观看99| 又大又黄又爽视频免费| 国产亚洲精品久久久com| 国产精品麻豆人妻色哟哟久久| 亚洲国产高清在线一区二区三| 丰满人妻一区二区三区视频av| 亚洲va在线va天堂va国产| 国产综合懂色| 国产乱人视频| 搡老乐熟女国产| 色视频www国产| 91狼人影院| 欧美日韩视频精品一区| 91久久精品国产一区二区成人| 亚洲久久久久久中文字幕| 成年版毛片免费区| xxx大片免费视频| 好男人在线观看高清免费视频| 中文在线观看免费www的网站| 九九在线视频观看精品| 爱豆传媒免费全集在线观看| 久久人人爽人人片av| 欧美高清成人免费视频www| 久久久午夜欧美精品| av在线天堂中文字幕| 婷婷色av中文字幕| 一级毛片黄色毛片免费观看视频| 各种免费的搞黄视频| 亚洲av免费高清在线观看| 亚洲丝袜综合中文字幕| 少妇人妻一区二区三区视频| 午夜福利视频精品| 免费看av在线观看网站| 天天一区二区日本电影三级| 国产精品无大码| 一边亲一边摸免费视频| 男人狂女人下面高潮的视频| 人体艺术视频欧美日本| 一级毛片电影观看| 亚洲国产成人一精品久久久| 亚洲欧美精品自产自拍| 欧美日韩视频高清一区二区三区二| 亚洲国产成人一精品久久久| 国产精品一区二区性色av| 欧美国产精品一级二级三级 | 最近手机中文字幕大全| 另类亚洲欧美激情| 亚洲av成人精品一二三区| 欧美成人一区二区免费高清观看| 欧美潮喷喷水| 婷婷色综合大香蕉| 亚洲四区av| 亚洲丝袜综合中文字幕| 男女下面进入的视频免费午夜| 王馨瑶露胸无遮挡在线观看| 人人妻人人看人人澡| 国产精品无大码| 毛片一级片免费看久久久久| 日本一本二区三区精品| 深夜a级毛片| 国产成人a∨麻豆精品| av在线app专区| 麻豆久久精品国产亚洲av| 久久精品综合一区二区三区| 精品99又大又爽又粗少妇毛片| 国产高清不卡午夜福利| av一本久久久久| 国产精品久久久久久精品古装| 亚洲色图av天堂| 一级毛片 在线播放| 在线免费十八禁| 国产精品不卡视频一区二区| 久久99热这里只有精品18| 国产精品熟女久久久久浪| 亚洲av中文字字幕乱码综合| 国产一区二区在线观看日韩| 午夜福利视频1000在线观看| 国产伦理片在线播放av一区| 国产成人91sexporn| 在现免费观看毛片| 国产成人精品福利久久| 国产精品三级大全| 国产精品蜜桃在线观看| 91久久精品电影网| 一个人观看的视频www高清免费观看| 老师上课跳d突然被开到最大视频| 国产免费一级a男人的天堂| 亚洲精品国产成人久久av| 亚洲av不卡在线观看| 最近中文字幕2019免费版| 啦啦啦在线观看免费高清www| 麻豆乱淫一区二区| 在线天堂最新版资源| 欧美少妇被猛烈插入视频| 婷婷色av中文字幕| 欧美日韩一区二区视频在线观看视频在线 | av一本久久久久| 国产精品蜜桃在线观看| 国产爱豆传媒在线观看| 国产免费视频播放在线视频| 99久久精品一区二区三区| 日韩欧美精品v在线| 王馨瑶露胸无遮挡在线观看| 国产精品不卡视频一区二区| 亚洲精品乱久久久久久| 如何舔出高潮| 欧美日韩一区二区视频在线观看视频在线 | 少妇人妻精品综合一区二区| 五月开心婷婷网| 亚洲精品色激情综合| 久久久久性生活片| 国产成人一区二区在线| 亚洲精品aⅴ在线观看| 欧美亚洲 丝袜 人妻 在线| 麻豆乱淫一区二区| 国产色爽女视频免费观看| 中文天堂在线官网| 中文精品一卡2卡3卡4更新| 下体分泌物呈黄色| 80岁老熟妇乱子伦牲交| 交换朋友夫妻互换小说| 日本猛色少妇xxxxx猛交久久| 久久精品夜色国产| 赤兔流量卡办理| 91aial.com中文字幕在线观看| 性插视频无遮挡在线免费观看| 91精品国产九色| 国产av码专区亚洲av| 一边亲一边摸免费视频| 国产乱来视频区| 久久久久国产网址| 亚洲怡红院男人天堂| 一区二区三区乱码不卡18| 亚洲在线观看片| 亚洲av欧美aⅴ国产| 嫩草影院精品99| 成人鲁丝片一二三区免费| 九草在线视频观看| 亚洲国产av新网站| 欧美日韩国产mv在线观看视频 | 久久久久久久大尺度免费视频| 国产男女超爽视频在线观看| 国产成人免费观看mmmm| 性色avwww在线观看| 久久久a久久爽久久v久久| 久久精品人妻少妇| av在线天堂中文字幕| 亚洲精品国产av蜜桃| 精品国产露脸久久av麻豆| 七月丁香在线播放| 国产亚洲av嫩草精品影院| 老师上课跳d突然被开到最大视频| 成人亚洲欧美一区二区av| 视频区图区小说| 黄色怎么调成土黄色| 成人二区视频| kizo精华| 国产午夜精品久久久久久一区二区三区| 男女国产视频网站| 国产 精品1| 成人美女网站在线观看视频| 日本免费在线观看一区| 亚洲av男天堂| 日韩欧美精品免费久久| videos熟女内射| 成人亚洲欧美一区二区av| 精品午夜福利在线看| 日本欧美国产在线视频| 国产白丝娇喘喷水9色精品| 男女边吃奶边做爰视频| 国产高清三级在线| 亚洲伊人久久精品综合| 国产乱来视频区| 亚洲第一区二区三区不卡| 一级毛片久久久久久久久女| 日韩一区二区视频免费看| 日韩制服骚丝袜av| 亚洲国产成人一精品久久久| 亚洲精品视频女| 亚洲国产最新在线播放| 黄色配什么色好看| 18禁在线播放成人免费| 天堂网av新在线| 欧美bdsm另类| 欧美丝袜亚洲另类| av网站免费在线观看视频| 18禁在线无遮挡免费观看视频| 中文资源天堂在线| 午夜激情久久久久久久| 中文字幕人妻熟人妻熟丝袜美| 蜜臀久久99精品久久宅男| 婷婷色综合www| 久久精品国产亚洲av涩爱| 精品久久久久久久久av| 国产日韩欧美在线精品| 亚洲综合色惰| 亚洲精品国产av蜜桃| 国产黄色免费在线视频| 好男人在线观看高清免费视频| 白带黄色成豆腐渣| 亚洲欧美日韩卡通动漫| 日韩国内少妇激情av| 有码 亚洲区| 中文字幕人妻熟人妻熟丝袜美| 少妇人妻一区二区三区视频| 国产精品一二三区在线看| 国产精品一区二区性色av| 一个人观看的视频www高清免费观看| 欧美老熟妇乱子伦牲交| 国产一区二区亚洲精品在线观看| 精品一区二区三卡| 大码成人一级视频| videos熟女内射| 婷婷色综合大香蕉| 国产一级毛片在线| 日本与韩国留学比较| 老司机影院成人| h日本视频在线播放| 我要看日韩黄色一级片| 看黄色毛片网站| 五月玫瑰六月丁香| 蜜桃久久精品国产亚洲av| 国产精品福利在线免费观看| 久久精品国产亚洲网站| 亚洲激情五月婷婷啪啪| 熟女人妻精品中文字幕| 亚洲av欧美aⅴ国产| 春色校园在线视频观看| av天堂中文字幕网| 王馨瑶露胸无遮挡在线观看| 国产欧美亚洲国产| 亚洲,一卡二卡三卡| av在线亚洲专区| 久久久久久国产a免费观看| 涩涩av久久男人的天堂| 菩萨蛮人人尽说江南好唐韦庄| 欧美丝袜亚洲另类| 日本与韩国留学比较| 欧美xxxx性猛交bbbb| 欧美三级亚洲精品| 日韩免费高清中文字幕av| 久久久久久久久久久免费av| 午夜福利视频精品| h日本视频在线播放| 久久精品久久精品一区二区三区| 久久精品国产亚洲av涩爱| 国产精品爽爽va在线观看网站| 男人添女人高潮全过程视频| 最近最新中文字幕免费大全7| 久久久久久久大尺度免费视频| 免费看a级黄色片| 亚洲人与动物交配视频| 日韩av在线免费看完整版不卡| 日韩一区二区视频免费看| 各种免费的搞黄视频| 亚洲国产精品国产精品| 欧美bdsm另类| 99热这里只有精品一区| 亚洲熟女精品中文字幕| 特大巨黑吊av在线直播| 99热国产这里只有精品6| 国产精品伦人一区二区| 五月伊人婷婷丁香| 欧美极品一区二区三区四区| 久久精品国产自在天天线| 婷婷色av中文字幕| av免费观看日本| 色视频www国产| 久久久午夜欧美精品| 日韩 亚洲 欧美在线| 亚洲精品aⅴ在线观看| 国内精品宾馆在线| 国产免费福利视频在线观看| 男人和女人高潮做爰伦理| 高清午夜精品一区二区三区| 最近手机中文字幕大全| 永久网站在线| 老司机影院成人| 亚洲欧美一区二区三区黑人 | 国产亚洲午夜精品一区二区久久 | 在线观看美女被高潮喷水网站| 亚洲精品一二三| 国产亚洲最大av| 可以在线观看毛片的网站| 午夜激情久久久久久久| 99久久精品一区二区三区| 一边亲一边摸免费视频| 久久久精品免费免费高清| 高清av免费在线| 亚洲三级黄色毛片| 国产黄a三级三级三级人| 国产探花在线观看一区二区| 九色成人免费人妻av| 亚洲在线观看片| xxx大片免费视频| 免费少妇av软件| 精品午夜福利在线看| 成人美女网站在线观看视频| 2018国产大陆天天弄谢| 80岁老熟妇乱子伦牲交| 国产精品精品国产色婷婷| 另类亚洲欧美激情| 日韩强制内射视频| 久久久久久久久久成人| 一二三四中文在线观看免费高清| 亚洲色图av天堂| 天天躁日日操中文字幕| 精品国产露脸久久av麻豆| 国语对白做爰xxxⅹ性视频网站| 插逼视频在线观看| 性色avwww在线观看| 国产有黄有色有爽视频| 国产片特级美女逼逼视频| 亚洲精品aⅴ在线观看| 简卡轻食公司| 免费高清在线观看视频在线观看| 免费看av在线观看网站| 久久精品久久久久久噜噜老黄| 午夜精品国产一区二区电影 | 婷婷色av中文字幕| 少妇丰满av| 亚洲美女搞黄在线观看| 日韩大片免费观看网站| 欧美潮喷喷水| 成人漫画全彩无遮挡| 伦精品一区二区三区| 久久精品久久精品一区二区三区| 人人妻人人澡人人爽人人夜夜| 大片免费播放器 马上看| 亚洲av二区三区四区| 久久午夜福利片| 国产老妇女一区| 高清欧美精品videossex| 国内精品美女久久久久久| 色婷婷久久久亚洲欧美| 男人舔奶头视频| 日本黄色片子视频| 成人国产麻豆网| 97在线人人人人妻| 少妇人妻 视频| 午夜激情福利司机影院| 久久精品国产亚洲av涩爱| 少妇高潮的动态图| 嫩草影院精品99| 波多野结衣巨乳人妻| 国产av不卡久久| 老师上课跳d突然被开到最大视频| 亚洲精品,欧美精品| 在线免费十八禁| 国产高清不卡午夜福利| 一个人观看的视频www高清免费观看| 97超碰精品成人国产| 2018国产大陆天天弄谢| 亚洲精品一区蜜桃| 少妇丰满av| 国产精品久久久久久精品古装| 永久网站在线| 国精品久久久久久国模美| av专区在线播放| 视频区图区小说| 国产乱来视频区| 成人毛片a级毛片在线播放| 国产免费一区二区三区四区乱码| 91狼人影院| 国产一区二区亚洲精品在线观看| 中文字幕人妻熟人妻熟丝袜美| 日本色播在线视频| 久久久久国产精品人妻一区二区| 三级男女做爰猛烈吃奶摸视频| 日本av手机在线免费观看| 久久这里有精品视频免费| 日韩三级伦理在线观看| 亚洲图色成人| 日本午夜av视频| 国产精品国产三级国产专区5o| 边亲边吃奶的免费视频| 久久久久九九精品影院| 色综合色国产| 美女xxoo啪啪120秒动态图| 亚洲av中文字字幕乱码综合| 日韩三级伦理在线观看| 亚洲美女搞黄在线观看| 久久精品国产自在天天线| 美女高潮的动态| 午夜日本视频在线| 国产老妇伦熟女老妇高清| av天堂中文字幕网| 国产亚洲精品久久久com| 中文字幕av成人在线电影| 亚洲成人一二三区av| 最新中文字幕久久久久| 国产免费又黄又爽又色| 偷拍熟女少妇极品色| 国产黄片美女视频| 国产淫片久久久久久久久| 日日啪夜夜撸| 高清毛片免费看| 日韩三级伦理在线观看| 日韩成人伦理影院| 男男h啪啪无遮挡| 天堂中文最新版在线下载 | 中文资源天堂在线| 亚洲婷婷狠狠爱综合网| 国产精品一区二区性色av| 日韩一区二区三区影片| 草草在线视频免费看| 国产色婷婷99| 午夜亚洲福利在线播放| 国产久久久一区二区三区| 亚洲伊人久久精品综合| 80岁老熟妇乱子伦牲交| 91久久精品国产一区二区成人| 亚洲国产高清在线一区二区三| 51国产日韩欧美| 国产 一区精品| 精品人妻一区二区三区麻豆| 男女国产视频网站| 看十八女毛片水多多多| 日韩亚洲欧美综合| 深夜a级毛片| 中文字幕人妻熟人妻熟丝袜美| 男女边吃奶边做爰视频| 黄色一级大片看看| 男人爽女人下面视频在线观看| 街头女战士在线观看网站| 亚洲av在线观看美女高潮| 日本猛色少妇xxxxx猛交久久| 内射极品少妇av片p| videossex国产| 中文字幕av成人在线电影| 18禁动态无遮挡网站| 国产一区亚洲一区在线观看| 人妻一区二区av| 久久午夜福利片| 中文字幕av成人在线电影| 久久久久久久大尺度免费视频| 国产一区亚洲一区在线观看| 亚洲,欧美,日韩| 亚洲va在线va天堂va国产| 婷婷色麻豆天堂久久| 亚洲图色成人| 日本爱情动作片www.在线观看| 熟女电影av网| 国产精品一二三区在线看| 亚洲综合色惰| 97超碰精品成人国产| 18禁在线无遮挡免费观看视频| 国内少妇人妻偷人精品xxx网站| 欧美日韩亚洲高清精品| 一级片'在线观看视频| 人人妻人人爽人人添夜夜欢视频 | 国产精品人妻久久久影院| 亚洲国产精品999| 一区二区三区免费毛片| 亚洲精品视频女| 亚洲不卡免费看| 欧美国产精品一级二级三级 | 麻豆成人av视频| 婷婷色综合www| 国产精品成人在线| 亚洲人成网站高清观看| 校园人妻丝袜中文字幕| 三级国产精品欧美在线观看| 美女视频免费永久观看网站| 又粗又硬又长又爽又黄的视频| 天天一区二区日本电影三级| 精品一区二区三区视频在线| 看非洲黑人一级黄片| 成年女人在线观看亚洲视频 | av天堂中文字幕网| 久久精品国产亚洲网站| av免费在线看不卡| 伦精品一区二区三区| 2021少妇久久久久久久久久久| 亚洲欧美精品专区久久| 成人美女网站在线观看视频| av国产精品久久久久影院| 两个人的视频大全免费| 国产精品蜜桃在线观看| 亚洲欧洲日产国产| 亚洲经典国产精华液单| 免费在线观看成人毛片| 高清日韩中文字幕在线| 欧美变态另类bdsm刘玥| 亚洲高清免费不卡视频| 日韩av在线免费看完整版不卡| 国产在线一区二区三区精| 在线观看国产h片| 精品国产露脸久久av麻豆| 成年免费大片在线观看| 黑人高潮一二区| 十八禁网站网址无遮挡 | 婷婷色综合www| 深夜a级毛片| 大又大粗又爽又黄少妇毛片口| 亚洲aⅴ乱码一区二区在线播放| 天天躁夜夜躁狠狠久久av| 久久精品国产a三级三级三级| 男女国产视频网站| 国产熟女欧美一区二区| 超碰97精品在线观看| 又粗又硬又长又爽又黄的视频| 岛国毛片在线播放| 亚洲欧美一区二区三区国产| 久久久久久久久大av| 中国美白少妇内射xxxbb| 亚洲最大成人手机在线| 在线观看一区二区三区| 国产在线男女| 久久久久九九精品影院| 亚洲成人中文字幕在线播放| 国产伦精品一区二区三区四那| 一区二区三区免费毛片| 亚洲成人精品中文字幕电影| 成人美女网站在线观看视频| 日日摸夜夜添夜夜爱| 亚洲av中文av极速乱| 新久久久久国产一级毛片| 美女xxoo啪啪120秒动态图| 在线看a的网站| 丝瓜视频免费看黄片| 国产高清国产精品国产三级 | 国产男女内射视频| 欧美最新免费一区二区三区| 三级国产精品片| 国产欧美另类精品又又久久亚洲欧美| www.av在线官网国产| 最近的中文字幕免费完整| 成人国产av品久久久| 国产免费视频播放在线视频| 午夜福利在线在线| 国产高清三级在线| 一二三四中文在线观看免费高清| 汤姆久久久久久久影院中文字幕| 免费少妇av软件| 美女内射精品一级片tv| 一区二区三区精品91| 午夜精品一区二区三区免费看| 91精品伊人久久大香线蕉| 蜜臀久久99精品久久宅男| 国产色婷婷99| 一区二区三区免费毛片| 97人妻精品一区二区三区麻豆| 久热这里只有精品99| 一本久久精品| 欧美bdsm另类| av在线亚洲专区| 乱系列少妇在线播放| 亚洲av男天堂| 岛国毛片在线播放|