楊曉炳 閆澤鵬 尹升華 楊航

摘 要：為探明粗骨料膏體的管輸阻力特性，開展了關(guān)于粗骨料膏體輸送特性的工業(yè)級 環(huán)管輸送試驗，利用白金漢方程對環(huán)管數(shù)據(jù)擬合處理后獲得膏體的實際流變參數(shù)，并建立了基于流動度測試的粗骨料膏體阻力方程.通過響應(yīng)面分析法（RSM-BBD）對膏體輸送行為進(jìn)行分析，分別得到單因素和多因素耦合對阻力損失的影響規(guī)律，并優(yōu)化了粗骨料膏體的管道 輸送參數(shù).結(jié)果表明：流變測試得到的流變參數(shù)普遍大于環(huán)管擬合參數(shù)，其中屈服應(yīng)力相差10.4%～12.2%、塑性黏度相差 21%～32.8%.通過響應(yīng)面分析得到了單因素對阻力損失影響的敏 感程度為質(zhì)量分?jǐn)?shù)>管徑>流量，并發(fā)現(xiàn)因素間的交互作用影響顯著，其中質(zhì)量分?jǐn)?shù)和管徑的交互影響起到?jīng)Q定性作用.以金川實際充填需求為背景采用響應(yīng)面優(yōu)化后的輸送參數(shù)為：管 徑 220mm、流量142.4 m3/h、質(zhì)量分?jǐn)?shù)73.2%.研究成果可為粗骨料膏體充填系統(tǒng)運行參數(shù)的選 取提供參考.

關(guān)鍵詞：粗骨料膏體;管輸阻力;環(huán)管試驗;流變特性;響應(yīng)面分析

中圖分類號：TD862.2 文獻(xiàn)標(biāo)志碼：A

Transportation Pressure Model and Optimization of Coarse Aggregate Paste Backfill Based on Pipe Loop Test

YANG Xiaobing1，2，3，YAN Zepeng1，3?，YIN Shenghua1，3，YANG Hang2

（1.School of Civil and Resources Engineering，University of Science and Technology Beijing，Beijing100083，China;

2.State Key Laboratory of Mineral Processing，Beijing102628，China;

3.Key Laboratory of High-Efficient Mining and Safety of Metal Mines of Ministry of Education（University of Science and Technology Beijing），Beijing100083，China）

Abstract：To determine the pipe pressure drop characteristics in the coarse aggregate paste backfill（CPB）pro-cess.An industrial-scale loop test on the transportation characteristics of CPB was carried out.The actual rheological parameters of CPB are obtained by fitting the Buckingham equation to the loop test data，and the resistance equationwas established based on the fluidity test.The response surface method（RSM-BBD）was used to analyze the trans-portation behavior of the CPB，and the effects of single-factor and multi-factor coupling on the pressure loss were ob-tained，respectively.Finally，RSM-BBD was used to optimize the pipeline transportation parameters of CPB.The re-sults show that the rheological parameters obtained by the rheological test were generally larger than those obtained by the loop test，where the differences in yield stress and plastic viscosity between the two are10.4%～12.2% and 21%～32.8%，respectively.Response surface analysis showed that the sensitivity of a single factor on pressure loss was as follows concentration > pipe diameter > flow rate.And it was found that the interaction between factors had a significant effect，in which the interaction between mass fraction and pipe diameter played a decisive role.Taking the actual filling demand of Jinchuan as the background，the transport parameters optimized by response surface were as follows：pipe diameter 220mm，flow rate142.4 m3/h，and concentration 73.2%.The research results are of great significance to the selection of the CPB system operating parameters.

Key words：coarse aggregate paste;pressure drop;loop test;rheological properties;response surface analysis

礦產(chǎn)資源開發(fā)利用過程中產(chǎn)生大量的廢石、尾 砂等固體廢物[1-2].據(jù)統(tǒng)計[3]，每產(chǎn)出1萬 t礦石，平均 要產(chǎn)出0.93萬 t尾礦、0.34萬 t廢石.礦山固體廢料的排放不僅占用大量的田地，還會帶來一系列環(huán)境問 題，如廢石場的泥石流、尾礦庫的滲漏和廢水對環(huán)境的污染等[4-6].經(jīng)過眾多專家學(xué)者的努力，粗骨料膏 體充填技術(shù)成為解決礦山固體廢物問題的有效方式，并在世界各地采礦實踐中越來越受重視[7-8].采用粗骨料膏體充填，通常的做法是在地面充填站將 粗骨料、尾砂、膠凝材料和水混合制成復(fù)合材料，然 后通過管道系統(tǒng)輸送至采空區(qū)[9].膏體輸送至采空區(qū)一段時間后形成具有強度的膠結(jié)充填體，這可以限制圍巖變形和提高礦石回收率.

管道輸送是粗骨料膏體充填的關(guān)鍵環(huán)節(jié)[10].由于膏體質(zhì)量分?jǐn)?shù)較高，管輸阻力大，導(dǎo)致管道輸送難度偏大.因此，有必要對粗骨料膏體的輸送特性進(jìn)行研究.流變特性是膏體輸送特性研究中的一個重要的分支，主要研究方法有L管法、傾斜管法和流變儀 測試等[11].其中L管法和傾斜管法由于流動距離、循 環(huán)時間等因素的影響導(dǎo)致測試結(jié)果存在較大偏差，流變儀測試成為比較普遍的研究方法[12].但是，相關(guān)研究表明流變儀測試所得流變參數(shù)與實際膏體參數(shù)存在較大差距.阻力損失也是表征粗骨料膏體可輸 送性能的一個重要依據(jù)[13]，國內(nèi)外對阻力損失的研究方法分為三種：經(jīng)驗公式法[14]、數(shù)值模擬法[15]和環(huán) 管試驗法[16].相比經(jīng)驗公式和數(shù)值模擬而言，環(huán)管試 驗充分考慮了各種因素對輸送的影響，所得試驗結(jié)果與工程實際的吻合度最高[17].例如：王洪江等[18]采用環(huán)管試驗法對全尾砂料漿的輸送特性進(jìn)行分析，并得到了流變參數(shù)與灰砂比和濃度相關(guān)聯(lián)的管道阻力預(yù)測方程.王勇等[19]揭示了料漿管輸實際值和理論計算值之間存在的關(guān)系.鄭伯坤等[20]研究了改性 全尾砂料漿在環(huán)管試驗條件下的輸送性能.但是，前 述研究主要側(cè)重于全尾砂漿，粗骨料膏體在物料組成與性質(zhì)上和全尾砂充填料漿存在著本質(zhì)區(qū)別，基于全尾砂料漿所得試驗結(jié)果并不能直接應(yīng)用于粗骨 料膏體.

因此，本文以金川二礦區(qū)粗骨料膏體充填系統(tǒng)為研究背景，基于工業(yè)級環(huán)管試驗系統(tǒng)對全尾砂-粗 骨料膏體的輸送特性進(jìn)行研究，分析了流變儀測試結(jié)果與環(huán)管試驗結(jié)果的區(qū)別，并建立考慮流動度的輸送阻力方程.通過響應(yīng)面設(shè)計法（Box Behnken De-sign）對粗骨料膏體的管道輸送參數(shù)進(jìn)行設(shè)計，研究了單因素和多因素耦合對阻力損失的影響，并確定 各因素對阻力損失影響程度的大小，最后對管道輸 送參數(shù)進(jìn)行優(yōu)化.研究成果可為粗骨料充填系統(tǒng)的設(shè)計提供理論依據(jù).

1試驗材料和方法

1.1試驗材料

進(jìn)行環(huán)管試驗的試驗材料均取自金川二礦區(qū)，其中細(xì)骨料為全尾砂，粗骨料為廢石和棒磨砂.

1）骨料.金川二礦區(qū)全尾砂密度為2.645 g/cm3，對其進(jìn)行XRF 元素測試分析（結(jié)果見表1），發(fā)現(xiàn)尾 礦中主要化學(xué)成分為SiO2、MgO.采用激光粒度儀測試得到全尾砂粒級分布曲線見圖1（a），尾砂中極細(xì) 顆粒（粒徑<20μm）質(zhì)量分?jǐn)?shù)為54.7%，粒徑<45 μm 顆粒質(zhì)量分?jǐn)?shù)為77.3%，粒徑<74 μm 顆粒質(zhì)量分?jǐn)?shù)為91%;計算得到尾砂的Cu 值為7.82，Cc 值為1.61，可知全尾砂級配良好.試驗所用廢石密度為2.876 g/ cm3，松散堆積密實度和密實孔隙率分別為0.582 4和31.57%.采用篩分法測試得到廢石的粒級分布曲線（見圖1（b）），可以看出廢石的主要粒徑分布集中在2～15mm，粒徑小于10mm的細(xì)顆粒質(zhì)量分?jǐn)?shù)超過80%.棒磨砂密度為2.794 g/cm3，松散堆積密實度和密實孔隙率分別為0.5576和40.23%.篩分結(jié)果表明 棒磨砂顆粒比廢石顆粒小，6 mm以下的顆粒質(zhì)量分?jǐn)?shù)為83.1%.

2）水泥.采用金川礦用非標(biāo)38.5水泥（于充填站水泥 倉處取得），測得密度為3.105 g/cm3，容 重為1.12 g/cm3.

1.2? 試驗裝置

1.2.1? 環(huán)管試驗系統(tǒng)

充填環(huán)管試驗以高可靠度的方式對充填料漿在礦山現(xiàn)場類似工況條件下的運行參數(shù)進(jìn)行模擬，為膏體充填工藝設(shè)計提供依據(jù).本研究采用的環(huán)管試 驗系統(tǒng)主要包括：配料機(jī)、計量秤、螺旋輸送機(jī)、皮帶給料機(jī)等設(shè)備，圖2給出了環(huán)管系統(tǒng)結(jié)構(gòu)示意圖.試 驗物料按照配比設(shè)計，通過配料系統(tǒng)進(jìn)行精確計量，然后進(jìn)入攪拌系統(tǒng)，同時根據(jù)設(shè)計濃度加入清水;攪 拌系統(tǒng)分為二段臥式攪拌，分別采用德國 BHS 臥式雙軸強力攪拌機(jī)和雙螺旋臥式攪拌機(jī).泵送系統(tǒng)采用德國普斯邁斯特公司（Putzmeister）的KOS1070P 型柱塞泵，額定流量60m3/h，最大工作壓力8.5mPa，泵送流量通過頻率來調(diào)節(jié).管道系統(tǒng)由3種尺寸的無縫鋼管組成.料漿在環(huán)管系統(tǒng)中完成測試后返回雙螺旋臥式攪拌機(jī)中，實現(xiàn)物料的閉路循環(huán).通過隔 膜式壓力傳感器（共12個，編號為P1～P12，未全部在示意圖中標(biāo)明）采集料漿輸送過程中管段內(nèi)的壓力值，數(shù)據(jù)采集頻率為1次/s，通過DCS 控制系統(tǒng)實時采集、記錄.

1.2.2 流變儀

選用產(chǎn)自丹麥的ICAR 混凝土流變儀對環(huán)管系統(tǒng)攪拌槽的取樣漿體進(jìn)行流變測試.ICAR的轉(zhuǎn)子形式為四片裝葉輪（轉(zhuǎn)速為0.001～0.6 r/s），轉(zhuǎn)子高12.7 cm、直徑為6.35 cm，可用來測量骨料尺寸范圍為6～32mm的粗骨料膏體.ICAR 混凝土流變儀的優(yōu)點是可以對穩(wěn)定階段的流變特性曲線進(jìn)行自動擬合，從而 計算出流變參數(shù)（初始剪切應(yīng)力τ0和黏度系數(shù)μ）.

流變測試流程如下：首先設(shè)置轉(zhuǎn)子轉(zhuǎn)速為40r/s保持100s，然后放置 200s，以消除隨機(jī)因素的干擾;然后設(shè)置轉(zhuǎn)子轉(zhuǎn)速在100s 內(nèi)從0r/s增加到150r/s，在變剪切速率的條件下獲取料漿的流變性能.

1.2.3水泥膠砂流動度儀

依據(jù)《水泥膠砂流動度測定方法》（GB/T 2419— 2005），選擇 NLD-3型水泥膠砂流動度儀（跳桌）對粗 骨料膏體的振動擴(kuò)散度進(jìn)行測試，進(jìn)而表征其流動性能.

1.3 環(huán)管試驗方案

由于環(huán)管試驗單次試驗物料消耗大，不能實現(xiàn)所有配比的環(huán)管試驗.為此，選擇具有代表性的配比進(jìn)行試驗，以達(dá)到預(yù)測相近配比膏體管輸特性的目的，試驗方案如表2所示（由于實際質(zhì)量分?jǐn)?shù)不易控 制，表2中僅給出質(zhì)量分?jǐn)?shù)的調(diào)控范圍）.環(huán)管試驗的具體步驟在文獻(xiàn)[16]和文獻(xiàn)[18]中進(jìn)行了詳細(xì)的敘述，不贅述.下面針對本次環(huán)管試驗過程中的工況參數(shù)調(diào)節(jié)過程進(jìn)行說明.

1）設(shè)備調(diào)試.系統(tǒng)開始使用前先做清水試驗，其 目的是檢測密封性、檢測儀器精度可靠性以及泵對流量的調(diào)節(jié)能力，檢測數(shù)據(jù)如圖3所示.

2）同配比料漿質(zhì)量分?jǐn)?shù)調(diào)節(jié).通過添加外部清水進(jìn)行濃度調(diào)節(jié)，試驗配料質(zhì)量分?jǐn)?shù)設(shè)定為3種，在74%～78%范圍內(nèi)調(diào)節(jié)，每改變一次濃度，必須經(jīng)過10～20min的連續(xù)循環(huán)混合，使膏體料漿混合均勻后 進(jìn)行數(shù)據(jù)讀取.

3）柱塞泵流量調(diào)節(jié).通過調(diào)節(jié)泵送控制系統(tǒng)流量，泵送頻率范圍為55%～100%，對應(yīng)理論流量范圍應(yīng)為24～58 m3/h，每種流量有效測試時間不少于5min.

2? 基于環(huán)管試驗的管輸阻力計算

2.1 環(huán)管試驗結(jié)果

經(jīng)過環(huán)管輸送試驗后，選取壓力表P8和P9所記 錄的水平管段的數(shù)據(jù)進(jìn)行處理.限于篇幅，本節(jié)內(nèi)容 僅給出 m 尾砂 ∶m 廢石 ∶m 棒磨砂為1∶0.5 ∶0.5的數(shù)據(jù)進(jìn)行分析，其余組均具有相似的變化規(guī)律，且實際質(zhì)量分?jǐn)?shù)以現(xiàn)場取樣測定為準(zhǔn)（實測質(zhì)量分?jǐn)?shù)分別為73.4%、75.2%和77.1%）.圖4為不同濃度條件下壓力數(shù)據(jù)和流量數(shù)據(jù)的變化情況，可以看出：壓力值隨著泵送流量的降低而降低，且隨著料漿濃度的增加，壓力表P8和P9所測得的壓力數(shù)據(jù)的差值越來越大.

2.2? 環(huán)管試驗結(jié)果擬合分析

由于壓力傳感器的敏感度較高，當(dāng)廢石、棒磨砂 粗顆粒撞擊傳感器觸頭時會造成數(shù)值的波動（見圖4），這將導(dǎo)致測試所得數(shù)據(jù)因受到“噪聲”的干擾而 較為離散、數(shù)據(jù)精度低.因此，對無效數(shù)據(jù)進(jìn)行剔除 后，取膏體平穩(wěn)運行段的數(shù)據(jù)進(jìn)行數(shù)據(jù)擬合，得到不同配比條件下膏體輸送阻力和流速之間的關(guān)系式，如表3所示.同時，因膏體管輸阻力與流速間呈現(xiàn)顯 著的線性關(guān)系，膏體充填料可近似視為賓漢塑性體，適用于管流的流變方程可用白金漢方程[21-22]描述：

將式（2）與表中擬合曲線方程進(jìn)行聯(lián)立求解，將 管道直徑代入后就可獲得實際條件下膏體的屈服應(yīng)力和塑性黏度，計算結(jié)果見表3.

將環(huán)管測試所得流變參數(shù)與流變測試所得結(jié)果繪制成柱狀圖進(jìn)行對比分析，如圖5所示.

圖5表明，流變儀測試所得流變參數(shù)普遍比環(huán) 管擬合參數(shù)要大，其中屈服應(yīng)力相差10.4%～12.2%、塑性黏度相差 21%～32.8%.分析認(rèn)為，在進(jìn)行流變儀 測試時采用的是槳式轉(zhuǎn)子，該轉(zhuǎn)子可有效地減少壁 面滑移效應(yīng)的發(fā)生，但是壁面滑移在實際輸送過程中無法避免.同時，由于實際輸送的時間遠(yuǎn)大于流變測試的時間，這就導(dǎo)致環(huán)管中料漿的觸變性對實際 流變參數(shù)的影響較大.由圖5 還可以看出屈服應(yīng)力和塑性黏度均隨著質(zhì)量分?jǐn)?shù)的增加而增大，且當(dāng)質(zhì)量分?jǐn)?shù)固定時料漿的屈服應(yīng)力隨著物料組分的增加而降低.這主要是因為粗骨料膏體會形成一種由全尾砂包裹粗顆粒的穩(wěn)定漿體.而由粗骨料粒徑分析可知，棒磨砂6 mm以下的顆粒質(zhì)量分?jǐn)?shù)為83.1%，相對廢石而言是一種細(xì)骨料.因此，當(dāng)骨料為三組分時，廢石和棒磨砂的級配更為合理，全尾砂漿體對二 者的懸浮包裹效果更佳，減少了顆粒間的運動阻力，表現(xiàn)為粗骨料膏體的屈服應(yīng)力減小.

式（1）表明傳統(tǒng)的阻力計算通常根據(jù)測出的流變參數(shù)并結(jié)合相應(yīng)的輸送參數(shù)（如管徑、流速等）進(jìn)行計算.然而礦山并不具備相應(yīng)測試條件，通常委托 試驗機(jī)構(gòu)進(jìn)行流變參數(shù)的獲取[23].但是圖5結(jié)果表明實驗室所得參數(shù)和實際參數(shù)存在較大的誤差.為了更加有效、準(zhǔn)確且便捷地評估粗骨料膏體的輸送 性能，在傳統(tǒng)的阻力計算公式中引入相對容易獲得的流動度參數(shù).將膏體的流動度值與屈服應(yīng)力、塑性黏度進(jìn)行擬合，進(jìn)而得到基于流動度的阻力計算新模型，流變參數(shù)擬合結(jié)果如圖6所示（R2 均大于0.9）.將圖6中的擬合方程代入式（1）即可得到不同骨料條件下的輸送阻力經(jīng)驗公式，如式（3）～式（5）.

3 管道輸送參數(shù)優(yōu)化分析

3.1? 響應(yīng)面法試驗水平設(shè)計

本節(jié)內(nèi)容以全尾砂-廢石-棒磨砂膏體為主要研究對象，選擇管道直徑、輸送流量和固體質(zhì)量分?jǐn)?shù)作為自變量，輸送阻力損失作為響應(yīng)值，結(jié)合礦山實際 充填需求確定本次試驗中管徑及輸送流量的取值范圍，管道直徑為140mm、180mm、220mm;輸送流量為110m3/h、130m3/h、150m3/h，為了便于計算膏體的質(zhì)量分?jǐn)?shù)，根據(jù)環(huán)管試驗時的參數(shù)進(jìn)行取整，膏體質(zhì)量分?jǐn)?shù)分別為73%、75%、77%.最后，采用Design Expert 軟件對方案進(jìn)行設(shè)計，進(jìn)而分析多因素交互 作用對阻力損失的影響，設(shè)計后的因素水平及編碼 如表4所示.試驗方案如表5所示，將試驗所對應(yīng)參數(shù)值代入式（5）中進(jìn)行計算，分別得到不同條件下的輸送阻力.

Design Export對表5中數(shù)據(jù)進(jìn)行多元回歸擬合后得到響應(yīng)面函數(shù)，如式（6）所示.對所建立的響應(yīng) 面回歸模型進(jìn)行可靠性檢驗（R2 檢驗），得到每個影 響因子的顯著程度（P 值），發(fā)現(xiàn)所建立響應(yīng)面模型的R2為0.9973且 P 值小于0.0001，說明響應(yīng)面回歸模型的回歸效果極好、可靠性高.

3.2 單因素對膏體輸送阻力的影響

圖7（a）為膏體質(zhì)量分?jǐn)?shù)為73%時管道直徑與阻力損失的關(guān)系.由圖7（a）可知，當(dāng)流量和質(zhì)量分?jǐn)?shù)不變時，隨著管道直徑的增大阻力損失呈下降的趨勢，但是下降幅度在逐漸降低.這是因為管徑增大時，管 壁剪切作用力減小，邊界層的滑移效果減弱，此時阻力損失主要受到管內(nèi)顆粒間摩擦作用的影響.管徑增大時各流層間的距離增大[22]，顆粒的碰撞行為減 少，導(dǎo)致阻力損失快速減小;但是管徑達(dá)到220mm 后臨界速度發(fā)生變化，料漿流動形態(tài)由層流向紊流 轉(zhuǎn)變，顆粒間的碰撞行為加劇，導(dǎo)致阻力損失的降低 幅度變緩.同時，結(jié)合表6發(fā)現(xiàn)管道直徑 X1的F 值較大，說明管徑對阻力損失的影響較為顯著.圖7（b）為膏體質(zhì)量分?jǐn)?shù)為73%時管道流量與阻力損失的關(guān)系.由圖7（b）可知，當(dāng)管徑和質(zhì)量分?jǐn)?shù)不變時，隨著 流量的增大，阻力損失呈上升的趨勢.同時，結(jié)合表6發(fā)現(xiàn)，管道流量的F 值最小，說明其對阻力損失的影響程度最低.流速的增大導(dǎo)致膏體物料所承受的作用力發(fā)生變化，極大地提高了顆粒間的碰撞概率，相應(yīng)的阻力損失增大.但是圖7（b）曲線增長緩慢和F 值均表明流量對阻力損失的影響明顯不如管徑及膏體質(zhì)量分?jǐn)?shù).圖7（c）為管徑110mm時質(zhì)量分?jǐn)?shù)與阻力損失的關(guān)系.由圖7（c）可知，膏體質(zhì)量分?jǐn)?shù)與阻力損失呈正相關(guān)關(guān)系，且隨著質(zhì)量分?jǐn)?shù)的增加，阻力損失的增幅變大.分析認(rèn)為，膏體中的水分起到潤滑的效果，當(dāng)質(zhì)量分?jǐn)?shù)增加時相應(yīng)的膏體含水量降低，導(dǎo)致漿體在流經(jīng)管道時所受的摩擦增大，從而提高了膏體輸送的阻力損失.在滿足膏體可輸送性能的同時適當(dāng)提高質(zhì)量分?jǐn)?shù)可以降低水泥的用量，有助于降低生產(chǎn)成本.同時，結(jié)合表6發(fā)現(xiàn)，質(zhì)量分?jǐn)?shù)的F 值最大，說明其對阻力損失的影響程度最大，最終得到管徑、流量和質(zhì)量分?jǐn)?shù)對阻力損失的敏感程度為：質(zhì)量分?jǐn)?shù)>管徑>流量.

3.3 因素交互作用對輸送阻力的影響

結(jié)合表6中模型參數(shù)的方差分析結(jié)果發(fā)現(xiàn)，X1X3的F 值最大，說明管徑和質(zhì)量分?jǐn)?shù)的交互作用對膏 體阻力損失的影響最為顯著.圖8（a）為膏體質(zhì)量分?jǐn)?shù)為73%時，管徑與流量間的交互作用對膏體管道 輸送阻力損失的影響.可以看出，當(dāng)流量較低（110m3/h）時，管徑由140mm增大至220mm，阻力損失降 低了54.9%，這主要是由于各流層間的距離增大所致;當(dāng)管徑為180mm，流量從110m3/h增大至150m3/h時，阻力損失提升了9.35%，可以看出隨著管徑的增大，阻力損失對流量的敏感度逐漸降低，這主要是由于流量提升導(dǎo)致顆粒碰撞的效果小于流層間距增大導(dǎo)致的減阻效果.綜上，在一定范圍內(nèi)同時增加管徑和流量有助于礦山充填在提升產(chǎn)能的同時降低 損耗.圖8（b）為流量110m3/h時，管徑和質(zhì)量分?jǐn)?shù)的交互作用對膏體管道輸送阻力損失的影響.圖8（b）表明當(dāng)管徑為140mm，膏體質(zhì)量分?jǐn)?shù)由 73%增加至77%時，阻力損失提高了163.2%;而當(dāng)質(zhì)量分?jǐn)?shù)為73%時，管徑由140mm增大至220mm，阻力損失下降了55.1%.同時可以看出，阻力損失對管徑的敏感度隨著質(zhì)量分?jǐn)?shù)的增加而降低.這主要是因為膏體的流動性能主要依賴于其內(nèi)部含水量，質(zhì)量分?jǐn)?shù)的增大導(dǎo)致含水量的減少，抵消了管徑增大所帶來的減阻效果.圖8（c）為管徑180mm時，流量和質(zhì)量分?jǐn)?shù)的交互作用對膏體管道輸送阻力損失的影響.圖8（c）表明當(dāng)流量為110m3/h，膏體質(zhì)量分?jǐn)?shù)由 73%增加至77%時，阻力損失上漲了177.9%;而當(dāng)質(zhì)量分?jǐn)?shù)為73%時，流量由110m3/h增大至150m3/h，阻力損失增加了9.89%.由于流量增大時會導(dǎo)致膏體內(nèi)部 顆粒的紊動進(jìn)而提高輸送阻力，而質(zhì)量分?jǐn)?shù)的提高導(dǎo)致膏體內(nèi)部固體顆粒增多，因此阻力損失對流量的敏感度隨著質(zhì)量分?jǐn)?shù)的增加而變得顯著.

3.4 響應(yīng)面模型對輸送參數(shù)的優(yōu)化分析

金川全尾砂-廢石膏體充填站設(shè)計充填能力200萬 m3/a（兩套充填管路），其中單套充填系統(tǒng)承擔(dān)能力100萬 m3/a，經(jīng)換算后每小時的充填能力不應(yīng)小于115m3，將其代入響應(yīng)面模型進(jìn)行優(yōu)化設(shè)計，響應(yīng)面 優(yōu)化結(jié)果如圖9和圖10所示.可以看出，在管道參數(shù)為管徑 220mm、流量142.4 m3/h、質(zhì)量分?jǐn)?shù)73.2%時，管道阻力損失最低為1.939kPa/m，此結(jié)果的可信度為0.996（最大值為1）.

4結(jié)論

1）利用白金漢方程對環(huán)管試驗數(shù)據(jù)擬合后得到膏體的實際流變參數(shù)，發(fā)現(xiàn)流變儀測試所得流變參數(shù)普遍比環(huán)管擬合參數(shù)要大，其中屈服應(yīng)力相差10.4%～12.2%、塑性黏度相差 21%～32.8%，這是由實 際輸送過程中壁面滑移和觸變性的影響所致.

2）由于流變測試所得參數(shù)和實際參數(shù)存在較大誤差，在傳統(tǒng)阻力計算公式中引入相對容易獲得的流動度參數(shù)，建立了基于流動度的阻力計算新模型.

3）基于RSM-BBD開展全尾砂-粗骨料膏體輸送阻力行為分析，研究了質(zhì)量分?jǐn)?shù)、管徑和流量及其交 互作用條件下管道輸送阻力的變化規(guī)律.得到了單 因素對阻力損失影響的敏感度大小為：質(zhì)量分?jǐn)?shù)>管 徑>流量，并發(fā)現(xiàn)因素間的交互作用影響顯著，其中質(zhì)量分?jǐn)?shù)和管徑的交互影響起決定性的作用.

4）以金川二礦區(qū)粗骨料膏體充填系統(tǒng)為背景，采用響應(yīng)面對輸送參數(shù)進(jìn)行優(yōu)化，最終得到最佳的輸送參數(shù)為：管徑 220mm、流量142.4 m3/h、質(zhì)量分?jǐn)?shù)73.2%.該成果對于粗骨料膏體充填系統(tǒng)運行參數(shù)的選取具有一定的指導(dǎo)意義.

參考文獻(xiàn)

[1]李夕兵，周健，王少鋒，等.深部固體資源開采評述與探索[J].中國有色金屬學(xué)報，2017，27（6）：1236-1262.

LI X B，ZHOU J，WANG S F，et al.Review and practice of deep mining for solid mineral resources [J].The Chinese Journal of Nonferrous Metals，2017，27（6）：1236-1262.（In Chinese）

[2]YANG L，YILMAZ E，LI J W，et al.Effect of superplasticizer

type and dosage on fluidity and strength behavior of cemented tail-ings backfill with different solid contents [J].Construction and Building Materials，2018，187：290-298.

[3]黃志偉，古德生.我國礦山無廢開采的現(xiàn)狀[J].礦業(yè)研究與開發(fā)，2002，22（4）：9-10.

HUANG Z W，GU D S.The state-of arts of wasteless mining in China[J].Mining Research and Development，2002，22（4）：9-10.（In Chinese）

[4]尹升華，劉家明，陳威，等.不同粗骨料對膏體凝結(jié)性能的影響及配比優(yōu)化[J].工程科學(xué)學(xué)報，2020，42（7）：829-837.

YIN S H，LIU J M，CHEN W，et al.Optimization of the effect and formulation of different coarse aggregates on performance of the paste backfill condensation[J].Chinese Journal of Engineering，2020，42（7）：829-837.（In Chinese）

[5]BENZAAZOUA M，BUSSI?RE B，DEMERS I，et al.Integrated

mine tailings management by combining environmental desulphu-rization and cemented paste backfill：application to mine Doyon，Quebec，Canada [J]. Minerals Engineering，2008，21（4） ：330-340.

[6]楊曉炳，肖柏林，高謙.非牛頓流體充填料漿的管輸摩擦阻力預(yù)測[J].湖南大學(xué)學(xué)報（自然科學(xué)版），2020，47（6）：125-131.YANG X B，XIAO B L，GAO Q.Prediction of friction pressure for non-Newtonian backfilling pipe-flow slurry [J].Journal of Hunan University（Natural Sciences），2020，47（6）：125-131.（In Chinese）

[7]楊天雨，喬登攀，王俊，等.廢石-風(fēng)砂高濃度料漿管道輸送數(shù)值模擬及管輸阻力新模型[J].中國有色金屬學(xué)報，2021，31（1）：234-244.

YANG T Y，QIAO D P，WANG J，et al.Numerical simulation and new model of pipeline transportation resistance of waste rock-aeolian sand high concentration slurry[J].The Chinese Journal of Nonferrous Metals，2021，31（1）：234-244.（In Chinese）

[8]甘德清，閆澤鵬，薛振林，等.考慮壁面滑移效應(yīng)的充填料漿管道輸送阻力研究[J].金屬礦山，2020（9）：26-32.

GAN D Q，YAN Z P，XUE Z L，et al.Study on the resistance of cemented paste backfill slurry pipeline transport considering wall slip effect[J].Metal Mine，2020（9）：26-32.（In Chinese）

[9]WU D，YANG B G，LIU Y C.Transportability and pressure drop

of fresh cemented coal gangue-fly ash backfill（CGFB）slurry in pipe loop[J].Powder Technology，2015，284：218-224.

[10]顏丙恒，李翠平，吳愛祥，等.膏體料漿管道輸送中粗骨料顆粒運動規(guī)律分析[J].中南大學(xué)學(xué)報（自然科學(xué)版），2019，50（1）：172-179.

YAN B H，LI C P，WU A X，et al.Analysis of law of movement of coarse aggregate particles in pipeline transportation of paste[J].Journal of Central South University（Science and Technology），2019，50（1）：172-179.（In Chinese）

[11]胡華，孫恒虎，黃玉誠，等.似膏體粘彈塑性流變模型與流變方程研究[J].中國礦業(yè)大學(xué)學(xué)報，2003，32（2）：119-122.

HU H，SUN H H，HUANG Y C，et al.Rheological model and equation of viscoelastic-plasticity of paste-like backfill pulp[J].Journal of China University of Mining & Technology，2003，32（2）：119-122.（In Chinese）

[12]苗苗，雪凱旺，苗芳，等.石灰石粉對水泥漿體水化特性及流變性能的影響[J].湖南大學(xué)學(xué)報（自然科學(xué)版），2018，45（12）：90-96.

MIAO M，XUE K W，MIAO F，et al.Influence of limestone pow-der on hydration characteristics and rheological properties of ce-ment paste[J].Journal of Hunan University（Natural Sciences），2018，45（12）：90-96.（In Chinese）

[13]張修香，喬登攀 .粗骨料高濃度充填料漿的管道輸送模擬及試驗[J].中國有色金屬學(xué)報，2015，25（1）：258-266.

ZHANG X X，QIAO D P.Simulation and experiment of pipeline transportation of high density filling slurry with coarse aggregates[J].The Chinese Journal of Nonferrous Metals，2015，25（1）： 258-266.（In Chinese）

[14]陳秋松，張欽禮，王新民，等.全尾砂似膏體管輸水力坡度計算模型研究[J].中國礦業(yè)大學(xué)學(xué)報，2016，45（5）：901-906.?? CHEN Q S，ZHANG Q L，WANG X M，et al.Pipeline hydraulic gradient model of paste-like unclassified tailings backfill slurry[J].Journal of China University of Mining & Technology，2016，45（5）：901-906.（In Chinese）

[15]WU A X，RUAN Z E，WANG Y M，et al.Simulation of long-distance pipeline transportation properties of whole-tailings paste with high sliming[J].Journal of Central South University，2018，25（1）：141-150.

[16]王少勇，吳愛祥，阮竹恩，等.基于環(huán)管實驗的膏體流變特性及影響因素[J].中南大學(xué)學(xué)報（自然科學(xué)版），2018，49（10）： 2519-2525.

WANG S Y，WU A X，RUAN Z E，et al.Rheological properties of paste slurry and influence factors based on pipe loop test [J].Journal of Central South University（Science and Technology），2018，49（10）：2519-2525.（In Chinese）

[17]YANG X B，XIAO B L，GAO Q，et al.Determining the pressure

drop of cemented Gobi sand and tailings paste backfill in a pipe flow [J]. Construction and Building Materials，2020，255：119371.

[18]王洪江，王小林，寇云鵬，等.全尾砂高濃度膠結(jié)充填的環(huán)管試驗[J].工程科學(xué)學(xué)報，2021，43（2）：215-222.

WANG H J，WANG X L，KOU Y P，et al.Loop test study on the high-concentration cemented filling of full tailings [J].Chinese Journal of Engineering，2021，43（2）：215-222.（In Chinese）

[19] 王勇，那慶，吳愛祥，等.粉煤灰-爐渣漿體環(huán)管實驗和理論計算阻力對比研究）[J/OL].中南大學(xué)學(xué)報（自然科學(xué)版）：1-11（2021-07-26）[2021-10-12].

https：//kns.cnki.net/kcms/detail/ 43.1426.N.20210723.1649.002.html.

WANG Y，NA Q，WU A X，et al.Comparison of resistance be-tween pipe loop test and theoretical calculations in fly ash?slag slurry [J/OL].

Journal of Central South University（Science and Technology）：1-11（2021-07-26）[2021-10-12].

https：//kns.cnki.net/kcms/detail/43.1426.N.20210723.1649.002.html.（In Chinese）

[20]鄭伯坤，姚維，黃騰龍，等.基于環(huán)管試驗的改性全尾砂充填料漿輸送性能[J].中國有色金屬學(xué)報，2021，31（2）：520-529. ZHENG B K，YAO W，HUANG T L，et al.Transportable perfor-mance of modified unclassed tailings filling slurry based on loop test[J].The Chinese Journal of Nonferrous Metals，2021，31（2）： 520-529.（In Chinese）

[21]吳愛祥，劉曉輝，王洪江，等.考慮時變性的全尾膏體管輸阻力計算[J].中國礦業(yè)大學(xué)學(xué)報，2013，42（5）：736-740.

WU A X，LIU X H，WANG H J，et al.Calculation of resistance in total tailings paste piping transportation based on time-varying be-havior[J].Journal of China University of Mining & Technology，2013，42（5）：736-740.（In Chinese）

[22]CHEN Q S，ZHANG Q L，WANG X M，et al.A hydraulic gradi-ent model of paste-like crude tailings backfill slurry transported by a pipeline system[J].Environmental Earth Sciences，2016，75（14）：1-9.

[23]侯永強，尹升華，戴超群，等.尾礦膏體流變特性和管輸阻力計算模型[J].中國有色金屬學(xué)報，2021，31（2）：510-519.

HOU Y Q，YIN S H，DAI C Q，et al.Rheological properties and pipeline resistance calculation model in tailings paste [J].The Chinese Journal of Nonferrous Metals，2021，31（2）：510-519.（In Chinese）