鮞狀赤鐵礦石深度還原過程中金屬鐵顆粒粒度預(yù)測模型

孫永升，韓躍新，高鵬，李艷軍

?

鮞狀赤鐵礦石深度還原過程中金屬鐵顆粒粒度預(yù)測模型

孫永升，韓躍新，高鵬，李艷軍

(東北大學(xué)資源與土木工程學(xué)院，遼寧沈陽，110819)

以湖北官店鮞狀赤鐵礦為研究對象，對其進(jìn)行深度還原試驗，利用光學(xué)顯微圖像分析技術(shù)對還原物料中金屬鐵顆粒粒度進(jìn)行測量，考察還原溫度和還原時間對鐵顆粒粒度的影響，并采用MATLAB軟件對試驗數(shù)據(jù)進(jìn)行擬合分析，建立鐵顆粒粒度與還原條件之間的數(shù)學(xué)模型。研究結(jié)果表明：不同還原條件下金屬鐵顆粒粒度累積特性曲線呈現(xiàn)出相同的變化規(guī)律；升高還原溫度或延長還原時間可使鐵顆粒粒度明顯增加；建立鐵顆粒粒度80與還原溫度和還原時間之間的預(yù)測模型；模型的計算值與試驗值具有良好的吻合性，可用于預(yù)估深度還原過程中金屬鐵顆粒的粒度；基于該模型，可通過調(diào)整溫度和時間以實現(xiàn)金屬鐵顆粒粒度的優(yōu)化與控制。

深度還原；鐵顆粒；粒度；擬合分析；數(shù)學(xué)模型

針對采用常規(guī)選礦方法難以實現(xiàn)鐵礦物有效富集的復(fù)雜難選鐵礦資源，研究人員提出了深度還原?磁選新技術(shù)，即在低于礦石熔化溫度下，以煤粉為還原劑將鐵礦物還原為金屬鐵，并促使金屬鐵聚集生長為一定粒度的鐵顆粒，還原物料經(jīng)磨礦磁選獲得可用于煉鋼的金屬鐵粉[1]。深度還原?磁選技術(shù)為復(fù)雜難選鐵礦資源的開發(fā)利用開辟了新途徑，成為選礦和冶金領(lǐng)域研究熱點之一。國內(nèi)外學(xué)者圍繞深度還原技術(shù)及理論進(jìn)行了大量試驗研究，取得了一些具有科學(xué)意義和應(yīng)用價值的研究成果。鮞狀赤鐵礦、高鋁鐵礦、鐵尾礦、白云鄂博鐵礦、赤泥等含鐵原料深度還原提鐵研究表明，在適宜的還原溫度、還原時間、還原劑用量條件下，可獲得金屬化率大于90%的還原物料，經(jīng)磁選可得到鐵品位85%~95%，鐵回收率大于90%的深度還原鐵粉[2?10]。深度還原包含鐵礦物還原和金屬鐵聚集生長2個主要過程[11]。深度還原過程中，鐵礦物按照Fe2O3→Fe3O4→FeO→Fe的順序逐級反應(yīng)，最終還原為金屬鐵，其還原歷程包含多個階段，各階段動力學(xué)模型不同[12?14]。由于復(fù)雜難選含鐵物料中Si和Al等雜質(zhì)含量較高，鐵礦物還原生成的低價鐵氧化物(FeO)易與SiO2和Al2O3反應(yīng)形成鐵橄欖石和尖晶石等復(fù)雜鐵化合物，隨著還原過程的進(jìn)行，復(fù)雜鐵化合物會被進(jìn)一步還原為金屬鐵[15?16]。金屬鐵顆粒的形成及生長特性是深度還原理論體系的重要組成部分，因為只有當(dāng)金屬鐵生長到一定粒度，才能實現(xiàn)金屬鐵相與渣相的良好解離，進(jìn)而實現(xiàn)金屬鐵的磁選富集。目前，僅少數(shù)文獻(xiàn)指出鐵顆粒粒度可以采用圖像分析技術(shù)進(jìn)行測量，提高還原溫度和延長還原時間有利于金屬鐵顆粒的聚集生長[17?20]。因此，針對復(fù)雜難選鐵礦深度還原過程中金屬顆粒粒度特性進(jìn)行系統(tǒng)研究十分必要。本文作者以我國湖北官店鮞狀赤鐵礦為研究對象，對其進(jìn)行深度還原試驗，采用圖像分析技術(shù)對還原樣品中金屬鐵顆粒粒度進(jìn)行檢測，得出不同試驗條件下金屬鐵顆粒粒度變化規(guī)律。在此基礎(chǔ)上，利用MATLAB軟件建立金屬鐵顆粒粒度80與還原溫度、時間之間的數(shù)學(xué)模型，并對模型進(jìn)行驗證，以期預(yù)測深度還原過程中鐵顆粒粒度的變化趨勢，從而實現(xiàn)鐵顆粒粒度的優(yōu)化與控制。

1 材料與方法

1.1 原料

所用礦樣為取自湖北官店的鮞狀赤鐵礦，主要化學(xué)成分見表1。由表1可知：礦石中TFe品位(質(zhì)量分?jǐn)?shù)，下同)為42.21%，SiO2為21.80%，Al2O3為5.47%，CaO為4.33%，有害元素S的含量較低，而P含量達(dá)1.31%。試驗采用的還原劑為產(chǎn)自吉林的普通煙煤，表2給出了其工業(yè)分析結(jié)果。由表2可以發(fā)現(xiàn)：煤的固定碳含量較高，達(dá)到67.83%(質(zhì)量分?jǐn)?shù)，下同)；揮發(fā)分和灰分含量相對較低，分別為18.45%和12.02%；有害元素S和P含量很低。礦樣和煤分別經(jīng)顎式破碎機(jī)粗碎、對輥破碎機(jī)細(xì)碎至粒度小于2 mm，然后混勻分樣用于深度還原試驗。

表1 礦樣的主要化學(xué)成分(質(zhì)量分?jǐn)?shù))

表2 煤的工業(yè)分析結(jié)果(質(zhì)量分?jǐn)?shù))

1.2 試驗方法

深度還原試驗在KSL?1400X箱式電阻爐中進(jìn)行。將礦樣和煤粉按質(zhì)量比100:29.2混合均勻。取80 g混勻樣品裝入體積為100 mL的陶瓷坩堝中，并在樣品頂部均勻鋪一層約2 mm厚的煤粉(2 g)用于保證坩堝內(nèi)的還原氣氛。電阻爐升溫至設(shè)定溫度后，迅速將裝有樣品的坩堝放入爐膛內(nèi)，并開始計時。當(dāng)還原時間達(dá)到試驗設(shè)計值時還原試驗完成，迅速將還原樣品取出，并通過水冷方式將其冷卻至室溫。冷卻樣品經(jīng)過濾后，在80 ℃下由真空干燥箱烘干，最終得到還原物料。

1.3 測試方法

深度還原物料中金屬鐵顆粒粒度可通過光學(xué)顯微圖像分析技術(shù)進(jìn)行有效測量[17, 21]。主要步驟為：1) 取適量的還原物料與膠木粉按照體積比1:5均勻混合，采用鑲嵌機(jī)將其壓制成直徑30 mm、高度10 mm的圓柱體，采用拋光機(jī)對圓柱體進(jìn)行拋光，使樣品中金屬相的表面暴露和平滑；2) 利用BX41M型金相顯微鏡(配備Moticam 2300型數(shù)字相機(jī))在拋光截面上隨機(jī)選取30個視域進(jìn)行圖像采集，獲得圖像存儲于計算機(jī)硬盤中；3) 采用Motic Image Advanced 3.2圖像分析軟件對圖像進(jìn)行手動分割、選定區(qū)域、邊緣強(qiáng)化、二值化處理，運用軟件中的測量模塊對目標(biāo)區(qū)域的個數(shù)和粒度進(jìn)行統(tǒng)計測量，獲得每個顆粒的粒度。需要指出的是：圖像分析技術(shù)測量的顆粒直徑為當(dāng)量表觀直徑，即與顆粒二維截面積相等的圓的直徑，故本文中測得的顆粒粒度為其當(dāng)量表觀粒度。

在礦物加工領(lǐng)域，顆粒粒度通常用顆粒質(zhì)量的累積產(chǎn)率分布加以描述。已有研究表明：金屬顆粒主要以球狀或類球狀顆粒形式存在深度還原物料中[9,17]。故此，本文將金屬鐵顆粒視為球體處理，則依據(jù)球體的體積公式可以得出鐵顆粒累積粒度特性曲線。

式中：()為粒度小于的負(fù)累計產(chǎn)率(質(zhì)量分?jǐn)?shù))，%；d為金屬鐵顆粒粒度，μm；為粒級邊界，μm；為粒度≤的鐵顆粒個數(shù)；為鐵顆?？倐€數(shù)。

2 金屬鐵顆粒粒度特性

不同試驗條件下還原物料中金屬鐵顆粒的粒度累積特性曲線如圖1所示。由圖1可知：不同還原條件下金屬鐵顆粒粒度累積曲線呈現(xiàn)出相同的變化趨勢，即隨著粒徑的增加，負(fù)累積產(chǎn)率明顯增加。同時，還可以發(fā)現(xiàn)還原溫度和還原時間對金屬顆粒粒度特性影響顯著。在同一還原溫度下，隨著還原時間的延長，鐵顆粒粒度累積曲線依次向右移動，表明鐵顆粒粒度隨還原時間的延長而逐漸增大。由圖1(c)可見：當(dāng)還原溫度為1 200 ℃，還原時間由10 min延長到120 min時，金屬鐵顆粒的特征粒徑50和80分別由148.77 μm和213.96 μm增加到531.22 μm和596.07 μm。在相同還原時間條件下，隨著還原溫度的升高，曲線也呈現(xiàn)出依次向右移動的變化規(guī)律。例如，當(dāng)還原時間為50 min，還原溫度從1 150℃升高到1 300 ℃時，金屬鐵顆粒的50和80分別由194.32 μm和257.47 μm增加到772.92 μm和946.95 μm。上述結(jié)果表明：隨著還原時間的延長和還原溫度的升高，還原物料中金屬顆粒粒度逐漸增大。這是由于還原時間越長，金屬鐵的含量就越高，金屬鐵顆粒生長越充分；還原溫度升高則使得金屬鐵的擴(kuò)散速率增加，進(jìn)而促進(jìn)了金屬鐵的聚集生長。

3 金屬鐵顆粒粒度模型

3.1 數(shù)據(jù)選擇及處理

還原溫度和還原時間是影響鐵顆粒生長的2個關(guān)鍵因素，通過改變還原溫度和還原時間，可以調(diào)控鐵顆粒的粒度。為實現(xiàn)深度還原過程中鐵顆粒粒度的調(diào)控，本文采用MATLAB數(shù)學(xué)軟件研究鐵顆粒粒度與還原溫度和還原時間之間的關(guān)系，建立鐵顆粒粒度預(yù)測數(shù)學(xué)模型。礦物加工領(lǐng)域通常用80來表示顆粒粒度的分布特征，故選取金屬鐵顆粒的80作為建模數(shù)據(jù)。由圖1可得出還原溫度為1 150~1 300 ℃，還原時間為10~120 min范圍內(nèi)的80，如表3所示。

由于在對某一變量進(jìn)行比較時，盡管具有相同的量綱，但其數(shù)值的變化幅度、變化趨勢存在著較大的差異，造成數(shù)值比較不便。為了消除這一影響，在建模之前需要采用標(biāo)準(zhǔn)差法對自變量的數(shù)值進(jìn)行標(biāo)準(zhǔn)化處理，其轉(zhuǎn)化公式為

式中：x為經(jīng)標(biāo)準(zhǔn)化處理后數(shù)據(jù)；為原始數(shù)據(jù)；為平均值；為標(biāo)準(zhǔn)差。經(jīng)過處理后，數(shù)據(jù)符合標(biāo)準(zhǔn)的正態(tài)分布，消除了原始數(shù)據(jù)因數(shù)量水平和變化幅度不同而造成的影響。還原溫度和還原時間標(biāo)準(zhǔn)化處理后的結(jié)果如表4所示。

表3 不同還原工藝條件下金屬鐵顆粒的D80

溫度/℃：(a) 1 150；(b) 1 175；(c) 1 200；(d) 1 225；(e) 1 250；(f) 1 275；(g) 1 300

表4 還原溫度和還原時間的標(biāo)準(zhǔn)化結(jié)果

注：平均值為1 225，標(biāo)準(zhǔn)差為50.40；平均值為56.67，標(biāo)準(zhǔn)差為35.24。

3.2 模型建立

以還原溫度和還原時間為自變量，鐵顆粒粒度80為因變量，采用MATLAB軟件對其進(jìn)行擬合分析，結(jié)果如圖2所示。圖2中擬合曲面的方程為

80(，)=481.3+221.8+188.1+63.52+72.3?10.12

(1 150≤≤1 300 ℃，10≤≤120 min) (3)

相關(guān)系數(shù)2為0.975 5，表明擬合效果良好。由圖2可以看出：試驗數(shù)據(jù)點基本全部分布在擬合曲面上，試驗值與模型值之間的差值較小，隨機(jī)散布在0?0平面兩側(cè)。同時還可發(fā)現(xiàn)：模型值隨還原溫度和還原時間的變化規(guī)律與試驗規(guī)律相同，均隨著還原溫度升高和還原時間延長而不斷增大。綜合上述分析，可以得出擬合曲面對原觀測值擬合良好，擬合方程可作為描述金屬鐵顆粒粒度80與還原溫度和還原時間之間關(guān)系的數(shù)學(xué)模型。

(a) 曲面投影到θ?t平面的等高線圖；(b) 模型主圖；(c) 殘差圖

3.3 模型驗證

由擬合方案計算出本文試驗條件下金屬鐵顆粒粒度的模型預(yù)測值，預(yù)測值與測量值之間的比較結(jié)果如圖3所示。從圖3可以發(fā)現(xiàn)：模型預(yù)測值數(shù)據(jù)點在直線=上下兩側(cè)均勻分布且呈現(xiàn)出良好的線性關(guān)系，說明預(yù)測值與試驗測量值之間具有良好的一致性。這一結(jié)果表明建立的數(shù)學(xué)模型精度較高，可以用于預(yù)測深度還原過程中金屬鐵顆粒的粒度。

為進(jìn)一步驗證該模型的準(zhǔn)確性與可行性，選取還原溫度1 220 ℃、還原時間45 min作為試驗點進(jìn)行驗證。經(jīng)測量該條件下的金屬鐵顆粒粒度累積特性曲線如圖4所示。

圖3 金屬鐵顆粒粒度模型預(yù)測值與試驗測量值對比

圖4 還原溫度1 220 ℃和還原時間45 min條件下金屬鐵顆粒粒度累積特性曲線

鐵顆粒的特征粒徑80為392.14 μm。將=1 220 ℃，=45 min標(biāo)準(zhǔn)化后代入模型80(，)= 481.3+221.8+188.1+63.52+72.3?10.12，可求得80(，)=398.89 μm。據(jù)此可知模型預(yù)測結(jié)果與實際測量結(jié)果較為接近，故本文建立的鐵顆粒粒度模型具有良好的穩(wěn)定性。結(jié)合該模型，通過調(diào)整和可實現(xiàn)對深度還原過程中鐵顆粒粒度的優(yōu)化與控制。

4 結(jié)論

1) 鮞狀赤鐵礦石深度還原過程中，金屬鐵顆粒粒度累積特性曲線在整體上呈現(xiàn)相同的變化規(guī)律，即隨著粒徑的增大，負(fù)累積產(chǎn)率逐漸增加。還原溫度和還原時間對鐵顆粒粒度影響顯著，隨著還原溫度升高和還原時間延長，鐵顆粒粒度累積曲線依次向右移動，說明鐵顆粒粒度不斷變大。

3) 提出的金屬鐵顆粒粒度模型可有效預(yù)測深度還原過程中金屬鐵顆粒粒度隨還原溫度和時間的變化趨勢，通過調(diào)節(jié)和從而實現(xiàn)對深度還原過程中鐵顆粒粒度的控制。然而該模型的建立是基于鮞狀赤鐵礦深度還原，對于其他難選鐵礦石深度還原過程中鐵顆粒粒度預(yù)測是否適用還需通過試驗進(jìn)一步研究。

[1] 韓躍新, 李艷軍, 劉杰, 等. 難選鐵礦石深度還原-高效分選技術(shù)[J]. 金屬礦山, 2011, 40(11): 1?4. HAN Yuexin, LI Yanjun, LIU Jie, et al. Deep reduction-efficient separation technology applied to refractory iron ore[J]. Metal Mine, 2011, 40(11): 1?4.

[2] LI Keqing, NI Wen, ZHU Ming, et al. Iron extraction from oolitic iron ore by a deep reduction process[J]. Journal of Iron and Steel Research, International, 2011, 18(8): 9?13.

[3] SUN Yongsheng, HAN Yuexin, GAO Peng, et al. Recovery of iron from high phosphorus oolitic iron ore using coal-based reduction followed by magnetic separation[J]. International Journal of Minerals, Metallurgy and Materials, 2013, 20(5): 411?419.

[4] CHUN Tiejun, LONG Hongming, LI Jiaxin. Alumina-Iron separation of high alumina iron ore by carbothermic reduction and magnetic separation[J]. Separation Science and Technology, 2015, 50(5): 760?766.

[5] SRIVASTAVA U, KAWATRA S K. Strategies for processing low-grade iron ore minerals[J]. Mineral Processing & Extractive Metallurgy Review, 2009, 30(4): 361?371.

[6] 賈巖, 倪文, 王中杰, 等. 拜耳法赤泥深度還原提鐵實驗[J]. 北京科技大學(xué)學(xué)報, 2011, 33(9): 1059?1064. JIA Yan, NI Wen, WANG Zhongjie, et al. Deep reduction experiments of Bayer red mud for iron recovery[J]. Journal of University of Science and Technology Beijing, 2011, 33(9): 1059?1064.

[7] 林海, 許曉芳, 董穎博, 等. 深度還原?弱磁選回收稀土尾礦中鐵的試驗研究[J]. 東北大學(xué)學(xué)報(自然科學(xué)版), 2013, 34(7): 1039?1044. LIN Hai, XU Xiaofang, DONG Yingbo, et al. Test research on recovery of iron from rare earth tailings by depth reduction roasting-low intensity magnetic separation[J]. Journal of Northeastern University (Natural Science), 2013, 34(7): 1039?1044.

以上測定均取3次獨立實驗結(jié)果的平均值，數(shù)據(jù)用平均值 ± 標(biāo)準(zhǔn)差(SD)表示。使用Origin 9.0軟件對數(shù)據(jù)進(jìn)行雙總體t檢驗。差異顯著性以P < 0.05表示。

[8] YU Wen, SUN Tichang, KOU Jue, et al. The function of Ca(OH)2and Na2CO3as additive on the reduction of high-phosphorus oolitic hematite-coal mixed pellets[J]. ISIJ International, 2013, 53(3): 427?433.

[9] LI Yongli, SUN Tichang, KOU Jue, et al. Study on phosphorus removal of high-phosphorus oolitic hematite by coal-based direct reduction and magnetic separation[J]. Mineral Processing & Extractive Metallurgy Review, 2014, 35(1): 66?73.

[10] KOU Jue, SUN Tichang, TAO D, et al. Coal-based direct reduction and magnetic separation of lump hematite ore[J]. Minerals & Metallurgical Processing, 2014, 31(3): 150?161.

[11] 韓躍新, 張成文, 孫永升, 等. Na2CO3促進(jìn)復(fù)雜難選鐵礦石深度還原的機(jī)理分析[J]. 東北大學(xué)學(xué)報(自然科學(xué)版), 2012, 33(11): 1633?1636. HAN Yuexin, ZHANG Chengwen, SUN Yongsheng, et al. Mechanism analysis on deep reduction of complex refractory iron ore promoted by Na2CO3[J]. Journal of Northeastern University (Natural Science), 2012, 33(11): 1633?1636.

[12] SUN Yongsheng, GAO Peng, HAN Yuexin, et al. Reaction behavior of iron minerals and metallic iron particles growth in coal-based reduction of an oolitic iron ore[J]. Industrial & Engineering Chemistry Research, 2013, 52(6): 2323?2329.

[13] LIU Guisu, STREZOV V, LUCAS J A, et al. Thermal investigations of direct iron ore reduction with coal[J]. Thermochemica Acta, 2004, 410(1/2): 133?140.

[14] SUN Yongsheng, HAN Yuexin, GAO Peng, et al. Investigation of kinetics of coal based reduction of oolitic iron ore[J]. Ironmaking and Steelmaking, 2014, 41(10): 763?768.

[15] 魏玉霞, 孫體昌, 寇玨, 等. 內(nèi)配煤用量對某難選鐵礦石壓塊直接還原焙燒的影響[J]. 中南大學(xué)學(xué)報(自然科學(xué)版), 2013, 44(4): 1305?1311. WEI Yuxia, SUN Tichang, KOU Jue, et al. Effect of coal dosage on direct reduction roasting of refractory iron ore briquettes[J]. Journal of Central South University (Science and Technology), 2013, 44(4): 1305?1311.

[16] 范敦城, 倪文, 李瑾, 等. 鐵尾礦再選粗精礦深度還原含鐵硅酸鹽礦物的生成與還原[J]. 中南大學(xué)學(xué)報(自然科學(xué)版), 2015, 46(6): 1973?1981. FAN Duncheng, NI Wen, LI Jin, et al. Generation and reduction mechanism of silicate minerals containing iron in deep reduction of rough concentrate from iron tailings[J]. Journal of Central South University (Science and Technology), 2015, 46(6): 1973?1981.

[17] GAO Peng, SUN Yongsheng, REN Duozheng, et al. Growth of metallic iron particles during coal-based reduction of a rare earths-bearing iron ore[J]. Minerals & Metallurgical Processing, 2013, 30(1): 75?80.

[18] SUN Yongsheng, HAN Yuexin, GAO Peng, et al. Growth kinetics of metallic iron phase in coal-based reduction of oolitic iron ore[J]. ISIJ International, 2016, 56(10): 1697?1704.

[19] 高鵬, 韓躍新, 李艷軍, 等. 基于數(shù)字圖像處理的鐵礦石深度還原評價方法[J]. 東北大學(xué)學(xué)報(自然科學(xué)版), 2012, 33(1): 133?136. GAO Peng, HAN Yuexin, LI Yanjun, et al. Evaluation on deep reduction of iron ore based on digital image processing techniques[J]. Journal of Northeastern University (Natural Science), 2012, 33(1): 133?136.

[20] 孫體昌, 秦曉萌, 胡學(xué)平, 等. 低品位鐵礦石直接還原過程鐵顆粒生長和解離特性[J]. 北京科技大學(xué)學(xué)報, 2011, 33(9): 1048?1052. SUN Tichang, QIN Xiaomeng, HU Xueping, et al. Grain growth and cleavage characteristics of metallic phase in direct reduction of a low-grade iron ore[J]. Journal of University of Science and Technology Beijing, 2011, 33(9): 1048?1052.

[21] SUN Yongsheng, HAN Yuexin, GAO Peng, et al. Particle Size Measurement of Metallic Iron in Reduced Materials Based on Optical Image Analysis[J]. Chemical Engineering & Technology, 2014, 37(12): 2030?2036.

(編輯 陳愛華)

Prediction model of metallic iron particle size during coal-based reduction of oolitic hematite ore

SUN Yongsheng, HAN Yuexin, GAO Peng, LI Yanjun

(College of Resources and Civil Engineering, Northeastern University, Shenyang 110819, China)

An oolitic iron ore taken from Guandian in Hubei province was reduced, and the size of metallic iron particles in reduced product was measured using optical image analysis. The effects of reduction temperature and time on the size of metallic iron particles were investigated. The experimental data were analyzed by MATLAB software, and the mathematic model of metallic iron particle size considering reduction conditions was proposed. The results indicate that the curves of size cumulative passing percentage of metallic iron particles present similar variation trend under different conditions. When the reduction temperature increases and reduction time extends, the size of metallic iron particles increases obviously. The prediction model for80of metallic iron particles considering reduction temperature and time was established. The calculation values determined by this model correlate well with the test results, indicating that it can be used to predict the particle size of metallic iron in coal-based reduction. Based on the model, the particle size of metallic iron can be optimized by means of adjusting reduction temperature and time.

coal-based reduction; metallic iron particles; particle size; fitting analysis; mathematic model

10.11817/j.issn.1672?7207.2017.02.002

TD925.7

A

1672?7207(2017)02?0282?07

2016?04?25；

2016?07?20

國家自然科學(xué)基金資助項目(51134002)；中央高?；究蒲袠I(yè)務(wù)費專項資金資助項目(N140108001)(Project(51134002) supported by the National Natural Science Foundation of China; Project(N140108001) supported by the Fundamental Research Funds for the Central Universities of China)

孫永升，博士，副教授，從事復(fù)雜難選礦產(chǎn)資源高效開發(fā)與利用研究；E-mail：neusunyongsheng@163.com