CN102901994A - Verification system and method for measuring sensor by using seismograph - Google Patents

Abstract

The present invention relates to a verification system and method for measuring a sensor by using a seismograph, which perform time zone bandpass filter and hilbert transform to seismological observation measurement data of a novel seismological observation instrument and an aging seismometer or shaking table testing data so as to raise a response spectrum amplitude and phase determination efficiency of the seismograph. According to the invention, the response spectrum of a seismic sensor is obtained through time zone bandpass filter and hilbert transform to the seismological observation measurement data of a novel seismological observation instrument and an aging seismometer or the shaking table testing data, thereby providing the effects of raising the amplitude and phase determination accuracy. Moreover, characteristics of the sensor can be provided irrelevant to the kind of the sensor to be measured, thereby improving the inconvenience of requirement of introducing sensors of the same kind.

Description

Measurement sensor verification system and method using seismograph

technical field

The present invention relates to a measurement sensor verification system and method using a seismograph, and in more detail, relates to a time zone band-pass filter for seismic observation measurement data or shaking table test data of novel seismic observation instruments and aging seismometers A measurement sensor verification system and method using a seismograph to improve the response spectrum amplitude and phase measurement accuracy of a seismograph using Bandpass filter and Hilbert transform.

Background technique

The earthquake monitoring system is a device that measures the voltage of a seismometer sensor (velocity, acceleration) that mainly uses a seismograph to measure ground movement. When purchasing a new seismograph, it is necessary to confirm the performance of the seismometer and sensor that the manufacturer suggested. There may be a need to verify the performance of older, aging seismic monitoring systems located at seismic observatories.

As a verification method of a seismograph, a seismograph sensor is installed on a vibrating table (vibration generating device), and the vibration of the seismograph sensor to be measured is measured by exciting (generating vibration) with reference to the absolute reference sensor.

At this time, there is a method of simply performing Fourier transform on the measured vibration value to transform it into the frequency domain, and then divide it by the absolute reference sensor measurement value to measure it.

In this case, there is a drawback of mastering performance and proven absolute seismometer performance.

In addition, when the excitation vibration value is very small, measurement errors may occur due to a lot of noise measurement.

Another method, as a method for verifying the performance of a seismometer, is a method in which two sensors are installed in an actual seismic observatory, and the frequencies of the sensors measured from the sensor data (information) having the measured response spectrum of the seismometer sensor are performed. Conversion, a method of measuring using a spectral ratio.

In the case of such a measurement, the part having the noise of the sensor itself contains a lot of noise in the part of the response spectrum of the measurement sensor.

Also, if only the above-mentioned frequency conversion is performed, there is only one measurement value for each frequency, so the error in the response spectrum can only be amplified.

Therefore, in terms of obtaining the seismometer sensor response spectrum of the vibration data measured by the seismometer, the spectrum of components lower than the noise of the seismometer sensor itself produces a portion that distorts the sensor response spectrum, and a method that can suppress the sensor response spectrum by suppressing A method to improve the accuracy of the vibration data spectrum of the seismometer sensor's own noise.

In addition, when carrying out seismic research, it is necessary to obtain precise seismic data, to correct the seismograph, or to lose the measurable information of the seismograph. In this case, a method for correcting the seismograph is required. .

As an example of the background art, it is disclosed in "Seismic Wave Detection System" as Korean Patent No. 0594625 (published on June 30, 2006).

After the seismic wave signal is input and filtered in the above background technology, the ratio of the short time average (Short Time average) to the long time average (Long Timeaverage) of the seismic wave signal is exchanged for the interval value of the specified time window (TimeWindow), and the determination is based on the verification Differentiate the measured analog signal, and control the substitution, storage, recording control, calculation, etc. of seismic wave variables, so that they can be used for differentiation. slope to control any two intervals.

However, the above-mentioned background technology cannot distinguish background noise and seismic wave amplitude (Amplitude) when generating micro earthquakes (Micro Earthquake), so that the performance of seismic wave detection is significantly reduced, so it has the disadvantage that it is impossible to accurately determine seismic factors. .

prior art literature

patent documents

(Patent Document 1) Korea Authorized Patent 0594625 (2006.06.30)

Contents of the invention

Therefore, the present invention is made in consideration of the problems of the prior art as described above. Bandpass filter and Hilbert transform are used to obtain the response spectrum of the seismic sensor, thereby improving the accuracy of amplitude and phase measurement.

Another object of the present invention is to provide sensor characteristics regardless of the type of sensor to be measured.

Another object of the present invention is that when two sensors are installed in a general seismic observation station, it is not necessary to bring the sensors to the laboratory shaking table, and the performance confirmation and verification of the sensors can be performed using the data obtained in the seismic observation.

Still another object of the present invention is to perform seismograph correction only by changing the channel of the seismograph regardless of the excitation source, so that not only the shaking table test but also the general seismic observation station can only change the channel to easily perform corrections.

In order to achieve the purpose of the problem to be solved by the present invention, a measurement sensor verification system utilizing a seismograph according to an embodiment of the present invention includes: a vibration excitation device 200, which excites a vibration table; a vibration table 300, which is configured to Have reference sensor 400 and measurement sensor 45, be connected with above-mentioned excitation device, and exciter vibration when starting vibration device; Signal line converter 600, it is used to measure the vibration signal of above-mentioned reference sensor and measurement sensor; Signal amplifier 700 , which amplifies the vibration signal transformed in the above-mentioned signal line converter, and transmits the amplified vibration signal to the seismic recorder; the seismic recorder 800, which is used to measure the amplified vibration signal transmitted from the above-mentioned signal amplifier; seismic data processing Mechanism 100, which calculates the response spectrum of the reference sensor and the measurement sensor measured by the above-mentioned seismograph, transforms the vibration signal of the reference sensor into the frequency domain, and corrects it into a response spectrum, then compares the above-mentioned calculated response spectrum with the corrected The response spectra are compared to judge whether the measuring sensor is normal or not. This solves the problem of the invention.

The measurement sensor verification system and method utilizing the seismograph of the present invention with the above structures and effects performs a time zone band-pass filter (Bandpass filter) on the seismic observation measurement data or shaking table test data of novel seismic observation instruments and aging seismometers. filter) and Hilbert transform to obtain the response spectrum of the seismic sensor, thereby providing the effect of improving the accuracy of amplitude and phase measurement.

Furthermore, since the characteristics of the sensor are provided regardless of the type of sensor to be measured, the inconvenience of having to introduce the same type of sensor can be improved.

In addition, when two sensors are installed in a general seismic observation station, it is not necessary to bring the sensors to the laboratory shaking table, and the performance confirmation and verification of the sensors can be performed using the data obtained in the seismic observation, thereby providing convenience in use.

In addition, it is possible to perform seismograph correction by simply changing the channel of the seismograph regardless of the excitation source, so that not only in shaking table tests but also in general seismic observation stations, only the channel can be changed. Perform correction work easily.

Description of drawings

FIG. 1 is a configuration diagram schematically constituting a measurement sensor verification system using a seismograph according to an embodiment of the present invention.

Fig. 2 is a block diagram of a seismic data processing mechanism of a measurement sensor verification system using a seismograph according to an embodiment of the present invention.

FIG. 3 is a block diagram of a response spectrum calculation unit of a measurement sensor verification system using a seismograph according to an embodiment of the present invention.

FIG. 4 is a flowchart showing a measurement sensor verification method using a seismograph according to an embodiment of the present invention.

FIG. 5 is a flowchart showing a processing procedure of a response spectrum calculation unit of a measurement sensor verification system using a seismograph according to an embodiment of the present invention.

6 to 7 are diagrams showing a test execution procedure diagram of a seismograph for correcting a measurement sensor verification system using a seismograph according to an embodiment of the present invention.

8 is a block diagram of a seismograph for modifying a measurement sensor verification system using a seismograph according to an embodiment of the present invention.

9 is a flowchart showing a seismograph correction method using the measurement sensor verification system for a seismograph according to an embodiment of the present invention.

Explanation of reference signs

100: Seismic data processing organization

200: Vibration device

300: Shaker

400: Reference sensor

450: Measuring sensor

600: signal line converter

700: signal amplifier

800: Seismic Recorder

Detailed ways

The measurement sensor verification system utilizing the seismic recorder according to an embodiment of the present invention for achieving the above object is characterized in that it includes: a vibration excitation device 200, which excites a vibration table; a vibration table 300, which is configured to Have reference sensor 400 and measurement sensor 45, be connected with above-mentioned excitation device, and exciter vibration when starting vibration device; Signal line converter 600, it is used to measure the vibration signal of above-mentioned reference sensor and measurement sensor; Signal amplifier 700 , which amplifies the vibration signal transformed in the above-mentioned signal line converter, and transmits the amplified vibration signal to the seismic recorder; the seismic recorder 800, which is used to measure the amplified vibration signal transmitted from the above-mentioned signal amplifier; seismic data processing Mechanism 100, which calculates the response spectrum of the reference sensor and the measurement sensor measured by the above-mentioned seismograph, transforms the vibration signal of the reference sensor into the frequency domain, and corrects it into a response spectrum, then compares the above-mentioned calculated response spectrum with the corrected The response spectra are compared to judge whether the measuring sensor is normal or not.

At this time, the above-mentioned seismic data processing mechanism 100 is characterized in that it includes: a reference sensor frequency domain conversion unit 110, which obtains the vibration signal of the reference sensor, and transforms it into the frequency domain by Fourier transform; a response spectrum correction unit 120, which uses The response spectrum is corrected by dividing the frequency spectrum converted by the above-mentioned reference sensor frequency domain conversion part by the theoretical response spectrum or the known response spectrum; the same type sensor determination part 130 is used to determine whether the reference sensor and the measurement sensor are of the same type The sensor; the unit conversion unit 140, which is used to convert the measurement sensor in the frequency domain using a differential operator or an integral operator so as to have the same measurement unit when the reference sensor and the measurement sensor are not the same type of sensor The response spectrum calculation part 150, which is used to calculate the response spectrum by bandpass filter filtering and Hilbert transform when the reference sensor and the measurement sensor are the same type of sensor; the measurement sensor is normal or not judgment part 160, which It is used to compare the response spectrum calculated by the above-mentioned response spectrum calculation part with the response spectrum corrected by the response spectrum correction part to judge whether the measurement sensor is normal or not; the result value output part 170 is used to output whether the above-mentioned measurement sensor is normal or not. The value of the result of judgment by the no judgment unit.

At this time, the above-mentioned response spectrum calculation part 150 is characterized in that it includes: a real number part and an imaginary part setting part 151, which uses the measurement sensor information and the reference sensor information to perform time-domain band-pass filter filtering in the set band-pass frequency band. , and after performing the Hilbert transform, the bandpass filter filtering information is set as the real number part, and the Hilbert transform information is set as the imaginary number part; the amplitude and phase calculation part 152 is used to total the real number part and the After the real number part and the imaginary number part set by the imaginary number part setting part, the amplitude and phase that change with time are calculated; the amplitude extraction part 153 is used for the time series of the reference sensor and the measurement sensor calculated by the above-mentioned amplitude and phase calculation part The amplitude is used to extract only the amplitude value higher than the noise of the reference sensor or the measurement sensor itself; the response amplitude spectrum ratio calculation part 154 is for the amplitude extracted by the above-mentioned amplitude extraction part, by summing up the time-series amplitude values of the reference sensor and the measurement sensor Calculate the amplitude ratio (measurement sensor/reference sensor) to calculate the response amplitude spectrum; the response phase spectrum calculation part 155 is used to calculate the phase difference between the reference sensor and the measurement sensor (measurement sensor phase-reference sensor phase), and calculate the phase difference The mean and standard deviation values, in the phase difference value, are calculated in the mean value for the mean value of the phase difference values within the standard deviation range to calculate the response phase spectrum.

On the other hand, the measurement sensor verification method using a seismograph according to the present invention is characterized in that it includes the following steps: a vibrating device vibration scanning step (S100) in which the vibrating table is vibrated by activating the vibrating device; The processing mechanism obtains the experimental information acquisition step (S200) of the reference sensor and the vibration signal of the measurement sensor measured by the seismograph; Reference sensor frequency domain conversion step (S300); the response spectrum correction unit 120 divides the frequency domain transformed by the reference sensor frequency domain conversion unit by the theoretical response spectrum correction of the reference sensor or the known response spectrum to correct the response of the response spectrum Spectrum correction step (S400); same-type sensor determination step (S500) in which the same-type sensor determination unit 130 determines whether the reference sensor and the measurement sensor are the same type of sensor; if the reference sensor and the measurement sensor are the same type of sensor, Response spectrum calculation part 150 utilizes band-pass filter filtering and Hilbert transform to calculate the response spectrum calculation step (S700) of response spectrum; The response spectrum corrected by the response spectrum correction unit is compared to determine whether the measurement sensor is normal or not. The measurement sensor is normal or not judging step (S800); the result value output unit 170 outputs the result of the result value judged by the measurement sensor normal or not judging unit. Value output step.

In this case, it is characterized in that, when the reference sensor and the measurement sensor are not the same type of sensor as a result of the determination in the above-mentioned same-type sensor determination step (S500), the unit conversion unit 140 performs a differential operation in the frequency domain using the measurement sensor as a reference. A unit conversion step (S600) of converting a symbol or an integral operator into the same measurement unit.

At this time, the above-mentioned response spectrum calculation step (S700) is characterized in that it includes the following steps: the real number part and imaginary number part setting part 151 performs time-domain band-pass filtering in the set band-pass frequency band using the measurement sensor information and the reference sensor information filter, and after performing the Hilbert transform, set the bandpass filter filtering information as the real number part, and set the Hilbert transform information as the imaginary part of the real number part and the imaginary part setting step (S710); and the phase calculation part 152 calculates the amplitude and phase calculation step (S720) of the amplitude and phase that change with time after the real number part and the imaginary number part set by the real number part and the imaginary number part setting part are totaled (S720); the amplitude extraction part 153 For the time-series amplitudes of the reference sensor and the measurement sensor calculated by means of the amplitude and phase calculation unit, only the amplitude extraction step (S730) of extracting the amplitude value higher than the noise of the reference sensor or the measurement sensor itself; the response amplitude spectrum ratio calculation unit 154 Response amplitude spectrum ratio calculation step (S740) for calculating the response amplitude spectrum ratio by summing the time-series amplitude values of the reference sensor and the measurement sensor to calculate the amplitude ratio (measurement sensor/reference sensor) for the amplitude extracted by the amplitude extraction unit (S740); The phase spectrum calculation unit 155 calculates the phase difference between the reference sensor and the measurement sensor (measurement sensor phase-reference sensor phase), and calculates the average value and standard deviation value of the phase difference. A response phase spectrum calculation step (S750) for calculating the response phase spectrum by averaging the phase difference values within the range.

On the other hand, the above-mentioned seismic data processing mechanism according to yet another embodiment of the present invention is characterized in that it further includes: a primary response spectrum ratio calculation unit 180, which connects the channel one based on the first excitation The frequency domain transformation of the vibration signal of the reference sensor and the measured value of the response spectrum are divided by the measured value of the frequency domain transformation and the response spectrum of the vibration signal of the measurement sensor connected to the channel two to calculate the response spectrum ratio; the secondary response spectrum Ratio calculation part 185, which divides the measured value of the frequency domain transformation and the response spectrum of the vibration signal including the vibration signal connected to the channel 1 based on the second excitation and the response spectrum by the vibration signal including the reference sensor connected to the channel 2 The frequency domain transformation of the signal and the measured value of the response spectrum are used to calculate the response spectrum ratio; the seismic recorder response acquisition unit 190, by comparing the response spectrum ratio calculated by the above-mentioned primary response spectrum ratio calculation unit and the secondary response spectrum ratio calculation unit Multiply and multiply by 1/2 squared to get (channel one seismograph response/channel two seismograph response).

In this case, it is characterized in that the value acquired by the above-mentioned seismograph response acquisition unit is multiplied by (the seismograph response of channel 1/the seismograph of channel 2 response) to correct the seismic recorder.

Hereinafter, a detailed description will be given by using an embodiment of a measurement sensor verification system and method based on the seismograph of the present invention.

FIG. 1 is a configuration diagram schematically constituting a measurement sensor verification system using a seismograph according to an embodiment of the present invention.

As shown in Figure 1, as the measuring sensor verification system utilizing seismograph of the present invention, it is characterized in that, comprises: Vibration device 200, it excites vibrating table; And measurement sensor 45, is connected with above-mentioned excitation device, and excites when starting vibration device; Signal line converter 600, it is used to measure the vibration signal of above-mentioned reference sensor and measuring sensor; Signal amplifier 700, it will be in The vibration signal converted in the above-mentioned signal line converter is amplified, and the amplified vibration signal is transmitted to the seismic recorder; the seismic recorder 800 is used to measure the amplified vibration signal transmitted from the above-mentioned signal amplifier; the seismic data processing mechanism 100, its Calculate the response spectrum of the reference sensor and measuring sensor measured by the above-mentioned seismograph, convert the vibration signal of the reference sensor into the frequency domain, correct it into a response spectrum, and compare the above-calculated response spectrum with the corrected response spectrum To judge whether the measuring sensor is normal or not.

The reference sensor and the measurement sensor are fully and well fastened to the upper surface of the vibrating table, so that the vibration generated on the vibrating table can be transmitted smoothly without distortion.

Referring to FIG. 1 , in order to test many kinds of reference sensors 400 and measurement sensors 450 , a signal line converter 600 , a sensor power supply device 600 for supplying power to the sensors, and amplifying the signals of sensors with low output voltages into The signal amplifier 700 suitable for the seismic recorder 800 is constituted.

In the above-mentioned signal line converter, each manufacturer uses different terminals, so wiring needs to be changed so that the signal line cable can be input to the seismic recorder. Another term can be defined as a signal line wiring converter.

At this time, for example, the response spectrum of the sensor is known in the reference sensor or the measurement sensor, or the experiment is performed on the assumption that the sensor behaves almost the same as the theoretical response spectrum.

FIG. 2 is a block diagram of a seismic data processing mechanism of a measurement sensor verification system using a seismograph according to an embodiment of the present invention.

As shown in Fig. 2 , the seismic data processing mechanism 100 includes: a reference sensor frequency domain conversion unit 110, which obtains the vibration signal of the reference sensor, and transforms it into a frequency domain by Fourier transform; a response spectrum correction unit 120, which is used to convert The frequency spectrum converted by the reference sensor frequency domain conversion unit is divided by the theoretical response spectrum or the known response spectrum to correct the response spectrum; the same sensor determination unit 130 is used to determine whether the reference sensor and the measurement sensor are the same type of sensor; The unit conversion unit 140 is used to convert the reference sensor and the measurement sensor using a differential operator or an integral operator in the frequency domain based on the measurement sensor so as to have the same measurement unit; the response spectrum Calculation section 150, which is used to calculate the response spectrum by band-pass filter filtering and Hilbert transform when the reference sensor and the measurement sensor are of the same type; The response spectrum calculated by the above-mentioned response spectrum calculation part is compared with the response spectrum corrected by the response spectrum correction part to determine whether the measurement sensor is normal or not; the result value output part 170 is used to output the normal or non-normality of the measurement sensor in the above-mentioned measurement sensor. Judgment result value.

The response spectrum correction unit 120 corrects the response spectrum by dividing the spectrum converted by the reference sensor frequency domain conversion unit by a theoretical response spectrum provided by the manufacturer or a known response spectrum.

The theoretical response spectrum described in the present invention refers to the inherent characteristics of the sensor manufactured by the seismograph sensor manufacturing company, and the response spectrum is provided by providing pole (pole) and zero (zero) values.

In addition, the response spectrum of each sensor directly measured by the manufacturer is provided to the user. The known response spectrum refers to the response spectrum supplied by the manufacturer or the response spectrum measured by absolute measurement, and can be used as an absolute reference sensor.

At this time, the above-mentioned unit conversion unit 140 performs a function that can be converted with the same measurement unit. For example, there may be an accelerometer for measuring the acceleration of the ground motion and a speedometer for measuring the speed. Therefore, the measurement of the above-mentioned accelerometer and the speedometer may be unified. unit.

In addition, the seismic data described in the specification of the present invention means data showing ground vibration values measured over time.

FIG. 4 is a flowchart showing a measurement sensor verification method using a seismograph according to an embodiment of the present invention.

As shown in Figure 4, the measuring sensor verification method using a seismograph is characterized in that it includes the following steps: a vibration scanning step (S100) of the vibration excitation device vibrating the vibration table by starting the vibration excitation device; Obtaining the experimental information acquisition step (S200) of obtaining the vibration signals of the reference sensor and the measuring sensor measured by the seismograph; the reference sensor frequency domain conversion unit 110 obtains the vibration signal of the reference sensor, and transforms it into a reference sensor in the frequency domain by using Fourier transform Frequency domain transformation step (S300); the response spectrum correction unit 120 divides the frequency domain transformed by the reference sensor frequency domain transformation unit by the theoretical response spectrum correction of the reference sensor or the known response spectrum to correct the response spectrum of the response spectrum Step (S400); same-type sensor determination unit 130 determines whether the reference sensor and the measurement sensor are the same type of sensor (S500); when the reference sensor and the measurement sensor are the same type of sensor, the response spectrum The calculation part 150 uses the band-pass filter to filter and the response spectrum calculation step (S700) of Hilbert transform to calculate the response spectrum; The response spectrum corrected by the correction part is compared to judge whether the measuring sensor is normal or not. The step of determining whether the measuring sensor is normal or not (S800); the result value output part 170 outputs the result value output of the result value judged by the measuring sensor normal or not judging part. step.

At this time, if the reference sensor and the measurement sensor are not the same type of sensor as a result of the determination in the above-mentioned same-type sensor determination step (S500), the unit conversion unit 140 performs a differential operator or an integral operation in the frequency domain using the measurement sensor as a reference. A unit conversion step (S600) in which symbols are converted into the same measurement unit.

Referring to Fig. 2 and Fig. 4, in order to verify the sensor, the seismic data processing mechanism 100 includes: a reference sensor frequency domain conversion unit 110, a response spectrum correction unit 120, a same sensor determination unit 130, a unit conversion unit 140, and a response spectrum calculation unit 150 . The measuring sensor is normal or not judging part 160 and the result value output part 170 .

To describe the working process, the vibration table is excited (S100) by activating the vibration excitation device, and the vibration signals of the reference sensor and the measurement sensor measured by the seismograph are obtained in the seismic data processing mechanism 100 (S200).

Thereafter, the reference sensor frequency domain conversion unit 110 of the seismic data processing unit acquires the vibration signal of the reference sensor, and transforms it into the frequency domain by Fourier transform ( S300 ).

At this time, when the vibrating table performs excitation, the excitation source includes frequency components in the vibration range to be obtained.

For example, when the vibration performs a vibration sweep (sweep) of 0.1 Hz to 50 Hz, excitation is performed for about 20 seconds.

After that, it is confirmed whether the test data is transmitted from the seismic recorder and obtained by the seismic data processing facility.

Thereafter, the reference sensor data is converted into the frequency domain by Fourier transform in the reference sensor frequency domain conversion unit.

Thereafter, the response spectrum correction unit 120 divides the frequency domain converted by the reference sensor frequency domain conversion unit by the theoretical response spectrum correction of the reference sensor or a known response spectrum to correct the response spectrum ( S400 ).

Thereafter, the same-type sensor determination unit 130 determines whether the reference sensor and the measurement sensor are of the same type (S500). If the reference sensor and the measurement sensor are not the same type of sensor, the unit conversion unit 140 performs A unit conversion step ( S600 ) of converting the frequency domain into the same measurement unit using a differential operator or an integral operator.

That is, the measurement sensor is used as a reference, and the differential operator (iw) or the integral operator (1/iw) is used in the frequency domain to convert to the same measurement unit.

In this case, w represents each frequency, and i represents a complex number.

In particular, differentiation is performed in the frequency domain of the reference sensor data when the reference sensor is velocity and the measurement sensor is acceleration, and integration is performed in the frequency domain of the reference sensor data when the reference sensor is acceleration and the measurement sensor is velocity .

Then, when the reference sensor and the measurement sensor are the same type of sensor, the response spectrum calculation unit 150 calculates the response spectrum by bandpass filtering and Hilbert transform ( S700 ).

This will be described in more detail with reference to FIGS. 3 and 5 .

Thereafter, the measurement sensor normality determination unit 160 compares the response spectrum calculated by the response spectrum calculation unit with the response spectrum corrected by the response spectrum correction unit to determine whether the measurement sensor is normal or not ( S800 ).

Thereafter, the result value output unit 170 outputs the result value determined by the measurement sensor normality determination unit, and if it outputs a conclusion that the measurement sensor is abnormal, correction or repair of the sensor is determined.

At this time, the result itself of the calculated response spectrum can also be obtained.

3 is a block diagram of a response spectrum calculation unit of a measurement sensor verification system using a seismograph according to an embodiment of the present invention.

As shown in FIG. 3 , the above-mentioned response spectrum calculation part 150 includes: a real number part and an imaginary number part setting part 151, which performs time-domain bandpass filter filtering using measurement sensor information and reference sensor information in the set bandpass frequency band, and After the Hilbert transform is performed, the bandpass filter filtering information is set as the real number part, and the Hilbert transform information is set as the imaginary number part; the amplitude and phase calculation part 152 is used to total the real number part and the imaginary number part After the real number part and imaginary number part set by the setting part, calculate the amplitude and phase that change with time; the amplitude extraction part 153, for the time series amplitude of the reference sensor and the measurement sensor calculated by the above-mentioned amplitude and phase calculation part, Extract only amplitude values that are higher than the noise of the reference sensor or the measurement sensor itself; the response amplitude spectrum ratio calculation unit 154 calculates the amplitude ratio by summing the time-series amplitude values of the reference sensor and the measurement sensor for the amplitude extracted by the above-mentioned amplitude extraction unit. (measurement sensor/reference sensor) to calculate the response amplitude spectrum ratio; the response phase spectrum calculation part 155 is used to calculate the phase difference between the reference sensor and the measurement sensor (measurement sensor phase-reference sensor phase), and calculate the average value of the phase difference and The standard deviation value, in the phase difference value, calculates the average value of the phase difference values within the standard deviation range in the average value to calculate the response phase spectrum.

FIG. 5 is a flowchart showing a processing procedure of a response spectrum calculation unit of a measurement sensor verification system using a seismograph according to an embodiment of the present invention.

Referring to Fig. 3 and Fig. 5, the above-mentioned response spectrum calculation step (S700) is described. As far as the response spectrum calculation step is concerned, it is characterized in that it includes the following steps: the real number part and the imaginary number part setting part 151 use in the set bandpass frequency band After the measurement sensor information and reference sensor information are filtered with a time-domain bandpass filter and Hilbert transform is performed, the bandpass filter information is set as the real part and the Hilbert transform information is set as the imaginary part. Real number part and imaginary number part setting step (S710); the amplitude and phase calculation part 152 calculates the amplitude and phase that change with time after summing the real number part and imaginary number part set by the real number part and imaginary number part setting part Amplitude and phase calculation step (S720); the amplitude extraction unit 153 extracts only the amplitude value higher than the amplitude value of the noise of the reference sensor or the measurement sensor itself for the time-series amplitudes of the reference sensor and the measurement sensor calculated by the amplitude and phase calculation unit Extraction step (S730): The response amplitude spectrum ratio calculation unit 154 calculates the response amplitude spectrum by summing the time-series amplitude values of the reference sensor and the measurement sensor to calculate the amplitude ratio (measurement sensor/reference sensor) for the amplitude extracted by the amplitude extraction unit Response amplitude spectrum ratio calculation step (S740); response phase spectrum calculation section 155 calculates the phase difference between the reference sensor and the measurement sensor (measurement sensor phase-reference sensor phase), and calculates the average value and standard deviation of the phase difference. The phase difference value is calculated from the mean value in a response phase spectrum calculation step ( S750 ) for calculating the response phase spectrum by the mean value of the phase difference values within the range of the standard deviation.

That is, the real number part and imaginary number part setting part 151 performs time-domain bandpass filter filtering using measurement sensor information and reference sensor information in the set bandpass frequency band, and after performing Hilbert transform, the bandpass filter The information is set to the real number part, and the Hilbert transform information is set to the imaginary number part ( S710 ).

As a procedure for creating a response spectrum by filtering with a bandpass filter (minimum phase delay filter or zero phase filter) and using Hilbert transform, the minimum frequency is calculated from 1/T (T: total recording time) to the maximum frequency 1/(2X×t), it is theoretically most appropriate to set the increase and decrease frequency of each frequency band to 1/T, but considering the calculation time, for example, when the band pass band is less than 1Hz, set the band to 1.5/T, When the frequency is above 1 Hz, the band is set to 0.25 Hz, and the steps of band-pass filter and Hilbert transformation are repeatedly performed from the lowest frequency to the highest frequency.

Thereafter, the amplitude and phase calculation unit 152 calculates the amplitude and phase that change with time after adding up the real number part and the imaginary number part set by the real number part and imaginary number part setting unit ( S720 ).

Thereafter, the amplitude extracting unit 153 extracts only amplitude values that are higher than the noise of the reference sensor or the measurement sensor itself from the time-series amplitudes of the reference sensor and the measurement sensor calculated by the amplitude and phase calculation unit (S730), and calculates the response amplitude spectrum ratio. The unit 154 calculates an amplitude ratio (measurement sensor/reference sensor) by summing up the time-series amplitude values of the reference sensor and the measurement sensor with respect to the amplitude extracted by the amplitude extraction unit ( S740 ).

Thereafter, the response phase spectrum calculation unit 155 calculates the phase difference between the reference sensor and the measurement sensor (measurement sensor phase−reference sensor phase), calculates the average value and standard deviation value of the phase difference, and calculates the standard deviation of the average value from the phase difference value calculating the response phase spectrum by averaging the phase difference values within the range (S750).

Through the steps described above, frequency-dependent response spectra are obtained for different frequencies.

The method of obtaining the sensor response spectrum as described above is also applicable to the noise and seismic observation seismic data of general seismic observation stations.

FIGS. 6 to 7 are diagrams of a test execution procedure diagram for correcting a seismograph correction of a measurement sensor verification system using a seismograph according to an embodiment of the present invention.

First of all, channel one is connected to the reference sensor, channel two is connected to the measurement sensor, the first excitation source is executed as shown in Figure 6, channel one is connected to the measurement sensor, channel two is connected to the reference sensor, and the second excitation source is executed as Figure 7.

8 is a block diagram of a seismograph for modifying a measurement sensor verification system using a seismograph according to an embodiment of the present invention.

9 is a flowchart showing a seismograph correction method using the measurement sensor verification system for a seismograph according to an embodiment of the present invention.

After connecting Channel 1 to the reference sensor, and Channel 2 to the measuring sensor, perform the first excitation (S1000). The frequency-domain transformation of the vibration signal and the measurement value of the response spectrum are divided by the measurement value of the frequency-domain transformation and the response spectrum of the vibration signal including the measurement sensor connected to the channel 2 to calculate a response spectrum ratio ( S1100 ).

Thereafter, vibration is performed after the sensor channels are mutually changed (S1200), and the secondary response spectrum ratio calculation unit 185 converts the frequency domain transform and the response spectrum of the vibration signal including the measurement sensor connected to channel one based on the second vibration The response spectrum ratio is calculated by dividing the measurement value of the frequency domain transformation and the response spectrum measurement value including the vibration signal of the reference sensor connected to the second channel ( S1300 ).

Thereafter, the seismograph response acquisition unit 190 multiplies the response spectrum ratio calculated by the primary response spectrum ratio calculation unit and the secondary response spectrum ratio calculation unit and performs 1/2 square acquisition (seismograph response of channel one/channel two seismograph response) (S1400).

This will be described in detail as follows.

usable $\frac{Z_{11}\left(w\right)}{Z_{\mathrm{twenty\; one}}\left(w\right)}=\frac{R_1{S}_A}{R_2{S}_B},$ $\frac{Z_{12}\left(w\right)}{Z_{\mathrm{twenty\; two}}\left(w\right)}=\frac{R_1{S}_B}{R_2{S}_A}$ export $\frac{Z_{11}\left(w\right)}{Z_{\mathrm{twenty\; one}}\left(w\right)}\×\frac{Z_{12}\left(w\right)}{Z_{\mathrm{twenty\; two}}\left(w\right)}={\left(\frac{R_1}{R_2}\right)}^2.$

At this time, Z ₁₁ (w)=R ₁ S _A U ₁ (w): the frequency domain value measured on channel 1 based on the first excitation

Z ₂₁ (w)=R ₂ S _B U ₁ (w): Frequency domain value measured in channel 2 based on the first excitation

Z ₁₂ (w)=R ₁ S _A U ₁ (w): Frequency domain value measured at channel 1 based on the second excitation

Z ₂₂ (w)=R ₂ S _B U ₁ (w): Frequency domain value measured in channel 2 based on the second excitation

But, R ₁ : Channel 1 recorder response, R ₂ : Channel 2 recorder response

S _A : A sensor response, S _B : B sensor response

U ₁ : first excitation, U ₂ : second excitation

At this time, the value acquired by the above-mentioned seismograph response acquisition unit is corrected by multiplying the vibration value measured from the sensor connected to channel 2 by (the seismograph response of channel 1/the seismograph response of channel 2) Seismograph.

That is, if the value derived above is multiplied by the recorded value of channel two (R ₁ /R ₂ ), a correction to the seismograph can be performed.

Specifically, when the seismograph value derived from channel two is multiplied by (the seismograph response of channel one / the seismograph response of channel two), then the sensor values obtained on channel one and channel two are calculated as the same channel A record response.

The recorder values recorded in each channel can be regarded as "(channel recorder response spectrum X sensor response spectrum) X (foundation vibration spectrum)" in the frequency domain.

Therefore, if the response spectrum ratio related to the channel recorder is obtained, and the recorder response of the actual channel 1 is obtained, the response spectrum of channel 2 can also be known.

Generally, most recorders provide voltage/bit (voltage/bit) or voltage/count (voltage/count), and if the response spectrum ratio is obtained, it is expressed as a specified constant value according to different frequencies.

Those of ordinary skill in the technical field to which the present invention pertains as described above can understand that the present invention can be implemented in other specific forms without changing the technical idea or essential features of the present invention. Therefore, it should be understood that the embodiments described above are illustrative in all respects and not limited thereto.

The scope of the present invention depends more on the appended claims than the above detailed description, and it should be construed that all changes or modifications derived from the meaning and scope of the claims and their equivalent concepts are included in the scope of the present invention. .

Claims (8)

1. a determination sensor verification system that utilizes seismograph is characterized in that, comprising:

Exciting device (200), it carries out exciting to shaking table;

Shaking table (300), it is provided with reference sensor (400) and determination sensor (45), be connected with above-mentioned exciting device, and when starting exciting device exciting;

Signal wire transducer (600), it is used for measuring the vibration signal of said reference sensor and determination sensor,

Signal amplifier (700), its vibration signal to conversion in above-mentioned signal wire transducer amplifies, and transmits the vibration signal that amplifies to seismograph;

Seismograph (800), it is used for measuring the amplification vibration signal that transmits from above-mentioned signal amplifier;

Seismic data processing mechanism (100), the reference sensor that its calculating is measured in above-mentioned seismograph and the response spectra of determination sensor, the vibration signal of reference sensor is transformed to frequency domain, and it is modified to after the response spectra, to the response spectra of the response spectra of above-mentioned calculating and correction compare to judge determination sensor normally whether.

2. the determination sensor verification system that utilizes seismograph according to claim 1 is characterized in that, above-mentioned seismic data processing mechanism (100) comprising:

Reference sensor frequency domain transformation component (110), the vibration signal that it obtains reference sensor utilizes Fourier transform that it is transformed to frequency domain;

Response spectra correction portion (120), it is used for by coming the frequency spectrum of conversion to revise response spectra divided by theoretical property response spectra or known response spectra by said reference sensor frequency domain transformation component;

Sensor judging part of the same race (130), it is used for the judgment standard sensor and whether determination sensor is the sensor of identical type;

Unit transformation section (140), it is used for when reference sensor and determination sensor are not the sensor of identical type, take determination sensor as benchmark, uses differentiate symbol or integral operation to accord with at frequency domain and carries out conversion, in order to possess identical analytical unit;

Response spectra calculating part (150), it is used for when reference sensor and determination sensor are the sensor of identical type, utilizes band-pass filter and Hilbert transform to come calculated response to compose;

Determination sensor is judging part (160) whether normally, its be used for to the response spectra that calculates by above-mentioned response spectra calculating part and response spectra by the correction of response spectra correction portion compare to judge determination sensor normally whether;

End value efferent (170), it is used for output in the end value of the whether normal judgement section judges of said determination sensor.

3. the determination sensor verification system that utilizes seismograph according to claim 2 is characterized in that, above-mentioned response spectra calculating part (150) comprising:

Real part and imaginary part configuration part (151), it utilizes determination sensor information and reference sensor information and executing time domain band-pass filter in the band passband of setting, and after the execution Hilbert transform, being real part with the band-pass filter information setting, is imaginary part with the Hilbert transform information setting;

Amplitude and phase calculation section (152), after it be used for to add up to the real part and imaginary part of setting by real part and imaginary part configuration part, amplitude and phase place that calculating changes in time;

Amplitude extraction section (153), it only extracts the amplitude of the noise that is higher than reference sensor or determination sensor itself for by above-mentioned amplitude and the reference sensor of phase calculation section calculating and the time series amplitude of determination sensor;

Response amplitude spectrum is than calculating part (154), and it is for the amplitude that extracts by above-mentioned amplitude extraction section, and the time series amplitude calculated amplitude ratio by aggregated basis sensor and determination sensor is that determination sensor/reference sensor comes the calculated response spectral amplitude;

Response phase spectrum calculating part (155), its phase differential that is used for calculating reference sensor and determination sensor is determination sensor phase place-reference sensor phase place, calculate mean value and the standard deviation value of phase differential, in phase difference value, come the calculated response phase spectrum at fall into a trap the get it right mean value of the phase difference value in the standard deviation scope of mean value.

4. the determination sensor verification system that utilizes seismograph according to claim 1 is characterized in that, above-mentioned seismic data processing mechanism also comprises:

A response spectra is than calculating part (180), its by will comprise with based on the first time exciting the frequency domain conversion of vibration signal of the reference sensor that is connected of channel one and the measured value of response spectra come calculated response spectrum ratio divided by the frequency domain conversion of the vibration signal that comprises the determination sensor that is connected with channel two and the measured value of response spectra;

Quadratic response spectrum is than calculating part (185), its by will comprise with based on the second time exciting the frequency domain conversion of vibration signal of the determination sensor that is connected of channel one and the measured value of response spectra come calculated response spectrum ratio divided by the frequency domain conversion of the vibration signal that comprises the reference sensor that is connected with channel two and the measured value of response spectra;

Seismograph response obtaining section (190), it will be by responding than multiplying each other and multiply by 1/2 square of seismograph of obtaining the seismograph response/channel two of channel one than the response spectra that calculating part calculates than calculating part and quadratic response spectrum by an above-mentioned response spectra.

5. the determination sensor verification system that utilizes seismograph according to claim 4, it is characterized in that, will above-mentioned seismograph response obtaining section obtain on duty be that the seismograph of the seismograph response/channel two of channel one responds to revise seismograph with the foundation vibration value of measuring from the sensor that is connected with channel two.

6. a determination sensor verification method that utilizes seismograph is characterized in that, may further comprise the steps:

Make the exciting device oscillating scanning step (S100) of shaking table vibration by starting exciting device;

Obtain the experiment information of the vibration signal of the reference sensor measured by seismograph and determination sensor at the seismic data processing mechanism and obtain step (S200);

Reference sensor frequency domain transformation component (110) is obtained the vibration signal of reference sensor, utilizes Fourier transform to make it be transformed to the reference sensor frequency domain shift step (S300) of frequency domain;

Response spectra correction portion (120) is by reference sensor frequency domain transformation component and the response spectra correction step (S400) that the frequency domain of conversion is revised response spectra divided by the correction of theoretical property response spectra or the known response spectra of reference sensor;

Whether sensor judging part of the same race (130) judgment standard sensor and determination sensor are the sensor determining step of the same race (S500) of the sensor of identical type;

Be that response spectra calculating part (150) utilizes the response spectra calculation procedure (S700) of band-pass filter and Hilbert transform calculated response spectrum in the situation of sensor of identical type at reference sensor and determination sensor;

Determination sensor is the normal determining step (S800) whether of the whether normal determination sensor of judging part (160) by the response spectra that calculates by the response spectra calculating part and response spectra by the correction of response spectra correction portion being compared judge determination sensor whether normally;

End value efferent (170) output is in the end value output step of the end value of the whether normal judgement section judges of determination sensor.

7. the determination sensor verification method that utilizes seismograph according to claim 6, it is characterized in that, result in above-mentioned sensor determining step of the same race (S500) judgement, when reference sensor and determination sensor were not the sensor of identical type, executable unit's transformation component (140) accorded with being transformed to the unit transformation step (S600) of identical analytical unit with differentiate symbol or integral operation in frequency domain take determination sensor as benchmark.

8. the determination sensor verification method that utilizes seismograph according to claim 6 is characterized in that, above-mentioned response spectra calculation procedure (S700) may further comprise the steps:

Real part and imaginary part configuration part (151) utilize determination sensor information and reference sensor information and executing time domain band-pass filter in the band passband of setting, and after the execution Hilbert transform, being real part with the band-pass filter information setting, is that real part and the imaginary part of imaginary part set step (S710) with the Hilbert transform information setting;

After amplitude and phase calculation section (152) add up to the real part set by real part and imaginary part configuration part and imaginary part, calculate the amplitude that changes in time and amplitude and the phase calculation step (S720) of phase place;

Amplitude extraction section (153) only extracts the amplitude extraction step (S730) of the amplitude of the noise that is higher than reference sensor or determination sensor itself for by amplitude and the reference sensor of phase calculation section calculating and the time series amplitude of determination sensor;

For the amplitude that extracts by amplitude extraction section, the time series amplitude calculated amplitude by aggregated basis sensor and determination sensor is than being that determination sensor/reference sensor comes the response amplitude of calculated response Amplitude spectrum ratio to compose than calculation procedure (S740) than calculating part (154) for response amplitude spectrum;

The phase differential that response phase spectrum calculating part (155) calculates reference sensor and determination sensor is determination sensor phase place-reference sensor phase place, and calculate mean value and the standard deviation value of phase differential, in phase difference value, fall into a trap to get it right at mean value and come the response phase spectrum calculation procedure (S750) of calculated response phase spectrum in the mean value of the phase difference value in the standard deviation scope.

Images