Linear engineering geology three-dimensional zoning method, equipment and medium

Abstract

The invention discloses a linear engineering geology three-dimensional zoning method and device and a medium, and relates to the technical field of rock mass engineering geomechanics. The method comprises the steps that plane engineering geology zoning is conducted on a working area, and a preliminary plane engineering geology zoning result is obtained; optimizing the preliminary plane engineering geological zoning result to obtain a final plane engineering geological zoning result; performing section engineering geological zoning on the working area to obtain a preliminary section engineering geological zoning result; optimizing the preliminary section engineering geological zoning result to obtain a final section engineering geological zoning result; and combining the final plane engineering geology zoning result and the final section engineering geology zoning result to obtain a linear engineering geology three-dimensional zoning image. According to the method, engineering geology three-dimensional zoning can be carried out, and the efficiency of linear engineering geology zoning work and the precision of zoning results are improved.

Description

Technical Field

[0001] The present application relates to the technical field of rock mass engineering geomechanics, and in particular to a linear engineering geological three-dimensional zoning method, equipment and medium. Background Art

[0002] Engineering geological zoning is the process of dividing the work area into sections according to the needs of different stages of engineering construction and the principle of similar engineering geological conditions and similar engineering geological problems. In recent years, the planning, design and construction of linear projects such as railways and highways in my country that take into account the layout of surface and underground projects are gradually being carried out. Especially in complex geological environments, the safety and stability of engineering geology are crucial.

[0003] However, there are many existing engineering geological plane zoning indicators and profile zoning indicators for long and long linear projects in complex mountainous areas, and the zoning indicator data are difficult to obtain, which restricts the efficiency of long and long linear engineering geological zoning work in complex mountainous areas and the accuracy of zoning results. It is urgent to select key indicators for engineering geological plane zoning and profile zoning of long and long linear projects in complex mountainous areas and establish optimization indicator systems respectively. However, there is currently a lack of effective methods for optimizing key indicators of engineering geological zoning, and the existing technology does not involve the development of three-dimensional engineering geological zoning. Summary of the invention

[0004] The purpose of this application is to provide a linear engineering geological three-dimensional zoning method, equipment and medium, which can carry out engineering geological three-dimensional zoning, improve the efficiency of linear engineering geological zoning work and the accuracy of zoning results.

[0005] To achieve the above objectives, this application provides the following solutions.

[0006] In a first aspect, the present application provides a linear engineering geological three-dimensional zoning method, and the linear engineering geological three-dimensional zoning method includes the following steps.

[0007] The work area is subjected to plane engineering geological zoning to obtain preliminary plane engineering geological zoning results; the preliminary plane engineering geological zoning results are the results obtained based on plane zoning indicators.

[0008] The preliminary planar engineering geological zoning results are optimized to obtain final planar engineering geological zoning results.

[0009] Perform engineering geological zoning of the work area to obtain preliminary engineering geological zoning results; the preliminary engineering geological zoning results are the results obtained based on the profile zoning indicators.

[0010] The preliminary cross-section engineering geological zoning results are optimized to obtain final cross-section engineering geological zoning results.

[0011] The final planar engineering geological zoning result and the final cross-sectional engineering geological zoning result are combined to obtain a linear engineering geological three-dimensional zoning image.

[0012] In a second aspect, the present application provides a computer device, comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the above-mentioned linear engineering geological three-dimensional zoning method.

[0013] In a third aspect, the present application provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the above-mentioned linear engineering geological three-dimensional zoning method.

[0014] According to the specific embodiments provided in this application, this application discloses the following technical effects.

[0015] The present application provides a linear engineering geological three-dimensional zoning method, equipment and medium, the method comprising: performing plane engineering geological zoning on the working area to obtain preliminary plane engineering geological zoning results; the preliminary plane engineering geological zoning results are the results obtained based on the plane zoning index; optimizing the preliminary plane engineering geological zoning results to obtain the final plane engineering geological zoning results; performing profile engineering geological zoning on the working area to obtain preliminary profile engineering geological zoning results; the preliminary profile engineering geological zoning results are the results obtained based on the profile zoning index; optimizing the preliminary profile engineering geological zoning results to obtain the final profile engineering geological zoning results; combining the final plane engineering geological zoning results with the final profile engineering geological zoning results to obtain a linear engineering geological three-dimensional zoning image. The present application optimizes the plane engineering geological zoning results and the profile engineering geological zoning results, and fuses them to form a linear engineering geological three-dimensional zoning image, which can more intuitively understand the disaster-pregnant background of linear engineering geological problems and their control over engineering geological conditions, thereby scientifically guiding the prevention and control of linear engineering geological risks that take into account both surface and underground engineering layouts. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0017] Figure 1 This is an application environment diagram of a linear engineering geological three-dimensional zoning method in one embodiment of the present application.

[0018] Figure 2A schematic flow chart of a linear engineering geological three-dimensional zoning method provided in one embodiment of the present application.

[0019] Figure 3 A schematic diagram of the overall process of a linear engineering geological three-dimensional zoning method provided in one embodiment of the present application.

[0020] Figure 4 A schematic diagram of the structure of a computer device provided in one embodiment of the present application. DETAILED DESCRIPTION

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

[0022] In order to make the above-mentioned objects, features and advantages of the present application more obvious and easy to understand, the present application is further described in detail below in conjunction with the accompanying drawings and specific implementation methods.

[0023] Carrying out three-dimensional engineering geological zoning, that is, considering the zoning evaluation of surface engineering plane engineering geological conditions and the segmented evaluation of underground engineering profile engineering geological conditions at the same time, can provide an engineering geological basis for the planning, design and construction of linear projects that take into account the layout of both surface and underground projects, thereby effectively preventing and controlling potential geological disasters and ensuring the safety and smooth progress of project construction.

[0024] The linear engineering geological three-dimensional zoning method provided in the embodiment of the present application can be applied to Figure 1In the application environment shown. Among them, the terminal 102 communicates with the server 104 through the network. The data storage system can store the data that the server 104 needs to process. The data storage system can be set up separately, integrated on the server 104, or placed on the cloud or other servers. The terminal 102 can send the work area to the server 104. After the server 104 receives the work area, for the work area, the server 104 performs a plane engineering geological zoning on the work area to obtain a preliminary plane engineering geological zoning result; the preliminary plane engineering geological zoning result is a result obtained based on the plane zoning index; the preliminary plane engineering geological zoning result is optimized to obtain the final plane engineering geological zoning result; the work area is subjected to a profile engineering geological zoning to obtain a preliminary profile engineering geological zoning result; the preliminary profile engineering geological zoning result is a result obtained based on the profile zoning index; the preliminary profile engineering geological zoning result is optimized to obtain the final profile engineering geological zoning result; the final plane engineering geological zoning result and the final profile engineering geological zoning result are combined to obtain a linear engineering geological stereoscopic zoning image. The server 104 can feed back the obtained linear engineering geological three-dimensional zoning image to the terminal 102. In addition, in some embodiments, the linear engineering geological three-dimensional zoning method can also be implemented by the server 104 or the terminal 102 alone, such as the terminal 102 can directly perform linear engineering geological three-dimensional zoning on the working area, or the server 104 can obtain the working area from the data storage system and perform linear engineering geological three-dimensional zoning on the working area.

[0025] The terminal 102 may be, but is not limited to, various desktop computers, laptop computers, smart phones, and tablet computers. The server 104 may be implemented as an independent server or a server cluster consisting of multiple servers, or may be a cloud server.

[0026] In an exemplary embodiment, Figure 2 and Figure 3 As shown, a linear engineering geological three-dimensional zoning method is provided. The method is executed by a computer device, and can be executed by a computer device such as a terminal or a server alone, or can be executed by a terminal and a server together. In the embodiment of the present application, the method is applied to Figure 1 The server 104 in the example is used for explanation, and the steps include the following steps S1 to S5.

[0027] S1: Performing plane engineering geological zoning on the work area to obtain preliminary plane engineering geological zoning results; the preliminary plane engineering geological zoning results are the results obtained based on plane zoning indicators.

[0028] S2: Optimizing the preliminary planar engineering geological zoning results to obtain final planar engineering geological zoning results.

[0029] S3: Performing cross-sectional engineering geological zoning on the work area to obtain preliminary cross-sectional engineering geological zoning results; the preliminary cross-sectional engineering geological zoning results are results obtained based on cross-sectional zoning indicators.

[0030] S4: Optimizing the preliminary cross-section engineering geological zoning results to obtain final cross-section engineering geological zoning results.

[0031] S5: combining the final planar engineering geological zoning result and the final cross-sectional engineering geological zoning result to obtain a linear engineering geological three-dimensional zoning image.

[0032] By implementing the above steps S1 to S5, three-dimensional engineering geological zoning can be carried out, thereby improving the efficiency of linear engineering geological zoning work and the accuracy of zoning results.

[0033] As an optional implementation, in step S1, the following contents are specifically included.

[0034] S11: Acquire engineering geological condition data of the work area; the engineering geological condition data of the work area includes: geological pictures, remote sensing images and geological survey information.

[0035] S12: Determine the plane zoning index according to the engineering geological condition data of the work area.

[0036] S13: Divide the working area into grids to obtain a plurality of first grid units of the same size; each of the first grid units has a value of a corresponding plane zoning index.

[0037] S14: performing standardization processing on each of the first grid units to obtain a standardized first unit.

[0038] S15: Based on the Pearson correlation coefficient in the SPSS software, determine the first correlation between the plane division indicators in the first unit after the standardization process.

[0039] S16: taking the plane division indicators whose first correlation is greater than the first preset threshold as preliminary screening plane division indicators.

[0040] S17: According to the engineering geological comprehensive zoning model, the comprehensive weight of each indicator in the preliminary screening plane zoning indicators is calculated; the engineering geological comprehensive zoning model is a model that introduces Euclidean distance to establish a hierarchical analysis method-entropy value method coupling.

[0041] S18: Based on the comprehensive weights of each indicator in the preliminary screening plane zoning indicators, the first engineering geological zoning index is calculated using the geographic information system spatial analysis method and the multi-indicator superposition method; the first engineering geological zoning index divides the engineering geological conditions into four levels: good, medium, poor and bad.

[0042] S19: According to the first engineering geological zoning index, a preliminary planar engineering geological zoning result is obtained by adopting natural breakpoint classification or expert experience method.

[0043] After obtaining the preliminary planar engineering geological zoning results, it is necessary to optimize the preliminary planar engineering geological zoning results. Therefore, in step S2, the following contents are specifically included.

[0044] S21: Randomly arrange and combine the preliminary screened plane zoning indicators and use them as input layers respectively; use the preliminary plane engineering geological zoning results as output layer results; and perform training based on the DNN model to obtain the first optimal indicator combination.

[0045] S22: Using the first optimal indicator combination as the input layer, respectively, to perform DNN simulation prediction results to obtain the final plane engineering geological zoning results.

[0046] In an exemplary embodiment, step 1: performing planar engineering geological zoning of the work area includes the following contents.

[0047] A1: Construct a planar engineering geological zoning index system, establish an engineering geological zoning model, and use the geographic information system to conduct planar engineering geological zoning of the work area.

[0048] A2: Based on the results of the planar engineering geological zoning in step A1, model training is performed with the help of machine learning to screen out a combination of key zoning indicators with general adaptability and high accuracy.

[0049] A3: Optimize the plane zoning results based on the combination of key indicators and machine learning.

[0050] The specific steps include the following.

[0051] Step 1.1: Preliminarily determine the plane zoning indicators and collect indicator data. With the help of ArcGIS software, the data of each indicator is initially processed, and the work area is gridded to obtain several grid units of the same size, that is, each grid unit has the value of each indicator. The grid within the plane work area is converted into point data based on the grid conversion tool. The data of each zoning indicator raster file is extracted into point data using the multi-value extraction to point tool. Excel tables or txt text data can be exported for easy processing and analysis. The data of each indicator needs to be standardized, and the correlation between the indicators is analyzed based on the Pearson correlation coefficient in SPSS software. When the absolute value of the correlation coefficient between two factors is greater than 0.5, it is considered to have a high correlation, and the zoning indicators can be preliminarily screened.

[0052] The hierarchical method is adopted to divide the planar engineering geological zoning index system into five levels: overall target layer, element layer (first-level indicator layer), indicator layer (second-level indicator layer), indicator classification layer and zoning object layer. The target layer is the target of the system analysis and the engineering geological zoning of the work area; the element layer is the influencing factor of the target layer, which can consider structural stability and surface stability; the indicator layer is the influencing factor of the element layer, which is the various factors affecting the regional crustal stability; the indicator classification layer is the numericalization of the abstract concept of the element layer, which is divided and assigned based on the weight calculation of the above element layer; the zoning object layer is the result of the engineering geological zoning of the target layer area.

[0053] The Euclidean distance was introduced to establish a comprehensive engineering geological zoning model that coupled the analytic hierarchy process and the entropy value method, and the comprehensive weights of each indicator were obtained. According to the expert scores, the judgment matrices of the zoning element layer and the indicator layer were constructed respectively. The analytic hierarchy process weights of each indicator were calculated by formula (1), and the relative weights of each factor in the judgment matrix were calculated and consistency tests were performed.

[0054] (1); in, is the AHP weight of each indicator, i.e., the eigenvector (weight coefficient), is the number of factors, is the factor matrix row, is the factor matrix column, For indicators With indicators A quantitative expression of the relative importance or priority between For indicators With indicators A quantitative expression of the relative importance or priority between them.

[0055] The objective weights of factor indicators are calculated using the entropy method. Based on MATLAB software, the decision matrix is ​​standardized and objective weights are assigned to the indicator factors, as shown in formula (2). The entropy value and difference value of the indicator are calculated and the entropy method weights of each indicator are obtained.

[0056] (2); in, is the entropy weight of each indicator, is the number of factors, is the sample number, is the indicator number, This is a sample dataset of indicators for the research area.

[0057] The Euclidean distance is introduced to couple the weights of each indicator obtained by the hierarchical analysis method and the entropy method to obtain the combined weighted value, as shown in formula (3).

[0058] (3); in, is the comprehensive weight of each indicator, α is the preference degree of the AHP weight of each indicator, that is, the preference degree of the subjective weight, β is the preference degree of the entropy weight of each indicator, that is, the preference degree of the objective weight, is the analytic hierarchy process weight of each indicator, is the entropy weight of each indicator.

[0059] Determine the difference degree equation by establishing the Euclidean distance function α and β .

[0060] (4); (5); (6); The engineering geological zoning index (Formula 7) is obtained by using the spatial analysis method of the geographic information system and the multi-index superposition method. The index classification is assigned based on the reclassification tool. The raster calculator tool is used to multiply each index layer with the comprehensive weight and superimpose the result layer. The calculation results of the plane engineering geological zoning of the work area are obtained by using methods such as natural breakpoint classification or expert experience. Based on the comprehensive analysis of the calculation results and the field geological survey data, the engineering geological conditions are divided into four levels: good, medium, poor, and poor engineering geological conditions. The typical engineering geological characteristics and major engineering geological problems in different zones are explained respectively, the possibility and spatial location of major engineering geological problems are evaluated and predicted, and suggestions for preventing engineering geological problems or comparing engineering routes are proposed.

[0061] (7); in, is the engineering geological zoning index, is the coupling weight of AHP-entropy method, Is the indicator index.

[0062] Step 1.2: Based on the torch library of Python language, construct a DNN model for training. Randomly arrange and combine the initial screened plane zoning indicators, and gradually increase the number of combinations from 1 to all indicators, as the input layer. Use the four levels of engineering geological conditions of the work area zoning results obtained in step A1, namely good, medium, poor, and bad, to assign values ​​1, 2, 3, and 4 respectively as the output layer results. The number of neurons is set to 100, 1 hidden layer, Relu is used as the activation function, the cross entropy loss function is used, the learning rate is set to 0.01, the ratio of training set to test set is 7:3, and the DNN model is trained.

[0063] According to the results of DNN model training, the accuracy, weighted recall rate, and F1 value (comprehensive accuracy and recall rate) of each indicator combination of plane zoning are obtained. By sorting and analyzing the above results from high to low, we can obtain an indicator combination with universal adaptability and high accuracy.

[0064] Step 1.3: According to the optimal indicator combination obtained in step 1.2, each combination is used as the input layer to perform DNN simulation prediction results. The parameters of the model remain unchanged and the simulation results are used as the final planar engineering geological zoning results of the work area.

[0065] As an optional implementation, in step S3, the following contents are specifically included.

[0066] S31: Acquire engineering geological condition data of the work area; the engineering geological condition data of the work area includes: geological pictures, remote sensing images and geological survey information.

[0067] S32: Determine profile zoning indicators based on the engineering geological conditions data of the work area.

[0068] S33: Divide the working area into grids to obtain a plurality of second grid units of the same size; each of the second grid units has a value corresponding to the profile zoning index.

[0069] S34: performing standardization processing on each of the second grid units to obtain a standardized second unit.

[0070] S35: Based on the Pearson correlation coefficient in the SPSS software, determine the second correlation between the profile zoning indicators in the second unit after the standardization process.

[0071] S36: Use the plane zoning index whose second correlation is greater than the second preset threshold as the preliminary screening section zoning index.

[0072] S37: Using the analytic hierarchy process, determine the weight of each indicator in the preliminary screening profile zoning indicator.

[0073] S38: According to the weight of each indicator in the preliminary screening profile zoning indicators, the second engineering geological zoning index is calculated by using the geographic information system spatial analysis method and the multi-indicator superposition method; the second engineering geological zoning index divides the engineering geological conditions into four levels: good, medium, poor and bad.

[0074] S39: According to the second engineering geological zoning index, a natural breakpoint classification or expert experience method is used to obtain a profile engineering geological zoning result.

[0075] Similarly, after obtaining the profile engineering geological zoning results, the profile engineering geological zoning results need to be optimized, so in step S4, the following contents are specifically included.

[0076] S41: Randomly arrange and combine the preliminary screening profile zoning indicators and use them as input layers respectively; the preliminary profile engineering geological zoning results are used as output layer results; and the second optimal indicator combination is obtained by training based on the DNN model.

[0077] S42: Using the second optimal indicator combination as the input layer, respectively, to perform DNN simulation prediction results to obtain the final profile engineering geological zoning results.

[0078] In an exemplary embodiment, step 2: performing cross-sectional engineering geological zoning of the work area includes the following contents.

[0079] B1: Construct a profile engineering geological zoning index system, establish an engineering geological zoning model, and use the geographic information system to conduct profile engineering geological zoning of the work area.

[0080] B2: Based on the results of the cross-section engineering geological zoning in step B1, model training is performed with the help of machine learning to screen out a combination of key zoning indicators with universal adaptability and high accuracy.

[0081] B3: Optimize the profile zoning results based on the combination of key indicators and machine learning.

[0082] The specific steps include the following.

[0083] Step 2.1: Based on the plane zoning results obtained in step A1, the profile zoning results are based on the engineering design elevation to obtain the segmented evaluation map of the engineering geological conditions of the recommended route profile. Preliminary determination of the profile zoning indicators, and collection of indicator data to construct an indicator system for the engineering geological zoning of key sections. The target layer is the engineering geological zoning of the profile; the element layer can consider factors such as geological structure, stratum lithology and rock and soil characteristics, adverse geological effects, and ground temperature; the indicator layer is the various factors that affect the engineering geological conditions; the indicator classification layer is the numericalization of the abstract concept of the element layer, which is divided and assigned based on the weight calculation of the above element layer; the zoning object layer is the result of the engineering geological zoning of the key section profile.

[0084] With the help of ArcGIS software, the data of each indicator of the profile are preliminarily processed. Similarly, the work area is gridded. The indicator grids within the work area are extracted into point data using the multi-value extraction to point tool. Excel tables or txt text data can be exported. The data of each indicator are standardized, and the correlation between the indicators is analyzed based on the Pearson correlation coefficient in SPSS software. When the absolute value of the correlation coefficient between two factors is greater than 0.5, it is considered to have a high correlation, and the profile zoning indicators can be preliminarily screened.

[0085] The weights of the engineering geological zoning indicators of the profile are determined by the analytic hierarchy process. According to the expert scores, the judgment matrices of the zoning element layer and the indicator layer are constructed respectively, and the weights of the profile zoning indicators are obtained according to the calculation. The engineering geological zoning index is obtained by means of the spatial analysis method of the geographic information system and the multi-indicator superposition method. The indicator classification is assigned based on the reclassification tool. The raster calculator tool is used to multiply each indicator layer with the comprehensive weight and superimpose them to obtain the result layer. The calculation results of the engineering geological zoning of the work area profile are obtained by using methods such as natural breakpoint classification or expert experience. Based on the comprehensive analysis of the calculation results and the field geological survey data, the engineering geological conditions are zoned, and the engineering geological conditions are divided into four levels: good, medium, poor, and poor.

[0086] Step 2.2: Based on the torch library of Python language, construct a DNN model for training. Randomly arrange and combine the zoning indicators of the initially screened profiles, and gradually increase the number of combinations from 1 to all indicators, as the input layer. Use the four levels of engineering geological conditions of the work area zoning results obtained in step 1, namely good, medium, poor, and bad, to assign values ​​1, 2, 3, and 4 respectively as the output layer results. The number of neurons is set to 100, 1 hidden layer, Relu is used as the activation function, the cross entropy loss function is used, the learning rate is set to 0.01, the ratio of training set to test set is 7:3, and the DNN model is trained.

[0087] According to the results of DNN model training, the accuracy, weighted recall rate, and F1 value (comprehensive accuracy and recall rate) of each group of indicators in the profile zoning are obtained. The above results are sorted and analyzed from high to low to obtain an indicator combination with universal adaptability and high accuracy.

[0088] Step 2.3: According to the optimal indicator combination obtained in step 2.2, each combination is used as the input layer to perform DNN simulation prediction results. The parameters of the model remain unchanged, and the simulation results are used as the final engineering geological zoning results of the working area.

[0089] As an optional implementation, in step S5, the following contents are specifically included.

[0090] S51: Smoothing the final planar engineering geological zoning result and the final cross-sectional engineering geological zoning result to obtain a smoothed planar engineering geological zoning result and a smoothed cross-sectional engineering geological zoning result.

[0091] S52: Based on ArcGIS software, the smoothed profile engineering geological zoning results are extracted in the form of points, the points are matched one by one with the smoothed plane engineering geological zoning results, and Kriging interpolation is used to generate the engineering geological zoning results of the profile railway line passing through the area and convert them into plane expression.

[0092] S53: Determine the engineering geological conditions expressed in a plane according to the smoothed section engineering geological zoning result, and obtain a planarized section zoning result.

[0093] S54: Using the mosaic to new grid tool, the planarized section zoning result is superimposed on the smoothed planar engineering geological zoning result to obtain a planar section combined zoning result of a key section.

[0094] S55: Based on the focus statistics tool in ArcGIS software, the horizontal section combined with the zoning results of the key section are smoothed and optimized to obtain a linear engineering geological three-dimensional zoning image.

[0095] In an exemplary embodiment, step 3: combining plan and section engineering geological zoning results, includes the following contents.

[0096] C1: Optimized plane and section engineering geological zoning results obtained after smoothing steps 1 and 2 C2: Combine the optimized plan and section zoning results in step 1; C3: Optimize and smooth the delineation result image.

[0097] The specific steps include the following.

[0098] Step 3.1: Smooth the optimized plane and section engineering geological zoning results obtained in step 1 and step 2. Based on the focus statistics tool in ArcGIS software, the statistical type is MAJORITY, and the plane and section results are smoothed respectively.

[0099] Step 3.2: Based on the range of not less than 1 km on both sides of the engineering line to be further compared and optimized, the smoothed section engineering geological zoning results are extracted in the form of points based on ArcGIS software, and the points are matched one by one with the plane. The engineering geological zoning results of the section railway line through the area are converted into plane expression by Kriging interpolation. According to the engineering geological zoning results of the section, the plane expression is divided into good, medium, poor and bad engineering geological conditions. The mosaic to new raster tool is used, and the mosaic operator selects MAXIMUM, that is, the principle of "higher rather than lower". The planarized section zoning results are superimposed on the plane zoning results to obtain the combined planar and section zoning results of the key section.

[0100] Step 3.3: The grid data of the three-dimensional zoning result map obtained in step 3.2 can be smoothed. Based on the focus statistics tool in ArcGIS software, the statistical type is MAJORITY, and the result is smoothed and optimized to obtain the final linear engineering geological three-dimensional zoning image.

[0101] This application proposes a method for optimizing linear engineering geological zoning indicators and optimizing engineering geological three-dimensional zoning results based on geographic information system and machine learning. The key indicators of engineering geological plane zoning and profile zoning of long and large linear projects in complex mountainous areas are selected, and optimization indicator systems are established respectively. Detailed steps are given for integrating the engineering geological plane zoning results and profile segmentation results to form the engineering geological three-dimensional zoning results and optimizing them. This method can more intuitively understand the disaster-prone background of linear engineering geological problems and their control over engineering geological conditions, thereby scientifically guiding the prevention and control of linear engineering geological risks that take into account both surface and underground engineering layouts.

[0102] The present application also provides an application scenario, which applies the above-mentioned linear engineering geological three-dimensional zoning method. Specifically: the linear engineering geological three-dimensional zoning method provided in this embodiment can be applied in the linear engineering geological three-dimensional zoning scenario. The linear engineering geological three-dimensional zoning scene includes: a plane engineering geological zoning link, a plane engineering geological zoning result optimization link, a profile engineering geological zoning link, a profile engineering geological zoning result optimization link and a combination link; first, a plane engineering geological zoning is performed on the working area to obtain a preliminary plane engineering geological zoning result; the preliminary plane engineering geological zoning result is a result obtained based on the plane zoning index; secondly, the preliminary plane engineering geological zoning result is optimized to obtain a final plane engineering geological zoning result; then, a profile engineering geological zoning is performed on the working area to obtain a preliminary profile engineering geological zoning result; the preliminary profile engineering geological zoning result is a result obtained based on the profile zoning index; thirdly, the preliminary profile engineering geological zoning result is optimized to obtain a final profile engineering geological zoning result; finally, the final plane engineering geological zoning result and the final profile engineering geological zoning result are combined to obtain a linear engineering geological three-dimensional zoning image.

[0103] In an exemplary embodiment, a computer device is provided. The computer device may be a server or a terminal. The internal structure diagram thereof may be as follows: Figure 4 As shown. The computer device includes a processor, a memory, an input / output interface (Input / Output, referred to as I / O) and a communication interface. The processor, the memory and the input / output interface are connected through a system bus, and the communication interface is connected to the system bus through the input / output interface. The processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The database of the computer device is used to store a workspace. The input / output interface of the computer device is used to exchange information between the processor and an external device. The communication interface of the computer device is used to communicate with an external terminal through a network connection. When the computer program is executed by the processor, a linear engineering geological three-dimensional zoning method is implemented.

[0104] Those skilled in the art will understand that Figure 4 The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components.

[0105] In an exemplary embodiment, a computer device is further provided, including a memory and a processor, wherein a computer program is stored in the memory, and the above-mentioned method embodiments are implemented when the processor executes the computer program.

[0106] In an exemplary embodiment, a computer-readable storage medium is provided, storing a computer program, and when the computer program is executed by a processor, the above-mentioned method embodiments are implemented.

[0107] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, stored data, displayed data, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of relevant data must comply with relevant regulations.

[0108] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to the memory, database or other medium used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetoresistive random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM may be in various forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM).

[0109] The database involved in each embodiment provided in this application may include at least one of a relational database and a non-relational database. The non-relational database may include a distributed database based on blockchain, etc., but is not limited thereto. The processor involved in each embodiment provided in this application may be a general-purpose processor, a central processing unit, a graphics processor, a digital signal processor, a programmable logic device, a data processing logic device based on quantum computing, etc., but is not limited thereto.

[0110] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0111] This article uses specific examples to illustrate the principles and implementation methods of this application. The description of the above embodiments is only used to help understand the method and core ideas of this application. At the same time, for those skilled in the art, according to the ideas of this application, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as limiting this application.

Claims

1. A linear engineering geological three-dimensional zoning method, characterized in that: The linear engineering geological three-dimensional zoning method includes: Performing plane engineering geological zoning on the work area to obtain preliminary plane engineering geological zoning results; the preliminary plane engineering geological zoning results are the results obtained based on the plane zoning indicators; Optimizing the preliminary planar engineering geological zoning results to obtain final planar engineering geological zoning results; Performing cross-sectional engineering geological zoning on the work area to obtain preliminary cross-sectional engineering geological zoning results; the preliminary cross-sectional engineering geological zoning results are results obtained based on cross-sectional zoning indicators; Optimizing the preliminary engineering geological zoning results of the profile to obtain final engineering geological zoning results of the profile; The final planar engineering geological zoning result and the final cross-sectional engineering geological zoning result are combined to obtain a linear engineering geological three-dimensional zoning image.

2. The linear engineering geological three-dimensional zoning method according to claim 1 is characterized in that: Carry out plane engineering geological zoning of the work area and obtain preliminary plane engineering geological zoning results, including: Obtaining engineering geological condition data of the work area; the engineering geological condition data of the work area includes: geological pictures, remote sensing images and geological survey information; Determine the plane zoning index according to the engineering geological conditions of the work area; Divide the working area into grids to obtain a plurality of first grid units of the same size; each of the first grid units has a value of a corresponding plane zoning index; Performing a standardization process on each of the first grid units to obtain a first unit after the standardization process; Based on the Pearson correlation coefficient in SPSS software, determining the first correlation between the plane division indicators in the first unit after the standardization process; Taking the plane division index whose first correlation is greater than the first preset threshold as the preliminary screening plane division index; According to the engineering geological comprehensive zoning model, the comprehensive weight of each index in the preliminary screening plane zoning index is calculated; the engineering geological comprehensive zoning model is a model that introduces Euclidean distance to establish a hierarchical analysis method-entropy value method coupling; According to the comprehensive weight of each index in the preliminary screening plane zoning index, the first engineering geological zoning index is calculated by using the geographic information system spatial analysis method and the multi-index superposition method; the first engineering geological zoning index divides the engineering geological conditions into four levels: good, medium, poor and bad; According to the first engineering geological zoning index, a preliminary planar engineering geological zoning result is obtained by using natural breakpoint classification or expert experience method.

3. The linear engineering geological three-dimensional zoning method according to claim 2 is characterized in that: The preliminary planar engineering geological zoning results are optimized to obtain the final planar engineering geological zoning results, which specifically include: The preliminary screening plane zoning indicators are randomly arranged and combined as input layers respectively; the preliminary plane engineering geological zoning results are used as output layer results; training is performed based on the DNN model to obtain the first optimal indicator combination; The first optimal indicator combination is used as the input layer to perform DNN simulation prediction results to obtain the final plane engineering geological zoning results.

4. The linear engineering geological three-dimensional zoning method according to claim 2 is characterized in that: The calculation formula for the comprehensive weight of each indicator is: ； in, is the comprehensive weight of each indicator, α is the preference degree of the AHP weight of each indicator, β is the preference degree of the entropy weight of each indicator, is the analytic hierarchy process weight of each indicator, is the entropy weight of each indicator.

5. The linear engineering geological three-dimensional zoning method according to claim 2 is characterized in that: The calculation formula of the engineering geological zoning index is: ； in, is the engineering geological zoning index, is the coupling weight of AHP-entropy method, Is the indicator index.

6. The linear engineering geological three-dimensional zoning method according to claim 1, characterized in that: Carry out cross-section engineering geological zoning of the work area and obtain preliminary cross-section engineering geological zoning results, including: Obtaining engineering geological condition data of the work area; the engineering geological condition data of the work area includes: geological pictures, remote sensing images and geological survey information; Determine the section zoning index according to the engineering geological conditions of the work area; Divide the working area into grids to obtain a plurality of second grid units of the same size; each of the second grid units has a value corresponding to the profile zoning index; Performing a standardization process on each of the second grid units to obtain a standardized second unit; Based on the Pearson correlation coefficient in SPSS software, the second correlation between the profile zoning indicators in the second unit after the standardization process is determined; Taking the plane division index whose second correlation is greater than the second preset threshold as the preliminary screening section division index; Using the analytic hierarchy process, the weight of each index in the preliminary screening profile zoning index is determined; According to the weight of each index in the preliminary screening profile zoning index, the second engineering geological zoning index is calculated by using the geographic information system spatial analysis method and the multi-index superposition method; the second engineering geological zoning index divides the engineering geological conditions into four levels: good, medium, poor and bad; According to the second engineering geological zoning index, the profile engineering geological zoning result is obtained by using natural breakpoint classification or expert experience method.

7. The linear engineering geological three-dimensional zoning method according to claim 6 is characterized in that: The preliminary cross-section engineering geological zoning results are optimized to obtain the final cross-section engineering geological zoning results, which specifically include: The preliminary screening profile zoning indicators are randomly arranged and combined as input layers; the preliminary profile engineering geological zoning results are used as output layer results; training is performed based on the DNN model to obtain the second optimal indicator combination; The second optimal indicator combination is used as the input layer to perform DNN simulation prediction results to obtain the final profile engineering geological zoning results.

8. The linear engineering geological three-dimensional zoning method according to claim 1 is characterized in that: The final planar engineering geological zoning result and the final cross-sectional engineering geological zoning result are combined to obtain a linear engineering geological three-dimensional zoning image, which specifically includes: Smoothing the final planar engineering geological zoning result and the final cross-sectional engineering geological zoning result to obtain a smoothed planar engineering geological zoning result and a smoothed cross-sectional engineering geological zoning result; Based on ArcGIS software, the smoothed cross-sectional engineering geological zoning results are extracted in the form of points, the points are matched one by one with the smoothed plane engineering geological zoning results, and the engineering geological zoning results of the cross-sectional railway line passing through the area are converted into plane expression by using Kriging interpolation; Determining the engineering geological conditions expressed in a plane according to the smoothed section engineering geological zoning result to obtain a planarized section zoning result; Using the mosaic to new grid tool, the planarized section zoning result is superimposed with the smoothed planar engineering geological zoning result to obtain a planar section combined zoning result of a key section; Based on the focus statistics tool in ArcGIS software, the horizontal section combined with the zoning results of the key section are smoothed and optimized to obtain a linear engineering geological three-dimensional zoning image.

9. A computer device comprising: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the linear engineering geological three-dimensional zoning method according to any one of claims 1 to 8.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the linear engineering geological three-dimensional zoning method described in any one of claims 1 to 8 is implemented.

Images