CN108508479B - 一种空地井立体重磁数据协同目标位置反演方法 - Google Patents
一种空地井立体重磁数据协同目标位置反演方法 Download PDFInfo
- Publication number
- CN108508479B CN108508479B CN201810185427.1A CN201810185427A CN108508479B CN 108508479 B CN108508479 B CN 108508479B CN 201810185427 A CN201810185427 A CN 201810185427A CN 108508479 B CN108508479 B CN 108508479B
- Authority
- CN
- China
- Prior art keywords
- gravity
- inversion
- data
- magnetic data
- anomaly
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 20
- 230000005484 gravity Effects 0.000 claims abstract description 39
- 239000000126 substance Substances 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 239000002356 single layer Substances 0.000 claims description 6
- 230000002159 abnormal effect Effects 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims 3
- 238000005070 sampling Methods 0.000 claims 1
- 230000001502 supplementing effect Effects 0.000 claims 1
- 230000009467 reduction Effects 0.000 abstract description 10
- 230000004044 response Effects 0.000 abstract description 5
- 230000009286 beneficial effect Effects 0.000 abstract description 3
- 239000013589 supplement Substances 0.000 description 4
- 230000008859 change Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Images
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变换可得到
优选的,所述步骤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变换可得到
优选的,所述步骤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和地质体位置反演;
S4、根据不同数据的精度分配不同数据权重W,并根据多空间重磁数据按照相应公式进行反演;
S5、最后将单层反演窗口缩小而利用其它空间数据形成补充,重复上述操作,根据不同数据的精度和根据精度设定相应的权重系数,进而组合为矩阵进行反演,最后得出地质体空间位置信息。
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为需反演的地质体位置坐标。
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 CN108508479A (zh) | 2018-09-07 |
CN108508479B true 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) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112462442B (zh) * | 2020-11-30 | 2022-04-08 | 山东大学 | 重磁位场场源位置估计方法、系统、介质及电子设备 |
CN114779352B (zh) * | 2022-05-31 | 2022-11-22 | 成都理工大学 | 一种基于空-地-井多维电磁法的地质体勘测系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103258128A (zh) * | 2013-05-08 | 2013-08-21 | 浙江大学 | 一种地磁场空间延拓算法的评估方法 |
CN104537714A (zh) * | 2015-01-07 | 2015-04-22 | 吉林大学 | 磁共振与瞬变电磁空间约束联合反演方法 |
Family Cites Families (3)
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 |
CN106291725B (zh) * | 2015-05-13 | 2018-11-30 | 核工业北京地质研究院 | 一种快速反演地下地质体空间位置的方法 |
-
2018
- 2018-03-07 CN CN201810185427.1A patent/CN108508479B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103258128A (zh) * | 2013-05-08 | 2013-08-21 | 浙江大学 | 一种地磁场空间延拓算法的评估方法 |
CN104537714A (zh) * | 2015-01-07 | 2015-04-22 | 吉林大学 | 磁共振与瞬变电磁空间约束联合反演方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108508479A (zh) | 2018-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111927552B (zh) | 基于天空地勘察技术的复杂山区长大深埋隧道勘察方法 | |
CN102927934B (zh) | 一种利用单个InSAR干涉对获取矿区地表三维形变场的方法 | |
Barnes et al. | Geological analysis of Martian rover‐derived digital outcrop models using the 3‐D visualization tool, Planetary Robotics 3‐D Viewer—Pro3D | |
Li et al. | An efficient and adaptive approach for modeling gravity effects in spherical coordinates | |
CN105938193B (zh) | 一种无需地面辅助的升降轨InSAR监测沉降区绝对地表形变方法 | |
CN110007357B (zh) | 一种航空tem和航空mt联合反演方法 | |
CN110045432A (zh) | 基于3d-glq的球坐标系下重力场正演方法及三维反演方法 | |
CN110187382B (zh) | 一种回折波和反射波波动方程旅行时反演方法 | |
Luo et al. | Full magnetic gradient tensor from triaxial aeromagnetic gradient measurements: Calculation and application | |
Vinnik et al. | Joint inversion of P-and S-receiver functions and dispersion curves of Rayleigh waves: the results for the Central Anatolian Plateau | |
Shockley | Ground vehicle navigation using magnetic field variation | |
CN106483559A (zh) | 一种地下速度模型的构建方法 | |
CN110361788B (zh) | 一种空-地联合三维重力数据特征分析及密度反演方法 | |
GB2481643A (en) | Gravity survey data processing | |
CN108508479B (zh) | 一种空地井立体重磁数据协同目标位置反演方法 | |
CN105629299A (zh) | 角度域叠前深度偏移的走时、角度表获取方法及成像方法 | |
CN102636809A (zh) | 一种传播角度域共成像点道集的生成方法 | |
Guo et al. | A hybrid positive-and-negative curvature approach for detection of the edges of magnetic anomalies, and its application in the South China Sea | |
CN113960532A (zh) | 一种基于假想源的二次定位计算的微地震定位方法 | |
CN113900069A (zh) | 一种基于干涉成像高度计的垂线偏差计算方法及其系统 | |
CN113552637A (zh) | 一种航空-地面-井中磁异常数据协同三维反演方法 | |
CN112596113A (zh) | 一种基于重力不同阶梯度特征值交点的场源位置识别方法 | |
CN106338760B (zh) | 误差补偿的起伏地表偏移方法 | |
CN111983668B (zh) | 一种获取地震参数估计的方法及系统 | |
CN108375794A (zh) | 基于对称观测的vsp缝洞绕射成像技术方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200911 |