Method for obtaining rock mechanical-geometric parameters and holographic scanning system

Abstract

The invention discloses a method for obtaining the geometrical and mechanical parameters of rock samples and a holographic scanning system thereof, wherein the system includes an observation mechanism, a multi-scale penetration mechanism, a grinding mechanism, a rock sample installation mechanism arranged on a three-axis precision motion platform, and an industrial computer controlling the operation mode of each mechanism of the platform Indentation/rotary penetration test, pulse echo signal acquisition, three-dimensional surface topography reconstruction, layer by layer grinding and repeated experiments are carried out. The geometric parameters and corresponding mechanical field parameters are obtained by spatial interpolation of the three-dimensional parameter lattice accumulated by several layers of single-layer rock parameters. The holographic scanning system and method can obtain the real spatial distribution of various media in rock samples. Combined with high performance numerical calculation method, it provides a more scientific method for the analysis of rock mechanical properties, failure and instability.

Claims

1. A method for obtaining rock mechanical and geometric parameters of rock samples, the steps being as follows: S1, acquiring an image of an upper surface of a rock sample by one of two cameras symmetrically and obliquely arranged above the rock sample; S2, dividing the upper surface of the rock sample into N.sub.1×N.sub.1 grids with high resolution by laser, and then using a Berkovich indenter for indentation test in each grid according to a preset loading rate w.sub.1 and an indentation depth L.sub.1, so as to obtain a dynamic signal of displacement changing with time D.sub.1(t) by a displacement measurement system during each indentation process, and a dynamic signal of load changing with time F.sub.1(t) by a pressure sensor which is sitting at a bottom of the rock sample, and then obtaining a three-dimensional mechanical parameter lattice of a corresponding rock sample including an elastic modulus E and a Poisson's ratio μ.sub.1 by processing obtained parameters; and/or, dividing the upper surface of the same rock sample into N.sub.2×N.sub.2 grids with low resolution by laser, and then using a twist drill for rotary cutting penetration test in each grid according to a preset drilling rate w.sub.2 and a drilling depth L.sub.2, so as to obtain a dynamic signal of feed force P(t) changing with time by a feed force sensor connected with the twist drill during each cutting process, a dynamic signal of torque M(t) changing with time by a torque sensor set at an adjacent side of the twist drill during each cutting process, and a dynamic signal of displacement D.sub.2(t) changing with time by a displacement measurement system during each cutting process, and then obtaining three-dimensional mechanical parameters including an elastic modulus E, a cohesive force C and an internal friction angle φ and a compressive strength R.sub.p; S3, acquiring a pulse echo signal of the rock sample by an ultrasonic sensor excitation device and an ultrasonic sensor receiving device which are symmetrically arranged on both sides of the rock sample, so as to obtain a longitudinal wave velocity V.sub.p and a transverse wave velocity V.sub.s of an ultrasonic wave passing through the rock sample and a density ρ of the rock sample, and then calculating a dynamic elastic modulus Ed and a dynamic Poisson's ratio μ.sub.d of the rock sample; $\begin{array}{cc}{E}_d=\frac{\rho V_s^2\left(3V_p^2-4V_s^2\right)}{V_p^2-V_s^2},& {\mu }_d=\frac{V_p^2-2V_s^2}{2\left(V_p^2-V_s^2\right)}\end{array};$ S4, moving the rock sample to a position identical to step S1 again, and then projecting a structured light coding pattern on the upper surface of the rock sample, so as to collect images of the upper surface of the rock sample by the two cameras which are symmetrically and obliquely arranged above the rock sample, and reconstruct a three-dimensional morphology of the upper surface of the rock sample further; S5, polishing the upper surface of the rock sample, and lubricating and cooling by running water until the upper surface of the rock sample returning to a smooth plane; S6, repeating steps S1˜S5 to grind the rock sample layer by layer and carrying out identical parameter acquisition experiment until the rock sample reaches a last layer of a preset depth, so as to acquire a geometry of a complete rock sample and a 3D parameter lattice of mechanical parameters from an accumulation of parameters of several layers of single-layer rock; S7, using an IDW spatial interpolation algorithm to process the 3D lattice of mechanical parameters obtained from step S6: $R\left(x\right)=\underset{i=0}{\overset{n}{.Math.}}\frac{{\omega }_i\left(x\right)u_i}{S},$ wherein, ${\omega }_i\left(x\right)=\frac{1}{\ln \left(d\left(x,x_i\right)\right)};S=\underset{j=0}{\overset{N}{.Math.}}{\omega }_j\left(x\right);$ x.sub.i is an i.sup.th sample point; x is an unsolved point; ω.sub.i(x) is a weight of x.sub.i; u.sub.i is an actual value of the i.sup.th sample point; R(x) is mechanical parameters of rock obtained from different mechanical tests (indentation test, rotary cutting penetration test); d(x,x.sub.i) is a distance between the sample point and the unsolved point:
d(x,x.sub.i)=√{square root over ((x−x.sub.i).sup.2+(y−y.sub.i).sup.2+(z−z.sub.i).sup.2)}, wherein, x.sub.i, y.sub.i, z.sub.i are coordinates of the i.sup.th sample point on X, Y, and Z axes, respectively; obtaining a virtual lattice of mechanical parameters by weighted averages using distances between interpolation points and the sample points as weight, and then obtaining three-dimensional mechanical parameter fields of the complete rock sample, including a compressive strength field, a cohesive force field, an internal friction angle field, an elastic modulus field and a Poisson's ratio field.

2. The method for obtaining the geometrical and mechanical parameters of rock samples according to claim 1, wherein when the indentation test is conducted, a value of N.sub.1 varies according to a shape of the Berkovich indenters, and ranges from 50×50 to 100×100, and an area of the Berkovich indenter is 1/200˜ 1/50 of an area of the upper surface of the rock sample, and the preset depth is from 0.5 mm to 2 mm with quasi-static loading.

3. The method for obtaining the geometrical and mechanical parameters of rock samples according to claim 1, wherein when the indentation test is conducted, the method to obtain the three-dimensional mechanical parameters is as follows: correlating the elastic modulus of the rock sample and the Berkovich indenter by a reduction of elastic modulus E.sub.r, and a formula for calculating an indentation hardness H.sub.I from the indentation test is: $H_I=\frac{P_{\max }}{A_c},$ wherein, a formula of the reduction of elastic modulus E.sub.r is: $E_r=\frac{\sqrt{\pi }}{2}\frac{S}{\sqrt{A_C}},$ wherein, a formula of elastic modulus E.sub.i is: $E_i=\left(1-{\upsilon }_i^2\right){\left(\frac{1}{E_r}-\frac{1-{\upsilon }_s^2}{E_s}\right)}^{-1},$ wherein, P.sub.max is a maximum indentation force; E.sub.s is an elastic modulus of the Berkovich indenter; v.sub.s is a Poisson's ratio of the indenter; u.sub.i is a Poisson's ratio of the rock sample; wherein, A.sub.c is a projected area of the Berkovich indenter:
A.sub.C=24.5h.sub.C.sup.2, wherein, h.sub.c is a pressure depth, which is obtained by a following formula: $h_C=L-\frac{3P_{\max }}{4S},$ wherein, S is an unloading stiffness, which is obtained from a load-displacement curves at an unloading stage, and a tangent slope at one third of a curve is taken as the unloading stiffness.

4. The method for obtaining the geometrical and mechanical parameters of rock samples according to claim 1, wherein when the rotary cutting penetration test is conducted, details of obtaining three-dimensional mechanical parameters of the rock sample are as follows: assuming a feed displacement per revolution being proportionate to a bit pressure P and a torque M, $\frac{1}{K_p}=\frac{1}{2}{\tau }_{f}R{\delta }_c\tan \left(\frac{\pi }{2}-\beta \right),$ $\frac{1}{K_m}=\frac{1}{2}{\tau }_f{R}^2{\delta }_c^2,$ wherein, K.sub.p is a slope of a feed displacement-bit pressure curve, which is obtained from signals P(t) and D.sub.2(t); K.sub.m is a slope of the feed displacement-bit torque curve, which is obtained from M(t) and D.sub.2(t); β is a mean friction angle of the bit; τ.sub.f is the tensile strength of the rock sample; R is a radius of the bit, and a cutting coefficient is: ${\delta }_c=\frac{\cos \beta }{\sin \varphi .Math.\cos \left(\varphi +\beta \right)},$ wherein, a mean friction angle of the bit β is: $\beta =\frac{\pi }{2}-\mathrm{arc}\tan \frac{R{K}_m}{K_p},$ wherein, a cutting angle is: $\varphi =45°-\frac{\beta }{2},$ wherein, a friction angle φ is:
φ=90°−2ϕ, wherein, a normal stress σ and a tangential stress τ are obtained from: $\sigma =\frac{P}{2\mathrm{aR}\sec \theta \tan \varphi },$ $\tau =\frac{P}{2\mathrm{aR}\sec \theta },$ wherein, θ is the bit angle, according to a Mohr-coulomb criterion, the cohesive force of the rock sample is obtained as follows: $C=\tau -\sigma \tan \phi =\frac{P}{2\mathrm{aR}\sec \theta }-\frac{P}{2\mathrm{aR}\sec \theta \tan \varphi },$ wherein, a compressive strength of the rock sample is: $R_p=\frac{2C\cos \phi }{1-\sin \phi },$ wherein, according to an energy conservation theory, a total work done by the bit is mainly used to convert into an internal energy W.sub.inter, a work done by rock cutting W.sub.cutting, and a work done by pressing into rock sample W.sub.press, and a corresponding formula is:
W.sub.total=W.sub.inter+W.sub.cutting+W.sub.press, wherein, according to a bit pressure-time curve P.sub.(t), a drilling speed v, a bit torque-time curve M.sub.(t), and a rotation speed of the bit ω, the internal energy W.sub.inter is: $W_{total}=\underset{t}{\int }{P}_{\left(t\right)}.Math.v.Math.\mathrm{dt}+\underset{t}{\int }{M}_{\left(t\right)}.Math.\omega .Math.\mathrm{dt},$ wherein, the drilling speed is obtained from: $V=\frac{L}{t},$ wherein, L is a predetermined depth, t is a drilling time; wherein, the internal energy W.sub.inter is obtained from: $W_{\mathrm{inter}}=\underset{\Omega 1}{\int \int \int }{T}_1\left(\Delta T\right)C_{1}dv+\int \underset{\Omega 2}{\int }\int T_2\left(\Delta T\right)C_2\mathrm{dv},$ wherein, T.sub.1(ΔT) is a temperature distribution function of the twist drill; T.sub.2(ΔT) is a temperature distribution function of the rock sample; C.sub.1 is a specific heat of the twist drill; C.sub.1 is a specific heat of the rock sample; Ω.sub.1 is a domain size of the twist drill; Ω.sub.2 is a domain size of the rock sample; T.sub.1(ΔT) is a temperature difference of the twist drill before and after drilling; T.sub.2(ΔT) is a temperature difference of the rock sample before and after drilling; wherein, the work done by cutting the rock sample is divided into two parts: a work done by a crushing torque, and a work done by a friction torque:
W.sub.Cutting=∫M.sub.Crush.Math.ωdt+∫M.sub.Friction.Math.ωdt, and the elastic modulus of the rock sample E is obtained from: $E=\frac{F_2^{2}L}{2\pi R^{2}C\left(W_{\mathrm{total}}-W_{\mathrm{inter}}-W_{\mathrm{cutting}}\right)}.$

5. The method for obtaining the geometrical and mechanical parameters of rock samples according to claim 1, wherein when the indentation test is conducted, the value of N.sub.2 varies according to a shape of drill bits, and ranges from 10×10 to 50×50, and a diameter of the twist drill is 1/50 to 1/10 of an area of the upper surface of the rock sample, and the drilling depth L2 is from 0.5 mm to 2 mm.

6. The method for obtaining the geometrical and mechanical parameters of rock samples according to claim 1, wherein in step S3, a wavelength emitted by the ultrasonic sensor excitation device is one tenth of a path length through the rock sample.

7. A holographic scanning system for realizing the method for obtaining the geometrical and mechanical parameters of rock samples according to claim 6, wherein an equipment includes an observation mechanism arranged on a three-axis precision motion platform, a multi-scale penetration mechanism, a grinding mechanism, a rock sample installation mechanism, and an industrial computer for controlling an operation of the three-axis precision motion platform, the observation mechanism, the multi-scale penetration mechanism, the grinding mechanism and the rock sample installation mechanism; wherein, the three-axis precision motion platform comprises a first motion mechanism which is configured to realize up and down motion, a second motion mechanism which is arranged on the first motion mechanism and is configured to realize left and right motion, a third motion mechanism which is arranged on the second motion mechanism and is configured to realize back and forth motion, and a cover plate which is located above the three motion mechanisms; wherein the observation mechanism, the multi-scale penetration mechanism, and the grinding mechanism are successively installed on a floor of the cover plate, and the rock sample installation mechanism is installed on the third moving mechanism; the observation mechanism comprises a micro projector arranged on a bottom surface of the cover plate and a first camera, and a second camera symmetrically arranged on both sides of the micro projector; wherein a lens of the first camera and a lens of the second camera are arranged in an oblique manner towards the micro projector; the multi-scale penetration mechanism comprises a rotary cutting penetration device at a front side, an indentation detection device at a rear side, and a displacement measurement system arranged between the rotary cutting penetration device and the indentation detection device; wherein, the rotary cutting penetration device comprises a first rigid servo motor, a first ball screw, a feed force sensor and a standard twist drill which are connected from top to bottom in turn, and a displacement measurement system arranged between the rotary cutting penetration device and the indentation detection device, the indentation detection device includes a second rigid servo motor, a second ball screw, a vertical rigid spindle and a triangular shell type indenter connected from top to bottom in turn; the rock sample installation mechanism comprises a square steel fixture for fixing the rock sample, an ultrasonic sensor exciting device and an ultrasonic sensor receiving device symmetrically arranged on both sides of the square steel fixture; a pressure sensor for detecting a pressure on the rock sample is also arranged on the third moving mechanism of the square steel fixture.

8. The holographic scanning system according to claim 7, wherein the grinding mechanism comprises a grinding wheel and a drying device; wherein the drying device is an annular disc-shaped structure sheathed on an outer side of the grinding wheel, on which a plurality of hot air pipes are arranged along a annular direction, and each hot air pipe is connected with a hot air circulation box through a pipeline, so that a hot air is vertically ejected from the upper surface of the rock sample through a plurality of hot air pipes.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) The accompanying drawings are included to provide a further understanding of the disclosure, and are incorporated in and constitute a part of this specification. The drawings illustrate exemplary embodiments of the disclosure and, together with the description, serve to explain the principles of the disclosure.

(2) FIG. 1 is a schematic of rock mechanical-geometric parameters and holographic scanning system;

(3) FIG. 2. is a main view of the three-axis precision motion platform in rock mechanical-geometric parameters and holographic scanning system;

(4) FIG. 3. is a top view of the three-axis precision motion platform in rock mechanical-geometric parameters and holographic scanning system;

(5) FIG. 4. is a left view of the three-axis precision motion platform in rock mechanical-geometric parameters and holographic scanning system;

(6) FIG. 5. is a schematic of observation mechanism in rock mechanical-geometric parameters and holographic scanning system;

(7) FIG. 6. is a flowchart of rock mechanical-geometric parameters and holographic scanning system.

DESCRIPTION OF THE EMBODIMENTS

(8) The method for obtaining rock mechanical and geometric parameters of the present invention will be described in detail below with reference to the embodiments and the accompanying drawings.

Embodiment 1

(9) As shown in FIG. 1, a holographic scanning system, which includes a observation mechanism 4 arranged on a three-axis precision motion platform 7, a multi-scale penetration mechanism 3, a grinding mechanism 2, a rock sample installation mechanism, and an industrial computer for controlling the operation of three-axis precision motion platform 7, the observation mechanism 4, the multi-scale penetration mechanism 3, the grinding mechanism 2 and the rock sample installation mechanism; wherein,

(10) As shown in FIG. 2, the three-axis precision motion platform 7 is purchased from common commercially available equipment, which includes the first support column a, the second support column B, the third support column C and the fourth support column D. between the first support column A and the second support column B, and between the third support column C and the fourth support column D, there is a horizontal longitudinal rail 13 on each of the two horizontal rails, Each longitudinal rail 13 is provided with a first sliding block which can move back and forth on the longitudinal rail, at the same time, a vertical rail 16 is respectively arranged on the side wall of each supporting column which is in contact with the horizontal longitudinal beam, and the two ends of each horizontal longitudinal beam are respectively provided with a second sliding block which is matched with the vertical rail, so that the horizontal longitudinal beam can move back and forth along the vertical rail 16 in the axial direction The sliding block is fixed with a first bearing platform 15, on which two transverse rails 14 are symmetrically erected along the direction perpendicular to the horizontal longitudinal beam, each transverse rail 14 is provided with a third sliding block matched with it, and the second bearing platform 15 is fixed on the two third sliding blocks, so that the second bearing platform 15 can use the third sliding block to slide back and forth along the two transverse rails 14; on each supporting column, each horizontal rail 14 is fixed with a second bearing platform 15 The longitudinal beam, the first bearing platform 15 and the second bearing platform 15 are also provided with a plurality of servo motors for controlling the movement of two first sliding blocks, four second sliding blocks and two third sliding blocks respectively, and each servo motor is connected with the corresponding sliding blocks through a ball screw to control the movement of the sliding blocks;

(11) As shown in FIG. 1, a horizontal cover plate is fixed on the top surface of the first support column a, the second support column B, the third support column C and the fourth support column D, so that the observation mechanism, the multi-scale penetration mechanism and the grinding mechanism are arranged on the bottom surface of the cover plate in turn along the axial direction of the horizontal longitudinal beam;

(12) As shown in FIG. 1, FIG. 2, and FIG. 5, the observation mechanism 4 comprises a micro projector 19 arranged on the bottom surface of the cover plate and a first camera 17, and a second camera 18 symmetrically arranged on both sides of the micro projector 19; wherein, the micro projector 19, the first high-precision digital camera 17 and the second high-precision digital camera 18 are spaced on the same straight line along the axial direction of the horizontal longitudinal beam, and the lenses of the first high-precision digital camera 17 and the second high-precision digital camera 18 are arranged towards the micro projector 19 in an inclined way of 45°; wherein, the micro projector 19 is used to project a number of points to the top surface of the rock sample 8 The first high-precision digital camera 17 and the second high-precision digital camera 18 are used to collect the top photos of the rock sample 8 with digital images, specifically the three-dimensional morphology of the indentation (groove) on the rock surface, so as to draw the three-dimensional image of the rock surface in the computer;

(13) As shown in FIG. 1, the multi-scale penetration mechanism 3 comprises a rotary cutting penetration device at the front side, an indentation detection device at the rear side, and a displacement measurement system arranged between the rotary cutting penetration device and the indentation detection device; wherein, the indentation detection device is used to test the micro scale of rock sample 8 at first, the rotary cutting penetration device is used to test the macro scale of rock sample 8 again, and the displacement measurement system is used to test the dynamic signal of displacement changing with time during drilling and indentation of bit or indenter;

(14) Wherein, the rotary cutting penetration device comprises a first rigid servo motor, a first ball screw, a feed force sensor and a standard twist drill which are connected from top to bottom in turn, and a displacement measurement system arranged between the rotary cutting penetration device and the indentation detection device; the standard twist bit is equipped with various diameter bits, including 0.5 mm diameter standard twist bit, 1 mm diameter standard twist bit, 2 mm diameter standard twist bit and 10 mm diameter standard twist bit; a power amplifier is set between the rigid servo motor and the ball screw to provide the required drilling conditions for the bit; the working principle is as follows: The rigid servo motor drives the ball screw and loads the bit through the power amplifier to realize the rotary cutting of the top surface of rock sample 8; at the same time, the feed force sensor is used to measure the axial WOB of the bit in real time, so as to realize the high precision loading of the bit; the torque sensor is connected with the drill bit through the anti-drying shielded cable, the two sensors can obtain the dynamic signals of the feed force F measured by the feed force sensor and the torque T measured by the torque sensor, which change with time, for subsequent data analysis;

(15) Wherein, the indentation detection device includes a second rigid servo motor, a second ball screw, a vertical rigid spindle and a triangular shell type indenter connected from top to bottom in turn; the triangular shell type indenter is equipped with replaceable indenters of various diameters, including the triangular indenter with side length of 0.5 mm, the triangular indenter with side length of 1 mm and the triangular indenter with side length of 2 mm; each triangular indenter is made of tungsten carbide; the second rigid servo motor is connected with the second driving ball screw through a power amplifier to realize the second rigid servo motor; the working principle of the indentation detection device is as follows: the vertical rigid spindle is driven by the motor to transmit the pressure to the triangular bayonet indenter, which makes indentation on the rock surface under the set load, and uses the displacement measurement system to measure the indentation depth in real time, so as to obtain each size triangle Load displacement curve of indenter;

(16) As shown in FIG. 1, the grinding mechanism 2 comprises a grinding wheel and a drying device; wherein the drying device is an annular disc-shaped structure sheathed on the outer side of the grinding wheel, on which a plurality of hot air pipes are arranged along the annular direction, and each hot air pipe is connected with the hot air circulation box through a pipeline, so that the hot air is vertically ejected from the top surface of the rock sample 8 through a plurality of hot air pipes;

(17) As shown in FIG. 1, FIG. 2, and FIG. 3, the rock sample installation mechanism comprises a square steel fixture 10 for fixing the rock sample 8, an ultrasonic sensor exciting device 9 and an ultrasonic sensor receiving device 11 symmetrically arranged on both sides of the square steel fixture; wherein, the twist bit of the rotary cutting penetration device matches the size of the shell type indenter of the indentation detection device, and the size of the second bearing platform 15 meets the setting of the fixed rock sample 8 with the size of 100 mm×100 mm×180 mm; the ultrasonic sensor exciting device 9 and the ultrasonic sensor receiving device 11 respectively contact with the side wall of the rock sample 8, and the contact surface is coated with a layer of Vaseline; the ultrasonic sensor exciting device 9 and the ultrasonic sensor receiving device 11 have the functions of ultra-low receiving noise and high performance pulse control, which can generate the excitation signal and receive the pulse echo signal; in order to meet the test requirements, the ultrasonic sensor excitation device 9 and the ultrasonic sensor receiving device 11 should ensure that the wavelength of the ultrasonic wave emitted and received by them meets the following requirements: the width of the rock sample is more than ten times of the wavelength, so as to meet the physical premise of an approximate infinite body; in addition, a pressure sensor for detecting the pressure on the rock sample is also arranged on the third moving mechanism of the square steel fixture 10;

(18) The industrial computer controls the three-axis precision motion with multiple driving mechanisms of the three-axis precision motion platform, two cameras and projection equipment of the observation mechanism, each servo motor and sensor of the multi-scale penetration mechanism, grinding wheel and drying device of the grinding mechanism, ultrasonic sensor excitation device 9, ultrasonic sensor receiving device 11 and pressure sensor on the rock sample installation mechanism The platform, observation mechanism, multi-scale penetration mechanism, grinding mechanism and rock sample installation mechanism operate normally according to the above operation mode, and receive the real-time data sent back by the corresponding equipment for subsequent technical personnel to process;

(19) When the holographic scanning system of rock geometric mechanical parameters is used, the NI control and testing system 12 can be installed to assist the three-axis precision motion platform 7. The rock sample installation mechanism with rock sample 8 is arranged in the order of observation mechanism 4, indentation detection device of multi-scale penetration mechanism 3, rotary cutting penetration device of multi-scale penetration mechanism 3, observation mechanism 4 and grinding mechanism 4 The grinding mechanism 2 performs in sequence and circulates to realize the step-by-step and layer by layer scanning of rock sample 8.

Embodiment 2

(20) As shown in FIG. 6, a method for obtaining rock mechanical and geometric parameters by using the holographic scanning system in Embodiment 1, the steps are as follows:

(21) S1. Placing a rock sample 8 taken from the target test rock on the second bearing platform 15, and the bottom of the rock sample 8 is consolidated, and the rock sample is moved to the lower part of the observation mechanism through the three-axis precision motion platform 7; specifically, the rock sample 8 coincides with the central axis of the micro projector 19; then, one of the two cameras symmetrically and obliquely arranged above the rock sample is used the image of the upper surface of the rock sample is captured by a camera;

(22) S2. Dividing the upper surface of the rock sample into 50×50 grids with high resolution by laser, and then using the triangle shell type Berkovich indenter with side length of 0.5 mm for indentation test in each grid according to the preset loading rate w.sub.1=10 mm/s (Quasi static state) and indentation depth L.sub.1=0.5 mm, so as to obtain the dynamic signal of displacement changing with time D.sub.1(t) by the displacement measurement system during each indentation process and the dynamic signal of load changing with time F.sub.1(t) by the pressure sensor which is sitting at the bottom of the rock sample;

(23) Specifically, because the indenter is not completely rigid, the mechanical parameters of rock can not be obtained directly from the load displacement data; therefore, correlating the elastic modulus of rock sample and indenter by the reduction of elastic modulus E.sub.r, and the formula for calculating indentation hardness H.sub.I from indentation test is:

(24) $H_I=\frac{P_{\max }}{A_c},$

(25) wherein, the formula of reduction of elastic modulus E.sub.r is:

(26) 0$E_r=\frac{\sqrt{\pi }}{2}\frac{S}{\sqrt{A_C}},$

(27) wherein, the formula of elastic modulus E.sub.i is:

(28) $E_i=\left(1-{\upsilon }_i^2\right){\left(\frac{1}{E_r}-\frac{1-{\upsilon }_s^2}{E_s}\right)}^{-1},$

(29) wherein, P.sub.max is the maximum indentation force; E.sub.s is the elastic modulus of indenter; v.sub.s is the Poisson's ratio of indenter; u.sub.i is the Poisson's ratio of rock sample;

(30) wherein, A.sub.c is the projected area of indenter:
A.sub.C=24.5h.sub.C.sup.2,

(31) wherein, h.sub.c is the pressure depth, which is obtained by the following formula:

(32) $h_C=L-\frac{3P_{\max }}{4S},$

(33) Wherein, S is the unloading stiffness, which is obtained from the load-displacement curves at the unloading stage, and the tangent slope at one third of the curve is taken as the unloading stiffness;

(34) And then, obtaining three-dimensional mechanical parameter lattice of the corresponding rock sample including elastic modulus E and Poisson's ratio μ.sub.i by processing the above parameters;

(35) S3. Acquiring a pulse echo signal of the rock sample by the ultrasonic sensor excitation device and ultrasonic sensor receiving device which are symmetrically arranged on both sides of the rock sample, so as to obtain the longitudinal wave velocity V.sub.p and transverse wave velocity V.sub.s of the ultrasonic wave passing through the rock sample and the density ρ of the rock sample, and then calculating the dynamic elastic modulus Ed and dynamic Poisson's ratio μ.sub.d of the rock sample;

(36) $\begin{array}{cc}{E}_d=\frac{\rho V_S^2\left(3V_p^2-4V_s^2\right)}{V_p^2-V_s^2},& {\mu }_d=\frac{V_p^2-2V_s^2}{2\left(V_p^2-V_s^2\right)}\end{array};$

(37) S4. Moving the rock sample to the same position as step S1 again, and then projecting the structured light coding pattern on the upper surface of the rock sample, so as to collect images of the upper surface of the rock sample by two cameras which are symmetrically and obliquely arranged above the rock sample, and reconstruct the three-dimensional morphology of the upper surface of the rock sample further;

(38) S5. Polishing the surface of the rock sample, and lubricating and cooling by running water until the upper surface of the rock sample returning to a smooth plane;

(39) S6. Repeating steps S1˜S5 to grind the rock sample layer by layer and carrying out the same parameter acquisition experiment until the rock sample reaches the last layer of the preset depth, so as to acquire the 3D parameter lattice of the complete rock sample geometry and mechanics from the accumulation of several layers of single-layer rock parameters;

(40) S7. Using IDW spatial interpolation algorithm to process the three-dimensional lattice of mechanical parameters obtained from step S6:

(41) $R\left(x\right)=\underset{i=0}{\overset{n}{.Math.}}\frac{{\omega }_i\left(x\right)u_i}{S},$

(42) Wherein,

(43) $\begin{array}{cc}{\omega }_i\left(x\right)=\frac{1}{\ln \left(d\left(x,x_i\right)\right)};& S=\underset{j=0}{\overset{N}{.Math.}}{\omega }_j\left(x\right);\end{array}$
x.sub.i is the i.sup.th sample point; x is the unsolved point; ω.sub.i(x) is the weight of x.sub.i; u.sub.i is the value of i.sup.th actual sample point; R(x) is the mechanical parameters of rock obtained from different mechanical tests (indentation test, rotary cutting penetration test); d(x,x.sub.i) is the distance between the measured point and unsolved point:
d(x,x.sub.i)=√{square root over ((x−x.sub.i).sup.2+(y−y.sub.i).sup.2+(z−z.sub.i).sup.2)},

(44) Wherein, x.sub.i, y.sub.i, z.sub.i are the coordinates of the ith sample point on the X, Y, and Z axes, respectively;

(45) Obtaining a virtual lattice of mechanical parameters by weighted averages using distances between interpolation points and the sample points as weight, and then obtaining the three-dimensional mechanical parameter fields of the whole rock sample, including compressive strength field, cohesive force field, internal friction angle field, elastic modulus field and Poisson's ratio field.

(46) In the same way, take two rock samples which are the same as those in the present example, and test them in the same step by using a 1 mm side length bayonet indenter and a 2 mm side length bayonet indenter respectively, so as to obtain the geometric mechanical parameters of rock samples under different load conditions. For example, when the size of the indenter is 0.5 mm, 1 mm and 2 mm, the mesh is divided into 50×50, 80×80 and 100×100.

Embodiment 3

(47) A method for obtaining rock mechanical and geometric parameters by using the holographic scanning system in Embodiment 1; wherein the steps are the same except that step S2 is different from embodiment 1.

(48) Specifically, the specific steps of step S2 in this embodiment are:

(49) Dividing the upper surface of the same rock sample into 40×40 grids with low resolution by laser, and then using the twist drill for rotary cutting penetration test in each grid according to a preset drilling rate w.sub.2=75 mm/min and drilling depth L.sub.2=5 mm, so as to obtain the dynamic signal of feed force P(t) changing with time by the feed force sensor connected with the twist drill during each cutting process, the dynamic signal of torque M(t) changing with time by the torque sensor set at the adjacent side of the twist drill during each cutting process, and the dynamic signal of displacement D.sub.2(t) changing with time by the displacement measurement system during each cutting process,

(50) Assuming the feed displacement per revolution is proportionate to the bit pressure P and torque M,

(51) $\frac{1}{K_p}=\frac{1}{2}{\tau }_{f}R{\delta }_c\tan \left(\frac{\pi }{2}-\beta \right),$$\frac{1}{K_m}=\frac{1}{2}{\tau }_f{R}^2{\delta }_c^2,$

(52) Wherein, K.sub.p is the slope of the feed displacement-bit pressure curve, which is obtained from signals P(t) and D.sub.2(t); K.sub.m is the slope of the feed displacement-bit torque curve, which is obtained from M(t) and D.sub.2(t); β is the mean friction angle of the bit; τ.sub.f is the tensile strength of rock sample; R is radius of the bit, and the cutting coefficient is:

(53) ${\delta }_c=\frac{\cos \beta }{\sin \varphi .Math.\cos \left(\varphi +\beta \right)},$

(54) wherein, the mean friction angle of the bit β is:

(55) $\beta =\frac{\pi }{2}-\mathrm{arc}\tan \frac{R{K}_m}{K_p},$

(56) Wherein, cutting angle is:

(57) $\varphi =45°-\frac{\beta }{2},$

(58) Wherein, the friction angle φ is:
φ=90°−2ϕ,

(59) Wherein, the normal stress σ and tangential stress τ is obtained from:

(60) 0$\sigma =\frac{P}{2\mathrm{aR}\sec \theta \tan \varphi },$$\tau =\frac{P}{2\mathrm{aR}\sec \theta },$

(61) Wherein, θ is the bit angle. According to the Mohr-coulomb criterion, the cohesive of rock sample is obtained as follows:

(62) $C=\tau -\sigma \tan \phi =\frac{P}{2\mathrm{aR}\sec \theta }-\frac{P}{2\mathrm{aR}\sec \theta \tan \varphi },$

(63) wherein, the compressive strength of rock sample is:

(64) $R_p=\frac{2C\cos \phi }{1-\sin \phi },$

(65) Wherein, according to the energy conservation theory, the total work done by bit is mainly used to convert into internal energy W.sub.inter, work done by rock cutting W.sub.cutting, and work done by pressing into rock sample W.sub.press, and the corresponding formula is:
W.sub.total=W.sub.inter+W.sub.cutting+W.sub.press,

(66) Wherein, according to the bit pressure-time curve P.sub.(t), the drilling speed v, the bit torque-time curve M.sub.(t), and rotation speed of bit ω, the internal energy W.sub.inter is:

(67) $W_{total}=\underset{t}{\int }{P}_{\left(t\right)}.Math.v.Math.\mathrm{dt}+\underset{t}{\int }{M}_{\left(t\right)}.Math.\omega .Math.\mathrm{dt},$

(68) Wherein, the drilling speed is obtained from:

(69) $V=\frac{L}{t},$

(70) Wherein, L is the predetermined depth, t is the drilling time;

(71) Wherein, the internal energy W.sub.inter is obtained from:

(72) $W_{inter}=\int \underset{\mathrm{\Omega 1}}{\int }\int T_1\left(\Delta T\right)C_{1}dv+\int \underset{\mathrm{\Omega 2}}{\int }\int T_2\left(\Delta T\right)C_2\mathrm{dv},$

(73) Wherein, T.sub.1(ΔT) is the temperature distribution function of drill; T.sub.2(ΔT) is the temperature distribution function of rock samples; C.sub.1 is the specific heat of drill; C.sub.1 is the specific heat of rock samples; Ω.sub.1 is the domain size of drill; Ω.sub.2 is the domain size of rock samples; T.sub.1(ΔT) is the temperature difference of drill before and after drilling; T.sub.2(ΔT) is the temperature difference of rock samples before and after drilling;

(74) Wherein, the work done by cutting rock samples can be divided into two parts: the work done by the crushing torque, and the work done by friction torque:
W.sub.Cutting=∫M.sub.Crush.Math.ωdt+∫M.sub.Friction.Math.ωdt,

(75) And elastic modulus of rock samples E is obtained from:

(76) $E=\frac{F_2^{2}L}{2\pi R^{2}C\left(W_{\mathrm{total}}-W_{\mathrm{inter}}-W_{\mathrm{cutting}}\right)},$

(77) And then obtaining the three-dimensional mechanical parameters including elastic modulus E, cohesive force C and internal friction angle φ and compressive strength R.sub.p.

(78) In the same way, take another three rock samples which are the same as those in this example, and test them in the same steps by using the twist drill bit with a diameter of 1 mm, the twist drill bit with a diameter of 2 mm and the twist drill bit with a diameter of 10 mm respectively, so as to obtain the geometric mechanical parameters of rock samples under different WOB conditions. For example, when the diameter of the twist drill is 0.5 mm, 1 mm, 2 mm and 10 mm, the mesh can be set as 10×10, 25×25, 35×35 and 50×50.

(79) It will be apparent to those skilled in the art that various modifications and variations can be made to the disclosed embodiments without departing from the scope or spirit of the disclosure. In view of the foregoing, it is intended that the disclosure covers modifications and variations provided that they fall within the scope of the following claims and their equivalents.