王少勇 ，吳愛祥 ，王洪江 ，尹升華 ，顧曉春

(1.北京科技大學(xué) 土木與環(huán)境工程學(xué)院，北京 100083；2.北京科技大學(xué) 金屬礦山高效開采與安全教育部重點實驗室，北京 100083；3.云南銅業(yè)集團公司，昆明 650051)

銅資源現(xiàn)狀為三多三少，即貧礦多、富礦少；共生、伴生礦多，單一礦少；難處理礦多，易處理礦少。品位低、復(fù)雜難處理是我國銅資源特點。由此可見，大力開發(fā)處理低品位、復(fù)雜難處理銅礦技術(shù)對增加銅產(chǎn)量、緩解銅的供需矛盾，增強我國在國際銅市場上的競爭能力是十分必要的。高含泥氧化銅礦屬于復(fù)雜難處理銅資源的典型代表，這類礦石的特點是易破碎、含泥量大，氧化率高。我國氧化銅礦資源中，高含泥氧化銅礦占 2/3。傳統(tǒng)選冶工藝無法從中經(jīng)濟地回收銅[1?3]，采用堆浸技術(shù)是一種比較合理的工藝。但是，由于含泥量大導(dǎo)致的堆浸過程中滲透性差的問題，直接影響堆浸工藝的應(yīng)用與推廣。堆浸技術(shù)的成功與否主要取決于兩個方面[4]：一方面是溶浸液對有用礦物的溶解能力；另一方面是溶浸液在礦堆中的滲透效果。高含泥氧化礦氧化率極高，具有良好的可浸性。因此，改善高泥礦堆滲透性就成為處理高含泥氧化礦的關(guān)鍵。制粒技術(shù)[5?7]是處理高泥氧化礦的有效手段，但是制粒工藝需要制粒設(shè)備，且要使用粘結(jié)劑；目前制粒粘結(jié)劑多為堿性，而高泥氧化礦最優(yōu)的浸出方式是在酸性條件下，因此，制粒工藝也受到很大的限制。在礦堆中添加疏松物得到了嘗試，但受到試劑來源、價格等因素的影響，大規(guī)模應(yīng)用還為時尚早。在布液方式方面，加壓逆流浸出技術(shù)和加負壓進行強化浸出技術(shù)目前還沒有得到工業(yè)應(yīng)用[8]。本文作者以云南某高含泥氧化銅礦為研究對象，通過實驗研究探明顆粒結(jié)構(gòu)、粒徑及泥質(zhì)對堆體滲透性的影響規(guī)律，提出水洗?分級筑堆的堆浸新工藝，并在云南某銅礦得到成功應(yīng)用，對處理高含泥氧化銅礦具有重要的指導(dǎo)意義。

1 實驗

1.1 實驗礦樣

1.1.1 礦樣的化學(xué)性質(zhì)

礦樣取自于云南某銅礦，主要元素分析結(jié)果如表1所列，銅品位在1.013%。銅礦銅物相分析結(jié)果如表2所列。通過表2所列銅物相來看，氧化率為63.02%，結(jié)合率為21.13%。說明該礦具有良好的可浸性。

表1 銅礦多元素分析Table 1 Chemical composition of copper ore (mass fraction, % )

表2 銅礦銅物相分析Table 2 Phase analysis of copper ore

1.1.2 礦樣的粒級組成

將礦樣倒入標準分析篩中，石子孔徑分別為10.0、5.0、2.0、1.0、0.5、0.2和0.1 mm，置于篩析機上震篩10～15 min，得該氧化礦的級配曲線如圖1所示。

圖1 實驗礦樣級配曲線圖Fig.1 Particle size distribution of ore sample

1.2 實驗原理

室內(nèi)滲透實驗按原理分有兩大類[9]：變水頭法和常水頭法。前者適用于滲透性較小的介質(zhì)，后者適用于透水性較大的介質(zhì)。根據(jù)取樣粒徑在?25 mm，屬于透水性較大的介質(zhì)，固滲透性實驗采用變水頭法測量滲透系數(shù)，變水頭法滲透實驗原理如圖2所示。在滲透試驗過程中時刻t0時，水頭量管中水頭為h0，t1時刻的水頭量管中水頭為 h1，在時間間隔Δt=t1?t2內(nèi)，量管中水量變化為

式中：Q為流量，m3/s；a為截面積，m2；h0和h1為時間t0與t1時刻的液位高度，m。

由Darcy定律可知，流經(jīng)樣品的流體流量為

式中：A為樣品斷面積，m2；L為樣品長度，m；h為水頭高度，m；K為滲透系數(shù)，m/s。

圖2 變水頭法滲透系數(shù)測試原理圖Fig.2 Schematic diagram of falling head permeability test

由式(1)和 (2)可得

由于本試驗裝置為統(tǒng)一管徑，a=A，故式(3)可變?yōu)?/p>

式中：L0為滲透柱上下表面之間的長度，m；L1為起始液面至滲透柱上表面長度，m；L2為最終液面至滲透柱上表面長度，m；tΔ為液面從L1到L2經(jīng)歷的時間，s。

1.3 實驗步驟

1)將原礦進行篩分分級為+1～?5 mm、+1～?2 mm、?1 mm、?0.6 mm、?0.28 mm、+1～?25 mm、+2～?25 mm、+5～?25 mm；

2)按照上述實驗方法進行滲透性測試；

3)對堆場中的+1～?5 mm 和+5～?25 mm 兩級水洗礦樣進行滲透性測試；

4)記錄實驗數(shù)據(jù)，對比原礦與分級礦樣、水洗礦樣滲透性變化規(guī)律。

2 結(jié)果分析

不同粒級滲透性測試結(jié)果見表3。由表3可見，?0.28 mm、?0.6 mm和?1 mm粒級礦樣的滲透系數(shù)低于未分級原礦的滲透系數(shù)，說明對于粒級在?1 mm以下的礦樣，通過分級來改善顆粒結(jié)構(gòu)的措施對提高浸堆的滲透性效果不明顯。分析原因主要是?1 mm粒級的顆粒為泥質(zhì)，其對溶浸液有黏滯作用。所以，將?1 mm粒級的礦樣不入堆進行堆浸，而進行槽浸和攪拌浸出。+1～?2 mm和+1～?5 mm粒級礦樣的滲透系數(shù)分別為原礦的3.2和9.4倍，+1～?25 mm、+2～?25 mm和+5～?25 mm粒級礦樣的滲透性也大幅度改善，其滲透系數(shù)分別是原礦的19.1、75.9和47.6倍，說明調(diào)節(jié)粒徑組合可明顯提高浸堆的滲透性。

表3 不同粒徑滲透系數(shù)記錄表Table 3 Log sheet of osmotic coefficient of different particle size

對顆粒粒徑與滲透性進行擬合分析(見圖3)，結(jié)果表明，顆粒直徑與其組成的堆體的滲透系數(shù)成二階指數(shù)關(guān)系。當粒徑小于10 mm時，粒徑改變對提高浸堆的滲透性無明顯作用；當顆粒粒徑大于10 mm時，顆粒粒徑變化對浸堆的滲透性影響明顯。

圖3 顆粒直徑與滲透系數(shù)的關(guān)系Fig.3 Relationship between particle size and permeability coefficient

3 高泥礦堆滲透性影響因素分析

3.1 黏土礦物對滲透性的影響

高泥礦堆中的黏土礦物與水相互作用時，在黏土顆粒周圍會產(chǎn)生一個電場[10]，呈強極性狀態(tài)的水分子和水溶液中的陽離子在土粒表面的靜電引力作用下，被土粒強烈地吸附，它們共同構(gòu)成了強結(jié)合水層(見圖4(a))。強結(jié)合水因受土粒表面的強吸附力控制而不服從靜水力學(xué)規(guī)律，不受重力控制而流動。隨著離土粒表面距離的加大，水化陽離子的濃度逐漸降低，直至達到孔隙中水溶液的正常濃度為止，水分子亦從定向排列逐漸變成正常的隨機狀態(tài)。從強結(jié)合水層表面直至陽離子濃度正常這個范圍稱為弱結(jié)合水層(見圖4(b))，弱結(jié)合水的黏滯度亦較正常狀態(tài)的大。強結(jié)合水層和弱結(jié)合水層構(gòu)成了擴散層，土粒表面的負電荷和擴散層合起來稱為雙電層，隨著土顆粒距離進一步加大，顆粒之間孔隙變大，可以存在自由流動的自由水(見如圖4(c))，自由水受重力控制，決定了整個礦堆的滲透性。

圖4 不同狀態(tài)下黏土顆粒表面水理性質(zhì)Fig.4 Water physical properties of clay particle surface in different states: (a)Solid or semisolid state; (b)Plastic state; (c)Flow state

黏土顆粒的膨脹性是其水理性質(zhì)的重要體現(xiàn)，其膨脹大致分兩個階段[11]：第一階段為干顆粒表面吸附單層水分子形成粒間膨脹；第二階段由于雙電層的形成，使黏?；蚓舆M一步推開，形成滲透膨脹。黏土礦物遇水時膨脹，使得土粒間孔隙更小，阻礙自由水的流動；其次黏土中的結(jié)合水對自由水有較強黏滯作用而產(chǎn)生很大的阻力，這也正是高泥礦堆滲透性差的一個重要原因。

3.2 顆粒結(jié)構(gòu)對滲透性的影響

浸堆的顆粒結(jié)構(gòu)是指入堆散體顆粒和孔隙的空間相互排列以及顆粒間的聯(lián)結(jié)特征的綜合，其中的孔隙是溶浸液滲流的通道，而孔徑大小及其分布決定著溶浸液的滲流路徑和滲流速度。因此，堆浸礦巖散體的結(jié)構(gòu)是影響溶液滲流的主要因素，并最終影響到浸出效果。按顆粒的排列及聯(lián)結(jié)情況可以分為下列4種[12]。

1)單粒結(jié)構(gòu)

單粒結(jié)構(gòu)是破碎制粒堆浸散體的特征，其特征是顆粒粒徑大小相似，顆粒間沒有聯(lián)結(jié)力，或聯(lián)結(jié)力非常微弱，可以忽略不計。按顆粒的相互排列，可以把單粒結(jié)構(gòu)分為疏松的和緊密的(見圖5)。

2)蜂窩狀結(jié)構(gòu)

當浸堆的粗顆粒(粒徑大于 5 mm)體積分數(shù)大于70%時，礦石顆粒以松散拱架形式堆積，細料填不滿粗料間的空隙，堆體整體以大顆粒構(gòu)成骨架承受載荷，小顆?；静怀惺茌d荷，浸堆具有大孔隙性、彈性，這些性質(zhì)與散體的強度和變形有密切關(guān)系(見圖6)。

3)包裹結(jié)構(gòu)

浸堆的粗顆粒體積分數(shù)小于30%時，粗顆粒礦塊被粉礦分割包裹，浸堆呈現(xiàn)黏土特性，一般孔隙性和滲透性較低(見圖7)。當粗顆粒體積分數(shù)介于 30%和70%之間時，浸堆的力學(xué)特性和滲流特性由粗細顆粒共同承擔。

4)A型和B型結(jié)構(gòu)

ABERG[13?15]根據(jù)粒級曲線的形狀，將顆粒的結(jié)構(gòu)分為兩種：A型和B型。A型顆粒集合體中微細顆粒比較多，顆粒級配曲線的下端相對較陡；而B型顆粒集合體中微細顆粒相對較少，因此，顆粒級配曲線的下端比較平緩(見圖8)，這類顆粒中的微細顆粒容易發(fā)生遷移，而導(dǎo)致孔隙堵塞。在A型結(jié)構(gòu)的顆粒集合體中，每個顆粒和周圍充分多的顆粒接觸，其位置是不變的。具有B型結(jié)構(gòu)的顆粒集合體中，僅有部分顆粒位置固定，他們組成土的骨架，其它顆粒和周圍的顆粒接觸較少，可在不擾動周圍其它顆粒的情況下在一定范圍內(nèi)自由移動(見圖9)。

圖5 浸堆體系的單粒結(jié)構(gòu)[12]Fig.5 Single particle structures of dump leaching system[12]:(a)Loose; (b)Closely packed

圖6 浸堆體系的蜂窩狀結(jié)構(gòu)[12]Fig.6 Cellular structure of dump leaching system[12]

圖7 浸堆體系的包裹結(jié)構(gòu)[12]Fig.7 Surrounded structure of dump leaching system [12]

圖8 A和B型顆粒結(jié)構(gòu)級配曲線Fig.8 Particle size distribution curves of particle structure types A and B

圖9 B型顆粒結(jié)構(gòu)堵塞示意圖Fig.9 Schematic diagram of blockage of particle structure type B

由圖1所示級配曲線可以得出，該高含泥氧化銅礦為的粗顆粒(+5 mm)在 63%左右，介于包裹結(jié)構(gòu)與蜂窩狀結(jié)構(gòu)之間，滲流特性由粗細顆粒共同承擔；其次，曲線下凹，為B型顆粒結(jié)構(gòu)，細顆粒容易在粗顆粒組成的骨架顆粒之間運移，在孔隙較小的地方容易發(fā)生拱堵現(xiàn)象，降低浸堆整體的滲透性，這也是高泥礦堆滲透性隨著浸出時間的增加逐漸降低的主要原因。

3.3 分級筑堆對滲透性的影響機理

分級筑堆的設(shè)計原理就是利用不同粒級物料的滲透性不同來設(shè)計筑堆。盡管粗粒級區(qū)的孔隙大，在噴淋強度很小的情況下，從粗粒級區(qū)流出的溶液所占的比率遠遠小于從細粒級區(qū)流出溶液所占的比率。這主要是因為細粒級區(qū)比表面積大，對溶液的吸附能力強，并且細粒級區(qū)孔隙和裂隙直徑較小，能產(chǎn)生較大的負孔隙水壓力，使溶液從粗粒級區(qū)進入細粒級區(qū)[16?17]，如圖10所示。隨著噴淋強度的加大，細粒級區(qū)逐漸達到飽和狀態(tài)，溶液開始作橫向流動，進入粗粒級區(qū)，粗粒級區(qū)出液率逐漸增大，并遠遠超過細粒級區(qū)。因此，在噴淋強度較大的情況下，液流優(yōu)先發(fā)生在粗粒級區(qū)。

圖10 粗、細顆粒區(qū)內(nèi)優(yōu)先液流示意圖Fig.10 Schematic diagrams of preferential flow in coarse and fine particles

4 工程應(yīng)用

4.1 存在問題

2007年9月，云南某銅礦建成濕法廠相應(yīng)的基礎(chǔ)設(shè)施，進行了堆浸工業(yè)試驗，首期入堆礦量10多萬t，主要為露天開采的氧化礦，礦石經(jīng)兩段破碎后，粒級在?25 mm，為入堆高含泥原礦，由自卸汽車運輸筑堆，并開始布液浸出實驗。至2007年11月底，經(jīng)過近3個月的堆浸試生產(chǎn)，浸出效果差，浸出率不到10%，堆場滲透性差，溶浸液無法正常滲透循環(huán)，浸堆表面徑流與積液現(xiàn)象嚴重(見圖11)，電銅廠堆浸因此停產(chǎn)。

圖11 水洗?分級前浸堆表面溶浸液積聚Fig.11 Solution gathering in surface of heap before washingclassification

4.2 水洗?分級工藝

根據(jù)以上分析研究，本文作者提出了基于洗礦機+旋流器的水洗?分級工藝。水洗指采用洗礦機將礦石進行預(yù)處理，?1 mm細砂礦約占處理礦量的20%左右，對堆體的滲透性具有惡化作用，不入堆進行堆浸，而改為槽浸和攪拌浸出，其工藝流程如圖12所示。首先將原礦礦石破碎至?25 mm，并進行篩分。+5～?25 mm進入洗礦機，其返砂粗粒級直接入堆，而+1～?5 mm粒級直接進入旋流器繼續(xù)分級，返砂+1～?5 mm粒級入堆，? 1 mm粒級進入攪拌槽進行槽浸。?0.28 mm粒級的泥質(zhì)進入濃密機，50%濃度底流進行攪拌浸出，酸度為50 g/L，液固比為3:1，浸出2 h后進行壓濾，濾液進行澄清后進入料液池，濾餅保證水分不超過15%，加石灰中和后進入堆渣場。

分級指將礦石根據(jù)粒級組成分為粗中細粒級，并對粗、細顆粒進行分區(qū)堆置，從而達到提高粒級均勻度的目的(見圖13)。塊狀礦石量較大，可堆放在堆場中央?yún)^(qū)域。細粒礦石產(chǎn)量較小，在堆場周邊區(qū)域堆放，杜絕均勻混合堆放，保證堆場主體部分能夠順利浸出。在分區(qū)堆放之后，采用不同的布液強度進行分區(qū)布液。塊狀礦石區(qū)采用較大的布液強度，細粒礦石區(qū)采用較小的布液強度，使整個堆場得到充分的浸潤。+5 mm的塊狀礦石堆放區(qū)域采用的布液強度為 20～100 L/(m2·h)，+1～?5 mm的細粒礦石堆放區(qū)域采用的布液強度為 5～20 L/(m2·h)。

4.3 效果分析

圖12 水洗?分級浸出工藝流程示意圖Fig.12 Schematic diagram of dump leaching with washing and classification

圖13 分級筑堆、分區(qū)布液示意圖Fig.13 Schematic diagram of classification heap and partition liquid arrangement

通過采用水洗?分級工藝，浸堆的滲透性明顯改善，消除了堆場表面的徑流與積液現(xiàn)象。圖14所示為采用此工藝后的堆場表面，較未采用此工藝堆場表面(見圖11)積液明顯減少。在浸堆處取兩種粒級的礦樣(+1～?5 mm和+5～?25 mm)，對其滲透性進行測試并與原礦進行對比，其測試結(jié)果見表4。由表4可見，原礦的滲透系數(shù)最低，且隨著浸出時間的增長，下降也比較快。+1～?5 mm粒級礦樣的滲透系數(shù)為原礦的8～10倍，+5～?25 mm粒級礦樣的滲透系數(shù)為原礦的50倍左右，且隨著浸出時間的增長，滲透系數(shù)下降緩慢且能滿足堆浸要求，說明水洗?分級工藝對改善浸堆滲透性效果顯著。

表4 現(xiàn)場堆樣滲透系數(shù)變化記錄表Table 4 Permeability coefficient of different size fraction in heap leaching

圖14 水洗?分級后堆場浸堆表面Fig.14 Surface of dump after washing-classification

此浸出工藝特點是堆浸為主，槽浸、攪浸為輔的堆浸?槽浸?攪拌浸聯(lián)合浸出；進入堆浸的礦量占70%左右，浸出率在59.3%左右；+0.28～?1 mm粒級的細砂進行槽浸，浸出率71.75%；將?0.28 mm粒級的微細粒泥礦進行攪拌浸出，浸出率在 43.5%；綜合浸出率達到了63.98%。

5 結(jié)論

1)分析了高含泥礦堆滲透性差的原因。首先，黏土礦物遇水膨脹導(dǎo)致顆粒間孔隙變小，阻礙了自由水的順利通過；其次，顆粒表面結(jié)合水對溶浸液自由水具有黏滯及吸收作用，增大溶浸液的滲流阻力。

2)堆體的顆粒粒徑與其滲透系數(shù)成二階指數(shù)關(guān)系，當粒徑小于10 mm時，其滲透系數(shù)隨著粒徑的增加而緩慢增加；當粒徑大于10 mm時，其滲透系數(shù)隨著粒徑的增加迅速增加。

3)?1 mm粒級物料對堆體滲透性影響最大，對堆體滲透性具有惡化作用，不宜入堆。以?1 mm粒徑作為堆浸與槽浸或攪浸的分界線。

4)提出了水洗?分級浸出新工藝。礦堆中+1～?5 mm和+5～?25 mm粒級物料占70%，采用分級筑堆、分區(qū)噴淋的堆浸方式；+0.28～?1 mm粒級物料視為細砂，進入進行槽浸；?0.28 mm視為泥質(zhì)，進行攪拌浸出。運用此工藝可使堆體的滲透性明顯提高，滲透系數(shù)為原礦的8～50倍，浸出率由原來的不到10%提高到63.98%。

[1]趙 勁.高含泥、高風化金礦石合理選礦工藝研究[J].云南冶金, 2006, 25(3): 291?296.ZHAO Jin.A study on the rational dressing technology of Au ore with high argillation & weathering intensity[J].Yunnan Metallurgy, 2006, 25(3): 291?296.

[2]孔勝武.某高含泥氧化銅礦石的可選性探討[J].有色金屬(選礦部分), 2002(4): 1?3.KONG Sheng-wu.Discussion of the wash ability of some high clay oxidized copper[J].Nonferrous Metals(Mineral Processing Section), 2002(4): 1?3.

[3]胡福成, 張健良.低品位高含泥鐵帽金礦提金工藝的研究[J].礦冶, 1995, 4(2): 51?56.HU Fu-cheng; ZhANG Jian-liang.Study on gold extraction from low-grade iron cap gold ore with high mud content[J].Mining and Metallurgy, 1995, 4(2): 51?56.

[4]李青松, 吳愛祥, 尹升華.礦堆中溶浸液滲流效果的控制機理研究[J].西部探礦工程, 2005(12): 73?75.LI Qing-song, WU Ai-xiang, YIN Shen-hua.Study of control mechanical on effect of seepage in dump leaching[J].West-china Exploration Engineering, 2005(12): 73?75.

[5]譚海明.低品位高含泥氧化銅礦的制粒堆浸新工藝研究[J]邵陽學(xué)院學(xué)報: 自然科學(xué)版, 2005, 2(1): 92?93.TAN Hai-ming.A study on granulated heap leaching of lowly graded and highly-mudded oxidized copper ore[J].Journal of Shaoyang University: Natural Science, 2005, 2(1): 92?93

[6]呂 萍.低品位高含泥氧化銅礦制粒堆浸新工藝的研究[J].礦業(yè)研究與開發(fā), 2001, 21(2): 32?34.Lü Ping.A study on granulated heap leaching of low graded and highly-mudded oxidized copper ore[J].Mining Research and Development, 2001,21(2): 32?34.

[7]黎湘虹, 黎澄宇, 王 卉.鑫泰含泥氧化銅礦制粒預(yù)處理堆浸工藝[J]．有色金屬，2009(1): 86?90.LI Xiang-hong, LI Cheng-yu, WANG Hui.Acidified granulation pretreatment-heap leaching processing of argillious oxidized copper ore from Xintai Co.Ltd [J].Nonferrous Metals, 2009(1):86?90.

[8]王洪江, 吳愛祥, 顧曉春, 張新普, 張 儀.高含泥氧化銅礦石分粒級筑堆技術(shù)及其應(yīng)用[J].黃金, 2011, 32(2): 46?50.WANG Hong-jiang, WU Ai-xiang, GU Xiao-cun, ZHANG Xin-pu, ZHANG Yi. Technique of classification heap-constructing of copper oxide ore with high-clay and its application[J].Gold, 2011, 32(2): 46?50.

[9]王貽明, 吳愛祥, 左 恒, 楊保華.微粒滲濾沉積作用對銅礦排土場滲流特性的影響[J].中國有色金屬學(xué)報, 2007, 17(12):2074?2078.WANG Yi-ming, WU Ai-xiang, ZUO Heng, YANG Bao-hua.Effect of particles sedimentation during leaching on seepage characteristic of copper dumps[J].The Chinese Journal of Nonferrous Metals, 2007, 17(12): 2074?2078.

[10]洪毓康.土質(zhì)學(xué)與土力學(xué)[M].2版.北京: 人民交通出版社，1993: 27?29.HONG Yu-kang.Soil science and soil mechanics[M].2nd Edition.Beijing: China Communications Press, 1993: 27?29.

[11]吳愛祥, 李青松, 尹升華, 姜立春.改善高泥礦堆滲透性的機理研究[J].湘潭礦業(yè)學(xué)院學(xué)報, 2003, 18(4): 1?5.WU Ai-xiang, LI Qing-song, YIN Sheng-hua, JIANG Li-chun.Study on the mechanism of improving the penetrability of the high-clay ore heaps[J].Journal of Xiangtan Mine Institute, 2003,18(4): 1?5.

[12]王貽明.應(yīng)力波強化堆浸滲流的理論與試驗研究[D].長沙:中南大學(xué), 2008: 84?86.WANG Yi-ming.Theory and test study on seepage reinforcement by stress wave in heap leaching[D].Changsha:Central South University, 2008: 84 ?86.

[13]ABERG B.Void ratio of non-cohesive soils and similar materials[J].Journal of Geotechnical Engineering, ASCE, 1992,118(9):1315?1334.

[14]ABERG B.Hydraulic conductivity of Non-cohesive soils[J].Journal of Geotechnical Engineering, ASCE, 1992, 118(9):1335?1347.

[15]ABERG B.Washout of grains from filtered sand and gravel materials[J].Journal of Geotechnical Engineering, ASCE, 1993,119(1): 36?53.

[16]WU Ai-xiang, YIN Sheng-hua, LIU Jin-zhi, YANG Bao-hua.Formative mechanism of preferential solution flow[J].Journal of Central South University of Technology, 2006, 13(5): 590?594.

[17]LUXMOORE R J, JARDINE P M, WILSON G V, ZELAZNY L W.Physical and chemical controls of preferred path flow through hill slope[J].Geodenna, 1990, 46:139?154.