METODA MERJENJA GEO-STATIČNIH NAPETOSTI, KI TEMELJI NA VGRAJENEM VRTALNEM IN OPTIČNO MIKROSKOPSKEM SLIKOVNEM SISTEMU Izvleček Konvencionalne metode merjenja geostatičnih napetosti so omejene zaradi pomanjkljivosti, vključno z izmerljivo globino, kompleksnostjo in dolgim trajanjem izvedbe. Da bi rešili te probleme in dosegli merjenje geostatičnih napetosti v globokih vodnjakih v pogojih kompleksnih visokih tlakov in visokih temperatur, je predlagana nova merilna metoda za geostatične napetosti na osnovi integriranega sistema vrtanja in optičnega mikroskopa. Njena inovativna integrirana konstrukcija odpravlja probleme, povezane s kompleksnimi postopki in globinskimi omejitvami, ter se izogiba lezenju kamnin zaradi dolgotrajnosti postopka, kar bistveno izboljša natančnost in obseg meritev. Deluje s pomočjo mikroskopskega slikanja in neposrednimi sondami, da zajame spremembe zunanje površine prereza vrtine pred in po zmanjšanju napetosti. Rezultirajoče slike se analizirajo z iskalnimi krogi, da se pridobijo položaji vrhov sond, ki so lahko vstavljeni v elipse, ki opisujejo obrise, in izračunajo stanje napetosti. Veljavnost in natančnost metode sta bili preverjeni z laboratorijskimi preizkusi in aplikacijami na terenu v vrtini ZK1. Rezultati kažejo, da: (1) se integrirani sistem lahko uporablja za merjenje mikrometerskih deformacij; (2) pristop iskalnega kroga lahko natančno pridobi položaje vrhov sond; (3) metoda merjenja napetosti, ki temelji na tem sistemu, je točna in izvedljiva. Jinchao Wang Chinese Academy of Sciences, Institute of Rock and Soil Mechanics, Key Laboratory of Geomechanics and Geotechnical Engineering Wuhan, Kitajska University of Chinese Academy of Sciences Beijing ,China Chuanying Wang Chinese Academy of Sciences, Institute of Rock and Soil Mechanics, Key Laboratory of Geomechanics and Geotechnical Engineering Wuhan, Kitajska Zengqiang Han (vodilni avtor) Chinese Academy of Sciences, Institute of Rock and Soil Mechanics, Key Laboratory of Geomechanics and Geotechnical Engineering Wuhan, Kitajska E-posta: zqhan@whrsm.ac.cn Yiteng Wang Chinese Academy of Sciences, Institute of Rock and Soil Mechanics, Key Laboratory of Geomechanics and Geotechnical Engineering Wuhan, Kitajska University of Chinese Academy of Sciences Beijing, Kitajska Xinjian Tang Chinese Academy of Sciences, Institute of Rock and Soil Mechanics, Key Laboratory of Geomechanics and Geotechnical Engineering Wuhan, Kitajska Ključne besede geostatična napetost; sonda; optična mikroskopija; merjenje med vrtanjem; iskalni krog 16. Acta Geotechnica Slovenica, 2018/1 A GEOSTRESS MEASUREMENT METHOD BASED ON AN INTEGRATED DRILLING AND OPTICAL MICROSCOPIC IMAGING SYSTEM Jinchao Wang Chinese Academy of Sciences, Institute of Rock and Soil Mechanics, Key Laboratory of Geomechanics and Geotechnical Engineering Wuhan, China University of Chinese Academy of Sciences Beijing ,China Chuanying Wang Chinese Academy of Sciences, Institute of Rock and Soil Mechanics, Key Laboratory of Geomechanics and Geotechnical Engineering Wuhan, China Zengqiang Han (corresponding author) Chinese Academy of Sciences, Institute of Rock and Soil Mechanics, Key Laboratory of Geomechanics and Geotechnical Engineering Wuhan, China E-mail: zqhan@whrsm.ac.cn Yiteng Wang Chinese Academy of Sciences, Institute of Rock and Soil Mechanics, Key Laboratory of Geomechanics and Geotechnical Engineering Wuhan, China University of Chinese Academy of Sciences Beijing, China Xinjian Tang Chinese Academy of Sciences, Institute of Rock and Soil Mechanics, Key Laboratory of Geomechanics and Geotechnical Engineering Wuhan, China Keywords geostress; probe; optical microscopy; measurement while drilling; search circle https://doi.org/10.18690/actageotechslov.15.1.17-27.2018 Abstract Conventional geostress measurement methods are limited by deficiencies including the measurable depth, the complexity, and the long duration of operation. To address these problems and achieve the measurement of geostress in deep wells under conditions of complex high pressures and high temperatures, we propose a new measurement method for geostress based on an integrated drilling and optical microscopy system. Its innovative integrated structure eliminates the problems associated with complex procedures and depth limits, and avoids rock creep caused by long delays, significantly improving the accuracy and range of the measurements. It works by using microscopic imaging and direct contact probes to capture the changes of a borehole's cross-sectional outlines before and after stress relief. The resulting images are analyzed with search circles to obtain the positions of probe apices, which can be fitted into ellipses that describe the outlines, and calculate the state of the stress. The validity and accuracy of the method was verified by in-door tests and field applications in the ZK1 borehole. The results show that: (1) the integrated system can be used to measure micrometer-grade deformations; (2) the search-circle approach can accurately obtain the positions of probe apices; and (3) the stress measurement method based on the system is accurate and feasible. 1 INTRODUCTION Geostress, while a universal phenomenon, has been difficult to assess and measure. As we continue our exploration into the deeper crusts, their stresses will play an increasingly greater role in the design and decisionmaking of our deep engineering projects. Work on the theory and practice of stress measurements has led to the creation of new techniques and instruments. The most typical such device was the USBM gauge, developed by the U.S. Bureau of Mines in the 1950s, and the stress-relief method it represents has remained the most commonly used approach [1,2,3]. The stress-relief method is based on the theory of elasticity. It assumes a homogeneous and continuous rock mass, and the same stress-strain function when both loading and unloading stress. It also considers the rock elements being removed as having negligible weight (gravity stress) compared to their load. The measurements are usually made with the following steps: 1) A large borehole is drilled to a predetermined depth close to the measurement position; 2) A smaller hole concentric with the borehole is drilled into the bottom, and three-axis strain gauges (for DOI Acta Geotechnica Slovenica, 2018/1 17. J. Wang et al.: A geostress measurement method based on an integrated drilling and optical microscopic imaging system measuring wall strains) or deformation transducers (for diameter deformations) are installed into the hole; 3) The larger borehole is then extended down, separating the now pipe-like core from the base rock to relieve the stress. This conventional method involves a series of complicated tasks: drilling the large well, leveling the bottom, shaping the pilot hole, drilling the measurement hole, installing transducers, and drilling the relief well. In this process, the drilling rod must be lifted several times from the well, have the drill replaced, and lower back to the bottom. The usefulness of the method in deeper wells is further impeded by the lengths of the transducer cables. For these reasons the utilization of the stressrelief method has often been restricted to the testing of shallow wells [4,5]. Another problem stems from the time delays in measurements. The duration from the installation of the transducers to the end of the stress relief lasts a few days at least, and the test process lasts several tens of minutes by itself. The delays are long enough for rock creep to significantly alter the body's mechanical parameters, deviating from the assumed elastic relationships [6,7]. The resulting stress data are seldom very accurate for this reason [8]. To address these common issues, this paper introduces an integrated drilling and optical microscopic measurement system. It incorporates an outer drill and an inner drill, respectively, for relief and measurement holes. A control module connects and disconnects the drills as needed, replacing the complicated traditional approach. An optical module is used to directly acquire deformation data from the walls. This system can be utilized in deep wells, mitigates the delay-induced errors by real-time measurements, and effectively solves the major technical issues encountered in stress relief. A stress-measurement method based on the system is also introduced. 2 INTEGRATED DRILLING AND OPTICAL MICROSCOPIC MEASUREMENT SYSTEM The system can be divided into the drilling and acquisition sub-systems, the former for drilling the measurement and relief holes, and the latter for acquiring and recording the borehole's diameter deformations. 2.1 Theory The integrated system is still based on the principles of stress relief. It measures the deformations of the rock body before and after the stress relief, and extrapolates the initial stress from them in the following steps: 1) The integrated system drills the measurement hole, and lowers the optical module to the measurement position (see Fig. 1(a)); 2) After the optical module arrives at the position, the control module separates the inner and outer drills, the inner drill stays at rest, while the outer drill resumes drilling (see Fig. 1(b)); 3) The outer drill overtakes the measurement hole, separating the pipe-like core from the rock mass (see Fig. 1(c)), which relieves the stress around the optical module; 4) The optical module detects the minuscule diameter changes, which is used to compute the deformation of the cross-section and its horizontal stress; 5) Measurements are made in three intersecting boreholes to determine the three-dimensional state of stress [9,10,11,12]. c) b) Figure 1. Measurement process: (a) drilling the measurement hole; (b) drilling the relief hole; (c) stress relief. 2.2 Structure of the drilling system The drilling system consists of the outer drill, the control module, and the inner drill, shown in Fig. 2. The inner drill is used to create the measurement hole, and lower the optical microscopic component to the measurement position. After the preset bottom of the measurement hole is reached, the control module will decouple the two drills, and the inner drill will be stationary relative to the surrounding rock body, allowing the acquisition system inside to perform its measurements, while the outer drill continues drilling downwards to create the relief well. As the relief well takes shape, stress can be completely relieved from the core, which is detected by the optical module. At the end of the entire operation, all the system components can be retrieved together. 18. Acta Geotechnica Slovenica, 2018/1 J. Wang et al.: A geostress measurement method based on an integrated drilling and optical microscopic imaging system (a) Outer drill (relief) (b) Control module (decoupling) Figure 2. Structure of the drilling system. (c) Inner drill (measurement) 2.3 Structure of the acquisition system The acquisition system consists of a number of probes, a camera, and an electronic compass, shown in Fig. 3. As the surrounding rock core becomes separated from the larger rock mass, it is relieved of the rock mass' stresses, causing the measurement hole to deform. Each probe is in direct contact with the walls of the hole via a prod protruding from the drill, allowing the deformation to be recorded in real time. Figure 4. Image data from the probes. Compass Camera A \ll\i Figure 3. Structure of the acquisition system. The changing positions of the probes are captured by the camera, consisting of a fixed-focus lens and a microscopic CCD component (see Fig. 4). Since the target area is small, a lens with a small FOV is chosen. Since the CCD image has a greater width than its height, the shorter side is used to determine its physical resolution, which is the physical length of the side divided by its number of pixels, e.g., 576px for a standard capturing component with a resolution of 768*576. The device used in this study captures a 500px-tall image from a 3.5mm-tall area, meaning the image has a physical resolution of 0.007mm. The directional data of the deformation are an important parameters that cannot be accurately recorded by the images. Hence, an electronic compass is used, with its north set to be the same direction as one of the probes. The direction of each probe can then be derived. 2.4 Advantages of the system The proposed system has the following standout features. Firstly, its deployment is vastly simpler than the conventional method. This is accomplished by using a control module to coordinate the movements of the outer and inner drills, integrating the functions of downward drilling, borehole shaping, and measurement in one device. Additionally, the measurement module can start working as soon as the measurement hole is formed, and continue its recording throughout the relief operation, generating more valid data. By eliminating the long delay that comes from replacing drills, the influence of rock creep can be minimized, and the resulting errors can be avoided. This means that the system is more efficient, and less vulnerable to external disturbances than conventional devices. Secondly, the system uses probes in direct contact with the walls, captures their deformation with microscopic imaging, calculates the borehole's deformation after relief from the images, and determines the stress from the deformation. The microscopic imaging allows accuracy at the micrometer level; the probes are made with Acta Geotechnica Slovenica, 2018/1 19. J. Wang et al.: A geostress measurement method based on an integrated drilling and optical microscopic imaging system materials resistant to abrasion, ensuring the system can be operational under unusually harsh conditions. These features give the system high accuracy and adaptability, and present a solution to the limited reach and accuracy in conventional stress-relief methods. 3 ACQUISITION OF DEFORMATION DATA From the image data, it is possible to directly identify the positions of the probes' sharp apices by computer. However, for our required accuracy, this is a process prone to errors. In this study, the images are pre-processed, allowing the probe edges to be detected. The precise coordinates of the apices are then determined as the intersections of these edge lines. 3.1 Pre-processing The inside of a borehole is not a homogeneous environment. The image data can be affected by the noises from external signals and a minor turbidity in the liquid between the camera and the probes, as seen in Fig. 5(a). Since the adaptive median filter (AMF) algorithm can preserve the edges in an image, it is chosen for the de-noise operation. AMF functions by adjusting the filter window based on noise density, using different methods to process noise and signal pixels, maintaining the gray-scale values of the signals, while applying a median filter to the noises. Let the grayscale of the pixel (x,y) be fxy, the current work window be Axy , and the preset largest window permitted be A, in which the minimum, median and maximum gray scales are respectively fmin , fmed and fmax. When fmin Figure 9. Deformation in the borehole coordinate system. Assuming the borehole exis ts in a 3-dimensional, infinite frame, the following relationship can be found between the deformation Ua and the 3-dimensiona 1 stress components: ,d U = <7x(-[(1 +2cos20)(1-v2) + dv2]) + E dvN (6) a, (- [0 - 2 cos 20)(1 - v2) + d v 2]) + az (-—) + (- 4(1 - v) sin 20) E E E where d is the borehole diametea; 9 as the bngle between the direction of the deformed diameter and e'; rotating counterclockwise (CCW) from x' to the direction of the said diameter is positive; E and v are, respectively, the elastic modulus and Poisson's ratio of the rock at the measurement position; and a'x, , a'z , Txb can be converted to the global coordinate system via Equation (5). 4.3 Equation of deformed cross-section After the stress is completed relieved, the borehole cross-section becomes elliptic in shape [16]. Under ideal conditions, the probe apices should also fall onto points on an ellipse, but since the borehole walls are not perfectly smooth, deviations and errors can still be present. Here, we approximate the ellipse with the least-squares method. Now the apices have coordinates (ea,ba) with a=1,2,..., N, Nbeing the total number of probes. The target function can be constructed: F(A,B,C,D,E) ip(xa2 + Axe ■ ba + Bba' + Cxa + Dba + E)2 (7) ail We need to find the parameters (A,B,C,D,E) that allow the apices (xa,ya) to fall onto an ellipse as much as possible. This problem can be represented as obtaining the minimum value of F(A,B,C,D,E). Based on the least-squares theory, this implies the following relationships: =d>F_ dA ~ dB ~ dC ~ dD ~ dE -i 0 (8) The following can be derived from Equation (8): N N N N N Z2 2 3 2 2 xa ■ ba P xa ■ ba P xa ■ ba P xa ■ ba P xa ■ ba a il a il ail a 1 a 1 N N N N N Z3 ^ . ^ 2^3 2 xa ■ ba P ba P xa ■ ba P ba P ba a il ail ail ail ail N N N N N P xa2 ■ ba P xa ■ ba2 P xa2 P xa ■ ba P xa a il ail ail ail ail N N N N N P xa ■ ba2 P ba3 P xa ■ ba P ba P ba ail ail ail ail ail N N N N P xa ■ ba P ba2 P xa P ba N P xa3 ■ ba i ii N Zt t xa ■ ba N P xa2 ■ ba (9) Solving the five-order linear system results in A, B, C, D and E, which are the parameters of an elliptical equation. For the ellipse, the center (x0,bo), major axis a, minor axis b, and angle 9 between the major axis and x can then be determined as the following: , ,2BC - AD 2D - AD. (x°'bo) i } b i 2(ACD - BC2 - D2 + 4BE - A2E) (A2 - 4 B) (B -J A2 + (1 - B)2 +1) (10) 2( ACD - BC2 - D2 + 4BE - A2E) (A2 - 4B)(B + yjA2 + (1 - B)2 +1) 9 i tan-1 a2 - b2B a2B - b2 a 22. Acta Geotechnica Slovenica, 2018/1 J. Wang et al.: A geostress measurement method based on an integrated drilling and optical microscopic imaging system 4.4 Obtaining the stress data The system can measure the stress in a cross-section perpendicular to the borehole axis. But only three equations can be obtained from one borehole, which is not enough to solve for the six components in a 3-dimensional system. To obtain the 3-dimensional stress data, the measurements must be made in three intersecting boreholes. The following uses one borehole as an example to show the process. As in Fig. 10(a), the outline of the cross-sectio n h as changed from the circle before the stress reltef to an ellipse. The deformation can be represented and compared using the outline's elliptical equation. From the origin o, three rays are drawn: the first along the positive x-axis, the second rotated CCW 60 degreeo, and the third rotated CCW 120 degrees. The three rays, respectively, intersect with the pre-relief circle at N1, N2 and N3, and with the post-relief ellipse at M1, M2 and M3, with the coordinates Nj(xni,yni) and M,(xm,,ym,), ¿=1,2,3. Therefore, the distance between each pair of points |MiNi| is: MiNil (xmi -Xni)2 o(ymi -yni)2 (i h 1,2,3) (11) (a) Intersections of rays and the outline pre- and post-deformation Ui*/2 ur/2 x' m u>72 (b) Directions of deformation and stress Figure 10. Deformation and state of stress in a cross-section perpendicular to the borehole. Along the directions of the three rays, three diameter deformations U', U2 and U3' are defined as 2|MNi|, 2\M2N2\ and 2|M3N3|, rotated 60 degrees from each other, as shown in Fig. 10(b). Let fi be the angle between the first deformation UJ and the maximum principal stress a', and a CCW rotation from U' to a' as positive, the state of stress, i.e., the amounts and directions of the stresses a' and a'2 can be computed from the three deformations: a = -[(U ' + U ' + U 3) +—]> 6 d(1 -n2) 1 2 2 U(u -U )2+(U -U1 U +(U3-U U V 1 2 2 3 -'l a ' =-+—2 [(U' +U' +U3 )-—]> 2 6d( 1 -n2) 1 2 3 2 p = e+n 2 U UU - UU )2+(U[- U[ ) 2 +( U-[/(2 (12) where v is the rock's Poisson's ratio and E is its elastic modulus. If we attach an axial strain gauge to the core obtained from the borehole, and apply a confining pressure to it, the ratio of its axial and radial strain is its Poisson's ratio. The elastic modulus can be obtained by using a borehole deformation gauge to measure its confining pressure-deformation curve, using the following equation: 4PnrR2 E = - 0 U(R2 - r2) (13) where P0 is the confining pressure; U is the average fine deformation caused by the confining pressure; R and r are, respectively, the outer and inner radii of the core. If the axis of the borehole is along the same direction as a principal stress, the amount of which is already known, e.g., the gravity stress in a perfectly vertical well, measuring one borehole will be enough to determine the 3-dimensional state of stress at the measurement point. 5 CASE STUDY_ We used the system to measure the stresses in the bottom outlets of a certain hydropower station. For the borehole ZK1 with diameter 168 mm, the system performed the measurement operation along a depth of 150 m, using a 91-mm outer drill, and a 40-mm inner drill. The cross-sectional stress data were obtained at 15 points, each separated by a depth interval of 10 m. As an example, data from the first three points are analyzed below. Acta Geotechnica Slovenica, 2018/1 23. J. Wang et al.: A geostress measurement method based on an integrated drilling and optical microscopic imaging system 5.1 Solving for stress First, the system conducted data acquisition from depths 10, 20 and 30 m, generating images of probes before and after the stress relief. The images were processed by computer and stored as image data with directional information. The data then went under pre-processing, edge detection, and determination of probe apices. With the positions of the probe apices known, they were fitted using least squares into circles and ellipses representing the statuses of the borehole before and after the stress relief. The circles or ellipses were described by equations, as shown in Fig. 11, Cartesian systems were established, and rays were drawn to calculate the three deformations U', U2 and U3' in each instance. * Niwi \ \ / \ \ / \ \ / K 0 MV3 t**. (a) Cross-section at 10m (b) Cross-section at 20m V —-, /NN f IV/Î, V Figure 11. Cross-sections at depths 10, 20 and 30m, pre- and post-stress relief. The deformations U', U2 and U3 can then be used to determine the stresses in each cross-section using Equation (12), while the elliptical equation in each instance can be used to determine the angle 9 between the major axis and the x-axis, thus obtaining the directions of two principal stresses for each measurement point. The results are shown in Table 2. 5.2 Analyses of the result We conducted tests to validate the feasibility and reliability of the integrated system and the method. First, a three-axis compressor was used to perform a series of compression tests. For each core, confining pressures were applied that have the same amount and direction as the cross-sectional principal stresses above. Standard instruments were then used to measure the diameters of each hole before and after applying the pressures. The instrument is an inside micrometer with a precision of 0.001mm. The deformations caused by applying the known stresses were compared with the integrated system's results (see Fig. 12). (a) Deformations of the 10 m core Since the calculation requires knowing the rocks' Poisson ratios and elastic moduli, a RMT-150C rock testing machine was used to measure these parameters, with results shown in Table 1. Table 1. Poissons ratios and elastic moduli for cores at depths 10, 20 and 30 m. (b) Deformations of the 20 m core Core location 10 m 20 m 30 m Elastic modulus E/GPa 25.1 37.5 34.0 Poisson's ratio v 0.27 0.27 0.27 (c) Deformations of the 30 m core (Blue: borehole measurements; Red: compression test results) Figure 12. Comparison of deformations in borehole measurements and compression tests. Table 2. Deformations and principal stresses along the cross-sections. Location Deformation Principal stresses along the cross-section U'/pm U'/pm U /pm amount direction (¡2 amount o'2 direction 10 m 282 318 282 9.1MPa N29.88°E 1.0MPa S60.12°E 20 m 264 328 284 14.3MPa N55.34°E 1.5MPa S44.66°E 30 m 287 280 326 10.7MPa N39.67°W 1.3MPa S50.33°W 24. Acta Geotechnica Slovenica, 2018/1 J. Wang et al.: A geostress measurement method based on an integrated drilling and optical microscopic imaging system The greatest variance between the two is 1.06% at U3' in Fig. 12(a), 2.44% at U2 in Fig. 12(b), and 2.44% at U' in Fig. 12(c). It can be seen that the general tendencies of the deformations are similar between the two groups of results, with errors below 2.5%. Therefore, the system's measurements can be deemed valid. The inner surfaces of a borehole might crack or peel off when the stress reaches a certain level. Such occurrences can also be used to determine the direction of the stress [17,18,19]. In our validation of the system, we used it to observe this borehole, and found the existence of cracks at 20 m (see Fig. 13(a)), and peeled off surfaces at 49 m and 58 m (see Figs. 13(b) and (c)). Their directions were consistent with calculations of the system, which were, respectively, along the west-east and north-south axes. For the measurements of ZK1, the borehole was along the vertical direction. If we consider the geostress to be a principal stress, one measurement is enough to determine a point's 3-dimensional state of stress. Hence, the measurements over the full lengths of the borehole can be used to describe the vertical distribution of the stresses, which is depicted as rose diagrams in Table 3. The circumferential direction represents the directions of the stresses, and rotates clockwise from true north. The radial direction represents the amount of stress, with weaker stresses closer to the center, and stronger ones farther away from the center. It can be seen in Table 3 that ZKl's maximum principal stress largely falls in the 10 MPa range, existing along the west-east direction, while the minimum principal Table 3. Summary of stress distribution in ZK1. Principal Stress m Rose Diagram Principal Stress <72 Rose Diagram Principal stress-depth correlations Blue: <71; Red: <72 (Ver.: MPa; Hor.: 10m) Acta Geotechnica Slovenica, 2018/1 25. J. Wang et al.: A geostress measurement method based on an integrated drilling and optical microscopic imaging system stress is around 1.3 MPa along the north-south direction. There is no linear correlation between the stress and depth; the stresses are affected by multiple factors, and vary greatly between the depths. The distribution of stress is at least partially stochastic. It can only be described using a unified equation when the distribution is normal or log-normal [20]. For this survey, linear interpolation is used to describe the stress within each 10 m segment between measurement points, and reduce the complexity of stochasticity. The stress at each depth value is assumed to be a linear function of the two measurement points before and after it. The results show that the stress varies considerably with depth. We chose the maximum values as an intuitive way to describe the borehole's stress characteristics. The maximum values of the principal stresses a1 and a2 are used as the reference values of the stress in this borehole, while the reference directions are obtained by observing the rose diagram. The measurements of ZK1 indicate the system can effectively acquire the borehole deformations along the well, avoiding changes caused by rock creep during the long periods of waiting, which leads to more accurate results. The system also enables more efficient measurements of greater depths. The method of stress determination based on the system is a good approach to describe the movements of the probes, which can be used to calculate the deformations and principal stresses. The resulting cross-sectional stresses from the depths can be used to provide an intuitive depiction of the borehole's stress variations using rose and contoured diagrams. 6 CONCLUSION This paper proposes a measurement method for geostress based on an integrated drilling and optical microscopy system. It employs microscopic imaging and direct contact probes to capture the changes of a borehole's cross-sections before and after the stress relief. The resulting images are analyzed with search circles to obtain the positions of the probe apices, which can be fitted into ellipses that describe the cross-sectional outlines, and calculate the states of stress. The validity and accuracy of the method was verified by tests and field application in the ZK1 borehole. The results show that: (1) the integrated system avoids inaccuracies caused by time delays; (2) the integrated system can be used to measure micrometer-grade deformations; (3) the search circle approach can accurately obtain the positions of probe apices; (4) the stress-measurement method based on the system is accurate and feasible; (5) the method can be used for stress measurements in deep positions and complex environments. REFERENCES [1] Liu, Y. F. 2000. Rock Body Geostress and Construction Engineering [in Chinese], Hubei Science & Technology, Wuhan. [2] Ge, X.R., Hou, M.X. 2012. Principle of in-situ 3D rock stress measurement with borehole wall stress relief method and its preliminary applications to determination of in-situ rock stress orientation and magnitude in Jinping hydropower station. Sci China Tech Sci. 55, 939-949. DOI: 10.1007/s11431-011-4680-x [3] Worotnicki, G., Walton, R. 1976. Triaxial "hollow inclusion" gauges for determination of rock stresses in-situ. Investigation of Stress in Rock - Advances in Stress Measuremnent, Proc, Int, Symp, Sydney,1-8 (supplement). [4] Nur, A., Slmmons, G. 1969. Stress-induced velocity anisotropy in rock: an experimental study. J. of Geophys. Res. 74(27), 667~6674. DOI: 10.1029/ JB074i027p06667 [5] AI Kai, Han Xiao-yu, Li Yong-song 2006. Regression Analysis of In-situ Stress for Wudongde Hydropower Station, Jinsha River [in Chinese]. Chinese Journal of Underground Space and Engneering 2, (26), 930-933. DOI: 10.3969/j.issn.1673-0836.2006.06.011 [6] Xu, Z.L. 2006. Elastic Mechanics [in Chinese], China Higher Education, Beijing. [7] Luo, C.H., Wang, W. 2003. Refined theory of generalized plane-stress problems in elasticity. Acta Mechanica Solida Sinica 24(02), 192-196. DOI: 10.3969/j.issn.0254-7805.2003.02.009 [8] Yoshikawa, S., Mogi, K. 1981. A new method for estimation of the crustal stress from cored rock samples: laboratory study in the case of uniaxial compression. Tectonophysics 74, 323-339. DOI: 10.1016/0040-1951(81)90196-7 [9] Cai, M.F., Qiao, L., Li, H.B. 1995. Theory and Techniques of Geostress Measurement [in Chinese]. China Science, Beijing. [10] Amadei, B., Stephansson, O. 1997. Rock stress and its measurement, Chapman and Hall, London. [11] Duncan Fama, M.E., Pender, M.J. 1980. Analysis of the hollow inclusion technique for measuring in-situ rock stress. Int. J. Rock Mech. Min. Sci. And Geomech. Abstr. 17, 137-146. DOI: 10.1016/0148-9062(80)91360-1 [12] Wang, L.J., Pan, L.Z. 1991. Geostress Measurement and Its Application in Engineering [in Chinese]. China Geology, Beijing. [13] Zhang, Z., Xu, C. 2014. Digital Image Processing and Machine Vision, China Posts & Telecom, Beijing. 26. Acta Geotechnica Slovenica, 2018/1 J. Wang et al.: A geostress measurement method based on an integrated drilling and optical microscopic imaging system [14] Brudy, M., Zoback, M.D., Fuchs, K. et al. 1997. Estimation of the complete stress tensor to 8 km depth in the KTB scientific drill holes: implications for crustal strength. Journal of Geophysical Research 102(B8), 18453-18475. DOI: 10.1029/96JB02942 [15] Sun, W.C., Min, H., Wang, C.Y. 2008. Three-dimensional geostress measurement and geomechanical analysis. Chinese Journal of Rock Mechanics and Engineering 27(z2), 3778-3784. DOI: 10.3321/j.issn:1000-6915.2008.z2.070 [16] Wang, C.Y., Han, Z.Q., Wang, J.C., Wang, Y.T. 2016. Study of borehole geometric shape features under plane stress state, Chinese Journal of Rock Mechanics and Engineering 35(z1), 2836-2846. DOI: 10.13722/j.cnki.jrme.2015.1045 [17] Quan Jiang, Xia-ting Feng, Jing Chen, Ke Huang, Ya-li Jiang 2013. Estimating in-situ rock stress from spalling veins: A case study. Engineering Geology 152(1), 38-47. DOI: 10.1016/j. enggeo.2012.10.010 [18] Jiang, Q., Feng, X.T., Xiang, T.B., Su, G.S. 2010. Rockburst characteristics and numerical simulation based on a new energy index: a case study of a tunnel at 2500 m depth. Bulletin of Engineering Geology and the Environment 69, 381-388. DOI: 10.1007/s10064-010-0275-1 [19] Quan Jiang, Xia-ting Feng, Yilin Fan, Qixiang Fan, Guofeng Liu, Shufeng Pei 2017. In situ experimental investigation of basalt spalling. Tunnelling and Underground Space Technology 68, 82-94. DOI: 10.1016/j.tust.2017.05.020 [20] Jiang Q., Zhong S., Cui J., Feng X.T., Song, L. 2016. Statistical Characterization of the Mechanical Parameters of Intact Rock Under Triaxial Compression_ An Experimental Proof of the Jinping Marble. Rock Mechanics & Rock Engineering 49 (12), 4631-4646. DOI: 10.1007/s00603-016-1054-5 Acta Geotechnica Slovenica, 2018/1 27.