Transformer substation noise coupling simulation analysis method and device based on sound ray simulation scene

Abstract

The invention provides a transformer substation noise coupling simulation analysis method and device based on a sound ray simulation scene, and the method comprises the steps: constructing a three-dimensional space scene based on the spatial position and geometric dimension information of a spatial entity in a transformer substation; controlling the noise source to emit a plurality of sound rays according to a preset angle, and determining a collision position of each sound ray with a space entity when the sound ray is propagated in the three-dimensional space scene and a reflection direction after each collision; on the basis of the propagation path, the collision position and the reflection direction of each sound ray, acoustic parameter changes generated by distance attenuation and collision reflection in the propagation process of each sound ray are determined, and an acoustic parameter change sequence is generated; coupling is carried out based on the target acoustic parameter change sequence of the propagation path passing through the spatial points, and a noise coupling result is determined; and traversing the noise coupling result of each spatial point to obtain the noise coupling distribution condition of the space where the transformer substation is located. According to the method, the dynamic propagation and coupling characteristics of the noise in a complex space scene are accurately reflected.

Description

Technical Field

[0001] This invention relates to the field of computer technology, and in particular to a method and apparatus for noise coupling simulation analysis of substations based on sound ray simulation scenarios. Background Technology

[0002] As a critical hub in the power system, substations generate noise during operation that not only disturbs the living environment of nearby residents but may also affect the stable operation of equipment. Therefore, accurate analysis of substation noise is of great significance for developing effective noise reduction measures and optimizing substation layout.

[0003] Currently, substation noise analysis primarily employs acoustic parameter analysis methods based on fixed monitoring points. This involves setting up several fixed monitoring points within and around the substation, using noise sensors to collect parameters such as sound pressure level and frequency at different times, and then combining these with empirical formulas or statistical models to assess the propagation characteristics and impact range of the noise. However, substations contain numerous spatial entities (equipment, walls, etc.), and noise propagation is subject to reflection and refraction, making it impossible for fixed monitoring points to capture attenuation and changes along the propagation path. Furthermore, noise from different sources couples and superimposes in space, and fixed monitoring points can only obtain local superposition results, failing to comprehensively reflect the overall coupling patterns. Moreover, when equipment operating conditions change dynamically (such as the start-up and shutdown of some equipment), the noise distribution changes in real time, making it difficult for fixed monitoring points to respond quickly to such dynamic changes. Therefore, existing methods struggle to accurately reflect the dynamic propagation and coupling characteristics of noise in complex spatial scenarios. Summary of the Invention

[0004] This invention provides a method and apparatus for substation noise coupling simulation analysis based on acoustic ray simulation scenarios, aiming to accurately reflect the dynamic propagation and coupling characteristics of noise in complex spatial scenarios.

[0005] In a first aspect, the present invention provides a substation noise coupling simulation analysis method based on acoustic ray simulation scenarios, including: A three-dimensional spatial scene is constructed based on the spatial location and geometric dimensions of spatial entities within the substation. For each noise source in the substation, multiple sound rays are emitted by the noise source according to a preset angle, and the collision position of each sound ray with the spatial entity when it propagates in the three-dimensional space scene and the reflection direction after each collision are determined until the sound ray propagates to the boundary of the three-dimensional space scene; each sound ray carries the initial acoustic parameters of the noise source. Based on the propagation path, collision location, and reflection direction of each sound ray, the changes in acoustic parameters of each sound ray during propagation due to distance attenuation and collision reflection are determined, generating an acoustic parameter change sequence for each sound ray; the acoustic parameter change sequence includes the acoustic parameter values ​​and corresponding spatial coordinates of each node on the propagation path of the sound ray. For each spatial point in the three-dimensional spatial scene, the noise coupling result of each spatial point is determined by coupling based on the sequence of changes in the target acoustic parameters of the target sound ray passing through the spatial point through the propagation path. By traversing the noise coupling results of each spatial point in the three-dimensional spatial scene, the noise coupling distribution of the space where the substation is located is obtained.

[0006] In a second aspect, the present invention also provides a substation noise coupling simulation analysis device based on a sound ray simulation scenario, applied to the substation noise coupling simulation analysis method based on a sound ray simulation scenario as described in the first aspect; the substation noise coupling simulation analysis device based on a sound ray simulation scenario includes: The spatial scene construction module is used to construct a three-dimensional spatial scene based on the spatial location and geometric dimensions of spatial entities within the substation. The sound ray path tracking module is used to control each noise source in the substation to emit multiple sound rays at a preset angle, and to determine the collision position of each sound ray with the spatial entity when it propagates in the three-dimensional space scene and the reflection direction after each collision, until the sound ray propagates to the boundary of the three-dimensional space scene; each sound ray carries the initial acoustic parameters of the noise source. The acoustic parameter monitoring module is used to determine the changes in acoustic parameters of each sound ray during propagation due to distance attenuation and collision reflection, based on the propagation path, collision location, and reflection direction of each sound ray, and to generate an acoustic parameter change sequence for each sound ray; the acoustic parameter change sequence includes the acoustic parameter values ​​and corresponding spatial coordinates of each node on the propagation path of the sound ray; A single spatial point noise coupling module is used to couple each spatial point in the three-dimensional spatial scene based on the sequence of changes in target acoustic parameters of the target sound ray passing through the spatial point through the propagation path, and to determine the noise coupling result of each spatial point. The full-space point noise coupling module is used to traverse the noise coupling results of each spatial point in the three-dimensional spatial scene to obtain the noise coupling distribution of the space where the substation is located.

[0007] Thirdly, the present invention also provides an electronic device, comprising: a memory for storing computer software programs; and a processor for reading and executing the computer software programs, thereby realizing the substation noise coupling simulation analysis method based on sound ray simulation scenarios as described above.

[0008] Fourthly, the present invention also provides a non-transitory computer-readable storage medium storing a computer software program, which, when executed by a processor, implements the substation noise coupling simulation analysis method based on a sound ray simulation scenario as described above.

[0009] Fifthly, the present invention also provides a computer program product, including a computer program that, when executed by a processor, implements the substation noise coupling simulation analysis method based on a sound ray simulation scenario as described above.

[0010] The substation noise coupling simulation analysis method based on acoustic ray simulation scenarios provided in this invention emits multiple acoustic rays and tracks their propagation paths. This allows for a complete record of noise reflection and attenuation during propagation, covering all possible propagation paths within the substation and eliminating limitations imposed by fixed points. Furthermore, by tracking the acoustic parameter changes of each acoustic ray carrying the initial acoustic parameters of the corresponding noise source at various nodes, the method determines the noise coupling result at each spatial point by coupling based on the acoustic parameter change sequences of all propagation paths passing through that spatial point. This accurately reflects the coupling results of different noise sources at that spatial point. Moreover, when the equipment's operating state changes dynamically, the acoustic parameters of the noise source can be updated, acoustic rays can be re-emitted, and noise coupling can be performed to quickly obtain a new noise coupling distribution. This allows for dynamic response to changes in equipment state. Therefore, this invention, through full-path acoustic ray tracking, acoustic parameter change recording, and full-spatial-point coupling, can accurately reflect the dynamic propagation and coupling characteristics of noise in complex spatial scenarios. Attached Figure Description

[0011] Figure 1 This is a flowchart illustrating the substation noise coupling simulation analysis method based on a sound ray simulation scenario provided in this embodiment of the invention. Figure 2 This is a schematic diagram of the substation noise coupling simulation analysis device based on a sound ray simulation scenario provided in an embodiment of the present invention; Figure 3 An embodiment diagram of the electronic device provided in this invention; Figure 4 An embodiment diagram of a computer-readable storage medium provided in accordance with the present invention. Detailed Implementation

[0012] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0013] In the description of this invention, the terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of indicated technical features. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more of the stated features. In the description of this invention, "a plurality of" means two or more, unless otherwise explicitly specified.

[0014] In the description of this invention, the term "for example" is used to mean "used as an example, illustration, or description." Any embodiment described as "for example" in this invention is not necessarily to be construed as being more preferred or advantageous than other embodiments. The following description is provided to enable any person skilled in the art to make and use the invention. Details are set forth in the following description for purposes of explanation. It should be understood that those skilled in the art will recognize that the invention can be made without using these specific details. In other instances, well-known structures and processes will not be described in detail to avoid obscuring the description of the invention with unnecessary detail. Therefore, the invention is not intended to be limited to the embodiments shown, but is consistent with the broadest scope of the principles and features disclosed herein.

[0015] Optional, see below Figure 1 , Figure 1 This is a flowchart illustrating the substation noise coupling simulation analysis method based on a sound ray simulation scenario provided by the present invention. In this embodiment of the invention, the execution entity of the substation noise coupling simulation analysis method based on a sound ray simulation scenario is a noise coupling simulation analysis device. Therefore, the substation noise coupling simulation analysis method based on a sound ray simulation scenario includes: Step 10: Construct a three-dimensional spatial scene based on the spatial location and geometric dimensions of the spatial entities within the substation.

[0016] Optionally, the noise coupling simulation analysis device collects the spatial location and geometric dimension information of all spatial entities within the substation. These spatial entities include electrical equipment (such as transformers and reactors), buildings (such as control rooms and duty rooms), and obstructions (such as walls and guardrails). Furthermore, based on the above information, the noise coupling simulation analysis device will construct a three-dimensional spatial scene in the computer system that is consistent with the actual spatial structure of the substation. During the construction process, each spatial entity needs to be accurately modeled, including the entity's shape, size, and position coordinates, to ensure that the three-dimensional spatial scene can realistically reflect the relative positions and spatial relationships of each spatial entity within the substation.

[0017] In one embodiment, the substation contains a transformer of model S11-2000 / 10 with spatial coordinates (10, 5, 0) and geometric dimensions of 5 meters long, 3 meters wide, and 2.5 meters high; a control room building with spatial coordinates (20, 10, 0) and geometric dimensions of 15 meters long, 8 meters wide, and 6 meters high; and a section of wall serving as a barrier with spatial coordinates from (0, 0, 0) to (30, 0, 2) and geometric dimensions of 30 meters long, 0.5 meters wide, and 2 meters high. Therefore, for the transformer, a cuboid model is constructed in the three-dimensional coordinate system with (10, 5, 0) as the center of the base, and dimensions of 5 meters long, 3 meters wide, and 2.5 meters high to represent the transformer. For the control room building, a cuboid model is constructed with (20, 10, 0) as the center of the base, and dimensions of 15 meters long, 8 meters wide, and 6 meters high. For the perimeter wall, a cuboid model is constructed from coordinates (0, 0, 0) to (30, 0, 2), with dimensions of 30 meters long, 0.5 meters wide, and 2 meters high to serve as a barrier. Through the modeling of spatial entities, a three-dimensional spatial scene of the substation is finally constructed, which can accurately present the position and geometric features of spatial entities such as the transformer, control room, and perimeter wall.

[0018] Step 20: For each noise source within the substation, control the noise source to emit multiple sound rays at a preset angle, and determine the collision position of each sound ray with spatial entities as it propagates in the three-dimensional space scene, as well as the reflection direction after each collision, until the sound ray propagates to the boundary of the three-dimensional space scene. Each sound ray carries the initial acoustic parameters of the noise source.

[0019] Furthermore, the noise coupling simulation analysis device identifies each noise source within the substation, which is typically electrical equipment within the substation, such as transformers and reactors.

[0020] Furthermore, for each noise source, the noise coupling simulation analysis device controls the noise source to emit multiple sound rays according to a preset angle (e.g., every 10 degrees is an angle interval). Each sound ray carries the initial acoustic parameters of the noise source, including sound pressure level and frequency spectrum.

[0021] During the propagation of sound rays, the collision between each sound ray and spatial entities in the three-dimensional space scene is monitored in real time to determine the collision location. The reflection direction after each collision is calculated based on the surface characteristics of the spatial entities and the collision angle, and this process continues until the sound ray propagates to the boundary of the three-dimensional space scene.

[0022] Continuing with the aforementioned three-dimensional spatial scenario, the transformer serves as the noise source, emitting sound rays at preset 10-degree intervals. The transformer's initial acoustic parameters are a sound pressure level of 85 dB and a frequency spectrum ranging from 50 Hz to 2000 Hz. A noise coupling simulation analysis device controls the transformer to emit sound rays at preset angles, for example, starting from 0 degrees and emitting one ray every 10 degrees, for a total of 36 sound rays. One of these sound rays is emitted from a 30-degree angle, with an initial acoustic parameter of 85 dB and a frequency spectrum ranging from 50 Hz to 2000 Hz. During its propagation, this sound ray first collides with the wall of the control room building at coordinates (20, 10, 0). The noise coupling simulation analysis device detects the collision location at coordinates (18, 8, 2) and calculates the reflection direction as 45 degrees to the normal direction based on the wall's material (assumed to be concrete with a reflection coefficient of 0.8) and the collision angle. The sound then continued to propagate, collided with the wall, determined the collision location and a new reflection direction, until it stopped at the boundary of the three-dimensional space scene (such as the virtual boundary coordinates (50, 50, 10) of the substation).

[0023] Step 30: Based on the propagation path, collision location, and reflection direction of each sound ray, determine the changes in acoustic parameters of each sound ray during propagation due to distance attenuation and collision reflection, and generate an acoustic parameter change sequence for each sound ray. The acoustic parameter change sequence includes the acoustic parameter values ​​and corresponding spatial coordinates of each node along the propagation path of the sound ray.

[0024] Furthermore, for each sound ray, the noise coupling simulation analysis device calculates the changes in acoustic parameters during propagation based on its propagation path, collision location, and reflection direction. The changes in acoustic parameters of each sound ray during propagation due to distance attenuation and collision reflection generate a sequence of acoustic parameter changes for each sound ray, as described in steps 301 to 304. The sequence of acoustic parameter changes includes the acoustic parameter values ​​and corresponding spatial coordinates of each node (such as emission point, collision point, boundary point, etc.) on the propagation path of the sound ray.

[0025] Step 40: For each spatial point in the three-dimensional spatial scene, couple the target acoustic parameter change sequence based on the target sound ray passing through the spatial point through the propagation path to determine the noise coupling result of each spatial point.

[0026] Furthermore, the noise coupling simulation analysis device selects each spatial point in the three-dimensional spatial scene and determines the target sound rays that all propagation paths pass through that spatial point.

[0027] Furthermore, for each target acoustic ray, the noise coupling simulation analysis device extracts the acoustic parameters corresponding to the spatial point in the target acoustic parameter change sequence, and couples the acoustic parameters of the target acoustic ray at the spatial point to obtain the noise coupling result of each spatial point, as detailed in steps 401 to 404.

[0028] Step 50: Traverse the noise coupling results of each spatial point in the three-dimensional spatial scene to obtain the noise coupling distribution of the space where the substation is located.

[0029] Furthermore, the noise coupling simulation analysis device traverses each spatial point in the three-dimensional spatial scene according to a certain rule (such as according to the order of grid division) and collects the noise coupling results of each spatial point. Further, the noise coupling simulation analysis device organizes and statistically analyzes the above results to form the noise coupling distribution of the space where the substation is located. Optionally, embodiments of the present invention can present this in the form of noise distribution cloud maps, data tables, etc., to intuitively display the noise magnitude and frequency distribution characteristics at different spatial locations.

[0030] Continuing with the above embodiments, the noise coupling simulation analysis device divides the constructed three-dimensional spatial scene of the substation into a grid of 1m x 1m x 0.5m, traversing the spatial points corresponding to each grid node. For example, the noise coupling result of spatial point (5, 3, 0.5) is a sound pressure level of 68 dB, the noise coupling result of spatial point (25, 10, 3) is a sound pressure level of 55 dB, and the noise coupling result of spatial point (40, 40, 5) is a sound pressure level of 40 dB, etc. The noise coupling simulation analysis device organizes the noise coupling results of all spatial points and generates a noise coupling distribution cloud map using computer software. Different colors in the cloud map represent different sound pressure level ranges, such as red indicating a sound pressure level above 70 dB, yellow indicating 60-70 dB, green indicating 50-60 dB, and blue indicating below 50 dB. Simultaneously, a data table is generated, recording the coordinates of each spatial point and its corresponding sound pressure level and frequency spectrum. These methods reveal the noise coupling distribution within the substation's space, clearly demonstrating the noise distribution pattern: areas near noise sources such as transformers experience higher noise levels, while areas far from noise sources and shielded by buildings experience lower noise levels.

[0031] This invention, by emitting multiple sound rays and tracking their propagation paths, can comprehensively record the reflection and attenuation of noise during propagation, covering all possible propagation paths within the substation, no longer limited to fixed point locations. Furthermore, by tracking the acoustic parameter change sequence of each sound ray carrying the initial acoustic parameters of the corresponding noise source as it propagates at various nodes, the noise coupling result at each spatial point is determined by coupling based on the acoustic parameter change sequence of all propagation paths passing through that spatial point, accurately reflecting the coupling result of different noise sources at that spatial point. Moreover, when the equipment's operating state changes dynamically, the acoustic parameters of the noise source can be updated, sound rays can be re-emitted, and noise coupling can be performed to quickly obtain a new noise coupling distribution, thus enabling dynamic response to changes in equipment state. Therefore, this invention, through full-path sound ray tracking, acoustic parameter change recording, and full-space-point coupling, can accurately reflect the dynamic propagation and coupling characteristics of noise in complex spatial scenarios.

[0032] In one embodiment, steps 301 to 304 include: Step 301: For each sound ray, based on the collision position and reflection direction of the sound ray, extract each node on its propagation path, and map and bind the three-dimensional spatial coordinates of each node and its ordinal position in the propagation path to obtain the coordinate mapping result of each node. Nodes include collision points and non-collision points.

[0033] Optionally, for each sound ray, the noise coupling simulation analysis device sorts out its complete propagation path based on the collision position and reflection direction of the sound ray, and extracts each node from the propagation path, as in steps 3011 to 3014. The nodes include collision points (the position where the sound ray collides with the spatial entity) and non-collision points (such as the emission point of the sound ray, the midpoint on the propagation path, the endpoint reaching the boundary of the three-dimensional spatial scene, etc.).

[0034] Furthermore, the noise coupling simulation analysis device maps and binds the three-dimensional spatial coordinates of each node to its ordinal identifier in the propagation path (i.e., the sequential number of the node in the propagation path) to form the coordinate mapping result of each node, thereby clarifying the positional relationship of each node in the sound propagation path.

[0035] Continuing with the example of the sound emitted from the S11-2000 / 10 transformer (coordinates (10, 5, 0)) at a 30-degree angle in step 20, its propagation path is as follows: starting from the transformer, it first collides with the wall of the control room building (coordinates (20, 10, 0)), then with the surrounding wall, and finally reaches the boundary of the three-dimensional space scene. Therefore, the nodes extracted by the noise coupling simulation analysis device include: the emission point (non-collision point), the collision point with the wall of the control room building (collision point), the collision point with the surrounding wall (collision point), and the endpoint reaching the boundary of the three-dimensional space scene (non-collision point). The three-dimensional space coordinates of each node are (10, 5, 0), (18, 8, 2), (25, 0, 1), and (50, 50, 10), respectively. The sequence identifiers are 1, 2, 3, and 4 in the order of propagation. Furthermore, the noise coupling simulation analysis device maps and binds the three-dimensional spatial coordinates of the above nodes with the sequence identifiers to obtain the coordinate mapping results: sequence identifier 1 corresponds to coordinates (10, 5, 0), sequence identifier 2 corresponds to coordinates (18, 8, 2), sequence identifier 3 corresponds to coordinates (25, 0, 1), and sequence identifier 4 corresponds to coordinates (50, 50, 10).

[0036] Step 302: For each collision point, based on the sound energy absorption coefficient and frequency response coefficient of the material of the collision interface passed through by the collision point, and combined with the dominant frequency and incident angle of each sound ray at the collision point, determine the reflection attenuation coefficient of each sound ray at the collision point. The dominant frequency is determined based on the frequency spectrum of each sound ray at the collision point, and the incident angle is determined based on the reflection direction and incident direction of each sound ray at the collision point.

[0037] Furthermore, for each collision point of each sound ray, the noise coupling simulation analysis device determines the material of the collision interface (such as concrete, steel, wood, etc.) that the collision point passes through, and obtains the sound energy absorption coefficient and frequency response coefficient of the material.

[0038] Furthermore, the noise coupling simulation analysis device determines the dominant frequency (usually the frequency with the highest energy in the frequency spectrum) at the collision point from the frequency spectrum carried by the sound ray, and calculates the incident angle (the angle between the incident direction and the normal direction of the collision interface) based on the incident direction and reflection direction of the sound ray at the collision point.

[0039] Furthermore, based on the obtained acoustic absorption coefficient, frequency response coefficient, dominant frequency, and incident angle, the noise coupling simulation analysis device determines the reflection attenuation coefficient of the sound ray at the collision point using acoustic calculation methods. The reflection attenuation coefficient characterizes the degree of attenuation of acoustic parameters after the sound ray is reflected upon collision. In this embodiment of the invention, the formula for calculating the reflection attenuation coefficient is as follows: ,in, Let be the reflection attenuation coefficient of the k-th collision node. Let m be the sound energy absorption coefficient of the collision interface material (range 0-1, where 0 represents complete reflection and 1 represents complete absorption). For the collision interface material m at the main frequency The frequency response coefficient (reflects the material's reflection characteristics at a specific frequency, and its value is usually between 0 and 1). Let be the incident angle of the sound ray at the k-th collision node (the angle between the incident direction and the normal direction of the collision interface).

[0040] Continuing with the example of the sound ray in step 301, the collision point (coordinates (18, 8, 2)) between the sound ray and the wall of the control room is traversed by a collision interface material of concrete m. According to research, the sound energy absorption coefficient of this concrete m material is 0.2, and its frequency response coefficient is 0.95 in the range of 50 Hz to 2000 Hz. In the frequency spectrum of this sound ray at the collision point, the energy is highest at 100 Hz, therefore the dominant frequency is 100 Hz. Through calculation, the angle between the incident direction of the sound ray at the collision point and the normal direction of the wall (incident angle) is 30 degrees. Therefore, the reflection attenuation coefficient at the collision node is... The value is (1-0.2)*0.95*cos30≈0.659.

[0041] Step 303: For each non-collision node, based on the straight-line distance between the non-collision node and its next node, and combined with the dominant frequency wavelength corresponding to each sound ray in each propagation path, determine the distance attenuation coefficient from the non-collision node to its next node.

[0042] Furthermore, for each non-collision node of each sound ray, the noise coupling simulation analysis device calculates the straight-line distance between the non-collision node and its next node (which can be a collision node or a non-collision node). Further, the dominant frequency of the sound ray is determined from the frequency spectrum of the non-collision node, and the corresponding dominant frequency wavelength is calculated based on the dominant frequency. The formula for calculating the dominant frequency wavelength in this embodiment is: λ n =c / f n , where λ n Let f be the wavelength (in meters) of the dominant frequency of the sound ray at the nth non-collision node, c be the speed of sound (in meters per second, taken as 340 meters per second), and f be the wavelength (in meters per second). n The dominant frequency (Hertz) of the sound ray at the nth non-collision node.

[0043] Furthermore, based on the calculated straight-line distance and dominant frequency wavelength, the distance attenuation coefficient from the non-collision node to its next node is determined using an acoustic formula. The distance attenuation coefficient reflects the attenuation of acoustic parameters due to distance along the propagation path of the sound ray. The formula for calculating the distance attenuation coefficient in this embodiment of the invention is: D nk =(λ n / (4πd nk ))², where Dnk Let d be the distance decay coefficient from the nth non-collision node to the kth node. nk It is the straight-line distance (in meters) from the nth non-collision node to the kth node.

[0044] Continuing with the example of the sound ray in step 301, the next node after its emission point (coordinates (10, 5, 0), non-collision node, sequence identifier 1) is the collision point with the wall of the control room building (coordinates (18, 8, 2), sequence identifier 2). The noise coupling simulation analysis device calculates the straight-line distance between these two nodes: according to the formula for the distance between two points in space, the distance d ≈ 8.77 meters. The dominant frequency of the sound ray in the frequency spectrum at the emission point is 100 Hz, and the wavelength of the dominant frequency λ = 340 m / s ÷ 100 Hz = 3.4 meters. Based on the straight-line distance of 8.77 meters and the dominant frequency wavelength of 3.4 meters, the noise coupling simulation analysis device calculates that the distance attenuation coefficient from the non-collision node to the next node is 0.38. For example, the node preceding the endpoint of the sound ray at the boundary of the three-dimensional space scene (coordinates (50, 50, 10), non-collision node, sequence identifier 4) is the collision point with the wall (coordinates (25, 0, 1), sequence identifier 3). The straight-line distance between the two nodes is calculated to be 56.62 meters. The main frequency at this node is still 100 Hz, and the wavelength is 3.4 meters. The distance attenuation coefficient is calculated to be 0.05.

[0045] Step 304: Based on the reflection attenuation coefficient, distance attenuation coefficient, and the initial acoustic parameters and coordinate mapping results of each node, an analysis is performed to generate a sequence of acoustic parameter changes for each sound ray.

[0046] Furthermore, the noise coupling simulation analysis device analyzes the initial acoustic parameters and coordinate mapping results of each node based on the reflection attenuation coefficient and distance attenuation coefficient, and generates a sequence of acoustic parameter changes for each sound ray, as described in steps 3041 to 3044.

[0047] The embodiments of the present invention can accurately track the changes in acoustic parameters of each sound ray during propagation and associate the changes in acoustic parameters with specific spatial locations to generate a complete and accurate sequence of acoustic parameter changes. This provides detailed basic data for noise coupling calculations at spatial points in subsequent steps, thereby ensuring the accuracy and reliability of the noise coupling results.

[0048] In one embodiment, steps 3011 to 3014 include: Step 3011: Based on the initial emission direction, the first entity boundary of the spatial entity, and the starting coordinates of the sound ray, determine the initial spatial domain that intersects with the first entity boundary during the initial propagation phase.

[0049] Optionally, for each sound ray, the noise coupling simulation analysis device determines the initial emission direction of the sound ray (i.e., the propagation direction vector of the sound ray when it originates from the noise source) and the starting point coordinates (i.e., the spatial coordinates of the noise source). Further, the noise coupling simulation analysis device traverses the first entity boundaries (boundaries of spatial entity surfaces, such as the outer shell surface of a transformer, the walls of a building, etc.) of all spatial entities in the three-dimensional spatial scene of the substation, and determines through geometric calculations which first entity boundaries the sound ray might intersect with during the initial propagation stage (the stage before any collision occurs after originating from the starting point). The spatial range of these intersections is the initial spatial domain. It should be noted that the initial spatial domain in this embodiment of the invention is a virtual spatial region used to define the range of entity boundaries where the sound ray may first collide.

[0050] Optionally, the acoustic parameter equation in this embodiment of the invention is as follows: ,in, Let be the spatial coordinate vector of the sound ray at propagation parameter t. Let t be the coordinate vector of the starting point of the sound ray, and t be the propagation parameter (t≥0). is the initial emission direction unit vector of the sound ray.

[0051] Continuing with the example of the sound ray emitted from the S11-2000 / 10 transformer (starting coordinates (10, 5, 0)) at a 30-degree angle in step 20, its initial emission direction vector is calculated to be (cos30°, sin30°, 0) ≈ (0.866, 0.5, 0). The noise coupling simulation analysis device traverses the spatial entities in the three-dimensional space scene of the substation and finds that the first entity boundaries that the sound ray may intersect in the initial propagation stage include the east wall of the control room building (coordinate range x=20, y=5 to 15, z=0 to 6) and the north wall of the perimeter wall (coordinate range y=0, x=0 to 30, z=0 to 2). The noise coupling simulation analysis device determines through geometric modeling that the initial spatial domain formed by the intersection area of ​​these two first entity boundaries and the sound ray is a spatial region extending along the initial emission direction with the starting point (10, 5, 0) as the vertex, containing the above two walls. The spatial points within this region satisfy the geometric intersection conditions of the sound ray propagation path and the wall boundaries.

[0052] Step 3012: Based on the collision position of each first entity boundary and the sound ray within the initial spatial domain, determine the theoretical collision distance between each first entity boundary and the sound ray, and based on the theoretical collision distance, concavity and convexity, and curvature of each first entity boundary, determine the collision priority of each first entity boundary.

[0053] Furthermore, for each first entity boundary within the initial spatial domain, the noise coupling simulation analysis device, based on the initial emission direction and starting coordinates of the sound ray, and combined with the geometric equations of the first entity boundary (such as plane equations, surface equations, etc.), calculates the theoretical collision distance (i.e., the straight-line distance from the starting coordinates to the collision point) of the sound ray to that first entity boundary by solving the intersection points of the spatial straight line and the entity boundary. In this embodiment of the invention, the theoretical collision distance is calculated as follows: For a planar first entity boundary, its plane equation is... , It is a plane normal vector. For constant terms, Let be the coordinate vector of any point on the plane, and be the propagation parameter at the intersection of the sound ray and the plane. The theoretical collision distance D = | |

[0054] Furthermore, the noise coupling simulation analysis device acquires the concavity (characterizing the degree of concavity and convexity of the boundary surface, with positive values ​​for convex surfaces, negative values ​​for concave surfaces, and 0 values ​​for planes) and curvature (characterizing the degree of bending of the boundary surface, with plane curvature being 0 and surface curvature being positive and the larger the value, the more obvious the bending) of each first entity boundary. Based on the theoretical collision distance, concavity and convexity, and curvature, the device determines the collision priority of each first entity boundary through the set priority calculation rules (such as the shorter the theoretical collision distance, the higher the priority; the convex boundary with moderate curvature has a higher priority than the concave boundary or the boundary with excessive curvature).

[0055] Optionally, the collision priority calculation in this embodiment of the invention is as follows: ,in, This represents the collision priority of the kth first entity boundary (the smaller the value, the higher the priority). , , Weighting coefficients ( + + =1), Let be the theoretical collision distance of the k-th first entity boundary. The concavity / convexity of the k-th first entity boundary (normalized value). The curvature of the k-th first entity boundary (normalized value).

[0056] Continuing with the above embodiment, the first entity boundary within the initial spatial domain is the east wall of the control room building and the north wall of the perimeter wall. For the east wall of the control room building, its planar equation is x=20, and the parametric equations for the sound ray are x=10+0.866t, y=5+0.5t, z=0 (t is a parameter representing the propagation time or distance ratio). Solving the simultaneous equations yields t=(20-10) / 0.866≈11.55, and the theoretical collision distance is t*1 (since the coefficient of parameter t is the magnitude of the direction vector, which is 1 here)≈11.55 meters. This wall surface is planar, with zero concavity and zero curvature. For the north wall of the enclosure, the plane equation is y=0. Solving the simultaneous acoustic parametric equations, we get 5+0.5t=0, t=-10 (negative values ​​indicate the area behind the starting point and have no practical significance). After correction, the effective theoretical collision distance is the distance from the starting point to the actual point of intersection on the north wall, which is recalculated to be 15 meters. This wall is planar with zero concavity and zero curvature. Based on the theoretical collision distance, the noise coupling simulation analysis device determines that the collision priority of the east wall of the control room building is higher than that of the north wall (because 11.55 meters < 15 meters), with collision priority values ​​of 1 (highest) and 2 respectively.

[0057] Step 3013: Based on the initial launch direction, starting point coordinates, and the first boundary normal vector and first position parameters of the first target entity with the highest collision priority, determine the first spatial coordinates of the first collision point.

[0058] Furthermore, the noise coupling simulation analysis device selects the first entity boundary with the highest collision priority from all first entity boundaries as the first target entity boundary.

[0059] Furthermore, the noise coupling simulation analysis device obtains the first boundary normal vector (a vector perpendicular to the surface of the entity boundary pointing outward) and the first position parameters (such as the plane equation parameters of the boundary, the coordinates of the curvature center of the surface, etc.) of the first target entity boundary.

[0060] Furthermore, the noise coupling simulation analysis device, by combining the initial emission direction and starting point coordinates of the sound ray with the intersection of the spatial straight line and the entity boundary, determines the first spatial coordinates of the first collision point between the sound ray and the first target entity boundary.

[0061] Continuing with the above embodiment, the boundary of the first target entity with the highest collision priority is the east wall of the control room building. The first boundary normal vector of this wall is (1, 0, 0) (perpendicular to the wall pointing outwards), and the first position parameters are the plane equations x=20, y∈[5, 15], z∈[0, 6]. The parametric equations of the sound ray are x=10+0.866t, y=5+0.5t, z=0. Substituting x=20 into the parametric equations of the sound ray, we get t=(20-10) / 0.866≈11.55. At this time, y=5+0.5*11.55≈5+5.775=10.775, z=0. The coordinates (20, 10.775, 0) are within the y and z coordinate range of the wall. Therefore, the first spatial coordinates of the first collision point are determined to be (20, 10.775, 0).

[0062] Step 3014: Based on the first boundary normal vector, the first spatial coordinates, and the reflection direction, determine each node of the sound ray on the propagation path.

[0063] Furthermore, the noise coupling simulation analysis device determines each node of the sound ray on the propagation path based on the first boundary normal vector, the first spatial coordinates, and the reflection direction, as specifically in steps 30141 to 30144.

[0064] The embodiments of the present invention can accurately track all nodes in the process of each sound ray from emission to propagation to the boundary of the three-dimensional spatial scene, and clarify the spatial coordinates of each node and its order relationship in the propagation path, laying a precise spatial position foundation for subsequent calculation of the acoustic parameter change sequence of the sound ray.

[0065] In one embodiment, the process of steps 30141 to 30144 includes: Step 30141: Based on the first boundary normal vector and the reflection direction, perform reflection reconstruction to obtain the reconstructed reflection direction at the first collision point.

[0066] Optionally, the noise coupling simulation analysis device performs reflection reconstruction based on the first boundary normal vector (the normal vector of the entity boundary at the first collision point) and the initially calculated reflection direction. Reflection reconstruction corrects the initial reflection direction, eliminating deviations caused by microscopic unevenness or material inhomogeneity on the entity boundary surface, ensuring that the reconstructed reflection direction strictly follows the law of reflection (the incident angle and the reflection angle are on opposite sides of the normal vector and are equal in angle, and all three are coplanar). Specifically, the noise coupling simulation analysis device verifies whether the initial reflection direction, the first boundary normal vector, and the incident direction satisfy the geometric relationship of the law of reflection. If deviations exist, adjustments are made to ultimately obtain the reconstructed reflection direction at the first collision point.

[0067] Optionally, the formula for verifying the reconstructed reflection direction in this embodiment of the invention is as follows: = ,in, Angle of incidence The angle of reflection, Let be the unit vector of the incident direction. To reconstruct the unit vector of the reflection direction, It is the unit vector of the first boundary normal.

[0068] Continuing with the above embodiment, the first collision point is (20, 10.775, 0), the first boundary normal vector is (1, 0, 0) (pointing outwards from the entity), the incident direction vector of the sound ray is (0.866, 0.5, 0), and the initially calculated reflection direction vector is (-0.866, 0.5, 0). The noise coupling simulation analysis device verifies whether this reflection direction satisfies the law of reflection: the cosine of the angle between the incident direction and the normal vector is (0.866*1+0.5*0+0*0)=0.866 (corresponding to 30°), and the cosine of the angle between the reflection direction and the normal vector is (-0.866*1+0.5*0+0*0)=-0.866 (absolute value corresponding to 30°), and the incident direction, reflection direction, and normal vector are coplanar (all in the z=0 plane), which conforms to the law of reflection. Therefore, the noise coupling simulation analysis device determines that the reflection direction after reconstruction is still (-0.866, 0.5, 0).

[0069] Step 30142: Based on the reconstructed reflection direction, the second entity boundary of the spatial entity that did not collide, and the first spatial coordinates, determine the secondary spatial domain where the reflected sound ray intersects with the second entity boundary.

[0070] Furthermore, the noise coupling simulation analysis device uses the reconstructed reflection direction as the new propagation direction and the first spatial coordinates (coordinates of the first collision point) as the new starting point. It traverses the second entity boundaries (i.e., entity boundaries other than the entity to which the first collision point belongs, such as walls, other building walls, etc.) of all non-collision spatial entities in the three-dimensional space scene of the substation. Through geometric calculations, it determines which second entity boundaries the reflected sound ray might intersect with; the spatial range of these intersections is the secondary spatial domain. The secondary spatial domain is used to define the range of entity boundaries that the sound ray might collide with during the propagation phase after the first collision.

[0071] Optionally, the boundary equation of the quadratic spatial domain in this embodiment of the invention is: for the planar second entity boundary, its intersection condition with the sound ray is: (t≥0), where, The second boundary normal vector, Let be the first spatial coordinate vector. To reconstruct the unit vector of the reflection direction, t is the constant term of the plane equation, and t is the propagation parameter.

[0072] Continuing with the above embodiment, the reconstructed reflection direction vector is (-0.866, 0.5, 0), and the first spatial coordinates are (20, 10.775, 0). The spatial entities that did not collide include the enclosure wall (the second entity boundary is its north wall, with the plane equation y=0, x∈[0, 30], z∈[0, 2]) and the outer surface of another reactor. The noise coupling simulation analysis device calculates t=(0-10.775) / 0.5=-21.55 (a negative value, which needs to be corrected to the effective intersection point in the direction of sound ray propagation) by simultaneously solving the acoustic ray parameter equation (x=20-0.866t, y=10.775+0.5t, z=0, t≥0) and the equation y=0 of the north wall of the enclosure wall. After recalculation, the effective intersection point is on the sound ray propagation path; therefore, the north wall of the enclosure wall belongs to the secondary spatial domain. The secondary spatial domain is the area extending from (20, 10.775, 0) along the direction of (-0.866, 0.5, 0), including the space of the north wall of the enclosure.

[0073] Step 30143: Based on the reconstructed reflection direction, the first spatial coordinates, the second boundary normal vector and the second position parameters of the second target entity boundary with the highest collision priority in the secondary spatial domain, determine the second spatial coordinates of the secondary collision point.

[0074] Furthermore, the noise coupling simulation analysis device selects the second target entity boundary with the highest collision priority from the second entity boundary in the secondary spatial domain according to the collision priority rules (theoretical collision distance, concavity and convexity, curvature) similar to those in step 3012, and obtains the second boundary normal vector (the normal vector of the boundary surface) and the second position parameters (such as plane equation parameters) of the second target entity boundary.

[0075] Furthermore, the noise coupling simulation analysis device determines the second spatial coordinates of the secondary collision point by combining the second boundary normal vector, the second position parameters, the reconstructed reflection direction, and the first spatial coordinates, and by solving the intersection point of the sound ray and the boundary of the second target entity.

[0076] Continuing with the above embodiment, the second entity boundary within the secondary spatial domain is the north wall of the enclosure, with its second boundary normal vector being (0, -1, 0) (pointing outward in the negative y-axis direction). The second position parameter is the plane equation y = 0, x ∈ [0, 30], z ∈ [0, 2]. The acoustic ray parameter equations are x = 20 - 0.866t, y = 10.775 + 0.5t, z = 0. Solving the equations simultaneously, we get 10.775 + 0.5t = 0 → t = -21.55 (negative values ​​indicate opposite directions, requiring correction to the propagation direction corresponding to the absolute value. The intersection point in the actual acoustic ray propagation direction is calculated as follows: along the reflection direction, the acoustic ray y-coordinate needs to be reduced, hence the effective t is 21.55). At this point, x = 20 - 0.866 * 21.55 ≈ 20 - 18.66 = 1.34, z = 0, and the coordinates (1.34, 0, 0) are within the x and z range of the north wall of the enclosure. The theoretical collision distance at this boundary is 21.55 meters (t value), with 0 concavity / convexity and 0 curvature, and the highest collision priority. Therefore, the noise coupling simulation analysis device determines the second spatial coordinates of the secondary collision point as (1.34, 0, 0).

[0077] Step 30144: Based on the spatial coordinates of the start point, end point, and non-collision points in the sound ray, combined with the first and second spatial coordinates, sort the start point, end point, first collision point, second collision point, and non-collision points according to spatial coordinate order to obtain the nodes of the sound ray on the propagation path. Non-collision points are the nodes between the start point and end point excluding the first and second collision points.

[0078] Furthermore, the noise coupling simulation analysis device first clarifies the starting point coordinates (noise source coordinates), ending point coordinates (three-dimensional spatial scene boundary points), and non-collision points (preset intermediate nodes between the starting point and the ending point, excluding collision points, such as sampling points set at certain intervals).

[0079] Furthermore, the noise coupling simulation analysis device combines the first spatial coordinates (first collision point) and the second spatial coordinates (secondary collision point), and sorts each point according to the order of their positions on the sound propagation path (by calculating the propagation parameter t value corresponding to each point; the larger the t value, the later the position). Finally, it obtains each node on the sound propagation path, including the starting point, all collision points, non-collision points, and the ending point.

[0080] Optionally, the node sorting parameters in this embodiment of the invention are calculated as: the propagation parameters corresponding to each node. satisfy (From the starting point to the first collision point) or (After the first collision point), among which, Let be the coordinate vector of the k-th node. The propagation parameters for the first collision point, This is the initial launch direction vector. This is the reconstructed reflection direction vector. Nodes are arranged according to... Sort by size from smallest to largest.

[0081] Continuing with the above embodiment, the starting coordinates of the sound ray are (10, 5, 0) (t=0). The non-collision points are preset to be points on the propagation path at t=5 and t=30, with corresponding coordinates as follows: at t=5, x=10+0.866*5≈14.33, y=5+0.5*5=7.5, z=0→(14.33, 7.5, 0); at t=30, propagating along the reflection direction, x=20-0.866*(30-11.55)≈20-0.866*18.45≈20-16≈4, y=10.775+0.5*(30-11.55)≈10.775+9.225≈20→(4, 20, 0) (t=11.55 is the propagation parameter corresponding to the first collision point). The ending coordinates are (50, 50, 10) (t=100, assumed value). The propagation parameter t values ​​for each point are as follows: starting point t=0, non-collision point 1 t=5, first collision point t=11.55, second collision point t=11.55+21.55=33.1, non-collision point 2 t=30 (actually after the first collision point and before the second collision point, the sorting logic needs to be corrected according to the actual path order), and ending point t=100. After correction, the order by t value is: starting point (10, 5, 0) → non-collision point 1 (14.33, 7.5, 0) → first collision point (20, 10.775, 0) → non-collision point 2 (4, 20, 0) → second collision point (1.34, 0, 0) → ending point (50, 50, 10).

[0082] The embodiments of the present invention can accurately determine the reflection direction of the sound ray after the first collision, the secondary collision point, and the order of all nodes, and completely construct the node sequence on the sound ray propagation path, ensuring that the spatial coordinates and order relationship of each node are accurate, providing a coherent and precise spatial path basis for the generation of subsequent acoustic parameter change sequences.

[0083] In one embodiment, the ancient city in steps 3041 to 3045 includes: Step 3041: Based on the initial sound pressure level and distance attenuation coefficient of each sound ray at the non-collision node, and combined with the angle between each sound ray in each propagation path and the reference axis in the three-dimensional space scene, determine the first sound pressure level attenuated by distance from the non-collision node to its next node.

[0084] Optionally, for each non-collision node of each sound ray, the noise coupling simulation analysis device obtains the initial sound pressure level of that node (i.e., the sound pressure level value carried by the node itself) and the distance attenuation coefficient from the non-collision node to its next node. Simultaneously, it calculates the angle between the sound ray's propagation path and a reference axis (such as the X, Y, or Z axis) in the three-dimensional spatial scene. This angle is used to correct for the effect of distance attenuation on the sound pressure level (since sound wave propagation attenuation may differ in different directions). Further, based on the initial sound pressure level, the distance attenuation coefficient, and the aforementioned angle, the noise coupling simulation analysis device calculates the first sound pressure level caused by distance attenuation from the non-collision node to its next node using a specific formula.

[0085] The formula for calculating the first sound pressure level in this embodiment of the invention is as follows: ,in, This is the first sound pressure level (decibels). The initial sound pressure level (decibels) for non-collision nodes. This is the distance attenuation coefficient. The angle between the sound ray propagation path and the i-th reference axis (i=1, 2, 3 correspond to the X-axis, Y-axis, and Z-axis, respectively).

[0086] Continuing with the above embodiment, taking a non-collision node A (coordinates (10, 5, 0)) of a certain sound ray as an example, its initial sound pressure level is 85 dB, and the distance attenuation coefficient from node A to the next node B (collision point) is 0.38. The sound ray makes an angle of 30 degrees with the X-axis, 60 degrees with the Y-axis, and 90 degrees with the Z-axis in this propagation path of the three-dimensional scene (because the propagation path is in the horizontal plane). The noise coupling simulation analysis device, combining these parameters, calculates the first sound pressure level attenuated by distance from node A to node B. The calculated first sound pressure level is 85 dB * 0.38 * cos30° ≈ 28.5 dB.

[0087] Step 3042: Based on the first sound pressure level of each sound ray before reaching the collision point, combined with the reflection attenuation coefficient and reflection angle of each sound ray at the collision point, determine the second sound pressure level after reflection at each collision point. The reflection angle is determined based on the reflection direction of each sound ray at the collision point.

[0088] Furthermore, for each collision point of each sound ray, the noise coupling simulation analysis device obtains the first sound pressure level before the sound ray reaches the collision point, and combines the reflection attenuation coefficient and reflection angle (the angle between the reflection direction and the normal direction of the collision interface) of the collision point. The reflection angle is used to correct the influence of the difference in sound energy distribution during reflection on the sound pressure level, and to determine the second sound pressure level of the collision point after reflection.

[0089] Optionally, the formula for calculating the second sound pressure level in this embodiment of the invention is: ,in, This is the second sound pressure level (decibels). The first sound pressure level (decibels) before reaching the point of impact. The reflection attenuation coefficient is... The reflection angle is in radians.

[0090] Continuing with the above embodiment, the first sound pressure level when the sound ray reaches the collision point B is 28.5 dB. The reflection attenuation coefficient of the collision point is 0.75, and the reflection angle is 45 degrees. Therefore, the second sound pressure level is 28.5 dB * 0.75 * sin45° ≈ 14.5 dB. This value is the second sound pressure level after reflection at the collision point B.

[0091] Step 3043: Based on the initial frequency spectrum and distance attenuation coefficient of each sound ray at the non-collision node, and combined with the angular frequency of each sound ray in each propagation path, determine the first frequency spectrum of each sound ray after distance attenuation from the non-collision node to its next node.

[0092] Furthermore, for each non-collision node of each sound ray, the noise coupling simulation analysis device acquires the initial frequency spectrum of that node (containing the sound energy distribution corresponding to different frequencies) and the distance attenuation coefficient from that non-collision node to its next node. Simultaneously, it calculates the angular frequency of the sound ray along this propagation path (angular frequency = 2π * frequency, reflecting the speed of sound wave vibration). Based on the initial frequency spectrum, distance attenuation coefficient, and angular frequency, the first frequency spectrum after distance attenuation from the non-collision node to its next node is calculated using a formula. The sound energy value corresponding to each frequency in the first frequency spectrum will decrease due to distance attenuation.

[0093] Optionally, the formula for calculating the acoustic energy at a certain frequency in the first frequency spectrum in this embodiment of the invention is as follows: ,in, The acoustic energy (in units) corresponding to frequency f in the first frequency spectrum. D is the acoustic energy (in units) corresponding to frequency f in the initial frequency spectrum, D is the distance attenuation coefficient, and f is the current frequency (Hertz). The dominant frequency (Hertz) of the propagation path.

[0094] Continuing with the example of non-collision node A, its initial frequency spectrum shows acoustic energy of 100 units at 50 Hz, 80 units at 100 Hz, and 50 units at 200 Hz. The distance attenuation coefficient from node A to node B is 0.38, the dominant frequency corresponding to this propagation path is 100 Hz, and the angular frequency is 2π*100≈628 radians / second. The noise coupling simulation analysis device calculates the first frequency spectrum as follows: acoustic energy at 50 Hz is 100 units*0.38*(50 / 100)≈100*0.38*0.5≈19 units; acoustic energy at 100 Hz is 80 units*0.38*(100 / 100)≈30.4 units; and acoustic energy at 200 Hz is 50 units*0.38*(200 / 100)≈38 units. Therefore, the first frequency spectrum is 50 Hz corresponding to 19 units, 100 Hz corresponding to 30.4 units, and 200 Hz corresponding to 38 units.

[0095] Step 3044: Based on the first frequency spectrum of each sound ray before reaching the collision point, combined with the reflection attenuation coefficient and reflection phase angle of each sound ray at the collision point, determine the second frequency spectrum after reflection at each collision point. The reflection phase angle is determined based on the waveform interference between the reflection direction and the incident direction of each sound ray at the collision point.

[0096] Furthermore, for each collision point of each sound ray, the noise coupling simulation analysis device acquires the first frequency spectrum of the sound ray before it reaches that collision point.

[0097] Furthermore, the noise coupling simulation analysis device combines the reflection attenuation coefficient and reflection phase angle of the collision point (determined by the waveform interference between the reflection direction and the incident direction, reflecting the influence of the phase change of the reflected sound wave on the frequency spectrum) to calculate the second frequency spectrum of the collision point after reflection using a formula. The acoustic energy value corresponding to each frequency in the second frequency spectrum will change due to reflection attenuation and phase interference.

[0098] Optionally, the formula for calculating the acoustic energy of a certain frequency in the second frequency spectrum in this embodiment of the invention is as follows: ,in, The acoustic energy (in units) corresponding to frequency f in the second frequency spectrum. Let f be the acoustic energy (in units) corresponding to frequency f in the first frequency spectrum, and R be the reflection attenuation coefficient. f is the reflection phase angle (in radians) and f is the current frequency (in Hertz). Main frequency (Hertz).

[0099] Continuing with the above embodiment, the first frequency spectrum before the sound ray reaches collision point B is 19 units for 50 Hz, 30.4 units for 100 Hz, and 38 units for 200 Hz. The reflection attenuation coefficient at this collision point is 0.75, and the reflection phase angle is 60 degrees. The noise coupling simulation analysis device calculates the second frequency spectrum: the sound energy corresponding to 50 Hz is 19 units * 0.75 * cos(60° * (50 / 100)) ≈ 12.3 units; the sound energy corresponding to 100 Hz is 30.4 units * 0.75 * cos(60° * (100 / 100)) ≈ 11.4 units; and the sound energy corresponding to 200 Hz is 38 units * 0.75 * cos(60° * (200 / 100)) ≈ -14.25 units (the negative value indicates opposite phase, and the actual sound energy is taken as the absolute value of 14.25 units). Therefore, the second frequency spectrum corresponds to 12.3 units at 50 Hz, 11.4 units at 100 Hz, and 14.25 units at 200 Hz.

[0100] Step 3045: Associate and bind the first sound pressure level, second sound pressure level, first frequency spectrum, second frequency spectrum and coordinate mapping results of each node in each sound ray to obtain the acoustic parameter change sequence of each sound ray.

[0101] Furthermore, for each sound ray, the noise coupling simulation analysis device associates and binds the first sound pressure level, second sound pressure level, first frequency spectrum, and second frequency spectrum of each node with the coordinate mapping results (the spatial coordinates and sequence identifiers of each node). That is, the sequence identifier of each node corresponds to a sound pressure level (first sound pressure level or second sound pressure level), a frequency spectrum (first frequency spectrum or second frequency spectrum), and corresponding spatial coordinates, generating a sequence of acoustic parameter changes for each sound ray that includes the acoustic parameter values ​​of each node and the corresponding spatial coordinates.

[0102] Continuing with the above embodiments, the nodes of the sound ray include non-collision node A (order 1, coordinates (10, 5, 0)), collision point B (order 2, coordinates (18, 8, 2)), non-collision node C (order 3, coordinates (25, 0, 1)), etc. The parameters of each node are associated and bound: Sequence 1 corresponds to the initial sound pressure level of 85 dB, the initial frequency spectrum (100 units at 50 Hz, 80 units at 100 Hz, 50 units at 200 Hz), and the coordinates (10, 5, 0); Sequence 2 corresponds to the first sound pressure level of 28.5 dB, the first frequency spectrum (19 units at 50 Hz, 30.4 units at 100 Hz, 38 units at 200 Hz), and the coordinates (18, 8, 2), as well as the second sound pressure level of 14.5 dB, the second frequency spectrum (12.3 units at 50 Hz, 11.4 units at 100 Hz, 14.25 units at 200 Hz), and the coordinates (18, 8, 2); Sequence 3 corresponds to the first sound pressure level (assumed to be 10 dB) after propagation from the collision point B, the first frequency spectrum (assumed to be 5 units at 50 Hz, 4 units at 100 Hz, 6 units at 200 Hz), and the coordinates (25, 0, 1), etc., generating a sequence of acoustic parameter changes.

[0103] The embodiments of the present invention can accurately calculate the sound pressure level and frequency spectrum changes of each sound ray at each node along the propagation path, and accurately correlate them with the spatial coordinates of the nodes to generate a continuous sequence of acoustic parameter changes. Therefore, the acoustic parameter change sequence comprehensively reflects the evolution of acoustic characteristics of the sound ray from emission to propagation to the boundary, providing accurate basic data for subsequent noise coupling calculations of each spatial point in the three-dimensional spatial scene.

[0104] In one embodiment, steps 401 to 404 include: Step 401: Extract the node pair closest to the spatial point in the target acoustic parameter change sequence, and determine the relative position parameters of the spatial point in the path segment based on the distance between the two nodes in the node pair, the distance between the spatial point and the first node in the node pair, and the propagation path direction of the path segment where the node pair is located.

[0105] Optionally, for each spatial point in the three-dimensional scene, the noise coupling simulation analysis device first filters out all target sound rays whose propagation paths pass through that spatial point, and extracts the target acoustic parameter change sequence for each target sound ray. For each target sound ray, the noise coupling simulation analysis device finds the two nodes (i.e., node pairs) closest to that spatial point in its target acoustic parameter change sequence. These two nodes are located on both sides of the spatial point (or one side is the start / end point), forming a path segment that includes that spatial point.

[0106] Furthermore, the noise coupling simulation analysis device calculates the straight-line distance between the two nodes in the node pair, the straight-line distance between the spatial point and the first node in the node pair, and the propagation path direction vector of the path segment where the node pair is located. Based on the above parameters, it calculates and determines the relative position parameters of the spatial point in the path segment. The relative position parameters are used to characterize the position ratio of the spatial point between the two nodes.

[0107] Optionally, the formula for calculating the relative position parameter in this embodiment of the invention is as follows: ,in, For relative position parameters, This is the distance (in meters) between the spatial point and the first node in the node pair. This represents the distance (in meters) between the two nodes in a node pair.

[0108] In one embodiment, a spatial point P (coordinates (15, 6, 1)) is selected in a three-dimensional spatial scene. A target acoustic parameter variation sequence for a target sound ray contains node A (coordinates (10, 5, 0), sequence 1) and node B (coordinates (20, 10, 2), sequence 2). These two nodes form the node pair closest to spatial point P. A noise coupling simulation analysis device calculates the distance between node A and node B. The distance is 1.36 meters. This is the distance between spatial point P and node A. The distance is 5.20 meters. The propagation path direction vector of a node to its corresponding path segment is the coordinate vector of node B minus the coordinate vector of node A, i.e., (10, 5, 2). Relative position parameters. 5.20 / 11.36≈0.458 indicates that spatial point P is on the path segment from node A to node B, and its distance from node A is approximately 45.8%.

[0109] Step 402: Interpolate the sound pressure level and frequency spectrum of the node pair based on the relative position parameters to obtain the target sound pressure level and target frequency spectrum of each target sound ray at the spatial point.

[0110] Furthermore, for each target sound ray, the noise coupling simulation analysis device interpolates the sound pressure level and frequency spectrum of the node pair based on the relative position parameters. Optionally, in the interpolation process, this embodiment of the invention determines the values ​​of the sound pressure level and frequency spectrum at the spatial point between the parameters of the two nodes based on the relative position parameters, thereby obtaining the target sound pressure level and target frequency spectrum of each target sound ray at the spatial point. For the sound pressure level, linear interpolation or nonlinear interpolation based on distance attenuation is used; for the frequency spectrum, interpolation is performed separately for each frequency component to ensure the continuity of frequency characteristics.

[0111] Optionally, the target sound pressure level interpolation formula in this embodiment of the invention is: ,in, The target sound pressure level (decibels) at a spatial point. , These are the sound pressure levels (decibels) of the two nodes in the node pair, respectively. This refers to the relative position parameter.

[0112] Optionally, the target frequency spectrum interpolation formula in this embodiment of the invention is: ,in, The target frequency spectrum acoustic energy (in units) corresponding to frequency f at a spatial point; The acoustic energy (in units) of the first node A in the node pair at frequency f. The acoustic energy (in units) of the second node B in the node pair at frequency f. is the relative position parameter of the spatial point on the path segment (value range 0-1); f is the current calculated frequency (Hertz); The dominant frequency (Hertz) of the path segment to which the node is located reflects the main frequency characteristics of the sound wave in that path segment.

[0113] Continuing with the above embodiment, the sound pressure level at node A is 85 dB, corresponding to 100 units at 50 Hz, 80 units at 100 Hz, and 50 units at 200 Hz in the frequency spectrum; the sound pressure level at node B is 28.5 dB, corresponding to 19 units at 50 Hz, 30.4 units at 100 Hz, and 38 units at 200 Hz in the frequency spectrum. Relative position parameters. =0.458. The noise coupling simulation analysis device interpolates the sound pressure level: target sound pressure level. =85 + (28.5 - 85) * 0.458 ≈ 59.2 dB. Interpolating the frequency spectrum, the target frequency spectrum sound energy is 62.9 units at 50 Hz; 57.3 units at 100 Hz; and 44.5 units at 200 Hz. Therefore, the target sound pressure level at spatial point P is 59.2 dB, with target frequency spectra of 62.9 units at 50 Hz, 57.3 units at 100 Hz, and 44.5 units at 200 Hz.

[0114] Step 403: Based on the propagation path direction and the direction of the propagation path and the direction of the spatial point normal, determine the azimuth angle of each target sound ray at the spatial point.

[0115] Furthermore, for each target sound ray, the noise coupling simulation analysis device determines its propagation path direction at the spatial point (i.e., the propagation direction vector of the target sound ray when it passes through the spatial point). Further, the noise coupling simulation analysis device determines the normal direction of the spatial point (determined based on the virtual or actual surface where the spatial point is located; if it is a free space point, the normal direction can be set to be perpendicular to the propagation path plane).

[0116] Furthermore, the noise coupling simulation analysis device determines the azimuth angle of each target sound ray at the spatial point by calculating the angle between the propagation path direction vector and the spatial point normal direction vector. The azimuth angle is used to characterize the directional characteristics of the sound ray when it reaches the spatial point.

[0117] Optionally, the azimuth calculation formula in this embodiment of the invention is as follows: ,in, It is the azimuth angle (in radians). The propagation path direction is the unit vector. Let be the unit vector of the normal direction of a point in space.

[0118] Continuing with the above embodiment, the propagation path direction vector of the target sound ray at spatial point P is the direction vector of the node pair (10, 5, 2), which, after normalization, is (10 / 11.36, 5 / 11.36, 2 / 11.36) ≈ (0.880, 0.440, 0.176). Assuming spatial point P is a free space point, its normal direction vector is set to (0, 0, 1) (perpendicular to the horizontal plane). The noise coupling simulation analysis device calculates the angle (azimuth) between the propagation path direction vector and the normal direction vector. According to the vector dot product formula: cos =0.176, then ≈79.8 That is, the azimuth angle of the target sound ray at point P in space is approximately 79.8 degrees.

[0119] Step 404: Based on the azimuth angle, target sound pressure level and target frequency spectrum of each target sound ray at the spatial point, couple the noise coupling results at each spatial point.

[0120] Furthermore, the noise coupling simulation analysis device couples each target sound ray at a spatial point based on the azimuth angle, target sound pressure level, and target frequency spectrum to determine the noise coupling result at each spatial point, as described in steps 4041 to 4044.

[0121] The embodiments of the present invention can accurately calculate the noise coupling result of each spatial point in a three-dimensional spatial scene. Taking into account the path position, acoustic parameters and directional characteristics of the target sound ray, the final noise coupling result accurately reflects the actual noise superposition of each spatial point, laying a precise numerical foundation for the subsequent generation of the noise coupling distribution of the substation space.

[0122] In one embodiment, steps 4041 to 4044 include: Step 4041: Determine the acoustic parameter coupling coefficient based on the azimuth angle of each target sound ray at the spatial point. The acoustic parameter coupling coefficient reflects the influence of the sound ray direction on the acoustic parameters received at the spatial point.

[0123] Optionally, for each target sound ray, the noise coupling simulation analysis device, based on the azimuth angle of the sound ray at a spatial point, uses a functional relationship... ,in, Let be the acoustic parameter coupling coefficient of the i-th target sound ray. Calculate the acoustic parameter coupling coefficient based on the azimuth angle (in radians) of the sound ray at the spatial point. The acoustic parameter coupling coefficient ranges from 0 to 1, and its magnitude is related to the azimuth angle. The smaller the azimuth angle (the closer the sound ray propagation direction is to the normal direction of the spatial point), the larger the coupling coefficient, indicating that the sound ray direction has a greater influence on the acoustic parameters received at the spatial point; conversely, the larger the azimuth angle, the smaller the coupling coefficient, and the smaller the influence. This is used to quantify the degree of influence of the sound ray direction on the acoustic parameters received at the spatial point.

[0124] Continuing with the above embodiment, spatial point P has two target sound rays. The azimuth angle of the first target sound ray is 79.8 degrees, and the azimuth angle of the second target sound ray is 30 degrees. The noise coupling simulation analysis device uses the azimuth angle combined with the functional relationship... Calculate the coupling coefficient of the acoustic parameters, specifically the coupling coefficient of the first acoustic ray. ≈0.176; Coupling coefficient of the second vocal ray ≈0.866. This indicates that because the second sound ray has a smaller azimuth angle, its direction has a greater impact on the acoustic parameters received by the spatial point P.

[0125] Step 4042: Modulate the target sound pressure level and target frequency spectrum based on the coupling coefficient of the acoustic parameters of each target sound ray to obtain the coupled sound pressure level and coupled frequency spectrum of each target sound ray after directional coupling.

[0126] Furthermore, for each target sound ray, the noise coupling simulation analysis device multiplies the acoustic parameter coupling coefficient with the target sound pressure level and target frequency spectrum of the sound ray at the spatial point, respectively, to complete the modulation of the target sound pressure level and target frequency spectrum, and obtain the coupled sound pressure level and coupled frequency spectrum of each target sound ray after directional coupling, which reflects the actual influence of the sound ray direction on the acoustic parameters received at the spatial point.

[0127] Continuing with the above embodiment, the target sound pressure level of the first target sound ray is 59.2 dB, the target frequency spectrum is 62.9 units at 50 Hz, 57.3 units at 100 Hz, and 44.5 units at 200 Hz, and the acoustic parameter coupling coefficient is 0.176. The noise coupling simulation analysis device calculates its coupled sound pressure level as 59.2 dB * 0.176 ≈ 10.42 dB; the coupled frequency spectrum is: 50 Hz 62.9 units * 0.176 ≈ 11.07 units, 100 Hz 57.3 units * 0.176 ≈ 10.08 units, and 200 Hz 44.5 units * 0.176 ≈ 7.83 units. The target sound pressure level of the second target sound ray is 65 dB, the target frequency spectrum is 50 units at 50 Hz, 70 units at 100 Hz, and 30 units at 200 Hz, and the acoustic parameter coupling coefficient is 0.866. The coupled sound pressure level is approximately 56.29 dB (65 dB * 0.866 ≈ 55 dB); the coupled frequency spectrum is approximately 43.3 units (50 Hz * 0.866 ≈ 40.62 Hz (100 Hz * 0.866 ≈ 60.62 Hz (70 Hz * 0.866 ≈ 60.62 Hz (30 Hz * 0.866 ≈ 25.98 ...

[0128] Step 4043: Based on the phase difference of each target sound ray in the frequency components, the coupled frequency spectra of each target sound ray are interferometrically superimposed to obtain the superimposed total frequency spectrum.

[0129] Furthermore, for all target sound rays at a spatial point, the noise coupling simulation analysis device determines the phase difference of each sound ray at each frequency component (the phase difference is determined by factors such as the difference in the propagation path length of the sound ray and the number of reflections, and can be calculated by the ratio of the path length to the wavelength).

[0130] Furthermore, the noise coupling simulation analysis device performs interference superposition on the acoustic energy of the corresponding frequency components in the coupled frequency spectrum of each sound ray based on the phase difference. Sound energies in the same phase are added together, and sound energies in opposite phases are subtracted together, finally obtaining the superimposed total frequency spectrum. This spectrum reflects the interference effect of different sound rays on each frequency component.

[0131] Optionally, the formula for the total frequency spectrum acoustic energy after superposition in this embodiment of the invention is as follows: ,in, The total frequency spectrum acoustic energy (in units) after superposition at frequency f. =2 The phase difference of the i-th target sound ray at frequency f ( For sound path difference, (where f is the wavelength, j is the imaginary unit, and n is the number of target sound rays).

[0132] Continuing with the above embodiment, the frequency spectra of the two target sound rays after coupling are as follows: the first spectrum is 11.07 units at 50 Hz, 10.08 units at 100 Hz, and 7.83 units at 200 Hz; the second spectrum is 43.3 units at 50 Hz, 60.62 units at 100 Hz, and 25.98 units at 200 Hz. Calculations show that the phase difference between the 50 Hz frequency components is 0 degrees (in phase), and the resulting acoustic energy after interference superposition is 11.07 units + 43.3 units = 54.37 units; the phase difference between the 100 Hz frequency components is 0 degrees, and the resulting acoustic energy after superposition is 10.08 units + 60.62 units = 70.7 units; the phase difference between the 200 Hz frequency components is 180 degrees (out of phase), and the resulting acoustic energy after superposition is |7.83 units - 25.98 units| = 18.15 units. Therefore, the total frequency spectrum after superposition is 54.37 units at 50 Hz, 70.7 units at 100 Hz, and 18.15 units at 200 Hz.

[0133] Step 4044: Integrate the superimposed total frequency spectrum and the sound pressure level after coupling of each target sound ray to obtain the noise coupling result of each spatial point.

[0134] Furthermore, the noise coupling simulation analysis device integrates the superimposed total frequency spectrum and the sound pressure level after coupling of each target sound ray to obtain the noise coupling result for each spatial point.

[0135] The embodiments of the present invention can comprehensively consider the influence of the direction of the sound ray on the acoustic parameters received at the spatial point, and combine the interference effect between the sound rays to finally obtain the accurate noise coupling result for each spatial point. Therefore, the total sound pressure level and total frequency spectrum contained in the noise coupling result fully reflect the noise superposition characteristics at the spatial point, providing accurate single-node data support for the subsequent generation of noise coupling distribution in the substation space.

[0136] The substation noise coupling simulation analysis device based on sound ray simulation scenario provided by the present invention will be described below. The substation noise coupling simulation analysis device based on sound ray simulation scenario described below can be referred to in correspondence with the substation noise coupling simulation analysis method based on sound ray simulation scenario described above.

[0137] Optional, refer to Figure 2 , Figure 2 This is a schematic diagram of the substation noise coupling simulation analysis device based on a sound ray simulation scenario provided by the present invention. The substation noise coupling simulation analysis device based on a sound ray simulation scenario includes: The spatial scene construction module 210 is used to construct a three-dimensional spatial scene based on the spatial location and geometric size information of spatial entities within the substation. The sound path tracking module 220 is used to control the noise source to emit multiple sound rays according to a preset angle for each noise source in the substation, and to determine the collision position of each sound ray with the spatial entity when it propagates in the three-dimensional space scene and the reflection direction after each collision, until the sound ray propagates to the boundary of the three-dimensional space scene; each sound ray carries the initial acoustic parameters of the noise source. The acoustic parameter monitoring module 230 is used to determine the changes in acoustic parameters of each sound ray during propagation due to distance attenuation and collision reflection based on the propagation path, collision position and reflection direction of each sound ray, and generate an acoustic parameter change sequence for each sound ray; the acoustic parameter change sequence includes the acoustic parameter values ​​and corresponding spatial coordinates of each node on the propagation path of the sound ray; The single spatial point noise coupling module 240 is used to couple each spatial point in the three-dimensional spatial scene based on the sequence of changes in the target acoustic parameters of the target sound ray passing through the spatial point through the propagation path, and to determine the noise coupling result of each spatial point. The full-space point noise coupling module 250 is used to traverse the noise coupling results of each spatial point in the three-dimensional spatial scene to obtain the noise coupling distribution of the space where the substation is located.

[0138] This invention, by emitting multiple sound rays and tracking their propagation paths, can comprehensively record the reflection and attenuation of noise during propagation, covering all possible propagation paths within the substation, no longer limited to fixed point locations. Furthermore, by tracking the acoustic parameter change sequence of each sound ray carrying the initial acoustic parameters of the corresponding noise source as it propagates at various nodes, the noise coupling result at each spatial point is determined by coupling based on the acoustic parameter change sequence of all propagation paths passing through that spatial point, accurately reflecting the coupling result of different noise sources at that spatial point. Moreover, when the equipment's operating state changes dynamically, the acoustic parameters of the noise source can be updated, sound rays can be re-emitted, and noise coupling can be performed to quickly obtain a new noise coupling distribution, thus enabling dynamic response to changes in equipment state. Therefore, this invention, through full-path sound ray tracking, acoustic parameter change recording, and full-space-point coupling, can accurately reflect the dynamic propagation and coupling characteristics of noise in complex spatial scenarios.

[0139] Please see Figure 3 , Figure 3 An embodiment diagram of an electronic device provided in accordance with the present invention. For example... Figure 3 As shown, an embodiment of the present invention provides an electronic device 300, including a memory 310, a processor 320, and a computer program 311 stored in the memory 310 and executable on the processor 320. When the processor 320 executes the computer program 311, it implements the processes of steps 10 to 50.

[0140] Please see Figure 4 , Figure 4 An embodiment diagram of a computer-readable storage medium provided in accordance with an embodiment of the present invention is shown. Figure 4 As shown, this embodiment provides a computer-readable storage medium 400 on which a computer program 311 is stored. When the computer program 311 is executed by a processor, it implements the processes of steps 10 to 50.

[0141] On the other hand, the present invention also provides a computer program product, which includes a computer program that can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer is able to execute the substation noise coupling simulation analysis method based on the sound ray simulation scenario provided by the above methods, which includes the process of steps 10 to 50.

[0142] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs. Those skilled in the art can understand and implement this without any creative effort.

[0143] Through the above description of the embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus necessary general-purpose hardware platforms, and of course, it can also be implemented by hardware. Based on this understanding, the above technical solutions, in essence or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute the methods described in the various embodiments or some parts of the embodiments.

[0144] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for noise coupling simulation analysis of substations based on acoustic ray simulation scenarios, characterized in that, include: A three-dimensional spatial scene is constructed based on the spatial location and geometric dimensions of spatial entities within the substation. For each noise source in the substation, multiple sound rays are emitted by the noise source according to a preset angle, and the collision position of each sound ray with the spatial entity when it propagates in the three-dimensional space scene and the reflection direction after each collision are determined until the sound ray propagates to the boundary of the three-dimensional space scene; each sound ray carries the initial acoustic parameters of the noise source. Based on the propagation path, collision location, and reflection direction of each sound ray, the changes in acoustic parameters of each sound ray during propagation due to distance attenuation and collision reflection are determined, generating an acoustic parameter change sequence for each sound ray; the acoustic parameter change sequence includes the acoustic parameter values ​​and corresponding spatial coordinates of each node on the propagation path of the sound ray. For each spatial point in the three-dimensional spatial scene, the noise coupling result of each spatial point is determined by coupling based on the sequence of changes in the target acoustic parameters of the target sound ray passing through the spatial point through the propagation path. By traversing the noise coupling results of each spatial point in the three-dimensional spatial scene, the noise coupling distribution of the space where the substation is located is obtained.

2. The substation noise coupling simulation analysis method based on acoustic ray simulation scenario according to claim 1, characterized in that, Acoustic parameters include sound pressure level and frequency spectrum; Based on the propagation path, collision location, and reflection direction of each sound ray, the acoustic parameter changes of each sound ray during propagation due to distance attenuation and collision reflection are determined, generating a sequence of acoustic parameter changes for each sound ray, including: For each sound ray, based on the collision position and reflection direction of the sound ray, each node on the propagation path is extracted, and the three-dimensional spatial coordinates of each node and its ordinal position in the propagation path are mapped and bound to obtain the coordinate mapping result of each node; the nodes include collision points and non-collision points. For each collision point, the reflection attenuation coefficient of each sound ray at the collision point is determined based on the sound energy absorption coefficient and frequency response coefficient of the material of the collision interface through which the collision point passes, combined with the dominant frequency and incident angle of each sound ray at the collision point; the dominant frequency is determined based on the frequency spectrum of each sound ray at the collision point, and the incident angle is determined based on the reflection direction and incident direction of each sound ray at the collision point. For each non-collision node, the distance attenuation coefficient from the non-collision node to its next node is determined based on the straight-line distance between the non-collision node and its next node, combined with the dominant frequency wavelength corresponding to each sound ray in each propagation path; the dominant frequency wavelength is determined based on the frequency spectrum of each sound ray at the non-collision node. Based on the analysis of the reflection attenuation coefficient, distance attenuation coefficient, initial acoustic parameters of each node, and coordinate mapping results, an acoustic parameter change sequence for each sound ray is generated.

3. The substation noise coupling simulation analysis method based on acoustic ray simulation scenario according to claim 2, characterized in that, The analysis, based on the reflection attenuation coefficient, distance attenuation coefficient, and the initial acoustic parameters and coordinate mapping results of each node, generates a sequence of acoustic parameter changes for each sound ray, including: Based on the initial sound pressure level and distance attenuation coefficient of each sound ray at a non-collision node, and combined with the angle between each sound ray's propagation path and the reference axis in the three-dimensional spatial scene, the first sound pressure level attenuated by distance from the non-collision node to its next node is determined. The second sound pressure level after reflection at each collision point is determined based on the first sound pressure level before each sound ray reaches the collision point, combined with the reflection attenuation coefficient and reflection angle of each sound ray at the collision point; the reflection angle is determined based on the reflection direction of each sound ray at the collision point. Based on the initial frequency spectrum and distance attenuation coefficient of each sound ray at the non-collision node, and combined with the angular frequency of each sound ray in each propagation path, the first frequency spectrum of each sound ray after distance attenuation from the non-collision node to its next node is determined. The second frequency spectrum after reflection at each collision point is determined based on the first frequency spectrum of each sound ray before it reaches the collision point, combined with the reflection attenuation coefficient and reflection phase angle of each sound ray at the collision point; the reflection phase angle is determined based on the waveform interference between the reflection direction and the incident direction of each sound ray at the collision point. By associating and binding the first sound pressure level, second sound pressure level, first frequency spectrum, second frequency spectrum, and coordinate mapping results of each node in each sound ray, a sequence of acoustic parameter changes for each sound ray is obtained.

4. The substation noise coupling simulation analysis method based on acoustic ray simulation scenario according to claim 2, characterized in that, The method of extracting each node along the propagation path based on the collision position and reflection direction of the sound ray includes: For each sound ray, based on the initial emission direction, the first entity boundary of the spatial entity, and the starting coordinates of the sound ray, the initial spatial domain that intersects with the first entity boundary during the initial propagation stage is determined. Based on the collision position of each first entity boundary and the sound ray within the initial spatial domain, the theoretical collision distance between each first entity boundary and the sound ray is determined, and based on the theoretical collision distance, concavity and convexity, and curvature of each first entity boundary, the collision priority of each first entity boundary is determined. Based on the initial launch direction, the starting point coordinates, and the first boundary normal vector and first position parameter of the first target entity with the highest collision priority, the first spatial coordinates of the first collision point are determined. Based on the first boundary normal vector, the first spatial coordinates, and the reflection direction, the nodes of the sound ray on the propagation path are determined.

5. The substation noise coupling simulation analysis method based on acoustic ray simulation scenario according to claim 4, characterized in that, The step of determining each node of the sound ray on the propagation path based on the first boundary normal vector, the first spatial coordinates, and the reflection direction includes: Based on the first boundary normal vector and the reflection direction, the reflection is reconstructed to obtain the reconstructed reflection direction at the first collision point; Based on the reconstructed reflection direction, the second entity boundary of the spatial entity that did not collide, and the first spatial coordinates, a secondary spatial domain is determined where the reflected sound ray intersects with the second entity boundary. Based on the reconstructed reflection direction, the first spatial coordinates, and the second boundary normal vector and second position parameters of the second target entity boundary with the highest collision priority in the secondary spatial domain, the second spatial coordinates of the secondary collision point are determined. Based on the spatial coordinates of the start point, end point, and non-collision point in the sound ray, combined with the first and second spatial coordinates, the start point, end point, first collision point, second collision point, and non-collision point are sorted according to the spatial coordinate order to obtain the nodes of the sound ray on the propagation path; the non-collision point is the node between the start point and the end point excluding the first and second collision points.

6. The substation noise coupling simulation analysis method based on acoustic ray simulation scenario according to claim 1, characterized in that, The coupling based on the target acoustic parameter change sequence of the target sound ray passing through spatial points along the propagation path determines the noise coupling result at each spatial point, including: Extract the node pair closest to the spatial point from the target acoustic parameter change sequence, and determine the relative position parameters of the spatial point in the path segment based on the distance between the two nodes in the node pair, the distance between the spatial point and the first node in the node pair, and the propagation path direction of the path segment where the node pair is located. Interpolating the sound pressure level and frequency spectrum of the node pair based on the relative position parameters yields the target sound pressure level and target frequency spectrum of each target sound ray at the spatial point. Based on the direction of the propagation path and the direction of the propagation path relative to the normal direction of the spatial point, the azimuth angle of each target sound ray at the spatial point is determined. The noise coupling result at each spatial point is determined by coupling based on the azimuth angle, sound pressure level, and frequency spectrum of each target sound ray at the spatial point.

7. The substation noise coupling simulation analysis method based on acoustic ray simulation scenario according to claim 6, characterized in that, The coupling based on the azimuth angle, sound pressure level, and frequency spectrum of each target sound ray at a spatial point determines the noise coupling result at each spatial point, including: The acoustic parameter coupling coefficient is determined based on the azimuth angle of each target sound ray at a spatial point; the acoustic parameter coupling coefficient reflects the influence of the sound ray direction on the acoustic parameters received at the spatial point; Based on the coupling coefficient of the acoustic parameters of each target sound ray, the target sound pressure level and target frequency spectrum are modulated to obtain the coupled sound pressure level and coupled frequency spectrum of each target sound ray after directional coupling; Based on the phase difference of each target sound ray in the frequency components, the frequency spectra of each target sound ray after coupling are interferometrically superimposed to obtain the superimposed total frequency spectrum; The total frequency spectrum after superposition and the sound pressure level after coupling of each target sound ray are integrated to obtain the noise coupling result of each spatial point.

8. A substation noise coupling simulation analysis device based on acoustic ray simulation scenarios, characterized in that, The method for substation noise coupling simulation analysis based on acoustic ray simulation scenarios as described in any one of claims 1 to 7; the substation noise coupling simulation analysis device based on acoustic ray simulation scenarios includes: The spatial scene construction module is used to construct a three-dimensional spatial scene based on the spatial location and geometric dimensions of spatial entities within the substation. The sound ray path tracking module is used to control each noise source in the substation to emit multiple sound rays at a preset angle, and to determine the collision position of each sound ray with the spatial entity when it propagates in the three-dimensional space scene and the reflection direction after each collision, until the sound ray propagates to the boundary of the three-dimensional space scene; each sound ray carries the initial acoustic parameters of the noise source. The acoustic parameter monitoring module is used to determine the changes in acoustic parameters of each sound ray during propagation due to distance attenuation and collision reflection, based on the propagation path, collision location, and reflection direction of each sound ray, and to generate an acoustic parameter change sequence for each sound ray; the acoustic parameter change sequence includes the acoustic parameter values ​​and corresponding spatial coordinates of each node on the propagation path of the sound ray; A single spatial point noise coupling module is used to couple each spatial point in the three-dimensional spatial scene based on the sequence of changes in target acoustic parameters of the target sound ray passing through the spatial point through the propagation path, and to determine the noise coupling result of each spatial point. The full-space point noise coupling module is used to traverse the noise coupling results of each spatial point in the three-dimensional spatial scene to obtain the noise coupling distribution of the space where the substation is located.

9. An electronic device, comprising: Memory, used to store computer software programs; A processor for reading and executing the computer software program, characterized in that, when the processor executes the computer software program, it implements the substation noise coupling simulation analysis method based on sound ray simulation scenario as described in any one of claims 1 to 7.

10. A non-transitory computer-readable storage medium, wherein a computer software program is stored therein, characterized in that, When the computer software program is executed by the processor, it implements the substation noise coupling simulation analysis method based on sound ray simulation scenario as described in any one of claims 1 to 7.

Images