CN108304618B - 一种重力数据与大地电磁数据联合反演方法 - Google Patents
一种重力数据与大地电磁数据联合反演方法 Download PDFInfo
- Publication number
- CN108304618B CN108304618B CN201810017213.3A CN201810017213A CN108304618B CN 108304618 B CN108304618 B CN 108304618B CN 201810017213 A CN201810017213 A CN 201810017213A CN 108304618 B CN108304618 B CN 108304618B
- Authority
- CN
- China
- Prior art keywords
- data
- magnetotelluric
- gravimetric
- inverting
- model
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种重力数据与大地电磁数据联合反演方法,包括数据采集步骤、数据筛选步骤、初始模型确定步骤和联合反演步骤;本发明以数据质量为标准,选取联合反演迭代时的初始模型,相对于现有技术中统一给定均匀半空间模型相比,本发明具有如下有益效果:提高了联合反演的效率,减少了反演计算时间;避免了由于数据误差导致的迭代不收敛问题,提高了联合反演结果的精度。
Description
技术领域
本发明涉及一种地球物理反演方法,尤其涉及一种重力数据与大地电磁数据联合反演方法。
背景技术
地球物理反演的基本流程为:先给定初始模型,根据算法搜索模型更新量,更新模型,再计算当前模型计算出的理论数据与实际数据的误差,再重新搜索模型更新量,如此循环,直至误差收敛或是达到最大循环次数。根据大地电磁数据可以反演地下的电阻率模型,重力数据可以反演出地下的密度模型。而所谓联合反演,是指同时对两种或两种以上的地球物理数据进行反演,求取一个使所有参加反演的地球物理数据都满足的统一的地质模型。大地电磁数据横向分辨率较低,但是其丰富的频率成分保证了其较高的纵向分辨率,而重力数据横向分辨率较高,将二者联合可以使两种方法优势互补。联合反演过程中,初始模型一般都是预先给定均匀半空间模型,例如给定地下电阻率为10Ω、密度为1g/cm3的均匀半空间模型,由反演算法自动搜索模型的更新量。这样做有以下弊端:首先,初始模型为均匀半空间,计算速度慢,有可能导致迭代次数过多,效率低下;其次,各种地球物理方法的数据质量不一样,若根据数据质量较差的数据,由均匀半空间模型出发进行第一次反演迭代,则在后续的迭代过程中,造成误差累计,给最终结果造成影响,甚至导致迭代过程不收敛。
发明内容
本发明旨在解决现有技术中存在的问题,提供一种重力数据与大地电磁数据联合反演方法,能够提高联合反演的效率和精度。
本发明的技术方案如下:一种重力数据与大地电磁数据联合反演方法,其特征在于,包括如下步骤:
步骤1),数据采集步骤:采集工作区所有预先选定的测量点的重力数据和大地电磁数据,其中,所述重力数据是指剩余重力数据,所述大地电磁数据是指电场数据和磁场数据;
步骤2),数据筛选步骤:分别计算出所述重力数据和大地电磁数据中跳变数据的数据量占各自数据总量的百分比r1和r2,其中,所述跳变数据是指与相邻采集点的数据误差均超过误差阈值的采集点的数据;
步骤3),初始模型确定步骤:比较r1与r2,若r1≥r2,则使用均匀半空间密度模型作为重力数据反演初始模型,使用第一次重力数据反演结果的密度模型所对应的电阻率模型作为大地电磁数据反演初始模型,并进行大地电磁数据的第一次反演;若r1<r2,则使用均匀半空间电阻率模型作为大地电磁反演初始模型,使用第一次大地电磁反演结果的电阻率模型所对应的密度模型作为重力数据反演的初始模型,并进行重力数据的第一次反演;
步骤4),反演步骤:进行重力数据和大地电磁的联合反演,每进行一次重力数据、大地电磁数据的反演计算并更新密度和电阻率模型后,计算重力数据拟合差Grms和大地电磁数据拟合差MTrms,判断是否同时小于对应的数据拟合差阈值Gth和MTth,若成立,则退出循环,若不成立,则判断迭代次数n是否达到最大的迭代次数Nmax,若成立,则退出循环,若不成立,则继续进行重力数据、大地电磁数据的迭代计算。
优选的,所述步骤1)中的测量点个数为100-2500个,相邻测量点之间的距离为50-500米。
优选的,所述步骤2)中的误差阈值为20%-40%。
优选的,所述步骤3)中的密度模型与电阻率模型之间的对应关系根据对所述工作区的岩石样本统计分析确定。
优选的,所述步骤4)的反演步骤具体为:按顺序分别进行重力数据和大地电磁数据的反演迭代计算,其中,从第二次反演迭代开始,利用最新的电阻率模型约束重力数据的反演,利用最新的密度模型约束大地电磁数据的反演。
优选的,所述约束为交叉梯度约束。
本发明以数据质量为标准,选取联合反演迭代时的初始模型,相对于现有技术中统一给定均匀半空间模型相比,本发明具有如下有益效果:提高了联合反演的效率,减少了反演计算时间;避免了由于数据误差导致的迭代不收敛问题,提高了联合反演结果的精度。
附图说明
图1是重力数据与大地电磁数据联合反演方法的流程图;
图2是步骤4)中反演步骤的流程图。
具体实施方式
下面结合附图对本发明的技术方案作进一步详细说明。
图1是重力数据与大地电磁数据联合反演方法的流程图,图2是步骤4)中反演步骤的流程图。如图1-2所示,本发明的联合反演方法包括:
步骤1),数据采集步骤:预先在工作区选取测量点100个,点距为50米,在测量点处分别采集重力数据和大地电磁数据,其中,重力数据是指剩余重力数据,大地电磁数据是指电场数据和磁场数据。
步骤2),数据筛选步骤:分别计算出所述重力数据和大地电磁数据中跳变数据的数据量占各自数据总量的百分比r1和r2,其中,所述跳变数据是指与相邻采集点的数据误差均超过误差阈值的采集点的数据;误差阈值为20%-40%。
步骤3),初始模型确定步骤:采集工作区的岩石样本至少100个,分别测量样本的密度和电阻率,使用三次样条函数拟合二者之间的关系,由此得到密度模型与电阻率模型之间的对应关系。比较r1与r2,若r1≥r2,则使用均匀半空间密度模型作为重力数据反演初始模型,使用第一次重力数据反演结果的密度模型所对应的电阻率模型作为大地电磁数据反演初始模型,并进行大地电磁数据的第一次反演;若r1<r2,则使用均匀半空间电阻率模型作为大地电磁反演初始模型,使用第一次大地电磁反演结果的电阻率模型所对应的密度模型作为重力数据反演的初始模型,并进行重力数据的第一次反演。
步骤4),反演步骤:进行重力数据和大地电磁的联合反演,按顺序分别进行重力数据和大地电磁数据的反演迭代计算,其中,从第二次反演迭代开始,利用最新的电阻率模型使用交叉梯度函数约束重力数据的反演,利用最新的密度模型使用交叉梯度函数约束大地电磁数据的反演。每进行一次重力数据、大地电磁数据的反演计算并更新密度和电阻率模型后,计算重力数据拟合差Grms和大地电磁数据拟合差MTrms,判断是否同时小于对应的数据拟合差阈值Gth和MTth,若成立,则退出循环,若不成立,则判断迭代次数n是否达到最大的迭代次数Nmax,若成立,则退出循环,若不成立,则继续进行重力数据、大地电磁数据的迭代计算。
表1-2给出了本发明与现有技术的结果对比:
表1r1=8%,r2=11%时本发明与现有技术对比
反演用时(h) | 重力数据拟合差(%) | 大地电磁数据拟合差(%) | |
本发明 | 20.2 | 1.5 | 2.4 |
现有技术 | 34.7 | 2.7 | 3.1 |
表2r1=8%,r2=26%时本发明与现有技术对比
反演用时(h) | 重力数据拟合差(%) | 大地电磁数据拟合差(%) | |
本发明 | 24.9 | 1.7 | 3.3 |
现有技术 | 不收敛 | — | — |
由表中可以看到,本发明以数据质量为标准,选取联合反演迭代时的初始模型,相对于现有技术中统一给定均匀半空间模型相比,本发明具有如下有益效果:提高了联合反演的效率,减少了反演计算时间;避免了由于数据误差导致的迭代不收敛问题,提高了联合反演结果的精度。
Claims (5)
1.一种重力数据与大地电磁数据联合反演方法,其特征在于,包括如下步骤:
步骤1),数据采集步骤:采集工作区所有预先选定的测量点的重力数据和大地电磁数据,其中,所述重力数据是指剩余重力数据,所述大地电磁数据是指电场数据和磁场数据;
步骤2),数据筛选步骤:分别计算出所述重力数据和大地电磁数据中跳变数据的数据量占各自数据总量的百分比r1和r2,其中,所述跳变数据是指与相邻采集点的数据误差均超过误差阈值的采集点的数据;
步骤3),初始模型确定步骤:采集工作区的岩石样本至少100个,分别测量样本的密度和电阻率,使用三次样条函数拟合二者之间的关系,由此得到岩石样本的密度模型与岩石样本的电阻率模型之间的对应关系;比较r1与r2,若r1≥r2,则使用均匀半空间密度模型作为重力数据反演初始模型,使用第一次重力数据反演结果的密度模型所对应的电阻率模型作为大地电磁数据反演初始模型,并进行大地电磁数据的第一次反演;若r1<r2,则使用均匀半空间电阻率模型作为大地电磁反演初始模型,使用第一次大地电磁反演结果的电阻率模型所对应的密度模型作为重力数据反演的初始模型,并进行重力数据的第一次反演;
步骤4),反演步骤:进行重力数据和大地电磁的联合反演,每进行一次重力数据、大地电磁数据的反演计算并更新密度和电阻率模型后,计算重力数据拟合差Grms和大地电磁数据拟合差MTrms,判断是否同时小于对应的数据拟合差阈值Gth和MTth,若成立,则退出循环,若不成立,则判断迭代次数n是否达到最大的迭代次数Nmax,若成立,则退出循环,若不成立,则继续进行重力数据、大地电磁数据的迭代计算;
所述步骤3)中的岩石样本的密度模型与岩石样本的电阻率模型之间的对应关系根据对所述工作区的岩石样本统计分析确定。
2.如权利要求1所述的重力数据与大地电磁数据联合反演方法,其特征在于,所述步骤1)中的测量点个数为100-2500个,相邻测量点之间的距离为50-500米。
3.如权利要求1所述的重力数据与大地电磁数据联合反演方法,其特征在于,所述步骤2)中的误差阈值为20%-40%。
4.如权利要求1所述的重力数据与大地电磁数据联合反演方法,其特征在于,所述步骤4)的反演步骤具体为:按顺序分别进行重力数据和大地电磁数据的反演迭代计算,其中,从第二次反演迭代开始,利用最新的电阻率模型约束重力数据的反演,利用最新的密度模型约束大地电磁数据的反演。
5.如权利要求4所述的重力数据与大地电磁数据联合反演方法,其特征在于,所述约束为交叉梯度约束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810017213.3A CN108304618B (zh) | 2018-01-05 | 2018-01-05 | 一种重力数据与大地电磁数据联合反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810017213.3A CN108304618B (zh) | 2018-01-05 | 2018-01-05 | 一种重力数据与大地电磁数据联合反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108304618A CN108304618A (zh) | 2018-07-20 |
CN108304618B true CN108304618B (zh) | 2019-05-24 |
Family
ID=62868823
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810017213.3A Active CN108304618B (zh) | 2018-01-05 | 2018-01-05 | 一种重力数据与大地电磁数据联合反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108304618B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111158059B (zh) * | 2020-01-08 | 2021-04-27 | 中国海洋大学 | 基于三次b样条函数的重力反演方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2062071B1 (en) * | 2006-09-04 | 2014-10-22 | Geosystem S.r.l. | Method for building velocity models for pre-stack depth migration via the simultaneous joint inversion of seismic, gravity and magnetotelluric data |
CN104216006A (zh) * | 2013-06-04 | 2014-12-17 | 中国石油化工股份有限公司 | 利用重磁电震同步联合参数模型进行提高成像质量的方法 |
CN104375195A (zh) * | 2013-08-15 | 2015-02-25 | 中国石油天然气集团公司 | 时频电磁的多源多分量三维联合反演方法 |
WO2017048445A1 (en) * | 2015-09-15 | 2017-03-23 | Exxonmobil Upstream Research Company | Accelerated occam inversion using model remapping and jacobian matrix decomposition |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104102814B (zh) * | 2014-06-11 | 2017-07-11 | 中国科学院地质与地球物理研究所 | 一种基于大地电磁数据反演电阻率和磁化率的方法及系统 |
CN104537714B (zh) * | 2015-01-07 | 2019-10-29 | 吉林大学 | 磁共振与瞬变电磁空间约束联合反演方法 |
CN106501867B (zh) * | 2016-10-19 | 2019-03-15 | 中国科学院电子学研究所 | 一种基于横向平滑约束的瞬变电磁反演方法 |
-
2018
- 2018-01-05 CN CN201810017213.3A patent/CN108304618B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2062071B1 (en) * | 2006-09-04 | 2014-10-22 | Geosystem S.r.l. | Method for building velocity models for pre-stack depth migration via the simultaneous joint inversion of seismic, gravity and magnetotelluric data |
CN104216006A (zh) * | 2013-06-04 | 2014-12-17 | 中国石油化工股份有限公司 | 利用重磁电震同步联合参数模型进行提高成像质量的方法 |
CN104375195A (zh) * | 2013-08-15 | 2015-02-25 | 中国石油天然气集团公司 | 时频电磁的多源多分量三维联合反演方法 |
WO2017048445A1 (en) * | 2015-09-15 | 2017-03-23 | Exxonmobil Upstream Research Company | Accelerated occam inversion using model remapping and jacobian matrix decomposition |
Non-Patent Citations (2)
Title |
---|
"基于GPU 的二维地震波全波形反演理论测试";崔丛越等;《中国地球科学联合学术年会2017》;20171015;正文第1810页第3节 |
"基于交叉梯度的重磁和大地电磁三维联合反演研究";付少帅;《中国优秀硕士学位论文全文数据库 基础科学辑》;20160115(第1期);正文第28-43页第4.2节,图4-3,4-4 |
Also Published As
Publication number | Publication date |
---|---|
CN108304618A (zh) | 2018-07-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104750896B (zh) | 一种缝洞型碳酸盐岩油藏数值模拟方法 | |
CN104537714B (zh) | 磁共振与瞬变电磁空间约束联合反演方法 | |
CN103064124B (zh) | 一种校正电磁勘探中地形影响的比值方法 | |
CN101769151B (zh) | 基于油水前缘的电阻率数据采集处理方法 | |
CN105114063B (zh) | 一种斜井双侧向测井曲线校正方法 | |
CN104727813B (zh) | 一种源储共生型地层的孔隙度测量方法 | |
CN109799540B (zh) | 基于地质信息约束下的火山岩型铀矿床磁化率反演方法 | |
CN104749617A (zh) | 一种多尺度裂缝储层正演模型建立方法 | |
CN106199733B (zh) | 用于评价地热富水地质异常体的方法及装置 | |
CN106650192B (zh) | 一种火山岩型铀矿床磁性界面反演方法 | |
CN106443774A (zh) | 一种提高起伏地表叠前深度偏移成像精度的方法 | |
CN109597136A (zh) | 一种矿井全空间瞬变电磁数据处理方法 | |
CN102338887B (zh) | 不规则尺寸空变网格层析成像静校正方法 | |
CN108304618B (zh) | 一种重力数据与大地电磁数据联合反演方法 | |
CN109655890B (zh) | 一种深度域浅中深层联合层析反演速度建模方法及系统 | |
CN110532507B (zh) | 一种提高水平井致密油储层钻遇率的方法 | |
CN103477248A (zh) | 用于估算碳氢化合物生产区域中至少一个参数值的计算方法,以在该区域规划和实施作业 | |
CN111965712B (zh) | 一种校正可控源音频大地电磁法静态效应的方法 | |
CN103543478A (zh) | 地质形态插值的km方法 | |
CN110765665B (zh) | 一种地学动态建模方法及系统 | |
CN105464650A (zh) | 一种随钻测井解释方法 | |
CN113516347B (zh) | 基于地球物理与水文数据融合的地下水水质评估方法 | |
CN104182558A (zh) | 一种缝洞野外露头水驱油数值模拟方法 | |
CN105242317B (zh) | 一种纵波速度的确定方法及装置 | |
Pirrone et al. | E2GPR-Edit your geometry, execute GprMax2D and plot the results! |
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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20221009 Address after: 276000 Hexi Village, Zhulu Town, Lingang Economic Development Zone, Linyi City, Shandong Province Patentee after: Linyi Lingang Water Affairs Co.,Ltd. Address before: 318001 room 1406, building 2, Boao business sea, Jiaojiang District, Taizhou, Zhejiang. Patentee before: TAIZHOU CHUANGXING ENVIRONMENTAL PROTECTION SCIENCE & TECHNOLOGY Co.,Ltd. |
|
TR01 | Transfer of patent right |