admin管理员组文章数量:1530845
2024年7月16日发(作者：)
MAGPICK-磁异常图及剖面数据处理
使 用 手 册
译者：badboy
2007-1
Mikhai lMikhail Tchernychev. 16/02/98
Last revised 5/13/2004
目 录
1MagPick的发展及特点 ................................................................................................................................................. 1
1.11.0版1996-1998 .................................................................................................................................................. 1
1.22.0版1998-1999 .................................................................................................................................................. 1
1.3Version 2.x - 2000 - 2004 ...................................................................................................................................... 1
2安装 ............................................................................................................................................................................... 2
2.1系统要求 ............................................................................................................................................................. 2
2.2编辑 ..................................................................................................................................................................... 3
2.3 MS Windows环境下安装 ................................................................................................................................... 4
3MagPick基本操作 ......................................................................................................................................................... 4
3.1 Magpick图形(网格)阅读器 ................................................................................................................................ 6
3.1.1Transforming equalized color palette into gradient. ................................................................................. 14
3.2Drawing and clipping on top of the map ............................................................................................................ 14
3.2.1Magpick vector formats: lines, polygons, points, clip area. ..................................................................... 14
3.2.2Bringing an AutoCAD DXF drawing into Magpick. ................................................................................ 15
3.2.3Using ArcInfo(TM) shape files in magpick .............................................................................................. 15
3.3Simple picking of magnetic anomalies ............................................................................................................... 16
3.3.1Simple pick. .............................................................................................................................................. 16
3.3.2 Pick export ............................................................................................................................................... 17
3.4Automatic pick .................................................................................................................................................... 18
3.5Profile information loading and viewing ............................................................................................................ 19
3.5.1Profile loading. .......................................................................................................................................... 19
3.5.2Easy profile load. ...................................................................................................................................... 21
3.5.3Operation with profiles on the map. ......................................................................................................... 23
3.5.4Parameters of the profile view .................................................................................................................. 23
3.6Printing ................................................................................................................................................................ 25
4Grid operations. ............................................................................................................................................................ 25
4.1Making grid out of profile data. .......................................................................................................................... 26
4.1.1Triangulation with linear interpolation ..................................................................................................... 26
4.1.2Gridding with splines. ............................................................................................................................... 26
4.2Save grid. ............................................................................................................................................................ 26
4.3Create clip path. .................................................................................................................................................. 26
4.4Clip the grid. ....................................................................................................................................................... 27
5Profile data transformations. ......................................................................................................................................... 28
5.1Saving and re-loading profile data. ..................................................................................................................... 28
5.2Profile data operations ........................................................................................................................................ 28
5.3Profile position operations. ................................................................................................................................. 28
5.3.1 Simple shift .............................................................................................................................................. 29
5.3.2 Spline smooth. ......................................................................................................................................... 29
5.3.3 Shift / dragging along profile. .................................................................................................................. 29
5.3.4 Dragging method details. ......................................................................................................................... 29
6UTM / Gauss-Kruger transformation in magpick ......................................................................................................... 29
7地球磁场模型(IGRF). ................................................................................................................................................. 30
8位场转换 ..................................................................................................................................................................... 30
8.1上延(Upward continuation) ............................................................................................................................... 31
9Estimation of magnetic sources. ................................................................................................................................... 32
9.1Magnetic modelling based on grid data. ............................................................................................................. 32
9. 2Magnetic modelling based on profile data. ........................................................................................................ 35
9.3Magpick worksheet - results of inversions. ........................................................................................................ 35
9.4Option of mass estimation .................................................................................................................................. 35
9.5Automatic estimation of magnetic sources. ........................................................................................................ 35
10Implementation of the inversion in magpick .............................................................................................................. 35
11Export to the CAD and other graphic programs ......................................................................................................... 35
12Tutorial指南 ............................................................................................................................................................... 35
12.1View grid files with MagPick. .......................................................................................................................... 35
12.2Load profile information into MagPick. ........................................................................................................... 39
12.3Putting things together: map and profile data. .................................................................................................. 44
12.4Detecting of 51
12.5Finding magnetic dipoles .................................................................................................................................. 53
12.6Finding magnetic pipes. .................................................................................................................................... 58
1MagPick的发展及特点
1.11.0版1996-1998
提供基本的异常解释操作，例如用网格数据最大和最小值调节异常图颜色；移动窗体显
示部分图形，自动改变颜色；在一般情况下，磁性体在异常的最大值和最小值之间，程序能
用鼠标取两点并将数据保存在文件中；选择最大和最小的区，能自动找到其精确位置。
1.22.0版1998-1999
没有被广泛的使用。
1.3Version 2.x - 2000 - 2004
-1-
根据用户需要，增加了以下功能：导出DXF文件功能；剖面数据编辑功能，样条曲线滤
波；网格数据编辑；网格数据转换，例如向上、下延续，梯度(倾斜度)估计；网格数据滤波，
样条曲线插值等；调色板导出GMT和SURFER格式；
2安装
2.1系统要求
MagPick是C++ 程序，运行在U?nix操作系统下。现行版本的MagPick在Linux (U?nix
clone for PC’s)发展和使用。由于MagPick是一个完全的视窗客户机程序，因此必须有X11
运行。它需要连接下列的X11公共库：
-2-
选择系统库与另外的操作系统不同，例如Solaris2.x。
MagPick使用netcdf库阅读netCDF文件，最近发表的CDF文件能从下面的网址获得：
ftp:///pub/netcdf/.Z
下面说明在微机上安装netCDF库。
为了执行制图用户界面(GUI)，magpick是用wxWin C++库(1.6x版),通过以下网址找
到：
//wxwin/
有不同运行wxWin方法，库建立在标准的工具包：xview (OpenLoolk) or Athena
widget set (Xt)。一般的MagPick能用这些平台进行编译。而正在工作的版本能用wx-Xt
库编译(版本166d)。下面说明安装wxWin库，在设置期间要确信Xlib被激活，你需要的wxgrid
库是wxWin的部件。从：
ftp:///pub/packages/wxwin/contrib/wxgridd/wxgrid
获得工作表。有必要使用test/demo程序检查设置。当设置完这些库和gcc C++编辑
器，你可以编辑MagPick。
2.2编辑
建立子目录并用gzip和tar程序下载文件。编辑文件要与
-3-
你的设置相匹配。需要改变下面的文件：
WXDIR指向安装的wxWin安装的根目录，
INC指向包含的文件夹，假如netcdf库没有在系统的位置，也要设置包含的路径。
THISDIR当前放置原始资料的位置。
LDFLAGS连接标记，那里所有库能被找到，应是wxWin和netcdf库的路径。
LDLIBS选择库，假如在你的系统wxWin中需要另外的库，应改变列表，在改变makefile
后，你能键入make-f 文件，使magpick xt成为可执行版本的程序。
2.3 MS Windows环境下安装
假如有文件，可能有的问题是库，需要正确的平台。下一步
是复制这个文件成为Windows系统文件，可能的目录是在NT下C：WINNTSYSTEM32，
对于95/98系统在C：WINDOWSSYSTEM。注意：该文件不能存放在与同
一位置。后来版本的magpick有一个安装程序。
3MagPick基本操作
MagPick用网格文件数据和剖面文件。网格文件：SURFER ASCII，NetCDF格式，剖
面文件，用空格/TAB分隔的ASCII文件。网格文件是必需的，另外一些数据如路线图/线路
数据放在网格的头部。
-4-
装载网格文件不是合成的，必须是能MagPick被阅读的文件。
剖面文件，必须注意剖面类型，MagPick能准备这种文件而很少能控制。不必装载图形
文件而得到剖面数据。
在剖面数据文件被装载后，可用内插法建立一个能被MagPick识别的文件。
在磁异常图上，能够设置标志，可以描绘你感兴趣的地方，MagPick能把它作为一个文
件保存并可再次调用。MagPick能导出其它格式SURFER和AutoCAD文件。
能估算磁性源的位置。能辨别两种磁性源：磁偶极(
Magnetic dipole
)，它是一种有效的
本地磁信号源；另一种是线性磁性源(
magnetic line
)，例如管道。在一些情况下，能辨别磁
性体的几何位置和强度。估算结果可作为工作表(worksheets)保存为文件并导出。在U?nix
启动MagPick下，用命令：magpick xt；在MS windows下：，MagPick
首先读文件(从MSWindows目录下- windows 95/98；用者剖面目录-
Windows NT)。这些文件控制程序设置，通过对话框能够转换，在MagPick退出时，这些
文件能保存。
如果文件没有找到，使用缺省的设置。如启动成功，在U?nix下将见到小窗体菜单(X不
支持MDI界面)；在window下出现空的多行文本文件。这是MagPick主窗体，将保持不
变(U?nix)。此外，主窗体有另外的子窗体，现在版本的MagPick有三种类型的子窗体。
(1)图形窗体：显示2-D彩色或黑白网格图形，图形与其它窗体没有关系；
(2)剖面浏览窗体：一组剖面文件装入程序，可浏览每一个2-D线的图形。窗体可以和图
-5-
形窗体连接；
(3)剖面倒置观察：可看到穿插母图的所有或部分剖面。仅有一个剖面可观察图形。在剖
面未装载前，剖面不能观察。当选择时，将产生一个新的空窗体。假如文件被装载，
然后选择时，能见到想观察的图形。
3.1 Magpick图形(网格)阅读器
启动Magpick选择建立新的窗体，单击选择netCDF或SURFER (ASCII
或binary)网格文件装载进入窗体。彩色图形将出现在窗体上。改变大小选择
/，可选择网格单元的水平大小(像素)，垂直和放大因数。注意大尺寸增加微
机内存。可用
Autoscale on zoom
，注意与最大最小值一致，可观察图形细节。
如显示的图形不满意，应调节颜色，网格线，轮廓线。另外可通过/
出现“common parameters”对话框，设置参数(参见图1)。其中：
1 Field direction表示区域方向；
2Data Min Data Max最小数据-最大数据：当“autoscale”未被选择时，可定义网格
的标图界限，调色板调节大数据使用大号的颜色。但选择“
Wrap color scale
”，最小数据-
最大数据有不同的意义；
3自动比例(Autoscale)：按最大最小数据自动调节；
4约束颜色比例：使用“color scale wrapping”，可见图形细节。如数据范围从0到200，
颜色比例在此范围；
-6-
图1：通用参数对话框。
5Color Map:彩图.调节彩图和消失颜色。这个控制被标记为(*)，单击检查另外的设置(对
话框)，在对话框中可选择彩图的几种模式，包括：
No color map
：没有颜色；
Simple color map
：数值用颜色表示，
Shaded relief
：设置方位，高度和角度(
azimuth
and
elevation
, degrees)，另外要设置
数据比例(
Data scal
)，一般为1—10，应使用相同的调色板(
Uniform palette
)设置。
Color illuminated map
：一个简单的彩图与阴影图结合效果更佳，尤其是地形图。
-7-
6改变调色板()；
7Show crosses；
8 Show profiles:显示剖面;
9等值线contours；
10Snap to grid:网格对奇;
:网格颜色；
12Show X Y grid：显示XY网格；
13Mesh steps.网格间距设置。
Contours被选择，弹出等值线设置见图2。
1Fixed values/Automatic：固定值和自动形成等值线；
2Total：当使用Automatic时，
3Min., Max：等值线最大、最小值；
4Interval：等值线间隔；
-8-
5Calculate：输入等值线最大、最小值、等值线间隔后，需按下“Calculate”按钮后有
效，等值线数值在列表框显示；
6Isolines：等值线列表，可以计算或从文件加载；
7 Read ：等值线数值的ASCII文件，结果在等值线列表框显示。
：定义等值线颜色；
9 Fast drawing：
10 Label isolines：等值线注记；
11 Decimals.：注记的小数位数；
12 Distance between labels, map units：控制绘制等值线，程序计算等值线距离，
当超出距离，则应在XY单位指定距离；
13 Label length, map units：注记文本长度。
-9-
图2：等值线参数对话框
输出剖面参数：
1Show profile direction：显示剖面方向；
2 Interval between signs.：建议100-200像素；
3 Plot stack profile：如要观看一组剖面，此项应被选择；
4 Datasets to draw：每个剖面文件有几种类型数据，
5 Data scale:
-10-
6 Preferred direction:
7 Type:
8 Fixed direction azimuth:
9 Fill type:
10 Fill color:
11 Clip stack plot
12 Min:, Max:
13 Middle line approximation
14 Stack labels
-11-
图3：剖面参数对话框
-12-
图4：建立剖面的不同方法
图5：剖面标注实例
图6：剖面标注对话框
虚拟网格文件格式：
-13-
DUMMY 文件头
100 100 列数和行数
62740 63200 X坐标极值
33130 33820 Y坐标极值
-9366.12 17963.5 虚拟场值的最小和最大值
3.1.1Transforming equalized color palette into gradient.
3.2Drawing and clipping on top of the map
3.2.1Magpick vector formats: lines, polygons, points, clip area.
4种类型信息能被输出：
1LINE：准备下列格式的文件：
> UTM GREY 0
0 0
10 10
-14-
> UTM RED 0.1
10 10
20 0
每一线起始用符号：>，便于与GMT系统兼容。UTM指示为米制；GEO指示是经纬度。
2Point：准备下列格式的文件：
POINTS 12 点 字体大小
UTM 5 5 10 RED Point num 1 点1 X Y 注记大小
UTM 10 10 10 BLACK Point num 2 点2
3Polygons：多边形
4Clip.
3.2.2Bringing an AutoCAD DXF drawing into Magpick.
导入DXF文件。
3.2.3Using ArcInfo(TM) shape files in magpick
-15-
3.3Simple picking of magnetic anomalies
3.3.1Simple pick.
有两种拾取的方式：
Simple pick Go /Simple pick
鼠标变为“+”，移动鼠标到第二个点，单击左键。如果Snap to grid选择为开则拾取的
结束点总在网格中心，
Square picking
矩形拾取，鼠标改变形状，拾取。
利用/ /保存文件(ASCII file。每一拾取表示
为一行，字川格式：
1 X Y 拾取部位最大值(3个数据)
2 X Y 拾取部位最小值(3个数据)
3 从最大值(第一个点)到最小值(第二个点)的方位+场方向值，
4 拾取点水平距离
-16-
5 拾取中心点坐标
6 总是0.000(默认)
当拾取产生错误时，用/Delete pick删除。
3.3.2 Pick export
在
/Export pick..
菜单导出*.DAT、BNA and DXF文件。有两种方式：
沿X Y轴线性导出，包括缩放比例和偏移；
UTM方式的转化，参考Universal Transverse Mercator / Gauss Kruger in MagPick。
应用于经纬度的情况。
-17-
图10为导出文件设置。
3.4Automatic pick
步骤：
1找到本地最大值、最小值的所有点；
2建立一对最大值和最小值，这些对值应十分大，最大值与最小值匹配，减少数据对数量
(见图11)；
-18-
3用户定义两点距离，进行处理；
4用户定义两点的场值差异；
5方位，例如北半球异常从最大值应该由最小定位到北方，因此方位范围应符合这一情况。
6中间点场值在端点场值之间。
在Edit / Auto pick，选择角度、距离和场值范围，
3.5Profile information loading and viewing
3.5.1Profile loading.
-19-
为任何ASCII文件，文件至少有三列：X Y 场值。为表征剖面特征，应建立用文件名、
列数或
Easy profile load
一个特殊文件，后者访问剖面简单但有限制，例如每一个文件仅有一
种数据设置，
可以用
Easy profile load
装载文件或作为文件列表保存，文件列表可作为模板。数据文件
每一行包含了文件装载信息，场值用空格分割，行格式为：
1“SPLIT=NN”，SPLIT指示数据设置必须被分离为行，NN可以是整数或一个带.的实数，
整数时，NN表示文件行数；MAGPICK假设当字串改变时为新行开始；实数时，被作为点间
的距离处理，分离数据成行，假如距离>现在点间距离，则为新行开始。
2 XFALSE=NN1和YFALSE=NN2：NN1和NN2为假设的北和东，文件中被加到XY坐标。
3剖面文件名，如不在当前目录则带路径。
4磁场图顶部不活动的图形颜色名。
5活动的剖面颜色名。
6当前文件数据量设置。一个文件能包含多个剖面。
7程序仅能读部分数据(如每次读10个数据)，为增加读数和后期处理速度，内部执行。
8XYZT数据顺序，全部数据为
the number of data sets
×
4
.，一般情况下Z轴向下，文
件中无Z，则指定一个常数H=Value代表Z列。
-20-
9下一个信息时可选的，反映了磁力仪排列特征，可以计算磁力仪探头高程，
文件实例：
当从一个数据文件装载几个剖面数据时，第一列是#number1_number2，number1剖
面数据起始字串号；number2剖面数据结束字串号，例如：
3.5.2Easy profile load.
(1)选择Go to
File / Profiles / Simple load
.
(2)弹出图12对话框：
-21-
在拾取数据前输入假设北、东。
(3)单击Add出现下部的对话框，设置分离线的方法，假设文件包含一行则不分离。
(4)键入
Decimation
，可装载部分数据，假设
decimation
=2，则仅每个第二个点装载。
(5)假如有场值显示的文件头，显示包含XY场值的列，否则显示列号。
(6)活动剖面和非活动剖面颜色，按OK确认。
-22-
(7)新的字串出现在列表中，通知MAGPICK如何处理文件。用Add按钮增加许多文件。
(8)拾取全部文件后，选择
Save as profile list
，保存文件。
3.5.3Operation with profiles on the map.
加载后，文件列表出现在程序中。每一个文件分为活动和非活动两种状态，活动文件高
亮显示。当选择时，出现对话框，观察选择分为三种：
网格观察；
剖面观察；
所有剖面观察。
文件装载后，Edit…菜单有效，可进行一下操作：
方向：建立一个坐标轴，单击左键显示坐标轴起点，再次单击显示方向，
拾取剖面：
设置标记：
3.5.4Parameters of the profile view
选择/，显示对话框(图14)：
-23-
剖面显示有两种模式：
1收缩模式：剖面收缩与最小、最大数据相一致状态；
2固定比例：
剖面参数对话框内容：
1场值限制：
2水平轴限制：
3数据限制：
4X轴类型：取决于建立剖面观察设置，有2—4种选择：
-24-
(1)无另外坐标轴的剖面倒置观察；
(2)沿剖面XY无另外坐标轴的剖面倒置观察；
(3)沿剖面和指示器；
(4)沿指示器
5收缩；
6网格；
7尺寸：选择5为适合指定尺寸。
8比例。
下一个重要菜单是，
3.6Printing
御览，
打印设置对话框Print setup，
4Grid operations.
-25-
4.1Making grid out of profile data.
在
Files / Interpolate grid
菜单和图形观察
Profiles / Interpolate grid
。基本步骤：
1加载文件数据，
2如主要目的是制图，在
Profiles / Interpolate grid
，拾取添加方法；
3添写参数添加；
4添加完成后，MAGPICK保存为一个文件，自动装载新的计算图进入当前图形观察点，
5调节观察参数；
6观察使用
File / New
和
Grid view
，
4.1.1Triangulation with linear interpolation
用线性插值分成三角形，用于离散点插值，形成规则网格。图15。
4.1.2Gridding with splines.
4.2Save grid.
4.3Create clip path.
-26-
在
Edit / Polygon
，鼠标变为笔，绘一个多边形，左键增加节点，右键删除最新的节点，
单击第一个点关闭多边形。
在
Operations / Save clip
，
文件保存后，
4.4Clip the grid.
-27-
5Profile data transformations.
5.1Saving and re-loading profile data.
5.2Profile data operations
5.3Profile position operations.
-28-
5.3.1 Simple shift
5.3.2 Spline smooth.
5.3.3 Shift / dragging along profile.
5.3.4 Dragging method details.
6UTM / Gauss-Kruger transformation in magpick
-29-
7地球磁场模型(IGRF).
地磁场计算，利用文件计算。
8位场转换
包括上延、下延，减少磁极、磁场分解，转换使用快速富里叶转换(FFT)。转换应使用剩
-30-
余磁场数据，在转换前，观测数据应减去地磁场数据，即IGRF数据。
8.1上延(Upward continuation)
F
(
u
,
v
) =
s
(
u
,
v
)×
f
(
u
,
v
) =)
F
(
u
,
v
) =
s
(
u
,
v
)1+
s
(
u
,
v
)2 ×
f
(
u
,
v
)
式中：
f
(
u
,
v
)为被转换的磁场频谱；
F
(
u
,
v
)为上延或下延频谱；
s
(
u
,
v
)为转换频谱；
为调整参数。
延拓计算在菜单“
Operations / Upward
”。
8.2水平垂直梯度合成(Synthetic horizontal and vertical gradients).
当上延计算时，垂直梯度能被计算出。
8.3 Reduction to the magnetic pole and pseudo gravity rans ormation
-31-
如果地磁场和磁化方向已知，则全场异常能被简化到地磁极异常。参数设置：
9Estimation of magnetic sources.
9.1Magnetic modelling based on grid data.
当磁场图装入后，进行处理包括以下步骤：
1决定磁极(偶极)的数目；
2借助缩放工具部分的选择子窗体，得到一个或更多的磁异常分区，
-32-
3单击磁极的大致位置一个或多的简单拾取，进行解释。然后调用
/Run..
，设置
参数运行。
4分析结果出现在工作表中，
5分析结果，决定是否重新输入参数
6重复上述步骤
7
/
.可移去或放弃结果。
为发现磁极，必须提供本地坐标位置和地磁场参数，这些参数必须由
Inversionparameters”加载，参数：
I—磁倾角；D—磁偏角；A—X轴与磁北夹角。单位体积的磁化系数(cgs-厘米克秒制)；
-33-
单位体积的剩余磁场，gauss单位。高程。
页边距：
Margin space.
反演次数；
反演误差；
起始深度；
密度；
J数值，
如推测物体仅有感应磁场，
仅是线性答案；
网格类型；
位置文件：推测结果除了在工作表外，另外保存位置文件，格式与拾取文件一致。
计算出的网格；
误差；
-34-
9. 2Magnetic modelling based on profile data.
9.3Magpick worksheet - results of inversions.
9.4Option of mass estimation
9.5Automatic estimation of magnetic sources.
10Implementation of the inversion in magpick
11Export to the CAD and other graphic programs
12Tutorial指南
12.1View grid files with MagPick.
以文件为例：
1Check initial state检查初始状态：选择菜单
Options/Size
，弹出对话框，确定
X-size,
cell
和
Y-size, cell
为1。选择
Options/Settings
菜单，弹出普通参数对话框，注意除了
Autoscale
和
olored map
外其它不选择，单击“
Change palette
”，选择
Uniform
，OK确定，
弹出
color palette
对话框，调色板号为34，选择“
Equalize colors
”。操作完成后，准备加载
数据。
2Pick and load the file. Inspect result拾取和装载文件
-35-
在
File/Open
菜单，打开文件，显示图形。用菜单
Options/Info
能观察网格
的几何参数。
3Change geometrical parameters of the view
启动普通参数对话框，设置观察几何参数。主窗体状态栏显示First和当前鼠标位置的XYZ
值。
4Zoom a picture.
观察图的细节使用：用菜单
Edit/Zoom
鼠标改变形态，鼠标左键和拖拽。注意，如果老
的窗体使活动的，假如
autoscale on zoom
被选择，则颜色比例自动调节。
5Change data color mapping.
启动普通参数对话框，如选择“
Autoscale
”，则Data Min和Data Max不可用。不选择
“
Autoscale
”，Data Min=1500，Data Max=2000，<1500为蓝色，>2500为红色。
6Color palette.
假如
Equalize colors
未被选，则颜色是线性的；否则，为非线性，每一种颜色在屏幕占
据相同面积，因此色彩丰富。
7Advanced mapping techniques.
高级设置。包括阴影和亮度。在普通参数设置选择“
Change palette
”，确定设置为
-36-
“
Uniform
”后在普通参数对话框单击“
Colored map
”选择框，出现“color map draw
mode”对话框，选择“
Shaded relief
”，选择光线角度315，高度25，数据比例5，将看到
黑白阴影图。
8Make your map nice - add contours (isolines).
图形美观，增加等值线。在普通参数设置单击“
contours
”复选框，弹出“
contours
”
对话框：
-37-
选择
Fast drawing
，输入Min和Max及间隔，单击“计算”则在“Contours”列表框显
示等值线值；不选择“
ast drawin
”，则“Label contours”可用(标注等值线)，“Decimals”
小数点位数，“Distance detween labels”沿等值线按设置的长度标注，“Label length”标
注长度。
9Printing your map.
打印图形。MAGPICK打印容量有限，打印图形总是与打印纸相匹配，因此需合理的选择
打印比例。打印时选择“
Options/Title
”菜单，弹出对话框：
-38-
输入题目，字体大小。关闭对话框，在“file”菜单选择“
Preview
”预览：
12.2Load profile information into MagPick.
剖面文件导入MAGPICK。需要的文件：和，通常由Paul Kot提
供。文件为Magmap Geometrics程序导出的计算结果，描述G-858磁力仪磁场数据。
文件为本地设置，用于GPS接收器，California地区。
1Loading one profile data from one file.
从一个文件加载一个剖面文件。用“
File/Profiles/Simple Load
.”菜单：
-39-
单击“Add…”按钮弹出普通对话框，选择文件“”文件，弹出对话框：
“”文件有一个各列名称的头字串：
X Y TOP_RDG BOTTOM_RDG VRT_GRAD TIME DATE
LINE MARK
200.00 99.58 49975.415 50032.320 113.810 12:02:09.60
05/26/99 50 0
-40-
……………………………….
如没有头字串则仅知道列数。应知道哪列表示XY和场值。数据分为9列，TOP_RDG为顶
部传感器磁场；BOTTOM_RDG为底部传感器磁场；VRT_GRAD为磁场插值。
在列表框中选择数据类型，设置完成后，单击OK对话框关闭，回到“Load XYZ Prifiles”
对话框，显示新的字串，单击OK，提示文件加载。
2Loading several profiles (lines) from one data file.
单击“Add…”按钮弹出普通对话框，选择文件“”文件，弹出对话框，在
“Split”列表框中，选择“distance”分离线。
如文件有一特殊的列表示线的改变，则可用该列分离线。此时在“Split”选择“Line”
则线分离并全部加载，可在file菜单的“Browse”浏览各线的曲线：
3Viewing of profile data.
观察剖面数据。在“
File/New
”，弹出的下面对话框选择“Profiles view”，单击OK，
-41-
在主窗体显示空的银屏和一垂直红线(右边)。缺省显示T(x)数据，但“”沿
着Y方向的线足够长，为改变观察到T(y)，可在
Parameters/Settings
选择类型，选择
Settings
菜单后弹出对话框：
在左下角列表框选择Along Y，单击OK，可见到改正后的曲线图。在状态栏可见到鼠标
所在位置的T(y)值。
如设置为“
Along profile
”，水平坐标为沿着线的距离(线可以不是直线)，距离总是正的
且从0开始，
如选择“
Shrink to canvas
”复选框，为缩放显示。
4Browsing profile data.
-42-
在file菜单的“Browse”浏览各线的曲线。如何缩放，选择菜单“
Parameters/Settings
”，
不选择“
shrink to canvas
”，
T scale=
1，
X scale=0.
1，则曲线比例变化。
5Add another file to the profile list
增加另一文件到文件列表。
用菜单选择“
File/Profiles/Simple Load
”，加载“”文件，数据格式：
X Y TOP_RDG BOTTOM_RDG VRT_GRAD TIME DATE
LINE MARK
6156055.370 2085444.39 49979.053 50016.776 75.446 12:02:09.30
05/26/99 50 1
……………………………………..
选择“Line”分割剖面，数据选择为底部传感器数据，OK确定，回到起始对话框，OK
确定。弹出：
选择“Yes”则剖面总数变为102个，创建新的视点观察。
-43-
6Saving file list for future use.
当用菜单选择“
File/Profiles/Simple Load
”时，对话框有一“
Save as profile list
”复
选框，选择则提示是否保存剖面分离的文件列表，如保存为列表文件，则再次启动MAGPICK
用“
File/Profiles/Load list
”可打开所有剖面。
12.3Putting things together: map and profile data.
1Load profile data.
启动MAGPICK，用“
File/Profiles/Simple Load
”加载文件，选择Line分
离线，“
Save as profile list.
”复选框为选择，
False Easting
和
False Northing
设为0，OK确
定，文件对话框弹出，选择“”文件，确定。
2Loaded grid data.
选“
File/New
”菜单，在弹出的对话框选择“
Profile view
”， 然后用菜单“
File/Interpolate
grid/Splines
”(即网格插值/样条曲线方式)，弹出“
Interpolation parameters
”(插值参数)
对话框：
-44-
“
X parameters
”选“Data limits”，表示X范围与X最小、最大坐标相匹配；
X间隔选择为2.5；
“Y
parameters
”选“Data limits”，表示Y范围与Y最小、最大坐标相匹配；
Y间隔选择为2.5；
“
Low constrain
”最小约束数据值；“
High constrain
”最大约束数据值，
“
Tension
”设为0. 25；“Format”设为“GS ASCII”(Golden Software ASCII)。
“
Iterations
”设为250(反复数)，大多数情况足够用。
“
Convergence limit
”(收敛极限限制)设为0；
-45-
单击“
”，输入文件名，保存网格文件。
3View interpolated map
观察图形。选“
File/New
”菜单，在弹出的对话框选择“
Grid view
”，出现空的窗体，
用“
Options/Size
”菜单，X Size和Y Size设为1，用“
Options/Settings
”菜单，进行以下
设置：
选择“
Autoscale.
”；
选择“
Show x grid
and
Show y grid
”，XY设为10；
不选“
Snap to grid
and
Contours
”；
选择“
Show profiles
”，在结构对话框“
Profile plot parameters
”不选“
Show profile
direction
and
Plot stack profile
”；
单击OK；
选择“
Color map
”，在结构对话框“
Color map draw mode
”选择为“
Drawing mod
”，
单击OK；
单击“
Change palette..
”，选择“Uniform”，OK后在弹出的对话框选择调色板33，选
择“
Equalize colors
”。
-46-
回到菜单
File/Open
，打开文件，在屏幕上有两个窗体，一个是剖面数据，
另一个是用垂直线表示剖面位置的磁场图。
4Simple profile - map interaction
现在设置剖面视点和水平轴设置。在剖面窗体上，在
Parameters/Settings
菜单，Type
设为Along Y，注意，当前的剖面线为高亮的红色。其它剖面用菜单“
File/Browse
”观看。
回到主图上，用菜单“
Edit/Pickprofile
”编辑，光标改变形状，在一些线上用鼠标左键，
点中的线变为红色，可直接观察剖面而不需要“
Profile browser
”对话框。
下面的操作时如何绑定点，用设标记的方法，即向上或向下的箭头。第一次是向上的箭
头，第二次则为向下箭头，第三次做时低一个箭头消失，新的向上箭头出现在新的地方…。
在剖面观察窗口，用菜单“
Edit/Set marker
”，光标变为+字形，单击曲线出现箭头；
在主图形窗体上，再次用菜单“
Edit/Set marker
”，光标变为⊙形，其圆心放在红线的
地方，单击出现向下箭头，在剖面窗体箭头随之变化出现在同样地方。
重复上述操作，注意在剖面端点可能看不到一些类型的标记，扩大图形，标记能被放置。
注意，只能在活动的剖面设置标记。
5Another approach to grid profile data.
样条曲线内插法可能形成较宽的异常，与实际不符。可以改变“
Tension
”参数，
Tension
变为“a minimum curvature”(最小曲率)，合理值为0.25-0.8。
-47-
有另外一种插值方法，在剖面窗体用菜单“
File/Interpolate grid/Triangulate
.”(三角形
插值)，弹出对话框：
在
X parameters
选择为
Data limits
，设置间隔0.25；
在
Y parameters
选择为
Data limits
，设置间隔0.25；
选择“
Filter data (recommended-
推荐
)
.”；
设置
Method
为
Average value
，
Format
为
GS ASCII
。
最后变文件名为，t表示由三角形而形成的。
现在，再次用“
File/New
”和“Grid view”等菜单观察图形。
6Draw stack profile.
-48-
绘一组剖面。上述观看的剖面是在分离的窗体上。可以在一张图上同时看到多条剖面。
用“
Options/Settings
”弹出普通对话框后，双击“
Show profiles
”，出现“
Profile plot
parameters
”对话框：
对话框中：
选择“
Show profile direction
”—每200个像素显示剖面方向；
选择“
Plot stack profiles.
”，
设置
Datasets to draw
为1，一般情况下为1；
设置“
Data scale: 100
.”；
选择“
Type
”，为“
Fixed direction
.”，即固定方向，意味着正异常绘到指定方向，在“
Fixed
-49-
direction azimuth
.”输入90，选择“
End points
”类型，意味着剖面端点作为计算方位的点；
选择“
Normal to path
”，意味着一组剖面的每一个点垂直于路线。
“
Preferred direction
”，正异常部分移动的方位；
选择“
Fill type:
”为
positive
，填充颜色为绿色；
现在不要选择“
Clip stack plot
”；
OK关闭对话框，关闭普通对话框，图形分布一组剖面。在“
Options/Settings
”不选彩
色图，可清楚地观看剖面。
7Advanced feature - how to add something on top of magnetic map
磁场图观察几何或地形信息。
另外的信息可来自USGS DLG、AutoCAD DXF、DLG, ArcInfo files等等。
用记事本新建一个数据文件，；
键入如下数据：
> UTM GRAY 5
10 10
-50-
10 70
70 70
70 10
10 10
然后保存文件。
在图形窗体下，用菜单“
Options/Drawings
”，弹出对话框，选择“Add”；弹出对话框，
选择“Lines”，OK确定。选择重新绘图，在图形左下角出现一个大的灰色的矩形。如果矩形
大，可改变文件> UTM GRAY 5为> UTM GREEN 0.，重新绘图。
UTM表示UTM坐标格式，GEO表示经纬度格式。
12.4Detecting of magnetic bodies.
探测磁性体。举例：226.5Kg炸弹，用装有上下两个传感器的g858磁力仪探测，网格为
17X20m，磁化系数推测为0.87cgs，剩余磁化0.24高斯。文件放在bomb子目录。
1Simple pick of magnetic object.
启动MAGPICK，加载文件，为35×41网格。用“size”菜单将XY设为10，网
格间隔设为2×2，进行不选择均等色彩设置，显示图形。下面进行定位：
-51-
选择菜单“
Edit/Simple pick
.”，光标变为+字，移到认为最大值的点击左键，出现一虚
线，然后移到认为最小值的地方再击左键，一个实心的黑色线段出现，线段的中心则为炸弹
的位置。
在普通对话框，如选择“
Snap to grid
”，则起点和终点被捕捉到网格单元中心，假如不
选择“
Snap to grid
”，则为停在所设置的位置。复选框“
Show crosses
”可使在拾取的端点
绘制一小的倾斜十字，
2More advanced: Rectangular pick.
矩形拾取。简单拾取存在最大值、最小值拾取的问题，矩形拾取可克服这一问题。
使用“
Edit/Rectangular pick
.”，光标改变形状，将鼠标移到6,12拖动鼠标到10,8释放，
搜索到最大场值；
移动鼠标并观察状态栏，状态栏出现“Second”，将鼠标移到6,14拖动鼠标到10,10释放，
搜索到最小场值，同时最大和最小场值位置连成一条线，即为拾取线。
3Removing the pick.
用菜单“
dit/Delete pick
”。
4Saving and loading of pick file.
用菜单
File/Save pick
，以“”名保存拾取。再次启动MAGPICK，装载
文件，用
File/Open pick
重新装载拾取。
-52-
“”文件，每一字串为一个拾取，
8.50 9.50 151.79 8.00 11.50 -55.41 346.0 2.06 8.25 10.50 0.000
-31.67
最后一个数总是0.000，在0.000前的两个数是拾取的线段中心X、Y。
5Notes about picking.
许多局部磁性体接近于球体磁场(偶极)，这个简单的磁性源有6个参数：X，Y，Z，Jx，
Jy，Jz，不仅考虑感磁，同时考虑剩磁，MAGPICK认为磁场表示为任意方向的磁化J。利用IRGF
估算本地地磁方向和数值。
6Pick export.
12.5Finding magnetic dipoles
MAGPICK估算全场磁矩，基于剩磁化强度Jrem和磁化系数κ，
德国汉堡地磁参数：
Ttotal
= 49011
nT
，
D
= −0.1，
I
= 68.3，
Jrem
= 0.24
cgs
，
0.87
gauss
，
mass
= 226.5
kg
。
=
1Estimation based on magnetic map.
基于磁异常图推测。打开“”文件，可见到典型的偶极异常图。Zoom图形从
4,4到4,16。
-53-
用矩形拾取的拾取，其拾取的中间点为偶极最初的位置。用菜单“
Inverse/Run
”，弹出
对话框(反演参数)：
Man filed inclination
：倾角；
Man filed declination
：偏角；
azimuth of X axis=90，
X轴方位，因X轴为东西向
；
total magnetic field
：全场值；
Susceptibility
：磁化系数；
-54-
Remanent
：剩余磁场；
Margin space
：；
Number of iterations
：反演次数=4000；
Limit to stop
：停止限制；=0；
Initial depth
：起始深度；=5；
Density
和
J scale
：密度和J；先不用设置，
Induced only
和
Linear only
不被选择。
Grid type
为
GS ASCII grd
，
改变位置文件名
Position file:
为；
Calculated grid:
计算网格为；
Difference
：；
对象类型
Object type
：为
Dipole(s)
偶极。
单击开始后，进度栏显示反演次数，误差逐渐减小，当不见小停留于某一数值，则停止，
反演结束。用+和一个数字表示偶极的水平位置。另一种表示形式为工作表，在菜单“
Inverse
-55-
/ Worksheet/ Dipoles
”，可见到反演的数值。
加载，比较计算场值和原场值的差别。
加载，观察剩余异常。
球体真正的位置：X=8.5m，Y=10m，Z=2m，从工作表中比较。
2Stability of the solution.
得到稳定的答案。有两种类型的标准：
检查观测和计算磁场的差别；
改变初始的位置和深度，看反演结果是否一致。如下表：
删除拾取，改变模式为简单拾取，在另外的地方拾取，启动反演。反演后，观察结果是
否趋于稳定。注意：初始深度不能为0。
3A solution based on profile data.
-56-
基于剖面数据的答案。用插值法形成的磁场图反演，有时不准确，现在用剖面计算：
装载“”文件；
用“文件“菜单“
Profiles/Load list
”装载“bomb_”，19条剖面被装载；
将窗体图形zoom从4,4到14,16；
新建一窗体，选择所有剖面，注意所有剖面线在窗体上变为红色，这个窗体有反演菜单；
设置剖面参数：
Parameters/Settings
，type设置为along Y，OK后出现各剖面曲线；
在菜单“
File/Browser
”，可见到所有选择的线，单击“
Deselect All
.”(取消所有选择)，
此时窗体为空，选择一条或多条线，曲线随之出现；
单击“
Deselect All
.”(取消所有选择)，回到图形窗体，用“
Edit/Pick profile
”拾取几
条穿过异常的剖面，在剖面窗体可见到曲线；
在异常图上，到菜单“
Edit/Direction
”，在图形点左键，以大约45度方向向北东移动，
再次点左键，出现箭头。回到剖面窗体，在“
Parameters/Settings
”，
Type
设为
along pointer
，
则剖面图上数据为设定方向上的数据，起点即图形窗体○符号点。
重新设置为along Y。
在许多情况下,需要裁减数据，图形窗体使用ZOOM方法，剖面窗体上有更灵活的方法，
在图形窗体选择菜单“
Edit/Polygon
”，画一多边形，
-57-
在图形窗体做一个拾取，作为起始位置，然后在剖面窗体反演
Inverse/Run
，在对话框中
改变文件
Positions file
Output file
，启动反演。在图形窗体见到×数字的偶极位置。工作表
中，P为剖面反演结果，剖面曲线理论值为虚线。
深度大约0.6m(提示：顶部、底部传感器的分离)，当起始深度<0.6m，反演结果不趋于
一致，起始深度应比推测的深度大。
基于剖面的反演比图形反演稳定性差，取决于剖面的位置。
4Worksheet operation
12.6Finding magnetic pipes.
模型的长度要比管的直径大许多。
1Start MagPick and load all information.
-58-
版权声明：本文标题：MagPick使用手册 内容由热心网友自发贡献，该文观点仅代表作者本人， 转载请联系作者并注明出处：https://m.elefans.com/dianzi/1721117475a861134.html， 本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌抄袭侵权/违法违规的内容，一经查实，本站将立刻删除。
发表评论