A linear engineering geological three-dimensional zoning method, device and medium

Abstract

The present application discloses a method, device and medium for linear engineering geological three-dimensional zoning, which relates to the technical field of rock mass engineering geology. The method includes: conducting plane engineering geological zoning on the working area to obtain a preliminary plane engineering geological zoning result; optimizing the preliminary plane engineering geological zoning result to obtain a final plane engineering geological zoning result; conducting profile engineering geological zoning on the working area to obtain a preliminary profile engineering geological zoning result; optimizing the preliminary profile engineering geological zoning result to obtain a final profile engineering geological zoning result; combining the final plane engineering geological zoning result and the final profile engineering geological zoning result to obtain a linear engineering geological three-dimensional zoning image. The present application can carry out engineering geological three-dimensional zoning, improving the efficiency of linear engineering geological zoning work and the accuracy of the zoning result.

Description

Technical Field

[0001] The present application relates to the technical field of rock mass engineering geology, and particularly to a method, device and medium for linear engineering geological three-dimensional zoning. Background Art

[0002] Engineering geological zoning is a process of dividing and segmenting the working area according to the needs of different stages of engineering construction, based on 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 China, which take into account both surface and underground engineering layouts, have been gradually carried out. Especially in complex geological environments, the engineering geological safety and stability are of crucial importance.

[0003] However, there are many engineering geological plane zoning indicators and profile zoning indicators for existing long and large linear projects in complex mountainous areas, and it is difficult to obtain the zoning indicator data, which restricts the efficiency of the engineering geological zoning work and the accuracy of the zoning results for long and large linear projects in complex mountainous areas. It is urgent to select the key indicators for the engineering geological plane zoning and profile zoning of long and large linear projects in complex mountainous areas and establish optimized indicator systems respectively. However, there is currently a lack of an effective method for selecting key indicators for engineering geological zoning, and the prior art does not involve carrying out engineering geological three-dimensional zoning. Summary of the Invention

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

[0005] To achieve the above purpose, the present application provides the following solutions.

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

[0007] Carry out plane engineering geological zoning on the working area to obtain a preliminary plane engineering geological zoning result; the preliminary plane engineering geological zoning result is a result based on plane zoning indicators.

[0008] Optimize the preliminary plane engineering geological zoning result to obtain a final plane engineering geological zoning result.

[0009] Carry out profile engineering geological zoning on the working area to obtain a preliminary profile engineering geological zoning result; the preliminary profile engineering geological zoning result is a result based on profile zoning indicators.

[0010] Optimize the preliminary profile engineering geological zoning result to obtain a final profile engineering geological zoning result.

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

[0012] In a second aspect, the present application provides a computer device, including: a memory, a processor, and a computer program stored on the memory and executable on the processor, where 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, on which a computer program is stored, and when the computer program is executed by a processor, it implements the above-mentioned linear engineering geological three-dimensional zoning method.

[0014] According to the specific embodiments provided by the present application, the following technical effects are disclosed.

[0015] The present application provides a linear engineering geological three-dimensional zoning method, device, and medium. The method includes: performing planar engineering geological zoning on a working area to obtain a preliminary planar engineering geological zoning result; the preliminary planar engineering geological zoning result is a result based on planar zoning indicators; optimizing the preliminary planar engineering geological zoning result to obtain a final planar engineering geological zoning result; performing sectional engineering geological zoning on the working area to obtain a preliminary sectional engineering geological zoning result; the preliminary sectional engineering geological zoning result is a result based on sectional zoning indicators; optimizing the preliminary sectional engineering geological zoning result to obtain a final sectional engineering geological zoning result; combining the final planar engineering geological zoning result and the final sectional engineering geological zoning result to obtain a linear engineering geological three-dimensional zoning image. After optimizing the planar engineering geological zoning result and the sectional engineering geological zoning result, the present application integrates them to form a linear engineering geological three-dimensional zoning image, which can more intuitively understand the disaster-forming background of linear engineering geological problems and their control effects on engineering geological conditions, thereby scientifically guiding the prevention and control of linear engineering geological risks that take into account the layout of surface and underground projects. Description of the Drawings

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

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

[0018] Figure 2Schematic flowchart of a linear engineering geological three-dimensional zoning method provided by an embodiment of the present application.

[0019] Figure 3 Overall flowchart of a linear engineering geological three-dimensional zoning method provided by an embodiment of the present application.

[0020] Figure 4 Schematic structural diagram of a computer device provided by an embodiment of the present application. Detailed implementation manners

[0021] Next, the technical solutions in the embodiments of the present application will be clearly and completely described in conjunction with the accompanying drawings in the embodiments of the present application. Obviously, the described embodiments are only a part of the embodiments of the present application, rather than all the embodiments. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present application without creative efforts shall fall within the protection scope of the present application.

[0022] To make the above objects, features, and advantages of the present application more obvious and understandable, the present application will be further described in detail below with reference to the accompanying drawings and specific implementation manners.

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

[0024] The linear engineering geological three-dimensional zoning method provided by the embodiments of the present application can be applied, for example, 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 separately, integrated on the server 104, 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 conducts a planar engineering geological zoning on the work area to obtain a preliminary planar engineering geological zoning result; the preliminary planar engineering geological zoning result is a result obtained based on planar zoning indicators; optimize the preliminary planar engineering geological zoning result to obtain a final planar engineering geological zoning result; conduct a sectional engineering geological zoning on the work area to obtain a preliminary sectional engineering geological zoning result; the preliminary sectional engineering geological zoning result is a result obtained based on sectional zoning indicators; optimize the preliminary sectional engineering geological zoning result to obtain a final sectional engineering geological zoning result; combine the final planar engineering geological zoning result and the final sectional engineering geological zoning result to obtain a linear engineering geological three-dimensional zoning image. The server 104 can feedback 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 independently by the server 104 or the terminal 102. For example, the terminal 102 can directly conduct a linear engineering geological three-dimensional zoning on the work area, or the server 104 can obtain the work area from the data storage system and conduct a linear engineering geological three-dimensional zoning on the work area.

[0025] Among them, the terminal 102 can be, but is not limited to, various desktop computers, laptop computers, smart phones and tablet computers. The server 104 can be implemented by an independent server or a server cluster composed of multiple servers, and can also be a cloud server.

[0026] In an exemplary embodiment, as Figure 2 and Figure 3 shown, a linear engineering geological three-dimensional zoning method is provided. This method is executed by a computer device, and can specifically be executed independently by a computer device such as a terminal or a server, or jointly executed by a terminal and a server. In the embodiments of the present application, taking this method applied to Figure 1 the server 104 in as an example for illustration, it includes the following steps S1 to S5.

[0027] S1: Conduct a planar engineering geological zoning on the work area to obtain a preliminary planar engineering geological zoning result; the preliminary planar engineering geological zoning result is a result obtained based on planar zoning indicators.

[0028] S2: Optimize the preliminary planar engineering geological zoning result to obtain a final planar engineering geological zoning result.

[0029] S3: Conduct profile engineering geological zoning for the working area to obtain the preliminary profile engineering geological zoning result; the preliminary profile engineering geological zoning result is the result obtained based on the profile zoning index.

[0030] S4: Optimize the preliminary profile engineering geological zoning result to obtain the final profile engineering geological zoning result.

[0031] S5: Combine the final plane engineering geological zoning result and the final profile engineering geological zoning result to obtain a linear engineering geological three-dimensional zoning image.

[0032] Implementing the above steps S1 to S5 can carry out engineering geological three-dimensional zoning, improving the efficiency of linear engineering geological zoning work and the accuracy of the zoning result.

[0033] As an alternative implementation method, in step S1, it specifically includes the following contents.

[0034] S11: Obtain the engineering geological condition data of the working area; the engineering geological condition data of the working area includes: geological pictures, remote sensing images, and geological exploration information.

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

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

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

[0038] S15: Based on the Pearson correlation coefficient in SPSS software, determine the first correlation between the plane zoning indices in the standardized first unit.

[0039] S16: Use the plane zoning indices with the first correlation greater than the first preset threshold as the preliminary screened plane zoning indices.

[0040] S17: According to the engineering geological comprehensive zoning model, calculate the comprehensive weights of the indices in the preliminary screened plane zoning indices; the engineering geological comprehensive zoning model is a model that couples the analytic hierarchy process and the entropy method by introducing the Euclidean distance.

[0041] S18: According to the comprehensive weights of each index in the preliminary screening of plane zoning indicators, using the spatial analysis method of geographic information system and the multi-index superposition method, calculate the first engineering geological zoning index; the first engineering geological zoning index divides the engineering geological conditions into 4 grades: good, medium, poor, and very poor.

[0042] S19: According to the first engineering geological zoning index, use the natural breakpoint grading or expert experience method to obtain the preliminary plane engineering geological zoning result.

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

[0044] S21: Randomly permute and combine the preliminary screening plane zoning indicators and use them as the input layer respectively; the preliminary plane engineering geological zoning result is used as the output layer result; train based on the DNN model to obtain the first optimal index combination.

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

[0046] In an exemplary embodiment, step 1: Conduct plane engineering geological zoning on the working area, including the following contents.

[0047] A1: Construct a plane engineering geological zoning index system, establish an engineering geological zoning model, and conduct plane engineering geological zoning on the working area with the help of geographic information system.

[0048] A2: Based on the results of the plane engineering geological zoning in step A1, use machine learning to train the model and screen out the key index combinations with universal adaptability and high accuracy.

[0049] A3: Optimize the plane zoning result based on the key index combination and machine learning.

[0050] The specific steps include the following contents.

[0051] Step 1.1: Initially determine the plane zoning indicators and collect the indicator data. With the help of the geographic information system ArcGIS software, initially process the data of each indicator, divide the working area into grids, and obtain several grid units of the same size, that is, each grid unit has the values of each indicator. Convert the raster within the plane working area into point data based on the raster to point tool, and use the multi-value extraction to point tool to extract the data in the zoning indicator raster file into the point data. An excel table or txt text data can be exported for convenient processing and analysis. It is necessary to standardize the data of each indicator, and analyze the correlation between each indicator based on the Pearson correlation coefficient in the SPSS software. When the absolute value of the correlation coefficient between two factors is greater than 0.5, it is considered that they have a high correlation, and the zoning indicators can be initially screened.

[0052] The plane engineering geological zoning indicator system is divided into 5 levels: the overall goal layer, the element layer (the first-level indicator layer), the indicator layer (the second-level indicator layer), the indicator grading layer, and the zoning object layer by using the hierarchical step-by-step method. The goal layer is the goal of the system analysis, which is the engineering geological zoning of the working area; the element layer is the influencing factor of the goal layer, and tectonic stability and surface stability can be considered; the indicator layer is the influencing factor of the element layer, which are the various factors affecting regional crustal stability; the indicator grading 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 regional engineering geological zoning of the goal layer.

[0053] Introduce the Euclidean distance to establish an engineering geological comprehensive zoning model coupling the analytic hierarchy process and the entropy method to obtain the comprehensive weight of each indicator. According to the expert scoring, construct the judgment matrices of the zoning element layer and the indicator layer respectively, and calculate the weights of each indicator by the analytic hierarchy process using formula (1), calculate the relative weights of each factor in the judgment matrix and conduct a consistency test.

[0054] (1);

[0055] Among them, is the weight of each indicator by the analytic hierarchy process, that is, the eigenvector (weight coefficient), is the number of factors, is the row of the factor matrix, is the column of the factor matrix, is the indicator and the indicator is the quantitative representation of the relative importance or priority between them, is the indicator and the indicator is the quantitative representation of the relative importance or priority between them.

[0056] The objective weights of factor indicators are calculated using the entropy value method. Based on MATLAB software, the decision matrix is standardized, and objective weight assignment is performed on the indicator factors. As shown in Equation (2), the entropy value and difference value of the indicators are calculated to obtain the entropy value method weights of each indicator.

[0057] (2);

[0058] Among them, is the entropy value method weight of each indicator, is the number of factors, is the sample number, is the indicator number, is the indicator sample data set of the research area.

[0059] The Euclidean distance is introduced to couple the weights of each indicator obtained by the analytic hierarchy process and the entropy value method to obtain the combined weight value, as shown in Equation (3).

[0060] (3);

[0061] Among them, is the comprehensive weight of each indicator, α is the preference degree of the analytic hierarchy process weight of each indicator, that is, the preference degree of the subjective weight, β is the preference degree of the entropy value method 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 value method weight of each indicator.

[0062] According to the Euclidean distance function, a difference degree equation is established to determine α and β .

[0063] (4);

[0064] (5);

[0065] (6);

[0066] The engineering geological zoning index (Formula 7) is obtained by means of the spatial analysis method of geographic information system and the multi-index superposition method. The index grading is assigned values based on the reclassification tool, and each index layer is multiplied by the comprehensive weight respectively using the raster calculator tool and then superimposed to obtain the result layer. Methods such as natural breakpoint grading or expert experience are used to obtain the calculation results of the plane engineering geological zoning in the working area. On the basis of comprehensively analyzing the calculation results and the field geological survey data, the engineering geological conditions are partitioned, and accordingly divided into 4 grades: good, medium, poor, and very poor engineering geological conditions. The typical engineering geological characteristics and main engineering geological problems in different partitions are described respectively, the possibility and spatial location of major engineering geological problems are evaluated and predicted, and suggestions for preventing engineering geological problems or suggestions for route selection of engineering lines are put forward.

[0067] (7);

[0068] Among them, is the engineering geological zoning index, is the coupling weight of analytic hierarchy process-entropy method, is the index index.

[0069] Step 1.2: Based on the torch library of Python language, a DNN model is constructed for training. The preliminary screened zoning indicators in each area of the plane are randomly arranged and combined, and the combination quantity gradually increases from 1 to all indicators, which are used as the input layer respectively. The 4 grades of good, medium, poor, and very poor engineering geological conditions of the zoning results of the working area obtained in Step A1 are assigned 1, 2, 3, and 4 respectively as the output layer results. The number of neurons is set to 100, with 1 hidden layer. The activation function uses Relu, the cross-entropy loss function is adopted, the learning rate is set to 0.01, and the ratio of the training set to the test set is 7:3 for DNN model training.

[0070] According to the results obtained from the DNN model training, the accuracy, weighted recall rate, and F1 value (comprehensive accuracy and recall rate) of each index combination of the plane zoning are obtained. The above results are sorted and analyzed from high to low, and an index combination with universal adaptability and high accuracy can be obtained.

[0071] Step 1.3: According to the optimal index combination obtained in Step 1.2, it is used as the input layer respectively to perform the DNN simulation prediction results. The parameters of the model remain unchanged, and the simulation results obtained are used as the final plane engineering geological zoning results of the working area.

[0072] As an optional implementation manner, in Step S3, it specifically includes the following contents.

[0073] S31: Obtain the engineering geological condition data of the working area; the engineering geological condition data of the working area includes: geological pictures, remote sensing images, and geological exploration information.

[0074] S32: Determine the profile zoning index according to the engineering geological condition data of the work area.

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

[0076] S34: Perform standardization processing on each of the second grid units to obtain the standardized second unit.

[0077] S35: Based on the Pearson correlation coefficient in SPSS software, determine the second correlation between the profile zoning indexes in the standardized second unit.

[0078] S36: Take the plane zoning index with the second correlation greater than the second preset threshold as the preliminary screening profile zoning index.

[0079] S37: Use the analytic hierarchy process to determine the weight of each index in the preliminary screening profile zoning index.

[0080] S38: According to the weight of each index in the preliminary screening profile zoning index, use the geographic information system spatial analysis method and the multi-index superposition method to calculate and obtain the second engineering geological zoning index; the second engineering geological zoning index divides the engineering geological conditions into 4 grades: good, medium, poor, and very poor.

[0081] S39: According to the second engineering geological zoning index, use the natural breakpoint grading or expert experience method to obtain the profile engineering geological zoning result.

[0082] Similarly, after obtaining the profile engineering geological zoning result, it is necessary to optimize the profile engineering geological zoning result. Therefore, in step S4, the following specific contents are included.

[0083] S41: Randomly arrange and combine the preliminary screening profile zoning indexes as the input layer respectively; the preliminary profile engineering geological zoning result is used as the output layer result; train based on the DNN model to obtain the second optimal index combination.

[0084] S42: Use the second optimal index combination as the input layer respectively to perform DNN simulation prediction results to obtain the final profile engineering geological zoning result.

[0085] In an exemplary embodiment, step 2: Perform profile engineering geological zoning on the work area, including the following contents.

[0086] B1: Construct an index system for profile engineering geological zoning, establish an engineering geological zoning model, and conduct profile engineering geological zoning for the working area with the aid of a geographic information system.

[0087] B2: Based on the results of the profile engineering geological zoning in step B1, conduct model training with the aid of machine learning to screen out a combination of key zoning indicators with universal adaptability and high accuracy.

[0088] B3: Optimize the profile zoning results based on the combination of key indicators with the aid of machine learning.

[0089] The specific steps include the following.

[0090] Step 2.1: Based on the planar zoning results obtained in step A1, with the engineering design elevation as the baseline for the profile zoning results, obtain a sectional evaluation map of the engineering geological conditions of the recommended line profile. Initially determine the profile zoning indicators, collect indicator data, and construct an index system for the profile engineering geological zoning of key sections. The target layer is the profile engineering geological zoning; the factor layer can consider elements such as geological structure, stratigraphic lithology and geotechnical properties, adverse geological effects, and ground temperature; the indicator layer is each factor affecting the engineering geological conditions; the indicator classification layer is the numericalization of the abstract concepts of the factor layer, which is divided and assigned based on the calculation of the weights of the above factor layers; the zoning object layer is the result of the profile engineering geological zoning of key sections.

[0091] With the aid of the geographic information system ArcGIS software, conduct preliminary processing of the profile indicator data. Similarly, divide the working area into grids, extract the indicator rasters within the working area to point data using the multi-value extraction to point tool, and export the data to an excel table or txt text. Standardize each indicator data, analyze the correlation between indicators 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 that they have a high correlation, and the profile zoning indicators can be initially screened.

[0092] Use the analytic hierarchy process to determine the weights of the profile engineering geological zoning indicators. According to expert scoring, construct judgment matrices for the zoning factor layer and indicator layer respectively. Through calculation, obtain the weights of the profile zoning indicators. With the aid of the spatial analysis method and multi-index overlay method of the geographic information system, obtain the engineering geological zoning index, assign values to the indicator classification based on the reclassification tool, and use the raster calculator tool to multiply each indicator layer by the comprehensive weight and overlay them to obtain the result layer. Use methods such as natural breakpoint classification or expert experience to obtain the calculation results of the profile engineering geological zoning of the working area. Based on the comprehensive analysis of the calculation results and field geological survey data, conduct engineering geological condition zoning, and correspondingly divide it into 4 grades: good, medium, poor, and very poor engineering geological conditions.

[0093] Step 2.2: Based on the torch library of the Python language, construct a DNN model for training. Randomly permute and combine the initially screened zoning indicators for each section. The number of combinations gradually increases from 1 to all indicators, which are used as the input layer respectively. The engineering geological conditions of the zoning results of the working area obtained in Step 1, namely good, medium, poor, and very poor, are assigned values of 1, 2, 3, and 4 respectively as the output layer results. The number of neurons is set to 100, with 1 hidden layer. The activation function uses Relu, the cross-entropy loss function is adopted, the learning rate is set to 0.01, and the ratio of the training set to the test set is 7:3 for DNN model training.

[0094] According to the results obtained from the DNN model training, calculate the accuracy, weighted recall rate, and F1 value (combining accuracy and recall rate) for each set of indicator combinations of the section zoning. Sort and analyze the above results from high to low to obtain indicator combinations with universal adaptability and high accuracy.

[0095] Step 2.3: Use the optimal indicator combinations obtained in Step 2.2 as the input layer respectively to perform DNN simulation predictions. Keep the parameters of the model unchanged, and use the obtained simulation results as the final section engineering geological zoning results of the working area.

[0096] As an alternative implementation, in Step S5, it specifically includes the following content.

[0097] S51: Smooth the final planar engineering geological zoning results and the final section engineering geological zoning results to obtain the smoothed planar engineering geological zoning results and the smoothed section engineering geological zoning results.

[0098] S52: Based on the ArcGIS software, extract the smoothed section engineering geological zoning results in the form of points, correspond each point to the smoothed planar engineering geological zoning results, and use Kriging interpolation to generate the engineering geological zoning results of the area through which the section railway line passes and convert them into a planar expression.

[0099] S53: Determine the engineering geological conditions of the planar expression according to the smoothed section engineering geological zoning results to obtain the planarized section zoning results.

[0100] S54: Use the Mosaic to New Raster tool to overlay the planarized section zoning results with the smoothed planar engineering geological zoning results to obtain the combined planar and section zoning results for key sections.

[0101] S55: Based on the Focal Statistics tool in the ArcGIS software, perform smoothing and optimization processing on the combined planar and section zoning results for key sections to obtain a linear engineering geological three-dimensional zoning image.

[0102] In an exemplary embodiment, Step 3: Combine the plane and profile engineering geological zoning results, including the following content.

[0103] C1: Smooth the optimized plane and profile engineering geological zoning results obtained in Steps 1 and 2

[0104] C2: Combine the optimized plane and profile zoning results in Step 1;

[0105] C3: Optimize and smooth the zoning result image.

[0106] The specific steps include the following content.

[0107] Step 3.1: Smooth the optimized plane and profile engineering geological zoning results obtained in Steps 1 and 2. Based on the Focal Statistics tool in ArcGIS software, with the statistical type being MAJORITY, smooth the plane and profile results respectively.

[0108] Step 3.2: According to the range on both sides of the engineering line to be further compared and optimized, with a width of not less than 1 km on each side, based on ArcGIS software, extract the smoothed profile engineering geological zoning results in the form of points, correspond the points to the plane one by one, use Kriging interpolation to generate the engineering geological zoning results of the area where the profile railway line passes and convert them into a plane expression. According to the engineering geological zoning results of the profile, divide the plane expression into good, medium, poor, and very poor engineering geological conditions, and use the Mosaic to New Raster tool, select the MAXIMUM mosaic operator, that is, the principle of "taking the higher value rather than the lower value", and overlay the planarized profile zoning results with the plane zoning results to obtain the combined plane and profile zoning results of the key section.

[0109] Step 3.3: The raster data of the three-dimensional zoning result map obtained in Step 3.2 can be smoothed. Based on the Focal Statistics tool in ArcGIS software, with the statistical type being MAJORITY, perform smoothing and optimization processing on the results to obtain the final linear engineering geological three-dimensional zoning image.

[0110] This application proposes a method for optimizing the selection of linear engineering geological zoning indicators and the results of engineering geological three-dimensional zoning based on geographic information systems and machine learning. It selects the key indicators for engineering geological plane zoning and profile zoning of long linear engineering in complex mountainous areas and establishes optimization index systems respectively. It gives the detailed steps for fusing the engineering geological plane zoning results and profile sectioning results to form the engineering geological three-dimensional zoning results and optimizing them, which can more intuitively understand the disaster-forming background of linear engineering geological problems and their control effects on engineering geological conditions, so as to scientifically guide the prevention and control of linear engineering geological risks that take into account the layout of surface and underground projects.

[0111] 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 to the linear engineering geological three-dimensional zoning scenario. The linear engineering geological three-dimensional zoning scenario includes: a planar engineering geological zoning link, an optimization link for the planar engineering geological zoning result, a sectional engineering geological zoning link, an optimization link for the sectional engineering geological zoning result, and a combination link; First, conduct a planar engineering geological zoning on the working area to obtain a preliminary planar engineering geological zoning result; The preliminary planar engineering geological zoning result is the result obtained based on the planar zoning index; Secondly, optimize the preliminary planar engineering geological zoning result to obtain the final planar engineering geological zoning result; Then, conduct a sectional engineering geological zoning on the working area to obtain a preliminary sectional engineering geological zoning result; The preliminary sectional engineering geological zoning result is the result obtained based on the sectional zoning index; Thirdly, optimize the preliminary sectional engineering geological zoning result to obtain the final sectional engineering geological zoning result; Finally, combine the final planar engineering geological zoning result and the final sectional engineering geological zoning result to obtain a linear engineering geological three-dimensional zoning image.

[0112] In an exemplary embodiment, a computer device is provided. The computer device can be a server or a terminal, and its internal structure diagram can be as Figure 4 shown. The computer device includes a processor, a memory, an input / output interface (Input / Output, abbreviated as I / O), and a communication interface. Among them, 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. Among them, 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 the working area. The input / output interface of the computer device is used to exchange information between the processor and external devices. The communication interface of the computer device is used to communicate with external terminals through a network connection. When the computer program is executed by the processor, it realizes a linear engineering geological three-dimensional zoning method.

[0113] Those skilled in the art can understand that Figure 4 the structure shown in

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

[0115] 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 method embodiments are implemented.

[0116] 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 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 need to comply with relevant regulations.

[0117] Those of ordinary skill in the art can understand that all or part of the processes of implementing the above method embodiments can be completed by instructing relevant hardware through a computer program. 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 above method embodiments. Among them, any reference to a memory, database, or other medium used in the embodiments provided in this application can include at least one of non-volatile and volatile memories. Non-volatile memories 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 memories can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM), etc.

[0118] In each of the embodiments provided in the present application, the database involved 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., without limitation. In each of the embodiments provided in the present application, the processor involved may be a general-purpose processor, a central processing unit, a graphics processing unit, a digital signal processor, a programmable logic device, a data processing logic device based on quantum computing, etc., without limitation.

[0119] The technical features of the above embodiments can be combined arbitrarily. For the sake of concise description, 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, it should be considered as the scope recorded in this specification.

[0120] Specific examples are used in this article to elaborate on the principles and implementation manners of the present application. The description of the above embodiments is only used to help understand the method and its core idea of the present application; at the same time, for those of ordinary skill in the art, according to the idea of the present application, there will be changes in the specific implementation manners and application scopes. In summary, the content of this specification should not be construed as a limitation to the present 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; Combining the final planar engineering geological zoning result with the final cross-sectional engineering geological zoning result to obtain a linear engineering geological three-dimensional zoning image; The preliminary planar engineering geological zoning results are optimized to obtain the final planar engineering geological zoning results, which specifically include: The preliminary screened 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 index combination is used as the input layer to perform DNN simulation prediction results to obtain the final plane engineering geological zoning results; 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 respectively; 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 index combination is used as the input layer to perform DNN simulation prediction results to obtain the final cross-section engineering geological zoning results; 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.

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 coupled analytic hierarchy process and entropy method; 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 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.

4. 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.

5. The linear engineering geological three-dimensional zoning method according to claim 1 is 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.

6. 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 5.

7. 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 5 is implemented.

Images