礦區(qū)沉陷預(yù)計(jì)程序與計(jì)算機(jī)仿真

前言：想要寫(xiě)出一篇引人入勝的文章？我們特意為您整理了礦區(qū)沉陷預(yù)計(jì)程序與計(jì)算機(jī)仿真范文，希望能給你帶來(lái)靈感和參考，敬請(qǐng)閱讀。

1礦區(qū)概況

1.1開(kāi)采區(qū)及建筑物分布情況

蒲白礦務(wù)局馬村礦南盤(pán)采區(qū)引起的開(kāi)采沉陷損害影響范圍較大，西固鎮(zhèn)的多數(shù)建筑物都受到不同程度的損害。開(kāi)采區(qū)煤系地層屬于石炭系山西組和太原群，包含5個(gè)煤層，5#煤層為該區(qū)主采煤層，基巖主要包括粗砂巖、細(xì)砂巖、砂頁(yè)巖互層、頁(yè)巖和泥質(zhì)砂巖。西固鎮(zhèn)約900戶、4000人和12個(gè)村落。全鎮(zhèn)的地面建筑物多為舊窯洞、瓦房、平房及少量樓房。受到開(kāi)采破壞的地面建筑物約有350間平房、300間樓房及部分磚拱窯洞[2]。

1.2擬定開(kāi)采方案

按照變條帶協(xié)調(diào)開(kāi)采的原則，結(jié)合開(kāi)采區(qū)域的煤層賦存條件、礦區(qū)地表建筑物的分布情況及煤柱的采留寬理論計(jì)算分析，擬定2種開(kāi)采方案：方案Ⅰ傾斜變條帶協(xié)調(diào)開(kāi)采；方案Ⅱ走向變條帶協(xié)調(diào)開(kāi)采。開(kāi)采參數(shù)如下：方案Ⅰ留寬a=25m，采寬b=18～32m，依次按2m的步距逐漸遞增；方案Ⅱ留煤柱寬度b=25m，采寬b=25m[2]。通過(guò)YHL-12預(yù)計(jì)軟件系統(tǒng)和YHL-12數(shù)值模擬計(jì)算分析，確定安全、合理、經(jīng)濟(jì)的開(kāi)采方案。

2地表移動(dòng)變形預(yù)計(jì)方法

2.1地表移動(dòng)變形預(yù)計(jì)數(shù)學(xué)模型

設(shè)計(jì)模擬計(jì)算模型時(shí)，為了提高預(yù)計(jì)計(jì)算準(zhǔn)確度，常使用多塊段劃分計(jì)算，可以保證滿足開(kāi)采方法、形狀和時(shí)間變化的要求，同時(shí)可以減小開(kāi)采厚度、高度、深度及開(kāi)采時(shí)間和順序因素所造成的影響。

2.2YLH-12預(yù)計(jì)評(píng)價(jià)系統(tǒng)介紹

目前，國(guó)內(nèi)外對(duì)礦區(qū)開(kāi)采引起的地表沉陷、破壞預(yù)計(jì)理論和方法做了大量的研究,取得了豐碩的成果。開(kāi)采沉陷涉及開(kāi)采方法、覆巖內(nèi)部結(jié)構(gòu)、巖性和破壞機(jī)理、礦區(qū)地質(zhì)構(gòu)造及受保護(hù)建筑設(shè)施的結(jié)構(gòu)等多重因素,是一個(gè)十分復(fù)雜的影響系統(tǒng)。YLH-12預(yù)計(jì)評(píng)價(jià)系統(tǒng)正是集地表位移變形預(yù)計(jì)、開(kāi)采方法優(yōu)選、建筑物與地表破壞程度圖形分區(qū)分類(lèi)評(píng)價(jià)及參數(shù)優(yōu)化選擇為一體的綜合評(píng)價(jià)系統(tǒng)。

2.3影響函數(shù)仿真模擬計(jì)算

2.3.1預(yù)計(jì)理論

西固村鎮(zhèn)開(kāi)采區(qū)域，濕陷性黃土覆蓋層厚度大，厚度大約在80m以上，屬于特厚黃土覆蓋層下厚煤層開(kāi)采類(lèi)型。煤層埋藏約為230m，地質(zhì)構(gòu)造較為復(fù)雜，開(kāi)采條件變化較大，開(kāi)采時(shí)間先后不同且開(kāi)采形狀不規(guī)則。為了提高預(yù)計(jì)計(jì)算的準(zhǔn)確度，減小開(kāi)采厚度、高度、深度及開(kāi)采時(shí)間和順序因素所造成的影響，在計(jì)算中采用多塊段劃分計(jì)算[5]，采用極坐標(biāo)閉合回路積分通用模型。

2.3.2確立預(yù)計(jì)參數(shù)

根據(jù)礦區(qū)條帶開(kāi)采的地表觀測(cè)數(shù)據(jù)，結(jié)合同類(lèi)采礦地質(zhì)條件下條帶開(kāi)采的觀測(cè)數(shù)據(jù)，以及南盤(pán)采區(qū)鉆孔柱狀巖性綜合分析，確立了概率積分法預(yù)計(jì)參數(shù)：水平移動(dòng)系數(shù)b=0.28～0.32；下沉系數(shù)η=0.89；最大下沉角φ=86°；主要影響角正切值tanβ=2.5，β=68.2°；拐點(diǎn)偏移距d=0.03H[6]。

2.3.3預(yù)計(jì)成果分析

采用預(yù)計(jì)評(píng)價(jià)軟件YLH-12，計(jì)算出2種方案開(kāi)采后地表的移動(dòng)變形值，再經(jīng)過(guò)圖像處理軟件制作出移動(dòng)變形等值線圖。根據(jù)受損地表建筑物涉及范圍內(nèi)的移動(dòng)變形值及建筑物的受破壞情況綜合分析，確定開(kāi)采損害對(duì)地表建筑物的影響，從而分析和判斷開(kāi)采方案的可行性，并對(duì)開(kāi)采方案進(jìn)行優(yōu)化設(shè)計(jì)，選擇最實(shí)用的開(kāi)采方案。采用變條帶協(xié)調(diào)開(kāi)采方案，地表下沉盆地呈偏態(tài)分布，位于采區(qū)邊界一側(cè)，下沉盆地拐點(diǎn)位置較正常情況下平緩得多，從而能夠控制和減小建筑物涉及范圍內(nèi)的地表移動(dòng)變形值，有效地保護(hù)了建筑物不受損害。方案Ⅰ中建筑物位置的移動(dòng)變形更小，較方案Ⅱ更適合西固鎮(zhèn)下煤層開(kāi)采條件，因此確定開(kāi)采方案為傾向變條帶開(kāi)采。這種開(kāi)采方案中，由東向西逐漸加寬開(kāi)采條帶，根據(jù)開(kāi)采地表建筑物受損情況可以及時(shí)調(diào)整開(kāi)采寬度，降低地表建筑物的損害程度，確保了建筑物的安全使用[7]。

2.4FLAC

3D數(shù)值模擬以方案Ⅰ中的參數(shù)進(jìn)行數(shù)值模擬計(jì)算，進(jìn)一步確定開(kāi)采方案。模擬煤層的埋深為238.7m，模擬模型全長(zhǎng)600m，高270m，寬550m，共劃分147400個(gè)單元。數(shù)值模擬9個(gè)開(kāi)采工作面，模擬沿井田邊界至上山方向450m的傾斜變條帶開(kāi)采。開(kāi)采結(jié)束后豎向位移的走向剖面圖見(jiàn)圖3，最大主應(yīng)力的走向剖面圖見(jiàn)圖4。計(jì)算結(jié)論表明，采用留煤柱開(kāi)采法，開(kāi)采后地表最大豎向位移量小于75mm，地表沉陷引起的水平變形遠(yuǎn)遠(yuǎn)小于1mm/m，不會(huì)引起地表建筑物的變形損害，因此在開(kāi)采過(guò)程中地表基本上保持穩(wěn)定。FLAC3D數(shù)值模擬與YHL-12影響函數(shù)的預(yù)計(jì)地表移動(dòng)變形參數(shù)和結(jié)果吻合較好，從而確定開(kāi)采方案為傾向變條帶開(kāi)采。

3結(jié)語(yǔ)

煤層開(kāi)采引起的地表移動(dòng)變形是西部礦區(qū)開(kāi)采的一個(gè)特殊問(wèn)題,采動(dòng)引起的開(kāi)采沉陷問(wèn)題具有一定的特殊性。借助于計(jì)算機(jī)技術(shù)，運(yùn)用YhL-12預(yù)計(jì)軟件對(duì)蒲白礦務(wù)局西固鎮(zhèn)下開(kāi)采的地表移動(dòng)變形進(jìn)行了預(yù)計(jì),并與數(shù)值模擬計(jì)算結(jié)果進(jìn)行了對(duì)比，分析了條帶開(kāi)采時(shí)地表沉陷的基本規(guī)律及特點(diǎn)，確定了開(kāi)采方案，對(duì)礦區(qū)綠色開(kāi)采、生態(tài)環(huán)境保護(hù)和可持續(xù)發(fā)展具有重要的指導(dǎo)意義。