Method for monitoring seismic intensity of ordinary speed railway and determining influence section and related equipment thereof

Abstract

The application discloses a kind of methods for monitoring earthquake intensity of ordinary speed railway and judging influence interval and related equipment thereof, the method comprises: obtaining seismic event information of seismic network;According to the seismic event information, determine whether there is reported intensity information along the line in the corresponding region of ordinary speed railway;If not, draw intensity distribution map according to the preset earthquake intensity attenuation relationship model;If yes, in the case that the reported intensity information has effective instrument intensity value, the intensity distribution map is drawn by updating the earthquake intensity attenuation relationship model using Kriging interpolation algorithm;From the epicenter to the influence range formed by the preset distance parameter, determine the mileage interval intersected with railway line under the encirclement of different intensity contour lines, and publish externally after creating the corresponding earthquake intensity influence distribution map.The application realizes the whole-process automatic processing of earthquake information monitoring and acquisition, intensity estimation, result display and information release, and improves the earthquake influence analysis and work efficiency of ordinary speed railway.

Description

Technical Field

[0001] This invention relates to the field of railway earthquake monitoring technology, and in particular to a method and related equipment for monitoring earthquake intensity and determining the impact zone of conventional railways. Background Technology

[0002] Earthquakes are a major natural disaster threatening railway safety. As of the end of 2022, the total operating mileage of railways in China reached 155,000 kilometers. For high-speed railways, the high speed of trains means that derailments could cause serious casualties under seismic activity, thus requiring a highly real-time earthquake early warning system. In addition to high-speed railways, my country still has over 100,000 kilometers of conventional railways. According to traffic control rules, conventional railways in the vicinity must be suspended after an earthquake until maintenance personnel conduct on-line inspections. Compared to high-speed railways, conventional railways are more numerous and denser. Even a small earthquake can lead to the suspension of multiple nearby conventional railway lines for inspection; currently, the focus is primarily on post-earthquake response and rapid recovery.

[0003] Currently, China mainly uses the seismic intensity attenuation formula to calculate the distribution of the seismic impact field to quickly assess the impact of earthquakes. However, conventional railways have large east-west or north-south spans and cover areas with significant differences in geological structure. Using a uniform attenuation law to calculate the seismic intensity distribution results in errors compared to the actual situation, which is not conducive to providing a favorable and rapid decision-making basis for the rapid resumption of railway operations. Summary of the Invention

[0004] This invention provides a method and related equipment for monitoring the seismic intensity and determining the impact zone of conventional railways, which solves the problem that existing technologies rely on manual calculation of the impact of earthquakes on railways. This process is cumbersome and requires a lot of manual effort, making it impossible to quickly determine the scope of action based on basic railway data.

[0005] This invention provides a method for monitoring seismic intensity and determining the impact zone of conventional railways, including:

[0006] Obtain earthquake event information pushed by the seismic network;

[0007] Based on the earthquake event information, determine whether there is reported intensity information along the conventional railway line in the corresponding area;

[0008] If not, then draw an intensity distribution map based on the preset earthquake intensity attenuation relationship model;

[0009] If so, and if the reported intensity information has a valid instrument intensity value, the Kriging interpolation algorithm is used to update the earthquake intensity attenuation relationship model and draw an intensity distribution map.

[0010] Within the influence range formed by the distance parameter from the epicenter, based on the railway line that intersects with the intensity contour lines of the intensity distribution map, the first location point closest to the epicenter on the railway line, the distance value, the seismic intensity value corresponding to the first location point, and the mileage intervals that intersect with the railway line under different intensity contour lines are calculated, and the corresponding seismic intensity influence distribution map is created and released to the public.

[0011] According to the present invention, a method for monitoring seismic intensity and determining the impact zone of conventional railways includes the following steps: when the reported intensity information has valid instrument intensity values, updating the seismic intensity attenuation relationship model using a Kriging interpolation algorithm and drawing an intensity distribution map, the steps include:

[0012] Obtain the intensity values ​​contained in the intensity information reported by the seismic intensity monitoring equipment deployed along the conventional railway line;

[0013] With the epicenter as the center, within a fan-shaped area at a preset angle, if any of the earthquake intensity monitoring devices reports an intensity value that meets the preset conditions and is a valid instrument intensity value, then it is determined that there is a valid instrument intensity value within the fan-shaped area.

[0014] If, within the fan-shaped area, the intensity values ​​reported by all the earthquake intensity monitoring devices do not meet the preset conditions for valid instrument intensity values, then it is determined that there are no valid instrument intensity values ​​within the fan-shaped area. Virtual intensity monitoring devices are then inserted along the center line of the fan-shaped area at preset intervals, and virtual instrument intensity values ​​are generated.

[0015] Based on the effective instrument intensity value and / or the virtual instrument intensity value, the corresponding intensity distribution map is drawn using the Kriging interpolation algorithm.

[0016] According to the present invention, a method for monitoring seismic intensity and determining the impact zone of conventional railways includes the following steps for determining that the intensity value is a valid instrument intensity value that meets preset conditions:

[0017] Obtain the earthquake event time contained in the earthquake event information;

[0018] The theoretical intensity value is calculated based on the distance between the latitude and longitude coordinates of the earthquake intensity monitoring equipment and the latitude and longitude coordinates of the epicenter;

[0019] When the absolute value of the difference between the intensity value and the theoretical intensity value is less than a preset first threshold, and the difference between the earthquake event time and the initial start time of the earthquake intensity monitoring device is less than a preset second threshold, it is determined to be a valid instrument intensity value that meets the preset conditions.

[0020] According to the present invention, a method for monitoring seismic intensity and determining the impact zone of conventional railways includes the step of drawing a corresponding intensity distribution map using a Kriging interpolation algorithm based on the effective instrument intensity value and / or the virtual instrument intensity value, comprising:

[0021] Construct a set, and use the second location points corresponding to the earthquake intensity monitoring equipment and / or the virtual intensity monitoring equipment that have obtained the effective instrument intensity value and / or the virtual instrument intensity value as elements of the set;

[0022] The spatial geographic location range for which the Kriging interpolation algorithm is used is determined, and each second location point in the set is marked as an intensity point within the spatial geographic location range. The first intensity value of the intensity point is the effective instrument intensity value or the virtual instrument intensity value.

[0023] Calculate the first semivariance between the first intensity values ​​of any two different intensity points;

[0024] Within the spatial geographic area, a grid is divided to obtain equally spaced grid points. For any one of the grid points, the distance between the grid point and any one of the intensity points is calculated, and the second semi-variance between the grid point and any one of the intensity points is calculated based on the distance.

[0025] Based on the first semi-variance and the second semi-variance, calculate the weight coefficient of each grid point for any intensity point, and calculate the second intensity value of each grid point based on the first intensity value of the intensity point and the weight coefficient;

[0026] Connect the grid points with the same second intensity value to obtain the intensity distribution map with intensity contour lines.

[0027] According to the present invention, a method for monitoring seismic intensity and determining the impact zone of a conventional railway includes the step of calculating the mileage intervals that intersect with the railway line under the coverage of different intensity contour lines, comprising:

[0028] For any intensity contour line, the railway line that intersects with the intensity contour line is divided into intervals, and the interval points are taken as discrete points;

[0029] Based on discrete points within the intensity contour lines, determine the first railway length within all discrete points;

[0030] Obtain the intersection point of the railway line and the intensity contour line, and calculate the first distance and the second distance between the discrete points at both ends of the intensity contour line and the nearest intersection point;

[0031] Based on the first distance and the second distance, the mileage interval marked with railway mileage K is generated.

[0032] According to the present invention, a method for monitoring seismic intensity and determining the impact zone of conventional railways includes the following steps: creating and publicly releasing a corresponding seismic intensity impact distribution map.

[0033] Mark the first location point, the distance value, the earthquake intensity value, and the mileage interval on the map interface, and then save the screenshot.

[0034] A textual description is provided for the first location point, the distance value, the seismic intensity value, and the mileage interval;

[0035] The screenshots and text descriptions are saved locally and linked to the website interface. The corresponding earthquake intensity impact distribution map is automatically generated in web page format and then published externally via a URL link.

[0036] The present invention also provides a system for monitoring the seismic intensity and determining the impact zone of a conventional railway. The system includes: a rapid reporting information processing center for seismic intensity of a conventional railway and seismic intensity monitoring equipment arranged along the conventional railway line. The seismic intensity monitoring equipment is communicatively connected to the rapid reporting information processing center for seismic intensity of a conventional railway.

[0037] The earthquake intensity monitoring equipment includes: a data acquisition unit, a data analysis unit, and a data transmission unit, wherein,

[0038] The data acquisition unit is used for real-time monitoring and data collection of seismic motion along conventional railway lines.

[0039] The data analysis unit is used to generate intensity information, state information, and waveform information, including intensity values, based on seismic motion monitoring data.

[0040] The data transmission unit is used to upload the intensity information to the conventional railway earthquake intensity rapid reporting information processing center;

[0041] The conventional railway earthquake intensity rapid reporting information processing center includes: an interface front-end server, an information processing server, a geographic information server, a database server, and an equipment management server.

[0042] The device manager is used for remote access and control of the earthquake intensity monitoring device;

[0043] The interface front-end server is used to receive the intensity information reported by the seismic intensity monitoring equipment and the seismic event information pushed by the seismic network.

[0044] The database server is used to provide data storage services, including storing the intensity information, the status information, the waveform information, and the earthquake event information.

[0045] The geographic information server is used to provide geospatial data and earthquake location information along the conventional railway lines.

[0046] The information processing server is equipped with a seismic intensity attenuation relationship model, which is used to draw an intensity distribution map based on the intensity information and the seismic event information, and to calculate the mileage intervals that intersect with conventional railway lines under the different intensity contour lines of the intensity distribution map based on the intensity distribution map, the geospatial data and the earthquake location information, and to publish the corresponding seismic intensity impact distribution map after it is created.

[0047] The present invention also provides a device for monitoring the seismic intensity and determining the impact zone of conventional railways, the device comprising:

[0048] The acquisition module is used to acquire earthquake event information pushed by the seismic network;

[0049] The intensity information determination module is used to determine whether there is reported intensity information along the conventional railway line in the corresponding area based on the earthquake event information;

[0050] The effective intensity value determination module is used to determine whether there is a valid instrument intensity value based on the reported intensity information.

[0051] The intensity distribution map drawing module is used to draw an intensity distribution map according to a preset seismic intensity attenuation relationship model when the intensity information is unavailable, or to update the seismic intensity attenuation relationship model and draw an intensity distribution map using a Kriging interpolation algorithm when the reported intensity information has a valid instrument intensity value.

[0052] The earthquake intensity impact distribution map creation module is used to determine the first location point closest to the epicenter, its distance value, and the earthquake intensity value corresponding to the first location point, based on the railway line intersecting with the intensity contour lines of the intensity distribution map within the influence range formed from the epicenter to a preset distance parameter. It also calculates the mileage intervals that intersect with the railway line under different intensity contour lines and publishes the corresponding earthquake intensity impact distribution map.

[0053] The present 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 when the processor executes the program, it implements the method for monitoring seismic intensity and determining the impact zone of conventional railways as described above.

[0054] The present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the method for monitoring seismic intensity and determining the impact zone of conventional railways as described above.

[0055] This invention provides a method and related equipment for monitoring seismic intensity and determining the impact zone of conventional railways. The method involves acquiring seismic event information pushed by a seismic network; determining whether there is reported intensity information along the conventional railway line within the corresponding area based on the seismic event information; if not, drawing an intensity distribution map based on a preset seismic intensity attenuation relationship model; if so, updating the seismic intensity attenuation relationship model and drawing an intensity distribution map using a Kriging interpolation algorithm, provided the reported intensity information has valid instrument intensity values; within the impact zone formed by the epicenter and preset distance parameters, determining the first location point closest to the epicenter on the railway line, its distance value, and the seismic intensity value corresponding to the first location point, based on the railway line intersecting with the intensity contour lines on the intensity distribution map; calculating the mileage intervals intersecting with the railway line under different intensity contour lines; and then publishing the corresponding seismic intensity impact distribution map. This invention enables timely acquisition of earthquake event information from the national earthquake network's published earthquake catalog and valid instrument intensity values ​​from measured instruments along the railway line after an earthquake. Based on basic railway data, it rapidly determines the K-mileage section of the railway under different intensities. Taking into account the fault characteristics of the epicenter, it quickly generates an earthquake intensity impact distribution map. This map visualizes the earthquake impact data and pushes it to subscribed users. It automates the entire process of earthquake information monitoring and acquisition, intensity estimation, result display, and information dissemination, improving the efficiency of earthquake impact analysis and work on conventional railways. It provides data support for railway management units to conduct real-time and efficient post-earthquake response and provides decision-making basis for the rapid resumption of railway operations. Attached Figure Description

[0056] To more clearly illustrate the technical solutions in this invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of this invention. For those skilled in the art, other drawings can be obtained from these drawings without creative effort.

[0057] Figure 1 This is one of the flowcharts illustrating a method for monitoring seismic intensity and determining the affected area of ​​a conventional railway, provided by an embodiment of the present invention.

[0058] Figure 2 A system architecture diagram of a conventional railway seismic intensity monitoring system provided in an embodiment of the present invention;

[0059] Figure 3An information flow diagram of a conventional railway seismic intensity monitoring system provided in an embodiment of the present invention;

[0060] Figure 4 This is the second flowchart illustrating a method for monitoring seismic intensity and determining the affected area of ​​a conventional railway, as provided in an embodiment of the present invention.

[0061] Figure 5 This is a schematic diagram illustrating the calculation of railway mileage K intervals under the enclosure of intensity contour lines, as provided in an embodiment of the present invention.

[0062] Figure 6 This is a schematic diagram of the structure of a conventional railway seismic intensity monitoring and impact zone determination device provided in an embodiment of the present invention;

[0063] Figure 7 This is a schematic diagram of the structure of an electronic device provided in an embodiment of the present invention.

[0064] Figure label:

[0065] 21: Acquisition Module; 22: Intensity Information Determination Module; 23: Effective Intensity Value Judgment Module; 24: Intensity Distribution Map Drawing Module; 25: Seismic Intensity Influence Distribution Map Creation Module. Detailed Implementation

[0066] To make the objectives, technical solutions, and advantages of this invention clearer, the invention will be described clearly and completely below with reference to specific embodiments and accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of this invention, and not all of them. All other embodiments obtained by those skilled in the art based on the embodiments of this invention without creative effort are within the scope of protection of this invention.

[0067] It should be noted that those skilled in the art will understand, explicitly and implicitly, that the embodiments described in this invention can be combined with other embodiments without conflict. Unless otherwise defined, the technical or scientific terms used in this invention should be understood in their ordinary sense by those skilled in the art. The terms "a," "an," "an," "the," etc., used in this invention do not indicate quantity limitation and can represent singular or plural. The terms "comprising," "including," "having," and any variations thereof used in this invention are intended to cover non-exclusive inclusion; the terms "first," "second," "third," etc., used in this invention are merely to distinguish similar objects and do not represent a specific ordering of objects.

[0068] Currently, the analysis of the impact of earthquakes on railway lines typically involves manual verification of the three earthquake elements (earthquake intensity, seismic intensity, and seismic velocity), manual location and calculation, and then delineation of the impact radius using offline maps and railway distribution maps. The covered railway mileage is then estimated based on the affected field. The data analysis process first requires manually inputting the earthquake catalog to locate the epicenter, then manually drawing a range at a certain distance according to regulations to find intersections with railway lines and estimate railway mileage. To understand the magnitude of the impact on conventional railways, it is necessary to find the nearest railway line to the epicenter and manually collect data on the map. The distance is then calculated based on the latitude and longitude coordinates of the epicenter and the nearest point, and finally, the attenuation relationship on the Chinese zoning map is used to calculate the acceleration value at the nearest point. This entire calculation process is cumbersome and requires significant manual effort, making it impossible to automate or conduct efficient post-earthquake response and recovery work in real time. Therefore, a method for monitoring seismic intensity and determining the impact zone for conventional railways is needed to improve the analysis of earthquake impacts on conventional railways and provide a basis for decision-making regarding the rapid resumption of railway operations.

[0069] The following detailed description, in conjunction with the accompanying drawings, of a method and equipment for monitoring seismic intensity and determining the impact zone of conventional railways provided by the present invention, through specific embodiments and application scenarios, will be provided in detail.

[0070] Example 1

[0071] Reference Figure 1 As shown in this embodiment, the method for monitoring seismic intensity and determining the impact zone of conventional railways includes:

[0072] Step S1: Obtain earthquake event information pushed by the seismic network;

[0073] Step S2: Determine whether there is reported intensity information along the conventional railway line in the corresponding area based on the earthquake event information;

[0074] If not, proceed to step S3: Draw an intensity distribution map based on the preset earthquake intensity attenuation relationship model;

[0075] If so, proceed to step S4: if the reported intensity information has a valid instrument intensity value, use the Kriging interpolation algorithm to update the seismic intensity attenuation relationship model and draw an intensity distribution map.

[0076] Step S5: Within the influence range formed by the distance parameters from the epicenter to the preset distance parameters, based on the railway line that intersects with the intensity contour lines of the intensity distribution map, determine the first location point closest to the epicenter on the railway line, the distance value, the seismic intensity value corresponding to the first location point, and calculate the mileage intervals that intersect with the railway line under different intensity contour lines. Then, release the corresponding seismic intensity influence distribution map to the public.

[0077] In this embodiment, step S4 specifically includes:

[0078] Step S41: Obtain the intensity values ​​contained in the intensity information reported by the seismic intensity monitoring equipment deployed along the conventional railway line;

[0079] Step S42: With the epicenter as the center, within a fan-shaped area at a preset angle, if any seismic intensity monitoring device reports an intensity value that meets the preset conditions and is a valid instrument intensity value, then it is determined that there is a valid instrument intensity value within the fan-shaped area.

[0080] Step S43: If the intensity values ​​reported by all earthquake intensity monitoring devices within the sector area do not meet the preset conditions for valid instrument intensity values, then it is determined that there are no valid instrument intensity values ​​within the sector area. Virtual intensity monitoring devices are then inserted along the center line of the sector area at preset intervals, and virtual instrument intensity values ​​are generated.

[0081] Step S44: Draw the corresponding intensity distribution map based on the effective instrument intensity value and / or virtual instrument intensity value.

[0082] In this embodiment, step S42, which involves determining the intensity value as a valid instrument intensity value that meets preset conditions, includes:

[0083] Step S421: Obtain the earthquake event time contained in the earthquake event information;

[0084] Step S422: Calculate the theoretical intensity value based on the distance between the latitude and longitude coordinates of the earthquake intensity monitoring equipment and the latitude and longitude coordinates of the epicenter;

[0085] Step S423: When the absolute value of the difference between the intensity value and the theoretical intensity value is less than the preset first threshold, and the difference between the earthquake event time and the initial movement time of the earthquake intensity monitoring equipment is less than the preset second threshold, it is determined to be a valid instrument intensity value that meets the preset conditions.

[0086] In this embodiment, step S44 specifically includes:

[0087] Step S441: Construct a set, taking the second location points corresponding to the seismic intensity monitoring equipment and / or virtual intensity monitoring equipment that have obtained valid instrument intensity values ​​and / or virtual instrument intensity values ​​as elements of the set;

[0088] Step S442: Determine the spatial geographic range using the Kriging interpolation algorithm, and mark each second location point in the set as an intensity point within the spatial geographic range. The first intensity value of the intensity point is either the effective instrument intensity value or the virtual instrument intensity value.

[0089] Step S443: Calculate the first semivariance between the first intensity values ​​of any two different intensity points;

[0090] Step S444: Divide the area into grids within the spatial geographic location range to obtain grid points with equal spacing. For any grid point, calculate the distance between the grid point and any intensity point, and calculate the second semivariance between the grid point and any intensity point based on the distance.

[0091] Step S445: Calculate the weight coefficient of each grid point for any intensity point based on the first semi-variance and the second semi-variance, and calculate the second intensity value of each grid point based on the first intensity value and the weight coefficient of the intensity point;

[0092] Step S446: Connect the grid points with the same second intensity value on the grid to obtain an intensity distribution map with intensity contour lines.

[0093] In this embodiment, step S5, which involves calculating the mileage intervals that intersect with the railway line under the coverage of contour lines of different intensities, includes:

[0094] Step S51: For any intensity contour line, divide the railway line that intersects with the intensity contour line into intervals, and take the interval points as discrete points;

[0095] Step S52: Based on the discrete points within the intensity contour lines, determine the first railway length within all discrete points;

[0096] Step S53: Obtain the intersection point of the railway line and the intensity contour line, and calculate the first distance and the second distance between the discrete points at both ends of the intensity contour line and the nearest intersection point;

[0097] Step S54: Based on the first distance and the second distance, generate the mileage interval marked with railway mileage K.

[0098] In this embodiment, step S5, which involves creating and publicly releasing a corresponding seismic intensity impact distribution map, includes:

[0099] Step S55: Mark the first location point, distance value, seismic intensity value, and mileage interval on the map interface, then take a screenshot and save it;

[0100] Step S56: Create a textual description of the first location point, distance value, seismic intensity value, and mileage interval;

[0101] Step S57: Save the screenshots and text descriptions locally and link them to the website interface. The corresponding earthquake intensity impact distribution map will be automatically generated in web page format and then published externally via a URL link.

[0102] In this embodiment, the specific application scenario of this method is based on a conventional railway seismic intensity monitoring system, the system architecture of which is as follows: Figure 2As shown. The seismic intensity monitoring system for conventional railways includes:

[0103] The conventional railway earthquake intensity rapid reporting information processing center and the earthquake intensity monitoring equipment deployed along the conventional railway line are connected in communication with the conventional railway earthquake intensity rapid reporting information processing center.

[0104] As an example, the conventional railway earthquake intensity rapid reporting information processing center includes one information processing server, two database servers, one geographic information server, one interface front-end server, one equipment management server, and several gateways. Several earthquake intensity monitoring devices are deployed along the conventional railway line. These devices support VPDN technology and can establish virtual communication links with the conventional railway earthquake intensity rapid reporting information processing center to achieve data communication.

[0105] Earthquake intensity monitoring equipment includes: a data acquisition unit, a data analysis unit, and a data transmission unit.

[0106] The data acquisition unit is used for real-time monitoring and data collection of seismic motion along conventional railway lines.

[0107] The data analysis unit is used to generate intensity information, state information, and waveform information, including intensity values, based on seismic ground motion monitoring data.

[0108] The data transmission unit is used to upload intensity information to the conventional railway earthquake intensity rapid reporting information processing center;

[0109] Specifically, the seismic intensity monitoring equipment includes a data acquisition unit, a data analysis unit, a data transmission unit, and a power supply unit. The data acquisition unit is used for real-time monitoring and acquisition of seismic ground motion data along conventional railway lines. The data analysis unit generates intensity information, waveform information, and status information containing intensity values ​​based on the seismic ground motion monitoring data. This information is then uploaded to the interface front-end server of the conventional railway seismic intensity rapid reporting information processing center via the data transmission unit. The seismic intensity monitoring equipment is deployed along conventional railway lines at certain intervals. Due to the long railway lines and wide distribution of the equipment, the data transmission unit uses a 4G wireless router to establish communication links, sending the various information processed and generated by the data analysis unit to the interface front-end server with a public IP address.

[0110] Meanwhile, the data transmission unit supports VPDN technology, which can establish a virtual communication link with the equipment management server in the conventional railway earthquake intensity rapid reporting information processing center, enabling the equipment management server to remotely access and control the earthquake intensity monitoring equipment.

[0111] The power supply unit adopts two power supply methods: solar power and national grid power, to meet the equipment's power needs in the field.

[0112] The conventional railway earthquake intensity rapid reporting information processing center includes: an interface front-end server, an information processing server, a geographic information server, a database server, and an equipment management server.

[0113] Device Manager is used for remote access and control of earthquake intensity monitoring equipment;

[0114] The interface front-end server is used to receive intensity information reported by earthquake intensity monitoring equipment and earthquake information pushed by the seismic network.

[0115] The database server is used to provide data storage services, including storing intensity information, status information, waveform information, and earthquake information.

[0116] Geographic information servers are used to provide geospatial data and earthquake location information along conventional railway lines.

[0117] The information processing server is equipped with a seismic intensity attenuation relationship model, which is used to draw intensity distribution maps based on intensity information and seismic event information, and to calculate the mileage intervals that intersect with conventional railway lines under the different intensity contour lines of the intensity distribution map based on the intensity distribution map, geospatial data and earthquake location information, and to publish the corresponding seismic intensity impact distribution map after it is created.

[0118] Specifically, the conventional railway earthquake intensity rapid reporting information processing center includes: an interface front-end server, an information processing server, a database server, a geographic information server, an equipment management server, and a gateway. The interface front-end server receives intensity information reported by earthquake intensity monitoring equipment along the conventional railway line and earthquake information pushed by the National Seismic Network, pushes it to the information processing server through a secure platform, and provides internet services. The information processing server performs intensity analysis and estimation. The equipment management server implements functions such as equipment status monitoring and equipment management. The database server provides data storage services, including intensity information, status information, and waveform information reported by earthquake intensity monitoring equipment. The geographic information server provides geospatial data and earthquake location information along the conventional railway line. Its information flow diagram is as follows: Figure 3 As shown.

[0119] See Figure 4 As shown, this paper details the methods for monitoring seismic intensity and determining the impact zone of conventional railways, using the conventional railway seismic intensity monitoring system. The specific steps are as follows:

[0120] 1) The interface front-end server captures earthquake event information pushed by the National Earthquake Network and pushes it to the information processing server through the security platform. The information processing server treats it as an earthquake event.

[0121] 2) The information processing server determines whether the seismic intensity monitoring equipment along the conventional railway line in the corresponding area has reported intensity information based on the earthquake information event;

[0122] If not, proceed to step 3; if yes, proceed to step 4.

[0123] 3) The information processing server uses an elliptical model as the earthquake intensity attenuation relationship model to calculate the intensity distribution. The ellipse angle is the tangent angle between the epicenter and the nearest point on the earthquake fault within 10 kilometers. The calculation formula of the model is as follows:

[0124] I = A + BM + Clg(R + R0) (1)

[0125] In the formula, A, B, C, and R0 are regression coefficients; M is the surface wave magnitude; I is the earthquake intensity; and R is the epicentral distance.

[0126] The values ​​of regression coefficients A, B, C, and R0 for different regional zones are shown in Table 1:

[0127] Table 1 Regression coefficients for earthquake intensity attenuation in each region

[0128]

[0129] 4) The information processing server receives the intensity information reported by the seismic intensity monitoring equipment deployed along the conventional railway line. Taking the epicenter as the center and due north as the 0-degree direction, it divides the area into fan-shaped regions of 90 degrees ([0°, 90°), [90°, 180°), [180°, 270°), [270°, 360°)) to determine whether there are valid instrument intensity values. The determination of valid instrument intensity values ​​is based on the following two conditions (a) and (b):

[0130] (a) Determine whether the deviation between the intensity value reported by the earthquake intensity monitoring device and the theoretical intensity value is less than a preset first threshold. Calculate the distance R based on the latitude and longitude coordinates of the earthquake intensity monitoring device and the epicenter, and substitute it into equation (1) to calculate the theoretical intensity value. If the absolute value of the difference between the actual intensity value generated by the earthquake intensity monitoring device and the theoretical intensity value is within the first threshold range, then the device is considered a valid earthquake intensity monitoring device, and the intensity value it generates is a valid instrument intensity value. Therefore, the intensity value generated by the earthquake intensity monitoring device is used when drawing the intensity distribution map. If the deviation is not within the preset first threshold range, then the intensity value generated by the device is considered invalid, and the intensity value generated by the device is not used when drawing the intensity distribution map.

[0131] (b) Obtain the earthquake event time contained in the earthquake information pushed by the National Earthquake Network and determine whether the deviation range between the time of the earthquake event and the initial motion time triggered by the earthquake intensity monitoring equipment is less than a preset second threshold. If the deviation is within the second threshold range, the intensity value generated by the earthquake intensity monitoring equipment is deemed valid; if the deviation is not within the second threshold range, it is not included in the intensity distribution map calculation.

[0132] The intensity values ​​reported by earthquake intensity monitoring equipment can only be considered valid instrument intensity values ​​if they meet the above two conditions.

[0133] 5) If any seismic intensity monitoring device within the sector reports an intensity value that meets the conditions for a valid instrument intensity value, then it is determined that there is a valid instrument intensity value within the sector, and no virtual seismic intensity monitoring device is inserted; if all seismic intensity monitoring devices within the sector report intensity values ​​that do not meet the conditions for a valid instrument intensity value, then it is determined that there is no valid instrument intensity value within the sector, and a virtual intensity monitoring device is inserted every 10km along the center line of the sector (i.e., the center line with angles of 45°, 135°, 225°, and 315°), and the virtual instrument intensity value is calculated using formula (1) based on the insertion location and the distance from the epicenter.

[0134] 6) Based on the set of location points corresponding to the obtained valid instrument intensity values ​​or virtual instrument intensity values, the seismic intensity attenuation relationship model is updated using the Kriging interpolation algorithm, and an intensity distribution map is drawn. The specific steps are as follows:

[0135] ① Substitute the epicenter latitude and longitude and intensity I=0 into equation (1) to obtain the theoretical radius value R of the earthquake's impact coverage area. With the epicenter as the center, select four points (R,R), (-R,R), (-R,-R), and (R,-R) and connect them to obtain a rectangular area as the spatial geographical location range of the Kriging interpolation algorithm.

[0136] ② Define the rectangular area as A and construct a set. The location points corresponding to the seismic intensity monitoring equipment and / or virtual intensity monitoring equipment that have obtained valid instrument intensity values ​​and / or virtual instrument intensity values ​​are used as elements of the set. Mark the location points (x, y) of each known intensity value in the set in A, where x and y are the longitude and latitude of that location point, respectively. Let z... i =z(x i ,y i Let be the first intensity value at location i, where i = 1, 2, 3, ..., b, and n be the number of locations in A;

[0137] Calculate the first semivariance γ between the first intensity values ​​at two different locations. ij Its calculation formula is shown in equation (2):

[0138]

[0139] ③ Calculate the distance d between each pair of locations using equation (3). ij :

[0140] d ij =r0*arccos(sinx) i *sinx j +cosx i *cosx j *cos(y i -y j ))(3)

[0141] ④ Then select the exponential function model to fit the distance d ij With semivariance γ ij The relationship between them:

[0142] γ(d)=C0+C1(1-e -d / a (4)

[0143] Where C0 and C1 are parameters to be determined, and a is a constant.

[0144] ⑤ Divide the area within the Kriging spatial geographic location into a grid, obtaining equally spaced grid points. For each grid point (x0, y0), calculate the distance d from each grid point to the known intensity point i according to equation (3). i0 Then, calculate the second semivariance γ between the grid point and the known intensity point i according to formula (4). i0 .

[0145] ⑥ Calculate the weight coefficient of each grid point according to equation (5):

[0146]

[0147] Where w1,…,w n The weighting coefficient for each grid point relative to known intensity points, where ε is a constant.

[0148] ⑦ Calculate the second intensity value of each grid point according to formula (6):

[0149]

[0150] ⑧ Connect the grid points with the same second intensity value on the grid with a smooth curve to obtain an intensity distribution map with intensity contour lines.

[0151] 7) Based on the distance parameters set by the system, such as 100km, 200km, etc., from the epicenter to the influence range formed by this distance parameter, find the first location point closest to the epicenter on the railway line and calculate the closest distance value between the railway line and the epicenter; determine the seismic intensity value of the first location point on the railway line according to the intensity distribution map; at the same time, it can also be determined which railway line the first location point is located on, and obtain the railway line information by outputting the latitude and longitude coordinates of the first location point and performing reverse address analysis.

[0152] 8) If the intensity contour lines on the intensity distribution map intersect with conventional railway lines, calculate the K-mileage interval where each intensity contour line intersects with the railway line. The calculation method is as follows:

[0153] ① Divide the railway line into 100m curve length intervals and take the interval points as discrete points, labeling them sequentially as K0, K100, K200... For each discrete point, use the ray method to determine whether the point is within the range corresponding to the intensity contour line. That is, draw a horizontal ray from the point. If the number of intersections with the polygon corresponding to the intensity contour line is odd, then the point is within the range corresponding to the intensity contour line; if the number is even, then the point is outside the range corresponding to the intensity contour line.

[0154] ② Form a set of discrete points within the contour line range (assume K100, K200, ..., K1800), such as Figure 5 The circle to the left of the "+" mark in the enlarged image is shown. Calculating the total number of discrete points n in the set, which is 18, we can find the length of the first railway line among all discrete points as (n-1)*100 = (18-1)*100 = 1700 (m).

[0155] ③ Mark the intersection points of the railway line and the intensity contour lines ( Figure 5 (Marked with a "+" dot).

[0156] ④ Calculate the first distance x between the discrete points at both ends of the intensity contour line and the nearest intersection point (e.g., Figure 5 (As shown in the enlarged diagram) and the second distance x' (not shown in the diagram). Then, the second railway length X of the railway line under this intensity contour line can be approximately equal to the first railway length of the railway line within the discrete points plus the first distance x and the second distance x' on both sides, that is, X = (n-1)*100 + x + x'.

[0157] ⑤ Since the discrete point is the point marking the K mileage, based on the first distance x and the second distance x', the railway line mileage interval under this intensity contour line can be generated as the K mileage interval ([K(100-x')—K(1800+x)]).

[0158] ⑥ Similarly, the railway mileage K intervals enclosed by the other intensity isopleths can be obtained.

[0159] 9) After marking the calculated intensity distribution map, the first location point closest to the epicenter on the railway line, the seismic intensity value, and the mileage intervals of the railway line where different intensity contour lines intersect on the system map interface, select an appropriate interface scale to take a screenshot and save it, ensuring that all calculation results are included in the screenshot. Explain the above calculation results in text form, along with the screenshot. Figure 1 Save it locally.

[0160] Examples of textual explanations are as follows:

[0161] "According to the Taiwan Network Center, an earthquake of magnitude xx occurred at xx:xx:xx on xx:xx in xx month xx year, with a focal depth of xx km. The nearest conventional railway line is the xx railway line, located at xx (its longitude is xx, latitude is xx), and the estimated earthquake intensity is xx."

[0162] "The railway lines affected by this earthquake include xx and xx. The mileage section of the xx railway line with an earthquake intensity ≥3 is Kxxx~Kxxx, the mileage section with an earthquake intensity ≥4 is Kxxx~Kxxx, ..., and the mileage section with an intensity ≥8 is Kxxx~Kxxx."

[0163] 10) Link the text description and screenshots to the website interface, automatically format and edit them into a webpage to generate the corresponding earthquake intensity impact distribution map. Once completed, push the earthquake intensity impact distribution map to the SMS sending server or relevant public platform as a URL link so that subscribers can receive the information.

[0164] In summary, this embodiment provides a method for monitoring seismic intensity and determining the impact zone of conventional railways. It performs real-time seismic intensity monitoring along the railway line, promptly acquiring seismic event information from the national earthquake network's published earthquake catalog and valid instrument intensity values ​​from measured instruments along the line after an earthquake. This generates an intensity distribution map with intensity contour lines. Based on railway baseline data, it quickly determines the railway mileage K intervals enclosed by different intensity contour lines. Taking into account the fault characteristics of the seismic source, it rapidly generates a seismic intensity impact distribution map marked with several seismic impact analysis data points. This map is then visualized and pushed to subscribed users. This method automates the entire process of seismic information monitoring and acquisition, intensity estimation, result display, and information dissemination, improving the efficiency of seismic impact analysis and work on conventional railways. It provides data support for railway management units to conduct real-time and efficient post-earthquake response and provides a decision-making basis for the rapid resumption of railway operations.

[0165] Example 2

[0166] Based on the same inventive concept as the above method, and referring to... Figure 6As shown, this embodiment provides a device for monitoring seismic intensity and determining the impact zone of conventional railways. The device includes:

[0167] Module 21 is used to acquire earthquake information pushed by the seismic network;

[0168] The intensity information determination module 22 is used to determine whether there is reported intensity information along the conventional railway line within the earthquake damage area based on the earthquake information;

[0169] The effective intensity value judgment module 23 is used to determine whether there is an effective intensity value based on the reported intensity information;

[0170] The intensity distribution map drawing module 24 is used to draw an intensity distribution map according to a preset earthquake intensity attenuation relationship model when there is no intensity information, or to draw an intensity distribution map according to a preset earthquake intensity attenuation relationship model when the intensity information has a valid intensity value.

[0171] The mileage relationship calculation module 25 is used to calculate the mileage relationship of at least one railway line covered by the intensity distribution map in different intensity intervals, and to create and publish the corresponding mileage relationship influence map.

[0172] Module 21 is used to acquire earthquake event information pushed by the seismic network;

[0173] The intensity information determination module 22 is used to determine whether there is reported intensity information along the conventional railway line in the corresponding area based on the earthquake event information;

[0174] The effective intensity value judgment module 23 is used to determine whether there is an effective instrument intensity value based on the reported intensity information;

[0175] The intensity distribution map drawing module 24 is used to draw an intensity distribution map according to a preset seismic intensity attenuation relationship model when the intensity information is unavailable, or to update the seismic intensity attenuation relationship model and draw an intensity distribution map by using a Kriging interpolation algorithm when the reported intensity information has a valid instrument intensity value.

[0176] The earthquake intensity impact distribution map creation module 25 is used to determine the first location point closest to the epicenter, the distance value, and the earthquake intensity value corresponding to the first location point on the railway line intersecting with the intensity contour lines of the intensity distribution map within the influence range formed from the epicenter to the preset distance parameters. It also calculates the mileage intervals that intersect with the railway line under different intensity contour lines and publishes the corresponding earthquake intensity impact distribution map.

[0177] It should be noted that the method for monitoring the seismic intensity and determining the impact zone of conventional railways provided in the method embodiments of the present invention can be executed by a device for monitoring the seismic intensity and determining the impact zone of conventional railways, or by a control module in the device for monitoring the seismic intensity and determining the impact zone of conventional railways for executing the method for monitoring the seismic intensity and determining the impact zone of conventional railways.

[0178] The implementation process of the functions and roles of each module in the above system is detailed in the implementation process of the corresponding steps in the above method. Therefore, relevant parts can be referred to in the description of the method embodiment, and will not be repeated here.

[0179] The system embodiments described above are merely illustrative. For example, the division of modules is only a logical functional division. In actual implementation, there may be other division methods. In the embodiments, each functional module can be integrated into a processor, or each module can be a separate device, or two or more modules can be integrated into a device. Each functional module in each embodiment can be implemented in hardware or in the form of hardware plus software functional units.

[0180] Example 3

[0181] Reference Figure 7 As shown, this embodiment provides an electronic device, which includes a processor 310, a communication interface 320, a memory 330, and a communication bus 340. The processor 310, communication interface 320, and memory 330 communicate with each other via the communication bus 340. The processor 310 can call logical instructions stored in the memory 330. The processor 310 executes the conventional railway seismic intensity monitoring and impact zone determination method described in the above embodiment, which includes:

[0182] Obtain earthquake event information pushed by the seismic network;

[0183] Determine whether there is reported intensity information along the conventional railway line in the corresponding area based on earthquake event information;

[0184] If not, then draw an intensity distribution map based on the preset earthquake intensity attenuation relationship model;

[0185] If so, and provided that the reported intensity information has valid instrument intensity values, the Kriging interpolation algorithm is used to update the earthquake intensity attenuation relationship model and draw an intensity distribution map.

[0186] Within the influence range formed by the distance parameters from the epicenter, based on the railway line that intersects with the intensity contour lines on the intensity distribution map, the first location point closest to the epicenter on the railway line, the distance value, the seismic intensity value corresponding to the first location point, and the mileage intervals that intersect with the railway line under different intensity contour lines are calculated, and the corresponding seismic intensity influence distribution map is created and released to the public.

[0187] Furthermore, the logical instructions in the aforementioned memory 330 can be implemented as software functional units and, when sold or used as independent products, can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, essentially, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0188] On the other hand, the present invention also provides a computer program product, which includes a computer program that can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer is able to execute the method for monitoring the seismic intensity and determining the impact zone of conventional railways as described in the above-described method embodiments. The method includes:

[0189] Obtain earthquake event information pushed by the seismic network;

[0190] Determine whether there is reported intensity information along the conventional railway line in the corresponding area based on earthquake event information;

[0191] If not, then draw an intensity distribution map based on the preset earthquake intensity attenuation relationship model;

[0192] If so, and provided that the reported intensity information has valid instrument intensity values, the Kriging interpolation algorithm is used to update the earthquake intensity attenuation relationship model and draw an intensity distribution map.

[0193] Within the influence range formed by the distance parameters from the epicenter, based on the railway line that intersects with the intensity contour lines on the intensity distribution map, the first location point closest to the epicenter on the railway line, the distance value, the seismic intensity value corresponding to the first location point, and the mileage intervals that intersect with the railway line under different intensity contour lines are calculated, and the corresponding seismic intensity influence distribution map is created and released to the public.

[0194] Example 4

[0195] This embodiment provides a non-transitory computer-readable storage medium storing a computer program thereon. When the computer program is executed by a processor, it implements the method for monitoring seismic intensity and determining the impact zone of conventional railways as described in the above-described method embodiment. The method includes:

[0196] Obtain earthquake event information pushed by the seismic network;

[0197] Determine whether there is reported intensity information along the conventional railway line in the corresponding area based on earthquake event information;

[0198] If not, then draw an intensity distribution map based on the preset earthquake intensity attenuation relationship model;

[0199] If so, and provided that the reported intensity information has valid instrument intensity values, the Kriging interpolation algorithm is used to update the earthquake intensity attenuation relationship model and draw an intensity distribution map.

[0200] Within the influence range formed by the distance parameters from the epicenter, based on the railway line that intersects with the intensity contour lines on the intensity distribution map, the first location point closest to the epicenter on the railway line, the distance value, the seismic intensity value corresponding to the first location point, and the mileage intervals that intersect with the railway line under different intensity contour lines are calculated, and the corresponding seismic intensity influence distribution map is created and released to the public.

[0201] The various embodiments in this invention are described in a progressive manner. Similar or identical parts between embodiments can be referred to mutually. Each embodiment focuses on describing the differences from other embodiments. In particular, the device and medium embodiments are relatively simple in description because they are fundamentally similar to the method embodiments; relevant parts can be referred to the descriptions in the method embodiments.

[0202] The devices, media, and methods provided in the embodiments of the present invention are one-to-one correspondences. Therefore, the devices and media also have similar beneficial technical effects as their corresponding methods. Since the beneficial technical effects of the methods have been described in detail above, the beneficial technical effects of the devices and media will not be repeated here.

[0203] It should also be noted that the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process method or product that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process method or product. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process method or product that includes that element.

[0204] The above are merely embodiments of the present invention and are not intended to limit the invention. Although the present invention has been described in detail above with general descriptions and specific embodiments, modifications or improvements can be made to it, which will be obvious to those skilled in the art. Therefore, all such modifications or improvements made without departing from the spirit of the present invention fall within the scope of protection claimed by the present invention.

Claims

1. A method for monitoring seismic intensity and determining the impact zone of conventional railways, characterized in that, include: Obtain earthquake event information pushed by the seismic network; Based on the earthquake event information, determine whether there is reported intensity information along the conventional railway line in the corresponding area; If not, then draw an intensity distribution map based on the preset earthquake intensity attenuation relationship model; If so, and if the reported intensity information has a valid instrument intensity value, the Kriging interpolation algorithm is used to update the earthquake intensity attenuation relationship model and draw an intensity distribution map. Within the influence range formed by the distance parameter from the epicenter to the preset distance, based on the railway line that intersects with the intensity contour lines of the intensity distribution map, determine the first location point on the railway line closest to the epicenter, the closest distance value between the railway line and the epicenter, the seismic intensity value corresponding to the first location point, and calculate the mileage intervals that intersect with the railway line under different intensity contour lines. Then, release the corresponding seismic intensity influence distribution map to the public. The step of calculating the mileage intervals that intersect with the railway line under different intensity isolines includes: for any intensity isoline, dividing the railway line that intersects with the intensity isoline into intervals, and taking the interval points as discrete points; based on the discrete points within the intensity isoline range, determining the first railway length within all discrete points; obtaining the intersection point of the railway line and the intensity isoline, and calculating the first distance and the second distance between the discrete points at both ends within the intensity isoline range and the nearest intersection point; and generating the mileage interval marked with railway mileage K based on the first distance and the second distance.

2. The method for monitoring seismic intensity and determining the impact zone of conventional railways according to claim 1, characterized in that, The step of updating the seismic intensity attenuation relationship model and drawing an intensity distribution map using the Kriging interpolation algorithm when the reported intensity information has valid instrument intensity values ​​includes: Obtain the intensity values ​​contained in the intensity information reported by the seismic intensity monitoring equipment deployed along the conventional railway line; With the epicenter as the center, within a fan-shaped area at a preset angle, if any of the earthquake intensity monitoring devices reports an intensity value that meets the preset conditions and is a valid instrument intensity value, then it is determined that there is a valid instrument intensity value within the fan-shaped area. If, within the fan-shaped area, the intensity values ​​reported by all the earthquake intensity monitoring devices do not meet the preset conditions for valid instrument intensity values, then it is determined that there are no valid instrument intensity values ​​within the fan-shaped area. Virtual intensity monitoring devices are then inserted along the center line of the fan-shaped area at preset intervals, and virtual instrument intensity values ​​are generated. Based on the effective instrument intensity value and / or the virtual instrument intensity value, the corresponding intensity distribution map is drawn using the Kriging interpolation algorithm.

3. The method for monitoring seismic intensity and determining the impact zone of conventional railways according to claim 2, characterized in that, The step of determining that the intensity value is a valid instrument intensity value that meets preset conditions includes: Obtain the earthquake event time contained in the earthquake event information; The theoretical intensity value is calculated based on the distance between the latitude and longitude coordinates of the earthquake intensity monitoring equipment and the latitude and longitude coordinates of the epicenter; When the absolute value of the difference between the intensity value and the theoretical intensity value is less than a preset first threshold, and the difference between the earthquake event time and the initial start time of the earthquake intensity monitoring device is less than a preset second threshold, it is determined to be a valid instrument intensity value that meets the preset conditions.

4. The method for monitoring seismic intensity and determining the impact zone of conventional railways according to claim 2, characterized in that, The step of drawing the corresponding intensity distribution map using the Kriging interpolation algorithm based on the effective instrument intensity value and / or the virtual instrument intensity value includes: Construct a set, and use the second location points corresponding to the earthquake intensity monitoring equipment and / or the virtual intensity monitoring equipment that have obtained the effective instrument intensity value and / or the virtual instrument intensity value as elements of the set; The spatial geographic range for which the Kriging interpolation algorithm is used is determined, and each second location point in the set is marked as an intensity point within the spatial geographic range. The first intensity value of the intensity point is the effective instrument intensity value or the virtual instrument intensity value. Calculate the first semivariance between the first intensity values ​​of any two different intensity points; Within the spatial geographic area, a grid is divided to obtain equally spaced grid points. For any one of the grid points, the distance between the grid point and any one of the intensity points is calculated, and the second semi-variance between the grid point and any one of the intensity points is calculated based on the distance. Based on the first semi-variance and the second semi-variance, calculate the weight coefficient of each grid point for any intensity point, and calculate the second intensity value of each grid point based on the first intensity value of the intensity point and the weight coefficient; Connect the grid points with the same second intensity value to obtain the intensity distribution map with intensity contour lines.

5. The method for monitoring seismic intensity and determining the impact zone of conventional railways according to claim 1, characterized in that, The steps of creating and releasing a corresponding earthquake intensity impact distribution map include: Mark the first location point, the distance value, the earthquake intensity value, and the mileage interval on the map interface, and then save the screenshot. A textual description is provided for the first location point, the distance value, the seismic intensity value, and the mileage interval; The screenshots and text descriptions are saved locally and linked to the website interface. The corresponding earthquake intensity impact distribution map is automatically generated in web page format and then published externally via a URL link.

6. A system for monitoring seismic intensity and determining the impact zone of conventional railways, characterized in that, The system includes: a conventional railway earthquake intensity rapid reporting information processing center and earthquake intensity monitoring equipment deployed along the conventional railway line, wherein the earthquake intensity monitoring equipment is communicatively connected to the conventional railway earthquake intensity rapid reporting information processing center; The earthquake intensity monitoring equipment includes: a data acquisition unit, a data analysis unit, and a data transmission unit, wherein, The data acquisition unit is used for real-time monitoring and data collection of seismic motion along conventional railway lines. The data analysis unit is used to generate intensity information, state information, and waveform information, including intensity values, based on seismic ground motion monitoring data. The data transmission unit is used to upload the intensity information to the conventional railway earthquake intensity rapid reporting information processing center; The conventional railway earthquake intensity rapid reporting information processing center includes: an interface front-end server, an information processing server, a geographic information server, a database server, and an equipment management server. The device manager is used for remote access and control of the earthquake intensity monitoring device; The interface front-end server is used to receive the intensity information reported by the seismic intensity monitoring equipment and the seismic event information pushed by the seismic network. The database server is used to provide data storage services, including storing the intensity information, the status information, the waveform information, and the earthquake event information. The geographic information server is used to provide geospatial data and earthquake location information along the conventional railway lines. The information processing server is equipped with a seismic intensity attenuation relationship model, which is used to draw an intensity distribution map based on the intensity information and the seismic event information, and to calculate the mileage intervals that intersect with conventional railway lines under the different intensity contour lines of the intensity distribution map based on the intensity distribution map, the geospatial data and the earthquake location information, and to publish the corresponding seismic intensity impact distribution map after it is created. The step of calculating the mileage intervals that intersect with the railway line under different intensity isolines includes: for any intensity isoline, dividing the railway line that intersects with the intensity isoline into intervals, and taking the interval points as discrete points; based on the discrete points within the intensity isoline range, determining the first railway length within all discrete points; obtaining the intersection point of the railway line and the intensity isoline, and calculating the first distance and the second distance between the discrete points at both ends within the intensity isoline range and the nearest intersection point; and generating the mileage interval marked with railway mileage K based on the first distance and the second distance.

7. A device for monitoring seismic intensity and determining the impact zone of conventional railways, characterized in that, The device includes: The acquisition module is used to acquire earthquake event information pushed by the seismic network; The intensity information determination module is used to determine whether there is reported intensity information along the conventional railway line in the corresponding area based on the earthquake event information; The effective intensity value determination module is used to determine whether there is a valid instrument intensity value based on the reported intensity information. The intensity distribution map drawing module is used to draw an intensity distribution map according to a preset seismic intensity attenuation relationship model when the intensity information is unavailable, or to update the seismic intensity attenuation relationship model and draw an intensity distribution map using a Kriging interpolation algorithm when the reported intensity information has a valid instrument intensity value. The earthquake intensity impact distribution map creation module is used to determine the first location point closest to the epicenter on the railway line intersecting with the intensity contour lines of the intensity distribution map within the influence range formed from the epicenter to a preset distance parameter, the closest distance value between the railway line and the epicenter, the earthquake intensity value corresponding to the first location point, and to calculate the mileage intervals that intersect with the railway line under different intensity contour lines. After creating the corresponding earthquake intensity impact distribution map, it is released to the public. The step of calculating the mileage intervals that intersect with the railway line under different intensity isolines includes: for any intensity isoline, dividing the railway line that intersects with the intensity isoline into intervals, and taking the interval points as discrete points; based on the discrete points within the intensity isoline range, determining the first railway length within all discrete points; obtaining the intersection point of the railway line and the intensity isoline, and calculating the first distance and the second distance between the discrete points at both ends within the intensity isoline range and the nearest intersection point; and generating the mileage interval marked with railway mileage K based on the first distance and the second distance.

8. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the program, it implements the method for monitoring the seismic intensity and determining the impact zone of conventional railways as described in any one of claims 1-5.

9. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the method for monitoring the seismic intensity and determining the impact zone of conventional railways as described in any one of claims 1-5.

Images