CN103827694A - Seismic intensity estimation device, seismic intensity estimation method, and computer-readable recording medium - Google Patents
Seismic intensity estimation device, seismic intensity estimation method, and computer-readable recording medium Download PDFInfo
- Publication number
- CN103827694A CN103827694A CN201280046946.0A CN201280046946A CN103827694A CN 103827694 A CN103827694 A CN 103827694A CN 201280046946 A CN201280046946 A CN 201280046946A CN 103827694 A CN103827694 A CN 103827694A
- Authority
- CN
- China
- Prior art keywords
- adjacent area
- seismic intensity
- unit
- data
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims description 23
- 230000001143 conditioned effect Effects 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 abstract description 3
- 230000000295 complement effect Effects 0.000 description 27
- 238000004891 communication Methods 0.000 description 19
- 238000010586 diagram Methods 0.000 description 6
- 230000015654 memory Effects 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 230000005055 memory storage Effects 0.000 description 5
- 238000003860 storage Methods 0.000 description 5
- 206010044565 Tremor Diseases 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- 238000003325 tomography Methods 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
- 230000005540 biological transmission Effects 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 239000012925 reference material Substances 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
Images
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 seismic intensity estimation device (10) is provided with: a seismic source determination unit (11) which, when seismographs (20) installed in respective regions along the boundary of a continental plate output data by which the amplitude of a seismic wave can be specified, determines whether the position of a seismic source is within a set region; an adjacent area specification unit (12) which, when the position of the seismic source is within the set region, specifies an adjacent area adjacent to an area where the seismograph that detected the seismic wave first is installed; an attenuation determination unit (13) which acquires data from the seismograph in the adjacent area and determines whether attenuation is occurring in the amplitude of the data; a fault rupture point setting unit (14) which, when the attenuation does not occur, sets a fault rupture estimation point on the set region on the basis of the position of the adjacent area; and an estimated seismic intensity calculation unit (15) which calculates an estimated seismic intensity with the fault rupture estimation point as the center.
Description
Technical field
The present invention relates to a kind of for the seismic intensity estimation unit of the seismic intensity in the each place of stage forecast and seismic intensity method of estimation in early days in the time that earthquake occurs, and a kind of computer readable recording medium storing program for performing that stores the program for realizing described device and described method.
Background technology
In recent years, in Japan, seismism recording geometry has been introduced into so that promptly to society's notice earthquake information.In the time that earthquake occurs, the observation data that the analysis of seismism recording geometry is immediately captured by near the seismograph being positioned at focus after occurring, estimate focus and earthquake magnitude, and the seismic intensity in the each place of the further estimation of the focus based on estimated and earthquake magnitude and the time of arrival (for example, seeing non-patent file 1) of main motion.Now, generally speaking, until as the observation data of analytic target being the P ripple initial motion in stage in early days in the time that earthquake occurs.
Then seismism recording geometry is issued estimated information (focus, earthquake magnitude, seismic intensity and time of arrival) widely via various media as the early stage alarm of earthquake.As a result, people can promptly take action to evacuate in factory, office and residence, and therefore reduce by the caused loss of 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, the meteorological Room, on July 29th, 2008,3-15 page
Summary of the invention
Problem to be solved by this invention
Incidentally, in non-patent file 1 above-mentioned disclosed seismism recording geometry to get focus be place and estimate seismic intensity centered by this place.For this reason, in non-patent file 1 above-mentioned disclosed seismism recording geometry near its focus, cause tomography and coverage relatively little earthquake be effectively, but there is following problem with respect to the subduction zone earthquake of adjoining land generation tomography on making on a large scale.
First, subduction zone earthquake causes the Fault Rupture starting from focus on may be on a large scale.In this case, seismic event can be propagated on the scope than wider around the concentric scope of focus.For this reason, compared with the scope of the seismic intensity that in non-patent file 1, disclosed seismism recording geometry predicts, macroseism is quivered and will on wider scope, be occurred.Therefore, not yet utilize the early stage alarm prediction of earthquake will get an unexpected blow to the resident in the region of seismic intensity, and loss may increase.
Exemplary purpose of the present invention is to solve foregoing problems, and provide a kind of seismic intensity estimation unit, seismic intensity method of estimation and computer readable recording medium storing program for performing, even if it is in the situation that subduction zone earthquake occurs, also can specify reliably the region of predicting trembling, and make it possible to estimate the seismic intensity in this region.
For the means of dealing with problems
In order to realize the object of stating above, seismic intensity estimation unit is according to an aspect of the present invention following seismic intensity estimation unit, and it comprises:
Focus determining unit, when the seismograph of installing in each region in multiple regions of setting on border along particular continent plate detects that while being used to specify the data of described earthquake wave amplitude by the caused seismic event of generation of earthquake and output, described focus determining unit determines whether the position of the focus of the earthquake occurring is in preset areas;
Adjacent area designating unit, if determine that by described focus determining unit the position of described focus is in described district, described adjacent area designating unit is specified the adjacent adjacent area in region being installed in the described seismograph that described seismic event first detected;
Decay determining unit, described decay determining unit is obtained the data of output thus from being arranged on described seismograph described adjacent area, and determines in the specified amplitude of the data by accessed whether decay based on pre-conditioned;
Fault Rupture site setting unit, if determined and not yet decayed by described decay determining unit, the place of wherein estimating to occur Fault Rupture is set in the position of described Fault Rupture site setting unit based on described adjacent area in described district; And
Estimate seismic intensity computing unit, described estimation seismic intensity computing unit calculates estimation seismic intensity centered by the place setting by described Fault Rupture site setting unit.
In addition, in order to realize the object of stating above, seismic intensity method of estimation is according to an aspect of the present invention following seismic intensity method of estimation, and it comprises following steps:
(a) when the seismograph of installing in each region in the multiple regions that set along the border of particular continent plate detects while being used to specify the data of described earthquake wave amplitude by the caused seismic event of generation of earthquake and output, determine whether the position of the focus of the earthquake occurring is in preset areas;
(b), if determine that in step (a) position of described focus is in described district, specify the adjacent adjacent area in region being installed in the described seismograph that described seismic event first detected;
(c) obtain the data of output thus from being arranged on described seismograph described adjacent area, and determine in the specified described amplitude of the data by accessed whether decay based on pre-conditioned;
(d), if determine and not yet decay in step (c), the place of wherein estimating to occur Fault Rupture is set in the position based on described adjacent area in described district; And
(e), centered by the place being set in step (d), estimation seismic intensity is calculated.
And in order to realize the target of stating above, computer readable recording medium storing program for performing is according to an aspect of the present invention following computer readable recording medium storing program for performing, it has program stored therein, and described routine package is containing being useful on the order that makes computing machine carry out following steps:
(a) when the seismograph of installing in each region in the multiple regions that set along the border of particular continent plate detects while being used to specify the data of described earthquake wave amplitude by the caused seismic event of generation of earthquake and output, determine whether the position of the focus of the earthquake occurring is in preset areas;
(b), if determine that in step (a) position of described focus is in described district, specify the adjacent adjacent area in region being installed in the described seismograph that described seismic event first detected;
(c) obtain the data of output thus from being arranged on described seismograph described adjacent area, and determine in the specified described amplitude of the data by accessed whether decay based on pre-conditioned;
(d), if determine and not yet decay in step (c), the place of wherein estimating to occur Fault Rupture is set in the position based on described adjacent area in described district; And
(e), centered by the place being set in step (d), estimation seismic intensity is calculated.
Beneficial effect of the present invention
As described above, according to the present invention, even if in the situation that subduction zone earthquake occurs, also can specify reliably the region of predicting trembling, and can estimate the seismic intensity in this region.
Accompanying drawing explanation
Fig. 1 illustrates the block diagram of the configuration of seismic intensity estimation unit according to an embodiment of the invention.
Fig. 2 is the figure of the example of the state that illustrates that wherein the seismograph that uses is in an embodiment of the present invention mounted.
Fig. 3 illustrates figure Japanese and continental-bearing plate around it.
Fig. 4 illustrates the process flow diagram of the operation of seismic intensity estimation unit according to an embodiment of the invention.
Fig. 5 A is illustrated in the figure that has determined the exemplarily seismic wave in the situation that of not yet there is decay in steps A 3, and Fig. 5 B is illustrated in the figure that has determined the exemplarily seismic wave in the situation that of there is decay in steps A 3.
Fig. 6 illustrates to implement the block diagram of the illustrative computer of seismic intensity estimation unit 10 according to an embodiment of the invention.
Embodiment
Embodiment
Hereinafter, with reference to Fig. 1 to 6, seismic intensity estimation unit, seismic intensity method of estimation and program are according to an embodiment of the invention described.
Device configuration
First, will describe with Fig. 1 to 3 according to the configuration of the seismic intensity estimation unit 10 of the present embodiment.Fig. 1 illustrates the block diagram of the configuration of seismic intensity estimation unit according to an embodiment of the invention.Fig. 2 is the figure of the example of the state that illustrates that wherein the seismograph that uses is in an embodiment of the present invention mounted.Fig. 3 illustrates figure Japanese and continental-bearing plate around it.
As shown in fig. 1, seismic intensity estimation unit 10 is connected to multiple shake instrument 20 via network 50.In the present embodiment, seismic intensity estimation unit 10 is also connected to seismism recording geometry 40 via network 50.
Described in background technology, in the time that earthquake occurs, seismism recording geometry 40 specifies focus and the focus based on specified to calculate earthquake magnitude.The seismism recording geometry 40 also earthquake magnitude based on calculated is estimated seismic intensity in each place and the time of arrival of main motion.
The region that seismograph 20 is installed in by them is divided in groups.Reference numeral 21a to 21d represents each group of seismograph 20.In the time that earthquake occurs, each seismograph 20 detects seismic event, and output is used to specify the data (hereinafter referred to as " seismic data ") of this earthquake wave amplitude.In the present embodiment, seismic data can be the data of only specifying earthquake wave amplitude, or can be to specify the amplitude of seismic event and the data of waveform.
In the example of Fig. 2, it is in the A to D of region that seismograph 20 is installed in four regions.Setting regions A to D is carried out on border 22 along continental-bearing plate.In each in region A, B, C and D, the seismograph 20 of respective sets is installed.Note, Japan and its continental-bearing plate are around shown in Fig. 3.
In addition, as shown in fig. 1, seismic intensity estimation unit 10 comprises focus determining unit 11, adjacent area designating unit 12, decay determining unit 13, Fault Rupture site setting unit 14 and estimates seismic intensity computing unit 15.
If the seismograph 20 being arranged in each region detects that the focus determining unit 11 among the constituent components of enumerating determines whether the position of the focus (see figure 2) of the earthquake having occurred is (hereinafter referred to as " setting district " in preset areas by the caused seismic event of generation and the output seismic data of earthquake in the above; See Fig. 2) in 30.
If determined in the Shi setting district, position 30 of focus 23 by focus determining unit 11, adjacent area designating unit 12 is specified the adjacent region (hereinafter referred to as " adjacent area "), region being installed in the seismograph that seismic event detected at first.
Decay determining unit 13 obtains from this seismograph the seismic data of being exported by the seismograph being arranged on specified adjacent area, and determines in the specified amplitude of the seismic data by accessed whether decay based on pre-conditioned.
In addition, if determine that by decay determining unit 13 decay not yet occurs, the position of the adjacent area of Fault Rupture site setting unit 14 based on specified, wherein estimates in the interior setting in setting district 30 place (hereinafter referred to as " Fault Rupture place ") 24 that Fault Rupture occurs.
Estimate that seismic intensity computing unit 15 calculates estimation seismic intensity centered by the Fault Rupture place 24 setting.In the present embodiment, estimate that the estimation seismic intensity in calculated each place is outputed to seismism recording geometry 40 by seismic intensity computing unit 15.In this case, seismism recording geometry 40 is via various media, not only issue widely the estimation seismic intensity that seismism recording geometry 40 is calculated centered by focus 23, but also issue the estimation seismic intensity of calculating by seismic intensity measurement mechanism 10.
Therefore, in the present embodiment, whether seismic intensity estimation unit 10 definitely seismic wave has been attenuated in the adjacent area adjacent with the region of close focus.In this case, if the Fault Rupture starting from focus occurs on a large scale, in adjacent area, do not decay.Therefore, seismic intensity estimation unit 10 can determine whether Fault Rupture occurs, and depend on determine result on a large scale, calculates estimation seismic intensity.As described above, adopt seismic intensity estimation unit 10, even if occur therein in the situation of subduction zone earthquake, also can specify reliably the region of predicting to trembling, and can estimate the seismic intensity in this region.
Now will describe in more detail according to the configuration of the seismic intensity estimation unit 10 of the present embodiment.As shown in fig. 1, seismic intensity estimation unit 10 further comprises communication unit 16.Communication unit 16 via network 30 to seismograph 20 and seismism recording geometry 40 transmission information and from seismograph 20 and seismism recording geometry 40 reception information.
In the present embodiment, in the time that seismograph 20 is exported seismic data, seismism recording geometry 40 specifies focus 23(to see Fig. 2), and send the positional information (for example, specifying the information of longitude and the latitude of focus with it) about focus 23.Therefore, the positional information (hereinafter referred to as " source location information ") that communication 16 receives about focus 23, and outputed to focus determining unit 11.In the present embodiment, focus determining unit 11 has been output in source location information to carry out under its condition and has determined and process.Note, source location information can be used as the early stage alarm of earthquake and is sent out.
In the present embodiment, as shown in Figures 2 and 3, preferably, with setting by focus determining unit 11 setting district 30 of being used in determining with shape of the border 22 along continental-bearing plate, and especially, so that comprise border 22.This is because generally speaking, Fault Rupture is to occur near the border 22 of continental-bearing plate.
Fig. 3 shows the possible setting district 30a to 30f in Japan.Note, setting district 30a to 30b shown in Fig. 3 is only example, and in practice, be near continental-bearing plate based on existing subject area, set described setting district about result of the simulation of the generation of the statistics of the focus of earthquake in the past, earthquake etc.Although near the region shown in the example in Fig. 3 Japan, the present embodiment is equally applicable to the country except Japan.
In the present embodiment, adjacent area designating unit 12 can also be specified new adjacent area by the adjacent area based on specified.In this case, the definite information 13 of decay is obtained the seismic data of the seismograph output being installed in new adjacent area from this seismograph.Then decay determining unit 13 decays in the specified amplitude of the seismic data by accessed based on pre-conditioned determining whether.
Further, if determine in new adjacent area and not yet decay, the position of the adjacent area based on new again, Fault Rupture site setting unit 14, sets on the border 22 of continental-bearing plate the place (Fault Rupture place) of wherein estimating to occur Fault Rupture.In this case, estimate that seismic intensity computing unit 15 further calculates estimation seismic intensity centered by the Fault Rupture place 24 being reset by Fault Rupture site setting unit 14.
According to above-described pattern, the region that first seismic intensity estimation unit 10 can be detected from seismic event, sequentially determine whether the decay of spot seismic wave in adjacent area.Therefore, even if being on a grand scale of Fault Rupture also can estimate seismic intensity.
Device operation
Next, will the operation of seismic intensity estimation unit 10 according to an embodiment of the invention be described with Fig. 4.Fig. 4 illustrates the process flow diagram of the operation of seismic intensity estimation unit according to an embodiment of the invention.In the following description, will depend on the circumstances with reference to Fig. 1 and 2.In the present embodiment, implement seismic intensity method of estimation by operation seismic intensity estimation unit 10.Therefore, replace according to the description of the seismic intensity method of estimation of the present embodiment with the following description of the operation of seismic intensity estimation unit 10.
Receive the step of seismic data and source location information
As shown in Figure 4, initial, in seismic intensity estimation unit 10, the reception seismic data of communication unit 16 from seismograph 20, and further receive source location information (steps A 1) from seismism recording geometry 40.Communication unit 16 also outputs to focus determining unit 12 by received source location information, and received seismic data is outputed to adjacent area designating unit 12.Note, if communication unit 16 not yet receives source location information and seismic data, communication unit 16 enters stand-by state.
Determine the step of source location
Next, the source location information of focus determining unit 11 based on receiving in steps A 1 determines whether the position of the focus 23 of the earthquake having occurred is in setting district 30 (steps A 2).If the definite result in steps A 2 is in the Bu setting district, position 30 of focus 23, the processing in focus estimation unit 10 finishes.On the other hand, if the definite result in steps A 2 is in the Shi setting district, position 30 of focus 23, focus determining unit 11 notifies adjacent area designating unit 12 with execution step A3.
Adjacent area given step
Next,, when receive notice from focus determining unit 11, adjacent area designating unit 12 is specified the adjacent adjacent area (steps A 3) in region being installed in the seismograph 20 of first exporting seismic data.For example, in Fig. 2, if first the seismograph 20 being arranged in the A of region exports seismic data, adjacent area designating unit 12 appointed area B are as adjacent area.
Determine the step of decay
Next, decay determining unit 13 is obtained seismic data via communication unit 16 from the seismograph 20 being arranged on the specified adjacent area of steps A 2.Decay determining unit 13 is then based on pre-conditioned determine whether decay in the specified amplitude of the seismic data by accessed (steps A 4).
Now use Fig. 5 A and 5B are described steps A 4.Fig. 5 A is illustrated in the figure that has wherein determined the exemplarily seismic wave in the situation that of not yet there is decay in steps A 3, and Fig. 5 B is illustrated in the figure that wherein determines the exemplarily seismic wave in the situation that of having there is decay in steps A 3.
In steps A 4, for example, if S wave amplitude 30 seconds or later Preset Time point (hereinafter referred to as " determining time point ") after the generation in earthquake are for example located or be more than or equal to afterwards threshold value, the determining unit that decays 13 determines and not yet decays, as shown in Figure 5 A.On the other hand, as shown in Figure 5 B, if S wave amplitude at definite time point place or be less than afterwards threshold value, the determining unit that decays 13 is determined and is decayed.Note, although threshold value is set to 80% of peak swing in Fig. 5 A and 5B, the present embodiment is not limited to this.The threshold value not only percentage based on peak swing is recently set, and can be set by the absolute value based on amplitude, or the combination of the absolute value of number percent that can be based on peak swing and amplitude is set.
Decay determining unit 13 can also be by the time as for decay, whether carried out definite coefficient having occurred.In this case, continue a setting-up time section if earthquake wave amplitude is more than or equal to threshold value, the determining unit that decays 13 determines that decay not yet occurs.
In the present embodiment, multiple seismographs 20 are installed in each region, as described above.But if the number of seismograph 20 is too large, the processing speed of seismic intensity estimation unit 20 may reduce.For this reason, in steps A 4, from being arranged on after seismograph 20 adjacent area gets seismic data, decay determining unit 13 can obtain the mean value of the amplitude specified by accessed seismic data.In this case, decay determining unit 13 can determine whether to decay with obtained mean value.
If the definite result in steps A 4 is to decay to occur, the processing in seismic intensity estimation unit 10 finishes.On the other hand, if the definite result in steps A 4 is to decay not yet to occur, Fault Rupture site setting unit 14 performs step A5.
Set the step in Fault Rupture place
In steps A 5, the position of the adjacent area of Fault Rupture site setting unit 14 based on specifying in steps A 2, in setting district 30, interior setting estimates that the Fault Rupture place 24(that Fault Rupture occurs is shown in Fig. 2).In the present embodiment, setting district 30 is set in the band shape on the border 22 of continental-bearing plate.Therefore, in steps A 5, for example, Fig. 2 is shown at the border of continental-bearing plate 22(in Fault Rupture site setting unit 14) upper specify the place close to adjacent area, and be Fault Rupture place 24 by this site setting.
Calculate the step of estimating seismic intensity
Next, estimate that seismic intensity computing unit 15 calculates estimation seismic intensity (steps A 6) centered by the Fault Rupture place 24 of setting in steps A 5.If not yet calculated before the execution of steps A 6 by the focus 23(of earthquake and see Fig. 2) centered by estimation seismic intensity, estimate seismic intensity computing unit 15 can also be in steps A 6 the focus 23(by earthquake see Fig. 2) centered by calculate estimation seismic intensity.
In the present embodiment, estimate that for calculating the method for seismic intensity comprises the disclosed computing method of non-patent file 1 above-mentioned, and for example, estimation seismic intensity computing unit 15 use expression formula 1 to 4 is below calculated and is estimated seismic intensity I.
In expression formula 1 below, PGV
600that Mw is instantaneous earthquake magnitude, and D is the degree of depth (if the degree of depth of focus is 10kM or less ,=0) of focus in the maximal rate [m/s] with reference under basis (the S wave velocity of 600m/s).In expression formula 1 below, x is the distance between the place that is installed in of focus 24 or Fault Rupture place 24 and seismograph 20.Note, as seismograph 20 in this case, depend on the circumstances to select optimal seismograph 20, such as being located at the seismograph at regional center place or close to the seismograph on the border 22 of continental-bearing plate.
In expression formula 2 below, Mjma is " Japan Meteorological Agency earthquake magnitude " (seeing the non-patent file 1 of mentioning above).In addition,, in expression formula 3 and 4, PGV is the maximal rate [m/s] on the ground surface at each place place.In expression formula 4, ARVi is the speed amplification coefficient on the ground surface at each place place.
[expression formula 1]
log(PCV
600)=0.58Mw+0.0038D-1.29-log(R
m+0.0028*10
0.50M)=0.002x
[expression formula 2]
Mw=Mjma-0.171
[expression formula 3]
PGV=ARVi×0.90×PGV
600
[expression formula 4]
I=2.68+1.721og(PGV)±0.21
The step of output estimation seismic intensity
Next, estimate that calculated estimation seismic intensity is outputed to seismism recording geometry 40(steps A 7 by seismic intensity computing unit 15).Therefore, calculated estimation seismic intensity is published to each place by seismism recording geometry 40.If there is Fault Rupture, described estimation seismic intensity is calculated in larger region, and therefore, seismism recording geometry 40 can be issued the estimation seismic intensity in larger region.
After steps A 7 is performed, steps A 3 is carried out again.In steps A 3 in this case, the adjacent area of adjacent area designating unit 12 based on specified specified new adjacent area.For example, if region B has been designated as adjacent area, specify the region C adjacent with region B.
Thereafter, steps A 4 is performed, and decay determining unit 13 determines whether to decay in the seismic event in the adjacent area of reassigning.If decay occurs, steps A 5, A6 and A7 are further carried out.
In other words, in the present embodiment, steps A 3 to A7 be repeatedly performed until in steps A 4 determine occurred decay.For this reason, the region first being detected from seismic event, in respective regions, one by one, determine whether to decay in seismic event.Therefore, according to the present embodiment, even if for being attended by the earthquake that extensive Fault Rupture occurs, also can suitably estimate seismic intensity.
Program
Be configured to make computing machine to carry out the steps A 1 being illustrated in Fig. 4 to A7 according to the program of the present embodiment.By by this installation on computers and carry out the program of installing and realize according to the seismic intensity estimation unit 10 of the present embodiment and seismic intensity method of estimation.In this case, the CPU(CPU (central processing unit) of computing machine) in being used as focus determining unit 11, adjacent area designating unit 12, decay determining unit 13, Fault Rupture site setting unit 14 and estimation seismic intensity computing unit 15, carry out processing.
In addition, being provided with according to the computing machine of the program of the present embodiment can be the computing machine that forms existing seismism recording geometry 40.In this case, Estimation of Seismetic Magnitude device 10 according to the present invention will be established as a part for seismism recording geometry 40.
To computing machine that realize seismic intensity estimation unit 10 according to the program of the present embodiment by carrying out be described with Fig. 6 here.Fig. 6 illustrates to realize the block diagram of the illustrative computer of seismic intensity estimation unit 10 according to an embodiment of the invention.
As shown in Figure 6, computing machine 110 is equipped with CPU111, primary memory 112, memory storage 113, input interface 114, display controller 115, data reader/write device 116 and communication interface 117.These unit are such that via bus 121 mode that can carry out data communication and are connected to each other.
CPU111 by primary memory 112 expansion in memory storage 113 storage according to the program of the present embodiment (code) and carry out these codes with predefined procedure and carry out various types of operations.Primary memory 12 is typically such as DRAM(dynamic RAM) volatile storage.And, be provided under the state being stored in computer readable recording medium storing program for performing 120 according to the program of the present embodiment.Note, according to the program of the present embodiment can be distributed to connect via communication interface 117 on the Internet.
Except hard disk, be presented the specific example as memory storage 113 such as the semiconductor storage of flash memory or analog.Input interface 114 is reconciled CPU111 and such as the data communication between the input media 118 of keyboard and mouse.Display controller 115 is connected to display device 119, and controls the demonstration in display device 119.Data reader/write device 116 is reconciled the data communication between CPU111 and recording medium, and carry out by program, reading and carrying out from recording medium 120 will be written to recording medium 120 by the handled result of computing machine 110.Communication interface 117 is reconciled the data communication between CPU111 and another computing machine.
Such as CF(Compact Flash(registered trademark)) or SD(Secure Digital) general semiconductor memory storage, such as the magnetic storage medium of flexible disk or such as CD-ROM(Compact Disk Read Only Memory) optical storage medium be presented the specific example as recording medium 120.
Above-described embodiment can fully or partly be represented by below described note 1 to 12, but be not limited to following description.
(complementary annotations 1)
A kind of seismic intensity estimation unit, it comprises:
Focus determining unit, when the seismograph of installing in each region in multiple regions of setting on border along particular continent plate detects that while being used to specify the data of described earthquake wave amplitude by the caused seismic event of generation of earthquake and output, described focus determining unit determines whether the position of the focus of the earthquake occurring is in preset areas;
Adjacent area designating unit, if determine that by described focus determining unit the position of described focus is in described district, described adjacent area designating unit is specified the adjacent adjacent area in region being installed in the described seismograph that described seismic event first detected;
Decay determining unit, described decay determining unit is obtained the data of output thus from being arranged on described seismograph described adjacent area, and determines in the specified amplitude of the data by accessed whether decay based on pre-conditioned;
Fault Rupture site setting unit, if determined and not yet decayed by described decay determining unit, the place of wherein estimating to occur Fault Rupture is set in the position of described Fault Rupture site setting unit based on described adjacent area in described district; And
Estimate seismic intensity computing unit, described estimation seismic intensity computing unit calculates estimation seismic intensity centered by the place setting by described Fault Rupture site setting unit.
(complementary annotations 2)
According in the seismic intensity estimation unit of complementary annotations 1, set described district with the band shape on the border along described continental-bearing plate.
(complementary annotations 3)
According in the seismic intensity estimation unit of complementary annotations 2,
Described adjacent area designating unit is specified new adjacent area based on described adjacent area,
Described decay determining unit is obtained the data of output thus from being arranged on described seismograph described new adjacent area, and pre-conditionedly determines in the specified amplitude of the data by accessed whether decay based on described;
If determine and not yet decay, the position of described Fault Rupture site setting unit based on described new adjacent area resets the place of wherein estimating to occur Fault Rupture on the border of described continental-bearing plate, and
Described estimation seismic intensity computing unit further calculates described estimation seismic intensity centered by the place resetting by described Fault Rupture site setting unit.
(complementary annotations 4)
According in the seismic intensity estimation unit of any one in complementary annotations 1 to 3,
In each region in described multiple regions, multiple seismographs are installed, and
Described decay determining unit is obtained described data from the each described seismograph being arranged on the adjacent area of appointment, obtain by the mean value of the specified described amplitude of the data of each acquisition, and use obtained mean value to determine whether to decay.
(complementary annotations 5)
According in the seismic intensity estimation unit of any one in complementary annotations 1 to 4,
If be more than or equal to setting threshold so that continued setting-up time section by the specified described amplitude of accessed data, described decay determining unit is definite not yet decays.
(complementary annotations 6)
A kind of seismic intensity method of estimation, it comprises following steps:
(a) when the seismograph of installing in each region in the multiple regions that set along the border of particular continent plate detects while being used to specify the data of described earthquake wave amplitude by the caused seismic event of generation of earthquake and output, determine whether the position of the focus of the earthquake occurring is in preset areas;
(b), if determine that in step (a) position of described focus is in described district, specify the adjacent adjacent area in region being installed in the described seismograph that described seismic event first detected;
(c) obtain the data of output thus from being arranged on described seismograph described adjacent area, and determine in the specified described amplitude of the data by accessed whether decay based on pre-conditioned;
(d), if determine and not yet decay in step (c), the place of wherein estimating to occur Fault Rupture is set in the position based on described adjacent area in described district; And
(e), centered by the place being set in step (d), estimation seismic intensity is calculated.
(complementary annotations 7)
According in the seismic intensity method of estimation of complementary annotations 5, set described district with the band shape on the border along described continental-bearing plate.
(complementary annotations 8)
According in the seismic intensity method of estimation of complementary annotations 7,
In step (b), specify new adjacent area based on described adjacent area,
In step (c), obtain from the described seismograph being arranged on described new adjacent area the data of being exported by this seismograph, and pre-conditionedly determine in the specified amplitude of the data by accessed whether decay based on described;
In step (d), if determine and not yet decay in step (c), the position of the adjacent area based on new resets the place of wherein estimating to occur Fault Rupture on the border of continental-bearing plate, and
In step (e), centered by the place resetting in step (d), estimation seismic intensity is further calculated.
(complementary annotations 9)
According in the seismic intensity method of estimation of any one in complementary annotations 6 to 8,
In each in described multiple regions, multiple seismographs are installed, and
In step (c), obtain data from the each described seismograph being arranged on specified adjacent area, obtain by the mean value of the specified amplitude of the data of each acquisition, and use obtained mean value to determine whether to decay.
(complementary annotations 10)
According in the seismic intensity method of estimation of any one in complementary annotations 6 to 9,
In step (c), if be more than or equal to setting threshold so that continue setting-up time section by the specified amplitude of accessed data, determine that decay not yet occurs.
(complementary annotations 11)
The computer readable recording medium storing program for performing having program stored therein, described routine package contains for making computing machine carry out the order of following steps:
(a) when the seismograph of installing in each region in the multiple regions that set along the border of particular continent plate detects while being used to specify the data of described earthquake wave amplitude by the caused seismic event of generation of earthquake and output, determine whether the position of the focus of the earthquake occurring is in preset areas;
(b), if determine that in step (a) position of described focus is in described district, specify the adjacent adjacent area in region being installed in the described seismograph that described seismic event first detected;
(c) obtain the data of output thus from being arranged on described seismograph described adjacent area, and determine in the specified described amplitude of the data by accessed whether decay based on pre-conditioned;
(d), if determine and not yet decay in step (c), the place of wherein estimating to occur Fault Rupture is set in the position based on described adjacent area in described district; And
(e), centered by the place being set in step (d), estimation seismic intensity is calculated.
(complementary annotations 12)
According in the computer readable recording medium storing program for performing of complementary annotations 11, set described district with the band shape on the border along described continental-bearing plate.
(complementary annotations 13)
According in the computer readable recording medium storing program for performing of complementary annotations 12,
In step (b), specify new adjacent area based on described adjacent area,
In step (c), obtain from the described seismograph being arranged on described new adjacent area the data of being exported by this seismograph, and pre-conditionedly determine in the specified amplitude of the data by accessed whether decay based on described;
In step (d), if determine and not yet decay in step (c), the position of the adjacent area based on new resets the place of wherein estimating to occur Fault Rupture on the border of continental-bearing plate, and
In step (e), centered by the place resetting in step (d), estimation seismic intensity is further calculated.
(complementary annotations 14)
According in the computer readable recording medium storing program for performing of any one in complementary annotations 11 to 13,
In each in described multiple regions, multiple seismographs are installed, and
In step (c), obtain data from the each described seismograph being arranged on specified adjacent area, obtain by the mean value of the specified amplitude of the data of each acquisition, and use obtained mean value to determine whether to decay.
(complementary annotations 15)
According in the computer readable recording medium storing program for performing of any one in complementary annotations 11 to 14,
In step (c), if be more than or equal to setting threshold so that continue setting-up time section by the specified amplitude of accessed data, determine that decay not yet occurs.
Although reference example has been described invention required for protection in the above, invention required for protection is not limited to above-described embodiment.It will be appreciated by those of skill in the art that, in scope of invention required for protection, can make various modifications to the configuration of invention required for protection and details.
The application based on and require the rights and interests of right of priority of the existing Japanese patent application No.2011-208608 submitting on September 26th, 2011, its whole disclosure is combined in herein by reference.
Industrial usability
According to the present invention, even if in the situation that subduction zone earthquake therein occurs, also can specify reliably the predicted region of arriving of the generation of wherein trembling, and can estimate the seismic intensity in this region.The present invention is applicable to need to issue as quickly as possible the system about the information of earthquake in the time that earthquake occurs.
Description of reference numerals
10 seismic intensity estimation units
11 focus determining units
12 adjacent area designating unit
13 decay determining units
14 Fault Rupture site setting unit
15 estimate seismic intensity computing unit
16 communication units
20 seismographs
21a to 21d seismograph group
The border of 22 continental-bearing plates
23 focus
24 Fault Rupture places
30 designation area
40 seismism recording geometrys
50 networks
110 computing machines
111?CPU
112 primary memorys
113 memory storages
114 input interfaces
115 display controllers
116 data readers/write device
117 communication interfaces
118 input medias
119 display device
120 recording mediums
121 buses
Claims (7)
1. a seismic intensity estimation unit, it comprises:
Focus determining unit, when the seismograph of installing in each region in multiple regions of setting on border along particular continent plate detects that while being used to specify the data of described earthquake wave amplitude by the caused seismic event of generation of earthquake and output, described focus determining unit determines whether the position of the focus of the earthquake occurring is in preset areas;
Adjacent area designating unit, if determine that by described focus determining unit the position of described focus is in described district, described adjacent area designating unit is specified the adjacent adjacent area in region being installed in the described seismograph that described seismic event first detected;
Decay determining unit, described decay determining unit is obtained the data of output thus from being arranged on described seismograph described adjacent area, and determines in the specified amplitude of the data by accessed whether decay based on pre-conditioned;
Fault Rupture site setting unit, if determined and not yet decayed by described decay determining unit, the place of wherein estimating to occur Fault Rupture is set in the position of described Fault Rupture site setting unit based on described adjacent area in described district; And
Estimate seismic intensity computing unit, described estimation seismic intensity computing unit, centered by the place setting by described Fault Rupture site setting unit, calculates estimation seismic intensity.
2. seismic intensity estimation unit according to claim 1, wherein,
Band shape with the border along described continental-bearing plate is set described district.
3. seismic intensity estimation unit according to claim 2, wherein,
Described adjacent area designating unit is specified new adjacent area based on described adjacent area,
Described decay determining unit is obtained the data of output thus from being arranged on described seismograph described new adjacent area, and pre-conditionedly determines in the specified amplitude of the data by accessed whether decay based on described;
If determine and not yet decay, the position of described Fault Rupture site setting unit based on described new adjacent area resets the place of wherein estimating to occur Fault Rupture on the border of described continental-bearing plate, and
Described estimation seismic intensity computing unit, centered by the place resetting by described Fault Rupture site setting unit, further calculates described estimation seismic intensity.
4. according to the seismic intensity estimation unit described in any one in claims 1 to 3, wherein,
In each region in described multiple regions, multiple seismographs are installed, and
Described decay determining unit is obtained described data from the each described seismograph of installing the adjacent area of specifying, obtain the mean value by the specified described amplitude of each obtained data, and use obtained mean value to determine whether to decay.
5. according to the seismic intensity estimation unit described in any one in claim 1 to 4,
If be more than or equal to setting threshold so that continued setting-up time section by the specified described amplitude of accessed data, described decay determining unit is definite not yet decays.
6. a seismic intensity method of estimation, it comprises following steps:
(a) when the seismograph of installing in each region in the multiple regions that set along the border of particular continent plate detects while being used to specify the data of described earthquake wave amplitude by the caused seismic event of generation of earthquake and output, determine whether the position of the focus of the earthquake occurring is in preset areas;
(b), if determine that in step (a) position of described focus is in described district, specify the adjacent adjacent area in region being installed in the described seismograph that described seismic event first detected;
(c) obtain the data of output thus from being arranged on described seismograph described adjacent area, and determine in the specified described amplitude of the data by accessed whether decay based on pre-conditioned;
(d), if determine and not yet decay in step (c), the place of wherein estimating to occur Fault Rupture is set in the position based on described adjacent area in described district; And
(e), centered by the place being set in step (d), estimation seismic intensity is calculated.
7. the computer readable recording medium storing program for performing having program stored therein, described routine package contains for making computing machine carry out the order of following steps:
(a) when the seismograph of installing in each region in the multiple regions that set along the border of particular continent plate detects while being used to specify the data of described earthquake wave amplitude by the caused seismic event of generation of earthquake and output, determine whether the position of the focus of the earthquake occurring is in preset areas;
(b), if determine that in step (a) position of described focus is in described district, specify the adjacent adjacent area in region being installed in the described seismograph that described seismic event first detected;
(c) obtain the data of output thus from being arranged on described seismograph described adjacent area, and determine in the specified described amplitude of the data by accessed whether decay based on pre-conditioned;
(d), if determine and not yet decay in step (c), the place of wherein estimating to occur Fault Rupture is set in the position based on described adjacent area in described district; And
(e), centered by the place being set in step (d), estimation seismic intensity is calculated.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2011208608 | 2011-09-26 | ||
| JP2011-208608 | 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 true CN103827694A (en) | 2014-05-28 |
| CN103827694B 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 (6)
| 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 |
| US20090033511A1 (en) * | 2007-07-31 | 2009-02-05 | Nec Corporation | Earthquake noticing server, earthquake noticing system, earthquake noticing method, and program |
| CN102156298A (en) * | 2011-03-11 | 2011-08-17 | 中国地震局地壳应力研究所 | Rapid seismic intensity assessment method based on non-statistical hypothesis test |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 CN CN201280046946.0A patent/CN103827694B/en not_active Expired - Fee Related
- 2012-08-21 JP JP2013536076A patent/JP5742951B2/en active Active
- 2012-08-21 WO PCT/JP2012/071065 patent/WO2013047008A1/en active Application Filing
Patent Citations (6)
| 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 |
| US20090033511A1 (en) * | 2007-07-31 | 2009-02-05 | Nec Corporation | Earthquake noticing server, earthquake noticing system, earthquake noticing method, and program |
| CN102156298A (en) * | 2011-03-11 | 2011-08-17 | 中国地震局地壳应力研究所 | Rapid seismic intensity assessment method based on non-statistical hypothesis test |
Non-Patent Citations (2)
| Title |
|---|
| 刘瑞丰 等: "中国地震台网震级的对比", 《地震学报》 * |
| 张杨 等: "华北地区烈度衰减模型建立及其用于震中区域和震级的定量估算", 《地震学报》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103827694B (en) | 2016-11-02 |
| JP5742951B2 (en) | 2015-07-01 |
| WO2013047008A1 (en) | 2013-04-04 |
| JPWO2013047008A1 (en) | 2015-03-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| WO2021033547A1 (en) | Seismic observation device, seismic observation method, and recording medium | |
| Meruane | Online sequential extreme learning machine for vibration-based damage assessment using transmissibility data | |
| CN110780347B (en) | Earthquake destructive power prediction device and method based on cyclic neural network | |
| US7392139B2 (en) | System and method for providing earthquake data | |
| JP6915279B2 (en) | Worker allocation device, worker allocation method and worker allocation program | |
| Lin et al. | Development of a tsunami early warning system for the South China Sea | |
| CN113076358A (en) | Report generation method, device, equipment and storage medium | |
| US20170186145A1 (en) | Method for determining reusability, apparatus for displaying boundary sample, and method for displaying boundary sample | |
| Kuge | Source modeling using strong-motion waveforms: toward automated determination of earthquake fault planes and moment-release distributions | |
| JP6981061B2 (en) | Production plan generation program, production plan generation method and production plan generation device | |
| Roumelioti et al. | Slip distribution of the 7 September 1999 Athens earthquake inferred from an empirical Green's function study | |
| CN103827694A (en) | Seismic intensity estimation device, seismic intensity estimation method, and computer-readable recording medium | |
| JP7406828B2 (en) | Earthquake estimation method, earthquake estimation program, and earthquake estimation device | |
| Rohmer et al. | Source characterisation by mixing long-running tsunami wave numerical simulations and historical observations within a metamodel-aided ABC setting | |
| US20170184488A1 (en) | Facility state analysis device, analysis method for facility state, storage medium, and facility management system | |
| JP5770051B2 (en) | Immediate estimation method of the rupture area of a large earthquake | |
| JP5157779B2 (en) | Earthquake early warning device, earthquake early warning method, and earthquake early warning program | |
| CN110545145B (en) | Method for performing radio wave communication in bimodal terrain through computer simulation | |
| TWI647475B (en) | An earthquake prediction method and system | |
| JP5760268B2 (en) | Earthquake motion duration prediction system | |
| CN103778057B (en) | LBS method of testings and device | |
| CN103782197A (en) | Magnitude estimation device, magnitude estimation method, and computer-readable recording medium | |
| CN110635856B (en) | Method for performing radio wave communication in forest land terrain through computer simulation | |
| Carvalho | The PoD curve for the detection of planar defects using a multi-channel ultrasonic system | |
| US10234374B2 (en) | Method for selecting boundary sample, apparatus for selecting boundary sample, storage medium, and system for selecting boundary sample |
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 |