CN108492204A - 一种高精度空地井不同空间重磁数据变换方法 - Google Patents
一种高精度空地井不同空间重磁数据变换方法 Download PDFInfo
- Publication number
- CN108492204A CN108492204A CN201810185429.0A CN201810185429A CN108492204A CN 108492204 A CN108492204 A CN 108492204A CN 201810185429 A CN201810185429 A CN 201810185429A CN 108492204 A CN108492204 A CN 108492204A
- Authority
- CN
- China
- Prior art keywords
- abnormal
- gravity
- precision
- continuation
- different spaces
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000005484 gravity Effects 0.000 title claims abstract description 16
- 230000002159 abnormal effect Effects 0.000 claims abstract description 43
- 238000007689 inspection Methods 0.000 claims description 15
- 230000008569 process Effects 0.000 claims description 7
- 238000005516 engineering process Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 3
- 230000002035 prolonged effect Effects 0.000 claims 2
- 230000000694 effects Effects 0.000 abstract description 6
- 230000008901 benefit Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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、然后利用泰勒展开式方程来计算异常的二阶和四阶垂直导数,并且利用测得的导数来和泰勒展开式计算向下延伸的初始异常数据
S3、然后根据公式计算下界面向上延伸地平面后的异常数据T1,并把T1与原始异常数据T进行比较得出误差数据ΔT1;
S4、如果ΔT1均方差大于ε时,则重复步骤2和步骤3来重复计算过程,直到测得的ΔTm均方差小于ε时输出下延伸异常数据
优选的,所述步骤S1中观测面上重磁异常之间的关系式为:
其中,T0(x,y,0)和T(x,y,h)表示高差为h的两个观测面上的异常,当h>0时表示向上延拓,当h<0时表示向下延拓,对公式(1)进行Fourier变换可得到
其中,表示异常T0(x,y,0)的Fourier谱,(kx,ky)是分别表示x和y方向上的波数,弥为延拓因子。
优选的,所述步骤S2中泰勒展开式计算深度为h界面上的异常表达式为:
其中,T(x,y,h)是观测面h上的异常,h为延拓高度,T(x,y,0)是观测面上异常,向上延拓运算的泰勒展开式可表示为:
其中,T(x,y,-h)为观测面-h上的异常,将式(3)和(4)相加后可得:
优选的,所述步骤S2中对初始异常数据进行技术的方程式为:
且采用Laplace方程来计算异常的二阶和四阶垂直导数:
优选的,所述步骤S3中计算T1与原始异常数据T的公式为:
ΔT1(x,y,0)=T(x,y,0)-T1(x,y,0) (8)
利用公式(6)将ΔT1向下延拓高度h,可以得到:
其中,Δ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、然后利用泰勒展开式方程来计算异常的二阶和四阶垂直导数,并且利用测得的导数来和泰勒展开式计算向下延伸的初始异常数据
S3、然后根据公式计算下界面向上延伸地平面后的异常数据T1,并把T1与原始异常数据T进行比较得出误差数据ΔT1;
S4、如果ΔT1均方差大于ε时,则重复步骤2和步骤3来重复计算过程,直到测得的ΔTm均方差小于ε时输出下延伸异常数据
具体的,所述步骤S1中观测面上重磁异常之间的关系式为:
其中,T0(x,y,0)和T(x,y,h)表示高差为h的两个观测面上的异常,当h>0时表示向上延拓,当h<0时表示向下延拓,对公式(1)进行Fourier变换可得到
其中,表示异常T0(x,y,0)的Fourier谱,(kx,ky)是分别表示x和y方向上的波数,称为延拓因子。
具体的,所述步骤S2中泰勒展开式计算深度为h界面上的异常表达式为:
其中,T(x,y,h)是观测面h上的异常,h为延拓高度,T(x,y,0)是观测面上异常,向上延拓运算的泰勒展开式可表示为:
其中,T(x,y,-h)为观测面-h上的异常,将式(3)和(4)相加后可得:
具体的,所述步骤S2中对初始异常数据进行技术的方程式为:
且采用Laplace方程来计算异常的二阶和四阶垂直导数:
具体的,所述步骤S3中计算T1与原始异常数据T的公式为:
ΔT1(x,y,0)=T(x,y,0)-T1(x,y,0) (8)
利用公式(6)将ΔT1向下延拓高度h,可以得到:
其中,Δ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 (6)
1.一种高精度空地井不同空间重磁数据变换方法,其特征在于,包括如下步骤:
S1、首先在需要勘测地点的地片面计算再地表面的原始异常数据T,并同时勘测向上延伸的异常数据T-h;
S2、然后利用泰勒展开式方程来计算异常的二阶和四阶垂直导数,并且利用测得的导数来和泰勒展开式计算向下延伸的初始异常数据Ti h;
S3、然后根据公式计算下界面向上延伸地平面后的异常数据T1,并把T1与原始异常数据T进行比较得出误差数据ΔT1;
S4、如果ΔT1均方差大于ε时,则重复步骤2和步骤3来重复计算过程,直到测得的ΔTm均方差小于ε时输出下延伸异常数据Ti h。
2.根据权利要求1所述的一种高精度空地井不同空间重磁数据变换方法,其特征在于:所述步骤S1中观测面上重磁异常之间的关系式为:
其中,T0(x,y,0)和T(x,y,h)表示高差为h的两个观测面上的异常,当h>0时表示向上延拓,当h<0时表示向下延拓,对公式(1)进行Fourier变换可得到
其中,表示异常T0(x,y,0)的Fourier谱,(kx,ky)是分别表示x和y方向上的波数,称为延拓因子。
3.根据权利要求1所述的一种高精度空地井不同空间重磁数据变换方法,其特征在于:所述步骤S2中泰勒展开式计算深度为h界面上的异常表达式为:
其中,T(x,y,h)是观测面h上的异常,h为延拓高度,T(x,y,0)是观测面上异常,向上延拓运算的泰勒展开式可表示为:
其中,T(x,y,-h)为观测面-h上的异常,将式(3)和(4)相加后可得:
4.根据权利要求1所述的一种高精度空地井不同空间重磁数据变换方法,其特征在于:所述步骤S2中对初始异常数据Ti h进行技术的方程式为:
且采用Laplace方程来计算异常的二阶和四阶垂直导数:
5.根据权利要求1所述的一种高精度空地井不同空间重磁数据变换方法,其特征在于:所述步骤S3中计算T1与原始异常数据T的公式为:
ΔT1(x,y,0)=T(x,y,0)-T1(x,y,0) (8)
利用公式(6)将ΔT1向下延拓高度h,可以得到:
其中,Δ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)。
6.根据权利要求1所述的一种高精度空地井不同空间重磁数据变换方法,其特征在于:所述步骤S4中计算过程,直至ΔTm(x,y,0)的均方差小于给定值ε,因此下界面h上的最终异常为:
Tm(x,y,h)=T(x,y,h)+ΔT1(x,y,h)+…+ΔTm(x,y,h) (11)。
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 true CN108492204A (zh) | 2018-09-04 |
CN108492204B 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) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110133749A (zh) * | 2019-05-30 | 2019-08-16 | 中国地质科学院 | 一种地质资源勘探中重磁数据处理方法及其系统 |
CN114167511A (zh) * | 2021-11-26 | 2022-03-11 | 兰州大学 | 一种基于连分式展开向下延拓的位场数据快速反演方法 |
Citations (3)
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 | 核工业北京地质研究院 | 一种快速反演地下地质体空间位置的方法 |
-
2018
- 2018-03-07 CN CN201810185429.0A patent/CN108492204B/zh not_active Expired - Fee Related
Patent Citations (3)
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 | 核工业北京地质研究院 | 一种快速反演地下地质体空间位置的方法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110133749A (zh) * | 2019-05-30 | 2019-08-16 | 中国地质科学院 | 一种地质资源勘探中重磁数据处理方法及其系统 |
CN114167511A (zh) * | 2021-11-26 | 2022-03-11 | 兰州大学 | 一种基于连分式展开向下延拓的位场数据快速反演方法 |
CN114167511B (zh) * | 2021-11-26 | 2023-08-11 | 兰州大学 | 一种基于连分式展开向下延拓的位场数据快速反演方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108492204B (zh) | 2021-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Witte et al. | Fracture evolution in oil-producing sills of the Rio Grande Valley, northern Neuquén Basin, Argentina | |
CN102636812B (zh) | 一种获得碳酸盐岩储层储集空间体积的方法 | |
CN105549082B (zh) | 超深碳酸盐储层三维地质力学场建立方法和系统 | |
CN109238229A (zh) | 一种基于静力水准仪监测的地表沉降值的温度补偿方法 | |
CN104537714A (zh) | 磁共振与瞬变电磁空间约束联合反演方法 | |
CN105223617A (zh) | 利用地震属性控制初始波阻抗的建模方法及系统 | |
CN105137482A (zh) | 一种沉积体古坡度的计算方法 | |
CN109799540B (zh) | 基于地质信息约束下的火山岩型铀矿床磁化率反演方法 | |
CN109458176A (zh) | 碳酸盐岩储层压力的预测方法及其应用 | |
CN108492204A (zh) | 一种高精度空地井不同空间重磁数据变换方法 | |
CN107784159A (zh) | 一种储层电阻率各向异性系数的确定方法 | |
CN109839670A (zh) | 一种热液型铀矿床基底界面反演方法 | |
CN103541718A (zh) | 用图版解释油井含水率的方法 | |
CN109490978A (zh) | 一种起伏地层的频率域快速高精度正演方法 | |
CN112114374A (zh) | 复杂地质体的多密度界面反演方法 | |
Zahorec et al. | The role of near topography and building effects in vertical gravity gradients approximation | |
Elmasry et al. | Exploration of geothermal resources utilizing geophysical and borehole data in the Abu Gharadig basin of Egypt’s Northern Western desert | |
Okubo et al. | Thermal and crustal structure of the Aso volcano and surrounding regions constrained by gravity and magnetic data, Japan | |
Khutorskoi et al. | Thermal tomography of the West Arctic basin | |
Brady et al. | Surface-gravity monitoring of the gas cap water injection project, Prudhoe Bay, Alaska | |
CN101598806A (zh) | 一种提高构造图准确度的偏差消除方法 | |
Atkinson et al. | Using precision gravity data in geothermal reservoir engineering modeling studies | |
CN104007436B (zh) | 基于机载LiDAR数据的公路勘测设计方法 | |
Meng et al. | Numerical Modeling and Stability Analysis of Surrounding Rock of Yuanjue Cave | |
Kawasaki et al. | Estimation of regional stress and stress distribution |
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 |