Reciprocity in electromagnetics: application to modelling marine magnetometric resistivity data

被引:8
作者
Chen, JP [1 ]
Oldenburg, DW
Haber, E
机构
[1] Univ British Columbia, Geophys Invers Facil, Dept Earth & Ocean Sci, Vancouver, BC V6T 1Z4, Canada
[2] Emory Univ, Dept Math & Comp Sci, Atlanta, GA 30322 USA
基金
加拿大自然科学与工程研究理事会;
关键词
Lorentz reciprocity; adjoint; magnetometric resistivity; electromagnetic induction; forward modelling; marine electromagnetics;
D O I
10.1016/j.pepi.2004.08.015
中图分类号
P3 [地球物理学]; P59 [地球化学];
学科分类号
0708 ; 070902 ;
摘要
The marine magnetometric resistivity (MMR) can be used to obtain resistivity structure beneath the ocean floor. Because of logistical reasons, many more transmitters than receivers are deployed in a survey. This makes it difficult to carry out 3-D forward modelling of magnetometric resistivity data, since each transmitter source requires a separate solution of Maxwell's equations in order to generate the fields. Two methods are presented to overcome this difficulty. The first is based upon the Lorentz reciprocity theorem. With this theorem, the magnetic field at a receiver, generated by a long vertical electrical bipole, is exactly the same as the normalized electromotive force induced in the transmitter wire generated by an artificial magnetic dipole located at the receiver position. The second is the adjoint method in which the magnetic field can be obtained by solving an adjoint equation with an artificial source at each receiver. We show that these two methods are eventually identical: the artificial source in both methods is a steady current in a loop, and the "measurement" is the voltage along the transmitter wire. However the adjoint algorithm is computationally more efficient and we use it in the 3-D marine MMR forward modelling. We verify the code with a synthetic 3-D example. Use of the reciprocity significantly reduces the computational load, making the practical marine MNM problem tractable. (c) 2004 Elsevier B.V. All rights reserved.
引用
收藏
页码:45 / 61
页数:17
相关论文
共 24 条
[1]   Imaging and time reversal in random media [J].
Borcea, L ;
Papanicolaou, G ;
Tsogka, C ;
Berryman, J .
INVERSE PROBLEMS, 2002, 18 (05) :1247-1279
[2]  
CHAVE AD, 1991, ELECTROMAGNETIC METH, V2, P931
[3]   Three-dimensional numerical modelling and inversion of magnetometric resistivity data [J].
Chen, JP ;
Haber, E ;
Oldenburg, DW .
GEOPHYSICAL JOURNAL INTERNATIONAL, 2002, 149 (03) :679-697
[4]   A nonlinear inversion method for 3D electromagnetic imaging using adjoint fields [J].
Dorn, O ;
Bertete-Aguirre, H ;
Berryman, JG ;
Papanicolaou, GC .
INVERSE PROBLEMS, 1999, 15 (06) :1523-1558
[5]   ON MEASURING THE ELECTRICAL-CONDUCTIVITY OF THE OCEANIC-CRUST BY A MODIFIED MAGNETOMETRIC RESISTIVITY METHOD [J].
EDWARDS, RN ;
LAW, LK ;
DELAURIER, JM .
JOURNAL OF GEOPHYSICAL RESEARCH, 1981, 86 (NB12) :1609-1615
[6]  
EDWARDS RN, 1991, ELECTROMAGNETIC METH, V2, P47
[7]   Hydrothermal circulation at the Cleft-Vance overlapping spreading center: Results of a magnetometric resistivity survey [J].
Evans, RL ;
Webb, SC ;
Jegen, M ;
Sananikone, K .
JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH, 1998, 103 (B6) :12321-12338
[8]   Crustal resistivity structure at 9°50′N on the East Pacific Rise:: Results of an electromagnetic survey -: art. no. 1082 [J].
Evans, RL ;
Webb, SC .
GEOPHYSICAL RESEARCH LETTERS, 2002, 29 (06)
[9]   Fast simulation of 3D electromagnetic problems using potentials [J].
Haber, E ;
Ascher, UM ;
Aruliah, DA ;
Oldenburg, DW .
JOURNAL OF COMPUTATIONAL PHYSICS, 2000, 163 (01) :150-171
[10]  
HABER E, 1943, SIAM J SCI COMPUT, V22