CN105549099A - 基于全空间正则化下延数据的视磁化强度三维反演方法 - Google Patents
基于全空间正则化下延数据的视磁化强度三维反演方法 Download PDFInfo
- Publication number
- CN105549099A CN105549099A CN201510931800.XA CN201510931800A CN105549099A CN 105549099 A CN105549099 A CN 105549099A CN 201510931800 A CN201510931800 A CN 201510931800A CN 105549099 A CN105549099 A CN 105549099A
- Authority
- CN
- China
- Prior art keywords
- value
- sigma
- inversion
- regularization
- alpha
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
- G01V3/26—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device
- G01V3/28—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device using induction coils
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Electromagnetism (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
技术领域
本发明属于数学三维模型技术领域,涉及基于全空间正则化下延数据的视磁化强度三维反演方法。
背景技术
传统的视磁化强度反演方法,是利用在观测面采集的实际磁力数据来反推地下异常体的磁性参数以及空间赋存状态的一种地球物理手段,这种方法的基本思想就是对地下待反演地质体进行网格剖分,通过不停的反演迭代直至达到精度要求,来计算出所有剖分网格单元的视磁化强度。虽然很容易找到满足数据拟合空间的解模型,但是由于重磁场的体积效应、观测数据有限性及反演问题的欠定性等几种因素的存在,反演结果往往产生多解性。为了减小这种多种性的影响,提高反演的精度,前人做了大量的研究并提出了一系列方法,如已知井资料定量约束、概率成像定性约束等等。这些方法虽然一定程度上克服了多解性的影响,提高了反演结果的精度,但并没有解决上文提到的导致反演多解性的根本问题——数据的有限性。这种有限性主要体现在两个方面:(1)磁化强度是地下不同深度场的累加效应,而在数据采集时,只能测量到在地表观测到的数据;(2)当场源埋深较大时,地表观测的异常信号较弱,不能有效突出异常信号。
发明内容
本发明的目的在于提供基于全空间正则化下延数据的视磁化强度三维反演方法,解决了目前三维反演过程中由于数据有限,导致反演精度低的问题。
本发明所采用的技术方案是按照以下步骤进行:
步骤1:采用三维正则化下延算子实现不同深度的延拓得到稳定的磁异常全空间场,采用的三维正则化算子为:
步骤2:求解目标方程的极小值,得到最小值对应的β、α、f1、f2,代入正则化下延算子从而计算出不同深度的正延场值Eh(x,y,z)=E0(x,y,z)*hα,E0为地表观测值,hα为Hα的反傅氏变换;
步骤3:构建全空间场反演目标函数;
选取反演目标函数为 式中,为全空间场值Ek(Ek=E0+Eh)的平均值,该观测场包含地表观测值E0和步骤1下延得到的不同深度的场值Eh,fk表示第k个观测点的观测值和正演计算值之间的差;为fk的平均值,其数学表达式如下:
步骤4:采用多维搜索黄金分割算法求解反演目标函数的极小值,提高反演速度。
进一步,所述步骤1中β取值在0.5-2.5之间,以0.1为单位。
本发明的有益效果是实现全深度稳定下延,能够获得精度较为可靠的全空间场。
附图说明
图1是本发明反演方法流程图;
图2展示了单个异常体模型的反演结果;
图2(a)模型地表磁异常正演等值线图;
图2(b)下延不同深度的磁异常分布图;
图2(c)常规的仅用地表数据反演结果;
图2(d)全空间数据反演结果;
图3展示了两个不等深地质体模型的反演结果;
图3(a)模型地表磁异常正演等值线图;
图3(b)下延不同深度的磁异常分布图;
图3(c)常规的仅用地表数据反演结果;
图3(d)全空间数据反演结果。
具体实施方式
下面结合具体实施方式对本发明进行详细说明。
如图1所示,本发明提出了一种新的视磁化强度反演方法——基于正则化下延方法的视磁化强度全空间反演。该方法的核心思想是通过将观测磁异常场向下延拓,得到某一深度内的磁异常全空间场,进而增大反演数据量;另外,在向下延拓的过程中,磁异常场不断向场源接近,可以有效突出异常信号。
得到精度较高的全空间场是实现本发明关键所在,是实现本发明方法的前提。然而许多位场下延方法在延拓至场源体顶面深度附近区域时,会因为高频振荡效应使场值发生畸变,这样获得的全空间场显然在场源附近是不准确的。本发明采用一种基于深度变化的磁场正则化下延方法,利用正则化理论对下延算子中的各参数进行了求取,经实践证明,该方法可以实现全深度稳定下延,能够获得精度较为可靠的全空间场。
步骤1:采用三维正则化下延算子实现不同深度的延拓得到稳定的磁异常全空间场,采用的三维正则化算子为 其中,α是正则化参数,β是滤波系数,f1、f2为通带范围,f为波数,z为深度,λx、λy为基波波长,u,v为二维频域坐标系。构建目标方程进行自动求取最佳的α、β、f1、f2值,目标方程为:
步骤2:由前人经验β取值在0.5-2.5之间效果最好,以0.1为单位,在0.5-2.5之间分别赋予β,不同的β时,求解目标方程的极小值,最小值对应的α、f1、f2即为所求。将得到的α、β、f1、f2代入正则化下延算子从而计算出不同深度的正延场值Eh(x,y,z)=E0(x,y,z)*hα,E0为地表观测值,hα为Hα的反傅氏变换。
步骤3:通过步骤1中得到的结果,构建合理的全空间场反演目标函数,大大提高了反演数据量,增加了先验信息的约束。
选取反演目标函数为 式中,为全空间场值Ek(Ek=E0+Eh)的平均值,该观测场包含地表观测值E0和步骤1下延得到的不同深度的场值Eh。fk表示第k个观测点的观测值和正演计算值之间的差;为fk的平均值,其数学表达式如下: 该目标函数的优点是模型本身不受正常场的影响,即当实际数据与正演计算值之间存在较大差时,目标函数同样可以取到最小。
步骤4:采用改进的多维搜索黄金分割算法求解反演目标函数的极小值,提高反演速度。
图2展示了单个异常体模型的反演结果,模型中异常体磁化强度为600×10-3A/m,磁倾角90°,磁偏角为0,模型为200m×200m×200m规则正方体,在X方向700~900m,在Y方向700~900m,在Z方向600~800m;由反演结果可以看出,图2(c)常规仅用地表数据反演出的地质体磁化强度最高值为530A/m,反演出的目标体埋度大于600,且反演结果在纵向上相对于实际位置被明显拉长,在横向上与实际位置基本吻合;图2(d)利用本文提出的全空间数据进行反演,得到的地质体磁化强度最高值为628A/m,反演目标体的埋深在600~850m之间,在横向上与实际位置基本吻合,但在纵向上其结果要明显好于常规的仅用地表数据反演的方法(图中的黑框为异常体实际位置);
图3展示了两个不等深地质体模型的反演结果,模型中异常体磁化强度为500×10-3A/m,磁倾角90°,磁偏角为0,两个异常体均为200m×200m×200m规则正方体,具体的位置分别是异常体1在X方向300~500m,在Y方向700~900m,在Z方向200~400m;异常体2在X方向1100~1300m,在Y方向700~900m,在Z方向300~500m。由反演结果可以看出,图3(c)仅利用地表观测值反演出的结果是两个大小不一,但磁化强度相近的两个异常体,而且还有多余构造。图3(d)是本文提出的利用全空间数据进行反演的结果,其反演结果是两个大小相似、磁化强度相近、埋深不同的两个异常体,与模型实际情况的吻合度要明显高于常规仅用地表观测值反演的方法(图中的黑框为异常体实际位置)。因此,本文提出的方法可以在一定程度上克服反演多解性的影响,可以更好地反应地下异常体的真实形态。
以上所述仅是对本发明的较佳实施方式而已,并非对本发明作任何形式上的限制,凡是依据本发明的技术实质对以上实施方式所做的任何简单修改,等同变化与修饰,均属于本发明技术方案的范围内。
Claims (2)
1.基于全空间正则化下延数据的视磁化强度三维反演方法,其特征在于按照以下步骤进行反演:
步骤1:采用三维正则化下延算子实现不同深度的延拓得到稳定的磁异常全空间场,采用的三维正则化算子为:
为偏差泛函,为包含误差的混合值,Smn为不含误差的有效值;
步骤2:求解目标方程的极小值,得到最小值对应的β、α、f1、f2,代入正则化下延算子从而计算出不同深度的正延场值Eh(x,y,z)=E0(x,y,z)*hα,E0为地表观测值,hα为Hα的反傅氏变换;
步骤3:构建全空间场反演目标函数;
选取反演目标函数为 式中,为全空间场值Ek(Ek=E0+Eh)的平均值,该观测场包含地表观测值E0和步骤1下延得到的不同深度的场值Eh,fk表示第k个观测点的观测值和正演计算值之间的差;为fk的平均值,其数学表达式如下:
步骤4:采用多维搜索黄金分割算法求解反演目标函数的极小值,提高反演速度。
2.按照权利要求1所述基于全空间正则化下延数据的视磁化强度三维反演方法,其特征在于:所述步骤1中β取值在0.5-2.5之间,以0.1为单位。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510931800.XA CN105549099A (zh) | 2015-12-11 | 2015-12-11 | 基于全空间正则化下延数据的视磁化强度三维反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510931800.XA CN105549099A (zh) | 2015-12-11 | 2015-12-11 | 基于全空间正则化下延数据的视磁化强度三维反演方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105549099A true CN105549099A (zh) | 2016-05-04 |
Family
ID=55828398
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510931800.XA Pending CN105549099A (zh) | 2015-12-11 | 2015-12-11 | 基于全空间正则化下延数据的视磁化强度三维反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105549099A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106023085A (zh) * | 2016-06-12 | 2016-10-12 | 西南科技大学 | 一种稀疏约束sar图像重建正则化参数的sure黄金分割自动搜索算法 |
CN106199742A (zh) * | 2016-06-29 | 2016-12-07 | 吉林大学 | 一种频率域航空电磁法2.5维带地形反演方法 |
CN106855642A (zh) * | 2016-11-21 | 2017-06-16 | 中国石油天然气股份有限公司 | 确定基底火山岩分布的方法和装置 |
CN110989002A (zh) * | 2019-11-07 | 2020-04-10 | 吉林大学 | 地空时域电磁系统浅部低阻异常体高精度数据解释方法 |
CN113484919A (zh) * | 2021-06-30 | 2021-10-08 | 长江大学 | 基于卷积神经网络的磁异常反演方法、系统、终端及介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101285896A (zh) * | 2008-06-13 | 2008-10-15 | 杨辉 | 一种地球物理勘探中的重磁数据处理方法 |
CN101661115A (zh) * | 2008-08-29 | 2010-03-03 | 中国石油集团东方地球物理勘探有限责任公司 | 基于标准格架的快速三维重力、磁力物性反演的方法 |
CN103399350A (zh) * | 2013-07-29 | 2013-11-20 | 中国人民解放军国防科学技术大学 | 一种基于积分迭代算法的航空重力向下延拓方法 |
CN103869376A (zh) * | 2014-03-20 | 2014-06-18 | 中国石油大学(华东) | 一种基于深度变化的三维重力位场正则化下延方法及应用 |
CN104866653A (zh) * | 2015-04-29 | 2015-08-26 | 中国地质科学院矿产资源研究所 | 一种获取地下三维密度结构的方法 |
-
2015
- 2015-12-11 CN CN201510931800.XA patent/CN105549099A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101285896A (zh) * | 2008-06-13 | 2008-10-15 | 杨辉 | 一种地球物理勘探中的重磁数据处理方法 |
CN101661115A (zh) * | 2008-08-29 | 2010-03-03 | 中国石油集团东方地球物理勘探有限责任公司 | 基于标准格架的快速三维重力、磁力物性反演的方法 |
CN103399350A (zh) * | 2013-07-29 | 2013-11-20 | 中国人民解放军国防科学技术大学 | 一种基于积分迭代算法的航空重力向下延拓方法 |
CN103869376A (zh) * | 2014-03-20 | 2014-06-18 | 中国石油大学(华东) | 一种基于深度变化的三维重力位场正则化下延方法及应用 |
CN104866653A (zh) * | 2015-04-29 | 2015-08-26 | 中国地质科学院矿产资源研究所 | 一种获取地下三维密度结构的方法 |
Non-Patent Citations (4)
Title |
---|
刘展 等: ""济阳坳陷花沟地区火成岩重磁成像解释方法"", 《中国石油大学学报(自然科学版)》 * |
张旺: ""基于深度变化的三维重力位场正则化下延方法研究"", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
李中芹: ""基于深度变化的磁场正则化下延方法研究"", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
管志宁: "《地磁场和磁力勘探》", 31 August 2005 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106023085A (zh) * | 2016-06-12 | 2016-10-12 | 西南科技大学 | 一种稀疏约束sar图像重建正则化参数的sure黄金分割自动搜索算法 |
CN106199742A (zh) * | 2016-06-29 | 2016-12-07 | 吉林大学 | 一种频率域航空电磁法2.5维带地形反演方法 |
CN106855642A (zh) * | 2016-11-21 | 2017-06-16 | 中国石油天然气股份有限公司 | 确定基底火山岩分布的方法和装置 |
CN110989002A (zh) * | 2019-11-07 | 2020-04-10 | 吉林大学 | 地空时域电磁系统浅部低阻异常体高精度数据解释方法 |
CN110989002B (zh) * | 2019-11-07 | 2021-01-05 | 吉林大学 | 地空时域电磁系统浅部低阻异常体数据解释方法 |
CN113484919A (zh) * | 2021-06-30 | 2021-10-08 | 长江大学 | 基于卷积神经网络的磁异常反演方法、系统、终端及介质 |
CN113484919B (zh) * | 2021-06-30 | 2022-11-04 | 长江大学 | 基于卷积神经网络的磁异常反演方法、系统、终端及介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100429531C (zh) | 目标最小化的三维电磁快速反演方法 | |
CN105549099A (zh) | 基于全空间正则化下延数据的视磁化强度三维反演方法 | |
WO2022227206A1 (zh) | 一种基于全卷积神经网络的大地电磁反演方法 | |
CN103077558B (zh) | 碳酸盐岩缝洞型油藏大型溶洞储集体分布模型的建模方法 | |
CN108873103A (zh) | 一种结构约束的二维重力梯度和大地电磁联合反演方法 | |
CN105277980A (zh) | 高精度空间和时间任意倍数可变网格有限差分正演方法 | |
CN112949134B (zh) | 基于非结构有限元方法的地-井瞬变电磁反演方法 | |
CN106199742A (zh) | 一种频率域航空电磁法2.5维带地形反演方法 | |
CN109917466B (zh) | 一种基于电磁场垂直分量的电阻率测量方法 | |
CN110286416B (zh) | 一种基于物性函数的快速二维密度反演方法 | |
CN105204078B (zh) | 一种计算三维大地电磁测深主阻抗的方法 | |
CN107305600A (zh) | 最小二乘法电阻率三维近似反演技术 | |
CN102565866A (zh) | 物探二维电测深资料层析法反演处理技术 | |
CN105738952A (zh) | 一种水平井区储层岩石相建模方法 | |
CN104991272A (zh) | 一种针对无井地震反演的地震速度扰动建模方法 | |
CN109490978A (zh) | 一种起伏地层的频率域快速高精度正演方法 | |
CN113406707A (zh) | 一种大地电磁多尺度、多时段探测方法 | |
KR101348788B1 (ko) | 대수적 재구성법을 이용한 3차원 자력역산 방법 | |
Feng et al. | Contrast between 2D inversion and 3D inversion based on 2D high-density resistivity data | |
CN102967889A (zh) | 一种矿床储量的测量方法 | |
CN104570104A (zh) | 一种基于两步法avf的纵横波地震品质因子提取方法 | |
TANG et al. | Determination of borders for resistive oil and gas reservoirs by deviation rate using the hole‐to‐surface resistivity method | |
Deng et al. | 3D finite element numerical simulation of advanced detection in roadway for DC focus method | |
CN103267980B (zh) | 地球电磁场场源修正的物探方法及测量装置 | |
Zheng et al. | Normalization processing and research of multi-field data in reservoir seepage investigation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160504 |
|
RJ01 | Rejection of invention patent application after publication |