應(yīng)力狀態(tài)對煤巷頂板錨固孔鉆進(jìn)速度的影響

劉少偉,劉棟梁,馮友良,李鑫濤,尚鵬翔

(河南理工大學(xué)能源科學(xué)與工程學(xué)院,河南焦作 454003)

應(yīng)力狀態(tài)對煤巷頂板錨固孔鉆進(jìn)速度的影響

劉少偉,劉棟梁,馮友良,李鑫濤,尚鵬翔

(河南理工大學(xué)能源科學(xué)與工程學(xué)院,河南焦作 454003)

采用有限元數(shù)值軟件對煤巷頂板常見巖石在不同應(yīng)力狀態(tài)下模擬鉆進(jìn),發(fā)現(xiàn)應(yīng)力狀態(tài)對煤巷頂板巖石錨固孔的鉆進(jìn)速度影響不明顯。從鉆進(jìn)位移和實時鉆進(jìn)速度兩方面對煤巷頂板不同類型巖石鉆進(jìn)結(jié)果進(jìn)行分析討論,認(rèn)為平均鉆進(jìn)速度可作為衡量巖石類型的一個較好指標(biāo)。通過各類巖石瞬時鉆進(jìn)速度的分析,鉆頭的回彈是不同巖石力學(xué)特性的反應(yīng),可通過對回彈次數(shù)、正面積和負(fù)面積等參數(shù)對不同巖層進(jìn)行識別。

煤巷頂板錨固孔;鉆進(jìn)速度;應(yīng)力狀態(tài)

自20世紀(jì)90年代以來,錨桿支護(hù)以其顯著的技術(shù)和經(jīng)濟優(yōu)越性在煤巷圍巖控制中獲得了廣泛應(yīng)用,是巷道支護(hù)技術(shù)的一場革命[1－2]。煤巷頂板錨固孔鉆進(jìn)是錨桿支護(hù)的一個重要環(huán)節(jié),錨固孔鉆進(jìn)速度不僅影響著錨固孔的成孔的速度和效率,還對煤巷的施工管理和工程進(jìn)度有著重要的影響。煤巷頂板為沉積層狀巖體,各層的巖性不同,多為砂巖、粉砂巖、砂質(zhì)泥巖、泥巖、砂質(zhì)頁巖及石灰?guī)r等,其中直接頂主要以粉砂巖、砂頁巖、泥巖為主,基本頂以砂巖、石灰?guī)r及砂礫巖為主[3－4]。頂板錨固孔現(xiàn)場鉆進(jìn)過程中,常出現(xiàn)時快時慢現(xiàn)象。鉆進(jìn)速度的變化影響因素較多,包括頂板巖層物理力學(xué)結(jié)構(gòu)(頂板組合巖層的尺寸、節(jié)理發(fā)育程度、巖石類型)及應(yīng)力狀態(tài)等。反之,鉆進(jìn)速度也能反映頂板的巖層結(jié)構(gòu)及應(yīng)力狀態(tài)。頂板巖層結(jié)構(gòu)及應(yīng)力狀態(tài)是錨桿支護(hù)煤巷冒頂?shù)脑?所以通過頂板錨固孔的鉆進(jìn)速度可間接反映煤巷頂板的穩(wěn)定性,如圖1所示。本文主要研究應(yīng)力狀態(tài)與頂板錨固孔鉆進(jìn)速度的關(guān)系,為頂板錨固孔鉆進(jìn)速度的重要影響因素的確定提供理論參考。

圖1 煤巷頂板錨固孔鉆進(jìn)速度與冒頂?shù)年P(guān)系Fig.1 Relationship between the drilling velocity of anchorhole in coal roadway roof and roof fall

1 鉆進(jìn)模型基本假設(shè)及建立

1.1 鉆進(jìn)模型基本假設(shè)

煤巷頂板錨固孔的鉆進(jìn)是一個復(fù)雜的過程,主要受頂板巖體力學(xué)參數(shù)和錨桿鉆機參數(shù)的影響。為便于分析,模型做以下基本假設(shè)[3,5－7]:

(1)錨固孔鉆進(jìn)前,巖體未發(fā)生塑性變形,處于彈性變化范圍內(nèi),即:巷道隨掘隨支。

(2)鉆進(jìn)過程中,鉆孔軌跡控制良好,鉆頭以垂直于巷道頂板的方式進(jìn)行鉆進(jìn),且不考慮鉆桿與孔壁的摩擦影響。

(3)由于錨桿鉆機的鉆頭相對于巷道頂板巖體具有很高的硬度和強度,因此假設(shè)鉆頭為剛體。

(4)巖石單元鉆進(jìn)失效后,不再考慮其被重復(fù)破碎的問題,破碎的巖石單元不再影響后續(xù)的鉆進(jìn)工作。

(5)現(xiàn)場錨桿鉆機常采用氣動式錨桿鉆機,風(fēng)壓是氣動錨桿鉆機的動力源,給鉆機提供軸向推力和旋轉(zhuǎn)切削力。利用CAE軟件建立?28 mm的氣動錨桿鉆機鉆頭的三維模型,如圖2所示。鉆頭的邊界條件:推力為6 000 N,轉(zhuǎn)速220 r/min,給鉆頭施加除軸向的其他方向的約束,即保持鉆孔垂直于巷道頂板。

圖2 錨桿鉆機鉆頭三維模型Fig.2 3D model of bore bit

(6)研究表明,當(dāng)巷道埋深H≥20R0(R0為巷道半徑)時,可以忽略巷道影響范圍(3～5倍R0)內(nèi)巖石的自重(圖3),水平原巖應(yīng)力可以簡化為均布的[3]。因此,為便于對結(jié)果的處理,對組成煤巷頂板的各巖層進(jìn)行單獨鉆進(jìn)分析,不同應(yīng)力狀態(tài)下對不同巖層的鉆進(jìn)時間統(tǒng)一為2 s。

圖3 深埋巷道力學(xué)特點Fig.3 Mechanical characters of deep roadway

1.2 鉆進(jìn)模型的建立

鉆進(jìn)模型選用沽源煤礦煤巷頂板巖石,通過實驗室實驗得出煤巷頂板各巖層力學(xué)參數(shù)(表1)。

表1 巖石力學(xué)參數(shù)Table 1 Rock mechanics parameters

鉆進(jìn)的巖體模型尺寸為300 mm×300 mm× 300 mm,滿足鉆孔中心距邊界大于5倍的鉆頭直徑要求[8]。巖石和鉆頭模型采用四面體自由網(wǎng)格劃分,為了保證結(jié)果精度及節(jié)省計算時間,對鉆頭鉆進(jìn)影響較大的區(qū)域盡量細(xì)化網(wǎng)格,而其他地方適當(dāng)增加網(wǎng)格尺寸[6,9－10]。巖石的屈服準(zhǔn)則較多,對于均質(zhì)的巖石類型,常采用的是Drucker－Prager準(zhǔn)則和Mohr－Coulomb準(zhǔn)則[6]。Drucker－Prager準(zhǔn)則是在Coulomb準(zhǔn)則和Mises準(zhǔn)則的基礎(chǔ)上擴展和推廣而得,且在一定情況下可與Mohr－Coulomb準(zhǔn)則相互轉(zhuǎn)化[9],能較好地反映巖石的屈服過程,故選用Drucker－Prager模型,巖石的破壞選用剪切損傷破壞。圖4為煤巷頂板巖石鉆進(jìn)試驗?zāi)Ｐ褪疽?利用該方案對沽源煤礦某煤巷頂板的各類型巖石進(jìn)行鉆進(jìn),參數(shù)見表1[11],研究應(yīng)力狀態(tài)對煤巷頂板錨固孔鉆進(jìn)速度的影響。

圖4 煤巷頂板沉積巖體鉆進(jìn)試驗方案示意Fig.4 Drilling test plan schematic of sedimentary rock in coal roadway roof

2 鉆進(jìn)方案的確定

有關(guān)研究表明[2,12],一般將側(cè)壓系數(shù)大于1.2的區(qū)域為高構(gòu)造應(yīng)力區(qū),小于1.2的區(qū)域為一般構(gòu)造應(yīng)力區(qū)。因此,鉆進(jìn)方案主要分兩類各3組。

沽源煤礦的巷道埋深約800 m,垂直應(yīng)力σv約為20 MPa,最小水平主應(yīng)力約為0.75倍的垂直應(yīng)力,即15 MPa,通過最大水平主應(yīng)力變化得出不同應(yīng)力狀態(tài)下的鉆進(jìn)方案(表2)。

表2 不同應(yīng)力狀態(tài)下鉆進(jìn)方案Table 2 Drilling programs under different stress states

3 模擬結(jié)果分析

3.1 錨固孔鉆進(jìn)特征分析

錨固孔的鉆進(jìn)是在氣動錨桿鉆機的作用下,通過風(fēng)壓來控制氣腿的推力和鉆頭轉(zhuǎn)速完成旋轉(zhuǎn)切削破巖。圖5為各類巖石在高構(gòu)造應(yīng)力區(qū)(λ=1.5)情況下鉆進(jìn)2 s時的鉆進(jìn)深度。泥巖、砂質(zhì)泥巖、粉砂巖、灰?guī)r的鉆進(jìn)深度分別為83.8,59.0,41.1,34.9 mm,不同巖石的鉆進(jìn)深度不同。不同巖石類型具有不同的平均鉆進(jìn)速度,不同平均鉆進(jìn)速度也正是不同類型巖石的反映。由此可知,平均鉆進(jìn)速度是評價巖石類型的較好指標(biāo),對于煤巷頂板錨固孔鉆進(jìn)過程的巖層類型識別以及預(yù)防煤巷冒頂有重要意義。

為具體分析錨固孔鉆進(jìn)過程,對鉆進(jìn)砂質(zhì)泥巖(鉆頭完全進(jìn)入)的瞬時速度進(jìn)行分析,如圖6所示。由圖可知,鉆頭在錨固孔的鉆進(jìn)整體呈“上臺階式”鉆進(jìn)過程,鉆頭存在回彈的現(xiàn)象,巖石鉆進(jìn)速度不均勻,呈波動性變化。

圖6 砂質(zhì)泥巖的鉆進(jìn)過程示意Fig.6 Illustration of the process of drilling sandy mudstone

3.2 應(yīng)力狀態(tài)對錨固孔鉆進(jìn)速度的影響

根據(jù)模擬結(jié)果得出不同應(yīng)力狀態(tài)下各巖石鉆進(jìn)時鉆頭位移,見表3。

表3 不同應(yīng)力狀態(tài)下各巖石鉆進(jìn)時鉆頭位移Table 3 Drilling bit displacement various types of rock under different stress statesmm

由表3可看出:①相同應(yīng)力狀態(tài)下,鉆進(jìn)各類型巖石的鉆頭位移不同,即其平均鉆進(jìn)速度也是不同的,主要與巖石的力學(xué)參數(shù)有一定的關(guān)系。②鉆進(jìn)各類型巖石鉆頭的位移隨最大水平主應(yīng)力σH(或側(cè)壓系數(shù)λ)的增加并沒有呈現(xiàn)出明顯的變化規(guī)律。③在不同應(yīng)力狀態(tài)下,鉆進(jìn)泥巖時鉆頭位移最大為105.0 mm,最小為102.0 mm,相差為3 mm;砂質(zhì)泥巖時最大位移為80.3 mm,最小為75.6 mm,相差為4.7 mm;粉砂巖時位移最大為62.7 mm,最小為59.8 mm,相差僅1.9 mm;灰?guī)r時位移最大為57.0 mm,最小為53.3 mm,相差為3.7 mm。通過鉆進(jìn)過程數(shù)據(jù)的實時采集(剔除初始鉆頭與巖石的距離20 mm,如圖6(a)中的鉆頭未接觸巖石區(qū)域),可以得出各巖石類型的平均鉆進(jìn)速度:泥巖的平均鉆進(jìn)速度最大為46.4 mm/s,最小為44.7 mm/s,相差僅1.7 mm/s;砂質(zhì)泥巖的平均鉆進(jìn)速度最大為30.3 mm/s,最小為32.9 mm/s,相差為2.6 mm/s;粉砂巖的平均鉆進(jìn)速度最大為23.3 mm/s,最小為21.7 mm/s,相差僅1.6 mm/s;灰?guī)r的平均鉆進(jìn)速度最大為20.5 mm/s,最小為18.0 mm/s,相差為2.5 mm/s。僅從鉆進(jìn)速度來說,應(yīng)力狀態(tài)對各巖石的平均鉆進(jìn)速度均影響不大,但是一定長度錨固孔的鉆進(jìn)是一個時間問題,不同的巖性之間更是有差別。以砂質(zhì)泥巖為例,若鉆進(jìn)2 400 mm的錨固孔,最短為74.3 s,最長為80.4 s,相差僅6.1 s,而鉆進(jìn)位移卻相差達(dá)180 mm,可見應(yīng)力狀態(tài)雖然對砂質(zhì)泥巖鉆進(jìn)速度的影響不明顯,但如果錨固孔長度較大時,例如長度較大的錨索孔(5～10 m)鉆進(jìn)時間存在較大影響。因此,不論高構(gòu)造應(yīng)力區(qū)域還是一般構(gòu)造應(yīng)力區(qū)域的煤巷對頂板錨固孔的鉆進(jìn)時間存在一定的影響,但對鉆進(jìn)速度影響不明顯。

3.3 錨固孔鉆進(jìn)速度的探討

由于應(yīng)力狀態(tài)對錨固孔鉆進(jìn)速度的影響不大,取一般構(gòu)造應(yīng)力區(qū)(λ=1.1)情況下鉆進(jìn)不同巖石類型時鉆頭的位移和瞬時鉆進(jìn)速度進(jìn)行分析(圖7,8)。

圖7 鉆進(jìn)不同巖石類型時鉆頭位移Fig.7 Drilling bit displacement of different types of rock

(1)從鉆進(jìn)位移看巖石的平均鉆進(jìn)速度。

圖8 鉆進(jìn)不同巖石類型時鉆頭瞬時鉆進(jìn)速度Fig.8 Instantaneous rock drilling velocity of different types of rock

分析圖7可知,不同巖石的位移鉆進(jìn)曲線在一定斜率的直線附近波動,即鉆進(jìn)速度可以通過直線的斜率近似的表示,稱為平均鉆進(jìn)速度。4種巖石的平均鉆進(jìn)速度大致有以下關(guān)系:泥巖＞砂質(zhì)泥巖＞粉砂巖＞灰?guī)r,灰?guī)r與粉砂巖的平均鉆進(jìn)速度差別不大,而灰?guī)r與泥巖的差別還是比較明顯的。由于煤巷頂板是不同巖性的巖層組合,其厚度不一、巖石力學(xué)參數(shù)也各不同,在實際煤巷頂板錨固孔鉆進(jìn)過程中,鉆進(jìn)的位移曲線是波動性的前進(jìn)過程,即使實時監(jiān)測,當(dāng)巖層發(fā)生變化時,易誤認(rèn)為仍在同一巖層中鉆進(jìn),尤其是對于巖石的力學(xué)參數(shù)比較接近的巖層。當(dāng)巖石鉆進(jìn)速度變化不太明顯的情況下,僅從鉆進(jìn)位移看鉆進(jìn)速度具有明顯的滯后性,因此,相關(guān)學(xué)者[13－16]對頂板錨固孔的鉆進(jìn)的研究,采用多因素對頂板巖層結(jié)構(gòu)進(jìn)行識別,如:鉆進(jìn)速度－轉(zhuǎn)速－推力。

(2)從實時鉆進(jìn)速度看平均鉆進(jìn)速度。

由圖8可知,4種巖石的鉆進(jìn)速度均呈波動狀,波形各不相同。圖中各巖石的回彈次數(shù)(即“0”下的振動次數(shù))不盡相同。鉆進(jìn)過程中鉆頭的回彈是正?，F(xiàn)象,根據(jù)回彈原理可測定巖石的硬度[17],巖石的硬度對錨固孔的鉆進(jìn)速度有一定的影響,不同的鉆進(jìn)速度在一定程度上可以反映出不同強度的巖石,可見回彈是不同巖石特性的反應(yīng),巖石在鉆進(jìn)過程中的回彈也正是不同巖石力學(xué)參數(shù)的反應(yīng)?；貜棿螖?shù)越多,回彈的時間就越長,鉆進(jìn)的位移就越小,平均鉆進(jìn)速度也就越小。對于速度曲線的速度線與時間軸所圍成的面積,即為路程,在圖中面積可分為正面積SZ和負(fù)面積SF,正面積為向上鉆進(jìn)路程,負(fù)面積向下回彈的路程,那么,平均鉆進(jìn)速度=(SZ－SF)/t。回彈次數(shù)越少,正面積越大,負(fù)面積越小,巖石的平均鉆進(jìn)速度就越大。圖中4種巖石間的鉆進(jìn)回彈次數(shù)之間的關(guān)系為:泥巖＜砂質(zhì)泥巖＜粉砂巖,泥巖＜砂質(zhì)泥巖＜灰?guī)r,粉砂巖(15次)和灰?guī)r(13次)間的回彈次數(shù)相近;且負(fù)面積間也有明顯的關(guān)系:泥巖＜砂質(zhì)泥巖＜粉砂巖＜灰?guī)r;由表3可知,4種巖石的平均鉆進(jìn)速度的關(guān)系為:泥巖＞砂質(zhì)泥巖＞粉砂巖＞灰?guī)r,而粉砂巖和灰?guī)r間的平均鉆進(jìn)速度較為接近。事實上,鉆頭回彈次數(shù)的監(jiān)測是方便易行的,實時鉆進(jìn)過程中采集鉆進(jìn)速度曲線也是具有可操作性的,故在鉆進(jìn)中,對回彈次數(shù)、正面積和負(fù)面積進(jìn)行實時分析,在巖層識別方面有著重要的意義。

4 結(jié) 論

(1)應(yīng)力狀態(tài)對煤巷頂板巖石錨固孔鉆進(jìn)速度的影響不明顯。因此,煤巷頂板錨固孔鉆進(jìn)過程中速度變化時,出現(xiàn)的偶發(fā)冒頂現(xiàn)象,應(yīng)力狀態(tài)對其影響不大,為實驗室鉆進(jìn)模擬模型簡化提供了理論基礎(chǔ)。

(2)從鉆進(jìn)位移和實時鉆進(jìn)速度兩個方面對煤巷頂板不同類型巖石鉆進(jìn)結(jié)果分析可知,平均鉆進(jìn)速度是衡量巖石類型的較好指標(biāo),為現(xiàn)場快速判別煤巷頂板錨固孔內(nèi)巖石類型提供依據(jù)。

(3)鉆頭的回彈是不同巖石力學(xué)特性的反應(yīng),利用回彈次數(shù)、正面積和負(fù)面積等參數(shù)可進(jìn)行巖石類型的判斷,為煤巷頂板錨固孔內(nèi)巖層特征識別和冒頂?shù)念A(yù)防具有重要的意義。

參考文獻(xiàn):

[1] 康紅普.煤巷錨桿支護(hù)理論與成套技術(shù)[M].北京:煤炭工業(yè)出版社,2007.

[2] 劉少偉.錨桿支護(hù)煤巷冒頂危險的應(yīng)力影響及工程應(yīng)用[J].采礦與安全工程學(xué)報,2007,24(2):239－242.

Liu Shaowei.Effect of stress on roof falling in coal drift supported by bolts and its engineering application[J].Journal of Mining&Safety Engineering,2007,24(2):239－242.

[3] 錢鳴高,石平五.礦山壓力與巖層控制[M].徐州:中國礦業(yè)大學(xué)出版社,2003.

[4] 車書成,張榮偉.煤礦地質(zhì)學(xué)[M].徐州:中國礦業(yè)大學(xué)出版社, 1996.

[5] 陳 勇,萬教育,趙子仁,等.PDC鉆頭破巖的動態(tài)數(shù)值模擬[J].新疆石油天然氣,2009,5(3):69－73.Chen Yong,Wan Jiaoyu,Zhao Ziren,et al.Dynamic numerical simulation of PDC bits rock breaking[J].Xinjiang Oil&Gas,2009,5 (3):69－73.

[6] 林鐵軍,練章華,孟英峰,等.空氣鉆井中動態(tài)破巖有限元仿真研究[J].巖石力學(xué)與工程學(xué)報,2008,27(S2):3592－3597.

Lin Tiejun,Lian Zhanghua,Meng Yingfeng,et al.Finite elements method study on dynamic rock breaking in air drilling[J].Chinese Journal of Rock Mechanics and Engineering,2008,27 (S2):3592－3597.

[7] 馬鳳清,唐 波,馮光通,等.三牙輪鉆頭與地層相互作用動力學(xué)數(shù)值模擬[J].石油礦場機械,2012,41(9):5－9.

Ma Fengqing,Tang Bo,Feng Guangtong,et al.Study of dynamical numerical simulation of the interaction between three-cone bit and formation[J].Oil Field Equipment,2012,41(9):5－9.

[8] 徐小荷,余 靜.巖石破碎學(xué)[M].北京:煤炭工業(yè)出版社,1984.

[9] 朱以文,蔡元奇,徐 晗.ABAQUS與巖土工程分析[M].香港:中國圖書出版社,2005.

[10] 陳衛(wèi)忠,伍國軍,賈善坡.ABAQUS在隧道及地下工程中的應(yīng)用[M].北京:中國水利水電出版社,2010.

[11] 河南理工大學(xué)巖石力學(xué)實驗室.沽源金牛能源有限公司巷道圍巖物理力學(xué)參數(shù)試驗報告[R].焦作:河南理工大學(xué),2010.

Laboratory of Rock Mechanics,Henan Polytechnic University.Test report of rock mechanical properties of Guyuan Jinniu Energy Co., Ltd[R].Jiaozuo:Henan Polytechnic University,2010.

[12] 蒲文龍,張國華.長溝峪煤礦地應(yīng)力測量及煤巷穩(wěn)定性分析與控制[J].煤炭工程,2006(10):68－70.

Pu Wenlong,Zhang Guohua.The ore crustal stress measure of Changgouyu Mine and analyse the stability of country rock and controlled[J].Coal Engineering,2006(10):68－70.

[13] 劉少偉,馬念杰,武建棟,等.錨桿支護(hù)煤巷冒頂危險區(qū)預(yù)測及工程應(yīng)用[M].北京:煤炭工業(yè)出版社,2009.

[14] 賈明魁.巖層組合劣化探測技術(shù)研究[J].中國礦業(yè)大學(xué)學(xué)報, 2006,35(1):44－48.

Jia Mingkui.Research on detecting technique for inferior strata combination[J].Journal of China University of Minging&Technology,2006,35(1):44－48.

[15] 張汝海,劉少偉.煤巷錨桿支護(hù)全線跟蹤設(shè)計方法[J].中國礦業(yè),2007,16(2):54－56.

Zhang Ruhai,Liu Shaowei.Design method of bolt spport tracking on whole line in roadway[J].China Mining Magazine,2007,16(2): 54－56.

[16] 鄧廣濤,賈明魁,馬念杰.煤巷頂板巖層結(jié)構(gòu)的探測與確定[J].煤炭工程,2006(2):48－50.

Deng Guangtao,Jia Mingkui,Ma Nianjie.Detection and determination of roof rock structure in mine roadway[J].Coal Engineering, 2006(2):48－50.

[17] 賴海輝.機械巖石破碎學(xué)[M].長沙:中南工業(yè)大學(xué)出版社, 1991.

Influence of stress states on drilling velocity of anchorage hole on coal roadway roof

LIU Shao-wei,LIU Dong-liang,FENG You-liang,LI Xin-tao,SHANG Peng-xiang
(School of Energy Science and Engineering,Henan Plytechnic University,Jiaozuo 454003,China)

The paper used finite element software to simulate drilling of common rocks in coal roadway roof under different stress states,and got different stress states on roadway roof rock have no obvious impact on anchor hole drilling velocity.From two aspects including drilling displacement and real time velocity,drilling results of various rocks in coal roadway roof were analyzed and discussed,meanwhile the conclusion that average drilling velocity can be considered as a good index to measure the type of rock is got.Through analyzing instantaneous drilling velocity of various rocks,it is obtained that the resilience of bore bit may react different rock mechanics characteristics,various rocks can be identified by different parameters such as resilience frequency,positive area,negative area and so on.

anchorage hole on coal roadway roof;drilling velocity;stress states

TD353

A

0253－9993(2014)04－0608－06

劉少偉,劉棟梁,馮友良,等.應(yīng)力狀態(tài)對煤巷頂板錨固孔鉆進(jìn)速度的影響[J].煤炭學(xué)報,2014,39(4):608－613.

10.13225/j.cnki.jccs.2013.1087

Liu Shaowei,Liu Dongliang,Feng Youliang,et al.Influence of stress states on drilling velocity of anchorage hole on coal roadway roof[J].Journal of China Coal Society,2014,39(4):608－613.doi:10.13225/j.cnki.jccs.2013.1087

2013－07－28 責(zé)任編輯:常 琛

國家自然科學(xué)基金面上基金資助項目(51274087);國家自然科學(xué)基金資助項目(51104055);河南省青年骨干教師資助項目(2010GGJS－053)

劉少偉(1977—),男,遼寧錦州人,副教授,碩士生導(dǎo)師,博士。Tel:0391－3987948,E－mail:lswxll@126.com