09重磁资料处理及解释新技术方法--徐世浙
重磁资料数据处理某些方法.
常的影响被削弱了年, ‘ 我们用上述方法编制了 , 币机算法语言程序 , 并作了模型试验 , , 后来又处理了实际资料闭有一定效果。
原始数据中存在有局部异常成分也是造成虚假异常的原因之一数据点随机分布叠加异常的划分趋势分析方法去处理规则格网数据吼 , 所以可以用不规则测网。
‘’ 至今仍是一个没有很好解决的问题、 , 还有许多工作要做二异常的反演 , 异常的反演一多模型最优化选择法、即根据观测异常求地质体的位置。
大小、产状和物性参数数据处理和解释工作中的一个不可缺少的环节欲反演的目标大体上有三种—是重磁资料求地质体用规则几何形体近似 , 的几何参数和物求一定范围物。
性参数性的空间分布的实质在于 , 。
求一定范围的物性分界面 , 反映某一地质层位的起伏目前采用的反演方法可分为两类直接法 , 。
根据观测整理的异常直接曲线拟合一选择法选择法 , 计算地质体的某些参数多用于解释简单的异常一次完成。
将实测异常曲线与一系列已知形状模型体产生的理论异常曲线进行比较当实测曲线与某一理论曲线符合为实际地质体的近似结果。
在给定的误差范围内选择法需迭代完成。
时 , 就将该理论曲线所对应的模型体作这里先介绍多模型最优化选择法在一个矿区、的应用情况 , 然后讨论此法的应用条件 , 年 , 武汉地质学院磁法组应用长方体组合模型采用了十五个模型 , 改进的马奎特法。
〕冀东对。
‘ 一区的地磁异常进行了反演川得到了各模型体的参数的理论曲线。
图表示反演得到的模型体的平面位置和由它算出理论曲线与实测曲线有些模型体的参数也。
一拟合得较好 , , 滋儡节静之火浓姆丫一之文侧气二么爹又一夕‘ 之毛‘ ‘’ 币—、、比较接近附近钻孔中的见矿情况月又卜女例如间见到体 , 孔位于第块之间 , 块和第米剩图一汤火在井深入米赤铁矿和 , 米磁性矿米磁块和第按计算结果是该处应有又如在第一、 , 性矿体块之间的孔 , 按反演推米磁性矿米到米米磁铁算在体 , 米左右应有实际钻探结果在当然此外 , 之间见到了矿区。
磁法数据处理过程
基点（日变站）的选择与基站T0的测定按照《地面高精度磁测技术规定》和《青海省当江地区地质矿产调查项目总体设计书》的要求，项目组于2008年7月15日和2009年07月29日开工前选择了日变站（分基点），日变站（分基点）位于平稳场内，且远离帐篷，高压线等干扰源，并在探头周围设置防护围栏，以防人畜干扰。
利用十子剖面法，在半径为2.5m的范围内,做了正南北,正东西两条剖面,经计算所得结果： 08年总场梯度变化为：1.8nT，符合规范及设计要求。
09年总场梯度变化为：北——南：ΔVi =1.4nT东——西：ΔVi=-1.5Nt总场梯度变化均小于总均方误差的1/2,符合设计及规范要求。
2、校正点的选择为了检验仪器每天工作是否正常，就必须建立校正点进行早晚基读数，仪器校正点选择在日变站或住地附近，并避开人文干扰，为了避免探头相互磁化造成影响，校正点离日变站点距大于20m设立，校正点打有木桩作标记。
测点观测采用总场测量方式，观测参数为地磁场总场强度，探头高度为1.70米（三节探杆）。
野外测量采用早校正点→观测点→晚校正点的闭合方式进行；每个测点均作两次读数，两次读数之差在2nT内方可储存；野外工作严格执行设计及规范要求，在外业生产前对磁力仪操作人员进行“去磁”，禁止磁力仪操作人员携带任何磁性物品；在磁力测量过程中，非操作人员远离测点。
对无GPS定位坐标的磁测数据或不符合设计规范要求的数据、干扰大的磁测数据经室内有关人员审核确定后，均予报废或次日补测处理。
每日野外测量工作结束后，对野外测量数据进行回放、日变改正和数据备份工作。
野外磁测资料整理的顺序是：第一、采用质子磁力仪实地采集野外数据；第二、将采集的数据回放到电脑中，根据现场记录和GPS行迹对数据进行编辑，前后对比检查实测数据是否合理，对于突变数据要求操作员进行重复观测，确保数据可靠，然后进行日变改正；第三、进行正常梯度改正、高度改正和总基改正，最后求出△T值。
带相位信息的一维大地电磁曲线对比反演法
（4） （5）
1 1 + dIog （ T ） （ T） / dIog T ｝ ! " 4｛ Pa
4
4 .1
模型实验
模型一
这是一个 4 层电阻率模型， 该模型是从多次被用于检验一维大地电磁反演的 COPROD 数据 中提炼出来的 （Meju et aI . ，1992） 。模型的电阻率和深度值如图 1b 中虚线所示， 实线是带相位 信息的曲线对比法反演的结果， 基本上正确反映了模型的电阻率值和层的位置， 反演模型的视 电阻率曲线和正演的视电阻率曲线拟合也比较好 （图 1a） 。迭代误差如图 1c 所示， 加入相位信 息的曲线对比法达到 25 次左右拟合， 其误差降至 1% 以下， 而不含相位信息的反演则需要 50 次 左右 （徐世浙， ， 可见收敛速度提高很多； 且迭代 100 次反演一条测深曲线仅需不到 1 min， 1995） 反演速度也得到提高。
Abstract
On the basis of the magnetoteiiuric curve comparison method inversion ，we introduce phase information into the inversion procedure . It means inserting phase correction parameter into apparent resistivity correction eiement during the inversion ，which makes the program more rapid and the procedure more robust . Comparing with other inversion methods of modei research，the resuit of using curve comparision method with phase information can represent the eiectricai property distribution of the modei correctiy，and can be used as the initiai modei of muiti - dimension inversion . K words Magnetoteiiuric，Curve comparison method，Phase information，Resistivity pseudo section
重磁资料处理及解释新技术方法
50000 950 45000 850 750 40000 650 550 35000 450 350 30000 250 150 25000 50 -50 20000 -150 -250 15000 -350 -450 10000 -550
y(m)
30000
25000
20000
15000
10000
0 0 5000 10000 15000 20000 25000 30000 35000 40000 45000 50000 55000
x(m) (a) Field continued from flying altitude to 3500m
55000
x(m) (b) Field re-continued from 3500m to flying altitude
55000
50000 2200 45000 2000 1800 40000 1600 1400 35000 1200 1000 30000 800 600 25000 400 20000 200 0 15000 -200 -400 10000 -600
50000 950 45000 850 750 40000 650 550 35000 450 350 30000 250 150 25000 50 -50 20000 -150 -250 15000 -350 -450 10000 -550
5150000
5150000
nT
5140000 220 200 5130000
升深4 宋深3 升深6 升深2-1 卫深4 升深2 升深101 卫深5 升深201 卫深501 升深3 芳深8
nT
5140000
180 5130000 160
实验四位场边缘识别程序设计实验
《重磁资料处理与解释》实验四位场边缘识别程序设计实验专业名称：地球物理学学生姓名：学生学号：指导老师：王万银、纪新林、纪晓琳、邱世灿提交日期：2016-1-31、基本原理地质目标体边缘时指断裂构造线、不同地质体的边界线，实际上是具有一定密度或磁性差异的地质体的边界线，在地质体的边缘附近，重、磁异常变化率较大，故所有的边缘识别方法均利用这一特点进行设计。
现在有重、磁位场边缘识别方法分为数理统计、数值计算以及其他三大类。
数值类边缘识别方法均利用极大值位置或零值位置确定地质体的边缘位置，其依据的理论基础是二度体铅垂台阶模型的重力异常特征。
在该模型边缘处重力异常总水平导数和解析信号振幅达到极大值、垂向导数达到零值。
故可以利用这些特征位置来确定二度体铅垂台阶的边缘位置，确定倾斜二度体、不规则二度体及三度体边缘位置的理论均为二度体铅垂台阶模型理论的推广，但确定的边缘位置和真实的位置有一定的偏差。
该偏差随着地质体边界形状、埋深、水平尺寸及物性差异等的变化而变化。
因此，边缘识别结果是一种定性或半定量解释结果，与定量解释结果有一定区别，识别结果可作为边缘位置定量反演的初值。
(1)垂向导数：垂向导数方法利用零值位置确定地质体的边缘位置，重力异常可以直接使用，对磁力异常必须转化为磁源重力异常或化极磁力异常才可以使VD R x y z：g(x, y, z)用VD Rx y去‘az ( 1.1)(2)解析信号振幅：解析信号振幅也是利用极大值位置来确定地质体的边缘位置适用于重、磁力异常ASM二THDR2 VDR2(3)总水平导数(THDR)2、输入/输出数据格式设计依据上述原理，现在对上述各种边缘识别方法进行程序设计。
(1-2)THDR(x, y,z)(1-3)2.1 输入数据格式设计2.3 参数文件数据格式设计将以上部分量保存在一个文件中，该文件名变量为 cmdfile ，字符串变量， 长度不超过 80，全路径名。
在该文件中保存的参数如下：输入数据文件名 input_file ，字符串变量，长度不超过 80；输出 vdr 数据文件名 output_file_vdr ，字符串变量，长度不超过 80； 输出 thdr 数据文件名 output_file_thdr ，字符串变量，长度不超过 80; 输出asm 数据文件名output_file_asm ，字符串变量，长度不超过 80factor_m ：扩边比例因子，实型变量（＞1）;本次实验给了正演的重力异常数据，为例如： .GRD 格式，均为实型变量DSAA201 -1000.000000 -1000.000000 5.549671E-01 5.549671E-01 5.897312E-012011000.000000 1000.000000 23.539846 5.634658E-01 5.987253E-015.721339E-01 5.808522E-016.078691E-01 6.171604E-012.2 输出数据格式设计 计算结果输出数据格式与输入格式对应，例如： 格式为.GRD 格式，均为实型变量DSAA201 -1000.000 -1000.000 -0.1465084 -3.6523044E-02 -3.2688729E-02 -3.2654848E-02-3.3138681E-02201 1000.000 1000.000 0.3190881 -3.3485338E-02 -3.2723978E-02-3.3061244E-02 -3.2748722E-02 -3.2787599E-02 -3.2927759E-023.总体设计此次程序采用IPO结构设计，首先通过读取cmd文件，得到相关输入参数：输入数据文件名gravity.grd、输出vdr文件名field_vrd.grd、输出thdr文件名field_thdr.grd、输出asm文件名field_asm.grd、扩边比例因子factor_m ;然后确定确定扩边网格的大小，扩边数据点号位置；再从观测面位场数据文件中读取数据。
重磁勘探课程教学实践与探索
重磁勘探课程教学实践与探索作者：肖锋许家姝来源：《黑龙江教育·高教研究与评估》 2019年第12期肖锋许家姝摘要：重磁勘探课程是地质类院校开设的应用地球物理系列课程的重要组成部分，是勘查技术与工程、地球物理学专业本科生的主干课程。
文章以重磁勘探相关课程为例，从助课制度、课程沿革、备课环节、教学方法、课程评价等方面对近年来的教学实践活动进行探讨，探索进一步完善重磁勘探课程建设的思路与途径。
关键词：重磁勘探；教学实践；地质类院校中图分类号：G642.0文献标识码：A文章编号：1002-4107（2019）12-0015-02重力与磁法勘探是两种历史悠久的地球物理工作方法，在地壳深部构造、区域地质调查、固体矿产勘查、石油天然气勘查等方面发挥了重要作用[1]。
重力勘探与磁法勘探作为专业课程一直是勘查技术与工程、地球物理学专业本科生的主干课程[2-3]。
仅以吉林大学勘查技术与工程专业为例，目前与重磁勘探有关的专业必修课程包括“重磁勘探原理与方法”（64学时）、“重磁勘探数据处理与解释”（32学时）。
选修课程包括“重磁勘探数据处理与解释实验”（32学时）、“重磁勘探新技术新方法”（16学时）等。
这说明重磁勘探系列课程在传统物探方法类课程中仍占有重要地位。
尽管地质类院校在开设重磁勘探相关课程方面已经具有丰富经验，但相应的课程建设与改革并没有举步不前。
有的优化教学内容，并在教学方法中将互动式、启发式、探究式多种方法相结合，为学生营造主动学习的氛围[4-5]；有的在“重磁勘探”中引入翻转课堂教学模式来提高教学效果[6]；有的以重磁异常为课程主线，以新方法和新知识为切入点，培养学生自主获取知识的能力[7]。
吉林大学重磁勘探课程的培养方案和教学大纲历经多次修订，逐渐趋于完善。
目前教学团队面临新老教师更替的挑战。
如何继承和发扬优良的教学传统，并将新颖的教学方法和手段引入课堂教学中，我们一直在思考这一问题，也在一步步探索和实践。
磁法数据处理_异常反演与解释的常用方法及常见问题探讨
磁法数据处理、异常反演与解释的常用方法及常见问题探讨张湖源（安徽省地质矿产勘查局313地质队）磁法勘探是最经典的物探方法，可广泛用于地质调查的各个阶段、工程地质及考古学等众多方面，尤其是在铁矿勘查中，更是必不可少的先行手段。
可以说没有其它的地球物理方法有如此广阔的应用范围，花费少而提供信息的丰富。
因磁参数多为矢量且常见干扰较多（与其它物探方法相比），使数据处理及异常推断解释变得较复杂，在实际应用中，产生不少使人困惑的问题。
结合笔者实际工作经验，对其进行初步探讨，供同行参考。
1磁法数据预处理常用方法对实测数据进行日变、基点、正常场等改正后，应注意消除异常数据的误差与干扰。
误差主要源于仪器的状态和操作及点位误差；干扰主要是指人文或地质因素的干扰。
在严格执行技术规范下，含有人文干扰的数据一般不作为成果；地质干扰通常指与勘探目标无关的地质因素引起的异常。
浅表局部的地质干扰体分为两类：①在空间上有一定分布规律的，如出露的岩石；②孤立的、无规律的，如滚石等。
孤立的地质干扰具有随机性，具有白噪声特征，而一些出露岩石虽然不具备随机特征，但往往具有相同的走向等特征，在某方向上具有有色噪声特点。
由于多数误差和干扰具有随机性特征，其均值为零，因此，可以通过小范围异常进行平均来消除这种误差和干扰。
平均圆滑能够有效地消除随机误差和干扰，但圆滑后有可能改变异常形态特征，给一些利用异常形态特征进行异常解释的工作带来困难，可采用多项式圆滑。
深部大型的地质干扰体引起异常特征为：磁性基底异常强度大，但相对平缓；岩浆岩（强磁性、较强磁性）有相当的强度，但有一定变化；火山岩或火山碎屑岩强度不大，但变化大。
对局部磁异常进行方差统计，方差较大被视为隐伏岩浆岩或火山岩异常，由此推断隐伏岩体的存在。
即可通过场的分离来剥离大型地质干扰体引起的区域异常。
剔处各种干扰后编绘图件时，成图数据位置最好使用xy坐标，以避免用点线号成图造成异常形态扭曲，尤其测线为斜线时更应注意。
2017重磁数据处理与解释复习题及答案
第九章重、磁异常的处理与转换1、什么是重磁异常的处理与转换？处理转换的目的及包括哪些内容？（1）定义：重、磁异常处理与转换的过程是根据重、磁异常的数学物理特征，对实测异常进行必要的数学加工处理，(2) 目的：提高信噪比，突出有用异常使实际异常满足或接近解释理论所要求的条件。
（3）内容：重、磁异常处理与转换的内容很广泛，它可以在空间域中进行，也可以在波数域中进行。
其主要内容包括：数据的网格化、异常的圆滑、异常的划分、异常的解析延拓、异常的分量转换、异常的导数计算、磁化方向的转换，异常的分离。
2、异常的预处理有哪些内容？数据的网格化、异常的圆滑、区域异常与局部异常的划分3、什么是数据的网格化，网格化的实质是什么？书中介绍了几种方法？（1）定义：将不规则的实测数据或数字化仪取出的数据换算成规则网格节点上的数据，这个过程就是数据的网格化。
（2）实质：数据网格化的实质问题就是对不规则的数据点进行插值。
（插值的方法很多，有拉格朗日多项式法、克里格法（Kriging）、最小二乘拟合法（多项式回归法）和加权平均法（近邻法）等等）（3）方法：拉格朗日多项式法，最小二乘拟合法4、什么是数据的圆滑，圆滑的实质是什么？书中介绍了几种方法？（1）定义：为了去掉数据中的误差或随机干扰对原始重、磁异常做圆滑处理。
（2）实质：其实质就是数学拟合。
（3）方法：1．徒手圆滑法；2．平均圆滑法；3．最小二乘圆滑法。
※异常的圆滑类型：剖面圆滑和平面圆滑5、趋势分析法的数学实质是什么？与最小二乘拟合、最小二乘圆滑法有什么区别？（1）实质：趋势分析法的实质是用一个多项式拟合区域场，是函数拟合。
（2）区别：趋势分析方法的实质与异常圆滑计算中的最小二乘圆滑是一样的，它们都属于函数拟合。
但趋势分析要一次性地利用全测区（或整条测线）中所有测点的异常数据，异常圆滑则是多次利用计算点附近一定范围内的数据。
也就是前者是整体拟合，后者是局部拟合。
正因为如此，圆滑计算时需要取计算点为坐标原点，计算点变化移动时坐标原点以及周围参与计算的已知点的异常值都随着变化，多次移动计算出多组多项式系数；而趋势分析计算时坐标原点必须固定且一次性求解出多项式的全部系数。
磁测资料的定性(半定量)解释推断.doc
学无止境物探学习参考资料磁测资料的定性（半定量）解释推断编写人：徐新忠目录前言高精度磁测是目前开展的战略性矿产远景调查项目中应用的主要物探方法之一。
磁测资料的解释推断是磁测野外工作的继续，更是应用磁测资料解决各种地质问题，直接或间接勘查矿产的必要步骤和关键环节。
为了帮助地质人员在应用综合方法进行综合研究开展地质找矿中利用好磁测资料，特编辑此学习资料。
本资料力求为地质人员提供简单、方便、实用、有效的磁测资料定性（半定量）解释方法，读者阅后对学习资料的批评、意见及建议请及时转告总工办。
一、磁测资料解释推断的基本任务和方法（一）磁测资料解释推断的基本任务磁测资料解释推断的基本任务，就是依据磁异常特征、岩矿石磁性资料和地质及其它物化探资料，正确判断引起磁异常的地质体的性质，并确定其空间位置和几何参数，同时结合地质规律，对地质构造和矿产分布做出相应的结论。
在区调和普查阶段，解释推断的首要任务是确定引起磁异常的磁性体的性质，其次是概略判断磁性体的空间位置和几何参数。
在详查和勘探阶段，则侧重于确定磁性体的空间位置和几何参数，以便为验证工程布置提供依据。
由于勘查对象的复杂性及勘查资料和解释推断工作者认识上的局限性，一个正确的完全符合客观实际的解释推断结果，往往需要多次经历解释、验证、再解释、再验证的过程才能完成。
（二）磁测资料解释推断的基本方法1、地质、物探资料对比方法将各种地质、物探资料综合起来，进行详细的对比和研究，按由已知到未知的原则，总结已知地质条件下的物探异常的特征和规律，然后利用这些特征和规律，结合解释地区的具体情况，对磁异常进行解释推断。
对比方法是区域调查和普、详查找矿工作中磁测资料解释的基本方法。
2、数学物理分析方法建立各种规则或不规则磁性体的物理模型或成矿模式，用数学物理方法求解模型周围空间磁场的分布及与模型各参量的关系。
根据这种分布和关系去分析待解释的异常，从而对引起异常的磁性地质体的赋存状态和磁化状态作出推断。
不同磁处理方法对铁皮石斛原球茎增殖生长的影响
不同磁处理方法对铁皮石斛原球茎增殖生长的影响徐忠传;曹秀;蔡国超;徐式近【摘要】研究了直接磁处理和磁处理水2种方法对铁皮石斛原球茎增殖生长的影响.结果表明,这2种磁处理方法在适宜条件下对铁皮石斛原球茎的增殖生长均具有促进作用,其中促进铁皮石斛原球茎增殖生长的最佳条件是:直接磁处理为25mT/80 min,磁处理水为0.8T.【期刊名称】《江苏农业科学》【年(卷),期】2013(041)005【总页数】3页(P229-231)【关键词】铁皮石斛;原球茎;直接磁处理;磁处理水;组织培养【作者】徐忠传;曹秀;蔡国超;徐式近【作者单位】常熟理工学院生物与食品工程学院,江苏常熟215500;常熟理工学院生物与食品工程学院,江苏常熟215500;常熟理工学院生物与食品工程学院,江苏常熟215500;苏州大学金螳螂建筑与城市环境学院,江苏苏州215123;常熟理工学院生物与食品工程学院,江苏常熟215500;苏州大学金螳螂建筑与城市环境学院,江苏苏州215123【正文语种】中文【中图分类】S567.239铁皮石斛(Dendrobium officinale Kimura et Migo)是传统名贵中药材，能养阴清热、益胃生津、解热镇痛，并有增强新陈代谢、抗衰老等作用[1]。
由于人们长期无节制地采挖铁皮石斛，同时也由于铁皮石斛的自然繁殖率很低[2]，使得其野生资源濒临枯竭，因此铁皮石斛的人工繁育显得十分重要，而铁皮石斛人工繁育的关键则在于其原球茎的增殖生长。
近年来，关于磁处理(包括直接磁处理和磁处理水)对植物生长影响的研究也成为热门话题。
直接磁处理能提高种子的发芽率，促进细胞生长。
范玲娟等研究发现，经直接磁处理的金鱼藻和水绵藻的光合强度均显著高于对照；经磁场处理后，小麦种子的发芽率和长势也有所提高[3]。
章继敏等在对原生动物(纤毛虫)生长速度的研究中发现，经磁场处理后，水密度会减小，水分子间的氢键遭到破坏，水分子的有序成分比例减小，改变了膜内外的渗透压平衡，有利于物质运输，促进细胞代谢[4]。
高精度重力资料的密度界面反演
高精度重力资料的密度界面反演肖鹏飞;陈生昌;孟令顺;杨金玉【摘要】Parker-Oldenburg法(PO)界面反演方法的核心算子是向下延拓,由于该算子的不稳定性,限制了PO反演法在重力密度界面反演应用中的发展.利用徐世浙提出的一种精度高、稳定性好的向下延拓的迭代法进行PO反演中的向下延拓,改进了PO反演中的稳定性.把改进的PO反演方法应用到理论模型和莺歌海基底反演中,均取得了较好的结果.【期刊名称】《物探与化探》【年(卷),期】2007(031)001【总页数】5页(P29-33)【关键词】重力资料;密度界面反演;Parker法【作者】肖鹏飞;陈生昌;孟令顺;杨金玉【作者单位】吉林大学,地球探测科学与技术学院,吉林,长春,130026;浙江大学,理学院,浙江,杭州,310027;吉林大学,地球探测科学与技术学院,吉林,长春,130026;浙江大学,理学院,浙江,杭州,310027;青岛海洋地质研究所,山东,青岛,266071【正文语种】中文【中图分类】P631由重力异常反演密度界面是位场处理解释中的一大问题。
其中地壳底界面即莫霍面及密度基底的反演就是典型应用。
很多学者提出了各种根据已知重力异常计算密度界面起伏的方法，包括利用等效棱柱体拟合密度界面进行反演。
该方法先计算每个棱柱体产生的异常，然后求和就能得到相应界面产生的重力异常，但该方法由于计算速度慢而没有广泛采用。
1973年，Parker在位场界面正演计算中引入快速傅里叶变换[1]，其方法具有速度快、功能强的特点，成为重磁异常定量解释中的工具，适合于含油气构造区的地球物理资料综合解释及联合反演。
其后Oldenburg在此基础上，提出了确定起伏密度界面的迭代反演法[2]。
由于该方法中包含了一个不稳定因子——向下延拓算子，导致了本方法的不稳定性，通常的做法在反演过程中加入一个低通滤波器，以限制高频振荡[3-4]。
但这种做法在消除高频假波的同时必然会对有效信号进行一定的压制，造成反演精度的降低。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
A
-550
5000
5000
0 0 5000 10000 15000 20000 25000 30000 35000 40000 45000 50000 55000
0 0 5000 10000 15000 20000 25000 30000 35000 40000 45000 50000 55000
x(m) (a) Field continued from flying altitude to 3500m
55000
x(m) (b) Field re-continued from 3500m to flying altitude
55000
50000 2200 45000 2000 1800 40000 1600 1400 35000 1200 1000 30000 800 600 25000 400 200 0 15000 -200 -400 10000 -600
5000
y(m)
延拓至4500m
45000
40000
35000
从4500m至 地形
5000
0 0 5000 10000 15000 20000 25000 30000 35000 40000 45000 50000 55000
0 0 5000 10000 15000 20000 25000 30000 35000 40000 45000 50000 55000
5050000 650000 660000 670000 680000 690000 700000 710000 720000
见见见见见
见见见见见
徐徐徐徐徐钓钓
徐徐徐徐徐徐徐徐 2km钓的的徐钓钓
2.位场曲化平的插值迭代法
航磁资料处理遇到的问题：曲化平；向下 延拓 现有的曲化平和向下延拓主要方法：等效 源法和傅立叶变换法。存在问题：延拓深度 <3-5点距；计算速度慢 我们开发的插值迭代法向下延拓的距离可 以达到20倍点距，计算深度：对1024*1024 数 据点的资料，平面向下延拓<60秒，曲面延拓 <600秒
35000
40000
45000
50000
55000
60000
再从1000m至200m
25000
资料点距50m
20000
15000
35000
40000
45000
50000
55000
60000
快速傅立叶变换法(FFT)的延拓效果 195m高度的航磁 从195m向下延拓至145m
从195m延拓至95m 资料点距50m
y(m)
0 0 5000 10000 15000 20000 25000 30000 35000 40000 45000 50000 55000
y(m)
0 0 5000 10000 15000 20000 25000 30000 35000 40000 45000 50000 55000
x(m) (c) Field continued from flying altitude to 4500m
50000 950 45000 850 750 40000 650 550 35000 450 350 30000 250 150 25000 50 -50 20000 -150 -250 15000 -350 -450 10000 -550
y(m)
30000
25000
20000
15000
10000
5150000
5150000
nT
5140000 220 200 5130000
升钓 4 宋钓 3 升钓 6升钓 2-1 卫钓 4 升钓 2 升钓 101 卫钓 5 升钓 201 卫钓 501 升钓 3
801 升钓 5
nT
5140000
180 5130000 160
a
5120000
D b
140 120 100 80 60
600 550 500 450
25000
400 350 300 250 200
20000
20000
150 100 50 0 -50
15000
15000
-100 -150 35000 40000 45000 50000 55000 60000 -200 -250 -300 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 50 0 -50 -100 -150 -200 -250 -300 -350 -400
5120000
B 芳钓
5110000
B 801 升钓5 芳钓
5110000
芳钓 8
升钓 4 升钓 6升钓 2-1 宋钓 3 卫钓 4 升钓 2 升钓 101 卫钓 5 升钓 201 卫钓 501 升钓 3
a
D b
芳钓 8
升钓 7
升钓 7
5100000
芳钓 3
芳钓 6 芳钓 7 徐钓 徐钓 6 1 芳钓 5 701 芳钓
50000
E
950 850 750 650 550
45000
45000
40000
40000
D350003Fra bibliotek000450
30000
30000
C
350 250 150
25000
25000
50 -50
20000
20000
B
-150 -250
地形上异常
15000
15000
-350 -450
10000
10000
准葛尔盆地陆东地区切割半径4km的局部异常 •见油、见气钻孔分布在局部异常的边缘； •显示火山岩的带状分布； •岩性与构造对油气藏的控制关系？
根据切割法成果推断的陆东地区火山岩分布图
实例3:徐家围子重磁异常解释
• 徐家围子的天然气富集受构造断裂与岩性 的共同控制，其中酸性火山岩（流纹岩）储层 产能很高，而中基性火山岩产能较低。如：徐 深1（53万方），升深2（31万方），肇深10 （24万方），汪深1（7万方）都在流纹岩中。 由于酸性流纹岩有一定的磁性，用磁法圈定流 纹岩的可能性是存在的。但是，在徐家围子地 区存在磁性较大的中基性火山岩，产生较大的 磁异常。流纹岩是酸性火山岩，磁性弱，产生 的异常低，常常被中基性火山岩的磁异常淹没， 在磁异常图上往往难以识别。因此，磁法解释 的主要任务是，如何从磁异常中，识别磁性较 低的流纹岩。
50000 950 45000 850 750 40000 350 300 35000 650 550 450 350 30000 250 150 25000 50 -50 20000 -150 -250 15000 -350 -450 10000 -550
y(m)
延拓至5500m
35000
200 30000 150 25000 100 50 20000 0 -50 15000 -100 -150 10000 -200
55000
x(m) (d) Field re-continued from 4500m to flying altitude
55000
50000 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 50 0 -50 -100 -150 -200 -250 -300
x(m) (e) Field continued from flying altitude to 5500m
55000
x(m) (f) Field re-continued from 5500m to flying altitude
55000
50000 500 45000 450 400 40000
迭代法平面向下延拓
迭代法延拓效果 200m
25000 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 50 0 -50 -100 -150 -200 -250 -300 -350 -400 -450
200m至1000m
重磁资料处理及解释新技术方法
浙江大学地球科学系 徐世浙 2007-9-27
1.重磁场分离的切割法； 2.重磁场曲化平的插值迭代法； 3.重磁场物性反演
1.重磁场分离的切割法
• 用切割法将浅源与深源 的磁异常,重力异常分离 开来
资料处理方法－ 资料处理方法－切割法
60 g(mGa) a 40 20 0 0 40 30 20 10 0 0 0 2000 4000 6000 II 8000 10000 0 10000 20000 30000 40000 I 10000 20000 Distance(m) b B c 30000 40000 A
5080000
肇钓 6 肇钓 8
-80 -100 -120 5060000 5070000
C
肇钓 4 肇钓 10
肇钓 6 肇钓 8
70 65 60 55 50 45 40 35 30 25 20 15 10 5 0 -5 -10 -15 -20 -25 -30 -35 -40 -45 -50
5050000 650000 660000 670000 680000 690000 700000 710000 720000
Flying altitude
Magnetic field at flying altitude
55000
y(m)
y(m)
曲面向下延拓 资料点距100m 地形
55000
50000 5400 5300 5200 5100 5000 4900 4800 4700 4600 4500 4400 4300 4200 4100 4000 3900 3800 3700 3600 3500