CN115201909B - Comprehensive judging method and system for risk coefficient of earthquake activity fault section - Google Patents
Comprehensive judging method and system for risk coefficient of earthquake activity fault section Download PDFInfo
- Publication number
- CN115201909B CN115201909B CN202210799016.8A CN202210799016A CN115201909B CN 115201909 B CN115201909 B CN 115201909B CN 202210799016 A CN202210799016 A CN 202210799016A CN 115201909 B CN115201909 B CN 115201909B
- Authority
- CN
- China
- Prior art keywords
- section
- earthquake
- weight
- fault
- sections
- 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.)
- Active
Links
- 230000000694 effects Effects 0.000 title claims abstract description 176
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000004088 simulation Methods 0.000 claims description 16
- 238000005259 measurement Methods 0.000 claims description 7
- 238000012544 monitoring process Methods 0.000 abstract description 2
- 230000000903 blocking effect Effects 0.000 description 11
- 238000010586 diagram Methods 0.000 description 10
- 238000004364 calculation method Methods 0.000 description 4
- 230000035939 shock Effects 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 2
- 238000011081 inoculation Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 206010027336 Menstruation delayed Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/642—Faults
Abstract
The embodiment of the invention relates to the technical field of seismic monitoring, and particularly discloses a method and a system for comprehensively judging risk coefficients of an earthquake activity fault section. The method comprises the steps of dividing a plurality of seismic activity fault sections into a strong-seismic fracture space section, a fault movement locking section, a middle-small seismic sparse section and a coulomb stress enhancement section; subdividing and setting weights of the strong earthquake breaking empty section, the fault movement locking section, the small and medium earthquake sparse section and the coulomb stress enhancement section respectively; and adding the weights corresponding to the seismic activity fault sections, and comprehensively obtaining the risk coefficients corresponding to the seismic activity fault sections. The method can divide the strong-earthquake fracture space section, the fault motion locking section, the middle-small earthquake sparse section and the coulomb stress enhancement section of the earthquake activity fracture section, subdivide the sections, set different weights, comprehensively add the weights to obtain the corresponding danger coefficients of the different earthquake activity fracture sections, and avoid the false report phenomenon of single-method identification.
Description
Technical Field
The invention belongs to the technical field of seismic monitoring, and particularly relates to a method and a system for comprehensively judging risk coefficients of a seismic active fault section.
Background
Strong earthquake is a serious threat affecting social development and welfare of common people due to strong destructiveness and irreversibility of disasters, and strong earthquake prediction has strong social demands. Of course, it is necessary to give definite three elements of "time, space and intensity" for the seismic complex prediction conclusion, and it is preferable to express the prediction result in a quantitative manner. For earthquake time prediction, aiming at different stages of the strong earthquake inoculation process, the time can be divided into: the time periods of the same earthquake, post earthquake, interearthquake and pre earthquake are divided into the following earthquake periods of the strong earthquake inoculation: early, medium and late.
The method for judging the late period between the earthquakes is summarized as the most representative method: a strong earthquake fracture empty section of earthquake geology, a fault movement locking section measured in the earth, a middle and small earthquake sparse section of earthquake activity, a coulomb stress enhancement significant section of numerical simulation and the like. From the past earthquake examples, the strong earthquake crack empty section, the fault movement locking section and the fault coulomb stress enhancement section are common earthquake location prediction methods, the earthquake examples of the middle and small earthquake sparse sections are relatively fewer, meanwhile, the strong earthquake identification results are all 'necessary insufficient' criteria for determining important dangerous areas of earthquakes, and the phenomenon of false report exists in different degrees through single method identification.
Disclosure of Invention
The embodiment of the invention aims to provide a comprehensive judging method and system for risk coefficients of an earthquake activity fault section, and aims to solve the problems in the background technology.
In order to achieve the above object, the embodiment of the present invention provides the following technical solutions:
a comprehensive judging method for risk coefficients of an earthquake activity fault section specifically comprises the following steps:
obtaining information identification results of earthquake geology, geodetic survey, earthquake activity and numerical simulation of main faults of the boundary zone of the movable block, and dividing a plurality of earthquake activity fault sections into strong earthquake fracture empty sections, fault movement locking sections, small and medium earthquake sparse sections and coulomb stress enhancement sections according to the information identification results;
subdividing the strong earthquake fracture empty section, the fault movement locking section, the middle and small earthquake sparse section and the coulomb stress enhancement section respectively to generate subdivision results, and setting different weights for a plurality of earthquake activity fault sections according to the subdivision results to generate weight setting results;
and adding weights corresponding to each seismic activity fault section according to the subdivision result and the weight setting result, and comprehensively obtaining a plurality of risk coefficients corresponding to the seismic activity fault sections.
As a further limitation of the technical scheme of the embodiment of the invention, the method for acquiring the information identification results of the earthquake geology, geodetic measurement, earthquake activity and numerical simulation of the main faults of the boundary zone of the movable block, and dividing the plurality of earthquake activity fault sections into strong earthquake breaking empty sections, fault movement locking sections, small and medium earthquake sparse sections and coulomb stress enhancement sections according to the information identification results specifically comprises the following steps:
acquiring information identification results of earthquake geology, geodetic survey, earthquake activity and numerical simulation of main faults of the boundary zone of the movable land block;
acquiring a plurality of earthquake activity fault sections according to the information identification result;
the plurality of seismic activity fault sections are divided into a strong-earthquake fracture space section, a fault movement locking section, a middle-small earthquake sparse section and a coulomb stress enhancement section.
As a further limitation of the technical solution of the embodiment of the present invention, the method further includes the steps of respectively subdividing the strong earthquake fracture space section, the fault motion locking section, the small and medium earthquake sparse section and the coulomb stress enhancement section to generate subdivision results, and performing different weight settings on the plurality of earthquake activity fault sections according to the subdivision results, where the generation of weight setting results specifically includes the following steps:
dividing the strong-earthquake breaking empty section into a high-release-rate empty section, a short-earthquake empty section and a non-empty section, determining a plurality of earthquake activity fault sections which are marked into the high-release-rate empty section, setting a first weight, determining a plurality of earthquake activity fault sections which are marked into the short-earthquake empty section, setting a second weight, determining a plurality of earthquake activity fault sections which are marked into the non-empty section, and setting a third weight;
dividing the fault movement locking section into a strong locking section, a partial locking section, an unblocking section and a non-data section, determining a plurality of seismic activity fault sections marked into the strong locking section, setting a fourth weight, determining a plurality of seismic activity fault sections marked into the partial locking section, setting a fifth weight, determining a plurality of seismic activity fault sections marked into the unblocking section, setting a sixth weight, determining a plurality of seismic activity fault sections marked into the non-data section, and setting a seventh weight;
dividing the small and medium-sized earthquake sparse sections into sparse sections and non-sparse sections, determining a plurality of earthquake activity fault sections marked into the sparse sections, setting eighth weights, determining a plurality of earthquake activity fault sections marked into the non-sparse sections, and setting ninth weights;
the coulomb stress enhancement section is divided into a significant enhancement section, an enhancement section and a non-significant enhancement section, a plurality of seismic-motion fault sections marked into the significant enhancement section are determined, a tenth weight is set, a plurality of seismic-motion fault sections marked into the enhancement section are determined, an eleventh weight is set, a plurality of seismic-motion fault sections marked into the non-significant enhancement section are determined, and a twelfth weight is set.
As a further limitation of the technical solution of the embodiment of the present invention, adding weights corresponding to each of the seismic activity fault sections according to the subdivision result and the weight setting result, and comprehensively obtaining risk coefficients corresponding to a plurality of the seismic activity fault sections specifically includes the following steps:
determining a plurality of set weights corresponding to each seismic activity fault section according to the subdivision result and the weight setting result;
and adding a plurality of set weights corresponding to each seismic activity fault section, and comprehensively obtaining a plurality of risk coefficients corresponding to the seismic activity fault sections.
As a further limitation of the technical solution of the embodiment of the present invention, the first weight is set to 3.0, the second weight is set to 2.0, the third weight is set to 0.0, the fourth weight is set to 3.0, the fifth weight is set to 1.5, the sixth weight is set to 0.0, the seventh weight is set to 0.0, the eighth weight is set to 1.0, the ninth weight is set to 0.0, the tenth weight is set to 3.0, the eleventh weight is set to 1.5, and the twelfth weight is set to 0.0.
As a further limitation of the technical scheme of the embodiment of the invention, the coulomb stress enhancement of the significant enhancement section is not less than 0.1MPa, the coulomb stress enhancement of the enhancement section is between 0.01 and 0.1MPa, and the coulomb stress enhancement of the non-significant enhancement section is less than 0.01MPa.
The system comprises an active fault dividing unit, a subdivision weight setting unit and a risk coefficient generating unit, wherein:
the movable fault dividing unit is used for acquiring information identification results of earthquake geology, geodetic measurement, earthquake activity and numerical simulation of main faults of the boundary zone of the movable land block, and dividing a plurality of earthquake movable fault sections into a strong earthquake breaking empty section, a fault movement locking section, a middle and small earthquake sparse section and a coulomb stress enhancement section according to the information identification results;
the subdivision weight setting unit is used for respectively subdividing the strong earthquake fracture space section, the fault movement locking section, the small and medium earthquake sparse section and the coulomb stress enhancement section to generate subdivision results, and carrying out different weight settings on the plurality of earthquake movement fault sections according to the subdivision results to generate weight setting results;
and the risk coefficient generation unit is used for adding the weights corresponding to each seismic activity fault section according to the subdivision result and the weight setting result, and comprehensively obtaining a plurality of risk coefficients corresponding to the seismic activity fault sections.
As a further limitation of the technical solution of the embodiment of the present invention, the active fault dividing unit specifically includes:
the recognition result acquisition module is used for acquiring information recognition results of earthquake geology, geodetic measurement, earthquake activity and numerical simulation of the main faults of the boundary zone of the movable land block;
the active fault acquisition module is used for acquiring a plurality of earthquake active fault sections according to the information identification result;
the movable fault dividing module is used for dividing the plurality of earthquake movable fault sections into a strong earthquake fracture empty section, a fault movement locking section, a middle and small earthquake sparse section and a coulomb stress enhancement section.
As a further limitation of the technical solution of the embodiment of the present invention, the subdivision weight setting unit specifically includes:
the first subdivision weight module is used for dividing the strong-earthquake breaking empty section into a high-release-rate empty section, a short-earthquake empty section and a non-empty section, determining a plurality of earthquake activity broken layer sections which are marked into the high-release-rate empty section, setting a first weight, determining a plurality of earthquake activity broken layer sections which are marked into the short-earthquake empty section, setting a second weight, determining a plurality of earthquake activity broken layer sections which are marked into the non-empty section, and setting a third weight;
the second subdivision weight module is used for dividing the fault movement locking section into a strong locking section, a partial locking section, an unblocking section and a non-data section, determining a plurality of seismic activity fault sections marked into the strong locking section, setting a fourth weight, determining a plurality of seismic activity fault sections marked into the partial locking section, setting a fifth weight, determining a plurality of seismic activity fault sections marked into the unblocking section, setting a sixth weight, determining a plurality of seismic activity fault sections marked into the non-data section, and setting a seventh weight;
the third subdivision weight module is used for dividing the small and medium-sized earthquake sparse sections into sparse sections and non-sparse sections, determining a plurality of earthquake activity fault sections marked into the sparse sections, setting eighth weights, determining a plurality of earthquake activity fault sections marked into the non-sparse sections, and setting ninth weights;
and the fourth subdivision weight module is used for dividing the coulomb stress enhancement section into a significant enhancement section, an enhancement section and a non-significant enhancement section, determining a plurality of seismic activity fault sections marked into the significant enhancement section, setting tenth weight, determining a plurality of seismic activity fault sections marked into the enhancement section, setting eleventh weight, determining a plurality of seismic activity fault sections marked into the non-significant enhancement section, and setting twelfth weight.
As a further limitation of the technical solution of the embodiment of the present invention, the risk coefficient generating unit specifically includes:
the weight determining module is used for determining a plurality of setting weights corresponding to each seismic activity fault section according to the subdivision result and the weight setting result;
and the weight adding module is used for adding a plurality of set weights corresponding to each seismic activity fault section and comprehensively obtaining a plurality of risk coefficients corresponding to the seismic activity fault sections.
Compared with the prior art, the invention has the beneficial effects that:
the method comprises the steps of dividing a plurality of seismic activity fault sections into a strong-seismic fracture space section, a fault movement locking section, a middle-small seismic sparse section and a coulomb stress enhancement section; subdividing and setting weights of the strong earthquake breaking empty section, the fault movement locking section, the small and medium earthquake sparse section and the coulomb stress enhancement section respectively; and adding the weights corresponding to the seismic activity fault sections, and comprehensively obtaining the risk coefficients corresponding to the seismic activity fault sections. The method can divide the strong-earthquake fracture space section, the fault motion locking section, the middle-small earthquake sparse section and the coulomb stress enhancement section of the earthquake activity fracture section, subdivide the sections, set different weights, comprehensively add the weights to obtain the corresponding danger coefficients of the different earthquake activity fracture sections, and avoid the false report phenomenon of single-method identification.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the following description will briefly introduce the drawings that are needed in the embodiments or the description of the prior art, and it is obvious that the drawings in the following description are only some embodiments of the present invention.
Fig. 1 shows a flowchart of a method provided by an embodiment of the present invention.
FIG. 2 shows a flow chart of the active fault identification partitioning in the method provided by the embodiment of the invention.
Fig. 3 shows a flowchart of subdivision and weight setting in the method provided by the embodiment of the invention.
Fig. 4 shows a flow chart of fault correspondence weight addition in the method provided by the embodiment of the invention.
Fig. 5 shows an application architecture diagram of a system provided by an embodiment of the present invention.
Fig. 6 shows a block diagram of an active fault dividing unit in the system according to an embodiment of the present invention.
Fig. 7 shows a block diagram of the structure of a subdivision weight setting unit in the system according to the embodiment of the present invention.
Fig. 8 is a block diagram showing a configuration of a risk coefficient generation unit in the system according to the embodiment of the present invention.
Fig. 9 shows a specific weight allocation rule diagram provided by an embodiment of the present invention.
Detailed Description
The present invention will be described in further detail with reference to the drawings and examples, in order to make the objects, technical solutions and advantages of the present invention more apparent. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention.
It will be appreciated that the methods available for inter-seismic late determination are, in summary, the most representative: a strong earthquake fracture empty section of earthquake geology, a fault movement locking section measured in the earth, a middle and small earthquake sparse section of earthquake activity, a coulomb stress enhancement significant section of numerical simulation and the like. From the past earthquake examples, the strong earthquake crack empty section, the fault movement locking section and the fault coulomb stress enhancement section are common earthquake location prediction methods, the earthquake examples of the middle and small earthquake sparse sections are relatively fewer, meanwhile, the strong earthquake identification results are all 'necessary insufficient' criteria for determining important dangerous areas of earthquakes, and the phenomenon of false report exists in different degrees through single method identification.
In order to solve the problems, the embodiment of the invention divides a plurality of earthquake activity fault sections into a strong earthquake crack empty section, a fault movement locking section, a middle and small earthquake sparse section and a coulomb stress enhancement section; subdividing and setting weights of the strong earthquake breaking empty section, the fault movement locking section, the small and medium earthquake sparse section and the coulomb stress enhancement section respectively; and adding the weights corresponding to the seismic activity fault sections, and comprehensively obtaining the risk coefficients corresponding to the seismic activity fault sections. The method can divide the strong-earthquake fracture space section, the fault motion locking section, the middle-small earthquake sparse section and the coulomb stress enhancement section of the earthquake activity fracture section, subdivide the sections, set different weights, comprehensively add the weights to obtain the corresponding danger coefficients of the different earthquake activity fracture sections, and avoid the false report phenomenon of single-method identification.
Fig. 1 shows a flowchart of a method provided by an embodiment of the present invention.
Specifically, the method for comprehensively judging the risk coefficient of the earthquake activity fault section specifically comprises the following steps:
step S101, obtaining information identification results of earthquake geology, geodetic survey, earthquake activity and numerical simulation of main faults of the boundary zone of the movable block, and dividing a plurality of earthquake activity fault sections into strong earthquake breaking empty sections, fault movement locking sections, small and medium earthquake sparse sections and coulomb stress enhancement sections according to the information identification results.
In the embodiment of the invention, based on earthquake geology, geodetic measurement, earthquake activity and numerical simulation data, the information identification result of main faults of the boundary zone of the movable land is obtained, and a plurality of earthquake activity fault sections are determined according to the information identification result, so that the plurality of earthquake activity fault sections are divided into strong earthquake breaking empty sections, fault movement locking sections, middle and small earthquake sparse sections and coulomb stress enhancement sections.
It can be understood that the strong-shock fracture empty section refers to a fracture section in which strong-shock fracture does not occur for a long time in a large-scale structural band such as a plate block boundary band or a movable block boundary band; the fault movement locking section refers to a fault section in which the fault rebound is staggered after the strain is accumulated to a critical point to generate earthquake fluctuation due to the fact that the relative movement of the plots at the two sides of the pre-earthquake san andersoi fault and the fault is in a locking state, so that the accumulation of elastic strain energy in the crust near the fault layer is caused; the middle and small earthquake sparse sections are fault sections with the activity degree of less than 0.4 according to the middle and small earthquake catalogues, the fault nodes are used as the centers, the number of earthquakes with more than 3 levels in the section is counted, and the activity degree normalization processing is carried out according to the boundary zone of each fault section; for the coulomb stress enhancement section, a common method for calculating the fault stress accumulation level is coulomb stress change numerical simulation on a fault plane, and many vibration examples show that the coulomb stress change caused by strong vibration can be used for long-term prediction of a strong vibration site.
Specifically, fig. 2 shows a flowchart of an active fault identification division in the method provided by the embodiment of the invention.
In the preferred embodiment provided by the invention, the method for acquiring the information identification results of the earthquake geology, geodetic measurement, earthquake activity and numerical simulation of the main faults of the boundary zone of the movable block, and dividing a plurality of earthquake activity fault sections into strong earthquake breaking empty sections, fault movement locking sections, small and medium earthquake sparse sections and coulomb stress enhancement sections according to the information identification results specifically comprises the following steps:
and S1011, acquiring information identification results of earthquake geology, geodetic survey, earthquake activity and numerical simulation of the main faults of the boundary zone of the movable land block.
Step S1012, acquiring a plurality of earthquake activity fault sections according to the information identification result.
Step S1013, dividing the plurality of seismic activity fault sections into a strong-earthquake fracture space section, a fault motion locking section, a middle-small earthquake sparse section and a fault coulomb stress enhancement section.
Further, the method for comprehensively judging the risk coefficient of the earthquake activity fault section further comprises the following steps:
step S102, the strong earthquake fracture space section, the fault movement locking section, the small and medium earthquake sparse section and the coulomb stress enhancement section are subdivided respectively to generate subdivision results, and different weight settings are carried out on a plurality of earthquake activity fault sections according to the subdivision results to generate weight setting results.
In the embodiment of the invention, as shown in fig. 9, a specific weight distribution rule diagram provided by the embodiment of the invention is shown, a strong-earthquake breaking air section is divided into a high-release-rate air section, a short-earthquake air section and a non-air section, a plurality of earthquake activity air sections which are marked into the high-release-rate air section are determined, a first weight is set to be 3.0, a plurality of earthquake activity air sections which are marked into the short-earthquake air section are determined, a second weight is set to be 2.0, a plurality of earthquake activity air sections which are marked into the non-air section are determined, and a third weight is set to be 0.0; dividing a fault movement locking section into a strong locking section, a partial locking section, an unblocking section and a non-data section, determining a plurality of earthquake activity fault sections marked into the strong locking section, setting a fourth weight to be 3.0, determining a plurality of earthquake activity fault sections marked into the partial locking section, setting a fifth weight to be 1.5, determining a plurality of earthquake activity fault sections marked into the unblocking section, setting a sixth weight to be 0.0, determining a plurality of earthquake activity fault sections marked into the non-data section, and setting a seventh weight to be 0.0; dividing a middle and small earthquake sparse section into a sparse section and a non-sparse section, determining a plurality of earthquake activity fault sections marked into the sparse section, setting an eighth weight to be 1.0, determining a plurality of earthquake activity fault sections marked into the non-sparse section, and setting a ninth weight to be 0.0; the coulomb stress enhancement section is divided into a significant enhancement section, an enhancement section and a non-significant enhancement section, a plurality of seismic-motion fault sections marked into the significant enhancement section are determined, a tenth weight is set to 3.0, a plurality of seismic-motion fault sections marked into the enhancement section are determined, an eleventh weight is set to 1.5, a plurality of seismic-motion fault sections marked into the non-significant enhancement section are determined, and a twelfth weight is set to 0.0.
Specifically, fig. 3 shows a flowchart of subdivision and weight setting in the method provided by the embodiment of the invention.
In the preferred embodiment provided by the invention, the strong earthquake breaking empty section, the fault movement locking section, the small and medium earthquake sparse section and the coulomb stress enhancement section are respectively subdivided to generate subdivision results, different weight settings are carried out on a plurality of earthquake activity fault sections according to the subdivision results, and the generation of the weight setting results specifically comprises the following steps:
step S1021, dividing the strong-earthquake breaking empty section into a high-release-rate empty section, a lack-earthquake empty section and a non-empty section, determining a plurality of earthquake activity broken layer sections marked into the high-release-rate empty section, setting a first weight, determining a plurality of earthquake activity broken layer sections marked into the lack-earthquake empty section, setting a second weight, determining a plurality of earthquake activity broken layer sections marked into the non-empty section, and setting a third weight.
In the embodiment of the invention, the strong shock fracture air section is divided into: (1) The high-release-rate empty section, the last strong vibration release time is close to or exceeds the fault section of the recurrence period; (2) The lack of earthquake empty section has no relevant research result of the evanescent rate and no fault section recorded by the history strong earthquake; (3) And a non-empty section, namely a strong shock fracture fault section divided by the two types.
Step S1022, dividing the fault motion blocking section into a strong blocking section, a partial blocking section, an unblocking section and a non-data section, determining a plurality of seismic activity fault sections marked into the strong blocking section, setting a fourth weight, determining a plurality of seismic activity fault sections marked into the partial blocking section, setting a fifth weight, determining a plurality of seismic activity fault sections marked into the unblocking section, setting a sixth weight, determining a plurality of seismic activity fault sections marked into the non-data section, and setting a seventh weight.
In the embodiment of the present invention, the locking degree results can be classified into 4 types, namely, strong locking, partial locking, non-locking and no data. The strong locking means that the locking depth of the broken layer section is close to the depth of the pregnant and vibrating layer in the region, and the locking rate is not lower than 0.8; the partial locking means that the locking depth of the broken layer section is smaller than the depth of the zone of the pregnant and vibrating layer but larger than 10km, and the locking rate is not lower than 0.8, or the locking depth of the broken layer section is close to the depth of the zone of the pregnant and vibrating layer, and the locking rate is between 0.6 and 0.8; unblocking means that the blocking depth of the broken layer section is less than 5km; the data-free means that the geodetic data of the area near the fault section is insufficient, and a reliable inversion result cannot be obtained.
Step S1023, dividing the small and medium-sized earthquake sparse section into a sparse section and a non-sparse section, determining a plurality of earthquake activity fault sections marked into the sparse section, setting an eighth weight, determining a plurality of earthquake activity fault sections marked into the non-sparse section, and setting a ninth weight.
In the embodiment of the invention, a calculation node is taken every 5km along the fault trend, 15 km in front of and behind each calculation node and 25 km on both sides of each calculation node are statistical windows, the number of times of earthquakes above 3.0 level in each calculation window is calculated, the logarithm is taken, a density scanning result is obtained, the normalization processing of the earthquake activity degree is carried out according to the boundary zone where each fault section is located, and if the normalized earthquake activity degree is less than 0.4, the fault section is judged to be a sparse section of medium and small earthquakes.
Step S1024, dividing the coulomb stress enhancement section into a significant enhancement section, an enhancement section and a non-significant enhancement section, determining a plurality of seismic-motion fault sections marked into the significant enhancement section, setting a tenth weight, determining a plurality of seismic-motion fault sections marked into the enhancement section, setting an eleventh weight, determining a plurality of seismic-motion fault sections marked into the non-significant enhancement section, and setting a twelfth weight.
In the embodiment of the invention, the coulomb stress enhancement of the significant enhancement section is not less than 0.1Mpa; the coulomb stress enhancement of the enhancement section is 0.01-0.1MPa; the coulomb stress enhancement of the non-significant enhancement segment is less than 0.01MPa.
Further, the method for comprehensively judging the risk coefficient of the earthquake activity fault section further comprises the following steps:
and step S103, adding weights corresponding to each seismic activity fault section according to the subdivision result and the weight setting result, and comprehensively obtaining a plurality of risk coefficients corresponding to the seismic activity fault sections.
In the embodiment of the invention, according to the subdivision result and the weight setting result, a plurality of setting weights corresponding to each seismic activity fault section are matched, the plurality of setting weights corresponding to each seismic activity fault section are added to obtain the risk coefficient corresponding to each seismic activity fault section, and the risk coefficients corresponding to all the seismic activity fault sections are comprehensively obtained.
Specifically, fig. 4 shows a flow chart of adding corresponding weights of faults in the method provided by the embodiment of the invention.
In the preferred embodiment of the present invention, according to the subdivision result and the weight setting result, adding weights corresponding to each of the seismic activity fault sections, and synthesizing to obtain risk coefficients corresponding to a plurality of seismic activity fault sections, the method specifically includes the following steps:
and step S1031, determining a plurality of setting weights corresponding to each seismic activity fault section according to the subdivision result and the weight setting result.
And S1032, adding a plurality of set weights corresponding to each seismic activity fault section, and comprehensively obtaining a plurality of risk coefficients corresponding to the seismic activity fault sections.
For example: a fault strong shock fracture empty section belongs to a high-release rate empty section, namely 3.0; the fault movement locking section belongs to partial locking, namely 1.5; the middle and small earthquake sparse section belongs to a non-sparse section, namely 0.0; the coulomb stress enhancement section belongs to 0.05MPa, i.e. 1.5, then the risk factor result is 3.0+1.5+0.0+1.5=6.0.
Further, fig. 5 shows an application architecture diagram of the system provided by the embodiment of the present invention.
In another preferred embodiment of the present invention, a system for comprehensively determining risk coefficients of an earthquake activity fault section includes:
the movable fault dividing unit 101 is configured to obtain information identification results of seismic geology, geodetic measurement, seismic activity and numerical simulation of a main fault in a boundary zone of a movable land block, and divide a plurality of seismic activity fault sections into a strong-earthquake breaking empty section, a fault movement locking section, a middle-small earthquake sparse section and a coulomb stress enhancement section according to the information identification results.
Specifically, fig. 6 shows a block diagram of the structure of the active fault dividing unit 101 in the system according to the embodiment of the present invention.
In a preferred embodiment of the present invention, the active fault dividing unit 101 specifically includes:
the recognition result obtaining module 1011 is configured to obtain information recognition results of seismic geology, geodetic, seismic activity, and numerical simulation of the main fault in the boundary zone of the movable land parcel.
The active fault acquisition module 1012 is used for acquiring a plurality of seismic active fault sections according to the information identification result.
The movable fault dividing module 1013 is configured to divide each of the plurality of seismic movable fault sections into a strong-earthquake breaking void section, a fault movement locking section, a middle-small earthquake sparse section and a coulomb stress enhancement section.
Further, the system for comprehensively judging the risk coefficient of the earthquake activity fault section further comprises:
the subdivision weight setting unit 102 is configured to subdivide the strong earthquake breaking empty section, the fault motion locking section, the small and medium earthquake sparse section, and the coulomb stress enhancement section, respectively, to generate subdivision results, and perform different weight settings on the plurality of earthquake activity fault sections according to the subdivision results, to generate weight setting results.
Specifically, fig. 7 shows a block diagram of the subdivision weight setting unit 102 in the system according to the embodiment of the present invention.
In a preferred embodiment of the present invention, the subdivision weight setting unit 102 specifically includes:
the first subdivision weight module 1021 is configured to divide the strong-earthquake breaking space into a high-release-rate space, a lack-earthquake space and a non-space, determine a plurality of earthquake activity breaking sections marked into the high-release-rate space, set a first weight, determine a plurality of earthquake activity breaking sections marked into the lack-earthquake space, set a second weight, determine a plurality of earthquake activity breaking sections marked into the non-space, and set a third weight.
A second subdivision weight module 1022, configured to divide the fault motion blocking section into a strong blocking section, a partial blocking section, an unblocking section, and a no-data section, determine a plurality of seismic-motion segments that are marked into the strong blocking section, set a fourth weight, determine a plurality of seismic-motion segments that are marked into the partial blocking section, set a fifth weight, determine a plurality of seismic-motion segments that are marked into the unblocking section, set a sixth weight, determine a plurality of seismic-motion segments that are marked into the no-data section, and set a seventh weight.
And a third subdivision weight module 1023, configured to divide the small and medium seismic sparse segments into sparse segments and non-sparse segments, determine a plurality of seismic activity fracture segments marked into the sparse segments, set an eighth weight, determine a plurality of seismic activity fracture segments marked into the non-sparse segments, and set a ninth weight.
A fourth subdivision weight module 1024, configured to divide the coulomb stress enhancement segment into a significant enhancement segment, an enhancement segment, and a non-significant enhancement segment, determine a plurality of seismic-motion intervals that are marked into the significant enhancement segment, and set a tenth weight, determine a plurality of seismic-motion intervals that are marked into the enhancement segment, and set an eleventh weight, determine a plurality of seismic-motion intervals that are marked into the non-significant enhancement segment, and set a twelfth weight.
Further, the system for comprehensively judging the risk coefficient of the earthquake activity fault section further comprises:
and the risk coefficient generating unit 103 is configured to add weights corresponding to each of the seismic activity fault segments according to the subdivision result and the weight setting result, and comprehensively obtain risk coefficients corresponding to a plurality of seismic activity fault segments.
Specifically, fig. 8 shows a block diagram of the risk coefficient generation unit 103 in the system according to the embodiment of the present invention.
In a preferred embodiment provided by the present invention, the risk factor generation unit 103 specifically includes:
the weight determining module 1031 is configured to determine a plurality of set weights corresponding to each seismic activity fault section according to the subdivision result and the weight setting result.
And the weight adding module 1032 is configured to add the plurality of set weights corresponding to each seismic activity fault section, and comprehensively obtain a plurality of risk coefficients corresponding to the seismic activity fault sections.
It should be understood that, although the steps in the flowcharts of the embodiments of the present invention are shown in order as indicated by the arrows, these steps are not necessarily performed in order as indicated by the arrows. The steps are not strictly limited to the order of execution unless explicitly recited herein, and the steps may be executed in other orders. Moreover, at least some of the steps in various embodiments may include multiple sub-steps or stages that are not necessarily performed at the same time, but may be performed at different times, nor do the order in which the sub-steps or stages are performed necessarily performed in sequence, but may be performed alternately or alternately with at least a portion of the sub-steps or stages of other steps or other steps.
Those skilled in the art will appreciate that all or part of the processes in the methods of the above embodiments may be implemented by a computer program for instructing relevant hardware, where the program may be stored in a non-volatile computer readable storage medium, and where the program, when executed, may include processes in the embodiments of the methods described above. Any reference to memory, storage, database, or other medium used in the various embodiments provided herein may include non-volatile and/or volatile memory. The nonvolatile memory can include Read Only Memory (ROM), programmable ROM (PROM), electrically Programmable ROM (EPROM), electrically Erasable Programmable ROM (EEPROM), or flash memory. Volatile memory can include Random Access Memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in a variety of forms such as Static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double Data Rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous Link DRAM (SLDRAM), memory bus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), among others.
The technical features of the above-described embodiments may be arbitrarily combined, and all possible combinations of the technical features in the above-described embodiments are not described for brevity of description, however, as long as there is no contradiction between the combinations of the technical features, they should be considered as the scope of the description.
The foregoing examples illustrate only a few embodiments of the invention and are described in detail herein without thereby limiting the scope of the invention. It should be noted that it will be apparent to those skilled in the art that several variations and modifications can be made without departing from the spirit of the invention, which are all within the scope of the invention. Accordingly, the scope of protection of the present invention is to be determined by the appended claims.
The foregoing description of the preferred embodiments of the invention is not intended to be limiting, but rather is intended to cover all modifications, equivalents, and alternatives falling within the spirit and principles of the invention.
Claims (8)
1. The comprehensive judging method for the risk coefficient of the earthquake activity fault section is characterized by comprising the following steps of:
obtaining information identification results of earthquake geology, geodetic survey, earthquake activity and numerical simulation of main faults of the boundary zone of the movable block, and dividing a plurality of earthquake activity fault sections into strong earthquake fracture empty sections, fault movement locking sections, small and medium earthquake sparse sections and coulomb stress enhancement sections according to the information identification results;
subdividing the strong earthquake fracture empty section, the fault movement locking section, the middle and small earthquake sparse section and the coulomb stress enhancement section respectively to generate subdivision results, and setting different weights for a plurality of earthquake activity fault sections according to the subdivision results to generate weight setting results;
and adding weights corresponding to each seismic activity fault section according to the subdivision result and the weight setting result, and comprehensively obtaining a plurality of risk coefficients corresponding to the seismic activity fault sections.
2. The method for comprehensively determining the risk coefficient of the seismic activity fault section according to claim 1, wherein the steps of respectively subdividing the strong-earthquake fault section, the fault motion locking section, the small-medium earthquake sparse section and the coulomb stress enhancement section to generate subdivision results, and setting different weights for the plurality of seismic activity fault sections according to the subdivision results, and generating the weight setting results specifically comprise the following steps:
dividing the strong-earthquake breaking empty section into a high-release-rate empty section, a short-earthquake empty section and a non-empty section, determining a plurality of earthquake activity fault sections which are marked into the high-release-rate empty section, setting a first weight, determining a plurality of earthquake activity fault sections which are marked into the short-earthquake empty section, setting a second weight, determining a plurality of earthquake activity fault sections which are marked into the non-empty section, and setting a third weight;
dividing the fault movement locking section into a strong locking section, a partial locking section, an unblocking section and a non-data section, determining a plurality of seismic activity fault sections marked into the strong locking section, setting a fourth weight, determining a plurality of seismic activity fault sections marked into the partial locking section, setting a fifth weight, determining a plurality of seismic activity fault sections marked into the unblocking section, setting a sixth weight, determining a plurality of seismic activity fault sections marked into the non-data section, and setting a seventh weight;
dividing the small and medium-sized earthquake sparse sections into sparse sections and non-sparse sections, determining a plurality of earthquake activity fault sections marked into the sparse sections, setting eighth weights, determining a plurality of earthquake activity fault sections marked into the non-sparse sections, and setting ninth weights;
the coulomb stress enhancement section is divided into a significant enhancement section, an enhancement section and a non-significant enhancement section, a plurality of seismic-motion fault sections marked into the significant enhancement section are determined, a tenth weight is set, a plurality of seismic-motion fault sections marked into the enhancement section are determined, an eleventh weight is set, a plurality of seismic-motion fault sections marked into the non-significant enhancement section are determined, and a twelfth weight is set.
3. The method for comprehensively determining risk coefficients of an earthquake activity fault section according to claim 1, wherein the step of adding weights corresponding to each earthquake activity fault section according to the subdivision result and the weight setting result to comprehensively obtain the risk coefficients corresponding to a plurality of earthquake activity fault sections specifically comprises the following steps:
determining a plurality of set weights corresponding to each seismic activity fault section according to the subdivision result and the weight setting result;
and adding a plurality of set weights corresponding to each seismic activity fault section, and comprehensively obtaining a plurality of risk coefficients corresponding to the seismic activity fault sections.
4. The seismic-motion interval risk coefficient comprehensive decision method according to claim 2, wherein the first weight is set to 3.0, the second weight is set to 2.0, the third weight is set to 0.0, the fourth weight is set to 3.0, the fifth weight is set to 1.5, the sixth weight is set to 0.0, the seventh weight is set to 0.0, the eighth weight is set to 1.0, the ninth weight is set to 0.0, the tenth weight is set to 3.0, the eleventh weight is set to 1.5, and the twelfth weight is set to 0.0.
5. The method of claim 2, wherein the significant enhancement segment coulomb stress enhancement is not less than 0.1MPa, the enhancement segment coulomb stress enhancement is at [0.01MPa,0.1 MPa), and the non-significant enhancement segment coulomb stress enhancement is less than 0.01MPa.
6. The system for comprehensively judging the risk coefficient of the earthquake activity fault section is characterized by comprising an activity fault dividing unit, a subdivision weight setting unit and a risk coefficient generating unit, wherein:
the movable fault dividing unit is used for acquiring information identification results of earthquake geology, geodetic measurement, earthquake activity and numerical simulation of main faults of the boundary zone of the movable land block, and dividing a plurality of earthquake movable fault sections into a strong earthquake breaking empty section, a fault movement locking section, a middle and small earthquake sparse section and a coulomb stress enhancement section according to the information identification results;
the subdivision weight setting unit is used for respectively subdividing the strong earthquake fracture space section, the fault movement locking section, the small and medium earthquake sparse section and the coulomb stress enhancement section to generate subdivision results, and carrying out different weight settings on the plurality of earthquake movement fault sections according to the subdivision results to generate weight setting results;
and the risk coefficient generation unit is used for adding the weights corresponding to each seismic activity fault section according to the subdivision result and the weight setting result, and comprehensively obtaining a plurality of risk coefficients corresponding to the seismic activity fault sections.
7. The seismic activity interval risk coefficient comprehensive decision system according to claim 6, wherein the subdivision weight setting unit specifically includes:
the first subdivision weight module is used for dividing the strong-earthquake breaking empty section into a high-release-rate empty section, a short-earthquake empty section and a non-empty section, determining a plurality of earthquake activity broken layer sections which are marked into the high-release-rate empty section, setting a first weight, determining a plurality of earthquake activity broken layer sections which are marked into the short-earthquake empty section, setting a second weight, determining a plurality of earthquake activity broken layer sections which are marked into the non-empty section, and setting a third weight;
the second subdivision weight module is used for dividing the fault movement locking section into a strong locking section, a partial locking section, an unblocking section and a non-data section, determining a plurality of seismic activity fault sections marked into the strong locking section, setting a fourth weight, determining a plurality of seismic activity fault sections marked into the partial locking section, setting a fifth weight, determining a plurality of seismic activity fault sections marked into the unblocking section, setting a sixth weight, determining a plurality of seismic activity fault sections marked into the non-data section, and setting a seventh weight;
the third subdivision weight module is used for dividing the small and medium-sized earthquake sparse sections into sparse sections and non-sparse sections, determining a plurality of earthquake activity fault sections marked into the sparse sections, setting eighth weights, determining a plurality of earthquake activity fault sections marked into the non-sparse sections, and setting ninth weights;
and the fourth subdivision weight module is used for dividing the coulomb stress enhancement section into a significant enhancement section, an enhancement section and a non-significant enhancement section, determining a plurality of seismic activity fault sections marked into the significant enhancement section, setting tenth weight, determining a plurality of seismic activity fault sections marked into the enhancement section, setting eleventh weight, determining a plurality of seismic activity fault sections marked into the non-significant enhancement section, and setting twelfth weight.
8. The seismic activity interval risk coefficient comprehensive decision system according to claim 6, wherein the risk coefficient generation unit specifically comprises:
the weight determining module is used for determining a plurality of setting weights corresponding to each seismic activity fault section according to the subdivision result and the weight setting result;
and the weight adding module is used for adding a plurality of set weights corresponding to each seismic activity fault section and comprehensively obtaining a plurality of risk coefficients corresponding to the seismic activity fault sections.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210799016.8A CN115201909B (en) | 2022-07-06 | 2022-07-06 | Comprehensive judging method and system for risk coefficient of earthquake activity fault section |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210799016.8A CN115201909B (en) | 2022-07-06 | 2022-07-06 | Comprehensive judging method and system for risk coefficient of earthquake activity fault section |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115201909A CN115201909A (en) | 2022-10-18 |
| CN115201909B true CN115201909B (en) | 2023-12-19 |
Family
ID=83580378
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210799016.8A Active CN115201909B (en) | 2022-07-06 | 2022-07-06 | Comprehensive judging method and system for risk coefficient of earthquake activity fault section |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115201909B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115903035B (en) * | 2022-11-17 | 2023-08-29 | 中国地震局地震预测研究所 | Earthquake triggering probability determining method and system based on geological parameters and coulomb stress |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA1280502C (en) * | 1987-03-11 | 1991-02-19 | Tom C. Holloway | Method for selection of mining and drilling sites using synthesized three dimensional seismic data |
| JPH09145849A (en) * | 1995-11-13 | 1997-06-06 | Serata Shiyoosei | Method and device for determination and prediction of direct-under ground type giant earthquake |
| KR101834783B1 (en) * | 2017-06-01 | 2018-03-06 | 이범근 | 3D the Earth crust modelling method and earthquake early warning system |
| CN111596354A (en) * | 2020-05-11 | 2020-08-28 | 同济大学 | Seismic reservoir prediction method considering space constraint under machine learning framework |
| CN113009574A (en) * | 2021-03-01 | 2021-06-22 | 应急管理部国家自然灾害防治研究院 | Method for evaluating earthquake risk of active fault |
| CN114048676A (en) * | 2021-11-02 | 2022-02-15 | 中国石油大学(华东) | Seismic fault and horizon interpretation method and device based on multitask deep learning |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7460436B2 (en) * | 2005-12-05 | 2008-12-02 | The Board Of Trustees Of The Leland Stanford Junior University | Apparatus and method for hydraulic fracture imaging by joint inversion of deformation and seismicity |
-
2022
- 2022-07-06 CN CN202210799016.8A patent/CN115201909B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA1280502C (en) * | 1987-03-11 | 1991-02-19 | Tom C. Holloway | Method for selection of mining and drilling sites using synthesized three dimensional seismic data |
| JPH09145849A (en) * | 1995-11-13 | 1997-06-06 | Serata Shiyoosei | Method and device for determination and prediction of direct-under ground type giant earthquake |
| KR101834783B1 (en) * | 2017-06-01 | 2018-03-06 | 이범근 | 3D the Earth crust modelling method and earthquake early warning system |
| CN111596354A (en) * | 2020-05-11 | 2020-08-28 | 同济大学 | Seismic reservoir prediction method considering space constraint under machine learning framework |
| CN113009574A (en) * | 2021-03-01 | 2021-06-22 | 应急管理部国家自然灾害防治研究院 | Method for evaluating earthquake risk of active fault |
| CN114048676A (en) * | 2021-11-02 | 2022-02-15 | 中国石油大学(华东) | Seismic fault and horizon interpretation method and device based on multitask deep learning |
Non-Patent Citations (5)
| Title |
|---|
| Qiu Qiang et al.coulomb stress perturbation after great earthquakes in the Sumatran subduction zone: Potential impacts in the surrounding region.JOURNAL OF ASIAN EARTH SCIENCE.2019,第180卷全部. * |
| 段虎荣 ; 周仕勇 ; 李闰 ; 陈胜雷 ; 闫全超 ; .芦山地震断层的滑动分布与汶川地震断层的关系.地球物理学报.2020,(第01期),全文. * |
| 罗艳 ; 赵里 ; 曾祥方 ; 高原 ; .芦山地震序列震源机制及其构造应力场空间变化.中国科学:地球科学.2015,(第04期),全文. * |
| 芦山地震序列震源机制及其构造应力场空间变化;罗艳;赵里;曾祥方;高原;;中国科学:地球科学(第04期);全文 * |
| 芦山地震断层的滑动分布与汶川地震断层的关系;段虎荣;周仕勇;李闰;陈胜雷;闫全超;;地球物理学报(第01期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115201909A (en) | 2022-10-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Abdelnaby | Fragility curves for RC frames subjected to Tohoku mainshock-aftershocks sequences | |
| Jena et al. | Seismic hazard and risk assessment: a review of state-of-the-art traditional and GIS models | |
| Abo-El-Ezz et al. | Seismic fragility assessment of low-rise stone masonry buildings | |
| CN115201909B (en) | Comprehensive judging method and system for risk coefficient of earthquake activity fault section | |
| CN104153814A (en) | Micro-seismic multi-parameter early warning method for rock burst | |
| Katsanos et al. | Prediction of inelastic response periods of buildings based on intensity measures and analytical model parameters | |
| Valentini et al. | Relaxing segmentation on the Wasatch fault zone: Impact on seismic hazard | |
| CN115440008B (en) | Landslide radar wave monitoring method and system | |
| CN115390130A (en) | Coal mining high-energy microseismic event prediction method and device | |
| Astriana et al. | Assessing seismic performance of moment resisting frame and frame-shear wall system using seismic fragility curve | |
| Sousa et al. | Including multiple IMTs in the development of fragility functions for earthquake loss estimation | |
| CN111538071B (en) | Quantitative prediction method for displacement of steep dip stratified rock mass cavern group high side wall | |
| Joughin | Dealing with uncertainty and risk in the design of deep and high stress mining excavations | |
| CN112230275A (en) | Seismic waveform identification method and device and electronic equipment | |
| Park et al. | Fragility assessment model of building structures using characteristics of artificial aftershock motions | |
| Pejović et al. | Dependence of RC high-rise buildings response on the earthquake intensity | |
| Kyriakides et al. | Vulnerability curves for RC substandard buildings | |
| Mignan et al. | Harmonized approach to stress tests for critical infrastructures against natural hazards (STREST) | |
| CN114896660A (en) | Durability design method for underground engineering anchoring and supporting system | |
| Abdelhamid et al. | Seismic fragility and uncertainty analysis of concrete gravity dams under near-fault ground motions | |
| CN113139681A (en) | Neural network rock burst prediction method based on time series data | |
| Sborshchikovi et al. | The Analysis of Dynamical Changes and Local Seismic Activity of the Enguri Arch Dam | |
| Bogdanov et al. | A probabilistic model of seismicity: Kamchatka earthquakes | |
| Tierney et al. | Seismic exclusions and re-entry from a risk perspective | |
| CN117708600B (en) | Bridge blasting shock hazard prediction and evaluation method |
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 |