CN108229028A - A kind of method for weighting Seismic net than generation based on ceiling capacity - Google Patents

A kind of method for weighting Seismic net than generation based on ceiling capacity Download PDF

Info

Publication number
CN108229028A
CN108229028A CN201810010153.2A CN201810010153A CN108229028A CN 108229028 A CN108229028 A CN 108229028A CN 201810010153 A CN201810010153 A CN 201810010153A CN 108229028 A CN108229028 A CN 108229028A
Authority
CN
China
Prior art keywords
node
earthquake
network
magnitude
represent
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201810010153.2A
Other languages
Chinese (zh)
Other versions
CN108229028B (en
Inventor
王超飞
任涛
李松威
王冉冉
张轶馨
李伟超
刘思邈
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northeastern University China
Original Assignee
Northeastern University China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Northeastern University China filed Critical Northeastern University China
Priority to CN201810010153.2A priority Critical patent/CN108229028B/en
Publication of CN108229028A publication Critical patent/CN108229028A/en
Application granted granted Critical
Publication of CN108229028B publication Critical patent/CN108229028B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • 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

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention belongs to complex network technical fields, and in particular to a kind of based on method of the ceiling capacity than generation weighting Seismic net.Technical solution is as follows：A kind of method for being weighted Seismic net than generation based on ceiling capacity, is included the following steps：Step 1：Build network analysis topological property；Step 2：According to magnitude, influence diameter and influence time relationship carries out network modelling；Step 3：Undirected Seismic net and oriented Seismic net are generated by simulation；Step 4：Define based on ceiling capacity than company's side right value method；Step 5：Draw out the distribution probability table of side right value in network；Step 6：Draw the probability distribution graph of side right value in network；Step 7：Draw the probability double-log distribution map of side right value in network.The method provided by the invention for weighting Seismic net than generation based on ceiling capacity, for Study of Seismic propagation path.

Description

A kind of method for weighting Seismic net than generation based on ceiling capacity
Technical field
The invention belongs to complex network technical fields, and in particular to one kind is based on ceiling capacity than generation weighting Seismic net Method.
Background technology
Through the ages, earthquake often brings serious harm to the mankind, and the annual whole world can all have about ten to be caused several times to the mankind The earthquake of grievous injury occurs.Earthquake is difficult to avoid that, but if can correctly understand its rule, scientifically grasps earthquake Affecting laws and countermeasure can effectively mitigate earthquake and loss.It can study ground by building Seismic net model Propagation path is shaken, existing weighting Seismic net is using common earthquake number or earthquake effect number as side right value.
Invention content
The present invention provides a kind of method for weighting Seismic net than generation based on ceiling capacity, and road is propagated for Study of Seismic Diameter.
Technical scheme is as follows：
A kind of method for being weighted Seismic net than generation based on ceiling capacity, is included the following steps：
Step 1：Network analysis topological property is built, network node is divided using the method for longitude and latitude partitioning site, by ground The geographic area that is included of shake data is divided into several sub-boxes, remove geographic area that incoherent seismic data includes and After earthquake magnitude is less than 1 grade of zonule, these regions are numbered；
Step 2：According to magnitude, influence diameter and influence time relationship carries out network modelling；
1) for network node i, earthquake magnitude is with influencing diameter and influencing the relational expression of duration：
log Ti=a1Mi+b1 (1)
logLi=a2Mi+b2 (2)
Wherein, MiExpression is happened at the magnitude of node i, TiRepresent that M occurs for node iiInfluence time after grade earthquake, LiRepresent that M occurs for node iiBiggest impact distance after earthquake, that is, influence diameter, a1、a2、b1、b2It is according to the normal of statistics gained Number；
2) judged between node i and node j with the presence or absence of even side by formula (3) and (4)：
Δ t=ti-tj (3)
Wherein, Δ t represents node earthquake interval, tiRepresent the earthquake moment of node i, tjWhen representing the earthquake of node j It carves, dijRepresent the spacing distance between node i and node j, R is earth radius, ilatRepresent that the longitude of earthquake occurs for node i, ilonRepresent that the latitude of earthquake, j occur for node ilatRepresent that the longitude of earthquake, j occur for node jlonRepresent that earthquake occurs for node j Latitude；The earthquake magnitude of acquiescence focus is more than its region earthquake magnitude influenced；
3) the company side between nodes is defined：In one earthquake event, if node i, j successively occur earthquake and Node j is within the coverage of node i and influence time, i.e. Δ t≤TiAnd dij≤Li, then exist between node i and j and connect Side；If earthquake has occurred in node i, j, v successively and node v is within the influence time and coverage of node i and j, There is even side between node v and node i and j；Weight side is excluded, while excludes the earthquake that node is occurred by itself affect, i.e., from ring；
Step 3：Undirected Seismic net and oriented Seismic net are generated by simulation；
Step 4：Define based on ceiling capacity than company's side right value method：All it is that primary energy discharges per secondary earthquake Process, and the relationship between energy and earthquake magnitude is as follows：
E=101.5*M+11.8 (5)
Wherein, E is the energy of earthquake, and M is earthquake magnitude；
Define the side right value s of m secondary earthquake event interior joints i and jijmFor：
Wherein, Eim, EjmThe energy of node i and j during m secondary earthquakes are represented respectively；
Formula (5) is substituted into formula (6), i.e.,
Wherein, MimRepresent the earthquake magnitude of focus node i during m secondary earthquakes, MjmRepresent the earthquake magnitude of node j during m secondary earthquakes；
Step 5：Draw out the distribution probability table of side right value in network；
Step 6：Draw the probability distribution graph of side right value in network；
Step 7：Draw the probability double-log distribution map of side right value in network.
The method for weighting Seismic net than generation based on ceiling capacity, wherein software environment are Windows7 systems System, simulated environment Matlab.
Beneficial effects of the present invention are：The present invention provides a kind of new method to define the side right value in Seismic net model, The method provides more efficient information for structure Seismic net model and prediction earthquake propagation path.
Description of the drawings
Fig. 1 is the oriented earthquake network of 0.1*0.1 longitudes and latitudes partitioning site generation；
Fig. 2 is based on undirected Seismic net side right value probability distribution graph of the ceiling capacity than generation；
Fig. 3 is based on oriented earthquake network edge weights probability distribution graph of the ceiling capacity than generation；
Fig. 4 is based on undirected Seismic net side right value probability double-log distribution map of the ceiling capacity than generation；
Fig. 5 is based on oriented earthquake network edge weights probability double-log distribution map of the ceiling capacity than generation.
Specific embodiment
A kind of method for weighting Seismic net than generation based on ceiling capacity, software environment are Windows7 systems, are emulated Environment is Matlab, is included the following steps：
Step 1：Network analysis topological property is built, network node is divided using the method for longitude and latitude partitioning site, by ground The geographic area that is included of shake data is divided into several sub-boxes, remove geographic area that incoherent seismic data includes and After earthquake magnitude is less than 1 grade of zonule, these regions are numbered；It is 579 to choose number of nodes, and number of edges is 3496；
Step 2：According to magnitude, influence diameter and influence time relationship carries out network modelling；
1) for network node i, earthquake magnitude is with influencing diameter and influencing the relational expression of duration：
log Ti=a1Mi+b1 (1)
logLi=a2Mi+b2 (2)
Wherein, MiExpression is happened at the magnitude of node i, TiRepresent that M occurs for node iiInfluence time after grade earthquake, LiRepresent that M occurs for node iiBiggest impact distance after earthquake, that is, influence diameter, a1、a2、b1、b2It is according to the normal of statistics gained Number；
2) judged between node i and node j with the presence or absence of even side by formula (3) and (4)：
Δ t=ti-tj (3)
Wherein, Δ t represents node earthquake interval, tiRepresent the earthquake moment of node i, tjWhen representing the earthquake of node j It carves, dijRepresent the spacing distance between node i and node j, R is earth radius, ilatRepresent that the longitude of earthquake occurs for node i, ilonRepresent that the latitude of earthquake, j occur for node ilatRepresent that the longitude of earthquake, j occur for node jlonRepresent that earthquake occurs for node j Latitude；The earthquake magnitude of acquiescence focus is more than its region earthquake magnitude influenced；
3) the company side between nodes is defined：In one earthquake event, if node i, j successively occur earthquake and Node j is within the coverage of node i and influence time, i.e. Δ t≤TiAnd dij≤Li, then exist between node i and j and connect Side；If earthquake has occurred in node i, j, v successively and node v is within the influence time and coverage of node i and j, There is even side between node v and node i and j；Weight side is excluded, while excludes the earthquake that node is occurred by itself affect, i.e., from ring；
Step 3：It is simulated by Matlab and generates undirected Seismic net and oriented Seismic net；Oriented Seismic net such as Fig. 1 It is shown；
Step 4：Define based on ceiling capacity than company's side right value method：All it is that primary energy discharges per secondary earthquake Process, and the relationship between energy and earthquake magnitude is as follows：
E=101.5*M+11.8 (5)
Wherein, E is the energy of earthquake, and M is earthquake magnitude；Earthquake mainly passes through seimic wave propagation energy；Seismic wave is divided into shear wave And longitudinal wave, longitudinal wave causes ground dipping and heaving, and shear wave carries main energy and transmits energy along seismic spread direction, so as to Cause close region that earthquake occurs；
Define the side right value s of m secondary earthquake event interior joints i and jijmFor：
Wherein, Eim, EjmThe energy of node i and j during m secondary earthquakes are represented respectively；
Formula (5) is substituted into formula (6), i.e.,
Wherein, MimRepresent the earthquake magnitude of focus node i during m secondary earthquakes, MjmRepresent the earthquake magnitude of node j during m secondary earthquakes；
Step 5：Draw out the distribution probability table of side right value in network；It is shown in Table 1：
Table 1 based on ceiling capacity than Seismic net side right value probability
As shown in Table 1, in undirected Seismic net, there is not the company side that weights are more than or equal to 1, and weights are not less than 0.5 Lian Bianzhan overall 18%；In oriented Seismic net, similary there is no the company side that weights are more than or equal to 1, weights are not less than 0.5 company side accounts for 21%；Because when earthquake occurs, close region spot can be influenced by only carrying the shear wave of most of energy Shake, therefore the company side that weights are more than or equal to 1 can not possibly occur；
Step 6：The probability distribution graph of side right value in network is drawn, as shown in Figure 2,3；
Step 7：Draw the probability double-log distribution map of side right value in network；It is taken by opposite side weights and its shared probability Double-log, and then study whether side right value meets power-law distribution characteristic；As shown in Figures 4 and 5, either in undirected Seismic net also It is in oriented Seismic net, side right value all disobeys power-law distribution；
In conclusion the weighting earthquake network model based on ceiling capacity than generation has prediction earthquake propagation path More efficient meaning to predict the generation of earthquake from now on, reduces earthquake and provides strong foundation to the harm that the mankind bring.

Claims (2)

  1. It is 1. a kind of based on method of the ceiling capacity than generation weighting Seismic net, which is characterized in that include the following steps：
    Step 1：Network analysis topological property is built, network node is divided using the method for longitude and latitude partitioning site, by earthquake number Several sub-boxes are divided into according to the geographic area included, remove geographic area and earthquake magnitude that incoherent seismic data includes Behind less than 1 grade of zonule, these regions are numbered；
    Step 2：According to magnitude, influence diameter and influence time relationship carries out network modelling；
    1) for network node i, earthquake magnitude is with influencing diameter and influencing the relational expression of duration：
    log Ti=a1Mi+b1 (1)
    log Li=a2Mi+b2 (2)
    Wherein, MiExpression is happened at the magnitude of node i, TiRepresent that M occurs for node iiInfluence time after grade earthquake, LiIt represents M occurs for node iiBiggest impact distance after earthquake, that is, influence diameter, a1、a2、b1、b2It is the constant according to statistics gained；
    2) judged between node i and node j with the presence or absence of even side by formula (3) and (4)：
    Δ t=ti-tj (3)
    Wherein, Δ t represents node earthquake interval, tiRepresent the earthquake moment of node i, tjRepresent the earthquake moment of node j, dijRepresent the spacing distance between node i and node j, R is earth radius, ilatRepresent that the longitude of earthquake, i occur for node ilonTable Show that the latitude of earthquake, j occur for node ilatRepresent that the longitude of earthquake, j occur for node jlonRepresent that the latitude of earthquake occurs for node j； The earthquake magnitude of acquiescence focus is more than its region earthquake magnitude influenced；
    3) the company side between nodes is defined：In one earthquake event, if earthquake and node successively occur for node i, j J is within the coverage of node i and influence time, i.e. Δ t≤TiAnd dij≤Li, then there is even side between node i and j；Such as Earthquake has occurred in fruit node i, j, v successively and node v is within the influence time and coverage of node i and j, then node v There is even side between node i and j；Weight side is excluded, while excludes the earthquake that node is occurred by itself affect, i.e., from ring；
    Step 3：Undirected Seismic net and oriented Seismic net are generated by simulation；
    Step 4：Define based on ceiling capacity than company's side right value method：All it is the process of primary energy release per secondary earthquake, And the relationship between energy and earthquake magnitude is as follows：
    E=101.5*M+11.8 (5)
    Wherein, E is the energy of earthquake, and M is earthquake magnitude；
    Define the side right value s of m secondary earthquake event interior joints i and jijmFor：
    Wherein, Eim, EjmThe energy of node i and j during m secondary earthquakes are represented respectively；
    Formula (5) is substituted into formula (6), i.e.,
    Wherein, MimRepresent the earthquake magnitude of focus node i during m secondary earthquakes, MjmRepresent the earthquake magnitude of node j during m secondary earthquakes；
    Step 5：Draw out the distribution probability table of side right value in network；
    Step 6：Draw the probability distribution graph of side right value in network；
    Step 7：Draw the probability double-log distribution map of side right value in network.
  2. It is 2. according to claim 1 based on method of the ceiling capacity than generation weighting Seismic net, which is characterized in that software Environment be Windows7 systems, simulated environment Matlab.
CN201810010153.2A 2018-01-05 2018-01-05 Method for generating weighted seismic network based on maximum energy ratio Expired - Fee Related CN108229028B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810010153.2A CN108229028B (en) 2018-01-05 2018-01-05 Method for generating weighted seismic network based on maximum energy ratio

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810010153.2A CN108229028B (en) 2018-01-05 2018-01-05 Method for generating weighted seismic network based on maximum energy ratio

Publications (2)

Publication Number Publication Date
CN108229028A true CN108229028A (en) 2018-06-29
CN108229028B CN108229028B (en) 2021-03-02

Family

ID=62642947

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810010153.2A Expired - Fee Related CN108229028B (en) 2018-01-05 2018-01-05 Method for generating weighted seismic network based on maximum energy ratio

Country Status (1)

Country Link
CN (1) CN108229028B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106468782A (en) * 2015-08-19 2017-03-01 中国石油化工股份有限公司 A kind of based on ceiling capacity than method crack prediction method
CN106533742A (en) * 2016-10-31 2017-03-22 天津大学 Time sequence mode representation-based weighted directed complicated network construction method
CN107220493A (en) * 2017-05-24 2017-09-29 王欣 Shale gas horizontal well network fracture modeling method based on micro-seismic event

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106468782A (en) * 2015-08-19 2017-03-01 中国石油化工股份有限公司 A kind of based on ceiling capacity than method crack prediction method
CN106533742A (en) * 2016-10-31 2017-03-22 天津大学 Time sequence mode representation-based weighted directed complicated network construction method
CN107220493A (en) * 2017-05-24 2017-09-29 王欣 Shale gas horizontal well network fracture modeling method based on micro-seismic event

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
SUMIYOSHI ABEA 等: "Small-worldstructure of earthquake network", 《PHYSICA A》 *
何璇 等: "基于时空影响域的地震网络构造方法", 《东北大学学报（自然科学版）》 *
徐久强 等: "基于时空影响域地震加权网络的研究", 《东北大学学报（自然科学版）》 *
赵海 等: "基于时空影响域的地震网络动力学演化特征分析", 《东北大学学报（自然科学版）》 *

Also Published As

Publication number Publication date
CN108229028B (en) 2021-03-02

Similar Documents

Publication Publication Date Title
Okal The quest for wisdom: lessons from 17 tsunamis, 2004–2014
Megawati et al. Response spectral attenuation relationships for Sumatran-subduction earthquakes and the seismic hazard implications to Singapore and Kuala Lumpur
CN109375252A (en) Consider the ground motion parameter evaluation method of different seismic structure maximum credible earthquakes
CN109375253A (en) Ground motion parameter evaluation method based on whole seismic structure maximum credible earthquakes
CN116203619B (en) Regional earthquake landslide simulation method based on spatial cross-correlation multi-seismic vibration parameters
Padgham et al. Planning for sandbagging as a response to flooding: A tool and case study
Tiwari Status of seismicity in the northeast India and earthquake disaster mitigation
Sato Fractal characteristics of phase spectrum of earthquake motion
CN117270047B (en) Earthquake risk evaluation method, system, equipment and medium
Karimzadeh et al. A study on fragility analyses of masonry buildings in Erzincan (Turkey) utilizing simulated and real ground motion records
CN108229028A (en) A kind of method for weighting Seismic net than generation based on ceiling capacity
Kemal et al. Evaluation of tsunami evacuation plans in Padang, Indonesia
Ladina et al. Feasibility study of an earthquake early warning system in eastern central Italy
Orsolya et al. Seismic risk and disaster management perspectives in Hungary, presented on a case study performed in Győr
Sorabella et al. Ground motion suite selection for Boston, Massachusetts
Sofyan et al. Estimation of return period of the destructive earthquake in aceh using the maximum likelihood method
Aulov et al. Seismic hazard analysis (SHA) considerations for South China Sea territories of Malaysia
Lemoine et al. Pligurian earthquake: Seismic and tsunami scenario modeling, from hazard to risk assessment towards evacuations planning
Fujiwara et al. Revision of probabilistic seismic hazard assessment for Japan after the 2011 Tohoku-oki Mega-thrust earthquake (M9. 0)
CN108229743B (en) Method for constructing weighting network based on maximum magnitude ratio
Fujiwara et al. Seismic hazard assessment for Japan after the 2011 Tohoku-oki megathrust earthquake (Mw9. 0)
Torrentira Jr et al. Framework of disaster preparedness among coastal communities in Davao City, Philippines
Termizi et al. Earthquake Threats in Ranau–From The Sources of Mensaban and Mesilou Fault
Bevilacqua Doubly stochastic models for volcanic vent opening probability and pyroclastic density current hazard at Campi Flegrei caldera
Igarashi et al. Building a tsunami simulation database for the tsunami warning system in the Philippines

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210302