Method and device for calculating elastic modulus of tension-compression dual-modulus rock and electronic equipment

Abstract

The invention relates to the technical field of rock mechanics and engineering, and discloses an elastic modulus calculation method and device for tension-compression dual-modulus rock and electronic equipment. The elastic modulus calculation method of the tension-compression dual-modulus rock comprises the following steps: constructing a mathematical model of rock stress based on a numerical model of rock tension-compression dual-modulus assumption; and calculating the elastic modulus of the rock based on the mathematical model of the rock stress. According to the method, the mathematical model of the rock stress is constructed based on the numerical model of the rock tension-compression dual-modulus assumption, so that the calculated stress can be more fit with the actual stress condition of the rock. And calculating the elastic modulus of the rock based on the mathematical model of the rock stress. Compared with a traditional stress calculation formula, the mathematical model of the stress is more real and more accurate, so that the elastic modulus calculated based on the mathematical model of the stress is more accurate and smaller in error, and the problems that in the prior art, the measurement precision of the tensile modulus and the compression modulus is not high, and the measurement error is large are solved.

Description

Technical Field

[0001] The present invention relates to the field of rock mechanics and engineering technology, and in particular to a method, a computing device and an electronic device for calculating the elastic modulus of a tensile and compressive dual modulus rock. Background Art

[0002] When rock is subjected to load, its primary response is deformation. As the load gradually increases, the deformation of the rock will become more and more significant until it reaches its ultimate bearing capacity and is destroyed. The deformation capacity of rock, as one of its core mechanical properties, plays a vital role in analyzing the stability of rock engineering and guiding engineering practice.

[0003] The Brazilian disk test, as a classic method for measuring the tensile strength of rocks, has been widely respected since its inception. The experiment is based on the isotropic elastic assumption, which simplifies the difficulty of analysis, but it implicitly assumes that the elastic properties of rocks under tension and compression are the same. However, this is not the case in reality. Most rocks exhibit different elastic responses under tension and compression, which is called the tensile-compression dual modulus characteristic. Specifically, the slope of the stress-strain curve of the rock during compression (i.e., the compression modulus) E c ) is often greater than its slope during stretching (i.e., the tensile modulus E t ). Different rock types E c / E t The ratio varies, with dense rocks being closer to 1:1, some sandstones being 2:1, and weaker rocks being as high as 10:1.

[0004] In view of the importance of the tensile and compressive dual modulus properties of rocks, some scholars have used Brazilian disk tests and other methods to measure the tensile and compressive elastic moduli of rocks. However, the accuracy of these methods still needs to be improved. The relative error of the compression modulus sometimes exceeds 10%, while the relative error of the tensile modulus is even more significant, exceeding 20%. Therefore, it is necessary to further explore more accurate measurement techniques to more accurately reveal the mechanical behavior of rocks. Therefore, how to improve the measurement accuracy of the tensile modulus and compression modulus of rocks and reduce the measurement error is an urgent problem that the industry needs to solve. Summary of the invention

[0005] The present invention provides a method, a computing device and an electronic device for calculating the elastic modulus of a tensile-compressive dual-modulus rock, which are used to solve the problems of low measurement accuracy and large measurement error of tensile modulus and compression modulus in the prior art.

[0006] The first aspect of the present invention provides a method for calculating the elastic modulus of a tensile and compressive dual modulus rock, comprising: Based on the numerical model of rock tension and compression dual modulus assumption, a mathematical model of rock stress is constructed; Based on the mathematical model of rock stress, the elastic modulus of the rock is calculated.

[0007] According to the elastic modulus calculation method of the tensile and compressive dual modulus rock provided by the present invention, in the step of calculating the elastic modulus of the rock based on the mathematical model of rock stress, the method comprises: The mathematical model based on rock stress, rock strain and Poisson's ratio under rock compression n c The elastic modulus of rock is calculated based on the relationship between them.

[0008] According to the method for calculating the elastic modulus of a tensile-compressive dual-modulus rock provided by the present invention, the elastic modulus includes the tensile modulus of the rock and / or the compression modulus of the rock.

[0009] According to the elastic modulus calculation method of tensile and compressive dual modulus rock provided by the present invention, the constitutive equation of the numerical model in the principal stress coordinate system is: Formula (1); Formula (1) is the stress tensor, is the strain tensor, is the fourth-order flexibility tensor in the principal stress coordinate system; Formula (2); Formula (3); In formula (2), E t represents the tensile modulus of rock, E c They represent the compression modulus of rock, n t represents the Poisson's ratio associated with the stretched state, n c represents the Poisson's ratio associated with the compressive state.

[0010] According to the elastic modulus calculation method of the tensile and compressive dual modulus rock provided by the present invention, in the step of constructing a mathematical model of rock stress based on the numerical model of the tensile and compressive dual modulus assumption of the rock, the method comprises: Based on the numerical model of rock tension and compression dual modulus assumption, parameter sensitivity analysis and data fitting are performed on the key parameters of rock to obtain parameter sensitivity analysis results and data fitting results; Based on the parameter sensitivity analysis results and the data fitting results, a mathematical model of rock stress is constructed.

[0011] According to the elastic modulus calculation method of tensile and compressive dual modulus rock provided by the present invention, in the step of constructing a mathematical model of rock stress based on the parameter sensitivity analysis result and the data fitting result, the method specifically includes: Based on the parameter sensitivity analysis result, constructing a first function of rock stress and the key parameter; Based on the data fitting result, obtaining a fitting function of the key parameters; Based on the first function and the fitting function, a mathematical model of the rock stress is constructed.

[0012] According to the elastic modulus calculation method of tensile and compressive dual modulus rock provided by the present invention, the step of constructing a first function of rock stress and the key parameter based on the parameter sensitivity analysis result includes: Based on the parameter sensitivity analysis results, construct a second function of a plurality of the key parameters; constructing a third function of rock stress and one of the key parameters; Based on the second function and the third function, a first function of rock stress and the key parameter is constructed.

[0013] According to the elastic modulus calculation method of the tensile and compressive dual modulus rock provided by the present invention, the key parameters include the Poisson's ratio of the rock under compression state. n c , rock strain, rock compression modulus E c and the tensile modulus of rock E t .

[0014] A second aspect of the present invention provides a device for calculating elastic modulus of tensile and compressive dual modulus rock, comprising: Model building module, used to build a mathematical model of rock stress based on the numerical model of rock tension and compression dual modulus assumption; The calculation module is used to calculate the elastic modulus of the rock based on the mathematical model of rock stress.

[0015] A third aspect of the present invention provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements a method for calculating the elastic modulus of tensile and compressive dual modulus rock as described in any one of the above items when executing the computer program.

[0016] The present invention provides a method, a computing device, and an electronic device for calculating the elastic modulus of a tensile and compressive dual modulus rock. By constructing a mathematical model of rock stress based on a numerical model of the tensile and compressive dual modulus assumption of rock, the calculated stress can be made to fit the actual stress condition of the rock more accurately and realistically. The elastic modulus of the rock is then calculated based on the mathematical model of rock stress. Because the mathematical model of stress is more realistic and accurate than the traditional stress calculation formula, the elastic modulus calculated based on the mathematical model of stress is also more accurate and has a smaller error, thereby solving the problem of low measurement accuracy and large measurement error of the tensile modulus and compression modulus in the prior art. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0018] Figure 1 It is one of the flow charts of the elastic modulus calculation method of the tensile and compressive dual modulus rock provided by the present invention.

[0019] Figure 2 This is the second flow chart of the elastic modulus calculation method for tensile and compressive dual modulus rock provided by the present invention.

[0020] Figure 3 This is a schematic diagram of the existing Brazilian disk experiment.

[0021] Figure 4 is a stress and strain sensitivity analysis diagram of the center point of the Brazilian disk of Example 1 of the present invention; wherein (a) E c and k = E c / E t The impact of n c =0.35); (b) k = E c / E t and n c The impact of E c =10GPa).

[0022] Figure 5 is the stress and e y / e x and ν c relationship diagram.

[0023] Figure 6 3 is a histogram of the calculation error of the elastic modulus of Example 1 of the present invention.

[0024] Figure 7 3 is a histogram of the calculation error of the elastic modulus of Comparative Example 1 of the present invention.

[0025] Figure 8 It is a structural schematic diagram of elastic modulus calculation of tensile and compressive dual modulus rock provided by the present invention.

[0026] Figure 9 It is a structural schematic diagram of the electronic device provided by the present invention. DETAILED DESCRIPTION

[0027] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention will be clearly and completely described below in conjunction with the drawings of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0028] In the embodiments of the present invention, "at least one" refers to one or more, and "more than one" refers to two or more. "And / or" describes the association relationship of associated objects, indicating that three relationships may exist. For example, A and / or B may represent: A exists alone, A and B exist at the same time, and B exists alone, where A and B may be singular or plural. In the textual description of the present invention, the character " / " generally indicates that the previous and next associated objects are in an "or" relationship. In addition, it should be noted that the serial numbers themselves, such as "first", "second", etc., used to distinguish the objects described in the present invention are only used to distinguish the objects described, and do not have any order or technical meaning.

[0029] Combine the following Figure 1 to Figure 9The elastic modulus calculation method, computing device and electronic device of the tensile and compressive dual modulus rock of the present invention are described in detail. Among them, the execution subject of the elastic modulus calculation method of the tensile and compressive dual modulus rock is an electronic device or a server. The electronic device can be a personal computer (PC), a portable device, a laptop, a smart phone, a tablet computer, a portable wearable device and other devices; the server can refer to a server, or a server cluster composed of multiple servers, a cloud server, etc. The present invention does not specifically limit the specific form of the electronic device or the server. Furthermore, the elastic modulus calculation method of the tensile and compressive dual modulus rock can also be applied to a computing device set in an electronic device or a server, and the computing device can be implemented by software, hardware or a combination of the two. The following takes the execution subject of the elastic modulus calculation method of the tensile and compressive dual modulus rock as a computing device as an example to describe the elastic modulus calculation method of the tensile and compressive dual modulus rock of the present invention.

[0030] In order to facilitate understanding of the elastic modulus calculation method for tensile and compressive dual modulus rock provided by the embodiment of the present invention, the elastic modulus calculation method for tensile and compressive dual modulus rock provided by the present invention is described in detail below through the following exemplary embodiments. It can be understood that the following exemplary embodiments can be combined with each other, and the same or similar concepts or processes may not be repeated in some embodiments.

[0031] like Figure 1 As shown, an embodiment of the first aspect of the present invention provides a method for calculating the elastic modulus of a tensile-compressive dual-modulus rock. The method for calculating the elastic modulus of a tensile-compressive dual-modulus rock comprises: S100, a numerical model based on the assumption of rock tension and compression dual modulus, to construct a mathematical model of rock stress; S200, based on the mathematical model of rock stress, calculate the elastic modulus of rock.

[0032] In this embodiment, by constructing a mathematical model of rock stress based on a numerical model based on the assumption of rock tension and compression dual modulus, the calculated stress can be more closely aligned with the actual stress condition of the rock and more accurate and true. Then, based on the mathematical model of rock stress, the elastic modulus of the rock is calculated. Because the mathematical model of stress is more realistic and accurate than the traditional stress calculation formula, the elastic modulus calculated based on the mathematical model of stress is also more accurate and has a smaller error, thus solving the problem of low measurement accuracy and large measurement error of tensile modulus and compression modulus in the prior art.

[0033] In some embodiments, the elastic modulus includes the tensile modulus of the rock and / or the compression modulus of the rock.

[0034] In some embodiments, the step of calculating the elastic modulus of rock based on the mathematical model of rock stress includes: Mathematical model based on rock stress, rock strain and Poisson's ratio under rock compression n c The elastic modulus of rock is calculated based on the relationship between them.

[0035] Specifically, the tensile modulus E t and compression modulus E c They can be calculated by formula (16) respectively.

[0036] Formula (16).

[0037] In formula (16), E t and E c Respectively represent the tensile modulus and compression modulus of rock; σ x and σ y are the horizontal stress component and the vertical stress component respectively; ε x and ε y is the strain component, which can be obtained through strain gauges; v c It represents the Poisson's ratio of rock under compression and can be obtained through uniaxial compression test.

[0038] In this embodiment, σ x and σ y They are calculated through the constructed mathematical model of stress.

[0039] In some embodiments, the constitutive equation of the numerical model in the principal stress coordinate system is: Formula (1); In formula (1), is the stress tensor, is the strain tensor, is the fourth-order flexibility tensor in the principal stress coordinate system; Formula (2); Formula (3); In formula (2), E t represents the tensile modulus of rock, E c They represent the compression modulus of rock, n t represents the Poisson's ratio associated with the stretched state, n c represents the Poisson's ratio associated with the compressive state.

[0040] like Figure 2 As shown, in some embodiments, the step of constructing a mathematical model of rock stress based on a numerical model of rock tension and compression dual modulus assumption includes: S110. Based on the numerical model of rock tension and compression dual modulus assumption, parameter sensitivity analysis and data fitting are performed on key parameters of rock to obtain parameter sensitivity analysis results and data fitting results; S120. Based on the parameter sensitivity analysis results and data fitting results, a mathematical model of rock stress is constructed.

[0041] In this embodiment, based on the numerical model, parameter sensitivity analysis and data fitting are performed on the key parameters of the rock; based on the parameter sensitivity analysis results and the data fitting results, a mathematical model of rock stress is constructed to obtain a stress distribution with high accuracy, solving the problem that the stress distribution at the center of the Brazilian disk of a tensile and compressive dual modulus rock cannot be accurately derived in the prior art. Based on the mathematical model of the stress, the elastic modulus is calculated, which can improve the measurement accuracy of the elastic modulus and reduce the measurement error.

[0042] It should be noted that the mathematical model of rock stress in this embodiment includes a mathematical model of horizontal stress and / or a mathematical model of vertical stress.

[0043] In some embodiments, the step of constructing a mathematical model of rock stress based on the parameter sensitivity analysis results and the data fitting results specifically includes: Based on the results of parameter sensitivity analysis, the first function of rock stress and key parameters is constructed; Based on the data fitting results, obtain the fitting function of key parameters; Based on the first function and the fitting function, a mathematical model of rock stress is constructed.

[0044] Exemplarily, the step of constructing a mathematical model of rock stress based on the parameter sensitivity analysis results and the data fitting results specifically includes: Based on the results of parameter sensitivity analysis, the first function of rock horizontal stress and key parameters is constructed; Based on the data fitting results, obtain the fitting function of key parameters; Based on the first function and the fitting function, a mathematical model of rock horizontal stress is constructed.

[0045] Exemplarily, the step of constructing a mathematical model of rock stress based on the parameter sensitivity analysis results and the data fitting results specifically includes: Based on the results of parameter sensitivity analysis, the first function of rock vertical stress and key parameters is constructed; Based on the data fitting results, obtain the fitting function of key parameters; Based on the first function and the fitting function, a mathematical model of rock vertical stress is constructed.

[0046] Furthermore, based on the parameter sensitivity analysis results, the step of constructing the first function of rock stress and key parameters includes: Based on the results of parameter sensitivity analysis, a second function of multiple key parameters is constructed; Construct a third function of rock stress and a key parameter; Based on the second function and the third function, the first function of rock stress and key parameters is constructed.

[0047] Preferably, the key parameters include Poisson's ratio of the rock under compression n c , rock strain, rock compression modulus E c and the tensile modulus of rock E t .

[0048] Furthermore, the step of constructing a second function of multiple key parameters based on the parameter sensitivity analysis results specifically includes: Based on the results of parameter sensitivity analysis, the relationship between parameter k and strain and Poisson's ratio under rock compression state is constructed. n c The second function, where the parameter k is the compression modulus E c and tensile modulus E t The absolute value of the ratio.

[0049] Exemplarily, based on the parameter sensitivity analysis results, the step of constructing a first function of rock stress and key parameters includes: Based on the results of parameter sensitivity analysis, the relationship between parameter k and strain and Poisson's ratio under rock compression state is constructed. n c The second function of Construct the third function of rock horizontal stress and parameter k; Based on the second function and the third function, the first function of rock horizontal stress and key parameters is constructed.

[0050] Exemplarily, based on the parameter sensitivity analysis results, the step of constructing a first function of rock stress and key parameters includes: Based on the results of parameter sensitivity analysis, the relationship between parameter k and strain and Poisson's ratio under rock compression state is constructed. n c The second function of Construct the third function of rock vertical stress and parameter k; Based on the second function and the third function, the first function of rock vertical stress and key parameters is constructed.

[0051] In some embodiments, in the step of performing parameter sensitivity analysis and data fitting on key parameters of rock and obtaining parameter sensitivity analysis results and data fitting results, the method of performing data fitting includes: A bivariate cubic polynomial function was used to fit the key parameters of the data.

[0052] Specifically, a cubic polynomial function is used to calculate the strain and Poisson's ratio of rock under compression. n c Perform data fitting.

[0053] Example 1

[0054] This embodiment provides a method for calculating the elastic modulus of a compressive dual-modulus rock, the method comprising the following steps: Step 1: The principle of the numerical model of the rock tension and compression dual modulus assumption in this embodiment is as follows: Based on the coupled constitutive model driven by stress based on the assumption of rock tension and compression dual modulus, the constitutive equation in the principal stress coordinate system is in the following form: Formula (1); Among them, i, j, k, l can be 1, 2 or 3 respectively.

[0055] In formula (1), is the stress tensor, is the strain tensor, is the fourth-order elastic tensor in the principal stress coordinate system; The value of the element is determined by the sign of the principal stress, specifically: Formula (2); Formula (3); In formula (2), E t represents the tensile modulus of rock, E c They represent the compression modulus of rock, n t represents the Poisson's ratio associated with the stretched state, n c represents the Poisson's ratio associated with the compressive state.

[0056] It is worth noting that according to the laws of thermodynamics, the following relationship can be derived: n t / E t = nc / E c .

[0057] Formula (1) represents the constitutive equation in the principal stress coordinate system. In order to transform this equation from the principal coordinate system O-α-β-γ to the global coordinate system Oxyz, the following relationship is applied:

[0058] The above formulas are solved using the finite element method. The elastic modulus and Poisson's ratio at each point in the model are not only determined by the basic properties of the material, but also affected by the current stress state at that point. Correspondingly, the stress and strain at each point are also affected by the local elastic constants. Therefore, there is a cross-coupling effect between the elastic constants and the stress state of the material. In order to ensure the convergence and accuracy of the numerical solution, an iterative stress transfer algorithm is implemented in the computational code. The convergence criterion is set as the Euclidean norm of the global elastic matrix reaching a preset tolerance threshold.

[0059] like Figure 4 The parameter sensitivity analysis results obtained by the above numerical model are given. The key parameters at this time include the Poisson's ratio under rock compression n c , the horizontal strain of rock e x , the vertical strain of rock e y , the compression modulus of rock E c and the tensile modulus of rock E t . For stress (including horizontal stress and vertical stress ) and strain variables were normalized to make it easier to compare the effects of different parameters. The analysis results show that stress (including horizontal stress and vertical stress ) changes significantly with the parameter k, however, it is not affected by E c and n c From a mathematical point of view, the horizontal stress and vertical stress can all be expressed as functions related only to the parameter k, in the following form: Formula (6); Formula (7).

[0060] also, e y / e x With the parameters k and n c changes with changes in E c Therefore, e y / e x It can be expressed as parameters k and n c Function: Formula (8).

[0061] It should be noted that for parameters k and n c For every unique combination of e y / e x The value corresponds to it, and vice versa. That is, e y / e x and n c Each unique combination of uniquely determines the parameter k. It reflects a one-to-one correspondence. Therefore, the parameter k can be expressed as e y / e x and n c There is an inverse function: Formula (9).

[0062] Substituting formula (9) into formula (6) and formula (7) respectively, we can obtain: Formula (10).

[0063] Among them, formula (10) is the function model of parameter sensitivity analysis results.

[0064] Step 3: Perform data fitting on the key parameters of rock based on the above numerical model.

[0065] Figure 5 Shows stress, e y / e x and n c3D scatter plot between . In this figure, discrete data points are presented in the form of solid circles and are marked as "value points" in the legend. After a detailed comparison and analysis of various curve fitting methods, a bivariate cubic polynomial function was selected to fit these data, and the specific form is as follows: Formula (11).

[0066] In formula (11) , , a i and b i are the fitting coefficients.

[0067] Determine the fitting coefficients through regression analysis based on discrete data a i and b i Formula (11) is the data fitting result.

[0068] Step 4: Based on the results of parameter sensitivity analysis and data fitting, a mathematical model of rock horizontal stress and a mathematical model of rock vertical stress are constructed.

[0069] Specifically, formula (11) is substituted into formula (10) to obtain the rock horizontal fitting stress and the rock numerical fitting stress, as shown in the following formula (12).

[0070] Formula (12).

[0071] Step 5: Calculate the elastic modulus of rock based on the mathematical model of rock stress.

[0072] In the Brazilian disk test, a thin circular rock disk is placed between two rigid platens and compressed so that the applied load acts along the diameter of the specimen, e.g. Figure 3 This test method induces a tensile stress state in the center area of ​​the disc specimen, which eventually causes the specimen to break and fail due to tension. It is generally believed that the horizontal force state and vertical stress state at the center point of a uniform isotropic Brazilian disc are: Formula (13-1); Formula (13-2).

[0073] In formula (13-1) and formula (13-2), P is the external load, D is the disk diameter, and L is the disk thickness.

[0074] According to the elastic theory, for the tensile-compressive dual modulus Brazilian disk, the strain constitutive equation is simplified as follows: In formula (14-1) and formula (14-2), E t and E c Respectively represent the tensile modulus and compression modulus of rock; σ x and σ y are the horizontal stress component and the vertical stress component, respectively, which can be calculated by formula (13-1) and formula (13-2). x and ε y is the strain component, which can be obtained through strain gauges; v c It represents the Poisson's ratio of rock under compression and can be obtained through uniaxial compression test.

[0075] After processing formula (13-1), formula (13-2), formula (14-1) and formula (14-2), we can get: Formula (16).

[0076] Combining equation (12) and equation (16), the known Poisson's ratio ν can be used c ,strain e x and e y And the fitting stress and The compression modulus and tensile modulus of the Brazilian specimens were determined using the formulas. The calculated values ​​of the compression modulus and tensile modulus obtained using these formulas are listed in Table 1. It can be seen that the calculated modulus closely matches the input modulus in all cases. Figure 6 A histogram of the relative error between the calculated modulus and the input modulus is given. In all cases, is always less than 0.5%, and Keep below 0.4% ( k = 1). This deviation occurs due to small differences between the fitted stress and the actual stress. k = 1. In other words, through parameter sensitivity analysis and stress fitting, we obtained a more accurate stress value at the center of the Brazilian disk, effectively solving the problem of large error in elastic modulus measurement.

[0077] Table 1 Calculation results obtained using formula (12)

[0078] Comparative Example 1

[0079] The center of the Brazilian disk of tensile and compressive dual modulus rock was selected as the calculation point to verify the above measurement method. Under the condition of keeping the applied load constant at 10kN, the strain value at the center of the disk was determined by existing numerical simulation methods. e x and e y The stress at the center of the disk is calculated using formula (13-1) and formula (13-2): s x and s y In order to verify the accuracy and reliability of this method, five different cases were considered, and all the calculation results are listed in detail in Table 2. In addition, histograms of the calculation errors of the tensile modulus and compression modulus were plotted, as shown in Table 2. Figure 7 As shown. It is found that there is a non-negligible difference between the calculated modulus and the input modulus. Compression modulus E c The maximum calculation error of the tensile modulus is as high as 18%. E t The maximum calculation error reaches 55%.

[0080] Table 2 Calculation results obtained using formula (13-1) and formula (13-2)

[0081] By comparing the calculation results of Example 1 with those of Comparative Example 1, it can be seen that the method of this embodiment can improve the measurement accuracy of the elastic modulus and reduce the measurement error.

[0082] like Figure 8 As shown, a specific embodiment of the second aspect of the present invention provides a device for calculating the elastic modulus of a tensile-compressive dual modulus rock. The device for calculating the elastic modulus of a tensile-compressive dual modulus rock includes a model building module and a calculation module; the model building module is used to build a mathematical model of rock stress based on a numerical model of the tensile-compressive dual modulus assumption of rock. The calculation module is used to calculate the elastic modulus of rock based on the mathematical model of rock stress.

[0083] Figure 9 An example of a physical structure diagram of an electronic device is shown in FIG. Figure 9 As shown, the electronic device may include: a processor 810, a communication interface 820, a memory 830 and a communication bus 840, wherein the processor 810, the communication interface 820 and the memory 830 communicate with each other through the communication bus 840. The processor 810 may call the logic instructions in the memory 830 to execute the elastic modulus calculation method of the tension-compression dual modulus rock, the method comprising: constructing a mathematical model of rock stress based on a numerical model of the tension-compression dual modulus assumption of rock; and calculating the elastic modulus of rock based on the mathematical model of rock stress.

[0084] In addition, the logic instructions in the above-mentioned memory 830 can be implemented in the form of a software functional unit and can be stored in a computer-readable storage medium when it is sold or used as an independent product. Based on this understanding, the technical solution of the present invention is essentially or the part that contributes to the prior art or the part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions to enable a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk and other media that can store program codes.

[0085] On the other hand, the present invention also provides a computer program product, which includes a computer program. The computer program can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer can execute the elastic modulus calculation method of tensile and compressive dual modulus rock provided by the above-mentioned methods. The method includes: constructing a mathematical model of rock stress based on a numerical model of the tensile and compressive dual modulus assumption of rock; and calculating the elastic modulus of rock based on the mathematical model of rock stress.

[0086] On the other hand, the present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, it is implemented to execute the elastic modulus calculation method of tensile and compressive dual modulus rock provided by the above-mentioned methods. The method includes: constructing a mathematical model of rock stress based on a numerical model of the tensile and compressive dual modulus assumption of rock; and calculating the elastic modulus of rock based on the mathematical model of rock stress.

[0087] The device embodiments described above are merely illustrative, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. Those of ordinary skill in the art may understand and implement it without creative work.

[0088] Through the description of the above implementation methods, those skilled in the art can clearly understand that each implementation method can be implemented by means of software plus a necessary general hardware platform, and of course, can also be implemented by hardware. Based on this understanding, the above technical solution is essentially or the part that contributes to the prior art can be embodied in the form of a software product, and the computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a disk, an optical disk, etc., including a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in each embodiment or some parts of the embodiments.

[0089] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for calculating the elastic modulus of tensile and compressive dual modulus rock, characterized in that: include: Based on the numerical model of rock tension and compression dual modulus assumption, a mathematical model of rock stress is constructed; Based on the mathematical model of rock stress, the elastic modulus of the rock is calculated.

2. The elastic modulus calculation method of tensile and compressive dual modulus rock according to claim 1 is characterized in that: The step of calculating the elastic modulus of rock based on the mathematical model of rock stress includes: The mathematical model based on rock stress, rock strain and Poisson's ratio under rock compression ν c The elastic modulus of rock is calculated based on the relationship between them.

3. The elastic modulus calculation method of tensile and compressive dual modulus rock according to claim 1 is characterized in that: The elastic modulus includes the tensile modulus of the rock and / or the compression modulus of the rock.

4. The elastic modulus calculation method of tensile and compressive dual modulus rock according to claim 1 is characterized in that: The constitutive equation of the numerical model in the principal stress coordinate system is: Formula (1): In formula (1), is the stress tensor, is the strain tensor, is the fourth-order flexibility tensor in the principal stress coordinate system; Formula (2): Formula (3): In formula (2), E t represents the tensile modulus of rock, E c They represent the compression modulus of rock, ν t represents the Poisson's ratio associated with the stretched state, ν c represents the Poisson's ratio associated with the compressive state.

5. The elastic modulus calculation method of tensile and compressive dual modulus rock according to any one of claims 1 to 4, characterized in that: The step of constructing a mathematical model of rock stress based on the numerical model of rock tension and compression dual modulus assumption includes: Based on the numerical model of rock tension and compression dual modulus assumption, parameter sensitivity analysis and data fitting are performed on the key parameters of rock to obtain parameter sensitivity analysis results and data fitting results; Based on the parameter sensitivity analysis results and the data fitting results, a mathematical model of rock stress is constructed.

6. The elastic modulus calculation method of tensile and compressive dual modulus rock according to claim 5 is characterized in that: The step of constructing a mathematical model of rock stress based on the parameter sensitivity analysis results and the data fitting results specifically includes: Based on the parameter sensitivity analysis result, constructing a first function of rock stress and the key parameter; Based on the data fitting result, obtaining a fitting function of the key parameters; Based on the first function and the fitting function, a mathematical model of the rock stress is constructed.

7. The elastic modulus calculation method of tensile and compressive dual modulus rock according to claim 6 is characterized in that: The step of constructing a first function of rock stress and the key parameter based on the parameter sensitivity analysis result includes: Based on the parameter sensitivity analysis results, construct a second function of a plurality of the key parameters; constructing a third function of rock stress and one of the key parameters; Based on the second function and the third function, a first function of rock stress and the key parameter is constructed.

8. The elastic modulus calculation method of tensile and compressive dual modulus rock according to claim 5, characterized in that: The key parameters include Poisson's ratio under rock compression ν c , rock strain, rock compression modulus E c and the tensile modulus of rock E t .

9. A device for calculating elastic modulus of tensile and compressive dual modulus rock, characterized in that: include: Model building module, used to build a mathematical model of rock stress based on the numerical model of rock tension and compression dual modulus assumption; The calculation module is used to calculate the elastic modulus of the rock based on the mathematical model of rock stress.

10. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the computer program, the method for calculating the elastic modulus of tensile and compressive dual modulus rock as described in any one of claims 1 to 8 is implemented.

Images