CN-121982869-A - Method for acquiring MEMS vibration data and receiving and releasing early warning information
Abstract
The application discloses a method for collecting MEMS vibration data, receiving and releasing early warning information, wherein each MEMS monitoring point is a miniature early warning information release node, so that a release network with the same density and wide coverage as a sensor network is formed, and the cost is effectively reduced; and the earthquake early warning center server makes a decision, and the early warning is issued to the monitoring point to the terminal user equipment within the preset range, so that an innovative framework of cloud decision and edge issue is constructed, the long path of early warning information from the center server to the user through the operator network is omitted, the end-to-end time from the decision to the touch is greatly shortened, and the most precious seconds or even sub-seconds is strived for the near-earthquake area. Even after the communication with the central server is interrupted, the authorized monitoring points can still act as independent early warning beacons, so that early warning autonomy of a specific area is realized.
Inventors
- WEI YONGXIANG
- WU CANJIN
- CHEN JIALIANG
- ZHOU YUEYONG
Assignees
- 厦门帝嘉科技股份有限公司
- 福建省地震局
Dates
- Publication Date
- 20260505
- Application Date
- 20260313
Claims (10)
- 1. The method for collecting and receiving and releasing the MEMS vibration data and the early warning information is characterized by comprising the following steps: Acquiring triaxial acceleration characteristic parameters of the earth surface through an MEMS vibration sensor arranged at least one monitoring point in a monitoring area, wherein the triaxial acceleration characteristic parameters comprise peak ground acceleration characteristic parameters, P-wave characteristic period characteristic parameters, accumulated absolute speed characteristic parameters, vertical and horizontal peak value characteristic parameters; Calculating the triaxial acceleration characteristic parameters in real time to obtain comprehensive characteristic parameters; When the comprehensive characteristic parameters are monitored to exceed a first preset trigger threshold, judging a suspected seismic event and generating seismic characteristic information, wherein the seismic characteristic information comprises coordinate information of monitoring points, triaxial acceleration characteristic parameters and comprehensive characteristic parameters; the seismic characteristic information is sent to a seismic early warning center server through a wireless communication module; The earthquake early warning center server generates earthquake early warning information according to the earthquake characteristic information of the preset number of monitoring points, wherein the earthquake early warning information comprises an earthquake early warning area, an estimated earthquake magnitude and an estimated earthquake intensity; The earthquake early warning center server issues early warning to the monitoring points in the earthquake early warning area according to the earthquake early warning information; and the monitoring point issues earthquake early warning to the terminal user equipment within a preset range.
- 2. The method for acquiring and issuing MEMS vibration data and early warning information according to claim 1, wherein the method for acquiring the triaxial acceleration characteristic parameters of the earth surface by the MEMS vibration sensor arranged at least one monitoring point in the monitoring area is as follows: Calculating the absolute value of the triaxial acceleration synthesis vector, and taking the maximum value in a preset time window as a peak ground acceleration characteristic parameter; In a P wave triggering time window, carrying out numerical integration on the triaxial acceleration time sequence data to obtain speed and displacement data, and solving the P wave characteristic period characteristic parameters according to a preset P wave characteristic period calculation formula; numerical integration is carried out on the absolute value of the triaxial acceleration synthesis vector in a preset time window, and an accumulated absolute speed characteristic parameter is obtained; and respectively calculating a vertical acceleration maximum value and a horizontal composite acceleration maximum value, and taking the ratio of the vertical acceleration maximum value to the horizontal composite acceleration maximum value as a vertical peak value and a horizontal peak value identification parameter.
- 3. The method for acquiring and receiving and issuing MEMS vibration data and early warning information according to claim 2, wherein the method for calculating the absolute value of the triaxial acceleration synthesis vector and taking the maximum value in a preset time window as the peak ground acceleration characteristic parameter is as follows: calculating a composite acceleration according to the formula ai=sqrt (Xi 2+yi2+zi 2), xi being an X-axis acceleration, yi being a Y-axis acceleration, zi being a Z-axis acceleration, ai being a composite acceleration; and searching the maximum value of Ai as a peak ground acceleration characteristic parameter in a preset time period.
- 4. The method for acquiring and distributing MEMS vibration data and early warning information according to claim 3, wherein the method comprises the following steps: The calculation formula of the P-wave characteristic period is that τc=2pi Sqrt ([ mean of Dx 2+mean of Dy 2+mean of Dz2 ]/[ mean of Vx 2+mean of Vy 2+mean of Vz2 ]), dx is X-axis displacement data, dy is Y-axis displacement data, dz is Z-axis displacement data, vx is velocity data obtained by integrating X-axis acceleration once, vy is velocity data obtained by integrating Y-axis acceleration once, and Vz is velocity data obtained by integrating Z-axis acceleration once.
- 5. The method for acquiring and receiving and issuing the early warning information with the MEMS vibration data according to claim 3, wherein the method for respectively calculating the maximum value of the vertical acceleration and the maximum value of the horizontal composite acceleration and taking the ratio of the maximum value of the vertical acceleration to the maximum value of the horizontal composite acceleration as the characteristic parameter of the ratio of the vertical peak value to the horizontal peak value is as follows: Acquiring the maximum value of the absolute value of the Z-axis acceleration; Calculating the horizontal combined acceleration through a formula hi=sqrt (Xi 2+yi 2), and obtaining the maximum value of the horizontal combined acceleration; and taking the ratio of the maximum value of the vertical acceleration to the maximum value of the horizontal composite acceleration as the vertical and horizontal peak value bit parameters.
- 6. The method for acquiring and issuing MEMS vibration data and early warning information according to claim 5, wherein the method for calculating the triaxial acceleration characteristic parameters in real time to obtain the comprehensive characteristic parameters is as follows: setting a scoring function and weight for each characteristic parameter of the peak ground acceleration characteristic parameter, the P wave characteristic period characteristic parameter, the accumulated absolute speed characteristic parameter, the vertical peak value characteristic parameter and the horizontal peak value characteristic parameter; Obtaining corresponding scores through the scoring function according to the peak ground acceleration characteristic parameters, the P wave characteristic period characteristic parameters, the accumulated absolute speed characteristic parameters and the vertical and horizontal peak value characteristic parameters which are calculated in real time; and multiplying the corresponding fractions obtained by the peak ground acceleration characteristic parameter, the P wave characteristic period characteristic parameter, the accumulated absolute speed characteristic parameter and the vertical and horizontal peak value characteristic parameters by weights, and then summing to obtain the comprehensive characteristic parameter.
- 7. The method for collecting and distributing MEMS vibration data according to claim 6, wherein the calculation formula of the comprehensive characteristic parameter is: T=(W 1 F 1 (PGA))+(W 2 F 2 (τc))+(W 3 F 3 (CAV))+(W 4 F 4 (V/H)),W 1 ,W 2 ,W 3 ,W 4 , wherein the weight sum of each characteristic parameter is 1;F 1 (PGA), F 2 (τc), F 3 (CAV), F 4 (V/H), and F 4 (V/H).
- 8. The method for acquiring and receiving and releasing the early warning information of the MEMS vibration data according to claim 7, wherein the method for generating the early warning information of the earthquake by the earthquake early warning center server according to the earthquake feature information of the preset number of monitoring points is as follows: Based on a space consistency principle, when the number of triggered monitoring points reaches a preset threshold value in a preset time and geographic window, confirming that an effective earthquake event occurs; Based on the converged seismic characteristic information, performing seismic positioning and magnitude estimation to respectively obtain a seism position and an estimated magnitude; calculating the estimated earthquake intensity of the peripheral area by utilizing an earthquake attenuation model and a site amplification effect model according to the estimated earthquake magnitude and the earthquake center position, and defining an area of which the estimated earthquake intensity reaches or exceeds a preset alarm threshold value as an earthquake early warning area; And integrating the earthquake early-warning area, the estimated magnitude and the estimated earthquake intensity information to generate a structured earthquake early-warning information data packet.
- 9. The method for acquiring and releasing the early warning information and the MEMS vibration data according to claim 8, wherein the method for performing earthquake positioning and earthquake level estimation based on the converged earthquake characteristic information to respectively obtain the earthquake center position and the estimated earthquake level comprises the following steps: Longitude according to the formula \u seism=Σ (longitude_i) W_i)/ΣW_i and latitude seismometer = Σ (latitude_i Performing seismic localization by using W_i)/sigma W_i, wherein the weight W_i is a P-wave characteristic period characteristic parameter reported by each monitoring point; determining the maximum value of the P-wave characteristic period characteristic parameters according to the P-wave characteristic period characteristic parameters reported by each monitoring point; According to the formula m=a The estimated magnitude is calculated by log10 (τc_max) +b, wherein M is the estimated magnitude, τc_max is the maximum value of the characteristic periodic parameters of the P wave, and a and b are constants.
- 10. The method of claim 9, wherein said estimating the magnitude and mid-shock location based on the MEMS vibration data acquisition, pre-warning information reception and distribution, the method for calculating the estimated earthquake intensity of the peripheral area by utilizing the earthquake attenuation model and the site amplification effect model and delineating the area of which the estimated earthquake intensity reaches or exceeds a preset alarm threshold as an earthquake early warning area comprises the following steps: calculating estimated earth surface peak ground acceleration of each grid point in the early warning area by using a preset earthquake damping relation model and a field amplification effect model based on the estimated earthquake magnitude and the earthquake center position; Converting the estimated peak ground acceleration of each grid point into estimated seismic intensity according to a preset instrument intensity conversion relation; comparing the estimated earthquake intensity with a preset minimum alarm intensity threshold value, and delineating a continuous geographic area formed by grid points with the estimated earthquake intensity larger than the minimum alarm intensity threshold value as an earthquake early warning area of the earthquake.
Description
Method for acquiring MEMS vibration data and receiving and releasing early warning information Technical Field The invention relates to the technical field of earthquake early warning, in particular to a method for acquiring MEMS vibration data, receiving and releasing early warning information. Background Earthquake early warning is a key technology for providing emergency risk avoidance time for a specific area before destructive earthquake waves arrive. The existing earthquake early warning is mainly based on an official station network system constructed by an expensive professional seismograph, and although the professional station network system has accurate data, the construction cost is extremely high, so that the station density is insufficient, and the early warning blind area is large. More importantly, after the early warning information is generated in the central server, the early warning information is pushed to the mobile phone of the user through a third party network such as a telecom operator. This process involves data aggregation, processing, decision making and multi-level distribution, with long internal delays, and in extreme disaster scenarios, the communication network may be congested or interrupted, resulting in the inability of the early warning information to reach the user effectively. Disclosure of Invention The application provides a method for acquiring MEMS vibration data and receiving and releasing early warning information, which solves the problems of long early warning information transmission link, high delay and unreliable communication link under extreme conditions caused by adopting a centralized decision and remote release mode in the existing early warning system. In order to solve the technical problems, the application provides a method for acquiring MEMS vibration data and receiving and releasing early warning information, which comprises the following steps: Acquiring triaxial acceleration characteristic parameters of the earth surface through an MEMS vibration sensor arranged at least one monitoring point in a monitoring area, wherein the triaxial acceleration characteristic parameters comprise peak ground acceleration characteristic parameters, P-wave characteristic period characteristic parameters, accumulated absolute speed characteristic parameters, vertical and horizontal peak value characteristic parameters; Calculating the triaxial acceleration characteristic parameters in real time to obtain comprehensive characteristic parameters; When the comprehensive characteristic parameters are monitored to exceed a first preset trigger threshold, judging a suspected seismic event and generating seismic characteristic information, wherein the seismic characteristic information comprises coordinate information of monitoring points, triaxial acceleration characteristic parameters and comprehensive characteristic parameters; the seismic characteristic information is sent to a seismic early warning center server through a wireless communication module; The earthquake early warning center server generates earthquake early warning information according to the earthquake characteristic information of the preset number of monitoring points, wherein the earthquake early warning information comprises an earthquake early warning area, an estimated earthquake magnitude and an estimated earthquake intensity; The earthquake early warning center server issues early warning to the monitoring points in the earthquake early warning area according to the earthquake early warning information; and the monitoring point issues earthquake early warning to the terminal user equipment within a preset range. In one embodiment, the method for acquiring the triaxial acceleration characteristic parameters of the earth surface by using the MEMS vibration sensor arranged at least one monitoring point in the monitoring area comprises the following steps: Calculating the absolute value of the triaxial acceleration synthesis vector, and taking the maximum value in a preset time window as a peak ground acceleration characteristic parameter; In a P wave triggering time window, carrying out numerical integration on the triaxial acceleration time sequence data to obtain speed and displacement data, and solving the P wave characteristic period characteristic parameters according to a preset P wave characteristic period calculation formula; numerical integration is carried out on the absolute value of the triaxial acceleration synthesis vector in a preset time window, and an accumulated absolute speed characteristic parameter is obtained; and respectively calculating a vertical acceleration maximum value and a horizontal composite acceleration maximum value, and taking the ratio of the vertical acceleration maximum value to the horizontal composite acceleration maximum value as a vertical peak value and a horizontal peak value identification parameter. In one embodiment, the method for calculating the