利用磁场水平调整方法实现航磁数据融合——以山东省齐河—禹城地区为例.pdf
《利用磁场水平调整方法实现航磁数据融合——以山东省齐河—禹城地区为例.pdf》由会员分享,可在线阅读,更多相关《利用磁场水平调整方法实现航磁数据融合——以山东省齐河—禹城地区为例.pdf(7页珍藏版)》请在咨信网上搜索。
1、 第 47 卷第 4 期物 探 与 化 探Vol.47,No.4 2023 年 8 月GEOPHYSICAL&GEOCHEMICAL EXPLORATION Aug.,2023doi:10.11720/wtyht.2023.1345吴成平,杨雪,于长春,等.利用磁场水平调整方法实现航磁数据融合 以山东省齐河禹城地区为例J.物探与化探,2023,47(4):1071-1077.http:/doi.org/10.11720/wtyht.2023.1345Wu C P,Yang X,Yu C C,et al.Aeromagnetic data merging based on magnetic da
2、ta leveling:A case study of the Qihe-Yucheng area,Shandong ProvinceJ.Geophysical and Geochemical Exploration,2023,47(4):1071-1077.http:/doi.org/10.11720/wtyht.2023.1345利用磁场水平调整方法实现航磁数据融合 以山东省齐河禹城地区为例吴成平1,杨雪1,于长春1,熊盛青1,范正国1,苏永军2,郝兴中3(1.中国自然资源航空物探遥感中心,北京 100083;2.中国地质调查局 天津地质调查中心,天津 300170;3.山东省地质调查院,
3、山东 济南 250104)摘 要:山东省齐河禹城地区于不同年份开展过高精度航磁测量,为了更好地实现老资料的合理再利用,对已有数据进行融合处理形成更大比例尺的高精度航磁资料是十分必要的。本文以 2011 年和 2017 年两次航磁测量数据为基础,经过数据改正、延拓、合并等处理,使用切割线调平、弱磁场区调平以及微调平等方法,融合形成了 1 2.5万高精度航磁数据,其疏密度为 25016.2 m,总精度为 1.3 nT。新数据具有更好的分辨力,包含更丰富的地球物理信息,填补了该地区 12.5 万航磁资料的空白,是后续数据处理和解释的重要基础资料。关键词:齐河禹城地区;数据融合;航空地球物理;磁场水平
4、调整中图分类号:P631 文献标识码:A 文章编号:1000-8918(2023)04-1071-07收稿日期:2022-07-13;修回日期:2022-08-19基金项目:ND工程区航空物探调查(DD20221640);青岛海洋科学与技术试点国家实验室“问海计划”项目(2021WHZZB0400);国家重点研发计划项目(2017YFC0602201);山东省地质勘查项目(鲁勘字(2020)63 号)第一作者:吴成平(1982-),男,正高级工程师,从事航空地球物理数据处理和综合解释工作。Email:chengpingwu 0 引言航空地球物理资料是深部地质矿产勘查十分有效的手段,而大比例航空
5、地球物理资料可以获取更加丰富的地质地球物理信息。然而,想要获取大比例的航空地球物理资料,往往需要付出更高成本。通过利用已有的航空地球物理资料,融合形成更大比例尺的资料是一种经济、有效和便捷的方法。山东省齐河禹城地区开展过多次航空地球物理测量,随着航空地球物理探测仪器和导航定位系统精度的不断提高,数据质量也不断提高,但是该地区缺少比例尺大于 1 2.5 万的高精度航空地球物理资料。2011 年和 2017 年分别在该地区开展了比例尺为 1 5万的航空地球物理测量,为融合形成 1 2.5万比例尺的航空地球物理资料提供了有利条件。在数据融合方面,前人主要利用数学方法或地球物理反演方法进行重磁数据融合
6、。支澳威等1以船磁为基准,在航磁、船磁水平位置和高程改正基础上,利用克里金插值方法实现两者的融合。刘福香等2采用空间域迭代法延拓、加权平均、回归分析等方法将不同观测手段获得的位场数据在同一标准下归算至同一平面或曲面,提出适用于多维多尺度重磁位场数据融合方法。郝燕玲等3基于移去恢复技术的融合方法构建高精度海洋重力图。Minty 等4利用加权最小二乘法使不同网格数据叠合区误差最小,达到数据融合、拼接的目的。而本文采用了不同的思路,在数据预处理基础上,使用磁场水平调整方法(简称调平)进行数据融合,以实测切割线磁场作为参考磁场,通过切割线调平、弱异常区调平方法5和微调平方法等进行齐河禹城地区航磁数据融
7、合,形成更大比例尺的航磁资料,更加真实地反映磁场强度大小及其变化趋势,具有更丰富的地球物理信息。物 探 与 化 探47 卷 1 数据准备齐河禹城地区航磁数据融合使用的航磁资料来源于 2011 年、2017 年在该地区完成的 1 5万航空地球物理测量工作6-7,测量方式均为沿地形起伏飞行测量。其中 2011 年数据(图 1a)基本覆盖研究区,但未覆盖研究区西北角部分地区,2017 年数据(图 1b)东侧边部不完整。上述两次航空地球物理测网交错布置,即在原15万比例尺数据(测线间距 500 m)基础上进行了加密,加密后测线间距为 250 m(图 2)。为了保证数据融合的效果,所使用的是未经过调平处
8、理的航磁数据,只进行了航磁正常场校正和磁日变校正a2011 年航磁 T 异常;b2017 年航磁 T 异常aaeromagnetic total filed acquired in 2011;baeromagnetic total filed acquired in 2017图 1 山东齐河禹城地区航磁异常(调平后)Fig.1 Aeromagnetic total field in Qihe-Yucheng area of Shandong Province(after leveling)图 2 山东齐河禹城地区航磁测网Fig.2 Aeromagnetic survey grid in Qih
9、e-Yucheng area of Shandong Province(图 3)。若使用调平处理的数据,因进行了滤波处理,在一定程度上融入了数据处理人员的主观认识,在局部区域会妨碍对原始磁场面貌的认识,增加后续数据融合的难度,导致更大误差。2 数据融合磁场水平调整是由于仪器观测、导航定位的点位、飞行高度以及日变观测等诸方面都存在误差,尤其在地形崎岖地区,不同测线与切割线之间飞行高度相差甚大,在交叉点形成很大的差别。在飞行后需利用测线与切割线交叉点测量观测的结果,按误差分配原理进行分配调整,使各测线统一到同一水平上8。切割线对磁场趋势及强度等起到有效地控制作用,本次数据融合采用包含切割线调平在内
10、的调平方法,能更加真实地反映磁场面貌,使用的其他调平方法包括弱磁场区调平方法(局部调平)以及微调平方法。2701 4 期昊成平等:利用磁场水平调整方法实现航磁数据融合 以山东省齐河禹城地区为例a2011 年航磁 T 异常;b2017 年航磁 T 异常aaeromagnetic total field in 2011;baeromagnetic total field in 2017图 3 正常场和日变校正后的航磁异常(未调平)Fig.3 Aeromagnetic field map after normal field reduction and diurnal correction(with
11、out leveling)2.1 数据合并由于不同年代获取的航磁数据飞行高度是有差异的,对飞行高度进行统计发现 2011 年数据的平均飞行高度为 152 m,2017 年数据的平均飞行高度为208 m,平均高度差为 56 m,本次以最新获取的数据为基准,对 2011 年数据进行向上延拓处理;同时,由于两者零值不统一,对 2011 年数据的零值进行了调整(调整量为 45 nT);最后完成两套数据的数据库合并,合并后的航磁面貌见图 4a。2.2 切割线调平切割线大致垂直于测线方向,用于联络测线磁场水平及航空地球物理测量质量检查。切割线调平是利用实测切割线网对测线的磁场水平进行调整,使研究区内所有测
12、线背景场的场值统一到同一水平或同一基点上9。在实际应用中,将测线与切割线交点处(交叉点)的磁场值作为调平参考值。切割线调平主要以最小二乘法原理为基础,使测线和切割线交点的差值平方和最小,实现磁场水平调整。针对齐河禹城地区 2011、2017 年两次获取的15万航磁资料,首先计算测线与切割线交叉点磁场差值,检查交叉点处磁场是否平静,同时去除磁场差值较大的交叉点,使其不参与调平过程;其次,对切割线磁场值进行修正9(也作统计调平),统一两次切割线磁场背景值 将一条切割线上各测点的磁场值减去此切割线与测线在平静交叉点处的磁场差值的算术平均值,并照此方法推广到研究区内所有切割线。最后,运用最小二乘法原理
13、,使所有交叉点上的差值平方和最小。切割线调平后,恢复了区域磁场的宏观面貌,消除了大部分磁场条带,对部分交叉点数据不足的测线或磁场复杂、短波长的磁场条带处理效果有限(图 4b)。2.3 弱磁场区调平方法弱磁场区调平方法5是一种局部调平方法,用于处理切割线调平处理后仍存在的沿测线方向的磁场局部条带,能使后续微调平等方法10-13达到更好的处理效果。其原理是根据测线及相邻测线磁场的空间变化趋势,在测线不同位置设置磁场水平调整量 在残存磁场条带两端调整量为零,中部根据磁场变化趋势设置调整量。设置好测线上的调整量后,依据点数或者距离将调整量分配至测线的每个数据点,完成测线的局部磁场水平调整(图 4c)。
14、可结合后续微调平结果,确定该阶段调整到何种程度以及重点调平区域。2.4 微调平方法经过上述调平处理后,航磁数据仍然存在一些微弱磁场条带,微调平方法利用频率域或空间域滤波方法去除此类磁场条带,分离地质信息与噪声10-13。空间域方法使用条件更为苛刻,而频率域方法应用更为广泛,其结果更加平滑,受局部异常或干扰异常的影响更小。空间域方法主要针对测线方3701物 探 与 化 探47 卷 a延拓和零值调整后的磁场异常;b切割线调平;c弱磁场区调平方法;d微调平amagnetic field after upward continuation and zero level altering;btie li
15、ne leveling;cleveling for the weak magnetic field;dmicroleveling图 4 磁场调平方法数据融合Fig.4 Data merging with leveling methods向与网格的 x 或 y 轴方向一致的情况,通过对测线方向进行空间域低通滤波,以及对切割线方向进行空间域高通滤波,获得调平误差,用原数据减去调平误差得到微调平处理结果。参数选择方面,低通滤波窗口应大于浅源局部异常宽度,高通滤波窗口选择 24 倍的测线间距。频率域方法则是对网格进行扩边和假值填充后,使用二维滤波进行处理,本次使用 Geosoft Oasis 软件实现
16、,其滤波方法为巴特沃斯高通滤波和方向余弦滤波组合,方向滤波器阶数取 2,巴特沃斯滤波器阶数取 6,波长选择 4 倍测线间距为宜。再通过 Naudy 滤波(波长通常为 5 倍测线间距)进一步从噪声中分离出地质信息,返回到调平结果中去。微调平基本消除了磁场条带,较好地保留了深部或浅部磁性体引起的磁场信息(图 4d)。3 数据融合效果3.1 数据质量评价对融合后的数据进行了质量评价,主要针对测网的实际分布情况以及数据精度,即测网疏密和磁4701 4 期昊成平等:利用磁场水平调整方法实现航磁数据融合 以山东省齐河禹城地区为例测总精度评价。测网疏密度反映了实际测线飞行质量,在研究区按一定间隔布置虚拟切割
17、线(通常按 35 倍测线距布置虚拟切割线),计算实际测线间距与设计的测线间距的差值,这些差值的平均值即为全区的测网疏密度。在研究区采用 1 km 间隔布置虚拟切割线,经计算测网疏密度为 25016.2 m,测线间距主要分布在 230270 m 之间(图 5)。航磁总精度采用测线与切割线交叉点残差的均方差进行评价,经计算融合形成的 1 2.5 万数据总精度为1.30 nT。3.2 数据融合前后对比分析为了便于对比分析数据融合前后的磁场细节特征,选取了研究区部分区域进行对比研究(图 6a)。经对比,融合后 1 2.5 万数据比原 1 5万数据突出了更多的细节信息。虽然高频异常多是浅表异常源引起的,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 利用 磁场 水平 调整 方法 实现 数据 融合 山东省 齐河 禹城 地区
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。