CN108492204B - 一种高精度空地井不同空间重磁数据变换方法 - Google Patents

一种高精度空地井不同空间重磁数据变换方法 Download PDF

Info

Publication number
CN108492204B
CN108492204B CN201810185429.0A CN201810185429A CN108492204B CN 108492204 B CN108492204 B CN 108492204B CN 201810185429 A CN201810185429 A CN 201810185429A CN 108492204 B CN108492204 B CN 108492204B
Authority
CN
China
Prior art keywords
anomaly
data
continuation
upward
abnormal
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
Application number
CN201810185429.0A
Other languages
English (en)
Other versions
CN108492204A (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 CN201810185429.0A priority Critical patent/CN108492204B/zh
Publication of CN108492204A publication Critical patent/CN108492204A/zh
Application granted granted Critical
Publication of CN108492204B publication Critical patent/CN108492204B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02—Agriculture; Fishing; Forestry; Mining

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Mining & Mineral Resources (AREA)
  • Agronomy & Crop Science (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Animal Husbandry (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Complex Calculations (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种高精度空地井不同空间重磁数据变换方法，所述高精度空地井不同空间重磁数据变换方法通过泰勒级数展开式证明了向下延拓计算可通过向上延拓和水平导数之和来完成，向上延拓和水平导数的计算是稳定的，受噪音影响较小，且为了弥补导数阶次为有效项造成的误差，采用一种异常向下延拓的水平导数迭代法，有效地增强了计算结果的稳定性和准确性，该发明通过泰勒级数展开式证明了向下延拓计算可通过向上延拓和水平导数之和来完成，向上延拓和水平导数的计算是稳定的，受噪音影响较小，且为了弥补导数阶次为有效项造成的误差，采用一种异常向下延拓的水平导数迭代法，有效地增强了计算结果的稳定性和准确性。

Description

一种高精度空地井不同空间重磁数据变换方法
技术领域
本发明属于地球科学技术领域，更具体地说，尤其涉及一种高精度空地井不同空间重磁数据变换方法。
背景技术
针对矿产、油气资源潜力区域，为了了解地下物质的不同深度上物质的分布，往往在一个地区开展航空、地面、井中重磁测量，从而利用不同空间的数据联合完成地下地质体的深度及物性反演工作，可有效地提高解释结果的精度和可靠性。但是在实际勘探中，不同空间数据的范围和比例尺是不一样的，从而大多时候需要利用其它空间的数据进行换算填补缺失的数据，所以不同空间重磁数据变换技术是进行空-地-井位场数据联合反演的前提，其变换精度也是制约反演精度的主要因素。
不同空间位场数据的获得是采用向上和向下延拓技术来实现，向上延拓采用频率域公式来进行，具有计算速度快、受噪声干扰下小的优点；向下延拓若采用频率域技术进行是不稳定的，且会明显增大噪声的干扰，使测得曲线发生强烈地跳动，从而造成测量精确度低，测得数据不准确的问题，越来越不能满足人们的需求，不利于广泛的推广和普及。
发明内容
本发明的目的是为了解决现有技术中存在的缺点，而提出的一种高精度空地井不同空间重磁数据变换方法。
为实现上述目的，本发明提供如下技术方案：一种高精度空地井不同空间重磁数据变换方法，包括如下步骤：
S1、首先在需要勘测地点的地片面计算再地表面的原始异常数据T，并同时勘测向上延伸的异常数据T-h；
S2、然后利用泰勒展开式方程来计算异常的二阶和四阶垂直导数，并且利用测得的导数来和泰勒展开式计算向下延伸的初始异常数据
Figure BSA0000160193750000021
S3、然后根据公式计算下界面向上延伸地平面后的异常数据T1，并把T1与原始异常数据T进行比较得出误差数据ΔT1；
S4、如果ΔT1均方差大于ε时，则重复步骤2和步骤3来重复计算过程，直到测得的ΔTm均方差小于ε时输出下延伸异常数据
Figure BSA0000160193750000022
优选的，所述步骤S1中观测面上重磁异常之间的关系式为：
Figure BSA0000160193750000023
其中，T0(x，y，0)和T(x，y，h)表示高差为h的两个观测面上的异常，当h＞0时表示向上延拓，当h＜0时表示向下延拓，对公式(1)进行Fourier变换可得到
Figure BSA0000160193750000024
其中，
Figure BSA0000160193750000025
表示异常T0(x，y，0)的Fourier谱，(kx，ky)是分别表示x和y方向上的波数，
Figure BSA0000160193750000026
弥为延拓因子。
优选的，所述步骤S2中泰勒展开式计算深度为h界面上的异常表达式为：
Figure BSA0000160193750000027
其中，T(x，y，h)是观测面h上的异常，h为延拓高度，T(x，y，0)是观测面上异常，向上延拓运算的泰勒展开式可表示为：
Figure BSA0000160193750000028
其中，T(x，y，-h)为观测面-h上的异常，将式(3)和(4)相加后可得：
Figure BSA0000160193750000029
Figure BSA0000160193750000031
优选的，所述步骤S2中对初始异常数据
Figure BSA0000160193750000032
进行技术的方程式为：
Figure BSA0000160193750000033
且采用Laplace方程来计算异常的二阶和四阶垂直导数：
Figure BSA0000160193750000034
优选的，所述步骤S3中计算T1与原始异常数据T的公式为：
ΔT1(x，y，0)＝T(x，y，0)-T1(x，y，0) (8)
利用公式(6)将ΔT1向下延拓高度h，可以得到：
Figure BSA0000160193750000035
其中，ΔT1(x，y，h)为异常ΔT1(x，y，0)下延h后异常，ΔT1(x，y，-h)为异常ΔT1(x，y，0)上延h后异常，ΔT1(x，y，h)作下界面异常的一个修正量，则下界面异常可改写为：
T1(x，y，h)＝T(x，y，h)+ΔT1(x，y，h) (10)。
优选的，所述步骤S4中计算过程，直至ΔTm(x，y，0)的均方差小于给定值ε，因此下界面h上的最终异常为：
Tm(x，y，h)＝T(x，y，h)+ΔT1(x，y，h)+…+ΔTm(x，y，h) (11)。
本发明的技术效果和优点：本发明提供的一种高精度空地井不同空间重磁数据变换方法，与传统的变换方法相比，本发明以通过泰勒级数展开式证明了向下延拓计算可通过向上延拓和水平导数之和来完成，向上延拓和水平导数的计算是稳定的，受噪音影响较小，且为了弥补导数阶次为有效项造成的误差，采用一种异常向下延拓的水平导数迭代法，有效地增强了计算结果的稳定性和准确性，该发明通过泰勒级数展开式证明了向下延拓计算可通过向上延拓和水平导数之和来完成，向上延拓和水平导数的计算是稳定的，受噪音影响较小，且为了弥补导数阶次为有效项造成的误差，采用一种异常向下延拓的水平导数迭代法，有效地增强了计算结果的稳定性和准确性。
附图说明
图1为本发明的向下延拓迭代过程工艺流程图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1所示的一种高精度空地井不同空间重磁数据变换方法，包括如下步骤：
S1、首先在需要勘测地点的地片面计算再地表面的原始异常数据T，并同时勘测向上延伸的异常数据T-h；
S2、然后利用泰勒展开式方程来计算异常的二阶和四阶垂直导数，并且利用测得的导数来和泰勒展开式计算向下延伸的初始异常数据
Figure BSA0000160193750000041
S3、然后根据公式计算下界面向上延伸地平面后的异常数据T1，并把T1与原始异常数据T进行比较得出误差数据ΔT1；
S4、如果ΔT1均方差大于ε时，则重复步骤2和步骤3来重复计算过程，直到测得的ΔTm均方差小于ε时输出下延伸异常数据
Figure BSA0000160193750000042
具体的，所述步骤S1中观测面上重磁异常之间的关系式为：
Figure BSA0000160193750000051
其中，T0(x，y，0)和T(x，y，h)表示高差为h的两个观测面上的异常，当h＞0时表示向上延拓，当h＜0时表示向下延拓，对公式(1)进行Fourier变换可得到
Figure BSA0000160193750000052
其中，
Figure BSA0000160193750000053
表示异常T0(x，y，0)的Fourier谱，(kx，ky)是分别表示x和y方向上的波数，
Figure BSA0000160193750000054
称为延拓因子。
具体的，所述步骤S2中泰勒展开式计算深度为h界面上的异常表达式为：
Figure BSA0000160193750000055
其中，T(x，y，h)是观测面h上的异常，h为延拓高度，T(x，y，0)是观测面上异常，向上延拓运算的泰勒展开式可表示为：
Figure BSA0000160193750000056
其中，T(x，y，-h)为观测面-h上的异常，将式(3)和(4)相加后可得：
Figure BSA0000160193750000057
具体的，所述步骤S2中对初始异常数据
Figure BSA0000160193750000058
进行技术的方程式为：
Figure BSA0000160193750000059
且采用Laplace方程来计算异常的二阶和四阶垂直导数：
Figure BSA00001601937500000510
具体的，所述步骤S3中计算T1与原始异常数据T的公式为：
ΔT1(x，y，0)＝T(x，y，0)-T1(x，y，0) (8)
利用公式(6)将ΔT1向下延拓高度h，可以得到：
Figure BSA0000160193750000061
其中，ΔT1(x，y，h)为异常ΔT1(x，y，0)下延h后异常，ΔT1(x，y，-h)为异常ΔT1(x，y，0)上延h后异常，ΔT1(x，y，h)作下界面异常的一个修正量，则下界面异常可改写为：
T1(x，y，h)＝T(x，y，h)+ΔT1(x，y，h) (10)。
具体的，所述步骤S4中计算过程，直至ΔTm(x，y，0)的均方差小于给定值ε，因此下界面h上的最终异常为：
Tm(x，y，h)＝T(x，y，h)+ΔT1(x，y，h)+…+ΔTm(x，y，h) (11)。
综上所述：本发明提供的一种高精度空地井不同空间重磁数据变换方法，与传统的变换方法相比，本发明以通过泰勒级数展开式证明了向下延拓计算可通过向上延拓和水平导数之和来完成，向上延拓和水平导数的计算是稳定的，受噪音影响较小，且为了弥补导数阶次为有效项造成的误差，采用一种异常向下延拓的水平导数迭代法，有效地增强了计算结果的稳定性和准确性，该发明通过泰勒级数展开式证明了向下延拓计算可通过向上延拓和水平导数之和来完成，向上延拓和水平导数的计算是稳定的，受噪音影响较小，且为了弥补导数阶次为有效项造成的误差，采用一种异常向下延拓的水平导数迭代法，有效地增强了计算结果的稳定性和准确性。
最后应说明的是：以上所述仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (1)

1.一种高精度空地井不同空间重磁数据变换方法，其特征在于,包括如下步骤：
S1、首先在需要勘测地点的地片面计算再地表面的原始异常数据T，并同时勘测向上延伸的异常数据T-h；
S2、然后利用泰勒展开式方程来计算异常的二阶和四阶垂直导数，并且利用测得的导数来和泰勒展开式计算向下延伸的初始异常数据Ti h；
S3、然后根据公式计算下界面向上延伸地平面后的异常数据T1，并把T1与原始异常数据T进行比较得出误差数据△T1；
S4、如果△T1均方差大于ε时，则重复步骤2和步骤3来重复计算过程，直到测得的△Tm均方差小于ε时输出下延伸异常数据Ti h；
所述步骤S1中观测面上重磁异常之间的关系式为：
Figure FDA0003017310650000011
其中，T0(x,y,0)和T(x,y,h)表示高差为h的两个观测面上的异常，当h>0时表示向上延拓，当h<0时表示向下延拓，对公式(1)进行Fourier变换可得到：
Figure FDA0003017310650000012
其中，
Figure FDA0003017310650000013
表示异常T0(x,y,0)的Fourier谱，kx,ky分别表示x和y方向上的波数，
Figure FDA0003017310650000014
称为延拓因子；
所述步骤S2中泰勒展开式计算深度为h界面上的异常表达式为：
Figure FDA0003017310650000021
其中，T(x,y,h)是观测面h上的异常，h为延拓高度，T(x,y,0)是观测面上异常，向上延拓运算的泰勒展开式可表示为：
Figure FDA0003017310650000022
其中，T(x,y,-h)为观测面-h上的异常，将式(3)和(4)相加后可得：
Figure FDA0003017310650000023
所述步骤S2中对初始异常数据Ti h进行技术的方程式为：
Figure FDA0003017310650000024
且采用Laplace方程来计算异常的二阶和四阶垂直导数：
Figure FDA0003017310650000025
所述步骤S3中计算T1与原始异常数据T的公式为：
ΔT1(x，y，0)＝T(x，y，0)-T1(x，y，0) (8)；
利用公式(6)将△T1向下延拓高度h，可以得到：
Figure FDA0003017310650000026
其中，ΔT1(x,y,h)为异常ΔT1(x,y,0)下延h后异常，ΔT1(x,y,-h)为异常ΔT1(x,y,0)上延h后异常，ΔT1(x,y,h)作下界面异常的一个修正量，则下界面异常可改写为：
T1(x，y，h)＝T(x，y，h)+ΔT1(x，y，h) (10)；
所述步骤S4中计算过程，直至ΔTm(x，y，0)的均方差小于给定值ε，因此下界面h上的最终异常为：
Tm(x，y，h)＝T(x，y，h)+ΔT1(x，y，h)+…+ΔTm(x，y，h) (11)。
CN201810185429.0A 2018-03-07 2018-03-07 一种高精度空地井不同空间重磁数据变换方法 Expired - Fee Related CN108492204B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810185429.0A CN108492204B (zh) 2018-03-07 2018-03-07 一种高精度空地井不同空间重磁数据变换方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810185429.0A CN108492204B (zh) 2018-03-07 2018-03-07 一种高精度空地井不同空间重磁数据变换方法

Publications (2)

Publication Number Publication Date
CN108492204A CN108492204A (zh) 2018-09-04
CN108492204B true CN108492204B (zh) 2021-06-01

Family

ID=63341759

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810185429.0A Expired - Fee Related CN108492204B (zh) 2018-03-07 2018-03-07 一种高精度空地井不同空间重磁数据变换方法

Country Status (1)

Country Link
CN (1) CN108492204B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110133749B (zh) * 2019-05-30 2021-06-22 中国地质科学院 一种地质资源勘探中重磁数据处理方法及其系统
CN114167511B (zh) * 2021-11-26 2023-08-11 兰州大学 一种基于连分式展开向下延拓的位场数据快速反演方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100610686B1 (ko) * 2006-02-27 2006-08-09 (주)인텔리지오 중력 또는 자력 탐사자료를 이용하여 탐사대상체의 2차원적위치를 구하는 방법.
KR20080023946A (ko) * 2006-09-12 2008-03-17 한국지질자원연구원 오일러 해를 이용한 지하공동의 3차원 중력역산 방법 및이를 이용한 3차원 영상화 방법
CN106291725A (zh) * 2015-05-13 2017-01-04 核工业北京地质研究院 一种快速反演地下地质体空间位置的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100610686B1 (ko) * 2006-02-27 2006-08-09 (주)인텔리지오 중력 또는 자력 탐사자료를 이용하여 탐사대상체의 2차원적위치를 구하는 방법.
KR20080023946A (ko) * 2006-09-12 2008-03-17 한국지질자원연구원 오일러 해를 이용한 지하공동의 3차원 중력역산 방법 및이를 이용한 3차원 영상화 방법
CN106291725A (zh) * 2015-05-13 2017-01-04 核工业北京地质研究院 一种快速反演地下地质体空间位置的方法

Also Published As

Publication number Publication date
CN108492204A (zh) 2018-09-04

Similar Documents

Publication Publication Date Title
CN114966856B (zh) 基于多频带地震资料的碳封存场址优选方法、系统和设备
US6885943B2 (en) Simultaneous resolution enhancement and dip correction of resistivity logs through nonlinear iterative deconvolution
CN107783187B (zh) 一种将测井速度和地震速度结合建立三维速度场的方法
Araffa et al. Implementation of magnetic and gravity methods to delineate the subsurface structural features of the basement complex in central Sinai area, Egypt
CN110954943B (zh) 被动源地震频率谐振勘探方法
Zhang et al. An improved and stable downward continuation of potential field data: The truncated Taylor series iterative downward continuation method
CN108492204B (zh) 一种高精度空地井不同空间重磁数据变换方法
CN112114374B (zh) 复杂地质体的多密度界面反演方法
WO2017135183A1 (ja) 重力偏差データの表層密度値推定方法
Zeng et al. Recursive impedance inversion of ground-penetrating radar data in stochastic media
CN108319566B (zh) 基于向上延拓的航空重力点对点向下延拓解析方法
CN110737018B (zh) Vsp地震资料各向异性建模方法
CN108957554A (zh) 一种地球物理勘探中地震反演方法
CN109490978B (zh) 一种起伏地层的频率域快速高精度正演方法
CN106970426B (zh) 一种基于测线差分与主成分分析的航空电磁数据调平方法
CN112241025B (zh) 一种井震联合地层压力确定方法及系统
CN115577616A (zh) 基于深度学习的碳酸岩缝洞地震刻画方法及装置
CN108508479A (zh) 一种空地井立体重磁数据协同目标位置反演方法
CN104007465A (zh) 基于弹性反演的纯纵波数据刻画储层发育细节方法
CN109901221B (zh) 一种基于动校正速度参数的地震资料各向异性建模方法
CN109061737B (zh) 一种基于合成地震记录的储层预测方法及装置
Wu et al. Study on the relationship between multi-stage strike-slip mechanism and basin evolution in Fangzheng fault depression
CN111142158B (zh) 一种提高断阶构造地震解释精度的标定方法
CN108549103A (zh) 饱和多孔介质叠后地震烃类检测方法
Guo et al. A wavenumber-domain iterative approach for apparent density mapping of an undulant layer and its application in central South China

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: 20210601

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