Earth surface dynamics real-time monitoring system

Abstract

The embodiment of the invention discloses a real-time monitoring system for earth surface dynamics, and the system comprises a message and data subsystem which is used for the management of the whole process from the access to the archiving of messages and data, and provides a message exchange and data acquisition platform for other subsystems; the data real-time processing subsystem is used for acquiring seismic data, performing earth surface dynamic parameter real-time inversion, landslide event discrimination and parameter determination, and returning a result to the message and data subsystem; the result comprises an earth surface dynamics real-time inversion result and landslide event information; and the visual display system is used for acquiring the messages and the data to realize visual display. According to the earth surface dynamics real-time monitoring system provided by the embodiment of the invention, earth surface dynamics parameter real-time inversion, landslide event discrimination, parameter determination and the like based on seismic data can be realized, results are visually displayed, and data support can be provided for the aspects of earthquake disaster prevention and control, earthquake prediction and prevention and the like.

Description

Technical Field

[0001] The invention relates to the technical field of surface dynamics research, and in particular to a real-time monitoring system for surface dynamics. Background Art

[0002] Surface dynamics monitoring is an important research and application field, which is mainly used to monitor and analyze the dynamic changes of the earth's surface.

[0003] Surface dynamics monitoring covers a variety of technologies and methods. For example, satellite remote sensing technology can be used to observe changes in surface topography, vegetation coverage, soil moisture, etc. over a large area and for a long time. Through a network of sensors installed on the ground, data such as ground displacement, vibration, and temperature can be obtained in real time.

[0004] The application scope of surface dynamics monitoring is very wide. In the prevention and control of geological disasters, for example, monitoring the surface movement of potential dangerous areas such as landslides and mud-rock flows, issuing early warnings, and protecting people's lives and property. In urban construction, surface dynamics monitoring is carried out on the foundations of large buildings or infrastructure to ensure their stability.

[0005] For example, when building a highway in a mountainous area, surface dynamics monitoring can promptly detect signs of mountain instability caused by construction or natural factors, and take corresponding reinforcement measures. Another example is that in earthquake-prone areas, continuous surface dynamics monitoring helps to understand the laws of crustal movement and provide data support for earthquake prediction and prevention.

[0006] In short, surface dynamics monitoring is of great significance for ensuring the safety of human activities, rationally utilizing natural resources, and deeply understanding the changing mechanisms of the earth's surface. Summary of the invention

[0007] The purpose of the embodiments of the present invention is to provide a real-time monitoring system for surface dynamics, which provides data support for earthquake disaster prevention, earthquake prediction and prevention, etc.

[0008] To achieve the above object, an embodiment of the present invention provides a real-time monitoring system for surface dynamics, comprising:

[0009] The message and data subsystem is responsible for the management of the entire process of message and data from access to archiving, and provides a message exchange and data acquisition platform for other subsystems;

[0010] A real-time data processing subsystem, used to obtain seismic data from the message and data subsystem, perform real-time inversion of surface dynamic parameters, landslide event identification and parameter determination based on the seismic data, and return the results to the message and data subsystem; the results include real-time inversion results of surface dynamics and landslide event information;

[0011] The visual display system is used to obtain messages and data from the message and data subsystem to achieve visual display.

[0012] As a specific implementation of the present application, the message and data subsystem includes:

[0013] A message exchange center is used to manage various time-series messages, including seismic data, real-time inversion results of surface dynamics returned by the real-time data processing subsystem, and landslide event information;

[0014] Data cache center, used to provide seismic data within the time validity range;

[0015] The first functional module is used to obtain seismic data from a seismic real-time data server or a seismic data file and store it in the message exchange center;

[0016] The second functional module is used to obtain seismic data from the message exchange center in real time and store it in a data cache center;

[0017] The third functional module is used to obtain the real-time inversion results of surface dynamics and landslide event information from the message exchange center and archive them.

[0018] As a specific implementation of the present application, the messages and data acquired by the message and data subsystem include access data, and the access data includes continuous waveform data, offline seismic data and historical data;

[0019] The processing of the continuous waveform data is as follows:

[0020] Continuous waveform data is obtained based on two different protocols, NetSeis / IP and SeedLink; the continuous waveform data is accessed in the form of original MSEED data packets;

[0021] Extracting network station information from the original MSEED data packet;

[0022] Using the network station information as the key, the original MSEED data packets are directly stored in different partitions of the message exchange center to achieve distributed storage of data and facilitate parallel computing of subsequent processing modules;

[0023] The processing method for offline seismic data is as follows:

[0024] The MSEED data and SAC data are connected to the message and data subsystem in the form of files, which facilitates the real-time data processing subsystem to calculate surface dynamic parameters, identify landslide events and estimate landslide parameters.

[0025] As a specific implementation of the present application, the access data also includes intermediate calculation results, which are processed as follows:

[0026] The intermediate calculation results are divided into spatial grid points; the inversion results at each of the spatial grid points are event continuous data;

[0027] The inversion results are stored in different partitions of the message exchange center using the grid point identifier as a key, and the storage format is JSON.

[0028] As a specific implementation of the present application, the message and data acquired by the message and data subsystem include the identified landslide event information, which is processed as follows:

[0029] The identified landslide event information is stored in the data cache center after parameter optimization; each landslide event information at least includes event ID, longitude and latitude, estimated volume, inverted force time series, stations and components used in inversion, data frequency bands used in inversion, signal-to-noise ratio of each component observation data, and correlation coefficient between observation data and theoretical data.

[0030] As a specific implementation of the present application, the real-time data processing subsystem includes:

[0031] Initialization module, used to initialize the system, including parameter file reading, spatial grid point division and Green's function calculation;

[0032] A preprocessing module, used to obtain seismic data from the message and data subsystem in real time, perform preprocessing and calculate the signal-to-noise ratio; the preprocessing includes sliding averaging, bandpass filtering and removing instrument response;

[0033] The surface dynamic parameter inversion module is used to traverse all spatial grid points, use the recorded data of the seismic stations associated with the grid points and the corresponding Green's function, and use the least squares method to complete the inversion of surface dynamic parameters in the time domain;

[0034] The landslide event identification module is used to obtain the intermediate results of the inversion from the message and data subsystem, determine whether it exceeds a given threshold, and at the same time, determine whether the event exceeding the threshold is a landslide through the mutual correlation coefficient between the observed data and the theoretical data.

[0035] As a specific implementation of the present application, the initialization module is specifically used for:

[0036] Take the parameter file and divide the space grid points by given parameters;

[0037] Bind the stations within a given distance range around each grid point to the grid point;

[0038] Calculate the Green's function for each grid point;

[0039] Initialize the instrument response at each station by given parameters;

[0040] The station information, grid point information and the binding information between them are stored in the data cache center of the information and data subsystem.

[0041] Furthermore, as a preferred implementation of the present application, the real-time data processing subsystem further includes a parameter optimization module, which is used to:

[0042] For the landslide events that are preliminarily identified, the spatial grid is refined, the Green's function at each grid point is calculated, and the dynamic parameter inversion is completed;

[0043] Traverse the inversion results and select the optimal position;

[0044] After obtaining the results, the empirical formula is used to obtain the important parameters of landslide volume.

[0045] As a specific implementation of the present application, the visual display subsystem includes:

[0046] Visualization display module, used to obtain and display inversion results from the message and data subsystem, and provide a visualization engine for the entire system;

[0047] The human-computer interaction module is used to receive user operations for the selected landslide event to realize re-inversion.

[0048] The visual display module is specifically used for:

[0049] Spatial display of intermediate inversion results: Obtain intermediate inversion results from the message exchange center and display the change of the value of each spatial grid point over time in real time;

[0050] Time series display of intermediate inversion results: For one or more given grid points, the intermediate inversion results are obtained from the message exchange center, and the evolution characteristics of the inversion results over time are displayed in real time in the form of waveforms; for multiple grid points, it has the function of sorting according to different attributes and scaling amplitude and time;

[0051] Final inversion result display: landslide events within a given time range are given in a list format.

[0052] The real-time surface dynamics monitoring system provided by the embodiment of the present invention can realize real-time inversion of surface dynamic parameters based on seismic data, identification of landslide events and parameter determination, etc., and visualize the results, which can provide data support for earthquake disaster prevention and control, earthquake prediction and prevention, etc. BRIEF DESCRIPTION OF THE DRAWINGS

[0053] In order to more clearly illustrate the specific implementation of the present invention or the technical solution in the prior art, the drawings required for use in the specific implementation or the description of the prior art will be briefly introduced below.

[0054] Figure 1 is a structural diagram of a real-time monitoring system for surface dynamics provided by an embodiment of the present invention;

[0055] Figure 2 is the time series diagram of the intermediate inversion results;

[0056] Figure 3 is the final inversion result diagram;

[0057] Figure 4 is the final inversion result - adjustment parameter map. DETAILED DESCRIPTION

[0058] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0059] It should be understood that when used in this specification and the appended claims, the terms "include" and "comprises" indicate the presence of described features, integers, steps, operations, elements and / or components, but do not exclude the presence or addition of one or more other features, integers, steps, operations, elements, components and / or combinations thereof.

[0060] Please refer to Figure 1 , is a real-time monitoring system for surface dynamics provided by an embodiment of the present invention, including a message and data subsystem, a real-time data processing subsystem and a visualization display subsystem.

[0061] The message and data subsystem is responsible for the management of the entire process of message and data from access to archiving, and provides a message exchange and data acquisition platform for other subsystems;

[0062] A real-time data processing subsystem, used to obtain seismic data from the message and data subsystem, perform real-time inversion of surface dynamic parameters, landslide event identification and parameter determination based on the seismic data, and return the results to the message and data subsystem; the results include real-time inversion results of surface dynamics and landslide event information;

[0063] The visual display system is used to obtain messages and data from the message and data subsystem to achieve visual display.

[0064] The following will introduce the composition and functions of the three subsystems respectively:

[0065] 1. Message and Data Subsystem

[0066] The subsystem consists of two centers and three functional modules, namely:

[0067] The message exchange center is responsible for the management of various time-series messages including seismic data;

[0068] Data Cache Center: used to provide earthquake data within the time validity range:

[0069] The first functional module is used to obtain earthquake data from an earthquake real-time data server or an earthquake data file and store it in a message exchange center;

[0070] The second functional module is used to obtain seismic data from the message exchange center in real time and store it in the data cache center;

[0071] The third functional module is used to obtain the real-time inversion results of surface dynamics and landslide event information from the message exchange center and archive them.

[0072] In this embodiment, the messages and data acquired by the message and data subsystem include access data, identified landslide event information, and other program running parameters. The system processes these types of data in different ways, as follows:

[0073] 1.1 Access Data

[0074] (1) Continuous waveform data

[0075] The message and data subsystem obtains continuous waveform data based on two different protocols, NetSeis / IP and SeedLink. To ensure the orderliness of seismic data, the original MSEED data packets are directly stored in different partitions of the message exchange center using the station code as the key, realizing distributed storage of data and facilitating parallel computing of subsequent processing modules.

[0076] At the same time, the data can also be transferred to the data cache center through the data distribution module.

[0077] For continuous waveform data, the message exchange center can realize sequential processing of data, while the data cache center provides data query and random access.

[0078] (2) Offline seismic data

[0079] Regarding offline seismic data, the system connects MSEED data and SAC data into the system in the form of files, completes the calculation of surface dynamic parameters, landslide event identification and landslide parameter estimation, and stores the results in designated files.

[0080] (3) Historical data playback

[0081] For given historical seismic data, the system reads offline seismic data and sends it to the solution service at the original speed or a specified multiple of the speed for real-time solution; the system simulates real-time seismic data stream processing, completes surface dynamic parameter calculation, landslide event identification and landslide parameter estimation, and stores the results in the specified file.

[0082] (4) Intermediate calculation results

[0083] The intermediate calculation results are divided into grid points in space, and the inversion results at each spatial grid point are time-continuous data. The inversion results are stored in different partitions of the message exchange center with the grid point identifier as the key, and the storage format is JSON. Since there is no need to query the intermediate results, there is no need to transfer them to the data cache center.

[0084] 1.2 Identified landslide events

[0085] For the identified landslide events, after parameter optimization, they are stored in the data cache center for future query. The storage information of each event includes at least the event ID, longitude and latitude, estimated volume, inverted force time series, stations and components used for inversion, data frequency bands used for inversion, signal-to-noise ratio of each component observation data, and correlation coefficient between observation data and theoretical data. The updated results after manual interaction are also stored in the data cache center. At the same time, all the above information is archived.

[0086] 1.3 Other parameters of program running

[0087] The station information (network code, station code, channel code, longitude and latitude coordinates, instrument response information), grid information (longitude and latitude, grid code), the association information between the station and the grid, the association information between the grid and the station, etc. are stored in the data cache center in JSON format to facilitate the use of other modules.

[0088] The station list, station coordinates, grid point list, grid point coordinates, etc. are all stored in redis, and other systems can obtain the coordinate information of all stations and grid points through these data.

[0089] It should be noted that the message and data subsystem, as the message and data hub of the entire system, undertakes important functions such as the import and distribution, storage, archiving and caching of messages and data. There is no direct communication between the computing and display modules. The message / data exchange between all modules is carried out through the message and data subsystem; the independence between modules is improved, the coupling between modules is reduced, and the entire system is easy to manage and expand. In order to ensure the stability, reliability, backup, fault tolerance and scalability of the system, this system needs to be supported by big data technology.

[0090] 2. Real-time Data Processing Subsystem

[0091] The subsystem consists of six functional modules, including:

[0092] Initialization module, used for system initialization;

[0093] A calculation module, used for calculating Green's function;

[0094] A preprocessing module, used for preprocessing seismic data;

[0095] The surface dynamic parameter inversion module is used to invert the surface dynamic parameters;

[0096] Landslide event identification module, used for landslide identification based on the intermediate results of inversion;

[0097] The parameter optimization module is used to optimize the inversion parameters for re-inversion.

[0098] 2.1 System Initialization

[0099] When the data real-time processing subsystem is started, it first reads the parameter file, divides the spatial grid points by given parameters, binds the stations within a given distance range around each grid point to the grid point, and calculates the Green's function of each grid point. The instrument response of each station is initialized by given parameters. The station information, grid point information, and the binding information between them are saved in the data cache center.

[0100] 2.2 Green’s function calculation

[0101] The real-time data processing subsystem obtains the velocity model from Crust1.0 for a given landslide and station location, calls the Qseis06 program to calculate the single-force Green's function, and converts the result into a rectangular coordinate system. This function is used to obtain the Green's function of each grid point to its bound station when the system is initialized, and after identifying an event, optimize the inversion result by refining the spatial position.

[0102] 2.3 Seismic data preprocessing

[0103] The real-time data processing subsystem obtains seismic data from the message and data subsystem in real time, and performs preprocessing operations including sliding averaging, bandpass filtering, and instrument response removal. At the same time, it calculates the signal-to-noise ratio of the data to provide a basis for selecting weights for subsequent inversion calculations. Preprocessing is performed in the time domain; therefore, it is necessary to save the intermediate state variables of each step of processing.

[0104] 2.4 Inversion of surface dynamic parameters

[0105] The real-time data processing subsystem traverses all spatial grid points, uses the recorded data of the seismic stations associated with the grid points and the corresponding Green's function, and uses the least squares method to complete the inversion of surface dynamic parameters in the time domain.

[0106] After the inversion is completed, the system converts the results into JSON format and sends them to the message exchange center of the message and data subsystem. At the same time, the correlation coefficient between the observed data and the theoretical data of each record is saved for subsequent evaluation.

[0107] 2.5 Landslide event identification

[0108] The real-time data processing subsystem obtains the intermediate results of the inversion from the message and data subsystem to determine whether they exceed the given threshold. At the same time, the correlation coefficient between the observed data and the theoretical data is used to determine whether the event exceeding the threshold is a landslide. The same event observed at multiple grid points at the same time is merged.

[0109] 2.6 Landslide parameter optimization

[0110] The real-time data processing subsystem refines the spatial grid for the initially determined landslide events, calculates the Green's function at each grid point, and completes the inversion of dynamic parameters. It traverses the inversion results and selects the optimal position. After obtaining the results, the important parameters of the landslide volume are obtained using empirical formulas. The data can be re-inverted by adjusting the landslide time and assigning different weights to seismic stations and components.

[0111] It should be noted that the real-time data processing subsystem undertakes the computing tasks of the entire system. Each of its submodules obtains input data from the message and data subsystem, completes the calculation, and then returns the result to the message and data subsystem. This system is composed of multiple functional modules in a loose mode, which do not directly depend on each other; therefore, it is very easy to add and delete modules. The system needs to be scalable to cope with the growing computing needs. The system can be implemented under the computing framework of real-time processing of big data, or it can be implemented through parallel computing technology.

[0112] 3. Visualization display subsystem

[0113] The subsystem consists of two functional modules:

[0114] Visualization display module, used to obtain and display inversion results from the message and data subsystem, and provide a visualization engine for the entire system;

[0115] The human-computer interaction module is used to receive user operations for the selected landslide event to realize re-inversion.

[0116] 3.1 Spatial display of intermediate inversion results

[0117] The page displays the change of the value of each space grid point of the page over time in real time in the form of a heat map. Select the timing time and refresh frequency to set the time to a specific point. The refresh frequency is used to obtain the heat map data from the background. Select different components to display the intermediate inversion results from different dimensions.

[0118] The visualization system obtains the calculated inversion result data through Kafka, and saves the inversion result data in the form of different components (Z, E, N) to MySQL for persistent storage after processing the calculation time and frequency. At the same time, the grid coordinate data has been stored in the cache. The grid ID in the inversion result and the grid coordinates in the cache are combined to form the heat map coordinates for display. The grid list and grid coordinates are stored in the Redis cache respectively, and the grids and grid coordinates are added at the JVM level as secondary cache according to the page usage format. When the data calculated by the real-time processing subsystem is transmitted through the message and data system through Kafka, the grid coordinates and calculation results are organized into the specified format in combination with the JVM secondary cache and returned to the front-end page for real-time display.

[0119] The page has three sets of built-in color scales that can be switched. Select a different built-in color scale and the front end will render the point color according to the intermediate inversion result value corresponding to the grid point.

[0120] Set the maximum value of the intermediate inversion result: After modifying the maximum value from the page, click the OK button to re-query the data and the front end will automatically reset the interval and corresponding color of the intermediate inversion result legend according to the set maximum value. At the same time, if the inversion result exceeds the maximum value, it will default to the maximum value.

[0121] Set the minimum value of the intermediate inversion result: After modifying the minimum value from the page, click the OK button to re-query the data and the front end will automatically reset the interval and corresponding color of the intermediate inversion result legend according to the set minimum value. At the same time, if the inversion result is less than the minimum value, it will default to the minimum value.

[0122] Set point size: After modifying the minimum value on the page, click the OK button to requery the data and the front end will re-render the size of the grid points on the map according to the set point size. When the point is too small, if the distance between the grid points is too large, the points cannot be fused into a heat map, and the intermediate inversion information of the grid points cannot be displayed intuitively. When the point is too large, if the distance between the grid points is too small, the grid points will be directly fused together and cannot be effectively displayed.

[0123] 3.2 Time series display of intermediate inversion results

[0124] For a given grid point (one or more), the intermediate inversion results are obtained from the message exchange center, and the evolution characteristics of the inversion results over time are displayed in real time in the form of waveforms. For multiple grid points, it has the function of sorting according to different attributes; it has the function of scaling amplitude and time.

[0125] Please refer to Figure 2 , select the grid through the fuzzy query box, and the grid list is loaded to the browser on the front end of the page. The fuzzy query front end page is directly implemented without interacting with the back end every time to save system resources and network traffic. The page displays the data in the given time interval of one or more grids in the form of a line graph in the form of different components (Z, E, N). The grid can select multiple components at the same time, and these components are directly displayed in multiple line graphs. The page can also cancel the selected grid by clicking the selected grid multi-select group. The backend visualization system obtains the calculated inversion result data through Kafka, and saves the inversion result data in the form of different components (Z, E, N) to MySQL for persistent storage after processing the calculation time and calculation frequency. If the data in the selected time interval is not in the same calculation cycle, it will automatically go to the results of the next calculation cycle for query. If a grid does not have a calculation result at a certain time point, the data defaults to 0. The backend caches the grid list in Redis and the JVM distribution of the visualization system. The Redis cache is under great pressure due to the use of multiple systems. The visualization system caches the grid list twice, which can respond faster when the grid list is loaded on the page.

[0126] 3.3 Final inversion results display

[0127] Please refer to Figure 3 , lists the landslide events within a given time range. For the selected event, displays the inverted force time series, the distribution of stations used for inversion, the contour lines of inversion residuals, the matching status of observed waveforms and theoretical waveforms, etc.

[0128] The real-time processing system synchronizes the landslide event data and the final inversion results of the calculation to the visualization system through the message and data system, and the visualization system saves the data. When the landslide event is selected on the page and the OK button is clicked, the station coordinates used for inversion and the coordinates of the landslide event will be displayed in the form of a map on the right side, and the coordinate information will be displayed when the mouse is placed on it.

[0129] Below it, the inverted stress time series is displayed in the form of different components, the correlation coefficients between the observed waveform and the theoretical waveform, VR, and the inversion residual are given in the form of a list. At the bottom, the theoretical waveform and the observed waveform of the station used for inversion are displayed in the form of line graphs of different components.

[0130] It should be noted that the relevant parameters and corresponding operations are as follows:

[0131] Landslide list: displays landslide events. Click the OK button to display the corresponding final inversion results.

[0132] Manual inversion status: not calculated, no manual parameter adjustment has been performed; inverting, the system is inverting after manual parameter adjustment; inversion completed, there is inversion completed data after manual parameter adjustment; final result, the result of manual inversion is saved on the landslide event.

[0133] Adjust parameters: This is used to manually adjust parameters for re-inversion. For details, please refer to the human-computer interaction optimization inversion results below.

[0134] Save inversion results: Save the data of the human-computer interaction optimized inversion results to the landslide event to replace the original landslide event inversion results.

[0135] View inversion results: If the manual inversion result is inversion completed, it means there are inversion results with manual parameter adjustment. Click this button to view the latest inversion data.

[0136] Delete: Permanently delete the landslide event and the associated final inversion results.

[0137] 3.4 Human-computer interaction optimization inversion results

[0138] Please refer to Figure 4 For a given landslide event, by adjusting the inversion parameters, assigning different weights to seismic stations and components, and re-inverting, more reliable results can be obtained. By selecting different empirical formulas or manually giving formula coefficients, different landslide volume estimation results can be obtained.

[0139] On the final inversion result page, if the final inversion result is not satisfactory, you can adjust the parameters to re-invert. In the data inversion list, the adjusted parameters and the original event inversion results are pushed to the data real-time processing system through the Kafka message for re-inversion, and the new inversion results are displayed on the page. At this time, the list shows that the data inversion is over. If the final inversion result is satisfactory, click the Save button to replace the final inversion result data on the landslide event with the latest inversion result. If you are not satisfied, click Inversion again, the data will be re-inverted, and the landslide event can be deleted to discard the landslide event.

[0140] From the above description, we can know that the visualization subsystem obtains and displays the inversion results from the message and data subsystem, providing a visualization engine for the entire system. It can not only display the spatial distribution of the inversion results over time, but also display the evolution of the inversion results over time for a given grid point; for a selected landslide event, it can also display the comparison between the observed waveform and the theoretical waveform, the distribution of seismic stations, etc. At the same time, it also has some human-computer interaction functions: for a selected landslide event, the user can re-invert by adjusting the inversion parameters, screening the seismic records, etc., to obtain more reliable results.

[0141] The real-time surface dynamics monitoring system provided by the embodiment of the present invention can realize real-time inversion of surface dynamic parameters based on seismic data, identification of landslide events and parameter determination, etc., and visualize the results, which can provide data support for earthquake disaster prevention and control, earthquake prediction and prevention, etc.

[0142] The above is only a specific embodiment of the present invention, but the protection scope of the present invention is not limited thereto. Any technician familiar with the technical field can easily think of various equivalent modifications or replacements within the technical scope disclosed by the present invention, and these modifications or replacements should be included in the protection scope of the present invention. Therefore, the protection scope of the present invention shall be based on the protection scope of the claims.

Claims

1. A real-time monitoring system for surface dynamics, characterized in that: include: The message and data subsystem is responsible for the management of the entire process of message and data from access to archiving, and provides a message exchange and data acquisition platform for other subsystems; A real-time data processing subsystem, used to obtain seismic data from the message and data subsystem, perform real-time inversion of surface dynamic parameters, landslide event identification and parameter determination based on the seismic data, and return the results to the message and data subsystem; the results include real-time inversion results of surface dynamics and landslide event information; The visual display system is used to obtain messages and data from the message and data subsystem to achieve visual display.

2. The real-time monitoring system for surface dynamics according to claim 1, characterized in that: The message and data subsystem includes: A message exchange center is used to manage various time-series messages, including seismic data, real-time inversion results of surface dynamics returned by the real-time data processing subsystem, and landslide event information; Data cache center, used to provide seismic data within the time validity range; The first functional module is used to obtain seismic data from a seismic real-time data server or a seismic data file and store it in the message exchange center; The second functional module is used to obtain seismic data from the message exchange center in real time and store it in a data cache center; The third functional module is used to obtain the real-time inversion results of surface dynamics and landslide event information from the message exchange center and archive them.

3. The real-time monitoring system for surface dynamics according to claim 1, characterized in that: The messages and data acquired by the message and data subsystem include access data, and the access data includes continuous waveform data, offline seismic data and historical data; The processing of the continuous waveform data is as follows: Continuous waveform data is obtained based on two different protocols, NetSeis / IP and SeedLink; the continuous waveform data is accessed in the form of original MSEED data packets; Extracting network station information from the original MSEED data packet; Using the network station information as the key, the original MSEED data packets are directly stored in different partitions of the message exchange center to achieve distributed storage of data and facilitate parallel computing of subsequent processing modules; The processing method of offline seismic data is as follows: The MSEED data and SAC data are connected to the message and data subsystem in the form of files, which facilitates the real-time data processing subsystem to calculate surface dynamic parameters, identify landslide events and estimate landslide parameters.

4. The real-time monitoring system for surface dynamics according to claim 3, characterized in that: The access data also includes intermediate calculation results, which are processed as follows: The intermediate calculation results are divided into spatial grid points; the inversion results at each of the spatial grid points are event continuous data; The inversion results are stored in different partitions of the message exchange center using the grid point identifier as a key, and the storage format is JSON.

5. The real-time monitoring system for surface dynamics according to claim 4, characterized in that: The messages and data acquired by the message and data subsystem include the identified landslide event information, which is processed as follows: The identified landslide event information is stored in the data cache center after parameter optimization; each landslide event information at least includes event ID, longitude and latitude, estimated volume, inverted force time series, stations and components used in inversion, data frequency bands used in inversion, signal-to-noise ratio of each component observation data, and correlation coefficient between observation data and theoretical data.

6. The real-time monitoring system for surface dynamics according to claim 1, characterized in that: The real-time data processing subsystem comprises: Initialization module, used to initialize the system, including parameter file reading, spatial grid point division and Green's function calculation; A preprocessing module, used to obtain seismic data from the message and data subsystem in real time, perform preprocessing and calculate the signal-to-noise ratio; the preprocessing includes sliding averaging, bandpass filtering and removing instrument response; The surface dynamic parameter inversion module is used to traverse all spatial grid points, use the recorded data of the seismic stations associated with the grid points and the corresponding Green's function, and use the least squares method to complete the inversion of surface dynamic parameters in the time domain; The landslide event identification module is used to obtain the intermediate results of the inversion from the message and data subsystem, determine whether it exceeds a given threshold, and at the same time, determine whether the event exceeding the threshold is a landslide through the mutual correlation coefficient between the observed data and the theoretical data.

7. The real-time monitoring system for surface dynamics according to claim 6, characterized in that: The initialization module is specifically used for: Take the parameter file and divide the space grid points by given parameters; Bind the stations within a given distance range around each grid point to the grid point; Calculate the Green's function for each grid point; Initialize the instrument response at each station by given parameters; The station information, grid point information and the binding information between them are stored in the data cache center of the information and data subsystem.

8. The real-time monitoring system for surface dynamics according to claim 6, characterized in that: The real-time data processing subsystem also includes a parameter optimization module, which is used to: For the landslide events that are preliminarily identified, the spatial grid is refined, the Green's function at each grid point is calculated, and the inversion of dynamic parameters is completed; Traverse the inversion results and select the optimal position; After obtaining the results, the empirical formula is used to obtain the important parameters of landslide volume.

9. The real-time monitoring system for surface dynamics according to claim 1, characterized in that: The visual display subsystem includes: Visualization display module, used to obtain and display inversion results from the message and data subsystem, and provide a visualization engine for the entire system; The human-computer interaction module is used to receive user operations for the selected landslide event to realize re-inversion.

10. The real-time monitoring system for surface dynamics according to claim 9, characterized in that: The visual display module is specifically used for: Spatial display of intermediate inversion results: Obtain intermediate inversion results from the message exchange center and display the change of the value of each spatial grid point over time in real time; Time series display of intermediate inversion results: For one or more given grid points, the intermediate inversion results are obtained from the message exchange center, and the evolution characteristics of the inversion results over time are displayed in real time in the form of waveforms; For multiple grid points, it has the function of sorting by different attributes and scaling amplitude and time; Final inversion result display: landslide events within a given time range are given in a list format.

Images