嚴(yán)鵬　陳紹民　黃敏　黃志國　張寶

摘要：合理的結(jié)構(gòu)參數(shù)對維護(hù)采場穩(wěn)定、提高單個采場產(chǎn)能、降低采切工程量具有重要意義。以蒙庫鐵礦為研究對象，利用巖體遙測與結(jié)構(gòu)分析系統(tǒng)Sirovision對采場頂板、側(cè)幫及上盤圍巖節(jié)理面進(jìn)行調(diào)查和數(shù)字化識別，獲取采場節(jié)理組的空間方位統(tǒng)計信息及優(yōu)勢節(jié)理產(chǎn)狀；基于Mathews穩(wěn)定圖解法對沿走向布置與垂直走向布置采場的頂板、側(cè)幫及上盤開展穩(wěn)定性分析，得到采場頂板、側(cè)幫及上盤的實(shí)際水力半徑和允許水力半徑的相關(guān)關(guān)系；采用數(shù)值分析方法，通過分析礦柱的位移、最大剪切應(yīng)變率及采場的主應(yīng)力等因素，確定了2種典型采場結(jié)構(gòu)礦柱的最小安全厚度。研究結(jié)果表明：采場沿礦體走向布置時，當(dāng)?shù)V體厚度小于20 m時，長度控制在40 m之內(nèi)的采場穩(wěn)定性較好；采場垂直礦體走向布置時，礦體厚度不超過30 m，寬度小于25 m的采場穩(wěn)定性較好；對于沿走向布置的40 m×15 m典型采場，礦柱最小安全厚度為6 m；對于垂直走向布置的30 m×25 m典型采場，礦柱最小安全厚度為8 m。研究結(jié)果可為蒙庫鐵礦采場結(jié)構(gòu)參數(shù)優(yōu)化提供依據(jù)，并能為國內(nèi)外同類礦山采場結(jié)構(gòu)設(shè)計提供參考。

關(guān)鍵詞：厚大礦體；Mathews穩(wěn)定圖解法；采場結(jié)構(gòu)參數(shù)；數(shù)值分析；礦柱厚度

中圖分類號：TD35文章編號：1001-1277（2023）09-0046-08

文獻(xiàn)標(biāo)志碼：Adoi：10.11792/hj20230908

引 言

中國經(jīng)濟(jì)從高速增長階段轉(zhuǎn)變?yōu)楦哔|(zhì)量發(fā)展階段，經(jīng)濟(jì)結(jié)構(gòu)的改變導(dǎo)致其對礦產(chǎn)資源的依賴程度增加，需要消耗大量的礦產(chǎn)資源來支撐［1］。然而，中國礦產(chǎn)資源整體開發(fā)利用率較低，存在一定程度的資源浪費(fèi)現(xiàn)象，制約礦產(chǎn)資源開發(fā)向規(guī)?；?、集約化方向發(fā)展［2-3］。蒙庫鐵礦東段礦體總體走向SE—NW，傾向SW，傾角63°～85°，整體呈不連續(xù)分布，8#勘探線—18#勘探線礦體厚10～30 m，是典型的急傾斜厚大礦體。因此，結(jié)合礦山工程地質(zhì)條件、工藝技術(shù)水平及現(xiàn)場工程布置等情況，從保障礦山生產(chǎn)能力、降低采礦成本、提高采礦效率及確保采場的穩(wěn)定性和安全性等方面出發(fā)，對厚大礦體合理采場結(jié)構(gòu)參數(shù)的確定及優(yōu)化開展研究工作具有重要意義。

目前，國內(nèi)外相關(guān)領(lǐng)域?qū)W者利用工程類比法、理論分析法和數(shù)值模擬法等方法，對采場結(jié)構(gòu)參數(shù)的確定及優(yōu)化進(jìn)行了豐富的研究［4-7］。其中，Mathews穩(wěn)定圖解法和數(shù)值模擬法因其實(shí)用性及準(zhǔn)確性得到了廣泛應(yīng)用［8- 9］。趙興東等［10］利用3種常見的巖體質(zhì)量分級方法進(jìn)行巖體質(zhì)量分級綜合評價，采用Mathews穩(wěn)定圖解法對采場結(jié)構(gòu)參數(shù)進(jìn)行初步計算，在此基礎(chǔ)上借助數(shù)值模擬軟件對采場穩(wěn)定性進(jìn)行分析，最終確定了最優(yōu)的采場結(jié)構(gòu)參數(shù)。劉建東等［11］針對礦山開采過程中出現(xiàn)的地壓災(zāi)害現(xiàn)象，通過理論計算、Mathews穩(wěn)定圖解法及數(shù)值模擬分析等方法對采場穩(wěn)定性進(jìn)行綜合評價，得到合理的采場結(jié)構(gòu)參數(shù)，確保高應(yīng)力條件下礦體高效開采的穩(wěn)定性。郭進(jìn)平等［12-13］為了解決礦體開采過程中采場頂板的穩(wěn)定性問題，應(yīng)用Mathews穩(wěn)定圖解法對采場結(jié)構(gòu)參數(shù)進(jìn)行初步選擇，借助Flac3D軟件對采場結(jié)構(gòu)參數(shù)進(jìn)行針對性的分析，為礦山提供合理的采場結(jié)構(gòu)參數(shù)。

本研究通過現(xiàn)場工程地質(zhì)調(diào)查獲取采場不同區(qū)域內(nèi)巖體分布特征，利用Q系統(tǒng)分級評價方法對礦（巖）體進(jìn)行綜合評價，在此基礎(chǔ)上，結(jié)合Mathews穩(wěn)定圖解法對不同賦存條件下的礦體采場結(jié)構(gòu)參數(shù)進(jìn)行初步確定，利用Flac3D數(shù)值模擬軟件對采場穩(wěn)定性進(jìn)行分析，最終確定合理的采場結(jié)構(gòu)參數(shù)，并為類似礦山采場結(jié)構(gòu)參數(shù)設(shè)計提供參考依據(jù)。

1 采場礦巖工程地質(zhì)調(diào)查

采用三維數(shù)字?jǐn)z影測量和巖體結(jié)構(gòu)分析系統(tǒng)（Sirovision），對深度為500 m左右的礦體頂板、側(cè)幫及下盤圍巖典型區(qū)域進(jìn)行節(jié)理裂隙掃描測量，統(tǒng)計節(jié)理裂隙分布規(guī)律并進(jìn)行優(yōu)勢節(jié)理劃分，如圖1所示。根據(jù)調(diào)查結(jié)果，采場頂板礦巖巖組中控制性的節(jié)理產(chǎn)狀主要發(fā)育有2組：143°∠83°，39°∠76°；采場頂板為近水平狀（傾角取0°），產(chǎn)狀為221°∠0°。采場側(cè)幫的礦巖巖組中控制性的節(jié)理產(chǎn)狀主要發(fā)育有2組：143°∠83°，39°∠76°；采場側(cè)幫為豎直狀，產(chǎn)狀為221°∠90°。采場上盤圍巖均為變粒巖，變粒巖巖組中主要發(fā)育3組控制性節(jié)理：324°∠85°，111°∠74°，318°∠58°；采場上盤產(chǎn)狀與礦體產(chǎn)狀類似，取221°∠72°（取礦體平均傾角72°）。

根據(jù)Q系統(tǒng)分級評價方法［13］，對礦巖和變粒巖進(jìn)行分級評價，Q值計算見式（1）：

式中：RQD為巖石質(zhì)量指標(biāo)；Jn為巖體的節(jié)理數(shù)；Jr為巖體的節(jié)理粗糙度系數(shù)；Ja為巖體的節(jié)理蝕變系數(shù)；Jw為節(jié)理滲水折減系數(shù)；SRF為應(yīng)力折減系數(shù)。

根據(jù)Q系統(tǒng)參數(shù)評分表，并結(jié)合礦巖節(jié)理裂隙的實(shí)際調(diào)查結(jié)果，礦巖和變粒巖的評價結(jié)果如表1所示。

2 計算采場穩(wěn)定性系數(shù)

本文基于Mathews穩(wěn)定圖解法設(shè)計采場結(jié)構(gòu)參數(shù)，Mathews穩(wěn)定圖解法是預(yù)測采場穩(wěn)定性的實(shí)用方法，在工程中得到了廣泛應(yīng)用。Mathews穩(wěn)定圖解法設(shè)計過程以穩(wěn)定系數(shù)和水力半徑的計算為基礎(chǔ)，根據(jù)這2個因子將采場劃分為穩(wěn)定區(qū)、無支護(hù)過渡區(qū)、支護(hù)穩(wěn)定區(qū)、支護(hù)過渡區(qū)和崩落區(qū)。

穩(wěn)定性系數(shù)反映了在一定應(yīng)力條件下巖體的自穩(wěn)能力，其計算見式（2）：

式中：N為穩(wěn)定性系數(shù)；Q′為修正的Q值；A為巖石應(yīng)力系數(shù)；B為節(jié)理方位系數(shù)；C為重力調(diào)整系數(shù)。

水力半徑考慮了單獨(dú)采場暴露表面的尺寸和形狀，水力半徑可用表面積與暴露面周長的比值來表示，計算見式（3）：

式中：HR為水力半徑（m）；l為采場的長度（m）；h為采場的寬度或高度（m）。

Mathews穩(wěn)定性系數(shù)與水力半徑之間的關(guān)系如圖2所示，通過該圖來確定容許的水力半徑。

1）Q′值的計算。Mathews穩(wěn)定圖解法使用修正后的NGI隧道質(zhì)量指標(biāo)值Q′，Q′值的計算見式（4）：

根據(jù)表1，可計算出礦巖和變粒巖的Q′值分別為30.67和19.78。

2）系數(shù)A、B、C的計算。巖石應(yīng)力系數(shù)A與完整巖體的單軸抗壓強(qiáng)度和平行開挖面的最大誘導(dǎo)應(yīng)力的比值呈線性關(guān)系。A值為0.1～1.0，其值可根據(jù)圖3-a）確定。結(jié)合室內(nèi)巖石力學(xué)試驗(yàn)及現(xiàn)場應(yīng)力測試得到礦巖及變粒巖的值分別為8.12，7.47，結(jié)合圖3-a），可以確定礦巖及變粒巖的巖石應(yīng)力系數(shù)A的取值分別為0.79，0.72。

節(jié)理方位系數(shù)B受不連續(xù)面的影響，其值可根據(jù)控制性節(jié)理與采場表面的相對方位確定（如圖3-b）所示）。根據(jù)前述的地質(zhì)調(diào)查結(jié)果，可以得到采場節(jié)理方位系數(shù)B的計算結(jié)果如表2所示。

結(jié)合最不利原則，進(jìn)行采場穩(wěn)定性分析時，頂板取B=0.91（礦巖）；側(cè)幫取B=0.20（礦巖）；上盤取B=0.70（變粒巖）。

重力調(diào)整系數(shù)C考慮了重力對采場暴露表面崩落、滑落等穩(wěn)定性的影響（如圖3-c）所示），重力調(diào)整系數(shù)和采場表面傾角的關(guān)系式為：

C=8-6cos α （5）

式中：α為采場工作面傾角（°）。

根據(jù)采場頂板、側(cè)幫及上盤的產(chǎn)狀，通過式（5）可計算出重力調(diào)整系數(shù)C的值分別為2.0，8.0，6.2。

根據(jù)式（2）及上述系數(shù)的取值，得到相關(guān)的參數(shù)及計算結(jié)果見表3。由表3可知，采場頂板、側(cè)幫及上盤的穩(wěn)定性系數(shù)N分別為45.91，38.66，61.81。

3 采場結(jié)構(gòu)參數(shù)設(shè)計

根據(jù)蒙庫鐵礦東段500 m以上采場實(shí)際布置情況，并結(jié)合礦體的賦存條件，采用Mathews穩(wěn)定圖解法分析采場穩(wěn)定性?，F(xiàn)場工程地質(zhì)調(diào)查結(jié)果表明，蒙庫鐵礦東段500 m以上礦體直接圍巖以變粒巖為主，夾少量石榴子石，礦體主要為磁鐵礦或塊狀石榴子石磁鐵礦。

礦體500 m中段采用階段空場采礦法開采，階段高度為50 m，當(dāng)?shù)V體厚度小于20 m時，采場沿礦體走向布置，當(dāng)?shù)V體厚度大于20 m時，采場垂直礦體走向布置。

3.1 沿礦體走向布置

當(dāng)采場沿礦體走向布置時，根據(jù)水力半徑的計算方法，得出不同采場結(jié)構(gòu)參數(shù)下采場頂板的水力半徑，以及采場側(cè)幫及上盤水力半徑，計算結(jié)果分別如圖4、表4及表5所示。

由圖4可知，當(dāng)?shù)V體厚度小于20 m，采場長度小于40 m時，采場頂板最大水力半徑僅為6.67 m，小于頂板穩(wěn)定區(qū)的允許水力半徑11.2 m，表明采場頂板的穩(wěn)定性較好；同理，根據(jù)表4、表5可以看出，采場側(cè)幫的最大水力半徑為6.67 m及上盤的最大水力半徑為10 m，也都小于它們穩(wěn)定區(qū)的允許水力半徑10.51 m和12.52 m。因此，當(dāng)?shù)V體厚度小于20 m時，采場沿礦體走向布置，當(dāng)采場長度控制在40 m之內(nèi)，穩(wěn)定性較好，能確保開采時采場安全。

3.2 垂直礦體走向布置

當(dāng)采場垂直礦體走向布置時，根據(jù)水力半徑的計算方法，得出不同采場結(jié)構(gòu)參數(shù)下采場頂板水力半徑，以及采場側(cè)幫及上盤水力半徑，計算結(jié)果如圖5、表6及表7所示。

根據(jù)Mathews穩(wěn)定圖解法采場頂板容許水力半徑的計算，穩(wěn)定區(qū)的水力半徑小于11.2 m，由圖5可知，當(dāng)?shù)V體厚度小于30 m時，寬度為25 m的采場對應(yīng)的水力半徑為6.82 m，表明采場頂板穩(wěn)定性較好；高度為40 m，礦體厚度小于30 m的采場側(cè)幫最大水力半徑為8.57 m，小于側(cè)幫穩(wěn)定區(qū)允許的水力半徑10.51 m；而寬度小于25 m的采場上盤最大水力半徑為7.69 m，小于上盤穩(wěn)定區(qū)允許的水力半徑12.52 m。因此，采場垂直礦體走向布置，長度小于25 m采場，其頂板、側(cè)幫及上盤的穩(wěn)定性均較好，能確保開采時采場安全。

4 礦柱寬度設(shè)計

在礦巖質(zhì)量分級基礎(chǔ)上采用Mathews穩(wěn)定圖解法，通過考慮采場頂板、側(cè)幫及上盤的穩(wěn)定性，得到了采場安全回采的結(jié)構(gòu)參數(shù)。Mathews穩(wěn)定圖解法是一種經(jīng)驗(yàn)的參數(shù)選擇方法，盡管沒有考慮應(yīng)力環(huán)境對采場結(jié)構(gòu)參數(shù)的影響，但為采場結(jié)構(gòu)參數(shù)選取提供了可供參考的范圍。本節(jié)基于Flac3D數(shù)值分析軟件，利用Mathews穩(wěn)定圖解法確定采場結(jié)構(gòu)參數(shù)，確定合理的礦柱寬度。

4.1 數(shù)值模擬方案

根據(jù)采場布置形式選擇典型的采場參數(shù)，分析采場間礦柱的寬度。對于沿礦體走向布置采場，選擇的典型采場參數(shù)為40 m×15 m，如圖6-a）所示；垂直走向布置的采場，選擇的典型采場參數(shù)為30 m×15 m，如圖6-b）所示。

運(yùn)用巖體工程地質(zhì)力學(xué)方法進(jìn)行工程地質(zhì)巖組劃分、巖體結(jié)構(gòu)分類，結(jié)合巖體分級結(jié)果及原有資料綜合選取巖體力學(xué)參數(shù)，得到蒙庫鐵礦東段500 m以上礦體及圍巖的物理力學(xué)參數(shù)，如表8所示。

礦柱是決定采場穩(wěn)定狀態(tài)的重要結(jié)構(gòu)單元，礦柱結(jié)構(gòu)破壞必然引發(fā)采場力學(xué)狀態(tài)發(fā)生變化，影響整個采空區(qū)穩(wěn)定性。本次數(shù)值模擬方案在固定跨度條件下，通過研究不同礦柱厚度對應(yīng)的采場及礦柱的穩(wěn)定性，進(jìn)而確定礦柱最小安全厚度，分析方案如表9所示。

4.2 沿走向布置采場礦柱最小安全厚度的確定

位移是反應(yīng)采場結(jié)構(gòu)穩(wěn)定性的重要參數(shù)，采場某些區(qū)域的位移過大預(yù)示著這些區(qū)域很可能發(fā)生局部失穩(wěn)。不同厚度礦柱的位移云圖如圖7所示。當(dāng)?shù)V柱厚度為8 m時，礦柱的最大位移為1.85 cm;礦柱厚度為6 m、4 m對應(yīng)的最大位移分別為2.46 cm、3.97 cm。這表明，隨著礦柱厚度的減小，礦柱的最大位移逐漸增加，并且礦柱厚度越小，礦柱最大位移的增量越大。

最大剪切應(yīng)變率則可反映采場中發(fā)生剪切破壞的區(qū)域，剪切應(yīng)變率大的區(qū)域則形成剪切破壞帶。不同厚度礦柱的最大剪切應(yīng)變率分布云圖如圖8所示。礦柱厚度為8 m時，礦柱表面絕大部分的剪切應(yīng)變率處于隨機(jī)分布的狀態(tài)，礦柱整體剪切變形比較均勻，只有礦柱中間部位出現(xiàn)了剪切應(yīng)變率較大的區(qū)域；礦柱厚度變?yōu)? m時，中間部位的剪切變形區(qū)域向礦柱頂?shù)撞繑U(kuò)展，但是沒有形成連續(xù)的剪切變形條帶；隨著礦柱厚度的進(jìn)一步降低，礦柱表面已經(jīng)形成了連續(xù)的剪切變形條帶，整個礦柱的剪切變形呈現(xiàn)明顯的不均勻狀態(tài)，礦柱容易沿著形成的剪切變形條帶發(fā)生剪切破壞。

最大主應(yīng)力和最小主應(yīng)力的分布對于判斷采場穩(wěn)定性是至關(guān)重要的。采場的最大主應(yīng)力分布云圖如圖9所示。由圖9可知：不同礦柱采場的最大主應(yīng)力大小幾乎相同，分布十分相似，盡管礦柱側(cè)幫及采場側(cè)壁處都出現(xiàn)了拉應(yīng)力，拉應(yīng)力的最大值約為0.2 MPa，但都小于圍巖及礦石的抗拉強(qiáng)度，因而采場不會出現(xiàn)拉伸破壞。隨著礦柱厚度的不斷減小，礦柱的最小主應(yīng)力呈逐漸增大的趨勢，如圖10所示。但6 m礦柱采場相較于8 m礦柱采場的最小主應(yīng)力變化不大，而4 m礦柱采場的最小主應(yīng)力卻有了較為明顯的增加，結(jié)合圖8-c）可知，4 m礦柱采場在較高的最小主應(yīng)力下發(fā)生壓剪破壞的可能性增大。

綜上分析，6 m礦柱采場與8 m礦柱采場的最大位移、最大主應(yīng)力及最小主應(yīng)力的大小基本相同、分布大致相似，沒有表現(xiàn)出明顯的差異。盡管6 m礦柱采場的最大剪切應(yīng)變率相較于8 m礦柱采場而言有了進(jìn)一步的增加，但增加的幅度不大，沒有形成貫穿的剪切帶，可以認(rèn)為6 m礦柱采場仍然處于穩(wěn)定狀態(tài)。而4 m礦柱采場的最大位移和最小主應(yīng)力相較于6 m礦柱采場有了較為明顯的不同，并且4 m礦柱采場形成貫穿的剪切帶，據(jù)此可以得出，對于40 m×15 m的采場而言，礦柱的最小安全厚度為6 m。

4.3 垂直走向布置采場礦柱最小安全厚度確定

垂直走向布置時采場的最大位移分布云圖如圖11所示，最大位移總體也呈現(xiàn)出隨礦柱厚度變小最大位移增大的規(guī)律，并且礦柱厚度越小，礦柱最大位移的增量越大，這與沿走向布置采場的位移變化規(guī)律一致。最大剪切應(yīng)變率的分布也表現(xiàn)出了由無序到規(guī)則的變化趨勢（如圖12所示），但6 m礦柱采場的剪切帶在礦柱頂?shù)滓呀?jīng)大致形成，雖未貫通，但剪切破壞的可能性很大。

采場最小主應(yīng)力、最大主應(yīng)力均隨礦柱厚度的減小呈現(xiàn)出增大的趨勢，并且增加的幅度較為明顯，如圖13、圖14所示。8 m礦柱采場的最大主應(yīng)力為0.5 MPa，但未達(dá)到礦柱或者圍巖的最大抗拉強(qiáng)度，而6 m礦柱采場的最大主應(yīng)力達(dá)到了0.6 MPa，非常接近礦柱或者圍巖的最大抗拉強(qiáng)度。結(jié)合圖12-c）分析可知，6 m礦柱采場很可能發(fā)生拉伸與剪切聯(lián)合破壞，處于不穩(wěn)定的臨界狀態(tài)。

綜上分析，垂直走向布置時，礦柱位移、最大剪切應(yīng)變率及采場的最大、最小主應(yīng)力都隨著礦柱厚度的變小而增大，但是6 m礦柱采場的頂部出現(xiàn)了較大范圍的剪切帶，并且采場周邊的最大拉應(yīng)力接近于圍巖或礦體的抗拉強(qiáng)度，在回采過程中受到其他外力擾動，很有可能發(fā)生剪切或拉伸破壞，據(jù)此可以得出，對于30 m×25 m的采場而言，礦柱的最小安全厚度為8 m。

5 結(jié) 論

1）采用巖體結(jié)構(gòu)分析系統(tǒng)對礦體頂板、側(cè)幫及下盤圍巖典型區(qū)域進(jìn)行節(jié)理裂隙掃描測量，得到采場頂板的2組控制性節(jié)理產(chǎn)狀：143°∠83°，39°∠76°；采場側(cè)幫的2組控制性節(jié)理產(chǎn)狀：143°∠83°，39°∠76°；采場上盤圍巖的3組控制性節(jié)理產(chǎn)狀：324°∠85°，111°∠74°，318°∠58°。

2）基于Mathews穩(wěn)定圖解法針對沿走向布置與垂直走向布置采場的頂板、側(cè)幫及上盤開展穩(wěn)定性分析，采場沿礦體走向布置時，礦體厚度小于20 m，采場長度控制在40 m之內(nèi)，穩(wěn)定性較好；采場垂直礦體走向布置時，礦體厚度不超過30 m，寬度小于25 m的采場，穩(wěn)定性較好。

3）通過分析礦柱的位移、最大剪切應(yīng)變率及采場的主應(yīng)力等因素，確定了2種典型采場結(jié)構(gòu)的礦柱最小安全厚度。對于沿走向布置的40 m×15 m的典型采場，礦柱最小安全厚度為6 m；對于垂直走向布置的30 m×25 m的典型采場，礦柱最小安全厚度為8 m。

［參 考 文 獻(xiàn)］

［1］王安建，高芯蕊.中國能源與重要礦產(chǎn)資源需求展望［J］.中國科學(xué)院院刊，2020，35（3）：338-344.

［2］薛小蒙，莫東旭，羅佳，等.大直徑深孔落礦工藝在低品位厚大礦床開采中的應(yīng)用［J］.礦業(yè)研究與開發(fā)，2020，40（4）：8-11.

［3］李啟月，劉愷，李夕兵.基于協(xié)同回采的深部厚大礦體分段充填采礦法［J］.工程科學(xué)學(xué)報，2016，38（11）：1 515-1 521.

［4］陳霖，黃明清，唐紹輝，等.大直徑深孔空場嗣后充填法采場結(jié)構(gòu)參數(shù)優(yōu)化及穩(wěn)定性分析［J］.金屬礦山，2022（11）：44-51.

［5］鄒平，賀超，李愛兵，等.某緩傾斜鐵礦采場結(jié)構(gòu)參數(shù)與開采順序優(yōu)化研究［J］.礦冶工程，2020，40（2）：37-42.

［6］羅先偉，韋軍，高峰，等.深部急傾斜厚大礦體采場寬度優(yōu)化［J］.金屬礦山，2018（6）：7-10.

［7］徐帥，安龍，李元輝，等.基于SOM的深埋厚大礦體采場結(jié)構(gòu)參數(shù)優(yōu)化研究［J］.采礦與安全工程學(xué)報，2015，32（6）：883-888.

［8］劉嘉偉，黃明清，譚偉.基于拓展的Mathews穩(wěn)定圖法的采場結(jié)構(gòu)參數(shù)優(yōu)化［J］.有色金屬工程，2023，13（1）：106-113.

［9］羅來和，夏國進(jìn)，王沉，等.基于FLAC3D的采礦方法優(yōu)選及采場結(jié)構(gòu)參數(shù)優(yōu)化［J］.礦冶工程，2021，41（6）：129-133.

［10］趙興東，周鑫，魏慧.三山島金礦分段空場嗣后充填采場結(jié)構(gòu)參數(shù)設(shè)計［J］.金屬礦山，2022（1）：101-106.

［11］劉建東，解聯(lián)庫，曹輝.大規(guī)模充填采礦采場穩(wěn)定性研究與結(jié)構(gòu)參數(shù)優(yōu)化［J］.金屬礦山，2018（12）：10-13.

［12］郭進(jìn)平，晏承園，張雯，等.分段礦房嗣后充填采礦采場結(jié)構(gòu)參數(shù)優(yōu)化研究［J］.礦業(yè)研究與開發(fā)，2021，41（10）：1-7.

［13］張小輝，姜春雨，程小輝.基于Q系統(tǒng)的巖體質(zhì)量快速評價優(yōu)化創(chuàng)新方法及應(yīng)用［J］.有色金屬（礦山部分），2022，74（4）：102-109.

Research on optimization of stope structure parameters of thick orebody

Yan Peng1，2 ，Chen Shaomin1，3，Huang Min1，4，Huang Zhiguo5，Zhang Bao4

（1.Zijin Mining Group Co.，Ltd.; 2.Guizhou Xinhengji Mining Industry Co.，Ltd.; 3.Fuyun Jinshan Mining & Metallurgy Co.，Ltd.; 4.Zijin （Changsha） Engineering Technology Co.，Ltd.;5.School of Earth Sciences，East China Institute of Technology）

Abstract：Reasonable structural parameters are of great significance to maintain stope stability，improving the productivity of a single stope，and reducing mining and cutting quantities.Taking Mengku Iron Mine as the research object，the rock mass telemetering and structural analysis system Sirovision are used to investigate and digitally identify the joint surface of the roof，side wall，and hanging wall of the stope，and obtain the spatial orientation statistical information and dominant joint occurrence of the stope joint group.Based on Mathews graphic method，the stability analysis is carried out for the roof，side wall and hanging wall of the stope arranged along the strike and perpendicular to strike，and the correlation between the actual hydraulic radius and the allowable hydraulic radius of the stope roof，side wall and hanging wall is obtained.The minimum safe thickness of pillar for two typical stope structures is determined by numerical analysis method through analyzing the displacement of pillar，maximum shear strain rate and main stress of stope.The analysis results show that when the stope is arranged along the strike of the ore body，the thickness of the orebody is less than 20 m，and the stope length is controlled within 40 m，the stability is good.When the stope is arranged perpendicular to the strike of the orebody，and the thickness of the ore body is not more than 30 m，the stope with a width of less than 25 m has good stability.For the 40 m×15 m typical stope arranged along the strike，the minimum safety thickness of the pillar is 6 m，and for the 30 m×25 m typical stope perpendicular to strike，the minimum safe thickness of the pillar is 8 m.The research results can provide a basis for the optimization of stope structure parameters of Mengku Iron Mine and can provide a reference for the design of stope structures of similar mines.

Keywords：thick orebody;Mathews graphic method;stope structural parameters;numerical analysis;pillar thickness

收稿日期：2023-03-25； 修回日期：2023-05-08

基金項(xiàng)目：福建省自然科學(xué)基金（2019J05039）；紫金礦業(yè)集團(tuán)股份有限公司2022年第二批科技計劃項(xiàng)目（5401KY2022100003，2674KY2022100001，2674KY2022100002）

作者簡介：嚴(yán) 鵬（1986—），男，正高級工程師，碩士，從事金屬礦山安全高效開采及技術(shù)標(biāo)準(zhǔn)化方向的研究工作;E-mail：yan.peng@zijinmining.com