CN103827694B - Earthquake intensity estimation unit, earthquake intensity method of estimation and computer readable recording medium storing program for performing - Google Patents
Earthquake intensity estimation unit, earthquake intensity method of estimation and computer readable recording medium storing program for performing Download PDFInfo
- Publication number
- CN103827694B CN103827694B CN201280046946.0A CN201280046946A CN103827694B CN 103827694 B CN103827694 B CN 103827694B CN 201280046946 A CN201280046946 A CN 201280046946A CN 103827694 B CN103827694 B CN 103827694B
- Authority
- CN
- China
- Prior art keywords
- unit
- adjacent area
- earthquake intensity
- earthquake
- decay
- 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
Links
- 238000000034 method Methods 0.000 title claims description 36
- 230000001143 conditioned effect Effects 0.000 claims description 11
- 230000007423 decrease Effects 0.000 claims description 5
- 238000003325 tomography Methods 0.000 claims description 3
- 238000009434 installation Methods 0.000 claims 1
- 230000000295 complement effect Effects 0.000 description 27
- 238000004891 communication Methods 0.000 description 17
- 238000003860 storage Methods 0.000 description 14
- 206010044565 Tremor Diseases 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000033001 locomotion Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 2
- 238000000205 computational method Methods 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 239000012925 reference material Substances 0.000 description 1
- 230000000284 resting effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- CCEKAJIANROZEO-UHFFFAOYSA-N sulfluramid Chemical group CCNS(=O)(=O)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)F CCEKAJIANROZEO-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/01—Measuring 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 kind of earthquake intensity estimation unit (10), including: focus determines unit (11), its for be installed in export for specifying the data of amplitude of seismic wave along the seismic detector (20) in each region on the border of particular continent plate time, determine whether the position of focus is in setting district;Adjacent area designating unit (12), if the position of focus is in setting district, then its specify with first it is detected that the adjacent adjacent area in the region that is installed in of the seismic detector of seismic wave;Decay determines unit (13), and it obtains data and determining from the seismic detector that is arranged on adjacent area and the most decays in the amplitude of described data;Fault Rupture site setting unit (14), if decay not yet occurs, then its position based on adjacent area sets the place wherein estimating Fault Rupture in described setting district;And estimate earthquake intensity computing unit (15), it calculates estimation earthquake intensity centered by the place estimating Fault Rupture.
Description
Technical field
The present invention relates to the ground of a kind of earthquake intensity in each place of stage forecast the most in early days
Shake earthquake intensity estimation unit and earthquake intensity method of estimation, and a kind of storage has for realizing described device and the journey of described method
The computer readable recording medium storing program for performing of sequence.
Background technology
In recent years, in Japan, seismism observation system has been introduced into promptly to notify earthquake information to society.?
When earthquake occurs, seismism observation system analysis is immediately captured by the seismic detector being positioned near focus after occurring
Observation data, estimate focus and earthquake magnitude, and estimate the earthquake in each place further based on estimated focus and earthquake magnitude
The time of advent (such as, seeing non-patent file 1) of earthquake intensity and main motion.Now, it is however generally that, to be used as the observation analyzing object
Data are the P ripple initial motions in stage the most in early days.
Seismism observation system then via various media widely by estimated information (focus, earthquake magnitude, earthquake intensity
And the time of advent) issue as earthquake alarm in early days.As a result, people can promptly take action in factory, office
Room and residence are evacuated, and therefore reduce by the loss caused by earthquake.
Reference listing
Non-patent file
Non-patent file 1: " Overview of earthquake early warning and technical
Reference material relating to processing techniques ", seismology and the volcano department of the Chinese Academy of Sciences, meteorological
The Room, on July 29th, 2008, the 3-15 page
Summary of the invention
Problem to be solved by this invention
Incidentally, the seismism observation system disclosed in non-patent file 1 above-mentioned takes focus for ground
Put and estimate earthquake intensity centered by this place.For this reason, in non-patent file 1 above-mentioned
Disclosed seismism observation system for cause near its focus tomography and coverage relatively small earthquake for
It is effective, but relative to making the subduction zone earthquake that tomography succeedingly occurs on a large scale then have the following problems.
First, subduction zone earthquake may on a large scale on cause the Fault Rupture started from focus.In this case,
Seismic wave can propagated than on the broader scope of concentric scope of focus.For this reason, with in non-patent file 1
The scope of the earthquake intensity that disclosed seismism observation system is predicted is compared, and macroseism is quivered on broader scope
Occur.Therefore, not yet utilize earthquake alarm prediction in early days will get an unexpected blow to the resident in the region of earthquake intensity, and damage
Lose and may increase.
The exemplary purpose of the present invention is to solve foregoing problems, and provides a kind of earthquake intensity estimation unit, earthquake strong
Degree method of estimation and computer readable recording medium storing program for performing, even if it is in the case of occurring in subduction zone earthquake, it is also possible to can
Specify the region being predicted trembling by ground, and make it possible to the earthquake intensity estimated in this region.
Means for solving the above
In order to realize the purpose being set forth above, earthquake intensity estimation unit according to an aspect of the present invention is followingly
Shake earthquake intensity estimation unit, comprising:
Focus determines unit, when in each region in the multiple regions set along the border of particular continent plate
The seismic detector installed detects by the seismic wave caused by the generation of earthquake and exports the amplitude for specifying described seismic wave
Data time, described focus determines that unit determines whether the position of the focus of occurred earthquake is in preset areas;
By described focus, adjacent area designating unit, if determining that unit determines that the position of described focus is in described district
In, the most described adjacent area designating unit specify with first it is detected that the district that is installed in of the described seismic detector of described seismic wave
The adjacent area that territory is adjacent;
Decay determine unit, described decay determine unit from the described seismic detector that is arranged on described adjacent area obtain by
This output data, and based on pre-conditioned determine the most send out in by the amplitude specified by accessed data
Raw decay;
Fault Rupture site setting unit, decays if determining unit to determine not yet by described decay, then described
Fault Rupture site setting unit position based on described adjacent area sets in described district wherein to be estimated Fault Rupture
Place;And
Estimating earthquake intensity computing unit, described estimation earthquake intensity computing unit is to set by described Fault Rupture place
To estimating that earthquake intensity calculates centered by place set by cell.
Additionally, in order to realize the purpose being set forth above, earthquake intensity method of estimation according to an aspect of the present invention is
Following earthquake intensity method of estimation, it includes following steps:
A () is when the ground installed in each region in the multiple regions set by the border along particular continent plate
Shake instrument detect by the seismic wave caused by the generation of earthquake and export for the data of the amplitude specifying described seismic wave time,
Whether the position determining the focus of occurred earthquake is in preset areas;
If b () determines that the position of described focus is in described district in step (a)., then specify and first it is detected that institute
State the adjacent area that region that the described seismic detector of seismic wave is installed in is adjacent;
C () obtains, from the described seismic detector being arranged on described adjacent area, the data thus exported, and based on presetting
Condition determines and the most decays in by the described amplitude specified by accessed data;
Not yet decay if d () determines in step (c), then position based on described adjacent area, in described district
Interior setting wherein estimates to occur the place of Fault Rupture;And
E estimation earthquake intensity, centered by the place being set in step (d), is calculated by ().
And, in order to realize the target being set forth above, computer readable recording medium storing program for performing according to an aspect of the present invention
Being following computer readable recording medium storing program for performing, it has program stored therein, and described program includes for making computer perform following step
Rapid order:
A () is when the ground installed in each region in the multiple regions set by the border along particular continent plate
Shake instrument detect by the seismic wave caused by the generation of earthquake and export for the data of the amplitude specifying described seismic wave time,
Whether the position determining the focus of occurred earthquake is in preset areas;
If b () determines that the position of described focus is in described district in step (a)., then specify and first it is detected that institute
State the adjacent area that region that the described seismic detector of seismic wave is installed in is adjacent;
C () obtains, from the described seismic detector being arranged on described adjacent area, the data thus exported, and based on presetting
Condition determines and the most decays in by the described amplitude specified by accessed data;
Not yet decay if d () determines in step (c), then position based on described adjacent area, in described district
Interior setting wherein estimates to occur the place of Fault Rupture;And
E estimation earthquake intensity, centered by the place being set in step (d), is calculated by ().
Beneficial effects of the present invention
As described above, according to the present invention, even if in the case of subduction zone earthquake occurs, it is also possible to reliably refer to
The fixed region being predicted trembling, and the earthquake intensity in this region can be estimated.
Accompanying drawing explanation
Fig. 1 is the block diagram illustrating the configuration of earthquake intensity estimation unit according to an embodiment of the invention.
Fig. 2 is the figure of the example illustrating the mounted state of the seismic detector that uses the most in an embodiment of the present invention.
Fig. 3 is the figure illustrating Japan with continental plate about.
Fig. 4 is the flow chart illustrating the operation of earthquake intensity estimation unit according to an embodiment of the invention.
Fig. 5 A is the figure being shown in and determining the exemplarily seismic wave in the case of not yet decaying in step A3, and
And Fig. 5 B is the figure being shown in and determining the exemplarily seismic wave in the case of decaying in step A3.
Fig. 6 is the frame illustrating and implementing the illustrative computer of earthquake intensity estimation unit 10 according to an embodiment of the invention
Figure.
Detailed description of the invention
Embodiment
Hereinafter, referring to figs. 1 to 6, earthquake intensity estimation unit, earthquake according to an embodiment of the invention will be described
Earthquake intensity method of estimation and program.
Device configures
First, the configuration that Fig. 1 to 3 will be used to describe the earthquake intensity estimation unit 10 according to the present embodiment.Fig. 1 is to show
Go out the block diagram of the configuration of earthquake intensity estimation unit according to an embodiment of the invention.Fig. 2 is to illustrate wherein in the reality of the present invention
Execute the figure of the example of the mounted state of the seismic detector used in example.Fig. 3 is the figure illustrating Japan with continental plate about.
As shown in fig. 1, earthquake intensity estimation unit 10 is connected to multiple shake instrument 20 via network 50.In this enforcement
In example, earthquake intensity estimation unit 10 is also connected to seismism observation system 40 via network 50.
As described in the background, when earthquake occurs, seismism observation system 40 specifies focus also
And calculate earthquake magnitude based on specified focus.It is each to estimate that seismism observation system 40 is additionally based upon calculated earthquake magnitude
Earthquake intensity in side and the time of advent of main motion.
Japan Meteorological Agency earthquake magnitude that then seismism observation system 40 by estimated focus, is calculated, estimated
The time of advent of earthquake intensity and main motion is published to various medium as earthquake alarm in early days.Note, " the shake at this moment calculated
Level " it is Japan Meteorological Agency earthquake magnitude.Japan Meteorological Agency earthquake magnitude is the earthquake magnitude announced by Japan Meteorological Agency, and its computational methods quilt
It is disclosed in non-patent literature 1.
The region that seismic detector 20 is installed in by them is divided in groups.Reference 21a to 21d represents seismic detector 20
Each group.When an earthquake occurs, each seismic detector 20 detects seismic wave, and exports the number being used for specifying the amplitude of this seismic wave
According to (hereinafter referred to as " seismic data ").In the present embodiment, seismic data can be only to specify shaking of seismic wave
The data of width, or may refer to determine the data of the amplitude of seismic wave and waveform.
In the figure 2 example, during seismic detector 20 is installed in four regions i.e. region A to D.Border along continental plate
22 carry out setting regions A to D.In each in region A, B, C and D, the seismic detector 20 of respective sets is installed.Note, day
This and continental plate about are as shown in Fig. 3.
Additionally, as shown in fig. 1, it is single that earthquake intensity estimation unit 10 includes that focus determines that unit 11, adjacent area are specified
Unit 12, decay determine unit 13, Fault Rupture site setting unit 14 and estimate earthquake intensity computing unit 15.
If installing seismic detector 20 in each area to detect by the seismic wave caused by the generation of earthquake and defeated
Going out seismic data, the focus among constituent components the most enumerated above determines that unit 11 determines the focus of the earthquake occurred
Whether the position of (see figure 2) is in preset areas (hereinafter referred to as " setting district ";See Fig. 2) in 30.
If determining that unit 11 determines that the position of focus 23 is in setting district 30 by focus, then adjacent area is specified single
Unit 12 specifies the region adjacent with the region that the seismic detector being initially detected seismic wave is installed in (hereinafter referred to as
" adjacent area ").
Decay determines that the seismic detector institute that unit 13 obtains by being arranged on specified adjacent area from this seismic detector is defeated
The seismic data gone out, and determine by the amplitude specified by accessed seismic data based on pre-conditioned
The most decay.
If additionally, determine that unit 13 determines that decay not yet occurs by decay, then Fault Rupture site setting unit 14
Position based on specified adjacent area, to set in setting district 30 and wherein estimates that the place of Fault Rupture generation is (under
Literary composition is referred to as " Fault Rupture place ") 24.
Estimate that earthquake intensity computing unit 15 calculates estimation earthquake centered by set Fault Rupture place 24 strong
Degree.In the present embodiment, the estimation earthquake intensity output in each place that earthquake intensity computing unit 15 will be calculated is estimated
To seismism observation system 40.In this case, seismism observation system 40, via various media, is the most only broadly sent out
The estimation earthquake intensity that cloth seismism observation system 40 is calculated centered by focus 23, but also issue and pass through earthquake intensity
The estimation earthquake intensity that estimation unit 10 is calculated.
Therefore, in the present embodiment, whether earthquake intensity estimation unit 10 seismic wave definitely is in the region with close focus
Adjacent adjacent area is attenuated.In this case, if the Fault Rupture started from focus is on a large scale
Occur, then decay in adjacent area.Therefore, earthquake intensity estimation unit 10 can determine that Fault Rupture is the most
Occur, and depend on determining that result calculates estimation earthquake intensity on a large scale.As described above, use earthquake strong
, even if in the case of there is subduction zone earthquake wherein, it is also possible to reliably specify trembling in degree estimation unit 10
The region of row prediction, and the earthquake intensity in this region can be estimated.
The configuration of earthquake intensity estimation unit 10 according to the present embodiment now be will be described in further detail.As shown in fig. 1, ground
Shake earthquake intensity estimation unit 10 farther includes communication unit 16.Communication unit 16 via network 50 to seismic detector 20 and seismism
Observation system 40 sends information and receives information from seismic detector 20 and seismism observation system 40.
Received seismic data, when receiving seismic data from seismic detector 20, is exported by communication unit 16
Adjacent area designating unit 12 and decay determine unit 13.Additionally, communication unit 16 will be by estimating earthquake intensity computing unit 15
The estimation earthquake intensity calculated is sent to seismism observation system 40.
In the present embodiment, when seismic detector 20 exports seismic data, seismism observation system 40 specifies focus 23
(see figure 2), and send the positional information (such as, using it to specify longitude and the information of latitude of focus) about focus 23.
Therefore, communication unit 16 receives the positional information (hereinafter referred to as " source position information ") about focus 23, and will
Its output determines unit 11 to focus.In the present embodiment, focus determines that unit 11 has been output at source position information
Perform under conditions of it to determine process.Noting, source position information can be sent as earthquake alarm in early days.
In the present embodiment, as shown in Figures 2 and 3, it is preferred that set with the band-shaped of border 22 along continental plate
Fixed determined the unit 11 setting district 30 used in determining by focus, and especially, so that including border 22.This be because of
For, it is however generally that, Fault Rupture is proximate to the border 22 of continental plate and occurs.
Fig. 3 shows the possible setting district 30a to 30f in Japan.Noting, the setting district 30a shown in Fig. 3 is extremely
30b is only example, and in practice, be based on the continental plate existed near subject area, relevant pass by earthquake focus
Statistical data, the result of simulation of generation etc. of earthquake set described setting district.Although shown in example in figure 3
Region near Japan, but the present embodiment is equally applicable to the country in addition to Japan.
In the present embodiment, adjacent area designating unit 12 can also specify new phase based on specified adjacent area
Neighbouring region.In this case, decay determines that information 13 is installed in the earthquake new adjacent area from the acquisition of this seismic detector
The seismic data of instrument output.Decay determines that unit 13 is then based on pre-conditioned determining whether by accessed earthquake
Amplitude specified by wave datum is decayed.
Further, if it is determined that not yet decay in new adjacent area, then Fault Rupture site setting unit
14 positions based on new adjacent area again, set on the border 22 of continental plate and wherein estimate Fault Rupture
Place (Fault Rupture place).In this case, estimate that earthquake intensity computing unit 15 is with by Fault Rupture site setting list
Further to estimating that earthquake intensity calculates centered by the Fault Rupture place 24 that unit 14 resets.
According to pattern described above, earthquake intensity estimation unit 10 can be opened from the region that seismic wave is detected first
Begin, sequentially determine whether the decay of spot seismic wave in adjacent area.Therefore, though being on a grand scale of Fault Rupture,
It also is able to estimate earthquake intensity.
Device operates
It follows that the operation that Fig. 4 will be used to describe earthquake intensity estimation unit 10 according to an embodiment of the invention.Figure
4 is the flow chart illustrating the operation of earthquake intensity estimation unit according to an embodiment of the invention.In the following description, will be regarding feelings
With reference to Fig. 1 and 2 depending on condition.In the present embodiment, earthquake intensity estimation side is implemented by operation earthquake intensity estimation unit 10
Method.Therefore, the earthquake intensity according to the present embodiment is replaced to estimate with the following description of the operation of earthquake intensity estimation unit 10
The description of method.
Receive seismic data and the step of source position information
As shown in Figure 4, initially, in earthquake intensity estimation unit 10, from seismic detector 20 one of communication unit 16
Receive seismic data, and receive source position information (step A1) from seismism observation system 40 further.Communication unit
Received source position information is also exported focus and determines unit 12 by unit 16, and by received seismic data
Output is to adjacent area designating unit 12.Note, if communication unit 16 not yet receives source position information and earthquake wave number
According to, then communication unit 16 enters resting state.
Determine the step of hypocentral location
Occur it follows that focus determines that unit 11 determines based on the source position information received in step A1
Whether the position of the focus 23 of earthquake is in setting district 30 (step A2).If the result of the determination in step A2 is focus 23
Position not in setting district 30, then the process in earthquake intensity estimation unit 10 terminates.On the other hand, if in step A2
The result determined be the position of focus 23 be in setting district 30, then focus determines that unit 11 notifies adjacent area designating unit 12
To perform step A3.
Adjacent area given step
It follows that when determining that unit 11 receives notice from focus, adjacent area designating unit 12 is specified with the most defeated
Go out the adjacent area (step A3) that region that the seismic detector 20 of seismic data is installed in is adjacent.Such as, in fig. 2, as
First the seismic detector 20 that fruit is installed in region a exports seismic data, then adjacent area designating unit 12 specifies region B conduct
Adjacent area.
Determine the step of decay
It follows that decay determine unit 13 via communication unit 16 from the adjacent area being arranged on specified by step A2
Seismic detector 20 obtains seismic data.Decay determines that unit 13 is then based on pre-conditioned determining whether by accessed
Amplitude specified by seismic data is decayed (step A4).
Now step A4 is described by use Fig. 5 A and 5B.Fig. 5 A is to be shown in which to determine not yet in step A3
The figure of the exemplarily seismic wave in the case of decaying, and Fig. 5 B is to be shown in which to determine in step A3 to decline
The figure of the exemplarily seismic wave in the case of subtracting.
In step A4, such as, if the amplitude of S ripple such as after the generation of earthquake 30 seconds or later default time
Between point (hereinafter referred to as " determining time point ") place or afterwards more than or equal to threshold value, then decay determines that unit 13 determines still
Decay, as shown in Figure 5 A.On the other hand, as shown in Figure 5 B, if the amplitude of S ripple is determining at time point or it
Afterwards less than threshold value, then decay determines unit 13 to determine to decay.Note, although threshold value is set in Fig. 5 A and 5B
The 80% of big amplitude, but the present embodiment is not limited to this.Threshold value can be based not only on the percentage ratio of peak swing and be set, and
And can be set by absolute value based on amplitude, or can the group of absolute value of percentage ratio based on peak swing and amplitude
Incompatible it is set.
Decay determines that the time can also be used as the most to have occurred the coefficient of carried out determination for decay by unit 13.?
In this case, if the amplitude of seismic wave sets the time period more than or equal to threshold value down to continuing one, then decay determines unit
13 determine that decay not yet occurs.
In the present embodiment, multiple seismic detector 20 is installed in each area, as described above.But,
If the number of seismic detector 20 is too big, then the processing speed of earthquake intensity estimation unit 20 will be likely to reduced.For this reason,
In step A4, after getting seismic data from the seismic detector 20 being arranged in adjacent area, decay determines unit 13
It is obtained in that by the meansigma methods of the amplitude specified by accessed seismic data.In this case, decay determines unit
13 can use obtained meansigma methods to determine whether decays.
If the result of the determination in step A4 is to decay to occur, then the knot of the process in earthquake intensity estimation unit 10
Bundle.On the other hand, if the result of the determination in step A4 is to decay not yet to occur, then Fault Rupture site setting unit 14 is held
Row step A5.
Set the step in Fault Rupture place
In step A5, Fault Rupture site setting unit 14 position based on the adjacent area specified in step A2,
In setting district 30, set the Fault Rupture place 24(estimating to occur Fault Rupture see Fig. 2).In the present embodiment, setting district 30
Be set at the border 22 along continental plate band-shaped in.Therefore, in step A5, such as, Fault Rupture site setting list
Fig. 2 sees at the border 22(of continental plate in unit 14) above specify the place closest to adjacent area, and by this site setting
For Fault Rupture place 24.
Calculate the step estimating earthquake intensity
It follows that estimate that earthquake intensity computing unit 15 comes centered by the Fault Rupture place 24 set in step A5
Calculate and estimate earthquake intensity (step A6).Figure is seen with the focus 23(of earthquake if not yet calculated before the execution of step A6
2) the estimation earthquake intensity centered by, then estimate that earthquake intensity computing unit 15 can also be with the focus 23 of earthquake in step A6
Estimation earthquake intensity is calculated centered by (see figure 2).
In the present embodiment, it is used for calculating and estimates that the method for earthquake intensity includes non-patent file 1 above-mentioned
Disclosed in computational methods, and such as, estimate that earthquake intensity computing unit 15 uses table below to reach formula 1 to 4 and calculates and estimate
Meter earthquake intensity I.
In following expression formula 1, PGV600It is with reference to the maximal rate [m/ under basis (the S wave velocity of 600m/s)
S], Mw is instantaneous earthquake magnitude, and D is the degree of depth (if the degree of depth of focus is 10kM or the least ,=0) of focus.At following table
Reach in formula 1, the distance between the place that x is focus 23 or Fault Rupture place 24 and seismic detector 20 is installed in.Note, make
For seismic detector 20 in this case, to depend on the circumstances and select optimal seismic detector 20, such as be positioned in the zone
Seismic detector at the heart or the seismic detector on the border 22 closest to continental plate.
In following expression formula 2, Mjma is " Japan Meteorological Agency earthquake magnitude " (the non-patent file 1 see mentioned above).
Additionally, in expression formula 3 and 4, PGV is the maximal rate [m/s] on the ground surface at each place.In expression formula 3,
ARVi is the speed amplification coefficient on the ground surface at each place.
[expression formula 1]
log(PCV600)=0.58Mw+0.0038D-1.29-log(Rm+0.0028*100.50M)=0.002x
[expression formula 2]
Mw=Mjma-0.171
[expression formula 3]
PGV=ARVi×0.90×PGV600
[expression formula 4]
I=2.68+1.721og(PGV)±0.21
The step of output estimation earthquake intensity
It follows that estimate that the estimation earthquake intensity calculated is exported seismism observation by earthquake intensity computing unit 15
System 40(step A7).Therefore, the estimation earthquake intensity calculated is published to each place by seismism observation system 40.
If there is Fault Rupture, then bigger region is calculated described estimation earthquake intensity, and therefore, seismism is observed
System 40 can issue the estimation earthquake intensity of larger area.
After step A7 is performed, step A3 is executed once again.In step A3 in this case, adjacent area refers to
Cell 12 specifies new adjacent area based on specified adjacent area.Such as, if region B has been designated as adjacent region
Territory, then specify the region C adjacent with region B.
Thereafter, step A4 is performed, and decay determines that unit 13 determines whether in the adjacent area reassigned
Seismic wave is decayed.If decay occurs, then step A5, A6 and A7 are further executed.
In other words, in the present embodiment, step A3 to A7 is repeatedly performed until determining in step A4 and having occurred
Till decay.For this reason, from the beginning of the region that seismic wave is detected first, in respective regions one by one
Ground, it is determined whether decay in seismic wave.Therefore, according to the present embodiment, extensive disconnected even for being attended by generation
The earthquake that layer ruptures, it is also possible to suitably estimate earthquake intensity.
Program
Program according to the present embodiment is configured to make computer perform the step A1 to A7 being illustrated in the diagram.Pass through
This program is installed on computers and performs installed program and realize the earthquake intensity estimation dress according to the present embodiment
Put 10 and earthquake intensity method of estimation.In this case, the CPU(CPU of computer) determining list as focus
Unit 11, adjacent area designating unit 12, decay determine unit 13, Fault Rupture site setting unit 14 and estimate earthquake intensity
Process is performed while computing unit 15.
Additionally, the computer being provided with the program according to the present embodiment can be to constitute existing seismism observation system
The computer of 40.In this case, seismism will be established as according to the earthquake intensity estimation unit 10 of the present invention to see
A part for examining system 40.
Here, earthquake intensity estimation dress is realized by execution according to the program of the present embodiment by using Fig. 6 to describe
Put the computer of 10.Fig. 6 is the example calculation illustrating and realizing earthquake intensity estimation unit 10 according to an embodiment of the invention
The block diagram of machine.
As shown in Figure 6, computer 110 is equipped with CPU111, main storage 112, storage device 113, input interface
114, display controller 115, data reader/write device 116 and communication interface 117.These unit via bus 121 by with
The mode allowing for data communication is connected to each other.
CPU111 is by extending the program according to the present embodiment of storage in storage device 113 in main storage 112
(code) and perform these codes to perform various types of operation with predefined procedure.Main storage 12 is the most such as
DRAM(dynamic random access memory) volatile storage.Further, it is being stored in meter according to the program of the present embodiment
It is provided under state in calculation machine readable medium recording program performing 120.Note, can be distributed to via logical according to the program of the present embodiment
Letter interface 117 connected on the Internet.
In addition to a hard disk, the semiconductor storage of such as flash memory or the like is given as storing device
The particular example of 113.Input interface 114 is reconciled the data between the input equipment 118 of CPU111 and such as keyboard and mouse and is led to
Letter.Display controller 115 is connected to display device 119, and controls the display in display device 119.Data reader/write
Enter device 116 and reconcile the data communication between CPU111 and record medium, and perform program from the reading recorded medium 120
And perform to be written to record medium 120 by the result handled by computer 110.Communication interface 117 reconcile CPU111 with
Data communication between another computer.
Such as CF(Compact Flash(registered trade mark)) or SD(Secure Digital) general semiconductor storage dress
Put, the magnetic storage medium of such as flexible disk or such as CD-ROM(Compact Disk Read Only Memory) optics deposit
Storage media is given as recording the particular example of medium 120.
Embodiment described above can fully or partly be represented to 12 by note 1 described below,
But it is not limited to following description.
(complementary annotations 1)
A kind of earthquake intensity estimation unit, comprising:
Focus determines unit, when in each region in the multiple regions set along the border of particular continent plate
The seismic detector installed detects by the seismic wave caused by the generation of earthquake and exports the amplitude for specifying described seismic wave
Data time, described focus determines that unit determines whether the position of the focus of occurred earthquake is in preset areas;
By described focus, adjacent area designating unit, if determining that unit determines that the position of described focus is in described district
In, the most described adjacent area designating unit specify with first it is detected that the district that is installed in of the described seismic detector of described seismic wave
The adjacent area that territory is adjacent;
Decay determine unit, described decay determine unit from the described seismic detector that is arranged on described adjacent area obtain by
This output data, and based on pre-conditioned determine the most send out in by the amplitude specified by accessed data
Raw decay;
Fault Rupture site setting unit, decays if determining unit to determine not yet by described decay, then described
Fault Rupture site setting unit position based on described adjacent area sets in described district wherein to be estimated Fault Rupture
Place;And
Estimating earthquake intensity computing unit, described estimation earthquake intensity computing unit is to set by described Fault Rupture place
To estimating that earthquake intensity calculates centered by place set by cell.
(complementary annotations 2)
In the earthquake intensity estimation unit according to complementary annotations 1, with the border along described continental plate band-shaped come
Set described district.
(complementary annotations 3)
In the earthquake intensity estimation unit according to complementary annotations 2,
Described adjacent area designating unit specifies new adjacent area based on described adjacent area,
Described decay determines what unit thus exported from the described seismic detector acquisition being arranged on described new adjacent area
Data, and based on described pre-conditioned determine the most decline in by the amplitude specified by accessed data
Subtract;
If it is determined that not yet decay, the most described Fault Rupture site setting unit is based on described new adjacent area
Position, resets the place wherein estimating to occur Fault Rupture on the border of described continental plate, and
Described estimation earthquake intensity computing unit is with the place reset by described Fault Rupture site setting unit
Centered by calculate further described estimation earthquake intensity.
(complementary annotations 4)
In the earthquake intensity estimation unit according to any one in complementary annotations 1 to 3,
Each region in the plurality of region is installed multiple seismic detector, and
Described decay determines that unit obtains described data from each described seismic detector being arranged on the adjacent area specified,
Obtain by the meansigma methods of the described amplitude specified by the data of each acquisition, and use obtained meansigma methods to determine to be
No decay.
(complementary annotations 5)
In the earthquake intensity estimation unit according to any one in complementary annotations 1 to 4,
If by the described amplitude specified by accessed data more than or equal to when setting threshold value so that persistently setting
Between section, the most described decay determines unit to determine not yet to decay.
(complementary annotations 6)
A kind of earthquake intensity method of estimation, it includes following steps:
A () is when the ground installed in each region in the multiple regions set by the border along particular continent plate
Shake instrument detect by the seismic wave caused by the generation of earthquake and export for the data of the amplitude specifying described seismic wave time,
Whether the position determining the focus of occurred earthquake is in preset areas;
If b () determines that the position of described focus is in described district in step (a)., then specify and first it is detected that institute
State the adjacent area that region that the described seismic detector of seismic wave is installed in is adjacent;
C () obtains, from the described seismic detector being arranged on described adjacent area, the data thus exported, and based on presetting
Condition determines and the most decays in by the described amplitude specified by accessed data;
Not yet decay if d () determines in step (c), then position based on described adjacent area, in described district
Interior setting wherein estimates to occur the place of Fault Rupture;And
E estimation earthquake intensity, centered by the place being set in step (d), is calculated by ().
(complementary annotations 7)
In the earthquake intensity method of estimation according to complementary annotations 6, with the border along described continental plate band-shaped come
Set described district.
(complementary annotations 8)
In the earthquake intensity method of estimation according to complementary annotations 7,
In step (b), specify new adjacent area based on described adjacent area,
In step (c), obtain from the described seismic detector being arranged on described new adjacent area and exported by this seismic detector
Data, and based on described pre-conditioned determine the most occur in by the amplitude specified by accessed data
Decay;
In step (d), if determining in step (c) and not yet decaying, then position based on new adjacent area,
The border of continental plate resets the place wherein estimating to occur Fault Rupture, and
In process step (e), centered by the place reset in step (d), to estimating that earthquake intensity is entered further
Row calculates.
(complementary annotations 9)
In the earthquake intensity method of estimation according to any one in complementary annotations 6 to 8,
Each in the plurality of region is provided with multiple seismic detector, and
In step (c), obtain data from each described seismic detector being arranged on specified adjacent area, it is thus achieved that logical
Cross the meansigma methods of amplitude specified by the data of each acquisition, and use obtained meansigma methods to determine whether to decline
Subtract.
(complementary annotations 10)
In the earthquake intensity method of estimation according to any one in complementary annotations 6 to 9,
In step (c), if by the amplitude specified by accessed data more than or equal to setting threshold value so that holding
The continuous setting time period, it is determined that decay not yet occurs.
(complementary annotations 11)
A kind of computer readable recording medium storing program for performing having program stored therein, described program comprises for making computer perform following
The order of step:
A () is when the ground installed in each region in the multiple regions set by the border along particular continent plate
Shake instrument detect by the seismic wave caused by the generation of earthquake and export for the data of the amplitude specifying described seismic wave time,
Whether the position determining the focus of occurred earthquake is in preset areas;
If b () determines that the position of described focus is in described district in step (a)., then specify and first it is detected that institute
State the adjacent area that region that the described seismic detector of seismic wave is installed in is adjacent;
C () obtains, from the described seismic detector being arranged on described adjacent area, the data thus exported, and based on presetting
Condition determines and the most decays in by the described amplitude specified by accessed data;
Not yet decay if d () determines in step (c), then position based on described adjacent area, in described district
Interior setting wherein estimates to occur the place of Fault Rupture;And
E estimation earthquake intensity, centered by the place being set in step (d), is calculated by ().
(complementary annotations 12)
In the computer readable recording medium storing program for performing according to complementary annotations 11, band-shaped with the border along described continental plate
Set described district.
(complementary annotations 13)
In the computer readable recording medium storing program for performing according to complementary annotations 12,
In step (b), specify new adjacent area based on described adjacent area,
In step (c), obtain from the described seismic detector being arranged on described new adjacent area and exported by this seismic detector
Data, and based on described pre-conditioned determine the most occur in by the amplitude specified by accessed data
Decay;
In step (d), if determining in step (c) and not yet decaying, then position based on new adjacent area,
The border of continental plate resets the place wherein estimating to occur Fault Rupture, and
In process step (e), centered by the place reset in step (d), to estimating that earthquake intensity is entered further
Row calculates.
(complementary annotations 14)
In the computer readable recording medium storing program for performing according to any one in complementary annotations 11 to 13,
Each in the plurality of region is provided with multiple seismic detector, and
In step (c), obtain data from each described seismic detector being arranged on specified adjacent area, it is thus achieved that logical
Cross the meansigma methods of amplitude specified by the data of each acquisition, and use obtained meansigma methods to determine whether to decline
Subtract.
(complementary annotations 15)
In the computer readable recording medium storing program for performing according to any one in complementary annotations 11 to 14,
In step (c), if by the amplitude specified by accessed data more than or equal to setting threshold value so that holding
The continuous setting time period, it is determined that decay not yet occurs.
Although describing invention required for protection above with reference to embodiment, but invention required for protection being not limited to
Above-described embodiment.It will be appreciated by those of skill in the art that, can be to claimed in the range of invention required for protection
The configuration of invention and details make various amendment.
The application based on and require in the existing Japanese patent application No.2011-208608 that JIUYUE in 2011 is submitted on the 26th
The rights and interests of priority, its entire disclosure is incorporated herein.
Industrial usability
According to the present invention, even if in the case of subduction zone earthquake wherein occurs, it is also possible to reliably specify and wherein tremble
The region being predicted to, and the earthquake intensity in this region can be estimated.The present invention is applicable to when earthquake occurs
Time need to issue as quickly as possible the system of the information about earthquake.
Description of reference numerals
10 earthquake intensity estimation units
11 focus determine unit
12 adjacent area designating unit
13 decay determine unit
14 Fault Rupture site setting unit
15 estimate earthquake intensity computing unit
16 communication units
20 seismic detectors
21a to 21d seismic detector group
The border of 22 continental plates
23 focus
24 Fault Rupture places
30 specify district
40 seismism observation systems
50 networks
110 computers
111 CPU
112 main storages
113 storage devices
114 input interfaces
115 display controllers
116 data readers/write device
117 communication interfaces
118 input equipment
119 display devices
120 record media
121 buses
Claims (6)
1. an earthquake intensity estimation unit, comprising:
Focus determines unit, when installing in each region in the multiple regions set along the border of particular continent plate
Seismic detector detect by the seismic wave caused by the generation of earthquake and the number that exports amplitude for specifying described seismic wave
According to time, described focus determines that unit determines whether the position of the focus of occurred earthquake is in preset areas;
Adjacent area designating unit, if determining that unit determines that the position of described focus is in described district by described focus,
The most described adjacent area designating unit specify with first it is detected that the region that is installed in of the described seismic detector of described seismic wave
Adjacent adjacent area;
Decay determines that unit, described decay determine that unit obtains the most defeated from the described seismic detector being arranged on described adjacent area
The data gone out, and based on pre-conditioned determine the most decline in by the amplitude specified by accessed data
Subtract;
Fault Rupture site setting unit, decays if determining unit to determine not yet by described decay, the most described tomography
Rupture site setting unit position based on described adjacent area in described district, set the ground wherein estimating to occur Fault Rupture
Point;And
Estimating earthquake intensity computing unit, described estimation earthquake intensity computing unit is with by described Fault Rupture site setting list
Centered by place set by unit, to estimating that earthquake intensity calculates.
Earthquake intensity estimation unit the most according to claim 1, wherein,
To set described district along the band-shaped of the border of described continental plate.
Earthquake intensity estimation unit the most according to claim 2, wherein,
Described adjacent area designating unit specifies new adjacent area based on described adjacent area,
Described decay determines that unit obtains, from the described seismic detector being arranged on described new adjacent area, the data thus exported,
And based on described pre-conditioned determine the most decay in by the amplitude specified by accessed data;
If it is determined that not yet decay, the most described Fault Rupture site setting unit position based on described new adjacent area
Put, the border of described continental plate resets the place wherein estimating to occur Fault Rupture, and
During described estimation earthquake intensity computing unit with the place reset by described Fault Rupture site setting unit is
The heart, calculates described estimation earthquake intensity further.
Earthquake intensity estimation unit the most according to claim 1, wherein,
Each region in the plurality of region is provided with multiple seismic detector, and
Described decay determines that unit, obtains to obtain described data from each described seismic detector of installation the adjacent area specified
By the meansigma methods of the described amplitude specified by each obtained data, and obtained meansigma methods must be used to determine whether
Decay.
Earthquake intensity estimation unit the most according to claim 1,
If by the described amplitude specified by accessed data more than or equal to setting threshold value so that persistently setting the time period,
The most described decay determines unit to determine not yet to decay.
6. an earthquake intensity method of estimation, it includes following steps:
A () is when the seismic detector installed in each region in the multiple regions set by the border along particular continent plate
Detect by the seismic wave caused by the generation of earthquake and export for the data of the amplitude specifying described seismic wave time, determine
Whether the position of the focus of the earthquake occurred is in preset areas;
If b () determines that in step (a) position of described focus is in described district, then specify with first it is detected that describedly
The adjacent area that region that the described seismic detector of seismic wave is installed in is adjacent;
C () obtains, from the described seismic detector being arranged on described adjacent area, the data thus exported, and based on pre-conditioned
Determine and the most decay in by the described amplitude specified by accessed data;
Not yet decay if d () determines in step (c), then position based on described adjacent area, sets in described district
The fixed place wherein estimating Fault Rupture;And
E estimation earthquake intensity, centered by the place being set in step (d), is calculated by ().
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011-208608 | 2011-09-26 | ||
JP2011208608 | 2011-09-26 | ||
PCT/JP2012/071065 WO2013047008A1 (en) | 2011-09-26 | 2012-08-21 | Seismic intensity estimation device, seismic intensity estimation method, and computer-readable recording medium |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103827694A CN103827694A (en) | 2014-05-28 |
CN103827694B true CN103827694B (en) | 2016-11-02 |
Family
ID=47995059
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201280046946.0A Expired - Fee Related CN103827694B (en) | 2011-09-26 | 2012-08-21 | Earthquake intensity estimation unit, earthquake intensity method of estimation and computer readable recording medium storing program for performing |
Country Status (3)
Country | Link |
---|---|
JP (1) | JP5742951B2 (en) |
CN (1) | CN103827694B (en) |
WO (1) | WO2013047008A1 (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101697227B1 (en) * | 2016-05-23 | 2017-01-17 | 한국지질자원연구원 | Method for determining epicenter location in order to prevent false alarms of an earthquake early warning system using forced association of adjacent observation station information |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002040152A (en) * | 2000-07-28 | 2002-02-06 | Nec Ocean Eng Ltd | Seabed observation system |
JP2006078297A (en) * | 2004-09-09 | 2006-03-23 | Railway Technical Res Inst | Early instrumental seismic intensity prediction method and apparatus for the same |
JP2007040949A (en) * | 2005-07-07 | 2007-02-15 | Taisei Corp | System for predicting seismic motion utilizing real-time seismic information |
JP2007108012A (en) * | 2005-10-13 | 2007-04-26 | Oki Electric Ind Co Ltd | Earthquake disaster prevention system |
CN102156298A (en) * | 2011-03-11 | 2011-08-17 | 中国地震局地壳应力研究所 | Rapid seismic intensity assessment method based on non-statistical hypothesis test |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4385233B2 (en) * | 2007-07-31 | 2009-12-16 | 日本電気株式会社 | Earthquake prediction information notification server, earthquake prediction information notification system, earthquake prediction information notification method and program |
JP2013007728A (en) * | 2011-06-22 | 2013-01-10 | Home Seismometer:Kk | Real time estimation method of hypocentral region of giant earthquake |
JP5770051B2 (en) * | 2011-09-05 | 2015-08-26 | 公益財団法人鉄道総合技術研究所 | Immediate estimation method of the rupture area of a large earthquake |
-
2012
- 2012-08-21 JP JP2013536076A patent/JP5742951B2/en active Active
- 2012-08-21 WO PCT/JP2012/071065 patent/WO2013047008A1/en active Application Filing
- 2012-08-21 CN CN201280046946.0A patent/CN103827694B/en not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002040152A (en) * | 2000-07-28 | 2002-02-06 | Nec Ocean Eng Ltd | Seabed observation system |
JP2006078297A (en) * | 2004-09-09 | 2006-03-23 | Railway Technical Res Inst | Early instrumental seismic intensity prediction method and apparatus for the same |
JP2007040949A (en) * | 2005-07-07 | 2007-02-15 | Taisei Corp | System for predicting seismic motion utilizing real-time seismic information |
JP2007108012A (en) * | 2005-10-13 | 2007-04-26 | Oki Electric Ind Co Ltd | Earthquake disaster prevention system |
CN102156298A (en) * | 2011-03-11 | 2011-08-17 | 中国地震局地壳应力研究所 | Rapid seismic intensity assessment method based on non-statistical hypothesis test |
Non-Patent Citations (2)
Title |
---|
中国地震台网震级的对比;刘瑞丰 等;《地震学报》;20080930;第29卷(第5期);第467-476页 * |
华北地区烈度衰减模型建立及其用于震中区域和震级的定量估算;张杨 等;《地震学报》;20090531;第31卷(第3期);第290-306页 * |
Also Published As
Publication number | Publication date |
---|---|
JP5742951B2 (en) | 2015-07-01 |
CN103827694A (en) | 2014-05-28 |
JPWO2013047008A1 (en) | 2015-03-26 |
WO2013047008A1 (en) | 2013-04-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Brown et al. | Development of the ElarmS methodology for earthquake early warning: Realtime application in California and offline testing in Japan | |
Wu et al. | Relationship between peak ground acceleration, peak ground velocity, and intensity in Taiwan | |
Limongelli | Seismic health monitoring of an instrumented multistory building using the interpolation method | |
CN111160125A (en) | Railway foreign matter intrusion detection method based on railway monitoring | |
CN103686065A (en) | Cloud mirror cluster control method and device of monitoring equipment based on GIS (geographic information system) interoperability | |
CN106033514A (en) | Suspicious process detection method and device | |
US20120134527A1 (en) | Hazard detection for asset management | |
CN108176050A (en) | Method for searching and device | |
US20220120607A1 (en) | Optical fiber sensing system, monitoring apparatus, monitoring method, and computer readable medium | |
CN105867899A (en) | Method and device for identifying device | |
Xu et al. | Reliability assessment on earthquake early warning: A case study from Taiwan | |
JP2006208249A (en) | Earthquake disaster mitigation system, earthquake disaster mitigation method, earthquake disaster mitigation program, and computer-readable storage medium | |
Pollino et al. | An open source GIS system for earthquake early warning and post-event emergency management | |
CN103827694B (en) | Earthquake intensity estimation unit, earthquake intensity method of estimation and computer readable recording medium storing program for performing | |
Premus | Modal scintillation index: A physics-based statistic for acoustic source depth discrimination | |
CN110392239A (en) | Specified area monitoring method and device | |
CN108304578B (en) | Map data processing method, medium, device and computing equipment | |
CN112396021B (en) | High-altitude dangerous object detection method, device and storage medium | |
Yamada et al. | The extended integrated particle filter method (IPFx) as a high‐performance earthquake early warning system | |
EP3407086A1 (en) | Mobile object measurement system and method of determining number of people in measurement area | |
US11841305B2 (en) | System, device and method for effective deployment of a dust accumulation sensor | |
CN111382643B (en) | Gesture detection method, device, equipment and storage medium | |
Perkins et al. | The annotation approach used for marine imagery impacts the detection of temporal trends in seafloor biota | |
CN106845214A (en) | Based on safety protecting method and system under virtualized environment | |
CN104349348B (en) | Problem takes place frequently the determination method, apparatus and server in region in drive test |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: DE Ref document number: 1195813 Country of ref document: HK |
|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: WD Ref document number: 1195813 Country of ref document: HK |
|
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20161102 Termination date: 20210821 |
|
CF01 | Termination of patent right due to non-payment of annual fee |