本发明涉及一种变压器套管内部电容芯子绝缘结构设计，尤其涉及一种基于kriging模型的套管电容芯子电场优化方法。
背景技术：
电力变压器是构成电网的主要设备，在电力系统中起着电压等级变换的作用。随着电力系统输电电压的提高，电力变压器的电压等级也不断提高，这就对变压器套管的绝缘特性提出了更高的要求。套管是变压器箱外的主要绝缘装置，它将变压器的高压和低压绕组的引出线引出到油箱外部，同时起固定引线的作用。目前较为常用的变压器套管有油浸式和电容式两种，其中电容式套管多用于100kv以上的电力变压器。
电容式套管的主要绝缘结构为电容芯子，它是将绝缘纸和铝箔加压力交替卷在导电杆上形成的。这样的设计相当于以导电杆为中心串联多个电容器，根据串联电路分压原理，导电杆对接地法兰的电压应等于各个电容器的电压之和，通过控制各层的电容量可使电场均匀的分布在套管内部。套管中的最大场强值直接关系到变压器的正常运行，而且50％以上的变压器事故都是由套管故障引起的，因此对套管内部电场进行最优化设计尤为重要。
套管电容芯子中电场分布的均匀程度直接影响到套管内部的最大场强值，而最大场强值则主要受电容芯子的间距和厚度影响。传统的电场计算方法很难得到电容芯子的最优参数，且得到的参数无法进行验证缺少可信度。
因此，寻找一种能够对电容芯子最优参数进行可靠预测的方法具有重要的实用价值。
技术实现要素：
发明的主要目的在于克服现有技术的缺点与不足，提供一种基于kriging模型的套管电容芯子电场优化方法。
该方法利用仿真软件计算出电容芯子和变压器油的最大场强值，通过kriging模型预测最优解，最后对最优解进行验证，预测结果可靠且计算量小。
为实现上述目的，本发明提供如下技术方案：
一种基于kriging模型的套管电容芯子电场优化方法，包括以下步骤：
(1)利用cad参数化建模功能,建立电容式变压器套管模型；
(2)选择电容芯子间距和厚度为标注约束参数，确定研究范围，得到电容芯子的结构参数；
(3)利用cad参数化建模功能依次改变电容芯子的结构参数，并将模型导入仿真软件comsol中，最终计算出各个参数下变压器油和电容芯子的最大场强值；
(4)利用步骤(3)得到的结果初步分析变压器油和电容芯子的最大场强值与电容芯子的间距和厚度之间的关系；
(5)构建kriging模型对电容芯子的结构参数进行优化处理；
(6)得到最优化的电容芯子结构参数，并将优化参数导入comsol进行验证。
进一步地，所述步骤(1)中cad参数化建模功能具体包括：将模型中各个部分设置为相互约束关系，包括相切、相等、垂直、同心约束关系，并选择电容芯子之间的距离和厚度为标注约束参数，通过改变标注约束参数值进而实现对整个模型的调整。
进一步地，所述步骤(2)中电容芯子的间距在1.0～2.35mm范围内等间隔选取10个数值进行研究，电容芯子厚度在0.2～0.35mm范围内等间隔选取10个数值进行研究。
进一步地，所述步骤(3)具体为：
(31)将改变标注约束参数值后得到的变压器套管模型导入comsol软件中；
(32)对模型进行条件加载，包括：各部分材料属性的添加、边界条件的设定及接地和电势的接入；
(33)运行软件计算出电容芯子内部的电场分布，最终得到电容芯子和油间隙中最大场强的数值及出现的位置。
进一步地，所述步骤(5)中构建kriging模型的具体步骤为：
(51)设n×p试验样本点x＝{x1，x2…xn},
y＝{y1(x),y2(x)…yn(x)}；
(52)据输出函数，未知参数β，σ2可估计为：
β，＝(ftr-1f)-1ftr-1y
式中：f为每个试验点f(x)估值向量；r为试验点相关的函数矩阵；
(53)由于β，σ2均与相关参数θ相关，因此需要先求得θ值，用最大似然估计法对表达式求最小值：
(54)获取θ后，通过响应最佳线性预测可得β，σ2值。
进一步地，所述步骤(6)中所述最优化的电容芯子结构参数即为使得套管内部变压器油和电容芯子的最大场强值最低时，电容芯子间距和厚度的参数值。
进一步地，所述步骤(6)中验证方法具体为：将所得最优化参数代入模型中，再将模型导入仿真软件comsol中计算出变压器油和电容芯子的最大场强值，并将所得结果与优化之前获得的结果进行比较。
本发明与现有技术相比，具有如下优点和有益效果：
1、本发明采用参数化建模技术，只需要设定各部分之间的约束关系，选择电容芯子之间的距离和厚度为标注约束参数，通过改变标注约束参数即可实现对模型的修改，无需重新建模；
2、本发明采用kriging模型对电容芯子内部电场进行优化处理，计算量小且结果可靠。
附图说明
图1是本发明实施例的流程示意图。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。
实施例1
如图1所示，基于kriging模型的套管电容芯子电场优化方法包括以下步骤：
s1、利用cad参数化建模功能,建立电容式变压器套管模型。
本步骤采用cad参数化建模技术，将模型的各个部分设置成各种约束关系,并对套管内部电容芯子结构进行细化，具体如：将模型中各个部分设置为相互约束关系，包括相切、相等、垂直、同心约束关系，并选择电容芯子之间的距离和厚度为标注约束参数，通过改变标注约束参数值进而实现对整个模型的调整。
影响电容芯子中最大场强值的因素有很多，如中心导电杆的粗细、均压球的大小、铝箔的层数、变压器油中水分和气泡含量、温度和承受电压类型等等。在实际运行中变压器油温度、水分和气泡等难以人为控制，且导电杆的粗细、均压球的大小都有变压器的相关电气特性作为设计依据，而电容芯子的间距和厚度是影响最大场强值的可控因素。因此步骤s1中选择电容芯子之间的距离和厚度为标注约束参数，通过改变标注约束参数值进而实现对整个模型的调整。
s2、选择电容芯子间距和厚度为标注约束参数，确定研究范围，得到电容芯子的结构参数。
本步骤重点研究电容芯子厚度和间距对套管内电场分布的影响，因此选择间距和厚度为标注约束参数。按照设计要求，选取电容芯子间距和厚度的研究范围分别为：1～2.35mm和0.2～0.35mm，等间距选取10组结构参数进行研究。
s3、利用cad参数化建模功能依次改变电容芯子的结构参数，并将模型导入仿真软件comsol中，最终计算出各个参数下变压器油和电容芯子的最大场强值，具体包括：
s31、将改变标注约束参数值后得到的变压器套管模型导入comsol软件中；
s32、对模型进行条件加载，包括：各部分材料属性的添加、边界条件的设定及接地和电势的接入；
s33、运行软件计算出电容芯子内部的电场分布，最终得到电容芯子和油间隙中最大场强的数值及出现的位置。
本步骤用有限元法建立变压器套管的电场计算模型，给模型各部分赋予相应的材料属性，并对电容芯子结构中的铝箔施加悬浮电位。在对模型进行剖分时，对剖分精度不做太高要求，可以减少计算量，对结果并无太大影响。在将中心导电杆接入电势，套管外侧屏蔽罩设置接地后，通过计算得到结果，绘制电场等势面图后分析电场分布。记录每一组参数下得到的电容芯子内部和变压器油中最大场强值，这样便得到一组完整的样本数据。
s4、利用步骤s3、得到的结果初步分析变压器油和电容芯子的最大场强值与电容芯子的间距和厚度之间的关系。
s5、建立kriging模型，对电容芯子的结构参数进行优化处理。
初步建立kriging模型，并利用样本数据进行训练，确定kriging模型结构。kriging方法为数学、地质学广泛使用，是基于随机过程的随机预测法，其根据区域内若干信息样品的某种特征数据对该区域同类特征未知数做线性无偏、最小方差估计，具有平滑效应及估计方差最小的统计特征，构建kriging模型具体包括步骤：
s51、设n×p试验样本点x＝{x1，x2…xn},
y＝{y1(x),y2(x)…yn(x)}；
s52、据输出函数，未知参数β，σ2可估计为：
β’＝(ftr-1f)-1ftr-1y
式中：f为每个试验点f(x)估值向量；r为试验点相关的函数矩阵；
s53、由于β，σ2均与相关参数θ相关，因此需要先求得θ值，用最大似然估计法对表达式求最小值：
s54、获取θ后，通过响应最佳线性预测可得β，σ2值。
s6、得到最优化的电容芯子结构参数，并将优化参数导入comsol进行验证。
将上述过程得到的最优解模型重新导入comsol中，按照上述过程进行仿真计算，将得到结果与预测结果进行比较。
本实施例采用有限元方法进行仿真计算，并利用kriging方法进行最优化求解，若得到广泛推广，必将降低电容式套管在运行时的故障率。
上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。