CN113779013B - 基于地质约束的小范围磁场数据补缺处理方法 - Google Patents
基于地质约束的小范围磁场数据补缺处理方法 Download PDFInfo
- Publication number
- CN113779013B CN113779013B CN202111097134.6A CN202111097134A CN113779013B CN 113779013 B CN113779013 B CN 113779013B CN 202111097134 A CN202111097134 A CN 202111097134A CN 113779013 B CN113779013 B CN 113779013B
- Authority
- CN
- China
- Prior art keywords
- magnetic
- data
- geological
- measurement
- extracting
- 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.)
- Active
Links
- 230000007812 deficiency Effects 0.000 title claims abstract description 9
- 238000003672 processing method Methods 0.000 title claims abstract description 8
- 238000005259 measurement Methods 0.000 claims abstract description 65
- 238000000034 method Methods 0.000 claims abstract description 10
- 238000013075 data extraction Methods 0.000 claims description 35
- 238000000605 extraction Methods 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 238000004519 manufacturing process Methods 0.000 claims description 3
- 230000006698 induction Effects 0.000 claims 1
- 238000011161 development Methods 0.000 abstract description 7
- 239000011435 rock Substances 0.000 abstract description 6
- 238000011835 investigation Methods 0.000 abstract description 5
- 230000007613 environmental effect Effects 0.000 abstract description 4
- 238000004364 calculation method Methods 0.000 abstract description 3
- 239000000284 extract Substances 0.000 description 4
- 230000001502 supplementing effect Effects 0.000 description 2
- YTAHJIFKAKIKAV-XNMGPUDCSA-N [(1R)-3-morpholin-4-yl-1-phenylpropyl] N-[(3S)-2-oxo-5-phenyl-1,3-dihydro-1,4-benzodiazepin-3-yl]carbamate Chemical compound O=C1[C@H](N=C(C2=C(N1)C=CC=C2)C1=CC=CC=C1)NC(O[C@H](CCN1CCOCC1)C1=CC=CC=C1)=O YTAHJIFKAKIKAV-XNMGPUDCSA-N 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/21—Design, administration or maintenance of databases
- G06F16/215—Improving data quality; Data cleansing, e.g. de-duplication, removing invalid entries or correcting typographical errors
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
Abstract
本发明涉及一种基于地质约束的小范围磁场数据补缺处理方法，通过收集工作区地质及磁性参数资料，对工作区磁测数据进行提取、处理，并提取磁测空白区范围。根据空白区内不同地层或岩体的分布及其磁性特征，建立基于地质约束的磁性体模型，结合磁测空白区边缘数据，根据磁性体模型开展基于边缘磁测数据约束的磁正演计算，然后逐条线拟合磁剖面曲线。提取拟合磁剖面数据，与各测区磁测数据网格进行拼接，形成完整磁测数据，作为进一步磁资料开发利用的基础资料。本发明解决了部分地区小范围缺少磁数据的问题，在不重新开展磁测的情况下提供一种填补磁数据空白的快速有效技术方法，从而节约成本并减少因实地勘查带来的环境破坏。
Description
技术领域
本发明涉及一种磁场数据补缺方法，具体地说是一种基于地质约束的小范围磁场数据补缺处理方法。
背景技术
一些地区开展过不同年代、不同测量比例尺以及不同采集方式的磁测工作，但是局部地区属于磁测工作的空白区，由于磁测空白区的存在，造成在区域磁资料综合利用开发时面临困难。由于实地采集磁测数据成本也较高，急需探索快速有效且成本低的磁数据处理方法技术。随着已有磁测数据越来越多，对磁资料的融合开发越来越迫切，进一步研究磁数据补缺的方法，对我国已有磁资料的开发、降低国家财政支出及减少因实地勘查带来的环境破坏具有重要意义。
发明内容
本发明的目的就是提供一种基于地质约束的小范围磁场数据补缺处理方法，以解决现在存在磁测数据空白区以及实地勘查成本高、会破坏环境的问题。
本发明是这样实现的：一种基于地质约束的小范围磁场数据补缺处理方法，包括以下步骤。
a.收集工作区的地质图，并对工作区地质图进行矢量化处理。
b.收集工作区及工作区周边的磁性参数资料，统计磁性参数的测量地点、测量数量、磁化率和剩余磁化率，并将各磁性参数的点坐标投影到步骤a中得到的矢量化地质图上。
c.对工作区磁测数据进行提取整理，并提取磁测空白区范围坐标。
d.在空白区布设数据提取线，将数据提取线投影到步骤b中得到的地质图上，分析经过数据提取线的不同地质体的分布状态，沿数据提取线提取地质节点，制作沿提取线的地质断面。
e.提取经过数据提取线的地质体的磁性参数，对地质断面上不同地质体进行磁性参数赋值，建立基于地质约束的磁性体断面模型。
f.沿数据提取线提取分布于空白区边缘的已有磁测数据，形成边缘磁数据，然后将边缘磁数据导入磁正反演软件中，导入步骤e中建立的基于地质约束的磁性体断面模型，开展磁性模型正演。
g.在基于地质约束的磁性体断面模型上，不断调整模型端面上不同地质体的深度及空间范围，并对不同地质体的磁性参数进行微调，直到基于地质约束的磁性体断面模型的拟合正演磁数据曲线的两端与边缘磁数据重合。
h.对所有的数据提取线逐条进行建模，并根据模型正演得到拟合磁剖面曲线数据。
i.提取每条数据提取线上的拟合磁剖面数据及其坐标，形成明码格式文件，并进行网格化。
j.对步骤c中的工作区磁测数据进行网格化，并与步骤i中得到的网格化数据进行拼接，形成完整的网格化磁测数据。
进一步地，在步骤c中，对收集的工作区磁测数据进行整理，提取并投影到同一坐标系上，形成可供进一步利用的明码文件；并提取前人开展的磁法项目工作区范围，编制工作程度图件，圈出磁测空白区范围并提取磁测空白区的坐标。
进一步地，在步骤d中，数据提取线的方向尽量与周围磁测数据的测线方向一致，并重复500m~1000m，数据提取线的线间距为200~500m。
进一步地，在步骤e中，进行赋值的磁性参数包括感应磁化率、剩余磁化率、磁偏角与磁倾角。
进一步地，在步骤i和j中，网格化间距取数据提取线间距的1/4。
进一步地，在步骤j中，采用缝合法对网格化数据进行拼接。
进一步地，在步骤b中，统计磁性参数的测量地点、测量数量、磁化率和剩余磁化率。
本发明通过收集工作区地质及磁性参数资料，对工作区磁测数据进行提取、处理，并提取磁测空白区范围。通过分析磁测空白区地质特征，根据空白区内不同地层或岩体的分布及其磁性特征，建立基于地质约束的磁性体模型，结合磁测空白区边缘数据，根据磁性体模型开展基于边缘磁测数据约束的磁正演计算，然后逐条线拟合磁剖面曲线。提取拟合磁剖面数据，形成网格文件，与各测区磁测数据网格进行拼接，形成完整磁测数据，作为进一步磁资料开发利用的基础资料。本发明解决了部分地区小范围缺少磁数据的问题，在不重新开展磁测的情况下提供一种填补磁数据空白的快速有效技术方法，从而节约成本并减少因实地勘查带来的环境破坏。
附图说明
图1是本发明的流程图。
具体实施方式
如图1所示，下面结合附图对本发明作进一步详细说明。
a.收集工作区的地质图，并对工作区地质图进行矢量化处理。
具体的地，所收集的地质图为大比例地质图，对地质图进行矢量化后，为后续提取空白区不同地层或岩体的界限坐标提供基础。
b.收集工作区及工作区周边的磁性参数资料，统计磁性参数的测量地点、测量数量、磁化率和剩余磁化率等参数，并将各磁性参数的点坐标投影到步骤a中得到的矢量化地质图上。
c.对工作区磁测数据进行提取整理，并提取磁测空白区范围坐标。
具体地，对收集的工作区磁测数据进行整理，提取并投影到同一坐标系上，形成可供进一步利用的明码文件；并提取前人开展的磁法项目工作区范围，编制工作程度图件，圈出磁测空白区范围并提取磁测空白区的坐标。
d.在空白区布设数据提取线，将数据提取线投影到步骤b中得到的地质图上，分析经过数据提取线的不同地质体的分布状态，沿数据提取线提取地质节点，制作沿提取线的地质断面。
其中，数据提取线的方向尽量与周围磁测数据的测线方向一致，并重复500m~1000m，数据提取线的线间距为200~500m，具体根据周边测量比例尺确定，保持大致一致的原则。
e.提取经过数据提取线的地质体的磁性参数，对地质断面上不同地质体进行磁性参数赋值，建立基于地质约束的磁性体断面模型。
其中，进行赋值的磁性参数包括感应磁化率、剩余磁化率、磁偏角与磁倾角。
f.沿数据提取线提取分布于空白区边缘外围的已有磁测数据，形成边缘磁数据，然后将边缘磁数据导入磁正反演软件中，导入步骤e中建立的基于地质约束的磁性体断面模型，开展磁性模型正演。
g.在基于地质约束的磁性体断面模型上，不断调整模型端面上不同地质体的深度及空间范围，并对不同地质体的磁性参数进行微调，直到基于地质约束的磁性体断面模型的拟合正演磁数据曲线的两端与边缘磁数据重合。
h.对所有的数据提取线逐条进行建模，并根据模型正演得到拟合磁剖面曲线数据。
i.提取每条数据提取线上的拟合磁剖面数据及其坐标，形成明码格式文件，并进行网格化，网格化间距取数据提取线间距的1/4。
j.对步骤c中的工作区磁测数据进行网格化，网格化间距取数据提取线间距的1/4，并与步骤i中得到的网格化数据进行拼接，采用缝合法对网格化数据进行拼接，形成完整的网格化磁测数据，作为进一步磁资料开发利用的基础。
本发明通过收集工作区地质及磁性参数资料，对工作区磁测数据进行提取、处理，并提取磁测空白区范围。通过分析磁测空白区地质特征，根据空白区内不同地层或岩体的分布及其磁性特征，建立基于地质约束的磁性体模型，结合磁测空白区边缘数据，根据磁性体模型开展基于边缘磁测数据约束的磁正演计算，然后逐条线拟合磁剖面曲线。提取拟合磁剖面数据，形成网格文件，与各测区磁测数据网格进行拼接，形成完整磁测数据，作为进一步磁资料开发利用的基础资料。本发明解决了部分地区小范围缺少磁数据的问题，在不重新开展磁测的情况下提供一种填补磁数据空白的快速有效技术方法，从而节约成本并减少因实地勘查带来的环境破坏。
Claims (2)
1.一种基于地质约束的小范围磁场数据补缺处理方法，其特征在于，包括以下步骤：
a.收集工作区的地质图，并对工作区地质图进行矢量化处理；
b.收集工作区及工作区周边的磁性参数资料，统计磁性参数的测量地点、测量数量、磁化率和剩余磁化率，并将各磁性参数的点坐标投影到步骤a中得到的矢量化地质图上；
c.对工作区磁测数据进行提取整理，并提取磁测空白区范围坐标：对收集的工作区磁测数据进行整理，提取并投影到同一坐标系上，形成可供进一步利用的明码文件；并提取前人开展的磁法项目工作区范围，编制工作程度图件，圈出磁测空白区范围并提取磁测空白区的坐标；
d.在空白区布设数据提取线，将数据提取线投影到步骤b中得到的地质图上，分析经过数据提取线的不同地质体的分布状态，沿数据提取线提取地质节点，制作沿提取线的地质断面，数据提取线的方向尽量与周围磁测数据的测线方向一致，并重复500m~1000m，数据提取线的线间距为200~500m；
e.提取经过数据提取线的地质体的磁性参数，对地质断面上不同地质体进行磁性参数赋值，建立基于地质约束的磁性体断面模型，进行赋值的磁性参数包括感应磁化率、剩余磁化率、磁偏角与磁倾角；
f.沿数据提取线提取分布于空白区边缘的已有磁测数据，形成边缘磁数据，然后将边缘磁数据导入磁正反演软件中，导入步骤e中建立的基于地质约束的磁性体断面模型，开展磁性模型正演；
g.在基于地质约束的磁性体断面模型上，不断调整模型端面上不同地质体的深度及空间范围，并对不同地质体的磁性参数进行微调，直到基于地质约束的磁性体断面模型的拟合正演磁数据曲线的两端与边缘磁数据重合；
h.对所有的数据提取线逐条进行建模，并根据模型正演得到拟合磁剖面曲线数据；
i.提取每条数据提取线上的拟合磁剖面数据及其坐标，形成明码格式文件，并进行网格化；
j.对步骤c中的工作区磁测数据进行网格化，并与步骤i中得到的网格化数据进行拼接，采用缝合法对网格化数据进行拼接，形成完整的网格化磁测数据。
2.根据权利要求1所述的基于地质约束的小范围磁场数据补缺处理方法，其特征在于，在步骤i和j中，网格化间距取数据提取线间距的1/4。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111097134.6A CN113779013B (zh) | 2021-09-18 | 2021-09-18 | 基于地质约束的小范围磁场数据补缺处理方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111097134.6A CN113779013B (zh) | 2021-09-18 | 2021-09-18 | 基于地质约束的小范围磁场数据补缺处理方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113779013A CN113779013A (zh) | 2021-12-10 |
| CN113779013B true CN113779013B (zh) | 2024-01-30 |
Family
ID=78852211
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111097134.6A Active CN113779013B (zh) | 2021-09-18 | 2021-09-18 | 基于地质约束的小范围磁场数据补缺处理方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113779013B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2978500A1 (en) * | 2015-03-04 | 2016-09-09 | Institute of Mineral Resources, Chinese Academy of Geological Sciences | Method for automatically extracting tectonic framework of potential field data |
| CN109799540A (zh) * | 2018-12-25 | 2019-05-24 | 核工业北京地质研究院 | 基于地质信息约束下的火山岩型铀矿床磁化率反演方法 |
| CN110657798A (zh) * | 2019-09-11 | 2020-01-07 | 广州海洋地质调查局 | 一种地磁日变数据补缺方法及处理终端 |
| CN113050182A (zh) * | 2021-03-12 | 2021-06-29 | 安徽省勘查技术院(安徽省地质矿产勘查局能源勘查中心) | 一种水域地磁场观测方法及系统 |
-
2021
- 2021-09-18 CN CN202111097134.6A patent/CN113779013B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2978500A1 (en) * | 2015-03-04 | 2016-09-09 | Institute of Mineral Resources, Chinese Academy of Geological Sciences | Method for automatically extracting tectonic framework of potential field data |
| CN109799540A (zh) * | 2018-12-25 | 2019-05-24 | 核工业北京地质研究院 | 基于地质信息约束下的火山岩型铀矿床磁化率反演方法 |
| CN110657798A (zh) * | 2019-09-11 | 2020-01-07 | 广州海洋地质调查局 | 一种地磁日变数据补缺方法及处理终端 |
| CN113050182A (zh) * | 2021-03-12 | 2021-06-29 | 安徽省勘查技术院(安徽省地质矿产勘查局能源勘查中心) | 一种水域地磁场观测方法及系统 |
Non-Patent Citations (1)
| Title |
|---|
| 基于凸集投影方法的重磁数据规则缺失重建;闫浩飞;刘国峰;薛典军;王林飞;;地球物理学进展(05);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113779013A (zh) | 2021-12-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102609986B (zh) | 数字矿山中矿体建模系统及其建模和检测方法 | |
| CN109508508B (zh) | 一种露天矿山治理勘查设计方法 | |
| CN102750413B (zh) | 一种输电线路塔位地形测量的数据处理及成图方法 | |
| CN110276732B (zh) | 一种顾及地形特征线要素的山区点云空洞修复方法 | |
| CN102609982B (zh) | 空间地质数据非结构化模式的拓扑发现方法 | |
| CN107103153A (zh) | 一种基于三维激光扫描技术的矿产资源消耗量评估方法 | |
| Wu et al. | Three-dimensional geological modeling and its application in digital mine | |
| CN111339691A (zh) | 一种基于voxler软件的智慧岩土工程参数三维分析与评价系统及方法 | |
| CN106709987B (zh) | 一种三维地质剖面模型动态构建方法 | |
| CN110070038B (zh) | 一种基于城市绿化覆盖面积测算城市绿地率的方法 | |
| CN101915570A (zh) | 一种基于灭点的地面移动测量影像线段自动提取分类方法 | |
| CN115560690B (zh) | 一种基于三维激光扫描技术的结构物整体变形分析方法 | |
| CN112906096B (zh) | 一种公路工程设计数据重构解析的方法及系统 | |
| CN101964009A (zh) | 一种基于insar制作3d产品的系统及方法 | |
| CN111191307A (zh) | 一种基于bim+gis技术的土方工程虚拟施工方法 | |
| CN113467315A (zh) | 一种基于bim技术的隧道工程自动化监测控制方法及系统 | |
| Ni et al. | Lineament length and density analyses based on the segment tracing algorithm: a case study of the gaosong field in gejiu tin mine, China | |
| CN110197013B (zh) | 基于Morphing的河床基岩面建模方法 | |
| CN113779013B (zh) | 基于地质约束的小范围磁场数据补缺处理方法 | |
| CN109711680A (zh) | 基于gis数字化技术的河流水电规划方法 | |
| CN105205311A (zh) | 一种地层剖面绘制的复数运算法 | |
| Lin et al. | A simple generation technique of complex geotechnical computational model | |
| Gao et al. | Mathematical Interpolation and Correction of Three-Dimensional Modelling of High-Speed Railway. | |
| CN114511240A (zh) | 一种矿山采掘计划生成方法及存储介质 | |
| CN114492093A (zh) | 一种城市可视化方法及系统 |
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 |