AU2013204299B2 - System and method to estimate compressional to shear velocity (Vp/Vs) ratio in a region remote from a borehole - Google Patents

Abstract

In some aspects of the disclosure, a method for creating three-dimensional images of non-linear properties and the compressional to shear velocity ratio in a region (130) remote from a borehole using a conveyed logging tool (140) is disclosed. In some aspects, the method includes arranging a first source (105) in the borehole (110) and generating a steered beam of elastic energy at a first frequency; arranging a second source (115) in the borehole (110) and generating a steerable beam of elastic energy at a second frequency, such that the steerable beam at the first frequency and the steerable beam at the second frequency intercept at a location (130) away from the borehole (110); receiving at the borehole (110) by a sensor (135) a third elastic wave, created by a three wave mixing process, with a frequency equal to a difference between the first and second frequencies and a direction of propagation towards the borehole (110).

Description

WO 2010/121202 PCT/US2010/031490 SYSTEM AND METHOD TO ESTIMATE COMPRESSIONAL TO SHEAR VELOCITY (Vp/Vs) RATIO IN A REGION REMOTE FROM A BOREHOLE CROSS-REFERENCE TO RELATED APPLICATIONS This application claims benefit from U.S. Patent Application Serial No. 5 12/463,796 filed on May 11, 2010, which claims priority to U.S. Provisional Patent Application number 61/170,070 filed on April 16, 2009, both of which are incorporated herein by reference in their entirety, under 35 U.S.C. § 119 (e). GOVERNMENT RIGHTS This invention was made with Government support under Cooperative 10 Research and Development Agreement (CRADA) Contract Number DE-AC52 06NA25396 awarded by the United States Department of Energy. The Government may have certain rights in this invention. FIELD The present invention relates generally to seismic interrogation of rock 15 formations and more particularly to creating three-dimensional images of non-linear properties and the compressional to shear velocity ratio in a region remote from a borehole using a combination of sources in a borehole configured to provide elastic energy, and receiving and analyzing a resultant third wave formed by a three wave mixing process. 20 BACKGROUND Acoustic interrogation of subsurface features tends to be limited by the size and power of practical sources, and in practice, the output of down hole acoustic transducers is limited by the power transmission capabilities of the wireline cable. High frequency signals have a relatively short penetration distance, while low 25 frequency signals generally require large sources, clamped to the borehole wall, to maximize energy transfer to the formation and minimize unwanted signals within the well bore. Currently, acoustic borehole tools are designed with acoustic sources in the borehole to detect returning acoustic waves that are propagating along the borehole walls or scattered by inhomogeneities of linear properties of rock formations 1 surrounding the borehole. U.S. Patent No. 7,301,852 by Leggett, III el al. discloses a Logging While Drilling tool, designed to detect rock formation boundaries, with two acoustic source arrays emitting from a borehole generating a third wave by assumed non-linear mixing in rocks at the location of intersection of the acoustic signals. The 5 third wave is scattered by heterogeneities in subsurface properties, and the scattered signal is detected by sensors in the logging tool. The source arrays are merely disclosed to be directional without any further description. Attempts have been made to characterize the non-linear properties of a formation in the area of oil and gas prospecting from boreholes, but each has its own 10 limitations. For example, U.S. Patent No. 5,521,882 by D'Angelo et al. discloses an acoustic tool designed to record with pressure receivers the non-linear waves that propagate along the borehole wall with limited penetration into the surrounding rock formation and refract back into the well bore fluid. U.S. Patent No. 6,175,536 by Khan discloses a method to estimate the degree of non-linearity of earth formations 15 from spectral analysis of seismic signals transmitted into the formation from one borehole and received in a second borehole. In light of these prior attempts, there is a need for an apparatus and method for generating three-dimensional images of non linear properties in a remote region from a borehole. Reference to any prior art in the specification is not, and should not be taken 20 as, an acknowledgment or any form of suggestion that this prior art forms part of the common general knowledge in Australia or any other jurisdiction or that this prior art could reasonably be expected to be ascertained, understood and regarded as relevant by a person skilled in the art. SUMMARY 25 In accordance with some aspects of the disclosure, there is provided a method for estimating compressional to shear velocity (Vp/Vs) ratio in a region remote from a borehole using a conveyed logging tool, the method comprising: arranging a first source in the borehole and generating elastic energy at a first frequency; arranging a second source in the borehole and generating elastic energy at a second frequency, 30 wherein the elastic energy at the first and the second frequency intersect at a mixing zone located away from the borehole; receiving, by an array of sensors, at the 2 borehole a third elastic wave, created by a three wave non-collinear mixing process in a non-linear mixing zone, with a frequency equal to a difference between the first and second frequencies and a direction of propagation towards the borehole; identifying a location of the three wave mixing zone based on the arrangement of the first and the 5 second sources, the direction of the third wave and selection rules governing non collinear mixing in acoustically non-linear media; and estimating Vp/Vs ratio using data recorded by repeating the generating, receiving and identifying at a plurality of azimuths, inclinations and longitudinal locations of the first source, the second source, or both within the borehole. 10 In accordance with some aspects of the disclosure, there is provided a computer implemented method arranged to perform processing of a third wave generated through a three wave mixing process of two non-collinear primary beams or waves in a remote non-linear region of a formation surrounding a borehole, and recorded in the borehole, in order to estimate Vp/Vs ratio of a formation surrounding the borehole, 15 the computer-implemented method comprising: analyzing a frequency content of a recorded third wave by spectral analysis and selecting a signal corresponding to a difference frequency to isolate the third wave signal generated by the non-linear mixing process; analyzing an amplitude of the recorded third wave as a function of frequency ratios of the primary mixing beams or waves and determining a mixing 20 location where the third wave signals originated, from selection rules of non-collinear mixing in non-linear media, wavenumbers of the first and second beams or waves and the third wave and locations of two primary sources and a sensor array; confirming a direction of the third wave impinging on the borehole and recorded the sensor array based on analyzing the third elastic wave detected at the sensory array; and estimating 25 the Vp/Vs ratio of the zone of interaction of the two non-collinear primary beams or waves from the amplitude and direction of the third elastic wave arrivals. In accordance with some aspects of the disclosure, a method for creating three-dimensional images of non-linear properties in a region remote from a borehole using a conveyed logging tool is disclosed. The method includes arranging with a 30 specific spatial configuration a first source in the borehole and generating a primary wave of elastic energy at a first frequency; arranging a second source in the borehole and generating a primary wave of elastic energy at a second frequency; receiving by 3 an array of three component sensors at the borehole the arrival of the third elastic wave created by a three wave mixing process, with a frequency equal to a difference between the first and second primary frequencies, that propagates back to the borehole; determining the propagation direction of the third wave from the signals 5 received by the sensor array; imaging the locus of the three wave mixing region based on the arrangement of the first and second sources and the propagation direction of the third wave; and creating three-dimensional images of the non-linear properties using data recorded by repeating the generating, receiving, determining and imaging steps at a plurality of azimuths, inclinations and longitudinal locations within the 10 borehole. In accordance with some aspects of the disclosure, further methods for creating three-dimensional images of non-linear properties in a region remote from a borehole using a conveyed logging tool are disclosed. These share the common configuration of two sources and an array of sensors in the borehole, but differ in that 15 the one or other of the sources may generate a steerable beam or a wave of elastic energy, and the sensor units in the array may be a combination of non-directional and three component devices. The method includes arranging with a specific spatial configuration a first source in the borehole and generating either a steerable primary beam of elastic energy or a primary wave of elastic energy at a first frequency; 20 arranging a second source in the borehole and generating either a steerable primary beam of elastic energy or a primary wave of elastic energy at a second frequency, such that the energy from the two sources mixes at locations away from the borehole; receiving by a sensor array at the borehole the direct arrival of the third elastic wave, created by a three wave mixing process, with a frequency equal to a difference 25 between the first and second primary frequencies, that propagates back to the borehole in a specific direction; locating the three wave mixing region based on the arrangement of the first and second sources and on properties of the third wave signal; and creating three-dimensional images of the non-linear properties using data recorded by repeating the generating, receiving and locating steps at a plurality of 30 azimuths, inclinations and longitudinal locations within the borehole. In accordance with some aspects of the disclosure, three dimensional images of the non-linear properties of the formations surrounding the borehole are 4 transformed to reservoir properties using appropriate relations between formation non-linearity and said properties. The images may be of properties at the time of logging, or may represent changes between two logging runs separated by the passage of time. 5 In accordance with some aspects of the disclosure, methods to create three dimensional images of the ratio of compressional to shear acoustic velocity of rocks surrounding the borehole are disclosed. These methods are variations of the methods for creating three dimensional images of non-linear properties discussed in paragraphs [0005] to [0008], above. 10 In accordance with some aspects of the disclosure, an apparatus for creating three-dimensional images of non-linear properties and the compressional to shear velocity ratio of the rock formations remote from a borehole using a conveyed logging tool is disclosed. The apparatus includes a first source arranged in the borehole and configured to generate a steerable beam or a wave of elastic energy at a 15 first frequency; a second source arranged in the borehole and configured to generate a steerable beam or a wave of elastic energy at a second frequency, such that the beams or waves at the first frequency and the second frequency intercept at a location away from the borehole; and a non-directional or three component sensor array configured to receive a third elastic wave if that the non-linear properties of the region of interest 20 result in the creation of the third elastic wave by a three wave mixing process having a frequency equal to a difference of the first and the second frequencies and a specific direction of propagation back to the borehole; a first processor arranged in the borehole to control source firing and recording of the third elastic wave; a device configured to transmit data up-hole through a wireline cable for a wireline tool; and a 25 second processor arranged to create the three-dimensional images based, in part, on properties of the received third wave and the arrangement of the first and second sources. These and other objects, features, and characteristics of the present invention, as well as the methods of operation and functions of the related elements of structure 30 and the combination of parts and economies of manufacture, will become more apparent upon consideration of the following description and the appended claims 5 with reference to the accompanying drawings, all of which form a part of this specification, wherein like reference numerals designate corresponding parts in the various Figures. It is to be expressly understood, however, that the drawings are for the purpose of illustration and description only and are not intended as a definition of 5 the limits of the invention. As used in the specification and in the claims, the singular form of "a", "an", and "the" include plural referents unless the context clearly dictates otherwise. As used herein, except where the context requires otherwise, the term "comprise" and variations of the term, such as "comprising", "comprises" and 10 "comprised", are not intended to exclude further additives, components, integers or steps. BRIEF DESCRIPTION OF THE DRAWINGS Figure I shows a configuration for creating three-dimensional images of non linear properties in a region remote from a borehole in accordance with various 15 aspects of the disclosure Figure 2 shows a configuration for creating three-dimensional images of non linear properties in a region remote from a borehole in accordance with aspects of the disclosure. Figure 3 shows a configuration for creating three-dimensional images of non 20 linear properties in a region remote from a borehole in accordance with aspects of the disclosure. Figure 4 shows a flow chart for creating three-dimensional images of non linear properties in a region remote from a borehole in accordance with various aspects of the disclosure. 5A WO 2010/121202 PCT/US2010/031490 Figures 5a, 5b and 5c shows a numerical simulation of selection rule 1 for a beam-beam interaction of Table 1 when the two primary waves are beams. Figure 6 illustrates the geometry of the generation of the difference frequency third wave by non-linear mixing of two primary acoustic waves as governed by the 5 non-linear mixing selection rule. Figure 7 shows an application of aspects of the disclosure for imaging using a beam and broad beam or plane wave. DETAILED DESCRIPTION Figure 1 shows one of several possible configurations for creating three 10 dimensional images of non-linear properties and the compressional to shear velocity ratio in a region remote from a borehole in accordance with various aspects of the disclosure. First source 105 is arranged in borehole 110 to generate a steerable primary beam of acoustic energy at a first frequency fi. Second source 115 is also arranged in borehole 110 to generate a steerable primary beam of acoustic energy at a 15 second frequency f2. By way of a non-limiting example, both first source 105 and second source 115 may be a phased array of sources and may be configured to generate either compressional or shear steerable beams. As shown in Figure 1, first source 105 is arranged on first tool body 120 and second source 115 is arranged on second tool body 125. However, the disclosure is 20 not so limiting as first tool body 120 and second tool body 125 may also be arranged together on a common tool body (not shown). Tool bodies 120 and 125 are arranged to be independently moveable within bore hole 110 in at least two degrees of freedom including translation along the longitudinal axis 150 of borehole 110 and rotation 155 in azimuth about the longitudinal axis of borehole 110. First source 105 may be 25 arranged above or below second source 115 in borehole 110. Tool bodies 120 and 125 may be arranged on a conveyed logging tool (not shown) within borehole 110. For a given azimuth orientation of first source 105 and second source 115, the beam generated by second source 115 and the beam generated by first source 105 are configured such that the beams converge and intercept in a mixing zones 130 remote 30 from borehole 110. By a combination of independently steering the beams and 6 WO 2010/121202 PCT/US2010/031490 changing the separation between the sources 105, 115, the mixing zones 130 move in the plane defined by the beams and the longitudinal borehole axis 150, while controlling the angle of interception. The distance of mixing zones 130 from borehole 110 can range from near the edge of borehole 110 to about 300 meters into the 5 surrounding subsurface rock formation. By way of a non-limiting example, the phase difference and/or time delays between adjacent elements in the source array referred to in paragraph [0018] may be modified to focus the acoustic energy of the primary beams at a particular mixing zone. The non-linear properties of the earth at the location between the two waves 10 result in the generation of a third elastic wave. The third elastic wave is a result of a three-wave mixing process that occurs in nonlinear materials, in this case, rock formations. In this process, two converging non-collinear waves of different frequencies, fi andf 2 , also called primary waves, mix to form additional waves at the harmonic and intermodulation frequencies f, - f2, fi + f2, 2xfj and 2xf 2 , etc. The 15 strength of the third wave is a function the non-linearity of the rocks in the mixing zones. By way of a non-limiting example, when a primary compressional (P) wave with a frequencyfl and a primary shear (SV) wave with a frequencyf 2 cross in a non linear medium, a third compressional (P) or shear (SV) wave is generated with a frequencyfi -f2. For further description, see Johnson et al. (1987) and Johnson and 20 Shankland (1989), which is hereby incorporated by reference in its entirety. As further discussed in paragraphs [0031] and [0033] below under propagation selection rules, the third wave propagation vector is co-planar with the propagation vectors of the two primary waves. Certain combinations of angle of intersection, fi/f2 ratio and compressional to shear velocity ratio result in a third elastic 25 wave with frequency fj - f2 propagating in a specific angle relative to the primary beams back to the borehole 110. Sensor or receiver array 135 is arranged at specific location in borehole 110 to detect the third wave returning to the borehole 110. In some aspects of the disclosure, as shown in the Figure, sensor array 135 comprises more than one sensor arranged as 30 an array of sensors on sensor tool body 140 and separate from tool bodies 120 and 125. Sensor 135 is configured to be independently moveable within bore hole 110 along the longitudinal axis 150 of borehole 110. In some aspects, sensor tool body 7 WO 2010/121202 PCT/US2010/031490 140 is arranged below tool bodies 120 and 125 or arranged above and below tool bodies 120 and 125. In some aspects, sensor tool body 140 is connected to either or both tool bodies 120 and 125. The third wave is detected at borehole 110 by sensor array 135. Figure 2 5 shows an arrangement similar to Figure 1, wherein receiver 135 is replaced by three component geophone 145 clamped to the borehole walls. The resultant signal is decomposed by processing into its inclination and azimuth in order to add redundancy to the system by determining the direction of the incoming third wave arrival. In some aspects, a first processor configured to execute machine-readable 10 instructions (not shown) may be arranged in borehole 110 to perform various processing tasks, such as controlling source firing and compressing or filtering the data recorded by sensor array 135. A second processor configured to execute machine-readable instructions (not shown) may be arranged outside borehole 110 to assist the first processor or perform different processing tasks than the first processor. 15 For example, the second processor may perform part or all processing activities in creating the three-dimensional images. A transmitter or transceiver (not shown) may be arranged in borehole 110 to transmit data up-hole through a wireline cable (not shown). At a given depth along the borehole of one of the sources 105, 115, sweeping 20 the beams in inclination at constant relative bearing to spatially scan the mixing zone in a plane passing through the borehole axis, rotating the sources azimuthally to rotationally scan the mixing region and moving the whole assembly along borehole 110, results in scanning a 3D volume of mixing zones around the borehole for non linear properties. With sources 105, 115 and sensor array 135 located on independent 25 tool bodies, high redundancy in the data can be obtained and the depth of investigation can be varied. In this way, a 3D volume of the rocks surrounding the borehole can be interrogated for non-linear properties and a 3D image of non-linear properties can be processed and computed from the returned signals. Figure 3 shows another arrangement for creating three-dimensional images of 30 non-linear properties in a region remote from a borehole in accordance with various aspects of the disclosure. The arrangement of Figure 3 is similar to the arrangement 8 WO 2010/121202 PCT/US2010/031490 in Figure 2, with the primary difference being that the sources are arranged in borehole 110 to produce elastic waves instead of steerable beams. With reference to Figure 3, first source 305 is arranged in borehole 110 on first tool body 320 to generate a first elastic wave of acoustic energy at a first frequencyfy. Second source 5 315 is arranged in borehole 110 on second tool body 325 to generate a second elastic wave of acoustic energy at a second frequency f2. First and second elastic waves produced by sources 305, 315 are arranged to intercept away from borehole 110 at various mixing zones 130. Receiver 145 is arranged within borehole 110 to receive a third wave that is produced in the mixing zones 130 by the three-wave mixing process 10 discussed above, and further discussed below. Since the waves produced by sources 305, 315 are essentially non-directional, mixing between the waves occurs simultaneously in the entire area of mixing zones 130, that also extends out of the plane of the Figure, and receiver 145 tends to have directional characteristics. By way of a non-limiting example, a three component geophone array may be used for this 15 purpose. The resultant signal is decomposed by processing into multiple arrival signals at a range of inclinations and azimuths and travel times. Given the locations of sources and the receivers, the travel times and directions of each decomposed directional arrival, there is sufficient information to apply the selection rules described below in paragraphs [0031] to [0033] to determine a unique mixing zone where the 20 third wave was generated. This unique mapping allows the construction of a three dimensional image from the properties of the received signal. Figure 4 shows a method for creating three-dimensional images of non-linear properties and the compressional to shear velocity ratio in a region remote from a borehole using a conveyed logging tool. The method begins at 405 where a first 25 source is arranged in the borehole to generate a steerable beam elastic energy at a first frequency and a second source is arranged in the borehole to generate a steerable beam of elastic energy at a second frequency. The steerable beams at the first and second frequency are arranged to intercept at a location away from the borehole. As such, the second beam is generated at the same azimuth as the first beam, but at a 30 different inclination relative to the longitudinal axis of the borehole. The method continues at 410 where a third elastic wave is received at the borehole by a sensor array. As discussed above, the third elastic wave is created by a three wave mixing process, with a frequency equal to a difference between the first and second 9 WO 2010/121202 PCT/US2010/031490 frequencies and a direction of propagation towards the borehole. At 415, a three wave mixing location away from the borehole is determined from the arrangement of the first and second sources and properties of the third wave, by recourse to the selection rules discussed below in paragraphs [0031] to [0033]. At 420, three-dimensional 5 images are created of the non-linear properties using data recorded by repeating the generating of step 405, the receiving of step 410 and the determining of step 415 at a plurality of azimuths, inclinations and longitudinal locations within the borehole. In cases of compressional-shear interaction the received signals are analyzed in step 425 for the compressional/shear velocity (Vp/Vs) ratio according to paragraph [0038]. At 10 430, the non-linear properties are transformed to physical reservoir properties such as fluid saturation, effective stress, fracture density and mineralogy. In some aspects of the present disclosure, the first and second sources may be beam or cylindrical or spherical wave sources, and the sensor array may be any combination of non-directional single component sensors and three component 15 geophones. Alternative permutations of the component parts offer different degrees of redundancy in signal processing and imaging. Experimental demonstrations of non-linear mixing of two acoustic waves in solid have been reported, for example by Rollins, Taylor and Todd (1964), Johnson et al. (1987) and Johnson and Shankland (1989), which are hereby incorporated by 20 reference in their entirety. In the special case where a primary compressional (P) wave with a frequency fi and a primary shear (S) wave with a frequencyf 2 cross each other, in a non-linear medium, a third P or S wave is generated with the frequency f2. If the primary P and S waves are beams with wave vectors k 1 and k 2 , respectively, and the non-linear formation property is uniform, the kinematics of wave interaction 25 requires the resulting third wave to be a plane wave with wave vector k 3 that obeys the selection rule k 1 - k 2 = k 3 . The selection rule imposes a very tight restriction on the permissible crossing angles for the primary waves and a specific propagation direction of the third wave. The general kinematic theory for non-linear mixing of two linear plane waves and the selection rules and amplitude responses have 30 contributions from Jones and Kobett (1963), Rollins, Taylor et al. (1964) and later by Korneev, Nihei and Myer (1998), all of which are hereby incorporated by reference in their entirety, who also provide specific relationships between non-linear parameters 10 WO 2010/121202 PCT/US2010/031490 of the mixing medium and the non-linear mixing signal strength. For example, Equation 53 and 54 of Korneev, Nihei and Myer show that the mixing strength of P and SV (vertically polarized shear) plane waves is proportional to a specific combination of non-linear parameters of the rocks. 5 The selection rules of Komeev, Nihei and Myer governing the nonlinear interaction of two elastic plane waves can be used as guidance for the interaction of two elastic beams. These plane wave selection rules dictate that the following six nonlinear interactions produce backscattered waves. Table I - Selection Rules Governing Non-Linear Interaction of Two Elastic 10 Plane Waves. In this table, and elsewhere in this document,fi is greater thanf 2 . Selection Rules 1" beam or wave 2 beam or wave Resultant 3d beam or wave from 1 't + 2 nd 1 P(f 1 ) SV(f 2 ) P(f; -f2) 2 P(f) SV(f2) SV(fl -f2) 3 P(fl) SH(f 2 ) SH(fl -f2) 4 P(1) SV(f 2 ) P(f +f 2 ) 5 SV(fl) SV(f2) P(fz +f 2 ) 6 SH(fl) SH(f 2 ) P(fl +f) Figures 5a, 5b and 5c shows a numerical simulation of selection rule 1 of Table 1 when the two primary waves are beams of a beam-beam interaction. A 25 kHz compressional beam, shown in Figure 5a, and a 18 kHz shear beam, shown in Figure 5b, mix to form a third beam, shown in Figure 5c, with frequency 7 kHz = 25 15 kHz - 18 kHz. In this example, in accordance with the plane wave predictions of Korneev, Nihei and Myer, a third back propagating P beam with frequency (fl -f2) at an angle of 133' to the P(fz) wave is generated by nonlinear mixing in the region where the P(fl) and SV(f2) beams overlap. 11 WO 2010/121202 PCT/US2010/031490 The kinematics of non-linear interactions of beams results in the generation of specific combinations of wave vectors and frequencies. The third wave returns at a specific travel time, and with specific frequenciesf 3 and wave vectors k 3 such as fs=fl - f2 and k 3 =kl - k 2 . For a combination offi, f2, k 2 and k 3 , there is a well-defined 5 propagation wave vector k 3 of the third wave in the same plane, defined by k 1 and k 2 . There is a direct correspondence between the signal detected at a particular receiver position and the location where the non-linear mixing of the two primary waves k, and k 2 takes place. The signal strength of the receiver would be proportional to the strength of the non-linearity of the rocks in the mixing zone, among other factors, and 10 reach a maximum for a receiver lying on vector k 3 . Therefore, the signal strength at the receivers can be geometrically mapped onto the non-linearity of the rocks along the beam trajectory as indicated by Figure 1. The geometrical theory of wave propagation indicates that the beam generated in each interaction zone would arrive at the borehole at a specific receiver defined by 15 the geometry of the three wave vectors k 1 , k 2 and k 3 , after a specific time delay. The strength of the returning signal at a specific location in the borehole at a particular time is dependent on the degree of non-linearity of the interaction location, and hence a time image of the relative strength of the non-linear properties of the rocks along the beam can be constructed. The amplitude magnitude of a returned signal at the 20 receivers is itself indicative of certain petrophysical properties of the mixing zone. If the beam and plane wave are scanned in azimuth and inclination while preserving the necessary convergence angle, a localized circumferential and radial 3D image of non linear properties of rocks surrounding the borehole can be obtained. By moving the entire assembly up and down the borehole, repeated 3D images of non-linear 25 properties of rocks surrounding the borehole are obtained. By making weighted stacks of these repeated images, a final image of non-linear properties of rocks surrounding the entire borehole can be constructed through subsequent computer processing. In addition, if the sources and the receivers are part of three separate tool bodies, one or two can be moved while the third one is fixed (for example, the sources 30 are fixed while the receiver tool body is moved up and down). Alternatively, several descents into the well may be made with different spacing between the tool bodies. 12 WO 2010/121202 PCT/US2010/031490 For non-linear mixing between an elastic beam and a broader beam (quasi plane wave), the selection rule is relaxed. Third waves of frequencyfi -f2, centered around the wave vector k 3 = k 1 - k 2 , are generated continuously along the primary beam if the beam width is about ten wavelengths of the third wave. The resulting 5 signal strength forfi =f; -f2 is a function of the average non-linear properties of the mixing region, the average ratio of velocity off, propagation and average velocity for f 2 propagation (noting thatfi andf 2 may be compressional or shear), the volume of the mixing zone and the geometry of the mixing. This function can be computed for various mixing modes. For example, the signal strength for a particular important 10 mixing mode such as compressional wave P forfi and SV forf 2 is given by U = 2) 2 1 2PSP A 1

B

2 fl 1 r2 )fVPSAP)FPSvP PS'P VP2V r where U is the displacement amplitude of the third wave received at the borehole, A, is the longitudinal polarization of the compressional wave and B 2 is the transverse polarization of the shear wave. 8 is a function of the A, B and C parameters of Landau 15 and Lifschitz representing the non-linearity of the rocks in the mixing zone. v is the volume of the mixing zone, r is the distance from mixing zone to the receiver. F is the geometric form factor of order 1 which is dependent on the geometry of the incident beams and can be numerically computed from the Korneev, Nihei, Myers theory for the particular geometry. A is a selection rule form factor which is a 20 numerically computable function of the wave vectors ki, k 2 and k 3 and is only significant if the interaction geometry honors the selection rules. The subscript PSvP in the formula refers to compressional-shear interaction generating a compressional wave. In accordance with certain aspects of this disclosure, an image of the 25 compressional to shear velocity ratio may be constructed as follows. When one of the sources generates a compressional wave (P-wave) with frequency fi and the other source generates an SV-wave with frequencyf2 and both waves are steered towards a specific mixing volume, the propagation direction of the third compressional wave (P wave) with difference frequencyfs =fl -f2 is controlled by the average in situ Vp/Vs 30 ratio of the rock in the mixing zone as governed by the selection rules as shown in Figure 6. From the measurements of the signal in the three component receiver array 13 WO 2010/121202 PCT/US2010/031490 145 on Figure 2 or Figure 3, the direction of this third wave can be determined and thereby, the in situ Vp/Vs of the mixing zone can be computed. If the beam and plane wave are scanned in azimuth and inclination while preserving the necessary convergence angle, a localized circumferential and radial 3D image of in situ Vp/Vs 5 ratio of rocks surrounding the borehole can be obtained. By moving the entire assembly up and down the borehole, repeated 3D images of in situ Vp/Vs of rocks surrounding the borehole may be obtained. By making weighted stack of these repeated images, a final image of in situ Vp/Vs of rocks surrounding the entire borehole can be constructed through subsequent computer processing. Alternatively, 10 several descents into the well may be made with different fixed spacing between the tool bodies. In some aspects of this disclosure, an alternative determination of Vp/Vs ratio is achieved through scanning the ratio of the frequencies fi to f2 of the primary beams. Figure 6 illustrates the geometry of the interaction of two beams such as those 15 generated in the configuration of Figure 1, that may be analyzed using the vector mathematics and trigonometry described above in paragraphs [0031] to [0034]. The lengths k, and k 2 of vectors k 1 and k 2 are defined by the ratio of their corresponding frequencies and velocities. As shown in Figure 6, the returning angle (p is a function of ft/f2, Vp/Vs ratio and the intersection angle 0 of the two primary beams. In 20 addition, the physical selection rules only permit the generation of a third wave at specific combinations of f/f2, Vp/Vs ratio and angle of interception 0, such as the example illustrated on Figure 5. Using the symbol r for the Vp/Vs ratio and the terms defined on Figure 6, the magnitude k 3 of vector k 3 is given by the vector sum of ki and -k 2 , that is 25 k |- k| and also by the cosine rule that states k = 191- kj- 2kzka mes. Combining the two equations, and substituting fi/Vp for k, andf 2 /Vs for k 2 , leads to a statement of the geometric conditions imposed by the selection rules. The quadratic equation r -2 cs; f - + 2 = a may be solved for r, the VpVs ratio of the mixing zone. This leads to a non-limiting alternative 30 method for measuring in situ Vp/Vs ratio of a particular mixing region by the 14 WO 2010/121202 PCT/US2010/031490 following sequence: a) record a standard sonic waveform log to determine Vp and Vs near the wellbore to acquire data to estimate the phase differences between adjacent elements in a phased source array to steer the beams at the approximate convergence angle for the geometry of the planned measurement; b) steer the P and SV sources to 5 converge at a controlled angle 0 and mix at a particular region in space surrounding the borehole; c) vary f2 while fixing fi and measure the amplitude of the received signal at the difference frequency fi -f2 at the sensors in the borehole; d) identify the frequency at which the signal each receiver in the array reaches a maximum amplitude strength; and e) determine angles 0 and p from the geometry of the sources 10 and receivers. By sweeping the beams in inclination, rotating in azimuth, and moving the entire assembly up and down the borehole and repeating the above procedure, the VpVs ratio of a 3D volume around the borehole is interrogated and thereby 3D images of in situ Vp/Vs ratio of rocks surrounding the borehole may be obtained. The methods described above offer an advantageous property in that the 15 frequency differences -f2 is very specific, allowing for spectral analysis to enhance the signal to noise ratio of the measurements. Moreover, if both frequenciesfi andf 2 are simultaneously chirped proportionally, the resulting difference frequency signal fi -f2 would also be a well defined chirped signal. The time-varying code may include one or more of a variation in amplitude, a variation in frequency, and/or a variation in 20 phase of the first, the second, or both the first and the second beams or waves. The third difference wave can be broad band if one of the primary frequencies is swept through a range of frequencies while their frequency ratio is fixed. Thus, the resulting third beam f2 - fi will be swept across a wide frequency range, while preserving the same direction. This allows for improvement in signal to noise by standard auto 25 correlation of the chirped or coded signal. Since the wave vector k 3 = ki - k 2 is well defined, the signal to noise discrimination of the recorded third wave from receivers 135 can be enhanced further by employing three-component receivers in the borehole. The signals from the three components can be tuned to specific directivity by a technique, such as, hodogram 30 analysis. In some aspects of the present disclosure, the signal to noise ratio can be improved by repeating the above steps with an inverse polarity (180 degrees out of 15 WO 2010/121202 PCT/US2010/031490 phase) and adding the results together. The returning difference frequency signal will add coherently as its amplitude is proportional to the product of the amplitudes of the two primary waves and therefore will not reverse polarity when the polarity of the primary source is reversed, while any linear noises generated by the primary waves in 5 the system will reverse polarity and cancel upon addition. Alternative methods can be devised with various non-exclusive combinations of beams and waves. By way of a non-limiting example, a method to generate images by computer processing of acoustic and seismic signals includes the follow steps. First, perform spectral analysis of the frequency content of the recorded third wave 10 and applicable selection rules of the difference frequency signal in order to isolate the third wave signal generated by the non-linear mixing process. In the case that the sensors include three component geophones, determine the direction of the third wave impinging on the borehole using orientation techniques. The method continues by analyzing the amplitude of the recorded third wave as a function of frequency ratios 15 of the primary mixing waves and determining the mixing location where the third wave signals originated, from the selection rules of non-collinear mixing in non-linear media, the wavenumbers of the first and second beams and the third wave and the locations of the two beam sources and the sensor array. The method continues by constructing seismograms determined by cross-correlation of the received signals 20 with chirped transmitter signals for each source-receiver combination. The method continues by performing three dimensional time or depth imaging to the entire data set, in order to obtain three dimensional images of the non-linear properties of the formation surrounding a borehole in either or both of time and distance. The methods for generating images from seismograms are known, for example, Hill et al., which is 25 hereby incorporated by reference, have provided the general methodology for the special case of imaging from beams. Another non-limiting alternative imaging method is illustrated on Figure 7, which shows the case of interactions of a narrow 705 and a broad (wide) beam 710. Given a smooth background model of Vp and Vs of the investigated volume, 30 application of the selection rules enables the geometric mapping of the energy detected at a receiver location 735 on to mixing zones 730 along the narrow beam. A time image of the non-linear property can thus be constructed along the narrow beam. By 16 WO 2010/121202 PCT/US2010/031490 rotating in azimuth and moving the assembly along the borehole, a three dimensional time image can be constructed of a volume centered on the borehole. Successive repetition of the measurement at different beam inclinations, and altering the f2/f, frequency ratio a yields a series of three dimensional time images. This redundancy in 5 imaging permits the further refinement of the smooth background model and a three dimensional spatial image. Non-linear parameters of rocks have been found to be related to a number of important hydrocarbon reservoir parameters, such as variations with gas, oil and water saturation, effective stress, fracture density and mineralogical content. For example, 10 see Ostrovsky and Johnson 2001, which is hereby incorporated by reference. In certain aspects of this disclosure, the 3D images of non-linear properties constructed by this method are transformed to provide quantitative information on the distribution of these properties around the borehole at the time of recording. In addition, sequential repetitions of this method are used to detect changes in reservoir properties 15 over time for reservoir monitoring purposes. The recordings of received waveforms are processed to generate an image of the non-linear characteristics of the formation. The directivity of the beam and the time of flight may fix the locations where scattered waves are generated, distinguishing this device from normal sonic imaging techniques using conventional 20 non-directional monopole and dipole sources. Although the invention has been described in detail for the purpose of illustration based on what is currently considered to be the most practical and preferred embodiments, it is to be understood that such detail is solely for that purpose and that the invention is not limited to the disclosed embodiments, but, on the 25 contrary, is intended to cover modifications and equivalent arrangements that are within the spirit and scope of the appended claims. As a further example, it is to be understood that the present invention contemplates that, to the extent possible, one or more features of any embodiment can be combined with one or more features of any other embodiment. 17

Claims (22)

1. A method for estimating compressional to shear velocity (Vp/Vs) ratio in a region remote from a borehole using a conveyed logging tool, the method comprising: arranging a first source in the borehole and generating elastic energy at a first 5 frequency; arranging a second source in the borehole and generating elastic energy at a second frequency, wherein the elastic energy at the first and the second frequency intersect at a mixing zone located away from the borehole; receiving, by an array of sensors, at the borehole a third elastic wave, created 10 by a three wave non-collinear mixing process in a non-linear mixing zone, with a frequency equal to a difference between the first and second frequencies and a direction of propagation towards the borehole; identifying a location of the three wave mixing zone based on the arrangement of the first and the second sources, the direction of the third wave and selection rules 15 governing non-collinear mixing in acoustically non-linear media; and estimating Vp/Vs ratio using data recorded by repeating the generating, receiving and identifying at a plurality of azimuths, inclinations and longitudinal locations of the first source, the second source, or both within the borehole. 20

2. The method according to claim 1, further comprising creating three dimensional images of the Vp/Vs ratio using data recorded by repeating the generating, receiving and identifying at a plurality of azimuths, inclinations and longitudinal locations within the borehole. 25

3. The method according to claim 1, wherein the first source and the second source are configured to generate elastic energy selected from the group consisting of: steerable beams, waves with limited directionality, and combinations thereof. 18 1000321052

4. The method according to claim 3, wherein the first source and the second source comprise an array of sources.

5. The method according to claim 3, wherein the first source, the second source, 5 and the sensor array are arranged either on a common tool body or on separate tool bodies of the conveyed logging tool.

6. The method according to claim 5, wherein the separate tool bodies can be moved independently along a longitudinal axis of the borehole. 10

7. The method according to claim 3, wherein the sensor array comprises three component sensors clamped to a wall of the borehole.

8. The method according to claim 3, comprising concentrating energy of the 15 wave at the first frequency and/or the wave at the second frequency at a selected azimuth and inclination, relative to a longitudinal axis of the borehole.

9. The method according to claim 3, comprising generating the wave at the second frequency at a range of frequencies f2 of the form f2 = a fl, and sweeping 20 values of a.

10. The method according to claim 9, comprising repeating the concentrating of energy of the wave at the second frequency at a selected azimuth and inclination, relative to a longitudinal axis of the borehole and the generating the wave at the 25 second frequency at a range of frequencies f2 of the form f2= a fl, and sweeping values of a after rotating the second source to a selected azimuth with respect to a longitudinal axis of the borehole. 19 1000321052

11. The method according to claim 9, comprising repeating the concentrating of energy of the wave at the second frequency at a selected azimuth and inclination, relative to a longitudinal axis of the borehole and the generating the wave at the second frequency at a range of frequencies f2 of the form f2 = a fl, and sweeping 5 values of a after rotating the second source to a selected azimuth with respect to a longitudinal axis of the borehole, after moving the sources and/or the sensor array along a longitudinal axis of the borehole.

12. The method according to claim 9, comprising repeating with inverse polarity 10 the concentrating of energy of the wave at the second frequency at a selected azimuth and inclination, relative to a longitudinal axis of the borehole and the generating the wave at the second frequency at a range of frequencies f2 of the form f2 = a fl, and sweeping values of a after rotating the second source to a selected azimuth with respect to a longitudinal axis of the borehole, after moving the sources and the sensor 15 array along a longitudinal axis of the borehole.

13. The method according to claim 9, comprising obtaining a plurality of radial scans by altering spacings between the first and/or the second source and/or the sensor array in order to change a depth of investigation and resolution of the three 20 dimensional images.

14. The method according to claim 9, comprising either chirping or coding or both chirping and coding either or both the acoustic energy transmitted by the first source and the second source. 25

15. The method according to claim 9, comprising modulating an amplitude, a phase, a period or any combination thereof of either the acoustic energy at the first frequency or the acoustic energy at the second frequency. 30

16. The method according to claim 1, further comprising: 20 1000321052 varying one of the first or the second frequency while fixing the other frequency and measuring an amplitude of the received third elastic wave; identifying a frequency at which the received signal reaches a maximum amplitude strength for each receiver in the sensor array; and 5 measuring in situ Vp/Vs ratio of a particular mixing region away from the borehole by determining the angle of intersection of the beams of the first frequency and the beams at the second frequency and a returning angle of the third elastic wave from the geometry of the first and the second sources, the sensor array and selection rules governing non-linear and non-collinear mixing properties. 10

17. A computer-implemented method arranged to perform processing of a third wave generated through a three wave mixing process of two non-collinear primary beams or waves in a remote non-linear region of a formation surrounding a borehole, and recorded in the borehole, in order to estimate Vp/Vs ratio of a formation 15 surrounding the borehole, the computer-implemented method comprising: analyzing a frequency content of a recorded third wave by spectral analysis and selecting a signal corresponding to a difference frequency to isolate the third wave signal generated by the non-linear mixing process; analyzing an amplitude of the recorded third wave as a function of frequency 20 ratios of the primary mixing beams or waves and determining a mixing location where the third wave signals originated, from selection rules of non-collinear mixing in non-linear media, wavenumbers of the first and second beams or waves and the third wave and locations of two primary sources and a sensor array; confirming a direction of the third wave impinging on the borehole and 25 recorded the sensor array based on analyzing the third elastic wave detected at the sensory array; and estimating the Vp/Vs ratio of the zone of interaction of the two non-collinear primary beams or waves from the amplitude and direction of the third elastic wave arrivals. 21 1000321052

18. The computer-implemented method according to claim 17, further comprising imaging, applied to the entire data set, in order to obtain three dimensional images both in time and in distance of the Vp/Vs ratio of the formation surrounding a 5 borehole.

19. The computer-implemented method according to claim 17, further comprising determining a location of the non-linear region where the three wave mixing process took place. 10

20. The computer-implemented method according to claim 17, wherein the first source and the second source are configured to generate elastic energy selected from the group consisting of: steerable beams, waves with limited directionality, and combinations thereof. 15

21. The computer-implemented method according to claim 17, further comprising determining the direction of the third wave impinging on the borehole if the sensors include three component geophones. 20

22. The computer-implemented method according to claim 17, wherein the recorded signals obtained from two consecutive opposite polarity firings of sources are added in order to suppress noise and to enhance the amplitude of the nonlinear third wave. 22