Highway line selection method and system based on space-air-ground internal intelligence integrated study and judgment

Abstract

The invention discloses a road line selection method and system based on space-air-ground internal intelligence integrated study and judgment, and relates to the technical field of road line selection, and the method comprises the steps: employing the InSAR and remote sensing image interpretation technology to carry out the preliminary recognition and screening of the geological hazard hidden danger of a corridor zone from a large scale range; an airborne laser radar technology is adopted to carry out detailed investigation on geological disasters of the corridor from a small-scale range; checking geological disasters around the preliminary line by adopting comprehensive survey means of geological survey, geophysical prospecting test and drilling verification; risk evaluation is carried out on geological disaster hidden danger points in the preliminary line corridor zone range, and a risk evaluation model is constructed; a plurality of suggested lines are given after the lines are optimized by a geological engineer and a line engineer, and an optimal planning line is calculated by adopting an intelligent expert research and judgment model. The technical problem that in the prior art, due to dependence on single geological disaster data, scientificity and objectivity of line selection are limited is solved.

Description

Technical Field

[0001] This invention relates to the field of highway route selection technology, and in particular to a highway route selection method and system based on integrated air-space-ground-indoor intelligent analysis. Background Technology

[0002] More and more highways in my country are extending into mountainous areas prone to geological disasters, and the focus of highway construction in Hubei Province is also shifting towards the geologically complex mountainous regions of western and northwestern Hubei. Due to the constraints of these mountainous areas' complex terrain and geological conditions, the variety of complex terrain and adverse geological phenomena that need to be considered in highway route selection will significantly increase. An unreasonable route will expose roadbed slopes, bridges, and tunnels to significant geological disaster risks and structural safety risks during construction, potentially leading to substantial economic losses.

[0003] Due to the complexity, diversity, and concealment of geological hazards, accurate identification is extremely difficult. Traditional identification methods mainly rely on manual on-site investigations and borehole exploration. These traditional single-hazard investigation methods fail to fully utilize the advantages of technologies such as remote sensing and drone aerial photography, resulting in huge investments of manpower and resources. Furthermore, they suffer from incomplete and insufficient investigations, making it difficult to meet the needs of geological hazard prevention and control. At the same time, existing route selection methods mostly depend on data from single geological hazards, and often suffer from independent decision-making and a lack of data complementarity. Route selection often depends on subjective human will, lacking objective scientific basis, leading to the underutilization of multi-source data and limiting the scientific rigor and rationality of route selection. Summary of the Invention

[0004] To address the aforementioned technical problems in existing technologies, this invention provides a highway route selection method and system based on integrated air-space-ground intelligent analysis. The technical solution is as follows:

[0005] On the one hand, a highway route selection method based on integrated air-space-ground-indoor intelligent analysis is provided. The method includes: preliminary zoning of the geological hazard susceptibility within the highway corridor based on the overall highway route planned during the engineering feasibility study phase, and preliminary marking based on historical geological hazard data to obtain a preliminary spatial distribution map of geological hazards along the highway; preliminary identification and screening of potential geological hazards in the corridor area based on InSAR and remote sensing image interpretation technology to obtain a large-scale spatial distribution map of potential geological hazards; detailed investigation of the large-scale spatial distribution map of potential geological hazards based on UAV-borne lidar and real-time dynamic positioning technology to obtain a small-scale spatial distribution map of potential geological hazards; and preliminary route design based on the small-scale spatial distribution map of potential geological hazards to obtain a preliminary... The design phase involves defining the local route corridor area; verifying surrounding geological hazards within the local route corridor area using comprehensive survey methods, and constructing a geological condition information database; these comprehensive survey methods include geological surveys, geophysical testing, and drilling verification; based on the geological condition information database, conducting local geological risk assessments of potential geological hazard points within the local route corridor area using qualitative methods, and obtaining a risk assessment model; based on the risk assessment model, making preliminary judgments on risk sources and risk zones along the local route corridor area, and adjusting and optimizing the route within the local route corridor area to obtain multiple routes to be planned; constructing an intelligent expert judgment model using a combination of qualitative and quantitative methods, and judging and comparing the multiple routes to be planned based on the intelligent expert judgment model to obtain the target planned route.

[0006] Optionally, based on the overall highway route planned during the feasibility study phase, the susceptibility to geological hazards within the highway corridor is preliminarily divided into zones, and preliminary marking is performed based on historical geological hazard data to obtain a preliminary spatial distribution map of geological hazards along the highway. This includes: classifying the susceptibility to geological hazards within the highway corridor into multiple regions according to the topographical features that easily induce geological hazards, based on the overall highway route planned during the feasibility study phase; the topographical features include: rivers, riverbanks, mountains, and karst areas; marking historical hazards in each region; the historical hazards include: landslides, debris flows, collapses, and unstable rock masses; generating image maps of the highway route within a predetermined range based on satellite data, and marking the historical hazards in the image maps to obtain a preliminary spatial distribution map of geological hazards along the highway.

[0007] Optionally, based on InSAR and remote sensing image interpretation technology, preliminary identification and screening of potential geological hazards in the corridor zone of the preliminary geological hazard spatial distribution map along the highway are performed to obtain a large-scale spatial distribution map of potential geological hazards. This includes: based on synthetic aperture radar interferometry, analyzing the interference effect between radar images of the target area acquired at different times to identify the displacement of the ground surface between two imaging sessions, obtaining multiple SAR images of the same area at different time points; the target area is the area corresponding to the corridor zone in the preliminary geological hazard spatial distribution map along the highway; registration, interferometry, filtering, and unwrapping calculations are performed on the multiple SAR images to obtain the deformation of the target area over the entire time period. The deformation time series is used to obtain surface deformation monitoring data; based on the surface deformation monitoring data, the three-dimensional coordinate information of the target area before and after deformation is obtained, and based on the three-dimensional coordinate information of the target area before and after deformation, the deformation characteristics, volume, and drop information of the unfavorable geological body are obtained; based on the deformation characteristics of the unfavorable geological body, the type of unfavorable geological body is identified; the types of unfavorable geological bodies include: landslides, debris flows, collapses, and unstable rock masses; based on the type of unfavorable geological body and the volume and drop information, the impact range of geological hazards is estimated; the impact range of geological hazards is marked in the multi-scene SAR imagery to obtain a large-scale spatial distribution map of geological hazard risks.

[0008] Optionally, a detailed investigation of the large-scale geological hazard spatial distribution map is conducted based on UAV-borne LiDAR and real-time dynamic positioning technology to obtain a small-scale geological hazard spatial distribution map. This includes: using UAV-borne LiDAR and UAV-borne camera aerial photography technology, flying close to the geological hazard points in the large-scale geological hazard spatial distribution map to obtain image data of unfavorable geological bodies within the geological hazard and 3D point cloud data acquired based on LiDAR; using SFM 3D image reconstruction technology to generate a 3D model; based on the 3D model, obtaining the attitude information of exposed structural surfaces of unfavorable geological bodies in different geological hazard areas; embedding the 3D model into the image map of the large-scale geological hazard spatial distribution map, and marking the attitude information of exposed structural surfaces of unfavorable geological bodies in the 3D model to obtain the small-scale geological hazard spatial distribution map.

[0009] Optionally, a geological condition information database is constructed by verifying the surrounding geological hazards within the local railway corridor area based on comprehensive exploration methods. This includes: obtaining first geological information of the local railway corridor area using geological survey methods; the first geological information includes: stratigraphic lithology, landform type, hydrology and meteorology, geological structure, vegetation cover information, and information on the impact of human activities; obtaining second geological information by conducting ground tests on the local railway corridor area using geophysical testing methods such as high-density electrical resistivity tomography, CSATM magnetotelluric method, micro-motion surface wave method, and 3D seismic method; obtaining third geological information by conducting internal geological exploration of the soil and rock at potential geological hazard points within the local railway corridor area using drilling verification methods such as internal drilling, comprehensive well logging, and deep displacement monitoring; and constructing a geological condition information database based on the first, second, and third geological information.

[0010] Optionally, based on the geological condition information database, a qualitative method is used to conduct a local geological risk assessment of the geological hazard points within the local route corridor, resulting in a risk assessment model. This model includes: constructing multiple risk indicators based on the geological condition information database, and assigning corresponding weights to each risk indicator; the risk indicators include: topographic indicators, geological structure indicators, stratigraphic lithology indicators, hydrogeological condition indicators, and adverse geological body indicators; and conducting expert panel scoring of the geological hazard points within the local route corridor based on the multiple risk indicators to construct the risk assessment model.

[0011] Optionally, based on the risk assessment model, a preliminary judgment is made on the risk sources and risk zones along the local route corridor, and the routes within the local route corridor are adjusted and optimized to obtain multiple routes to be planned. This includes: based on the risk assessment model, a preliminary judgment is made on the risk sources and risk zones along the local route corridor to obtain the risk level of multiple risk zones; the multiple risk zones are ranked based on the risk level of each risk zone; the routes within the local route corridor are adjusted and optimized based on the ranking of the multiple risk zones to obtain multiple routes to be planned; wherein the optimization objective is to minimize the sum of the risk levels of the routes passing through the risk zones.

[0012] Optionally, an intelligent expert judgment model is constructed using a combination of qualitative and quantitative methods. Based on the intelligent expert judgment model, the multiple routes to be planned are judged and compared to obtain the target planned route. This includes: determining multiple impact indicators based on the risk assessment model; the multiple image indicators include multiple technical indicators, multiple geological indicators, and multiple economic indicators; screening multiple sensitive indicators based on the multiple impact indicators; constructing an intelligent expert judgment model using a combination of qualitative and quantitative methods based on the multiple sensitive indicators; and judging and comparing the multiple routes to be planned based on the intelligent expert judgment model to obtain the target planned route.

[0013] Optionally, a combined qualitative and quantitative approach is used to construct an intelligent expert judgment model, including: constructing an initial model structure based on an optimal route layer, an influencing indicator layer, and a proposed route layer; establishing a first comparison matrix for the influencing indicator layer; the first comparison matrix includes... , Indicator Factors and Compared to the importance of what is obtained, The first factor influencing highway route selection is defined as: Based on the first comparison matrix, a first importance matrix is ​​constructed relative to the optimal route layer O, whereby the influencing index layer is the i-th factor. This first importance matrix includes: Based on the first importance matrix, the largest eigenvalue and the corresponding first eigenvector are obtained, and normalized to obtain a normalized first eigenvector. Based on the magnitude of the elements in the normalized first eigenvector, the importance ranking of each factor in the influencing index layer is obtained; each element in the normalized first eigenvector represents the weight of the influencing index layer relative to each factor in the optimal route layer; a second comparison matrix for the proposed route layer is established; the second comparison matrix includes... , Indicator Factors With factors Compared to the importance of what is obtained, For the proposed first i One alternative route; based on the second comparison matrix, the proposed route layer is obtained relative to the first influencing index layer. i The second importance matrix of the factors; the second importance matrix includes: Based on the second importance matrix, the largest eigenvalues ​​and corresponding second eigenvectors are obtained and normalized to obtain normalized second eigenvectors. Each element in the normalized second eigenvector represents the weight of the proposed route layer relative to each factor in the optimal route layer. Based on the magnitude of the elements in the normalized second eigenvector, the comprehensive weight of each factor in the proposed route layer to the optimal route layer is obtained. The comprehensive weight includes: In the formula, n i Let n be the comprehensive weight of the i-th factor in the proposed route layer to the optimal route layer, and n be the total number of factors in the proposed route layer. i,j m is an element in the normalized second eigenvector. j The elements are those in the normalized first feature vector; based on the comprehensive weight, the order of each candidate route in the proposed route layer is determined.

[0014] On the other hand, a highway route selection system based on integrated air-space-ground-indoor intelligent analysis is also provided, applied to the highway route selection method based on integrated air-space-ground-indoor intelligent analysis provided in the embodiments of the present invention. The system includes: a preliminary zoning module, a first screening module, a second screening module, a preliminary design module, a construction module, an evaluation module, an optimization module, and an analysis module. The preliminary zoning module is used to perform preliminary zoning of the geological hazard susceptibility within the highway corridor based on the macro-route of the entire highway planned in the engineering feasibility study stage, and to perform preliminary marking based on historical geological hazard data to obtain a preliminary spatial distribution map of geological hazards along the highway. The first screening module is used to perform preliminary identification and screening of corridor geological hazard hazards in the preliminary spatial distribution map of geological hazards along the highway based on InSAR and remote sensing image interpretation technology to obtain a large-scale spatial distribution map of geological hazard hazards. The second screening module is used to conduct a detailed investigation of the large-scale spatial distribution map of geological hazard hazards based on UAV-borne lidar and real-time dynamic positioning technology to obtain a small-scale spatial distribution map of geological hazard hazards. The system comprises the following modules: a spatial distribution map of geological hazards; a preliminary design module for preliminary route design based on the small-scale spatial distribution map of geological hazards, resulting in a local route corridor range for the preliminary design stage; a construction module for verifying the surrounding geological hazards within the local route corridor range using comprehensive exploration methods, including geological surveys, geophysical testing, and drilling verification; an evaluation module for conducting local geological risk assessments of geological hazard points within the local route corridor range using qualitative methods based on the geological condition information database, resulting in a risk assessment model; an optimization module for making preliminary judgments on risk sources and risk zones along the local route corridor range based on the risk assessment model, and adjusting and optimizing the route within the local route corridor range to obtain multiple routes to be planned; and a judgment module for constructing an intelligent expert judgment model using a combination of qualitative and quantitative methods, and judging and comparing the multiple routes to be planned based on the intelligent expert judgment model to obtain the target planned route.

[0015] This invention provides a method and system for highway route selection based on integrated air-space-ground-indoor intelligent analysis. By integrating multi-source information such as InSAR data, airborne lidar data, geological surveys, geophysical tests, and drilling verification, a risk assessment model and an intelligent expert analysis model are constructed. This solves the problems of data bias and strong subjectivity when selecting routes based on a single factor. By comprehensively considering multiple influencing factors and comparing the importance of any two factors to calculate their weights, the scientific nature and timeliness of highway route selection are improved. Attached Figure Description

[0016] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0017] Figure 1 This is a flowchart of a highway route selection method based on integrated air-space-ground-indoor intelligence analysis provided by an embodiment of the present invention; Figure 2 This is a preliminary spatial distribution map of geological hazards along a highway provided by an embodiment of the present invention; Figure 3 This is a large-scale spatial distribution map of potential geological hazards provided by an embodiment of the present invention; Figure 4 This is a small-scale spatial distribution map of potential geological hazards provided by an embodiment of the present invention; Figure 5 This is a schematic diagram of the structure of an intelligent expert judgment model provided in an embodiment of the present invention; Figure 6 This is a schematic diagram of a highway route selection system based on integrated air-space-ground-indoor intelligent analysis provided in an embodiment of the present invention. Detailed Implementation

[0018] The technical solution of the present invention will now be described with reference to the accompanying drawings.

[0019] In embodiments of the present invention, words such as "exemplarily," "for example," etc., are used to indicate that something is an example, illustration, or description. Any embodiment or design described as "exemplary" in the present invention should not be construed as being more preferred or advantageous than other embodiments or designs. Specifically, the use of the word "exemplary" is intended to present the concept in a concrete manner. Furthermore, in embodiments of the present invention, the meaning expressed by "and / or" can be both, or either one.

[0020] To make the technical problems, technical solutions and advantages of the present invention clearer, a detailed description will be given below in conjunction with the accompanying drawings and specific embodiments.

[0021] Figure 1 This is a flowchart of a highway route selection method based on integrated air-space-ground-indoor intelligent analysis according to an embodiment of the present invention. Figure 1 As shown, the method specifically includes the following steps:

[0022] Step S102: Based on the overall highway route planned in the feasibility study stage, the geological hazard susceptibility within the highway corridor is preliminarily divided into zones, and preliminary marking is performed based on historical geological hazard data to obtain a preliminary spatial distribution map of geological hazards along the highway.

[0023] Step S104: Based on InSAR and remote sensing image interpretation technology, preliminary identification and screening of potential geological hazards in the corridor zone of the preliminary geological hazard spatial distribution map along the highway are carried out to obtain a large-scale spatial distribution map of potential geological hazards.

[0024] Step S106: Based on UAV-borne lidar and real-time dynamic positioning technology, a detailed investigation of the spatial distribution map of large-scale geological hazard risks is conducted to obtain a spatial distribution map of small-scale geological hazard risks.

[0025] Step S108: Based on the spatial distribution map of small-scale geological hazard risks, a preliminary route design is carried out to obtain the local route corridor range in the preliminary design stage.

[0026] Step S110: Verify the surrounding geological hazards within the local corridor area based on comprehensive exploration methods and construct a geological condition information database; the comprehensive exploration methods include geological survey, geophysical testing and drilling verification.

[0027] Step S112: Based on the geological condition information database, a qualitative method is used to conduct a local geological risk assessment of the geological hazard points within the local route corridor area, and a risk assessment model is obtained.

[0028] Step S114: Based on the risk assessment model, a preliminary judgment is made on the risk sources and risk zones along the local route corridor, and the routes within the local route corridor are adjusted and optimized to obtain multiple routes to be planned.

[0029] Step S116: A smart expert judgment model is constructed using a combination of qualitative and quantitative methods. Based on the smart expert judgment model, multiple routes to be planned are judged and compared to obtain the target planned route.

[0030] Specifically, step S102 further includes the following steps: Step S1021: Based on the overall highway route planned in the feasibility study stage, the terrain and geomorphological features that are prone to geological disasters are classified into multiple areas within the highway corridor. The terrain and geomorphological features include: rivers, riverbanks, mountains, and karst areas. Step S1022: Mark historical disasters in each area; historical disasters include: landslides, debris flows, collapses, and unstable rock masses; Step S1023: Generate an image map of the entire highway's macroscopic route within a predetermined range based on satellite data, and mark historical disasters on the image map to obtain a preliminary spatial distribution map of geological disasters along the highway, such as... Figure 2 As shown. Among them, Figure 2 This is a preliminary spatial distribution map of geological hazards along a highway, provided according to an embodiment of the present invention.

[0031] Among them, historical disasters refer to disasters that have occurred in the history of the highway corridor.

[0032] Specifically, step S104 further includes the following steps: Step S1041: Based on the synthetic aperture radar interferometry (InSAR) technology, analyze the interference effect between radar images of the target area acquired at different times to identify the displacement of the ground surface between the two imaging, and obtain multiple SAR images of the same area at different time points; the target area is the area corresponding to the corridor zone in the preliminary spatial distribution map of geological hazards along the highway. Step S1042: Register, interferometrically, filter, and unwrap calculations are performed on multiple SAR images to obtain the deformation time series of the target area for the entire time period, and surface deformation monitoring data are obtained based on the deformation time series. Step S1043: Based on the surface deformation monitoring data, obtain the three-dimensional coordinate information of the target area before and after deformation, and based on the three-dimensional coordinate information of the target area before and after deformation, obtain the deformation characteristics, volume and drop information of the unfavorable geological body. Step S1044: Based on the deformation characteristics of the unfavorable geological body, identify the type of unfavorable geological body; the types of unfavorable geological bodies include: landslides, debris flows, collapses, and unstable rock masses; Step S1045: Based on the type, volume, and elevation information of the adverse geological body, estimate the impact range of the geological hazard; Step S1046: Mark the area affected by geological disasters in multiple SAR images to obtain a large-scale spatial distribution map of potential geological disaster hazards.

[0033] Figure 3 This is a large-scale spatial distribution map of potential geological hazards provided by an embodiment of the present invention.

[0034] Specifically, step S106 further includes the following steps: Step S1061: Based on UAV-borne LiDAR and UAV-borne camera aerial photography technology, fly close to the geological hazard points in the large-scale geological hazard spatial distribution map to obtain image data of unfavorable geological bodies in the geological hazard and three-dimensional point cloud data based on LiDAR. Use SFM three-dimensional image reconstruction technology to generate a three-dimensional model. Step S1062: Based on the three-dimensional model, obtain the attitude information of exposed structural surfaces of unfavorable geological bodies in different geological hazard risks; Step S1063: The three-dimensional model is embedded into the image map of the large-scale spatial distribution map of geological hazards, and the occurrence information of exposed structural surfaces of adverse geological bodies is marked in the three-dimensional model to obtain the small-scale spatial distribution map of geological hazards.

[0035] Figure 4 This is a small-scale spatial distribution map of potential geological hazards provided by an embodiment of the present invention.

[0036] Specifically, step S110 further includes the following steps: Step S1101: Using geological survey methods, obtain the first geological information of the local route corridor zone; the first geological information includes: stratigraphic lithology, landform type, hydrology and meteorology, geological structure, vegetation cover information and information on the impact of human activities; Step S1102: Using geophysical testing methods, ground tests are conducted on the local corridor area using high-density electrical resistivity tomography, CSATM magnetotelluric method, micro-motion surface wave method, and three-dimensional seismic method to obtain second geological information. Step S1103: Using the drilling verification method, internal geological exploration of the soil and rock is carried out by means of internal drilling, comprehensive logging and deep displacement monitoring of potential geological disaster points in the local line corridor area to obtain third geological information. Step S1104: Based on the first geological information, the second geological information, and the third geological information, construct a geological condition information database containing multi-source heterogeneous information from the ground and interior.

[0037] Specifically, step S112 further includes the following steps: Step S1121: Based on the geological condition information database, construct multiple risk indicators and assign corresponding weights to each risk indicator; the risk indicators include: topographic and geomorphological indicators, geological structure indicators, stratigraphic and lithological indicators, hydrogeological condition indicators, and adverse geological body indicators. Preferably, the topographic and geomorphological indicators include: altitude, relative elevation, and slope. Geological structural indicators include: active faults, dead faults, and tectonic stress; Stratigraphic lithological indicators include: loose deposits, fractured rock masses, and intact rock masses; Hydrogeological indicators include: the distribution range of groundwater, the distribution range of surface water, and the corrosivity of water to structures; Unfavorable geological indicators include: landslides, debris flows, collapses, and unstable rock masses; Step S1122: Based on multiple risk indicators, an expert panel scores the geological hazard points within the local corridor area and constructs a risk assessment model.

[0038] Specifically, the risk assessment model includes: Based on expert experience, the weights of topographic and geomorphological indicators, geological structure indicators, stratigraphic lithology indicators, hydrogeological condition indicators, and adverse geological body indicators were obtained respectively, and they conform to the relationship formula (1): (1) In the formula: Weights for topographic and geomorphological indicators Weights of geological structural indicators Weights of stratigraphic lithology indicators The weights of hydrogeological condition indicators The weights of the indicators for unfavorable geological bodies; Then, based on expert experience, the weights of altitude, relative elevation, and slope gradient are obtained, and they conform to the relationship formula (2): (2) In the formula: Weighting for altitude Weighting relative elevation The weight of the slope gradient; Based on expert experience, the weights of active faults, dead faults, and tectonic stresses are obtained, and they conform to the relationship formula (3): (3) In the formula: Weight of active faults, Weights of dead faults The weight of the structural stress; Based on expert experience, the weights of loose sedimentary layers, fractured rock masses, and intact rock masses are obtained, and they conform to the relationship formula (4): (4) In the formula: For the weight of loosely packed layers, For the weight of the fractured rock mass, The weight of the intact rock mass; Based on expert experience, the weights of groundwater distribution range, surface water distribution range, and water corrosivity to structures were obtained, and they conformed to the relationship formula (5): (5) In the formula: Weights for groundwater distribution range Weights for the distribution range of surface water The weight of the corrosiveness of water to a structure; Based on expert experience, the weights of landslides, debris flows, collapses, and unstable rock masses are obtained, and they conform to the relationship formula (6): (6) In the formula: For the weight of landslides, Weight of debris flow For the collapse of weight, The weight of the unstable rock mass; Experts scored the altitude, relative elevation, and slope gradient to obtain scores for topographic indicators. (7) In the formula: Fractions of altitude Fractions for relative elevation It represents the fraction of the slope. By having experts score active faults, dead faults, and tectonic stress, scores for geological structural indicators are obtained: (8) In the formula: For the fraction of active faults, The fraction of dead faults, This represents the fraction of structural stress. By having experts score the loose sedimentary layers, fractured rock masses, and intact rock masses, scores for the stratigraphic lithology indicators are obtained: (9) In the formula: For the fraction of loosely packed layers, The fraction of the fractured rock mass The fraction of intact rock mass; By having experts score the distribution range of groundwater, the distribution range of surface water, and the corrosiveness of water to structures, scores for hydrogeological conditions are obtained: (10) In the formula: Fractions of the distribution range of groundwater Fractions representing the distribution range of surface water The score representing the corrosiveness of water to the structure; Experts scored landslides, debris flows, collapses, and unstable rock masses to obtain scores for unfavorable geological indicators: (11) In the formula: For the landslide score, For the fraction of debris flow, For the collapsed score, The fraction representing the unstable rock mass; Based on the scores of topographic and geomorphological indicators, geological structure indicators, stratigraphic lithology indicators, hydrogeological condition indicators, and adverse geological body indicators, and their respective weights, a risk assessment model for a specific area is obtained: (12) Specifically, step S114 further includes the following steps: Step S1141: Based on the risk assessment model, make a preliminary judgment on the risk sources and risk zones along the local line corridor and obtain the risk level of multiple risk zones. Step S1142: Sort multiple risk partitions based on the risk level of each risk partition; Step S1143: Based on the ranking of multiple risk zones, the routes within the local route corridor are adjusted and optimized to obtain multiple routes to be planned; the optimization objective is to minimize the sum of the risk levels of the routes passing through the risk zones.

[0039] Specifically, step S116 further includes the following steps: Step S1161: Based on the risk assessment model, determine multiple impact indicators; these multiple image indicators include: multiple technical indicators, multiple geological indicators, and multiple economic indicators. Preferably, the technical specifications include: mileage, horizontal alignment, longitudinal profile, cross profile, burial depth, and slope; Geological indicators include: topographic and geomorphological indicators, geological structure indicators, stratigraphic and lithological indicators, hydrogeological condition indicators, and indicators of unfavorable geological bodies; Economic indicators include: project cost, demolition cost, and environmental impact losses; Step S1162: Based on multiple influencing indicators, screen multiple sensitive indicators; Step S1163: Based on multiple sensitive indicators, construct an intelligent expert judgment model using a combination of qualitative and quantitative methods; Step S1164: Based on the intelligent expert judgment model, multiple routes to be planned are judged and compared to obtain the target planned route.

[0040] Figure 5 This is a schematic diagram of the structure of an intelligent expert judgment model provided according to an embodiment of the present invention. Figure 5 As shown, step S1163 further includes the following steps:

[0041] Step S1: Construct the initial model structure based on the optimal route layer, the influencing index layer, and the proposed route layer; Specifically, the optimal path layer, using This indicates the final selected highway route; Influence indicator layer ,use This indicates the factors that influence highway route selection; Proposed line layer ,use This indicates the initial list of n alternative routes. Step S2, establish the first comparison matrix affecting the indicator layer; the first comparison matrix includes , Indicator Factors and Compared to the importance of what is obtained, This is the i-th factor influencing highway route selection; Specifically, M is a positively reciprocal matrix. The degree of importance follows these rules: When factors and When equally important, ; When factors Compare When it is slightly important, ; When factors Compare When it is obviously important, ; When factors Compare When it is extremely important, ; When factors Compare When extremely important and strong, ; When factors Compare Compared to when it is in the middle of the above adjacent importance levels, they are respectively , , , ; When factors With factors The importance of the comparison is At that time, then factors With factors The importance of the comparison is ; Step S3: Based on the first comparison matrix, construct the first importance matrix of the influence index layer relative to the optimal line layer O; the first importance matrix includes: Step S4: Based on the first importance matrix, obtain the largest eigenvalue and the corresponding first eigenvector, and normalize them to obtain the normalized first eigenvector. Based on the weights in the normalized first eigenvector, obtain the importance ranking of each factor influencing the indicator layer; where the first eigenvector... The elements in the layer are the influencing indicators. Relative to the optimal line layer Various factors Corresponding weights .

[0042] Step S5: Establish the second comparison matrix for the proposed route layer; the second comparison matrix includes... , Indicator Factors With factors Compared to the importance of what is obtained, This is the proposed i-th alternative route; Wherein, the second comparison matrix N is a positive reciprocal matrix. The importance of satisfying the rules and The rules for determining their importance are the same; Step S6: Based on the second comparison matrix, obtain the second importance matrix of the proposed line layer relative to the i-th factor in the influencing index layer; the second importance matrix (a total of m matrices) includes: It should be noted that the elements at the same position within different matrices are different, such as b in M1. 1,2 and b in M2 1,2 They are not necessarily equal. i,j This represents any two elements combined randomly.

[0043] Step S7: Based on the second importance matrix, obtain each maximum eigenvalue and the corresponding second eigenvector, and normalize them to obtain the normalized second eigenvector. Based on the weights in the normalized second eigenvector, obtain the comprehensive weights of each factor in the proposed line layer for the optimal line layer. Specifically, the second feature vector The elements in the diagram represent the proposed circuit layer. n factors affect the indicator layer Factors in The weight is , ... ( ),in To influence the indicator layer The number of a certain factor; Step S8: Based on the comprehensive weight, determine the ranking of each candidate route in the proposed route layer.

[0044] Specifically, according to weight and weight The proposed circuit layer is obtained. The i-th factor For the optimal route Overall weighting: ; That is, the overall weight of each factor is: according to Sort by size to obtain the proposed circuit layer The sorting, The corresponding route is the optimal route.

[0045] Optionally, this embodiment of the invention also provides a specific embodiment of an intelligent expert judgment model as follows: To facilitate calculation and reduce the number of matrices and their elements, two factors were selected from each of the technical, geological, and economic indicators, for a total of six factors for calculation; the proposed routes were calculated as three routes.

[0046] Technical specifications include: planar alignment and burial depth; Geological indicators include: geological structural indicators and indicators of unfavorable geological bodies; Economic indicators include: project cost and demolition cost; Optimal path layer, using This indicates the selection of the final highway route; Influence indicator layer Including: planar linear ( ), burial depth ( ), geological structural indicators ( ), adverse geological indicators ( ), project cost ( ), demolition cost ( ); Proposed line layer ,include: , , ; First comparison matrix ,use It is represented as a positively reciprocal matrix. Indicates any two factors and Compared to the importance of what is obtained, that is The degree of importance follows these rules: When factors and When equally important, ; When factors Compare When it is slightly important, ; When factors Compare When it is obviously important, ; When factors Compare When it is extremely important, ; When factors Compare When extremely important and strong, ; When factors Compare Compared to when it is in the middle of the above adjacent importance levels, they are respectively , , , ; When factors With factors The importance of the comparison is At that time, then factors With factors The importance of the comparison is ; Based on impact indicators and The importance rules are used to determine the influence indicator layer. Relative to the optimal line layer Importance matrix: According to the matrix To obtain its largest eigenvalue and the corresponding normalized feature vector eigenvectors The elements in the layer are the influencing indicators. Relative to the optimal line layer Various factors Corresponding weights Based on the weights, the influencing indicator layer is obtained. Various factors Ranking by importance; Second comparison matrix ,use It is represented as a positively reciprocal matrix. Indicates any two factors With factors Compared to the importance of what is obtained, The importance of satisfying the rules and The rules for determining their importance are the same; Based on the proposed route, the proposed route layer is obtained. Compared to the influence indicator layer Importance matrix, including the proposed circuit layer Compared to the influence indicator layer Medium factors Importance matrix (6 matrices in total): , , , , ; Based on 6 matrices The largest eigenvalues ​​are obtained. and the corresponding normalized feature vector ,include: , ; , ; , ; , ; , ; , ; According to each matrix The corresponding normalized feature vector The total eigenvector is obtained. : Feature vector The elements in the diagram represent the proposed circuit layer. n factors affect the indicator layer Factors in The weight is , , ( ), To influence the indicator layer The number of a certain factor; According to weight and weight The proposed circuit layer is obtained. The i-th factor For the optimal route The weights are: ; That is, the weight of each factor is: according to Sort by size to obtain the proposed circuit layer The sorting, that is The corresponding route is the optimal route.

[0047] As described above, the embodiments of the present invention provide a highway route selection method based on integrated air-space-ground-indoor intelligent analysis. By integrating multi-source information such as InSAR data, airborne lidar data, geological surveys, geophysical tests, and drilling verification, a risk assessment model and an intelligent expert analysis model are constructed. This solves the problems of data bias and strong subjectivity when selecting routes based on a single factor. By comprehensively considering multiple influencing factors and comparing the importance of any two factors to calculate their weights, the scientificity and timeliness of highway route selection are improved.

[0048] Figure 6 This is a schematic diagram of a highway route selection system based on integrated air-space-ground-indoor intelligent analysis according to an embodiment of the present invention, applied to the highway route selection method based on integrated air-space-ground-indoor intelligent analysis provided in the embodiment of the present invention. Figure 6 As shown, the system includes: a preliminary partitioning module 10, a first screening module 20, a second screening module 30, a preliminary design module 40, a construction module 50, an evaluation module 60, an optimization module 70, and a judgment module 80.

[0049] Specifically, the preliminary zoning module 10 is used to perform preliminary zoning of the geological hazard susceptibility within the highway corridor based on the macro-route of the entire highway planned in the engineering feasibility study stage, and to perform preliminary marking based on historical geological hazard data to obtain a preliminary spatial distribution map of geological hazards along the highway. The first screening module 20 is used to perform preliminary identification and screening of potential geological hazards in the corridor zone of the preliminary spatial distribution map of geological hazards along the highway based on InSAR and remote sensing image interpretation technology, so as to obtain a large-scale spatial distribution map of potential geological hazards. The second screening module 30 is used to conduct a detailed investigation of the spatial distribution map of large-scale geological hazard hazards based on UAV-borne lidar and real-time dynamic positioning technology, and obtain a spatial distribution map of small-scale geological hazard hazards. Preliminary design module 40 is used to conduct preliminary route design based on a small-scale spatial distribution map of geological hazard risks, and to obtain the local route corridor range in the preliminary design stage; Module 50 is used to verify the surrounding geological hazards within the local corridor area based on comprehensive exploration methods and to build a geological condition information database; the comprehensive exploration methods include geological survey, geophysical testing and drilling verification; Assessment module 60 is used to conduct local geological risk assessment of potential geological hazards within the local corridor area based on a geological condition information database and using qualitative methods to obtain a risk assessment model; The optimization module 70 is used to make a preliminary judgment on the risk sources and risk zones along the local route corridor based on the risk assessment model, and to adjust and optimize the route within the local route corridor to obtain multiple routes to be planned. The analysis module 80 is used to construct an intelligent expert analysis model using a combination of qualitative and quantitative methods, and to analyze and compare multiple routes to be planned based on the intelligent expert analysis model to obtain the target planned route.

[0050] The above description is merely a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the technical scope disclosed in the present invention should be included within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims.

Claims

1. A highway route selection method based on space-earth-integrated wisdom judgment, characterized in that, The method comprises: Step S102, according to the highway macroscopic line planned in the engineering feasibility study stage, the geological disaster susceptibility in the highway corridor is preliminarily zoned, and the historical geological disaster data is preliminarily marked to obtain a preliminary geological disaster spatial distribution map along the highway; Step S104, based on InSAR and remote sensing image interpretation technology, the corridor geological disaster hidden danger in the preliminary geological disaster spatial distribution map along the highway is preliminarily identified and screened to obtain a large-scale geological disaster hidden danger spatial distribution map; Step S106, based on the unmanned aerial vehicle airborne laser radar and real-time dynamic positioning technology, the large-scale geological disaster hidden danger spatial distribution map is further investigated to obtain a small-scale geological disaster hidden danger spatial distribution map; Step S108, based on the small-scale geological disaster hidden danger spatial distribution map, a preliminary route design is carried out to obtain a local route corridor range in the preliminary design stage; Step S110, based on the comprehensive survey means, the surrounding geological disasters of the local route corridor range are checked, and a geological condition information database is constructed; the comprehensive survey means include geological investigation, geophysical prospecting test and drilling verification; Step S112, based on the geological condition information database, a qualitative method is used to carry out local geological risk assessment on the geological disaster hidden danger points in the local route corridor range to obtain a risk assessment model; Step S114, based on the risk assessment model, the on-line risk sources and risk partitions of the local route corridor range are preliminarily judged, and the route in the local route corridor range is adjusted and optimized to obtain a plurality of to-be-planned routes; Step S116, an intelligent expert judgment model is constructed by using a combination of qualitative and quantitative methods, and the plurality of to-be-planned routes are judged and compared based on the intelligent expert judgment model to obtain a target planning route.

2. The method of claim 1, wherein, According to the highway macroscopic line planned in the engineering feasibility study stage, the geological disaster susceptibility in the highway corridor is preliminarily zoned, and the historical geological disaster data is preliminarily marked to obtain a preliminary geological disaster spatial distribution map along the highway, comprising: According to the highway macroscopic line planned in the engineering feasibility study stage, the geological disaster susceptibility in the highway corridor is preliminarily zoned, and the historical geological disaster data is preliminarily marked to obtain a preliminary geological disaster spatial distribution map along the highway, comprising: According to the highway macroscopic line planned in the engineering feasibility study stage, the geological disaster susceptibility in the highway corridor is preliminarily zoned, and the historical geological disaster data is preliminarily marked to obtain a preliminary geological disaster spatial distribution map along the highway, comprising: According to the highway macroscopic line planned in the engineering feasibility study stage, the geological disaster susceptibility in the highway corridor is preliminarily zoned, and the historical geological disaster data is preliminarily marked to obtain a preliminary geological disaster spatial distribution map along the highway, comprising:

3. The method of claim 1, wherein, According to the highway macroscopic line planned in the engineering feasibility study stage, the geological disaster susceptibility in the highway corridor is preliminarily zoned, and the historical geological disaster data is preliminarily marked to obtain a preliminary geological disaster spatial distribution map along the highway, comprising: Based on the synthetic aperture radar interferometry technology, the interference effect between radar images of a target region acquired at different times is analyzed to identify the displacement of the ground surface between the two imaging times, and multi-scene SAR images of the same region at different time points are obtained; the target region is a corresponding region of a corridor in a preliminary geological disaster spatial distribution map along the highway; The multi-scene SAR images are respectively registered, interfered, filtered, and unwrapped to obtain a deformation time sequence of the target region in the entire time period, and based on the deformation time sequence, surface deformation monitoring data is obtained; Based on the surface deformation monitoring data, surface three-dimensional coordinate information of the target region before and after deformation is obtained, and based on the surface three-dimensional coordinate information of the target region before and after deformation, deformation characteristics, volume, and fall information of the adverse geological body are obtained; Based on the deformation characteristics of the adverse geological body, the type of the adverse geological body is identified; the type of the adverse geological body includes: landslide, debris flow, collapse, and dangerous rock mass; Based on the type of the adverse geological body and the volume and fall information, the influence range of the geological disaster is estimated; The influence range of the geological disaster is marked in the multi-scene SAR images to obtain a large-scale geological disaster hidden danger spatial distribution map.

4. The method of claim 1, wherein, Based on the unmanned aerial vehicle (UAV) airborne laser radar and real-time dynamic positioning technology, the large-scale geological disaster hidden danger spatial distribution map is further investigated to obtain a small-scale geological disaster hidden danger spatial distribution map, including: Based on the UAV airborne laser radar and UAV airborne camera photography aerial photography technology, the UAV flies close to the geological disaster hidden danger points in the large-scale geological disaster hidden danger spatial distribution map to obtain image data and three-dimensional point cloud data of the adverse geological body in the geological disaster hidden danger, and a three-dimensional model is generated by using the SFM three-dimensional image reconstruction technology; Based on the three-dimensional model, the occurrence information of the exposed structural surface of the adverse geological body in different geological disaster hidden dangers is obtained; The three-dimensional model is implanted into the image map of the large-scale geological disaster hidden danger spatial distribution map, and the occurrence information of the exposed structural surface of the adverse geological body is marked in the three-dimensional model to obtain a small-scale geological disaster hidden danger spatial distribution map.

5. The method of claim 1, wherein, Based on comprehensive investigation means, the surrounding geological disasters in the local line corridor range are checked to construct a geological condition information database, including: First geological information of the local line corridor range is obtained by using a geological survey method; the first geological information includes: stratum lithology, landform type, hydrology and meteorology, geological structure, vegetation coverage information, and information of human activity influence; Second geological information is obtained by using a geophysical prospecting test method, by testing the local line corridor range on the ground by using high-density electrical method, CSATM magnetotelluric method, micro surface wave method, and three-dimensional seismic method; Third geological information is obtained by using a drilling verification method, by carrying out internal drilling, comprehensive well logging, and deep displacement monitoring for the internal geology of the local line corridor range; Based on the first geological information, the second geological information, and the third geological information, a geological condition information database is constructed.

6. The method of claim 1, wherein, Based on the geological condition information base, a local geological risk assessment is performed on the geological disaster hidden danger points in the local line corridor range by using a qualitative method, and a risk assessment model is obtained, including: Based on the geological condition information base, a plurality of risk indicators are constructed, and each risk indicator is assigned a corresponding weight; the risk indicators include: topography and geomorphology indicators, geological structure indicators, stratum lithology indicators, hydrogeological condition indicators, and adverse geological body indicators; Based on the plurality of risk indicators, an expert group scoring is performed on the geological disaster hidden danger points in the local line corridor range, and a risk assessment model is constructed.

7. The method of claim 1, wherein, Based on the risk assessment model, a preliminary judgment is made on the line risk sources and risk partitions along the local line corridor range, and the line in the local line corridor range is adjusted and optimized to obtain a plurality of to-be-planned lines, including: Based on the risk assessment model, a preliminary judgment is made on the line risk sources and risk partitions along the local line corridor range, and the risk degrees of the plurality of risk partitions are obtained; Based on the risk degrees of each risk partition, the plurality of risk partitions are sorted; Based on the sorting of the plurality of risk partitions, the line in the local line corridor range is adjusted and optimized to obtain a plurality of to-be-planned lines; the optimization goal is to minimize the sum of the risk degrees of the risk partitions passed through.

8. The method of claim 1, wherein, An intelligent expert judgment model is constructed by using a combination of qualitative and quantitative methods, and based on the intelligent expert judgment model, the plurality of to-be-planned lines are judged and compared to obtain a target planning line, including: Based on the risk assessment model, a plurality of influence indicators are determined; the plurality of influence indicators include: a plurality of technical indicators, a plurality of geological indicators, and a plurality of economic indicators; Based on the plurality of influence indicators, a plurality of sensitive indicators are screened; Based on the plurality of sensitive indicators, an intelligent expert judgment model is constructed by using a combination of qualitative and quantitative methods; Based on the intelligent expert judgment model, the plurality of to-be-planned lines are judged and compared to obtain a target planning line.

9. The method of claim 8, wherein, An intelligent expert judgment model is constructed by using a combination of qualitative and quantitative methods, including: Based on the optimal route layer, the influence indicator layer, and the proposed route layer, an initial model structure is constructed; Establish a first comparison matrix for the influence index layer; the first comparison matrix includes , Indicator Factors and Compared to the importance of what is obtained, The first factor affecting highway alignment i One factor; Based on the first comparison matrix, a first importance degree matrix of the influence indicator layer relative to the optimal line layer O is constructed; the first importance degree matrix includes: According to the first importance degree matrix, the maximum eigenvalue and the corresponding first eigenvector are obtained, and the normalized first eigenvector is obtained by normalization, and according to the size of the elements in the normalized first eigenvector, the importance degree sorting of each factor in the influence indicator layer is obtained; each element in the normalized first eigenvector is the weight corresponding to each factor in the influence indicator layer relative to the optimal line layer; establishing a second contrast matrix of the tentative route layer; the second contrast matrix comprises , representing the importance degree obtained by comparing the factor with the factor , is the tentative first i alternative route; Based on the second comparison matrix, the proposed line layer is obtained relative to the influence index layer. i The second importance matrix of the factors; the second importance matrix includes: According to the second importance degree matrix, maximum eigenvalues and corresponding second eigenvectors are obtained, and normalization is performed to obtain normalized second eigenvectors, wherein each element of the normalized second eigenvectors is a weight of each factor in the proposed route layer corresponding to each factor in the optimal route layer, and according to the size of the elements in the normalized second eigenvectors, a comprehensive weight of each factor in the proposed route layer to the optimal route layer is obtained; the comprehensive weight comprises: , wherein n i is the comprehensive weight of the ith factor in the proposed route layer to the optimal route layer, n is the total number of factors in the proposed route layer, n i,j is an element in the normalized second eigenvectors, m j is an element in the normalized first eigenvectors. Based on the comprehensive weight, the sorting of each alternative route in the proposed route layer is determined.

10. A highway route selection system based on space-earth-integrated wisdom judgment, characterized in that, The application is applied to the highway route selection method based on space-earth-ground integrated judgment, and belongs to the technical field of highway route selection. The preliminary zoning module is configured to preliminarily zone the geological disaster proneness in a highway corridor based on a macro route of the whole highway planned in the engineering feasibility study stage, and preliminarily mark the geological disaster proneness based on historical geological disaster data to obtain a preliminary spatial distribution map of geological disasters along the highway. The first screening module is configured to preliminarily identify and screen the geological disaster hidden danger in the corridor based on the preliminary spatial distribution map of geological disasters along the highway by using InSAR and remote sensing image interpretation technology to obtain a large-scale spatial distribution map of geological disaster hidden dangers. The second screening module is configured to perform detailed investigation on the large-scale spatial distribution map of geological disaster hidden dangers by using an unmanned aerial vehicle (UAV) and real-time kinematic (RTK) technology to obtain a small-scale spatial distribution map of geological disaster hidden dangers. The preliminary design module is configured to perform preliminary route design based on the small-scale spatial distribution map of geological disaster hidden dangers to obtain a local route corridor range in the preliminary design stage. The construction module is configured to check the surrounding geological disasters of the local route corridor range by using comprehensive survey means, and construct a geological condition information database. The evaluation module is configured to perform local geological risk assessment on the geological disaster hidden danger points in the local route corridor range by using a qualitative method based on the geological condition information database to obtain a risk assessment model. The optimization module is configured to preliminarily judge the risk sources and risk zones along the local route corridor range based on the risk assessment model, and adjust and optimize the route in the local route corridor range to obtain a plurality of to-be-planned routes. The judgment module is configured to construct an intelligent expert judgment model by using a combination of qualitative and quantitative methods, and judge and compare the plurality of to-be-planned routes based on the intelligent expert judgment model to obtain a target planned route.

Images