趙麗娟,田 震,王 野

(1.遼寧工程技術(shù)大學(xué)機(jī)械工程學(xué)院,遼寧阜新 123000;2.沈陽(yáng)航天新光集團(tuán)有限公司,遼寧沈陽(yáng) 110043)

采煤機(jī)外噴霧系統(tǒng)數(shù)值模擬研究

趙麗娟1,田 震1,王 野2

(1.遼寧工程技術(shù)大學(xué)機(jī)械工程學(xué)院,遼寧阜新 123000;2.沈陽(yáng)航天新光集團(tuán)有限公司,遼寧沈陽(yáng) 110043)

為提高采煤機(jī)外噴霧系統(tǒng)性能,改善工作面作業(yè)環(huán)境,基于Pro/E與GAMBIT聯(lián)合建立了采煤工作面模型,選用標(biāo)準(zhǔn)k-ε湍流模型構(gòu)造采煤機(jī)外流場(chǎng)數(shù)學(xué)模型,并通過FLUENT仿真得到巷道內(nèi)風(fēng)流平均速度為1.8 m/s左右,符合《煤礦安全規(guī)程》的規(guī)定;對(duì)采煤機(jī)外噴霧系統(tǒng)進(jìn)行了數(shù)值模擬,通過改變噴嘴口徑、噴霧壓力和擴(kuò)散角等參數(shù),得到了不同噴霧參數(shù)下的霧滴粒徑和霧化濃度的變化規(guī)律:隨著噴霧壓力的增大,霧滴粒徑減小,噴霧濃度分布更加均勻;小口徑噴嘴霧化濃度較高,大口徑噴嘴霧化區(qū)域濃度均勻、捕集面積大;小霧化擴(kuò)散角下的霧粒速度相對(duì)較大,大擴(kuò)散角下的霧流濃度分布均勻,但其穿透力差。分析結(jié)果為提高采煤機(jī)噴霧降塵效率,降低采煤工作面粉塵濃度提供了依據(jù)。

采煤機(jī);外噴霧;粒徑;霧化濃度;數(shù)值模擬

近年來(lái),高濃度煤礦粉塵引起的事故造成了大量的人員傷亡,如2005年11月,黑龍江東風(fēng)煤礦發(fā)生特大煤塵爆炸事故,造成171人死亡,48人受傷[1]; 2006年7月,山西藺家莊煤礦發(fā)生一起特別重大煤塵爆炸事故,共有53人遇難[2];2007年12月,山西新窯煤礦由于煤塵爆炸繼而引起瓦斯爆炸,造成105人遇難[3];2012年5月新疆洪山煤礦粉塵爆炸事故造成6名礦工死亡[4]。高濃度粉塵不僅使煤礦的安全開采面臨嚴(yán)峻的考驗(yàn),還將會(huì)導(dǎo)致煤礦工人患上職業(yè)塵肺病,危害煤礦工人的身心健康[5-6]。據(jù)相關(guān)數(shù)據(jù)統(tǒng)計(jì),我國(guó)塵肺病患者2010年已累計(jì)超過67萬(wàn)人,2011年,我國(guó)煤礦檢出塵肺病人達(dá)14 000例,呈上升趨勢(shì)[7]。因此,降低煤礦井下工作面粉塵濃度和減少呼吸性粉塵的危害已成為煤礦安全生產(chǎn)的主要研究課題。

在綜采工作面,采煤機(jī)截割過程產(chǎn)生的粉塵占整個(gè)工作面粉塵總量的60%～85%[8]。在采煤機(jī)開采過程中,工作面粉塵量過大、噴嘴易堵等問題使采煤機(jī)作業(yè)環(huán)境粉塵的沉降難以達(dá)到煤礦安全規(guī)程的要求[9]。本文基于流體動(dòng)力學(xué)仿真技術(shù),對(duì)采煤機(jī)外噴霧系統(tǒng)進(jìn)行數(shù)值模擬,得到不同噴嘴口徑、噴霧壓力和擴(kuò)散角條件下噴霧系統(tǒng)的工作效果,找出了不同噴霧參數(shù)下的霧滴粒徑和霧化濃度的變化規(guī)律;將數(shù)值模擬技術(shù)應(yīng)用到煤礦粉塵控制的研究中,通過調(diào)整噴霧參數(shù),達(dá)到提高噴霧降塵效果的目的,為改善工作面作業(yè)環(huán)境和提高煤礦生產(chǎn)安全提供了參考。

1 模型的建立

綜采工作面較多的設(shè)備(如采煤機(jī)、刮板輸送機(jī)、液壓支架等)導(dǎo)致巷道內(nèi)部通風(fēng)區(qū)域形狀極為復(fù)雜,建立準(zhǔn)確模型十分困難,需對(duì)計(jì)算模型進(jìn)行適當(dāng)?shù)暮?jiǎn)化[10-11],如圖1所示,其中工作面尺寸如下:進(jìn)風(fēng)區(qū)8.5 m(長(zhǎng))×4.17 m(寬)×1.3 m(高);采空區(qū)17.5 m(長(zhǎng))×4.8 m(寬)×1.3 m(高);工作面高度為1.3 m,液壓支柱間距1.5 m。

圖1 綜采工作面Fig.1 Simplified diagram of fully mechanized coal face

利用Pro/E建立簡(jiǎn)化后的采煤機(jī)模型(其中在搖臂殼體輸出端頸部上下各分布3個(gè)外噴霧噴嘴),將其導(dǎo)入到GAMBIT中,并在此基礎(chǔ)之上建立綜采工作面的整體模型。在GAMBIT中劃分計(jì)算網(wǎng)格,考慮到數(shù)值模擬準(zhǔn)確度及仿真時(shí)間問題,網(wǎng)格精度應(yīng)根據(jù)物理量場(chǎng)在模型區(qū)域中的變化情況酌情而定。模型選擇Tet/Hybrid網(wǎng)格單元進(jìn)行劃分,劃分類型為TGrid[12],劃分好的網(wǎng)格模型如圖2所示。

圖2 模型網(wǎng)格劃分Fig.2 Mesh division ofmodel

為了便于對(duì)噴霧顆粒進(jìn)行采樣分析,研究噴霧參數(shù)變化對(duì)噴霧效果的影響,在滾筒與煤壁距離100～600 mm范圍內(nèi),沿滾筒軸線方向依次建立與煤壁平行的A～F六個(gè)采樣平面,如圖3所示。

圖3 采樣截面位置Fig.3 Sampling section position

2 外流場(chǎng)數(shù)值模擬

針對(duì)采煤機(jī)外流場(chǎng)的數(shù)值模擬,采用標(biāo)準(zhǔn)k-ε湍流模型,在湍動(dòng)能k方程的基礎(chǔ)上,引入湍動(dòng)耗散率ε方程建立封閉的數(shù)學(xué)模型,其中湍流動(dòng)能方程[13]為

式中,Gk為由平均速度梯度引起的湍動(dòng)能的產(chǎn)生項(xiàng); Gb為由浮力引起的湍動(dòng)能的產(chǎn)生項(xiàng);YM為可壓縮湍流中脈動(dòng)擴(kuò)張的貢獻(xiàn);C1ε,C2ε和C3ε為經(jīng)驗(yàn)常數(shù);σk和σε分別為與湍動(dòng)能和耗散率對(duì)應(yīng)的Prandtl數(shù);Sk和Sε為用戶定義的源項(xiàng)。

模型選擇Pressure-Swirl Atomizer(壓力旋轉(zhuǎn)霧化噴嘴)作離散相噴射源,用液滴的最可幾分布直徑以及分布指數(shù)來(lái)確定液滴的Rosin-Rammler分布[14], Rosin-Ramm ler分布是以累積體積分布來(lái)表達(dá)的,可用式(3)描述。

式中,Vc為粒徑在D以下的所有顆粒的體積與總體積的比值;D為霧滴粒徑;n為均勻度指數(shù);c為常數(shù)。

為了與實(shí)驗(yàn)數(shù)據(jù)更好的擬合,Rizk和Lefebvre等[15]對(duì)Rosin-Rammler分布進(jìn)行了修正,修正后的表達(dá)式為

求解得到整個(gè)工作面通風(fēng)流場(chǎng)的速度矢量,如圖4所示。風(fēng)流進(jìn)入工作面后,受到液壓支柱的阻礙作用,在其附近呈現(xiàn)繞流流動(dòng),導(dǎo)致風(fēng)流速度減小,模型中距離風(fēng)流入口最遠(yuǎn)處的液壓支柱附近風(fēng)流速度減小到了0.2 m/s左右。在采煤工作面內(nèi),采煤機(jī)阻礙了風(fēng)流運(yùn)動(dòng),并且占據(jù)了一部分截割段巷道空間內(nèi)的體積,促使風(fēng)流速度有所增加,在采煤機(jī)機(jī)身附近風(fēng)流速度大約為2.64 m/s。由于風(fēng)流在截割段沖擊煤壁和采煤機(jī)滾筒,致使?jié)L筒附近產(chǎn)生了回流區(qū)和渦流區(qū),整個(gè)巷道內(nèi)風(fēng)流平均速度為1.8 m/s左右,符合《煤礦安全規(guī)程》[16]中規(guī)定的0.25～4.0 m/s之內(nèi)。

圖4 通風(fēng)流場(chǎng)的速度矢量Fig.4 Velocity vectorgraph of ventilation flow field

3 噴霧參數(shù)對(duì)霧滴粒徑的影響分析

3.1 噴霧壓力對(duì)霧滴粒徑的影響

當(dāng)噴嘴口徑為1.5 mm、噴霧擴(kuò)散角為60°時(shí),對(duì)其分別施加3,6和8 MPa的壓力。利用采樣截面對(duì)噴霧顆粒進(jìn)行采樣統(tǒng)計(jì),得到霧滴平均粒徑Y(jié)向分布情況,如圖5所示。由圖5可見,壓力越大,噴嘴噴出的霧滴粒徑越小,不同壓力下霧滴粒徑沿程分布變化趨勢(shì)相同,且壓力越高,霧滴粒徑沿程分布越均勻,霧滴粒徑隨著距離煤壁越來(lái)越近而逐漸變大,低壓下霧滴粒徑沿程增大趨勢(shì)要大于高壓噴霧下的情況。

選定F截面作為采樣截面,對(duì)截面上霧滴粒徑大小進(jìn)行統(tǒng)計(jì),統(tǒng)計(jì)結(jié)果如圖6所示。當(dāng)壓力為3,6和8 MPa時(shí),在直徑60～160μm范圍內(nèi)的霧粒分別占采樣數(shù)據(jù)總數(shù)的78%,82%和74%。從粒徑分布的集中區(qū)域及其所占的百分比可以看出,隨著噴霧壓力的增加,霧粒直徑分布集中區(qū)域有逐漸向小粒度范圍集中的趨勢(shì)。

圖5 不同壓力下霧滴粒徑沿Y向分布Fig.5 Diameter distribution along the Y direction of different pressures

圖6 F截面上霧粒直徑的統(tǒng)計(jì)分布Fig.6 Diameter’s statistical distribution of spray particles at F section

3.2 噴嘴口徑對(duì)霧滴粒徑的影響

對(duì)口徑分別為1.0,1.2和1.5 mm的噴嘴(噴霧擴(kuò)散角均為60°)分別施加3,4,5和6 MPa的噴霧壓力,統(tǒng)計(jì)得到B截面在不同噴嘴口徑下的霧滴粒徑見表1。其中大口徑噴嘴的霧滴粒徑較大,但在相同壓力下,隨著噴嘴口徑的增大,霧滴粒徑的增幅逐漸減小;隨著噴霧壓力的提高,霧滴粒徑隨口徑增大而增大的幅度也逐漸減緩。

表1 不同噴嘴口徑下的霧滴粒徑Table 1 Diameter of spray par ticles of d ifferent nozzle caliber

3.3 霧化擴(kuò)散角對(duì)霧滴粒徑的影響

對(duì)口徑為1.0 mm、噴霧擴(kuò)散角分別為30°,60°和75°的噴嘴均施加3 MPa的壓力。選取B采樣平面進(jìn)行采樣統(tǒng)計(jì),得到在不同噴霧擴(kuò)散角度下的霧滴平均粒徑以及噴霧顆粒沿Y軸的平均速度見表2。霧滴粒徑隨著噴霧擴(kuò)散角度的增加而減小,噴霧擴(kuò)散角的增大,加大了噴嘴霧化的擴(kuò)散區(qū)域,同時(shí)擴(kuò)大了多個(gè)噴嘴之間的噴霧重疊區(qū),使更多的噴霧顆粒發(fā)生碰撞產(chǎn)生碎裂,進(jìn)而導(dǎo)致了噴霧粒徑的減小。隨著噴霧擴(kuò)散角的增大,在整體坐標(biāo)系下霧滴沿Y軸的運(yùn)動(dòng)速度減小,但是減小的幅度很小。由于Y軸方向并不是噴霧軸線的方向,因此,在B采樣平面處,隨霧化角的增大,霧滴沿Y軸方向的運(yùn)動(dòng)速度變化不是很明顯。

表2 不同噴霧擴(kuò)散角下的霧滴粒徑及速度Table2 Diameter and velocity of spray particles under different diffusion angles

4 噴霧參數(shù)對(duì)霧化濃度的影響分析

4.1 噴霧壓力對(duì)霧化濃度的影響

對(duì)口徑為1.0 mm、噴霧擴(kuò)散角60°的噴嘴,分別施加3和6 MPa的壓力,得到噴嘴霧化濃度如圖7所示。由圖7可見,采煤機(jī)搖臂外噴霧可以較好的將粉塵抑制在產(chǎn)塵源處,當(dāng)壓力為3 MPa時(shí),霧化區(qū)域周邊濃度較低,只有4.56 g/m3,但區(qū)域內(nèi)的最大濃度達(dá)到了91.2 g/m3,整體濃度主要集中在27.3～63.8 g/m3;當(dāng)壓力達(dá)到6 MPa時(shí),噴霧區(qū)域的霧化濃度更加的均勻,其濃度分布規(guī)律與3 MPa時(shí)大致相同,都是中間部分濃度高,周邊濃度低,但隨著壓力的增高,噴霧區(qū)域的整體濃度也隨之增加,整體濃度的大小主要集中在36～72 g/m3。

4.2 噴嘴口徑對(duì)霧化濃度的影響

對(duì)口徑分別為1.0和1.2 mm的噴嘴(噴霧擴(kuò)散角均為60°)施加3 MPa的壓力。選取F截面作為觀察平面,得到不同噴嘴口徑下的霧化濃度,如圖8所示。由圖8可見,當(dāng)壓力及霧化擴(kuò)散角相同時(shí),口徑為1.0 mm的噴嘴其霧化濃度高于口徑為1.2 mm的噴嘴。噴嘴口徑為1.0 mm時(shí),霧流較為緊密,噴霧區(qū)的中心處出現(xiàn)了高濃度區(qū),且噴霧區(qū)的整體濃度變化梯度比較明顯;當(dāng)噴嘴口徑為1.2 mm時(shí),霧流較為分散,噴霧區(qū)的濃度分布相對(duì)均勻,沒有出現(xiàn)明顯的梯度變化。這主要是由于大口徑噴嘴的霧化液滴具有較大粒徑,在噴霧過程中沒有進(jìn)一步破碎,進(jìn)而降低了噴霧區(qū)的濃度。

4.3 噴霧擴(kuò)散角對(duì)霧化濃度的影響

圖7 不同壓力下噴嘴霧化濃度Fig.7 Atomizing concentration of different pressures

圖8 不同噴嘴口徑在F截面處的霧化濃度Fig.8 Atomizing concentration of different calibers at F section

對(duì)口徑為1.0 mm、噴霧擴(kuò)散角分別為30°和75°的噴嘴施加3 MPa的壓力,得到噴嘴霧化濃度如圖9所示。由圖9可見,隨著噴霧擴(kuò)散角度的增大,噴霧區(qū)域體積增加,霧化濃度逐漸均勻。噴霧擴(kuò)散角為30°時(shí),霧流較細(xì),在噴嘴出口較遠(yuǎn)處各噴嘴的霧流開始重疊,搖臂上側(cè)與下側(cè)霧化區(qū)域中心處濃度分布較高,最大濃度達(dá)到98.2 g/m3;當(dāng)霧化角度為75°時(shí),霧流在離開噴嘴不遠(yuǎn)處開始重疊,霧滴發(fā)生激烈碰撞進(jìn)而破碎,使霧化區(qū)域濃度分布較為均勻,且霧化區(qū)域體積明顯增大,同時(shí)由于霧滴破碎成了較小液滴,在霧化區(qū)域周邊有部分液滴顆粒產(chǎn)生了比較明顯的飄散現(xiàn)象。

5 結(jié) 論

圖9 不同擴(kuò)散角下的霧化濃度Fig.9 Atom izing concentration of different diffusion angles

(1)由于受到液壓支柱和采煤機(jī)等設(shè)備的阻礙,巷道內(nèi)的風(fēng)流呈現(xiàn)出不均勻分布,其中在出風(fēng)口液壓支柱處風(fēng)速減小到了0.2 m/s左右,采煤機(jī)處風(fēng)速為2.64 m/s左右,在滾筒附近產(chǎn)生了回流區(qū)和渦流區(qū);整個(gè)巷道內(nèi)風(fēng)流平均速度在1.8 m/s左右,符合《煤礦安全規(guī)程》的規(guī)定。

(2)隨著噴霧壓力增加,霧粒分布集中區(qū)域出現(xiàn)逐漸向小粒度范圍集中的趨勢(shì),同時(shí)噴霧區(qū)域也會(huì)增大,且分布更加均勻。小口徑噴嘴因霧滴粒徑較小而使霧化濃度較高,但大口徑噴嘴所形成的霧化區(qū)域濃度更加均勻,捕集面積更大,粒徑增大的幅度隨口徑的增加逐漸趨于平緩。小擴(kuò)散角下的霧流較細(xì),霧粒速度相對(duì)較大;大擴(kuò)散角下的液膜破碎更加完全,使得其捕集范圍較大,霧流濃度分布均勻,但霧滴速度較小。在實(shí)際的生產(chǎn)中,根據(jù)不同噴霧參數(shù)下霧滴粒徑和霧化濃度的變化規(guī)律,對(duì)相關(guān)參數(shù)進(jìn)行合理匹配,可以提高噴霧降塵效率。

(3)基于本文提出的研究方法,可為采煤機(jī)外噴霧系統(tǒng)的優(yōu)化設(shè)計(jì)和性能評(píng)價(jià)提供技術(shù)支持,同時(shí)也為提高外噴霧降塵效率、降低采煤工作面粉塵濃度、減少呼吸性粉塵的危害提供參考。

[1] 梁 冬,曹霽陽(yáng),劉 錚.七臺(tái)河礦難系重大責(zé)任事故[N].新華每日電訊,2005-12-02(1).

[2] 呂曉宇.山西靈石藺家莊煤礦發(fā)生礦難[N].人民日?qǐng)?bào),2006-07-17(2).

[3] 車艷輝,王 軍.爆炸事故中受傷人員的心理應(yīng)激狀態(tài)分析及心理干預(yù)[J].中國(guó)傷殘醫(yī)學(xué),2011,19(7):99-100.

Che Yanhui,Wang Jun.Psychological stressanalysisand psychological intervention of injured in the explosion accident[J].Chinese Journal of Trauma and Disability Medicine,2011,19(7):99-100.

[4] 新疆煤礦爆炸事故續(xù):6人死亡因被困時(shí)間過長(zhǎng)[EB/N].http://news.sina.com.cn,2012-05-16.

[5] 聶百勝,李祥春,楊 濤,等.工作面采煤期間PM2.5粉塵的分布規(guī)律[J].煤炭學(xué)報(bào),2013,38(1):33-37.

Nie Baisheng,Li Xiangchun,Yang Tao,et al.Distribution of PM2.5 dust duringmining operation in coal workface[J].Journal of China Coal Society,2013,38(1):33-37.

[6] 時(shí)訓(xùn)先,蔣仲安,周姝嫣,等.綜采工作面粉塵分布規(guī)律的實(shí)驗(yàn)研究[J].煤炭學(xué)報(bào),2008,33(10):1117-1121.

Shi Xunxian,Jiang Zhongan,Zhou Shuyan,et al.Experimental study on dust distribution regularity of fully mechanized mining face[J].Journal of China Coal Society,2008,33(10):1117-1121.

[7] 武曉娟.塵肺病危害生命猛于事故[N].中國(guó)能源報(bào),2013-01-21(16).

[8] 安愛英.采煤機(jī)噴霧系統(tǒng)噴嘴噴霧仿真分析研究[D].太原:太原理工大學(xué),2010:1-4.

An Aiying.The nozzle simulation and analysis for shearer spary system[D].Taiyuan:Taiyuan University of Technology,2010:1-4.

[9] 馬素平,寇子明.噴霧降塵效率及噴霧參數(shù)匹配研究[J].中國(guó)安全科學(xué)學(xué)報(bào),2006,16(5):84-88.

Ma Suping,Kou Ziming.Study on efficiency of dust suppression by mist spray and its matched parameters[J].China Safety Science Journal,2006,16(5):84-88.

[10] Mohammad Masroor M lam.An integrated approach to dust control in coalmining face areas of a continuousminer and its computational fluid dynamics modeling[D].America:Southern Illinois University Carbondale,2006.

[11] 陳清華,張國(guó)樞.采煤機(jī)噴霧嘴內(nèi)流場(chǎng)數(shù)值模擬[J].煤炭科學(xué)技術(shù),2006,34(5):25-26,29.

Chen Qinghua,Zhang Guoshu.Numerical simulation on flow field in spray nozzle of coal shearer[J].Coal Science and Technology, 2006,34(5):25-26,29.

[12] 郝志勇,劉 偉.礦用局部通風(fēng)機(jī)改進(jìn)前后空氣動(dòng)力性能對(duì)比試驗(yàn)[J].煤炭學(xué)報(bào),2012,37(5):874-877.

Hao Zhiyong,Liu Wei.Aerodynamic performance contrast test before and after improvement of local fans for mines[J].Journal of China Coal Society,2012,37(5):874-877.

[13] 王曉珍,蔣仲安,王善文,等.煤巷掘進(jìn)過程中粉塵濃度分布的數(shù)值模擬[J].煤炭學(xué)報(bào),2007,32(4):386-390.

Wang Xiaozhen,Jiang Zhongan,Wang Shanwen,et al.Numerical simulation of distribution regularities of dust concentration during the ventilation process of coal roadway driving[J].Journal of China Coal Society,2007,32(4):386-390.

[14] Rosin P,Rammler E.The laws governing the fineness of powdered coal[J].Inst.Fuel,1933(7):29-36.

[15] Rizk N K,Lefebvre A H.Drop-size distribution characteristics of spill-return atomizers[J].AIA-A Propulsion and Power,1985, 1(3):16-22.

[16] 國(guó)家安全生產(chǎn)監(jiān)督管理總局,國(guó)家煤礦安全監(jiān)察局.煤礦安全規(guī)程[M].北京:煤炭工業(yè)出版社,2011.

Numerical simulation of shearer external spray system

ZHAO Li-juan1,TIAN Zhen1,WANG Ye2
(1.College ofMechanical Engineering,Liaoning Technical University,Fuxin 123000,China;2.Shenyang Aerospace Xinguang Group Co.,Ltd.,Shenyang 110043,China)

In order to improve the shearer external spray system performance and environment ofworking face,amodel ofmechanized coal face was established by the union of Pro/E and GAMBIT,a standard k-εturbulence flow model was selected to build the external flow field of shearer,the average speed of ventilation in roadway was about1.8 m/ s,itwas comply with the provisions of the Safety Standard of Coal Mine.The numerical simulation of shearer external spray situation wasestablished by the software FLUENT,the law of the influence about the effectof parameter variation to spray was obtained through the changes of nozzle caliber,spray pressure and diffusion angle.As the spray pressure increase,the diameter decreases and the spray concentration distribution ismore uniform;Small-bore nozzle’s concentration is higher,large-bore nozzle’s concentration is uniform and trapping area is large;the fog particles velocity with small diffusion angle is opposite bigger,the spray concentration distribution with large diffusion angle ismore uniform, but its penetration is poor.The result provides foundation for improving the efficiency of shearer spraying dust and reducing the coalface dust concentration.

shearer;external spray;diameter;atomizing concentration;numerical simulation

TD421.6

A

0253-9993(2014)06-1172-05

趙麗娟,田 震,王 野.采煤機(jī)外噴霧系統(tǒng)數(shù)值模擬研究[J].煤炭學(xué)報(bào),2014,39(6):1172-1176.

10.13225/j.cnki.jccs.2013.0981

Zhao Lijuan,Tian Zhen,Wang Ye.Numerical simulation of shearer external spray system[J].Journal of China Coal Society,2014,39(6): 1172-1176.doi:10.13225/j.cnki.jccs.2013.0981

2013-07-10 責(zé)任編輯:許書閣

中國(guó)煤炭工業(yè)科技計(jì)劃資助項(xiàng)目(MTKJ2009-264)

趙麗娟(1964—),女,遼寧阜新人,教授,博士生導(dǎo)師。Tel:0418-3352413,E-mail:zlj2120@163.com。通訊作者:田 震(1987—),男,安徽界首人,博士研究生。E-mail:lntutian2008@126.com