專利名稱：變速恒頻風(fēng)力發(fā)電機(jī)動態(tài)最優(yōu)能量捕獲方法
技術(shù)領(lǐng)域：
本發(fā)明涉及一種變速恒頻風(fēng)力發(fā)電機(jī)動態(tài)最優(yōu)能量捕獲方法。主要應(yīng) 用于現(xiàn)代水平軸升力型變速恒頻風(fēng)力發(fā)電機(jī)，可以提高風(fēng)力發(fā)電系統(tǒng)在低 于額定風(fēng)速情況下的從湍流風(fēng)中捕獲的能量，增加風(fēng)力發(fā)電機(jī)的年發(fā)電量。
背景技術(shù)：
現(xiàn)代風(fēng)力發(fā)電機(jī)一般為水平軸升力型，也就是說利用發(fā)電機(jī)葉片的特 殊設(shè)計的氣動外形，在風(fēng)速和風(fēng)輪轉(zhuǎn)速共同作用下在葉片的每一段上產(chǎn)生 升力，在升力的作用下，產(chǎn)生旋轉(zhuǎn)力矩，帶動發(fā)電機(jī)發(fā)出電能。相對于阻 力型風(fēng)力發(fā)電機(jī)，現(xiàn)代水平軸升力型風(fēng)力發(fā)電機(jī)具有效率高的特點(diǎn)。為了 能夠在不同的風(fēng)速下，獲得最大的能量捕獲，現(xiàn)代大型水平軸升力型風(fēng)力 發(fā)電機(jī)， 一般采用變速恒頻方式操作風(fēng)力發(fā)電機(jī)。因此，會在不同的風(fēng)速 下采用不同的控制策略。根據(jù)參考文獻(xiàn)l,在風(fēng)速低于額定風(fēng)速下工作時， 風(fēng)力發(fā)電機(jī)主要是工作在欠功率狀態(tài)，即葉片的槳距角被調(diào)整到最佳槳距角位置，調(diào)節(jié)風(fēng)輪轉(zhuǎn)速，達(dá)到較高的能量捕獲；在風(fēng)速高于額定風(fēng)速情況 下，保持轉(zhuǎn)速不變，調(diào)節(jié)槳距角，減少風(fēng)輪能量捕獲，進(jìn)入恒功率工作狀 態(tài)。圖1是IEC標(biāo)準(zhǔn)中的三類風(fēng)場的風(fēng)速分布，我國大部分風(fēng)場介于二類風(fēng)場和三類風(fēng)場之間，具有很強(qiáng)的代表性。以1.5兆瓦變速恒頻風(fēng)力發(fā)電 機(jī)為例，風(fēng)輪直徑可以達(dá)到70-77米，額定風(fēng)速在12米/秒左右。風(fēng) 力發(fā)電機(jī)主要年工作小時數(shù)，從切入風(fēng)速(4m/s左右)到額定風(fēng)速(12m/s 左右)為大約6000小時，而額定風(fēng)速以上的年工作小時數(shù)大約為1381小 時。這樣，提高在風(fēng)力發(fā)電機(jī)工作在低于額定風(fēng)速情況下的發(fā)電效率，將 對于提高風(fēng)力發(fā)電機(jī)年發(fā)電量產(chǎn)生巨大的影響。對于兆瓦級變速恒頻風(fēng)力發(fā)電機(jī)來說，為了能夠優(yōu)化地捕獲風(fēng)中的能 量，主要依靠通過靜態(tài)仿真得出的優(yōu)化轉(zhuǎn)矩轉(zhuǎn)速曲線進(jìn)行開環(huán)控制。如圖 2所示的黑色曲線(摘自參考文獻(xiàn)l)。而一般的控制系統(tǒng)框圖如圖3所示。 根據(jù)靜態(tài)曲線控制有以下問題.-1、 靜態(tài)曲線考慮的是正面吹來的風(fēng)作用。而實際上，由于風(fēng)向時刻改 變，而發(fā)電機(jī)的偏航控制反應(yīng)較慢，導(dǎo)致偏航誤差，風(fēng)會在左右方向和上 下方向與風(fēng)輪所成的夾角，使風(fēng)力發(fā)電機(jī)處于非最優(yōu)能量捕獲狀態(tài)，降低 了能量捕獲。2、 由于靜態(tài)曲線根據(jù)穩(wěn)態(tài)風(fēng)得出，而實際的風(fēng)場為三維湍流風(fēng)場。這 樣，就造成風(fēng)輪獲得的升力與靜態(tài)情況有所不同，導(dǎo)致能量捕獲下降。3、 靜態(tài)曲線一般根據(jù)葉片數(shù)據(jù)，有常見的載荷仿真軟件，如GHBladed 和Flex5，根據(jù)葉片數(shù)據(jù)直接計算出來，由于實際葉片和葉片數(shù)據(jù)模型之 間存在誤差，經(jīng)過計算得出的曲線本身就存在誤差，使用這個曲線控制時， 也會降低能量捕獲。發(fā)明內(nèi)容本發(fā)明的目的是針對上述現(xiàn)有技術(shù)的不足，克服上述恒頻變速風(fēng)力發(fā) 電機(jī)在低于額定風(fēng)速下的靜態(tài)非最優(yōu)能量捕獲的問題，提供一種變速恒頻 風(fēng)力發(fā)電機(jī)動態(tài)最優(yōu)能量捕獲方法。本發(fā)明的目的是這樣實現(xiàn)的 一種變速恒頻風(fēng)力發(fā)電機(jī)動態(tài)最優(yōu)能量捕 獲方法，其特征在于所述方法包括風(fēng)輪氣動轉(zhuǎn)矩測量計算方法和動態(tài)最優(yōu) 能量捕獲調(diào)節(jié)方法兩部分，所述的風(fēng)輪氣動轉(zhuǎn)矩測量計算方法是在低速軸兩端安裝扭轉(zhuǎn)角度測 量儀，該測量儀測量低速軸在工作狀態(tài)下，實時測量出低速軸的扭轉(zhuǎn)角度， 扭轉(zhuǎn)角度以弧度為單位，低速軸在鍛造的時候通過做實驗的方法將低速軸 的扭轉(zhuǎn)剛度和扭轉(zhuǎn)阻尼系數(shù)測量出來，根據(jù)低速軸的這兩個參數(shù)以及實時 測量的低速軸的扭轉(zhuǎn)角度，可以實時測量計算出低速軸的扭轉(zhuǎn)力矩，根據(jù) 低速軸的扭轉(zhuǎn)力矩和風(fēng)輪轉(zhuǎn)速的變化計算出瞬時風(fēng)輪的氣動轉(zhuǎn)矩，所述動態(tài)最優(yōu)能量捕獲調(diào)節(jié)方法是根據(jù)測量計算出的瞬時風(fēng)輪的氣動 轉(zhuǎn)矩以及功率的變化梯度計算出最優(yōu)能量捕獲轉(zhuǎn)矩，施加到發(fā)電機(jī)軸上， 該發(fā)電機(jī)軸上轉(zhuǎn)矩再通過齒輪箱施加到低速軸，從而調(diào)節(jié)風(fēng)輪轉(zhuǎn)速，使風(fēng) 輪達(dá)到最優(yōu)能量捕獲狀態(tài)，達(dá)到動態(tài)最優(yōu)能量捕獲的目的，具體計算過程 如下，所述瞬時風(fēng)輪的氣動轉(zhuǎn)矩的計算公式如下-r"=力ty二+r =+尺,e+w ， 其中，r。為風(fēng)輪獲得的氣動轉(zhuǎn)矩，r為低速軸與齒輪箱鏈接處的扭矩，單 位為牛頓米，,為風(fēng)輪的轉(zhuǎn)動慣量，單位是公斤米平方，w;為風(fēng)輪轉(zhuǎn)速的導(dǎo)數(shù)，6為低速軸兩端扭轉(zhuǎn)角度，單位是弧度，^為低速軸的扭轉(zhuǎn)角度的導(dǎo) 數(shù)，《,是低速軸的扭轉(zhuǎn)剛度，單位是牛頓米/弧度，5、為低速軸的扭轉(zhuǎn)阻 尼系數(shù)，單位是牛頓米/(弧度/秒)。所述動態(tài)最優(yōu)能量捕獲調(diào)節(jié)方法的控制公式為7; (/t+1)=厶("—-1)) / r、.+^+s, c—維—1)) / r、 + 7(《w - d i)) / ;其中，7;為發(fā)電機(jī)的電轉(zhuǎn)矩，單位是牛頓米，(t+i)， (w，和(/t-i)分 別表示下一時刻，當(dāng)前時刻，和前一時刻對應(yīng)變量的值，令7 = /《，《為 測量的電功率。"作為歸一化的調(diào)節(jié)參數(shù)，"的調(diào)節(jié)范圍為0.1-0.3。


圖1是以往IEC標(biāo)準(zhǔn)中的三類風(fēng)場的風(fēng)速分布圖。圖2是優(yōu)化轉(zhuǎn)矩轉(zhuǎn)速曲線圖。圖3是通常的控制系統(tǒng)框圖。圖4是本發(fā)明變速恒頻風(fēng)力發(fā)電機(jī)動態(tài)最優(yōu)能量捕獲方法的方案框圖。圖5是本發(fā)明的風(fēng)力發(fā)電機(jī)簡化模型圖。
具體實施方式
本發(fā)明所涉及的變速恒頻風(fēng)力發(fā)電機(jī)動態(tài)最優(yōu)能量捕獲方法具體的方案框圖4所示，方案由兩部分組成風(fēng)輪氣動轉(zhuǎn)矩測量計算方法和動態(tài)最 優(yōu)能量捕獲調(diào)節(jié)方法。所述風(fēng)輪氣動轉(zhuǎn)矩測量計算所述的風(fēng)輪氣動轉(zhuǎn)矩測量計算方法包括測 量和計算兩部分。第一部分、風(fēng)輪氣動轉(zhuǎn)矩測量風(fēng)輪氣動轉(zhuǎn)矩測量是通過設(shè)置在低速軸兩端的扭轉(zhuǎn)角度測量儀，該測 量儀測量低速軸在工作狀態(tài)下，實時直接測量出低速軸的扭轉(zhuǎn)角度入以 弧度為單位。再根據(jù)低速軸的扭轉(zhuǎn)剛度和扭轉(zhuǎn)阻尼系數(shù)(低速軸在鍛造的 時候通過做實驗的方法將低速軸的扭轉(zhuǎn)剛度和扭轉(zhuǎn)阻尼系數(shù)測量出來)，根 據(jù)低速軸的這兩個參數(shù)以及實時測量的低速軸的扭轉(zhuǎn)角度，可以實時測量 出低速軸的扭轉(zhuǎn)力矩。根據(jù)低速軸的扭轉(zhuǎn)力矩和風(fēng)輪角速度的變化計算出 瞬時風(fēng)輪的氣動轉(zhuǎn)矩。風(fēng)輪的氣動轉(zhuǎn)矩直接跟風(fēng)的湍流密度，偏航誤差， 以及風(fēng)的垂直入射角度有關(guān)。因此，直接測量的風(fēng)輪捕獲的氣動轉(zhuǎn)矩考慮 到了風(fēng)速湍流密度，偏航誤差，以及風(fēng)的垂直入射角度。根據(jù)測算出的風(fēng) 輪捕獲的氣動轉(zhuǎn)矩，疊加功率的梯度信息，設(shè)計了動態(tài)優(yōu)化跟蹤調(diào)節(jié)算法。風(fēng)輪氣動轉(zhuǎn)矩計算風(fēng)輪捕獲的氣動轉(zhuǎn)矩直接反映了風(fēng)力發(fā)電機(jī)的能量捕獲狀態(tài)，因此如 何獲得風(fēng)輪捕獲的氣動轉(zhuǎn)矩，是最優(yōu)能量捕獲算法的基礎(chǔ)。根據(jù)參考文獻(xiàn) [1]，可知風(fēng)輪捕獲的氣動轉(zhuǎn)矩與風(fēng)速和風(fēng)輪轉(zhuǎn)速以及葉片特性相關(guān)7;=上一3~^2 (i)其中，7;為風(fēng)輪獲得的氣動轉(zhuǎn)矩，單位為牛頓米，p為空氣密度，單 位公斤/立方米，冗為圓周率，/ 為風(fēng)輪半徑，單位為米， 〃為風(fēng)輪平面上 的等效風(fēng)速，單位為米/秒，G為功率捕獲系數(shù)，入為葉尖速比。由于風(fēng)力發(fā)電機(jī)的風(fēng)輪平面很大，例如1.5MW的風(fēng)力發(fā)電機(jī)風(fēng)輪平面面積至少為 3800平方米，因此測量整個風(fēng)輪平面的等效風(fēng)速非常困難，所以這個公式在實際的控制中無法使用。根據(jù)風(fēng)力發(fā)電機(jī)的簡化模型(如圖5)，可以得出如下關(guān)系<formula>formula see original document page 9</formula>其中，r。為風(fēng)輪獲得的氣動轉(zhuǎn)矩，r為低速軸與齒輪箱鏈接處的扭矩，z;為 發(fā)電機(jī)電轉(zhuǎn)矩，轉(zhuǎn)矩的單位為牛頓米，入和^分別為風(fēng)輪和發(fā)電機(jī)的轉(zhuǎn)動 慣量，單位是公斤米平方，A和^分別是風(fēng)輪和電動機(jī)的轉(zhuǎn)速，單位是弧 度/秒，w;和《分別是風(fēng)輪和電動機(jī)的轉(zhuǎn)速的導(dǎo)數(shù)，e為低速軸兩端扭轉(zhuǎn)角 度，單位是弧度，伊為低速軸的扭轉(zhuǎn)角度的導(dǎo)數(shù)，K、是低速軸的扭轉(zhuǎn)剛度， 單位是牛頓米/弧度，A為低速軸的扭轉(zhuǎn)阻尼系數(shù)，單位是牛頓米/(弧度/ 秒)。"為齒輪箱變速比，無單位。有文章采用上述公式作為自適應(yīng)控制的基礎(chǔ)，而在實際控制過程中， 由于齒輪箱內(nèi)部的齒輪嚙合過程中，轉(zhuǎn)速的傳遞會出現(xiàn)非線性，也就是說用公式(4)后半部分公式^((^-丄^)心+^^-丄"」，計算出的r誤差較大，從而導(dǎo)致控制誤差，不能達(dá)到最優(yōu)的能量捕獲的目的。如果能夠精 確測量和計算出風(fēng)輪獲得的氣動轉(zhuǎn)矩，才能更加準(zhǔn)確地控制風(fēng)力發(fā)電機(jī)。 因此，本發(fā)明提出了在低速軸的兩端安裝軸的扭轉(zhuǎn)測量裝置，直接測量出 低速軸的扭轉(zhuǎn)角度《及其變化，根據(jù)公式(4)中K" + S^'，就可以得出精確的低速軸上的扭轉(zhuǎn)力矩，這樣，就可以通過低速軸的扭轉(zhuǎn)力矩和風(fēng)輪速度的改變，計算出瞬時風(fēng)輪從風(fēng)中獲取的氣動轉(zhuǎn)矩<formula>formula see original document page 10</formula>。測速裝置可以參見文獻(xiàn)2。第二部分、動態(tài)最優(yōu)能量捕獲方法 瞬時電功率C7;.^，電功率對時間的導(dǎo)數(shù)為-<formula>formula see original document page 10</formula>根據(jù)參考文獻(xiàn)3，基于梯度的最優(yōu)化方法是一種常用的優(yōu)化方法，本 發(fā)明為了能夠優(yōu)化地獲得最大功率，將梯度法應(yīng)用于風(fēng)力發(fā)電機(jī)動態(tài)最優(yōu) 能量捕獲，即沿梯度(&)方向改變電功率，經(jīng)過多次的迭代就會得到最優(yōu)的電功率捕獲。由于整個傳動系統(tǒng)的轉(zhuǎn)動慣量比較大，速度很難快速改變，所以電功 率的改變直接和電轉(zhuǎn)距的改變成正比# ^.^。因此，電轉(zhuǎn)矩的改變 梯度和功率改變梯度是一致的。電功率的梯度表示了電功率的變化方向， 如果電功率的梯度為正數(shù)，則說明電功率增加，如果電功率的梯度為負(fù)數(shù)， 電功率變小。為了獲得最大功率，就可以按照電功率變化梯度來調(diào)節(jié)電轉(zhuǎn) 距大小，在每步的控制步驟中，按照功率梯度方向調(diào)整功率。也就是說， 為了能夠動態(tài)跟蹤，達(dá)到最優(yōu)能量捕獲，在電氣轉(zhuǎn)矩上疊加一個電功率的 梯度項，艮P，<formula>formula see original document page 10</formula>其中;/為調(diào)節(jié)參數(shù)，是正常數(shù)，其中，^， e和《是可測量的，而《;,， 和，是可以計算的，因此可以直接進(jìn)行控制。公式(6)就是動態(tài)最優(yōu)能量捕獲的控制公式，但是，計算機(jī)控制系統(tǒng)為離散控制系統(tǒng)，所以將^,e'，和^離散化，得到<formula>formula see original document page 11</formula>其中，z;為控制步長，單位為秒。由公式(7)可以推出實際控制公式<formula>formula see original document page 11</formula> 其中，(* + 1)，("，和(^-l)分別表示下一時刻，當(dāng)前時刻，和前一時刻對 應(yīng)變量的值。如果;/越大，對功率變化放大越大，動態(tài)跟蹤越快，太大則 會導(dǎo)致系統(tǒng)不穩(wěn)定；如果; 越小，對功率變化放大越小，動態(tài)跟蹤越慢， 系統(tǒng)越穩(wěn)定。7值與風(fēng)力發(fā)電機(jī)的葉片和當(dāng)前功率有關(guān)系，其中由于變速 恒頻風(fēng)力發(fā)電機(jī)的功率范圍很寬，所以功率的影響更大，也就是說在小功 率時"值應(yīng)該稍大一些，在大功率時"值應(yīng)該稍小一點(diǎn)。因此，令;7 = "/《。 "作為歸一化的調(diào)節(jié)參數(shù)。把公式(8)用軟件實現(xiàn)，作為外部控制器進(jìn)行仿真研究。仿真條件是 在IEC3類風(fēng)場中，將傳統(tǒng)控制器和本發(fā)明實現(xiàn)的控制器進(jìn)行比對試驗，得 出表格l.由仿真得出a的調(diào)節(jié)范圍為0.1-0.3。從表格中的數(shù)據(jù)可以看出 本方法優(yōu)于傳統(tǒng)方法，每年可額外產(chǎn)生12947元的額外收益(風(fēng)電上網(wǎng)電 價按0.5元/千瓦時計算)。參考文獻(xiàn)[1] Tony Burton, David Sharpe， Nick Jenkins, Ervin Bossanyi. 化WJohn Wiley & Sons, Ltd, West Sussex, P019 IUD，England. 2001。[2]喻洪麟，朱傳新，楊張利，光柵扭矩動態(tài)測量系統(tǒng)設(shè)計及實現(xiàn)。應(yīng) 用光學(xué)，2006年05期.[3]袁亞湘，最優(yōu)化理論與方法，科學(xué)出版社。表I能量捐l獲方法對比表平均風(fēng)速 (m/s)4. 005.006. 007. 008. 00傳統(tǒng)方法 捕獲能量 (kW)40. 50106.20208.10348.70535.90本發(fā)明方 法捕獲能 量(kW)41.07107.40211.93358. 89551.21a取值0. 300. 250. 220. 180. 10風(fēng)速分布 時間/年780. 28860. 42885. 82862. 27799. 63多發(fā)電度 數(shù)(跳)產(chǎn)生效益 (元)多發(fā)電能 (跳)443. 381029. 923392. 248788.9312239. 6525894. 1212947. 0權(quán)利要求
1、一種變速恒頻風(fēng)力發(fā)電機(jī)動態(tài)最優(yōu)能量捕獲方法，其特征在于所述方法包括風(fēng)輪氣動轉(zhuǎn)矩測量計算方法和動態(tài)最優(yōu)能量捕獲調(diào)節(jié)方法兩部分，所述的風(fēng)輪氣動轉(zhuǎn)矩測量計算方法是在低速軸兩端安裝扭轉(zhuǎn)角度測量儀，該測量儀測量低速軸在工作狀態(tài)下，實時測量出低速軸的扭轉(zhuǎn)角度，扭轉(zhuǎn)角度以弧度為單位，低速軸在鍛造的時候通過做實驗的方法將低速軸的扭轉(zhuǎn)剛度和扭轉(zhuǎn)阻尼系數(shù)測量出來，根據(jù)低速軸的這兩個參數(shù)以及實時測量的低速軸的扭轉(zhuǎn)角度，可以實時測量計算出低速軸的扭轉(zhuǎn)力矩，根據(jù)低速軸的扭轉(zhuǎn)力矩和風(fēng)輪轉(zhuǎn)速的變化計算出瞬時風(fēng)輪的氣動轉(zhuǎn)矩，所述動態(tài)最優(yōu)能量捕獲調(diào)節(jié)方法是根據(jù)測量計算出的瞬時風(fēng)輪的氣動轉(zhuǎn)矩以及功率的變化梯度計算出最優(yōu)能量捕獲轉(zhuǎn)矩，施加到發(fā)電機(jī)軸上，該發(fā)電機(jī)軸上轉(zhuǎn)矩再通過齒輪箱施加到低速軸，從而調(diào)節(jié)風(fēng)輪轉(zhuǎn)速，使風(fēng)輪達(dá)到最優(yōu)能量捕獲狀態(tài)，具體計算過程如下，所述瞬時風(fēng)輪的氣動轉(zhuǎn)矩的計算公式如下Ta＝Jrω′+T＝Jrωr′+Ksθ+Bsθ′，其中，Ta為風(fēng)輪獲得的氣動轉(zhuǎn)矩，T為低速軸與齒輪箱鏈接處的扭矩，單位為牛頓米，Jr為風(fēng)輪的轉(zhuǎn)動慣量，單位是公斤米平方，ωr′為風(fēng)輪轉(zhuǎn)速的導(dǎo)數(shù)，θ為低速軸兩端扭轉(zhuǎn)角度，單位是弧度，θ′為低速軸的扭轉(zhuǎn)角度的導(dǎo)數(shù)，Ks是低速軸的扭轉(zhuǎn)剛度，單位是牛頓米/弧度，Bs為低速軸的扭轉(zhuǎn)阻尼系數(shù)，單位是牛頓米/(弧度/秒)，所述動態(tài)最優(yōu)能量捕獲調(diào)節(jié)方法的控制公式為Te(k+1)＝Jr(ωr(k)-ωr(k-1))/Ts+Ksθ+Bs(θ(k)-θ(k-1))/Ts+η(Pe(k)-Pe(k-1))/Ts其中，Te為發(fā)電機(jī)的電轉(zhuǎn)矩，單位是牛頓米，(k+1)，(k)，和(k-1)分別表示下一時刻，當(dāng)前時刻，和前一時刻對應(yīng)變量的值，令η＝α/Pe，其中Pe為測量的電功率。α作為歸一化的調(diào)節(jié)參數(shù)，α的調(diào)節(jié)范圍為0.1-0.3。
全文摘要
本發(fā)明涉及一種變速恒頻風(fēng)力發(fā)電機(jī)動態(tài)最優(yōu)能量捕獲方法，包括風(fēng)輪氣動轉(zhuǎn)矩測量計算方法和動態(tài)最優(yōu)能量捕獲調(diào)節(jié)方法兩部分，風(fēng)輪氣動轉(zhuǎn)矩測量計算方法是在低速軸兩端安裝扭轉(zhuǎn)角度測量儀，實時測量出低速軸的扭轉(zhuǎn)角度，將低速軸的扭轉(zhuǎn)剛度和扭轉(zhuǎn)阻尼系數(shù)測量出來，根據(jù)低速軸的這兩個參數(shù)以及實時測量的低速軸的扭轉(zhuǎn)角度，可以實時測量計算出低速軸的扭轉(zhuǎn)力矩，根據(jù)低速軸的扭轉(zhuǎn)力矩和風(fēng)輪轉(zhuǎn)速的變化計算出瞬時風(fēng)輪的氣動轉(zhuǎn)矩，所述動態(tài)最優(yōu)能量捕獲調(diào)節(jié)方法是根據(jù)測量計算出的瞬時風(fēng)輪的氣動轉(zhuǎn)矩以及功率的變化梯度計算出最優(yōu)能量捕獲轉(zhuǎn)矩，施加到發(fā)電機(jī)軸上，該發(fā)電機(jī)軸上轉(zhuǎn)矩再通過齒輪箱施加到低速軸，從而調(diào)節(jié)風(fēng)輪轉(zhuǎn)速，使風(fēng)輪達(dá)到最優(yōu)能量捕獲狀態(tài)。
文檔編號H02P9/04GK101252334SQ20081002043
公開日2008年8月27日 申請日期2008年3月6日 優(yōu)先權(quán)日2008年3月6日
發(fā)明者雷 張 申請人:江陰遠(yuǎn)景能源科技有限公司