A method for quickly calculating the impact fields of a large number of earthquake events

Abstract

The present invention relates to a method for rapidly calculating the impact fields of a large number of earthquake events. According to the theoretical method of seismic motion relationship attenuation, in the same earthquake attenuation zone, earthquakes with the same fault direction and the same magnitude cause the same seismic motion impact field, while the seismic motion impact field caused by earthquakes with different fault directions can be obtained by rotation. The beneficial effect of the present invention is that the present invention establishes a set of indexes and impact field templates of seismic motion relationship attenuation, thereby achieving the purpose of rapidly evaluating the impact of a set of earthquake events on a benchmark.

Description

Technical Field

[0001] The invention belongs to the field of earthquake catastrophe models, and in particular relates to a method for rapidly calculating the impact fields of a large number of earthquake events. Background Art

[0002] In earthquake catastrophe models, in order to assess the losses caused by medium- and long-term earthquakes, it is often necessary to simulate hundreds of millions of earthquake events. Assessing the impact of such a large number of earthquakes on targets (especially a large number of targets) requires a lot of storage space, computing time, and computing power. Summary of the invention

[0003] The method for constructing the index of the attenuation of the earthquake relationship and the influence field template provided by the present invention comprises the following steps:

[0004] Discretize the magnitude to get the earthquake sequence;

[0005] Establishing a two-dimensional array based on the earthquake sequence and earthquake attenuation partitions;

[0006] Acquire attenuation relationship parameters according to the magnitude and the earthquake attenuation partition;

[0007] The attenuation partition influence field reference data of the earthquake attenuation partition corresponding to the magnitude and attenuation relationship parameters are calculated by the following formula;

[0008] lgY=A+BM+Clg(R+De EM )

[0009] Among them, Y is the seismic parameter, R is the epicenter distance, M is the magnitude, and A, B, C, D, and E are all regression coefficients;

[0010] The reference data of the influence field of each attenuation partition is stored in a data file DataFile;

[0011] Storing the starting position and the ending position of each attenuation partition influence field reference data into the two-dimensional array;

[0012] The two-dimensional array is stored in the index file IndexFile.

[0013] The reference data of the influence field of the example may be organized as shown in the following table. It should be noted that there are multiple storage methods for the influence field. For example, for the sake of intuitiveness, the present invention is explained in text format.

[0014]

[0015] The present invention also provides a method for quickly calculating an influence field based on an index of earthquake motion relationship attenuation and an influence field template, comprising the following steps:

[0016] Obtain the spatial position of the target and the earthquake;

[0017] Load the index file and read the file content into a two-dimensional array;

[0018] According to the spatial position of the earthquake, the earthquake attenuation zone to which the earthquake belongs is determined; then, according to the magnitude of the earthquake and the earthquake attenuation zone, the row number and the column number are obtained;

[0019] Read the values ​​in the two-dimensional array according to the row number and the column number to obtain the starting position and the ending position of the impact field data;

[0020] Open the data file and read the data between the starting position and the ending position;

[0021] Taking the position coordinates of the earthquake epicenter as the origin, the position coordinates of the target are converted into relative position coordinates;

[0022] According to the strike angle of the earthquake fault, the target position is rotated in the opposite direction to obtain a new target position;

[0023] The seismic motion parameters are read from the template influence field according to the position, that is, the seismic motion of the target is affected by the earthquake.

[0024] The present invention also provides a computer device, comprising a memory, a processor and a computer program stored in the memory, wherein the processor executes the computer program to implement the steps of the above method.

[0025] The present invention also provides a computer-readable storage medium having a computer program / instruction stored thereon, which implements the steps of the aforementioned method when the computer program / instruction is executed by a processor.

[0026] The present invention also provides a computer program product, comprising a computer program / instruction, which implements the steps of the above method when executed by a processor.

[0027] The beneficial effect of the present invention is that the present invention establishes a set of indexes and influence field templates of earthquake motion relationship attenuation, thereby achieving the purpose of quickly evaluating the impact of earthquake event sets on the benchmark. BRIEF DESCRIPTION OF THE DRAWINGS

[0028] Figure 1 Schematic diagram of earthquake motion attenuation relationship zoning;

[0029] Figure 2 A sequence of earthquakes with magnitudes of 5 to 9 with a magnitude interval of 0.1;

[0030] Figure 3 Schematic diagram with earthquake sequences as rows and earthquake attenuation zones as columns;

[0031] Figure 4 Schematic diagram of a two-dimensional array;

[0032] Figure 5 Schematic diagram of the assignment of two-dimensional array and storage of data file DataFile;

[0033] Figure 6 Flow chart of the retrieval method of the embodiment; DETAILED DESCRIPTION

[0034] The following specific implementations can disclose many technical details, but these specific implementations cannot be used to limit the protection scope of the present invention.

[0035] First, the construction method of the index of the attenuation of the ground motion relationship and the impact field template includes the following steps:

[0036] Discretize the magnitude to get the earthquake sequence;

[0037] Establishing a two-dimensional array based on the earthquake sequence and earthquake attenuation partitions;

[0038] Acquire attenuation relationship parameters according to the magnitude and the earthquake attenuation partition;

[0039] The attenuation partition influence field reference data of the earthquake attenuation partition corresponding to the magnitude and attenuation relationship parameters are calculated by the following formula;

[0040] lgY=A+BM+Clg(R+De EM )

[0041] Among them, Y is the seismic parameter, R is the epicenter distance, M is the magnitude, and A, B, C, D, and E are all regression coefficients;

[0042] The reference data of the influence field of each attenuation partition is stored in a data file DataFile;

[0043] Storing the starting position and the ending position of each attenuation partition influence field reference data into the two-dimensional array;

[0044] The two-dimensional array is stored in the index file IndexFile.

[0045] In some embodiments, the earthquake magnitudes are discretized from small to large to obtain an earthquake sequence, and a two-dimensional array is established with the earthquake sequence as rows and the earthquake attenuation partitions as columns. In other embodiments, the discrete magnitude interval in the earthquake sequence is 0.1.

[0046] Among them, "magnitude" includes but is not limited to Richter magnitude; "earthquake attenuation zone" can be defined according to the theoretical method of earthquake motion relationship attenuation, for example, it can refer to "China Earthquake Motion Parameter Zoning Map (GB 18306-2015)", and can also be determined based on the corresponding data of other geographical regions. The term "establishing a two-dimensional array" can be understood as establishing a storage area of ​​a two-dimensional array on a storage medium. "Attenuation relationship parameters" are selected from regression coefficients including A, B, C, D, E, etc. When the attenuation relationship equation changes, increasing or decreasing the regression coefficient is within the concept of the present invention. The present invention does not specifically limit the method of obtaining the attenuation relationship parameters, including manual or computer reference to known data, actual measurement, software estimation, etc.

[0047] In the same earthquake attenuation zone, earthquakes with the same fault orientation and magnitude will cause the same seismic motion impact field, while the seismic motion impact field caused by earthquakes with different fault orientations can be obtained by rotation.

[0048] In some embodiments of the method for quickly calculating the influence field based on the index of the attenuation of the ground motion relationship and the influence field template, the following steps are included:

[0049] Obtain the spatial position of the target and the earthquake;

[0050] Load the index file and read the file content into a two-dimensional array;

[0051] According to the spatial position of the earthquake, the earthquake attenuation zone to which the earthquake belongs is determined; then, according to the magnitude of the earthquake and the earthquake attenuation zone, the row number and the column number are obtained;

[0052] Read the values ​​in the two-dimensional array according to the row number and the column number to obtain the starting position and the ending position of the impact field data;

[0053] Open the data file and read the data between the starting position and the ending position;

[0054] Taking the position coordinates of the earthquake epicenter as the origin, the position coordinates of the target are converted into relative position coordinates;

[0055] According to the strike angle of the earthquake fault, the target position is rotated in the opposite direction to obtain a new target position;

[0056] The seismic motion parameters are read from the template influence field according to the position, that is, the seismic motion of the target is affected by the earthquake.

[0057] The following further describes implementation methods of some embodiments of the present invention.

[0058] In some embodiments, reference Figure 1As shown, the distribution of each earthquake attenuation zone is exemplified. According to the theoretical method of seismic motion relationship attenuation, in the same earthquake attenuation zone, earthquakes with the same fault direction and the same magnitude will cause the same seismic motion influence field, while the seismic motion influence field caused by earthquakes with different fault directions can be obtained by rotation.

[0059] In some specific embodiments, the earthquake magnitude is discretized from small to large to obtain an earthquake sequence. In a specific embodiment, the precision, i.e., discretization, needs to be set during the implementation process, such as Figure 2 As shown in the figure, taking the magnitude interval of 0.1, 5 to 9 earthquake as an example, the following steps are continued:

[0060] 1. Create a two-dimensional array with earthquake sequences as rows and earthquake attenuation partitions as columns. Please refer to Figure 3 and Figure 4 .

[0061] 2. According to the magnitude and earthquake zone, the parameters of the attenuation relationship are obtained. In this specific embodiment, see Table 1, and the parameters of the attenuation relationship can be obtained from the table. For a magnitude 7 earthquake, the distance to the medium-strong earthquake zone, the parameters of the long axis are shown in Table 2, and the parameters of the short axis are shown in Table 3:

[0062] Table 1 Earthquake attenuation coefficients for each sub-area

[0063]

[0064] Table 2: Parameters of the major axis

[0065] M A B C D E 7 3.706 0.298 -2.079 2.802 0.295

[0066] Table 3: Parameters of the short axis

[0067] M A B C D E 7 2.690 0.321 -1.723 1.295 0.331

[0068] 3. Based on the magnitude and parameters obtained in step 4, the corresponding impact field can be calculated using the following formula.

[0069] lgY=A+BM+Clg(R+De EM )

[0070] Among them, Y is the seismic parameter, R is the epicenter distance, M is the magnitude, and A, B, C, D, and E are all regression coefficients;

[0071] 4. Store the calculated reference data of the influence field of each attenuation partition into the data file DataFile, refer to Figure 5 shown.

[0072] 5. Store the starting and ending positions of the data in the above two-dimensional array, refer to Figure 5 shown.

[0073] 6. Repeat steps 2 to 5 to calculate each value of the above two-dimensional array.

[0074] 7. Store this two-dimensional array in a file (index file IndexFile).

[0075] Search Method

[0076] refer to Figure 6 As shown, the position of the target and the spatial position of the earthquake are obtained;

[0077] Load the index file IndexFile and read the file content into the two-dimensional array IDX[,];

[0078] According to the spatial position of the earthquake, the earthquake attenuation zone to which the earthquake belongs is determined; then according to the magnitude of the earthquake and the earthquake attenuation zone, the row number r and the column number c are obtained;

[0079] Read the value IDX[r,c] in the two-dimensional array IDX[,] according to the row number r and the column number c to obtain the starting position As and the ending position Ae of the impact field data;

[0080] Open the data file DataFile and read the data between the starting position As and the ending position Ae;

[0081] Taking the location coordinates (x0, y0) of the earthquake epicenter as the origin, convert the location coordinates (x, y) of the target into relative location coordinates (x′, y′):

[0082]

[0083] According to the strike angle θ of the earthquake fault, the target position (x′, y′) is rotated in the opposite direction by θ to obtain the new target position (x″, y″)

[0084]

[0085] The seismic motion parameters are read from the template influence field according to the position (x″, y″), that is, the seismic motion of the target is affected by the earthquake.

[0086] Storage space comparison test example

[0087] In the method of this specific embodiment, in the same earthquake attenuation zone, earthquakes of the same magnitude (this method regards earthquakes with magnitudes less than a threshold value as having the same magnitude, and the threshold value is set to 0.1 m / s in the above example) 2 ) earthquake uses one impact field data. Therefore, in this specific embodiment, only 164 impact field data need to be stored (the number of impact fields that need to be stored and calculated will be different depending on the threshold, such as the threshold is set to 0.01m / s 2, 1640 impact fields need to be stored). In the above implementation, there are 41 magnitudes to be calculated for earthquakes of magnitude 5 to 9 with an interval of 0.1, and each magnitude needs to be calculated for 4 attenuation partitions, so 4*41=164 impact field data need to be calculated. By analogy, when the interval is 0.01, 4*410=1640 impact field data need to be calculated. Traditional methods generally save one impact field data for each earthquake. The insurance industry's catastrophe model needs to use simulated events of 100,000 or even 1 million years for medium- and long-term earthquake loss assessment (taking 1 million years as an example, China alone has more than 60 million simulated events), so it can save about 300,000 times the storage space (the threshold is set to 0.01m / s 2 When the number of earthquake events increases, the storage space is saved by 30,000 times, and so on). The query speed of this method does not increase with the increase of the number of earthquake events, and the time complexity is O(1).

[0088] The embodiments and functional operations of the subject matter described in this specification may be implemented in: digital electronic circuitry, tangibly implemented computer software or firmware, computer hardware, including the structures disclosed in this specification and their structural equivalents, or a combination of more than one of the above. The embodiments of the subject matter described in this specification may be implemented as one or more computer programs, i.e., one or more modules of computer program instructions encoded on one or more tangible non-transitory program carriers, for execution by a data processing device or to control the operation of a data processing device.

[0089] A computer program (which may also be referred to or described as a program, software, software application, module, software module, script or code) may be written in any form of programming language, including compiled or interpreted languages ​​or declarative or procedural languages, and the computer program may be deployed in any form, including as a stand-alone program or as a module, component, subroutine or other unit suitable for use in a computing environment. A computer program may, but need not, correspond to a file in a file system. A program may be stored in a portion of a file that stores other programs or data, for example, one or more scripts stored in the following: in a markup language document; in a single file dedicated to a related program; or in multiple collaborative files, for example, files storing one or more modules, subroutines or code portions. A computer program may be deployed to execute on one or more computers, the computers being located in one place or distributed to multiple locations and interconnected by a communication network.

[0090] The processes and logic flows described in this specification may be performed by one or more programmable computers that execute one or more computer programs to perform functions by operating on input data and generating output. The processes and logic flows may also be performed by special purpose logic circuits, such as FPGAs (field programmable gate arrays) or ASICs (application specific integrated circuits), and the apparatus may also be implemented as special purpose logic circuits.

[0091] The computer suitable for implementing the computer program includes and exemplarily can be based on a general-purpose microprocessor or a special-purpose microprocessor or both of the above-mentioned processors, or any other kind of central processing unit. Generally, the central processing unit will receive instructions and data from a read-only memory or a random access memory or both. The main element of the computer is a central processing unit for running or executing instructions and one or more memory devices for storing instructions and data. Generally, the computer will also include or be operably coupled to receive data from one or more large-capacity storage devices for storing data or transfer data to the large-capacity storage device, or receive and transfer both, and the large-capacity storage device is, for example, a magnetic disk, a magneto-optical disk or an optical disk. However, the computer does not necessarily have such a device. In addition, the computer can be embedded in another device, for example, a mobile phone, a personal digital assistant (PDA), a mobile audio or video player, a game console, a global positioning system (GPS) receiver or a removable storage device, for example, a universal serial bus (USB) flash drive, etc.

[0092] Computer readable media suitable for storing computer program instructions and data include all forms of nonvolatile memory, media and memory devices, including, as examples: semiconductor memory devices, such as EPROM, EEPROM and flash memory devices; magnetic disks, such as internal hard disks or removable disks; magneto-optical disks; CD-ROM and DVD-ROM disks. The processor and memory may be supplemented by or incorporated into a special purpose logic circuit.

[0093] Although this specification contains many specific implementation details, these should not be interpreted as limitations on the scope of any invention or the scope of what can be claimed, but rather as an illustration of features that can concretize a particular embodiment of a particular invention. The specific features described in this specification in the context of independent embodiments can also be implemented in combination with a single embodiment. Conversely, the various features described in the context of a single embodiment can also be implemented independently in multiple embodiments, or in any suitable sub-combination. In addition, although the features can be described above as a combination effect and even initially required, one or more features from the required combination can be removed from the combination in some cases, and the required combination can be turned to a sub-combination or a variation of the sub-combination.

[0094] Similarly, although operations are described in the accompanying drawings in a particular order, it should not be understood that: in order to achieve the desired results, such operations are required to be performed in the particular order shown or in a sequential order, or that all illustrated operations are performed. In certain cases, multitasking and parallel processing can be advantageous. In addition, the separation of various system modules and components in the above-mentioned embodiments should not be understood as requiring such separation in all embodiments, and it should be understood that program components and systems can generally be integrated into a single software product or packaged into multiple software products.

[0095] Specific implementations of the subject matter have been described. Other implementations are within the scope of the following claims. For example, the activities recited in the claims can be performed in a different order and still achieve the desired results. As an example, the processes described in the accompanying drawings do not necessarily require the specific order or sequential order shown in order to achieve the desired results. In certain implementations, multitasking and parallel processing may be advantageous.

Claims

1. A method for constructing an index of earthquake motion relationship attenuation and an influence field template, characterized in that: The steps include: Discretizing the magnitude to obtain an earthquake sequence, the specific steps are: discretizing the earthquake magnitude from small to large to obtain an earthquake sequence, wherein the discrete magnitude interval in the earthquake sequence is 0.1; A two-dimensional array is established with the earthquake sequence and the earthquake attenuation partitions, and the specific steps are: a two-dimensional array is established with the earthquake sequence as rows and the earthquake attenuation partitions as columns; Acquire attenuation relationship parameters according to the magnitude and the earthquake attenuation partition; The attenuation partition influence field reference data of the earthquake attenuation partition corresponding to the magnitude and attenuation relationship parameters are calculated by the following formula; lgY=A+BM+Clg(R+De EM , Among them, Y is the seismic parameter, R is the epicenter distance, M is the magnitude, and A, B, C, D, and E are all regression coefficients; The reference data of the influence field of each attenuation partition is stored in a data file DataFile; Storing the starting position and the ending position of each attenuation partition influence field reference data into the two-dimensional array; The two-dimensional array is stored in the index file IndexFile; In the same earthquake attenuation zone, earthquakes with the same fault orientation and magnitude will cause the same seismic motion impact field, while the seismic motion impact field caused by earthquakes with different fault orientations can be obtained by rotation.

2. A method for quickly calculating an influence field based on the index of the attenuation of the ground motion relationship and the influence field template according to claim 1, characterized in that: The following steps are included: Obtain the spatial position of the target and the earthquake; Load the index file and read the file content into a two-dimensional array; According to the spatial position of the earthquake, the earthquake attenuation zone to which the earthquake belongs is determined; then, according to the magnitude of the earthquake and the earthquake attenuation zone, the row number and the column number are obtained; Read the values ​​in the two-dimensional array according to the row number and the column number to obtain the starting position and the ending position of the impact field data; Open the data file and read the data between the starting position and the ending position; Taking the position coordinates of the earthquake epicenter as the origin, the position coordinates of the target are converted into relative position coordinates; According to the strike angle of the earthquake fault, the target position is rotated in the opposite direction to obtain a new target position; The seismic motion parameters are read from the template influence field according to the position, that is, the seismic motion of the target is affected by the earthquake.

3. A method for quickly calculating an influence field based on the index of the attenuation of the ground motion relationship and the influence field template of claim 1, characterized in that: The following steps are included: Obtain the spatial position of the target and the earthquake; Load the index file IndexFile and read the file content into the two-dimensional array IDX[,]; According to the spatial position of the earthquake, the earthquake attenuation zone to which the earthquake belongs is determined; then according to the magnitude of the earthquake and the earthquake attenuation zone, the row number r and the column number c are obtained; Read the value IDX[r,c] in the two-dimensional array IDX[,] according to r and c to obtain the starting position As and the ending position Ae of the impact field data; Open the data file DataFile and read the data between the starting position As and the ending position Ae; Taking the location coordinates (x0, y0) of the earthquake epicenter as the origin, convert the location coordinates (x, y) of the target into relative location coordinates (x', y'): According to the strike angle θ of the earthquake fault, the target position (x', y') is rotated in the opposite direction by θ to obtain the new target position (x", y"). Read the seismic motion parameters from the template influence field according to the position (x", y"), that is, the seismic motion of the target is affected by the earthquake.

4. A computer device comprising a memory, a processor and a computer program stored in the memory, characterized in that: The processor executes the computer program to implement the steps of the method according to any one of claims 1 to 3.

5. A computer-readable storage medium having a computer program / instruction stored thereon, characterized in that: When the computer program / instructions are executed by a processor, the steps of the method according to any one of claims 1 to 3 are implemented.

6. A computer program product comprising a computer program / instructions, characterized in that When the computer program / instructions are executed by a processor, the steps of the method according to any one of claims 1 to 3 are implemented.

Images