Rapid inversion of electromagnetic reconnaissance survey data
원문보기
IPC분류정보
국가/구분
United States(US) Patent
등록
국제특허분류(IPC7판)
G01V-001/40
G01V-003/00
G01R-023/16
G06G-007/58
출원번호
US-0375188
(2007-07-30)
등록번호
US-8538699
(2013-09-17)
국제출원번호
PCT/US2007/016986
(2007-07-30)
§371/§102 date
20090126
(20090126)
국제공개번호
WO2008/033184
(2008-03-20)
발명자
/ 주소
Jing, Charlie
Willen, Dennis E.
Carazzone, James J.
Pavlov, Dmitriy A.
출원인 / 주소
ExxonMobil Upstream Research Company
대리인 / 주소
ExxonMobil Upstream Research Company Law Dept.
인용정보
피인용 횟수 :
1인용 특허 :
126
초록▼
Method for rapid inversion of data from a controlled-source electromagnetic survey of a subterranean region. Selected (51) common-receiver or common-source gathers of the data are reformed into composite gathers (52) by summing their data. Each composite gather is forward modeled (in the inversion p
Method for rapid inversion of data from a controlled-source electromagnetic survey of a subterranean region. Selected (51) common-receiver or common-source gathers of the data are reformed into composite gathers (52) by summing their data. Each composite gather is forward modeled (in the inversion process) with multiple active source locations (53). Computer time is reduced in proportion to the ratio of the total number of composite gathers to the total number of original common-receiver or common-source gathers. The data may be phase encoded to prevent data cancellation. Methods for mitigating loss of far offset information by data overlap in the summing process are disclosed.
대표청구항▼
1. A computer implemented method for rapid inversion of electromagnetic data from a controlled source electromagnetic survey of a subterranean region, comprising: (a) summing measured electromagnetic survey data, using a computer, to form one or more composite gathers, each composite gather being fo
1. A computer implemented method for rapid inversion of electromagnetic data from a controlled source electromagnetic survey of a subterranean region, comprising: (a) summing measured electromagnetic survey data, using a computer, to form one or more composite gathers, each composite gather being formed from at least two selected ordinary gathers, said ordinary gathers being either common-receiver or common-source gathers; wherein loss of sensitivity to a subsurface resistivity anomaly due to data overlap or data cancellation, or both, in the one or more composite gathers is mitigated; wherein said loss of sensitivity is defined relative to sequential inversion, meaning inverting one ordinary gather at a time;wherein data overlap is mitigated by one or more of the following measures when forming the one or more composite gathers: (i) summing only a subset of the ordinary gathers in the survey, depending upon one or more of depth range of interest, earth conductivity, frequency of the data, and other parameters of the survey; (ii) preserving far-offset data by summing only selected components of the electromagnetic field from neighboring receivers; andwherein data cancellation is mitigated by phase encoding the ordinary gathers in forming the one or more composite gathers; and(b) inverting the composite gathers, using a computer, for resistivity or another subsurface electrical-property parameter. 2. The method of claim 1, further comprising inverting at least one ordinary gather with the composite gathers. 3. The method of claim 1 wherein the ordinary gathers are common-receiver gathers and summing measured electromagnetic survey data to form a composite gather comprises: (i) selecting a set of transmitter (source) locations for forming the composite gather;(ii) selecting at least one electromagnetic field data component for each selected ordinary gather; and(iii) at every selected transmitter location, summing all selected data components from all selected ordinary gathers , thereby forming the composite gather. 4. The method of claim 1, wherein the ordinary gathers are common-source gathers and summing measured electromagnetic survey data to form a composite gather comprises: (i) selecting a set of receiver locations for forming the composite gather;(ii) selecting at least one electromagnetic field data component for each selected ordinary gather; and(iii) at every selected receiver location, summing all selected data components from all selected ordinary gathers , thereby forming the composite gather. 5. The method of claim 3 or claim 4, wherein the selection of ordinary gathers and electromagnetic field components for the composite gather is governed by consideration of factors including (i) reducing degree of data overlap between neighboring ordinary gathers to preserve contributions from longer offsets and (ii) enhancing inversion efficiency by increasing the number of ordinary gathers in one composite gather. 6. The method of claim 5, wherein ordinary gathers are selected to have receiver spacing in the case of common-receiver gathers or source spacing in the case of common-source gathers great enough to reduce data overlap to a pre-determined level. 7. The method of claim 5, wherein at least two different electromagnetic field components are selected alternately from neighboring ordinary gathers to reduce data overlap. 8. The method of claim 5 wherein an additional factor considered is a pre-selected depth range of interest in the subsurface region. 9. The method of claim 5, wherein the method is performed in frequency domain, and the selection of ordinary gathers to reduce data overlap is frequency dependent. 10. The method of claim 1, wherein said inverting the composite gathers comprises: (i) constructing an initial resistivity model;(ii)forward modeling at least one electromagnetic field component, wherein Maxwell's electromagnetic field equations are numerically solved assuming the resistivity model and with all sources associated with the composite gather assumed to be simultaneously active, thus generating simulated data;(iii) calculating data mismatch between the simulated data and the measured survey data for the corresponding composite gathers; and(iv) adjusting the resistivity model to reduce the data mismatch repeating steps (ii) to (iv) until a preset tolerance level or other stopping point is met. 11. The method of claim 10, wherein said adjusting the resistivity model is done by minimizing an objective function, said objective function being selected so as to be a measure of the data mismatch. 12. The method of claim 1, wherein the summing and inverting steps are performed in time domain. 13. The method of claim 1, wherein the summing and inverting steps are performed in frequency domain, wherein time dependence of measured survey data are transformed to frequency dependence. 14. The method of claim 1, wherein common-receiver gathers are used to form the composite gathers if the survey has more source positions (as determined by data binning strategy in the case of a moving source) than receiver positions, and otherwise common-source gathers are used. 15. The method of claim 1, wherein the survey has more source positions (as determined by data binning strategy in the case of a moving source) than receiver positions, and common-receiver gathers are used to form the composite gathers, and wherein source and receiver roles are exchanged (reciprocity invoked) for forward-modeling computations in the inverting step. 16. The method of claim 1, wherein phase encoding is used in forming composite gathers and in forward-modeling computations associated with the inverting step to avoid signal cancellation from summing two electromagnetic field values of opposite phase. 17. The method of claim 1, wherein the measured survey data to be summed includes magnetic field data. 18. The method of claim 17, wherein reciprocity of source and receiver is invoked in forward modeling calculations in the inverting step, and wherein computational magnetic sources (corresponding to actual magnetic receivers) are each modeled as an effective single-turn circular antenna with core material of the same magnetic permeability as background material at the actual receiver's location. 19. The method of claim 10, wherein the forward modeling uses a method selected from a group consisting of (a) finite difference method; (b) finite element method, and (c) integral equation method. 20. A method for producing hydrocarbons from a subsurface region, comprising: (a) performing a controlled-source electromagnetic survey of the subsurface region;(b) obtaining a resistivity model of the subsurface region produced by inverting electromagnetic data measured in the survey, said inversion having been performed by using a method of claim 1;(c) drilling a well into a zone indicated as an anomaly in the resistivity model; and(d) producing hydrocarbons from the well. 21. The method of claim 16, wherein the phase encoding is performed using encoding functions that are designed to reduce the signal cancellation.
Hayashi,Takeshi; Sonobe,Kouichi; Yoshimura,Takahiro; Suzuki,Masakazu, Automatic exposure control method of image and automatic exposure control system using the method.
Schoen, Juergen S.; Fanini, Otto N.; Georgi, Daniel, Combined characterization and inversion of reservoir parameters from nuclear, NMR and resistivity measurements.
Alumbaugh,David; Chen,Jiuping; Morrison,Huntly F., Computing values for surveying a subterranean structure based on measurements according to different electromagnetic survey techniques.
Tasci M. Tahsin (Denver CO) Jordan John M. (Golden CO), Electromagnetic imaging device and method for delineating anomalous resistivity patterns associated with oil and gas tra.
Ramamoorthy, Raghu; Srivastava, Ashok; Essawi, Amr; Gibson, William Dean, Generating and displaying a virtual core and a virtual plug associated with a selected piece of the virtual core.
King Thomas C. (Rte. 2 Box 59 San Saba TX 76877) Arnason Byron T. (4104 Honeycomb Rock Cl. Austin TX 78731), Geophysical prospecting with collimated magnetotelluric fields.
Rosenquist Mark L. ; Dwan Fa S. ; Bryndzia L. Taras ; Pelletier John H., Method and apparatus for detecting shallow water flow sands using converted shear waves.
Minerbo Gerald N. ; Rosthal Richard A., Method and apparatus for evaluating the resistivity of formations with high dip angles or high-contrast thin layers.
Srnka Leonard J. (Houston TX), Method and apparatus for offshore electromagnetic sounding utilizing wavelength effects to determine optimum source and.
Freedman Robert (Houston TX) Minerbo Gerald N. (Missouri City TX), Method and apparatus for producing a more accurate resistivity log from data recorded by an induction sonde in a borehol.
Omeragic, Dzevat; Habashy, Tarek M.; Esmersoy, Cengiz, Method for calculating a distance between a well logging instrument and a formation boundary by inversion processing measurements from the logging instrument.
Neff Dennis B. ; Runnestrand Scott A. ; Butler Edgar L. ; Vienot Michael E., Method for determining distribution of reservoir permeability, porosity and pseudo relative permeability.
Gregory Joseph Jorgensen ; Jerry Lee Kisabeth ; Alan Royce Huffman ; John B. Sinton ; David W. Bell, Method for gravity and magnetic data inversion using vector and tensor data with seismic imaging and geopressure prediction for oil, gas and mineral exploration and production.
Gregory Joseph Jorgensen ; Jerry Lee Kisabeth ; Alan Royce Huffman ; John B. Sinton ; David W. Bell, Method for integrating gravity and magnetic inversion data with model based seismic data for oil, gas and mineral exploration and production.
Gregory Joseph Jorgensen ; Jerry Lee Kisabeth ; Alan Royce Huffman, Method for integrating gravity and magnetic inversion with geopressure prediction for oil, gas and mineral exploration and production.
Guerillot Dominique,FRX ; Roggero Frederic,FRX, Method for predicting, by means of an inversion technique, the evolution of the production of an underground reservoir.
Safinya Kambiz A. (Garches CT FRX) Habashy Tarek M. (Danbury CT) Beren Jeffrey A. (Westport CT), Method of evaluating a geological formation using a logging tool including slot antenna having two nonparallel elements.
Laws,Robert; Morice,Stephen Patrick, Method of seismic surveying, a marine vibrator arrangement, and a method of calculating the depths of seismic sources.
Edington B. Leon (Houston TX) Khan Tawassul A. (Houston TX), Method of source coding and harmonic cancellation for vibrational geophysical survey sources.
Druskin Vladimir L'vovich ; Lee Ping ; Knizhnerman Leonid Aronovich,RUX, Method, apparatus, and article of manufacture for solving 3D Maxwell equations in inductive logging applications.
Kosloff Dan D. (Ramat-Aviv ILX) Meshbey Oleg V. (Lod ILX) Koren Zvi (Ramat-Gan ILX) Litvin Alexander (London GBX), Migration process using a model based aperture technique.
Chavarria,Juan Andres; Routh,Partha Sarathi; Kisabeth,Jerry Lee; Jorgensen,Gregory Joseph, Modeling gravity and tensor gravity data using poisson's equation for airborne, surface and borehole applications.
Routh, Partha Sarathi; Kisabeth, Jerry Lee; Jorgensen, Gregory Joseph, Nonlinear constrained inversion method to determine base of salt interface from gravity and gravity tensor data.
Fanini, Otto N.; Withjack, Eric; Treitel, Sven, Reservoir monitoring in a laminated reservoir using 4-D time lapse data and multicomponent induction data.
Thomsen,Leon; Allegar,Norman C.; Dellinger,Joseph A.; Jilek,Petr; Johnson,Daniel; Xia,Ganyuan, System and method for using time-distance characteristics in acquisition, processing, and imaging of t-CSEM data.
Richard A. Mollison ; Otto N. Fanini ; Berthold Kriegshauser ; Milomir Pavlovic, Water saturation and sand fraction determination from borehole resistivity imaging tool, transverse induction logging and a tensorial dual water saturation model.
Mollison, Richard A.; Schoen, Juergen S.; Fanini, Otto N.; Kriegshauser, Berthold; Pavlovic, Milomir, Water saturation and sand fraction determination from borehole resistivity imaging tool, transverse induction logging and a tensorial water saturation model.
Richard A. Mollison ; Juergen S. Schoen AT; Otto N. Fanini ; Berthold Kriegshauser ; Milomir Pavlovic, Water saturation and sand fraction determination from borehole resistivity imaging tool, transverse induction logging and a tensorial water saturation model.
Clark Brian (Redding CT) Wu Peter T. (Sugar Land TX) Grijalva Victor E. (Greenwich CT), Well logging apparatus and method with borehole compensation including multiple transmitting antennas asymmetrically dis.
※ AI-Helper는 부적절한 답변을 할 수 있습니다.