Attribute information-oriented tunnel three-dimensional geological modeling and updating method and system

Abstract

The present application relates to railway tunnel construction and three-dimensional geological modeling technical field, provide the tunnel three-dimensional geological modeling and update method and system for attribute information, including: obtaining profile data and surface elevation data, after establishing tunnel structure model, attribute model is established with structure model constraint; the attribute model is divided into several grids, and attribute assignment is carried out, and attribute grid is obtained; in response to incremental data, attribute grid is constructed without attribute grid, and attribute filling is carried out according to the proportion of attribute grid in attributeless grid according to the principle of weighted filling. Real-time control of the engineering geological conditions of the surrounding rock of the tunnel is realized, and the three-dimensional geological model construction of the whole process of tunnel construction and excavation is suitable for various complex geological conditions.

Description

Technical Field

[0001] This invention belongs to the field of railway tunnel construction and three-dimensional geological modeling technology, and in particular relates to a method and system for three-dimensional geological modeling and updating of tunnels based on attribute information. Background Technology

[0002] The statements in this section are merely background information related to the present invention and do not necessarily constitute prior art.

[0003] More and more deep-buried railway tunnels are being built in mountainous areas with complex geological conditions. During the design and construction of railway tunnels, it is necessary to understand the relevant properties of the surrounding rock to provide parameters for tunnel support design. Geological objects possess information in three aspects: geometry, topology, and attributes. In existing 3D geological modeling processes, most focus only on constructing the geometric structure model of 3D geology while neglecting the attribute model of 3D geology. This results in a model that looks good but has little practical use and is not very practical.

[0004] Moreover, existing three-dimensional geological models are all static models, while tunnel construction is a dynamic process. As the tunnel advances, new geological information is constantly added, and this information has not been well utilized.

[0005] In summary, most existing 3D geological modeling is static modeling, and when building models, they often consider the structural information of geological bodies while ignoring the numerous attribute information of the geological bodies being modeled, resulting in models with weak applicability. Summary of the Invention

[0006] To address the technical problems mentioned above, this invention provides a method and system for three-dimensional geological modeling and updating of tunnels based on attribute information. This method enables real-time control of the engineering geological conditions of the surrounding rock of the tunnel and is applicable to the construction of three-dimensional geological models throughout the entire tunnel construction and excavation process under various complex geological conditions.

[0007] To achieve the above objectives, the present invention adopts the following technical solution: The first aspect of the present invention provides a method for three-dimensional geological modeling and updating of tunnels based on attribute information, comprising: After acquiring profile data and surface elevation data, and establishing a tunnel structure model, an attribute model constrained by the structure model is established. The attribute model is divided into several grids, and attribute values ​​are assigned to obtain the attribute grid; In response to incremental data, an attribute-free grid is constructed, and attribute filling is performed according to the proportion of attribute grids to attribute-free grids, following the principle of weighted filling.

[0008] Furthermore, after regularizing the profile data and surface elevation data, a tunnel structure model is established.

[0009] Furthermore, the profile data includes borehole data and profile data.

[0010] Furthermore, based on the borehole data, the profile data is converted into borehole data and virtual borehole data to establish a tunnel structure model.

[0011] Furthermore, the attribute model is divided into several regular hexahedral meshes, and the attribute value of each mesh is assigned to the body center of the regular hexahedron.

[0012] Furthermore, the attribute assignment employs three-dimensional spatial attribute interpolation.

[0013] Furthermore, the three-dimensional spatial attribute interpolation employs radial basis functions, Kriging and its variants, inverse distance weighting, spline interpolation, or Thiessen polygon interpolation methods.

[0014] A second aspect of the present invention provides a tunnel 3D geological modeling and updating system oriented towards attribute information, comprising: The model building module is configured to: acquire profile data and surface elevation data, build a tunnel structure model, and then build an attribute model constrained by the structure model. The model partitioning module is configured to: divide the attribute model into several grids, assign attribute values, and obtain the attribute grid; The model update module is configured to: in response to incremental data, construct an attribute-free grid and fill it with attributes according to the proportion of attribute grids to attribute-free grids, following the principle of weighted filling.

[0015] A third aspect of the present invention provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps in the attribute-information-oriented tunnel three-dimensional geological modeling and updating method described above.

[0016] A fourth aspect of the present invention provides a computer device including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the steps in the attribute-information-oriented tunnel three-dimensional geological modeling and updating method described above.

[0017] Compared with the prior art, the beneficial effects of the present invention are: This invention provides a method for three-dimensional geological modeling and updating of tunnels based on attribute information. By assigning values ​​to the structural model and adding incremental data to update the model, it achieves real-time control of the engineering geological conditions of the surrounding rock of the tunnel, especially the parameters related to support design. It is applicable to the construction of three-dimensional geological models for the entire process of tunnel construction and excavation under various complex geological conditions.

[0018] This invention provides a method for three-dimensional geological modeling and updating of tunnels based on attribute information. After determining the attribute values ​​of each newly created block and completing the attribute filling, the overall newly created block model will be smoother than the previous attribute data model, and the transition between two different attributes will be more gradual. Attached Figure Description

[0019] The accompanying drawings, which form part of this invention, are used to provide a further understanding of the invention. The illustrative embodiments of the invention and their descriptions are used to explain the invention and do not constitute an improper limitation of the invention.

[0020] Figure 1 This is a flowchart of the tunnel 3D geological modeling and updating method oriented towards attribute information according to Embodiment 1 of the present invention; Figure 2 This is a diagram showing the numerical attribute update and filling situation during model update in case A of Embodiment 1 of the present invention; Figure 3 This is a diagram showing the numerical attribute update and filling situation during model update in case B of Embodiment 1 of the present invention; Figure 4 This is a diagram showing the numerical attribute update and filling situation during model update in case C of Embodiment 1 of the present invention; Figure 5 This is a diagram showing the numerical attribute update and filling situation during model update in case D of Embodiment 1 of the present invention. Detailed Implementation

[0021] The present invention will be further described below with reference to the accompanying drawings and embodiments.

[0022] It should be noted that the following detailed description is illustrative and intended to provide further explanation of the invention. Unless otherwise specified, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention pertains.

[0023] Example 1 This embodiment provides a method for three-dimensional geological modeling and updating of tunnels based on attribute information.

[0024] The attribute-information-oriented tunnel 3D geological modeling and updating method provided in this embodiment focuses on the geological structure, rock type, stratum thickness, water-bearing capacity, wave velocity, uniaxial compressive strength, and shear strength surrounding the tunnel. c / φ The description and expression of parameters such as ) enable unified modeling of the geometric objects, topological relationships, and attribute information of geological bodies, providing parameter basis for subsequent tunnel support design.

[0025] The tunnel 3D geological modeling and updating method provided in this embodiment uses incremental data to dynamically update the established model in real time.

[0026] The tunnel 3D geological modeling and updating method based on attribute information provided in this embodiment, such as Figure 1 As shown, the process includes defining the modeling scope, modeling the tunnel structure, constraining the structure model, dividing the tunnel mesh, assigning tunnel attributes, updating the model, and visualizing the model.

[0027] (1) Determine the scope of tunnel modeling.

[0028] (101) The top surface of the model is based on the ground elevation DEM data. The top surface reflects the surface undulation, the boundary lines of different lithologies and their topological relationships. (102) The bottom surface of the model is a plane at a certain elevation that is a certain distance below the tunnel track surface elevation. (103) The left and right sides of the model are parallel to the tunnel direction and are a certain distance away from the tunnel. (104) The model length needs to include the tunnel entrance and exit. This length is not the tunnel mileage length, but should be greater than the tunnel mileage length.

[0029] In this embodiment, the determined modeling range is a large area reflecting the entire area around the tunnel, which is a regional model.

[0030] (2) Tunnel structure modeling.

[0031] (201) Data regularization processing: Regularize the profile data (including borehole data and profile data) and surface elevation (DEM) data; (202) Based on borehole data, other profile data are converted into borehole data and virtual borehole data to establish a tunnel structure model. The profile includes stratigraphic boundaries and other boundaries. To convert it into virtual borehole data, it is only necessary to set a series of vertical straight lines at a certain sampling distance, extract the coordinates of the intersection of the straight lines and the profile boundaries, and extract the lithological information between two adjacent coordinates.

[0032] (203) If the data conflict with each other during modeling, the data from different sources should be prioritized before modeling, with the borehole data having the highest priority.

[0033] These borehole control points, distributed across various strata, are subjected to Delaunay triangulation to form a borehole stratigraphic control triangular grid covering the entire modeling area. This is used to determine stratigraphic lithological boundaries.

[0034] In this embodiment, the tunnel geometry modeling yields a tunnel region model.

[0035] (3) Structural model constraints.

[0036] The property distribution around the tunnel is based on the structural model, and a property model constrained by the structural model is constructed. The structural constraints of the property model are generated according to the top and bottom positions of each geological unit, and then the property model volume mesh is constructed. The structural constraints of the property model refer to the framework of the model, mainly including the top and bottom surfaces of the model, as well as the boundaries of the lithology of each stratum.

[0037] (4) Tunnel mesh generation. Divide the attribute model into several regular hexahedral meshes.

[0038] In the horizontal direction, a regular grid is divided according to the modeling accuracy requirements. The size of each grid represents the length and width of the volume element, which determines the modeling accuracy in the horizontal direction. Generally, a regular hexahedral grid is used to construct the attribute model grid, with the center of the regular hexahedron as the calculation point.

[0039] (5) Assigning tunnel attributes. Using the body center of the regular hexahedron as the calculation point, three-dimensional spatial attribute interpolation is performed to obtain the attribute mesh.

[0040] The body center of the regular hexahedral mesh and its corresponding stratigraphic unit are determined. The number of boreholes involved in the interpolation calculation and their associated attributes are also determined. Attributes are assigned using three-dimensional spatial attribute interpolation. Methods for three-dimensional spatial attribute interpolation include radial basis functions, Kleekin and its variants, inverse distance weighting, spline interpolation, and Thiessen polygon interpolation. The series of attributes includes rock type, water abundance, wave velocity, uniaxial compressive strength, and shear strength. c / φ ) and other parameters.

[0041] In this embodiment, assigning attributes to the mesh is assigning attributes to the rock mass surrounding the tunnel wall, which is an outcrop-scale attribute model.

[0042] (6) Model update. In response to incremental data, construct an attribute-free grid and fill it with attributes according to the proportion of attribute grids to attribute-free grids, following the principle of weighted filling.

[0043] The model is updated and corrected using incremental data obtained from tunnel excavation and exploration. The incremental data includes information on lithology and stratum thickness revealed by excavation, as well as TSP and induced polarization data for advanced geological prediction.

[0044] After the update, the following information will be displayed: geological structure, rock type, stratum thickness, water-bearing capacity, wave velocity, uniaxial compressive strength, and shear strength. c / φ Detailed information on parameters such as )

[0045] When filling the newly constructed layer of attribute-free mesh blocks with attributes, a weighted filling principle will be followed, that is, filling according to the proportion of each established regular cube to the attribute blocks. For ease of analysis, the three-dimensional spatial problem will be transformed into a two-dimensional planar spatial problem for a detailed explanation of the attribute filling principle. In the two-dimensional plane, four cases A, B, C, and D will occur. Specifically, when updating numerical attributes, the filling process will be divided into these four cases: A, B, C, and D. Scenario A: such as Figure 2 The newly created block is completely enclosed by the attribute grid. Therefore, it will be determined that the newly created block will fully display the attributes within the attribute grid. Figure 2 In this context, assuming the attribute value of the attribute grid is "5", the attribute of the newly created block will also be represented as "5". Case B: such as Figure 3 As shown, the newly created block appears in two attribute grids. Assuming both attribute grids have the same attribute and a value of "5", then the newly created block will also have an attribute of "5". If... Figure 3 As shown, the attribute grid has attributes "5" and "6". Therefore, it is necessary to consider using a weighted method to calculate the attribute value of this newly created block. The weighted formula is as follows:

[0046] In the formula, Y n This represents the attribute value of the target newly created block; A This represents the attribute value of the first attribute grid occupied by the newly created block. N 1 represents the weight of the newly created block in the first attribute grid; B This represents the attribute value of the second attribute grid occupied by the newly created block. N 2 represents the weight of the newly created block in the second attribute grid; Case C: such as Figure 4 As shown, the newly created block appears in three attribute grids. Assuming the three attribute grids have the same attribute and all have an attribute value of "5", then the newly created block will also have an attribute of "5". Figure 4As shown, the attributes of the attribute grid are "4", "5", and "6" respectively. Therefore, it is necessary to consider using a weighted method to calculate the attribute values ​​of this newly created block. The weighted formula is as follows:

[0047] In the formula, Y n This represents the attribute value of the target newly created block; A This represents the attribute value of the first attribute grid occupied by the newly created block. N 1 represents the weight of the newly created block in the first attribute grid; B This represents the attribute value of the second attribute grid occupied by the newly created block. N 2 represents the weight of the newly created block in the second attribute grid; C This represents the attribute value of the third attribute grid occupied by the newly created block. N 3 represents the weight of the newly created block in the third attribute grid; Case D: such as Figure 5 As shown, the newly created block appears in four attribute grids. Assuming these four attribute grids have the same attribute and all four attribute grids have an attribute value of "5", then the newly created block will also have an attribute of "5". If... Figure 5 As shown, the attributes of the attribute grid are "3", "4", "5" and "6". Therefore, a weighted approach needs to be used to calculate the attribute values ​​of this newly created block. The weighted formula is as follows:

[0048] In the formula, Y n This represents the attribute value of the target newly created block; A This represents the attribute value of the first attribute grid occupied by the newly created block. N 1 represents the weight of the newly created block in the first attribute grid; B This represents the attribute value of the second attribute grid occupied by the newly created block. N 2 represents the weight of the newly created block in the second attribute grid; C This represents the attribute value of the third attribute grid occupied by the newly created block. N 3 represents the weight of the newly created block in the third attribute grid; D This represents the attribute value of the fourth attribute grid occupied by the newly created block. N 4 represents the weight of the newly created block in the fourth attribute grid.

[0049] After determining the attribute values ​​of each newly created block and completing the attribute filling, the overall newly created block model will be smoother than the previous attribute data model, and the transition between two different attributes will be more gradual.

[0050] After determining the attribute values ​​for each grid cell, these values ​​are assigned to the volume center of the spatial cube. The cube is then divided using these data values, given a range of variation. Since the x, y, and z coordinates of each grid cell's centroid are assumed to be three data values, even a blank cubic grid can be divided into regions based on the coordinate value range. Similarly, a certain data index is first interpolated across the entire cubic grid to fill each cell with that index value, and then regions are divided based on the distribution range.

[0051] (7) Model visualization; Different attributes are displayed using different layers, and different colors are used to render attributes within the same layer to represent differences in attribute size, level, etc., and the display can be zoomed in / out, moved, sliced, and roamed.

[0052] This embodiment provides a method for 3D geological modeling and updating of tunnels based on attribute information. This includes determining the tunnel modeling range; converting various multi-source exploration data into borehole data for geometric structure modeling; performing attribute modeling under the constraints of the structural model; assigning and updating attribute values ​​to the model mesh; and visualizing the model. It fully utilizes exploration data to solve the problem of 3D geological models being "visually appealing but not practical." Through 3D geological modeling and visualization of the tunnel excavation model, it quickly and accurately obtains geological attribute information around the tunnel and integrates this information into a single 3D geological model. This model can provide information on geological structure, rock type, stratum thickness, water content, wave velocity, uniaxial compressive strength, and shear strength. c / φ Detailed information on parameters such as ) will enable more intuitive model data to be provided for the design, construction, and exploration layout of tunnel projects, and to provide scientific and reliable engineering geological basis and parameters for tunnel construction.

[0053] The tunnel 3D geological modeling and updating method provided in this embodiment only requires data regularization to establish an initial structural model during the modeling process. Then, the structural model is assigned values ​​and incremental data is added to update the model. This enables real-time control of the engineering geological conditions of the surrounding rock of the tunnel, especially the parameters related to support design. It is applicable to the construction of 3D geological models for the entire tunnel construction and excavation process under various complex geological conditions.

[0054] Example 2 This embodiment provides a tunnel 3D geological modeling and updating system oriented towards attribute information, which specifically includes: The model building module is configured to: acquire profile data and surface elevation data, build a tunnel structure model, and then build an attribute model constrained by the structure model. The model partitioning module is configured to: divide the attribute model into several grids, assign attribute values, and obtain the attribute grid; The model update module is configured to: in response to incremental data, construct an attribute-free grid and fill it with attributes according to the proportion of attribute grids to attribute-free grids, following the principle of weighted filling.

[0055] It should be noted that each module in this embodiment corresponds one-to-one with each step in Embodiment 1, and their specific implementation processes are the same, so they will not be repeated here.

[0056] Example 3 This embodiment provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the steps in the attribute-information-oriented tunnel 3D geological modeling and updating method described in Embodiment 1 above.

[0057] Example 4 This embodiment provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, it implements the steps in the tunnel three-dimensional geological modeling and updating method oriented towards attribute information as described in Embodiment 1 above.

[0058] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of hardware embodiments, software embodiments, or embodiments combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage and optical storage) containing computer-usable program code.

[0059] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0060] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0061] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0062] Those skilled in the art will understand that all or part of the processes in the above embodiments can be implemented by a computer program instructing related hardware. The program can be stored in a computer-readable storage medium, and when executed, it can include the processes of the embodiments of the above methods. The storage medium can be a magnetic disk, optical disk, read-only memory (ROM), or random access memory (RAM), etc.

[0063] The above description is merely a preferred embodiment of the present invention and is not intended to limit the invention. Various modifications and variations can be made to the present invention by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

Claims

1. A method for three-dimensional geological modeling and updating of tunnels based on attribute information, characterized in that, include: After acquiring profile data and surface elevation data, and establishing a tunnel structure model, an attribute model constrained by the structure model is established. The attribute model is divided into several grids, and attribute values ​​are assigned to obtain the attribute grid; In response to incremental data, an attribute-free grid is constructed, and attribute filling is performed according to the proportion of attribute grids to attribute-free grids, following the principle of weighted filling. include When a newly created block is completely enclosed by the attribute grid, it is determined that the newly created block fully displays the attributes in the attribute grid. When a newly created block appears in two attribute grids, if the attributes of the two attribute grids are different, then the weighted formula is applied. Calculate the attribute values ​​of the newly created block. Y n This represents the attribute value of the newly created target block; A This represents the attribute value of the first attribute grid occupied by the newly created block. N 1 represents the weight of the newly created block in the first attribute grid; B This represents the attribute value of the second attribute grid occupied by the newly created block. N 2 represents the weight of the newly created block in the second attribute grid; When a newly created block appears in three attribute grids, if the attributes of the three attribute grids are different, then the weighted formula is applied. Calculate the attribute values ​​of the newly created block, where C This represents the attribute value of the third attribute grid occupied by the newly created block. N 3 represents the weight of the newly created block in the third attribute grid; When a newly created block appears in one of the four attribute grids, it is calculated using the weighted formula. Calculate the attribute values ​​of the newly created block, where D This represents the attribute value of the fourth attribute grid occupied by the newly created block. N 4 represents the weight of the newly created block in the fourth attribute grid; The attribute model is divided into several regular hexahedral meshes, and the attribute value of each mesh is assigned to the body center of the regular hexahedron.

2. The method for three-dimensional geological modeling and updating of tunnels based on attribute information as described in claim 1, characterized in that, After regularizing the profile data and surface elevation data, a tunnel structure model was established.

3. The method for three-dimensional geological modeling and updating of tunnels based on attribute information as described in claim 1, characterized in that, The profile data includes borehole data and profile data.

4. The method for three-dimensional geological modeling and updating of tunnels based on attribute information as described in claim 3, characterized in that, Based on borehole data, the profile data is converted into borehole data and virtual borehole data to establish a tunnel structure model.

5. The method for three-dimensional geological modeling and updating of tunnels based on attribute information as described in claim 1, characterized in that, The attribute assignment uses three-dimensional spatial attribute interpolation.

6. The method for three-dimensional geological modeling and updating of tunnels based on attribute information as described in claim 5, characterized in that, The three-dimensional spatial attribute interpolation adopts radial basis function, Kriging and its variants, inverse distance weight, spline interpolation or Thiessen polygon interpolation method.

7. A tunnel 3D geological modeling and updating system oriented towards attribute information, the system being used to execute the tunnel 3D geological modeling and updating method oriented towards attribute information as described in any one of claims 1-6, characterized in that, include: The model building module is configured to: acquire profile data and surface elevation data, build a tunnel structure model, and then build an attribute model constrained by the structure model. The model partitioning module is configured to: divide the attribute model into several grids, assign attribute values, and obtain the attribute grid; The model update module is configured to: in response to incremental data, construct an attribute-free grid and fill it with attributes according to the proportion of attribute grids to attribute-free grids, following the principle of weighted filling.

8. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the program is executed by the processor, it implements the steps in the method for tunnel 3D geological modeling and updating based on attribute information as described in any one of claims 1-6.

9. A computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the program, it implements the steps in the tunnel three-dimensional geological modeling and updating method oriented towards attribute information as described in any one of claims 1-6.

Images