Zond-IP Eng
Zond-IP Eng
ZONDIP1D
User manual
1
U r Rm J 0 mr dm .
r 0
K
к 100 % .
K 1 / 100
Value definitions
U
Ro_a – apparent resistivity, a G , where G is array coefficient.
I
U ВП
Eta_a – apparent polarizability, a 100 % , in percent.
U ПР
Separation: AB/2 for Schlumberger, Wenner and symmetrical array; AМ for two-electrode
array; AO for three-electrode array; O’O for dipole-dipole array.
Pseudodepth – approximate investigation depth connected with array coefficient (array
separation in this case).
To start up «ZondIP1D» it is necessary to create data file of certain format which contains
information on acquisition geometry and measured values of apparent resistivity and polarizability.
One profile data usually corresponds to one file.
In case of using multi-electrode arrays it must be remembered that sounding curve is
considered as data element in «ZondIP1D» program. That’s why data should be represented as a set
of sounding curves.
Text data files of «ZondIP1D» format have «*.zlf» and «*.zlp» extensions. See data file
format for details. The program also supports text files of IPI2WIN program.
For correct running of the program data file must not contain:
curves with less than 3 separations;
incorrect symbols of records separator (TAB and SPACE use only);
absurd data values.
Desirably, total record number should be no more than 5000 in one file.
Model of previous interpretation loads together with field data while opening «*.zlp» file.
Program presents universal data format which contains information on array geometry,
coordinates and relative elevations of sounding stations. «ZondIP1D» data files have *.zlf
extension.
First three lines contain information on measurement parameters which are the same for the
whole profile.
First line – must contain the following records separated by space or tab.
First record – is reserved by «ZondHED1d» program for frequency or seabed operations (8
values). Second record – (0-5) – type of array:“0” – Schlumberger, Wenner (length of potential
circuit must be specified for it), three-electrode array; “1” – dipole-axial array; “2” – two-electrode
array; “3” – three-electrode multi-purpose array; “4” – four-electrode multi-purpose array; “5” –
any multi-purpose array (pic. 2).
The dialog is meant for new data entry and can be found in the main menu of the program
File/VEZ Notepad. Picture 3 shows dialog window of VEZ Notepad option.
Copy
Paste
Start inversion.
Use XYZ field to specify coordinates of sounding station and its name (field «name1» in
pic. 3).
Table in the left part of the window contains the following graphs. First one or two columns
(depending on array type) are meant for array geometry specification. For Schlumberger array:
Use button to create new sounding site. When all parameters are set, press button
to start inversion.
Profile interpretation result is hold in ZondIP1D project files (extension *.zlp). This file
contains name of file with field data, fitted parameters and parameters variation limits for each
sounding. Use button in the tool bar or corresponding menu option to save interpretation
result. If autosave option is on interpretation result is saved automatically in selected periods of
time. Theoretical curves can be also saved in *.zlf format.
There are a few options to save data:
Project data file Save measured values and current subsurface model.
Calculated data file Save calculates values.
XYZ model file Save subsurface model as table file.
Section file Save current model in Section format. This file type can be used as
background.
Model with calculated Save current model and calculated values.
Model in columns Save current model as parameters table.
Data export
«ZONDIP1D» program allows exporting data to Microsoft Excel. This function can be
found in File/Export to tab. It is possible to export two types of data. Function Export to / Excel
report creates file which contains information on each sounding site, measured and calculated
values, its coordinates, models of parameters, etc. Function Export to / Excel map creates file
which contains names of stations, their coordinates, resistivity values (calculated), and positions of
When «*.zlf» file is created, load it using button or corresponding menu item. After
successful loading of the file functional buttons to work with data are activated in the main window
toolbar and short information on first sounding station appears in the right section of status panel.
While moving the cursor in created windows, coordinates which correspond to axes in every
window are displayed in the left section of status panel. Use the first system button of the
window, marked in plan as “objects setup”, to open context menu which allows running setup
dialog for each object in the window. Methods of editing are described in detail in the following
chapters (Model editor, Section, Parameters table, Pseudosection, Profile plot, Graph of
theoretical and observed curves).
Second and third buttons serve to proceed from one sounding station to another.
Style sheet
The toolbar is meant for quick run of the most frequently used functions. It contains the
following functional buttons (from left to right):
Add layer. User-selected layer is divided into two layers with equal
thicknesses (in logarithmic scale).
Delete selected layer.
Run inversion process for current station [left click] or for following
stations [right click] – profile inversion.
Copy model at current stake to buffer.
Paste model from buffer and use it at current station [left click] or at
all following stations [right click].
Cancel previous step.
The following table lists items found in the Menu with their corresponding functions:
“Hot” keys
The dialog serves to adjust autosaving options and specifics of starting model definition, set
default values and inversion parameters. Use button in the toolbar or corresponding menu
item (File/Program setup) to run it.
Tab Options
Tab Defaults
This tab serves to set different parameters used while working with model. Rows correspond
to types of layer parameters (properties and thickness), columns correspond to options (pic. 6).
First row Value – default values (used for starting model).
When starting model parameters are set use button to apply them to current project without
closing parameters setup window.
Pic. 10. Models created using different inversion algorithms (Standart – black, Focused – blue,
Smooth – red)
Option Threshold sets maximum contrast value of adjacent layers. When this value is
achieved parameters of adjacent layers are not averaged (thus boundary is considered to exist
between layers). It’s an empiric value (0.001-1). Low values of this parameter lead to algorithm
discrepancy (it’s recommended to increase it). High values result in smooth distribution.
Picture 11 shows models created using two different values of parameter Threshold (0.01 –
black, 0.1 – red).
Additional parameter - offset to prevent P-effect is assigned to each curve in the window.
This parameter is minimized for all curves during fitting and thus it significantly decreases P-effect
influence.
Current algorithm differs from standard inversion in additional parameters and construction
of smoothing operator. Parameters which model P-effect have less weight compared to the others.
Main characteristics of the algorithm are given below.
• Fitting is conducted simultaneously for all curves in the window and central point has the
greatest weight for discrepancy calculation.
• P-effect of each curve is selected during inversion process.
• Own model with common bottom part and variable top part corresponds to each curve.
Pic. 15. Results of 1.5D inversion algorithm (A) compared to results of standard 1D inversion (B)
«ZondIP1D» program allows solving 1D forward and inverse problems of VES. Program
automatically selects resistivity values and layer thicknesses. It is also possible to fix parameters
and specify their variation limits and define weights of certain measurements.
Linear filtering algorithm is used for forward modeling and Newton's method is used for
inverse modeling.
Interpretation mode is available when data file is read.
Data inversion
This option is used for quick data inversion of one or all profile stations. Newton’s method
is used for minimization of least-square deviation between calculated and observed curves. Left
click on button to invert current sounding station data, right click on it to invert the whole
profile. Information on current relative divergence between calculated and observed curves is
shown in the status bar. Inversion stops when current relative divergence between calculated and
observed curves or number of iterations is equal to specified beforehand.
Inversion gives non-unique solution since calculated curves can be identical for different
models. So it is necessary to take into account aprior information on geoelectrical section when
specifying starting model. If aprior information is reliable it is advisable to fix known parameters or
limit their change range in starting model and to add geological boundaries.
Program objects
It serves to display observed and calculated sounding curves and to visualize and edit
calculated model parameters curve.
Model editor window is shown in picture 16. Red and blue lines with filled circles are
experimental curves of apparent resistivity (left red axis) and polarizability (right blue axis) versus
separation in meters (bottom horizontal axis). Calculated model is displayed in red solid line.
model curve (cursor shape must change at that moment) and with left button pressed drag
selected part of curve. Green circle on selected curve indicates active layer.
Alteration of vertical curve pieces positions corresponds to alteration of model geometry (that
is thicknesses [right button] and depths of layers tops [left button]).
Alteration of horizontal curve pieces positions corresponds to alteration of model layer
parameters.
Double click in object area to run context menu which contains the following options:
It is displayed if Interpretation style is selected (tab Window). Window can contain one,
two or three similar blocks – graphs. On default left axis usually corresponds to red curves and
right one corresponds to blue ones. If only one measurement type is used in the program then right
axis is missing (pic.18).
The window consists of two parts. Editing curve is displayed in the left part. Separations
and values of editing parameter are shown in the right one. Right click to delete point on the curve.
Deleted point will be highlighted in grey in the table then. Left click to recover point. Use scroll
wheel with [ALT] button pressed to delete points. Change limits of deletion using scroll wheel.
Table editing window is used to change model parameters via keyboard. Table contains 3 or
4 columns. Each table row contains properties of one layer only.
First column contains resistivity values of layers, the second one (if present) – apparent
polarizability values, the third – thickness, and the last one – depth to layer top with sounding
station elevation taken into account. The cell is light grey on default (color can be specified) if
parameter variation limit is set. The cell is dark grey on default (color can be specified) if layer
parameter is fixed.
Right click on table cells to run context menu (pic. 21). If you right click on the first row
(table heading) then selected menu option will be applied to specified parameter of all layers (cell
is changed into col then). If you right click on the first column then selected menu option will be
applied to all parameters of specified layer.
Tab Limites
This tab serves to set parameter variation limits (pic. 22). Toolbar buttons are used to select
parameter type to specify limits to. [Auto] button sets limits to all parameters of selected type
automatically, according to their values and specified divergence percent. Model of selected
parameter (black), its lower (red) and upper (blue) variation limits are shown in the graph.
Use Equalence plot option to run equivalence model cloud plot window. It implements
quite resource-intensive algorithm of equivalence model access scan within specified error range
using curve-fitting method (pic. 29).
calculation is chosen then (Resistivity). After that access scan can be started . Algorithm result in
set of model curves. Use scrolling to proceed from one curve to another.
Table graphical parameters setup dialog (Options)
Option Lock sets fixed parameter cell color.
Option Range sets cell color for parameter with specified variation limits.
Option Free sets cell color for parameter without specified variation limits.
Option Min sets color of parameter minimum limit.
Option Max sets color of parameter maximum limit.
Option Active sets frame color of active cell.
Option Font sets cell font.
Option Cell height sets cell height.
Section
Current object visualizes geoelectrical section variations along profile. Graph is created in
profile coordinate versus depth axes. Color bar represents correlation between displayed parameter
value and color (pic. 30).
If mouse cursor is located within geological section then use slider to select and highlight
layer which is located below it. Cursor shape changes when it is reaching boundary of layers which
can be edited then. To do that drag selected boundary with left button pressed. If right button is
Double click in the object area runs context menu with the following options:
Bad data interpolation Interpolate models (with large fitting error) between two
selected sounding stations.
Pseudosection
Current object visualizes variations of observed values along profile as contour map
(pic. 36).
Double click in object area runs context menu with the following options:
Setup Run object parameters setup dialog.
Print preview Print preview dialog (see for details).
Select Options/Set lines in the main menu to set a number of profile lines. Line settings
window appears then (pic. 40). It displays areal positions of sounding stations.
A few profile lines can be specified in the program simultaneously. When profile and all
points around it are set, all included in profile points are displayed in blue. Left click to
include/exclude profile point. If profile line does not cross points then position of sounding station
projection is displayed in green.
Right click on necessary point to view and edit sounding station coordinates. Coordinates
which can be edited are displayed in opening window.
Press button to start data interpretation when profile is set. Main window of the
program appears then. Use window in the main window toolbar to switch
profiles.
Option Plane data (Options/Plane data) is used to create areal distribution maps of
selected parameter (resistivity, apparent resistivity, elevations, etc.) depending on depth or time.
Picture 41 shows example of apparent resistivity areal distribution map.
When plotting model parameters there are two windows in the toolbar which are used to
select one of parameters and specify depth to display it at. Depth counting method is specified in the
right window: Depth from topo – depth values are counted off the ground, Absolute depth –
absolute depth values are used, Layer index – contour map is generated for specified layer.
When plotting contours of measured parameters, layer number corresponds to initial data
time grid.
Field Axis scales is meant for axes scale definition. Scales can only be set with pressed
Aprior data (well data) availability helps to increase reliability of received geoelectrical
sections significantly. «ZONDIP1D» program has integral module which is used to visualize aprior
data graphically in sections.
Choose Options/Create/Edit borehole data option in the main menu to create stratigraphy
file. Dialog window of the module Add borehole data (pic. 45) appears.
Delete well.
Refresh window.
Advanced options.
The main toolbar also contains fields used to specify horizontal coordinate of well (offset
from profile beginning) - horizontal position and elevation - Position from surface. Vertical and
Press button in the toolbar to start creating lithology file. New table appears in Data
window. Use button to set number of layers, then edit table by specifying thickness and layer
bottom depth and select pattern type according to lithology. Double left click on C column to run
pattern setup dialog Pattern Color Editor (pic. 46). There is a wide range of lithology pattern in
the program. Use option Color to select pattern color.
When well data is entered, press button to display well graphically. Then specify
horizontal and vertical coordinate of well in kilometers in the toolbar and after that well is displayed
according to its coordinates. Active well is displayed in red in graphics window.
Palette can be saved to make work with large number of wells more convenient. In order to
create palette select pattern in pattern column of Data window, then right click on pattern domain
in the main toolbar. Set of patterns can be created in this way. To save it press button and
select Save default palette. Saved set of patterns can be used while creating new lithology and
Function Set percent, run via button , serves to change scale of well data graphical
representation.
Pic. 47. Display of lithology data in geoelectrical section (A) and in model editor (B)
Use Options menu to run boundaries setup dialog Set boundaries. It allows taking into
account aprior geological information. Menu with the following buttons appear when this tab is
selected:
Pic. 48. Example of geoelectrical section after inversion without (A) and with (B) regard to
geological boundaries
The dialog is meant for object palette setup and can be run using [Palette] button (pic. 49).
The dialog allows choosing one of default palettes (forward and inverse rainbow, grey scale, etc.)
or creating own one. Right click with [Ctrl] button pressed to add slider to the scale. Use [Delete]
button to remove slider. Use button to save scale or load already existent scale using
button.
Axes editor
Use File/Print preview function of the program main menu to run print preview dialog.
Double click on any program object also runs it. Two options are available when running via main
menu:
Option Station serves to print sounding curves and current station model (pic. 63). Model
parameters are displayed as table which contains number, resistivity value and depth of top layer
boundary.
- print setup button. Size and page orientation, print properties, number of
pages per sheet and other parameters can be set in opening window.
- press this button to start printing when necessary parameters are specified.
Lithologic columns are hold in certain file formats. First type of files has txt extension.
The following file structure is used to create lithology data file:
First column contains depth (from ground surface) of layer boundary. Second column
contains zeros. Third column defines layer color for visualization, forth – type of pattern
Lift of first 23 patterns which can be used for lithologic column creation is given below.
0 1 0 13 Top of layer 1
4 1 0 13 Bottom of layer 1
4 1 0 19 Top of layer 2
11 1 0 19 Bottom of layer 2
11 1 0 27 Top of layer 3
16 1 0 27 Bottom of layer 3
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 data visualization for any number of wells is
described below.
2280.txt First line – logging or lithology data file
скв2280 Second line – Well name (is displayed on well)
18 2 2 1 0 1 0 0 Third line contains control parameters.
Data record 18 – well coordinate on profile.
2 – image width (in percents to profile length, usually 1 - 20).
2 – type of data visualization 0 - 3.
0 - logging data (as graph); carot1.crt
1 - logging data (interpolated colour column), section colour scale is used for visualization;
carot2.crt
2 - lithologic column; strati.crt
3 - logging data (colour column), colours for data visualization correspond to model colour scale,
column colours are selected in compliance with model colour scale;;