维普资讯 http://www.cqvip.com 第12卷第4期 燃烧科学与技术 V01.12 No.4 2006年8月 Journal of Combustion Science and Technology Aug.2006 基于SVM的燃煤电站锅炉飞灰含碳量预测 蔡杰进,马晓茜 (华南理工大学电力学院,广州510640) 摘要:将支持向量机方法引人燃煤电站锅炉飞灰含碳量预测领域.该预测方法很好地建立了燃煤电站锅炉飞灰 含碳量特性与运行参数之间的复杂关系模型,并考虑到运行参数之间的耦合性,具有预测能力强、全局最优及泛化 性好等优点.将该方法应用于某300 MW燃煤电站锅炉中,经过训练后的SVM模型对检验样本飞灰含碳量进行预 测,均方根误差和平均相对误差分别为1.39%和1.30%,相当于BP网络模型的22.20%和21.07%.应用结果表明, 支持向量机方法优于多层BP神经网络法,能很好地满足预测要求. 关键词:锅炉;飞灰含碳量;支持向量机 中图分类号:TK223 文献标志码:A 文章编号:1006.8740(2006)04—0312—06 Forecasting Unburned Carbon Content in the Fly Ash from Coal-Fired Utility Boilers Based on SVM CAI Jie-jin,MA Xiao—qian (Electric Power CoHege,South China University of Technology,Guangzhou 510640,China) Abstract:A new algorithmforforecastingthe unburned carbon contentintheflyashfrom coal—ifred utiliyt boilers based OH hte support vector machine(SVM)method is presented.This forecasting method establishes a model to relfect hte compli— cated relations between the unburned carbon content characteristics in the fly ash and the operating parameters. tll the coupling performance of every parameter considered,It has the advantages of hish forecasting accuracy,global optima prop— erty,and more generalized performance.Apphed to a 300 MW coal—burning utiliyt boiler,the SVM model which had been trained forecasted hte unburned carbon in the ny ash in the test smaples set,and got the mena square root error nad the mean relaitve error of 1.39%,and 1.30%,resepctively,which are 22.20%and 21.07%of BP newtork mode1.These re— suhs show that SVM method is more accurate than the BP neural network,and Can satisfy the forecastign demand wel1. Keywords:utiilyt boiler;unburned carbon content;suppo ̄vector machine 实时、在线、准确地检测飞灰含碳量将有利于监测 较多等问题,且存在造价高、维修保养要求高等缺点. 锅炉燃烧,降低煤耗,提高机组运行的经济性和安全 而传统的软测量技术都是基于人工神经网络方法来建 性.但燃煤锅炉飞灰含碳量受到多种因素的影响(如 模,由于神经网络的学习算法采用经验风险最小化原 煤种、设计安装水平及锅炉运行操作水平等),这就给 理,缺乏定量分析与机理完备的理论结果,因而基于人 飞灰含碳量的实时检测带来了困难.目前,国内外有较 工神经网络燃煤电站锅炉飞灰含碳量软测量方法存在 多关于燃煤电站锅炉飞灰含碳量在线监测的研究,主 很大的局限性. 要有锅炉飞灰测碳仪技术Ll 和软测量技术 应用 针对这一现状,笔者将支持向量机方法应用于燃 于在线预测.飞灰测碳仪的在线测量,在技术上遇到了 煤电站锅炉飞灰含碳量预测领域.支持向量机是Vap— 样品管堵灰、附加设备复杂和对测量精度的影响因素 nik等人根据统计学理论提出的一种新的机器学习方 收稿日期:2005—10—10. 作者简介:蔡杰进(1977一 ),男,博士研究生.联系人:马晓茜,epxqraa@scut.edu.ca 维普资讯 http://www.cqvip.com 2006年8月 蔡杰进等:基于SVM的燃煤电站锅炉飞灰含碳量预测 ・313・ 法,它是建立在统计学习理论的VC维理论和结构风 险最小原理基础上的,能较好地解决小样本、非线性、 高维数和局部极小点等实际问题【7-10 ̄.目前,支持向 量机已被成功应用于电网、天然气等领域的负荷预 测 .燃煤电站锅炉飞灰含碳量预测也可看作是一 种对飞灰含碳量特性及其影响因子间的复杂的非线性 函数关系的逼近问题,笔者尝试将支持向量机理论引 入燃煤锅炉飞灰含碳量预测建模之中. 任何函数只要满足Mercer条件¨ ,都可作为核 函数,较常用的核函数有线性函数多项式函数、高斯函 数和Sigmoid函数等.本文中选取高斯函数作为核 函数,即 唧(- ) (6) 式中, 是高斯核函数的宽度参数. 2基于SVM的燃煤电站锅炉飞灰含碳量预 1支持向量机 支持向量机通过引入核函数用一非线性映射将输 入空间映射到高维的特征空间,利用一个定义在高维 空间中的线性函数集来进行回归估计,目标函数除经 验风险外还包括一个用来控制回归函数复杂性的正则 化部分 H . 给定样本集{(Xi,Y。)} (Xi∈R 为输入变量的 值,Y ∈R为相应的输出值,Ⅳ为样本个数,I"t为输入 变量的维数).SVM所用拟合函数为 Y= )=∑ 咖 +b (1) {((b( ),Y )} 是从输入空间到高维特征空间的非 线性映射,系数{ } 和b由最小化方程式来估算, 即 l 1+c (), 置)) (2)i =1 y )= 一 s (3) 式中:第1项÷II I I为正则化部分;第2项是经验 风险,由式(3)给出的不敏感损失函数 来度量;C为 用来平衡两者的权重参数,称为正则化参数. 为了得到系数 和b,引入松弛变量 和 ,将 式(2)转化为 w II z+C・ ( 一 rYi—W。咖( )一6≤ + i (4) s.t.{W・咖( )+b—Y ≤ + L ’≥0, ≥0 最后,引入拉格朗日乘子,回归函数式(1)可变成精确 形式为 )=∑( i—Ixi*)K(置, )+b (5) 式中: i、 为对偶参数. 测模型 2.1输入、输出参数 燃煤锅炉飞灰含碳量特性受到如煤种、运行参数 和锅炉设计、制作、安装等因素的影响,关系较为复杂. 但如将锅炉考虑成黑箱模型,由于锅炉已经建成运行, 其设计和安装参数均已确定,因此将锅炉负荷、炉膛出 口烟气含氧量、各二次风挡板开度、燃尽风挡板开度、 燃料风挡板开度、煤种特性及各磨煤机给煤量、炉膛与 风箱差压、一次风总风压、燃烧器摆角等运行参数作为 模型的输入量.由于现役机组中存在各种不同型号的 燃煤锅炉,它们之间的一些参数,如磨煤机数目、燃烧 器数目等,存在差异,即使是同一型号的锅炉,不同生 产厂家制造,这些参数也会存在一定的差异,所以在实 际应用工程中,输入参数还需根据具体的实际锅炉模 型进一步细化.飞灰含碳量作为模型的输出. 2.2算法参数 利用SVM方法进行燃煤电站锅炉飞灰含碳量预 测,为了得到准确的预测精度,还必须精心选择SVM 的超参数.在采用高斯核函数的SVM中,超参数包括 以下3项: 1)s不敏感损失函数中s的值越大,支撑向量的 数目就越少,因而解的表达就越稀疏,然而大的 会降 低回归函数的逼近精度,即它被用来平衡解的稀疏程 度和回归函数的逼近精度; . 2)正则化参数c决定着模型复杂度与训练误差 在目标函数中的比重; 3)高斯核函数的宽度参数艿隐含地定义了从输 入空间到高维特征空间的非线性映射,决定了特征空 间的结构,因而控制了最终解的复杂性. 2.3步骤 应用SVM进行燃煤电站锅炉飞灰含碳量预测的 步骤如下: 1)读取数据,并对数据进行平滑处理(剔除奇异 值)和归一化处理(防止计算溢出); 2)形成、更新训练样本集,预测样本集; 3)选取合适的结构参数 、C和艿; 维普资讯 http://www.cqvip.com 燃 烧 科 学 与 技 术 第12卷第4期 尽次 次 风风 风 燃二 一 4)在训练样本上建立形如式(4)的目标函数; 5)利用序列极小化方法求解目标函数; 6)将预测样本的输入变量和求解目标函数所得 的参数代人式(5),用于预测该燃煤电站锅炉某一工 况下的飞灰含碳量. l口口 回 口 回△ = ∑l” 1 l Ll×100% (8) fh l , ,v l c 一 l 式中,C 为C 的估计值. 3.3输入参数细化 一一一一一一一一该300 MW锅炉采用滑压运行方式,采用发电负 3实例分析 荷和人炉总空气量描述锅炉负荷对飞灰含碳量特性的 影响;5台磨煤机对应5个磨煤机给煤量,用以描述一 本文将SVM算法模型应用于某300 MW燃煤电 站锅炉实际运行过程中的飞灰含碳量预测,并和BP 算法的预测进行比较研究. 3.1 对象简介 该机组锅炉是上海锅炉厂生产的SG一1025/18.3一 M317型亚临界中间再热控制循环汽包锅炉,采用正压 直吹式制粉系统;配有5台RP863型中速磨煤机,炉 膛尺寸为14 022 mm x 12 330 mm.炉膛四角布置切向 摆动燃烧器,燃烧器为一二次风间隔布置,燃烧器可在 上、下方向摆动以控制炉膛出口烟温及调节再热蒸汽 温度,采用同心反切燃烧系统,二次风偏离一次风25。 送人炉膛.满负荷条件下投用4层一次风,其余一层一 次风备用.锅炉燃烧器布置如图1所示.图中,OFA代 表顶部燃尽风,该风室高度440 mln,A、B、C、D、E分别 表示各层燃料风,AA、AB、BC、CD、DE和EF分别表示 各层辅助风,括号内数据为风室高度,单位为mm. . j / 八/ ,’蓉 J● / 、 ! 图1燃烧器布置示意 3.2评价指标 采用了以下两个统计量来评价模型的拟合程度和 诊断效果,这两个统计量是均方根误差△ (root mean square eITOr)、平均相对误差△MRE(mean relative eITor) ARMsE 次风粉量沿炉高分配对炉内燃烧的影响;5台磨煤机 对应5个煤粉细度 帅,用以描述煤粉细度对飞灰含碳 量特性的影响;二次风共投用6层,同层联动,底层AA 层二次风门开度固定,因此有5个二次风挡板开度值, 以描述二次风配风方式对燃烧的影响;1个燃尽风挡 板开度值以描述燃尽风的影响;取炉膛与风箱差压测 点的平均值作为1个输人参数以考察其对二次风速的 影响;取炉膛出口烟气含氧量平均值作为1个输人参 数,描述燃烧氧量的影响;煤种特性取发热量Q 、挥 发分 、灰分A 和水分 共4个参数,以描述煤种 特性的影响;各层燃烧器联锁摆动,对应1个燃烧器摆 角;各燃料风开度联锁,对应1个燃料风开度;取1个 一次风总风压以考虑一次风压对一次风速的影响.共 有27个输人参数.按等截面法在磨煤机出口每根煤粉 管道上取样,测得煤粉细度R帅;其他运行参数均可由 电厂DCS中获得. 于是,可确定样本集{(置,Y )} 。的结构,Y 是某 一工况下的飞灰含碳量;置是该工况下锅炉的其他运 行参数. 3.4热态测试结果构造样本集 在该300 MW燃煤锅炉上进行了热态多工况实炉 测试,对影响锅炉燃烧特性的各种运行参数:锅炉负 荷、磨煤机组合、燃尽风(OFA)配风方式、二次风配风 方式、煤种和燃烧器摆角等,进行变工况试验,具体工 况及测得的飞灰含碳量见表1及表2. 取表1和表2中工况1、工况3~5、工况7~10、工 况12~13、工况15~23对应的19个样本构成16个训 练样本;工况2、工况6、工况l1及工况14对应的4个 样本作为检验样本. 3.5算法参数选择 用SVM算法模型进行燃煤电站锅炉飞灰含碳量 预测时,采用高斯函数作为核函数,参数选择分别为 =5,c=10 000和 =0.000 1;用BP神经网络模型 进行燃煤电站锅炉飞灰含碳量预测时,采用标准的3 层BP网络,选用27个输入节点,1个输出节点,32个 隐节点. 维普资讯 http://www.cqvip.com 2006年8月 蔡杰进等:基于SVM的燃煤电站锅炉飞灰含碳量预测 ・315・ 表1现场热态试验工况表(I) 负荷/ 入炉总风量G/ 磨煤机给煤量q/(t・h ) 磨煤机煤粉细度R90/% 风箱炉膛压差 炉膛出口烟气 工况 (W・MW ) (t・h ) A B C D E A B C D E z ̄/Pa 氧量[02]/% 1 3o2 1 145 26 26 26 25 0 25.0 26.0 27.6 19.4 0 496.5 3.3 2 304 1 16o 25 25 25 23 0 25.9 26.8 27.6 2O.7 0 516.5 3.05 3 30o 1 O11 24 24 24 23 0 25.9 26.8 28.6 2O.0 0 506 3.85 4 293 1 141 24 24 24 24 O 25.9 26.8 28.6 2O.0 0 542.5 4.7 5 273 762 22 22 22 23 0 26.8 27.7 29.6 2O.7 0 516 3.15 6 270 915 22 22 21 22 0 26.8 27.7 29.6 21.4 0 7l1.5 3.75 7 271 958 23 23 22 21 0 26.8 26.8 29.6 21.4 0 599 4.35 8 239 656 2O 2O 2O 20 0 27.8 28.7 3O.8 21.4 0 441 3.15 9 241 810 2O 2O 2O 2O 0 27.8 27.7 29.6 21.4 0 489 4.o5 1O 239 850 2O 2O 2O 20 0 26.8 27.7 3O.8 21.4 0 556 4.7 11 292 912 26 26 26 24 0 25.0 27.5 32.0 31.0 0 388.5 3.05 12 288 992 25 26 24 24 O 25.O 26.O 27.2 14.0 0 756.5 3.15 13 289 945 25 25 24 24 0 26.8 29.7 27.6 19.4 0 758.5 2.8 14 292 1 014 25 24 24 O 24 26.8 29.7 28.6 O 33.2 746 2.85 15 293 988 25 24 24 O 24 25.9 26.8 27.6 0 18.0 793 2.9 16 298 99o 25 25 25 0 24 25.9 26.8 27.6 0 22.6 874.5 2.9 17 297 968 25 25 25 O 24 25.9 26.O 27.6 O 22.6 917.5 2.8 18 295 946 24 24 24 O 23 25.O 26.O 26.7 O 21.9 911 3.1 19 295 964 25 25 25 24 0 25.0 26.0 26.7 19.4 0 640 3.15 2O 280 945 25 25 25 24 0 25.O 26.O 27.6 19.4 0 726 3.35 21 294 949 26 26 25 25 0 24.2 25.2 26.7 18.8 O 666 2.7 22 296 988 26 26 25 25 O 24.2 25.2 26.7 18.8 0 835.5 2.85 23 29o 988 26 26 25 25 0 24.2 25.2 26.7 l8.8 0 794.5 3.25 表2现场热态试验工况表(II) 燃尽风开度 煤种特性 二次风门开度 2/% 1,o/ ̄/o 燃烧器 燃料风 一次风 飞灰含 工况 摆角 开度 总压 碳量 肘 A t,O/AB BC CD DE EF 0FA o1% 1 p/kPa Cr. . (kJ・kg ) l 11.7 l6.74 29.2l 22 919 l00 l0o l00 l0o 5 l0o 20 59 7.75 0.84 2 l1.7 13.42 28.13 23 264 10o 10o 10o 10o 5 10o 2O 62 8.15 1.17 3 l1.7 13.42 28.13 23 264 10o 70 70 10o 5 10o 2O 64 8.25 1.47 4 l1.7 13.42 28.13 23 264 10o 70 70 10o 5 91 2O 64 8.15 O.62 5 l1.99 17.59 28.O3 22 349 70 70 10o 10o 5 64 5O 56 7.40 1.28 6 l1.99 17.59 28.O3 22 349 70 70 10o 10o 5 6O 10o 56 7.35 O.62 7 l1.99 17.59 28.O3 22 349 70 70 70 70 5 61 10o 56 7.40 0.70 8 12.1 17.3O 28.34 22 421 70 70 70 70 5 19 10o 44 6.65 1.27 9 l2.1 17.3O 28.34 22 421 10o 10o 10o 10o 5 2l 10o 45 6.65 0.84 1O l2.1 l7.3O 28.34 22 421 10o 10o 10o 10o 5 19 5O 45 6.55 0.93 l1 l3 15.13 29.69 22 692 10o 5O 5O 10o 5 89 10o 53 7.45 1.12 12 l3 15.13 29.69 22 692 10o 5O 5O 10o 5 84 5O 53 7.35 O.5O 13 l3.69 16.64 29.37 21 672 5O 5O 10o 10o 5 85 5O 53 7.5O 0.45 14 l3.69 16.64 29.37 21 672 10o 10o 10o 5 10o 89 10o 59 7.75 0.86 15 l3.69 16.64 29.37 21 672 10o 10o 10o 5 10o 91 5O 52 7.3O 0.43 16 l2.O5 21.O5 28.26 21 055 10o 5O 5O 5 10o 97 10o 53 7.40 O.31 17 l2.05 21.05 28.26 21 O55 1O0 50 5O 5 10o 96 10o 40 7.1O 0.42 18 12.05 21.05 28.26 21 055 1oo 1oo 5O 5 5o 93 1oo 51 7.25 0.43 19 9.63 16.11 27.61 23 964 10o 10o 10o 10o 5 93 5O 6o 7.85 3.20 2O 9.63 16.11 27.61 23 964 10o 5O 5O 10o 5 73 10o 62 7.95 1.44 21 9.63 16.11 27.61 23 964 5O 5O 10o 10o 5 92 5O 61 7.95 2.97 22 9.63 16.11 27.61 23 964 10o 10o 5O 5O 5 95 10o 61 8.O5 1.24 23 9.63 16.11 27.61 23 964 10o 5O 10o 5O 5 87 10o 61 8.o5 2.48 注:由于现场热态试验是在锅炉燃烧及制粉系统性能优化调整刚完成不久进行的,所以该炉的运行状态较优,飞灰含碳量较低。 维普资讯 http://www.cqvip.com 燃 烧 科 学 与 技 术 第12卷第4期 3.6结果及分析 SVM模型和BP网络模型两个模型训练学习过程 及其误差曲线如图2所示,预测输出的相对误差分别 为0.001 5%及0.001 7%,由于两种模型的训练样本输 出和实测值几乎完全一致,所以图中预测值中几乎看 不到“+”号;4个检验样本的预测输出及相对误差如 表3所示,均方根误差△ MsE分别为1.39%和6.26%,平 均相对误差△MRE分别为1.30%和6.17%. (a)训练与预测过程曲线 2 \ 械 (b)相对误差曲线 图2模型训练过程与误差曲线 表3 SVM与BP网络法预测比较 % BP网络法 SVM方法 工况点 实测值 相对误差 相对误差 预测值 绝对值 预测值 绝对值 2 1.17 1.1l2 4.957 1.152 1.538 6 0.62 0.664 7.097 0.624 0.645 , 11 1.12 1.182 5.536 1.142 1.964 14 O.86 0.799 7.O93 O.851 1.047 平均 6.170 1.3oo 均方根误差 6.26 1.39 由以上结果可知,两种预测模型训练过程输出误 差均很小,说明均获得了正确的输入一输出相应行为, 建立了锅炉运行参数和飞灰含碳量关系的神经网络模 型和支持向量机模型.在验证模型正确性的过程中, SVM模型预测误差均比BP网络方法小,均方根误差 △ 眦和平均相对误差△ 分别是BP网络方法的 22.20%和21.07%.说明用SVM方法来预测电站锅炉 的飞灰含碳量,其数据吻合程度相当好,明显优于BP 网络方法,而且精度能够满足实际使用的需要.这主要 是因为SVM模型能充分利用训练样本的分布特性,根 据部分训练样本构建判别函数,不需过多的先验信息 和使用技巧,并最终转化为二次寻优问题,得到的将是 全局最优解,且SVM模型的小样本处理能力较强,在 样本数不多的情况下也能很好地进行统计学习. 4结论 (1)由于支持向量机方法以统计学习理论作为其 理论基础,采用结构风险最小化原则,具有精度高、全 局最优及泛化性好等特点,所以本文提出了基于SVM 的燃煤电站锅炉飞灰含碳量预测模型.该模型可将电 站锅炉飞灰含碳量特性与各个运行参数之间的复杂关 系及参数之间的耦合关系简单地表示出来,可准确、有 效地对飞灰含碳量进行预报. (2)将SVM飞灰含碳量预测方法应用于300 MW 燃煤电站锅炉中.经过训练后的SVM模型对检验样本 飞灰含碳量进行预报,均方根误差△ 和平均相对误 差△MRE分别为1.39%、1.30%,相当于BP网络模型的 22.20%和21.07%.表明SVM方法有更高的预测精度, 更好的泛化能力,具有很好的推广应用前景. (3)SVM模型的建立,为电站锅炉飞灰含碳量的 实时测量提供了一种更为准确的手段,为实时监测电 站锅炉燃烧经济性提供便利.如结合全局寻优算法,可 在燃用任何煤种条件下寻找出最优的操作参数,以获 得高的燃烧经济性. (4)此作为一种新的尝试,SVM方法不仅为电站 锅炉飞灰含碳量预测带来了新的思路和方法,在电力 系统中其他带有更加复杂条件的预测分析方面也有非 常重要的意义. 参考文献: [1]Ouazzane A K,Castagner J L,Jones A R,et 1a.Design of an optical instrument to measure the carbon content of lfy ash [J].Fuel,2002,81(15):1907—1911. [2]方湘涛,宣益民,冯长青,等.MCE-1型烟道式飞灰含碳 量在线测量系统的研制[J].工程热物理学报,2003,24 (3):523—.527. Fang Xinagtao,Xuan Yimin,Feng Changqing,et a1.Devel・ opment of MCE・-1 carbon・-in・-sah measuring system on・-line non-extranctive[J].Journal ofEngineering Thermophysics, 维普资讯 http://www.cqvip.com 2006年8月 蔡杰进等:基于SVM的燃煤电站锅炉飞灰含碳量预测 2003,14(5):993—1009. ・317・ 2003,24(3):523—527(in Chinese). [3] 刘鸿,周克毅.锅炉飞灰测碳仪的技术现状及发展趋 [9]王锡淮,朱思锋.基于支持向量机的船舶电力负荷预测 [J].中国电机工程学报,2004,24(10):36—39. Wang Xihuai,Zhu Sifeng.Ship power load forecasting using 势[J].锅炉技术,2004,35(2):65-_68. Liu Hong,Zhou Xeyi.Recent research and devolopment of unburned carbon in fly ash detector[J].Boiler Technology, 2004,35(2):65—68(in Chinese). support vector machine[J].Proceedings of te CSEE,h 2004,24(10):36—39(in Chinese). [4] 周吴,朱洪波,曾庭华,等.基于人工神经网络的大型 [10]李元诚,方廷健,于尔铿.短期负荷预测的支持向量机 方法研究[J].中国电机工程学报,2003,23(6):55— 电厂锅炉飞灰含碳量建模[J].中国电机工程学报, 2002,22(6):96—100. Zhou Hao,Zhu Hongbo,Zeng TinghHa,et a1.Artiifcial neural network modelling on the unburned carbon in fly ash from utility boilers[J].Proceedings ofthe CSEE,2002,22 (6):96—100(in Chinese). [5] 方湘涛,叶念渝.基于BP神经网络的电厂锅炉飞灰含碳 量预测[J].华中科技大学学报:自然科学版,2003,31 (12):75—77. Fang Xiangtao,Ye Nianyu.A systemforforecasting the an- burned carbon of hte fly ash from utiilyt boilers based on BP artiifcial neural netwoks[J].Jouma/Huazhong Uruverstiy ofScience&Technology:Nature Sciecne Edition,2003,31 (12):75—77(in Chinese). [6] Marta Sebastia,Inaki Fernandez Olmo,Angel Irabien.Neu. ral network prediction of unconfined compressive strength of coal lfy sah.cement mixtures『J].Cement and Concrete Re. search,20o3,33(8):Il37一l145. [7] Lin Kuanming,Lin Chihen.A study on reduced support vector machines『J].IEEE Transactoim on Neural Net. works,2003,14(6):1449—1459. [8] nAsuit ̄D,Boni A,Ridella S.A digital architecture for support vector machines:Theory,algorithm,and FPGA im- plementation[J].IEEE Transactoins on Neural Networks, 59. Li Yuancheng,Fang Tingjina,Yu Erkeng. Study of sup- oprt vector machines for sh0rt-term load forecasting[J]. rPoceedigns foteh CSEE,2003,23(6):55—59(in Chi- nese). [11]Zhu Sifeng,Wang Xihua1.A SVM approach to ship power load forecasting based on RBF kernel[c].IEEE rPoceed- igns of teh 5th Worm Congress OFt Intelligent Control and Automation,Hangzhou,China:2004:1 824--1 828. [12]Pal Pin ̄eng,Hong Weichang.Support vector macihnes with simulated annealing lagorithms in electriciyt load fore— casting[J].Eneryg Conversion and Management,2005, 46(17):2669--2688. [13] Liu Han,Liu Ding,Liang Yamning,et a1.Research on natural gas load forecasting based on least squares support vector machine[J].IEEE Proceedigns ofteh Third Inter- natoinal Conference On Machine Learrlign and Cybeenrtics, 2004(5):26—29. [14] Vapnik V N.The Nature of Statistical Learnign Theory [M].New York:Springer,1999. [15]Smola A,Scholkopf B.A tutorila on support vector regres. sion[R].London:Rc’yal Holloway College,1998.