A method and apparatus for correcting the cubic law of seepage through rough single-jointed channels in granite.

Abstract

This invention provides a method and apparatus for correcting the cubic law of seepage through rough single-jointed channels in granite. The method includes: dividing the target rock mass into several measurement intervals and determining the roughness coefficient of the actual joints in the target rock mass; determining the actual permeability of a single fracture in the target rock mass based on the cubic law correction, and correcting the initial width of the fracture to obtain the final width after fracture displacement; establishing a calculation model of the single-fracture rock mass under triaxial stress conditions, and obtaining the displacement data of the target rock mass through the calculation model. This invention makes the simulation results more accurate, optimizes the permeability calculation formula, and solves the problem of poor accuracy in existing studies of seepage characteristics in fractured rock masses.

Description

Technical Field

[0001] This invention relates to the field of rock mass seepage characteristics technology, and in particular to a method and apparatus for correcting the cubic law of seepage through rough single-jointed channels in granite. Background Technology

[0002] Many geological and engineering hazards related to seepage in fractured rock masses occur frequently, making seepage in fractured rock masses a major research topic. Currently, scholars both domestically and internationally have conducted extensive research on the seepage characteristics of single-jointed channels and have achieved a series of practically significant research results.

[0003] Regarding the seepage characteristics of fractured rock masses under stress-seepage coupling, some scholars have established a multi-fracture network model for coupled analysis of the seepage field and stress field of the rock mass by classifying the fractures in the rock mass. Calculation results from actual engineering examples demonstrate the necessity of incorporating the coupling effect into engineering design. In terms of stress-strain-seepage coupling models, some scholars have conducted theoretical analysis and laboratory experiments on the permeability characteristics of single fractures in granite under triaxial stress, establishing a coupling model of fracture permeability and triaxial stress. The fracture permeability coefficient exhibits an exponential relationship with triaxial stress, and it is proposed that the lateral stress influence coefficient can be used to evaluate the impact of lateral stress on the permeability coefficient. Regarding shear-seepage coupling, some scholars have conducted seepage tests on split granite samples, studying the influence of normal stress on the nonlinear seepage characteristics of single fractures under shear dislocation conditions. They found that as the normal stress increases, the effective contact area of ​​the fracture surface increases, making the seepage path more complex and enhancing the nonlinear seepage capacity of the fracture.

[0004] Although many scholars at home and abroad have conducted extensive research on the seepage characteristics of single-jointed channels, the seepage test results still show a significant gap between the actual seepage situation and the real situation. There is a lack of in-depth research on the seepage characteristics of single-jointed fractures in rough rocks, and the accuracy of the tests is not good.

[0005] There is currently no effective solution to the problem of poor accuracy in existing studies on the seepage characteristics of fractured rock masses. Summary of the Invention

[0006] This invention provides a method and apparatus for correcting the cubic law of seepage in coarse single-jointed channels of granite, in order to solve the defect of poor accuracy in the study of seepage characteristics of fractured rock masses in the prior art.

[0007] In a first aspect, the present invention provides a method for correcting the cubic law of seepage through rough single-jointed channels in granite, comprising:

[0008] The target rock mass is divided into several measurement intervals, and the roughness coefficient of the actual joints of the target rock mass is determined;

[0009] The permeability of a single fracture in the target rock mass is determined by correcting the cubic law, and the initial width of the fracture in the target rock mass is corrected to obtain the fracture width after the fracture displacement.

[0010] A calculation model of a single-fracture rock mass under triaxial stress conditions is established, and the displacement data of the target rock mass is obtained through the calculation model of the single-fracture rock mass.

[0011] Based on the stress state of the target rock mass and the displacement data of the target rock mass, the fracture width of the target rock mass is determined.

[0012] By combining the fracture width of the target rock mass and the actual permeability of a single fracture in the target rock mass, the roughness coefficient of the actual joints of the target rock mass is introduced, and the relationship between the roughness coefficient and the permeability is corrected and fitted to obtain an initial empirical relationship.

[0013] The confining pressure of the target rock mass is compensated and corrected, and the corrected confining pressure is then substituted into the initial empirical relation to obtain the target empirical relation.

[0014] According to the present invention, a method for correcting the cubic law of seepage through rough single joints in granite is provided to determine the roughness coefficient of the actual joints in the target rock mass, including:

[0015] Intermediate parameters are obtained based on the number of segments in the measurement interval of the target rock mass and the protrusion height of each measuring point;

[0016] The intermediate parameters are substituted into a pre-set regression equation to determine the roughness coefficient of the actual joints of the target rock mass.

[0017] According to the present invention, a method for correcting the cubic law of seepage through rough single-jointed channels in granite is provided. Based on the correction of the cubic law, the actual permeability of a single fracture in the target rock mass is determined, including:

[0018] Based on the cubic law, the relationship between the amount of seepage through the fracture surface and the fracture width of the target rock mass under steady water flow is determined.

[0019] The actual permeability of a single fracture in the target rock mass is determined by correcting the relationship between the permeability of the target rock mass through the fracture surface and the fracture width under stable water flow.

[0020] According to the present invention, a method for correcting the cubic law of seepage through rough single-jointed channels in granite is provided, which corrects the initial width of the fractures in the target rock mass to obtain the fracture width after the fracture displacement in the target rock mass, including:

[0021] The fracture deformation of the target rock mass is determined based on the hyperbolic deformation formula of the normal deformation stress change of the fracture structure surface.

[0022] The initial width of the cracks in the target rock mass is corrected based on the crack deformation of the target rock mass to obtain the crack width after the crack displacement.

[0023] According to the present invention, a method for correcting the cubic law of seepage through rough single-jointed channels in granite is provided, which establishes a calculation model of a single-fracture rock mass under triaxial stress conditions, and obtains the displacement data of the target rock mass through the single-fracture rock mass calculation model, including:

[0024] Construct the calculation model of the single-fracture rock mass;

[0025] The rock blocks on both sides of the fracture surface of the target rock mass are regarded as linear elastic bodies. The total displacement, fracture displacement and displacement of the rock blocks on both sides of the single fracture rock mass in the target direction are obtained by the single fracture rock mass calculation model.

[0026] According to the present invention, a method for correcting the cubic law of seepage through rough single-jointed channels in granite is provided, which determines the fracture width of the target rock mass based on the stress state of the target rock mass and the displacement data of the target rock mass, including:

[0027] Based on the stress state of the target rock mass and the displacement data of the target rock mass, the strain increment of the single-fracture rock mass and the rock blocks on both sides of the target rock mass along the target direction is obtained.

[0028] The fracture width of the target rock mass is determined based on the strain increment of the single fracture rock mass and the rock blocks on both sides along the target direction.

[0029] According to the present invention, a method for correcting the cubic law of seepage through rough single joint channels in granite is provided. This method combines the fracture width of the target rock mass and the actual permeability of a single fracture in the target rock mass, introduces the roughness coefficient of the actual joints in the target rock mass, and corrects and fits the relationship between the roughness coefficient and the permeability to obtain an initial empirical formula, including:

[0030] The quantitative relationship between the flow rate of a single joint in the target rock mass and its roughness, confining pressure, and seepage pressure is obtained. The roughness coefficient of the actual joint in the target rock mass is introduced to correct the relationship between the roughness coefficient and the seepage rate.

[0031] The image of the permeability and roughness coefficient of the target rock mass is obtained, and the relationship between the permeability and roughness coefficient of a single joint fracture is fitted to obtain the fitting coefficient.

[0032] The initial empirical relation is determined based on the fitting coefficients.

[0033] Secondly, the present invention also provides a device for correcting the cubic law of seepage in rough single-jointed channels of granite, comprising:

[0034] The coefficient determination module is used to divide the target rock mass into several measurement intervals and determine the roughness coefficient of the actual joints of the target rock mass;

[0035] The width correction module is used to correct and determine the actual permeability of a single fracture in the target rock mass based on the cubic law, and to correct the initial width of the fracture in the target rock mass to obtain the fracture width after the fracture displacement in the target rock mass.

[0036] The data processing module is used to establish a calculation model of a single-fracture rock mass under triaxial stress conditions, and to obtain the displacement data of the target rock mass through the single-fracture rock mass calculation model.

[0037] The width determination module is used to determine the crack width of the target rock mass based on the stress state of the target rock mass and the displacement data of the target rock mass.

[0038] The correction fitting module is used to combine the fracture width of the target rock mass and the actual permeability of a single fracture in the target rock mass, introduce the roughness coefficient of the actual joint of the target rock mass, and correct and fit the relationship between the roughness coefficient and the permeability to obtain an initial empirical relationship.

[0039] The confining pressure correction module is used to compensate and correct the confining pressure of the target rock mass, and then input the corrected confining pressure into the initial empirical relation to obtain the target empirical relation.

[0040] Thirdly, the present invention also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the cubic law correction method for seepage through rough single jointed channels of granite as described in the first aspect above.

[0041] In a fourth aspect, the present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the method for correcting the cubic law of seepage through rough single-jointed channels in granite as described in the first aspect above.

[0042] Fifthly, the present invention also provides a computer program product, including a computer program that, when executed by a processor, implements the cubic law correction method for seepage through rough single-jointed channels in granite as described in the first aspect above.

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

[0044] This invention provides a method for correcting the cubic law of seepage in granite rough single-joint channels. In this method, the initial fracture width of the fractured rock mass is corrected based on the displacement of rock blocks on both sides under different stress states. A relationship between the permeability of rough joints and the roughness of the single joint, the seepage pressure, and the confining pressure is proposed. Based on the proposed permeability calculation formula, different stress states are substituted into the formula for back-calculation. It is found that the formula shows better agreement under different roughness coefficients and different seepage pressures. To improve the applicability of this method under different confining pressures, a second correction is performed for different confining pressures, resulting in a new permeability formula. Through the above process, the simulation results are more accurate, a new empirical relationship between the permeability of rough joints and the roughness of the single joint, the seepage pressure, and the confining pressure is proposed, and the permeability calculation formula is optimized, solving the problem of poor accuracy in existing studies of seepage characteristics in fractured rock masses. Attached Figure Description

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

[0046] Figure 1 This is a flowchart of the method for correcting the cubic law of seepage through rough single-jointed channels in granite provided by the present invention;

[0047] Figure 2 This is a calculation model diagram of a single-fracture rock mass in an embodiment of the present invention;

[0048] Figure 3 This is a schematic diagram of the relationship between the permeation amount Q and JRC obtained in an embodiment of the present invention;

[0049] Figure 4 This is a comparative analysis chart of the experimental and corrected values ​​of permeation in an embodiment of the present invention;

[0050] Figure 5 This is a structural block diagram of the granite rough single-jointed channel seepage cubic law correction device provided by the present invention;

[0051] Figure 6 This is a schematic diagram of the structure of the electronic device provided by the present invention. Detailed Implementation

[0052] To make the objectives, technical solutions, and advantages of this invention clearer, the technical solutions of this invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of this invention. All other embodiments obtained by those skilled in the art based on the embodiments of this invention without creative effort are within the scope of protection of this invention.

[0053] This invention provides a method for correcting the cubic law of seepage through rough single-jointed channels in granite. Figure 1 This is a flowchart of the method for correcting the cubic law of seepage in coarse single-jointed channels of granite provided by the present invention, as shown below. Figure 1 As shown, the method includes the following steps:

[0054] Step S101: Divide the target rock mass into several measurement intervals and determine the roughness coefficient of the actual joints of the target rock mass;

[0055] Step S102: Based on the cubic law, the actual permeability of a single fracture in the target rock mass is determined, and the initial width of the fracture in the target rock mass is corrected to obtain the fracture width after the fracture displacement in the target rock mass.

[0056] Step S103: Establish a calculation model of a single-fracture rock mass under triaxial stress conditions, and obtain the displacement data of the target rock mass through the single-fracture rock mass calculation model;

[0057] Step S104: Determine the fracture width of the target rock mass based on the stress state of the target rock mass and the displacement data of the target rock mass.

[0058] Step S105: Combining the fracture width of the target rock mass and the actual permeability of a single fracture in the target rock mass, the roughness coefficient of the actual joints of the target rock mass is introduced, and the relationship between the roughness coefficient and the permeability is corrected and fitted to obtain the initial empirical relationship.

[0059] Step S106: Compensate and correct the confining pressure of the target rock mass, and substitute the corrected confining pressure into the initial empirical relation to obtain the target empirical relation.

[0060] For example, firstly, the roughness coefficient of the target rock mass specimen with a length of m cm is calculated, and then the roughness coefficient of the actual joint (JRC, Joint Roughness Coefficient) is obtained through an empirical formula. Next, the permeability of the actual single fracture is calculated based on a correction to the cubic law, and the initial fracture width of the target rock mass is corrected. Finally, the fracture width after displacement is obtained. Then, a calculation model of the single-fracture rock mass under triaxial stress conditions is established, and the permeability along the target direction is obtained through this model. Displacement data of a single-fracture rock mass (direction). The displacement data includes total displacement. Crack displacement and displacement of rock blocks on both sides Next, based on the stress state of the fractured rock mass, the strain increment of the single fractured rock mass and the rock blocks on both sides of the fracture along the Z direction is calculated using the calculation formulas for the fracture width. Then, the fracture width is calculated by inverse calculation using the formula for calculating the fracture width. Next, the quantitative relationship between the flow rate of a single joint fracture and roughness, confining pressure, and permeability pressure is obtained. The expression for the fracture width is substituted into the expression for the permeability, and the JRC value is introduced. Based on the permeability test results of different JRC values ​​under the same stress state, the relationship is corrected and fitted. Finally, the fitted coefficients are substituted into the expression to obtain the initial empirical relationship between the average permeability of the fracture and the roughness coefficient JRC, confining pressure, and permeability pressure. Finally, for a certain joint segment, the confining pressure is compensated and corrected based on the test results. The corrected confining pressure is substituted into the initial empirical relationship between the average permeability of the fracture and the roughness coefficient JRC, confining pressure, and permeability pressure to obtain the final target empirical relationship.

[0061] In this method, the initial fracture width of the fractured rock mass is corrected based on the displacement of rock blocks on both sides under different stress states. A relationship between the permeability of rough joints and the roughness of a single joint, the seepage pressure, and the confining pressure is proposed. Based on the proposed permeability calculation formula, different stress states are substituted into the formula for inverse calculation. It is found that the formula shows better agreement under different roughness coefficients and different seepage pressures. To improve the applicability of this method under different confining pressures, a second correction is made for different confining pressures, resulting in a new permeability formula. Through the above process, the simulation results are more accurate, a new empirical relationship between the permeability of rough joints and the roughness of a single joint, the seepage pressure, and the confining pressure is proposed, the permeability calculation formula is optimized, and the problem of poor accuracy in existing studies of the seepage characteristics of fractured rock masses is solved.

[0062] In some embodiments, step S101, determining the roughness coefficient of the actual joint of the target rock mass, includes: obtaining intermediate parameters based on the number of segments in the measurement interval of the target rock mass and the protrusion height of each measuring point; and substituting the intermediate parameters into a pre-set regression equation to determine the roughness coefficient of the actual joint of the target rock mass.

[0063] For example, in this invention, m=10, and the following calculations are performed:

[0064] First, the first derivative of the square root. The intermediate parameters are determined by the following formula:

[0065]

[0066] in, M The number of segments in the measurement section on the joint profile. Indicates the first i The height of the protrusion at each measuring point Indicates the first i The x-coordinates of each measuring point are obtained. After obtaining the intermediate parameters, they are substituted into the regression equation between the intermediate parameters and JRC:

[0067]

[0068] Finally, the roughness coefficient of the actual joint is calculated using an empirical formula, as follows:

[0069]

[0070] in, Represents the roughness coefficient of the actual joint. This represents the roughness coefficient of a target rock mass specimen with a length of 10 cm. It is 10cm. The value represents the actual joint length (cm), which is taken as 5cm in this embodiment. The calculation results of JRC in this embodiment are shown in Table 1:

[0071] Table 1 Results of actual joint roughness coefficient

[0072]

[0073] In some embodiments, step S102, which corrects and determines the actual permeability of a single fracture in the target rock mass based on the cubic law, includes: determining the relationship between the permeability of the target rock mass through the fracture surface and the fracture width under steady water flow based on the cubic law; and correcting and determining the actual permeability of a single fracture in the target rock mass based on the relationship between the permeability of the target rock mass through the fracture surface and the fracture width under steady water flow.

[0074] The cubic law assumes that the fluid is incompressible and derives the relationship between the amount of water seeping through a fracture surface and the fracture width under steady-state water flow. The expression is as follows:

[0075]

[0076] in, q This indicates the amount of permeation through the fracture surface. g Represents gravitational acceleration. v Represents the kinematic viscosity coefficient of a fluid. J Indicates the hydraulic gradient along the direction of the parallel plate. b This represents the crack width. According to Darcy's law, we have:

[0077]

[0078] in, K Let represent the permeability coefficient of a fracture in a smooth flat plate. The formula for calculating the permeability of an actual single fracture, corrected based on the cubic law, is as follows:

[0079]

[0080] in, w This indicates the transverse length of the crack.

[0081] Based on this embodiment, the initial width of the crack in the target rock mass is corrected to obtain the crack width after the crack displacement. This includes: determining the crack deformation of the target rock mass based on the hyperbolic deformation formula of the normal deformation stress change of the crack structure surface; and correcting the initial width of the crack in the target rock mass according to the crack deformation to obtain the crack width after the crack displacement.

[0082] Based on the hyperbolic deformation formula for the change of normal deformation stress on the fracture structure surface, the initial width of the fracture is corrected, and the fracture deformation is then... S is:

[0083]

[0084]

[0085] in, This represents the initial fracture width of a single-fracture rock mass, which can be calculated from the results of preliminary seepage experiments. This represents the normal stiffness coefficient of the fractured structural surface. This represents the initial normal stiffness coefficient of the fractured structural surface. Let be the increment of the normal stress on the crack surface relative to the initial state. Then the crack width after crack displacement is:

[0086]

[0087] in, This represents the normal strain increment of the fracture structure surface in a single-fracture rock mass.

[0088] In some embodiments, step S103, establishing a single-fracture rock mass calculation model under triaxial stress conditions, and obtaining displacement data of the target rock mass through the single-fracture rock mass calculation model, includes: constructing a single-fracture rock mass calculation model; treating the rock blocks on both sides of the fracture surface of the target rock mass as linear elastic bodies, and obtaining the total displacement of the single-fracture rock mass in the target direction through the single-fracture rock mass calculation model. Crack displacement and displacement of rock blocks on both sides .

[0089] For example, establish as Figure 2The single-fracture rock mass calculation model shown is as follows: Figure 2 This is a calculation model diagram of a single-fracture rock mass according to an embodiment of the present invention. In the diagram, Let be the stress in the Z direction experienced by the rock blocks on both sides of the fracture surface, b be the fracture width, and h be the rock block thickness. Then along... Total displacement of a single-fracture rock mass Crack displacement and displacement of rock blocks on both sides The calculation formula is as follows:

[0090]

[0091]

[0092]

[0093]

[0094] in, For single-fracture rock mass along Strain increment in direction For the rock blocks on both sides along Strain increment in direction For the fracture structure surface along Strain increment in the direction.

[0095] In some embodiments, step S104, determining the fracture width of the target rock mass based on the stress state of the target rock mass and the displacement data of the target rock mass, includes: obtaining the strain increment of the single fracture rock mass and the rock blocks on both sides of the target rock mass along the target direction based on the stress state of the target rock mass and the displacement data of the target rock mass; and determining the fracture width of the target rock mass based on the strain increment of the single fracture rock mass and the rock blocks on both sides of the target rock mass along the target direction.

[0096] Considering the confining pressure and seepage water pressure during the seepage process of a single-fracture granite sample, the strain increment of the single-fracture rock mass and the rock blocks on both sides of the fracture along the Z direction can be expressed as:

[0097]

[0098] in, This represents the strain increment of a single-fractured rock mass and the rock blocks on both sides of the fracture along the Z direction, where E represents the elastic modulus of the rock blocks on both sides of the fracture surface. This represents the stress in the Z-direction experienced by the rock blocks on both sides of the fractured structural plane. This represents the stress in the X direction experienced by the rock blocks on both sides of the fractured structural plane. Let be the stress in the Y direction experienced by the rock blocks on both sides of the fracture surface. Simplifying the above equation, we get:

[0099]

[0100] in, For the fracture structure surface along The strain increment in the direction, b is the fracture width, and h is the rock block thickness. The elastic modulus of the fractured rock mass. Poisson's ratio for fractured rock masses Poisson's ratio of the rock blocks on both sides of the fracture surface; and The calculation formula is as follows:

[0101]

[0102]

[0103] in, This represents the normal stiffness coefficient of the fractured structural surface.

[0104] In this embodiment, the stress state of the single-fracture rock mass is simplified. During the seepage process, it is only subjected to confining pressure and seepage water pressure. Only the normal stress acting on the fracture surface of the single-fracture rock mass is considered, and the effect of shear stress is not considered. Therefore, the stress state of the single-fracture rock mass in this embodiment is as follows:

[0105]

[0106] in, For confining pressure, Let be the osmotic pressure. Combining the above formulas, we can obtain:

[0107]

[0108] Simplifying the above equation, we get:

[0109]

[0110] in, The initial fracture width is denoted as . In this embodiment, the selected target rock mass sample has an elastic modulus E of 14.57 GPa, a Poisson's ratio of 0.17, and a fracture surface normal stiffness coefficient of 4.324 TPa·m. -1 .

[0111] Based on the above embodiments, step S105, combining the fracture width of the target rock mass and the actual permeability of a single fracture in the target rock mass, introduces the roughness coefficient of the actual joint of the target rock mass, and corrects and fits the relationship between the roughness coefficient and the permeability to obtain an initial empirical formula. This includes: obtaining the quantitative relationship between the flow rate of a single joint fracture in the target rock mass and roughness, confining pressure, and seepage pressure; introducing the roughness coefficient of the actual joint of the target rock mass; correcting the relationship between the roughness coefficient and the permeability; obtaining an image of the permeability and roughness coefficient of the target rock mass, and fitting the relationship between the permeability of a single joint fracture and the roughness coefficient to obtain a fitting coefficient; and determining the initial empirical formula based on the fitting coefficient.

[0112] For example, a JRC value is introduced, and it is corrected based on the results of permeation tests with different JRC values ​​under the same stress state, as shown in the following formula:

[0113]

[0114] in, , is the fitted curve of the roughness coefficient, where a, m and c are all fitting coefficients. Figure 3 This is a schematic diagram of the relationship between the permeation amount Q and JRC obtained in an embodiment of the present invention, as shown in the figure. Figure 3 As shown, the relationship between the permeability of a single joint fracture and the JRC can be fitted by the following formula:

[0115]

[0116] in, Let J be the seepage rate (m³ / s) of a single-jointed fracture under confining pressure σ = 5 MPa and seepage water pressure P = 0.3 MPa. J is the hydraulic gradient, and g is the acceleration due to gravity, taken as 9.8 m / s². 2 w is the transverse length of the fracture (m), b is the corrected fracture width (m), and v is the kinematic viscosity of water, which is taken as 1.007 × 10⁻⁶ in this embodiment. -6 m / s 2 a, m, and c are all fitting coefficients, and JRC is the roughness coefficient of a single joint. The specific values ​​are shown in Table 1.

[0117] The correlation coefficient R obtained from the fitting 2 =0.9868, indicating a relatively ideal fitting effect. The obtained nonlinear coefficients are a=2.565, m=0.07627, and c=1.525.

[0118] In summary, the initial empirical relationships between the average permeability of the fracture and the roughness coefficient JRC, confining pressure, and seepage pressure are as follows:

[0119]

[0120] Where Q represents the permeability of a single joint fracture. J This indicates the hydraulic gradient along the direction of the parallel plate.

[0121] Because the displacement generated by the structural planes on both sides of the fractured rock mass has a nonlinear relationship with the applied confining pressure, the degree of correction to the confining pressure is too large when the confining pressure is too small or too large. Therefore, it is necessary to compensate for the effect of the confining pressure.

[0122] This invention takes a certain joint segment to correct the confining pressure. For example, when taking When the confining pressure is at that time, the following formula is used to correct it:

[0123]

[0124] in, This represents the corrected confining pressure. The correlation coefficient R... 2 =0.9962, indicating a relatively good fit. The rationality of the correction function is verified by calculating the average permeability of the fracture surface using the corrected confining pressure. The calculation results are shown in Table 2 below:

[0125] Table 2 Comparison of average permeability of single fracture under different confining pressures

[0126]

[0127] As shown in the table above, the errors are 0.005%, 5.71%, 4.53%, and 0.945%, respectively, all around 5%, indicating a relatively good fitting effect. Figure 4 This is a comparative analysis chart of the experimental and corrected values ​​of permeability in embodiments of the present invention. Figure 4 It can be seen that adopting a shape like It is reasonable to correct the confining pressure in the equation, where A, B, and C are the influence coefficients of fracture roughness, obtained from experimental data. In summary, considering the corrected confining pressure, the target empirical relationship between the final average fracture permeability and the roughness coefficient JRC, confining pressure, and permeable water pressure is:

[0128]

[0129] in, This indicates the corrected confining pressure.

[0130] The present invention also provides a device for correcting the cubic law of seepage in coarse single-jointed channels of granite. The following is a description of the device for correcting the cubic law of seepage in coarse single-jointed channels of granite provided by the present invention. The device for correcting the cubic law of seepage in coarse single-jointed channels of granite described below can be referred to in correspondence with the method for correcting the cubic law of seepage in coarse single-jointed channels of granite described above. Figure 5This is a structural block diagram of the granite rough single-jointed channel seepage cubic law correction device provided by the present invention, as shown in the figure. Figure 5 As shown, the device includes:

[0131] The coefficient determination module 501 is used to divide the target rock mass into several measurement intervals and determine the roughness coefficient of the actual joints of the target rock mass;

[0132] The width correction module 502 is used to correct and determine the actual permeability of a single fracture in the target rock mass based on the cubic law, and to correct the initial width of the fracture in the target rock mass to obtain the fracture width after the fracture displacement in the target rock mass.

[0133] Data processing module 503 is used to establish a calculation model of a single-fracture rock mass under triaxial stress conditions and obtain displacement data of the target rock mass through the single-fracture rock mass calculation model;

[0134] The width determination module 504 is used to determine the crack width of the target rock mass based on the stress state of the target rock mass and the displacement data of the target rock mass.

[0135] The correction fitting module 505 is used to combine the fracture width of the target rock mass and the actual permeability of a single fracture in the target rock mass, introduce the roughness coefficient of the actual joint of the target rock mass, and correct and fit the relationship between the roughness coefficient and the permeability to obtain the initial empirical relationship.

[0136] The confining pressure correction module 506 is used to compensate and correct the confining pressure of the target rock mass, and to input the corrected confining pressure into the initial empirical relation to obtain the target empirical relation.

[0137] For example, firstly, the coefficient determination module 501 calculates the roughness coefficient of the target rock mass specimen with a length of m cm, and then obtains the joint roughness coefficient (JRC) of the actual joint using an empirical formula. Next, the width correction module 502 calculates the actual permeability of a single fracture based on the cubic law correction, and corrects the initial fracture width of the target rock mass, finally obtaining the fracture width after displacement. The data processing module 503 then establishes a single-fracture rock mass calculation model under triaxial stress conditions, and obtains the permeability along the target direction (…). Displacement data of a single-fracture rock mass (direction). The displacement data includes total displacement. Crack displacement and displacement of rock blocks on both sides The width determination module 504 then calculates the strain increment of the single-fracture rock mass and the rock blocks on both sides of the fracture along the Z direction by substituting the stress state of the fractured rock mass into the calculation formulas of the displacement data mentioned above. The fracture width is then calculated using the inverse formula. Next, the correction and fitting module 505 obtains the quantitative relationship between the flow rate of a single joint fracture and its roughness, confining pressure, and permeability pressure. The expression for the fracture width is substituted into the expression for the permeability, and the JRC value is introduced. Based on permeability test results with different JRC values ​​under the same stress state, the module corrects and fits the relationship. Finally, the fitting coefficients are substituted into the expression to obtain the initial empirical relationship between the average permeability of the fracture and the roughness coefficient JRC, confining pressure, and permeability pressure. Finally, the confining pressure correction module 506 compensates and corrects the confining pressure for a certain joint segment based on test results. The corrected confining pressure is substituted into the initial empirical relationship between the average permeability of the fracture and the roughness coefficient JRC, confining pressure, and permeability pressure to obtain the final target empirical relationship.

[0138] During operation, this device corrects the initial fracture width of the fractured rock mass based on the displacement of rock blocks on both sides under different stress states, and proposes a relationship between the permeability of rough joints and the roughness of a single joint, seepage pressure, and confining pressure. Based on the proposed permeability calculation formula, different stress states are substituted into the formula for back-calculation. It is found that the formula shows better agreement under different roughness coefficients and seepage pressures. To improve the applicability of this method under different confining pressures, a second correction is made for different confining pressures, resulting in a new permeability formula. Through the above process, the simulation results are more accurate, a new empirical relationship between the permeability of rough joints and the roughness of a single joint, seepage pressure, and confining pressure is proposed, the permeability calculation formula is optimized, and the problem of poor accuracy in existing studies of seepage characteristics of fractured rock masses is solved.

[0139] Figure 6 An example is a schematic diagram of the physical structure of an electronic device, such as... Figure 6 As shown, the electronic device may include: a processor 601, a communication interface 602, a memory 603, and a communication bus 604. The processor 601, communication interface 602, and memory 603 communicate with each other via the communication bus 604. The processor 601 can call logical instructions from the memory 603 to execute a method for correcting the cubic law of seepage through rough single-jointed channels in granite. This method includes:

[0140] The target rock mass is divided into several measurement intervals, and the roughness coefficient of the actual joints of the target rock mass is determined;

[0141] The permeability of a single fracture in the target rock mass is determined by correcting the cubic law, and the initial width of the fracture in the target rock mass is corrected to obtain the fracture width after the fracture displacement.

[0142] A calculation model of a single-fracture rock mass under triaxial stress conditions was established, and the displacement data of the target rock mass was obtained through the calculation model of the single-fracture rock mass.

[0143] Based on the stress state of the target rock mass and the displacement data of the target rock mass, the fracture width of the target rock mass is determined.

[0144] By combining the fracture width of the target rock mass and the actual permeability of a single fracture in the target rock mass, the roughness coefficient of the actual joints of the target rock mass is introduced, and the relationship between the roughness coefficient and the permeability is corrected and fitted to obtain the initial empirical relationship.

[0145] The confining pressure of the target rock mass is compensated and corrected, and the corrected confining pressure is then substituted into the initial empirical relation to obtain the target empirical relation.

[0146] Furthermore, the logical instructions in the aforementioned memory 603 can be implemented as software functional units and, when sold or used as independent products, can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods of the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0147] On the other hand, the present invention also provides a computer program product, which includes a computer program that can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer is able to execute the cubic law correction method for seepage through rough single jointed channels of granite provided by the above methods, the method comprising:

[0148] The target rock mass is divided into several measurement intervals, and the roughness coefficient of the actual joints of the target rock mass is determined;

[0149] The permeability of a single fracture in the target rock mass is determined by correcting the cubic law, and the initial width of the fracture in the target rock mass is corrected to obtain the fracture width after the fracture displacement.

[0150] A calculation model of a single-fracture rock mass under triaxial stress conditions was established, and the displacement data of the target rock mass was obtained through the calculation model of the single-fracture rock mass.

[0151] Based on the stress state of the target rock mass and the displacement data of the target rock mass, the fracture width of the target rock mass is determined.

[0152] By combining the fracture width of the target rock mass and the actual permeability of a single fracture in the target rock mass, the roughness coefficient of the actual joints of the target rock mass is introduced, and the relationship between the roughness coefficient and the permeability is corrected and fitted to obtain the initial empirical relationship.

[0153] The confining pressure of the target rock mass is compensated and corrected, and the corrected confining pressure is then substituted into the initial empirical relation to obtain the target empirical relation.

[0154] In another aspect, the present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the method for correcting the cubic law of seepage through rough single-jointed channels in granite provided by the methods described above, the method comprising:

[0155] The target rock mass is divided into several measurement intervals, and the roughness coefficient of the actual joints of the target rock mass is determined;

[0156] The permeability of a single fracture in the target rock mass is determined by correcting the cubic law, and the initial width of the fracture in the target rock mass is corrected to obtain the fracture width after the fracture displacement.

[0157] A calculation model of a single-fracture rock mass under triaxial stress conditions was established, and the displacement data of the target rock mass was obtained through the calculation model of the single-fracture rock mass.

[0158] Based on the stress state of the target rock mass and the displacement data of the target rock mass, the fracture width of the target rock mass is determined.

[0159] By combining the fracture width of the target rock mass and the actual permeability of a single fracture in the target rock mass, the roughness coefficient of the actual joints of the target rock mass is introduced, and the relationship between the roughness coefficient and the permeability is corrected and fitted to obtain the initial empirical relationship.

[0160] The confining pressure of the target rock mass is compensated and corrected, and the corrected confining pressure is then substituted into the initial empirical relation to obtain the target empirical relation.

[0161] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs. Those skilled in the art can understand and implement this without any creative effort.

[0162] Through the above description of the embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus necessary general-purpose hardware platforms, and of course, it can also be implemented by hardware. Based on this understanding, the above technical solutions, in essence or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute the methods described in the various embodiments or some parts of the embodiments.

[0163] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for modifying the cubic law for flow through a rough single-jointed channel in a granite, characterized by, The method comprises the following steps: dividing a target rock mass into several measurement intervals and determining the roughness coefficient of the actual joints of the target rock mass; correcting the permeability of the actual single fissure of the target rock mass based on the cubic law and correcting the initial width of the fissure of the target rock mass to obtain the fissure width after displacement of the target rock mass; establishing a single fissure rock mass calculation model under three-dimensional stress conditions and obtaining displacement data of the target rock mass through the single fissure rock mass calculation model; determining the fissure width of the target rock mass according to the stress state of the target rock mass and in combination with the displacement data of the target rock mass; combining the fissure width of the target rock mass and the permeability of the actual single fissure of the target rock mass, introducing the roughness coefficient of the actual joints of the target rock mass, correcting and fitting the relationship between the roughness coefficient and the permeability to obtain an initial empirical relationship; compensating and correcting the confining pressure of the target rock mass and bringing the corrected confining pressure into the initial empirical relationship to obtain a target empirical relationship.

2. The method of claim 1, wherein the method is a method of modifying the cubic law for flow through a rough single-jointed channel in a granite, characterized by, The method for determining the roughness coefficient of the actual joints of the target rock mass comprises the following steps: obtaining an intermediate parameter according to the number of segments of the measurement intervals of the target rock mass and the protrusion height of each measurement point; bringing the intermediate parameter into a pre-set regression equation to determine the roughness coefficient of the actual joints of the target rock mass.

3. The method for correcting the cubic law of seepage through rough single-jointed channels in granite according to claim 1, characterized in that, The method for correcting and determining the permeability of the actual single fissure of the target rock mass based on the cubic law comprises the following steps: determining a relationship between the permeability through a fissure face and the fissure width of the target rock mass under stable water flow based on the cubic law; correcting and determining the permeability of the actual single fissure of the target rock mass according to the relationship between the permeability through a fissure face and the fissure width of the target rock mass under stable water flow.

4. The method for correcting the cubic law of seepage through rough single-jointed channels in granite according to claim 1, characterized in that, The method for correcting the initial width of the fissure of the target rock mass to obtain the fissure width after displacement of the target rock mass comprises the following steps: determining the deformation of the fissure of the target rock mass based on a hyperbolic deformation formula of the normal deformation stress change of the fissure structure surface; correcting the initial width of the fissure of the target rock mass according to the deformation of the fissure of the target rock mass to obtain the fissure width after displacement of the target rock mass.

5. The method for correcting the cubic law of seepage through rough single-jointed channels in granite according to claim 1, characterized in that, The method for establishing a single fissure rock mass calculation model under three-dimensional stress conditions and obtaining displacement data of the target rock mass through the single fissure rock mass calculation model comprises the following steps: constructing the single fissure rock mass calculation model; regarding the rock masses on both sides of the fissure face of the target rock mass as linear elastic bodies and obtaining the total displacement of the single fissure rock mass of the target rock mass in a target direction, the fissure displacement and the displacement of the rock masses on both sides through the single fissure rock mass calculation model.

6. The method for correcting the cubic law of seepage through rough single-jointed channels in granite according to claim 1, characterized in that, The method for determining the fissure width of the target rock mass according to the stress state of the target rock mass and in combination with the displacement data of the target rock mass comprises the following steps: obtaining the strain increment of the single fissure rock mass and the rock masses on both sides of the target rock mass along a target direction according to the stress state of the target rock mass and in combination with the displacement data of the target rock mass; determining the fissure width of the target rock mass according to the strain increment of the single fissure rock mass and the rock masses on both sides of the target rock mass along a target direction.

7. The method for correcting the cubic law of seepage through rough single-jointed channels in granite according to claim 1, characterized in that, The roughness coefficient of the actual joint of the target rock mass is introduced to correct and fit the relationship between the roughness coefficient and the permeation amount, and an initial empirical relationship is obtained, including: The quantitative relationship between the single-joint fissure flow of the target rock mass and the roughness, confining pressure and permeation water pressure is obtained, the roughness coefficient of the actual joint of the target rock mass is introduced, and the relationship between the roughness coefficient and the permeation amount is corrected. An image between the permeation amount and the roughness coefficient of the target rock mass is obtained, and the relationship between the permeation amount and the roughness coefficient of the single-joint fissure is fitted to obtain a fitting coefficient. The initial empirical relationship is determined according to the fitting coefficient.

8. A device for correcting the cubic law of seepage in rough single-jointed channels of granite, characterized in that, The coefficient determination module is configured to divide the target rock mass into a plurality of measurement intervals and determine the roughness coefficient of the actual joint of the target rock mass. The width correction module is configured to correct the permeation amount of the actual single fissure of the target rock mass based on the cubic law and correct the initial width of the fissure of the target rock mass to obtain the fissure width after displacement of the target rock mass. The data processing module is configured to establish a single-fissure rock mass calculation model under a three-way stress condition and obtain displacement data of the target rock mass through the single-fissure rock mass calculation model. The width determination module is configured to determine the fissure width of the target rock mass according to the stress state of the target rock mass and the displacement data of the target rock mass. The correction and fitting module is configured to introduce the roughness coefficient of the actual joint of the target rock mass to correct and fit the relationship between the roughness coefficient and the permeation amount based on the fissure width of the target rock mass and the permeation amount of the actual single fissure of the target rock mass. The confining pressure correction module is configured to compensate and correct the confining pressure of the target rock mass, and introduce the corrected confining pressure into the initial empirical relationship to obtain a target empirical relationship. The processor executes the program to implement the granite rough single-joint channel seepage cubic law correction method of any one of claims 1 to 7.

9. An electronic device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, The computer program is executed by the processor to implement the granite rough single-joint channel seepage cubic law correction method of any one of claims 1 to 7. 10.A non-transitory computer-readable storage medium having stored thereon a computer program, characterized in that, ​

Images