发明名称 一种变压器维修风险决策方法
摘要 一种变压器维修风险决策方法,所述方法使用模糊评价和精确评价相结合的方法构建初始维修决策矩阵,根据正交表的试验次数组成决策方案布点集,计算每个维修方案与正负理想解之间的马氏距离,设定维修决策影响因子,进而构建正负前景价值矩阵,最后计算各维修方案的综合前景值并排序选择最优方案。本发明采用定量和定性相结合的方式,避免了将模糊数精确化过程造成信息的损失,同时引入灰色累积前景理论,提高了决策方法的有效性。
申请公布号 CN105303331A 申请公布日期 2016.02.03
申请号 CN201510829987.2 申请日期 2015.11.25
申请人 国网江西省电力科学研究院;国家电网公司 发明人 曾伟;孙旻;范瑞祥;曹蓓
分类号 G06Q10/06(2012.01)I;G06Q50/06(2012.01)I 主分类号 G06Q10/06(2012.01)I
代理机构 南昌市平凡知识产权代理事务所 36122 代理人 姚伯川
主权项 一种变压器维修风险决策方法,其特征在于,所述方法使用模糊评价和精确评价相结合的方法构建初始维修决策矩阵,根据正交表的试验次数组成决策方案布点集,计算每个维修方案布点集与正负理想解之间的马氏距离,进而构建正负前景价值矩阵,最后计算各维修方案的综合前景值并排序选择最优方案;所述方法包括以下步骤:(1)构建变压器维修决策指标属性集合I={I<sub>1</sub>,I<sub>2</sub>,I<sub>3</sub>,I<sub>4</sub>,I<sub>5</sub>,I<sub>6</sub>},I<sub>1</sub>为维修需要的技术水平,I<sub>2</sub>为维修方案可能带来的风险,I<sub>3</sub>为维修方案对用户负荷产生的影响,I<sub>4</sub>为故障对电网安全的威胁,I<sub>5</sub>为维修方案的综合维修费用,I<sub>6</sub>为维修方案的综合效益;每一个指标的评价值用E(I<sub>i</sub>),i=1,2,3,4,5,6来表示;对I<sub>2</sub>,I<sub>3</sub>,I<sub>4</sub>,I<sub>5</sub>采用精确评估方法;E(I<sub>1</sub>)表示不同的维修方案需要不同级别的维修技术,E(I<sub>2</sub>)为不同维修方案在实施中可能发生的风险,E(I<sub>3</sub>)为不同方案维修后对变电站的影响,E(I<sub>4</sub>)表示在不同的维修方案下,设备的故障引起电力系统损失的期望值,E(I<sub>5i</sub>(t))表示不同维修方案的总体维修费用的期望值,E(I<sub>6</sub>)表示维修附加效益值;(2)设有n个维修决策方案A={A<sub>1</sub>,A<sub>2</sub>,...,A<sub>n</sub>},衡量每个方案优劣的指标属性有m个I={I<sub>1</sub>,I<sub>2</sub>,...,I<sub>m</sub>},每一个备选方案在指标下对应的区间数指标值为<img file="FDA0000857515870000011.GIF" wi="247" he="79" /><img file="FDA0000857515870000012.GIF" wi="222" he="77" />1≤i≤m,1≤j≤n;由s<sub>ij</sub>构成初始区间决策矩阵S=(s<sub>ij</sub>)<sub>m×n</sub>;对初始区间决策矩阵进行规范化处理,得到规范化区间矩阵<img file="FDA0000857515870000013.GIF" wi="238" he="101" /><img file="FDA0000857515870000014.GIF" wi="245" he="85" />根据指标属性进行规范化操作;对决策矩阵<img file="FDA0000857515870000019.GIF" wi="34" he="70" />进行加权处理,<img file="FDA0000857515870000015.GIF" wi="795" he="86" />其中ω<sub>j</sub>为各指标的权重;(3)寻找正负理想解,正理想解X<sup>+</sup>通过下式计算得到:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msup><mi>X</mi><mo>+</mo></msup><mo>=</mo><mrow><mo>(</mo><mo>&lsqb;</mo><msubsup><mi>x</mi><mrow><mi>i</mi><mn>1</mn></mrow><mrow><mi>l</mi><mo>*</mo><mo>+</mo></mrow></msubsup><mo>,</mo><msubsup><mi>x</mi><mrow><mi>i</mi><mn>1</mn></mrow><mrow><mi>u</mi><mo>*</mo><mo>+</mo></mrow></msubsup><mo>&rsqb;</mo><mo>,</mo><mo>...</mo><mo>,</mo><mo>&lsqb;</mo><msubsup><mi>x</mi><mrow><mi>i</mi><mi>m</mi></mrow><mrow><mi>l</mi><mo>*</mo><mo>+</mo></mrow></msubsup><mo>,</mo><msubsup><mi>x</mi><mrow><mi>i</mi><mi>m</mi></mrow><mrow><mi>u</mi><mo>*</mo><mo>+</mo></mrow></msubsup><mo>&rsqb;</mo><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mo>&lsqb;</mo><munder><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mi>i</mi></munder><msubsup><mi>x</mi><mrow><mi>i</mi><mn>1</mn></mrow><mi>l</mi></msubsup><mo>,</mo><munder><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mi>i</mi></munder><msubsup><mi>x</mi><mrow><mi>i</mi><mn>1</mn></mrow><mi>u</mi></msubsup><mo>&rsqb;</mo><mo>,</mo><mo>...</mo><mo>,</mo><mo>&lsqb;</mo><munder><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mi>i</mi></munder><msubsup><mi>x</mi><mrow><mi>i</mi><mi>m</mi></mrow><mi>l</mi></msubsup><mo>,</mo><munder><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mi>i</mi></munder><msubsup><mi>x</mi><mrow><mi>i</mi><mi>m</mi></mrow><mi>u</mi></msubsup><mo>&rsqb;</mo><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000857515870000016.GIF" wi="1398" he="118" /></maths>其中,<img file="FDA0000857515870000017.GIF" wi="155" he="94" />表示取矩阵列向量最大值;<img file="FDA0000857515870000018.GIF" wi="198" he="77" />表示正理想解的第1个区间数,<img file="FDA0000857515870000021.GIF" wi="78" he="74" />为较小值,<img file="FDA0000857515870000022.GIF" wi="85" he="69" />为较大值;<img file="FDA0000857515870000023.GIF" wi="197" he="71" />表示正理想解的第m个区间数,<img file="FDA0000857515870000024.GIF" wi="77" he="71" />为较小值,<img file="FDA0000857515870000025.GIF" wi="85" he="70" />为较大值;<img file="FDA0000857515870000026.GIF" wi="154" he="71" />分别表示决策矩阵<img file="FDA0000857515870000027.GIF" wi="39" he="63" />中第1列区间数的较小值和较大值;<img file="FDA0000857515870000028.GIF" wi="181" he="71" />分别表示决策矩阵<img file="FDA0000857515870000029.GIF" wi="39" he="62" />中第m列区间数的较小值和较大值;负理想解X<sub>‑</sub>通过下式计算得到:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msup><mi>X</mi><mo>-</mo></msup><mo>=</mo><mrow><mo>(</mo><mo>&lsqb;</mo><msubsup><mi>x</mi><mrow><mi>i</mi><mn>1</mn></mrow><mrow><mi>l</mi><mo>*</mo><mo>-</mo></mrow></msubsup><mo>,</mo><msubsup><mi>x</mi><mrow><mi>i</mi><mn>1</mn></mrow><mrow><mi>u</mi><mo>*</mo><mo>-</mo></mrow></msubsup><mo>&rsqb;</mo><mo>,</mo><mo>...</mo><mo>,</mo><mo>&lsqb;</mo><msubsup><mi>x</mi><mrow><mi>i</mi><mi>m</mi></mrow><mrow><mi>l</mi><mo>*</mo><mo>-</mo></mrow></msubsup><mo>,</mo><msubsup><mi>x</mi><mrow><mi>i</mi><mi>m</mi></mrow><mrow><mi>u</mi><mo>*</mo><mo>-</mo></mrow></msubsup><mo>&rsqb;</mo><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mo>&lsqb;</mo><munder><mi>min</mi><mi>i</mi></munder><msubsup><mi>x</mi><mrow><mi>i</mi><mn>1</mn></mrow><mi>l</mi></msubsup><mo>,</mo><munder><mi>min</mi><mi>i</mi></munder><msubsup><mi>x</mi><mrow><mi>i</mi><mn>1</mn></mrow><mi>u</mi></msubsup><mo>&rsqb;</mo><mo>,</mo><mo>...</mo><mo>,</mo><mo>&lsqb;</mo><munder><mi>min</mi><mi>i</mi></munder><msubsup><mi>x</mi><mrow><mi>i</mi><mi>m</mi></mrow><mi>l</mi></msubsup><mo>,</mo><munder><mi>min</mi><mi>i</mi></munder><msubsup><mi>x</mi><mrow><mi>i</mi><mi>m</mi></mrow><mi>u</mi></msubsup><mo>&rsqb;</mo><mo>)</mo></mrow></mrow>]]></math><img file="FDA00008575158700000210.GIF" wi="1373" he="111" /></maths>其中,<img file="FDA00008575158700000211.GIF" wi="150" he="89" />表示取矩阵列向量最小值;<img file="FDA00008575158700000212.GIF" wi="198" he="70" />表示负理想解的第1个区间数,<img file="FDA00008575158700000213.GIF" wi="76" he="70" />为较小值,<img file="FDA00008575158700000214.GIF" wi="85" he="69" />为较大值;<img file="FDA00008575158700000215.GIF" wi="198" he="70" />表示负理想解的第m个区间数,<img file="FDA00008575158700000216.GIF" wi="78" he="70" />为较小值,<img file="FDA00008575158700000217.GIF" wi="85" he="70" />为较大值;<img file="FDA00008575158700000218.GIF" wi="150" he="71" />分别表示决策矩阵<img file="FDA00008575158700000219.GIF" wi="39" he="70" />中第1列区间数的较小值和较大值;<img file="FDA00008575158700000220.GIF" wi="174" he="71" />分别表示决策矩阵<img file="FDA00008575158700000221.GIF" wi="40" he="61" />中第m列区间数的较小值和较大值;(4)设计正交试验表,结合步骤(2)每个方案的区间数变量<img file="FDA00008575158700000222.GIF" wi="173" he="79" />水平1表示选择<img file="FDA00008575158700000223.GIF" wi="79" he="78" />水平2表示选择<img file="FDA00008575158700000224.GIF" wi="78" he="71" />得到每个维修方案A<sub>j</sub>的布点集T<sub>j</sub>;根据正交试验表,结合步骤(3)的正负理想解,得到正负理想解的布点集<img file="FDA00008575158700000225.GIF" wi="72" he="70" />和<img file="FDA00008575158700000226.GIF" wi="93" he="62" />(5)分别计算步骤(4)的每个维修方案布点集T<sub>j</sub>到正负理想解<img file="FDA00008575158700000227.GIF" wi="70" he="70" />和<img file="FDA00008575158700000228.GIF" wi="63" he="62" />的马氏距离,用ZMD<sub>j×k</sub>和FMD<sub>j×k</sub>表示,其中1≤k≤10,1≤j≤n;(6)按如下公式计算步骤(5)中马氏距离的灰色关联系数,分别用<img file="FDA00008575158700000229.GIF" wi="78" he="67" />和<img file="FDA00008575158700000230.GIF" wi="77" he="69" />表示:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>&xi;</mi><mrow><mi>j</mi><mo>&times;</mo><mi>k</mi></mrow><mo>-</mo></msubsup><mo>=</mo><mfrac><mrow><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mi>j</mi></munder><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mi>k</mi></munder><mo>|</mo><msub><mi>FMD</mi><mrow><mi>j</mi><mo>&times;</mo><mi>k</mi></mrow></msub><mo>|</mo><mo>+</mo><mi>&rho;</mi><munder><mi>max</mi><mi>j</mi></munder><munder><mi>max</mi><mi>k</mi></munder><mo>|</mo><msub><mi>FMD</mi><mrow><mi>j</mi><mo>&times;</mo><mi>k</mi></mrow></msub><mo>|</mo></mrow><mrow><mo>|</mo><msub><mi>FMD</mi><mrow><mi>j</mi><mo>&times;</mo><mi>k</mi></mrow></msub><mo>|</mo><mo>+</mo><mi>&rho;</mi><munder><mi>max</mi><mi>j</mi></munder><munder><mi>max</mi><mi>k</mi></munder><mo>|</mo><msub><mi>FMD</mi><mrow><mi>j</mi><mo>&times;</mo><mi>k</mi></mrow></msub><mo>|</mo></mrow></mfrac></mrow>]]></math><img file="FDA00008575158700000231.GIF" wi="835" he="182" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>&xi;</mi><mrow><mi>j</mi><mo>&times;</mo><mi>k</mi></mrow><mo>+</mo></msubsup><mo>=</mo><mfrac><mrow><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mi>j</mi></munder><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mi>k</mi></munder><mo>|</mo><msub><mi>ZMD</mi><mrow><mi>j</mi><mo>&times;</mo><mi>k</mi></mrow></msub><mo>|</mo><mo>+</mo><mi>&rho;</mi><munder><mi>max</mi><mi>j</mi></munder><munder><mi>max</mi><mi>k</mi></munder><mo>|</mo><msub><mi>ZMD</mi><mrow><mi>j</mi><mo>&times;</mo><mi>k</mi></mrow></msub><mo>|</mo></mrow><mrow><mo>|</mo><msub><mi>ZMD</mi><mrow><mi>j</mi><mo>&times;</mo><mi>k</mi></mrow></msub><mo>|</mo><mo>+</mo><mi>&rho;</mi><munder><mi>max</mi><mi>j</mi></munder><munder><mi>max</mi><mi>k</mi></munder><mo>|</mo><msub><mi>ZMD</mi><mrow><mi>j</mi><mo>&times;</mo><mi>k</mi></mrow></msub><mo>|</mo></mrow></mfrac></mrow>]]></math><img file="FDA00008575158700000232.GIF" wi="822" he="181" /></maths>其中,ρ为分辨系数,一般取0.5,min(·)为取最小值操作;(7)根据步骤(6)的灰色关联系数,得到各维修方案的前景价值矩阵,分别用<img file="FDA00008575158700000233.GIF" wi="68" he="68" />和<img file="FDA00008575158700000234.GIF" wi="74" he="63" />表示,其矩阵元素分别按如下公式计算得到:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>v</mi><mrow><mi>j</mi><mo>&times;</mo><mi>k</mi></mrow><mo>-</mo></msubsup><mo>=</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msubsup><mi>&xi;</mi><mrow><mi>j</mi><mo>&times;</mo><mi>k</mi></mrow><mo>-</mo></msubsup><mo>)</mo></mrow><mi>&alpha;</mi></msup></mrow>]]></math><img file="FDA0000857515870000031.GIF" wi="301" he="86" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msubsup><mi>v</mi><mrow><mi>j</mi><mo>&times;</mo><mi>k</mi></mrow><mo>+</mo></msubsup><mo>=</mo><mo>-</mo><mi>&theta;</mi><msup><mrow><mo>(</mo><mo>-</mo><mo>(</mo><mrow><msubsup><mi>&xi;</mi><mrow><mi>j</mi><mo>&times;</mo><mi>k</mi></mrow><mo>+</mo></msubsup><mo>-</mo><mn>1</mn></mrow><mo>)</mo><mo>)</mo></mrow><mi>&beta;</mi></msup></mrow>]]></math><img file="FDA0000857515870000032.GIF" wi="421" he="102" /></maths>其中,参数α,β一般取值0.88,θ取值2.25;(8)根据步骤(7)的前景价值矩阵,按如下公式计算前景权重函数值:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>Q</mi><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msup><mi>Q</mi><mi>+</mi></msup><mo>=</mo><mfrac><msup><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>q</mi></msub><mo>)</mo></mrow><msup><mi>&gamma;</mi><mo>+</mo></msup></msup><msup><mrow><mo>(</mo><mrow><msup><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>q</mi></msub><mo>)</mo></mrow><msup><mi>&gamma;</mi><mi>+</mi></msup></msup><mi>+(1-</mi><msub><mi>&omega;</mi><mi>q</mi></msub><msup><mi>)</mi><msup><mi>&gamma;</mi><mo>+</mo></msup></msup></mrow><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><msup><mi>&gamma;</mi><mo>+</mo></msup></mrow></msup></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><msup><mi>Q</mi><mo>-</mo></msup><mo>=</mo><mfrac><msup><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>q</mi></msub><mo>)</mo></mrow><msup><mi>&gamma;</mi><mo>-</mo></msup></msup><msup><mrow><mo>(</mo><mrow><msup><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>q</mi></msub><mo>)</mo></mrow><msup><mi>&gamma;</mi><mo>-</mo></msup></msup><mo>+</mo><msup><mrow><mi>(1-</mi><msub><mi>&omega;</mi><mi>q</mi></msub><mi>)</mi></mrow><msup><mi>&gamma;</mi><mo>-</mo></msup></msup></mrow><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><msup><mi>&gamma;</mi><mo>-</mo></msup></mrow></msup></mfrac></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000857515870000033.GIF" wi="645" he="463" /></maths>一般取参数γ<sup>+</sup>=0.61,γ<sup>‑</sup>=0.69,权重ω<sub>q</sub>=1/10;其中,Q<sup>+</sup>表示正理想解的前景权重函数值,Q<sup>‑</sup>为负理想解的前景权重函数值;(9)结合步骤(7)和步骤(8)计算各方案的综合前景值:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>V</mi><mi>j</mi></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mn>10</mn></munderover><mrow><msubsup><mi>v</mi><mrow><mi>j</mi><mo>&times;</mo><mi>k</mi></mrow><mo>+</mo></msubsup><msup><mi>Q</mi><mo>+</mo></msup></mrow><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mn>10</mn></munderover><mrow><msubsup><mi>v</mi><mrow><mi>j</mi><mo>&times;</mo><mi>k</mi></mrow><mo>-</mo></msubsup><msup><mi>Q</mi><mo>-</mo></msup></mrow></mrow>]]></math><img file="FDA0000857515870000034.GIF" wi="532" he="134" /></maths>在计算出各维修方案综合最优前景值后,对方案前景值进行排序,最大值对应的方案确定为最优方案。
地址 330096 江西省南昌市民营科技园民强路88号