CN112346139B - 一种重力数据多层等效源延拓与数据转换方法 - Google Patents
一种重力数据多层等效源延拓与数据转换方法 Download PDFInfo
- Publication number
- CN112346139B CN112346139B CN202011104563.7A CN202011104563A CN112346139B CN 112346139 B CN112346139 B CN 112346139B CN 202011104563 A CN202011104563 A CN 202011104563A CN 112346139 B CN112346139 B CN 112346139B
- Authority
- CN
- China
- Prior art keywords
- gravity
- data
- gravity field
- equivalent source
- 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
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
- G01V7/02—Details
- G01V7/06—Analysis or interpretation of gravimetric records
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
本发明提供一种重力数据等效源延拓与数据转换方法,包括:输入起伏观测曲面上的重力场数据d0,根据观测面所在区域的地形信息,建立地形起伏曲面;根据起伏观测曲面以及设定的反演最大深度,确定网格剖分范围,并对所述范围进行网格剖分以确定模型空间;根据重力背景场,基于模型空间对重力场数据d0进行带约束项及规整化项的积分方程三维反演计算,得到重力异常体的等效源模型;根据延拓面的位置确定灵敏度矩阵G′,然后利用等效源模型进行基于积分方程的重力场正演计算,得到异常体产生的延拓后的重力场数据Us;对Us进行梯度张量转换。本发明能对重力异常体进行自适应且快速高效准确地生成所需要延拓与数据转换后的数据。
Description
技术领域
本发明涉及重力勘探技术领域,尤其涉及一种重力数据多层等效源延拓与数据转换方法。
背景技术
在重力探测中,人们通常对重力场的铅垂一次导数进行探测,但在实际的数据解释中,又往往需要将这些数据转换成所需要的参数和类型,比如不同观测高度(延拓)、梯度张量数据的换算等,重力场数据类型转换的主要挑战在于航测的起伏测量曲面以及数据的不规则测量位置(非网格化测量),传统的数据类型转换以及延拓计算难以处理这些数据,而采用传统单层、双层、三层等效源的数据转换方法,处理结果精度较低,无法满足应用需求。
现有文献1“Dampney,C.N.G.THE EQUIVALENT SOURCE TECHNIQUE[J].geophysics,1969,34(1):39.”首次提出了等效源方法,利用虚拟场源模拟实测异常,可用于位场数据的空间延拓(包括曲面延拓)、梯度计算以及分量转换等;选择单层等效源且将其布置于近地表是等效源方法的主要特点,比如,文献2“Lyrio J C S DO.2011.Equivalent source:A natural choice for gridding scatter gravity data[J].In 12th International Congress of the Brazilian Geophysical Society,pp.661-665,doi:10.1190/sbgf2011-136.”使用了单层等效源进行重力场数据的分析,文献3“Martinez C,Li Y G.2016.Denoising of gravity gradient data using anequivalent source technique[J].Geophysics,81(4):G67-G79,doi:10.1190/geo2015-0379.1.”使用单层等效源对重力梯度张量数据进行了去噪;为了在保证计算效率的同时高精度重构位场,多层等效源方法是一个相对合理的选择,文献4“李端;陈超;杜劲松;梁青;黎海龙;基于多层等效源的重力数据处理方法[C]//2017中国地球科学联合学术年会.0.”提出了一种将地下等效源分为浅、中、深三层来实现重力场数据的转换,文献5“陈涛,张贵宾.利用重力异常计算重力梯度的等效源技术[J].地球物理学进展.”提出了一种多层等效源进行重力梯度计算的方法,其在垂向方向上将剖分区域分为三个子区域。
这些技术的应用使得人们可以通过等效源技术来进行地面重力数据的化极与数据类型转换运算,但存在以下问题:1)上述方法构建的等效源层数通常小于等于三层,且网格剖分不连续;2)需单独对每一个等效源的深度位置进行估算,然后单独放置;3)传统延拓方法对于数据的延拓距离有限制(一般小于6倍的数据点距),大于限制距离的延拓计算精度较差。
发明内容
有鉴于此,本发明提供了一种重力数据多层等效源延拓与数据转换方法,采用连续的结构化非均匀网格剖分,并基于积分方程正反演理论框架,引入深度规整化因子,在反演过程中直接确定多层等效源的深度和分布,包括以下步骤:
S1、输入起伏观测曲面上的重力场数据d0,并根据起伏观测曲面所在区域的地形高度信息,建立地形起伏曲面;
S2、根据起伏观测曲面以及设定的反演最大深度,确定网格剖分的空间范围,并根据S1建立的地形起伏曲面对所述空间范围进行连续的结构化非均匀网格剖分,进一步确定反演模型求解空间;
S3、根据重力场的背景场,基于S2反演模型求解空间对重力场数据d0进行带深度规整化因子、正值约束项以及规整化项的积分方程三维反演计算,得到重力异常体的多层等效源模型;
S4、设定延拓面,并根据延拓面得到灵敏度矩阵G′,然后利用S3得到的多层等效源模型进行基于积分方程的重力场正演计算,得到重力异常体产生的延拓后的重力场数据Us;
S5、基于S4所求的重力场数据Us,利用重力场梯度张量转换公式,进行重力场数据的转换。
进一步地,S2中,所述网格剖分的空间范围包括上顶面和下底面,其中,所述上顶面为地形起伏曲面的最大高度所确定的平面,所述下底面为设定的反演最大深度所确定的平面。
进一步地,S2中,根据地形起伏曲面的最低点对网格剖分的空间范围进行划分,其中,对所述最低点以上的空间范围进行均匀网格剖分得到精细网格,对所述最低点以下的空间范围进行非均匀网格剖分得到扩展网格;所述地形起伏曲面至下底面之间的空间范围构成反演模型求解空间。
进一步地,所述扩展网格的垂直边以精细网格的垂直边的α1倍速度增长,且设定其最大增速为α2,其中,α2>α1>1。
进一步地,S3中所述多层等效源模型的模型深度面的数目大于3层。
进一步地,S3中,积分方程三维反演计算的目标函数为:
其中,φ表示优化目标即误差,m表示待求解的多层等效源模型的密度矩阵,m≥0;F(·)表示基于积分方程的重力场正演操作,G表示灵敏度矩阵,与起伏观测曲面的位置以及网格剖分相关;β表示根据实际需求添加的规整化因子;mref表示参考模型的密度矩阵,Wr表示深度规整化因子。
进一步地,所述深度规整化因子Wr为:
其中,z表示等效源到地形起伏曲面的距离,z0表示起伏观测曲面的高度,r表示深度系数。
进一步地,S4得到重力异常体产生的延拓后的重力场数据Us=F(m,G′)=G′m,其中,m为密度矩阵,F(·)表示基于积分方程的重力场正演操作,G′为由延拓面确定的灵敏度矩阵。
进一步地,S5中重力场数据的转换如下:
本发明提供的技术方案带来的有益效果是:本发明提出的技术方案能对复杂环境中的地下重力异常体产生的重力场进行基于积分方程的非均匀网格的多层等效源延拓与数据类型转换运算,采用连续网格剖分,数目通常在3层以上,精度更高;同时加入深度规整化因子,不需要单独估计等效层的深度和范围,能对重力异常体进行自适应且快速高效准确地生成所需要的延拓与数据类型转换后的数据,具有更高的稳定性和精度。
附图说明
图1是本发明一种重力数据多层等效源延拓与数据转换方法的流程图;
图2是本发明一种重力数据多层等效源延拓与数据转换方法的非均匀网格剖分的示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地描述。
请参考图1,本发明的实施例提供了一种重力数据多层等效源延拓与数据转换方法,包括以下步骤:
S1、输入起伏观测曲面上的重力场数据d0,并根据观测面所在区域的地形高度信息,建立地形起伏曲面;所述重力场数据d0可以是重力异常铅垂数据或重力梯度张量数据,本实施例以重力异常铅垂数据为例;
S2、根据地形起伏曲面以及设定的反演最大深度,确定网格剖分的空间范围,并根据地形起伏曲面对所述空间范围进行连续的结构化非均匀网格剖分,进一步确定反演模型求解空间;
网格剖分的空间范围包括上顶面和下底面,其中,根据起伏观测曲面的最大高度确定上顶面,然后基于现有探测技术或实际经验估计出重力异常体可能存在的最大深度(即反演最大深度),并据此确定下底面;
在确定网格剖分的空间范围之后,以地形起伏曲面的最低点对所述空间范围进行划分,对最低点以上的空间范围进行均匀网格剖分得到精细网格,对最低点以下的空间范围进行非均匀的扩展网格剖分得到扩展网格,所述地形起伏曲面至下底面之间的空间范围构成反演模型求解空间;优选地,若所述精细网格的垂直边为1长度单位(该长度单位的具体数值可根据观测区域空间大小进行设定,比如100m为1长度单位),则所述扩展网格的垂直边以精细网格垂直边的1.2倍速增长,且设定最大增速为1.5倍,由此在保证一定反演精度的基础上降低计算量。请参考图2,其为非均匀网格剖分的结果示意图,图中坐标轴分别为东向(Easting)、北向(Northing)、垂向(Depth)。
S3、根据重力场的背景场,基于S2反演模型求解空间对重力场数据d0进行带深度规整化因子、正值约束项以及规整化项的积分方程三维反演计算,得到重力异常体的多层等效源模型,所述多层等效源模型具体是指反演模型求解空间中包含多个模型深度面,基于上述方案,可求解得到的模型深度面数目通常大于3层。
基于积分方程的三维反演计算的目标函数为:
式中,φ表示优化目标,即误差,表示目标函数的数值约束,表示目标函数的模型约束,m表示待求解的多层等效源模型的密度矩阵,考虑到物体的物理性质导致密度必须为正值,故进行正值约束,即m≥0;F(·)表示基于积分方程的重力场正演操作,G表示灵敏度矩阵,与起伏观测曲面的位置以及网格剖分相关;β表示根据实际需求添加的规整化因子,若不需要可令β=1;mref表示参考模型的密度矩阵,Wr表示深度规整化因子,计算公式如下:
其中,z表示等效源到地形起伏曲面的距离,z0表示起伏观测曲面的高度,r表示深度系数,本实施例中取r=3。
S4、设定延拓面,并根据延拓面得到灵敏度矩阵G′,然后利用步骤S3得到的多层等效源模型进行基于积分方程的重力场正演计算,得到重力异常体产生的延拓后的重力场数据Us,即:
Us=F(m,G′)=G′m,
其中,m表示S3输出的密度矩阵;
通过下述公式从重力场数据Us转换为重力异常数据:
其中,d为铅垂方向的重力异常数据;
S5、基于所求的重力场数据Us,利用下述的重力场梯度张量转换公式,进行重力场数据的转换:
上面结合附图对本发明的实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可做出很多形式,这些均属于本发明的保护之内。
Claims (2)
1.一种重力数据多层等效源延拓与数据转换方法,其特征在于,包括以下步骤:
S1、输入起伏观测曲面上的重力场数据d0,并根据起伏观测曲面所在区域的地形高度信息,建立地形起伏曲面;
S2、根据起伏观测曲面以及设定的反演最大深度,确定网格剖分的空间范围,并根据S1建立的地形起伏曲面对所述空间范围进行连续的结构化非均匀网格剖分,进一步确定反演模型求解空间;根据地形起伏曲面的最低点对网格剖分的空间范围进行划分,其中,对所述最低点以上的空间范围进行均匀网格剖分得到精细网格,对所述最低点以下的空间范围进行非均匀网格剖分得到扩展网格;所述地形起伏曲面至下底面之间的空间范围构成反演模型求解空间;
S3、根据重力场的背景场,基于S2反演模型求解空间对重力场数据d0进行带深度规整化因子、正值约束项以及规整化项的积分方程三维反演计算,得到重力异常体的多层等效源模型;所述多层等效源模型的模型深度面的数目大于3层;
积分方程三维反演计算的目标函数为:
其中,φ表示优化目标,即误差,m表示待求解的多层等效源模型的密度矩阵,m≥0;F(·)表示基于积分方程的重力场正演操作,G表示灵敏度矩阵,与起伏观测曲面的位置以及网格剖分相关;β表示根据实际需求添加的规整化因子;mref表示参考模型的密度矩阵,Wr表示深度规整化因子;所述深度规整化因子Wr为:
其中,z表示等效源到地形起伏曲面的距离,z0表示起伏观测曲面的高度,r表示深度系数;
S4、设定延拓面,并根据延拓面得到灵敏度矩阵G′,然后利用S3得到的多层等效源模型进行基于积分方程的重力场正演计算,得到重力异常体产生的延拓后的重力场数据Us;重力场数据Us=F(m,G′)=G′m,其中,m为密度矩阵,F(·)表示基于积分方程的重力场正演操作,G′为由延拓面确定的灵敏度矩阵;
S5、基于S4所求的重力场数据Us,利用重力场梯度张量转换公式,进行重力场数据的转换;重力场数据的转换如下:
2.根据权利要求1所述的一种重力数据多层等效源延拓与数据转换方法,其特征在于,S2中,所述网格剖分的空间范围包括上顶面和下底面,其中,所述上顶面为地形起伏曲面的最大高度所确定的平面,所述下底面为设定的反演最大深度所确定的平面。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011104563.7A CN112346139B (zh) | 2020-10-15 | 2020-10-15 | 一种重力数据多层等效源延拓与数据转换方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011104563.7A CN112346139B (zh) | 2020-10-15 | 2020-10-15 | 一种重力数据多层等效源延拓与数据转换方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112346139A CN112346139A (zh) | 2021-02-09 |
CN112346139B true CN112346139B (zh) | 2021-12-21 |
Family
ID=74360903
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011104563.7A Active CN112346139B (zh) | 2020-10-15 | 2020-10-15 | 一种重力数据多层等效源延拓与数据转换方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112346139B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113627051B (zh) * | 2021-07-23 | 2024-01-30 | 中国地质科学院地球物理地球化学勘查研究所 | 一种重力异常场分离方法、系统、存储介质和电子设备 |
CN114167511B (zh) * | 2021-11-26 | 2023-08-11 | 兰州大学 | 一种基于连分式展开向下延拓的位场数据快速反演方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10242126B2 (en) * | 2012-01-06 | 2019-03-26 | Technoimaging, Llc | Method of simultaneous imaging of different physical properties using joint inversion of multiple datasets |
CN104316972A (zh) * | 2014-10-16 | 2015-01-28 | 中国海洋石油总公司 | 一种磁源重力视强度反演成像方法 |
CN105549106B (zh) * | 2016-01-07 | 2018-06-08 | 中国科学院地质与地球物理研究所 | 一种重力多界面反演方法 |
CN110045432B (zh) * | 2018-12-05 | 2020-06-30 | 中南大学 | 基于3d-glq的球坐标系下重力场正演方法及三维反演方法 |
CN109375280B (zh) * | 2018-12-10 | 2020-03-31 | 中南大学 | 一种球坐标系下重力场快速高精度正演方法 |
CN111323830B (zh) * | 2020-01-14 | 2021-06-25 | 东华理工大学 | 一种基于大地电磁和直流电阻率数据的联合反演方法 |
-
2020
- 2020-10-15 CN CN202011104563.7A patent/CN112346139B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN112346139A (zh) | 2021-02-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112363236B (zh) | 一种基于pde的重力场数据等效源延拓与数据类型转换方法 | |
CN111856599B (zh) | 一种基于pde的磁测数据等效源化极与类型转换方法 | |
CN112346139B (zh) | 一种重力数据多层等效源延拓与数据转换方法 | |
CN112147709B (zh) | 一种基于部分光滑约束的重力梯度数据三维反演方法 | |
CN111856598B (zh) | 一种磁测数据多层等效源上延拓与下延拓方法 | |
CA2713019A1 (en) | Gridless geological modeling | |
AU2013299551A1 (en) | Method for improving the accuracy of rock property values derived from digital images | |
AU2012336262B2 (en) | Method of generating and combining multiple horizons to determine a seismic horizon and its uncertainty | |
CN109884710B (zh) | 针对激发井深设计的微测井层析成像方法 | |
GB2474740A (en) | Gridless geological modeling of a structural framework | |
CN112116708B (zh) | 一种获取三维地质实体模型的方法及系统 | |
CN112150582B (zh) | 一种面向多模态数据的地质剖面图近似表达方法 | |
CN109085643A (zh) | 早至波的分步联合反演方法 | |
CN114943178A (zh) | 一种三维地质模型建模方法、装置及计算机设备 | |
CN109444973B (zh) | 一种球坐标系下重力正演加速方法 | |
CN111859251B (zh) | 一种基于pde的磁测数据等效源上延拓与下延拓方法 | |
EP3281044B1 (en) | Method for estimating elastic parameters of subsoil | |
CN111045098B (zh) | 一种拾取地下深部构造信息的方法 | |
CN111880236A (zh) | 一种构建多层等效源模型计算化极与数据类型转换的方法 | |
CN114200541B (zh) | 一种基于余弦点积梯度约束的三维重磁联合反演方法 | |
CN113591030B (zh) | 基于多gpu的重力梯度数据灵敏度矩阵压缩及调用方法 | |
CN113970732A (zh) | 一种三维频率域探地雷达双参数同步反演方法 | |
CN111508075B (zh) | 一种闭孔泡沫铝的三维真实有限元模型建模方法 | |
Yin et al. | A fast 3D gravity forward algorithm based on circular convolution | |
CN112612935A (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 |