US7769572B2 - Method of imaging subsurface formations using a virtual source array - Google Patents

Method of imaging subsurface formations using a virtual source array Download PDF

Info

Publication number
US7769572B2
US7769572B2 US10/228,453 US22845302A US7769572B2 US 7769572 B2 US7769572 B2 US 7769572B2 US 22845302 A US22845302 A US 22845302A US 7769572 B2 US7769572 B2 US 7769572B2
Authority
US
United States
Prior art keywords
source
virtual
radial
bipole
dipole
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active, expires
Application number
US10/228,453
Other languages
English (en)
Other versions
US20030050759A1 (en
Inventor
Leonard J. Srnka
James J. Carazzone
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
ExxonMobil Upstream Research Co
Original Assignee
ExxonMobil Upstream Research Co
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
Application filed by ExxonMobil Upstream Research Co filed Critical ExxonMobil Upstream Research Co
Priority to US10/228,453 priority Critical patent/US7769572B2/en
Assigned to EXXONMOBIL UPSTREAM RESEARCH COMPANY reassignment EXXONMOBIL UPSTREAM RESEARCH COMPANY ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: CARAZZONE, JAMES J., SRNKA, LEONARD J.
Publication of US20030050759A1 publication Critical patent/US20030050759A1/en
Priority to US12/796,374 priority patent/US8078439B2/en
Application granted granted Critical
Publication of US7769572B2 publication Critical patent/US7769572B2/en
Active legal-status Critical Current
Adjusted expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/08Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
    • G01V3/083Controlled source electromagnetic [CSEM] surveying
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/12Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves

Definitions

  • This invention relates to the field of controlled-source electromagnetic surveys for geophysical applications. More specifically, this invention relates to a method of imaging subsurface formations using a virtual source array.
  • Electromagnetic sources used in geophysical surveying generally produce unfocused regions of subsurface excitation. Unfocused excitations result both from the diffusive nature of low frequency electromagnetic wave propagation in the earth and from the use of relatively simple source configurations in the surveys. Unfocused excitations are a severe disadvantage when a target-oriented survey is needed, such as when it is desirable to probe the electrical properties of a prospective formation in a location where a complex geologic background may be present.
  • a method of source focusing is disclosed in U.S. patent application Ser. No. 09/656,191, filed Sep. 6, 2000 (hereinafter referred to as “Srnka”).
  • the method referred to as the remote reservoir resistivity mapping method (“R3M”) uses a spatially extended grounded electromagnetic source array that is positioned on the earth's surface, or near the seafloor, to focus electromagnetic energy on the subsurface target.
  • the source array consists of two concentric electrode rings whose radii are optimally determined from theoretical modeling to focus the subsurface excitation at the approximate depth of the subsurface target.
  • a limitation of the Srnka method is that for typical reservoir depths, which range from 500 meters to greater than 5000 meters, the total length of electrode wire that is required in the concentric electrode rings ranges from about 23.5 km to greater than 235 km.
  • Electrode wire raise numerous operational challenges. For example, difficulties in handling and deploying the electrode wires, and problems in ensuring accurate positioning of the wires in the desired electrode ring locations may be involved, and there is the need to avoid obstacles in the area of the desired ring locations, and the requirement that the electrodes be uniformly grounded along their lengths. Other problems influence the economics of the survey, including the large cost of the electrode wire, some of which is likely to be damaged or lost in each survey, and the large labor and logistics costs associated with mobilization, demobilization, deployment, and recovery of the wire.
  • FIG. 1 depicts both the R3M concentric ring method and this radial bipole alternative.
  • Concentric electrode rings 4 and 5 lie on surface 1 of earth 2 and have radii a and b, respectively, and a center 7 generally above reservoir 3 .
  • these concentric rings are replaced by radial electrode bipole sources 6 .
  • Srnka discloses that a minimum of six radial bipole sources 6 is preferable.
  • a signal equivalent to the signal that a receiver array would receive from the concentric ring source array is produced mathematically by summing the signals received by the receiver array from each of the bipole sources 6 (For convenience, FIG. 1 does not depict the location of the receiver array).
  • This alternative method which allows the simulation of a concentric ring source array, is referred to as the virtual concentric ring source method, or as the radial bipole method, and substantially reduces the required length of electrode wire as compared to the R3M physical concentric ring source method.
  • both methods have other inherent disadvantages, in particular for certain aspects of the subsurface resistivity-imaging problem. Specifically, both methods produce maximum responses near a pre-specified subsurface excitation focus depth (in other words, at the approximate depth of the target), but provide little differential sensitivity to features within an appreciable vertical range of that depth. More specifically, the ability to image separately features above or below the target a distance of less than five percent of the target's depth below the surface is limited. This limitation greatly hinders the imaging of stacked reservoirs. In addition, amplifier saturation effects can influence data quality and hinder analysis of such data.
  • a third issue is that both methods produce primarily vertical subsurface currents in the region of the target. These currents are not adequate for imaging reservoir electrical macro-anisotropy (vertical electrical isotropy), such as exists in many reservoirs due to the presence of shale interbeds. Imaging this anisotropy is essential for estimating reservoir net and gross volumes, and for determining hydrocarbon pore volumes.
  • This invention relates generally to a method of simulating the signal of an electromagnetic source using one or more dipole sources.
  • a dipole source is located at an excitation location corresponding to a segment of the electromagnetic source to be simulated.
  • the dipole source is activated, and an electromagnetic signal recorded at one or more receiver locations. This process is repeated for additional excitation locations corresponding to additional segments of the source being simulated.
  • the data from the sequence of dipole source excitation locations is processed to determine the simulated signal of the electromagnetic source.
  • FIG. 1 depicts a perspective view of the Srnka method.
  • FIG. 2 is a flow chart of the steps for implementation of a first embodiment of the present invention.
  • FIG. 3 depicts a perspective view of a first embodiment of the present invention.
  • FIG. 4 depicts a plan view of a second embodiment of the present invention.
  • FIG. 5 depicts a data acquisition procedure for a marine application of an embodiment of the present invention.
  • FIG. 6 depicts a plan view of a method for acquiring data for the marine application of FIG. 5 .
  • FIG. 7 depicts a plan view of an embodiment of the present invention in which variations in the locations of the target depth and focal point can be made.
  • FIG. 8 depicts a model target for use with a bipole source processing example.
  • FIG. 9 depicts the total electric field calculated for the bipole source processing example of FIG. 8 .
  • FIG. 10 depicts a model target for use with an embodiment of the present invention in which dipole sources are used to simulate the data acquired by the bipole source example of FIG. 8 .
  • FIG. 11 depicts the total electric field calculated according to the embodiment of the present method depicted in FIG. 10 .
  • the present invention is a method of simulating the advantageous multicomponent electromagnetic response of a resistive reservoir to excitation by a concentric ring source array, or its radial bipole equivalent, with a significantly reduced amount of electrode wire.
  • the method uses one or more short dipole sources, each located sequentially at source excitation locations corresponding to a location of a radial bipole to be simulated, to obtain measurement data during field surveys in the region of the subsurface target to be imaged.
  • Phase and amplitude receiver data are synchronized with data from the source locations and summed using a Green's function integration to determine the response of the target.
  • the method has an inherent ability to create a large number of virtual source configurations by varying the number and positions of the source excitation locations that are integrated in the receiver data, thereby facilitating much higher spatial resolution of the near surface, as well as for other depths, than can be attained by previously proposed methods.
  • Additional virtual source configurations can be created by use of Green's function reciprocity conditions which can be used to relate electric field recordings made with exchanged transmitter and detector locations.
  • the present method allows for the recording of signals from a large number of linear electric sources. Integration of receiver data from these linear sources facilitates the imaging of subsurface resistivity outside of the region that would be the focal point of a concentric ring source array. This advantage allows for measurement of the subsurface horizontal resistivity response over a range of depths, and allows for detection of reservoir resistivity anisotropy.
  • An advantage of the use of a multiplicity of sources is production of a multiplicity of electric current probes of the subsurface target, which can have positive effects on data redundancy.
  • the method of the present invention can be more fully described with reference to the flow chart in FIG. 2 and the embodiment depicted in FIG. 3 .
  • the initial step, FIG. 2-step 100 requires calculation of the dimensions of the R3M method's concentric electrode rings for the subsurface target of interest.
  • subsurface target 3 within earth 2 is at depth D below the earth's surface 1
  • concentric electrode rings 4 ′ and 5 ′ have diameters a and b respectively.
  • Rings 4 ′ and 5 ′ are not physical electrode rings, but rather are solely used in the calculations required in the present method (As used herein, the addition of a prime to a reference number indicates that the referenced item is part of a virtual source, i.e.
  • a source to be simulated as opposed to a physical source.
  • electrode rings 4 and 5 are physical sources
  • electrode rings 4 ′ and 5 ′ are virtual sources.
  • Srnka also discloses that a minimum ratio b/a ⁇ 4 provides a sufficient source excitation at depth D.
  • FIG. 2-step 102 the number, length L, and location of radial bipole sources for Srnka's virtual concentric ring source method must be determined.
  • the use of these bipole sources provide a mechanism for creating an electromagnetic source focal point at the depth D of target 3 .
  • the radial bipole sources e.g. the radial bipole sources 6 of FIG. 1
  • FIG. 2-step 102 results in the calculation of the number, length and location of the radial bipole sources to be simulated. For simplicity, the location of only one simulated radial bipole source 6 ′ is shown in FIG. 3 .
  • Length l is restricted to l ⁇ a so that the radius of the inner virtual concentric ring 4 ′ is at least one radial dipole source length, for ease of data integration following the steps described below.
  • the number N may range from about 3 to more than 300.
  • FIG. 2-step 106 the position of receiver array 12 is determined. Factors to be considered in designing and determining a location for the receiver array, which will be understood to those skilled in the art, are further discussed in Srnka.
  • FIG. 2-step 108 a first radial dipole source is placed at a first source excitation location corresponding to a first bipole source to be simulated.
  • the first bipole source to be simulated is 6 ′
  • radial dipole source 8 of length l is placed at first source excitation location 10 1 along a radius 9 of the virtual outer concentric ring 5 ′.
  • first source excitation location 10 1 extends from virtual outer concentric ring 5 ′ along radius 9 inward toward center 7 .
  • the present method is not limited to this location of first source excitation location 10 1 .
  • radial dipole source 8 should be grounded at each end over a distance of at least ten percent of length l, and, as will be understood to those skilled in the art, should be powered by variable frequency (preferably 10 ⁇ 4 -10 4 Hz), high current (preferably 10 2 -10 6 Amperes) electrical power sources and controllers. Further details on sources for use in electromagnetic surveying are provided by Srnka.
  • FIG. 2-step 110 the physical location of the radial dipole source 8 at the source excitation location 10 1 must be determined. This location can be determined by geodetic or other methods well known to those skilled in the art.
  • the required positioning accuracies for each end of the radial dipole source are preferably 1 meter on land and 5 to 10 meters offshore.
  • FIG. 2-step 112 the radial dipole source is energized and the amplitude and phase of the source signal is recorded.
  • the source signal at each source excitation location is preferably held constant, with variations in phase limited to 0.1 degree and in amplitude to 0.1%, as further described in Srnka.
  • step 114 the amplitude and phase of the response at the receiver array 12 are recorded.
  • FIG. 2-step 116 The sequence encompassed by FIG. 2-steps 108 through 114 is repeated, FIG. 2-step 116 , for a radial dipole source located at a second source excitation location.
  • the second source excitation location is indicated as 10 2 .
  • the process of steps 108 through 114 is repeated for additional source excitation locations 10 3 , 10 4 , . . . 10 N ⁇ 1 until measurements have been recorded along the entire length of the radial bipole source 6 ′ being simulated. In FIG. 3 , this length corresponds to the portion of radius 9 extending from virtual electrode ring 5 ′ to virtual electrode ring 4 ′.
  • This process completes the source excitations and receiver response recordings necessary to simulate radial bipole source 6 ′.
  • the mathematical calculations of the present method require the normalizing of the multicomponent receiver responses to the source electric dipole signal (magnitude and phase) recorded for each source excitation location. Normalization assures that unwanted changes in receiver signals due to variations in either source coupling or power source current are removed from the receiver data.
  • Source phase and amplitude control, as described above, together with recording of the source output with accuracies of 0.1 degree in phase and 0.1% in amplitude is preferred.
  • Additional processing steps can include manipulation of the current in the transmitting electric dipole to achieve noise suppression and focusing effects. Manipulation of the dipole current can be achieved during data recording or by computational means after data is recorded.
  • FIG. 2-step 118 The sequence encompassed by FIG. 2-steps 108 through 116 is then repeated, FIG. 2-step 118 , for the second and each subsequent bipole source locations, as specified in FIG. 2-step 102 (the second and subsequent bipole sources to be simulated are not shown in FIG. 3 ). This process completes the recording of the data required in the present method to compute the response of target 3 in FIG. 3 .
  • FIG. 2-step 120 the received signals for dipole source excitations corresponding to each bipole source are integrated in either the frequency or time domain using received signal amplitude and phase control relative to the source signal.
  • This process is referred to as amplitude and phase normalized Green's function integration.
  • the integration results in the creation of a virtual radial bipole source at each of the locations specified in FIG. 2-step 102 .
  • the virtual concentric ring source method of Srnka can be applied to determine the characteristics of the target.
  • a second novel aspect of the present invention involves the ability to create a large number of virtual concentric ring source array configurations (in other words, to obtain simulated data for virtual concentric rings with a range of values of a and b).
  • This advantage derives from the ability to vary the number and locations of source excitations that are summed in the receiver data, and allows accurate target characterizations over a broader range of depths than can be characterized by the R3M or virtual concentric ring source array methods.
  • FIG. 4 demonstrates this aspect of the present invention.
  • FIG. 4 depicts a sequence of virtual outer concentric electrode rings 25 1 , 25 2 , 25 3 , . . . 25 m , with diameters b 1 , b 2 , b 3 , . . . b m , and a sequence of virtual inner concentric electrode rings 24 1 , 24 2 , 24 3 , . . . 24 m , with diameters a 1 , a 2 , a 3 , . . . a m . All rings are centered at 27 .
  • a first virtual bipole radial source to be simulated, extending from ring 24 1 to ring 25 1 is shown as 26 1 .
  • Dipole source excitation locations 30 1 , 30 2 , 30 3 , . . . 30 N extend inward from ring 25 1 , along radius 29 toward center 27 . (For convenience in discussing FIG. 4 , the prime reference, which otherwise signifies the distinction between a virtual and a physical source, has been omitted).
  • the diameters b 1 and a 1 of virtual concentric electrode rings 25 i and 24 1 and length l of the dipole sources to be used at each dipole source excitation location 30 1 are determined as discussed above.
  • the constraint l ⁇ a noted above will determine the radius for the innermost, inner concentric electrode ring 24 m . This constraint thereby establishes the total number of virtual concentric electrode rings in this embodiment of the present invention. Since N ⁇ 1 radial dipole positions exist which exclude one electrode end at zero radius, (in other words, for which the constraint l ⁇ a is satisfied), it follows from the well-known binomial theorem of mathematics that ((N ⁇ 1)!)/(2!(N ⁇ 3)! virtual ring source combinations exist.
  • These virtual ring source combinations cover a range of aspect ratios of the radius of the virtual outer concentric ring to the radius of the virtual inner concentric ring (in other words, of the ratio b 1 /a 1 ), including aspect ratios near unity that produce narrow (annular) and shallow subsurface excitations.
  • a broad range of virtual bipole sources 26 i can be simulated, resulting in a characterization of the subsurface region over a much broader range of depths than has previously been possible.
  • additional virtual source array configurations can be generated using Green's function symmetries.
  • Source and receiver data is recorded for each source excitation location extending along the entire length of radius 29 .
  • these source excitation locations extend from outer concentric ring 25 1 along radius 29 to the innermost inner concentric ring 24 m .
  • This data recording process continues for each radius 29 along which a virtual radial bipole source is to be simulated (only a first radius 29 is depicted in FIG. 4 ).
  • a minimum of six virtual bipole sources is preferred, separated by no more than 60 degrees.
  • the processing of data in this embodiment also follows from the discussion of FIGS. 2 and 3 above, except that the additional amount of data acquired in this embodiment adds additional processing steps.
  • This data processing sequence is repeated for the source excitation locations that correspond to each additional pair of concentric rings 25 i and 24 i , each pair of which facilitates the characterization of the subsurface at depth D 1 .
  • data may be processed for source excitation locations extending along radius 29 from ring 25 1 to 24 2 , thereby creating a distinct virtual bipole to be simulated, and a distinct subsurface target depth. Therefore, in the present method, data can be acquired once over the areal region on the surface of the earth overlying the target of interest, with multiple processing sequences used to increase the ability to characterize the entire subsurface below that areal region.
  • a preferred embodiment for land-based surveys uses only one electromagnetic dipole source that is moved sequentially in a radial direction inward, outward, or in an alternating sequence in both directions.
  • Use of only one dipole source reduces the total amount of source electrode wire from tens or hundreds of kilometers to only a few hundred meters, as is determined by the specified dipole length l.
  • the dipole length l for land-based surveys is preferably 500-1000 meters, to minimize the number of necessary source locations, since deployment of grounded source dipoles on land is relatively difficult.
  • the source dipole occupies positions 30 1 , 30 2 , 30 3 , and so on corresponding to a first bipole source to be simulated ( FIG. 4 ). The process is repeated for each bipole source to be simulated.
  • FIGS. 5 and 6 show sequential source dipole movement in the marine case.
  • the marine implementation of the present invention has advantages over the land-based survey implementation, including the ease of re-positioning the source dipole and the improved electrical coupling of the source electrodes to the seafloor and/or water so that the effort to ground the ends of the electrode is unnecessary.
  • FIG. 5 shows a survey vessel 35 at the surface 60 of a body of water and an electric source dipole 63 at or near, but preferably within 50 meters of, the bottom 61 of the body of water, such as the seafloor.
  • Source dipole 63 has a length l of typically 100-500 meters in the marine case, so as to minimize the difficulties of controlling the source remotely and accurately measuring its position from the survey vessel.
  • the survey is carried out by the vessel either: 1) sequentially moving the source from a first source excitation location 64 1 to a second source excitation location 64 2 , where sequential source excitation locations are centered at a distance separated by the length l of the source dipole, and then remaining stationary during the excitation of the source, with the process repeated for the source excitation locations required to simulate the bipole source of interest; or 2) continuously towing the dipole at low speed (typically 1 to 4 knots, approximately 0.5 to 2.0 meters/sec) while the source is under excitation.
  • the source towing speed In the towed dipole source case, the source towing speed must be sufficiently slow so that distortion of the source excitation waveform due to the Doppler effect is negligible.
  • v s 1 meter/sec
  • f 0.5 Hertz
  • 1 ⁇ -m
  • ⁇ f 0.045%
  • M values also referred to as data stacks
  • M M values
  • M the data stack value M ⁇ 50.
  • FIG. 6 A preferred method of acquiring data for the marine case is depicted in FIG. 6 .
  • Vessel 35 tows an electric dipole source 36 of length l along a first diameter 37 1 of virtual concentric ring 38 1 .
  • vessel 35 Upon completion of data acquisition along the entire length of diameter 37 1 , vessel 35 turns and commences data acquisition along a second diameter 37 2 , separated by angle ⁇ from first diameter 37 1 .
  • a series of diameters 37 1 (where the subscript i refers to each of the individual diameters along which data is to be acquired), preferably from 3 to 16, separated by equal angles ⁇ , allows formation of a series of virtual concentric source rings of diameter ratios b 1 /a i , as described above.
  • Model studies have shown that 3 to 16 equally spaced diameters are preferable to approximate concentric ring sources.
  • the present method allows for the simulation of a large number of linear sources. These simulations derive from the integration of the data using other combinations of signals resulting from the excitations by the short dipoles, and thus can be carried out at no additional survey cost.
  • This aspect of the invention allows for imaging of subsurface resistivity outside of the central focusing zones of the concentric ring sources, such as near the location of the outer concentric electrode ring which is simulated in the processing described above. For example, in FIG. 7 , a set of full-diameter bipole sources 40 1 , 40 2 , . . . 40 D , corresponding to virtual concentric rings 44 1 , 44 2 , . . .
  • each of these full-diameter bipole sources can be simulated using the method of the present invention, as described above. During the data acquisition phase of the present method, it is preferable that data are acquired between the innermost inner concentric ring 44 M and center 47 . This preference results in source and receiver data recordings for the entire length of each full-diameter bipole source.
  • the dipole source excitation locations at a first end of the simulated bipole source may be excluded from the Green's function integration, thereby shifting the symmetry point from center 47 along the simulated bipole source in a direction away from the end at which the sources were excluded.
  • This exclusion also reduces the effective depth of excitation, and allows for a different central focusing zone to be probed, in other words for a different subsurface region both in depth and in source focal point to be evaluated.
  • the source focal point will be directly below the midpoint of diameter 40 1 , in other words directly below center 47 .
  • the data corresponding to dipole source excitation location 40 1A is excluded from the Green's function integration. This exclusion shifts the source focal point to a point which is nevertheless below diameter 40 1 , but which is shifted horizontally by the distance 1 ⁇ 2 away from location 40 1A .
  • a second exclusion of the data corresponding to dipole source excitation location 40 1B will shift the source focal point an additional distance of 1 ⁇ 2 in that direction. Each of these exclusions will also reduce the depth of the source focal point.
  • This process provides a method for developing a broader characterization of the subsurface.
  • the horizontal electric field excitations resulting from exclusion of these individual dipole source locations in the integration when combined with the primarily vertical excitation due to the virtual concentric ring sources, allows for detection of reservoir resistivity anisotropy, as discussed in Srnka.
  • a hypothetical model reservoir target 50 is a 5000 meter ⁇ 5000 meter square slab 20 meters thick of 100 ⁇ -m (ohm-meter) resistive rock that is buried in a uniform background of 1 ⁇ -m conductive rock 51 .
  • the top of target 50 is buried 1000 meters below the seafloor, which in turn is covered by 1000 meters of uniform 0.25 ⁇ -m seawater.
  • An insulating air layer 5000 meters high is placed above the seawater.
  • a uniform grid of 5,329 electromagnetic receivers 52 spaced 250 meters apart in both horizontal (X and Y) directions, is placed at the seafloor, uniformly distributed within the coordinates ⁇ 9000 ⁇ Y ⁇ 9000 and ⁇ 9000 ⁇ X ⁇ 9000.
  • Bipole 53 is assumed to be placed in the seawater one meter above the seafloor and emits 100 Amperes of electrical current at a sinusoidal frequency of 1.0 Hertz.
  • FIG. 9 shows the magnitude of the total vertical electric field in Volts/meter at the receiver positions, calculated at a frequency of 1.0 Hertz, for the single bipole model of FIG. 8 .
  • the response is smoothly varying in X and Y, except for large, mathematically infinite, values at the two ends of the bipole where current is injected.
  • Persons skilled in the art will understand the source of these singularities in the numerical calculation procedure, and will further understand the decrease of the total vertical electric field with distance from source bipole 53 .
  • the ends of adjacent dipoles 60 are separated by one meter in the Y direction for more stable calculation in the numerical simulation (This separation is not depicted at the grid scale of this figure).
  • FIG. 11 shows the calculated response for the model of FIG. 10 using the present method, again at a frequency of 1.0 Hertz.
  • the response in FIG. 10 is substantially identical to the response of FIG.
  • FIG. 11 depicts the response calculated according to the present method for a model target of specified characteristics.
  • Data processing techniques will be required to determine the resistivity and depth characteristics of a target for which data have been acquired according to the present method. These processing techniques generally involve one of two methods. Curve matching techniques, in which modeled responses, such as in FIG. 11 , are compared to actual measured data to thereby infer the actual characteristics of the data acquisition target, are a first such method.
  • a second method involves use of an imaging algorithm such as described by Srnka in which a wave equation based algorithm is used to invert the acquired data into a model that contains the characteristics of the target.

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Electromagnetism (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)
US10/228,453 2001-09-07 2002-08-27 Method of imaging subsurface formations using a virtual source array Active 2028-12-30 US7769572B2 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US10/228,453 US7769572B2 (en) 2001-09-07 2002-08-27 Method of imaging subsurface formations using a virtual source array
US12/796,374 US8078439B2 (en) 2001-09-07 2010-06-08 Method of imaging subsurface formations using a virtual source array

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US31808201P 2001-09-07 2001-09-07
US10/228,453 US7769572B2 (en) 2001-09-07 2002-08-27 Method of imaging subsurface formations using a virtual source array

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US12/796,374 Continuation US8078439B2 (en) 2001-09-07 2010-06-08 Method of imaging subsurface formations using a virtual source array

Publications (2)

Publication Number Publication Date
US20030050759A1 US20030050759A1 (en) 2003-03-13
US7769572B2 true US7769572B2 (en) 2010-08-03

Family

ID=23236560

Family Applications (2)

Application Number Title Priority Date Filing Date
US10/228,453 Active 2028-12-30 US7769572B2 (en) 2001-09-07 2002-08-27 Method of imaging subsurface formations using a virtual source array
US12/796,374 Expired - Fee Related US8078439B2 (en) 2001-09-07 2010-06-08 Method of imaging subsurface formations using a virtual source array

Family Applications After (1)

Application Number Title Priority Date Filing Date
US12/796,374 Expired - Fee Related US8078439B2 (en) 2001-09-07 2010-06-08 Method of imaging subsurface formations using a virtual source array

Country Status (12)

Country Link
US (2) US7769572B2 (ru)
EP (1) EP1433090A1 (ru)
CN (1) CN1276377C (ru)
AU (1) AU2002329890B2 (ru)
BR (1) BR0211611A (ru)
CA (1) CA2458170C (ru)
EA (1) EA007644B1 (ru)
MA (1) MA26217A1 (ru)
MX (1) MXPA04002105A (ru)
MY (1) MY141136A (ru)
NO (1) NO20040964L (ru)
WO (1) WO2003025803A1 (ru)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080089526A1 (en) * 2006-10-16 2008-04-17 Siemens Audiologische Technik Gmbh Method for estimating an interference field for a coil
US20100017156A1 (en) * 2006-08-24 2010-01-21 Anton Ziolkowski Reduction of noise in electrical field measurements
WO2015038102A3 (en) * 2013-09-10 2015-11-19 Halliburton Energy Services Inc. Surface calibration of a wellbore resistivity logging tool

Families Citing this family (41)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2390904B (en) * 2002-07-16 2004-12-15 Univ Southampton Electromagnetic surveying for hydrocarbon reservoirs
GB2404444B (en) * 2003-07-28 2006-11-29 Statoil Asa Transmitter antena
GB2409900B (en) 2004-01-09 2006-05-24 Statoil Asa Processing seismic data representing a physical system
MXPA06009118A (es) * 2004-02-13 2007-01-26 Exxonmobil Upstream Res Co Sistema y metodo para remolcar antena vertical submarina.
GB2412741B (en) * 2004-04-03 2009-02-25 Statoil Asa Electromagnetic data processing
US7347271B2 (en) 2004-10-27 2008-03-25 Schlumberger Technology Corporation Wireless communications associated with a wellbore
US7477160B2 (en) 2004-10-27 2009-01-13 Schlumberger Technology Corporation Wireless communications associated with a wellbore
GB2420855B (en) 2004-12-02 2009-08-26 Electromagnetic Geoservices As Source for electromagnetic surveying
GB2422673B (en) * 2005-02-01 2010-03-24 Electromagnetic Geoservices As Optimum signal for sea bed logging
BRPI0607222A2 (pt) 2005-02-16 2009-08-18 Exxonmobil Upstream Res Co estimativa de ruìdo em uma freqüência por amostragem de ruìdo em outras freqüências
RU2284555C1 (ru) * 2005-06-01 2006-09-27 Николай Иванович РЫХЛИНСКИЙ Способ морской геоэлектроразведки с фокусировкой электрического тока (варианты)
US7894989B2 (en) * 2005-06-09 2011-02-22 Exxonmobil Upstream Research Co. Method for determining earth vertical electrical anisotropy in marine electromagnetic surveys
NO323889B3 (no) 2005-11-03 2007-07-16 Advanced Hydrocarbon Mapping As Framgangsmåte for kartlegging av hydrokarbonreservoarer samt apparat for anvendelse ved gjennomføring av framgangsmåten
GB2435693A (en) * 2006-02-09 2007-09-05 Electromagnetic Geoservices As Seabed electromagnetic surveying
CA2643057C (en) 2006-02-21 2015-05-26 Exxonmobil Upstream Research Company Method for electromagnetic air-wave suppression by active cancellation and shielding
US8014988B2 (en) * 2006-04-06 2011-09-06 Exxonmobil Upstream Research Co. Method for obtaining resistivity from controlled source electromagnetic data
AU2007248882B2 (en) * 2006-05-04 2011-01-27 Exxonmobil Upstream Research Company Time lapse analysis with electromagnetic data
NO327007B1 (no) * 2006-05-24 2009-04-06 Norsk Hydro As Fremgangsmate for elektromagnetisk geofysisk kartlegging av undersjoiske bergartsformasjoner
GB2439378B (en) 2006-06-09 2011-03-16 Electromagnetic Geoservices As Instrument for measuring electromagnetic signals
GB2442749B (en) 2006-10-12 2010-05-19 Electromagnetic Geoservices As Positioning system
GB2445582A (en) 2007-01-09 2008-07-16 Statoil Asa Method for analysing data from an electromagnetic survey
US7859942B2 (en) * 2007-03-01 2010-12-28 Christof Stork Measuring and modifying directionality of seismic interferometry data
WO2009006464A1 (en) * 2007-07-03 2009-01-08 Shell Oil Company Method of determining electrical anisotropy in a subsurface formation
US8358562B2 (en) * 2008-09-03 2013-01-22 Christof Stork Measuring and modifying directionality of seismic interferometry data
US8901931B2 (en) * 2009-03-13 2014-12-02 Schlumberger Technology Corporation Electromagnetic surface-to-borehole look around systems and methods of monitoring in horizontal wells
US9075163B2 (en) 2009-04-17 2015-07-07 Westerngeco L.L.C. Interferometric seismic data processing
US8554482B2 (en) * 2009-05-05 2013-10-08 Baker Hughes Incorporated Monitoring reservoirs using array based controlled source electromagnetic methods
AU2010262878B2 (en) * 2009-06-19 2014-06-12 Shell Internationale Research Maatschappij B.V. Multi-mode electromagnetic surveying
EP2550552A4 (en) * 2010-03-23 2013-09-11 Pan Geo Subsea Inc METHOD USING STATIONARY STAR-FORMED ANTENNA FOR HANDLING FOCUSED BEAM-FORMED FIELDS AND BEAM-ORIENTED ELECTROMAGNETIC SIGNALS FROM VERY FINE SEDIMENTARY STRATIGRAPHIC FORMATIONS DEEPLY LOCATED IN EARTH
CN102466822B (zh) * 2010-11-04 2013-09-04 中国石油天然气集团公司 一种海洋电磁勘探四极互组合布极方法
US20140153365A1 (en) * 2012-11-30 2014-06-05 Chevron U.S.A. Inc. System and method for producing local images of subsurface targets
US9268052B2 (en) * 2012-12-04 2016-02-23 Chevron U.S.A. Inc. Method and system of using a data weighted electromagnetic source for time-lapse monitoring of a reservoir production or hydraulic fracturing
US9377552B2 (en) * 2013-02-28 2016-06-28 Chevron U.S.A. Inc. System and method for detecting a fracture in a rock formation using an electromagnetic source
US9274241B2 (en) * 2013-03-14 2016-03-01 Pgs Geophysical As Method and system for suppressing swell-induced electromagnetic noise
US9611736B2 (en) * 2013-08-29 2017-04-04 Saudi Arabian Oil Company Borehole electric field survey with improved discrimination of subsurface features
US20140191760A1 (en) * 2014-01-17 2014-07-10 Bentsion Zinger Method and apparatus for suppression of the airwave in subsea exploration
US10126456B2 (en) * 2015-01-07 2018-11-13 The Regents Of The University Of California System and method for groundwater detection and evaluation
US10878142B2 (en) * 2015-10-26 2020-12-29 Pgs Geophysical As Bipole source modeling
EP3413092B1 (en) * 2017-06-08 2022-06-01 Total Se Method for evaluating a geophysical survey acquisition geometry over a region of interest, related process, system and computer program product
US11073637B2 (en) * 2018-10-04 2021-07-27 Saudi Arabian Oil Company Data structure for fast invasion percolation modeling software
CN117310824B (zh) * 2023-06-20 2024-07-12 北京自动化控制设备研究所 电偶源激发可控源电磁场时间序列正演建模方法及系统

Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3551797A (en) * 1966-09-29 1970-12-29 Schlumberger Prospection High frequency electromagnetic well logging methods and apparatus
US3727231A (en) 1971-09-24 1973-04-10 Collins Radio Co Collinear dipole antenna
US4247821A (en) 1978-07-20 1981-01-27 Commonwealth Scientific And Industrial Research Organization Transient electromagnetic prospecting apparatus possessing noise supression features
US4446434A (en) 1978-12-20 1984-05-01 Conoco Inc. Hydrocarbon prospecting method with changing of electrode spacing for the indirect detection of hydrocarbon reservoirs
US4535293A (en) 1979-07-09 1985-08-13 Compagnie Generale De Geophysique Method for geophysical prospecting using impulsine currents
US4617518A (en) 1983-11-21 1986-10-14 Exxon Production Research Co. Method and apparatus for offshore electromagnetic sounding utilizing wavelength effects to determine optimum source and detector positions
US4633182A (en) 1983-03-03 1986-12-30 Instytut Gornictwa Naftowego I Gazownictwa Method and system for direct prospecting of hydrocarbon deposits
US4875015A (en) 1987-07-20 1989-10-17 University Of Utah Research Institute Multi-array borehole resistivity and induced polarization method with mathematical inversion of redundant data
US4968940A (en) * 1987-10-30 1990-11-06 Schlumberger Technology Corporation Well logging apparatus and method using two spaced apart transmitters with two receivers located between the transmitters
US5563513A (en) 1993-12-09 1996-10-08 Stratasearch Corp. Electromagnetic imaging device and method for delineating anomalous resistivity patterns associated with oil and gas traps
US5594343A (en) * 1994-12-02 1997-01-14 Schlumberger Technology Corporation Well logging apparatus and method with borehole compensation including multiple transmitting antennas asymmetrically disposed about a pair of receiving antennas
RU2084929C1 (ru) 1993-03-24 1997-07-20 Владимир Сергеевич Могилатов Способ геоэлектроразведки
WO1998028638A1 (fr) 1995-03-06 1998-07-02 Alexandr Zhakferovich Nizamov Procede de recherche de produits petroliers dans le sol
WO2000013046A1 (en) 1998-08-28 2000-03-09 Den Norske Stats Oljeselskap A.S. Method and apparatus for determining the nature of subterranean reservoirs
US6304086B1 (en) * 1999-09-07 2001-10-16 Schlumberger Technology Corporation Method and apparatus for evaluating the resistivity of formations with high dip angles or high-contrast thin layers

Patent Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3551797A (en) * 1966-09-29 1970-12-29 Schlumberger Prospection High frequency electromagnetic well logging methods and apparatus
US3727231A (en) 1971-09-24 1973-04-10 Collins Radio Co Collinear dipole antenna
US4247821A (en) 1978-07-20 1981-01-27 Commonwealth Scientific And Industrial Research Organization Transient electromagnetic prospecting apparatus possessing noise supression features
US4446434A (en) 1978-12-20 1984-05-01 Conoco Inc. Hydrocarbon prospecting method with changing of electrode spacing for the indirect detection of hydrocarbon reservoirs
US4535293A (en) 1979-07-09 1985-08-13 Compagnie Generale De Geophysique Method for geophysical prospecting using impulsine currents
US4633182A (en) 1983-03-03 1986-12-30 Instytut Gornictwa Naftowego I Gazownictwa Method and system for direct prospecting of hydrocarbon deposits
US4617518A (en) 1983-11-21 1986-10-14 Exxon Production Research Co. Method and apparatus for offshore electromagnetic sounding utilizing wavelength effects to determine optimum source and detector positions
US4875015A (en) 1987-07-20 1989-10-17 University Of Utah Research Institute Multi-array borehole resistivity and induced polarization method with mathematical inversion of redundant data
US4968940A (en) * 1987-10-30 1990-11-06 Schlumberger Technology Corporation Well logging apparatus and method using two spaced apart transmitters with two receivers located between the transmitters
RU2084929C1 (ru) 1993-03-24 1997-07-20 Владимир Сергеевич Могилатов Способ геоэлектроразведки
US5563513A (en) 1993-12-09 1996-10-08 Stratasearch Corp. Electromagnetic imaging device and method for delineating anomalous resistivity patterns associated with oil and gas traps
US5594343A (en) * 1994-12-02 1997-01-14 Schlumberger Technology Corporation Well logging apparatus and method with borehole compensation including multiple transmitting antennas asymmetrically disposed about a pair of receiving antennas
WO1998028638A1 (fr) 1995-03-06 1998-07-02 Alexandr Zhakferovich Nizamov Procede de recherche de produits petroliers dans le sol
WO2000013046A1 (en) 1998-08-28 2000-03-09 Den Norske Stats Oljeselskap A.S. Method and apparatus for determining the nature of subterranean reservoirs
US6304086B1 (en) * 1999-09-07 2001-10-16 Schlumberger Technology Corporation Method and apparatus for evaluating the resistivity of formations with high dip angles or high-contrast thin layers

Non-Patent Citations (55)

* Cited by examiner, † Cited by third party
Title
"Analysis of multiple source bipole-quadripole resistivity surveys using the apparent resistivity tensor". H.M. Bibby. Geophysics Apr. 1986. pp. 972-983. *
Alumbaugh, D. L. et al, (1997) "Three-Dimensional Massively Parallel Electromagnetic Inversion-II. Analysis of a Crosswell Electromagnetic Experiment", Geophys. J. Int., v. 128, pp. 355-363.
Alumbaugh, D. L. et al, (1998) "3D EM Imaging From a Single Borehole; A Numerical Feasibility Study", SEG Annual Meeting Expanded Abstracts, v. 1, pp. 448-451.
Alumbaugh, D. L. et al, (1998) "Image Appraisal for 2D & 3D Electromagnetic Inversion", SEG Annual Mtg. Expanded Abstracts, v. 1, pp. 456-459.
Bing et al. "Explicit Expressions and Numerical Calculations for the Frechet and second derivatives in 2.5D Helmholtz equation inversion" Jul. 1999. *
Buselli, G. et al, (1996) "Robust Statistical Methods for Reducing Sferics Noise Contaminating Transient Electromagnetic Measurements", Geophysics, v. 61, pp. 1633-1646.
Caldwell et al. "Instantaneous apparent resistivity tensor: a visualization scheme for LOTEM electric field measurements".1998. *
Caldwell, T. G. et al, (1998) "The Instantaneous Apparent Resistivity Tensor: A Visualization Scheme for LOTEM Electric Field Measurements", Geophys. J. Int., v. 135, pp. 817-834.
Chave, A. D., Flosadottir, A. H. and Cox, C. S. (1990) "Some Comments on Seabed Propagation of ULF/ELF Electromagnetic Fields", Radio Science, vol. 25, No. 5, pp. 825-836.
Cheesman, S., Edwards, R. and Chave, A. (1987) "On the Theory of Sea-Floor Conductivity Mapping Using Transient Electromagnetic Systems", Geophysics, vol. 52, No. 2, pp. 204-217.
Constable, S. and Cox, C. S. (1996) "Marine Controlled-Source Electromagnetic Sounding 2. The PEGASUS Experiment", Mar. 1996, J. Geophys. Res., vol. 101, No. B3, pp. 5519-5530.
de Hoop, Adrianus T. (1995) Handbook of Radiation and Scattering of Waves, Academic Press, ISBN 0-12-208655-4, to include Chapters 24, 25 & 26.
Dey, A. et al, (1975) "Electric Field Response of Two-Dimensional Inhomogeneities to Unipolar and Bipolar Electrode Configurations", Geophysics, v. 40, No. 4, pp. 630-640.
Edwards, R. N. (1997) "On the Resource Evaluation of Marine Gas Hydrate Deposits Using Sea-Floor Transient Electric Dipole-Dipole Methods", Geophysics, vol. 62, No. 1, pp. 63-74.
Egbert, G. D., (1997) "Robust Multiple-Station Magnetotelluric Data Processing", Geophys. J. Int., v. 130, pp. 475-496.
Eiedesmo, T. et al, (2002) "Remote Detection of Hydrocarbon Filled Layers Using Marine Controlled Source Electromagnetic Sounding", Z-99, EAGE, 64th Conf., Florence, Italy, May 27-30, 2002, 4 pgs.
Ellingsrud, S. et al., (2002) "Remote Sensing of Hydrocarbon Layers by Seabed Logging (SBL): Results from a Cruise Offshore Angola", The Leading Edge, Oct. 2002, pp. 972-982.
Eswarappa et al., "Mixed boundary semicircular and 120 degrees-sectoral microstrip antennas." IEEE Antennas & Prop. Soc. Int. Symp. pp. 1688-1690. Jun. 1989.
Feynman, Richard P. et al (1964) The Feynman Lectures on Physics, The Electromagnetic Field, Addison-Wesley Pub. Co., Library of Congress No. 63-20717, 3 pgs.
Flosadottir, A. H. and Constable, S. (1996) "Marine Controlled-Source Electroseismic Sounding, Modeling and Experimental Design", Mar. 1996, J. Geophys. Res., v. 101, No. B3, pp. 5507-5517.
Garg, N. and Keller, G. (1984) "Synthetic Electric Sounding Surveys Over Known Oil Fields", Geophysics, Nov. 1984, vol. 49, No. 11, pp. 1959-1967.
Greaves, R. J., Beydoun, W. B. and Spies, B. R. (1991) "New Dimensions in Geophysics for Reservoir Monitoring", SPE Formation Evaluation, Jun. 1991, pp. 141-150.
Gupta, R. N. et al, (1963) "Unipole Method of Electrical Profiling", Geophysics, v. 28, pp. 608-616.
Hoversten, G. M. and Nichols, E. (1992) "Seaborne Electromagnetic Sub-Salt Exploration", Abstracts, American Geophysical Union, p. 313.
Hoversten, G. M. et al, (1998) "Electromagnetics 2: Modeling for Petroleum and Mining Applications", SEG Annual Meeting Expanded Abstracts, v. 1, pp. 425-428.
Hoversten, G. M. et al, (1998) "Marine Magnetotellurics for Petroleum Exploration, Part II: Numerical Analysis of Subsalt Resolution", Geophysics, v. 63, pp. 826-840.
Jupp, D. et al, (1997) "Resolving Anisotrophy in Layered Media by Joint Inversion", Geophys. Prospecting, v. 25, pp. 460-470.
Kaufman, A. A. and Keller, G. V. (1983) "Frequency and Transient Soundings" pub. Elsevier, pp. 285-313.
Kaufman, A. A., Goldman, M., Lee, D. and Keller, G. (1982) "Marine Electromagnetic Prospecting System", Apr. 1982, Geophysics, vol. 47, p. 431.
Leland, "Estimation of boundary value processes applied to shape determination of a circular antenna from observations on the boundary." IEEE Proc. 30th Cond. Decision & Control. Sections 4-5. Dec. 1991.
Li, X. and Pedersen, L. B. (1991) Geophysics, v. 56, No. 9, pp. 1456-1461.
Lockwood et al. "Optimizing the Radiation Pattern of Sparse Periodic Two-Dimensional Arrays". Geoffrey Lockwood and F. Stuart Foster. IEEE. 1996, Abstract only. *
MacGregor, L. M. and Sinha, M. C. (2000) "Use of Marine Controlled-Source Electromagnetic Sounding for Sub-Basalt Exploration", Geophys. Prosp., v. 48, pp. 1091-1106.
MacGregor, L., Constable, S. and Sinha, M. (1998) "The RAMESSES Experiment—III. Controlled-Source Electromagnetic Sounding of the Reykjanes Ridge at 57deg. 45′N", Jun. 1998, Geophys. J. Int., v. 135, pp. 773-789.
MacGregor, L., Constable, S. and Sinha, M. (1998) "The RAMESSES Experiment-III. Controlled-Source Electromagnetic Sounding of the Reykjanes Ridge at 57deg. 45'N", Jun. 1998, Geophys. J. Int., v. 135, pp. 773-789.
Maurer, H. et al, (1998) "Optimized and Robust Experimental Design: A Non-Linear Application to EM Sounding", Geophysics J. Int., v. 132, pp. 458-468.
Mogilatov, V. S. et al, (1996) "A New Method of Geoelectrical Prospecting by Vertical Electric Current Soundings", J. Appl. Geophys., v. 36, pp. 31-41.
Nabighian, M. N. (ed.) (1988) "Electromagnetic Methods in Applied Geophysics", vols. 1 & 2, SEG Investigations in Geophysics No. 3.
Nekut, A. G. and Spies, B. R., (1989) "Petroleum Exploration Using Controlled-Source Electromagnetic Methods", Proceedings IEEE, v. 77, pp. 338-362.
Newman, G. A. et al, (1997) "Three-Dimensional Massively Parallel Electromagnetic Inversion -I. Theory", Report SAND96-0582, Sandia Nat'l Labs, 1996 and Geophys. J. Int., v. 128, pp. 345-354.
Peters, L. J. and Bardeen, J. (1932) "Some Aspects of Electrical Prospecting Applied in Locating Oil Structures", Early Geophysical Papers of the Society of Exploration Geophysicists-vol. II, pp. 145-164.
Rinvert. Schlumberger . Scientific Software Group. 1998. *
Sinha, M. C., Navin, D. A., MacGregor, L. M., Constable, S., Peirce, C., White, A., Heinson, G. and Inglis, M. A. (1997) "Evidence for Accumulated Melt Beneath the Slow-Spreading Mid-Atlantic Ridge", Phil. Trans. R. Soc. Lond., vol. 355, pp. 233-253.
Smith, B. W. and Dzwinel, J. (1984) "WEGA-D System", WEGA-D Geophysical Research Ltd.
Spies, B. R. (1989) "Depth of Investigation in Electromagnetic Sounding Methods", Geophysics v. 54, pp. 872-888.
Strack, K-M (1992) "Exploration with Deep Transient Electromagnetics", Elsevier.
Tarkhov, A. G., (1957) "On Electric Geophysical Exploration Methods of a Pure Anomaly", Bull. Acad. Sci. U.S.S.R., Geophys. Ser., No. 8, pp. 11-22.
U.S. Appl. No. 09/656,191, Sep. 6, 2000, Srnka.
Verma, S. K. et al, (1995) "Focused Resolution of Thin Conducting Layers by Various Dipole EM Systems", Geophysics, v. 60, pp. 381-389.
Wikipedia. Electrostatics, Oct. 28, 2006. *
Zhdanov, M. et al, (1995) "Resistivity Imaging by Time Domain Electromagnetic Migration (TDEMM)", Exploration Geophysics, v. 26, pp. 186-194.
Zhdanov, M. et al, (1996) "Three-Dimensional Quasi-Linear Electromagnetic Inversion", Radio Science, v. 31, pp. 741-754.
Zhdanov, M. et al, (1997) "Time-Domain Electromagnetic Migration in the Solution of Inverse Problems", Geophys. J. Int., v. 131, pp. 293-309.
Zhdanov, M. et al, (1998) "Preconditioned Time Domain Electromagnetic Migration", SEG Annual Meeting Expanded Abstracts, v. 1, pp. 461-468.
Zhdanov, M. S. and Keller, G. V. (1994) "The Geoelectrical Methods in Geophysical Exploration", Elsevier.

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100017156A1 (en) * 2006-08-24 2010-01-21 Anton Ziolkowski Reduction of noise in electrical field measurements
US20080089526A1 (en) * 2006-10-16 2008-04-17 Siemens Audiologische Technik Gmbh Method for estimating an interference field for a coil
WO2015038102A3 (en) * 2013-09-10 2015-11-19 Halliburton Energy Services Inc. Surface calibration of a wellbore resistivity logging tool
GB2533711A (en) * 2013-09-10 2016-06-29 Halliburton Energy Services Inc Surface calibration of a wellbore resistivity logging tool
US10359536B2 (en) 2013-09-10 2019-07-23 Halliburton Energy Services Inc. Surface calibration of a wellbore resistivity logging tool

Also Published As

Publication number Publication date
CA2458170C (en) 2013-11-19
EA007644B1 (ru) 2006-12-29
MY141136A (en) 2010-03-15
BR0211611A (pt) 2004-08-24
AU2002329890B2 (en) 2009-06-04
MXPA04002105A (es) 2004-06-07
CN1276377C (zh) 2006-09-20
US20030050759A1 (en) 2003-03-13
EP1433090A1 (en) 2004-06-30
WO2003025803A1 (en) 2003-03-27
EA200400405A1 (ru) 2004-06-24
CA2458170A1 (en) 2003-03-27
CN1552035A (zh) 2004-12-01
NO20040964L (no) 2004-05-06
US8078439B2 (en) 2011-12-13
MA26217A1 (fr) 2004-07-01
US20100250211A1 (en) 2010-09-30

Similar Documents

Publication Publication Date Title
US7769572B2 (en) Method of imaging subsurface formations using a virtual source array
AU2002329890A1 (en) Method of imaging subsurface formations using a virtual source array
CA2383931C (en) Remote reservoir resistivity mapping
US4617518A (en) Method and apparatus for offshore electromagnetic sounding utilizing wavelength effects to determine optimum source and detector positions
US7203599B1 (en) Method for acquiring transient electromagnetic survey data
US9103927B2 (en) Providing a tow cable having plural electromagnetic receivers and one or more electromagnetic sources
US7894989B2 (en) Method for determining earth vertical electrical anisotropy in marine electromagnetic surveys
CA2650105C (en) Time lapse analysis with electromagnetic data
US8030934B2 (en) Method for hydrocarbon reservoir mapping and apparatus for use when performing the method
US7340348B2 (en) Method for acquiring and interpreting seismoelectric and electroseismic data
Myer et al. Marine CSEM of the Scarborough gas field, Part 1: Experimental design and data uncertainty
RU2397512C1 (ru) Способ электромагнитных геофизических исследований подводных пластов пород
EA018145B1 (ru) Способ определения ориентации датчика для электромагнитной разведки
EA022910B1 (ru) Способ проведения электромагнитной разведки
US20120191352A1 (en) CDP Electromagnetic Marine Data Acquisition and Processing
CN105158808A (zh) 一种浅海瞬变电磁海空探测及其解释方法
US20100235100A1 (en) Method for determining resistivity anisotropy from earth electromagnetic responses
RU2706456C1 (ru) Способ получения результатов измерений в ходе сканирования скрытых или открытых объектов
D. Tereshchenko et al. Radial component of the super-low-frequency magnetic field of a horizontal electric dipole at low conductivity of the ground along the propagation path
Ákos et al. VESZ ADATOK ÉRTELMEZÉSÉNEK ÚJ ELJÁRÁSA: A 1.5-D SZIMULTÁN INVERZIÓS MÓDSZER
NO323241B1 (no) Fremgangsmate for hydrokarbonprospektering i et marint miljo

Legal Events

Date Code Title Description
AS Assignment

Owner name: EXXONMOBIL UPSTREAM RESEARCH COMPANY, TEXAS

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:SRNKA, LEONARD J.;CARAZZONE, JAMES J.;REEL/FRAME:013236/0153

Effective date: 20020821

STCF Information on status: patent grant

Free format text: PATENTED CASE

FPAY Fee payment

Year of fee payment: 4

MAFP Maintenance fee payment

Free format text: PAYMENT OF MAINTENANCE FEE, 8TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1552)

Year of fee payment: 8

MAFP Maintenance fee payment

Free format text: PAYMENT OF MAINTENANCE FEE, 12TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1553); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY

Year of fee payment: 12