CN113297759A - GIL telescopic joint acoustic transmission modeling method - Google Patents
GIL telescopic joint acoustic transmission modeling method Download PDFInfo
- Publication number
- CN113297759A CN113297759A CN202110507493.8A CN202110507493A CN113297759A CN 113297759 A CN113297759 A CN 113297759A CN 202110507493 A CN202110507493 A CN 202110507493A CN 113297759 A CN113297759 A CN 113297759A
- Authority
- CN
- China
- Prior art keywords
- gil
- finite element
- wave
- attenuation
- simulation model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
The invention relates to a sound wave transmission modeling method, belongs to the technical field of power monitoring, and particularly relates to a GIL telescopic joint sound wave transmission modeling method. The method combines the test and numerical simulation modes, can consider the noise interference when the GIL ultrasonic fault positioning online monitoring system works and the influence of the sensitivity of the sensor on the result, calculates the attenuation and time delay of sound waves passing through different types of GIL telescopic joints, and has strong universality and wider engineering application prospect.
Description
Technical Field
The invention relates to a sound wave transmission modeling method, belongs to the technical field of power monitoring, and particularly relates to a GIL telescopic joint sound wave transmission modeling method.
Background
The GIL is a power transmission line with large transmission capacity, high reliability and small environmental interaction, the overall length of the GIL is generally long, and the GIL is usually installed in a segmented modular manner. When an insulation fault occurs, how to quickly and accurately locate the fault position so as to timely replace the damaged pipe section, and other pipe sections do not need to be disassembled, so that the hot topic in GIL operation and maintenance is achieved.
When a telescopic joint structure exists between a sensor and a fault point, the existing GIL ultrasonic fault location online monitoring method has the problems of low location result accuracy, poor system reliability and the like, because a large number of gas-solid interfaces exist in the telescopic joint, sound waves transmitted along a GIL shell generate complex refraction and reflection and waveform conversion phenomena in the telescopic joint, and relatively serious time delay and attenuation are generated relative to a straight pipe section. Therefore, the propagation characteristic of the sound wave passing through the GIL telescopic joint needs to be mastered, so that a basis is provided for the improvement of the GIL ultrasonic fault positioning online monitoring method.
The propagation condition of sound waves passing through the telescopic link is complex, analysis and calculation are difficult to perform in a theoretical mode, and no relevant empirical formula exists. The numerical simulation method can solve the problems of large theoretical calculation difficulty, complicated test and the like faced by sound propagation in a complex structure. However, the numerical simulation results under ideal conditions, the actual GIL ultrasonic fault location online monitoring system is often interfered by field noise when working, and the sensitivity of the used sensors is different. When insulation faults occur, effective wave heads of time domain waveforms which can be measured by ultrasonic fault positioning online monitoring systems installed on different GILs are different, and action thresholds set by the systems are different.
Therefore, it is necessary to consider the influence of the above factors in the results obtained by numerical simulation so that the numerical simulation results better guide the practice. On the other hand, GIL expansion joints are of various types, and the number of corrugated pipes in expansion joints used in different GILs is different. Therefore, the invention provides a calculation method for reducing attenuation and delay of sound waves passing through an expansion joint according to the joint under the condition of considering the sensitivity of an actual ultrasonic positioning online monitoring system and the interference of external noise by combining tests and numerical simulation.
Disclosure of Invention
The following presents a simplified summary of one or more aspects in order to provide a basic understanding of such aspects. This summary is not an extensive overview of all contemplated aspects, and is intended to neither identify key or critical elements of all aspects nor delineate the scope of any or all aspects. Its sole purpose is to present some concepts of one or more aspects in a simplified form as a prelude to the more detailed description that is presented later.
The invention mainly aims to solve the technical problems in the prior art and provides a GIL telescopic joint acoustic wave transmission modeling method. The method combines the test and numerical simulation modes, can consider the noise interference when the GIL ultrasonic fault positioning online monitoring system works and the influence of the sensitivity of the sensor on the result, calculates the attenuation and time delay of sound waves passing through different types of GIL telescopic joints, and has strong universality and wider engineering application prospect.
In order to solve the problems, the scheme of the invention is as follows:
a GIL telescopic joint acoustic wave transmission modeling method comprises the following steps:
the method comprises the following steps of (1) testing acoustic transmission characteristics, namely knocking a shell of a GIL transmission system comprising an expansion joint, and recording a knocking force waveform signal and an actually measured vibration time domain waveform signal transmitted by sound waves in the GIL;
establishing a first finite element simulation model and a second finite element simulation model representing the geometric form of the GIL transmission system, wherein elements and the geometric dimensions in the first finite element simulation model are the same as those of the GIL transmission system, and the second finite element simulation model is obtained by replacing an expansion joint corrugated pipe in the first finite element simulation model with a GIL straight pipe with the same length as that of the expanded corrugated pipe;
a finite element numerical simulation step, fitting according to the actually measured knocking force waveform signal to obtain an excitation force waveform; applying the excitation force waveform to the first finite element simulation model and the second finite element simulation model in a total force mode, and recording simulated vibration time domain waveform signals obtained by simulation in each finite element model;
analyzing an effective wave head of the simulated vibration time domain waveform signal and a time delay calculation reference point on the wave head, and calculating the attenuation and time delay characteristics of the sound wave after passing through the first finite element simulation model and the second finite element simulation model;
and acoustic transmission modeling, namely establishing a mathematical model of the attenuation and delay of the GIL expansion joint reduced by sections based on the attenuation and delay characteristics.
Preferably, in the sound transmission characteristic test step, a standard knocking sound source is generated by knocking the GIL housing by a piezoelectric impact hammer, and two sensors mounted in front of and behind the GIL telescopic joint unit are selected to record vibration time domain waveform signals transmitted by sound waves in the GIL.
Preferably, in the above method for modeling GIL telescopic joint acoustic wave transmission, the finite element model modeling step adds an acoustic-structure coupling boundary condition, a support constraint boundary condition, a medium damping characteristic, and infinite length boundary conditions at both ends of the model to the finite element simulation model.
Preferably, in the above method for modeling GIL telescopic joint acoustic wave transmission, in the finite element numerical simulation step, the excitation position at which the excitation force waveform is applied to the GIL finite element simulation model and the measurement position at which the simulated vibration time domain waveform signal is recorded correspond to the tapping position and the recording position in the acoustic transmission characteristic test step.
Preferably, in the finite element numerical simulation step, the excitation force waveform is obtained by fitting the actually measured knocking force waveform signal in a manner of splicing a quadratic function and a gaussian pulse function in sections.
Preferably, in the attenuation and delay analysis step, in a wave group of the horizontal-vertical aliasing wave, a region between a start point of a rising part of each wave head in a wave group of a wave head of the waveform and a peak point of the wave head is selected as a distribution region of a delay calculation reference point, a wave velocity range corresponding to each region is calculated, an effective wave head in the wave group of the waveform head of the simulated vibration time domain waveform signal is determined based on an actual measurement wave velocity obtained by the actual measurement vibration time domain waveform signal calculation so as to determine a position of the delay calculation reference point, and an attenuation amount and a delay amount of the sound wave after passing through the first finite element simulation model and the second finite element simulation model are calculated based on the position of the reference point.
Preferably, in the GIL telescopic joint sound wave transmission modeling method, based on the attenuation a1 and the delay T1 of sound waves passing through the first finite element simulation model and the attenuation a2 and the delay T2 before and after passing through the second finite element simulation model, the additional attenuation a1-a2 and the delay T1-T2 caused by the telescopic joint structure are calculated; calculating the joint attenuation amount and the time delay amount of a single telescopic joint based on the extra attenuation amount A1-A2 and the time delay amount T1-T2 caused by the telescopic joint structure; and establishing a mathematical model of the attenuation and the time delay of the GIL expansion joint reduced according to the sections based on the section attenuation and the section time delay.
Therefore, compared with the prior art, the invention has the following advantages:
(1) an excitation sound source is generated by knocking the GIL shell through a force hammer, so that an excitation signal can be recorded, a large number of tests can be rapidly carried out, and damage to the GIL due to the fact that the sound source is punctured through an actual alternating current withstand voltage test can be effectively avoided;
(2) the sensor of the ultrasonic positioning on-line monitoring system arranged on the GIL shell is used for testing, so that the test arrangement is convenient, the test cost is reduced, and the test result has stronger pertinence;
(3) the influence of noise interference and sensor sensitivity on the result in an actual system is considered, and an effective wave head of a numerical simulation result waveform and a time delay calculation reference point on the wave head are determined according to a test result, so that the numerical simulation result can be suitable for different GIL ultrasonic positioning online monitoring systems;
(4) the algorithm is strong in universality, and can accurately and quantitatively calculate the attenuation and delay of sound waves passing through various GIL telescopic joints with different structures.
Drawings
The accompanying drawings, which are incorporated herein and form a part of the specification, illustrate embodiments of the present invention and, together with the description, further serve to explain the principles of the invention and to enable a person skilled in the pertinent art to make and use the disclosure.
FIG. 1 is a flowchart of a method for calculating attenuation and delay of sound waves passing through a GIL telescopic joint according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a GIL telescopic joint unit acoustic transmission characteristic test arrangement provided in an embodiment of the present invention;
fig. 3 is a schematic diagram of a process of straightening a corrugated pipe in a GIL telescopic joint into a straight pipe according to an embodiment of the present invention, where the length of the telescopic joint is L0, and the length of the straightened straight pipe is L;
fig. 4 is an impact force waveform output by the piezoelectric impact hammer according to the embodiment of the present invention and a waveform obtained by segment-wise splicing a quadratic function and a gaussian pulse function;
fig. 5 is a schematic diagram of an acoustic vibration time domain waveform on a GIL acoustic-structural coupling finite element simulation model surface probe and a method for selecting a reference point distribution region for calculating a head wave group delay of the waveform according to an embodiment of the present invention, where (a) is an acoustic vibration time domain waveform on a finite element simulation model probe, and (b) is a selection of a reference point distribution region for calculating a head wave group delay of an acoustic vibration time domain waveform.
Embodiments of the present invention will be described with reference to the accompanying drawings.
Detailed Description
Examples
The invention is further elucidated with reference to the drawings and the embodiments.
As shown in FIG. 1, the invention provides a GIL telescopic joint acoustic transmission modeling method, which comprises the following steps:
step 1: as shown in fig. 2, the GIL expansion joint unit sound transmission characteristic test is carried out on a section of GIL containing expansion joints, and comprises the following steps:
step 1-1: the method comprises the steps that a standard knocking sound source is generated in a mode that a piezoelectric type impact hammer knocks a GIL shell, and the fact that the GIL shell is made of aluminum alloy and is soft in texture is considered, so that the impact hammer uses a rubber hammer head and is connected to an oscilloscope through a signal conditioner to output an excitation force signal;
step 1-2: a sensor of an ultrasonic positioning online monitoring system arranged on a GIL shell outputs a sound wave vibration time domain waveform signal, two sensors arranged in front of and behind a GIL expansion joint unit are selected for testing, and the sensors are fixed on the GIL shell through a steel rolled strip.
The purpose of the test is to obtain a vibration time domain waveform actually-measured signal acquired by the GIL ultrasonic positioning online monitoring system, the signal contains the influence of system sensitivity and background noise on a measurement result, the signal is compared with a vibration time domain waveform simulation signal obtained by a subsequent simulation model to determine a calculation reference point in a simulation waveform, so that the simulation result is more practical, and meanwhile, optimization can be performed on different GIL ultrasonic positioning online monitoring systems.
Step 2: establishing an acoustic-structure coupling finite element simulation model containing a GIL expansion joint and a straight pipe unit, comprising the following steps of:
step 2-1: establishing a corresponding geometric model according to the sizes of structural units such as a GIL expansion joint, a flange, a basin-type insulator, a three-pillar insulator, a movable support, a fixed support and the like, and establishing a section of GIL three-dimensional simulation model containing the expansion joint in a module combination mode, wherein the integral structure of the model is consistent with the GIL section tested in the sound transmission characteristic test of the GIL expansion joint unit, and the model is called as a first finite element simulation model;
step 2-2: calculating the length L of the whole unit after the corrugated pipe is straightened in the expansion joint unit, wherein the straightening process is shown in figure 3, the expansion joint unit in the first finite element simulation model established in the step 2-1 is replaced by a section of GIL straight pipe unit with the length of L, the rest part of the model is kept unchanged, and the newly established model is called a second finite element simulation model;
step 2-3: adding corresponding boundary conditions including an acoustic-structure coupling boundary condition, a support constraint boundary condition, a medium damping characteristic, infinite length boundary conditions at two ends of the model and the like to the established first and second finite element simulation models;
the acoustic-structure coupling boundary condition means that at the acoustic-solid interface, fluid pressure acts on the solid with normal unit area load; acceleration of the solid along the normal to the interface acts on the fluid, corresponding to normal acceleration, to create an acoustic source. The support constraint boundary condition means that a fixed constraint is set at the joint of the support and the foundation, namely the displacement of the joint is considered to be 0; the medium damping characteristic boundary condition is that Rayleigh damping is added to a solid domain in a finite element simulation model, and a damping coefficient is 0.0005;
the infinite boundary condition refers to the condition that perfect matching layers, low reflection and spherical wave radiation boundary are added at two ends of a finite element simulation model, so that the boundary at the two ends of the model can completely absorb the transmitted sound waves without reflecting the sound waves.
Step 2-4: and (3) carrying out proper mesh subdivision on the model according to the frequency of the excitation sound source, wherein the number of mesh units per unit wavelength is generally 5. Because the transmission of sound waves in the GIL is a transient process changing with time, a transient (time-dependent) solver is required to solve, and an appropriate solution time step is set for the transient process, generally a fixed time step is used, and the step is 1/60 of the excitation sound source period.
And step 3: carrying out finite element numerical simulation based on the established first and second finite element simulation models to couple two physical fields of pressure acoustics and solid mechanics to obtain a numerical simulation result of the acoustic vibration signal under the excitation of the standard knocking sound source on the GIL shell, and comprising the following steps of:
step 3-1: according to an impact force waveform signal output when a piezoelectric impact hammer strikes the GIL shell, an excitation force waveform is fitted in a mode of splicing a quadratic function and a Gaussian pulse function in a segmented mode, an excitation sound source generated when the GIL shell is struck is simulated, a fitting result is shown in FIG. 4, and a mathematical model of a fitting excitation force waveform function f (t) is as follows, wherein A is a fitting function peak value, t is a time independent variable, and t isaIs the pulse start time, tpAt the time of the peak of the pulse, tbAt the end of the pulse, f0Gaussian pulse bandwidth, e is the natural constant:
applying the excitation force waveform to the surface of the GIL finite element simulation model in a total force mode, keeping the applied position consistent with the knocking point of the impact force hammer in the test of the step 1, and calculating the sound vibration signal of the GIL surface under the excitation of the sound source;
step 3-2: and (3) according to the installation position of the sensor on the GIL in the test in the step 1, arranging probes at corresponding positions on the surface of the first finite element simulation model, and recording the acoustic vibration time-domain waveform on each probe.
And 4, step 4: determining an effective wave head of a numerical simulation result waveform and a time delay calculation reference point on the wave head under the condition of considering the sensitivity of an actual ultrasonic positioning online monitoring system and external noise interference according to the test result obtained in the step 1, and analyzing the attenuation and time delay characteristics of a sound wave passing through an expansion joint and a straight pipe unit with the same length as the expansion joint unit corrugated pipe after being expanded, wherein the method comprises the following steps:
step 4-1: the acoustic vibration time domain waveform recorded by the surface probe of the first finite element simulation model is shown in fig. 5, because the propagation speed of longitudinal waves in a solid is higher than that of transverse waves, and the excitation mode mainly takes excitation transverse waves as a main mode, the acoustic vibration time domain waveform is composed of two parts, namely direct longitudinal waves with lower amplitude at the front and transverse-longitudinal aliasing waves with higher amplitude at the back; the method comprises the steps that a head wave of a direct longitudinal wave with a lower amplitude value is used as a time delay calculation reference point, the obtained wave velocity is 5000m/s and is very close to a theoretical value (5027m/s) of the wave velocity of the longitudinal wave in a pipe wall, but the actual GIL ultrasonic fault positioning online monitoring system cannot distinguish the direct longitudinal wave waveform with a lower amplitude value buried in a noise signal due to the sensitivity of a used sensor, field noise interference and the like, the aliasing wave with a higher amplitude value at the back is usually selected as the reference point of time delay calculation, and the set reference wave velocity is about 2000 m/s; in order to enable the numerical simulation result to better reflect the test result and further better guide the improvement of the actual system, the factors are considered to determine the effective wave head of the waveform of the numerical simulation result and the time delay calculation reference point on the wave head; therefore, the area between the rising part starting point and the wave head peak point of each wave head in the wave group of the wave form head wave is selected as the distribution area of the time delay calculation reference points and respectively recorded as area1, area2 and area3, the starting point and the ending point of each area are recorded as k11, k12, k21, k22, k31 and k32, the corresponding time points are recorded as t11, t12, t21, t22, t31 and t32, and the time point when the force hammer strikes the GIL shell is recorded as t 0;
step 4-2: calculating the wave velocity range corresponding to each region according to a formula v ═ Δ s/Δ t, wherein Δ s is the distance from the knocking point to the probe point in the model, Δ t is the time difference t11-t0, t12-t0, t21-t0, t12-t0 & gtbetween the corresponding time of the starting point and the ending point of each region and the knocking time of the force hammer, and the calculated wave velocity ranges are v 11-v 12, v 21-v 22 and v 31-v 32 & gtrespectively;
step 4-3: calculating the actually measured wave speed between a knocking point and a sensor according to the actually measured waveform of the vibration time domain obtained in the step 1, wherein the sensor used in the test is a sensor used by the GIL ultrasonic positioning on-line monitoring system, so that the influence of the sensitivity of the system and the interference of external noise is considered in the test result; and (3) comparing the actually measured wave speed obtained by testing in the step (1) with the wave speed range corresponding to each wave head time delay calculation reference point distribution area obtained by calculating in the step (4-2), determining the effective wave head in the wave group of the wave head of the waveform of the surface probe of the model, continuously reducing the range of the time delay calculation reference point distribution area on the wave head, finally determining the position of the time delay calculation reference point, and recording the percentage of the longitudinal coordinate value of the point to the amplitude value of the wave head where the point is located, wherein the percentage is used as the basis for determining the numerical simulation waveform time delay calculation reference point. In the previous step 4-2, the wave velocity range corresponding to each wave head region of the simulation waveform is obtained. In step 4-3, it is determined which of the plurality of wave velocity ranges is met according to the measured wave velocity (only one value) obtained by the test, and the wave head in which the wave velocity range including the measured wave velocity is located is the effective wave head. This step relates the test results to the simulation results.
Step 4-4: and respectively arranging probes on two sides of the straight pipe unit with the same length after the expansion joint of the first finite element simulation model and the bellows of the expansion joint unit are straightened, determining the effective wave head of the wave form of the probes and the position of the time delay calculation reference point on the wave head, and calculating the front and back attenuation and time delay of the sound wave passing through the expansion joint and the straight pipe unit with the same length after the bellows of the expansion joint unit is straightened according to the amplitude of the effective wave head and the time corresponding to the time delay calculation reference point.
And 5: based on the attenuation and time delay characteristics of the sound wave passing through the expansion joint and the straight pipe unit with the same length as the expansion joint unit corrugated pipe after being straightened, a mathematical model of the attenuation and time delay of the GIL expansion joint calculated according to the joint is established, and the method comprises the following steps:
step 5-1: based on the attenuation A1 and the delay T1 of the sound wave passing through the expansion joint and the attenuation A2 and the delay T2 of the sound wave passing through the front and back straight pipe units with the same length as the expansion joint unit corrugated pipe after being straightened, the extra attenuation and delay caused by the expansion joint structure can be calculated and obtained to be A1-A2 and T1-T2 respectively, wherein the A2 comprises the diffusion attenuation and medium absorption attenuation of the sound wave propagation process;
step 5-2: calculating the joint attenuation amount a of the expansion joint, namely (A1-A2)/k, and the joint delay amount T, namely (T1-T2)/k, wherein k is the number of corrugated joints of the expansion joint unit corrugated pipe;
step 5-3: the mathematical model of the attenuation of the GIL expansion joint reduced according to the joints is as follows: a is ka + A0, wherein k is the number of telescopic joints, a is the joint attenuation of the telescopic joints, A0 is the attenuation of the sound wave before and after the sound wave passes through the straight pipe unit with the same length as the telescopic joint unit corrugated pipe after being straightened, and the propagation distance of the sound wave and the sound wave accord with the exponential attenuation law;
step 5-4: a mathematical model of the GIL telescopic joint delay metric, which is calculated according to the joints, is as follows: t is kt + T0, where k is the number of expansion joints, T is the time delay of the expansion joint, and T0 is the time delay before and after the sound wave passes through the straight pipe unit with the same length as the expansion joint unit corrugated pipe after being straightened, and is proportional to the length of the straight pipe unit.
The invention provides a node-by-node reduction algorithm of GIL telescopic joint attenuation and time delay under the consideration of the influence of noise interference and sensor sensitivity on results in an actual system, which can calculate the attenuation and time delay of sound waves passing through different types of GIL telescopic joints, thereby providing a basis for the improvement of a GIL ultrasonic fault positioning online monitoring method, and has strong universality and wide engineering application prospect.
It is noted that references in the specification to "one embodiment," "an example embodiment," "some embodiments," etc., indicate that the embodiment described may include a particular feature, structure, or characteristic, but every embodiment may not necessarily include the particular feature, structure, or characteristic. Moreover, such phrases are not necessarily referring to the same embodiment. Further, when a particular feature, structure, or characteristic is described in connection with an embodiment, it is submitted that it is within the knowledge of one skilled in the art to effect such feature, structure, or characteristic in connection with other embodiments whether or not explicitly described.
The previous description of the disclosure is provided to enable any person skilled in the art to make or use the disclosure. Various modifications to the disclosure will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other variations without departing from the spirit or scope of the disclosure. Thus, the disclosure is not intended to be limited to the examples and designs described herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
Claims (7)
1. A GIL telescopic joint acoustic wave transmission modeling method is characterized by comprising the following steps:
the method comprises the following steps of (1) testing acoustic transmission characteristics, namely knocking a shell of a GIL transmission system comprising an expansion joint, and recording a knocking force waveform signal and an actually measured vibration time domain waveform signal transmitted by sound waves in the GIL;
establishing a first finite element simulation model and a second finite element simulation model representing the geometric form of the GIL transmission system, wherein elements and the geometric dimensions in the first finite element simulation model are the same as those of the GIL transmission system, and the second finite element simulation model is obtained by replacing an expansion joint corrugated pipe in the first finite element simulation model with a GIL straight pipe with the same length as that of the expanded corrugated pipe;
a finite element numerical simulation step, fitting according to the actually measured knocking force waveform signal to obtain an excitation force waveform; applying the excitation force waveform to the first finite element simulation model and the second finite element simulation model in a total force mode, and recording simulated vibration time domain waveform signals obtained by simulation in each finite element model;
analyzing an effective wave head of the simulated vibration time domain waveform signal and a time delay calculation reference point on the wave head, and calculating the attenuation and time delay characteristics of the sound wave after passing through the first finite element simulation model and the second finite element simulation model;
and acoustic transmission modeling, namely establishing a mathematical model of the attenuation and delay of the GIL expansion joint reduced by sections based on the attenuation and delay characteristics.
2. The GIL telescopic joint sound wave transmission modeling method according to claim 1, wherein in the sound transmission characteristic testing step, a standard knocking sound source is generated in a mode that a piezoelectric impact hammer knocks a GIL shell, and two sensors arranged in front of and behind a GIL telescopic joint unit are selected to record vibration time domain waveform signals transmitted by sound waves in the GIL.
3. The GIL telescopic joint acoustic wave transmission modeling method as claimed in claim 1, wherein the finite element model modeling step adds an acoustic-structure coupling boundary condition, a bracket constraint boundary condition, a medium damping characteristic and an infinite length boundary condition at both ends of the model to the finite element simulation model.
4. The GIL telescopic joint acoustic wave transmission modeling method as claimed in claim 1, wherein in said finite element numerical simulation step, an excitation position at which said excitation force waveform is applied to the GIL finite element simulation model and a measurement position at which a simulated vibration time domain waveform signal is recorded correspond to a tapping position and a recording position in said acoustic transmission characteristic test step.
5. The GIL telescopic joint acoustic wave transmission modeling method as claimed in claim 1, wherein in the finite element numerical simulation step, an excitation force waveform is obtained by fitting the actually measured knocking force waveform signal in a manner of piecewise splicing of a quadratic function and a Gaussian pulse function.
6. The GIL telescopic node acoustic wave transmission modeling method as claimed in claim 1, wherein in the attenuation and time delay analysis step, in a wave group of a horizontal-longitudinal aliasing wave, a region between a start point of a rising part of each wave head in a wave group of a wave head of the wave shape and a peak point of the wave head is selected as a distribution region of a time delay calculation reference point, a wave velocity range corresponding to each region is calculated, an effective wave head in the wave group of the wave shape of the simulated vibration time domain waveform signal is determined based on an actually measured wave velocity calculated by the actually measured vibration time domain waveform signal to determine a position of the time delay calculation reference point, and an attenuation amount and a time delay amount of the acoustic wave after passing through the first finite element simulation model and the second finite element simulation model are calculated based on the position of the reference point.
7. The GIL telescopic joint sound wave transmission modeling method of claim 1, wherein an additional attenuation A1-A2 and a delay T1-T2 caused by the telescopic joint structure are calculated based on an attenuation A1 and a delay T1 of sound waves passing through the first finite element simulation model and an attenuation A2 and a delay T2 before and after passing through the second finite element simulation model; calculating the joint attenuation amount and the time delay amount of a single telescopic joint based on the extra attenuation amount A1-A2 and the time delay amount T1-T2 caused by the telescopic joint structure; and establishing a mathematical model of the attenuation and the time delay of the GIL expansion joint reduced according to the sections based on the section attenuation and the section time delay.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110507493.8A CN113297759B (en) | 2021-05-10 | 2021-05-10 | GIL telescopic joint acoustic transmission modeling method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110507493.8A CN113297759B (en) | 2021-05-10 | 2021-05-10 | GIL telescopic joint acoustic transmission modeling method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113297759A true CN113297759A (en) | 2021-08-24 |
| CN113297759B CN113297759B (en) | 2022-04-15 |
Family
ID=77321327
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110507493.8A Active CN113297759B (en) | 2021-05-10 | 2021-05-10 | GIL telescopic joint acoustic transmission modeling method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113297759B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114910753A (en) * | 2022-04-21 | 2022-08-16 | 武汉朗德电气有限公司 | Acoustic wave transmission device and detection method for GIL fault location |
| CN116699339A (en) * | 2023-08-04 | 2023-09-05 | 武汉朗德电气有限公司 | GIL arc fault positioning method based on bimodal acoustic guided wave |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130145842A1 (en) * | 2010-03-09 | 2013-06-13 | Cidra Corporate Services Inc. | Dispersion Compensation Technique for Differential Sonar Measurement - Density Meter. |
| CN108830025A (en) * | 2018-08-21 | 2018-11-16 | 广西电网有限责任公司电力科学研究院 | A kind of simulating analysis of guided wave propagation characteristic in energy-stored spring |
| CN110703051A (en) * | 2019-10-15 | 2020-01-17 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | Precise positioning method for fault point of GIL (Gate in L) voltage withstand test based on acoustic wave detection method |
| CN111460706A (en) * | 2020-03-19 | 2020-07-28 | 国网江苏省电力有限公司电力科学研究院 | High-voltage alternating current GI L pipe gallery online monitoring and temperature state distinguishing method, system and storage medium thereof |
| CN111896848A (en) * | 2020-06-24 | 2020-11-06 | 国网江苏省电力有限公司电力科学研究院 | Internal discharge positioning method, terminal, system and storage medium of power transmission line |
-
2021
- 2021-05-10 CN CN202110507493.8A patent/CN113297759B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130145842A1 (en) * | 2010-03-09 | 2013-06-13 | Cidra Corporate Services Inc. | Dispersion Compensation Technique for Differential Sonar Measurement - Density Meter. |
| CN108830025A (en) * | 2018-08-21 | 2018-11-16 | 广西电网有限责任公司电力科学研究院 | A kind of simulating analysis of guided wave propagation characteristic in energy-stored spring |
| CN110703051A (en) * | 2019-10-15 | 2020-01-17 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | Precise positioning method for fault point of GIL (Gate in L) voltage withstand test based on acoustic wave detection method |
| CN111460706A (en) * | 2020-03-19 | 2020-07-28 | 国网江苏省电力有限公司电力科学研究院 | High-voltage alternating current GI L pipe gallery online monitoring and temperature state distinguishing method, system and storage medium thereof |
| CN111896848A (en) * | 2020-06-24 | 2020-11-06 | 国网江苏省电力有限公司电力科学研究院 | Internal discharge positioning method, terminal, system and storage medium of power transmission line |
Non-Patent Citations (2)
| Title |
|---|
| 刘云鹏等: "特高压GIL故障定位超声衰减特性及试验研究", 《电网技术》 * |
| 宋超然: "高压交流GIL多物理场耦合分析及在线监测系统设计", 《中国优秀硕士学位论文全文数据库 (工程科技Ⅱ辑)》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114910753A (en) * | 2022-04-21 | 2022-08-16 | 武汉朗德电气有限公司 | Acoustic wave transmission device and detection method for GIL fault location |
| CN116699339A (en) * | 2023-08-04 | 2023-09-05 | 武汉朗德电气有限公司 | GIL arc fault positioning method based on bimodal acoustic guided wave |
| CN116699339B (en) * | 2023-08-04 | 2023-11-17 | 武汉朗德电气有限公司 | GIL arc fault positioning method based on bimodal acoustic guided wave |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113297759B (en) | 2022-04-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113297759B (en) | GIL telescopic joint acoustic transmission modeling method | |
| CN103344934B (en) | The detection check method and system of Partial Discharge in Power Transformer sonac | |
| CN102706317B (en) | Online monitoring device for thermal expansion amount of pressure-containing member of power station boiler | |
| CN106556647A (en) | A kind of impact echo data processing method | |
| CN104537251B (en) | A kind of shock load recognition methods of fan blade | |
| CN112162034B (en) | Steel-concrete combined beam damage identification method applying structural noise | |
| Setshedi et al. | Estimation of rail properties using semi-analytical finite element models and guided wave ultrasound measurements | |
| CN109299539A (en) | A method of utilizing curvilinear corner index evaluation cantilever beam rigidity weak link | |
| JP3719740B2 (en) | Nondestructive testing method and apparatus using ultrasonic waves | |
| CN111157065A (en) | Acoustic time delay measuring method in ultrasonic signal transmission loop of gas ultrasonic flowmeter | |
| TW201310053A (en) | On-site instant earthquake analysis system and method, and storage medium thereof | |
| Mu et al. | Echo signal envelope fitting based signal processing methods for ultrasonic gas flow-meter | |
| CN113932942B (en) | Method and device for detecting temperature field inside cable by ultrasonic wave | |
| CN114910753B (en) | Acoustic wave transmission device and detection method for GIL fault location | |
| CN108362431A (en) | Non-intervention type pressure detection method based on time delay spacing between adjacent longitudinal wave and measuring system | |
| WO2021057288A1 (en) | Pipe creep measurement system and method | |
| CN113624395A (en) | Bolt axial force detection coefficient calibration method based on ultrasound | |
| JP2004061361A (en) | Piping breakage investigating apparatus | |
| JP3421412B2 (en) | Pipe thinning measurement method and equipment | |
| CN115356142B (en) | Hydraulic rock drill performance detection system and method | |
| CN115420407B (en) | CFRP plate load identification method based on FBG monitoring response information | |
| CN105548356A (en) | Method for detecting defects of small-bending-radius bend pipe with girth joint based on T-mode guided waves | |
| CN114397072B (en) | Site damage assessment method for static test and thermal strength test of structure | |
| CN212275668U (en) | Scanning type impact echo method concrete structure reinforcing steel bar nondestructive testing equipment | |
| CN109855536B (en) | Oil and gas pipeline blockage detection method based on strain measurement |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |