JP5064946B2 - Predicted seismic intensity calculation device for alarm, earthquake alarm notification system - Google Patents

Description

  The present invention relates to a technique for calculating a predicted value of a measured seismic intensity with high accuracy from an acceleration of ground motion during an earthquake.

  When an earthquake warning is issued in real time during an earthquake, a mechanism may be considered in which the magnitude of the earthquake motion is measured with an earthquake motion index for determining whether there is an alarm, and an alarm is issued if a reference value (e.g., equivalent to a seismic intensity of less than 5) is exceeded. At this time, an easy-to-understand earthquake motion index, that is, a highly universal one is required. Of the various indices representing the strength of seismic motion, the most universal is considered to be the Japan Meteorological Agency seismic intensity class.

This JMA seismic intensity class is determined by the measured seismic intensity. The measured seismic intensity is a value calculated by performing a calculation process such as a filter process according to the method disclosed in Non-Patent Document 1 based on the acceleration waveform from the beginning to the end of the ground motion (for example, (See Patent Document 1.)
In addition, as a mechanism for measuring the magnitude of ground motion during an earthquake and issuing a seismic warning when the reference value is exceeded, 5 Hz acceleration (vertical motion, horizontal composite motion, 3 components) exemplified in Non-Patent Document 2, The magnitude of seismic motion is measured by the real-time SI value exemplified in Patent Document 3 and the product-sum average value (DI value) of speed and acceleration exemplified in Non-Patent Document 4, and an alarm is issued when a certain value is exceeded. The mechanism is being used. In addition, about the method using 5 Hz acceleration, there exists an example applied to the Shinkansen and a conventional line along a seismometer. Moreover, there exists an example used in gas supply infrastructure about the method using real-time SI value. Further, the technique using the DI value is also disclosed in Patent Documents 1 to 3.

In addition, as a mechanism for predicting the final seismic intensity (DI value) from the very first part of the ground motion, the PI value is applied to the initial part of the ground motion using the DI value calculation method exemplified in Non-Patent Document 5. There are examples. In addition, there is an example used for the seismometer of the “Shinkansen early earthquake detection system” as a mechanism for predicting the PI value.
JP 2000-292547 A JP 2001-147273 A JP 2006-10664 A 1996 Japan Meteorological Agency announcement No.4 "Earthquake warning system in JNR" Yutaka Nakamura, 1985, Railway Technology 42-10, pages 371-376 “Development of Intelligent Seismic Sensors” Furukawa et al., 1999, “Saving Review” vol. 17 "Examination of rational seismic intensity index values" Yutaka Nakamura, 2003, JSCE Earthquake Engineering Papers "Approximation method of measured seismic intensity equivalent value suitable for real-time sequential processing" Kotokaku, Makoto Aoi, Hiromitsu Nakamura, Hiroyuki Fujiwara, Shigeki Adachi (National Institute for Disaster Prevention) Autumn Meeting of the Seismological Society of Japan, 2006, p. 128

  However, in the technique disclosed in Non-Patent Document 1, the measured seismic intensity can be calculated only after the earthquake ends, and cannot be calculated during the earthquake. In addition, since the calculation process of the measured seismic intensity is complicated, it is required that the apparatus for calculating the measured seismic intensity be equipped with a high-performance, that is, expensive arithmetic unit and complicated software.

  Further, the technique disclosed in Non-Patent Document 2 has a problem that an error occurs between the calculated measured seismic intensity and the measured seismic intensity obtained by the Japan Meteorological Agency because the frequency characteristics are different even with the same acceleration value. Specifically, an error of about ± 0.3 occurs (see Non-Patent Document 4). For this reason, an unnecessary alarm may be generated.

  Further, since the technology disclosed in Non-Patent Document 3 has characteristics closer to the seismic intensity measured by the Japan Meteorological Agency than the technology disclosed in Non-Patent Document 2, the calculated seismic intensity and the seismic intensity measured by the Japan Meteorological Agency Although the error is smaller than when the technique disclosed in Non-Patent Document 2 is used, there is still a problem that an error occurs. Specifically, an error of about ± 0.3 occurs even with the same SI value (see Non-Patent Document 4). For this reason, an unnecessary alarm may be generated.

  The technology disclosed in Non-Patent Document 4 has characteristics close to the seismic intensity measured by the Japan Meteorological Agency, but there is a problem that an error occurs between the calculated seismic intensity calculated by the Japan Meteorological Agency. Specifically, an error of about ± 0.2 occurs even with the same DI value (see Non-Patent Document 4). For this reason, an unnecessary alarm may be generated.

  The technique disclosed in Non-Patent Document 5 aims to predict the later seismic intensity by applying the DI value calculation method to the initial motion part of the ground motion. Research results have not been published. In addition, about this technique, it is thought that a correlation with a seismic intensity is inferior to the case where DI value is used in principle, and there is a possibility that many unnecessary alarms may be generated for that reason.

  The present invention has been made in view of such problems, and the object of the present invention is to calculate the predicted value of the measured seismic intensity in real time from the acceleration of the earthquake motion during the earthquake without waiting for the end of the earthquake. Is to provide.

According to the alarm predictive seismic intensity calculation device according to claim 1 made to solve the above-described problems, the acquisition means (51: In this section, in order to facilitate understanding of the invention, the “invention The sign used in the “Best Mode for Implementation” column is attached, but this sign does not mean that the scope of claims is limited.) However, the acceleration of the vertical motion of the earthquake motion measured by an acceleration sensor or the like. The value is acquired, and the correction unit corrects the acceleration value of the vertical motion of the earthquake motion acquired by the acquisition unit. Then, the calculation means calculates the alarm predicted seismic intensity Iap, which is a predicted value of the measured seismic intensity, with respect to the vertical acceleration value corrected by the correcting means , the maximum value between the arrival of the seismic wave and the sampling time (i), Afu [ n] and predicted alarm seismic intensity Iap using the following relational expression (1).

Relational expression (1): Iap = 2.0αlog (Afu [n]) + 2.0β + 0.94
However, the coefficient α and the coefficient β are regression coefficients calculated in advance by regression analysis.
According to the seismic intensity prediction apparatus of the present invention configured as described above, the predicted value of the measured seismic intensity can be calculated in real time from the acceleration of the seismic motion measured during the occurrence of the earthquake without waiting for the end of the earthquake.

  Further, as illustrated in FIG. 5, the earthquake alarm based on the predicted seismic intensity for warning includes a slight error from the final seismic intensity. On the other hand, as illustrated in FIG. Although it is the frequency of occurrence of alarms, the alarm can be issued earliest, which can contribute to ensuring safety during an earthquake.

  In this case, as in claim 2, it is conceivable that the correcting means corrects the acceleration value of the vertical motion of the seismic motion acquired by the acquiring means using a bandpass filter. As an example, the correcting means corrects the acceleration value of the vertical motion of the earthquake motion using the following relational expression (2) (Claim 3).

Equation (2): y [n] = b 1 x [n] + b 2 x [n-1] + b 3 x [n-2] -a 2 y [n-1] -a 3 y [n-2 ]
However, x [n] is the acceleration amplitude value at the sampling time (i), x [n−1] is the acceleration amplitude value at the sampling time (i−1), and x [n−2] is the sampling time ( i-2) is an acceleration amplitude value, and y [n] is an output at a sampling time iΔt, that is, an acceleration value after correction, and y [n]
n−1] is the output at the sampling time (i−1) Δt, that is, the corrected acceleration value, and y [n−2] is the output at the sampling time (i−2) Δt, that is, after the correction. The coefficient b ₁ , the coefficient b ₂ , the coefficient b ₃ , the coefficient a _2, and the coefficient a ₃ are coefficients set in advance through experiments or the like.

  In this way, in the method on the public notice, the entire seismic waveform is Fourier transformed to the frequency domain, subjected to a predetermined filtering process, and then subjected to the inverse Fourier transform to return the seismic waveform to the time domain. On the other hand, in the present invention, by correcting the acceleration value of the vertical motion of the seismic motion using a bandpass filter, the three processes in the above-mentioned method of the announcement can be performed by the process by one bandpass filter, The time required for processing can be shortened and the burden on the processing apparatus can be reduced.

  By the way, the earthquake alarm notification system according to claim 4 made to solve the above-described problem includes a plurality of alarm predictive seismic intensity calculating devices according to any one of claims 1 to 3, a central processing unit, The plurality of alarm predicted seismic intensity calculation devices and the central processing unit can transmit and receive data to and from each other.

  Among them, each of the plurality of alarm predicted seismic intensity calculation devices includes a device-side communication unit that transmits / receives various information to / from the central processing unit, and the alarm predicted seismic intensity Iap and the alarm predicted seismic intensity Iap calculated by the calculating unit. Control device side communication means for controlling alarm-side seismic intensity data including at least information for specifying the time of measurement acceleration value used for calculation and information for specifying the alarm-predicted seismic intensity calculating device of the transmission source to the central processing unit Communication control means for transmitting.

  On the other hand, the central processing unit, center side communication means for transmitting and receiving various information between each of the plurality of alarm predicted seismic intensity calculation devices, storage means for storing the alarm predicted seismic intensity data received by the center side communication means, Analyzes one or more predicted earthquake seismic intensity data stored in the storage means, and based on the analysis result, determines the earthquake alarm to be notified for each alarm predicted seismic intensity calculation device, and the contents of the earthquake alarm Earthquake warning communication means for reporting the earthquake warning determined by the determination means.

  According to the earthquake alarm notification system of the present invention configured as described above, the central processing unit grasps the predicted earthquake intensity for warning transmitted from each alarm predicted earthquake intensity calculator, and calculates the predicted earthquake intensity for each alarm. Determine the earthquake alarm to be communicated for each device and contact it. As a result, if an earthquake occurs, the seismic intensity can be predicted in real time, and earthquake warnings can be issued early in various places as needed.

Embodiments of the present invention will be described below with reference to the drawings.
[First embodiment]
FIG. 1 is a schematic configuration diagram of a seismometer 1 of the first embodiment.

[Description of configuration of seismometer 1]
As shown in FIG. 1, the seismometer 1 as a warning predicted seismic intensity calculation device includes an acceleration sensor 10, an A / D conversion unit 20, a timer control unit 30, a GPS clock 40, a signal processing unit 50, a mode transition unit 60, and An output control unit 70.

The acceleration sensor 10 measures accelerations of three components (vertical, north-south, east-west) of the ground motion, and outputs the measured acceleration values to the A / D converter 20 as analog signals.
The A / D conversion unit 20 converts the analog signal output from the acceleration sensor 10 into a digital signal, adds the time information output from the timer control unit 30 to the signal as a time stamp, and the signal processing unit 50. To the normal processing unit 51.
The time information output from the timer control unit 30 is added to the signal as a time stamp.

  The GPS clock 40 receives an incoming radio wave from a GPS satellite (not shown) by a GPS antenna (not shown), and extracts time information included in the received radio wave. Then, the GPS clock 40 outputs the extracted time information to the timer control unit 30.

  The timer control unit 30 controls the GPS clock 40 and outputs time information output from the GPS clock 40 to the A / D conversion unit 20, the signal processing unit 50, and the mode transition unit 60.

  The signal processing unit 50 includes a normal processing unit 51 and a warning seismic intensity and warning predicted seismic intensity calculation unit (hereinafter, calculation unit) 52. The normal processing unit 51 corresponds to an acquisition unit, and the calculation unit 52 corresponds to a correction unit and a calculation unit.

  Among them, the normal processing unit 51 is a well-known CPU, ROM, RAM, I / O as an input / output circuit, a bus line connecting these components, a signal processing circuit for processing a signal from an input operation, and a calculation unit 52. And a signal output circuit for controlling each unit such as the timer control unit 30. The CPU performs control according to control programs and data stored in the ROM and RAM. The ROM has a program storage area and a data storage area. A control program is stored in the program storage area, and data necessary for the operation of the control program is stored in the data storage area. In addition, the control program operates on the RAM in a form in which the work memory is a work area.

The normal processing unit 51 configured as described above acquires the acceleration value of the seismic motion after digital conversion output from the A / D conversion unit 20 and outputs the acceleration value to the calculation unit 52 inside the signal processing unit 50.
Similarly to the normal processing unit 51, the calculation unit 52 includes a well-known CPU, ROM, RAM, I / O that is an input / output circuit, a bus line that connects these components, and the like. Performs processing to calculate seismic intensity and predicted earthquake intensity for warning.

  The mode transition unit 60 includes a nonvolatile memory and is used to store various data such as a flag indicating whether an earthquake has occurred. The types of flags include an earthquake mode flag indicating that an earthquake has occurred and a normal mode flag indicating that an earthquake has not occurred.

  Similar to the normal processing unit 51, the output control unit 70 includes a known CPU, ROM, RAM, I / O that is an input / output circuit, a bus line that connects these components, and the like. The alarm measured seismic intensity and the alarm predicted seismic intensity calculated by the unit 52 are output to the outside.

[Calculation of measured seismic intensity for warning and predicted seismic intensity for warning]
(A) The present applicant created the following recursive filter (Formula 1) having characteristics close to those defined by the Japan Meteorological Agency by superimposing narrow bandpass filters. That is, the filtered acceleration value y [n] is the acceleration amplitude value x [n] of the time series data of the seismic waveform when the seismic motion is constantly observed, the acceleration amplitude value x [n−1] at the previous sampling timing, Based on x [n-2] and post-filter acceleration values y [n-1] and y [n-2] at the previous sampling timing, the calculation can be performed by the following equation (1). In addition, the characteristic of the filter represented by Formula (1) is shown in FIG. This calculation is performed for each of the acceleration values of the three components (vertical, north-south, east-west).

y [n] = b ₁ x [n] + b ₂ x [n−1] + b ₃ x [n−2] −a ₂ y [n−1] −a ₃ y [n−2] (1)

However, x [n] is the acceleration amplitude value at the sampling time (i), x [n−1] is the acceleration amplitude value at the sampling time (i−1), and x [n−2] is the sampling time ( i−2) is an acceleration amplitude value, y [n] is an output at the sampling time iΔt, that is, an acceleration value after filtering, and y [n−1] is at the sampling time (i−1) Δt. The output, that is, the acceleration value after filtering, and y [n−2] is the output at the sampling time (i−2) Δt, that is, the acceleration value after filtering, and the coefficient b
₁ , coefficient b ₂ , coefficient b ₃ , coefficient a _2, and coefficient a ₃ are coefficients set in advance through experiments or the like.

In the present embodiment, the coefficient b ₁ is set to the numerical value “0.06089422”, the coefficient b ₂ is set to the numerical value “0”, the coefficient b ₃ is set to the numerical value “−0.06089422”, and the coefficient a ₂ is set to a numerical value “−1.88625943”, and the coefficient a ₃ is set to a numerical value “0.898850964”.

(B) Then, the vector composite value v [n] of the post-filter acceleration values y [n] of the three components (up / down, north / south, east / west) obtained until a certain time (sampling time) is v [n] in the calculation method announced by the Japan Meteorological Agency. ], The procedure for obtaining the measured seismic intensity Ia based on v [n] can conform to the published procedure. That is, v [n] is rearranged in descending order, and v [n] that is 30th (in the case of 100 Hz sampling) from the maximum value is amax, and the measured seismic intensity Ia can be obtained by the following equation (2).

Ia = 2.0 log (amax) +0.94 (2)

In addition, FIG. 3 is explanatory drawing (1) which shows the result of having compared the alarm predicted seismic intensity and the measured seismic intensity. In addition, in order to distinguish from the measured seismic intensity Ijma calculated by the method as announced, the measured seismic intensity calculated by the method of this case will be referred to as “measured seismic intensity for warning”. As shown in FIG. 3, with respect to past earthquakes (5403 seismic waveform, magnitude 6.0 or more), the average value of the difference between the alarm seismic intensity I calculated by this method and the measured seismic intensity Ijma was compared. Is 0.083 and the standard deviation is 0.058. Since the effective number of measured seismic intensity is the second decimal place, it can be seen that the alarm seismic intensity reproduces the measured seismic intensity with the necessary and sufficient accuracy.

(C) By the way, the first part of the ground motion is the P wave, and it is well known that the vertical motion is superior to the horizontal motion in the P wave. From this, the applicant thinks that it is better to focus on vertical motion in order to predict the final seismic intensity from the first part of the ground motion, and predicts the final measured seismic intensity from the maximum value of vertical motion. A correlation formula was created. That is, for the three component acceleration values, the alarm seismic intensity filter processing is performed, and the maximum value of the vertical motion acceleration (hereinafter referred to as Afu) is used as the explanatory variable (x), and the 0.3 second rule is applied to the three component composite values. A regression coefficient by power approximation is obtained with the maximum value amax later as y (see FIG. 4), and Expression (3) is created based on this regression coefficient and Expression (2). Equation (3) is an equation for calculating the seismic intensity equivalent value from the maximum value of Afu [n]. By expanding Afu [n] to an arbitrary time, the predicted value Iap (hereinafter “alarm for alarm”) is expanded. This is an equation for calculating the “predicted seismic intensity”.

Iap = 2.0αlog (Afu [n]) + 2.0β + 0.94 (3)

However, the coefficient α and the coefficient β are regression coefficients calculated by regression analysis, and the vertical motion acceleration Afu [n] is obtained when the alarm measurement seismic intensity filter processing is performed on the three component acceleration values at the sampling time i. Vertical acceleration. In the present embodiment, the coefficient α is set to a numerical value “0.9931”, and the coefficient β is set to a numerical value “0.3456”.

  (D) Note that in order to advance the timing for notifying an earthquake alarm, it is only necessary to lower the threshold value for any index, but on the other hand, no alarm is generated and no damage occurs as a result. (Hereinafter referred to as “missing warning”) will be increased. Therefore, when considering the superiority or inferiority of judgment indicators for earthquake warnings in terms of early warning, set a threshold at which the frequency of occurrence of an alarm is the same, and the early warning of an earthquake that requires an alarm at that threshold. Should be considered. Therefore, the threshold of the earthquake warning based on the predicted seismic intensity for warning was set so that it would be the same level as the frequency of the idling warning at 5 Hz acceleration of 40 gal, which has been used for many years as an index of earthquake warnings along railway railway seismometers. FIG. 5 is an explanatory diagram showing the speed of the earthquake warning based on the predicted seismic intensity for warning, and shows that the warning is issued at an earlier timing as the shape projects to the upper left on the graph. As a reference, the same analysis as the predicted seismic intensity for warning was performed for the existing technology 5 Hz acceleration, real-time SI value (see Non-Patent Document 2), and PI value (calculated by the method estimated based on Non-Patent Document 3). . From this, it is understood that the predicted seismic intensity for warning is an index that can notify the warning earlier even when compared with the existing technology.

[Explanation of warning seismic intensity calculation process]
Next, the alarm predicted seismic intensity calculation processing executed by the signal processing unit 50 of the seismometer 1 will be described with reference to the flowchart of FIG. This process is executed when the seismometer 1 is in an energized state.

  First, an acceleration value of earthquake motion is acquired (S105). Specifically, the acceleration sensor 10 constantly measures acceleration values of three components (upper and lower, north and south, east and west), and the A / D converter 20 digitally converts the measured acceleration values. Then, the normal processing unit 51 of the signal processing unit 50 acquires the acceleration value of the seismic motion after digital conversion output from the A / D conversion unit 20 and outputs the acceleration value to the calculation unit 52 inside the signal processing unit 50.

  Subsequently, the calculation unit 52 of the signal processing unit 50 executes acceleration waveform filter processing on the input acceleration value using the above-described equation (1) (S110). With this process, the acceleration value after the filter process is calculated.

Subsequently, the calculation unit 52 calculates the alarm measured seismic intensity Ia and the alarm predicted seismic intensity Iap in each sample using the above-described formulas (2) and (3) (S115).
Subsequently, the normal processing unit 51 refers to the flag stored in the mode transition unit 60 to determine whether or not the normal mode indicating the normal state (S120). If it is determined that the current mode is the normal mode (S120: YES), it is determined whether or not the acceleration value measured in S105 is larger than the reference value (S125). The reference value is set in advance based on the analysis result regarding past earthquakes. If it is determined that the acceleration value is greater than the reference value (S125: YES), it is determined that an earthquake is occurring and the flag stored in the mode transition unit 60 is changed from the normal mode to the earthquake mode (S130). On the other hand, if it is determined that the acceleration value is equal to or less than the reference value (S125: NO), the process directly returns to S105.

  On the other hand, when it is determined that the current mode is not the normal mode but the earthquake mode (S120: NO), the normal processing unit 51 determines whether or not the alarm predicted seismic intensity Iap is larger than the reference value (S135). ). The reference value is set in advance based on the analysis result regarding past earthquakes. When it is determined that the alarm predicted seismic intensity Iap is larger than the reference value (S135: YES), the normal processing unit 51 controls the output control unit 70 to issue an early warning (S140), and the process proceeds to S145. To do. On the other hand, if it is determined that the predicted earthquake intensity Iap for warning is equal to or less than the reference value (S135: NO), the process proceeds to S145 as it is.

  Subsequently, the normal processing unit 51 determines whether or not the alarm measured seismic intensity Ia is larger than the reference value (S145). The reference value is set in advance based on the analysis result regarding past earthquakes. When it is determined that the alarm seismic intensity Ia is greater than the reference value (S145: YES), the normal processing unit 51 controls the output control unit 70 to issue a seismic intensity alarm (S150), and the process proceeds to S155. To do. On the other hand, when it is determined that the alarm seismic intensity Ia is equal to or less than the reference value (S145: NO), the process proceeds to S155 as it is.

  Subsequently, the normal processing unit 51 determines whether or not the acceleration value measured in S105 is smaller than the reference value (S155). The reference value is set in advance based on the analysis result regarding past earthquakes. When it is determined that the acceleration value is smaller than the reference value (S155: YES), the normal processing unit 51 determines that the earthquake has ended and changes the flag stored in the mode transition unit 60 from the earthquake mode to the normal mode. Change (S160) and return to S105. On the other hand, when it is determined that the acceleration value is greater than or equal to the reference value (S155: NO), the process returns to S105 as it is.

[Effect of the first embodiment]
(1) As described above, according to the seismometer 1 of the first embodiment, the acceleration value of the measured ground motion is corrected, the alarm seismic intensity Ia is calculated using the above-described equation (2), and the measurement is performed. The predicted earthquake intensity Iap for warning, which is the predicted value of seismic intensity, is calculated using the above-described equation (3) established between the corrected acceleration value Afu [n] and the predicted earthquake intensity Iap for warning. Therefore, the predicted value of the measured seismic intensity can be calculated in real time from the acceleration of the ground motion measured during the occurrence of the earthquake without waiting for the end of the earthquake.

  As shown in Fig. 4, the alarm seismic intensity is almost the same as the seismic intensity after the end of the earthquake motion. There is an effect of providing an accurate alarm.

  Furthermore, the calculation processing is much simpler than any existing technology, so that the calculation device mounted on the calculation device can be made inexpensive, and the effort required for software development, that is, the development cost can be reduced. . Thereby, the capital investment cost for obtaining an earthquake motion index substantially equivalent to the measured seismic intensity can be reduced.

  Further, as illustrated in FIG. 5, the earthquake alarm based on the predicted seismic intensity for warning includes a slight error from the final seismic intensity. On the other hand, as illustrated in FIG. Although it is the frequency of occurrence of alarms, the alarm can be issued earliest, which can contribute to ensuring safety during an earthquake.

  (2) Further, according to the seismometer 1 of the first embodiment, the acceleration value of the vertical motion of the ground motion is corrected using the bandpass filter of the above formula (1). As a result, in the method on the public notice, the entire seismic waveform is subjected to Fourier transform in the frequency domain, subjected to a predetermined filtering process, and then subjected to computation to return the seismic waveform to the time domain by Fourier inverse transform. By correcting the acceleration value of the vertical motion of the earthquake motion using a bandpass filter, the three processes in the above-mentioned method of announcement can be performed with a single bandpass filter, reducing the time required for the process. And the burden on the processing apparatus can be reduced.

[Other Embodiments]
As mentioned above, although one Embodiment of this invention was described, this invention is not limited to the said embodiment, It is possible to implement in various aspects as follows.

  (1) As illustrated in FIG. 10A, the present invention includes a plurality of seismometers 1 and a central processing unit 3, and the plurality of seismometers 1 and the central processing unit 3 exchange data with each other. You may implement | achieve as the possible earthquake warning alerting | reporting system 100. FIG. Specifically, the seismometer 1 is a device-side communication unit that transmits and receives various types of information to and from the central processing unit 3, and the alarm seismic intensity Iap and the acceleration used when calculating the alarm predicted seismic intensity Iap. Communication control means for controlling the apparatus side communication means to transmit to the central processing unit 3 alarm predicted seismic intensity data including at least information for specifying the value measurement time and the source seismometer 1. On the other hand, the central processing unit 3 includes a center-side communication unit that transmits / receives various information to / from each of the plurality of seismometers 1, a storage unit that stores alarm predicted seismic intensity data received by the center-side communication unit, and a storage unit Analyzes one or more predicted seismic intensity data stored in the alarm, and based on the result of the analysis, the earthquake alarm determination means for determining the earthquake alarm to be communicated for each seismometer 1 and the earthquake warning content determination means And an earthquake alarm communication means for communicating an earthquake alarm.

  In this case, the earthquake warning determination means determines a safety confirmation procedure as an earthquake warning according to the measured seismic intensity. For example, as illustrated in FIG. 10B, when the measured seismic intensity is “0.0 to 3.9”, “no regulation” is determined as the safety confirmation procedure, and the measured seismic intensity is “ If it is “4.0-4.4”, “slow driving” is determined as the safety confirmation procedure, and if the measured seismic intensity is “4.5-4.9”, “partial” If the seismic intensity is “5.0 to 5.4”, “All-ground patrol and partial structure inspection” is determined as the safety confirmation procedure, and the seismic intensity is “5.5. In the case of "-", "all-line ground patrol and structure inspection" is determined as the safety confirmation procedure.

  In this way, the central processing unit 3 grasps the predicted seismic intensity for warning transmitted from the seismometers 1 at each location and determines a safety confirmation procedure as an earthquake warning to be communicated to each seismometer 1 at each location. Contact me. As a result, if an earthquake occurs, the seismic intensity can be predicted in real time, and safety confirmation procedures can be issued early in various places as needed.

(2) In the above embodiment, time information is added to the acceleration value of the earthquake motion at the time of digital conversion. Specifically, the A / D conversion unit 20 converts an analog signal output from the acceleration sensor 10 into a digital signal, and adds the time information output from the timer control unit 30 to the signal as a time stamp. And output to the normal processing unit 51 of the signal processing unit 50. That is, time information is added to the acceleration value of the earthquake motion at the time of digital conversion. On the other hand, time information may be added to the acceleration value of the seismic motion after digital conversion. For example, the normal processing unit 51 acquires the acceleration value of the seismic motion after the digital conversion output from the A / D conversion unit 20, and provides time information with respect to the acquired acceleration value of the seismic motion after the digital conversion. For example, it is added and output to the calculation unit 52. Further, instead of the normal processing unit 51, the calculation unit 52 may add time information to the acceleration value of the seismic motion after digital conversion.
[Second Embodiment]
FIG. 8 is a schematic configuration diagram of the seismometer 2 of the second embodiment.

[Description of configuration of seismometer 2]
As shown in FIG. 8, the seismometer 2 as a warning predicted seismic intensity calculation device includes an acceleration sensor 10, an A / D conversion unit 20, a timer control unit 30, a GPS clock 40, a signal processing unit 80, a mode transition unit 60, and An output control unit 70. In addition, about the acceleration sensor 10, the A / D conversion part 20, the timer control part 30, the GPS timepiece 40, the mode transition part 60, and the output control part 70, it is the same as each structure with which the seismometer 1 of 1st embodiment is provided. Therefore, detailed description thereof is omitted here.

  The signal processing unit 80 includes a normal processing unit 51, a calculation unit 52, and various conventional index calculation units 81. In addition, about the normal process part 51 and the calculation part 52, since it is the same as each structure with which the seismometer 1 of 1st embodiment is provided, the detailed description is abbreviate | omitted here.

  Conventionally, the various index calculation units 81 use the methods disclosed in Non-Patent Documents 2 to 4 above, 5 Hz acceleration (vertical movement, horizontal composite movement, three components), real-time SI value, and product-sum average of speed and acceleration. The magnitude of seismic motion is measured by the value (DI value).

[Explanation of warning seismic intensity calculation process]
Next, the alarm predicted seismic intensity calculation process executed by the signal processing unit 80 of the seismometer 2 will be described with reference to the flowchart of FIG. This process is executed when the seismometer 2 is in an energized state.

  The difference between this process and the alarm predicted seismic intensity calculation process executed by the signal processing unit 50 of the seismometer 1 is that S210 is executed after S115 is executed, and when an affirmative determination is made when S150 is executed In addition, S220 and S230 are executed. Hereinafter, it demonstrates in order.

  First, in S210, various conventional index calculation units 81 perform 5 Hz acceleration (vertical movement, horizontal composite movement, three components), real-time SI value, speed and acceleration by the methods disclosed in Non-Patent Documents 2 to 4 described above. The magnitude of seismic motion is measured by the product-sum average (DI value).

  In S220, the normal processing unit 51 determines whether or not the conventional various indexes are larger than the reference value. The reference value is set in advance based on the analysis result regarding past earthquakes. When it is determined that the various conventional indexes are larger than the reference value (S220: YES), the normal processing unit 51 controls the output control unit 70 to transmit a seismic intensity alarm (S230), and the process proceeds to S155. . On the other hand, when it is determined that the conventional various indexes are below the reference value (S220: NO), the process proceeds to S155 as it is.

[Effects of Second Embodiment]
As described above, according to the seismometer 2 of the second embodiment, the necessity of generating a seismic intensity alarm is determined based on various conventional indexes, so that the accuracy of alarm generation when an earthquake occurs can be increased.

[Other Embodiments]
As mentioned above, although one Embodiment of this invention was described, this invention is not limited to the said embodiment, It is possible to implement in various aspects as follows.

  (1) The present invention is realized as an earthquake alarm notification system 100 that includes a plurality of seismometers 2 and a central processing unit 3 and that can transmit and receive data to and from each other. Also good. Specifically, in the earthquake alarm notification system 100 illustrated in FIG. 10A, the seismometer 2 is installed instead of the seismometer 1.

  In this way, the central processing unit 3 grasps the predicted seismic intensity for warning transmitted from the seismometers 2 at each location, and determines a safety confirmation procedure as an earthquake warning to be communicated to each seismometer 2 at each location. Contact me. As a result, if an earthquake occurs, the seismic intensity can be predicted in real time, and safety confirmation procedures can be issued early in various places as needed.

It is a schematic block diagram of the seismometer of 1st embodiment. It is explanatory drawing which shows the characteristic of the filter used when correcting an acceleration value. It is explanatory drawing (1) which shows the result of having compared the prediction seismic intensity for warning and the measured seismic intensity. It is explanatory drawing which shows correlation with 3D maximum value y [gal] (final time) after a 0.3 second rule application, and maximum value Afu [gal] (maximum time) of the vertical motion acceleration after IIR. It is explanatory drawing (2) which shows the result of having compared the prediction seismic intensity for warning and the measured seismic intensity. It is explanatory drawing which shows the result of having compared the alarm timing according to parameter | index. It is a flowchart explaining the prediction seismic intensity calculation process for warning which the seismometer of 1st embodiment performs. It is a schematic block diagram of the seismometer of 2nd embodiment. It is a flowchart explaining the prediction earthquake intensity calculation process for warning which the seismometer of 2nd embodiment performs. (A) It is a schematic block diagram of an earthquake warning system, (b) is explanatory drawing which shows the correspondence of a measured seismic intensity and a safety confirmation procedure.

Explanation of symbols

DESCRIPTION OF SYMBOLS 1, 2 ... Seismometer, 3 ... Central processing unit, 10 ... Acceleration sensor, 20 ... A / D conversion part, 30 ... Timer control part, 40 ... GPS clock, 50, 80 ... Signal processing part, 51 ... Normal processing part 52 ... Measured seismic intensity for warning and predicted seismic intensity calculation unit for alarm, 60 ... mode transition unit, 70 ... output control unit, 80 ... signal processing unit, 81 ... various conventional index calculation unit, 100 ... earthquake alarm notification system

Claims (4)

An acquisition means for acquiring an acceleration value of the vertical motion of the measured earthquake motion;
Correction means for correcting the acceleration value of the vertical motion of the earthquake motion acquired by the acquisition means;
The alarm predicted seismic intensity Iap, which is the predicted value of the measured seismic intensity, is the maximum value Afu [n] between the arrival of the seismic wave and the sampling time (i) with respect to the vertical motion acceleration value corrected by the correcting means and the alarm. A calculating means for calculating using the following relational expression (1) established between the predicted seismic intensity Iap for use:
A predictive seismic intensity calculation device for warning, comprising:
Relational expression (1): Iap = 2.0αlog (Afu [n]) + 2.0β + 0.94
However, the coefficient α and the coefficient β are regression coefficients calculated in advance by regression analysis. In the alarm prediction seismic intensity calculation device according to claim 1,
The alarm means predictive seismic intensity calculation device, wherein the correction means corrects the acceleration value of the vertical motion of the earthquake motion acquired by the acquisition means using a bandpass filter. In the alarm prediction seismic intensity calculation device according to claim 2,
The said correction | amendment means correct | amends the acceleration value of the vertical motion of the seismic motion acquired by the said acquisition means using the following relational expression (2), The warning predicted seismic intensity calculation apparatus characterized by the above-mentioned.
Equation (2): y [n] = b 1 x [n] + b 2 x [n-1] + b 3 x [n-2] -a 2 y [n-1] -a 3 y [n-2 ]
However, x [n] is the acceleration amplitude value at the sampling time (i), x [n−1] is the acceleration amplitude value at the sampling time (i−1), and x [n−2] is the sampling time ( i-2) is an acceleration amplitude value, y [n] is an output at a sampling time iΔt, that is, an acceleration value after correction, and y [n−1] is an output at a sampling time (i−1) Δt. That is, the corrected acceleration value, and y [n−2] is the output at the sampling time (i−2) Δt, that is, the corrected acceleration value, and the coefficient b ₁ , coefficient b ₂ , coefficient b ₃ , coefficient a ₂ and coefficient a ₃ are coefficients set in advance through experiments or the like. A plurality of alarm predicted seismic intensity calculating devices according to any one of claims 1 to 3 and a central processing unit, wherein the plurality of alarm predictive seismic intensity calculating devices and the central processing unit mutually exchange data. An earthquake alarm notification system capable of transmitting and receiving,
Each of the plurality of predicted earthquake intensity calculating devices for warning,
A device-side communication means for transmitting and receiving various information to and from the central processing unit;
Information for identifying the alarm predicted seismic intensity Iap calculated by the calculating means, the measurement time of the measured acceleration value used when calculating the alarm predicted seismic intensity Iap, and the source alarm predicted seismic intensity calculating device Communication control means for controlling the apparatus-side communication means and transmitting the predicted earthquake intensity data for alarm including at least to the central processing unit,
The central processing unit is
Center-side communication means for transmitting and receiving various information to and from each of the plurality of alarm predicted seismic intensity calculation devices,
Storage means for storing predicted earthquake intensity data for warning received by the center side communication means;
Analyzing one or more predicted earthquake intensity data for warning stored in the storage means, and an earthquake warning determining means for determining an earthquake alarm to be contacted for each predicted earthquake intensity calculation device for each alarm based on the analysis result;
An earthquake warning notification system comprising: an earthquake warning contact means for contacting the earthquake warning determined by the earthquake warning content determination means.