Data processing system for mapping fracture length using downhole ground penetrating radar

Abstract

Ground penetrating radar (GPR) measurements from a downhole well tool in a wellbore are obtained to identify length of fractures adjacent the wellbore. A ground penetrating radar transmitter of the downhole tool emits an electromagnetic pulse. The electromagnetic wave of the ground penetrating radar is diffracted on encountering an end or tip of a fracture, which acts as a secondary source. The diffracted signal is then collected by downhole receiver(s) of the downhole tool. Length of the fracture is determined based on the time of travel of the electromagnetic wave from its emission until its collection as a diffracted signal by the downhole receiver(s).

Claims

1. A data processing system for determining location of end portions of a fracture in a subsurface formation penetrated by a well, comprising: (a) an uphole telemetry system receiving, from a well logging tool in the well, an uphole signal comprising measures of sensed time of receipt of diffracted ground penetrating radar electromagnetic energy sensed by a radar receiver in the well logging tool after emission of electromagnetic energy by a ground penetrating radar electromagnetic energy transmitter in the well logging tool; (b) the uphole telemetry system further receiving, from the uphole signal, sensed arrival of a first break pulse indicating spacing of the radar transmitter and radar receiver, and sensed arrival of a second break pulse indicating presence of a fracture as a secondary electromagnetic energy source; and (c) a processor computer processing the received uphole signal to determine, from the measures of the sensed time of receipt of the sensed electromagnetic energy, the location of end portions of the fracture in the subsurface formation.

2. The data processing system of claim 1, wherein the processor computer further determines the location of end portions of the fracture by determining the location of end portions of the fracture from an arrival time of the second break pulse.

3. The data processing system of claim 1 wherein the fracture in the formation is formed by fracturing the subsurface formation.

4. The data processing system of claim 1 wherein the well is formed by horizontal drilling of the subsurface formation.

5. The data processing system of claim 1 further including: a memory storing the determined location of end portions of the fracture in the subsurface formation.

6. The data processing system of claim 1 further including: an output display forming an image of the determined location of end portions of the fracture in the subsurface formation.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a schematic diagram, partly in cross-section, of a ground penetrating radar system according to the present invention deployed in a subsurface formation penetrated by a horizontal wellbore.

(2) FIG. 2 is a schematic diagram of fracture length determination according to the present invention with the ground penetrating radar system of FIG. 1.

(3) FIG. 3 is schematic diagram of a laboratory model used in test mapping of fracture length with a ground penetrating radar system according to the present invention.

(4) FIGS. 4, 5, 6, 7, 8 and 9 are plots of received radar signals as functions of time obtained according to the present invention during testing of the laboratory model of FIG. 3.

(5) FIGS. 10A, 10B, 11A, 11B, 12A and 12B are comparative plots of mapped fractures of different lengths obtained according to the present invention during testing of the laboratory model of FIG. 3 for different ground penetrating radar frequencies.

(6) FIGS. 13A, 13B and 13C are comparative plots of mapped fractures obtained according to the present invention during testing of the laboratory model of FIG. 3 with different reservoir conductivities.

(7) FIG. 14 is a schematic diagram of the data processing system of the ground penetrating radar system of FIG. 1.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

(8) According to the present invention, a method is provided to mapping fracture length with downhole ground penetrating radar waves that eliminates the need for nearby monitoring wells and surface sensors. The present invention takes advantage of what is known in physics as the knife-edge effect. The knife-edge effect is based on the Huygens-Fresnel principle, which is that a well-defined obstruction to an electromagnetic wave acts as a secondary source of an electromagnetic wave. The obstruction acting as a secondary source creates a new electromagnetic wavefront appearing to emanate from the obstruction.

(9) In the drawings, a ground penetrating radar well logging L is shown in FIG. 1. The ground penetrating radar well logging L includes a downhole ground penetrating radar Tx which emits when activated a pulse of electromagnetic wave energy that travels into the earth formation. A portion of the emitted pulse of electromagnetic wave energy travels directly from the transmitter Tx to the receiver Rx and is detected as a first break pulse having a first travel time. When there is a fracture 20 present in a subsurface formation in which the well logging tool is present, a portion of the emitted pulse of electromagnetic wave energy is diffracted by a tip 22 of fracture 20. The diffracted signal is then sensed by one or more downhole receiver(s) Rx. Length of the fracture 20 or distance to the tip 22 is then determined, as will be explained, from the travel time of the diffracted ground penetrating radar electromagnetic pulse sensed at the receiver(s) Rx.

(10) The well logging tool L includes a sonde or housing body 30 which is suspended for movement in a wellbore 32 for movement by a wireline logging cable 34. The well logging tool T is moved in the wellbore 32 to well depths of interest in a formation 36 which is of interest for characterization of the length of fracture 20. As shown in FIG. 1, the wellbore 32 for an oil or gas shale formation 36 of an unconventional reservoir has transitioned from vertical near the earth surface to a generally horizontal plane by directional drilling.

(11) Measures of travel time of the diffracted ground penetrating radar electromagnetic pulse sensed at the receiver(s) Rx of the downhole logging tool L are telemetered by a downhole telemetry system 38 in the downhole sonde 30 over the wireline logging cable 34 to an uphole telemetry and data processing system D (FIGS. 1 and 14). Travel time measurements from the logging tool are received by the uphole telemetry and data processing system D in a data input unit 42 as functions of borehole depth or length of extent in wellbore 32.

(12) A preprocessing unit 44 and a processor computer 46 receive and process the measures of travel time of the diffracted ground penetrating radar electromagnetic pulses so that length of the fracture 20 present adjacent the wellbore 32 at a depth of interest may be determined. The determined fracture lengths at depths of interest are then stored or recorded in memory 48 and displayed by action of a display 50 in the desired manner, which is usually a plot of the parameters of interest as a function of depth at which they are determined within the wellbore 32.

(13) A surface depth measurement system, such as a depth measure sheave wheel 52 and associated circuitry is provided to indicate depth of the logging tool L in the wellbore 32. The travel time data from the downhole subassembly 10 are recorded as functions of borehole depth in the data memory 48. Once recorded, the travel time measurements are transferred as needed into the preprocessing unit 44 and processor computer 46 of the data processing system D. The travel time measurements are after transfer subjected to conventional preprocessing in the preprocessing unit 44 and then transferred to computer 46. The processed data measurements in computer 46 are then available for processing to determine fracture length of fractures such as shown at 20 according to the present invention in a manner to be set forth below.

(14) A series of numerical simulations were performed to demonstrate capabilities of the present invention. A numerical two-dimensional or 2-D model (FIGS. 2 and 3) was created representing a section of an unconventional reservoir), with a section of horizontal wellbore 60 and a fracture 62 that is perpendicular (FIG. 2) to the wellbore 60. The length of the fracture 62 to a tip 64 was configured in the numerical model to be variable to confirm the sensitivity of the present invention to the fracture length. A location 66 for transmitter Tx and a location 68 for receiver Rx were configured in the numerical model to be located inside the wellbore 60, spaced equidistant on either side of the fracture 62. The positions 66 and 68 of transmitter Tx and receiver Rx were moved for several different simulation scenarios to determine the effect of transmitter and receiver position on the received ground penetrating radar electromagnetic pulse signal.

(15) In an initial test scenario (FIG. 3) for the numerical model, a 1 cm wide and 5 m long fracture 62 was simulated. The fracture was mapped according to the present invention using a 2 GHz electromagnetic pulse. The transmitter Tx location 66 and the receiver location 68 in the numerical model were 10 m apart from each other. Values for subsurface formation electromagnetic field parameters .sub.r, , and .sub.r in the numerical model are indicated in FIG. 3. An electromagnetic pulse was then generated at by numerical modeling at location 66 and its propagation was monitored in the numerical model.

(16) The received signal from numerical modeling is shown as a plot 70 in FIG. 4. The received signal plot 70 clearly indicates two signal breaks, a first signal break 72 at 129 nanoseconds (ns) and a second break 74 at 180 ns. The time when the first break 72 occurs correlates with the distance in the numerical model between transmitter location 66 and receiver location 68 (10 m in this case), while the second break 74 correlates with the length of the path an emitted ground penetrating radar electromagnetic pulse would travel if it was being diffracted by a fracture tip as if it were a secondary behaving as described above with regard to FIG. 2.

(17) To confirm that the second break travel time corresponds to the fracture length, second and third numerical modeling scenarios were simulated with fracture 62 having lengths of 10 m and 15 m, respectively. The received second breaks 76 (FIG. 5) and 78 (FIG. 6) show that as fracture length increases, the time at which the second break occurs becomes greater (FIGS. 5 and 6). Further, the indications of travel time by second breaks 76 and 78 are in agreement with what it would be if the pulse was propagating as a result of secondary electromagnetic waves as illustrated schematically in FIG. 2. The results in FIGS. 4, 5 and 6 of numerical modeling simulations suggest that the fracture tip 62 is indeed acting as a secondary electromagnetic source. These results confirm that there is a relationship between the length of the fracture and the time at which the second break of the emitted ground penetrating radar electromagnetic wave pulse occurs.

(18) To investigate whether the present invention methodology is sensitive at lower frequencies, the numerical modeling scenario for a 10 m long fracture was repeated but this time the pulse was centered at 100 MHz, rather than 2 GHz. Results are plotted in FIG. 7. The received signal clearly shows a first break 80 corresponding to the distance between transmitter and receiver, and a second break 82 related to the fracture length. A notable difference between using 2 GHz and 100 MHz, aside from the amplitude, is the lack of ringing in the 100 MHz signal. This is apparently explained by the fact that at 2 GHz the ground penetrating radar signal can in effect penetrate or see inside the fracture, while at 100 MHz the ground penetrating radar signal cannot penetrate or see inside the fracture because the wavelength at the lower 100 MHz frequency is not small enough. In both cases, however, the present invention is capable to determine the fracture length with ground penetrating radar.

(19) A summary of the results obtained by numerical modeling and plotted in FIGS. 4, 5, 6 and 7 is contained in Table 1.

(20) TABLE-US-00001 TABLE 1 Fracture Path length Frequency Direct length Calculated Simulation (m) (GHz) shot (ns) (m) time (ns) time (ns) 5 2 129 14.1 180 ~180 10 2 129 22.4 290 ~290 10 0.1 129 22.4 290 ~290 15 2 129 31.6 410 ~410

(21) Table 1 above shows that second breaks from these simulations are in accord with the time a second break should occur if indeed the tip of the fracture is acting as a secondary source. This indicates that the present invention can be used to estimate the length of a fracture that is perpendicular to the wellbore.

(22) To demonstrate of the present invention, a further set of numerical simulations were performed. In the numerical simulations, the transmitter position 66 and receiver position 68 were located at different spacings from each other in the wellbore for a fixed fracture length (10 m) and with a 2 GHz ground penetrating radar pulse emitted. Initially, results were obtained with the transmitter position 66 and receiver position 68 located 10 m apart from each other in the numerical model, which is the configuration for the modeling depicted in FIG. 7. The results of numerical modeling of 10 m spacing are shown in the travel times displayed in FIG. 5.

(23) For the second and third scenario of the further set, both the transmitter position 66 and receiver position 68 were moved closer to the fracture 60. First the transmitter position 66 and receiver position 68 were placed 3 m away from the fracture 60 on either side. The results of numerical modeling for 6 m spacing between the transmitter position 66 and receiver position 68 are shown in the travel times displayed in FIG. 8. The transmitter position 66 and receiver position 68 were then moved to 1 m away from the fracture 60 on either side. The results of numerical modeling of 2 m spacing between the transmitter position 66 and receiver position 68 are shown in the travel times displayed in FIG. 9. The received signals plotted in FIGS. 8 and 9 again show that the time when the second breaks 84 (FIG. 8) and 86 (FIG. 9) occur correlates with the length of the path the signal would have to travel if the fracture tip 64 is acting as a secondary source.

(24) Table 2 presents a summary of the results obtained in the further set of simulations.

(25) TABLE-US-00002 TABLE 2 Distance Fracture Fre- between Direct Path length quency Tx- shot length Calculated Simulation (m) (MHz) Rx (m) (ns) (m) time (ns) time (ns) 10 100 10 129 22.4 290 ~290 10 100 6 77 20.9 270 ~270 10 100 2 26 20.1 260 ~260

(26) The numerical modeling results plotted in FIGS. 7, 8 and 9 and summarized in Table 2 indicate that the present invention is sensitive to spacing of ground penetrating radar transmitter and receiver position from the fracture.

(27) Next, further numerical modeling simulations were performed to study the effect of ground penetrating radar frequency on the ability to map fractures with different lengths. Fractures with three different lengths (5 m, 10 m, and 15 m) were simulated and surveyed using 10 MHz and 100 MHz pulses for each of three different fracture lengths.

(28) FIGS. 10A, 11A and 12A are plots of the numerical modeling of fractures mapped with ground penetrating radar emitting a pulse at 20 MHz and fracture lengths of 5 m, 10 m and 15 m, respectively. FIGS. 10B, 11B and 12B are plots of the numerical modeling of fractures mapped with ground penetrating radar emitting a pulse at 100 MHz and fracture lengths of 5 m, 10 m and 15 m, respectively.

(29) The numerical modeling results plotted in FIGS. 10A through 12B show that the lower the ground penetrating radar frequency, the longer the fracture needs to be to be able to properly map its length. While lower ground penetrating radar frequencies offer increased penetration, a tradeoff between penetration and ability to map a fracture's length exists.

(30) Additionally, the effect of the conductivity of the reservoir in the received signal was analyzed. FIGS. 13A, 13B and 13C are plots of numerical models of received ground penetrating radar pulses as a function of time for a fracture having a length of 10 m. FIG. 13A is plot of a received signal with reservoir conductivity assumed to have a zero value. FIG. 13B is plot of a received signal with a reservoir conductivity having a value of 0.003 S/m, or Siemens per meter. FIG. 13C is plot of a received signal with a reservoir conductivity having a value of 0.03 S/m. FIGS. 13A through 13C confirm that higher conductivity in a reservoir formation results in increased attenuation of the ground penetrating radar waves. Therefore the present invention is more preferably suited mapping fractures in low conductivity reservoir formation environments.

(31) The present invention with ground penetrating radar electromagnetic waves offers advantages in contrast to prior seismic waves to map fractures, in addition to electromagnetic waves and seismic waves being different types of wave energy. The present invention does not require monitoring wells in proximity to the wellbore (typically within hundreds of meters for microseismic monitoring). The present invention is also less complicated because it does not require a receiver array at the surface, and the data is less complex to interpret. The present invention can be used for deep reservoirs as opposed to microseismic monitoring. As has been mentioned, microseismic monitoring is limited to shallow reservoirs (typically no deeper than 3 km) due to signal attenuation. The present invention is also not affected by acoustic noise. It can be seen that the present invention provides an improved method to map hydraulic fractures.

(32) The invention has been sufficiently described so that a person with average knowledge in the field of reservoir modeling and simulation may reproduce and obtain the results mentioned herein described for the invention. Nonetheless, any skilled person in the field of technique, subject of the invention herein, may carry out modifications not described in the request herein, to apply these modifications to a determined structure and methodology, or in the use and practice thereof, requires the claimed matter in the following claims; such structures and processes shall be covered within the scope of the invention.

(33) It should be noted and understood that there can be improvements and modifications made of the present invention described in detail above without departing from the spirit or scope of the invention as set forth in the accompanying claims.