Regional rock mass quality assessment method and related components considering unloading fracture effect

Abstract

The present invention discloses a method for assessing the quality of regional rock mass considering the unloading fracture effect and related components. The method includes obtaining all sample joint surface data in the sample area; based on all sample joint surface data, a discrete fracture network model; obtaining microseismic data in the monitoring area, and calculating the geometric parameters of the new fracture surface based on the microseismic data; importing the geometric parameters of the new fracture surface into the discrete fracture network model to obtain a three-dimensional discrete fracture network calculation model; performing a simulated drilling operation in the three-dimensional discrete fracture network calculation model to obtain the intersection of each discrete fracture surface and the simulated drilling, and calculating the target RQD value based on the intersection. The three-dimensional discrete fracture network calculation model obtained by the method constructs the rock mass in the sample area with a three-dimensional model, and then dynamically simulates the state of the rock mass, so as to more quickly and effectively evaluate the quality of the rock mass, that is, the target RQD value can be quickly and accurately obtained, and the accuracy of the rock mass quality evaluation can be improved.

Description

Technical Field

[0001] The present invention relates to the field of data modeling, and in particular to a regional rock mass quality assessment method considering unloading fracture effect and related components. Background Art

[0002] At present, with the urgent demand for infrastructure in the development of the national economy, the number of underground engineering construction is gradually increasing, which often inevitably requires the excavation of rock and soil. During the excavation process, rock unloading has a great influence on the stability of the rock mass of underground engineering. Therefore, it is of great significance to evaluate the quality of rock mass in the underground engineering area under the fracture effect of excavation unloading. Among them, accurately grasping the geometric characteristics of the rock mass structural surface is crucial for rock mass quality evaluation and stability analysis.

[0003] Under natural conditions, the structural surface parameters of rock mass can be measured by traditional manual on-site contact measurement methods. With tape measure and compass as the main tools, the structural surface information is manually measured one by one, investigated and collected, such as line measurement method, window measurement method, etc. At the same time, during the excavation process, new fracture surfaces induced by unloading will be generated in the rock mass, and the parameters of these fracture surfaces are difficult to obtain using traditional line measurement methods. Microseismic monitoring is a non-destructive real-time passive monitoring technology that studies rock fracture. There are countless rock fracture processes during rock excavation and construction. Therefore, microseismic monitoring technology can monitor the entire process of structural surface incubation in rock under the effect of unloading fracture. Through real-time microseismic monitoring, accurate and effective parameters of new fracture surfaces can be grasped in time, which is conducive to further rock mass quality assessment. For a long time, people have mainly relied on limited drilling holes and limited rock exposure surfaces to infer rock mass quality, accumulated a lot of experience and achieved rich research results. Among them, the most representative descriptions of rock mass quality are the RMR rock mass engineering classification system and the Q classification system. As an important evaluation parameter in these classification systems, the rock quality index RQD has been unanimously recognized by the majority of engineers and is widely used in rock engineering. At the same time, although the traditional drilling line measurement method is more cumbersome, it can still accurately and effectively measure and calculate the RQD value of the rock mass in its natural state.

[0004] However, due to the unloading fracture effect of the excavated rock mass, new fracture surfaces will be generated inside it, and the traditional drilling line measurement method cannot accurately and effectively calculate the RQD value of the rock mass. Summary of the invention

[0005] The purpose of the present invention is to provide a regional rock mass quality assessment method, device and related components taking into account the unloading fracture effect, aiming to solve the problem that the existing drilling line measurement method cannot accurately and effectively calculate the RQD value of the rock mass.

[0006] In order to solve the above technical problems, the object of the present invention is achieved through the following technical solutions: providing a method for assessing the quality of regional rock mass taking into account the unloading fracture effect, which comprises:

[0007] Based on the selected sample area, obtaining sample joint surface data corresponding to all natural rock mass structural surfaces in the sample area;

[0008] Based on all sample joint surface data, a discrete fracture network model is established to characterize the spatial distribution of joints.

[0009] Based on the preset ESG microseismic monitoring system, the microseismic data of the monitoring area is obtained, and the geometric parameters of the new fracture surface are calculated based on the microseismic data;

[0010] Importing the geometric parameters of the newly formed fracture surface into the discrete fracture network model to obtain a three-dimensional discrete fracture network calculation model;

[0011] A simulated drilling operation is performed in the three-dimensional discrete fracture network calculation model to obtain the intersection points between each discrete fracture surface and the simulated drilling hole, and the target RQD value is calculated based on the intersection points.

[0012] In addition, the technical problem to be solved by the present invention is to provide a regional rock mass quality assessment device taking into account the unloading fracture effect, which comprises:

[0013] A sample joint surface data acquisition unit, used to acquire sample joint surface data corresponding to all natural rock mass structural surfaces in the sample area based on the selected sample area;

[0014] A discrete fracture network model building unit is used to build a discrete fracture network model that characterizes the spatial distribution of joints based on all sample joint surface data;

[0015] A geometric parameter calculation unit of a newly formed fracture surface, used to obtain microseismic data of a monitoring area based on a preset ESG microseismic monitoring system, and calculate geometric parameters of a newly formed fracture surface based on the microseismic data;

[0016] A three-dimensional discrete fracture network calculation model establishment unit, used for importing the geometric parameters of the new fracture surface into the discrete fracture network model to obtain a three-dimensional discrete fracture network calculation model;

[0017] The RQD value calculation unit is used to perform a simulated drilling operation in the three-dimensional discrete fracture network calculation model to obtain the intersection points of each discrete fracture surface and the simulated drilling, and calculate the target RQD value based on the intersection points.

[0018] In addition, an embodiment of the present invention provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the method for assessing regional rock quality considering the unloading fracture effect described in the first aspect above is implemented.

[0019] In addition, an embodiment of the present invention further provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the processor executes the regional rock quality assessment method considering the unloading fracture effect as described in the first aspect above.

[0020] The embodiment of the present invention discloses a method, device and related components for assessing the quality of regional rock mass considering the unloading fracture effect, wherein the method comprises: based on the selected sample area, obtaining the sample joint surface data corresponding to all the natural rock mass structural surfaces in the sample area; based on all the sample joint surface data, establishing a discrete fracture network model that characterizes the spatial distribution of joints; based on the preset ESG microseismic monitoring system, obtaining the microseismic data of the monitoring area, and calculating the geometric parameters of the new fracture surface based on the microseismic data; importing the geometric parameters of the new fracture surface into the discrete fracture network model to obtain a three-dimensional discrete fracture network calculation model; performing a simulated drilling operation in the three-dimensional discrete fracture network calculation model to obtain the intersection of each discrete fracture surface and the simulated drilling, and calculating the target RQD value based on the intersection. The three-dimensional discrete fracture network calculation model obtained by the method is to construct the rock mass in the sample area with a three-dimensional model, and then dynamically simulate the state of the rock mass, so as to more quickly and effectively assess the quality of the rock mass, that is, to quickly and accurately obtain the target RQD value, and improve the accuracy of the rock mass quality evaluation. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0022] Figure 1 A schematic flow chart of a method for assessing regional rock mass quality taking into account unloading fracture effects provided by an embodiment of the present invention;

[0023] Figure 2 A joint surface pole diagram of a regional rock mass quality assessment method considering unloading fracture effect provided by an embodiment of the present invention;

[0024] Figure 3A joint pole isodensity map of a regional rock mass quality assessment method considering unloading fracture effect provided by an embodiment of the present invention;

[0025] Figure 4 A joint parameter histogram and a fitting curve of a regional rock mass quality assessment method considering unloading fracture effect provided by an embodiment of the present invention;

[0026] Figure 5 A three-dimensional discrete fracture network calculation model for a regional rock mass quality assessment method considering unloading fracture effects provided by an embodiment of the present invention;

[0027] Figure 6 A schematic diagram of the structure of a simulated borehole in a three-dimensional discrete fracture network calculation model of a regional rock mass quality assessment method considering unloading fracture effects provided by an embodiment of the present invention;

[0028] Figure 7 A two-dimensional cross-sectional diagram of a simulated borehole in a three-dimensional discrete fracture network calculation model of a regional rock mass quality assessment method considering unloading fracture effects provided by an embodiment of the present invention;

[0029] Figure 8 A schematic block diagram of a regional rock mass quality assessment device considering unloading fracture effect provided by an embodiment of the present invention;

[0030] Fig. 9 A schematic block diagram of a computer device provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0031] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments 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.

[0032] It should be understood that when used in this specification and the appended claims, the terms "include" and "comprises" indicate the presence of described features, integers, steps, operations, elements and / or components, but do not exclude the presence or addition of one or more other features, integers, steps, operations, elements, components and / or combinations thereof.

[0033] It should also be understood that the terms used in this specification of the present invention are only for the purpose of describing specific embodiments and are not intended to limit the present invention. As used in the specification of the present invention and the appended claims, unless the context clearly indicates otherwise, the singular forms "a", "an" and "the" are intended to include plural forms.

[0034] It should be further understood that the term "and / or" used in the present description and the appended claims refers to and includes any and all possible combinations of one or more of the associated listed items.

[0035] See also Figure 1 , Figure 1 A schematic flow chart of a method for assessing regional rock mass quality taking into account unloading fracture effects provided by an embodiment of the present invention;

[0036] like Figure 1 As shown, the method includes steps S101 to S105.

[0037] S101, based on the selected sample area, obtaining sample joint surface data corresponding to all natural rock mass structural surfaces in the sample area;

[0038] S102, based on all sample joint surface data, establishing a discrete fracture network model that characterizes the spatial distribution of joints;

[0039] S103, based on a preset ESG microseismic monitoring system, obtaining microseismic data of the monitoring area, and calculating geometric parameters of the newly formed rupture surface based on the microseismic data;

[0040] S104, importing the geometric parameters of the newly formed fracture surface into the discrete fracture network model to obtain a three-dimensional discrete fracture network calculation model;

[0041] S105, performing a simulated drilling operation in the three-dimensional discrete fracture network calculation model to obtain the intersection points between each discrete fracture surface and the simulated drilling, and calculating the target RQD value based on the intersection points.

[0042] In this embodiment, through manual on-site contact measurement, with a tape measure and a compass as the main tools, the sample joint surface data corresponding to the natural rock structural surface of the sample area is investigated and collected. It should be noted that the natural rock in each sample area has multiple joint surfaces, and the sample joint surface data of each joint surface needs to be collected. Then, based on the collected sample joint surface data, a discrete fracture network model characterizing the spatial distribution of joints is established.

[0043] It should be noted that the ESG microseismic monitoring system of the present application includes an ESG microseismic monitoring data acquisition box, a host (usually a computer) and multiple sensors. During installation, the local rock mass of the underground project is used as the monitoring area (i.e., the sample area of ​​the present application), and the sensors equipped with the ESG microseismic monitoring system are installed in the rock mass area whose quality is to be evaluated. At least 6 sensors are installed, and the sensors are arranged in a manner that enables the sensors to form a network structure in space and cover the monitoring area. Each sensor is connected to the ESG microseismic monitoring data acquisition box, and then the data acquisition box is connected to the host for processing data signals. The excavation process of the underground project will induce microseismic events in the monitoring area. During the excavation process of the underground project, the monitoring area is monitored by the ESG microseismic monitoring system, and the microseismic data of the microseismic events generated in the monitoring area are measured. The microseismic data include the source location and the time of occurrence of the microseism. It should be noted that excavation unloading will induce the appearance of new fracture surfaces in the rock.

[0044] After obtaining the microseismic data corresponding to the sample area, the geometric parameters of the corresponding new fracture surface are calculated, and then the geometric parameters of the new fracture surface are input into the discrete fracture network model to obtain a three-dimensional discrete fracture network calculation model; finally, the three-dimensional discrete fracture network calculation model is used to simulate the drilling operation to obtain the intersection points of each discrete fracture surface and the simulated drilling, and the target RQD value is calculated based on the intersection points.

[0045] The three-dimensional discrete fracture network calculation model obtained in the present application constructs the rock mass in the sample area using a three-dimensional model, and then dynamically simulates the state of the rock mass, so as to more quickly and effectively evaluate the rock mass quality, that is, the target RQD value can be quickly and accurately obtained, thereby improving the accuracy of rock mass quality evaluation.

[0046] Microseismic monitoring technology is used to capture the geometric parameters of new rock fracture surfaces induced by excavation unloading, and the rock mass in a specific area is constructed into a three-dimensional model through 3DEC software to dynamically simulate the rock mass state, so as to facilitate a faster and more effective assessment of rock mass quality.

[0047] In a specific embodiment, step S101 includes the following steps:

[0048] S10, delineating a sample area, and obtaining sample parameters of the sample area, wherein the sample parameter data includes a half-trace length parameter, a spacing parameter, a strike parameter, a dip parameter, and a dip angle parameter;

[0049] S11, converting the sample parameters into a corresponding sample format, and drawing a corresponding joint surface pole diagram and a joint pole isodensity diagram based on the converted sample data;

[0050] S12, based on a preset region selection rule, selecting and grouping all corresponding regions in the joint surface pole map and the joint pole equal density map to obtain a corresponding target region;

[0051] S13, calculating the disk diameter parameter and the joint surface quantity parameter of the disk corresponding to each target area.

[0052] In this embodiment, it should be noted that the strike parameter, dip parameter and dip angle parameter are also called the occurrence data of the joint surface. When actually obtaining the sample parameters of the sample area, the strike parameter, dip parameter and dip angle parameter need to be input into an Excel table in the order of strike, dip and dip angle. Then, the strike parameter, dip parameter and dip angle parameter are converted into the dip angle / dip angle format required by the Dips (Development Information Processing System) data conversion program using the Dips (Development Information Processing System) data conversion program. After the format conversion, the corresponding joint surface pole diagram (such as Figure 2 ) and joint pole isodensity maps (as shown in Figure 3 As shown in the figure, the joint surface pole diagram and the joint pole isodensity diagram can be used to intuitively determine the areas where the occurrence distribution is relatively concentrated. These areas are the target areas of this application. Figure 2 and Figure 3 As shown, there are 3 groups of target area settings.

[0053] It is also necessary to conduct statistical analysis on the dip parameters, inclination parameters, spacing values ​​and half-trace length parameters of all joint surfaces in each target area, and draw a histogram and fitting curve, as shown in the attached figure. Figure 4 As shown in Table 1, the probability distribution parameters obeyed by each joint parameter (strike parameter, dip parameter and inclination parameter) are finally determined according to the fitting curve, as shown in Table 1 below:

[0054] Table 1:

[0055]

[0056]

[0057] In a specific embodiment, the Bqecher model is used to assume that the joint surface is a disk model, the joint trace length of the joint surface is regarded as the chord length of the disk model, and the number of joints is controlled by the joint density, wherein the number of joints is obtained by step S24, and the step S13 includes the following steps:

[0058] S21, obtaining all half-track length parameters corresponding to all joint surfaces in each target area, and based on the half-track length parameters, calculating the mean value a of the half-track length parameters in each target area according to the following formula: h :

[0059]

[0060] Among them, a g is the overall mean of the total trace length of all joint surfaces, a is the half trace length parameter of the joint surface, μ is the probability density of the trace length of the joint surface, and F(a) is the cumulative probability distribution of the trace length of the joint surface;

[0061] In this embodiment, the cumulative probability distribution of the trace length of the joint surface is obtained by fitting a probability density function curve by trace length data, and then the probability distribution is obtained by integrating the probability density function;

[0062] S22, calculate the average length of the whole track in each target area according to the following formula

[0063]

[0064] Wherein, b is the radius of the disk model;

[0065] It should be noted that the mean value of the full trace length is the average chord length of the disk, because in the Bqecher model, the trace is often regarded as the chord length of the disk. The mean value of the full trace length obtained in step S21 is equal to twice the mean value of the half trace length. The relationship between the average chord length of the disk and the radius is obtained in step S22. Since the mean value of the full trace length is the average chord length of the disk, the relationship between the radius and the mean value of the full trace length can be obtained in this way, that is, there is a positive relationship between the mean value of the full trace length and the mean value of the diameter. relationship.

[0066] S23, calculate the volume density (c v ) j :

[0067]

[0068]

[0069] Among them, (c i ) j is the linear density of the jth group of joints, E(d 2 ) is the second moment of the disk diameter distribution; e i is the average vector direction of the i-th group of joints; e i is the number of joint surfaces in the i-th group;

[0070] It should be noted that the second moment of the disk diameter distribution is obtained by the weighted average of the squares of the joint surface diameters.

[0071] S24. Obtain the number of joint surfaces n according to the following formula:

[0072] n=Vc v

[0073] Where V represents the simulated space volume.

[0074] It should be noted that the above calculation process is to calculate the disk diameter parameter based on the half-trace length data measured on site, and to calculate the joint surface quantity parameter based on the spacing parameters measured on site.

[0075] In a specific embodiment, step S102 includes the following steps:

[0076] S30. Input the dip distribution parameter, dip angle distribution parameter, disk diameter distribution parameter and the number of joint surfaces into 3DEC software in sequence to establish a discrete fracture network model that characterizes the spatial distribution of joints, wherein the strike parameter, dip parameter and dip angle parameter adopt a normal distribution, the disk diameter adopts a log-normal distribution, and the disk diameter is set based on a preset diameter range.

[0077] In this embodiment, the minimum diameter of the disk model is 0.1m and the maximum diameter is 1m. It should be noted that the center distribution of the joint disk mainly includes uniform distribution, Gaussian distribution and other custom distribution forms. In order to simplify the model and increase calculation efficiency, the joint position is regarded as a uniform distribution here.

[0078] In a specific embodiment, step S103 includes the following steps:

[0079] S40, obtaining a point source corresponding to the earthquake source based on a preset point source setting rule;

[0080] It should be noted that when the observation distance from the earthquake source to the sensor and the wavelength of the seismic wave are much larger than the earthquake source rupture scale, the earthquake source can be assumed to be a point source, the observation distance is far field compared to the wavelength, and the rock mass between the sensor and the earthquake source is assumed to be a continuous, uniform and isotropic infinite space medium.

[0081] S41, the long-range P-wave displacement field s of the point source as described by the following formula:

[0082] s=F′T=kQT

[0083] in, ρ represents the density of the rock mass, v represents the P-wave velocity, R represents the observation distance, and Q represents the excitation matrix;

[0084] It should be noted that the long-range P-wave displacement field (s) can be expressed as the product of the spatial derivative of the Green’s function (F) and the second-order moment tensor (T);

[0085] S42, wherein the excitation matrix Q is calculated as follows:

[0086] Q=α o α u

[0087] Where α is the direction cosine vector between the point source and the sensor;

[0088] It should be noted that o and u can be 1, 2, or 3, respectively. When o and u are 1, they represent the x-axis direction of the coordinate system where the rock mass is located; when o and u are 2, they represent the y-axis direction of the coordinate system where the rock mass is located; and when o and u are 3, they represent the z-axis direction of the coordinate system where the rock mass is located.

[0089] S43. Calculate the microseismic moment tensor of the initial P wave motion as follows:

[0090]

[0091] Among them, R n represents the observation distance from the point source to the nth sensor, represents the direction cosine vector between the point source and the nth sensor; T represents the moment tensor;

[0092] In this embodiment, R represents the observation distance from the source to the sensor. n represents the observation distance to the nth sensor, This is the calculation of the excitation matrix Q in the above formula. The superscript represents the sensor number, the subscripts 1, 2, and 3 represent the x-axis direction, y-axis direction, and z-axis direction of the coordinate system of the region, respectively, and a is the direction cosine vector between the point source (rupture source) and the sensor.

[0093] S44. Calculate the eigenvalue of the moment tensor T as follows:

[0094]

[0095] Where s represents the displacement of the fracture surface in the direction of movement, A represents the surface area of ​​the new fracture surface, and μ represents the Lame constant. represents the movement direction of the new fracture surface, Indicates the normal direction of the new fracture surface;

[0096] S45. Calculate the eigenvector corresponding to the eigenvalue of the moment tensor according to the following formula:

[0097]

[0098] Vector and The angle between and The angle and and The angle between them is β / 2, which can be calculated as follows:

[0099]

[0100] S46. Calculate the normal direction of the fracture surface according to the following formula and direction of movement

[0101]

[0102]

[0103] S47. According to the spatial position relationship of the newly formed fracture surface, the fracture surface occurrence data of the newly formed fracture surface geometric parameters are calculated as follows:

[0104]

[0105] Where: θ is the strike of the new fracture surface, γ is the inclination angle of the new fracture surface;

[0106] S48. Calculate the radius of the fracture surface as follows:

[0107]

[0108] Where: z s represents the S wave velocity; p s Indicates the S-wave angular frequency.

[0109] In a specific embodiment, step S104 includes the following steps:

[0110] S50, inputting the geometric dimensions of the rock mass into the 3DEC (3 Dimension Distinct Element Code) software;

[0111] S51, inputting the radius parameter of the new fracture surface circle into the 3DEC software;

[0112] S52, inputting the strike parameters and dip parameters of the new fracture surface into the 3DEC software;

[0113] S53, inputting the earthquake source coordinate parameters and the position of the new fracture surface circle into the 3DEC software to obtain a three-dimensional discrete fracture network calculation model.

[0114] In a specific embodiment, step S105 includes the following steps:

[0115] S60, performing a simulated drilling operation in the three-dimensional discrete fracture network calculation model to obtain the intersection points of each discrete fracture surface and the simulated drilling, and calculating the target RQD value based on the intersection points, including:

[0116] S61, using the geometry function of 3DEC to generate numerical simulation drilling and obtain all corresponding target holes;

[0117] S62, calculating the intersection points of each target hole and the joint surface to obtain a corresponding intersection point set;

[0118] S63, calculating the spacing value between adjacent intersection points in the intersection point set;

[0119] S64, traversing and determining whether each spacing value is greater than a preset spacing threshold, if the current spacing value is greater than the preset spacing threshold, taking the current spacing value as a target spacing value;

[0120] It should be noted that the spacing threshold of this application is set to 10cm;

[0121] S65. After all target spacing values ​​are obtained, all target spacing values ​​are added together to obtain a total length;

[0122] S66. Calculate the sub-target RQD value as follows:

[0123]

[0124] Among them, L i represents the i-th spacing value, and L represents the total number of spacing values ​​in the intersection set;

[0125] S67, after obtaining the sub-target RQD values ​​corresponding to all the target holes, averaging all the sub-target RQD values ​​to obtain the target RQD value of the rock mass

[0126] In this embodiment, 3DEC software is used to establish a three-dimensional discrete fracture network calculation model of the rock mass in the local area of ​​the underground project including the excavation, and it is ensured that the calculation model and the actual underground area rock mass are in the same coordinate system, such as Figure 5 As shown in the figure, the calculation model is simulated by the parameters of the structural surface. The structural surface is composed of a disk model. The horizontal width × vertical height × longitudinal length of the model is 5m × 5m × 5m. The geometry function of 3DEC is used to generate numerical simulation boreholes. Several boreholes that penetrate the rock mass are formed on the assumed three-dimensional model. All boreholes are arranged parallel to each other in the rock mass, as shown below Figure 6 and Figure 7 As shown, find the intersections of each discrete fracture surface and the simulated borehole, the number of intersections is N, and count the total length L of the line segments with an intersection spacing greater than or equal to 10 cm. This embodiment makes three groups of simulated boreholes, and their basic parameters are as shown in Tables 3 and 4:

[0127] Table 3:

[0128] Drilling 1 Drilling 2 Drilling 3 Intersection 1 -14 -80 -186 Intersection 2 29 -77 -55 Intersection 3 38 -23 -12 Intersection 4 106 20 -9 Intersection 5 222 85 Intersection 6 124

[0129] Table 4:

[0130]

[0131] According to the calculation formula of RQD: The RQD value of each simulated drill hole was calculated.

[0132] The calculation example is as follows:

[0133] Drilling 1:

[0134] Drill 2:

[0135] Drilling 3:

[0136] Take the average value of the RQD values ​​of several boreholes to get the average RQD value of the rock mass

[0137]

[0138] It should be noted that the structural surface that appears in the present application is the joint surface or the fracture surface, and the new fracture surface is the new fracture surface generated by vibration.

[0139] The embodiment of the present invention also provides a regional rock mass quality assessment device considering the unloading fracture effect, and the regional rock mass quality assessment device considering the unloading fracture effect is used to execute any embodiment of the regional rock mass quality assessment method considering the unloading fracture effect. Figure 8 , Figure 8 It is a schematic block diagram of a regional rock mass quality assessment device taking into account unloading fracture effects provided by an embodiment of the present invention.

[0140] like Figure 8 As shown, a regional rock mass quality assessment device 500 considering unloading fracture effect comprises:

[0141] The sample joint surface data acquisition unit 501 is used to acquire the sample joint surface data corresponding to all natural rock mass structural surfaces in the sample area based on the selected sample area;

[0142] A discrete fracture network model building unit 502 is used to build a discrete fracture network model that characterizes the spatial distribution of joints based on all sample joint surface data;

[0143] A geometric parameter calculation unit 503 of a newly formed fracture surface is used to obtain microseismic data of a monitoring area based on a preset ESG microseismic monitoring system, and calculate geometric parameters of a newly formed fracture surface based on the microseismic data;

[0144] A three-dimensional discrete fracture network calculation model establishing unit 504 is used to import the geometric parameters of the new fracture surface into the discrete fracture network model to obtain a three-dimensional discrete fracture network calculation model;

[0145] The RQD value calculation unit 505 is used to perform a simulated drilling operation in the three-dimensional discrete fracture network calculation model to obtain the intersection points between each discrete fracture surface and the simulated drilling, and calculate the target RQD value based on the intersection points.

[0146] In a specific embodiment, the sample joint surface data acquisition unit includes:

[0147] A sample area delineation unit is used to delineate a sample area and obtain sample parameters of the sample area, wherein the sample parameter data includes a half-trace length parameter, a spacing parameter, a strike parameter, a dip parameter and a dip angle parameter;

[0148] A drawing unit, used for converting the sample parameters into a corresponding sample format, and drawing a corresponding joint surface pole diagram and a joint pole isodensity diagram based on the converted sample data;

[0149] Selecting a target area unit, for selecting and grouping all corresponding areas in the joint surface pole map and the joint pole equal density map based on a preset area selection rule, to obtain a corresponding target area;

[0150] A parameter calculation unit, used for calculating a disk diameter parameter and a joint surface quantity parameter of a disk corresponding to each target area;

[0151] Wherein, the parameter calculation unit includes:

[0152] A half-trace length parameter mean value calculation unit is used to obtain all half-trace length parameters corresponding to all joint surfaces in each target area, and based on the half-trace length parameters, calculate the mean value a of the half-trace length parameters in each target area according to the following formula: h :

[0153]

[0154] Among them, a g is the overall mean of the total trace length of all joint surfaces, a is the half trace length parameter of the joint surface, μ is the probability density of the trace length of the joint surface, and F(a) is the cumulative probability distribution of the trace length of the joint surface;

[0155] The full trace length average unit is used to calculate the full trace length average in each target area according to the following formula

[0156]

[0157] Wherein, b is the radius of the disk model;

[0158] The volume density unit of the joint surface is used to calculate the volume density (c v ) j :

[0159]

[0160]

[0161] Among them, (c i ) h is the linear density of the jth group of joints, E(d 2 ) is the second moment of the disk diameter distribution; e i is the average vector direction of the i-th group of joints; e i is the number of joint surfaces in the i-th group;

[0162] The joint surface number calculation unit is used to obtain the joint surface number n according to the following formula:

[0163] n=Vc v

[0164] Where V represents the simulated space volume.

[0165] The three-dimensional discrete fracture network calculation model obtained by the device constructs the rock mass in the sample area using a three-dimensional model, and then dynamically simulates the state of the rock mass to facilitate a faster and more effective assessment of the rock mass quality, that is, the target RQD value can be quickly and accurately obtained, thereby improving the accuracy of rock mass quality evaluation.

[0166] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the above-described devices and units can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0167] The above-mentioned regional rock mass quality assessment device considering the unloading fracture effect can be realized in the form of a computer program. The computer program can be used in Fig. 9 Runs on the computer device shown.

[0168] See also Fig. 9 , Fig. 9 11 is a schematic block diagram of a computer device provided by an embodiment of the present invention. The computer device 1100 is a server, which may be an independent server or a server cluster composed of multiple servers.

[0169] See also Fig. 9 The computer device 1100 includes a processor 1102 , a memory, and a network interface 1105 connected via a system bus 1101 , wherein the memory may include a non-volatile storage medium 1103 and an internal memory 1104 .

[0170] The non-volatile storage medium 1103 may store an operating system 11031 and a computer program 11032. When the computer program 11032 is executed, the processor 1102 may execute a regional rock mass quality assessment method considering unloading fracture effects.

[0171] The processor 1102 is used to provide computing and control capabilities to support the operation of the entire computer device 1100 .

[0172] The internal memory 1104 provides an environment for the operation of the computer program 11032 in the non-volatile storage medium 1103. When the computer program 11032 is executed by the processor 1102, the processor 1102 can execute a regional rock mass quality assessment method that takes into account the unloading fracture effect.

[0173] The network interface 1105 is used for network communication, such as providing data information transmission, etc. Those skilled in the art will understand that Fig. 9 The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present invention, and does not constitute a limitation on the computer device 1100 to which the solution of the present invention is applied. The specific computer device 1100 may include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components.

[0174] Those skilled in the art will understand that Fig. 9 The embodiments of the computer device shown in the figure do not constitute a limitation on the specific composition of the computer device. In other embodiments, the computer device may include more or fewer components than shown in the figure, or combine certain components, or arrange the components differently. For example, in some embodiments, the computer device may only include a memory and a processor. In such embodiments, the structure and function of the memory and the processor are the same as those of the embodiment of the present invention. Fig. 9 The embodiments shown are consistent and will not be described again here.

[0175] It should be understood that in the embodiment of the present invention, the processor 1102 may be a central processing unit (CPU), and the processor 1102 may also be other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. Among them, the general-purpose processor may be a microprocessor or the processor may also be any conventional processor, etc.

[0176] In another embodiment of the present invention, a computer-readable storage medium is provided. The computer-readable storage medium may be a non-volatile computer-readable storage medium. The computer-readable storage medium stores a computer program, wherein when the computer program is executed by a processor, a method for assessing regional rock mass quality considering unloading fracture effect according to an embodiment of the present invention is implemented.

[0177] The storage medium is a physical, non-transient storage medium, for example, it can be a USB flash drive, a mobile hard disk, a read-only memory (ROM), a magnetic disk or an optical disk, etc., which can store program codes.

[0178] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the above-described equipment, devices and units can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0179] The above is only a specific embodiment of the present invention, but the protection scope of the present invention is not limited thereto. Any technician familiar with the technical field can easily think of various equivalent modifications or replacements within the technical scope disclosed by the present invention, and these modifications or replacements should be included in the protection scope of the present invention. Therefore, the protection scope of the present invention shall be based on the protection scope of the claims.

Claims

1. A method for assessing regional rock mass quality considering unloading fracture effects, characterized in that: include: Based on the selected sample area, obtaining sample joint surface data corresponding to all natural rock mass structural surfaces in the sample area; Based on all sample joint surface data, a discrete fracture network model is established to characterize the spatial distribution of joints. Based on the preset ESG microseismic monitoring system, the microseismic data of the monitoring area is obtained, and the geometric parameters of the new fracture surface are calculated based on the microseismic data; Importing the geometric parameters of the newly formed fracture surface into the discrete fracture network model to obtain a three-dimensional discrete fracture network calculation model; Performing a simulated drilling operation in the three-dimensional discrete fracture network calculation model to obtain intersections between each discrete fracture surface and the simulated drilling, and calculating a target RQD value based on the intersections; Wherein, based on the selected sample area, obtaining sample parameters corresponding to all natural rock mass structural surfaces in the sample area includes: Delimiting a sample area, and obtaining sample parameters of the sample area, wherein the sample parameter data includes a half-trace length parameter, a spacing parameter, a strike parameter, a dip parameter, and a dip angle parameter; Converting the sample parameters into a corresponding sample format, and drawing a corresponding joint surface pole diagram and a joint pole isodensity diagram based on the converted sample data; Based on a preset region selection rule, all corresponding regions are selected and grouped in the joint surface pole map and the joint pole equal density map to obtain a corresponding target region; Calculating a disk diameter parameter and a joint surface quantity parameter of a disk corresponding to each target area; The step of calculating the disk diameter parameter and the joint surface quantity parameter of the disk corresponding to each target area includes: Obtain all the half-track length parameters corresponding to all the joint surfaces in each of the target areas, and based on the half-track length parameters, calculate the mean value a of the half-track length parameters in each target area according to the following formula: h : Among them, a g is the overall mean of the total trace length of all joint surfaces, a is the half trace length parameter of the joint surface, μ is the probability density of the trace length of the joint surface, and F(a) is the cumulative probability distribution of the trace length of the joint surface; The mean value of the total trace length in each target area is calculated as follows: Wherein, b is the radius of the disk model; The volume density (c v ) j : Among them, (c i ) j is the linear density of the jth group of joints, E(d 2 ) is the second moment of the disk diameter distribution; e i is the average vector direction of the i-th group of joints; e i is the number of joint surfaces in the i-th group; The number of joint surfaces n is obtained as follows: n=yc v Where V represents the simulated space volume.

2. The method for assessing regional rock mass quality considering unloading fracture effect according to claim 1, characterized in that: The discrete fracture network model characterizing the spatial distribution of joints is established based on all sample joint surface data, including: The dip distribution parameters, dip angle distribution parameters, disk diameter distribution parameters and the number of joint surfaces are sequentially input into the 3DEC software to establish a discrete fracture network model characterizing the spatial distribution of joints, wherein the strike parameters, dip parameters and dip angle parameters adopt a normal distribution, the disk diameter adopts a log-normal distribution, and the disk diameter is set based on a preset diameter range.

3. The method for assessing regional rock mass quality considering unloading fracture effect according to claim 1, characterized in that: The calculating geometric parameters of the newly formed fracture surface based on the microseismic data comprises: Based on the preset point source setting rules, a point source corresponding to the earthquake source is obtained; The long-range P-wave displacement field s of a point source is given by: s=FT=kQT in, ρ represents the density of the rock mass, v represents the P-wave velocity, R represents the observation distance, and Q represents the excitation matrix; Among them, the excitation matrix Q is calculated as follows: 0＝α o α u Where α is the direction cosine vector between the point source and the sensor; The microseismic moment tensor of the initial P wave motion is calculated as follows: Among them, R n represents the observation distance from the point source to the nth sensor, represents the direction cosine vector between the point source and the nth sensor, and T represents the moment tensor; The eigenvalue of the moment tensor T is calculated as follows: Where s represents the displacement of the fracture surface in the direction of movement, A represents the surface area of ​​the new fracture surface, and μ represents the Lame constant. represents the movement direction of the new fracture surface, Indicates the normal direction of the new fracture surface; The eigenvector corresponding to the eigenvalue of the moment tensor is calculated as follows: Vector and The angle between and The angle and and The angle between them is β / 2, which can be calculated as follows: The normal direction of the fracture surface is calculated as follows: and direction of movement According to the spatial position relationship of the new fracture surface occurrence, the fracture surface occurrence data of the new fracture surface geometric parameters are calculated as follows: Where: θ is the strike of the new fracture surface, γ is the inclination angle of the new fracture surface; Calculate the fracture surface radius parameter as follows: Where: z s represents the S wave velocity; p s Indicates the S-wave angular frequency.

4. The method for assessing regional rock mass quality considering unloading fracture effect according to claim 3, characterized in that: The step of importing the geometric parameters of the newly formed fracture surface into the discrete fracture network model to obtain a three-dimensional discrete fracture network calculation model comprises: Inputting the radius parameter of the new fracture surface circle into the 3DEC software; Inputting the strike parameters and dip parameters of the new fracture surface into the 3DEC software; The earthquake source coordinate parameters and the position of the new fracture surface circle are input into the 3DEC software to obtain a three-dimensional discrete fracture network calculation model.

5. The method for evaluating regional rock mass quality considering unloading fracture effect according to claim 4, characterized in that: The method of performing a simulated drilling operation in the three-dimensional discrete fracture network calculation model to obtain the intersection points of each discrete fracture surface and the simulated drilling hole, and calculating the target RQD value based on the intersection points, includes: Use the geometry function of 3DEC to generate numerical simulation drilling and obtain all corresponding target holes; Calculate the intersection points between each target hole and the joint surface to obtain a corresponding intersection point set; Calculating the spacing values ​​between adjacent intersection points in the intersection point set; Traversing and determining whether each spacing value is greater than a preset spacing threshold, if the current spacing value is greater than the preset spacing threshold, taking the current spacing value as the target spacing value; After obtaining all target spacing values, add all target spacing values ​​together to obtain the total length; The sub-target RQD value is calculated as follows: Among them, L i represents the i-th spacing value, and L represents the total number of spacing values ​​in the intersection set; After obtaining the sub-target RQD values ​​corresponding to all the target holes, all the sub-target RQD values ​​are averaged to obtain the target RQD value of the rock mass.

6. A regional rock mass quality assessment device considering unloading fracture effect, characterized in that: include: A sample joint surface data acquisition unit, used to acquire sample joint surface data corresponding to all natural rock mass structural surfaces in the sample area based on the selected sample area; A discrete fracture network model building unit is used to build a discrete fracture network model that characterizes the spatial distribution of joints based on all sample joint surface data; A geometric parameter calculation unit of a newly formed fracture surface, used to obtain microseismic data of a monitoring area based on a preset ESG microseismic monitoring system, and calculate geometric parameters of a newly formed fracture surface based on the microseismic data; A three-dimensional discrete fracture network calculation model establishment unit, used for importing the geometric parameters of the new fracture surface into the discrete fracture network model to obtain a three-dimensional discrete fracture network calculation model; An RQD value calculation unit is used to perform a simulated drilling operation in the three-dimensional discrete fracture network calculation model to obtain the intersection points of each discrete fracture surface and the simulated drilling hole, and calculate the target RQD value based on the intersection points; Wherein, the sample joint surface data acquisition unit comprises: A sample area delineation unit is used to delineate a sample area and obtain sample parameters of the sample area, wherein the sample parameter data includes a half-trace length parameter, a spacing parameter, a strike parameter, a dip parameter and a dip angle parameter; A drawing unit, used for converting the sample parameters into a corresponding sample format, and drawing a corresponding joint surface pole diagram and a joint pole isodensity diagram based on the converted sample data; Selecting a target area unit, for selecting and grouping all corresponding areas in the joint surface pole map and the joint pole equal density map based on a preset area selection rule, to obtain a corresponding target area; A parameter calculation unit, used for calculating a disk diameter parameter and a joint surface quantity parameter of a disk corresponding to each target area; Wherein, the parameter calculation unit includes: A half-trace length parameter mean value calculation unit is used to obtain all half-trace length parameters corresponding to all joint surfaces in each target area, and based on the half-trace length parameters, calculate the mean value a of the half-trace length parameters in each target area according to the following formula: h : Among them, a g is the overall mean of the total trace length of all joint surfaces, a is the half trace length parameter of the joint surface, μ is the probability density of the trace length of the joint surface, and F(a) is the cumulative probability distribution of the trace length of the joint surface; The full trace length average unit is used to calculate the full trace length average in each target area according to the following formula Wherein, b is the radius of the disk model; The volume density unit of the joint surface is used to calculate the volume density (c v ) j : Among them, (c i ) j is the linear density of the jth group of joints, E(d 2 ) is the second moment of the disk diameter distribution; e i is the average vector direction of the i-th group of joints; e i is the number of joint surfaces in the i-th group; The joint surface number calculation unit is used to obtain the joint surface number n according to the following formula: n=You v Where V represents the simulated space volume.

7. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the method for assessing regional rock mass quality considering unloading fracture effects according to any one of claims 1 to 5 is implemented.

8. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, which, when executed by a processor, causes the processor to execute the method for assessing regional rock mass quality taking into account unloading fracture effects as claimed in any one of claims 1 to 5.

Images