CN111045090A - 一种磁异常网格快速缝合方法 - Google Patents

一种磁异常网格快速缝合方法 Download PDF

Info

Publication number
CN111045090A
CN111045090A CN201911407381.4A CN201911407381A CN111045090A CN 111045090 A CN111045090 A CN 111045090A CN 201911407381 A CN201911407381 A CN 201911407381A CN 111045090 A CN111045090 A CN 111045090A
Authority
CN
China
Prior art keywords
grid
stitching
gridb
grida
correction
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201911407381.4A
Other languages
English (en)
Other versions
CN111045090B (zh
Inventor
喻翔
张东辉
武鼎
杨龙泉
胡小婕
汪硕
梁永顺
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Research Institute of Uranium Geology
Original Assignee
Beijing Research Institute of Uranium Geology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Research Institute of Uranium Geology filed Critical Beijing Research Institute of Uranium Geology
Priority to CN201911407381.4A priority Critical patent/CN111045090B/zh
Publication of CN111045090A publication Critical patent/CN111045090A/zh
Application granted granted Critical
Publication of CN111045090B publication Critical patent/CN111045090B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
    • G01V3/081—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices the magnetic field is produced by the objects or geological structures
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Electromagnetism (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明属于铀矿地球物理勘查技术领域，具体涉及一种磁异常网格快速缝合方法。本发明包括以下步骤：步骤1、在铀资源勘查区域内，划分A、B两区分别开展地面高精度磁测，在两区开展不同比例尺，但A、B区有部分的磁数据重合，分别形成网格图；步骤2、对步骤1中的进行化极处理；步骤3、对步骤2中两张网格图进行趋势面去除处理；步骤4、对步骤3中趋势面去除后的网格开始进行缝合；步骤5、选取网格图校正的权重；步骤6、拓展网格图的校正宽度，得到校正后网格图；步骤7、将步骤6中校正后的两张网格图进行插值，重复步骤3‑步骤7直至形成最终缝合后的网格图。本发明能够准确、迅速地拼接整合同一测区的磁异常数据。

Description

一种磁异常网格快速缝合方法
技术领域
本发明属于铀矿地球物理勘查技术领域，具体涉及一种磁异常网格快速缝合方法。
背景技术
在常规铀资源地球物理勘查中，航空磁法、地面高精度磁法测量是铀矿勘查较为有效和直接的探测方法。同一勘查区，由于面积较大，往往需要长时间(数年)才能完成全区的面积性测量任务，或者同一勘查区由多支地球物理测量队伍共同完成，磁异常网格可能为不同时期、不同测量间距、不同投影下产生，磁异常网格也可能由不同插值算法计算形成，为了及时掌握勘查区全区的磁异常特征，快速的整合区内磁测数据则显得尤为重要，对后续开展其他地球物理方法及钻探工程具有重要的指导意义。
在铀资源地球物理勘查评价工作中，以往并没有明确的方法或没有明确的针对磁异常网格快速缝合的方法流程，直接影响后续地质地球物理找矿综合评价效果，因此亟需提供一种磁异常网格快速缝合方法。
发明内容
本发明的目的在于提供一种磁异常网格快速缝合方法，能够准确、迅速地拼接整合同一测区的磁异常数据。
实现本发明目的的技术方案：
一种磁异常网格快速缝合方法，包括以下步骤：
步骤1、在铀资源勘查区域内，划分A、B两区分别开展地面高精度磁测，在两区开展不同比例尺，但A、B区有部分的磁数据重合，分别形成magA和magB网格图；步骤2、对步骤1中的magA和magB进行化极处理，化极后形成gridA和gridB网格图；步骤3、对步骤2中gridA和gridB两张网格图进行趋势面去除处理，得到gridA-1和gridB-1网格图；步骤4、对步骤3中趋势面去除后的gridA-1和gridB-1网格开始进行缝合，缝合路径可在位于gridA-1和gridB-1网格图的重叠区域自定义3点连接而成；步骤5、选取gridA-1和gridB-1网格图校正的权重，0值表示所有校正都应用到网格magA中；值1表示所有校正都应用到网格magB中；值0.5表示两个网格图之间平均共享校正；步骤6、拓展gridA-1和gridB-1网格图的校正宽度，得到校正后网格图gridA-1-1和gridB-1-1；步骤7、将步骤6中校正后的gridA-1-1和gridB-1-1两张网格图进行插值，插值方法选取Akima样条插值方法，最大值间隔设置为0，所有间隔将被插值，最终形成缝合后的网格图gridC，查看两张网格图的缝合效果，若网格缝合处还出现空白数据，则回到步骤3，将去除趋势面参数选择2或3阶，重复步骤3-步骤7直至网格缝合处异常平滑，形成最终缝合后的网格图gridC。
所述步骤1中，插值方法采用最小曲率法，插值间距选取1/4个网格大小。
所述步骤2中，化极处理需查找地球基本磁场参数，通过球谐模型计算出该勘查区背景磁场强度值、磁倾角和磁偏角。
所述步骤3中，趋势面移除处理中可选择零阶或静态校正从整个网格图移除单一的平均值；或选择一阶校正移除最佳拟合的平面或斜率；或选择二阶或三阶校正拟合测点的二次或三次多项式，校正计算使用gridA和gridB网格图重叠的测点。
所述步骤4中，两张网格图沿缝合线的测点值通过调整路径两侧的网格图进行校正以使其相互匹配；与路径点相邻的点可能进行调整，以生成两个网格图之间的平滑过渡。
所述步骤4中，在每个缝合点将需要的校正在半径等于定义校正宽度的圆圈向外发散，这种方法减少了缝合路径的拐角和弯曲附近出现的失真数量，也保证了运算的速度。
所述步骤7中，Akima样条插值方法：任取6个实测点，设用i(i＝1,2,3,4,5,6)表示这6个实测数据点的序号，其坐标为(xi,yi),差值点(x，y)位于第3和第4个实测点之间，即x3<x<x4,则插值可用下式计算：
f(x，y)＝p0+p1(x-X3)+p2(x-X3)2+p3(x-X3)3 (1)
其中，p0＝y3，p1＝t3，p2＝[3(y4-y3)/(X4-X3)-2t3-t4]/(X4-X3),p3＝[t3+t4-2(y4-y3)/(X4-X3)]/(X4-X3)2,t3，t4分别是第3和第4号测点要素的斜率，他们分别用1,2,3,4,5和2,3,4,5,6点上的实测值表示。
本发明的有益技术效果在于：
(1)本发明提供的一种磁异常网格快速缝合方法，适用于全国范围铀资源地球物理预测或勘查区，可操作性强、有效性高。
(2)本发明提供的一种磁异常网格快速缝合方法，利用去除趋势面、自定义缝合线等将同一测区部分重合网格进行缝合拼接，方法流程操作简明，具有高效性，从而减少原始数据剔除、调整、改变投影方式等预处理工作量，达到快速地了解区内磁异常信息的功能。
具体实施方式
下面结合实施例对本发明作进一步详细说明。
以江西相山铀矿田为例，本发明涉及一种磁异常网格快速缝合方法，依次包括以下步骤：
步骤1、在选定铀资源勘查区域内，划分出A、B两区分别开展地面高精度磁测，A区比例尺为1:1万，B区比例尺为1:5000，且在两区布设测线方向、点距均不一致，但A、B区有部分重合的磁法测量数据，分别形成magA和magB网格图，插值方法采用最小曲率法，插值间距选取1/4个网格大小；
步骤2、对步骤1中的magA和magB进行化极处理，化极后形成gridA和gridB网格图；
所述步骤2中化极处理需查找地球基本磁场参数，通过利用http://www.ngdc.noaa.gov/seg/geomag/jsp/IGRF.jsp网站提供的球谐模型计算出该勘查区背景磁场强度值、磁倾角和磁偏角；
步骤3、对步骤2中gridA和gridB两张网格图进行趋势面去除处理，得到gridA-1和gridB-1网格图；
其中趋势面移除处理中可选择零阶或静态校正从整个网格图移除单一的平均值；或选择一阶校正移除最佳拟合的平面或斜率；或选择二阶或三阶校正拟合测点的二次或三次多项式，校正计算使用gridA和gridB网格图重叠的测点；
步骤4、对步骤3中趋势面去除后的gridA-1和gridB-1网格开始进行缝合，缝合路径可在位于gridA-1和gridB-1网格图的重叠区域自定义3点连接而成；
在上述步骤4中，其中两张网格图沿缝合线的测点值通过调整路径两侧的网格图进行校正以使其相互匹配，例如如果在缝合路径上gridA-1网格图的值比gridB-1网格图的值大0.1，可通过求取平均来消除这种差异。然后与路径点相邻的点可能进行调整，以生成两个网格图之间的平滑过渡。此种缝合方法不用考虑缝合路径所穿过特征的幅值和波长，并且确保了从gridA-1网格图到gridB-1的转换保持平滑；
在上述步骤4中，在每个缝合点将需要的校正在半径等于定义“校正宽度”的圆圈向外发散，这种方法减少了缝合路径的拐角和弯曲附近出现的失真数量，也保证了运算的速度；
步骤5、选取gridA-1和gridB-1网格图校正的权重，0值表示所有校正都应用到网格magA中；值1表示所有校正都应用到网格magB中；值0.5表示两个网格图之间平均共享校正，本发明实例中选取权重值为0.5；
步骤6、拓展gridA-1和gridB-1网格图的校正宽度，得到校正后网格图gridA-1-1和gridB-1-1，其中选取校正宽度值为32,校正宽度值越高则缝合需要时间越长，因为受影响的网格面积更大；
步骤7、将步骤6中校正后的gridA-1-1和gridB-1-1两张网格图进行插值，插值方法选取Akima样条插值方法，最大值间隔设置为0，所有间隔将被插值，最终形成缝合后的网格图gridC，查看两张网格图的缝合效果，若网格缝合处还出现空白数据，则回到步骤3，将去除趋势面参数选择2或3阶，重复步骤3-步骤7直至网格缝合处异常平滑，形成最终缝合后的网格图gridC。
其中Akima样条插值方法：任取6个实测点，设用i(i＝1,2,3,4,5,6)表示这6个实测数据点的序号，其坐标为(xi,yi),差值点(x，y)位于第3和第4个实测点之间，即x3<x<x4,则插值可用下式计算：
f(x，y)＝p0+p1(x-X3)+p2(x-X3)2+p3(x-X3)3 (1)
其中，p0＝y3，p1＝t3，p2＝[3(y4-y3)/(X4-X3)-2t3-t4]/(X4-X3),p3＝[t3+t4-2(y4-y3)/(X4-X3)]/(X4-X3)2,t3，t4分别是第3和第4号测点要素的斜率，他们分别用1,2,3,4,5和2,3,4,5,6点上的实测值表示。
上面实施例对本发明作了详细说明，但是本发明并不限于上述实施例，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下做出各种变化。本发明中未作详细描述的内容均可以采用现有技术。

Claims (7)

1.一种磁异常网格快速缝合方法，其特征在于：包括以下步骤：
步骤(1)、在铀资源勘查区域内，划分A、B两区分别开展地面高精度磁测，在两区开展不同比例尺，但A、B区有部分的磁数据重合，分别形成magA和magB网格图；步骤(2)、对步骤(1)中的magA和magB进行化极处理，化极后形成gridA和gridB网格图；步骤(3)、对步骤(2)中gridA和gridB两张网格图进行趋势面去除处理，得到gridA-1和gridB-1网格图；步骤(4)、对步骤(3)中趋势面去除后的gridA-1和gridB-1网格开始进行缝合，缝合路径可在位于gridA-1和gridB-1网格图的重叠区域自定义3点连接而成；步骤(5)、选取gridA-1和gridB-1网格图校正的权重，0值表示所有校正都应用到网格magA中；值1表示所有校正都应用到网格magB中；值0.5表示两个网格图之间平均共享校正；步骤(6)、拓展gridA-1和gridB-1网格图的校正宽度，得到校正后网格图gridA-1-1和gridB-1-1；步骤(7)、将步骤(6)中校正后的gridA-1-1和gridB-1-1两张网格图进行插值，插值方法选取Akima样条插值方法，最大值间隔设置为0，所有间隔将被插值，最终形成缝合后的网格图gridC，查看两张网格图的缝合效果，若网格缝合处还出现空白数据，则回到步骤(3)，将去除趋势面参数选择2或3阶，重复步骤(3)-步骤(7)直至网格缝合处异常平滑，形成最终缝合后的网格图gridC。
2.根据权利要求1所述的一种磁异常网格快速缝合方法，其特征在于：所述步骤(1)中，插值方法采用最小曲率法，插值间距选取1/4个网格大小。
3.根据权利要求1所述的一种磁异常网格快速缝合方法，其特征在于：所述步骤(2)中，化极处理需查找地球基本磁场参数，通过球谐模型计算出该勘查区背景磁场强度值、磁倾角和磁偏角。
4.根据权利要求1所述的一种磁异常网格快速缝合方法，其特征在于：所述步骤(3)中，趋势面移除处理中可选择零阶或静态校正从整个网格图移除单一的平均值；或选择一阶校正移除最佳拟合的平面或斜率；或选择二阶或三阶校正拟合测点的二次或三次多项式，校正计算使用gridA和gridB网格图重叠的测点。
5.根据权利要求4所述的一种磁异常网格快速缝合方法，其特征在于：所述步骤(4)中，两张网格图沿缝合线的测点值通过调整路径两侧的网格图进行校正以使其相互匹配；与路径点相邻的点可能进行调整，以生成两个网格图之间的平滑过渡。
6.根据权利要求4所述的一种磁异常网格快速缝合方法，其特征在于：所述步骤(4)中，在每个缝合点将需要的校正在半径等于定义校正宽度的圆圈向外发散，这种方法减少了缝合路径的拐角和弯曲附近出现的失真数量，也保证了运算的速度。
7.根据权利要求1所述的一种磁异常网格快速缝合方法，其特征在于：所述步骤(7)中，Akima样条插值方法：任取6个实测点，设用i(i＝1,2,3,4,5,6)表示这6个实测数据点的序号，其坐标为(xi,yi),差值点(x，y)位于第3和第4个实测点之间，即x3<x<x4,则插值可用下式计算：
f(x，y)＝p0+p1(x-X3)+p2(x-X3)2+p3(x-X3)3 (1)
其中，p0＝y3，p1＝t3，p2＝[3(y4-y3)/(X4-X3)-2t3-t4]/(X4-X3),p3＝[t3+t4-2(y4-y3)/(X4-X3)]/(X4-X3)2,t3，t4分别是第3和第4号测点要素的斜率，他们分别用1,2,3,4,5和2,3,4,5,6点上的实测值表示。
CN201911407381.4A 2019-12-31 2019-12-31 一种磁异常网格快速缝合方法 Active CN111045090B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911407381.4A CN111045090B (zh) 2019-12-31 2019-12-31 一种磁异常网格快速缝合方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911407381.4A CN111045090B (zh) 2019-12-31 2019-12-31 一种磁异常网格快速缝合方法

Publications (2)

Publication Number Publication Date
CN111045090A true CN111045090A (zh) 2020-04-21
CN111045090B CN111045090B (zh) 2022-05-20

Family

ID=70241546

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911407381.4A Active CN111045090B (zh) 2019-12-31 2019-12-31 一种磁异常网格快速缝合方法

Country Status (1)

Country Link
CN (1) CN111045090B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070159524A1 (en) * 2006-01-09 2007-07-12 Samsung Electronics Co., Ltd. Method and apparatus for providing panoramic view with high speed image matching and mild mixed color blending
CN103592691A (zh) * 2013-10-30 2014-02-19 中国石油天然气集团公司 一种反向剩余磁性磁异常识别方法及装置
US20160125555A1 (en) * 2014-10-31 2016-05-05 Larisa V. Branets Handling Domain Discontinuity With The Help of Grid Optimization Techniques
CN106981093A (zh) * 2017-05-17 2017-07-25 上海同岩土木工程科技股份有限公司 一种分区约束耦合的三维地层并行建模方法
EP3246871A1 (en) * 2016-05-20 2017-11-22 Shenyang Neusoft Medical Systems Co., Ltd. Image splicing
CN109544447A (zh) * 2018-10-26 2019-03-29 广西师范大学 一种图像拼接方法、装置及存储介质
CN109636913A (zh) * 2018-12-04 2019-04-16 山东理工大学 基于Delaunay剖分的三角网格增量拓扑拼接方法
CN110211234A (zh) * 2019-05-08 2019-09-06 上海索辰信息科技有限公司 一种网格模型缝合系统和方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070159524A1 (en) * 2006-01-09 2007-07-12 Samsung Electronics Co., Ltd. Method and apparatus for providing panoramic view with high speed image matching and mild mixed color blending
CN103592691A (zh) * 2013-10-30 2014-02-19 中国石油天然气集团公司 一种反向剩余磁性磁异常识别方法及装置
US20160125555A1 (en) * 2014-10-31 2016-05-05 Larisa V. Branets Handling Domain Discontinuity With The Help of Grid Optimization Techniques
EP3246871A1 (en) * 2016-05-20 2017-11-22 Shenyang Neusoft Medical Systems Co., Ltd. Image splicing
CN106981093A (zh) * 2017-05-17 2017-07-25 上海同岩土木工程科技股份有限公司 一种分区约束耦合的三维地层并行建模方法
CN109544447A (zh) * 2018-10-26 2019-03-29 广西师范大学 一种图像拼接方法、装置及存储介质
CN109636913A (zh) * 2018-12-04 2019-04-16 山东理工大学 基于Delaunay剖分的三角网格增量拓扑拼接方法
CN110211234A (zh) * 2019-05-08 2019-09-06 上海索辰信息科技有限公司 一种网格模型缝合系统和方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
尹航 等: "编制１：500万中国海陆磁力异常图关键技术", 《地球物理学进展》 *
杨柳 等: "Geosoft环境下两种网格拼接的方法", 《数字技术与应用》 *
陈上及 等编著: "《海洋数据处理分析方法及其应用》", 31 August 1991, 海洋出版社 *

Also Published As

Publication number Publication date
CN111045090B (zh) 2022-05-20

Similar Documents

Publication Publication Date Title
US10459117B2 (en) Extended subspace method for cross-talk mitigation in multi-parameter inversion
US20170176613A1 (en) Method To Design Geophysical Surveys Using Full Wavefield Inversion Point- Spread Function Analysis
US20140278120A1 (en) Methods and systems for locating seismic events
EP2877880B1 (en) System and method for migration velocity modeling
CN104237940A (zh) 一种基于动力学特征的绕射波成像方法及装置
CN105319589A (zh) 一种利用局部同相轴斜率的全自动立体层析反演方法
CN105353405A (zh) 一种全波形反演方法和系统
CN103499836A (zh) 空变-多时窗融合高精度剩余静校正方法
CN106959467A (zh) 地震波速度反演方法和装置
CN116774292A (zh) 一种地震波走时确定方法、系统、电子设备及存储介质
CN113671570B (zh) 一种地震面波走时和重力异常联合反演方法与系统
CN111045090B (zh) 一种磁异常网格快速缝合方法
WO2024078134A1 (zh) 基于多参数约束和结构校正的随掘巷道全波形反演方法
CN109581494B (zh) 叠前偏移方法及装置
CN106054252A (zh) 一种叠前时间偏移的方法及装置
Hardy High resolution tomographic MVA with automation
CN111435172A (zh) 层析静校正的方法和装置
CN114428341A (zh) 振幅保真地震成像方法、装置、电子设备及介质
CN111751871B (zh) 三维地震沿层相干切片处理方法、装置及电子设备
CN104166162B (zh) 基于迭代三参数小波变换的缝洞发育带检测方法
CN113267810B (zh) 地震勘探全深度速度建模方法及装置
CN116626751B (zh) 基于多目标函数的黏弹性参数同步反演方法、装置和设备
CN112130198B (zh) 基于时差校正的多期次三维地震资料融合方法及装置
CN113514887B (zh) 确定黄土山地区地震资料共中心点参考面的方法及装置
US7295929B2 (en) Process, device and software package for locally determining the shape of geological horizons

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant