阅读说明：本技术 一种利用振幅谱自相关的高铁运行速度估计方法 (It is a kind of to utilize the autocorrelative high-speed rail speed of service estimation method of amplitude spectrum ) 是由 王晓凯 陈文超 师振盛 于 2019-08-16 设计创作，主要内容包括：本发明公开了一种利用振幅谱自相关的高铁运行速度估计方法,在单检波器地震数据中截取高铁经过时所激发的信号；对截取的信号做傅里叶变换,得到其振幅谱；计算前一步骤截取的信号的能量谱累积函数,确定能量频率区间；计算信号振幅谱的自相关函数；在自相关函数中寻找最大峰值、第二峰值以及其对应的频率,进而获得最大峰值与第二峰值之间的频率间隔；利用最大峰值与第二峰值之间的频率间隔及高铁列车单个车厢长度,获得高铁列车运行速度的估计。(The autocorrelative high-speed rail speed of service estimation method of amplitude spectrum is utilized the invention discloses a kind of, 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 energy spectrum cumulative function for calculating the signal of previous step interception, determines energy frequency section；Calculate the auto-correlation function of signal amplitude spectrum；Peak-peak, the second peak value and its corresponding frequency are found in auto-correlation function, and then obtain the frequency interval between peak-peak and the second peak value；Using the frequency interval and the single car length of high-speed rail train between peak-peak and the second peak value, the estimation of high-speed rail train running speed is obtained.)

1. a kind of utilize the autocorrelative high-speed rail speed of service estimation method of amplitude spectrum, which comprises the following steps:

S1, the signal excited when high-speed rail is passed through is intercepted in single detector seismic data；

S2, Fourier transformation is done to the signal of interception, obtains its amplitude spectrum；

S3, the energy spectrum cumulative function for calculating step S2 intercept signal, determine energy frequency section；

S4, the auto-correlation function for calculating signal amplitude spectrum；

S5, peak-peak, the second peak value and its corresponding frequency are found in auto-correlation function, so obtain peak-peak with Frequency interval between second peak value；

S6, the frequency interval and the single car length of high-speed rail train between peak-peak and the second peak value, acquisition high-speed rail column are utilized The estimated value v of the vehicle speed of service_final。

2. according to claim 1 utilize the autocorrelative high-speed rail speed of service estimation method of amplitude spectrum, which is characterized in that step In rapid S2, it is assumed that the caused signal of the high-speed rail operation intercepted out is y (t), does Fourier transformation to the signal, obtains its amplitude Spectrum | Y (ω) | are as follows:

Wherein, [t₁,t₂] it is time range corresponding to the useful signal of wave detector acquisition, y (t) is the high-speed rail operation intercepted out Caused signal, j are that imaginary unit ω is frequency, and t is speed.

3. according to claim 1 utilize the autocorrelative high-speed rail speed of service estimation method of amplitude spectrum, which is characterized in that step Rapid S3 specifically:

The gross energy E of intercept signal is obtained first_y, the upper frequency range of intercept signal is then calculated, to amplitude square by frequency Zero starts to add up, until accumulated value and gross energy E_yRatio be higher than setting retention rate η₁, corresponding frequency is upper frequency range ω_max；It when obtaining lower frequency border, adds up by upper frequency range, until accumulated value and gross energy E_yRatio be higher than setting Retention rate η₂, corresponding frequency is lower frequency border ω at this time_min, the frequency separation where determining 95% energy is [ω_min, ω_max]。

4. according to claim 3 utilize the autocorrelative high-speed rail speed of service estimation method of amplitude spectrum, which is characterized in that cut The signal gross energy E taken_yAre as follows:

Upper frequency range ω_maxAre as follows:

Lower frequency border ω_minAre as follows:

5. according to claim 1 utilize the autocorrelative high-speed rail speed of service estimation method of amplitude spectrum, which is characterized in that step In rapid S4, in frequency separation [ω_min,ω_max] in range, calculate signal amplitude spectrum | Y (ω) | auto-correlation function Corr (u) Are as follows:

Wherein, the value range of auto-correlation function Corr (u) is [0,1].

6. according to claim 1 utilize the autocorrelative high-speed rail speed of service estimation method of amplitude spectrum, which is characterized in that step In rapid S5, peak-peak, the second peak value and its corresponding frequency are found in auto-correlation function, obtain peak-peak and second Frequency interval u between peak value_interval。

7. according to claim 1 utilize the autocorrelative high-speed rail speed of service estimation method of amplitude spectrum, which is characterized in that step In rapid S6, high-speed rail train running speed v_finalAre as follows:

v_final=Lu_interval

Wherein, L is the length of single-unit high-speed rail compartment, u_intervalFor the frequency interval between peak-peak and the second peak value.

Technical field

The invention belongs to exploration geophysics fields, in particular to a kind of to utilize the autocorrelative high-speed rail speed of service of amplitude spectrum Estimation method.

Background technique

So far, Chinese high-speed rail revenue kilometres are up to 3.1 ten thousand kilometers, close to the 70% of the total revenue kilometres of world's high-speed rail. From when 10 days zero July in 2019, china railway implements new route map of train, and starting motor train unit train daily is more than 3000 pairs. High-speed rail train high-speed cruising is on high-speed rail route, and the speed of service is the important parameter for not reflecting train operating safety still, together When be also key that later use high-speed rail train causes vibration signal.The method for obtaining high-speed rail train running speed existing at present Mainly include:

Summary of the invention

In view of the above-mentioned deficiencies in the prior art, the technical problem to be solved by the present invention is that providing a kind of using amplitude Autocorrelative high-speed rail speed of service estimation method is composed, the collected earthquake number of the outer single detector of high-speed rail route isolated area is used only According to estimating train running speed using the amplitude spectrum auto-correlation for collecting data, provided for subsequent judgement train operation state Data.

The invention adopts the following technical scheme:

It is a kind of to utilize the autocorrelative high-speed rail speed of service estimation method of amplitude spectrum, comprising the following steps:

S1, the signal excited when high-speed rail is passed through is intercepted in single detector seismic data；

S2, Fourier transformation is done to the signal of interception, obtains its amplitude spectrum；

S3, the energy spectrum cumulative function for calculating step S2 intercept signal, determine energy frequency section；

S4, the auto-correlation function for calculating signal amplitude spectrum；

S5, peak-peak, the second peak value and its corresponding frequency are found in auto-correlation function, and then obtain maximum peak Frequency interval between value and the second peak value；

S6, the frequency interval and the single car length of high-speed rail train between peak-peak and the second peak value, acquisition height are utilized The estimated value v of iron train running speed_final。

Specifically, in step S2, it is assumed that the caused signal of the high-speed rail operation intercepted out is y (t), is in Fu to the signal Leaf transformation obtains its amplitude spectrum | Y (ω) | are as follows:

Wherein, [t₁, t₂] it is time range corresponding to the useful signal of wave detector acquisition, y (t) is the high-speed rail intercepted out The caused signal of operation, j is that imaginary unit ω is frequency, and t is speed.

Specifically, step S3 specifically:

The gross energy E of intercept signal is obtained first_y, then calculate intercept signal upper frequency range, to amplitude square by The zero of frequency starts to add up, until accumulated value and gross energy E_yRatio be higher than setting retention rate η₁, corresponding frequency is frequency Rate upper bound ω_max；It when obtaining lower frequency border, adds up by upper frequency range, until accumulated value and gross energy E_yRatio be higher than The retention rate η of setting₂, corresponding frequency is lower frequency border ω at this time_min, the frequency separation where determining 95% energy is [ω_min, ω_max]。

Further, the signal gross energy E of interception_yAre as follows:

Upper frequency range ω_maxAre as follows:

Lower frequency border ω_minAre as follows:

Specifically, in step S4, in frequency separation [ω_min, ω_max] in range, calculate signal amplitude spectrum | Y (ω) | from Correlation function Corr (u) are as follows:

Wherein, the value range of auto-correlation function Corr (u) is [0,1].

Specifically, peak-peak, the second peak value and its corresponding frequency are found in auto-correlation function in step S5, Obtain the frequency interval u between peak-peak and the second peak value_interval。

Specifically, in step S6, high-speed rail train running speed v_finalAre as follows:

v_final=Lu_interval

Wherein, L is the length of single-unit high-speed rail compartment, u_intervalBetween frequency between peak-peak and the second peak value Every.

Compared with prior art, the present invention at least has the advantages that

The present invention be using amplitude spectrum auto-correlation carry out the estimation of the high-speed rail speed of service a kind of method, merely with one every The estimation of high-speed rail train running speed can be realized in geophone data outside from area.The present invention calculates wave detector first and is connect The signal amplitude spectrum that the high-speed rail train operation received is excited, then calculates the auto-correlation function of the amplitude spectrum, then using certainly The frequency interval between peak-peak and the second largest peak value in correlation function finally estimates high-speed rail train using the frequency interval The speed of service.It is compared to conventional high-speed rail train speed estimation method, present invention depends only upon the earthquake inspections outside isolated area Wave device data are the speed of service for being convenient to obtain high-speed rail train.

Further, the signal excited when high-speed rail is passed through is intercepted in single detector seismic data, and it is every to be conducive to estimation The speed by wave detector train is plowed, while advantageously reducing the operand of subsequent calculated amplitude spectrum.

Further, the energy spectrum cumulative function of the signal of interception and determining energy frequency section, advantageously reduce participation The frequency range that cross-correlation function calculates, while being conducive to the noise immunity of improvement method.

Further, the auto-correlation function for calculating signal amplitude spectrum, is conducive to subsequent determining peak-peak and the second peak value Interval；.

Further, peak-peak, the second peak value and its corresponding frequency are found in auto-correlation function, and then are obtained Frequency interval between peak-peak and the second peak value can get variation the closest in train running speed and signal spectrum Amount.

Further, using the frequency interval and the single car length of high-speed rail train between peak-peak and the second peak value, The estimation for obtaining high-speed rail train running speed, help to obtain the accurate estimation of high-speed rail train running speed.

In conclusion the present invention can effectively and quickly realize the high-speed rail train merely with a geophone data Speed of service estimation, used method are that calculated amplitude composes auto-correlation function, have the characteristics that high reliablity, while being detection High-speed rail train running speed provides a method independently of equipment in vehicle-mounted/isolated area.

Below by drawings and examples, technical scheme of the present invention will be described in further detail.

Detailed description of the invention

Fig. 1 is flow chart of the present invention；

Fig. 2 is one of high-speed rail focus seismic signal that single detector receives when train 1 passes through；

Fig. 3 is the amplitude spectrum for one of high-speed rail focus seismic signal that single detector receives when train 1 passes through；

Fig. 4 is the auto-correlation function of the spectrum of signal amplitude shown in Fig. 3.

Specific embodiment

Autocorrelative high-speed rail speed of service estimation method is composed using signal amplitude the present invention provides a kind of, merely with one The estimation of high-speed rail train running speed can be realized in geophone data outside isolated area.Wave detector is calculated first to be received To high-speed rail train operation excited signal amplitude spectrum, then calculate the auto-correlation function of the amplitude spectrum, then use from phase The peak-peak in function and the frequency interval between the second largest peak value are closed, is finally transported using frequency interval estimation high-speed rail train Scanning frequency degree.It is compared to conventional high-speed rail train speed estimation method, present invention depends only upon a seismic detections outside isolated area Device data are the speed of service for being convenient to obtain high-speed rail train.

Referring to Fig. 1, the present invention is a kind of using the autocorrelative high-speed rail speed of service estimation method of signal amplitude spectrum, including with Lower step:

S1, the signal excited when high-speed rail is passed through is intercepted in single detector seismic data；

The embedding detector outside high-speed rail route isolated area is intercepted from the signal that wave detector receives high when high-speed rail is passed through The signal excited when iron passes through, time range corresponding to the useful signal of acquisition are [t₁, t₂]。

S2, Fourier transformation is done to the signal of interception, obtains its amplitude spectrum；

Assuming that signal caused by the high-speed rail operation intercepted out is y (t), Fourier transformation is done to the signal, obtains its vibration Width spectrum | Y (ω) | are as follows:

S3, the cumulative function for calculating this signal energy spectrum determine the frequency separation where its most of energy；

The gross energy E of the signal is obtained first_y:

Then the upper frequency range for calculating the signal, to square adding up by the zero of frequency for amplitude, until accumulated value With gross energy E_yRatio be higher than setting retention rate η₁, corresponding frequency is required upper frequency range ω at this time_max:

Similarly, it when obtaining lower frequency border, can add up by upper frequency range, until accumulated value and gross energy E_yRatio Value is higher than the retention rate η of setting₂, corresponding frequency is required lower frequency border ω at this time_min:

By η₁And η₂It is set to 0.975 and 0.95, the frequency separation where thereby determining that 95% energy is [ω_min, ω_max]。

S4, the auto-correlation function for calculating signal amplitude spectrum；

In frequency separation [ω_min, ω_max] in range, calculate signal amplitude spectrum | Y (ω) | auto-correlation function Corr (u) Are as follows:

The value range of the auto-correlation function is between [0,1].

S5, peak-peak, the second peak value and its corresponding frequency are found in auto-correlation function, and then obtain maximum peak Frequency interval between value and the second peak value；

According to the characteristic of auto-correlation function, the position of peak-peak is sought in auto-correlation function generally in zero frequency The second peak value is looked for, recording its position is u_interval, then the frequency interval between peak-peak and the second peak value is u_interval。

S6, the frequency interval u between peak-peak and the second peak value is utilized_intervalAnd the single car length of high-speed rail train, Obtain the estimation v of high-speed rail train running speed_final。

High-speed rail compartment be usually L meters (Chinese high-speed rail compartment length is usually 25 meters) and peak-peak with Frequency interval u between second peak value_intervalIt can estimate high-speed rail train running speed v_final:

v_final=Lu_interval (6)

Wherein, L is the length of single-unit high-speed rail compartment, and unit is rice.

In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is A part of the embodiment of the present invention, instead of all the embodiments.The present invention being described and shown in usually here in attached drawing is real The component for applying example can be arranged and be designed by a variety of different configurations.Therefore, below to the present invention provided in the accompanying drawings The detailed description of embodiment be not intended to limit the range of claimed invention, but be merely representative of of the invention selected Embodiment.Based on the embodiments of the present invention, those of ordinary skill in the art are obtained without creative efforts The every other embodiment obtained, shall fall within the protection scope of the present invention.

A kind of high-speed rail speed of service estimation method using velocity scanning of the present invention, away from high-speed rail route when being passed through with high-speed rail For signal received by the single low-frequency detector of 75m.

Table 1 is the high-speed rail train running speed estimated when 8 trains pass through using this method

Referring to Fig. 2, Fig. 2 is vibration letter caused by high-speed rail focus received by single detector when train 1 passes through Number, sampling interval 5ms shares 3001 sampled points.Fig. 3 and Fig. 4 are please referred to, Fig. 3 is caused vibration letter when train 1 passes through Number amplitude spectrum, Fig. 4 is the amplitude spectrum auto-correlation function obtained using Fig. 3, between available peak-peak and the second peak value Distance be 3.3396 hertz, when car length is 25 meters, the speed of service of available this train be 83.49 meter per seconds (i.e. 300.56 thousand ms/h).Data received by single detector are divided using this patent method when passing through to 8 trains Analysis, 8 times obtained train speed estimated values are as shown in table 1, coincide with the commercial operation speed of China's high-speed rail.

The above content is merely illustrative of the invention's technical idea, and this does not limit the scope of protection of the present invention, all to press According to technical idea proposed by the present invention, any changes made on the basis of the technical scheme each falls within claims of the present invention Protection scope within.