Method and system for determining parameters of stratified jointed rock mass model

Abstract

The invention relates to the technical field of layered jointed rock mass, and provides a parameter determination method and system based on a layered jointed rock mass model.The method comprises the steps that a cylinder with a measuring line as the axis is arranged in a geological measuring window, the total number of cracks of the cylinder is calculated, and the ratio of the total number of cracks to the height of the cylinder serves as the linear density of a dominant joint group; based on the linear density of the dominant joint group, the fracture density is calculated; and obtaining the dominant orientation and fracture equivalent diameter of the dominant joint group, combining the fracture density to form a tensor multiplication form, obtaining a second-order fracture tensor and a second-order damage tensor, and determining the shear strength index of the main level layer and the rock matrix. The practicability of the layered joint model can be improved, and the practical engineering value is achieved.

Description

Technical Field

[0001] The invention belongs to the technical field of layered jointed rock mass, and in particular relates to a parameter determination method and system based on a layered jointed rock mass model. Background Art

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

[0003] The layered rock mass commonly seen in engineering usually has special structural characteristics: the primary planes of oriented stratification such as bedding and foliation and the secondary structural planes such as the tangent joint group and random cracks in the layer cut each other to form a "secondary structure" layered jointed rock mass containing the primary and secondary structural planes. This "secondary structure" layered jointed rock mass often shows more complex anisotropic characteristics. Generally speaking, the primary plane belongs to the primary structural plane, which has a high degree of development, strong ductility and good connectivity, and is considered to be the main factor controlling the deformation and destruction of the surrounding rock of underground caverns. The secondary structural planes such as the dominant joint group belong to the secondary or tectonic structural plane, which has a transformative effect on the rock mass. Its mechanical effect is mainly reflected in two aspects: first, it has different occurrence characteristics, such as orientation, spacing and length, which destroys the integrity of the layered rock mass and the connectivity of the layers at a microscopic level, that is, the structural effect; second, it deteriorates the mechanical parameters of the layered rock mass at a macroscopic level, that is, the deterioration effect.

[0004] In this regard, the prior art has proposed a mechanical model that can reflect the structural effects and degradation effects of layered jointed rock masses. However, since the model contains three types of rock / body materials, namely, rock matrix, primary layer surface and secondary structural surface, and involves more than ten key mechanical parameters, the solution process of many parameters is relatively complicated.

[0005] The key rock mass mechanical parameters contained in the layered jointed rock mass are mainly divided into three categories. One category is the rock matrix related mechanical parameters: elastic modulus E r , Poisson's ratio ρ, internal friction angle and cohesion c r ; Another type of main level related mechanical parameters: normal stiffness KN J , tangential stiffness KS J , internal friction angle and cohesion c J ; The third category is the relevant mechanical parameters of the secondary structural surface: internal friction angle Cohesion JDFN , the second-order crack tensor F ij and the second-order damage tensor ω ij Among them, the second-order crack tensor F ij and the second-order damage tensor ω ijThere is no accurate method to measure and confirm this, which in turn affects the determination of effective shear strength indices for the rock matrix and primary levels, and is therefore an urgent problem to be solved in actual engineering tests. Summary of the invention

[0006] In order to solve the technical problems existing in the above-mentioned background technology, the present invention provides a parameter determination method and system based on a layered jointed rock model. Through the geological window measurement method, the equivalent mechanical parameters of the main level surface and the rock matrix under the influence of the secondary structural surface are derived from the two aspects of structural effect and degradation effect, and the key mechanical parameters required in the layered joint model are accurately and efficiently obtained, which helps to improve the practicability of the layered joint model and has practical engineering value.

[0007] In order to achieve the above object, the present invention adopts the following technical solution:

[0008] A first aspect of the present invention provides a parameter determination method based on a layered jointed rock mass model, comprising:

[0009] In the geological survey window, a cylinder with the survey line as the axis is set, and the total number of cracks in the cylinder is calculated. The total number of cracks and the ratio of the cylinder height are used as the linear density of the dominant joint group. Based on the linear density of the dominant joint group, the fracture density is calculated.

[0010] The dominant orientation of the dominant joint group and the equivalent diameter of the crack are obtained, combined with the crack density, and combined into a tensor multiplication form to obtain the second-order crack tensor and the second-order damage tensor, and then the shear strength index of the main level and rock matrix is ​​determined.

[0011] Furthermore, the second-order crack tensor is expressed as:

[0012]

[0013] Where ρ is the crack density, f(r) represents the distribution function of the crack equivalent diameter, r is the crack equivalent diameter, and r max represents the maximum equivalent diameter, n is the normal vector of the crack, and n i and n j is the crack normal cosine, E(n) represents the probability density function of the crack distribution orientation, and dΩ represents the microelement.

[0014] Furthermore, the steps for determining the shear strength index of the main level layer include: respectively obtaining the internal friction angle and cohesion of the secondary structural surface and the main level layer, calculating the damage tensor of the main level layer based on the second-order damage tensor, and based on the damage tensor, internal friction angle and cohesion of the main level layer, combining the internal friction angle and cohesion of the secondary structural surface, obtaining the shear strength index of the main level layer through the connectivity weighted average method.

[0015] Furthermore, the step of determining the shear strength index of the rock matrix includes: obtaining the internal friction angle and cohesion of the rock matrix, calculating the overall average damage level of the rock mass caused by the secondary structural surface based on the second-order damage tensor, combining the internal friction angle and cohesion of the rock matrix, and obtaining the effective shear strength index of the rock matrix through the connectivity weighted average method.

[0016] A second aspect of the present invention provides a parameter determination system based on a layered jointed rock mass model, comprising:

[0017] The fracture density determination module is configured to: in the geological survey window, set a cylinder with the survey line as the axis, calculate the total number of cracks in the cylinder, take the total number of cracks and the ratio of the cylinder height as the dominant joint group line density, and calculate the fracture density based on the dominant joint group line density;

[0018] The shear strength determination module is configured to obtain the dominant orientation of the dominant joint group and the equivalent diameter of the crack, combine them with the crack density, combine them into a tensor multiplication form, obtain the second-order crack tensor and the second-order damage tensor, and then determine the shear strength index of the main level and the rock matrix.

[0019] Furthermore, the second-order crack tensor is expressed as:

[0020]

[0021] Where ρ is the crack density, f(r) represents the distribution function of the crack equivalent diameter, r is the crack equivalent diameter, and r max represents the maximum equivalent diameter, n is the normal vector of the crack, and n i and n j is the crack normal cosine, E(n) represents the probability density function of the crack distribution orientation, and dΩ represents the microelement.

[0022] Furthermore, the steps for determining the shear strength index of the main level layer include: respectively obtaining the internal friction angle and cohesion of the secondary structural surface and the main level layer, calculating the damage tensor of the main level layer based on the second-order damage tensor, and based on the damage tensor, internal friction angle and cohesion of the main level layer, combining the internal friction angle and cohesion of the secondary structural surface, obtaining the shear strength index of the main level layer through the connectivity weighted average method.

[0023] Furthermore, the step of determining the shear strength index of the rock matrix includes: obtaining the internal friction angle and cohesion of the rock matrix, calculating the overall average damage level of the rock mass caused by the secondary structural surface based on the second-order damage tensor, combining the internal friction angle and cohesion of the rock matrix, and obtaining the effective shear strength index of the rock matrix through the connectivity weighted average method.

[0024] A third aspect of the present invention provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps in the parameter determination method based on a layered jointed rock mass model as described above.

[0025] The fourth aspect of the present invention provides a computer device, including a computer-readable storage medium, a processor, and a computer program stored on the computer-readable storage medium and executable on the processor, wherein when the processor executes the program, the steps in the parameter determination method based on a layered jointed rock model as described above are implemented.

[0026] Compared with the prior art, the present invention has the following beneficial effects:

[0027] The present invention uses the geological window measurement method to derive the equivalent mechanical parameters of the main level surface and the rock matrix under the influence of the secondary structural surface from the two aspects of structural effect and degradation effect, accurately and efficiently obtains the key mechanical parameters required in the layered joint model, helps to improve the practicability of the layered joint model, and has practical engineering value. BRIEF DESCRIPTION OF THE DRAWINGS

[0028] The accompanying drawings in the specification, which constitute a part of the present invention, are used to provide a further understanding of the present invention. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute improper limitations on the present invention.

[0029] Figure 1 It is a flow chart of a parameter determination method based on a layered jointed rock mass model according to the first embodiment of the present invention;

[0030] Figure 2 This is a schematic diagram of the distribution of the cylindrical inner structural surface of the geological measuring window according to the first embodiment of the present invention;

[0031] Figure 3 Schematic diagram of a layered jointed rock mass in a secondary structure according to Example 1 of the present invention;

[0032] Figure 4 It is a schematic diagram comparing the uniaxial compressive strength results of the numerical simulation test and the indoor test of muddy siltstone at different layer inclination angles according to the first embodiment of the present invention;

[0033] Figure 5 It is a schematic diagram comparing the elastic modulus results of the numerical simulation test and the indoor test of muddy siltstone at different layer inclination angles in Example 1 of the present invention;

[0034] Figure 6 It is a structural diagram of a computer device according to a fourth embodiment of the present invention. DETAILED DESCRIPTION

[0035] To make the purpose, technical solution and advantages of the embodiments of the present invention more clear, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the accompanying drawings in the embodiments of the present invention.

[0036] It should be noted that the following detailed descriptions are all illustrative and intended to provide further explanation of the present invention. Unless otherwise specified, all technical and scientific terms used herein have the same meanings as those commonly understood by those skilled in the art to which the present invention belongs.

[0037] Embodiment 1

[0038] This embodiment provides a parameter determination method based on a layered jointed rock mass model.

[0039] This embodiment provides a parameter determination method based on a layered jointed rock mass model to solve the problem that the prior art is difficult to measure and confirm the second-order fracture tensor F ij and the second-order damage tensor ω ij technical issues.

[0040] The parameter determination method based on a layered jointed rock mass model provided in this embodiment is used to determine the second-order fracture tensor F of the secondary structural surface of the layered jointed rock mass model. ij and the second-order damage tensor ω ij .

[0041] This embodiment provides a parameter determination method based on a layered jointed rock mass model, such as Figure 1 As shown, the following steps are included:

[0042] Step S11: The second-order crack tensor and the second-order damage tensor are expressed as crack scale, crack density and orientation vector and combined in the form of tensor multiplication.

[0043] Specifically, the crack tensor in step S11 combines the crack scale, density and orientation vector into a tensor multiplication form, as shown in formula (1), to characterize the geometric structural characteristics of the spatial distribution of the crack system. The crack tensor is an even-order dimensionless symmetric matrix, which realizes the single parameter characterization of the geometric structural characteristics of the crack space and provides a theoretical basis for the evaluation of the spatial deformation characteristics of the crack system. Assuming that the orientation of the secondary structural surface is independent of its effective diameter, the continuous form of the second-order crack tensor can be written as follows:

[0044]

[0045] Its discrete formula is:

[0046]

[0047] Where ρ is the volume density of the crack Defined as the number of fractures per unit volume of rock mass, m V is the number of cracks in a representative unit with volume V, f(r) represents the distribution function of the equivalent diameter of the crack, r is the equivalent diameter of the crack, r max represents the maximum equivalent diameter in V, n is the normal vector of the crack, n i and n j is the crack normal cosine, E(n) represents the probability density function of the crack distribution orientation, Ω represents a representative unit body, and dΩ represents a microelement.

[0048] In engineering, the secondary structural surface that affects the stability of the rock mass is mainly the dominant joint group. The crack density ρ of the dominant joint group in formula (1) is calculated as follows.

[0049] Step S12: In the geological survey window, a cylinder with the survey line T as the axis is set; along the survey line T, the linear density λ of the dominant joint group is calculated. T ; The crack density is obtained based on the line density.

[0050] Step S121: Figure 2 In the geological window shown, a normal line with direction n is set. j The measuring line T is a cylinder with a diameter of r and a height of L, and its upper and lower planes are composed of multiple structural surfaces with direction vectors of n and equivalent diameter of r. Its cross-sectional area is equal to The volume is Assuming L is long enough, the total number of structural faces in the cylinder can be written as

[0051] Step S122: For any microelement dΩ in the cylinder, calculate the number of cracks dN (k) , thus obtaining the total number of cracks N.

[0052] For any microelement dΩ, the number of cracks dN (k) It can be written as:

[0053]

[0054] The total number of cracks in the cylinder is:

[0055]

[0056] Step S123: along the measuring line T, calculate the dominant joint group linear density λ by the total number of cracks and the height of the cylinder T .

[0057] The linear density of the dominant joint group λ can be calculated along the survey line T T for:

[0058]

[0059] Where, d Τ Represents the gap spacing.

[0060] Step S124: Based on the linear density of the dominant joint group, the crack density ρ is obtained, that is, the structural surface volume density ρ:

[0061]

[0062] Step S13: Collect the dominant orientation (normal vector n) and effective diameter of the dominant joint group, and bring the crack density into the tensor multiplication form of the second-order crack tensor and the second-order damage tensor to obtain the second-order crack tensor F ij and the second-order damage tensor ω ij The value of .

[0063] Specifically, substitute formula (6) into formula (1) to obtain the second-order crack tensor of the secondary structural surface:

[0064]

[0065] In this embodiment, the dominant orientation of the dominant joint group and the effective diameter of the fracture in formula (7) are determined by on-site geological survey. It should be noted that the effective diameter here refers to the diameter of a disk equal to the actual fracture area.

[0066] In summary, in actual engineering, it is necessary to obtain the crack spacing d on site. Τ , the dominant orientation of the dominant joint group (normal vector n of the structural surface), and the effective diameter of the crack, the second-order crack tensor F can be calculated ij .

[0067] Joint second-order damage tensor ω ij Although the meanings of the second-order crack tensor Fij are different, their continuous and discrete expressions are completely consistent, as shown in formulas (1) and (2). In this embodiment, when the second-order crack tensor Fij is obtained, the second-order damage tensor ω can be confirmed at the same time. i j, that is:

[0068] The following is a specific calculation example using formula (7). Figure 3 As shown in the figure, the main level surface in the horizontal layered jointed rock mass is a highly continuous penetrating structural surface with a dip angle of α J , tends to be β J , the layer spacing is d T J ; A group of dominant joints, namely secondary structural planes, are developed in the rock mass, with an effective diameter of r, which follows a uniform distribution U(r J1 ,r J2 ), the inclination angle is α P, tends to be β P , the distance between the structural surfaces is d T P , r J1 and r J2 are the upper and lower limit parameters of the uniform distribution function, respectively. It should be noted that, for the sake of convenience, it is assumed that the main geometric parameters of the primary and secondary structural surfaces are uniformly distributed. In actual engineering, it is necessary to select a reasonable probability distribution function based on the specific statistical characteristics of each parameter. Since different probability density functions have no essential difference in the basic formulas for calculating the second-order crack tensor and the second-order damage tensor, it does not affect the rationality of this implementation case.

[0069] The second-order crack tensor and second-order damage of the secondary structural surface are calculated using formula (7):

[0070]

[0071] Step S15: Based on the second-order damage tensor ω obtained above ij It is used to confirm the effective internal friction angle and effective cohesion of the primary level of the layered jointed rock model, the effective internal friction angle and effective cohesion of the rock matrix, and then obtain the shear strength index of the secondary structure facing the primary level and the effective shear strength index of the rock matrix.

[0072] The yield criteria of the rock matrix and the main level in the layered jointed rock mass model all adopt the Mohr-Coulomb yield criterion with tensile cutoff, in which the tensile strength is positive and the compressive strength is negative. The secondary structural surfaces such as the dominant structural surface in the layered jointed rock mass are regarded as the damage of the rock mass and it is assumed that the reduction of the effective area is the main cause of the damage. This embodiment further establishes a second-order damage tensor describing the joint cracks based on the crack tensor in steps S11-S14 to reflect the direction of the discontinuity surface in the jointed rock mass and the change of its effective area.

[0073] The mechanical meaning of the second-order damage tensor of joints and fissures corresponds to the connectivity of the rock mass. The connectivity of joints in the rock mass is defined as the ratio of the area occupied by the crack part to the area occupied by the intact part on a certain section. The second-order damage tensor can be considered as a tensor expression of the connectivity. Therefore, the effective shear strength index of the rock mass matrix and the main level can be obtained by using the weighted average method of the connectivity.

[0074] The main plane is the inherent bedding plane of the layered jointed rock mass. It is assumed that all the main planes completely penetrate the representative unit with a volume of V, that is, the main planes in the numerical model are distributed throughout each calculation unit. The damage of the secondary structure to the strength of the layered rock mass is simplified to its damage to the main plane and the rock matrix.

[0075] (1) Damage tensor ω for the primary leveln It can be expressed as the normal unit vector n of the section (main level) and the second-order damage tensor ω ij A continuous function of :

[0076] ω n =n·ω ij ·n (8)

[0077] (2) Determine the internal friction angle of the secondary structural surface through indoor or field tests and cohesion c JDFN ; Internal friction angle of the main level and cohesion c J ; Internal friction angle of rock matrix and cohesion c r .

[0078] (3) The shear strength index of the secondary structure facing the primary level is obtained by using the weighted average method of connectivity:

[0079]

[0080] in, represents the effective internal friction angle of the primary layer, c eff J Represents the effective cohesion at the primary level; are the internal friction angle and cohesion of the secondary structural surface, and c J They are the internal friction angle and cohesion of the main level respectively. These strength parameters can be determined through indoor or field tests.

[0081] (4) Damage effect of secondary structure on rock matrix Simplified first invariant of second-order damage tensor I 1 ω The equivalent degradation of its strength parameters, where the first invariant I 1 ω It represents the overall average damage level of the rock mass caused by the secondary structural surface, and its expression is:

[0082] I 1 ω =ω 1 +ω 2 +ω 3 (10)

[0083] In the formula, ω 1 ,ω 2 ,ω 3 are the maximum damage amount, the intermediate damage amount and the minimum damage amount, respectively. In this embodiment, they are obtained by the following formula:

[0084] |ω ij -ωδij |=0 (11)

[0085] In the formula, ω represents the loss amount; || is the absolute value symbol; δ ij is the Kronecker function. It should be noted that formula (11) is a formula for finding the eigenvalue of a second-order matrix. After being solved, it is a cubic function with three solutions.

[0086] (5) The effective shear strength index of the rock matrix is ​​obtained by using the connectivity weighted average method:

[0087]

[0088] in, represents the effective internal friction angle of the rock matrix, c eff r Represents the effective cohesion of the rock matrix; and c r They are the internal friction angle and cohesion of the rock matrix, respectively, and can be determined through indoor or field tests.

[0089] The following is a detailed description based on the data:

[0090] 1. Methods for determining key mechanical parameters of rock matrix and main level.

[0091] First, prepare standard cylindrical specimens with different main-level surface inclination angles and main-level surface angles with the loading direction: 0°, 15°, 30°, 45°, 60°, 75° and 90°. (diameter × height); then, indoor uniaxial compression tests were carried out on samples with different inclination angles to obtain the key mechanical parameters of standard rock block samples containing rock matrix and main-level layers, including comprehensive elastic modulus, comprehensive internal friction angle, comprehensive cohesion and comprehensive tensile strength; finally, numerical simulation uniaxial compression tests corresponding to indoor uniaxial compression tests with different main-level inclination angles were carried out, and the key deformation and strength mechanical parameters of rock matrix and main-level layers were obtained by inversion, rock matrix: elastic modulus E r , Poisson's ratio ρ, internal friction angle Cohesion r ; Main level: Normal stiffness KN J , tangential stiffness KS J , internal friction angle Cohesion J .

[0092] Table 1 summarizes the average results of the main mechanical parameters of the uniaxial compression test of muddy siltstone at different layer inclination angles in a pumped storage power station project area. In general, the uniaxial test results have obvious transverse isotropy characteristics, which are specifically manifested in: its uniaxial compressive strength basically shows a "U"-shaped feature, and there is a minimum value when the loading direction and the layer are at an angle of about 30°. According to Jaeger's single weak surface theory (self-locking effect), the friction angle within the layer should satisfy the following formula:

[0093]

[0094] From this, we can infer that the internal friction angle of the layer

[0095] Table 1. Basic mechanical parameters of indoor rock blocks

[0096]

[0097] Further combined with the uniaxial compression test of saturated muddy siltstone cylindrical specimens at different layer inclinations, the mechanical parameters corresponding to the small-scale rock blocks were determined. The relevant comparison results are statistically analyzed in Figure 4 and Figure 5 The mechanical parameters of the rock matrix obtained by inversion are listed in Table 2.

[0098] Table 2. Range of mechanical parameters of small-scale argillaceous siltstone cylindrical specimens derived from inversion

[0099]

[0100] 2. Method for determining key mechanical parameters of secondary structural surfaces.

[0101] Conduct direct shear tests on secondary structural surfaces, select 3 to 5 normal stress levels, and obtain the internal friction angle of the structural surface. Cohesion JDFN .

[0102] The method for determining the mechanical parameters of the layered jointed rock mass model proposed in the present invention mainly includes the second-order crack tensor F of the secondary structural surface. ij and the second-order damage tensor ω ij , an indicator of the effective shear strength of the rock matrix and the main level.

[0103] Compared with the traditional constitutive model reflecting the anisotropy of layered rock mass, such as the ubiquitous joint model, the layered joint model can more effectively reflect the influence of the secondary structural surface on the mechanical properties of the layered jointed rock mass, but it is difficult to determine the second-order crack tensor F of the secondary structural surface. ij and the second-order damage tensor ω ijTo solve this problem, the present invention uses the geological window measurement method to derive the equivalent mechanical parameters of the main level surface and the rock matrix under the influence of the secondary structural surface from the two aspects of structural effect and degradation effect, accurately and efficiently obtains the key mechanical parameters required in the layered joint model, improves the practicability of the layered joint model, and has practical engineering value.

[0104] Embodiment 2

[0105] This embodiment provides a parameter determination system based on a layered jointed rock mass model, which specifically includes:

[0106] The fracture density determination module is configured to: in the geological survey window, set a cylinder with the survey line as the axis, calculate the total number of cracks in the cylinder, take the total number of cracks and the ratio of the cylinder height as the dominant joint group line density, and calculate the fracture density based on the dominant joint group line density;

[0107] The shear strength determination module is configured to obtain the dominant orientation of the dominant joint group and the equivalent diameter of the crack, combine them with the crack density, combine them into a tensor multiplication form, obtain the second-order crack tensor and the second-order damage tensor, and then determine the shear strength index of the main level and the rock matrix.

[0108] Among them, the second-order crack tensor is expressed as:

[0109]

[0110] Where ρ is the crack density, f(r) represents the distribution function of the crack equivalent diameter, r is the crack equivalent diameter, and r max represents the maximum equivalent diameter, n is the normal vector of the crack, and n i and n j is the crack normal cosine, E(n) represents the probability density function of the crack distribution orientation, and dΩ represents the microelement.

[0111] Among them, the steps for determining the shear strength index of the main level include: obtaining the internal friction angle and cohesion of the secondary structural surface and the main level respectively, calculating the damage tensor of the main level based on the second-order damage tensor, and based on the damage tensor, internal friction angle and cohesion of the main level, combined with the internal friction angle and cohesion of the secondary structural surface, the shear strength index of the main level is obtained by the weighted average method of connectivity rate.

[0112] Among them, the steps for determining the shear strength index of the rock matrix include: obtaining the internal friction angle and cohesion of the rock matrix, calculating the overall average damage level of the rock mass caused by the secondary structural surface based on the second-order damage tensor, and combining the internal friction angle and cohesion of the rock matrix to obtain the effective shear strength index of the rock matrix through the connectivity weighted average method.

[0113] It should be noted here that each module in this embodiment corresponds to each step in Example 1 one by one, and the specific implementation process is the same, which will not be repeated here.

[0114] Embodiment 3

[0115] This embodiment provides a computer-readable storage medium on which a computer program is stored. When the program is executed by a processor, the steps in the parameter determination method based on a layered jointed rock mass model as described in the first embodiment above are implemented.

[0116] Embodiment 4

[0117] This embodiment provides a computer device, such as Figure 6 As shown, it includes a display device, an input device, a computer-readable storage medium (volatile memory and non-volatile storage medium), a processor, a communication interface (i.e., a network interface), and a computer program stored on the computer-readable storage medium and executable on the processor, wherein the processor, the communication interface, and the computer-readable storage medium can be connected via a bus or other means. The communication interface is used to receive and send data, and when the processor executes the program, the steps in the parameter determination method based on a layered jointed rock model as described in the first embodiment above are implemented.

[0118] Among them, any reference to memory, storage, database or other media provided by the present application and used in the embodiments may include non-volatile and / or volatile memory. Non-volatile memory may include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM) or flash memory. Volatile memory may include random access memory (RAM) or external cache memory. As an illustration and not limitation, RAM is available in various forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), dual-speed data rate SDRAM (SSRSDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).

[0119] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0120] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0121] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0122] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. A parameter determination method based on a layered jointed rock mass model, characterized in that: include: In the geological survey window, a cylinder with the survey line as the axis is set, and the total number of cracks in the cylinder is calculated. The total number of cracks and the ratio of the cylinder height are used as the linear density of the dominant joint group. Based on the linear density of the dominant joint group, the fracture density is calculated. The dominant orientation of the dominant joint group and the equivalent diameter of the crack are obtained, combined with the crack density, and combined into a tensor multiplication form to obtain the second-order crack tensor and the second-order damage tensor, and then the shear strength index of the main level and rock matrix is ​​determined.

2. The parameter determination method based on a layered jointed rock mass model according to claim 1, characterized in that: The second-order crack tensor is expressed as: Where ρ is the crack density, f(r) represents the distribution function of the crack equivalent diameter, r is the crack equivalent diameter, and r max represents the maximum equivalent diameter, n is the normal vector of the crack, and n i and n j is the crack normal cosine, E(n) represents the probability density function of the crack distribution orientation, and dΩ represents the microelement.

3. The parameter determination method based on a layered jointed rock mass model according to claim 1, characterized in that: The steps for determining the shear strength index of the main level layer include: obtaining the internal friction angle and cohesion of the secondary structural surface and the main level layer respectively, calculating the damage tensor of the main level layer based on the second-order damage tensor, and obtaining the shear strength index of the main level layer by a connectivity weighted average method based on the damage tensor, internal friction angle and cohesion of the main level and the internal friction angle and cohesion of the secondary structural surface.

4. The parameter determination method based on a layered jointed rock mass model according to claim 1, characterized in that: The step of determining the shear strength index of the rock matrix includes: obtaining the internal friction angle and cohesion of the rock matrix, calculating the overall average damage level of the rock mass caused by the secondary structural surface based on the second-order damage tensor, and combining the internal friction angle and cohesion of the rock matrix to obtain the effective shear strength index of the rock matrix through a connectivity weighted average method.

5. A parameter determination system based on a layered jointed rock mass model, characterized in that: include: The fracture density determination module is configured to: in the geological survey window, set a cylinder with the survey line as the axis, calculate the total number of cracks in the cylinder, take the total number of cracks and the ratio of the cylinder height as the dominant joint group line density, and calculate the fracture density based on the dominant joint group line density; The shear strength determination module is configured to obtain the dominant orientation of the dominant joint group and the equivalent diameter of the crack, combine them with the crack density, combine them into a tensor multiplication form, obtain the second-order crack tensor and the second-order damage tensor, and then determine the shear strength index of the main level and the rock matrix.

6. The parameter determination system based on a layered jointed rock mass model according to claim 5, characterized in that: The second-order crack tensor is expressed as: Where ρ is the crack density, f(r) represents the distribution function of the crack equivalent diameter, r is the crack equivalent diameter, and r max represents the maximum equivalent diameter, n is the normal vector of the crack, and n i and n j is the crack normal cosine, E(n) represents the probability density function of the crack distribution orientation, and dΩ represents the microelement.

7. The parameter determination system based on a layered jointed rock mass model according to claim 5, characterized in that: The steps for determining the shear strength index of the main level layer include: obtaining the internal friction angle and cohesion of the secondary structural surface and the main level layer respectively, calculating the damage tensor of the main level layer based on the second-order damage tensor, and obtaining the shear strength index of the main level layer by a connectivity weighted average method based on the damage tensor, internal friction angle and cohesion of the main level and the internal friction angle and cohesion of the secondary structural surface.

8. The parameter determination system based on a layered jointed rock mass model according to claim 5, characterized in that: The step of determining the shear strength index of the rock matrix includes: obtaining the internal friction angle and cohesion of the rock matrix, calculating the overall average damage level of the rock mass caused by the secondary structural surface based on the second-order damage tensor, and combining the internal friction angle and cohesion of the rock matrix to obtain the effective shear strength index of the rock matrix through a connectivity weighted average method.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the program is executed by a processor, the steps in the parameter determination method based on a layered jointed rock model as described in any one of claims 1 to 4 are implemented.

10. A computer device comprising a computer-readable storage medium, a processor, and a computer program stored in the computer-readable storage medium and executable on the processor, characterized in that: When the processor executes the program, the steps in the parameter determination method based on a layered jointed rock model as described in any one of claims 1 to 4 are implemented.

Images