Seismic intensity determination method and system based on data correction, storage medium and equipment

Abstract

The embodiment of the invention discloses a seismic intensity determination method based on data correction. The method comprises the following steps: acquiring initial seismic parameters representing a seismic event; calculating seismic oscillation parameters PGA of a plurality of geographic positions in a target area based on a preset seismic oscillation attenuation relation model and the initial seismic parameters; performing terrain amplification correction on the seismic oscillation parameter based on the elevation data of the target area to obtain a first correction parameter PGAtopo; correcting the first correction parameter PGA < topo > based on the field geological condition data of the target area to obtain a seismic oscillation parameter PGA < site >; and based on the seismic oscillation parameter PGAsite and a preset conversion relationship, determining seismic intensity values corresponding to the plurality of geographic positions. According to the invention, a dual correction mechanism of elevation data and site geological condition data is introduced, so that the accuracy of seismic intensity assessment is remarkably improved.

Description

Technical Field

[0001] This invention relates to the field of earthquake intensity monitoring technology, and in particular to earthquake intensity determination methods, systems, storage media and devices based on data correction. Background Technology

[0002] Seismic intensity is a key indicator for assessing the impact of earthquakes on the Earth's surface and engineering structures, and is crucial for post-earthquake emergency response and disaster assessment.

[0003] Traditional earthquake intensity assessment methods mainly rely on post-earthquake field investigations and limited instrument monitoring, which have prominent problems such as poor timeliness, incomplete spatial coverage, and high labor costs, making it difficult to meet the needs of modern earthquake emergency response for rapid and accurate assessment. Summary of the Invention

[0004] Therefore, it is necessary to propose a data-corrected method for determining seismic intensity to address the above problems.

[0005] A method for determining seismic intensity based on data correction, the method comprising the following steps:

[0006] Obtain initial seismic parameters characterizing the seismic event;

[0007] Based on the preset seismic motion attenuation relationship model and the initial seismic parameters, the seismic motion parameters PGA of several geographical locations within the target area are calculated;

[0008] Based on the elevation data of the target area, the seismic motion parameters are subjected to topographic amplification correction to obtain the first correction parameter PGA. topo ;

[0009] Based on the site geological condition data of the target area, the first correction parameter PGA is adjusted. topo Correction is performed to obtain the ground motion parameter PGA. site ;

[0010] Based on the aforementioned ground motion parameter PGA site Based on the preset conversion relationship, the earthquake intensity values ​​corresponding to the several geographical locations are determined.

[0011] In the above scheme, the calculation of the ground motion parameters PGA at several geographical locations within the target area based on the preset ground motion attenuation relationship model and the initial seismic parameters specifically includes:

[0012] The initial earthquake parameters include magnitude and epicentral distance;

[0013] The ground motion parameters PGA at several geographical locations within the target area are determined based on the following ground motion attenuation relationship model:

[0014]

[0015] Where Y is the ground motion parameter PGA, M is the magnitude, R is the epicentral distance, S is the site type parameter, and C1, C2, C3, C4, and C5 are the first regression coefficient, the second regression coefficient, the third regression coefficient, the fourth regression coefficient, and the fifth regression coefficient, respectively.

[0016] In the above scheme, the ground motion parameters are subjected to topographic amplification correction based on the elevation data of the target area to obtain the first correction parameter PGA. topo Specifically, it includes:

[0017] The terrain shape ratio of the target area is determined using the elevation data of the target area. Where H is the terrain height and L is the terrain length;

[0018] The terrain magnification factor A is determined based on the terrain shape ratio of the target area. h And according to the terrain magnification factor A h The ground motion parameters are subjected to topographic amplification correction to obtain the first correction parameter PGA. topo .

[0019] In the above scheme, the terrain units of the target area are divided into mountain tops, ridges and depressions according to the elevation data, and the different terrain units correspond to matching terrain magnification factors.

[0020] In the above scheme, the terrain magnification factor A is determined based on the terrain shape ratio of the target area. h And according to the terrain magnification factor A h The ground motion parameters are subjected to topographic amplification correction to obtain the first correction parameter PGA. topo Specifically, it includes:

[0021] Determine the terrain magnification factor A using the following formula. h :

[0022]

[0023] in, A represents the terrain shape ratio of the target area. h This is the terrain magnification factor;

[0024] The first correction parameter PGA is determined according to the following formula. topo :

[0025] PGA topo =PGA×A h

[0026] Among them, A hPGA is the terrain amplification factor, and PGA is the ground motion parameter. topo This is the first correction parameter.

[0027] In the above scheme, the first correction parameter PGA is adjusted based on the site geological condition data of the target area. topo Correction is performed to obtain the ground motion parameter PGA. site Specifically, it includes:

[0028] The ground motion parameter PGA is determined using the following formula. site :

[0029]

[0030] Among them, PGA site For ground motion parameters, PGA topo V is the first correction parameter. s30 This is data on the site's geological conditions.

[0031] In the above scheme, the step based on the seismic motion parameter PGA site Based on a preset conversion relationship, the seismic intensity values ​​corresponding to the aforementioned geographical locations are determined, specifically including:

[0032]

[0033] Among them, PGA site Here, I represents the seismic motion parameter, and I represents the seismic intensity value.

[0034] This application also proposes a seismic intensity determination system based on data correction, the system comprising: an initial seismic parameter acquisition unit, a ground motion parameter PGA acquisition unit, a correction unit, and a seismic intensity value calculation unit;

[0035] The initial seismic parameter acquisition unit is used to acquire initial seismic parameters characterizing seismic events;

[0036] The ground motion parameter PGA acquisition unit is used to calculate the ground motion parameter PGA of several geographical locations within the target area based on a preset ground motion attenuation relationship model and the initial earthquake parameters.

[0037] The correction unit is used to perform topographic amplification correction on the seismic motion parameters based on the elevation data of the target area to obtain a first correction parameter PGA. topo Based on the site geological condition data of the target area, the first correction parameter PGA is adjusted. topo Correction is performed to obtain the ground motion parameter PGA. site ;

[0038] The earthquake intensity calculation unit is used to calculate the earthquake motion parameter PGA.site Based on the preset conversion relationship, the earthquake intensity values ​​corresponding to the several geographical locations are determined.

[0039] This application also proposes a readable storage medium storing a computer program, which, when executed by a processor, causes the processor to perform the following steps:

[0040] Obtain initial seismic parameters characterizing the seismic event;

[0041] Based on the preset seismic motion attenuation relationship model and the initial seismic parameters, the seismic motion parameters PGA of several geographical locations within the target area are calculated;

[0042] Based on the elevation data of the target area, the seismic motion parameters are subjected to topographic amplification correction to obtain the first correction parameter PGA. topo ;

[0043] Based on the site geological condition data of the target area, the first correction parameter PGA is adjusted. topo Correction is performed to obtain the ground motion parameter PGA. site ;

[0044] Based on the aforementioned ground motion parameter PGA site Based on the preset conversion relationship, the earthquake intensity values ​​corresponding to the several geographical locations are determined.

[0045] This application also proposes a computer 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 causes the processor to perform the following steps:

[0046] Obtain initial seismic parameters characterizing the seismic event;

[0047] Based on the preset seismic motion attenuation relationship model and the initial seismic parameters, the seismic motion parameters PGA of several geographical locations within the target area are calculated;

[0048] Based on the elevation data of the target area, the seismic motion parameters are subjected to topographic amplification correction to obtain the first correction parameter PGA. topo ;

[0049] Based on the site geological condition data of the target area, the first correction parameter PGA is adjusted. topo Correction is performed to obtain the ground motion parameter PGA. site ;

[0050] Based on the aforementioned ground motion parameter PGA site Based on the preset conversion relationship, the earthquake intensity values ​​corresponding to the several geographical locations are determined.

[0051] The embodiments of this invention have the following beneficial effects: By introducing a multi-level data correction mechanism, this invention significantly improves the accuracy and reliability of seismic intensity assessment. First, based on a preset seismic motion attenuation relationship model and initial seismic parameters, the ground motion parameters PGA at various geographical locations within the target area are calculated, providing basic data for subsequent correction. Second, the ground motion parameters PGA are corrected for topographic amplification using elevation data to obtain the first correction parameter, effectively eliminating the influence of topographic undulations on seismic wave propagation, especially under complex terrain conditions such as mountainous areas, significantly improving the accuracy of the calculation results. Third, the first correction parameter is corrected for site amplification effect based on site geological condition data, further correcting the differences in seismic response caused by site factors such as soil type and geological structure, making the corrected ground motion parameters closer to the actual seismic impact. Finally, through a preset conversion relationship between PGA and seismic intensity, the corrected ground motion parameters are converted into seismic intensity values, realizing a scientific mapping from ground motion parameters to intensity values. The entire process, through a dual correction mechanism of topography and site, significantly reduces the assessment bias caused by neglecting topography and site effects in traditional methods, improves the accuracy and applicability of seismic intensity assessment, and provides more reliable technical support for earthquake emergency response and disaster assessment. Attached Figure Description

[0052] To more clearly illustrate the technical solutions in the embodiments of the present invention 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 the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0053] in:

[0054] Figure 1 This is a schematic diagram of a data-corrected method for determining seismic intensity in one embodiment;

[0055] Figure 2 This is a schematic diagram of terrain height and feature length parameters in one embodiment. Detailed Implementation

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

[0057] In the following description, numerous specific details are set forth in order to provide a more thorough understanding of the invention; however, it will be apparent to those skilled in the art that the invention may be practiced without one or more of these details; in other instances, certain technical features well-known in the art have not been described in order to avoid confusion with the invention. It should be understood that the invention can be practiced in different forms and should not be construed as limited to the embodiments set forth herein; rather, these embodiments are provided to make the disclosure thorough and complete and to fully convey the scope of the invention to those skilled in the art.

[0058] The terminology used herein is for the purpose of describing particular embodiments only and is not intended to limit the invention. When used herein, the singular forms “a,” “an,” and “the” are also intended to include the plural forms, unless the context clearly indicates otherwise. The terms “comprising” and / or “including,” when used in this specification, identify the presence of said features, integers, steps, operations, elements, and / or components, but do not exclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and / or groups. When used herein, the term “and / or” includes any and all combinations of the associated listed items.

[0059] To fully understand the present invention, a detailed structure will be presented in the following description in order to illustrate the technical solution proposed by the present invention; optional embodiments of the present invention are described in detail below, however, in addition to these detailed descriptions, the present invention may have other embodiments.

[0060] like Figure 1 As shown, in one embodiment, a seismic intensity determination method based on data correction is provided. This method includes steps S101 to S105, which are detailed below:

[0061] S101. Obtain the initial seismic parameters characterizing the seismic event;

[0062] Specifically, initial earthquake parameters include magnitude and epicentral distance. These can be obtained automatically using web crawling technology to retrieve real-time earthquake information released by seismic networks or by receiving manually input parameters. These parameters are prerequisites for calculating seismic motion parameters and subsequent corrections; their accuracy and timeliness directly affect the reliability and efficiency of the entire assessment process.

[0063] Furthermore, the capture frequency is set as a configurable parameter, with a default value of 5 minutes. System administrators can dynamically adjust it to a range of 1-30 minutes based on emergency response needs. For the captured seismic data, the following three-level filtering verification is performed:

[0064] The first level is data integrity verification, ensuring that the magnitude and epicentral distance parameters are complete and in the correct format;

[0065] The second level is the geographical scope filtering, where the system presets the boundary coordinates of the study area and retains only earthquake events located within the target area;

[0066] The third level is magnitude threshold filtering, when the magnitude M ≥ M threshold (M threshold The subsequent analysis process is triggered when the default setting is 3.0. The threshold parameter can be adjusted according to the seismic activity characteristics and emergency response requirements of different regions.

[0067] S102. Based on the preset seismic motion attenuation relationship model and initial seismic parameters, calculate the seismic motion parameters PGA at several geographical locations within the target area;

[0068] The ground motion attenuation model was used to obtain the ground motion parameters (PGA) at various geographical locations. This model maps macroscopic parameters such as magnitude and distance to microscopic ground motion intensity using scientific formulas, laying a scientific foundation for subsequent correction and intensity determination.

[0069] In some embodiments, based on a preset seismic attenuation relationship model and initial seismic parameters, the ground motion parameters PGA at several geographical locations within the target area are calculated, specifically including:

[0070] Initial earthquake parameters include magnitude and epicentral distance;

[0071] The ground motion parameters PGA at several geographical locations within the target area are determined based on the following ground motion attenuation relationship model:

[0072]

[0073] Where Y is the ground motion parameter PGA, M is the magnitude, R is the epicentral distance, S is the site type parameter, and C1, C2, C3, C4, and C5 are the first regression coefficient, the second regression coefficient, the third regression coefficient, the fourth regression coefficient, and the fifth regression coefficient, respectively.

[0074] Furthermore, the values ​​of regression coefficients C1 to C5 are fitted using a two-step regression method. These regression coefficients have different values ​​in the major and minor axis directions to reflect the anisotropic characteristics of the geological structure.

[0075] Specifically, the first step is to fix the fourth regression coefficient C4 and fit the remaining coefficients using the least squares method; the second step is to optimize the fourth regression coefficient C4 using the grid search method based on the results of the first step, so as to maximize the overall goodness of fit.

[0076] S103. Based on the elevation data of the target area, perform terrain amplification correction on the seismic motion parameters to obtain the first correction parameter PGA. topo ;

[0077] Topography has a significant impact on seismic wave propagation, especially in mountainous areas and other regions with dramatic topographic relief, where seismic waves may be amplified due to the terrain amplification effect. This step, by introducing elevation data and calculating the terrain amplification factor, corrects for the PGA (Programme for Gaining Apparent Radical Aggregate), significantly improving the accuracy of the assessment results under complex terrain conditions. This correction mechanism effectively compensates for the shortcomings of traditional methods that neglect topographic factors, making the results closer to the actual impact of earthquakes.

[0078] In some embodiments, the ground motion parameters are subjected to topographic amplification correction based on the elevation data of the target area to obtain the first correction parameter PGA. topo Specifically, it includes:

[0079] Determine the terrain shape ratio of the target area using elevation data of the target area. Where H is the terrain height and L is the terrain length;

[0080] Determine the terrain magnification factor A based on the terrain shape ratio of the target area. h And based on the terrain magnification factor A h The ground motion parameters are corrected for topographic amplification to obtain the first correction parameter PGA. topo .

[0081] Specifically, based on DEM data (elevation data), the terrain height H and characteristic length L of each grid point in the study area are calculated, and the terrain shape ratio is determined. The calculation range is the ratio of the vertical distance from the top of the mountain to the horizontal distance from the bottom of the mountain.

[0082] Among them, the terrain magnification factor A h The value ranges from 1.1 to 2.1, with the magnification effect being most significant in the mountain top and ridge areas. For the concave terrain and foot areas, the magnification factor is set to 1.0 or a random number between 0.99 and 1.0.

[0083] In some embodiments, the terrain units of the target area are divided into mountain peaks, ridges, and depressions based on elevation data, and different terrain units correspond to matching terrain magnification factors.

[0084] like Figure 2 The diagram shows the terrain height and feature length parameters. In some embodiments, the terrain parameters are calculated based on 30-meter resolution DEM data, and the terrain features are extracted using a moving window algorithm.

[0085] For example, to calculate the terrain height H, the system sets a search window with a radius of 500 meters around each grid point, and calculates the difference between the maximum and minimum elevation values ​​within the window: H = h max -h minThe terrain length L is calculated using the slope change rate method. By analyzing the gradient field of the DEM data, the characteristic length of the terrain undulations is identified. in, The Laplace operator for elevation.

[0086] Furthermore, the terrain classification algorithm divides the study area into different types based on terrain shape ratio and local slope:

[0087] Mountain top area: And the local slope is less than 5°;

[0088] Ridge area: And the local slope is ≥5°;

[0089] Hillside area:

[0090] Plains area:

[0091] Depression area: Identified through terrain curvature analysis, curvature > 0.01m -1 .

[0092] Furthermore, the calculation strategy for the terrain magnification factor differs for different terrain types:

[0093] Mountain tops and ridges: The value range is 1.3-2.1;

[0094] Hillside area: The value range is 1.1-1.5;

[0095] Plains area: A h =1.0;

[0096] Depressed area: A h =0.95+0.05×random, introducing random perturbations to simulate complex terrain effects.

[0097] In some embodiments, the terrain magnification factor A is determined based on the terrain shape ratio of the target area. h And based on the terrain magnification factor A h The ground motion parameters are corrected for topographic amplification to obtain the first correction parameter PGA. topo Specifically, it includes:

[0098] Determine the terrain magnification factor A according to the following formula. h :

[0099]

[0100] in, A represents the terrain shape ratio of the target area. hThis is the terrain magnification factor;

[0101] The first correction parameter PGA is determined according to the following formula. topo :

[0102] PGA topo =PGA×A h

[0103] Among them, A h PGA is the terrain amplification factor, and PGA is the ground motion parameter. topo This is the first correction parameter.

[0104] S104. Adjusting the first correction parameter PGA based on the site geological condition data of the target area. topo Correction is performed to obtain the ground motion parameter PGA. site ;

[0105] Site geological conditions (such as soil type and shear wave velocity) have a significant impact on seismic response, and the intensity of seismic ground motion varies considerably under different site conditions. This step further improves the accuracy of the assessment by incorporating site geological data to perform a secondary correction on the topographically corrected PGA. This site amplification effect correction mechanism allows the final seismic ground motion parameters to more accurately reflect the actual dynamic response of the surface and subsurface structures.

[0106] In some embodiments, the first correction parameter PGA is adjusted based on site geological condition data of the target area. topo Correction is performed to obtain the ground motion parameter PGA. site Specifically, it includes:

[0107] The ground motion parameter PGA is determined using the following formula. site :

[0108]

[0109] Among them, PGA site For ground motion parameters, PGA topo V is the first correction parameter. s30 This is data on the site's geological conditions.

[0110] Specifically, Vs30 data for each grid point in the study area were obtained. The Vs30 data ranged from 180 m / s to 760 m / s, with 760 m / s used as the reference value for bedrock velocity structure.

[0111] In some embodiments, the site correction calculation takes into account the nonlinear response characteristics of the PGA value, and the complete site correction formula is:

[0112] PGA site =PGA topo ×exp(FLIN +F NL )

[0113] The linear part includes: Nonlinear component (considered when PGA > 0.1g): PGA site For ground motion parameters, PGA topo Vs30 is the first correction parameter, and Vs30 is the site geological condition data.

[0114] S105, Based on seismic motion parameters PGA site Based on the preset conversion relationship, determine the earthquake intensity values ​​corresponding to several geographical locations.

[0115] Seismic intensity is a macroscopic indicator describing the degree of impact of an earthquake and has a corresponding relationship with seismic motion parameters. This step maps the corrected PGA value to a seismic intensity value through a pre-defined conversion relationship, achieving a scientific conversion from physical quantity to the degree of disaster impact. This conversion relationship is usually based on historical earthquake data and engineering experience, and has broad applicability and authority, ensuring the scientific validity and practicality of the intensity assessment results.

[0116] In some embodiments, based on ground motion parameters PGA site Based on the preset conversion relationship, determine the seismic intensity values ​​corresponding to several geographical locations, specifically including:

[0117]

[0118] Among them, PGA site Here, I represents the seismic motion parameter, and I represents the seismic intensity value.

[0119] Preferably, the corrected PGA value (Geostationary Ground Motion Parameter PGA) is used. site The simulation is converted into seismic intensity values ​​and a spatial distribution map of intensity is generated. Based on the spatial distribution map of intensity, an intensity contour map and analysis report are generated, supporting PDF and Word format output. Considering the influence of active fault distribution on the spatial distribution of intensity, accuracy requirements are imposed on the generated intensity distribution results, and the deviation between the simulated intensity and the intensity monitored by the instrument does not exceed 0.5 degrees.

[0120] Specifically, the influence of active faults is considered using a distance attenuation function, which amplifies the intensity values ​​near the fault zone (distance less than 5 km):

[0121]

[0122] Where, d fault f is the distance to the nearest active fault. activity This is the fracture activity coefficient.

[0123] To ensure the spatial continuity of intensity distribution, the Laplace smoothing algorithm can be used to optimize the initial calculation results:

[0124]

[0125] Where α = 0.1 is the smoothing coefficient;

[0126] Furthermore, the rationality of the calculation results was verified by statistical comparison with historical earthquake intensity distributions.

[0127] For an earthquake of magnitude M, its intensity attenuation should satisfy the empirical relationship:

[0128] I(R) = I0 - 3.5log 10 (R)-0.006R

[0129] Where I0 = 1.5M - 1.4 is the empirical estimate of the epicenter intensity, and R is the epicenter distance (km); when the calculated result deviates from the empirical value by more than 0.8 degrees, it is automatically marked as an outlier and a manual review process is triggered.

[0130] In some embodiments, a cartographic knowledge rule base is constructed, and map templates of different types and scales are pre-made; satellite imagery, topographic maps, and annotated maps are automatically loaded as base maps; multi-layer overlay display is output, including intensity contour lines, active faults, historical earthquakes, and infrastructure distribution; and WebGIS functions are provided, the WebGIS functions including at least the following interactive operations: spatial query, distance measurement, and area calculation.

[0131] To achieve fast response, a multi-level caching and parallel computing optimization strategy is adopted. This strategy includes:

[0132] Pre-calculate and store the base attenuation matrix under different magnitude ranges within the study area, and cache commonly used terrain correction coefficients and site correction parameters to reduce real-time computational load;

[0133] The study area was divided into multiple computing units according to a geographic grid, and a MapReduce distributed computing model was adopted, with each computing node processing the PGA calculation and correction process in different areas in parallel:

[0134] T total =max(T) i )+T m

[0135] Among them, T i Let T be the processing time of the i-th computing node. m For the time of result merging;

[0136] When a new seismic event occurs, the radius of influence R is recalculated only for grid points within the affected area. i (Unit: km) Determined based on magnitude M: R i =10 0.5M+1.5 ;

[0137] By using memory mapping technology to directly access DEM and Vs30 data, redundant disk I / O operations are avoided, and the data access time is optimized from O(n) to O(1).

[0138] In summary, this invention achieves high accuracy and reliability in seismic intensity assessment through a multi-level data correction mechanism. First, a preliminary PGA calculation is performed based on a ground motion attenuation model, providing fundamental data for subsequent corrections. Second, terrain amplification correction effectively eliminates the influence of terrain undulations on seismic wave propagation, significantly improving assessment accuracy, especially under complex terrain conditions. Third, site amplification effect correction further corrects ground motion response deviations caused by differences in site conditions, making the results closer to the actual earthquake impact. Finally, through a scientific conversion relationship, the corrected PGA is mapped to seismic intensity values, achieving a reasonable mapping from ground motion parameters to intensity values.

[0139] The entire process from S101 to S105, through a dual correction mechanism of topography and site, significantly reduces the assessment bias caused by neglecting topography and site effects in traditional methods, improves the accuracy, applicability, and stability of seismic intensity assessment, and provides more reliable technical support for earthquake emergency response, disaster assessment, and earthquake prevention and mitigation.

[0140] This application also proposes a seismic intensity determination system based on data correction, which includes: an initial seismic parameter acquisition unit, a ground motion parameter PGA acquisition unit, a correction unit, and a seismic intensity value calculation unit;

[0141] The initial seismic parameter acquisition unit is used to acquire the initial seismic parameters characterizing the seismic event;

[0142] The ground motion parameter PGA acquisition unit is used to calculate the ground motion parameter PGA of several geographical locations within the target area based on a preset ground motion attenuation relationship model and initial earthquake parameters.

[0143] The correction unit is used to perform topographic amplification correction on the seismic motion parameters based on the elevation data of the target area to obtain the first correction parameter PGA. topo The first correction parameter PGA is adjusted based on the site geological condition data of the target area. topo Correction is performed to obtain the ground motion parameter PGA. site ;

[0144] Earthquake intensity value calculation unit, used for calculating seismic motion parameters PGA siteBased on the preset conversion relationship, determine the earthquake intensity values ​​corresponding to several geographical locations.

[0145] This application also proposes a readable storage medium storing a computer program, which, when executed by a processor, causes the processor to perform the following steps:

[0146] Obtain initial seismic parameters characterizing the seismic event;

[0147] Based on the preset seismic motion attenuation relationship model and the initial seismic parameters, the seismic motion parameters PGA of several geographical locations within the target area are calculated;

[0148] Based on the elevation data of the target area, the seismic motion parameters are subjected to topographic amplification correction to obtain the first correction parameter PGA. topo ;

[0149] Based on the site geological condition data of the target area, the first correction parameter PGA is adjusted. topo Correction is performed to obtain the ground motion parameter PGA. site ;

[0150] Based on the aforementioned ground motion parameter PGA site Based on the preset conversion relationship, the earthquake intensity values ​​corresponding to the several geographical locations are determined.

[0151] This application also proposes a computer 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 causes the processor to perform the following steps:

[0152] Obtain initial seismic parameters characterizing the seismic event;

[0153] Based on the preset seismic motion attenuation relationship model and the initial seismic parameters, the ground motion parameters PGA of several geographical locations within the target area are calculated;

[0154] Based on the elevation data of the target area, the seismic motion parameters are subjected to topographic amplification correction to obtain the first correction parameter PGA. topo ;

[0155] Based on the site geological condition data of the target area, the first correction parameter PGA is adjusted. topo Correction is performed to obtain the ground motion parameter PGA. site ;

[0156] Based on the aforementioned ground motion parameter PGA site Based on the preset conversion relationship, the earthquake intensity values ​​corresponding to the several geographical locations are determined.

[0157] Those skilled in the art will understand that implementing all or part of the processes in the above embodiments can be accomplished by instructing related hardware through a computer program. The program can be stored in a non-volatile computer-readable storage medium, and when executed, it can include the processes of the embodiments of the above methods. Any references to memory, storage, databases, or other media used in the embodiments provided in this application can include non-volatile and / or volatile memory. Non-volatile 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 various forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), dual data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), Rambus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.

[0158] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0159] The embodiments described above are merely illustrative of 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 application's patent. 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. The embodiments disclosed above are merely preferred embodiments of the present invention and should not be construed as limiting the scope of the present invention. Therefore, equivalent variations made according to the claims of this invention are still within the scope of this invention.

Claims

1. A method for determining seismic intensity based on data correction, characterized in that, The method includes: Obtain initial seismic parameters characterizing the seismic event; Based on the preset seismic motion attenuation relationship model and the initial seismic parameters, the seismic motion parameters PGA of several geographical locations within the target area are calculated; Based on the elevation data of the target area, the seismic motion parameters are subjected to topographic amplification correction to obtain the first correction parameter PGA. topo ; Based on the site geological condition data of the target area, the first correction parameter PGA is adjusted. topo Correction is performed to obtain the ground motion parameter PGA. site ; Based on the aforementioned ground motion parameter PGA site Based on the preset conversion relationship, the earthquake intensity values ​​corresponding to the several geographical locations are determined.

2. The seismic intensity determination method based on data correction according to claim 1, characterized in that, The calculation of ground motion parameters (PGA) at several geographical locations within the target area, based on a preset ground motion attenuation relationship model and the initial seismic parameters, specifically includes: The initial earthquake parameters include magnitude and epicentral distance; The ground motion parameters PGA at several geographical locations within the target area are determined based on the following ground motion attenuation relationship model: Where Y is the ground motion parameter PGA, M is the magnitude, R is the epicentral distance, S is the site type parameter, and C1, C2, C3, C4, and C5 are the first regression coefficient, the second regression coefficient, the third regression coefficient, the fourth regression coefficient, and the fifth regression coefficient, respectively.

3. The seismic intensity determination method based on data correction according to claim 1, characterized in that, The ground motion parameters are corrected by topographic amplification based on the elevation data of the target area to obtain the first correction parameter PGA. topo Specifically, it includes: The terrain shape ratio of the target area is determined using the elevation data of the target area. Where H is the terrain height and L is the terrain length; The terrain magnification factor A is determined based on the terrain shape ratio of the target area. h And according to the terrain magnification factor A h The ground motion parameters are subjected to topographic amplification correction to obtain the first correction parameter PGA. topo .

4. The seismic intensity determination method based on data correction according to claim 3, characterized in that, Based on the elevation data, the terrain units of the target area are divided into mountain peaks, ridges, and depressions, and the different terrain units correspond to matching terrain magnification factors.

5. The seismic intensity determination method based on data correction according to claim 3, characterized in that, The terrain magnification factor A is determined based on the terrain shape ratio of the target area. h And according to the terrain magnification factor A h The ground motion parameters are subjected to topographic amplification correction to obtain the first correction parameter PGA. topo Specifically, it includes: Determine the terrain magnification factor A using the following formula. h : in, A represents the terrain shape ratio of the target area. h This is the terrain magnification factor; The first correction parameter PGA is determined according to the following formula. topo : PGA topo =PGA×A h Among them, A h PGA is the terrain amplification factor, and PGA is the ground motion parameter. topo This is the first correction parameter.

6. The seismic intensity determination method based on data correction according to claim 5, characterized in that, The first correction parameter PGA is adjusted based on the site geological condition data of the target area. topo Correction is performed to obtain the ground motion parameter PGA. site Specifically, it includes: The ground motion parameter PGA is determined using the following formula. site : Among them, PGA site For ground motion parameters, PGA topo V is the first correction parameter. s30 This refers to the site's geological conditions data.

7. The seismic intensity determination method based on data correction according to claim 6, characterized in that, The seismic motion parameter PGA site Based on a preset conversion relationship, the seismic intensity values ​​corresponding to the aforementioned geographical locations are determined, specifically including: Among them, PGA site Here, I represents the seismic motion parameter, and I represents the seismic intensity value.

8. A seismic intensity determination system based on data correction, characterized in that, The system includes: an initial seismic parameter acquisition unit, a ground motion parameter PGA acquisition unit, a correction unit, and a seismic intensity value calculation unit; The initial seismic parameter acquisition unit is used to acquire initial seismic parameters characterizing seismic events; The ground motion parameter PGA acquisition unit is used to calculate the ground motion parameter PGA of several geographical locations within the target area based on a preset ground motion attenuation relationship model and the initial earthquake parameters. The correction unit is used to perform topographic amplification correction on the seismic motion parameters based on the elevation data of the target area to obtain a first correction parameter PGA. topo Based on the site geological condition data of the target area, the first correction parameter PGA is adjusted. topo Correction is performed to obtain the ground motion parameter PGA. site ; The earthquake intensity calculation unit is used to calculate the earthquake motion parameter PGA. site Based on the preset conversion relationship, the earthquake intensity values ​​corresponding to the several geographical locations are determined.

9. A readable storage medium storing a computer program that, when executed by a processor, causes the processor to perform the steps of the method as claimed in any one of claims 1 to 7.

10. A computer device comprising a memory and a processor, the memory storing 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.

Images