发明名称 一种变电站接地网腐蚀速率的模糊预测方法
摘要 本发明公开了一种变电站接地网腐蚀速率的模糊预测方法,包括依据土壤理化指标与变电站接地网腐蚀相关关系对隶属度的确定、对不同土壤理化指标对接地网腐蚀贡献度的确定以及不同地域变电站接地网间贴近度的计算。本发明采用模糊预测模型以解决接地网腐蚀速率预测中存在的涉及因素较多,关系错综复杂,以及所涉及因素中的不确定性和模糊性,同时避免了因小样本引起的模型预测精度差,对地域上邻近或土壤理化指标相近的接地网站点,在工程上能够快速而有效的预测其变电站接地网腐蚀速率。
申请公布号 CN103163064B 申请公布日期 2015.01.28
申请号 CN201310082265.6 申请日期 2013.03.14
申请人 陕西电力科学研究院;国家电网公司 发明人 闫爱军;杜京义;冯拉俊;李娜;唐小华;刘磊
分类号 G01N17/00(2006.01)I 主分类号 G01N17/00(2006.01)I
代理机构 西安通大专利代理有限责任公司 61200 代理人 蔡和平
主权项 一种变电站接地网腐蚀速率的模糊预测方法,其特征在于,包括以下步骤:步骤一、获取m个不同地域变电站接地网腐蚀速率及其附近n种土壤理化指标;步骤二、确定各土壤理化指标与腐蚀速率相关系数;步骤三、根据相关系数确定隶属度和贡献度;步骤四、计算m个变电站接地网土壤与待预测变电站接地网土壤之间贴近度;步骤五、对计算贴近度结果进行由大到小排序;步骤六、选择贴近度最大的前三变电站接地网腐蚀数据作为参考,对待预测变电站接地网腐蚀速率进行预测;步骤二具体包括:通过公式2计算n个土壤理化指标与腐蚀速率之间的相关系数R;R=[r<sub>1n+1</sub>,r<sub>2n+1</sub>,…,r<sub>nn+1</sub>]其中r<sub>yn+1</sub>表示第y种土壤理化指标与腐蚀速率之间的相关系数;y=1…n;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>r</mi><mi>ij</mi></msub><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mrow><mo>(</mo><msub><mi>x</mi><mi>ik</mi></msub><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>x</mi><mi>jk</mi></msub><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>j</mi></msub><mo>)</mo></mrow></mrow><mrow><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>ik</mi></msub><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>jk</mi></msub><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>j</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000554034850000011.GIF" wi="1083" he="241" /></maths>其中:r<sub>ij</sub>表示第i种特征因素与第j种特征因素的相关系数;x<sub>ik</sub>为第k个地域的第i种特征因素;i≤n时的特征因素为土壤理化指标,i=n+1时的特征因素为接地网腐蚀速率;x<sub>jk</sub>为第k个地域的第j种特征因素;j≤n时的特征因素为土壤理化指标,j=n+1时的特征因素为接地网腐蚀速率;<img file="FDA0000554034850000012.GIF" wi="49" he="61" />为第i种土壤理化指标的平均值,<img file="FDA0000554034850000013.GIF" wi="56" he="67" />为第j种土壤理化指标的平均值,k=1…m,i=1…n+1,j=1…n+1;步骤三具体包括:当第y种土壤理化指标与腐蚀速率相关系数r<sub>yn+1</sub>>0时,通过公式3,当r<sub>yn+1</sub><0时通过公式4将土壤理化指标特征因素转换成[0,1]区间上的数值s<sub>yo</sub>,即为隶属度,m个不同地域的n个土壤理化指标特征因素得到隶属度矩阵S:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>S</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>T</mi><mn>1</mn></msub></mtd><mtd><msub><mi>T</mi><mn>2</mn></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>T</mi><mi>o</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>s</mi><mn>11</mn></msub></mtd><mtd><msub><mi>s</mi><mn>12</mn></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>s</mi><mrow><mn>1</mn><mi>o</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>s</mi><mn>21</mn></msub></mtd><mtd><msub><mi>s</mi><mn>22</mn></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>s</mi><mrow><mn>2</mn><mi>o</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>s</mi><mrow><mi>y</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>s</mi><mrow><mi>y</mi><mn>2</mn></mrow></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>s</mi><mi>yo</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000554034850000021.GIF" wi="1401" he="313" /></maths>公式3、公式4分别为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>s</mi><mi>yo</mi></msub><mo>=</mo><mi>u</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>yo</mi></msub><mo>,</mo><mi>a</mi><mo>,</mo><mi>b</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mi>x</mi><mo>&le;</mo><mi>a</mi></mtd></mtr><mtr><mtd><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>g</mi><mo>[</mo><mrow><mo>(</mo><mfrac><mrow><mi>a</mi><mo>+</mo><mi>b</mi></mrow><mn>2</mn></mfrac><mo>-</mo><msub><mi>x</mi><mi>yo</mi></msub><mo>)</mo></mrow><mfrac><mn>1</mn><mrow><mrow><mo>(</mo><mi>a</mi><mo>-</mo><msub><mi>x</mi><mi>yo</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><mi>b</mi><mo>-</mo><msub><mi>x</mi><mi>yo</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>]</mo></mrow></msup></mrow></mfrac></mtd><mtd><mi>a</mi><mo>&lt;</mo><mi>x</mi><mo>&lt;</mo><mi>b</mi></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mi>x</mi><mo>&GreaterEqual;</mo><mi>b</mi></mtd></mtr></mtable></mfenced><mo>,</mo><msub><mi>r</mi><mrow><mi>yn</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>></mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000554034850000022.GIF" wi="1464" he="372" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>s</mi><mi>yo</mi></msub><mo>=</mo><mi>u</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>yo</mi></msub><mo>,</mo><mi>a</mi><mo>,</mo><mi>b</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mi>x</mi><mo>&le;</mo><mi>a</mi></mtd></mtr><mtr><mtd><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mi>g</mi><mo>[</mo><mrow><mo>(</mo><mfrac><mrow><mi>a</mi><mo>+</mo><mi>b</mi></mrow><mn>2</mn></mfrac><mo>-</mo><msub><mi>x</mi><mi>yo</mi></msub><mo>)</mo></mrow><mfrac><mn>1</mn><mrow><mrow><mo>(</mo><mi>a</mi><mo>-</mo><msub><mi>x</mi><mi>yo</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><mi>b</mi><mo>-</mo><msub><mi>x</mi><mi>yo</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>]</mo></mrow></msup></mrow></mfrac></mtd><mtd><mi>a</mi><mo>&lt;</mo><mi>x</mi><mo>&lt;</mo><mi>b</mi></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>x</mi><mo>&GreaterEqual;</mo><mi>b</mi></mtd></mtr></mtable></mfenced><mo>,</mo><msub><mi>r</mi><mrow><mi>yn</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>&lt;</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000554034850000023.GIF" wi="1504" he="381" /></maths>其中x<sub>yo</sub>为第o个地域土壤样品的第y个土壤理化指标,y=1…n,o=1…m,g为1‑20之间的正整数,a为第y个土壤理化指标的最小值,b为第y个土壤理化指标的最大值,s<sub>yo</sub>为第o个地域的第y个土壤理化指标的隶属度;根据第y种土壤理化指标与腐蚀速率之间的相关系数r<sub>yn+1</sub>的绝对值的大小,对土壤理化指标特征因素对接地网腐蚀速率的贡献度划分;令|r<sub>yn+1</sub>|≥0.5贡献度为1,0.4≤|r<sub>yn+1</sub>|≤0.49贡献度为0.9,0.3≤|r<sub>yn+1</sub>|≤0.39贡献度为0.7.0.2≤|r<sub>yn+1</sub>|≤0.29贡献度为0.5,0≤|r<sub>yn+1</sub>|≤0.19贡献度为0.3;得到n个贡献度大小为:w<sub>1</sub>,w<sub>2</sub>,…,w<sub>n</sub>;步骤四具体包括:对待预测变电站接地网附近土壤理化指标检测,并确定其各个指标隶属度s<sub>h测</sub>;检测的待预测变电站接地网附近土壤理化指标的种类与步骤一中检测的每个变电站接地网附近土壤理化指标的种类相同;通过公式6计算待预测变电站接地网处土壤与步骤一中m个不同地域变电站接地网附近土壤的贴进度:<img file="FDA0000554034850000031.GIF" wi="1259" he="128" />其中<img file="FDA0000554034850000032.GIF" wi="1296" he="99" />p=1…m,h=1…n;步骤六具体包括:将选出的三个与待预测变电站接地网处土壤贴进度最高的变电站接地网处土壤与待预测变电站接地网处土壤的贴进度,以及三个选出的变电站接地网腐蚀速率带入公式7中,计算出待预测变电站接地网的腐蚀速率;<img file="FDA0000554034850000033.GIF" wi="1833" he="129" />其中,λ为系数,取值范围为0.8‑1.6。
地址 710054 陕西省西安市友谊东路308号