A method and device for rapid earthquake event retrieval

Abstract

The present invention relates to a method and device for rapid retrieval of seismic events. The beneficial effect of the present invention is that, in view of the impact characteristics of earthquakes, a set of rapid retrieval methods are designed, which can quickly retrieve seismic events that may cause losses. This method first needs to create a spatial index of seismic events. Although traditional indexes can quickly retrieve the locations of data, since the relevant data are often stored in different places, it is impossible to read them in one piece at a time. Frequent reading operations reduce the speed of data retrieval. The present invention reorganizes the storage order and location of data according to the spatial positions of earthquakes and ensures that the time complexity of retrieval is O(1).

Description

Technical Field

[0001] The present invention belongs to the field of earthquake catastrophe models, and particularly relates to a method and device for quickly retrieving earthquake events. Background Art

[0002] In earthquake catastrophe models, in order to evaluate the losses caused by medium- and long-term earthquakes, it is often necessary to simulate tens of millions of earthquake events. Taking China and its surrounding areas as an example, these earthquake events are distributed throughout China, and their influence ranges are different according to the different epicenter positions and magnitudes of the earthquakes. When evaluating the medium- and long-term impact of an earthquake on a certain target, it is often not necessary to calculate each earthquake. According to the position of the target and the earthquake epicenter, as well as the earthquake magnitude, quickly screening out the earthquake events that may affect the target is a key technology that people expect but has not been reported, which affects the calculation speed of earthquake catastrophe models. Summary of the Invention

[0003] One aspect of the present invention provides a method for constructing a spatial index of earthquake events, including the following steps:

[0004] Dividing the earthquake distribution range into several grids;

[0005] According to the epicenter position of each earthquake event within the earthquake distribution range, attributing the earthquake data corresponding to the epicenter position to each grid;

[0006] Storing the earthquake data in the same grid continuously;

[0007] Establishing a Hash index for the grids, where the Hash index includes a unique identifier key and a value value corresponding to each grid, and the value value is the storage position and length of the earthquake data corresponding to the grid.

[0008] Another aspect of the present invention further provides a method for quickly retrieving earthquake events based on the aforementioned spatial index of earthquake events, including the following steps:

[0009] Selecting several nearby grids according to the target position;

[0010] Querying the storage position and length of the earthquake data in each grid in the Hash index of each grid according to the grid key;

[0011] Reading all the earthquake events corresponding to each grid in full;

[0012] Analyzing the impact of each earthquake event on the target.

[0013] Another aspect of the present invention further provides a computer program product, including a computer program, which when executed by a processor implements the steps of the aforementioned method.

[0014] Another aspect of the present invention also provides a computer device, the device comprising at least one processor; and

[0015] The memory stores instructions, and when the instructions are executed by at least one processor, the steps according to the above method are implemented.

[0016] The beneficial effect of the present invention is that a set of rapid retrieval methods is designed according to the impact characteristics of earthquakes, and the method can quickly retrieve earthquake events that may cause losses. The method first needs to create a spatial index of earthquake events. Although the traditional index can quickly retrieve the location of data, it is impossible to read the entire block at one time because related data is often stored in different places. Frequent reading operations reduce the speed of data retrieval. The present invention reorganizes the storage order and position of data according to the spatial position of the earthquake, and ensures that the time complexity of retrieval is O(1). BRIEF DESCRIPTION OF THE DRAWINGS

[0017] Figure 1 Assign earthquake data to each grid diagram according to the location of the earthquake epicenter;

[0018] Figure 2 Schematic diagram of continuous storage of earthquake and its impact field data in the same grid;

[0019] Figure 3 Schematic diagram of establishing a Hash index for a grid; DETAILED DESCRIPTION

[0020] Some embodiments relate to a method for constructing a spatial index of seismic events, comprising the steps of:

[0021] Divide the earthquake distribution range into several grids;

[0022] According to the earthquake epicenter position of each earthquake event within the earthquake distribution range, the earthquake data corresponding to the earthquake epicenter position are assigned to each grid;

[0023] Continuously store the seismic data of the same grid;

[0024] A Hash index of the grid is established, wherein the Hash index includes a corresponding unique identifier key and a value value of each grid, wherein the value value is the storage location and length of the corresponding seismic data in the grid.

[0025] The term "seismic distribution range" refers to, for example, the distribution range of a set of seismic random events. The set of seismic random events refers to a set of events simulated based on potential source regions and their seismic activity parameters. The GIS technology can be used to mark the locations of all earthquakes using coordinates (or longitude and latitude). The essence of the vector elements formed by these seismic events is the seismic distribution range referred to herein.

[0026] The term "division" refers to, for example, using the visualization interface provided by a computer to separate the seismic events in the target geographical area into several retrievable or locatable grids, or, for another example, establishing a mapping between the location information and grid information of the target geographical area in a Geographic Information System (GIS).

[0027] The term "grid" refers to the smallest area of the target geographical area characterized by location information.

[0028] The term "epicenter location" refers to a set of data indicating the location of the epicenter in a seismic event. For example, each set of data includes longitude and latitude.

[0029] The term "assigned to each grid" refers to, for example, assigning a seismic event to the grid that contains the epicenter location of the earthquake.

[0030] The term "seismic data" includes but is not limited to magnitude, epicenter location of the earthquake, fault strike, focal depth, and the resulting influence field, etc.

[0031] The term "influence field data" refers to, for example, the ground motion parameters at each location around the earthquake after the earthquake occurs. The influence field data has various storage and presentation forms. Commonly, it is in the form of a raster file, with the ground motion parameters within each equal-area spatial range.

[0032] The term "continuous storage" means that the storage locations of data on a readable storage medium are continuous based on a certain storage rule, achieving the storage of spatially adjacent data together.

[0033] The term "Hash index" is the hash index. It is implemented based on a hash table, and only queries that precisely match all columns of the index are effective. For each row of data, the storage engine calculates a hash code for all index columns. The hash code is a relatively small value, and the hash codes calculated for rows with different key values are also different. The hash index stores all the hash codes in the index and simultaneously saves pointers to each data row in the hash table.

[0034] The term "storage location and length" refers to the address for storing seismic data and its influence field data, and the length occupied by the seismic data and its influence field data on the storage medium;

[0035] In some other embodiments, the method further includes the steps of:

[0036] obtaining the grid size Rm; and then,

[0037] dividing the seismic distribution range into different grids according to the grid size Rm, where

[0038] the grid size Rm is calculated in terms of distance or longitude and latitude. Preferably, the distance is 100 km and the longitude and latitude are 1°.

[0039] In some other embodiments, the seismic data in each of the grids is called simultaneously.

[0040] In some other embodiments, the seismic data does not include the influence field data, and the data is obtained by reading seismic information and performing real-time calculation or query. In some specific embodiments, after reading the seismic data, the real-time calculation or query method is adopted. The overall process is consistent with the foregoing optional embodiments, and only the logic of real-time calculation or query needs to be added.

[0041] In some other embodiments, the unique identifier key is obtained through the lower left coordinates (long, lat) by a formula: .

[0042] The present invention will be further explained below with reference to the accompanying drawings.

[0043] 1. Select the grid size Rm according to the project requirements. The setting of this parameter can refer to the influence range calculated according to the attenuation relationship of the earthquake (please refer to GB18306-2015, and the detailed description is omitted). It is usually expressed in terms of distance or longitude and latitude, such as 100 km or 1° (in the following examples, 1° is used as the grid size for illustration, and the actual parameter determination should consider the attributes of the object to be analyzed, zoning, site soil characteristics, etc.).

[0044] 2. Divide the seismic distribution range into different grids according to the grid size Rm.

[0045] 3. According to the epicenter position of the earthquake, attribute the seismic data to each grid, see Figure 1 .

[0046] 4. Continuously store the earthquakes and their influence field data in the same grid, see Figure 2 .

[0047] 5. Establish a Hash index for the grid: The key is the unique identifier of the above grid, and this identifier can be obtained through the lower left coordinates (long, lat) by a formula: ; the value is the address of the seismic data stored in the grid and the length occupied by the seismic data. As shown in the following table, specifically see Figure 3 .

[0048]

[0049] Retrieval

[0050] 1. Select the nearby grids according to the target location.

[0051] 2. Query the storage locations and lengths of the seismic data and impact field data in the grid in the Hash index of the grid according to the grid key.

[0052] 3. Read all the seismic events in these grids.

[0053] 4. Analyze the impact on the target for each event based on the results of step 3.

[0054] It should be noted that the present invention does not limit the specific analysis method, and various analysis methods for the impact of various targets by seismic events are within the selection of the present invention. For example, an analysis method including the following steps:

[0055] 1. Calculate the seismic motion parameters at the location of the target according to the fault trend of the seismic event, the attenuation zone where the epicenter of the earthquake is located, and the distance between the target and the epicenter of the earthquake, in accordance with the seismic motion parameter attenuation relationship model (GB18306 - 2015);

[0056] 2. Simulate the damage suffered by the target when the above - mentioned seismic motion parameters occur according to the type of the target building structure.

[0057] The embodiments and functional operations of the subject matter described in this specification can be implemented in: digital electronic circuits, tangibly implemented computer software or firmware, computer hardware, including the structures disclosed in this specification and structural equivalents thereof, or a combination of one or more of the above. The embodiments of the subject matter described in this specification can 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 being executed by a data processing device or for controlling the operation of a data processing device.

[0058] A computer program (which may also be referred to or described as a program, software, software application, module, software module, script, or code) can be written in any form of programming language, including compiled languages or interpreted languages or declarative or procedural languages, and the computer program can 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. The program can be stored in a part of a file that holds other programs or data, e.g., stored in one or more scripts in: a markup language document; a single file dedicated to the relevant program; or multiple co-operating files, e.g., files that hold one or more modules, subroutines, or portions of code. The computer program can be deployed to execute on one computer or on multiple computers, which may be located at one site, or distributed across multiple sites and interconnected by a communication network.

[0059] Although this specification contains many specific implementation details, these should not be construed as limitations on the scope of any invention or on the scope of what can be claimed, but rather as illustrations of features that may be specific to particular embodiments of a particular invention. The specific features described in the context of separate embodiments in this specification can also be implemented in combination in a single embodiment. Conversely, the various features described in the context of a single embodiment can also be implemented independently

[0060] in multiple embodiments, or in any suitable sub-combination. Moreover, although the features may be described above as acting in combination and even initially claimed as such, one or more features from a claimed combination can in some cases be removed from that combination, and the claimed combination can be directed to a sub-combination or variation of a sub-combination.

[0061] Similarly, although operations are depicted in the figures in a particular order, it should not be understood that such operations are required to be performed in the particular order shown or in sequential order in order to achieve the desired results, or that all illustrated operations are to be performed. In certain circumstances, multitasking and parallel processing may be advantageous. Additionally, the separation of various system modules and components in the above embodiments should not be understood as required in all embodiments, and it should be understood

[0062] that program components and systems can generally be integrated in a single software product or packaged into multiple software products.

[0063] Specific embodiments of the subject matter have been described. Other embodiments are within the scope of the following claims. For example, the acts recited in the claims can be performed in a different order and still achieve the desired result. As one example, in order to achieve the desired result, the processes described in the figures do not necessarily require the particular order or sequential order shown. In certain implementations, multitasking and parallel processing may be advantageous.

Claims

1. A method for retrieving fast seismic events based on a spatial index of seismic events, characterized in that, Including the following steps: Select several grids near the target according to the target location; According to the unique identifier key of each grid, query the storage location and length of seismic data and the storage location and length of influence field data in the Hash index of each grid; Read all seismic events corresponding to each grid; Analyze the influence of each seismic event on the target; Among them, the method for constructing the spatial index of the seismic event includes the following steps: Obtain the grid size Rm; then, divide the seismic distribution range into several different grids according to the grid size Rm, where the grid size Rm is calculated by distance or longitude and latitude; the seismic distribution range is the distribution range of the seismic random event set; the seismic random event set refers to the event set simulated based on the potential seismic source area and its seismic activity parameters, and the positions of all earthquakes are calibrated using coordinates or longitude and latitude; According to the epicenter positions of each seismic event within the seismic distribution range, assign the seismic data corresponding to the epicenter positions to each grid; the epicenter positions include longitude, latitude, and azimuth; Store the seismic data and its influence field data of the same grid continuously; the seismic data includes any one of magnitude, epicenter position of the earthquake, fault strike, and focal depth; the influence field data refers to the ground motion parameters at each position around after the earthquake occurs; Build a Hash index for the grid. The Hash index includes a corresponding unique identifier key and value for each grid. The value is the storage location and length of the corresponding seismic data in the grid. The unique identifier key is obtained by calculating according to the formula using the longitude and latitude of the lower left coordinate. The formula is: where long represents longitude and lat represents latitude.

2. The method according to claim 1, characterized in that, The distance is 100 km and the longitude and latitude are 1°.

3. The method according to claim 1, wherein The seismic data in each of the grids is called simultaneously.

4. The method according to claim 1, wherein The seismic data is data obtained by real-time calculation or query of seismic information.

5. A computer program product, comprising a computer program, characterized in that, When the computer program is executed by a processor, it implements the steps of the method described in claims 1-4.

6. A computer device, characterized in that, The device includes at least one processor; and a memory that stores instructions, and when the instructions are executed by at least one processor, the method described in any one of claims 1-4 is implemented.

Images