MARE2DEM can model isotropic to triaxially anisotropic conductivity but tilted or generally anisotropic media is not supported. Switching between the various types of anisotropy is easily accomplished using the anisotropy menu in the Mamba2D model building assistant. For most modeling purposes, you will likely want to use simple isotropic conductivity, but for modeling a stack of interbedded marine sediments the transversely isotropic vertical (TIZ) option can be useful. The other options are for more specialized anisotropies. In general if you have no idea about anisotropy or which setting you should use, go for the isotropic default.
Here's a guide to the five possible options:
isotropic - conductivity is the same in all direction of electric current flow. The conductivity tensor has the form:
triaxial - conductivity varies in the direction of the three coordinate axes. The conductivity tensor has the form:
where are the 2D model axes.
There are three other possible anisotropies where the conductivity along one axis is different than the other two. This type of anisotropy is generally referred to as transverse isotropy, where conductivity is symmetric about an axis normal to the plane of isotropy. MARE2DEM uses the abbreviation TI for transverse isotropy and the third letter denotes the transverse axis. The three possible cases are:
TIX - transversely isotropic perpendicular to the x axis. The conductivity tensor has the form:
where denotes conductivity along the x axis, and is conductivity in the transverse plane.
TIY - transversely isotropic perpendicular to the y axis. The conductivity tensor has the form:
TIZ - transversely isotropic perpendicular to the z axis. The conductivity tensor has the form:
This is also known in the exploration community as vertically transverse isotropy (VTI).
Note that when you specify the resistivity for anisotropic parameters using Mamba2D.m, the edit box expects a specific order for the components, as shown below for the TIZ anisotropy commonly used for marine sediments. In this case, the resistivity menu expects the resistivity for the transverse component (vertical resistivity ) to be listed before the horizontal resistivity. In the example below, the vertical to horizontal resistivity anisotropy factor is 3 since the vertical resistivity is 6 and the horizontal resistivity is 2.
For triaxial resistivity, the edit menu expects resistivity in order of :
Receiver rotations are described by three angles: θ, α and β, as shown in the image below. These rotate from the model coordinates to the fully rotated measurement coordinates . For the rotation matrices associated with these angles, see Key and Lockwood (2010). The θ angle describes the horizontal rotation of the receiver's axisfrom the 2D model strike direction . For a typical CSEM data set , the data should be rotated so that θ = 0. Therefore, for an array of receivers inline to a horizontal electric dipole source (with transmitter azimuth of 90°), the inline electric field will be the component. There angles α and β describe tilts of the receiver. α is the vertical tilt of the component. If α =0, then β describes the vertical tilt of the , otherwise β describes a slanted angle to .
MARE2DEM allows for arbitrary orientations of transmitter dipoles. The orientation of electric and magnetic dipole transmitters are described by the two angles of the dipole axis: the horizontal azimuth and the vertical dip angle.
For MT receivers, the x coordinate of the receivers is ignored by MARE2DEM since the 2D MT fields are strike invariant. For controlled-source electromagnetic receivers, there are a few considerations you need to mind about the accuracy of MARE2DEM, which uses a total field implementation for the source and a wavenumber domain transformation to compute the spatially varying EM fields:
The receiver x component should be close the to transmitter x position. Receivers don't need to be perfectly inline, but if they are more than 1-2 km down-strike from the transmitter (gray regions below), the default settings for the wavenumber domain transforms in MARE2DEM may break down due to highly oscillatory kernel functions.
Due to the source singularity in the 2.5D wavenumber domain, receivers located down strike from the transmitter (pink region below) may be inaccurate unless they are positioned much deeper than the source.
If modeling point dipoles, remember that receivers closer than a few real dipole lengths from the transmitter will be inaccurate due to finite length dipole effects (red region below). You can instead specify a finite length dipole in the data file, but this comes at an increased numerical effort for MARE2DEM (more on that in another tutorial).
Note: UTM referencing is only used for plotting purposes. MARE2DEM completely ignores this for the 2D EM calculations.
To help with post-inversion plotting of well logs and SEG-Y seismic overlays on your 2D inversion models, the data file has a line for specifying the UTM zone, origin and rotation of the 2D model coordinates:
The 2D strike angle is the geographic angle corresponding to the x direction of the 2D model, clockwise from north. The example below has a 2D strike of 20º.