CN110286416B - 一种基于物性函数的快速二维密度反演方法 - Google Patents

一种基于物性函数的快速二维密度反演方法 Download PDF

Info

Publication number
CN110286416B
CN110286416B CN201910393070.0A CN201910393070A CN110286416B CN 110286416 B CN110286416 B CN 110286416B CN 201910393070 A CN201910393070 A CN 201910393070A CN 110286416 B CN110286416 B CN 110286416B
Authority
CN
China
Prior art keywords
inversion
physical property
function
gravity
property function
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
CN201910393070.0A
Other languages
English (en)
Other versions
CN110286416A (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 CN201910393070.0A priority Critical patent/CN110286416B/zh
Publication of CN110286416A publication Critical patent/CN110286416A/zh
Application granted granted Critical
Publication of CN110286416B publication Critical patent/CN110286416B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种基于物性函数的快速二维密度反演方法，包括以下步骤：S1，野外实测采集重力数据；S2，依据重力数据或已有资料确定反演地区的物性函数阶次；S3，依据物性函数分布，将地下半空间中反演区域横向剖分；S4，建立相应的核函数矩阵；S5，对参数反演结果进行正演，检验反演拟合程度。本发明引入物性函数，减少了运算时的内存占用并提高了运算速度；在引入物性函数时对地下半空间密度模型的构建方式；反演速度提高；计算过程中占用更少的内存空间。

Description

一种基于物性函数的快速二维密度反演方法
技术领域
本发明涉及二维密度反演方法，具体涉及一种基于物性函数的快速二维密度反演方法。
背景技术
目前，我国矿产资源接替基地面临的主要找矿难题是：老矿山深部和各类隐伏区的探矿难度大，急需先进、高效的理论和技术方法指导深部找矿。我国大部分金属矿山位于地形条件相对较好的地区，探查和开采深度均停留在500米以上范围。
随着卫星、航空等勘探手段的综合应用，以及数据精度和勘探精度要求的不断提高，重磁勘探在深部探测方面的作用显著增强。但传统的密度反演方法运算量大、计算速度慢，难以满足当前大数据量下高精度深部反演的需求。
重力勘探以勘探目标与围岩间密度差异为物性基础，通过测量重力异常来研究地下空间地质构造特征。因其具有经济、勘探深度大及快速获得面积上信息等优点，在地球深部构造探测、区域地质构造单元划分、沉积盆地圈定、固体矿产及油气资源勘查等领域获得广泛应用。近年来，随着硬件设备的升级和计算机技术的快速发展，重力勘探装备精度与效率都有了很大提高，同时数据处理与解释方法也由传统的定性解释逐步向定量解释发展。
密度反演是根据观测的异常数据对地下地质体密度分布的定量计算，进而可以用来推测地质体的空间分布体积，另外还能够为地震解释和钻井工程在前期提供有利的密度分布信息，进一步提高地球物理解释的精度。
现有技术对于大数据量的高精度反演速度慢，并且对于计算机性能有较高的要求。
需要发明一种基于物性函数的快速二维密度反演方法，其能够快速的进行密度反演并得到高精度的解释结果。
发明内容
本发明的主要目的在于提供一种基于物性函数的快速二维密度反演方法。
本发明采用的技术方案是：一种基于物性函数的快速二维密度反演方法，包括以下步骤：
S1，野外实测采集重力数据；
S2，依据重力数据或已有资料确定反演地区的物性函数阶次；
S3，依据物性函数分布，将地下半空间中反演区域横向剖分；
S4，建立相应的核函数矩阵；
S5，对参数反演结果进行正演，检验反演拟合程度。
进一步地，所述步骤S1具体为：
野外实测采集重力数据，均匀采样后得到适当间距的重力异常数据
Figure 677985DEST_PATH_IMAGE001
；
更进一步地，所述步骤S2具体为：
依据重力数据或已有资料确定反演地区的物性函数阶次
Figure 245232DEST_PATH_IMAGE002
，物性函数表示为：
Figure 668123DEST_PATH_IMAGE003
其中：
Figure 586401DEST_PATH_IMAGE004
点处的密度；
Figure 487361DEST_PATH_IMAGE005
均为物性函数参数；
更进一步地，所述步骤S3具体为：
依据物性函数分布，将地下半空间中反演区域横向剖分成
Figure 440273DEST_PATH_IMAGE006
个块体单元，
Figure 768486DEST_PATH_IMAGE007
点的重力异常等效于各块体单元的物性函数的各项在该点处重力效应的总和；
更进一步地，所述步骤S4具体为：
建立相应的核函数矩阵
Figure 675525DEST_PATH_IMAGE008
，其中：
Figure 114596DEST_PATH_IMAGE009
，其构建方法如下：
对于物性函数，剖分方式取
Figure 922015DEST_PATH_IMAGE010
个点对其进行拟合计算：
令：
Figure 686709DEST_PATH_IMAGE011
，则各点密度可表示为：
Figure 313999DEST_PATH_IMAGE012
其中：
Figure 822341DEST_PATH_IMAGE013
Figure 251311DEST_PATH_IMAGE014
Figure 921327DEST_PATH_IMAGE015
Figure 35913DEST_PATH_IMAGE016
Figure 82366DEST_PATH_IMAGE017
Figure 864378DEST_PATH_IMAGE018
令：
Figure 970874DEST_PATH_IMAGE019
Figure 808642DEST_PATH_IMAGE020
则各点的重力异常表示为：
Figure 393207DEST_PATH_IMAGE021
Figure 307123DEST_PATH_IMAGE023
Figure 396301DEST_PATH_IMAGE024
更进一步地，所述步骤S5具体为：
对参数反演结果进行正演，检验反演拟合程度。
本发明的优点：
本发明引入物性函数，减少了运算时的内存占用并提高了运算速度；在引入物性函数时对地下半空间密度模型的构建方式；反演速度提高；计算过程中占用更少的内存空间。
本发明通过引入物性函数减少了所需反演的未知量个数，进而减少了反演过程中的数据量，相应的提高了反演的速度，并一定程度的降低了对计算机性能的要求。
本发明创造性的利用物性函数对改进了密度反演方法，极大地提高了反演的效率，在实际地区密度反演中得到了良好的应用。
除了上面所描述的目的、特征和优点之外，本发明还有其它的目的、特征和优点。下面将参照图，对本发明作进一步详细的说明。
附图说明
构成本申请的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。
图1是本发明的反演流程图；
图2是本发明的地下空间离散化正演网格；
图3是本发明的拟合随迭代变化图；
图4是本发明的密度反演成果图；
图5是本发明的正演结果与原始重力异常对比图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
参考图1至图5，如图1至图5所示，一种基于物性函数的快速二维密度反演方法，包括以下步骤：
S1，野外实测采集重力数据；
S2，依据重力数据或已有资料确定反演地区的物性函数阶次；
S3，依据物性函数分布，将地下半空间中反演区域横向剖分；
S4，建立相应的核函数矩阵；
S5，对参数反演结果进行正演，检验反演拟合程度。
进一步地，所述步骤S1具体为：
野外实测采集重力数据，均匀采样后得到适当间距的重力异常数据
Figure 784557DEST_PATH_IMAGE001
；
更进一步地，所述步骤S2具体为：
依据重力数据或已有资料确定反演地区的物性函数阶次
Figure 777046DEST_PATH_IMAGE002
，物性函数表示为：
Figure 225345DEST_PATH_IMAGE003
其中：
Figure 536241DEST_PATH_IMAGE004
点处的密度；
Figure 728188DEST_PATH_IMAGE005
均为物性函数参数；
更进一步地，所述步骤S3具体为：
依据物性函数分布，将地下半空间中反演区域横向剖分成
Figure 73718DEST_PATH_IMAGE006
个块体单元，
Figure 692919DEST_PATH_IMAGE007
点的重力异常等效于各块体单元的物性函数的各项在该点处重力效应的总和；
更进一步地，所述步骤S4具体为：
建立相应的核函数矩阵
Figure 258154DEST_PATH_IMAGE008
，其中：
Figure 253792DEST_PATH_IMAGE009
，其构建方法如下：
对于物性函数，剖分方式取
Figure 453829DEST_PATH_IMAGE010
个点对其进行拟合计算：
令：
Figure 243931DEST_PATH_IMAGE011
，则各点密度可表示为：
Figure 529419DEST_PATH_IMAGE012
其中：
Figure 328748DEST_PATH_IMAGE013
Figure 383291DEST_PATH_IMAGE014
Figure 857478DEST_PATH_IMAGE015
Figure 630261DEST_PATH_IMAGE016
Figure 967702DEST_PATH_IMAGE017
Figure 876752DEST_PATH_IMAGE018
令：
Figure 274235DEST_PATH_IMAGE019
Figure 534315DEST_PATH_IMAGE020
则各点的重力异常表示为：
Figure 911332DEST_PATH_IMAGE021
Figure 940468DEST_PATH_IMAGE025
Figure 243274DEST_PATH_IMAGE026
Figure 990650DEST_PATH_IMAGE027
更进一步地，所述步骤S5具体为：
对参数反演结果进行正演，检验反演拟合程度。
本发明的具体实施方案：
1）采集某地区重力剖面数据，并对其进行均匀采样；
2）采用2阶物性函数对其进行拟合计算，并按照如图2所示方式建立物性函数模型空间；
3）代入数据，令
Figure DEST_PATH_IMAGE028
构建核函数矩阵；
4）迭代计算并在拟合差足够小时停止迭代；
5）绘制反演所得密度空间分布图像（如图3所示）；
6）
Figure 201051DEST_PATH_IMAGE029
时，反演所得密度空间分布图像（如图4所示）。
下表为本发明与传统物性反演方法，在采用同等反演精度下（相同采样点距、相同成图精度、相同的截止条件），反演相同数据时采用同一的计算机的同一软件反演时的运行时间以及内存占用率对比表：
运行时间 内存占用率
本发明 37.9219s 25%
传统方法 990.2031s 31%
本发明引入物性函数，减少了运算时的内存占用并提高了运算速度；在引入物性函数时对地下半空间密度模型的构建方式；反演速度提高；计算过程中占用更少的内存空间。
本发明通过引入物性函数减少了所需反演的未知量个数，进而减少了反演过程中的数据量，相应的提高了反演的速度，并一定程度的降低了对计算机性能的要求。
本发明创造性的利用物性函数对改进了密度反演方法，极大地提高了反演的效率，在实际地区密度反演中得到了良好的应用。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (1)

1.一种基于物性函数的快速二维密度反演方法，其特征在于，包括以
下步骤：
S1，野外实测采集重力数据；
S2，依据重力数据或已有资料确定反演地区的物性函数阶次；
S3，依据物性函数分布，将地下半空间中反演区域横向剖分；
S4，建立相应的核函数矩阵；
S5，对参数反演结果进行正演，检验反演拟合程度;
所述步骤S1具体为：
野外实测采集重力数据，均匀采样后得到适当间距的重力异常数据
Figure 3842DEST_PATH_IMAGE001
；
所述步骤S2具体为：
依据重力数据或已有资料确定反演地区的物性函数阶次
Figure 186561DEST_PATH_IMAGE002
，物性函数表示为：
Figure 497457DEST_PATH_IMAGE003
其中：
Figure 423825DEST_PATH_IMAGE004
点处的密度；
Figure 566093DEST_PATH_IMAGE005
均为物性函数参数；
所述步骤S3具体为：
依据物性函数分布，将地下半空间中反演区域横向剖分成
Figure 654135DEST_PATH_IMAGE006
个块体单元，
Figure 452326DEST_PATH_IMAGE007
点的重力异常等效于各块体单元的物性函数的各项在该点处重力效应的总和；
所述步骤S4具体为：
建立相应的核函数矩阵
Figure 683850DEST_PATH_IMAGE008
，其中：
Figure 680625DEST_PATH_IMAGE009
，其构建方法如下：
对于物性函数，剖分方式取
Figure 205147DEST_PATH_IMAGE010
个点对其进行拟合计算：
令：
Figure 225056DEST_PATH_IMAGE011
，则各点密度可表示为：
Figure 289963DEST_PATH_IMAGE012
其中：
Figure 344507DEST_PATH_IMAGE013
Figure 541396DEST_PATH_IMAGE014
Figure 314179DEST_PATH_IMAGE015
Figure 917199DEST_PATH_IMAGE016
Figure 91829DEST_PATH_IMAGE017
Figure 958153DEST_PATH_IMAGE018
令：
Figure 952654DEST_PATH_IMAGE019
Figure 562627DEST_PATH_IMAGE020
则各点的重力异常表示为：
Figure 889965DEST_PATH_IMAGE021
Figure 192771DEST_PATH_IMAGE022
Figure 408989DEST_PATH_IMAGE024
Figure DEST_PATH_IMAGE026
所述步骤S5具体为：
对参数反演结果进行正演，检验反演拟合程度。
CN201910393070.0A 2019-05-13 2019-05-13 一种基于物性函数的快速二维密度反演方法 Expired - Fee Related CN110286416B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910393070.0A CN110286416B (zh) 2019-05-13 2019-05-13 一种基于物性函数的快速二维密度反演方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910393070.0A CN110286416B (zh) 2019-05-13 2019-05-13 一种基于物性函数的快速二维密度反演方法

Publications (2)

Publication Number Publication Date
CN110286416A CN110286416A (zh) 2019-09-27
CN110286416B true CN110286416B (zh) 2020-12-22

Family

ID=68001595

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910393070.0A Expired - Fee Related CN110286416B (zh) 2019-05-13 2019-05-13 一种基于物性函数的快速二维密度反演方法

Country Status (1)

Country Link
CN (1) CN110286416B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112861065B (zh) * 2021-02-10 2022-07-08 湖北省地震局(中国地震局地震研究所) 一种自适应的震源区重力密度变化实时反演计算方法
CN113311494B (zh) * 2021-05-26 2022-04-26 中南大学 一种卫星重力场反演方法
CN113514899A (zh) * 2021-07-12 2021-10-19 吉林大学 一种重力梯度的自适应剖分反演方法
CN113486591B (zh) * 2021-07-13 2022-04-19 吉林大学 一种卷积神经网络结果的重力多参量数据密度加权反演方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003023447A2 (en) * 2001-09-07 2003-03-20 Conocophillips Company A nonlinear constrained inversion method to determine base of salt interface from gravity and gravity tensor data
CN101661115A (zh) * 2008-08-29 2010-03-03 中国石油集团东方地球物理勘探有限责任公司 基于标准格架的快速三维重力、磁力物性反演的方法
WO2017131848A1 (en) * 2016-01-26 2017-08-03 Saudi Arabian Oil Company Gravity modeling a rifted continental margin
CN107561592A (zh) * 2017-09-11 2018-01-09 中南大学 一种密度为多项式的任意多面体的重力场正演方法
CN108229082A (zh) * 2018-04-12 2018-06-29 吉林大学 一种基于数据空间快速计算的联合反演方法
CN108490496A (zh) * 2018-03-26 2018-09-04 中国石油化工股份有限公司 基于拟径向基函数神经网络的重力场密度反演方法
CN108873103A (zh) * 2018-09-14 2018-11-23 吉林大学 一种结构约束的二维重力梯度和大地电磁联合反演方法
CN109471190A (zh) * 2018-11-12 2019-03-15 吉林大学 一种不同高度重力数据和井中重力数据联合反演方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8571842B2 (en) * 2010-06-21 2013-10-29 Schlumberger Technology Corporation Method of determining parameter from sparse measurement data
CN103913774B (zh) * 2014-04-02 2017-03-01 西南石油大学 基于微地震事件的储层地质力学参数反演方法
CN104866653B (zh) * 2015-04-29 2017-10-24 中国地质科学院矿产资源研究所 一种获取地下三维密度结构的方法
CN107491411B (zh) * 2017-06-23 2020-07-17 中国海洋大学 基于n阶多项式密度函数的重力异常反演方法
CN107748399B (zh) * 2017-09-12 2019-12-31 中国石油化工股份有限公司 利用重力界面反演识别山前带深部构造层方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003023447A2 (en) * 2001-09-07 2003-03-20 Conocophillips Company A nonlinear constrained inversion method to determine base of salt interface from gravity and gravity tensor data
CN101661115A (zh) * 2008-08-29 2010-03-03 中国石油集团东方地球物理勘探有限责任公司 基于标准格架的快速三维重力、磁力物性反演的方法
WO2017131848A1 (en) * 2016-01-26 2017-08-03 Saudi Arabian Oil Company Gravity modeling a rifted continental margin
CN107561592A (zh) * 2017-09-11 2018-01-09 中南大学 一种密度为多项式的任意多面体的重力场正演方法
CN108490496A (zh) * 2018-03-26 2018-09-04 中国石油化工股份有限公司 基于拟径向基函数神经网络的重力场密度反演方法
CN108229082A (zh) * 2018-04-12 2018-06-29 吉林大学 一种基于数据空间快速计算的联合反演方法
CN108873103A (zh) * 2018-09-14 2018-11-23 吉林大学 一种结构约束的二维重力梯度和大地电磁联合反演方法
CN109471190A (zh) * 2018-11-12 2019-03-15 吉林大学 一种不同高度重力数据和井中重力数据联合反演方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
"Compact gravity inversion";Last B J等;《Geophysics》;19831231;第46卷(第6期);第713-721页 *
"基于Born敏感核函数的速度、密度双参数全波形反演";杨积忠 等;《地球物理学报》;20160331;第59卷(第31期);第1082-1094页 *
"重力及其梯度张量数据高精度成像和密度反演方法";王泰涵;《中国博士学位论文全文数据库 基础科学辑》;20190115(第01（2019）期);第A011-159页 *
"重力资料聚焦反演成像改进及其应用";栗学磊 等;《世界地质》;20130630;第32卷(第2期);第366-371页 *

Also Published As

Publication number Publication date
CN110286416A (zh) 2019-09-27

Similar Documents

Publication Publication Date Title
CN110286416B (zh) 一种基于物性函数的快速二维密度反演方法
Han et al. A 3D fracture network model for the undisturbed rock mass at the Songta dam site based on small samples
CN112949134B (zh) 基于非结构有限元方法的地-井瞬变电磁反演方法
CN107748399B (zh) 利用重力界面反演识别山前带深部构造层方法
KR20160024232A (ko) 지반정보에 기초한 3차원 공간 모델링 방법
Li et al. Part II: A demonstration of integrating multiple-scale 3D modelling into GIS-based prospectivity analysis: A case study of the Huayuan-Malichang district, China
CN111221048B (zh) 基于跨孔电阻率ct多尺度反演的孤石边界识别与成像方法
Koike et al. Incorporation of fracture directions into 3D geostatistical methods for a rock fracture system
CN109459787B (zh) 基于地震槽波全波形反演的煤矿井下构造成像方法及系统
CN110119586B (zh) 轴向电导率各向异性瞬变电磁三分量三维fdtd正演方法
CN110133748B (zh) 一种集成碱交代型铀矿深部成矿信息的方法
CN115292890A (zh) 基于多源辅助数据开发的场地土壤污染物浓度三维空间预测方法
Li et al. Uncertainty visualisation of a 3D geological geometry model and its application in GIS-based mineral resource assessment: a case study in Huayuan District, Northwestern Hunan Province, China
Amigun et al. Euler 3-D deconvolution of analytical signal of magnetic anomalies over iron ore deposit in Okene, Nigeria
CN113406707A (zh) 一种大地电磁多尺度、多时段探测方法
CN109490978A (zh) 一种起伏地层的频率域快速高精度正演方法
CN112882124B (zh) 一种三维铀成矿有利部位的圈定方法
CN105242317B (zh) 一种纵波速度的确定方法及装置
CN105572732A (zh) 一种属性变化率逐次晋级的裂缝发育带检测方法
CN114280688A (zh) 一种基于无人机平台的地下人工空洞多物理场综合探测系统及其运行方法
CN110632665B (zh) 一种砂体构型量化分析方法
CN113934806A (zh) 一种土壤地球化学测量数据精细化处理方法
CN116068663A (zh) 基于磁震联合低频建模的火成岩波阻抗反演方法
Al-Hadithi et al. Using source parameter imaging technique to the aeromagnetic data to estimate the basement depth of Tharthar Lake and surrounding area in Central Iraq
CN113031076A (zh) 一种在电阻率法超前探测中考虑围岩各向异性影响的方法

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

Termination date: 20210513

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