1.一种开关类设备操作强电磁骚扰特性分析方法，其特征在于，包括以下步骤：步骤1：暂态强电磁骚扰的测量，包括利用倒置式电容分压器测量暂态过电压、利用罗氏线圈和无源RC积分器测量暂态大电流、利用电场测量系统测量脉冲电场、利用磁场测量系统测量脉冲磁场；

步骤2：研究AIS和GIS变电站开关类设备操作下传导、耦合和辐射暂态电磁骚扰的产生机理、传播途径和分布规律；

步骤3：通过大量实测，研究AIS及GIS变电站高压场域二次设备安装关键位置电磁骚扰时域、频域分布规律，给出变电站关键位置骚扰的具体特征；

开展AIS和GIS开关类设备操作强电磁暂态过程仿真计算：建立不同电站瞬态电磁环境三维空间仿真模型，静电场和感应电场建模分开建模，以模拟装置实际现场布置为依据，利用COMSOL搭建三维瞬态电磁场有限元模型，仿真中，考虑设备的最小尺寸，采用自由四面体作为剖分单元，多次细化剖分后发现结果未改变，表明剖分网格已满足计算需求，根据具体导线截面尺寸和测量点与导线的距离，决定是否在仿真中忽略导线截面尺寸与集肤效应带来的影响，基于该模型研究空间瞬态电磁场与电压电流波形的幅值频率间的关系，以及空间电磁场分布特性；

步骤4：研究开关类设备操作强电磁骚扰对典型二次设备致损的作用机理，以及在此工况下二次设备敏感元器件的易损性评估及失效敏感度分析，考虑实际情况下对不同二次设备致损机理所掌握的程度不同，将不同的二次设备分类并发展相应的易损性评估方法；

所述步骤4中将不同的二次设备分三类进行易损性评估，第一类设备，其强电磁骚扰作用下的物理机理与统计特征均已知，易损性模型可以通过完整的数学表达式表示，进而采用不同的统计方法进行模型参数回归，最终以比例优势模型进行评估；第二类设备，只有部分的致损机理已知但无法表示出完整的效应评估过程，评估中采用由描述致损机理的白箱模型与假设的黑箱模型组成的混合模型用于描述效应易损性模型，采用贝叶斯混合模型评估；第三类设备，考虑复杂二次设备，其致损机理与阈值概率特征均未知，无法用数学模型加以描述，评估中采用数据驱动的非参模型来描述，采用基于支持向量机和高斯过程回归的非参模型评估；所述第三类设备的评估过程是：将多种环境参数作为评价指标评估强电磁干扰对二次电气设备的易损性，包括脉冲辐射电场幅值、端口耦合电流幅值、端口耦合电流能量积分，用于评估分析的试验数据可表示为：n

其中，xi∈R，是n维评价指标组成向量，i{0,1,…,N}；yi{0,1,…,m}，0代表正常状态，m为失效状态的等级数，则总的输出维数为m+1维，m≥1且m∈N；

分类结果的发生概率为在给定训练数据时求：

s.t.P0(xi)+P1(xi)+…+Pm(xi)＝1分别针对各个等级进行评估计算给出各等级的中间分类指标，进而比较各等级的分类指标得到其发生概率，在具体计算中间分类指标时，先采用SVM方法将训练数据从离散量转化为连续的分类指数；基于此，采用GPR方法对评估空间进行分类指标的回归计算，这样便可以综合SVM和GPR各自在的优势，实现模型总体上的无参化；

1)采用SVM方法将训练数据从离散量转化为连续的分类指数将所研究的某等级数据视为正类，用“1”表示，其余所有数据视为负类，用“‑1”表示，那么针对某等级数据，训练数据集可以被变换为：l

其中，l是所研究的等级，l∈{0,1,…,m}；yi∈{+1,‑1}；

如果训练数据线性可分，其分离超平面可以表示为：

其中，ωl为法向量，bl为截距；

而相应的分类决策函数可以表示为：

如果 则样本点分类结果为正类；如果 则样本点分类

结果为负类；如果 则样本点在分界面上；分类超平面(l,bl)关于样本点l(xi,yji)的几何间隔可以表示分类预测的准确程度，几何间隔由γ i表示，l其中样本点为xi，其类标记为yi＝+1；

进一步定义分类超平面(l,bl)关于训练数据集Tl的几何间隔为超平面(l,bl)关于Tl中所有样本点的几何间隔最小值：l

相对应的，分类超平面(ωl,bl)关于样本点(xi,yi)和训练数据集Tl的函数间隔可以分别表示为：而它们与几何间隔之间的关系可以表示为：

用约束最优化问题来表示：

l

即使得分类超平面ωl,bl)关于训练数据集Tl的几何间隔γ 最大化，同时约束条件表示l分类超平面(ωl,bl)关于每个样本点的几何间隔至少是γ ；

先基于几何间隔与函数间隔的关系，即式(11)，式(12)将改写为：为了保证解是唯一的，需要对 做一些限制，取 即将全局的函数间隔定义为1，也即是将离超平面最近的点的距离定义为1/||ωl||，求1/||ωl||的最大值相当于求||ωl||

2/2的最小值，因此改写后结果为：

根据式求得分类超平面的法向量ωl和截距bl；

2)采用GPR方法对评估空间进行分类指标的回归计算

a)基于高斯过程的先验模型建立

在GPR建模过程中，首先取模型先验为高斯过程模型，通过对均值函数和协方差函数的确定来表示，其中均值函数表示对目标函数的代理模型，协方差函数表示模型的不确定度，其中，m(x)＝E[f(x)]为均值函数； 为协方差函数，先验均值函数表示为一组基函数线性组合：

T

m(x)＝s(x)β       (16)

T

其中，s(x)＝[s1(x),s2(x),…,sp(x)]为基函数组成的向量；β为p×1的回归参数向量；

实际上，为了运算方便，先验均值函数取0；

协方差函数确定着样本间的相互依赖关系，预测数据与已知数据的距离越大，其协方差就越大，结果就越不确定，协方差函数的选取对回归结果有较大的影响，是高斯过程模型中重要的组成部分，采用核函数表示，通用的表达式为：2

其中，σ为方差参数； 为相关核，θ为参数向量；

b)基于联合高斯分布的后验模型回归

采用SVM对分类指数的计算，可以得到用于建立后验模型的训练数据，即：考虑到在测量环节中引入的噪声，实际所得的训练数据应表示为模型实际值与噪声之和，即：2

其中，ε为均值为0，方差为σn的高斯噪声；

那么基于高斯过程先验的定义，设先验均值函数为0，先验协方差函数为 则训练数据X＝[x1,x2,…,xN]和预测数据x的概率分布分别表示为：fl(X)～N(0,Kl+Δl)                     (20)fl(x)～N(0,kl(x,x))                     (21)其中，

基于(20)和(21)，二者的联合高斯分布可以表示为：

T

其中，kl(x)＝kl(x,x)，kl(x)＝[kl(x,xi)]i＝1,…,N；

由联合高斯分布的边缘分布性质可得预测数据x服从分布：

其中：

c)基于“一对其他”策略的多等级效应概率评估方法

将SVM和GPR结合可以得到对预测空间任一点处分类指数分布的估计，即(23)～(25)，得到针对类别l进行的两分类回归预测，采用“一对其他”的策略，重复以上的步骤，便可以得到m+1个分类指数分布模型，即：假设分类指数越大的效应状态，那么结果将最终表现为该状态，通过对各个效应状态分类指数的积分，就可以得到在x点处状态i的发生概率：同时，各状态发生概率的和为1，即：

2.根据权利要求1所述的一种开关类设备操作强电磁骚扰特性分析方法，其特征在于，上述步骤1中为了提高暂态过电压测量准确度，在进行暂态测量过程，应将稳态基准在同一时间段同时测量，即要求高速采集变换器有较大的数据存储量，同时记录稳态和暂态波形数据。

3.根据权利要求1所述的一种开关类设备操作强电磁骚扰特性分析方法，其特征在于，上述步骤1中所述暂态大电流的测量具体是：通过罗氏线圈测量电流，得到与电流成比例的微分电压信号，再利用无源RC积分器对微分电压信号进行积分，所述罗氏线圈的线圈骨架选用聚四氟乙烯材料，罗氏线圈设有三层铝合金制作的屏蔽外壳，对杂散磁场和电场进行屏蔽，提高罗氏线圈的测量结果准确性，所述无源RC积分器包括穿心电容器和无引线金属膜电阻、同轴全屏蔽结构，电阻选用高频特性较好的金属膜电阻，去掉电阻引线，导电直接用电阻端部的帽盖金属固件，以减少引线带来的电感，电容选用穿心电容器；穿心电容器结构只有上下两个平板电极，本身残余电感极小，环形电容器无引脚线，可以将引线电感减到最低，导电采用上下环形导电面，导电面分别与同轴外壳及导电片连接，这种接线方式相比普通引脚线电容器具有接线电感小、导电面积大、接触电阻小的优点，电阻和电容器安装在一个全封闭的同轴金属壳体内，输入和输出采用同轴连接器，同轴结构可降低导体的传导电感，金属全封闭结构可增强抗电磁干扰性能。

4.根据权利要求1所述的一种开关类设备操作强电磁骚扰特性分析方法，其特征在于，上述步骤1中的所述电场测量系统由电场接收天线、光发射机与光接收机组成，所述磁场测量系统由B‑dot天线、光发射机与光接收机组成。

5.根据权利要求4所述的一种开关类设备操作强电磁骚扰特性分析方法，其特征在于，上述步骤1中测量脉冲磁场的过程是：由宽频带的B‑dot天线接收空间中的脉冲磁场信号，该信号经光发射机进行电光转换后利用光纤进行远距离传输，在接收端利用光接收机进行光电转换并最终将接收信号输入示波器。

6.根据权利要求1所述的一种开关类设备操作强电磁骚扰特性分析方法，其特征在于，上述步骤3中建立不同电站瞬态电磁环境三维空间仿真模型后，要将试验中有限的测量点的电磁场与模型中对应点的瞬态电磁场波形进行对比以验证模型的准确性，在验证模型准确性后，才可基于该模型研究，验证模型的准确性的方法是：将实测瞬态电流时域信号转化到频域，即信号分解为不同频率正弦波叠加，对于某一个特定点，将每一个频率的瞬态电流分量进行电磁场计算，获得该点不同频率下瞬态电磁场的幅值和相位，再通过快速傅里叶逆变换将结果转化为时域，即可获得该点的电磁场波形。