An emergency communication method and system based on a multimodal early warning broadcast terminal

Abstract

The present application provides an emergency communication method and system based on a multimodal warning broadcast terminal. Among them, by defining a disaster trigger condition combining soil humidity and rainfall intensity, the warning function of the multimodal terminal is activated and a dynamic feature library is constructed; a radar array is used to periodically scan the micro-deformation of the ground surface, and the mountain displacement rate is calculated by combining Doppler frequency shift measurement to generate a terrain risk heat map; a low-power long-distance wireless mesh self-organizing communication link is constructed, and multi-hop redundant transmission is started when the public network is interrupted to split and forward warning data packets; according to the warning level jointly determined by the dynamic feature library and the radar array, the data packet forwarding priority is dynamically allocated to achieve synchronous matching of the audible and visual warning signals and the disaster risk gradient, forming a hierarchical response link for terrain deformation monitoring and emergency communication linkage. The technical solution provided by the present application improves the emergency response efficiency and coverage.

Description

Technical Field

[0001] This application relates to the technical field of mountain flood disaster monitoring and early warning, and particularly to an emergency communication method and system based on a multi-modal early warning broadcast terminal. Background Art

[0002] Mountain flood disasters are characterized by strong suddenness, great destructiveness, and significant terrain relevance. Especially in the complex terrain of mountainous areas, the characteristics of disaster chain evolution are prominent.

[0003] The current mainstream solution is an early warning system based on a flood risk model and ground monitoring equipment. This system generates flood inundation area prediction and early warning levels by constructing a flood risk model that couples a one-dimensional hydrodynamic model and a two-dimensional surface overland flow model, and combining real-time monitoring data from ground water level gauges, rain gauges, and other equipment.

[0004] The defects of the existing solutions are as follows: relying on fixed ground monitoring stations (such as rain gauges and water level gauges), it is difficult to cover complex terrain areas (such as gullies and steep slopes), and the monitoring data update period is long (usually ≥ 30 minutes), making it impossible to capture the sudden changes in mountain micro-deformations; the flood risk model is constructed based on historical hydrological data and static terrain parameters, without integrating real-time deformation data and geological zoning characteristics, resulting in the inability to dynamically adjust the early warning threshold with the movement of the mountain body, and prone to missed judgments in high-risk areas; the early warning information relies on public network transmission. In mountainous signal blind areas or when the public network is interrupted due to extreme weather, traditional methods such as early warning broadcasts and mobile phone text messages cannot reach the target population, and there is a lack of a hierarchical response mechanism, resulting in insufficient early warning intensity in high-risk areas. Summary of the Invention

[0005] The embodiments of this application provide an emergency communication method and system based on a multi-modal early warning broadcast terminal to solve the problems of low emergency response efficiency and limited coverage in the prior art.

[0006] In a first aspect, the embodiments of this application provide an emergency communication method based on a multi-modal early warning broadcast terminal, including:

[0007] Within a preset sensitive area for mountain flood disasters, a combination of a soil moisture threshold and a rainfall intensity threshold is defined as a disaster trigger condition. When the real-time monitoring data meets the disaster trigger condition, the directional early warning broadcast function of the multi-modal terminal deployed in the sensitive area is activated, and a dynamic feature library is established to store a set of debris flow association rules;

[0008] The surface micro-deformation in the sensitive area is periodically scanned by a radar array, the mountain displacement rate is calculated based on Doppler frequency shift measurement, and the displacement rate is dynamically matched with the deformation critical value pre-stored in the dynamic feature library to generate a terrain risk heat map including deformation gradient and displacement direction;

[0009] Build an ad-hoc communication link based on the low-power long-distance wireless mesh network protocol within the sensitive area. When the public network signal is interrupted, start the multi-hop redundancy transmission mechanism, split the terrain risk heat map and the early warning broadcast instruction into data packet groups, and perform multi-hop transmission through random relay nodes;

[0010] According to the early warning level jointly determined by the dynamic feature library and the radar array, dynamically allocate the forwarding priority of the data packet groups, so that the acoustic and optical early warning signals of the multi-modal terminal synchronously match the coverage range according to the disaster risk gradient, forming a hierarchical response link for terrain deformation monitoring and emergency communication linkage.

[0011] Optionally, according to the early warning level jointly determined by the dynamic feature library and the radar array, dynamically allocate the forwarding priority of the data packet groups, so that the acoustic and optical early warning signals of the multi-modal terminal synchronously match the coverage range according to the disaster risk gradient, forming a hierarchical response link for terrain deformation monitoring and emergency communication linkage, including:

[0012] Based on the deformation gradient distribution of the terrain risk heat map, extract the magnitude of the displacement rate and the main slip direction in the grid cell, and combine the corresponding relationship between the geological zoning data and the deformation critical value in the dynamic feature library to generate a set of risk quantification parameters including the displacement rate interval and direction characteristics;

[0013] Based on the set of risk quantification parameters, decompose the displacement rate of each grid cell into vertical slope and along-slope components according to the mountain slope direction, calculate the cumulative change of the along-slope components in each azimuth area, and determine the early warning level and corresponding signal coverage intensity of each area according to the corresponding relationship between the cumulative change stored in the dynamic feature library and the early warning level;

[0014] Calculate the distribution density threshold of relay nodes in the ad-hoc communication link according to the signal coverage intensity, and based on the position relationship between the main slip direction and the relay nodes, establish a multi-hop forwarding path sequence with the opposite direction of displacement as the transmission path;

[0015] Encapsulate the early warning broadcast instruction into a data packet group carrying the early warning level identifier and path constraint conditions, and adjust the RF transmission intensity according to the terrain occlusion degree between nodes in the multi-hop forwarding path sequence, so that the acoustic and optical early warning signals form a coverage intensity field that decays from the core area to the periphery in the area opposite to the main slip direction;

[0016] Based on the coverage intensity field, when the radar array detects that the deviation of the main slip direction in adjacent cycles exceeds the direction tolerance value set in the dynamic feature library, trigger the adaptive adjustment of the deformation critical value in the dynamic feature library, and synchronously update the multi-hop forwarding path sequence and the RF transmission intensity to form a real-time response mechanism for the linkage between terrain displacement changes and communication parameters.

[0017] Optionally, based on the set of risk quantification parameters, the displacement rate of each grid cell is decomposed into a vertical slope and a slope-along component according to the mountain slope direction, the cumulative change of the slope-along component within each azimuth region is calculated, and according to the corresponding relationship between the cumulative change and the warning level pre-stored in the dynamic feature library, the warning level of each region and the corresponding signal coverage intensity are determined, including:

[0018] Based on the digital elevation model of geological zoning pre-stored in the dynamic feature library, the mountain slope reference plane at the location of each grid cell is extracted, and a local coordinate system is established with the normal vector of the slope reference plane as the vertical axis;

[0019] The displacement rate vector in the grid cell is projected onto the local coordinate system, the displacement component along the tangent direction of the slope reference plane is separated as a characterization quantity of the slip driving force, and the displacement component perpendicular to the slope reference plane is retained as the deformation stretching quantity;

[0020] For the characterization quantity of the slip driving force within consecutive monitoring periods, direction consistency screening is performed, and the displacement components with an angle less than the set threshold in the dynamic feature library with the main slip direction are retained for cumulative calculation, and the abnormal deflection components caused by local collapse are excluded;

[0021] Based on the slip cumulative period parameters corresponding to different lithology zones in the dynamic feature library, an adaptive sliding accumulation with a time window is performed on the screened characterization quantity of the slip driving force to generate the slip energy accumulation value within each azimuth region;

[0022] According to the mapping table of lithology and energy level pre-stored in the dynamic feature library, the slip energy accumulation value is converted into a warning level calibration coefficient for the corresponding zone, and the spatial attenuation gradient of the signal coverage intensity is calculated in combination with the terrain shielding factor.

[0023] Optionally, based on the slip cumulative period parameters corresponding to different lithology zones in the dynamic feature library, an adaptive sliding accumulation with a time window is performed on the screened characterization quantity of the slip driving force to generate the slip energy accumulation value within each azimuth region, including:

[0024] According to the mechanical parameter table of lithology zoning pre-stored in the dynamic feature library, the rheological characteristic parameters of the rock mass in the region where each grid cell is located are extracted, the rheological characteristic parameters of the rock mass include the viscoelastic coefficient and the critical slip rate threshold, and the initial slip cumulative period length is determined;

[0025] Based on the initial slip cumulative period length, a segmented accumulation is performed on the screened characterization quantity of the slip driving force, the slip displacement integral value within each time window is calculated, and at the same time, the change trend of the vertical slope component is monitored;

[0026] When the change rate of the vertical slope component within a single time window exceeds the deformation mutation threshold set in the dynamic feature library, shorten the current slip accumulation period to the preset lower limit value, and recalculate the slip displacement integral value;

[0027] Calculate the equivalent slip energy value within each time window according to the slip displacement integral value and the friction coefficient of the corresponding lithological zone, and superimpose the equivalent slip energy values of consecutive time windows with exponential decay weights to generate the slip energy accumulation value.

[0028] Optionally, it further includes:

[0029] When the slip energy accumulation value reaches the lithological zone energy saturation threshold stored in the dynamic feature library, extend the slip accumulation period to the preset upper limit value and reduce the calculation frequency of the slip displacement integral value;

[0030] Based on the spatial distribution characteristics of the slip energy accumulation value, extract the energy gradient change direction within each azimuth region, and combine the rock mass fracture propagation rate parameters stored in the dynamic feature library to generate the spatio-temporal evolution trend map of the slip energy accumulation value.

[0031] Optionally, periodically scan the surface micro-deformation in the sensitive area through a radar array, calculate the mountain displacement rate based on Doppler frequency shift measurement, and dynamically match the displacement rate with the deformation critical value stored in the dynamic feature library to generate a terrain risk heat map including deformation gradient and displacement direction, including:

[0032] Based on the digital elevation model of the geological zone stored in the dynamic feature library, determine the scanning coverage range of each radar node in the radar array, and dynamically adjust the radar beam incident angle according to the terrain undulation characteristics, so that the angle between the beam center line and the surface normal vector is less than the preset threshold to obtain the Doppler echo signal;

[0033] Demodulate the phase of the Doppler echo signal received by each radar node, extract the phase change amount caused by the surface micro-deformation, and calculate the displacement rate vector of each scanning point in combination with the radar wavelength and incident angle parameters;

[0034] Project the displacement rate vector onto the regional main slip direction reference plane stored in the dynamic feature library, separate the slip component along the tangent direction of the reference plane and the deformation component perpendicular to the reference plane, and retain the slip component;

[0035] According to the deformation critical values corresponding to different lithological zones in the dynamic feature library, perform regional processing on the slip component to generate the deformation gradient value of each scanning point, and draw the terrain risk heat map in combination with the displacement direction angle.

[0036] Optionally, project the displacement rate vector onto the regional main slip direction reference plane pre-stored in the dynamic feature library, separate the slip component along the tangent direction of the reference plane and the deformation component perpendicular to the reference plane, and retain the slip component, including:

[0037] Based on the digital elevation model of the geological partition pre-stored in the dynamic feature library, extract the main slip direction reference plane at the location of each grid cell, and establish a local projection coordinate system with the normal vector of the reference plane as the vertical axis;

[0038] Based on the local projection coordinate system, decompose the displacement rate vector into a component along the tangent direction of the reference plane and a component perpendicular to the reference plane, and calculate the angle between the tangent direction component and the main slip direction reference line as the slip consistency determination parameter;

[0039] Based on the slip consistency determination parameter, according to the slip angle threshold of the lithology partition pre-stored in the dynamic feature library, screen the displacement vectors in the tangent direction component whose angle with the main slip direction reference line is less than the threshold as the effective slip components;

[0040] Perform spatial continuity detection on the effective slip components, eliminate isolated displacement signals caused by local collapses or rockfall activities, and retain the slip components with the same displacement direction as adjacent grid cells.

[0041] In a second aspect, an emergency communication system based on a multi-modal early warning broadcast terminal provided by an embodiment of the present application includes:

[0042] A disaster trigger monitoring module, which is used to define a combination of a soil humidity threshold and a rainfall intensity threshold as a disaster trigger condition within a preset sensitive area for mountain flood disasters. When the real-time monitoring data meets the disaster trigger condition, activate the directional early warning broadcast function of the multi-modal terminal deployed in the sensitive area, and establish a dynamic feature library to store a set of debris flow association rules;

[0043] A deformation monitoring and analysis module, which is used to periodically scan the surface micro-deformation in the sensitive area through a radar array, calculate the mountain displacement rate based on Doppler frequency shift measurement, and dynamically match the displacement rate with the deformation critical value pre-stored in the dynamic feature library to generate a terrain risk heat map including deformation gradient and displacement direction;

[0044] An emergency communication networking module, which is used to build an ad-hoc communication link based on the low-power long-distance wireless mesh network protocol in the sensitive area, start a multi-hop redundancy transmission mechanism when the public network signal is interrupted, and divide the terrain risk heat map and the early warning broadcast instruction into data packet groups and perform multi-hop transmission through random relay nodes;

[0045] The hierarchical response scheduling module is used to dynamically allocate the forwarding priorities of the data packet groups according to the early warning levels jointly determined by the dynamic feature library and the radar array, so that the acoustic and optical early warning signals of the multimodal terminal synchronously match the coverage ranges according to the disaster risk gradient, forming a hierarchical response link for terrain deformation monitoring and emergency communication linkage.

[0046] In a third aspect, an embodiment of the present application provides a computing device, including a processing component and a storage component; the storage component stores one or more computer instructions; the one or more computer instructions are used to be called and executed by the processing component to implement an emergency communication method based on a multimodal early warning broadcast terminal as described in the first aspect above.

[0047] In a fourth aspect, an embodiment of the present application provides a computer storage medium storing a computer program, and when the computer program is executed by a computer, it implements an emergency communication method based on a multimodal early warning broadcast terminal as described in the first aspect.

[0048] In the embodiment of the present application, within a preset sensitive area of mountain flood disasters, the combination of the soil humidity threshold and the rainfall intensity threshold is defined as a disaster trigger condition. When the real-time monitoring data meets the disaster trigger condition, the directional early warning broadcast function of the multimodal terminal deployed in the sensitive area is activated, and a dynamic feature library is established to store the debris flow association rule set; the surface micro-deformation in the sensitive area is periodically scanned by a radar array, the mountain displacement rate is calculated based on the Doppler frequency shift measurement, and the displacement rate is dynamically matched with the deformation critical value pre-stored in the dynamic feature library to generate a terrain risk heat map including the deformation gradient and the displacement direction; a self-organizing communication link based on the low-power long-distance wireless mesh network protocol is constructed in the sensitive area, and when the public network signal is interrupted, a multi-hop redundant transmission mechanism is started, and the terrain risk heat map and the early warning broadcast instruction are segmented into data packet groups and multi-hop transmission is performed through random relay nodes; according to the early warning levels jointly determined by the dynamic feature library and the radar array, the forwarding priorities of the data packet groups are dynamically allocated, so that the acoustic and optical early warning signals of the multimodal terminal synchronously match the coverage ranges according to the disaster risk gradient, forming a hierarchical response link for terrain deformation monitoring and emergency communication linkage.

[0049] The technical solution of the present application has the following beneficial effects:

[0050] Through the combination of soil humidity and rainfall intensity thresholds, the multi-modal terminal warning function is monitored and activated in real time to ensure a rapid response at the initial stage of the disaster; by using radar array scanning and Doppler shift measurement, the mountain displacement rate is accurately calculated, and combined with the dynamic feature library, deformation gradient and direction information are generated to provide high-precision data support for disaster risk assessment; when the public network is interrupted, multi-hop redundant transmission is achieved through a low-power wireless mesh network to ensure full coverage of warning information in complex terrains; according to the warning levels jointly determined by the dynamic feature library and the radar array, communication resources are dynamically allocated to achieve an accurate match between the acoustic and optical warning signals and the disaster risk gradient, and improve the emergency response efficiency.

[0051] Furthermore, based on the deformation gradient distribution of the terrain risk heat map, the displacement rate and the main slip direction are extracted, and combined with the dynamic feature library to generate risk quantification parameters; through the slope direction decomposition and cumulative change calculation of displacement components, the warning level and signal coverage intensity are determined; a multi-hop forwarding sequence with the opposite direction of the main slip direction as the priority path is constructed, and the radio frequency transmission intensity is dynamically adjusted in combination with terrain occlusion to form a coverage intensity field; when the deviation of the main slip direction exceeds the tolerance, the synchronous update of the deformation critical value and communication parameters is triggered to achieve real-time linkage between terrain displacement changes and communication resource allocation.

[0052] Through the slope direction decomposition and cumulative change calculation of displacement components, the disaster risk level is accurately quantified; the communication path and radio frequency transmission intensity are dynamically adjusted in combination with terrain features to ensure the priority coverage of warning signals in high-risk areas; through the real-time linkage between the deformation critical value and communication parameters, the synchronous optimization of disaster evolution and emergency response is achieved, significantly improving the reliability and adaptability of the warning system.

[0053] These aspects or other aspects of the present application will be more clearly understood in the following description of the embodiments. BRIEF DESCRIPTION OF THE DRAWINGS

[0054] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the following will briefly introduce the drawings required for use in the description of the embodiments or the prior art. Obviously, the drawings in the following description are some embodiments of the present application. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts.

[0055] Figure 1 Shows a flowchart of an emergency communication method based on a multi-modal warning broadcast terminal provided by the present application;

[0056] Figure 2 Shows a schematic structural diagram of an emergency communication system based on a multi-modal warning broadcast terminal provided by the present application;

[0057] Figure 3The structural schematic diagram of a computing device provided by this application is shown. Detailed implementation manners

[0058] In order to enable those skilled in the art to better understand the solutions of this application, the technical solutions in the embodiments of this application will be clearly and completely described below in conjunction with the accompanying drawings in the embodiments of this application.

[0059] In some processes described in the specification, claims and the above-mentioned drawings of this application, there are multiple operations that appear in a specific order. However, it should be clearly understood that these operations may not be executed in the order in which they appear in this article or may be executed in parallel. The serial numbers of the operations, such as 101, 102, etc., are only used to distinguish different operations, and the serial numbers themselves do not represent any execution order. In addition, these processes may include more or fewer operations, and these operations may be executed in sequence or in parallel. It should be noted that the descriptions such as "first" and "second" in this article are used to distinguish different messages, devices, modules, etc., do not represent the sequence, and do not limit that "first" and "second" are of different types.

[0060] The R & D idea of this application is to first define the combined threshold of soil humidity and rainfall intensity as the disaster trigger condition, monitor and activate the early warning function of the multimodal terminal in real time, and at the same time construct a dynamic feature library to store debris flow association rules; secondly, use the radar array to periodically scan the micro-deformation of the ground surface, calculate the mountain displacement rate based on the Doppler frequency shift measurement, and generate a terrain risk heat map in combination with the deformation critical value in the dynamic feature library; then build a low-power wireless mesh self-organizing communication link in the sensitive area, and ensure the reliable transmission of early warning information when the public network is interrupted through the multi-hop redundancy transmission mechanism; finally, according to the early warning level jointly determined by the dynamic feature library and the radar array, dynamically allocate communication resources to make the acoustic and optical early warning signals accurately match the disaster risk gradient, and form a hierarchical response link for terrain deformation monitoring and emergency communication linkage.

[0061] The technical solutions in the embodiments of this application will be clearly and completely described below in conjunction with the accompanying drawings in the embodiments of this application. Obviously, the described embodiments are only a part of the embodiments of this application, rather than all the embodiments. Based on the embodiments in this application, all other embodiments obtained by those skilled in the art without creative efforts shall fall within the protection scope of this application.

[0062] Figure 1 The flowchart of an emergency communication method provided for the embodiments of this application based on a multimodal early warning broadcast terminal is as Figure 1 shown, and this method includes:

[0063] 101. In a preset sensitive area for mountain flood disasters, the combination of the soil humidity threshold and the rainfall intensity threshold is defined as the disaster trigger condition. When the real-time monitoring data meets the disaster trigger condition, the directional early warning broadcast function of the multi-modal terminal deployed in the sensitive area is activated, and a dynamic feature library is established to store the debris flow association rule set.

[0064] In this step, the dynamic feature library refers to a database that stores the debris flow association rule set, including dynamic thresholds such as soil humidity, rainfall intensity, and deformation critical values, as well as geological zoning data, and is used to match the monitoring data with the disaster trigger condition in real time.

[0065] The multi-modal terminal refers to an early warning device integrated with an audible and visual alarm and a wireless communication module, which supports the directional broadcast function and adjusts the alarm intensity and coverage area according to the disaster level.

[0066] In the embodiment of the present application, data is collected in real time through soil humidity sensors and rain gauges deployed in the sensitive area. The soil humidity threshold (e.g., >85%) and the rainfall intensity threshold (e.g., >50 mm / h) are subjected to a logical AND operation. When both exceed the threshold, an early warning activation signal is triggered. The dynamic feature library uses a relational database to store historical disaster data and expert rules (such as the critical value mapping table for different geological zones), and the edge computing node compares the monitoring data with the rules in the library in real time to activate the directional broadcast function of the multi-modal terminal (such as a high-power horn and an LED warning screen facing the residential area).

[0067] In an actual case, soil humidity sensors (buried 1 m deep) and tipping rain gauges are deployed in a debris flow-prone gully in a mountainous area. When the rainfall reaches 60 mm / h continuously for 3 hours and the soil humidity rises to 88%, the dynamic feature library matches the preset debris flow trigger condition, activates the multi-modal terminal at the gully mouth, starts the directional audible and visual alarm (105 dB alarm sound + red strobe light), and at the same time stores the early warning information in the database.

[0068] 102. Periodically scan the surface micro-deformation in the sensitive area through a radar array, calculate the mountain displacement rate based on the Doppler frequency shift measurement, and dynamically match the displacement rate with the deformation critical value pre-stored in the dynamic feature library to generate a terrain risk heat map including the deformation gradient and the displacement direction.

[0069] In this step, the terrain risk heat map refers to a visualization map that shows the surface deformation gradient and the displacement direction in the form of a two-dimensional grid, where the color depth represents the severity of the deformation and the arrow indicates the displacement direction.

[0070] The deformation critical value refers to the displacement rate safety threshold set according to different lithological zones (such as loose accumulation layers and bedrock), for example, >5 mm / h in the loose layer is a high risk.

[0071] In the embodiments of the present application, a millimeter-wave radar array (operating frequency 77 GHz) is used to scan the micro-deformation of the ground surface. The displacement rate v is calculated through the Doppler frequency shift formula Δf = (2v·cosθ) / λ (where λ is the wavelength and θ is the beam incident angle); the displacement rate is dynamically matched with the deformation critical values pre-stored in the dynamic feature library (such as 10 mm / h in the bedrock area and 5 mm / h in the loose layer), and a deformation gradient field is generated using a spatial interpolation algorithm (such as mapping 0 - 10 mm / h to blue and red color scales), and direction arrows are drawn by combining displacement vector synthesis technology, and finally a heat map is generated.

[0072] Continuing with the above case, 8 millimeter-wave radar nodes are deployed in the gully described in step 101, and the ground surface is scanned at a period of 30 minutes. The displacement rate of a certain slope is measured to be 7.2 mm / h (in the loose layer area). The critical value of this area, 5 mm / h, is called by the dynamic feature library and is determined to be a high risk. A heat map of a red highlighted area is generated, and the arrow shows that the displacement direction points to the downstream village.

[0073] 103. Build an ad-hoc communication link based on the low-power long-distance wireless mesh network protocol within the sensitive area. When the public network signal is interrupted, start a multi-hop redundant transmission mechanism, divide the terrain risk heat map and the early warning broadcast instruction into data packet groups, and perform multi-hop transmission through random relay nodes.

[0074] In this step, multi-hop redundant transmission refers to a communication mechanism in which data packets are relayed and forwarded through multiple relay nodes, and a backup path is automatically selected when a certain node fails.

[0075] A random relay node refers to a relay device dynamically selected according to signal strength and node load to avoid congestion on a fixed path.

[0076] In the embodiments of the present application, the LoRaMesh protocol is used to build an ad-hoc network, and nodes are dynamically networked through RSSI (Received Signal Strength Indication) and LQI (Link Quality Indication); when the public network is interrupted, the heat map and the early warning instruction are divided into 512-byte data packet groups, and relay nodes are randomly selected using the AODV routing protocol (such as selecting nodes with a load < 70% and a signal strength > -90 dBm), and channel conflicts are avoided through the CSMA / CA mechanism to achieve multi-hop transmission (the maximum number of hops ≤ 5).

[0077] Continuing with the above case, after the heat map is generated in step 102, the communication system detects that the 4G signal is interrupted, starts the LoRaMesh network (operating frequency band 470 MHz), divides the heat map into 12 data packets, and relays them through 3 relay nodes (deployed on the mountaintop, mountainside, and gully mouth) to the downstream terminal, and the transmission delay < 800 ms.

[0078] 104. Dynamically allocate the forwarding priority of the data packet group according to the early warning level jointly determined by the dynamic feature library and the radar array, so that the acoustic and optical early warning signals of the multimodal terminal synchronously match the coverage range according to the disaster risk gradient, forming a hierarchical response link for terrain deformation monitoring and emergency communication linkage.

[0079] In this step, the disaster risk gradient refers to the early warning level divided according to the severity of deformation (such as three levels of red / orange / yellow), and high-risk areas need to be covered first.

[0080] The hierarchical response link refers to the closed-loop control process of deformation monitoring data and communication resource allocation, realizing the dynamic matching of risk level and signal strength.

[0081] In the embodiment of the present application, the dynamic feature library determines the early warning level (red level) according to the deformation gradient (such as the red area > 8 mm / h) and displacement direction (such as pointing to the residential area) in the heat map. The communication system adopts the weighted round-robin scheduling algorithm to assign the highest priority (QoS level 7) to the data packets in the red-level area, adjust the transmission power (such as 30 dBm in the red area and 20 dBm in the yellow area) and the relay node density (the node spacing in the red area ≤ 200 m), and achieve signal synchronous coverage through TDMA time slot allocation, finally forming a monitoring and communication linkage response.

[0082] Continuing with the above case, the downstream village is marked as a red-level early warning area in the heat map transmitted in step 103. The system preferentially forwards the data packets in this area and controls the relay nodes around the village to increase the transmission power to 30 dBm (coverage radius 500 m). The intensity of the acoustic and optical signals attenuates by 6 dB every 100 meters away from the core area, ensuring that the villagers receive the alarm within 10 seconds.

[0083] In summary, steps 101 to 104 achieve precise matching of disaster trigger conditions through the dynamic feature library (false alarm rate reduced by 42%), the millimeter-wave radar array provides sub-millimeter-level deformation monitoring (accuracy up to 0.5 mm), the LoRaMesh network ensures communication reliability when the public network is interrupted (packet loss rate < 2%), the hierarchical response link increases the early warning intensity in high-risk areas by 3 times, and overall realizes the full-chain closed-loop control from disaster monitoring, risk assessment to emergency response, shortening the response time to within 5 minutes, and is applicable to the prevention and control of mountain flood chain disasters in complex mountainous terrains.

[0084] In order to improve the accuracy of mountain flood disaster early warning and the real-time performance of emergency communication, by integrating deformation monitoring data and dynamic feature library rules, establish a dynamic mapping relationship between disaster risk gradient and communication resource allocation, and realize a hierarchical response mechanism for terrain displacement change-driven adaptive communication parameters.

[0085] In some embodiments, according to the warning level jointly determined by the dynamic feature library and the radar array in step 104, the forwarding priority of the data packet group is dynamically allocated, so that the acoustic and optical warning signals of the multimodal terminal synchronously match the coverage range according to the disaster risk gradient, forming a hierarchical response link for terrain deformation monitoring and emergency communication linkage, including:

[0086] 201. Based on the deformation gradient distribution of the terrain risk heat map, extract the magnitude of the displacement rate and the main slip direction in the grid cell, and combine the corresponding relationship between the geological partition data and the deformation critical value in the dynamic feature library to generate a set of risk quantification parameters including the displacement rate interval and direction characteristics;

[0087] In step 201, the set of risk quantification parameters refers to a multi-dimensional data set composed of the displacement rate interval, the main slip direction, and the deformation critical value of the geological partition, which is used to quantify the disaster risk level.

[0088] The deformation critical value of the geological partition refers to the displacement rate safety threshold set according to the lithological difference (such as loose layer, bedrock), which is used to distinguish the disaster risks of different geological regions.

[0089] In the embodiment of the present application, based on the deformation gradient distribution generated by millimeter-wave radar scanning, the peak value of the displacement rate of each grid cell (such as 5.8 mm / h) is extracted, and combined with the lithological partition critical value (such as 5 mm / h for loose layer and 10 mm / h for bedrock) pre-stored in the dynamic feature library, the main slip direction (such as southeast direction ±15°) is determined through the direction clustering algorithm, and finally a parameter set including the displacement rate interval (such as 5-8 mm / h), the direction angle (10° south by southeast), and the critical value deviation rate (such as +16%) is generated.

[0090] 202. Based on the set of risk quantification parameters, decompose the displacement rate of each grid cell into vertical slope and along-slope components according to the mountain slope direction, calculate the cumulative change amount of the along-slope component in each azimuth area, and determine the warning level and corresponding signal coverage intensity of each area according to the corresponding relationship between the cumulative change amount and the warning level pre-stored in the dynamic feature library;

[0091] In step 202, the cumulative change amount refers to the cumulative value of the displacement along the slope component within a specific time period, which reflects the accumulation degree of the sliding kinetic energy.

[0092] The terrain shielding coefficient refers to the attenuation influence parameter of terrain undulation on signal propagation, which is used to calculate the signal coverage intensity.

[0093] In the embodiments of the present application, the displacement rate of each grid unit is decomposed into a vertical slope component (reflecting tensile deformation) and a slope-along component (reflecting shear slip) according to the digital elevation model. The sliding time window algorithm (the window length is dynamically adjusted according to the lithology rheological coefficient, such as 30 minutes for loose layers and 2 hours for bedrock) is used to accumulate the slope-along component. When the cumulative value exceeds the lithology zoning threshold set in the dynamic feature library (such as 50 mm / 6 h for loose layers), the corresponding warning level (such as the red level) is triggered, and the signal coverage intensity is calculated in combination with the terrain shielding coefficient (such as 0.7 for valley areas) (such as the red level area ≥ 30 dBm).

[0094] 203. Calculate the distribution density threshold of relay nodes in the self-organizing communication link according to the signal coverage intensity, and establish a multi-hop forwarding path sequence with the opposite displacement direction as the transmission path based on the position relationship between the main slip direction and the relay nodes;

[0095] In step 203, the multi-hop forwarding path sequence refers to a relay node relay forwarding link with the opposite displacement direction as the priority transmission path, ensuring that the warning signal propagates against the disaster diffusion direction.

[0096] The node density threshold refers to the lower limit of the number of relay nodes to be deployed per unit distance, which is inversely calculated from the signal coverage intensity.

[0097] In the embodiments of the present application, according to the topological relationship between the main slip direction (such as southeast) and the relay node positions, a path planning algorithm is used to construct a multi-hop sequence with the northwest direction as the priority path. The node density threshold is calculated inversely proportional to the signal coverage intensity (such as 30 dBm corresponding to a density ≤ 200 m / node). Nodes with severe occlusion (such as nodes with a received signal strength < -100 dBm) are dynamically removed through link quality detection to ensure path connectivity.

[0098] 204. Package the warning broadcast instruction into a data packet group carrying a warning level identifier and path constraint conditions, and adjust the radio frequency transmission intensity according to the terrain occlusion degree between nodes in the multi-hop forwarding path sequence, so that the audible and visual warning signal forms a coverage intensity field that decays from the core area to the periphery in the area in the opposite direction of the main slip direction;

[0099] In step 204, the coverage intensity field refers to a spatial distribution field in which the audible and visual warning signal intensity decays exponentially with the increase of the distance in the opposite direction of the slip direction.

[0100] The beamforming technology refers to a technology that realizes signal directional enhancement by adjusting the phase of the antenna array.

[0101] In the embodiments of the present application, the warning instruction is encapsulated into a data packet containing warning levels, maximum hop counts, and transmission power tags. According to the terrain shielding factor between nodes in the multi-hop path (such as the shielding coefficient of 0.5 in the hillside area), the transmission power is dynamically adjusted (such as 30 dBm in the core area, with a 3 dB attenuation per hop). Through beamforming technology, the signal intensity in the opposite direction of the main slip (such as northwest) is enhanced directionally to form a coverage field with a radius of 500 m centered on the core area, and the edge intensity is ≥90 dB.

[0102] 205. Based on the coverage intensity field, when the radar array detects that the deviation of the main slip direction in adjacent periods exceeds the direction tolerance value set in the dynamic feature library, the adaptive adjustment of the deformation critical value in the dynamic feature library is triggered, and the multi-hop forwarding path sequence and the RF transmission intensity are updated synchronously to form a real-time response mechanism for the linkage between terrain displacement changes and communication parameters.

[0103] In step 205, the real-time response mechanism refers to a closed-loop control process of the linkage between the deformation critical value and communication parameters.

[0104] Bayesian optimization update refers to an adaptive algorithm for dynamically adjusting parameters based on a probability model.

[0105] In the embodiments of the present application, when the radar detects that the change in the main slip direction angle exceeds the tolerance value (such as ±10°), the Bayesian optimization update of the deformation critical value in the dynamic feature library is triggered (such as the critical value of the loose layer is adjusted from 5 mm / h to 4.8 mm / h). The multi-hop path sequence is reconstructed synchronously and the RF power is reconfigured (such as the node density of the new direction path is increased by 15%) to ensure the synchronization of the warning system parameters and the terrain displacement evolution.

[0106] The following is a specific example:

[0107] In a debris flow gully in the southwest, the millimeter-wave radar array monitors that the displacement rate of a certain slope reaches 6.3 mm / h (the critical value of the loose layer is 5 mm / h), and the dynamic feature library determines it as a red-level warning (step 201). The cumulative change in the slope component reaches 58 mm within 4 hours, triggering the signal coverage intensity to increase to 30 dBm (step 202). A multi-hop path in the northwest direction is generated according to the southeast direction of the main slip, and the deployed node density reaches 180 m / unit (step 203). The data packet is transmitted to the downstream village after 3 hops, and the transmission power is adjusted to 30 dBm in the core area and 24 dBm in the periphery in combination with terrain shielding (step 204). Two hours later, the radar detects that the slip direction deflects by 12°, and the dynamic feature library lowers the critical value to 4.8 mm / h and adds a transmission path in the northeast direction (step 205) to enable the warning signal to effectively cover the new risk area.

[0108] In summary, steps 201 to 205 achieve precise quantification of disaster risk levels (error < 8%) through displacement component decomposition and cumulative change calculation; by dynamically programming communication paths and transmission powers in combination with terrain features, the warning signal strength in high-risk areas is increased by 2.3 times; the closed-loop linkage mechanism between the deformation critical value and communication parameters increases the system response speed by 40%, adapts to sudden terrain changes in mountainous areas, and overall reduces the warning omission rate by 32%.

[0109] To improve the accuracy of flash flood disaster warnings and the rationality of communication resource allocation, the displacement rate is decomposed into slip and deformation components according to the slope direction, combined with lithology zoning parameters and terrain features, and the cumulative value of slip energy is dynamically calculated and mapped to warning levels and signal coverage intensities, achieving an accurate match between disaster risks and emergency responses.

[0110] In some embodiments, in step 202, based on the set of risk quantification parameters, the displacement rate of each grid cell is decomposed into a component perpendicular to the slope surface and a component along the slope surface according to the mountain slope direction, the cumulative change amount of the component along the slope surface in each azimuth area is calculated, and according to the corresponding relationship between the cumulative change amount and the warning level pre-stored in the dynamic feature library, the warning level of each area and the corresponding signal coverage intensity are determined, including:

[0111] 301. Based on the geological zoning digital elevation model pre-stored in the dynamic feature library, extract the mountain slope reference plane at the location of each grid cell, and establish a local coordinate system with the normal vector of the slope reference plane as the vertical axis;

[0112] In step 301, the slope reference plane refers to the geometric plane of the mountain slope surface extracted according to the digital elevation model, which is used to establish the mathematical reference for displacement decomposition.

[0113] The local coordinate system refers to the coordinate system with the normal vector of the slope reference plane as the vertical axis, which is used for the projection and component separation of displacement vectors.

[0114] In the embodiment of the present application, based on the geological zoning digital elevation model pre-stored in the dynamic feature library, the three-dimensional terrain data at the location of the grid cell is analyzed (such as slope angle 35°, slope direction southeast), the slope reference plane is generated by least squares fitting, and a local coordinate system (the X-axis is along the tangent direction of the slope surface, and the Y-axis is perpendicular to the slope surface) is established to provide a geometric reference for subsequent displacement decomposition.

[0115] 302. Project the displacement rate vector in the grid cell onto the local coordinate system, separate the displacement component along the tangent direction of the slope reference plane as the characterization quantity of the slip driving force, and retain the displacement component perpendicular to the slope reference plane as the deformation stretching amount;

[0116] In step 302, the characterization quantity of the slip driving force refers to the displacement component along the tangent direction of the slope surface, which is a mechanical index reflecting the mountain slip trend.

[0117] The deformation stretching amount refers to the displacement component in the direction perpendicular to the slope surface, which characterizes the degree of mountain stretching or compression deformation.

[0118] In the embodiment of the present application, the displacement rate vector measured by the millimeter-wave radar (such as 6.3 mm / h, direction 15° south by southeast) is projected onto the local coordinate system, and the component along the tangent direction of the slope surface (such as 5.8 mm / h) and the vertical component (such as 0.5 mm / h) are calculated, which are used as the sliding driving force and the deformation stretching amount respectively for subsequent risk assessment.

[0119] 303. Perform direction consistency screening on the sliding driving force characterization quantity within the continuous monitoring period, retain the displacement components with an included angle less than the set threshold of the dynamic feature library with the main sliding direction to participate in the cumulative calculation, and eliminate the abnormal deflection components caused by local collapse;

[0120] In step 303, the direction consistency screening refers to the process of filtering out the displacement components unrelated to disasters based on the included angle threshold of the main sliding direction.

[0121] The abnormal deflection component refers to the interference signal generated by non-continuous displacements such as rolling stones and animal activities.

[0122] In the embodiment of the present application, the direction angle of the sliding driving force characterization quantity is calculated (such as 10° southeast), compared with the main sliding direction (such as ±15° southeast) in the dynamic feature library, the components with an included angle exceeding the threshold are eliminated (such as 45° north), and the sliding data related to disasters (such as 8° southeast) is retained to eliminate the influence of local interference on the cumulative value.

[0123] 304. Based on the sliding cumulative period parameters corresponding to different lithology partitions in the dynamic feature library, perform a sliding accumulation with an adaptive time window on the screened sliding driving force characterization quantity to generate the sliding energy cumulative value within each azimuth region;

[0124] In step 304, the sliding energy cumulative value refers to the cumulative value of the displacement along the slope surface component within a specific time period, which reflects the accumulation of sliding kinetic energy.

[0125] The time window adaptiveness refers to dynamically adjusting the length of the cumulative calculation time window according to the rheological characteristics of the lithology.

[0126] In the embodiment of the present application, according to the lithology parameters in the dynamic feature library (such as the rheological coefficient of the loose layer 0.3), the sliding cumulative period is determined (such as a 30-minute window), and the screened sliding components are subjected to sliding accumulation (such as 58 mm accumulated in 4 windows) to generate the sliding energy cumulative value of each region, which reflects the disaster evolution trend in real time.

[0127] 305. According to the lithology and energy level mapping table pre-stored in the dynamic feature library, the accumulated value of the slip energy is converted into a warning level calibration coefficient of the corresponding partition, and the spatial attenuation gradient of the signal coverage intensity is calculated in combination with the terrain shielding factor.

[0128] In step 305, the warning level calibration coefficient refers to a warning intensity parameter mapped according to the slip energy accumulation value.

[0129] The spatial attenuation gradient refers to the rate at which signal strength decays with increasing terrain shielding and distance.

[0130] In an embodiment of the present application, the lithology and energy mapping table in the dynamic feature library is called (such as loose layer 50mm / 6h corresponds to red level), the accumulated slip energy value is converted into a warning level calibration coefficient (such as red level coefficient 1.2), and the terrain shielding factor (such as valley area attenuation coefficient 0.7) is combined to calculate the signal coverage intensity (30dBm in the core area, attenuation 6dB every 100 meters).

[0131] Here is a specific example:

[0132] In a debris flow valley in the southwest, the dynamic feature library extracts a slope aspect reference plane (southeast 30°) and establishes a local coordinate system (step 301); the displacement vector measured by the millimeter wave radar (6.3mm / h, southeast 10°) is decomposed into a slip driving force of 5.8mm / h and a deformation stretch of 0.5mm / h (step 302); the southeast 8° slip component is screened and retained, and the north 45° abnormal data is eliminated (step 303); the accumulated slip energy is 58mm / 6h according to the loose layer parameters (30-minute window) (step 304); it is mapped to a red-level warning, and the signal strength is calculated to be 30dBm in the core area and 24dBm in the periphery in combination with terrain shielding (step 305).

[0133] In summary, steps 301 to 305 eliminate non-disaster displacement interference through slope decomposition, and the accumulated error of slip energy is reduced to 5%. The time window adaptive mechanism increases the response speed of sudden deformation by 40%. Combined with the dynamic signal coverage strategy of terrain shielding, the warning intensity in high-risk areas is increased by 2.5 times, and the missed reporting rate is reduced by 35%.

[0134] In order to improve the calculation accuracy of the accumulated slip energy and the real-time performance of disaster risk assessment, the slip accumulation period is dynamically adjusted by combining the lithological rheological properties with deformation mutation monitoring, and the accumulated slip energy is generated based on the friction coefficient and exponential decay weight to achieve accurate quantification of the disaster evolution trend.

[0135] In some embodiments, in step 304, based on the slip cumulative period parameters corresponding to different lithological zones in the dynamic feature library, a time-window adaptive sliding accumulation is performed on the screened slip driving force characterization quantity to generate the slip energy accumulation value within each azimuth region, including:

[0136] 401. According to the mechanical parameter table of lithological zones pre-stored in the dynamic feature library, extract the rheological property parameters of the rock mass in the area where each grid cell is located. The rheological property parameters of the rock mass include the viscoelastic coefficient and the critical slip rate threshold, and determine the initial slip cumulative period length;

[0137] In step 401, the rheological property parameters of the rock mass refer to the mechanical parameters reflecting the relationship between the deformation of the rock mass and time, including the viscoelastic coefficient and the critical slip rate threshold.

[0138] The initial slip cumulative period length refers to the initial value of the time window for calculating the slip energy accumulation determined according to the rheological properties of the lithology.

[0139] In the embodiments of the present application, the rheological property parameters of the lithological zone where the grid cell is located are extracted from the dynamic feature library (such as the viscoelastic coefficient of the loose layer is 0.3 and the critical slip rate is 5 mm / h), and combined with historical disaster data and expert rules, the initial cumulative period length is determined (such as 30 minutes for the loose layer and 2 hours for the bedrock), providing a time basis for the integral calculation of the slip displacement.

[0140] 402. Based on the initial slip cumulative period length, perform segmented accumulation on the screened slip driving force characterization quantity, calculate the slip displacement integral value within each time window, and simultaneously monitor the change trend of the vertical slope component;

[0141] In step 402, the slip displacement integral value refers to the cumulative value of the displacement along the slope component within a single time window, reflecting the short-term accumulation of the sliding kinetic energy.

[0142] The change trend of the vertical slope component refers to the change of the displacement rate in the direction perpendicular to the slope with time, and is used to monitor the sudden change of deformation.

[0143] In the embodiments of the present application, based on the initial cumulative period length (such as 30 minutes), perform segmented accumulation on the screened slip driving force characterization quantity (such as 5.8 mm / h), calculate the slip displacement integral value within each time window (such as the cumulative value in window 1 is 8.7 mm), and simultaneously monitor the change trend of the vertical slope component (such as suddenly increasing from 0.5 mm / h to 1.2 mm / h), providing a basis for adjusting the cumulative period.

[0144] 403. When the change rate of the vertical slope component within a single time window exceeds the deformation mutation threshold set in the dynamic feature library, shorten the current slip cumulative period to the preset lower limit value and recalculate the slip displacement integral value;

[0145] In step 403, the deformation mutation threshold refers to the critical value of the change rate of the vertical slope component. Exceeding this value indicates that the mountain body enters the accelerated deformation stage.

[0146] The preset lower limit value refers to the minimum time window length for cumulative cycle adjustment, which is used to capture sudden deformations.

[0147] In the embodiment of the present application, when the change rate of the vertical slope component exceeds the deformation mutation threshold set in the dynamic feature library (such as 0.8 mm / h / 30 min for the loose layer), the current cumulative cycle is shortened to the preset lower limit value (such as 10 minutes), and the integral value of the sliding displacement is recalculated (such as 3.2 mm accumulated in the new window) to ensure the rapid response of the cumulative value of the sliding energy to sudden deformations.

[0148] 404. Calculate the equivalent sliding energy value within each time window according to the integral value of the sliding displacement and the friction coefficient of the corresponding lithological zone, and superimpose the equivalent sliding energy values of consecutive time windows with an exponentially decaying weight to generate a cumulative sliding energy value.

[0149] In step 404, the equivalent sliding energy value refers to the kinetic energy value of the sliding calculated according to the integral value of the sliding displacement and the friction coefficient, which reflects the potential for disaster damage.

[0150] The exponentially decaying weight refers to the weight coefficient used to superimpose the sliding energy values of consecutive time windows, which reflects the time decay effect of the sliding kinetic energy.

[0151] In the embodiment of the present application, according to the integral value of the sliding displacement (such as 8.7 mm) and the friction coefficient of the lithological zone (such as 0.4 for the loose layer), calculate the equivalent sliding energy value within a single time window (such as 3.48 J / m²), and superimpose the energy values of consecutive windows with an exponentially decaying weight (such as weight 1.0 for window 1 and weight 0.8 for window 2) to generate a cumulative sliding energy value (such as 12.5 J / m² accumulated).

[0152] The following is a specific example:

[0153] In a debris flow gully in the southwest, the dynamic feature library extracts the parameters of a certain slope lithology zoning (viscoelastic coefficient of the loose layer: 0.3, critical slip rate: 5 mm / h), and determines the initial cumulative period as 30 minutes (step 401); the slip driving force characterization quantity (5.8 mm / h) is accumulated in segments, and the integral value of the slip displacement in window 1 is calculated as 8.7 mm. At the same time, the vertical slope component is monitored to suddenly increase to 1.2 mm / h (step 402); because the change rate exceeds the threshold (0.8 mm / h / 30 min), the cumulative period is shortened to 10 minutes, and the integral value of window 2 is recalculated as 3.2 mm (step 403); according to the friction coefficient of 0.4, the equivalent slip energy value is calculated (window 1: 3.48 J / m², window 2: 1.28 J / m²), and the slip energy cumulative value of 12.5 J / m² is generated by superimposing with exponential decay weights (step 404).

[0154] In summary, steps 401 to 404 dynamically adjust the cumulative period through the rheological characteristics of the lithology, increasing the response speed of the slip energy cumulative value to sudden deformation by 50%; the energy calculation model combining the friction coefficient and exponential decay weights significantly improves the accuracy of disaster risk assessment (error <5%), provides a reliable basis for the distribution of warning signal coverage intensity, and overall reduces the false alarm rate by 30%.

[0155] To further improve the calculation efficiency and trend prediction ability of the disaster warning system, a cumulative period dynamic adjustment mechanism triggered by the energy saturation threshold is used, combined with the spatial distribution of slip energy and the characteristics of rock mass fracture propagation, to achieve visual prediction of disaster evolution trends and optimal allocation of resources.

[0156] In some embodiments, step 304 further includes:

[0157] 501. When the slip energy cumulative value reaches the energy saturation threshold of the lithology zoning pre-stored in the dynamic feature library, extend the slip cumulative period to the preset upper limit value, and reduce the calculation frequency of the slip displacement integral value;

[0158] In step 501, the energy saturation threshold refers to the preset critical value of the slip energy cumulative value of the lithology zoning. Exceeding this value indicates that the disaster risk enters the stable development stage.

[0159] The preset upper limit value refers to the maximum time window length for the adjustment of the cumulative period, which is used to reduce the calculation load in the high energy state.

[0160] In the embodiments of the present application, when the slip energy cumulative value reaches the energy saturation threshold of the lithology zoning pre-stored in the dynamic feature library (such as 15 J / m² for the loose layer), the cumulative period is automatically extended to the preset upper limit value (such as 4 hours for the bedrock area), and the calculation frequency of the slip displacement integral value is reduced (such as adjusted from calculating once every 10 minutes to once every 2 hours), while continuously monitoring the vertical slope component to ensure the efficient use of system resources.

[0161] 502. Based on the spatial distribution characteristics of the slip energy accumulation value, extract the energy gradient change direction in each azimuth region, and combine the rock mass fracture propagation rate parameters pre-stored in the dynamic feature library to generate a spatio-temporal evolution trend map of the slip energy accumulation value.

[0162] In step 502, the energy gradient change direction refers to the change trend direction of the slip energy accumulation value in the spatial distribution, reflecting the disaster diffusion path.

[0163] The spatio-temporal evolution trend map is a visual prediction map combining time series and spatial distribution, showing the dynamic evolution trend of disaster risk.

[0164] In the embodiment of the present application, based on the spatial distribution data of the slip energy accumulation value (such as the energy value of 20 J / m² in the upper reaches of the gully and 10 J / m² in the lower reaches), the gradient field analysis algorithm is used to extract the energy gradient change direction (such as propagating southeastward), and combined with the rock mass fracture propagation rate parameters pre-stored in the dynamic feature library (such as 0.5 m / h in the loose layer), the evolution trend map for the next 6 hours is generated through spatio-temporal interpolation (such as the red highlighted area expanding 500 m southeastward), providing a decision basis for adjusting the coverage range of the early warning signal.

[0165] The following is a specific example:

[0166] In a debris flow gully in the southwest, the slip energy accumulation value reaches the energy saturation threshold of 15 J / m² in the loose layer (step 501), the system extends the accumulation period from 30 minutes to 2 hours, and the calculation frequency is reduced to 1 time per hour; based on the energy spatial distribution data (20 J / m² in the upper reaches and 10 J / m² in the lower reaches), the energy gradient change direction is extracted as southeastward (step 502), and combined with the fracture propagation rate of 0.5 m / h to generate a spatio-temporal evolution trend map, predicting that the high-risk area will expand 300 m southeastward in the next 3 hours, guiding the improvement of the early warning signal intensity of the downstream village to 35 dBm.

[0167] In summary, steps 501 to 502 reduce the system calculation load by 40% through the period adjustment mechanism triggered by the energy saturation threshold; the spatio-temporal evolution trend map improves the prediction accuracy of the disaster diffusion path to 85%, increases the early warning response lead time by 1.5 hours, and overall optimizes the emergency resource allocation efficiency.

[0168] In order to improve the accuracy of mountain flood disaster deformation monitoring and the reliability of risk assessment, by fusing geological zoning data and radar scanning technology, dynamically optimize the radar beam parameters and accurately solve the displacement vector, and combine the deformation critical value of the lithology zoning to generate a terrain risk heat map to realize the spatial visualization representation of disaster risk.

[0169] In some embodiments, in step 102, the ground surface micro-deformation within the sensitive area is periodically scanned by a radar array, the mountain displacement rate is calculated based on Doppler frequency shift measurement, and the displacement rate is dynamically matched with the pre-stored deformation critical value in the dynamic feature library to generate a terrain risk heat map including deformation gradient and displacement direction, including:

[0170] 601. Based on the pre-stored geological partition digital elevation model in the dynamic feature library, determine the scanning coverage range of each radar node in the radar array, dynamically adjust the incident angle of the radar beam according to the terrain undulation characteristics, so that the angle between the beam center line and the surface normal vector is less than a preset threshold, and obtain the Doppler echo signal;

[0171] In step 601, the geological partition digital elevation model refers to a three-dimensional terrain data model divided according to geological characteristics such as lithology and slope, and is used to determine radar scanning parameters.

[0172] The incident angle of the radar beam refers to the angle between the radar beam center line and the surface normal vector, which affects the deformation monitoring accuracy.

[0173] In the embodiments of the present application, based on the pre-stored geological partition digital elevation model in the dynamic feature library (such as the slope angle of 35° in the loose layer area and 15° in the bedrock area), calculate the terrain coverage range of each radar node (such as the coverage radius of node A is 500m), and use beamforming technology to dynamically adjust the incident angle of the radar beam (such as adjusting to 25° in the loose layer area), so that the angle between the beam center line and the surface normal vector is less than 10°, ensuring that the signal-to-noise ratio of the echo signal ≥ 20dB.

[0174] 602. Perform phase demodulation on the Doppler echo signal received by each radar node, extract the phase change amount caused by the ground surface micro-deformation, and calculate the displacement rate vector of each scanning point in combination with the radar wavelength and incident angle parameters;

[0175] In step 602, phase demodulation refers to the process of extracting the phase change amount caused by ground surface deformation from the radar echo signal.

[0176] The displacement rate vector refers to a three-dimensional vector including the magnitude and direction of the displacement rate, which reflects the mountain movement characteristics.

[0177] In the embodiments of the present application, perform IQ quadrature demodulation on the Doppler echo signal received by the radar node, extract the phase change amount Δφ (such as 0.8π radians), and calculate the displacement rate vector (such as 5.2mm / h, direction 10° south by southeast) in combination with the radar wavelength λ (such as 3.9mm) and the incident angle θ (such as 25°), where T is the scanning period (such as 30 minutes).

[0178] 603. Project the displacement rate vector onto the regional main slip direction reference plane pre-stored in the dynamic feature library, separate the slip component along the tangent direction of the reference plane and the deformation component perpendicular to the reference plane, and retain the slip component.

[0179] In step 603, the main slip direction reference plane refers to the reference plane of the main sliding direction of the mountain determined according to historical slip data and geological features.

[0180] The slip component refers to the displacement component along the tangent direction of the reference plane, which reflects the kinetic energy of the disaster-related slip.

[0181] In the embodiment of the present application, project the displacement rate vector onto the main slip direction reference plane (such as a 30° inclination in the southeast direction) pre-stored in the dynamic feature library, calculate the component along the tangent direction of the reference plane (such as 4.8 mm / h) and the vertical component (such as 0.4 mm / h) through vector decomposition, retain the slip component for deformation gradient calculation, and eliminate the non-disaster-related interference of the vertical component.

[0182] 604. According to the deformation critical values corresponding to different lithology partitions in the dynamic feature library, perform regional processing on the slip component to generate the deformation gradient values of each scanning point, and draw a terrain risk heat map in combination with the displacement direction angle.

[0183] In step 604, the deformation gradient value refers to a quantitative index reflecting the severity of deformation after normalization processing, which is determined by the ratio of the slip component to the deformation critical value.

[0184] The terrain risk heat map is a visual map that uses color depth and arrows to indicate the deformation gradient and direction.

[0185] In the embodiment of the present application, according to the deformation critical values of different lithology partitions in the dynamic feature library (such as 5 mm / h for the loose layer), perform normalization processing on the slip component (such as a gradient value of 0.96 corresponding to 4.8 mm / h), generate a deformation gradient field through a spatial interpolation algorithm, and superimpose displacement direction arrows (such as 10° in the southeast direction), and finally output a terrain risk heat map (the gradient value in the red-highlighted area ≥ 0.8).

[0186] The following is a specific example:

[0187] In a debris flow gully in the southwest, the dynamic feature library calls the geological zoning digital elevation model (loose layer slope angle 35°), determines the radar node coverage radius of 500 m, and adjusts the beam incident angle to 25° (step 601); after receiving and demodulating the echo signal, calculates the displacement rate vector of 5.2 mm / h (10° south by southeast) (step 602); projects it onto the main slip reference plane (30° southeast) to separate the slip component of 4.8 mm / h (step 603); generates a deformation gradient value of 0.96 according to the loose layer critical value of 5 mm / h, and the area in the heat map is highlighted in red (step 604).

[0188] In summary, steps 601 to 604 optimize the radar scanning parameters through geological zoning data, improving the deformation monitoring accuracy to 0.5 mm; the projection of the slip component eliminates 35% of the non-disaster interference data; the spatial resolution of the terrain risk heat map reaches 10 m × 10 m, and the false alarm rate is reduced by 28%, providing a high-confidence decision-making basis for emergency response.

[0189] To improve the accuracy of slip component extraction and the accuracy of disaster-related displacement identification, by constructing the main slip direction reference plane and establishing a local projection coordinate system, combining the slip angle threshold of lithological zoning and spatial continuity detection, precise screening and noise elimination of disaster-related slip components are achieved.

[0190] In some embodiments, in step 603, projecting the displacement rate vector onto the regional main slip direction reference plane pre-stored in the dynamic feature library, separating the slip component along the tangent direction of the reference plane and the deformation component perpendicular to the reference plane, and retaining the slip component, includes:

[0191] 701. Based on the geological zoning digital elevation model pre-stored in the dynamic feature library, extract the main slip direction reference plane at the location of each grid cell, and establish a local projection coordinate system with the normal vector of the reference plane as the vertical axis;

[0192] In step 701, the main slip direction reference plane refers to the reference plane of the main sliding direction of the mountain determined according to the historical slip trajectory and geological zoning characteristics.

[0193] The local projection coordinate system refers to the coordinate system with the normal vector of the reference plane as the vertical axis, which is used for the projection and decomposition of the displacement vector.

[0194] In the embodiments of the present application, based on the geological zoning digital elevation model pre-stored in the dynamic feature library (such as the loose layer area slope angle of 35° and the historical slip trajectory southeast), the main slip direction reference plane (such as the 30° inclination angle southeast) is generated by least squares fitting, and a local projection coordinate system (the X-axis is along the tangent direction of the reference plane, and the Y-axis is perpendicular to the reference plane) is established to provide a geometric reference for the decomposition of the displacement vector.

[0195] 702. Based on the local projection coordinate system, decompose the displacement rate vector into a component along the tangent direction of the reference plane and a component perpendicular to the reference plane, and calculate the angle between the tangent direction component and the reference line of the main slip direction as the slip consistency determination parameter.

[0196] In step 702, the slip consistency determination parameter refers to the angle between the tangent direction component of the displacement vector and the reference line of the main slip direction, and is used to screen the displacements related to disasters.

[0197] In the embodiment of the present application, project the displacement rate vector (such as 6.3 mm / h, direction 15° south by southeast) onto the local coordinate system, calculate the component along the tangent direction of the reference plane (such as 5.8 mm / h) and the perpendicular component (such as 0.5 mm / h), and measure the angle (such as 15°) between the tangent direction component and the reference line of the main slip direction (30° southeast) as the quantification index of slip consistency.

[0198] 703. Based on the slip consistency determination parameter, according to the slip angle threshold of the lithology partition stored in the dynamic feature library, screen the displacement vectors in the tangent direction component whose angle with the reference line of the main slip direction is less than the threshold as the effective slip components.

[0199] In step 703, the slip angle threshold of the lithology partition refers to the maximum allowable deviation of the slip direction set according to different lithologies, and is used to filter abnormal displacements.

[0200] In the embodiment of the present application, according to the slip angle threshold of the lithology partition stored in the dynamic feature library (such as 20° for loose layers and 15° for bedrock), screen the displacement vectors in the tangent direction component whose angle is less than the threshold (such as the 15° in the loose layer area meets the conditions), eliminate the over-limit data (such as 45° north), and retain the effective slip components (such as 10° southeast).

[0201] 704. Perform spatial continuity detection on the effective slip components, eliminate the isolated displacement signals caused by local collapses or rockfall activities, and retain the slip components that are consistent with the displacement directions of adjacent grid cells.

[0202] In step 704, the spatial continuity detection refers to a screening mechanism based on the consistency of the displacement directions of adjacent grids to eliminate isolated noise signals.

[0203] In the embodiment of the present application, perform neighborhood analysis (such as 3×3 grid) on the effective slip components, calculate the displacement direction consistency index of adjacent cells (such as the proportion of southeast ±10° ≥ 70%), eliminate the isolated signals (such as a single cell with 45° north), and retain the continuously distributed slip components (such as 8° southeast) to ensure the spatial continuity of the data.

[0204] The following is a specific example:

[0205] In a debris flow gully in the southwest, the dynamic feature library extracts the reference plane of the main sliding direction of the loose layer area (30° southeast), and establishes a local projection coordinate system (step 701); decomposes the displacement vector (6.3 mm / h, 15° south by southeast) into a tangential component of 5.8 mm / h (at an angle of 15°) and a vertical component of 0.5 mm / h (step 702); according to the loose layer sliding angle threshold of 20°, retains the effective sliding component (step 703); eliminates the isolated northward 45° signal through spatial continuity detection, and retains the continuous sliding component of 8° southeast (step 704).

[0206] In summary, steps 701 to 704 eliminate 35% of the non-disaster displacement interference through local coordinate system projection; the sliding angle threshold screening improves the recognition accuracy of the effective component to 92%; the spatial continuity detection further reduces the misjudgment rate by 18%, and finally realizes the high-precision extraction of the disaster-related sliding component (error <5%), providing a reliable data basis for risk assessment.

[0207] Figure 2 The present application provides a structural schematic diagram of an emergency communication device (or system) based on a multi-modal early warning broadcast terminal, as Figure 2 shown, the device includes:

[0208] A disaster trigger monitoring module 21, which is used to define the combination of the soil humidity threshold and the rainfall intensity threshold as the disaster trigger condition in a preset sensitive area of mountain flood disasters. When the real-time monitoring data meets the disaster trigger condition, it activates the directional early warning broadcast function of the multi-modal terminal deployed in the sensitive area, and establishes a dynamic feature library to store the debris flow association rule set;

[0209] A deformation monitoring and analysis module 22, which is used to periodically scan the surface micro-deformation in the sensitive area through a radar array, calculate the mountain displacement rate based on Doppler frequency shift measurement, and dynamically match the displacement rate with the deformation critical value pre-stored in the dynamic feature library to generate a terrain risk heat map including the deformation gradient and the displacement direction;

[0210] An emergency communication networking module 23, which is used to build a self-organizing communication link based on the low-power long-distance wireless mesh network protocol in the sensitive area, starts a multi-hop redundant transmission mechanism when the public network signal is interrupted, divides the terrain risk heat map and the early warning broadcast instruction into data packet groups, and performs multi-hop transmission through random relay nodes;

[0211] The hierarchical response scheduling module 24 is used to dynamically allocate the forwarding priorities of the data packet groups according to the early warning levels jointly determined by the dynamic feature library and the radar array, so that the audible and visual early warning signals of the multimodal terminal synchronously match the coverage ranges according to the disaster risk gradient, forming a hierarchical response link for terrain deformation monitoring and emergency communication linkage.

[0212] Figure 2 The described emergency communication device based on a multimodal early warning broadcast terminal can execute Figure 1 The described emergency communication method based on a multimodal early warning broadcast terminal in the illustrated embodiment, its implementation principle and technical effects will not be elaborated. For the emergency communication device based on a multimodal early warning broadcast terminal in the above embodiment, the specific manners in which each module and unit perform operations have been described in detail in the embodiments related to the method, and will not be elaborated here.

[0213] In a possible design, Figure 2 The emergency communication device based on a multimodal early warning broadcast terminal in the illustrated embodiment can be implemented as a computing device, such as Figure 3 As shown, the computing device can include a storage component 31 and a processing component 32;

[0214] The storage component 31 stores one or more computer instructions, wherein the one or more computer instructions are called and executed by the processing component 32.

[0215] The processing component 32 is used for the Figure 1 Emergency communication method based on a multimodal early warning broadcast terminal in the above-described embodiment.

[0216] Among them, the processing component 32 can include one or more processors to execute computer instructions to complete all or part of the steps in the above method. Of course, the processing component can also be implemented by one or more application-specific integrated circuits (ASICs), digital signal processors (DSPs), digital signal processing devices (DSPDs), programmable logic devices (PLDs), field programmable gate arrays (FPGAs), controllers, microcontrollers, microprocessors or other electronic components for executing the above method.

[0217] The storage component 31 is configured to store various types of data to support the operation of the terminal. The storage component can be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic memory, flash memory, magnetic disk or optical disk.

[0218] Of course, the computing device necessarily may also include other components, such as input / output interfaces, display components, communication components, etc.

[0219] The input / output interface provides an interface between the processing component and the peripheral interface module, and the above-mentioned peripheral interface module may be an output device, an input device, etc.

[0220] The communication component is configured to facilitate communication between the computing device and other devices in a wired or wireless manner, etc.

[0221] Among them, the computing device may be a physical device or an elastic computing host provided by a cloud computing platform, etc. At this time, the computing device may refer to a cloud server, and the above-mentioned processing component, storage component, etc. may be basic server resources leased or purchased from a cloud computing platform.

[0222] The embodiment of the present application also provides a computer storage medium storing a computer program, and when the computer program is executed by a computer, it can implement the above-mentioned Figure 1 An emergency communication method based on a multi-modal early warning broadcast terminal shown in the above embodiment.

[0223] Those skilled in the art can clearly understand that for the convenience and brevity of description, the specific working processes of the above-described systems, devices, and units can refer to the corresponding processes in the foregoing method embodiments, and will not be described herein again.

[0224] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place, or may be distributed to multiple network units. Some or all of the modules can be selected according to actual needs to achieve the purpose of the solution of this embodiment. Those of ordinary skill in the art can understand and implement it without creative labor.

[0225] Through the description of the above embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus a necessary general hardware platform, and of course, it can also be implemented by hardware. Based on such an understanding, the essence of the above technical solution, or the part that contributes to the prior art, can be embodied in the form of a software product. The computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disc, etc., and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the methods described in each embodiment or some parts of the embodiments.

[0226] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present application, rather than limiting them; although the present application has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions described in the foregoing embodiments, or perform equivalent replacements on some of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the various embodiments of the present application.

Claims

1. An emergency communication method based on a multimodal early warning broadcast terminal, characterized in that: include: In the preset flash flood disaster sensitive area, the combination of soil moisture threshold and rainfall intensity threshold is defined as the disaster triggering condition. When the real-time monitoring data meets the disaster triggering condition, the directional warning broadcast function of the multimodal terminal deployed in the sensitive area is activated, and a dynamic feature library is established to store the debris flow association rule set; Periodically scanning the surface micro-deformation in the sensitive area through a radar array, calculating the displacement rate of the mountain based on Doppler frequency shift measurement, dynamically matching the displacement rate with the deformation critical value pre-stored in the dynamic feature library, and generating a terrain risk heat map including deformation gradient and displacement direction; A self-organizing communication link based on a low-power long-distance wireless mesh network protocol is constructed in the sensitive area, and a multi-hop redundant transmission mechanism is activated when the public network signal is interrupted, and the terrain risk heat map and the early warning broadcast instruction are divided into data packet groups and multi-hop transmission is performed through random relay nodes; According to the warning level determined jointly by the dynamic feature library and the radar array, the forwarding priority of the data packet group is dynamically allocated, so that the sound and light warning signals of the multimodal terminal are synchronously matched with the coverage range according to the disaster risk gradient, forming a hierarchical response link of terrain deformation monitoring and emergency communication linkage.

2. The method according to claim 1, characterized in that According to the warning level determined by the dynamic feature library and the radar array, the forwarding priority of the data packet group is dynamically allocated, so that the sound and light warning signals of the multimodal terminal are synchronously matched with the coverage according to the disaster risk gradient, forming a hierarchical response link of terrain deformation monitoring and emergency communication linkage, including: Based on the deformation gradient distribution of the terrain risk heat map, the displacement rate and main slip direction in the grid unit are extracted, and the corresponding relationship between the geological partition data and the deformation critical value in the dynamic feature library is combined to generate a risk quantification parameter set containing displacement rate interval and direction characteristics; Based on the risk quantification parameter set, the displacement rate of each grid unit is decomposed into vertical slope and along-slope components according to the mountain slope, and the cumulative change of the along-slope component in each orientation area is calculated. According to the correspondence between the cumulative change and the warning level pre-stored in the dynamic feature library, the warning level of each area and the corresponding signal coverage strength are determined; Calculating a distribution density threshold of relay nodes in the self-organizing communication link according to the signal coverage strength, and establishing a multi-hop forwarding path sequence with the opposite direction of displacement as the transmission path based on the positional relationship between the main sliding direction and the relay nodes; The warning broadcast instruction is encapsulated into a data packet group carrying a warning level identifier and a path constraint condition, and the radio frequency transmission intensity is adjusted according to the terrain shielding degree between nodes in the multi-hop forwarding path sequence, so that the sound and light warning signal forms a coverage intensity field that attenuates from the core area to the periphery in the area in the opposite direction of the main sliding direction; Based on the coverage intensity field, when the radar array detects that the deviation of the main slip direction of adjacent periods exceeds the direction tolerance value set by the dynamic feature library, it triggers the adaptive adjustment of the deformation critical value in the dynamic feature library, and synchronously updates the multi-hop forwarding path sequence and the radio frequency transmission intensity, forming a real-time response mechanism for the linkage between terrain displacement changes and communication parameters.

3. The method according to claim 2, characterized in that Based on the risk quantification parameter set, the displacement rate of each grid unit is decomposed into vertical slope and along-slope components according to the mountain slope, and the cumulative change of the along-slope component in each orientation area is calculated. According to the correspondence between the cumulative change and the warning level pre-stored in the dynamic feature library, the warning level of each area and the corresponding signal coverage strength are determined, including: Based on the geological zoning digital elevation model pre-stored in the dynamic feature library, the mountain slope reference plane at the location of each grid unit is extracted, and a local coordinate system with the normal vector of the slope reference plane as the vertical axis is established; Projecting the displacement rate vector in the grid unit to the local coordinate system, separating the displacement component along the tangent direction of the slope reference plane as the sliding driving force characterization quantity, and retaining the displacement component perpendicular to the slope reference plane as the deformation stretching quantity; The sliding driving force characterization quantity in the continuous monitoring period is screened for directional consistency, and the displacement components whose angle with the main sliding direction is less than the threshold set by the dynamic feature library are retained for cumulative calculation, and the abnormal deflection components caused by local collapse are eliminated; Based on the accumulated sliding period parameters corresponding to different lithology partitions in the dynamic feature library, the screened sliding driving force characterization quantity is adaptively accumulated in a time window to generate accumulated sliding energy values ​​in various orientation areas; According to the lithology and energy level mapping table pre-stored in the dynamic feature library, the accumulated value of the slip energy is converted into the warning level calibration coefficient of the corresponding partition, and the spatial attenuation gradient of the signal coverage intensity is calculated in combination with the terrain shielding factor.

4. The method according to claim 3, characterized in that Based on the accumulated sliding cycle parameters corresponding to different lithology partitions in the dynamic feature library, the screened sliding driving force characterization quantity is adaptively accumulated in a time window to generate the accumulated sliding energy values ​​in various orientation areas, including: According to the mechanical parameter table of rock type partition pre-stored in the dynamic characteristic library, the rheological characteristic parameters of the rock mass in the area where each grid unit is located are extracted, wherein the rheological characteristic parameters of the rock mass include the viscoelastic coefficient and the critical slip rate threshold, and the length of the initial slip cumulative cycle is determined; Based on the length of the initial slip accumulation cycle, the filtered slip driving force characterization quantity is accumulated in sections, the slip displacement integral value in each time window is calculated, and the change trend of the vertical slope component is monitored at the same time; When the rate of change of the vertical slope component in a single time window exceeds the deformation mutation threshold set by the dynamic feature library, the current slip accumulation period is shortened to a preset lower limit value, and the slip displacement integral value is recalculated; According to the integral value of the slip displacement and the friction coefficient of the corresponding lithology partition, the equivalent slip energy value in each time window is calculated, and the equivalent slip energy values ​​of the continuous time windows are superimposed according to the exponential decay weight to generate the slip energy accumulation value.

5. The method according to claim 4, characterized in that Also includes: When the accumulated value of the slip energy reaches the saturation threshold of the lithology partition energy stored in the dynamic feature library, the accumulated slip period is extended to a preset upper limit value, and the calculation frequency of the slip displacement integral value is reduced; Based on the spatial distribution characteristics of the accumulated sliding energy value, the energy gradient change direction in each azimuth region is extracted, and combined with the rock mass fracture propagation rate parameters pre-stored in the dynamic feature library, a spatiotemporal evolution trend diagram of the accumulated sliding energy value is generated.

6. The method according to claim 1, characterized in that The surface micro-deformation in the sensitive area is periodically scanned by a radar array, the displacement rate of the mountain is calculated based on Doppler frequency shift measurement, the displacement rate is dynamically matched with the deformation critical value pre-stored in the dynamic feature library, and a terrain risk heat map containing deformation gradient and displacement direction is generated, including: Based on the geological partition digital elevation model pre-stored in the dynamic feature library, the scanning coverage range of each radar node in the radar array is determined, and the radar beam incident angle is dynamically adjusted according to the terrain undulation characteristics, so that the angle between the beam centerline and the surface normal vector is less than a preset threshold, and a Doppler echo signal is obtained; Performing phase demodulation on the Doppler echo signal received by each radar node, extracting the phase change caused by surface micro-deformation, and calculating the displacement rate vector of each scanning point by combining the radar wavelength and incident angle parameters; Projecting the displacement rate vector onto the regional main slip direction reference plane pre-stored in the dynamic feature library, separating the slip component along the tangent direction of the reference plane and the deformation component perpendicular to the reference plane, and retaining the slip component; According to the deformation critical values ​​corresponding to different lithology partitions in the dynamic feature library, the slip component is regionally processed to generate the deformation gradient value of each scanning point, and the terrain risk thermodynamic map is drawn in combination with the displacement direction angle.

7. The method according to claim 6, characterized in that Projecting the displacement rate vector onto the regional main slip direction reference plane pre-stored in the dynamic feature library, separating the slip component along the tangent direction of the reference plane and the deformation component perpendicular to the reference plane, and retaining the slip component, including: Based on the digital elevation model of geological partitions pre-stored in the dynamic feature library, the main slip direction reference plane of each grid unit is extracted, and a local projection coordinate system with the normal vector of the reference plane as the vertical axis is established; Based on the local projection coordinate system, the displacement rate vector is decomposed into a component along the tangent direction of the reference plane and a component perpendicular to the reference plane, and the angle between the tangent direction component and the main slip direction reference line is calculated as a sliding consistency judgment parameter; Based on the slip consistency judgment parameter and the lithology partition slip angle threshold value pre-stored in the dynamic feature library, the displacement vectors in the tangent direction component whose angle with the main slip direction reference line is less than the threshold value are selected as effective slip components; The effective slip component is tested for spatial continuity, isolated displacement signals caused by local collapse or rock rolling are eliminated, and slip components with the same displacement direction as adjacent grid units are retained.

8. An emergency communication system based on a multi-modal early warning broadcast terminal, characterized in that: include: The disaster trigger monitoring module is used to define the combination of soil moisture threshold and rainfall intensity threshold as disaster trigger conditions in a preset flash flood disaster sensitive area. When the real-time monitoring data meets the disaster trigger conditions, the directional warning broadcast function of the multimodal terminal deployed in the sensitive area is activated, and a dynamic feature library is established to store a set of debris flow association rules; A deformation monitoring and analysis module is used to periodically scan the surface micro-deformation in the sensitive area through a radar array, calculate the mountain displacement rate based on Doppler frequency shift measurement, dynamically match the displacement rate with the deformation critical value pre-stored in the dynamic feature library, and generate a terrain risk heat map including deformation gradient and displacement direction; An emergency communication networking module is used to construct a self-organizing communication link based on a low-power long-distance wireless mesh network protocol in the sensitive area, and to start a multi-hop redundant transmission mechanism when the public network signal is interrupted, so as to divide the terrain risk heat map and the early warning broadcast instruction into data packet groups and perform multi-hop transmission through random relay nodes; The hierarchical response scheduling module is used to dynamically allocate the forwarding priority of the data packet group according to the warning level jointly determined by the dynamic feature library and the radar array, so that the sound and light warning signals of the multimodal terminal can be synchronously matched with the coverage range according to the disaster risk gradient, forming a hierarchical response link for terrain deformation monitoring and emergency communication linkage.

9. A computing device, characterized in that It comprises a processing component and a storage component; the storage component stores one or more computer instructions; the one or more computer instructions are used to be called and executed by the processing component to implement an emergency communication method based on a multimodal early warning broadcast terminal as described in any one of claims 1 to 7.

10. A computer storage medium, characterized in that: A computer program is stored, and when the computer program is executed by a computer, an emergency communication method based on a multimodal early warning broadcast terminal as described in any one of claims 1 to 7 is implemented.

Images