CN109061719B - Method for predicting earthquake - Google Patents

Method for predicting earthquake Download PDF

Info

Publication number
CN109061719B
CN109061719B CN201810731412.0A CN201810731412A CN109061719B CN 109061719 B CN109061719 B CN 109061719B CN 201810731412 A CN201810731412 A CN 201810731412A CN 109061719 B CN109061719 B CN 109061719B
Authority
CN
China
Prior art keywords
earthquake
monitoring
monitoring points
taking
trough
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810731412.0A
Other languages
Chinese (zh)
Other versions
CN109061719A (en
Inventor
张崇耀
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to CN201810731412.0A priority Critical patent/CN109061719B/en
Publication of CN109061719A publication Critical patent/CN109061719A/en
Application granted granted Critical
Publication of CN109061719B publication Critical patent/CN109061719B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/01Measuring or predicting earthquakes

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

A method for predicting an earthquake, comprising the steps of: arranging monitoring points on a map at intervals of 100km in a determinant manner, monitoring parameter change of the monitoring points before an earthquake in a fixed time period, drawing a parameter change curve by taking time as an abscissa, comparing each peak and each valley in the curve, and recording peaks and valleys which continuously appear in the curve and are similar in shape; recording a connecting line of two monitoring points with the largest peak or trough area and the smallest area as m, making a parallel line n of m through the other monitoring point, and taking the other three monitoring points as e and f; the straight lines m, n, e and f are intersected to form a parallelogram, and the position of the parallelogram is the epicenter position; the number of wave crests or wave troughs in the curve is marked as N, and N +4 is the magnitude of the epicenter region. The method can predict the epicenter position, the onset time and the magnitude of the earthquake between the first unit time and the last six unit times or between the last nine unit times and the last sixteen unit times when the last wave crest or wave trough is finished.

Description

Method for predicting earthquake
Technical Field
The invention relates to a method for predicting earthquake.
Background
Earthquake is also called earthquake and earth vibration, and is a natural phenomenon that earthquake waves are generated during the vibration caused in the process of quickly releasing energy from the earth crust.
The direct earthquake disaster is a disaster caused by the primary phenomenon of earthquake, such as earthquake fault dislocation and ground vibration caused by earthquake waves. Mainly comprises the following steps: destruction of the ground, destruction of buildings and structures, destruction of natural objects such as mountains (such as landslides and debris flows), tsunamis, and burnings of the ground.
The current earthquake prediction mode can only predict approximate time and azimuth, and cannot accurately predict the earthquake starting time, the epicenter position and the magnitude.
Disclosure of Invention
The technical problem to be solved by the invention is to provide a method capable of predicting the epicenter position, the onset time and the magnitude.
In order to solve the technical problem, the invention provides a method for predicting earthquake, which comprises the following steps:
(1) the method comprises the following steps of establishing a planar rectangular coordinate system on a map by taking the center of a region to be monitored as an origin, the direction from west to east as an x-axis forward direction and the direction from south to north as a y-axis forward direction, arranging monitoring points in a determinant manner by taking the origin as the center and taking the x-axis direction and the y-axis direction as intervals, monitoring changes of physical quantities before an earthquake by the monitoring points in a fixed time period, and drawing a physical quantity curve by taking time as a horizontal coordinate, wherein the monitored physical quantities comprise: atmospheric pressure, air temperature, ground stress, earth electricity, geomagnetism, altitude, gravity acceleration, water level, water temperature, radon content in water, carbon dioxide content in water, nitrite ion content in water and rainfall;
(2) comparing the shapes of each peak and each trough in the physical quantity curve drawn in the step (1), and recording the peaks and the troughs which continuously appear in each physical quantity curve and are similar in shape;
(3) taking three monitoring points near an area to be detected, comparing the monitoring results of the monitoring points on the same physical quantity in the same time period, marking the connection line of the two monitoring points with the maximum peak or trough area and the minimum area on a map as a straight line m, and then marking the parallel line passing through the other monitoring point as the straight line m as a straight line n;
(4) taking the other three monitoring points near the area to be detected, recording the connecting line of the two monitoring points with the maximum peak or trough area and the minimum area on the map as a straight line e for the monitoring results of different physical quantities in the same time period in the step (3), and then recording the parallel line passing through the other monitoring point to be the straight line e as a straight line f;
(5) the straight lines m, n, e and f are intersected to form a parallelogram, and the position of the parallelogram is the epicenter position;
(6) recording the physical quantities monitored by the monitoring points in the step (3) and the step (4), recording the number of wave crests or wave troughs in a physical quantity curve as N, wherein N +4 is the magnitude of the earthquake in the earthquake area;
(7) the origin time interval: and the unit time is between the unit time before the last peak or trough is finished and six unit times later, or between the nine unit times after the last peak or trough is finished and sixteen unit times later, wherein the unit time is the division value of the abscissa of the corresponding physical quantity curve.
Preferably, the value range of N is more than or equal to 1 and less than or equal to 6, wherein N is a positive integer.
For simplicity of explanation, the method for predicting earthquake according to the present invention is simply referred to as the present method below.
The method has the advantages that: the method can accurately predict the epicenter position, the onset time and the magnitude of the earthquake without expensive instruments, strives for time for people to evacuate from the earthquake area and rescue substances, has simple prediction mode and no need of large amount of calculation, can adopt an unattended mode at each monitoring point, has low labor cost and is suitable for popularization and use.
Drawings
FIG. 1 is a line graph of acceleration of gravity at monitoring point A, B, C over time.
Fig. 2 is a line graph of the water level at monitoring point D, E, F over time.
FIG. 3 is a epicenter position indication diagram.
Fig. 4 shows a grid pattern for recording peaks and valleys of the physical quantity curve.
FIG. 5 is a pressure day-to-day line from 2 months to 3 months from a Paschen stage earthquake in 1966.
Detailed Description
Referring to fig. 1-3, a method of predicting an earthquake, comprising the steps of:
(1) the method comprises the following steps of establishing a planar rectangular coordinate system on a map by taking the center of a region to be monitored as an origin, the direction from west to east as an x-axis forward direction and the direction from south to north as a y-axis forward direction, arranging monitoring points in a determinant manner by taking the origin as the center and taking the x-axis direction and the y-axis direction as intervals, monitoring changes of physical quantities before an earthquake by taking one day as a time period and drawing a physical quantity change curve by taking time as a horizontal coordinate, wherein the monitored physical quantities comprise: atmospheric pressure, air temperature, ground stress, earth electricity, geomagnetism, altitude, gravity acceleration, water level, water temperature, radon content in water, carbon dioxide content in water, nitrite ion content in water and rainfall;
(2) comparing the shapes of each peak and each trough in the curve drawn in the step (1), and recording the peaks and the troughs which continuously appear in each curve and are similar in shape;
(3) taking three monitoring points near the area to be detected, respectively recording as a monitoring point A, a monitoring point B and a monitoring point C, and comparing the gravity acceleration (the unit of ordinate is m/s) of each monitoring point to the gravity acceleration in the same time period-2) According to the monitoring result, a connecting line of two monitoring points with the maximum peak area (monitoring point A) and the minimum peak area (monitoring point B) on a map is marked as a straight line m, and then a parallel line passing through the other monitoring point (monitoring point C) and making the straight line m is marked as a straight line n;
(4) taking the other three monitoring points near the area to be detected, respectively recording the other three monitoring points as a monitoring point D, a monitoring point E and a monitoring point F, recording the connection line of the two monitoring points with the maximum peak area (the monitoring point D) and the minimum peak area (the monitoring point E) on the map as a straight line E for the monitoring result of the water level (the unit of the vertical coordinate is m) in the same time period in the step (3), and then, taking the parallel line of the straight line E passing through the other monitoring point (the monitoring point F) as a straight line F;
(5) the straight lines m, n, e and f are intersected to form a parallelogram, and the position of the parallelogram is the epicenter position;
(6) recording the number of wave crests in the physical quantity curve monitored by the monitoring points in the step (3) and the step (4), wherein the number of the wave crests is 4, and the seismic level of the obtained epicenter area is 8;
(7) the origin time interval: the time between one day before the end of the last wave crest or wave trough and six days after the end of the last wave crest or wave trough, or the time between nine days after the end of the last wave crest or wave trough and sixteen days after the end of the last wave crest or wave trough.
Specifically, the shape of the continuous peaks or valleys in the recorded physical quantity curve may be recorded using the M-table grid in fig. 4.
For the first peak or trough, in the region of the number a 1; recording in the region marked with the number b1 on the upper side of the first time if the position of the second-order peak or trough is higher than the first time, in the region marked with the number b1 on the upper side of the first time if the second-order peak or trough is substantially equal to the first time, and in the region marked with the number b1 on the lower side of the first time if the second-order peak or trough is lower than the first time; the region labeled with the number c1 on the upper side of the second time if the third-order peak or trough position is higher than the second time, the region labeled with the number c1 on the upper side of the second time if the third-order peak or trough position is substantially equal to the second time, the region labeled with the number c1 on the lower side of the second time if the third-order peak or trough position is lower than the second time, and so on.
Example verification:
FIG. 5 shows the average line of the gas pressure day from 2 months to 3 months of the Paschen platform earthquake in 1966, three peaks are found, and the ending time of the last peak is 3 months and 5 days, so according to the calculation of the method, the earthquake grade is 7, the earthquake starting time is 3 months and 4 days to 3 months and 11 days or 3 months and 14 days to 3 months and 21 days, the actual earthquake starting time is 3 months and 8 days, and the earthquake grade is 6.8 grades, which is consistent with the prediction.

Claims (2)

1. A method of predicting an earthquake, comprising the steps of:
(1) the method comprises the following steps of establishing a planar rectangular coordinate system on a map by taking the center of a region to be monitored as an origin, the direction from west to east as an x-axis forward direction and the direction from south to north as a y-axis forward direction, arranging monitoring points in a determinant manner by taking the origin as the center and taking the x-axis direction and the y-axis direction as intervals, monitoring changes of physical quantities before an earthquake by the monitoring points in a fixed time period, and drawing a physical quantity curve by taking time as a horizontal coordinate, wherein the monitored physical quantities comprise: atmospheric pressure, air temperature, ground stress, earth electricity, geomagnetism, altitude, gravity acceleration, water level, water temperature, radon content in water, carbon dioxide content in water, nitrite ion content in water and rainfall;
(2) comparing the shapes of each peak and each trough in the physical quantity curve drawn in the step (1), and recording the peaks and the troughs which continuously appear in each physical quantity curve and are similar in shape;
(3) taking three monitoring points near an area to be detected, comparing the monitoring results of the monitoring points on the same physical quantity in the same time period, marking the connection line of the two monitoring points with the maximum peak or trough area and the minimum area on a map as a straight line m, and then marking the parallel line passing through the other monitoring point as the straight line m as a straight line n;
(4) taking the other three monitoring points near the area to be detected, recording the connecting line of the two monitoring points with the maximum peak or trough area and the minimum area on the map as a straight line e for the monitoring results of different physical quantities in the same time period in the step (3), and then recording the parallel line passing through the other monitoring point to be the straight line e as a straight line f;
(5) the straight lines m, n, e and f are intersected to form a parallelogram, and the position of the parallelogram is the epicenter position;
(6) recording the physical quantities monitored by the monitoring points in the step (3) and the step (4), recording the number of wave crests or wave troughs in a physical quantity curve as N, wherein N +4 is the magnitude of the earthquake in the earthquake area;
(7) the origin time interval: and the unit time is between the unit time before the last peak or trough is finished and six unit times later, or between the nine unit times after the last peak or trough is finished and sixteen unit times later, wherein the unit time is the division value of the abscissa of the corresponding physical quantity curve.
2. A method of predicting an earthquake as recited in claim 1, wherein: the value range of N is more than or equal to 1 and less than or equal to 6, wherein N is a positive integer.
CN201810731412.0A 2018-07-05 2018-07-05 Method for predicting earthquake Active CN109061719B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810731412.0A CN109061719B (en) 2018-07-05 2018-07-05 Method for predicting earthquake

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810731412.0A CN109061719B (en) 2018-07-05 2018-07-05 Method for predicting earthquake

Publications (2)

Publication Number Publication Date
CN109061719A CN109061719A (en) 2018-12-21
CN109061719B true CN109061719B (en) 2021-06-18

Family

ID=64819394

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810731412.0A Active CN109061719B (en) 2018-07-05 2018-07-05 Method for predicting earthquake

Country Status (1)

Country Link
CN (1) CN109061719B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021102675A1 (en) * 2019-11-26 2021-06-03 谢镕键 System and method for performing joint prediction of earthquake on the basis of network

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1166207A (en) * 1995-10-06 1997-11-26 阿莫科公司 Method and appts. for seismic signal processing and exploration
WO2000023821A1 (en) * 1998-10-16 2000-04-27 Strm, Llc Method for 4d permeability analysis of geologic fluid reservoirs
GB2319611B (en) * 1996-11-18 2000-06-21 Western Atlas Int Inc Marine seismic surveying
CN1445561A (en) * 2002-03-18 2003-10-01 王风田 Method for forecasting earthquake
CN105899971A (en) * 2013-11-27 2016-08-24 界标制图有限公司 Method and system for gradiational seismic volume classification
CN106251585A (en) * 2016-06-25 2016-12-21 占行波 Ground water temperature alarm point geometirc graphical calculation prediction earthquake location and the method for scope

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1166207A (en) * 1995-10-06 1997-11-26 阿莫科公司 Method and appts. for seismic signal processing and exploration
GB2319611B (en) * 1996-11-18 2000-06-21 Western Atlas Int Inc Marine seismic surveying
WO2000023821A1 (en) * 1998-10-16 2000-04-27 Strm, Llc Method for 4d permeability analysis of geologic fluid reservoirs
CN1445561A (en) * 2002-03-18 2003-10-01 王风田 Method for forecasting earthquake
CN105899971A (en) * 2013-11-27 2016-08-24 界标制图有限公司 Method and system for gradiational seismic volume classification
CN106251585A (en) * 2016-06-25 2016-12-21 占行波 Ground water temperature alarm point geometirc graphical calculation prediction earthquake location and the method for scope

Also Published As

Publication number Publication date
CN109061719A (en) 2018-12-21

Similar Documents

Publication Publication Date Title
US20190271680A1 (en) Determination risk of natural disaster based on moisture content information
CN111105600B (en) Cutting slope stability dynamic monitoring and early warning system and method based on rainfall condition
Megawati et al. Ground‐motion attenuation relationship for the Sumatran megathrust earthquakes
Press Earthquake prediction
CN104732104A (en) Method for calculating extreme high water levels in different reappearance periods under insufficient long-term tide level data condition
CN113283802A (en) Landslide risk assessment method for complex and difficult mountain area
JP6450129B2 (en) Slope failure prediction method and slope failure prediction device
CN104615847A (en) Tsunami dangerousness forecasting method for the South China Sea based on probabilistic method
Vojtko et al. Late Quaternary fault activity in the Western Carpathians: evidence from the Vikartovce Fault (Slovakia)
CN109242234B (en) Automatic initial assessment method for underground nuclear power plant site based on area discrete method
Lee et al. Development and application of landslide susceptibility analysis techniques using geographic information system (GIS)
Salami et al. Trend analysis of hydro-meteorological variables in the coastal area of Lagos using Mann-Kendall trend and Standard Anomaly Index methods
Wang et al. Historical change and future scenarios of sea level rise in Macau and adjacent waters
CN109061719B (en) Method for predicting earthquake
CN104602297B (en) Radio sensor network monitoring method and system based on compressed sensing
Svalova et al. Complex environmental monitoring in Russia and India
JP4880440B2 (en) Snow accretion prediction method and snow accretion prediction program
Kim et al. Seismoacoustic explosion yield and depth estimation: Insights from the large surface explosion coupling experiment
Lermo et al. Local earthquakes of the Mexico basin in Mexico City: κ, Q, source spectra, and stress drop
CN115687448A (en) Short-rainfall forecasting method and device, electronic equipment and storage medium
Partono et al. Seismic Microzonation of Semarang, Indonesia, based on Probabilistic and Deterministic Combination Analysis
Parise Evaluating rockfall hazard from carbonate slopes in the Sele Valley, Southern Italy
Amakuni et al. Basic analysis on building damages by tsunami due to the 2011 Great East Japan Earthquake Disaster using GIS
Eneva et al. Application of radar Interferometry to Detect subsidence and Uplift at the Heber Geothermal Field, southern california
Dermadi Analysis of Numerical Model Result To Estimate Tsunami Damage Based On Inundation Data

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant