Complex non-uniform rock mass equivalent deformation modulus calculation method and system and medium

Abstract

The invention discloses an equivalent deformation modulus calculation method and system for a complex uneven rock mass and a medium. Relates to the technical field of rock engineering. The method comprises the following steps: partitioning according to geological information of rock masses and counting proportions, then respectively carrying out sound wave testing on the rock masses in different regions to obtain single-hole sound wave velocities, and calculating deformation moduli of the rock masses in each region by utilizing a correlation relationship between the single-hole wave velocities and the deformation moduli; a mixed equivalent model of strain coordination and stress coordination is constructed based on the rock mass strain energy proportion so as to calculate the equivalent deformation modulus of the special uneven rock mass; a new thought is provided for obtaining deformation modulus parameters containing multiple types of altered rocks, meanwhile, a targeted investigation and test scheme is provided, and an important basis is provided for design of rock mass slopes, caverns and the like.

Description

Technical Field

[0001] The present invention relates to the technical field of rock engineering, and particularly relates to a method, a system and a medium for calculating the equivalent deformation modulus of a complex and non-uniform rock mass. Background Art

[0002] Alteration refers to the geological process in which rocks or minerals undergo chemical or mineralogical composition changes under the action of fluids (such as hydrothermal fluids, groundwater, and meteoric water). Alteration is usually accompanied by changes in temperature, pressure, and fluid composition, resulting in the replacement of minerals in the original rock by new minerals, forming new mineral assemblages and structures. In recent years, with the construction development in China, a large number of projects have encountered the influence of altered rocks. For example, the clay-ification alteration causes slope instability when it meets water and softens, the altered rock zone has low strength leading to tunnel collapses, and the kaolin-ification alteration has low bearing capacity resulting in uneven settlement of building foundations, etc.

[0003] Due to the differences in fluid properties, original rock types, and tectonic environments during the alteration process, changes in mineral composition, structure, and cementation state will occur, thus forming many different types of altered rocks. The physical and mechanical properties of rock masses with different alteration types vary greatly, and they are often unevenly distributed, with various shapes and wrapped around each other, which will cause the unevenness of the rock mass properties and form geological defects.

[0004] Accurately obtaining the deformation modulus of such non-uniform rock masses containing multiple types of altered rocks is of great significance for the stability of engineering slopes, underground chambers, and foundations. Conventional indoor tests such as uniaxial compression tests are mostly applicable to relatively intact and uniform rock masses, but such rock masses are extremely non-uniform and it is difficult to prepare representative specimens. Moreover, large-scale in-situ deformation tests on site are costly and time-consuming. At the same time, there is no constitutive model for such special rock masses. Therefore, a new method is needed to efficiently and accurately obtain the deformation modulus of non-uniform rock masses containing multiple types of altered rocks. Summary of the Invention

[0005] The technical problem to be solved by the present invention is as follows: Inhomogeneous rock masses containing multiple types of altered rocks are difficult to be systematically divided into independent engineering geological units due to uneven distribution, various shapes, mutual wrapping of the altered rocks, and often significant differences in physical and mechanical properties. Limited by engineering geological testing techniques, methods, and sampling and sample preparation conditions, it is difficult to directly obtain relatively accurate deformation moduli through relevant specifications and tests. The purpose of the present invention is to provide a method, system, and medium for calculating the equivalent deformation modulus of complex inhomogeneous rock masses. The rock masses are divided into regions according to geological information and the proportion is statistically calculated. Then, the single-hole acoustic wave velocities of the rock masses in different regions are obtained through acoustic wave testing, and the deformation moduli of the rock masses in each region are calculated using the correlation between the single-hole wave velocity and the deformation modulus. A hybrid equivalent model of strain coordination and stress coordination is constructed based on the proportion of rock mass strain energy to calculate the equivalent deformation modulus of such special inhomogeneous rock masses, providing a new idea for obtaining the deformation modulus parameters of rock masses containing multiple types of altered rocks, and at the same time proposing a targeted investigation and test plan, providing an important basis for the design of slopes, caverns, etc. of such rock masses.

[0006] The present invention is realized through the following technical solutions:

[0007] This solution provides a method for calculating the equivalent deformation modulus of complex inhomogeneous rock masses, including:

[0008] Obtain the geological information of the rock masses in the target area; the geological information of the rock masses includes the types of altered rock masses and the content of each type of altered rock mass;

[0009] Based on the geological information of the rock masses and combined with the engineering geological properties, conduct engineering geological division of the types of altered rock masses in the target area to obtain the original rock area and n types of altered rock areas, and the statistical proportion of each type of altered rock area;

[0010] Conduct wave velocity testing on the rock masses in the original rock area and each altered rock area to obtain the single-hole acoustic wave velocity values of various rock masses, and calculate the deformation moduli of various rock masses according to the single-hole acoustic wave velocity values;

[0011] Input the deformation moduli of various rock masses into the constructed equivalent deformation modulus model to calculate the equivalent deformation modulus of the target area.

[0012] A further optimized solution is that the target area simultaneously meets the following conditions:

[0013] The target area simultaneously includes altered rock masses and original rock masses; among them, the altered rock masses include at least two types;

[0014] The distribution of the altered rock masses in the target area is uneven;

[0015] In the target area: there are differences in strength between the altered rock masses and the original rock masses, and there are differences in strength between different types of altered rock masses;

[0016] The proportion of altered rock masses in the target area is ≥ 30%.

[0017] A further optimization plan is to conduct engineering geological classification of altered rock mass types in the target area based on the geological information of rock masses and combined with engineering geological properties; the method includes: obtaining the content of each mineral in the geological information of rock masses in the target area.

[0018] The target area is divided into multiple sub - areas according to the mineral with the highest content, and the corresponding mineral is used as the representative mineral.

[0019] When the content of the representative mineral in each sub - area exceeds the content threshold, the representative mineral is used as a type of altered rock mass.

[0020] Statistical proportion of each type of altered rock mass.

[0021] A further optimization plan is to calculate the deformation modulus of various rock masses based on the single - hole acoustic wave velocity value; the method includes:

[0022] Calculate the deformation modulus of various rock masses according to the following formula:

[0023]

[0024]

[0025] Among them, E s represents the deformation modulus of the altered rock mass, V ps represents the single - hole wave velocity of the altered rock mass, E r represents the deformation modulus of the original rock mass, V pr represents the single - hole wave velocity of the original rock mass.

[0026] A further optimization plan is that the construction method of the equivalent deformation modulus model includes:

[0027] Taking the strain of all components in the material as the same when stressed, and the stress is distributed according to the elastic modulus of each region as a strain - uniform equivalent model; taking the stress of all components in the material as the same when stressed, and the strain is distributed according to the elastic modulus of each region as a stress - uniform equivalent model.

[0028] Obtain the first equivalent deformation modulus E t V of the target area based on the strain - uniform equivalent model, and obtain the second equivalent deformation modulus of the target area based on the stress - uniform equivalent model And introduce a weight coefficient α to obtain the equivalent deformation modulus model:

[0029]

[0030] A further optimization plan is that the first equivalent deformation modulus includes:

[0031]

[0032] The second equivalent deformation modulus includes:

[0033]

[0034] Wherein, P r represents the area ratio of the original rock on the research section; E r represents the deformation modulus of the original rock; P i represents the area ratio of the i-th type of altered rock mass on the research section; i = 1, 2, 3, …, n; E si represents the deformation modulus of the i-th type of altered rock mass.

[0035] A further optimized scheme is that the method for obtaining the weight coefficient α includes:

[0036] Obtain the total strain energy W V under the strain-uniform equivalent model, and the total strain energy W R under the stress-uniform equivalent model;

[0037] Calculate the weight coefficient α according to the following formula:

[0038] A further optimized scheme is that the total strain energy W V under the strain-uniform equivalent model is calculated according to the following formula:

[0039]

[0040] The total strain energy W R under the stress-uniform equivalent model is calculated according to the following formula:

[0041]

[0042] Wherein, P r represents the area proportion of the original rock on the research section; E r represents the deformation modulus of the original rock; P i represents the area proportion of the i-th type of altered rock mass on the research section; i = 1, 2, 3, …, n; E si represents the deformation modulus of the i-th type of altered rock mass; ε represents the strain of the rock mass; ρ represents the density of the rock mass; g represents the acceleration due to gravity; h represents the buried depth of the rock mass.

[0043] This solution also provides an equivalent deformation modulus calculation system for complex non-uniform rock masses, which is used to implement the equivalent deformation modulus calculation method for complex non-uniform rock masses described above; the system includes:

[0044] A collection module for obtaining the rock mass geological information of a target area; the rock mass geological information includes the types of altered rock masses and the content of each type of altered rock mass.

[0045] A division module for performing engineering geological division of the altered rock mass types on the target area based on the rock mass geological information and combining the engineering geological properties, to obtain the original rock area and n types of altered rock areas, and the proportion statistics of each type of altered rock area.

[0046] A first calculation module for performing wave velocity tests on the rock masses in the original rock area and each altered rock area to obtain the single-hole acoustic wave velocity values of various rock masses, and calculating the deformation modulus of various rock masses according to the single-hole acoustic wave velocity values.

[0047] A second calculation module for inputting the deformation modulus of various rock masses into the pre-constructed equivalent deformation modulus model to calculate the equivalent deformation modulus of the target area.

[0048] This solution also provides a computer-readable medium, on which a computer program is stored, and the computer program can be executed by a processor to implement the equivalent deformation modulus calculation method for complex non-uniform rock masses as described above.

[0049] Compared with the prior art, the present invention has the following advantages and beneficial effects:

[0050] The present invention provides an equivalent deformation modulus calculation method, system and medium for complex non-uniform rock masses; the area is divided according to the rock mass geological information and the proportion is statistically calculated, then the acoustic wave test is respectively carried out on the rock masses in different areas to obtain the single-hole acoustic wave velocity, the deformation modulus of the rock masses in each area is calculated by using the correlation between the single-hole wave velocity and the deformation modulus, and a hybrid equivalent model of strain coordination and stress coordination is constructed based on the proportion of rock mass strain energy to calculate the equivalent deformation modulus of such special non-uniform rock masses; it provides a new idea for obtaining the deformation modulus parameters of rock masses containing multiple types of altered rocks, and at the same time proposes a targeted investigation and test plan, providing an important basis for the design of such rock mass slopes, caverns, etc. BRIEF DESCRIPTION OF THE DRAWINGS

[0051] In order to more clearly illustrate the technical solutions of the exemplary embodiments of the present invention, the drawings required for use in the embodiments will be briefly introduced below. It should be understood that the following drawings only show some embodiments of the present invention, and therefore should not be regarded as limiting the scope. For those of ordinary skill in the art, other related drawings can be obtained based on these drawings without creative efforts. In the drawings:

[0052] Figure 1 It is a schematic flow chart of the equivalent deformation modulus calculation method for complex non-uniform rock masses;

[0053] Figure 2 It is a schematic diagram of the engineering geological division result of the altered rock mass types. Specific implementation manners

[0054] To make the objectives, technical solutions and advantages of the present invention clearer and more understandable, the present invention will be further described in detail below in conjunction with embodiments and the accompanying drawings. The illustrative implementation manners of the present invention and their descriptions are only used to explain the present invention and are not intended to limit the present invention.

[0055] For the heterogeneous rock mass containing multiple types of altered rocks, due to the uneven distribution of alteration, various shapes, mutual wrapping, and often large differences in physical and mechanical properties, it is difficult to divide these different types of altered rocks into independent engineering geological units. Limited by engineering geological testing techniques, methods, and sampling and sample preparation conditions, it is difficult to directly obtain relatively accurate deformation moduli through relevant specifications and tests. In view of this, the following embodiments are provided in this solution to solve the above technical problems.

[0056] Embodiment 1

[0057] This embodiment provides a method for calculating the equivalent deformation modulus of a complex heterogeneous rock mass, as Figure 1 shown, including:

[0058] Step 1: Obtain the rock mass geological information of the target area; the rock mass geological information includes the types of altered rock masses and the content of each type of altered rock mass;

[0059] The target area simultaneously meets the following conditions:

[0060] The target area simultaneously contains altered rock masses and unaltered original rock masses; among them, the altered rock masses include at least 2 types; such as chloritization alteration, argillization alteration, sericitization alteration, etc.;

[0061] The distribution of the altered rock masses in the target area is uneven; the uneven distribution described in this solution means various shapes and mutual wrapping, and it cannot be divided into independent engineering geological units for classification research;

[0062] In the target area: there are differences in strength between the altered rock masses and the original rock masses, and there are differences in strength between different types of altered rock masses;

[0063] The proportion of the altered rock masses in the target area ≥ 30%;

[0064] The rock mass geological information also includes the color, shape, material composition, genetic mechanism, development and distribution law, etc. of the altered rocks.

[0065] Step 2: Based on the rock mass geological information, combine the engineering geological properties to conduct engineering geological classification of the types of altered rock masses in the target area, obtain the original rock area and n types of altered rock areas, and the statistical proportion of each type of altered rock area;

[0066] This step specifically includes the method:

[0067] Obtain the contents of various minerals in the rock mass geology of the target area;

[0068] Divide the target area into multiple sub-areas according to the mineral with the highest content, and use the corresponding mineral as the representative mineral;

[0069] When the content of the representative mineral in each sub-area exceeds the content threshold, use the representative mineral as a type of altered rock mass;

[0070] Statistical the proportion of each type of altered rock mass.

[0071] Specifically, the following types of altered rock masses are divided:

[0072] The altered rock with sericite as the main constituent mineral, part of quartz and a small amount of feldspar, and the content of sericite mineral is greater than 50%, is divided into sericitization alteration;

[0073] The altered rock with chlorite as the main constituent mineral, part of carbonate minerals such as calcite and dolomite, and a small amount of epidote and quartz, and the content of chlorite mineral is greater than 50%, is divided into chloritization alteration;

[0074] The altered rock with clay minerals such as montmorillonite and kaolinite as the main constituent minerals, part of carbonate minerals such as calcite and dolomite, and a small amount of quartz, and the content of clay minerals is greater than 50%, is divided into argillization altered rock;

[0075] The altered rock with epidote as the main constituent mineral, part of quartz, calcite, feldspar and a small amount of other minerals, and the content of epidote mineral is greater than 50%, is divided into epidotization alteration.

[0076] The schematic diagram of the division result is as Figure 2 shown, obtaining Sr, S1, S2,..., S n area, where the Sr area is the original rock mass area, and S1 - S n are the corresponding areas of different types of altered rock masses, and n is the number of divided alteration types.

[0077] Step 3: Conduct wave velocity tests on the rock masses in the original rock area and each altered rock area to obtain the single-hole acoustic wave velocity values of various rock masses, and calculate the deformation modulus of various rock masses according to the single-hole acoustic wave velocity values; in this step, the methods for calculating the deformation modulus of various rock masses according to the single-hole acoustic wave velocity values include:

[0078] Calculate the deformation modulus of various rock masses according to the following formula:

[0079]

[0080] Among them, E s represents the deformation modulus of the altered rock mass, and V psDenote the single-hole wave velocity of the altered rock mass as E r Denote the deformation modulus of the original rock mass as V pr Denote the single-hole wave velocity of the original rock mass.

[0081] Step 4: Input the deformation modulus and other required parameters of various types of rock masses into the constructed equivalent deformation modulus model to calculate the equivalent deformation modulus of the target area.

[0082] Since the complex inhomogeneous rock mass contains a relatively intact and strongly bonded unaltered original rock part and an altered rock mass part that has been altered and weakly bonded, the former cannot undergo relative slip due to strong interface bonding and has strain coordination when stressed, while the latter has some weak interfaces and has stress coordination when stressed. Based on the above analysis, for such special rock masses, it is impossible to accurately express them through a single equivalent model. Therefore, a hybrid model based on the strain uniformity assumption and the stress uniformity assumption needs to be considered simultaneously to improve the accuracy.

[0083] The construction method of the equivalent deformation modulus model includes:

[0084] Taking the strain of all components in the material to be the same when stressed, and the stress being distributed according to the elastic modulus of each region as the strain-uniform equivalent model; taking the stress of all components in the material to be the same when stressed, and the strain being distributed according to the elastic modulus of each region as the stress-uniform equivalent model; specifically, the strain-uniform equivalent model includes:

[0085]

[0086] The stress-uniform equivalent model includes:

[0087]

[0088] Among them, P r Denote the area ratio of the original rock on the research cross-section; E r Denote the deformation modulus of the original rock; P i Denote the area ratio of the i-th type of altered rock mass on the research cross-section; i = 1, 2, 3,..., n; E si Denote the deformation modulus of the i-th type of altered rock mass.

[0089] Obtain the first equivalent deformation modulus of the target area based on the strain-uniform equivalent model Obtain the second equivalent deformation modulus of the target area based on the stress-uniform equivalent model In order to comprehensively reflect the actual characteristics of the inhomogeneous structure of the rock mass, the weighted average method is adopted, and the weight coefficient α is introduced to reconcile the two extreme models to make the calculation of the equivalent modulus more accurate and closer to the actual situation. The equivalent deformation modulus model is obtained by introducing the weight coefficient α:

[0090]

[0091] The method for obtaining the weight coefficient α includes: obtaining the total strain energy W under the strain-uniform equivalent model V , and the total strain energy W under the stress-uniform equivalent model R ; calculating the weight coefficient α according to the following formula:

[0092] Strain energy is the elastic energy stored in materials, reflecting the contribution of different components to the overall deformation during the stress process. For non-uniform rock masses, the strain energy stored in regions with different stiffnesses is different. Within the linear elastic range, the calculation method for the strain energy of a rock mass is as follows:

[0093]

[0094] Under the strain-uniform equivalent model, the strains of each component in the rock mass are the same, but the stresses are different. The part with a higher deformation modulus will store more strain energy. At this time, the total strain energy W V is calculated according to the following formula:

[0095]

[0096] Under the stress-uniform equivalent model, the stresses of each component in the rock mass are the same, but the strains are different. The part with a lower deformation modulus will store more strain energy. At this time, the total strain energy W R is calculated according to the following formula:

[0097]

[0098] where P r represents the area ratio of the original rock on the research section; E r represents the deformation modulus of the original rock; P i represents the area ratio of the i-th type of altered rock mass on the research section; i = 1, 2, 3,..., n; E si represents the deformation modulus of the i-th type of altered rock mass; ε represents the strain of the rock mass; ρ represents the density of the rock mass; g represents the acceleration due to gravity; h represents the buried depth of the rock mass.

[0099] To a certain extent, the strain energy ratio of the two models quantifies the tendency of the rock mass structure and directly reflects the mechanical contributions of different components.

[0100] Example 2

[0101] This example provides an equivalent deformation modulus calculation system for complex non-uniform rock masses, which is used to implement the equivalent deformation modulus calculation method for complex non-uniform rock masses described in Example 1; the system includes:

[0102] A collection module, which is used to obtain the geological information of the rock mass in the target area; the geological information of the rock mass includes the types of altered rock masses and the content of each type of altered rock mass;

[0103] A dividing module, configured to perform engineering geological division of altered rock mass types on a target area based on rock mass geological information and combined with engineering geological properties, obtain a protolith area and n types of altered rock mass areas, and perform proportion statistics of each type of altered rock mass area;

[0104] A first calculation module, configured to perform wave velocity tests on the rock masses in the protolith area and each altered rock mass area, obtain the single-hole acoustic wave velocity values of various rock masses, and calculate the deformation moduli of various rock masses according to the single-hole acoustic wave velocity values;

[0105] A second calculation module, configured to input the deformation moduli of various rock masses into a constructed equivalent deformation modulus model to calculate the equivalent deformation modulus of the target area.

[0106] Embodiment 3

[0107] This embodiment provides a computer-readable medium, on which a computer program is stored, and the computer program, when executed by a processor, can implement the method for calculating the equivalent deformation modulus of a complex inhomogeneous rock mass as described in Embodiment 1.

[0108] Step 1: Obtain the rock mass geological information of the target area; the rock mass geological information includes the types of altered rock masses and the content of each type of altered rock mass;

[0109] Step 2: Based on the rock mass geological information, perform engineering geological division of altered rock mass types on the target area in combination with engineering geological properties to obtain a protolith area and n types of altered rock mass areas;

[0110] Step 3: Perform wave velocity tests on the rock masses in the protolith area and each altered rock mass area, obtain the single-hole acoustic wave velocity values of various rock masses, and calculate the deformation moduli of various rock masses according to the single-hole acoustic wave velocity values;

[0111] Step 4: Input the deformation moduli of various rock masses into a constructed equivalent deformation modulus model to calculate the equivalent deformation modulus of the target area.

[0112] Embodiment 4

[0113] This embodiment takes an inhomogeneous rock mass containing altered rock in a hydropower station project in a southwestern canyon area as an example for detailed description.

[0114] S1. Conduct geological investigation and testing work on the rock masses in the study area. The geological investigation and testing show that the rock masses contain altered rock and unaltered protolith. The types of alteration include epidotization alteration, argillization alteration, and chloritization alteration. Moreover, the three types of altered rock are unevenly distributed, with different shapes, wrapped around each other, and significant differences in strength. The proportion of altered rock mass is 45%, so such rock masses meet the requirements of the present invention. The rock mass density is 2700 kg / m 3 , and the rock mass burial depth is 150 m.

[0115] S2. According to the geological survey and test results, combined with the engineering geological properties, the regional division of 3 different types of altered rock masses is carried out, namely S1 (epidotization altered area), S2 (argillization altered area), and S3 (chloritization altered area).

[0116] S3. According to the divided altered types and regions, the area ratios P of the unaltered protolith and the 3 types of altered rocks on the research section are respectively counted r = 0.55, P1 = 0.21, P2 = 0.15, P3 = 0.09.

[0117] S4. According to the geological division results of different altered types, the wave velocity tests are respectively carried out on the protolith and the altered rock masses within the ranges of each altered type area to obtain the single-hole acoustic wave velocity values V of the protolith and different types of altered rock masses in each divided area pr = 5120 m / s, V ps1 = 4236 m / s, V ps2 = 3551 m / s, V ps3 = 3896 m / s. The deformation moduli of the protolith and each type of altered rock are calculated to be E r = 19.5 GPa, E s1 = 5.6 GPa, E s2 = 2.6 GPa, E s3 = 3.8 GPa.

[0118] S5. Calculate the deformation moduli under the strain uniform equivalent model and the stress uniform equivalent model respectively, and the results are E t V = 12.6 GPa, E t R = 6.8 GPa.

[0119] S6. Determination of the weight coefficient α. Based on the rock mass density and sampling buried depth parameters, the stress σ = 4 MPa is calculated, and the total strain energy W V = 636.6 J / m 3 under the condition of the strain uniform equivalent model, and the total strain energy W R = 1176.7 J / m 3 under the condition of the stress uniform equivalent model, then the weight coefficient α = 0.35.

[0120] S7. Calculation of the equivalent deformation modulus. Using the deformation moduli under the strain uniform equivalent model and the stress uniform equivalent model and the weight coefficient, the equivalent deformation modulus E of this type of inhomogeneous rock mass is obtained t= 8.8 GPa. The deformation modulus of such rock masses was measured by large-scale in-situ deformation tests to be 9.8 GPa. It can be seen from this that the calculation results of the equivalent model using the strain uniformity assumption are on the high side, the calculation results of the equivalent model using the stress uniformity assumption are on the low side, and the results calculated by the hybrid model constructed based on the proportion of rock mass strain energy are more in line with the actual situation.

[0121] In order to obtain the equivalent deformation modulus of complex inhomogeneous rock masses relatively accurately and efficiently, based on the characteristics of such special rock masses, according to geological investigation and testing, they are divided into unaltered rock masses and various types of altered rock masses (such as chloritization alteration, argillization alteration, sericitization alteration, etc.) and the proportion is counted. Then, sonic tests are carried out on the rock masses in different regions to obtain the single-hole sonic wave velocity, and the deformation modulus of the rock masses in each region is calculated by using the correlation between the single-hole wave velocity and the deformation modulus. Finally, a hybrid equivalent model of strain coordination and stress coordination is constructed based on the proportion of rock mass strain energy to calculate the equivalent deformation modulus of such special inhomogeneous rock masses; it provides a new idea for obtaining the deformation modulus parameters of rock masses with multiple types of inhomogeneous alterations, and at the same time proposes a targeted investigation and test plan, providing an important basis for the design of slopes, caverns, etc. of such rock masses.

[0122] The specific embodiments described above have further elaborated on the purpose, technical solutions, and beneficial effects of the present invention. It should be understood that the above description is only the specific embodiments of the present invention and is not used to limit the protection scope of the present invention. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present invention shall be included in the protection scope of the present invention.

Claims

1. A method for calculating the equivalent deformation modulus of a complex inhomogeneous rock mass, characterized in that, including: Obtain the rock mass geological information of the target area; the rock mass geological information includes the altered rock mass types and the content of each altered rock mass type; Based on the rock mass geological information, combined with the engineering geological properties, conduct an engineering geological division of the altered rock mass types in the target area to obtain the original rock area and n types of altered rock areas, and the proportion statistics of each type of altered rock area; Conduct wave velocity tests on the rock masses in the original rock area and each altered rock area to obtain the single-hole acoustic wave velocity values of various rock masses, and calculate the deformation moduli of various rock masses according to the single-hole acoustic wave velocity values; Input the deformation moduli of various rock masses into the constructed equivalent deformation modulus model to calculate the equivalent deformation modulus of the target area.

2. The method for calculating the equivalent deformation modulus of a complex inhomogeneous rock mass according to claim 1, characterized in that, The target area simultaneously satisfies the following conditions: The target area simultaneously contains altered rock masses and original rock masses; among which, the altered rock masses include at least two types; The distribution of the altered rock masses in the target area is uneven; In the target area: there are differences in strength between the altered rock masses and the original rock masses, and there are differences in strength between different types of altered rock masses; The proportion of the altered rock masses in the target area ≥ 30%; 3. The method for calculating the equivalent deformation modulus of a complex inhomogeneous rock mass according to claim 2, characterized in that, The engineering geological division of the altered rock mass types in the target area based on the rock mass geological information, combined with the engineering geological properties; includes Method: Obtain the content of each mineral in the rock mass geology of the target area; Divide the target area into multiple sub-areas according to the mineral with the highest content, and use the corresponding mineral as the representative mineral; When the content of the representative mineral in each sub-area exceeds the content threshold, use the representative mineral as one type of altered rock mass; Statistical proportion of each type of altered rock mass.

4. The method for calculating the equivalent deformation modulus of a complex inhomogeneous rock mass according to claim 2, wherein The calculation of the deformation moduli of various rock masses according to the single-hole acoustic wave velocity values; includes Method: Calculate the deformation moduli of various rock masses according to the following formula: Among them, E s represents the deformation modulus of the altered rock mass, V ps represents the single-hole wave velocity of the altered rock mass, E r represents the deformation modulus of the original rock mass, V pr represents the single-hole wave velocity of the original rock mass.

5. The method for calculating the equivalent deformation modulus of a complex inhomogeneous rock mass according to claim 2, characterized in that, The construction method of the equivalent deformation modulus model includes: Taking the model with the same strain of all components in the material when stressed and the stress distributed according to the elastic modulus of each area as the strain-uniform equivalent model; taking the model with the same stress of all components in the material when stressed and the strain distributed according to the elastic modulus of each area as the stress-uniform equivalent model; Obtain the first equivalent deformation modulus of the target area based on the strain uniform equivalent model Obtain the second equivalent deformation modulus of the target area based on the stress uniform equivalent model And introduce a weight coefficient α to obtain an equivalent deformation modulus model:

6. The method for calculating the equivalent deformation modulus of a complex inhomogeneous rock mass according to claim 5, characterized in that The first equivalent deformation modulus includes: Among them, P r represents the area ratio of the original rock on the research section; E r represents the deformation modulus of the original rock; P i represents the area ratio of the i-th type of altered rock mass on the research section; i = 1, 2, 3, …, n; E si represents the deformation modulus of the i-th type of altered rock mass; The second equivalent deformation modulus includes: Among them, P r represents the area ratio of the original rock on the research section; E r represents the deformation modulus of the original rock; P i represents the area ratio of the i-th type of altered rock mass on the research section; i = 1, 2, 3, …, n; E si represents the deformation modulus of the i-th type of altered rock mass.

7. The method for calculating the equivalent deformation modulus of a complex inhomogeneous rock mass according to claim 5, characterized in that, The method for obtaining the weight coefficient α includes: Obtain the total strain energy \(W\) under the strain-uniform equivalent model V and the total strain energy \(W\) under the stress-uniform equivalent model R ; Calculate the weight coefficient α according to the following formula:

8. The method for calculating the equivalent deformation modulus of a complex inhomogeneous rock mass according to claim 7, characterized in that, The total strain energy W under the strain-uniform equivalent model V Calculated according to the following formula: The total strain energy W under the stress uniform equivalent model R It is calculated according to the following formula: Among them, P r represents the area ratio of the original rock on the research section; E r represents the deformation modulus of the original rock; P i represents the area ratio of the i-th type of altered rock mass on the research section; i = 1, 2, 3, …, n; E si represents the deformation modulus of the i-th type of altered rock mass; ε represents the strain of the rock mass; ρ represents the density of the rock mass; g represents the acceleration due to gravity; h represents the depth of burial of the rock mass.

9. The equivalent deformation modulus calculation system for complex inhomogeneous rock masses is characterized in that For implementing the equivalent deformation modulus calculation method of the complex and uneven rock mass described in any one of claims 1-8; the system includes: A collection module, used to obtain the rock mass geological information of the target area; the rock mass geological information includes the altered rock mass types and the content of each altered rock mass type; A division module, used to conduct an engineering geological division of the altered rock mass types in the target area based on the rock mass geological information, combined with the engineering geological properties, to obtain the original rock area and n types of altered rock areas, and the proportion statistics of each type of altered rock area; A first calculation module, used to conduct wave velocity tests on the rock masses in the original rock area and each altered rock area to obtain the single-hole acoustic wave velocity values of various rock masses, and calculate the deformation moduli of various rock masses according to the single-hole acoustic wave velocity values; A second calculation module, used to input the deformation moduli of various rock masses into the constructed equivalent deformation modulus model to calculate the equivalent deformation modulus of the target area.

10. A computer-readable medium having a computer program stored thereon, characterized in that, The computer program, when executed by a processor, can implement the method for calculating the equivalent deformation modulus of a complex inhomogeneous rock mass as described in any one of claims 1-8.

Images