Underground engineering three-dimensional geological model dynamic visualization method, system, device and storage medium

Abstract

The present application relates to the technical field of geological model, especially to a three-dimensional geological model dynamic visualization method, system, equipment and storage medium for underground engineering. The technical scheme comprises the following steps: data acquisition and integration: obtaining geological data through various exploration methods, and performing standardized and fusion processing to obtain unified three-dimensional coordinate geological data; three-dimensional geological model initialization construction: constructing a stratigraphic framework based on the integrated data, and accurately embedding geological structures therein to form a complete three-dimensional geological model framework; model dynamic visualization rendering: setting specific lighting and material properties for the constructed three-dimensional geological model. The present application constructs a high-precision model through multi-source data fusion, optimizes the visualization effect and interactivity, realizes real-time monitoring and continuous optimization of the model, provides accurate geological basis for underground engineering, enhances decision accuracy, ensures safety, improves construction quality and economic benefits.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of geological model, and particularly relates to a three-dimensional geological model dynamic visualization method, system, device and storage medium for underground engineering. BACKGROUND

[0002] With the continuous development of underground engineering construction, such as subway tunnel excavation, mining, underground parking lot and large underground building construction, accurate understanding and visualization of underground geological conditions become increasingly critical. Traditional geological exploration and analysis methods gradually reveal many limitations when facing complex underground engineering needs.

[0003] In the early stage, geological exploration mainly relies on single drilling method. Although drilling can directly obtain underground core samples to determine the lithology, mineral composition and other information of the stratum, drilling information is only limited to the drilling position, and the regional geological conditions between the drill holes can only be speculated, and it is difficult to build a complete and continuous three-dimensional geological model. This leads to a large number of blind areas of geological information in the process of underground engineering planning and construction, which greatly increases the engineering risk, for example, when tunneling, unexpected faults, soft strata and other unfavorable geological bodies may be encountered, causing collapse, water inrush and other engineering accidents.

[0004] In order to make up for the limitations of drilling, geophysical exploration methods such as seismic wave exploration have been gradually applied. Seismic wave exploration transmits and receives seismic waves to infer stratum interfaces and geological structures according to the propagation characteristics of seismic waves in different strata. However, the interpretation of seismic wave exploration data has certain ambiguity and multiple solutions, and the geological model constructed by relying on seismic wave exploration data alone has limited accuracy, and it is difficult to effectively integrate with the fine lithology data obtained by drilling, and the advantages of the two kinds of data cannot be fully utilized.

[0005] At the same time, traditional geological model visualization is mostly static display, that is, even if a three-dimensional geological model is constructed, only fixed geological information at a specific stage is presented, and it cannot be dynamically updated as the underground engineering construction advances. In actual construction process, the stress state, hydrological conditions and other conditions of underground geological bodies will change, which may lead to stratum deformation, geological structure activation and other conditions, but the existing static visualization model cannot timely reflect these dynamic changes, making it difficult for engineering personnel to adjust the construction scheme according to the real-time geological conditions, which is not conducive to the safe and efficient advancement of underground engineering.

[0006] In addition, the existing visualization technology also has deficiencies in interactivity. When engineers view the geological model, they can only perform simple zooming and panning operations, and it is difficult to in-depth query detailed information of a specific geological body, which cannot meet the needs of full-range, multi-angle and in-depth analysis of geological information in the decision-making process of complex underground engineering. For example, when dealing with large underground complex engineering, the development of underground space with different purposes involves fine analysis of geological conditions in different regions, and the traditional visualization method cannot provide effective support.

[0007] In summary, it has extremely important practical significance to develop an underground engineering three-dimensional geological model dynamic visualization method capable of integrating multi-source geological data, constructing a high-precision and dynamically updated three-dimensional geological model, and having strong interactive functions. It will provide strong technical support for the design, construction, safety monitoring and management of underground engineering, effectively improve the quality and efficiency of underground engineering construction, and reduce engineering risks and costs. Therefore, the present application proposes an underground engineering three-dimensional geological model dynamic visualization method, system, device and storage medium. SUMMARY

[0008] The purpose of the present application is to solve the problem that there is no three-dimensional geological model capable of integrating multi-source geological data and constructing a high-precision and dynamically updated three-dimensional geological model in the background art, and proposes an underground engineering three-dimensional geological model dynamic visualization method, system, device and storage medium.

[0009] In a first aspect, the present application provides an underground engineering three-dimensional geological model dynamic visualization method, comprising the following steps:

[0010] Data acquisition and integration: obtain geological data through various exploration methods, and perform standardization and fusion processing to obtain unified three-dimensional coordinate geological data;

[0011] Three-dimensional geological model initialization construction: based on the integrated data, construct a stratigraphic framework, and accurately embed geological structures therein to form a complete three-dimensional geological model framework;

[0012] Model dynamic visualization rendering: set specific lighting and material properties for the constructed three-dimensional geological model, and develop a visualization program to realize dynamic display and interactive operation of the model;

[0013] Model dynamic updating and optimization: during the construction of underground engineering, monitor the changes of geological bodies with the aid of a sensor network, update and optimize the model according to the monitoring data, and ensure that the model is highly consistent with the actual geological conditions.

[0014] Optionally, in the data acquisition and integration, the multi-source data acquisition is specifically:

[0015] Using high-precision drilling equipment, drilling according to the predetermined exploration grid, drilling depth interval is 0.5 meters, to obtain the lithology, mineral composition, porosity and depth information of the core sample;

[0016] Using seismic wave tomography exploration, transmitting different frequency seismic waves with 2m x 2m measuring point spacing, constructing the wave velocity model of the underground geological body according to the reflection wave characteristics, and inferring the stratigraphic interface and geological structure;

[0017] Carrying out ground geological survey, recording the surface stratum lithology and geological structure characteristics, and correlating and calibrating with drilling and seismic wave data.

[0018] Optionally, in the data collection and integration, the data standardization and fusion specifically includes the following steps:

[0019] Converting the drilling data, seismic wave data and ground geological survey data into three-dimensional coordinate data format with the geodetic coordinate system as the reference, wherein the x and y coordinates are accurate to 0.01 meters, and the z coordinate is accurate to 0.1 meters;

[0020] For stratigraphic interface data, weighted average method is used for fusion, assuming that the stratigraphic interface depth determined by drilling is d d , the stratigraphic interface depth inferred by seismic wave is d s , their weights are w d =0.6 and w s =0.4 respectively, and the fused stratigraphic interface depth d f is calculated by the formula d f =w d ×d d +w s ×d s .

[0021] Optionally, in the three-dimensional geological model initialization construction, the method for constructing the stratigraphic framework includes:

[0022] Based on the fused stratigraphic interface data, a Delaunay triangulation algorithm is used to construct a stratigraphic triangular facet model, adjacent stratigraphic interfaces are filled with triangular prismatic units, and the mismatch error of the common surface of adjacent triangular prismatic units is checked, and if the error is greater than 0.05 square meters, local adjustment and optimization are performed;

[0023] The volume and centroid of the stratigraphic model are calculated, the volume V t of the stratigraphic model is obtained by summing the volumes of the triangular prisms, the volume of the triangular prism is calculated by vector cross product method, assuming that the vertices of the triangular prism are A(x a ,y a ,z a ), B(x b ,y b ,z b ), and C(x c ,y cz c ), high h, then its volume:

[0024]

[0025] where, [A-C, B-C] represents the cross product of vector A-C and B-C, formation centroid coordinates (x g ,y g ,z g ) are calculated according to the volume weighted average of each triangular prism, the formula is:

[0026]

[0027] n is the number of triangular prisms in the formation.

[0028] Optionally, in the three-dimensional geological model initialization construction step, the geological structure embedding manner is:

[0029] For fault modeling, according to the fault strike angle α, the tendency angle β, the dip angle γ and the fault throw Δh, the fault surface is constructed by the parametric equation at the corresponding formation position, assuming that the point P (x, y, z) on the fault surface, taking the known point P0 (x0, y0, z0) as the base, then x = x0 + t × cos α × cos γ, y = y0 + t × sin α × cos γ, z = z0 + t × sin γ, t is the parameter, and the formation on both sides of the fault surface is moved according to the fault throw;

[0030] For fold structure, according to the fold axis equation L: ax + by + cz + d = 0 and the fold amplitude A f , the stratum model is deformed by using the surface deformation algorithm, assuming that the distance from the point Q (x q ,y q ,z q ) on the stratum model to the fold axis is d q , the coordinates Q' (x' q ,y' q ,z' q ) after deformation are calculated by vertical axis displacement, and the displacement amount is T is the fold wavelength, which is determined according to geological survey.

[0031] Optionally, in the model dynamic visual rendering, the features of light and material setting are:

[0032] The ambient light, parallel light and point light source are set for the three-dimensional geological model, the ambient light intensity I a is 0.3, the parallel light direction vector is determined according to the observation direction of the underground engineering, the parallel light intensity I p is 0.6, and the point light source is 5 meters above the center of the model, and the intensity I sis 0.8, point light attenuation coefficient k a = 0.05, k b = 0.002, k c = 0.0001;

[0033] Different lithologic strata and geological structures are assigned with material properties, hard rock strata have a diffuse reflection color of (0.5, 0.5, 0.5), a specular reflection color of (0.2, 0.2, 0.2), and a glossiness of 30; soft rock strata have a diffuse reflection color of (0.7, 0.7, 0.3), a specular reflection color of (0.1, 0.1, 0.1), and a glossiness of 10; a fault has a diffuse reflection color of (1, 0, 0), a specular reflection color of (0.5, 0, 0), a glossiness of 5, and a fault transparency of 0.3.

[0034] Optionally, the dynamic display and interaction are manifested as:

[0035] A visualization program is developed based on a graphical programming interface, and a frame rate of 60 frames per second is set to realize real-time rendering and dynamic display of the model;

[0036] User interaction functions are set, the left mouse button is used to rotate the model around the y-axis with a rotation angle step of 0.5 degrees, the right mouse button is used to translate the model with a step of 0.05 meters, the scroll wheel is used to scale the model with a scale factor of 1.1 times per scale, and shortcut keys are set, i.e., the “F1” key is used to control the display or hiding of the stratum model, the “F2” key is used to control the display or hiding of the geological structure model, and the “F3” key is used to switch the lighting mode.

[0037] In a second aspect, the application provides a dynamic visualization system for a three-dimensional geological model of an underground engineering, which is used to realize the dynamic visualization method for the three-dimensional geological model of the underground engineering according to the first aspect, and includes:

[0038] A data acquisition and integration module: geological data are obtained through various exploration methods, and standardized and fused to obtain unified three-dimensional coordinate geological data;

[0039] A three-dimensional geological model initialization construction module: a stratum framework is constructed based on the integrated data, and geological structures are accurately embedded therein to form a complete three-dimensional geological model framework;

[0040] A model dynamic visualization rendering module: specific lighting and material properties are set for the constructed three-dimensional geological model, and a visualization program is developed to realize dynamic display and interactive operation of the model;

[0041] A model dynamic updating and optimization module: during the construction of the underground engineering, geological body changes are monitored by means of a sensor network, and the model is updated and optimized according to the monitoring data to ensure that the model is highly consistent with the actual geological conditions.

[0042] In a third aspect, the application provides a device for dynamic visualization of a three-dimensional geological model of underground works, comprising:

[0043] a processor;

[0044] a memory having executable instructions of the processor stored therein;

[0045] wherein the processor is configured to execute the steps of the method for dynamic visualization of a three-dimensional geological model of underground works according to the first aspect via execution of the executable instructions.

[0046] In a fourth aspect, the application provides a computer-readable storage medium for storing a program, which, when executed by a processor, implements the steps of the method for dynamic visualization of a three-dimensional geological model of underground works according to the first aspect.

[0047] Compared with the prior art, the application has at least one of the following beneficial technical effects:

[0048] The integration of multi-source geological data and accurate fusion can construct a high-precision three-dimensional geological model, accurately depict the details of strata and geological structures, reduce the uncertainty of geological information, and provide reliable basis for engineering planning.

[0049] Through unique lighting and material settings, different geological bodies can be visually distinguished in the model, enhancing the understanding and cognition of the underground geological environment and improving the accuracy of decision-making.

[0050] Through the high-frame-rate visualization program combined with convenient mouse and keyboard interaction and shortcut functions, the model can be observed from multiple angles and with different levels of detail, geological information can be efficiently obtained, and work efficiency and analysis comprehensiveness can be improved.

[0051] In construction, the sensor network monitors the changes of geological bodies in real time, triggers model updating when the threshold is exceeded, provides safety warnings, and ensures the safety of construction personnel and the smooth progress of the project.

[0052] Based on the monitoring data, the model parameters are inverted and optimized to keep the model highly consistent with the actual geological conditions, continuously provide reliable support for construction decisions, and improve construction quality and economic benefits.

[0053] The application constructs a high-precision model through multi-source data fusion, optimizes the visualization effect and interactivity, realizes real-time monitoring and continuous optimization of the model, provides accurate geological basis for underground works, enhances decision-making accuracy, ensures safety, and improves construction quality and economic benefits. BRIEF DESCRIPTION OF DRAWINGS

[0054] Figure 1 Flowchart of the method for dynamic visualization of a three-dimensional geological model of underground works; DETAILED DESCRIPTION

[0055] Following, the embodiments of the present application are illustrated by specific, concrete examples, and those skilled in the art can easily understand other advantages and effects of the present application from the disclosure of the present application. The present application can also be implemented or applied by other different specific embodiments, and various modifications or changes can be made to the details in the present application without departing from the spirit of the present application. It should be noted that the embodiments in the present application and the features in the embodiments can be combined with each other without conflict.

[0056] Embodiment 1

[0057] As shown in the figure, the underground engineering three-dimensional geological model dynamic visualization method provided by the present application comprises the following steps: Figure 1

[0058] I. Data acquisition and integration

[0059] 1. Multi-source data acquisition

[0060] High-precision drilling equipment is used to drill according to the predetermined exploration grid, to obtain core samples at different depths, analyze the lithology, mineral composition, and porosity parameters, and record the corresponding depth information. The drilling depth interval is set to 0.5 meters to ensure the fineness of the data.

[0061] Seismic wave tomography exploration is implemented, different frequency seismic waves are emitted, reflected wave signals are received, and the wave velocity model of the underground geological body is constructed by analyzing the time delay, amplitude change and other characteristics of the reflected waves, so as to infer the position and shape of the stratigraphic interface and geological structure. The point spacing of the seismic wave exploration is 2 meters x 2 meters.

[0062] Ground geological survey is carried out, and the lithology and geological structure characteristics of the exposed strata on the ground are recorded in detail, including the location, strike, and dip angle of the fault outcrop, and are associated and calibrated with the drilling and seismic wave exploration data.

[0063] 2. Data standardization and fusion

[0064] Drilling data, seismic wave data, and ground geological survey data are uniformly converted into three-dimensional coordinate data format with the geodetic coordinate system as the reference. Among them, the x and y coordinates are accurate to 0.01 meters, and the z coordinate (depth) is accurate to 0.1 meters.

[0065] For the same geological feature from different data sources, a data fusion algorithm is established. For stratigraphic interface data, a weighted average method is used for fusion, and the stratigraphic interface depth determined by drilling is d d , and the stratigraphic interface depth inferred by seismic wave is d s , and their weights are w d = 0.6 and w s ​= 0.4, depth of the merged stratigraphic interface d f The calculation formula is d f = w d × d d + w s × d s By integrating multi-source data such as high-precision drilling, seismic wave tomography exploration, and surface geological survey, the strengths of each data source are fully utilized. Drilling data provides accurate lithology and depth information, seismic wave data helps to infer large-scale stratigraphic interfaces and geological structures, and surface geological survey supplements key information from surface outcrops. Using data fusion algorithms such as weighted averaging, the fusion processing of stratigraphic interface data makes the constructed three-dimensional geological model more accurate and complete in spatial distribution and geological feature description, reducing errors and uncertainties that may be caused by single data sources, and providing reliable geological basis for the preliminary planning of underground engineering.

[0066] Initialization and construction of three-dimensional geological model

[0067] 1. Stratigraphic framework construction

[0068] Based on the fused stratigraphic interface data, a triangular facet model of the strata is constructed using the Delaunay triangulation algorithm. For adjacent stratigraphic interfaces, a three-prism unit is constructed to fill the space between the strata. During the construction process, the connectivity of adjacent three-prism bodies is checked. If the mismatch error of the common face of adjacent three-prism bodies exceeds 0.05 square meters, local adjustment and optimization are performed.

[0069] 2. Calculate the volume and centroid of the stratigraphic model. The calculation of the volume V t of the strata is obtained by summing the volumes of all three-prism bodies. The volume of a three-prism body is calculated using the vector cross product method. Let the coordinates of the three vertices of the three-prism body be A(x a ,y a ,z a ), B(x b ,y b ,z b ), and C(x c ,y c ,z c ), and the height be h, then the volume is:

[0070]

[0071] where [A-C, B-C] represents the cross product of vectors A-C and B-C, and the centroid coordinates (x g ,y g ,z g ) of the strata are calculated by weighted average of the volumes of each three-prism body, with the formula:

[0072]

[0073] n is the number of prisms in the formation.

[0074] 3. Geological structure embedding

[0075] For fault modeling, a fault plane is constructed at the corresponding formation position according to the geometric parameters and the fault throw Ah of the fault, the geometric parameters including the strike angle a, the trend angle b, and the dip angle g. The fault plane is represented by a parametric equation, and a known point P0(x0, y0, z0) on the fault plane is taken as the reference. Then x = x0 + t x cos a x cos g, y = y0 + t x sin a x cos g, and z = z0 + t x sin g, where t is the parameter. The strata on both sides of the fault plane are moved according to the fault throw to ensure the geometric consistency of the strata and the fault plane.

[0076] For fold structure, the strata model is deformed using a surface deformation algorithm according to the fold axis equation: a straight line equation L: ax + by + cz + d = 0 and the fold amplitude A f . Let the distance from a point Q(x q , y q , z q ) on the strata model to the fold axis be d q , and the deformed coordinates Q'(x' q , y' q , z' q ) are calculated by displacing the Q point along the direction perpendicular to the axis. The displacement amount T is the fold wavelength, which is determined according to geological surveys. In the initialization and construction process of the three-dimensional geological model, advanced algorithms are used to construct the strata framework and embed geological structures based on the fused data. For example, the Delaunay triangulation algorithm is used to construct the strata triangular facet model and fill the prismatic elements, which can well adapt to complex strata morphology and accurately depict the geometric characteristics of the strata. The modeling method for faults and folds and other geological structures is constructed according to their geometric parameters and deformation characteristics, so that the model can accurately reflect the true situation of the geological structure, whether it is the continuity of the strata, the spatial form of the geological structure, or the mutual relationship between them, which can help engineers to predict the possible geological problems in underground engineering construction, such as the treatment of fault fracture zones and the stability analysis of fold regions.

[0077] III. Dynamic visualization rendering of the model

[0078] 1. Light and material settings

[0079] Set the environment light, parallel light, and point light sources for the three-dimensional geological model. The environment light intensity I a is set to 0.3, the parallel light direction vector According to the main observation direction of the underground project, represents the downward irradiation from the top, the parallel light intensity I p = 0.6, the point light source is located 5 meters above the center of the model, the intensity I s = 0.8, the attenuation coefficient k a = 0.05, k b = 0.002, k c = 0.0001, used to calculate the attenuation of the point light source illumination intensity with distance.

[0080] Assign material properties to different lithological formations and geological structures. The diffuse reflection color of hard rock formation is (0.5, 0.5, 0.5), the specular reflection color is (0.2, 0.2, 0.2), and the glossiness is 30; the diffuse reflection color of soft rock formation is (0.7, 0.7, 0.3), the specular reflection color is (0.1, 0.1, 0.1), and the glossiness is 10; the diffuse reflection color of fault is (1, 0, 0), the specular reflection color is (0.5, 0, 0), the glossiness is 5, and the transparency of fault is set to 0.3. Carefully set lighting and material properties significantly improve the visualization effect of the model. Different lithological formations and geological structures have unique material characteristics such as diffuse reflection color, specular reflection color, glossiness, and transparency. For example, hard rock formation and soft rock formation can be clearly distinguished visually, and fault is highlighted with bright red color and specific transparency, which enables engineers to intuitively identify different geological bodies when viewing the model, as if they were in the underground geological environment, enhancing their understanding and cognition of the geological situation, and helping to improve the accuracy and efficiency of decision-making.

[0081] 2. Dynamic display and interaction

[0082] Develop a visualization program based on a graphical programming interface (including OpenGL or DirectX) to realize real-time rendering and dynamic display of the model. Set the frame rate to 60 frames per second to ensure smooth display of the model.

[0083] The user interaction function is realized, and the rotation, translation and scaling of the model are controlled by mouse and keyboard operation. The left mouse button is pressed and dragged to realize the rotation of the model around the y-axis, and the rotation angle step is 0.5 degrees; the right mouse button is pressed and dragged to realize the translation of the model, and the translation step is 0.05 meters; the mouse wheel is rolled to realize the scaling of the model, and the scaling factor is 1.1 times / wheel scale. At the same time, shortcut keys are set, such as "F1" key to display or hide the stratigraphic model, "F2" key to display or hide the geological structure model, and "F3" key to switch different lighting modes. The developed visualization program has high frame rate (60 frames / second) and rich interactive functions. Through mouse and keyboard operation, flexible rotation, translation and scaling of the model are realized, and the operation step is carefully designed, such as rotation angle step of 0.5 degrees, translation step of 0.05 meters, and scaling factor of 1.1 times / wheel scale. At the same time, convenient shortcut keys are set, such as "F1", "F2" and "F3" to control the display and hiding of the stratigraphic model, the geological structure model and the switching of the lighting mode respectively. These interactive functions facilitate engineers to observe and analyze the model from any angle and any detail level, quickly focus on the area or geological body of interest, efficiently obtain the required geological information, and greatly improve the work efficiency and the comprehensiveness of analysis.

[0084] Four, model dynamic updating and optimization

[0085] 1. Data monitoring and updating trigger

[0086] In the process of underground engineering construction, a sensor network is arranged, including displacement sensors, stress sensors, underground water level sensors, etc., to monitor the deformation, stress change and hydrogeological condition change of the geological body in real time. The measurement accuracy of the displacement sensor is 0.001 meters, and the measurement accuracy of the stress sensor is 0.1 MPa.

[0087] When the change amplitude of the monitoring data exceeds the preset threshold, the model updating is triggered. For example, when the stratigraphic displacement change exceeds 0.1 meters or the stress change exceeds 1 MPa, the model updating program is started. In the process of underground engineering construction, the sensor network arranged can monitor the deformation, stress change and hydrogeological condition change of the geological body in real time, with accuracy of 0.001 meters for displacement sensor and 0.1 MPa for stress sensor. Once the change amplitude of the monitoring data exceeds the preset threshold, such as the stratigraphic displacement change exceeding 0.1 meters or the stress change exceeding 1 MPa, the model updating program is triggered in time. This provides a real-time safety monitoring and early warning mechanism for underground engineering construction, and engineers can know in advance the occurrence trend of potential geological risks such as stratigraphic collapse and underground water inflow, so as to take corresponding preventive measures in time, ensuring the life safety of construction personnel and the smooth progress of the project.

[0088] 2. Model optimization and calibration

[0089] According to the new monitoring data, the parameters of the geological model are optimized by using a back analysis algorithm. For example, the elastic modulus and Poisson's ratio of the stratum are inverted by using the monitored stratum displacement data. Let the actual displacement of the displacement monitoring point be The displacement calculated by the model is By minimizing the objective function:

[0090]

[0091] k is the number of monitoring points, and the gradient descent method is used to adjust the elastic modulus E and Poisson's ratio v of the stratum. During the optimization process, the adjustment step of the elastic modulus is 0.1 GPa, and the adjustment step of the Poisson's ratio is 0.01.

[0092] The optimized model is compared and calibrated with the new monitoring data, and the error index of the model and the actual data is calculated:

[0093]

[0094] wherein d j is the jth data value calculated by the model, d mj is the jth monitoring data value, and l is the number of comparison data. If the error index is greater than 0.1, the model optimization is continued until the error index meets the requirements.

[0095] Based on the new monitoring data, the parameters of the geological model are optimized by using a back analysis algorithm, for example, the inversion and adjustment of the elastic modulus and Poisson's ratio of the stratum. By minimizing the objective function and using optimization algorithms such as the gradient descent method, combined with reasonable parameter adjustment steps (the elastic modulus adjustment step is 0.1 GPa, and the Poisson's ratio adjustment step is 0.01), the model can continuously adapt to the actual changes of the geological body and always remain highly consistent with the underground geological conditions. After the optimized model is compared and calibrated with the new monitoring data, it is ensured that the error index meets the requirements, thereby providing continuous and reliable basis for subsequent construction decisions, such as reasonably adjusting the construction progress, optimizing the support scheme, effectively improving the construction quality and economic benefits of underground engineering.

[0096] Example 2

[0097] The embodiment provides a kind of underground engineering three-dimensional geological model dynamic visualization system, for realizing the underground engineering three-dimensional geological model dynamic visualization method described in example 1, it includes:

[0098] Data acquisition and integration module: obtain geological data by various exploration means, and carry out standardization and fusion processing, obtain unified three-dimensional coordinate geological data;

[0099] A three-dimensional geological model initialization construction module: based on the integrated data, a stratigraphic framework is constructed, and geological structures are accurately embedded therein to form a complete three-dimensional geological model framework;

[0100] A model dynamic visualization rendering module: specific lighting and material properties are set for the constructed three-dimensional geological model, and a visualization program is developed to realize dynamic display and interactive operation of the model;

[0101] A model dynamic updating and optimization module: during the construction of underground engineering, the sensor network is used to monitor the change of geological body, and the model is updated and optimized according to the monitoring data, so as to ensure that the model is highly consistent with the actual geological conditions.

[0102] Embodiment 3

[0103] The embodiment provides a kind of underground engineering three-dimensional geological model dynamic visualization equipment, comprising:

[0104] Processor;

[0105] Memory, wherein the executable instruction of the processor is stored;

[0106] Wherein, the processor is configured to execute the executable instruction to execute the steps of the underground engineering three-dimensional geological model dynamic visualization method described in embodiment 1.

[0107] Embodiment 4

[0108] The embodiment provides a kind of computer readable storage medium, for storing program, the program is executed by processor when realizing the steps of the underground engineering three-dimensional geological model dynamic visualization method described in embodiment 1.

[0109] The above specific embodiments are only several optional embodiments of the present application, based on the technical scheme of the present application and the related inspiration of the above embodiment, the person skilled in the art can make various alternative improvements and combinations to the above specific embodiments.

Claims

1. A method for dynamic visualization of a three-dimensional geological model of an underground structure, characterized in that, The method comprises the following steps: Data acquisition and integration: Obtain geological data through various exploration methods, and perform standardization and fusion processing to obtain unified three-dimensional coordinate geological data; Three-dimensional geological model initialization construction: Based on the integrated data, construct a stratigraphic framework, and accurately embed geological structures therein to form a complete three-dimensional geological model framework; Model dynamic visualization rendering: Set specific lighting and material properties for the constructed three-dimensional geological model, and develop a visualization program to realize dynamic display and interactive operation of the model; Model dynamic updating and optimization: During the construction of underground projects, monitor the changes of geological bodies with the aid of a sensor network, update and optimize the model according to the monitoring data, and ensure that the model is highly consistent with the actual geological conditions; The parameters of the geological model are optimized by using a back analysis algorithm, and the elastic modulus and Poisson's ratio of the stratum are inverted by using the monitored stratum displacement data. The actual displacement of the displacement monitoring point is , the displacement calculated by the model is , and the objective function is minimized by ； For monitoring the number of points, the gradient descent method is used to adjust the elastic modulus of the formation And Poisson's ratio In the optimization process, the adjustment step of the elastic modulus is 0.1 GPa, and the adjustment step of Poisson's ratio is 0.01; Compare and calibrate the optimized model and new monitoring data to calculate the error index of the model and the actual data: ； wherein, the first data value calculated for the model, the first data value calculated for the model, the first monitoring data value, the first monitoring data value, the number of comparison data, if the error index is greater than 0.1, the model optimization is continued until the error index meets the requirements.

2. The method of claim 1, wherein, In the data acquisition and integration, the multi-source data acquisition specifically comprises the following steps: High-precision drilling equipment is used to drill at a predetermined exploration grid, the drilling depth interval is 0.5 meters, and the lithology, mineral composition, porosity and depth information of the core sample are obtained; Seismic wave tomography exploration is used, different frequency seismic waves are emitted at a 2m x 2m measurement point interval, and a wave velocity model of the underground geological body is constructed according to the reflection wave characteristics to deduce the stratigraphic interface and geological structure; Ground geological survey is carried out to record the surface stratigraphic lithology and geological structure characteristics, and the data are associated and calibrated with the drilling and seismic wave data.

3. The method of claim 2, wherein, In the data acquisition and integration, the data standardization and fusion specifically comprises the following steps: The drilling data, seismic wave data and ground geological survey data are converted into three-dimensional coordinate data format with the geodetic coordinate system as the reference, wherein the x and y coordinates are accurate to 0.01 meters, and the z coordinate is accurate to 0.1 meters; For the stratum interface data, the weighted average method is used for fusion, the stratum interface depth determined by drilling is , the stratum interface depth inferred by the seismic wave is , the weights are and respectively, the fused stratum interface depth is , and the calculation formula is .

4. The method of claim 1, wherein, In the three-dimensional geological model initialization construction, the method for constructing the stratigraphic framework comprises the following steps: Based on the integrated stratigraphic interface data, a stratigraphic triangular facet model is constructed using the Delaunay triangulation algorithm, adjacent stratigraphic interfaces are filled with triangular prism units, and the mismatch error of the common surface of adjacent triangular prisms is checked; if the error is greater than 0.05 square meters, local adjustment and optimization are performed; Computing the stratigraphic model volume and centroid, stratigraphic volume Summing the volumes of the triangular prisms, the volume of the triangular prism is computed using the cross product of vectors, let the triangular prism have vertices , , , and height , then its volume is: ； wherein, denotes the cross product of the vectors and the formation centroid coordinates are calculated as a volume-weighted average of the individual prisms, according to the formula: ； ； ； is the number of prisms in the formation.

5. The method of claim 4, wherein, In the three-dimensional geological model initialization construction step, the geological structure is embedded in the following manner: For fault modeling, according to the strike angle , the trend angle , the dip angle and the fault throw , the fault plane is constructed by parametric equation at the corresponding stratum position, and the point on the fault plane is set. With the known point as the base, the parameter , , , is obtained, and the strata on both sides of the fault plane are dislocated according to the fault throw. For the fold structure, according to the fold axis equation and the fold amplitude , the stratum model is deformed by the surface deformation algorithm, and the point on the stratum model is set as the distance to the fold axis , the coordinate after deformation is calculated by the vertical axis displacement, and the displacement amount , is the fold wavelength, which is determined according to the geological survey.

6. The method of claim 1, wherein, In the model dynamic visualization rendering, the lighting and material setting features are as follows: Ambient light, parallel light and point light source are set for the three-dimensional geological model, the ambient light intensity is 0.3, the parallel light direction vector is determined according to the observation direction of the underground engineering, the parallel light intensity is 0.6, the point light source is 5 meters above the center of the model, the intensity of the point light source is 0.8, and the point light source attenuation coefficient , , ; For different lithological strata and geological structure to define material properties, hard rock strata diffuse reflection color is , mirror reflection color is , glossiness 30; soft rock strata diffuse reflection color is , mirror reflection color is , glossiness 10; fault diffuse reflection color , mirror reflection color is , glossiness 5, and fault transparency 0.

3.

7. The method of claim 6, wherein, The dynamic display and interaction are as follows: A visualization program is developed based on a graphical programming interface, and the frame rate is set to 60 frames per second to realize real-time rendering and dynamic display of the model; User interaction functions are set, the left mouse button is used to rotate the model around the y-axis with a rotation angle step of 0.5 degrees, the right mouse button is used to translate the model with a step of 0.05 meters, and the scroll wheel is used to scale the model with a scale factor of 1.1 times per scale; and shortcut keys are set, "F1" is used to control the display and hiding of the stratigraphic model, "F2" is used to control the display and hiding of the geological structure model, and "F3" is used to switch the lighting mode.

8. A system for dynamic visualization of a three-dimensional geological model of an underground structure, for implementing the method for dynamic visualization of a three-dimensional geological model of an underground structure according to claim 1, characterized in that, It comprises: A data acquisition and integration module: geological data is obtained through various exploration methods, and standardized and fused to obtain unified three-dimensional coordinate geological data; A three-dimensional geological model initialization construction module: based on the integrated data, a stratigraphic framework is constructed, and geological structures are accurately embedded therein to form a complete three-dimensional geological model framework; Model dynamic visualization rendering module: set specific light and material properties for the constructed three-dimensional geological model, and develop a visualization program to realize dynamic display and interactive operation of the model; Model dynamic updating and optimization module: during the construction of underground engineering, monitor the changes of geological bodies with the help of sensor network, update and optimize the model according to the monitoring data, and ensure that the model is highly consistent with the actual geological conditions.

9. An underground engineering three-dimensional geological model dynamic visualization device, characterized in that, Comprise: a processor; a memory, wherein executable instructions of the processor are stored; wherein the processor is configured to execute the steps of the underground engineering three-dimensional geological model dynamic visualization method of any one of claims 1-7 by executing the executable instructions.

10. A computer readable storage medium for storing a program, characterized in that, The program is executed by the processor to realize the steps of the underground engineering three-dimensional geological model dynamic visualization method of any one of claims 1-7.

Images