CN109061719B - A method of predicting earthquakes - Google Patents

A method of predicting earthquakes 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
monitoring
trough
earthquake
monitoring points
taking
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.)
Expired - Fee Related
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
Expired - Fee Related 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

一种预测地震的方法,步骤为:在地图上以100km为间隔,以行列式排布监测点,监测点以固定的时间周期监测地震前的参数变化并以时间为横坐标绘制参数变化曲线,对比曲线中每个波峰、波谷,对曲线中连续出现的形状相似波峰、波谷进行记录;将波峰或波谷面积最大和面积最小的两个监测点连线记为m,过另一监测点做m的平行线n,取另外三个监测点,做e、f;直线m、n、e、f相交形成一个平行四边形,该平行四边形所在的位置即为震中位置;曲线中的将波峰或波谷数量记为N,N+4即为震中区的震级。最后一个波峰或波谷结束前一个至后六个单位时间之间,或后九个至后十六个单位时间之间,本方法能够预测震中位置、起震时间以及震级。

Figure 201810731412

A method for predicting earthquakes. The steps are: arranging monitoring points in a determinant pattern at intervals of 100 km on a map, monitoring the parameter changes before the earthquake with a fixed time period, and drawing a parameter change curve with time as the abscissa, Compare each peak and trough in the curve, and record the continuous peaks and troughs with similar shapes in the curve; record the line connecting the two monitoring points with the largest and smallest areas of the peak or trough as m, and pass another monitoring point as m The parallel line n of , take the other three monitoring points as e, f; the straight lines m, n, e, f intersect to form a parallelogram, and the position of the parallelogram is the epicenter position; the number of crests or troughs in the curve Denoted as N, N+4 is the magnitude of the epicenter. The method can predict the epicenter location, onset time, and magnitude between one and six units before the end of the last wave crest or trough, or between nine and sixteen units after the end.

Figure 201810731412

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 A method of predicting earthquakes Expired - Fee Related CN109061719B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810731412.0A CN109061719B (en) 2018-07-05 2018-07-05 A method of predicting earthquakes

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810731412.0A CN109061719B (en) 2018-07-05 2018-07-05 A method of predicting earthquakes

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 Expired - Fee Related CN109061719B (en) 2018-07-05 2018-07-05 A method of predicting earthquakes

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
CN108492236B (en) Multiple current Tsunami disaster appraisal procedure based on Monte Carlo stochastic simulation
CN103090831B (en) The determination methods of icing region electric power line ice-covering thickness
CN104615847A (en) Tsunami dangerousness forecasting method for the South China Sea based on probabilistic method
CN115457739B (en) Geological disaster early warning method and device, electronic equipment and storage medium
CN112085270A (en) Storm surge extreme water level prediction method based on random statistical model and hydrodynamic model coupling
Megawati et al. Ground‐motion attenuation relationship for the Sumatran megathrust earthquakes
Vojtko et al. Late Quaternary fault activity in the Western Carpathians: evidence from the Vikartovce Fault (Slovakia)
CN111563135A (en) GM(1,3) Model Urban Land Subsidence Prediction Method Combining Terrain Factors and Neural Networks
JP4880440B2 (en) Snow accretion prediction method and snow accretion prediction program
Svalova et al. Complex environmental monitoring in Russia and India
CN109061719B (en) A method of predicting earthquakes
CN104602297B (en) Radio sensor network monitoring method and system based on compressed sensing
CN105654488A (en) Method for estimating silt interception amount of silt dam based on GIS
CN119494539A (en) A collapse geological disaster vulnerability model and risk assessment method and system
Lermo et al. Local earthquakes of the Mexico basin in Mexico City: κ, Q, source spectra, and stress drop
CN117709711A (en) Ecological slope protection instability early warning method and system based on TRIGRS model
CN117724158A (en) Method for evaluating ground stress based on seismic activity parameter b value
JP5770051B2 (en) Immediate estimation method of the rupture area of a large earthquake
Amakuni et al. Basic analysis on building damages by tsunami due to the 2011 Great East Japan Earthquake Disaster using GIS
Toprak et al. Seismic damage probabilities for segmented buried pipelines in liquefied soils
Anyadiegwu et al. Assessment of flood hazard in Owerri West, Imo State, Nigeria using Analytical Hierarchical Process and Index Based Approach
Parise Evaluating rockfall hazard from carbonate slopes in the Sele Valley, Southern Italy
JPH09237262A (en) Damage estimating method at the time of earthquake
CN117828384B (en) Geological disaster risk point monitoring and emergency response system and method
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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210618