CN104008425B - 基于引力搜索的水火电系统多目标调峰方法 - Google Patents
基于引力搜索的水火电系统多目标调峰方法 Download PDFInfo
- Publication number
- CN104008425B CN104008425B CN201410198218.2A CN201410198218A CN104008425B CN 104008425 B CN104008425 B CN 104008425B CN 201410198218 A CN201410198218 A CN 201410198218A CN 104008425 B CN104008425 B CN 104008425B
- Authority
- CN
- China
- Prior art keywords
- power station
- individual
- thermal power
- vol
- output
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 230000005484 gravity Effects 0.000 title abstract 3
- 238000005457 optimization Methods 0.000 claims abstract description 16
- 238000010845 search algorithm Methods 0.000 claims abstract description 6
- 238000003860 storage Methods 0.000 claims description 29
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 11
- 230000000694 effects Effects 0.000 claims description 8
- 230000005611 electricity Effects 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 8
- 230000009897 systematic effect Effects 0.000 claims description 8
- 230000001133 acceleration Effects 0.000 claims description 7
- 239000000203 mixture Substances 0.000 claims description 7
- 230000005619 thermoelectricity Effects 0.000 claims description 6
- 238000011144 upstream manufacturing Methods 0.000 claims description 6
- 101100390771 Danio rerio fitm1l gene Proteins 0.000 claims description 3
- 238000005520 cutting process Methods 0.000 description 4
- 238000000205 computational method Methods 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 240000002853 Nelumbo nucifera Species 0.000 description 1
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种水火电系统多目标调峰的引力搜索优化方法，涉及电力系统调度技术领域。该方法首先建立水火电系统多目标调峰模型，模型为具有多个约束条件的目标函数，然后采用以Pareto优化理论为基础的多目标引力搜索算法优化求解多目标调峰模型，获得水火电系统对目标调峰的非劣方案集，得到同时满足系统运行费用最小与水电调量最大两个目标的最优解集，供调度决策部门选择。
Description
技术领域
本发明涉及电力系统调度技术领域，尤其涉及一种在调度过程中智能优化算法。
背景技术
社会生产生活中每天各时段对电力的需求是有差异的，存在峰谷差。水火电力系统调峰就是依据水电站与火电厂的不同特点通过安排各电厂各时段的输出功率的在使电网各时段的输出功率与社会需求负荷相等。
现有电力系统调度技术一般按照电量平衡原则通过切负荷的方式确定水电站和火电厂出力，调峰目标单一。在处理多目标调峰问题时通过目标加权，也将多目标转化为单目标，一次求解只能得到一个调度方案，缺乏优化机制。
下列方法为现有的方法：
(1)火电厂运行费用计算方法：
其中：为为阀点效应造成的叠加值；
参考文献：覃晖，周建中.基于多目标文化差分进化算法的水火电力系统优化调度.电力系统保护与控制，2011,39-22
(2)引力搜索算法(GSA)中关于引力常数G的计算方法，同时借鉴了算法中关于位置矢量X、引力引起的加速度a，和速度向量V的概念；
参考文献：徐遥，王士同.引力搜索算法的改进.计算机工程与应用，2011,47(35)
(3)水电站输出功率计算方法；
(4)水火电系统联合运行过程中的约束处理方法；
(5)Pareto非劣解集理论；
(6)非劣解集拥挤距离截断方法。
发明内容
本发明的目的是，提供基于引力搜索的水火电系统多目标调峰方法，该方法首先建立水火电系统多目标调峰模型，模型为具有多个约束条件的目标函数，然后采用多目标引力搜索算法优化求解多目标调峰模型，获得水火电系统对目标调峰的非劣方案集，供调度决策部门选择。
为了达到上述目的，本发明采用如下技术方案：
基于引力搜索的水火电系统多目标调峰方法，该方法包含下列步骤：
a、水火电系统多目标调峰模型的建立
a1、目标函数一
其中：minF1为水火电系统运行费用最小的目标函数，T为调度期的时段数，Ns为火电厂的个数，为第t时段第i个火电厂的输出功率，为火电运行费用函数，Ei,t为阀点效应函数；
a2、目标函数二
其中：minF2为水火电系统水电调峰量最大的目标函数，Nh为水电站的个数，PD,t为第t时段的水火电系统的负荷需求，为第i个水电站在第t时段的输出功率；
a3、约束条件
a31、系统电力负荷平衡约束
其中：PD,t为各时段的预测负荷，Ns为火电厂的个数，为第t时段第i个火电厂的输出功率，t＝1,2，…，T，i＝1,2，…，Ns；Nh为水电站的个数，为第t时段第i个水电站的输出功率，t＝1,2,…,T，i＝1,2,…,Nh；
a32、火电厂与水电站的输出功率限制
其中：和分别第为第i个火电厂最小与最大输出功率，和分别为第i个水电站最小与最大输出功率；
a33、水电站的发电流量约束
Qi,min≤Qi,t≤Qi,max (6)
其中：Qi,t为第i个水电站第t时段的发电流量，Qi,min与Qi,max分别为第i个水电站最小与最大发电流量；
a34、水电站的库容限制
VOLi,min≤VOLi,t≤VOLi,max (7)
其中：VOLi,min与VOLi,max分别为第i个水电站最小与最大库容，VOLi,t为第i个水电站第t时段的库容；
a35、水电站调度期的始末库容约束
其中：VOLi,B与VOLi,E分别为第i个水电站调度期始末库容限制，VOLi,0为第i个水电站调度期计算开始时的库容，VOLi,T为第i个水电站调度期计算结束时的库容；
a36、火电输出功率的爬坡速度约束
其中：为第i个火电厂第t时段的输出功率，为第i个火电厂第t-1时段的输出功率μi，为第i个火电厂的爬坡速度；
b、用多目标引力搜索算法优化求解步骤a中的多目标调峰模型，步骤如下：
b1、的初始值为满足约束条件的随机数，i＝1,2,…,Ns，t＝1,2,…,T；Qi,t的初始值为满足Qi,min≤Qi,t≤Qi,max约束条件的随机数，i＝1,2,…,Nh，t＝1,2,…,T，初始化设定阈值ε和最大迭代次数Max_it；
b2、按照步骤b1的要求初始化第i个水电站第t时段的发电流量Qi,t，i＝1,2,…,Nh，t＝1,2,…,T和第i个火电厂第t时段的输出功率i＝1,2,…,Ns，t＝1,2,…,T，将初始值按如下方式排列：
组成的行向量称为个体；
b3、重复步骤b2操作产生N个个体组成的集合称为一个群体X，其中第i个个体用Xi表示,初始化个体Xi的速度向量Vi为零向量；
b4、将个体Xi中第i个水电站第t时段的发电流量代入公式 其中为上游第h个水电站的发电流量，τh为上游第h个水电站的流达时间，nh为上游水电站水电个数；Ii,t为第i个水电站第t时段的自然入库流量。求出第i个水电站第t时段的水库库容VOLi,t，i＝1,2,…,Nh，t＝1,2,…,T，对发电流量Qi,t和水电站库容VOLi,t进行约束处理；进而通过查询水库库容与水头关系曲线得出第i个水电站第t时段的水头Hi,t，i＝1,2,…,Nh，t＝1,2,…,T；将Qi,t与Hi,t代入公式其中：Ki为第i个水电站的出力系数，求出第i个水电站第t时段的输出功率i＝1,2,…,Nh，t＝1,2,…,T；将代入目标函数2中计算其函数值fit2；
b5、对个体Xi中第i个火电厂第t时段的输出功率进行系统电力负荷平衡约束处理，i＝1,2,…,Ns，t＝1,2,…,T；代入目标函数1中计算其函数值fit1；
b6、重复按照步骤b4与步骤b5求出群体X中所有个体的两个目标函数值；
b7、根据群体X中个体的两个目标函数值筛选出非劣个体，当非劣个体数量超过NP时采用拥挤距离方法进行截断，直至非劣个体数到达NP，将这NP个非劣个体的集合称为外部档案集Y，其中第j个个体用Yj表示；
b8、计算群体中个体Xi受外部档案集Y吸引产生的加速度，更新个体速度和位置，计算方法如下：
b8.1、求出外部档案集Y中第j个个体Yj，j＝1,2,…,NP，到群体中个体Xi的欧式距离Disti,j，j＝1,2,…,NP，其中设最大欧式距离为worsti，最小欧式距离为besti；
b8.2、根据公式与求出外部档案集Y中个体Yj相对于群体中个体Xi的惯性质量Mi,j，j＝1,2,…,NP；
b8.3、根据公式求出群体中个体Xi受外部档案集Y的吸引产生的加速度ai，rand为(0,1)之间的一个随机数，G为引力常数；
b8.4、根据公式Vi＝rand·Vi+ai更新个体Xi速度向量，由公式Xi＝Xi+Vi更新个体Xi；
b8.5、重复步骤b8.1—步骤b8.5直到群体X中所有个体在外部档案集Y的吸引下全部更新；
b9、重复步骤b6，求出更新后群体X全部个体的两个目标函数值，筛选出外部档案集Y与更新后群体共同的非劣个体，如果非劣个体数量超过NP，采用拥挤距离方法进行截断，直至非劣个体数到达NP，这个非劣个体集合组成新的外部档案集Y；
b10、在按照设定的最大迭代次数内，重复步骤b8和步骤b9进行迭代操作；
b11、每迭代10次计算出第n次迭代产生外部档案集与第n-10次迭代产生的外部档案集对应个体的两个目标函数值的差值Δfit1j和Δfit2j，j＝1,2,..,NP，根据公式计算出相对差值Δfitj；求出相对差值的平均值，如果该平均值小于设定的阈值ε停止迭代，否则重复步骤10直到达到最大迭代次数Max_it；
b12、计算结束，输出外部档案集作为优化结果，外部档案集中的个体代表不同的水火电联合调峰调度方案；如果决策者更偏好系统运行费用的调度目标，则选择系统运行费用最小调度方案，如果偏好水电调峰目标，则选择水电调峰量最大方案。
与现有技术相比，本发明具有以下优点和效果：
1、已有技术多按照电量平衡原则通过切负荷的方式确定水电站和火电厂出力，缺乏优化机制，本发明通过智能优化方法优化调峰目标函数，能按照设定目标实现水火电优化协调调峰。
2、已有技术多是考虑水电调峰量单一目标的水火电系统运行方法，本发明是结合了运营成本与水电调峰能力的多目标调峰运行方法。且采用多目标优化算法在优化求解，一次求解可得到满足Pareto最优理论的多个非劣解，获得多种优化调峰运行方案，有助于提高决策效率。
3、已有技术在处理多目标调峰问题时通过目标加权将多目标转化为单目标，一次求解只能得到一个调度方案。本发明利用Pareto最优理论求解得到同时满足系统运行费用最小与水电调峰量最大两个目标的最优解集，得到多个调度方案供决策者选择。
具体实施方式
以下结合实施例，对本发明作进一步的说明。
实施例一
基于引力搜索的水火电系统多目标调峰方法，该方法包含下列步骤：
步骤1，数据准备
调度期为1天，调度期时段T＝24；各时段系统负荷PD,t；水电站个数Nh，各水电站最大出力值最小出力值各水电站最大库容值VOLi,max、最小库容值VOLi,min，各水电站调度期初始库容值VOLi,B和期末库容值VOLi,E，各水电站最大发电流量值Qi,max，最小发电流量Qi,min，各水电站水位库容关系曲线，各水电站出力系数Ki，各水电站各时段自然入库径流Ii,t；火电厂的个数Ns，火电厂运行费用函数阀点效应函数Ei,t，各火电厂最大出力值最小出力值各火电厂出力爬坡速度限制值μi；
步骤2，建立多目标调峰模型
(2.1)系统运行费用最小目标：
其中：为为阀点效应造成的叠加值；为t时段第i个火电厂的出力，ai,bi和ci为运行费用函数系数。
研究表明，忽略阀点效应会影响运行费用的求解精度，本发明在计算系统运行成本最小目标时考虑了阀点效应。
(2.2)水电调峰量最大目标：
其中：为第i个水电站在t时段的出力。调峰电量最大模型的目标就是使控制期内的最大余荷最小。
所描述的目标函数需满足如下约束条件：
(2.3)系统电力负荷平衡约束
其中：PD,t为各时段的预测负荷，Ns为火电厂的个数，为第t时段第i个火电厂的输出功率，t＝1,2，…，T，i＝1,2，…，Ns；Nh为水电站的个数，为第t时段第i个水电站的输出功率，t＝1,2，…,T，i＝1,2，…，Nh；
(2.4)火电厂与水电站的输出功率限制
其中：和分别第为第i个火电厂最小与最大输出功率，和分别为第i个水电站最小与最大输出功率；
(2.5)水电站的发电流量约束
Qi,min≤Qi,t≤Qi,max (14)
其中：Qi,t为第i个水电站第t时段的发电流量，Qi,min与Qi,max分别为第i个水电站最小与最大发电流量；
(2.6)水电站的库容限制
VOLi,min≤VOLi,t≤VOLi,max (15)
其中：VOLi,min与VOLi,max分别为第i个水电站最小与最大库容，VOLi,t为第i个水电站第t时段的库容；
(2.7)水电站调度期的始末库容约束
其中：VOLi,B与VOLi,E分别为第i个水电站调度期始末库容限制，VOLi,0为第i个水电站调度期计算开始时的库容，VOLi,T为第i个水电站调度期计算结束时的库容；
(2.8)火电输出功率的爬坡速度约束
其中：为第i个火电厂第t时段的输出功率，为第i个火电厂第t-1时段的输出功率μi，为第i个火电厂的爬坡速度；
步骤3、用多目标调峰的引力搜索优化方法求解步骤2的调峰模型，步骤如下：
(3.1)的初始值为满足约束条件的随机数，i＝1,2,…,Ns，t＝1,2,…,T；Qi,t的初始值为满足Qi,min≤Qi,t≤Qi,max约束条件的随机数，i＝1,2,…,Nh，t＝1,2,…,T，初始化设定阈值ε和最大迭代次数Max_it；
(3.2)按照步骤b1的要求初始化第i个水电站第t时段的发电流量Qi,t，i＝1,2,…,Nh，t＝1,2,…,T和第i个火电厂第t时段的输出功率2,…,T，将初始值按如下方式排列：
组成的行向量称为个体；
(3.3)重复步骤(3.2)操作产生N个个体组成的集合称为一个群体X，其中第i个个体用Xi表示,初始化个体Xi的速度向量Vi为零向量。
(3.4)将个体Xi中所有水电站每个时段的发电流量代入公式 (τh为上游水电站到下游水电站的流达时间)中求出所有水电站的每个时段的水库库容VOLi,t，i＝1,2,…,Nh，t＝1,2,…,T，对发电流量Qi,t和水电站库容VOLi,t进行约束处理。其中水电站调度始末库容约束处理具体操作如下：
(3.4.1)设定约束处理最大循环次数Ny1和阀值α；
(3.4.2)已根据Qi,t计算出调度期末库容求出违反约束量ΔVOL；
(3.4.3)将ΔVOL平均分配到Qi,t中即Qi,t＝Qi,t+ΔVOL/Nh；
(3.4.4)重新计算ΔVOL；
(3.4.5)重复步骤(3.4.2)至(3.4.4)的操作直到循环次数达到Ny1或ΔVOL<α。
(3.4.6)根据上步求出的各时段水电站库容VOLi,t查询水电站库容与水头关系曲线图得出所有水电站每个时段的水头Hi,t。将Qi,t与Hi,t代入公式中(Ki为水电站i的出力系数)，求出公式所有水电站各时段的输出功率i＝1,2,…,Nh，t＝1,2,…,T。将代入目标函数2中计算其函数值fit2。
(3.5)对个体Xi中所有火电厂每个时段的输出功率进行系统电力负荷平衡约束处理，其中系统电力负荷平衡约束处理方法如下：
(3.5.1)设定约束处理最大循环次数Ny2和阀值β；
(3.5.2)计算系统t时刻所有水电站与火电厂出力之和与系统负荷需求PD,t之间差值为ΔP；
(3.5.3)将ΔP平均分配到中，即：
(3.5.4)重新计算ΔP；
(3.5.5)重复步骤(3.5.2)至(3.5.4)的操作直到循环次数达到Ny2或ΔP<β。
将约束处理后得到的代入目标函数1中计算其函数值fit1。
(3.6)重复按照步骤(3.4)、(3.5)求出群体X中所有个体的两个目标函数值。
(3.7)根据群体X中个体的两个目标函数值筛选出非劣个体。当非劣个体数量超过NP时采用拥挤距离方法进行截断，直至非劣个体数到达NP，将这个非劣个体的集合称为外部档案集Y，其中第j个个体用Yj表示。
(3.8)计算群体中个体Xi受外部档案集Y吸引产生的加速度，更新个体速度和位置，具体计算如下：
(3.8.1)求出外部档案集Y中全部个体Yj，j＝1,2,…,NP，到群体中个体Xi的欧式距离Disti,j，j＝1,2,…,NP，其中设最大欧式距离为worsti，最小欧式距离为besti；
(3.8.2)根据公式与求出外部档案集Y中个体Yj相对于群体中个体Xi的惯性质量Mi,j；
(3.8.3)重复按照(3.8.2)操作求出外部档案集中所有个体相对于群体中个体Xi的惯性质量Mi,j，j＝1,2,…,NP；
(3.8.4)根据公式求出群体中个体Xi受外部档案集的吸引产生的加速度ai，rand为(0,1)之间的一个随机数，G为引力常量；
(3.8.5)根据公式Vi＝rand·Vi+ai更新个体Xi速度向量，由公式Xi＝Xi+Vi更新个体Xi；
(3.8.6)重复步骤(3.8.1)—(3.8.5)直到群体X中所有个体在外部档案集Y的吸引下全部更新。
(3.9)重复步骤(3.6)操作求出更新后群体X全部个体的两个目标函数值，筛选出外部档案集Y与更新后群体共同的的非劣个体，如果非劣个体数量超过NP，采用拥挤距离方法进行截断，直至非劣个体数到达NP，这个非劣个体集合组成新的外部档案集。
(3.10)在按照设定的最大迭代次数内，重复步骤(3.7)和步骤(3.8)进行迭代优化搜索。
(3.11)每迭代10次计算出第n次迭代产生的外部档案集与第n-10次迭代产生的外部档案集对应个体的两个目标函数值的差值Δfit1j和Δfit2j，j＝1,2,..,NP，根据公式计算出相对差值Δfitj。求出相对差值的平均值，如果该平均值小于设定的阈值ε停止迭代，否则重复步骤(3.10)直到达到最大迭代次数Max_it。
(3.12)计算结束，输出外部档案集作为优化结果，外部档案集中个体代表不同水火电联合调峰调度方案。如果决策者更偏好系统运行费用的调度目标，则选择系统运行费用最小调度方案，如果偏好水电调峰目标，则选择水电调峰量最大方案。
其中，非劣关系判断，设置外部档案集与截断操作是以Pareto理论为基础的多目标优化技术典型特征。
Claims (1)
1.一种水火电系统多目标调峰的引力搜索优化方法，其特征在于，该方法包含下列步骤：
a、水火电系统多目标调峰模型的建立
a1、目标函数一
其中：minF1为水火电系统运行费用最小的目标函数，T为调度期的时段数，Ns为火电厂的个数，为第t时段第i个火电厂的输出功率，为火电运行费用函数，Ei，t为阀点效应函数；
a2、目标函数二
其中：minF2为水火电系统水电调峰量最大的目标函数，Nh为水电站的个数，PD,t为第t时段的水火电系统的负荷需求，为第i个水电站在第t时段的输出功率；
a3、约束条件
a31、系统电力负荷平衡约束
其中：Ns为火电厂的个数，为第t时段第i个火电厂的输出功率，t＝1,2，…，T，i＝1,2，…，Ns；Nh为水电站的个数，为第t时段第i个水电站的输出功率，t＝1,2，…,T，i＝1,2，…，Nh；
a32、火电厂与水电站的输出功率限制
其中：和分别第为第i个火电厂最小与最大输出功率，和分别为第i个水电站最小与最大输出功率；
a33、水电站的发电流量约束
Qi,min≤Qi,t≤Qi,max (5)
其中：Qi,t为第i个水电站第t时段的发电流量，Qi,min与Qi,max分别为第i个水电站最小与最大发电流量；
a34、水电站的库容限制
VOLi,min≤VOLi,t≤VOLi,max (6)
其中：VOLi,min与VOLi,max分别为第i个水电站最小与最大库容，VOLi,t为第i个水电站第t时段的库容；
a35、水电站调度期的始末库容约束
其中：VOLi,B与VOLi,E分别为第i个水电站调度期始末库容限制，VOLi,0为第i个水电站调度期计算开始时的库容，VOLi,T为第i个水电站调度期计算结束时的库容；
a36、火电输出功率的爬坡速度约束
其中：为第i个火电厂第t时段的输出功率，为第i个火电厂第t-1时段的输出功率，μ i 为第i个火电厂的爬坡速度；
b、用多目标引力搜索算法优化求解步骤a中的多目标调峰模型，步骤如下：
b1、的初始值为满足约束条件的随机数，i＝1、2、…、Ns，t＝1、2、…、T；Qi,t的初始值为满足Qi,min≤Qi,t≤Qi,max约束条件的随机数，i＝1、2、…、Nh，t＝1、2、…、T，初始化设定阈值ε和最大迭代次数Max_it；
b2、按照步骤b1的要求初始化第i个水电站第t时段的发电流量Qi,t，i＝1、2、…、Nh，t＝1、2、…、T和第i个火电厂第t时段的输出功率i＝1、2、…、Ns，t＝1、2、…、T，将初始值按如下方式排列：
组成的行向量称为个体；
b3、重复步骤b2操作产生N个个体组成的集合称为一个群体X，其中第i个个体用Xi表示,初始化个体Xi的速度向量Vi为零向量；
b4、将个体Xi中第i个水电站第t时段的发电流量代入公式其中为上游第h个水电站的发电流量，τh为上游第h个水电站的流达时间，nh为上游水电站水电个数；Ii,t为第i个水电站第t时段的自然入库流量，求出第i个水电站第t时段的水库库容VOLi，t，i＝1、2、…、Nh，t＝1、2、…、T，对发电流量Qi,t和水电站库容VOLi，t进行约束处理；进而通过查询水库库容与水头关系曲线得出第i个水电站第t时段的水头Hi,t，i＝1、2、…、Nh，t＝1、2、…、T；将Qi,t与Hi,t代入公式其中：Ki为第i个水电站的出力系数，求出第i个水电站第t时段的输出功率i＝1、2、…、Nh，t＝1、2、…、T；将代入目标函数2中计算其函数值fit2；
b5、对个体Xi中第i个火电厂第t时段的输出功率进行系统电力负荷平衡约束处理，i＝1、2、…、Ns，t＝1、2、…、T；代入目标函数1中计算其函数值fit1；
b6、重复按照步骤b4与步骤b5求出群体X中所有个体的两个目标函数值；
b7、根据群体X中个体的两个目标函数值筛选出非劣个体，当非劣个体数量超过NP时采用拥挤距离方法进行截断，直至非劣个体数到达NP，将这NP个非劣个体的集合称为外部档案集Y，其中第j个个体用Yj表示，NP为外部档案集设定规模；
b8、计算群体中个体Xi受外部档案集Y吸引产生的加速度，更新个体速度和位置，计算方法如下：
b8.1、求出外部档案集Y中第j个个体Yj，j＝1，2，3…NP，到群体中个体Xi的欧式距离Disti,j，j＝1，2，3…NP，其中设最大欧式距离为worsti，最小欧式距离为besti；
b8.2、根据公式与求出外部档案集Y中个体Yj相对于群体中个体Xi的惯性质量Mi,j，j＝1，2，3…NP；
b8.3、根据公式求出群体中个体Xi受外部档案集Y的吸引产生的加速度ai，rand为(0,1)之间的一个随机数，G为引力常数；
b8.4、根据公式Vi＝rand*Vi+ai更新个体Xi速度向量，由公式Xi＝Xi+Vi更新个体Xi；
b8.5、重复步骤b8.1—步骤b8.5直到群体X中所有个体在外部档案集Y的吸引下全部更新；
b9、重复步骤b6，求出更新后群体X全部个体的两个目标函数值，筛选出外部档案集Y与更新后群体共同的非劣个体，如果非劣个体数量超过NP，采用拥挤距离方法进行截断，直至非劣个体数到达NP，这个非劣个体集合组成新的外部档案集Y；
b10、在按照设定的最大迭代次数内，重复步骤b8和步骤b9进行迭代操作；
b11、每迭代10次计算出第n次迭代产生外部档案集与第n-10次迭代产生的外部档案集对应个体的两个目标函数值的差值Δfit1j和Δfit2j，j＝1,2,..,NP，根据公式计算出相对差值Δfitj；求出相对差值的平均值，如果该平均值小于设定的阈值ε停止迭代，否则重复步骤10直到达到最大迭代次数Max_it；
b12、计算结束，输出外部档案集作为优化结果，外部档案集中的个体代表不同的水火电联合调峰调度方案；如果决策者更偏好系统运行费用的调度目标，则选择系统运行费用最小调度方案，如果偏好水电调峰目标，则选择水电调峰量最大方案。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410198218.2A CN104008425B (zh) | 2014-05-12 | 2014-05-12 | 基于引力搜索的水火电系统多目标调峰方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410198218.2A CN104008425B (zh) | 2014-05-12 | 2014-05-12 | 基于引力搜索的水火电系统多目标调峰方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104008425A CN104008425A (zh) | 2014-08-27 |
| CN104008425B true CN104008425B (zh) | 2017-05-17 |
Family
ID=51369073
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410198218.2A Active CN104008425B (zh) | 2014-05-12 | 2014-05-12 | 基于引力搜索的水火电系统多目标调峰方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104008425B (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104504474A (zh) * | 2014-12-23 | 2015-04-08 | 上海电机学院 | 一种微网能量多目标优化方法 |
| CN104533701B (zh) * | 2014-12-23 | 2017-03-08 | 华中科技大学 | 一种水轮机调速系统控制参数的自动整定方法 |
| CN104538992A (zh) * | 2015-01-14 | 2015-04-22 | 云南电力调度控制中心 | 一种大小水电及风电协调优化调度方法 |
| CN105868841B (zh) * | 2016-03-21 | 2019-07-16 | 广西电网有限责任公司电力科学研究院 | 一种基于风电优先上网的风水火联合调度方法 |
| CN107451562B (zh) * | 2017-07-31 | 2020-04-24 | 湖北工业大学 | 一种基于混沌二进制引力搜索算法的波段选择方法 |
| CN108345998B (zh) * | 2018-02-09 | 2021-04-27 | 湖北工业大学 | 一种基于杂交水稻算法的水火电经济调度方法及系统 |
| CN110598919B (zh) * | 2019-08-28 | 2022-04-01 | 华中科技大学 | 一种梯级水电站群生态调控方法与系统 |
| CN110601233B (zh) * | 2019-09-30 | 2023-02-21 | 国家电网公司西北分部 | 一种电力系统中储能电站的调峰调度方法 |
| CN111461478B (zh) * | 2020-01-17 | 2022-04-22 | 华中科技大学 | 一种大规模水光能源互补调度方法与系统 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103049800A (zh) * | 2012-12-17 | 2013-04-17 | 上海大学 | 一种存储时间受限的自动化立体仓库调度多目标优化方法 |
| CN103489131A (zh) * | 2013-08-29 | 2014-01-01 | 国家电网公司 | 一种基于光柴储供电系统的运行调度方法 |
| CN103728881A (zh) * | 2014-01-02 | 2014-04-16 | 东南大学 | 一种多楼宇冷热电联供系统的优化运行方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2013026731A1 (en) * | 2011-08-22 | 2013-02-28 | Abb Technology Ag | System and method to optimize operation of a water network |
-
2014
- 2014-05-12 CN CN201410198218.2A patent/CN104008425B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103049800A (zh) * | 2012-12-17 | 2013-04-17 | 上海大学 | 一种存储时间受限的自动化立体仓库调度多目标优化方法 |
| CN103489131A (zh) * | 2013-08-29 | 2014-01-01 | 国家电网公司 | 一种基于光柴储供电系统的运行调度方法 |
| CN103728881A (zh) * | 2014-01-02 | 2014-04-16 | 东南大学 | 一种多楼宇冷热电联供系统的优化运行方法 |
Non-Patent Citations (1)
| Title |
|---|
| 李超顺 等.基于引力搜索核聚类算法的水电机组振动故障诊断.《中国电机工程学报》.2013,第33卷(第2期), * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104008425A (zh) | 2014-08-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104008425B (zh) | 基于引力搜索的水火电系统多目标调峰方法 | |
| CN109523059B (zh) | 一种梯级水电站水库生态调度智能优化方法及系统 | |
| CN106682810B (zh) | 巨型水电站动态投产下跨流域梯级水电站群长期运行方法 | |
| CN105574605A (zh) | 一种梯级泵站输水系统日优化调控方法 | |
| CN109345019B (zh) | 一种基于改进粒子群算法的微电网经济调度优化策略 | |
| CN110598919B (zh) | 一种梯级水电站群生态调控方法与系统 | |
| CN109447393A (zh) | 一种电力系统经济负荷分配的改进的粒子群算法 | |
| CN107609683B (zh) | 一种基于萤火虫算法的梯级水库群调度优化方法 | |
| CN106169109A (zh) | 一种基于混沌差分粒子群算法的水电站优化调度方法 | |
| CN109345010A (zh) | 一种梯级泵站的多目标优化调度方法 | |
| CN104636821A (zh) | 基于动态惯性权重粒子群的火电机组负荷优化分配方法 | |
| CN113471989B (zh) | 基于灰狼优化法的微电网智能二次频率控制方法 | |
| Rahi et al. | Optimization of hydro power plant design by particle swarm optimization (PSO) | |
| CN104504455B (zh) | 一种梯级蓄能控制下水电站群长期优化调度方法 | |
| CN104252132A (zh) | 基于自适应遗传算法的行星际轨道控制优化方法 | |
| CN116757446B (zh) | 基于改进粒子群算法的梯级水电站调度方法及系统 | |
| CN104965954B (zh) | 一种基于全特性空间曲面建模的梯级负荷经济分配方法 | |
| CN104216383B (zh) | 一种小水电站机组运行效率优化方法 | |
| CN114254571B (zh) | 一种抽水蓄能电站极端工况下机组控制规律优化决策方法 | |
| CN109214110B (zh) | 一种长距离提水工程优化调度方法 | |
| CN111612170A (zh) | 一种输电设备的检修计划制定方法及系统 | |
| CN111461478B (zh) | 一种大规模水光能源互补调度方法与系统 | |
| CN117610828A (zh) | 一种步长逐密多目标梯级水电站群优化调度方法 | |
| CN103088783A (zh) | 一种面向生态的水库调度方案生成方法 | |
| CN109409569B (zh) | 一种考虑直流输电约束的水库群中长期调度的离散梯度逐步优化算法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |