Volume 44 Number 6
  • E-ISSN: 1365-2478



Hydrocarbon exploration in the Papuan fold belt is made extremely difficult by mountainous terrain, equatorial jungle and thick karstified Miocene limestones at the surface. The high‐velocity karstified limestones at or near the surface often render the seismic technique useless for imaging the subsurface. In such areas magnetotellurics (MT) provides a valuable capability for mapping subsurface structure. The main structural interface which can be mapped with MT, due to the large electrical contrast, is the contact between the resistive Darai limestone and the underlying conductive sediments of the Ieru Formation. In some areas the base of the Darai can be mapped with reasonable accuracy by fitting 1D models to the observed MT data. However, in many cases where 2D and 3D effects are severe, 1D interpretations can yield dramatically incorrect results. Numerical and field data examples are presented which demonstrate the severity of the 1D errors and the improvements in accuracy which can be achieved using a 2D inverse solution.

Two MT lines over adjacent anticlines, both with well control and seismic data, are used to demonstrate the application of 1D and 2D inversions for structural models. In both cases the seismic data provide no aid in the interpretations. The example over the Hides anticline illustrates a situation where 1D inversion of either TE or TM mode provides essentially the same depth to base of Darai as 2D inversion of both TE and TM. Both models provide base Darai depth estimates which are within 10% of that measured in the well. The example over the Angore anticline illustrates the inadequacy of 1D inversion in structurally complex geology complicated by electrical statics. The TE mode fits a 1D Darai thickness of 800 metres while the TM mode fits a 1D Darai thickness of 3500 metres, bracketing the thickness of 2450 metres observed in the well. The final 2D inversion model provides a depth estimate of 2250 metres. Four MT lines along the Angore anticline have been interpreted using 2D inversion. A high degree of correlation exists between lineaments observed on an airborne radar image and zones of low resistivity within the high‐resistivity material interpreted as Darai limestone. These low‐resistivity zones are interpreted as fault zones.

Three‐dimensional modelling has been used to simulate 3D statics in an otherwise 2D earth. These data were used to test the Groom‐Bailey (GB) decomposition for possible benefits in reducing static effects and estimating geoelectric strike in the Papua New Guinea (PNG) field data. It has been found that the GB decomposition can provide improved regional 2D strike estimates in 3D contaminated data. However, in situations such as PNG, where the regional 2D strike is well established and hence can be fixed, the GB decomposition provides apparent resistivities identical to those simply rotated to strike.


Article metrics loading...

Loading full text...

Full text loading...


  1. BahrK.1988. Interpretation of the magnetotelluric impedance tensor; regional induction and local telluric distortion. Journal of Geophysics62, 119–127.
    [Google Scholar]
  2. BennettJ.K. and ParishM.1993. The use of magneto‐telluric surveying as a risk reduction tool for hydrocarbon exploration in the Papuan fold belt. Petroleum Exploration in Papua New Guinea. Proceedings of 2nd PNG petroleum convention, Port Moresby, pp. 335–350.
  3. BillingsA.J. and ThomasJ.H.1990. The use and limitations of non‐seismic geophysics in the Papuan thrust belt: Petroleum Exploration in Papua New Guinea. Proceedings of the 1st PNG petroleum convention, Port Moresby, pp. 51–62.
  4. ChristophersonK.R.1991. Applications of magnetotellurics to petroleum exploration in Papua New Guinea; a model for frontier areas. The Leading Edge10, 21–27.
    [Google Scholar]
  5. DeGroot‐HedlinC. and ConstableS.1990. Occam's inversion to generate smooth, two‐dimensional models from magnetotelluric data. Geophysics55, 1613–1624.
    [Google Scholar]
  6. GroomR.W. and BaileyR.C.1991. Analytic investigations of the effects of near‐surface three‐dimensional galvanic scatters on MT tensor decompositions. Geophysics56, 496–518.
    [Google Scholar]
  7. Hill, K.C.1991. Structure of the Papuan fold belt, Papua New Guinea. AAPG Bulletin56, 857–872.
    [Google Scholar]
  8. MaddenT.R. and MackieR.L.1989. Three‐dimensional magnetotelluric modeling and inversion. Proceedings of the IEEE77, 318–333.
    [Google Scholar]
  9. OrangeA. S.1989. Magnetotelluric exploration for hydrocarbons. Proceedings of the IEEE77, 287–317.
    [Google Scholar]
  10. ReddigR.1991. WSE Associates Inc. Electrical Methods Course. Unpublished.
  11. RossL.1993. Evolution of structural models of two anticlines in the Papuan thrust belt by application of magnetotellurics. Petroleum Exploration and Development in Papua New Guinea. Proceedings of the 2nd PNG petroleum convention, Port Moresby, pp. 351–370.
  12. SternbergB.K., WashburneJ.C. and PellerinL.1988. Correction for the static shift in magnetotellurics using transient electromagnetic soundings. Geophysics53, 1459–1468.
    [Google Scholar]
  13. Torres‐VerdinC. and Bostick, F.X.Jr.1992. Principles of spatial surface electric field filtering in magnetotellurics; Electromagnetic array profiling (EMAP). Geophysics57, 603–622.
    [Google Scholar]
  14. WannamakerP.E., HohmannG.W. and Ward, S.H.1984. Magnetotelluric responses of threedimensional bodies in layered eafihs. Geophysics49, 1517–1533.
    [Google Scholar]
  15. WannamakerP.E., StodtJ.A. and RijoL.1987. A stable finite‐element solution for twodimensional magnetotelluric modelling. Geophysical Journal of the Royal Astronomical Society88, 277–296.
    [Google Scholar]
  • Article Type: Research Article
This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error