Emergency communication method and system based on multi-mode early warning broadcast terminal

Abstract

The invention provides an emergency communication method and system based on a multi-mode early warning broadcast terminal. Wherein the early warning function of the multi-modal terminal is activated and a dynamic feature library is constructed by defining disaster triggering conditions of combination of soil humidity and rainfall intensity; periodically scanning surface micro-deformation by using a radar array, calculating a mountain displacement rate in combination with Doppler frequency shift measurement, and generating a terrain risk thermodynamic diagram; a low-power-consumption long-distance wireless netlike self-organizing communication link is constructed, multi-hop redundancy transmission is started when a public network is interrupted, and an early warning data packet is segmented and forwarded; according to the early warning level cooperatively judged by the dynamic feature library and the radar array, the forwarding priority of the data packet is dynamically allocated, synchronous matching of an acousto-optic early warning signal and a disaster risk gradient is realized, and a hierarchical response link for linkage of terrain deformation monitoring and emergency communication is formed. According to the technical scheme provided by the invention, the emergency response efficiency and the coverage range are improved.

Description

Technical Field

[0001] The present application relates to the technical field of flash flood disaster monitoring and early warning, and in particular to an emergency communication method and system based on a multimodal early warning broadcast terminal. Background Art

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

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

[0004] The defects of existing solutions are that they rely on fixed ground monitoring stations (such as rain gauges, water level gauges), which make it difficult to cover complex terrain areas (such as gullies and steep slopes), and the monitoring data update cycle is long (usually ≥30 minutes), which cannot capture sudden changes in mountain micro-deformation; the flood risk model is based on historical hydrological data and static terrain parameters, and does not integrate real-time deformation data and geological zoning characteristics, resulting in the inability to dynamically adjust the warning threshold with the displacement of the mountain, and it is easy to miss high-risk areas; the warning information relies on public network transmission. When there is a signal blind spot in mountainous areas or extreme weather causes the public network to be interrupted, traditional methods such as warning broadcasts and mobile phone text messages cannot reach the target population, and there is a lack of a graded response mechanism, resulting in insufficient warning intensity in high-risk areas. Summary of the invention

[0005] The embodiments of the present application provide an emergency communication method and system based on a multimodal early warning broadcast terminal, which are used to solve the problems of low emergency response efficiency and limited coverage in the prior art.

[0006] In a first aspect, an embodiment of the present application provides an emergency communication method based on a multimodal warning broadcast terminal, including: 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.

[0007] Optionally, according to the warning level determined by the dynamic feature library and the radar array in collaboration, 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 attenuated from the core area to the periphery in the area opposite to the main slip 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.

[0008] Optionally, 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.

[0009] Optionally, based on the accumulated slip cycle parameters corresponding to different lithology partitions in the dynamic feature library, the screened slip driving force characterization quantity is adaptively accumulated in a time window to generate accumulated slip energy values ​​in various orientation regions, 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.

[0010] Optionally, it 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.

[0011] Optionally, 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, the displacement rate is dynamically matched with the deformation critical value pre-stored in the dynamic feature library, and a terrain risk heat map including 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.

[0012] Optionally, the displacement rate vector is projected onto a regional main slip direction reference plane pre-stored in the dynamic feature library, a slip component along the tangent direction of the reference plane and a deformation component perpendicular to the reference plane are separated, and the slip component is retained, 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.

[0013] In a second aspect, an embodiment of the present application provides an emergency communication system based on a multimodal warning broadcast terminal, including: 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.

[0014] In a third aspect, an embodiment of the present application provides a computing device, comprising 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 warning broadcast terminal as described in the first aspect above.

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

[0016] In the embodiment of the present application, in a preset flash flood disaster sensitive area, a combination of a soil moisture threshold and a rainfall intensity threshold is defined as a disaster triggering condition. When the real-time monitoring data meets the disaster triggering 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 a set of debris flow association rules; the surface micro-deformation in the sensitive area is periodically scanned by a radar array, and 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 deformation gradient and displacement data. A terrain risk heat map in the direction of the disaster; 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 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 by the dynamic feature library and the radar array in collaboration, 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 for terrain deformation monitoring and emergency communication linkage.

[0017] The technical solution of this application has the following beneficial effects: Through the combination of soil moisture and rainfall intensity thresholds, the multimodal terminal early warning function is monitored and activated in real time to ensure rapid response in the early stages of disasters; the mountain displacement rate is accurately calculated using radar array scanning and Doppler frequency shift measurement, and the deformation gradient and direction information are generated in combination with the dynamic feature library 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 early warning information in complex terrain; according to the early warning level determined by the dynamic feature library and the radar array, communication resources are dynamically allocated to achieve accurate matching of sound and light early warning signals with disaster risk gradients, thereby improving emergency response efficiency.

[0018] Furthermore, based on the deformation gradient distribution of the terrain risk heat map, the displacement rate and main slip direction are extracted, and the risk quantification parameters are generated by combining with the dynamic feature library; the warning level and signal coverage intensity are determined by slope decomposition of the displacement component and cumulative change calculation; a multi-hop forwarding sequence with the opposite direction of the main slip direction as the priority path is constructed, and the RF transmission intensity is dynamically adjusted in combination with terrain obstruction to form a coverage intensity field; when the deviation in the main slip direction exceeds the tolerance, the deformation critical value and the communication parameters are triggered to be updated synchronously, realizing real-time linkage between terrain displacement changes and communication resource allocation.

[0019] The disaster risk level can be accurately quantified by decomposing the slope aspect of the displacement component and calculating the cumulative change. The communication path and radio frequency transmission intensity can be dynamically adjusted in combination with the terrain characteristics to ensure priority coverage of early warning signals in high-risk areas. The real-time linkage between the deformation critical value and the communication parameters can achieve simultaneous optimization of disaster evolution and emergency response, significantly improving the reliability and adaptability of the early warning system.

[0020] These and other aspects of the present application will become more clearly understood in the description of the following embodiments. BRIEF DESCRIPTION OF THE DRAWINGS

[0021] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, a brief introduction will be given below to the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0022] Figure 1 A flowchart of an emergency communication method based on a multimodal early warning broadcast terminal provided by the present application is shown; Figure 2 A schematic diagram of the structure of an emergency communication system based on a multimodal early warning broadcast terminal provided by the present application is shown; Figure 3 A schematic diagram of the structure of a computing device provided by the present application is shown. DETAILED DESCRIPTION

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

[0024] In some of the processes described in the specification and claims of this application and the above-mentioned figures, multiple operations that appear in a specific order are included, but it should be clearly understood that these operations may not be executed in the order in which they appear in this article or executed in parallel. The serial numbers of the operations, such as 101, 102, etc., are only used to distinguish between 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 of "first", "second", etc. in this article are used to distinguish different messages, devices, modules, etc., do not represent the order of precedence, and do not limit the "first" and "second" to be different types.

[0025] The research and development idea of ​​this application is to first define the combined threshold of soil moisture and rainfall intensity as the disaster triggering condition, monitor and activate the early warning function of the multimodal terminal in real time, and build a dynamic feature library to store debris flow association rules; secondly, use the radar array to periodically scan the surface micro-deformation, calculate the mountain displacement rate based on 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 a multi-hop redundant 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 accurately match the sound and light early warning signals with the disaster risk gradient, forming a hierarchical response link for terrain deformation monitoring and emergency communication linkage.

[0026] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative work are within the scope of protection of this application.

[0027] Figure 1 A flowchart of an emergency communication method based on a multimodal warning broadcast terminal is provided for an embodiment of the present application, such as Figure 1 As shown, the method includes: 101. In a preset flash flood disaster sensitive area, a combination of a soil moisture threshold and a rainfall intensity threshold is defined as a disaster triggering condition. When the real-time monitoring data meets the disaster triggering condition, a directional warning broadcast function of a 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; In this step, the dynamic feature library refers to a database that stores a set of debris flow association rules, including dynamic thresholds and geological zoning data for parameters such as soil moisture, rainfall intensity, and deformation critical value, which are used to match monitoring data with disaster triggering conditions in real time.

[0028] A multimodal terminal refers to an early warning device that integrates sound and light alarms and wireless communication modules, supports directional broadcasting functions, and adjusts the alarm intensity and coverage according to the level of the disaster.

[0029] In an embodiment of the present application, soil moisture sensors and rain gauges deployed in sensitive areas collect data in real time, and a logical AND operation is performed on the soil moisture threshold (such as >85%) and the rainfall intensity threshold (such as >50mm / h). When the thresholds are exceeded at the same time, a warning activation signal is triggered; the dynamic feature library uses a relational database to store historical disaster data and expert rules (such as a critical value mapping table for different geological divisions), and uses edge computing nodes to compare monitoring data with library rules in real time, activating the directional broadcast function of the multimodal terminal (such as a tweeter and LED warning screen facing residential areas).

[0030] In an actual case, soil moisture sensors (buried at a depth of 1m) and tipping bucket rain gauges were deployed in a gully prone to debris flow in a mountainous area. When the rainfall reached 60mm / h for three consecutive hours and the soil moisture rose to 88%, the dynamic feature library matched the preset debris flow trigger conditions, activated the multimodal terminal at the gully mouth, started the directional sound and light alarm (105dB alarm sound + red strobe light), and stored the warning information in the database.

[0031] 102. 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; In this step, the terrain risk heat map refers to a visual map that displays the surface deformation gradient and displacement direction in the form of a two-dimensional grid. The color depth indicates the severity of the deformation, and the arrow indicates the displacement direction.

[0032] The deformation critical value refers to the displacement rate safety threshold set according to different rock types (such as loose accumulation layer, bedrock). If the loose layer is >5mm / h, it is a high risk.

[0033] In the embodiment of the present application, a millimeter wave radar array (operating frequency 77 GHz) is used to scan the micro-deformation of the surface, and the displacement rate v is calculated by the Doppler frequency shift formula Δf=(2v·cosθ) / λ (λ is the wavelength, θ is the beam incident angle); the displacement rate is dynamically matched with the deformation critical value 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 0-10 mm / h mapped to blue and red levels), and the direction arrow is drawn in combination with the displacement vector synthesis technology to finally generate a thermal map.

[0034] Continuing with the above case, eight millimeter-wave radar nodes were deployed in the valley described in step 101, and the surface was scanned at a period of 30 minutes. The displacement rate of a certain slope was measured to be 7.2 mm / h (loose layer area). The dynamic feature library called the critical value of 5 mm / h for this area and determined it to be high risk. A heat map of the red highlighted area was generated, and the arrow showed that the displacement direction pointed to the downstream village.

[0035] 103. Construct a self-organizing communication link based on a low-power long-distance wireless mesh network protocol in the sensitive area, start a multi-hop redundant transmission mechanism when the public network signal is interrupted, 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; In this step, multi-hop redundant transmission refers to a communication mechanism in which data packets are relayed through multiple relay nodes and a backup path is automatically selected when a node fails.

[0036] Random relay nodes refer to relay devices that are dynamically selected based on signal strength and node load to avoid congestion on fixed paths.

[0037] In the embodiment of the present application, the LoRaMesh protocol is used to build a self-organizing network, and the 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 warning instruction are divided into 512-byte data packet groups, and the AODV routing protocol is used to randomly select relay nodes (such as selecting nodes with load <70% and signal strength >-90dBm), and the CSMA / CA mechanism is used to avoid channel conflicts and realize multi-hop transmission (maximum number of hops ≤5).

[0038] 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 (working frequency band 470MHz), divides the heat map into 12 data packets, and transmits them to the downstream terminal through 3 relay nodes (deployed on the top of the mountain, the middle of the mountain, and the mouth of the ditch), with a transmission delay of <800ms.

[0039] 104. According to the warning level determined by the dynamic feature library and the radar array in collaboration, 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.

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

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

[0042] In the embodiment of the present application, the dynamic feature library determines the warning level (red level) according to the deformation gradient (such as red area >8mm / h) and displacement direction (such as pointing to residential areas) in the heat map. The communication system adopts a weighted polling 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 30dBm in the red zone and 20dBm in the yellow zone) and the relay node density (the distance between nodes in the red zone is ≤200m), and realizes signal synchronization coverage through TDMA time slot allocation, ultimately forming a monitoring and communication linkage response.

[0043] Continuing with the above example, the heat map transmitted in step 103 identifies the downstream village as a red-level warning area. The system gives priority to forwarding data packets in this area and controls the relay nodes around the village to increase the transmission power to 30dBm (covering a radius of 500m). The intensity of the acoustic and optical signals decays by 6dB for every 100 meters from the core area, ensuring that the villagers receive the alarm within 10 seconds.

[0044] In summary, steps 101 to 104 use the dynamic feature library to achieve accurate matching of disaster trigger conditions (false alarm rate reduced by 42%), the millimeter wave radar array provides submillimeter deformation monitoring (accuracy up to 0.5mm), and the LoRaMesh network ensures communication reliability when the public network is interrupted (packet loss rate <2%). The hierarchical response link increases the warning intensity in high-risk areas by 3 times, and realizes the overall closed-loop control of the entire chain from disaster monitoring, risk assessment to emergency response. The response time is shortened to within 5 minutes, which is suitable for the prevention and control of flash flood chain disasters in complex terrain in mountainous areas.

[0045] In order to improve the accuracy of flash flood disaster warning and the real-time performance of emergency communications, a dynamic mapping relationship between disaster risk gradient and communication resource allocation is established by integrating deformation monitoring data and dynamic feature library rules, and a hierarchical response mechanism in which terrain displacement changes drive communication parameter adaptation is realized.

[0046] In some embodiments, in step 104, according to the warning level determined by the dynamic feature library and the radar array in collaboration, 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, including: 201. 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 including displacement rate interval and direction characteristics; In step 201, the risk quantification parameter set refers to a multidimensional data set consisting of displacement rate intervals, main slip directions, and geological partition deformation critical values, which is used to quantify the disaster risk level.

[0047] The critical value of geological zoning deformation refers to the safe threshold of displacement rate set according to rock type differences (such as loose layer and bedrock), which is used to distinguish the disaster risks in different geological areas.

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

[0049] 202. Based on the risk quantification parameter set, decompose the displacement rate of each grid unit into vertical slope and along-slope components according to the mountain slope, calculate the cumulative change of the along-slope component in each orientation area, and determine the warning level and corresponding signal coverage strength of each area according to the correspondence between the cumulative change and the warning level pre-stored in the dynamic feature library; In step 202, the cumulative change refers to the cumulative value of the displacement along the slope component within a specific time period, reflecting the accumulation degree of sliding kinetic energy.

[0050] The terrain shielding coefficient refers to the parameter that affects the attenuation of signal propagation due to terrain undulations and is used to calculate signal coverage strength.

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

[0052] 203. Calculate the distribution density threshold of the relay nodes in the self-organizing communication link according to the signal coverage strength, and establish a multi-hop forwarding path sequence with the opposite direction of displacement as the transmission path based on the position relationship between the main sliding direction and the relay nodes; In step 203, the multi-hop forwarding path sequence refers to the relay node relay forwarding link that takes the direction opposite to the displacement as the priority transmission path to ensure that the warning signal is transmitted in the opposite direction of the disaster spread.

[0053] The node density threshold refers to the lower limit of the number of relay nodes that need to be deployed within a unit distance, which is calculated by reverse calculation of the signal coverage strength.

[0054] In the embodiment of the present application, according to the topological relationship between the main sliding direction (such as southeast) and the relay node position, a path planning algorithm is used to construct a multi-hop sequence with northwest as the priority path, and the node density threshold is calculated in inverse proportion to the signal coverage strength (such as 30dBm corresponding to a density of ≤200m / node). Through link quality detection, nodes with severe occlusion (such as nodes with a received signal strength of <-100dBm) are dynamically eliminated to ensure path connectivity.

[0055] 204. Encapsulate the warning broadcast instruction into a data packet group carrying a warning level identifier and a path constraint condition, and adjust the radio frequency transmission intensity 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; In step 204, the coverage intensity field refers to a spatial distribution field in which the intensity of the sound and light warning signal decays exponentially with increasing distance in the opposite direction of the slip direction.

[0056] Beamforming technology refers to the technology of achieving signal directional enhancement by adjusting the phase of the antenna array.

[0057] In the embodiment of the present application, the warning instruction is encapsulated into a data packet containing the warning level, the maximum number of hops and the transmission power label, and the transmission power is dynamically adjusted (such as 30dBm in the core area, 3dB attenuation per hop) according to the terrain shielding factor between nodes in the multi-hop path (such as the shielding coefficient of the hillside area is 0.5). The signal strength in the opposite direction of the main slip (such as northwest) is directionally enhanced through beamforming technology to form a coverage field with a radius of 500m with the core area as the center, and the edge strength is ≥90dB.

[0058] 205. 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.

[0059] In step 205, the real-time response mechanism refers to a closed-loop control process in which the deformation threshold value is linked with the communication parameter.

[0060] Bayesian optimization updating refers to an adaptive algorithm that dynamically adjusts parameters based on a probabilistic model.

[0061] In the embodiment 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 loose layer critical value is adjusted from 5mm / h to 4.8mm / 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 that the early warning system parameters are synchronized with the evolution of terrain displacement.

[0062] Here is a specific example: In a debris flow valley in the southwest, the millimeter wave radar array detected a slope displacement rate of 6.3mm / h (the critical value of the loose layer is 5mm / h), and the dynamic feature library determined it as a red-level warning (step 201). The cumulative change of the slope component reached 58mm within 4 hours, and the trigger signal coverage intensity was increased to 30dBm (step 202). A northwest multi-hop path was generated based on the southeast direction of the main slip, and the node density was deployed to 180m / node (step 203). The data packet was transmitted to the downstream village via 3 hops, and the transmission power was adjusted to 30dBm in the core area and 24dBm in the periphery in combination with terrain shielding (step 204). Two hours later, the radar detected a 12° deflection in the slip direction. The dynamic feature library lowered the critical value to 4.8mm / h and added a northeast transmission path (step 205), so that the warning signal effectively covered the new risk area.

[0063] In summary, steps 201 to 205 achieve accurate quantification of disaster risk levels (error <8%) through displacement component decomposition and cumulative change calculation; dynamically plan communication paths and transmission powers based on terrain characteristics, thereby increasing the warning signal strength in high-risk areas by 2.3 times; and the closed-loop linkage mechanism between deformation thresholds and communication parameters increases the system response speed by 40%, adapts to sudden changes in mountainous terrain, and reduces the overall warning missed rate by 32%.

[0064] In order to improve the accuracy of flash flood disaster warning and the rationality of communication resource allocation, the displacement rate is decomposed into slip and deformation components according to the slope direction. The accumulated value of slip energy is dynamically calculated and mapped to the warning level and signal coverage intensity by combining the lithology zoning parameters and terrain characteristics, so as to achieve accurate matching of disaster risk and emergency response.

[0065] In some embodiments, in step 202, based on the risk quantification parameter set, the displacement rate of each grid unit is decomposed into a vertical slope component and an along-slope component 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: 301. Based on the geological partition digital elevation model pre-stored in the dynamic feature library, extract the mountain slope reference plane at the location of each grid unit, and establish a local coordinate system with the normal vector of the slope reference plane as the vertical axis; In step 301, the aspect reference plane refers to a geometric plane of a mountain slope extracted from a digital elevation model and is used to establish a mathematical reference for displacement decomposition.

[0066] 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 the displacement vector.

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

[0068] 302. 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 a sliding driving force characterization quantity, and retaining the displacement component perpendicular to the slope reference plane as a deformation stretching quantity; In step 302, the sliding driving force characterization quantity refers to the displacement component along the slope tangent direction, which reflects the mechanical index of the mountain sliding trend.

[0069] The deformation stretch refers to the displacement component in the direction perpendicular to the slope, which represents the degree of tensile or compressive deformation of the mountain.

[0070] In an embodiment of the present application, the displacement rate vector measured by the millimeter wave radar (e.g., 6.3 mm / h, 15° southeast) is projected to the local coordinate system, and the component along the slope tangent direction (e.g., 5.8 mm / h) and the vertical direction component (e.g., 0.5 mm / h) are calculated as the slip driving force and deformation stretching, respectively, for subsequent risk assessment.

[0071] 303. Perform directional consistency screening on the slip driving force characterization quantity during the continuous monitoring period, retain the displacement components whose angle with the main slip direction is less than the threshold set by the dynamic feature library to participate in the cumulative calculation, and eliminate the abnormal deflection components caused by local collapse; In step 303, directional consistency screening refers to the process of filtering non-disaster related displacement components based on the main slip direction angle threshold.

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

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

[0074] 304. Based on the accumulated sliding period parameters corresponding to different lithology partitions in the dynamic feature library, adaptively accumulate the sliding driving force characterization quantity after screening in a time window to generate accumulated sliding energy values ​​in various azimuth regions; In step 304, the accumulated sliding energy value refers to the accumulated displacement value of the slope component within a specific time period, reflecting the accumulated sliding kinetic energy.

[0075] Time window adaptation refers to dynamically adjusting the length of the time window for cumulative calculation according to the lithological rheological characteristics.

[0076] In the embodiment of the present application, the cumulative slip period (such as a 30-minute window) is determined based on the lithological parameters in the dynamic feature library (such as a loose layer rheological coefficient of 0.3), and the sliding accumulation of the screened slip components is performed (such as a cumulative total of 58 mm in 4 windows) to generate the cumulative slip energy value of each area, thereby reflecting the evolution trend of the disaster in real time.

[0077] 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.

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

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

[0080] 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).

[0081] Here is a specific example: 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).

[0082] 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%.

[0083] 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.

[0084] In some embodiments, in step 304, based on the accumulated sliding period parameters corresponding to different lithological 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 regions, including: 401. Extracting rheological characteristic parameters of rock mass in the area where each grid unit is located according to the mechanical parameter table of rock type partition pre-stored in the dynamic characteristic library, wherein the rheological characteristic parameters of rock mass include viscoelastic coefficient and critical slip rate threshold, and determining the length of the initial slip cumulative cycle; In step 401, the rheological characteristic parameters of the rock mass refer to mechanical parameters reflecting the relationship between rock mass deformation and time, including viscoelastic coefficient and critical slip rate threshold.

[0085] The length of the initial slip accumulation period refers to the initial value of the time window for slip energy accumulation calculation determined according to the lithological rheological properties.

[0086] In the embodiment of the present application, the rheological characteristic parameters of the lithological partition where the grid unit is located (such as the viscoelastic coefficient of the loose layer of 0.3 and the critical slip rate of 5 mm / h) are extracted from the dynamic feature library, and the initial cumulative cycle length (such as 30 minutes for loose layer and 2 hours for bedrock) is determined by combining historical disaster data and expert rules, providing a time reference for the calculation of the slip displacement integral.

[0087] 402. Based on the initial slip cumulative cycle length, 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; In step 402, the sliding displacement integral value refers to the cumulative value of the displacement along the slope component in a single time window, reflecting the short-term accumulation of sliding kinetic energy.

[0088] The vertical slope component change trend refers to the change of the displacement rate in the vertical slope direction over time, which is used to monitor sudden deformation changes.

[0089] In the embodiment of the present application, based on the initial cumulative cycle length (such as 30 minutes), the screened slip driving force characterization quantity (such as 5.8 mm / h) is segmented and accumulated, and the slip displacement integral value in each time window is calculated (such as 8.7 mm accumulated in window 1), and the change trend of the vertical slope component is monitored (such as a sudden increase from 0.5 mm / h to 1.2 mm / h) to provide a basis for adjusting the cumulative cycle.

[0090] 403. 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; In step 403, the deformation mutation threshold refers to the critical value of the vertical slope component change rate. Exceeding this value indicates that the mountain has entered the accelerated deformation stage.

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

[0092] In an embodiment of the present application, when the rate of change of the vertical slope component exceeds the deformation mutation threshold set by the dynamic feature library (such as 0.8mm / h / 30min for a loose layer), the current accumulation period is shortened to a preset lower limit (such as 10 minutes), and the slip displacement integral value is recalculated (such as 3.2mm accumulated in a new window) to ensure that the slip energy accumulation value responds quickly to sudden deformation.

[0093] 404. Calculate the equivalent slip energy value in each time window according to the slip displacement integral value and the friction coefficient of the corresponding lithology partition, and superimpose the equivalent slip energy values ​​of consecutive time windows according to the exponential decay weight to generate a slip energy accumulation value.

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

[0095] The exponential decay weight refers to the weight coefficient used to superimpose the sliding energy values ​​of continuous time windows, reflecting the time decay effect of the sliding kinetic energy.

[0096] In the embodiment of the present application, the equivalent slip energy value (such as 3.48 J / m²) in a single time window is calculated based on the integral value of the slip displacement (such as 8.7 mm) and the friction coefficient of the lithology zoning (such as 0.4 for the loose layer), and the energy values ​​of the continuous windows are superimposed according to the exponential decay weights (such as 1.0 for window 1 and 0.8 for window 2) to generate the cumulative slip energy value (such as a cumulative 12.5 J / m²).

[0097] Here is a specific example: In a debris flow valley in the southwest, the dynamic feature library extracts the rock zoning parameters of a slope (loose layer viscoelastic coefficient 0.3, critical slip rate 5mm / h), and determines the initial accumulation period to be 30 minutes (step 401); the slip driving force characterization quantity (5.8mm / h) is accumulated in sections, and the slip displacement integral value of window 1 is calculated to be 8.7mm. At the same time, it is monitored that the vertical slope component suddenly increases to 1.2mm / h (step 402); because the change rate exceeds the threshold (0.8mm / h / 30min), the accumulation period is shortened to 10 minutes, and the integral value of window 2 is recalculated to be 3.2mm (step 403); the equivalent slip energy value is calculated based on the friction coefficient of 0.4 (window 1: 3.48J / m², window 2: 1.28J / m²), and the slip energy accumulation value of 12.5J / m² is generated by exponential decay weight superposition (step 404).

[0098] In summary, steps 401 to 404 dynamically adjust the accumulation period through the lithological rheological characteristics, so that the response speed of the accumulated slip energy value to sudden deformation is increased by 50%; the energy calculation model combining the friction coefficient and the exponential decay weight is significantly improved to significantly improve the accuracy of disaster risk assessment (error <5%), providing a reliable basis for the allocation of early warning signal coverage intensity, and reducing the overall underreporting rate by 30%.

[0099] In order to further improve the computational efficiency and trend prediction capability of the disaster early warning system, a dynamic adjustment mechanism of the cumulative period triggered by the energy saturation threshold is used, combined with the spatial distribution of slip energy and the propagation characteristics of rock rupture, to achieve visual prediction of disaster evolution trends and optimal resource allocation.

[0100] In some embodiments, step 304 further includes: 501. When the accumulated value of the sliding energy reaches the lithology partition energy saturation threshold value pre-stored in the dynamic feature library, the sliding accumulation period is extended to a preset upper limit value, and the calculation frequency of the sliding displacement integral value is reduced; In step 501, the energy saturation threshold refers to the preset critical value of the accumulated sliding energy of the lithology partition. Exceeding this value indicates that the disaster risk has entered a stable development stage.

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

[0102] In an embodiment of the present application, when the accumulated value of the slip energy reaches the energy saturation threshold of the rock type partition pre-stored in the dynamic feature library (such as 15 J / m² for the loose layer), the accumulation period is automatically extended to the preset upper limit value (such as 4 hours in the bedrock area), and the calculation frequency of the slip displacement integral value is reduced (such as from once every 10 minutes to once every 2 hours), while maintaining continuous monitoring of the vertical slope component to ensure efficient use of system resources.

[0103] 502. 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.

[0104] 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.

[0105] The spatiotemporal evolution trend chart is a visual prediction chart that combines time series and spatial distribution, showing the dynamic evolution trend of disaster risks.

[0106] In the embodiment of the present application, based on the spatial distribution data of the accumulated value of slip energy (such as the energy value of 20J / m² in the upstream area of ​​the valley and 10J / m² in the downstream area), a gradient field analysis algorithm is used to extract the direction of energy gradient change (such as southeast propagation), combined with the rock rupture propagation rate parameters pre-stored in the dynamic feature library (such as 0.5m / h in the loose layer), and a time-space interpolation is used to generate an evolution trend map for the next 6 hours (such as the red highlighted area extending 500m to the southeast), providing a decision-making basis for adjusting the coverage range of the early warning signal.

[0107] Here is a specific example: In a debris flow valley in the southwest, the accumulated sliding energy reaches the loose layer energy saturation threshold of 15 J / m² (step 501). The system extends the accumulation period from 30 minutes to 2 hours and reduces the calculation frequency to once per hour. Based on the energy spatial distribution data (20 J / m² upstream and 10 J / m² downstream), the energy gradient change direction is extracted as southeast (step 502). Combined with the rupture propagation rate of 0.5 m / h, a spatiotemporal evolution trend map is generated. It is predicted that the high-risk area will expand 300 m to the southeast in the next 3 hours, guiding the downstream villages to increase the warning signal strength to 35 dBm.

[0108] In summary, steps 501 to 502 reduce the system computing load by 40% through the periodic adjustment mechanism triggered by the energy saturation threshold; the spatiotemporal evolution trend diagram improves the prediction accuracy of the disaster diffusion path to 85%, increases the advance time of the early warning response by 1.5 hours, and optimizes the overall efficiency of emergency resource allocation.

[0109] In order to improve the accuracy of deformation monitoring and the reliability of risk assessment of flash flood disasters, the geological zoning data and radar scanning technology are integrated to dynamically optimize the radar beam parameters and accurately solve the displacement vector. The terrain risk heat map is generated by combining the critical value of rock zoning deformation to achieve spatial visualization of disaster risk.

[0110] In some embodiments, in step 102, 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 including deformation gradient and displacement direction is generated, including: 601. Based on the geological partition digital elevation model pre-stored in the dynamic feature library, determine the scanning coverage of each radar node in the radar array, dynamically adjust the radar beam incident angle 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 obtain a Doppler echo signal; In step 601, the geological zoning digital elevation model refers to a three-dimensional terrain data model divided according to geological features such as lithology and slope, and is used to determine radar scanning parameters.

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

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

[0113] 602. Perform phase demodulation on the Doppler echo signal received by each radar node, extract the phase change caused by surface micro-deformation, and calculate the displacement rate vector of each scanning point by combining radar wavelength and incident angle parameters; In step 602, phase demodulation refers to the process of extracting the phase variation caused by the ground surface deformation from the radar echo signal.

[0114] The displacement rate vector refers to a three-dimensional vector that contains the magnitude and direction of the displacement rate, reflecting the movement characteristics of the mountain.

[0115] In an embodiment of the present application, the Doppler echo signal received by the radar node is IQ orthogonal demodulated, the phase change Δφ (e.g., 0.8π radians) is extracted, and the radar wavelength λ (e.g., 3.9 mm) and the incident angle θ (e.g., 25°) are combined to calculate the displacement rate vector (e.g., 5.2 mm / h, 10° southeast), where T is the scanning period (e.g., 30 minutes).

[0116] 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; In step 603, the main sliding direction reference plane refers to the reference plane of the main sliding direction of the mountain determined based on historical sliding data and geological characteristics.

[0117] The slip component refers to the displacement component along the tangent direction of the base surface, reflecting the disaster-related sliding kinetic energy.

[0118] In an embodiment of the present application, the displacement rate vector is projected onto the main slip direction reference plane (such as a 30° inclination in the southeast direction) pre-stored in the dynamic feature library, and the component along the tangent direction of the reference plane (such as 4.8 mm / h) and the vertical direction component (such as 0.4 mm / h) are calculated by vector decomposition. The slip component is retained for deformation gradient calculation, and non-disaster-related interference of the vertical component is eliminated.

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

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

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

[0122] In the embodiment of the present application, according to the deformation critical values ​​of different lithology partitions in the dynamic feature library (such as 5mm / h for loose layers), the slip component is normalized (such as 4.8mm / h corresponds to a gradient value of 0.96), the deformation gradient field is generated through a spatial interpolation algorithm, and the displacement direction arrows are superimposed (such as 10° to the southeast), and finally the terrain risk heat map is output (the gradient value of the red highlighted area is ≥0.8).

[0123] Here is a specific example: In a debris flow valley 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 500m, and adjusts the beam incidence angle to 25° (step 601); after receiving the echo signal demodulation, the displacement rate vector 5.2mm / h (southeast by south 10°) is calculated (step 602); projected to the main slip reference plane (southeast 30°) to separate the slip component 4.8mm / h (step 603); based on the loose layer critical value of 5mm / h, a deformation gradient value of 0.96 is generated, and the area in the thermal map is highlighted in red (step 604).

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

[0125] In order to improve the accuracy of slip component extraction and disaster-related displacement identification, the main slip direction reference plane is constructed and a local projection coordinate system is established. The slip angle threshold of lithology zoning and spatial continuity detection are combined to achieve accurate screening and noise elimination of disaster-related slip components.

[0126] In some embodiments, in step 603, the displacement rate vector is projected onto the regional main slip direction reference plane pre-stored in the dynamic feature library, the slip component along the tangent direction of the reference plane and the deformation component perpendicular to the reference plane are separated, and the slip component is retained, including: 701. Based on the geological partition digital elevation model pre-stored in the dynamic feature library, extract the main slip direction reference plane of each grid unit, and establish a local projection coordinate system with the reference plane normal vector as the vertical axis; In step 701, the main sliding direction reference plane refers to the reference plane of the main sliding direction of the mountain determined according to the historical sliding trajectory and the characteristics of the geological division.

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

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

[0129] 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 main sliding direction reference line as a sliding consistency determination parameter; In step 702, the sliding consistency judgment parameter refers to the angle between the tangent direction component of the displacement vector and the main sliding direction reference line, which is used to screen the disaster-related displacement.

[0130] In an embodiment of the present application, the displacement rate vector (e.g., 6.3 mm / h, 15° southeast) is projected to the local coordinate system, the component along the tangent direction of the reference plane (e.g., 5.8 mm / h) and the vertical component (e.g., 0.5 mm / h) are calculated, and the angle (e.g., 15°) between the tangent direction component and the main slip direction reference line (30° southeast) is measured as a quantitative indicator of slip consistency.

[0131] 703. Based on the slip consistency judgment parameter and the lithology partition slip angle threshold 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 are selected as effective slip components; In step 703, the lithology partition slip angle threshold refers to the maximum allowable slip direction deviation set according to different lithologies, and is used to filter abnormal displacement.

[0132] In an embodiment of the present application, according to the slip angle threshold of the lithology zoning pre-stored in the dynamic feature library (such as 20° for loose layer and 15° for bedrock), the displacement vectors with angles less than the threshold in the tangent direction component are screened (such as 15° in the loose layer area meets the conditions), the excessive data (such as 45° to the north) are eliminated, and the effective slip components (such as 10° to the southeast) are retained.

[0133] 704. Perform spatial continuity detection on the effective slip component, eliminate isolated displacement signals caused by local collapse or rock rolling, and retain the slip component with the same displacement direction as the adjacent grid unit.

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

[0135] In the embodiment of the present application, a neighborhood analysis (such as a 3×3 grid) is performed on the effective slip component, the consistency index of the displacement direction of adjacent units is calculated (such as the southeast ±10° accounts for ≥70%), isolated signals (such as a single unit 45° north) are eliminated, and the continuously distributed slip component (such as southeast 8°) is retained to ensure the spatial continuity of the data.

[0136] Here is a specific example: In a debris flow valley in the southwest, the dynamic feature library extracts the main slip direction reference plane of the loose layer area (southeast 30°) and establishes a local projection coordinate system (step 701); the displacement vector (6.3mm / h, southeast 15°) is decomposed into a tangent component of 5.8mm / h (angle 15°) and a vertical component of 0.5mm / h (step 702); according to the loose layer slip angle threshold of 20°, the effective slip component is retained (step 703); through spatial continuity detection, the isolated north 45° signal is eliminated, and the southeast 8° continuous slip component is retained (step 704).

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

[0138] Figure 2 A structural diagram of an emergency communication device (or system) based on a multimodal warning broadcast terminal is provided for an embodiment of the present application, such as Figure 2 As shown, the device comprises: The disaster trigger monitoring module 21 is used to define the combination of soil moisture threshold and rainfall intensity threshold as a disaster trigger condition in a preset flash flood disaster sensitive area. When the real-time monitoring data meets the disaster trigger 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 a set of debris flow association rules; The deformation monitoring and analysis module 22 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; The emergency communication networking module 23 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 24 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 are 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.

[0139] Figure 2 The emergency communication device based on the multi-modal early warning broadcast terminal can perform Figure 1 The implementation principle and technical effect of the emergency communication method based on a multimodal warning broadcast terminal described in the illustrated embodiment will not be repeated. The specific manner in which each module and unit performs operations in the emergency communication device based on a multimodal warning broadcast terminal in the above embodiment has been described in detail in the embodiment of the method, and will not be elaborated here.

[0140] In one possible design, Figure 2An emergency communication device based on a multimodal warning broadcast terminal of the embodiment shown can be implemented as a computing device, such as Figure 3 As shown, the computing device may include a storage component 31 and a processing component 32; 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 .

[0141] The processing component 32 is used for the above Figure 1 The embodiment provides an emergency communication method based on a multimodal early warning broadcast terminal.

[0142] The processing component 32 may 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 may 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 to perform the above method.

[0143] The storage component 31 is configured to store various types of data to support operations at 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.

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

[0145] The input / output interface provides an interface between the processing component and the peripheral interface module, which may be an output device, an input device, etc.

[0146] The communication component is configured to facilitate, among other things, wired or wireless communications between the computing device and other devices.

[0147] Among them, the computing device can be a physical device or an elastic computing host provided by a cloud computing platform, etc. In this case, the computing device can refer to a cloud server, and the above-mentioned processing components, storage components, etc. can be basic server resources rented or purchased from the cloud computing platform.

[0148] The present application also provides a computer storage medium storing a computer program, wherein the computer program can achieve the above-mentioned Figure 1The illustrated embodiment is an emergency communication method based on a multimodal early warning broadcast terminal.

[0149] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the systems, devices and units described above can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0150] The device embodiments described above are merely illustrative, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. Ordinary technicians in this field can understand and implement it without paying creative labor.

[0151] Through the description of the above implementation methods, those skilled in the art can clearly understand that each implementation method can be implemented by means of software plus a necessary general hardware platform, and of course, can also be implemented by hardware. Based on this understanding, the above technical solution is essentially or the part that contributes to the prior art can be embodied in the form of a software product, and the computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a disk, an optical disk, etc., including a number of instructions for a computer device (which can 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.

[0152] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present application, rather than to limit it. Although the present application has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the 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