Rock mass viscosity inversion method and device based on sound wave detection

Abstract

The invention relates to the technical field of geotechnical engineering, in particular to a rock mass viscosity inversion method and device based on sound wave detection, and the method comprises the steps: obtaining sound wave detection data of engineering exploration, and obtaining the elevation distribution rule of rock mass sound waves according to the sound wave detection data; constructing a rock mass viscosity inversion equation set according to the elevation distribution rule of the rock mass sound waves, and solving the rock mass viscosity inversion equation set to invert rock mass viscosity parameters; and determining a distribution rule of the rock mass viscosity along the elevation based on the rock mass viscosity parameter so as to determine a rock mass comprehensive viscosity coefficient of the rock mass in the target range. Therefore, the problems that the viscosity data of the rock mass cannot be obtained before field construction and the viscosity of the rock mass in the whole construction area cannot be reflected due to the fact that a rock mass viscosity testing method in the related technology needs to be combined with field production and reflects the viscosity characteristics of the rock mass in a certain fixed range are solved.

Description

Technical Field

[0001] The present application relates to the field of geotechnical engineering technology, and in particular to a rock viscosity inversion method and device based on acoustic wave detection. Background Art

[0002] The viscosity parameters of rock masses are fundamental for studying the changing patterns of stress wave propagation attenuation and rock deformation under dynamic loads. Currently, measurements of rock viscosity primarily involve indoor specimen testing and in-situ testing. In-situ testing is generally based on the propagation attenuation of stress waves in rock masses, using energy dissipation or attenuation coefficients to invert rock viscosity. This method requires an artificial excitation source, which increases field testing costs. Furthermore, the rock mass tested is limited to the rock mass between the excitation source and the measurement point, failing to reflect the viscosity characteristics of the rock mass throughout the construction area.

[0003] In related technologies, many engineering practitioners and scholars have made innovative contributions to rock viscosity. For example, a wave-based test method for the viscosity coefficient of jointed rock has been developed. This method proposes a method for inverting rock viscosity based on the amplitude attenuation law of wavelets of different frequencies as they propagate through jointed rock. By decomposing the vibration waveform of a measuring point into wavelets, the energy loss between measuring points is calculated, and the viscosity coefficient of the rock mass is then calculated. This method improves the stability of frequency-domain inversion of rock viscosity to a certain extent.

[0004] However, the test data in the testing methods of related technologies require on-site in-situ blasting tests, which often need to be combined with on-site production. This means that rock viscosity data cannot be obtained before on-site construction. On the other hand, the use of fluctuation test data reflects the viscosity characteristics of rock within a fixed range, which is a comprehensive parameter. For the more complex rock composition in actual projects, it still cannot reflect the viscosity of the rock in the entire construction area, and urgently needs to be improved. Summary of the Invention

[0005] The present application provides a rock viscosity inversion method and device based on acoustic wave detection to solve the problem that the rock viscosity testing method in the related art needs to be combined with on-site production and reflects the viscosity characteristics of the rock mass in a fixed range, which leads to the inability to obtain the viscosity data of the rock mass before on-site construction and the inability to reflect the viscosity of the rock mass in the entire construction area.

[0006] The first aspect of the present application provides a rock viscosity inversion method based on acoustic wave detection, comprising the following steps: obtaining acoustic wave detection data from engineering exploration, and obtaining the elevation distribution law of rock sound waves based on the acoustic wave detection data; constructing a rock viscosity inversion equation group based on the elevation distribution law of the rock sound waves, and solving the rock viscosity inversion equation group to invert rock viscosity parameters; determining the distribution law of rock viscosity along the elevation based on the rock viscosity parameters to determine the comprehensive rock viscosity coefficient of the rock mass in the target range.

[0007] Through the above technical solution, the embodiment of the present application can carry out acoustic detection by drilling acoustic logging in the rock mass to be tested, and construct a rock viscosity inversion equation group based on the detection data. Compared with the existing viscosity coefficient inversion method, it can make full use of engineering geological exploration data in the early stage without the need to carry out on-site in-situ blasting tests; and can further determine the distribution law of rock viscosity along the elevation based on the rock viscosity parameters to determine the comprehensive viscosity coefficient of the rock mass in the target range, taking into account the viscosity difference of the rock mass at different elevations, and calculating the comprehensive viscosity coefficient of the rock mass in different ranges. Not only can the rock viscosity coefficient be obtained before construction, reducing the workload of subsequent tests, but it can also adapt to the rock viscosity parameters required for construction in different parts, and provide a scientific basis for rock viscosity in subsequent construction and scientific research.

[0008] Optionally, in one embodiment of the present application, the rock viscosity inversion equations are determined according to the relationship between the acoustic wave velocity and the elastic wave velocity of the rock formation, and are in the form of:

[0009]

[0010] Among them, Cp-so1 represents the distribution law of the first group of frequency sound wave velocity with elevation, Cp-so2 represents the distribution law of the second group of frequency sound wave velocity with elevation, and ω1 and ω2 represent the two groups of sound wave test frequencies respectively.

[0011] Through the above technical solution, the embodiment of the present application can reflect the relationship between the elevation distribution law of the acoustic wave velocity, the law of change of the longitudinal wave viscosity coefficient of the rock mass with elevation, the acoustic wave test frequency, and the law of change of the elastic wave velocity of the rock mass with elevation based on the relationship between the acoustic wave velocity and the elastic wave velocity of the rock mass, and thus can solve unknown parameters when other parameters are known.

[0012] Optionally, in one embodiment of the present application, the elastic wave velocity is obtained from the static physical and mechanical parameters of the rock sample, wherein the calculation formula of the static physical and mechanical parameters of the rock sample is:

[0013]

[0014] Among them, C pIt represents the law of change of rock mass elastic wave velocity with elevation, λ and μ are the law of change of rock mass elastic Lame constant with elevation, and ρ is the law of change of rock mass density with elevation.

[0015] Through the above technical solution, the embodiment of the present application can determine the law of change of the elastic wave velocity of the rock mass with elevation, and use the determined law of change of the elastic wave velocity of the rock mass with elevation as a known formula combined with the rock viscosity inversion equation group to calculate the law of change of the longitudinal wave viscosity coefficient of the rock mass with elevation, so as to solve the comprehensive viscosity coefficient of the rock mass in the target range.

[0016] Optionally, in one embodiment of the present application, the calculation formula for the comprehensive viscosity coefficient of the rock mass in the target range is:

[0017]

[0018] Among them, η ps It represents the comprehensive longitudinal wave viscosity coefficient of the rock mass, and h1 and h2 represent the range of the rock mass to be calculated.

[0019] Through the above technical solution, the embodiment of the present application can calculate the comprehensive viscosity coefficient of the rock mass in different parts, take into account the viscosity differences of the rock mass at different elevations, and well reflect the viscosity of each layer of rock mass, which can be applied to the viscosity inversion of various complex rock masses in actual engineering.

[0020] The second aspect of the present application provides a rock viscosity inversion device based on acoustic wave detection, including: an acquisition module for acquiring acoustic wave detection data of engineering exploration, and obtaining the elevation distribution law of rock sound waves based on the acoustic wave detection data; a construction module for constructing a rock viscosity inversion equation group according to the elevation distribution law of the rock sound waves, and solving the rock viscosity inversion equation group to invert the rock viscosity parameters; an inversion module for determining the distribution law of rock viscosity along the elevation based on the rock viscosity parameters to determine the comprehensive rock viscosity coefficient of the rock in the target range.

[0021] Through the above technical solution, the embodiment of the present application can carry out acoustic detection by drilling acoustic logging in the rock mass to be tested, and construct a rock viscosity inversion equation group based on the detection data. Compared with the existing viscosity coefficient inversion method, it can make full use of engineering geological exploration data in the early stage without the need to carry out on-site in-situ blasting tests; and can further determine the distribution law of rock viscosity along the elevation based on the rock viscosity parameters to determine the comprehensive viscosity coefficient of the rock mass in the target range, taking into account the viscosity difference of the rock mass at different elevations, and calculating the comprehensive viscosity coefficient of the rock mass in different ranges. Not only can the rock viscosity coefficient be obtained before construction, reducing the workload of subsequent tests, but it can also adapt to the rock viscosity parameters required for construction in different parts, and provide a scientific basis for rock viscosity in subsequent construction and scientific research.

[0022] Optionally, in one embodiment of the present application, the rock viscosity inversion equations are determined according to the relationship between the acoustic wave velocity and the elastic wave velocity of the rock formation, and are in the form of:

[0023]

[0024] Among them, Cp-so1 represents the distribution law of the first group of frequency sound wave velocity with elevation, Cp-so2 represents the distribution law of the second group of frequency sound wave velocity with elevation, and ω1 and ω2 represent the two groups of sound wave test frequencies respectively.

[0025] Through the above technical solution, the embodiment of the present application can reflect the relationship between the elevation distribution law of the acoustic wave velocity, the law of change of the longitudinal wave viscosity coefficient of the rock mass with elevation, the acoustic wave test frequency, and the law of change of the elastic wave velocity of the rock mass with elevation based on the relationship between the acoustic wave velocity and the elastic wave velocity of the rock mass, and thus can solve unknown parameters when other parameters are known.

[0026] Optionally, in one embodiment of the present application, the elastic wave velocity is obtained from the static physical and mechanical parameters of the rock sample, wherein the calculation formula of the static physical and mechanical parameters of the rock sample is:

[0027]

[0028] Among them, C p It represents the law of change of rock mass elastic wave velocity with elevation, λ and μ are the law of change of rock mass elastic Lame constant with elevation, and ρ is the law of change of rock mass density with elevation.

[0029] Through the above technical solution, the embodiment of the present application can determine the law of change of the elastic wave velocity of the rock mass with elevation, and use the determined law of change of the elastic wave velocity of the rock mass with elevation as a known law combined with the rock mass viscosity inversion equation group to calculate the law of change of the longitudinal wave viscosity coefficient of the rock mass with elevation, so as to solve the comprehensive viscosity coefficient of the rock mass.

[0030] Optionally, in one embodiment of the present application, the calculation formula for the comprehensive viscosity coefficient of the rock mass in the target range is:

[0031]

[0032] Among them, η ps It represents the comprehensive longitudinal wave viscosity coefficient of the rock mass, and h1 and h2 represent the range of the rock mass to be calculated.

[0033] Through the above technical solution, the embodiment of the present application can calculate the comprehensive viscosity coefficient of the rock mass in different parts, take into account the viscosity differences of the rock mass at different elevations, and well reflect the viscosity of each layer of rock mass, which can be applied to the viscosity inversion of various complex rock masses in actual engineering.

[0034] The third aspect of the present application 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 executes the program to implement the rock viscosity inversion method based on acoustic wave detection as described in the above embodiment.

[0035] A fourth aspect of the present application provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the above-mentioned rock viscosity inversion method based on acoustic wave detection.

[0036] A fifth aspect of the present application provides a computer program product, which stores a computer program that, when executed by a processor, implements the above-mentioned rock viscosity inversion method based on acoustic wave detection.

[0037] The embodiment of the present application can obtain the elevation distribution law of rock mass acoustic waves through the acoustic wave detection data obtained from on-site geological exploration, so as to make full use of on-site engineering survey data without the need to carry out other additional tests; and then construct a rock mass viscosity inversion equation group, solve the equation group to realize the inversion of rock mass viscosity parameters, determine the distribution law of rock mass viscosity along the elevation, and finally determine the comprehensive viscosity coefficient of the rock mass according to the research scope, providing a scientific basis for rock mass viscosity in subsequent construction and scientific research. Thus, it solves the problem that the rock mass viscosity test method of the related technology needs to be combined with on-site production and reflects the viscosity characteristics of the rock mass in a fixed range, which leads to the problem that the viscosity data of the rock mass cannot be obtained before on-site construction and the viscosity of the rock mass in the entire construction area cannot be reflected.

[0038] Additional aspects and advantages of the present application will be given in part in the description below, and in part will become apparent from the description below, or will be learned through practice of the present application. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] The above and / or additional aspects and advantages of the present application will become apparent and easily understood from the following description of the embodiments in conjunction with the accompanying drawings, in which:

[0040] Figure 1 A flowchart of a rock viscosity inversion method based on acoustic wave detection according to an embodiment of the present application;

[0041] Figure 2 A schematic diagram of an engineering layered rock mass according to a specific embodiment of the present application;

[0042] Figure 3 Schematic diagram of the structure of a rock viscosity inversion device based on acoustic wave detection according to an embodiment of the present application;

[0043] Figure 4 The figure is a schematic diagram of the structure of an electronic device provided according to an embodiment of the present application.

[0044] Reference numerals:

[0045] 1- explosion source, 2- rock layer 1, 3- rock layer 2, 4- rock layer 3, 5- ground surface, 6- rock layer interface;

[0046] 10- Rock viscosity inversion device based on acoustic wave detection: 100- acquisition module, 200- construction module and 300- inversion module; 401- memory, 402- processor and 403- communication interface. DETAILED DESCRIPTION

[0047] The following describes in detail embodiments of the present application, examples of which are shown in the accompanying drawings, wherein the same or similar reference numerals throughout represent the same or similar elements or elements having the same or similar functions. The embodiments described below with reference to the accompanying drawings are exemplary and are intended to be used to explain the present application, and should not be construed as limiting the present application.

[0048] The following describes the rock viscosity inversion method and device based on acoustic wave detection of an embodiment of the present application with reference to the accompanying drawings. In view of the fact that the rock viscosity testing method of the related technology mentioned in the background technology center needs to be combined with on-site production and reflects the viscosity characteristics of the rock in a fixed range, which leads to the problem that the viscosity data of the rock cannot be obtained before on-site construction and the viscosity of the rock in the entire construction area cannot be reflected, the present application provides a rock viscosity inversion method based on acoustic wave detection. In this method, the elevation distribution law of the rock sound wave can be obtained through the acoustic wave detection data obtained from on-site geological exploration, so as to make full use of the on-site engineering survey data without the need to carry out additional tests; and then construct a rock viscosity inversion equation group, solve the equation group to realize the inversion of rock viscosity parameters, determine the distribution law of rock viscosity along the elevation, and finally determine the comprehensive viscosity coefficient of the rock according to the research scope, so as to provide a scientific basis for rock viscosity in the subsequent construction and scientific research process. This solves the problem that the rock viscosity testing method in the relevant technology needs to be combined with on-site production and reflects the viscosity characteristics of the rock mass in a fixed range, which leads to the inability to obtain rock viscosity data before on-site construction and the inability to reflect the viscosity of the rock mass in the entire construction area.

[0049] Specifically, Figure 1 A schematic flow chart of a rock viscosity inversion method based on acoustic wave detection provided in an embodiment of the present application.

[0050] like Figure 1 As shown in FIG, the rock viscosity inversion method based on acoustic wave detection includes the following steps:

[0051] In step S101, acoustic wave detection data of engineering exploration is obtained, and the elevation distribution law of rock mass acoustic waves is obtained based on the acoustic wave detection data.

[0052] It is understandable that the acoustic wave detection data includes but is not limited to waveform records, time series data, acoustic wave velocity, amplitude information, etc. For example, the amplitude information can reflect the amplitude change of the acoustic wave signal, and thus reflect the absorption or scattering characteristics of the medium to the acoustic wave.

[0053] During the actual implementation process, the embodiment of the present application can obtain acoustic wave detection data by conducting acoustic wave detection on acoustic wave logging drilled in the rock mass to be tested. Specifically, the embodiment of the present application can carry out two groups of acoustic wave tests at different frequencies, and the frequencies can be 20kHz and 30kHz respectively; those skilled in the art should understand that some other embodiments can also select different frequencies for acoustic wave testing, which can be specifically set in combination with actual scenarios.

[0054] Through the above technical solution, the embodiment of the present application can obtain acoustic wave detection data by performing acoustic wave detection on the rock mass to be tested, and reflect the elevation distribution law of the rock mass acoustic waves based on the acquired acoustic wave detection data, so as to better construct the rock mass viscosity inversion equation group in step S102, which can make full use of on-site engineering survey data without the need for additional tests, and can adapt to the complex composition of the rock mass, providing technical support for subsequent rock mass-related research.

[0055] In step S102, a rock viscosity inversion equation group is constructed according to the elevation distribution law of rock mass acoustic waves, and the rock viscosity inversion equation group is solved to invert the rock mass viscosity parameters.

[0056] Among them, after carrying out two groups of acoustic wave tests at different frequencies respectively, the embodiment of the present application can divide the elevation distribution law of rock acoustic waves into the first group of frequency acoustic wave velocity distribution law with elevation and the second group of frequency acoustic wave velocity distribution law with elevation, and construct a rock viscosity inversion equation group based on the first group of frequency acoustic wave velocity distribution law with elevation and the second group of frequency acoustic wave velocity distribution law with elevation combined with other parameter characteristics.

[0057] As a feasible approach, the construction process of the rock viscosity inversion equation group is as follows:

[0058] Each layer of rock mass is considered as a three-dimensional Kelvin viscoelastic body, including the longitudinal wave viscosity coefficient η p and shear wave viscosity coefficient η s The two parameters, the longitudinal wave viscosity coefficient and the transverse wave viscosity coefficient of the rock mass, satisfy the relationship described in formula (1):

[0059]

[0060] Among them, η p is the longitudinal wave viscosity coefficient, which can reflect the change law of the longitudinal wave viscosity coefficient of the rock mass with elevation.

[0061] The elastic wave velocity of each rock layer is calculated according to the static physical and mechanical parameters of the rock sample according to formula (2):

[0062]

[0063] Among them, C p It represents the law of change of rock mass elastic wave velocity with elevation, λ and μ are the law of change of rock mass elastic Lame constant with elevation, and ρ is the law of change of rock mass density with elevation.

[0064] The acoustic wave velocity and elastic wave velocity of each rock layer satisfy the relationship described in formula (3):

[0065]

[0066] Among them, Cp-so1 represents the distribution law of the first group of frequency sound wave velocity with elevation, Cp-so2 represents the distribution law of the second group of frequency sound wave velocity with elevation, ω1 and ω2 represent two groups of sound wave test frequencies, which can be determined according to the transducer parameters.

[0067] According to equations (1) to (3), the rock viscosity inversion equations are constructed:

[0068] When the two different acoustic wave test frequencies are 20kHz and 30kHz respectively, the rock viscosity inversion equations can be expressed as:

[0069]

[0070] Through the above technical solution, the embodiment of the present application can use the acoustic wave detection data obtained from on-site geological exploration to stratify the rock mass, and establish a set of inversion equations for the viscosity of each layer of rock mass. By solving the set of equations, the viscosity of each layer of rock mass is inverted. Compared with the existing technology that uses wave test data to reflect the viscosity characteristics of rock mass in a fixed range, it can be applicable to more complex rock mass compositions in actual engineering projects, reflect the viscosity of rock mass in the entire construction area, and provide scientific support for rock viscosity in subsequent construction and scientific research processes.

[0071] In step S103, the distribution pattern of rock mass viscosity along the elevation is determined based on the rock mass viscosity parameters to determine the comprehensive rock mass viscosity coefficient of the rock mass in the target range.

[0072] The target range can be understood as the specified range of the rock mass to be calculated, and the rock mass composite viscosity coefficient can be understood as a parameter that describes the combined viscosity and elastic properties of the rock mass under long-term loads or rheological conditions (such as crustal movement, slope creep, and tunnel convergence). It reflects the rock mass's ability to resist shear deformation or flow and is an important indicator in rock mass rheology.

[0073] Specifically, the embodiment of the present application can calculate the comprehensive viscosity coefficient of the rock mass by formula (5):

[0074]

[0075] Among them, η ps It represents the comprehensive longitudinal wave viscosity coefficient of the rock mass, and h1 and h2 represent the range of the rock mass to be calculated.

[0076] The embodiment of the present application can calculate the comprehensive viscosity coefficient of the rock mass in the corresponding range by determining the range of the rock mass to be calculated. It can not only obtain the rock mass viscosity coefficient before construction and reduce the workload of subsequent tests, but also adapt to the rock mass viscosity parameters required for construction in different parts.

[0077] In order to enable those skilled in the art to more clearly understand the rock viscosity inversion method based on acoustic wave detection of the present application, the rock viscosity inversion method based on acoustic wave detection of the present application is further explained below with a specific embodiment.

[0078] like Figure 2 As shown, Figure 2 The overlying rock mass of an underground project consists of three layers. The viscosity calculation process can be achieved through the following steps:

[0079] Step S1: Acoustic wave measured data is obtained under 20kHz and 30kHz transducer tests, and the acoustic wave velocities of the three rock layers are set to Cp-so1(h) and Cp-so2(h) in sequence.

[0080] Step S2: Obtain the static parameters of the three-layer rock mass according to the rock sample static test data, and calculate the elastic wave velocity thereof as Cp(h) according to formula (2).

[0081] Step S3: constructing a rock mass viscosity inversion equation group;

[0082]

[0083] Step S4: Solve the equations to obtain the longitudinal wave viscosity coefficient η of each layer of rock mass p (h), and the rock mass shear wave viscosity coefficient η is obtained according to formula (1): s (h).

[0084] Step S5: Set the blast source elevation and the surface measurement point elevation of the construction area to h1 and h2 respectively.

[0085] Step S6: Substitute the calculated viscosity coefficient distribution law into formula (5) to calculate the comprehensive viscosity coefficient of the rock mass under this blasting.

[0086] Optionally, in one embodiment of the present application, the rock viscosity inversion equations are determined according to the relationship between the acoustic wave velocity and the elastic wave velocity of the rock formation, and are in the form of:

[0087]

[0088] Among them, Cp-so1 represents the distribution law of the first group of frequency sound wave velocity with elevation, Cp-so2 represents the distribution law of the second group of frequency sound wave velocity with elevation, and ω1 and ω2 represent the two groups of sound wave test frequencies respectively.

[0089] During the actual implementation process, the embodiment of the present application can carry out acoustic wave detection on the acoustic wave logging drilled in the rock mass to be tested, obtain measured acoustic wave data, and then derive two sets of acoustic wave velocity elevation distribution laws based on the measured acoustic wave data.

[0090] The embodiment of the present application can reflect the relationship between the acoustic wave velocity elevation distribution law, the law of change of the longitudinal wave viscosity coefficient of the rock mass with elevation, the acoustic wave test frequency, and the law of change of the elastic wave velocity of the rock mass with elevation based on the relationship between the acoustic wave velocity and the elastic wave velocity of the rock mass in the target range, and thus can solve the unknown parameters when other parameters are known.

[0091] Optionally, in one embodiment of the present application, the elastic wave velocity is obtained from the static physical and mechanical parameters of the rock sample, wherein the calculation formula of the static physical and mechanical parameters of the rock sample is:

[0092]

[0093] Among them, C p It represents the law of change of rock mass elastic wave velocity with elevation, λ and μ are the law of change of rock mass elastic Lame constant with elevation, and ρ is the law of change of rock mass density with elevation.

[0094] In actual implementation, the embodiment of the present application can determine the variation pattern of the elastic wave velocity of the rock mass with elevation based on the variation pattern of the elastic Lame constant of the rock mass with elevation and the variation pattern of the rock mass density with elevation.

[0095] The embodiment of the present application can determine the law of change of the elastic wave velocity of the rock mass with elevation, and use the determined law of change of the elastic wave velocity of the rock mass with elevation as a known variable in combination with the rock mass viscosity inversion equation group to calculate the law of change of the longitudinal wave viscosity coefficient of the rock mass with elevation, so as to solve the comprehensive viscosity coefficient of the rock mass.

[0096] Optionally, in one embodiment of the present application, the calculation formula for the comprehensive viscosity coefficient of the rock mass in the target range is:

[0097]

[0098] Among them, η ps It represents the comprehensive longitudinal wave viscosity coefficient of the rock mass, and h1 and h2 represent the range of the rock mass to be calculated.

[0099] During actual implementation, the embodiment of the present application can calculate the comprehensive longitudinal wave viscosity coefficient of the rock mass through the calculation formula of the comprehensive viscosity coefficient of the rock mass according to the range of the rock mass to be calculated and the law of change of the longitudinal wave viscosity coefficient with elevation.

[0100] The embodiment of the present application can calculate the comprehensive viscosity coefficient of the rock mass at different parts, take into account the viscosity differences of the rock mass at different elevations, and well reflect the viscosity of each layer of rock mass. It can be applied to the viscosity inversion of various complex rock masses in actual engineering.

[0101] According to the rock viscosity inversion method based on acoustic wave detection proposed in the embodiment of the present application, the elevation distribution pattern of rock acoustic waves can be obtained through the acoustic wave detection data obtained from on-site geological exploration, so as to make full use of on-site engineering survey data without the need to carry out other additional tests; and then a rock viscosity inversion equation group is constructed, and the equation group is solved to realize the inversion of rock viscosity parameters, determine the distribution pattern of rock viscosity along the elevation, and finally determine the comprehensive viscosity coefficient of the rock according to the research scope, providing a scientific basis for rock viscosity in subsequent construction and scientific research. This solves the problem that the rock viscosity test method of the related technology needs to be combined with on-site production and reflects the viscosity characteristics of the rock in a fixed range, which leads to the problem that the viscosity data of the rock cannot be obtained before on-site construction and the viscosity of the rock in the entire construction area cannot be reflected.

[0102] Next, refer to the attached Figure 3 The present invention describes a rock viscosity inversion device based on acoustic wave detection proposed in an embodiment of the present application.

[0103] Figure 3 It is a block diagram of a rock viscosity inversion device based on acoustic wave detection according to an embodiment of the present application.

[0104] like Figure 3 As shown, the rock viscosity inversion device 10 based on acoustic wave detection includes: an acquisition module 100, a construction module 200, and an inversion module 300.

[0105] The acquisition module 100 is used to acquire acoustic wave detection data of engineering exploration, and obtain the elevation distribution law of rock mass acoustic waves based on the acoustic wave detection data.

[0106] The construction module 200 is used to construct a rock viscosity inversion equation group according to the elevation distribution law of the rock mass acoustic wave, and solve the rock viscosity inversion equation group to invert the rock mass viscosity parameters.

[0107] The inversion module 300 is used to determine the distribution pattern of rock mass viscosity along the elevation based on the rock mass viscosity parameters, so as to determine the comprehensive rock mass viscosity coefficient of the rock mass in the target range.

[0108] Optionally, in one embodiment of the present application, the rock viscosity inversion equations are determined according to the relationship between the acoustic wave velocity and the elastic wave velocity of the rock formation, and are in the form of:

[0109]

[0110] Among them, Cp-so1 represents the distribution law of the first group of frequency sound wave velocity with elevation, Cp-so2 represents the distribution law of the second group of frequency sound wave velocity with elevation, and ω1 and ω2 represent the two groups of sound wave test frequencies respectively.

[0111] Optionally, in one embodiment of the present application, the elastic wave velocity is obtained from the static physical and mechanical parameters of the rock sample, wherein the calculation formula of the static physical and mechanical parameters of the rock sample is:

[0112]

[0113] Among them, C p It represents the law of change of rock mass elastic wave velocity with elevation, λ and μ are the law of change of rock mass elastic Lame constant with elevation, and ρ is the law of change of rock mass density with elevation.

[0114] Optionally, in one embodiment of the present application, the calculation formula for the comprehensive viscosity coefficient of the rock mass in the target range is:

[0115]

[0116] Among them, η ps It represents the comprehensive longitudinal wave viscosity coefficient of the rock mass, and h1 and h2 represent the range of the rock mass to be calculated.

[0117] It should be noted that the above explanation of the embodiment of the rock viscosity inversion method based on acoustic wave detection is also applicable to the rock viscosity inversion device based on acoustic wave detection in this embodiment, and will not be repeated here.

[0118] According to the rock viscosity inversion device based on acoustic wave detection proposed in the embodiment of the present application, the elevation distribution law of rock sound waves can be obtained through the acoustic wave detection data obtained from on-site geological exploration, so as to make full use of on-site engineering survey data without the need to carry out other additional tests; and then a rock viscosity inversion equation group is constructed, and the equation group is solved to realize the inversion of rock viscosity parameters, determine the distribution law of rock viscosity along the elevation, and finally determine the comprehensive viscosity coefficient of the rock according to the research scope, providing a scientific basis for rock viscosity in subsequent construction and scientific research. This solves the problem that the rock viscosity test method of the related technology needs to be combined with on-site production and reflects the viscosity characteristics of the rock in a fixed range, which leads to the problem that the viscosity data of the rock cannot be obtained before on-site construction and the viscosity of the rock in the entire construction area cannot be reflected.

[0119] Figure 4 This is a schematic diagram of the structure of an electronic device provided in an embodiment of the present application. The electronic device may include:

[0120] Memory 401 , processor 402 , and computer programs stored in the memory 401 and executable on the processor 402 .

[0121] When the processor 402 executes the program, the rock viscosity inversion method based on acoustic wave detection provided in the above embodiment is implemented.

[0122] Furthermore, the electronic device further includes:

[0123] The communication interface 403 is used for communication between the memory 401 and the processor 402 .

[0124] The memory 401 is used to store computer programs that can be run on the processor 402 .

[0125] The memory 401 may include a high-speed RAM memory, and may also include a non-volatile memory (non-volatile memory), such as at least one disk memory.

[0126] If the memory 401, the processor 402, and the communication interface 403 are implemented independently, the communication interface 403, the memory 401, and the processor 402 can be connected to each other via a bus and communicate with each other. The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus, or an Extended Industry Standard Architecture (EISA) bus. The bus can be divided into an address bus, a data bus, a control bus, etc. For ease of representation, Figure 4 Only one thick line is used in the diagram, but this does not mean that there is only one bus or one type of bus.

[0127] Optionally, in a specific implementation, if the memory 401 , the processor 402 and the communication interface 403 are integrated on a chip, the memory 401 , the processor 402 and the communication interface 403 can communicate with each other through an internal interface.

[0128] The processor 402 may be a central processing unit (CPU), an application specific integrated circuit (ASIC), or one or more integrated circuits configured to implement the embodiments of the present application.

[0129] An embodiment of the present application further provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the above-mentioned rock viscosity inversion method based on acoustic wave detection.

[0130] An embodiment of the present application further provides a computer program product, which stores a computer program. When the program is executed by a processor, the above-mentioned rock viscosity inversion method based on acoustic wave detection is implemented.

[0131] In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present application. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described can be combined in any one or N embodiments or examples in a suitable manner. In addition, those skilled in the art can combine and combine different embodiments or examples described in this specification and features of different embodiments or examples without contradiction.

[0132] Furthermore, the terms "first" and "second" are used for descriptive purposes only and should not be understood to indicate or imply relative importance or implicitly specify the number of technical features indicated. Thus, a feature specified as "first" or "second" may explicitly or implicitly include at least one such feature. In the description of this application, "N" means at least two, for example, two, three, etc., unless otherwise specifically defined.

[0133] Any process or method description in a flowchart or otherwise described herein may be understood to represent a module, fragment or portion of code comprising one or N executable instructions for implementing a custom logical function or process step, and the scope of the preferred embodiments of the present application includes alternative implementations in which functions may be performed in a different order than shown or discussed, including performing functions in a substantially simultaneous manner or in a reverse order depending on the functions involved, which should be understood by those skilled in the art to which the embodiments of the present application pertain.

[0134] The logic and / or steps represented in the flowcharts or otherwise described herein, for example, can be considered as a sequenced list of executable instructions for implementing the logical functions, and can be embodied in any computer-readable medium for use by, or in conjunction with, an instruction execution system, apparatus, or device (e.g., a computer-based system, a system including a processor, or other system that can fetch and execute instructions from an instruction execution system, apparatus, or device). For purposes of this specification, a "computer-readable medium" can be any device that can contain, store, communicate, propagate, or transport a program for use by, or in conjunction with, an instruction execution system, apparatus, or device. More specific examples (a non-exhaustive list) of computer-readable media include the following: an electrical connection with one or N wires (electronic devices), a portable computer disk cartridge (magnetic device), random access memory (RAM), read-only memory (ROM), erasable and programmable read-only memory (EPROM or flash memory), fiber optic devices, and a portable compact disc read-only memory (CDROM). In addition, the computer-readable medium may even be paper or other suitable medium on which the program is printed, since the program can be obtained electronically by optically scanning the paper or other medium and then editing, interpreting or processing it in other suitable ways as necessary, and then storing it in a computer memory.

[0135] It should be understood that various parts of the present application can be implemented using hardware, software, firmware, or a combination thereof. In the above embodiment, the N steps or methods can be implemented using software or firmware stored in a memory and executed by a suitable instruction execution system. If implemented using hardware, as in another embodiment, it can be implemented using any one or a combination of the following technologies known in the art: a discrete logic circuit having a logic gate circuit for implementing a logic function on a data signal, an application-specific integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

[0136] Those skilled in the art will understand that all or part of the steps in the method of the above embodiment can be completed by instructing related hardware through a program, and the program can be stored in a computer-readable storage medium. When the program is executed, it includes one or a combination of the steps of the method embodiment.

[0137] In addition, the functional units in the various embodiments of the present application may be integrated into a processing module, or each unit may exist physically separately, or two or more units may be integrated into a module. The above-mentioned integrated module may be implemented in the form of hardware or in the form of a software functional module. If the integrated module is implemented in the form of a software functional module and sold or used as an independent product, it may also be stored in a computer-readable storage medium.

[0138] The storage medium mentioned above may be a read-only memory, a magnetic disk, or an optical disk, etc. Although the embodiments of the present application have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting the present application. Persons skilled in the art may make changes, modifications, substitutions, and variations to the above embodiments within the scope of the present application.

Claims

1. A rock viscosity inversion method based on acoustic wave detection, characterized in that: The following steps are involved: Acquiring acoustic wave detection data from engineering exploration, and obtaining the elevation distribution law of rock mass acoustic waves based on the acoustic wave detection data; Constructing a rock mass viscosity inversion equation group according to the elevation distribution law of the rock mass acoustic wave, and solving the rock mass viscosity inversion equation group to invert the rock mass viscosity parameters; The distribution law of rock mass viscosity along the elevation is determined based on the rock mass viscosity parameters to determine the comprehensive rock mass viscosity coefficient of the rock mass in the target range.

2. The method according to claim 1, characterized in that The rock viscosity inversion equations are determined based on the relationship between the acoustic wave velocity and the elastic wave velocity of the rock layer, and are in the form of: Among them, Cp-so1 represents the distribution law of the first group of frequency sound wave velocity with elevation, Cp-so2 represents the distribution law of the second group of frequency sound wave velocity with elevation, and ω1 and ω2 represent the two groups of sound wave test frequencies respectively.

3. The method according to claim 2, characterized in that The elastic wave velocity is obtained from the static physical and mechanical parameters of the rock sample, wherein the calculation formula of the static physical and mechanical parameters of the rock sample is: Among them, C p It represents the law of change of rock mass elastic wave velocity with elevation, λ and μ are the law of change of rock mass elastic Lame constant with elevation, and ρ is the law of change of rock mass density with elevation.

4. The method according to any one of claims 1 to 3, characterized in that The calculation formula of the comprehensive viscosity coefficient of the rock mass in the target range is: Among them, η ps It represents the comprehensive longitudinal wave viscosity coefficient of rock mass, h1 and h2 represent the range of rock mass to be calculated.

5. A rock viscosity inversion device based on acoustic wave detection, characterized in that: include: An acquisition module is used to acquire acoustic wave detection data of engineering exploration and obtain the elevation distribution law of rock mass acoustic waves based on the acoustic wave detection data; A construction module is used to construct a rock viscosity inversion equation group according to the elevation distribution law of the rock mass acoustic wave, and solve the rock viscosity inversion equation group to invert the rock mass viscosity parameters; The inversion module is used to determine the distribution law of rock viscosity along the elevation based on the rock viscosity parameters to determine the comprehensive viscosity coefficient of the rock mass in the target range.

6. The device according to claim 5, characterized in that The rock viscosity inversion equations are determined based on the relationship between the acoustic wave velocity and the elastic wave velocity of the rock layer, and are in the form of: Among them, Cp-so1 represents the distribution law of the first group of frequency sound wave velocity with elevation, Cp-so2 represents the distribution law of the second group of frequency sound wave velocity with elevation, and ω1 and ω2 represent the two groups of sound wave test frequencies respectively.

7. The device according to claim 5, characterized in that The elastic wave velocity is obtained from the static physical and mechanical parameters of the rock sample, wherein the calculation formula of the static physical and mechanical parameters of the rock sample is: Among them, C p It represents the law of change of rock mass elastic wave velocity with elevation, λ and μ are the law of change of rock mass elastic Lame constant with elevation, and ρ is the law of change of rock mass density with elevation.

8. An electronic device, characterized in that: include: 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 rock viscosity inversion method based on acoustic wave detection according to any one of claims 1 to 4.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: The program is executed by a processor to implement the rock viscosity inversion method based on acoustic wave detection as described in any one of claims 1 to 4.

10. A computer program product comprising a computer program, characterized in that The computer program is executed to implement the rock viscosity inversion method based on acoustic wave detection as described in any one of claims 1 to 4.

Images