TWI849508B - Method of determining earthquake by intersection of time series and energy of a plurality of seismographs - Google Patents

Abstract

An earthquake detection method for determining earthquakes by timing intersection of multiple seismographs includes the following steps: receiving a first signal triggered in response to a first vibration status at a main detection location and a second signal triggered in response to a second vibration status at a first auxiliary detection location; first determination rule: determining whether the first signal and the second signal both show in triggered status at an intersection time point; second determination rule: determining when the intersection time point occurs, whether triggered time points of the first and the second signals satisfies a timeliness relationship; third determination rule: determining within an energy determination period after the intersection time point occurs, whether an energy relationship between a first vibration intensity measurement value related to the first vibration status and a second vibration intensity measurement value related to the second vibration status does not exceed a specific value range; and when the timeliness relationship satisfies and the energy relationship does not exceed the specific value range, continuing to determine whether the first vibration status is a true seismic event.

Description

Method of judging earthquakes by the intersection of time sequence and energy of multiple seismometers

The present invention relates to an earthquake detection system, an earthquake detection device, and an accurate earthquake detection method, and in particular to an earthquake detection system with multiple detection locations, an earthquake detection device, and an accurate earthquake detection method.

Since the velocity of longitudinal waves is about 6-8 kilometers per second on the surface, and the velocity of transverse waves is about half of that of longitudinal waves on the surface, it is obvious that there is a considerable gap between the arrival times of longitudinal waves and transverse waves at a certain distance from the epicenter. For example, the time gap between the arrival of longitudinal waves and transverse waves at a distance of about 10 kilometers from the epicenter is 3 seconds. If there is a device that can effectively determine whether it is a real earthquake event based on the measurement data of longitudinal waves within this time, there will be an opportunity to take safety measures in time or eliminate non-earthquake events, reducing the equipment and personnel losses that may be caused by earthquakes to a certain limit.

However, earthquake detectors installed in different areas may face different detection environments. For example, earthquake detectors installed near factories are more likely to detect vibrations caused by the operation of machines in the factory, while earthquake detectors installed near homes or office buildings are more likely to be affected by people walking around, or detect vibrations caused by opening and closing doors or knocking on doors. For another example, when a building construction requires laying the foundation, nearby earthquake detectors are more likely to detect larger vibrations and may be judged as earthquakes. For another example, the vibration of the railroad is not a real earthquake. To avoid these interferences, the location of the sensor needs to be a certain distance away from the railroad. Since these vibrations are not caused by real earthquakes, seismic detectors must exclude these vibrations caused by human factors to increase the reliability of earthquake warnings.

In particular, due to different early warning requirements in different application scenarios, for example, a few false alarms for general schools, general residences, and corporate offices will not cause serious economic losses. However, for technology companies (such as wafer factories) that mass-produce products through automated production lines, the shutdown of production lines due to false alarms may cause significant unnecessary economic losses. Therefore, improving the accuracy of earthquake early warnings is particularly important for technology companies.

The Republic of China Patent Publication No. I553327 provides an earthquake detection system, which includes an earthquake data receiving module, a threshold value setting module and an earthquake detection device, which is used to determine whether the new earthquake data is an earthquake event when receiving new earthquake data based on the plurality of earthquake data and the earthquake threshold value, so as to generate a judgment result.

Although the threshold setting module in the above patent can adjust the earthquake threshold value by using the ratio of the short-term average value to the long-term average value to automatically adjust the earthquake threshold value, this method is only applicable to the same detection location with a fixed vibration mode. For example, during the day, the movement of people at work causes the seismometer to frequently detect larger vibrations, while at night, the vibrations detected by the seismometer are smaller and less frequent. In this case, the earthquake threshold value can be automatically increased during the day, and the earthquake threshold value can be automatically lowered at night. However, this method can only reduce the misjudgment caused by noise, and there is still room for improvement in whether it can quickly or accurately judge as a real earthquake event.

In view of the shortcomings of the known technology, it is expected to propose a reliable detection system, detection method, and detection device for judging earthquake events, and to provide the required warnings for different objects.

The present invention improves the reliability of earthquake warning by deploying earthquake detection devices at different earthquake detection locations. It uses the triggering number of the earthquake detection device to determine whether an earthquake has occurred, which can solve the problem of unreliable warnings or false alarms.

The present invention also uses the longitudinal waves detected by the seismic detection devices at different seismic detection sites to trigger the monitoring signal, and further judges whether an earthquake has occurred based on the difference in the time point of the triggering time of the monitoring signal. At the same time, the energy of the longitudinal waves detected by the seismic detection devices at different seismic detection sites can also be used to calculate the energy difference, and further judge whether an earthquake has occurred or whether an early warning alarm should be issued for different objects based on the energy difference. This method can be applied to the accuracy and reliability when low-intensity earthquakes are required, such as precision factories that need accurate predictions of low-intensity earthquakes to shut down operations in advance to avoid product losses and increase the time cost of downtime due to false alarms. The earthquake detection system and earthquake detection device using this method are more complex, and in some cases, the deployment cost may be higher, but they are also more accurate. The present invention can also first use the majority rule to determine whether it is a potential earthquake event based on the triggering conditions of the monitoring signals detected at different locations, and then directly determine whether it is a real earthquake event based on the longitudinal wave energy difference detected at different locations.

According to the above concept, the present invention provides an earthquake detection method, comprising the following steps: receiving a first signal triggered in response to a first vibration state of a main detection location and a second signal triggered in response to a second vibration state of a first auxiliary detection location; judgment rule 1: judging whether the first signal and the second signal are both displayed as triggered at an intersection time point; judgment rule 2: judging when the intersection time point occurs, the time when the first signal and the second signal are triggered Whether the time sequence relationship between the first and second vibration points is satisfied; Judgment rule three: Determine whether an energy relationship between a first vibration intensity measurement value related to the first vibration state and a second vibration intensity measurement value related to the second vibration state does not exceed a specific value range during an energy judgment period after the intersection time point occurs; and when the time sequence relationship is satisfied and the energy relationship does not exceed the specific value range, continue to judge whether the first vibration state is a real earthquake event.

According to the above concept, the present invention provides an earthquake detection method, comprising the following steps: receiving a first signal triggered in response to a first vibration state of a main detection location and at least one second signal triggered in response to at least one second vibration state of at least one auxiliary detection location; judging rule 1: judging whether a quantity displayed as a trigger in the at least one second signal satisfies a quantity condition and Whether the first signal and the second signal satisfying the quantity condition are both displayed as triggered at an intersection time point; judgment rule 2: whether an interval between the time point of the triggering of the first signal and the time point of the triggering of one of the second signals satisfying the quantity condition is within a specific period, and when the interval is within the specific period, continue to judge whether the first vibration state is the real earthquake event. And judgment rule 3: within an energy judgment period after the intersection time point occurs, if the vibration intensity of the first vibration state and the second vibration state corresponding to the second signal for judging whether the quantity condition is satisfied satisfies an energy relationship, continue to judge whether the first vibration state is a real earthquake event. The method for determining whether the energy relationship is satisfied includes: capturing a first vibration intensity measurement value related to the first vibration state during the energy determination period, and capturing at least one second vibration intensity measurement value in the second vibration state corresponding to the second signal for determining whether the quantity condition is satisfied; obtaining an energy calculation value between the first vibration intensity measurement value and the at least one second vibration intensity measurement value; and determining whether the energy calculation value does not exceed a specific value range. The order of the steps of the second and third determination rules can be interchanged.

Based on the above concept, the present invention provides an equipment, device, system, or framework using the earthquake detection method of the present invention.

The earthquake detection method and device proposed by the present invention for judging earthquakes by the time sequence intersection of multiple seismometers can be used in areas where earthquakes can occur to reduce the damage to businesses caused by earthquakes through timely and reliable early warnings, and has industrial applicability.

Trig_M: First event trigger signal

Trig_S1: Second event trigger signal

Trig_S2: The third event trigger signal

TrigC: event judgment signal

Trig_W: Energy judgment signal

This case can be explained in detail through the following diagrams for a deeper understanding:

Figure 1 A: Schematic diagram of the monitoring signal of the earthquake detection system of an embodiment of the present invention.

Figure 1 B: Schematic diagram of the monitoring signal of the earthquake detection system of an embodiment of the present invention.

Figure 2: A schematic diagram of the process of an earthquake detection method according to an embodiment of the present invention.

Figure 3: Schematic diagram of another embodiment of the earthquake detection method of the present invention.

Figure 4: A schematic diagram of the process of another earthquake detection method of a preferred embodiment of the present invention.

Figure 5: Schematic diagram of the digital waveforms of the monitoring signal, the time determination signal, and the energy determination signal of an embodiment of the present invention.

Figure 6: A schematic diagram of the digital waveforms of the judgment formulas D and E in another preferred embodiment of the present invention.

Figure 7: A schematic diagram of a digital waveform using time judgment and energy judgment in an embodiment of the present invention.

Figure 8: Schematic diagram of vibration waves for energy determination in a preferred embodiment of the present invention.

Figure 9: A schematic diagram of another preferred embodiment of the present invention using a vibration wave for energy determination.

Figure 10: Schematic diagram of the specific process of the judgment formula E of the preferred embodiment of the present invention.

Figure 11: A schematic diagram of the specific process of determining formula F in an embodiment of the present invention.

Figure 12: A schematic diagram of a digital waveform determined as a non-seismic event by an embodiment of the present invention.

Figure 13: A schematic diagram of a digital waveform determined as a non-seismic event by an embodiment of the present invention.

Figure 14: A schematic diagram of a digital waveform determined as a non-seismic event by an embodiment of the present invention.

Figure 15: A schematic diagram of a digital waveform determined as a non-seismic event by an embodiment of the present invention.

Figure 16: A schematic diagram of a digital waveform determined as a non-seismic event by an embodiment of the present invention.

Figure 17: A schematic diagram of a digital waveform determined as a non-seismic event by an embodiment of the present invention.

Figure 18: A schematic diagram of a digital waveform determined as a non-seismic event by an embodiment of the present invention.

Figure 19: Schematic diagram of an earthquake detection method for judging earthquakes by the intersection of time sequence and energy of multiple seismometers according to an embodiment of the present invention.

Figure 20: Schematic diagram of another preferred embodiment of the earthquake detection method of the present invention.

Please refer to the attached drawings of this specification to read the following detailed description, wherein the attached drawings of this specification introduce various different embodiments of the present invention by way of example, and provide an understanding of how to implement the present invention. The embodiments of the present invention provide sufficient content for a person skilled in the art to implement the embodiments disclosed by the present invention, or to implement embodiments derived from the contents disclosed by the present invention. It should be noted that these embodiments are not mutually exclusive, and some embodiments can be appropriately combined with one or more other embodiments to form new embodiments, that is, the implementation of the present invention is not limited to the embodiments disclosed below. In addition, for the sake of simplicity and clarity, the relevant details will not be excessively disclosed in each embodiment. Even if specific details are disclosed, they are only given as examples to make the reader understand. The relevant specific details in each embodiment are not used to limit the disclosure of this case.

Please refer to the first figure A and B, which are schematic diagrams of monitoring signals related to earthquakes in the preferred embodiment of the present invention. The horizontal axis represents time, and the vertical axis represents the monitoring signals Trig_M and Trig_S1 with digital signal characteristics and when they are triggered. It should be noted that the seismometers (sensors) used in the embodiment of the present invention are configured at different locations in the same area. When an earthquake occurs, the first longitudinal wave will successively trigger each seismometer or sensor when passing through the detection location where these seismometers or sensors are configured, and then successively trigger the signal. The monitoring signals in the embodiment of the present invention can be triggered in a positive triggering or negative triggering manner, but it is not limited here.

According to one aspect of the present invention, by observing the presence or absence of these monitoring signals and comparing the time between these monitoring signals, it is possible to determine whether an earthquake has occurred. The embodiment in the first figure A is an example in which the time determination signal Trig_C is triggered, while the embodiment in the first figure B is an example in which the time determination signal Trig_C is not triggered. As can be seen from the first figure A, the interval between the triggering time point TGM of the first signal Trig_M and the triggering time point TGS1 of the second signal Trig_S1 is within the time determination period DT1, so the time determination signal Trig_C is triggered, that is, it is triggered at the triggering time point TGC to determine that the first vibration state is a type of earthquake event, that is, a potential earthquake event can be determined to be this type of earthquake event.

One of the purposes of setting the specific period DT1 is to avoid as much as possible the situation where both the first signal Trig_M of the main detection location and the second signal Trig_S1 of the auxiliary detection location are triggered but are actually non-earthquake events (for example, the main detection location and the auxiliary detection location each have unpredictable local environmental vibration factors or the sensor at the auxiliary detection location is abnormally triggered due to a malfunction, etc.), but are judged by the system as a false alarm of an earthquake event. For example, from the first figure B, it can be seen that the interval between the time point TGM when the first signal Trig_M is triggered and the time point TGS1 when the second signal Trig_S1 is triggered is not within the time determination period DT1, so the time determination signal Trig_C is not triggered, that is, if it is not triggered, it can be determined that the potential earthquake event is a non-earthquake event.

In any embodiment of the present invention, the number of monitoring signal triggers detected in real time can be used to determine whether there is a potential earthquake event according to the majority decision rule, and then enter the time difference judgment. If there are two monitoring signals including the main detection location and the auxiliary detection location, it can be reasonably set that the monitoring signals of the main detection location and the auxiliary detection location must be detected as triggered (in compliance with the majority decision, if the number of monitoring signal triggers is equal to the number of monitoring signals that are not triggered, it does not meet the majority decision rule), and then enter the subsequent judgment. According to one embodiment, if there is no monitoring signal trigger at the local end of the main detection location, it can be judged as a non-earthquake event regardless of whether it meets the majority decision rule. If there is a signal trigger at the local end of the main detection location, but the triggered signal changes to a non-triggered state during the specific period, it can also be judged as a non-earthquake event.

In any embodiment of the present invention, if the first signal Trig_M is not triggered, it is judged as a non-earthquake event. If the first signal Trig_M is triggered but does not exceed the preset time length (for example, DT1), it means that the maintenance time of the first vibration wave is insufficient. Although it can be judged whether it is an earthquake event based on the data, the estimated value of the earthquake intensity cannot be calculated for data that is less than the preset time length, so it is also judged as a non-earthquake event.

Generally speaking, earthquakes usually occur in a relatively large area, such as more than ten kilometers, and will not only appear in an area with a few buildings. The specific situation depends on the earthquake. For example, the velocity of longitudinal waves is about 6~8km/sec. Therefore, if the main detection point and the auxiliary detection point are set within a radius of 1 kilometer, the time difference between the longitudinal wave reaching the main detection point and the auxiliary detection point will not exceed 0.2 seconds. Therefore, it is practical and reliable to set the specific period in the range of 3 to 6 seconds as a judgment time for judging earthquake events. In addition, the embodiment of the present invention does not need to estimate the distance between the earthquake source and the detection site in advance, does not need to consider the transmission speed of longitudinal waves in different geological characteristics, and does not need to detect the occurrence time of multiple events of longitudinal waves. Instead, it uses the longitudinal waves detected by the main detection site and multiple auxiliary detection sites within a shorter and simultaneously accurate and effective judgment time to make a comparison, such as the time difference of detecting longitudinal waves, to shorten the judgment time.

The embodiments of the present invention in the first figure A and B are relatively simple examples of earthquake signal detection, which use a main detection point and an auxiliary detection point at different locations, and the time difference between different earthquake longitudinal waves detected by the main detection point and the auxiliary detection point to increase the reliability of earthquake judgment. However, if there are multiple auxiliary detection points, it can be further improved to more effectively judge earthquake events or non-earthquake events. The method of adding multiple auxiliary detection points can be applied to rapid detection.

In the case where the first signal Trig_M of the main detection site is triggered among multiple auxiliary detection sites, a majority decision method can be used to determine whether a potential earthquake event has occurred. For example, in an embodiment where there are only three detection sites including the main detection site, it is determined whether the first, second, and third signals Trig_M, Trig_S1, and Trig_S2 have an intersection in time. If there is an intersection, it is determined how many signals are in the time period of the intersection. When the number of signals triggered in the time period of the intersection is greater than half of the total number of monitoring signals (equivalent to at least one of the second/third signals Trig_S1/Trig_S2 being triggered), it can be determined as a potential earthquake event. The second/third signal Trig_S1/Trig_S2 can come from different auxiliary detection locations respectively.

In one embodiment of the present invention, the main detection point and a set of auxiliary detection points are respectively arranged at different detection points in the same area for a period of time. During this period of time, a plurality of monitoring signals Trig_M, Trig_S1, Trig_S2 can be firstly counted, and according to the historical data during this period of time, it can be judged whether errors or false alarms occur frequently. For example, Often, multiple monitoring signals Trig_M, Trig_S1, and Trig_S2 intersect and are judged as potential earthquake events. This may be due to the environmental vibration, which causes the seismic instruments in the primary/auxiliary detection sites to detect too much noise, which can easily lead to misjudgment. In this case, the majority decision condition can be adaptively adjusted to a full decision condition to avoid misjudgment with more stringent conditions. Or in areas with relatively quiet environments and less human-induced vibrations, a fixed number decision method can be used for judgment, rather than a majority decision method.

In an embodiment of the present invention, when the interval between the time point of the intersection signal triggering the second signal Trig_S1 and the third signal Trig_S2 (e.g., the intersection triggering time point TGS1S2 in the sixth figure) and the triggering time point TGM of the first signal Trig_M is within the time determination period DT1, a time determination signal Trig_C is triggered to determine that the first vibration state corresponding to the first signal Trig_M is a quasi-earthquake event, otherwise it is a non-earthquake event.

In an embodiment of the present invention, for example, the time determination period DT1 in the first figure A and B can be adjusted to DT1' according to the situation. In this case, the time point TGS1 triggered by the second signal Trig_S1 in the first figure A or the time point TGM triggered by the first signal Trig_M in the first figure B is not necessarily the starting point of the specific period DT1', and the time length can also be adjusted. The time determination periods DT1 and DT1' can also be appropriately adjusted so that potential earthquake events can be quickly determined.

Please refer to the second figure, which is a schematic diagram of the process of an earthquake detection method S10 of an embodiment of the present invention. The earthquake detection method S10 mainly includes a trigger quantity judgment step S101 and a time judgment step S102, and the judgment method is as follows. Step S101 (judgment rule 1) judges whether a quantity displayed as a trigger in Trig_S1 and Trig_S2 meets a quantity condition and whether Trig_M and Trig_S1 or Trig_S2 that meets the quantity condition are simultaneously displayed as a trigger at an intersection time (judgment formula D). If yes, it can be judged that a potential earthquake event is established and recorded. Then enter step S102. Step S102 (Judgment Rule 2) determines whether the interval between the time point when the intersection occurs and Trig_M is triggered and the time point when Trig_S1 or Trig_S2 that meets the quantity condition is triggered is within a time judgment period DT1 (Judgment Formula E). If yes, the time judgment signal Trig_C is triggered to judge whether the first vibration state corresponding to the first signal Trig_M is a quasi-earthquake event (or a potential earthquake event is a real earthquake event). Otherwise, proceed to step S103 to judge it as a non-earthquake event.

Please refer to the third figure, which is a schematic diagram of an earthquake detection method S20 of an embodiment of the present invention. The process of the earthquake detection method in the third figure can be summarized as follows. Step S201, in response to a first vibration state of a main detection location and a second vibration state of a first auxiliary detection location, a first signal Trig_M and a second signal Trig_S1 are triggered respectively. Step S202, receiving the first signal Trig_M and the second signal Trig_S1. Step S203 (determination rule 1): Determine whether the first signal Trig_M and the second signal Trig_S1 are both displayed as triggered at an intersection time point. Step S204 (Judgment Rule 2): When the intersection time point occurs and the interval between the time points TGM and TGS1 when the first and second signals Trig_M and Trig_S1 are triggered is within a time judgment period DT1, a time judgment signal Trig_C is triggered, so as to judge the first vibration state as a type of earthquake event.

In addition to increasing the reliability of earthquake warning by using the time difference factor between the main detection site and multiple auxiliary detection sites at different locations, the energy difference factor between the main detection site and multiple auxiliary detection sites at different locations can also be added as a judgment basis, which can further enhance the accuracy and reliability of earthquake warning. Adding energy judgment can improve the accuracy and reliability of earthquake warning. For example, precision factories need accurate predictions to shut down in advance to avoid product losses and increase downtime costs due to false alarms. Earthquake detection systems and earthquake detection devices using this method are more complex and have higher deployment costs, but they are also more accurate.

Please refer to the fourth figure, which is a schematic diagram of the process of another embodiment of the earthquake detection method S30 of the present invention. The earthquake detection method S30 in the fourth figure can continue the step S102 in the second figure and enter the step S304 of energy judgment. Step S304 (Judgment Rule 3): Determine whether Abs(I_M-I_S)≦1 within an energy judgment period after the conditions such as Trig_M=1 and the intersection triggering time point of Trig_S1 and Trig_S2 intersection is within the time judgment period DT1 are met (i.e., within the energy judgment period after the time judgment signal Trig_C is triggered). (Judgment formula F) (i.e., whether the magnitude difference is within the range of positive or negative 1), where I_M represents the magnitude value of the vibration wave associated with the first vibration state (if an earthquake event actually occurs, it may be the first seismic longitudinal wave PW1), and I_S represents the magnitude value of the vibration wave associated with the second vibration state (if the earthquake event actually occurs, it may be the second seismic longitudinal wave PW2). This embodiment may cover not only the case of receiving monitoring signals Trig_S1 and Trig_S2 from two auxiliary detection locations, but it is assumed that at least one of Trig_S1 and Trig_S2 must be in the triggering state to establish a potential earthquake event.

The judgment formula F in the fourth figure only lists a simpler calculation, and it is not necessary to use the magnitude value for calculation. For example, I_M and I_S can be the first and second vibration intensity measurement values and are respectively related to at least one of the maximum surface acceleration value, maximum velocity value, maximum displacement value, magnitude and spectral magnitude detected at the corresponding primary detection and auxiliary detection locations. The difference between the two (or addition and subtraction calculation value) can also be a proportional value, square difference and other calculation values, but is not limited to this. The proportional value calculation can be that the difference between the first and second vibration intensity measurement values is within 20% of the first vibration intensity measurement value.

In addition, the order of the steps in the second judgment rule and the third judgment rule can be interchanged.

Please refer to the fifth and sixth figures, which are schematic diagrams of the digital waveforms of the monitoring signals Trig_M, Trig_S1, Trig_S2, the time determination signal Trig_C, and the energy determination signal Trig_W of the preferred embodiment of the present invention. The horizontal axis represents time, and the vertical axis represents the monitoring signals Trig_M, Trig_S1, Trig_S2, the time determination signal Trig_C, and the energy determination signal Trig_W with digital signal characteristics. When ig_W is triggered, TGM, TGS1, TGS2, TGC, and TGW represent the triggering time points of the monitoring signals Trig_M, Trig_S1, Trig_S2, the time judgment signal Trig_C, and the energy judgment signal Trig_W respectively, and TGS1S2 represents the triggering time point of the intersection Trig_S1∩Trig_S2 of the monitoring signals Trig_S1 and Trig_S2 of the auxiliary detection site. This embodiment can cover the situation of receiving more than just the monitoring signals Trig_S1 and Trig_S2 of two auxiliary detection sites, but it is assumed that at least Trig_S1 and Trig_S2 must be in the triggered state to establish a potential earthquake event, that is, a quantitative condition is equal to 2 or greater than 2. In another embodiment, if there are other N auxiliary signals from auxiliary detection locations, the ratio of the monitoring signals Trig_S1, Trig_S2 to one of the other N-2 auxiliary signals can be set to be equal to 2/N or greater than 2/N.

In one embodiment, for example, the sampling frequency of the longitudinal wave is 200 Hz, and the longitudinal wave is sampled every 0.005 seconds. Therefore, the number of sampling points in the time determination period DT1 in the first figure A and B will be 600 points. The number of sampling points is not absolute here, and it can be adjusted appropriately, for example, with the sampling frequency of the longitudinal wave, and the sampling frequency of the longitudinal wave can be adaptively adjusted with the sophistication of the measuring instrument, or with the progress of the measurement and determination method.

In the fifth figure, the computing host (not shown) monitors all the first, second, and third signals Trig_M, Trig_S1, Trig_S2 in real time. When the second signal (first auxiliary signal) Trg_S1 is detected to be triggered at the trigger time point TGS1, the trigger quantity is only 1, and it is judged that the quantity condition (equal to 2 or greater than 2) is not met. When the computing host monitors the first signal Trig_M in real time and triggers at the trigger time point TGM, the total number of triggers of Trig_S1 and Trig_S2 is still 1, and it is judged that the quantity condition is not met. When the third signal (second auxiliary signal) Trig_S2 is detected in real time and is triggered at the trigger time point TGS2, the total number of Trig_S1 and Trig_S2 triggers is 2, and it is judged that the quantity condition is met, so it is judged as the potential earthquake event. The above embodiment can also refer to the judgment formula D of step S101 in the second figure, or steps S201~S203 in the third figure.

In one embodiment, if the number of auxiliary signal triggers can satisfy the quantity condition, then until the quantity condition can no longer be satisfied, even if a new auxiliary signal may be triggered, the triggering status of other auxiliary signals can be temporarily stopped from being monitored (taking the quantity condition as 2 as an example, when the signal Trig_S2 is triggered, the total number of auxiliary signal triggers is 2, which can satisfy the quantity condition, so the triggering status of other auxiliary signals can be temporarily stopped). In another embodiment, the triggering status of the auxiliary signal can also be continuously monitored (taking the quantity condition as greater than or equal to 2 as an example, when the signal Trig_S2 is triggered, the total number of auxiliary signal triggers is 2. Although the quantity condition has been met, the triggering status of other auxiliary signals can still be continuously monitored. If a new auxiliary signal is triggered, the number of auxiliary signal triggers is greater than 2).

Please refer to the sixth figure, which is a schematic diagram of the digital waveform of the preferred embodiment of the present invention related to the judgment formula D, E. In the sixth figure, after monitoring that the first signal Trig_M and the intersection signal Trig_S1∩Trig_S2 (or the third signal Trig_S2) of the second and third signals are all displayed as triggered at an intersection time point TGS1S2 (or TGS2), it can be judged as a potential earthquake event, and then the time judgment factor is added. If the interval between the triggering time point TGS1S2 of Trig_M and Trig_S1∩Trig_S2 (or the time point TGS2 of Trig_S2 triggering) is within the time determination period DT1, then the time determination signal Trig_C is triggered at the intersection time point TGS1S2 (or TGS2) to determine that the first vibration state corresponding to the first signal Trig_M is a quasi-earthquake event. The above-mentioned embodiment can also refer to the judgment formula D of step S101 and the judgment formula E of step S102 in the second figure, or steps S201~S204 in the third figure.

Please refer to Figure 8, which is a schematic diagram of the longitudinal wave using energy judgment in the preferred embodiment of the present invention, which is an example of a real earthquake event. The horizontal axis represents time and the vertical axis represents the surface acceleration value of the longitudinal wave.

In addition to time judgment, energy judgment can also be added. Please refer to Figures 7 and 8 together. Energy judgment is added within the 50-point sampling time (i.e., energy judgment period) after the intersection time point TGS1S2 (or TGS2) occurs or the time judgment signal Trig_C is triggered. In the embodiment of the present invention, the energy difference between the vibration waves detected by the main detection point and the auxiliary detection point is judged after 0.25 seconds, with a sampling interval of 0.005 seconds for a total of 50 times. Among them, the ground acceleration components of the vibration wave related to the first vibration state (if an earthquake event actually occurs, it can be the first seismic longitudinal wave PW1) detected by the detector (not shown) in the east-west axis, north-south axis, and vertical axis are represented by P_MX, P_MY, and P_MZ respectively, and the peak ground acceleration (peak ground acceleration, PGA) of the vibration wave related to the first vibration state in the vertical axis component is represented by PGA_MZ. The ground acceleration component of the vibration wave related to the second vibration state (if an earthquake event actually occurs, it can be the second seismic longitudinal wave PW2) in the vertical axis is represented by PGA_S1Z, and the peak ground acceleration of the vibration wave related to the second vibration state in the vertical axis component is represented by PGA_S1Z. The PGA_MZ of the vertical axial component of the vibration wave associated with the first vibration state and the PGA_S1Z of the vertical axial component of the vibration wave associated with the second vibration state are respectively converted into a first vibration intensity measurement value (e.g., seismic intensity value) and a second vibration intensity measurement value (e.g., seismic intensity value). For example, in the eighth figure, PGA_MZ is 40.5gal, which is equivalent to a magnitude of 4, and PGA_S1Z in the eighth figure is 84.4gal, which is equivalent to a magnitude of 5. The magnitude of the two magnitudes subtracted is less than or equal to the range of positive and negative magnitudes. After sampling the vibration wave to obtain PGA_MZ and PGA_S1Z during the energy judgment period DT2, the energy judgment signal Trig_W is triggered after the energy judgment period DT2 (as shown in the seventh figure), based on the difference between the two magnitudes being less than or equal to a specific value range (assuming a magnitude of 1 in this embodiment), and the earthquake event can be determined to be a true earthquake event. The energy judgment period DT2 in Figure 8 can be adjusted appropriately. The energy judgment period DT2 does not need to be related to the period of the vibration wave and is usually shorter than the aging judgment period DT1.

Please refer to Figure 9, which is a schematic diagram of another preferred embodiment of the present invention using energy judgment of vibration waves, which is an example of a non-seismic event, the horizontal axis represents time, and the vertical axis represents the ground acceleration value of the vibration wave. Please refer to Figure 9, within the 50-point sampling time after the intersection time point TGS1S2 (or TGS2) occurs or the time judgment signal Trig_C is triggered, energy judgment begins to be added. In the embodiment of the present invention, a total of 50 samples are sampled every 0.005 seconds, that is, after 0.25 seconds, the energy difference between the vibration wave of the first vibration state and the vibration wave of the second vibration state detected by the main detection point and the auxiliary detection point is judged. For example, in the ninth figure, PGA_MY is 23.1gal, which is equivalent to a magnitude of 3, and PGA_S1Z is 0.4gal, which is equivalent to a magnitude of 0. The difference between the two magnitude values is greater than the range of positive and negative magnitude values. After sampling the vibration wave to obtain PGA_MY and PGA_S1Z during the energy judgment period DT2, the energy judgment signal Trig_W will not be triggered after the energy judgment period DT2, and it can be judged as a non-earthquake event, based on the difference between the two magnitudes being greater than the specific value range (assuming a magnitude of 1 in this embodiment).

In one embodiment, if there are vibration waves from multiple auxiliary detection locations, then within an energy judgment period DT2 after the intersection time point occurs or the time judgment signal Trig_C is triggered, if the auxiliary signal related to the first vibration state (or corresponding to Trig_M) and the auxiliary signal related to judging whether the quantity condition is satisfied (for example, taking the quantity condition equal to 2 as an example, the auxiliary signal related to judging whether the quantity condition is satisfied) is equal to 2, then the auxiliary signal related to judging whether the quantity condition is satisfied is equal to 2. The auxiliary signals that meet the condition are Trig_S1 and Trig_S2; taking the condition that the quantity is greater than or equal to 2 as an example, if the auxiliary signals related to determining whether the condition is met include other auxiliary signals in addition to Trig_S1 and Trig_S2, and may also include other auxiliary signals), and the vibration intensities of the vibration states corresponding to the auxiliary signals meet an energy relationship, then it can be further determined whether the first vibration state is a real earthquake event. In one embodiment, the method for determining whether the energy relationship is satisfied includes capturing a vibration intensity meter value related to the first vibration state within the energy determination period DT2, and capturing a vibration intensity meter value related to the auxiliary signal (e.g., Trig_S1 and Trig_S2) for determining whether the quantity condition is satisfied; obtaining an energy operation value between the first vibration intensity meter value and the at least one second vibration intensity meter value; and determining whether the energy operation value does not exceed a specific value range. The vibration intensity meter value may be related to at least one of the maximum surface acceleration value, the maximum velocity value, the maximum displacement value, the magnitude, and the spectral magnitude detected at the detection site. The energy operation value may be an addition and subtraction operation, a proportional operation, a square difference operation, etc., but is not limited thereto. The proportional calculation may be that the difference between the first and second vibration intensity measurement values is within 20% of the first vibration intensity measurement value.

In any embodiment of the present invention, the main detection point only uses the maximum ground acceleration value on any axis and the maximum ground acceleration value on any axis of the auxiliary detection point for calculation, rather than fixedly using the maximum ground acceleration value on the vertical axis of the main detection point and the maximum ground acceleration value on the vertical axis of the auxiliary detection point for calculation. Because under normal circumstances, if the maximum acceleration values of the main detection point and the auxiliary detection point are converted into magnitudes due to human factors, the subtraction will be greater than the range of positive and negative magnitudes, so the problem of false alarms or misjudgments caused by human factors can be eliminated. On the contrary, when an earthquake event does occur, because longitudinal waves will transfer energy over a large range, the maximum surface acceleration value at the primary detection site and the maximum surface acceleration value at the auxiliary detection site will be very close, that is, the magnitude difference is within one level, based on which it can be judged as a real earthquake event.

Since the PGA value and maximum magnitude of the first vibration state, as well as the PGA value and maximum magnitude of the second vibration state, are added as judgment factors, the judgment of earthquake events is more accurate.

Please refer to FIG. 10, which is a schematic diagram of a specific process of another embodiment of the present invention, judgment formula E. Please refer to FIG. 2 and FIG. 10, step S102, when judging whether an interval of the time point of the intersection signal Trig_S1∩Trig_S2 of the first signal Trgi_M, the second signal (first auxiliary signal) Trgi_S1 and the third signal (second auxiliary signal) Trgi_S2 is within a time judgment period DT1, if yes, then enter step S104, judge it as a quasi-earthquake event; if not, then enter step S103, judge it as a non-earthquake event. If the first signal Trig_M is not triggered or the interval between the triggering time points of Trig_M and Trig_S1∩Trig_S2 is not within the time determination period DT1, then the process proceeds to step S103 and determines that it is a non-earthquake event.

Please refer to FIG. 11, which is a schematic diagram of the specific process of the judgment formula F of the preferred embodiment of the present invention. Please also refer to step S304 of FIG. 4 and steps S300-S303 of FIG. 11. Step S300 captures the PGA (Peak Ground Acceleration) values of the surface acceleration components P_MX, P_MY, and P_MZ related to the vibration wave of the first vibration state (if an earthquake event actually occurs, it can be the first seismic longitudinal wave PW1) during the energy judgment period DT2 after the time judgment signal Trig_C is triggered, that is, the PGA values of each P_MX, P_MY, and P_MZ, and then takes the largest PGA value among the three. Step S301, determine whether the PGA value in the energy judgment period DT2 after the time judgment signal Trig_C is triggered is ≧ a threshold value, such as 8gal. If so, proceed to the next step S302. In step S302, in addition to obtaining the PGA of the surface acceleration components P_MX, P_MY, and P_MZ of the vibration wave related to the first vibration state in the energy judgment period DT2 after the time judgment signal Trig_C is triggered in step S300, the PGA value of the surface acceleration component P_S1Z of the vibration wave related to the second vibration state (if an earthquake event actually occurs, it can be the second seismic longitudinal wave) is also obtained. In step S303, the largest PGA value among P_MX, P_MY, and P_MZ is converted into a first intensity value I_M, and the PGA value of P_S1 is converted into a second intensity value I_S. If there are multiple auxiliary detection sites and multiple vibration waves, for example, the vertical axial component P_S1Z of the vibration wave in the second vibration state and the vertical axial component P_S2Z of the vibration wave in the third vibration state (if an earthquake event actually occurs, it can be the third earthquake longitudinal wave) are detected at two auxiliary detection sites respectively, and it is assumed that at least Trig_S1 and Trig_S2 must be in the triggering state to establish a potential earthquake event (that is, the quantitative condition is equal to 2 or greater than or equal to 2), then the relevant information for determining whether the quantitative condition is The maximum PGA value of the vertical axial components P_S1Z and P_S2Z of the vibration wave of the vibration state corresponding to the satisfied auxiliary signal (for example, taking the quantity condition equal to 2 as an example, the auxiliary signals related to judging whether the quantity condition is satisfied are Trig_S1 and Trig_S2; taking the quantity condition greater than or equal to 2 as an example, the auxiliary signals related to judging whether the quantity condition is satisfied include other auxiliary signals in addition to Trig_S1 and Trig_S2, and may also include the other auxiliary signals) is converted into the second intensity value I_S. Step S304, determine whether the interval between the time points of the triggering of Trig_M and Trig_S1∩Trig_S2 is within the energy judgment period DT2 after the time judgment period DT1 (i.e., determine whether Abs(I_M-I_S) is less than or equal to C, where C represents a constant or a threshold value (determination formula F), i.e., the shock The intensity difference is within a specific value range, such as the intensity difference within the positive and negative C level), where I_M represents the intensity value of the first vibration state, and I_S represents the intensity value related to the second vibration state. If multiple auxiliary detection locations detect multiple vibration states when the potential earthquake event is triggered, I_S represents the intensity value of the vibration state of the auxiliary detection location with the largest PGA. If the judgment result of the judgment formula F in step S304 is yes, then enter step S305 and judge it as a real earthquake event; if not, enter step S103 and judge it as a non-earthquake event.

Please refer to Figures 12, 13, 14, 15, 16, and 17, which are schematic diagrams of digital waveforms of non-earthquake events determined by the preferred embodiments of the present invention. These embodiments can cover not only the situation of receiving monitoring signals Trig_S1 and Trig_S2 of two auxiliary detection locations, but it is assumed that at least Trig_S1 and Trig_S2 must be in the triggering state before continuing to determine whether it is a real earthquake event, that is, a quantitative condition is 2 or greater than or equal to 2. In FIG. 12, all the first, second, and third signals Trig_M, Trig_S1, and Trig_S2 are monitored in real time. When the second signal Trg_S1 is detected to be triggered at the trigger time point TGS1, the number of Trig_S1 and Trig_S2 triggered is only 1, and it is judged that the number condition (equal to 2 or greater than 2) is not met. When the first signal Trig_M is detected to be triggered at the trigger time point TGM, the total number of Trig_S1 and Trig_S2 triggered is still 1, and it is judged that the number condition is not met. When the real-time monitoring detects that the third signal Trig_S2 is triggered at the trigger time point TGS2, the total number of Trig_S1 and Trig_S2 triggers is 2, and it is judged that the quantity condition is met, so it can be judged as the potential earthquake event. Then, the time judgment factor is added. Since the interval between the trigger time point TGM of the first signal Trig_M and the trigger time point TGS1S2 of the intersection signal of the second signal Trig_S1 and the third signal Trig_S2 is not within the time judgment period DT1, it does not meet the condition of judgment formula E, so it is judged as a non-earthquake event. In addition, in the above embodiment, it is assumed that the time judgment step is performed before the energy judgment step. Therefore, when the relevant conditions of the time judgment step are not met, it can be judged as a non-earthquake event without the need to perform the energy judgment step (for example, under this assumption, if Trig_C is not triggered, Trig_W will not be triggered). However, in practical applications, the order of the time judgment step and the energy judgment step can be interchanged.

In FIG. 13 , all the first, second, and third signals Trig_M, Trig_S1, and Trig_S2 are monitored in real time. When the first signal Trg_M is detected to be triggered at the trigger time point TGM, both Trig_S1 and Trig_S2 are displayed as triggered. Therefore, the total number of Trig_S1 and Trig_S2 triggers is 2, and it is determined that the quantity condition is met, so it can be determined as the potential earthquake event. Then, the time judgment factor is added. Since the interval between the triggering time point TGM of the first signal Trig_M and the triggering time point TGS1S2 of the intersection signal of the second signal Trig_S1 and the third signal Trig_S2 is not within the time judgment period DT1, it does not meet the conditions of judgment formula E, so it is judged as a non-earthquake event.

In Figure 14, all the first, second, and third signals Trig_M, Trig_S1, and Trig_S2 are monitored in real time. When the second signal Trg_S1 is detected to be triggered at the trigger time point TGS1, the number of triggers is only 1, and it is judged that the quantity condition is not met. When the third signal Trig_S2 is detected to be triggered at the trigger time point TGS2, the total number of triggers is 2, and it is judged that the quantity condition is met. However, since the first signal Trig_M has never been triggered, it does not meet the condition of judgment formula D, so it is judged as a non-earthquake event.

In FIG. 15, although the number of Trig_S1 and Trig_S2 triggers is 1 at the trigger time point, since one of the third signals Trig_S2 related to the auxiliary detection location is never triggered, it is determined that the number of triggers does not meet the number condition, and does not meet the condition of judgment formula D, so it is determined to be a non-earthquake event. However, if the number condition in the embodiment of FIG. 15 is set to be equal to 1 or greater than 1, it may be the potential earthquake event. If the time judgment factor is added, since the interval between the triggering time point TGM of the first signal Trig_M and the triggering time point TGS1 of the second signal Trig_S1 is within the time judgment period DT1, the condition of judgment formula E is met, so it can be judged as a quasi-earthquake event.

In Figure 16, at the triggering time point TGM or TGS2, the number of Trig_S1 and Trig_S2 triggers is only 1, which is judged to be not meeting the quantity condition, so it is judged to be a non-earthquake event.

Please refer to Figure 17 and Figure 9 together. In Figure 17, although the time judgment signal Trig_C has been triggered, that is, the judgment formula E has been satisfied, but within the energy judgment period DT2 after the time point when the time judgment signal Trig_C is triggered, if the judgment formula F is not satisfied, that is, if Abs(I_M-I_S) is less than or equal to 1, it is judged as a non-earthquake event. For example, in the embodiment of Figure 9

Please refer to Figure 18, which is a waveform diagram of another preferred embodiment of the present invention that is judged as a non-seismic event. Although the first signal Trig_M and one of the time points of the intersection signal of the second signal Trig_S1 and the third signal Trig_S2 are triggered within the time determination period DT1, the first signal Trig_M is in a triggered state that fails to last for the time length of DT1. Even if the time point TGS1S2 of the intersection signal of the second signal Trig_S1 and the third signal Trig_S2 is triggered within the time determination period DT1, it can still be judged as a non-seismic event.

Please refer to FIG. 19, which is a schematic diagram of a method S40 for earthquake detection of a plurality of seismographs for determining earthquakes by intersection of timing and energy according to a preferred embodiment of the present invention. Step S401: In response to a first vibration state of a primary detection location and a second vibration state of a first auxiliary detection location, a first signal and a second signal are triggered respectively. Step S402: Receive the first signal and the second signal of the auxiliary detection location. Step S403 (determination rule 1): Determine whether the first and second signals are both displayed as triggered at an intersection time point. Step S404 (Judgment Rule 2): Judgment whether a timing relationship is satisfied between the time points at which the first signal and the second signal are triggered when the intersection time point occurs. Step S405 (Judgment Rule 3): Judgment whether an energy calculation value between a first vibration intensity measurement value related to the first vibration state and a second vibration intensity measurement value related to the second vibration state does not exceed a specific value range during an energy judgment period after the intersection time point occurs. S406: When the timing relationship is satisfied and the energy calculation value does not exceed the specific value range, continue to judge whether the first vibration state is a real earthquake event. Among them, the order of steps in the second judgment rule and the third judgment rule can be interchanged.

In any embodiment of the present invention, the time relationship is that the time point at which the first signal Trig_M is triggered and the time point at which the second signal Trig_S1 is triggered are separated by a specific period DT1.

In any embodiment of the present invention, the earthquake detection method S40 further includes receiving a third signal Trig_S2 triggered in response to a third vibration state of a second auxiliary detection location, wherein the time point of triggering the third signal Trig_S2 may be earlier than or later than the time point of triggering the second signal Trig_S1.

In any embodiment of the present invention, at the intersection time point, the second signal Trig_S1 and the third signal Trig_S2 are both displayed as triggered, and the earthquake detection method S40 further includes capturing a first, a second, and a third maximum vibration intensity measurement value related to the first, second, and third vibration states, respectively, in the energy judgment period DT2, taking the maximum value of the second and third maximum vibration intensity measurement values as a second operation value I_S, the first maximum vibration intensity measurement value as a first operation value I_M, and taking the difference between the first operation value I_M and the second operation value I_S as the energy operation value.

Please refer to FIG. 20, which is a schematic diagram of another preferred embodiment of the present invention, an earthquake detection method S70. The method S70 comprises: Step S701: In response to a first vibration state of a primary detection location and at least a second vibration state of at least one auxiliary detection location, a first signal Trig_M and at least a second signal Trig_S1, Trig_S2 are triggered respectively. Step S702: Receive the first signal Trig_M and the at least a second signal Trig_S1, Trig_S2. Step S703 (determination rule 1): Determine whether a quantity displayed as a trigger in the at least one second signal Trig_S1, Trig_S2 satisfies a quantity condition (for example, equal to 2 or greater than or equal to 2) and whether the first signal Trig_M and the second signal (for example, Trig_S2) that meets the quantity condition are both displayed as triggers at an intersection time point (for example, TGS1S2 or TGS2). Step S704 (Judgment Rule 3): Within an energy judgment period DT2 after the intersection time point occurs, if an energy relationship is satisfied between the vibration intensity of the second vibration state corresponding to the second signal Trig_S1 and Trig_S2 related to the first vibration state and the second vibration state corresponding to the second signal Trig_S1 and Trig_S2 related to judging whether the quantity condition is satisfied, continue to judge whether the first vibration state is a real earthquake event. In addition, any embodiment of the present invention may also include: judgment rule 2: judging whether an interval between the time point TGM at which the first signal Trig_M is triggered and the time point (e.g., TGS2) at which one of the second signals that meets the quantity condition (e.g., Trig_S2) is triggered is within a specific period DT1, wherein the step sequence of the judgment rule 2 and the judgment rule 3 can be interchanged; and when the interval is within the specific period, continue to judge whether the first vibration state is the real earthquake event.

Any embodiment of the present invention may include the following steps: providing an intersection signal (e.g., Trig_S1∩Trig_S2), triggering the intersection signal when the quantity meets the quantity condition; when the first signal Trig_M and the intersection signal Trig_S1∩Trig_S2 are both displayed as triggered, determining that the intersection time point occurs (e.g., TGS1S2 or TGS2). In addition, the at least one second signal has a first quantity displayed as triggered and a second quantity displayed as untriggered, and the quantity condition may be that the first quantity is greater than or equal to the second quantity.

Any embodiment of the present invention may include the following steps: providing an energy determination signal Trig_W, wherein: the energy determination signal Trig_W is triggered when the energy relationship is satisfied during the energy determination period DT2.

In any embodiment of the present invention, the method for determining whether the energy relationship is satisfied includes: capturing a first vibration intensity measurement value (e.g., first intensity value I_M) related to the first vibration state within the energy determination period DT2, and capturing at least one second vibration intensity measurement value (e.g., second intensity value I_S) in the second vibration state corresponding to the second signal for determining whether the quantity condition is satisfied, wherein the vibration intensity measurement values are corresponding to the detection location Abs(I_M-I_S) is less than or equal to 1).

In any embodiment of the present invention, the auxiliary detection location includes a first auxiliary detection location and a second auxiliary detection location, and the at least one second signal includes a first auxiliary signal Trig_S1 related to a first vibration state and a second auxiliary signal Trig_S2 related to a second vibration state and comes from the first auxiliary detection location and the second auxiliary detection location respectively, and the method further includes: receiving the first auxiliary signal Trig_S1; receiving the second auxiliary signal Trig_S2, wherein the first auxiliary signal The two auxiliary signals Trig_S2 are transmitted via a remote network; the quantity condition is greater than or equal to 2; at the intersection time point, the first auxiliary signal Trig_S1 and the second auxiliary signal Trig_S2 are both displayed as triggered; and the time relationship is that the time point TGM of the triggering of the first signal Trig_M and the later time point (for example, TGS2) of the triggering of the first auxiliary signal Trig_S1 and the second auxiliary signal Trig_S2 are separated within a specific period DT1. If the time relationship is not satisfied, it is determined to be a non-earthquake event. In addition, during the energy judgment period DT2, a first maximum magnitude value I_M of the first vibration state and a second maximum magnitude value I_S between the first auxiliary and second auxiliary vibration states can be captured, the second maximum magnitude value is taken as a second calculation value, and the first maximum magnitude value is taken as a first calculation value. When the difference between the first calculation value and the second calculation value does not exceed the specific value range (for example, Abs(I_M-I_S) is less than or equal to 1), the earthquake event is judged to be a real earthquake event.

Professionals in this field can understand that the various embodiments of the present invention described above can be implemented using electronic information equipment, devices, systems, or architectures, such as servers or computers, with appropriate software, hardware, or firmware.

Many of the variations and other embodiments of the present disclosure set forth herein will be understood by those skilled in the art with the benefit of the teachings presented in the above description and the associated drawings. Therefore, it should be understood that the present disclosure is not limited to the specific embodiments disclosed, and variations and other embodiments are intended to be included within the scope of the following patent application.

S40: Method for judging earthquakes by the intersection of time sequence and energy of multiple seismometers

S401~S406: Steps of the method for judging earthquakes by the intersection of time sequence and energy of multiple seismometers

Claims (10)

A method for earthquake detection by using the intersection of the time sequence and energy of multiple seismometers to determine earthquakes includes the following steps: receiving a first signal triggered in response to a first vibration state of a primary detection location and a second signal triggered in response to a second vibration state of a first auxiliary detection location; judging rule 1: judging whether the first signal and the second signal are both displayed as triggered at an intersection time point; judging rule 2: judging whether the first signal and the second signal are both displayed as triggered at an intersection time point; judging rule 3: judging whether the first signal and the second signal are both displayed as triggered at an intersection time point; judging rule 4: judging whether the first signal and the second signal are both displayed as triggered at an intersection time point; judging rule 5: judging whether the first signal and the second signal are both displayed as triggered at an intersection time point; judging rule 6: judging whether the first signal and the second signal are both displayed as triggered at an intersection time point; judging rule 7: judging whether the first signal and the second signal are both displayed as triggered at an intersection time point; judging rule 8: judging whether the first signal and the second signal are both displayed as triggered at an intersection time point; judging rule 9: judging whether the first signal and the second signal are both displayed as triggered at an intersection time point; judging rule 10: judging whether the first signal and the second signal are both displayed as triggered at an intersection time point; judging rule 11: judging whether the first signal and the second signal are both displayed as triggered at an intersection time point; judging rule 12: judging whether the first signal and the second signal are both displayed as triggered at an intersection time point; judging rule 13: judging whether the first signal and the second signal are both displayed as triggered at an intersection time point; judging rule 14: judging whether the first signal and the second signal are both displayed as triggered at an intersection time point; judging rule 15: judging whether the first signal and the second signal are both displayed as triggered at an intersection time point; judging rule 16: judging whether the first signal and the second signal are both displayed as triggered at an intersection time point; judging rule 17: judging whether the first signal and the second signal are both displayed as triggered at an intersection time point; judging rule 18: judging whether the first signal and Whether the time points of the triggering satisfy a time sequence relationship; Judgment rule three: Determine whether an energy relationship between a first vibration intensity measurement value related to the first vibration state and a second vibration intensity measurement value related to the second vibration state does not exceed a specific range during an energy judgment period after the intersection time point occurs; and when the time sequence relationship is satisfied and the energy relationship does not exceed the specific range, continue to judge whether the first vibration state is a real earthquake event. The earthquake detection method as described in claim 1, wherein the timing relationship is that the time point at which the first signal is triggered and the time point at which the second signal is triggered are separated by a specific period. The earthquake detection method as described in claim 1, wherein the step sequence of the judgment rule 2 and the judgment rule 3 can be interchanged. An earthquake detection method as described in claim 1, wherein: the energy relationship is one of an addition-subtraction relationship, a proportional relationship, and a square difference relationship between the first and second vibration intensity measurement values; and the first and second vibration intensity measurement values are related to at least one of the maximum surface acceleration value, maximum velocity value, maximum displacement value, seismic intensity, and spectral seismic intensity detected at the corresponding primary detection and auxiliary detection locations. A method for earthquake detection comprises the following steps: receiving a first signal triggered in response to a first vibration state of a main detection location and at least one second signal triggered in response to at least one second vibration state of at least one auxiliary detection location; judging rule 1: judging whether a quantity indicated as a trigger in the at least one second signal satisfies a quantity condition and whether the first signal and the quantity satisfying the quantity condition are equal. Whether the second signal is displayed as a trigger at an intersection time point; and judgment rule 2: within an energy judgment period after the intersection time point occurs, if an energy relationship is satisfied between the vibration intensity of the second vibration state corresponding to the second signal related to the first vibration state and the second vibration state corresponding to the second signal related to judging whether the quantity condition is satisfied, continue to judge whether the first vibration state is a real earthquake event. The earthquake detection method as described in claim 5, wherein the method for determining whether the energy relationship is satisfied includes: capturing a first vibration intensity measurement value related to the first vibration state during the energy determination period, and capturing at least one second vibration intensity measurement value in the second vibration state corresponding to the second signal for determining whether the quantity condition is satisfied, wherein the vibration intensity measurement values are related to at least one of the maximum surface acceleration value, maximum velocity value, maximum displacement value, magnitude and spectral magnitude detected at the corresponding detection location; obtaining an energy calculation value between the first vibration intensity measurement value and the at least one second vibration intensity measurement value; and determining whether the energy calculation value does not exceed a specific range. The earthquake detection method as described in claim 6 further comprises the following steps: providing an intersection signal, wherein: the intersection signal is triggered when the quantity meets the quantity condition; when the first signal and the intersection signal are both displayed as triggered, it is determined that the intersection time point occurs; and the at least one second signal has a first quantity displayed as triggered and a second quantity displayed as untriggered, and the quantity condition is that the first quantity is greater than or equal to the second quantity. The earthquake detection method as described in claim 6 further comprises the following steps: providing an energy judgment signal, wherein: the energy judgment signal is triggered when the energy relationship is satisfied during the energy judgment period; and the energy calculation value is one of the addition and subtraction operations, proportional operations, and square difference operations between the first and second vibration intensity measurement values. The earthquake detection method as described in claim 5 further includes: judgment rule 2: judging whether an interval between the time point of the triggering of the first signal and the time point of the triggering of one of the second signals satisfying the quantity condition is within a specific period, wherein the step sequence of the judgment rule 2 and the judgment rule 3 can be interchanged; and when the interval is within the specific period, continuing to judge whether the first vibration state is the real earthquake event. An apparatus, device, system, or structure for an earthquake detection method for determining earthquakes by using the time sequence and energy intersection of multiple seismometers as in claim 1-9.

Images