CN105510962A - Method for simulation of scope of aftershock - Google Patents
Method for simulation of scope of aftershock Download PDFInfo
- Publication number
- CN105510962A CN105510962A CN201510432096.3A CN201510432096A CN105510962A CN 105510962 A CN105510962 A CN 105510962A CN 201510432096 A CN201510432096 A CN 201510432096A CN 105510962 A CN105510962 A CN 105510962A
- Authority
- CN
- China
- Prior art keywords
- aftershock
- axis
- point
- coordinate system
- coverage
- 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.)
- Pending
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a method for simulation of the scope of aftershocks, and the method comprises the steps of establishing a coordinate system with the aftershock point longitude as an X axis, the aftershock point latitude as a Y axis and the aftershock point magnitude as a Z axis; taking a plane of the X axis and the Y axis as a projection plane, and transforming the coordinate system into a projection coordinate system; obtaining a relationship mode of the distance among aftershock points through fitting by using a correlation analysis method of point pattern analysis, decomposing the relationship mode into different scale-free intervals in which a numerical value is less than a numerical value of the next adjacent scale-free interval; performing nuclear density simulation analysis of points in the scale-free intervals, and obtaining corresponding pattern layers to the scale-free intervals respectively; determining a shielding content between adjacent pattern layers Li and Li+1, removing the corresponding shielding content of the Li+1, then combining the Li and the Li+1, processing all the pattern layers in sequence, and outputting a pattern layer Ri after the combination. The invention can quickly evaluate the scope of the aftershocks based on point data of aftershocks released publicly.
Description
Technical field
The present invention relates to seismic surveying technical field, particularly a kind of analogy method of aftershock coverage.
Background technology
Earthquake effect scope simulation is one of element task of earthquake disaster comprehensive assessment.Earthquake disaster is by force sudden, and predictability is poor, and frequency is higher in recent years, and the serious secondary disaster of association is very large to social influence.In this context, study based on open, retrievable Aftershocks point sequence data, carry out earthquake effect scope Fast simulation, its realistic price is apparent.
The existing method to the assessment of earthquake disaster coverage, be main mainly with parameters such as Maximum earthquake magnitude, ground peak accelerators, implement the assessment of earthquake disaster coverage subregion, this method needs more seismic physical parameter, therefore process is loaded down with trivial details, and operability is not strong.Based on the point data of the Aftershocks published, estimate earthquake effect scope fast, Data Source is simple, and cost is lower; Based on the data volume of real-time release, repeatedly can assess and approach, real-time and ease of handling stronger.
Summary of the invention
The present invention is intended at least solve one of problems of the prior art.
One of object that the present invention will reach is the point data based on the Aftershocks published, and can realize fast estimating earthquake effect scope.
For solving the problems of the technologies described above, the invention provides a kind of analogy method of aftershock coverage, its step comprises:
Set up with aftershock point longitude be X-axis, the aftershock point latitude coordinate system that is Z axis for Y-axis and aftershock point earthquake magnitude;
With described X-axis and face, Y-axis place for projecting plane, by described coordinate system transformation in projected coordinate system;
Utilize some pattern analysis association analysis method, matching obtain aftershock dot spacing from relation schema, be decomposed into different dimensionless interzone S
1, S
2..., S
n, wherein n is aftershock point number and n is more than or equal to 2, wherein dimensionless interzone S
iinterior numerical value is less than dimensionless interzone S
i+1interior numerical value, wherein 1≤i≤n-1;
In described dimensionless interzone, carry out cuclear density sunykatuib analysis a little, obtain described dimensionless interzone S respectively
1, S
2..., S
ncorresponding layer L
1, L
2..., L
n;
Judge adjacent described layer L
iwith L
i+1between block content, remove L
i+1corresponding block content after, L
iwith L
i+1combine, process All Layers L successively
1, L
2..., L
n, export the layer R after combination
i.
Further, its step also comprises according to described layer R
iidentify aftershock coverage.
Further, for making the map of earthquake effect scope.
Beneficial effect of the present invention there are provided a kind of analogy method of aftershock coverage, the present invention is based on the point data of the Aftershocks published, can rapid evaluation earthquake effect scope, and Data Source is simple, and cost is lower; Meanwhile, the present invention based on the data volume of real-time release, repeatedly can assess and approaches, real-time and ease of handling stronger.
Accompanying drawing explanation
Figure 1 shows that the process flow diagram of the analogy method of a kind of aftershock coverage of the embodiment of the present invention.
Figure 2 shows that the single yardstick of the embodiment of the present invention is the Density Estimator figure of the aftershock of 320km.
Figure 3 shows that the Density Estimator figure of the aftershock of embodiment of the present invention Multiscale Fusion.
Embodiment
Hereafter will describe the present invention in detail in conjunction with specific embodiments.It should be noted that the combination of technical characteristic or the technical characteristic described in following embodiment should not be considered to isolated, they can mutually be combined thus be reached better technique effect.
As described in Figure 1, the invention provides a kind of analogy method of aftershock coverage, its step comprises:
100: set up with aftershock point longitude be X-axis, the aftershock point latitude coordinate system that is Z axis for Y-axis and aftershock point earthquake magnitude;
200: with described X-axis and face, Y-axis place for projecting plane, by described coordinate system transformation in projected coordinate system;
300: the association analysis method utilizing some pattern analysis, matching obtain aftershock dot spacing from relation schema, be decomposed into different dimensionless interzone S1, S2 ..., Sn, wherein n is aftershock point number and n is more than or equal to 2, numerical value wherein in dimensionless interzone Si is less than the numerical value in dimensionless interzone Si+1, wherein 1≤i≤n-1;
400: in described dimensionless interzone, carry out cuclear density sunykatuib analysis a little, obtain respectively described dimensionless interzone S1, S2 ..., layer L1 that Sn is corresponding, L2 ..., Ln;
500: judge to block content between adjacent described layer Li and Li+1, remove Li+1 corresponding block content after, Li and Li+1 combines, successively process All Layers L1, L2 ..., Ln, export combination after layer Ri.
600: identify aftershock coverage according to described layer Ri.
The analogy method of a kind of aftershock coverage of the present invention can also be applied to the map making earthquake effect scope.
The embodiment of the present invention is regional for objective for implementation with Wenchuan:
1. share center (national Seismological science data shares center) at national Seismological science data, download the aftershock data after Wenchuan earthquake, being entered in EXCEL and arranging, data attribute comprises longitude (X), latitude (Y) and earthquake magnitude (Z); In the ArcMap of ArcGIS, select the AddXYdata under Toos, import the EXCEL data put in order, set up with aftershock point longitude be X-axis, aftershock point latitude is for Y-axis and the aftershock point earthquake magnitude coordinate system that is Z axis and carry out projection and transform, and is finally projected as square grid; Load the vector border that Wenchuan aftershock relates to city, county;
2. in EXCEL, carry out correlation dimension analysis to aftershock point, its analytical formula is: C (r)=(2N
r)/N (N-1), (R<r), wherein R is the distance between certain two pairs of earthquake, and r is the scale being used for weighing earthquake spacing chosen, N
rfor distance R is less than the earthquake logarithm of r, N is the total logarithm of earthquake selected; For selected a series of r values, make double logarithmic curve lg [C (r)]-lg (r), curve is found out that section of near linear, carry out matching D by least square method
2=lg [C (r)]/lg (r); R adopts equidistant increment method, increment is 500m, span is 1km-324km, obtains final fitting result and shows: aftershock o'clock linked character used in 13.5-20km, 30-43km, 66-82km, 225-236km, 317-321.5km, five intervals is remarkable;
3. on the basis of Spatial analysis, choose the intermediate value 16km of dimensionless interzone, 36km, 74km, 230km, 320km are as Density Estimator bandwidth;
4. in the ArcMap of ArcGIS, open KernelDensity instrument, inputpoint selects Wenchuan aftershock point data, and populationfield selects the earthquake magnitude field of Wenchuan aftershock point data, and first searchradius is set to 320km, exports layer 320;
5. do not affect experimental result, under the condition that in order to be effective figure is attractive in appearance, filter out the outermost layer of 320.SpatialAnalyst, RasterCalculator, 320> layer minimum value, obtains calculation;
6. in the ArcMap of ArcGIS, open ExtractbyMask instrument, inputraster selects calculation, inputrasterorfeaturemaskdata to select 320, exports c-320;
7.SpatialAnalyst, RasterCalculator, [c-320] * [320] obtains calculation2, and int arranges, and obtains calculation3, as shown in Figure 2;
8. calculate aftershock coverage when Density Estimator bandwidth is 230km, 74km, 36km, 16km respectively according to above-mentioned steps; The result of minimum bandwidth 16km is divided into two-stage, and remaining bandwidth result is one-level, as shown in Figure 3.
The invention provides a kind of analogy method of aftershock coverage, the present invention is based on the point data of the Aftershocks published, can rapid evaluation earthquake effect scope, Data Source is simple, and cost is lower; Meanwhile, the present invention based on the data volume of real-time release, repeatedly can assess and approaches, real-time and ease of handling stronger.
Although give some embodiments of the present invention, it will be understood by those of skill in the art that without departing from the spirit of the invention herein, can change embodiment herein.Above-described embodiment is exemplary, should using embodiment herein as the restriction of interest field of the present invention.
Claims (3)
1. an analogy method for aftershock coverage, is characterized in that, its step comprises:
Set up with aftershock point longitude be X-axis, the aftershock point latitude coordinate system that is Z axis for Y-axis and aftershock point earthquake magnitude;
With described X-axis and face, Y-axis place for projecting plane, by described coordinate system transformation in projected coordinate system;
Utilize some pattern analysis association analysis method, matching obtain aftershock dot spacing from relation schema, be decomposed into different dimensionless interzone S
1, S
2..., S
n, wherein n is aftershock point number and n is more than or equal to 2, wherein dimensionless interzone S
iinterior numerical value is less than dimensionless interzone S
i+1interior numerical value, wherein 1≤i≤n-1;
In described dimensionless interzone, carry out cuclear density sunykatuib analysis a little, obtain described dimensionless interzone S respectively
1, S
2..., S
ncorresponding layer L
1, L
2..., L
n;
Judge adjacent described layer L
iwith L
i+1between block content, remove L
i+1corresponding block content after, L
iwith L
i+1combine, process All Layers L successively
1, L
2..., L
n, export the layer R after combination
i.
2. the analogy method of aftershock coverage according to claim 1, is characterized in that, its step also comprises according to described layer R
iidentify aftershock coverage.
3. the analogy method of aftershock coverage according to claim 1 and 2, is characterized in that, for making the map of earthquake effect scope.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510432096.3A CN105510962A (en) | 2015-07-22 | 2015-07-22 | Method for simulation of scope of aftershock |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510432096.3A CN105510962A (en) | 2015-07-22 | 2015-07-22 | Method for simulation of scope of aftershock |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN105510962A true CN105510962A (en) | 2016-04-20 |
Family
ID=55719070
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510432096.3A Pending CN105510962A (en) | 2015-07-22 | 2015-07-22 | Method for simulation of scope of aftershock |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105510962A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2018044775A (en) * | 2016-09-12 | 2018-03-22 | 戸田建設株式会社 | Forecasting method for expansion of damages from aftershocks, and forecasting system for the same |
| CN113687422A (en) * | 2021-08-27 | 2021-11-23 | 湖北省地震局(中国地震局地震研究所) | Aftershock sequence deleting method based on fault buffer zone |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5877967A (en) * | 1996-03-28 | 1999-03-02 | The United States Of America As Represented By The Secretary Of The Navy | Site and workspaces layout process employing MDS and a PDI formula in which density is calculated using a unit lattice superposed over circumscribing-convex-hulls |
| US6012036A (en) * | 1996-03-28 | 2000-01-04 | The United States Of America As Represented By The Secretary Of The Navy | Site and workspaces layout process employing MDS a PDI formula in which density is calculated using measured span of circumscribing-convex-hulls |
| CN102565855A (en) * | 2012-01-02 | 2012-07-11 | 吉林大学 | Ground micro-seismic data processing method of oil field fracturing |
| CN104199097A (en) * | 2014-09-17 | 2014-12-10 | 西南石油大学 | Novel quantitative judgment method for potential fractures of extrusion structural system |
-
2015
- 2015-07-22 CN CN201510432096.3A patent/CN105510962A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5877967A (en) * | 1996-03-28 | 1999-03-02 | The United States Of America As Represented By The Secretary Of The Navy | Site and workspaces layout process employing MDS and a PDI formula in which density is calculated using a unit lattice superposed over circumscribing-convex-hulls |
| US6012036A (en) * | 1996-03-28 | 2000-01-04 | The United States Of America As Represented By The Secretary Of The Navy | Site and workspaces layout process employing MDS a PDI formula in which density is calculated using measured span of circumscribing-convex-hulls |
| CN102565855A (en) * | 2012-01-02 | 2012-07-11 | 吉林大学 | Ground micro-seismic data processing method of oil field fracturing |
| CN104199097A (en) * | 2014-09-17 | 2014-12-10 | 西南石油大学 | Novel quantitative judgment method for potential fractures of extrusion structural system |
Non-Patent Citations (1)
| Title |
|---|
| 董丞妍 等: "汶川及芦山地震余震分布的空间尺度效应", 《地震学报》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2018044775A (en) * | 2016-09-12 | 2018-03-22 | 戸田建設株式会社 | Forecasting method for expansion of damages from aftershocks, and forecasting system for the same |
| CN113687422A (en) * | 2021-08-27 | 2021-11-23 | 湖北省地震局(中国地震局地震研究所) | Aftershock sequence deleting method based on fault buffer zone |
| CN113687422B (en) * | 2021-08-27 | 2022-07-22 | 湖北省地震局(中国地震局地震研究所) | Aftershock sequence deleting method based on fault buffer zone |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Schorlemmer et al. | Earthquake likelihood model testing | |
| Zuo | Exploring the effects of cell size in geochemical mapping | |
| CN102147479B (en) | Modelling method of reservoir space physical property parameters | |
| CN102841385B (en) | Local geomagnetic chart constructing method based on multi-fractal Krigin method | |
| CN104778369A (en) | Method and system for decision making and early warning based on ground subsidence monitoring | |
| CN106054214B (en) | A kind of GNSS data processing method for weakening multipath effect | |
| CN104330433B (en) | Method and device for obtaining transverse relaxation time distribution of target reservoir | |
| Wang et al. | A spatial econometric modeling of online social interactions using microblogs | |
| Rosindell et al. | Universal scaling of species‐abundance distributions across multiple scales | |
| Irwansyah et al. | Earthquake hazard zonation using peak ground acceleration (PGA) approach | |
| Causse et al. | Calibrating median and uncertainty estimates for a practical use of empirical Green’s functions technique | |
| CN103823951A (en) | Method for quantifying characteristics of road network and effect of characteristics of road network on land utilization | |
| He et al. | Earthquake networks based on space–time influence domain | |
| CN105510962A (en) | Method for simulation of scope of aftershock | |
| CN106556877A (en) | A kind of earth magnetism Tonghua method and device | |
| CN105389631B (en) | More Sea-Crossing Passage distributed areas traffic accessibility analysis methods | |
| Jiang et al. | Effect of removing the common mode errors on linear regression analysis of noise amplitudes in position time series of a regional GPS network & a case study of GPS stations in Southern California | |
| Telesca et al. | Fisher–Shannon analysis of seismograms of tsunamigenic and non-tsunamigenic earthquakes | |
| CN115186037B (en) | Map display method, map display device, electronic device and computer readable medium | |
| Liu et al. | A new risk probability calculation method for urban ecological risk assessment | |
| Li et al. | Application of simulated annealing particle swarm optimization in underwater acoustic positioning optimization | |
| Mascaró et al. | Uncertainty analysis along the ecological quality status of water bodies: the response of the Posidonia oceanica multivariate index (POMI) in three Mediterranean regions | |
| CN102759747A (en) | Method for building seismic data matching pursuit common frequency body | |
| CN102880753B (en) | Based on the land utilization space characteristic dimension conversion method of fractal dimension | |
| CN109100165B (en) | Bridge operation modal analysis method, device, terminal and computer readable storage medium |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160420 |
|
| RJ01 | Rejection of invention patent application after publication |