An intelligent assignment method for rock-soil parameters of a non-continuous numerical model based on coding

Abstract

The present application relates to the technical field of intelligent design of geotechnical engineering, and discloses a method for intelligently assigning geotechnical parameters based on a coded discontinuous numerical model, which comprises the following steps: obtaining parameter values for discontinuous numerical simulation of geotechnical engineering; importing a three-dimensional geological model with attributes and an excavation contour design model to generate an initial block model; converting the initial block model into a discontinuous numerical calculation model; coding each block and the contact surface between blocks of the discontinuous numerical calculation model, with the coding values of each block and each contact surface being unique; simulating joints of the discontinuous numerical calculation model, cutting the blocks in the discontinuous numerical calculation model according to the geometric information of the joints, and synchronously updating the block coding and the contact surface coding; then determining the type of the main body to be assigned, and assigning values according to the type of the main body to be assigned; and the present application achieves intelligent assignment by independently coding the blocks and the contact surfaces in the discontinuous numerical model and associating the coding with geological information.

Description

Technical Field

[0001] This invention relates to the field of intelligent design technology for geotechnical engineering, specifically to an intelligent assignment method for geotechnical parameters based on coding in discontinuous numerical models. Background Technology

[0002] Numerical simulation technology, first applied to geotechnical engineering practice in the 1980s, is the earliest high-end digital technology put into production application in the field, primarily serving as a research tool for safety evaluation of geotechnical engineering design schemes. For a long time, due to the lagging digitalization levels in geological exploration and geotechnical design, digital numerical simulation technology lacked interaction with these processes. With significant improvements in the digitalization levels of these processes, establishing interaction between them can change the application of numerical simulation technology; intelligent cloud computing is the latest development direction. In engineering design and construction, a three-dimensional digital model containing geological and engineering conditions is directly used as the input condition for the numerical model, automatically completing the entire process of modeling, setting, assigning values, calculation, and result analysis without involving other numerical calculation software. This process is defined as intelligent cloud computing. Intelligent parameter assignment is a crucial step. Discontinuous numerical simulation requires input geotechnical parameters including the physical and mechanical parameters of strata, faults, and joint groups. The number of parameters related to the constitutive and strength criteria used in the calculation is large; therefore, intelligent cloud computing often employs database-based categorized storage and management. A three-dimensional discontinuous numerical calculation model can contain tens of thousands of blocks and hundreds of thousands of contact surfaces, and each block and contact surface needs to be assigned different geotechnical parameter values.

[0003] In the prior art, such as the patent document with publication number CN118228451A, an engineering numerical cloud computing method and system are disclosed. It takes advantage of the fact that current numerical simulation software is basically driven by command stream files. Its innovation lies in constructing a command stream file generation and software driving method applicable to different numerical simulation software. However, it does not involve the specific implementation technology of intelligent assignment. When there are a large number of blocks and contact surfaces in the model, it cannot quickly associate with the database.

[0004] Therefore, based on the existing patented technologies mentioned above, it is necessary to design a technical solution that focuses on the basic work required for intelligent assignment of parameters of discontinuous numerical models in the process of intelligent cloud computing services, and focuses on the intelligent parameter assignment process, so as to efficiently associate the database when there are a large number of blocks and contact surfaces in the model. Summary of the Invention

[0005] To address the shortcomings of existing technologies that do not address the specific implementation techniques for intelligent parameter assignment in discontinuous numerical models, and the inability to quickly connect to the database when a large number of blocks and contact surfaces appear in the model, a coding-based intelligent parameter assignment method for discontinuous numerical models is provided, comprising the following steps:

[0006] The parameter values ​​for discontinuous numerical simulation of geotechnical engineering are obtained in three-dimensional geological design software; the parameter values ​​include the physical and mechanical parameters of stratigraphic groups, fault groups, and joint groups;

[0007] Import the attribute-containing 3D geological model and the excavation outline design model to generate an initial block model; convert the initial block model into a discontinuous numerical calculation model;

[0008] Each block and the contact surface between the blocks in the discontinuous numerical calculation model are encoded separately, and the encoded value of each block and each contact surface remains unique;

[0009] The discontinuous numerical calculation model is simulated with joints, and the blocks in the discontinuous numerical calculation model are cut according to the geometric information of the joints, and the block codes and contact surface codes are updated synchronously; the geometric information of the joints includes joint spacing, trace length, connectivity, dip direction and dip angle;

[0010] After completing joint cutting, block and contact surface encoding updates, determine the type of the subject to be assigned a value, and assign a value according to the type of the subject to be assigned a value; the types of the subject to be assigned a value include blocks and contact surfaces;

[0011] When the subject to be assigned is a block, the stratum where the centroid of the subject to be assigned is located is obtained, the code of the subject to be assigned is assigned with the stratum partition code, and then the physical and mechanical parameters of the stratum are assigned.

[0012] When the subject to be assigned is a contact surface, the model information of the contact surface is obtained, the type of the model information is determined, and the corresponding parameters are assigned to the subject to be assigned. The model information of the contact surface includes auxiliary surfaces, excavation contours, strata, faults, and joint information required for modeling.

[0013] Furthermore, when the type of model information is stratigraphy, fault, and joint information, it is associated with the physical and mechanical parameters. After assigning the code of the subject to be assigned with the stratigraphic zone code, fault zone code, and joint zone code, the corresponding physical and mechanical parameters of the stratigraphic group, fault group, and joint group are assigned.

[0014] When the type of model information is auxiliary surface and excavation contour required for modeling, it cannot be associated with the physical and mechanical parameters, so the artificial surface parameters of the subject to be assigned are given.

[0015] Furthermore, the parameter values ​​are stored and managed using a relational database, which is stored and managed in the form of database forms, including stratigraphic lithology forms, fault forms, and joint set forms.

[0016] Furthermore, the stratigraphic lithology form associates stratigraphic physical and mechanical parameters with stratigraphic codes, which include stratigraphic codes, density, deformation modulus, Poisson's ratio, cohesion, and internal friction angle.

[0017] Furthermore, the fault form associates fault physical and mechanical parameters with fault numbers, which include fault number, filling type, tangential stiffness, normal stiffness, cohesion, and friction coefficient.

[0018] Furthermore, the joint group form associates joint physical and mechanical indicators with the joint group number. The joint physical and mechanical indicators include joint group number, dip direction, dip angle, tangential stiffness, normal stiffness, cohesion, and friction coefficient.

[0019] Furthermore, the attribute-containing three-dimensional geological model includes the ground surface, stratigraphic interfaces, and faults, and the excavation profile design model includes the excavation slope profile.

[0020] Furthermore, the contact surface is defined as a common edge or common surface between blocks, and the model information of the contact surface is attached when the contact surface is defined.

[0021] Furthermore, the discontinuous numerical calculation model is generated by transforming the excavation outline, surface, strata, faults, and selected calculation boundaries using a topological closed intersection technique.

[0022] Furthermore, when the discontinuous numerical calculation model is generated using the topological closed intersection technique, a stratigraphic model is also generated based on the surface, strata, and selected calculation boundaries. The stratigraphic model divides the initial block model into 6 material partitions based on the stratigraphic profile lines of the stratigraphic interface, which serve as the basis for determining the stratigraphic parameter values.

[0023] The beneficial effects of this invention are:

[0024] (1) When creating a calculation model based on 3D geological design software, the geological classification information in the model is automatically inherited. It can automatically mark the strata to which each block belongs and the fault or joint group to which each contact surface belongs in the calculation model, and associate it with the geotechnical parameter database to complete the intelligent assignment of parameters. It realizes standardized digital operation + cross-professional data interoperability, forming a dual guarantee of data and system collaboration.

[0025] (2) Independent coding is adopted for blocks and contact surfaces in discontinuous calculation models. The unique code is associated with geological information, and intelligent assignment of geotechnical parameters of the calculation model is realized through relational database. This greatly reduces tedious and repetitive manual operations and reduces the probability of calculation result distortion caused by incorrect input conditions during numerical calculation.

[0026] (3) Based on the geological model with attributes and the design model, the discontinuous numerical calculation model is generated by the cutting and closing technology, and the calculation model construction and parameter assignment are completed quickly. This opens up the key link of intelligent cloud computing and solves the industry problem of three-dimensional interaction between geology, design and calculation. Attached Figure Description

[0027] Figure 1 This is a flowchart of the intelligent assignment method for geotechnical parameters in a non-continuous numerical model based on coding provided by the present invention;

[0028] Figure 2 This is a schematic diagram of a relational database provided by the present invention;

[0029] Figure 3 This is a schematic diagram of the attribute-containing three-dimensional geological model and excavation contour design model provided by the present invention;

[0030] Figure 4 This invention provides a discontinuous numerical calculation model. Figure 4 (a) is the initial block model generated based on the excavation profile, surface, strata, faults, and selected computational boundaries. Figure 4 (b) A stratigraphic model generated based on the surface, strata, and selected computational boundaries;

[0031] Figure 5 This is a schematic diagram of the coding and definition method for blocks / contact surfaces provided by the present invention;

[0032] Figure 6 This is a schematic diagram of the coding update method for blocks and contact surfaces when further cutting using joints, provided by the present invention;

[0033] Figure 7 This is a schematic diagram of the method for determining the value of the subject parameter to be assigned provided by the present invention;

[0034] Explanation of reference numerals in the attached drawings: Stratigraphic designations T1-T6; Fault numbers F1-F2; Joint group numbers J1-J4. Detailed Implementation

[0035] The technical solution of the present invention is further described below, but the scope of protection is not limited to what is described.

[0036] This invention provides an intelligent method for assigning geotechnical parameters to discontinuous numerical models based on coding, such as... Figure 1As shown, it includes the following steps:

[0037] Step S100: Obtain parameter values ​​for discontinuous numerical simulation of geotechnical engineering in three-dimensional geological design software; the parameter values ​​include physical and mechanical parameters of stratigraphic groups, fault groups, and joint groups;

[0038] The parameter values ​​are stored and managed using a relational database, which is stored and managed in the form of database forms. The relational database includes stratigraphic lithology forms, fault forms, and joint set forms, containing stratigraphic lithology, faults, joint sets, and their corresponding physical and mechanical parameters.

[0039] The stratigraphic lithology form associates stratigraphic physical and mechanical parameters with stratigraphic codes. These parameters include stratigraphic code, density, deformation modulus, Poisson's ratio, cohesion, and internal friction angle. The stratigraphic codes are T1-T6, numbered from top to bottom according to the stratigraphic structure.

[0040] The fault form associates fault physical and mechanical parameters with fault numbers. The fault physical and mechanical parameters include fault number, filling type, tangential stiffness, normal stiffness, cohesion, and friction coefficient. The fault numbers are F1-F2, and are numbered according to a preset order. The preset order includes left to right, right to left, top to bottom, and bottom to top.

[0041] The joint group form associates joint physical and mechanical indicators with joint group numbers. The joint physical and mechanical indicators include joint group number, dip direction, dip angle, tangential stiffness, normal stiffness, cohesion, and friction coefficient. The joint group numbers are J1-J4, and are numbered according to the preset order.

[0042] In summary, the relational database includes three main categories: stratigraphy, joint sets, and faults, such as... Figure 2 As shown.

[0043] Step S200: Import the attribute-containing 3D geological model and excavation outline design model to generate an initial block model; convert the initial block model into a discontinuous numerical calculation model; the attribute-containing 3D geological model and excavation outline design model are as follows: Figure 3 As shown.

[0044] Specifically, based on the excavation outline, surface, strata, faults, and selected computational boundaries (lateral and bottom boundaries), the initial block model is divided into a set of several two-dimensional or three-dimensional blocks to generate an initial discontinuous computational model. Furthermore, a stratigraphic model is generated based on the intersection and closure of the surface, strata, and computational boundaries. Each block in the stratigraphic model represents a set of stratigraphic partitions, which are stored separately as the basis for intelligent assignment of stratigraphic parameters.

[0045] The attributed three-dimensional geological model includes the ground surface, stratigraphic interfaces, and faults, and the excavation contour design model includes the excavation slope contour; both are associated with the relational database through attribute codes.

[0046] The discontinuous numerical calculation model is generated by transforming the excavation outline, surface, strata, faults, and selected calculation boundaries using a topological closed intersection technique.

[0047] When generating a discontinuous computational model, common edges (two-dimensional) or common surfaces (three-dimensional) between blocks are defined as contact surfaces. When defining contact surfaces, model information (excavation outline, strata, faults, etc.) is attached.

[0048] When a fault in a 3D geological model does not penetrate the surface or the boundary of the calculation range, an auxiliary surface is added. This auxiliary surface is collinear with the unpenetrated fault location, shares nodes, and penetrates the boundary, thus maintaining the complete and realistic shape of the fault. After adding the auxiliary surface, a new closed intersection is performed to generate blocks and define the contact surface. The auxiliary surface type is added with the information when defining the contact surface.

[0049] A discontinuous numerical computation model of a two-dimensional block set is transformed using a topological closed intersection technique, such as... Figure 4 As shown, Figure 4 (a) is the initial block model generated based on the excavation profile, surface, strata, faults, and selected computational boundaries (lateral and bottom boundaries). Figure 4 (b) is a stratigraphic model generated based on the surface, strata and selected computational boundaries; wherein the stratigraphic model divides the initial block model into 6 material partitions (T1-T6) based on the stratigraphic profile of the stratigraphic interface, which are used to determine the subordinate relationship between the block and the strata in the subsequent process.

[0050] Step S300: Encode each block and the contact surface between the blocks in the discontinuous numerical calculation model, and keep the encoding value of each block and each contact surface unique;

[0051] When generating blocks and contact surfaces using the sectioning and closure technique, each initial block and contact surface is assigned a unique code. The assigned code is fixed and does not change with newly generated blocks and contact surfaces.

[0052] Figure 5 The diagram shows the coding and definition method of blocks / contact surfaces. In this embodiment, each block in the initial discontinuous calculation model is first assigned an independent and unique code as the initial code, and the contact surfaces (i.e. common edges) between blocks are assigned an independent and unique code as the initial code. The contact surface code is accompanied by model information (excavation outline, strata, faults, auxiliary surfaces, etc.).

[0053] Step S400: When the computational model needs to further simulate joints, the discontinuous numerical computational model is simulated with joints. The blocks in the discontinuous numerical computational model are cut according to the geometric information of the joints, and the block codes and contact surface codes are updated synchronously. The joint geometric information includes joint spacing, trace length, connectivity, dip direction, and dip angle.

[0054] like Figure 6 The diagram shows a block / contact surface encoding update method when further cutting using joint groups. The blocks in the calculation model are further cut according to the joint geometry information, and the block encoding and contact surface encoding are updated.

[0055] The joint group dip and dip angle are associated with a database. Spacing, trace length, and connectivity are obtained from survey data and then defined by the user. When the connectivity is less than 100%, it means that the joint surface does not penetrate the boundary. The processing method is the same as in step S200, which involves adding auxiliary surfaces to ensure the accuracy and reliability of the joint connectivity.

[0056] Blocks and contact surfaces unaffected by cutting retain their original codes. Blocks split in two by a single joint are designated as the "old" block, retaining their original code, while the other block is designated as the "new" block, receiving a new code based on its existing code. Contact surfaces split in two by a single joint are handled in the same way. Furthermore, newly generated contact surfaces due to joint cutting receive a new code based on their existing code, and the model information for the contact surface is included when defining it.

[0057] Step S500: After completing joint cutting, block and contact surface encoding updates, determine the type of the subject to be assigned values, and assign values ​​according to the type of the subject to be assigned values; the types of the subject to be assigned values ​​include blocks and contact surfaces; the method for determining the assignment of parameters of the subject to be assigned values ​​is as follows: Figure 7 As shown.

[0058] When the subject to be assigned is a block, the stratum where the centroid of the subject to be assigned is located is obtained, the code of the subject to be assigned is assigned with the stratigraphic partition code, and then the physical and mechanical parameters of the stratigraphic group are assigned.

[0059] In this embodiment, after further cutting using four sets of joint surfaces, the number of blocks in the model is 5810. Specifically, the relative positional relationship between the centroid of the block and the stratigraphic model can be calculated using the ray method. Test rays are emitted along the Z direction, and the intersections of the rays with all stratigraphic contours are detected. The intersections are traversed, and the number of times the stratigraphic region is crossed is recorded. The stratigraphic region crossed an odd number of times is the stratigraphic region to which the centroid of the block belongs. The block is then coded and assigned a stratigraphic zoning code, which is used to associate it with the stratigraphic attributes in the geotechnical parameter database.

[0060] When the subject to be assigned is a contact surface, the model information of the contact surface is obtained, the type of the model information is determined, and the corresponding parameters are assigned to the subject to be assigned. The model information of the contact surface includes auxiliary surfaces, excavation contours, strata, faults, and joint information required for modeling.

[0061] In this embodiment, when the type of model information is strata, fault, and joint information, it is associated with the physical and mechanical parameters. After assigning the code of the subject to be assigned with the strata partition code, fault zone code, and joint zone code, the corresponding physical and mechanical parameters of the strata, fault, and joint group are assigned.

[0062] When the type of model information is auxiliary surface and excavation contour required for modeling, it cannot be associated with the physical and mechanical parameters, so the artificial surface parameters of the subject to be assigned are given.

[0063] In this embodiment, after further cutting using four sets of joint surfaces, the number of contact surfaces in the model is 31,319. Based on the model information attached to the contact surfaces (strata, faults, joints, auxiliary surfaces, excavation outlines), the corresponding contact surfaces are coded and assigned type codes. The fault and joint codes are then associated with the fault and joint groups in the geotechnical parameter database, thereby linking them to the fault and joint attributes in the database. For contact surfaces (auxiliary surfaces, excavation outlines, etc.) that cannot be associated with the geotechnical parameter database, they are assigned the same constitutive model (artificial surface parameters) as faults and joints, increasing the strength parameters to ensure strata continuity.

[0064] The above-disclosed embodiments are merely specific examples of the present invention. However, the present invention is not limited thereto, and any variations that can be conceived by those skilled in the art should fall within the protection scope of the present invention.

Claims

1. A method for intelligent assignment of geotechnical parameters based on a coded non-continuous numerical model, characterized in that, The method comprises the following steps: Obtaining parameter values of geotechnical engineering discontinuous numerical simulation in three-dimensional geological design software; the parameter values include physical and mechanical parameters of stratum group, fault group and joint group; Importing a three-dimensional geological model with attributes and an excavation contour design model to generate an initial block model; converting the initial block model into a discontinuous numerical calculation model; the three-dimensional geological model with attributes includes a ground surface, stratum interfaces and faults; the excavation contour design model includes an excavation slope contour; the three-dimensional geological model with attributes and the excavation contour design model are both associated with a relational database through attribute codes; the discontinuous numerical calculation model is converted and generated by using a topological closed intersection technology based on the excavation contour, the ground surface, the strata, the faults and selected calculation boundaries; Encoding each block and the contact surface between blocks of the discontinuous numerical calculation model respectively; the encoding values of each block and each contact surface remain unique; Simulating joints of the discontinuous numerical calculation model; cutting the blocks in the discontinuous numerical calculation model according to the geometric information of the joints; synchronously updating the block encoding and the contact surface encoding; the joint geometric information includes joint spacing, trace length, connectivity rate, inclination and dip angle; After completing the joint cutting, the block and contact surface encoding updating, judging the type of the subject to be assigned; assigning values according to the type of the subject to be assigned; the type of the subject to be assigned includes blocks and contact surfaces; When the subject to be assigned is a block, obtaining the stratum where the centroid of the subject to be assigned is located; assigning the encoding of the subject to be assigned with the stratum partition code; and assigning the physical and mechanical parameters of the stratum; When the subject to be assigned is a contact surface, obtaining the model information of the contact surface; judging the type of the model information; and assigning corresponding parameters to the subject to be assigned; the model information of the contact surface includes auxiliary surfaces, excavation contours, strata, fault and joint information required for modeling.

2. The intelligent assignment method of rock-soil parameters based on the coded non-continuous numerical model according to claim 1, characterized in that, When the type of the model information is stratum, fault and joint information, associating with the physical and mechanical parameters; assigning the encoding of the subject to be assigned with the stratum partition code, the fault area code and the joint area code; and assigning the corresponding physical and mechanical parameters of the stratum group, the fault group and the joint group; When the type of the model information is auxiliary surfaces and excavation contours required for modeling, the physical and mechanical parameters cannot be associated, and then artificial surface parameters are assigned to the subject to be assigned.

3. The intelligent assignment method of rock-soil parameters based on the coded non-continuous numerical model according to claim 1, characterized in that, The parameter values are stored and managed by using a relational database; the relational database stores and manages in the form of database tables; the relational database includes a stratum lithology table, a fault table and a joint group table.

4. The intelligent assignment method of rock-soil parameters based on the coded non-continuous numerical model according to claim 3, characterized in that, The stratum lithology table is associated with stratum physical and mechanical index parameters through stratum codes; the stratum physical and mechanical index parameters include stratum codes, density, deformation modulus, Poisson's ratio, cohesion and internal friction angle.

5. The intelligent assignment method of rock-soil parameters based on the coded non-continuous numerical model according to claim 3, characterized in that, The fault table is associated with fault physical and mechanical index parameters through fault numbers; the fault physical and mechanical index parameters include fault numbers, filling types, tangential stiffness, normal stiffness, cohesion and friction coefficient.

6. The intelligent assignment method of rock-soil parameters based on the coded non-continuous numerical model according to claim 3, characterized in that, The joint group table is associated with joint physical and mechanical indexes through joint group numbers, and the joint physical and mechanical indexes include joint group numbers, inclination, dip angle, tangent stiffness, normal stiffness, cohesive force and friction coefficient.

7. The intelligent assignment method of rock-soil parameters based on the coded non-continuous numerical model according to claim 1, characterized in that, The contact surface is defined as a common edge or a common surface between blocks, and model information of the contact surface is defined.

8. The intelligent assignment method of rock-soil parameters based on the coded non-continuous numerical model according to claim 1, characterized in that, When the non-continuous numerical calculation model is converted by using the topological closed intersection technology, a stratum model is also generated based on the ground surface, strata and selected calculation boundaries, the stratum model divides the initial block model into six material partitions based on profile stratum lines of stratum interfaces, and serves as a judgment basis for stratum parameter assignment.

Images