1 / 51

SAGE MT

SAGE MT. George R. Jiracek San Diego State University. "Understanding is More Important Than Knowledge". THE INPUT. THE OUTPUT. MT DATA. LIGHTNING. SOLAR WIND. BLACK BOX EARTH. MT Data Collection. Marlborough, New Zealand. Southern Alps, New Zealand. Southern Alps, New Zealand.

trixie
Télécharger la présentation

SAGE MT

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. SAGE MT George R. Jiracek San Diego State University "Understanding is More Important Than Knowledge"

  2. THE INPUT THEOUTPUT MT DATA LIGHTNING SOLAR WIND BLACK BOX EARTH

  3. MT Data Collection

  4. Marlborough, New Zealand

  5. Southern Alps, New Zealand

  6. Southern Alps, New Zealand

  7. Southern Alps, New Zealand

  8. Taupo, New Zealand 2010-12

  9. Southern Alps, New Zealand The “Banana”

  10. Southern Alps, New Zealand (Jiracek et al., 2007)

  11. Southern Alps, New Zealand New Zealand Earthquakes vs. Resistivity in Three-Dimensions

  12. Three-Dimensional MT Taupo Volcanic Geothermal Field, New Zealand (Heise et al. , 2008)

  13. MT Phase Tensor Plot at 0.67s Period from the Taupo Volcanic Field

  14. Magnetotellurics (MT) • Low frequency (VLF to subHertz) • Natural source technique • Energy diffusion governed by ρ(x,y,z) (Ack. Paul Bedrosian, USGS) Techniques - MT

  15. Magnetotelluric Signals (Ack. Paul Bedrosian, USGS) Techniques - MT

  16. Quasi-static approx, σ >> εω Always Must SatisfyMaxwell’s Equations rfis free charge density (Ack. Paul Bedrosian, USGS) Magnetotellurics

  17. Quasistatic Approximation d is skin depth (Ack. Paul Bedrosian, USGS)

  18. Graphical Description of Skin Depth, d

  19. Magnetotelluric Impedance After Fourier transforming the E(t) and H(t) data into the frequency domain the MT surface impedance is calculated from: Ex(w) = Z(w) Hy(w)

  20. Note, that since Ex(w) = Z(w) Hy(w) is a multiplication in the frequency domain, it is a convolution in the time domain. Therefore, this is a filtering operation, i.e., Ex(t) Hy(t) Z(t)

  21. Apparent resistivity, ra and phase, f Apparent resistivity is the resistivity of an equivalent, but fictitious, homogeneous, isotropic half-space Phase is phase of the impedance f = tan-1 (Im Z/Re Z)

  22. The goal of MT is the resistivity distribution, r(x,y,z), of the subsurface as calculated from the surface electromagnetic impedance, Zs r7 r1 r2 Dimensionality: r3 • One-Dimensional • Two-Dimensional • Three-Dimensional r4 r5 r6

  23. Geoelectric Dimensionality 3-D 1-D 2-D

  24. x Shallow Resistive Layer Intermediate Conductive Layer Deep Resistive Layer y ra Log z Log Period(s) 1-D MT Sounding Curve raa |Z2|

  25. 104 Ex 100 103 Hy 500 102 Ohm-m 30 101 Apparent resistivity 80 Impedance Phase 60 Degrees 40 20 0 10-2 100 102 104 Period (s) Layered (1-D) Earth 1000 Longer period  deeper penetration ( )m Using a range of periods a depth sounding can be obtained (Ack., Paul Bedrosian, USGS)

  26. MT “Screening” of Deep Conductive Layer by Shallow Conductive Layer (Ack., Martyn Unsworth, Univ. Alberta)

  27. When the Earth is either 2-D or 3-D: Ex(w) = Z(w) Hy(w) Now Ex(w) = Zxx(w) Hx(w) + Zxy(w) Hy(w) Ey(w) = Zyx(w) Hx(w) + Zyy(w) Hy(w) This defines the tensor impedance, Z(w)

  28. 3-D MT Tensor Equation

  29. 1-D, 2-D, and 3-D Impedance • 1-D • 2-D • Assumes geoelectric strike • 3-D • No geoelectric assumptions [ ] is Tensor Impedance (Ack., Paul Bedrosian, USGS)

  30. 3- D MT Data Measure time variations of electric (E) and magnetic (H) fields at the Earth‘s surface. Estimate transfer functions of the E and H fields. Subsurface resistivity distribution recovered through modeling and inversion. Impedance Tensor: App Resistivity & Phase: (Ack. Paul Bedrosian, USGS) Techniques - MT

  31. x y ra Log z Log Period (s) 2-D MT (Tensor Impedance reduces to two off- diagonal elements) Geoelectric Strike raa |Z2|

  32. TM TE Boundary Conditions • E-Fields parallel to the geoelectric strike are continuous (called TE mode) • E-Fields perpendicular to the geoelectric strike are discontinuous (called TM mode) Map View E-Parallel Log ra E- Perpendicular Log Period (s)

  33. TE (Transverse Electric) and TM (Transverse Magnetic) Modes MT2 • 2-D Earth structure • Different results at MT1 (Ex and Hy) • and MT2 (Ey and Hx) MT1 TRANSVERSE ELECTRIC MODE (TE) TRANSVERSE MAGNETIC MODE (TM) Visualizing Maxwell’s Curl Equations (Ack., Martyn Unsworth, Univ. Alberta)

  34. MT Phase Tensor Described as “elegant” by Berdichevsky and Dmitriev (2008) and a “major breakthrough” by Weidelt and Chave (2012) “Despite its deceiving simplicity, students attending the SAGE program often have problems grasping the essence of the MT phase tensor” (Jiracek et al., 2014) The MT Phase Tensor and its Relation to MT Distortion (Jiracek Draft, June, 2014)

  35. MT Phase Tensor • X and Y are the real and imaginary parts of impedance tensor Z, i.e., Z = X + iY • Ideal 2-D, β=0 • Recommended β <3° for ~ 2-D by Caldwell et al., (2004)

  36. MT Phase Tensor Ellipse • Ellipses are traced out at every period by the multiplication of • the real 2 x 2 matrix from a MT phase tensor, F(f) and • arotating, family of unit vectors, c(w), that describe a unit circle. 2-D Tensor Ellipse p2D(w) is: http://www-rohan.sdsu.edu/~jiracek/DAGSAW/Rotation_Figure/

  37. Phase Tensor Example for Single MT Sounding at Taupo Volcanic Field, New Zealand (Bibby et al., 2005) 1-D TP Tc 2-D TP Tc 2-D TP

  38. Tc Tc 1-D TP 2-D TP 2-D TP Phase Tensor Determinations of Dimensionality (1-D. 2-D), Transition Periods (TP), and Threshold Periods (Tc)

  39. SAGE MT Caja Del Rio

  40. Geoelectric Section From Stitched 1-D TE Inversions (MT Sites Indicated by Triangles) E W Conductive Basin Elevation (m) Resistive Basement Distance (m) W E Basin Elevation (m) Basement

  41. 2-D MT Inversion/Finite-Difference Grid • M model parameters, N surface measurements, M>>N • A regularized solution narrows the model subspace • Introduce constraints on the smoothness of the model (Ack. Paul Bedrosian, USGS) Techniques - MT

  42. Geoelectric Section From 2-D MT Inversion (MT Sites Indicated by Triangles) W E Conductive Basin Elevation (m) Resistive Basement Distance (m) W E Basin Elevation (m) Basement

  43. SAGE – Rio Grande Rift, New Mexico (Winther, 2009)

  44. Resistivity Values of Earth Materials

  45. MT Interpretation Geology Well Logs

  46. SAGE – Rio Grande Rift, New Mexico (Winther, 2009)

  47. MT-Derived Midcrustal Conductor Physical State Eastern Great Basin (EGB), Transition Zone (TZ), and Colorado Plateau (CP) (Wannamaker et al., 2008)

  48. Field Area Now The Future?

  49. References Bibby, H. M., T. G. Caldwell, and C. Brown, 2005, Determinable and non-determinable parameters of galvanic distortion in magnetotellurics, Geophys. J. Int., 163, 915 -930. Caldwell, T. G., H. M.Bibby, and C. Brown, 2004, The magnetotelluric phase tensor, Geophys. J. Int., 158, 457- 469. Heise, W., T. G. Caldwell, H. W. Bibby, and C. Brown, 2006, Anisotropy and phase splits in magnetotellurics, Phys. Earth. Planet. Inter., 158, 107-121. Jiracek, G.R., V. Haak, and K.H. Olsen, 1995, Practical magnetotellurics in continental rift environments, in Continental rifts: evolution, structure, and tectonics, K.H. Olsen, ed., 103-129. Jiracek, G. R., V. M Gonzalez, T. G. Caldwell, P. E. Wannamaker, and D. Kilb, 2007, Seismogenic, Electrically Conductive, and Fluid Zones at Continental Plate Boundaries in New Zealand, Himalaya, and California-USA, in Tectonics of A Continental Transform Plate Boundary: The South Island, New Zealand, Amer. Geophys. Un. Mono. Ser. 175, 347-369.

More Related