CN109782136A - A method of contact net defective locations are determined based on longitudinal wave shear wave velocity difference - Google Patents
A method of contact net defective locations are determined based on longitudinal wave shear wave velocity difference Download PDFInfo
- Publication number
- CN109782136A CN109782136A CN201811615587.1A CN201811615587A CN109782136A CN 109782136 A CN109782136 A CN 109782136A CN 201811615587 A CN201811615587 A CN 201811615587A CN 109782136 A CN109782136 A CN 109782136A
- Authority
- CN
- China
- Prior art keywords
- wave
- shear wave
- duration
- value
- defective locations
- 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
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The invention discloses a kind of method for determining contact net defective locations based on longitudinal wave shear wave velocity difference, includes the following steps: S1, determines the time difference of longitudinal wave and shear wave arrival test point at defective locations;S2, the defective locations are calculated at a distance from test point.Specifically, by choosing the P-wave amplitude maximum value in multiple certain period of times, corresponding longitudinal wave repeatedly matches the time point of shear wave, obtain the propagation time difference of longitudinal wave and shear wave, so that it is determined that defective locations are at a distance from test point, the accurate positioning for realizing defective locations, improves efficiency, reduces costs.
Description
Technical field
The present invention relates to contact net technical fields, determine contact net based on longitudinal wave shear wave velocity difference more particularly, to one kind
The method of defective locations.
Background technique
Contact net is to set up along rail overhead " it " font in electric railway, take the high pressure of stream defeated for pantograph
Electric wire.Contact net is the main truss of railway electrification project, is that powering to electric locomotive for downline overhead erection is special
The transmission line of electricity of form.The quality and working condition of contact net will directly affect the transport capacity of electric railway, therefore,
It no matter in any condition to seek contact net, can guarantee to supply electric power locomotive electric energy well, guarantee electric locomotive on the line
Safety, high-speed cruising, and in the case where meeting above-mentioned requirements, it reduces investment outlay as much as possible, is structurally reasonable, maintenance simplicity, leading to
Often in the position of contact net anchorage, the hard of contact line can be monitored by installing vibrating sensor according to the detection data of vibrating sensor
The defects of point, crackle, but determining defect is typically only capable in which anchor section, and can not accurately determine the position of defect.It analyzes latent
Defect may, but can not determine the accurate location of defect, when service worker's site examining and repairing, to devote a tremendous amount of time
Confirm the position of defect, it is not easy to look for, or even be also possible to can not find, waste a large amount of manpower and material resources.
Summary of the invention
The object of the present invention is to provide a kind of methods for determining contact net defective locations based on longitudinal wave shear wave velocity difference, lead to
It crosses the longitudinal wave for determining same vibration source and shear wave reaches the defective locations time difference, determine defective locations, improve efficiency, reduce manpower
The investment of material resources.
Foregoing invention purpose of the invention has the technical scheme that
A method of contact net defective locations are determined based on longitudinal wave shear wave velocity difference, are included the following steps:
S1, the time difference for determining longitudinal wave and shear wave propagation arrival test point at defective locations;
S2, the defective locations are calculated at a distance from test point.
The present invention is further arranged to: in step S1, being included the following steps:
Sampled data in K A1, selection the first durations, the P-wave amplitude average value and shear wave sought in each first duration shake
Width average value;
A2, the P-wave amplitude average value maximum value present position is chosen as first position i, the correspondence moment was denoted as the first moment
t1;
A3, at the first position i, seek the first P-wave amplitude average value in M continuous second duration and first horizontal
Wave-amplitude average value successively obtains the ratio of the first P-wave amplitude average value and the first shear wave mean value of amplitude, and to all ratios
Value seeks mean square deviation δ 1;
A4, after moving forward or moving back one second duration, the second shear wave mean value of amplitude of M continuous second durations is sought again,
The ratio of the first P-wave amplitude average value and the second shear wave mean value of amplitude is successively obtained, and mean square deviation δ is sought to all ratios
2;
A5, step A4 is repeated, is more than the entire anchor Duan Suoxu of sound transmission until continuous second duration of all M corresponding time
Time;
Minimum Mean Square Error δ min in A6, all mean square deviations of selection;
A7, when Minimum Mean Square Error δ min is less than setting value, then it is assumed that shear wave this moment t2 at at the first moment t1
Longitudinal wave be the same vibration source vibration wave;
The time difference of A8, the first moment t1 and moment t2, the propagation time difference d of the as described longitudinal wave and the shear wave.
The present invention is further arranged to: the value range of the M is less than 200, and the value range of the K is less than 100.
The present invention is further arranged to: the value of the M is equal to 50.
The present invention is further arranged to: first duration is equal to second duration.
The present invention is further arranged to: being chosen M continuous second durations at the first position, is in first position
Front and back respectively choose M/2 the second durations.
The present invention is further arranged to: first duration or the second duration are no more than 0.3 second.
The present invention is further arranged to: when Minimum Mean Square Error δ min is more than or equal to setting value, then it is assumed that selected hits
According to for interference signal and abandoning, sampled data is chosen again and is analyzed.
The present invention is further arranged to: at the moment t2, being chosen third duration and is repeated step A4-A7, the third
Duration is less than second duration, improves positioning accuracy.
The present invention is further arranged to: in step S2, defective locations and test point distance S are obtained by following formula:
S=v1*v2*d/(v1-v2)
Wherein, v1 indicates the spread speed of longitudinal wave;V2 indicates the spread speed of shear wave, and d indicates the propagation time of longitudinal wave and shear wave
Difference.
In conclusion advantageous effects of the invention are as follows:
1. the present invention by the time difference to same vibration source longitudinal wave and shear wave propagation, accurately determines defective locations, improves efficiency,
The investment for reducing manpower and material resources, reduces cost;
2. further, for the present invention by comparing to sampled signal in time domain, acquisition shear wave and longitudinal wave correct time are poor,
Improve the accuracy of defective locations judgement.
3. further, the present invention by after coarse positioning further progress be accurately positioned, not only reduce calculation amount, but also guarantee
Positioning accuracy.
Detailed description of the invention
Fig. 1 is the method flow schematic diagram of a specific embodiment of the invention.
Fig. 2 is the determination longitudinal wave and shear transit time difference flow diagram of a specific embodiment of the invention.
Fig. 3 is the longitudinal wave and shear wave waveform diagram of a specific embodiment of the invention.
Specific embodiment
Below in conjunction with attached drawing, invention is further described in detail.
A kind of method that contact net defective locations are determined based on longitudinal wave shear wave velocity difference of the present invention, flow chart such as Fig. 1
It is shown, include the following steps:
S1, the time difference for determining longitudinal wave and shear wave propagation arrival test point at defective locations;
S2, the defective locations are calculated at a distance from test point.
When vibration sampling, while the shear wave and longitudinal wave of vibration are acquired, shear wave refers to the vibration perpendicular to conductive line surfaces, and longitudinal wave refers to
It is parallel to the vibration in conducting wire direction;Assuming that the spread speed of longitudinal wave is v1, propagation time t1;The spread speed of shear wave is v2,
Propagation time is t2;Wherein longitudinal wave is main related with material, temperature with the spread speed of shear wave, can measure in advance.
Propagation time t1 and t2 can not be measured directly, but due to the longitudinal wave same defective locations corresponding with shear wave, i.e. longitudinal wave
With shear wave generate at the time of be it is identical, their time difference d can be measured, t2=t1+d.
For a certain defect, " amplitude over time " characteristic of shear wave and longitudinal wave is similar, as shown in figure 3, in shear wave and indulging
Wave number finds the same vibration source (i.e. defective locations) in, as shown in Fig. 2, including the following steps:
Sampled data in K A1, selection the first durations, the P-wave amplitude average value and shear wave sought in each first duration shake
Width average value;
A2, the P-wave amplitude average value maximum value present position is chosen as first position i, the correspondence moment was denoted as the first moment
t1;
A3, at the first position i, seek the first P-wave amplitude average value in M continuous second duration and first horizontal
Wave-amplitude average value successively obtains the ratio of the first P-wave amplitude average value and the first shear wave mean value of amplitude, and to all ratios
Value seeks mean square deviation δ 1;
A4, after moving back one second duration, the second shear wave mean value of amplitude of M continuous second durations is sought again, successively
Mean square deviation δ 2 is sought to the ratio of the first P-wave amplitude average value and the second shear wave mean value of amplitude, and to all ratios;
A5, step A4 is repeated, is more than the entire anchor Duan Suoxu of sound transmission until continuous second duration of all M corresponding time
Time;
Minimum Mean Square Error δ min in A6, all mean square deviations of selection;
A7, when Minimum Mean Square Error δ min is less than setting value, then it is assumed that shear wave this moment t2 at at the first moment t1
Longitudinal wave be the same vibration source vibration wave;
The time difference of A8, the first moment t1 and moment t2, the propagation time difference d of the as described longitudinal wave and the shear wave.
In a specific embodiment of the invention, the sampled data in vibration sampling database in K period is chosen, i.e.,
It opens a window to sampling database, i.e. first window, the time tk of each first window is identical, and first window is in sampling database
In be dispersion, seek the P-wave amplitude average value in each first window, be respectively labeled as f1, f2 ... fk.Therefrom choose
The maximum position i of P-wave amplitude average value is t1 at the time of corresponding, then takes M continuous second windows in total before and after the t1 moment
Mouthful, the time tm of each second window is identical, seeks the P-wave amplitude average value fm of sampled data and cross in this M window respectively
Wave-amplitude average value hm, is respectively labeled as: fm (1), fm (2) ... fm (m), hm (1), hm (2) ... hm (m) is successively right
The P-wave amplitude average value and shear wave mean value of amplitude of same sequence number carry out division arithmetic, obtain longitudinal and shear wave amplitude scale p1, point
Not Wei p11=fm (1)/hm (1), p12=fm (2)/hm (2) ... p1m=fm (m)/hm (m), then to it is all in length and breadth
Wave-amplitude ratio p1 carries out mean square deviation operation, obtains the first mean square deviation δ 1.
The second window of M connection is chosen at the t1 moment, it is preferable that choose M/2 connection respectively before and after the t1 moment
The second window.
The sampled data of shear wave is moved into a second window time tm backward, recalculates shear wave mean value of amplitude, point
Biao Ji not be (2) ... hm (m), hm (m+1), successively at the P-wave amplitude average value at t1 moment and t1+tm moment
Shear wave mean value of amplitude carries out division arithmetic, obtains longitudinal and shear wave amplitude scale p2, and respectively p21=fm (1)/hm (2), p22=
Fm (2)/hm (3) ... p2m=fm (m)/hm (m+1), then mean square deviation operation is carried out to all longitudinal and shear wave amplitude scale p2,
Obtain the second mean square deviation δ 2;
The rest may be inferred, calculates q mean square deviation operation, and the value for obtaining q mean square deviation δ q, q is propagated in entire anchor section by sound
Required time T is determined, if the length of anchor section is G meters, the minimum speed of sound transmission is 3000 meter per seconds, then sound is in anchor section
T=G/3000 the time required to propagating, thus obtains q=T/tm.
The smallest mean square deviation δ min is sought in q mean square deviation δ 1, δ 2 ... δ q, when mean square deviation δ min is less than setting value, then
The longitudinal wave of the shear wave and t1 moment of thinking this moment is generated by the same vibration source, is remembered at the time of mean square deviation δ min is corresponded to
For t2, time difference d=t2-t1 of longitudinal wave and shear wave arrival test point.
Preferably, the length of anchor section is G=1500 meter, then T is up to 0.5 second.Such as choosing tm is 10 milliseconds, is carried out aforementioned
It calculates, then positioning accuracy is about 75 meters, is at this time coarse positioning.
After coarse positioning, window time tm is reduced, 1 millisecond is such as chosen for, carries out above-mentioned meter near the t2 moment above-mentioned
It calculates, it can be achieved that positioning accuracy is about 7.5 meters.
When mean square deviation δ min is more than or equal to setting value, then it is assumed that sampled data at this time is interference signal, is abandoned, again
Sampled data is chosen to be analyzed.
Preferably, the value range of M is less than 200, and the value range of the K is less than 100.
Further, the value of M is equal to 50.
After obtaining the propagation time difference d of longitudinal wave and shear wave, the distance S of calculating defective locations to test point:
S=v1*t1
S=v2*t2
Solution obtains:
S=v1*v2*d/(v1-v2)。
Thus formula can be seen that the range accuracy of defective locations and test point by longitudinal wave, shear wave spread speed measurement
Precision and the measurement accuracy of time difference d determine.
The embodiment of present embodiment is presently preferred embodiments of the present invention, not limits protection of the invention according to this
Range, therefore: the equivalence changes that all structures under this invention, shape, principle are done, should all be covered by protection scope of the present invention it
It is interior.
Claims (10)
1. a kind of method for determining contact net defective locations based on longitudinal wave shear wave velocity difference, it is characterised in that: including walking as follows
It is rapid:
S1, the time difference for determining longitudinal wave and shear wave arrival test point at defective locations;
S2, the defective locations are calculated at a distance from test point.
2. according to the method described in claim 1, it is characterized by: including the following steps: in step S1
Sampled data in K A1, selection the first durations, the P-wave amplitude average value and shear wave sought in each first duration shake
Width average value;
A2, the maximum value present position for choosing the P-wave amplitude average value are first position i, when the correspondence moment is denoted as first
Carve t1;
A3, at the first position i, seek the first P-wave amplitude average value in M continuous second duration and first horizontal
Wave-amplitude average value successively obtains the ratio of the first P-wave amplitude average value and the first shear wave mean value of amplitude, and to all ratios
Value seeks mean square deviation δ 1;
A4, after moving forward or moving back one second duration, the second shear wave mean value of amplitude of M continuous second durations is sought again,
The ratio of the first P-wave amplitude average value and the second shear wave mean value of amplitude is successively obtained, and mean square deviation δ is sought to all ratios
2;
A5, step A4 is repeated, is more than the entire anchor Duan Suoxu of sound transmission until continuous second duration of all M corresponding time
Time;
Minimum Mean Square Error δ min in A6, all mean square deviations of selection;
A7, when Minimum Mean Square Error δ min is less than setting value, then it is assumed that shear wave this moment t2 at at the first moment t1
Longitudinal wave be the same vibration source vibration wave;
The time difference of A8, the first moment t1 and moment t2, the propagation time difference d of the as described longitudinal wave and the shear wave.
3. according to the method described in claim 2, it is characterized by: the value range of the M less than 200, the value model of the K
It encloses less than 100.
4. according to the method described in claim 3, it is characterized by: the value of the M is equal to 50.
5. according to the method described in claim 2, it is characterized by: first duration is equal to second duration.
6. according to the method described in claim 2, it is characterized by: when choosing M continuous second at the first position
It is long, it is that M/2 the second durations are respectively chosen before and after first position.
7. according to the method described in claim 2, it is characterized by: first duration or the second duration are no more than 0.3 second.
8. according to the method described in claim 2, it is characterized by: when Minimum Mean Square Error δ min be more than or equal to setting value when, then
Think that selected sampled data is interference signal and abandons, chooses sampled data again and analyzed.
9. according to the method described in claim 2, it is characterized by: choosing third duration at the moment t2 and repeating step
A4-A7, the third duration improve positioning accuracy less than the second duration.
10. according to the method described in claim 1, it is characterized by: in step S2, defective locations and test point distance S, by
Following formula obtains:
S=v1*v2*d/(v1-v2)
Wherein, v1 indicates the spread speed of longitudinal wave;V2 indicates the spread speed of shear wave, and d indicates the propagation time of longitudinal wave and shear wave
Difference.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811615587.1A CN109782136B (en) | 2018-12-27 | 2018-12-27 | Method for determining defect position of contact net based on longitudinal wave and transverse wave speed difference |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811615587.1A CN109782136B (en) | 2018-12-27 | 2018-12-27 | Method for determining defect position of contact net based on longitudinal wave and transverse wave speed difference |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109782136A true CN109782136A (en) | 2019-05-21 |
CN109782136B CN109782136B (en) | 2021-04-09 |
Family
ID=66497894
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811615587.1A Active CN109782136B (en) | 2018-12-27 | 2018-12-27 | Method for determining defect position of contact net based on longitudinal wave and transverse wave speed difference |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109782136B (en) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001021607A (en) * | 1999-07-12 | 2001-01-26 | Sumitomo Electric Ind Ltd | Device for detecting partial discharge generating position |
CN102707190A (en) * | 2012-01-10 | 2012-10-03 | 成都唐源电气有限责任公司 | Direct-current-side short-circuit fault distance measuring device and method of metro tractive power supply system |
CN105203917A (en) * | 2015-09-09 | 2015-12-30 | 昆明理工大学 | Through type cophase traction power supply system traction network fault interval traveling wave positioning method |
CN105738766A (en) * | 2016-02-23 | 2016-07-06 | 武汉大学 | Overhead contact system single-ended traveling wave fault positioning device |
US20180088164A1 (en) * | 2016-09-26 | 2018-03-29 | KW Associates LLC | Estimation of arc location in three dimensions |
US20180136269A1 (en) * | 2016-11-11 | 2018-05-17 | Schweitzer Engineering Laboratories, Inc. | Traveling wave based single end fault location |
CN108761273A (en) * | 2018-05-25 | 2018-11-06 | 国家电网公司 | A kind of traveling wave fault positioning method not influenced by velocity of wave |
-
2018
- 2018-12-27 CN CN201811615587.1A patent/CN109782136B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001021607A (en) * | 1999-07-12 | 2001-01-26 | Sumitomo Electric Ind Ltd | Device for detecting partial discharge generating position |
CN102707190A (en) * | 2012-01-10 | 2012-10-03 | 成都唐源电气有限责任公司 | Direct-current-side short-circuit fault distance measuring device and method of metro tractive power supply system |
CN105203917A (en) * | 2015-09-09 | 2015-12-30 | 昆明理工大学 | Through type cophase traction power supply system traction network fault interval traveling wave positioning method |
CN105738766A (en) * | 2016-02-23 | 2016-07-06 | 武汉大学 | Overhead contact system single-ended traveling wave fault positioning device |
US20180088164A1 (en) * | 2016-09-26 | 2018-03-29 | KW Associates LLC | Estimation of arc location in three dimensions |
US20180136269A1 (en) * | 2016-11-11 | 2018-05-17 | Schweitzer Engineering Laboratories, Inc. | Traveling wave based single end fault location |
CN108761273A (en) * | 2018-05-25 | 2018-11-06 | 国家电网公司 | A kind of traveling wave fault positioning method not influenced by velocity of wave |
Non-Patent Citations (1)
Title |
---|
张礼得: "接触网单端故障行波定位方法研究", 《武汉大学学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN109782136B (en) | 2021-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP7411655B2 (en) | Method and system for monitoring track sections | |
CN105783799B (en) | A kind of fragment-free track slab gap depth lossless detection method and equipment based on vibration | |
CN102798804B (en) | High-voltage power cable fault on-line positioning device | |
CN102288883B (en) | Oscillation wave partial discharge identifying and positioning method for asynchronous double-end power cable | |
JP3524431B2 (en) | measuring device | |
CN102809715A (en) | On-line fault location device for high-voltage power cable | |
CN105571830B (en) | The method for measuring super-narrow line width laser device laser line width | |
CN106043356A (en) | Urban rail train rim vertex circle diameter online detection method and system | |
CN104631516B (en) | A kind of portable reinforced earth retaining walls facial disfigurement monitoring device and monitoring method | |
CN108519444A (en) | A kind of accurate measurement method of contact line defective locations | |
CN104931218A (en) | System of modal vibration mode of non-contact measuring flexible structure | |
CN101762347A (en) | Method for measuring rope force of multi-span steel stay rope by using half-wave method | |
Tanaka et al. | Practical application of portable trolley for the continuous measurement of rail surface roughness for rail corrugation maintenance | |
CN109855771B (en) | Method for detecting temperature based on vertical acceleration power spectrum density of steel rail | |
CN110470422B (en) | Fundamental frequency optimization method based on inhaul cable vibration spectrum analysis | |
CN204903118U (en) | System for non -contact measures flexible construction's mode mode of vibration | |
Jiang et al. | Assessment of pantograph-catenary interaction in a railway overlap section via a novel optical-based method | |
US20090320598A1 (en) | Method for testing the microstructure of a welded joint | |
Marendić et al. | Measurement of bridge dynamic displacements and natural frequencies by RTS | |
CN109693686A (en) | A kind of track alignment irregularity measurement method based on double antenna GNSS/INS | |
CN109782136A (en) | A method of contact net defective locations are determined based on longitudinal wave shear wave velocity difference | |
Song et al. | Current Collection Quality of High-speed Rail Pantograph-catenary Considering Geometry Deviation at 400 km/h and Above | |
CN107192763B (en) | Utilize the test method for the rail vibration attenuation rate that running train is excitation | |
CN111044848B (en) | Cable fault high-precision positioning method based on characteristic parameter adjustment wave speed | |
CN103123303B (en) | Quantifying and online monitoring method of bridge girder safe reliability |
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 |