ZondGM2D Manual
ZondGM2D Manual
ZondGM2D Manual
ZONDGM2D
PROGRAM FUNCTIONALITY __________________________________________________ 2
SYSTEM REQUIREMENTS ____________________________________________________ 4
PROGRAM INSTALLATION AND DEINSTALLATION _____________________________ 4
PROGRAM REGISTRATION ___________________________________________________ 5
DENSITY OF THE ROCKS _____________________________________________________ 5
MAGNETIC SUSCEPTIBILITY OF ROCKS AND ORES ____________________________ 6
DESCRIPTION OF THE PROGRAM FUNCTIONAL _______________________________ 8
The toolbar of the main window ______________________________________________________ 8
Functions of the main window menu __________________________________________________ 9
Hot keys _______________________________________________________________________ 14
Status panel ______________________________________________________________________ 14
Selection of medium model in inversion (modes of regular cells grid, polygonal and randomly
layered medium) __________________________________________________________________ 40
MODELING ________________________________________________________________ 43
Modes of polygonal modeling and modeling using arbitrary layered medium _______________ 45
PROGRAM FUNCTIONALITY
ZondGM2D is a computer program for 2D interpretation of profile and areal multilevel data obtained from magnetic and gravity surveys. User-friendly interface and rich options
for data presentation allows solving existing geological problems with maximum efficiency.
The traditional exploration method for iron-bearing formations is magnetic survey. It
studies magnetic fields of objects, which contain ferromagnetic minerals. Physical interrelation
between data measured at the surface and magnetic properties at depths allows assuming
presence of magnetic bodies.
In magnetic exploration, total magnetic field is measured. It is comprised of the Earth
normal field, anomalous fields caused by magnetized bodies and magnetic field variations,
mostly related to the solar activity. A useful component associated with a studying area is
anomalous field, which can be identified, taking into account the normal field and measuring
magnetic variations in the vicinity of a survey area.
The magnetic field on the Earth surface can be represented as a vector sum:
, where
and
magnetic
.
, which, to a first-order approximation, is the magnetic field of the
Earth, is a field of a homogeneous magnetized sphere. The difference between the dipole field
(calculated) and areal and normal field (measured by satellite) is a non-dipole part
of the
normal field. It is often named as a residual field or a field of continental anomalies (sizes of
2
these anomalies are commensurate with the sizes of the continents). Its maxima does not exceed
30% of a dipole field value.
values gradually increase from 33,000 nT at the equator to 68000 nT near the poles, a
vertical component of a normal field in the North Pole area reaches 60,000 nT, changes sign to a
negative when passing through the equator, and gradually decreases from 0at the equator to -68
000 nT at the south pole of the Earth . Horizontal component has the maximum near the equator
(33,000 nT) and decreases to 0 at the poles. A gradient of the normal magnetic field is about 5
nT per kilometer.
Sources of the anomalous magnetic field
surface of the Earth. Maximum depth of magnetic rocks is about 50 kilometers. At greater
depths, magnetic properties of rocks disappear due to high temperatures.
A value of the field
Anomaly, where the field
is about 10% of
contribution to the total field is less than 1%. The most important variations can be divided into
periodic (solar-diurnal and short-periodic variations) and aperiodic (bay disturbance and
magnetic storms).
All variations are induced by external (relative to the Earth) sources: interaction of
charged particles with the ionosphere; electric currents propagating in stripe areas of high
latitude at an altitude of 100-150 km, etc. In general, an accounting for magnetic variations is
important in magnetic surveys because they are cause significant distortions in the observed data.
Program ZondGM2D allows solving forward and inverse problems of magnetic and
gravimetric surveys (recovery of anomalous magnetic susceptibility, density and geometry of
magnetic objects).
In the program, magnetic susceptibility is set in SI system (n*10-5), density in g/sm3,
measured values in nanotesla or milligals.
To solve an inverse problem (inversion) various versions of deconvolution methods and
Newton's method with focusing regularization are used.
A W
T
WA C T RC m AT W T f C T RCm
where A is the matrix of partial derivatives of measured data with respect to model
parameters (Jacobian), C smoothing operator, W matrix of relative measurement errors, m the
3
model parameter vector, - regularizing parameter, f vector of residuals between the observed
and calculated values, R focusing operator.
During the development of inverse problem, special attention was given to a priori
information usage (weights of individual measurements, ranges of parameters).
ZondGM2D uses a simple and clear data file format.
Program allows importing and visualizing data using other methods which makes data
interpretation process more integrated.
ZondGM2D is an easy-to-use instrument for automatic and interactive multilevel data
interpretation of magnetic and gravity surveys, and can be used on IBM-PC compatible PCs with
Windows operating system.
SYSTEM REQUIREMENTS
ZondGM2D can be installed on a PC with OS Windows 98 and higher. Recommended
system parameters are processor P IV-2 GHz, memory 512 Mb, screen resolution 1024 X 768,
color mode True color (screen resolution change is not recommended while working with
data).
As far as the program is actively using the registry, it is recommended to launch it as
administrator (right click on program shortcut run as administrator), when using higher than
Windows XP systems.
PROGRAM REGISTRATION
For registration click Registration file item of the main menu of the program. When a dialog
appears, select registration file name, and save it. Created file is transmitted to specified in the
contract address. After that user receives unique password which depends on HDD serial
number. Input this password in Registration field. The second option is to use the program
with supplied SenseLock key inserted in a USB-port while working.
Rock
Density (g/sm3)
Oil
0,8 -1,0
Coal
1,0
Water
1,1 2
Soil
1,13 - 2,0
Sand
1,4 2
Clay
2 - 2,2
Sandstone
1,8 - 2,8
Limestone
2,3 - 3,0
Salt
2,1 - 2,4
Granite
2,4 - 2,9
Gneiss
2,6 - 2,9
Gabbro
2,8 - 3,1
Basalt
2,7 - 3,3
Peridotite
2,8 - 3,4
Copperpyrite
4,1 - 4,3
Magnetite, hematite
4,9 - 5,2
mine salt, marble, oil, ice, graphite, gold, silver, lead, copper, etc. are diamagnetic [Hmelevskoj,
1997].
Paramagnetic rocks have positive magnetic susceptibility with low values. The majority
of minerals, sedimentary, metamorphic, and volcanic rocks are paramagnetic.
Ferromagnetic minerals (e.g. magnetite, titaniferous magnetite, ilmenite, pyrrhotite) have
very large values (up to several millions of 10-5 SI units).
Magnetic susceptibility of the majority of rocks depends largely on the presence and
percentage of ferromagnetic minerals in their composition.
Table below gives values of some rock-forming minerals and rocks. It can be seen from
the table that ferromagnetic minerals are strongly magnetic. Ultrabasic and base rocks are
characterized by the highest magnetic susceptibility values among other volcanic rocks, acidic
rocks are magnetic and weakly magnetic. Magnetic susceptibility of metamorphic rocks is lower
than the one of volcanic rocks. Sedimentary rocks with the exception of some sandstones and
clays are almost non-magnetic.
10-5 (SI units)
Mineral, rock
range
average
Quartz
10
Calcite
7 12
Gypsum
12
Coal
25
Sphalerite
750
Hematite
500 - 50000
6000
Pyrrhotine
103-107
150000
Ilmenite
5*105 -5*106
106
Magnetite
106-107
5*106
Limestone
25 - 3500
300
Sandstone
0 - 20000
400
7
Gneiss
100 - 20000
Granite
0 - 40000
2000
Diabase
1000 - 15000
5000
Gabbro
Basalt
30 - 150000
Peridotite
Sedimentary (average)
0 - 5000
1000
Metamorphic (average)
0 - 75000
50000
60 - 120000
60000
8000
30000
To convert values in CGS system which is used in the program, magnetic susceptibility
in SI units is divided by 4.
Magnetic susceptibility of para- and ferromagnetic rocks decrease with increasing
temperature and almost disappears at Curie temperature which ranges from +400 to +700C for
different minerals. Maximum depth of magnetic survey investigation is approximately 25-50 km.
Temperature is higher than Curie point at greater depths and all occurring there rocks become
essentially non-magnetic.
Prospecting geological structures and ores with magnetic susceptibility
are embedded
by rocks with magnetic susceptibility . This is the reason why similarly to gravity survey
redundant or effective magnetic susceptibility is of interest. Value of can vary and be
positive or negative. Magnetic anomalies appear due to variation of from zero [Hmelevskoj,
1997].
File/Open file
File/Save file
File/Print Preview
File/Edit file
File/Recent
File/
File/English
File/Reg file
File/Exit
Exit program.
Options/Mesh constructor
Options/Project information
Options/Program setup
Options/Data editor
Options/
Modules/Downward
continuation
Options/Modules/3D
sections plot
Options/Modules/Graphics
map
Options/Modules/Geologica
l editor
Options/Inversion/Set
boundaries
Options/Inversion/Resolutio
n
Options/Inversion/reference
model
model cells.
Options/Inversion/Focused
filter
Options/Model/Block
section
10
Options/Model/Smooth
section
Options/Model/Contour
section
Options/Model/Infinite
borders
Options/Model/X=Z scale
Make scales along the vertical and horizontal axes the same.
Options/Model/Extend
bottom
Options/Model/Show
modeler tools
Options/Topography/Topo
coefficient 1
Options/Topography/Import
topography
Options/Topography/Remov Remove topography data (all data will be assigned zero height).
e topography
Options/Topography/Restor
e topography
Options/Topography/Edit
topography
Options/Topography/Denois
e topography
Options/Topography/Smoot
Smooth topography.
h topography
Options/Topography/Import
3D topography
Options/Topography/Spline
intermediate
Options/Topography/Absolu Turn the absolute height mode on (when the mode is turned off,
te coordinates
Options/Topography/Extend
ed nodes
Options/Boreholes/Load
Open and display a file with borehole log data and stratigraphic
borehole data
columns (more).
Options/Boreholes/
columns width
Options/Import/Export/Imp
/Import model/graph
Options/Import/Export
Remove from the project a plot derived from the imported data.
/Remove data
Options/Import/Export
/Save selection
Options/Import/Export
Open file with the selected cells, and paste it at the current cursor
/Load selection
position.
Options/Import/Export
/Extract 1d log
Options/Import/Export
/Load 1d log
Options/Import/Export
Save 2D model settings file that can later be used in others Zond
programs.
Options/Import/Export
programs.
Options/Import/Export
Save 1D model settings file that can later be used in others Zond
programs.
Options/Import/Export
programs.
background image
Options/Import/Export/Rem
w background image
Options/Import/Export/Resi
ze background image
Options/Import/Export/Exp
ort to Excel
Options/Import/Export/Exp
w model in Surfer
Buffer
Operations/Line operations
Operations/Subtract
field
Operations/Subtract median Subtract median value.
value
Operations/Subtract
value
Operations/ Subtract med in Subtract an average value in a current window (useful when
window
Operations/ Subtract field Set and subtract a specified value if the current value exceeds it.
from height
Operations/Reduce to pole
Operations/Reduce to plane
Operations/Denoise data
Operations/Return
original
Operations/Areal operations
Help/About
Help/Context
Run manual.
^ File of imported data should consist of two columns: a coordinate of measurement, measured
value. A vertical axis of the imported data is right axis of a profile graphs plan. An imported
model must be retained in a project file Zond 2D.Whileloading it, a new window containing a
model appears. Moving the cursor in a modeling frame, a position of an active cell is displayed
on the model of the imported file.
Hot keys
Cursor pad /cursor in the model editor
Delete
Insert
Delete selection.
Space
Enter
Status panel
Status panel is divided into a few sections which contain different information:
Coordinates of cursor and active cell.
Parameters of active cell
Additional information.
14
WORKING GUIDELINES
Preparing and opening a data file
To begin working with ZONDGM2D, it is necessary to have a data file of specific
format which contains information about measurement positions and measured values.
One file usually contains data obtained from one profile. Text data files created in
ZONDGM2D format have *.GM extension. See data file format section for more details.
For correct work of the program data file must not contain:
unconventional separation symbols (use TABs and SPACEs only);
invalid values of measured data
It is desirable to not have more than 5000 observed data values in a single file.
16
17
During cursor movement in the model editor, cursor will be displayed in all other
sections, according to the size of the current cell.
18
Using the option Save/Load selection, it is possible to save or load a model tile. To save
the tile, it is necessary to switch on the model displaying with blocks mode (Blocks section) and
use Selection tool to extract a tile. Then, finally click Save selection.
A model tile can be loaded in the following way: Select a small area of the current model.
The upper-left corner of the selected area will be considered as the place from which the tile will
be placed. Call option Load selection and select a file name. If the selection is not available,
then the fragment will be pasted from the upper-left corner of the model.
In order to save or load a vertical parameter profile for a horizontal coordinate, use the
option Extract 1d log/Load 1d log. When saving the vertical profile it is necessary to set X
coordinate in the dialog box. When loading a vertical profile it is necessary to specify the X axis
range. This option can be used, for example, for incorporation of logging data, or for inspection
of profile intersections.
If a priori information is available, the option Import/export/Background image allows
using it (as an underlying image for the model editor). This can, for instance, be a geological,
electrical or seismic sections, or a section from an adjacent profile.
The image can be loaded from a *.sec file in scale.
* .sec file has the following structure:
1st line image file name;
2d line four coordinates X1 Y1 X2 Y2 of top left and bottom right corners of the image
are given with space separators
sect.emf
0 0 152.4 53.3
There is also an ability to load any image file. In this case, a dialog Set rectangle
appears, allowing you to set the boundaries of the imported graphic file.
19
In the Blocks section mode, cells with values different from the background medium will
be shown. Thus, it is possible to model anomalous objects on top of an underlying image.
In the Smooth section mode, colors of an underlying image and the current model will be
mixed and it is possible to see features of two sections simultaneously.
20
, if Automatic option is
ON. Otherwise the same image as it is displayed on the screen will be saved.
The option Export to Excel allows saving the current model parameters in a file with
Microsoft Excel format. Data in the obtained file are organized by rows, the first two rows are
distance and height in meters, other rows contain parameter data depending on a model cell.
The option Export model to DXF is intended for exporting a model to an AutoCAD
compatible file format.
item. In the dialog box, it is also possible to select the data format for saving calculated data for
the current model values of magnetic and gravity fields, and also image (Model, WorkSheet) in
*.BMP format.
Zond project data
calculated data
Worksheet
Model
Grid file
Section file
or respective
menu item. Upon successful file loading, the window measurement parameters selection appears
and dialog of initial model settings, in which it is possible to select survey layout parameters.
This dialog is also available in the Options item of the main menu. Window of measurement
parameters selection is divided into two parts: Magnetic survey and Gravity survey,
corresponding to settings for the magnetic and gravity fields.
22
Anomalous magnetic field and directional cosines of the normal magnetic field (Cx, Cz)
are calculated using following formulae:
Shown below are isoline maps of declination, inclination and total magnetic field in the
year of 2005. The current values for specific latitude and longitude can be found on specialized
sites or using GIS packages.
23
Area X nodes contains options, which allow setting the horizontal model grid parameters.
Minimum set the minimum coordinate of the modelling area.
Maximum set the maximum coordinate of modelling area.
Nodes number if the option In points is OFF, it sets a number of uniformly distant
nodes of the horizontal grid; otherwise nodes are set at measurement points.
24
Regular mesh set a uniform partition of a grid between measurement points (if the
option in points is ON).
Area Z nodes includes options, which allow setting parameters of a model vertical grid
(active for the borehole measurements).
StartH set the first layer thickness. This value need to meet the required resolution.
Factor set a ratio between the adjacent layers thicknesses. Values of the parameter are
generally selected in the range from 1 to 2.
Nodes number sets a number of layers.
on the
main panel, and is intended to set parameters related to the solution of forward and inverse
problems.
25
A W WA C C m A W
T
As seen from the equation, inversion model contrast is not minimized. This algorithm
allows reaching minimum misfit values. It is recommended to use it in the initial interpretation,
in most cases.
Occam - inversion by least squares method with using a smoothing operator and
additional contrast minimization [Constable, 1987]. As a result of using this algorithm, we get
the more smooth parameters distribution.
26
A W
T
WA C T C m AT W T f C T Cm
A W
T
WA I m AT W T f
Careless using of this inversion modification may lead to getting unstable result or to
increasing root-mean-square deviation, that is algorithm divergence.
It is best to use the Marquardt method as refining (to reduce misfit), after carrying out
the inversion with algorithms Smoothness constrained or Occam.
Focused - inversion by least squares method using a smoothing operator, and additional
contrast focusing [Portniaguine, 2000]. As a result of this algorithm a piecewise smooth
parameters distribution is obtained, that is a model, which consists of blocks with a constant
parameter value.
27
A W
T
WA C T RC m AT W T f C T RCm
Using this inversion variant requires careful Threshhold parameter selection. This
parameter defines the threshold contrast value for adjacent cells. When reaching this value the
cell parameters are not averaged between each other (that is, it is believed that a boundary
between cells exists). The dependence of the degree (or weight) of two neighboring cells
averaging Ri from the contrast threshold e, and contrast between these cells relooks as follow:
Ri
e2
e 2 ri
Blocks - selection of separate area parameters differing in the parameter value. Areas
with the same parameter values are treated as united blocks.
The matrix equation for this inversion algorithm is identical to the Marquardt algorithm:
A W
T
WA I m AT W T f
It is better to use the algorithm after results from previous methods were obtained
(Focused is recommended). To this end, cells should be preliminarily merged using the option
Cell summarization. To selecting individual blocks manually should use the model editor and
extract specific areas with different parameters. Some blocks will be marked with border during
work with this dialog.
Smoothing factor - sets correlation between measurements misfit minimization and
model misfit (for example, smoothness). For data with high noise levels or to get more smooth
28
and stable distribution of parameters, relatively large values of smoothing parameter (0.5 2)
should be selected; for high quality data values : 0.005 - 0.01 can be used. As a rule, large values
of data misfit are obtained for large values of smoothing parameter. It is used in inversion
algorithms Occam and Focused.
The figures below show results of inversion with smoothing factor values 0.01 and 1. The
resulting misfit for the first model - 4.5% and for the second 6%.
29
The figures below show the results of the focusing inversion with Threshold parameters
0.1 and 10.
Robust weighting scheme - this option should be switched ON if data contains some
spikes related to the systematic measurement errors. If amount of invalid data is comparable with
amount of valid measurements, the algorithm may not produce feasible results.
The second inlay Options is used for adjustment of additional inversion parameters.
Area Model
Min susceptibility, Max susceptibility - sets the limits of model parameters variation for
inversion. The values of the magnetic susceptibility need to be specified in the GHS system
(n*10-5).
30
Min density, Max density - sets the limits of variation of the model parameters for
inversion. Density values need to be specified in g / cm3.
Smoothness ratio defines a ratio of smoothing in horizontal and vertical directions. For
horizontally layered structures use value <1, for vertically layered >1. Usually, the parameter
values vary between 0.2 and 1.
The figures below show the inversion results for the two values of the parameter
Smoothness ratio: 1 and 0.3.
Cell grouping - use this option for large models. The option activates a table that allows
to combine adjacent cells and to get smaller number of unknown inversion parameters. When the
option is used, number of cells for the forward problem remains the same, but number of cells
for the inversion is reduced. Ideally, number of unknown parameters should be close to the
number of data measurements.
The table contains three columns: in the first (Layer) number of initial model layer is
specified; in the third (ID) number of inversion layer grid is set; in the second (Width) it is
necessary to set number of cells (in horizontal direction) contained in each cell of the inversion
grid for a current layer. Inversion grid will be displayed in the model editor during its
configuration. Double-clicking the left mouse button on a cell in the column Width allows
merging cells in the horizontal direction for a current layer, where as clicking the right button
does the same for the current and all deeper layers.
Double-clicking the left mouse button on a cell in the ID column allows merging cells in
the vertical direction, and clicking the right button - for current and all deeper layers.
31
Below you can see examples of three inversion grids: in the first, inversion grid
corresponds to model, in the second, from the second layer, cells are merged in groups of two,
and in the third from the second layer on, cells are grouped by four cells.
Program ZondGM2D also allow working with log data. In this case the window looks as
follow:
Fig. 4 The main program window when working in the log data mode.
All editing options are also available in this mode.
Model editor
Model editor allows changing parameter values of individual model cells using the
mouse. To the right from the editing model area a color scale is located, which links a color to
the parameter value. To select the current value, click on the scale using the right mouse button,
and its value will be displayed below the color scale.
The size and position of the color scale can be changed by pulling the slider, with the left
or right mouse button pressed.
Working with model cells is similar to the pixels editing in a raster image editors. When
moving cursor in the model area, on the low status panel of the main program window
coordinates and parameters of the active cell in which the cursor is located are displayed. The
currently active cell is highlighted by a rectangle.
The selected or fixed cell is marked in the middle by point.
33
The model can be displayed with cells Options/ Model/Block-section (Fig.5), with a
smooth interpolated palette Options/Model/Smooth-section (fig.5B) and as the contour section
Options/Model/Contour-section (fig.5C).
Double-clicking the mouse in the area of the model editor context menu appears with the
following options:
34
objects Indicates
whether
to
display objects
border
border.
Setup
Color bar
Zoom&Scroll
Print preview
Setup
Set minimum
Set maximum
Vertical axis
Log scale
Colors as histogram
Set maximum
Redivide
Thick mesh
Thin mesh
Horizontal axis
Redivide
the
unique
positions
of
measurement points.
Thick mesh
Remove
every
second
node
of
the
35
When right button clicking in the model editing area the context menu appears with the options:
Edit mode
Selection\Free form
Highlight set of cells within editing area using mouse. Area bounds is
selection
specified by a user.
Selection\Rectangular
Highlight set of cells within editing area using mouse. Field has
selection
rectangular form.
Selection\Elliptical
Highlight set of cells within edit area using mouse. Field has elliptic
selection
form.
Selection\Magic wand
Highlight set of cells within edit area using mouse. Active cell and
adjacent cells whose parameters are close to its parameter are
highlighted. Proximity is specified in model parameters setup dialog.
Selection\Remove
Delete selection.
selection
Mesh options\add
column /row
new boundary.
Mesh options\remove
column /row
Mesh options\resize
column /row
Clear model
Additionally, it is possible to use the option Show modeler tools in the menu Options/Model.
After selecting this option, the editing toolbar appears:
Switch on model editing mode.
36
Click on the cell with CTRL button pressed to drag selected set of cells within edit area
using mouse. While dragging the selection with left button pressed content of selected cells is
copied to new place. While dragging the selection with right button pressed content of selected
cells is moved to new domain.
37
Parameter alteration - determines value of augment to parameters of the selected cells (in
percentage terms relative to value of the parameter) when working in mode Edit with the bottom
Shift pressed.
Font button calls font settings dialog.
Colors tab
38
39
This option is more useful solving inverse problem (automatic rebuild density or magnetic
susceptibility sections). Result of inversion is a good initial point for further geological and
geophysical interpretation. Program options available in this mode are described in detail in the
previous sections.
40
41
window.
This model type is useful for data interpretation in regions with quasi-layered rocks
bedding. Details about the set of options are described in the part Modeling. Before starting to
work, it is necessary to set starting model parameters on the tab Modeling/Layered media
setup, where the following settings are available:
Layered numbers - number of layers
Boundary subnodes - number of nodes on each border
Build create the specified medium model.
Parameters (density, magnetic susceptibility and thickness) for selected number of layers
are set in the table.
42
The program solves forward and inverse problems for random layered medium and
allows getting densities and magnetic susceptibilities, and also geometry of layer boundaries.
When clicking button
, a pop-up
For this model type joint gravity and magnetic data fitting is available.
MODELING
Modeling is an important tool prior to the field work. It allows choosing optimal
parameters for measuring system in order to solve existing geological problems. Using initial
information about the survey target, interpreter can model different geological situations while
planning geophysical works.
43
On the first stage measuring system parameters should be selected taking into
consideration specifics and depth of burial of the object of interest, work conditions (relief,
grounding conditions), and equipment capability. Use File/Create synthetic survey option to run
measuring system parameters setup dialog:
Gravity and Magnetic options allow switching on or off calculation of the magnetic or
gravity fields.
Creating of a medium model is done in the model editor - bottom graphic part of the
program window during the switched on mode block-section. Read more in the chapter on
Model editor.
The main goals of mathematical modeling are assessment of signal level and assessment
of selected surveying system resolution. In order to calculate responses from a given medium
model (to solve forward problem), press the button
on the toolbar.
To recover the underlying model, that is, to solve the inverse problem for calculated from
the given model data, save theoretical signals by selecting a file filter - Zond project with
calculated data, and then open as observed. To compare restored and initial models use the
option Import model/data (part Importing and Exporting Data), after the initial model was
saved - Zond model with calculated.
window. In this mode, anomalous objects are replaced by a set of polyhedron with an infinite
spread in a direction perpendicular to profile. Each polyhedron has its own geometry, magnetic
susceptibility and density, for which forward problem is solved. Polygon modeling mode is used
in the second stage of data interpretation. After getting a model using one of the inversion
algorithms, it is recommended to switch over to polygonal modeling mode and define polygons
based on the obtained model and known geology, and perform inversion again.
Modeling mode using polygonal-layered medium is activated by pressing button
in
the main program menu. This mode is recommended to use for regions with concordant rocks
bedding. Before starting to work, it is necessary to specify initial model parameters in the tab
Modeling / Layered medium setup.
In order to do inversion in polygonal modeling mode and modeling using a horizontally
layered medium model use the button
Show background
Remove background
Remove polygons
Save polygons
Save polygons
Load polygons
Export to CAD
Display ColorScale
46
Polygons creating and editing is done using the control panel that appears when the
button
is clicked.
The mode Add polygon is called by clicking the button
mode is intended to create a local polygon. To add a new point to a polygon click the left mouse
button. To close the polygon (connecting the last polygon point to the first) click the right mouse
button. If it is impossible to create a local polygon (that is any faces intersect, or other object
exists in a polygon) program does not allow closing the polygon. When creating a polygon try to
not place points too close to each other.
The mode Delete polygon is called by clicking the button
mode is intended to remove a polygon. To delete polygon click the right mouse button on it.
Color of a polygon contour changes to red when the cursor is inside.
The mode Create a joined polygon is called by clicking the button
at the control
panel. This mode is intended to create a polygon (an additional part of a polygon), joined to the
existing polygons, model borders or relief, that is to any coupled model area. Coupled area
means polygons set and model borders with common faces. To add a new polygon point, click
the left mouse button. To close the polygon (connecting the last polygon point to the first, across
a joined area border) click the right mouse button. It is necessary that the first and last (closing)
polygon points are placed on the outside border (which changes color to red when the cursor is
47
approaching) of a coupled area. If it is impossible to create a local polygon (that is any faces
intersect, or other object exists in a polygon) program does not allow user to close the polygon
and remove all created points. Note that polygons joined to the left, right and bottom edges of the
model have infinite strike in those directions (that is they extend beyond the model).
The mode Disconnect polygon is called by clicking the button
This mode is intended to disconnect a polygon from a set of adjoint polygons or model edges.
Polygon disconnected from the model edges loses its infinite strike (it will be limited by the
model). To disconnect a polygon click the right mouse button on it. Color of the polygon contour
changes to red when the cursor is inside. Further, using the button Move polygon, it is possible
to move disconnected polygon part from the main polygon.
The mode Split polygon is used to create two parts inside a polygon. It is called by
clicking the button
on the control panel. The mode is intended to split the polygon into two
new connected polygons. The interface is defined by two points at borders or nodes of the
polygon, which is split. To select the first border point click the left mouse button. To select the
second point and to split polygon click the right mouse button. If the operation is impossible
(that is any faces intersect, or border is outside the polygon) program does not allow user to split
the polygon and remove the created border. Color of the borders and points of a polygon changes
to red when the cursor is approaching.
The mode Move polygon is called by clicking the button
on the control
panel. This mode is intended to move a polygon and all connected with it. To select a polygon
being moved, click the left mouse button; after that the connected area moves with the cursor. To
fix new position of polygons click the right mouse button. If the operation is impossible (that is
any borders intersect, or a polygon exists in other polygon) program does not allow user to move
48
the polygons and return them to original positions. Color of the polygon contour changes to red
when the cursor is inside.
The mode Add point is called by clicking the button
intended to add a new point on a border of an existing polygon. To add a polygon point, click the
right mouse button on its border. Color of the polygon borders change to red when the cursor is
approaching.
The mode Remove point is called by clicking the button
mode is intended to remove new point of existing polygon. To remove a polygon point, click the
right mouse button on it. The operation is impossible in the next cases: borders intersect, a
polygon is in other one or number of polygon points is less than three. Color of polygon points
changes to red when the cursor is approaching.
The mode Merge points is called by clicking the button
mode is intended to merge two points into one by joining points to the border of another
polygon, or to the model edges. Selection of the first merging point.
Brush calls polygon filling setup dialog. If this function is on, specified parameters are
used in all polygons of the model.
Gradient calls polygon gradient filling setup dialog. If this function is on, specified
parameters are used in all polygons of the model.
sets polygon density value.
sets polygon magnetic susceptibility value.
The following option specifies type of label displayed on a polygon. If this function is on,
specified type is used in all polygons of the model.
Value None there is no label on a polygon.
Value Density polygon density value is displayed on a polygon.
Value Susceptibility polygon magnetic susceptibility value is displayed on polygon.
Value Density & Susceptibility polygon density and magnetic susceptibility values are
displayed on a polygon.
Value User text value of the following field is displayed on polygon. The following option
specifies type of label displayed on a polygon. If this function is on, specified type is used in all
polygons of the model.
50
The area Window is used for window calculation parameters setup. Windows Minimum
and Maximum length are intended to set minimum and maximum width of the window,
respectively. The field Active allows selecting windows to be displayed on the model.
The area Index allows selecting value of structural index, which points will be displayed.
Structural indexes related to degree of field attenuation from a studied type source. For example,
for a linear mass N = 1, for points mass and dipole lines N = 2, for dipoles N = 3. The value N =
0.5 in magnetic survey occurs in a model with an inclined step model object, and also in gravity
survey for models consisting of thin layers, including dikes and sills. Points corresponding to
different indexes will be displayed with different icons.
Display with error show fitting error
Iteration number the number of iterations
Button Calculate to calculate position of singular points with selected parameters
Button Repaint update image
51
Downward continuation
Field downward continuation is a type of transformation intended for extraction of highfrequency and low-frequency field components. Downward continuation of the measured field Downward
continuation/Standard,
and
vertical
derivative
of
field
Downward
A)
B)
C)
Fig. 9 a) Plot of g, b) the result of downward continuation of the measured field, c) Downward
continuation of the vertical derivative.
52
When using the Blocks mode, depending on a contrast parameter (Contrast factor), cells
aggregation with similar parameters in area with a constant value are done. Start layers option
specifies a layer index from which the operation is applied.
In the Smooth mode, depending on a smoothing factor (Smooth factor), model cell
parameters are averaged. End layers option specifies the index of the layer up to which this
operation is applied.
Button
After getting one of models, click on the Buffer tab of the main menu and select Model 1.
The program offers to set name of a current model, in which different parameters can be stored.
Then, the current model will be saved, and a new inversion can be done with different program
settings and save as a Model 2. It is possible to save up to five different models. To open the
window of saved models, select Buffer/Open in the program main menu.
54
Geological texture selected by an interpreter is drawn on these local objects and layers. The
module allows also displaying borehole data, which simplifies the process of model construction.
Thus, the main goal of the module is a rapid construction of geological sections using
geophysical results and its export in CAD systems.
Before starting, it is necessary to select carefully a section type and its graphic settings.
The best variant is to display a section as isolines.
Next, Geological editor is run and actual work with a section begins. Bodies and the
layers are set as closed and unclosed polygons, color and texture are set to them.
Buttons for editing polygons and lines are located on the Geological editor toolbar.
The mode Add polygon is called by clicking the button
mode is intended to create a local polygon. To add a new point to a polygon, click the left mouse
button. To close the polygon (connecting the last polygon point to the first), click the right
mouse button. If it is impossible to create a local polygon (that is any faces intersect, or other
object exists in a polygon) program does not allow closing the polygon. When creating a
polygon try to not place points too close to each other.
The mode Delete polygon is called by clicking the button
mode is intended to remove a polygon. To delete polygon click the right mouse button on it.
Color of a polygon contour changes to red when the cursor is inside.
The mode Create a joined polygon is called by clicking the button
at the control
panel. This mode is intended to create a polygon (an additional part of a polygon), joined to the
existing polygons, model borders or relief, that is to any coupled model area. Coupled area
means polygons set and model borders with common faces. To add a new polygon point, click
the left mouse button. To close the polygon (connecting the last polygon point to the first, across
55
a joined area border), click the right mouse button. It is necessary that the first and last (closing)
polygon points are placed on the outside border (which changes color to red when the cursor is
approaching) of a coupled area. If it is impossible to create a local polygon (that is any faces
intersect, or other object exists in a polygon) program does not allow user to close the polygon
and remove all created points. Note that polygons joined to the left, right and bottom edges of the
model have infinite strike in those directions (that is they extend beyond the model).
The mode Disconnect polygon is called by clicking the button
This mode is intended to disconnect a polygon from a set of adjoint polygons or model edges.
Polygon disconnected from the model edges loses its infinite strike (it will be limited by the
model). To disconnect a polygon click the right mouse button on it. Color of the polygon contour
changes to red when the cursor is inside. Further, using the button Move polygon, it is possible
to move disconnected polygon part from the main polygon.
The mode Split polygon is intended to create two parts inside a polygon. It is called by
clicking the button
on the control panel. The mode is intended to split the polygon into two
new connected polygons. The interface is defined by two points at borders or nodes of the
polygon, which is split. To select the first border point, click the left mouse button. To select the
second point and to split polygon, click the right mouse button. If the operation is impossible
(that is any faces intersect, or border is outside the polygon) program does not allow user to split
the polygon and remove the created border. Color of the borders and points of a polygon changes
to red when the cursor is approaching.
The mode Move polygon is called by clicking the button
on the control
panel. This mode is intended to move a polygon and all connected with it. To select a polygon
being moved, click the left mouse button; after that the connected area moves with the cursor. To
fix new position of polygons click the right mouse button. If the operation is impossible (that is
56
any borders intersect, or a polygon exists in other polygon) program does not allow user to move
the polygons and return them to original positions. Color of the polygon contour changes to red
when the cursor is inside.
The mode Add point is called by clicking the button
mode is intended to add a new point on a border of an existing polygon. To add a polygon point,
click the right mouse button on its border. Color of the polygon borders change to red when the
cursor is approaching.
The mode Remove point is called by clicking the button
mode is intended to remove new point of existing polygon. To remove a polygon point, click the
right mouse button on it. The operation is impossible in the next cases: borders intersect, a
polygon is in other one or number of polygon points is less than three. Color of polygon points
changes to red when the cursor is approaching.
The mode Merge points is called by clicking the button
mode is intended to merge two points into one by joining points to the border of another
polygon, or to the model edges. Selection of the first merging point
The mode Disconnect points is called by clicking the button
mode is intended to move polygon point. To select moving point, click the left mouse button;
after that polygon point moves with the cursor. To fix position of the new point click the right
mouse button. If the operation is not possible (that is faces intersect) program does not allow you
to move a point and returns it to its original position. Points located on model boundaries are
moved only along appropriate edges. Color of polygon points changes to red when the cursor is
approaching.
Polygon parameters setup dialog
This dialog is used to adjust different polygon parameters. Double click on necessary
polygon to run it.
57
Add line
Move node
Remove node
Add node
Remove line
Create polygon of two lines
Create polygon of two lines
Move line
Save line
58
59
setting and selecting points along it, all included profile points will be displayed in blue. It is also
possible to exclude/include a point in the profile by pressing the left mouse button. If profile line
does not go through points, than position of measuring point projection on profile will be
displayed in green.
To display and edit coordinates of the measuring points, click the right mouse button on
an interesting point. The coordinates, which can be edited, will be displayed in the appeared
window show.
After profile setting it is necessary to click the button to call the data interpretation mode,
and then the main program window appears. To switch between profiles use Active line index
window
Option Graphics map (Fig.13) allows displaying graphs plan of the selected parameter
(magnetic or gravity field) for all profiles. To change scales of plots, click the left mouse button
in the graph area, and use the scroll wheel. It is also possible to switch on/off displaying of each
profile in the legend window on the right. To select one plot and switch off other, press Shift +
60
left-click on the legend. To move to other plots use scroll wheel. To setup the axes parameter,
click the right mouse button while pressing Shift (more).
62
Area Axis scales is intended to set axes scale. The scales are set only while the button
is pressed on the toolbar.
64
Windows of horizontal coordinate setting (distance from the profile start) are placed on
the main panel - horizontal position, and elevation - Position from surface. The horizontal and
vertical coordinates are set in meters. Empty square windows are used to create a fill set.
The module contains two main windows. On the left is Data window, containing a table
with the following columns: N - layer index, H - layer thickness in meters, Z layer bottom
depth in meters, C filling type. Borehole data are displaying on the right window in graphical
form.
To begin creating a lithology file click the button
appears in the Data Window. To set number of layers, click the button
. Then it is necessary
to edit the table: set thickness or bottom depth for every layer, and also choose the filling type
according to the lithology. To call filling dialog Pattern Color Editor, click double the left mouse
button in column C of the data window (Fig. 18). The program offers a rich variety of lithology
fills. The option Color allows choosing the fill color.
65
the graphics window. After that it is necessary to set the horizontal and vertical borehole
coordinates on the toolbar in kilometers, after that the borehole will be displayed in accordance
with its coordinates. In the graphics window active borehole is displayed in red.
For comfortable working with a large number of boreholes it is possible to create a
palette in the program. To create a palette, select a fill on the fills column of the Data window
and then click the right mouse button within fills area in the main program panel. In this way,
fills set can be created and saved. To do this, click the button
palette. Saved fills set can be used when creating a new data file of lithology and logging (
66
To load the available logging data of boreholes prepare a file format described in the part
"Logging and lithology data file format" and load it using the button
Enable/Disable editing
boundaries mode
Add new boundary
Delete boundary
Inclusion of the prior geological boundaries in the inverse problem is one of the most
important methods to improve quality of interpretation. On the one hand it increases the problem
67
stability, on the other - reduces the ambiguity and allows getting more plausible structure. In
those model areas where the parameters have very low sensitivity - it is practically the only way
to get acceptable result.
When boundaries exist, it is best to use the Occam algorithm. Normally 1-2 boundaries
are used and one should not forget that geological boundaries do not always coincide with the
petrophysical boundaries.
The boundaries on the resulting density section or section of the magnetic susceptibility
should be plotted taking into account the borehole data, or based on a priori information about
the structure of the studied area. To set boundaries click the left mouse button when the edit
boundaries mode is switched on. To close boundaries click the right mouse button. When
plotting boundaries do not use a lot of nodes. It is desirable that the boundaries are as smooth as
possible and go near the nodes of the inversion grid.
After boundaries are created, the inversion needs to be restarted, which will be calculated
with the actual boundaries (Fig. 20 B).
Fig.20 A section example of the inversion results without (A) and with (B) geological
boundaries.
mouse button with pressing Ctrl. To remove a slider, press the key Delete. It is also possible to
save a custom palette. To do this, click the button
Axes editor
Editor is used for adjusting graphic and scale axes parameters. Right click on necessary
axis with SHIFT button pressed to run it. Pop-up menu with two fields (Options and Default)
appears. The first one runs dialog, the second sets values on default.
First tab of Scales dialog contains options for axes scale parameters setup.
69
Option Auto defines how minimum and maximum axis values are chosen. If this option is ON
axis limits are set automatically. Otherwise Minimum and Maximum values specified by user are
selected.
Option Visible shows/hides selected axis.
Option Inverted defines axis orientation.
Button Increment change runs dialog for axis label step definition.
Option Logarithmic selects logarithmic or linear axis scale. In case of sign-changing scale
additionally use options from LinLog options field.
Option Base sets logarithm base for logarithmic axis.
Field LinLog options contains options for linear-logarithmic axis adjustment. Linearlogarithmic scale allows representing sign-changing or zero containing data in logarithmic scale.
Option Dec Shift sets indent (in logarithmic decades) relative to maximum axis limit modulo to
zero. Minimum decade (prezero) has linear scale, others have logarithmic.
Option Min dec sets and fixes minimum (prezero) decade value if option is ON.
Option Rounded limits defines whether it is necessary to round minimum and maximum axis
values or not.
Fields Minimum and Maximum contain options for axis limits adjustment.
Option Auto defines whether axis limit is selected automatically or using Change button.
Option Offset sets percentage axis limit shift relative to its actual value.
Tab Text:
Button Font calls header font setup dialog.
Button Outline calls dialog for header letters outline adjustment.
Option Inter-char spacing sets letter spacing in axis header.
Button Gradient calls gradient fill setup dialog for header text.
Option Outline gradient specifies if gradient fill is used in outline or interior of letters.
Button Shadow calls axis header shadow setup dialog.
70
Tab Text:
Button Font calls label font setup dialog.
Button Outline calls dialog for label letters outline adjustment.
Option Inter-char spacing sets letter spacing in label text.
Button Gradient calls label gradient fill setup dialog.
Option Outline gradient specifies whether gradient fill is used in outline or interior of letters.
Button Shadow calls label shadow setup dialog.
Tab Title contains options for axis header adjustment.
Tab Style:
Option Title sets axis header text.
Option Angle sets header text rotation angle.
Option Size sets header text indent. If 0 value is specified it is selected automatically.
Option Visible shows/hides axis header.
Tab Text:
Button Font calls header font setup dialog.
Button Outline calls dialog for header letters outline adjustment.
Option Inter-char spacing sets letter spacing in axis header.
Button Gradient calls gradient fill setup dialog for header text.
Option Outline gradient specifies if gradient fill is used in outline or interior of letters.
Button Shadow calls axis header shadow setup dialog.
71
Tab Text:
Button Font calls label font setup dialog.
Button Outline calls dialog for label letters outline adjustment.
Option Inter-char spacing sets letter spacing in label text.
Button Gradient calls label gradient fill setup dialog.
Option Outline gradient specifies whether gradient fill is used in outline or interior of letters.
Button Shadow calls label shadow setup dialog.
Button Minor calls internal intermediate axis tick line setup dialog. Option Len sets its length
Option Count sets number of intermediate ticks between main ones.
Tab Position defines axis size and position.
Option Position % sets axis indent relative to its standard position on graph (in percent to graph
size or in screen units depending on selected option Units).
Option Start % sets axis start indent relative to its standard position on graph (in percent to
graph size).
Option End % sets axis end indent relative to its standard position on graph (in percent to graph
size).
Graphs editor
Graphics editor is intended for graphic interface adjustment. Right click with SHIFT
button pressed on necessary graphic to run it. Pop-up menu with two fields (Options and
Default) appears. The first one runs dialog, the second sets values on default.
Tab Style.
Option Visible is used to show/hide plot point marking.
Option Draw every allows plotting every second, third and so on marking depending on
selected value.
Option Angle sets point marking rotation angle.
Option Clipped defines whether point marking is plotted or not if it is located beyond
graphic borders.
74
The squares at the top of the sheet are used for making seals, stamps or company logos.
Click the right mouse button on the square and in the appeared window select a raster image to
paste. Sizes of the square can be changed with the mouse.
At the bottom of the sheet editable table is located. To add text, click the right mouse
button in the table and in the created window select required text. It is also possible to save all
the comments in the table files with clicking on the button
clicking on the button
LONGITUDE
GRAV
479
-3.49.38.31598
-40.49.47.04385
156.645
977996.7625
478
-3.50.14.23385
-40.49.34.86847
176.749
977993.118
477
-3.50.41.81796
-40.49.12.96376
158.598
977996.16
476
-3.50.56.86075
-40.48.43.78614
153.596
977995.8746
475
-3.51.29.57703
-40.48.33.59533
163.494
977992.9616
474
-3.51.34.51079
-40.47.58.78146
150.784
977994.8251
76
473
-3.51.57.72425
-40.47.37.80829
173.923
977989.1426
472
-3.52.13.61573
-40.47.09.33508
151.016
977993.4726
471
-3.52.13.05078
-40.46.50.40375
164.658
977990.0839
470
-3.52.27.34988
-40.46.39.50728
170.948
977988.923
469
-3.52.44.12399
-40.46.27.26992
154.136
977992.3446
468
-3.52.54.85479
-40.46.09.14952
161.646
977990.9782
467
-3.53.10.34938
-40.45.59.98908
162.519
977990.9039
466
-3.53.27.27992
-40.46.01.57194
168.091
977990.0246
465
-3.53.45.00000
-40.46.02.08070
172.989
977989.2686
464
-3.53.50.90568
-40.45.46.82738
175.239
977988.88
118.3035394 0
126.9002384 0
1.5
123.4170888 0
116.1519574 0
2.5
117.240884
111.9424174 0
3.5
142.0405875 0
125.3686538 0
4.5
521.0730567 0
735.5232592 0
5.5
707.7315998 0
706.3561614 0
6.5
725.9945623 0
722.433627
77
7.5
717.0991126 0
716.9836552 0
8.5
725.5024012 0
722.3551713 0
9.5
731.5717173 0
10
723.5097884 0
10.5
726.8844987 0
11
725.962034
11.5
743.2485878 0
12
726.4061156 0
12.5
734.399887
13
727.9166309 0
13.5
116.1921851 0
14
517.9613065 0
14.5
125.3706264 0
15
111.2952478 0
15.5
131.911879
16
107.9217309 0
16.5
114.9327361 0
17
134.0939196 0
17.5
138.4457143 0
18
129.1165104 0
78
The first column contains measure point depth (from surface) of lithology horizon. The
second column should be completed by zeroes. The third column is the layer color on the
lithology scale. The forth column is texture type on the lithology scale.
Second type of files has *.crt extension; these are control files which specify type of data
and way of visualization. Structure of CRT file for lithology and logging data visualization for
any quantity of boreholes is described below.
2280.txt
brhl2280
18 2 2 1 0 1 0 0
3246.txt
brhl3246
102 2 2 1 0 1 0 0
80