CN110514862A - A kind of high-speed rail speed of service estimation method using velocity scanning - Google Patents

A kind of high-speed rail speed of service estimation method using velocity scanning Download PDF

Info

Publication number
CN110514862A
CN110514862A CN201910760699.4A CN201910760699A CN110514862A CN 110514862 A CN110514862 A CN 110514862A CN 201910760699 A CN201910760699 A CN 201910760699A CN 110514862 A CN110514862 A CN 110514862A
Authority
CN
China
Prior art keywords
speed
amplitude spectrum
signal
cross
train
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910760699.4A
Other languages
Chinese (zh)
Other versions
CN110514862B (en
Inventor
王晓凯
陈文超
师振盛
刘璞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201910760699.4A priority Critical patent/CN110514862B/en
Publication of CN110514862A publication Critical patent/CN110514862A/en
Application granted granted Critical
Publication of CN110514862B publication Critical patent/CN110514862B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P3/00Measuring linear or angular speed; Measuring differences of linear or angular speeds
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/303Analysis for determining velocity profiles or travel times

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Train Traffic Observation, Control, And Security (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

The invention discloses a kind of high-speed rail speed of service estimation methods using velocity scanning, design amplitude spectrum stencil function according to the canonical parameter of the speed of service and China's high-speed rail train;The signal excited when high-speed rail is passed through is intercepted in single detector seismic data;Fourier transformation is done to the signal of interception, obtains its amplitude spectrum;The cumulative function for calculating this signal energy spectrum determines the frequency separation where its most of energy;Amplitude spectrum stencil function and actual amplitude spectrum are done into cross-correlation;Maximum cross-correlation coefficient and its corresponding speed are found in all cross-correlation coefficients;Using this speed as basic point, velocity variations interval is reduced, amplitude spectrum stencil function is re-established;The cross-correlation calculation of amplitude spectrum stencil function and actual amplitude spectrum is carried out again;Speed corresponding to maximum cross-correlation coefficient is found, which is denoted as to the final estimating speed v of train operationfinal

Description

High-speed rail running speed estimation method utilizing speed scanning
Technical Field
The invention belongs to the technical field of exploration geophysics, and particularly relates to a high-speed rail running speed estimation method utilizing speed scanning.
Background
In 1964, commercial operation of the new main line in japan opened a new era of world high-speed rail development. Subsequently, countries such as france, germany, canada, italy, sweden, and korea strived to build high-speed rails and opened commercial operations. The first commercial operation high-speed rail line, Jingjin intercity railway, was opened in China in 2008, 8 months and 1 day. At present, the business mileage of the Chinese high-speed rail reaches 3.1 kilometres, which is close to 70 percent of the total business mileage of the world high-speed rail. And from 7 month and 10 day zero in 2019, implementing a new train operation diagram on the railways in China, and running the motor train unit trains by 3310 pairs every day. The method is one of important means for monitoring the running safety of the train, and is also the key for subsequently utilizing the high-speed train to cause a vibration signal. The existing method for acquiring the running speed of the high-speed train mainly comprises the following steps:
prior art 1: acquiring the train running speed by utilizing the vehicle-mounted equipment: train operating speed can be obtained directly by using a tachometer on the train, but the speed of the train when passing a certain position cannot be determined. In addition, a GPS device carried by a high-speed train can provide train speed and real-time train location. The required equipment is installed on the train, and therefore, the permission of the high-speed railway department is required.
Prior art 2: the method comprises the following steps of installing video, optical, radar and other equipment in a high-speed rail line isolation area: common external velocity measurement systems and methods include a camera-based velocity estimation system, an optical sensor-based or two-vibration sensor-based train velocity estimation method, a radar velocity measurement method using the doppler effect, and a wheel count-based velocity estimation method. The above method requires installation in a position where the rails can be seen, and permission to access the isolation zone and install equipment in the isolation zone.
Disclosure of Invention
The technical problem to be solved by the present invention is to provide a method for estimating the operating speed of a high-speed rail by using speed scanning, which aims at overcoming the defects in the prior art, and uses seismic data acquired by a single detector outside a high-speed rail line isolation region to estimate the operating speed of a train by using the speed scanning method, so as to provide data for subsequently judging the operating state of the train.
The invention adopts the following technical scheme:
a high-speed rail running speed estimation method utilizing speed scanning comprises the following steps:
s1, designing an amplitude spectrum template function | F (omega, v) | according to the running speed and typical parameters of the high-speed train;
s2, intercepting the signal excited by the passing of the high-speed rail in the seismic data of the single detector, wherein the time range corresponding to the effective signal acquired by the detector is [ t [ [ t ]1,t2];
S3, carrying out Fourier transform on the intercepted signal to obtain an amplitude spectrum of the intercepted signal;
s4, calculating an energy spectrum accumulation function of the signal intercepted in the step S3, and determining an energy frequency interval;
s5, performing cross-correlation calculation on the amplitude spectrum template function and the actual amplitude spectrum;
s6, finding the maximum cross correlation coefficient in all the cross correlation coefficients of the step S5, and the corresponding speed vr
S7, taking the speed determined in the step S6 as a base point, reducing the speed change interval, and reestablishing an amplitude spectrum template function;
s8, performing cross-correlation calculation of the amplitude spectrum template function and the actual amplitude spectrum again;
s9, finding the speed corresponding to the maximum cross correlation coefficient, and recording the speed as the final estimated speed v of train operationfinal
Specifically, in step S1, the amplitude spectrum template function | F (ω, v) | is:
v=vinitial+mΔv1
wherein D is the length of each car body, v is the preset train speed, vinitialIs the initial scanning speed, m is an index of the scanning speed, Δ v1The initial velocity scan interval.
Further, the number N of carriages is set to be 16, the length D of each car body is 25 m, vinitialSet to 30 km/h, the value range of m is [0,370 ]],Δv1At 1 km/h, a total of 371 functions were included in a series of amplitude spectra template functions.
Specifically, in step S3, assuming that the signal resulting from the high-speed rail operation is truncated is Y (t), the signal is fourier-transformed, and the amplitude spectrum | Y (ω) | is obtained as:
wherein, [ t ]1,t2]And y (t) is a time range corresponding to the effective signal acquired by the detector, the intercepted signal caused by the running of the high-speed rail, j is an imaginary unit omega and is frequency, and t is speed.
Specifically, in step S4, the total energy E of the signal is obtained firstyThen, the upper frequency bound of the signal is calculated, and the square of the amplitude is accumulated from zero value of the frequency until the accumulated value is equal to the total energy EyIs higher than the set retention η1The corresponding frequency is the upper frequency bound omegamax(ii) a When acquiring the lower frequency bound, the accumulation is started from the upper frequency bound until the accumulated value and the total energy EyIs higher than the set retention η2The corresponding frequency is the lower bound omega of the frequencyminDetermining that 95% of energy is in a frequency interval [ omega ]minmax]。
Further, the total energy E of the intercepted signaly
Upper frequency bound omegamax
Lower frequency bound omegamin
Specifically, in step S5, cross-correlation is performed on the existing amplitude spectrum template function and the actually received data amplitude spectrum, and a normalized cross-correlation coefficient corr (v) corresponding to the cross-correlation function is obtained as follows:
wherein, | Y (ω) | is the amplitude spectrum of the vibration signal caused by the high-speed rail, d ω is the frequency infinitesimal, and the value range of the cross-correlation coefficient is between [0,1 ].
Specifically, in step S6, the corresponding speed vrComprises the following steps:
vr=argmaxvCorr(v)
wherein v isrIs a rough estimate of the velocity v.
Specifically, in step S7, [ v ]r-Δv1,vr+Δv1]Within range, the velocity search interval is reduced to Δ v2Reconstructing the amplitude spectrum template function | F (ω, v) | from the amplitude spectrum template function | F (ω, v) |, where v ═ v |r-Δv1+mΔv2
Specifically, in step S9, [ v ]r-Δv1,vr+Δv1]Within the range, the corresponding speed v determined according to step S6rFinding the speed corresponding to the maximum cross correlation coefficient, and recording the speed as the final estimated speed v of train operationfinal
Compared with the prior art, the invention has at least the following beneficial effects:
the invention relates to a method for estimating the running speed of a high-speed rail by using speed scanning, which can realize the estimation of the running speed of the high-speed rail train by only using geophone data outside an isolated area. The method comprises the steps of firstly generating a series of stress function amplitude spectrums for the structural parameters of the high-speed train and the running speed of a preset train, then calculating the cross-correlation function between the stress function amplitude spectrums at various speeds and the signal amplitude spectrum received by a single detector, and finally selecting the speed corresponding to the maximum cross-correlation function value as the estimated value of the running speed of the train. Compared with the conventional high-speed train speed estimation method, the method has the advantage that the running speed of the high-speed train can be obtained only by relying on one geophone data outside the isolation area.
Furthermore, the amplitude spectrum template function is designed according to the running speed and the typical parameters of the high-speed train, so that the subsequent calculation of the cross-correlation function of the amplitude spectrum template function and the actual signal is facilitated.
Furthermore, signals excited when high-speed rails pass through are intercepted from seismic data of the single detector, so that the speed of the train passing through the detector in each pass can be estimated, and the calculation amount of the subsequent calculation of the amplitude spectrum can be reduced.
Furthermore, the energy spectrum accumulation function of the intercepted signal and the energy frequency interval are determined, so that the frequency range participating in the cross-correlation function calculation is favorably reduced, and the noise resistance of the method is favorably improved.
Furthermore, the amplitude spectrum template function and the actual amplitude spectrum are subjected to cross-correlation calculation, and a data basis is provided for the subsequent search of the maximum correlation coefficient.
Further, the largest cross correlation coefficient and the corresponding speed thereof are found from all the cross correlation coefficients in step S5, which provides a coarser speed estimation, and is beneficial to improving the accuracy of the subsequent speed estimation.
Further, the velocity determined in step S6 is used as a base point, the velocity change interval is narrowed, the amplitude spectrum template function is reestablished, and the cross-correlation calculation between the amplitude spectrum template function and the actual amplitude spectrum is performed again, which is beneficial to improving the accuracy of subsequent velocity estimation.
Further, the speed corresponding to the maximum cross correlation coefficient is found and recorded as the final estimated speed of train operation, which is beneficial to obtaining the train speed estimation with higher precision.
In conclusion, the method can effectively and quickly realize the high-speed train running speed estimation by only using one geophone data, the adopted method is speed scanning, has high reliability and the like, and meanwhile, the method independent of equipment in a vehicle-mounted/isolation area is provided for detecting the high-speed train running speed.
The technical solution of the present invention is further described in detail by the accompanying drawings and embodiments.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 shows a seismic signal of a high-speed rail source received by a single detector when the train 1 passes by;
FIG. 3 is an amplitude spectrum of a seismic signal of a high-speed rail source received by a single detector when the train 1 passes by;
FIG. 4 is an amplitude spectrum template function generated at 300 km/h;
fig. 5 is a cross-correlation function of a series of amplitude spectrum template functions with the amplitude spectrum of an actual signal.
Detailed Description
The invention provides a high-speed rail running speed estimation method by using speed scanning, which can realize the estimation of the running speed of a high-speed rail train by only using geophone data outside an isolated area. The method comprises the steps of firstly generating a series of stress function amplitude spectrums for the structural parameters of the high-speed train and the running speed of a preset train, then calculating the cross-correlation function between the stress function amplitude spectrums at various speeds and the signal amplitude spectrum received by a single detector, and finally selecting the speed corresponding to the maximum cross-correlation function value as the estimated value of the running speed of the train.
Referring to fig. 1, the method for estimating the operating speed of a high-speed rail by using speed scanning according to the present invention includes the following steps:
s1, designing an amplitude spectrum template function according to the running speed and the typical parameters of the high-speed train;
the parameters of the high-speed train are relatively fixed, and typical train parameters are selected as follows: the number N of the carriages is 16, and the length D of each carriage body is 25 meters. At a certain speed interval deltav1Continuously varying the preset train speed v (which may vary from 30 km/h to 400 km/h, for example, at 1 km/h intervals) yields a series of speed-dependent amplitude spectrum template functions | F (ω, v) |:
v=vinitial+mΔv1 (2)
if v is to beinitialThe value of (a) is set to 30 km/h, and the value range of m is limited to [0,370 ]],Δv1At 1 km/h, a total of 371 functions are included in the series of amplitude spectrum template functions.
S2, intercepting the excited signal when the high-speed rail passes through in the seismic data of the single detector;
embedding a detector outside the high-speed rail circuit isolation area, intercepting a signal excited when the high-speed rail passes from signals received by the detector when the high-speed rail passes, wherein the time range corresponding to the obtained effective signal is [ t [ [ t ]1,t2]。
S3, carrying out Fourier transform on the intercepted signal to obtain an amplitude spectrum of the intercepted signal;
assuming that the signal caused by the intercepted high-speed rail operation is Y (t), performing Fourier transform on the signal to obtain an amplitude spectrum | Y (omega) | of the signal:
s4, calculating the accumulation function of the signal energy spectrum, and determining the frequency interval where most of the energy is located;
first, the total energy E of the signal is obtainedy
Then, the upper frequency bound of the signal is calculated, and the square of the amplitude is accumulated from the zero value of the frequency until the accumulated value is equal to the total energy EyIs higher than the set retention η1The corresponding frequency is the upper limit omega of the frequencymax
Similarly, when acquiring the lower frequency bound, the accumulation may be started from the upper frequency bound until the accumulated value and the total energy EyIs higher than the set retention η2The corresponding frequency is the lower bound omega of the frequencymin
Will eta1And η2Set to 0.975 and 0.95, respectively, thus determining that 95% of the energy is in the frequency interval [ omega ]minmax]。
S5, performing cross correlation on the amplitude spectrum template function and the actual amplitude spectrum;
in the frequency interval [ omega ]minmax]In the range, cross-correlation is carried out on a series of existing amplitude spectrum template functions and actually received data amplitude spectra, and a normalized cross-correlation coefficient Corr (v) corresponding to the cross-correlation coefficient is obtained and is as follows:
the value range of the cross correlation coefficient is between [0,1], and the closer to 1, the stronger the correlation between the amplitude spectrum template function under the current velocity v and the actual data amplitude spectrum is.
S6, finding the maximum cross correlation coefficient in all the cross correlation coefficients and the corresponding speed;
finding the largest cross-correlation coefficient among all cross-correlation coefficients and finding the corresponding velocity vrComprises the following steps:
vr=argmaxvCorr(v) (8)
wherein v isrIs a rough estimate of the velocity v.
S7, taking the speed determined in the step S6 as a base point, reducing the speed change interval, and reestablishing an amplitude spectrum template function;
velocity vrAfter determination is good, at [ v ]r-Δv1,vr+Δv1]Within range, the velocity search interval is reduced to Δ v2(which may take values of 0.01 km/h), reconstructing the amplitude spectrum template function | F (ω, v) | according to equation (1), wherein
v=vr-Δv1+mΔv2 (9)
Wherein,
s8, performing cross-correlation calculation of the amplitude spectrum template function and the actual amplitude spectrum again;
the cross-correlation coefficient corr (v) between the new template function and the actual signal is calculated using equation (7).
S9, finding the speed corresponding to the maximum cross correlation coefficient, and recording the speed as the final estimated speed v of train operationfinal
In [ v ]r-Δv1,vr+Δv1]Within the range, the speed corresponding to the maximum cross correlation coefficient is found by using the formula (8), and the speed is recorded as the final estimated speed v of train operationfinal
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some embodiments of the present invention, but not all embodiments. The components of the embodiments of the present invention generally described and illustrated in the figures herein may be arranged and designed in a wide variety of different configurations. Thus, the following detailed description of the embodiments of the present invention, presented in the figures, is not intended to limit the scope of the invention, as claimed, but is merely representative of selected embodiments of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without any inventive step, are within the scope of the present invention.
The invention relates to a high-speed rail running speed estimation method by using speed scanning, which takes a signal received by a single low-frequency detector at a distance of 75m from a high-speed rail passing time as an example.
Table 1 shows the running speed of the high-speed train estimated by the method when 8 trains pass
Referring to fig. 2, fig. 2 shows a vibration signal caused by a high-speed rail seismic source received by a single detector when the train 1 passes by, the sampling interval is 5ms, and 3001 sampling points are provided. Referring to fig. 3, 4 and 5, fig. 3 is an amplitude spectrum of a seismic signal caused by the passage of the train 1, fig. 4 is an amplitude spectrum template function generated at 300km/h, and fig. 5 is a cross-correlation function of a series of amplitude spectrum template functions with an actual signal amplitude spectrum. From fig. 5, it can be obtained that the maximum value of the cross-correlation coefficient corresponds to a velocity of 83.58 m/s (i.e., 300.89 km/h). The data received by a single detector when 8 trains pass by is analyzed by the method, and the obtained estimated speed value of the 8 trains is shown in table 1 and is consistent with the commercial operation speed of high-speed rails in China.
The above-mentioned contents are only for illustrating the technical idea of the present invention, and the protection scope of the present invention is not limited thereby, and any modification made on the basis of the technical solution according to the technical idea proposed by the present invention falls within the protection scope of the claims of the present invention.

Claims (10)

1. A high-speed rail running speed estimation method utilizing speed scanning is characterized by comprising the following steps:
s1, designing an amplitude spectrum template function | F (omega, v) | according to the running speed and typical parameters of the high-speed train;
s2, intercepting the signal excited by the passing of the high-speed rail in the seismic data of the single detector, wherein the time range corresponding to the effective signal acquired by the detector is [ t [ ]1,t2];
S3, carrying out Fourier transform on the intercepted signal to obtain an amplitude spectrum of the intercepted signal;
s4, calculating an energy spectrum accumulation function of the signal intercepted in the step S3, and determining an energy frequency interval;
s5, performing cross-correlation calculation on the amplitude spectrum template function and the actual amplitude spectrum;
s6, finding the largest cross correlation coefficient in all the cross correlation coefficients of the step S5 and the corresponding speed vr
S7, taking the speed determined in the step S6 as a base point, reducing the speed change interval, and reestablishing an amplitude spectrum template function;
s8, performing cross-correlation calculation of the amplitude spectrum template function and the actual amplitude spectrum again;
s9, finding the speed corresponding to the maximum cross correlation coefficient, and recording the speed as the final estimated speed v of the train operationfinal
2. The method of estimating speed of operation of a high-speed railway according to claim 1, wherein in step S1, the amplitude spectrum template function | F (ω, v) | is:
v=vinitial+mΔv1
wherein D is the length of each car body, v is the preset train speed, vinitialIs the initial scanning speed, m is an index of the scanning speed, Δ v1The initial velocity scan interval.
3. High speed rail operation speed estimation with speed scanning according to claim 2The method is characterized in that the number N of carriages is set to be 16, the length D of each carriage is 25 m, vinitialSet to 30 km/h, the value range of m is [0,370 ]],Δv1At 1 km/h, a total of 371 functions were included in a series of amplitude spectra template functions.
4. The method according to claim 1, wherein in step S3, assuming that the signal resulting from the intercepted high-speed rail operation is Y (t), the signal is fourier-transformed to obtain an amplitude spectrum | Y (ω) | of the signal:
wherein, [ t ]1,t2]And y (t) is a time range corresponding to the effective signal acquired by the detector, the intercepted signal caused by the running of the high-speed rail, j is an imaginary unit omega and is frequency, and t is speed.
5. The method of claim 1, wherein the total energy E of the intercepted signal is obtained first in step S4yThen, the upper frequency bound of the signal is calculated, and the square of the amplitude is accumulated from zero value of the frequency until the accumulated value is equal to the total energy EyIs higher than the set retention η1The corresponding frequency is the upper frequency bound omegamax(ii) a When acquiring the lower frequency bound, the accumulation is started from the upper frequency bound until the accumulated value and the total energy EyIs higher than the set retention η2When the corresponding frequency is the lower frequency bound omegaminDetermining that 95% of energy is in a frequency interval [ omega ]minmax]。
6. Method for estimating the speed of operation of a high-speed railway according to claim 5, characterized in that the total energy E of the intercepted signaly
Upper frequency bound omegamax
Lower frequency bound omegamin
7. The method as claimed in claim 1, wherein in step S5, the template function of the existing amplitude spectrum is cross-correlated with the amplitude spectrum of the data actually received, and the normalized cross-correlation coefficient corr (v) corresponding to the template function is obtained as follows:
wherein, | Y (ω) | is the amplitude spectrum of the vibration signal caused by the high-speed rail, d ω is the frequency infinitesimal, and the value range of the cross-correlation coefficient is between [0,1 ].
8. The method of claim 1, wherein the corresponding speed v is determined in step S6rComprises the following steps:
vr=argmaxvCorr(v)
wherein v isrIs a rough estimate of the velocity v.
9. The method of estimating a speed of operation of a high speed railway using a speed sweep as claimed in claim 1, wherein in step S7, [ v ] vr-Δv1,vr+Δv1]Within range, the velocity search interval is reduced to Δ v2Reconstructing according to the amplitude spectrum template function | F (ω, v) |Establishing an amplitude spectrum template function | F (ω, v) |, wherein v ═ v |r-Δv1+mΔv2
10. The method of estimating a speed of operation of a high speed railway using a speed sweep as claimed in claim 1, wherein in step S9, [ v ] vr-Δv1,vr+Δv1]Within the range, the corresponding speed v determined according to step S6rFinding the speed corresponding to the maximum cross correlation coefficient, and recording the speed as the final estimated speed v of train operationfinal
CN201910760699.4A 2019-08-16 2019-08-16 High-speed rail running speed estimation method utilizing speed scanning Active CN110514862B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910760699.4A CN110514862B (en) 2019-08-16 2019-08-16 High-speed rail running speed estimation method utilizing speed scanning

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910760699.4A CN110514862B (en) 2019-08-16 2019-08-16 High-speed rail running speed estimation method utilizing speed scanning

Publications (2)

Publication Number Publication Date
CN110514862A true CN110514862A (en) 2019-11-29
CN110514862B CN110514862B (en) 2020-10-27

Family

ID=68625499

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910760699.4A Active CN110514862B (en) 2019-08-16 2019-08-16 High-speed rail running speed estimation method utilizing speed scanning

Country Status (1)

Country Link
CN (1) CN110514862B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112230019A (en) * 2020-10-10 2021-01-15 西安交通大学 High-speed rail train running acceleration estimation method using multiple geophones
CN113341175A (en) * 2021-06-04 2021-09-03 西安交通大学 High-speed rail running acceleration estimation method and system based on single detector

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101858924A (en) * 2010-05-14 2010-10-13 北京理工大学 Method and device for measuring vehicle speed by utilizing weak magnetic signal correlation analysis
CN102053166A (en) * 2009-10-07 2011-05-11 罗伯特.博世有限公司 Method for determining speed
US8711337B2 (en) * 2008-12-18 2014-04-29 University Of Zurich Passive translational velocity measurement from optical information
JP2015108563A (en) * 2013-12-05 2015-06-11 公益財団法人鉄道総合技術研究所 Train speed measurement method and train position specifying method, and device thereof
CN108313089A (en) * 2017-01-18 2018-07-24 扬州立鼎恒新微电子科技有限公司 A kind of train real-time location method based on MEMS shock sensors
CN108596946A (en) * 2018-03-21 2018-09-28 中国航空工业集团公司洛阳电光设备研究所 A kind of moving target real-time detection method and system
CN108802811A (en) * 2017-04-28 2018-11-13 中国石油化工股份有限公司 A kind of normal-moveout spectrum automatic pick method and device
CN108959176A (en) * 2018-06-20 2018-12-07 西南交通大学 A kind of distributed vehicle speed estimation method based on adaptive volume Kalman filtering

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8711337B2 (en) * 2008-12-18 2014-04-29 University Of Zurich Passive translational velocity measurement from optical information
CN102053166A (en) * 2009-10-07 2011-05-11 罗伯特.博世有限公司 Method for determining speed
CN101858924A (en) * 2010-05-14 2010-10-13 北京理工大学 Method and device for measuring vehicle speed by utilizing weak magnetic signal correlation analysis
JP2015108563A (en) * 2013-12-05 2015-06-11 公益財団法人鉄道総合技術研究所 Train speed measurement method and train position specifying method, and device thereof
CN108313089A (en) * 2017-01-18 2018-07-24 扬州立鼎恒新微电子科技有限公司 A kind of train real-time location method based on MEMS shock sensors
CN108802811A (en) * 2017-04-28 2018-11-13 中国石油化工股份有限公司 A kind of normal-moveout spectrum automatic pick method and device
CN108596946A (en) * 2018-03-21 2018-09-28 中国航空工业集团公司洛阳电光设备研究所 A kind of moving target real-time detection method and system
CN108959176A (en) * 2018-06-20 2018-12-07 西南交通大学 A kind of distributed vehicle speed estimation method based on adaptive volume Kalman filtering

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王晓凯 等: "高铁震源地震信号的挤压时频分析应用", 《地球物理学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112230019A (en) * 2020-10-10 2021-01-15 西安交通大学 High-speed rail train running acceleration estimation method using multiple geophones
CN113341175A (en) * 2021-06-04 2021-09-03 西安交通大学 High-speed rail running acceleration estimation method and system based on single detector
CN113341175B (en) * 2021-06-04 2022-10-25 西安交通大学 High-speed rail running acceleration estimation method and system based on single detector

Also Published As

Publication number Publication date
CN110514862B (en) 2020-10-27

Similar Documents

Publication Publication Date Title
CN110658005B (en) Method for identifying rail corrugation diseases based on vehicle body acceleration
CN110514862B (en) High-speed rail running speed estimation method utilizing speed scanning
Lindveld et al. Evaluation of online travel time estimators and predictors
CN105424388B (en) A kind of Railway wheelset bearing fault transient feature detection method based on parametrization Doppler's transient model
DE102014226612B4 (en) Method for determining the speed of a rail vehicle
CN110514861B (en) High-speed rail running speed estimation method utilizing amplitude spectrum autocorrelation
CN109116048A (en) A kind of rotating speed of automobile engine detection and control method
López-Valcarce et al. Estimation of road vehicle speed using two omnidirectional microphones: A maximum likelihood approach
CN110530509B (en) High-speed motor train unit axle box vibration dominant frequency prediction method based on maximum entropy spectrum analysis
CN110285877A (en) Train real time location tracking and speed calculation method based on Spark Streaming
CN102914364B (en) Dynamic identifying method for instant contact loss of wheel rail under impact effect of wheel rail
CN109696609A (en) A kind of acquisition methods and its acquisition system of novel bow net current carrying quality
IWATA et al. Improved algorithms of seismic parameters estimation and noise discrimination in earthquake early warning
CN115758289B (en) Rail wave mill identification method based on multitask learning neural network
DE102012219111A1 (en) Method for localization of rail vehicle in rail network, involves matching corresponding period of reference signal of route section of rail network with selected portion of actual time signal
CN106404159B (en) A kind of determining method of Rocket sled test continuous shaking transmitting spectrum
CN115905811A (en) Railway track irregularity detection method based on FIR-wavelet transform
Lu et al. Information fusion and image processing based arc detection and localization in pantograph-catenary systems
CN108528475B (en) A kind of track transition fault alarm method based on multi-level fusion
Toyoda et al. Traffic monitoring with ad-hoc microphone array
Kojima et al. Fault detection of railway track by multi-resolution analysis
CN104977357B (en) Steel rail crack sound emission signal extracting and denoising method based on kalman filtering
CN106650198A (en) Rail transit vehicle velocity calculating method and system based on information fusion
CN112230019A (en) High-speed rail train running acceleration estimation method using multiple geophones
CN113341175B (en) High-speed rail running acceleration estimation method and system based on single detector

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