基于交叉熵重要抽样的元件可靠性参数优化解析方法与流程

1.本发明属于电力系统的模型优化技术领域，涉及基于交叉熵重要抽样的元件可靠性参数优化解析方法。


背景技术：

2.元件可靠性参数优化是指通过优选部分元件采取某些物理强化措施降低被优化元件的故障发生概率，从而降低系统失负荷风险。尽管优化元件可靠性参数可以显著降低元件的故障概率，但是对所有系统元件进行可靠性优化的成本是十分高昂的。因此可靠性参数优化需要考虑可靠性与经济性的平衡，选择合适的元件进行优化。
3.目前可靠性参数优化的研究主要分为两个方面，考虑极端条件下电网的加固规划以及针对正常运行条件下的元件可靠性参数优化。
4.在可靠性参数优化研究中，由于系统可靠性指标并非以元件可靠性参数的解析函数给出，整个优化模型只能采用智能算法求解，并且针对智能算法产生的每个候选解都需要执行一次完整的可靠性评估来评判解的优劣，计算极其耗时。
5.为了解决可靠性参数优化问题在求解过程中需要对变化的可靠性参数进行反复的系统可靠性评估导致计算效率极低这一问题，需要避免因可靠性参数变化而需要多次进行可靠性评估的弊端。因此，如何建立新的元件可靠性参数优化模型是目前本领域技术人员亟待解决的技术问题。


技术实现要素：

6.为解决现有技术中的不足，本技术提供基于交叉熵重要抽样的元件可靠性参数优化解析方法，建立基于交叉熵重要抽样的元件可靠性参数优化解析模型并进行求解，可提升样本代表性和抽样效率，本发明的解析计算精度与传统可靠性优化方法十分接近，在保证计算精度的前提下，极大的提高了可靠性参数优化问题的计算效率。
7.为了实现上述目标，本发明采用如下技术方案：
8.基于交叉熵重要抽样的元件可靠性参数优化解析方法，包括以下步骤：
9.步骤1：输入系统的电气参数、初始可靠性参数，设置系统的元件总个数；
10.步骤2：采用交叉熵重要抽样产生系统失负荷事件，并进行系统可靠性评估，存储失负荷事件、削负荷量、似然比，并累计所有最优抽样事件的eens，得出系统eens指标；
11.步骤3：基于步骤2得到的失负荷事件信息，建立基于交叉熵重要抽样的系统可靠性指标解析模型；
12.步骤4：通过构建辅助变量将系统可靠性指标解析模型融入元件可靠性参数优化问题中，建立可靠性参数优化解析模型；
13.步骤5：将可靠性参数优化解析模型分解为上层问题和下层问题，所述上层问题为决策变量为模型系数和辅助变量的线性规划问题，所述下层问题为决策变量为被优化元件的投资费用的非线性规划问题；
14.步骤6：基于失负荷事件、削负荷量和似然比，针对系统中的所有元件求解上层问题得到解析模型系数，并对系统元件进行重要度排序；
15.步骤7：根据元件重要度排序选取被优化元件集合，依次求解上层问题与下层问题，得到被优化元件投资费用及被优化后的不可用率；
16.步骤8：输出被优化元件集合及其中各元件的投资费用及被优化后的不可用率。
17.本发明进一步包括以下优选方案：
18.优选地，步骤2具体包括：
19.步骤2.1：交叉熵重要抽样初始化，包括设置方差系数收敛阈值β0，每次迭代抽样次数n，初始化迭代次数k＝1，设置初始化最优重要抽样概率密度函数；
20.步骤2.2：预抽样：采用步骤2.1设置的初始化最优重要抽样概率密度函数进行抽样，生成n个系统状态样本；
21.步骤2.3：对每个样本进行可靠性评估，确定削负荷量h(xi)，xi表示第i个样本的系统状态，i＝1,
…
,n，通过公式获得下一次迭代中各元件的最优不可用度q
k,j
，得到下一次迭代采用的最优重要抽样概率密度函数，并进行下一次迭代；
22.当前后两次迭代中各元件不可用度变化的最大值小于或等于0.01时，认为迭代过程收敛，输出最后一次迭代所得各元件的不可用度即为最优不可用度；
23.其中，q
k,j
表示第k次迭代的元件j的最优不可用度的估计值；
24.w(xi；p；q
k-1
)表示第k次迭代的似然比：
[0025][0026]
式中，w(xi；p；q
k-1
)表示第k次迭代的似然比，x
ij
表示系统状态样本xi中元件j的状态，q
k-1,j
为第k-1次迭代中元件j的不可用率；
[0027]
步骤2.4：循环可靠性评估：采用步骤2.3最后一次迭代使用的最优重要抽样概率密度函数进行正式抽样，抽取一个事件x，对抽取的样本进行系统状态分析，计算事件x的失负荷量为h(x)；
[0028]
计算由正式抽样产生的全部随机事件对应的eens指标方差系数β，若β＜β0，转步骤2.5，否则，抽取下一个系统状态进行可靠性评估，直至跳出循环；
[0029]
步骤2.5：存储失负荷事件、削负荷量、似然比：
[0030]
设每个元件有两个状态：正常状态用0表示，故障状态用1表示，则任意失负荷事件对应的元件状态用一个按照元件编号排序后的二进制数表示，记失负荷事件的总数为n
l
；
[0031]
将正式抽样得到的全部失负荷事件对应的二进制状态序列存储在n
l
×
ne维矩阵b中，并将对应的削负荷量、似然函数值分别存储到n
l
维列向量h与w中；
[0032]
其中，b的每一行对应一个失负荷事件的二进制状态序列。
[0033]
步骤2.6：累计所有最优抽样事件的eens，得出系统eens指标。
[0034]
累计所有最优抽样事件的eens，得出系统eens指标：
[0035]
优选地，步骤2.1中，设置初始化最优重要抽样概率密度函数为原始概率密度函数，即：i(x；q0)＝f(x；p)；
[0036]
i(x；q0)为初始化最优重要抽样概率密度函数，用于预抽样；
[0037]
f(x；p)为系统状态向量x的原始概率密度函数；
[0038]
设重要抽样概率密度函数i(x；q)与原始概率密度函数f(x；p)具有相同的形式，其联合概率密度函数表示为：用于求解似然比；
[0039]
其中，q0表示初始抽样概率密度分布的可靠性参数，q表示重要抽样概率密度分布的可靠性参数，qj为元件j的不可用度，x
ij
表示系统状态样本xi中元件j的状态，j＝1,
…
,ne。
[0040]
优选地，步骤3包括以下步骤：
[0041]
步骤3.1：采用状态枚举法得到系统eens指标解析模型；
[0042]
步骤3.2：在步骤3.1解析模型推导基础上，结合步骤2得到的系统eens指标，得到基于交叉熵重要抽样的系统可靠性指标解析模型。
[0043]
优选地，步骤3.1中，记采用状态枚举法得到的系统eens指标为：
[0044][0045]
假设被优化元件的个数为nh，每个元件只有工作和失效两种状态，且各元件之间的故障相互独立，则被优化元件的状态组合共有个；
[0046]
设代表nh个被优化元件处在第j个组合状态的系统事件集合，则系统失效事件集合ψ被分为m个子集：
[0047]
被优化元件的m个状态用ej(j＝1,2,
…
,m)表示；
[0048]
显然，这m个状态彼此互斥且构成一个完备事件组；
[0049]
则根据任意失负荷事件x对应的被优化元件所处的组合状态，将所述采用状态枚举法得到的系统eens指标公式拆分为m项之和，即：
[0050][0051]
将公式(1)改写为：
[0052][0053]
结合条件概率的定义，将公式(2)写为：
[0054][0055]
其中，p(ej)是被优化元件组合状态ej的发生概率；
[0056]
p(x|ej)是给定被优化元件组合状态为ej时系统失负荷事件x发生的条件概率；
[0057]
最终，将公式(3)简写为：
[0058][0059]
其中，对任意j＝1,2,...,m
[0060][0061][0062]
其中，在系统拓扑结构、电气参数和运行参数不变的前提下，p(x|ej)取决于n
e-nh个未被优化元件的可靠性水平；
[0063]
假设n
e-nh个其他元件的可靠性水平不变，则k
1-km为恒定值，与被优化元件的可靠性参数变化无关；
[0064]
φ
ju
和φ
jd
分别是组合状态ej中处于工作和故障状态的被优化元件集合；
[0065]
a对应被优化元件的可用率。
[0066]
优选地，步骤3.2中，根据任意失负荷事件x对应的被优化元件所处的组合状态，将步骤2得到的系统eens指标拆分为m项之和，即：
[0067]
[0068]
当抽样样本足够大时，采用状态枚举法和基于交叉熵重要抽样得到的eens指标是相等的，对比公式(4)与公式(8)，得到：
[0069][0070]
则式(9)求得系数kj：
[0071][0072]
通过式(10)可得基于交叉熵重要抽样的系统eens指标关于被优化元件可用率的解析函数，即为基于交叉熵重要抽样的系统可靠性指标解析模型：
[0073][0074]
优选地，步骤4中，以最小化系统eens指标为优化目标建立可靠性参数优化解析模型：
[0075][0076]
s.t.
[0077][0078]
[0079][0080]
上述可靠性参数优化解析模型为混合整数非线性优化，在模型中，矩阵b、g和p，向量h、w、b、和τ
all
均为已知量；
[0081]
其中，b、h、w是通过执行步骤2得到的失负荷事件、削负荷量、似然函数值输出结果；为系统全部元件的初始不可用率，为全部元件不可用率的取值下限，为全部元件的投资费用系数，τ
all
是对应于全部元件的常系数向量，和τ
all
均为ne维行向量；
[0082]
决策变量为：代表优化决策的0-1变量矩阵y和被优化元件i的投资费用ci；
[0083]
索引i、j和l的取值范围分别为：其中，为系统被优化元件集合；为系统被优化元件的所有组合状态集合，为系统被优化元件的所有组合状态集合，为系统元件集合；
[0084]
y为ne×
nh维表征系统元件是否被选择优化的二进制变量矩阵，其中，ne为系统元件总个数，nh为系统被优化元件个数；若元件l被优化，则矩阵y中第l行中的任意一个元素为1；否则，第l行元素均为0；式(15)是矩阵y自身的逻辑约束；
[0085]
式(13)用于求取优化解析模型系数kj，计算公式如式(13.a)中所示，其中，αj为n
l
维列向量，用于线性相加失负荷事件的失负荷信息；
[0086]
式(13.b)为辅助变量向量αj的下限约束，式(13.b)中构造了nh×ncs
维常数矩阵g、p与n
cs
维常数向量b；
[0087]
g矩阵用于存储nh个被优化元件的n
cs
种组合状态，每个状态可以用一个按照元件编号排序后的二进制数的序列表示，g的每一列对应一个组合状态的二进制数序列；
[0088]
矩阵p的维度与g相同，若g的第i行第j列的元素g
i,j
＝1，则p
i,j
＝1；若g
i,j
＝0，则p
i,j
＝-1，p的第j列记为pj；
[0089]
b为n
cs
维行向量；向量b的第j个元素bj的计算公式如式(13.e)所示；
[0090]
式(13.c)为辅助变量向量αj中每个元素的取值范围，式(13.d)与(14.h)分别为被优化元件组合状态ej的初始概率p
rd
(ej)与元件不可用率变化后的概率p(ej)的计算公式；
[0091]
式(14)表示被优化元件不可用率ui与投资费用ci之间的函数关系，式(14.b)-(14.d)分别表示通过矩阵相乘提取被优化元件对应的投资费用系数c
rd
、不可用率初始值u
rd
和不可用率下限u
min
，其中，，其中，
[0092]
式(14.e)为被优化元件不可用率u(u＝u1,
…
,uh)的上下限约束，由于元件可靠性参数优化只会降低元件不可用率，因此不可用率的上限即为初始值，式(14.e)中τ是个常系数向量，式(14.g)为元件优化费用的总预算约束。
[0093]
优选地，步骤5中，分解得到上层问题为：
[0094][0095]
s.t.
[0096][0097][0098]
上述上层问题用于求取解析模型系数kj；
[0099]
分解得到下层问题为：
[0100][0101]
s.t.
[0102][0103]
优选地，步骤6中，根据系统可靠性指标对各元件投资费用的灵敏度对系统元件进行重要度排序，系统可靠性指标对各元件投资费用的灵敏度计算公式为：
[0104][0105]
优选地，步骤7中，依次求解上层问题与下层问题，得到被优化元件的投资费用，依
据投资费用与元件不可用率的函数关系得到被优化元件不可用率。
[0106]
本技术所达到的有益效果：
[0107]
本技术首先为了提升样本代表性和抽样效率，采用交叉熵重要抽样产生系统失负荷事件，在此基础上建立可靠性指标解析模型。
[0108]
其次，通过构建辅助变量将可靠性指标解析模型融入传统可靠性参数优化问题中，实现了可靠性参数优化问题的解析表达，得到可靠性参数优化解析模型。
[0109]
鉴于所提可靠性参数优化解析模型为混合整数非线性规划问题，本发明进一步将其分解为双层问题进行求解，上层根据目标函数对元件投资费用的解析灵敏度，通过元件重要度排序选择待优化元件；下层的非线性规划问题进一步分解为一个线性规划与一个小型非线性规划，依次求解两个子问题，最终得到被优化元件的投资费用和不可用率。
[0110]
应用rbts、ieee-rts79和ieee-rts96系统测试了本发明所建立模型和求解方法的有效性，并与现有可靠性参数优化方法进行对比，结果表明，本发明所提解析模型及求解方法，不仅可以选择被优化元件，还能确定被优化元件的不可用率及投资费用，在保证计算精度的前提下，极大的提高了可靠性参数优化问题的计算效率。
附图说明
[0111]
图1为本发明方法流程图；
具体实施方式
[0112]
下面结合附图对本技术作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本技术的保护范围。
[0113]
如图1所示，本发明的基于交叉熵重要抽样的元件可靠性参数优化解析方法，包括以下步骤：
[0114]
本发明具体实施时，可靠性指标是指eens指标，可靠性参数是指可用率。
[0115]
步骤1：输入系统的电气参数、初始可靠性参数等，设置系统的元件总个数ne；
[0116]
具体可参见下文的实施例部分，本发明采用电力系统可靠性领域公认的三个标准测试系统：rbts、rts79和rts96来验证本发明的正确性和有效性，系统的电气参数(系统拓扑、负荷水平、机组容量)、初始可靠性参数(机组、线路可用率)、元件总数等均可从网络上获得，本领域技术人员公知。
[0117]
步骤2：采用交叉熵重要抽样产生系统失负荷事件，并进行系统可靠性评估，存储失负荷事件、削负荷量、似然比，并累计所有最优抽样事件的eens，得出系统eens指标，具体包括：
[0118]
步骤2.1：交叉熵重要抽样初始化，包括设置方差系数收敛阈值β0，每次迭代抽样次数n，初始化迭代次数k＝1，设置初始化最优重要抽样概率密度函数；
[0119]
具体实施时，设置初始化最优重要抽样概率密度函数为原始概率密度函数，即：
[0120]
为初始化最优重要抽样概率密度函数，用于预抽样；
[0121]
f(x；p)为系统状态向量x(元件随机故障、分级负荷等)的原始概率密度函数；
[0122]
设重要抽样概率密度函数i(x；q)与原始概率密度函数f(x；p)具有相同的形式，其联合概率密度函数表示为：用于下文求解似然比；
[0123]
其中，q0表示初始抽样概率密度分布的可靠性参数，q表示重要抽样概率密度分布的可靠性参数，qj为元件j的不可用度，x
ij
表示系统状态样本xi中元件j的状态，j＝1,
…
,ne。
[0124]
步骤2.2：预抽样：采用步骤2.1设置的初始化最优重要抽样概率密度函数进行抽样，生成n个系统状态样本；
[0125]
这里并非抽取失负荷事件，事件是否失负荷由之后的可靠性评估过程判断；
[0126]
步骤2.3：对每个样本进行可靠性评估，确定削负荷量h(xi)，i＝1,
…
,n，xi表示第i个样本的系统状态，通过公式获得下一次迭代中各元件的最优不可用度q
k,j
，得到下一次迭代采用的最优重要抽样概率密度函数，并进行下一次迭代；
[0127]
一个样本经过可靠性评估过程后即可获得削负荷量。本发明的重点在于解析模型，可靠性评估就是采用传统的最优削负荷模型进行优化，得到系统状态xi的削负荷量。
[0128]
当前后两次迭代中各元件不可用度变化的最大值小于或等于0.01时，认为迭代过程收敛，输出最后一次(即迭代过程收敛时的当前次)迭代所得各元件的不可用度即为最优不可用度；
[0129]
其中，q
k,j
表示第k次迭代的元件j的最优不可用度的估计值；
[0130]
w(xi；p；q
k-1
)表示第k次迭代的似然比：
[0131][0132]
式中，w(xi；p；q
k-1
)表示第k次迭代的似然比，x
ij
表示系统状态样本xi中元件j的状态，q
k-1,j
为第k-1次迭代中元件j的不可用率；
[0133]
进一步的，预抽样以各元件的原始不可用度作为初始值，此时似然比w(xi；p；q0)＝1。
[0134]
步骤2.4：循环可靠性评估：采用步骤2.3最后一次迭代使用的最优重要抽样概率密度函数进行正式抽样，抽取一个事件x，对抽取的样本进行系统状态分析，计算事件x的失负荷量为h(x)，令n＝n 1；
[0135]
计算由正式抽样产生的全部随机事件对应的eens指标方差系数β，若β＜β0，转步骤2.5，否则，抽取下一个系统状态进行可靠性评估，直至跳出循环；
[0136]
步骤2.4为可靠性评估过程，为一个循环过程，每次循环抽取一个系统状态，首先抽取系统状态x1，判断系统状态x1是否失负荷，然后计算eens的方差系数，如果满足收敛条件，则跳出循环，如果不满足收敛条件，则进行下一次循环，再抽取一个系统状态x2。直至跳
出循环。
[0137]
步骤2.5：存储失负荷事件、削负荷量、似然比：
[0138]
设每个元件有两个状态：正常状态用0表示，故障状态用1表示，则任意失负荷事件对应的元件状态用一个按照元件编号排序后的二进制数表示，记失负荷事件的总数为n
l
；
[0139]
将正式抽样得到的全部失负荷事件对应的二进制状态序列存储在n
l
×
ne维矩阵b中，并将对应的削负荷量、似然函数值分别存储到n1维列向量h与w中；
[0140]
其中，b的每一行对应一个失负荷事件的二进制状态序列。
[0141]
步骤2.6：累计所有最优抽样事件的eens，得出系统eens指标。
[0142]
累计所有最优抽样事件的eens，得出系统eens指标：
[0143]
步骤3：基于步骤2得到的失负荷事件信息，建立基于交叉熵重要抽样的系统可靠性指标解析模型，包括以下步骤：
[0144]
步骤3.1：采用状态枚举法得到系统eens指标解析模型；
[0145]
步骤3.1中，记采用状态枚举法得到的系统eens指标为：
[0146][0147]
假设被优化元件的个数为nh，每个元件只有工作和失效两种状态，且各元件之间的故障相互独立，则被优化元件的状态组合共有个；
[0148]
设代表nh个被优化元件处在第j个组合状态的系统事件集合，则系统失效事件集合ψ被分为m个子集：
[0149]
被优化元件的m个状态用ej(j＝1,2,
…
,m)表示；
[0150]
显然，这m个状态彼此互斥且构成一个完备事件组；
[0151]
则根据任意失负荷事件x对应的被优化元件所处的组合状态，将所述采用状态枚举法得到的系统eens指标公式拆分为m项之和，即：
[0152][0153]
将公式(1)改写为：
[0154][0155]
结合条件概率的定义，将公式(2)写为：
[0156][0157]
其中，p(ej)是被优化元件组合状态ej的发生概率；
[0158]
p(x|ej)是给定被优化元件组合状态为ej时系统失负荷事件x发生的条件概率；
[0159]
最终，将公式(3)简写为：
[0160][0161]
其中，对任意j＝1,2,...,m
[0162][0163]
在系统拓扑结构、电气参数和运行参数不变的前提下，p(x|ej)取决于n
e-nh个未被优化元件的可靠性水平。同时，h(x)的取值只与系统的负荷水平、网络结构等有关。由于本发明假设n
e-nh个其他元件的可靠性水平不变，因此k
1-km应为恒定值，与被优化元件的可靠性参数变化无关。由于假设各元件之间的故障相互独立，p(ej)的计算公式为：
[0164][0165]
φ
ju
和φ
jd
分别是组合状态ej中处于工作和故障状态的被优化元件集合。a对应被优化元件的可用率；
[0166]
将公式(5)和公式(6)代入公式(4)得到采用状态枚举法得到的系统eens指标解析模型，具体为：
[0167][0168]
步骤3.2：在步骤3.1解析模型推导基础上，结合步骤2得到的系统eens指标，得到基于交叉熵重要抽样的系统可靠性指标解析模型。
[0169]
步骤3.2中，在步骤3.1解析模型推导基础上，同样假设被优化的元件的个数为nh，则被优化元件的组合状态一共有个。
[0170]
可根据任意失负荷事件x对应的被优化元件所处的组合状态，将步骤2得到的系统eens指标拆分为m项之和，即：
[0171][0172]
由于公式(8)中没有显含事件的发生概率，故无法从公式中直接提取p(fj)。然而对于模拟法，理论上当抽样样本足够大时，采用状态枚举法和基于交叉熵重要抽样得到的eens指标是相等的，对比公式(4)与公式(8)，得到：
[0173][0174]
则式(9)求得系数kj：
[0175][0176]
通过式(10)可得基于交叉熵重要抽样的系统eens指标关于被优化元件可用率的解析函数，即为基于交叉熵重要抽样的系统可靠性指标解析模型：
[0177][0178]
需要指出的是，上式中的解析模型系数kj需要经过一次完整的可靠性评估计算才能从系统失负荷信息与系统状态信息中得到，并且一旦求得kj，当被优化元件的不可用率发生变化后，无需重复进行系统可靠性评估计算，只需将改变后的不可用率代入解析式即可直接解析求取系统eens指标，极大的提高了计算效率。即式(11)用于式(12)，传统可靠性参数优化问题因为系统eens指标不是元件可靠性参数的解析函数，只能通过一次可靠性评估过程得到系统eens指标，无法像式(12)一样写为类似表达式，更无法像式(12)、(13)一样直接给出优化模型，传统可靠性参数优化均采用智能算法求解，过于费时。
[0179]
步骤4：通过构建辅助变量将系统可靠性指标解析模型融入元件可靠性参数优化问题中，建立可靠性参数优化解析模型；
[0180]
具体实施时，以最小化系统eens指标为优化目标建立可靠性参数优化解析模型：
[0181][0182]
s.t.
[0183][0184][0185][0186]
上述可靠性参数优化解析模型为混合整数非线性优化，在模型中，矩阵b、g和p，向量h、w、b、和τ
all
均为已知量；
[0187]
其中b、h、w是通过执行步骤2.5得到的输出结果。为系统全部元件的初始不可用率，为全部元件不可用率的取值下限，为全部元件的投资费用系数，τ
all
是对应于全部元件的常系数向量，和τ
all
均为ne维行向量。
[0188]
决策变量为：代表优化决策的0-1变量矩阵y和被优化元件i的投资费用ci。在下文的公式描述中，如无特别说明，索引i、j和l的取值范围分别为：其中，为系统被优化元件集合；为系统被优化元件的所有组合状态集合，所有组合状态集合，为系统元件集合。
[0189]
y为ne×
nh维表征系统元件是否被选择优化的二进制变量矩阵，其中ne为系统元件总个数，nh为系统被优化元件个数。若元件l被优化，则矩阵y中第l行中的任意一个元素为1；否则，第l行元素均为0。式(15)是矩阵y自身的逻辑约束。
[0190]
式(13)用于求取本发明所提解析模型系数kj，计算公式如式(13.a)中所示，其中，αj为n
l
维列向量，用于线性相加失负荷事件的失负荷信息，从而避免逻辑运
算。
[0191]
式(13.b)为辅助变量向量αj的下限约束，式(13.b)中构造了nh×ncs
维常数矩阵g、p与n
cs
维常数向量b；
[0192]
g矩阵用于存储nh个被优化元件的n
cs
种组合状态，每个状态可以用一个按照元件编号排序后的二进制数的序列表示，g的每一列对应一个组合状态的二进制数序列。例如，若nh＝2，则矩阵g可表示为：
[0193][0194]
矩阵p的维度与g相同，若g的第i行第j列的元素g
i,j
＝1，则p
i,j
＝1；若g
i,j
＝0，则p
i,j
＝-1。p的第j列记为pj。
[0195]
b为n
cs
维行向量；向量b的第j个元素bj的计算公式如式(13.e)所示；
[0196]
式(13.c)为辅助变量向量αj中每个元素的取值范围，式(13.d)与(14.h)分别为被优化元件组合状态ej的初始概率p
rd
(ej)与元件不可用率变化后的概率p(ej)的计算公式；
[0197]
式(14)表示被优化元件不可用率ui与投资费用ci之间的函数关系，式(14.b)-(14.d)分别表示通过矩阵相乘提取被优化元件对应的投资费用系数c
rd
、不可用率初始值u
rd
和不可用率下限u
min
，其中，，其中，式(14.e)为被优化元件不可用率u(u＝u1,
…
,uh)的上下限约束，由于元件可靠性参数优化只会降低元件不可用率，因此不可用率的上限即为初始值。式(14.e)中τ是个常系数向量，式(14.g)为元件优化费用的总预算约束。
[0198]
由于本发明所提可靠性参数优化解析模型是一个混合整数非线性规划问题，很难直接求解。但因为本发明所提模型是解析的，针对任意被优化元件可以快速求取目标函数对于该元件投资费用的灵敏度。因此本发明采用基于灵敏度的薄弱环节辨识方法选取被优化元件，消除0-1变量，将混合整数非线性规划问题转化为连续非线性规划问题后，通过内点法求解，具体为步骤5-8。
[0199]
步骤5：将可靠性参数优化解析模型分解为上层问题和下层问题；
[0200]
具体实施时，分解得到上层问题为：
[0201][0202]
s.t.
[0203][0204]
[0205]
上层问题用于求取本发明所提解析模型系数kj；
[0206]
式(17)即式(13)；
[0207]
约束式中，αj为n
l
维列向量，用于线性相加失负荷事件的失负荷信息，从而避免逻辑运算。
[0208]
第二条为辅助变量向量αj的下限约束，其中引入了nh×ncs
维常数矩阵g、p与n
cs
维常数向量b。g矩阵用于存储nh个被优化元件的n
cs
种组合状态，每个状态可以用一个按照元件编号排序后的二进制数的序列表示，g的每一列对应一个组合状态的二进制数序列。
[0209]
矩阵p的维度与g相同，若g的第i行第j列的元素g
i,j
＝1，则p
i,j
＝1；若g
i,j
＝0，则p
i,j
＝-1。
[0210]
b为n
cs
维行向量。向量b的第j个元素bj的计算公式如第五条约束所示。
[0211]
第四条约束为被优化元件组合状态ej的初始概率p
rd
(ej)的计算公式。
[0212]
分解得到下层问题为：
[0213][0214]
s.t.
[0215][0216]
式(20)即式(14)；
[0217]
下层问题的约束式表示被优化元件不可用率ui与投资费用ci之间的函数关系，第二至第四条约束分别表示通过矩阵相乘提取被优化元件对应的投资费用系数c
rd
、不可用率初始值u
rd
和不可用率下限u
min
，第六条约束为被优化元件不可用率的上下限约束，由于元件可靠性参数优化只会降低元件不可用率，因此不可用率的上限即为初始值。第五条约束中τ是个常系数向量。第七条约束为元件优化费用的总预算约束；第八条约束为元件不可用率变化后的概率p(ej)的计算公式。
[0218]
步骤6：基于失负荷事件、削负荷量和似然比，针对系统中的所有元件求解上层问题得到解析模型系数，并对系统元件进行重要度排序；
[0219]
具体实施时，根据系统可靠性指标对各元件投资费用的灵敏度对系统元件进行重要度排序，系统可靠性指标对各元件投资费用的灵敏度计算公式为：
[0220][0221]
步骤7：根据元件重要度排序选取被优化元件集合，依次求解上层问题与下层问题，得到被优化元件投资费用及被优化后的不可用率；
[0222]
具体实施时，采用正序排序，从前往后取，取多少个由主观给定。
[0223]
根据步骤6得到基于灵敏度的元件重要度排序后，即可确定被优化元件集合，式(13)—(15)中的优化决策矩阵y也随之确定，降低了模型求解难度，同时因为变量kj仅由优化决策矩阵y唯一决定，不受其他变量影响，消除了模型中的整数变量，将可靠性参数优化模型进一步简化，依次求解上层问题与下层问题，即可得到各元件被优化后的不可用率。
[0224]
进一步的，本发明提出了元件可靠性参数与系统可靠性指标的解析函数，并将该函数嵌入可靠性参数优化问题中，提出了可靠性参数优化解析模型。以往的可靠性参数优化模型由于可靠性评估环节并未解析化，所以只能通过智能算法求解，求解过程中需要多次(甚至1000次)评估系统可靠性，费时费力，本发明将可靠性评估解析化后，可将可靠性参数优化问题直接放入商用优化求解器求解，效率提升几百至1000倍。
[0225]
具体实施时，上层问题的决策变量为kj与αj，为线性规划问题；
[0226]
下层问题的决策变量为被优化元件的投资费用ci，为非线性规划问题，本发明采用内点法进行求解。
[0227]
在确定被优化元件集合后，依次求解上层问题与下层问题即可求得被优化元件的投资费用，依据投资费用与元件不可用率的函数关系得到被优化元件投资费用及被优化后的不可用率。
[0228]
投资费用与元件不可用率的函数关系即式(14.1)：
[0229]
具体实施时，本发明上层模型和下层模型均可采用成熟的商用求解器进行求解。
[0230]
步骤8：输出被优化元件集合及其中各元件的投资费用及被优化后的不可用率。
[0231]
实施例
[0232]
将本发明提出的基于交叉熵重要抽样的元件可靠性参数优化解析模型应用于rbts、ieee-rts79及ieee-rts96系统，具体执行过程如下：
[0233]
第一步，为了提升样本代表性和抽样效率，采用交叉熵重要抽样产生系统失负荷事件，见步骤1—步骤2；
[0234]
第二步，建立系统可靠性指标关于被优化元件可靠性参数的解析模型，见步骤3；
[0235]
第三步，通过构建辅助变量将可靠性指标解析模型融入传统可靠性参数优化问题中，建立可靠性参数优化解析模型，见步骤4；
[0236]
第四步，本发明将混合整数非线性规划问题其分解为双层问题进行求解，上层根据目标函数对元件投资费用的解析灵敏度，通过元件重要度排序选择待优化元件；下层的
非线性规划问题进一步分解为一个线性规划与一个小型非线性规划，依次求解两个子问题，最终得到被优化元件的投资费用和不可用率，见步骤5—步骤8；
[0237]
第五步，为评估本发明的有效性，将利用本发明得到的rbts、ieee-rts79和ieee-rts96系统计算结果与现有的可靠性参数优化方法进行对比，三个测试系统的负荷均采用峰荷数据。算例测试在intel core i5-7500、3.4ghz、8gb的计算机中进行。
[0238]
本发明中被优化元件包括机组和线路，输电线路的c
rd
取10k$/km,机组的c
rd
取10k$/mw，rbts系统和rts79系统的总投资预算均取100,000k$,而rts96取300,000k$。工程实际中，由于元件不可能被优化到完全可靠，假设元件不可用率的下限u
min
为1e-5。
[0239]
实施结果见表1-8，表1和表2为交叉熵重要抽样法相较于标准蒙特卡洛抽样法的计算精确度和计算效率对比表；
[0240]
表3、表4、表5和表6为优化元件数目不同时的rts-79系统优化结果对比表；
[0241]
表7为优化元件数目为8时的rts-96系统优化结果对比表；
[0242]
表8为传统可靠性优化与本发明所提解析模型计算时间对比表。
[0243]
表1两种方法得到的系统eens指标对比
[0244][0245]
表2两种方法的计算效率对比
[0246][0247]
表3被优化元件数目为2时的rbts系统优化结果
[0248][0249][0250]
表4被优化元件数目为4时的rbts系统优化结果
[0251]
[0252]
表5被优化元件数目为4时的rts-79系统优化结果
[0253][0254]
表6被优化元件数目为6时的rts-79系统优化结果
[0255][0256]
表7被优化元件数目为8时的rts-96系统优化结果
[0257][0258]
表8传统可靠性优化与本发明所提解析模型计算时间比较
[0259][0260]
结果表明，本发明所提解析模型及算法，不仅可以选择被优化元件，还能确定被优
化元件的不可用率及投资费用，在保证计算精度的前提下，极大的提高了可靠性参数优化问题的计算效率。
[0261]
本发明所设计的基于交叉熵重要抽样的元件可靠性参数优化解析模型，在求解过程中只需要执行一次基于交叉熵重要抽样的可靠性评估计算，避免了因可靠性参数变化而需要多次进行可靠性评估的弊端。本发明基于交叉熵重要抽样法建立可靠性指标解析模型，并将可靠性指标解析模型嵌入到传统可靠性参数优化问题中，实现可靠性参数优化模型的解析表达，大幅提高了可靠性参数优化问题的求解效率。
[0262]
综上所述，本技术提供了一种可靠性参数优化解析模型。首先为了提升样本代表性和抽样效率，采用交叉熵重要抽样产生系统失负荷事件，在此基础上建立可靠性指标解析函数。其次，通过构建辅助变量将可靠性指标解析模型融入传统可靠性参数优化问题中，实现了可靠性参数优化问题的解析表达。鉴于所提模型为混合整数非线性规划问题，本发明将其分解为双层问题进行求解，上层根据目标函数对元件投资费用的解析灵敏度，通过元件重要度排序选择待优化元件；下层的非线性规划问题进一步分解为一个线性规划与一个小型非线性规划，依次求解两个子问题，最终得到被优化元件的投资费用和不可用率。应用rbts、ieee-rts79和ieee-rts96系统测试了所提模型和算法的有效性，并与现有可靠性参数优化方法进行对比，结果表明，本发明所提解析模型及算法，不仅可以选择被优化元件，还能确定被优化元件的不可用率及投资费用，在保证计算精度的前提下，极大的提高了可靠性参数优化问题的计算效率。
[0263]
本发明申请人结合说明书附图对本发明的实施示例做了详细的说明与描述，但是本领域技术人员应该理解，以上实施示例仅为本发明的优选实施方案，详尽的说明只是为了帮助读者更好地理解本发明精神，而并非对本发明保护范围的限制，相反，任何基于本发明的发明精神所作的任何改进或修饰都应当落在本发明的保护范围之内。