WO2013177340A1 - Gravity transducer and application to hydrocarbon exploration - Google Patents

Gravity transducer and application to hydrocarbon exploration Download PDF

Info

Publication number
WO2013177340A1
WO2013177340A1 PCT/US2013/042315 US2013042315W WO2013177340A1 WO 2013177340 A1 WO2013177340 A1 WO 2013177340A1 US 2013042315 W US2013042315 W US 2013042315W WO 2013177340 A1 WO2013177340 A1 WO 2013177340A1
Authority
WO
WIPO (PCT)
Prior art keywords
signal
method
gravity
recording
comprises
Prior art date
Application number
PCT/US2013/042315
Other languages
French (fr)
Inventor
George LISZICASZ
Xiang GUI
Original Assignee
Nxt Energy Solutions, Inc.
Forest, Carl A.
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 to US201261650413P priority Critical
Priority to US61/650,413 priority
Application filed by Nxt Energy Solutions, Inc., Forest, Carl A. filed Critical Nxt Energy Solutions, Inc.
Publication of WO2013177340A1 publication Critical patent/WO2013177340A1/en
Priority claimed from US14/550,779 external-priority patent/US9817153B2/en
Priority claimed from US15/727,157 external-priority patent/US10078152B2/en

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS
    • G01V7/00Measuring gravitational fields or waves; Gravimetric prospecting or detecting
    • G01V7/005Measuring gravitational fields or waves; Gravimetric prospecting or detecting using a resonating body or device, e.g. string

Abstract

A gravity transducer includes a particle system characterized by internal vibrations relating to its de Broglie wave, a resonant cavity for trapping the particle without holding it in a lattice structure; a source of a phonon wave, wherein the de Broglie wave and the phonon wave interact over a junction area; a power source for applying electrical power across the junction; a sensing system for measuring resistance, voltage, or current across the junction and for producing a sensed signal; and a recording system for recording the sensed signal. The transducer is used in a method of detecting potential hydrocarbon deposits, the method comprising: providing the transducer for sensing a change in a spatial orientation of gravity; flying the gravity transducer across the hydrocarbon deposit; sensing a change in spatial orientation of gravity to produce a signal indicative of geologic subsurface features, generally associated with hydrocarbon deposit; and recording the signal.

Description

GRAVITY TRANSDUCER AND APPLICATION TO

HYDROCARBON EXPLORATION

FIELD OF THE INVENTION

This invention relates in general to gravity transducers, and more

particularly to the detection of minute variation of the gravitational field induced by subsurface stress perturbations, and even more particularly to an application of such a transducer to hydrocarbon exploration. BACKGROUND OF THE INVENTION

Gravity variations have been detected using suspended large masses on the order of kilograms. It was believed the masses had to be large because the force of gravity is the weakest known force. More recently, atom interferometry, and in particular, accelerometers utilizing interferometry have been proposed as gravitational probes. See, for example, B. M. Anderson, J. M. Taylor and V. M. Galitski, "Interferometry with synthetic gauge fields", Physical Review A, vol. 83, 031602(R), 201 1 .

These devices have practical limitations. The devices using large masses can only detect large variations in the gravitational field. There are a number of complex technical challenges to overcome so as to make interferometers practical for field operations. Some of these requirements include cooling and keeping atoms near absolute zero and removing unwanted inertial effects. Thus, it would be highly desirable to have a device that could detect relatively small gravity fluctuations, such as those due to subsurface stresses relating to geologic features, and which could be used conveniently at reasonable temperatures, such as room temperature.

BRIEF DESCRIPTION OF THE DRAWINGS FIG. 1 is a diagram illustrating a typical geologic feature that the invention may be useful in distinguishing;

FIG. 2 is a block diagram illustrating the main components of a gravity transducer and detection system that is capable of identifying potential

hydrocarbon features; FIG. 3 is a detailed block diagram of the system of FIG. 2 showing the parts and cabling between the parts;

FIG. 4 is a sketch based on a photograph showing the system of FIGS. 2 and 3 in a Cessna Citation 560 aircraft;

FIG. 5 is a sketch of an open sensor case and an associated signal conditioning unit which are employed in the system of FIG. 4;

FIG. 6A shows side view of a stress field detector (SFD) sensor showing typical dimensions;

FIG. 6B shows an end plane view of the SFD sensor of FIG. 6A showing typical dimensions;

FIG. 7A is a side perspective view of another embodiment of an SFD sensor;

FIG. 7B is an end view of the SFD sensor of FIG. 7A;

FIG. 8 is a sketch showing a microscopic view of an SFD sensor illustrating how atoms may be trapped forming a sensor junction;

FIG. 9 is a sketch illustrating one embodiment of the invention comprising a multi-junction embodiment, each comprising a superposition of the wave of a trapped particle and the phonon wave generated in the junction cavity;

FIG. 10 is a sketch illustrating a standing wave of the Trapped Particle and a phonon wave front in a highly constricted region between a metal pin and a metal bead FIGS. 6A through 7B;

FIGS. 1 1 - 13 illustrate examples of signals generated by the system of FIGS. 2 through 6B on test flights that verify the abilities of the system, and more specifically:

FIG. 1 1 is a signal for a test flight line crossing Ladyfern gas field, BC,

Canada along a maximum stress line in the plains;

Fig. 12 shows a signal from a second test flight line Burnt Timber gas field,

AB, Canada using the same sensor as used in FIG. 1 1 , and along a maximum stress line, but in the foothills;

Fig. 13 shows signals from a third test flight crossing Ladyfern gas field, BC,

Canada, using the same sensor as used for FIGS. 1 1 and 12, but in this case the flight is along a minimum stress line; FIG. 14 shows a survey flight line map for a flight on April 14, 2012 over the La Jarilla oil field in Comodoro-Rivadiva, Argentina;

FIG. 15 shows a signal developed along the flight line of FIG. 14 taken off of sensor "Pearl";

FIG. 16 shows another signal developed along the flight line of FIG. 14 taken off of sensor "String", a different sensor from sensor "Pearl";

FIG. 17 is a sensor signal illustrating the concept of "signal amplitude";

FIG. 18 is a sensor signal in which the amplitude is averaged;

FIG. 19 is a sensor signal illustrating the concept of a geologic edge;

FIG. 20 is a sensor signal illustrating the concept of signal "frequency";

FIG. 21 is a sensor signal illustrating various frequency trends;

FIG. 22 is a sensor signal illustrating the concept of signal pattern change;

FIG. 23 is a sensor signal illustrating the concept of signal character change;

FIG. 24 is a sensor signal illustrating the concept of sequence identification;

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT

The present application discloses a quantum scale device for airborne geophysical exploration, which is sufficiently stable and for which the results are consistent. Thus, the new device and system have the ability to provide standalone reliability.

Embodiments of a practical gravitational field detector are disclosed herein. The specific application of the embodiment described herein is to the detection of potential hydrocarbon deposits. As will be seen in detail below, the geologic stresses in the region of hydrocarbon deposits in the earth's crust are different from in regions of solid rock. The detector described herein detects those differences. In the course of detecting the potential hydrocarbon deposits, it has been found that the embodiments disclosed also detect other geologic features within a sedimentary basin, such as traps and faults. It will also be evident to one skilled in the art that the stress field detector and the specific embodiments described has and will have many additional features and applications. Fig. 1 is a diagram illustrating a typical geologic feature 100 that the embodiments described herein may be useful in distinguishing. An exemplary area 1 10 of a sedimentary basin having a plurality of sedimentary layers 102, 104, 106, 108 and 1 10 is illustrated, which area includes areas having shale sequence 1 12 surrounding a porous sandstone formation 1 14, which may contain fluid, such as, for example, gas, oil, or water. In this example, the porous sandstone and the shale have the same density. Equal density implies equal gravitational

acceleration. Therefore, mass-based gravimeters yield the same measurement value over both. However, the two areas will have different stress regimes as suggested by 1 18, 126, 127, and 120. The shale rock 1 12, distributes stress differently from porous sandstone 1 14. The difference in the elastic (shear) properties between the shale and the porous sandstone containing significant fluid creates stress anomalies oblivious to conventional gravity detection devices, but detectible by the sensors described herein, due to changes in gravity field vibration associated with the subsurface stress perturbations.

I. DETAILED DESCRIPTION OF THE ELECTRONICS

FIG. 2 is a block diagram illustrating the main components of one embodiment of the gravity transducer and detection system capable of identifying potential hydrocarbon features. As shown in Fig. 2, the system integrates, stress field detection sensors, 154, signal conditioning systems 156, aircraft tracking and communication systems 160, data recording systems 164 and 166, quality control systems 170 including data post-processing systems. All the equipment is self- contained and is connected together through appropriate cabling 155, 161 , 158, 157, 167 and 168.

FIG. 3 shows a block circuit diagram showing the individual electronic components of an embodiment of the system 150 incorporated into a survey aircraft. Embodiment 150 includes a first Stress Field Sensor (SFD) system 202, a second SFD system 204, a main data acquisition computer system 214, a realtime display computer system 216, and a pilot navigation system 218. First SFD sensor system 202 includes an SFD sensor array including six SFD sensors and 2 accelerometers. These are contained in a sensor box 220. The structure of the sensors and the sensor boxes will be discussed below. First sensor system 202 also includes an SFD signal conditioning system 222 and a battery 225. Grounds throughout the system 150 are shown by ground symbols, such as 226. Second SFD system 204 includes a sensor box 230 containing eight SFD sensors and six accelerometers. Second system 204 also includes an SFD signal conditioning system 232 and a battery 235.

One embodiment of main computer 214 includes a QNX 4.25 data acquisition computer 290 together with cards and peripherals including keyboard 240 coupled to a keyboard IN port 299 via coupler 242. Computer 290 is powered from a power source 294, which is powered from the airplane, which power is input into 1 10 V AC power in unit 291 . Data acquisition computer 290 may also include data acquisition card 300 connected to the sensor data input port 250, data acquisition card 304 which is connected to accelerometer data input port 252, GPS serial port 306 which is connected to GPS Com 1 port 274, NMEA / qtalk serial port 308 which is connected to GPS Com 2 port 276, pilot display output 310, and USB port 314 which provides output to an external hard drive. Coupler 242 may be a cable or a wireless device and permits the keyboard to be moved about the airplane. Computer 214 may include: cable interface 248, M6259 2kHz data input interface 250, with up to 16 channels use by the SFD sensor; M6225 200 Hz data acquisition interface 252 for the accelerometers; a remote start output 256 for synchronizing the mainframe with the other computers; system one input port 258 connected to system one SFD signal conditioner 222; system two input port 260 connected to system two signal conditioner 232; and accelerometer/battery monitor input 264. Main computer 214 also includes a GPS module 268 which may include: GPS card 268; GPS antenna 270; Com 1 port 274; Com 2 port 276; DC power in port 280; GPS battery 282; and 15 volt dc source 286. GPS module 268 runs off battery 282, which is charged from the airplane generator when the airplane is flying.

Pilot navigation system 218 includes 5 volt DC power supply 320, pilot navigation interface 324 which receives the output from the pilot display output 310 of computer 290, and pilot navigation display 326 which is projected onto the windshield of the airplane.

Real-time display computer 216 includes a laptop 330, which may be a Panasonic CF30 with an OziExplorer moving map utility to track survey flights, and a WinDaq data acquisition system to show the SFD waveform during survey and also serve for backup data recording. Laptop 330 may also include

communication port 334 connected to the GPS communication ports 304 and 276 of computer 214, a USB port 326, RAM DC power in card 338 connected to 28 V DC power source 346, and USB port 340 connected to external hard drive 344. Real-time display computer 216 also may include mapping module 350 which may include data acquisition card 350, WinDaq input ports 356 and 358, and USB port 354 that connects to USB port 336 on the laptop 330. The WinDaq ports 356 connect to the SFD signal conditioning units 222 and 232, respectively.

Main computer 214 operates to collect and digitize sensor data,

accelerometer data, and GPS data. The data acquisition cards mentioned above include analog-digital convertors which digitize the analog outputs of SFD sensors and accelerometers. The accelerometer data serves a quality control function. If the impact of acceleration on the SDF sensors is too high, for example, 1 g or greater, extra care needs to be exercised with respect to data quality. This is performed by correlating accelerometer data with the SFD sensor data by the main computer 214. Main computer 214 also correlates the SFD data with the airplane position data output by the GPS module 268. Real time display computer 216 provides signal outputs during surveys that show the SFD output as a function of time, which waveforms will be discussed in detail below. These waveforms may be output in real time on the display of laptop 330. The time is also correlated to position via the GPS output enabling the system to provide real time display of flights which correlate time with position. Examples of such maps will be shown and discussed below.

Fig. 4 is a sketch of a photograph showing the sensor survey equipment in place aboard a Cessna Citation 560 aircraft 380. Real time display laptop 330 sits at the top of the nearest rack 384, while main computer 290 is directly below it. Mapping module 350 is at the bottom of this rack. Cable interface 248 can be seen at the bottom right of the rack 384. System one SFD sensor case 230 can be seen at the top of second rack 388 while system two sensor cases 230 is below it. System two signal conditioners 232 can be seen on the rack level below the sensor case 230. Navigation module 218 can be seen in the cockpit 390 at the far end of the airplane 380. Fig. 5 shows one embodiment of an open sensor case 230 and a multiple- channel signal conditioning unit 232. The gravity field sensors are individual units mounted inside a pelican hard case 400 with vibration and shock isolation. In the embodiment of FIG. 5 the vibration and shock isolation comprises a silicone gel mat 410. Sixteen SFD sensors 404, 406, 408, 414, 416, 418, 420, 422, 424, 426, 428, 430, 434, 435, 436 and 437 are mounted on the mat 410. Three

accelerometers, 440, 441 , and 442 are mounted on top of SFD sensors 434, 435, and 437, respectively, and thus are supported on the mat 410. These

accelerometers measure acceleration in each of three orthogonal directions. In addition, three accelerometers 444, 445 and 446 are mounted directly to case 400. These measure acceleration in each of three airframe orthogonal directions. Many other configurations of SFD sensors and accelerometers are possible.

Redundancy of SFD sensors and accelerometers is often a feature of

embodiments. In one embodiment, signals from fourteen SFD sensors in two separate systems are recorded simultaneously in differential mode to reduce uncertainty of identifying subsurface geologic features and eliminate unwanted noise. In another embodiment, twenty-four SFD sensors are contained in three separate systems.

The gravity field sensor according to the invention are quantum devices that respond to structural, stratigraphic trapped fluids caused by variations primarily in the horizontal stress directions in the subsurface, hence the term Stress Field Detectors. The preferred embodiment of the sensor is described in reference to FIGS. 6A - 10 below. Depending on the design, the sensor is a two or three terminal transducer that produces a time-varying analog electrical signal, which is fed through a multiple channel signal conditioning unit. The output, in the form of an analog voltage signal, flows into the data acquisition cards where it is converted to a digital signal for display on a computer system. Data may be recorded at 2,000 samples per second.

FIG. 6 shows side (FIG. 6A) and plan (FIG. 6B) views showing details of one embodiment of a sensor junction 416 according to the invention. The junction comprises a pin 450, having a coating 462 and a bead 454 comprising a first half- sphere 456, a disc 458 and second half-sphere 460. Pin 450 may be made of gold or copper, though it may be made of other materials. Coating 462 may be made of silver oxide, i.e., Ag20, though it may be made of other oxides or other dielectrics. The thickness of the coating 462 may be thirty nanometers or less; more preferably, twenty nanometers or less; and most preferably ten nanometers or less. The circular center disc portion 458 of bead 454 may be made of copper, though it may be made of many other metals. Exemplary dimensions of the various parts are shown in FIGS. 6A and 6B. Half-spheres 456 and 460 may be 2mm long, and may be made of alloys, such as alloys of tin, nickel, lead and silver. In one embodiment, first half-sphere 456 may be 63%Sn and 37%Pb, while second half-sphere 460 may be 96.3%Sn and 3.7%Ag. Disc 458 may be 1 .75 mm wide. Pin 450 may be 25 mm long and 1 .45 mm in diameter. The circular diameter of half-spheres 456 and 460 and disc 458 may be 3.5 mm. In operation, pin 450 is held at the ground voltage via contact 504, while bead 454 is held at the positive VDC voltage via contact 508. Some other materials for the sensor components detailed above may also be used such as nickel, carbon, platinum, others and their alloys. Metal oxides that may be used to coat pin 450 include tantalum oxide, aluminum oxide, tin oxide, nickel oxide and others.

FIGS. 7A and 7B are side and plane views showing perspective of another embodiment of a sensor junction 470. This embodiment includes pin 474, and bead 476 disc 478, disc 480, and disc 482. Pin 474 is coated with coating 472. Pin 474 and disc 480 may be made of a good conductor, such as gold or copper.

Discs 478 and 482 may be made of a metal alloy, such as alloys of tin, nickel, lead and silver. Coating 472 may be made of a metal oxide, such as silver oxide or other oxides mentioned above. A pin contact 490 may be made of copper, or other materials mentioned above, and may be grounded at 494. Disc 482 may be connected to the positive VDC voltage. The positive VDC voltage may be typically between 6 and 9 volts. In operation, the position of the pin 450, 474 is adjusted in the bore of the bead 454, 476 using a precision micro-mechanical control, which may be adjusted via a piezo-driven nanopositioner.

II. OPERATIONAL PRINCIPALS

While the system described herein is largely tested against empirical evidence over approximately 15 years, and steady improvement of the system has been made based on the empirical evidence, a discussion in terms of operational principals may be useful in understanding the system. FIG. 8 is a sketch showing a microscopic view of the sensor junction 416 illustrating how electrically neutral or charged particle(s) 512, 514, may be trapped forming an active sensor junction 416. As in FIGS. 6A and 6B, sensor 416 includes a bead 454 and a metal pin 450. Bead 454 includes portions 456 and 460 which may be made of metal alloy and portion 458 which may be made of metal, such as copper. Pin 450 is coated with a coating 462 which may be made of a metal oxide. Pin 450 is grounded at 504 and bead 454 is connected to the +VDC voltage at 508. Pin 450 is adjusted in directions orthogonal to its length, as shown at 520. As pin 450 is adjusted, electrically neutral or charged particle(s), such as 512 and 514 may become suitably trapped between pin 450 and bead 454. When the +VDC voltage is applied across the pin450 and bead 454, electrical current produces heat, which generates phonons. A particle, such as 514, can become suitably trapped in a phonon resonance forming at electronic junction of dissimilar metals.

FIG. 9 is a sketch illustrating one embodiment comprising a multi-junction embodiment, each comprising a superposition of the wave, such as 515, of a trapped particle, such as 514, and the phonon wave, such as 530, generated in the junction cavity. As shown, adjustment 526 traps particle 514 in a standing phonon wave 530, and adjustment 536 traps particle 533 in a phonon standing wave 534. The trapped particle is confined in the gap 500 between oxide coating of pin 450 and the inner surface 460 of bead 454. FIG. 10 is a sketch illustrating a trap 520, 501 of a particle 514 between metal pin 450 and a metal bead 460.

The phonon-wave controlled, active interface, resonance transducer is a geophysical sensor system which is capable of detecting subtle variations in the earth's gravitational field. Due to its compact and highly portable configuration, the sensor is easily employed in airborne survey applications, which ideally position it as a rapid and cost-effective tool for wide-area mineral exploration. The present survey system is designed specifically for oil and gas prospecting. Additional future applications of the system include metallic mineral exploration, kimberlite (diamond) prospecting, locating shallow fresh water beds, the prediction of earthquakes and volcanic eruptions, identifying and testing quantum mechanical interactions, (phenomena and processes), investigating quantum mechanical events in the domain of general relativity, exploring the understanding of the origin of mass, gravity and inertia from the most fundamental laws of physics, and testing general relativity postulates, and potential astronomical and astrophysical observations. However, these applications may require modifications to the current designs and material selection, added control functions and the selection of suitable particles.

The basic operating premise of the active interface is based on the fundamental relationship between stress and gravity.

The vertical stress component of the earth is dominated by gravity acting on the overburden. In contrast, the maximum and minimum horizontal stress components significantly affect horizontal gravity. Accordingly, the majority of gravitational anomalies with orientation changes will occur in the horizontal direction. Hayes showed that stress fields not only will introduce changes in the magnitude of gravity but also couple in the tensor directions (T. J. Hayes, K. F. Tiampo, J. Fernandez, and J. B. Rundle, "A gravity gradient method for

characterizing the post-seismic deformation field for a finite fault" Geophysical

Journal International, vol. 173, pp. 802-805, 2008). By using extensive calculations and models of Joshua Tree Earthquake, Hayes showed that gravity gradiometry will exhibit similar shape as stress model over a finite strike slip motion in horizontal plane. One exception where stress does not follow the gravity gradient is in trapped fluids. If a solid is in contact with a fluid, the shear stress at the interface must be zero (fluids cannot support shear) but the normal component of the stress (i.e., the traction) must be continuous. The effects of seismic shear waves are identical for trapped fluids. Stress fields of trapped fluids, therefore, differ significantly from those of solid rocks and also cause unique gravitational perturbations.

According to general theory of relativity, all forms of energy are sources for the gravitational field, and the addition or subtraction of energy alters the gravitational field potential of an object. Consequently, gravitating bodies will develop real perturbations in the geometry of their gravitational field in response to applied stress. Stressors may include heat, rotation, compression or tension, and the origin of such stressors may be internal or external to the body.

The transducer of the invention is a mesoscopic scale non-linear device designed (1 ) to selectively respond to stress induced gravitational anomalies; (2) to achieve a dynamical signal integration over time; (3) to reject constant or gradually varying gravitational sources; (4) to perform these functions from high altitudes . e.g. 1000 to 20,000 feet, and at a high sample rate (e.g., 2,000 samples per second); (5) to allow the interaction of classical and quantum mechanical domains through a specially designed interface. The gravity transducer system exploits the stress-gravity relationship to infer changes in stress patterns within a sedimentary basin and/or in the earth's shallow crust by focusing on the detection of subtle variations in the form of gravity field vibrations coming from stress perturbations. To be clear, the gravity transducer system according to the invention is designed to detect anomalies, as variations in vibrations, in the gravity-energy frequency domain, as opposed to measuring the magnitude of acceleration which is what the state-of-the-art conventional gravimeters and gravity gradiometers rely on. As such, the system of the invention represents a substantial departure from known technologies employed in the field geophysics.

A large 'test mass' is unsuitable to measure small distortions in gravity on large scale interactions. Furthermore, large mass has short de Broglie

wavelength. In contrast, particle scale mass has long de Broglie wavelength. Energy density is related to mass, and mass could be defined in terms of de Broglie frequency (J.W. G. Wignall, "Proposal for an absolute, atomic definition of mass", Physical Review Letters, vol. 68, pp. 5-8, 1992; S.-Y. Lan, P.-C. Kuan, B. Estey, D. English, J. Brown, M. Hohensee, and H. Muller, "A clock directly linking time to a particle's mass", Science, vol. 339, pp. 554-557, 2013). Therefore, all vibrations experienced by a mass are fundamentally associated with its de Broglie frequency.

As the sensor is flown though the gravitational field at a high speed, its vibrating active interface is tuned to a specific frequency of interest. This will allow the sensor to interact with the localized gravitational energy density within a specific geometry of spacetime relating to a specific subsurface stress

perturbation. The perturbation with an energy density comparable to that of the sensor will bring about localized vibrations interfering with the pre-set vibrations of the sensor and initiate changes in the phonon population within the resonating junction cavity. Thus the resulting electrical waveform output can be used to identify signal attributes indicative of trapped fluids and other geologic features. In addition, orientation effects are induced by subsurface in-situ anisotropic "stress states" associated with subsurface discontinuities, which can be detected as the sensor is flown along straight path. Stress perturbations associated with buoyancy (an upward force exerted by a fluid that opposes the weight of an immersed object) can be identified as well, based on stress energy patterns and frequency shifts.

Long wavelength gravity field vibrations, such as those produced by stress perturbations, can be transferred to particles with the corresponding wavelength of, for instance, the detection element consisting of the phonons and the trapped particle. The process described above allows an efficient energy-momentum transfer and the amplification of the resulting effects of small field vibrations and distortions. SFD probes gravitational energy density spectrum in order to investigate and discern stress changes in the subsurface associated with trapped fluids and other geological features. The present invention exploits a novel approach in identifying and describing gravitational effects that were undetectable before.

The device describe herein may be a phonon standing wave controlled transducer that consists of a resonant cavity or a low resistance, weakly interacting metal-oxide-metal junction (LONCI Junction) with highly constricted spatial regions. Between the junction terminals a particular non-lattice bound matter structure having its associated wave, i.e., a particle element may be trapped. The particle element will be called the Trapped Particle hereafter. With no excitation applied, the Trapped Particle would naturally oscillate about its equilibrium position.

However, with proper excitation, the Trapped Particle may behave as an Active Interface in conjunction with the phonon standing wave in the detection of localized gravitational vibrations or distortions. The coupling of the Trapped Particle matter- wave and the phonon standing wave constitutes the gravity detection mechanism.

Other quantum mechanical devices employ "passive junction interfaces" in which the junction barrier is composed of gas/air or thin film or some other dielectrics, etc., characterized by high resistance and usually operates at extremely low temperatures. In contrast, an active interface sensor may utilize low junction resistance in which particles may move through the junction as waves, as opposed to high resistance junctions of about 4.1 ΚΩ or more where particles moves through the junction as particles. Further, the described sensor operates at room temperature and employs Active Interface in a form of a vibrating Trapped Particle matter-wave, which is coupled directly with the phonon standing wave system. The resulting arrangement is used in the detection of highly localized minute vibrational/directional changes in the gravity field.

The LONCI junction is biased by DC current. The passage of the current through the low resistance junction produces Joule heat. When phonons are absorbed in the junction they are rapidly thermalized which raises junction temperature. As a result, the sensor system undergoes excitation and phonons are generated. This results in electron-phonon scattering, junction impedance change, Coulomb blockade formation, and initiates particle tunneling. The system is tuned electromechanically (1 ) to attain the desired resonance by modifying the cavity dimensions and/or changing junction capacitance, (2) to produce coherent standing phonon wave within the cavity out of two travelling phonon waves in opposite directions, (3) to allow a particle wave complex to couple to the gravity field locally, and (4) to amplify signals, and (5) to exploit stochastic resonance. By the assistance of ubiquitous Johnson-Nyquist electronic noise (due to thermal agitation of charge carriers) generated in the sensor junction, amplification and optimization of extremely weak signals resulting from gravity interactions may be achieved. It also affects sensor mode operations.

As the sensor is moved along a survey path, the matter-wave of the

Trapped Particle coupled with the phonon standing wave will stretch, compress or relax in response to field variations. The changing wave system continually affects junction impedance, which is observed as a voltage signal output against time. This is partly due to electron-phonon scattering with controlling electric charge distribution across the junction.

While the impedance change is evident, its cause is not yet fully

understood. Some of the processes governing the impedance change are potentially: superposition of coupling of fields and waveforms, phonon-phonon coupling, electron-phonon coupling, Coulomb Blockade formation, and particle tunneling. As bias voltage is applied, charge carriers can tunnel through the junction. Conductance varies according to the charging energy of junction, thermal energy of the system, conductance quantization and spin polarization, as charge carriers move across the oxide layer. The current transducer system is calibrated to focus on an effective observation window of approximately 0.5 to 1 .5 radial kilometers. At this scale, empirical evidence suggests that the majority of gravity effects arise from the interplay of in-situ tectonic stresses acting upon local geologic features. As a result, the inferred stress patterns are often indicative of prospective hydrocarbon trap and reservoir conditions in sedimentary basins including, for example, the presence of fault blocks, channelized reservoirs, or carbonate reefs. Certain other stress patterns may indicate geologic domain changes, such as fault systems, fracture belts, lithologic boundaries, and over-pressure.

IV. EXAMPLES OF HYDROCARBON SENSING RESULTS

To verify that the system according to the invention works, flights were made over several distinct potential hydrocarbon locations and the results compared to previous knowledge of the fields. The first of these fields is the Canadian onshore carbonate field referred to as the Ladyfern Gas Field, in BC, Canada. Found in 2000, the Ladyfern field was estimated to contain 1 .5 Tcf of natural gas reserves. Three flight lines were flown in different directions to confirm the relation of porous areas to the sensor signal. The Ladyfern gas field is enclosed within a carbonate barrier reef with porosity preserved through the process of leaching and dolomitization. Faults and fracture systems on the boundaries are strike-slip in nature. Based only on seismic, the dolomitization and trapping is not readily visible, although the sag features are quite clear.

Ladyfern is a natural gas accumulation hosted in a stratigraphic trap - Devonian slave point reef. The productive reservoir fades have developed via burial digenesis called hydrothermal dolomitization ("HTD"). Splays and antithetic faults associated with the nearby Hay River dextral transform system are believed to have acted as conduits for the hydrothermal fluids. The aggressive nature of HTD digenesis has created a series of 'sag' or 'collapse' features in and around the reef, which are readily visible via seismic studies.

Figs. 1 1 - 13 show the sensor signals obtained from two flights flown over the Ladyfern field and one flight flown over the Burnt Timber gas field, AB, Canada. Each of FIGS. 1 1 through 13 shows signal output in volts versus flight time. The output voltage is normalized to fall between zero and 1 volt. The resultant sensor signals indicate sub-surface geological changes and anomalies which are indicators of potential hydrocarbon bearing reservoirs. A very strong correlation to the distribution of the trapped reservoir was found. In this field, it was difficult to identify the carbonate pockets using only seismic, and the sensing system described herein provided a good corroboration of reservoir opportunities.

FIG. 1 1 shows a sensor signal 600 developed along a maximum stress line over the British Columbia plains. At 602 the signal begins to rise showing an approaching geological change. Sharp peaks at 604 indicate stress fractures and a geologic change indicated at arrow 606. The curve can best be understood by viewing it as a whole. It is seen that there is a strong anomaly at 614 bracketed by similar signals 608 and 620 on either side. Following the interpretation guidelines discussed below, we interpret this as a hydrocarbon reserve in the area 615 with typical approaches at 610 and 619. At 618 the signal shows the rise out of the anomaly. The area 615 corresponds to an area of drilled wells. We also can see a minor anomaly at 625 indicated by the bar 626. This is an undrilled area. At 630, we see the decline of the signal indicating the plane is leaving the anomalous area.

FIG. 12 is a signal 630 developed with the same sensor as used in FIG. 1 1 , but in the foothills of the Rocky Mountains and also along the maximum stress line. Here a geological change is indicated by the sharp peak at 632, which is marked by the arrow 633. The signal begins to rise at 634 along the direction 636. The bracketing structures are at 642 and 654 in the areas 640 and 658, respectfully. An anomaly 650 is indicated by the dip in the signal at 646 in the area 644. The signal falls off at 664 indicating a geological change indicated by arrow 660. Overall, this signal is much more fractured, such as shown at 638 and 646 which indicate that the terrain is more fractured.

FIG. 13 is a signal 670 developed over the same anomaly as the anomaly 614 of FIG. 1 1 , except along the line of minimum stress. The same anomaly shows up at 672 in the area 674.

FIGS. 14, 15 and 16 show signals developed in an R&D test flight over the "La Jarilla" oil field in Comodoro-Rivadiva, Argentina. This particular series is shown to illustrate that different sensors will not show identical signals. FIG. 14 shows the display 700 as shown by the mapping computer 216. Display 700 shows the "La Jarilla" field at 707 and several smaller fields at 720 and 722. The flight identification number is shown at 712 and the flight direction is shown at 714. The area surveyed is shown at 710, and the time line is shown at 714. FIGS. 15 and 16 show the same area and time covered, but with different sensors.

Comparing FIGS. 15 and 16, it is seen that a local geological change 734 in FIG. 15 and 750 in FIG. 16 and the start of an anomalous area show up differently. In FIG. 15 it is indicted by a sharp peak at 730, and in FIG. 16, it is indicated by a signal buildup at 752. A local geological change 736 in FIG. 15 and 758 in FIG. 16, show up similarly as minor peaks 738 in FIG. 15 and 755 in FIG. 16. The anomaly is shown by a dip at 741 and 742 in area 766 FIG. 15 and a single dip at 762 at area 740 in FIG. 16. The end of the anomalous area is shown at 746 in FIG. 15 by a frequency change and amplitude shift. This is also shown by a frequency change and amplitude shift at 770 in FIG. 16, though the change is not as pronounced. These curves show that what is important is the changes in the signals, not the absolute values of the signals. This will be made clearer in the signal interpretation section below.

V. OPERATION OF THE SYSTEM

1 . Aircraft Tracking and Communication Systems

The aircraft tracking and communication systems consist of two parts. The first part is a "flight-path following system" that provides the capability for the recording operator in the aircraft and the ground crew to follow the flight path and ensure the aircraft is on its intended course. The system updates the position of the aircraft each minute. This ensures that the ground crew has access to the plane's location with a high degree of accuracy. The main purpose of the tracking system is to conduct safe flights and to allow the R&D engineers, and, in the future, potential clients, government and military agencies to track the aircraft accurately during the survey flights. The system is satellite-based through the iridium network which enables global access to position data and enables voice and data communications outside of cellular coverage areas.

The second tracking system is connected to the GPS navigation and gives a real-time- map based update to the sensor operator who is monitoring the onboard quality control system. The onboard Novatel GPS System is supplemented by a subscription service to Omnistar for wide area differential GPS Corrections. This setup provides a very high degree of data positioning accuracy within +/- 2 meters.

2. Data Recording Systems

Two computer systems are used to record the data independently and simultaneously: the QNX system with 16-bit precision and the WinDaq system with 14-bit precision. All systems are controlled via the main QNX computer. Specifically, the QNX computer maintains automatic synchronization of the redundant data acquisition by automatically signaling the start and stop of recording within 0.5 seconds. Furthermore, the QNX computer also records the GPS position and elevation information, visual ground references and any other relevant in-flight data required for proper quality control at the end of the flight. GPS elevation data is utilized along with USGS-provided DEM data at the end of the flight to calculate the height of the plane above ground. This ensures that the flight, preferably done at 10,000 feet, was conducted within the altitude tolerance limits of the sensors. However, flights can be conducted within a range of 1 ,000 feet and 20,000 feet.

3. Data Quality Control Systems

The quality control system incorporates the raw digital output from QNX and the graphical display of the output from the WinDaq acquisition system. Since this is similar to the signal that is used for interpretation, determination of the quality of the signal can be made in-flight. The result from the real-time GPS output and the planned fight are compared to ensure a minimal deviation. The cross track, an indicator of how far off the planned course the aircraft is located, is monitored in real time during the flight both by the operator and the pilots. Normally, a cross track distance within 10 meters is maintained.

In addition, the quality control system contains power supply monitors to ensure the maintenance of proper voltage levels throughout the data acquisition process. Acce I ero meters are also used to monitor the acceleration of the aircraft on the vertical and two horizontal axes during the data acquisition process.

4. Data Post-Processing Systems

At the end of the flight, the QNX system automatically transfers all data onto an external disk. The raw data is checked to ensure the correct record lengths. In particular, the data is checked to ensure that both QNX and WinDaq acquisition systems have recorded the raw data with no problems.

The raw data in hexadecimal format is processed to check for consistency of data and is then reformatted for input into the DADiSP engineering worksheet. DADiSP is a commercial interactive graphics worksheet, a visually oriented software package for the display, management, analysis and presentation of scientific and technical data. It can be used to manipulate, edit, reduce, transform, display and analyze technical data. DADiSP can handle and simplify the large volumes of data output by the sensor system. As mentioned above, digitized sensor output contains 2,000 points per second and flights can average about

10,000 seconds, resulting in a 1 GB per flight data volume or more in a raw form.

Once the data is loaded into DADiSP, it may average more than 6 GB per flight.

This large volume of data from two recording systems can be handled by DADiSP in an efficient manner.

5. Flight-Planning Process

The overall sensor survey flight planning process starts by researching and understanding the following:

1 ) Maps - topographic and aviation maps of the target area are prepared in both paper and digital forms.

2) Sensor technology requirements - certain parameters must be followed to obtain reliable, quality data. These include long straight segments, large turning radius and low turbulence. As a rule of thumb, in one embodiment, the minimum straight length is 100 km, the minimum turning radius is 25 km, and the acceleration is not greater than 1 g.

3) Base of operations for the aircraft - selected based on the following

criteria: proximity to the survey area, runway length (minimum 5000 feet), fuel availability (jet A1 with Prist), flight planning facilities, security, accommodations and communications facilities.

4) Other flight restrictions - these include both regulatory agency and

geographic considerations as follows:

• Air traffic control rules that include departure and entry vectors, permit approval processes and IFR flight rules.

• Air Force permits required for restricted and prohibited air space. • Military restrictions that include hostile regions, border buffer zones, and altitude restrictions.

• Geographic restrictions including mountains and borders.

• Prevailing weather trends including precipitation, winds, storms, and clouds.

The above information is utilized to begin designing each individual target flight line. Flight lines are designed as a sequence of straight segment waypoints and accurately controlled turn waypoints.

In one embodiment, to get ready to fly, there are six processes that are performed:

1 ) The waypoints are loaded into the aircraft's flight management system.

This system guides the aircrafts autopilot mode over the planned route to a high degree of accuracy.

2) The waypoints also are loaded into the sensor flight tacking system. This enables the onboard sensor operator to track the progress of the aircraft to ensure that the correct flight path is followed. This system also provides accurate location information for the data that is being collected.

3) Sensors are loaded onto the aircraft and calibrated.

4) Pilots file the flight plan with the relevant administration of air traffic control. 5) Pilots and the sensor operator board the aircraft and perform the planned flight.

6) After completing the flight, the data then is evaluated to ensure that it meets the prescribed quality standards.

In this process, it is important to realize that gradual changes in altitude (in the vertical direction) do not affect the sensors adversely; however, a slight aircraft rotation along the horizontal plane does. The noise inherent within the sensor signal largely originates from rotational aircraft movements.

6. Data Interpretation Process

The sensor signal interpretation process is generally performed at the head office, and in one embodiment, involves three major steps. The first step is the single line sensor by sensor interpretation. This is followed by the multiline sensor signal integration. Then the results go through a final ranking and mapping. The identification and characterization of sensor signals is the end goal of the sensor interpretation process. The process of sensor signal interpretation involves the recognition the anomalous patterns within the waveform of a sensor signal from its background. The key challenge is to recognize a section of anomalous changes in signal character from a noisy background signal with confidence. Sensor signal interpretation employs several methods to improve the ability to identify an anomaly including spatial confirmation of anomalies across multiple sensors, recognition of the specific sequence of signal character changes, and independent identification of anomalies by multiple teams. For interpretation purposes the signal is continuously recorded from the beginning of a stable flight path until just before landing. In one embodiment, it is important to the

interpretation process to capture the relative changes in the "baseline" signal characteristics for each sensor over the entire flight. This enables the recognition of the normal background waveforms and the specific pattern or mode changes.

For a potential anomaly to be subject to interpretation it should show up on a plurality of sensors. The use of multiple sensors plays a significant role in quality control. For a potential anomaly to become further analyzed, it has to be independently located in the same spatial location by at least two separate teams of interpreters. Additionally, in all instances, it needs to occur in a specific sequence.

A key part of the interpretation process is the recognition of the trends and character changes of the background signal and the identification of the sensor "mode". A mode is a long term signal trend which is recognizable for a particular sensor.

A sensor signal character is considered anomalous over an area when it conforms to an idealized pattern or sequence of responses that are known to correlate to geologic changes in the subsurface. These changes include the development of fluid charged reservoir or the existence of structural or

stratigraphic trapping mechanisms. Although the sensor signal can indicate the presence of trapped fluid, it currently cannot distinguish between hydrocarbons and water.

Anomalies can exist as point changes in the signal patterns and are often indicative of a significant change in geology. Area anomalies consist of signal changes with a defined beginning and end and will have a recognizable shape, sequence, character change or pattern. Anomalies in the gravity sensor signal are often recognized by frequency changes over different rolling window sizes, trend changes over different window lengths, and magnitude and slope changes over an area.

The correct recognition of "changes" is an intricate process. Fundamentally, the process is complicated by the existence of several background modes of each sensor. The anomaly character development will depend on the background mode. Additionally, the character in the signal is relative; that is, the initial conditions will affect the subsequent sensor responses. For example, crossing a large fault can mask a smaller anomaly that closely follows it.

The output of the data sensor interpretation process is a prospectivity map, which shows recommendation boxes that describe the spatial extent of the reservoir portion of the geological anomaly within the area of the sensor survey.

The following section covering FIGS. 17 through 24 is a basic description of the most important concepts involved in SFD data interpretation using ideal or near- ideal sensor response examples. Each of FIGS. 17 through 24 shows the signal voltage on the ordinate as a function of time along the abscissa. This description is intended to introduce one skilled in the art to the sensor data interpretation process. The most common anomaly attributes and sequences are presented. All major anomalies may be identified by one of the following generalized

interpretation concepts, and many less significant anomalies may be identified by either an extension, corollary or grouping of these concepts. Because sensors may exhibit stochastic behaviour, this description may not apply to all sensors. This description focuses on the concepts applicable to a few of the primary sensors. The concepts can be extended with some variations to other sensors. However, it should be understood that it may be advisable to learn the response of each sensor individually by empiric data. Many subtle variations of these examples exist which may or may not alter the anomaly identification. To serve the above purpose, the emphasis of this document will be placed on the detailed visual inspection of the waveforms, although various signal processing techniques may be applied in the interpretation process.

A. Signal Attributes Used in SFD Interpretation The following signal attributes are utilized in this disclosure for a majority of the interpretation process: Signal amplitude; Signal frequency; Signal character and pattern. The current representation of the SFD sensor output is a simple 2D voltage-versus-time plot, which therefore implies that each of the above mentioned attributes can be represented by a particular behavior of the sensor voltage output as a function of time. This section of the document therefore attempts to describe each of these as applicable to the SFD interpretation process. It is important to note that SFD attributes are not defined in absolutes. Each of the following terms/attributes refers to a sequence observed over a floating time period, nominally of the order of 10-50 seconds.

I) Amplitude

The term "amplitude" as applied to SFD implies a change in the observed mean output voltage of the SFD sensor during a certain interval. There is no fixed definition of the time interval being used to define the mean output voltage nor is it required that the time intervals be of equal duration. FIG. 17 highlights the output voltages and their perceived mean values in a series of time intervals. In sections A and B the mean values shown by arrows 802 and 803 do not appear to have changed; however the two areas may be distinguished on the basis of the remaining signal attributes: frequency and character. In Section C the amplitude shown by arrow 806 has changed. This section may be classified as a transition period from A and B to D. The term amplitude, and quantities and concepts derived from it such as amplitude change, voltage baseline etc., would therefore refer to the voltage value represented by the arrows 802, 804, 806 and 808 which are 0.7V in sections A and B, 0.65V in area C and 0.55V in area D, as opposed to the discrete time amplitude value such as 0.69V at timestamp 6475. Amplitude and Amplitude trends are therefore directly realized by the mean magnitude of the signal voltage in any given time window and play an important role in the identification and ranking of the anomaly. With reference to the previous example shown in FIG. 17, calculating a moving average using 5 second time intervals over the same signal duration yields the curve 810 of FIG. 18. While such a moving average is useful in identifying the major trend intervals, in D and E it tends to obscure the trend of the baseline voltages within the time window. These trends often play a role in determining the ranking of the SFD anomalies. Thus, in this embodiment, moving averages are not the best way to analyze the data.

FIG. 19 shows an example of an edge 820 on a SFD signal 812. The edge is an example of a sudden amplitude change between a first area 814 and a second area 818 and it may or may not be accompanied by a mode/background pattern change. An edge is often used to interpret the transition of the sensor from one setting to another - identification of these points is generally used to determine the start/end of signal sequences that are used to define anomalies.

II) Frequency

The term "frequency" as applied to SFD implies the number of observed significant peak-to-peak (Vmax-Vmin) values during a certain time interval. The significance is determined based on the comparison with the waveforms prior to and after the time interval of interest. There is no fixed definition of the time interval being used to define the mean output voltage nor is it required that the time intervals be of equal duration. In most cases an identified frequency change or shift will need to be a length longer than 10 seconds (based on the current flying speed of the aircraft) to be considered applicable for SFD interpretation purposes. Consider signal 830 in FIG. 20. The areas A and F will fit the definition of the high frequency characteristics described previously while areas C and D will be considered low frequency areas. Area G can be considered as an example of a tight mode signal and in this case the notion of frequency as applied to SFD is not relevant. Anomaly identification in this mode will be driven primarily by amplitude and pattern characteristics. In particular note area E in which four sub areas a, b, c and d are defined. Within these sub areas, the Vmax-Vmin differences and the number of times the particular change described by the Vmax-Vmin event which happens during the transitions from a-b, b-c and c-d will fit the high frequency definition and more appropriately can be seen as a frequency transition from A-B- C-D to E. Hence, only area E will be considered to have high frequency relevant for anomaly identification purposes.

Signal 840 of FIG. 21 shows a transition from a low frequency area (A) to a high frequency area (B) and finally to a tight-mode area (C). It is quite apparent that area B is the densest part of the signal. The event described by the change in the mean Vmax-Vmin value can therefore be simplified by noting the following two items: it represents a change from a low Vmax-Vmin to a high Vmax-Vmin along the transition from A to B; and, the frequency of the high Vmax-Vmin event is much higher in area B. It may be noted that the reverse event (high to low Vmax-Vmin) transition is equally applicable to interpretation and anomaly identification. This is a direct consequence of the original observation that singular time events and changes on SFD are generally meaningless and only hold value if part of an overall signal development. Frequency and Frequency trends are therefore directly recognized by the occurrence and regularity (or lack thereof) of a significant standard deviation change of the signal voltage in any given time window and play an important role in the identification and ranking of the anomaly.

Ill) Character & Pattern

A character change of an SFD signal can imply the occurrence of large and diverse number of events. Character changes can be very loosely tied to mean output voltage changes but this is not always the case. In general a pattern change on a SFD sensor's output is observed as a shift or change in a sequence of regularly occurring signal events. These events are a combination of both amplitude and frequency shifts as described above. Most of these discernible pattern changes are derived from historical surveys and template crossings._A character change in general determines the overall shape of the signal as characterized by amplitude, frequency, envelope evolution etc. A particular waveform consisting of these character changes is visually noticeable on a SFD signal. Some SFD signal events associated with character changes include: A change or shift in the amplitude trend of the signal; A change or shift in the character of the signal development - for example a signal showing only positive deviations about the perceived baseline moves to a state where it shows both positive and negative deviations; A change or shift in the signal pattern - for example a signal showing a high frequency baseline trend moves to a low frequency baseline trend or a varying frequency trend; A change or shift in the signal envelope - for example a signal exhibiting a cyclic/periodic behaviour to a more discrete development. Signal character changes are therefore directly realized by a sequence driven combination of the frequency and/or amplitude characteristics as described previously in any given time window and play an important role in the determination and ranking of the anomaly. Some signal pattern changes are illustrated in signal 850 in FIG. 22. At the arrow 856, it is seen that the baseline voltages move up from area 854 to area 860. Also at this point the response type changes to high frequency with a U-shaped waveform development, specifically U's 861 and 862, which are accompanied by high frequency such as at 863. Further, when comparing the signals on either side of arrow 856, it is seen that the response range changes from mostly above the baseline to either side of the baselines.

FIG. 23 illustrates several signal character changes that are relevant to signal interpretation. The baseline voltage of the signal at 867 is shown at 868. The baseline voltage of the signal at 869 is shown at 870. It is evident that there is a character change in that the baseline voltages move down. It can be seen that the downward development in the area 867 changes to high frequency with only upwards development in the region 869 which forms a dip in the area 866. It is also evident in FIG. 23, the response changes from mostly below the baseline at the start to either side of the baseline in the region 869 and then back to the high frequency downward trending signal in the region 872. These are all examples of signal character changes.

B. Sequence Identification

The initial point of the SFD interpretation process is the determination of the start and end of the applicable time windows i.e. - the timestamps at which the change in the frequency, amplitude or character/pattern can be observed &/or established. The identification of these time windows may be useful in establishing the overall signal quality sequences that will eventually be used for anomaly identification. The determination of the anomaly and its subsequent strength will be determined by the type and sequence of attribute changes observed within & at the immediate periphery of these time windows. Consider the signal sequence shown in FIG. 24. It has been divided into a succession of time windows - A, B and C which have been differentiated based on a change or shift in one of the primary SFD signal attributes. The transition from area A to area B is defined by the following events. As shown by comparing the levels 876 and 872 there is a baseline amplitude change and a signal character change in that the signal changes both trend, i.e., upward trend to flat trend, and pattern change, i.e., the (directionality of the deviations about the baseline. There is also a signal frequency change as comparing signal 874 to signal 880. Similarly the shift from area B to area C is defined by the following events: a Baseline amplitude change as shown by arrows 872 and 886; and a signal frequency change as by comparing the signals at 880 and 884. The determination of these sequences is useful in the interpretation of the signal with respect to anomalous vs. non-anomalous waveforms, which then would relate to the overall ranking of the area/areas.

It is particularly noted that the signal waveform is a series of continuous changes where the same attributes that determine the anomalous areas may also determine the non-anomalous areas. Signal development throughout the length of the signal is composed of a series of changes of the individual attributes - however, the primary difference between an area and waveform marked as anomalous as compared to an area or waveform not considered anomalous is the change sequence for the individual attributes. Signal development over the course of a spatial interval will generally be composed of a number of possible sequential changes in the attributes and the starting point is usually the identification of the start and end of a particular development subsection. This is referred to as the determination of the "time windows".

A feature of one embodiment described is that it includes dynamic signal integration while in motion. All present gravimeters and gravity gradiometers take static measurements, i.e., sensor is not required to be in motion and, therefore, these instruments (1 ) are incapable of performing a dynamic signal integration or detecting the energy accumulation over time; (2) have no ability to selectively respond to subsurface stress induced energy changes (perturbations from each subsurface discontinuity); (3) are incapable of sampling vibration signals of the gravitational field at a high sample rate (e.g., 2,000 samples per second); and (4) have no ability to identify the presence of localized trapped fluids in the subsurface.

In one embodiment, there is a gravity transducer system comprising: a junction comprising a first metal and a second metal different than the first metal, the junction forming a resonant cavity; a particle characterized by a de Broglie wave, the particle contained within the resonant cavity; a phonon wave within the resonant cavity, the de Broglie wave coupled with the phonon wave; a power source for applying electrical power across the junction; a sensor capable of sensing resistance, voltage, or current, the sensor electrically connected with the junction, the sensor producing a sensed signal; and a recording system for recording the sensed signal. In one alternative, the particle is not held in a lattice structure. In another alternative, the system further comprises an accelerometer. In one alternative, the recording system comprises a geographical positioning system (GPS). In one embodiment, the system further comprises a vehicle capable of moving the transducer system. In one alternative, the vehicle is an airplane. In one alternative the recording system comprises a mapping system for mapping the location of the vehicle. In one embodiment, the junction comprises: a conductive bead having longitudinal bore through it; and a conductive pin having a dielectric coating, the pin located in the longitudinal bore to form the electrical junction. In one alternative, the bead and pin are made of metal, and the dielectric comprises a metal oxide. In one embodiment, the metal oxide comprises silver oxide. In one embodiment, the bead comprises a sandwich comprising a copper or gold layer formed between first and second alloy layers. In another embodiment, the power source is a voltage source and the voltage is between 3 volts and 9 volts.

In another embodiment, there is a method of detecting orientation changes or localized variations of a gravitational field associated with sub-surface stress, the method comprising: creating a resonant cavity; trapping a particle in the resonant cavity; applying an electrical field across the particle in the resonant cavity to create an electrical junction; exposing the junction to the orientation changes or localized variations of a gravitational field; sensing the resistance, voltage, or current across the junction to produce a sensed signal characteristic of the orientation changes or localized variations of a gravitational field associated with sub-surface stress; and recording the sensed signal; and analyzing the signal to detect the orientation changes or localized variations of a gravitational field. 11 n one embodiment, the method further comprises, moving the trapped particle in the resonant cavity through the gravitational field. In one alternative, the moving comprises moving the trapped particle in the resonant cavity in a vehicle. In another embodiment, the moving comprises moving the trapped particle in the resonant cavity in an airplane. In one alternative, the recording comprises recording the sensed signal as a function of time. In another embodiment, the recording comprises recording the sensed signal as a function of geographical position of the trapped particle in the resonant cavity. In one embodiment, the analyzing comprises analyzing the signal amplitude, signal frequency, or signal character and pattern of the signal.

In another embodiment, there is a method of detecting a hydrocarbon deposit, the method comprising: providing a gravity transducer for sensing a change in a gravity field; moving the gravity transducer through the atmosphere or on land above the hydrocarbon deposit; sensing a change in the gravity field to produce a gravity change signal indicative of the hydrocarbon deposit; recording the gravity change signal; and analyzing the gravity change signal to detect the hydrocarbon deposit. In one alternative, the moving is performed by placing the gravity transducer in a vehicle and moving the vehicle. In one embodiment, the method further comprises tracking the position of the vehicle and the recording comprises recording the position of the vehicle. In one alternative, the analyzing comprises correlating the recorded gravity change signal with the position of the vehicle. In one embodiment, the moving comprises placing the transducer in an aircraft and flying the aircraft over the deposit. In one embodiment, the method further comprises tracking the geographical position of the aircraft and recording the tracked position. In another embodiment, the method further comprises displaying the real time position of the aircraft. In another embodiment, the displaying comprises displaying a planned flight line and deviations from the flight line. In another embodiment, the displaying comprises displaying the planned flight line and the deviations on the windshield of the aircraft. In another

embodiment, the method further comprises monitoring rotation of the aircraft along a horizontal plane to provide a rotation signal, the recording comprises recording of the rotation signal, and the analyzing comprises using the rotation signal to improve the quality of the gravity change signal. In one embodiment, the flying comprises flying the aircraft at an altitude of between 1 ,000 feet and 20,000 feet. In another embodiment, the method comprises integrating the gravity change signal dynamically over time. In another embodiment, the method further comprises detecting the energy accumulation in the signal over time. In another embodiment, the signal is sampled at a sample rate of 1000 samples per second or more.

There is also a method of detecting a subterranean fluid trapped beneath a layer of earth, the method comprising: providing a subterranean fluid transducer capable of detecting the trapped subterranean fluid; moving the transducer through the atmosphere or on land above the trapped fluid; sensing the trapped fluid to produce a trapped fluid signal indicative of the trapped fluid; recording the trapped fluid signal; and analyzing the trapped fluid signal to detect the trapped fluid. In one embodiment, the moving is performed by placing the transducer in a vehicle and moving the vehicle. In another embodiment, the moving comprises placing the transducer in an aircraft and flying the aircraft over the deposit. In another embodiment, the method further comprises tracking the geographical position of the aircraft and recording the tracked position. In another embodiment, the method further comprises monitoring rotation of the aircraft along a horizontal plane to provide a rotation signal, the recording comprises recording of the rotation signal, and the analyzing comprises using the rotation signal to improve the quality of the gravity change signal. In another embodiment, the comprises flying the aircraft at an altitude of between 1 ,000 feet and 20,000 feet. In one embodiment, the trapped fluid is a hydrocarbon.

There has been described a practical gravitational detector, which detector is stable, compact, rugged and can be operated at room temperature.

Any of the parts of any one of the embodiments described herein may be combined with any of the parts of any of the other embodiments. Equivalent structures and processes may be substituted for the various structures and processes described; the sub-processes of the inventive method may, in some instances, be performed in a different order; or a variety of different materials and elements may be used. Consequently, the invention is to be construed as embracing each and every novel feature and novel combination of features present in and/or possessed by the gravity transducer apparatus and methods described.

Claims

CLAIMS We claim:
1 . A gravity transducer system comprising:
a junction comprising a first metal and a second metal different than said first metal, said junction forming a resonant cavity;
a particle characterized by a de Broglie wave, said particle contained within said resonant cavity;
a phonon wave within said resonant cavity, said de Broglie wave coupled with said phonon wave;
a power source for applying electrical power across said junction;
a sensor capable of sensing resistance, voltage, or current, said sensor electrically connected with said junction, said sensor producing a sensed signal; and
a recording system for recording said sensed signal.
2. A gravity transducer system as in claim 1 wherein said particle is not held in a lattice structure.
3. A gravity transducer system as in claim 1 and further comprising an accelerometer.
4. A gravity transducer system as in claim 1 wherein said recording system comprises a geographical positioning system (GPS).
5. A gravity transducer system as in claim 1 and further comprising a vehicle capable of moving said transducer system.
6. A gravity transducer as in claim 1 wherein said vehicle is an airplane.
7. A gravity transducer system as in claim 5 wherein said recording system comprises a mapping system for mapping the location of said vehicle.
8. A gravity transducer system as in claim 1 wherein said junction comprises
a conductive bead having longitudinal bore through it; and
a conductive pin having a dielectric coating, said pin located in said longitudinal bore to form said electrical junction.
9. A gravity transducer system as in claim 8 wherein said bead and pin are made of metal, and said dielectric comprises a metal oxide.
10. A gravity transducer system as in claim 9 wherein said metal oxide comprises silver oxide.
1 1 . A gravity transducer system as in claims 8 wherein said bead comprises a sandwich comprising a copper or gold layer formed between first and second alloy layers.
12. A gravity transducer system as in claim 1 wherein said power source is a voltage source and said voltage is between 3 volts and 9 volts.
13. A method of detecting orientation changes or localized variations of a gravitational field associated with sub-surface stress, said method comprising: creating a resonant cavity;
trapping a particle in said resonant cavity;
applying an electrical field across said particle in said resonant cavity to create an electrical junction;
exposing said junction to said orientation changes or localized variations of a gravitational field;
sensing the resistance, voltage, or current across said junction to produce a sensed signal characteristic of said orientation changes or localized variations of a gravitational field;
recording said sensed signal; and analyzing said signal to detect said orientation changes or localized variations of a gravitational field associated with sub-surface stress.
14. A method as in claim 13 and further comprising moving said trapped particle in said resonant cavity through said gravitational field.
15. A method as in claim 14 wherein said moving comprises moving said trapped particle in said resonant cavity in a vehicle.
16. A method as in claim 15 wherein said moving comprises moving said trapped particle in said resonant cavity in an airplane.
17. A method as in claim 13 wherein said recording comprises recording said sensed signal as a function of time.
18. A method as in claim 13 wherein said recording comprises recording said sensed signal as a function of geographical position of said trapped particle in said resonant cavity.
19. A method as in claim 13 wherein said analyzing comprises analyzing the signal amplitude, signal frequency, or signal character and pattern of said signal.
20. A method of detecting a hydrocarbon deposit, said method comprising:
providing a gravity transducer for sensing a change in a gravity field;
moving said gravity transducer through the atmosphere or on land above said hydrocarbon deposit;
sensing a change in said gravity field to produce a gravity change signal indicative of said hydrocarbon deposit;
recording said gravity change signal; and
analyzing said gravity change signal to detect said hydrocarbon deposit.
21 . A method as in claim 20 wherein said moving is performed by placing said gravity transducer in a vehicle and moving said vehicle.
22. A method as in claim 21 and further comprising tracking the position of said vehicle and said recording comprises recording said position of said vehicle.
23. A method as in claim 22 wherein said analyzing comprises correlating said recorded gravity change signal with said position of said vehicle.
24. A method as in claim 21 wherein said moving comprises placing said transducer in an aircraft and flying said aircraft over said deposit.
25. A method as in claim 24 and further comprising tracking the geographical position of said aircraft and recording said tracked position.
26. A method as in claim 25 and further comprising displaying the real time position of said aircraft.
27. A method as in claim 26 wherein said displaying comprising displaying a planned flight line and deviations from said flight line.
28. A method as in claim 27 wherein said displaying comprises displaying said planned flight line and said deviations on the windshield of said aircraft.
29. A method as in claim 24 and further comprising monitoring rotation of said aircraft along a horizontal plane to provide a rotation signal, said recording comprises recording of said rotation signal, and said analyzing comprises using said rotation signal to improve the quality of said gravity change signal.
30. A method as in claim 24 wherein said flying comprises flying said aircraft at an altitude of between 1 ,000 feet and 20,000 feet.
31 . A method as in claim 20 and further comprising integrating said gravity change signal dynamically over time.
32. A method as in claim 20 and further comprising detecting the energy accumulation in said signal over time.
33. A method as in claim 20 wherein said signal is sampled at a sample rate of 1000 samples per second or more
34. A method of detecting a subterranean fluid trapped beneath a layer of earth, said method comprising:
providing a subterranean fluid transducer capable of detecting said trapped subterranean fluid;
moving said transducer through the atmosphere or on land above said trapped fluid;
sensing said trapped fluid to produce a trapped fluid signal indicative of said trapped fluid;
recording said trapped fluid signal; and
analyzing said trapped fluid signal to detect said trapped fluid.
35. A method as in claim 34 wherein said moving is performed by placing said transducer in a vehicle and moving said vehicle.
36. A method as in claim 34 wherein said moving comprises placing said transducer in an aircraft and flying said aircraft over said deposit.
37. A method as in claim 34 and further comprising tracking the geographical position of said aircraft and recording said tracked position.
38. A method as in claim 36 and further comprising monitoring rotation of said aircraft along a horizontal plane to provide a rotation signal, said recording comprises recording of said rotation signal, and said analyzing comprises using said rotation signal to improve the quality of said gravity change signal.
39. A method as in claim 36 wherein said flying comprises flying said aircraft at an altitude of between 1 ,000 feet and 20,000 feet.
A method as in claim 34 wherein said trapped fluid is a hydrocarbon.
PCT/US2013/042315 2012-05-22 2013-05-22 Gravity transducer and application to hydrocarbon exploration WO2013177340A1 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US201261650413P true 2012-05-22 2012-05-22
US61/650,413 2012-05-22

Applications Claiming Priority (10)

Application Number Priority Date Filing Date Title
BR112014029304A BR112014029304A2 (en) 2012-05-22 2013-05-22 transducer gravity and application for hydrocarbon exploration
CA2874450A CA2874450C (en) 2012-05-22 2013-05-22 Gravity transducer
MX2014014309A MX350913B (en) 2012-05-22 2013-05-22 Gravity transducer and application to hydrocarbon exploration.
CN201380039132.9A CN104685380B (en) 2012-05-22 2013-05-22 Converter and gravity of hydrocarbon exploration applications
EP13794555.6A EP2852854A4 (en) 2012-05-22 2013-05-22 Gravity transducer and application to hydrocarbon exploration
JP2015514168A JP6175131B2 (en) 2012-05-22 2013-05-22 Gravity transducer
RU2014151766A RU2608632C2 (en) 2012-05-22 2013-05-22 Gravity transducer and application to hydrocarbon exploration
US14/550,779 US9817153B2 (en) 2012-05-22 2014-11-21 Gravity transducer system and method including a junction with a first metal and a second metal
IN2885/KOLNP/2014A IN2014KN02885A (en) 2012-05-22 2014-12-10 Gravity transducer and application to hydrocarbon exploration
US15/727,157 US10078152B2 (en) 2012-05-22 2017-10-06 Gravity transducer system and method including junctions with a first metal and a second metal

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US14/550,779 Continuation-In-Part US9817153B2 (en) 2012-05-22 2014-11-21 Gravity transducer system and method including a junction with a first metal and a second metal

Publications (1)

Publication Number Publication Date
WO2013177340A1 true WO2013177340A1 (en) 2013-11-28

Family

ID=49624328

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2013/042315 WO2013177340A1 (en) 2012-05-22 2013-05-22 Gravity transducer and application to hydrocarbon exploration

Country Status (9)

Country Link
EP (1) EP2852854A4 (en)
JP (1) JP6175131B2 (en)
CN (1) CN104685380B (en)
BR (1) BR112014029304A2 (en)
CA (1) CA2874450C (en)
IN (1) IN2014KN02885A (en)
MX (1) MX350913B (en)
RU (1) RU2608632C2 (en)
WO (1) WO2013177340A1 (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104618033A (en) * 2015-01-08 2015-05-13 东南大学 Multi-layer self-adapting morphological filtering gravity signal noise inhibition method
WO2016077915A1 (en) * 2014-11-21 2016-05-26 Nxt Energy Solutions, Inc. Gravity transducer
US9817153B2 (en) 2012-05-22 2017-11-14 Nxt Energy Solutions, Inc. Gravity transducer system and method including a junction with a first metal and a second metal
US10078152B2 (en) 2012-05-22 2018-09-18 Nxt Energy Solutions, Inc. Gravity transducer system and method including junctions with a first metal and a second metal

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4920313A (en) * 1979-01-11 1990-04-24 Constant James N Gravitational mass detector
US4992656A (en) 1987-10-26 1991-02-12 Clauser John F Rotation, acceleration, and gravity sensors using quantum-mechanical matter-wave interferometry with neutral atoms and molecules
US20020092350A1 (en) * 2000-11-28 2002-07-18 Business Arts Inc. Gravity gradiometry
US20040000910A1 (en) * 2002-06-28 2004-01-01 Tryggvason Bjarni V. System and method for surveying underground density distributions
US20050197773A1 (en) 2004-03-05 2005-09-08 Bell Geospace Inc. Method and system for evaluating geophysical survey data
US20070010946A1 (en) 2005-07-08 2007-01-11 Bell Geospace Inc. Method and system for geophysical data acquisition on an airship
US20080020935A1 (en) * 2006-07-22 2008-01-24 Boris Volfson Phonon maser
US20080163682A1 (en) * 2005-10-06 2008-07-10 Technological Resources Pty. Limited Gravity Gradiometer
US20110012693A1 (en) * 2009-06-08 2011-01-20 Stmicroelectronics S.A. Bulk-mode resonator

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4920313A (en) * 1979-01-11 1990-04-24 Constant James N Gravitational mass detector
US4992656A (en) 1987-10-26 1991-02-12 Clauser John F Rotation, acceleration, and gravity sensors using quantum-mechanical matter-wave interferometry with neutral atoms and molecules
US20020092350A1 (en) * 2000-11-28 2002-07-18 Business Arts Inc. Gravity gradiometry
US20040000910A1 (en) * 2002-06-28 2004-01-01 Tryggvason Bjarni V. System and method for surveying underground density distributions
US20050197773A1 (en) 2004-03-05 2005-09-08 Bell Geospace Inc. Method and system for evaluating geophysical survey data
US20070010946A1 (en) 2005-07-08 2007-01-11 Bell Geospace Inc. Method and system for geophysical data acquisition on an airship
US20080163682A1 (en) * 2005-10-06 2008-07-10 Technological Resources Pty. Limited Gravity Gradiometer
US20080020935A1 (en) * 2006-07-22 2008-01-24 Boris Volfson Phonon maser
US20110012693A1 (en) * 2009-06-08 2011-01-20 Stmicroelectronics S.A. Bulk-mode resonator

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
B. M. ANDERSON; J. M. TAYLOR; V. M. GALITSKI: "Interferometry with synthetic gauge fields", PHYSICAL REVIEW A, vol. 83, 2011, pages 031602
MOODY M VOL ET AL.: "Review of Scientific Instruments", vol. 73, 11 November 2002, AIP, article "Three-axis superconducting gravity gradiometer for sensitive gravity experiments", pages: 3957 - 3974
R COLELLA ET AL.: "Observation of gravitationally induced quantum interference", PHYSICAL REVIEW LETTERS, 9 June 1975 (1975-06-09)
See also references of EP2852854A4 *
T. J. HAYES ET AL.: "A gravity gradient method for characterizing the post-seismic deformation field for a finite fault", GEOPHYSICAL JOURNAL INTERNATIONAL, vol. 173, 1 June 2008 (2008-06-01), pages 802 - 805, XP055327234, DOI: doi:10.1111/j.1365-246X.2008.03795.x

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9817153B2 (en) 2012-05-22 2017-11-14 Nxt Energy Solutions, Inc. Gravity transducer system and method including a junction with a first metal and a second metal
US10078152B2 (en) 2012-05-22 2018-09-18 Nxt Energy Solutions, Inc. Gravity transducer system and method including junctions with a first metal and a second metal
WO2016077915A1 (en) * 2014-11-21 2016-05-26 Nxt Energy Solutions, Inc. Gravity transducer
CN104618033A (en) * 2015-01-08 2015-05-13 东南大学 Multi-layer self-adapting morphological filtering gravity signal noise inhibition method
CN104618033B (en) * 2015-01-08 2015-09-09 东南大学 A multilayer gravity adaptive morphological filtering signal noise suppression method

Also Published As

Publication number Publication date
IN2014KN02885A (en) 2015-05-08
CN104685380B (en) 2018-04-13
BR112014029304A2 (en) 2019-01-02
MX350913B (en) 2017-09-25
CN104685380A (en) 2015-06-03
RU2014151766A (en) 2016-07-10
EP2852854A4 (en) 2017-01-18
JP6175131B2 (en) 2017-08-02
EP2852854A1 (en) 2015-04-01
CA2874450A1 (en) 2013-11-28
MX2014014309A (en) 2015-09-21
RU2608632C2 (en) 2017-01-23
JP2015520858A (en) 2015-07-23
CA2874450C (en) 2017-08-01

Similar Documents

Publication Publication Date Title
Hilterman Three-dimensional seismic modeling
CA2488511C (en) System and method for surveying underground density distributions
Artru et al. Acoustic waves generated from seismic surface waves: propagation properties determined from Doppler sounding observations and normal-mode modelling
Wu et al. Magnitude estimation using the first three seconds P‐wave amplitude in earthquake early warning
Gomberg et al. Stress/strain changes and triggered seismicity at The Geysers, California
US8566036B2 (en) Method and system for passive electroseismic surveying
Liu et al. Comparison of phase velocities from array measurements of Rayleigh waves associated with microtremor and results calculated from borehole shear-wave velocity profiles
Chávez-García et al. An alternative approach to the SPAC analysis of microtremors: exploiting stationarity of noise
Bell et al. Gravity gradiometry resurfaces
Paterson et al. Applications of gravity and magnetic surveys: The state-of-the-art in 1985
CA2252594C (en) Gravity meter
Langbein et al. High‐rate real‐time GPS network at Parkfield: Utility for detecting fault slip and seismic displacements
Telford et al. Applied geophysics
Graizer Tilts in strong ground motion
US7069780B2 (en) Gravity techniques for drilling and logging
Ivanova et al. Monitoring and volumetric estimation of injected CO2 using 4D seismic, petrophysical data, core measurements and well logging: a case study at Ketzin, Germany
Nabighian et al. Historical development of the gravity method in explorationHistorical Development of Gravity Method
Priestley et al. The African upper mantle and its relationship to tectonics and surface geology
AU778333B2 (en) A method for monitoring seafloor subsidence and for gravity monitoring an underground hydrocarbon reservoir
US2963641A (en) Exploration of earth formations associated with petroleum deposits
Zadro et al. Measurements and interpretations of tilt–strain gauges in seismically active areas
Jekeli A review of gravity gradiometer survey system data analyses
US20060036367A1 (en) Method and system for processing geophysical survey data
Mandea et al. Geomagnetic jerks: rapid core field variations and core dynamics
US7554883B2 (en) Fault filter for seismic discontinuity data

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 13794555

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase in:

Ref document number: 2874450

Country of ref document: CA

Ref document number: 2015514168

Country of ref document: JP

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: MX/A/2014/014309

Country of ref document: MX

NENP Non-entry into the national phase in:

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 2013794555

Country of ref document: EP

ENP Entry into the national phase in:

Ref document number: 2014151766

Country of ref document: RU

Kind code of ref document: A

REG Reference to national code

Ref country code: BR

Ref legal event code: B01A

Ref document number: 112014029304

Country of ref document: BR

REG Reference to national code

Ref country code: BR

Ref legal event code: B01E

Ref document number: 112014029304

Country of ref document: BR

REG Reference to national code

Ref country code: BR

Ref legal event code: B01E

Ref document number: 112014029304

Country of ref document: BR

Free format text: EM ADITAMENTO A EXIGENCIA PUBLICADA NA RPI 2428 DE 18/07/2017, A ENTRADA NA FASE NACIONAL SE DEU EM DISSONANCIA COM O PEDIDO INTERNACIONAL QUE CONSTAVA COMO DEPOSITANTES ?NXT ENERGY SOLUTIONS, INC? E ?CARL A FOREST?. TENDO EM VISTA A RESPOSTA DE EXIGENCIA, ATRAVES DA PETICAO NO 870170068510, A CESSAO DE PCT APRESENTADA ESTA IRREGULAR POIS POSSUI DATA POSTERIOR A ENTRADA NA FASE NACIONAL. ISTO POSTO: 1) APRESENTE CESSAO DO PEDIDO INTERNACIONAL PCT COM DATA ANTERIOR A ENTRADA NA FASE NACIONAL; OU 2) RETIFIQUE A ENTRADA NA FASE NACIONAL, CONSTANDO AMBOS COMO DEPOSITANTES (APRESENTANDO A PROCURACAO DE ?CARL A FOREST?), E SOLICITE A TRANSFERENCIA DO PEDIDO COM GRU E RITO PROPRIOS.

ENP Entry into the national phase in:

Ref document number: 112014029304

Country of ref document: BR

Kind code of ref document: A2

Effective date: 20141124