煤礦井下爆破作業(yè)頻帶能量分布規(guī)律*

肖　鑄,張鵬翔,3,李青松,3,4,朱權(quán)潔

(1.貴州省礦山安全科學(xué)研究院,貴陽 550025;2.貴州省煤礦設(shè)計(jì)研究院,貴陽 550025;3.貴州省煤礦瓦斯防治工程技術(shù)研究中心,貴陽 550025；4.貴州安和礦業(yè)科技工程股份有限公司, 貴陽 550025)

?

煤礦井下爆破作業(yè)頻帶能量分布規(guī)律*

肖鑄1,2,張鵬翔1,2,3,李青松1,2,3,4,朱權(quán)潔1,2

(1.貴州省礦山安全科學(xué)研究院,貴陽 550025;2.貴州省煤礦設(shè)計(jì)研究院,貴陽 550025;3.貴州省煤礦瓦斯防治工程技術(shù)研究中心,貴陽 550025；4.貴州安和礦業(yè)科技工程股份有限公司, 貴陽 550025)

摘要:為了有效降低煤礦井下爆破作業(yè)擾動(dòng),減少爆破誘發(fā)的礦震和煤與瓦斯突出等次生災(zāi)害,需要對(duì)煤礦井下爆破作業(yè)能量變化規(guī)律進(jìn)行研究。采用錨桿外掛式方法安裝傳感器,優(yōu)化了微震監(jiān)測系統(tǒng)安裝方法。對(duì)煤礦井下爆破作業(yè)過程的微震信號(hào)進(jìn)行采集,并應(yīng)用小波包分解方法分析其頻譜,發(fā)現(xiàn)煤礦井下爆破作業(yè)信號(hào)的能量頻帶分布范圍廣,各個(gè)頻帶內(nèi)爆破微震信號(hào)的能量百分比不斷變化。爆破微震信號(hào)92%以上的能量集中在0～500 Hz的低頻帶范圍內(nèi)。隨傳播距離的增加,微震信號(hào)的高頻部分衰減速度快。在煤礦現(xiàn)場應(yīng)用低爆速、低密度、小直徑的炸藥,降低最大一段的炸藥藥量,采用分散布藥與不耦合裝藥方式有利于降低地震效應(yīng)。

關(guān)鍵詞:煤層；爆破；微震；頻譜分析

煤礦井下爆破作業(yè)時(shí),能量較大,擾動(dòng)影響大,對(duì)含瓦斯煤巖體破壞嚴(yán)重。井下系統(tǒng)環(huán)境中瓦斯運(yùn)移場、地應(yīng)力場受到破壞后失衡,導(dǎo)致煤巖體受力不均衡,瓦斯運(yùn)移過程突變,導(dǎo)致發(fā)生頂板冒落、煤與瓦斯突出等事故。根據(jù)多年來煤礦事故統(tǒng)計(jì)分析,發(fā)現(xiàn)所有的作業(yè)過程中,爆破作業(yè)最易誘發(fā)頂板事故、煤與瓦斯突出事故的作業(yè)。尤其是貴州中小煤礦多數(shù)煤礦采用炮掘作業(yè)方式進(jìn)行掘進(jìn),因此,本文引入微震手段對(duì)煤礦井下煤巖體狀態(tài)進(jìn)行監(jiān)測監(jiān)控,目的是研究煤巖體受爆破作業(yè)影響后的變化、失衡、破壞過程,有效降低爆破后的振動(dòng)效應(yīng),減少煤礦事故的發(fā)生。

1概況

微震監(jiān)測系統(tǒng)主要包括信號(hào)采集、信號(hào)接收系統(tǒng)、信號(hào)處理系統(tǒng)三大部分。微震監(jiān)測技術(shù)最早在井下工程中應(yīng)用于貴重金屬礦山開采誘發(fā)的地震。截止目前為止,在國內(nèi)礦山主要應(yīng)用在監(jiān)測礦柱殘存地壓、沖擊地壓、煤與瓦斯突出、巖爆監(jiān)測預(yù)警、巖體穩(wěn)定性等各個(gè)方面,并取得較好的效果。相對(duì)于傳統(tǒng)的監(jiān)測手段,微震技術(shù)可以監(jiān)測煤層內(nèi)部微震事件,追蹤煤層深部受擾動(dòng)地質(zhì)體裂隙發(fā)育、發(fā)展和突出的變化過程和受采動(dòng)影響區(qū)域內(nèi)的地應(yīng)力的重新分布過程。

井下巷道掘進(jìn)爆破時(shí),大部分的爆破能量在破壞和粉碎煤巖體的過程中被耗散,小部分能量傳播出去[1-5]。傳播的過程是一個(gè)能量不斷衰減的過程。由于煤巖體物理力學(xué)性質(zhì)不均勻,且含有不連續(xù)的結(jié)構(gòu)面,因此,在介質(zhì)中產(chǎn)生的波動(dòng)現(xiàn)象十分復(fù)雜。一方面,煤巖體介質(zhì)中波的能量隨著反射和投射減少,另一方面,隨著內(nèi)摩擦現(xiàn)象而被吸收,并使煤巖體積蓄一部分彈性能。爆破在煤巖體中發(fā)生后,產(chǎn)生沖擊波、應(yīng)力波、地震波。其中,沖擊波使煤巖體破碎,應(yīng)力波使煤巖體指點(diǎn)發(fā)生位移和應(yīng)變、擴(kuò)張和應(yīng)變,地震波引起裂隙的發(fā)育。沖擊波破壞作用最大,傳播距離最短；應(yīng)力波、地震波作用范圍較大,時(shí)間較長。經(jīng)筆者對(duì)煤礦井下現(xiàn)場爆破作業(yè)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)表明,與露天臺(tái)階式爆破不同,煤礦井下爆破微震波信號(hào)時(shí)間為250 ms以上。微震監(jiān)控傳感器布置于應(yīng)力波、地震波的作用范圍內(nèi)。

2煤礦井下爆破能量

煤與瓦斯突出事故相關(guān)資料統(tǒng)計(jì)分析研究結(jié)果表明,煤礦井下爆破作業(yè)是煤與瓦斯突出的最主要誘導(dǎo)方式,掘進(jìn)工作面是煤與瓦斯突出事故發(fā)生的主要地點(diǎn)區(qū)域。煤礦井下爆破震動(dòng)的危害程度受到爆破震動(dòng)幅值、頻譜、持續(xù)時(shí)間等多參數(shù)的作用影響。而且與爆破次數(shù)、爆破前應(yīng)力狀態(tài)、周圍煤巖體特性等相關(guān)。井下爆破作業(yè)釋放的能量大部分在破壞、拋擲煤巖體的過程中耗散掉,僅有一小部分在傳播介質(zhì)中轉(zhuǎn)化為煤巖體、空氣等介質(zhì)中傳送的地震波。地震波能量占總能量的百分比因傳播介質(zhì)不同而迥異[6],大部分介質(zhì)中能量占比在3%～20%之間[7]。在均勻介質(zhì)體中,爆破振動(dòng)地震波傳播過程中能量消耗較非均勻介質(zhì)少,引起的振動(dòng)效應(yīng)影響大。

煤礦井下掘進(jìn)巷道爆破作用在炮孔中的最大壓力與爆破炸藥的性質(zhì)、密度、爆速以及爆破時(shí)的裝藥結(jié)構(gòu)有關(guān),在時(shí)間歷程上,先按照指數(shù)規(guī)律升壓然后再降壓。影響爆破振動(dòng)效應(yīng)的因素包括炸藥量、炮孔個(gè)數(shù)、炮孔布置方式、單段最大藥量、傳播距離、煤巖體性質(zhì)等[8]。爆破時(shí),起爆藥量、藥包形狀、炸藥特性、煤巖體性質(zhì)都會(huì)影響爆破地震波的頻譜和幅值特性。

地震波在煤巖體介質(zhì)中傳播過程中,距離越長,發(fā)射、折射次數(shù)越多,衰減幅度越大,頻率越少越小,相應(yīng)的幅值也變小。地震波速度V與最大段藥量Q、水平距離R,地形、地質(zhì)條件、煤巖體性質(zhì)等有關(guān)的系數(shù)K、a的相互間關(guān)系可以利用薩道夫斯基公式V=K(Q1/3/R)a表示。煤礦井下爆破掘進(jìn)過程中多采用延長藥包法,裝藥結(jié)構(gòu)采用炮孔法,起爆方法采用微差爆破。微差爆破能有效降低地震效應(yīng)。根據(jù)利文斯頓彈性變形方程,以煤巖體在炸藥包臨界深度破壞時(shí),藥包臨界埋深Le與煤巖體彈性變形系數(shù)Eb、臨界藥量Qe的關(guān)系為Le=Eb(Qe)1/3。大部分煤體抗壓強(qiáng)度為1～25 MPa、縱波波速在1500～2400 km/s、彈性模量為1.5～4.5 GPa,孔隙率為2.5%～8.0%。與巖體相比,煤體抗壓強(qiáng)度弱、波速慢、彈性模量小、孔隙率高。炸藥的密度與爆速高時(shí),對(duì)煤巖體的振動(dòng)強(qiáng)度增大。在煤礦井下選用低爆速、低密度、小直徑的炸藥,有利于降低地震效應(yīng)。降低爆破中最大一段的炸藥藥量,可有效降低振速。另外,分散布藥與不耦合裝藥可控制振動(dòng)能量的幅值。

3爆破作業(yè)微震信號(hào)的小波包分析

3.1微震傳感器布設(shè)方式

為有效采集煤礦井下巷道爆破掘進(jìn)過程中的微震信號(hào),根據(jù)義忠煤礦生產(chǎn)布局以及井下采掘巷道開挖情況,設(shè)計(jì)將微震監(jiān)測傳感器布設(shè)于11112掘進(jìn)工作面迎頭后方。煤體中微震信號(hào)傳播速度為2～5 m/ms,傳感器信號(hào)采集頻率為4000 Hz,單次爆破作業(yè)微震信號(hào)長度為200～800 ms之間,因此,將微震傳感器布設(shè)間距設(shè)計(jì)為10 m以上,在巷道兩幫依次布設(shè)。

以往傳感器安裝時(shí),采用紙杯裝錨固劑將傳感器黏結(jié)在安裝鉆孔中。安裝方式存在費(fèi)時(shí)、費(fèi)力、工程量大,安裝過程中易導(dǎo)致傳感器損壞,以及回收困難等弊端。在現(xiàn)場應(yīng)用過程中,克堅(jiān)攻難,對(duì)存在的問題采用專家頭腦風(fēng)暴法等現(xiàn)代科研方法進(jìn)行了技術(shù)攻關(guān),不斷改進(jìn)微震系統(tǒng)布設(shè)方法、埋設(shè)方式等。在義忠煤礦微震監(jiān)測試驗(yàn)過程中自主研發(fā)出了錨桿式和可伸張埋入式兩種安裝方式,在現(xiàn)場應(yīng)用中效果良好。傳感器錨桿式安裝示意圖及連接器見圖1。

3.2分析方法的選取

信號(hào)波形分析中有傅里葉變換、小波變換、小波包變換等各類不同的方法。傅里葉變換能有效表征信號(hào)全局的性質(zhì),但對(duì)局部細(xì)節(jié)內(nèi)容無法深入。小波變換的優(yōu)點(diǎn)在于能對(duì)信號(hào)尤其是低頻信號(hào)進(jìn)行較好的表述,但對(duì)高頻部分的細(xì)節(jié)無法繼續(xù)分解。小波包分解方法能對(duì)高頻、低頻部分分別進(jìn)行進(jìn)一步精細(xì)的分解,能對(duì)波形進(jìn)行有效的時(shí)頻分辨。井下爆破微震波信號(hào)頻率豐富,既含有高頻成分,又含有低頻成分,因此,選用小波包分析方法對(duì)井下爆破微震信號(hào)進(jìn)行研究。

常用的小波函數(shù)中,Daubechies小波具有時(shí)域有限、整數(shù)位移后政教歸一等特點(diǎn),小波隨階次增大,消失矩階數(shù)變高,越有利于更好的劃分頻帶,但計(jì)算效率降低,時(shí)域緊支撐性變?nèi)?。結(jié)合煤礦井下爆破微震信號(hào)的波形特點(diǎn),選擇Db12小波包進(jìn)行分析。小波包分解到n層,有2^n種分解方式,選用度量信息規(guī)律性的Shannon熵、對(duì)數(shù)能量熵、SURE熵等作為指導(dǎo)小波包分解的特征函數(shù),降低運(yùn)算量,挑出最優(yōu)小波包樹。

3.3煤礦井下爆破作業(yè)微震信號(hào)小波包分析

井下巷道爆破作業(yè)時(shí),產(chǎn)生含有豐富頻率的微震波,在煤巖體中傳播一段距離后,能被微震傳感器捕捉。通過采用小波包分析方法對(duì)微震信號(hào)進(jìn)行處理、分析,可以得到爆破微震信號(hào)的統(tǒng)計(jì)規(guī)律[9-12]。選用加拿大ESG微震監(jiān)測監(jiān)控系統(tǒng)對(duì)煤礦井下爆破作業(yè)微震信號(hào)進(jìn)行采集。爆破微震信號(hào)采集采用寬頻帶高效率加速度傳感器,采集最小工作頻率為50 Hz。微震信號(hào)的采樣頻率為4000 Hz,根據(jù)Nyquist-Shannon采樣定理,可知其奈奎斯頻率為2000 Hz。利用小波包對(duì)信號(hào)進(jìn)行4層分解,第4層小波包分解樹具有16個(gè)節(jié)點(diǎn),微震爆破信號(hào)在頻域范圍內(nèi)被劃分為16個(gè)子頻帶。各個(gè)頻帶寬間隔125 Hz,最低頻帶為0～125 Hz,最高頻帶為1875～2000 Hz。根據(jù)各種小波基的緊支性、光滑性和近似對(duì)成性,選擇有利于分析爆破微震非平穩(wěn)信號(hào)的Daubechies函數(shù),為更有利于劃分頻帶,選擇db12序列作為微震信號(hào)小波包分析基函數(shù)。

在煤礦井下炮掘巷道進(jìn)行了10 d 30次的現(xiàn)場爆破作業(yè)試驗(yàn)。從現(xiàn)場試驗(yàn)采集到的爆破微震信號(hào)中選取具有代表性的25個(gè)進(jìn)行4層小波包分解。運(yùn)用matlab編程,對(duì)第4層小波包分解系數(shù)進(jìn)行平方求和,得出各個(gè)節(jié)點(diǎn)的能量值,并進(jìn)行歸一化處理,得出不同頻帶的能量占總能量的比例。隨著掘進(jìn)工作面的不斷向前推進(jìn),主要頻帶的能量變化如圖2所示。

由圖2可以看出煤礦井下爆破微震信號(hào)的頻率豐富,在0～2000 Hz的各個(gè)頻帶范圍內(nèi)均有分布。低頻信號(hào)能量略為分散,其中125～250 Hz、375～500 Hz兩個(gè)頻帶范圍內(nèi)信號(hào)能量占比較大。高頻信號(hào)能量較為集中,主要分布在750～875 Hz、1500～1625 Hz的頻帶范圍內(nèi)?？梢?不同于一般信號(hào),可以采用單個(gè)主頻來反應(yīng)整個(gè)頻域信號(hào)的特征,在微震爆破信號(hào)分析的過程中,需要將不同頻帶的微震能量綜合考慮進(jìn)去。

由于井下爆破作業(yè)震源以及傳播路徑的復(fù)雜和不重復(fù)性,每次爆破產(chǎn)生的微震信號(hào)各個(gè)頻帶的能量占總能量的百分比也起伏不定。由圖2可以看出,能量分布中百分比最大的兩個(gè)頻帶為125～250 Hz、375～500 Hz頻帶,爆破產(chǎn)生微震信號(hào)的頻帶能量百分比在不斷變化。綜上可以看出,在不同的震源和傳播路徑條件下,爆破微震信號(hào)的主要頻率及其能量也在不斷發(fā)生變化。

爆破微震信號(hào)能量頻帶分布廣,但小于500 Hz的頻帶能量占比在92%以上,大于500 Hz的頻帶范圍內(nèi)信號(hào)能量占比在8%以下,可見義忠煤礦井下一次正向爆破中,微震波能量大部分集中在0～500 Hz的低頻頻帶范圍內(nèi)。

500 Hz以上的頻帶范圍中,微震能量百分比在8%以下,能量占比較小,主要集中在750～875 Hz、1500～1625 Hz的頻帶范圍內(nèi)。隨著掘進(jìn)巷道的不斷推進(jìn),爆破微震信號(hào)的傳播距離逐漸增大,高頻段信號(hào)能量百分比迅速減小,表明煤體中爆破震動(dòng)信號(hào)的高頻部分衰減速度快。

4小結(jié)與展望

(1)煤礦井下作業(yè)工序中,爆破作業(yè)釋放能量大,對(duì)煤礦井下采掘系統(tǒng)、煤巖層等造成劇烈的擾動(dòng),易引發(fā)頂板垮落,誘發(fā)煤與瓦斯突出等,需要對(duì)爆破作業(yè)及其前后的環(huán)境變化進(jìn)行實(shí)時(shí)監(jiān)測監(jiān)控,嚴(yán)防煤礦重大事故的發(fā)生,減少人員傷亡。

(2)煤礦井下煤巖體爆破過程中產(chǎn)生的能量大部分耗散在破碎、拋擲煤巖體的過程中,僅有約3%～20%的能量轉(zhuǎn)化為地震波進(jìn)行擴(kuò)散。地震波的傳播效應(yīng)受到炸藥量、爆破參數(shù)、煤巖體性質(zhì)、傳播距離、地質(zhì)構(gòu)造環(huán)境等因素的制約。在井下選用低爆速、低密度、小直徑的炸藥,降低最大一段的炸藥藥量,采用分散布藥與不耦合裝藥方式有利于降低地震效應(yīng)。

(3)通過合理設(shè)定傳感器安裝距離等方法優(yōu)化了傳感器布置方案；自主研發(fā)了錨桿外掛式、可伸張埋入式等不同傳感器安裝方法,有效減少了工程量,降低了傳感器損耗,保證了傳感器的回收和循環(huán)利用；使監(jiān)測監(jiān)控系統(tǒng)能有效采集煤礦井下爆破作業(yè)信號(hào)。

(4)煤礦井下爆破作業(yè)信號(hào)的能量頻帶分布范圍廣,各個(gè)頻帶內(nèi)爆破微震信號(hào)的能量百分比不斷變化。爆破微震信號(hào)92%以上的能量集中在0～500 Hz的低頻帶范圍內(nèi)。隨傳播距離的增加,微震信號(hào)的高頻部分衰減速度快。

參考文獻(xiàn)(References)

[1]薛振舉.煤礦斜井掘進(jìn)爆破振動(dòng)監(jiān)測與分析[J].爆破,2013,30(1):70-74.

[1]XUE Zhen-ju.Monitoring and analysis of blasting vibration induced by inclined shaft advancement in coal mine[J].Blasting,2013,30(1):70-74.(in Chinese)

[2]盧文波,張樂,周俊汝,等.爆破振動(dòng)頻率衰減機(jī)制和衰減規(guī)律的理論分析[J].爆破,2013,30(2):1-6,11.

[2]LU Wen-bo,ZHANG Le,ZHOU Jun-ru,et al.Theoretical analysis on decay mechanism and law of blasting vibration frequency[J].Blasting,2013,30(2):1-6,11.(in Chinese)

[3]宗琦,汪海波,周勝兵.爆破地震效應(yīng)的監(jiān)測和控制技術(shù)研究[J].巖石力學(xué)與工程學(xué)報(bào),2008(5):938-945.

[3]ZONG Qi,WANG Hai-bo,ZHOU Sheng-bing.Research on monitoring and conterolling techniquesconsidering effects of seismic shock[J].Chinese Journal of Rock Mechanics and Engineering,2008(5):938-945.(in Chinese)

[4]褚懷保.煤體爆破作用機(jī)理及試驗(yàn)研究[D].焦作市：河南理工大學(xué),2011.

[5]朱權(quán)潔,姜福興,于正興,等.爆破震動(dòng)與巖石破裂微震信號(hào)能量分布特征研究[J].巖石力學(xué)與工程學(xué)報(bào),2012(4):723-730.

[5]ZHU Quan-jie,JIANG Fu-xing,YU Zheng-xing,et al.Study on energy distribution characters about blasting vibration and rock fracture microseismic signal[J].Chinese Journal of Rock Mechanics and Engineering,2012(4):723-730.(in Chinese)

[6]石崇.爆破地震效應(yīng)分析與安全評(píng)價(jià)[D].青島：山東科技大學(xué),2005.

[7]馮叔瑜,朱節(jié)忠,朱乃耀.土巖工程大量爆破及其設(shè)計(jì)[M].北京：冶金工業(yè)出版社,1968：382,384-385.

[8]范孝鋒,周傳波,陳國平.爆破震動(dòng)影響因素的灰關(guān)聯(lián)分析[J].爆破,2005,22(2)：100-102.

[8]FAN Xiao-feng,ZHOU Chuan-bo,CHEN Guo-ping.The influential factors of blasting vibration by grey correlation analysis[J].Blasting,2005,22(2):100-102.(in Chinese)

[9]丁幼亮,李愛群,繆長青.基于小波包能量譜的結(jié)構(gòu)損傷預(yù)警方法研究[J].工程力學(xué),2006(8):42-48.

[9]DING You-liang,LI Ai-qun,MIAO Chang-qing.Investigation on the structural damage alarming method based on wavelet packet energy spectrum[J].Engineering Mechanics,2006(8):42-48.(in Chinese)

[10]李明愛,馬建勇,楊金福.基于小波包和熵準(zhǔn)則的最優(yōu)頻段提取方法[J].儀器儀表學(xué)報(bào),2012(8):1721-1728.

[10]LI Ming-ai,MA Jian-yong,YANG Jin-fu.Optimal frequency band extraction method based on wavelet packet and entropy criterion[J].Chinese Journal of Scientific Instrument,2012(8):1721-1728.(in Chinese)

[11]杜立志,殷琨,張曉培,等.基于最優(yōu)小波包基的降噪方法及其應(yīng)用[J].工程地球物理學(xué)報(bào),2008(1):25-29.

[11]DU Li-zhi,YIN Kun,ZHANG Xiao-pei,et al.The signal denoising algorithm based on optimization wavelet packet and its application[J].Engineering Journal of Geophysics,2008(1):25-29.(in Chinese)

[12]和雪松.用小波包識(shí)別地震和礦震[D].北京市：中國地震局地球物理研究所,2005.

Frequency Band Energy Distribution of Coal Mine Blasting

XIAOZhu1,2,ZHANGPeng-xiang1,2,3,LIQing-song1,2,3,4,ZHUQuan-jie1,2

(1.Guizhou Mine Safety Research Institute,Guiyang 550025,China;2.Guizhou Coalmine Design & Research Institute,Guiyang 550025,China;3.Guizhou Province Engineering Technology Research Center for Coal Gas Prevention ﹠ Control,Guiyang 550025,China;4.Guizhou AnHe mining technology and engineering Co Ltd,Guiyang 550025,China)

Abstract:In order to reduce the disturbance caused by coal mine underground blasting,reduce the secondary disasters such as mine earthquake and coal and gas outburst,need to study the energy variation law of coal mine blasting operation.Install the sensor using bolt plug-in method,installation method to optimize the microseismic monitoring system.Coal mine blasting process of microseismic signal acquisition and application of wavelet packet decomposition method for the analysis of the spectrum,it is found that the coal mine blasting signal frequency band energy distribution of a wide range of percentage of energy in each frequency band of blasting seismic signals constantly changing.More than 92% of the energy concentration of blasting seismic signals in the 0～500 Hz frequency band.With the increase of propagation distance,the high frequency part of microseismic signal attenuation speed.In the coal mine field,the low detonation velocity,the low density,the small diameter of the explosive,and the lower the maximum amount of explosive charge,and the distribution of the drug and the non coupling charging method are used to reduce the seismic effect.

Key words:coal seam；blasting；microseismic；spectrum analysis

doi:10.3963/j.issn.1001-487X.2016.02.015

收稿日期:2016-04-06

作者簡介:肖鑄(1965-)，男，高級(jí)工程師、學(xué)士，主要從事煤礦設(shè)計(jì)、瓦斯防治、安全工程與技術(shù)等方面的工程技術(shù)服務(wù)和科研工作，(E-mai)349635236@qq.com。 通訊作者:張鵬翔(1987- )，男，工程師、碩士，主要從事瓦斯防治、煤礦設(shè)計(jì)與咨詢、安全技術(shù)及工程等方面的工程技術(shù)服務(wù)和科研工作，(E-mail)244492525@qq.com。

基金項(xiàng)目：國家“十二五”科技支撐計(jì)劃項(xiàng)目(2012BAK04B07)；貴州省社會(huì)發(fā)展攻關(guān)計(jì)劃項(xiàng)目(黔科合社G字[2011]4003、黔科合SY字[2010]3036)

中圖分類號(hào):TD713

文獻(xiàn)標(biāo)識(shí)碼:A

文章編號(hào):1001-487X(2016)02-0078-05