上向分段充填開采覆巖移動(dòng)變形及影響范圍預(yù)測(cè)

上向分段充填開采覆巖移動(dòng)變形及影響范圍預(yù)測(cè)

——以貴州某磷礦為例摘要：為進(jìn)一步揭示上向分段充填開采覆巖移動(dòng)變形規(guī)律,以貴州A磷礦為例,采用理論計(jì)算與數(shù)值仿真分析相結(jié)合的方法,對(duì)上向分段充填開采后地表變形及巖體移動(dòng)對(duì)相鄰礦山(B、C、D磷礦)的影響范圍進(jìn)行了分析預(yù)測(cè),結(jié)果表明:a、b礦層整體開采后的最大地表傾斜、水平變形和曲率值均未達(dá)到標(biāo)準(zhǔn)規(guī)定的地表變形臨界值;各分段依次開采后的地表總體最大位移量約為84mm、地表最大傾斜變形值為0.12mm/m、最大水平變形值為0.0193mm/m,現(xiàn)階段A磷礦開采后地表整體上不會(huì)產(chǎn)生較大的移動(dòng)變形;隨著+900m、+920m、+937m、+951m分段的逐步上向開采,巖體移動(dòng)對(duì)相鄰礦山的影響范圍逐漸擴(kuò)大,對(duì)C磷礦影響范圍最大、B磷礦次之、D磷礦相對(duì)較小。建議礦山在今后開采過程中應(yīng)重點(diǎn)加強(qiáng)因A磷礦開采巖體移動(dòng)對(duì)C磷礦影響的預(yù)測(cè)預(yù)警及礦山壓力管理工作,同時(shí)加強(qiáng)A磷礦對(duì)礦區(qū)范圍礦層露頭一帶高陡邊坡的巡查及地表變形監(jiān)測(cè)工作,以確保礦山安全開采。關(guān)鍵詞：上向分段充填開采;巖體移動(dòng);概率積分法;地表變形;影響范圍;數(shù)值仿真0引言磷礦作為非碎屑沉積巖,是一種重要的非金屬礦產(chǎn)資源,主要用于化工和農(nóng)業(yè)生產(chǎn)。貴州磷礦資源儲(chǔ)量豐富,居全國(guó)第二位,且賦存較為集中,主要分布在開陽、息烽、甕安及福泉等地,多礦井相鄰集中開采的情況比較常見。已有研究表明,井工礦山開采過程中形成的巖體移動(dòng)范圍會(huì)對(duì)鄰近礦山開采造成影響,如削弱礦區(qū)及鄰近范圍內(nèi)地表邊坡的穩(wěn)定性、誘發(fā)地表移動(dòng)變形與破壞等,對(duì)自身及周邊礦井的安全產(chǎn)生重大影響。為此,國(guó)家礦山安全監(jiān)察局在2022年印發(fā)的《金屬非金屬礦山重大事故隱患判定標(biāo)準(zhǔn)》中指出,相鄰礦山開采巖體移動(dòng)范圍存在交叉重疊等相互影響時(shí),未按設(shè)計(jì)留設(shè)保安礦(巖)柱或者采取其他措施的視為重大事故隱患。因此,金屬、非金屬地下礦山建設(shè)和生產(chǎn)過程中必須分析并圈定礦山開采巖體移動(dòng)范圍,從而為留設(shè)保安礦(巖)柱或者采取其他措施提供依據(jù)。近年來,一些學(xué)者對(duì)井工礦山開采地表沉陷規(guī)律展開了研究,結(jié)合下沉、水平移動(dòng)、傾斜、曲率等參數(shù)對(duì)地表沉陷規(guī)律進(jìn)行了分析,并提出了針對(duì)性的措施和建議。如程立年等對(duì)“三下”開采移動(dòng)帶內(nèi)地表穩(wěn)定性進(jìn)行了分析,采用概率積分法計(jì)算出了礦山安全開采深度;]研究發(fā)現(xiàn)地表開采沉陷在受自重應(yīng)力作用影響的同時(shí),還受構(gòu)造應(yīng)力作用和礦層傾角大小、頂板巖體節(jié)理發(fā)育程度與結(jié)構(gòu)形式、開采方法與放礦規(guī)律等的影響;通過對(duì)不同充填體灰砂質(zhì)量比和充填率影響下的地表移動(dòng)進(jìn)行了模擬,得出膏體充填開采下覆巖移動(dòng)角會(huì)隨著覆巖水平高度的增加而逐漸增大;分析了塌陷坑移動(dòng)角和遠(yuǎn)近地表圍巖隨開采深度的變化規(guī)律。以上學(xué)者從安全開采深度、構(gòu)造應(yīng)力作用、不同充填率等方面對(duì)地表移動(dòng)變形破壞進(jìn)行了研究,為礦山地表沉陷、移動(dòng)變形等問題提供了有益參考。然而,現(xiàn)有研究大多集中在礦井開采引起的地表移動(dòng)變形、沉陷分析及監(jiān)測(cè)等方面,針對(duì)礦山開采巖體移動(dòng)范圍對(duì)相鄰礦山的影響鮮有提及。因此,本文以貴州A磷礦為例,采用理論計(jì)算與數(shù)值仿真分析相結(jié)合的方法,對(duì)上向分段充填開采后的地表移動(dòng)變形及巖體移動(dòng)對(duì)相鄰礦山的影響范圍進(jìn)行了分析和預(yù)測(cè),以期為類似條件下多個(gè)磷礦相鄰開采的優(yōu)化提供參考。1礦井概況1.1地質(zhì)概況貴州A磷礦礦區(qū)面積0.4273km2,開采深度為+1300～+900m標(biāo)高,礦區(qū)內(nèi)總體地形為一東高西低、向北西傾斜的斜坡,礦區(qū)環(huán)境較復(fù)雜,同時(shí)與B、C、D磷礦相鄰(見圖1)。礦區(qū)含磷巖組為陡山沱組,由磷塊巖、白云巖、硅質(zhì)巖等組成,礦體平均傾角83°,為急傾斜礦體。礦層由兩層磷塊巖礦層組成(見圖2),a層礦平均厚度18.17m,b層礦平均厚度15.54m,總體趨勢(shì)是東厚西薄。圖1礦區(qū)相鄰位置關(guān)系Fig.1Locationdiagramofadjacentminingareas圖2礦層特征圖Fig.2Characteristicsoforebed1.2開采現(xiàn)狀及接續(xù)計(jì)劃設(shè)計(jì)中將A磷礦劃分為一采區(qū)和二采區(qū)兩個(gè)采區(qū),其中一采區(qū)已開采完畢,形成了采空區(qū);二采區(qū)為當(dāng)前開采區(qū)域,礦體開采標(biāo)高為+980～+900m,按照礦體賦存條件及礦山實(shí)際情況,將二采區(qū)劃分成4個(gè)分段開采,即+900m、+920m、+937m、+951m分段,分段垂高14～20m。采用上向分段充填采礦法開采,礦房長(zhǎng)度為10～25m。首采分段為+900m分段,礦房?jī)?nèi)開采順序是先采b層礦,再采a層礦。2地表移動(dòng)變形理論計(jì)算2.1概率積分法概率積分法又稱隨機(jī)介質(zhì)理論法,是一種以隨機(jī)介質(zhì)理論為基礎(chǔ)的開采沉陷預(yù)計(jì)方法,是目前較為成熟、應(yīng)用最廣泛的地表移動(dòng)變形預(yù)計(jì)方法。該方法是由我國(guó)學(xué)者在隨機(jī)介質(zhì)理論的基礎(chǔ)上逐步完善發(fā)展而來的。依據(jù)概率積分法原理,地表最大傾斜變形imax、最大水平變形εmax、最大曲率Kmax與地表最大下沉值Wmax、開采深度H及主要影響角正切值tanβ之間存在如下關(guān)系:Wmax=mqcosα

,(1)(2)(3)(4)(5)式中:m表示礦層厚度,m;q表示下沉系數(shù),根據(jù)礦區(qū)實(shí)際測(cè)量資料確定;α表示礦層傾角,(°);r表示采空區(qū)邊界地表主要影響范圍半徑,m;tanβ為主要影響角的正切值;b表示水平移動(dòng)系數(shù),按礦區(qū)實(shí)際測(cè)量資料確定。2.2計(jì)算結(jié)果分析根據(jù)該礦地質(zhì)資料及實(shí)際生產(chǎn)條件,本文分析的開采區(qū)域礦體主要由17#勘探線控制,取下沉系數(shù)q為0.04、主要影響角正切值tanβ為1.5、水平移動(dòng)系數(shù)b為0.25,與其余相關(guān)參數(shù)值一并代入式(1)-式(5),計(jì)算得出地表下沉量、傾斜值、曲率值、水平變形值等(見表1)。表1地表移動(dòng)變形計(jì)算結(jié)果Table1Calculationresultsofgroundmovementanddeformation由表1可知,A磷礦a層礦、b層礦開采后地表移動(dòng)變形產(chǎn)生的最大區(qū)域地表傾斜值、水平變形值和曲率值均未達(dá)到GB50771-2012《有色金屬采礦設(shè)計(jì)規(guī)范》規(guī)定的地表變形臨界值,因此,預(yù)測(cè)A磷礦在a層礦、b層礦開采完畢后,整體上不會(huì)在地表產(chǎn)生大的移動(dòng)變形。3地表移動(dòng)變形數(shù)值仿真分析3.1數(shù)值仿真模型建立采用FLAC3D建立數(shù)值仿真模型,開展A磷礦上向分段充填開采對(duì)地表影響范圍的模擬分析。為使模型最大程度地接近真實(shí)情況,依據(jù)地質(zhì)資料,并結(jié)合A磷礦與各勘探線位置關(guān)系建立數(shù)值模型(見圖3)。模型x方向與勘探線平行,長(zhǎng)度1500m;y方向與勘探線方向垂直,長(zhǎng)度1000m;z方向?yàn)樨Q直方向,最低標(biāo)高為600m。模型共劃分為41175個(gè)單元和50128個(gè)節(jié)點(diǎn)。圖3三維數(shù)值仿真模型Fig.3Three-dimensionalnumericalsimulationmodel通過開展井下工程地質(zhì)調(diào)查及巖石力學(xué)室內(nèi)試驗(yàn),并采用Hoek-Brown經(jīng)驗(yàn)公式對(duì)各項(xiàng)參數(shù)進(jìn)行折減,同時(shí)參考區(qū)域內(nèi)其他礦山數(shù)據(jù),綜合確定模型巖石物理力學(xué)參數(shù)(見表2)。表2巖石物理力學(xué)參數(shù)Table2Thephysicalandmechanicalparametersofrock3.2分析計(jì)算分析過程中采用建(構(gòu))筑物損壞等級(jí)評(píng)判標(biāo)準(zhǔn)來驗(yàn)證移動(dòng)范圍內(nèi)的沉降變形。目前,礦山回采作業(yè)面主要集中于+900～+920m水平,采用充填法處理采空區(qū)。結(jié)合圖3對(duì)各分段進(jìn)行模擬開采,同時(shí)為便于分析A磷礦各分段開采后對(duì)地表的影響,分別在模型地表位置布置一號(hào)測(cè)線(x方向)和二號(hào)測(cè)線(y方向)(見圖4)。圖4模型地表測(cè)線布置圖Fig.4Groundsurveylayoutofmodel根據(jù)各監(jiān)測(cè)點(diǎn)的水平位移和豎直位移數(shù)值仿真數(shù)據(jù),可計(jì)算得出相鄰監(jiān)測(cè)點(diǎn)的水平變形值和傾斜變形值,計(jì)算公式如下:(6)(7)(8)式中:in為地表n點(diǎn)的傾斜值,mm/m;Tn+1、Tn分別為地表n點(diǎn)和n+1點(diǎn)的下沉值,mm;Kn為地表n點(diǎn)的曲率值,mm/m2;Dn～n+1為地表n點(diǎn)和n+1點(diǎn)的水平距離,m;εn為地表n點(diǎn)的水平變形值,mm/m;Hn+1、Hn分別為地表n點(diǎn)和n+1點(diǎn)的水平移動(dòng)值,mm。3.3模擬結(jié)果分析將一號(hào)、二號(hào)測(cè)線測(cè)得的位移代入式(6)-式(8)計(jì)算得到+900、+920、+937、+951m分段開采后的地表位移變化(見圖5)。由圖5可以看出:沿y方向的傾斜值比沿x方向的傾斜值大,最大值0.12mm/m為+900m分段開采后的地表傾斜值;各分段沿y方向傾斜值按開采順序依次減小,沿x方向傾斜最大值出現(xiàn)在+937m分段,為2.74×10-2

mm/m;沿x、y方向水平變形最大值均出現(xiàn)在+951m分段,最大值為1.93×10-2

mm/m,且變形值呈逐漸增大的趨勢(shì)。圖5各分段開采后地表位移變化圖Fig.5Thegrounddisplacementchangeduponeachsublevelmining因文章篇幅所限,以下僅給出+951m分段開采后的地表位移變形云圖(見圖6、圖7)。圖6+951m分段開采后地表垂直位移變形云圖Fig.6Cloudchartfortheverticalgrounddisplacementafterminingatlevel+951m圖7+951m分段開采后地表位移與礦界相對(duì)位置示意圖Fig.7Relativepositiondiagramforgrounddisplacementandminingboundaryatlevel+951m從圖中7可以看出,地表整體位移變形主要集中在礦井范圍內(nèi),緊鄰C磷礦,地表沿x方向的水平位移最大值約為9.9mm,沿y方向的水平位移最大值約為0.5mm,沿z方向的垂直位移最大值為83.5mm,地表總體最大位移量約為84mm。計(jì)算出+951m分段開采后該區(qū)域下沉盆地的地表傾斜值、曲率、水平變形值,從x、y方向算得最大地表傾斜值、曲率、水平變形值分別為:x方向傾斜值1.36×10-2

mm/m、曲率5.77×10-5

mm/m2、水平變形值1.1×10-2

mm/m;y方向傾斜值6.02×10-2

mm/m、曲率1.2×10-4

mm/m2、水平變形值1.93×10-2

mm/m。根據(jù)式(6)-式(8),并結(jié)合數(shù)值模擬結(jié)果計(jì)算得出不同開采時(shí)期的地表變形值(見表3)。由表3可知,A磷礦各分段開采后的地表最大傾斜值為0.12mm/m、最大水平變形值為0.0193mm/m,均未達(dá)到GB50771-2012規(guī)定的地表變形臨界值。因此,A磷礦現(xiàn)階段開采引起的地表變形整體上較小,初步預(yù)測(cè)對(duì)相鄰的B磷礦、C磷礦及D磷礦開采影響程度較弱。表3不同開采時(shí)期的地表變形值Table3Grounddeformationindifferentminingphases4巖體移動(dòng)對(duì)相鄰礦山影響范圍的預(yù)測(cè)分析由于巖體長(zhǎng)期強(qiáng)度和采空區(qū)穩(wěn)定性的時(shí)間效應(yīng),現(xiàn)階段地表位移較大的區(qū)域有可能成為礦山長(zhǎng)期開采潛在的地表移動(dòng)范圍,對(duì)比數(shù)值仿真地表移動(dòng)變形區(qū)域,并結(jié)合位移變化可預(yù)測(cè)得出A磷礦各分段在上向充填開采后巖體移動(dòng)對(duì)相鄰礦山的影響范圍(見圖8)。當(dāng)開采+900m分段礦體時(shí),礦體埋深大于400m,地下開采對(duì)地表影響程度較小,潛在的開采移動(dòng)影響范圍主要在16#-17#勘探線之間,影響范圍呈類橢圓形[見圖8(a)],長(zhǎng)軸方向長(zhǎng)約547m,短軸方向長(zhǎng)約434m,擴(kuò)展到C磷礦礦界內(nèi)約175m。由該礦地質(zhì)資料可知,該區(qū)域?yàn)榈V井主要開采部分,受采掘影響,上覆老采空區(qū)發(fā)生下沉,且采空區(qū)上方為高陡山體邊坡,為潛在的開采移動(dòng)影響范圍。當(dāng)開采+920m、+937m分段礦體時(shí),下分段采空區(qū)已完成充填,此時(shí)C磷礦受到的影響范圍較大,該區(qū)域地形為高邊坡并向A磷礦礦井延伸,因此隨著開采范圍的擴(kuò)大,形成的潛在開采移動(dòng)影響范圍逐漸向C磷礦擴(kuò)展,由175m增加到189m,而長(zhǎng)軸方向由530m增加到601m,整體影響范圍分別呈類橢圓形[見圖8(b)]、類三角形[見圖8(c)],對(duì)B磷礦的影響逐漸變小。當(dāng)開采+951m分段礦體時(shí),上覆地層所受擾動(dòng)范圍和程度亦增大,位移變化最大值處向南移動(dòng),同時(shí)影響范圍再次擴(kuò)展到相鄰的C磷礦,整體影響范圍呈類三角形[見圖8(d)],長(zhǎng)軸方向長(zhǎng)約683m,短軸方向長(zhǎng)約400m,下部影響范圍擴(kuò)展到距礦界約189m處。礦井所處區(qū)域地勢(shì)為北高南低、東高西低,主要開采范圍處于該磷礦左部邊界處,此處地表為高陡邊坡并向北延伸至相鄰礦井,因此潛在的開采移動(dòng)影響范圍向B磷礦擴(kuò)展。開采+951m分段礦體時(shí),因礦脈傾斜向下,老采空區(qū)位于礦井右翼邊界處,此處地表由D磷礦一路延伸而來形成高陡邊坡,因此潛在的開采移動(dòng)影響范圍再次向D磷礦擴(kuò)展。圖8各分段依次開采后巖體移動(dòng)對(duì)相鄰礦山影響范圍預(yù)測(cè)Fig.8Thepredictionofinfluencingrangeofrockmassmovementafterminingateachsublevelonadjacentblocks

5開采建議綜合理論計(jì)算和數(shù)值仿真分析可知,相比于其他礦井而言,A磷礦開采后對(duì)C磷礦的影響范圍最大,B磷礦次之,D磷礦相對(duì)較小。因此,提出如下開采建議:a.后期開采過程中要重點(diǎn)加強(qiáng)因A磷礦開采巖體移動(dòng)對(duì)C磷礦影響的預(yù)測(cè)預(yù)警工作。b.相鄰礦井應(yīng)重視本礦山處于A磷礦開采巖體移動(dòng)影響范圍內(nèi)的礦山壓力管理工作,嚴(yán)格按照設(shè)計(jì)留設(shè)相應(yīng)保護(hù)礦(巖)柱,加強(qiáng)影響區(qū)域采場(chǎng)和巷道支護(hù)。c.由于A磷礦+980m以上存在采用空?qǐng)?/p>