基于多物理场耦合的变压器绕组振动声纹分析方法及系统与流程

1.本技术涉及变压器绕组振动噪声技术领域，尤其涉及一种基于多物理场耦合的变压器绕组振动声纹分析方法及系统。


背景技术：

2.作为电力系统的核心设备之一，电力变压器的安全性与可靠性对电网的可靠运行至关重要。据统计数据表明，变压器的机械故障主要由长期累积形成，包括变压器的主件和部分构件，其中由于绕组变形所引起的故障占绝大比例。
3.长期以来，人们一直将变压器发出的声音视为噪声，但却忽略其内在价值，其运行声音包含着变压器运行状态信息。变压器的噪声主要是由铁芯磁致伸缩以及绕组受洛伦兹力振动产生。在变压器短路状态下，由于电压比较小，铁芯磁致伸缩的影响很小，所产生的声音与绕组振动产生的声音存在数量级上的差异，因此在这种运行状态下，只需要考虑变压器绕组振动所产生的声音。
4.目前国内外在变压器的声音信号的产生机理、降噪方法、声信号特征提取以及识别方法中已经取得了大量的成果，但在变压器正常及故障状态下的声音信号特性分析方面仍未有深入研究。这使得有效的特征量提取方法以及故障识别算法缺乏依据。


技术实现要素：

5.有鉴于此，本技术提供了一种基于多物理场耦合的变压器绕组振动声纹特性分析方法，包括如下步骤：
6.步骤一：根据实际尺寸建立变压器几何模型，并对变压器几何模型进行网格划分；
7.步骤二：根据变压器实际运行情况，仿真变压器电磁场分布情况；
8.步骤三：根据变压器绕组中磁场与结构力场的耦合，仿真变压器绕组受力以及振动情况；
9.步骤四：根据变压器绕组结构力场与压力声场耦合，仿真变压器外围空间声场；
10.步骤五：依据变压器绕组常见三种机械故障，仿真变压器绕组在故障状态下的声纹特性；
11.步骤六：进行外围声场仿真结果分析以及实验验证分析。
12.进一步地，所述步骤一包括：
13.根据实际参数对变压器进行几何建模，对绕组上下两端施加固定约束，模拟构件加紧作用，铁芯采用硅钢片叠压而成，绕组采用均匀多匝形式，并对变压器有限元模型进行有效的网格划分。
14.进一步地，所述步骤二中，
15.变压器的磁场包括主磁场和漏磁场，漏磁场是变压器绕组振动的根源。
16.进一步地，变压器磁感应强度计算方法包括：
17.当绕组之中流过电流时，由于交变漏磁场的存在，使得绕组上产生电磁力，假设稳
态运行下流过绕组的电流i为：
18.i＝imsin(ωt φ0)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
19.其中im为绕组电流幅值，ω为绕组中电流的角频率，φ0为绕组电流的初相位；
20.简化磁感应强度计算公式，在静态条件下，磁感应强度可表示为：
[0021][0022]
其中为磁感应强度与电流之间的电动力系数。
[0023]
进一步地，所述步骤三包括：
[0024]
计算变压器绕组所受的电动力为：
[0025][0026]
变压器绕组等效为多自由度线性弹簧质量系统，其固体力学微分方程可以表示为：
[0027][0028]
式中：m为模型质量矩阵；c为模型阻尼矩阵；k为模型弹性系数矩阵；z(t)、式中：m为模型质量矩阵；c为模型阻尼矩阵；k为模型弹性系数矩阵；z(t)、分别为绕组模型的形变位移、形变速度和形变加速度；f(t)为绕组所受的电动力大小；g为重力加速度。
[0029]
变压器绕组模型的振动形变位移计算如下：
[0030]
将式(3)所推导出的电动力公式代入(4),可得：
[0031][0032]
由式(5)可知，m、c、k皆为常数，其为二阶常系数微分方程，对于其齐次部分：
[0033][0034]
其通解可以表示为：
[0035]
z＝e
αt
(c
1 cosβt c
2 sinβt)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(7)
[0036]
式中：c1、c2为常数，其大小由初始条件确定；
[0037][0038][0039]
由式(5)可知，m、c、k、r、φ0等皆为常数，其为二阶常系数微分方程，求解方程可得：
[0040][0041]
其中：a、θ为常数且有初始条件决定；
[0042][0043][0044][0045][0046]
所述步骤四中，声压的波动方程表示为
[0047][0048]
其中：
[0049]
p(r,t)为声压，其为传播半径r和时间t的函数；c0为声速.
[0050]
进一步地，所述步骤四中，
[0051]
利用comsol中的声-结构边界模块，耦合固体力学与压力声学模块，在绕组表面与空气介质接触面向外辐射噪声，将固体力学仿真数据作为声场分析的基础，分别计算变压器前后以及左右两个面的辐射噪声。
[0052]
本发明还公开了一种基于前述方法的变压器绕组振动声纹特性分析系统，包括变压器建模模块、变压器电磁场分布仿真模块、变压器绕组受力及振动仿真模块、变压器外围空间声场仿真模块、变压器绕组声纹特性仿真模块；
[0053]
所述变压器建模模块根据实际尺寸建立变压器几何模型，并对变压器几何模型进行网格划分；
[0054]
所述变压器电磁场分布仿真模块根据变压器实际运行情况，仿真变压器电磁场分布；
[0055]
所述变压器绕组受力及振动仿真模块根据变压器绕组中磁场与结构力场的耦合，仿真变压器绕组受力以及振动情况；
[0056]
所述变压器外围空间声场仿真模块根据变压器绕组结构力场与压力声场耦合，仿真变压器外围空间声场；
[0057]
所述变压器绕组声纹特性仿真模块依据变压器绕组常见三种机械故障，仿真变压器绕组在故障状态下的声纹特性。
[0058]
进一步优选，
[0059]
所述变压器绕组振动声纹特性分析系统还包括外围声场仿真结果分析以及实验验证分析模块，分析变压器外围声场状态，对于变压器噪声进行分析，对比仿真结果，验证模型有效性。
[0060]
本发明同时公开了一种计算机系统，包括处理器及存储介质；
[0061]
所述存储介质用于存储指令；
[0062]
所述处理器用于根据所述指令进行操作以执行根据基于多物理场耦合的变压器绕组振动声纹特性分析方法的步骤。
[0063]
一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现基于多物理场耦合的变压器绕组振动声纹特性分析方法的步骤。
[0064]
本发明研究了三相油浸式变压器的振动噪声特性，对s13-m-200/10型号的油浸式变压器进行短路试验，并测量了表面的振动信号和噪声信号。根据真实的试验条件，建立基于有限元的变压器绕组振动噪声模型，并考虑变压器油对声学传播过程的影响。变压器的振动及声压仿真结果与实际测量结果误差很小，验证了计算模型及方法的有效性。在模型的基础上建立变压器机械故障的仿真模型，分析变压器绕组故障状态下的噪声特性，为变压器绕组故障诊断提供了依据。
附图说明
[0065]
图1是本发明的方法流程图。
[0066]
图2是变压器网格划分模型正视图。
[0067]
图3是变压器网格划分模型俯视图。
[0068]
图4是变压器外接仿真电路图图。
[0069]
图5是变压器绕组感应电流图。
[0070]
图6是变压器绕组磁通密度模图。
[0071]
图7是变压器铁芯磁通密度模图。
[0072]
图8是变压器绕组辐向受力图。
[0073]
图9是高压绕组受力形变图。
[0074]
图10是低压绕组受力形变图。
[0075]
图11是变压器绕组轴向受力图。
[0076]
图12是10kv变压器声信号测点图。
[0077]
图13是绕组松动故障模拟图。
[0078]
图14是绕组变形故障模拟图。
[0079]
图15是绕组绝缘脱落故障模拟图。
[0080]
图16是变压器正面声场图。
[0081]
图17是变压器侧面声场图。
[0082]
图18是变压器正面实验图。
[0083]
图19是变压器侧面实验图。
具体实施方式
[0084]
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例的附图，对本发明实施例的技术方案进行清楚、完整地描述。显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。基于所描述的本发明的实施例，本领域普通技术人员所获得的所有其他实施例，都属于本发明保护的范围。
[0085]
如图1所示，本发明提供了一种基于多物理场耦合的变压器绕组振动声纹特性分析方法，包括如下步骤：
[0086]
步骤一：以s13-m-200/10型号的油浸式变压器仿真原型，建立变压器几何模型，并对变压器模型进行网格划分。
[0087]
变压器结构复杂，由铁芯、绕组、油箱及其他辅助构件组成，本发明根据实际参数对变压器进行几何建模。本仿真主要构件包括铁芯、绕组、绝缘油、油箱壁、上夹件、下夹件、木制垫块、铁芯夹件和上顶板等构成。对绕组上下两端施加固定约束，模拟构件加紧作用，铁芯采用硅钢片叠压而成，绕组采用均匀多匝形式。
[0088]
对变压器有限元模型进行合理有效的网格划分，对于结构精细的部分，单元划分也更加精细，从而减少计算量，提高计算精度。
[0089]
图2示出了变压器网格划分模型正视图，图3示出了变压器网格划分模型俯视图。
[0090]
步骤二：根据变压器实际运行情况，仿真变压器绕组电磁场分布情况。
[0091]
变压器的磁场主要包括主磁场和漏磁场，漏磁场是变压器绕组振动的根源。当绕组之中流过电流时，由于交变漏磁场的存在，使得绕组上产生电磁力。假设稳态运行下流过绕组的电流为：
[0092]
i＝i
m sin(ωt φ0)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0093]
其中im为绕组电流幅值，ω为绕组中电流的角频率，φ0为绕组电流的初相位。
[0094]
在漏磁通流经非铁磁材料时，主要经过主气隙、绕组、压紧结构或是油箱闭合，并且主要是轴向磁密分量。简化磁感应强度计算公式，在静态条件下，磁感应强度bv
t
可表示为：
[0095][0096]
其中为磁感应强度与电流之间的电动力系数。
[0097]
变压器电路模型采用dyn11接线方式，其高压侧电压源的电压为400v,频率为50hz，且相位相差120
°
，外接仿真电路如图4所示，变压器绕组感应电流如图5所示，高压侧感应出的电流有效值为11.5a,低压侧感应出的电流有效值为288.7a,仿真结果与实际相符。
[0098]
图6是变压器绕组磁通密度模图，图7是变压器铁芯磁通密度模图，如图6-7所示，漏磁通主要集中于高压绕组与低压绕组的气隙之中，最大值为0.0028t，变压器的主磁通主要集中铁芯之中，其最大值为1.16t，可知在激励电流的作用下，磁通密度较小，没有达到饱和磁通，且主磁通约占全部磁通的99％以上，漏磁通仅占全部磁通的1％以下，与实际情况相符。
[0099]
步骤三：根据变压器绕组中磁场与结构力场的耦合，仿真变压器绕组受力以及振动情况。
[0100]
计算变压器绕组所受的电动力为：
[0101][0102]
式中，指磁感应强度，r指线圈半径，i指稳态运行下流过绕组的电流，im为绕组电流幅值，为磁感应强度与电流之间的电动力系数。
[0103]
根据上式可知，绕组线圈的振动角频率为2ω，对于基波频率为50hz的电网来说，100hz是变压器绕组振动的主要频率。由于电网中存在其他谐波电流分量以及在变压器中非线性材料的作用下，变压器绕组振动信号中也存在着少量的其他频率的分量。
[0104]
变压器绕组是一个典型的复杂的多自由度结构，可以等效为多自由度线性弹簧质量系统，其固体力学微分方程可以表示为：
[0105][0106]
式中：m为模型质量矩阵；c为模型阻尼矩阵；k为模型弹性系数矩阵；z(t)、量矩阵；c为模型阻尼矩阵；k为模型弹性系数矩阵；z(t)、分别为绕组模型的形变位移、形变速度和形变加速度；f(t)为绕组所受的电动力大小；g为重力加速度。
[0107]
将式(3)所推导出的电动力公式代入(4)，可得：
[0108][0109]
由式(5)可知，m、c、k皆为常数，其为二阶常系数微分方程，对于其齐次部分：
[0110][0111]
其通解可以表示为：
[0112]
z＝e
αt
(c
1 cosβt c
2 sinβt)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(7)
[0113]
式中：c1、c2为常数，其大小由初始条件确定；
[0114][0115][0116]
由式(5)可知，m、c、k、r、φ0等皆为常数，其为二阶常系数微分方程，求解方程可得：
[0117][0118]
其中：a、θ为常数且有初始条件决定；
[0119][0120][0121][0122][0123]
由上式可知，绕组的振动位移大小取决于质量m、阻尼c、弹性系数k及绕组的几何结构，而对于确定的变压器，其质量以及阻尼式固定不变的，因此其振动位移大小主要取决
于绕组的弹性系数及绕组的几何结构。
[0124]
变压器绕组在电磁力的作用下，受到辐向以及轴向两个方向的力。辐向力使得变压器绕组收到向内或者向外的牵引力作用，而轴向的作用力，则使得绕组受到向上或者向下的牵引力作用。
[0125]
变压器绕组辐向受力如图8所示，b相和c相绕组延y轴方向受力也相应较大，因此在沿y轴方向上的形变程度也更为明显，而沿x轴方向上的受力相对较小。可知绕组所受的辐向力，在同一时刻的不同方向是不相同的，其大小由漏磁场的强度决定。如图9和10所示，高压绕组收到向外的牵引力，绕组对应有向外膨胀的趋势，而低压绕组受到向内的牵引力，绕组对应有向内收缩的趋势。并且越靠近绕组中心位置形变的程度越大。
[0126]
如图11所示，b相和c相绕组上表面受到较大的向上的牵引力，而其他位置的轴向力相对来说较小，但是由于绕组上下表面施加了相对应的约束，因此没有发生大的形变，可知在生产以及使用过程中，对于绕组进行有效的固定，有助于防止绕组发生严重的形变。
[0127]
步骤四：变压器绕组压力声学仿真建模：根据变压器绕组结构力场与压力声场耦合，仿真变压器外围空间声场。
[0128]
声波在介质中传播，但是在不同介质之中，声波的传播特性以及波的振动形式都是不同的，根据物理学三个基本定律：质量守恒定律、能量守恒定律和动量守恒定律，推导出流体之中的三个基本方程：连续性定律、状态方程和运动方程。根据声波振动的条件，并通过对三个方程的消元，可获得声压的波动方程：
[0129][0130]
其中：p(r,t)为声压，其为传播半径r和时间t的函数；c0为声速。
[0131]
变压器绕组是变压器主要的振动源之一，其他辅助构件的振动噪声很小，本发明在研究绕组特性的仿真中忽略变压器其他辅助构件的振动噪声。利用comsol中的声-结构边界模块，耦合固体力学与压力声学模块，在绕组表面与空气介质接触面向外辐射噪声，将固体力学仿真数据作为声场分析的基础。分别计算变压器前后以及左右两个面的辐射噪声，将仿真测点布置离变压器30cm厘米处，高度为变压器油箱高度的一半。振动及声音信号仿真测点布置如图12所示。
[0132]
步骤五：仿真变压器绕组三种机械故障，仿真变压器绕组在故障状态下的声纹特性。
[0133]
变压器绕组的松动程度即变压器绕组预紧力大小，与垫块及压紧结构件的材料特性相关。仿真几何模型在保证变压器绕组圈数不变的情况下，对变压器绕组松动、绕组变形及绕组绝缘脱落等故障进行模拟仿真。
[0134]
垫块以及压紧结构件的材料非线性，导致其弹性模量随预紧力大小变化而变化，本发明通过改变变压器b相绕组弹性模量大小模拟b相绕组的松动故障，仿真模型如图13所示。
[0135]
变压器绕组的变形通常是由绕组轴向力产生，尤其是变压器绕组中流过短路电流或者发生励磁涌流时，此时绕组中流过的大电流对绕组损耗尤为严重。由上述仿真可知，变
压器高压绕组受到向外的牵引力，使得绕组向外膨胀，而变压器低压绕组向内收缩，使得低压绕组向内扭曲变形。本发明通过将圆形绕组变成具有明显凸起的六边形圆柱体，来模拟变压器绕组变形，仿真模型如图14所示。
[0136]
变压器绕组绝缘脱落会对变压器安全稳定运行造成巨大影响，本发明通过改变变压器b相绕组的匝数来模拟变压器绕组绝缘脱落故障仿真模型如图15所示。
[0137]
步骤六：变压器外围声场仿真实验验证
[0138]
研究变压器周围的声场情况，本仿真简化为绕组的接触面在电磁场中运动而产生。
[0139]
搭建实验室变压器数据测量平台，并分别对如调压器、采集仪等持续存在的噪声以及汽车声，脚步声等随机出现的噪声进行数据采集。完成对变压器短路运行情况下的数据采集，并采集不同位置的变压器噪声信号，变压器正面以及侧面的声场如图16-17所示，实验测点布置如18-19所示，分析变压器外围声场状态，对于变压器噪声进行分析，对比仿真结果，验证模型有效性。
[0140]
基于本发明的精神，本领域技术人员能够容易想到基于前述变压器绕组振动声纹特性分析方法可以得到一种计算机程序产品。计算机程序产品可以包括计算机可读存储介质，其上载有用于使处理器实现本公开的各个方面的计算机可读程序指令。即本技术还包括一种计算机系统，包括处理器及存储介质；
[0141]
所述存储介质用于存储指令；
[0142]
所述处理器用于根据所述指令进行操作以执行根据前述变压器绕组振动声纹特性分析方法的步骤。
[0143]
一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现前述变压器绕组振动声纹特性分析方法的步骤。
[0144]
计算机可读存储介质可以是可以保持和存储由指令执行设备使用的指令的有形设备。计算机可读存储介质例如可以是――但不限于――电存储设备、磁存储设备、光存储设备、电磁存储设备、半导体存储设备或者上述的任意合适的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括：便携式计算机盘、硬盘、随机存取存储器(ram)、只读存储器(rom)、可擦式可编程只读存储器(eprom或闪存)、静态随机存取存储器(sram)、便携式压缩盘只读存储器(cd-rom)、数字多功能盘(dvd)、记忆棒、软盘、机械编码设备、例如其上存储有指令的打孔卡或凹槽内凸起结构、以及上述的任意合适的组合。这里所使用的计算机可读存储介质不被解释为瞬时信号本身，诸如无线电波或者其它自由传播的电磁波、通过波导或其它传输媒介传播的电磁波(例如，通过光纤电缆的光脉冲)、或者通过电线传输的电信号。
[0145]
这里所描述的计算机可读程序指令可以从计算机可读存储介质下载到各个计算/处理设备，或者通过网络、例如因特网、局域网、广域网和/或无线网下载到外部计算机或外部存储设备。网络可以包括铜传输电缆、光纤传输、无线传输、路由器、防火墙、交换机、网关计算机和/或边缘服务器。每个计算/处理设备中的网络适配卡或者网络接口从网络接收计算机可读程序指令，并转发该计算机可读程序指令，以供存储在各个计算/处理设备中的计算机可读存储介质中。
[0146]
用于执行本公开操作的计算机程序指令可以是汇编指令、指令集架构(isa)指令、
机器指令、机器相关指令、微代码、固件指令、状态设置数据、或者以一种或多种编程语言的任意组合编写的源代码或目标代码，所述编程语言包括面向对象的编程语言—诸如smalltalk、c 等，以及常规的过程式编程语言—诸如“c”语言或类似的编程语言。计算机可读程序指令可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络—包括局域网(lan)或广域网(wan)—连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。在一些实施例中，通过利用计算机可读程序指令的状态信息来个性化定制电子电路，例如可编程逻辑电路、现场可编程门阵列(fpga)或可编程逻辑阵列(pla)，该电子电路可以执行计算机可读程序指令，从而实现本公开的各个方面。以上实例仅用于说明本发明的设计思想和特点，其目的在于使本领域内的技术人员能够了解本发明的内容并据以实施，本发明的保护范围不限于上述实例，所以，凡依据本发明所揭示的原理、设计思路所作的等同变化或修饰，例如改变评价指标等，均在本发明的保护范围之内。
[0147]
申请人结合说明书附图对本发明的实施例做了详细的说明与描述，但是本领域技术人员应该理解，以上实施例仅为本发明的优选实施方案，详尽的说明只是为了帮助读者更好地理解本发明精神，而并非对本发明保护范围的限制，相反，任何基于本发明的发明精神所作的任何改进或修饰都应当落在本发明的保护范围之内。