Mon Taj Induced Polarization
Mon Taj Induced Polarization
Mon Taj Induced Polarization
www.geosoft.com
The software described in this manual is furnished under license and may
only be used or copied in accordance with the terms of the license.
Written by, Nancy Whitehead and Chris Musselman. Please send comments
or questions to info@geosoft.com
Geosoft Incorporated
8th Floor
85 Richmond St. W.
Toronto, Ontario
M5H 2C9
Canada
Tel: (416) 369-0111
Fax: (416) 369-9599
Processing Sequence 6
Array Channels 9
Database Channels 10
References 70
1
1. DEFINITIONS:
In this Agreement:
"Licensed Program(s)" means the actual copy of all or any portion of Geosoft’s proprietary software technology, computer
software code, components, dynamic link libraries (DLLs) licensed through the Geosoft license server, including any
modifications, improvements or updates provided by GEOSOFT.
“Effective Date” is the date the Geosoft license is installed. This date is recorded by the Geosoft License server when the
Licensed Program(s) is installed.
"Services" means the Services described on Section 4.
"Termination" means the occurrences contemplated by Section 6 and 7.
2. LICENSE:
GEOSOFT grants to me a non-transferable and non-exclusive license to use the Licensed Program(s) for my own purposes
whereby the Licensed Program(s) are being used only by myself, on one computer, at any one time.
Title and all intellectual property rights in and to the License Program(s), including, without limitation, copyright, trade secrets
and trade marks, shall remain with GEOSOFT. I agree to refrain from raising any objection or challenge to such intellectual
property rights, or from assisting or causing or permitting other(s) to do so, during the term of the Agreement and thereafter
I may not assign this Agreement or any part thereof or sub-license the rights granted herein, or lend, rent, time-share, sell or
lease the software without the prior written consent of GEOSOFT.
I may not attempt to reverse engineer, de-compile or disassemble the software.
I may not make any attempt to circumvent the License Manager that controls the access to the software use.
3. TERM:
The Term of this Agreement shall commence on the Effective Date and shall continue until termination, as described in Section
6.
4. SERVICES:
(i) According to the terms of my initial purchase, GEOSOFT shall make available to me, without additional fees such
corrections and improvements to the Licensed Program(s) as may be generally incorporated into the Licensed Program(s) by
GEOSOFT. (Normally this will be for a period of twelve (12) months).
(ii) GEOSOFT has a strong commitment to customer service and product support. GEOSOFT offers me, subject to applicable
Service Charge(s), continuing support in the form of email or telephone advice and other assistance in problem diagnosis and
the correction of errors or faults in the Licensed Program(s) during the life of this License. When a problem occurs which
appears to be related to errors or faults in the Licensed Program(s), I may contact GEOSOFT and GEOSOFT will make an
honest effort to solve the problem. However, GEOSOFT cannot guarantee service results or represent or warrant that all errors
or program defects will be corrected. Also it is to be noted that each Licensed Program is designed to operate on a Windows
NT (sp 6 or later), Windows 2000 or Windows XP platform.
(iii) Further, if I request service relating to the modification of the Licensed Program(s) to meet a particular need or to conform
with a particular operating environment, GEOSOFT may, at its discretion, modify the Licensed Program(s) to meet these
particular needs, subject to applicable Services Charge(s). However, all intellectual property or other rights which may arise
from such modifications shall reside with GEOSOFT.
6. TERMINATION:
This agreement shall terminate upon the termination date, if any, specified in your purchase agreement with Geosoft.
This agreement may be terminated only upon thirty-days prior written notice to GEOSOFT.
GEOSOFT may terminate this Agreement upon prior written notice effective immediately if I fail to comply with any of the terms
and conditions of this Agreement.
This Agreement shall terminate automatically upon the institution, or consenting to the institution of proceedings in insolvency
or bankruptcy, or upon a trustee in bankruptcy or receiver being appointed for me/us for all or a substantial portion of my/our
assets.
8. WARRANTY:
GEOSOFT does not warrant that the functions contained in the Licensed Program will meet my requirements or will operate in
the combinations which may be selected for use by me, or that the operation of the Licensed Program will be uninterrupted or
error free or that all program defects will be corrected.
Each Licensed Program shall be furnished to me in accordance with the terms of this Agreement. No warranties, either express
or implied, are made to me regarding the Licensed Program.
THE FOREGOING WARRANTIES ARE IN LIEU OF ALL OTHER WARRANTIES, EXPRESSED OR IMPLIED, INCLUDING,
BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OR MERCHANTABILITY AND FITNESS FOR A PARTICULAR
PURPOSE.
9. LIMITATION OF REMEDIES
I agree to accept responsibility for the use of the programs to achieve my intended results, and for the results obtained from
use of said Program(s). I therefore accept complete responsibility for any decision made based on my use of the
aforementioned Licensed Program(s).
In no event shall GEOSOFT be liable for any damages arising from performance or non-performance of the Licensed
Program(s), or for any lost profits, lost savings or other consequential damages, even if GEOSOFT has been advised of the
possibility of such damages, or for any claim against me by any other party.
10. GENERAL:
I agree that this Agreement is a complete and exclusive statement of the agreement with GEOSOFT.
This Agreement supersedes all previous Agreements with respect to the Licensed Programs, with the exception of a current
signed Technical Service Agreements.
GEOSOFT is not responsible for failure to fulfill its obligations under the Agreement due to causes beyond its control.
Should any part of This Agreement for any reason be declared invalid, such declaration shall not affect the remaining portion
which shall remain in full force and effect as if this Agreement had been executed without the invalid portion thereof.
The relationship between the parties is that of independent contractors. Nothing contained in this Agreement shall be deemed
to constitute or create a partnership, association, joint venture or agency.
The provision of this Agreement shall be binding upon me and GEOSOFT and my respective successors and permitted
assigns.
This Agreement will be governed by the laws of the Province of Ontario and applicable laws of Canada.
The best way to find information in this system is to use the Search tab to perform a
full-text search across all help topics. If you still cannot find the information you are
looking for, the Online Books help system contains complete Geosoft manuals and
tutorials in Adobe PDF format.
The inherent differences between time domain and frequency domain are minimized
in the montaj Induced Polarization system by providing a common database
environment and processing tools for both survey types. This means that interpreters
can perform the same import and quality control procedures regardless of the type of
data. By processing and presenting data in a similar way, it is possible to reduce the
time required to learn the system and produce results.
The montaj Induced Polarization system imports data from a variety of instrument
dump files, IP files, and array configurations, including:
Multiple values obtained for the same measurement point (time windows or
frequency responses) are displayed as a curve (array) in a single cell of the database.
This way, the system handles data channels containing both single values and value
arrays.
Data values for duplicate samples can be edited interactively for each line and station
in the database. The Quality Control tool in the montaj Induced Polarization system
enables you to visualize the relative values for each sample and to specify which
samples are used for calculations, export, and plotting. The system uses a channel to
identify the status of each sample, so the original data are not modified.
Processing Sequence
The general processing steps for time domain and frequency domain data are the
same. The following diagram identifies the main processing tasks and describes what
the software does during each step. The procedures for these tasks are described in
Chapter 2: Quick Start Tutorial on page 17.
Chapter 1: System Capabilities and IP Concepts 7
The following image identifies the elements of the IP spreadsheet and profile panes
for a database containing time-domain survey data on an N-S line.
Single channel of windowed Apparent
Station Station Position
Line Number data (from array channel IP) Metal Factor Resistivity
Number (calculated) (same as X)
(N-S line Survey QC flag
(N-S line IP Average Induced
direction) Elevation Receiver
direction) Current
Survey Line Time Window Positions Self-Potential
“n” value
Number (array channel)
Fiducials
(survey
reading Spreadsheet
numbers) Pane
Duplicate
samples
Profile Data
Range Profile
Pane
Currently
Selected
Channel
The spreadsheet displays a database containing all the survey data. The spreadsheet
organises information for a survey line using rows and columns. The current survey
line number is located in the top left cell of the spreadsheet. For IP data, each
spreadsheet row contains a fiducial number. The fiducial number is incremented for
each survey station where a reading was taken.
Spreadsheet columns can contain different types of information. The second and third
columns contain location information. For IP data, these columns contain line and
station numbers. Data in columns of a spreadsheet are usually referred to as channels.
A channel is a column of data that contains measurements or values for each location
along a survey line.
Profiles are displayed in the Profile pane for the channel selected to be profiled. In the
sample image above, the Profile pane displays the IP array channel in red. The
selected data point or cell in the spreadsheet is indicated by a square in the Profile
pane. The numbers at the left indicate the range of data values. By changing these
numbers, you can modify the scale of the profile pane and alter the appearance of the
profile lines.
Array Channels
The first thing you will notice about an array channel is that each cell in the
spreadsheet column contains a curve instead of a number. The reason for this is that
an array channel contains more than one channel or column of data. For example, in
an IP time domain survey, a curve in a spreadsheet cell of an array channel would
represent a decay curve for a single survey location.
By representing data in an array channel, all the readings for a single location can be
put into one column of the spreadsheet instead of having several channels for the
multiple readings at each survey location. The array channels do not display the data
numerically, but represent the data as a curve. The numerical data stored in an array
channel may be displayed in sub-channels. For example, if the decay curve contains
256 time slices, you have the option of looking at any slice of that curve in its own
sub-channel.
Important! As with all cells in the database, when a cell in an array channel is
selected (highlighted), pressing a key overwrites the data.
10 Chapter 1: System Capabilities and IP Concepts
You can create a mask channel for any single array channel or a pair thereof. This
enables you to turn “off” individual data window values in the Average duplicate
samples dialog (IPDUP GX). The dialog creates a “mirror” array channel of type
“byte”, with the same number of columns as the selected channel. By default, all
these values are “1”, which indicates that all the values are selected to be used. IP QC
Tool recognizes that an array channel has been assigned a mask channel, and checks
the Mask check box. You can edit the database array channel mask values using
virtual channels. For instance, if the array mask channel is named IP_MASK, typing
IP_MASK[2] in the channel header would display column 2 (numbered from 0) of
the IP_MASK array mask channel, and then the individual mask values can be edited.
Two channels can be selected to use the same mask channel. This is useful, for
instance, for frequency domain systems where channels typically occur in pairs for
In-Phase and Quadrature, or Amplitude and Phase. Normally, if one of these values is
bad, the corresponding value in the other channel is also bad, and both will be
disabled at the same time. At present, the array mask channels are only used in the
Average duplicate samples dialog. They have no effect on plotting or windowing
functions. The QC channel still defines whether any single row of data is included.
Database Channels
The table below provides a list of the IP database channels you will find in the
columns of the Spreadsheet pane. The channels in this list may vary depending on the
instrument used to collect the data. You may not see some of these channels, or may
have imported other channels not listed here.
Channel Description
X Line number Note: This example assumes an E-W line direction. For a
database with an N-S line direction, these channels
Y Station
are swapped.
Z Estimated depth (plotted depth).
IP Time window (an array channel). For information on array channels, see
Array Channels on page 9. Note that the scale of the curve in each array
channel cell is unique.
IP_Avg IP Average (calculated by the Induced Polarization system during import).
N The “n” value, which is not defined for pole-pole or gradient surveys
(calculated by the Induced Polarization system during import).
I Transmitted Current.
MF Metal Factor.
QC Flag created using the IP|Quality Control The flag values are:
1 Accept: The value is used in calculations and plotted.
2 Flag: The value is not used in calculations, exports, and averaging, but
is plotted with square brackets on pseudo-sections if no “better” values
exist for the same station.
Chapter 1: System Capabilities and IP Concepts 11
The screenshot below depicts the IP Quality Control dialog – the GUI element
through which you perform the IP QC operations.
The data are represented graphically in right-hand section. The above example shows
data from a magnitude array channel in a frequency domain survey dataset from the
CRPD132.FLD sample file.
The following table describes the elements of the IP Quality Control dialog.
12 Chapter 1: System Capabilities and IP Concepts
Element Description
Select a regular or array channel from the drop-down list.
Data are displayed for a single line. Select a line from the
drop-down list.
For the given line and channel, the data are sorted, first by the
station value (taken from the “Stn” channel) and then by the
N value (taken from the “N” channel). If the N value is
undefined, the data are sorted by the station number only.
Use the arrow keys to move through the data, one (Station, N)
value at a time.
The double arrow buttons go to the first and last (Station, N)
values.
Only “raw” data values (the value in the Type channel is 0)
are displayed. “Averaged” data values (the value in the Type
channel is 1) are not displayed. However, the average of the
currently displayed channel with the Accept status is
displayed in green in the data section.
Displays the number of rows of data sharing the current
(Station, N) values, and the number of the sample currently
highlighted in red in the data section.
Use the arrow buttons to step through the samples.
The Acceptance option buttons are linked to the QC channel
values as follows:
• Accept – QC = 1 – Values are used in all calculations and
export operations.
• Flag – QC = 2 – Values are ignored in calculations and
export operations but are plotted as labels on pseudo-
sections, enclosed in square brackets (if there are no
“accepted” values at their (Station, N) position).
• Reject – QC = * (dummy) – Values are always ignored in
plotting, averaging, and export.
Chapter 1: System Capabilities and IP Concepts 13
Data value ranges These values, plotted to the left of the data section beside the
top and bottom edges, display the current range of data values
plotted. These are for display purposes only and are not
... editable.
The time domain (transient) IP method identifies areas in the ground where the
voltage decay takes a noticeable period of time after the current is removed. This
technique measures the bulk or average polarization of a volume of rock. Polarization
is the ability of the ground to store electrical energy, similar to a galvanic battery.
The frequency domain (or variable frequency) method uses the magnitude and
phase shift of the frequency to calculate the apparent resistivity, which is the bulk or
average resistivity of a volume of rock. The resistivity value of a rock mainly depends
on the porosity and salinity (that are translated to resistivity) of the solution filling the
pore spaces. The method is used to search for areas where the resistivity of the rocks
decreases with the increase of the frequency of the applied current.
The relationship between the frequency domain and the time domain methods is
defined by the following equation (the LaPlace Transform theory):
fe
m= ≅ fe ,
I + fe
Where the Chargeability parameter (m) measured in the time domain method is
exactly equivalent to the frequency effect parameter (fe) used in the frequency
domain method, and the current is represented by I.
IP surveys require two pairs of electrodes to introduce electric current into the ground
and then measure the response voltage. The current electrode pair introduces current
into the ground. The potential electrode pair measures the potential difference at
regularly spaced positions along each survey line. The distance between the potential
electrodes is generally represented by “a”. As this distance increases, the survey
depth also increases, meaning that a greater volume of rock is sampled. The “n” value
is an integer that refers to the distance and order of the potential electrodes. The n
value increases with the distance of the potential electrodes from the current or
transmitting electrode. The greater the n value, the deeper the penetration and the
larger the sample location (in general).
16 Chapter 1: System Capabilities and IP Concepts
In most array configurations, the IP and resistivity results are plotted as profiles along
a measurement line. For each survey line, a mathematical model is used to match
each reading to a theoretical position below the surface. The data (actual readings)
represent the true path of the current projected onto the surface. The model projects
these data to a probable path and assigns a new location for each data point. The
result is a two-dimensional vertical slice of data for each survey line, referred to as a
“pseudo-section”. The vertical scale increases downwards, and represents the “n”
value – a rough measure of depth.
The resulting data points describe the relative positioning of anomalies inherent in the
data. A sense of depth to individual anomalies is conferred by the vertical axis. This,
however, is not a true measure of depth. Data inversion is necessary to obtain a more
realistic depth measurement.
Field surveys use the different electrode arrays. The montaj Induced Polarization
system supports four array types for both time domain and frequency domain
surveys:
• Dipole-Dipole
• Pole-Dipole
• Pole-Pole
• Gradient
Chapter 2: Quick Start Tutorials 17
The system enables you to access files anywhere but it is a good strategy to carefully
organize your data (project information and files) before carrying out any processing.
To start this tutorial, create a working directory called D:\Tutorial. A general rule to
follow while working with Geosoft applications is to avoid working in the Geosoft
directory. In the tutorials, we follow this rule by keeping all the working data, found
in C:\Program Files\Geosoft\Oasis montaj\data\ip\, in your working directory
D:\Tutorial.
The time domain examples are based on dipole-dipole survey data in the 1P6.DMP
file. This file is an instrument dump file from an Iris Elrec-6 IP instrument.
All the frequency domain examples in this tutorial are based on the data from a
pole-dipole CRPD132.FLD data file (in the Zonge FLD format).
Other sample data from different instruments are also included. You can use the
sample files to try the various processing scenarios with the system. A description of
each sample data file is provided in the Appendix: Sample Data. The descriptions
include the instrument/file type, survey domain, and array configuration for each file.
Chapter 2: Quick Start Tutorials 19
The project also controls your working directory. Projects are saved as (*.gpf) files. If
you open an existing project from a directory, the system assumes that all your
project files are located in the same directory. To streamline your work, as well as to
keep it organized, make sure that your project file is in the same directory as the other
files you want to use. We recommend that each project you work on have its own
project (*.gpf) file. If you use a number of applications or add-on tools in Oasis
montaj that have different menus, you can use the project to display only the menus
you require.
The Project Explorer tool enables you to browse and open project items. The Project
Explorer pane has two tabbed sections. The Data section displays all data files
included in the project, and the Tools section organizes and maintains the project
tools. To access the Tools section, you click the Tools bar at the bottom of the Project
Explorer pane. To return to the Data pane, you click the Data bar at the top of the
Project Explorer pane.
TO CREATE A PROJECT :
Note: Oasis montaj assumes that your data is in the directory containing this project
(D:\Tutorial).
3. Specify a name and directory for the project. For example, name the project IP
and place it in the working directory D:\Tutorial.
4. Click the [Save] button.
The system saves the project and indicates that it is open by adding menus to the
menu bar, adding buttons to the toolbar, and by displaying the Project Explorer
pane. These are visual clues indicating that you are ready to start working with
the system.
Before you can start working with the montaj Induced Polarization system, you
have to load the IP menu in your project. If you require more detailed information on
setting menus, refer to the Oasis montaj Online Help System (Help >Help Topics).
1. From the GX menu, select Load Menu or click the Load Menu icon ( ) on the
toolbar.
Chapter 2: Quick Start Tutorials 21
The system displays the IP menu on the main Oasis montaj menu bar.
5. Click [Next].
Chapter 2: Quick Start Tutorials 23
6. Specify Windows for Averaged Total IP, Resistively Units, Metal factor formula,
Metal factor multiplier, Minimum N spacing, Maximum N spacing, N spacing
increment, Pseudo-section filter, Filter weights and Adjust app. res. for
topography? defaults.
7. Click [Next].
8. Create a “boilerplate” text that is plotted on all maps you produce using the
montaj Induced Polarization system.
9. Click [Finish].
24 Chapter 2: Quick Start Tutorials
Before you can process data, you must create a database and import the data into it.
The montaj Induced Polarization system can import data from instrument dump
files directly into an Oasis montaj database. In this example, we import an
instrument dump file from an Iris Elrec-6 IP instrument (1P6.DMP).
The new database (iptime) is displayed in the background. The Import IP Data
dialog is displayed.
Chapter 2: Quick Start Tutorials 25
4. From the Import from drop-down list, select the format of the file to import (in
this example, “Instrument Dump File”).
5. From the Import mode drop-down list, select the import mode (in this example,
“Overwrite”).
Note: If you select "Overwrite", all existing data are deleted from the database
before the import is performed. If you select "Append", the newly imported
data are appended to the existing data, line by line.
6. Click [OK].
7. Using the […] button, locate the Data file to import (IP6.DMP).
Note: You may need to change the Files of Type in the Data File dialog to (*.dmp),
then choose the file to import (1P6.DMP) and click [Open].
8. From the IP Instrument drop-down list, select the instrument as “Iris Elrec-6”.
9. Click [OK].
The Iris Elrec-6 Parameters dialog is displayed. Note that the dialog varies based
on the selected IP instrument.
The data are imported into the database and displayed in the spreadsheet pane.
After importing data into a database, it is always a good idea to save your changes to
avoid having to import the data again if something goes wrong.
1. From the Data menu, select Save Database Changes or click the Save Database
2. Click [Yes].
Note: To revert back to the last saved version of the database, you can use the Data
> Discard Database Changes menu option or click the button on the
main toolbar.
Chapter 2: Quick Start Tutorials 27
IP Quality Control Tool provides an interface for quality control operations on “raw”
samples in the current IP database. With this tool, you can evaluate duplicate samples
by navigating through the data by station and line number.
The data in the database are not changed during the quality control (QC) process.
Instead, the process sets the value of the QC channel to a number that indicates which
data points to include when plotting, exporting, and averaging duplicate samples.
Quality control can be performed on both single-value and array channels.
The following procedure guides you through the process of evaluating duplicate
samples for the IP array channel of a time domain database. The procedure uses the
iptime.gdb database you created earlier, in the Importing Time Domain Data
section.
Since the IP channel is an array channel, the data appears in the data pane as a
connected series of data values.
Each data point represents the value of the data for a specific time window. The
connecting line shows the signal decay over time. For details, see the IP Quality
Control Tool section.
4. Select the different data values by clicking the squares representing the values in
the data pane. Optionally, modify the individual data values by entering a new
value for the selected data point in the Data Value field. For example, to change
data value for window 8 of 10 at station 262.5, first select the data point with the
cursor. In the Data Value field, replace the number 13.3 with 7 and press [Tab]
on the keyboard. A new data point is displayed with a value of 7.
5. To specify which data sample you would like to use for calculation, use the
Acceptance option buttons.
6. To switch between the duplicate samples, use the Sample arrow buttons. For
example, to reject the second sample for station 262.5, switch to sample 2 of 2
and select the Reject option button.
7. Click [OK].
Chapter
apter 2: Quick Start Tutorials 29
In the QC channel of the database, an asterisk (*) indicates a dummy value for the
second sample at station 262.5.
This sample will not be used in any calculations, exports, or plots. You can also
edit the acceptance status of the data points directly by changing the numbers in
the QC channel.
Use the windowing utility to create a new channel of IP values to store one or more of
the time windows in an array channel column. This utility is useful when you want to
plot a single channel (from an array) in a pseudo-section.
When plotting an array channel, the system averages the channel values and creates a
new channel with the averaged values.
Note: For time domain array data, multiple channels are summed using a weighted
average based on the width of each time window. For frequency domain data,
values are simply averaged.
TO WINDOW AN ARRAY CHANNEL :
2. From the Array Channel drop-down list, select the channel to window (“IP”).
3. In the New channel name field, type the name of the channel for the windowed
values as (“IP_Window”). Note that, if this channel exists, the existing data will
be overwritten.
4. In the Windows/Frequencies to include field, enter the range of the time windows
to be averaged in the new windowed channel as “0:9”.
Ranges are allowed, indicated with a "-" or ":", and individual windows are
delimited with commas or spaces. If the field is left blank, all the windows are
used. If there are N windows, the first window is "0" and the last is "N-1". To
specify multiple channels, use the following notations:
0-3,5-6,8,10 Channels 0, 1, 2, 3, 5, 6, 8, and 10
0:9 All windows from 0 to 9
0-9 All windows from 0 to 9
3- Windows 3 and up
-5 Windows up to and including 5
<Blank> All windows
5. Click [OK].
A new channel is created in the database containing the average of the windowed
values you specified.
You can display topographic information related to the survey readings on your maps.
The topographic elevations are displayed as offsets next to the survey readings.
Adding topographical information does not affect any of the data in the IP database or
any derived calculations – the elevations are for display only. In the montaj Induced
Polarization system, topographical information can be imported either from a grid or
from an ASCII file. If no topography information is defined or imported, the
topography for all survey points is assumed to be 0.
Chapter 2: Quick Start Tutorials 31
The surface elevation for each survey point is imported into a new database
channel called Topo. Values in the Z channel are recalculated as follows:
⎛ n×a⎞
z = Topo − ⎜ ⎟
⎝ 2 ⎠
4. Save changes you made to the database. For procedure, see the Saving Database
Changes section.
TO IMPORT TOPOGRAPHIC INFORMATION FROM AN ASCII FILE :
1. From the IP menu, select Import > Topography > ASCII File.
2. Using the […] button, select the ASCII file. The selected ASCII (CSV) file must
contain topographic information with either commas or spaces separating values.
It must have three columns of values.
3. Click [OK].
32 Chapter 2: Quick Start Tutorials
The surface elevation for each survey point is imported into a new database
channel called Topo. Values in the Z channel are recalculated as follows:
⎛ n×a⎞
z = Topo − ⎜ ⎟
⎝ 2 ⎠
4. Save changes you made to the database. For procedure, see the Saving Database
Changes section.
If you are ready to import data into a new database, you can close your current
database.
2. Click [Yes].
Before you can process data, you must create a database and import the data into it.
The montaj Induced Polarization processing system can import data from
instrument dump files directly into an Oasis montaj database. All the frequency
domain examples in this tutorial are based on data from the CRPD132.FLD pole-
dipole data file (in the Zonge FLD format).
4. From the Import from drop-down list, select the format of the file you want to
import as “Zonge FLD Format”.
5. From the Import mode drop-down list, select “Overwrite”.
Note: If you select "Overwrite", all existing data are deleted from the database
before the import is performed. If you select "Append", the newly imported
data are appended to the existing data, line by line.
6. Click [OK].
7. Use the […] button to specify the Zonge FLD file to import as “CRPD132.FLD”.
8. From the Station Scaling drop-down list, select the scaling method (“None (1 to
1)”).
9. Click [OK].
The data are imported into the database and displayed in the spreadsheet pane. If
you are importing frequency data with more than one base frequency, an
amplitude and phase channel sre created in the database for each frequency. For
example, in the CRPD132.FLD data file, readings were taken using 3 base
frequencies: 0.125 Hz, 0.016Hz, and 0.063 Hz. From these frequencies, 6
channels are created in the database: F125Amp, F125Phs, F63Amp, F63Phs,
F16Amp, and F16Phs.
Depending on the survey instrument, the survey data might also have additional
specialized channels not listed in the Database Channels table, on page 10.
Chapter 2: Quick Start Tutorials 35
After importing data into a database, it is always a good idea to save your changes to
avoid having to import the data again if something goes wrong.
1. From the Data menu, select Save Database Changes or click the Save Database
2. Click [Yes].
Note: To revert back to the last saved version of the database, you can use the Data
> Discard Database Changes menu option or click the button on the
main toolbar.
This utility enables you to edit distant electrode positions for a selected line in the
database.
1. To select a line to edit, right-click the top left cell in the spreadsheet, select the
List option from the pop-up menu, and select a line from the list.
2. From the IP menu, select Distant Electrodes.
36 Chapter 2: Quick Start Tutorials
Two upper (read-only) fields display the current survey line and array type.
3. Specify the following settings (as appropriate for your array type):
• Current Electrode X (T1X, T1Y, T1Z) – Near current electrode positions.
These are requested for Gradient surveys only. If X or Y is left blank, the
electrode is considered to be at infinity for the purpose of calculating apparent
resistivity. If the Z value is left blank, no topographic correction can be
applied for calculated resistivities.
• Current Electrode Y (T2X, T2Y, T2Z) – Distant current electrode positions.
These are requested for Pole-Dipole, Pole-pole, and Gradient surveys. If X or
Y is left blank, the electrode is considered to be at infinity for the purpose of
calculating apparent resistivity. If the Z value is left blank, no topographic
correction can be applied for calculated resistivities.
• Potential Electrode Z (R2X, R2Y, R2Z) – Distant potential electrode
positions. These are requested for Pole-Pole surveys. If X or Y is left blank,
the electrode is considered to be at infinity for the purpose of calculating
apparent resistivity. If the Z value is left blank, no topographic correction can
be applied for calculated resistivities.
4. Optionally, click [Export to lines]. This creates channels for the distant electrode
locations (TX2, TY2, etc.), and fills them with the existing values. From that
point on, the distant electrode values are set in (and taken from) these channels.
This makes it possible to have different distant electrodes for each line.
5. Click [OK].
IP Quality Control Tool provides an interface for quality control operations on “raw”
samples in the current IP database. Using this tool, you can evaluate duplicate
samples by navigating through the data by station and line number.
Chapter 2: Quick Start Tutorials 37
The data in the database are not changed during the quality control (QC) process.
Instead, the process sets the value of the QC channel to a number that indicates which
data points to include when plotting, exporting, and averaging duplicate samples.
Quality control can be performed on both single-value and array channels.
The following procedure guides you through the process of evaluating duplicate
samples in a frequency domain database. The procedure uses the ipfreq.gdb database
you created earlier, in the Importing Frequency Domain Data section.
The crosshair cursor in the data section indicates the currently selected data
point. The green line showing the average of all the samples is adjusted to fit the
remaining samples. For details, see the IP Quality Control Tool section.
38 Chapter 2: Quick Start Tutorials
3. Select the required data value or sample by clicking the left mouse button on
another data point. You can also move through the samples by using the Sample
arrow buttons.
4. To reject a selected sample, select the Reject option button in the Acceptance box.
5. Select another sample. Note that the rejected sample changes to grey to indicate
that it will no longer be used in calculations.
6. Use the Channel, Line, and Station/Transmitter parameters to select other samples
to evaluate.
7. When done, click [OK].
Note: The samples you rejected are indicated with an asterisk in the QC channel of
the database. These samples will not be used in calculations, exports, or plots.
You can also edit the acceptance status of the data points directly in the
spreadsheet pane by changing the numbers in the QC channel.
Samples are considered duplicates when two or more rows of data exist for the same
station and N value. If N is a dummy value, all rows of data with the same station are
considered duplicates.
The Average Duplicate Samples utility averages every channel in the database except
for the QC and Type channels. The averaged data are used for plotting and exporting
data. If the duplicate samples in a database have not been averaged, the last raw
number is used when exporting or creating pseudo-sections.
The Type channel is used to record whether a sample is raw data (0) or averaged data
(1). Only raw data are used in the duplicate averaging process. In addition, only those
samples whose value is 1 in the QC channel are included in the average.
Averages are performed on all data that share the same X, Y, and N (if provided)
values.
The Average Duplicate Samples utility (the IPDUP GX) runs the IPRECALC GX at
its conclusion, so as not to average values such as those found in the ResCalc channel
but calculate them directly from the averaged I (current) and Vp (initial voltage)
values.
Chapter 2: Quick Start Tutorials 39
The montaj Induced Polarization system creates three types of plots: pseudo-
section plots, stacked section maps, and 3D views of stacked section maps.
In most array configurations, the IP and resistivity results are plotted as profiles along
a measurement line. For each survey line, a mathematical model is used to match
each reading to a theoretical position below the surface. The data (actual readings)
represent the true path of the current projected onto the surface. The model projects
these data to a probable path and assigns a new location for each data point. The
result is a two-dimensional vertical slice of data for each survey line, referred to as a
“pseudo-section”. The vertical scale increases downwards, and represents the “n”
value – a rough measure of depth.
The resulting data points describe the relative positioning of anomalies inherent in the
data. A sense of depth to individual anomalies is conferred by the vertical axis.
The following picture presents a sample pseudo-section with explanations of its main
elements.
Channel Cell
Profiles showing
channel values
along the
section line
Legend
Pseudo-sections by channel
Chapter 2: Quick Start Tutorials 41
Pseudo-section map names are constructed from the line number and the channel
specified in the control file, e.g., L1300_IP_Avg.map.
3. If the example IP Plot Control file (psplot.ipi) is not selected, use the […] button
to locate the file in your working directory.
4. In the Map Name Tag field, type “Tutorial”. The map name tag is used to avoid
overwriting previous maps for different pseudo-section plots created from the
same line(s). If no map name tag is specified, the map name is created based on
the line name; e.g., “L1300.map”.
5. From the Plot drop-down list, select the lines to plot as “Displayed line”. A new
map is created for each line that you select.
6. Click the [Modify] button to change the options in the IP Plot Control file.
Note: If you click the [OK] button, the pseudo-section is plotted immediately with
the parameters defined in the IP Plot Control file.
42 Chapter 2: Quick Start Tutorials
The Section Layout page is used to set general pseudo-section plot parameters,
such as Display, Station Scale Line, Label Sizes, and Line and Station Labels.
7. In the Display section, check the Plot colour bars check-box. A vertical color bar
is plotted, based on the color zones detected in the plotted pseudo-section grids.
Note: Space must be defined for the Plot Right Channel Label Space (cm) in the
Horizontal space for labels section on the Page Layout page.
Chapter 2: Quick Start Tutorials 43
8. Leave the default values in the rest of the fields and select the Data tab.
10. Define/edit the section parameters as required. For parameter descriptions, click
the [Help] button in the dialog. For this tutorial, leave the default values for all
parameters.
11. Click [OK].
In this dialog, you can opt to plot the selected channel’s data as a profile using the
Filter, N Value, or SP Profile option. For parameter descriptions, click the [Help]
button in the dialog.
13. Define/edit the profile parameters as required. For this tutorial, leave the default
values for all parameters.
14. Click [OK].
This page contains those parameters that determine how the pseudo-sections are
plotted on the map, such as Plot Size, Margins, and Horizontal space for labels.
Most of these parameters use intelligent defaults, so you should not modify them
unless necessary. For example, the scale is determined automatically when
plotting.
Chapter 2: Quick Start Tutorials 47
On this page, you define parameters for the profiles plotted in a graph above the
pseudo-sections. A profile window is only plotted if a specific data profile has
been assigned to it (see the Data and Additional Profiles pages).
17. For this tutorial, leave the default values in all the fields.
48 Chapter 2: Quick Start Tutorials
19. Define contouring options that apply to all sections. Contour levels for individual
sections are defined on the Data page by clicking the [Define] buttons.
Chapter 2: Quick Start Tutorials 49
21. Check the Plot Legend check box to enable parameter editing.
22. Using the [Browse] button, select the MAPPLOT Control File as “C:\Program
Files\Geosoft\Oasis montaj\etc”. This file contains commands for plotting a
legend, boilerplate, date, company logo, etc.
50 Chapter 2: Quick Start Tutorials
On the Additional Profiles page, you can select channels to plot as profiles even if
you are not plotting these channels as pseudo-sections. You can also select grids
as profile data sources. The grid data for the profile are extracted from the grid
along the line path.
24. To define an additional profile:
• From the Channel or Grid drop-down list, select the required channel/grid.
• Click the adjacent [Define] button.
• In the <Name> channel profile parameters or <Name> grid profile
parameters dialog that is displayed, define the additional profile parameters.
For this tutorial, do not define additional profiles.
25. Click [OK].
Chapter 2: Quick Start Tutorials 51
The montaj Induced Polarization system enables you to plot “stacked sections” for
a specific. The stacked section map shows a pseudo-section on each survey line for a
specified channel. The distances between the lines are displayed accurately to scale.
The idea of the map is to provide a pseudo cross-section of the survey area.
The stacked section map names are constructed from the line number and the channel
specified in the control file, for example, L1300_ResCalc.map.
5. Using the […] button, select the IP Plot Control file as “ipstak.ipi”.
The ipstak.ipi control file contains the settings that describe how the finished
map will look.
6. To modify the control file settings, click the [Modify] button.
You can use an existing IP Plot Control file (e.g., ipstak.ipi) or create a new
control file by clicking the [New] button.
7. Click [Next].
Chapter 2: Quick Start Tutorials 53
The Section parameters dialog contains the parameters that determine how the
sections are plotted on the map. Most of these parameters use intelligent defaults,
so you should not modify them unless necessary. For example, the scale is
determined automatically when plotting.
8. Click [Next].
54 Chapter 2: Quick Start Tutorials
11. Specify parameters that define how the contours are drawn on the sections.
13. Defines the various line attributes for the default contour and drawing lines.
15. Define text to be plotted in the title block of the map. The annotations are saved
in the mapplot control file. The default file for a staked section plot is psstak.con.
16. Click [Finish].
The requested stacked section map is created and displayed in your current
project.
The montaj Induced Polarization system enables you to create 3D views of your
stacked sections. This one-step process enables you to view the individual sections
oriented in real (3D) space. You initiate this process (the IP3DVIEW GX) by
selecting the Create 3D View in Stacked Section Map option from the IP menu.
The process was designed principally for the stacked section maps. It also works on
the single-line pseudo-section maps. However, because each of the sections in those
maps is at the same location, only one is plotted.
A new 3D View will be placed in the centre of the currently displayed map window.
You may move and scale the 3D View just as you would any other view on the map.
3. Click [OK].
4. Move, scale, and adjust the 3D view as required using 3D Tool. For procedure,
see the Online Help system and/or the Oasis montaj Quick Start Tutorials.
58 Chapter 2: Quick Start Tutorials
The IPEXPIPDATA GX is used to export one or more lines from the current database
to individual Geosoft IPDATA files.
In this section, we will use the ipfreq.gdb database. Make sure this database is open
and selected.
1. Make sure the current database contains the data you want to export.
2. From the IP menu, select Export > Geosoft IPDATA Format.
3. From the Channel drop-down list, select the regular or array channel you want to
export.
4. From the Export drop-down list, select the line or lines to export to individual
Geosoft IPDATA format files. The file names are created in the working directory
from the line numbers with the extension ".dat". One file is created for each line,
for example L1320.dat, L1330.dat, etc.
Note: Previously created files in the working directory with the same names will be
overwritten.
Chapter 2: Quick Start Tutorials 59
5. In the Title field, type title text to appear as the first line in the IPDATA files.
6. Click [OK].
You can export pseudo-sections created from the current database to Geosoft IPRED
format files. The IPRED format was introduced for the DOS IP system, and was
produced by running the IPRED program on IPDATA files.
The file names are created from the line numbers, with the extension appropriate to
the user-specified data type. One file is created for each line.
Note: Previously created files in the working directory with the same names will be
overwritten.
In this section, we will use the ipfreq.gdb database. Make sure this database is open
and selected.
1. Make sure the current database contains the data you want to export.
2. From the IP menu, select Export > Geosoft IPRED Format.
3. From the Export drop-down list, select the line or lines to export to individual
Geosoft IPRED files. One file is created for each line.
4. In the Title field, type title text to appear as the first line in the IPRED files.
60 Chapter 2: Quick Start Tutorials
5. From the Channel drop-down list, select the regular or array channel you want to
export.
6. In the Filename suffix (type) field, type an extension for the resulting IPRED files.
7. From the Filter type drop-down list, select “Pant-Leg”.
8. In the Filter weights field, type the filter weights to apply to each N expansion
starting with n=1. Filter weights must be listed in the w1,w2,w3... format,
separated with spaces or commas. Weights for separations greater than the N for
the last specified weight are set to 1.0. If this field is left blank, all weights are set
to 1.0.
9. Specify positions for the first and last stations. If the first station is left blank, the
minimum will be used; if the last station is left blank, the maximum will be used.
If both are left blank, the entire range is used (default setting).
10. In the Maximum n spacing field, type the maximum array expansion factor to
include in a pseudo-section. The default is 10.
11. Click [OK].
This export option enables you to export lines from the current database to a Scintrex
Geophysical Data Format (SGDF) file.
In this section, we will use the ipfreq.gdb database. Make sure this database is open
and selected.
1. Make sure the current database contains the data you want to export.
2. From the IP menu, select Export > Scintrex Geophysical Data Format.
3. Use the […] button to specify the name of the export file you want to create.
Chapter 2: Quick Start Tutorials 61
4. From the Export drop-down list, select the line(s) you want to export.
5. From the IP/Amplitude Channel drop-down list, select a channel to export:
• For time domain data, select the channel containing the IP data
• For frequency domain data, select the amplitude channel
6. For frequency domain data, from the Phase Channel drop-down list, select the
phase channel. If this field is left blank, the phase is not included in the output.
For time domain, this field must be left blank.
7. Click [OK].
This option is used to export one or more lines from the current database to individual
Interpex I2X files. The file names are created from the line numbers with the
extension ".i2x". One file is created for each line.
Note: Previously created files in the working directory with the same names will be
overwritten.
In this section, we will use the ipfreq.gdb database. Make sure this database is open
and selected.
1. Make sure the current database contains the data you want to export.
2. From the IP menu, select Export > Interpex I2X Format.
4. Click [OK].
This option is used to export data in the format required by the UBCDCINV2D and
UBCIPINV2D inversion programs developed at the University of British Columbia.
In this section, you use the ipfreq.gdb database. Make sure this database is open and
selected.
1. Make sure the current database contains the data you want to export.
2. From the IP menu, select Export > UBCDCINV2D/IPINV2D Format.
3. From the Data type drop-down list, select “RES” for DCINV2D or “IP” for
IPINV2D.
4. For the “RES” data type, from the Voltage channel drop-down list, select the
channel that stores the initial voltage.
5. For the “RES” data type, from the Current channel drop-down list, select the
channel that stores the current values.
6. For the “IP” data type, from the IP channel drop-down list, select “IP”.
7. From the Line selection drop-down list, select “Selected”, “Displayed”, or “All”.
8. To change the output directory, use the […] button next to the Output directory
field.
9. To output an export control file for (INP) each output line, check the Create
control files check box.
10. To define RES-specific or IP-specific export parameters, click [More].
11. For the “RES” data type, on the RES page, define parameters for the DCINV2D
control file.
64 Chapter 2: Quick Start Tutorials
12. For the “IP” data type, select the IP tab and define parameters for the IPINV2D
control file on the IP page.
This export option enables you to subset your data based on the various criteria
including displayed channels, all channels, a list of channels, or a masking channel.
The process automatically creates and displays a new database that includes the the X
and Y channels, as well as the selected channels.
In this section, we will use the ipfreq.gdb database. Make sure this database is open
and selected.
1. Make sure the current database contains the data you want to export.
2. From the IP menu, select Export > New X,Y Database.
Chapter 2: Quick Start Tutorials 65
3. Use the […] button to specify the Output Database as “IP Data.gdb”.
4. From the Channels to output drop-down list, select an option to specify whether
you want to export all channels, displayed channels, or selected channels in a list.
5. From the Lines to output drop-down list, select the line(s) to export.
6. Optionally, from the Mask channel drop-down list, select a mask channel to
restrict the data you export. For procedure used for creating a mask channels, see
the Oasis montaj Quick Start Tutorials or the Oasis montaj Online Help system.
7. In the N value field, type the N value to use. Depending on the N value you enter,
the program will function as follows:
• If the field is left blank (default), all the data are copied to the new database.
• If an N value of 0 is entered, instead of copying rows directly to the new
database, the pseudo-section function is called to calculate the Fraser-filtered
values of the requested parameters. The filter and its weights are taken from
the IP Default settings.
8. Click [OK].
If you selected the “Selected channels” option, a dialog appears showing 5 fields
for channel selection channels.
9. Select the required channels.
10. Click the [More] button to select more channels. The X and Y channels are
automatically included and do not appear in the list.
11. Click [OK].
The Create a New Database dialog is displayed with the name of the database to
be created.
12. Optionally, change the size of the database.
The new database is created, and the selected data are copied to the new database.
Appendix: Sample IP Dump and Input Data Files 67
The IP data directory contains a collection of files from a variety of instruments and
data formats. Some formats can be imported directly into an Oasis montaj database.
Others require some preparation involving setting default parameters, because the
files do not contain all the information required during import. For the default setting
procedure, see the Setting the montaj Induced Polarization System Defaults section.
If you do not edit the system default settings, the default values will be used (Dipole-
dipole array, E-W line direction, distant electrode locations are assumed to be at
infinity).
Station numbering in Zonge files does not always correspond directly to station
position in meters. In many cases, however, the station numbers may be converted
directly into locations by multiplication with the A-spacing. For example, if A = 50,
and the stations are numbered in Zonge as 0, 1, 2, and 3, the actual station locations
will be 0, 50, 100, and 150. The conversion method is specified during the import
process. Failure to set the conversion correctly results in nonsensical "N" and
ResCalc (calculated resistivity) values in the database.
References
John S. Milsom, John Wiley & Sons, 2003, Metal factor. Field Geophysics, Third
Edition, p.122.
Yanzhang Luo, Guiqing Zhang, 1997, Theory and Application of Spectral Induced
Polarization. Society of Exploration, 171 pp.