Method for calculating seismic intensity scales of surface using table-based calculation

Abstract

A method for calculating seismic intensity according to a disclosed embodiment comprises: receiving information on the magnitude and location of the epicenter of an earthquake from an external source; determining the epicentral distance based on the received location information of the epicenter based on a first table containing an earthquake ID (IDentification) pre-classified according to the location information and a pre-stored data file; determining the maximum bedrock acceleration based on the determined epicentral distance based on a second table containing the maximum bedrock acceleration pre-calculated according to the magnitude and epicentral distance of the earthquake; calculating a site amplification ratio by applying the determined maximum bedrock acceleration based on a third table containing a plurality of coefficient values; and calculating the seismic intensity at the surface of the earth according to the earthquake based on the calculated site amplification ratio.

Description

Technology Field

[0001] The present invention relates to a method for calculating progress based on table operations and a device for calculating progress applying the same. Background Technology

[0002] The existing method for calculating seismic intensity determines the magnitude and epicenter location (latitude and longitude) of the earthquake at the time of occurrence, and calculates the distance from the epicenter to each grid cell in the domestic inland. Furthermore, the existing method proceeded with seismic intensity calculation based on a bedrock motion prediction equation that uses the calculated distance to the grid, as well as the earthquake's magnitude and location, as parameters. Additionally, the existing method calculated the surface seismic intensity by determining the seismic intensity in the bedrock from the aforementioned bedrock motion prediction equation and determining the site amplification at each location as a ground parameter.

[0003] Conventionally, a program called Shake Map was used to calculate the intensity of an earthquake at each point on the Korean Peninsula based on the location and depth of the epicenter and the magnitude of the earthquake. However, since the existing method was utilized for early earthquake warning, fast computation speeds were required, and low-density grid information was used for this purpose. Consequently, the existing method suffered from a problem where the accuracy of intensity results for each administrative district was low. The problem to be solved

[0004] According to one disclosed embodiment, the invention relates to a seismic intensity calculation method and a seismic intensity calculation device applying the same, wherein seismic intensity can be displayed in a short time even on a high-density grid by calculating the seismic intensity of the surface based on a table, thereby shortening the announcement time for an earthquake early warning. means of solving the problem

[0005] A method for calculating seismic intensity according to a disclosed embodiment comprises: receiving information on the magnitude and location of the epicenter of an earthquake from an external source; determining the epicentral distance based on the received location information of the epicenter based on a first table containing an earthquake ID (IDentification) pre-classified according to the location information and a pre-stored data file; determining the maximum bedrock acceleration based on the determined epicentral distance based on a second table containing the maximum bedrock acceleration pre-calculated according to the magnitude and epicentral distance of the earthquake; calculating a site amplification ratio by applying the determined maximum bedrock acceleration based on a third table containing a plurality of coefficient values; and calculating the seismic intensity at the surface of the earth according to the earthquake based on the calculated site amplification ratio.

[0006] Determining the epicentral distance may include determining the ID of the earthquake that occurred based on the first table; loading a data file corresponding to the determined earthquake ID from a database storing data files classified by the earthquake ID; and determining the epicentral distance according to the location information of the received epicenter based on the loaded data file.

[0007] The above data file may include the epicentral distance from the epicenter included in the earthquake ID to a plurality of grids classified at intervals of 0.01° according to latitude and longitude.

[0008] Calculating the site amplification ratio may include determining at least three coefficient values ​​predetermined according to the IDs of the plurality of grids in the third table, and calculating the site amplification ratio based on the determined coefficient values ​​and the bedrock maximum acceleration.

[0009] Calculating the seismic intensity at the above-mentioned surface may include calculating the maximum acceleration at the surface based on the calculated site amplification ratio and the maximum acceleration of the bedrock, and calculating the seismic intensity at the surface due to the earthquake based on the calculated maximum acceleration at the surface.

[0010] It may further include calculating the maximum seismic value in the administrative district based on the maximum acceleration of the surface; and generating a seismic distribution map based on the maximum acceleration of the surface.

[0011] The first table, the second table, and the third table may be characterized by being converted into N-dimensional vector values ​​and stored in the database.

[0012] A seismic intensity calculation device according to another disclosed embodiment comprises: a communication unit that receives information on the magnitude of an earthquake and the location of the epicenter from an external source; a processor that calculates the seismic intensity of the earthquake on the surface based on the received information on the magnitude of the earthquake and the location of the epicenter; and a database that stores a first table, a second table, a third table, and a data file necessary for calculating the seismic intensity. The processor determines the epicentral distance according to the received location information of the epicenter based on the first table and the data file, which include an earthquake ID pre-classified according to location information; determines the maximum bedrock acceleration according to the determined epicentral distance based on the second table, which includes a maximum bedrock acceleration pre-calculated according to the magnitude of the earthquake and the epicentral distance; calculates a site amplification ratio by applying the determined maximum bedrock acceleration based on the third table, which includes a plurality of coefficient values; and calculates the seismic intensity on the surface of the earthquake that occurred on the surface based on the calculated site amplification ratio.

[0013] The processor can determine the ID of the earthquake that occurred based on the first table; load a data file corresponding to the determined earthquake ID from a database storing data files classified by the earthquake ID; and determine the epicentral distance based on the location information of the received epicenter based on the loaded data file.

[0014] The above data file may include the epicentral distance from the epicenter included in the earthquake ID to a plurality of grids classified at intervals of 0.01° according to latitude and longitude.

[0015] The processor can determine at least three coefficient values ​​predetermined according to the IDs of the plurality of grids in the third table, and calculate the site amplification ratio based on the determined coefficient values ​​and the bedrock maximum acceleration.

[0016] The processor can calculate the maximum surface acceleration based on the calculated site amplification ratio and the maximum bedrock acceleration, and calculate the seismic intensity based on the calculated maximum surface acceleration.

[0017] After the seismic intensity resulting from the earthquake is calculated, the processor can calculate the maximum seismic intensity value in an administrative district based on the maximum acceleration of the surface and generate a seismic intensity distribution map based on the maximum seismic intensity value based on the maximum acceleration of the surface.

[0018] The first table, the second table, and the third table may be characterized by being converted into N-dimensional vector values ​​and stored in the database.

[0019] A computer-readable recording medium storing a program for executing a control method of an electronic device according to another disclosed embodiment stores a first table, a second table, a third table, and a data file necessary for calculating the seismic intensity of an earthquake in advance; receives information on the magnitude of an earthquake and the location of the epicenter from an external source; determines the epicentral distance according to the received location information of the epicenter based on the first table and the previously stored data file, which include an earthquake ID classified in advance according to the location information; determines the maximum bedrock acceleration according to the determined epicentral distance based on the second table, which includes the maximum bedrock acceleration calculated in advance according to the magnitude of the earthquake and the epicentral distance; calculates a site amplification ratio by applying the determined maximum bedrock acceleration based on a third table including a plurality of coefficient values; and calculates the seismic intensity of an earthquake that occurs in advance based on the calculated site amplification ratio. Effects of the invention

[0020] The seismic intensity calculation method and seismic intensity calculation device applying the same according to one aspect disclosed can calculate the seismic intensity of the surface of the earth based on a table, thereby enabling seismic intensity display within a short time even on a high-density grid, and thus shortening the announcement time for earthquake early warnings. Brief explanation of the drawing

[0021] Figure 1 is a diagram for schematically illustrating the main players of a system to which the disclosed progress calculation method is applied. FIG. 2 is a control block diagram of a seismic calculation device (100) that performs the disclosed seismic calculation method. Figure 3 is a flowchart illustrating a disclosed method for calculating progress. Figure 4 is a table to explain the first table. Figure 5 is a diagram illustrating a data file. Figure 6 is a table to explain the second table. Figure 7 is a table to explain the third table. Figure 8 is a diagram illustrating the effect of the disclosed seismic calculation method. Specific details for implementing the invention

[0022] Throughout the specification, the same reference numerals refer to the same components. This specification does not describe all elements of the embodiments, and general content in the art to which the invention pertains or content that overlaps between embodiments is omitted.

[0023] Throughout the specification, when a part is described as being "connected" to another part, this includes not only cases where they are directly connected but also cases where they are indirectly connected, and indirect connections include connections made via a wireless communication network.

[0024] Furthermore, when it is stated that a part "includes" a certain component, this means that, unless specifically stated otherwise, it does not exclude other components but may include additional components.

[0025] Singular expressions include plural expressions unless there is an obvious exception in the context.

[0026] In addition, terms such as "~part," "~unit," "~block," "~part," and "~module" may refer to a unit that processes at least one function or operation. For example, the above terms may refer to at least one piece of hardware such as an FPGA (field-programmable gate array) or ASIC (application specific integrated circuit), at least one piece of software stored in memory, or at least one process processed by a processor.

[0027] The symbols attached to each step are used to identify each step and do not indicate the order of the steps relative to one another; the steps may be performed differently from the specified order unless a specific order is clearly indicated in the context.

[0028] Hereinafter, an embodiment regarding a method for calculating advance according to one aspect will be described in detail with reference to the attached drawings.

[0029] Figure 1 is a diagram for schematically illustrating the main players of a system to which the disclosed progress calculation method is applied.

[0030] Referring to FIG. 1, a system (1) according to a disclosed embodiment includes a meteorological agency (10) that detects the occurrence of an earthquake and a seismic intensity calculation device (100) that calculates seismic intensity on the surface based on the magnitude and epicenter location information of the earthquake collected by the meteorological agency (10). The disclosed system (1) transmits an earthquake early warning generated using the seismic intensity calculated by the seismic intensity calculation device (100) to a user terminal (11) so as to reduce damage caused by the earthquake in advance.

[0031] Specifically, the Meteorological Agency (10) identifies not only the magnitude of the earthquake but also the location where the earthquake occurred, i.e., the location of the epicenter, from seismic waves recorded at earthquake observation stations installed in each region. Here, the location information of the epicenter includes latitude and longitude.

[0032] The Korea Meteorological Administration (10) transmits information regarding the magnitude and location of the epicenter of an earthquake to the seismic intensity calculation device (100). The seismic intensity calculation device (100) calculates the seismic intensity on the surface based on the received information regarding the magnitude and location of the epicenter of the earthquake, and based on this, produces the maximum seismic intensity and seismic intensity distribution map for each administrative district. The seismic intensity calculation device (100) transmits the calculated maximum seismic intensity and seismic intensity distribution map back to the Korea Meteorological Administration (10).

[0033] The magnitude calculation device (100) and the user terminal (11) may be implemented as a computer or portable terminal capable of connecting to the meteorological agency (10) or a communication network (2). Here, the computer may include, for example, a laptop, desktop, laptop, tablet PC, slate PC, etc. equipped with a web browser, and the portable terminal may include, for example, all types of handheld-based wireless communication devices such as PCS (Personal Communication System), GSM (Global System for Mobile communications), PDC (Personal Digital Cellular), PHS (Personal Handyphone System), PDA (Personal Digital Assistant), IMT (International Mobile Telecommunication)-2000, CDMA (Code Division Multiple Access)-2000, W-CDMA (W-Code Division Multiple Access), WiBro (Wireless Broadband Internet) terminal, smartphone, etc., as wireless communication devices that ensure portability and mobility.

[0034] The meteorological agency (10) can warn the user terminal (11) of the occurrence of an earthquake through a communication network (2) and transmit the magnitude of the calculated earthquake, thereby transmitting an alarm or information to minimize damage that the user of the user terminal (11) may suffer from the earthquake.

[0035] Meanwhile, the conventional method for calculating seismic intensity was performed using a grid in which the inland area was divided into 0.1° units to speed up the seismic intensity calculation. However, this method had the problem of not matching the administrative units of the inland area, and in the event of an earthquake, it was not possible to provide accurate maximum seismic intensity or seismic intensity distribution maps to each local government. Accordingly, the disclosed seismic intensity calculation method calculates the seismic intensity at the surface at the same or faster speed than the conventional method (within about 2 seconds after receiving the location information of the earthquake) even though the inland area is divided into 0.01° units, thereby ensuring rapid warning from the Meteorological Agency (10) or local governments and rapid evacuation and safety of users. A detailed explanation of the disclosed seismic intensity calculation method will be described later through other drawings below.

[0036] FIG. 2 is a control block diagram of a seismic calculation device (100) that performs the disclosed seismic calculation method.

[0037] Referring to FIG. 2, the seismic intensity calculation device (100) may include an input unit (101) that receives an input command from a user who wishes to calculate seismic intensity, a communication unit (102) that receives information on the magnitude and location of the epicenter of an earthquake from a communication network (2), a processor (103) that calculates seismic intensity on the surface of the earth based on the received information on the magnitude and location of the epicenter of an earthquake, a database (104) that stores a first table, a second table, a third table, and data files necessary for calculating seismic intensity, and an output unit (105) that displays the calculated seismic intensity and a generated seismic intensity distribution map, or outputs various user interfaces necessary for calculating seismic intensity.

[0038] Specifically, the input unit (101) may include hardware devices such as a keyboard, mouse, switch, pedal, trackball, various levers, handle, or stick for receiving various input commands required for progress calculation. Additionally, the input unit (101) may include a GUI (Graphical User Interface), i.e., a software device such as a touch pad for receiving input commands. The touch pad may be implemented as a touch screen panel (TSP) and form a layered structure with the output unit (105).

[0039] The communication unit (102) may include at least one of a short-range communication module, a wired communication module, and a wireless communication module that communicates with a communication network (2).

[0040] The short-range communication module may include various short-range communication modules that transmit and receive signals using a wireless communication network at short range, such as a Bluetooth module, an infrared communication module, an RFID (Radio Frequency Identification) communication module, a WLAN (Wireless Local Access Network) communication module, an NFC communication module, and a Zigbee communication module.

[0041] Wired communication modules may include various wired communication modules such as Controller Area Network (CAN) communication modules, Local Area Network (LAN) modules, Wide Area Network (WAN) modules, or Value Added Network (VAN) modules, as well as various cable communication modules such as Universal Serial Bus (USB), High Definition Multimedia Interface (HDMI), Digital Visual Interface (DVI), RS-232 (recommended standard 232), power line communication, or plain old telephone service (POTS).

[0042] In addition to Wi-Fi modules and WiBro (Wireless broadband) modules, the wireless communication module may include wireless communication modules that support various wireless communication methods such as GSM (global System for Mobile Communication), CDMA (Code Division Multiple Access), WCDMA (Wideband Code Division Multiple Access), UMTS (universal mobile telecommunications system), TDMA (Time Division Multiple Access), and LTE (Long Term Evolution).

[0043] The processor (103) is configured to control the entire seismic intensity calculation device (100), and the disclosed seismic intensity calculation method can be programmed and implemented through the processor (103). Specifically, the processor (103) calls a lookup table (hereinafter referred to as "table") that is stored separately in the database (104), and in particular, can quickly determine the epicentral distance to all inland grids according to the location information of the epicenter through a large-capacity data file stored in the database (104). In addition, the processor (103) determines the maximum bedrock acceleration calculated in advance in the second table according to the magnitude of the received earthquake and the determined epicentral distance, calculates the site amplification ratio by applying the maximum bedrock acceleration determined based on the third table containing multiple coefficient values, and calculates the seismic intensity on the surface of the earth due to the earthquake based on the calculated site amplification ratio.

[0044] A processor (103) may refer to a data processing device embedded in hardware having a physically structured circuit to perform a function expressed by code or instructions included in a program. Examples of such data processing devices embedded in hardware may include, but are not limited to, processing devices such as a microprocessor, a central processing unit (CPU), a processor core, a multiprocessor, an application-specific integrated circuit (ASIC), a field programmable gate array (FPGA), a Graphics Processing Unit (GPU), and a neural Processing Unit (NPU). The progress calculation device (100) may include at least one processor (103) in hardware.

[0045] The database (104) may be implemented in at least one of a non-volatile memory device such as a cache, ROM (Read Only Memory), PROM (Programmable ROM), EPROM (Erasable Programmable ROM), EEPROM (Electrically Erasable Programmable ROM), and flash memory, a volatile memory device such as RAM (Random Access Memory), or a storage medium such as a hard disk drive (HDD) and CD-ROM, but is not limited thereto. The database (104) may be a memory implemented as a separate chip from the aforementioned processor (103), or it may be implemented as a single chip with the processor (103).

[0046] The output section (105) is configured to display the calculated progress or progress distribution and may be provided with a Digital Light Processing (DLP) panel, a Plasma Display Panel, a Liquid Crystal Display (LCD) panel, an Electro Luminescence (EL) panel, an Electrophoretic Display (EPD) panel, an Electrochromic Display (ECD) panel, a Light Emitting Diode (LED) panel, or an Organic Light Emitting Diode (OLED) panel, but is not limited thereto.

[0047] Meanwhile, the magnitude calculation device (100) in FIG. 2 may additionally include other configurations other than the aforementioned configuration, and the configuration mentioned in FIG. 2 may be omitted as needed.

[0048] FIG. 3 is a flowchart illustrating a disclosed progress calculation method. FIG. 4 is a table illustrating a first table, FIG. 5 is a diagram illustrating a data file, FIG. 6 is a table illustrating a second table, and FIG. 7 is a table illustrating a third table. To avoid redundant descriptions, they are described together below.

[0049] Referring first to FIG. 3, the processor (103) of the disclosed seismic intensity calculation method, or the seismic intensity calculation device (100) performing the seismic intensity calculation method, receives information on the magnitude of the earthquake and the location of the epicenter (110).

[0050] Here, the processor (103) receives information on the magnitude and location of the epicenter of an earthquake from an external source, for example, the processor (103) can receive information on the magnitude and location of the epicenter of an earthquake through various entities such as a meteorological agency (10), an earthquake observation station, or a communication network (2).

[0051] The processor (103) determines the ID of the earthquake classified according to location information based on the first table (location table, see FIG. 4) (120).

[0052] Here, the first table is a lookup table that pre-identifies an ID (POID) that distinguishes earthquakes according to the location information (latitude and longitude) of the epicenter. That is, the disclosed seismic intensity calculation method pre-identifies an earthquake ID for each epicenter where an earthquake may occur and pre-stores it in the database (104) as the first table.

[0053] When the processor (103) receives location information of the epicenter, that is, the longitude and latitude information of the epicenter, it determines the ID of the earthquake according to the received location information by referring to the first table.

[0054] Referring again to FIG. 3, the processor (103) loads a data file (see FIG. 5) according to the ID of the determined earthquake (130).

[0055] Here, the data file refers to a lookup table prepared in multiple folders by organizing earthquake IDs. That is, the data file is divided into multiple files for each earthquake ID included in the first table and stored in advance in the database (104). The processor (103) loads the data file corresponding to the earthquake ID determined in step 120 among the multiple folders stored in advance.

[0056] Meanwhile, the data file illustrated in Fig. 5 is an example of the epicenter distance according to the location of each region (inland). That is, EVID is the data labeling, and {'P009190': 12.27, 'P012945': 17.82, 'P015198'...} are the epicenter distances (Repi) according to each grid location (POID).

[0057] The processor (103) determines the epicentral distance (GRID: Repi) to all inland grids according to the location information of the epicenter based on the loaded data file (140).

[0058] A data file distinguished for each earthquake ID contains data that pre-stores the epicentral distance from the location of the hypocenter that may affect the country to each grid inland. Here, each grid inland is divided into 0.01° units, and the data file pre-stores the epicentral distance from the location of the hypocenter that may occur to each grid as an N-dimensional vector value. Therefore, the disclosed seismic intensity calculation method can have the same effect as calculating all epicentral distances to each grid based on the received hypocenter location information by loading all data files.

[0059] The processor (103) is the second table (PGA rock Based on Table (see Fig. 6), the maximum bedrock acceleration according to the epicentral distance is determined (141).

[0060] Here, the maximum bedrock acceleration (PGA) rock ) is m / s 2It represents the maximum acceleration of bedrock movement when an earthquake occurs, in units of . The processor (103) can quickly determine the maximum acceleration of bedrock in each grid according to the epicentral distance (Repi) determined from the second table and the magnitude (ML) of the received earthquake.

[0061] The processor (103) calculates the site amplification ratio with the maximum bedrock acceleration applied based on the third table (AF Table, see FIG. 7) (142).

[0062] The third table includes coefficient values ​​(SiteAmp_C0, SiteAmp_C1, SiteAmp_C2) required for calculating the site amplification ratio. The processor (103) determines coefficient values ​​that are distinguished in advance according to the ID of each grid. Once the coefficient values ​​are determined by the third table, the processor (103) calculates the site amplification ratio (AF) through the following mathematical formula 1. site ) can be produced (14).

[0063]

[0064] Here, C o , C1, C2 represent SiteAmp_C0, SiteAmp_C1, and SiteAmp_C2 included in the third table.

[0065] The processor (103) calculates the seismic intensity on the surface of the earth based on the calculated site amplification ratio (150).

[0066] Specifically, the processor (103) calculates the site amplification ratio (AF site By substituting ) into Equation 2 below, finally the maximum acceleration at the surface (PGA surface ) produces (15).

[0067]

[0068] Maximum bedrock acceleration (PGA) rock ) is the bedrock maximum acceleration determined in step 141, and the processor (103) calculates the maximum acceleration at the surface (PGA surface) is applied to the seismic conversion formula to finally calculate the seismic intensity at the surface.

[0069] Meanwhile, the disclosed processor (103) can perform updates to the first table, data file, second table, and third table included in the database (104) by a user using the disclosed progress calculation method. For example, if the country or inland area to which the disclosed progress calculation method applies changes, the processor (103) can change the data included in each table. In addition, the processor (103) can also update the data included in each table and data file even if external factors, such as changes in administrative districts, occur.

[0070] Figure 8 is a diagram illustrating the effect of the disclosed seismic calculation method.

[0071] Conventional seismic intensity calculation methods calculated the seismic intensity from the epicenter for each grid segment divided into 0.1° units in the inland area. However, while this method can calculate earthquakes at a rapid pace, the warning messages provided by local governments, including the maximum seismic intensity, did not match the actual seismic intensity due to the grid segments being divided only by 0.1°. Therefore, the disclosed seismic intensity calculation method enables seismic intensity calculation for grid segments further subdivided to 0.01°, and allows for seismic intensity calculation and seismic intensity distribution map generation within approximately 2 seconds, similar to the conventional method using 0.1°.

[0072] Referring to FIG. 8, the disclosed seismic intensity calculation method, or the processor (103) executing the seismic intensity calculation method, may take about 0.284 seconds to load each data (130) from the database (104) from the time it receives the seismic magnitude (Magnitude) and epicenter location information (Latitude, Longitude).

[0073] Here, each data includes determining the earthquake ID from a first table, determining the epicentral distance by loading a data file corresponding to the determined earthquake ID, determining the bedrock maximum acceleration from the epicentral distance determined according to a second table, and determining the coefficient value required for calculating the site amplification ratio from a third table.

[0074] When each data is loaded, the processor (103) may take 0.380 seconds to calculate the maximum acceleration (PGA) and seismic intensity at the surface for every grid separated by 0.01° from the epicenter.

[0075] When the progress is calculated for each grid, the processor (103) may take about 0.225 seconds to calculate the maximum progress in the pre-stored administrative districts, i.e., city, county, and district units.

[0076] The processor (103) calculates the maximum intensity for each administrative district and then generates an intensity distribution map based on this (160).

[0077] Specifically, the processor (103) may take about 1.158 seconds to calculate the intensity distribution map based on the maximum intensity calculated for each administrative district, and may take 0.276 seconds to save the produced intensity distribution map as an XML file (170).

[0078] Each execution time mentioned in Fig. 8 is the CPU Intel TM This was done after setting the i7 and RAM to 16GB. If the specifications of the progress calculation device (100) to which the disclosed progress calculation method is applied are high performance, the time required for progress calculation may be further shortened.

[0079] Meanwhile, the disclosed embodiments may be implemented in the form of a recording medium that stores instructions executable by a computer. The instructions may be stored in the form of program code and, when executed by a processor, may generate a program module to perform the operation of the disclosed embodiments. The recording medium may be implemented as a computer-readable recording medium.

[0080] Computer-readable recording media include all types of recording media that store instructions that can be decoded by a computer. Examples include ROM (read-only memory), RAM (random access memory), magnetic tape, magnetic disk, flash memory, optical data storage devices, etc.

[0081] As described above, the disclosed embodiments have been explained with reference to the attached drawings. Those skilled in the art will understand that the present invention may be practiced in forms different from the disclosed embodiments without changing the technical spirit or essential features of the invention. The disclosed embodiments are illustrative and should not be interpreted restrictively. Explanation of the symbols

[0082] 1: System 2: Communication Network 10: Korea Meteorological Administration 11: User terminal 100: Progression calculation device 101: Input section 102: Communications Department 103: Processor 104: Database 105: Output section

Claims

Claim 1 The method comprises: receiving information on the magnitude and location of the epicenter of an earthquake from an external source; determining the epicentral distance based on the received location information of the epicenter based on a first table containing an earthquake ID (IDentification) classified in advance according to the location information and a pre-stored data file; determining the maximum bedrock acceleration based on the determined epicentral distance based on a second table containing a maximum bedrock acceleration calculated in advance according to the magnitude and epicentral distance of the earthquake; calculating a site amplification ratio by applying the determined maximum bedrock acceleration based on a third table containing a plurality of coefficient values; and calculating the seismic intensity at the surface of the earth according to the earthquake based on the calculated site amplification ratio; wherein determining the epicentral distance includes determining the ID of the earthquake based on the first table; and loading a data file corresponding to the determined earthquake ID from a database storing data files classified by the earthquake ID. A method for calculating seismic intensity, comprising: determining the epicentral distance based on the location information of the received epicenter based on the loaded data file; wherein the data file includes epicentral distances from the epicenter included in the earthquake ID to a plurality of grids classified at intervals of 0.01° according to latitude and longitude, and wherein the first table, the second table, and the third table are converted into N-dimensional vector values ​​and stored in the database. Claim 2 delete Claim 3 delete Claim 4 A method for calculating seismic intensity according to claim 1, wherein calculating the site amplification ratio comprises determining at least three coefficient values ​​predetermined according to the IDs of the plurality of grids in the third table, and calculating the site amplification ratio based on the determined coefficient values ​​and the maximum acceleration of the bedrock. Claim 5 A method for calculating seismic intensity according to claim 1, wherein calculating seismic intensity at the surface of the earth comprises calculating a maximum acceleration at the surface of the earth based on the calculated site amplification ratio and the maximum acceleration of the bedrock, and calculating seismic intensity at the surface of the earth due to the earthquake based on the calculated maximum acceleration at the surface of the earth. Claim 6 A method for calculating seismic intensity according to claim 5, further comprising: calculating a maximum seismic value in an administrative district based on the maximum acceleration of the surface; and generating a seismic intensity distribution map based on the maximum acceleration of the surface. Claim 7 delete Claim 8 A communication unit that receives information on the magnitude and location of the epicenter of a generated earthquake from an external source; and a processor that calculates the intensity of the generated earthquake at the surface based on the received magnitude and location of the epicenter of the earthquake. and a database storing a first table, a second table, a third table, and a data file necessary for calculating the intensity; wherein the processor determines the epicentral distance according to the location information of the received epicenter based on the first table and the data file containing an earthquake ID pre-classified according to location information, determines the maximum bedrock acceleration according to the determined epicentral distance based on the second table containing a maximum bedrock acceleration pre-calculated according to the magnitude and epicentral distance of the earthquake, calculates a site amplification ratio by applying the determined maximum bedrock acceleration based on the third table containing a plurality of coefficient values, and calculates the intensity at the surface of the earthquake that occurred continuously based on the calculated site amplification ratio, and the processor determines the ID of the earthquake that occurred based on the first table; and loads a data file corresponding to the determined earthquake ID from a database storing data files classified by the earthquake ID; An earthquake seismic intensity calculation device characterized by determining the epicentral distance based on the location information of the received epicenter based on the loaded data file, wherein the data file includes epicentral distances from the epicenter included in the earthquake ID to a plurality of grids classified at intervals of 0.01° according to latitude and longitude, and wherein the first table, the second table, and the third table are converted into N-dimensional vector values ​​and stored in the database. Claim 9 delete Claim 10 delete Claim 11 In claim 8, the processor determines at least three coefficient values ​​predetermined according to the IDs of the plurality of grids in the third table, and calculates the site amplification ratio based on the determined coefficient values ​​and the bedrock maximum acceleration. Claim 12 In claim 8, the processor calculates a maximum surface acceleration based on the calculated site amplification ratio and the maximum bedrock acceleration, and calculates the seismic intensity based on the calculated maximum surface acceleration. Claim 13 In claim 12, the processor is a seismic calculation device that, after the seismic intensity due to the earthquake is calculated, calculates a maximum seismic value in an administrative district based on the maximum acceleration of the surface and generates a seismic distribution map based on the maximum seismic value based on the maximum acceleration of the surface. Claim 14 delete Claim 15 Pre-storing a first table, a second table, a third table, and a data file necessary for calculating the seismic intensity of an earthquake; receiving information on the magnitude of an earthquake and the location of the epicenter from an external source; determining the epicentral distance according to the received location information of the epicenter based on a first table containing an earthquake ID classified in advance according to the location information and a pre-storing data file; determining the maximum bedrock acceleration according to the determined epicentral distance based on a second table containing a maximum bedrock acceleration calculated in advance according to the magnitude of the earthquake and the epicentral distance; calculating a site amplification ratio by applying the determined maximum bedrock acceleration based on a third table containing a plurality of coefficient values; and calculating the seismic intensity of the earthquake that occurred in the event of an event based on the calculated site amplification ratio; wherein determining the epicentral distance includes determining the ID of the earthquake that occurred based on the first table; and loading a data file corresponding to the determined earthquake ID from a database storing data files classified by the earthquake ID. A computer-readable recording medium storing a program for executing a control method of an electronic device, characterized by comprising: determining the epicentral distance according to the location information of the received epicenter based on the loaded data file; wherein the data file includes epicentral distances from the epicenter included in the earthquake ID to a plurality of grids classified at intervals of 0.01° according to latitude and longitude, and wherein the first table, the second table, and the third table are converted into N-dimensional vector values ​​and stored in the database.

Images