發(fā)電機(jī)非線性狀態(tài)PI勵(lì)磁汽門協(xié)調(diào)控制器的研究

1、82001年5月繼電器RELAY 第29卷 第5期發(fā)電機(jī)非線性狀態(tài)PI勵(lì)磁汽門協(xié)調(diào)控制器的研究康忠健,陳學(xué)允(哈爾濱工業(yè)大學(xué)電氣工程系,黑龍江哈爾濱 150001)摘要:針對(duì)電力系統(tǒng)的強(qiáng)非線性和不確定性的特點(diǎn),直接對(duì)其非線性不確定對(duì)象設(shè)計(jì)了非線性狀態(tài)PI勵(lì)磁汽門協(xié)調(diào)控制器。避免了一般反饋線性化勵(lì)磁汽門綜合控制由于數(shù)學(xué)模型的誤差而影響控制器性能的缺點(diǎn)。所得的控制規(guī)律與系統(tǒng)運(yùn)行點(diǎn)和網(wǎng)絡(luò)結(jié)構(gòu)完全無(wú)關(guān),同時(shí)控制器本身結(jié)構(gòu)簡(jiǎn)單,僅需發(fā)電機(jī)局部信息,易于在多機(jī)系統(tǒng)中實(shí)現(xiàn)分散協(xié)調(diào)控制。仿真計(jì)算表明非線性狀態(tài)PI勵(lì)磁汽門協(xié)調(diào)控制器能夠有效地提高系統(tǒng)的暫態(tài)穩(wěn)定性。關(guān)鍵詞:電力系統(tǒng); 發(fā)電機(jī)協(xié)調(diào)控制; 非線性狀態(tài)

2、PI控制; 魯棒控制中圖分類號(hào):TM76 文獻(xiàn)標(biāo)識(shí)碼:A 文章編號(hào):1003 4897(2001)05 0008 051 引言發(fā)電機(jī)的勵(lì)磁控制和汽門控制是提高電力系統(tǒng)穩(wěn)定性的重要措施,只有發(fā)電機(jī)的勵(lì)磁控制而沒(méi)有良好的汽門控制器則不能從根本上提高電力系統(tǒng)的暫態(tài)穩(wěn)定性;相反,只有發(fā)電機(jī)的汽門控制而沒(méi)有性能良好的勵(lì)磁控制器的共同作用,也不能達(dá)到理想的控制效果。同步發(fā)電機(jī)的勵(lì)磁和汽門的協(xié)調(diào)控制一直是電力系統(tǒng)控制工作者研究的一個(gè)熱點(diǎn)問(wèn)題。1 2文針對(duì)近似線性系統(tǒng)模型將自適應(yīng)控制理論用于電力系統(tǒng)多機(jī)系統(tǒng)勵(lì)磁和汽門的協(xié)調(diào)控制,文3獻(xiàn)用擴(kuò)展直接反饋線性化的方法,綜合考慮了勵(lì)磁系統(tǒng)和調(diào)速系統(tǒng)的作用,設(shè)計(jì)了發(fā)電機(jī)

3、非線性綜4合控制器。文在直接反饋線性化的基礎(chǔ)上,設(shè)計(jì)好勵(lì)磁控制器后考慮勵(lì)磁控制器的調(diào)節(jié)作用設(shè)計(jì)汽5門控制器,文在直接反饋線性化的基礎(chǔ)上運(yùn)用變結(jié)構(gòu)理論設(shè)計(jì)了發(fā)電機(jī)的非線性變結(jié)構(gòu)綜合控制器。基于反饋線性化理論的發(fā)電機(jī)非線性綜合控制解決了發(fā)電機(jī)功率特性中的強(qiáng)非線性問(wèn)題,但由于模型參數(shù)和反饋?zhàn)兞空`差等不確定因素的存在,反饋線性化后的布魯諾夫斯基標(biāo)準(zhǔn)型中也存在一定的誤差,嚴(yán)重時(shí)還有可能形成一個(gè)不穩(wěn)定系統(tǒng)。本文在 參數(shù)動(dòng)態(tài)確定法 設(shè)計(jì)出的非線性擴(kuò)張5 6狀態(tài)觀測(cè)器基礎(chǔ)上,將自抗擾控制技術(shù)引入狀態(tài)7PI調(diào)節(jié)器的設(shè)計(jì),設(shè)計(jì)出適用于多輸入多輸出非線性系統(tǒng)的非線性狀態(tài)PI解耦控制器。并直接對(duì)電力系統(tǒng)非線性不確定模

4、型進(jìn)行設(shè)計(jì),得到同步發(fā)電機(jī)的勵(lì)磁和汽門的協(xié)調(diào)控制器的設(shè)計(jì),以解決電力系統(tǒng)的運(yùn)行點(diǎn)和網(wǎng)絡(luò)結(jié)構(gòu)的非線性不確定問(wèn)題,進(jìn)一步提高勵(lì)磁汽門綜合控制器的控制品質(zhì),改善2 非線性狀態(tài)PI解耦控制器的設(shè)計(jì)2.1 非線性擴(kuò)張狀態(tài)觀測(cè)器對(duì)于受未知外擾作用的不確定非線性系統(tǒng)x(n)(t)=fx(t)!,x(n-1)式中:fx(t)!,x(1)(t)為未知函數(shù),w(t)為未(t)+w(t)+bu(t)(n-1)(n-1)知外擾,u(t)為控制量,x(t),!,x(t)是系統(tǒng)(1)的狀態(tài)變量,x(t)為輸出。對(duì)系統(tǒng)(1)可以采用 參數(shù)動(dòng)態(tài)確定法 設(shè)計(jì)出非線性擴(kuò)張狀態(tài)觀測(cè)器,其描述為:k1z1=z2- g(z1-x(t)

5、g#(z1-x(t)k2z2=z3- g(z1-x(t)g#(z1-x(t)(2)knzn=zn+1- g(z1-x(t)+bu(t)g#(z1-x(t)kn+1zn+1=-g#(z-x(t)g(z1-x(t)其中,a|z|sign(z) |z| 0g(z)=fal(z, , 0)=z,|z|< 0dg(z)g#(z)=dz(n-1)由于式(2)可以z1,z2,!,zn和x(t),!,x(t)的誤差狀態(tài)方程進(jìn)行動(dòng)態(tài)補(bǔ)償,使其成為一個(gè)線(3)x1 x2 xn xn+=-k1-k2 -kn-kn+110 0001 00! !00 10x1 x2 xn xn+0 0-a#(t)(發(fā)電機(jī)非線性狀

6、態(tài)PI勵(lì)磁汽門協(xié)調(diào)控制器的研究9式中,a#(t)是a(t)的導(dǎo)數(shù), x1=z1-x1, x2=z2-x2,!, xn=zn-xn, xn+1=zn+1-a(t)。-k1-k2 -kn-kn+110 0001 00! !00為補(bǔ)償矩陣。則10有 自抗擾 控制、狀態(tài)反饋?zhàn)顑?yōu)控制和常規(guī)PID的優(yōu)點(diǎn)。2.3 非線性狀態(tài)PI解耦控制不失一般性,設(shè)有一雙輸入雙輸出不確定非線性系統(tǒng),其表達(dá)式為:x1 x2(n)A=(t)=f1x1(t)!,x1(n-1)(t),x2(t)!,x2(n-1)(t)+k1,k2,!,kn+1的取值必須保證補(bǔ)償矩陣的特征根全部落在復(fù)平面的左半平面上且其實(shí)部充分的負(fù)。式(2)中的g

7、(z)滿足條件g(0)=0且g#(z)%0。研究表明:選取不同的g(z)就可以得到不同形式的非線性擴(kuò)張狀態(tài)觀測(cè)器,擴(kuò)張狀態(tài)觀測(cè)器的性能只與補(bǔ)償矩陣的極點(diǎn)配置有關(guān)。(n-1)令a(t)=fx(t)!,x(t)+w(t),只要k1,k2,!,kn+1的取值合理,z1,z2,!,zn就能很好地跟蹤系統(tǒng)的狀態(tài)變量x(t),!,x能實(shí)時(shí)估計(jì)出a(t)的值,即zn+1(t)&a(t)(n-1)w1(t)+b11u1(t)+b12u2(t)(n)(t)=f2x1(t)!,x1(n-1)(t),x2(t)!,x2(n-1)(t)+w2(t)+b21u1(t)+b22u2(t)(7)式中,f1、f2分別

8、為通道1和通道2的不確定函數(shù),w1(t)、w2(t)為未知外擾,u1(t)、u2(t)是系統(tǒng)的兩個(gè)輸入。利用擴(kuò)張狀態(tài)觀測(cè)器能對(duì)系統(tǒng)不確定因素和外擾實(shí)時(shí)估計(jì)的特點(diǎn),將通道2的輸入對(duì)通道1的輸出的影響視為通道1的外擾,將通道1的輸入對(duì)通道2的輸出的影響視為通道2的外擾。分別設(shè)計(jì)通道1和通道2的非線性自抗擾狀態(tài)PI控制器,在控制過(guò)程中各自對(duì)其進(jìn)行動(dòng)態(tài)補(bǔ)償,從而實(shí)現(xiàn)系統(tǒng)的解耦控制。令式(6)中b分別為b11、b22,擴(kuò)張狀態(tài)觀測(cè)器中u(t)為通道1或通道2的輸入,那么非線性狀態(tài)PI控制器將把a(bǔ)1(t)=f1x1(t)!,x1(n-1)(t),zn+1就(4)2.2 非線性狀態(tài)PI控制器令u(t)=u0

9、(t)-zn+1(t),將式(4)其帶b入式(1)可得:(n)x(t)=u0(t)(5)式(5)表示的系統(tǒng)是一個(gè)單輸入單輸出線性系7統(tǒng),用常規(guī)狀態(tài)PI控制器設(shè)計(jì)方法對(duì)其進(jìn)行設(shè)計(jì)可以得到其控制方程如式(6)所示。k1z1=z2- g(z1-x1)g#(z1-x1)k2z2=z3- g(z1-x1)g#(z1-x1)knzn=zn+1-g#(z-x)g(z1-x1)+bu(t)11kn+1zn+1=- g(z1-x1)g#(z1-x1)u(t)=u0(t)-zn+1 bTu0(t)=(DZ-v0)+T(z1-v0)dti(6)T式中,v0是輸入,Z=z1z2!zn。擴(kuò)張狀態(tài)觀(t),x2(t)!,

10、x2(n-1)(t)+w1(t)+b12u2(t)(n-1)(n-1)和a2(t)=f2x1(t)!,x1(t),x2(t)!,x2(t)+w2(t)+b21u1(t)看成各自的外擾而進(jìn)行補(bǔ)償,從而實(shí)現(xiàn)解耦控制。兩個(gè)非線性狀態(tài)PI控制器可以單獨(dú)進(jìn)行設(shè)計(jì)。3 非線性狀態(tài)PI勵(lì)磁汽門協(xié)調(diào)控制器3.1 數(shù)學(xué)模型單機(jī)無(wú)窮大系統(tǒng)線路圖如圖1所示,其勵(lì)磁和汽門協(xié)調(diào)控制的數(shù)學(xué)模型為:=!0!=(Pm-Pe)-!HH1E #q=(-Eq+uf)T#d0P m=-PmCPm0+TTT(9)測(cè)器的作用是提供控制器所需的系統(tǒng)的n個(gè)狀態(tài)變量和系統(tǒng)不確定因素及外擾的估計(jì)值。由于在調(diào)節(jié)過(guò)程中,對(duì)系統(tǒng)的不確定因素及外擾進(jìn)行

11、了實(shí)時(shí)補(bǔ)償,從而使控制器具有很強(qiáng)的適應(yīng)性和魯棒性。由于常規(guī)狀態(tài)PI控制器具有二次型指標(biāo)最優(yōu)和具有常規(guī),(8)Pe=E#qIq+(xq-x#d)IdIq(,102001年5月繼電器RELAY 第29卷 第5期定系統(tǒng),運(yùn)用非線性狀態(tài)PI解耦控制對(duì)其進(jìn)行設(shè)計(jì)可以得到勵(lì)磁和汽門的協(xié)調(diào)控制規(guī)律為:uf=u=11zdt(d11z11+d12z12+d13z13)+Ti2111-z14 b1!為發(fā)電機(jī)轉(zhuǎn)子角速度增量,取標(biāo)幺值;0為系統(tǒng)同步角速度;D為發(fā)電機(jī)阻尼系數(shù);H為發(fā)電機(jī)組轉(zhuǎn)子的慣性時(shí)間常數(shù),單位s;u為汽門開(kāi)度控制;Pm為調(diào)節(jié)系統(tǒng)對(duì)應(yīng)的機(jī)械功率(只調(diào)節(jié)高壓缸時(shí),Pm=PH,只調(diào)節(jié)中、低壓缸時(shí),Pm=P

12、ML);Pe為發(fā)電機(jī)輸出的電磁功率;T為等效時(shí)間常數(shù)(含蒸汽容積時(shí)間常數(shù)和油動(dòng)機(jī)時(shí)間常數(shù));C為功率系數(shù)。xd為發(fā)電機(jī)d軸電抗;x#d為發(fā)電機(jī)d軸暫態(tài)電抗;xq為發(fā)電機(jī)的q軸電抗;Id、Iq分別為發(fā)電機(jī)定子電流的d軸和q軸分量,以上各量均取系統(tǒng)的標(biāo)幺值。11zdt-z24 b2(d21z21+d22z22+d23z23)+Ti222 2(12)取: 1=0.014, Ti1=71.2184,D1=d11 d12 d13=1 0.8059 0.0136。 2=-0.0011, Ti2=400.5911,D2=d21 d22 d23=1 0.5907 0.162擴(kuò)張狀態(tài)觀測(cè)器為:100zi1=zi

13、2- g(zi1-zi)g#(zi1-zi)TT圖1 單機(jī)(無(wú)窮大系統(tǒng)zi2=zi3-3750g(zi1-zi)g#(zi1-zi)3.2 非線性狀態(tài)PI勵(lì)磁汽門協(xié)調(diào)控制器的設(shè)計(jì)勵(lì)磁與汽門的協(xié)調(diào)控制策略是:在暫態(tài)過(guò)程中,汽門控制目標(biāo)始終是控制發(fā)電機(jī)的頻率;而勵(lì)磁控制在暫態(tài)過(guò)程開(kāi)始時(shí)的控制目標(biāo)是控制系統(tǒng)的功角,不考慮發(fā)電機(jī)機(jī)端電壓調(diào)節(jié);在暫態(tài)過(guò)程后期勵(lì)磁控制目標(biāo)需要同時(shí)考慮發(fā)電機(jī)的功角和電壓調(diào)節(jié)特性,以實(shí)現(xiàn)勵(lì)磁與汽門的協(xié)調(diào)控制。在暫態(tài)過(guò)程開(kāi)始時(shí)勵(lì)磁控制不考慮電機(jī)機(jī)端電壓調(diào)節(jié)特性,勵(lì)磁控制目標(biāo)是! =0,汽門控制目標(biāo)是!=0。令z1=! 0,z2=!,則式(9)可轉(zhuǎn)化為一個(gè)2輸入2輸出非線性系統(tǒng):

14、!P z1=m-Pe-)z1HH(10)!P)z2=m-Pe-)z2HH由于P e=E Iq+(xq-xd#)q#Iq+Eq# IdIq+Id Iq, P)m=故z1=!62500zi3=zi4- g(zi1-zi)-biufg#(zi1-zi)zi4=-390625g(zi1-zi)g#(zi1-zi)i=1,2(13)其中,zi1、zi2、zi3和zi4是式(13)表示的擴(kuò)張狀態(tài)觀測(cè)器的輸出,補(bǔ)償矩陣的極點(diǎn)配置為P1,2,3,4=-25。非線性函數(shù)g(z)為:|z|sign(z) |z|g(z)=fal(z, , )= |z|<令 a1(t)=PmCPm0xq-xd#Eq-+-( I

15、dIq+THTHHTd0#HIqIqE#qDCIdI q)-)z1- Iq-(+b1)uf+uHHTHHTd0#PmCPm+,-TTTTPmCPm0xq-xd#Eq-+-( IdIq+THTHHTd0#HIqE#q)z1- Iq-uf+uHHTHHTd0#IqId Iq)-!2P)ePmCPm a2(t)=-H)z-H+-+HT-HTHT,則z14、z24可以動(dòng)態(tài)地跟蹤a1(t)、a2(t)的HT值。在暫態(tài)過(guò)程后期時(shí)勵(lì)磁控制需考慮發(fā)電機(jī)端電壓特性,令式(13)中的z1為:z1=k!dt+kv!Vt(14)式中,!Vt=Vt0-Vt,Vt、Vt0分別是發(fā)電機(jī)機(jī)端電壓及其設(shè)定值,k、kv是發(fā)電機(jī)轉(zhuǎn)

16、子轉(zhuǎn)速偏差和機(jī)端電壓偏差的加權(quán)系數(shù)。通過(guò)數(shù)值仿真k和kv的值確定為:當(dāng)暫態(tài)過(guò)程開(kāi)始時(shí),k=7,kv=0;當(dāng)進(jìn)入暫k7,。P)PmCPmez=-z2-+-+uHHHTHTHTHT(11)令b1=,b2=,其中Iq0為故障前IqHTHTd0#(-Iq0發(fā)電機(jī)非線性狀態(tài)PI勵(lì)磁汽門協(xié)調(diào)控制器的研究11由系統(tǒng)的運(yùn)行點(diǎn)變化、網(wǎng)絡(luò)結(jié)構(gòu)變化、勵(lì)磁控制和汽門控制的相互影響等不確定因素引起的變化都可以通過(guò)擴(kuò)張狀態(tài)觀測(cè)器的輸出z14和z24估計(jì)出來(lái),并在式(12)表示的勵(lì)磁和汽門協(xié)調(diào)控制律中分別給予動(dòng)態(tài)補(bǔ)償,因此,對(duì)系統(tǒng)運(yùn)行點(diǎn)和網(wǎng)絡(luò)結(jié)構(gòu)的變化具有極強(qiáng)的適應(yīng)性和魯棒性,同時(shí)也實(shí)現(xiàn)了勵(lì)磁和汽門的協(xié)調(diào)控制。同時(shí),勵(lì)磁汽

17、門協(xié)調(diào)控制律只需反饋所調(diào)節(jié)發(fā)電機(jī)的轉(zhuǎn)子轉(zhuǎn)速偏差信號(hào)和機(jī)端電壓偏差信號(hào),易于在多機(jī)系統(tǒng)中實(shí)現(xiàn)分散協(xié)調(diào)控制,由于減少了反饋量,從而提高了控制器本身的可靠性。非線性狀態(tài)PI勵(lì)磁汽門協(xié)調(diào)控制器的結(jié)構(gòu)如圖2所示:圖3 仿真條件1下的系統(tǒng)響應(yīng)圖2 非線性狀態(tài)PI勵(lì)磁汽門協(xié)調(diào)控制器框圖圖4 仿真條件2下的系統(tǒng)響應(yīng)4 數(shù)值仿真及分析單機(jī)無(wú)窮大系統(tǒng)如圖1,系統(tǒng)參數(shù)為:Xd=0.1,X#d=0.031,Xq=0.1437,H=42,Pe=7.0,Vt=1.0,XL1=XL2=0.1705,XT=0.034(以上各量均為標(biāo)幺值),Td0=10.2s。仿真列出了非線性狀態(tài)PI、非線性PID和常規(guī)控制三種勵(lì)磁和汽門協(xié)調(diào)

18、控制器的結(jié)果。曲線1、曲線2和曲線3分別表示非線性狀態(tài)PI、非線性PID和常規(guī)控制三種勵(lì)磁汽門協(xié)調(diào)控制器的響應(yīng)曲線。仿真條件為:(1)系統(tǒng)于t=0.1s在F點(diǎn)發(fā)生三相短路故障,t=0.25s故障消除,仿真結(jié)果如圖3。(2)系統(tǒng)于t=0.1s在F點(diǎn)發(fā)生三相短路故障,t=0.31s故障消除,仿真結(jié)果如圖4。(3)條件同(1),發(fā)電機(jī)運(yùn)行點(diǎn)由 0=67.3變?yōu)?0=18.3,仿真結(jié)果如圖5。從圖3、圖4和圖5中可以看出,非線性狀態(tài)PI勵(lì)磁汽門協(xié)調(diào)控制器比非線性PID勵(lì)磁汽門協(xié)調(diào)控,圖5 仿真條件3下的系統(tǒng)響應(yīng)線性狀態(tài)PI勵(lì)磁汽門協(xié)調(diào)控制器比非線性PID勵(lì)磁汽門綜合控制器對(duì)系統(tǒng)運(yùn)行點(diǎn)的變化具有良好的魯

19、棒性。5 結(jié)論非線性狀態(tài)PI勵(lì)磁汽門協(xié)調(diào)控制器直接對(duì)電力系統(tǒng)不確定對(duì)象進(jìn)行設(shè)計(jì),避免了一般反饋線性化勵(lì)磁汽門協(xié)調(diào)控制由于數(shù)學(xué)模型的誤差而影響控制器性能的缺點(diǎn),因而具有更強(qiáng)的強(qiáng)適應(yīng)性和強(qiáng)魯,122001年5月繼電器RELAY 第29卷 第5期Transientstabilityenhancementandvoltageregulationofpow ersystem.IEEETransactiononPowerSystem,1993,8(2):620 626.5 韓京清.非線性狀態(tài)誤差反饋控制律(NLSEF.控制與決策,1995,10(3):221 225.6 韓京清.自抗擾控制器及其應(yīng)用.控制

20、與決策,1998,13(1):19 23.7 宋文忠,蔣國(guó)平.一種新型狀態(tài)PI調(diào)節(jié)器的研究.自動(dòng)化學(xué)報(bào),1997,23(6):786 792.收稿日期: 2000 10 18作者簡(jiǎn)介: 康忠健(1971-),男,博士生,主要研究方向?yàn)殡娏ο到y(tǒng)穩(wěn)定分析及非線性控制、微機(jī)自動(dòng)控制; 陳學(xué)允(1934-),男,博士生導(dǎo)師,主要研究方向?yàn)殡娏ο到y(tǒng)穩(wěn)定分析及控制。結(jié)構(gòu)無(wú)關(guān),控制器本身結(jié)構(gòu)簡(jiǎn)單,僅需發(fā)電機(jī)局部信息,易于在多機(jī)系統(tǒng)中實(shí)現(xiàn)分散協(xié)調(diào)控制。數(shù)值仿真結(jié)果表明:非線性狀態(tài)PI汽門控制器能有效地提高電力系統(tǒng)的暫態(tài)穩(wěn)定特性。參考文獻(xiàn):1 劉偉,余貽鑫,田樹(shù)苞.變結(jié)構(gòu)模型參考自適應(yīng)勵(lì)磁與調(diào)速綜合控制器(上)

21、.電力系統(tǒng)自動(dòng)化,1998,22(4):8 12.2 PahalawaththaNC,HopeGS,MalikOP.Multivariableself tuningpowersystemstabilizersimulationandimplement atioinstudies.IEEETransactiononEnergyConversion,1991,6(2):310 318.3 WangYouyi,DavidJHill.Robustnonlinearcoordinatedcon trolofpowersystems.Automatic,1996,32(4):611 618.4 WangY

22、ouyi,DavidJHill,RichardHMiddleton,LongGao.TheresearchonnonlinearstatePIcoordinatedcontrollerofgenerator+sexcitationandvalveKANGZhong jian,CHENXue yun(HarbinInstituteofTechnology,Harbin150001,China)Abstract: Asforthestrongnonlinearanduncertaincharacteristicofthepowersystem,thenonlinearstatePIcoordina

23、tedcontrollerofexcitationandvalveisproposedinthispaper.Thecontrollerisdesignedfortheuncertainmodelobjectandtheshortcomingofgeneralfeedbacklin earizationcoordinatedcontrolofexcitationandvalve,whichitsperformanceisaffectedbythemathematicalmodelerror,isavoided.Thecontrollawiscompletelyindependentontheo

24、perationconditionandthenetstructureofthepowersystem.Thecontrollerstructureissimpleandonlylocalinformationofgeneratorisneeded.Thedecentralizedcontrollercanbeaccomplishedeasilyinmultimachinepowersystem.ThesimulationshowsthatthenonlinearstatePIcoordinatedcontrollerofexcitationandvalveissatisfiedforimprovementofpowersystem.Keywords: powersystem; coordinatedcontrollerofexcitationandvalve; nonlinearstatePIcontrol; robustcontrol(上接第7頁(yè))7 盧繼平,葉一麟.一種用于平行雙回線的故障點(diǎn)精確定位方法.重慶大學(xué)學(xué)報(bào),1988,(4).8 陳臘生,孫淑信.采用實(shí)數(shù)對(duì)稱分量方法準(zhǔn)確探測(cè)平行雙回輸電線路故障距離.全國(guó)高校電自專業(yè)第四屆年會(huì)論文集,1988.收稿日期: 2000 10 16作者簡(jiǎn)介: 粟小華(1961-),男,工學(xué)碩士,工程師,從事高