Methods, apparatus and equipment for analyzing rock damage-permeability under explosive loading

Abstract

This invention provides a method, apparatus, and equipment for analyzing rock damage and permeability under explosive loading, relating to the field of geotechnical engineering. The invention utilizes finite element simulation to realistically reflect the propagation and response of explosive dynamic loads in rock mass, obtaining key principal strain data and establishing a first functional relationship between damage variables and the sum of tensile principal strains, thus intrinsically linking mechanical response with material damage. Furthermore, by constructing a second functional relationship between permeability and the sum of tensile principal strains, centered on undetermined key parameters, the invention incorporates damage variables and measured permeability observations into the inversion calculation, achieving objective and accurate calibration of key parameters. Finally, based on the undetermined key parameters, a quantitative coupling relationship between damage and permeability is established, enabling quantitative analysis of the permeability characteristics of rock mass after blasting, improving prediction accuracy, and providing a more scientific and precise analytical tool for engineering practices such as blasting permeability enhancement design and rock mass stability assessment.

Description

Technical Field

[0001] This invention relates to the field of geotechnical engineering technology, and in particular to a method, apparatus and equipment for analyzing rock damage-permeability under explosive loading. Background Technology

[0002] In geotechnical engineering projects such as mining, tunneling, in-situ leaching for uranium mining, and deep oil and gas resource development, blasting technology is one of the core methods for fracturing rock masses, forming desired engineering structures, or improving formation permeability. When explosive impact loads propagate through rock masses, they induce complex dynamic mechanical responses, leading to numerous microcracks, macroscopic fissures, and even destruction within the rock mass—a process collectively known as blast damage. The evolution of damage not only significantly degrades the mechanical strength and integrity of the rock mass but also fundamentally alters its permeability, directly impacting the safety and functionality of the project. For example, in uranium leaching, blasting is needed to create efficient seepage channels within the ore layer; in oil and gas reservoir enhancement, precise control of the fracture network formed by blasting is required to improve conductivity; and in tunnel excavation, the impact of blast damage zones on surrounding rock stability and groundwater seepage needs to be assessed. Therefore, accurately predicting the distribution of damage zones and the evolution of permeability in rock masses under explosive loads is a crucial theoretical basis for optimizing blasting design, evaluating engineering effectiveness, and mitigating disaster risks.

[0003] Currently, numerous numerical simulation methods and theoretical models have been developed for the response analysis of rock masses under dynamic loads such as explosions. In mechanical damage analysis, models based on continuum damage mechanics, fracture mechanics, or microstatistics can effectively describe crack initiation, propagation, and the deterioration process of the macroscopic strength of the rock mass. In seepage analysis, models based on Darcy's law or non-Darcy flow can be used to describe fluid flow in fractured media.

[0004] Most studies treat mechanical damage analysis and seepage characteristic assessment as relatively independent steps. A common approach is to first simulate the damage or fracture field caused by blasting, and then substitute it as a static condition into seepage calculations. This method fails to fully reflect the two-way coupling between the dynamic development of damage and real-time changes in permeability. In reality, crack initiation and propagation instantly alter local permeability, while changes in the seepage field (such as pore pressure) can also react on the crack tip, affecting its propagation behavior. This lack of coupling makes it difficult for models to accurately reflect the true physical picture of the spatiotemporal evolution of rock mass permeability characteristics after blasting.

[0005] Therefore, current methods for analyzing the permeability characteristics of rock masses after blasting suffer from low prediction accuracy. Summary of the Invention

[0006] This invention provides a method, apparatus, and equipment for analyzing rock damage and permeability under explosive loading, which solves the problem of low accuracy in predicting the permeability characteristics of rock mass after blasting.

[0007] In a first aspect, the present invention provides a method for analyzing rock damage and permeability under explosive loading. The method includes: obtaining mechanical response data of the target rock mass at various measuring points under explosive loading through finite element simulation of blasting; the mechanical response data includes principal strains in three mutually perpendicular directions; based on the mechanical response data of each measuring point, analyzing the mechanical condition and damage development state of each measuring point, and constructing a first functional relationship between the damage variable of the target rock mass and the sum of the tensile principal strains; using undetermined key parameters as independent variables and the permeability of the target rock mass as the dependent variable, establishing a second functional relationship between permeability and the sum of the tensile principal strains; based on the first and second functional relationships, as well as the damage variable and permeability observations at each measuring point, performing inversion calculations to determine the undetermined key parameters; and based on the undetermined key parameters, determining the quantitative coupling relationship between the damage variable and permeability.

[0008] Secondly, embodiments of the present invention provide an analysis device for rock damage-permeability under explosive loading. The device includes a communication module and a processing module. The communication module is used to acquire mechanical response data of the target rock mass at various measuring points under explosive loading through finite element simulation of blasting. The mechanical response data includes principal strains in three mutually perpendicular directions. The processing module is used to analyze the mechanical condition and damage development state of each measuring point based on the mechanical response data, and to construct a first functional relationship between the damage variable of the target rock mass and the sum of the tensile principal strains. Using undetermined key parameters as independent variables and the permeability of the target rock mass as the dependent variable, a second functional relationship between permeability and the sum of the tensile principal strains is established. Based on the first and second functional relationships, as well as the damage variable and permeability observations at each measuring point, inversion calculations are performed to determine the undetermined key parameters. Based on the undetermined key parameters, a quantitative coupling relationship between the damage variable and permeability is determined.

[0009] Thirdly, embodiments of the present invention provide an electronic device including a memory and a processor. The memory stores a computer program, and the processor is used to call and run the computer program stored in the memory to perform the steps of the first aspect and any possible implementation of the first aspect described above.

[0010] Fourthly, embodiments of the present invention provide a computer-readable storage medium storing a computer program, characterized in that, when the computer program is executed by a processor, it implements the steps of the method as described in the first aspect and any possible implementation thereof.

[0011] This invention provides a method, apparatus, and equipment for analyzing rock damage and permeability under explosive loading. The invention utilizes finite element simulation to realistically reflect the propagation and response of explosive dynamic loads in rock mass, obtaining key principal strain data and establishing a first functional relationship between damage variables and the sum of tensile principal strains, thus intrinsically linking macroscopic mechanical response with material damage. Furthermore, by constructing a second functional relationship between permeability and the sum of tensile principal strains, centered on undetermined key parameters, and incorporating damage variables and measured permeability observations into the inversion calculation, the invention achieves objective and accurate calibration of key parameters, avoiding errors caused by human experience-based values. Finally, based on the undetermined key parameters, a quantitative coupling relationship between damage and permeability is established, enabling quantitative analysis of the permeability characteristics of rock mass after blasting, improving prediction accuracy, and providing a more scientific and precise analytical tool for engineering practices such as blasting permeability enhancement design and rock mass stability assessment. Attached Figure Description

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

[0013] Figure 1 This is a schematic flowchart of a method for analyzing rock damage and permeability under explosive loading, provided in an embodiment of the present invention.

[0014] Figure 2 This is a schematic diagram of a method for constructing a rock damage-permeability coupling model according to an embodiment of the present invention;

[0015] Figure 3 This is a schematic flowchart of a method for determining key parameters to be determined according to an embodiment of the present invention;

[0016] Figure 4 This is a schematic diagram of the structure of a rock damage-permeability analysis device under explosive loading provided in an embodiment of the present invention;

[0017] Figure 5 This is a schematic diagram of the structure of an electronic device provided in an embodiment of the present invention. Detailed Implementation

[0018] In the following description, specific details such as particular system architectures and techniques are set forth for illustrative purposes and not for limitation, in order to provide a thorough understanding of the embodiments of the invention. However, those skilled in the art will understand that the invention can be implemented in other embodiments without these specific details. In other instances, detailed descriptions of well-known systems, apparatuses, circuits, and methods are omitted so as not to obscure the description of the invention with unnecessary detail.

[0019] In the description of this invention, unless otherwise stated, " / " means "or". For example, A / B can mean A or B. The term "and / or" in this document is merely a description of the relationship between related objects, indicating that three relationships can exist. For example, A and / or B can represent: A alone, A and B simultaneously, and B alone. Furthermore, "at least one" and "more than one" refer to two or more. The terms "first," "second," etc., do not limit the quantity or order of execution, and "first," "second," etc., do not necessarily imply differences.

[0020] In the embodiments of this application, the terms "exemplary" or "for example" are used to indicate that something is an example, illustration, or description. Any embodiment or design that is described as "exemplary" or "for example" in the embodiments of this application should not be construed as being more preferred or advantageous than other embodiments or design. Specifically, the use of terms such as "exemplary" or "for example" is intended to present the relevant concepts in a specific manner to facilitate understanding.

[0021] Furthermore, the terms "comprising" and "having," and any variations thereof, used in the description of this application are intended to cover non-exclusive inclusion. For example, a process, method, system, product, or device that includes a series of steps or modules is not limited to the steps or modules listed, but may optionally include other steps or modules not listed, or may optionally include other steps or modules inherent to such process, method, product, or device.

[0022] To make the objectives, technical solutions, and advantages of the present invention clearer, the following description will be provided in conjunction with the accompanying drawings and specific embodiments.

[0023] like Figure 1 As shown, this embodiment of the invention provides a method for analyzing rock damage and permeability under explosive loading. The method includes steps S101-S105.

[0024] S101. By finite element simulation of blasting, mechanical response data of the target rock mass at each measuring point under explosive load are obtained.

[0025] In this embodiment of the application, the mechanical response data includes principal strains in three mutually perpendicular directions.

[0026] For example, embodiments of the present invention can simulate the dynamic response process of rock mass under the action of explosive shock wave using finite element analysis software (such as LS-DYNA), and extract the three principal strains (ε1, ε2, ε3) of each integration point or element in the calculation model during the load history. For each principal strain, the tensile principal strains (i.e., principal strains with values ​​greater than zero) are selected. All selected tensile principal strains are algebraically summed to obtain the sum of the tensile principal strains at that point, θ, calculated using the formula: θ = Σ max(ε i , 0), where i = 1, 2, 3.

[0027] S102. Based on the mechanical response data of each measuring point, analyze the mechanical condition and damage development state of each measuring point, and construct the first functional relationship between the damage variable and the sum of the tensile principal strain of the target rock mass.

[0028] As one possible implementation, step S102 can be specifically implemented as steps S1021-S1025.

[0029] S1021. Based on the mechanical response data of each measuring point, calculate the sum of the principal strains of tension at each measuring point.

[0030] In some embodiments, the sum of the principal strains of tension is the sum of the principal strains of tension in three mutually perpendicular directions. The calculation method is to sum the positive values ​​of the three principal strains at each measuring point.

[0031] S1022. Based on the sum of the principal strains at each measuring point and the preset damage threshold, determine the damage development state of the target rock mass.

[0032] S1023. If the target rock mass is in a damaged state, the crack density of the target rock mass is determined based on the sum of the tensile principal strains at each measuring point and the preset damage threshold.

[0033] S1024. Based on the functional relationship between crack density and damage variables, and the crack density of the target rock mass, determine the damage variables at each measuring point of the target rock mass.

[0034] S1025. Based on the sum of the damage variables and the principal strain at each measuring point, establish the first functional relationship.

[0035] In some embodiments, the damage threshold θc is a key material parameter, which is the ratio of the dynamic tensile strength Tc of the rock to its elastic modulus E. When θ ≤ θc, the material is considered to be in an undamaged or damage-initiated state; when θ > θc, material damage is considered to have begun to initiate and develop as θ increases.

[0036] For example, such as Figure 2As shown, this invention provides a schematic flowchart of a method for constructing a rock damage-permeability coupling model. When θ ≤ θc, the crack density C is defined. d = 0, and correspondingly, the damage variable D = 0. When θ > θc, the crack density C d It is determined by the sum of the principal strains θ and its threshold θc. Among them, the crack density C... d The relationship with the sum of the principal strains θ during tensioning is as follows:

[0037] ;

[0038] Where m is a material constant (take m = 6), and the <> operator is <θ θc>= h(θ) θc)×(θ θc), h(θ) θc) is the Heaviside step function. When θ When θc is less than or equal to zero, h(θ) θc)=0; when θ When θc is greater than zero, h(θ) θc)=1. K IC R represents the fracture toughness of the undamaged material (a mechanical property index reflecting the material's ability to resist unstable crack propagation), where ρ is density and c is wave velocity. max The maximum strain rate experienced by the material is constant, and γk is combined into a single parameter. .

[0039] parameter The expression for a fixed Poisson's ratio v is: ;

[0040] in, 0 indicates that the parameter is 0 when v = 0. The value, 0 expression is ;

[0041] K is the stress intensity factor, and Ns is the crack shape parameter. For coin-shaped cracks, Ns = 1.12.

[0042] When θ When θc is less than zero, the material is considered to be in a state of no damage or damage not yet initiated; therefore, the crack density C... d The relationship with the sum of the principal strains θ during tensioning can be expressed as follows: ;

[0043] When θ When θc is greater than zero, it is considered that material damage begins to initiate and develops with increasing θ. Therefore, the crack density C dThe relationship with the sum of the principal strains θ during tensioning can be expressed as:

[0044] ;

[0045] Based on the above formula, the crack density C can be obtained. d Represented as ;

[0046] Damage variable D and crack density C d The relationship is ;

[0047] Therefore, the damage variable D is expressed as the relationship between the sum of the tensile principal strains θ, that is, the first functional relationship is: This ensures that the damage variable D ranges from 0 (no damage) to 1 (complete destruction), and has a clear physical meaning.

[0048] S103. Using the undetermined key parameters as independent variables and the permeability of the target rock mass as the dependent variable, establish a second functional relationship between permeability and the sum of tensile principal strain.

[0049] As one possible implementation, embodiments of the present invention can construct a second functional relationship based on the following formula.

[0050] ;

[0051] in, For penetration rate, The initial inherent permeability of the target rock mass. Key parameters to be determined The sum of the main strains of tension, To preset the damage threshold, This represents the critical value of tensile strain for permeability growth.

[0052] S104. Based on the first and second functional relationships, as well as the damage variables and permeability observations at each measuring point, inversion calculations are performed to determine the key parameters to be determined.

[0053] As one possible implementation, step S104 can be specifically implemented as steps S1041-S1043.

[0054] S1041. Based on the first and second functional relationships, construct a third functional relationship with the damage variable at each measuring point as the independent variable and the permeability calculation value as the dependent variable.

[0055] S1042. Construct an objective function that minimizes the overall error between the observed and calculated permeability values ​​at each measuring point.

[0056] S1043. Based on the damage variables and permeability observations at each measuring point, the third function relationship, and the objective function, iterative optimization is performed to obtain the undetermined key parameters that minimize the overall error.

[0057] For example, such as Figure 3 As shown, after determining the rock damage-permeability coupling model, the method for determining the key parameter C to be determined is an important part of this invention, and is determined through the following inversion process.

[0058] 1. Obtain experimental baseline data: Acquire measured permeability data of rock mass in a laboratory after being subjected to known explosive loads. This data should constitute a dataset. ;

[0059] Where N is the number of measurement points. The sum of the principal strains under tension is calculated through numerical simulation at the j-th measuring point. Let be the permeability observation value corresponding to the j-th measuring point.

[0060] 2. Construct an objective function to quantify the overall error between the model's predicted values ​​and the experimental data. Using the least squares method, construct the objective function J(C): ;

[0061] in, The permeability calculation value is derived from a theoretical model and depends on the undetermined key parameter C. This theoretical model is also known as the rock damage-permeability coupling model, or the third function relationship.

[0062] 3. Find the undetermined key parameter C that minimizes the objective function J(C) through inversion analysis. The specific process is as follows:

[0063] (1) Initialization: Set the initial guess value C0 of the undetermined key parameter C and the convergence threshold. And the maximum number of iterations Mmax.

[0064] (2) Single iteration calculation: For the k-th iteration, the current parameter value is C. k :

[0065] a. Model prediction: C k Substituting into the theoretical model, i.e., the third function relationship, calculate the permeability values ​​for all measuring points. .

[0066] b. Error Assessment: Calculate the current total error J(C) based on the objective function formula. k ).

[0067] c. Error assessment: Check J(C) k Is it less than the convergence threshold? Check if the iteration count k has reached Mmax. If either condition is met, terminate the iteration and let C = C k .

[0068] d. Parameter update: If the convergence condition is not met, a new parameter estimate C is generated. k+1 Then return to step a for the next iteration.

[0069] The obtained optimal undetermined key parameter C is substituted back into the third function relationship, and a complete numerical simulation is run. The inversion results are then finally verified by visually comparing the predicted permeability distribution with the spatial trend of experimental measurements. When the model prediction curve and the experimental data points achieve an optimal fit, the value of C can be finally confirmed.

[0070] S105. Based on the undetermined key parameters, determine the quantitative coupling relationship between damage variables and permeability.

[0071] As one possible implementation, embodiments of the present invention can determine the quantitative coupling relationship between the damage variables and permeability of the target rock mass based on undetermined key parameters and a third functional relationship.

[0072] This invention provides a method for analyzing rock damage and permeability under explosive loading. It utilizes finite element simulation to realistically reflect the propagation and response of explosive dynamic loads in rock mass, obtaining key principal strain data and establishing a first functional relationship between damage variables and the sum of tensile principal strains, thus intrinsically linking macroscopic mechanical response with material damage. Furthermore, this invention constructs a second functional relationship between permeability and the sum of tensile principal strains, centered on undetermined key parameters, and incorporates damage variables and measured permeability observations into the inversion calculation. This achieves objective and accurate calibration of key parameters, avoiding errors caused by human experience-based values. Finally, based on the undetermined key parameters, a quantitative coupling relationship between damage and permeability is established, enabling quantitative analysis of the permeability characteristics of rock mass after blasting. This improves prediction accuracy and provides a more scientific and precise analytical tool for engineering practices such as blasting permeability enhancement design and rock mass stability assessment.

[0073] Optionally, after step S105, the embodiments of the present invention can predict the spatial distribution of permeability of the target rock mass after the application of a set explosive load based on the quantitative coupling relationship; and determine the location result of the blasting permeability enhancement area based on the predicted spatial distribution of permeability.

[0074] Thus, this invention directly transforms the theoretical coupling model into an engineering prediction tool, which can quantitatively output the spatial distribution map of rock permeability after blasting, thereby accurately identifying areas with advantages in permeability enhancement. This provides an intuitive and reliable decision-making basis for the layout of uranium leaching injection holes and the optimization of oil and gas production enhancement blasting schemes, achieving a seamless connection from damage-seepage mechanism research to engineering design and effect prediction.

[0075] Optionally, before step S102, this embodiment of the invention may obtain the stress-strain curve of the target rock mass under dynamic load; based on the stress-strain curve, determine the dynamic tensile strength and dynamic elastic modulus of the target rock mass; calculate the ratio between the dynamic tensile strength and dynamic elastic modulus of the target rock mass, and determine it as a preset damage threshold.

[0076] Thus, this invention obtains dynamic tensile strength and elastic modulus through measured stress-strain curves, and defines a threshold based on their ratio, making the damage criterion have a clear physical meaning and experimental basis, significantly improving the rigor and universality of model parameter determination, and avoiding the arbitrariness of empirical values.

[0077] Optionally, the rock damage-permeability analysis method under explosive loading provided in this embodiment of the invention further includes steps A1-A5.

[0078] A1. Based on the quantitative coupling relationship, input the parameters of the proposed blasting scheme into the finite element model to simulate the explosion dynamic response and obtain simulation data.

[0079] In some embodiments, the simulation data includes the sum of the principal tensile strains of each computational unit in the finite element model during the explosion process.

[0080] A2. Substitute the sum of the tensile principal strains of each calculation unit during the explosion process into the quantitative coupling relationship, and iteratively calculate to obtain the permeability prediction value of each calculation unit.

[0081] A3. Based on the permeability prediction values ​​of each calculation unit, generate a spatial distribution map of the rock permeability after blasting.

[0082] A4. Evaluate the permeability enhancement effect of blasting based on the spatial distribution map.

[0083] A5. Based on the blasting permeability enhancement effect, adjust the blasting parameters, repeat the steps of finite element model simulation, permeability prediction, and blasting permeability enhancement effect evaluation, and perform iterative optimization until a permeability distribution that meets the set engineering requirements is obtained.

[0084] For example, after determining the quantitative coupling relationship between damage variables and permeability in step S105, the present invention can be used to predict the permeability enhancement effect of blasting in uranium ore formations and optimize the pore layout. The implementation process of the present invention is described in detail below with reference to a specific embodiment of blasting permeability enhancement design in uranium ore in-situ leaching mining.

[0085] Taking a specific mining area of ​​a sandstone-type uranium mine as the target rock mass, a finite element model reflecting its geological characteristics was established based on the geological exploration data and mechanical parameters of the mining area. The parameters of the proposed blasting scheme were input into the model, including: explosive type, single-hole charge amount, charge structure, spatial location, diameter, and depth of the blast holes, and the detonation method and delay sequence. Subsequently, the calibrated model parameters were input, including the undetermined key parameter C determined through inversion calculation in the previous steps, as well as the material's inherent initial permeability K0, damage threshold θc, and critical tensile strain θ for permeability increase. s , etc., are used as known constants to analyze the system.

[0086] Run finite element analysis software (such as LS-DYNA) to calculate the dynamic response under explosive loading. Extract the time history data of the three principal strains for all elements throughout the entire explosive process. For each element, calculate the sum of its tensile principal strains over the entire time history: θ = Σ max(ε i , 0), where i = 1, 2, 3. Where ε i This indicates that positive values ​​are taken (i.e., summation is only performed with respect to tensile strain).

[0087] Substitute the θ value of each unit into the established quantitative coupling relationship. Iteratively calculate the predicted permeability values ​​of all units in the model. Use post-processing software to generate a spatial distribution cloud map of the rock mass permeability after blasting.

[0088] Based on this result, it is possible to quantitatively assess whether the current blasting scheme can form a continuous and effective high-permeability zone in the uranium mine. If the target is not met, the blasting parameters can be adjusted in real time in the software and the simulation can be repeated until the optimal scheme is obtained. Finally, the spatial coordinates of the permeability enhancement zone output by the model can be directly used to guide the precise layout of injection wells and production wells, or to provide a scientific basis for the grouting reinforcement range of the surrounding rock of the roadway. This enables quantitative prediction and optimization of the blasting permeability enhancement effect before construction, and completely changes the traditional extensive design mode that relies on experience.

[0089] Thus, this invention deeply integrates the theoretical coupling model with the engineering blasting design process, enabling quantitative, visual prediction and dynamic optimization of the infiltration enhancement effect of any given blasting scheme. Through a closed-loop operation of parameter input, simulation calculation, distribution map generation and iterative adjustment, this method completely changes the traditional extensive design mode that relies on engineering experience, significantly improving the scientificity, reliability and economy of blasting infiltration enhancement schemes.

[0090] It should be understood that the sequence number of each step in the above embodiments does not imply the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present invention.

[0091] The following are device embodiments of the present invention. For details not described in detail, please refer to the corresponding method embodiments described above.

[0092] Figure 4 A schematic diagram of a rock damage-permeability analysis device under explosive loading, provided by an embodiment of the present invention, is shown. The analysis device 200 includes a communication module 201 and a processing module 202.

[0093] The communication module 201 is used to acquire mechanical response data of the target rock mass at various measuring points under explosive load through finite element simulation of blasting; the mechanical response data includes principal strains in three mutually perpendicular directions.

[0094] The processing module 202 is used to analyze the mechanical condition and damage development state of each measuring point based on the mechanical response data of each measuring point, and to construct the first functional relationship between the damage variables and the sum of the tensile principal strains of the target rock mass; using the undetermined key parameters as independent variables and the permeability of the target rock mass as the dependent variable, to establish the second functional relationship between the permeability and the sum of the tensile principal strains; based on the first and second functional relationships, as well as the damage variables and permeability observations of each measuring point, to perform inversion calculations to determine the undetermined key parameters; and based on the undetermined key parameters, to determine the quantitative coupling relationship between the damage variables and permeability.

[0095] Figure 5 This is a schematic diagram of the structure of an electronic device provided in an embodiment of the present invention. The electronic device 300 includes: a processor 301, a memory 302, and a computer program 303 stored in the memory 302 and executable on the processor 301. When the processor 301 executes the computer program 303, it implements the steps in the above-described method embodiments. Alternatively, when the processor 301 executes the computer program 303, it implements the functions of each module / unit in the above-described device embodiments.

[0096] For example, the computer program 303 may be divided into one or more modules / units, which are stored in the memory 302 and executed by the processor 301 to complete the present invention. The one or more modules / units may be a series of computer program instruction segments capable of performing a specific function, which describe the execution process of the computer program 303 in the electronic device 300.

[0097] The processor 301 may be a Central Processing Unit (CPU), or other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or any conventional processor.

[0098] The memory 302 can be an internal storage unit of the electronic device 300, such as a hard disk or memory of the electronic device 300. The memory 302 can also be an external storage device of the electronic device 300, such as a plug-in hard disk, smart media card (SMC), secure digital card (SD), flash card, etc., equipped on the electronic device 300. Furthermore, the memory 302 can include both internal and external storage units of the electronic device 300. The memory 302 is used to store the computer program and other programs and data required by the terminal. The memory 302 can also be used to temporarily store data that has been output or will be output.

[0099] The above-described embodiments are only used to illustrate the technical solutions of the present invention, and are not intended to limit it. 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. Such 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, and should all be included within the protection scope of the present invention.

Claims

1. A method for analyzing rock damage-permeability under explosive loading, characterized in that, include: The mechanical response data of the target rock mass at various measuring points under explosive load were obtained by finite element simulation of blasting. The mechanical response data includes principal strains in three mutually perpendicular directions; Based on the mechanical response data of each measuring point, the mechanical condition and damage development state of each measuring point are analyzed, and the first functional relationship between the damage variables and the sum of the tensile principal strains of the target rock mass is constructed. Using the undetermined key parameters as independent variables and the permeability of the target rock mass as the dependent variable, a second functional relationship between permeability and the sum of tensile principal strain is established. Based on the first and second functional relationships, as well as the damage variables and permeability observations at each measuring point, an inversion calculation is performed to determine the undetermined key parameters. This includes: constructing a third functional relationship based on the first and second functional relationships, with the damage variables at each measuring point as independent variables and the calculated permeability as dependent variables; constructing an objective function with the goal of minimizing the overall error between the observed and calculated permeability values ​​at each measuring point; and iteratively optimizing the solution based on the damage variables and permeability observations at each measuring point, the third functional relationship, and the objective function to obtain the undetermined key parameters that minimize the overall error. Based on the undetermined key parameters, the quantitative coupling relationship between damage variables and permeability is determined, including: based on the undetermined key parameters and the third functional relationship, the quantitative coupling relationship between damage variables and permeability of the target rock mass is determined.

2. The method for analyzing rock damage-permeability under explosive loading according to claim 1, characterized in that, Based on the mechanical response data from each measuring point, the mechanical condition and damage development state of each measuring point are analyzed to construct the first functional relationship between the damage variables and the sum of the tensile principal strains of the target rock mass, including: Based on the mechanical response data of each measuring point, the sum of the principal strains under tension at each measuring point is calculated; Based on the sum of the principal strains at each measuring point and the preset damage threshold, the damage development state of the target rock mass is determined. If the target rock mass is in a damaged state, the crack density of the target rock mass is determined based on the sum of the tensile principal strains at each measuring point and the preset damage threshold. Based on the functional relationship between crack density and damage variables, and the crack density of the target rock mass, the damage variables of each measuring point in the target rock mass are determined. The first functional relationship is established based on the sum of the damage variables and the principal strain at each measuring point.

3. The method for analyzing rock damage-permeability under explosive loading according to claim 1, characterized in that, The establishment of a second functional relationship between permeability and the sum of tensile principal strain, using undetermined key parameters as independent variables and the permeability of the target rock mass as the dependent variable, includes: The second functional relationship is constructed based on the following formula; ； in, For penetration rate, The initial inherent permeability of the target rock mass. Key parameters to be determined The sum of the main strains of tension, To preset the damage threshold, This represents the critical value of tensile strain for permeability growth.

4. The method for analyzing rock damage-permeability under explosive loading according to claim 1, characterized in that, After determining the quantitative coupling relationship between damage variables and permeability based on the undetermined key parameters, the process further includes: Based on the quantitative coupling relationship, the spatial distribution of permeability of the target rock mass after a set explosive load is predicted. Based on the predicted spatial distribution of permeability, the location of the blasting-enhanced permeability zone is determined.

5. The method for analyzing rock damage-permeability under explosive loading according to claim 1, characterized in that, Before analyzing the mechanical condition and damage development state of each measuring point based on the mechanical response data of each measuring point, and constructing the first functional relationship between the damage variables and the sum of the tensile principal strains of the target rock mass, the following steps are also included: Obtain the stress-strain curve of the target rock mass under dynamic loading; Based on the stress-strain curve, the dynamic tensile strength and dynamic elastic modulus of the target rock mass are determined; The ratio between the dynamic tensile strength and the dynamic elastic modulus of the target rock mass is calculated and determined as the preset damage threshold.

6. The method for analyzing rock damage-permeability under explosive loading according to claim 1, characterized in that, The method further includes: Based on the quantitative coupling relationship, the parameters of the proposed blasting scheme are input into the finite element model to simulate the explosion dynamic response and obtain simulation data. The simulation data includes the sum of the tensile principal strains of each calculation unit in the finite element model during the explosion process. Substitute the sum of the tensile principal strain of each calculation unit during the explosion process into the quantitative coupling relationship, and iteratively calculate to obtain the permeability prediction value of each calculation unit. Based on the permeability prediction values ​​of each calculation unit, a spatial distribution map of rock permeability after blasting is generated; Based on the aforementioned spatial distribution map, the permeability enhancement effect of blasting was evaluated; Based on the blasting permeability enhancement effect, adjust the blasting parameters, repeat the steps of finite element model simulation, permeability prediction, and blasting permeability enhancement effect evaluation, and perform iterative optimization until a permeability distribution that meets the set engineering requirements is obtained.

7. An analytical apparatus for rock damage-permeability under explosive loading, characterized in that, include: The communication module is used to acquire mechanical response data of the target rock mass at various measuring points under explosive load through finite element simulation of blasting; The mechanical response data includes principal strains in three mutually perpendicular directions; The processing module is used to analyze the mechanical condition and damage development state of each measuring point based on the mechanical response data of each measuring point, and construct a first functional relationship between the damage variables and the sum of the tensile principal strains of the target rock mass; using the undetermined key parameters as independent variables and the permeability of the target rock mass as the dependent variable, a second functional relationship between permeability and the sum of the tensile principal strains is established; based on the first and second functional relationships, as well as the damage variables and permeability observations of each measuring point, inversion calculations are performed to determine the undetermined key parameters; based on the undetermined key parameters, the quantitative coupling relationship between the damage variables and permeability is determined; The processing module is specifically used to construct a third functional relationship based on the first and second functional relationships, with the damage variables of each measuring point as independent variables and the calculated permeability as dependent variables; to construct an objective function with the goal of minimizing the overall error between the observed and calculated permeability values ​​of each measuring point; and to perform iterative optimization based on the damage variables and observed permeability values ​​of each measuring point, the third functional relationship, and the objective function to obtain the undetermined key parameters that minimize the overall error. The processing module is specifically used to determine the quantitative coupling relationship between the damage variables and permeability of the target rock mass based on the undetermined key parameters and the third functional relationship.

8. An electronic device, characterized in that, The electronic device includes a memory and a processor, the memory storing a computer program, and the processor being used to invoke and run the computer program stored in the memory to perform the method as described in any one of claims 1 to 6.

Images