陳志英，周小娜，盧超龍

(1.福建省高電壓技術(shù)重點(diǎn)實(shí)驗(yàn)室，福建 廈門 361024；2.廈門國(guó)毅科技有限公司，福建 廈門 361006 )

?

高壓斷路器合(分)閘線圈電流在線監(jiān)測(cè)系統(tǒng)的研制

陳志英1，周小娜1，盧超龍2

(1.福建省高電壓技術(shù)重點(diǎn)實(shí)驗(yàn)室，福建 廈門 361024；2.廈門國(guó)毅科技有限公司，福建 廈門 361006 )

針對(duì)目前高壓斷路器在線監(jiān)測(cè)系統(tǒng)的研究現(xiàn)狀與不足，為及時(shí)發(fā)現(xiàn)運(yùn)行中的高壓斷路器故障，提高運(yùn)行可靠性，提出了一種基于DSP的高壓斷路器合(分)閘線圈電流在線監(jiān)測(cè)系統(tǒng)的設(shè)計(jì)方案.該方案以TMS320C28346為系統(tǒng)控制核心，利用AD7606實(shí)時(shí)采集斷路器合(分)閘線圈電流信號(hào)，采用包絡(luò)均值法去除波形中的噪聲和干擾，提取線圈電流波形中的關(guān)鍵時(shí)間點(diǎn)和電流最大值，計(jì)算合(分)閘時(shí)間，根據(jù)電流最大值與合(分)閘時(shí)間是否在正常范圍內(nèi)進(jìn)行故障判斷，并將電流最大值、合(分)閘時(shí)間、日期和運(yùn)行狀態(tài)顯示在LCD上.實(shí)驗(yàn)結(jié)果證明，該方案可實(shí)現(xiàn)斷路器合(分)閘線圈電流的實(shí)時(shí)監(jiān)測(cè)，可及時(shí)發(fā)現(xiàn)故障，具有實(shí)時(shí)性能好、精度高等優(yōu)點(diǎn)，合閘線圈電流最大值誤差在1.29%，合閘時(shí)間誤差在0.26%；分閘線圈電流最大值誤差在1.61%，分閘時(shí)間誤差在0.7%.

高壓斷路器；合(分)閘線圈電流；在線監(jiān)測(cè)；TMS320C28346；包絡(luò)均值法

在電力系統(tǒng)的一次設(shè)備中，高壓斷路器是僅次于發(fā)電機(jī)和變壓器的大型電力設(shè)備.當(dāng)系統(tǒng)正常運(yùn)行時(shí)，它能切斷和接通線路及各種電氣設(shè)備；當(dāng)系統(tǒng)發(fā)生故障時(shí)，它和繼電器保護(hù)配合，迅速切斷故障電流，防止事故范圍擴(kuò)大.若高壓斷路器無(wú)法正常運(yùn)行，會(huì)造成線路、設(shè)備受損，電量損失，用戶大面積停電，影響正常的社會(huì)生產(chǎn)與生活[1].因此，高壓斷路器作為電力系統(tǒng)中的重要控制與保護(hù)設(shè)備，對(duì)其進(jìn)行維護(hù)與檢修以確保它正常運(yùn)行至關(guān)重要[2-3].

對(duì)斷路器在線監(jiān)測(cè)可以及時(shí)發(fā)現(xiàn)故障，提高運(yùn)行可靠性.合(分)閘線圈是斷路器操作機(jī)構(gòu)的重要組成部分，線圈電流波形中包含了斷路器運(yùn)行狀態(tài)的關(guān)鍵信息，如線圈卡滯、輔助觸點(diǎn)狀態(tài)等[4-7].因此，通過(guò)對(duì)高壓斷路器的合(分)閘線圈電流的實(shí)時(shí)監(jiān)測(cè)，分析波形中的關(guān)鍵信息，可判斷斷路器是否出現(xiàn)故障.

目前，現(xiàn)有的高壓斷路器在線監(jiān)測(cè)系統(tǒng)多采用高速采集卡采集信號(hào)，并將采集的數(shù)據(jù)通過(guò)串口傳輸給上位機(jī)進(jìn)行故障診斷，成本較高，由于需傳輸?shù)臄?shù)據(jù)量較大，實(shí)時(shí)性也受到影響.針對(duì)上述缺點(diǎn)，本文提出了一種以TMS320C28346為控制核心的合(分)閘線圈電流在線監(jiān)測(cè)系統(tǒng)設(shè)計(jì)方案，包括傳感器的選型、調(diào)理電路和AD采樣電路的設(shè)計(jì)以及軟件設(shè)計(jì)，斷路器的合(分)閘線圈電流信號(hào)經(jīng)過(guò)傳感器、調(diào)理電路以及A/D轉(zhuǎn)換后輸入到控制核心進(jìn)行分析、處理，實(shí)現(xiàn)電流信號(hào)的實(shí)時(shí)監(jiān)測(cè)、合閘時(shí)間等參數(shù)的計(jì)算以及故障判斷.

1 斷路器合(分)閘線圈電流信號(hào)處理與故障診斷

1.1 合(分)閘線圈電流

斷路器的合(分)閘線圈電流信號(hào)是判斷斷路器是否正常運(yùn)行的重要依據(jù)[8-9]，通過(guò)分析線圈電流波形可得到觸頭合(分)閘動(dòng)作啟動(dòng)時(shí)刻、合(分)閘時(shí)間等參數(shù)，其典型波形如圖1所示.圖1中線圈電流變化總共經(jīng)歷了4個(gè)階段：t0～t1為鐵芯觸發(fā)階段，線圈電流不斷增大，t0代表斷路器合(分)閘動(dòng)作的觸發(fā)時(shí)刻；t1～t2為鐵芯運(yùn)動(dòng)階段，電流逐漸減小，可根據(jù)這段波形判斷鐵芯的運(yùn)動(dòng)狀態(tài)是否正常；t2～t3為觸頭合(分)閘動(dòng)作階段，開(kāi)始時(shí)刻在t2附近，線圈電流在t3時(shí)達(dá)到穩(wěn)定值，此時(shí)斷路器的動(dòng)觸頭動(dòng)作到位；t3～t4為電流切斷階段，輔助開(kāi)關(guān)觸頭間的電弧被拉長(zhǎng)，電弧電壓快速升高，迫使電流迅速減小直至熄滅[10].因此，t0、t1、t2、t3、t4、I3代表了斷路器在動(dòng)作過(guò)程中線圈電流信號(hào)的關(guān)鍵點(diǎn)，提取這些關(guān)鍵點(diǎn)，并判斷電流的最大值與動(dòng)作時(shí)間是否在正常范圍內(nèi)，即可判斷該斷路器操動(dòng)機(jī)構(gòu)是否正常運(yùn)行[11-13].

1.2 包絡(luò)均值法

高壓斷路器在運(yùn)行過(guò)程中存在一些無(wú)法避免的干擾信號(hào)，為提高在線監(jiān)測(cè)系統(tǒng)的精度，正確提取合(分)閘線圈電流的特征信號(hào)，需對(duì)采集到的波形進(jìn)行濾波、平滑等處理.包絡(luò)均值法在保證波形不失真的前提下可以有效濾除干擾信號(hào)，其過(guò)程如下：

1)求極大值點(diǎn)對(duì)采集到的數(shù)據(jù)I(t)進(jìn)行處理，首先找出數(shù)據(jù)中所有的極大值點(diǎn)(X)以及極大值點(diǎn)所對(duì)應(yīng)的時(shí)間(Y)，組成一組新的節(jié)點(diǎn)，如式(1)所示：

(1)

2)用三次樣條插值函數(shù)擬合形成原數(shù)據(jù)的上包絡(luò)線在每個(gè)分段區(qū)間(Xi，Xi+1)內(nèi)求解三次多項(xiàng)式函數(shù)Si(x)，其一階導(dǎo)數(shù)和二階導(dǎo)數(shù)在區(qū)間(X0，Xn)內(nèi)是連續(xù)的，因此n個(gè)三次多項(xiàng)式的表達(dá)式如式(2)所示：

(2)

其中，ai、bi、ci、di代表4n個(gè)未知系數(shù)，表達(dá)式如(3)所示，i=0,1,…,n-1(已標(biāo)注的除外).

(3)

在求解mi時(shí)，其方程僅有n-1個(gè)，而未知數(shù)有n+1個(gè)，因此還需對(duì)X0和Xn兩端點(diǎn)的微分加些限制，增加兩個(gè)方程方可解出所有mi的值.常見(jiàn)的限制條件有自由邊界、固定邊界以及非節(jié)點(diǎn)邊界，不同的端點(diǎn)條件會(huì)得到不同的樣條曲線.

3)求極小值點(diǎn)，其過(guò)程與1)相同.

4)用三次樣條插值函數(shù)擬合形成原數(shù)據(jù)的下包絡(luò)線，其過(guò)程與2)相同.

5)取上包絡(luò)線和下包絡(luò)線的均值，得到新的波形.

1.3 斷路器動(dòng)作觸發(fā)時(shí)刻t0的提取方法

高壓斷路器動(dòng)作的觸發(fā)時(shí)刻t0為斷路器接到合(分)閘指令的時(shí)刻，如圖1所示，當(dāng)檢測(cè)到有上升沿信號(hào)時(shí)，即認(rèn)為斷路器接到指令.系統(tǒng)選用±10V作為采樣芯片AD7606的輸入范圍，其最低有效位為305.175μV，精度為0.03%，即：采樣最小電壓變化范圍ΔU=1V×0.03%=0.003V，而系統(tǒng)選用的霍爾電流傳感器CSM010GT的精度為0.7%，即：轉(zhuǎn)換最小電壓變化范圍ΔU=1V×0.07%=0.007V.因此，當(dāng)線圈電流經(jīng)轉(zhuǎn)換后的電壓變化大于0.007V時(shí)，認(rèn)為線圈有電流流過(guò)，并將此時(shí)刻定義為斷路器動(dòng)作觸發(fā)時(shí)刻t0，系統(tǒng)開(kāi)始保存采樣數(shù)據(jù).

1.4 剛合(分)點(diǎn)的提取方法

高壓斷路器的合(分)閘時(shí)間是斷路器的一個(gè)重要參數(shù)，其定義為斷路器合(分)閘動(dòng)作觸發(fā)時(shí)刻到動(dòng)靜觸頭剛剛接觸(或分離)的時(shí)間，而剛分、剛合點(diǎn)的提取一直是斷路器在線監(jiān)測(cè)分析的難點(diǎn).目前，剛合(分)點(diǎn)的提取方法主要有振動(dòng)信號(hào)分析法、觸頭加速度分析法等[14-17]，上述方法精確度較高，但實(shí)現(xiàn)過(guò)程復(fù)雜.本文采用輔助開(kāi)關(guān)的開(kāi)斷確定剛分、剛合點(diǎn)[18]，在斷路器動(dòng)作過(guò)程中，主動(dòng)軸帶動(dòng)主觸頭旋轉(zhuǎn)到指定位置，輔助開(kāi)關(guān)的狀態(tài)發(fā)生變化，輔助開(kāi)關(guān)狀態(tài)變化時(shí)刻，即為剛分、剛合時(shí)刻.1.5 故障診斷方法

同種型號(hào)的高壓斷路器的線圈電流波形具有很高的重復(fù)性，通過(guò)多次實(shí)驗(yàn)，可確定斷路器正常工作時(shí)的合(分)閘時(shí)間與電流最大值的范圍.當(dāng)線圈出現(xiàn)卡滯、匝間短路等故障時(shí)，合(分)閘時(shí)間與電流最大值會(huì)出現(xiàn)明顯變化.因此，可采用判斷合(分)閘時(shí)間與電流最大值是否超出正常范圍的方法進(jìn)行故障診斷.通過(guò)斷路器合(分)閘線圈電流在線監(jiān)測(cè)系統(tǒng)實(shí)時(shí)測(cè)量線圈電流，經(jīng)分析計(jì)算后，若線圈合(分)閘時(shí)間或電流最大值不在正常范圍內(nèi)時(shí)，則判斷斷路器出現(xiàn)故障.

2 系統(tǒng)硬件設(shè)計(jì)

斷路器合(分)閘線圈電流在線監(jiān)測(cè)系統(tǒng)的控制原理如圖2.系統(tǒng)上電后對(duì)運(yùn)行中斷路器的合(分)閘電流進(jìn)行不間斷實(shí)時(shí)監(jiān)測(cè)，合(分)閘線圈電流信號(hào)經(jīng)霍爾電流傳感器將電流信號(hào)轉(zhuǎn)換成電壓信號(hào)，通過(guò)調(diào)理電路處理后，將傳感器輸出電壓的范圍轉(zhuǎn)換到AD采樣芯片的工作電壓范圍內(nèi)，AD采樣芯片將模擬信號(hào)轉(zhuǎn)換為數(shù)字信號(hào)，并輸入到控制核心TMS320C28346，經(jīng)計(jì)算分析后，將監(jiān)測(cè)結(jié)果顯示在LCD上.其中，電源電路為控制核心以及其他部件提供穩(wěn)定電源，通過(guò)按鍵電路可選擇顯示界面及內(nèi)容.

2.1 信號(hào)調(diào)理電路

系統(tǒng)選用CSM010GT型霍爾電流傳感器將線圈電流轉(zhuǎn)換成電壓信號(hào)以便采集，該傳感器體積小、精度高、線性度好，若僅給傳感器供電，原邊無(wú)電流流過(guò)時(shí)，傳感器輸出基準(zhǔn)電壓2.5V.根據(jù)傳感器的技術(shù)參數(shù)，得到輸入輸出關(guān)系為：

(4)

其中：Iin為輸入電流，即線圈電流；Vout為傳感器輸出電壓.

高壓斷路器的合(分)閘線圈電流一般為0～5 A，則由式(4)可知傳感器的輸出電壓范圍為2.5～3 V.為滿足AD采樣芯片的分辨率與采樣電壓范圍，需將傳感器輸出電壓通過(guò)調(diào)理電路進(jìn)行放大，其原理圖如圖3所示.該放大電路的輸入輸出關(guān)系為：

(5)

2.2 AD采樣電路

為實(shí)現(xiàn)更高精度的采樣，增加準(zhǔn)確性，系統(tǒng)選用高速、低功耗16位的AD7606作為控制核心的外擴(kuò)采樣芯片，可對(duì)8個(gè)模擬輸入通道進(jìn)行同步采樣，每個(gè)通道的吞吐率為200 Kbit/s，且芯片內(nèi)置濾波器，傳感器輸出的電壓可直接接入芯片的模擬信號(hào)輸入引腳，采樣電路如圖4所示.

3 系統(tǒng)軟件設(shè)計(jì)

系統(tǒng)軟件主要包括：AD采樣子程序、合(分)閘時(shí)間計(jì)算子程序、按鍵子程序、LCD12864顯示子程序、RTC(實(shí)時(shí)時(shí)鐘芯片)時(shí)間讀取子程序，其主程序流程圖如圖5所示.上電后系統(tǒng)初始化，調(diào)用AD采樣子程序，不間斷地實(shí)時(shí)采集線圈電流；將采集到的波形經(jīng)包絡(luò)均值法濾波、平滑后，計(jì)算合(分)閘電流的最大值，并提取線圈電流的關(guān)鍵時(shí)間點(diǎn)，調(diào)用合(分)閘時(shí)間計(jì)算子程序，若超出正常范圍，則發(fā)出警報(bào)，調(diào)用RTC時(shí)間讀取子程序記錄合(分)閘動(dòng)作發(fā)生時(shí)間；等待按鍵按下，若有按鍵按下時(shí)，調(diào)用相應(yīng)按鍵子程序與LCD12864顯示子程序.

4 實(shí)驗(yàn)結(jié)果與分析

為驗(yàn)證系統(tǒng)的功能與正確性，在實(shí)驗(yàn)室搭建了斷路器合(分)閘線圈電流在線監(jiān)測(cè)的模擬實(shí)驗(yàn)平臺(tái)，以ABB生產(chǎn)的VD4真空斷路器為實(shí)驗(yàn)對(duì)象，采用220 V直流脫扣器，當(dāng)接到合(分)閘指令后，延時(shí)20 ms后合(分)閘線圈才開(kāi)始有電流，因此，實(shí)際合(分)閘時(shí)間應(yīng)加上20 ms.以廈門國(guó)毅科技有限公司生產(chǎn)的機(jī)械特性儀GYT_07(精度為0.5%)測(cè)量的數(shù)據(jù)為標(biāo)準(zhǔn)對(duì)比數(shù)據(jù)，系統(tǒng)上電后對(duì)斷路器進(jìn)行合閘操作，運(yùn)行一段時(shí)間后，對(duì)斷路器進(jìn)行分閘操作，采集到的合(分)閘線圈電流經(jīng)包絡(luò)均值法濾波，濾波前后的對(duì)比圖如圖6所示.重復(fù)合閘與分閘操作各100次，發(fā)現(xiàn)合(分)閘線圈電流變化趨勢(shì)基本不變，特征值也在小范圍內(nèi)變化，任取其中10次實(shí)驗(yàn)結(jié)果，如表1和表2所示.同時(shí)，通過(guò)模擬斷路器常見(jiàn)的幾種機(jī)械故障，記錄故障時(shí)的合(分)閘線圈電流的監(jiān)測(cè)結(jié)果，脫扣器卡滯和脫扣器卡死的電流最大值分別為1.14 A和1.19 A，分閘時(shí)間分別為59.8 ms和大于100 ms.

實(shí)驗(yàn)結(jié)果表明，包絡(luò)均值法能有效濾除干擾，平滑波形，且保證波形不失真；當(dāng)斷路器正常運(yùn)行時(shí)，合(分)閘線圈電流最大值、合(分)閘時(shí)間在一定的范圍內(nèi)波動(dòng)，與GYT_07測(cè)量的數(shù)據(jù)比較，誤差在可接受范圍內(nèi).以斷路器分閘動(dòng)作為例，當(dāng)出現(xiàn)脫扣器卡滯現(xiàn)象時(shí)，會(huì)造成分閘時(shí)間的增加；當(dāng)出現(xiàn)脫扣器卡死故障時(shí)，斷路器無(wú)法完成分閘動(dòng)作.

表1 斷路器合閘線圈電流與合閘時(shí)間在線監(jiān)測(cè)結(jié)果

表2 斷路器分閘線圈電流與分閘時(shí)間在線監(jiān)測(cè)結(jié)果

5 結(jié)語(yǔ)

本文提出了一種基于DSP的高壓斷路器合(分)閘線圈電流在線監(jiān)測(cè)系統(tǒng)的設(shè)計(jì)方案，通過(guò)多次實(shí)驗(yàn)進(jìn)行驗(yàn)證，得到以下結(jié)論：1)能夠?qū)崟r(shí)監(jiān)測(cè)斷路器合(分)閘線圈電流與合(分)閘時(shí)間，合閘線圈電流最大值誤差在1.29%，合閘時(shí)間誤差在0.26%，分閘線圈電流最大值誤差在1.61%，分閘時(shí)間誤差在0.7%；2)采用包絡(luò)均值法對(duì)電流信號(hào)進(jìn)行處理，實(shí)現(xiàn)波形的濾波與平滑且不失真；3)采用了一種簡(jiǎn)單且易于實(shí)現(xiàn)的剛分、剛合點(diǎn)分析方法，通過(guò)輔助開(kāi)關(guān)的狀態(tài)變化計(jì)算得到合(分)閘時(shí)間；4)根據(jù)合(分)閘時(shí)間和電流最大值是否在正常范圍內(nèi)進(jìn)行故障診斷，可及時(shí)發(fā)現(xiàn)故障；5)記錄動(dòng)作發(fā)生的時(shí)間，并將監(jiān)測(cè)結(jié)果(合(分)閘時(shí)間、電流最大值、日期、是否故障)顯示在LCD上.該方案具有實(shí)時(shí)性好、精度高等優(yōu)點(diǎn)，具有廣闊的市場(chǎng)應(yīng)用前景.

[1]林萃.現(xiàn)代高壓電器技術(shù)[M].北京:機(jī)械工業(yè)出版社,2002:1-7.

[2]張彥飛.斷路器狀態(tài)檢修策略分析[J].中國(guó)高新技術(shù)企業(yè),2014(23):125-126,127.

[3]趙洋,劉漢宇,曾慶軍.高壓真空斷路器機(jī)械特性在線監(jiān)測(cè)系統(tǒng)研制[J].高壓電器,2009,45(6):91-95.

[4]卞和穎,代克杰.基于分合閘線圈信號(hào)的斷路器狀態(tài)在線監(jiān)測(cè)系統(tǒng)設(shè)計(jì)[J].現(xiàn)代電子技術(shù),2014,37(3):100-104.[5]ZEINELDIN H H，EL-SAADANY E F，SALAMA M M A,et al.High voltage circuit breaker modeling for online model-based monitoring and diagnosis [J].IEEE Transactions on Power Delivery,2008,3(4):317-322.

[6]孫霞.高壓斷路器狀態(tài)監(jiān)測(cè)技術(shù)研究[J].中國(guó)科技信息,2014(9):210-211.

[7]張衛(wèi)正,李永麗,姚創(chuàng).基于最小二乘支持向量機(jī)的高壓斷路器故障診斷[J].高壓電器,2015(12):79-83.

[8]陳辰,厲小潤(rùn).基于線圈電流分析的操動(dòng)機(jī)構(gòu)故障診斷方法[J].浙江大學(xué)學(xué)報(bào)(工學(xué)版),2016,50(3):527-535．[9]張永奎,趙智忠,馮旭基,等.基于分合閘線圈電流信號(hào)的高壓斷路器機(jī)械故障診斷[J].高壓電器,2013,49(2):37- 42.

[10]牟濤,周水斌,王曉峰.斷路器機(jī)構(gòu)狀態(tài)監(jiān)測(cè)數(shù)據(jù)采集的實(shí)現(xiàn)[J].低壓電器,2013(3):26-30.

[11]金曉明,邵敏艷,王曉華.基于脫扣線圈電流的斷路器機(jī)械狀態(tài)預(yù)測(cè)算法研究[J].高壓電器,2010,46(4):47-51.

[12]彭鑫霞,關(guān)永剛,高凱,等.應(yīng)用形態(tài)學(xué)濾波方法去除高壓斷路器線圈電流信號(hào)的噪聲[J].高壓電器,2015,51(7):1-7.

[13]孫銀山,張文濤,張一茗,等.高壓斷路器分合閘線圈電流信號(hào)特征提取與故障判別方法研究[J].高壓電器,2015,51(9):134-139.

[14]廖姍姍,黃煒昭,呂啟深,等.真空斷路器在線狀態(tài)監(jiān)測(cè)技術(shù)[J].自動(dòng)化應(yīng)用,2015(12):108-109,112.

[15]黃新波,陶晨,劉斌.智能斷路器機(jī)械特性在線監(jiān)測(cè)技術(shù)和狀態(tài)評(píng)估[J].高壓電器,2015,51(3):129-134.

[16]駱虎,戴冬云,殷曉剛,等.一種基于觸頭受力監(jiān)測(cè)實(shí)現(xiàn)斷路器斷口判定的方法[J].高壓電器,2015,51(9):116-122.

[17]李勁彬,全江濤,陳雋,等.基于分合閘線圈電流的某換流站開(kāi)關(guān)故障分析[J].高壓電器,2015,51(2):141-145.

[18]李大衛(wèi),徐黨國(guó),孫云生,等.基于分合閘線圈電流的斷路器缺陷診斷及試驗(yàn)方法研究[J].高壓電器,2015,51(8):114-118.

(責(zé)任編輯 雨 松)

Design of Online Coil Current Monitoring System forHigh Voltage Circuit Breaker

CHEN Zhiying1，ZHOU Xiaona1，LU Chaolong2

(1.Key High-Voltage Laboratory of Fujian Province，Xiamen 361024，China;2.Xiamen Guoyi Technology Co.Ltd.,Xiamen 361024，China)

Presented in this paper is the design of an online monitoring system for coil current of high voltage circuit breaker based on TMS320C28346，aiming to discover faults of high voltage circuit breaker in time and improve operational reliability.With TMS320C28346 at the core of control，the design used AD7606 to acquire the real-time data of coil current and used envelope average method to filter noise and interference signals.It then extracted the key time and maximum value of the current waves，and calculated the closing-opening time and the maximum value to judge whether the breaker was at fault while the online-monitoring results were displayed on LCD screen.Experimental results show that the design is capable of real-time monitoring of closing/opening coil current of high voltage circuit breaker and of timely detection of failure and achieves good real-time performance and accuracy，with the maximum error of closing coil current standing around 1.29%，closing time error around 0.26%，the maximum error of opening coil current around 1.61% and opening time error around 0.7%.

high voltage circuit breaker;coil current;online monitoring;TMS320C28346;envelope average method

2016-07-06

2016-08-11

福建省教育廳省屬高校專項(xiàng)課題(JK2013036)

陳志英(1978-)，女，副教授，碩士，研究方向?yàn)殡娖髦悄芑夹g(shù)及應(yīng)用.E-mail：chzy207@163.com

TM561

A

1673-4432(2016)05-0006-07