CN109408929A - 确定高压直流电阻标准分压器均压环的优化设计方法 - Google Patents
确定高压直流电阻标准分压器均压环的优化设计方法 Download PDFInfo
- Publication number
- CN109408929A CN109408929A CN201811195074.XA CN201811195074A CN109408929A CN 109408929 A CN109408929 A CN 109408929A CN 201811195074 A CN201811195074 A CN 201811195074A CN 109408929 A CN109408929 A CN 109408929A
- Authority
- CN
- China
- Prior art keywords
- design variable
- value
- grading ring
- direct current
- design
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本申请提供一种确定高压直流电阻标准分压器均压环的优化设计方法，具体步骤包括：S1、列出主、辅均压环的考察指标，确定各设计变量间相互影响作用表，确定试验阶梯数据方案；S2、建立高压直流电阻标准分压器仿真模型；S3、将S1中的数据代入到S2的仿真模型中，得出电场数据并汇总；S4、根据S3中的数据计算各设计变量试验总和、各设计变量试验总和平均值、极差；S5、通过S4中极差的数值确定各设计变量间相互影响作用的大小；S6、结合S4和S5的数据确定各设计变量的取值，缩小范围；S7、在S6得出的高压直流电阻标准分压器各均压环尺寸和位置的取值，缩小实验阶梯数据范围，迭代重复S1至S7，直到满足工程要求为止。该方法能节省实验时间并获得更优解。
Description
技术领域
本发明属于电力技术领域，特别涉及一种确定高压直流电阻标准分压器均压环的优化设 计方法。
背景技术
直流高电压量值传递和量值溯源的基础是直流电压比例标准的建立，而高精度直流电阻 标准分压器可作为直流电压比例标准。目前，大多数国家仅针对300kV及以下电压等级开展 直流分压器分压比的量值溯源研究，中国和澳大利亚等国家开展了在较高电压等级下的研究 工作。我国于2015年通过了800kV直流电压比例标准的建标考核，大幅提升了电压等级和 比例量值的不确定度水平，电压等级从500kV提升为800kV，比例量值不确定度从2×10-4提 升到5×10-5。在开展800kV以上直流电压互感器现场校准实验时，确定其分压比的方法是在 800kV下用分节标定，这可能导致准确度降低。目前我国特高压直流输电工程的最高电压等 级已达1100kV，如果不在直流电阻分压器周围采取措施均匀电场分布，那么直流电阻分压器 的电压等级越高，其周围电场分布将愈加不均匀，额定电压下分压器电阻的局部发热量越大， 分压器内部的温升越高，分压器中电阻元件的阻值变化越大，分压比的变化量越大，导致分 压器的准确度降低，分压器内部温升过大时甚至可能导致电阻元件的热击穿进而损坏。此外， 随着电压等级的提高，分压器的绝缘难度随之增加，而且电压等级的提高会使分压器的电晕 电流和泄漏电流增大，导致分压器的准确度降低，而均压环能有效改善分压器周围的电场分 布。
发明内容
为了改进由于我国仅通过800kV直流电压比例标准的建标考核，导致800kV以上的高 压直流电阻标准分压器周围均压环尺寸和位置需要重新考虑，而传统方法确定该尺寸和位置 效率低，重复工作量大的问题，本申请提出了一种确定高压直流电阻标准分压器均压环的优 化设计方法，主要是基于软件搭建高压直流电阻标准分压器仿真模型，结合数学分析方法迭 代逼近均压环尺寸和位置的最优值；以此方法计算均压环尺寸和位置的最优值可有效提高计 算效率，数学分析法等可以减小试验次数、减少计算时间，多次迭代可以使得结果快速收敛 于最优值附近，符合工程要求。
本发明是这样实现的：
一种确定高压直流电阻标准分压器均压环的优化设计方法，其具体步骤如下：
S1、确定主均压环和辅均压环的设计变量，对各设计变量多次取值，确定包括误差对比 列在内的设计变量相互影响作用表，确定阶梯数据计算方案；
S2、建立高压直流电阻标准分压器仿真模型；
S3、将S1中各设计变量的取值数据逐一代入到S2中建立的高压直流电阻标准分压器仿 真模型中，得出高压直流电阻标准分压器的周围电场数据yi，汇总至阶梯数据计算方案中；
S4、根据S3中的数据计算各设计变量的指标，具体包括：
S41、计算第j列上第i个设计变量取值的计算结果总和Kij；
S42、计算总和Kij的平均值其中t为第j列上设计变量号i出现的次数；
S43、计算第j列极差Rj，
S5、通过S4中极差的数值，计算各设计变量之间的效应，具体包括以下步骤：
S51、计算各设计变量之间的联合效应[ab]ij；
S52、计算各设计变量之间的交互效应(ab)ij；
S53、判断各设计变量对周围电场数据yi的影响；
S6、结合S4和S5的数据确定各设计变量的取值；
S7、在S6得出的高压直流电阻标准分压器各均压环尺寸和位置的取值后，缩小各设计 变量之间差值，迭代重复S1至S7，使分压器周围最大场强小于等于分压器的起晕场强。
优选地，步骤S1中设计变量具体包括：辅均压环的内环半径A、辅均压环圆心到对称 轴的距离B、辅均压环圆心到地面的距离C、主均压环的内环半径D、主均压环圆心到对称轴的距离E、主均压环圆心到地面的距离F，各设计变量的取值Ai、Bi、Ci、Di、Ei、Fi，各 设计变量取r个值，i＝1～r，各设计变量每次取值的差值为ΔQA＝|A1-A2|、ΔQB＝|B1-B2|、 ΔQC＝|C1-C2|、ΔQD＝|D1-D2|、ΔQE＝|E1-E2|、ΔQF＝|F1-F2|。
优选地，步骤S51和S52中计算的具体步骤为：
计算两个设计变量之间的联合效应[ab]ij，[ab]ij＝μij-μ，其中μij为Ai与Bj的组合下的 均值，总均值为μ，该数值可反应Ai与Bj的组合对计算结果的总效应；
计算两个设计变量不同i值之间的交互效应(ab)ij，(ab)ij＝[ab]ij-ai-bj，该数值可 反应Ai与Bj的组合对计算结果的交互效应。
优选地，步骤S53中的具体步骤为：
计算结果的总偏差平方和ST，其中n＝每列同一设计变量的计算次数 与设计变量的个数相乘，T为数据总和，
计算第j列设计变量的偏差平方和Sj，若设计变量A安排在表中第 j列上，则SA＝Sj；
计算误差对比列的偏差平方和Se，Se为所有误差对比列的偏差平方和之和；
计算总偏差平方和ST的度fT、Sj的度fj、Se的度fe，fT＝n-1，fj＝r-1，Fe为fT与各 设计变量的度之和；
计算检验第j列设计变量对试验结果是否有显著影响的统计量值Fj，若由观察值计算出的Fj≥F1-α(fj，fj)，则认为该列设计变量对试验结果影响显著，否则，认为效果不显著，其中α为给定的显著性水平。
优选地，步骤S6中的具体步骤为：
取设计变量A和设计变量B，计算设计变量A、B间的相互作用A×B，若设计变量A、 B以及A×B对周围电场数据yi的影响均大于系统误差对周围电场数据yi的影响，则按照工程要求选择[ab]ij中最大或最小者；
若设计变量B与设计变量A×B对周围电场数据yi的影响均大于系统误差对周围电场数 据yi的影响，设计变量A对周围电场数据yi的影响小于系统误差对周围电场数据yi的影响，则 按照工程要求选择(ab)ij及bj中最大或最小者。
优选地，步骤S7中的具体步骤为：
选取新一轮迭代的各设计变量间差值为ΔQA′、ΔQB′、ΔQC′、ΔQD′、ΔQE′、ΔQF′，确定 新一轮迭代各设计变量的值为ΔQ′A、ΔQ′B、ΔQ′C、ΔQ′D、ΔQ′E、ΔQ′F，要求ΔQA′<ΔQA、 ΔQB′<ΔQB、ΔQC′<ΔQC、ΔQD′<ΔQD、ΔQE′<ΔQE、ΔQF′<ΔQF，A′k＝Ak±ΔQ′A、 B′k＝Bk±ΔQ′B、Ck′＝Ck±ΔQ′C、D′k＝Dk±ΔQ′D、E′k＝Ek±ΔQ′E、Fk′＝Fk±ΔQ′F，k为迭 代次数；
把新一轮迭代值继续带入S1中的设计变量相互影响作用表中，形成修正后设计变量相互 影响作用表。
与现有技术相比，本发明具有以下有益效果：
1、结合数学分析方法可以更加科学且快速地确定主均压环和辅均压环的半径和位置的 取值范围，大幅度减少传统枚举法需要的计算时间，提升效率；
2、多次迭代可以更加准确地收敛至最优值附近，突破一次计算无法达到工程所需的精 确值的瓶颈，改善分压器周围的电场分布，提高分压器的准确度。
附图说明
图1为本发明确定高压直流电阻标准分压器均压环的优化设计方法的流程图；
图2为高压直流电阻标准分压器仿真模型。
附图标记：分压器1；辅均压环2；主均压环3。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整的描 述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的 实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属 于本发明保护的范围。
本发明提供一种确定高压直流电阻标准分压器均压环的优化设计方法，具体步骤包括：
S1、列出1100kV高压直流电阻标准分压器的主均压环和辅均压环需要考察的六个设计 变量A、B、C、D、E、F，分别为辅均压环的内环半径、辅均压环圆心到对称轴的距离、辅均压环圆心到地面的距离、主均压环的内环半径、主均压环圆心到对称轴的距离、主均压环圆心到地面的距离，各设计变量的取值Ai、Bi、Ci、Di、Ei、Fi，各设计变量取r个值，令r＝3， 在本实施例中i取值为1、2、3，如表1所示；对阶梯数据计算方案设计如下：各设计变量每 次取值的差值为ΔQA＝|A1-A2|＝|190-200|＝10、ΔQB＝|B1-B2|＝20、ΔQC＝|C1-
C2|＝100、ΔQD＝|D1-D2|＝10、ΔQE＝|E1-E2|＝100、ΔQF＝|F1-F2|＝100，给出各设计变 量相互影响作用表，并安排试验，在众多计算条件中按一定规则选出少数代表性强的计算方 案，确定各设计变量相互影响作用表，其中设计变量A与B间相互影响为A×B、设计变量A 与C间相互影响为A×C、设计变量B与C间相互影响为B×C，包括误差对比列，如图2所 示，表2内部的1、2、3分别为各设计变量在表1中的“取值”。传统枚举法需要做36＝729次仿真试验，按此表格科学安排试验，可将仿真次数减小至27次，即试验次数m＝27。给出设计变量间的相互影响作用表的表头，如表1：
表1设计变量间的相互影响作用表范例
| 列号 | A | B | A×B | A×B | C | A×C | A×C | B×C | D | 误差对比 | B×C | E | F | 数据yi |
| 1 | ||||||||||||||
| 2 | ||||||||||||||
| 3 | ||||||||||||||
| 4 | ||||||||||||||
| 5 | ||||||||||||||
| 6 | ||||||||||||||
| 7 | ||||||||||||||
| 8 | ||||||||||||||
| 9 | ||||||||||||||
| 10 |
该范例中有六个需要考察的设计变量，A、B、C，分别为内环半径、圆心到对称轴的距离、 圆心到地面的距离，D、E、F，为主均压环结构参数，分别为内环半径、圆心到对称轴的距 离、圆心到地面的距离，一列误差对比列，每个设计变量有3个待考察值(即图中白色格子 内标号i＝1、2、3)。传统枚举法需要做36＝729次仿真试验，按此表格科学安排试验，可将 仿真次数减小至27次。
S2、建立高压直流电阻标准分压器仿真模型；
S3、将S1中的数据，逐一代入到S2中建立的仿真模型中，得出高压直流电阻标准分压 器周围电场数据yi，汇总至S1阶梯数据计算方案中；
S4、根据S3中的数据计算各设计变量的指标，具体包括：
S41、计算总和Kij，即第j列上第i个水平的计算结果总和；
S42、计算总和平均值其中t为第j列上设计变量号i出现的次数；
S43、计算极差Rj，即第j列的极差；
S5、通过S4中极差的数值，计算各设计变量之间的效应，具体包括以下步骤：
S51、计算某两个设计变量之间的联合效应[ab]ij，[ab]ij＝μij-μ，其中μij为Ai与Bj的 组合下的均值，总均值为μ，该数值可反应Ai与Bj的组合对计算结果的总效应；
S52、计算某两个水平之间的交互效应(ab)ij，(ab)ij＝[ab]ij-ai-bj，该数值可反应Ai与Bj的组合对计算结果的交互效应；
S53、判断各设计变量影响是否显著，具体确定方法如下：
S531、计算结果的总偏差平方和ST，其中n＝每列同一设计变量的计 算次数与设计变量的个数相乘，T为数据总和，
S532、计算第j列设计变量的偏差平方和Sj，若设计变量A安排在 表中第j列上，则SA＝Sj；
S533、计算ST的度fT、Sj的度fj、Se的度fe，fT＝n-1，fj＝r-1，fe＝fT-各水平(包 括交互作用)的度之和，Se为所有误差对比列的偏差平方和之和；
S534、计算检验第j列设计变量对试验结果是否有显著影响的统计量的值Fj，若由观察值计算出的Fj≥F1-α(fj，fj)则认为该列设计变量对试验结果影响 显著，否则，认为效果不显著，其中α为给定的显著性水平，α可以取包括0.005、0.001在内 的数值。
给出方差分析表表头如表2所示：
表2方差分析表范例
S6、结合S4和S5的数据确定各设计变量A～F的取值，具体确定方法如下：
S61、以设计变量A和设计变量B为例，若设计变量A、B以及A×B对周围电场数据yi的影 响均大于系统误差对周围电场数据yi的影响，即设计变量A、B以及A×B的Sj均大于误差对 比列的Sj，则按照工程要求选择[ab]ij中最大或最小者；
S62、若设计变量B与设计变量A×B对周围电场数据yi的影响均大于系统误差对周围电 场数据yi的影响，设计变量A对周围电场数据yi的影响小于系统误差对周围电场数据yi的影 响，即设计变量A及A×B的影响显著，但设计变量B的影响不显著，则按照工程要求选择 (ab)ij及bj中最大或最小者；
S7、在S6得出的高压直流电阻标准分压器各均压环尺寸和位置的取值后，缩小各设计 变量之间差值为ΔQA′<ΔQA、ΔQB′<ΔQB、ΔQC′<ΔQC、ΔQD′<ΔQD、ΔQE′<ΔQE、 ΔQF′<ΔQF，迭代重复S1至S7，直到满足工程要求为止，具体迭代过程如下：
S71、选取新一轮迭代的各设计变量间差值ΔQA′、ΔQB′、ΔQC′、ΔQD′、ΔQE′、ΔQF′；
S72、确定新一轮迭代各设计变量的值A′k＝Ak±ΔQ′A、B′k＝Bk±ΔQ′B、Ck′＝Ck±ΔQ′C、 D′k＝Dk±ΔQ′D、E′k＝Ek±ΔQ′E、F′k＝Fk±ΔQ′F；
S73、把新一轮迭代值继续带入S1中的表中，制定新一轮方案。
下面以1100kV高压直流电阻标准分压器为研究对象。如图1所示，为本发明所提供的 基于数理统计方法和迭代法的高压直流电阻标准分压器均压环的设计方法，具体步骤包括：
S1、列出1100kV高压直流电阻标准分压器的主均压环和辅均压环需要考察的六个设计 变量A、B、C、D、E、F，分别为辅均压环的内环半径、辅均压环圆心到对称轴的距离、 辅均压环圆心到地面的距离、主均压环的内环半径、主均压环圆心到对称轴的距离、主均压环圆心到地面的距离，各设计变量的取值Ai、Bi、Ci、Di、Ei、Fi，各设计变量取r个值， 令r＝3，则在本实施例中i取值为1、2、3，如表1所示，各设计变量每次取值的差值为 ΔQA＝|A1-A2|＝|190-200|＝10、ΔQB＝|B1-B2|＝20、ΔQC＝|C1-C2|＝100、 ΔQD＝|D1-D2|＝10、ΔQE＝|E1-E2|＝100、ΔQF＝|F1-F2|＝100，给出各设计变量相互影响 作用表，并安排试验，在众多计算条件中按一定规则选出少数代表性强的计算方案，确定各 设计变量相互影响作用表，其中设计变量A与B间相互影响为A×B、设计变量A与C间相 互影响为A×C、设计变量B与C间相互影响为B×C，包括误差对比列，如图2所示，表2 内部的1、2、3分别为各设计变量在表1中的“取值”。传统枚举法需要做36＝729次仿真试 验，按此表格科学安排试验，可将仿真次数减小至m＝27次。
表1各设计变量的取值
表2各设计变量相互影响作用表
| 列号 | A | B | A×B | A×B | C | A×C | A×C | B×C | D | 误差对比 | B×C | E | F |
| 试验号 | |||||||||||||
| 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 |
| 2 | 1 | 1 | 1 | 1 | 2 | 2 | 2 | 2 | 2 | 2 | 2 | 2 | 2 |
| 3 | 1 | 1 | 1 | 1 | 3 | 3 | 3 | 3 | 3 | 3 | 3 | 3 | 3 |
| 4 | 1 | 2 | 2 | 2 | 1 | 1 | 1 | 2 | 2 | 2 | 3 | 3 | 3 |
| 5 | 1 | 2 | 2 | 2 | 2 | 2 | 2 | 3 | 3 | 3 | 1 | 1 | 1 |
| 6 | 1 | 2 | 2 | 2 | 3 | 3 | 3 | 1 | 1 | 1 | 2 | 2 | 2 |
| 7 | 1 | 3 | 3 | 3 | 1 | 1 | 1 | 3 | 3 | 3 | 2 | 2 | 2 |
| 8 | 1 | 3 | 3 | 3 | 2 | 2 | 2 | 1 | 1 | 1 | 3 | 3 | 3 |
| 9 | 1 | 3 | 3 | 3 | 3 | 3 | 3 | 2 | 2 | 2 | 1 | 1 | 1 |
| 10 | 2 | 1 | 2 | 3 | 1 | 2 | 3 | 1 | 2 | 3 | 1 | 2 | 3 |
| 11 | 2 | 1 | 2 | 3 | 2 | 3 | 1 | 2 | 3 | 1 | 2 | 3 | 1 |
| 12 | 2 | 1 | 2 | 3 | 3 | 1 | 2 | 3 | 1 | 2 | 3 | 1 | 2 |
| 13 | 2 | 2 | 3 | 1 | 1 | 2 | 3 | 2 | 3 | 1 | 3 | 1 | 2 |
| 14 | 2 | 2 | 3 | 1 | 2 | 3 | 1 | 3 | 1 | 2 | 1 | 2 | 3 |
| 15 | 2 | 2 | 3 | 1 | 3 | 1 | 2 | 1 | 2 | 3 | 2 | 3 | 1 |
| 16 | 2 | 3 | 1 | 2 | 1 | 2 | 3 | 3 | 1 | 2 | 2 | 3 | 1 |
| 17 | 2 | 3 | 1 | 2 | 2 | 3 | 1 | 1 | 2 | 3 | 3 | 1 | 2 |
| 18 | 2 | 3 | 1 | 2 | 3 | 1 | 2 | 2 | 3 | 1 | 1 | 2 | 3 |
| 19 | 3 | 1 | 3 | 2 | 1 | 3 | 2 | 1 | 3 | 2 | 1 | 3 | 2 |
| 20 | 3 | 1 | 3 | 2 | 2 | 1 | 3 | 2 | 1 | 3 | 2 | 1 | 3 |
| 21 | 3 | 1 | 3 | 2 | 3 | 2 | 1 | 3 | 2 | 1 | 3 | 2 | 1 |
| 22 | 3 | 2 | 1 | 3 | 1 | 3 | 2 | 2 | 1 | 3 | 3 | 2 | 1 |
| 23 | 3 | 2 | 1 | 3 | 2 | 1 | 3 | 3 | 2 | 1 | 1 | 3 | 2 |
| 24 | 3 | 2 | 1 | 3 | 3 | 2 | 1 | 1 | 3 | 2 | 2 | 1 | 3 |
| 25 | 3 | 3 | 2 | 1 | 1 | 3 | 2 | 3 | 2 | 1 | 2 | 1 | 3 |
| 26 | 3 | 3 | 2 | 1 | 2 | 1 | 3 | 1 | 3 | 2 | 3 | 2 | 1 |
| 27 | 3 | 3 | 2 | 1 | 3 | 2 | 1 | 2 | 1 | 3 | 1 | 3 | 2 |
S2、建立高压直流电阻标准分压器仿真模型，如图2所示；
S3、将S1中的数据，如A1＝190，A2＝200，A3＝210逐一代入到S2中建立的仿真 模型中，得出高压直流电阻标准分压器周围电场数据yi，汇总至S1阶梯数据计算方案中如 表3所示：
表3电场数据yi
| 列号 | A | B | A×B | A×B | C | A×C | A×C | B×C | D | 误差对比 | B×C | E | F | 数据yi |
| 试验号 | ||||||||||||||
| 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | -1781.74 |
| 2 | 1 | 1 | 1 | 1 | 2 | 2 | 2 | 2 | 2 | 2 | 2 | 2 | 2 | -1686.23 |
| 3 | 1 | 1 | 1 | 1 | 3 | 3 | 3 | 3 | 3 | 3 | 3 | 3 | 3 | -1602.63 |
| 4 | 1 | 2 | 2 | 2 | 1 | 1 | 1 | 2 | 2 | 2 | 3 | 3 | 3 | -1592.75 |
| 5 | 1 | 2 | 2 | 2 | 2 | 2 | 2 | 3 | 3 | 3 | 1 | 1 | 1 | -1821.18 |
| 6 | 1 | 2 | 2 | 2 | 3 | 3 | 3 | 1 | 1 | 1 | 2 | 2 | 2 | -1769.53 |
| 7 | 1 | 3 | 3 | 3 | 1 | 1 | 1 | 3 | 3 | 3 | 2 | 2 | 2 | -1599.48 |
| 8 | 1 | 3 | 3 | 3 | 2 | 2 | 2 | 1 | 1 | 1 | 3 | 3 | 3 | -1599.45 |
| 9 | 1 | 3 | 3 | 3 | 3 | 3 | 3 | 2 | 2 | 2 | 1 | 1 | 1 | -1885.04 |
| 10 | 2 | 1 | 2 | 3 | 1 | 2 | 3 | 1 | 2 | 3 | 1 | 2 | 3 | -1607.42 |
| 11 | 2 | 1 | 2 | 3 | 2 | 3 | 1 | 2 | 3 | 1 | 2 | 3 | 1 | -1653.05 |
| 12 | 2 | 1 | 2 | 3 | 3 | 1 | 2 | 3 | 1 | 2 | 3 | 1 | 2 | -1818.13 |
| 13 | 2 | 2 | 3 | 1 | 1 | 2 | 3 | 2 | 3 | 1 | 3 | 1 | 2 | -1665.26 |
| 14 | 2 | 2 | 3 | 1 | 2 | 3 | 1 | 3 | 1 | 2 | 1 | 2 | 3 | -1612.26 |
| 15 | 2 | 2 | 3 | 1 | 3 | 1 | 2 | 1 | 2 | 3 | 2 | 3 | 1 | -1724.47 |
| 16 | 2 | 3 | 1 | 2 | 1 | 2 | 3 | 3 | 1 | 2 | 2 | 3 | 1 | -161058 |
| 17 | 2 | 3 | 1 | 2 | 2 | 3 | 1 | 1 | 2 | 3 | 3 | 1 | 2 | -1751.3 |
| 18 | 2 | 3 | 1 | 2 | 3 | 1 | 2 | 2 | 3 | 1 | 1 | 2 | 3 | -1657.97 |
| 19 | 3 | 1 | 3 | 2 | 1 | 3 | 2 | 1 | 3 | 2 | 1 | 3 | 2 | -1568.54 |
| 20 | 3 | 1 | 3 | 2 | 2 | 1 | 3 | 2 | 1 | 3 | 2 | 1 | 3 | -1672.67 |
| 21 | 3 | 1 | 3 | 2 | 3 | 2 | 1 | 3 | 2 | 1 | 3 | 2 | 1 | -1762.23 |
| 22 | 3 | 2 | 1 | 3 | 1 | 3 | 2 | 2 | 1 | 3 | 3 | 2 | 1 | -1659.56 |
| 23 | 3 | 2 | 1 | 3 | 2 | 1 | 3 | 3 | 2 | 1 | 1 | 3 | 2 | -157736 |
| 24 | 3 | 2 | 1 | 3 | 3 | 2 | 1 | 1 | 3 | 2 | 2 | 1 | 3 | -1715.01 |
| 25 | 3 | 3 | 2 | 1 | 1 | 3 | 2 | 3 | 2 | 1 | 2 | 1 | 3 | -1617.63 |
| 26 | 3 | 3 | 2 | 1 | 2 | 1 | 3 | 1 | 3 | 2 | 3 | 2 | 1 | -1696.39 |
| 27 | 3 | 3 | 2 | 1 | 3 | 2 | 1 | 2 | 1 | 3 | 1 | 3 | 2 | -1654.53 |
S4、根据S3中的数据计算各设计变量的指标，具体包括：
S41、计算总和Kij，即第j列上第i个水平的计算结果总和；
S42、计算总和平均值其中t为第j列上设计变量号i出现的次数；
S43、极差Rj，即第j列的极差；
S5、通过S4中极差的数值，计算各设计变量之间的效应，具体包括以下步骤：
S51、计算某两个设计变量之间的联合效应[ab]ij，[ab]ij＝μij-μ，其中μij为Ai与Bj的 组合下的均值，总均值为μ，该数值可反应Ai与Bj的组合对计算结果的总效应；
S52、计算某两个水平之间的交互效应(ab)ij，(ab)ij＝[ab]ij-ai-bj，该数值可反应Ai与Bj的组合对计算结果的交互效应；
S53、判断各设计变量影响是否显著，具体确定方法如下：
S531、计算结果的总偏差平方和ST，其中n＝每列同一设计变量的计 算次数与设计变量的个数相乘，T为数据总和，
S532、计算第j列设计变量的偏差平方和Sj，若设计变量A安排在 表中第j列上，则SA＝Sj；
S533、计算ST的度fT、Sj的度fj、Se的度fe，fT＝n-1，fj＝r-1，fe＝fT-各水平(包 括交互作用)的度之和，Se为所有误差对比列的偏差平方和之和。如表4所示。
S534、计算检验第j列设计变量对试验结果是否有显著影响的统计量的值Fj，若由观察值计算出的Fj≥F1-α(fj，fj)，则认为该列设计变量对试验结果影响显著，否则，认为效果不显著，其中α为给定的显著性水平，可取包括0.005、0.001在内 的数值，如表5所示。
表4各计算指标结果
表5方差分析表
S6、结合S4和S5的数据确定各设计变量的最终取值A*、B*、C*、D*、E*、F*分别为210、640、8700、390、1200、7900，此时yi为1568.54，具体确定方法如下：
S61、若某两个设计变量A、B以及A×B的影响都显著，则按照工程要求选择[ab]ij中最大或最小者；
S62、若设计变量A及A×B的影响显著，但设计变量B的影响不显著，则按照工程要求选择(ab)ij及bj中最大或最小者，如表6所示；
表6汇总表
| [ab]ij表 | (ab)i表 | |||||||
| B1 | B2 | B3 | B1 | B2 | B3 | |||
| A1 | -10.1115 | -47.7315 | -14.5681 | A1 | 17.5637 | -21.7519 | 4.188148 | |
| A2 | -12.7781 | 12.75852 | 6.805185 | A2 | -11.5019 | 12.33926 | -0.83741 | |
| A3 | 12.27519 | 29.44519 | 23.90519 | A3 | -6.06185 | 9.412593 | -3.35074 | |
| [bc]ij表 | (bc)ij表 | |||||||
| C1 | C2 | C3 | C1 | C2 | C3 | |||
| B1 | 22.09852 | -22.1981 | -72.3115 | B1 | -20.7896 | -24.3163 | -16.6907 | |
| B2 | 52.33519 | 7.885185 | -53.4348 | B2 | 7.751481 | 4.071481 | 0.49037 | |
| B3 | 64.84519 | 31.28185 | -30.5015 | B3 | 13.03815 | 20.24481 | 16.20037 |
S7、在S6得出的高压直流电阻标准分压器各均压环尺寸和位置的取值后，缩小各设计 变量之间差值为ΔQA′、ΔQB′、ΔQC′、ΔQD′、ΔQE′、ΔQF′，迭代重复S1至S7，直到满足工程 要求为止，具体迭代过程如下：
S71、选取新一轮迭代的各设计变量间差值ΔQ′A、ΔQ′B、ΔQ′C、ΔQ′D、ΔQ′E、ΔQ′F，要求 ΔQA′＜ΔQA、ΔQB′＜ΔQB、ΔQC′＜ΔQC、ΔQD′＜ΔQD、ΔQE′＜ΔQE、ΔQF′＜ΔQF，分别取 ΔQA＝5、ΔQB＝10、ΔQC＝50、ΔQD＝5、ΔQE＝50、ΔQF＝50。
S72、确定新一轮迭代各设计变量的值Ak′、Bk′、Ck′、Dk′、Ek′、Fk′，A′k＝Ak±ΔQ′A、B′k＝Bk±ΔQ′B、C′k＝Ck±ΔQ′C、D′k＝Dk±ΔQ′D、E′k＝Ek±ΔQ′E、F′k＝Fk±ΔQ′F，Ak、 Bk、Ck、Dk、Ek、Fk为上一轮设计变量值；
S73、把新一轮迭代值继续带入S1中的表中，制定新一轮方案，最终使分压器周围最大 场强小于等于分压器的起晕场强。
本实施例可证明，在工程实际中此方式不仅能缩短设计高压直流电阻标准分压器均压环 所需时间，且多次迭代能使得结果更趋近于最优值，具有实际应用价值。
最后应说明的是：以上所述的各实施例仅用于说明本发明的技术方案，而非对其限制； 尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然 可以对前述实施例所记载的技术方案进行修改，或者对其中部分或全部技术特征进行等同替 换；而这些修改或替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (6)
1.一种确定高压直流电阻标准分压器均压环的优化设计方法，其特征在于，包括以下步骤：
S1、确定主均压环和辅均压环的设计变量，对所述设计变量多次取值，确定包括误差对比列在内的设计变量相互影响作用表，确定阶梯数据计算方案；
S2、建立高压直流电阻标准分压器仿真模型；
S3、将S1中各设计变量的取值数据逐一代入到S2中建立的高压直流电阻标准分压器仿真模型中，得出高压直流电阻标准分压器的周围电场数据yi，汇总至阶梯数据计算方案中；
S4、根据S3中的数据计算各设计变量的指标，具体包括：
S41、计算第j列上第i个设计变量取值的计算结果总和Kij；
计算第j列上各设计变量在i值下的叠加总和Kij
S42、计算总和Kij的平均值 其中t为第j列上设计变量号i出现的次数；
S43、计算第j列极差Rj，
S5、通过S4中极差的数值，计算各设计变量之间的效应，具体包括以下步骤：
S51、计算各设计变量之间的联合效应[ab]ij；
S52、计算各设计变量之间的交互效应(ab)ij；
S53、判断各设计变量对周围电场数据yi的影响；
S6、结合S4和S5的数据确定各设计变量的取值；以及
S7、在S6得出的高压直流电阻标准分压器各均压环尺寸和位置的取值后，缩小各设计变量之间差值，迭代重复S1至S7，使分压器周围最大场强小于等于分压器的起晕场强。
2.根据权利要求1所述的确定高压直流电阻标准分压器均压环的优化设计方法，其特征在于，步骤S1中设计变量具体包括：辅均压环的内环半径A、辅均压环圆心到对称轴的距离B、辅均压环圆心到地面的距离C、主均压环的内环半径D、主均压环圆心到对称轴的距离E、主均压环圆心到地面的距离F，各设计变量的取值Ai、Bi、Ci、Di、Ei、Fi，各设计变量取r个值，i＝1～r，各设计变量每次取值的差值为ΔQA＝|A1-A2|、ΔQB＝|B1-B2|、ΔQC＝|C1-C2|、ΔQD＝|D1-D2|、ΔQE＝|E1-E2|、ΔQF＝|F1-F2|。
3.根据权利要求1所述的确定高压直流电阻标准分压器均压环的优化设计方法，其特征在于，步骤S51和S52中计算的具体步骤为：
计算两个设计变量之间的联合效应[ab]ij，[ab]ij＝μij-μ，其中μij为Ai与Bj的组合下的均值，总均值为μ，该数值可反应Ai与Bj的组合对计算结果的总效应；
计算两个设计变量不同i值之间的交互效应(ab)ij，(ab)ij＝[ab]ij-ai-bj，该数值可反应Ai与Bj的组合对计算结果的交互效应。
4.根据权利要求1所述的确定高压直流电阻标准分压器均压环的优化设计方法，其特征在于，步骤S53中的具体步骤为：
计算结果的总偏差平方和ST，其中n＝每列同一设计变量的计算次数与设计变量的个数相乘，T为数据总和，
计算第j列设计变量的偏差平方和Sj，若设计变量A安排在表中第j列上，则SA＝Sj；
计算误差对比列的偏差平方和Se，Se为所有误差对比列的偏差平方和之和；
计算总偏差平方和ST的度fT、Sj的度fj、Se的度fe，fT＝n-1，fj＝r-1，fe为fT与各设计变量的度之和；
计算检验第j列设计变量对试验结果是否有显著影响的统计量值Fj，若由观察值计算出的Fj≥F1-α(fj，fj)，则认为该列设计变量对试验结果影响显著，否则，认为效果不显著，其中α为给定的显著性水平。
5.根据权利要求1或2所述的确定高压直流电阻标准分压器均压环的优化设计方法，其特征在于，步骤S6中的具体步骤为：
取设计变量A和设计变量B，计算设计变量A、B间的相互作用A×B，若设计变量A、B以及A×B对周围电场数据yi的影响均大于系统误差对周围电场数据yi的影响，则选择[ab]ij中最大或最小者；
若设计变量B与设计变量A×B对周围电场数据yi的影响均大于系统误差对周围电场数据yi的影响，设计变量A对周围电场数据yi的影响小于系统误差对周围电场数据yi的影响，则选择(ab)ij及bj中最大或最小者。
6.根据权利要求1或2所述的确定高压直流电阻标准分压器均压环的优化设计方法，其特征在于，步骤S7中的具体步骤为：
选取新一轮迭代的各设计变量间差值为ΔQA′、ΔQB′、ΔQC′、ΔQD′、ΔQE′、ΔQF′，确定新一轮迭代各设计变量的值为ΔQ′A、ΔQ′B、ΔQ′C、ΔQ′D、ΔQ′E、ΔQ′F，要求ΔQA′＜ΔQA、ΔQB′＜ΔQB、ΔQC′＜ΔQC、ΔQD′＜ΔQD、ΔQE′＜ΔQE、ΔQF′＜ΔQF，Ak＝Ak±ΔQA、Bk＝Bk±ΔQB、C′k＝Ck±ΔQ′C、D′k＝Dk±ΔQ′D、E′k＝Ek±ΔQ′E、F′k＝Fk±ΔQ′F，k为迭代次数；
把新一轮迭代值继续带入S1中的设计变量相互影响作用表中，形成修正后设计变量相互影响作用表。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811195074.XA CN109408929B (zh) | 2018-10-15 | 2018-10-15 | 确定高压直流电阻标准分压器均压环的优化设计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811195074.XA CN109408929B (zh) | 2018-10-15 | 2018-10-15 | 确定高压直流电阻标准分压器均压环的优化设计方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109408929A true CN109408929A (zh) | 2019-03-01 |
| CN109408929B CN109408929B (zh) | 2022-12-13 |
Family
ID=65467896
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811195074.XA Active CN109408929B (zh) | 2018-10-15 | 2018-10-15 | 确定高压直流电阻标准分压器均压环的优化设计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109408929B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110470264A (zh) * | 2019-09-11 | 2019-11-19 | 国网冀北电力有限公司电力科学研究院 | 直流电阻分压器的均压环尺寸确定方法及装置 |
| CN110514159A (zh) * | 2019-09-11 | 2019-11-29 | 国网冀北电力有限公司电力科学研究院 | 直流电阻分压器的均压环尺寸及位置确定方法及装置 |
| CN110795870A (zh) * | 2019-09-25 | 2020-02-14 | 中国电力科学研究院有限公司 | 一种直流高压发生器均压环的优化设计方法 |
| CN112528534A (zh) * | 2020-11-19 | 2021-03-19 | 南方电网科学研究院有限责任公司 | 直流分压器的表面最大电场强度的获取方法、系统及装置 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103683269A (zh) * | 2013-12-24 | 2014-03-26 | 国家电网公司 | 一种考虑外绝缘污秽的直流分压器均压环优化设计方法 |
| CN104200096A (zh) * | 2014-08-29 | 2014-12-10 | 中国南方电网有限责任公司超高压输电公司昆明局 | 基于微分进化算法和bp神经网络的避雷器均压环优化法 |
-
2018
- 2018-10-15 CN CN201811195074.XA patent/CN109408929B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103683269A (zh) * | 2013-12-24 | 2014-03-26 | 国家电网公司 | 一种考虑外绝缘污秽的直流分压器均压环优化设计方法 |
| CN104200096A (zh) * | 2014-08-29 | 2014-12-10 | 中国南方电网有限责任公司超高压输电公司昆明局 | 基于微分进化算法和bp神经网络的避雷器均压环优化法 |
Non-Patent Citations (5)
| Title |
|---|
| YANZHEN ZHAO等: "Calculation of electric field and voltage distribution along glass insulator string and optimum design of its grading rings for 1000kV a.c. UHV transmission lines", 《2008 WORLD AUTOMATION CONGRESS》 * |
| 张宗喜等: "特高压±1600kV直流发生器空间电场分布仿真分析", 《中国电力》 * |
| 樊阳文等: "±1100 kV直流电子式电压互感器绝缘特性研究", 《高压电器》 * |
| 章述汉等: "±800kV换流站直流电压互感器现场校准试验", 《高电压技术》 * |
| 袁渊等: "±1200kV标准直流分压器设计及测量不确定度评定", 《高压电器》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110470264A (zh) * | 2019-09-11 | 2019-11-19 | 国网冀北电力有限公司电力科学研究院 | 直流电阻分压器的均压环尺寸确定方法及装置 |
| CN110514159A (zh) * | 2019-09-11 | 2019-11-29 | 国网冀北电力有限公司电力科学研究院 | 直流电阻分压器的均压环尺寸及位置确定方法及装置 |
| CN110795870A (zh) * | 2019-09-25 | 2020-02-14 | 中国电力科学研究院有限公司 | 一种直流高压发生器均压环的优化设计方法 |
| CN112528534A (zh) * | 2020-11-19 | 2021-03-19 | 南方电网科学研究院有限责任公司 | 直流分压器的表面最大电场强度的获取方法、系统及装置 |
| CN112528534B (zh) * | 2020-11-19 | 2024-05-28 | 南方电网科学研究院有限责任公司 | 直流分压器的表面最大电场强度的获取方法、系统及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109408929B (zh) | 2022-12-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109408929A (zh) | 确定高压直流电阻标准分压器均压环的优化设计方法 | |
| CN108920751B (zh) | 基于拓扑优化的电力变压器绕组变形状态逆向求解方法 | |
| CN104679959B (zh) | 1000kV特高压电容分压器宽频等效电路建模方法 | |
| CN111896906A (zh) | 对互感器的误差校准数据进行修正的方法及装置 | |
| CN111262246A (zh) | 一种变电站接地网快速模型建立与腐蚀准确诊断方法 | |
| CN108761241A (zh) | 一种射频同轴连接器的非线性效应建模方法 | |
| CN105258847A (zh) | 一种压力传感器校准的方法和装置 | |
| CN106021189A (zh) | 一种适应多种边界条件的多回电缆稳态温升获取方法 | |
| CN103632235A (zh) | 基于改进前推回代算法的电网低压台区线损信息处理方法 | |
| US20240012968A1 (en) | Performance evaluation method of lng ambient air vaporizer | |
| CN110110406B (zh) | 一种基于Excel计算平台实现LS-SVM模型的边坡稳定性预测方法 | |
| CN104965134B (zh) | 基于回声状态网络的变电站接地网地表电位分布计算方法 | |
| CN106018970A (zh) | 一种双层土壤中接地网接地电阻的测量方法 | |
| CN108593964A (zh) | 一种七孔气流探针的标定方法 | |
| CN113283784B (zh) | 基于区间优化算法的可靠性逆问题模型及其求解方法 | |
| CN204925189U (zh) | 一种10kV直流电压比例标准装置 | |
| CN117235950B (zh) | 基于牛顿迭代法的天然气管网稳态仿真方法、介质及设备 | |
| CN105426656B (zh) | 一种误差计算方法及系统 | |
| CN108877969A (zh) | 核功率理论模型建立及验证方法、系统及终端设备 | |
| CN107505536B (zh) | 计及多根金属管道的变电站接地网地表电位分布计算方法 | |
| CN111159936A (zh) | 一种基于广义时域有限差分的电缆接头热场计算方法 | |
| CN107122829B (zh) | 一种利用虚拟样本训练神经网络诊断变压器故障的方法 | |
| CN108763830A (zh) | 半导体器件的闪烁噪声模型及其提取方法 | |
| KR101169936B1 (ko) | 비보정인자와 위상오차를 측정하는 전압변성기 측정시스템, 이를 이용한 평가방법, 및 그 기록매체 | |
| CN103149472B (zh) | 一种基于二端子激励测量模式的电压-电流映射构造方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |