CN102628955A - Method for acquiring earthquake early warning magnitude - Google Patents

Method for acquiring earthquake early warning magnitude Download PDF

Info

Publication number
CN102628955A
CN102628955A CN2012101202459A CN201210120245A CN102628955A CN 102628955 A CN102628955 A CN 102628955A CN 2012101202459 A CN2012101202459 A CN 2012101202459A CN 201210120245 A CN201210120245 A CN 201210120245A CN 102628955 A CN102628955 A CN 102628955A
Authority
CN
China
Prior art keywords
earthquake
magnitude
warning
station
seismic event
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.)
Granted
Application number
CN2012101202459A
Other languages
Chinese (zh)
Other versions
CN102628955B (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.)
CHENGDU MEIHUAN TECHNOLOGY Co Ltd
INSTITUTE OF CARE-LIFE
Original Assignee
CHENGDU MEIHUAN TECHNOLOGY Co Ltd
INSTITUTE OF CARE-LIFE
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 CHENGDU MEIHUAN TECHNOLOGY Co Ltd, INSTITUTE OF CARE-LIFE filed Critical CHENGDU MEIHUAN TECHNOLOGY Co Ltd
Priority to CN201210120245.9A priority Critical patent/CN102628955B/en
Publication of CN102628955A publication Critical patent/CN102628955A/en
Application granted granted Critical
Publication of CN102628955B publication Critical patent/CN102628955B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention discloses a method for acquiring earthquake early warning magnitude, and the objective of the invention is to provide a novel method combining two earthquake early warning magnitude acquisition means of prediction and actual measurement to overcome shortcomings of existing early warning magnitude calculation methods. The technical scheme of the invention is that the method comprises the following steps: 1 P waves and S waves of earthquake are acquired by using a plurality of earthquake wave acquisition stations which are close to the epicenter; 2 feature information of first n seconds of the earthquake waves acquired by each of the earthquake wave acquisition stations after the earthquake is extracted, and earthquake early warning magnitude of each of the stations is calculated by using an earthquake magnitude prediction formula; 3 earthquake magnitude of each of the stations is calculated by using an actual measurement formula; 4 if the earthquake magnitude of one station calculated by using the earthquake magnitude actual measurement formula is greater than the earthquake early warning magnitude calculated by using the earthquake magnitude prediction formula or the end of the earthquake waves is monitored by the station, then the earthquake magnitude calculated by using the earthquake magnitude actual measurement formula is selected as the earthquake early warning magnitude for the station; otherwise the predicted earthquake magnitude is selected as the early warning magnitude; 5 the average value of the earthquake early warning magnitude of each of the earthquake wave acquisition stations is selected as comprehensive earthquake early warning magnitude of the earthquake early warning.

Description

A kind of earthquake pre-warning earthquake magnitude acquisition methods
Technical field
The present invention relates to the earthquake pre-warning technology, especially a kind of acquisition methods of earthquake pre-warning earthquake magnitude.
Background technology
Earthquake pre-warning (also can be described as the urgent speed newspaper of earthquake, for example in Japan) is meant after earthquake takes place, and gets the jump on seismic wave propagation before fortified defence area, shifts to an earlier date several seconds to tens of seconds to fortified defence area and gives the alarm, to reduce local loss.P Bob S velocity of wave propagation was fast after earthquake took place, and arrived fortified defence area faster than the S ripple.
In the prior art; Usually can in the more incident area of earthquake some seismic events be installed in advance and detect the station; Seismic event monitor station standing-meeting after the earthquake sends to earthquake pre-warning server (perhaps being called early warning center, early warning publisher server) in real time with the P ripple that collects, S ripple, and the earthquake pre-warning server obtains this magnitude of earthquake through extraction P ripple, S ripple information calculations.
The earthquake pre-warning earthquake magnitude computing formula that this area exists at present is a predictor formula.
Wherein predictor formula can but to be not limited to be that
Figure 2012101202459100002DEST_PATH_IMAGE001
or
Figure 215344DEST_PATH_IMAGE002
calculates magnitude of earthquake; Wherein is magnitude of earthquake;
Figure 817970DEST_PATH_IMAGE004
is the peak displacement of former seconds after the earthquake P ripples;
Figure 2012101202459100002DEST_PATH_IMAGE005
is epicentral distance, and is the eigenperiod of former seconds after the earthquake P ripples.
The principle of magnitude of earthquake predictor method is to obtain predictor formula according to historical seismic wave waveform through statistics, and what this predictor formula reacted is the statistical relationship between earthquake P ripple first motion characteristic (in general being preceding 3 seconds waveform character) and the magnitude of earthquake.Because the result who always exists individual difference, predictor formula to calculate is always big than actual measurement formula result calculated deviation.
In addition, the formula that is used to calculate magnitude of earthquake after earthquake finishes is an earthquake magnitude actual measurement formula.The actual measurement formula can but be not limited to body wave magnitude computing formula
Figure 2012101202459100002DEST_PATH_IMAGE007
; Wherein
Figure 429135DEST_PATH_IMAGE003
is magnitude of earthquake;
Figure 243507DEST_PATH_IMAGE008
is the earthquake displacement, and its value is the arithmetic mean of two levels to maximum earthquake displacement;
Figure 190603DEST_PATH_IMAGE005
is epicentral distance;
Figure 2012101202459100002DEST_PATH_IMAGE009
is calibration functions; Calibration functions is when on different epicentral distance observation stations, measuring earthquake magnitude with the particle movement maximum amplitude; Because of seismic event with range attenuation the corrected value that must add, the earthquake magnitude value of sympathetic earthquakes when its numerical value is equivalent on this distance, to record particle movement speed and is 1 μ m/s.The actual measurement formula is after earthquake is calmed down, and is used to calculate the formula of magnitude of earthquake, also is the earthquake magnitude computing formula that earthquake department uses now.
Summary of the invention
The objective of the invention is provides a kind of new earthquake pre-warning earthquake magnitude acquisition methods to the deficiency that existing early warning earthquake magnitude computing method exist; Need to prove; The actual measurement formula is after earthquake is calmed down, just to be applied to earthquake magnitude calculate in the prior art; And in the present invention, this formula is used for before earthquake calms down by novelty.Concrete steps are following:
The earthquake pre-warning earthquake magnitude acquisition methods that the present invention proposes may further comprise the steps:
Step 1: utilize the some seismic events nearer to gather P ripple and the S ripple that the station is gathered earthquake from the earthquake centre;
Step 2: extract the characteristic information of preceding n seismic event second of each seismic event collection station collection after the earthquake, and utilize magnitude of earthquake predictor formula to calculate magnitude of earthquake;
Step 3: utilize magnitude of earthquake actual measurement formula to calculate magnitude of earthquake;
Step 4: judge that each seismic wave gathers the earthquake magnitude whether earthquake magnitude that the magnitude of earthquake actual measurement formula of the station calculates calculates greater than magnitude of earthquake predictor formula:, then magnitude of earthquake is surveyed the earthquake magnitude that formula calculates and gather the earthquake pre-warning earthquake magnitude of the station as this seismic event if the earthquake magnitude that the magnitude of earthquake at certain seismic event picking platform station actual measurement formula calculates is big or this seismic event is gathered station and monitored seismic event and calm down; Otherwise, the earthquake magnitude that magnitude of earthquake predictor formula calculates is gathered the earthquake pre-warning earthquake magnitude of the station as this seismic event;
Step 5: the earthquake pre-warning earthquake magnitude to each seismic event picking platform station is averaged; Said mean value is the comprehensive earthquake pre-warning earthquake magnitude of secondary earthquake early warning for this reason.
Preferably; The method that the earthquake pre-warning earthquake magnitude of the calculating various places seismic wave collection station is averaged in the step 5 is: the variance of directly the earthquake pre-warning earthquake magnitude at each seismic event picking platform station being asked the earthquake pre-warning earthquake magnitude at arithmetic mean value or elder generation each seismic event picking platform station of calculating; Remove earthquake pre-warning earthquake magnitude then, at last remaining earthquake pre-warning earthquake magnitude is averaged greater than 2 times or 3 times variances.
Preferably, gathering station place at seismic event adopts strong-motion instrument or acceleration transducer to gather seismic signal.
Preferably; In the said step 2; Predictor formula is that
Figure 208369DEST_PATH_IMAGE001
or
Figure 469586DEST_PATH_IMAGE002
calculates magnitude of earthquake; Wherein
Figure 699317DEST_PATH_IMAGE003
is magnitude of earthquake;
Figure 754997DEST_PATH_IMAGE004
be after the earthquake before the peak displacement of n P second ripple;
Figure 571644DEST_PATH_IMAGE005
is epicentral distance,
Figure 370972DEST_PATH_IMAGE006
be after the earthquake before the eigenperiod of n P second ripple;
In the said step 3; Said actual measurement formula is
Figure 458139DEST_PATH_IMAGE007
; Wherein
Figure 12617DEST_PATH_IMAGE003
is magnitude of earthquake;
Figure 316560DEST_PATH_IMAGE008
is the earthquake displacement, and its value is seismic shock acceleration that said strong-motion instrument or acceleration transducer collect is carved into current magnitude of earthquake calculation time when earthquake takes place a quadratic integral;
Figure 591683DEST_PATH_IMAGE005
is epicentral distance, and
Figure 858323DEST_PATH_IMAGE009
is calibration functions.
Preferably, utilization is gathered the preceding after the earthquake n of the station collections P ripple and the S ripple of second from earthquake centre 40 kilometer range with interior some seismic events in the said step 1.
Preferably, said n is 3.
The invention allows for another kind of earthquake pre-warning earthquake magnitude acquisition methods, may further comprise the steps: step 1: utilize the some seismic events in 45 kilometer range of earthquake centre to gather P ripple and the S ripple that the station is gathered earthquake;
Step 2: extract seismic event after the earthquake and gather the characteristic information in preceding n seismic event second that the station gathers, and utilize magnitude of earthquake predictor formula to calculate magnitude of earthquake;
Step 3: utilize magnitude of earthquake actual measurement formula to calculate magnitude of earthquake;
Step 4: judge seismic waveshape whether peak value cross or seismic event is calmed down: the seismic event peak value is crossed or seismic event is calmed down if certain seismic event picking platform station detects, and then magnitude of earthquake is surveyed the earthquake magnitude that formula calculates and gathers the earthquake pre-warning earthquake magnitude of the station as this seismic event; Otherwise, the earthquake magnitude that magnitude of earthquake predictor formula calculates is gathered the earthquake pre-warning earthquake magnitude of the station as this seismic event;
Step 5: the earthquake pre-warning earthquake magnitude to each seismic event picking platform station is averaged; Said mean value is the comprehensive earthquake pre-warning earthquake magnitude of secondary earthquake early warning for this reason.
Preferably; The method that the earthquake pre-warning earthquake magnitude of the calculating various places seismic wave collection station is averaged in the step 5 is: the variance of directly the earthquake pre-warning earthquake magnitude at each seismic event picking platform station being asked the earthquake pre-warning earthquake magnitude at arithmetic mean value or elder generation each seismic event picking platform station of calculating; Remove earthquake pre-warning earthquake magnitude then, at last remaining earthquake pre-warning earthquake magnitude is averaged greater than 2 times or 3 times variances.
Preferably, gathering station place at seismic event adopts strong-motion instrument or acceleration transducer to gather seismic signal; Said n is 3.
Preferably; In the said step 2; Predictor formula is that
Figure 255807DEST_PATH_IMAGE001
or
Figure 47045DEST_PATH_IMAGE002
calculates magnitude of earthquake; Wherein
Figure 751958DEST_PATH_IMAGE003
is magnitude of earthquake;
Figure 312252DEST_PATH_IMAGE004
is the peak displacement of preceding after the earthquake 3 seconds P ripples; is epicentral distance, and
Figure 657707DEST_PATH_IMAGE006
is the eigenperiod of preceding after the earthquake 3 seconds P ripples;
In the said step 3; Said actual measurement formula is
Figure 602529DEST_PATH_IMAGE007
; Wherein is magnitude of earthquake;
Figure 959878DEST_PATH_IMAGE008
is the earthquake displacement, and its value is seismic shock acceleration that said strong-motion instrument or acceleration transducer collect is carved into current magnitude of earthquake calculation time when earthquake takes place a quadratic integral;
Figure 460130DEST_PATH_IMAGE005
is epicentral distance, and is calibration functions.
Last the present invention proposes a kind of earthquake pre-warning earthquake magnitude acquisition methods, may further comprise the steps:
Step 1: utilize the some seismic events nearer to gather P ripple and the S ripple that the station is gathered earthquake from the earthquake centre;
Step 2: the characteristic information in preceding n seismic event second of extraction each seismic event collection station collection after the earthquake, and utilize magnitude of earthquake predictor formula to calculate magnitude of earthquake;
Step 3: utilize magnitude of earthquake actual measurement formula to calculate magnitude of earthquake;
Step 4: judge that each seismic event gathers the earthquake magnitude whether earthquake magnitude that the magnitude of earthquake actual measurement formula of the station calculates calculates greater than magnitude of earthquake predictor formula:, then magnitude of earthquake is surveyed the earthquake magnitude that formula calculates and gather the earthquake pre-warning earthquake magnitude of the station as this seismic event if the earthquake magnitude that the magnitude of earthquake at certain seismic event picking platform station actual measurement formula calculates is big or this seismic event is gathered station and monitored seismic event and calm down; Otherwise, judge whether peak value mistake of seismic waveshape that this seismic event gathers the station: if peak value mistake, then magnitude of earthquake is surveyed the earthquake magnitude that formula calculates and gather the earthquake pre-warning earthquake magnitude of the station as this seismic event; Otherwise, the earthquake magnitude that magnitude of earthquake predictor formula calculates is gathered the earthquake pre-warning earthquake magnitude of the station as this seismic event;
Step 5: the earthquake pre-warning earthquake magnitude to each seismic event picking platform station is averaged; Said mean value is the comprehensive earthquake pre-warning earthquake magnitude of secondary earthquake early warning for this reason.
Preferably; The method that the earthquake pre-warning earthquake magnitude of the calculating various places seismic wave collection station is averaged in the step 5 is: the variance of directly the earthquake pre-warning earthquake magnitude at each seismic event picking platform station being asked the earthquake pre-warning earthquake magnitude at arithmetic mean value or elder generation each seismic event picking platform station of calculating; Remove earthquake pre-warning earthquake magnitude then, at last remaining earthquake pre-warning earthquake magnitude is averaged greater than 2 times or 3 times variances.
Preferably, said n is 3, and gathers station place at seismic event and adopt strong-motion instrument to gather seismic signal.
Above-mentioned epicentral distance
Figure 41732DEST_PATH_IMAGE005
is the distance that the liftoff seismic wave in earthquake centre is gathered the station.
In sum, owing to adopted technique scheme, the invention has the beneficial effects as follows:
Earthquake magnitude acquisition methods disclosed by the invention can obtain magnitude of earthquake more accurately estimating utilization measurement method calculating earthquake magnitude on the means.
Adopt strong-motion instrument or acceleration transducer to gather seismic event, effectively avoided causing sseismograph saturated because of magnitude of earthquake is excessive in the prior art, further guaranteed the precision that the earthquake pre-warning earthquake magnitude calculates when violent earthquake takes place.
Traditionally, the measurement method of magnitude of earthquake is peak displacement or a speed calculation of utilizing seismic waveshape after waiting earthquake to finish.And the measurement method that adopts among the present invention proposes in early-warning parameters computation process, and with the peak displacement or the velocity estimation early warning earthquake magnitude that have measured, the early warning earthquake magnitude among the present invention is called as earthquake pre-warning earthquake magnitude measured value.The measurement method that adopts among the present invention need not be waited until after earthquake finishes just can calculate earthquake magnitude, has better ageing with respect to existing earthquake measurement method.Simultaneously, at the initial stage that the early warning earthquake magnitude calculates, that utilizes earthquake P ripple estimates the early warning earthquake magnitude, helps the earthquake pre-warning earthquake magnitude and rapidly converges to end value.
Embodiment
Disclosed all characteristics in this instructions, or the step in disclosed all methods or the process except mutually exclusive characteristic and/or the step, all can make up by any way.
Disclosed arbitrary characteristic in this instructions (comprising any accessory claim, summary) is only if special narration all can be replaced by other equivalences or the alternative features with similar purpose.That is, only if special narration, each characteristic is an example in a series of equivalences or the similar characteristics.
Earthquake pre-warning earthquake magnitude acquisition methods disclosed by the invention is estimated and is surveyed two kinds of methods of obtaining magnitude of earthquake and organically combine of the prior art.
Embodiment 1: a kind of earthquake pre-warning earthquake magnitude acquisition methods may further comprise the steps:
Step 1: utilize from earthquake centre 45 kilometer range with interior several, for example 1 ~ 6, seismic event is gathered P ripple and the S ripple that the station is gathered earthquake;
Step 2: extract the characteristic information of preceding 3 seconds P ripples of each seismic event collection station collection after the earthquake; And utilize magnitude of earthquake predictor formula
Figure 902984DEST_PATH_IMAGE001
to calculate magnitude of earthquake; Wherein
Figure 641264DEST_PATH_IMAGE003
is magnitude of earthquake;
Figure 512178DEST_PATH_IMAGE004
is the peak displacement of preceding after the earthquake 3 seconds P ripples, and
Figure 635992DEST_PATH_IMAGE005
is epicentral distance;
Step 3: utilize magnitude of earthquake actual measurement formula
Figure 982660DEST_PATH_IMAGE007
to calculate magnitude of earthquake; Wherein
Figure 958968DEST_PATH_IMAGE003
is magnitude of earthquake;
Figure 111601DEST_PATH_IMAGE008
is the earthquake displacement, and its value is seismic shock acceleration that said strong-motion instrument or acceleration transducer collect is carved into current magnitude of earthquake calculation time when earthquake takes place a quadratic integral;
Figure 355501DEST_PATH_IMAGE005
is epicentral distance, and
Figure 106026DEST_PATH_IMAGE009
is calibration functions;
Step 4: judge that each seismic event gathers the earthquake magnitude whether earthquake magnitude that the magnitude of earthquake actual measurement formula of the station calculates calculates greater than magnitude of earthquake predictor formula:, then magnitude of earthquake is surveyed the earthquake magnitude that formula calculates and gather the earthquake pre-warning earthquake magnitude of the station as this seismic event if the earthquake magnitude that the magnitude of earthquake at certain seismic event picking platform station actual measurement formula calculates is big or this seismic event is gathered station and monitored seismic event and calm down; Otherwise, the earthquake magnitude that magnitude of earthquake predictor formula calculates is gathered the earthquake pre-warning earthquake magnitude of the station as this seismic event;
Step 5: the earthquake pre-warning earthquake magnitude to each seismic event picking platform station is asked the arithmetic mean value, and said arithmetic mean value is the comprehensive earthquake pre-warning earthquake magnitude of secondary earthquake early warning for this reason.
Embodiment 2:
A kind of earthquake pre-warning earthquake magnitude acquisition methods may further comprise the steps:
Step 1: utilize from earthquake centre 45 kilometer range with interior several, for example, 1 ~ 6, seismic event is gathered P ripple and the S ripple that station strong-motion instrument is gathered earthquake;
Step 2: extract the characteristic information of preceding 4 seconds P ripples of each seismic event collection station collection after the earthquake; And utilize magnitude of earthquake predictor formula
Figure 68165DEST_PATH_IMAGE002
to calculate magnitude of earthquake, is the eigenperiod of preceding after the earthquake 4 seconds P ripples;
Step 3: utilize magnitude of earthquake actual measurement formula
Figure 748994DEST_PATH_IMAGE007
to calculate magnitude of earthquake; Wherein
Figure 171885DEST_PATH_IMAGE003
is magnitude of earthquake;
Figure 119856DEST_PATH_IMAGE008
is the earthquake displacement, and its value is seismic shock acceleration that said strong-motion instrument or acceleration transducer collect is carved into current magnitude of earthquake calculation time when earthquake takes place a quadratic integral;
Figure 614291DEST_PATH_IMAGE005
is epicentral distance, and
Figure 567203DEST_PATH_IMAGE009
is calibration functions;
Step 4: judge that each seismic event gathers the earthquake magnitude whether earthquake magnitude that the magnitude of earthquake actual measurement formula of the station calculates calculates greater than magnitude of earthquake predictor formula:, then magnitude of earthquake is surveyed the earthquake magnitude that formula calculates and gather the earthquake pre-warning earthquake magnitude of the station as this seismic event if the earthquake magnitude that the magnitude of earthquake at certain seismic event picking platform station actual measurement formula calculates is big or this seismic event is gathered station and monitored seismic event and calm down; Otherwise, the earthquake magnitude that magnitude of earthquake predictor formula calculates is gathered the earthquake pre-warning earthquake magnitude of the station as this seismic event;
Step 5: the earthquake pre-warning earthquake magnitude to each seismic event picking platform station is averaged, and said mean value is the comprehensive earthquake pre-warning earthquake magnitude of secondary earthquake early warning for this reason.
Embodiment 3:
A kind of earthquake pre-warning earthquake magnitude acquisition methods may further comprise the steps: utilize from earthquake centre 45 kilometer range and gather P ripple and the S ripple that station strong-motion instrument is gathered earthquake with 1 ~ 6 interior tableland seismic wave;
Step 2: the characteristic information in the preceding 5 seconds P ripples of extraction each seismic event collection station collection after the earthquake; And utilize magnitude of earthquake predictor formula or to calculate magnitude of earthquake; Wherein is magnitude of earthquake;
Figure 655803DEST_PATH_IMAGE004
is the peak displacement of preceding after the earthquake 5 seconds P ripples;
Figure 187541DEST_PATH_IMAGE005
is epicentral distance, and
Figure 345990DEST_PATH_IMAGE006
is the eigenperiod of preceding after the earthquake 5 seconds P ripples;
Step 3: utilize magnitude of earthquake actual measurement formula
Figure 119910DEST_PATH_IMAGE007
to calculate magnitude of earthquake; Wherein
Figure 47415DEST_PATH_IMAGE003
is magnitude of earthquake;
Figure 747125DEST_PATH_IMAGE008
is the earthquake displacement, and its value is seismic shock acceleration that said strong-motion instrument or acceleration transducer collect is carved into current magnitude of earthquake calculation time when earthquake takes place a quadratic integral;
Figure 455186DEST_PATH_IMAGE005
is epicentral distance, and
Figure 767219DEST_PATH_IMAGE009
is calibration functions;
Step 4: judge whether seismic waveshape had occurred peak value (being maximum point) or seismic event is calmed down: the seismic event of the station peak value occurred or seismic event is calmed down if certain seismic event is gathered, and then magnitude of earthquake was surveyed the earthquake pre-warning earthquake magnitude of the earthquake magnitude of formula calculating as this seismic event collection station; Otherwise, the earthquake magnitude that magnitude of earthquake predictor formula calculates is gathered the earthquake pre-warning earthquake magnitude of the station as this seismic event;
Step 5: the earthquake pre-warning earthquake magnitude to each seismic event picking platform station is averaged, and said mean value is the comprehensive earthquake pre-warning earthquake magnitude of secondary earthquake early warning for this reason.
Embodiment 4:
A kind of earthquake pre-warning earthquake magnitude acquisition methods may further comprise the steps: step 1: utilize from earthquake centre 45 kilometer range and gather P ripple and the S ripple that station strong-motion instrument is gathered earthquake with 10 ~ 20 interior tableland seismic waves;
Step 2: the characteristic information in the preceding 3 seconds P ripples of extraction each seismic event collection station collection after the earthquake; And utilize magnitude of earthquake predictor formula
Figure 316274DEST_PATH_IMAGE001
or
Figure 626033DEST_PATH_IMAGE002
to calculate magnitude of earthquake; Wherein
Figure 759074DEST_PATH_IMAGE003
is magnitude of earthquake;
Figure 874798DEST_PATH_IMAGE004
is the peak displacement of preceding after the earthquake 3 seconds P ripples; is epicentral distance, and
Figure 257554DEST_PATH_IMAGE006
is the eigenperiod of preceding after the earthquake 3 seconds P ripples;
Step 3: utilize magnitude of earthquake actual measurement formula
Figure 110848DEST_PATH_IMAGE007
to calculate magnitude of earthquake; Wherein
Figure 30262DEST_PATH_IMAGE003
is magnitude of earthquake;
Figure 724549DEST_PATH_IMAGE008
is the earthquake displacement, and its value is seismic shock acceleration that said strong-motion instrument or acceleration transducer collect is carved into current magnitude of earthquake calculation time when earthquake takes place a quadratic integral;
Figure 438427DEST_PATH_IMAGE005
is epicentral distance, and
Figure 483743DEST_PATH_IMAGE009
is calibration functions;
Step 4: judge that each seismic event gathers the earthquake magnitude whether earthquake magnitude that the magnitude of earthquake actual measurement formula of the station calculates calculates greater than magnitude of earthquake predictor formula:, then magnitude of earthquake is surveyed the earthquake magnitude that formula calculates and gather the earthquake pre-warning earthquake magnitude of the station as this seismic event if the earthquake magnitude that the magnitude of earthquake at certain seismic event picking platform station actual measurement formula calculates is big or this seismic event is gathered station and monitored seismic event and calm down; Otherwise, execution in step 5;
Step 5: judge whether peak value mistake of seismic waveshape that this seismic event gathers the station: if peak value mistake, then magnitude of earthquake is surveyed the earthquake magnitude that formula calculates and gather the earthquake pre-warning earthquake magnitude of the station as this seismic event; Otherwise, the earthquake magnitude that magnitude of earthquake predictor formula calculates is gathered the earthquake pre-warning earthquake magnitude of the station as this seismic event;
Step 6: consider that the make mistakes seismic event substantial deviation that can cause its collection of seismic event picking platform station equipment sometimes is actual; We can calculate the variance of the earthquake pre-warning earthquake magnitude at each seismic event picking platform station earlier for this reason; Remove earthquake pre-warning earthquake magnitude then, at last remaining earthquake pre-warning earthquake magnitude is averaged greater than 2 times or 3 times variances;
Said average is the comprehensive earthquake pre-warning earthquake magnitude of secondary earthquake early warning for this reason.
If need, can the above earthquake pre-warning earthquake magnitude that calculates be transformed to the earthquake magnitude of other form according to the seismology common method, for example, surface wave magnitude etc.
What be worth explanation is: to gather the station all be to adopt sseismograph to gather seismic waveshape to seismic event in the prior art, and what sseismograph detected is speed, and when running into the big earthquake of earthquake magnitude, sseismograph is saturated easily, thereby can not obtain seismic signal accurately.When Wenchuan earthquake, the sseismograph in area, most of Sichuan differs greatly with actual earthquake magnitude because the magnitude of earthquake that the sseismograph in the saturated and feasible application Sichuan Province measures is 7.3 grades.The present invention adopts the strong-motion instrument or the acceleration transducer of sense acceleration, can well avoid " saturated " phenomenon, has further guaranteed earthquake magnitude reliability and the precision obtained.
The present invention is not limited to aforesaid embodiment.The present invention expands to any new feature or any new combination that discloses in this manual, and the arbitrary new method that discloses or step or any new combination of process.

Claims (10)

1. an earthquake pre-warning earthquake magnitude acquisition methods is characterized in that, may further comprise the steps:
Step 1: utilize the some seismic events nearer to gather P ripple and the S ripple that the station is gathered earthquake from the earthquake centre;
Step 2: extract the characteristic information of preceding n seismic event second of each seismic event collection station collection after the earthquake, and utilize magnitude of earthquake predictor formula to calculate magnitude of earthquake;
Step 3: utilize magnitude of earthquake actual measurement formula to calculate magnitude of earthquake;
Step 4: the earthquake magnitude of surveying formula calculating as if the magnitude of earthquake at certain seismic event picking platform station monitors seismic event than big or this seismic event collection station and calms down, then magnitude of earthquake is surveyed the earthquake pre-warning earthquake magnitude of the earthquake magnitude of formula calculating as this seismic event collection station; Otherwise, the earthquake magnitude that magnitude of earthquake predictor formula calculates is gathered the earthquake pre-warning earthquake magnitude of the station as this seismic event;
Step 5: the earthquake pre-warning earthquake magnitude to each seismic event picking platform station is averaged; Said mean value is the comprehensive earthquake pre-warning earthquake magnitude of secondary earthquake early warning for this reason.
2. a kind of earthquake pre-warning earthquake magnitude acquisition methods according to claim 1 is characterized in that, gathers station place at seismic event and adopts strong-motion instrument or acceleration transducer to gather seismic signal.
3. a kind of earthquake pre-warning earthquake magnitude acquisition methods according to claim 2; It is characterized in that; In the said step 2; Predictor formula is that
Figure 2012101202459100001DEST_PATH_IMAGE002
or
Figure 2012101202459100001DEST_PATH_IMAGE004
calculates magnitude of earthquake; Wherein
Figure 2012101202459100001DEST_PATH_IMAGE006
is magnitude of earthquake;
Figure 2012101202459100001DEST_PATH_IMAGE008
be after the earthquake before the peak displacement of n P second ripple; is epicentral distance,
Figure 2012101202459100001DEST_PATH_IMAGE012
be after the earthquake before the eigenperiod of n P second ripple;
In the said step 3; Said actual measurement formula is
Figure 2012101202459100001DEST_PATH_IMAGE014
; Wherein
Figure 742504DEST_PATH_IMAGE006
is magnitude of earthquake;
Figure 2012101202459100001DEST_PATH_IMAGE016
is the earthquake displacement, and its value is seismic shock acceleration that said strong-motion instrument or acceleration transducer collect is carved into current magnitude of earthquake calculation time when earthquake takes place a quadratic integral;
Figure 472694DEST_PATH_IMAGE010
is epicentral distance, and
Figure 2012101202459100001DEST_PATH_IMAGE018
is calibration functions.
4. a kind of earthquake pre-warning earthquake magnitude acquisition methods according to claim 3 is characterized in that, utilizes in the said step 1 from earthquake centre 45 kilometer range and gathers preceding after the earthquake 3 seconds P ripple and the S ripple of station collection with interior some seismic events.
5. a kind of earthquake pre-warning earthquake magnitude acquisition methods according to claim 4; It is characterized in that the method for in the said step 5 the earthquake pre-warning earthquake magnitude at each seismic event picking platform station being averaged is directly the earthquake pre-warning earthquake magnitude at each seismic event picking platform station to be asked the arithmetic mean value; Perhaps calculate the variance of the earthquake pre-warning earthquake magnitude at each seismic event picking platform station earlier, remove earthquake pre-warning earthquake magnitude then, at last remaining earthquake pre-warning earthquake magnitude is averaged greater than 2 times or 3 times variances.
6. an earthquake pre-warning earthquake magnitude acquisition methods is characterized in that, may further comprise the steps:
Step 1: utilize the some seismic events in 45 kilometer range of earthquake centre to gather P ripple and the S ripple that the station is gathered earthquake;
Step 2: the characteristic information in preceding n seismic event second of extraction each seismic event picking platform station collection after the earthquake, and utilize magnitude of earthquake predictor formula to calculate magnitude of earthquake;
Step 3: utilize magnitude of earthquake actual measurement formula to calculate magnitude of earthquake;
Step 4: judge seismic waveshape whether peak value cross or seismic event is calmed down: the seismic event peak value is crossed or seismic event is calmed down if certain seismic event picking platform station detects, and then magnitude of earthquake is surveyed the earthquake magnitude that formula calculates and gathers the earthquake pre-warning earthquake magnitude of the station as this seismic event; Otherwise, the earthquake magnitude that magnitude of earthquake predictor formula calculates is gathered the earthquake pre-warning earthquake magnitude of the station as this seismic event;
Step 5: the earthquake pre-warning earthquake magnitude to each seismic event picking platform station is averaged; Said mean value is the comprehensive earthquake pre-warning earthquake magnitude of secondary earthquake early warning for this reason.
7. a kind of earthquake pre-warning earthquake magnitude acquisition methods according to claim 6 is characterized in that, gathers station place at seismic event and adopts strong-motion instrument or acceleration transducer to gather seismic signal; Said n is 3.
8. a kind of earthquake pre-warning earthquake magnitude acquisition methods according to claim 7; It is characterized in that; In the said step 2; Predictor formula is that
Figure 314748DEST_PATH_IMAGE002
or
Figure 21542DEST_PATH_IMAGE004
calculates magnitude of earthquake; Wherein is magnitude of earthquake;
Figure 720693DEST_PATH_IMAGE008
is the peak displacement of preceding after the earthquake 3 seconds P ripples;
Figure 987727DEST_PATH_IMAGE010
is epicentral distance, and
Figure 999676DEST_PATH_IMAGE012
is the eigenperiod of preceding after the earthquake 3 seconds P ripples;
In the said step 3; Said actual measurement formula is
Figure 730872DEST_PATH_IMAGE014
; Wherein
Figure 724236DEST_PATH_IMAGE006
is magnitude of earthquake;
Figure 790150DEST_PATH_IMAGE016
is the earthquake displacement, and its value is seismic shock acceleration that said strong-motion instrument or acceleration transducer collect is carved into current magnitude of earthquake calculation time when earthquake takes place a quadratic integral;
Figure 589478DEST_PATH_IMAGE010
is epicentral distance, and
Figure 112863DEST_PATH_IMAGE018
is calibration functions.
9. an earthquake pre-warning earthquake magnitude acquisition methods is characterized in that, may further comprise the steps:
Step 1: utilize the some seismic events nearer to gather P ripple and the S ripple that the station is gathered earthquake from the earthquake centre;
Step 2: the characteristic information in preceding n seismic event second of extraction each seismic event collection station collection after the earthquake, and utilize magnitude of earthquake predictor formula to calculate magnitude of earthquake;
Step 3: utilize magnitude of earthquake actual measurement formula to calculate magnitude of earthquake;
Step 4: if the seismic waveshape that the earthquake magnitude that the magnitude of earthquake at certain seismic event picking platform station actual measurement formula calculates is big or this seismic event gathers that station monitors that seismic event has been calmed down or this seismic event is gathered station peak value mistake whether, then magnitude of earthquake is surveyed the earthquake magnitude that formula calculates and gather the earthquake pre-warning earthquake magnitude of the station as this seismic event; Otherwise, the earthquake magnitude that magnitude of earthquake predictor formula calculates is gathered the earthquake pre-warning earthquake magnitude of the station as this seismic event;
Step 5: the earthquake pre-warning earthquake magnitude to each seismic event picking platform station is averaged; Said mean value is the comprehensive earthquake pre-warning earthquake magnitude of secondary earthquake early warning for this reason.
10. a kind of earthquake pre-warning earthquake magnitude acquisition methods according to claim 9 is characterized in that said n is 3, and gathers station place at seismic event and adopt strong-motion instrument to gather seismic signal.
CN201210120245.9A 2012-04-24 2012-04-24 Method for acquiring earthquake early warning magnitude Active CN102628955B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210120245.9A CN102628955B (en) 2012-04-24 2012-04-24 Method for acquiring earthquake early warning magnitude

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210120245.9A CN102628955B (en) 2012-04-24 2012-04-24 Method for acquiring earthquake early warning magnitude

Publications (2)

Publication Number Publication Date
CN102628955A true CN102628955A (en) 2012-08-08
CN102628955B CN102628955B (en) 2014-12-17

Family

ID=46587247

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210120245.9A Active CN102628955B (en) 2012-04-24 2012-04-24 Method for acquiring earthquake early warning magnitude

Country Status (1)

Country Link
CN (1) CN102628955B (en)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103680069A (en) * 2012-09-12 2014-03-26 李宗洪 Earthquake prediction machine, earthquake prediction system and earthquake-precursor monitoring system
CN105676286A (en) * 2016-01-25 2016-06-15 西南交通大学 Real-time earthquake magnitude estimation method for earthquake early warning system
CN106094021A (en) * 2016-06-01 2016-11-09 北京科技大学 A kind of microseism focus method for rapidly positioning based on arrival time difference data base
CN106291662A (en) * 2016-09-09 2017-01-04 吉林大学 Fracturing causes early warning and the Forecasting Methodology that microseism occurs in hot dry rock
CN107741601A (en) * 2017-11-21 2018-02-27 高梧 Novel intelligent earthquake monitoring device and its control method
CN107798826A (en) * 2017-09-28 2018-03-13 北京交通大学 A kind of Forecasting scale of seismic magnitude method for earthquake pre-warning
CN107949796A (en) * 2016-03-14 2018-04-20 欧姆龙株式会社 Seismological sensor and earthquake detection method
CN108089225A (en) * 2017-12-02 2018-05-29 桂林理工大学 A kind of earthquake magnitude Method of fast estimating based on separate unit station first arrival P ripples
CN109901220A (en) * 2019-03-28 2019-06-18 夏仕安 A kind of method of discrimination of the magnitude for earthquake early-warning system
CN110412646A (en) * 2018-04-27 2019-11-05 曾玉婷 A kind of seismic precursor parameter acquiring method and prediction about imminent earthquakes method
CN110824560A (en) * 2019-11-11 2020-02-21 宜宾凯翼汽车有限公司 Earthquake monitoring system and method based on Internet of vehicles
CN114355438A (en) * 2022-01-11 2022-04-15 西南交通大学 Method, device and equipment for estimating seismic magnitude and readable storage medium
CN116594057A (en) * 2023-05-17 2023-08-15 中国地震局地球物理研究所 Earthquake early warning method and device based on deep learning and edge calculation

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
《华南地震》 20100930 张红才等 福建省"九五"地震监测台网台站震级校正及单台震级应用研究 , 第03期 *
YIH-MIN WU ET AL: "Magnitude estimation using the first three seconds P-wave amplitude in earthquake early warning", 《GEOPHYSICAL RESEARCH LETTERS》 *
张红才等: "福建省"九五"地震监测台网台站震级校正及单台震级应用研究", 《华南地震》 *
马强: "地震预警技术研究及应用", 《中国地震局工程力学研究所2008年度工学博士学位论文》 *

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103680069A (en) * 2012-09-12 2014-03-26 李宗洪 Earthquake prediction machine, earthquake prediction system and earthquake-precursor monitoring system
CN105676286B (en) * 2016-01-25 2017-11-14 西南交通大学 A kind of Real-time Seismic Estimation of Seismetic Magnitude method for earthquake early-warning system
CN105676286A (en) * 2016-01-25 2016-06-15 西南交通大学 Real-time earthquake magnitude estimation method for earthquake early warning system
CN107949796A (en) * 2016-03-14 2018-04-20 欧姆龙株式会社 Seismological sensor and earthquake detection method
CN106094021B (en) * 2016-06-01 2017-11-17 北京科技大学 A kind of microseism focus method for rapidly positioning based on arrival time difference database
CN106094021A (en) * 2016-06-01 2016-11-09 北京科技大学 A kind of microseism focus method for rapidly positioning based on arrival time difference data base
CN106291662A (en) * 2016-09-09 2017-01-04 吉林大学 Fracturing causes early warning and the Forecasting Methodology that microseism occurs in hot dry rock
CN107798826A (en) * 2017-09-28 2018-03-13 北京交通大学 A kind of Forecasting scale of seismic magnitude method for earthquake pre-warning
CN107741601A (en) * 2017-11-21 2018-02-27 高梧 Novel intelligent earthquake monitoring device and its control method
CN108089225A (en) * 2017-12-02 2018-05-29 桂林理工大学 A kind of earthquake magnitude Method of fast estimating based on separate unit station first arrival P ripples
CN108089225B (en) * 2017-12-02 2019-08-13 桂林理工大学 A kind of earthquake magnitude Method of fast estimating based on separate unit station first arrival P wave
CN110412646B (en) * 2018-04-27 2024-02-23 曾玉婷 Earthquake precursor parameter acquisition method and earthquake prediction method
CN110412646A (en) * 2018-04-27 2019-11-05 曾玉婷 A kind of seismic precursor parameter acquiring method and prediction about imminent earthquakes method
CN109901220A (en) * 2019-03-28 2019-06-18 夏仕安 A kind of method of discrimination of the magnitude for earthquake early-warning system
CN110824560A (en) * 2019-11-11 2020-02-21 宜宾凯翼汽车有限公司 Earthquake monitoring system and method based on Internet of vehicles
CN114355438A (en) * 2022-01-11 2022-04-15 西南交通大学 Method, device and equipment for estimating seismic magnitude and readable storage medium
CN116594057A (en) * 2023-05-17 2023-08-15 中国地震局地球物理研究所 Earthquake early warning method and device based on deep learning and edge calculation
CN116594057B (en) * 2023-05-17 2024-01-12 中国地震局地球物理研究所 Earthquake early warning method and device based on deep learning and edge calculation

Also Published As

Publication number Publication date
CN102628955B (en) 2014-12-17

Similar Documents

Publication Publication Date Title
CN102628955B (en) Method for acquiring earthquake early warning magnitude
CN102841835B (en) The method and system of hardware performance evaluation and test
Koo et al. Damage detection of shear buildings using deflections obtained by modal flexibility
TW201233986A (en) Method for analyzing structure safety
CN102227189A (en) Method and apparatus for fall detection
Cheynet et al. Damping estimation of large wind-sensitive structures
KR102209749B1 (en) Method and system for generating earthquake acceleration time history
CN103994817A (en) Vibration source identification method based on long-distance optical fiber frequent occurring events
CN102147232A (en) Equipment and method for monitoring multipoint displacement of landslide
JP5591759B2 (en) Epicenter distance estimation method for single station processing
CN108089225A (en) A kind of earthquake magnitude Method of fast estimating based on separate unit station first arrival P ripples
JP6644970B2 (en) Observation system
JP2014169960A (en) Method of predicting arrival time of principal shock of earthquake
CN109489931A (en) A kind of abnormal impact real-time detection method
CN107102342A (en) Gps coordinate time series discontinuity based on common-mode error supplies method
CN104865602A (en) Method of determining epicenter location and earthquake start time through multiple data sources
Park et al. Assessment of vibration serviceability for steel cable-stayed bridges using GNSS data
JP6099123B2 (en) Earthquake identification device, earthquake identification system and earthquake identification method using the same
CN105404708A (en) Actually measured acceleration-based outside period motivation and evaluation method for long-span bridge structures
CN103994816A (en) Identification method based on optical fiber multiple events
JP2013007728A (en) Real time estimation method of hypocentral region of giant earthquake
JP5946074B1 (en) Absolute speed response computing device, absolute speed response computing system using the same, and absolute speed response computing method
CN204666188U (en) A kind ofly measure the land movement of side slope body and the device of water cut
KR101386255B1 (en) Apparatus for measuring earthquake magnitude and distance and method for measuring the same
CN103294851B (en) Based on the bridge structure sensor arrangement method of disaster-ridden harmful sensitivity coefficient index

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant