CN113221413B - 一种变电站可听噪声预测的声波干涉修正方法及*** - Google Patents
一种变电站可听噪声预测的声波干涉修正方法及*** Download PDFInfo
- Publication number
- CN113221413B CN113221413B CN202110508682.7A CN202110508682A CN113221413B CN 113221413 B CN113221413 B CN 113221413B CN 202110508682 A CN202110508682 A CN 202110508682A CN 113221413 B CN113221413 B CN 113221413B
- Authority
- CN
- China
- Prior art keywords
- sound
- transformer substation
- interference
- equipment
- substation
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 76
- 238000012937 correction Methods 0.000 title claims abstract description 34
- 238000004364 calculation method Methods 0.000 claims abstract description 60
- 230000005855 radiation Effects 0.000 claims abstract description 27
- 238000004088 simulation Methods 0.000 claims abstract description 23
- 238000009826 distribution Methods 0.000 claims description 28
- 238000004422 calculation algorithm Methods 0.000 claims description 18
- 239000003990 capacitor Substances 0.000 claims description 18
- 238000006073 displacement reaction Methods 0.000 claims description 18
- 230000000694 effects Effects 0.000 claims description 12
- 230000002068 genetic effect Effects 0.000 claims description 12
- 238000005259 measurement Methods 0.000 claims description 12
- 239000012088 reference solution Substances 0.000 claims description 6
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000011076 safety test Methods 0.000 claims description 3
- 239000007787 solid Substances 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 2
- 229910052739 hydrogen Inorganic materials 0.000 claims description 2
- 239000001257 hydrogen Substances 0.000 claims description 2
- 125000004435 hydrogen atom Chemical class [H]* 0.000 claims description 2
- 238000009472 formulation Methods 0.000 abstract 1
- 239000000203 mixture Substances 0.000 abstract 1
- 230000006870 function Effects 0.000 description 11
- 238000010586 diagram Methods 0.000 description 9
- 238000004590 computer program Methods 0.000 description 7
- 238000012360 testing method Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 5
- 238000001816 cooling Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000001427 coherent effect Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 230000001755 vocal effect Effects 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 239000010977 jade Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
Images
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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/10—Noise analysis or noise optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明提供一种变电站可听噪声预测的声波干涉修正方法及***,包括a、获取变电站主设备声源参数,以及获取变电站设备、建筑物布局的空间几何尺寸,并构建变电站三维仿真模型用于声场计算;b、设置对应的声环境参数及边界条件,仿真计算出不考虑干涉因素的变电站厂界噪声;c、计算每台变电站设备的辐射声场;d、利用区间随机模型摄动法,获取单独设备的相对相位差,及其对应的单独设备的干涉声纹;e、采用步骤d中获取的干涉声纹,对步骤b中求取的不考虑干涉因素的变电站厂界噪声给予干涉修正,得到变电站厂界噪声。本发明有助于指导变电站厂界噪声预测计算和厂界达标治理工作,为噪声技改方案的制定提供算法上的支持。
Description
技术领域
本发明涉及变电站可听噪声预测领域,具体是一种变电站可听噪声预测的声波干涉修正方法及***。
背景技术
变电站环境噪声的计算主要依据ISO9613-2规定的通用工程算法,以及基于该算法开发的仿真软件,如Soundplan等。这些软件在对强干涉声源的辐射噪声预测时,并未计及相干声源的低频谐波特征和声波干涉现象。基于上述软件的预测误差普遍在5dB左右,甚至有多个误差10dB以上的超标厂界,给设计部门出具治理方案带来了很大困扰。
目前,能够考虑干涉现象的声场计算方法有以边界元法为代表的LMS VirtualLab、以圆锥束法和三棱锥束法为基础的Raynoise,但是,在变电站噪声预测这种大尺度范围的声场计算场合,难以兼顾计算量和计算精度。
在变电站噪声的预测方面,当前成果多集中于变压器、电抗器以及电容器的单独声源辐射特征的仿真计算与实测研究,例如周兵、王玉兴等人利用振动法-边界元法对高压电抗器的辐射声功率进行了计算求解;张鸣等人通过实验测量分析变电站的噪声特征及分布规律;林丽妲等人通过采用一种基于双球面法且测点绕垂直地面半径所在轴旋转测量的方式,设计了电容器声功率的测试方法。
对于变电站厂界上的噪声预测,实质上是针对户外声环境的声场计算问题。在该特定声环境中,主要声源为变压器、电抗器、电容器等电力设备。变电站的地表状况和其他建筑则在声波的传播过程中起到吸收、反射、折射的作用。最终厂界上某敏感点的声压可以看作是各声源产生的声波在经历传播衰减后,于敏感点处的能量叠加,势必需要考虑声源的干涉效应。
发明内容
本发明提出一种变电站可听噪声预测的声波干涉修正方法及***,提高变电站厂界噪声的预测精度。
本发明的技术方案:
一种变电站可听噪声预测的声波干涉修正方法,包括:
a、获取变电站主设备声源参数,以及获取变电站设备、建筑物布局的空间几何尺寸,并构建变电站三维仿真模型用于声场计算;
b、以a中获取的变电站主设备声源参数中的声功率为源强参数,设置对应的声环境参数及边界条件,仿真计算出不考虑干涉因素的变电站厂界噪声;
c、根据a中获取的变电站主设备声源参数中主设备表面振动或近场声压,计算每台变电站设备的辐射声场;
d、根据c中每台变电站设备的辐射声场,利用区间随机模型摄动法,获取单独设备的相对相位差,及其对应的单独设备的干涉声纹;
e、采用步骤d中获取的干涉声纹,对步骤b中求取的不考虑干涉因素的变电站厂界噪声给予干涉修正,得到变电站厂界噪声。
所述步骤a中实测获取安全区域内设备表面振动或近场声压数据,对设备安全区域外进行振动参数反演方法获取设备声源参数,其中所述的安全区域振动、声压数据为变电站主设备2/3高度以下的安全测试范围内获取的变压器、电抗器表面振动数据、近场声压数据。
所述步骤a中对设备安全区域外进行振动参数反演方法获取设备声源参数具体是采用安全区域内的设备表面振动或近场声压,预测无法测量的安全区域外的设备表面振动或近场声压数据,具体为:
根据设备图纸或实测结构尺寸,采用有限元法计算复杂表面的模态向量Φn(x),获取模态导纳Yn(ω),n=1,2,...,N,在M点的振动位移PSD谱和相对相位通过现场设备的表面振动测量得到:
[W]M×1=[W(x1),W(x2),...,W(xM)]T
然后,求解模态力向量f(x):
变电站设备的模态力展开为:
模态位移展开为:
主变、高抗顶部套管位置,以及侧面安全范围以外的振动位移表示为:
模态位移的重建精度取决于特征矩阵与实测已知点振动位移的准确性,采用直接法求解模态位移:
采用结构有限元法计算上述模态力激发的变电站主设备振动分布,采用ABAQUS程序,将模态力考虑到变压器、电抗器设备的流固耦合模型中,计算外壳表面的振动速度。
所述步骤a中对设备安全区域外进行振动参数反演方法获取设备声源参数还包括电容器等效声源的参数,电容器等效声源参数的反演采用现场测量物理可观声压数据,以遗传算法,作为演化算法进行群体搜索,使描述可观测物理量的计算值与观测值差的目标函数最小,具体为:
电容器塔的每层电容器看成n2列×n1行源强相等、相位相同的等效点声源的组合,总声压p为:
式中:A为距离每个等效声源1米处的声压幅值;ri为受声点和声源的距离,单位为m;ω为圆角频率,单位为rad/s;k为声波的波数,单位为m-1,
各个等效点声源的相位表示为:
ri=r1+(i1-1)l1 sinθ+(i2-1)l2 cosθ
其中,1≤i1≤n1,1≤i2≤n2,
等效源的总声压可写为:
在使用遗传算法进行等效点源的相位获取时,构造目标函数为实测2D声辐射指向性和总声压计算得到的声指向性之差f(Φ),遗传算法的搜索空间([Ql,Qr])将在获取初始数据后给出,在该空间中随机产生Np个种群,并计算每个点声源的目标函数值,选择适应度最大的个体作为全局最优个体,即为电容器等效声源的参数。
所述步骤a中的三维变电站仿真模型为DXF格式的地形文件和STP格式的CAD模型,这两种模型分别被用于步骤b中Soundplan软件的无干涉效应声场计算和步骤c中LMSVirtual lab软件的干涉声场计算;所述步骤b中采用的仿真软件是基于Soundplan的噪声仿真软件,所述步骤c中单独设备声场的计算采用LMS Virtual lab软件进行,在单独设备的声场计算时,已知边界取自步骤a中主设备的表面振动或近场声压数据;变电站数模取自步骤a中的CAD模型。
所述步骤d中单独设备的相对相位差是指不同单独设备之间存在的相位差,同一个单独设备的相位关系是确定的,且已经在骤a中获取。
所述步骤d中区间随机模型摄动法的具体实施步骤,
根据每台变电站设备的辐射声场,将等效声源的负荷波动、等效位置误差作为区间变量;
采用随机区间摄动法计算等效声源相位的期望与标准差变化范围;以现场测量的变电站主设备辐射声场观察点的实际相位数据作为参考解,验证评估随机区间摄动法的计算结果;
利用单独声源相位数据的参考解,调整每个声源独自运行时的声压分布,得到对应的单独设备的干涉声纹。
所述步骤e中干涉修正的步骤是:
将相位调整后的每个设备的干涉声纹归一化;
将归一化后的单独设备的干涉声纹加权求和,得到考虑了干涉效应的声压分布模式;
把步骤b中得到的不考虑干涉因素的变电站厂界噪声换算成声波能量后,减去考虑了干涉效应的声压分布模式,得到综合修正后的变电站噪声分布和厂界敏感点噪声,完成变电站厂界噪声的修正。
一种变电站可听噪声预测的声波干涉修正***,包括主设备声源参数获取模块,变电站设备、建筑物布局的空间几何尺寸获取模块,变电站三维仿真模型构建模块,不考虑干涉因素的变电站厂界噪声计算模块,辐射声场计算模块,干涉声纹获取模块,变电站厂界噪声生成模块,
所述主设备声源参数获取模块采用实测安全区域振动、声压数据,对设备其他区域进行振动参数反演的方法,获取安全区域外的设备表面振动或近场声压数据;
所述变电站设备、建筑物布局的空间几何尺寸获取模块用以获取变电站设备、建筑物布局的空间几何尺寸;
所述变电站三维仿真模型构建模块根据变电站设备、建筑物布局的空间几何尺寸获取模块获取的变电站设备、建筑物布局的空间几何尺寸构建变电站三维仿真模型;
所述不考虑干涉因素的变电站厂界噪声计算模块以设备声源参数模型获取模块的变电站主设备声功率为源强参数,设置对应的声环境参数及边界条件,仿真计算不考虑干涉因素的变电站厂界噪声;
所述辐射声场计算模块以设备声源参数模型获取模块的主设备表面振动或近场声压作为已知条件,单独计算每台变电站设备的辐射声场;
所述干涉声纹获取模块利用区间随机模型摄动法,获取单独设备的相对相位差,及其对应的单独设备的干涉声纹;
所述变电站厂界噪声生成模块采用干涉声纹获取模块获取的干涉声纹,对不考虑干涉因素的变电站厂界噪声计算模块中基于ISO9613-2下风向声压级计算法则求取的变电站噪声分布给予干涉修正,得到最终的变电站厂界噪声。
所述声波干涉修正***还包括安全区域振动、声压数据的测量装置。
与现有技术相比,本发明的有益效果是:通过考虑变电站设备辐射噪声的干涉效应,修正基于ISO9613-2下风向声压级计算法则求取的变电站噪声分布,提高变电站厂界噪声的预测精度。该方法可以在主要声源设备的表面或近场测量安全范围内的少数样本测点,反演不可测量设备区域的声源参数,在保证测试安全的前提下,实现声源的最佳等效模型。
附图说明
图1为本发明的方法流程图。
图2为考虑干涉效应的修正效果图。
图3为不考虑干涉效应的声场结果图。
图4为本发明的***框图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
声波在变电站内传播会遇到防火墙、站界围墙和站内建筑物的遮挡,因而产生反射、透射、折射等声学现象。通常噪声源辐射出的声波到达接收点,不只是由声源经过空气传播直接到达接收点,还可以经过地面反射后到达接收点,经楼房、围墙、隔声屏等建筑物后产生的反射与衍射到达接收点。本专利提出采用实测站内声学参照物的方法,建立站内声场的声传播模型。要求在变电站内部的建筑物、变电站功能性结构、围墙、防火墙等都需要实际测量其长度、宽度、高度以及相对位置,以保证计算模型的准确性。测量结束后,开展变电站模型建模工作,需要建立三维实体模型、有限元模型以及边界元计算模型;有限元模型建模时需要按照声学边界元计算规则确定网格尺寸。
特高压主变压器和高压电抗器表面的振动测量按国标《GBT 16539-1996声学振速法测定噪声源声功率级用于封闭机器的测量》的要求布置测点;测试时至少采用2个传感器。要求必须有1个传感器作为参考传感器,位于被测点中部,另外几个传感器用于表面振动的测量,每相测试都是由左侧顶部位置开始从上往下依次测量。测试时既需要记录测点的振动加速度幅值还需记录各测点与参考点的相位关系。测试结果作为计算边界输入到声学边界元预测模型中,开展变电站主设备干涉声场预测。
利用区间随机模型摄动法,获取相干声源相位、空间指向性等涉及变电站声场计算的关键参数,对ISO9613-2下风向声压级计算法则求取的变电站噪声分布给予适当的干涉修正。具体操作过程如下
首先,实际测量变电站主设备的近场噪声,利用等效声源法求解变电站内各设备的声功率,以声功率为计算输入,通过国标ISO9613-2下风向声压级计算法求解变电站内部声场分布。利用过程3介绍的分析方法,以遗传算法(Genetic Algorithm)作为演化算法进行群体搜索,使描述可观测物理量的计算值与观测值差的目标函数最小,从而获得变压器、电容器、冷却风机等主要设备声源的真实模型参数,开展变电站声场预测计算。
然后,实际测量变电站各主设备的表面振动加速度,同时利用边界元法建模,求解各主设备的近场声纹分布结果,并对计算结果开展归一化处理;实际测量各声源附近样本区域内的声场分布,对比样本区间声场分布计算结果与实测结果,确定声场分布的加权因子。
最后,利用归一化后的声波干涉特性及加权因子,对变电站声场分布特性开展干涉修正,完成变电站厂区干涉声场分布计算。
请参阅图1-3,一种变电站可听噪声预测的声波干涉修正方法,包括:
a、获取变电站主设备声源参数,以及获取变电站设备、建筑物布局的空间几何尺寸,并构建变电站三维仿真模型用于声场计算;
b、以a中获取的变电站主设备声源参数中的声功率为源强参数,设置对应的声环境参数及边界条件,仿真计算出不考虑干涉因素的变电站厂界噪声;
c、根据a中获取的变电站主设备声源参数中主设备表面振动或近场声压,计算每台变电站设备的辐射声场;
d、根据c中每台变电站设备的辐射声场,利用区间随机模型摄动法,获取单独设备的相对相位差,及其对应的单独设备的干涉声纹;
e、采用步骤d中获取的干涉声纹,对步骤b中求取的不考虑干涉因素的变电站厂界噪声给予干涉修正,得到变电站厂界噪声。
所述步骤a中实测获取安全区域内设备表面振动或近场声压数据,对设备安全区域外进行振动参数反演方法获取设备声源参数,其中所述的安全区域振动、声压数据为变电站主设备2/3高度以下的安全测试范围内获取的变压器、电抗器表面振动数据、近场声压数据。
所述步骤a中对设备安全区域外进行振动参数反演方法获取设备声源参数具体是采用安全区域内的设备表面振动或近场声压,预测无法测量的安全区域外的设备表面振动或近场声压数据,具体为:
根据设备图纸或实测结构尺寸,采用有限元法计算复杂表面的模态向量Φn(x),获取模态导纳Yn(ω),n=1,2,...,N,在M点的振动位移PSD谱和相对相位通过现场设备的表面振动测量得到:
[W]M×1=[W(x1),W(x2),...,W(xM)]T
然后,求解模态力向量f(x):
变电站设备的模态力展开为:
模态位移展开为:
主变、高抗顶部套管位置,以及侧面安全范围以外的振动位移表示为:
模态位移的重建精度取决于特征矩阵与实测已知点振动位移的准确性,采用直接法求解模态位移:
采用结构有限元法计算上述模态力激发的变电站主设备振动分布,采用ABAQUS程序,将模态力考虑到变压器、电抗器设备的流固耦合模型中,计算外壳表面的振动速度。
所述步骤a中的振动参数反演方法还包括电容器、冷却风机等效声源的参数建模,电容器等效声源参数的反演采用现场测量物理可观声压数据,以遗传算法,作为演化算法进行群体搜索,使描述可观测物理量的计算值与观测值差的目标函数最小,具体为:
电容器塔的每层电容器看成n2列×n1行源强相等、相位相同的等效点声源的组合,总声压p为:
式中:A为距离每个等效声源1米处的声压幅值;ri为受声点和声源的距离,单位为m;ω为圆角频率,单位为rad/s;k为声波的波数,单位为m-1,
各个等效点声源的相位表示为:
ri=r1+(i1-1)l1 sinθ+(i2-1)l2 cosθ
其中,1≤i1≤n1,1≤i2≤n2,
等效源的总声压可写为:
在使用遗传算法进行等效点源的相位获取时,构造目标函数为实测2D声辐射指向性和总声压计算得到的声指向性之差f(Φ),遗传算法的搜索空间([Ql,Qr])将在获取初始数据后给出,在该空间中随机产生Np个种群,并计算每个点声源的目标函数值,选择适应度最大的个体作为全局最优个体。
冷却风机作为变压器、电抗器的噪声附件并不单独运行。风机噪声的宽频特性使得它在设备本体噪声的谐波频率附近具有相对较小的幅值。在本专利中,将其用等效声功率近似为与风机个数相同的点声源组合。而各等效点源在工频谐波处的声波相位信息则由类似的遗传算法确定。
所述步骤a中的三维变电站仿真模型为DXF格式的地形文件和STP格式的CAD模型,这两种模型分别被用于步骤b中Soundplan软件的无干涉效应声场计算和步骤c中LMSVirtual lab软件的干涉声场计算;所述步骤b中采用的仿真软件是基于Soundplan的噪声仿真软件,所述步骤c中单独设备声场的计算采用LMS Virtual lab软件进行,在单独设备的声场计算时,已知边界取自步骤a中主设备的表面振动或近场声压数据;变电站数模取自步骤a中的CAD模型。
所述步骤d中单独设备的相对相位差是指不同单独设备之间存在的相位差,同一个单独设备的相位关系是确定的,且已经在骤a中获取。
所述步骤d中区间随机模型摄动法的具体实施步骤,
根据每台变电站设备的辐射声场,将等效声源的负荷波动、等效位置误差作为区间变量;
采用随机区间摄动法计算等效声源相位的期望与标准差变化范围;以现场测量的变电站主设备辐射声场观察点的实际相位数据作为参考解,验证评估随机区间摄动法的计算结果;
利用单独声源相位数据的参考解,调整每个声源独自运行时的声压分布,得到对应的单独设备的干涉声纹。
所述步骤e中干涉修正的步骤是:
将相位调整后的每个设备的干涉声纹归一化;
将归一化后的单独设备的干涉声纹加权求和,得到考虑了干涉效应的声压分布模式;
把步骤b中得到的不考虑干涉因素的变电站厂界噪声换算成声波能量后,减去考虑了干涉效应的声压分布模式,得到综合修正后的变电站噪声分布和厂界敏感点噪声,完成变电站厂界噪声的修正。
请参阅图4,一种变电站可听噪声预测的声波干涉修正***,包括主设备声源参数获取模块1,变电站设备、建筑物布局的空间几何尺寸获取模块2,变电站三维仿真模型构建模块3,不考虑干涉因素的变电站厂界噪声计算模块4,辐射声场计算模块5,干涉声纹获取模块6,变电站厂界噪声生成模块7,
所述主设备声源参数获取模块1采用实测安全区域振动、声压数据,对设备其他区域进行振动参数反演的方法,获取安全区域外的设备表面振动或近场声压数据;
所述变电站设备、建筑物布局的空间几何尺寸获取模块2用以获取变电站设备、建筑物布局的空间几何尺寸;
所述变电站三维仿真模型构建模块3根据变电站设备、建筑物布局的空间几何尺寸获取模块2获取的变电站设备、建筑物布局的空间几何尺寸构建变电站三维仿真模型;
所述不考虑干涉因素的变电站厂界噪声计算模块4以设备声源参数模型获取模块1的变电站主设备声功率为源强参数,设置对应的声环境参数及边界条件,仿真计算不考虑干涉因素的变电站厂界噪声;
所述辐射声场计算模块5以设备声源参数模型获取模块1的主设备表面振动或近场声压作为已知条件,单独计算每台变电站设备的辐射声场;
所述干涉声纹获取模块6利用区间随机模型摄动法,获取单独设备的相对相位差,及其对应的单独设备的干涉声纹;
所述变电站厂界噪声生成模块7采用干涉声纹获取模块6获取的干涉声纹,对不考虑干涉因素的变电站厂界噪声计算模块4中基于ISO9613-2下风向声压级计算法则求取的变电站噪声分布给予干涉修正,得到最终的变电站厂界噪声。
所述声波干涉修正***还包括安全区域振动、声压数据的测量装置。
本领域内的技术人员应明白,本申请的实施例可提供为方法、***、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求保护范围之内。
Claims (7)
1.一种变电站可听噪声预测的声波干涉修正方法,其特征在于,包括:
a、获取变电站主设备声源参数,以及获取变电站设备、建筑物布局的空间几何尺寸,并构建变电站三维仿真模型用于声场计算;
b、以a中获取的变电站主设备声源参数中的声功率为源强参数,设置对应的声环境参数及边界条件,仿真计算出不考虑干涉因素的变电站厂界噪声;
c、根据a中获取的变电站主设备声源参数中主设备表面振动或近场声压,计算每台变电站设备的辐射声场;
d、根据c中每台变电站设备的辐射声场,利用区间随机模型摄动法,获取单独设备的相对相位差,及其对应的单独设备的干涉声纹;
e、采用步骤d中获取的干涉声纹,对步骤b中求取的不考虑干涉因素的变电站厂界噪声给予干涉修正,得到变电站厂界噪声;
所述步骤a中实测获取安全区域内设备表面振动或近场声压数据,对设备安全区域外进行振动参数反演方法获取设备声源参数,其中所述的安全区域振动、声压数据为变电站主设备2/3高度以下的安全测试范围内获取的变压器、电抗器表面振动数据、近场声压数据;
所述步骤a中对设备安全区域外进行振动参数反演方法获取设备声源参数具体是采用安全区域内的设备表面振动或近场声压,预测无法测量的安全区域外的设备表面振动或近场声压数据,具体为:
根据设备图纸或实测结构尺寸,采用有限元法计算复杂表面的模态向量Φn(x),获取模态导纳Yn(ω),n=1,2,...,N,在M点的振动位移PSD谱和相对相位通过现场设备的表面振动测量得到:
[W]M×1=[W(x1),W(x2),...,W(xM)]T
然后,求解模态力向量f(x):
变电站设备的模态力展开为:
模态位移展开为:
主变、高抗顶部套管位置,以及侧面安全范围以外的振动位移表示为:
模态位移的重建精度取决于特征矩阵与实测已知点振动位移的准确性,采用直接法求解模态位移:
采用结构有限元法计算上述模态力激发的变电站主设备振动分布,采用ABAQUS程序,将模态力考虑到变压器、电抗器设备的流固耦合模型中,计算外壳表面的振动速度;
所述步骤a中对设备安全区域外进行振动参数反演方法获取设备声源参数还包括电容器等效声源的参数,电容器等效声源参数的反演采用现场测量物理可观声压数据,以遗传算法,作为演化算法进行群体搜索,使描述可观测物理量的计算值与观测值差的目标函数最小,具体为:
电容器塔的每层电容器看成n2列×n1行源强相等、相位相同的等效点声源的组合,总声压p为:
式中:A为距离每个等效声源1米处的声压幅值;ri为受声点和声源的距离,单位为m;ω为圆角频率,单位为rad/s;k为声波的波数,单位为m-1,
各个等效点声源的相位表示为:
ri=r1+(i1-1)l1 sinθ+(i2-1)l2 cosθ
其中,1≤i1≤n1,1≤i2≤n2,
等效源的总声压可写为:
在使用遗传算法进行等效点源的相位获取时,构造目标函数为实测2D声辐射指向性和总声压计算得到的声指向性之差f(Φ),遗传算法的搜索空间[Ql,Qr]将在获取初始数据后给出,在该空间中随机产生Np个种群,并计算每个点声源的目标函数值,选择适应度最大的个体作为全局最优个体,即为电容器等效声源的参数。
2.根据权利要求1所述的一种变电站可听噪声预测的声波干涉修正方法,其特征在于,所述步骤a中的三维变电站仿真模型为DXF格式的地形文件和STP格式的CAD模型,这两种模型分别被用于步骤b中Soundplan软件的无干涉效应声场计算和步骤c中LMS Virtual lab软件的干涉声场计算;所述步骤b中采用的仿真软件是基于Soundplan的噪声仿真软件,所述步骤c中单独设备声场的计算采用LMS Virtual lab软件进行,在单独设备的声场计算时,已知边界取自步骤a中主设备的表面振动或近场声压数据;变电站数模取自步骤a中的CAD模型。
3.根据权利要求1所述的一种变电站可听噪声预测的声波干涉修正方法,其特征在于,所述步骤d中单独设备的相对相位差是指不同单独设备之间存在的相位差,同一个单独设备的相位关系是确定的,且已经在骤a中获取。
4.根据权利要求1所述的一种变电站可听噪声预测的声波干涉修正方法,其特征在于,所述步骤d中区间随机模型摄动法的具体实施步骤,
根据每台变电站设备的辐射声场,将等效声源的负荷波动、等效位置误差作为区间变量;
采用随机区间摄动法计算等效声源相位的期望与标准差变化范围;以现场测量的变电站主设备辐射声场观察点的实际相位数据作为参考解,验证评估随机区间摄动法的计算结果;
利用单独声源相位数据的参考解,调整每个声源独自运行时的声压分布,得到对应的单独设备的干涉声纹。
5.根据权利要求4所述的一种变电站可听噪声预测的声波干涉修正方法,其特征在于,所述步骤e中干涉修正的步骤是:
将相位调整后的每个设备的干涉声纹归一化;
将归一化后的单独设备的干涉声纹加权求和,得到考虑了干涉效应的声压分布模式;
把步骤b中得到的不考虑干涉因素的变电站厂界噪声换算成声波能量后,减去考虑了干涉效应的声压分布模式,得到综合修正后的变电站噪声分布和厂界敏感点噪声,完成变电站厂界噪声的修正。
6.一种变电站可听噪声预测的声波干涉修正***,其特征在于:用以实现如权利要求1-5任一所述的方法,***包括主设备声源参数获取模块(1),变电站设备、建筑物布局的空间几何尺寸获取模块(2),变电站三维仿真模型构建模块(3),不考虑干涉因素的变电站厂界噪声计算模块(4),辐射声场计算模块(5),干涉声纹获取模块(6),变电站厂界噪声生成模块(7),
所述主设备声源参数获取模块(1)采用实测安全区域振动、声压数据,对设备其他区域进行振动参数反演的方法,获取安全区域外的设备表面振动或近场声压数据;
所述变电站设备、建筑物布局的空间几何尺寸获取模块(2)用以获取变电站设备、建筑物布局的空间几何尺寸;
所述变电站三维仿真模型构建模块(3)根据变电站设备、建筑物布局的空间几何尺寸获取模块(2)获取的变电站设备、建筑物布局的空间几何尺寸构建变电站三维仿真模型;
所述不考虑干涉因素的变电站厂界噪声计算模块(4)以设备声源参数模型获取模块(1)的变电站主设备声功率为源强参数,设置对应的声环境参数及边界条件,仿真计算不考虑干涉因素的变电站厂界噪声;
所述辐射声场计算模块(5)以设备声源参数模型获取模块(1)的主设备表面振动或近场声压作为已知条件,单独计算每台变电站设备的辐射声场;
所述干涉声纹获取模块(6)利用区间随机模型摄动法,获取单独设备的相对相位差,及其对应的单独设备的干涉声纹;
所述变电站厂界噪声生成模块(7)采用干涉声纹获取模块(6)获取的干涉声纹,对不考虑干涉因素的变电站厂界噪声计算模块(4)中基于ISO9613-2下风向声压级计算法则求取的变电站噪声分布给予干涉修正,得到最终的变电站厂界噪声。
7.根据权利要求6所述的一种变电站可听噪声预测的声波干涉修正***,其特征在于:所述声波干涉修正***还包括安全区域振动、声压数据的测量装置。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110508682.7A CN113221413B (zh) | 2021-05-10 | 2021-05-10 | 一种变电站可听噪声预测的声波干涉修正方法及*** |
EP22732381.3A EP4116857A4 (en) | 2021-05-10 | 2022-05-09 | SOUND WAVE INTERFERENCE CORRECTION METHOD AND SYSTEM FOR PREDICTING AUDIBLE NOISE OF A TRANSFORMER SUBSTATION |
PCT/CN2022/091600 WO2022237708A1 (zh) | 2021-05-10 | 2022-05-09 | 一种变电站可听噪声预测的声波干涉修正方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110508682.7A CN113221413B (zh) | 2021-05-10 | 2021-05-10 | 一种变电站可听噪声预测的声波干涉修正方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113221413A CN113221413A (zh) | 2021-08-06 |
CN113221413B true CN113221413B (zh) | 2023-02-17 |
Family
ID=77094432
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110508682.7A Active CN113221413B (zh) | 2021-05-10 | 2021-05-10 | 一种变电站可听噪声预测的声波干涉修正方法及*** |
Country Status (3)
Country | Link |
---|---|
EP (1) | EP4116857A4 (zh) |
CN (1) | CN113221413B (zh) |
WO (1) | WO2022237708A1 (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113221413B (zh) * | 2021-05-10 | 2023-02-17 | 国网湖北省电力有限公司电力科学研究院 | 一种变电站可听噪声预测的声波干涉修正方法及*** |
CN114152877A (zh) * | 2021-11-26 | 2022-03-08 | 国网山东省电力公司电力科学研究院 | 一种同步调相机辅机噪声测试方法及*** |
CN115099115B (zh) * | 2022-08-24 | 2022-12-30 | 国网江西省电力有限公司电力科学研究院 | 一种户内变电站通风降噪深度强化学习设计方法 |
CN117192166A (zh) * | 2022-12-20 | 2023-12-08 | 国网河南省电力公司电力科学研究院 | 传声器支承设备、可听噪声测量装置 |
CN116108672B (zh) * | 2023-02-17 | 2024-01-23 | 南京声远声学科技有限公司 | 一种基于地理信息***的户外声传播预测模型构建方法 |
CN116400627A (zh) * | 2023-04-13 | 2023-07-07 | 深圳市丰源升科技有限公司 | 一种基于5g的工业远程控制***及方法 |
CN116992197B (zh) * | 2023-09-27 | 2024-01-23 | 国网北京市电力公司 | 一种变压器稳态运行噪声计算方法 |
CN117371255B (zh) * | 2023-12-06 | 2024-02-20 | 华中科技大学 | 一种旋转声源辐射噪声预测模型的构建方法及应用 |
CN117872273B (zh) * | 2024-03-11 | 2024-05-31 | 厦门市盛迅信息技术股份有限公司 | 基于人工智能的多环境声场声线识别方法和*** |
CN118052168B (zh) * | 2024-04-16 | 2024-07-26 | 山东职业学院 | 一种高速列车通过噪声的计算方法 |
CN118310627B (zh) * | 2024-06-07 | 2024-08-16 | 国网浙江省电力有限公司金华供电公司 | 一种城市变电站噪声***性控制方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105095609A (zh) * | 2015-09-21 | 2015-11-25 | 武汉大学 | 基于有限元法的变压器电磁振动噪声计算方法 |
CN112818553A (zh) * | 2021-02-07 | 2021-05-18 | 国网河北省电力有限公司沧州供电分公司 | 一种基于等效源的变压器辐射声场计算方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107860469B (zh) * | 2017-11-22 | 2019-09-10 | 重庆大学 | 一种基于正交多项式拟合的变电站噪声预测方法 |
CA2994312A1 (en) * | 2018-02-08 | 2019-08-08 | Electro-Kinetic Solutions Inc. | A control system and an adaptive electrokinetic dewatering system incorporating same for dewatering tailings |
CN110243464B (zh) * | 2019-07-19 | 2022-06-07 | 国家电网有限公司 | 一种适用于变电站现场的高抗噪声声功率级测试修正方法 |
CN111753405A (zh) * | 2020-06-05 | 2020-10-09 | 国网湖北省电力有限公司电力科学研究院 | 一种变电站声源设备源强求解方法及*** |
CN112685928B (zh) * | 2020-12-21 | 2024-01-30 | 中国电力科学研究院有限公司 | 一种基于三相电抗器声源模型的噪声预测方法及*** |
CN112347705B (zh) * | 2021-01-07 | 2021-04-02 | 中国电力科学研究院有限公司 | 一种建立变电站厂界噪声模型的方法和*** |
CN113221413B (zh) * | 2021-05-10 | 2023-02-17 | 国网湖北省电力有限公司电力科学研究院 | 一种变电站可听噪声预测的声波干涉修正方法及*** |
-
2021
- 2021-05-10 CN CN202110508682.7A patent/CN113221413B/zh active Active
-
2022
- 2022-05-09 WO PCT/CN2022/091600 patent/WO2022237708A1/zh unknown
- 2022-05-09 EP EP22732381.3A patent/EP4116857A4/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105095609A (zh) * | 2015-09-21 | 2015-11-25 | 武汉大学 | 基于有限元法的变压器电磁振动噪声计算方法 |
CN112818553A (zh) * | 2021-02-07 | 2021-05-18 | 国网河北省电力有限公司沧州供电分公司 | 一种基于等效源的变压器辐射声场计算方法 |
Non-Patent Citations (3)
Title |
---|
基于有限元法的变压器电磁振动噪声分析;胡静竹等;《电工技术学报》;20160810(第15期);全文 * |
特高压变电站噪声特性、预测及调控技术研究;宋凯;《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》;20170615;第1-5、11-27、46-60页 * |
特高压变电站电抗器复杂声场的研究与控制;蔡萱等;《噪声与振动控制》;20201031;第40卷(第5期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
WO2022237708A1 (zh) | 2022-11-17 |
CN113221413A (zh) | 2021-08-06 |
EP4116857A4 (en) | 2023-08-02 |
EP4116857A1 (en) | 2023-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113221413B (zh) | 一种变电站可听噪声预测的声波干涉修正方法及*** | |
CN109670257B (zh) | 一种换流站声场仿真方法及*** | |
CN106529055A (zh) | 一种基于应变模态振型相关性的模型修正方法 | |
CN105224741B (zh) | 电动汽车驱动***电磁辐射测试规划方法 | |
CN105260568B (zh) | 基于离散型卡尔曼滤波的超高层建筑风荷载反分析方法 | |
CN102967189A (zh) | ***冲击波超压时空场重建方法 | |
CN104239646A (zh) | 微组装组件振动疲劳寿命预测仿真模型验证方法和*** | |
CN105466661A (zh) | 基于改进卡尔曼滤波的超高层建筑风荷载反分析方法 | |
CN109670202A (zh) | 一种基于云模型的仿真可信度评估方法 | |
CN117709207B (zh) | 多波束测量船的测量布线设计方法、装置、设备和介质 | |
CN105004416B (zh) | 基于逆边界元法机械噪声远场声压预测方法 | |
Hou et al. | Estimation of virtual masses for structural damage identification | |
Du et al. | Modeling of multiply connected sound spaces by the surface coupling approach | |
CN117332562A (zh) | 一种基于环境动态仿真的变电站噪声强度计算方法和设备 | |
CN113804285B (zh) | 确定电力变压器组低频谐波噪声传播声场的方法及*** | |
CN114492220B (zh) | 基于OpenFOAM平台的海洋大气边界层流动特性的预测方法及*** | |
CN113312759A (zh) | 一种计算换流变压器组噪声源强的方法及*** | |
Meissner | Prediction of low-frequency sound field in rooms with complex-valued boundary conditions on walls | |
Gallovic et al. | High-frequency directivity in strong ground motion modeling methods | |
Froli et al. | Experimental static and dynamic tests on a large-scale free-form Voronoi grid shell mock-up in comparison with finite-element method results | |
CN113299027A (zh) | 一种隧道火灾监测方法及***、终端设备、存储介质 | |
CN113076600B (zh) | 一种直升机漂浮稳定性验证方法 | |
CN114611270A (zh) | 基于梯度下降算法的变电站内干涉声场修正方法及*** | |
CN111046542A (zh) | 一种评估二十面体的le网格计算特性方法及离散方法 | |
CN110108425A (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 |