一种基于SPA和TOPSIS耦合的工业园区负荷灵活度量化方法

一种基于spa和topsis耦合的工业园区负荷灵活度量化方法
技术领域
1.本发明涉及工业园区负荷领域，具体涉及一种基于spa和topsis耦合的工业园区负荷灵活度量化方法。


背景技术：

2.能源是人类生存和发展必不可少的物质基础，是全部人类社会发展和进步的最根本的物质保障。目前，中国已经成为世界上最大的能源消耗国。但同时，大量化石燃料的燃烧带来的环境问题已经严重制约着我国的经济、能源和环境之间的协调发展，社会生产力的跨越式发展给传统的能源供应模式带来了巨大的挑战。积极有效开发利用清洁能源成为我国建设清洁、安全和高效的现代能源体系和实现非化石能源发展目标的重要措施。
3.近年来，风电的大量接入造成了系统调峰资源的紧张，传统的调度模式已经无法接纳全部的风电，寻求并合理利用除常规电源以外的调峰资源成为亟待解决的问题。在大容量用户与电网之间形成有效的互动机制，是解决问题的可行思路。
4.需求响应能力对负荷调节资源非常关键，用电客户里有非常多潜在需求响应资源还没有被挖掘。作为电力系统需求侧的重要组成部，工业园区负荷管理，利用用户的用电灵活性,来缓解负荷高峰时的供电紧张状况,有利于电力系统的安全经济运行,削弱电力市场中市场势力的影响,抑制价格尖峰。
5.通过研究发现，影响工业园区负荷灵活度的因素包括政策层面上和技术层面上和的因素。在政策层面上，主要影响因素包括电价、相关补贴政策等；在技术层面上，主要影响因素包括系统的调峰能力、调频能力、备用水平、负荷特性、电网的输电容量等。但大部分研究着重于供给侧的各种因素对可再生能源消纳的制约与影响。实践研究表明，电网具有不确定性及黑箱性，即电网内部复杂多变，影响因素繁多，只研究需求侧的影响因素未免过于片面。
6.目前，工业园区负荷占总负荷比例越来越重，对工业园区的研究成为热点。研究发现工业园区负荷具有足够的灵活性，能够在很短的时间内，通过调节电压对不断增加或减少能耗的要求做出反应，且不影响产品的质量。因此，它是一种很有前途的资源，以满足工业和电网的需求。一系列工业园区负荷灵活度的评估模型层出不穷。相对的，大量关于评估工业园区负荷灵活度的文献也随之出现，但目前的研究主要集中在定性评估工业园区负荷的灵活性，评估结果并不准确、不够直观。如何量化评估工业园区负荷的灵活度的文献较少，并且大多只应用需求侧的线性指标进行评估，面对复杂多变的大电网，这样的评估指标并不能有效的揭示工业园区负荷的灵活度。


技术实现要素：

7.发明目的
8.本发明提出一种基于spa和topsis耦合的工业园区负荷灵活度量化方法，可有效解决电网企业在量化评估工业园区负荷灵活度时评估结果并不准确、不够直观的问题。
9.技术方案
10.一种基于spa和topsis耦合的工业园区负荷灵活度量化方法，其特征在于：按照以下步骤执行：
11.步骤一：采集工业负荷参与调度后大电网的供给侧与需求侧的电气数据，包括电气运行参数、电能产出量和电能需求量；
12.步骤二：根据影响工业园区负荷灵活度的因素，得出评价指标，构建工业园区灵活度评估体系，将步骤一中采集的数据带入到评估体系中，得出不同评价指标下的运行模拟结果；
13.步骤三：将步骤二中的各评价指标下的运行模拟结果导入到spa和topsis耦合的工业园区负荷灵活度量化评估模型中，得到标准化矩阵；
14.步骤四：根据步骤三中得到的标准化矩阵，用改进的组合赋权法确定各评价指标的权重；
15.步骤五：将步骤四中确定的各评价指标的权重与步骤三得到的标准化矩阵相乘得到决策矩阵；
16.步骤六：将步骤五得到的决策矩阵代入topsis模型得到决策矩阵的最优解、中间解和最劣解；
17.步骤七：运用spa理论中的多元联系度，计算决策矩阵中的每一列要素与第六步中所得的最优解、中间解和最劣解的同一度、差异度和对立度；
18.步骤八：根据步骤七中的同一度、差异度和对立度计算指标之间的联系度，得到工业园区负荷灵活度的量化评估模型。
19.在步骤二所述的评估指标方面包括：经济性、优质性、安全性、高效性；经济性包含非线性的电价敏感度和用电成本；优质性包含可调节深度、响应偏差率和故障率；高效性包含最大可调节速率、最大可调节时间和能耗率；安全性包含用电安全记录、用电违规记录和调节电压对电网稳定性的影响，量化工业园区负荷的灵活度。
20.提出非线性电价敏感度，并且综合考虑建成成本、运维成本和折损成本的用电成本，量化工业园区负荷的灵活度；
21.指标c1：非线性的电价敏感度：
[0022][0023]
式中：a为工业园区负荷i能够有所转移最小电价差，b为工业园区负荷i极限转移率所对应的最大电价差，λ
gd
为该负荷的极限转移率，x为实时电价差，k为转移系数，取2；
[0024]
指标c2：用电成本：
[0025]
[0026]
式中：t
l
为工业园区负荷i的平均设备寿命期；为工业园区负荷i在第t年的设备利用小时数；c
i,t
为工业园区负荷i在第t年的装机规模；分别为工业园区负荷i的建成成本、运维成本和折损成本。
[0027]
在步骤四中所述的改进的组合赋权法，引入灰色关联系数对熵权法进行改进，考虑指标之间的关联程度，并且将改进后的熵权法和层次分析法相组合，以得到各个指标的权重；
[0028]
在步骤四中所述的改进的组合赋权法如下：
[0029]
①
改进熵权法确定工业园区负荷灵活度的客观权重；
[0030]
数据选取后，用初始矩阵a表示，针对具有不同的量纲和数量级的评价指标体系，数据进行规范化处理，以实现实际值到指标评价值的转化；采用初值化、均值化和区间相对化的方法；
[0031][0032]
对初始矩阵a进行正向化、无量纲、规范化处理后，得到矩阵b：
[0033][0034]
参考矩阵取：
[0035]
b0＝{max(b
11
,b
12
...b
1n
),max(b
21
,b
22
...b
2n
),...max(b
m1
,b
m2
...b
mn
)}
t
ꢀꢀꢀ
(5)
[0036]
灰色关联系数的计算：
[0037][0038]
式中：i＝1,2,
…
,m；为m个评价指标，j＝1,2,
…
，n；为n个工业园区负荷用户，ρ取0.5；
[0039]
对矩阵c进行标准化的矩阵z
[0040][0041]
计算第j项指标下第i个样本值占该指标的比重：
[0042]
[0043]
计算各项指标的熵：
[0044][0045]
式中，k＝1/ln(n)>0
[0046]
计算各项指标的权重：
[0047][0048][0049][0050]
式中：为所有e
j
不等于1的平均值，ω
i
为第i个指标对应的权重，ω
i1
，ω
i2
为第i个指标的分权重；
[0051]
②
层次分析法确定权重：
[0052][0053]
式中：b
ij
为因素指标b
i
对因素指标b
j
的重要性，υ
i
为层次分析法所确定的第i个指标所对应的权重；
[0054]
最终的各个指标的组合权重：
[0055][0056]
式中：w
i
为第i个指标的组合权重；
[0057]
根据上述公式，可以得到各项评价指标的综合权重向量：
[0058]
w＝(w1,w2,...,w
m
)
t
ꢀꢀꢀ
(15)
[0059]
在步骤六中将决策矩阵代入topsis模型得到决策矩阵的最优解、中间解和最劣解；
[0060]
采用改进组合赋权法的spa和topsis耦合后的工业园区负荷灵活度量化模型，将所有评价指标正向化后，构建正向化矩阵：
[0061][0062]
式中：x
ij
为第j个工业园区负荷的第i个评价指标；
[0063]
将正向化矩阵标准化，构建标准化矩阵z，z中的每一个元素为：
[0064][0065]
构建加权决策矩阵：将无量纲矩阵z
ij
与各指标对应的权重相乘，得到加权决策矩阵：
[0066]
v
ij
＝w
i
×
z
ij
，(i＝1，2，
…
，m；j＝1，2，...，n)
ꢀꢀꢀ
(18)
[0067]
确定评价体系的正理想解s
max
、负理想解s
min
和中间理想解s
mid
；s
max
为m个工业园区负荷各评价指标的最优值组成的集合，s
min
为m个工业园区负荷各评价指标的最劣值组成的集合，s
mid
为m个工业园区负荷各评价指标数据的平均值；
[0068]
计算工业园区负荷a的灵活度；计算工业园区负荷a
j
与理想解s
max
、中间理想解s
mid
和最劣理想解s
min
的欧氏距离；
[0069][0070]
式中：d
jmax
、d
jmid
、d
jmin
为决策矩阵中的第j个工业园区负荷与理想解s
max
、中间理想解s
mid
和最劣理想解s
min
的欧氏距离。
[0071]
在步骤八中，引入spa对topsis模型进行改进，从同、异、反三个方面考虑工业园区负荷的灵活度；
[0072]
在步骤八中所述的基于spa和topsis耦合的工业园区负荷灵活度量化模型为：
[0073]
u(a
j
,a0)＝a+bi
j
+cj
ꢀꢀꢀ
(20)
[0074]
[0075][0076]
式中j＝
‑
1，u(a
j
,a0)为第j个工业园区负荷的灵活度。
[0077]
优点及效果：
[0078]
1、考虑了非线性的电价敏感度和用电成本。传统灵活度模型评估模型主要是评估误差大的主要原因在于采用线性的需求侧指标，然而多种柔性负荷的园区配电网内部参数复杂，非线性因素繁多，具有复杂的非线性关系。且传统的量化评估模型仅从需求侧考虑了可调节容量、可维持时长、经济性等指标，电网内部复杂多变，复杂因素繁多，柔性条件的改变能够引起电网内部众多的连锁反应，导致需求侧与供给侧的数据变化并不是简单的线性关系，故仅考虑需求侧，而进行量化评估并不能有效的揭示出工业园区负荷的灵活度。本发明将非线性的电价敏感度和用电成本考虑到影响工业园区负荷的灵活度的评价指标内，能够更加准确地量化工业园区负荷的灵活度。
[0079]
2、建立了基于spa和topsis耦合的工业园区负荷灵活度量化模型。传统topsis模型在计算待评价对象与正、负理想解的距离时，如果评价对象与正、负理想解的距离都很近，则不能准确真实地反应各评价对象的灵活度。本发明引入集对分析理论(spa)和中间解。从同、异、反三个角度准确的反映出评价对象与各解之间的联系度，使得工业园区负荷灵活度量化评估方法更加客观，结果更加直观，更具有说服力。
附图说明
[0080]
附图1为工业园区负荷灵活度量化方法总体流程图；
[0081]
附图2为工业园区负荷灵活度评价指标图；
[0082]
附图3为改进组合赋权法总体流程图；
[0083]
附图4为灵活度求解逻辑关系图。
具体实施方式
[0084]
本方法是一种基于spa和topsis耦合的工业园区负荷灵活度量化模型。首先，确定了工业园区负荷的灵活度的评价指标和约束条件；其次，利用改进的组合赋权算法对各归一化后的评价指标进行赋权；再次；利根据各个归一化后的指标和权重得到对工业园区负荷的灵活度决策矩阵；最后，基于spa和topsis耦合模型得到对工业园区负荷的灵活度量化模型。
[0085]
本方法特征在于：该方法的步骤如下：
[0086]
步骤一：采集工业负荷参与调度后大电网的供给侧与需求侧的电气运行参数、电能产出量和电能需求量等电气数据；
[0087]
步骤二：根据影响工业园区负荷灵活度的因素，得出评价指标，构建工业园区灵活度评估体系，将第一步中采集的数据带入到评估体系中，得出不同评价指标下的运行模拟结果；
[0088]
步骤三：将第二步中的各评价指标下的运行模拟结果导入到spa和topsis耦合的工业园区负荷灵活度量化评估模型中得到标准化矩阵；
[0089]
步骤四：根据步骤三中得到的工业园区负荷的灵活度量化评估模型中得到的标准化矩阵，用改进的组合赋权法确定各评价指标的权重；
[0090]
步骤五：将第四步中的确定的各评价指标的权重与第三步得到的标准化矩阵相乘得到决策矩阵；
[0091]
步骤六：将第五步得到的决策矩阵代入topsis模型得到决策矩阵的最优解中间解和最劣解；
[0092]
步骤七：运用spa理论中的多元联系度计算决策矩阵中的每一列要素与第六步中所得的最优解、中间解和最劣解的同一度、差异度和对立度；
[0093]
步骤八：根据第七步中的同一度、差异度和对立度计算指标之间的联系度，得到工业园区负荷灵活度的量化评估模型。
[0094]
步骤二所述的评估指标包括经济、优质、安全、高效这四个方面，其中经济性包含电价敏感度和用电成本；优质性包含可调节深度、响应偏差率和故障率；高效性包含最大可调节速率、最大可调节时间和能耗率；安全性包含用电安全记录、用电违规记录和调节电压对电网稳定性的影响，量化工业园区负荷的灵活度。在本方法中首次提出非线性的电价敏感度和用电成本两个指标，可有效解决在电网企业在选择工业园区负荷用户作为峰谷时段市场备用时的评估问题的公约园区负荷的灵活度量化方法。
[0095]
所述的四个方面的评价指标和对应的约束条件的定义如下：
[0096]
经济性包含电价敏感度和用电成本。
[0097]
指标c1：非线性的电价敏感度：
[0098][0099]
式中：a为工业园区负荷i能够有所转移最小电价差，b为工业园区负荷i极限转移率所对应的最大电价差，λ
gd
为该负荷的极限转移率。x为实时电价差。
[0100]
指标c2：用电成本：
[0101][0102]
式中：t
l
为工业园区负荷i的平均设备寿命期；为工业园区负荷i在第t年的设备利用小时数；c
i,t
为工业园区负荷i在第t年的装机规模。分别为工业园区负荷i的建成成本、运维成本和折损成本。
[0103]
优质性包含可调节深度、响应偏差率和故障率。
[0104]
指标c3：可调节深度：
[0105][0106]
式中：为调节前t时刻的标准化后负荷用电功率；为调节后t时刻的标准化负荷用电功率
[0107]
指标c4：响应偏差率：
[0108][0109]
式中：p
i,need
为第i个工业园区负荷所需要响应的量，p
i,sj
为第i个工业园区负荷所实际响应的量。
[0110]
指标c5：故障率：
[0111][0112]
式中：j
i
为工业园区负荷i的故障停机时间，t
i,z
为工业园区负荷i的额定运行时间。
[0113]
高效性包含负荷的最大可调节速率、最大可调节时间、能耗率。
[0114]
指标c6：最大可调节速率：
[0115][0116]
式中；δ
i
为工业园区负荷i的调节速率；λ＝1表示负荷进行上调节，λ＝0表示负荷进行下调节；p
rated
为可调节负荷的额定值；p
max
、p
min
分别为负荷可调节能力上、下限。
[0117]
指标c7：最大可调节时间：
[0118]
t
adj,max
＝t
ret
‑
t
on
＝max(t
ret,i
‑
t
on,i
)
[0119]
式中：t
on
为初始调节时间；t
ret
为返回时间；t
ret,i
t
ret,i
分别为工业园区负荷i的返回时间、初始调节时间。
[0120]
指标c8：能耗率：
[0121][0122]
式中：w
z
为第i个工业园区负荷的能源消耗量，w
c
为第i个工业园区负荷的产量。
[0123]
安全性包含用电安全记录、用电违规记录和调节电压对电网稳定性的影响。
[0124]
指标c9：参与电网需求响应的完成率：
[0125][0126]
式中：r
i,n1
为第i个工业园区负荷完成电网需求响应的次数，r
i,n
为第i个工业园区负荷参与电网需求响应的总次数。
[0127]
指标c
10
：参与电网需求响应的违规率：
[0128]
r
i,v
＝1
‑
r
i,se
[0129]
式中：r
i,se
为第i个工业园区负荷参与电网需求响应的完成率。
[0130]
指标c
11
：调节电压对电网稳定性的影响：
[0131][0132]
式中：ui为线路过负荷程度，ui值越高，电压冲击性越强，对电力系统安全运行影响越大；u
i
为第i个工业园区负荷的调节电压，为第i个工业园区负荷的调节电压的平均值。
[0133]
约束条件如下：
[0134]
(1)电压限制
[0135]
u
min
≤u
t
≤u
max
[0136]
u
max
代表可调节电压的最大值；u
min
代表可调节电压的最小值。
[0137]
(2)热平衡约束
[0138]
热平衡是提供灵活性关键的问题,如果温度过高，冰晶石会融化到罐壁上，腐蚀罐。为了保证生产过程的安全，需要将温度限制在一定范围内。热平衡约束可以表示为：
[0139][0140]
式中：为t
keep
小时维持温度所需的最小输入能量。q
0i,t'
为时间段时第i个工业园区负荷的功率。
[0141]
(3)调节持续时长约束：
[0142]
t
lh,min
≤t
keep
≤t
lh,max
[0143]
t
lh,min
、t
lh,max
分别为可调节的允许的调节持续最小、最大时长。
[0144]
(4)可调节负荷功率约束：
[0145]
0≤δq
lh
≤δq
lhm
[0146]
式中：δq
lhmax
为调节工业园区负荷最大上调功率。
[0147]
(5)环保约束：
[0148]
化石能源产生的各类排放满足区域相关环保目标要求：
[0149][0150]
式中：为工业园区负荷i在第t年的单位碳排放系数；为工业园区负荷i在第t年的度电耗能系数；为第t年的碳排放上限。
[0151]
在步骤四中所述的改进的组合赋权法如下：
[0152]
1)层次分析法确定权重：
[0153]
构造判断矩阵：根据三标尺尺度法对评价指标的重要性进行评分，并构建判断矩阵b。
[0154][0155]
式中：b
ij
为因素(指标)b
i
对因素(指标)b
j
的重要性，b
ij
>0，b
ij
＝1/b
ji
。
[0156]
权重向量为：
[0157][0158]
一致性检验，判断矩阵最大n次根的计算：
[0159][0160]
一致性检验：ri通过查表可得，cr＜0.1则可以通过一致性检验。
[0161]
2)改进熵权法确定工业园区负荷灵活度的客观权重
[0162]
数据选取后，用初始矩阵a表示，由于评价指标体系具有不同的量纲和数量级，因此需要对数据进行规范化处理，以实现实际值到指标评价值的转化。通常可采用初值化、均值化和区间相对化等方法。
[0163][0164]
对初始矩阵a进行正向化、无量纲、规范化处理后，得到矩阵b：
[0165][0166]
参考矩阵取：
[0167]
b0＝{max(b
11
,b
12
...b
1n
),max(b
21
,b
22
...b
2n
),...max(b
m1
,b
m2
...b
mn
)}
t
[0168]
灰色关联系数的计算：
[0169][0170]
式中：i＝1,2,
…
,m；为m个评价指标，j＝1,2,
…
，n；为n个工业园区负荷用户。ρ取0.5。
[0171]
对矩阵c进行标准化的矩阵z
[0172][0173]
计算第j项指标下第i个样本值占该指标的比重：
[0174][0175]
计算各项指标的权重：
[0176][0177]
式中，k＝1/ln(n)>0
[0178][0179]
式中：为所有e
j
不等于1的平均值。
[0180][0181][0182]
最终的各个指标的组合权重：
[0183][0184]
根据上述公式，可以得到各项评价指标的综合权重向量：
[0185]
w＝(w1,w2,...,w
m
)
t
[0186]
以往对于工业园区负荷的灵活度的分析居多，但没有建立一个直观的定量的标准进行评估，且对于各个评价指标的权重采用人为，主观赋权，没有客观性，改进的组合赋权法解决了这一缺点，使得在对各评价指标进行赋权时既考虑了主观经验，又考虑了数据的客观性，同时引入了灰熵算法，考虑了指标之间的相关性，使得评价指标的权重更客观，准确。
[0187]
在步骤六中将决策矩阵代入topsis模型得到决策矩阵的最优解、中间解和最劣
解：
[0188]
采用改进组合赋权法的spa和topsis耦合后的工业园区负荷灵活度量化模型，将所有评价指标正向化后，构建正向化矩阵：
[0189][0190]
式中：x
ij
为第j个工业园区负荷的第i个评价指标。
[0191]
将正向化矩阵标准化，构建标准化矩阵z，z中的每一个元素为：
[0192][0193]
构建加权决策矩阵。将无量纲矩阵z
ij
与各指标对应的权重相乘，得到加权决策矩阵：
[0194]
v
ij
＝w
i
×
z
ij
(i＝1,2,...,m；j＝1,2,...,n)
[0195]
确定评价体系的正理想解s
max
、负理想解s
min
和中间理想解s
mid
。s
max
为m个工业园区负荷各评价指标的最优值组成的集合，s
min
为m个工业园区负荷各评价指标的最劣值组成的集合，s
mid
为m个工业园区负荷各评价指标数据的平均值。
[0196]
计算工业园区负荷a的灵活度。计算工业园区负荷a
j
与理想解s
max
、中间理想解s
mid
和最劣理想解s
min
的欧氏距离。
[0197][0198]
在步骤八中所述的基于spa和topsis耦合的工业园区负荷灵活度量化模型为：
[0199]
u(a
j
,a0)＝a+bi
j
+cj
[0200]
[0201][0202]
式中j＝
‑
1，u(a
j
,a0)为第j个工业园区负荷的灵活度。
[0203]
本方法提出一种基于spa和topsis耦合的工业园区负荷灵活度量化方法，改善了传统的量化模型中权重和指标线性加权的缺点，更加准确地量化工业园区负荷的灵活度。
[0204]
下面结合附图，对本发明加以具体描述；
[0205]
为实现上述发明目的，本发明一种基于spa和topsis耦合的工业园区负荷灵活度量化方法，其特征在于，包括以下步骤：
[0206]
如附图1所示，首先，确定了工业园区负荷的灵活度的评价指标和约束条件；其次，利用改进的组合赋权算法对各归一化后的评价指标进行赋权；再次；利根据各个归一化后的指标和权重得到对工业园区负荷的灵活度决策矩阵；最后，基于spa和topsis耦合模型得到对工业园区负荷的灵活度量化模型。
[0207]
具体实施方案如下：
[0208]
如图2所示，首先确定评价指标和约束条件：
[0209]
(1)评价指标：
[0210]
经济性包含电价敏感度和用电成本。
[0211]
指标c1：非线性的电价敏感度：
[0212][0213]
式中：a为工业园区负荷i能够有所转移最小电价差，b为工业园区负荷i极限转移率所对应的最大电价差，λ
gd
为该负荷的极限转移率。x为实时电价差。
[0214]
指标c2：用电成本：
[0215][0216]
式中：t
l
为工业园区负荷i的平均设备寿命期；为工业园区负荷i在第t年的设备利用小时数；c
i,t
为工业园区负荷i在第t年的装机规模。分别为工业园区负荷i的建成成本、运维成本和折损成本。
[0217]
优质性包含可调节深度、响应偏差率和故障率。
[0218]
指标c3：可调节深度：
[0219]
[0220]
式中：为调节前t时刻的标准化后负荷用电功率；为调节后t时刻的标准化负荷用电功率
[0221]
指标c4：响应偏差率：
[0222][0223]
式中：p
i,need
为第i个工业园区负荷所需要响应的量，p
i,sj
为第i个工业园区负荷所实际响应的量。
[0224]
指标c5：故障率：
[0225][0226]
式中：j
i
为工业园区负荷i的故障停机时间，t
i,z
为工业园区负荷i的额定运行时间。
[0227]
高效性包含负荷的最大可调节速率、最大可调节时间、能耗率。
[0228]
指标c6：最大可调节速率：
[0229][0230]
式中；δ
i
为工业园区负荷i的调节速率；λ＝1表示负荷进行上调节，λ＝0表示负荷进行下调节；p
rated
为可调节负荷的额定值；p
max
、p
min
分别为负荷可调节能力上、下限。
[0231]
指标c7：最大可调节时间：
[0232]
t
adj,max
＝t
ret
‑
t
on
＝max(t
ret,i
‑
t
on,i
)
[0233]
式中：t
on
为初始调节时间；t
ret
为返回时间；t
ret,i
t
ret,i
分别为工业园区负荷i的返回时间、初始调节时间。
[0234]
指标c8：能耗率：
[0235][0236]
式中：w
z
为第i个工业园区负荷的能源消耗量，w
c
为第i个工业园区负荷的产量。
[0237]
安全性包含用电安全记录、用电违规记录和调节电压对电网稳定性的影响。
[0238]
指标c9：参与电网需求响应的完成率：
[0239][0240]
式中：r
i,n1
为第i个工业园区负荷完成电网需求响应的次数，r
i,n
为第i个工业园区负荷参与电网需求响应的总次数。
[0241]
指标c
10
：参与电网需求响应的违规率：
[0242]
r
i,v
＝1
‑
r
i,se
[0243]
式中：r
i,se
为第i个工业园区负荷参与电网需求响应的完成率。
[0244]
指标c
11
：调节电压对电网稳定性的影响：
[0245][0246]
式中：ui为线路过负荷程度，ui值越高，电压冲击性越强，对电力系统安全运行影响越大；u
i
为第i个工业园区负荷的调节电压，为第i个工业园区负荷的调节电压的平均值。
[0247]
(2)约束条件：
[0248]
(1)电压限制
[0249]
u
min
≤u
t
≤u
max
[0250]
u
max
代表可调节电压的最大值；u
min
代表可调节电压的最小值。
[0251]
(2)热平衡约束
[0252]
热平衡是提供灵活性关键的问题,如果温度过高，冰晶石会融化到罐壁上，腐蚀罐。为了保证生产过程的安全，需要将温度限制在一定范围内。热平衡约束可以表示为：
[0253][0254]
式中：为t
keep
小时维持温度所需的最小输入能量。q
0i,t'
为时间段时第i个工业园区负荷的功率。
[0255]
(3)调节持续时长约束：
[0256]
t
lh,min
≤t
keep
≤t
lh,max
[0257]
t
lh,min
、t
lh,max
分别为可调节的允许的调节持续最小、最大时长。
[0258]
(4)可调节负荷功率约束：
[0259]
0≤δq
lh
≤δq
lhmax
[0260]
式中：δq
lhmax
为调节工业园区负荷最大上调功率。
[0261]
(5)环保约束：
[0262]
化石能源产生的各类排放满足区域相关环保目标要求：
[0263][0264]
式中：为工业园区负荷i在第t年的单位碳排放系数；为工业园区负荷i在第t年的度电耗能系数；为第t年的碳排放上限。
[0265]
其次，利用改进的组合赋权算法对各归一化后的评价指标进行赋权，如图3所示。
[0266]
1)层次分析法确定权重：
[0267]
构造判断矩阵：根据三标尺尺度法对评价指标的重要性进行评分，并构建判断矩阵b。
[0268][0269]
式中：b
ij
为因素(指标)b
i
对因素(指标)b
j
的重要性，b
ij
>0，b
ij
＝1/b
ji
。
[0270]
权重向量为：
[0271][0272]
一致性检验，判断矩阵最大n次根的计算：
[0273][0274]
一致性检验：ri通过查表可得，cr＜0.1则可以通过一致性检验。
[0275]
2)改进熵权法确定工业园区负荷灵活度的客观权重
[0276]
数据选取后，用初始矩阵a表示，由于评价指标体系具有不同的量纲和数量级，因此需要对数据进行规范化处理，以实现实际值到指标评价值的转化。通常可采用初值化、均值化和区间相对化等方法。
[0277][0278]
对初始矩阵a进行正向化、无量纲、规范化处理后，得到矩阵b：
[0279][0280]
参考矩阵取：
[0281]
b0＝{max(b
11
,b
12
...b
1n
),max(b
21
,b
22
...b
2n
),...max(b
m1
,b
m2
...b
mn
)}
t
[0282]
灰色关联系数的计算：
[0283][0284]
式中：i＝1,2,
…
,m；为m个评价指标，j＝1,2,
…
，n；为n个工业园区负荷用户。ρ取0.5。
[0285]
对矩阵c进行标准化的矩阵z
[0286][0287]
计算第j项指标下第i个样本值占该指标的比重：
[0288][0289]
计算各项指标的权重：
[0290][0291]
式中，k＝1/ln(n)>0
[0292][0293]
式中：为所有e
j
不等于1的平均值。
[0294][0295][0296]
最终的各个指标的组合权重：
[0297][0298]
根据上述公式，可以得到各项评价指标的综合权重向量：
[0299]
w＝(w1,w2,...,w
m
)
t
[0300]
以往对于工业园区负荷的灵活度的分析居多，但没有建立一个直观的定量的标准进行评估，且对于各个评价指标的权重采用人为，主观赋权，没有客观性，改进的组合赋权法解决了这一缺点，使得在对各评价指标进行赋权时既考虑了主观经验，又考虑了数据的客观性，同时引入了灰熵算法，考虑了指标之间的相关性，使得评价指标的权重更客观，准确。
[0301]
最后，基于spa和topsis耦合模型得到对工业园区负荷的灵活度量化模型。
[0302]
如附图4所示，采用改进组合赋权法的spa和topsis耦合后的工业园区负荷灵活度量化模型，将所有评价指标正向化后，构建正向化矩阵：
[0303][0304]
式中：x
ij
为第j个工业园区负荷的第i个评价指标。
[0305]
将正向化矩阵标准化，构建标准化矩阵z，z中的每一个元素为：
[0306][0307]
构建加权决策矩阵。将无量纲矩阵z
ij
与各指标对应的权重相乘，得到加权决策矩阵：
[0308]
v
ij
＝w
i
×
z
ij
(i＝1,2,...,m；j＝1,2,...,n)
[0309]
确定评价体系的正理想解s
max
、负理想解s
min
和中间理想解s
mid
。s
max
为m个工业园区负荷各评价指标的最优值组成的集合，s
min
为m个工业园区负荷各评价指标的最劣值组成的集合，s
mid
为m个工业园区负荷各评价指标数据的平均值。
[0310]
计算工业园区负荷a的灵活度。计算工业园区负荷a
j
与理想解s
max
、中间理想解s
mid
和最劣理想解s
min
的欧氏距离。
[0311][0312]
u(a
j
,a0)＝a+bi
j
+cj
[0313]
[0314][0315]
式中j＝
‑
1，u(a
j
,a0)为第j个工业园区负荷的灵活度。
[0316]
本方法提出一种基于spa和topsis耦合的工业园区负荷灵活度量化方法，改善了传统的量化模型中权重和指标线性加权的缺点，更加准确地量化工业园区负荷的灵活度。
[0317]
以上技术特征构成了本发明的实施例，其具有较强的适应性和实施效果，可根据实际需要增减非必要的技术特征，来满足不同情况的需求。