地磁作业(转换、延拓)
用解析延拓求取剩余磁异常来圈定浅部隐伏磁铁矿脉
用解析延拓求取剩余磁异常来圈定浅部隐伏磁铁矿脉李建华;李春鹏【摘要】Residual anomaly calculation is one high-pass filtering method by extracting shallow media induced anomaly by pole reduction and continuation. For magnetite veins surrounded by rocks with strong magnetic, the anomaly induced by shallow magnetite veins always covered by huge magnetic surrounding rock. Magnetic field strength is inversely proportional to the square (two-dimensional) or cubic(three-dimensional) of the distance between the station and the source. The high-frequency anomalies induced by shallow magnetic sources fade faster than low-frequency anomalies induced by deep sources. Therefore, when the field with analytic upward continuation is removed from the measured magnetic field, the residual anomalies can reflect shallow magnetic bodies with high-frequency anomalies. The residual anomaly is calculated by analytical continuation to extract the anomaly induced by shallow magnetite veins for exploring insidious magnetite veins in this paper. Furthermore, several insidious magnetite veins are found by this method on Morritos magnetite deposit in republic of Peru. It is validated to use this method exploring insidious magnetite veins surrounding huge magnetic rock.%剩余磁异常求取是一种对磁场强度进行化极,解析延拓后提取浅部磁性介质异常的高通滤波方法.对于处于强磁性围岩的磁铁矿脉来说,浅部磁铁矿脉的异常往往被大规模的磁性围岩所引起的异常所掩盖.磁场强度与场源到测点距离的平方(二度体)或者立方(三度体)成反比,进行向上延拓的时候,由浅部场源引起的高频异常衰减快,深部场源引起的低频异常衰减的慢,因此地面实测磁场去掉向上解析延拓后的磁场值,可突出反映近地表磁性体引起的高频异常.这里通过解析延拓求取剩余异常,实现浅部磁铁矿脉异常信息的提取,为寻找隐伏的磁铁矿脉提供依据.并且利用此方法在秘鲁共和国Morritos铁矿圈定隐伏磁铁矿脉,验证了剩余异常在大规模的高磁性围岩中寻找隐伏磁铁矿脉的应用效果.【期刊名称】《物探化探计算技术》【年(卷),期】2012(034)003【总页数】7页(P288-294)【关键词】剩余磁异常;解析延拓;隐伏磁铁矿脉;秘鲁;Morritos铁矿【作者】李建华;李春鹏【作者单位】中南大学地球科学与信息物理学院,湖南长沙410083;中南大学地球科学与信息物理学院,湖南长沙410083【正文语种】中文【中图分类】P631.2磁异常的各种处理转换方法，实际上就是各种滤波方法及其组合。
磁法勘探及应用实例2 (1)
3）磁性围岩的影响 矽卡岩型矿床产于岩浆岩和沉积岩的接触带内。 一般来说，沉积岩无磁性或磁性微弱，而岩浆岩磁 性较强。作为围岩的岩浆岩的异常与矿体异常相叠 加，往往使矿体异常畸变。对于热液型矿床，形成 矿体时的热液作用可导致围岩蚀变，在近矿围岩中 伴生出铁磁性矿物，这些矿物的磁场也会使矿体异 常发生畸变。 图2.4.3为山东某矿区的Za曲线。矿体位于大理 岩与闪长岩的接触带中，闪长岩的磁化强度约 1000l0-6GSM，整个异常为矿体异常与闪长岩导常 的叠加。
目前在空间域进行向下延拓使用的是近似方法,如多项式插值、级数正则化、调和分析等。 如使用逐次延拓逐次圆滑(即向下延拓一次就对曲线圆滑一次)的方法,常常可以取得较好的效果 。延拓中还应注意延拓深度不能太大，尤其不能超过磁体的埋藏深度。因为当观测面进入磁体 内部后，情况会空得十分复杂，这时若运用只适宜磁体外部的向下延拓公式，将导致极大的错 误。 认真研究通过向下延拓所获得的不同深度异常的特征，有助于判断磁异常性质，还可粗略估 计磁体埋深及其分布范围。
含有误差和随机干扰的异常剖面
2) 磁异常的划分 磁异常的划分指将叠加在一起的区域场与局部场分离开来。 如插值是常用的处理方法。
3) 磁异常数据的网格化 在实际工作中，由于这样或那样的原因,实际测点往往呈不规则分布，但在进行数据处理时， 总是要求数据按规则网格分布。因此就需要由不规则网格上的实际场值换算出规则网格节点上 的场值，这个过程就是数据网格化。显然，数据网格化的问题实际上是插值问题，即用不规则 分布的插值节点上的值来计算规则网格节点上的值． 插值方法很多，但通常采用拉格朗日插值的方式。
图2.4.5a为一个无限延深柱体位于山坡的情况。柱体左 侧的测点都位于其顶面之下，由于顶面负磁极在这些测点 处的总磁场方向都指向它自身，故垂直分量是向上的，使 得曲线左支出现明显的负值。柱体右侧的测点位于顶面之 上，故相对于水平面而言，磁体的埋深在这些点都增大了， 因而负磁极的影响相应地减小，致使曲线右支变缓。若不 了解这一情况,仍按水平面上磁异常的特点进行分析,就会 得出磁体向右倾斜的错误结论。当柱体位于山脊时,曲线 的畸变更加厉害,在正异常的两侧均出现了负值,这就很可 能作出磁体向下廷深很小的错误判断。如地形也具有磁性, 情况就更加复杂。因为磁体异常除受地形起伏的影响外, 还增加了一个地形磁性的影响,更增加了解释的难度。
解析延拓法在山阳磁法数据解释中的应用
解析延拓法在山阳磁法数据解释中的应用孙海龙;吕伟星;陈鑫;郑伟【摘要】依据磁法数据处理中的解析延拓原理,对陕西省商洛市山阳县小河口镇工区的磁法勘探数据分别进行了上、下延拓:0.5、2、8、10 m及200、240、260、280、360 m的向上延拓,2、4、6、8 m的向下延拓.纵观磁异常等值线平面图可以发现,该区存在2处主要高值剩余磁异常区域,异常值分别介于80~170 nT与70~90 nT,规模分别为100 m*350 m(异常4)和230 m*250 m(异常3).结合该范围地质填图成果,分析认为异常4是由位于约250m深处的矽卡岩脉引起的,而浅部干扰导致了异常3.该区解释延拓的应用表明向上延拓达到压制浅部干扰、突出深部异常,向下延拓压制深部干扰、相对突出浅部异常的目的.%Based on the analytical continuation principles in magnetometric data processing, have carried out upward (0.5, 2, 8, 10, and 200m, 240, 260, 280, 360m levels) and downward (2, 4, 6, 8m levels) analytical continuations respectively for magnetic prospecting da-ta from the Xiaohekou work area, Shanyang County, Shangluo City, Shaanxi Province. The comprehensive surveying of magnetic anoma-ly isogram plan can find that in the area have 2 mainhigh residual magnetic anomalies with outliers between 80~170nT and70~90nT and sizes 100m*350m (anomaly 4) and 230m*250m (anomaly 3) respectively. Combined with geological mapping in the area consid-ered that the anomaly 4 is caused by skarn dike at depth about 250m, while the anomaly 3 by shallow interference. The application of analytical continuation on the interpretation in the area has shown that the upward continuation can achieve the aim to suppress shallow interference andhighlight deep anomalies;while the downward continuation suppress deep interference and highlight shallow anoma-lies.【期刊名称】《中国煤炭地质》【年(卷),期】2017(029)002【总页数】7页(P76-82)【关键词】解析延拓;磁异常等值线;磁法勘探【作者】孙海龙;吕伟星;陈鑫;郑伟【作者单位】中煤航测遥感集团有限公司,西安 710199;武警黄金部队第五支队,西安 710100;中水东北勘测设计研究有限责任公司,长春 130061;中煤航测遥感集团有限公司,西安 710199【正文语种】中文【中图分类】P613.2磁法勘探是在地质学、地磁学、岩石的磁性理论以及数学理论的基础上，研究地磁场和相关磁性异常体所引起的磁异常在空间分布特征，从而达到找矿，研究地质构造等目的[1-2]。
地球物理勘查名词术语
中华人民共和国国家标准GB XXXX--XX地球物理勘查名词术语Terms Of geophysical exploration1 主题内容及适用范围本标准规定了地球物理勘查(包括重力勘查、磁勘查、电勘查、地震勘查、测井及核物探)中常用的、主要的、本学科专有的名词术语。
本标准适用于地球物理勘查工作的语言和文字交流。
2 基本术语2.1 地球物理勘查geophysical exploration运用物理学的原理、方法和仪器以研究地质情况或寻查埋藏物的一类勘查。
同义词物探；地球物理勘探：(勘探地球物理；地球物理探矿)注：1.取决于使用场合，该术语可附加后缀“法”或“学”。
2.根据具体情况，可以使用“航空物探”，“海洋物探”，“地面物探”，“地下物探”，“深部物探”，“区域物探”，“工程物探”，“环境物探”，……等术语。
2.2 正常场normal field物理场的相对平稳部分。
2.3 异常anomaly物理场对正常场的偏离。
2.3.1 理论异常theoretical anomaly正演所获得的异常。
同义词计算异常2.4 物性physical properties岩(矿)石或其它探测对象的物理性质。
2.5 异向性系数coefficient of anisotropy描述介质垂直层理(片理、节理等)方向与平行层理方向的物性差异的一种参数。
同义词(各向异性系数;非各向同性系数)2.6 地球物理正演geophysical direct problem根据地质体或其它探测对象的几何参数和物理参数计算地球物理场值。
同义词物探正演2.7 地球物理反演geophysical inversion根据地球物理场值，计算地质体或其它探测对象的几何参数和物性参数。
同义词物探反演国家技术监督局XXXX—XX—XX批准 XXXX—XX—XX实施2.7.1 交互解释interactive modelling直接在计算机的显示器上反复修改地质体或其它探测对象模型的参数，使模型的理论异常逐步逼近于实测异常的一种反演方法。
《地球物理勘探》例题
《地球物理勘探》例题一、填空题1．沉积岩密度值主要取决于岩石中孔隙度大小，干燥的岩石随孔隙度减少密度呈线性__________。
2．按照导电机制可将固体矿物分为三种类型：金属导体、半导体和__________。
3．电子导电矿物或矿石的电阻率随温度增高而上升，但__________岩石的电阻率随温度增高而降低。
4．在岩（矿）石的主要物理性质中，__________的变化范围是最大。
5．由不同地质体接触处由岩石的固相骨架与充满空隙空间的液相接触处的电荷自然产生的电动势的物理- 化学过程称为__________。
6．地球的重力场可分为正常重力场、重力随时间的变化及重力异常三部分，其中地球的正常重力是由赤道向两极逐渐__________。
7. 在重力测量重，由于负地形部分相对水准面缺少一部分物质，空缺物质产生的引力可以认为是负值，其垂直分量也是向上的，使仪器读数__________。
8．电法勘探是根据所测得的地下__________________的分布规律来查明地下地质构造和寻找有用矿产的一种常用物探方法。
9．电阻率法是__________类电法勘探方法之一，它是建立在地壳中各种岩矿石之间具有导电性差异的基础上，通过观测和研究与这些差异有关的天然电场或人工电场的分布规律，达到查明地下地质构造或寻找矿产资源之目的。
10．根据地质任务的不同，重力勘探可分为预查、普查、详查和__________四个阶段。
二、名词解释1、磁化率2、视电阻率3、有效磁化强度4、二度体异常三、简答题1、简述地球物理学的组成及研究内容2、什么是地球物理正演（正问题）和地球物理反演（反问题）?3、什么是重磁场的解析延拓，向上和向下延拓分别有什么作用？4、地磁要素主要有哪几部分组成？请画图说明。
5、影响视电阻率的主要因素有哪些？四、分析说明题1、如下图（a）和（b）所示两种剖面的地质条件，请分别说明沿地面能否观测到有相对变化的重力异常？为什么？（a）（b）2、（a）定性画出下面地质体的△ g 曲线，（b）定性画出下面磁性体的 Za 曲线。
磁法延拓在地磁背景场匹配中的应用研究
磁 场值
， ｙ ， ） 为（ ， ， ） 坐标 点待求 的磁 场值 。
２ ． ２ 磁 场延拓 计 算方 法
从式 （ ３ ） 中 可 以看 出 ， ｆ （ ， Ｙ ， ｚ ） 为ｆ （ ， Ｙ ， ０ ） 与 权 函数 关于变 量 （ ， ｙ ） 的 二 维 褶
区域 飞行高度 的地 磁 背景 场 。然后 将 飞 机磁 测 数据
飞机尾 部 的磁探仪 在该 区域探测 目标 时 ， 发现有 些地 磁背景 场 的信号特 征和 目标 的信 号特征 相似 ， 这 类背 景场信 号 也 能 触 发 ＭＡ Ｄ信 号 识 别 输 出 ， 产生虚警。 为 了减 小复杂地磁 背景场导致 的高虚警 率 ， 本文 尝试
使用 海上磁测 数据 ， 通过 向上 延拓 的方法 ， 得 到该
（ ｕ ， １ ２ ， ｚ ） 分别为延拓前 、 后 的频谱 ； （ “ ， ） 为权 函数相
２ 磁 场 延拓原 理与 计算 方法
２ ． １ 磁 场 延 拓 原 理
应 的频率响应 函数 ， 也 称延拓 因子 ， 可以 由下式求得 ：
磁 场 向上延 拓 的定义 ： 已知 高度 ＝ ０ （ 或某 一 固
ห้องสมุดไป่ตู้２
。
基金项 目：国家 ８ ６ ３计划 （ ２ ０ １ ３ Ａ Ａ ８ １ １ ２ ０ ０ ５ ； ２ ０ １ ４ Ａ Ａ ０ ６ Ａ ６ １ １ ） 。 作者简介 ： 王景然 （ １ ９ ８ ４ ． ） ， 女， 河北石家庄人 ， 工程师 ， 硕士 ， 主要 从事磁法勘探数据处理与信号识别相关研究。
磁法勘探06磁异常的处理与换算资料
第一节 磁异常的处理与换算的目的意义
应当指出，磁异常处理和转换时，有两个问题必须要明确： 1.应当合理的选择处理和转换的方法。由于转换、处理方法 较多，具有各自的特点、作用、适用条件，不应盲目的对各 种方法都使用一遍。应当认真分析磁异常特征、测区内地质、 物性情况及所要解决的地质问题，根据各个方法的功能和适 用条件，合理的选择若干种处理方法; 2.磁异常的处理和转换只是一种数学加工处理，它能使资料 中某些信息更加突出和明显。但不能获得在观测数据中不包 含的信息。数学变换只能改变异常的信噪比，而不能提供新 信息；因此，在应用各个方法时必须要注意到实际资料的精
15
第二节 磁异常的处理
1.剖面网格化
16
17
第三节 磁异常的空间转换
延拓是把原观测面的磁异常通过一定的数学方法换算到高 于或低于原观测面上,分为向上延拓与向下延拓。向上延拓 是一种常用的处理方法，它的主要用途是削弱局部干扰异 常，反映深部异常。我们知道，磁场随距离的衰减速度与 具磁性的地质体体积有关。体积大，磁场衰减慢；体积小， 磁场衰减快。对于同样大小的地质体，磁场随距离衰减的 速度与地质体埋深有关。埋深大，磁场衰减慢；埋深小， 磁场衰减快。因此小而浅的地质体磁场比大而深的地质体 磁场随距离衰减要快得多。这样就可以通过向上延拓来压 制局部异常的干扰，反映出深部大的地质体。
是很重要的。随着磁测量精度的不断提高，实测异常中所包含 的可靠信息也不断增加。如何有效地提取和利用这些信息，就 成为磁异常解释理论研究的重要课题。早在20世纪40、50年代， 诸如导数异常的计算，磁场解析延拓，化磁极等处理方法已相 继问世。到60、70年代，由于电子计算机的广泛应用，使磁异 常的处理和转换容易实现，从而其理论和方法得到了迅速的发 展，并不断得到完善。由于在实践中磁异常的转换和处理对提 高磁方法解决问题的能力和改善地质效果起到了应有的作用， 因此它已成为当今磁异常推断解释中不可缺少的重要环节。
地球的磁场及地磁倒转-物理小论文
地球的磁场及地磁倒转摘要：本文首先介绍了地磁场及其成因，进而引出地磁倒转现象，最后总结了对于地磁倒转成因的一些推论。
关键词：地磁场，地磁倒转1、关于地磁场的描述地球的磁性, 是地球内部的物理性质之一。
地球是一个大磁体, 在其周围形成磁场, 即表现出磁力作用的空间, 称作地磁场。
它和一个置于地心的磁偶极子的磁场很近似, 这是地磁场的最基本特性。
地磁场强度很弱, 这是地磁场的另一特性, 在最强的两极其强度不到10-4(T), 平均强度约为0.6x10-4(T), 而它随地点或时间的变化就更小, 因此常用(γ), 即10 -9(T)做为磁场强度单位地球结构，可分成地壳、地幔和地核三大圈层。
地核又分为内地核(固体)和外地核(液态物质)。
液核处于外地核中间部分，它的形成和存在，使其两侧部分具备了独立运动的可能性。
由予地球外壳(即地壳和地幔)自转速度的变化，导致了相对运动的产生。
液核通过液体的粘滞性摩擦阻力实现对其两侧部分的制约作用。
当外壳的自转速度加快时，由于液核的存在，内地核不能同时加速，但在液核内摩擦力的作用下会使内地核转速逐步地加快；反之，当外壳自转减速时，内地核转速又会逐渐地减慢下来。
不论内地核加速还是减速，都滞后于外壳转速的变化，其结果必然在外壳与内地核之间形成一个摩擦面。
由于转速的不同，摩擦面两侧的物质发生相对运动，从而导致了地球磁场的形成。
2、关于地磁倒转的概述古地磁的研究发现，在漫长的地质历史时期中，地球磁场的正负极性不断地发生180。
的来回倒转，并且是全球同步的．这一现象的发现是20世纪古地磁学的重大成果．在地质学中，它不仅可以通过岩石原位剩余磁性的研究及与岩石年龄测定对比，建立起地磁极性倒转——地质年代对应关系，有助于对地层进行更为可靠的时代划分与对比，而且可以据之进一步深入研究地磁场的起因和变化，以及地球的内部结构与自身运动的特点，因而引起国内外学者的极大关注．．C S考克斯于1969年提出了近4．5／Via的地磁极性年表，并被普遍应用，其简化形式如图1图一从图1可以看出，地球磁极倒转有不太规则的长短周期性，长周期约1Ma,短周期约0.2-0.3Ma,中间还有小的反复。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
应用地磁学实验报告实验2——磁异常转换计算学号： 10105218 姓名：朱占升一、实验目的1、掌握水平圆柱体磁场异常分布；2、用Matlab实现水平圆柱体的磁异常场正演计算；3、利用正演结果进行磁异常分量之间的换算；4、通过程序换算认知测点间距即采样点数对换算效果的影响；5、加深对磁法勘探的理解认识；二、程序代码%磁法异常换算%剖面为北向A’=0度，则有I=is。
%所测数据均在同一水平面，柱体深30m，半径8m,测点数为100,%测点间距依次选取2 4 6 8mclcclearfprintf('\n柱体深30m，半径8m,测点数为100,点距依次选取2 4 6 8m\n');for b=2:2:8fprintf('\n点距取%dm\n',b);figure('color','w','NumberTitle','off ','name','za-→ha');x2=1:b:100*b; %点距为bmx1=x2-100*b/2;h=30; %柱深R=8; %柱体半径 ms=pi*(R^2);%柱体截面积k=0.2; %磁铁矿磁化率u=4*pi*10^(-7); %磁导率B=50000; %nT磁感应强度H=B/u ; %磁化场强度M=k*H; %磁化强度m=M*s; %磁矩a=0; %剖面为北向A’=0度.I=90/180*pi; %倾斜角is=atan(tan(I)*csc(pi/2-a));hold onza=u*m*((h.^2-x1.^2)*sin(is)-2*h*x1.*cos(is))./(2*pi*(x1.^2+h.^2).^2);hax=-1*u*m*((h.^2-x1.^2)*cos(is)+2*h*x1.*sin(is))./(2*pi*(x1.^2+h.^2).^2);plot(x2, za,'.-m');plot(x2, hax,'.-g');title('za转换为hax')xlabel(' X剖面走向/m');ylabel('磁异常nT');c=[0.4268 0.1749 0.1103 0.0813 0.06450.0536 0.0458 0.0400 0.0355 0.1759];n=length(c); %转换系数个数m=length(za);for i=(n+1):(m-n)haxz=0;for j=1:nhaxz=haxz+c(j)*(za(i+j)-za(i-j));endhaxzh(i)=haxz;endfor i=1:80haxzh1(i)=haxzh(i+10);endx=((n+1)*b):b:(100*b-n*b);plot(x,haxzh1,'.-')legend('za','hax','za→hax');end三、输出图形及解释通过图形可以知道当点距大于6m是转换出的H a与正演的H a有明显的图形失真（假异常），随着测点间距的加大，失真越厉害。
分析其原因：随着测点间距加大，单位距离内的采样点数逐渐减少，即采样密度逐渐降低。
获取的有效信息不足造成转换后的图形失真，而且推测当采样间距达到一定值之后，转换将会产生不稳定的情况。
测点间距即采样密度与地下的场源密切相关，它与场源体的几何参数、物性参数有着直接的关系。
本次实验中选用了水平圆柱体，其几何物性参数如程序中所设，在这种条件下他的最佳采样点距为2～5m。
当超过这个范围则出现失真、假异常。
四、实验小结本次试验通过Matlab程序实现对水平圆柱体正演后，对其引起的磁异常场进行了分量变换。
通过设定不同的测量点距，认识了采样密度不同的情况下，磁异常场转换失真的特点。
本次实验加深了对空间域磁异常的处理与转换的理解，对公式有了更深入的理解。
对于通过运用计算机实现磁场问题的模拟有了更深刻的理解，也加强了对于模型设置等程序语言的锻炼。
通过实验自己分析问题解决问题的能力不断加强，在实践中知识得到了不断的丰富。
应用地磁学实验报告实验3——磁异常频率域延拓学号： 10105218 姓名：朱占升一、实验目的1、掌握水平圆柱体磁场异常三维分布；2、用Matlab实现水平圆柱体的磁异常场延拓计算；3、利用程序分析延拓效果与高度的关系；4、通过柱体三维延拓了解延拓的失真问题；5、加深对磁法勘探的理解认识；二、程序代码1、水平圆柱体二维向上延拓（浅部异常在右侧）% 10105218 朱占升%剖面为北向A’=0度，则有I=is。
%场源为不同深度的两水平圆柱体，延拓高度选取了50、80、100m三种高度clc;clear;% 测点分布范围b=5; % X方向测点间距X=-200:b:300;nx=length(X);a=0; %剖面为北向A’=0度.%浅部球体1参数I=90/180*pi; %倾斜角i=atan(tan(I)*csc(pi/2-a));%有效磁化倾角isR1=20; % 球体半径 ms1=pi*R1^2;u=4*pi*10^(-7); %磁导率k=0.2;B=50000;H=B/u;M=k*H ; %磁化强度 A/mm1=M*s1; %磁矩D1=30; % 球体埋深 m% 深部球体2参数I=90/180*pi; %倾斜角i=atan(tan(I)*csc(pi/2-a));%有效磁化倾角isR2=50; % 球体半径 mu=4*pi*10^(-7); %磁导率k=0.3;B=50000;H=B/u;M2=k*H ; %磁化强度 A/mm2=M2*s2; %磁矩D2=100; % 球体埋深 m%总Za磁异常x1=X-100;x2=X;Za1=u*m1*((D1.^2-x1.^2)*sin(i)-2*D1*x1.*cos(i))./(2*pi*(x1.^2+D1.^2).^2);Za2=u*m2*((D2.^2-x2.^2)*sin(i)-2*D2*x2.*cos(i))./(2*pi*(x2.^2+D2.^2).^2);Za=Za1+Za2;%总Hax磁异常Hax1=-1*u*m1*((D1.^2-x1.^2)*cos(i)+2*D1*x1.*sin(i))./(2*pi*(x1.^2+D1.^2).^2);Hax2=-1*u*m2*((D2.^2-x2.^2)*cos(i)+2*D2*x1.*sin(i))./(2*pi*(x1.^2+D2.^2).^2);Hax=Hax1+Hax2;figure('color','g','NumberTitle','off','name','单个柱体Za');plot(X,Za1,'.-',X,Za2,'.-');xlabel('x(m)'),ylabel('y(m)'),title('单个柱体Za'); xlabel('x剖面方向(m)'),ylabel('y(nT)');figure('color','g','NumberTitle','off','name','两柱体Za');plot(X,Za,'.-r'),title('总Za'); xlabel('x剖面方向(m)'),ylabel('y(nT)');%延拓部分figure('color','g','NumberTitle','off','name','延拓后Za');n=512; % 傅立叶变换充零后的点数Hu=50; % 向上延拓高度50mZayt=up(Za,b,nx,n,Hu); %延拓function函数upplot(X,Zayt,'.-');xlabel('x剖面方向(m)'),ylabel('y(nT)'),title('上延50、80、100m 后Za');Hu=80; % 向上延拓高度80mZayt=up(Za,b,nx,n,Hu);%延拓function函数uphold onplot(X,Zayt,'.-m');Hu=100; % 向上延拓高度100mZayt=up(Za,b,nx,n,Hu);%延拓function函数upplot(X,Zayt,'.-r');legend('上延50m','上延80m','上延100m');stem(0,-100)stem(100,-20)ylabel('Y坐标/m');x=[0 100];y=[-100 20];s=['柱体'];text(x,y,s,'color','g')2、水平柱体三维上延（浅部异常在左侧）clc;clear;% 测点分布范围dx=5; % X方向测点间距dy=5; % Y方向测点间距nx=81; % X方向测点数ny=81; % Y方向测点数xmin=-200; % X方向起点ymin=-200; % Y方向起点x=xmin:dx:(xmin+(nx-1)*dx); % X方向范围y=ymin:dy:(ymin+(ny-1)*dy); % Y方向范围[X,Y]=meshgrid(x,y); % 转化为排列a=0;% 浅部柱体1参数I=90/180*pi; %倾斜角i=atan(tan(I)*csc(pi/2-a));%有效磁化倾角is R1=20; % 柱体半径 ms1=pi*R1^2;u=4*pi*10^(-7); %磁导率k=0.2;B=50000;H=B/u;M=k*H ; %磁化强度 A/mm1=M*s1; %磁矩D1=30; % 柱体埋深 mx1=X+100;x2=X;% 深部柱体2参数I=90/180*pi; %倾斜角i=atan(tan(I)*csc(pi/2-a));%有效磁化倾角is R2=50; % 柱体半径 ms2=pi*R2^2;u=4*pi*10^(-7); %磁导率k=0.3;B=50000;H=B/u;M2=k*H ; %磁化强度 A/mm2=M2*s2; %磁矩D2=100; % 柱体埋深 m% 两柱体Za总磁异常Za1=u*m1*((D1.^2-x1.^2)*sin(i)-2*D1*x1.*cos(i))./(2*pi*(x1.^2+D1.^2).^2);Za2=u*m2*((D2.^2-x2.^2)*sin(i)-2*D2*x2.*cos(i))./(2*pi*(x2.^2+D2.^2).^2);Za=Za1+Za2;% 两柱体Hax总磁异常Hax1=-1*u*m1*((D1.^2-x1.^2)*cos(i)+2*D1*x1.*sin(i))./(2*pi*(x1.^2+D1.^2).^2);Hax2=-1*u*m2*((D2.^2-x2.^2)*cos(i)+2*D2*x1.*sin(i))./(2*pi*(x1.^2+D2.^2).^2);Hax=Hax1+Hax2;subplot(221)surf(X,Y,Za),shading interp,xlabel('x(m)'),ylabel('y(m)'),title('两柱体总Za异常');colorbar;axis autosubplot(222)surf(X,Y,Hax),shading interp,xlabel('x(m)'),ylabel('y(m)'),title('两柱体总Hax异常');colorbar;axis auton=512; % 傅立叶变换充零后的点数Hu=80; % 向上延拓高度mZau=upcon(Za,dx,nx,ny,n,Hu); % Za向上延拓axis autosubplot(223)surf(X,Y,Zau),shading interp,xlabel('x(m)'),ylabel('y(m)'),title('上延80m两柱体Za 总磁异常');colorbar;axis autosubplot(224)surf(X,Y,Za2),shading interp,xlabel('x(m)'),ylabel('y(m)'),title('深部柱体Za磁异常');colorbar;axis auto三、输出图形及解释1、二维向上延拓（浅部异常在右侧）在0处100米深有一水平柱体，100处20米深有一柱体。