CN110286416A - 一种基于物性函数的快速二维密度反演方法 - Google Patents
一种基于物性函数的快速二维密度反演方法 Download PDFInfo
- Publication number
- CN110286416A CN110286416A CN201910393070.0A CN201910393070A CN110286416A CN 110286416 A CN110286416 A CN 110286416A CN 201910393070 A CN201910393070 A CN 201910393070A CN 110286416 A CN110286416 A CN 110286416A
- Authority
- CN
- China
- Prior art keywords
- inversion
- physical property
- function
- fast
- 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.)
- Granted
Links
- 230000000704 physical effect Effects 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 32
- 230000006870 function Effects 0.000 claims abstract description 52
- 239000011159 matrix material Substances 0.000 claims abstract description 8
- 230000005484 gravity Effects 0.000 claims description 31
- 238000005259 measurement Methods 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 5
- 238000005070 sampling Methods 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 5
- 238000004088 simulation Methods 0.000 abstract 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 3
- 239000011707 mineral Substances 0.000 description 3
- 230000006872 improvement Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
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具体为:
野外实测采集重力数据,均匀采样后得到适当间距的重力异常数据()。
根据权利要求1所述的基于物性函数的快速二维密度反演方法,其
特征在于,所述步骤S2具体为:
依据重力数据或已有资料确定反演地区的物性函数阶次,物性函数表示为:
其中:为点处的密度;,均为物性函数参数。
根据权利要求1所述的基于物性函数的快速二维密度反演方法,其
特征在于,所述步骤S3具体为:
依据物性函数分布,将地下半空间中反演区域横向剖分成个块体单元,点的重力异常等效于各块体单元的物性函数的各项在该点处重力效应的总和。
根据权利要求1所述的基于物性函数的快速二维密度反演方法,其
特征在于,所述步骤S4具体为:
建立相应的核函数矩阵,其中:,其构建方法如下:
对于物性函数,依据图二所示剖分方式取个点对其进行拟合计算:
令:则各点密度可表示为:
其中:
令:
则各点的重力异常表示为:
其中:
其中,为万有引力常数;;;为第个异常数据的横坐标;;;距离单位均为。
更进一步地,所述步骤S5具体为:
对参数反演结果进行正演,检验反演拟合程度,并进行多次迭代反演,直至得到高精度的反演结果。
本发明的优点:
本发明引入物性函数,减少了运算时的内存占用并提高了运算速度;在引入物性函数时对地下半空间密度模型的构建方式;反演速度提高;计算过程中占用更少的内存空间。
本发明通过引入物性函数减少了所需反演的未知量个数,进而减少了反演过程中的数据量,相应的提高了反演的速度,并一定程度的降低了对计算机性能的要求。
本发明创造性的利用物性函数对改进了密度反演方法,极大地提高了反演的效率,在实际地区密度反演中得到了良好的应用。
除了上面所描述的目的、特征和优点之外,本发明还有其它的目的、特征和优点。下面将参照图,对本发明作进一步详细的说明。
附图说明
构成本申请的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。
图1是本发明的反演流程图;
图2是本发明的地下空间离散化正演网格;
图3是本发明的拟合随迭代变化图;
图4是本发明的密度反演成果图;
图5是本发明的正演结果与原始重力异常对比图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
参考图1至图5,如图1至图5所示,一种基于物性函数的快速二维密度反演方法,包括以下步骤:
S1,野外实测采集重力数据;
S2,依据重力数据或已有资料确定反演地区的物性函数阶次;
S3,依据物性函数分布,将地下半空间中反演区域横向剖分;
S4,建立相应的核函数矩阵;
S5,对参数反演结果进行正演,检验反演拟合程度。
进一步地,所述步骤S1具体为:
野外实测采集重力数据,均匀采样后得到适当间距的重力异常数据( )。
更进一步地,所述步骤S2具体为:
依据重力数据或已有资料确定反演地区的物性函数阶次(一般为取2),物性函数表示为:
其中:为点处的密度;,均为物性函数参数。
参考图2,如图2所示,,所述步骤S3具体为:
依据物性函数分布,将地下半空间中反演区域横向剖分成个块体单元,点的重力异常等效于各块体单元的物性函数的各项在该点处重力效应的总和。
所述步骤S4具体为:
建立相应的核函数矩阵,其中:,其构建方法如下:
对于物性函数,依据图二所示剖分方式取个点对其进行拟合计算:
令:则各点密度可表示为:
其中:
令:
则各点的重力异常表示为:
其中:
其中,为万有引力常数;;;为第个异常数据的横坐标;;;距离单位均为。
所述步骤S5具体为:
对参数反演结果进行正演,检验反演拟合程度,并进行多次迭代反演,直至得到高精度的反演结果。
本发明的具体实施方案:
1)采集某地区重力剖面数据,并对其进行均匀采样;
2)采用2阶物性函数对其进行拟合计算,并按照如图2所示方式建立物性函数模型空间;
3)代入数据,令构建核函数矩阵;
4)迭代计算并在拟合差足够小时停止迭代;
5)绘制反演所得密度空间分布图像(如图3所示);
6)时,反演所得密度空间分布图像(如图4所示)。
下表为本发明与传统物性反演方法,在采用同等反演精度下(相同采样点距、相同成图精度、相同的截止条件),反演相同数据时采用同一的计算机的同一软件反演时的运行时间以及内存占用率对比表:
运行时间 | 内存占用率 | |
本发明 | 37.9219s | 25% |
传统方法 | 990.2031s | 31% |
本发明引入物性函数,减少了运算时的内存占用并提高了运算速度;在引入物性函数时对地下半空间密度模型的构建方式;反演速度提高;计算过程中占用更少的内存空间。
本发明通过引入物性函数减少了所需反演的未知量个数,进而减少了反演过程中的数据量,相应的提高了反演的速度,并一定程度的降低了对计算机性能的要求。
本发明创造性的利用物性函数对改进了密度反演方法,极大地提高了反演的效率,在实际地区密度反演中得到了良好的应用。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种基于物性函数的快速二维密度反演方法,其特征在于,包括以
下步骤:
S1,野外实测采集重力数据;
S2,依据重力数据或已有资料确定反演地区的物性函数阶次;
S3,依据物性函数分布,将地下半空间中反演区域横向剖分;
S4,建立相应的核函数矩阵;
S5,对参数反演结果进行正演,检验反演拟合程度。
2.根据权利要求1所述的基于物性函数的快速二维密度反演方法,其
特征在于,所述步骤S1具体为:
野外实测采集重力数据,均匀采样后得到适当间距的重力异常数据()。
3.根据权利要求1所述的基于物性函数的快速二维密度反演方法,其
特征在于,所述步骤S2具体为:
依据重力数据或已有资料确定反演地区的物性函数阶次,物性函数表示为:
其中:为点处的密度;,均为物性函数参数。
4.根据权利要求1所述的基于物性函数的快速二维密度反演方法,其
特征在于,所述步骤S3具体为:
依据物性函数分布,将地下半空间中反演区域横向剖分成个块体单元,点的重力异常等效于各块体单元的物性函数的各项在该点处重力效应的总和。
5.根据权利要求1所述的基于物性函数的快速二维密度反演方法,其
特征在于,所述步骤S4具体为:
建立相应的核函数矩阵,其中:,其构建方法如下:
对于物性函数,依据图二所示剖分方式取个点对其进行拟合计算:
令:则各点密度可表示为:
其中:
令:
则各点的重力异常表示为:
其中:
其中,为万有引力常数;;;为第个异常数据的横坐标;;;距离单位均为。
6.根据权利要求1所述的基于物性函数的快速二维密度反演方法,其
特征在于,所述步骤S5具体为:
对参数反演结果进行正演,检验反演拟合程度,并进行多次迭代反演,直至得到高精度的反演结果。
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 true CN110286416A (zh) | 2019-09-27 |
CN110286416B 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) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112861065A (zh) * | 2021-02-10 | 2021-05-28 | 湖北省地震局(中国地震局地震研究所) | 一种自适应的震源区重力密度变化实时反演计算方法 |
CN113311494A (zh) * | 2021-05-26 | 2021-08-27 | 中南大学 | 一种卫星重力场反演方法 |
CN113486591A (zh) * | 2021-07-13 | 2021-10-08 | 吉林大学 | 一种卷积神经网络结果的重力多参量数据密度加权反演方法 |
CN113514899A (zh) * | 2021-07-12 | 2021-10-19 | 吉林大学 | 一种重力梯度的自适应剖分反演方法 |
Citations (13)
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 | 中国石油集团东方地球物理勘探有限责任公司 | 基于标准格架的快速三维重力、磁力物性反演的方法 |
US20110313737A1 (en) * | 2010-06-21 | 2011-12-22 | Schlumberger Technology Corporation | Mehtod of determining parameter from sparse measurement data |
CN103913774A (zh) * | 2014-04-02 | 2014-07-09 | 西南石油大学 | 基于微地震事件的储层地质力学参数反演方法 |
CN104866653A (zh) * | 2015-04-29 | 2015-08-26 | 中国地质科学院矿产资源研究所 | 一种获取地下三维密度结构的方法 |
WO2017131848A1 (en) * | 2016-01-26 | 2017-08-03 | Saudi Arabian Oil Company | Gravity modeling a rifted continental margin |
CN107491411A (zh) * | 2017-06-23 | 2017-12-19 | 中国海洋大学 | 基于n阶多项式密度函数的重力异常反演方法 |
CN107561592A (zh) * | 2017-09-11 | 2018-01-09 | 中南大学 | 一种密度为多项式的任意多面体的重力场正演方法 |
CN107748399A (zh) * | 2017-09-12 | 2018-03-02 | 中国石油化工股份有限公司 | 利用重力界面反演识别山前带深部构造层方法 |
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 | 吉林大学 | 一种不同高度重力数据和井中重力数据联合反演方法 |
-
2019
- 2019-05-13 CN CN201910393070.0A patent/CN110286416B/zh not_active Expired - Fee Related
Patent Citations (13)
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 | 中国石油集团东方地球物理勘探有限责任公司 | 基于标准格架的快速三维重力、磁力物性反演的方法 |
US20110313737A1 (en) * | 2010-06-21 | 2011-12-22 | Schlumberger Technology Corporation | Mehtod of determining parameter from sparse measurement data |
CN103913774A (zh) * | 2014-04-02 | 2014-07-09 | 西南石油大学 | 基于微地震事件的储层地质力学参数反演方法 |
CN104866653A (zh) * | 2015-04-29 | 2015-08-26 | 中国地质科学院矿产资源研究所 | 一种获取地下三维密度结构的方法 |
WO2017131848A1 (en) * | 2016-01-26 | 2017-08-03 | Saudi Arabian Oil Company | Gravity modeling a rifted continental margin |
CN107491411A (zh) * | 2017-06-23 | 2017-12-19 | 中国海洋大学 | 基于n阶多项式密度函数的重力异常反演方法 |
CN107561592A (zh) * | 2017-09-11 | 2018-01-09 | 中南大学 | 一种密度为多项式的任意多面体的重力场正演方法 |
CN107748399A (zh) * | 2017-09-12 | 2018-03-02 | 中国石油化工股份有限公司 | 利用重力界面反演识别山前带深部构造层方法 |
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 (6)
Title |
---|
LAST B J等: ""Compact gravity inversion"", 《GEOPHYSICS》 * |
YUNUS L E.: ""2D focusing inversion of gravity data with the use of parameter variation as a stopping criterion"", 《JOURNAL OF THE BALKAN GEOPHYSICAL SOCIETY》 * |
杨积忠 等: ""基于Born敏感核函数的速度、密度双参数全波形反演"", 《地球物理学报》 * |
栗学磊 等: ""重力资料聚焦反演成像改进及其应用"", 《世界地质》 * |
王泰涵: ""重力及其梯度张量数据高精度成像和密度反演方法"", 《中国博士学位论文全文数据库 基础科学辑》 * |
苏超: ""基于重力数据的地质体埋深及物性三维反演方法研究"", 《中国硕士学位论文全文数据库 基础科学辑》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112861065A (zh) * | 2021-02-10 | 2021-05-28 | 湖北省地震局(中国地震局地震研究所) | 一种自适应的震源区重力密度变化实时反演计算方法 |
CN112861065B (zh) * | 2021-02-10 | 2022-07-08 | 湖北省地震局(中国地震局地震研究所) | 一种自适应的震源区重力密度变化实时反演计算方法 |
CN113311494A (zh) * | 2021-05-26 | 2021-08-27 | 中南大学 | 一种卫星重力场反演方法 |
CN113311494B (zh) * | 2021-05-26 | 2022-04-26 | 中南大学 | 一种卫星重力场反演方法 |
CN113514899A (zh) * | 2021-07-12 | 2021-10-19 | 吉林大学 | 一种重力梯度的自适应剖分反演方法 |
CN113486591A (zh) * | 2021-07-13 | 2021-10-08 | 吉林大学 | 一种卷积神经网络结果的重力多参量数据密度加权反演方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110286416B (zh) | 2020-12-22 |
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 | |
Dannemiller et al. | A new method for determination of magnetization direction | |
Huang et al. | Bidimensional empirical mode decomposition (BEMD) for extraction of gravity anomalies associated with gold mineralization in the Tongshi gold field, Western Shandong Uplifted Block, Eastern China | |
CN105426619B (zh) | 一种确定盾构隧道纵向弯曲刚度的分析方法 | |
CN113488117B (zh) | 一种具有深度学习能力的深部金矿床成矿找矿方法 | |
Guo et al. | Generation and verification of three-dimensional network of fractured rock masses stochastic discontinuities based on digitalization | |
CN112949134B (zh) | 基于非结构有限元方法的地-井瞬变电磁反演方法 | |
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 | |
CN107748399B (zh) | 利用重力界面反演识别山前带深部构造层方法 | |
Surrette et al. | Regional evaluation of hydraulic properties in variably fractured rock using a hydrostructural domain approach | |
KR20160024232A (ko) | 지반정보에 기초한 3차원 공간 모델링 방법 | |
Koike et al. | Incorporation of fracture directions into 3D geostatistical methods for a rock fracture system | |
CN110133644B (zh) | 基于插值尺度函数法的探地雷达三维正演方法 | |
CN108646298B (zh) | 多井多参数约束下隧道波速反演成像及围岩分级的方法 | |
CN105738952A (zh) | 一种水平井区储层岩石相建模方法 | |
CN115327663A (zh) | 深部矿产资源勘查用空-地-井立体地球物理探测方法 | |
CN115292890A (zh) | 基于多源辅助数据开发的场地土壤污染物浓度三维空间预测方法 | |
Du et al. | Distance-gradient-based variogram and Kriging to evaluate cobalt-rich crust deposits on seamounts | |
Amigun et al. | Euler 3-D deconvolution of analytical signal of magnetic anomalies over iron ore deposit in Okene, Nigeria | |
CN109490978A (zh) | 一种起伏地层的频率域快速高精度正演方法 | |
CN112596113A (zh) | 一种基于重力不同阶梯度特征值交点的场源位置识别方法 | |
CN112882124B (zh) | 一种三维铀成矿有利部位的圈定方法 | |
CN105242317B (zh) | 一种纵波速度的确定方法及装置 | |
CN114280688A (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 |