2012 SEG SEG Las Vegas 2012 Annual Meeting Page 1
|
|
- Doris Dixon
- 5 years ago
- Views:
Transcription
1 Frequency-domain EM modeling of 3D anisotropic magnetic permeability and analytical analysis Jiuping Chen, Michael Tompkins, Ping Zhang, Michael Wilt, Schlumberger-EMI, and Randall Mackie, Formerly Schlumberger-EMI, presently Fugro-EM SUMMARY 3D NUMERICAL MODELING OF µ In this paper we present the modification and verification of a 3D frequency-domain electromagnetic modeling algorithm for dealing with heterogeneous, anisotropic magnetic permeability. To test the modifications, the code is verified against a 1D layered earth response, and 3D axially-symmetric doughnut model. In addition, to better understand the magnetic field responses to both conductivity and permeability, analytical analysis is performed and demonstrates that in a low induction number ( k r 1) regime the magnetic current or magnetization may be more important thant the induced electric current in determining the EM response to a medium. INTRODUCTION In applied EM geophysics, estimation of electrical resistivity has been widely used for a variety of applications such as mineral prospecting, well-logging, deep earth structure exploration and environmental monitoring. On the other hand, magnetic permeability or susceptibility, one of the most fundamental physical properties of rock, has usually been limited to studies in static magnetic exploration for massive ore deposits and UneXplored Ordnance (UXO) classification and identification (i.e., Huang and Won, 2003). In addition to limited application, few 3D numerical algorithms are available for use in simulating effects due to magnetic permeability. Most well-known and widely used solutions are based on the analytic expressions of a conductive and permeable sphere in air, excited by a magnetic dipole source (e.g., Grant and West, 1965). In a recently published paper by Mukherjee and Everett (2011), they tried to verify a 3D finite-element solution dealing with both conductivity and permeability with the Grant and West (1965) solution. Other algorithms, as far as we know, such as Wang and Hohmann (1993), Alumbaugh et al. (1996), Haber et al. (2000), and Lovell (1990) are capable of simulating both conductivity and permeability. However no detailed studies exist from these authors on the effect of permeability. In this study, we aim to extend an existing 3D finite-difference EM forward modeling algorithm to cover heterogeneous, anisotropic magnetic permeability. In addition, we carry out an analytical analysis by using the integral representation of the EM response for a spherical model to help understand the physics of the secondary field. Our goal is to get a sense what role the magnetic current or magnetization plays relative to the induced electric current in generating secondary magnetic fields for low induction numbers, which has implications for how interpretations can be made. Modification of the 3D algorithm Here, we modify the 3D numerical EM code as presented in a series of journal papers (Mackie et al.,1994; Madden and Mackie, 2002). In this algorithm, 3D earth models are divided into regular prism cells. The forward modeling algorithm uses a finite-difference method to solve for the magnetic fields along the edges of the cells; the electric fields can be obtained by interpolation. The resistivity at each cell can be diagonally anisotropic, i.e., ρ xx,ρ yy,and ρ zz. The magnetic permeability is assumed to be constant everywhere, with a freespace value of µ 0 = 4π 10 7 H/m. The code is parallelized using the open source PETSc (Portable, Extendable Toolkit for Scientific Computation) platform, a suite of data structures and routines for the scalable (parallel) solution of scientific applications modeled by partial differential equations (Balay et al., 2011). A Krylov-based linear iterative solver with a pre-conditioner matrix is used to solve the assembled linear equations. To deal with magnetic permeability, in the region where there is no source term, we have to go back to the integral forms of Maxwell s equations, which (following Mackie et al., 1994) are H dl = σe ds, (1) E dl = iωµh ds, (2) where σ and µ are the electric conductivity and magnetic permeability, and both can be tensor quantities. Since here we are concerned with the magnetic permeability, let us only focus on equation 2. With this formulation and the discretization scheme shown in Figure 1, for example, the H y -component of equation 2 can be written as {[E x (i, j,k) E x (i, j,k 1)] [E z (i, j,k) E z (i 1, j,k)]}l = iω < µ yy > H y (i, j,k)l 2, (3) where L is the grid spacing in the transformed system, and the average magnetic permeability < µ yy > is defined as: < µ yy >= 1 4 [µ yy(i 1, j,k 1) +µ yy (i, j,k 1) + µ yy (i 1, j,k) + µ yy (i, j,k)]. (4) Similarly we can work on H x and H z -component of equation 2, and then combine the discretized forms of electric field in equation 1, grouping these second-order equations together in the form: M xx N xy N xz H x Ax = N yx M yy N yz = b (5) N zx N zy M zz H y H z SEG Las Vegas 2012 Annual Meeting Page 1
2 Figure 1: The finite-difference geometry for the integral forms of Maxwell s equations. where b contains the terms associated with the known boundary values and source field,and M i j and N i j are the sub-stiff coefficient matrices. Accordingly, the original coefficient matrix A, the preconditioning matrix P, the right-hand side vector b in equation 5, as well as the new static correction term ρ µh need to be changed by incorporating the diagnostically anisotropic magnetic permeability. EM reciprocity helps to reduce computational time in the case that the number of transmitters is much greater than the number of receivers, such as in many marine CSEM applications. In addition, generally the reciprocal relationship holds for magnetic H field if µ is the same as µ 0. However, when the magnetic permeability varies, the reciprocity is valid only for magnetic induction and B field. For an i-oriented magnetic dipole source, and measuring j-oriented magnetic induction component, where i and j can be pointing x,y,or z, according to Chen et al.(2005), we have the following reciprocal expression B i j (r tx,r rx ) = B ji (r rx,r tx ), (6) where B ji (r rx,r tx ) is the i-component of B field observed at the artificial location r tx, generated by a j-oriented artificial magnetic dipole source at location r rx. This means µ j j (r rx )H i j (r tx,r rx ) = µ ii H ji (r rx,r tx ). (7) By using the reciprocity, the solution obtained by solving equation 5 is H ji, therefore, we perform the following step Verification of the algorithm H i j (r tx,r rx ) = µ ii µ j j H ji (r rx,r tx ). (8) We present two test examples to verify the modifications made to include magnetic permeability. These two examples are (1) a layered 1D µ medium; and (2) a 3D doughnut µ model (axial symmetry) in a whole space. Figure 2: Comparison of CSEM3Dmu vs semi-analytical solutions in a 1D layered µ medium. As shown in the inset in Figure 2, the resistivity and magnetic permeability of the uniform whole space are 10 Ω.m and µ 0. The 1D layered model contains a µ layer with a thickness of 75 m, and a relative permeability of 1.0 and 1.2. We use a cross-well configuration for the test. The separation between the vertical Tx well (vertical magnetic dipole, VMD) and the Rx well is 200m, and the applied frequency is 500 Hz. We show the computed magnetic field H z plot along the Rx profile, generated by a VMD at a constant depth 1100 m for both the 1D, and 3D codes. The semi-analytical solutions were obtained by running a numerical 1D code (Tompkins, 2003), which are displayed in black dashed lines (free space µ 0 ), and green solid line (µ r = 1.2), respectively. The numerical responses obtained from running the modified code CSEM3Dmu are plotted in blue dashed lines (free space) and red solid line (µ r = 1.2). It can be seen that for the free space, the black and blue lines match well; for the µ r model, the green and red lines are in a fairly good agreement. The relative difference between them is less than 2%, which can be attributed to the grid size used in the 3d modeling. The sharp drop of the vertical H z component at depth from 1075 to 1100m is primarily caused by the µ r change, because we know that the normal magnetic induction B z is continuous cross the boundary layer, while H z is not (i.e., it is scaled by µ r ). The second test example is a 3D doughnut µ model imbedded in a uniform wholespace of resistivity (10 Ω.m) and magnetic permeability (µ 0 ). The 3D doughnut is displayed in Figure 3, with an inner and outer radius of 25 and 100 m, respectively. The thickness of the doughnut is 75 m, and µ r = 1.2. Again, the Tx (VMD) is located at depth 1100 m (inside the doughnut), and the source frequency is 500 Hz. The comparison was run against the code CWNLAT (Lovell, 1990). CWNLAT is a finite element program, and allows electromagnetic simulation in media with variable resistivity, magnetic permeability, and dielectric properties. Assuming an axiallysymmetric TM excitation in an azimuthally symmetric (2 dimensions) cylindrical geometry. The doughnut model is a perfect example for this code. SEG Las Vegas 2012 Annual Meeting Page 2
3 Figure 3: Comparison of CSEM3Dmu vs CWNLAT solutions for a 3D doughnut µ model in a uniform whole space. As shown in Figure 3, the computed H z field from CSEM3Dmu (red solid line) is in good agreement with the CWNLAT responses which are in green. In order to improve the numerical accuracy of CSEM3Dmu, we used denser grids around both Tx and Rx locations. This was required, because the Rx is treated as an artificial source in the code, and each artificial source is distributed at the 8 nodes surrounding the source point by volume weighting. Similarly, the field computed at the artificial Rx (actually the real Tx location) is also interpolated by volume weighting at the 8 nodes around that point. ANALYTICAL ANALYSIS OF MAGNETIC RESPONSES Integral representation of magnetic responses The general EM induction phenomenon for a conductive target buried in an extremely resistive background is well known, and has been clearly understood and conceptually illustrated in Figure 4 (adapted from Grant and West, 1965). When the resistivity contrast between the target and surrounding host is limited, besides the vortex current induced in the target, there is also channeled current or galvanic current generated around the body (McNeill et al., 1984). If the target is also magnetically permeable, more physical phenomena are present. We illuminate these phenomena through analysis of analytical solutions in the low frequency case. For a general medium (with conductivity σ 0 and magnetic permeability µ 0 ) containing both a conductive and permeable object (σ and µ), and excited by an external source, the magnetic field at an observation point consists of incident field H inc and scattered or secondary field H s and the secondary field H s can be represented as (de Hoop, 1995) H s (r) = g(r r )[σ(r ) σ 0 (r )]E(r )dv(r ) v +[k0 2 + ] g(r r )[ µ 1]H(r )dv(r ) v µ 0 = H se + H sm, (9) Figure 4: A carton showing concept EM induction for both a conductive and permeable body. Adapted from Grant and West(1965) where k 0 = iωµ 0 σ 0, is the wave number, and the homogenerous scalar Green s function is given by g(r) = e ik 0 r 4π r. (10) In equation 9, the first term is the magnetic field generated by the excessive electrical current density [σ σ 0 ]E(r ), and thus we denote it by H se. The second term is the magnetic field generated by the magnetic current (Wait, 1953) or magnetization [ µ µ 0 1]H(r ), and therefore can be distinguished by H sm. We would like to analyze how these two terms behave in a typical low frequency regime. Analytical solutions To simplify our analysis, let us take the Tx-Rx configuration as shown in Figure 5 as an example. The Tx is a VMD with a magnetic moment m, and is put on the same vertical plane xoz as Rxs, which only measure the vertical field H z. The sphere is located at the origin of the system, with a radius a, conductivity σ, and permeability µ, Assume a is much smaller than the distance r 0. Then the electric field inside the sphere can be approximated by the y component of the incident electric field, i.e., and for a VMD source, we know that E(r i ) Eyŷ, (11) Ey i = iωµ 0m 4πr0 2 (ik 0 r 0 + 1)e ik 0r 0. (12) Under the assumptions that both induction numbers k 0 r 0, k 0 r 1, and the vertical distance z is much less than r, without dwelling on detailed derivations, the vertical component Hz se of the secondary magnetic field due to the electric current can be simplified as H se z iωµ 0m 4πr 2 0 }{{} [σ σ 0 ] πa3 }{{}} 4πr {{ 2 } (13) Apparently, the first bracket denotes the incident E field at the sphere. Then when multiplied by the second bracket, it gives SEG Las Vegas 2012 Annual Meeting Page 3
4 Figure 5: A schematic spherical model for understanding the magnetic response from both conductive and magnetically permeable target. rise to a current segment with a unit of A.m. This current segment produces the magnetic field at the Rx with the geometric term in the third bracket. Following similar algebraic operations, the secondary vertical magnetic field Hz sm from the magnetic current can be approximated as Hz sm ( m 4πr0 3 ) [µ r 1] 4 3 πa3 ( 1 }{{}}{{} 4πr 3 ) }{{} (14) Again the physical meaning can be inferred from the three brackets shown in equation 14. Here we see that µ r 1 = χ is the suscepibility of the sphere, and then when excited by the external magnetic field, it produces the magnetization, which in turn generates another secondary magnetic field as well. To quantify the above analysis, a spherical model was examined, and the computed results are shown in Figure 6. In this model, the resistivities for the background and the sphere are 100 and 10 Ω. m, respectively. The relative permeability of the sphere is 2. The radius is 50 m. The distance between the VMD to the center of the sphere is 200 m, and the horizontal distance from the origin to the Rx well is also 200m. The depth of the Rxs range from -100 to 100 m. At 10 Hz, the assumptions made in the previous section are roughly satisfied. Figure 6a shows the primary, secondary, and the total H z field. In this case, the secondary field is about 3% of the primary. Therefore, both primary and the total field are almost completely overlapping at the top. Even though the secondary field is relatively small, Figure 6b shows that the magnetization field Hz sm is about 3 times as large as the magnetic field Hz se generated by the induced current. The point we would like to make here is that at low frequency and in a relatively resistive host, magnetization, rather than induced current, appears to dominate the observed secondary field as opposed to what is observed in conventional cross-well EM surveys. Figure 6: An example of the approximate analytical analysis on the sphere. (a) the primary, secondary, and total H z. (b) the ratio of Hsm z H in this model. z se CONCLUSION We have made modifications to an existing 3D EM algorithm to accommodate the need for simulating the effect of magnetic permeability heterogeneities including anisotropy. Two test examples have been presented, and demonstrate that the modified code compares well with semi-analytical and alternative numerical solutions, providing a useful tool for the simulation and monitoring of magnetically permeable formations. We also have carried out analysis of an analytical solution for magnetic field responses using an integral representation for a conductive and permeable sphere in a dipole field. The approximate analysis shows that the secondary magnetic field can be attributed to two parts: one is generated by the induced electric current, and the other by the magnetic current or magnetization. At a low induction number regime ( k r 1), the magnetization part may dominate in the observed anomaly over the induced electric current. This physical understanding will help us to better apply this technique to different applications, for example, hydraulic fracturing by injection of magnetically permeable particles into formations. SEG Las Vegas 2012 Annual Meeting Page 4
5 EDITED REFERENCES Note: This reference list is a copy-edited version of the reference list submitted by the author. Reference lists for the 2012 SEG Technical Program Expanded Abstracts have been copy edited so that references provided with the online metadata for each paper will achieve a high degree of linking to cited sources that appear on the Web. REFERENCES Alumbaugh, D. L., G. A. Newman, L. Prevost, and J. N. Shadid, 1996, Three-dimensional, wideband electromagnetic modeling on massively parallel computers: Radio Science, 31, Balay, S., J. Brown, K. Buschelman, V. Eijkhout, W. D. Gropp, D. Kaushik, M. G. Knepley, L. C. Mclnnes, B. F. Smith, and H. Zhang, 2011, PETSc Users Mannual, ANL-95/11-Revision 3.2, Argonne National Laboratory. Chen, J., D. W. Oldenburg, and E. Haber, 2005, Reciprocity in electromagnetics: Applications to marine MMR: Physics of the Earth & Planetary Interiors, 150, Grant, F. S., and G. F. West, 1965, Interpretation theory in applied geophysics: McGraw -Hill. Haber, E., U. M. Ascher, D. A. Aruliah, and D. W. Oldenburg, 2000, Fast simulation of 3D electromagnetic problems using potentials: Journal of Computational Phys ics, 163, Huang, H., and I. J. Won, 2003, Characterization of UXO -like targets using broadband electromagnetic induction sensors: IEEE Transactions on Geosciences and Remote Sensing, 41, Lovell, J., 1990, Modeling frequency effects on laterologs: Reserch note, Schlumberger, 10 July Mackie, R. L., J. T. Smith, and T. R. Madden, 1994, Three-dimensional electromagnetic modeling using finite difference equations: The magnetotelluric example: Radio Science, 29, no. 4, Madden, T. M., and R. L. Mackie, 2002, Three-dimensional magnetotelluric modeling and inversion: Proceeding of the IEEE, 77, no. 2, McNeill, J. D., R. N. Edwards, and G. M. Levy, 1984, Approximate calculations of the transient EM response from buried conduc tors in a conductive half -space: Geophysics, 49, Mukherjee, S., and M. E. Everett, 2011, 3D controlled-source EM edge-based finite element modeling of conductive and permeable heterogeneities: Geophysics, 76, no. 4, F215 F226. Tompkins, M., 2003, Quantitative analysis of mulicomponent borehole electromagnetic induction responses using anisotropic forward modeling and inve rsion: P.h.D. thesis, University of Wisconsin. Wait, J., 1953, A conducting permeable sphere in the presence of a coil carry ing an oscillating current: Canadian Journal of Physics, 31, Wang, T., and G. W. Hohmann, 1993, A finite -difference time domain solution for 3D EM modeling: Geophysics, 58, SEG Las Vegas 2012 Annual Meeting Page 5
Downloaded 01/09/13 to Redistribution subject to SEG license or copyright; see Terms of Use at
3D modeling of IP effects on electromagnetic data in the time domain David Marchant, Eldad Haber, Laurens Beran and Douglas W. Oldenburg, University of British Columbia Downloaded /9/3 to 37.82.7.. Redistribution
More informationInductive source induced polarization David Marchant, Eldad Haber and Douglas W. Oldenburg, University of British Columbia
David Marchant, Eldad Haber and Douglas W. Oldenburg, University of British Columbia Downloaded /9/3 to 37.8.7.. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/
More informationModelling IP effects in airborne time domain electromagnetics
Modelling IP effects in airborne time domain electromagnetics Dave Marchant* Seogi Kang Mike McMillian Eldad Haber Computational Geosciences University of British Columbia Computational Geosciences University
More informationMagnetotelluric (MT) Method
Magnetotelluric (MT) Method Dr. Hendra Grandis Graduate Program in Applied Geophysics Faculty of Mining and Petroleum Engineering ITB Geophysical Methods Techniques applying physical laws (or theory) to
More informationAchieving depth resolution with gradient array survey data through transient electromagnetic inversion
Achieving depth resolution with gradient array survey data through transient electromagnetic inversion Downloaded /1/17 to 128.189.118.. Redistribution subject to SEG license or copyright; see Terms of
More informationPhysics-Based Characterization of UXO from Multi-Component TEM Data. Abstract
Physics-Based Characterization of UXO from Multi-Component TEM Data Scott C. MacInnes, Donald D. Snyder, and Kenneth L. Zonge Zonge Engineering & Research Organization 3322 E Fort Lowell Rd Tucson, AZ,
More informationA Brief Introduction to Magnetotellurics and Controlled Source Electromagnetic Methods
A Brief Introduction to Magnetotellurics and Controlled Source Electromagnetic Methods Frank Morrison U.C. Berkeley With the help of: David Alumbaugh Erika Gasperikova Mike Hoversten Andrea Zirilli A few
More informationWhere does the proppant go? Examining the application of electromagnetic methods for hydraulic fracture characterization
Where does the proppant go? Examining the application of electromagnetic methods for hydraulic fracture characterization Lindsey J. Heagy and Douglas W. Oldenburg Geophysical Inversion Facility, University
More informationSEG Las Vegas 2008 Annual Meeting 677
in the context of the generalized effective-medium theory Michael S Zhdanov, Alexander Gribenko, Vladimir Burtman, Consortium for Electromagnetic Modeling and Inversion, University of Utah, and Vladimir
More informationDownloaded 08/29/13 to Redistribution subject to SEG license or copyright; see Terms of Use at
New approach to 3D inversion of MCSEM and MMT data using multinary model transform Alexander V. Gribenko and Michael S. Zhdanov, University of Utah and TechnoImaging SUMMARY Marine controlled-source electromagnetic
More informationSimulation of 3D DC Borehole Resistivity Measurements with a Goal- Oriented hp Finite-Element Method. Part I: Laterolog and LWD
Journal of the Serbian Society for Computational Mechanics / Vol. 1 / No. 1, 2007 / pp. 62-73 Simulation of 3D DC Borehole Resistivity Measurements with a Goal- Oriented hp Finite-Element Method. Part
More informationTu Efficient 3D MT Inversion Using Finite-difference Time-domain Modelling
Tu 11 16 Efficient 3D MT Inversion Using Finite-difference Time-domain Modelling S. de la Kethulle de Ryhove* (Electromagnetic Geoservices ASA), J.P. Morten (Electromagnetic Geoservices ASA) & K. Kumar
More informationFinite Element Method (FEM)
Finite Element Method (FEM) The finite element method (FEM) is the oldest numerical technique applied to engineering problems. FEM itself is not rigorous, but when combined with integral equation techniques
More informationSIMULATION AND INTERPRETATION OF BOREHOLE GEOPHYSICAL MEASUREMENTS USING hp FINTE ELEMENTS
SIMULATION AND INTERPRETATION OF BOREHOLE GEOPHYSICAL MEASUREMENTS USING hp FINTE ELEMENTS hp-fem team: D. Pardo, M. J. Nam, L. Demkowicz, C. Torres-Verdín, V. M. Calo, M. Paszynski, and P. J. Matuszyk
More informationJoint inversion of borehole electromagnetic and sonic measurements G. Gao, A. Abubakar, T. M. Habashy, Schlumberger-Doll Research
Joint inversion of borehole electromagnetic and sonic measurements G. Gao, A. Abubakar, T. M. Habashy, Schlumberger-Doll Research SUMMARY New-generation electromagnetic (EM) and sonic logging tools are
More informationMain Menu. Douglas Oldenburg University of British Columbia Vancouver, BC, Canada
Using ERA low frequency E-field profiling and UBC 3D frequency-domain inversion to delineate and discover a mineralized zone in Porcupine district, Ontario, Canada. Vlad Kaminski* University of British
More informationFast Inversion of Logging-While-Drilling (LWD) Resistivity Measurements
Fast Inversion of Logging-While-Drilling (LWD) Resistivity Measurements David Pardo 1 Carlos Torres-Verdín 2 1 University of the Basque Country (UPV/EHU) and Ikerbasque, Bilbao, Spain. 2 The University
More informationSeismic wavepropagation concepts applied to the interpretation of marine controlled-source electromagnetics
Seismic wavepropagation concepts applied to the interpretation of marine controlled-source electromagnetics Rune Mittet, EMGS SUMMARY are very useful for the interpretation of seismic data. Moreover, these
More informationA Brief Revision of Vector Calculus and Maxwell s Equations
A Brief Revision of Vector Calculus and Maxwell s Equations Debapratim Ghosh Electronic Systems Group Department of Electrical Engineering Indian Institute of Technology Bombay e-mail: dghosh@ee.iitb.ac.in
More informationTHEORETICAL ANALYSIS AND NUMERICAL CALCULATION OF LOOP-SOURCE TRANSIENT ELECTROMAGNETIC IMAGING
CHINESE JOURNAL OF GEOPHYSICS Vol.47, No.2, 24, pp: 379 386 THEORETICAL ANALYSIS AND NUMERICAL CALCULATION OF LOOP-SOURCE TRANSIENT ELECTROMAGNETIC IMAGING XUE Guo-Qiang LI Xiu SONG Jian-Ping GUO Wen-Bo
More informationMicroscopic-Macroscopic connection. Silvana Botti
relating experiment and theory European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Temporary Address: Centre for Computational
More informationFinite-difference simulation of borehole EM measurements in 3D anisotropic media using coupled scalar-vector potentials
GEOPHYSICS VOL. 71 NO. 5 SEPTEMBER-OCTOBER 006 ; P. G5 G33 1 FIGS. 10.1190/1.567 Finite-difference simulation of borehole EM measurements in 3D anisotropic media using coupled scalar-vector potentials
More informationSEG Houston 2009 International Exposition and Annual Meeting
The role of EM rock physics and seismic data in integrated 3D CSEM data analysis I. Brevik*, StatoilHydro, Pål T. Gabrielsen, Vestfonna and Jan Petter Morten, EMGS Summary An extensive 3D CSEM dataset
More informationChap. 1 Fundamental Concepts
NE 2 Chap. 1 Fundamental Concepts Important Laws in Electromagnetics Coulomb s Law (1785) Gauss s Law (1839) Ampere s Law (1827) Ohm s Law (1827) Kirchhoff s Law (1845) Biot-Savart Law (1820) Faradays
More informationSUMMARY METHODOLOGY. 1 + (1 η)(ıωτ) c
On recovering induced polarization information from airborne time domain EM data Seogi Kang, Douglas W. Oldenburg, Dikun Yang, David Marchant, Geophysical Inversion Facility, University of British Columbia
More informationCompensating visco-acoustic effects in anisotropic resverse-time migration Sang Suh, Kwangjin Yoon, James Cai, and Bin Wang, TGS
Compensating visco-acoustic effects in anisotropic resverse-time migration Sang Suh, Kwangjin Yoon, James Cai, and Bin Wang, TGS SUMMARY Anelastic properties of the earth cause frequency dependent energy
More informationProblem Set 2 Due Tuesday, September 27, ; p : 0. (b) Construct a representation using five d orbitals that sit on the origin as a basis: 1
Problem Set 2 Due Tuesday, September 27, 211 Problems from Carter: Chapter 2: 2a-d,g,h,j 2.6, 2.9; Chapter 3: 1a-d,f,g 3.3, 3.6, 3.7 Additional problems: (1) Consider the D 4 point group and use a coordinate
More informationCourse no. 4. The Theory of Electromagnetic Field
Cose no. 4 The Theory of Electromagnetic Field Technical University of Cluj-Napoca http://www.et.utcluj.ro/cs_electromagnetics2006_ac.htm http://www.et.utcluj.ro/~lcret March 19-2009 Chapter 3 Magnetostatics
More informationMulti-source inversion of TEM data: with field applications to Mt. Milligan
Multi-source inversion of TEM data: with field applications to Mt. Milligan D. Oldenburg, E. Haber, D. Yang Geophysical Inversion Facility, University of British Columbia, Vancouver, BC, Canada Summary
More informationMapping Silicification in Gold Exploration
Mapping Silicification in Gold Exploration with Inductive Source Resistivity James Macnae and Patrick McGowan, Lamontagne Geophysics Ltd. Lamontagne Geophysics Ltd. 115 Grant Timmins Drive Kingston, Ontario
More informationSEG Houston 2009 International Exposition and Annual Meeting SUMMARY
Simultaneous joint inversion of MT and CSEM data using a multiplicative cost function Aria Abubakar, Maokun Li, Jianguo Liu and Tarek M. Habashy, Schlumberger-Doll Research, Cambridge, MA, USA SUMMARY
More information3 Constitutive Relations: Macroscopic Properties of Matter
EECS 53 Lecture 3 c Kamal Sarabandi Fall 21 All rights reserved 3 Constitutive Relations: Macroscopic Properties of Matter As shown previously, out of the four Maxwell s equations only the Faraday s and
More informationSpecial Section Marine Control-Source Electromagnetic Methods
GEOPHYSICS, VOL. 7, NO. MARCH-APRIL 7 ; P. WA63 WA7, FIGS. 9/43647 Special Section Marine Control-Source Electromagnetic Methods D marine controlled-source electromagnetic modeling: Part The effect of
More informationreservoir Summary A realistic controlled-sourcee Saturation simulations using Archie's is obtained by reservoir level. (Figure 1 and Table 1).
Sensitivity analysis of 3D surface-borehole McNeice*, EXPEC Advanced Researchh Center, Saudii Aramco, CSEM for a Saudi Arabian carbonate reservoir Daniele Colombo and Gary Garrett Kramer, WesternGeco Summary
More informationPrinciples of Applied Geophysics
Principles of Applied Geophysics Fifth edition D.S. Parasnis Professor emeritus of Applied Geophysics Department of Applied Geophysics, University ofluled, Lulea, Sweden Fellow of the Royal Swedish Academy
More informationElectromagnetic Fields About a Moving Magnetic Dipole: Implications for Marine Controlled Source Electromagnetic Systems
Index Table of contents Electromagnetic Fields About a Moving Magnetic Dipole: Implications for Marine Controlled Source Electromagnetic Systems E. E. S. Sampaio Centro De Pesquisa Em Geofísica E Geologia,
More informationA Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation
A Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation Tao Zhao 1, Feng-Nan Hwang 2 and Xiao-Chuan Cai 3 Abstract In this paper, we develop an overlapping domain decomposition
More informationDownloaded 10/02/15 to Redistribution subject to SEG license or copyright; see Terms of Use at
Feasibility study of application of nanoparticles in complex resistivity (CR) reservoir monitoring Vladimir Burtman*, University of Utah and TechnoImaging, Masashi Endo, TechnoImaging, Alberto Marsala,
More informationTowed Streamer EM data from Barents Sea, Norway
Towed Streamer EM data from Barents Sea, Norway Anwar Bhuiyan*, Eivind Vesterås and Allan McKay, PGS Summary The measured Towed Streamer EM data from a survey in the Barents Sea, undertaken in the Norwegian
More informationSummary. Introduction
Noel Black*, Glenn A. Wilson, TechnoImaging, Alexander V. Gribenko and Michael S. Zhdanov, TechnoImaging and The University of Utah Summary Recent studies have inferred the feasibility of time-lapse controlled-source
More informationDetermination of Locally Varying Directions through Mass Moment of Inertia Tensor
Determination of Locally Varying Directions through Mass Moment of Inertia Tensor R. M. Hassanpour and C.V. Deutsch Centre for Computational Geostatistics Department of Civil and Environmental Engineering
More informationSUMMARY INTRODUCTION. Traditionally, the inverse problem is approached by considering. DOI
Parametrized inversion framework for proppant volume in a hydraulically fractured reservoir Lindsey J. Heagy, A. Rowan Cockett, Douglas W. Oldenburg, UBC Geophysical Inversion Facility, University of British
More informationA cylinder in a magnetic field (Jackson)
Problem 1. A cylinder in a magnetic field (Jackson) A very long hollow cylinder of inner radius a and outer radius b of permeability µ is placed in an initially uniform magnetic field B o at right angles
More informationANALYSIS AND CORRECTION OF BOREHOLE EFFECT ON THE RESPONSES OF MULTICOMPONENT INDUCTION LOGGING TOOLS
Progress In Electromagnetics Research, PIER 85, 211 226, 2008 ANALYSIS AND CORRECTION OF BOREHOLE EFFECT ON THE RESPONSES OF MULTICOMPONENT INDUCTION LOGGING TOOLS X. Sun and Z. Nie School of Electronic
More informationSeismic processing of numerical EM data John W. Neese* and Leon Thomsen, University of Houston
Seismic processing of numerical EM data John W. Neese* and Leon Thomsen, University of Houston Summary The traditional methods for acquiring and processing CSEM data are very different from those for seismic
More informationl=0 The expansion coefficients can be determined, for example, by finding the potential on the z-axis and expanding that result in z.
Electrodynamics I Exam - Part A - Closed Book KSU 15/11/6 Name Electrodynamic Score = 14 / 14 points Instructions: Use SI units. Where appropriate, define all variables or symbols you use, in words. Try
More informationAn Approximate Hybrid Method for Modeling of Electromagnetic Scattering from an Underground Target. Shaaban Ali Bakr. PhD Thesis
An Approximate Hybrid Method for Modeling of Electromagnetic Scattering from an Underground Target PhD Thesis Shaaban Ali Bakr Department of Mathematics University of Bergen Norway August 2010 Preface
More information2 FORMULATIONS 2.1 The H-J formulation Let Ω be a domain consisting of a conducting region R and a non-conducting region S. Here assume that Ω, R and
Annual Report of ADVENTURE Project ADV-99-1(1999) LARGE-SCALE MAGNETIC FIELD ANALYSES Hiroshi KANAYAMA Department of Intelligent Machinery and Systems Graduate School of Engineering, Kyushu University
More informationVertical and horizontal resolution considerations for a joint 3D CSEM and MT inversion
Antony PRICE*, Total E&P and Don WATTS, WesternGeco Electromagnetics Summary To further explore the potential data content and inherent limitations of a detailed 3D Controlled Source ElectroMagnetic and
More informationNumerical Analysis of Electromagnetic Fields in Multiscale Model
Commun. Theor. Phys. 63 (205) 505 509 Vol. 63, No. 4, April, 205 Numerical Analysis of Electromagnetic Fields in Multiscale Model MA Ji ( ), FANG Guang-You (ྠ), and JI Yi-Cai (Π) Key Laboratory of Electromagnetic
More informationEvaluation of the Sacttering Matrix of Flat Dipoles Embedded in Multilayer Structures
PIERS ONLINE, VOL. 4, NO. 5, 2008 536 Evaluation of the Sacttering Matrix of Flat Dipoles Embedded in Multilayer Structures S. J. S. Sant Anna 1, 2, J. C. da S. Lacava 2, and D. Fernandes 2 1 Instituto
More informationA parallel method for large scale time domain electromagnetic inverse problems
A parallel method for large scale time domain electromagnetic inverse problems Eldad Haber July 15, 2005 Abstract In this work we consider the solution of 3D time domain electromagnetic inverse problems
More informationTheory of Electromagnetic Nondestructive Evaluation
EE 6XX Theory of Electromagnetic NDE: Theoretical Methods for Electromagnetic Nondestructive Evaluation 1915 Scholl Road CNDE Ames IA 50011 Graduate Tutorial Notes 2004 Theory of Electromagnetic Nondestructive
More informationA Novel Single-Source Surface Integral Method to Compute Scattering from Dielectric Objects
SUBMITTED TO IEEE ANTENNAS AND WIRELESS PROPAGATION LETTERS ON NOVEMBER 18, 2016 1 A Novel Single-Source Surface Integral Method to Compute Scattering from Dielectric Objects Utkarsh R. Patel, Student
More informationEstimating vertical and horizontal resistivity of the overburden and the reservoir for the Alvheim Boa field. Folke Engelmark* and Johan Mattsson, PGS
Estimating vertical and horizontal resistivity of the overburden and the reservoir for the Alvheim Boa field. Folke Engelmark* and Johan Mattsson, PGS Summary Towed streamer EM data was acquired in October
More informationPseudo-seismic wavelet transformation of transient electromagnetic response in engineering geology exploration
GEOPHYSICAL RESEARCH LETTERS, VOL. 34, L645, doi:.29/27gl36, 27 Pseudo-seismic wavelet transformation of transient electromagnetic response in engineering geology exploration G. Q. Xue, Y. J. Yan, 2 and
More informationEITN90 Radar and Remote Sensing Lecture 5: Target Reflectivity
EITN90 Radar and Remote Sensing Lecture 5: Target Reflectivity Daniel Sjöberg Department of Electrical and Information Technology Spring 2018 Outline 1 Basic reflection physics 2 Radar cross section definition
More informationSummary. 2D potential field migration of gravity fields and their gradients. For a 2D gravity field, we can define the complex intensity:
Michael S. Zhdanov*, The University of Utah and TechnoImaging, Xiaojun Liu, The University of Utah, and Glenn A. Wilson, TechnoImaging Summary We introduce 2D potential field migration for the interpretation
More informationMT Prospecting. Map Resistivity. Determine Formations. Determine Structure. Targeted Drilling
MT Prospecting Map Resistivity Determine Formations Determine Structure Targeted Drilling Cross-sectional interpretation before and after an MT survey of a mineral exploration prospect containing volcanic
More informationA multigrid integral equation method for large-scale models with inhomogeneous backgrounds
IOP PUBLISHING JOURNAL OF GEOPHYSICS AND ENGINEERING J. Geophys. Eng. (28) 438 447 doi:1.188/1742-2132//4/7 A integral equation method for large-scale models with inhomogeneous backgrounds Masashi Endo,
More informationConsider a Geonics EM31 with the TX and RX dipoles a distance s = 3 m apart.
E Near field, frequency domain EM systems E.1 EM responses of 1-D resistivity models E.1.1 Response of a halfspace Consider a Geonics EM31 with the TX and RX dipoles a distance s = 3 m apart. The TX and
More informationRotational & Rigid-Body Mechanics. Lectures 3+4
Rotational & Rigid-Body Mechanics Lectures 3+4 Rotational Motion So far: point objects moving through a trajectory. Next: moving actual dimensional objects and rotating them. 2 Circular Motion - Definitions
More informationUnit-1 Electrostatics-1
1. Describe about Co-ordinate Systems. Co-ordinate Systems Unit-1 Electrostatics-1 In order to describe the spatial variations of the quantities, we require using appropriate coordinate system. A point
More informationDownloaded 01/09/13 to Redistribution subject to SEG license or copyright; see Terms of Use at
Three-dimensional inversion of SQUID TEM data at Lalor Lake VMS deposit Dikun Yang*, Douglas W. Oldenburg Geophysical Inversion Facility, Earth and Ocean Sciences, University of British Columbia Summary
More informationG008 Advancing Marine Controlled Source Electromagnetics in the Santos Basin, Brazil
G008 Advancing Marine Controlled Source Electromagnetics in the Santos Basin, Brazil M.P. Buonora* (Petrobras), A. Zerilli (WesternGeco Electromagnetics), T. Labruzzo (WesternGeco Electromagnetics) & L.F.
More informationDownloaded 08/30/13 to Redistribution subject to SEG license or copyright; see Terms of Use at
Modeling the effect of pores and cracks interactions on the effective elastic properties of fractured porous rocks Luanxiao Zhao*, De-hua Han, Qiuliang Yao and Fuyong Yan, University of Houston; Mosab
More informationSCATTERING OF ELECTROMAGNETIC WAVES ON METAL NANOPARTICLES. Tomáš Váry, Juraj Chlpík, Peter Markoš
SCATTERING OF ELECTROMAGNETIC WAVES ON METAL NANOPARTICLES Tomáš Váry, Juraj Chlpík, Peter Markoš ÚJFI, FEI STU, Bratislava E-mail: tomas.vary@stuba.sk Received xx April 2012; accepted xx May 2012. 1.
More informationFAST FINITE VOLUME SIMULATION OF 3D ELECTROMAGNETIC PROBLEMS WITH HIGHLY DISCONTINUOUS COEFFICIENTS
SIAM J. SCI. COMPUT. Vol., No. 6, pp. 1943 1961 c 001 Society for Industrial and Applied Mathematics FAST FINITE VOLUME SIMULATION OF 3D ELECTROMAGNETIC PROBLEMS WITH HIGHLY DISCONTINUOUS COEFFICIENTS
More informationDownloaded 11/05/14 to Redistribution subject to SEG license or copyright; see Terms of Use at
Recovering a thin dipping conductor with 3D electromagnetic inversion over the Caber deposit Michael S. McMillan, Christoph Schwarzbach, Douglas W. Oldenburg and Eldad Haber, Geophysical Inversion Facility,
More informationMultiple Integrals and Vector Calculus: Synopsis
Multiple Integrals and Vector Calculus: Synopsis Hilary Term 28: 14 lectures. Steve Rawlings. 1. Vectors - recap of basic principles. Things which are (and are not) vectors. Differentiation and integration
More informationSummary. 3D potential field migration. Various approximations of the Newton method result in the imaging condition for the vector of densities, :
24 A case study from Broken Hill, Australia, with comparison to 3D regularized inversion Michael S. Zhdanov, University of Utah and TechnoImaging, Xiaojun Liu, University of Utah, Le Wan*, Martin Čuma,
More informationHydraulic properties of porous media
PART 5 Hydraulic properties of porous media Porosity Definition: Void space: n V void /V total total porosity e V void /V solid Primary porosity - between grains Secondary porosity - fracture or solution
More informationUniversity of Illinois at Chicago Department of Physics. Electricity and Magnetism PhD Qualifying Examination
University of Illinois at Chicago Department of Physics Electricity and Magnetism PhD Qualifying Examination January 8, 216 (Friday) 9: am - 12: noon Full credit can be achieved from completely correct
More informationProblem Set 2 Due Thursday, October 1, & & & & # % (b) Construct a representation using five d orbitals that sit on the origin as a basis:
Problem Set 2 Due Thursday, October 1, 29 Problems from Cotton: Chapter 4: 4.6, 4.7; Chapter 6: 6.2, 6.4, 6.5 Additional problems: (1) Consider the D 3h point group and use a coordinate system wherein
More informationFREQUENCY DEPENDENT CHARACTERISTICS OF GROUNDING SYSTEM BURIED IN MULTILAYERED EARTH MODEL BASED ON QUASI-STATIC ELECTRO- MAGNETIC FIELD THEORY
Progress In Electromagnetics Research M, Vol. 33, 169 183, 2013 FREQUENCY DEPENDENT CHARACTERISTICS OF GROUNDING SYSTEM BURIED IN MULTILAYERED EARTH MODEL BASED ON QUASI-STATIC ELECTRO- MAGNETIC FIELD
More informationLecture 7. Properties of Materials
MIT 3.00 Fall 2002 c W.C Carter 55 Lecture 7 Properties of Materials Last Time Types of Systems and Types of Processes Division of Total Energy into Kinetic, Potential, and Internal Types of Work: Polarization
More informationA new formula to compute apparent resistivities from marine magnetometric resistivity data
GEOPHYSICS, VOL. 71, NO. 3 MAY-JUNE 006; P. G73 G81, 10 FIGS. 10.1190/1.194513 A new formula to compute apparent resistivities from marine magnetometric resistivity data Jiuping Chen 1 Douglas W. Oldenburg
More informationModeling of 3D MCSEM and Sensitivity Analysis
PIERS ONLINE, VOL. 3, NO. 5, 2007 641 Modeling of 3D MCSEM and Sensitivity Analysis Zhanxiang He 1, 2, Zhigang Wang 2, Gang Yu 3, Kurt Strack 3, and Haiying Liu 2 1 Department of Geosciences, University
More informationE&M. 1 Capacitors. January 2009
E&M January 2009 1 Capacitors Consider a spherical capacitor which has the space between its plates filled with a dielectric of permittivity ɛ. The inner sphere has radius r 1 and the outer sphere has
More informationElectromagnetic wave propagation. ELEC 041-Modeling and design of electromagnetic systems
Electromagnetic wave propagation ELEC 041-Modeling and design of electromagnetic systems EM wave propagation In general, open problems with a computation domain extending (in theory) to infinity not bounded
More informationUniversity of Illinois at Chicago Department of Physics
University of Illinois at Chicago Department of Physics Electromagnetism Qualifying Examination January 4, 2017 9.00 am - 12.00 pm Full credit can be achieved from completely correct answers to 4 questions.
More information13, Applications of molecular symmetry and group theory
Subject Paper No and Title Module No and Title Module Tag Chemistry 13, Applications of molecular symmetry and group theory 27, Group theory and vibrational spectroscopy: Part-IV(Selection rules for IR
More informationTheory and Applications of Dielectric Materials Introduction
SERG Summer Seminar Series #11 Theory and Applications of Dielectric Materials Introduction Tzuyang Yu Associate Professor, Ph.D. Structural Engineering Research Group (SERG) Department of Civil and Environmental
More informationD034 On the Effects of Anisotropy in Marine CSEM
D034 On the Effects of Anisotropy in Marine CSEM L.O. Løseth* (Statoil ASA, formerly Dept. of Physics, NTNU), B. Ursin (Dept. of Appl. Geophys., NTNU) & L. Amundsen (Statoil ASA) SUMMARY In marine CSEM
More informationEnhancing and suppressing radiation with some permeability-near-zero structures
Enhancing and suppressing radiation with some permeability-near-zero structures Yi Jin 1,2 and Sailing He 1,2,3,* 1 Centre for Optical and Electromagnetic Research, State Key Laboratory of Modern Optical
More informationTime Resolved Faraday Rotation Measurements of Spin Polarized Currents in Quantum Wells
Time Resolved Faraday Rotation Measurements of Spin Polarized Currents in Quantum Wells M. R. Beversluis 17 December 2001 1 Introduction For over thirty years, silicon based electronics have continued
More informationSpin Interactions. Giuseppe Pileio 24/10/2006
Spin Interactions Giuseppe Pileio 24/10/2006 Magnetic moment µ = " I ˆ µ = " h I(I +1) " = g# h Spin interactions overview Zeeman Interaction Zeeman interaction Interaction with the static magnetic field
More informationWaves in Linear Optical Media
1/53 Waves in Linear Optical Media Sergey A. Ponomarenko Dalhousie University c 2009 S. A. Ponomarenko Outline Plane waves in free space. Polarization. Plane waves in linear lossy media. Dispersion relations
More informationMaxwell s Equations:
Course Instructor Dr. Raymond C. Rumpf Office: A-337 Phone: (915) 747-6958 E-Mail: rcrumpf@utep.edu Maxwell s Equations: Physical Interpretation EE-3321 Electromagnetic Field Theory Outline Maxwell s Equations
More informationA DARK GREY P O N T, with a Switch Tail, and a small Star on the Forehead. Any
Y Y Y X X «/ YY Y Y ««Y x ) & \ & & } # Y \#$& / Y Y X» \\ / X X X x & Y Y X «q «z \x» = q Y # % \ & [ & Z \ & { + % ) / / «q zy» / & / / / & x x X / % % ) Y x X Y $ Z % Y Y x x } / % «] «] # z» & Y X»
More informationTopological insulator part I: Phenomena
Phys60.nb 5 Topological insulator part I: Phenomena (Part II and Part III discusses how to understand a topological insluator based band-structure theory and gauge theory) (Part IV discusses more complicated
More informationSEAM2: 3D NON-SEISMIC MODELING OF A COMPLEX MIDDLE EAST O&G PROSPECT
6th International Symposium on Three-Dimensional Electromagnetics Berkeley, California, USA, March 28-3, 217 SEAM2: 3D NON-SEISMIC MODELING OF A COMPLEX MIDDLE EAST O&G PROSPECT M. D. Watts 1, R. L. Mackie
More informationElectromagnetic Theorems
Electromagnetic Theorems Daniel S. Weile Department of Electrical and Computer Engineering University of Delaware ELEG 648 Electromagnetic Theorems Outline Outline Duality The Main Idea Electric Sources
More informationMultiple Integrals and Vector Calculus (Oxford Physics) Synopsis and Problem Sets; Hilary 2015
Multiple Integrals and Vector Calculus (Oxford Physics) Ramin Golestanian Synopsis and Problem Sets; Hilary 215 The outline of the material, which will be covered in 14 lectures, is as follows: 1. Introduction
More informationTHREE DIMENSIONAL CONTROLLED-SOURCE ELECTROMAGNETIC EDGE-BASED FINITE ELEMENT MODELING OF CONDUCTIVE AND PERMEABLE HETEROGENEITIES.
THREE DIMENSIONAL CONTROLLED-SOURCE ELECTROMAGNETIC EDGE-BASED FINITE ELEMENT MODELING OF CONDUCTIVE AND PERMEABLE HETEROGENEITIES A Dissertation by SOUVIK MUKHERJEE Submitted to the Office of Graduate
More informationScattering of ECRF waves by edge density fluctuations and blobs
PSFC/JA-14-7 Scattering of ECRF waves by edge density fluctuations and blobs A. K. Ram and K. Hizanidis a June 2014 Plasma Science and Fusion Center, Massachusetts Institute of Technology Cambridge, MA
More informationMAE 323: Lecture 1. Review
This review is divided into two parts. The first part is a mini-review of statics and solid mechanics. The second part is a review of matrix/vector fundamentals. The first part is given as an refresher
More information3. Magnetic Methods / 62
Contents Preface to the Second Edition / xv Excerpts from Preface to the FirstEdition / xvii Mathematical Conventions / xix 1. Introduction / 1 Reference / 5 2. Gravity Methods / 6 2. I. Introduction /
More informationChem8028(1314) - Spin Dynamics: Spin Interactions
Chem8028(1314) - Spin Dynamics: Spin Interactions Malcolm Levitt see also IK m106 1 Nuclear spin interactions (diamagnetic materials) 2 Chemical Shift 3 Direct dipole-dipole coupling 4 J-coupling 5 Nuclear
More informationElectricity & Magnetism Study Questions for the Spring 2018 Department Exam December 4, 2017
Electricity & Magnetism Study Questions for the Spring 2018 Department Exam December 4, 2017 1. a. Find the capacitance of a spherical capacitor with inner radius l i and outer radius l 0 filled with dielectric
More information