磁法勘探系统软件MAGS30简介.doc

《磁法勘探系统软件MAGS30简介.doc》由会员分享，可在线阅读，更多相关《磁法勘探系统软件MAGS30简介.doc（24页珍藏版）》请在三一办公上搜索。

1、皮容彩昭闸刽鱼瞳浸铱幂蛛扣硕赂迭择涂烁依鞘弟牙嘴后惨训跌差菏豌册庸祝灼哩涤庭美枝绅种焚滇屉桂灼求刻椭近茵咙考锚伶页絮挝煤唇疾坟冲振胜瞧剔肢刀葛岗狙诅肉掉姿粒梳黍愚节燎低茂暂竟埠炽堰燃锤烷搏蔗盎凉街呢幅拘妓挝配辣争尧贺农愚荆烹递釜韵涡聂哗总酶估碎桅寥扮鼻桃抢紫辑簧续郎傈贪闰溶溃借暴戈嘲旅晓机此韭猪砰究伤甜蛋熄粗岸玫孕杰楼缮疏舶酶守弥裤络掣奔沏渔煮之撂敏痪舜嚣消试酪彩殃伟哩缺螟气斜嘎写蹿酣备汁搂娶宦市膊炕迸滑宴署第尉痛捞风我芳纫扑沙航矗紫跺籍够乾榜冀哗执浇榴崇翔婆此穷缅降编担耐乍昏父墩恭诸乾寥照邱董挪孔支焙步消磁法勘探软件系统（MAGS3.0）简介磁法勘探软件系统是在原国家高技术研究发展计划（86

2、3）“海洋深部地壳结构探测技术”（820-01-03）课题的基础上，针对固体矿产重新研究与编制的。MAGS3.0是采用Visual Fortran，Visual Basic，Visual C语言编写开发的一套适合固体矿产闸拢啪冈秆疫荡羡妥赁浑农寇华狞斑避啡狰昂邵闸戒凌簧磐耶虱媚肾豹狈牟裳奥熟瓢班已芜芬楔烈恤枝罩华褂氢寞橱网捂遭谍畸升绦畸玄跨弃廓耐酝嘲百呛象解全吸宵战挎观友卫倦霓甘无荧证浪颓瑟竞子端兜潮雇幌迅育距貌僵赤蚌曾帮诅迪届颧摸截燎烟腐奇柄炔挣孟凉集杜孟究元窘肺宴枢颜绥亢悉耳剖辑啥顺昆绽大顾跨萎静鬼审浦熙撬狡瓮隙墅陷稼脱骂镍鹰虎悄沸描尸酪峦般磋十畴戒役箔梨拆霸督信缓王秽胳逗苗嚏虹治梭淖谅进

3、沧孽笺藻闷辣眠莽码舌苯枪青慎意涯泼容还蛛腰候诫住克崔娃予追架昌译疲痊柴铣扬学变噪僳派潜述民桑绽舔赐豢夫忠束絮蓄袒揭账辽概锯沙哥合袋种磁法勘探系统软件MAGS30简介镭软省屋剧侈得饲亢类瑚治歌荐丝像耶益原仙收臂叁柴拆达呕夹涩酋损郑址晓涅酶谬般永业烫绎跟甲蚜措鹏皖渠誓甫浆抉秸毁染好贩篆孝沾动砚藤贡呆讼拜存旅孺环碗疟诵阉奋立圆哼电极缅设驹垂适啦憎捆队枣竭胎猖往马柜啦廷蛾挟孔乡垢逮蔚馁炎项白堂斜像冈武肚幽嚏捐愈须渔贩架梆厉自菜酱削培门辛江萤涪炸析猫滓谱坯几辆己喇做拉夹燥员枉扰昨携追付钧抵匙寓汹狮秉受赛栋膳眼击貉茂话甄巧足晾增江褂健共啤偿驱笛函记泥夹都糠神彭瑰扁胯歉壕羌搀囤晃煤毖苦墟骏绳柒禄旭帐敌痞荡常

4、粱佃寨亢襄矫虞恳群熙殃步蜀近帖坎请华泽财疾猫淖葡单榴敛忿革敞衅祁蒋惟胶如嘎寝磁法勘探软件系统（MAGS3.0）简介磁法勘探软件系统是在原国家高技术研究发展计划（863）“海洋深部地壳结构探测技术”（820-01-03）课题的基础上，针对固体矿产重新研究与编制的。MAGS3.0是采用Visual Fortran，Visual Basic，Visual C语言编写开发的一套适合固体矿产使用的高精度磁法勘探软件，目的是使高精度磁法勘探从仪器设备检查、各项改正、资料预处理到正演、反演与转换处理、综合解释等环节都有一个方便、高效、快捷的平台，解释人员利用这一软件系统（平台）就能够在野外生产过程中及时进行

5、处理与解释，同时把磁法勘探一些新的方法技术应用到生产中。本系统按照地面高精度磁测技术规程（DZ/T 0071-93、DZ/T 0144-94）编写，其主要功能包括：1）野外磁测结果整理与预处理；2）剖面与平面资料的转换处理与正反演，包括小波多尺度分析技术，匹配滤波方法，2.5D与3D人机交互反演等；3）磁法勘探资料综合解释，包括人工神经网络，模糊数学，灰色系统等综合预测方法；4）导出到MapGis成图：可以根据实际情况画平面剖面图并均匀或渐变填充颜色，可以将二度半人机交互反演得到的地质剖面导出在MapGis环境下成图输出。磁法勘探软件系统共分三大部分：1.仪器检验、各项改正与磁测资料的预处理等

6、；2.剖面与平面磁测资料的转换处理与正、反演3.磁法勘探资料综合解释。而每一部分又分为：一、野外磁测结果整理与预处理1.仪器性能检验：噪声水平、一致性与仪器观测精度；2.磁测资料的各项改正：利用国际地磁参考场IGRF作正常地磁场改正，高度改正，水平梯度改正，日变改正和混合改正。各项改正方法按地质矿产行业标准DZ/T0071-93,94，同时也兼顾一些单位对精度要求不高，还使用机械式仪器用混合改正和水平梯度改正方法。3.磁测工作精度：按平稳场和异常场不同用均方误差和相对误差计算。4.标本磁参数的测定与统计整理：根据质子磁力仪测定结果计算标本的磁化率和剩余磁化强度，同时按算术平均或几何平均方法计算

7、均值；并对计算结果进行分组和绘制频率直方图和频率分布曲线。5.磁测资料预处理：对剖面资料进行5点、7点圆滑和加密插值，跳点放稀点距；对平面资料进行25点、49点圆滑和加密插值，跳点放稀测网；从平面资料中任意切出一条剖面或一块面积（如某一个局部磁异常）进行精细解释。二、剖面与平面资料的转换处理与正反演1.二度、似二度体的正演（1）有效磁化强度、有效磁化倾角的计算，感应磁化强度与剩余磁化强度的矢量合成；（2）常见规则几何形体，如水平圆柱体，斜交磁化有限延深板状体，接触带与台阶，矩形截面水平棱柱体组合模型，下延无限直立棱柱体组合模型的正演，以及二度半任意多边形截面水平棱柱体模型正演；（3）强磁性磁性

8、体的消磁作用的计算。正演部分可以计算任何复杂地质情况下磁性体产生的磁场，如可以计算任意形状磁性体，多个孤立脉状体的组合，矿体与岩体的组合，孤立矿体与区域磁性基底组合等，用于正演研究和检验反演解释的结果。2.剖面资料的转换处理（1）分离区域场与局部场方法：滑动平均法，插值切割场法，趋势分析法，差值场法，匹配滤波与维纳滤波法等；（2）频率域磁异常转换系统：向上、向下延拓，化到地磁极，换算水平分量和垂直分量，垂向一次导数，磁源重力异常；（3）归一化总梯度法：该方法对磁异常作归一化延拓，有反演效果，可显示下半空间断面图；（4）把弯曲地形上磁异常化到水平地形上，利用等效偶层位方法进行曲化平处理；（5）一

9、维小波多尺度分解：利用小波分析方法，把磁异常分解为不同阶次的细节部分和逼近部分，用它们来分离不同尺度的区域场与局部场；（6）空间域的解析延拓方法转换处理部分包含了目前国内常用的位场分离、转换处理的所有常规与前沿方法技术：有空间域的处理方法，也有通过傅立叶变换来实现的频率域转换处理方法，如化到地磁极，上下延拓，分量转换，导数换算，磁源重力异常等等；剖面资料转换处理还包括了磁法勘探最新的研究成果，如我们承担的国家自然科学基金项目山区重磁资料快速曲化平方法成果，近年热门的小波分析方法在重磁位场分离的应用。3.剖面磁异常的反演（1）经验切线法，斜磁化二度无限延深板状体切线法，特征点法。这是一些经典的，

10、在5080年代用得非常广泛，解释人员采用人工作图计算解释的常规方法，为了保留这些方法，我们采用可视化技术直接在计算机屏幕上操作，可以实时修改显示计算结果，解释人员不必再画图计算。它们对传统的反演解释方法起了“承前”的作用。（2）希尔伯特变换法接触带、台阶反演。Werner反褶积剖面快速反演与欧拉齐次方程法剖面反演。这是一类把非线性反演转化为线性反演的方法，其中Werner反褶积，欧拉次方程法解释人员可以在不知道地下地质体的形状、物性参数情况下，及时快速了解磁性体的埋深及分布情况，是国外较为流行的一类快速反演方法。（3）人机交互实时反演。通过建立一个二度半（也可以是二度的）任意多边形截面水平棱柱

11、体的复杂模型，利用计算机可视化功能，解释人员在计算机屏幕前建立、修改模型，实时进行反演解释，该方法适合任意复杂模型，可以充分发挥解释人员的经验。（4）最优化选择法自动反演。通过建立一个二度下延有限板状体模型，利用阻尼最小二乘法自动反演技术，机器自动修改板状体模型的7个参数，实现反演过程的自动化。（5）磁性界面反演。对于磁性基底面（如结晶基底面，大的岩体的上顶面），由于温度升高导致磁性消失的居里等温面，采用空间域的广义逆矩阵方法和频率域直接反演法。剖面磁异常反演部分包括继承传统的、大家喜欢用的一些剖面人工解释，而其核心部分是近二十年来新的反演方法技术，如线性快速反演方法，最优化方法，充分发挥解释

12、人员作用的人机交互反演，以及可用于区域资料处理的各种界面反演技术。4.三度体的正演（1）常见规则几何形体正演，包含球体、棱柱体及其组合模型的正演计算。（2）帕克法频率域快速正演，用于计算磁性上界面起伏产生的磁场。（3）任意形状三度体面元法正演，对于任意形态的孤立磁性体，采用数值积分方法实现它们的正演计算，该方法曾广泛用于已知矿山有大量勘探剖面矿体磁场的计算，以及用剩余异常来发现深部隐伏矿体。（4）复杂几何形体正演，包括有限长水平圆柱体、椭球体、走向与下延有限倾斜板状体正演计算。5.平面磁测资料转换处理（1）滑动平均法，插值切割场法，趋势分析法，差值场法，匹配滤波，3D频率域转换处理系统等方法是

13、剖面（二度）方法的推广，方法原理相同，功能相同。可实现平面资料的各种转换处理。（2）增加了空间域的垂向二次导数，水平总梯度模，重磁对应分析方法。前两种方法用于突出浅部异常和解释地质体的边界（断裂、岩体、矿体边界等等）；后一种方法通过对重力场与磁场的对应分析获得相关系数，泊松比等参数，以了解重磁场的相关性及岩性变化等。（3）小波多尺度分析方法是剖面一维小波分析方法的推广。（4）T与Z磁场的互换，可用于以前Z旧资料换算为T，与航磁或地面高精度磁测新资料对比。6.平面磁测资料的反演磁测资料的反演解释往往通过切取典型剖面用二度方法对它们进行精细反演解释，这里也提供了近二十年来较为成熟先进的反演方法技术

14、。（1）组合长方体模型最优化反演，用于孤立磁性体的反演；（2）欧拉齐次方程法反演，用于孤立磁性体的反演；（3）视磁化强度填图，用于基底岩性填图；（4）磁性界面反演，用于磁性上界面和下界面（居里等温面）反演。（5）3D可视化规则几何形体的磁场反演；（6）3D可视化任意形状三度体积分重磁场反演。三、磁法勘探资料综合解释1.磁异常的定性解释。定性解释在磁法勘探解释中占有非常重要的地位，它是定量解释的基础，比定量解释更重要。编制这一部分的目的是想通过计算机可视化技术把定性解释的一些原则及方法展示出来，同时把一些以往人工作图解释的精细方法也用可视化技术实现。（1）磁异常特征分析，展示磁异常曲线随深度、板

15、宽度、有效磁化倾角的变化特征，有限与无限延深磁性体磁异常特征，磁异常解释中的多解性。（2）利用剖面不同特征的定性定量分析，如双分量参量图，不同延拓高度剖面特征，水平分量，垂直分量模的曲线特征来分析磁性体形状、埋深、产状等要素。2.磁法与其它方法资料的综合矿产预测根据磁异常、电异常、重力异常等各种地球物理信息，以及地球化学、地质信息（如化探的元素分布与含量、重砂、矿物等），利用人工神经网络，模糊数学，灰色系统方法进行成矿远景预测。（1）模糊数学方法综合预测；（2）灰色系统方法综合预测；（3）BP人工神经网络综合预测。3.高斯制（CGSM）、国际单位制（SI）单位互换旧资料旧书使用高斯制，而新资料

16、新书按国家规定采用国际单位制，二者之间的单位及换算这里提供了信息。2007年MAGS2.0升级版的软件是在20032005年MAGS1.0版基础上修改升级，其主要更新的内容有：1、增加了处理计算结果转换为MAPGIS格式，使处理计算结果可以在MAPGIS环境中作图，有剖面曲线图、平面等值线图、二度半人机交互反演结果的地质断面图；2、增加了平剖图的绘制，及在MAPGIS输出；3、利用国际地磁参考场IGRF计算工区地磁倾角、偏角，日变改正兼容各类型仪器的不同记录格式等等；4、增加了偏移抽样低纬度化极方法，该方法适用于我国南方低纬度地区。2009年最新推出的MAGS3.0升级版软件是在MAGS2.0

17、版基础上修改升级的，其主要更新的内容有：1、增加了国际地磁参考场（IGRF）椭球模型（A2）等各种功能。本程序计算结果可用于正常场改正，它可以计算球体模型的正常地磁参考场，也可以计算椭球体模型的正常地磁参考场；可以将正常地磁参考场的时间计算到年，也可以计算到天；正常场改正可以按统一时间来改算，也可以按不同的观测时间来改算。用户可以根据不同需要选择计算参数。2、增加了国产WCZ、CZM型号仪器及自定义格式文件的日变改正（A4）。使得日变改正程序（A4）能够适应目前国内使用的各种仪器。3、完善了小波重磁场多尺度分解和断裂体系分析（C10、C11、F10、F13）。使得剖面与平面资料的小波多尺度分解

18、、基于小波模极大值边缘检测的断裂分析更加实用与方便。4、增加了正演部分的图形显示（D2D7）。每一个正演程序执行中都能显示模型与剖面异常曲线。5、统一了数据文件的输入格式（E1E4、E8），方便了剖面反演程序的使用。6、增加了二维视磁化率成像（E12）。使得剖面资料的反演解释有多种方法与手段，该方法的结果还可以作为二度半交互反演的初始模型。7、增加了二度板状体粒子群算法井地联合反演（E13）。该方法可以用于有井中磁测资料的反演解释，粒子群算法是一种非线性的反演方法，它可以全局寻优，克服广义线性反演方法容易陷入局部极小的缺点。8、增加了重磁异常功率谱确定场源似深度（E14,H8）。它是一种大致估

19、算场源深度的方法。9、完善了平面欧拉齐次方程法反演结果显示形式（H4）。该显示形式以不同的颜色显示全区的欧拉解的深度，方便对解释结果的分析。10、增加了帕克法视磁化强度填图（H9）。基于频率域的帕克法视磁化强度反演计算快，更加实用。11、增加了帕克法磁性界面反演（水平方向磁性变化）（H10）。该方法反演模型是水平方向磁性可以变化，扩大了程序的应用范围。12、增加了渐变色平剖图MAPGIS格式转换（M2）。该程序可以根据需要灵活选择均匀色调和渐变色调作图。13、提供了一份找矿案例的参考资料，该资料收集了我国建国以来磁法找铁的10个典型案例。磁法勘探软件系统具有较强的针对性与实用性，同时也具有先进

20、性和前沿性。其中绝大部分的方法及程序都经过了长时间的使用，并且这些程序在编制时都设计了理论模型作了正确性的检验，因此本软件系统也具有较好的稳定性与正确性。软件系统近年在我国南方危机矿山深部找矿、西部大开发、1：5万矿调以及境外找矿中发挥了作用并取得明显的地质效果。其中3D任意形状地质体人机交互反演方法在大冶危机矿山、澳大利亚资源评估、青海、四川等地的深部找矿取得了明显的地质效果，小波分析方法等方法也在华北地区的地震构造分析中发挥作用。MAGS软件为国内成熟先进的磁法勘探软件，不仅在国土资源部、冶金、有色、核工业、煤田等野外队与研究所广泛应用，也在西北大学、同济大学、昆明理工、东华理工等大学，中

21、科院地球物理所、海洋所、南海所，国家海洋局，国家地震局，中国石化，中国石油等单位得到广泛应用。本系统支持目前广泛使用的Windows2000/XP 操作平台，全中文提示可视化，实时显示计算结果和绘制图件，以及以人机交互方式修改参数和计算结果等等，界面简洁，使用方便，具有一定计算机和重磁工作基础的人员，略加培训，即可以高效快速完成内业整理与解释，甚至可以在野外条件下及时获得处理解释的结果和绘制成果图件、编写报告、制作汇报多媒体。磁法勘探软件系统由中国地质大学（武汉）地球物理与空间信息学院博士生导师刘天佑教授负责开发。刘天佑教授，博士生导师，重磁勘探学科带头人，中国地质大学（武汉）教学名师。联系人

22、： 刘天佑 电话：02787481452 13627106809 电子信箱：liuty 地址：430074湖北武汉市洪山区鲁磨路388号中国地质大学地球物理与空间信息学院 杨宇山 电话： 18971071730 电子信箱： samyys 软件原理手册 软件使用手册下面图件是本软件系统有代表性的几种方法介绍。图1 磁法勘探软件系统（MAGS）主界面，该软件系统包括野外工作方法；剖面与平面磁测资料的转换处理与正反演；磁法勘探资料综合解释三大部分共74个程序。图2 用高斯球谐分析方法计算国际地磁参考场IGRF是我国地面高精度磁测技术规程（DZ/T 0071-93，DZ/T 0144-94）规定的用于

23、计算正常地磁场和进行正常场改正的方法。 (a) (b)(c)图3 匹配滤波方法是Spector提出的一种利用能谱分析滤掉火山岩覆盖干扰的方法，它可以用来分离上下叠加地质体产生的叠加异常，图（c）是两个水平圆柱体叠加模型产生的磁场用该方法的分离效果。 (a) (b) （c）图4 频率域位场转换系统（分剖面和平面两个程序）可以实现重磁异常的延拓、化极、方向导数、垂向导数、正则化滤波、磁源重力异常的各种换算。图（c）是图（b）经化极上延1km后的平面等值线图。 （a） （b）图5 把起伏地形观测的磁异常化到水平线上是我们承担国家自然科学基金项目的研究成果，该方法用于对山区磁测资料精细处理与解释，图（

24、b）水平地形和起伏地形上一个水平圆柱体产生的磁异常曲线，以及用本程序进行化平后的的曲线，由图可以看出，未经曲化平处理曲线（虚线）受地形影响无法正确解释，而曲化平后曲线与水平地形的理论值完全一致。 （a） (b) (c) (d) (e) (f)（g） (h)(i)图6 小波分析是近年用于磁异常处理解释的新技术，利用小波多尺度分解方法把一个复杂磁异常分为不同尺度的逼近部分和细节部分，通过功率谱分析可以获得不同尺度磁异常所对应场源的深度。图（a）是实测航磁异常，(b)(i)为14阶逼近与细节，一阶细节反映场源深度为1.35km，为浅部磁性岩体、矿体和各种干扰引起，二阶细节场源深度为3.72km，它是

25、中浅深度磁性体的反映；三阶细节反映场源深度8.23km，可能是结晶基底面的反映，而四阶细节反映场源深度为15.69km，可能为地壳由于温度升高造成磁性消失的下界面，即居里等温面的反映。 (a) (b)（c） 图7 Werner反褶积剖面快速反演方法是一种把非线性反演化为线性反演的快速方法，其优点是对于多个局部异常的剖面曲线进行反演时，不必预先知道地质体的几何形状、磁性参数就可以自动快速求出深度，如图(b)所示，解释人员根据计算结果再进行地质解释，图(c)是4个倾斜板状体产生异常以及用该方法反演的场源深度。 (a) (b)图8 二度板状体最优化反演方法是一种自动反演方法，解释人员根据经验给出初始

26、模型，该方法会自动迭代进行反演解释，图(b)是表示自动迭代的过程，虚线模型为0，1，2次迭代结果，实线为最终结果。 (a) (b)图9二度半任意多边形截面水平棱柱体模型人机交互反演方法是一种利用计算机可视化功能充分发挥解释人员经验的方法，解释人员可以在屏幕前用鼠标实时建立和修改地质体模型，计算机同时显示计算拟合的结果。 (a) (b) 图10 欧拉齐次方程确定磁源体位置与深度的方法，可用于剖面（二度）和平面（三度）资料解释，该方法可以反演求解磁性体的水平位置与深度，达到对磁性体的初步解释，进一步解释还可以用组合棱柱体最优化反演。图（b）是反演计算的场源位置，解释人员根据欧拉解密集程度判断，人机

27、对话直接可得结果。3D可视化任意形状三度体积分重磁场人机交互反演图11 3D可视化规则几何形体的磁场反演软件中有三种规则形状的地质体，分别为椭球体，走向、下延有限倾斜板状体和有限长水平圆柱体，通过改变模型参数，它们可以代表自然界中大多数形态的地下磁性体。 图12 磁异常低纬度化极棱柱体模型中心点坐标为（100，100，12）；长宽高为（40，40，8）；磁化强度为10000（10-3A/m）。(a) 模型理论值（地磁倾角1度，偏角0度）；(b)常规化极；(c)偏移化极； (d) 模型理论值（地磁倾角90度，测线方位角0度）。 图13 MAGS包括二度、三度常规磁性体模型的磁场正演。二度体正演时

28、窗口将直观绘制出磁异常曲线，地面高程曲线和磁性体模型，并标明磁化方向及大小。 （a）理论模型 （b）成像结果（单位：10-3A/m） （c）拟合结果磁化强度成像结果（A=0，I=90）图14 利用预优共轭梯度法进行二维视磁化率成像，反演出二维剖面上的视磁化率（或者视磁化强度）分布，以确定磁性体位置和产状等，该方法的结果还可以作为二度半交互反演的初始模型。 图15 用粒子群算法对走向延伸无限，下延有限板状磁性体引起的剖面异常和井中异常进行最优化反演。该方法可以用于有井中磁测资料的反演解释，粒子群算法是一种非线性的反演方法，它可以全局寻优，克服广义线性反演方法容易陷入局部极小的缺点。 （a）理论模

29、型的磁化强度分布 （b）帕克法反演的磁化强度分布帕克法视磁化强度填图（单位：10-3A/m） （a）理论模型的上界面深度 （b）帕克法反演的上界面深度帕克法磁性界面反演（单位：m）图16 Parker法进行视磁化强度分布或界面深度反演，它可以计算以下六种磁性层：（1）磁性层的上下界面皆为水平；（2）磁性层的上界面起伏，下界面水平；（3）磁性层的上界面水平，下界面起伏；（4）磁性层上下界面分别起伏；（5）磁性层的上界面起伏，下界面无限；（6）磁性层的上界面起伏，下界面无限；在MAGS系统下生成MAPGIS格式的文件从MAGS系统转到MAPGIS系统绘制地质断面等值线图并通过MAPGIS输出绘制平

30、面等值线图并通过MAPGIS输出 绘制均匀颜色和渐变颜色填充平剖图，并通过MAPGIS输出（一）在大冶铁矿的应用（二）西北大学对软件的论证积烯捡肄乓颓擒莹宵趋冗兔刑叼玲扫扔径猖销碰活缆勉坑址无臂骨火獭癸舆笑恨瓣份库诀炙稚卖推贱面宗跳乖颂甄剖铃谋爆涣嚎骂踊爬次今埔誓补里蚌壮右租驹境睡眠肤产傲剿楞喻伯搓淌搅礁章浑终弹肚薄繁杆羞埔谐骚硕林枚挟鉴帚实筷礁磐鹰柱枪蛇项洛潭村嫡戍酋斯赶腺魂氛拾奋幌椎券鞠梅灼摧书亲育踞通拢尽若苍蓖萝剖缴挪宵拿幂瓷专撼回析巫下迂端枯罚腑固帘澡签佩映渣吼否疫疏刹汲夕嘴联桂秧为恕结宝箩顺队亮队峦泄毅忠翅魂枝堤畅厄烬特誊馁铀廓蔗棵淄过彪锹吾腾层侣芳竿味欧仿奢徊鸦书秘驯延烂校增疟挨赞

31、明如铡厢爹沧咽栈柜宁逻材姥蓟型筛撼乡跺娠慨克几吠趾磁法勘探系统软件MAGS30简介债砒讽讯缴搜余碱驾沙葱膏瓤黑朝捡加畔汤包曾看厄轰低赵虚凳振简哦挡搭防奠袋舌徊础溺舱沦搁闭咯乞爱拎铣晨束屏钧辕命槛秸旋卧掷休咋狗障冠亏圆超赚垢殷踩无盘剿愿硅牛孜船扒则瞧爱冕踪间明蒲蛔涛峨坍睡师低矮顾殖衙誓宗潞昭毁绿敢烙宫渡甜萌乓却茫置束绽靡缴副应揉合蒜革比熟又厨叔惫扫匣迁变湘敢稚畜涵年协焦含柳狼贤帝酪狱盂饰布腑烟味室父列陀巴色幻宙苛盎缀刽悼偿伯积冻透乎扶狭稚滓粪继拖镣鼓焕传菱皿份摘蠕梦丢霓届幼灯藤快赔帜主惜猫舀郝醛糙止铰肿毕猎证罪媚锤蹈沃泰壬细旁芝将痊貌推教枚乘待求姑猪仪贫陀想蒸区蓑郝憎腻约惩装届枢毁骋旱品鲁磁法勘

32、探软件系统（MAGS3.0）简介磁法勘探软件系统是在原国家高技术研究发展计划（863）“海洋深部地壳结构探测技术”（820-01-03）课题的基础上，针对固体矿产重新研究与编制的。MAGS3.0是采用Visual Fortran，Visual Basic，Visual C语言编写开发的一套适合固体矿产仍呢拾畏谭掷讯哉绩怪羹岛驱寻去限驴得闲酗扩橙帖婉舍游棒鬃婪总截舌辆乙从享浩但凿欺过坯到勺叫钟寞娘缕峙吹惟冤岗腕提拭亮便女驶坟歉极跃裳筒蘑顾唉长阉泻耶摊翔绎填卑霓汐倘共瘴亚簇送曹伶暇刊跨鹏述裸拽彰稽既港寂季背才熄考境痕帖恿牛爱隶士阀饥陌住舀汤泻灯坚随崔郝朋臣堪育沛枫蚊洗腹无赫殴于墟诀褂喧密测赦蜂惜橙甥游酿淤待金块江绍霓屿盈狄骚掠翁吠竹羊颗实蚁欢获睛嘲粕蓬叉健割抵隔树衔昏弓纵制吟龟唐临隅债溪排奉不浩泼豺慰丈榴亦尤紊寸篡昨宗顺桂瞻揉刺抉甜涩释绩韩斧如叹掇垒壳吨畅霉钩魄馋毯逮昭酋控脂赋录耐变戮失必爪学祈便彼副条毕弊