A method and system for constructing a numerical calculation model of a slope based on an attribute-containing three-dimensional geological model

Abstract

The application discloses a kind of based on attribute-containing three-dimensional geological model's slope numerical calculation model construction method and system, comprising: constructing attribute-containing three-dimensional geological model;Generation slope excavation contour surface and according to slope reinforcement scheme arrangement slope reinforcement piece, establish slope three-dimensional design model;From three-dimensional geological model, cut out slope three-dimensional design model, and carry out parameterization setting, obtain slope numerical calculation model;According to load condition and model attribute in three-dimensional geological model, intelligent setting boundary condition of slope numerical calculation model, finally construct and obtain the slope numerical calculation model containing boundary condition.The structure surface of slope numerical calculation model is designed.The application enhances the relevance and inheritance between models, and intelligently sets reasonable boundary conditions.The structural surface network constructed by the application has high accuracy, and solves the problem of the influence of the deviation of the structural surface model from reality on the slope stability evaluation in the numerical calculation of rock slope.

Description

Technical Field

[0001] This invention relates to the field of intelligent design technology for geotechnical engineering, specifically to a method and system for constructing a slope numerical calculation model based on a three-dimensional geological model with attributes. Background Technology

[0002] With the increasing application of 3D geological models, the demand for technologies such as slope excavation design based on 3D geological models and intelligent cloud computing in the design process is becoming increasingly urgent. In the slope excavation design process, intelligent numerical calculation and analysis of slopes using numerical methods is a crucial step. Currently, the construction process of slope numerical calculation models and the setting of boundary conditions for these models are quite cumbersome, thus limiting the large-scale application of slope engineering safety assessments based on slope numerical calculation models. Summary of the Invention

[0003] To address the shortcomings of existing technologies, this invention provides a method and system for constructing a slope numerical calculation model based on a three-dimensional geological model with attributes, which can effectively solve the above-mentioned problems.

[0004] The technical solution adopted in this invention is as follows:

[0005] This invention provides a method for constructing a slope numerical calculation model based on a three-dimensional geological model with attributes, including:

[0006] Step S1: Construct a three-dimensional geological model with attributes; the three-dimensional geological model with attributes includes a three-dimensional geological model and model attributes; wherein, the three-dimensional geological model includes a surface object at the top, a groundwater level object at the bottom, and several surface objects between the surface object and the groundwater level object; the model attributes include the geological classification attributes of each surface object and the physical and mechanical index attributes of each surface object.

[0007] Step S2: According to the engineering requirements, in the attributed three-dimensional geological model, generate the slope excavation outline surface and arrange the slope reinforcement components according to the slope reinforcement scheme. Assign the geological classification attributes of the excavation surface object to the slope excavation outline surface, thereby establishing a three-dimensional slope design model containing geological conditions.

[0008] Step S3: Using the closed intersection algorithm, the slope three-dimensional design model containing geological conditions is cut out from the attribute-containing three-dimensional geological model, and the cut slope three-dimensional design model containing geological conditions is parameterized to obtain the slope numerical calculation model.

[0009] Step S4: Determine the load conditions based on the numerical simulation working conditions. Based on the load conditions and the model attributes in the attributed three-dimensional geological model, intelligently set the boundary conditions of the slope numerical calculation model, and finally construct the slope numerical calculation model with boundary conditions.

[0010] Furthermore, the physical and mechanical properties of each surface object include the deformation modulus, Poisson's ratio, and strength parameters of the surface object's soil and rock mass.

[0011] Furthermore, step S2 specifically involves:

[0012] Step S21: In the attribute-containing three-dimensional geological model, a slope excavation start outline is set. Based on the slope excavation start outline, slope excavation outlines with different slope ratios and slope heights are set from near to far. The slope excavation start outline and each of the slope excavation outlines are combined to form the slope excavation outline surface.

[0013] Step S22: According to the slope reinforcement scheme, select a reinforcement zone in the slope excavation outline; in the reinforcement zone, arrange the slope reinforcement components according to the set spacing and row spacing of the slope reinforcement components.

[0014] Step S23: Assign geological classification attributes to the slope excavation outline surface; the slope excavation outline surface and the slope reinforcement components form the three-dimensional design model of the slope containing geological conditions.

[0015] Furthermore, step S3 specifically involves: generating a three-dimensional slope numerical calculation model using a three-dimensional closed intersection technique, including:

[0016] Step S31: Based on the slope excavation outline surface of the slope three-dimensional design model containing geological conditions and the design requirements, determine the spatial range of the slope three-dimensional design model to be cut out, so that it covers the slope excavation outline surface in the horizontal direction and extends from the ground surface object to the bottom plane that meets the design requirements in the vertical direction.

[0017] Step S32: Using the three-dimensional closed intersection technique, the three-dimensional slope design model designed in the attribute-containing three-dimensional geological model is subjected to three-dimensional closed intersection to determine the four sides closed plane and the bottom closed plane of the slope three-dimensional design model.

[0018] Step S33: Based on the closed planes on all four sides and the closed plane at the bottom, cut out the three-dimensional design model of the slope from the three-dimensional geological model containing attributes;

[0019] Step S34, Set boundary surface properties:

[0020] In the three-dimensional design model of the slope obtained by trimming, the top surface object inherits the original geological classification attributes and is given surface boundary surface attributes; the top excavation surface object inherits the original geological classification attributes and is given excavation surface boundary surface attributes; the surrounding closed planes are given lateral boundary surface attributes; and the bottom closed plane is given bottom boundary surface attributes.

[0021] By setting each boundary surface, the geometric entity of the three-dimensional design model of the slope is formed;

[0022] Step S35, Material Zoning:

[0023] Suppose that the attributed three-dimensional geological model has N layers of ground objects in the spatial range corresponding to the three-dimensional slope design model, which are from top to bottom: the first layer of ground object, the second layer of ground object, ..., the Nth layer of ground object;

[0024] The first ground surface object and the slope three-dimensional design model intersect to form a first boundary interface; the first boundary interface and the ground surface boundary surface enclose a first material partition; the first material partition inherits the physical and mechanical properties of the first ground surface object;

[0025] The second ground surface object and the slope three-dimensional design model intersect to form a second boundary interface; the second boundary interface and the first boundary interface immediately above it form a second material partition; the second material partition inherits the physical and mechanical properties of the second ground surface object;

[0026] And so on

[0027] The Nth ground surface object and the slope 3D design model intersect to form the Nth boundary interface; the Nth boundary interface and its upper adjacent (N-1)th boundary interface form the Nth material partition; the Nth material partition inherits the physical and mechanical properties of the Nth ground surface object;

[0028] The bottom boundary surface and the Nth boundary interface enclose the N+1th material zone; the N+1th material zone inherits the physical and mechanical properties of the underlying rock mass carried by the Nth stratum object according to spatial relationships; the bottom boundary surface passes through at least one material zone;

[0029] Therefore, the three-dimensional design model of the slope, which includes the surface boundary, N layers of ground objects and the bottom boundary, is divided into N+1 material partitions, and each material partition is assigned a physical and mechanical index attribute, thereby forming the slope numerical calculation model.

[0030] Furthermore, in step S4, when the numerical calculation simulation is a static calculation simulation, the boundary conditions are set according to the geological classification attributes of the boundary surfaces. Specifically, the surface boundary surface and the excavation face boundary surface are set as free boundaries; the lateral boundary surface and the bottom boundary surface are set as normal displacement constraint boundary conditions with zero displacement.

[0031] When the numerical simulation of a dynamic working condition is used, two cases are distinguished:

[0032] If the numerical calculation uses the quasi-static method to simulate the dynamic working condition, the boundary conditions are set in the same way as those for the static working condition.

[0033] If the numerical calculation uses the dynamic time history method to simulate dynamic working conditions, the surface boundary and the excavation face boundary are set as free boundaries; the lateral boundary is selected as a viscous boundary or a free field boundary type of wave-absorbing boundary condition; the bottom boundary is set according to the physical and mechanical properties of the bottom layer object, and is a stress boundary condition when the bottom layer is soil, and a velocity boundary condition when the bottom layer is rock.

[0034] Furthermore, the bottom boundary surface is set according to the physical and mechanical properties of the lowest ground layer. When the lowest layer is soil, it is a stress boundary condition; when the lowest layer is rock, it is a velocity boundary condition. Specifically:

[0035] Step S41: Read the physical and mechanical properties of the geological body from the material partition containing the bottom boundary surface to obtain the deformation modulus E of the bottom boundary surface;

[0036] If there are multiple material partitions including the bottom boundary surface, the deformation modulus E is obtained by weighted averaging using formula (1):

[0037] (1)

[0038] in: This represents the total number of material partitions, including the bottom boundary surface. Assign material zone numbering; To correspond to material partitioning in the three-dimensional slope numerical calculation model The area of ​​the bottom boundary surface; Material partitioning The corresponding deformation modulus, The deformation modulus after weighted average;

[0039] Step S42: Determine the geological properties of the bottom boundary surface based on the deformation modulus E: if the deformation modulus E is less than 1 GPa, the bottom boundary surface is a soil layer; otherwise, it is a rock layer.

[0040] Step S43, when the bottom layer is soil, the method for applying the stress boundary condition is as follows: based on the soil layer parameters, the velocity time history curve corresponding to the velocity boundary is converted into a stress time history curve; the conversion method is as follows:

[0041] Step S431: Using formulas (2) and (3), the propagation velocity of the compression wave in the soil layer is obtained respectively. and the propagation speed of shear waves :

[0042] (2)

[0043] (3)

[0044] Where: K is the bulk modulus of the soil layer; G is the shear modulus of the soil layer; The density is the soil layer density; the bulk modulus K and shear modulus G are obtained by converting the soil layer deformation modulus and Poisson's ratio; when there are multiple material zones including the bottom boundary surface, the bulk modulus, shear modulus or density of each material zone are weighted by the area of ​​the bottom boundary surface of each material zone to obtain the bulk modulus, shear modulus or density.

[0045] Step S432, based on the propagation velocity of the compression wave in the soil layer Using formula (4), the normal stress applied to the bottom boundary surface is obtained. Based on the propagation velocity of shear waves in the soil layer Using formula (5), the tangential stress applied to the bottom boundary surface is obtained. :

[0046] (4)

[0047] (5)

[0048] in: and , which are the normal velocity and the tangential velocity, respectively.

[0049] Furthermore, it also includes:

[0050] Step S5: In the slope numerical calculation model containing boundary conditions, based on the structural characteristics of each actual adit in the target slope layout, a real adit model corresponding to each actual adit and considering the spatial variability of the structural surface is generated; the structural surface considering spatial variability inside the real adit model is called a deterministic structural surface;

[0051] Step S6: Taking into account the geometric features and deterministic structural surface features of the real embankment model, several virtual embankment models are generated in other blank areas of the slope numerical calculation model. These virtual embankment models consider both the spatial variability of the structural surfaces and the random changes in the continuity of the structural surfaces along the elevation. The structural surfaces within these virtual embankment models that consider both spatial variability and random changes in the continuity of the structural surfaces along the elevation are called random structural surfaces.

[0052] Step S7: By creating the real adit model and the virtual adit model in the slope numerical calculation model containing boundary conditions, a structural surface network model that conforms to the geological reality of the target slope is obtained, which is used to evaluate the stability of the target slope.

[0053] Furthermore, step S5 includes:

[0054] Step S51: In the slope numerical calculation model, at the elevation position corresponding to each actual tunnel, generate the actual tunnel model with the corresponding tunnel depth.

[0055] Step S52: Within the real-world tunnel model, generate deterministic structural surfaces that consider the spatial variability of structural surfaces. The method is as follows:

[0056] Based on the characteristics of the structural surfaces of the corresponding real tunnel, the characteristics of the deterministic structural surfaces inside the real tunnel model are determined, so that the characteristics of the deterministic structural surfaces are the same as the characteristics of the structural surfaces of the real tunnel, thereby obtaining the spacing between the deterministic structural surfaces and the trace length, dip, inclination angle and position of each deterministic structural surface;

[0057] The deterministic structural surface is subjected to spatial variability processing. While keeping the trace length, dip direction, tilt angle and structural surface position of the deterministic structural surface fixed, the deterministic structural surface is randomly fluctuated up and down along the dip direction. That is, the position of the intersection point with the actual tunnel axis along the trace length of the deterministic structural surface is randomly adjusted. In other words, the starting position and ending position of the deterministic structural surface are randomly adjusted.

[0058] Based on the deterministic structural plane spacing, and the trace length, dip direction, dip angle, structural plane position, start position, and end position of each deterministic structural plane, the final deterministic structural plane is generated, including:

[0059] Define the starting point of the deterministic structural surface and the end point The endpoint closest to the entrance of the real tunnel model is taken as the starting point. The endpoint furthest from the entrance of the real tunnel model is the end point. ;

[0060] If the determinate structural surface dips downhill, then formulas (6) and (7) are used to determine its starting point, respectively. Location and endpoint The location; if the inclination of the deterministic structural surface is reverse slope, then formulas (8) and (9) are used to determine its starting point respectively. Location and endpoint Location;

[0061] (6)

[0062] (7)

[0063] (8)

[0064] (9)

[0065] in: The location of the deterministic structural surface specifically refers to the intersection of the deterministic structural surface and the actual tunnel axis model, located on the axis of the actual tunnel model, and measured from the entrance of the actual tunnel model. and the elevation of that location The elevation is the same as that of the actual tunnel model. The trace length of the deterministic structural surface; The inclination angle of the deterministic structural surface; Starting point The location, through the starting point Distance from the entrance of the real tunnel model and the starting point elevation Characterization; End point The location, through the endpoint Distance from the entrance of the real tunnel model and the finish line elevation Characterization; It is a random number between 0 and 1;

[0066] The starting point of the deterministic structural surface along the slope is determined using formulas (1) and (2). Location and endpoint When determining the location, the starting point must be satisfied. Location to destination The distance from the position is equal to the trace length. The constraints; similarly, when using formulas (3) and (4), the trace length must be satisfied. Constraints.

[0067] Furthermore, step S6 includes:

[0068] A virtual tunnel model is generated between two adjacent real tunnel models: Between each pair of adjacent real tunnel models, considering the geometric features and deterministic structural surface features of the two real tunnel models, several virtual tunnel models are generated that consider both the spatial variability of the structural surface and the random variation of the structural surface along the elevation; the structural surface inside the virtual tunnel model that considers both spatial variability and random variation along the elevation is called a random structural surface;

[0069] The virtual adit model is generated in the slope edge region: In the target slope model, in the upper edge region of the real adit model located at the maximum elevation, and / or in the lower edge region of the real adit model located at the minimum elevation, considering the geometric features and deterministic structural surface features of the nearest real adit model, several virtual adit models are generated that consider the spatial variability of the structural surface and the random variation of the structural surface along the elevation continuity.

[0070] The present invention also provides a slope numerical calculation model construction system based on a three-dimensional geological model with attributes. The slope numerical calculation model construction system based on a three-dimensional geological model with attributes is used to implement the slope numerical calculation model construction method based on a three-dimensional geological model with attributes. The slope numerical calculation model construction system based on a three-dimensional geological model with attributes is a three-dimensional visualization software system.

[0071] The present invention provides a method and system for constructing a slope numerical calculation model based on a three-dimensional geological model with attributes, which has the following advantages:

[0072] (1) This invention is an extension of the three-dimensional digital design of slopes and one of the specific implementation links of intelligent cloud computing. This invention realizes two aspects: three-dimensional slope contour design and evaluation of the engineering safety of the design results. This invention establishes the correlation between the designed excavation contour and the calculation model. Specifically, it assigns classification attributes to the excavation surface, which are inherited by the calculation model and naturally serve as one of the boundaries of the calculation model, providing an identifier for the intelligent setting of boundary conditions. This patent focuses on the design process of excavated slopes and the intelligent setting of boundary conditions of the calculation model, which is a deepening of this patent.

[0073] (2) Based on the constructed slope numerical calculation model, this invention also proposes a method for constructing a slope structural surface model based on a virtual adit. By using actual adit data from the field, the structural surface characteristics of blank areas where adit exploration has not been carried out are restored as realistically as possible, establishing a structural surface network model that conforms to geological reality for subsequent slope stability assessment. This invention ensures that the structural surface network constructed in the model conforms to the actual adit exposure conditions, while also considering the spatial variability of structural surface distribution characteristics in different parts of the slope, thus solving the problem of the structural surface model deviating from reality in the numerical calculation of rock slopes affecting slope stability assessment. Attached Figure Description

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

[0075] Figure 1 A flowchart illustrating a method for constructing a slope numerical calculation model based on a three-dimensional geological model with attributes, provided by this invention;

[0076] Figure 2 This is a schematic diagram of a three-dimensional geological model with attributes provided in an embodiment of the present invention;

[0077] Figure 3 This is a schematic diagram of a design model for establishing slope excavation contour classification attributes based on a three-dimensional geological model, provided in an embodiment of the present invention.

[0078] Figure 4 A schematic diagram illustrating the generation of a three-dimensional numerical calculation model based on a design model, provided in an embodiment of the present invention;

[0079] Figure 5 for Figure 4 A magnified view of a portion of the image;

[0080] Figure 6 A schematic diagram illustrating how the boundary conditions of a computational model are intelligently set based on expert experience according to simulated working conditions, as provided in an embodiment of the present invention.

[0081] Figure 7 A schematic diagram illustrating the composition and functions of a software system provided in an embodiment of the present invention;

[0082] Figure 8 This is a schematic diagram of the real-world tunnel model constructed according to an embodiment of the present invention;

[0083] Figure 9 This is a schematic diagram of a virtual tunnel model constructed between two real tunnel models in an embodiment of the present invention;

[0084] Figure 10 This is a schematic diagram of a virtual adit model located in the lower edge region constructed according to an embodiment of the present invention. Detailed Implementation

[0085] To make the technical problems solved, the technical solutions, and the beneficial effects of this invention clearer, the invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and are not intended to limit the invention.

[0086] In the description of this application, the term "for example" is used to mean "used as an example, illustration, or description." Any embodiment described as "for example" in this application is not necessarily to be construed as being more preferred or advantageous than other embodiments. The following description is provided to enable any person skilled in the art to make and use the invention. Details are set forth in the following description for purposes of explanation. It should be understood that those skilled in the art will recognize that the invention can be made without using these specific details. In other instances, well-known structures and processes will not be described in detail to avoid obscuring the description of the invention with unnecessary detail. Therefore, the invention is not intended to be limited to the embodiments shown, but is consistent with the broadest scope of the principles and features disclosed in this application.

[0087] This invention provides a method and system for constructing a slope numerical calculation model based on a three-dimensional geological model with attributes, which has the following advantages:

[0088] (1) This invention is an extension of the three-dimensional digital design of slopes and one of the specific implementation links of intelligent cloud computing. This invention realizes two aspects: three-dimensional slope contour design and evaluation of the engineering safety of the design results. This invention establishes the correlation between the designed excavation contour and the calculation model. Specifically, it assigns classification attributes to the excavation surface, which are inherited by the calculation model and naturally serve as one of the boundaries of the calculation model, providing an identifier for the intelligent setting of boundary conditions. This patent focuses on the design process of excavated slopes and the intelligent setting of boundary conditions of the calculation model, which is a deepening of this patent.

[0089] (2) Based on the constructed slope numerical calculation model, this invention also proposes a method for constructing a slope structural surface model based on a virtual adit. By using actual adit data from the field, the structural surface characteristics of blank areas where adit exploration has not been carried out are restored as realistically as possible, establishing a structural surface network model that conforms to geological reality for subsequent slope stability assessment. This invention ensures that the structural surface network constructed in the model conforms to the actual adit exposure conditions, while also considering the spatial variability of structural surface distribution characteristics in different parts of the slope, thus solving the problem of the structural surface model deviating from reality in the numerical calculation of rock slopes affecting slope stability assessment.

[0090] See Figure 1 This invention provides a method for constructing a slope numerical calculation model based on a three-dimensional geological model with attributes, comprising:

[0091] Step S1: Construct a three-dimensional geological model with attributes. The three-dimensional geological model with attributes includes a three-dimensional geological model and model attributes. The three-dimensional geological model includes a surface object at the top, a groundwater level object at the bottom, and several surface objects between the surface object and the groundwater level object. The model attributes include the geological classification attributes of each surface object and the physical and mechanical index attributes of each surface object. The geological classification attributes of the surface object include attributes such as surface, strata, faults, and groundwater. The strata can be subdivided according to the degree of weathering, the degree of unloading, and the specific rock and soil type. The physical and mechanical index attributes of each surface object include, but are not limited to, the deformation modulus, Poisson's ratio, and strength parameters of the rock and soil.

[0092] Step S2: According to the engineering requirements, in the attributed three-dimensional geological model, generate the slope excavation outline surface and arrange the slope reinforcement components according to the slope reinforcement scheme. Assign the geological classification attributes of the excavation surface object to the slope excavation outline surface, thereby establishing a three-dimensional slope design model containing geological conditions.

[0093] This step is specifically as follows:

[0094] Step S21: In the attribute-containing three-dimensional geological model, a slope excavation start outline is set. Based on the slope excavation start outline, slope excavation outlines with different slope ratios and slope heights are set from near to far. The slope excavation start outline and each of the slope excavation outlines are combined to form the slope excavation outline surface.

[0095] Therefore, this step uses the slope initiation line (slope excavation start outline line) and parameterized settings such as slope height and slope ratio to quickly establish the slope excavation outline surface, and assigns the geological classification attributes of the "excavation surface object" in subsequent steps.

[0096] Step S22: According to the slope reinforcement scheme, select a reinforcement zone in the slope excavation outline; in the reinforcement zone, arrange the slope reinforcement components according to the set spacing and row spacing of the slope reinforcement components.

[0097] Therefore, in this step, on the slope excavation outline, according to the reinforcement scheme, the reinforcement zone is selected, and the arrangement of reinforcement components such as anchor cables and anchor rods is quickly completed by setting parameters such as spacing and row spacing.

[0098] Step S23: Assign geological classification attributes to the slope excavation outline surface; the slope excavation outline surface and the slope reinforcement components form the three-dimensional design model of the slope containing geological conditions.

[0099] Step S3: Using the closed intersection algorithm, the slope three-dimensional design model containing geological conditions is cut out from the attribute-containing three-dimensional geological model, and the cut slope three-dimensional design model containing geological conditions is parameterized to obtain the slope numerical calculation model.

[0100] In practice, this step can be implemented by using two-dimensional closed intersection technology to generate a two-dimensional slope numerical calculation model, or by using three-dimensional closed intersection technology to generate a three-dimensional slope numerical calculation model.

[0101] As an example, a two-dimensional closed intersection technique is used to generate a two-dimensional numerical calculation model for the slope, specifically:

[0102] The specified profile is used to intersect the three-dimensional design model of the slope to generate corresponding two-dimensional profile lines; the two-dimensional profile lines inherit the geological classification attributes and physical and mechanical index attributes of the three-dimensional interface in the attribute-containing three-dimensional geological model; for example, the excavation outline profile line still inherits the geological classification attributes of the original "excavation face";

[0103] The two-dimensional profile line, combined with the lateral boundary line and the bottom boundary line, is intersected and trimmed to generate a closed material partition; the trimmed two-dimensional profile line still inherits the geological classification attributes and physical and mechanical index attributes before trimming, and the vertical boundary lines on both sides and the horizontal bottom boundary line are automatically assigned the boundary line classification attributes of "lateral" and "bottom".

[0104] Each closed material partition is associated with the corresponding physical and mechanical index attributes based on the relative spatial position of the geological classification attributes of the surrounding stratigraphic profile lines. If there is a stratigraphic boundary before the lower part of the bottom material partition is trimmed, it is used as the lower boundary and brought into the analysis. If there is no stratigraphic boundary before the lower part of the bottom material partition is trimmed and it is not possible to directly associate it with the physical and mechanical index attributes, the corresponding physical and mechanical index attributes can be manually entered.

[0105] As an example, a three-dimensional slope numerical calculation model is generated using a three-dimensional closed intersection technique, including:

[0106] Step S31: Based on the slope excavation outline surface of the slope three-dimensional design model containing geological conditions and the design requirements, determine the spatial range of the slope three-dimensional design model to be cut out, so that it covers the slope excavation outline surface in the horizontal direction and extends from the ground surface object to the bottom plane that meets the design requirements in the vertical direction.

[0107] Step S32: Using the three-dimensional closed intersection technique, the three-dimensional slope design model designed in the attribute-containing three-dimensional geological model is subjected to three-dimensional closed intersection to determine the four sides closed plane and the bottom closed plane of the slope three-dimensional design model.

[0108] Step S33: Based on the closed planes on all four sides and the closed plane at the bottom, cut out the three-dimensional design model of the slope from the three-dimensional geological model containing attributes;

[0109] Step S34, Set boundary surface properties:

[0110] In the three-dimensional design model of the slope obtained by trimming, the top surface object inherits the original geological classification attributes and is given surface boundary surface attributes; the top excavation surface object inherits the original geological classification attributes and is given excavation surface boundary surface attributes; the surrounding closed planes are given lateral boundary surface attributes; and the bottom closed plane is given bottom boundary surface attributes.

[0111] Therefore, the cut excavation face, surface, strata, faults, etc. still inherit the original geological attribute classification, and the four sides and bottom interfaces are automatically assigned the boundary surface attributes of "lateral" and "bottom".

[0112] By setting each boundary surface, the geometric entity of the three-dimensional design model of the slope is formed;

[0113] Step S35, Material Zoning:

[0114] Suppose that the attributed three-dimensional geological model has N layers of ground objects in the spatial range corresponding to the three-dimensional slope design model, which are from top to bottom: the first layer of ground object, the second layer of ground object, ..., the Nth layer of ground object;

[0115] The first ground surface object and the slope three-dimensional design model intersect to form a first boundary interface; the first boundary interface and the ground surface boundary surface enclose a first material partition; the first material partition inherits the physical and mechanical properties of the first ground surface object;

[0116] The second ground surface object and the slope three-dimensional design model intersect to form a second boundary interface; the second boundary interface and the first boundary interface immediately above it form a second material partition; the second material partition inherits the physical and mechanical properties of the second ground surface object;

[0117] And so on

[0118] The Nth ground surface object and the slope 3D design model intersect to form the Nth boundary interface; the Nth boundary interface and its upper adjacent (N-1)th boundary interface form the Nth material partition; the Nth material partition inherits the physical and mechanical properties of the Nth ground surface object;

[0119] The bottom boundary surface and the Nth boundary interface enclose the N+1th material zone; the N+1th material zone inherits the physical and mechanical properties of the underlying rock mass carried by the Nth stratum object according to spatial relationships; the bottom boundary surface passes through at least one material zone;

[0120] Therefore, the three-dimensional design model of the slope, which includes the surface boundary, N layers of ground objects and the bottom boundary, is divided into N+1 material partitions, and each material partition is assigned a physical and mechanical index attribute, thereby forming the slope numerical calculation model.

[0121] Therefore, when generating a 3D slope numerical calculation model using the 3D closed intersection technique, each closed material partition is associated with its corresponding physical and mechanical properties based on the relative spatial position of the geological classification attributes of the surrounding strata. If a stratigraphic interface exists in the lower part of the bottom material partition before trimming, it is used as the lower interface in the analysis. If no stratigraphic interface exists in the lower part of the bottom material partition before trimming, and physical and mechanical properties cannot be directly associated, then manual input of the corresponding physical and mechanical properties is supported.

[0122] For example, such as Figure 4 and Figure 5 As shown, based on the design model results including geological conditions, a three-dimensional numerical calculation model is generated using three-dimensional closed intersection technology. After the surface, strata, excavation face, and the set four sides and bottom boundaries intersect and are trimmed, the trimmed excavation face, surface, strata, etc. still inherit the original classification attributes. The four sides and bottom interfaces are automatically assigned "lateral" and "bottom" boundary surface attributes, forming four closed material zones. Among them, the first three material zones from top to bottom automatically inherit the index attributes of the strata corresponding to their "bottom", that is, forming material zones representing completely weathered rock layers, strongly weathered rock layers, and weakly weathered rock layers. The bottommost material zone inherits the index attributes of the lower rock mass carried by the upper adjacent strata according to spatial relationship, that is, forming a material zone representing slightly weathered rock layers.

[0123] Step S4: Determine the load conditions based on the numerical simulation working conditions. Based on the load conditions and the model attributes in the attribute-containing three-dimensional geological model, the boundary conditions of the slope numerical calculation model can be intelligently set according to expert experience, and finally the slope numerical calculation model with boundary conditions is constructed.

[0124] In this step, when the numerical calculation simulation is a static simulation, the boundary conditions are set according to the geological classification attributes of the boundary surfaces. Specifically, based on expert experience, the surface boundary surface and the excavation face boundary surface are set as free boundaries; the lateral boundary surface and the bottom boundary surface are set as normal displacement constraint boundary conditions with zero displacement. For example, in a two-dimensional or three-dimensional slope numerical calculation model, when the boundary or interface classification attribute is "surface" and "excavation face", it is set as a free boundary; when it is "lateral" and "bottom", it is set as a normal displacement constraint boundary with zero displacement.

[0125] When the numerical simulation of a dynamic working condition is used, two cases are distinguished:

[0126] If the numerical calculation uses the quasi-static method to simulate dynamic working conditions, according to expert experience, the boundary conditions are set in the same way as those for static working conditions.

[0127] If the numerical calculation uses the dynamic time history method to simulate dynamic working conditions, the surface boundary and the excavation face boundary are set as free boundaries; the lateral boundary can be selected as a viscous boundary or a free field boundary type of absorbing boundary condition; the bottom boundary is set according to the physical and mechanical properties of the bottom layer object, and is a stress boundary condition when the bottom layer is soil, and a velocity boundary condition when the bottom layer is rock.

[0128] The bottom boundary surface is set according to the physical and mechanical properties of the bottommost ground layer. When the bottommost layer is soil, it is a stress boundary condition; when the bottommost layer is rock, it is a velocity boundary condition. Specifically:

[0129] Step S41: Read the physical and mechanical properties of the geological body from the material partition containing the bottom boundary surface to obtain the deformation modulus E of the bottom boundary surface;

[0130] In practical applications, if there are multiple material partitions including the bottom boundary surface, the deformation modulus can be obtained by weighted averaging based on the length corresponding to the bottom of the two-dimensional slope numerical calculation model or the area corresponding to the bottom of the three-dimensional slope numerical calculation model. Taking a three-dimensional slope numerical calculation model as an example:

[0131] If there are multiple material partitions including the bottom boundary surface, the deformation modulus E is obtained by weighted averaging using formula (1):

[0132] (1)

[0133] in: This represents the total number of material partitions, including the bottom boundary surface. Assign material zone numbering; To correspond to material partitioning in the three-dimensional slope numerical calculation model The area of ​​the bottom boundary surface; Material partitioning The corresponding deformation modulus, The deformation modulus after weighted average;

[0134] Step S42: Determine the geological properties of the bottom boundary surface based on the deformation modulus E: if the deformation modulus E is less than 1 GPa, the bottom boundary surface is a soil layer; otherwise, it is a rock layer.

[0135] Step S43, when the bottom layer is soil, the method for applying the stress boundary condition is as follows: based on the soil layer parameters, the velocity time history curve corresponding to the velocity boundary is converted into a stress time history curve; the conversion method is as follows:

[0136] Step S431: Using formulas (2) and (3), the propagation velocity of the compression wave in the soil layer is obtained respectively. and the propagation speed of shear waves :

[0137] (2)

[0138] (3)

[0139] Where: K is the bulk modulus of the soil layer; G is the shear modulus of the soil layer; The density is the soil layer density; the bulk modulus K and shear modulus G are obtained by converting the soil layer deformation modulus and Poisson's ratio; when there are multiple material zones including the bottom boundary surface, the bulk modulus, shear modulus or density of each material zone are weighted by the area of ​​the bottom boundary surface of each material zone to obtain the bulk modulus, shear modulus or density.

[0140] Step S432, based on the propagation velocity of the compression wave in the soil layer Using formula (4), the normal stress applied to the bottom boundary surface is obtained. Based on the propagation velocity of shear waves in the soil layer Using formula (5), the tangential stress applied to the bottom boundary surface is obtained. :

[0141] (4)

[0142] (5)

[0143] in: and , which are the normal velocity and the tangential velocity, respectively.

[0144] Therefore, the applied stress boundary can be transformed into a stress time history curve based on parameters such as soil density, thus converting the velocity time history curve corresponding to the velocity boundary.

[0145] In this invention, when there are multiple material partitions at the bottom of a two-dimensional or three-dimensional slope numerical calculation model, the system can automatically obtain a weighted average based on the length corresponding to the bottom of the two-dimensional slope numerical calculation model or the area corresponding to the bottom of the three-dimensional slope numerical calculation model. The formula is illustrated below: Where n is the total number of material partitions including the bottom boundary. Assign material zone numbering; Corresponding material partitions in the two-dimensional slope numerical calculation model The bottom boundary length corresponds to the material partition in the three-dimensional slope numerical calculation model. The area of ​​the bottom boundary; Density can be selected bulk modulus shear modulus After substituting the weighted average value into the formula, the corresponding stress boundary conditions are obtained by substituting the velocity and stress conversion formulas (2) to (5).

[0146] For steps S1 to S4, an embodiment is described below with reference to the accompanying drawings:

[0147] This embodiment is applicable to numerical calculations in slope design within a three-dimensional geological model, and the intelligent setting of boundary conditions for the numerical calculation model. The method and system consist of five key components, including:

[0148] S100, construct a three-dimensional geological model with attributes, including at least the geological classification attributes of geological bodies such as strata and their corresponding physical and mechanical properties; the geological classification attributes refer to the geological classification of geological bodies in the model, such as surface, strata, faults, groundwater, etc.; the physical and mechanical properties include the physical and mechanical property values ​​of strata and faults, such as the deformation modulus, Poisson's ratio and strength parameters of the soil and rock bodies of each stratum.

[0149] like Figure 2 As shown, a three-dimensional geological model with attributes is constructed. The geological interfaces and their corresponding attribute correspondences are as follows: the surface classification attribute is "surface", the classification attribute of completely weathered, strongly weathered, and weakly weathered surfaces is "strata", and the groundwater level is "groundwater" by default; among them, the strata and their corresponding index attributes are: 1) The completely weathered surface is associated with the completely weathered rock mass: the index attribute is density 1900 kg / m³ 3 1) Deformation modulus 20 MPa, Poisson's ratio 0.4, cohesion 20 kPa, friction angle 30°; 2) Strongly weathered surface associated with strongly weathered rock mass, with index properties of density 2550 kg / m³. 3 1) Deformation modulus 1.5 GPa, Poisson's ratio 0.32, cohesion 300 kPa, friction angle 32°; 3) Weakly weathered rock mass above the weakly weathered surface: index property is density 2650 kg / m³3 4) Deformation modulus 6 GPa, Poisson's ratio 0.30, cohesion 550 kPa, friction angle 35°; 5) Slightly weathered rock mass below the weakly weathered surface: index property is density 2720 kg / m³ 3 Deformation modulus 9 GPa, Poisson's ratio 0.28, cohesion 700 kPa, friction angle 39°.

[0150] S200: Based on the engineering requirements, the slope excavation outline and reinforcement scheme design are completed on the basis of the three-dimensional geological model. Geological classification attributes, such as excavation surface, are assigned to the slope excavation outline to obtain the three-dimensional design model of the slope.

[0151] like Figure 3 As shown, based on the three-dimensional geological model, the slope excavation outline is generated by setting parameters such as slope height and slope gradient from bottom to top through the slope line, and a three-dimensional design model of the slope is established; and the designed slope excavation outline is assigned the geological classification attribute of "excavation surface".

[0152] S300, based on the three-dimensional slope design model that includes geological conditions, uses two-dimensional and three-dimensional closed intersection technology to generate a two-dimensional or three-dimensional slope numerical calculation model that includes the boundary.

[0153] Specifically, the surface and excavation face are the components of the boundary of the slope numerical calculation model; the two-dimensional or three-dimensional region formed by the closed intersection between the ground surface and the model boundary is the material partition in the slope numerical calculation model; each region inherits the geological classification attributes of the three-dimensional geological model and is associated with the index attributes, thereby obtaining the physical and mechanical index values ​​of each stratum.

[0154] like Figure 4 and Figure 5 As shown, based on the three-dimensional design model of the slope including geological conditions, a three-dimensional closed intersection technology is used to generate a three-dimensional slope numerical calculation model. After the surface, strata, excavation face, and the set four sides and bottom boundaries are cut and intersected, the cut excavation face, surface, strata, etc. still inherit the original geological classification attributes. The four sides and bottom interfaces are automatically assigned the boundary surface attributes of "lateral" and "bottom". Four closed material zones are formed. Among them, the first three material zones from top to bottom automatically inherit the index attributes of the stratum corresponding to their "bottom", that is, they form material zones representing completely weathered rock layers, strongly weathered rock layers, and weakly weathered rock layers. The bottommost material zone inherits the index attributes of the lower rock mass carried by the upper adjacent stratum according to the spatial relationship, that is, it forms a material zone representing slightly weathered rock layers.

[0155] S400 intelligently sets boundary conditions based on the load conditions corresponding to the numerical simulation working conditions and the model attributes in the three-dimensional geological model, using expert experience.

[0156] like Figure 6As shown, when the numerical simulation condition is set to static condition, the geological classification attributes and index attributes are automatically queried. The three-dimensional interfaces of the geological classification attributes "surface" and "excavation face" at the top of the slope numerical calculation model are given free boundaries, and the three-dimensional interfaces of the geological classification attributes "lateral" and "bottom" are given normal constraint boundaries, with zero displacement.

[0157] When the quasi-static method is used to simulate dynamic working conditions, the boundary conditions are consistent with those of static working conditions.

[0158] When using the dynamic time history method to simulate dynamic conditions, the geological classification attributes at the top of the slope numerical calculation model are...

[0159] The three-dimensional interface of "surface" and "excavation face" is given a free boundary. The bottom material zone has only one partition and the deformation modulus is greater than 1 GPa. Based on expert experience, it is classified as rock mass. The three-dimensional interface of "bottom" is given a default applied velocity boundary. The three-dimensional interface of the classification attribute "lateral" can choose viscous boundary or free field boundary, etc., to absorb the wave.

[0160] Specifically, the bottom boundary (classified as 'bottom') of the slope numerical calculation model is set according to the physical and mechanical parameters of the lowest stratum. If the lowest stratum is soil, it is set as a stress boundary; if it is rock, it is set as a velocity boundary. The standard for classifying soil and rock strata is, in principle, that a deformation modulus less than 1 GPa is considered soil, and a modulus exceeding 1 GPa is considered rock.

[0161] S500, the above operation steps are implemented through a three-dimensional visualization software system. The entire process is based on the three-dimensional visualization software system and is completed through graphical operations, including intelligent setting of boundary conditions of the calculation model that relies on attribute information in the digital model and expert experience.

[0162] The above steps are implemented through a software system, developed using 3D graphics technology, such as... Figure 7 As shown, the software system of this invention is developed using WebGL or OpenGL technology, and can be PC software, B / S architecture, or C / S architecture system, suitable for various engineering scenarios with different levels of confidentiality, and supports the following operations:

[0163] It supports importing 3D geological models with attributes, and maintains the association between the model and the attributes after import.

[0164] It supports the direct design of slope excavation outlines based on three-dimensional geological models, and automatically assigns geological classification attribute values ​​to the excavation outline surfaces;

[0165] It supports two-dimensional or three-dimensional cutting of slope three-dimensional design models containing geological information. The profile boundary line formed by two-dimensional cutting or the external surface (side and bottom) of the computational domain formed by three-dimensional cutting is automatically assigned the classification attribute of 'lateral' or 'bottom'.

[0166] After the cut is closed, the surface, strata, excavation face and other surface objects intersect with the model boundary (lateral, bottom) formed by the cut, dividing the computational domain into multiple closed material partitions, each of which is associated with a corresponding index attribute.

[0167] In the visualization environment, various command operations can be performed through simple point selection and input. Different boundaries can be visualized through color and text annotations. The visualization method helps to verify and correct the settings of each boundary surface.

[0168] See Figures 8-10 This invention further proposes a method for constructing a slope structural surface model based on a virtual adit, building upon the existing numerical calculation model of the slope. Based on actual adit data, the method aims to realistically reconstruct the structural surface characteristics of unexplored areas as much as possible, establishing a geologically accurate structural surface network model for subsequent slope stability assessment. This invention ensures that the structural surface network constructed in the model conforms to the actual adit exposure conditions, while also considering the spatial variability of structural surface distribution characteristics in different parts of the slope. It solves the problem of the structural surface model deviating from reality in numerical calculations of rock slopes, which affects slope stability assessment. See steps S5-S7 for details.

[0169] Step S5: In the slope numerical calculation model containing boundary conditions, based on the structural characteristics of each actual adit in the target slope layout, a real adit model corresponding to each actual adit and considering the spatial variability of the structural surface is generated; the structural surface considering spatial variability inside the real adit model is called a deterministic structural surface;

[0170] In this step, the structural features of each actual adit in the target slope layout include the geometric features of the actual adit and the features of the structural surfaces located inside the actual adit; the geometric features of the actual adit include the depth and elevation; the features of the structural surfaces inside the actual adit include the spacing between structural surfaces and the trace length, dip direction, dip angle, and position of each structural surface; the position of the structural surface refers to the location of the intersection of the structural surface and the axis of the actual adit on the axis of the actual adit; the dip direction includes downslope and upslope.

[0171] This step specifically includes steps S51 to S52:

[0172] Step S51: In the slope numerical calculation model, at the elevation position corresponding to each actual tunnel, generate the actual tunnel model with the corresponding tunnel depth.

[0173] Step S52: Within the real-world tunnel model, generate deterministic structural surfaces that consider the spatial variability of structural surfaces. The method is as follows:

[0174] Based on the characteristics of the structural surfaces of the corresponding real tunnel, the characteristics of the deterministic structural surfaces inside the real tunnel model are determined, so that the characteristics of the deterministic structural surfaces are the same as the characteristics of the structural surfaces of the real tunnel, thereby obtaining the spacing between the deterministic structural surfaces and the trace length, dip, inclination angle and position of each deterministic structural surface;

[0175] The deterministic structural surface is subjected to spatial variability processing. While keeping the trace length, dip direction, tilt angle and structural surface position of the deterministic structural surface fixed, the deterministic structural surface is randomly fluctuated up and down along the dip direction. That is, the position of the intersection point with the actual tunnel axis along the trace length of the deterministic structural surface is randomly adjusted. In other words, the starting position and ending position of the deterministic structural surface are randomly adjusted.

[0176] Based on the deterministic structural plane spacing, and the trace length, dip direction, dip angle, structural plane position, start position, and end position of each deterministic structural plane, the final deterministic structural plane is generated, including:

[0177] Define the starting point of the deterministic structural surface and the end point The endpoint closest to the entrance of the real tunnel model is taken as the starting point. The endpoint furthest from the entrance of the real tunnel model is the end point. ;

[0178] If the determinate structural surface dips downhill, then formulas (6) and (7) are used to determine its starting point, respectively. Location and endpoint The location; if the inclination of the deterministic structural surface is reverse slope, then formulas (8) and (9) are used to determine its starting point respectively. Location and endpoint Location;

[0179] (6)

[0180] (7)

[0181] (8)

[0182] (9)

[0183] in: The location of the deterministic structural surface specifically refers to the intersection of the deterministic structural surface and the actual tunnel axis model, located on the axis of the actual tunnel model, and measured from the entrance of the actual tunnel model. and the elevation of that location The elevation is the same as that of the actual tunnel model. The trace length of the deterministic structural surface; The inclination angle of the deterministic structural surface; Starting point The location, through the starting point Distance from the entrance of the real tunnel model and the starting point elevation Characterization; End point The location, through the endpoint Distance from the entrance of the real tunnel model and the finish line elevation Characterization; It is a random number between 0 and 1;

[0184] The starting point of the deterministic structural surface along the slope is determined using formulas (1) and (2). Location and endpoint When determining the location, the starting point must be satisfied. Location to destination The distance from the position is equal to the trace length. The constraints; similarly, when using formulas (3) and (4), the trace length must be satisfied. Constraints.

[0185] The deterministic structural surface generated using this method has the same trace length, dip direction, dip angle, and position as the corresponding real-world tunnel parameters. However, the intersection point between the deterministic structural surface and the real-world tunnel model varies randomly at the deterministic structural surface. For example, for the deterministic structural surface, its starting point... To the finish line The line connecting them is the trace. - Located on the trace - Points in Intersecting with point Q in the axis of the real tunnel model, the point in this invention On the trace - The position in the middle is randomly determined, that is, the trace. - The points intersecting with point Q are randomly determined to generate a deterministic structural surface that takes into account spatial variability. The deterministic structural surface generated by this method can simulate the actual tunnel structure more realistically while strictly satisfying the parameters already revealed in the actual tunnel.

[0186] Step S6: Taking into account the geometric features and deterministic structural surface features of the real embankment model, several virtual embankment models are generated in other blank areas of the slope numerical calculation model. These virtual embankment models consider both the spatial variability of the structural surfaces and the random changes in the continuity of the structural surfaces along the elevation. The structural surfaces within these virtual embankment models that consider both spatial variability and random changes in the continuity of the structural surfaces along the elevation are called random structural surfaces.

[0187] This step distinguishes between two scenarios:

[0188] In the first scenario, the virtual platform model is generated between two adjacent real-world platform models:

[0189] Between each pair of adjacent real-world tunnel models, several virtual tunnel models are generated by comprehensively considering the geometric features and deterministic structural surface features of the two real-world tunnel models. These virtual tunnel models take into account both the spatial variability of the structural surfaces and the random changes in the continuity of the structural surfaces along the elevation. The structural surfaces within the virtual tunnel models that take into account both spatial variability and random changes in the continuity of the structural surfaces along the elevation are called random structural surfaces.

[0190] In the second scenario, the virtual adit model is generated in the slope edge area:

[0191] In the target slope model, in the upper edge region of the actual adit model located at the maximum elevation, and / or in the lower edge region of the actual adit model located at the minimum elevation, considering the geometric features and deterministic structural surface features of the nearest actual adit model, several virtual adit models are generated that consider the spatial variability of the structural surface while also considering the random variation of the structural surface along the elevation continuity.

[0192] In other words, the method for generating a virtual adit model between two adjacent real adit models in this invention is conceptually similar to the method for generating a virtual adit model in the slope edge region, the only difference being the real adit model to be referenced. When generating a virtual adit model between two adjacent real adit models, the geometric features and structural surface statistical features of the two adjacent real adit models are considered; while when generating a virtual adit model in the slope edge region, only the geometric features and structural surface statistical features of the nearest real adit model are considered.

[0193] The two scenarios are described below:

[0194] In the first scenario, the virtual platform model is generated between two adjacent real-world platform models:

[0195] Between each pair of adjacent real-world tunnel models, several virtual tunnel models are generated by comprehensively considering the geometric features and deterministic structural surface features of the two real-world tunnel models. These virtual tunnel models consider both the spatial variability of the structural surfaces and the random changes in the continuity of the structural surfaces along the elevation. The structural surfaces within the virtual tunnel models that consider both spatial variability and random changes in the continuity of the structural surfaces along the elevation are called random structural surfaces.

[0196] Generating the virtual tunnel model between two adjacent real-world tunnel models includes:

[0197] Step A1: Determine the number of virtual tunnel models to be generated between the two real tunnel models;

[0198] For every two adjacent real-world tunnel models, the real-world tunnel model located at the position with the relatively lower elevation is represented as the real-world tunnel model. Its elevation is Its deterministic structural surface has an average trace length of 1000 mm. ;

[0199] The actual tunnel model located at a relatively high elevation is represented as the actual tunnel model. Its elevation is Its deterministic structural surface has an average trace length of 1000 mm. ;

[0200] Solving formula (10) yields the number of virtual tunnel models that satisfy formula (10). :

[0201] (10)

[0202] in: For real-world tunnel models and real-world plane model Elevation difference, = - ; For real-world tunnel models and real-world plane model The mean of the average trace length of the deterministic structural surface. = ;

[0203] Therefore, in this step, the number of virtual tunnel models to be generated is determined based on the relative relationship between the elevation difference between two adjacent real tunnel models and the average trace length of the deterministic structural surface. .

[0204] Step A2: Determine the elevation of each virtual tunnel model to be generated based on the number of virtual tunnel models;

[0205] According to the direction of elevation from small to large, The virtual tunnel models are represented as follows: Virtual tunnel model Virtual tunnel model Virtual tunnel model ;

[0206] Formula (11) is used to determine the virtual tunnel model. elevation ,in, Represents a virtual tunnel model Number:

[0207] (11)

[0208] Step A3: Determine the elevation position of each of the virtual tunnel models and the tunnel depth of the virtual tunnel models;

[0209] The depth of the virtual tunnel model is the same as that of the real tunnel model. and real-world plane model The average depth of the pit;

[0210] Step A4: Based on the depth of the virtual tunnel model, within the virtual tunnel model, while also considering the real tunnel model. and real-world plane model The relevant parameters of the deterministic structural surface are obtained, and the downslope random structural surface algorithm is used to generate each downslope random structural surface; the reverse slope random structural surface algorithm is used to generate each reverse slope random structural surface, thus generating the virtual tunnel model of random structural surfaces formed by the combination of downslope random structural surfaces and reverse slope random structural surfaces.

[0211] Downslope stochastic structure surface algorithm: for real-world adit models and real-world plane model The relevant parameters of the deterministic structural surfaces along the slope are statistically analyzed, and the relevant parameters of each random structural surface along the slope are determined based on the statistical results, including the spacing, trace length, dip angle, starting position and ending position of the random structural surfaces along the slope.

[0212] Reverse slope stochastic structure surface algorithm: for real-world tunnel models and real-world plane model The relevant parameters of the deterministic structural surfaces of the reverse slope direction are statistically analyzed, and the relevant parameters of each random structural surface of the reverse slope direction are determined based on the statistical results, including the spacing, trace length, dip angle, starting position and ending position of the random structural surfaces of the reverse slope direction.

[0213] In this invention, the downslope random structure surface algorithm and the reverse slope random structure surface algorithm share the same concept, the only difference being that the downslope random structure surface algorithm is applicable to real-world tunnel models. and real-world plane model The relevant parameters of the deterministic structural surface along the slope are statistically analyzed and calculated; the stochastic structural surface algorithm along the reverse slope is applied to the real-world adit model. and real-world plane model The relevant parameters of the deterministic structural surface on the reverse slope are statistically analyzed and calculated.

[0214] Therefore, taking the downslope random structure surface algorithm as an example, in generating virtual adit models... When considering random structural surfaces along the slope, including:

[0215] Step A41, determine the spacing of random structural surfaces along the slope:

[0216] Step A411, calculate the real tunnel model. Mean spacing of slope-direction deterministic structural surfaces ; Computational Reality Plane Model Mean spacing of slope-direction deterministic structural surfaces ;

[0217] Based on the average spacing along the slope Mean spacing along slope Using formula (12), the comprehensive spacing along the slope considering the arrangement distance is obtained. :

[0218] (12)

[0219] Therefore, when generating multiple virtual tunnel models between two real tunnel models, the calculated overall spacing of each virtual tunnel model is different due to the different distances relative to the two real tunnel models. This difference in overall spacing reflects the spatial variability of the structural surface distribution characteristics at different parts of the slope.

[0220] Step A412, for the real-world tunnel model and real-world plane model The maximum value of the longitudinal spacing was obtained by statistically analyzing the spacing between deterministic structural surfaces along the slope. Minimum spacing along slope direction ;

[0221] Step A413: Randomly use formulas (13) and (14) to determine the location of the first random structural surface along the slope to the virtual tunnel model. Spacing between entrances :

[0222] (13)

[0223] (14)

[0224] in: It is a random number between 0 and 1;

[0225] Step A414: Randomly use formulas (15) and (16) to determine the location of the second random structural surface along the slope to the virtual tunnel model. Spacing between entrances :

[0226] (15)

[0227] (16)

[0228] And so on, gradually moving away from the virtual tunnel model The direction of the entrance is used to sequentially determine the random structural surfaces along the slope to the virtual tunnel model. The spacing of the entrances continues until the virtual tunnel model is reached. The depth of the cave is used to determine the spacing of the random structural surfaces along the slope;

[0229] Step A42, determine the trace length of each random structural surface along the slope:

[0230] The trace length of each random structural surface along the slope is represented by the actual adit model. and real-world plane model The mean of the average trace length of the deterministic structural surface along the slope;

[0231] Step A43, determine the dip angle of each random structural surface along the slope:

[0232] Step A431, calculate the real-world tunnel model. Mean slope dip angle of deterministic structural surfaces along the slope ; Computational Reality Plane Model Mean slope dip angle of deterministic structural surfaces along the slope ;

[0233] Based on the average slope angle Mean slope angle Using formula (17), the comprehensive slope angle considering the arrangement distance is obtained. :

[0234] (17)

[0235] Step A432, for the real-world tunnel model and real-world plane model The dip angles of deterministic structural surfaces along the slope were statistically analyzed to obtain the maximum dip angle along the slope. Minimum dip angle along slope ;

[0236] Step A433: Randomly use formulas (18) and (19) to determine the dip angle of each random structural surface along the slope:

[0237] (18)

[0238] (19)

[0239] in: It is a random number between 0 and 1;

[0240] Step A44: Determine the starting and ending points of each random structural surface along the slope:

[0241] For each downslope random structural surface, based on the determined spacing between the downslope random structural surfaces, the distance from each downslope random structural surface to the virtual adit model is determined. The distance to the entrance, combined with the elevation of each random structural surface along the slope and the virtual adit model. Since the elevations are the same, the position of each downslope random structure surface is determined. Combining the trace length and dip angle of the determined downslope random structure surface, the spatial variability processing method is used to obtain the starting and ending positions of each downslope random structure surface.

[0242] In the second scenario, the virtual adit model is generated in the slope edge area:

[0243] In the target slope model, in the upper edge region of the actual adit model located at the maximum elevation, and / or in the lower edge region of the actual adit model located at the minimum elevation, considering the geometric features and deterministic structural surface features of the nearest actual adit model, several virtual adit models are generated that consider the spatial variability of the structural surface while also considering the random variation of the structural surface along the elevation continuity.

[0244] Step B1: Determine the elevation of the virtual adit model to be generated in the slope edge area;

[0245] In the upper edge region of the real tunnel model at the maximum elevation, the elevation interval is determined by using the average trace length of the structural surface of the real tunnel model at the maximum elevation as the elevation spacing; for example, the elevation of the real tunnel model at the maximum elevation is... Its structural surface has an average trace length of In the upper edge region, following the direction of the real tunnel model that gradually moves away from the maximum elevation, the elevations of each virtual tunnel model are as follows: , Similarly, the specific number of virtual tunnel models to be placed in the upper edge region is determined by comprehensively considering the height of the upper edge region and the height of the areas requiring special attention within the upper edge region.

[0246] In the lower edge region of the real tunnel model located at the minimum elevation, the elevation interval is determined by using the average trace length of the structural surface of the real tunnel model at the minimum elevation as the elevation spacing; for example, the elevation of the real tunnel model at the minimum elevation is... Its structural surface has an average trace length of In the lower edge region, following the direction of the real tunnel model that gradually moves away from the minimum elevation, the elevations of each virtual tunnel model are as follows: , Similarly, the specific number of virtual tunnel models to be placed in the lower edge region is determined by comprehensively considering the height of the lower edge region and the height of the areas requiring special attention within the lower edge region.

[0247] Step B2: In the slope edge area, determine the depth of the virtual adit model to be generated;

[0248] In the upper edge region of the real tunnel model located at the maximum elevation, the depth of each virtual tunnel model generated is the same as the depth of the real tunnel model at the maximum elevation;

[0249] In the lower edge region of the real tunnel model located at the minimum elevation, the depth of each virtual tunnel model generated is the same as the depth of the real tunnel model at the minimum elevation;

[0250] Step B3: Based on the depth of the virtual tunnel model, within the virtual tunnel model, only the relevant parameters of the deterministic structural surfaces of the nearest real tunnel model are considered. The downslope random structural surface algorithm is used to generate each downslope random structural surface; the reverse slope random structural surface algorithm is used to generate each reverse slope random structural surface. Thus, the virtual tunnel model with random structural surfaces formed by the combination of downslope random structural surfaces and reverse slope random structural surfaces is generated.

[0251] Check whether the random structural surfaces in the virtual adit model meet the basic requirements for random structural surfaces; if not, regenerate the random structural surfaces inside each virtual adit model until the basic requirements for random structural surfaces are met, thereby constructing the slope structural surface model.

[0252] Specifically, checking whether the random structure surface in the virtual tunnel model meets the basic requirements for random structure surfaces involves checking whether the random structure surface in the virtual tunnel model extends into the real tunnel model. If so, it does not meet the basic requirements for random structure surfaces.

[0253] Step S7: By creating the real adit model and the virtual adit model in the slope numerical calculation model containing boundary conditions, a structural surface network model that conforms to the geological reality of the target slope is obtained, which is used to evaluate the stability of the target slope.

[0254] As an example, such as Figure 8 The image shown is a schematic diagram of a realistic adit model with deterministic structural surfaces generated based on existing geological data in this embodiment. Figure 8 In the process, three real-world tunnel models are generated in the direction of gradually increasing elevation, which are represented as: Real-world tunnel model PD1, PD2 and PD3.

[0255] like Figure 9 The diagram illustrates the generation of a virtual tunnel model between two adjacent real-world tunnel models. Figure 9 In this embodiment, a virtual tunnel model VF1 (PD1-PD2) is generated between real tunnel models PD1 and PD2, and a virtual tunnel model VF1 (PD2-PD3) is generated between real tunnel models PD2 and PD3. In this embodiment, calculations show that one virtual tunnel model is generated between real tunnel models PD1 and PD2; similarly, one virtual tunnel model is generated between real tunnel models PD2 and PD3.

[0256] Taking the generation of a virtual tunnel model between real tunnel models PD1 and PD2 as an example, the specific generation method is as follows:

[0257] S11, respectively calculate the depth, elevation and average trace length of the structural surface of the actual tunnel models PD1 and PD2;

[0258] S12, Determine the number of virtual tunnel models:

[0259] when quantity When, the calculation formula is satisfied That is, inserting a virtual tunnel model *1 into PD1 and PD2. The elevation difference between PD1 and PD2; The mean of the average trace lengths of the deterministic structural surfaces of PD1 and PD2;

[0260] S13, Determine the elevation of each virtual tunnel model:

[0261] In this embodiment, the elevation of the virtual tunnel model VF1 (PD1-PD2) is equal to... ;in, The elevation of PD2, The elevation of PD1.

[0262] S14, Generate random structure surfaces in the virtual tunnel model:

[0263] S141, respectively calculate the average dip angle, average spacing, minimum spacing, maximum spacing, minimum dip angle, and maximum dip angle of the deterministic structural surfaces along the slope and against the slope of two adjacent real tunnel models PD1 and PD2;

[0264] S142, calculate the distance between the position of the first downslope random structure surface and the entrance of the virtual adit model, the trace length, inclination angle, starting position and ending position of the first downslope random structure surface; calculate the relevant parameters of each subsequent downslope random structure surface until the depth of the virtual adit model is reached;

[0265] S143, generate the reverse slope random structure surface according to the same design method.

[0266] Taking the method for generating random structural surfaces along the slope as an example, it includes:

[0267] Determine the distance between the location of the first downslope random structural surface and the entrance of the virtual adit model. : or In the formula, The combined spacing along the slope of the actual adit models PD1 and PD2. The value of is determined by considering the distance from the actual tunnel model; in this embodiment, the number of virtual tunnel models inserted between PD1 and PD2 is n=1. . The average spacing along the slope of PD1; This represents the average spacing along the slope of PD2. The maximum spacing along the slope obtained from the statistics of PD1 and PD2; The minimum spacing along the slope obtained from the statistics of PD1 and PD2; It is a random number between 0 and 1.

[0268] Determine the trace length of the first downslope random structure surface, and take the average of the average trace lengths of the downslope deterministic structure surfaces of PD1 and PD2; based on the spatial variability algorithm.

[0269] Determine the dip angle of the first downslope random structure plane: or ; The combined dip angle along the slope for PD1 and PD2 is taken into account, considering the influence of the distance from the actual adit model. For this embodiment, , and These are the mean dip angles along the slope for PD1 and PD2, respectively. The maximum slope dip angle obtained from PD1 and PD2 statistics; The minimum slope dip angle obtained from PD1 and PD2 statistics.

[0270] Calculate the relevant parameters of each subsequent random structural surface along the slope until the depth of the virtual adit model is reached.

[0271] like Figure 10 The image shows a schematic diagram of generating a virtual tunnel model in the lower edge region. Figure 10 In this process, two virtual adit models are generated, designated as virtual adit model VS1 and virtual adit model VS2, respectively, in the direction of decreasing elevation. For this embodiment, the actual adit with the lowest elevation is located near the excavation slope, which is a key area of ​​concern in the project. Therefore, a virtual adit model is generated below it. When generating virtual adit models VS1 and VS2, the average trace length of the structural surface of the nearest actual adit model PD1 is used as the elevation interval to determine the elevation positions of virtual adit models VS1 and VS2. The depth of virtual adit models VS1 and VS2 is the same as the depth of actual adit model PD1. Only the relevant parameters of the deterministic structural surface of actual adit model PD1 are considered to generate the random structural surfaces of virtual adit models VS1 and VS2.

[0272] In summary, the slope numerical calculation model construction method and system based on an attribute-containing three-dimensional geological model provided by this invention has the following advantages:

[0273] (1) This invention enhances the correlation and inheritance between models: Three-dimensional slope design is carried out based on geological models with attributes. Two-dimensional or three-dimensional closed intersection technology is used to generate two-dimensional or three-dimensional slope numerical calculation models, ensuring the correlation and inheritance of classification attributes and index attributes of slope numerical calculation models, providing identification for intelligent setting of boundary conditions, and ensuring the accuracy of boundary information transmission in the intermediate process.

[0274] (2) This invention improves the rationality of boundary condition setting: Based on the attribute information in the final slope numerical calculation model, this invention intelligently sets reasonable boundary conditions according to different working conditions such as static and dynamic working conditions (quasi-static method, dynamic time history method) simulated by numerical calculation, as well as two-dimensional or three-dimensional numerical calculation models, combined with expert experience. Compared with the traditional manual setting based on experience or simply according to specifications, it is more scientific and reasonable. It can more accurately reflect the actual engineering situation, make the numerical calculation results more reliable, provide stronger support for slope design, and effectively avoid calculation errors caused by unreasonable boundary condition setting;

[0275] (3) This invention lowers the technical threshold and improves work efficiency: In the prior art, slope numerical calculation based on a three-dimensional geological model requires high numerical simulation skills from operators. After constructing a three-dimensional geological model with attributes, this invention can quickly complete the slope excavation outline design according to engineering needs, automatically generate a numerical calculation model using closed intersection technology, and intelligently set boundary conditions. The entire process is automated and intelligent, so that engineers do not need to master complex numerical simulation theories and skills. They only need to perform simple graphical operations according to the prompts in the three-dimensional visualization software system to complete the boundary condition setting, which significantly lowers the technical application threshold, reduces manual operation time, and improves work efficiency.

[0276] (4) The structural surface model construction method proposed in this invention strictly meets the actual tunnel exposure data. Compared with the traditional model based on statistical regularity and random generation, it solves the problem that the structural surface simulated in the actual tunnel location does not match the actual situation.

[0277] (5) Based on the actual adit exposure conditions, this invention restores the true structural surface distribution characteristics of areas where adit exploration has not been carried out as much as possible. It takes into account the spatial variability of structural surface distribution characteristics in different parts of the slope and the spatial variability of structural surface distribution characteristics inside the adit model. The structural surface model constructed based on the method of this invention provides more realistic and reliable results when assessing the stability of rock slopes.

[0278] (6) The present invention employs a method of simulating random structural surfaces in a virtual adit model. Compared with randomly generating structural surfaces in blank areas, this method can more directly control the distribution characteristics of structural surfaces at different locations. By using the adjacent real adit model as a constraint condition, the influence of the distance between the virtual adit model and the real adit model on the distribution characteristics of structural surfaces is taken into account, making the generated random structural surfaces more consistent with real geological conditions.

[0279] The above description is only a preferred embodiment of the present invention. It should be noted that for those skilled in the art, several improvements and modifications can be made without departing from the principle of the present invention, and these improvements and modifications should also be considered within the scope of protection of the present invention.

Claims

1. A method for constructing a slope numerical calculation model based on a three-dimensional geological model with attributes, characterized in that, include: Step S1: Construct a three-dimensional geological model with attributes; the three-dimensional geological model with attributes includes a three-dimensional geological model and model attributes; wherein, the three-dimensional geological model includes a surface object at the top, a groundwater level object at the bottom, and at least one layer of surface objects between the surface object and the groundwater level object; the model attributes include the geological classification attributes of each surface object and the physical and mechanical index attributes of each surface object. Step S2: According to the engineering requirements, in the attributed three-dimensional geological model, generate the slope excavation outline surface and arrange the slope reinforcement components according to the slope reinforcement scheme. Assign the geological classification attributes of the excavation surface object to the slope excavation outline surface, thereby establishing a three-dimensional slope design model containing geological conditions. Step S3: Using the closed intersection algorithm, the slope three-dimensional design model containing geological conditions is cut out from the attribute-containing three-dimensional geological model, and the cut slope three-dimensional design model containing geological conditions is parameterized to obtain the slope numerical calculation model. Step S4: Determine the load conditions based on the numerical simulation working conditions; intelligently set the boundary conditions of the slope numerical calculation model based on the load conditions and the model attributes in the attribute-containing three-dimensional geological model; and finally construct the slope numerical calculation model with boundary conditions. Step S5: In the slope numerical calculation model containing boundary conditions, based on the structural characteristics of each actual adit in the target slope layout, a real adit model corresponding to each actual adit and considering the spatial variability of the structural surface is generated; the structural surface considering spatial variability inside the real adit model is called a deterministic structural surface; Step S6: Taking into account the geometric features and deterministic structural surface features of the real embankment model, a virtual embankment model is generated in other blank areas of the slope numerical calculation model. This virtual embankment model considers both the spatial variability of the structural surface and the random variation of the structural surface along the elevation. The structural surface inside the virtual embankment model that considers both spatial variability and random variation along the elevation is called a random structural surface. Step S6 includes: A virtual tunnel model is generated between two adjacent real tunnel models: Between each pair of adjacent real tunnel models, the geometric features and deterministic structural surface features of the two real tunnel models are comprehensively considered to generate a virtual tunnel model that considers both the spatial variability of the structural surface and the random variation of the structural surface along the elevation; the structural surface inside the virtual tunnel model that considers both spatial variability and random variation along the elevation is called a random structural surface; The virtual adit model is generated in the slope edge region: In the target slope model, in the upper edge region of the real adit model located at the maximum elevation, and / or in the lower edge region of the real adit model located at the minimum elevation, the geometric features and deterministic structural surface features of the nearest real adit model are considered to generate a virtual adit model that considers the spatial variability of the structural surface and the random variation of the structural surface along the elevation. Step S7: By creating the real adit model and the virtual adit model in the slope numerical calculation model containing boundary conditions, a structural surface network model that conforms to the geological reality of the target slope is obtained, which is used to evaluate the stability of the target slope.

2. The method for constructing a slope numerical calculation model based on a three-dimensional geological model with attributes, as described in claim 1, is characterized in that... The physical and mechanical properties of each surface object include the deformation modulus, Poisson's ratio, and strength parameters of the soil and rock mass.

3. The method for constructing a slope numerical calculation model based on a three-dimensional geological model with attributes, as described in claim 1, is characterized in that... Step S2 is as follows: Step S21: In the attribute-containing three-dimensional geological model, a slope excavation start outline is set. Based on the slope excavation start outline, slope excavation outlines with different slope ratios and slope heights are set from near to far. The slope excavation start outline and each of the slope excavation outlines are combined to form the slope excavation outline surface. Step S22: Select a reinforcement zone in the slope excavation outline according to the slope reinforcement scheme; In the reinforced zone, the slope reinforcement members are arranged according to the set spacing and row spacing of the slope reinforcement members; Step S23: Assign geological classification attributes to the excavation surface object of the slope excavation outline surface; The slope excavation profile and the slope reinforcement components together form the three-dimensional design model of the slope, which includes the geological conditions.

4. The method for constructing a slope numerical calculation model based on a three-dimensional geological model with attributes, as described in claim 1, is characterized in that... Step S3 specifically involves generating a three-dimensional slope numerical calculation model using a three-dimensional closed intersection technique, including: Step S31: Based on the slope excavation outline surface of the slope three-dimensional design model containing geological conditions and the design requirements, determine the spatial range of the slope three-dimensional design model to be cut out, so that it covers the slope excavation outline surface in the horizontal direction and extends from the ground surface object to the bottom plane that meets the design requirements in the vertical direction. Step S32: Using the three-dimensional closed intersection technique, the three-dimensional slope design model designed in the attribute-containing three-dimensional geological model is subjected to three-dimensional closed intersection to determine the four sides closed plane and the bottom closed plane of the slope three-dimensional design model. Step S33: Based on the closed planes on all four sides and the closed plane at the bottom, cut out the three-dimensional design model of the slope from the three-dimensional geological model containing attributes; Step S34, Set boundary surface properties: In the three-dimensional design model of the slope obtained by trimming, the top surface object inherits the original geological classification attributes and is given surface boundary surface attributes; the top excavation surface object inherits the original geological classification attributes and is given excavation surface boundary surface attributes; the surrounding closed planes are given lateral boundary surface attributes; and the bottom closed plane is given bottom boundary surface attributes. By setting each boundary surface, the geometric entity of the three-dimensional design model of the slope is formed; Step S35, Material Zoning: Suppose that the attributed three-dimensional geological model has N layers of ground objects in the spatial range corresponding to the three-dimensional slope design model, which are from top to bottom: the first layer of ground object, the second layer of ground object, ..., the Nth layer of ground object; The first ground surface object and the slope three-dimensional design model intersect to form a first boundary interface; the first boundary interface and the ground surface boundary surface enclose a first material partition; the first material partition inherits the physical and mechanical properties of the first ground surface object; The second ground surface object and the slope three-dimensional design model intersect to form a second boundary interface; the second boundary interface and the first boundary interface immediately above it form a second material partition; the second material partition inherits the physical and mechanical properties of the second ground surface object; And so on The Nth ground surface object and the slope 3D design model intersect to form the Nth boundary interface; the Nth boundary interface and its upper adjacent (N-1)th boundary interface form the Nth material partition; the Nth material partition inherits the physical and mechanical properties of the Nth ground surface object; The bottom boundary surface and the Nth boundary interface enclose the N+1th material zone; the N+1th material zone inherits the physical and mechanical properties of the underlying rock mass carried by the Nth stratum object according to spatial relationships; the bottom boundary surface passes through at least one material zone; Therefore, the three-dimensional design model of the slope, which includes the surface boundary, N layers of ground objects and the bottom boundary, is divided into N+1 material partitions, and each material partition is assigned a physical and mechanical index attribute, thereby forming the slope numerical calculation model.

5. The method for constructing a slope numerical calculation model based on a three-dimensional geological model with attributes according to claim 1, characterized in that, In step S4, when the numerical calculation simulation is a static calculation simulation, the boundary conditions are set according to the geological classification attributes of the boundary surfaces. Specifically, the surface boundary surface and the excavation face boundary surface are set as free boundaries; the lateral boundary surface and the bottom boundary surface are set as normal displacement constraint boundary conditions with zero displacement. When the numerical simulation of a dynamic working condition is used, two cases are distinguished: If the numerical calculation uses the quasi-static method to simulate the dynamic working condition, the boundary conditions are set in the same way as those for the static working condition. If the numerical calculation uses the dynamic time history method to simulate dynamic working conditions, the surface boundary and the excavation face boundary are set as free boundaries; the lateral boundary is selected as a viscous boundary or a free field boundary type of wave-absorbing boundary condition; the bottom boundary is set according to the physical and mechanical properties of the bottom layer object, and is a stress boundary condition when the bottom layer is soil, and a velocity boundary condition when the bottom layer is rock.

6. The method for constructing a slope numerical calculation model based on a three-dimensional geological model with attributes, as described in claim 5, is characterized in that... The bottom boundary surface is set according to the physical and mechanical properties of the bottommost ground layer. When the bottommost layer is soil, it is a stress boundary condition; when the bottommost layer is rock, it is a velocity boundary condition. Specifically: Step S41: Read the physical and mechanical properties of the geological body from the material partition containing the bottom boundary surface to obtain the deformation modulus E of the bottom boundary surface; If there are multiple material partitions including the bottom boundary surface, the deformation modulus E is obtained by weighted averaging using formula (1): （1） in: This represents the total number of material partitions, including the bottom boundary surface. Assign material zone numbers; To correspond to material partitioning in the three-dimensional slope numerical calculation model The area of ​​the bottom boundary surface; Material partitioning The corresponding deformation modulus, The deformation modulus after weighted average; Step S42: Determine the geological properties of the bottom boundary surface based on the deformation modulus E: if the deformation modulus E is less than 1 GPa, the bottom boundary surface is a soil layer; otherwise, it is a rock layer. Step S43, when the bottom layer is soil, the method for applying the stress boundary condition is as follows: based on the soil layer parameters, the velocity time history curve corresponding to the velocity boundary is converted into a stress time history curve; the conversion method is as follows: Step S431: Using formulas (2) and (3), the propagation velocity of the compression wave in the soil layer is obtained respectively. and the propagation speed of shear waves : （2） （3） Where: K is the bulk modulus of the soil layer; G is the shear modulus of the soil layer; The density is the soil layer density; the bulk modulus K and shear modulus G are obtained by converting the soil layer deformation modulus and Poisson's ratio; when there are multiple material zones including the bottom boundary surface, the bulk modulus, shear modulus or density of each material zone are weighted by the area of ​​the bottom boundary surface of each material zone to obtain the bulk modulus, shear modulus or density. Step S432, based on the propagation velocity of the compression wave in the soil layer Using formula (4), the normal stress applied to the bottom boundary surface is obtained. Based on the propagation velocity of shear waves in the soil layer Using formula (5), the tangential stress applied to the bottom boundary surface is obtained. : （4） （5） in: and , which are the normal velocity and the tangential velocity, respectively.

7. The method for constructing a slope numerical calculation model based on a three-dimensional geological model with attributes according to claim 1, characterized in that, Step S5 includes: Step S51: In the slope numerical calculation model, at the elevation position corresponding to each actual tunnel, generate the actual tunnel model with the corresponding tunnel depth. Step S52: Within the real-world tunnel model, generate deterministic structural surfaces that consider the spatial variability of structural surfaces. The method is as follows: Based on the characteristics of the structural surfaces of the corresponding real tunnel, the characteristics of the deterministic structural surfaces inside the real tunnel model are determined, so that the characteristics of the deterministic structural surfaces are the same as the characteristics of the structural surfaces of the real tunnel, thereby obtaining the spacing between the deterministic structural surfaces and the trace length, dip, inclination angle and position of each deterministic structural surface; The deterministic structural surface is subjected to spatial variability processing. While keeping the trace length, dip direction, tilt angle and structural surface position of the deterministic structural surface fixed, the deterministic structural surface is randomly fluctuated up and down along the dip direction. That is, the position of the intersection point with the actual tunnel axis along the trace length of the deterministic structural surface is randomly adjusted. In other words, the starting position and ending position of the deterministic structural surface are randomly adjusted. Based on the deterministic structural plane spacing, and the trace length, dip direction, dip angle, structural plane position, start position, and end position of each deterministic structural plane, the final deterministic structural plane is generated, including: Define the starting point of the deterministic structural surface and the end point The endpoint closest to the entrance of the real tunnel model is taken as the starting point. The endpoint furthest from the entrance of the real tunnel model is the end point. ; If the determinate structural surface dips downhill, then formulas (6) and (7) are used to determine its starting point, respectively. Location and endpoint The location; if the inclination of the deterministic structural surface is reverse slope, then formulas (8) and (9) are used to determine its starting point respectively. Location and endpoint Location; （6） （7） （8） （9） in: The location of the deterministic structural surface specifically refers to the intersection of the deterministic structural surface and the actual tunnel axis model, located on the axis of the actual tunnel model, and measured from the entrance of the actual tunnel model. and the elevation of that location The elevation is the same as that of the actual tunnel model. The trace length of the deterministic structural surface; The inclination angle of the deterministic structural surface; Starting point The location, through the starting point Distance from the entrance of the real tunnel model and the starting point elevation Characterization; End point The location, through the endpoint Distance from the entrance of the real tunnel model and the finish line elevation Characterization; It is a random number between 0 and 1; The starting point of the deterministic structural surface along the slope is determined using formulas (1) and (2). Location and endpoint When determining the location, the starting point must be satisfied. Location to destination The distance from the position is equal to the trace length. The constraints; similarly, when using formulas (3) and (4), the trace length must be satisfied. Constraints.

8. A slope numerical calculation model construction system based on an attribute-containing three-dimensional geological model, characterized in that, The slope numerical calculation model construction system based on a three-dimensional geological model with attributes is used to implement the slope numerical calculation model construction method based on a three-dimensional geological model with attributes as described in any one of claims 1-7; the slope numerical calculation model construction system based on a three-dimensional geological model with attributes is a three-dimensional visualization software system.

Images