Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
AU2011250753B2 - Method and apparatus for measuring in-situ stress of rock using thermal crack - Google Patents
[go: Go Back, main page]

AU2011250753B2 - Method and apparatus for measuring in-situ stress of rock using thermal crack - Google Patents

Method and apparatus for measuring in-situ stress of rock using thermal crack Download PDF

Info

Publication number
AU2011250753B2
AU2011250753B2 AU2011250753A AU2011250753A AU2011250753B2 AU 2011250753 B2 AU2011250753 B2 AU 2011250753B2 AU 2011250753 A AU2011250753 A AU 2011250753A AU 2011250753 A AU2011250753 A AU 2011250753A AU 2011250753 B2 AU2011250753 B2 AU 2011250753B2
Authority
AU
Australia
Prior art keywords
coolant
crack
borehole
stress
wall
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Ceased
Application number
AU2011250753A
Other versions
AU2011250753A1 (en
Inventor
Bung-Hee Choi
Loui Porathur John
Chang-Ha Ryu
Dong-Woo Ryu
Joong-Ho Synn
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Korea Institute of Geoscience and Mineral Resources KIGAM
Original Assignee
Korea Institute of Geoscience and Mineral Resources KIGAM
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from AU2007314874A external-priority patent/AU2007314874B2/en
Application filed by Korea Institute of Geoscience and Mineral Resources KIGAM filed Critical Korea Institute of Geoscience and Mineral Resources KIGAM
Priority to AU2011250753A priority Critical patent/AU2011250753B2/en
Publication of AU2011250753A1 publication Critical patent/AU2011250753A1/en
Application granted granted Critical
Publication of AU2011250753B2 publication Critical patent/AU2011250753B2/en
Ceased legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

Abstract Disclosed is a method and apparatus for measuring in-situ stress in rock using a thermal crack. The method involves forming a borehole, cooling a wall of the borehole, 5 applying tensile thermal stress, forming a crack in the borehole wall, and measuring temperature and cracking point. Afterwards, the borehole wall is heated to close the formed crack, the borehole wall is cooled again to re-open the crack, and temperature is measured when the crack is re-opened. The in-situ stress of the rock is calculated using a first cracking temperature at which the crack is formed and a second cracking temperature at which the 10 crack is re-opened. Further, the apparatus cools, heats and re-cools the borehole wall, thereby measuring the first cracking temperature, the second cracking temperature, and the cracking point.

Description

AUSTRALIA Patents Act COMPLETE SPECIFICATION (ORIGINAL) Class Int. Class Application Number: Lodged: Complete Specification Lodged: Accepted: Published: Priority Related Art: Name of Applicant: Korea Institute of Geoscience and Mineral Resources Actual Inventor(s): Bung-Hee Choi, Chang-Ha Ryu, Loui Porathur John, Joong-Ho Synn, Dong-Woo Ryu Address for Service and Correspondence: PHILLIPS ORMONDE FITZPATRICK Patent and Trade Mark Attorneys 367 Collins Street Melbourne 3000 AUSTRALIA Invention Title: METHOD AND APPARATUS FOR MEASURING IN-SITU STRESS OF ROCK USING THERMAL CRACK Our Ref: 928593 POF Code: 487366/489284 The following statement is a full description of this invention, including the best method of performing it known to applicant(s): 600eq [DESCRIPTION] [Invention Title] METHOD AND APPARATUS FOR MEASURING IN-SITU STRESS OF ROCK USING THERMAL CRACK 5 The present application is a divisional application from Australian Patent Application No. 2007314874, the entire disclosure of which is incorporated herein by reference. [Technical Field] The present invention relates to a method and apparatus for measuring in-situ stress in rock, and more particularly, to a method and apparatus for measuring in-situ stress in rock, in 10 which the in-situ stress in rock is measured by applying heat to the rock using a cryogenic coolant to thereby generate cracks. [Background Art] The term "in-situ stress in rock" refers to stress that exist in the interior of rock, including gravitational stress, tectonic stress, and residual stress. Here, gravitational stress 15 indicates stress generated by the rock's own weight. Tectonic stress indicates stress generated by movement of the earth's crust. Residual stress indicates stress remaining after removal of its original cause, such as expansion or heating of the rock or a past surface load since removed by surface erosion. In the design or safety analysis of a large-scale structure in rock such as a tunnel or an 20 oil storage tank, calculation of in-situ stress in the rock is proving to be very important. This is because no underground structure can be designed and constructed in a stable and economical manner inside rock until the direction and magnitude of stresses acting on the rock have been accurately measured. For example, in tunnel construction, if pressure applied to surrounding rock is isotropic, the tunnel cross section is generally circular. On condition that there is strong 25 transverse pressure on the rock such as by a surface load, though the tunnel cross section is elliptical, the tunnel keep safe from collapse. For instance, when excavating a tunnel without exactly measuring the in-situ stress in the rock, the rock may become overstressed due to stress concentration on an excavated surface and may collapse or become unstable due to expansion of existing cracks. Thus, in 30 order to install a structure in rock, an accurate measurement of in-situ stress in rock is required.
]A
Methods for measuring in-situ stress in rock include hydraulic fracturing (Fairhurst, 1964), stress opening such as overcoring (Leeman and Hayes, 1966; Merrill, 1967; USPTO 4491022), an lB indirect method such as acoustic emission and so forth. Among these methods, hydraulic fracturing and overcoring are frequently used. A conventional hydraulic fiacturing system is illustrated in FIG. 1. Refening to FIG. 1, a borehole h is forced in the rock in which an underground structure will be constructed, and then packers I arc 5 installed at upper and lower ends of a test section a, thereby sealing the test section. Afterwards, a high-pressure pump (not shown) applies hydraulic pressure, thereby injecting fluid into test section a through a pipe p. When hydraulic pressure continuously increases to reach 'initial fracturing pressure', cracks are formed in a borehole wall w. When fluid is continuously injected, the cracks are gradually widened and pressure is lowered. When the injection is stopped, the cracks am closed and 10 the pressure reaches 'crack closing pressure'. Here, since the crack closure pressure is equal to or slightly greater than the press required to maintain the cracks, it represents the minimum principal stress, which acts in a direction perpendicular to the cracked surface. Further, while the above mentioned processes are repeated (second cycle), the cracks are re-opened (crack re-opening pressure). Thus, after crack re-opening pressure is measured, the maximum horizontal principal is stress can be calculated using the measured pressure. Among the methods for measuring in-situ stress in rock, hydraulic fracturing has an advantage in that it can be applied to deep underground rock as long as the borehole can be drilled, but has a disadvantage in that it is not easily applied to specific types of rock, for instance, sedimentary rock with a stratified structure. Further, when pressure is applied using packers 1, the 20 borehole may collapse or the packers may jam. This complicates withdrawal of equipment, especially at greater depths. Further, since the size of the equipment must be increased in order to apply high hydraulic pressure, there are difficulties in using the conventional hydraulic fracturing system. As described above, hydraulic featuring has a limitation in that only two parts of principal 25 stress can be measured. In other words, hydraulic fracturing has a limitation in that, since the vertical stress is set to the surface load (density * mass * height), the maximum and minimum principal stress can be measured in only a horizontal direction, and thus in-situ stress cannot be accurately measured. Furthermore, many researchers who have extensively studied hydraulic fracturing cast a doubt on the accuracy of crack re-opening pressure as well as the variation in pore hydraulic pressure in cracks (Ito 2 et al., 2001). Unlike hydraulic fracturing, overcoring has an advantage in that six stress components existing on three dimensions can be provided. However, since overcoring is based on the measurement of strain, it requires elaborate dual coring work. Hence, when applied, 5 overcoring is problematic in that it is highly complex as well as restricted to the depth of the borehole. A reference herein to a patent document or other matter which is given as prior art is not to be taken as an admission that that document or matter was known or that the information it contains was part of the common general knowledge as at the priority date of 10 any of the claims. [Disclosure] [Technical Problem] Accordingly, the present invention has been made in an effort to ameliorate the problems occurring in the related art, and an aspect of the present invention is to provide a 15 method for measuring in-situ stress in rock, which can be applied to rock at a great depth to measure in-situ stress by directly forming cracks in rock by means of heat, thereby enabling precise and easy measurement without elaborate equipment. Another aspect of the present invention is to provide an apparatus for measuring in-situ stress in rock, which can be applied to rock at a great depth to measure in-situ stress by 20 directly forming cracks in rock by means of heat, thereby enabling precise and easy measurement without elaborate equipment. [Technical Solution] In order to achieve the above aspects, according to one aspect of the present invention, there is provided an aspect of the present invention, there is provided an apparatus for 25 measuring in-situ stress in rock using a thermal crack by drilling a borehole into the ground and forming a crack by heating the wall of the borehole. The apparatus comprises: a coolant container, in which an annular containing space where coolant can be contained, and an inlet for the coolant to flow into and out of the containing space are formed; a close-contact means for bringing an outer surface of the coolant container into close contact with the borehole wall; 30 at least one temperature sensor which is installed on an outer surface of the coolant container so as to measure temperature of the borehole wall; and a crack detecting means for detecting cracks formed in the borehole wall by heat transmission between the coolant contained in the coolant container and the borehole wall; wherein the crack detecting means detects a burst 3 sound occurring when a crack is formed from the rock with the borehole wall, and comprises a plurality of acoustic emission sensors attached to an outer circumference of the coolant container. 5 [Advantageous Effects] The apparatus for measuring in-situ stress in rock using a thermal crack according to the present invention can measure the in-situ stress in rock at a great depth, and can be easily used due to simple configuration thereof. [Description of Drawings] 10 The above and other features and advantages of the present invention will become more apparent after a reading of the following detailed description when taken in conjunction with the drawings, in which: FIG. 1 is a schematic view for explaining an in-situ stress measuring method and apparatus using a conventional hydraulic fracturing method; 15 FIG.2 is a flowchart illustrating a method for measuring in-situ stress in rock using a thermal crack according to an exemplary embodiment of the present invention; FIG. 3 is a schematic view illustrating a borehole for explaining the in-situ stress measuring 4 method of FIG. 2; FIG. 4 is a schematic view illustrating the configuration of an apparatus for measuring in-situ stress in rock using a thermal crack according to an exemplary embodiment of the present invention, wherein FIG. 4A shows the state in which a coolant container is expanded, and FIG. 4B shows the 5 state in which a coolant container is contracted; FIG. 5 is a schematic cross-sectional view taken along line V-V of FIG. 4A; FIG. 6 is an exploded perspective view illustrating the coolant container, the main part of the in-situ stress measuring appartus illustrated in FIG. 4; FIG. 7 is a schematic view for explaining the detection of cracks by acoustic emission 10 sensors; FIG. 8 illustrates a resultant screen in which two-dimensional (plane) modeling is performed on the formation of a thermal crack from a borehole wall in the first cooling step using a numerical analysis program of "FLAC3D", wherein FIGS. 8A through 8J are resultant screens showing a thermal crack growing with the lapse of time; 15 FIG. 9 is a graph comparing the results of the modeling of FIG. 8 with theoretical values; FIG. 10 illustates a resultant screen in which two-dimensional (plane) modeling is performed on the formation of a thermal crack from a borehole wall in the second cooling step using a numerical analysis program of "FLAC3D", wherein FIGS. IOA through 1OH are resultant screens showing a thermal crack growing with the lapse of time; 20 FIG. 11 is a graph comparing the results ofthe modeling of FIG. 10 with theoretical values; FIG. 12 illustrates a resultant screen in which three-dimensional modeling is performed on the formation of a thermal crack from a borehole wall in the first cooling step using a numerical analysis program of "FLAC3D", wherein FIGS. 12A through 12H are resultant screens showing a thermal crack growing with the lapse of time; 25 FIG. 13 is a graph comparing the results of the modeling of FIG. 12 with theoretical values; FIG. 14 illustrates a resultant screen in which three-dimensional modeling is performed on the formation of a thermal crack from a borehole wall in the second cooling step using a numerical analysis program of "FLAC3D", wherein FIGS. 14A through 14F are resultant screens showing a thermal crack growing with the lapse of time; and 5 FIG.15 is a graph comparing the results ofthe modeling of FIG. 14 with theoretical values; [Best Model In order to achieve the above objects, according to one aspect of the present invention, ther is provided a method for measuring in-situ stress in rock using a thermal crack The method 5 comprises: a borehole forming step of forming a borehole in a target rock for measuring in-situ stress; a first cooling step of cooling a wall of the borehole, applying tensile thermal stress to the borehole wall, forming and growing the crack in the borehole wall, and measuring a first cracking temperature of the borehole wall when the crack occurs and a cracking point at which the crack is formed; a heating step of heating the borehole wall cooled in the first cooling step and closing the formed crack; 10 a second cooling step of cooling the borehole wall again, applying tensile terminal stss to the borehole wall, re-opening the crack closed in the heating step, and measuring a second cracking temperature of the borehole wall when the crack is re-opened; and a calculating step of calculating the in-situ stress of the rock using the first cracking temperatume of the borehole wall and the cracking point measured in the first cooling step, and using the second cracking temperature of the borehole 15 wall measured in the second cooling step. According to one aspect of the present invention, in-situ stress includes maximum and minimum horizontal principal stresses, which act in directions perpendicular to each other on a plane perpendicular to an axis of the borehole; and the maximum horizontal principal stress (ai) and the minimum horizontal principal stress (a 2 ) can be obtained using the following first equation. 20 <First Equation> (1- 2 cos26) a +(1+2cos2O)a 2+ a,=C( ta-t) where a, is the maximum horizontal principal stress, o2 is the minimum horizontal principal stress, 0 is the rotating angle measured in a counterclockwise direction from a point, on which the maximum horizontal principal stress acts, to the cracking point centered around a central point of the 25 borehole on the plane perpendicular to the axis of the borehole, atis the tensile strength of the rock, t. is the temperature of the rock before cooling, and t is one of the first and second cracking temperatures, and where C is the bi-axial thermo-elastic constant of the rock in which the borehole wall is 6 formed and is expressed by C = Ea/(l-v), where E is the elastic coefficient (Young's module) of the rock, a is the linear thermal expansion coefficient of the rock, and v is the Poisson's ratio of the rock. According to one aspect of the present invention, in-situ stress includes a vertical stress, which acts in directions perpendicular to directions of the maximum and minimum horizontal 5 principal stresses; and the vertical stress is obtained using the following second equation. <Second Equation>
C
3 =C( t,,- t 3 )- at where a 3 is the vertical stress, and t 3 is the temperature of the rock when a transverse crack perpendicular to an axial direction of the borehole is formed. 10 According to one aspect of the present invention, the tensile strength (or) of the rock is set to 0 when the first equation is established using the second cracking temperature measured by re opening the crack in the second cooling step. According to one aspect of the present invention, t3 is set by averaging the cracking start and end points when a transverse crack is formed in a circular shape along the borehole wall and when the 15 temperatures of the borehole wall are different from each other at the cracking start and end points. According to one aspect of the present invention, maximum and minimum horizontal principal stresses are decided by: forming a plurality of cracks at different points in the borehole wall, creating a plurality of first equations, and calculating a plurality of solutions of maximum and minimum horizontal principal stresses; and performing a least square method using the plurality of 20 solutions of maximum and minimum horizontal principal stresses. According to one aspect of the present invention, the heating step comprises of introducing external air into the borehole to heat the borehole. According to another aspect of the present invention, there is provided an apparatus for measuring in-situ stress in rock using a thermal crack by drilling a borehole into the ground and 25 fonning a crack by heating the wall of the borehole. The apparatus comprises: a coolant container, which has an annular containing space where coolant can be contained, and an inlet for coolant to flow into and out of the containing space; a close-contact means for bringing an outer surface of the coolant container into close contact with the borehole wall; at least one temperature sensor which is 7 installed on an outer surface of the coolant container so as to measure temperature of the borehole wall; and a crack detecting means for detecting cracks formed in the borehole wall by heat transmission between the coolant contained in the coolant container and the borehole wall. According to another aspect of the present invention, the coolant container is provided 5 therein with a filling space, which is surrounded by the inner wall, the upper cover, and the lower cover of the coolant container so as to fill and discharge a fluid; The close-contact means is a pump connected with the filling space so as to fill and discharge the fluid into and from the filling space; and the outer and inner walls of the coolant container are made of an elastic material which can expand and contract as the fluid is filled into and discharged from the filling space. 10 According to another aspect of the present invention, the apparatus further comprises a coolant injecting means for injecting coolant into the containing space. The coolant injecting means comprises: a cylinder, which is coupled to the coolant container such that a coolant chamber, in which the coolant is contained, is formed between the coolant container and the cylinder; a piston, which has a piston head inserted into the cylinder, and a piston rod having a bar shape and fixed to the 15 piston head, which is installed so as to be able to be reciprocated in the cyinder, and which presses the coolant toward the coolant container; a driver, which reciprocates the piston; a valve, which is installed on the coolant container so as to mutually communicate and shut the coolant chamber and the containing space. According to another aspect of the present invention, the crack detecting means detects a 20 burst sound which occurs when a crack is formed in the rock with the borehole wall, and includes at least one acoustic emission sensor attached to the outer circumference of the coolant container. [Mode for Invention] Hereinafter, an apparatus for measuring in-situ stress in rock using a thermal crack according to an exemplary embodiment of the present invention will be described with reference to the 25 accompanying drawings, and than a measuring method using the apparatus will be described in greater detail. FIG. 4 is a schematic view illustrating an apparatus for measuring in-situ stress in rock using a thermal crack according to an exemplary embodiment of the present invention, wherein FIG. 4A shows the state in which a coolant container is expanded, and FIG. 4B shows the state in which a 8 coolant container is contracted. FIG. 5 is a schematic cross-sectional view taken along line V-V of FIG. 4A. FIG. 6 is an exploded perspective view illustrating the coolant container, the main part of the in-situ stress measuring apparatus illustrated in FIG. 4. Referring to FIGS. 4 through 6, an appamtus 100 for measuring in-situ stress in rock using a 5 thermal crack according to an exemplary embodiment of the present invention bores a borehole h into the ground in which an underground structure such as a tunnel will be installed, and then applies cryogenic heat to a wall w of the borehole so as to generate cracks, thereby measuring the in-situ stress in the rock. The measuring apparatus comprises a coolant container 10, a close-contact means, a temperature sensor 40, a coolant injecting means, and a crack detecting means. 10 The coolant container 10 is a case that contains coolant r for applying the cryogenic heat to the borehole wall w. In this embodiment, the coolant container 10 has the structure of a double wall, and thus forms an annular containing space 15. To this end, the coolant container 10 includes an inner wall 11, an outer wall 12, an upper cover 13 and a lower cover 14. The inner wall II is formed in an approximately circular shape, and is disposed at a preset is height. The outer wall 12 is formed in a circular shape so as to correspond to a shape of borehole h, and is disposed outside the inner wall I I while spaced apart from an outer circumference of the inner wall 11. he outer wall 12 has a height equal to that of the inner wall 11. The inner wall l and the outer wall 12 are made of elastic material, which can be contracted and expanded. Further, the inner wall 11 and the outer wall 12 are made of high heat resistant material so as to be able to withstand 20 when coming into contact with a low-temperature coolant The upper cover 13 and the lower cover 14 are coupled to upper and lower sides of the inner wall 11 and the outer wall 12, respectively. In this case, the annular space, i.e., the containing space 15, surrounded by the inner wall 11, the outer wall 12, the upper cover 13, and the lower cover 14 is defined. The containing space 15 is a space in which coolant r is injected and contained. In this 25 embodiment, liquefied nitrogen (LN 2 ) is used as the coolant. The lower cover 14 is provided with an inlet 18a so as to allow coolant r to flow into and out of the containing space 15. To this end, the inlet 18a passes through the lower cover 14. The inlet 18a allows the containing space 15 to mutually communicate with a coolant chamber c, which will be -described below. Further, an inlet valve 18b for shutting and opening the inlet 18a is installed in the 9 inlet 1ga. Further, the upper cover 13 is provided with an outlet 19a for the coolant to flow out of the containing space 15. To this end, the outlet 19a passes through the upper cover 13. A pressure decreasing valve l9b for shutting and opening the outlet 19a is installed in the outlet 19a. Meanwhile, as described above, when the upper cover 13 and the lower cover 14 are coupled 5 to the upper and lower sides of the inner wall I I and the outer wall 12 respectively, a filling space 16 surrounded by the inside of the coolant container 10, that is, the inner wall 11, the upper cover 13 and the lower cover 14 is defined. The close-contact means functions to carry out smooth heat transmission between the coolant r and the borehole wall w by bringing the outer wall 12 of the coolant container 10 into close contact 10 with the borehole wall w in the state where the coolant r is contained in the coolant container 10. In this embodiment, as the close-contact means, a pump is used that can bring the inner wall 11 and the outer wall 12 of the coolant container 10, which are made of elastic material, into close contact with the borehole wall w by filling a fluid into the falling space 16. The pump can be varied depending on the fluid being filled into the filling space 16. Since the filing space 16 is filled with gas in this 15 embodiment, a pneumatic pump (not shown) is used. More specifically, the upper cover 13 is provided with a gas injecting hole 17a, which is connected with the filling space 16, at the center thereof. A pipe 17b, which is connected with the pneumatic pump, is installed through the gas injecting hole 17a, so that the pneumatic pump is mutually connected with the filling space 16. The pneumatic pump is operated be t lo fill and discharge the gas into and from the filling space 16, 20 In other words, when the filling space 16 is filled with the gas, the outer wall 12 of the coolant container 10 expands to come into close contact with the borehole wall w. In contrast, when the gas is discharged, the outer wall 12 of the coolant container 10 contracts and separates from the borehole wall w. The coolant injecting means functions to inject the coolant r into the containing space 15 of 25 the coolant container 10, and comprises a cylinder 21, a piston24, a driver, and the inlet valve 18b. The cylinder 21 has a cylindrical shape, and is fixed to the lower portion of the coolant container 10, particularly to the lower cover 14. When the cylinder 21 is coupled to the coolant container 10, the coolant chamber c is defined between the lower cover 14 of the coolant container 10 and the cylinder 21. The coolant r is temporally contained in the coolant chamber c before it is 10 injected into the coolant container 10. The piston 24 is installed on the cylinder 21, and includes a piston head 22 and a piston rod 23. In other words, the piston head 22, which has a disc shape and a diameter substantially equal to that of the cylinder 21, is fitted into the cylinder 21. The piston rod 23 having a bar shape is fixed to a lower end of the piston head 22, and is coupled to a driver, i.e., a 5 motor 30, which will be described below. This piston rod 23 is provided with a male thread 23a on an outer circumference thereof. he motor 30 functions to reciprocate the piston 24, makes use of a known reversible motor, and includes a stator (not shown) and a rotor 31. The rotor 31 is a member that is rotated when power is applied like the known motor, is formed in a hollow shape, and is provided with a female thread 32 10 in an inner circumference thereof. The piston rod 23 is fitted inside the rotor 31 of the motor 30. As described above, because the male thread 23a is formed on the outer circumference of the piston rod 23, the male thread 23a of the piston rod 23 is screwed together with the female thread 32 of the rotor 31. Thus, when the rotor 31 of the motor 30 is rotated in a forward or reverse direction, the piston rod 23 moves in an upward or downward direction. 15 When the rotor 31 of the motor is rotated in a forward direction, the piston rod 23 moves in an upward direction, and the piston head 22 pushes the coolant r in the coolant chamber c towards the coolant container 10, and thereby the coolant r is injected into the containing space 15 through the inlet valve 18b. In contrast, when the rotor 31 of the motor is rotated in a reverse direction, the piston rod 23 moves in a downward direction, and the coolant r in the containing space 15 is pushed into the 20 coolant chamber c through the inlet valve l8b. Further, the piston rod 23 extends through the piston head 22. In other words, the piston rod 23 extends through the coolant chamber c as well as the filling space 16 and the upper cover 13 of the coolant container 10. Thus, the piston rod 23 protrudes upwards from the coolant container 10 at an upper end thereof, and protrudes downwards from the rotor 31 of the motor 30 at a lower end thereof 25 Further, the piston rod 23 is provided therein with a plurality of holes 26,27 and 28 in a longitudinal direction (i.e. in an axial direction of the borehole). At least one of the holes 26, 27 and 28 acts as an air injecting hole 26, which introduces air in a downward direction, passes through a portion between the upper and lower ends of the piston rod 23 and reaches a lower part of the borehole h. The other holes are connecting holes 27 and 28, in which electric cables (not shown) connected with the motor l 1 or various sensors, which will be described below, are installed. Air, which is injected through an upper end of the air injecting hole 26 and then is discharged in a downward direction can function to heat the borehole wall w cooled by the coolant r while moving along the borehole h in an upward direction. In the present invention, as will be described below, the process of fonning cracks in the 5 rock using the coolant in order to measure in-situ stress in the rock is repeatedly perfonned. Thus, the borehole wall w cooled by the coolant r must be restored to room-temperature. Here, thawing of the borehole wall w can be accelerated by the air injecting. The temperature sensor 40 functions to measure temperature of the borehole wall w. As will be described below, in order to measure in-situ stress in rock, temperature of the borehole wall w 0 o must be accurately obtained at a point of time when cracks are generated in the borehole wall w due to the cooling of the coolant r. To this end, a plurality of temperature sensors 40 are installed on the outer circumference, for instance, the outer wall 12, of the coolant container 10, and continuously measure the temperature ofthe borehole wall w at the respective points. These temperature sensors 40 are electrically connected with the connecting cables (not is shown), which are installed between the inner wall 11 and the outer wall 12 of the coolant container 10. The connecting cables are connected with the electric cables installed in the inner connecting holes 27 and 28 of the piston rod 23. Thus, the temperature of the borehole wall w, which is measured by the temperature sensors 40, is transmitted to a calculating unit outside the borehole h through both connecting and electric cables. 20 The crack detecting means functions to detect whether or not cracks are generated in the borehole wall w, at which position of the borehole wall w the cracks are generated, and in which diction the cracks are widened at the borehole wall w when the borehole wall w is cooled by the coolant r. In this embodiment, the method of crack detection employs a plurality of acoustic emission sensors 51, 52 and 53. The acoustic emission sensors 51, 52 and 53 are known as detectors that 25 detect acoustic emission generated between about 20 kHz and about 2 MHz. Each acoustic emission sensor detects a sound wave emitted when a minute internal crack is formed due to external shock or friction of a stmucture, converts the detected sound wave into an electrical signal, and transmits the converted electrical signal. FIG. 7 shows three acoustic emission sensors 51, 52 and 53 attached to the coolant container 12 10. As illustrated in FIG. 7, cryogenic heat caused by the coolant is applied to the borehole wall, and thereby a crack is formed at a point indicated by a symbol s to a point indicated by a symbol f. In this case, the three acoustic emission sensors 51, 52 and 53 detect a micro burst sound when the crack is formed. More specifically, when the crack is formed at a start point s and thus the burst sound occurs, 5 the acoustic emission sensor 53 located at a distance bl nearest the start point s detects the burst sound first, and then the acoustic emission sensor 52 located at a near distance eI next to the distance bl, and the acoustic emission sensor 51 located at a farthest distance dI sequentially detect the burst sound at an interval of time. Further, the acoustic emission sensors 51, 52 and 53 detect the burst sound at the final point f of the crack depending on distance b2, e2 and d2 from the final point f The 10 speed of the sound wave traveling at the borehole wall w is constant so that the position at which the sound wave occurs can be found using a difference of between the times detected by the acoustic emission sensors. Further, when a crack is widened, its widening direction can be found in the same method. Consequently, the three acoustic emission sensors 51, 52 and 53 detect micro burst sounds, so that they can measure whether or not a crack is formed, the position at which the crack is formed in 15 the borehole wall, and the direction(longitudinal direction, transverse direction, etc.) in which the crack is widened at the borehole wall w. Hereinafter, a method for measuring in-situ stress in rock using a thermal crack according to an exemplary embodiment of the present invention will be described in greater detail with reference to the accompanying drawings. 20 FIG. 2 is a flowchart illustrating a method for measuring in-situ stress in rock using a thermal crack according to an exemplary embodiment ofthe present invention. Referring to FIG. 2, a method Ml00 for measuring in-situ stress in rock using a thermal crack according to an exemplary embodiment of the present invention comprises a borehole forming step M10, an apparatus installing step M20, a first cooling step M30, a heating step M40, a second 25 cooling step M50, and a calculating step M60. The borehole forming step M10 is scheduled to construct a structure such as an oil storage tank, and thus forms a borehole h in the rock, which is aimed at measuring in-situ stress, at a preset depth. The borehole h can be formed in an inclined direction, but generally in a vertical direction. When the borehole h is formed, the apparatus installing step M20 is performed. To this end, 13 the in-situ stress measuring apparatus 100 of the present invention which is illustrated in FIGS. 4 through 6 is inserted into a measuring section i in the borehole h. When the in-situ stress measuring apparatus 100 is inserted, no pneumatic pressure is applied to the filling space 16. Thus, as illustrated in FIG. 4B, the outer circumference (i.e. the outer wall 12) of the coolant container 10 is kept without 5 being in close contact with the borehole wall w. Further, the coolant r is contained in the coolant chamber c of the cylinder 21. In this state, when the motor 30 is powered on, the rotor 31 of the motor 30 is rotated in one direction. The piston rod 23, which is screwed together with the rotor 31, linearly moves in an upward direction, thereby pushing the coolant r of the coolant chamber c in an upward direction. 10 When the electrically controlled inlet valve 18b of the coolant container 10 is opened, the coolant r evaporates while introduced into the containing space 15 ofthe coolant container 10, and thus is filled in the containing space 15. When a preset amount of coolant r is filled in the containing space 15, the inlet valve 18b is shut to stop filling the coolant r. At this time, liquefied nitrogen remains in the coolant chamber c for the purpose of secondary cooling, which will be described below. 15 In this state, when the pneumatic pump (not shown) is operated to inject air into the filling space 16 through the pipe 17b, pneumatic pressure is applied to the filling space 16, the inner wall 11 and the outer wall 12 of the coolant container 10 expands outwards, as illustrated in FIG. 4A. As a result, the outer wall 12 of the coolant container 10 comes into close contact with the borehole wall w. As described above, when the coolant container 10 is in close contact with the borehole wall 20 w, the fist cooling step M30 is initiated. In this closely contacted state, the coolant r (liquefied nitrogen) about -196 *C, primarily cools the borehole wall w. As illustrated in FIG. 3, the borehole wall w in the test section i is subjected to tensile stress so that a crack occurs at a point in the borehole wall w where tangential compressive stress is weakest. The crack grows in proportion to cooling time, and similar cracks occur at other places in the borehole wall w. In other words, after a crack is 25 formed, the crack gradually grows from within the rock, and thus the size of the crack, i.e., a gap caused by the crack, gradually increases. Further, when tensile stress is continuously added, cracks occur at other points other than the point of origin. In brief, when rock is cooled, tensile stress caused by heat is increased. And if this stress is greater than the compressive stress at any given point in the borehole wall, a crack will occur. 14 Each tempemture sensor 40 continuously measures the temperature of the borehole wall w, and transmits a temperature signal to an external calculating unit (not shown) through the electric cable installed in the connecting hole 27 in the piston rod 23. Thus, the temperature (first crack temperature) of the point when the crack occurs can be accurately obtained. The acoustic emission 5 sensors 51, 52 and 53 transmit signals detecting micro burst sound to an external calculating unit (not shown) through the electric cable installed in the connecting hole 27. The calculating unit yields the position (cracking point) where the crack is formed and the direction in which the crack is expanded on the basis ofthe transmitted signals. Thereby, the first cooling step M30 is completed. When the first cooling step M30 is completed, the heating step M40 of heating the cooled 10 borehole wall w is performed. The pressure decreasing valve 19b is opened to discharge the coolant r, which is contained in the coolant container 10, from the container. The pneumatic pressure of the filling space 16 is released to contract the coolant container 10 as illustrated in FIG. 4 B. Further, room-temperature air (or heated air) is introduced into the borehole h through the air injecting hole 26, which axially extends through the piston rod 23. The room-temperature air injected into the is borehole h heats the cooled borehole wall w while moving along the borehole h in an upward direction. When the borehole wall w is heated at room temperature due to the injection of the air, tensile stress applied to the borehole wall w is released, and thus the crack becomes smaller. In other wonds, the gap widened by the crack is reduced in size, and thus only an extremely minute gap remains. However, once a crack is formed in the rock, the crack cannot be restored to its original 20 state, and thus the crack is still present In the conventional hydraulic fracturing, hydraulic pressure is primarily applied to generate cracks, and then is released to make the crack smaller. This concept is equal to that of closing the crack by way of the heating step in the present invention. When checked from each temperature sensor 40 in the heating step M40 that the borehole wall w and the borehole h am sufficiently restored to room temperture, the second cooling step M50 25 is performed. More specifically, the coolant r remaining in the coolant chamber c is again filled in the containing space 15 of the coolant container 10 by the same process as the first cooling step M30, and pneumatic pressure is again applied to the filling space 16. Thereby, the outer wall 12 of the coolant container 10 comes into close contact with the borehole wall w. The borehole wall w is re-cooled and subjected to tensile stress. The crack, which is generated and widened in the first cooling step 15 M30 and then is closed in the heating step M40, is re-opened. This means that the cack closed in the heating step M40 is widened again However, as described above, since the crack itself is still present in the heating step M40, the crack is re-opened but at a relatively higher temperature than in the first cooling step M30. In other words, in order to generate the crack in the first cooling step M30, the 5 tensile stress caused by heat at the original cracking point must be more than both compressive stess and tensile strength of the rock combined. However, because the crack already occurred in the first cooling step M30, the tensile stress required to re-open the crack in the second cooling step M50 (the tensile stress generated by the coolant) has only to be more than the value of the compressive stress at that point. In other words, since the crack still remains, the force corresponding to the tensile strength 10 of the rock is excluded in the second cooling step M50, compared to the first cooling step M30. To be brief, the tensile stress required to re-open the crack in the second cooling step M50 is detennined by subtracting the tensile strength of the rock from the tensile stress required to re-open the crack in the first cooling step M30. When the second cooling step M50 proceeds, the crack generated in the first cooling step 15 M30 is re-opened. More specifically, at the start, the cracks are opened in a direction (longitudinal direction) parallel to the axial direction of the borehole I. These longitudinal cracks are located along the inner circumference of the borehole wall at several points. Further, when cooling is continued, the cracks are also opened in a direction (transverse direction) perpendicular to the axial direction of the borehole h. These transverse cracks are located along the axial direction of the borehole h at 20 several points. While the cracks are opened by cooling, each temperature sensor 40 continuously measures the temperature (second crack temperature) when each crack is re-opened, and the acoustic emission sensors 51, 52 and 53 continuously measure the positions in which the cracks are re-opened and the direction in which the cracks are widened, and then transmit the measurements to an extemal calculating unit. 25 In the calculating step M60, in-situ stress of the rock is calculated using the data, for instance, the first crack temperature, the second crack temperature, and the cracking point, which are obtained in the first cooling step M30 and the second cooling step M50. Here, in-situ stress in the test section is calculated using the signals transmitted from the temperature sensors 40 and the signals transmitted from the acoustic emission sensors 51, 52 and 53. 16 The in-situ stress calculated by the calculating unit includes two stresses, i.e., maximum horizontal principal stress and minimum horizontal principal stress, which act in directions perpendicular to each other on a plane (xy plane of FIG. 3) perpendicular to the axis of the borehole. Further, in the exemplary embodiment of the present invention, the in-situ stress also includes vertical 5 stress. A vertical stress acts in a direction (z-axial direction) perpendicular to each of the maximum horizontal principal stress and the minimum horizontal principal stess. First, the maximum horizontal principal stress, the minimum horizontal principal stress, and the vertical stress of the borehole wall will be described, and then the calculating step M60 of the in situ stress measuring method of the present invention will be described. 10 In the past, the stresses around the borehole h (stresses of the borehole wall w) were studied by many researchers, and can be obtained using the following Equations 1, 2 and 3, which are disclosed in a reference (Jaegar J. C., Cook N. G. W., 1976, Fundamentals of Rock Mechanics, 2'" Edn., Chapman and Hall, London). A tangential stress, a radial stress, and a shear stress around the borehole h in a homogeneous medium having an anisotropic in-situ stress on the plane (xy plane) 15 perpendicular to the z-axial direction z of the borehole h are as follows. (01+02) a2 (01-2) 3a 2 ' 2 (1+ 2)- (1+ 2 )cos26......(1) r 2 r (GI+G2) a 2 (11-G2) 4a 2 3a 2 a 2 (1-~- (2- 2 + 2 )cos20......(2) r 2 r r ( C-a2) 2 a2 3 a2 (1+ 2 - 2 )sin2o......(3) 2 r r where r is the radius of the borehole, a, and a2 are the maximum and minimum principal 20 stresses on an infinite plane perpendicular to the axial direction of the borehole (in which the maximum and minimum principal stresses am perpendicular to each other), and 0 is the angle measured in a counterclockwise direction starting at the maximum principal stress ar. When three dimensional stress is analyzed, the equations are resultant equations under a condition of plane strain, and are valid on the assumption that a circular opening is oriented in one of principal stress 25 directions. 17 At the borehole wall, the tangential stresses cao.0 and aso-go at 6 = 0* and 0 = 90* are as follows: ao-=3a 2 - a 1 ...... (4) ae-o= 3a, -a 2 ...... (5) 5 a,= 0, and -r0 = 0 at the borehole wall r = a,. The compressive force in a direction of the principal stress is taken as a positive value. In most of the cases, the tangential stresses are compressive in nature. Therefore, de-stressing of the borehole requires the induction of tensile stresses in the form of thermal stresses by lowering the temperature. As disclosed in a reference (Bmo A. B., Jerome H. W., 1985, Theory of Thermal Stresses, Dover Publications, Minneola, New 10 York), if it is shown through thermal stress analysis that the thermal stress does not cause the shear stress, the stress analysis can be simplified. Meanwhile, the tangential thermal stress aeA at the borehole wall w is given by the following Equation 6, as disclosed in the above-mentioned reference (Bruno A. B., Jerome H. W., 1985, Theory of Thermal Stresses, Dover Publications, Minneola, New York). E ar E,= a At=C(t- t1) 1s(1-v) .........(6) where E is the elastic coefficient (Young's module) of the rock, v is the Poisson's ratio of the rock, a is the linear thermal expansion coefficient of the rock, At is the temperature change fiom the temperature t before cooling to the cracking temperature t (first cracking tempemture or second cracking temperature) of the borehole wall, and C is the bi-axial thermo-elastic constant of the rock 20 and is expressed by C= Ea/(1-v). It can be predicted that the first point, at which the crack is formed from the borehole wall w on the plane (xy plane of FIG. 3) perpendicular to the borehole h, is a point (intersection of the x axis and the borehole wall) on which the maximum principal stress aj (assuming that the direction when M-0 is the direction of the principal stress) acts. In other words, since the tangential stress (acting in a 25 direction perpendicular to the direction of the maximum horizontal principal stress) at the point on which the maximum horizontal principal stress acts is shown on the borehole wall to the minimum extent, the tangential stress at that point is the fst to become equal to the tensile stress caused by the 18 coolant. Thus, when the following condition is met, the crack is formed from the borehole wall w on the plane perpendicular to the axis of the borehole. The following Equation 7 shows the concept that the value of adding the tensile strength at of the rock to Equaton 4 (value of adding the tangential stress and the tensile strength at the point where 0 =0*) is equal to Equation 6 (tensile strss caused 5 by heat). More strictly speaking, when the right side of the equality sign is even slightly greater than the left side, a crack is formed. 3; 2s 1 + (, =C Q ..... (7) where at is the tensile strength of the rock, and ti is the temperature at which the crack is forced. 10 Additional cooling of the borehole wall cannot only generate secondary cracks at various angles in a counterclockwise direction on the basis of the direction of the maximum principal stress, but also widen the existing cracks. When the following condition is met, a crack is formed in a direction perpendicular to the maximum principal stress (at a point on which the minimum principal stress acts, or at a point on the y axis). The following Equation 8 shows the concept that the value of 15 adding the tensile strength a, of the rock to Equation 5 (value of adding the tangential stress and the tensile strength at the point where 0 = 90*) is equal to Equation 6 (tensile stress caused by heat). More strictly speaking, when the right side of the equality sign is even slightly greater than the left side, a crack is formed. 3a- 2+ 0 I=C ta 2) ......... ....... (8) 20 where t 2 is the temperature of the borehole wall while the crack is formed in the direction of the minimum principal stress. However, it can be predicted that a crack is not easily generated at the point on which the minimum principal stress acts. This is because the tangential stress (acting in a diretion perpendicular to the minimum principal stress) at the point on which the minimum principal stress acts is shown on the borehole wall to the maximum extent. In the case of actual cooling, a 25 crack does not often occur at the point on which the minimum principal stress acts. Meanwhile, a crack (transverse crack) perpendicular to the axial direction of the borehole can be generated on the following condition. 19 ......... (9) where t3 is the terry e w heere the crack is forced in a direction perpendicular to the axial direction of the borehole, and a3 is the vertical stess and acts in the axial direction of the borehole. Solving Equations 7, 8 and 9, the principal stresses can be decided as follows. C ar a t= [ 3( ta- t 2)+( IC,-t 1)]-2 5 ......... (10) aJ2 C.[3( t,-t I )+( la- 12))- a g 2, (y3=C( 'a- ts3- (T 2 When the material constants E, v, a and at are given, the maximum horizontal principal stress ai, the minimum horizontal principal stress 02, and the vertical stress a can be decided. These 10 physical properties (the material constants) of rock are evaluated by testing core samples obtained at a test position. Meanwhile, when an equation is made using data obtained from the crack re-opening test (second cooling step), the tensile strength at can be removed from Equations 10, 11 and 12. This is because the crack that has already been generated is re-opened in the second cooling step and thus the ts tensile strength of the rock is removed. The thermo-elastic constant C can be obtained through a separate test instead of deciding the material constants E, y and a. The temperature (or sunounding temperature) t, before cooling can be measured, and the temperatures ti, t 2 and t 3 at which various cracks are generated are values measured using the in-situ 20 stmss measuring apparatus 100. When only the tensile strength at and the thermo-elastic constant C of the rock which are the unknowns are found, the three principal stresses 0, 02 and a3 can be decided. Since the tensile strength at and the thermo-elastic constant C of the rock are physical property values of the rock which can be measured at a laboratory, the three principal stresses can be theoretically calculated using Equations 10, 11 and 12. 20 In order to obtain in-situ stress using Equations 10, 11 and 12, a crack must be generated at the point on which the minimum horizontal principal stress acts. However, as described above, because a crack is rarely generated at the point on which the minimum horizontal principal stress acts, in-situ stress cannot be obtained using Equations 10, 11 and 12 if a crack is not generated at that point. 5 Further, although it is predicted in the first cooling step M30 and the second cooling step M50 that a crack will be formed first at the point on which the maximum horizontal principal stress acts and generated last at the point on which the minimum horizontal principal stress acts, the results may vary under real circumstances. For this reason, in the present invention, in-situ strss is obtained using a more generalized equation, which will be described below in detail. 10 Meanwhile, the vertical borehole will be theoretically described. From the worldwide stress measurement data, it can be said that a =yH where y is the unit weight, H is the depth (Brown E. T., Hoek E., 1978, Trends in Relationships between Measured In Situ Stresses and Depth, Int. . Rock Mech. Min. Sci. & 15 Geomech. Abst. 15,211-215; Amadei B, Stephansson 0, 1997, Rock Stress and Its Measurement l' Edn. Chapman and Hall, London). Thus, in the case of the vertical borehole, since it can be assumed that the axis of the boxhole is parallel to the vertical stress a 3 , ca can be considered as the already known value. Arranging the thermo-elastic constant C by substituting the above equation into as of Equation 9, the thermo-elastic constant C is as in the following Equation 13. C= =+ (T t H+aI 20 t (- t 3 t ~ f3 ........... Substituting Equation 13 into Equations 7 and 8, the following results are obtained. (yH+ a,) 3 - C I+ (= at3) ( t t 21 (yH+ a,) 3 t a- t2+ /,a(t- t2) (5 The components S(aI) and Sh(a2) of the horizontal principal stress can be obtained from Equations 14 and 15 as follows. (YH+at) s H= 3 t 12H t" t10" 8(ta-t13) 2 (16) (yH+a,) t, sh= _ 3( a to-t I ta- t21-- 8( ta- t3) 2 ......... (17) 'Ihus, in the first cooling step in the vertical borehole, only the value of the tensile strength, which can be obtained from the laboratory test, is required. Accordingly, the in-situ stress in the vertical borehole can be obtained when only the tensile strength of the rock is found. Moreover, as described above, the re-opening of the crack (second 10 cooling step) is performed on the vertical borehole, the constant at becomes 0, and thus the laboratory test can be omitted. In other words, in the second cooling step, Equations 16 and 17 can be simplified as follows. yH S H= [H3(ta- t2)+( t a- ti)] 8(ta- 13) (8 SH=............. ......(18) Sh= ta- t IH g(. ............. (19) 15 With respect to the vertical borehole, the value of in-situ stress can be decided in the above mentioned method. However, since Equations 18 and 19 are also made using Equations 10, 11 and 12, it is assumed that a crack is formed at the point on which the minimum horizontal principal stress acts, as 22 described above. Hence, it is necessary to use the generalized equation that can obtain in-situ stress although the crack is not generated at the point on which the minimum horizontal principal stress acts. Further, in order to obtain in-situ stress using Equations 10, 11 and 12 or Equations 18 and 19, the tensile strength and the thermo-elastic constant of the rock must be found. The tensile strength at 5 and the themo-elastic constant C are the physical property values of the rock, and thus can be obtained by testing the samples taken on the spot However, a slight deviation may exist between the values obtained at the laboratory and the real environment. As such, it is necessary to reduce any errors resulting from this deviation by obtaining as much data as possible. Hereinafter, the calculating process performed in the calculating step M60 of the present 10 invention will be described in detail. In order to minimize the range of error of in-situ stress measured and calculated according to the present invention, the following Equation 20, i.e., fust equation, is used as the general equation in which the maximum and minimum principal stresses a, and r 2 , the crack position 0, the temperature t,the them-elastic constant C, the tensile strength at etc., am included as variables. In the first is equation, the term of the right side ofthe equality sign is the thermal stress (tensile stress) expressed in Equation 6, and the term of the left side is the tangential stress (compressive stress) at an arbitrary angle 0. In the first equation, the variables am equal to those described in the above equations. The maximum horizontal principal stress ac and the minimum horizontal principal stress a2 can be 20 obtained with respect to the borehole wall on the plane perpendicular to the axis of the borehole using the following first equation. <First Equation> (1- 2 cos 2 O) a +(1+ 2 cos2O)a 2 + a,=C( 1,-t) In Equation 20, the maximum horizontal principal stress ar and the minimum horizontal 25 principal stress q2, which are to be measured, are the unknowns. The tensile strength a, and the thenno-elastic constant C of the rock are the physical property values of the rock, and thus can be obtained by testing at a laboratory, but may be treated as unknowns. The crack position 0, the temperatures t and t are the values measured by the in-situ stress measuring apparatus 100. In the 23 case in which the tensile strength and the thermo-elastic constant of the rock are treated as the unknowns, a total of four unknowns are included in the first equation. Thus, a crack is formed at four or more points in the borehole wall by cooling, and then four or more first equations are made and solved simultaneously. Thereby, the maximum horizontal principal stress and the minimum 5 horizontal principal stress can be obtained. However, in the case in which the tensile strength at and the thenmo-elastic constant C of the rock are separately measured at a laboratory by taking rock samples, only two unknowns are included in the first equation. Thus, a crack is frned at two or more points in the borehole wall, and then two or more first equations ae made and solved simultaneously. Thereby, the maximum horizontal 10 principal stress and the minimum horizontal principal stress can be obtained. Similarly, when only one of the tensile strength atand the themo-elastic constant C of the rock is obtained at a laboratory, three unknowns remain in the fust equation. Thus, a crack is forced at three or more points in the borehole wall, and th ree or more first equations are made and solved simultaneously. Tbereby, the maximum horizontal principal stress and the minimum 15 horizontal principal stress can be calculated. However, in order to reduce the enors of the values of the maximum and minimum horizontal principal stresses calculated by the first equations, as much data must be secured by applying the first equations to the cracks generated at several points in the borehole wall. When a least square method is performed on a plurality of solutions of the maximum and minimum 20 horizontal principal stresses obtained by the plurality of first equations secured from the several points where the cracks are generated, least square solutions are calculated with respect to the maximum and minimum horizontal principal stresses. Thereby, the error range can be minimized. The first equation can be made using the cracks generated in the first cooling step M30 or the cracks re-opened in the second cooling step M50, When the equation is established in the first 25 cooling step, 6 will be the point at which the crack occurs, and t will be the temperature at which the crack occurs. When the equation is established in the second cooling step, 0 will be the point at which the crack is re-opened, and t will be the temperature at which the crack is re-opened. The re opened crack in the second cooling step is one that has already been generated in the first cooling step, so 0 is the same as in the first cooling step. Further, in the case in which the first equation is 24 established using the re-opened crack in the second cooling step, the tensile strength a, of the rock is set to 0 in the first equation, because the crack has already been generated in the first cooling step M30, and thus tensile strength of the rock is removed when the crack is re-opened in the second cooling step M50. As described above, the first equation can be established in the first cooling step 5 M30 and the second cooling step M50. The plurality of first equations is preferably formed by using the data, i.e., the point and temperature at which the crack is re-opened. This is because, in the first cooling step M30, a release of stress can occur around the crack, which is generated in the initial stage of the first cooling step M30, and which has an influence on stress concentration required to subsequently generate cracks. If so, the accurate value cannot be calculated, and thus errors can 10 occur. Meanwhile, the vertical stress as acting in a vertical direction is independent not of ac and 02 but of the thermo-elastic constant C and the tensile strength at 'he vertical stress 03 must be obtained by deciding the thermo-elastic constant C and the tensile strength amusing the equations based on the first equation (Equation 20), on a condition that the transverse crack is generated. 15 <Second Equation>
;
3 =C( t,- t 3 )- aI In general, transverse cracks can be initiated on the plane perpendicular to the axis of the borehole in directions of 0* and 1800 with respect to the direction of the maximum principal stress. These cracks will gradually be propagated while forming complete circular cracks along the 20 circumference of the borehole wall. In this case, unlike the case where the longitudinal cracks are generated, a great difference may exist between the temperature t 3 i at which each crack is initiated and the temperature tk at which each crack is completed. Thus, the temperature t 3 at which the transverse crack is formed has only to be calculated by a simple average of the temperature t 3 ; at which each crack is initiated and the temperature t3e at which each crack is completed, as in the following 25 equation. t 3 =( t3i+ t3e)/2 Meanwhile, when the equations based on Equation 13 are obtained, preferably they are obtained only when the initial crack is formed in the fust cooling step, i.e., only the point at which the 25 crack position 6 is 0*, and the equations for the subsequent cracks are obtained in the second cooling step M50. The reason why the equations for the subsequent cracks are not obtained in the first cooling step is because stress release occurs around the initial crack in the first cooling step and because the initial crack has a serious influence on the stress concentration required to generate the 5 subsequent cracks, and thus because the accurate values cannot be calculated. As described above, the borehole wall w is subjected to primary and secondary cooling using the apparatus according to the present invention, and thereby many cracks are generated. The equation is established for each crack using Equation 13, and thereby many solutions (maximum and minimum horizontal principal stresses, and thenno-elastic constant) can be obtained. The least 10 square method is applied to these solutions, and thus the least square solution is calculated. Thereby, the solution minimizing the enor range can be obtained. The vertical stress also can be obtained by minimizing the error range in the above-mentioned method. Numerical modeling studies are conducted using so-called 'FLAC3D' program with thennal option (FLAC3D manual, 1997) to check the validity of the above-mentioned method. The crack Is development is modeled by making the elements as null zones when the tensile stresses exceed the tensile strength of the rock. Though several other methods were also suggested by various researchers in the past for the modeling of generation and propagation of cracks in brittle materials, this method is simple and easy to use. This technique has also been used in the past with finite element method for predicting the fragment formation while blasting (Saharan M. R., 2004, Dynamic Modeling of Rock 20 Fracturing by Destress Blasting, Ph. D. Thesis, University of McGill, Canada). When the direction of generation and propagation of cracks is not known prior to the modeling, this is a very effective method for crack simulation, provided the finite element/difference mesh is fine enough. Large element size near the borehole wall can result in significant mechanical stress unbalance while cracking, and therefore should be avoided. 25 In order to simulate the size of a so-called NX borehole, a borehole having a diameter of 74 mm is made in a plane stain state using a single strip model. Only one quadrant is created using a condition of rotation symmetry. In order to simulate conditions of the vertical borehole, the direction of the borehole is set to the z-axial direction. The maximum horizontal principal stress SH and the minimum horizontal principal stress are considered parallel to the x axis and the y axis, 26 respectively. At the inside of the borehole, the coolant container (made of low stiffness material) is modeled, which would hold the cryogenic fluid (LN 2 ) at a temperture of -196 0 C. Transient heat flow analysis is performed while heat conducts through the coolant container to the borehole wall and into s the rock. Coupled thermal and mechanical stepping is performed until the tensile stresses in the rock (borehole wall) exceed the tensile strength of the rock. A so-called 'FISH' program checks the stress level of each element in the mesh during every calculation cycle. For any element, when the tensile stress exceeds the tensile strength of the rock, the element is changed to a "null element". This will result in significant unbalance of forces. Therefore, the 10 thermal cycling is stopped and the model is brought to mechanical equilibrium by performing only the mechanical cycles. Coupled thermal-mechanical cycles are resumed thereafter. In other words, when one process is completed in this way, the thermodynamics interction analysis is re-started. This process is continued until fractures (cracks) develop through out the periphery of the borehole wall The temperature of the borehole wall at the time of the first appearance of crack (usually along 15 the Su direction) is taken as ti and the appearance of subsequent fractures at angles, E 2, 0 3, E 4 with respect to the SH direction is taken as t 2 , t, t 4 , respectively. The examples (Example 1) performed here have the following general input data: Thermal and mechanical properties for the surrounding rock are as follows: Linear thermal expansion coefficient a=6.64e-6 PC ; Thermal conductivity k=2.63 W/mPC; 20 Specific heat at constant pressure C,=710 J/kg/C; Density p -2500 kg/m 3 ; a 1 =5 MPa; E=20GPa; v=0.25 Thermal and mechanical properties for the coolant container are as follows; Linear thermal expansion coefficient a=2.0e-7 / *C; Thermal conductivity k=3.0 W/m/ ; C,,=500 J/kg/"C ; Density =1000 kg/m 3 ; E=0.02GPa; v=0.25 25 A plane strain condition is created using a single strip model of a borehole of diameter 74mm (corresponding to NX size). Only a quarter of the borehole and the surrounding rock is created by using planes of symmetry passing through the borehole axis and along x and y axes. The borehole is oriented in the z axis direction to simulate the condition of a vertical borehole. The maximum horizontal principal stss (SH) is considered parallel to the x axis and the minimum horizontal 27 principal stress (Sb) parallel to they axis. Initial stress level taken for this study is as follows: ac=10 Maag=7 MPa, an=20 MPa, a,=O MPa, ay/=O MPa, a=O MPa Fig. 8 shows various stages of the thenno-mechanical modeling of the borehole cooling and 5 fracturing (cracking). The sequence of fracturing (cracking) during the initial cooling cycle can be visualized from Fig. 8A-8J). As exeted, the cracks el initiate at the borehole peiphery and in the direction of x axis, which is the maximum horizontal principal stress (SH) direction. Subsequently, second and third cracks c2 initiate at 9.50 and 17.5* respectively with the SH direction. This can be seen fiom Fig. 8B and 8C, respectively. All this above fractures can be treated as radial fractures 10 parallel to the borehole axis. In general, it is seen from the numerical modeling that the cracks initiate at intervals of 7-10. As new cracks are generated, it is also observed that the earlier cracks deepen to 2-3 elements' depth (Fig.8D-81). Seven cracks are observed at angles 0, 9.5, 17.5. 25.5, 32.5, 41.5 and 53.50. During the numerical modeling analysis, the wall temperature was also monitored in tandem with the fractum monitoring. The above seven crack initiations are found at wall temperatures is -56.85, -66.12, -73.6, -85.16, -89.76, -102.9 and -113.6*C, respectively. In this particular study, the out-off-plane stress o. is 2OMPa, which is less than the tangential stress at 0 =90 * (23MPa), therefore the crack perpendicular to the SN direction could not be studied, as cracks (transverse cracks) perpendicular to the borehole axis initiated prior to the fracture along y axis. It requires a more elaborate three-dimensional modelling to study fracturing in all the 20 three directions. In order to analyze the results of numerical modeling for potentials and verification of a new method for measuring in-situ stress, the theoretical stress values at the borehole periphery is compared with the numerical estimates of thermal stresses during the cracking process at various angles. 25 Fig. 9 shows the theoretical tangential stresses existing through out the borehole periphery at various angles with the x axis. These theoretical values are calculated using equation (1) for a=r. The numerical estimates for the seven crack angles are calculated from the wall temperature at the time of crack initiation using equation (6) and plotted in Fig. 9 after subtracting the tensile strength factor. From Fig. 9 the numerically estimated tensile thermal stress value (-oih) required for fracturing at 28 e=o, i.e. for the first crack cl, matches well with the theoretical tangential stress value. However, for the subsequent cracks c2, the numerical estimates are found over predicting the existing stress levels. This discrepancy is caused due to the stress release near the earlier cracks, which significantly affects the stress concentration needed for the subsequent crack initiations. Therefore, with the first 5 cycle of cooling(first cooling step), only one equation, viz. equation (7) can be obtained accurately. It is further speculated that when the tensile strength of the rock is lower compared to the ambient stress levels, this discrepancy may reduce. The depth of these fractures are in the order of a few millimeters, which is much lower compared to the hole diameter, and therefore in real testing these cracks may not significantly alter the stress picture of the surrounding rock 10 To simulate the second cycle of cooling (second cooling step) and crack reopening, it is necessary to thaw the borehole and bring it back to its original temperature state. This can be practically done by withdrawing the cryogenic fluid from the coolant container and blowing the borehole periphery with ambient temperatre air for a few minutes and at the same time monitoring the temperature, until it is stabilized. This will allow the earlier cracks to close and the compressive 15 stresses to act across them. Further re-cooling of the borehole generates tensile stresses, which will cause the earlier cracks to reopen at earlier angles at specific wall temperatures. Since the cracks are tension free surfaces, the tensile strength factor automatically disappear from the stress equation (equation (20)). This enables to obtain the solution of the remaining two unknown principal stresses using any two equations or to obtain a best-fit using multi-linear regression employing least square 20 analysis. The crack reopening, second-cycle cooling test is simulated using the following procedure. a) Mapping of the cracks caused by the fist-cycle cooling test. b) Changing the crack elements back to "elastic" and flagging them as "no-tension"material. c) Stabilizing the temperature back to the ambient and thus nullifying all the thermal stresses 25 in the medium. d) Re-application of -I 96*C at the coolant container periphery and start of the second-cycle cooling. e) Monitoring the wall temperature as well as crack reopening using a FISH program. 29 Since the crack elements are declared "no-tension" materials, the cracking is expected to initiate for those elements as the pm-existing compressive stresses am nullified by the tensile thermal stresses. The sequence of crack reopening during the second cycle of cooling is given in Fig. 10. The first crack ci reopening can be seen in Fig. lOB along the x axis. The wall temperature of reopening is 5 -30.5*C, which is lower compared to the first-cycle cooling test. Subsequent failure and crack extension can be seen at the angles, 9.5,17.5,25.5,32.5,41.5 and 53.5*, as shown in Fig. IOC- 10H. The wall temperatures during the crack reopening along the above angles are observed from the history plots as -32.97,-37.6, -41.1, -50.1, -62.17 and -75.3*C, respectively. The numerical estimates of thermal stresses at crack reopening are plotted in Fig. 11 along to with the pre-existing theoretical tangential stresses. One point to be noted is that for crack reopening, the tensile strength is taken as zero. From Fig. 11 it can be seen that the numerical estimates are very close to the theoretical values. Therefore, the crack reopening temperatures at any two angles can be used to resolve the two principal stresses, SH and S.,. However, the solution gives the principal stresses in terms of the thermo 15 elastic constant C. As stated easier, for a vertical borehole, C can be estimated as C= yH (t, - t3) where t-, is the wall temperature of reopening of a horizontal crack perpendicular to the borehole axis. But for all other conditions, it has to be estimated in the laboratory. In a two-dimensional analysis, the crack generation perpendicular to the borehole axis cannot be studied; this requires a three dimensional analysis. 20 If we know the value of C, which is equal to 0.1771 Wa/"C in this particular case, the principal stresses can be estimated. Best fit solutions are obtained by least square analysis for the principal stresses by considering considering 7 crack angles, which gives q1=10.17 MPa and or=7.06 MPa with and
R
2 =0.99. These values are very close to the applied stress values of a. and ai, in the numerical 25 model. The first example, being a two dimensional analysis, was unable to model the thermal cracks in planes perpendicular to the borehole axis. Therefore, a three-dimensional model is constructed in FLAC3D by taking a considerable length of the borehole section, where the low temperature is to be 30 applied. The procedure adopted for the 3D modelling is otherwise similar to that described in Example 1. The initial stress levels taken for this example are:a.=5 MPa, ag2 MPa, an=3 MPa, a=O MPa, ay=O MPa, a.=O MPa. 5 No shear stresses are applied in the Cartesian planes, which imply that the principal stresses are aligned along the axes, and the borehole is aligned along one of the principal stress direction viz., the z axis. During the initial cooling test run, the tensile cracks initiated through out the modeled section of the borehole, in a plane parallel to the borehole axis and at zero degrees to the x axis, which is the 10 maximum horizontal principal stress direction. This initial crack formation is shown in Fig. 12A. Similar cracks in planes parallel to the borehole axis and at various angles to the x-axis can be seen in the subsequent figures Fig. 12B-12H. At the same time, cracks perpendicular to the borehole are also found to originate from the direction of x axis (maximum principal stress direction) and propagate angularly along the borehole 15 wall until it reaches the direction of y axis, depicting a full-circle perpendicular crack. This crack formation can be seen from Fig. 12B-12R. Such perpendicular cracks are found to occur at intervals of 5-8 mm along the length of the borehole. To attribute an average wall temperature for the crack formation, a simple average of the temperature of crack initiation and crack completion is taken as, 13 =(t3 +t3,)/2......(21) 20 If the principal stresses in the other directions viz., o and o ae not widely different, it is expected that these cracks occur within a brief span of time and wall temperature. In this particular case, the average temperature of perpendicular crack formation is -15.3"C. The corresponding stress value for the perpendicular crack formation, computed using equation (12) is 2.14 MPa. ui of Equation (12) and (20) can be eliminated from multiple crack re-opening tests. Ths 25 makes it possible it decide the principal stresses using only the thermo-elastic constant C. From Fig. 12, six crack angles can be observed starting from 0* to 59.25* with the x axis. The wall temperatures corresponding to crack initiations at 0, 11.25, 23.25, 35.25, 47.25 and 59.250 are 10.7, -16.07, -27.8, -45.96, -69.63 and -85.03*C, respectively. Fig. 13 shows the comparison between 31 the theoretical tangential stress values and the predicted numerical estimates. For the initial crack formation, as in the case of example 1, the numerical prediction agrees well with theoretical values, but for the subsequent angles, the numerical values seem to over predict the tangential stress levels. A multi-linear regression analysis of the numerically predicted data 5 gives cy=6.8MPa, c 2 =2.3MPa and at=6.47MPa. These values are within 15-35 % error. A second-cycle cooling(second cooling step) for crack reopening modeling is conducted using similar procedure followed for the previous example. The cracks recorded during the first cycle are mapped, with some approximations. Fig.14A-14F shows the different stages of crack reopening. Referring to Fig.14, the temperatures corresponding to the radial 10 cracks reopening at angles 0, 11.25, 23.25, 35.25, 47.25 and 59.25' are 22.57, 17.1, 8.75, -5.5, -23.0, -39.5'C, respectively. The perpendicular cracks reopening initiated at an approximate wall temperature of 17'C (Fig.14B), which completed a full circle at about 2.4 0 C, averaging about 9.7'C. The ca level corresponding to this average temperature is 2.7MPa, which is within 15% error. Fig. 15 shows a comparison of the theoretical pre-existing tangential stress 15 levels through out the borehole periphery and the numerically predicted values with respect to the crack reopening thermal stresses. As in the case of example 1, the numerical predictions are in good agreement with the theoretical values for the reopening test. A multi-linear regression using least square analysis gives the best-fit values for aYx and yy as 5.78 and 2.13 MPa, respectively. These predictions estimated the principal stresses within 15% error. 20 The preliminary theoretical and numerical analysis conducted to validate the novel concept shows that the method is capable of accurate estimation of the in-situ stresses by monitoring tensile crack formation at borehole wall and the wall temperature at the time of crack initiation. The time taken for cooling the borehole periphery for the test will only be a few seconds. The method can be applied to boreholes of any directions. Moreover, all three 25 principal stresses can be obtained by monitoring crack formation at various angles along the borehole axis as well as across the borehole axis. The errors obtained in determination may be those associated with the numerical modeling procedure. With better and more sophisticated crack modeling procedures, it is expected that the error levels may come down. Where the terms "comprise", "comprises", "comprised" or "comprising" are used in 30 this specification (including the claims) they are to be interpreted as specifying the presence of the stated features, integers, steps or components, but not precluding the presence of one or more other features, integers, steps or components, or group thereto. 32

Claims (11)

1. An apparatus for measuring in-situ stress in rock using a thermal crack by drilling a borehole into the ground, and by forming a crack caused by heat in a wall of the borehole, the apparatus comprises: 5 a coolant container, in which an annular containing space where a coolant can be contained, and an inlet for the coolant to flow into and out of the containing space are formed; at least one temperature sensor, which is installed on an outer surface of the coolant container so as to measure temperature of the borehole wall; and a crack detecting means for detecting cracks formed in the borehole wall by heat 10 transmission between the coolant contained in the coolant container and the borehole wall; wherein the crack detecting means detects a burst sound occurring when a crack is formed from the rock with the borehole wall, and comprises a plurality of acoustic emission sensors attached to an outer circumference of the coolant container.
2. The apparatus according to claim 1, further comprising a close-contact means for 15 bringing an outer surface of the coolant container into close contact with the borehole wall.
3. The apparatus according to claim 2, wherein: the coolant container comprises an inner wall, an outer wall disposed outside the inner wall while spaced apart from an outer circumference of the inner wall, an upper cover coupled to the upper sides of the inner and outer walls, and a lower cover coupled to the lower sides of 20 the inner and outer walls; and the annular containing space is surrounded and defined by the inner and outer walls and the upper and lower covers.
4. The apparatus according to claim 3, wherein: the coolant container is provided therein with a filling space, which is surrounded by 25 the inner wall, the upper cover, and the lower cover of the coolant container so as to fill and discharge a fluid; the close-contact means uses a pump connected to the filling space so as to fill and discharge the fluid into and from the filling space; and the outer and inner walls of the coolant container are made of an elastic material, 30 which expand and contract as the fluid is filled into and discharged from the filling space.
5. The apparatus according to claim 1, wherein the coolant contained in the coolant container is liquefied nitrogen. 33
6. The apparatus according to claim 1, further comprising a coolant injecting means for injecting the coolant into the containing space.
7. The apparatus according to claim 6, wherein the coolant injecting means comprises: a cylinder, which is coupled to the coolant container such that a coolant chamber, in which the 5 coolant is contained, is formed between the coolant container and the cylinder; a piston, which has a piston head inserted into the cylinder, and a piston rod having a bar shape and fixed to the piston head, which is installed so as to be able to be reciprocated in the cylinder, and which pushes the coolant towards the coolant container; a driver, which reciprocates the piston; a valve, which is installed on the coolant container so as to mutually communicate and shut the 10 coolant chamber and the containing space.
8. The apparatus according to claim 7, wherein: the driver comprises a reversible motor, which includes a rotor having a hollow shape and a female thread on an inner circumference thereof; and the piston rod has a male thread on an outer circumference thereof, is screwed into the 15 rotor, and is linearly reciprocated when the rotor is rotated in forward and reverse directions.
9. The apparatus according to claim 8, wherein: the piston rod extends through the piston head, the coolant chamber of the cylinder and the coolant container, and protrudes from the coolant container at one end thereof and from the motor at the other end thereof; and 20 an air injection hole passes through the opposite ends of the piston rod in a longitudinal direction of the piston rod.
10. The apparatus according to claim 1, further comprising a pressure decreasing valve, which is installed on the coolant container in order to open and shut the containing space of the coolant container and to discharge the coolant of the containing space. 25
11. An apparatus for measuring in-situ stress in rock using a thermal crack substantially as hereinbefore described with reference to any one of the embodiments shown in the drawings. 34
AU2011250753A 2006-10-31 2011-11-11 Method and apparatus for measuring in-situ stress of rock using thermal crack Ceased AU2011250753B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
AU2011250753A AU2011250753B2 (en) 2006-10-31 2011-11-11 Method and apparatus for measuring in-situ stress of rock using thermal crack

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
KR10-2006-0106208 2006-10-31
KR10-2007-0101418 2007-10-09
KR10-2007-0101406 2007-10-09
AU2007314874A AU2007314874B2 (en) 2006-10-31 2007-10-30 Method and apparatus for measuring in-situ stress of rock using thermal crack
AU2011250753A AU2011250753B2 (en) 2006-10-31 2011-11-11 Method and apparatus for measuring in-situ stress of rock using thermal crack

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
AU2007314874A Division AU2007314874B2 (en) 2006-10-31 2007-10-30 Method and apparatus for measuring in-situ stress of rock using thermal crack

Publications (2)

Publication Number Publication Date
AU2011250753A1 AU2011250753A1 (en) 2011-12-08
AU2011250753B2 true AU2011250753B2 (en) 2012-12-20

Family

ID=45465628

Family Applications (1)

Application Number Title Priority Date Filing Date
AU2011250753A Ceased AU2011250753B2 (en) 2006-10-31 2011-11-11 Method and apparatus for measuring in-situ stress of rock using thermal crack

Country Status (1)

Country Link
AU (1) AU2011250753B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111678994A (en) * 2020-07-01 2020-09-18 山西宏安翔科技股份有限公司 Acoustic emission sensing system and method
CN113155589B (en) * 2021-05-17 2025-08-05 中南大学 A method for studying the drilling process and the development of cracks in the surrounding rock of deep high-temperature and high-pressure formations

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100089585A1 (en) * 2006-10-13 2010-04-15 Kaminsky Robert D Method of Developing Subsurface Freeze Zone

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100089585A1 (en) * 2006-10-13 2010-04-15 Kaminsky Robert D Method of Developing Subsurface Freeze Zone

Also Published As

Publication number Publication date
AU2011250753A1 (en) 2011-12-08

Similar Documents

Publication Publication Date Title
AU2007314874B2 (en) Method and apparatus for measuring in-situ stress of rock using thermal crack
Pakhotina et al. Diagnosing multistage fracture treatments with distributed fiber-optic sensors
AU2020217344A1 (en) Methods for estimating hydraulic fracture surface area
US11125061B2 (en) Simultaneous injection and fracturing interference testing
WO2022056194A1 (en) Method and system for determining energy-based brittleness
MX2013014709A (en) Systems and methods for measuring parameters of a formation.
US8408296B2 (en) Methods for borehole measurements of fracturing pressures
CN113626972B (en) Method and device for evaluating dynamic sealing performance of gas storage trap using stress field
Wiktorski et al. Experimental study of temperature effects on wellbore material properties to enhance temperature profile modeling for production wells
Bohloli et al. Determination of the fracture pressure from CO2 injection time-series datasets
Xiang et al. Numerical investigation of borehole breakout and rock spalling based on strain energy criteria
US10550687B2 (en) Methods for analyzing formation tester pretest data
AU2011250753B2 (en) Method and apparatus for measuring in-situ stress of rock using thermal crack
Seagraves et al. Fundamental wormhole characteristics in acid stimulation of perforated carbonates
Read et al. Technical summary of AECLs Mine-by Experiment phase I: Excavation response.
CN120764183B (en) A method for predicting easily leaky formations of drilling fluid using well logging data
TerHeege et al. Discrete element modelling of wellbore integrity in high temperature geothermal reservoirs
Zillur et al. The role of geomechanical earth modeling in the unconsolidated Pre-Khuff field completion design for Saudi Arabian gas wells
CN115541393B (en) Experimental system and method for determining dynamic closure characteristics of cracks
Khaksar et al. Use of combined analytical and numerical sand production prediction methods to define simple cut-off for sand control and well completion selection
Guglielmi et al. A new in-situ approach for hydromechanical characterization of mesoscale fractures: the High-Pulse Poroelasticity Protocol (HPPP)
Pei et al. Constraining in-situ stresses at BETA by analysis of borehole images and downhole pressure data
Xiaodong et al. Calculation of fluid-producing profiles from mine fiber-optic acoustic data: Model and application
Ramos et al. Development and Testing of Advanced Inter-Well Pressure Pulse Analysis for Fracture Diagnostics in Tight Gas Reservoirs
Islam et al. Calibration of rock strength and yielding by using hollow cylinder experimental data

Legal Events

Date Code Title Description
FGA Letters patent sealed or granted (standard patent)
MK14 Patent ceased section 143(a) (annual fees not paid) or expired