严钦渝 YAN Qin-yu
（成都理工大学地球物理学院，成都 610059）
摘要： 本文是用解析延拓的方法对某地区铜铁矿磁异常数据处理，并做了相应的解释。利用场的空间换算就可以划分不同深度的地质体的异常。也可以在目标场源空间深度最大限度地突出目标引起的磁异常, 更好地确定其几何形态。
教育期刊网 http://www.jyqkw.com
关键词 ： 延拓；数据处理；铜铁矿
中图分类号：P631 文献标识码：A 文章编号：1006-4311（2015）03-0323-02
作者简介：严钦渝（1987-），男，四川南充人，硕士，研究方向为油气与矿产地球物理。
1 区域地质概况
区域出露地层简单，以晚古生代海相碎屑岩，碳酸盐岩，第三纪陆相火山岩为主。岩浆较发育，主要为喜山期石英单斜辉石苏长岩、石英二长闪长岩、黑云花岗闪长岩、黑云二长花岗岩等侵入体及大量中酸、中基性岩脉。区域构造主要表现为一系列东西向构造挤压带，其次是北东、北西、南北向小规模次级断裂。
2 地球物理特征
1?誜100万航磁ΔT平面等值线图显示，异常发育，分布较复杂，以正异常为主，正负相间产出。测区159异常呈椭圆状元宝形，两端略北翘，形态规整，呈北西西向分布，75nT等值线所圈定的范围为：长约45km，宽约15km，南翼等值线较北翼稍稀，最高为159nT。在正异常北偏东有明显的负异常，峰值为105nT，正负异常中心相距约30km。在测区南东还分别有一个正负异常中心。
岩石地球物理性质的差异是地球物理勘探的前提，充分了解区内各种岩（矿）石的物性参数，是解释物探异常场必不可少的条件，也是对测量数据进行深入研究的基础。但是由于测区范围大、海拔高、地形复杂，通信和交通极为不便，地面物性测量工作程度低。为此，笔者收集了前人工作的资料，结合项目组野外实际采集的标本的测量结果，现将该区各种岩（矿）石的磁性参数叙述如表1。
3 数据处理方法及过程
3．1 圆滑滤波
在野外实测的磁异常，总是包含有各种观测误差和浅部磁性不均匀引起的异常，这些异常叠加在有用的异常之上，使异常发生不规律的跳动以至难以利用和解释。为了能利用有用异常，就需要对实测异常进行圆滑。圆滑方法的好坏直接关系到后面对磁场数据进行延拓的效果，一般主要根据算区异常特征选用合适的圆滑公式。
3．2 积分插值延拓
重磁位场的空间换算实质上就是由一个已知观测平面上的位场来计算不包含场源的不同高度的位场。
计算可得上延公式，即：
是用所有的圆环半径作为同心圆量板，然后将欲要计算的点与量板圆心重合，求得与该计算点有关的各个圆环上场的平均值
，再乘上相应的各环的延拓系数，最后累加求和便可求得某高度平面上的延拓值。
计算可得下延计算公式，即：
4 成果解释与推断
为了研究本区的地质构造和磁异常与磁性岩体的关系，笔者对地面磁测资料进行了转换处理，做了上延400、300、200、100及下延80、60、40、20米平面等值线图。
磁场随距离的衰减速度是与磁性体的形态体积有关系的，即体积大，磁场衰减慢，体积小磁场衰减快。由上延平面等值线图1、图2、图3和图4对比可看出上延高度越大，一些由深部岩体和地表花岗闪长岩体或中基性岩体引起的异常明显减少了很多。由图1可看出主要磁异常就是深部大的磁性体和大量磁铁矿露头产生的磁场。要是上延高度再大一点有可能这种效果会更为明显。
从下延图5看出一些重要的异常浅部不太明显了，可能是由于下延深度太大。又出现了许多新的异常，很可能是地下深部的磁性岩体引起的，图7效果最为理想，清楚的反映出了一些低缓异常，这对于区分是磁铁矿引起的异常还是其它干扰引起的异常很有帮助。
对上面各个延拓等值线平面图对比分析，根据磁异常的走向和磁异常梯度带变化特征，划分了两组构造破碎带，一组为近东西向的构造，共6条；另一组为近南北向的构造，共6条；圈定具有代表性的磁性岩体7处。见图8。
图9为该地区磁性基底等深线图，该等深线图反映了全区磁性基底埋深的变化特征，也可以用来定性分析隐伏磁性体的分布。从基底等深线图可以看出，在仲青、沙松南和桌布日阿岩体单元，磁性基底埋深较浅，而对应于该区在上延400米等值线平面图上为正磁场区，说明该区存在隐伏的磁性体分布。从木质顶到阿木弄一线，磁性基底埋深相对较深，而地面磁测却存在局部的高磁异常，而对应于该区在下延40米等值线平面图上为低缓磁场区，说明引起该区高磁异常很可能是由近地表的磁铁矿产生的，而磁铁矿的下延深度可能并不深。
5 结论与建议
通过本次高精度磁测工作，我们可以得出如下结论：
①推断解释了隐伏地质构造和磁性块体的分布，分析了全区磁性基底埋深的变化特征。特别是测区西南日阿岩体引起的磁异常，异常范围较大，存在隐伏的磁性体分布的可能性大，是最具有意义的成矿远景区。这些对今后的工作具有指导性的意义。
②根据对本次磁测工作的质量评价表明，本次磁测内外业工作质量可靠，所有资料符合规范要求，解析延拓方法在该地区数据处理解释中效果明显。
教育期刊网 http://www.jyqkw.com
参考文献：
[1]胡金柱.FORTRAN77与应用软件开发方法[M].化学工业出版社，1991.
[2]邓巍巍，王越男．Visual FORTRAN编程指南[M]．人民邮电出版社，2000．
[3]马幼捷，问虎龙，周雪松，李季，杨海珊.采用延拓法的风电系统稳定模型二维参数分岔边界的计算与研究[J].中国电机工程学报，2010（19）.