葉進霞

(山東省魯南地質(zhì)工程勘察院，山東 兗州 272100)

安全穩(wěn)定性約束下的排土場優(yōu)化設(shè)計

葉進霞

(山東省魯南地質(zhì)工程勘察院，山東 兗州 272100)

在保證露天礦排土場邊坡安全穩(wěn)定的情況下，優(yōu)化排土場幾何參數(shù)，達到在有限排土場范圍內(nèi)排土量最大的目的。根據(jù)排放物的物理性質(zhì)，選擇邊坡平面滑動和圓弧滑動的安全穩(wěn)定性系數(shù)為安全評價指標(biāo)，利用排土場底面內(nèi)接半徑把排土場的邊坡傾角、臺階高度與寬度聯(lián)系到一起，分析三個指標(biāo)之間的制約關(guān)系，將臺階寬度轉(zhuǎn)化到臺階高度和邊坡傾角的函數(shù)，分別對邊坡傾角和臺階高度進行分段求取，多次模擬計算，得出一組接近于最優(yōu)解的值作為排土場設(shè)計的最佳模型。實驗結(jié)果表明，和其他方法比較，本文算法考慮了安全穩(wěn)定性約束條件下，實現(xiàn)排土場設(shè)計盡量多排土、少占地的目的。

安全穩(wěn)定性系數(shù)；分段求解；最優(yōu)模型；排土場容積

露天采煤的第一個任務(wù)是剝離表土，需要剝離煤層上部的巖土，排棄大量的巖石和土壤，尤其較深的露天礦，往往占用較多的農(nóng)田。排土場應(yīng)選擇在靠近采礦場，盡量回填采空區(qū)，少占農(nóng)田，有條件的應(yīng)放置在山谷、洼地處，注意環(huán)境保護和造田、還田。排土場占用大量的農(nóng)田，不僅破壞環(huán)境，而且需要賠付耕地補償，很多礦山企業(yè)都對排土場進行了復(fù)墾，但是排土場的結(jié)構(gòu)并不適合種植，投入大量的資金需要改善排土場，并且復(fù)墾土地的經(jīng)濟效益不如原有土地[1]。

目前，對排土場優(yōu)化設(shè)計的相關(guān)研究中，任高峰等提出了優(yōu)化設(shè)計露天礦的邊坡空間形狀，最大限度的減少巖土剝離量，以達到減少占地的目的[2]；郭成等提出了基于安全穩(wěn)定性的大容積排土場堆置要素優(yōu)化，但是沒有將給出具體的體積計算方法[3]；衛(wèi)博等提出了基于圓錐模型的節(jié)地復(fù)墾排土場設(shè)計模型，但是沒有考慮排土場的安全穩(wěn)定性，只分析了引起排土場面積和容積的變化關(guān)系，沒有提出確定最佳模型的方法[4]。因此，本文在排土場各種物理屬性確定的情況下，分析人為控制的排土場因素與排土場邊坡安全穩(wěn)定性、排土場容積之間的關(guān)系，保證排土場邊坡安全穩(wěn)定性的情況下[5-8]，利用分段模擬法，使設(shè)計范圍固定的排土場的最大容積，減少排土場占用土。

1 排土場邊坡穩(wěn)定性分析

排土場穩(wěn)定性是影響排土場占地面積的關(guān)鍵因素，排土場參數(shù)的選擇關(guān)系到排土場設(shè)計的合理性。影響露天礦排土場邊坡穩(wěn)定性的因素很多，排土場邊坡滑動的類型主要包括平面滑動和圓弧滑動。

1.1 無確定拉張裂縫的平面滑動穩(wěn)定性系數(shù)

邊坡破壞時，滑動面為平面或近似平面時的滑動為平面滑動。露天礦排土場的土質(zhì)相對較均勻，沒有確定的滑面，也沒有確定的拉張裂縫，通過分析可以求得排土場的最危險的滑面和最可能出現(xiàn)拉張裂縫的位置[9-10]。圖1(a)為平面滑動示意圖。

平面法分析邊坡穩(wěn)定性時，滑面傾角小于邊坡傾角，滑體在受到自身重力作用下，沿著滑面向下滑動，邊坡穩(wěn)定性系數(shù)計算公式見式(1)。

(1)

式中：Fs平為邊坡穩(wěn)定性系數(shù)；γ為排棄物的容重；C為滑面上的單位粘聚力；φ為內(nèi)摩擦角；H為排土場臺階高度；Z為垂直裂面高度；α為排土場邊坡傾角；β為滑面傾角。其中γ、C、φ由排土場排棄物的物理屬性決定，最危險滑面的傾角β=(α+φ)/2，H和α根據(jù)實際情況設(shè)計，根據(jù)Fs平的極小值可求出拉張裂縫的臨界高度Z。

根據(jù)圖1(a)可知，在臺階上部最可能出現(xiàn)平面滑動的位置到臺階邊緣的距離D1=H(cosβ-cosα)，D1為平面滑動最危險距離，在設(shè)計排土臺階時，臺階寬度應(yīng)大于最危險距離D1。

1.2 圓弧滑動的穩(wěn)定性系數(shù)

圓弧滑動是指排棄物的滑面呈圓滑形狀的巖土滑動。圓弧滑動的基本假設(shè)條件：均質(zhì)粘性土坡滑動時，其滑動面常近似為圓弧形狀，假定滑動面以上的土體為剛性體，即設(shè)計中不考慮滑動土體內(nèi)部的相互作用力，假定土坡穩(wěn)定屬于平面應(yīng)變問題。圖1(b)為圓弧滑動示意圖。

圓弧滑動的穩(wěn)定性系數(shù)的基本公式為：取圓弧滑動面以上滑動體為脫離體，土體繞圓心O下滑的滑動力矩為Ms=Wa，阻止土體滑動的力是滑弧AFD上的抗滑力，其值等于土的抗剪強度與滑弧長度的乘積，穩(wěn)定系數(shù)Fs圓=抗滑力矩/滑動力矩，見式(2)。

(2)

(3)

滑動力矩由滑動圓弧內(nèi)的排棄物的作用力引起，滑體ABDF為三角體ABD和弓形體AEDF之和，而弓形體AEDF為扇形體AODF與三角體AOD之差，其各自的力矩見式(4)～(7)。則滑動力矩可表達為式(8)。

(4)

(5)

(6)

(7)

(8)

式中：H為排土場臺階高度；C為沿弧面的粘聚力；ω為弧心角之半；α為排土場邊坡傾角；β為弧弦的傾角；R為圓弧滑面半徑。其圓弧滑面半徑的計算公式見式(9)。沿弧面的粘聚力的計算公式見式(10)。

(9)

(10)

2 排土場容積推導(dǎo)

在排土車或者輸送皮帶將排棄物運送到排土場并直接傾倒在排土場時，排棄物在重力的作用下，根據(jù)排棄物的顆粒大小以及密度等性質(zhì)，排土場各個方向的邊緣都接近于圓弧形狀，臺階式排土場的各個臺階也就近似圓錐臺[11]，如圖2所示。

圖2 臺階設(shè)計示意圖

由于排土場的基底是不規(guī)則的面，處于排土場底部最高點以下的容積不能用視為圓錐臺，處于基底最高點下部的體積用不規(guī)則三角網(wǎng)TIN即可計算得到?；矫嬉陨系牟糠置恳患壟_階分為一個圓錐臺體，每一個臺階柱體的體積計算公式見式(11)。

(11)

式中：Si下和分別為第i個臺階上下面的面積；hi下和hi上分別表示第i個臺階上下面所在圓錐臺的高。由于排土場的各個臺階的側(cè)邊緣接近于圓弧形狀，再根據(jù)三角形相似性的性質(zhì)，排土場臺階的上下兩個面的面積滿足式(12)、式(13)。

(12)

(13)

式中：S上、S下分別代表排土場臺階上下面的面積；r為臺階下底面的內(nèi)接半徑；d為臺階寬度；H為臺階高度。上下兩個臺階的下底面的內(nèi)接半徑關(guān)系見式(14)。將式(12)和式(13)代入式(11)得到這個排土場的容積，見式(15)。

ri+1=ri-Hictgα-d

(14)

(15)

式中：V0表示基平面下部的容積，利用排土場底面的三維點構(gòu)建的TIN求得。

3 最佳模型模擬設(shè)計

排土場設(shè)計時，排棄物的物理性質(zhì)由礦區(qū)地質(zhì)環(huán)境決定，因此要提高排土場安全性和排土場容積只能調(diào)整排土場邊坡傾角α、臺階高度Hi與寬度di，結(jié)合排土場實際情況和排棄物的性質(zhì)，協(xié)調(diào)三個量使排土場保證安全的情況下達到最大的容積，減少排土場占地面積。排土場平面滑動和圓弧滑動的安全穩(wěn)定系數(shù)以及排土場容積的計算公式都涉及到α、Hi、di三個變量，各個排土臺階的底面內(nèi)接半徑ri也由這三個變量決定，其余各個量由排棄物的物理性質(zhì)或排土場的環(huán)境決定，不能人為地改變。選定排土場范圍之后，排土場基底面的內(nèi)接半徑r0底確定，確定邊坡傾角和臺階高度之后，平面滑動和圓弧滑動的最可能出現(xiàn)滑動的位置都相應(yīng)確定，為了削弱上下臺階之間的相互影響，設(shè)置臺階寬度為兩種滑動中的最危險距離，即d=max(D1,D2)，由此減少一個變量。

在兩個安全穩(wěn)定性的約束下，α、Hi增加到一定值時，受到安全因素的影響，繼續(xù)增大α、Hi的值，臺階級數(shù)必然減少，V的值也減小。由于Fs平、Fs圓、V與變量α、Hi的關(guān)系都比較復(fù)雜，無法用常規(guī)的方法進行規(guī)劃并求V的最大值，因此采用分段求最值的方法，確定一個值的情況下，求另一個值使得V最大，分別對α、Hi進行兩次分段求V的最大值，求得的α、Hi使V的值接近于最大值。

4 驗證實驗

研究區(qū)域范圍的面積為60492m2，內(nèi)接半徑為437m，區(qū)域內(nèi)的最高點和最低點的高程分別為478m和418m，排土場基底最高點以下部分的容積為72.39萬m3。由于排棄物的顆粒大小、質(zhì)地等不一樣，選擇各個參數(shù)的平均值作為整個區(qū)域的參數(shù)值。其中，內(nèi)摩擦角為46°，容重1.75kN/m3，單位粘聚力為1.50MPa。

邊坡傾角從0°到45°，以3°為間隔，保證平面滑動穩(wěn)定性系數(shù)和圓弧滑動穩(wěn)定性系數(shù)都大于1.2的情況下，分別計算得到對應(yīng)的臺階高度與寬度、臺階級數(shù)、排土場容積，計算結(jié)果見表1。

從表2 中可以看出，在所有的模擬值中，當(dāng)邊坡傾角為33°，對應(yīng)的臺階高度為35.47m，該排土場的容積最大，可知該區(qū)域的最佳模型的邊坡傾角在33°附近。因此，在臺階高度的基礎(chǔ)上，對33°的左右兩個分段區(qū)間再進行分段計算，邊坡傾角30°和36°對應(yīng)的最佳臺階高度為37.59m和33.58m，則從33.5m到38m按照0.5m進行分段，計算結(jié)果見表2。

按照分段計算得到臺階高度為36m以及臺階寬度為32m時，對應(yīng)邊坡傾角為32.53°，研究區(qū)域的排土容積最大，即該區(qū)域內(nèi)排土臺階設(shè)計的最佳模型為：當(dāng)邊坡傾角32.53°，臺階高度36m，臺階寬度32m，臺階級數(shù)4級時，基平面上方的容積最大，為134.73萬m3。

表1 按邊坡傾角的分段結(jié)果

表2 按臺階高度的分段結(jié)果

5 結(jié) 論

本文提出的模型是將排土場三維表面模擬為圓錐形，對于一些特殊情況可以把排土場表面簡化近視為圓錐。在保證邊坡平面滑動和圓弧滑動的安全穩(wěn)定性系數(shù)大于1.2的情況下，協(xié)調(diào)邊坡傾角、臺階高度和寬度三者之間的關(guān)系，使得給定范圍內(nèi)的排土場內(nèi)排土容積達到最大值，減少排土場占用土地面積。臺階式排土場基底面的內(nèi)接半徑將邊坡傾角、臺階高度、臺階寬度以及排土場高度聯(lián)系到了一起，形成彼此之間相互限制的關(guān)系，將最有模型轉(zhuǎn)化到邊坡傾角和臺階高度的二元函數(shù)之后，兩次分段查找體積最大的點，該點對應(yīng)的三個參數(shù)值接近于排土場設(shè)計的理想模型，可認為是臺階式排土場設(shè)計的最佳模型。在實際操作過程中，很難控制三個參數(shù)和設(shè)計值完全一致，因此保證三個參數(shù)值在一定容差范圍內(nèi)即可。

[1] 馬從安,才慶祥,韓可琦.露天礦生產(chǎn)與生態(tài)重建適宜性評價專家系統(tǒng)[J].中國礦業(yè)大學(xué)學(xué)報,2006,35(2):231-234.

[2] 任高峰，張世雄.深凹露天礦邊坡空間形狀優(yōu)化設(shè)計[J].采礦與安全工程學(xué)報,2006,23(4):472-475.

[3] 郭成,馬萃林.露天礦排土場堆置要素的優(yōu)化[J].金屬礦山,2011(1):42-45.

[4] 衛(wèi)博,付梅臣,白中科,等.基于節(jié)地的露天礦排土場設(shè)計[J].農(nóng)業(yè)工程學(xué)報,2006,22(6):230-232.

[5] 藍航,李鳳明,姚建國.露天煤礦排土場邊坡下采動沉陷規(guī)律研究[J].中國礦業(yè)大學(xué)學(xué)報,2007,36(4):482-486.

[6] 石建勛,劉新榮,廖紹波,等.礦區(qū)排土場堆載對邊坡穩(wěn)定性影響的分析[J].采礦與安全工程學(xué)報,2011,28(2)：258-264.

[7] 李小春，任偉，王少泉，等.論金屬礦山排土場設(shè)計規(guī)范中邊坡極限平衡計算方法的選取[J].巖石力學(xué)與工程學(xué)報，2011，30(S2)：4136-4142.

[8] 汪海濱，李小春，米子軍，等.排土場空間效應(yīng)及其穩(wěn)定性評價方法研究[J].巖石力學(xué)與工程學(xué)報，2011，30(10)：2103-2111.

[9] 蔣斌松,蔡美峰,都浩.平面滑動邊坡穩(wěn)定性的解析計算[J].巖石力學(xué)與工程學(xué)報,2004,23(1)：91-94.

[10] 吳侃,李亮,敖建鋒,等.開采沉陷引起地表土體裂縫極限深度探討[J].煤炭科學(xué)技術(shù),2010,38(6):108-112.

[11] 蔡利平,李鋼,史文中.增地節(jié)地型露天礦排土場優(yōu)化設(shè)計[J].煤炭學(xué)報，2013，38(12)：2208-2214.

Optimized design of dump based on the security and stability

YE Jin-xia

(Shandong Lunan Geo-engineering Exploration Institute，Yanzhou 272100，China)

The dump geometric parameters were optimized to achieve the maximum volume of dump within the limited range under the condition of the safety and stability of the open-pit dump.The security and stability factor of plane and circular sliding were selected to evaluate the safety based on the physical properties of emissions.Relationship between three indexes，slope angle，step height and width，were analyzed，followed by the computation of step height and slope inclination segment by segment.The values close to the optimal solutions were obtained via simulations.The experiment shows that，when the method takes into account the Security and Stability Constraint，the proposed method can pile up more solid waste and occupy less land than other methods.

security and stability factor；segmented solving；optimal model；dump volume

2014-04-03

葉進霞(1977-)，男，山東東平人，工程師，主要從事水文地質(zhì)、工程地質(zhì)、環(huán)境地質(zhì)工作。E-mail：yejinxia1977@163.com。

TD 216

A

1004-4051(2015)05-0153-04