CN112685928B - 一种基于三相电抗器声源模型的噪声预测方法及系统 - Google Patents

一种基于三相电抗器声源模型的噪声预测方法及系统 Download PDF

Info

Publication number
CN112685928B
CN112685928B CN202011517182.1A CN202011517182A CN112685928B CN 112685928 B CN112685928 B CN 112685928B CN 202011517182 A CN202011517182 A CN 202011517182A CN 112685928 B CN112685928 B CN 112685928B
Authority
CN
China
Prior art keywords
sound
sound source
noise
model
line
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.)
Active
Application number
CN202011517182.1A
Other languages
English (en)
Other versions
CN112685928A (zh
Inventor
倪园
郭兆枫
周兵
王延召
胡静竹
张建功
陈传敏
干喆渊
张业茂
谢辉春
刘震寰
赵军
路遥
李妮
刘兴发
刘健犇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
North China Electric Power University
Electric Power Research Institute of State Grid Hunan Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
North China Electric Power University
Electric Power Research Institute of State Grid Hunan Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, China Electric Power Research Institute Co Ltd CEPRI, North China Electric Power University, Electric Power Research Institute of State Grid Hunan Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202011517182.1A priority Critical patent/CN112685928B/zh
Publication of CN112685928A publication Critical patent/CN112685928A/zh
Application granted granted Critical
Publication of CN112685928B publication Critical patent/CN112685928B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明公开了一种基于三相电抗器声源模型的噪声预测方法及系统,属于电力系统设备噪声控制与预测技术领域。本发明方法,包括:对测点处的噪声声强或声压进行测量,获取每个测点处的声强值或声压值;根据声强值或声压值,确定线声源的声功率;搭建三相电抗器声源模型的井型立体线等效模型,搭建防火墙有限元模型,搭建边界元模型;对声压数据和声压级三维立体数据绘制二维多切面云图,根据二维多切面云图,预测线声源的近、远场噪声的分布情况。本发明对声源模型的精准化建立及声传播预测提供了新的理论依据与方法指导,利用本发明,可以用建立三相电抗器精准的声源模型,从而为声传播预测与新建变电站方案设计提供准确的仿真数据。

Description

一种基于三相电抗器声源模型的噪声预测方法及系统
技术领域
本发明涉及电力系统设备噪声控制与预测技术领域,并且更具体地,涉及一种基于三相电抗器声源模型的噪声预测方法及系统。
背景技术
随着社会的不断发展,各个行业的用电量持续增加,但是我国能源分布为西多东少,北多南多,70%左右的电力消费集中在东部以及中部地区,因此为满足远距离大容量输电的需求,我国已建成了多个1000kV特高压输变电工程。在工程投运后的近几年里,虽已有不少企业借助Cadna/A、SoundPlan等常用噪声预测软件对特高压变电站开展了变电站噪声的评估与治理研究,但在多次现场实测后发现,1000kV特高压变电站并联电抗器所产生的噪声并不是按照几何发散规律随距离的增加而单调的衰减,而是在附近空间上出现若干极大、极小值点,这说明特高压变电站并联电抗器噪声的声场存在明显干涉现象。
导致常用噪声预测软件所得的计算结果与实测结果存在一定差距,尤其是三相电抗器等设备附近区域计算值与实测值差别相对较大,其主要原因噪声预测计算所建立的声源模型与实际运行设备模型存在一定差距,且预测的源强参数不准确,传播过程中的干涉效应无法得到有效计算,导致预测结果与实测结果误差较大。
建立准确的三相电抗器预测模型,开展特高压变电站噪声预测技术研究,准确且有效的预测变电站噪声,是提高变电站规划设计合理性的必然要求,也是保障工程运行后厂界环境噪声排放达标的重要手段。
目前,已有研究将三相电抗器声学模型简化为以下几种:
1、点声源模型:将三相电抗器等效为一个点声源,基于点声源户外传播理论对高抗产生的噪声进行计算。此方法简单方便且容易计算,但是使用条件十分有限,仅在预测高抗远场噪声时较为准确,即当预测点距离高抗足够远(一般大于高抗尺寸的3倍时)可将其等效为一个点声源。当预测点与高抗距离较近时,噪声计算误差较大,特别是在特高压变电站中,三相电抗器体积较大,且具有三个声源,将其等效为一个点声源不太合理。
2、面声源模型:现阶段在电力设计院或科研高校,在计算三相电抗器噪声时,均是将高抗等效为一个面声源,或为高抗的横截面,或为高抗的立截面,此模型更接近于三相电抗器的实际情况,计算简单且方便,但在计算近场区时,噪声预测仍然不准确。
3、体声源模型:将三相电抗器简化为一个长方体,并忽略高抗顶和箱底振动的辐射噪声,认为噪声全部来源于整体的振动,计算基于Helmholtz积分公式,通过高抗整体表面法向振动加速度的相关计算来分析三相电抗器的声场声压,这种方法比简化为点声源模型时更接近于三相电抗器的实际情况,但计算得到的变压器声场声压仍然不是十分精确。
4、扬声器阵列:此模型是在将三相电抗器看作面声源的基础上,用一系列的扬声器来代替各个面,通过设置不同的扬声器幅值和初始相角的组合以得到任何需要得到的声场,实现声场重建。此方法十分复杂,难以实现扬声器的设置与三相电抗器原始声场一致。
5、有限元模型:建立三相电抗器铁芯和绕组的有限元模型,对其进行多物理场耦合分析、电磁分析、结构分析以及声场分析得到三相电抗器辐射声场的分布情况。根据有限元理论基础,计算分析了三相电抗器声场分布情况,利用试验中测得三相电抗器表面法向振动加速度分量数据进行了计算。将其作为求解声场分布的边界条件,声学计算基于Helmholtz积分公式,是一种理论计算与试验数据相结合的方法。但是由于计算量过大,计算结果不容易收敛,且随着预测范围的扩大,计算量将成倍增加,导致无法计算远场噪声。
发明内容
针对上述问题,本发明提出了一种基于三相电抗器声源模型的噪声预测方法,包括:
对三相电抗器进行井型立体拆分,获取线声源,并确定线声源的长短和位置;针对线声源,根据线声源的长短及位置,以预设的间隔设置测点,对测点处的噪声声强或声压进行测量,获取每个测点处的声强值或声压值;
根据声强值或声压值,确定线声源的声功率;
搭建三相电抗器声源模型的井型立体线等效模型,并将符合预设条件的声功率输入等效模型;搭建防火墙有限元模型,并针对防火墙有限元模型设置声固耦合边界条件;搭建边界元模型,并确定模型的噪声预测点及预测范围;
针对输入声功率的等效模型、防火墙有限元模型及边界元模型,使用统一坐标系,根据声学波动方程及波动方程对噪声预测点及预测范围进行噪声计算,获取边界元范围内的声压数据和声压级三维立体数据,对声压数据和声压级三维立体数据绘制二维多切面云图,根据二维多切面云图,预测线声源的近、远场噪声的分布情况。
可选的,对三相电抗器进行井型立体拆分,包括:
根据电抗器结构,在电抗器垂直方向三等分,在长水平方向三等分,在宽水平方向二等分,按等分线进行切分,将每相电抗器拆分成18个长方体,并获得45条线声源。
可选的,根据声强值或声压值,确定线声源的声功率,具体为:
对获取的每个测点处的声强值或声压值,根据等效线声源的个数进行平均分组或不平均分组;
对每个等效线声源组中测点的噪声声强级合成值和测量表面面积的附加声强级,进行叠加计算,获取各个等效线声源的声功率。
可选的,预设条件为声功率相差120°的相位。
可选的,噪声预测点具有不同方位,在不同方位内的预测点,分别确定相应方位内的预测点,使用统一坐标系,以三相电抗器A相或C相的任一底角为坐标原点,确定每个线声源及预测点的坐标。
本发明还提出了一种基于三相电抗器声源模型的噪声预测系统,包括:
测量模块,对三相电抗器进行井型立体拆分,获取线声源,并确定线声源的长短和位置;针对线声源,根据线声源的长短及位置,以预设的间隔设置测点,对测点处的噪声声强或声压进行测量,获取每个测点处的声强值或声压值;
计算模块,根据声强值或声压值,确定线声源的声功率;
模型搭建模块,搭建三相电抗器声源模型的井型立体线等效模型,并将符合预设条件的声功率输入等效模型;搭建防火墙有限元模型,并针对防火墙有限元模型设置声固耦合边界条件;搭建边界元模型,并确定模型的噪声预测点及预测范围;
预测模块,针对输入声功率的等效模型、防火墙有限元模型及边界元模型,使用统一坐标系,根据声学波动方程及波动方程对噪声预测点及预测范围进行噪声计算,获取边界元范围内的声压数据和声压级三维立体数据,对声压数据和声压级三维立体数据绘制二维多切面云图,根据二维多切面云图,预测线声源的近、远场噪声的分布情况。
可选的,对三相电抗器进行井型立体拆分,包括:
根据电抗器结构,在电抗器垂直方向三等分,在长水平方向三等分,在宽水平方向二等分,按等分线进行切分,将每相电抗器拆分成18个长方体,并获得45条线声源。
可选的,根据声强值或声压值,确定线声源的声功率,具体为:
对获取的每个测点处的声强值或声压值,根据等效线声源的个数进行平均分组或不平均分组;
对每个等效线声源组中测点的噪声声强级合成值和测量表面面积的附加声强级,进行叠加计算,获取各个等效线声源的声功率。
可选的,预设条件为声功率相差120°的相位。
可选的,噪声预测点具有不同方位,在不同方位内的预测点,分别确定相应方位内的预测点,使用统一坐标系,以三相电抗器A相或C相的任一底角为坐标原点,确定每个线声源及预测点的坐标。
本发明弥补了常用商业软件无法计算干涉效应的缺点,同时考虑到三相电抗器的三个相各不相同且存在防火墙,将高抗的声功率设置三种不同的相位,对不同相位的高抗之间设立声固耦合边界,同时考虑到声场的三维性,对每个相的高抗都采用立体井型线声源进行建立,使仿真噪声具有立体性,具有更高的真实还原度及预测精度。
本发明对声源模型的精准化建立及声传播预测提供了新的理论依据与方法指导,利用本发明,可以用建立三相电抗器精准的声源模型,从而为声传播预测与新建变电站方案设计提供准确的仿真数据。
附图说明
图1为本发明方法流程图;
图2为本发明方法实施例流程图;
图3为本发明立体井型线声源在电抗器上的分布示意图;
图4为本发明基于FEM声固耦合及FEM-BEM耦合方法的三相电抗器立体井型线声源模型噪声仿真计算出的声压级云图及切片图;
图5为本发明基于基于FEM声固耦合及FEM-BEM耦合方法的三相电抗器立体井型线声源模型噪声仿真计算出的声压云图及切片图
图6为本发明系统结构图。
具体实施方式
现在参考附图介绍本发明的示例性实施方式,然而,本发明可以用许多不同的形式来实施,并且不局限于此处描述的实施例,提供这些实施例是为了详尽地且完全地公开本发明,并且向所属技术领域的技术人员充分传达本发明的范围。对于表示在附图中的示例性实施方式中的术语并不是对本发明的限定。在附图中,相同的单元/元件使用相同的附图标记。
除非另有说明,此处使用的术语(包括科技术语)对所属技术领域的技术人员具有通常的理解含义。另外,可以理解的是,以通常使用的词典限定的术语,应当被理解为与其相关领域的语境具有一致的含义,而不应该被理解为理想化的或过于正式的意义。
本发明提出了一种基于三相电抗器声源模型的噪声预测方法,如图1所示,包括:
对三相电抗器进行井型立体拆分,获取线声源,并确定线声源的长短和位置;针对线声源,根据线声源的长短及位置,以预设的间隔设置测点,对测点处的噪声声强或声压进行测量,获取每个测点处的声强值或声压值;
根据声强值或声压值,确定线声源的声功率;
搭建三相电抗器声源模型的井型立体线等效模型,并将符合预设条件的声功率输入等效模型;搭建防火墙有限元模型,并针对防火墙有限元模型设置声固耦合边界条件;搭建边界元模型,并确定模型的噪声预测点及预测范围;
针对输入声功率的等效模型、防火墙有限元模型及边界元模型,使用统一坐标系,根据声学波动方程及波动方程对噪声预测点及预测范围进行噪声计算,获取边界元范围内的声压数据和声压级三维立体数据,对声压数据和声压级三维立体数据绘制二维多切面云图,根据二维多切面云图,预测线声源的近、远场噪声的分布情况。
其中,对三相电抗器进行井型立体拆分,包括:
根据电抗器结构,在电抗器垂直方向三等分,在长水平方向三等分,在宽水平方向二等分,按等分线进行切分,将每相电抗器拆分成18个长方体,并获得45条线声源。
其中,根据声强值或声压值,确定线声源的声功率,具体为:
对获取的每个测点处的声强值或声压值,根据等效线声源的个数进行平均分组或不平均分组;
对每个等效线声源组中测点的噪声声强级合成值和测量表面面积的附加声强级,进行叠加计算,获取各个等效线声源的声功率。
其中,预设条件为声功率相差120°的相位。
其中,噪声预测点具有不同方位,在不同方位内的预测点,分别确定相应方位内的预测点,使用统一坐标系,以三相电抗器A相或C相的任一底角为坐标原点,确定每个线声源及预测点的坐标。
下面结合实施例对本发明进行进一步的说明:
如图2所示,步骤101,电抗器拆分单元根据电抗器的模型结构,将其进行井型立体拆分,获取线声源的长短和位置;
根据电抗器结构模块,在垂直方向三等分,在长水平方向三等分,在宽水平方向二等分,按等分线进行切分,将每相电抗器拆分成18个长方体,并获得45条线声源。
步骤102,测点测量单元在获取的线声源处每隔0.5m设置一个测点,对测点处的噪声声强或声压进行测量,以获得每个测点处的声强值或声压值,每个侧面具有多个测点。
将测量轮廓线布置在电抗器基准发射机面外1m,采用声级计或声强计分别测量线声源处各测点的声压或声强,计算每个点的声功率。
步骤103,声功率级计算单元根据声功率级计算方法,分别对多组测点的噪声数据进行计算与统计,获取每段线声源线段上得到相应的声功率。
将测得的与每条线所对应的多个测点的噪声数据根据等效线声源的个数进行平均分组,或者将测得的与每条线所对应的多个测点的噪声数据根据所述等效线声源个数进行不平均分组。
例如,选取的某个面的等效线声源个数为12,将所测得的48个测点对的噪声数据不平均分为12组,其中水平方向平均分为6组,每组包含5个测点的噪声数据,分组情况为1/3高度处的左边5个测点为一组,中间5个测点为一组,右边5个测点为一组,2/3高度处的左边5个测点为一组,中间5个测点为一组,右边5个测点为一组;垂直方向平均分为6组,每组包含3个测点,分组情况为左上方3个测点为一组,左侧中间3个测点为一组,左下方3个测点为一组,右上方3个测点为一组,右侧中间3个测点为一组,右下方3个测点为一组。
测量表面积根据S=(H+2)*L计算,H为电抗器高度,L是测量轮廓线的长度。
例如,利用分得的5组噪声数据分别计算相应的线声源声功率级,电抗器高度为h=4.5m,每组的测量轮廓线为L=3m,获得1个点声源的声功率级。
根据每个等效线声源组中测点的噪声声强级合成值和测量表面面积的附加声强级,分别叠加计算得到各个等效线声源的声功率。
步骤104,AFEM模型建立单元基于声学有限元(AFEM)理论,建立针对三相电抗器井型立体线声源等效模型,分别在三相电抗器等效线声源处输入三相声功率。
将计算得到的声功率以相差120°的相位将其分别输入到三相高抗的有限元声源模型中。
步骤105,SFEM模型建立与FEM声固耦合单元基于固体力学有限元(SFEM)理论,建立防火墙有限元模型,并根据FEM声固耦合方法,设置声固耦合边界条件。
保证每相高抗两侧均有防火墙,共设置4座防火墙,并添加声固耦合边界条件。
步骤106,BEM模型建立单元基于边界元(BEM)理论,建立边界元模型,确定预测点及预测范围。
根据实际变电站占地面积及所需预测噪声声场体积,分别在x,y,z轴定义坐标、最小值与最大值,从而确定预测点与预测范围。
步骤107,噪声计算单元基于FEM-BEM耦合方法,使用统一坐标系,根据声学波动方程及波动方程对预测点及预测范围处噪声进行计算。优选地,对于不同方位内的预测点,分别确定相应方位内的预测点,使用统一坐标系,以三相电抗器A相或C相的某一底角为坐标原点,即i=0(i=x,y,z),确定每个线声源及预测点的坐标,基于FEM声固耦合及FEM-BEM耦合方法,根据声学波动方程及波波动方程对预测点处噪声进行计算。
步骤108,后处理单元获取边界元范围内的声压和声压级三维立体数据并绘制二维多切面云图,预测近、远场噪声分布情况。
根据所需预测的具体情况,分别设置x,y,z轴的值,确定声压和声压级切片的位置。
如图3所示,图3中说明了根据本发明中立体井型线声源在电抗器上的分布示意图;
图3所示高抗的西侧和南侧,取高抗分析侧面垂直中轴线上的点作为计算与对比分析点,建立统一坐标系,计算各噪声衰减分量后,得到预测点处噪声计算与实测对比结果如表1和表2所示,表1为高抗西面噪声衰减预测结果,表2为高抗南面噪声衰减预测结果。
表1
表2
距变压器外边界距离(m) 仿真计算结果(dB) 实测数据(dB) 误差(dB)
1 95.1 95.9 -0.8
3 92.3 93.6 -1.3
5 91.0 91.5 -0.5
7 90.3 91.3 -1.0
9 89.5 88.0 1.5
11 88.8 88.1 0.7
13 88.3 88.3 0.0
15 87.8 85.0 2.8
17 87.3 85.4 1.9
19 87.0 84.3 2.7
21 86.8 84.9 1.9
23 86.7 85.1 1.6
25 86.6 84.9 1.7
通过上述表1和表2可以看出,本发明的声源模型计算得到的高抗噪声数据与实测结果差距控制在3dB以内,误差较小。
如图4所示,通过声压级云图可以明显观察到变电站噪声的传播具有指向性,声压级在各个方向上都呈现出非对称分布,且传播更加参差,极大与极小值的差值更大。其中,在xy平面的传播呈现出“高X”型分布特点。
如图5所示,通过声压云图可以明显看到,声波的干涉效应在FEM声固耦合及FEM-BEM耦合方法下得到很好的体现。
本发明还提出了一种基于三相电抗器声源模型的噪声预测系统200,如图6所示,包括:
测量模块201,对三相电抗器进行井型立体拆分,获取线声源,并确定线声源的长短和位置;针对线声源,根据线声源的长短及位置,以预设的间隔设置测点,对测点处的噪声声强或声压进行测量,获取每个测点处的声强值或声压值;
计算模块202,根据声强值或声压值,确定线声源的声功率;
模型搭建模块203,搭建三相电抗器声源模型的井型立体线等效模型,并将符合预设条件的声功率输入等效模型;搭建防火墙有限元模型,并针对防火墙有限元模型设置声固耦合边界条件;搭建边界元模型,并确定模型的噪声预测点及预测范围;
预测模块204,针对输入声功率的等效模型、防火墙有限元模型及边界元模型,使用统一坐标系,根据声学波动方程及波动方程对噪声预测点及预测范围进行噪声计算,获取边界元范围内的声压数据和声压级三维立体数据,对声压数据和声压级三维立体数据绘制二维多切面云图,根据二维多切面云图,预测线声源的近、远场噪声的分布情况。
其中,对三相电抗器进行井型立体拆分,包括:
根据电抗器结构,在电抗器垂直方向三等分,在长水平方向三等分,在宽水平方向二等分,按等分线进行切分,将每相电抗器拆分成18个长方体,并获得45条线声源。
其中,根据声强值或声压值,确定线声源的声功率,具体为:
对获取的每个测点处的声强值或声压值,根据等效线声源的个数进行平均分组或不平均分组;
对每个等效线声源组中测点的噪声声强级合成值和测量表面面积的附加声强级,进行叠加计算,获取各个等效线声源的声功率。
其中,预设条件为声功率相差120°的相位。
其中,噪声预测点具有不同方位,在不同方位内的预测点,分别确定相应方位内的预测点,使用统一坐标系,以三相电抗器A相或C相的任一底角为坐标原点,确定每个线声源及预测点的坐标。
本发明弥补了常用商业软件无法计算干涉效应的缺点,同时考虑到三相电抗器的三个相各不相同且存在防火墙,将高抗的声功率设置三种不同的相位,对不同相位的高抗之间设立声固耦合边界,同时考虑到声场的三维性,对每个相的高抗都采用立体井型线声源进行建立,使仿真噪声具有立体性,具有更高的真实还原度及预测精度。
本发明对声源模型的精准化建立及声传播预测提供了新的理论依据与方法指导,利用本发明,可以用建立三相电抗器精准的声源模型,从而为声传播预测与新建变电站方案设计提供准确的仿真数据。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。本申请实施例中的方案可以采用各种计算机语言实现,例如,面向对象的程序设计语言Java和直译式脚本语言JavaScript等。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
尽管已描述了本申请的优选实施例,但本领域内的技术人员一旦得知了基本创造性概念,则可对这些实施例作出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本申请范围的所有变更和修改。
显然,本领域的技术人员可以对本申请进行各种改动和变型而不脱离本申请的精神和范围。这样,倘若本申请的这些修改和变型属于本申请权利要求及其等同技术的范围之内,则本申请也意图包含这些改动和变型在内。

Claims (10)

1.一种基于三相电抗器声源模型的噪声预测方法,所述方法包括:
对三相电抗器进行井型立体拆分,获取线声源,并确定线声源的长短和位置;针对线声源,根据线声源的长短及位置,以预设的间隔设置测点,对测点处的噪声声强或声压进行测量,获取每个测点处的声强值或声压值;
根据声强值或声压值,确定线声源的声功率;
搭建三相电抗器声源模型的井型立体线等效模型,并将符合预设条件的声功率输入等效模型;搭建防火墙有限元模型,并针对防火墙有限元模型设置声固耦合边界条件;搭建边界元模型,并确定模型的噪声预测点及预测范围;
针对输入声功率的等效模型、防火墙有限元模型及边界元模型,使用统一坐标系,根据声学波动方程及波动方程对噪声预测点及预测范围进行噪声计算,获取边界元范围内的声压数据和声压级三维立体数据,对声压数据和声压级三维立体数据绘制二维多切面云图,根据二维多切面云图,预测线声源的近、远场噪声的分布情况。
2.根据权利要求1所述的方法,所述对三相电抗器进行井型立体拆分,包括:
根据电抗器结构,在电抗器垂直方向三等分,在长水平方向三等分,在宽水平方向二等分,按等分线进行切分,将每相电抗器拆分成18个长方体,并获得45条线声源。
3.根据权利要求1所述的方法,所述根据声强值或声压值,确定线声源的声功率,具体为:
对获取的每个测点处的声强值或声压值,根据等效线声源的个数进行平均分组或不平均分组;
对每个等效线声源组中测点的噪声声强级合成值和测量表面面积的附加声强级,进行叠加计算,获取各个等效线声源的声功率。
4.根据权利要求1所述的方法,所述预设条件为声功率相差120°的相位。
5.根据权利要求1所述的方法,所述噪声预测点具有不同方位,在不同方位内的预测点,分别确定相应方位内的预测点,使用统一坐标系,以三相电抗器A相或C相的任一底角为坐标原点,确定每个线声源及预测点的坐标。
6.一种基于三相电抗器声源模型的噪声预测系统,所述系统包括:
测量模块,对三相电抗器进行井型立体拆分,获取线声源,并确定线声源的长短和位置;针对线声源,根据线声源的长短及位置,以预设的间隔设置测点,对测点处的噪声声强或声压进行测量,获取每个测点处的声强值或声压值;
计算模块,根据声强值或声压值,确定线声源的声功率;
模型搭建模块,搭建三相电抗器声源模型的井型立体线等效模型,并将符合预设条件的声功率输入等效模型;搭建防火墙有限元模型,并针对防火墙有限元模型设置声固耦合边界条件;搭建边界元模型,并确定模型的噪声预测点及预测范围;
预测模块,针对输入声功率的等效模型、防火墙有限元模型及边界元模型,使用统一坐标系,根据声学波动方程及波动方程对噪声预测点及预测范围进行噪声计算,获取边界元范围内的声压数据和声压级三维立体数据,对声压数据和声压级三维立体数据绘制二维多切面云图,根据二维多切面云图,预测线声源的近、远场噪声的分布情况。
7.根据权利要求6所述的系统,所述对三相电抗器进行井型立体拆分,包括:
根据电抗器结构,在电抗器垂直方向三等分,在长水平方向三等分,在宽水平方向二等分,按等分线进行切分,将每相电抗器拆分成18个长方体,并获得45条线声源。
8.根据权利要求6所述的系统,所述根据声强值或声压值,确定线声源的声功率,具体为:
对获取的每个测点处的声强值或声压值,根据等效线声源的个数进行平均分组或不平均分组;
对每个等效线声源组中测点的噪声声强级合成值和测量表面面积的附加声强级,进行叠加计算,获取各个等效线声源的声功率。
9.根据权利要求6所述的系统,所述预设条件为声功率相差120°的相位。
10.根据权利要求6所述的系统,所述噪声预测点具有不同方位,在不同方位内的预测点,分别确定相应方位内的预测点,使用统一坐标系,以三相电抗器A相或C相的任一底角为坐标原点,确定每个线声源及预测点的坐标。
CN202011517182.1A 2020-12-21 2020-12-21 一种基于三相电抗器声源模型的噪声预测方法及系统 Active CN112685928B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011517182.1A CN112685928B (zh) 2020-12-21 2020-12-21 一种基于三相电抗器声源模型的噪声预测方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011517182.1A CN112685928B (zh) 2020-12-21 2020-12-21 一种基于三相电抗器声源模型的噪声预测方法及系统

Publications (2)

Publication Number Publication Date
CN112685928A CN112685928A (zh) 2021-04-20
CN112685928B true CN112685928B (zh) 2024-01-30

Family

ID=75449661

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011517182.1A Active CN112685928B (zh) 2020-12-21 2020-12-21 一种基于三相电抗器声源模型的噪声预测方法及系统

Country Status (1)

Country Link
CN (1) CN112685928B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113283131B (zh) * 2021-04-22 2022-10-28 杭州申弘智能科技有限公司 一种适用于变电站的火势蔓延预测方法
CN113221413B (zh) * 2021-05-10 2023-02-17 国网湖北省电力有限公司电力科学研究院 一种变电站可听噪声预测的声波干涉修正方法及系统
CN113804285B (zh) * 2021-05-13 2024-04-05 中国电力科学研究院有限公司 确定电力变压器组低频谐波噪声传播声场的方法及系统
CN113656997B (zh) * 2021-07-09 2024-02-27 中国电力科学研究院有限公司 一种用于变压器油箱结构低噪声优化的方法及系统
CN114216635B (zh) * 2021-11-24 2023-08-22 广西电网有限责任公司电力科学研究院 并联电抗器拉紧螺杆影响振动和噪声的试验装置、方法
CN116151169B (zh) * 2022-11-30 2024-03-12 西安西电变压器有限责任公司 预测电抗器噪声强度的方法和装置
CN118296914A (zh) * 2024-05-29 2024-07-05 邢台职业技术学院 一种基于fem-bem点线复合声源模型的调相机厂房噪声预测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002365127A (ja) * 2001-06-06 2002-12-18 Mitsubishi Heavy Ind Ltd 音源の音響パワーレベル推定方法および音源の音響パワーレベル推定装置
CN106096184A (zh) * 2016-06-24 2016-11-09 中国电力科学研究院 一种基于声强法的变压器多点声源模型的噪声计算方法及系统
CN110319920A (zh) * 2019-05-15 2019-10-11 中国船舶重工集团公司第七一九研究所 一种基于设备等效声源点的噪声定量预报方法及系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002365127A (ja) * 2001-06-06 2002-12-18 Mitsubishi Heavy Ind Ltd 音源の音響パワーレベル推定方法および音源の音響パワーレベル推定装置
CN106096184A (zh) * 2016-06-24 2016-11-09 中国电力科学研究院 一种基于声强法的变压器多点声源模型的噪声计算方法及系统
CN110319920A (zh) * 2019-05-15 2019-10-11 中国船舶重工集团公司第七一九研究所 一种基于设备等效声源点的噪声定量预报方法及系统

Also Published As

Publication number Publication date
CN112685928A (zh) 2021-04-20

Similar Documents

Publication Publication Date Title
CN112685928B (zh) 一种基于三相电抗器声源模型的噪声预测方法及系统
CN106096184B (zh) 一种基于声强法的变压器多点声源模型的噪声计算方法及系统
CN111159928B (zh) 一种基于多线声源模型的变压器噪声计算方法及系统
Gamper Head-related transfer function interpolation in azimuth, elevation, and distance
CN102194027B (zh) 变电站噪声控制优化分析方法及系统
WO2022237708A1 (zh) 一种变电站可听噪声预测的声波干涉修正方法及系统
CN105223480B (zh) 天线阵列时差法定位变电站局部放电源的定位误差仿真方法
WO2012172670A1 (ja) 電波伝搬環境計測装置、無線ネットワーク構築システムおよび電波伝搬環計測方法
CN106646645B (zh) 一种重力正演加速方法
CN104811968B (zh) 一种传播模型的校正方法及装置
CN110139359B (zh) 一种干扰源定位处理方法及装置
CN111261188A (zh) 高压变压器噪声频谱确定方法及装置
CN105118091A (zh) 一种构建多精度非均匀地质网格曲面模型的方法和系统
CN108242962B (zh) 一种基于测量报告的室内信号传播损耗计算方法及装置
CN112818553A (zh) 一种基于等效源的变压器辐射声场计算方法
JP2018032939A (ja) 品質推定装置及び品質推定方法
CN106658538B (zh) 一种基于泰森多边形的手机基站信号覆盖范围模拟方法
US20140018084A1 (en) Method and apparatus for antenna placement of wireless basestations
CN112230277B (zh) 基于柱坐标系的隧道地震波传播数值模拟方法及系统
CN107995630B (zh) 一种无线网络覆盖的信号数据计算方法及装置
CN113537724A (zh) 一种电网设施抗震安全风险评估方法及装置
CN107566051A (zh) 一种mimo ota最大三维测试区域大小的确定方法及装置
CN117725769A (zh) 基于多尺度有限元的复杂模型航空电磁三维快速正演方法
CN113312759A (zh) 一种计算换流变压器组噪声源强的方法及系统
CN108154557A (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