Method and device for determining seismic theoretical intensity of attention area

Abstract

The embodiment of the specification provides a method and device for judging seismic theoretical intensity of a concerned area, wherein the method comprises the following steps: deducing long and short axis sizes of different intensity values under different magnitudes in the concerned area through a seismic influence field model of the concerned area; determining a seismic influence field direction of the concerned area; determining an elliptical model according to the long and short axis sizes and the seismic influence field direction, and obtaining new coordinates (X', Y') after rotating the elliptical model; bringing the elliptical model into the converted new coordinates (X', Y'), calculating a final value, comparing the final value with 1, and determining an intensity circle where the key concerned area is located according to a comparison result.

Description

Technical Field

[0001] This document relates to the field of computer technology, and in particular to a method and apparatus for determining the theoretical intensity of earthquakes in a region of interest. Background Technology

[0002] The distribution of earthquake energy near the epicenter, affecting surrounding buildings and facilities to varying degrees, is called the earthquake impact field. Its distribution is influenced by many objective factors, such as the earthquake's source parameters (including magnitude, location, and depth) and the effects of the propagation medium (e.g., geological structures and topography around the epicenter). The rapid determination method for the earthquake impact field is the theoretical basis for initially assessing the affected area after an earthquake. In the immediate aftermath of an earthquake, before information about the situation at the scene can be reported or collected, the rapid determination of the earthquake impact field provides government departments with the initial theoretical basis and guidance for earthquake relief efforts.

[0003] When an earthquake occurs, earthquake officials typically input the three key earthquake parameters into a local seismic intensity attenuation model immediately. By calculating the theoretical major and minor axis radii for different intensity values, they then use GIS tools to plot the earthquake impact field model onto a thematic base map of earthquake elements, based on the epicenter location and the major and minor axis radii. The earthquake emergency response commander or other personnel then use this thematic map information to advance subsequent earthquake emergency response efforts. Figure 1 This is a schematic diagram of the impact field assessment of the 2014 Ludian 6.5 earthquake in existing technology.

[0004] The seismic intensity attenuation model is: I = A + BM + Clg(R + R0), where I is the seismic intensity; M is the surface wave magnitude; R is the epicentral distance; and A, B, C, and R0 are regression coefficients. Generally, an elliptical model is used to ensure that the intensity is equal along the major and minor axes when R = 0, while maintaining the difference in intensity along the major and minor axes at intermediate distances, and also to make the isoseismal lines circular in the far field. To reflect the disappearance of the seismogenic tectonic influence in the far field and the tendency of the attenuation shape to become circular, the radius of the felt area is taken as the far-field control point, with a felt intensity value of III-IV.

[0005] When regressing the attenuation relationship, an elliptical major and minor axis joint attenuation model was adopted, and statistical regression was performed using the least squares method to obtain the attenuation relationship of seismic intensity in various regions of my country, as detailed in Table 1.

[0006] Table 1. Relationship coefficients between earthquake intensity attenuation

[0007]

[0008] In order to ensure that the impact of earthquake damage is close to the local characteristics, each region also determined the earthquake intensity attenuation relationship by measuring the length and short axis data of the isoseismal lines of each earthquake and using the above-mentioned least squares regression method, according to the corresponding earthquake parameters of different intensity values ​​V-VIII in the region.

[0009] When a perceptible earthquake occurs, a common method is to determine the size of the major and minor axes of the elliptical attenuation model based on the parameters of the elliptical attenuation model for different regions or the earthquake intensity attenuation relationship in Table 1 for all regions of the country. Then, GIS is used to draw the earthquake influence field according to the direction of the active fault. Finally, based on the geographical location of the key areas of concern (target buildings), GIS is used to determine whether they are located in the corresponding intensity zone.

[0010] In summary, while the commonly used method is relatively intuitive, it is complex, requires a high level of professional knowledge, and necessitates specialized GIS processing software such as ArcGIS or QGIS, as well as various basic geographic data (administrative divisions). This limits its widespread adoption. Therefore, there is an urgent need for a simpler, more user-friendly method for calculating earthquake intensity that does not require advanced professional knowledge. Summary of the Invention

[0011] The purpose of this invention is to provide a method and apparatus for determining the theoretical intensity of earthquakes in a region of interest, thereby addressing the aforementioned problems in the prior art.

[0012] This invention provides a method for determining the theoretical intensity of earthquakes in a region of interest, comprising:

[0013] By using the earthquake influence field model of the area of ​​interest, the magnitudes of the major and minor axes of different intensity values ​​under different magnitudes in the area of ​​interest are derived.

[0014] Determine the direction of the seismic influence field in the area of ​​interest;

[0015] The elliptical model is determined based on the magnitude of the major and minor axes and the direction of the seismic influence field, and the elliptical model is rotated to obtain the new coordinates (X', Y').

[0016] Substitute the elliptical model into the new coordinate system (X', Y') after transformation, calculate the final value, compare it with 1, and determine the intensity zone of the key concern area based on the comparison result.

[0017] This invention provides a device for determining the theoretical intensity of earthquakes in a region of interest, comprising:

[0018] The derivation module is used to derive the major and minor axes of different intensity values ​​under different magnitudes in the area of ​​interest from the earthquake influence field model of the area of ​​interest;

[0019] The determination module is used to determine the direction of the seismic influence field in the area of ​​interest;

[0020] The ellipse model module is used to determine the ellipse model based on the size of the major and minor axes and the direction of the seismic influence field, and to rotate the ellipse model to obtain a new coordinate (X', Y').

[0021] The judgment module is used to input the elliptical model into the new coordinates (X', Y') after transformation, calculate the final value and compare it with 1, and determine the intensity circle of the key concern area based on the comparison result.

[0022] This invention also provides an electronic device, including: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the computer program, when executed by the processor, implements the steps of the above-mentioned method for determining the theoretical seismic intensity of the area of ​​interest.

[0023] This invention also provides a computer-readable storage medium storing an information transmission implementation program, which, when executed by a processor, implements the steps of the above-mentioned method for determining the theoretical seismic intensity of the area of ​​interest.

[0024] Using the embodiments of the present invention, the theoretical intensity value of the earthquake impact on key areas (landmarks) can be determined without the need for GIS software processing. It has the advantage of not being limited by professional software and being able to quickly determine the earthquake intensity. Attached Figure Description

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

[0026] Figure 1 This is a schematic diagram of the impact field assessment of the 2014 Ludian 6.5 earthquake in existing technology;

[0027] Figure 2 This is a flowchart of the method for determining the theoretical seismic intensity of the area of ​​interest according to an embodiment of the present invention;

[0028] Figure 3 This is a schematic diagram of the earthquake theoretical intensity determination device for the area of ​​interest according to an embodiment of the present invention;

[0029] Figure 4 This is a schematic diagram of an electronic device according to an embodiment of the present invention. Detailed Implementation

[0030] To enable those skilled in the art to better understand the technical solutions in one or more embodiments of this specification, the technical solutions in one or more embodiments of this specification will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of this specification, and not all of the embodiments. Based on one or more embodiments of this specification, all other embodiments obtained by those skilled in the art without creative effort should fall within the protection scope of this document.

[0031] Method Implementation Examples

[0032] According to an embodiment of the present invention, a method for determining the theoretical intensity of earthquakes in a region of interest is provided. Figure 2 This is a flowchart of the method for determining the theoretical seismic intensity of the area of ​​interest according to an embodiment of the present invention, such as... Figure 2 As shown, the method for determining the theoretical seismic intensity of the area of ​​interest according to an embodiment of the present invention specifically includes:

[0033] Step S201 involves deriving the major and minor axes of different intensity values ​​at different magnitudes in the area of ​​interest using the earthquake influence field model of the area of ​​interest; specifically including:

[0034] By substituting magnitudes from 4.0 to 8.0 and intensity from 5 to 8 degrees into the seismic influence field model of the area of ​​interest, the lengths of the major and minor semi-axes of the elliptical attenuation model under different magnitudes and influence intensities can be derived.

[0035] Step S202 involves determining the direction of the seismic influence field in the area of ​​interest; specifically including:

[0036] The main active faults in the area of ​​interest are decomposed into polyline objects, and the center value of each segment is taken as the latitude and longitude coordinates of the fault. The attribute information is then assigned.

[0037] The data from the epicenter and the fault point are calculated and compared. The distance between the two points is calculated, and the fault with the minimum distance is taken as the fault with the highest probability of inducing an earthquake. The direction of the earthquake influence field is determined according to the direction of the rupture of the fault.

[0038] Step S203 involves determining an elliptical model based on the magnitudes of its major and minor axes and the direction of the seismic influence field, and then rotating the elliptical model to obtain new coordinates (X', Y'); specifically including:

[0039] Set the latitude and longitude coordinates of the earthquake epicenter as (Lon 震中 Lat 震中 The coordinates of the region of interest are (Lon 关注 Lat 关注 The magnitude is M, the distance of the area of ​​concern from the epicenter is D, the intensity is N, and the semi-major axis of intensity N is R.N长 The minor half-axis is R N短 The fault strike is T, and it is assumed that the distance of 1 degree of longitude or 1 degree of latitude is 111 kilometers;

[0040] The standard ellipse model is obtained according to formula 1-3:

[0041] The semi-major axis radian ρ of intensity N 长 =R N长 / 111;Formula 1;

[0042] minor semi-axis radian ρ 短 =R N短 / 111;Formula 2;

[0043] Fault strike converted to radian θ 断层 =π·T / 180; Formula 3;

[0044] The standard elliptical model is translated and rotated in polar coordinates, with the epicenter as the center and the fault strike as the direction of rotation. The new coordinates after rotation are (X', Y'):

[0045] X′=cos(θ 断层 )×(Lat 关注 -Lat 震中 )+din(θ 断层 )×(Lon 关注 -Lon 震中 );

[0046] Y′=-sin(θ 断层 )×(Lat 关注 -Lat 震中 )+cos(θ 断层 )×(Lon 关注 -Lon 震中 ).

[0047] Step S204 involves substituting the elliptical model into the transformed new coordinate system (X', Y'), calculating the final value, comparing it with 1, and determining the intensity zone of the key concern area based on the comparison result. Specifically, this includes:

[0048] Step 1, when Then it is determined that the area of ​​concern is located within the intensity zone, and its intensity value is N degrees;

[0049] Step 2, when Then it is determined that the area of ​​concern is located on the boundary of the intensity circle, and its intensity value is N degrees;

[0050] Step 3, when If the area of ​​concern is not within the intensity zone, the intensity value is reduced to N-1, and steps 1-3 are repeated for the same determination.

[0051] The overall solution of this invention is to use a series of mathematical models to solve for a method to quickly determine the theoretical intensity value of the earthquake impact on key areas (or landmarks), and the specific processing is as follows;

[0052] 1. Determine the major and minor axes.

[0053] Based on the earthquake influence field model for this region, the major and minor axes of different intensities at different magnitudes are derived. For example, the earthquake influence field model for Hubei Province is as follows:

[0054] Major axis direction: I = 4.924 + 1.083M - 1.337ln(Ra + 12.0)

[0055] Minor axis direction: I = 3.936 + 1.024M - 1.133ln(Rb + 6.0)

[0056] In the formula, I represents the earthquake intensity; M represents the surface wave magnitude; and R represents the epicentral distance (Ra is the major semi-axis; Rb is the minor semi-axis). By substituting earthquake magnitudes from 4.0 to 8.0 and intensity values ​​from V to VIII (5 degrees to 8 degrees) into the earthquake influence field model, the lengths of the major and minor semi-axis of the elliptical attenuation model under different magnitudes and influence intensities can be derived. Table 2 shows the major and minor axis radii of the earthquake influence field in Hubei Province.

[0057] Table 2

[0058]

[0059]

[0060] 2. Determination of the direction of the earthquake influence field

[0061] The direction of the earthquake's impact field is generally determined based on the strike of active faults near the epicenter. For example, the major active faults in Hubei Province can be decomposed into broken line objects, and the center value of each segment can be taken as the latitude and longitude coordinates of the fault. Attribute information can then be assigned values, as shown in Table 3. Table 3 contains fault data, with the Qianjiang Fault as an example.

[0062] Table 3

[0063]

[0064]

[0065] By comparing data from the epicenter and fault points, the distance between the two points is calculated, and the fault with the minimum distance is considered to have a higher probability of inducing earthquakes. The direction of its seismic influence field is defined as the direction of the fault rupture (i.e., strike).

[0066] 3. Reset the ellipsoid model

[0067] The latitude and longitude coordinates of the earthquake epicenter can be set as (Lon 震中 Lat 震中 The coordinates of the area (or feature) of interest are (Lon 关注 Lat 关注 The magnitude is M, the distance of the area of ​​concern from the epicenter is D, the intensity is N, and the semi-major axis of intensity N is R. N长 The minor half-axis is R N短 The fault strikes T, and it is assumed that the distance of 1 degree of longitude or 1 degree of latitude is 111 kilometers.

[0068] The semi-major axis radian ρ of intensity N 长 =R N长 / 111;

[0069] minor semi-axis radian ρ 短 =R N短 / 111;

[0070] Fault strike converted to radian θ 断层 =π·T / 180;

[0071] A standard elliptical model is translated and rotated in polar coordinates, with the epicenter as the center and the fault strike as the direction of rotation. The new coordinates after rotation are (X', Y'):

[0072] X′=cos(θ 断层 )×(Lat 关注 -Lat 震中 )+sin(θ 断层 )×(Lon 关注 -Lon 震中 )

[0073] Y′=-sin(θ 断层 )×(Lat 关注 -Lat 震中 )+cos(θ 断层 )×(Lon 关注 -Lon 震中 )

[0074] 4. Determine the theoretical intensity value of the area of ​​concern.

[0075] By using the new elliptical model formula generated after translation and rotation, and substituting the transformed new coordinates (X', Y'), the final value is calculated and compared with 1 to determine the intensity circle of the key area of ​​concern.

[0076] when Then it is determined that the area of ​​concern is located within the intensity zone, and its intensity is N degrees;

[0077] when Then it is determined that the area of ​​concern is located on the boundary of the intensity circle, and its intensity value is N degrees;

[0078] when If the area of ​​concern is not within the intensity zone, the intensity can be reduced to N-1, and the above steps can be repeated.

[0079] In summary:

[0080] 1. According to the mathematical formula of the earthquake influence field model, the major and minor axis values ​​of the influence range of the standard ellipse model under different intensities N are derived in this embodiment of the invention, which can be directly used as a two-dimensional array and stored in the database.

[0081] 2. In this embodiment of the invention, a fault data table is created to record the fault strike and location. By calculating the distance from the epicenter, the nearest fault is determined to be the most likely to have triggered the earthquake, that is, the fault strike is determined to be the rotation direction of the elliptical model.

[0082] 3. This invention, through rapid numerical calculation, can determine the theoretical intensity value of the earthquake impact on the area of ​​interest or key landmarks.

[0083] The technical solution of this invention eliminates the need for specialized GIS processing software, reducing the step of having professionals draw earthquake impact distribution maps. Furthermore, mathematical calculations can be used to determine the major and minor axes of each intensity value in the standard theoretical earthquake impact field ellipse attenuation formula, storing this information in a two-dimensional array in a database. When an earthquake occurs, the theoretical standard ellipse model is translated in the local earthquake impact field attenuation model's ellipsoid formula, rotating around the epicenter with the nearest fault strike as the rotation direction. The intensity of the area of ​​interest is then determined in polar coordinates. This invention is computationally simple and requires minimal computer equipment, including memory, processor, and computer application skills. For areas of interest (such as key water conservancy projects like dams, the Wuhan Greenland Center, large stadiums, and schools), their latitude and longitude are determined, allowing enterprises and institutions to independently determine the theoretical intensity values, providing a scientific basis for earthquake resistance assessment of key construction projects.

[0084] For example, if a magnitude 5.5 earthquake occurs with its epicenter located at (30.95°N, 110.76°E) and a focal depth of 10km, the distance from the epicenter to a key water conservancy project dam is approximately 57km. Calculations determine the strike of the nearest fault to be 40°. Using a seismic influence field model for the Hubei region, the major semi-axis of the ellipse model is 69.3km and the minor semi-axis is 50.4km. Substituting the latitude and longitude of the key water conservancy project dam (provided by Baidu) (30.766°N, 111.259°E) into the new ellipse formula yields a result of 1.298. This indicates that the key water conservancy project dam is not within circle V of the theoretical intensity of this earthquake, and its intensity value is below V.

[0085] If the nearest fault strike is determined to be 330° using the distance formula, then substituting the latitude and longitude coordinates of a key water conservancy project dam into the new ellipse formula will yield a result of 0.990. This indicates that the theoretical intensity of the dam under the influence of this earthquake is V.

[0086] Device Example 1

[0087] According to an embodiment of the present invention, a device for determining the theoretical intensity of earthquakes in a region of interest is provided. Figure 3 This is a schematic diagram of the earthquake theoretical intensity determination device for the area of ​​interest according to an embodiment of the present invention, as shown below. Figure 3 As shown, the earthquake theoretical intensity determination device for the area of ​​interest according to an embodiment of the present invention specifically includes:

[0088] Derivation module 30 is used to derive the major and minor axes of different intensity values ​​at different magnitudes in the area of ​​interest from the seismic influence field model of the area of ​​interest; specifically, it is used for:

[0089] By substituting magnitudes from 4.0 to 8.0 and intensity from 5 to 8 degrees into the seismic influence field model of the area of ​​interest, the lengths of the major and minor semi-axes of the elliptical attenuation model under different magnitudes and influence intensities can be derived.

[0090] Module 32 is used to determine the direction of the seismic influence field in the area of ​​interest; specifically, it is used for:

[0091] The main active faults in the area of ​​interest are decomposed into polyline objects, and the center value of each segment is taken as the latitude and longitude coordinates of the fault. The attribute information is then assigned.

[0092] The data from the epicenter and the fault point are calculated and compared. The distance between the two points is calculated, and the fault with the minimum distance is taken as the fault with the highest probability of inducing an earthquake. The direction of the earthquake influence field is determined according to the direction of the rupture of the fault.

[0093] Ellipse model module 34 is used to determine an ellipse model based on the magnitude of the major and minor axes and the direction of the seismic influence field, and to rotate the ellipse model to obtain a new coordinate system (X', Y'); specifically used for:

[0094] Set the latitude and longitude coordinates of the earthquake epicenter as (Lon 震中 Lat 震中 The coordinates of the region of interest are (Lon 关注 Lat 关注 The magnitude is M, the distance of the area of ​​concern from the epicenter is D, the intensity is N, and the semi-major axis of intensity N is R. N长 The minor half-axis is R N短 The fault strike is T, and it is assumed that the distance of 1 degree of longitude or 1 degree of latitude is 111 kilometers;

[0095] The standard ellipse model is obtained according to formula 1-3:

[0096] The semi-major axis radian ρ of intensity N 长 =R N长 / 111;Formula 1;

[0097] minor semi-axis radian ρ 短 =R N短 / 111;Formula 2;

[0098] Fault strike converted to radian θ 断层 =π·T / 180; Formula 3;

[0099] The standard elliptical model is translated and rotated in polar coordinates, with the epicenter as the center and the fault strike as the direction of rotation. The new coordinates after rotation are (X', Y'):

[0100] X′=cos(θ 断层 )×(Lat 关注 -Lat 震中 )+sin(θ 断层 )×(Lon 关注 -Lon 震中 );

[0101] Y′=-sin(θ 断层 )×(Lat 关注 -Lat 震中 )+cos(θ 断层 )×(Lon 关注 -Lon 震中 );

[0102] Module 36 is used to input the elliptical model into the transformed new coordinates (X', Y'), calculate the final value, compare it with 1, and determine the intensity circle of the key concern area based on the comparison result. Specifically, it includes:

[0103] The first judgment submodule is used when... Then it is determined that the area of ​​concern is located within the intensity zone, and its intensity value is N degrees;

[0104] The second judgment submodule is used when... Then it is determined that the area of ​​concern is located on the boundary of the intensity circle, and its intensity value is N degrees;

[0105] The third judgment submodule is used when... If the area of ​​interest is not within the intensity zone, the intensity value is reduced to N-1, and the first, second, and third judgment submodules are called repeatedly for judgment.

[0106] The embodiments of the present invention are device embodiments corresponding to the above method embodiments. The specific operation of each module can be understood with reference to the description of the method embodiments, and will not be repeated here.

[0107] Device Example 2

[0108] This invention provides an electronic device, such as... Figure 4 As shown, it includes: a memory 40, a processor 42, and a computer program stored in the memory 40 and executable on the processor 42, wherein the computer program, when executed by the processor 42, performs the steps as described in the method embodiment.

[0109] Device Example 3

[0110] This invention provides a computer-readable storage medium storing an information transmission implementation program, which, when executed by a processor 42, performs the steps described in the method embodiment.

[0111] The computer-readable storage media described in this embodiment include, but are not limited to, ROM, RAM, disk, or optical disk.

[0112] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some or all of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for determining the theoretical intensity of earthquakes in a region of interest, characterized in that, include: By using the earthquake influence field model of the area of ​​interest, the magnitudes of the major and minor axes of different intensity values ​​under different magnitudes in the area of ​​interest are derived. Determine the direction of the seismic influence field in the area of ​​interest; The elliptical model is determined based on the magnitude of the major and minor axes and the direction of the seismic influence field, and the elliptical model is rotated to obtain the new coordinates (X', Y'). Substitute the elliptical model into the new coordinate system (X', Y') after transformation, calculate the final value, compare it with 1, and determine the intensity zone of the key concern area based on the comparison result.

2. The method according to claim 1, characterized in that, Based on the earthquake influence field model of the area of ​​interest, the specific sizes of the major and minor axes at different intensity values ​​under different magnitudes in the area of ​​interest are derived as follows: By substituting magnitudes from 4.0 to 8.0 and intensity from 5 to 8 degrees into the seismic influence field model of the area of ​​interest, the lengths of the major and minor semi-axes of the elliptical attenuation model under different magnitudes and influence intensities can be derived.

3. The method according to claim 1, characterized in that, Determining the direction of the seismic influence field in the area of ​​interest specifically includes: The active faults in the area of ​​interest are decomposed into polyline objects, and the center value of each segment is taken as the latitude and longitude coordinates of the fault, and its attribute information is assigned. The data from the epicenter and the fault point are calculated and compared. The distance between the two points is calculated, and the fault with the minimum distance is taken as the fault that may induce an earthquake. The direction of the earthquake influence field is determined based on the direction of the rupture of the fault.

4. The method according to claim 1, characterized in that, Based on the magnitudes of the major and minor axes and the direction of the seismic influence field, an elliptical model is determined, and the elliptical model is rotated to obtain a new coordinate system (X', Y'). Specifically, this includes: Set the latitude and longitude coordinates of the earthquake epicenter as (Lon 震中 Lat 震中 The coordinates of the area of ​​interest are (Lon 关注 Lat 关注 The magnitude is M, the distance of the area of ​​concern from the epicenter is D, the intensity is N, and the semi-major axis of intensity N is R. N长 The minor half-axis is R N短 The fault strike is T, and it is assumed that the distance of 1 degree of longitude or 1 degree of latitude is 111 kilometers; The standard ellipse model is obtained according to formula 1-3: The semi-major axis radian ρ of intensity N 长 =R N长 / 111; Formula 1; minor semi-axis radian ρ 短 =R N短 / 111; Formula 2; Fault strike converted to arc system 断层 =π·T / 180; Formula 3; The standard elliptical model is translated and rotated in polar coordinates, with the epicenter as the center and the fault strike as the direction of rotation. The new coordinates after rotation are (X', Y'): ； 。 5. The method according to claim 4, characterized in that, Substitute the elliptical model into the transformed new coordinate system (X', Y'), calculate the final value, and compare it with 1. Based on the comparison result, determine the specific intensity zones of the key areas of concern, including: Step 1, when If so, the area of ​​concern is determined to be within the intensity zone, and its intensity value is N degrees; Step 2, when If so, the area of ​​concern is located on the boundary of the intensity circle, and its intensity value is N degrees. Step 3, when If the area of ​​concern is not within the intensity zone, the intensity value is reduced to N-1, and steps 1-3 are repeated for the same determination.

6. A device for determining the theoretical intensity of an earthquake in a region of interest, characterized in that, include: The derivation module is used to derive the major and minor axes of different intensity values ​​under different magnitudes in the area of ​​interest from the earthquake influence field model of the area of ​​interest; The determination module is used to determine the direction of the seismic influence field in the area of ​​interest; The ellipse model module is used to determine the ellipse model based on the size of the major and minor axes and the direction of the seismic influence field, and to rotate the ellipse model to obtain a new coordinate (X', Y'). The judgment module is used to input the elliptical model into the new coordinates (X', Y') after transformation, calculate the final value and compare it with 1, and determine the intensity circle of the key concern area based on the comparison result.

7. The apparatus according to claim 6, characterized in that, The derivation module is specifically used for: By substituting magnitudes from 4.0 to 8.0 and intensity from 5 to 8 degrees into the seismic influence field model of the area of ​​interest, the lengths of the major and minor semi-axes of the elliptical attenuation model under different magnitudes and influence intensities can be derived. The determining module is specifically used for: The active faults in the area of ​​interest are decomposed into polyline objects, and the center value of each segment is taken as the latitude and longitude coordinates of the fault, and its attribute information is assigned. The data from the epicenter and the fault point are calculated and compared. The distance between the two points is calculated, and the fault with the minimum distance is taken as the fault that may induce an earthquake. The direction of the earthquake influence field is determined based on the direction of the rupture of the fault.

8. The apparatus according to claim 6, characterized in that, The ellipse model module is specifically used for: Set the latitude and longitude coordinates of the earthquake epicenter as (Lon 震中 Lat 震中 The coordinates of the area of ​​interest are (Lon 关注 Lat 关注 The magnitude is M, the distance of the area of ​​concern from the epicenter is D, the intensity is N, and the semi-major axis of intensity N is R. N长 The minor half-axis is R N短 The fault strike is T, and it is assumed that the distance of 1 degree of longitude or 1 degree of latitude is 111 kilometers; The standard ellipse model is obtained according to formula 1-3: The semi-major axis radian ρ of intensity N 长 =R N长 / 111; Formula 1; minor semi-axis radian ρ 短 =R N短 / 111; Formula 2; Fault strike converted to arc system 断层 =π·T / 180; Formula 3; The standard elliptical model is translated and rotated in polar coordinates, with the epicenter as the center and the fault strike as the direction of rotation. The new coordinates after rotation are (X', Y'): ； ； The judgment module specifically includes: The first judgment submodule is used when... If so, the area of ​​concern is determined to be within the intensity zone, and its intensity value is N degrees; The second judgment submodule is used when... If so, the area of ​​concern is located on the boundary of the intensity circle, and its intensity value is N degrees. The third judgment submodule is used when... If the area of ​​interest is not within the intensity zone, the intensity value is reduced to N-1, and the first, second, and third judgment submodules are called repeatedly for judgment.

9. An electronic device, characterized in that, include: The memory, the processor, and the computer program stored in the memory and executable on the processor, wherein the computer program, when executed by the processor, implements the steps of the method for determining the theoretical seismic intensity of the area of ​​interest as described in any one of claims 1 to 5.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores an information transmission implementation program, which, when executed by a processor, implements the steps of the method for determining the theoretical seismic intensity of the area of ​​interest as described in any one of claims 1 to 5.

Images