A landslide-debris flow disaster chain evaluation method based on a graph neural network

Abstract

The application discloses a landslide-mudslide disaster chain evaluation method based on a graph neural network, generates a composite ground object unit through a composite curvature field watershed segmentation and hydrological analysis, and abstracts the composite ground object unit into a geographic node; multi-source environmental factors are optimized through mutual information screening; a directed weighted graph is constructed in combination with spatial similarity; an improved GraphSAGE model is adopted, neighbor weight, self-loop and residual information are introduced, and random sampling is cancelled; a sample set suitable for a small data area is constructed through a unit splitting-feature matching-terrain fitting mechanism; and model training and risk prediction are completed. The method solves the problems of single disaster evaluation, strong sample dependence and high omission rate in the prior art, has excellent performance indexes, effectively reduces the omission risk of a high-risk area, and provides precise and efficient technical support for disaster chain risk evaluation of a small watershed in a mountainous area with less data.

Description

TECHNICAL FIELD

[0001] The present application belongs to the technical field of geological disaster detection, and particularly relates to a landslide-debris flow disaster chain evaluation method based on a graph neural network. BACKGROUND

[0002] Landslides and debris flows are common geological disasters in mountainous areas, and there is a close correlation between the two in the process of occurrence. In some areas, loose deposits generated after landslides often become an important source of debris flow, and under the action of environmental factors such as precipitation and topography, these loose sources are easily transported by water flow, thereby inducing debris flow and forming a landslide-debris flow disaster chain. This type of disaster chain will have an impact on the ecological environment, infrastructure and personnel safety in the region, so it is of practical significance to carry out risk assessment on landslide-debris flow disaster chains.

[0003] At present, there has been some research on risk prediction technology for landslides and debris flows, but there are still many deficiencies. Firstly, existing technologies are mostly focused on risk evaluation of a single disaster, and there is no systematic quantitative evaluation framework for landslide-debris flow disaster chains, making it difficult to fully depict the causal transmission relationship between landslides and debris flows. Secondly, traditional evaluation methods rely too much on high-quality training samples, and in areas where the number of samples is insufficient or the quality is not high, the model prediction performance may decrease. Thirdly, the spatial feature expression ability of most methods is limited, and they cannot effectively capture the potential relevance between similar geographical environments in distant areas, making it difficult for the evaluation results to cover all risk characteristics in the region. Fourthly, existing methods are often limited by the computational complexity of modeling the entire graph, and it is difficult to achieve efficient operation when carrying out risk evaluation in large-scale regions (such as cross-city, cross-province basins), limiting the popularization and application of the technology in actual scenarios. SUMMARY

[0004] In order to solve the problems existing in the prior art, the present application provides a landslide-debris flow disaster chain evaluation method based on a graph neural network, which aims to provide a systematic technical path for disaster chain risk evaluation by constructing a risk evaluation system for landslide-debris flow disaster chains.

[0005] The technical solution adopted by the present application is as follows:

[0006] In a first aspect, the present application provides a landslide-debris flow disaster chain evaluation method based on a graph neural network, comprising the following steps:

[0007] Step 100: generating a composite ground feature unit by composite curvature field watershed segmentation and hydrological analysis, and abstracting the composite ground feature unit as a geographic node;

[0008] Step 200: constructing node features corresponding to the geographic nodes, and verifying the independence of each environmental factor in the node features using a mutual information method;

[0009] Step 300: a directed weighted graph is constructed according to the environmental factors and spatial similarity, and the directed weighted graph is used to reflect the spatial transmission relationship of the landslide-debris flow disaster chain;

[0010] Step 400: an improved GraphSAGE model is designed, and information in the directed weighted graph is transmitted and node representation is optimized by the improved GraphSAGE model;

[0011] Step 500: the improved GraphSAGE model is trained and verified, and risk prediction and spatial mapping of the landslide-debris flow disaster chain are completed.

[0012] With reference to the first aspect, in a first implementation manner of the first aspect, the process of generating the composite ground object unit in the step 100 comprises: first obtaining a slope unit through composite curvature field watershed segmentation, and then introducing a gully buffer zone on the basis of the slope unit, and further forming the composite ground object unit.

[0013] With reference to the first implementation manner of the first aspect, in a second implementation manner of the first aspect, after the composite ground object unit is formed in the step 100, the operation of quantifying the source area features and flow-through area features corresponding to the composite ground object unit is further included.

[0014] With reference to the first aspect, in a third implementation manner of the first aspect, the process of constructing the node features in the step 200 comprises: collecting multi-source heterogeneous environmental factors, and taking the multi-source heterogeneous environmental factors as initial input factors of the node features; and the step of verifying the independence of each environmental factor by using the mutual information method comprises: calculating the correlation degree between each initial input factor by using a normalized mutual information method.

[0015] With reference to the third implementation manner of the first aspect, in a fourth implementation manner of the first aspect, after the correlation degree between each initial input factor is calculated in the step 200, the operation of removing a strongly dependent factor according to the correlation degree is further included, and the strongly dependent factor is an initial input factor with a correlation degree exceeding a preset threshold.

[0016] With reference to the first aspect, in a fifth implementation manner of the first aspect, the process of designing the improved GraphSAGE model in the step 400 comprises: introducing a neighbor weight mechanism, retaining self-loop information, retaining residual information, and canceling a random sampling operation for the GraphSAGE model.

[0017] In combination with the fifth implementation manner of the first aspect, the application provides a sixth implementation manner of the first aspect, and the neighbor weight mechanism is specifically introducing that different weights are assigned to different neighbor nodes according to the influence of each geographic node on the geographic node.

[0018] In combination with the fifth implementation manner of the first aspect, the application provides a seventh implementation manner of the first aspect, and the self-loop information is specifically reserving the feature information of each geographic node itself when the improved GraphSAGE model performs information transmission, and the self-loop information of the geographic node is not excluded; and the residual error information is specifically connecting the node representation result of the previous round with the node representation result of the current round in residual error when the improved GraphSAGE model performs node representation optimization, and the historical representation information is reserved.

[0019] In combination with the first aspect, the application provides an eighth implementation manner of the first aspect, and the step 400 further includes an operation of constructing a disaster chain sample set before the improved GraphSAGE model performs information transmission and node representation optimization, and the disaster chain sample set is constructed by using a hierarchical sampling strategy to count positive and negative samples according to environment units to ensure balanced sample distribution.

[0020] In combination with the eighth implementation manner of the first aspect, the application provides a ninth implementation manner of the first aspect, and the process of constructing the disaster chain sample set further includes an operation of splitting disaster chain labels, and the data diversity of the disaster chain sample set is enhanced by the label splitting.

[0021] The application has the following beneficial effects:

[0022] (1) The application constructs a risk assessment system for a landslide-debris flow disaster chain, abstracts disasters and environmental factors into a graph structure, realizes modeling and causal relationship quantization of the landslide-debris flow disaster chain, makes up for the limitation of single disaster research in the prior art, and provides a systematic technical path for disaster chain risk assessment.

[0023] (2) The application uses an improved GraphSAGE model and a transfer learning mechanism, reduces the dependence on high-quality training samples, maintains good model performance in a sample-scarce area, improves the generalization ability of the model in a regional expansion scenario, and helps to expand the application scope of the technology.

[0024] (3) The application enhances the spatial feature expression ability by constructing a directed weighted graph according to environmental factors and spatial similarity, effectively captures the potential correlation between distant and geographically similar areas, and makes the risk assessment result more comprehensive and more in line with the actual disaster distribution characteristics.

[0025] (4) The GraphSAGE model is improved, neighbor weight mechanism is introduced, self-loop and residual information are reserved, and random sampling is cancelled, so that the information transmission and node representation process are optimized, the calculation complexity of whole graph modeling is reduced, efficient risk evaluation in large-scale region is facilitated, and the practical application value of the technology is improved;

[0026] (5) The present application optimizes the redundancy of multi-source heterogeneous environmental factors by mutual information method, quantifies the independence of factors and eliminates strong dependent factors, improves the robustness of input data; at the same time, through the hierarchical sampling strategy, the disaster chain sample set is constructed and the label splitting is supported, the sample distribution is balanced, the data diversity is enhanced, a more reliable data foundation is provided for model training, and the stability of risk prediction result is further improved. BRIEF DESCRIPTION OF DRAWINGS

[0027] Figure 1 is the schematic diagram of information transmission of the improved GraphSAGE model in the embodiment of the present application;

[0028] Figure 2 is the schematic diagram of node representation optimization of the improved GraphSAGE model in the embodiment of the present application;

[0029] Figure 3 is the schematic diagram of risk prediction mapping in the embodiment of the present application;

[0030] Figure 4 is the schematic diagram of risk evaluation of the target region in the embodiment by the random forest method in the prior art;

[0031] Figure 5 is the schematic diagram of risk evaluation of the target region in the embodiment by the support vector machine method in the prior art;

[0032] Figure 6 is the schematic diagram of risk evaluation of the target region in the embodiment by the multilayer perception machine method in the prior art;

[0033] Figure 7 is the schematic diagram of risk evaluation of the target region in the embodiment by the method in the embodiment of the present application. DETAILED DESCRIPTION

[0034] The present application will be further explained in combination with the drawings and specific embodiments.

[0035] To make the purposes, technical solutions, and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be described clearly and completely below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only some of the embodiments of the present application, rather than all the embodiments. The components of the embodiments of the present application described and shown in the drawings can be arranged and designed in various different configurations.

[0036] Therefore, the detailed description of the embodiments of the present application provided below in the drawings is not intended to limit the scope of the claimed present application, but only represents selected embodiments of the present application. All other embodiments obtained by a person of ordinary skill in the art without making creative efforts based on the embodiments in the present application are within the scope of protection of the present application. Embodiments

[0037] The scenario to be solved by the present embodiment is risk assessment of landslide-debris flow disaster chain in mountainous small watershed. The terrain in such areas is undulating, and gullies are densely distributed. Landslides are easy to induce debris flow due to the supply of loose deposits to gullies, but the existing assessment steps can be further refined in terms of adapting to the topographic characteristics, environmental factor specificity, and model information transmission of small watersheds, so as to better fit the disaster chain occurrence regularity of small watersheds.

[0038] The present embodiment discloses a landslide-debris flow disaster chain assessment method based on a graph neural network, comprising the following steps:

[0039] Step 100: For the topographic characteristics of mountainous small watersheds, first, the compound curvature field is used to divide the watershed to obtain slope units that adapt to the slope gradient and aspect characteristics of the region, and then, the gully buffer zone covering the main confluence area is introduced based on the slope units through hydrological analysis in combination with the gully strike and confluence path of the small watershed, to generate compound ground units that can completely contain the potential landslide occurrence area and potential debris flow flow-through area in the small watershed, and then each compound ground unit is abstracted as a geographic node with geographic spatial attributes;

[0040] Step 200: Collect multi-source heterogeneous environmental factors related to landslide-debris flow disaster chain in mountainous small watersheds, use these environmental factors as initial input of node features corresponding to the geographic nodes, and then use the mutual information method to calculate the correlation degree between each initial input factor, verify the independence of different environmental factors in the small watershed scenario, and select factors that are more targeted in affecting the disaster chain in the region;

[0041] Step 300: According to the environmental factors screened in step 200, the disaster influence relationship between each geographical node in the mountainous small watershed is analyzed due to factors such as terrain and material source, and the weight of the influence relationship between each geographical node is determined by combining the geographical environment similarity of different regions in the small watershed, and a directed weighted graph reflecting the spatial transmission relationship of the landslide-mudslide disaster chain in the small watershed is constructed.

[0042] Step 400: According to the characteristics of the number of geographical nodes and the correlation density between nodes in the mountainous small watershed, an improved GraphSAGE model is designed, a neighbor weight mechanism based on the disaster influence degree of the small watershed is introduced in the model, the feature information and historical representation information of each geographical node are retained, and random sampling is cancelled to ensure the information capture of the key associated nodes in the small watershed. The improved GraphSAGE model is used to efficiently transmit the information in the directed weighted graph and optimize the node representation.

[0043] Step 500: Collect historical landslide-mudslide disaster chain data in the mountainous small watershed, construct a disaster chain sample set suitable for this region, train the improved GraphSAGE model using the sample set, adjust the model parameters during the training process in combination with the disaster occurrence regularity of the small watershed, and then test the prediction effect of the model through independent verification data. After the test is passed, the model is used to complete the risk prediction and spatial mapping of the landslide-mudslide disaster chain in the mountainous small watershed.

[0044] Further, the embodiment is aimed at the characteristics of the small watershed, i.e. the narrow terrain, the significant spatial difference of the slope gradient, the alternating distribution of gentle slope section and steep slope section, the easy occurrence of landslide in the steep slope section, the easy accumulation of loose material in the gentle slope section, etc. The specific operations of obtaining the slope unit through the compound curvature field watershed segmentation, introducing the channel buffer zone to form the compound ground object unit, and quantifying the characteristics of the material source area and the flow area are optimized, and the scheme content is as follows:

[0045] Obtaining the slope unit through the compound curvature field watershed segmentation

[0046] 1. Composition of compound curvature field

[0047] The compound curvature field is a terrain feature field that combines plane curvature and profile curvature, which is used to accurately identify the boundary contour and internal uniformity of the slope in the small watershed.

[0048] Plane curvature: reflects the bending shape of the slope in the small watershed in the horizontal direction. The positive value corresponds to the ridge line, which is the natural boundary of the slope unit and can block the transport of material sources across the unit. The negative value corresponds to the valley line, which is the transition zone between the slope unit and the channel, and is the key path for the material source to enter the channel.

[0049] Profile curvature: reflects the bending shape of the small watershed slope in the vertical slope direction, negative value corresponds to convex slope, i.e. the transition zone from gentle slope to steep slope in the small watershed, loose material is easy to slide, which is the potential starting area of landslide, positive value corresponds to concave slope, i.e. the bottom of gentle slope or steep slope in the small watershed, which is easy to accumulate loose material, which is the source reserve area;

[0050] Combined with the slope characteristics of the small watershed, the composite curvature field is calculated:

[0051] According to the slope classification of the small watershed, the slope of the gentle slope section is less than 15°, the slope of the medium slope section is 15° to 35°, and the slope of the steep slope section is greater than 35°. Different slope sections are assigned different weights for plane curvature and profile curvature.

[0052] The profile curvature weight of the medium slope section and the steep slope section is increased to 0.5, i.e. the plane curvature weight is 0.5;

[0053] The plane curvature weight of the gentle slope section (source accumulation area) is increased to 0.6, i.e. the profile curvature weight is 0.4, to ensure that the segmentation logic is preferentially adapted to the topographic characteristics of the potential landslide area.

[0054] 2. Based on the composite curvature field, the steps of watershed segmentation to obtain the slope unit

[0055] (1) Small watershed DEM data preprocessing

[0056] Get the small watershed high-precision digital elevation model data, first fill the pit, i.e. eliminate local concave areas in DEM to avoid distortion in subsequent slope and curvature calculation, then extract the global slope data of the small watershed based on DEM, and complete the slope classification according to gentle slope section, medium slope section and steep slope section to form the small watershed slope classification map.

[0057] (2) Composite curvature field calculation

[0058] Based on the preprocessed DEM and slope classification map, use terrain analysis tools to calculate the plane curvature and profile curvature of each grid cell:

[0059] For medium slope section and steep slope section grid: composite curvature value = 0.5 x absolute value of plane curvature + 0.5 x absolute value of profile curvature;

[0060] For gentle slope section grid: composite curvature value = 0.6 x absolute value of plane curvature + 0.4 x absolute value of profile curvature;

[0061] Through the above calculation, the small watershed composite curvature field grid map is obtained, the high curvature value area in the grid map with curvature value greater than 2.0 corresponds to the slope boundary, and the low curvature value area with curvature value less than 0.8 corresponds to the uniform area inside the slope.

[0062] (3) Watershed segmentation and slope unit generation

[0063] First step: Generate the composite curvature field gradient map, that is, calculate the gradient of the composite curvature field through spatial analysis tools, and the high gradient value area corresponds to the slope boundary, and the low gradient value area corresponds to the slope interior;

[0064] Second step: Screen the segmentation seed points, and only in the low gradient value area (composite curvature value <0.8) of the slope section and steep slope section in the small watershed to screen the seed points, exclude the area without landslide potential in the gentle slope section, and ensure that the seed points correspond to the core area of the landslide potential area;

[0065] Third step: Perform watershed segmentation, take the gradient map as the segmentation boundary basis, and take the seed point as the unit core basis, complete the segmentation through the terrain analysis tool to obtain the initial slope unit raster map;

[0066] Fourth step: Unit post-processing, merge the initial slope units with an area <0.02 km² in the small watershed, avoid excessive unit calculation redundancy, delete the initial slope units with a slope <10°, and finally obtain the slope unit vector map adapted to the slope gradient characteristics of the small watershed.

[0067] Definition, data acquisition and generation logic of the gully buffer zone

[0068] 1. Definition of the gully buffer zone

[0069] The gully buffer zone in the small watershed refers to the area extending a certain width to the slope on both sides with the natural gully (including seasonal gully) in the small watershed as the center. This area is the transition zone of landslide material sources in the slope unit entering the gully, and is also the peripheral area of the core channel of the debris flow. Its width needs to be dynamically adjusted in combination with the slope gradient of the small watershed to avoid the lack of regional adaptability caused by fixed width.

[0070] 2. Data acquisition of the gully buffer zone

[0071] Basic data sources: small watershed 1:10,000 DEM data, high-resolution remote sensing images, and small watershed field hydrological investigation data, which contain gully direction, perennial flow / seasonal flow attributes, etc.;

[0072] Gully extraction: Based on DEM data, calculate the flow accumulation through hydrological analysis tools, select the grid with flow accumulation >2000 as the initial river network, and then superimpose high-resolution remote sensing images and field investigation data to modify the initial river network, that is, delete artificial drainage channels and supplement seasonal gullies not recognized by DEM, and then obtain the actual gully vector data of the small watershed;

[0073] Buffer width determination: Combined with the slope gradient classification of small watershed, the width of the buffer zone is set. The width of the buffer zone corresponding to the steep slope section with a gradient greater than 35° is 10-20 m. The source transport distance of the steep slope section is short, and the buffer zone does not need to be too wide. The width of the buffer zone corresponding to the medium slope section with a gradient ranging from 15° to 35° is 20-30 m. The width of the buffer zone corresponding to the gentle slope section with a gradient less than 15° is 30-40 m. The source transport distance of the gentle slope section is long, and the coverage range of the buffer zone needs to be expanded.

[0074] 3. Generation of channel buffer zone

[0075] Based on the corrected small watershed channel vector data, the channel buffer zone vector diagram is generated by the buffer analysis function of the spatial analysis tool according to the dynamic width, and then the small watershed boundary vector is superimposed to cut off the buffer zone part exceeding the small watershed range, thereby obtaining the final small watershed channel buffer zone data.

[0076] III. Quantitative description of the characteristics of the source area and the flow area

[0077] 1. Quantification of source area characteristics

[0078] Taking the slope unit as the basic unit and combining the slope gradient characteristics of the small watershed, the source amount potential and the source stability are quantified from two dimensions:

[0079] Source amount potential: The slope unit is assigned by slope classification, the steep slope section is assigned 3, the medium slope section is assigned 2, and the gentle slope section is assigned 1. Then, the relative thickness of the loose deposits in the unit is combined, that is, the vegetation coverage degree is inversely calculated based on the remote sensing image: the unit with vegetation coverage degree less than 30% is assigned 3, the unit with vegetation coverage degree ranging from 30% to 60% is assigned 2, and the unit with vegetation coverage degree greater than 60% is assigned 1. The product of the two is taken as the source amount potential index. The higher the index, the more landslide sources that can be provided in the unit.

[0080] Source stability: The slope unit is assigned by profile curvature type, the convex slope is assigned 3, the straight slope is assigned 2, and the concave slope is assigned 1. The more negative the curvature value, the worse the source stability, and the more prone to landslides and source replenishment.

[0081] 2. Quantification of flow area characteristics

[0082] Taking the channel buffer zone as the basic unit, the source transport capacity and the flow channel smoothness are quantified from two dimensions:

[0083] Source transport capacity: The channel buffer zone is classified by the corresponding flow accumulation. The flow accumulation greater than 5000 is assigned 3, the flow accumulation ranging from 3000 to 5000 is assigned 2, and the flow accumulation less than 3000 is assigned 1. The higher the flow accumulation, the stronger the water flow transport capacity.

[0084] Flowing channel smoothness: the channel buffer zone corresponding to the channel longitudinal slope gradient classification, namely longitudinal slope gradient > 5 ‰, 2 ‰-5 ‰, <2 ‰, the greater the longitudinal slope gradient, the steeper the channel slope, the more smooth the debris flow, and it is not easy to accumulate.

[0085] Four, the final formation of the composite land unit

[0086] The slope unit of the above-obtained adaptive small watershed slope gradient is spatially superimposed and analyzed with the channel buffer zone of dynamic width, the area where the slope unit and the channel buffer zone overlap and the overlapping area accounts for ≥30% of the slope unit area is retained, the slope unit and the channel buffer zone in the overlapping area are merged to form a composite land unit containing the landslide potential source area (slope unit) and the debris flow passing area (channel buffer zone), each composite land unit is attached with the above-mentioned quantified source area characteristic index and the flowing area characteristic index.

[0087] Further, for the node characteristics constructed in the above step 200 and the verification of the independence of environmental factors, optimization deepening is carried out, through specific mutual information calculation, independence quantization and redundancy elimination logic, the problems of low model training efficiency and prediction accuracy interference caused by multi-source factor redundancy are solved, and the composite land unit is connected.

[0088] The embodiment further clarifies the whole process of screening and optimization in view of the characteristics of small watershed multi-source environmental factor type, specific steps are as follows:

[0089] Step 201: Around the formation mechanism of landslide-debris flow disaster chain in small watershed, collect multi-source heterogeneous environmental factors corresponding to the composite land unit, and classify according to the influence dimension:

[0090] Topographic factors: average slope of composite land unit, relative height difference, profile curvature, plane curvature;

[0091] Geological factors: rock-soil type in the unit, fault structure density, weathering layer thickness grade;

[0092] Hydrological factors: multi-year average precipitation in the unit, runoff accumulation intensity, groundwater depth grade;

[0093] Vegetation factors: vegetation coverage in the unit, vegetation type (tree / shrub / herb);

[0094] Human activity factors: road density in the unit, engineering activity disturbance range proportion.

[0095] Factor pretreatment: unify different dimension factors to the same numerical interval to eliminate the interference of dimension difference on mutual information calculation;

[0096] Outlier treatment: remove extreme values caused by data acquisition errors, and keep effective data that conforms to the terrain and geological regularities according to the distribution characteristics of the same factor in small watersheds;

[0097] Factor gridding: convert all environmental factors into raster data with the same spatial resolution as the composite ground feature units, ensuring that each geographic node can accurately match the corresponding factor value.

[0098] Step 202: Calculate the correlation degree between any two environmental factors using normalized mutual information. Normalized mutual information quantifies the correlation degree as a value between 0 and 1, avoiding the distortion of original mutual information caused by different factor value ranges.

[0099] Mutual information matrix construction: take all collected environmental factors as matrix rows and columns, calculate the normalized mutual information value between each pair of factors, such as average slope and relative height difference, vegetation coverage and weathering layer thickness grade, and form a factor-factor mutual information matrix.

[0100] Matrix element definition: let the matrix be M n×n , n is the total number of factors, and element M ij represents the normalized mutual information value between the i th factor and the j th factor. M ij =M ji , the diagonal element M ii =1 .

[0101] Step 203: Calculate the comprehensive independence index for each environmental factor to quantify its independence degree with all other factors. The formula logic is as follows:

[0102] Average correlation degree of a single factor: calculate the mean value of all non-diagonal elements in the corresponding mutual information matrix row (or column) of the factor, denoted as NMI i , NMI i The larger the value, the stronger the overall correlation between the factor and other factors, and the weaker the independence.

[0103] Comprehensive independence index: I i =1-NMI i

[0104] I i The value range is 0-1, I iThe closer to 1, the stronger the comprehensive independence of the factor; the closer to 0, the higher the redundancy of the factor.

[0105] Independence classification: according to the characteristics of the small watershed disaster chain influencing factor, the comprehensive independence index I i is divided into three levels:

[0106] High independence: I i ≥0.7, the factor is weakly associated with other factors, and the information is unique;

[0107] Medium independence: 0.4≤ I i <0.7, the factor has weak association with some factors, but still has effective information);

[0108] Low independence: I i <0.4, the factor has strong association with other multiple factors, and the redundant information accounts for a high proportion.

[0109] Step 204: Combined with mutual information matrix and comprehensive independence index, double determination of strong dependent factors:

[0110] Rule 1: If the normalized mutual information value between two factors M ij ≥ preset adaptive threshold, which is calibrated based on small watershed historical disaster data, to ensure that truly redundant factors are removed, then the two are determined as a strong dependent factor pair;

[0111] Rule 2: In the strong dependent factor pair, the factor with I i Smaller comprehensive independence index is preferentially removed, and the factor with stronger independence and more unique information is retained;

[0112] Rule 3: If a single factor has M ij ≥ preset adaptive threshold, regardless of its I i Size, it is directly determined as a strong dependent factor and removed, and the factor has too much redundant information, which is limited to model gain.

[0113] Iterative removal process:

[0114] First round of removal: identify and remove strong dependent factors according to the above rules to obtain a preliminary screened factor set;

[0115] Second round of verification: rebuild the mutual information matrix for the preliminary screened factor set, calculate the comprehensive independence index of each factor, and if there are still strong dependent factor pairs, repeat the removal operation;

[0116] Termination condition: in the screened factor set, any two factors M ij < Preset adaptation threshold, and the comprehensive independence index of all factors I i ≥ 0.4, stop iteration.

[0117] Step 205: The adaptability verification of the screened factor set will finally screen the optimized factor set, and spatially match the composite ground object unit (geographic node) generated in step 100 to ensure that each geographic node can correspond to complete optimized factor data; at the same time, combined with the key influence mechanism of the small watershed landslide-mudflow disaster chain, it is verified whether the screened factor set covers the core dimension of the disaster chain formation, so as to avoid the loss of key information due to excessive elimination.

[0118] Through the above steps, a node feature set without redundancy, high independence and strong pertinence is finally constructed for each geographic node.

[0119] Further, for step 300, the embodiment provides further preferred limitations, and the specific optimization is as follows:

[0120] Step 301: Based on the function attributes of the source area slope unit and the flow-through area channel buffer of the composite ground object unit, combined with the conduction path of the landslide-mudflow disaster chain, the directed association between geographic nodes is determined to avoid meaningless spatial association.

[0121] The association direction rule is that the source area node points to its spatially adjacent flow-through area node to reflect the conduction direction of source supply to debris flow flow-through; the flow-through area node points to the downstream flow-through area node to reflect the conduction direction of flow-through continuation; the flow-through area node is not allowed to point to the source area node, and the downstream flow-through area node is not allowed to point to the upstream flow-through area node, which conforms to the irreversible conduction characteristics of the disaster chain.

[0122] The spatial adjacency determination is based on the vector boundary of the composite ground object unit. If the boundaries of two geographic nodes overlap or the distance does not exceed the preset threshold value based on the small watershed channel density calibration, such as 500 meters, it is determined to be spatially adjacent and has the basis for association.

[0123] Step 302: Based on the source area factor subset, the flow-through area factor subset and the common factor subset screened in step 200, the disaster influence intensity of each node on the associated node is quantified, and the weight size directly reflects the contribution degree of the factor to the conduction of the disaster chain.

[0124] The factor weight distribution follows a professional priority orientation. The influence weight of the material source area node is based on the distribution proportion of loose deposits, average slope, and rock-soil shear strength as core factors, with weight proportions of 30%, 25%, and 20% respectively. Vegetation coverage and fault structure density are auxiliary factors, with weight proportions of 15% and 10% respectively. The higher the factor value, such as a high proportion of loose deposits and steep slope, the greater the influence weight of the material source area node on the downstream flow area node. The influence weight of the flow area node is based on the channel longitudinal gradient and the convergence cumulative intensity as core factors, with weight proportions of 35% and 30% respectively. Channel width and rainstorm intensity grade are auxiliary factors, with weight proportions of 20% and 15% respectively. The higher the factor value, such as a large longitudinal gradient and strong convergence, the greater the influence weight of the flow area node on the downstream flow area node.

[0125] In the calculation of the comprehensive influence weight of the node, each factor value in the optimized factor set of each geographic node is standardized in the range of 0-1 to eliminate dimensional differences. Then the weighted sum is calculated according to the above weight distribution to obtain the comprehensive influence weight of the node, which is in the range of 0-1 and serves as a core component of the subsequent edge weight.

[0126] Step 303: Spatial similarity quantification focuses on the similarity of disaster chain transmission potential rather than simply geographical distance. It combines the environmental attributes of composite ground units with the similarity calculation logic defined by disaster chain links. The similarity evaluation index selects common factors strongly related to disaster chains, namely fault structure density and multi-year average precipitation, as well as exclusive factors for corresponding functional areas. For material source area nodes, rock-soil types are added, and for flow area nodes, riverbed material composition is added.

[0127] The similarity calculation method is to calculate the Euclidean distance of two associated nodes, such as material source area A and flow area B, and flow area B and flow area C, in each similarity index. Then, according to the importance of the index to the disaster chain, the weight is distributed, with a structure density of 30%, a precipitation of 30%, and an exclusive factor of 40%. The weighted sum is the index difference. The spatial similarity coefficient is 1 minus the normalized index difference, with a value range of 0-1. The closer the coefficient is to 1, the more similar the disaster chain influence conditions of the two nodes, and the stronger the correlation of the transmission potential.

[0128] At the same time, the spatial similarity coefficient is distance decay corrected. If the distance between the associated nodes is less than 500 meters, the correction coefficient is 1, i.e. no attenuation. If the distance exceeds 500 meters, the correction coefficient decreases with the distance gradient, such as 0.8 for 500-1000 meters and 0.6 for 1000-1500 meters, ensuring that the influence of spatial distance on transmission efficiency conforms to reality, i.e. the transmission potential decreases at a long distance.

[0129] Step 304: Based on the results of the above steps, a complete directed weighted graph is constructed. Each element of the graph corresponds to the professional logic of the disaster chain.

[0130] The node definition abstracts the geographical node after step 100 as the vertex of the graph, and each vertex is attached with node type, source area or flow area, comprehensive influence weight, optimization factor set data, spatial coordinates and other attributes, and the function positioning of the node in the disaster chain is clear.

[0131] The direction of the edge follows the association direction rule of step 301, and only the directed edges from the source area to the flow area and from the upstream flow area to the downstream flow area are drawn, and the arrow direction of the edge directly corresponds to the conduction path of the disaster chain.

[0132] The edge weight is the product of the node comprehensive influence weight and the modified spatial similarity coefficient, and the calculation formula is edge weight W ij The geometric mean of the comprehensive influence weight of node i and the comprehensive influence weight of node j is multiplied by the modified spatial similarity coefficient of node i and j. The geometric mean of the comprehensive influence weight of two nodes is taken to ensure the balance of two-way influence, and the similarity coefficient is multiplied to reflect the enhancement effect of similar environmental conditions on the conduction intensity.

[0133] When the graph structure is simplified, the directed edges with edge weight less than 0.1 are removed, and the core associated edges with weight not less than 0.1 are retained to avoid redundancy of the graph structure.

[0134] Step 305: Verify whether the graph structure fits the actual spatial association of the landslide-mudflow disaster chain to ensure no logical deviation. When verifying the conduction path, superimpose the historical landslide-mudflow disaster path graph of the small watershed to check whether the directed edges from the source area to the flow area and then to the downstream flow area in the graph are consistent with the actual disaster spread path. If there is deviation, such as a historical disaster path not reflected in the graph, adjust the spatial threshold or similarity index of node association.

[0135] When verifying the weight rationality, compare the edge weight distribution of high-risk and low-risk disaster areas to ensure that the node associated edge weight in the high-risk disaster area is significantly higher than that in the low-risk area, such as the average edge weight in the high-risk area not less than 0.4 and the low-risk area not more than 0.2, which reflects the positive correlation between weight and disaster risk. When verifying the function area association, check whether the source area node is associated with at least one flow area node to ensure that the source can flow into the flow area, and whether the flow area nodes form a continuous association chain from upstream to downstream to ensure that the debris flow can flow completely, avoiding isolated nodes without associated edges. If there are isolated nodes, re-determine their spatial adjacency.

[0136] Through the above refinement steps, the directed and weighted graph constructed not only has the characteristics of mathematical graph structure, but also deeply integrates the functional attributes of complex ground object units and the conduction rules of landslide-mudflow disaster chain, providing a structured input that fits the professional scenario for improving the information transmission of the GraphSAGE model in step 400.

[0137] Further, with reference to Figure 1 and Figure 2 the following preferred solutions are proposed for step 400.

[0138] Step 401: Explicitly improve the input adaptation logic of the GraphSAGE model, taking the directed and weighted graph constructed in step 300 as the core input carrier, and taking the optimized factor set of the geographical node as the initial feature vector of the model. At the same time, import the node type source area or flow area attribute, node comprehensive influence weight and edge weight data, ensure the deep matching of the model input and the spatial correlation characteristics of the landslide-debris flow disaster chain, the node function attribute, and provide data basis for subsequent targeted improvement.

[0139] Step 402: Cancel the original random sampling mechanism of the GraphSAGE model, and adopt a deterministic neighbor selection strategy based on the conduction intensity of the disaster chain.

[0140] For each target node in the directed and weighted graph, first filter the effective neighbor nodes according to the directionality of the edge, and only retain the downstream flow area neighbors of the source area nodes and the downstream flow area neighbors of the flow area nodes, excluding the reverse nodes that have no disaster conduction significance.

[0141] Then sort the effective neighbor nodes according to the edge weight from large to small, and dynamically determine the number of neighbor nodes according to the node type, that is, the source area node selects the top N neighbors, and the flow area node selects the top M neighbors, where the values of N and M are based on the calibration of the influence range of the small watershed disaster chain, to ensure that the selected neighbors are all core related nodes that significantly contribute to the disaster risk conduction of the target node, and to avoid missing key conduction information or interference from irrelevant nodes caused by random sampling.

[0142] Step 403: Introduce a neighbor weight mechanism that fits the disaster chain influence mechanism, and convert the edge weight calculated in step 300 into the aggregated weight of the neighbor nodes. For each neighbor node of the target node, its weight value directly adopts the weight of the corresponding edge, and the larger the edge weight, the stronger the disaster conduction influence of the neighbor node on the target node, and the higher the contribution ratio in information aggregation.

[0143] At the same time, adjust the weight distribution combined with the neighbor node type, if the neighbor is a source area node, superimpose a correction coefficient corresponding to the loose deposit distribution ratio on the basis of the weight, to highlight the enhancement effect of the source supply capacity on the weight; if the neighbor is a flow area node, superimpose a correction coefficient corresponding to the accumulation intensity, to strengthen the influence of the flow conduction capacity on the weight. In the aggregation process, a weighted average aggregation function is adopted to weight and sum the feature vectors of each neighbor node according to the corresponding weight, to obtain the neighbor aggregated features, and to ensure that the aggregation result accurately reflects the differentiated disaster conduction contribution of different neighbors to the target node.

[0144] Step 404: Retain the self-loop information of the model and optimize its integration logic, highlighting the fundamental impact of the geographical node's own properties on the disaster chain. Add a self-loop edge to each geographical node, and the self-loop weight is determined according to the node type and the comprehensive impact weight. The self-loop weight of the source area node is set to 1.2 times its comprehensive impact weight, highlighting the core role of the source area's own stability in landslide initiation. The self-loop weight of the flow-through area node is set to 1.1 times its comprehensive impact weight, emphasizing the fundamental support of the flow-through area's own topographic and hydrological conditions on debris flow transport. During information aggregation, multiply the target node's own initial feature vector by the self-loop weight to obtain the self-loop feature vector, and process it in parallel with the neighbor aggregated feature vector to avoid representation bias caused by the model ignoring the node's own properties.

[0145] Step 405: Introduce residual information connection mechanism to solve the information decay problem in deep layer training, while adhering to the continuity feature of disaster chain transmission. After each layer of aggregation and feature conversion in the model, construct a residual connection path to add the output feature vector of the current layer and the input feature vector of the previous layer element-wise.

[0146] To adapt to the information retention needs of different links of the disaster chain, set the residual weight coefficient, which is higher for source area nodes than for flow-through area nodes, ensuring that the initial topographic and geological features of source area nodes and historical transmission information of flow-through area nodes are not weakened in deep network. Through residual connection, node representation contains both current layer aggregated neighbor information and historical level key features, achieving cross-layer continuous transmission of disaster chain information.

[0147] Step 406: Design a multi-stage aggregation and feature conversion structure to gradually optimize node representation effects.

[0148] The model is set to K layers of aggregation structure, with K value determined according to the transmission level of the small watershed disaster chain, ensuring coverage of the complete chain from source initiation to flow-through spread.

[0149] The first layer of aggregation only integrates the target node's direct neighbor information and self-loop information, outputting a first-order representation vector. Each subsequent layer of aggregation is based on the current layer representation vector, combined with the corresponding layer representation vector of the neighbor nodes, repeating the neighbor weight aggregation, self-loop information integration, and residual connection process to achieve progressive integration of multi-stage disaster transmission information. After each layer of aggregation, a linear conversion layer and an activation function are connected to perform dimension conversion and non-linear mapping of the aggregated features, enhancing the model's ability to capture complex correlation features of the disaster chain while avoiding computational redundancy caused by feature dimension expansion.

[0150] Step 407: Perform the full flow operation of information transmission and node characterization optimization, start model calculation according to the disaster conduction order of the directed weighted graph, first perform feature aggregation and characterization update on all source area nodes, then process upstream flow area nodes and downstream flow area nodes in turn, to ensure that the information transmission direction is consistent with the actual disaster chain transmission path.

[0151] After each layer operation is completed, the output characterization vector is standardized to ensure that the characterization vectors of each node are in the same numerical interval, eliminating the influence of dimension difference on subsequent calculations. After K layers of aggregation and optimization, the final characterization vector of each geographic node is obtained, which integrates its own attributes, multi-order neighbor disaster conduction information, and cross-layer residual information, fully describing the risk characteristics of the node in the landslide-mudflow disaster chain.

[0152] Step 408: Verify the adaptability of the node characterization results by comparing the final characterization vectors of disaster nodes and non-disaster nodes to test whether there is a significant distinction between the two, ensuring that the characterization results can effectively reflect the disaster risk differences of the nodes; at the same time, analyze the dimension characteristics of the characterization vector in combination with the node type to verify whether the characterization vectors of the source area nodes and the flow area nodes present feature differences matching their own functional attributes. If the distinction is insufficient or the feature deviation is large, adjust the number of neighbor selection, the self-loop weight coefficient, or the residual weight coefficient, and re-execute the characterization optimization process until the results meet the risk assessment requirements of the landslide-mudflow disaster chain, providing high-quality node characterization data for model training and risk prediction in step 500.

[0153] Further, with reference to Figure 3 , for step 500, to achieve sample construction and model training based on unit splitting feature matching terrain fitting in a small data area, the following is specified:

[0154] Step 501: Split the complex ground feature units of the existing disaster area to extract fine-grained unit-level basic samples.

[0155] Multiple regional data of landslide-mudflow disasters that have occurred at home and abroad need to be collected, which can be limited to complete regions or only small range data around disaster points.

[0156] According to the complex ground feature unit segmentation logic in step 100, each data is split into the smallest granularity of source area units and flow area units, with the source area unit area not exceeding 0.05 square kilometers and the flow area unit area not exceeding 0.1 square kilometers. Ensure that each fine-grained unit only contains a single terrain feature and a single disaster chain link, such as a steep slope source unit corresponding only to the landslide initiation link and a narrow ditch flow unit corresponding only to the mudflow transport link. For each fine-grained unit, three types of core information are extracted simultaneously.

[0157] The first type is the topographic feature of the unit itself, including slope, rock-soil type, relative height difference, and channel longitudinal gradient, which is consistent with the optimization factor set in step 200.

[0158] The second type is the correlation coefficient between the unit and adjacent units, which is based on the edge weight calculation logic in step 300, and the correlation strength and direction data between the unit and upstream and downstream units are retained.

[0159] The third type is the disaster label of the unit, indicating whether the unit has experienced a disaster and the link of the disaster chain, such as landslide in the source unit and debris flow in the flow unit. Through the above operations, a unit-level basic sample library containing thousands of fine units is formed.

[0160] Step 502: Establish a unit feature topographic type mapping index to prepare for target area fitting. For each fine unit in the unit-level basic sample library, perform double classification and labeling according to topographic dimensions and disaster chain links.

[0161] The topographic dimensions are divided into 27 topographic types according to slope levels, rock-soil stability, and convergence conditions. The disaster chain links are divided into four types according to source initiation, source supply, flow transport, and flow deposition.

[0162] Based on the classification results, a mapping index table is constructed, with each row corresponding to a fine unit and the columns containing topographic type codes, disaster link codes, unit feature vectors, adjacent unit correlation coefficients, and disaster labels. This index table supports quick retrieval and matching of samples according to topographic types and disaster links, solving the unit screening efficiency problem in subsequent target area terrain adaptation.

[0163] Step 503: Analyze the target area terrain partition to determine the target template for unit fitting.

[0164] The target area is divided into topographic zones according to the composite ground units generated in step 100. First, the source area and flow area are split into large units according to the function in step 100, and then each large unit is matched with the 27 topographic types in step 502 according to its slope, rock-soil, and convergence three-dimensional characteristics, forming the target area terrain partition template.

[0165] For example, the topographic features of a large unit in a source area of the target area are steep slope greater than 35 degrees, low rock-soil stability, and medium convergence condition, corresponding to a specific topographic type code; the features of a large unit in a flow area are medium slope 15-35 degrees, medium rock-soil stability, and strong convergence condition, corresponding to another topographic type code. At the same time, the correlation between units in the target area is recorded, such as the adjacency of source area A and flow area B, and the upstream and downstream connection of flow area B and flow area C, forming a target area unit correlation topology graph.

[0166] Step 504: Generate target region exclusive sample set by terrain matching and associated adaptive fitting. Based on the target region's terrain zoning template and associated topological graph, filter and combine the fine-grained units from the unit-level basic sample library.

[0167] In the single-unit matching link, for each minimum granularity small unit split by step 100 in the target region, retrieve the fine-grained units with the same code from the sample library according to its terrain type code. If the number of completely matched units is not less than 5, randomly select 3 to 5 as candidate units. If there is no completely matched unit, select units that match the two core features of slope and rock-soil body and have similar confluence conditions, adjust the confluence-related feature values through linear interpolation, such as adjusting the sample unit confluence intensity 0.8 to the target unit required confluence intensity 0.6, to ensure feature adaptation.

[0168] In the associated relationship fitting link, refer to the target region unit associated topological graph, and supplement the correlation coefficients of the adjacent units of the matched candidate units from the sample library. If target units A and B are adjacent, and the candidate unit of A comes from unit X in the sample library and the candidate unit of B comes from unit Y in the sample library, calculate the historical correlation coefficient of units X and Y. Then, combined with the actual distance between A and B in the target region, calibrate and correct the correlation coefficient of A and B according to the spatial adjacency threshold of step 301, to ensure that the correlation of the fitted units conforms to the conduction law of the actual terrain of the target region.

[0169] In the combination verification link, combine the fitted units according to the associated topological graph of the target region to form the target region fitted sample set. Check whether the combined sample set covers all terrain types of the target region, such as steep slope source area and gentle slope flow-through area. At the same time, confirm that there is no logical contradiction in the correlation of each unit, such as no reverse correlation of downstream unit pointing to upstream unit. If there is a lack or contradiction, repeat the single-unit matching and associated relationship fitting steps to supplement and adjust.

[0170] Step 505: Split the disaster chain tags of the fitted sample set to enhance data diversity and adaptability. For the target region fitted sample set generated by fitting, according to the unit function and correlation of step 503, split the traditional binary tags of disaster occurrence or non-occurrence into two-dimensional tags of link intensity.

[0171] Link tags label the disaster chain link to which the unit belongs according to the source start, source recharge, flow transport, and flow deposition, such as labeling the fitted source area unit as source start and the flow-through area unit as flow transport.

[0172] The intensity label refers to the disaster intensity of the candidate unit in the sample library, such as landslide scale, debris flow flow, and combines the topographic features of the target area unit, such as the intensity of steep slope unit is higher than that of gentle slope unit, and labels the intensity of three levels of strong, medium and weak. For example, the target area steep slope source unit is labeled as strong source start, and the gentle slope flow unit is labeled as weak flow transport. Through label splitting, the fitted sample set can cover the characteristics of different links and intensities of the disaster chain even if the number of units is limited, and the learning ability of the model to the details of the disaster chain of the target area is improved.

[0173] Step 506: Based on the phased model training of the fitted sample set, the adaptability of the target area is strengthened.

[0174] The first stage carries out unit-level pre-training, taking a single unit in the fitted sample set as the training unit, and inputting the unit feature vector and link label. The improved GraphSAGE model designed in step 400 is used to learn the mapping relationship from the topographic features of the unit to the disaster link, and the feature capture ability of the model to a single unit is optimized. When the unit-level prediction accuracy is stable at more than 75%, the pre-training is stopped.

[0175] The second stage carries out correlation-level fine-tuning, introduces the unit correlation index in the fitted sample set, and constructs a unit correlation subgraph according to the correlation topology graph of the target area. The pre-trained model parameters are used as the initial value, and the correlation subgraph is used to train the model to learn the rules of correlation between units to disaster transmission, such as strong start of source unit plus high correlation index corresponding to strong transport of flow unit. The residual connection and neighbor weight mechanism of step 400 are retained during fine-tuning to ensure that the correlation information does not decay. When the F1 score of the correlation-level prediction is stable at more than 70%, the fine-tuning is stopped.

[0176] The third stage carries out regional-level verification, randomly selects 20% of the unit combinations from the fitted sample set as a verification sub-region to simulate the complete topography of the target area. The fine-tuned model is used to predict the disaster chain risk of the verification sub-region, and the difference between the prediction result and the label of the fitted sample set is compared. If the difference rate exceeds 15%, return to step 504 to adjust the unit fitting logic, such as supplementing candidate units of specific topographic types, until the difference rate does not exceed 15%, to ensure that the model adapts to the overall topographic features of the target area.

[0177] Step 507: Perform risk prediction and spatial mapping of target area, consistent with sample fitting logic. Split the target area into composite land units generated in step 100, according to the terrain zoning template in step 503. Input the trained improved GraphSAGE model, which follows the process of single-unit feature extraction, inter-unit correlation calculation, and multi-link risk determination, to output the disaster chain risk index for each unit, which integrates the comprehensive value of link label and intensity label. Divide the units into four risk levels: extremely high, high, medium, and low, according to the risk index. Combine the unit correlation topology graph of the target area to generate a risk spatial map. Label the unit function type source area or flow area and core disaster-causing factors during mapping, such as extremely high-risk steep slope source units with strong activation, to ensure that the risk map not only reflects unit-level risk but also demonstrates the correlation of disaster chain material sources, flow, and transmission, providing precise spatial guidance for subsequent disaster prevention deployment.

[0178] Step 508: Dynamically update unit sample library and model parameters to continuously optimize prediction effect. Regularly collect new environmental data in the target area, such as annual vegetation changes and rainfall data, and new disaster unit data from surrounding areas. Split the new disaster units into fine units according to step 501 and supplement them to the unit-level basic sample library. At the same time, repeat steps 503 to 506 based on the new data to fine-tune the fitting sample set and model parameters of the target area, ensuring that the model can continuously adapt to environmental changes in the target area and maintain the accuracy of risk prediction in the long term.

[0179] This embodiment provides an application example: landslide-debris flow disaster chain risk assessment in a small watershed in a mountainous area

[0180] Select a small watershed in a mountainous area in southwest China as the evaluation area. The area is about 80 square kilometers, with mainly medium and steep slopes, dense gullies, and scarce historical disaster data, with only 3 small-scale landslide-debris flow events recorded, which meets the evaluation scenario of few data regions.

[0181] Based on the technical process of this application, the watershed is first split into 1260 composite land units according to step 100, 12 core environmental factors are selected through step 200, a directed weighted graph is constructed in step 300, and the improved GraphSAGE model in step 400 is trained. The sample set is constructed according to the unit splitting-feature matching-terrain fitting mechanism in step 500, and the fine unit data of five different terrain regions are fused to generate the fitting.

[0182] At the same time, referring to Figures 4-7 , traditional machine learning models MLP, RF, and SVM are used as a comparison, and parallel testing is carried out under the same sample set and evaluation indicators. The performance indicators are compared as follows:

[0183] Precision: 0.895 for the method of the present application, 0.823 for RF, 0.801 for SVM, 0.786 for MLP;

[0184] Precision: 0.858 for the method of the present application, 0.812 for RF, 0.793 for SVM, 0.765 for MLP;

[0185] Recall: 0.885 for the method of the present application, 0.784 for RF, 0.752 for SVM, 0.731 for MLP;

[0186] F1 score: 0.820 for the method of the present application, 0.728 for RF, 0.705 for SVM, 0.689 for MLP;

[0187] AUC value: 0.854 for the method of the present application, 0.796 for RF, 0.778 for SVM, 0.753 for MLP;

[0188] Kappa coefficient: 0.816 for the method of the present application, 0.735 for RF, 0.712 for SVM, 0.690 for MLP.

[0189] The actual evaluation result shows that 42 high-risk units are identified in the small watershed. Among them, the method of the present application only misses 2 high-risk units, while RF misses 7, SVM misses 9, and MLP misses 11. The missed high-risk units are mostly in the complex terrain of steep slope source area and the junction of channel buffer zone. The method of the present application effectively reduces the risk of missing judgment in such areas by strengthening the capture of unit correlation characteristics and fitting the disaster chain links, providing more accurate spatial guidance for disaster prevention and monitoring deployment in the watershed.

[0190] The present application is not limited to the above-mentioned optional embodiments, and anyone can derive other various forms of products under the inspiration of the present application. The above specific embodiments should not be understood as limiting the scope of protection of the present application, and the scope of protection of the present application should be defined by the claims, and the specification can be used to explain the claims.

Claims

1. A landslide-debris flow disaster chain evaluation method based on a graph neural network, characterized in that, The method comprises the following steps: Step 100: generating a composite land unit by composite curvature field watershed segmentation and hydrological analysis, and abstracting the composite land unit into a geographic node; Step 200: constructing a node feature corresponding to the geographic node, and verifying the independence of each environmental factor in the node feature by using a mutual information method; Step 300: constructing a directed weighted graph according to the environmental factor and spatial similarity, wherein the directed weighted graph is used to reflect the spatial transmission relationship of the landslide-debris flow disaster chain; Step 400: designing an improved GraphSAGE model, and performing information transmission and node representation optimization on the information in the directed weighted graph by using the improved GraphSAGE model; Step 500: training and verifying the improved GraphSAGE model, and completing risk prediction and spatial mapping of the landslide-debris flow disaster chain; The process of designing the improved GraphSAGE model in step 400 comprises: introducing a neighbor weight mechanism, reserving self-loop information, reserving residual information, and canceling a random sampling operation for the GraphSAGE model; The neighbor weight mechanism is specifically introduced as follows: different neighbor nodes of each geographic node are assigned with differentiated weights according to the influence degree of the neighbor nodes on the geographic node; The self-loop information is specifically reserved as follows: when the improved GraphSAGE model performs information transmission, the feature information of each geographic node itself is reserved, and the information contribution of the geographic node to itself is not excluded; and the residual information is specifically reserved as follows: when the improved GraphSAGE model performs node representation optimization, the residual connection is performed between the node representation result of the previous round and the node representation result of the current round, and the historical representation information is reserved.

2. The landslide-debris flow disaster chain evaluation method based on a graph neural network according to claim 1, characterized in that, The process of generating the composite land unit in step 100 comprises: first obtaining a slope unit by composite curvature field watershed segmentation, and then introducing a gully buffer zone on the basis of the slope unit to form the composite land unit.

3. The landslide-debris flow disaster chain evaluation method based on a graph neural network according to claim 2, characterized in that, After the composite land unit is formed in step 100, an operation of quantifying the source area features and flow channel features corresponding to the composite land unit is further included.

4. The landslide-debris flow disaster chain evaluation method based on a graph neural network according to claim 1, characterized in that, The process of constructing the node feature in step 200 comprises: collecting multi-source heterogeneous environmental factors, and taking the multi-source heterogeneous environmental factors as initial input factors of the node feature; and the independence of each environmental factor is verified by using a mutual information method, and the correlation degree between each initial input factor is calculated by using a normalized mutual information method.

5. The landslide-debris flow disaster chain evaluation method based on a graph neural network according to claim 4, characterized in that, After the correlation degree between each initial input factor is calculated in step 200, an operation of removing a strongly dependent factor according to the correlation degree is further included, wherein the strongly dependent factor is an initial input factor with a correlation degree exceeding a preset threshold.

6. The landslide-debris flow disaster chain evaluation method based on a graph neural network according to claim 1, characterized in that, Before the improved GraphSAGE model performs information transmission and node representation optimization in step 400, an operation of constructing a disaster chain sample set is further included, and the construction of the disaster chain sample set adopts a hierarchical sampling strategy, and positive and negative samples are counted according to environmental units to ensure balanced sample distribution.

7. The method of claim 6, wherein the method is based on a graph neural network. The process of constructing the disaster chain sample set further includes an operation of splitting disaster chain labels, and the data diversity of the disaster chain sample set is enhanced through the label splitting.

Images