CN108508479A - 一种空地井立体重磁数据协同目标位置反演方法 - Google Patents

一种空地井立体重磁数据协同目标位置反演方法 Download PDF

Info

Publication number
CN108508479A
CN108508479A CN201810185427.1A CN201810185427A CN108508479A CN 108508479 A CN108508479 A CN 108508479A CN 201810185427 A CN201810185427 A CN 201810185427A CN 108508479 A CN108508479 A CN 108508479A
Authority
CN
China
Prior art keywords
gravity
data
magnetic data
magnetic
inverting
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
CN201810185427.1A
Other languages
English (en)
Other versions
CN108508479B (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.)
Jilin University
Original Assignee
Jilin University
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 Jilin University filed Critical Jilin University
Priority to CN201810185427.1A priority Critical patent/CN108508479B/zh
Publication of CN108508479A publication Critical patent/CN108508479A/zh
Application granted granted Critical
Publication of CN108508479B publication Critical patent/CN108508479B/zh
Expired - Fee Related 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
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging

Landscapes

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

Abstract

本发明公开了一种空地井立体重磁数据协同目标位置反演方法，所述空地井立体重磁数据协同目标位置反演方法通过采用同一目标函数，建立数据间的关系方程，依据不同空间数据的精度设定反演的权系数，从而获得地下物质的准确空间分布；获取同一地质目标的不同高度重磁响应的基础上，将立体重磁数据协同正反演，能有效地降低反演的多解性；该发明设计合理，思路清晰，操作方便，反演结果准确度高可有效地降低勘探成本，采用空地井协同正反演可有效地降低结果的多解性，从而提供更加准确的结果，准确地划定下一步勘探的靶区，大大地降低勘探成本，安全稳定，适用范围广，有利于推广和普及。

Description

一种空地井立体重磁数据协同目标位置反演方法
技术领域
本发明属于地球科学技术领域，更具体地说，尤其涉及一种空地井立体重磁数据协同目标位置反演方法。
背景技术
针对难以进入区域的地球物理勘探工作，以及提高地质体分辨率和降低解释结果多解性的目标，在已有地面重磁勘探的基础上，近几年航勘探和井中勘探成为新型重磁勘探手段，因此实现了空-地-井立体重磁数据的获取，重磁勘探具有大范围、速度快、经费低的优势，是所有地球物理工作中首先勘探的，其数据的解释结果为下一步地震、井等精细、昂贵地球物理勘探方法提供勘探靶区，其结果的准确性直接决定资源靶区的准确性，如果结果的不准确，会造成精细勘探的浪费，以往地面重磁勘探采用单一数据进行地下目标的反演，结果具有明显的多解性和不确定行。针对不同高度的重磁数据，人们往往采用的数据比值来进行反演，或者建立各自的反演方程，不能很好地综合利用多数据的优势，不利于广泛的推广和普及。
发明内容
本发明的目的是为了解决现有技术中存在的缺点，而提出的一种空地井立体重磁数据协同目标位置反演方法。
为实现上述目的，本发明提供如下技术方案：一种空地井立体重磁数据协同目标位置反演方法，包括如下步骤：
S1、首先在需要勘测地点的地平面计算地面重磁数据，并利用泰勒展开式方程来计算数据精度ε2；
S2、然后利用无人机计算航空重磁数据和下井测试井中重磁数据，并分别利用泰勒展开式方程来计算数据精度ε1和ε3；
S3、再对单一空间重磁数据进行地质体位置反演，并且根据不同的高度和深入重复步骤S2和步骤S3；
S4、根据不同数据的质量分配不同数据权重W，并根据多空间重磁数据按照相应公式进行反演；
S5、最后将单层反演窗口缩小而利用其它空间数据形成补充，重复上述操作根据不同数据的精度和根据精度设定相应的权重系数，进而组合为矩阵进行反演，最后得出地质体空间位置信息。
优选的，所述步骤S1和S2中勘测重磁数据的泰勒展开式方程关系式为：
其中，T0(x，y，0)和T(x，y，h)表示高差为h的两个观测面上的异常，当h＞0时表示向上延拓，当h＜0时表示向下延拓，对公式进行Fourier变换可得到
其中，表示异常T0(x，y，0)的Fourier谱，(kx，ky)是分别表示x和y方向上的波数，称为延拓因子。
优选的，所述步骤S3计算单一空间重磁数据进行地质体位置反演时的原理是重磁数据及其导数满足欧拉齐次方程：
(x-x0)·fx+(y-y0)·fy+(z-z0)·fz＝-N(f(x-x0，y-y0，z-z0)+A)
其中f为原始重磁异常，fx，fy，fz分别为异常f在不同方向的导数，x，y，z为观测点坐标，A为背景异，N为构造指数(人为给定)，x0，y0，z0为需反演的地质体位置坐标。
优选的，所述步骤S4中对多空间重磁数据采用方程进行位置反演时的表达式为：
其中，h为代表观测高度。
优选的，所述步骤S5中计算不同数据的精度，根据精度设定相应的权重系数，进而组合为矩阵进行反演的公式为：
本发明的技术效果和优点：本发明提供的一种空地井立体重磁数据协同目标位置反演方法，与传统的反演方法相比，本发明以采用同一目标函数，建立数据间的关系方程，依据不同空间数据的精度设定反演的权系数，从而获得地下物质的准确空间分布；获取同一地质目标的不同高度重磁响应的基础上，将立体重磁数据协同正反演，能有效地降低反演的多解性；空地井立体重磁数据反演技术主要考虑多空间数据联合的优势以及方式，同时在反演时可以将单层反演窗口缩小而利用其它空间数据形成补充，不会降低参与约束反演的数据量，窗口的缩小增强正下方目标体的响应，有效地避免了周围地质体的干扰，从而提高对地质体的分辨率和准确性；此外，窗口的缩小可有效地提高水平叠加异常的分辨率，可有效地提高现有反演技术的分辨率和准确性；该发明设计合理，思路清晰，操作方便，反演结果准确度高可有效地降低勘探成本，采用空地井协同正反演可有效地降低结果的多解性，从而提供更加准确的结果，准确地划定下一步勘探的靶区，大大地降低勘探成本，安全稳定，适用范围广，有利于推广和普及。
附图说明
图1为本发明的不同空间重磁数据的联合反演流程图；
图2为本发明的空地井联合反演窗口变化对比图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1-2所示的一种空地井立体重磁数据协同目标位置反演方法，包括如下步骤：
S1、首先在需要勘测地点的地平面计算地面重磁数据，并利用泰勒展开式方程来计算数据精度ε2；
S2、然后利用无人机计算航空重磁数据和下井测试井中重磁数据，并分别利用泰勒展开式方程来计算数据精度ε1和ε3；
S3、再对单一空间重磁数据进行地质体位置反演，并且根据不同的高度和深入重复步骤S2和步骤S3；
S4、根据不同数据的质量分配不同数据权重W，并根据多空间重磁数据按照相应公式进行反演；
S5、最后将单层反演窗口缩小而利用其它空间数据形成补充，重复上述操作根据不同数据的精度和根据精度设定相应的权重系数，进而组合为矩阵进行反演，最后得出地质体空间位置信息。
优选的，所述步骤S1和S2中勘测重磁数据的泰勒展开式方程关系式为：
其中，T0(x，y，0)和T(x，y，h)表示高差为h的两个观测面上的异常，当h＞0时表示向上延拓，当h＜0时表示向下延拓，对公式进行Fourier变换可得到
其中，表示异常T0(x，y，0)的Fourier谱，(kx，ky)是分别表示x和y方向上的波数，称为延拓因子。
优选的，所述步骤S3计算单一空间重磁数据进行地质体位置反演时的原理是重磁数据及其导数满足欧拉齐次方程：
(x-x0)·fx+(y-y0)·fy+(z-z0)·fz＝-N(f(x-x0，y-y0，z-z0)+A)
其中f为原始重磁异常，fx，fy，fz分别为异常f在不同方向的导数，x，y，z为观测点坐标，A为背景异，N为构造指数(人为给定)，x0，y0，z0为需反演的地质体位置坐标。
优选的，所述步骤S4中对多空间重磁数据采用方程进行位置反演时的表达式为：
其中，h为代表观测高度。
优选的，所述步骤S5中计算不同数据的精度，根据精度设定相应的权重系数，进而组合为矩阵进行反演的公式为：
综上所述：本发明提供的一种空地井立体重磁数据协同目标位置反演方法，与传统的反演方法相比，本发明以采用同一目标函数，建立数据间的关系方程，依据不同空间数据的精度设定反演的权系数，从而获得地下物质的准确空间分布；获取同一地质目标的不同高度重磁响应的基础上，将立体重磁数据协同正反演，能有效地降低反演的多解性；空地井立体重磁数据反演技术主要考虑多空间数据联合的优势以及方式，同时在反演时可以将单层反演窗口缩小而利用其它空间数据形成补充，不会降低参与约束反演的数据量，窗口的缩小增强正下方目标体的响应，有效地避免了周围地质体的干扰，从而提高对地质体的分辨率和准确性；此外，窗口的缩小可有效地提高水平叠加异常的分辨率，可有效地提高现有反演技术的分辨率和准确性；该发明设计合理，思路清晰，操作方便，反演结果准确度高可有效地降低勘探成本，采用空地井协同正反演可有效地降低结果的多解性，从而提供更加准确的结果，准确地划定下一步勘探的靶区，大大地降低勘探成本，安全稳定，适用范围广，有利于推广和普及。
最后应说明的是：以上所述仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (5)

1.一种空地井立体重磁数据协同目标位置反演方法，其特征在于，包括如下步骤：
S1、首先在需要勘测地点的地平面计算地面重磁数据，并利用泰勒展开式方程来计算数据精度ε2；
S2、然后利用无人机计算航空重磁数据和下井测试井中重磁数据，并分别利用泰勒展开式方程来计算数据精度ε1和ε3；
S3、再对单一空间重磁数据进行地质体位置反演，并且根据不同的高度和深入重复步骤S2和步骤S3；
S4、根据不同数据的质量分配不同数据权重W，并根据多空间重磁数据按照相应公式进行反演；
S5、最后将单层反演窗口缩小而利用其它空间数据形成补充，重复上述操作根据不同数据的精度和根据精度设定相应的权重系数，进而组合为矩阵进行反演，最后得出地质体空间位置信息。
2.根据权利要求1所述的一种空地井立体重磁数据协同目标位置反演方法，其特征在于：所述步骤S1和S2中勘测重磁数据的泰勒展开式方程关系式为：
其中，T0(x，y，0)和T(x，y，h)表示高差为h的两个观测面上的异常，当h＞0时表示向上延拓，当h＜0时表示向下延拓，对公式进行Fourier变换可得到
其中，表示异常T0(x，y，0)的Fourier谱，(kx，ky)是分别表示x和y方向上的波数，称为延拓因子。
3.根据权利要求1所述的一种空地井立体重磁数据协同目标位置反演方法，其特征在于：所述步骤S3计算单一空间重磁数据进行地质体位置反演时的原理是重磁数据及其导数满足欧拉齐次方程：
(x-x0)·fx+(y-y0)·fy+(z-z0)·fz＝-N(f(x-x0，y-y0，z-z0)+A)
其中f为原始重磁异常，fx，fy，fz分别为异常f在不同方向的导数，x，y，z为观测点坐标，A为背景异，N为构造指数(人为给定)，x0，y0，z0为需反演的地质体位置坐标。
4.根据权利要求1所述的一种空地井立体重磁数据协同目标位置反演方法，其特征在于：所述步骤S4中对多空间重磁数据采用方程进行位置反演时的表达式为：
其中，h为代表观测高度。
5.根据权利要求1所述的一种空地井立体重磁数据协同目标位置反演方法，其特征在于：所述步骤S5中计算不同数据的精度，根据精度设定相应的权重系数，进而组合为矩阵进行反演的公式为：
CN201810185427.1A 2018-03-07 2018-03-07 一种空地井立体重磁数据协同目标位置反演方法 Expired - Fee Related CN108508479B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810185427.1A CN108508479B (zh) 2018-03-07 2018-03-07 一种空地井立体重磁数据协同目标位置反演方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810185427.1A CN108508479B (zh) 2018-03-07 2018-03-07 一种空地井立体重磁数据协同目标位置反演方法

Publications (2)

Publication Number Publication Date
CN108508479A true CN108508479A (zh) 2018-09-07
CN108508479B CN108508479B (zh) 2020-09-11

Family

ID=63376198

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810185427.1A Expired - Fee Related CN108508479B (zh) 2018-03-07 2018-03-07 一种空地井立体重磁数据协同目标位置反演方法

Country Status (1)

Country Link
CN (1) CN108508479B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112462442A (zh) * 2020-11-30 2021-03-09 山东大学 重磁位场场源位置估计方法、系统、介质及电子设备
CN114779352A (zh) * 2022-05-31 2022-07-22 成都理工大学 一种基于空-地-井多维电磁法的地质体勘测系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100610686B1 (ko) * 2006-02-27 2006-08-09 (주)인텔리지오 중력 또는 자력 탐사자료를 이용하여 탐사대상체의 2차원적위치를 구하는 방법.
US20130018588A1 (en) * 2011-07-11 2013-01-17 Technolmaging, Llc. Method of real time subsurface imaging using gravity and/or magnetic data measured from a moving platform
CN103258128A (zh) * 2013-05-08 2013-08-21 浙江大学 一种地磁场空间延拓算法的评估方法
CN104537714A (zh) * 2015-01-07 2015-04-22 吉林大学 磁共振与瞬变电磁空间约束联合反演方法
CN106291725A (zh) * 2015-05-13 2017-01-04 核工业北京地质研究院 一种快速反演地下地质体空间位置的方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100610686B1 (ko) * 2006-02-27 2006-08-09 (주)인텔리지오 중력 또는 자력 탐사자료를 이용하여 탐사대상체의 2차원적위치를 구하는 방법.
US20130018588A1 (en) * 2011-07-11 2013-01-17 Technolmaging, Llc. Method of real time subsurface imaging using gravity and/or magnetic data measured from a moving platform
CN103258128A (zh) * 2013-05-08 2013-08-21 浙江大学 一种地磁场空间延拓算法的评估方法
CN104537714A (zh) * 2015-01-07 2015-04-22 吉林大学 磁共振与瞬变电磁空间约束联合反演方法
CN106291725A (zh) * 2015-05-13 2017-01-04 核工业北京地质研究院 一种快速反演地下地质体空间位置的方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
ALAN B. REID,ER AL.: "The structural index in gravity and magnetic interpretation: Errors, uses, and abuses", 《GEOPHYSICS》 *
RAYMOND JOHN DURRHEIM,ET AL.: "EULDEP: a program for the Euler deconvolution of magnetic and gravity data", 《COMPUTERS & GEOSCIENCES 》 *
REZA TOUSHMALANI，ET AL.: "Euler deconvolution of 3D gravity data interpretation: New approach", 《JOURNAL OF APPLIED SCIENCE AND AGRICULTURE》 *
刘双 等，: "井-地-空磁数据协同反演研究", 《中国地球科学联合学术年会2016》 *
周文月: "重力及梯度数据联合欧拉反褶积方法研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *
李丽丽 等，: "归一化总水平导数法在位场数据解释中的应用", 《地球物理学报》 *
王明 等，: "欧拉反褶积在重磁位场中应用与发展", 《物探与化探》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112462442A (zh) * 2020-11-30 2021-03-09 山东大学 重磁位场场源位置估计方法、系统、介质及电子设备
CN112462442B (zh) * 2020-11-30 2022-04-08 山东大学 重磁位场场源位置估计方法、系统、介质及电子设备
CN114779352A (zh) * 2022-05-31 2022-07-22 成都理工大学 一种基于空-地-井多维电磁法的地质体勘测系统
CN114779352B (zh) * 2022-05-31 2022-11-22 成都理工大学 一种基于空-地-井多维电磁法的地质体勘测系统

Also Published As

Publication number Publication date
CN108508479B (zh) 2020-09-11

Similar Documents

Publication Publication Date Title
CN111927552B (zh) 基于天空地勘察技术的复杂山区长大深埋隧道勘察方法
Quigley et al. Previously unknown fault shakes New Zealand's South Island
CN106094029A (zh) 利用偏移距矢量片地震数据预测储层裂缝的方法
CN102927934A (zh) 一种利用单个InSAR干涉对获取矿区地表三维形变场的方法
CN102096944A (zh) 地质体结构面三维激光扫描点云识别方法
Valkaniotis et al. The Mw= 5.6 Kanallaki earthquake of 21 March 2020 in west Epirus, Greece: Reverse fault model from InSAR data and seismotectonic implications for Apulia-Eurasia collision
Nakatsuka et al. Reduction of magnetic anomaly observations from helicopter surveys at varying elevations
Bossennec et al. Structural architecture and permeability patterns of crystalline reservoir rocks in the Northern Upper Rhine Graben: Insights from surface analogues of the Odenwald
CN108508479A (zh) 一种空地井立体重磁数据协同目标位置反演方法
Zhang et al. Three-dimensional P-wave velocity structure of the Zhuxi ore deposit, South China revealed by control-source first-arrival tomography
Xiong et al. Lithofacies palaeogeography restoration of the Middle Devonian Guanwushan Formation in the Longmenshan area
Menichetti et al. Sentinel-1 Interferometry and UAV Aerial Survey for Mapping Coseismic Ruptures: Mts. Sibillini vs. Mt. Etna Volcano
Dueholm et al. Reservoir analog studies using multimodel photogrammetry: A new tool for the petroleum industry
Eisner et al. Effective VTI anisotropy for consistent monitoring of microseismic events
CN111045087A (zh) 一种隐伏伟晶状白岗岩型铀矿的勘查方法
Qi et al. Three-Dimensional Geological–Geophysical Modeling and Prospecting Indications of the Ashele Ore Concentration Area in Xinjiang Based on Irregular Sections
Guang-Ming Geology and petroleum potential of northwestern China
Wu et al. Application of double GPS multi-rotor UAV in the investigation of high slope perilous rock-mass in an open pit iron mine
Yang et al. Application of Induced Polarization Method in Mineral Resource Exploration
BIAN et al. Applying Two‐Step Iterative Least Square Approach to Determine the Geometry and Physical Parameters of Magnetic Sphere Sources
García-Sellés et al. Fracture analog of the sub-Andean Devonian of southern Bolivia: Lidar applied to Abra Del Condor
Wang et al. Comprehensive Geophysical Measurement in Seismic Safety Evaluation: A Case Study of East Lake High-Tech Development Zone, China
CN113536693B (zh) 一种基于井中岩石物性约束的航空-地面-井中磁异常数据联合反演方法
CN109581521A (zh) Tti各向异性的局部层析方法及系统
Bjornstad et al. Geothermal Resource Exploration, NAF El Centro–Superstition Mountain Area, Imperial Valley, California

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200911

CF01 Termination of patent right due to non-payment of annual fee