Earthquake intensity value calculation method, system and equipment and storage medium

Abstract

The embodiment of the invention discloses a seismic intensity value calculation method, system and device, and a storage medium, and relates to the technical field of seismic intensity monitoring, the influence of a topographic condition on seismic oscillation propagation can be fully considered by constructing a topographic amplification coefficient based on DEM data in combination with a slope correction factor, so that the first PGA distribution is corrected by using the topographic amplification coefficient, and the seismic intensity value calculation accuracy is improved. The defect that a traditional method is insufficient in consideration of topographic conditions is overcome; the field shear wave velocity Vs30 data and the basin effect correction coefficient are used for realizing refined characterization of the seismic oscillation amplification effects of different soil layer structures, so that the field amplification coefficient is calculated by using the field shear wave velocity Vs30 data and the basin effect correction coefficient, and the second PGA distribution is corrected by using the field amplification coefficient; and the precision of the seismic intensity value under the complex geological and topographic conditions is obviously improved.

Description

Technical Field

[0001] This invention relates to the field of earthquake intensity monitoring technology, and in particular to a method, system, device and storage medium for calculating earthquake intensity values. Background Technology

[0002] Seismic intensity is a crucial indicator for measuring the impact of earthquakes on the Earth's surface and engineering structures, serving as a vital basis for earthquake emergency response, disaster assessment, and seismic design of engineering projects. Existing rapid seismic intensity assessment techniques primarily rely on empirical attenuation relationships, estimating intensity distribution using parameters such as magnitude and epicentral distance. However, these methods have the following shortcomings: First, they do not adequately consider the influence of topographic conditions on seismic motion propagation; complex mountainous terrain significantly amplifies or attenuates seismic waves. Second, they lack refined processing of site conditions; different soil structures exhibit vastly different amplification effects on seismic motion, particularly in mountainous areas with complex terrain and diverse geological conditions, leading to significant errors in the calculated seismic intensity values. Summary of the Invention

[0003] In view of this, the present invention provides a method, system, device and storage medium for calculating earthquake intensity values.

[0004] The specific technical solution of the first embodiment of the present invention is as follows: a method for calculating seismic intensity values, the method comprising: acquiring seismic event information and site shear wave velocity (Vs30) data of a target area; obtaining a first PGA distribution of the target area based on the seismic event information and the site shear wave velocity (Vs30) data; the first PGA distribution being a peak ground acceleration (PGA) value; acquiring DEM data of the target area; obtaining a topographic amplification factor of the target area based on the DEM data and a preset slope correction factor; correcting the first PGA distribution based on the topographic amplification factor to obtain a second PGA distribution; obtaining a site amplification factor of the target area based on the site shear wave velocity (Vs30) data and a preset basin effect correction factor; correcting the second PGA distribution based on the site amplification factor to obtain a third PGA distribution; and converting the third PGA distribution into a seismic intensity value of the target area based on a preset correspondence between PGA and seismic intensity.

[0005] Preferably, the first PGA distribution is obtained using the following formula:

[0006]

[0007] Among them, PGA base The first PGA distribution is defined as follows: C1, C2, C3, C4, and C5 are preset attenuation coefficients; M is the earthquake magnitude in the earthquake event information; R is the epicenter distance in the earthquake event information; S is the site type parameter of the target area; and V is the ground type parameter of the target area.s30 The data refers to the site shear wave velocity Vs30.

[0008] Preferably, obtaining the terrain magnification factor of the target area based on the DEM data and the preset slope correction factor includes: obtaining the elevation of a target point in the target area based on the DEM data; the target point is any point in the target area; obtaining the terrain shape ratio of the target point based on the elevation of the target point, the horizontal distance component from the target point along the maximum slope direction to the foot of the slope, and the vertical distance component from the target point along the maximum slope direction to the foot of the slope; and obtaining the terrain magnification factor of the target area based on the terrain shape ratio and the preset slope correction factor.

[0009] Preferably, the terrain magnification factor is obtained using the following formula:

[0010] A h =1 + 0.4 × arctan(4SR) × F s

[0011] Among them, A h Where SR is the terrain magnification factor, F is the terrain shape ratio, and F is the terrain magnification factor. s This is the preset slope correction factor.

[0012] Preferably, the second PGA distribution is obtained using the following formula:

[0013] PGA topo =PGA base ×A h ×(1+0.1×cos(θ-φ))

[0014] Among them, PGA topo This is the second PGA distribution, PGA base For the first PGA distribution, A h φ is the terrain magnification factor, θ is the azimuth angle from the epicenter to the target point, and φ is the mountain orientation angle of the target point.

[0015] Preferably, after obtaining the seismic intensity value of the target area, the method further includes: obtaining the fault zone enhancement amplitude and the distance from the target point to the fault zone in the target area; and correcting the seismic intensity value based on the fault zone enhancement amplitude and the distance from the target point to the fault zone to obtain a corrected seismic intensity value.

[0016] Preferably, the corrected seismic intensity value is obtained using the following formula:

[0017]

[0018] Among them, I faultThe corrected seismic intensity value is denoted as I, the seismic intensity value is ΔI, the fault zone enhancement amplitude is d, the distance from the target point to the fault zone is d, and the preset influence range parameter is σ.

[0019] The specific technical solution of the second embodiment of the present invention is as follows: a seismic intensity value calculation system, the system comprising: a data acquisition module, a first distribution calculation module, a first coefficient acquisition module, a first correction module, a second coefficient acquisition module, a second correction module, and an intensity value output module; the data acquisition module is used to acquire seismic event information and site shear wave velocity Vs30 data of a target area; the first distribution calculation module is used to obtain a first PGA distribution of the target area based on the seismic event information and the site shear wave velocity Vs30 data; the first PGA distribution is the peak ground acceleration (PGA) value; the first coefficient acquisition module is used to acquire DEM data of the target area. The system obtains the topographic amplification factor of the target area based on the DEM data and a preset slope correction factor; the first correction module corrects the first PGA distribution based on the topographic amplification factor to obtain a second PGA distribution; the second coefficient acquisition module obtains the site amplification factor of the target area based on the site shear wave velocity Vs30 data and a preset basin effect correction factor; the second correction module corrects the second PGA distribution based on the site amplification factor to obtain a third PGA distribution; and the intensity value output module converts the third PGA distribution into the seismic intensity value of the target area based on a preset correspondence between PGA and seismic intensity.

[0020] The specific technical solution of the third embodiment of the present invention is as follows: an earthquake intensity value calculation device, including a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the processor performs the steps of the method as described in any one of the first embodiments of this application.

[0021] The specific technical solution of the fourth embodiment of the present invention is as follows: a computer-readable storage medium storing a computer program, wherein when the computer program is executed by a processor, the processor performs the steps of the method as described in any one of the first embodiments of this application.

[0022] Implementing the embodiments of the present invention will have the following beneficial effects:

[0023] This invention constructs a terrain amplification factor based on DEM data and a slope correction factor, which can fully consider the influence of terrain conditions on earthquake propagation. Therefore, the terrain amplification factor is used to correct the first PGA distribution, making up for the shortcomings of traditional methods in not considering terrain conditions. The site shear wave velocity Vs30 data and the basin effect correction factor enable a refined characterization of the earthquake amplification effect of different soil structures. Therefore, the site amplification factor is calculated using the site shear wave velocity Vs30 data and the basin effect correction factor, and the site amplification factor is used to correct the second PGA distribution, which significantly improves the accuracy of earthquake intensity values ​​under complex geological and topographical conditions. Attached Figure Description

[0024] To more clearly illustrate the technical solutions in the embodiments of this application or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0025] Figure 1 This is a flowchart illustrating the steps involved in calculating earthquake intensity values.

[0026] Figure 2 This is a schematic diagram of the earthquake intensity calculation system.

[0027] Figure 3 This is a diagram of the internal structure of a computer device.

[0028] Among them, 201 is the data acquisition module; 202 is the first distribution calculation module; 203 is the first coefficient acquisition module; 204 is the first correction module; 205 is the second coefficient acquisition module; 206 is the second correction module; and 207 is the intensity value output module. Detailed Implementation

[0029] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, and not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those of ordinary skill in the art without creative effort are within the scope of protection of this application.

[0030] The terms "first," "second," etc., used in the specification, claims, and drawings of this application are used to distinguish different objects, not to describe a specific order. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover non-exclusive inclusion. For example, a process, method, system, product, or apparatus that includes a series of steps or modules is not limited to the listed steps or modules, but may optionally include steps or modules not listed, or may optionally include other steps or modules inherent to such processes, methods, products, or apparatus.

[0031] In this document, the term "embodiment" means that a particular feature, structure, or characteristic described in connection with an embodiment may be included in at least one embodiment of this application. The appearance of this phrase in various places throughout the specification does not necessarily refer to the same embodiment, nor is it a separate or alternative embodiment mutually exclusive with other embodiments. It will be explicitly and implicitly understood by those skilled in the art that the embodiments described herein can be combined with other embodiments.

[0032] Please see Figure 1 This is a flowchart illustrating the steps of a seismic intensity value calculation method in the first embodiment of this application, which improves the accuracy of seismic intensity values ​​under complex geological and topographical conditions. The method includes:

[0033] Step 101: Obtain seismic event information and site shear wave velocity Vs30 data for the target area;

[0034] Specifically, earthquake event information may include three elements: epicenter location coordinates, time of occurrence, and magnitude. The epicenter location coordinates include longitude and latitude information, and the magnitude is either Richter magnitude or moment magnitude. High-precision digital elevation model (DEM) data of the target area, site shear wave velocity (Vs30) data, active fault spatial distribution data, and historical earthquake event database are loaded. The DEM data resolution is not less than 30 meters, and the Vs30 data represents the average shear wave velocity within 30 meters below the surface.

[0035] The earthquake event information is acquired as follows: real-time earthquake information released by the China Earthquake Networks Center is automatically obtained through web crawling technology, with the crawling frequency set to once every 5 minutes. The acquired earthquake information undergoes format standardization and validity verification, eliminating duplicate records and abnormal data. First, a data structure analysis of the target area is established to identify key HTML tags and JSON data fields containing earthquake information. For example, for data from the China Earthquake Networks Center, XML format data containing information such as epicenter latitude and longitude, time of occurrence, magnitude, and focal depth is parsed. For USGS data, GeoJSON format earthquake catalog data is processed. To ensure the stability of data acquisition, a retry mechanism and fault tolerance strategy are adopted. When an access anomaly occurs in a data source, acquisition is automatically retried, and an anomaly log is recorded. Simultaneously, a data source priority mechanism is established; when the primary data source is unavailable, it automatically switches to a backup data source to ensure the continuity of earthquake information acquisition. When an earthquake event with a magnitude exceeding a preset threshold is detected, a rapid intensity calculation process is automatically triggered. Manual input of the three earthquake elements is supported for historical earthquake reproduction or hypothetical scenario analysis. Seismic networks may use different coordinate systems, such as WGS84, Beijing54, and CGCS2000. Through various coordinate system transformation algorithms, the coordinate system type of the input data can be automatically identified and uniformly converted to the WGS84 coordinate system. The transformation process adopts a seven-parameter transformation model to ensure that the transformation accuracy meets the requirements of subsequent calculations.

[0036] The time of an earthquake may be expressed in UTC time, local time, or other time zones. Based on the time zone settings of the data source, all time information is uniformly converted to UTC time and stored in the ISO8601 standard format.

[0037] Standardization of magnitude data includes magnitude type identification and numerical correction. Different seismic networks use different magnitude scales, such as Richter magnitude (ML), moment magnitude (Mw), and surface wave magnitude (Ms). A magnitude conversion database is established to uniformly convert different types of magnitudes into moment magnitude (Mw). The conversion relationship adopts the standard formula recommended by the International Association for Seismology and Physics of the Earth's Interior (IASPEI).

[0038] Step 102: Based on the earthquake event information and the site shear wave velocity Vs30 data, obtain the first PGA distribution of the target area; the first PGA distribution is the peak ground acceleration (PGA) value;

[0039] Specifically, based on the ground motion attenuation relationship model of the target area, the peak ground acceleration (PGA) values ​​of each spatial grid point in the target area are calculated to obtain the first PGA distribution. The ground motion attenuation relationship model comprehensively considers the influence of magnitude, epicentral distance and site conditions on the ground motion intensity.

[0040] Step 103: Obtain the DEM data of the target area, and obtain the terrain magnification factor of the target area based on the DEM data and the preset slope correction factor;

[0041] Step 104: Correct the first PGA distribution according to the terrain magnification factor to obtain the second PGA distribution;

[0042] Specifically, the terrain magnification factor of the target area is calculated based on high-precision DEM data. The terrain magnification factor is then used to correct the terrain effect of the first PGA distribution of each grid point to obtain the terrain-corrected second PGA distribution.

[0043] Step 105: Obtain the site amplification factor of the target area based on the site shear wave velocity Vs30 data and the preset basin effect correction factor;

[0044] Step 106: Correct the second PGA distribution according to the site magnification factor to obtain the third PGA distribution;

[0045] Specifically, the site amplification factor is obtained by analyzing the soil structure characteristics of the site based on Vs30 data. The site effect is then corrected on the second PGA distribution after topographic correction using the site amplification factor to obtain the third PGA distribution after comprehensive correction.

[0046] Step 107: Based on the preset correspondence between PGA and seismic intensity, convert the third PGA distribution into the seismic intensity value of the target area.

[0047] Specifically, based on the correspondence between PGA and seismic intensity, the corrected third PGA distribution is converted into seismic intensity values, and the spatial distribution of seismic intensity in the target area is generated.

[0048] Intensity conversion and spatial distribution calculations include: establishing the third PGA distribution. site The conversion relationship between seismic intensity I and seismic intensity I is represented by a piecewise linear model:

[0049]

[0050] Different PGA ranges correspond to different intensity levels. The determination of the conversion factor takes into account the following factors: the degree of building damage corresponding to different intensity levels; people's perception and reaction under different intensity levels; surface damage phenomena under different intensity levels; and the characteristics of seismic motion parameters under different intensity levels.

[0051] First paragraph (PGA) site<35gal) corresponds to intensity III-V. Earthquakes in this range are relatively weak, mainly manifested as slight tremors and very little building damage. The conversion factor of 1.5 reflects a strong positive correlation between the intensity and the logarithm of PGA in this range.

[0052] Second paragraph (35≤PGA) site <100gal) corresponds to intensity V-VII, within which significant building damage begins to appear, especially for buildings with poor seismic performance; the conversion factor of 2.0 indicates that the intensity is more sensitive to changes in PGA within this range.

[0053] Third paragraph (100≤PGA) site <300 gal) corresponds to an intensity of VII-IX, which would cause severe building damage and is a key focus of earthquake hazard assessment; the reduction of the conversion factor to 1.2 reflects the complexity of the conversion relationship in the high-intensity range.

[0054] Fourth paragraph (PGA) site ≥300 gal corresponds to an intensity of IX or higher, and earthquakes of this intensity will cause devastating damage; the further reduction of the conversion factor to 0.8 reflects the saturation effect of intensity increase under extremely strong earthquakes.

[0055] The method in this embodiment, based on DEM data and combined with slope correction factors, constructs a terrain amplification coefficient that fully considers the influence of terrain conditions on earthquake propagation. Therefore, the terrain amplification coefficient is used to correct the first PGA distribution, making up for the shortcomings of traditional methods in not considering terrain conditions. The site shear wave velocity Vs30 data and basin effect correction coefficient enable a refined characterization of the earthquake amplification effect of different soil structures. Therefore, the site amplification coefficient is calculated using the site shear wave velocity Vs30 data and basin effect correction coefficient, and the site amplification coefficient is used to correct the second PGA distribution, significantly improving the accuracy of seismic intensity values ​​under complex geological and topographical conditions.

[0056] In a specific embodiment, the first PGA distribution is obtained using the following formula:

[0057]

[0058] Among them, PGA base The first PGA distribution is defined as follows: C1, C2, C3, C4, and C5 are preset attenuation coefficients; M is the earthquake magnitude in the earthquake event information; R is the epicenter distance in the earthquake event information; S is the site type parameter of the target area; and V is the ground type parameter of the target area. s30 The data refers to the site shear wave velocity Vs30.

[0059] Specifically, C1, C2, C3, C4, and C5 are attenuation coefficients obtained by regression of strong earthquake records for the target area. Among them, C1 is a constant term, C2 is the magnitude influence coefficient, C3 is the distance attenuation coefficient, C4 is the near-field saturation distance parameter, and C5 is the site influence coefficient.

[0060] The method for determining the attenuation coefficient includes: collecting and organizing historical strong earthquake records and seismic intensity survey data within the target area, with a data sample of no less than 100 strong earthquake records. An elliptical attenuation model is established, considering the attenuation differences of seismic waves along and perpendicular to the fault strike; the ratio of the major and minor axes of the ellipse is determined based on the regional tectonic stress field.

[0061] The attenuation coefficients were obtained by fitting the parameters using a two-step least squares method. The values ​​of each attenuation coefficient are: C1∈[-2.0,2.0], C2∈[0.5,1.5], C3∈[-2.5,-1.0], C4∈[5,20], and C5∈[-0.5,0.5].

[0062] The attenuation relationships along the major and minor axes of the ellipse are expressed as follows:

[0063]

[0064] Among them, R major and R manor , respectively, are the equivalent distances along the major and minor axes of the ellipse, with subscripts a and b corresponding to the regression coefficients along the major and minor axes, respectively, and S being the area of ​​the ellipse.

[0065] For areas lacking measured Vs30 data, an empirical estimation method based on topography is used. Studies show a correlation between topographic gradient and Vs30; Vs30 values ​​are typically higher in mountainous and hilly areas, and lower in plains and basins. The empirical formula is as follows:

[0066] log(Vs30) = a + b × log(slope) + c × log(elevation); where slope is the terrain slope, elevation is the altitude, and a, b, and c are regression coefficients determined based on the regional geological conditions.

[0067] The PGA calculation uses a regular gridding method, which divides the target area into regular grid cells and calculates the PGA value at the center point of each grid. The choice of grid resolution needs to balance the calculation accuracy and calculation efficiency, and a grid spacing of 500 meters × 500 meters or 1 kilometer × 1 kilometer is used.

[0068] The calculation area is determined based on the earthquake magnitude and user requirements, using either a circular or rectangular area. A regular grid is generated within the determined area, with each grid point having unique latitude and longitude coordinates. The distance from each grid point to the epicenter is calculated, taking into account the influence of the Earth's curvature, using the great circle distance formula. The corresponding Vs30 value is extracted for each grid point to determine the site type parameter S.

[0069] The PGA value for each grid point is calculated based on the attenuation relationship model. To optimize computational efficiency, a parallel computing strategy is employed. Since the PGA calculations for each grid point are independent, they are well-suited for parallel processing. Based on the availability of computing resources, grid points are allocated to different computation threads, significantly improving computational speed. Distance calculations consider the influence of the Earth's curvature, employing a more accurate great circle distance formula.

[0070] R = R earth ×arccos(sin(φ1)sin(φ2)+cos(φ1)cos(φ2)cos(Δλ)); where, R earth φ is the Earth's radius (6371 km), φ1 and φ2 are the latitudes of the epicenter and the calculation point, respectively, and Δλ is the difference in longitude.

[0071] In a specific embodiment, obtaining the terrain magnification coefficient of the target area based on the DEM data and a preset slope correction factor includes: obtaining the elevation of a target point in the target area based on the DEM data; the target point is any point in the target area; obtaining the terrain shape ratio of the target point based on the elevation of the target point, the horizontal distance component from the target point along the maximum slope direction to the foot of the slope, and the vertical distance component from the target point along the maximum slope direction to the foot of the slope; and obtaining the terrain magnification coefficient of the target area based on the terrain shape ratio and the preset slope correction factor.

[0072] Specifically, based on DEM data, the terrain parameters of target points within the target area are calculated, and the formula for calculating the terrain height H is as follows: Among them, Z p The elevation of the grid point corresponding to the target point. The average elevation within a 1000-meter radius around this point; the method for calculating the terrain feature length L is as follows:

[0073] Where Δx and Δy are the horizontal distance components from the toe of the slope along the direction of maximum slope.

[0074] The formula for calculating the terrain shape ratio (SR) is: The terrain magnification factor of the target area is obtained by using the terrain shape ratio SR and the preset slope correction factor.

[0075] In a specific embodiment, the terrain magnification factor is obtained using the following formula:

[0076] A h =1 + 0.4 × arctan(4SR) × F s

[0077] Among them, A h Where SR is the terrain magnification factor, F is the terrain shape ratio, and F is the terrain magnification factor. s Let F be the preset slope correction factor, where F s Slope correction factor:

[0078]

[0079] In a specific embodiment, the second PGA distribution is obtained using the following formula:

[0080] PGA topo =PGA base ×A h ×(1+0.1×cos(θ-φ))

[0081] Among them, PGA topo This is the second PGA distribution, PGA base For the first PGA distribution, A h φ is the terrain amplification factor, θ is the azimuth angle from the epicenter to the target point, and φ is the mountain orientation angle of the target point, used to consider the influence of terrain orientation on earthquake motion amplification.

[0082] In a specific embodiment, the shear wave velocity Vs30 data of the target area is obtained, and the site magnification factor F is... t The calculation uses a piecewise function model: F t =F LIN +F NL F LIN For linear amplification terms: F NL For nonlinear amplification terms: In the formula, c3 = 0.03g, PGA rock The PGA value for the bedrock site;

[0083] The third PGA distribution after site correction site The calculation formula is:

[0084] PGA site =PGA topo ×exp(F t )×a b ; where a b Here is the basin effect correction factor, and the formula for the basin effect correction factor is:

[0085] In a specific embodiment, after obtaining the seismic intensity value of the target area, the method further includes: obtaining the fault zone enhancement amplitude and the distance from the target point to the fault zone in the target area; and correcting the seismic intensity value based on the fault zone enhancement amplitude and the distance from the target point to the fault zone to obtain a corrected seismic intensity value.

[0086] In a specific embodiment, the corrected seismic intensity value is obtained using the following formula:

[0087]

[0088] Among them, I fault The corrected seismic intensity value is denoted as I, the seismic intensity value is ΔI, the fault zone enhancement amplitude is d, the distance from the target point to the fault zone is d, and the preset influence range parameter is σ.

[0089] Specifically, active fault zones have a significant impact on the distribution of seismic intensity, often exhibiting anomalies in intensity near fault zones. The physical mechanisms underlying this impact include: fracturing rocks near fault zones, reducing seismic wave propagation speed; softer filling materials in fault zones, amplifying seismic motion; and the influence of fault zone geometry on seismic wave propagation paths.

[0090] Considering the influence of active fault zones on intensity distribution, the intensity I of the fault zone... fault It is necessary to correct for the earthquake intensity level I by using the fault zone enhancement amplitude: Where σ is the influence range parameter, which is taken as 5-10km.

[0091] In a specific embodiment, the method further includes determining the intensity I of the fault zone. fault It generates seismic intensity contour maps and thematic analysis reports, supporting data output and graphical display in multiple formats; it automatically loads multi-source base map data, including satellite remote sensing imagery, topographic shading maps, and administrative division annotation layers; it supports multi-layer overlay display, with layers including intensity contour lines, active fault distribution, historical earthquake epicenters, and infrastructure distribution; and it provides WebGIS-based online interactive functions, including spatial extent query, distance measurement, area statistics, and attribute information viewing.

[0092] The method in this embodiment also includes the calculation of multi-source combined effects. When multiple earthquake sources exist simultaneously, a probabilistic seismic hazard analysis method is used. When n sources exist simultaneously, the PGA distribution of the superimposed ground motion parameters from multiple sources is calculated. total The calculation formula is: Among them, PGA i Let w be the ground motion parameters generated by the i-th source.i The weighting coefficients are determined based on magnitude and distance: Considering the time correlation, a Poisson process model is adopted. The probability that an earthquake will occur at source i within time t is: P i (t)=1-exp(-λ i t), where λ i Let i be the annual average seismicity rate of earthquake source i;

[0093] The formula for calculating the maximum seismic intensity under multiple source combined action is:

[0094] I max =max{I1,I2,…I n}+βlog(n); where β is the multi-source superposition coefficient, with a value range of 0.1-0.3.

[0095] The continuous intensity distribution is generated using the Kriging interpolation method, and the variogram model is as follows:

[0096] Where q0 is the nugget effect, q1 is the arch height, a1 is the range, and h1 is the spatial lag distance;

[0097] The intensity distribution accuracy assessment uses a cross-validation method, and the accuracy index is:

[0098] The required RMSE is ≤ 0.5 degrees.

[0099] In a specific embodiment, please refer to Figure 2This is a schematic diagram of a seismic intensity value calculation system according to a second embodiment of this application. The system includes: a data acquisition module 201, a first distribution calculation module 202, a first coefficient acquisition module 203, a first correction module 204, a second coefficient acquisition module 205, a second correction module 206, and an intensity value output module 207. The data acquisition module 201 is used to acquire seismic event information and site shear wave velocity Vs30 data of the target area. The first distribution calculation module 202 is used to obtain a first PGA distribution of the target area based on the seismic event information and the site shear wave velocity Vs30 data. The first PGA distribution is the peak ground acceleration (PGA) value. The first coefficient acquisition module 203 is used to acquire the intensity value of the target area. The DEM data of the target area is used to obtain the topographic amplification factor of the target area based on the DEM data and a preset slope correction factor; the first correction module 204 is used to correct the first PGA distribution based on the topographic amplification factor to obtain a second PGA distribution; the second coefficient acquisition module 205 is used to obtain the site amplification factor of the target area based on the site shear wave velocity Vs30 data and a preset basin effect correction factor; the second correction module 206 is used to correct the second PGA distribution based on the site amplification factor to obtain a third PGA distribution; the intensity value output module 207 is used to convert the third PGA distribution into the seismic intensity value of the target area according to the preset correspondence between PGA and seismic intensity.

[0100] In this embodiment, the system constructs a terrain amplification coefficient based on DEM data and a slope correction factor, which can fully consider the influence of terrain conditions on earthquake propagation. Therefore, the terrain amplification coefficient is used to correct the first PGA distribution, making up for the shortcomings of traditional methods in not considering terrain conditions. The site shear wave velocity Vs30 data and the basin effect correction coefficient enable a refined characterization of the earthquake amplification effect of different soil structures. Therefore, the site amplification coefficient is calculated using the site shear wave velocity Vs30 data and the basin effect correction coefficient, and the site amplification coefficient is used to correct the second PGA distribution, which significantly improves the accuracy of earthquake intensity values ​​under complex geological and topographical conditions.

[0101] In a specific embodiment, the third embodiment of this application provides an earthquake intensity value calculation device, including a memory and a processor. The memory stores a computer program, and when the computer program is executed by the processor, the processor performs the steps of the method as described in any one of the first embodiments of this application.

[0102] In a specific embodiment, the fourth embodiment of this application provides a computer-readable storage medium storing a computer program, which, when executed by a processor, causes the processor to perform the steps of the method as described in any one of the first embodiments of this application.

[0103] Figure 3 An internal structural diagram of a computer device in one embodiment is shown. This computer device can specifically be a terminal or a server. Please refer to... Figure 3 The computer device includes a processor, memory, etc., connected via a system bus. The memory includes a non-volatile storage medium and internal memory. The non-volatile storage medium stores an operating system and may also store a computer program. When executed by the processor, this computer program causes the processor to implement the method described in this embodiment. The internal memory may also store a computer program, which, when executed by the processor, causes the processor to perform the method described in this embodiment. Those skilled in the art will understand that... Figure 3 The structure shown is merely a block diagram of a portion of the structure related to the present application and does not constitute a limitation on the computer device to which the present application is applied. Specific computer devices may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.

[0104] The above embodiments merely illustrate several implementation methods of this application, and their descriptions are relatively specific and detailed. However, they should not be construed as limiting the scope of this patent application. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of this application, and these all fall within the protection scope of this application. Therefore, the protection scope of this patent application should be determined by the appended claims.

[0105] The above description is merely a preferred embodiment of the present invention and is not intended to limit the present invention in any other way. Any person skilled in the art may make changes or modifications to the above-disclosed technical content to create equivalent embodiments for application in other fields. However, any simple modifications, equivalent changes, and modifications made to the above embodiments based on the technical essence of the present invention without departing from the scope of the present invention shall still fall within the protection scope of the present invention.

Claims

1. A method for calculating earthquake intensity values, characterized in that, The method includes: Acquire seismic event information and site shear wave velocity Vs30 data for the target area; Based on the earthquake event information and the site shear wave velocity Vs30 data, the first PGA distribution of the target area is obtained; the first PGA distribution is the peak ground acceleration (PGA) value. Obtain the DEM data of the target area, and obtain the terrain magnification factor of the target area based on the DEM data and a preset slope correction factor; The first PGA distribution is corrected according to the terrain magnification factor to obtain the second PGA distribution; The site amplification factor of the target area is obtained based on the site shear wave velocity Vs30 data and the preset basin effect correction factor. The second PGA distribution is corrected based on the site magnification factor to obtain the third PGA distribution; Based on the preset correspondence between PGA and seismic intensity, the third PGA distribution is converted into the seismic intensity value of the target area.

2. The method for calculating earthquake intensity as described in claim 1, characterized in that, The first PGA distribution is obtained using the following formula: Among them, PGA base The first PGA distribution is defined as follows: C1, C2, C3, C4, and C5 are preset attenuation coefficients; M is the earthquake magnitude in the earthquake event information; R is the epicenter distance in the earthquake event information; S is the site type parameter of the target area; and V is the ground type parameter of the target area. s30 The data refers to the site shear wave velocity Vs30.

3. The method for calculating earthquake intensity as described in claim 1, characterized in that, The step of obtaining the terrain magnification factor of the target area based on the DEM data and the preset slope correction factor includes: The altitude of the target point in the target area is obtained based on the DEM data; the target point can be any point in the target area. The terrain shape ratio of the target point is obtained based on the elevation of the target point, the horizontal distance component from the target point along the maximum slope direction to the foot of the slope, and the vertical distance component from the target point along the maximum slope direction to the foot of the slope. The terrain magnification factor of the target area is obtained based on the terrain shape ratio and the preset slope correction factor.

4. The method for calculating earthquake intensity as described in claim 3, characterized in that, The terrain magnification factor is obtained using the following formula: THE h =1+0.4×arctan(4SR)×F s Among them, A h Where SR is the terrain magnification factor, F is the terrain shape ratio, and F is the terrain magnification factor. s This is the preset slope correction factor.

5. The method for calculating earthquake intensity as described in claim 3, characterized in that, The second PGA distribution is obtained using the following formula: PGA topo =PGA base ×A h ×(1+0.1×cos(θ-φ)) Among them, PGA topo This is the second PGA distribution, PGA base For the first PGA distribution, A h φ is the terrain magnification factor, θ is the azimuth angle from the epicenter to the target point, and φ is the mountain orientation angle of the target point.

6. The method for calculating earthquake intensity as described in claim 1, characterized in that, After obtaining the seismic intensity value of the target area, the method further includes: Obtain the fracture zone enhancement amplitude and the distance from the target point to the fracture zone in the target area; The seismic intensity value is corrected based on the magnitude of the fault zone enhancement and the distance from the target point to the fault zone to obtain the corrected seismic intensity value.

7. The method for calculating earthquake intensity as described in claim 6, characterized in that, The corrected seismic intensity value is obtained using the following formula: Among them, I fault The corrected seismic intensity value is denoted as I, the seismic intensity value is ΔI, the fault zone enhancement amplitude is d, the distance from the target point to the fault zone is d, and the preset influence range parameter is σ.

8. A system for calculating earthquake intensity values, characterized in that, The system includes: a data acquisition module, a first distribution calculation module, a first coefficient acquisition module, a first correction module, a second coefficient acquisition module, a second correction module, and an intensity value output module; The data acquisition module is used to acquire seismic event information and site shear wave velocity Vs30 data for the target area. The first distribution calculation module is used to obtain the first PGA distribution of the target area based on the seismic event information and the site shear wave velocity Vs30 data; the first PGA distribution is the peak ground acceleration (PGA) value; The first coefficient acquisition module is used to acquire the DEM data of the target area, and obtain the terrain magnification coefficient of the target area based on the DEM data and a preset slope correction factor; The first correction module is used to correct the first PGA distribution according to the terrain magnification factor to obtain the second PGA distribution; The second coefficient acquisition module is used to obtain the site amplification coefficient of the target area based on the site shear wave velocity Vs30 data and the preset basin effect correction coefficient; The second correction module is used to correct the second PGA distribution according to the site magnification factor to obtain a third PGA distribution; The intensity value output module is used to convert the third PGA distribution into the seismic intensity value of the target area according to the preset correspondence between PGA and seismic intensity.

9. A device for calculating earthquake intensity values, comprising a memory and a processor, characterized in that, The memory stores a computer program that, when executed by the processor, causes the processor to perform the steps of the method as described in any one of claims 1 to 7.

10. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by a processor, it causes the processor to perform the steps of the method as described in any one of claims 1 to 7.

Images