地磁日变改正中磁湾磁钩等磁扰对海洋磁力资料的影响.pdf
《地磁日变改正中磁湾磁钩等磁扰对海洋磁力资料的影响.pdf》由会员分享,可在线阅读,更多相关《地磁日变改正中磁湾磁钩等磁扰对海洋磁力资料的影响.pdf(8页珍藏版)》请在咨信网上搜索。
1、 12 海洋科学/2023 年/第 47 卷/第 6 期 地磁日变改正中磁湾磁钩等磁扰对海洋磁力资料的影响 任相宇1,2,3,付永涛1,2,4,周章国1,2(1.中国科学院海洋研究所,山东 青岛 266071;2.中国科学院海洋地质与环境重点实验室,山东 青岛 266071;3.中国科学院大学,北京 100049;4.海洋油气勘探国家工程研究中心,北京 100028)摘要:在海洋磁力测量过程中,磁扰形态变化复杂的特性致使地磁日变的改正仍存在着许多困难。本文根据黄海、东海和南海北部十余年采集的超过 104 km 的海洋和陆地地磁资料,利用对比分析的方法,总结了 3 种磁扰现象对海洋磁力资料的影响
2、规律:1)微扰的海陆日变值持续时间短,幅度比值接近11,对海洋磁力资料影响很小,利用低通滤波可以完全消除其影响;2)磁钩海陆日变幅值比为 12,持续时间在 1 h 以内,校正后的地磁异常值仍有较大的误差,利用最小曲率法可以基本消除其影响;3)磁湾的持续时间长,海陆日变幅值变化大,无法准确消除其影响,采用最小曲率法等插值方法处理后,仍残留假异常,需要额外关注。关键词:磁湾;磁钩;日变改正;海洋磁力测量 中图分类号:P631.2+24 文献标识码:A 文章编号:1000-3096(2023)6-0012-08 DOI:10.11759/hykx20210428002 地磁日变改正是海洋磁力测量过程
3、中最基础的工作之一,但同时也是海洋磁力成果资料中误差的主要来源,若不进行准确校正,会导致虚假的磁异常,造成解释上的错误1-5。地磁日变主要与太阳辐射对地球电离层的影响有关,具有极强的时效性,每时、每日的观测结果都不相同。因此,如何提高日变改正的精度,一直是海洋磁力测量中需要优先考虑的问题。从 20 世纪60 年代开始,国内外众多学者针对海陆日变幅值、时差、基值等问题6-11,以及深远海日变测量技术12-17和地磁日变的校正方法18-23开展了大量工作,保障了近远海磁力测量工作,测网精度达到技术规范的标准。磁扰是地磁日变改正中典型的难点,对于较大的脉动等磁扰的校正,通常采用边刚3提出的校正方法:
4、将磁扰日期间的静日变化和扰动变化分离,分别对其进行改正。但对于相对较弱的磁钩(钩扰,magnetic crochet)、磁湾(湾扰,magnetic bay)等磁扰还缺乏充分的讨论。在实际工作当中,对于强度较小、持续时间较短的磁扰,通常认为是随机噪声的影响,利用平滑等方法来消除其影响,但仍可能会残留扰动,降低地磁日变改正的精度。大量实测资料经分析显示,磁湾、磁钩等磁扰现象广泛存在,校正后的地磁异常值出现了显著的假磁异常。本文以微扰、磁湾和磁钩等磁扰现象为研究对象,对比分析海洋磁力资料和陆地地磁资料,评价微扰、磁钩和磁湾等磁扰现象对海洋磁力成果资料的影响。1 数据来源数据来源 数据来源于中国科学
5、院海洋研究所 20052018 年的南黄海、东海和南海北部采集的海洋磁力资料和陆地地磁日变资料。海洋磁力资料大多数是在油气勘探任务中采集,作业船只以约 9 km/h 航速匀速航行,使用的仪器是Marine Magnetics 生产的海洋磁力仪,仪器精度为0.2 nT24-25。地磁日变观测在工区同纬度陆地区域进行,使用的仪器有 Sentinel 以及 WCZ-3S 等多种磁力仪,其中 WCZ-3S 在平静日测量时的精度约为0.3 nT,Sentinel 基站磁力仪的精度为0.1 nT,满足高精度地磁日变观测的需求26-27。2 数据分析数据分析 在海洋磁力测量期间,严格意义上的地磁平静日 收稿
6、日期:2021-04-28;修回日期:2022-06-26 基金项目:中国科学院海洋地质与环境重点实验室开放基金项目(MGE2022KG11)Foundation:Open Fund of Key Laboratory of Marine Geology and Envi-ronment,Chinese Academy of Sciences,No.MGE2022KG11 作者简介:任相宇(1996),男,山东潍坊人,硕士研究生,地质工程专业,现从事海洋重磁测量研究工作,E-mail:;付永涛(1970),通信作者,男,湖南绥宁人,副研究员,从事海洋重磁测量研究工作,E-mail: Marin
7、e Sciences/Vol.47,No.6/2023 13 天数很少,一个月内经常少于 3 d。而磁扰日相对较多,这些磁扰有梯度相对平缓、幅值在 10 nT 以内的微扰,也有幅值较大的磁钩、磁湾等磁扰,磁暴则相对较少。2.1 陆地地磁日变数据处理陆地地磁日变数据处理 陆地地磁日变数据处理主要包括五点平滑和基值处理两部分。地磁基值理论上为台站所在位置的基本地磁场值。本文中确定基值的方法是平均值法:选择航次期间某几天磁静日 24 h 地磁日变观测数据的平均值作为该地磁日变站的地磁场总强度值,以此为地磁日变基值28。用五点平滑后的地磁日变观测数据减去基值即可得到地磁日变值。2.2 微扰型磁扰特征分
8、析微扰型磁扰特征分析 微扰在本文中特指幅值在 10 nT 以内、持续时间在几分钟或十几分钟的磁扰以及突然脉动(一种突然的脉冲现象;水平分量和总场值突增或突降,幅度为 1 100 nT;时间持续 110 min)。微扰的海陆日变幅值很小,幅度比值接近 1,可以利用低通滤波方法消除其影响。琼东南盆地的测线13QDN01微扰型磁扰特征明显(图 1)。从图 1 中可以看到,测线 13QDN01 在世界时 16 时、17 时和 19 时左右的陆地地磁日变值与海测值减去正常场值同时出现有幅值约 5 nT、持续时间约 15 min 的磁扰。两者之间的幅值基本相等,因此校正后的地磁异常值(黑色)在对应时间段内
9、表现为相对平缓,只有一些幅值很低的残留扰动。本文共统计了 25 条存在微扰型磁扰的测线,这些微扰型地磁扰动对地磁日变改正的影响很小,经校正后的地磁异常值幅值不超过 1 nT。图 1 微扰型磁扰示意图 Fig.1 Perturbation schematic 注:为方便进行曲线比对,将陆地地磁日变值加上 30 nT,海测值-正常场值减去 10 nT 2.3 磁钩型磁扰特征分析磁钩型磁扰特征分析 钩扰是一种持续时间一般在 1 h 以内,幅度变化超过 10 nT 的磁扰,起始较急,形状类似于一个钩子,常常伴随有短波消失4。即便测量船只以 1.8 km/h航速航行,也往往会跨越 1020 km,产生显
10、著的峰值或谷值(幅度超过 10 nT)(图 2)。钩扰与太阳耀斑同时发生,是一种由太阳辐射引起的地磁效应,故又称太阳耀斑效应。在太阳耀斑发生时,紫外线辐射会突然增强,引起电离层的电离浓度突然增高,从而在电离层中形成一个短暂的电流体系,这就是钩扰的成因29。因此,磁钩只限于日照半球。图 2 磁钩型磁扰及低通滤波和最小曲率法处理前后的地磁异常值 Fig.2 Magnetic crochet schematic diagram and geomagnetic anomalies before and after low-pass filtering and mi-nimum curvature 注:
11、为方便进行曲线比对,将地磁异常值减去 10 nT,低通滤波处理后地磁异常值减去 20 nT,最小曲率法处理后地磁异常值减去 30 nT 以琼东南盆地测线 18QDN01 为例,在图 2 中世界时约 3 时出现了第一次钩扰,陆地地磁日变值有大约 10 nT 的幅值变化,海测值与正常场值的差有大约 20 nT 的幅值变化。世界时 4 时至 5 时之间出现了第二次钩扰。陆地地磁日变值的变化幅度约 20 nT,海测值减去正常场值的变化幅度约 35 nT。两者变化幅度的差异导致了地磁异常值在第二次钩扰时间段内出现波动,幅度约 15 nT。三条曲线在钩扰时间段内同时出现了变化形态相同但幅值不同的波动。这也
12、说明了传统的地磁日变校正未能从海洋磁测原始数据中完全消除钩扰的影响,校正后的地磁异常值仍存在高达 15 nT 的误差。对于这种海陆日变幅值不同的现象,有两种可靠的解释:一是由于受太阳辐射的影响,电离层的电离浓度突然增高,在海水和陆地上激发感生电动势和电流,而海洋电导率比陆地电导率高,这就导致了钩扰在海洋日变数据中的幅值要比陆地日变数 14 海洋科学/2023 年/第 47 卷/第 6 期 据大;二是由于地质体引起的真实异常。本文统计了36 条存在磁钩的测线,结果显示:这些测线的海陆地磁日变数据在钩扰时间段内都出现了幅值比在12 的波动。大量数据证明这种现象在海洋磁力资料中广泛存在,可以排除地质
13、体引起的异常。钩扰对海洋磁力资料造成的误差虽然可以被大洋中高梯度值(超过 100 nT)掩盖,但对于重磁震联合反演、测线网调差等工作的影响是显著的。在重磁震联合反演中,一般认为这种磁钩引起的假异常由地质体或地层磁性差异引起,这将导致错误的反演结果。图 2 显示了经过低通滤波和最小曲率法处理后的地磁异常值。由于钩扰的幅值相对较高,校正后的地磁异常值仍有较大的误差,利用低通滤波也不能消除。因此,本文利用最小曲率法来消除其影响。从图 2 中可以看到,低通滤波处理后的数据虽然比原始数据更加平滑,但钩扰处的地磁异常值仍存在假异常。而经过最小曲率法处理后的地磁异常值在钩扰处的幅值明显降低,基本消除了钩扰对
14、海洋磁力资料的影响。图 3 是 2018 年琼东南盆地两条相交测线(18QDN02 和 18QDN03)的交点情况,从图 3(a)中可以看到有磁钩假异常的测线段与无磁钩假异常的测线段的交点差在14.2 nT,这会显著提高测线的平均交点差值,破坏测线与所有相交测线之间的线性特征,从而扰乱该测线的平差参数,降低海洋磁力成果资料的质量。图 3(b)是经过最小曲率法处理后的交点差,经过最小曲率法处理后的两条测线的地磁异常值几乎相等,交点差明显降低,基本消除了钩扰对测线交点差的影响。表 1 是测线 18QDN02与相交测线的交点差情况,最小曲率法处理后与测线 18QDN03 交点差降低为 0.08 nT
15、,测线准确度由3.55 nT 提升为0.85 nT,充分证明了该方法的可行性。图 3 最小曲率法处理前后的测线交点差 Fig.3 Crossing error before and after minimum curvature 表 1 测线 18QDN02 交点差统计 Tab.1 Crossing error statistics of line 18QDN02 相交测线名称 原始交点差 改正值 改正后交点差08QDN01 0.34 0.34 07QDN01 0.05 0.05 06QDN01 1.43 1.43 18QDN03 14.20 14.28 0.08 11QDN01 0.42 0
16、.42 06QDN02 1.94 1.94 08QDN02 1.97 1.97 11QDN02 1.75 1.75 11QDN03 0.43 0.43 2.4 磁湾型磁扰特征分析磁湾型磁扰特征分析 磁湾是亚暴的磁场表现,持续时间约半小时到几小时,幅值可超过 10 nT,随时间变化形态像海湾,有正湾和负湾之分。在高纬度地磁台,子夜前常记录正湾,子夜后为负湾,表明分别有东向电集流和西向电集流流过头顶5,30。以琼东南盆地 2005 年的测线 05QDN01 为例,在图 4 中世界时 12 时至 14 时内出现了一次持续时间在 1.5 h 左右的湾扰。湾扰形态变化复杂,幅值变化大,出现多个波峰或波谷
17、。陆地地磁日变值、海测值与正常场值的差、地磁异常值在湾扰时间段内的 Marine Sciences/Vol.47,No.6/2023 15 波动的频率和相位虽然一致,但振幅却有很大的差异。若不对其进行仔细处理,会导致时间较长的假磁异常。与磁钩相同,磁湾对海洋磁力资料的影响也是由于海陆电导率不同而造成的,但磁湾持续时间长,海陆日变幅值变化大,可能会与地质磁异常混淆,因此在判断磁湾时需要仔细检查相邻测线磁异常、交叉测线在交点处附近的磁异常。图 4 磁湾型磁扰及低通滤波和最小曲率法处理前后的地磁异常值 Fig.4 Magnetic bay schematic and geomagnetic anom
18、alies before and after low-pass filtering and minimum curvature 注:为方便进行曲线比对,将低通滤波处理后地磁异常值减去10 nT,最小曲率法处理后地磁异常值减去 20 nT 图 4 显示了经过低通滤波和最小曲率法处理后的地磁异常值。低通滤波削弱了湾扰时间段内幅值较低的磁异常,但仍残留幅值较大的磁异常。经过最小曲率法处理后的地磁异常值在湾扰处的幅值虽然降低了 10 nT左右,但由于时间长达 1.5 h,在测线上跨度长达 15 km,会对真实的地磁异常进行插值,导致最小曲率法处理后的地磁异常并不准确。为了验证最小曲率法处理对磁湾型磁扰
19、的有效性,本文以 2005 年的测线 05QDN01 及两条相交测线为例。从图 5(a)中可以看到,有磁湾假异常的05QDN01 测线与无磁湾假异常的测线 05QDN06 和测线 11QDN04 的交点差均在10 nT 左右。图 5(b)为经过最小曲率法处理后,测线 05QDN01 与测线05QDN06 的交点差为0.12 nT。但由于磁湾的幅值变化大,测线 05QDN01 与测线 11QDN04 的交点差仍有 3 nT 左右的误差。本文统计了 15 条存在磁湾的测线,结果表明:经最小曲率法处理后的大部分交点差仍残留较大的误差。因此,最小曲率法并不能彻底消除磁湾对海洋磁力资料的影响。表 2 是
20、测线 05QDN01 与相交测线交点差的交点情况,最小曲率法处理后的交点差分别降低为0.12 nT 和 3.35 nT,测线准确度由2.93 nT 提升为1.1 nT,得到显著改善。图 5 最小曲率法处理前后的测线交点差 Fig.5 Crossing error before and after minimum curvature 相较于磁钩而言,磁湾对测线之间的线性误差、平差参数以及校正后的地磁异常值都会产生更严重的影响,给磁力资料处理解释带来困惑。因此,需仔细分析磁湾,尽可能降低其对海洋磁力资料的影响。3 讨论讨论 如前所述,地磁日变改正是海洋磁力测量中最主要的误差来源之一,而关于磁扰的处
21、理则是地磁日变改正当中最复杂的一部分。因此,国内外学者 16 海洋科学/2023 年/第 47 卷/第 6 期 表 2 测线 05QDN01 交点差统计(单位:nT)Tab.2 Crossing error statistics of line 05QDN01(unit:nT)相交测线 原始交点差 改正值 改正后交点差05QDN02 0.15 0.15 05QDN03 1.45 1.45 05QDN04 0.48 0.48 05QDN05 0.01 0.01 11QDN04 6.07 9.42 3.35 05QDN06 10.38 10.26 0.12 20QDN01 1.39 1.39 11
22、QDN05 2.4 2.4 05QDN07 0.59 0.59 开始寻找磁扰日的改正方法,吴水根等31对比当地时法、小时平均法、北京时法共三种方法对典型磁暴日、全区磁力交点均方根误差计算,结果表明按照北京时法处理得到的均方根误差最小,其实质是对原始数据进行世界时校正。而磁暴、磁扰往往是全球同时发生,在这种情况下磁暴、磁扰的变幅远远超过了磁静日变幅,因此取得了较好的效果。李兴康等32依据大洋实测地磁资料,注意到 530 km 大洋岛屿台站与工区锚系潜标日变站的地磁日变值在磁暴期间的相位和幅值是吻合的,磁暴期间的一条磁力测线,经远端洋岛台站地磁日变资料校正后,与 2 条平静期测线的交点差分别为3.
23、20 nT 和0.66 nT,因此需要适当关注磁暴期间的海洋地磁资料的处理。刘天佑等33提出应用维纳滤波和三次样条函数法对磁扰和磁暴日变曲线进行改正。理论模型试验表明,该方法优于目前使用的地方时法、小时平均法和北京时法,可取得较好的效果。边刚等3提出利用扰日法对磁扰日的地磁日变数据进行校正。该方法充分利用磁扰前后的磁静日变观测资料,实现了磁扰日的静日变化与扰动变化分离,结合各自的时间特点,分别对其进行改正,经实例验证表明,该方法可以提高磁扰日地磁日变改正精度,并且其日变改正精度与时间差(经度差)不存在相关性。傅里叶分析在磁场建模中具有广泛应用,可以运用傅里叶级数展开的方法把一天的太阳静日变化曲
24、线表示成各次谐波的叠加。张连伟等21基于傅里叶级数建立了日变数据处理谐波分析模型;应用于大批量实测日变数据处理,确定了谐波截断次数为5,实现了日变基值、平静日变改正和磁扰改正的合理分离,有助于磁扰和磁暴期间的海洋磁力测量数据处理。刘帆等20采用最小二乘法确定各次谐波系数建立太阳静日变化模型,以此进行地磁观测数据中太阳静日变化和磁扰的分离,并应用于海洋磁力测量的日变数据处理。海洋调查规范 第 8 部分:海洋地质与地球物理调查28中规定:必须准确记录初动、持续、消失的时间,并及时通知测量船,磁暴或磁扰持续时间较长时,应进行补测。但在实际工作中,对于磁湾、磁钩等持续时间较短的磁扰,无法做到补测。综上
25、,国内外学者关于磁扰日的改正主要着重于如何分离静日变化和扰动变化,结合各自的时间特点,然后分别对其进行校正。此类方法虽然方便快速,但无法彻底分离所有磁扰,也无法确定其改正的精度。而本文通过分析大量实测资料,总结了微扰、磁钩和磁湾型磁扰的海陆日变幅值变化规律,以此为基础,通过观察各条测线人工分离了这些磁扰。但在海洋磁力数据处理中,数据量通常较大,为了减轻磁扰校正工作导致的负担,需要提出合理的算法去分离微扰、磁湾和磁钩等磁扰现象,然后分别对其校正。此外,本文根据微扰和磁钩型磁扰对海洋磁力成果资料的影响规律提出了利用低通滤波和最小曲率方法。结果显示:经低通滤波和最小曲率法处理后的测线交点差显著降低,
- 配套讲稿:
如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。