CN111159928A - Transformer noise calculation method and system based on multi-line sound source model - Google Patents

Transformer noise calculation method and system based on multi-line sound source model Download PDF

Info

Publication number
CN111159928A
CN111159928A CN201911173488.7A CN201911173488A CN111159928A CN 111159928 A CN111159928 A CN 111159928A CN 201911173488 A CN201911173488 A CN 201911173488A CN 111159928 A CN111159928 A CN 111159928A
Authority
CN
China
Prior art keywords
transformer
model
predicted
sound source
noise
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.)
Granted
Application number
CN201911173488.7A
Other languages
Chinese (zh)
Other versions
CN111159928B (en
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 CN201911173488.7A priority Critical patent/CN111159928B/en
Publication of CN111159928A publication Critical patent/CN111159928A/en
Application granted granted Critical
Publication of CN111159928B publication Critical patent/CN111159928B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

The invention discloses a transformer noise calculation method and a system based on a multiline sound source model, wherein the method comprises the following steps: determining the relative position of a point to be predicted and a transformer; selecting one or more multi-line sound source equivalent models corresponding to the side faces of the transformers having influences on the points to be predicted according to the relative positions of the points to be predicted; inputting the positions of points to be predicted into the one or more multi-line sound source equivalent models, and obtaining the noise of one or more influencing transformer sides at the positions of the predicted points; superposing one or more noises at the position of the predicted point to obtain the comprehensive noise of the position to be predicted; the method and the system are divided into a plurality of small cubic modules by a finite element method, the small cubic modules are equivalent to a plurality of corresponding line sound sources, and a multi-line sound source model is established by a finite element-boundary element coupling method; the noise generated by the transformer can be accurately predicted for different prediction positions, and the problem of inaccurate near field noise prediction is solved.

Description

Transformer noise calculation method and system based on multi-line sound source model
Technical Field
The invention relates to the technical field of electric power, in particular to a transformer noise calculation method and system based on a multiline sound source model.
Background
When planning design and environment influence evaluation are carried out on an extra-high voltage alternating current power transmission and transformation project, the sound environment influence of a transformer substation needs to be predicted and calculated. Currently, sound environment of a transformer substation is often predicted and calculated by a design unit and an environmental impact evaluation unit by using soundPLAN or Cadna/A commercial noise prediction software, but the existing commercial noise prediction software has no calculation module specially aiming at transformer substation noise. The calculated result obtained by adopting the industrial noise module has a certain difference with the actual measured result, and particularly, the calculated value and the actual measured value of the area near the equipment such as a transformer, a reactor and the like have relatively large difference. The main reasons are that a certain difference exists between a sound source model established by noise prediction calculation and an actual operation equipment model, and the predicted source strength parameters are inaccurate, so that the error between a calculation result and an actual measurement result is large.
The method has the advantages that an accurate main equipment prediction model is established, the research on the ultra-high voltage transformer substation noise prediction technology is developed, the transformer substation noise is accurately and effectively predicted, the inevitable requirement for improving the rationality of the planning and design of the transformer substation is met, and the important means for ensuring that the environmental noise emission of the factory boundary reaches the standard after the engineering operation is also provided.
At present, the acoustic model of the transformer has been simplified into the following ones:
point sound source model: the transformer is equivalent to a point sound source, and the noise generated by the transformer is calculated based on the outdoor propagation theory of the point sound source. The method is simple, convenient and easy to calculate, but the use condition is very limited, and the method is accurate only when the far-field noise of the transformer is predicted, namely when the predicted point is far enough away from the transformer (generally larger than 3 times of the size of the transformer), the method can be equivalent to a point sound source. When the prediction point is closer to the transformer, the noise calculation error is larger, and especially in an extra-high voltage transformer substation, the main transformer equipment has larger volume and is not reasonable to be equivalent to a point sound source.
A surface sound source model: at the present stage, in an electric power design institute or a scientific research college, when the noise of a transformer substation is calculated, main equipment is equivalent to a plane sound source, or a cross section of the main equipment, or a vertical section of the main equipment, and the model is closer to the actual situation of a transformer and is simple and convenient to calculate, but the noise prediction is still inaccurate when a near-field area is calculated.
The body sound source model: the transformer oil tank is simplified into 6 surfaces of a cuboid, the radiation noise of the vibration of the top and the bottom of the transformer is ignored, the noise is considered to be all from the vibration of 4 tank walls, the calculation is based on a Helmholtz integral formula, the sound field and the sound pressure of the transformer are analyzed through the related calculation of the normal vibration acceleration of the surfaces of the 4 tank walls of the transformer, the method is closer to the actual situation of the transformer than the method when the method is simplified into a point sound source model, but the calculated sound field and the sound pressure of the transformer are still not very accurate.
A speaker array: the model is characterized in that on the basis that a transformer is regarded as a plane sound source, a series of loudspeakers are used for replacing each plane, and any sound field needing to be obtained is obtained by setting different loudspeaker amplitude values and initial phase angles, so that the sound field reconstruction is realized. The method is very complex, and the setting of the loudspeaker is difficult to be consistent with the original sound field of the transformer.
In conclusion, the prediction precision of the existing transformer acoustic model on noise cannot meet the prediction requirement under the extra-high voltage environment.
Disclosure of Invention
In order to solve the problem that the prediction precision of the existing transformer acoustic model on noise cannot meet the prediction requirement under the extra-high voltage environment in the background technology, the invention provides a transformer noise calculation method and a transformer noise calculation system based on a multiline sound source model, wherein the multiline sound source model is established by a finite element-boundary element coupling method through the method and the system, and the noise at the peripheral position of a transformer can be accurately predicted; the transformer noise calculation method based on the multiline sound source model comprises the following steps:
determining the relative position of a point to be predicted and a transformer;
selecting one or more multi-line sound source equivalent models corresponding to the side faces of the transformers which have influences on the points to be predicted according to the relative positions of the points to be predicted;
inputting the positions of the points to be predicted into the one or more multi-line sound source equivalent models, and obtaining noise of one or more influencing transformer sides at the positions of the predicted points;
and superposing the one or more noises at the position of the predicted point to obtain the comprehensive noise of the position of the point to be predicted.
Further, the selecting one or more multi-line sound source equivalent models corresponding to the side surfaces of the transformer having an influence on the point to be predicted includes:
the influence range of noise equivalently generated by each side face of the transformer is within 180 degrees of the vertical and outward direction of the side face;
determining whether each side surface has the influence of noise on the point to be predicted according to the relative position of the point to be predicted and the transformer; and selecting the equivalent model of the multi-line sound source corresponding to the side face which can generate noise influence.
Further, the method for establishing the equivalent model of the multi-line sound source corresponding to the side surface of each transformer comprises the following steps:
generating a geometric model in a COMSOL system according to the outline dimension of the transformer; dividing the geometric model into a plurality of cube small modules according to the structure of the geometric model;
determining material properties around the transformer, the material properties including air density and sound speed;
taking the geometric model as a domain of finite elements, and taking an infinite airspace outside the geometric model as a domain of boundary elements; defining the outer boundary of the set model as a boundary element-finite element coupling boundary;
calculating the corresponding sound power of the side surface of the transformer to be modeled;
inputting the divided geometric model, the material properties, the set finite element definition domain, the boundary element-finite element coupling boundary and the sound power into a COMSOL system, and establishing a multi-line sound source equivalent model corresponding to the side surface of the transformer to be modeled through the COMSOL system.
Further, the inputting the sound power level corresponding to the side surface of the transformer to be modeled includes:
receiving noise sound pressure measured at measuring points preset around the transformer; the preset measuring points are arranged at equivalent line sound source positions around the transformer; the equivalent line sound source of each small cube module is a line segment of a horizontal line on the small cube module, which is used for dividing the geometric model into a plurality of small cube modules, on the surface of the transformer;
calculating the sound power corresponding to each section of line sound source according to the acquired noise sound pressure and the sound power calculation method;
and summarizing the total sound power calculated by superposing the plurality of linear sound sources on the side surface of each transformer.
Further, the distance between the measuring points is fixed, and the number of the measuring points on each line sound source is determined according to the length of the line sound source;
the distance of the measuring point from the surface of the transformer is fixed.
Further, the calculation formula for calculating the acoustic power corresponding to each section of line sound source according to the collected noise sound pressure and acoustic power calculation method is as follows:
LP=LW-10lg(4πr2)
LW=10lg(W/W0)
wherein W is the acoustic power, W0Is a preset reference sound power; l isPIs the sound pressure level, LWAnd r is the distance from the sound pressure measuring point to the surface of the transformer.
The transformer noise calculation system based on the multiline sound source model comprises:
the model selection unit is used for determining the relative position of a point to be predicted and a transformer and selecting one or more multi-line sound source equivalent models corresponding to the side faces of the transformer, which have influences on the point to be predicted, according to the relative position of the point to be predicted;
the model calculation unit is used for inputting the positions of the points to be predicted to the one or more multi-line sound source equivalent models and obtaining the noise of one or more influencing transformer sides at the positions of the predicted points;
the model calculation unit is used for superposing the one or more noises at the position of the predicted point to obtain the comprehensive noise of the position of the point to be predicted.
Further, the model selecting unit is used for determining whether each side surface has the influence of noise on the point to be predicted according to the relative position of the point to be predicted and the transformer; selecting a multi-line sound source equivalent model corresponding to the side face which can generate noise influence;
each side of the transformer equivalently generates noise within an influence range of 180 degrees of a vertical and outward direction of the side.
Further, the system also comprises a model establishing unit;
the model establishing unit is used for generating a geometric model in a COMSOL system according to the overall dimension of the transformer; dividing the geometric model into a plurality of cube small modules according to the structure of the geometric model;
the model establishing unit is used for determining material properties around the transformer, wherein the material properties comprise air density and sound velocity;
the model establishing unit is used for setting the geometric model as a defined domain of a finite element and setting an infinite airspace outside the geometric model as a defined domain of a boundary element; setting the outer boundary of the set model as a boundary element-finite element coupling boundary;
the model establishing unit is used for calculating the acoustic power corresponding to the side surface of the transformer to be modeled;
the model establishing unit is used for inputting the divided geometric model, the material attribute, the set finite element definition domain, the boundary element-finite element coupling boundary and the sound power into a COMSOL system, and establishing a multi-line sound source equivalent model corresponding to the side face of the transformer to be modeled through the COMSOL system.
Further, the system also comprises a model establishing unit;
the model establishing unit is used for generating a geometric model in a COMSOL system according to the overall dimension of the transformer; dividing the geometric model into a plurality of cube small modules according to the structure of the geometric model;
the model establishing unit is used for determining material properties around the transformer, wherein the material properties comprise air density and sound velocity;
the model establishing unit is used for setting the geometric model as a defined domain of a finite element and setting an infinite airspace outside the geometric model as a defined domain of a boundary element; setting the outer boundary of the set model as a boundary element-finite element coupling boundary;
the model establishing unit is used for calculating the acoustic power corresponding to the side surface of the transformer to be modeled;
the model establishing unit is used for inputting the divided geometric model, the material attribute, the set finite element definition domain, the boundary element-finite element coupling boundary and the sound power into a COMSOL system, and establishing a multi-line sound source equivalent model corresponding to the side face of the transformer to be modeled through the COMSOL system.
Further, the distance between the measuring points is fixed, and the number of the measuring points on each line sound source is determined according to the length of the line sound source;
the distance of the measuring point from the surface of the transformer is fixed.
Further, the model building unit calculates the acoustic power corresponding to each section of line sound source according to the collected noise sound pressure and acoustic power calculation method by the following calculation formula:
LP=LW-10lg(4πr2)
LW=10lg(W/W0)
wherein W is the acoustic power, W0Is a preset reference sound power; l isPIs the sound pressure level, LWAnd r is the distance from the sound pressure measuring point to the surface of the transformer.
The invention has the beneficial effects that: the technical scheme of the invention provides a transformer noise calculation method and a system based on a multi-line sound source model, wherein the method and the system are divided into a plurality of small cubic modules by a finite element method, the small cubic modules are equivalent to a plurality of corresponding line sound sources, and the multi-line sound source model is established in a COMSOL system by a finite element-boundary element coupling method; according to the method and the system, an accurate sound source model is established, the noise generated by the transformer can be accurately predicted for different prediction positions, and the problem of inaccurate noise prediction of a near field region is solved.
Drawings
A more complete understanding of exemplary embodiments of the present invention may be had by reference to the following drawings in which:
FIG. 1 is a flow chart of a method for calculating transformer noise based on a multiline sound source model according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of the division of a small cube module and a line sound source on a geometric model of a transformer according to an embodiment of the present invention;
fig. 3 is a block diagram of a transformer noise calculation system based on a multiline sound source model according to an embodiment of the present invention.
Detailed Description
The exemplary embodiments of the present invention will now be described with reference to the accompanying drawings, however, the present invention may be embodied in many different forms and is not limited to the embodiments described herein, which are provided for complete and complete disclosure of the present invention and to fully convey the scope of the present invention to those skilled in the art. The terminology used in the exemplary embodiments illustrated in the accompanying drawings is not intended to be limiting of the invention. In the drawings, the same units/elements are denoted by the same reference numerals.
Unless otherwise defined, terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. Further, it will be understood that terms, such as those defined in commonly used dictionaries, should be interpreted as having a meaning that is consistent with their meaning in the context of the relevant art and will not be interpreted in an idealized or overly formal sense.
FIG. 1 is a flow chart of a method for calculating transformer noise based on a multiline sound source model according to an embodiment of the present invention; as shown in fig. 1, the method includes:
step 110, determining the relative position of a point to be predicted and a transformer;
in order to predict the noise influence on the environment after the transformer is built, simulating the relative position of a point to be predicted and the transformer; in the embodiment, a space coordinate system is established by taking the transformer as an origin; and transforming the relative position of the point to be predicted and the transformer into the coordinate of the space coordinate system so as to realize the clear expression of the relative position.
Step 120, selecting one or more multi-line sound source equivalent models corresponding to the side faces of the transformers which have influences on the points to be predicted according to the relative positions of the points to be predicted;
the influence range of noise equivalently generated by each side face of the transformer is within 180 degrees of the vertical and outward direction of the side face; that is, a spatial region covered in a direction perpendicular to and outward from a plane in which the side surface is located belongs to a range in which the side surface may be affected by noise.
Determining whether each side surface has the influence of noise on the point to be predicted according to the relative position of the point to be predicted and the transformer; and selecting the equivalent model of the multi-line sound source corresponding to the side face which can generate noise influence.
In the coordinate system established by taking the transformer as the origin in the embodiment, in order to improve the accuracy of the calculation of the near field region, the volume of the transformer cannot be ignored, that is, the geometric center point of the transformer is taken as the origin or a certain bottom angle is taken as the origin; each surface of the transformer can be displayed in a coordinate system, and the specific area covered by one or more planes of the point to be measured is determined by the noise influence determination method;
for example, a point to be predicted is right in front of a side surface of a transformer at the same level (i.e., the point to be predicted is taken as a perpendicular line to the side surface, and the perpendicular point is in the side surface), and at this time, the point to be predicted is only affected by noise generated by the side surface.
Step 130, inputting the position of the point to be predicted to the one or more multi-line sound source equivalent models, and obtaining noise of one or more influencing transformer sides at the position of the predicted point;
further, before step 130, a multiline sound source equivalent model of each side of the transformer needs to be established; the method comprises the following specific steps:
step 131, generating a geometric model in a COMSOL system according to the overall dimension of the transformer; dividing the geometric model into a plurality of cube small modules according to the structure of the geometric model;
as shown in fig. 2, fig. 2 is a schematic diagram of a small cube module division and a line sound source on a geometric model of a transformer according to an embodiment of the present invention; in fig. 2, a small cube module division method is illustrated: the length of the simulated transformer geometric model is 9m, the width and the length of one side are 4m, the width and the length of the other side are 4.5m, and the height of the other side is 4.5m, at the moment, the structure of the transformer model consists of 6 cuboids with the length of 3m, the width of 2.5m, the height of 1.5m, 3 cuboids with the length of 2.5m, the width of 2.5m, the height of 1.5m, 3 cuboids with the length of 3m, the width of 2m and the height of 1.5m, and the total number of the cuboids is 15.
The equivalent line sound source of each small cube module is a line segment of a horizontal line on the small cube module, which is used for dividing the geometric model into a plurality of small cube modules, on the surface of the transformer; namely, the upper and lower outer boundaries of the middle layer structure block are used as line sound sources (thickened lines in the figure), wherein the structure blocks are rectangles with different sizes, so that the length and the position of the obtained line sound sources are different.
And for the line sound sources on the same side, all the line sound sources on the changed side are obtained.
Step 132, determining material properties around the transformer, wherein the material properties comprise air density and sound velocity;
step 133, using the geometric model as a domain of finite elements, and using an infinite airspace outside the geometric model as a domain of boundary elements; defining the outer boundary of the set model as a boundary element-finite element coupling boundary;
in order to observe the sound pressure and the sound pressure level change conveniently, the definition domain of the boundary element is defined as a fixed region outside the set model;
step 134, calculating the acoustic power corresponding to the side surface of the transformer to be modeled;
the specific steps of calculating the acoustic power are as follows:
step 1341, receiving noise sound pressure measured at measurement points preset around the transformer; the preset measuring points are arranged at equivalent line sound source positions around the transformer;
the distance between the measuring points is fixed, and the number of the measuring points on each line sound source is determined according to the length of the line sound source; the distance between the measuring points may also be not fixed, but the line sound sources with different lengths and the corresponding number of the measuring points should be matched.
The distance of the measuring point from the surface of the transformer is fixed.
Step 1342, calculating the corresponding sound power of each section of line sound source according to the collected noise sound pressure and the sound power calculation method;
the calculation formula is as follows:
LP=LW-10lg(4πr2)
LW=10lg(W/W0)
wherein W is the acoustic power, W0Is a preset reference sound power; l isPIs the sound pressure level, LWAnd r is the distance from the sound pressure measuring point to the surface of the transformer.
Step 1343, summarizing the total sound power calculated by overlapping the plurality of line sound sources on the side surface of each transformer.
And 135, inputting the divided geometric model, the material attribute, the set finite element definition domain, the boundary element-finite element coupling boundary and the sound power into a COMSOL system, and establishing a multi-line sound source equivalent model corresponding to the side surface of the transformer to be modeled through the COMSOL system.
And respectively defining the determined contents in a COMSOL system to generate a simulated multiline sound source equivalent model.
And 140, superposing the one or more noises at the position of the point to be predicted to obtain the comprehensive noise of the position of the point to be predicted.
In this embodiment, the noise generated by any one side surface may be represented in a space coordinate system in a form of a vector, and when a plurality of noises are superimposed at the predicted point position, a vector sum is obtained for a plurality of vectors corresponding to the plurality of noises, that is, the sum is affected by the comprehensive noise at the predicted point position.
FIG. 3 is a block diagram of a transformer noise calculation system based on a multiline sound source model according to an embodiment of the present invention; as shown in fig. 3, the system includes:
the model selecting unit 310 is configured to determine a relative position of a point to be predicted and a transformer, and select one or more multi-line sound source equivalent models corresponding to transformer sides having influences on the point to be predicted according to the relative position of the point to be predicted;
the model selecting unit 310 is configured to determine whether each side surface has an influence on noise generated by a point to be predicted according to a relative position of the point to be predicted and the transformer; selecting a multi-line sound source equivalent model corresponding to the side face which can generate noise influence;
each side of the transformer equivalently generates noise within an influence range of 180 degrees of a vertical and outward direction of the side.
A model calculating unit 320, wherein the model calculating unit 320 is configured to input the position of the point to be predicted to the one or more multiline sound source equivalent models, and obtain noise of one or more influencing transformer sides at the position of the predicted point;
the model calculating unit 320 is configured to superimpose the one or more noises at a to-be-predicted point to obtain a comprehensive noise at the to-be-predicted point.
The system further comprises a model building unit 330;
the model establishing unit 330 is configured to generate a geometric model in the COMSOL system according to the external dimensions of the transformer; dividing the geometric model into a plurality of cube small modules according to the structure of the geometric model;
the model building unit 330 is configured to determine material properties around the transformer, where the material properties include air density and sound velocity;
the model establishing unit 330 is configured to set the geometric model as a domain of finite elements, and set an infinite airspace outside the geometric model as a domain of boundary elements; setting the outer boundary of the set model as a boundary element-finite element coupling boundary;
the model establishing unit 330 is configured to calculate an acoustic power corresponding to a side surface of the transformer to be modeled;
the model establishing unit 330 is configured to receive noise sound pressure measured at preset measurement points around the transformer; the preset measuring points are arranged at equivalent line sound source positions around the transformer; the equivalent line sound source of each small cube module is a line segment of a horizontal line on the small cube module, which is used for dividing the geometric model into a plurality of small cube modules, on the surface of the transformer;
the distance between the measuring points is fixed, and the number of the measuring points on each line sound source is determined according to the length of the line sound source;
the distance of the measuring point from the surface of the transformer is fixed.
The model establishing unit 330 is configured to calculate the acoustic power corresponding to each section of line sound source according to the collected noise sound pressure and the acoustic power calculation method, and sum up the total acoustic power calculated by superimposing the plurality of line sound sources on the side surface of each transformer.
The model establishing unit 330 calculates the acoustic power corresponding to each section of line sound source according to the collected noise sound pressure and acoustic power calculation method by the following calculation formula:
LP=LW-10lg(4πr2)
LW=10lg(W/W0)
wherein W is the acoustic power, W0Is a preset reference sound power; l isPIs the sound pressure level, LWAnd r is the distance from the sound pressure measuring point to the surface of the transformer.
The model establishing unit 330 is configured to input the divided geometric model, the material property, the set finite element definition domain, the boundary element-finite element coupling boundary, and the acoustic power into a COMSOL system, and establish a multi-line acoustic source equivalent model corresponding to a side surface of the transformer to be modeled by the COMSOL system.
In the description provided herein, numerous specific details are set forth. However, it is understood that embodiments of the disclosure may be practiced without these specific details. In some instances, well-known methods, structures and techniques have not been shown in detail in order not to obscure an understanding of this description.
Those skilled in the art will appreciate that the modules in the device in an embodiment may be adaptively changed and disposed in one or more devices different from the embodiment. The modules or units or components of the embodiments may be combined into one module or unit or component, and furthermore they may be divided into a plurality of sub-modules or sub-units or sub-components. All of the features disclosed in this specification (including any accompanying claims, abstract and drawings), and all of the processes or elements of any method or apparatus so disclosed, may be combined in any combination, except combinations where at least some of such features and/or processes or elements are mutually exclusive. Each feature disclosed in this specification (including any accompanying claims, abstract and drawings) may be replaced by alternative features serving the same, equivalent or similar purpose, unless expressly stated otherwise. Reference to step numbers in this specification is only for distinguishing between steps and is not intended to limit the temporal or logical relationship between steps, which includes all possible scenarios unless the context clearly dictates otherwise.
Moreover, those skilled in the art will appreciate that while some embodiments described herein include some features included in other embodiments, rather than other features, combinations of features of different embodiments are meant to be within the scope of the disclosure and form different embodiments. For example, any of the embodiments claimed in the claims can be used in any combination.
Various component embodiments of the disclosure may be implemented in hardware, or in software modules running on one or more processors, or in a combination thereof. The present disclosure may also be embodied as device or system programs (e.g., computer programs and computer program products) for performing a portion or all of the methods described herein. Such programs implementing the present disclosure may be stored on a computer-readable medium or may be in the form of one or more signals. Such a signal may be downloaded from an internet website or provided on a carrier signal or in any other form.
It should be noted that the above-mentioned embodiments illustrate rather than limit the disclosure, and that those skilled in the art will be able to design alternative embodiments without departing from the scope of the appended claims. The word "comprising" does not exclude the presence of elements or steps not listed in a claim. The word "a" or "an" preceding an element does not exclude the presence of a plurality of such elements. The disclosure may be implemented by means of hardware comprising several distinct elements, and by means of a suitably programmed computer. In the unit claims enumerating several systems, several of these systems may be embodied by one and the same item of hardware.
The foregoing is directed to embodiments of the present disclosure, and it is noted that numerous improvements, modifications, and variations may be made by those skilled in the art without departing from the spirit of the disclosure, and that such improvements, modifications, and variations are considered to be within the scope of the present disclosure.

Claims (12)

1. A method of transformer noise calculation based on a multiline sound source model, the method comprising:
determining the relative position of a point to be predicted and a transformer;
selecting one or more multi-line sound source equivalent models corresponding to the side faces of the transformers which have influences on the points to be predicted according to the relative positions of the points to be predicted;
inputting the positions of the points to be predicted into the one or more multi-line sound source equivalent models, and obtaining noise of one or more influencing transformer sides at the positions of the predicted points;
and superposing the one or more noises at the position of the predicted point to obtain the comprehensive noise of the position of the point to be predicted.
2. The method according to claim 1, wherein said selecting one or more equivalent models of the multi-line sound source corresponding to the side of the transformer having an influence on the point to be predicted comprises:
the influence range of noise equivalently generated by each side face of the transformer is within 180 degrees of the vertical and outward direction of the side face;
determining whether each side surface has the influence of noise on the point to be predicted according to the relative position of the point to be predicted and the transformer; and selecting the equivalent model of the multi-line sound source corresponding to the side face which can generate noise influence.
3. The method according to claim 1, wherein the method for establishing the equivalent model of the multi-line sound source corresponding to each transformer side comprises:
generating a geometric model in a COMSOL system according to the outline dimension of the transformer; dividing the geometric model into a plurality of cube small modules according to the structure of the geometric model;
determining material properties around the transformer, the material properties including air density and sound speed;
taking the geometric model as a domain of finite elements, and taking an infinite airspace outside the geometric model as a domain of boundary elements; defining the outer boundary of the set model as a boundary element-finite element coupling boundary;
calculating the corresponding sound power of the side surface of the transformer to be modeled;
inputting the divided geometric model, the material properties, the set finite element definition domain, the boundary element-finite element coupling boundary and the sound power into a COMSOL system, and establishing a multi-line sound source equivalent model corresponding to the side surface of the transformer to be modeled through the COMSOL system.
4. The method of claim 3, wherein inputting the corresponding acoustic power level of the transformer side to be modeled comprises:
receiving noise sound pressure measured at measuring points preset around the transformer; the preset measuring points are arranged at equivalent line sound source positions around the transformer; the equivalent line sound source of each small cube module is a line segment of a horizontal line on the small cube module, which is used for dividing the geometric model into a plurality of small cube modules, on the surface of the transformer;
calculating the sound power corresponding to each section of line sound source according to the acquired noise sound pressure and the sound power calculation method;
and summarizing the total sound power calculated by superposing the plurality of linear sound sources on the side surface of each transformer.
5. The method of claim 4, wherein:
the distance between the measuring points is fixed, and the number of the measuring points on each line sound source is determined according to the length of the line sound source;
the distance of the measuring point from the surface of the transformer is fixed.
6. The method according to claim 4, wherein the calculation formula for calculating the sound power corresponding to each line sound source according to the collected noise sound pressure and sound power calculation method is as follows:
LP=LW-10lg(4πr2)
LW=10lg(W/W0)
wherein W is the acoustic power, W0Is a preset reference sound power; l isPIs the sound pressure level, LWAnd r is the distance from the sound pressure measuring point to the surface of the transformer.
7. A transformer noise calculation system based on a multiline sound source model, the system comprising:
the model selection unit is used for determining the relative position of a point to be predicted and a transformer and selecting one or more multi-line sound source equivalent models corresponding to the side faces of the transformer, which have influences on the point to be predicted, according to the relative position of the point to be predicted;
the model calculation unit is used for inputting the positions of the points to be predicted to the one or more multi-line sound source equivalent models and obtaining the noise of one or more influencing transformer sides at the positions of the predicted points;
the model calculation unit is used for superposing the one or more noises at the position of the point to be predicted to obtain the comprehensive noise of the position of the point to be predicted.
8. The system of claim 7, wherein:
the model selecting unit is used for determining whether each side surface generates noise influence on the point to be predicted according to the relative position of the point to be predicted and the transformer; selecting a multi-line sound source equivalent model corresponding to the side face which can generate noise influence;
each side of the transformer equivalently generates noise within an influence range of 180 degrees of a vertical and outward direction of the side.
9. The system of claim 1, wherein: the system further comprises a model building unit;
the model establishing unit is used for generating a geometric model in a COMSOL system according to the overall dimension of the transformer; dividing the geometric model into a plurality of cube small modules according to the structure of the geometric model;
the model establishing unit is used for determining material properties around the transformer, wherein the material properties comprise air density and sound velocity;
the model establishing unit is used for setting the geometric model as a defined domain of a finite element and setting an infinite airspace outside the geometric model as a defined domain of a boundary element; setting the outer boundary of the set model as a boundary element-finite element coupling boundary;
the model establishing unit is used for calculating the acoustic power corresponding to the side surface of the transformer to be modeled;
the model establishing unit is used for inputting the divided geometric model, the material attribute, the set finite element definition domain, the boundary element-finite element coupling boundary and the sound power into a COMSOL system, and establishing a multi-line sound source equivalent model corresponding to the side face of the transformer to be modeled through the COMSOL system.
10. The system of claim 9, wherein:
the model establishing unit is used for receiving noise sound pressure measured at preset measuring points around the transformer; the preset measuring points are arranged at equivalent line sound source positions around the transformer; the equivalent line sound source of each small cube module is a line segment of a horizontal line on the small cube module, which is used for dividing the geometric model into a plurality of small cube modules, on the surface of the transformer;
the model establishing unit is used for calculating the corresponding sound power of each section of linear sound source according to the collected noise sound pressure and the sound power calculating method, and summarizing the total sound power calculated by overlapping a plurality of linear sound sources on the side surface of each transformer.
11. The system of claim 10, wherein:
the distance between the measuring points is fixed, and the number of the measuring points on each line sound source is determined according to the length of the line sound source;
the distance of the measuring point from the surface of the transformer is fixed.
12. The system of claim 10, wherein: the model establishing unit calculates the acoustic power corresponding to each section of line sound source according to the collected noise sound pressure and acoustic power calculation method by the following calculation formula:
LP=LW-10lg(4πr2)
LW=10lg(W/W0)
wherein W is the acoustic power, W0Is a preset reference sound power; l isPIs the sound pressure level, LWAnd r is the distance from the sound pressure measuring point to the surface of the transformer.
CN201911173488.7A 2019-11-26 2019-11-26 Transformer noise calculation method and system based on multi-line sound source model Active CN111159928B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911173488.7A CN111159928B (en) 2019-11-26 2019-11-26 Transformer noise calculation method and system based on multi-line sound source model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911173488.7A CN111159928B (en) 2019-11-26 2019-11-26 Transformer noise calculation method and system based on multi-line sound source model

Publications (2)

Publication Number Publication Date
CN111159928A true CN111159928A (en) 2020-05-15
CN111159928B CN111159928B (en) 2024-07-02

Family

ID=70556076

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911173488.7A Active CN111159928B (en) 2019-11-26 2019-11-26 Transformer noise calculation method and system based on multi-line sound source model

Country Status (1)

Country Link
CN (1) CN111159928B (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111753405A (en) * 2020-06-05 2020-10-09 国网湖北省电力有限公司电力科学研究院 Transformer substation sound source equipment source intensity solving method and system
CN111985113A (en) * 2020-08-28 2020-11-24 国网青海省电力公司电力科学研究院 Method and device for predicting noise of power transmission line
CN112380668A (en) * 2020-10-10 2021-02-19 中国电力科学研究院有限公司 Method and device for noise test based on transformer reduced scale model
CN113720441A (en) * 2021-08-26 2021-11-30 国网上海市电力公司 Transformer substation perimeter noise distribution reproduction method
CN113804285A (en) * 2021-05-13 2021-12-17 中国电力科学研究院有限公司 Method and system for determining low-frequency harmonic noise propagation sound field of power transformer bank
CN114048644A (en) * 2021-10-21 2022-02-15 国网新源控股有限公司 Method for constructing acoustic diagnosis database of dry-type transformer
CN114117629A (en) * 2021-11-08 2022-03-01 北京理工大学 Equivalent sound source simulation method of loudspeaker
CN114152824A (en) * 2021-11-15 2022-03-08 国网重庆市电力公司电力科学研究院 Real-time monitoring and optimal control method for environmental noise of transformer substation
WO2023103468A1 (en) * 2021-12-10 2023-06-15 国网河北省电力有限公司经济技术研究院 Method for determining transformer multi-sound-source noise equivalent model, terminal, and storage medium
CN118296914A (en) * 2024-05-29 2024-07-05 邢台职业技术学院 Camera factory building noise prediction method based on FEM-BEM dotted line composite sound source model

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11102154A (en) * 1997-02-27 1999-04-13 Sekisui Chem Co Ltd Noise simulation device and method
CN102722607A (en) * 2012-05-25 2012-10-10 山东大学 Method for optimizing parameters of secondary sound source in active noise reduction of power transformer
CN104239742A (en) * 2014-09-29 2014-12-24 国家电网公司 Transformer far-field noise prediction method and system
CN106096184A (en) * 2016-06-24 2016-11-09 中国电力科学研究院 A kind of noise calculation method and system of transformator multiple spot Source Model based on sound intensity technique

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11102154A (en) * 1997-02-27 1999-04-13 Sekisui Chem Co Ltd Noise simulation device and method
CN102722607A (en) * 2012-05-25 2012-10-10 山东大学 Method for optimizing parameters of secondary sound source in active noise reduction of power transformer
CN104239742A (en) * 2014-09-29 2014-12-24 国家电网公司 Transformer far-field noise prediction method and system
CN106096184A (en) * 2016-06-24 2016-11-09 中国电力科学研究院 A kind of noise calculation method and system of transformator multiple spot Source Model based on sound intensity technique

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
HU JINGZHU ET AL: "Analysis of substation noise prediction based on equivalent source method", 2014 IEEE PES ASIA-PACIFIC POWER AND ENERGY ENGINEERING CONFERENCE (APPEEC). PROCEEDINGS, 30 March 2015 (2015-03-30), pages 1 - 5 *
孙涛等: "特高压变压器噪声源模型及仿真分析", 高电压技术, vol. 40, no. 9, 30 September 2014 (2014-09-30), pages 2750 - 2756 *
章彬等: "电力变压器低频噪声的模拟与等效", 武汉大学学报, vol. 47, no. 5, 31 October 2014 (2014-10-31), pages 642 - 647 *
胡静竹等: "超高压变电站噪声特性分析与预测模型研究", 高压电器, vol. 53, no. 2, 16 February 2017 (2017-02-16), pages 32 - 38 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111753405A (en) * 2020-06-05 2020-10-09 国网湖北省电力有限公司电力科学研究院 Transformer substation sound source equipment source intensity solving method and system
CN111985113A (en) * 2020-08-28 2020-11-24 国网青海省电力公司电力科学研究院 Method and device for predicting noise of power transmission line
CN111985113B (en) * 2020-08-28 2024-04-16 国网青海省电力公司电力科学研究院 Method and device for predicting noise of power transmission line
CN112380668A (en) * 2020-10-10 2021-02-19 中国电力科学研究院有限公司 Method and device for noise test based on transformer reduced scale model
CN113804285A (en) * 2021-05-13 2021-12-17 中国电力科学研究院有限公司 Method and system for determining low-frequency harmonic noise propagation sound field of power transformer bank
CN113804285B (en) * 2021-05-13 2024-04-05 中国电力科学研究院有限公司 Method and system for determining low-frequency harmonic noise propagation sound field of power transformer group
CN113720441B (en) * 2021-08-26 2023-09-12 国网上海市电力公司 Substation perimeter noise distribution reproduction method
CN113720441A (en) * 2021-08-26 2021-11-30 国网上海市电力公司 Transformer substation perimeter noise distribution reproduction method
CN114048644A (en) * 2021-10-21 2022-02-15 国网新源控股有限公司 Method for constructing acoustic diagnosis database of dry-type transformer
CN114048644B (en) * 2021-10-21 2024-04-26 国网新源控股有限公司 Method for constructing acoustic diagnosis database of dry-type transformer
CN114117629A (en) * 2021-11-08 2022-03-01 北京理工大学 Equivalent sound source simulation method of loudspeaker
CN114152824A (en) * 2021-11-15 2022-03-08 国网重庆市电力公司电力科学研究院 Real-time monitoring and optimal control method for environmental noise of transformer substation
WO2023103468A1 (en) * 2021-12-10 2023-06-15 国网河北省电力有限公司经济技术研究院 Method for determining transformer multi-sound-source noise equivalent model, terminal, and storage medium
CN118296914A (en) * 2024-05-29 2024-07-05 邢台职业技术学院 Camera factory building noise prediction method based on FEM-BEM dotted line composite sound source model

Also Published As

Publication number Publication date
CN111159928B (en) 2024-07-02

Similar Documents

Publication Publication Date Title
CN111159928B (en) Transformer noise calculation method and system based on multi-line sound source model
CN106096184B (en) Noise calculation method and system of transformer multipoint sound source model based on sound intensity method
CN112685928B (en) Noise prediction method and system based on three-phase reactor sound source model
US9658272B2 (en) Method and device for determining the defect type of a partial discharge
Lakshmynarayanana et al. Application of a two-way partitioned method for predicting the wave-induced loads of a flexible containership
JP2017003574A (en) System and method for managing corrosion
US9117041B2 (en) Magnetic property analyzing apparatus and method
JP2016020841A (en) Noises and/or vibrations monitoring method and monitoring system
JP6027121B2 (en) Measuring position indicating device, measuring position indicating method
Shi et al. Damage detection in shear buildings using different estimated curvature
CN115758202A (en) Transformer noise signal identification method based on multi-physical field coupling simulation
US10268785B2 (en) Noise detection device
US20120101748A1 (en) Recording medium for acoustic analysis program, acoustic analysis method, and acoustic analysis device
Godinho et al. 3D numerical modelling of acoustic horns using the method of fundamental solutions
US10375501B2 (en) Method and device for quickly determining location-dependent pulse responses in signal transmission from or into a spatial volume
Siano et al. Sensitivity analysis and correlation Experimental/Numerical FEM-BEM for Noise Reduction assessment of an engine beauty cover
CN102270249B (en) Method for identifying characteristic frequency of parts
US20220187483A1 (en) Earthquake estimation method, non-transitory computer readable medium, and earthquake estimation device
JP2015215698A (en) Electromagnetic field simulation program
JP2020085631A (en) Internal state estimation system
JP5390447B2 (en) Method for estimating noise countermeasure effect
CN108154557A (en) A kind of tri patch method in the mute region based on domestic environment
CN110569526B (en) Method for analyzing sound radiation characteristics of power transformer in multiple operating states
JP2018197965A (en) Design support device of small speaker and design support method of speaker
CN106599395A (en) Numerical simulation calculation method for noise of oil immersed transformer

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