CN104316972A - 一种磁源重力视强度反演成像方法 - Google Patents
一种磁源重力视强度反演成像方法 Download PDFInfo
- Publication number
- CN104316972A CN104316972A CN201410549240.7A CN201410549240A CN104316972A CN 104316972 A CN104316972 A CN 104316972A CN 201410549240 A CN201410549240 A CN 201410549240A CN 104316972 A CN104316972 A CN 104316972A
- Authority
- CN
- China
- Prior art keywords
- magnetic
- centerdot
- gravity anomaly
- delta
- magnetic body
- 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
Landscapes
- Measuring Magnetic Variables (AREA)
Abstract
本发明涉及一种磁源重力视强度反演成像方法,其包括以下步骤:将磁异常转换成空间域的磁源重力异常;选取不同的切割半径,采用插值切割法对空间域的磁源重力异常进行不同深度空间层的磁源重力异常分离,获取各空间层的剩余磁源重力异常;采用迭代法将各空间层的剩余磁源重力异常向下延拓到各空间层的磁性体的顶界面上,获取磁性体顶界面的磁源重力异常;对磁性体顶界面的磁源重力异常的视强度参数进行反演,获取各空间层的视强度参数。本发明具有反演精度高且能够降低反演多解性的优点,本发明可以广泛应用于深部矿产勘探和深部地质构造研究中。
Description
技术领域
本发明涉及一种反演成像方法,特别是关于应用地球物理学领域中的一种磁源重力视强度反演成像方法。
背景技术
目前,国内外用于磁场物性(磁化强度)反演的方法分为两类:线性反演和非线性反演。具体的算法大致分为广义最小二乘法(Gauss法)、最速下降法、阻尼最小二乘法、奇异值分解法、人工神经网络法BP算法、遗传算法、模拟退火算法和基于Radon变换、小波多尺度边缘的算法等。上述反演方法在理论创新及小面积实际资料反演应用中都取得了一定的成绩,但在用于大面积的三维物性反演的实际生产当中没有发挥到应有的作用。其主要存在问题:1)物性参数与几何参数全局反演时,需要解大型线性代数方程组,计算时间长;2)求解方程组数学病态十分严重,解的稳定性很差;3)磁源偶极正负异常相互叠加效应导致异常反演解精度差、分辨率低。
发明内容
针对上述问题,本发明的目的是提供一种反演精度高且能够降低反演多解性的磁源重力视强度反演成像方法。
为实现上述目的,本发明采取以下技术方案:一种磁源重力视强度反演成像方法,其包括以下步骤:1)将磁异常△T(x,y)转换成空间域的磁源重力异常△g(x,y);2)选取不同的切割半径,采用插值切割法对空间域的磁源重力异常△g(x,y)进行不同深度空间层的磁源重力异常分离,获取各空间层的剩余磁源重力异常△gj(x,y,0);3)采用迭代法将各空间层的剩余磁源重力异常△gj(x,y,0)向下延拓到各空间层的磁性体的顶界面上,获取磁性体顶界面的磁源重力异常△g(x,y,hj);4)对磁性体顶界面的磁源重力异常△g(x,y,hj)的视强度参数进行反演,获取各空间层的视强度参数。
所述步骤1)中,将磁异常△T(x,y)转换成空间域的磁源重力异常△g(x,y),其包括以下步骤:(I)将地面观测到的磁异常△T(x,y)转换成频率域的磁异常即
式中,u和v分别表示x方向和y方向的波数;(II)将频率域的磁异常转换成化极磁异常即
式中,i表示虚分量,(P0,Q0,R0)表示磁性体所在地层的方向余弦,(P1,Q1,R1)表示磁性体在地表的方向余弦,(P2,Q2,R2)表示磁性体在磁倾角为90°的磁化方向下的方向余弦;(III)将化极磁异常的偶极磁源作类重力的单极源处理,即对沿磁性体南北极方向积分,得到频率域的磁源重力异常即
式中,M表示磁化强度矢量,G表示万有引力常量,ρ表示磁性体的剩余密度;(IV)对频率域的磁源重力异常进行反傅里叶变换,得到空间域的磁源重力异常△g(x,y),即
所述步骤2)中,获取各空间层的剩余磁源重力异常△gj(x,y,0),其具体包括以下步骤:(I)选取不同的切割半径r1,r2,r3…rn,其中,n表示切割半径的个数,采用插值切割方法对磁源重力异常△g(x,y)进行不同深度空间层的磁源重力异常分离,获取不同切割半径下的浅源异常:
L1(x,y),L2(x,y),L3(x,y)…Ln(x,y),
(II)根据不同切割半径下的浅源异常,得到各空间层的剩余磁源重力异常:
式中,△gj(x,y,0)表示第j层磁性体在地面反映的磁源重力异常,即剩余磁源重力异常。
所述步骤3)中,采用迭代法向下延拓,获取磁性体顶界面的磁源重力异常△g(x,y,hj),其具体包括以下步骤:(I)将地面ΓA上的异常值△g(x,y,0)中第s点的异常值GAs放到平面ΓB的垂直投影点s上,作为平面ΓB上的异常值即并通过傅里叶变换得到ΓB上异常值的频谱(II)当地面ΓA与平面ΓB之间没有场源时,重力异常满足拉普拉斯方程,采用向上延拓公式:
将平面ΓB上的异常值的频谱代入上式求出地面ΓA上s点的异常值即:
式中,表示地面ΓA(z=0)上的磁源重力异常值△g(x,y,0)的傅里叶变换;F-1表示反傅里叶变换,h表示延拓高度;(III)利用GAs与的差对平面ΓB上的异常值进行修正,得到新的异常值即:
式中,l表示步长;(IV)重复步骤(II)和步骤(III),得到如下迭代公式:
当 时,由式(16)可知 ε表示一个给定的趋近于零的数,则迭代次数n取20~50次;(V)采用与步骤(I)~(IV)相同的方法,计算得到平面ΓB上全部点的磁源重力异常值并利用向上延拓公式:
作反傅里叶变换得到△g(x,y,0)向下延拓hj深度的磁源重力异常△g(x,y,hj)。
所述步骤4)中,获取各空间层的视强度参数,其具体包括以下步骤:(I)将引起磁源重力异常△g(x,y,hj)的非规则的磁性体分割为M*N个大小规则、质心位置不同的细长垂直棱柱体;设棱柱体的质心点坐标为(x0k,y0k),其剩余密度为ρk,棱柱体的底面边长分别为a和b,棱柱体高度为dh,则在频率域中各棱柱体在顶界面(x,y)点产生的重力异常△gk(x,y)的频谱为:
式中,e表示自然常数,k表示棱柱体序号,k=1,2,3,…,M*N;(II)将各棱柱体在顶界面(x,y)点产生的重力异常的频谱进行叠加,得到整个深度层在顶界面(x,y)点产生的重力异常△g(x,y)的频谱为:
令
则磁源重力异常在顶界面的频谱响应为:
则频域中空间层的剩余磁源重力异常顶界面的视密度为:
(III)对频域中空间层的剩余磁源重力异常顶界面的视密度进行反傅里叶变换,得到该深度层顶面各点(x,y)的剩余密度值为:
(IV)将剩余密度值ρ(x,y)代入式
计算得到磁异常的视磁化率κ(x,y);将计算得到的磁异常的视磁化率κ(x,y)代入式
M=κ(x,y)*T0,
计算得到磁异常的视磁化强度M;式中,κ(x,y)为磁异常的视磁化率,M为磁异常的视磁化强度,ρ(x,y)为磁源重力的视密度,T0为正常磁场强度,κ'和ρ'分别是泊松公式转换时需要的已知磁化率值和密度值。
本发明由于采取以上技术方案,其具有以下优点:1、现有反演方法是基于空间域的优化算法,通常是物理参数及几何参数全局反演,反演过程复杂,计算精度低;与现有反演方法相比,本发明的反演方法只对物性强度参数进行反演,大大提高了反演精度,降低了反演的多解性。2、现有反演方法需求解大型线性方程组,计算机内存消耗巨大,反演时间长,因此不能进行大面积快速的三维成像反演;本发明的反演方法是在频率域进行的,运算速度快,无需求解大型线性方程组,计算机内存消耗小,因此能进行大范围、快速的三维视磁化强度反演成像。3、现有技术中偶极磁源采用传统的类重力反演方法(即用重力异常反演方法直接反演磁异常)时,会反演出两个不同的正负磁化强度参数;而本发明的磁源重力视强度反演成像方法直观,且能够反映磁性体的中心位置。4、现有反演方法大多是针对单一模型的,而对于多源模型反演,则因磁异常相互叠加,异常的等效性使得反演更加困难,反演精度差;而本发明的磁源重力视强度反演成像方法采用多源偶极磁源,能进行类似重力异常的单极叠加异常分离,并能进行基于磁源外(含磁源顶界面)任意深度空间磁异常信息获取的高分辨率视磁化强度反演求解。5、本发明的反演方法能够建立三维空间的视磁化强度属性分布。基于以上优点,本发明可以广泛应用于深部矿产勘探和深部地质构造研究中。
附图说明
图1是本发明磁源重力视强度反演成像方法的流程图;
图2是采用迭代法向下延拓的示意图;
图3是理论模型空间分布及地面投影图;其中,图(a)是理论模型空间分布图,图(b)是理论模型的地面投影图;
图4是理论模型地面磁异常示意图;
图5是地面磁异常化向地磁极异常示意图;
图6是模型地面磁源重力异常示意图;
图7是0~200米空间层剩余磁源重力异常图;
图8是200~500米空间层剩余磁源重力异常图;
图9是200米深顶界面磁源重力异常图;
图10是500米深顶界面磁源重力异常图;
图11是200米深顶界面磁源重力视密度反演图;
图12是500米深顶界面磁源重力视密度反演图;
图13是200米深顶界面磁源重力视磁化强度反演图;
图14是500米深顶界面磁源重力视磁化强度反演图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1所示,本发明磁源重力视强度反演成像方法包括以下步骤:
1)将磁异常△T(x,y)转换成空间域的磁源重力异常△g(x,y);
地下地质体间存在密度差异或磁性差异时,地面上可观测到其相应的重力异常△g(x,y)或磁异常△T(x,y),由于重磁场满足拉普拉斯位场理论,因此数学上可做相互转换,即重力异常可以用梯度的形式换算成重力源磁异常、磁异常也可以用积分的形式转换成磁源重力异常。下面具体说明磁异常转换成磁源重力异常的实现过程:
(I)将地面观测到的磁力异常△T(x,y)转换成频率域的磁异常
对地面观测到的磁异常△T(x,y)进行傅里叶变换,得到地面磁异常△T(x,y)的频谱即
式中,u和v分别表示x方向和y方向的波数。
(II)将频率域的磁异常换算成具有重力方向的磁异常,即转换成化极磁异常
设定地表实测磁力异常为△T(x,y)的磁性体其初始磁倾角为I1,磁偏角为D1,该磁性体在地表的方向余弦为(P1,Q1,R1);该磁性体所在地层的地磁倾角为I0,地磁偏角为D0,该磁性体所在地层的方向余弦为(P0,Q0,R0);该磁性体的磁化方向从磁倾角为I1转换到磁倾角为I2时,磁偏角为D2,该磁性体在该磁化方向下的方向余弦为(P2,Q2,R2),则将磁化方向从磁倾角为I1转换成到磁倾角为I2时,磁化方向转化公式为:
式中,i表示虚分量,即i2=-1。该磁性体在地表的方向余弦(P1,Q1,R1)为:
该磁性体所在地层的方向余弦(P0,Q0,R0)为:
该磁性体在磁倾角为I2的磁化方向下的方向余弦(P2,Q2,R2)为:
当该磁性体的磁化方向从磁倾角为I1转换到垂直磁化方向(化向地磁极)时,即I2=90°,将I2=90°代入式(3)中,则该磁性体在磁倾角为I2的磁化方向下的方向余弦(P2,Q2,R2)为:
将式(4)带入式(2)中,得到化向地磁极的频谱为:
(III)将化极磁异常的偶极磁源作类重力的单极源处理,即对沿磁性体南北极方向积分,得到频率域的磁源重力异常
由于重磁位梯度值分别对应梯度方向的磁场强度和重力场强度值,即
且磁位U与引力位V之间的关系满足泊松公式,即
式中,M表示磁化强度矢量,G表示万有引力常量,ρ表示磁性体的剩余密度。
由式(5)~式(7)计算得到频率域的磁源重力异常为:
(IV)对频率域的磁源重力异常进行反傅里叶变换,得到空间域的磁源重力异常△g(x,y),即
2)选取不同的切割半径,采用插值切割法对空间域的磁源重力异常△g(x,y)进行不同深度空间层的磁源重力异常分离,获取各空间层的剩余磁源重力异常△gj(x,y,0),其具体包括以下步骤:
(I)选取不同的切割半径r1,r2,r3…rn,其中,n表示切割半径的个数,采用插值切割方法对磁源重力异常△g(x,y)进行不同深度空间层的磁源重力异常分离,获取不同切割半径下的浅源异常:
L1(x,y),L2(x,y),L3(x,y)…Ln(x,y) (10)
(II)根据不同切割半径下的浅源异常,得到各空间层的剩余磁源重力异常:
式中,△gj(x,y,0)表示第j(j=0,1,2,3…n)层磁性体在地面反映的磁源重力异常,即剩余磁源重力异常。
3)采用迭代法将各空间层的剩余磁源重力异常△gj(x,y,0)向下延拓到各空间层的磁性体的顶界面上,获取磁性体顶界面的磁源重力异常△g(x,y,hj),其中,△g(x,y,hj)表示将j层剩余磁源重力异常△gj(x,y,0)向下延拓至地面以下hj深度,获取磁性体顶界面的磁源重力异常,深度hj为:
采用迭代法向下延拓,获取磁性体顶界面的磁源重力异常△g(x,y,hj)的具体步骤包括:
(I)如图2所示,将地面ΓA上的异常值△g(x,y,0)中第s点的异常值GAs放到平面ΓB的垂直投影点s上,作为平面ΓB上的异常值即并通过傅里叶变换得到ΓB上异常值的频谱
(II)当地面ΓA与平面ΓB之间没有场源时,重力异常满足拉普拉斯方程,采用向上延拓公式:
将平面ΓB上的异常值的频谱代入上式求出地面ΓA上s点的异常值即:
式(13)中,表示地面ΓA(z=0)上的磁源重力异常值△g(x,y,0)的傅里叶变换;F-1表示反傅里叶变换,h表示延拓高度。
(III)利用GAs与的差对平面ΓB上的异常值进行修正,得到新的异常值即:
式中,l表示步长,一般取值为1。
(IV)重复步骤(II)和步骤(III),得到如下迭代公式:
当 时,由式(16)可知 ε表示一个给定的趋近于零的数,则
一般迭代次数n取20~50次。
(V)采用与步骤(I)~(IV)相同的方法,计算得到平面ΓB上全部点的磁源重力异常值并利用向上延拓公式(13)得到△g(x,y,0)向下延拓hj深度的磁源重力异常△g(x,y,hj)。
4)对顶界面磁源重力异常△g(x,y,hj)的视强度参数进行反演,获取各空间层的视强度参数,其具体包括以下步骤:
(I)将引起磁源重力异常△g(x,y,hj)的非规则的磁性体分割为M*N个大小规则、质心位置不同的细长垂直棱柱体。
设棱柱体的质心点坐标为(x0k,y0k),其剩余密度为ρk(k=1,2,3,…,M*N),棱柱体的底面边长分别为a和b,棱柱体高度为dh,则在频率域中各棱柱体在顶界面(x,y)点产生的重力异常△gk(x,y)的频谱为:
式中,e表示自然常数,k表示棱柱体序号。
(II)将各棱柱体在顶界面(x,y)点产生的重力异常的频谱进行叠加,得到整个深度层在顶界面(x,y)点产生的重力异常△g(x,y)的频谱为:
令
则磁源重力异常在顶界面的频谱响应为:
由式(20)得到频域中空间层剩余磁源重力异常顶界面的视密度为:
(III)对频域中空间层剩余磁源重力异常顶界面的视密度进行反傅里叶变换,得到该深度层顶面各点(x,y)的剩余密度值为:
(IV)根据泊松公式,磁化率与密度存在如下线性关系:
M=κ(x,y)*T0 (24)
式中,κ(x,y)为磁异常的视磁化率,M为磁异常的视磁化强度,ρ(x,y)为磁源重力的视密度,T0为正常磁场强度,κ'和ρ'分别是泊松公式转换时需要的已知磁化率值和密度值。
将计算得到的剩余密度值ρ(x,y),即磁源重力的视密度,代入式(23)中,计算得到磁异常的视磁化率κ(x,y)。将计算得到的磁异常的视磁化率κ(x,y)代入式(24)中,计算得到磁异常的视磁化强度M。
实施例1:为了详细描述本发明反演成像方法的过程及技术上的优势,创建一个如图3所示的柱体模型组,各模型柱的正演几何参数和物理参数如表1所示。
表1 模型的具体参数
采用本发明反演成像方法对该模型表示的地质体进行反演,其具体实施过程为:
(1)获取如图4所示的地面磁场值,并视其为地面观测数据;
(2)对地面观测磁异常数据做重力方向转换,获取如图5所示的化极磁异常;
(3)对化极偶极磁源做单极源转换,获取如图6所示的磁源重力异常;
(4)分别利用200米和500米的切割半径对地面磁源重力进行异常分离,获取如图7所示的0~200米空间层及如图8所示的200~500米空间层的地面剩余磁源重力异常;
(5)采用迭代法对0~200米空间层剩余磁源重力进行向下延拓,获取如图9所示的200米深磁性体顶界面磁源重力异常,采用迭代法对200~500米空间层剩余磁源重力进行向下延拓,获取如图10所示的500米深磁性体顶界面磁源重力异常;
(6)分别对200米深磁性体顶界面磁源重力异常和500米深磁性体顶界面磁源重力做视密度反演,获取如图11所示的200米深顶界面磁源重力视密度反演图和如图12所示的500米深顶界面磁源重力视密度反演图;
(7)分别对200米深视密度和500米深磁性体的视密度做视磁化强度转换,获取如图13所示的200米深顶界面磁源重力视强度反演图和如图14所示的500米深顶界面磁源重力视强度反演图。
从图4~图14可以看出,采用本发明磁源重力视强度成像方法,反演得到的如图13和图14所示的参数与如图3所示的模型参数相比,有很好的精确度。其中,深源反演磁化强度为75,与理论数据80比较接近,浅源磁化强度反演值为100,与理论值100吻合,这说明方法理论可行。在实际应用中,由于选取切割半径时存在一定的随意性,因此分离出来的层间剩余异常不能真实反映实际磁性体的空间位置,采用迭代法向下延拓后的“磁性体顶界面”可能与真实磁性体的顶界面并不重合,使得反演出来的强度参数与真实参数值存在一定的差异。然而,反演参数并不影响横向多个磁性体的磁性定性比较,因此,本发明将所反演的参数叫做视强度参数反演。
上述各实施例仅用于说明本发明,其中各部件的结构、连接方式和方法步骤等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。
Claims (5)
1.一种磁源重力视强度反演成像方法,其包括以下步骤:
1)将磁异常△T(x,y)转换成空间域的磁源重力异常△g(x,y);
2)选取不同的切割半径,采用插值切割法对空间域的磁源重力异常△g(x,y)进行不同深度空间层的磁源重力异常分离,获取各空间层的剩余磁源重力异常△gj(x,y,0);
3)采用迭代法将各空间层的剩余磁源重力异常△gj(x,y,0)向下延拓到各空间层的磁性体的顶界面上,获取磁性体顶界面的磁源重力异常△g(x,y,hj);
4)对磁性体顶界面的磁源重力异常△g(x,y,hj)的视强度参数进行反演,获取各空间层的视强度参数。
2.如权利要求1所述的一种磁源重力视强度反演成像方法,其特征在于:所述步骤1)中,将磁异常△T(x,y)转换成空间域的磁源重力异常△g(x,y),其包括以下步骤:
(I)将地面观测到的磁异常△T(x,y)转换成频率域的磁异常即
式中,u和v分别表示x方向和y方向的波数;
(II)将频率域的磁异常转换成化极磁异常即
式中,i表示虚分量,(P0,Q0,R0)表示磁性体所在地层的方向余弦,(P1,Q1,R1)表示磁性体在地表的方向余弦,(P2,Q2,R2)表示磁性体在磁倾角为90°的磁化方向下的方向余弦;
(III)将化极磁异常的偶极磁源作类重力的单极源处理,即对沿磁性体南北极方向积分,得到频率域的磁源重力异常即
式中,M表示磁化强度矢量,G表示万有引力常量,ρ表示磁性体的剩余密度;
(IV)对频率域的磁源重力异常进行反傅里叶变换,得到空间域的磁源重力异常△g(x,y),即
3.如权利要求1或2所述的一种磁源重力视强度反演成像方法,其特征在于:所述步骤2)中,获取各空间层的剩余磁源重力异常△gj(x,y,0),其具体包括以下步骤:
(I)选取不同的切割半径r1,r2,r3…rn,其中,n表示切割半径的个数,采用插值切割方法对磁源重力异常△g(x,y)进行不同深度空间层的磁源重力异常分离,获取不同切割半径下的浅源异常:
L1(x,y),L2(x,y),L3(x,y)…Ln(x,y),
(II)根据不同切割半径下的浅源异常,得到各空间层的剩余磁源重力异常:
式中,△gj(x,y,0)表示第j层磁性体在地面反映的磁源重力异常,即剩余磁源重力异常。
4.如权利要求1或2所述的一种磁源重力视强度反演成像方法,其特征在于:所述步骤3)中,采用迭代法向下延拓,获取磁性体顶界面的磁源重力异常△g(x,y,hj),其具体包括以下步骤:
(I)将地面ΓA上的异常值△g(x,y,0)中第s点的异常值GAs放到平面ΓB的垂直投影点s上,作为平面ΓB上的异常值即并通过傅里叶变换得到ΓB上异常值的频谱
(II)当地面ΓA与平面ΓB之间没有场源时,重力异常满足拉普拉斯方程,采用向上延拓公式:
将平面ΓB上的异常值的频谱代入上式求出地面ΓA上s点的异常值即:
式中,表示地面ΓA(z=0)上的磁源重力异常值△g(x,y,0)的傅里叶变换;F-1表示反傅里叶变换,h表示延拓高度;
(III)利用GAs与的差对平面ΓB上的异常值进行修正,得到新的异常值即:
式中,l表示步长;
(IV)重复步骤(II)和步骤(III),得到如下迭代公式:
当时,由式(16)可知ε表示一个给定的趋近于零的数,则迭代次数n取20~50次;
(V)采用与步骤(I)~(IV)相同的方法,计算得到平面ΓB上全部点的磁源重力异常值并利用向上延拓公式:
作反傅里叶变换得到△g(x,y,0)向下延拓hj深度的磁源重力异常△g(x,y,hj)。
5.如权利要求1或2所述的一种磁源重力视强度反演成像方法,其特征在于:所述步骤4)中,获取各空间层的视强度参数,其具体包括以下步骤:
(I)将引起磁源重力异常△g(x,y,hj)的非规则的磁性体分割为M*N个大小规则、质心位置不同的细长垂直棱柱体;
设棱柱体的质心点坐标为(x0k,y0k),其剩余密度为ρk,棱柱体的底面边长分别为a和b,棱柱体高度为dh,则在频率域中各棱柱体在顶界面(x,y)点产生的重力异常△gk(x,y)的频谱为:
式中,e表示自然常数,k表示棱柱体序号,k=1,2,3,…,M*N;
(II)将各棱柱体在顶界面(x,y)点产生的重力异常的频谱进行叠加,得到整个深度层在顶界面(x,y)点产生的重力异常△g(x,y)的频谱为:
令
则磁源重力异常在顶界面的频谱响应为:
则频域中空间层的剩余磁源重力异常顶界面的视密度为:
(III)对频域中空间层的剩余磁源重力异常顶界面的视密度进行反傅里叶变换,得到该深度层顶面各点(x,y)的剩余密度值为:
(IV)将剩余密度值ρ(x,y)代入式
计算得到磁异常的视磁化率κ(x,y);将计算得到的磁异常的视磁化率κ(x,y)代入式
M=κ(x,y)*T0,
计算得到磁异常的视磁化强度M;
式中,κ(x,y)为磁异常的视磁化率,M为磁异常的视磁化强度,ρ(x,y)为磁源重力的视密度,T0为正常磁场强度,κ'和ρ'分别是泊松公式转换时需要的已知磁化率值和密度值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410549240.7A CN104316972A (zh) | 2014-10-16 | 2014-10-16 | 一种磁源重力视强度反演成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410549240.7A CN104316972A (zh) | 2014-10-16 | 2014-10-16 | 一种磁源重力视强度反演成像方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104316972A true CN104316972A (zh) | 2015-01-28 |
Family
ID=52372224
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410549240.7A Pending CN104316972A (zh) | 2014-10-16 | 2014-10-16 | 一种磁源重力视强度反演成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104316972A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105334542A (zh) * | 2015-10-23 | 2016-02-17 | 中南大学 | 任意密度分布复杂地质体重力场快速、高精度正演方法 |
CN106405664A (zh) * | 2016-08-25 | 2017-02-15 | 中国科学院地质与地球物理研究所 | 一种磁异常化极方法 |
CN112114374A (zh) * | 2020-08-13 | 2020-12-22 | 天津市地球物理勘探中心 | 复杂地质体的多密度界面反演方法 |
CN112346139A (zh) * | 2020-10-15 | 2021-02-09 | 中国地质大学(武汉) | 一种重力数据多层等效源延拓与数据转换方法 |
CN112964248A (zh) * | 2021-02-09 | 2021-06-15 | 自然资源部第二海洋研究所 | 基于重力梯度张量数据的目标体实时定位方法及*** |
CN113279688A (zh) * | 2021-04-25 | 2021-08-20 | 朝阳上为软件科技有限公司 | 水平定向钻控向***及其定位方法 |
CN113486591A (zh) * | 2021-07-13 | 2021-10-08 | 吉林大学 | 一种卷积神经网络结果的重力多参量数据密度加权反演方法 |
CN114167511A (zh) * | 2021-11-26 | 2022-03-11 | 兰州大学 | 一种基于连分式展开向下延拓的位场数据快速反演方法 |
CN115220117A (zh) * | 2022-07-14 | 2022-10-21 | 招商局重庆交通科研设计院有限公司 | 异常体方位预测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6615139B1 (en) * | 2002-03-28 | 2003-09-02 | Council Of Scientific & Industrial Research | Digitally implemented method for automatic optimization of gravity fields obtained from three-dimensional density interfaces using depth dependent density |
US20030220739A1 (en) * | 2000-04-14 | 2003-11-27 | Lockheed Martin Corp. | Method of Determining Boundary Interface Changes in a Natural Resource Deposit |
US20040000910A1 (en) * | 2002-06-28 | 2004-01-01 | Tryggvason Bjarni V. | System and method for surveying underground density distributions |
CN101661115A (zh) * | 2008-08-29 | 2010-03-03 | 中国石油集团东方地球物理勘探有限责任公司 | 基于标准格架的快速三维重力、磁力物性反演的方法 |
CN101937106A (zh) * | 2010-07-27 | 2011-01-05 | 浙江大学 | 一种海底大起伏测线磁测数据的处理方法 |
-
2014
- 2014-10-16 CN CN201410549240.7A patent/CN104316972A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030220739A1 (en) * | 2000-04-14 | 2003-11-27 | Lockheed Martin Corp. | Method of Determining Boundary Interface Changes in a Natural Resource Deposit |
US6615139B1 (en) * | 2002-03-28 | 2003-09-02 | Council Of Scientific & Industrial Research | Digitally implemented method for automatic optimization of gravity fields obtained from three-dimensional density interfaces using depth dependent density |
US20040000910A1 (en) * | 2002-06-28 | 2004-01-01 | Tryggvason Bjarni V. | System and method for surveying underground density distributions |
CN101661115A (zh) * | 2008-08-29 | 2010-03-03 | 中国石油集团东方地球物理勘探有限责任公司 | 基于标准格架的快速三维重力、磁力物性反演的方法 |
CN101937106A (zh) * | 2010-07-27 | 2011-01-05 | 浙江大学 | 一种海底大起伏测线磁测数据的处理方法 |
Non-Patent Citations (3)
Title |
---|
余海龙: "位场异常三维视物性快速反演", 《中国博士学位论文全文数据库 基础科学辑》 * |
刘光海: "位场转换的谱分析方法及程序", 《中国地质科学院 矿床地质研究所所刊》 * |
吴兴等: "磁源重力技术及其应用", 《物探化探计算技术》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105334542A (zh) * | 2015-10-23 | 2016-02-17 | 中南大学 | 任意密度分布复杂地质体重力场快速、高精度正演方法 |
CN106405664A (zh) * | 2016-08-25 | 2017-02-15 | 中国科学院地质与地球物理研究所 | 一种磁异常化极方法 |
CN112114374A (zh) * | 2020-08-13 | 2020-12-22 | 天津市地球物理勘探中心 | 复杂地质体的多密度界面反演方法 |
CN112114374B (zh) * | 2020-08-13 | 2023-07-18 | 天津市地球物理勘探中心 | 复杂地质体的多密度界面反演方法 |
CN112346139A (zh) * | 2020-10-15 | 2021-02-09 | 中国地质大学(武汉) | 一种重力数据多层等效源延拓与数据转换方法 |
CN112964248A (zh) * | 2021-02-09 | 2021-06-15 | 自然资源部第二海洋研究所 | 基于重力梯度张量数据的目标体实时定位方法及*** |
CN113279688A (zh) * | 2021-04-25 | 2021-08-20 | 朝阳上为软件科技有限公司 | 水平定向钻控向***及其定位方法 |
CN113279688B (zh) * | 2021-04-25 | 2024-01-30 | 北京上为星运科技有限公司 | 水平定向钻控向***及其定位方法 |
CN113486591A (zh) * | 2021-07-13 | 2021-10-08 | 吉林大学 | 一种卷积神经网络结果的重力多参量数据密度加权反演方法 |
CN114167511A (zh) * | 2021-11-26 | 2022-03-11 | 兰州大学 | 一种基于连分式展开向下延拓的位场数据快速反演方法 |
CN114167511B (zh) * | 2021-11-26 | 2023-08-11 | 兰州大学 | 一种基于连分式展开向下延拓的位场数据快速反演方法 |
CN115220117A (zh) * | 2022-07-14 | 2022-10-21 | 招商局重庆交通科研设计院有限公司 | 异常体方位预测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104316972A (zh) | 一种磁源重力视强度反演成像方法 | |
Caratori Tontini et al. | Rapid 3‐D forward model of potential fields with application to the Palinuro Seamount magnetic anomaly (southern Tyrrhenian Sea, Italy) | |
CN105334542B (zh) | 任意密度分布复杂地质体重力场快速、高精度正演方法 | |
EP3418778A1 (en) | Systems and methods to build sedimentary attributes | |
CN105954802B (zh) | 一种岩性数据体的转换方法及装置 | |
CN106405664B (zh) | 一种磁异常化极方法 | |
CN107329171A (zh) | 深度域储层地震反演方法及装置 | |
CN102636809B (zh) | 一种传播角度域共成像点道集的生成方法 | |
Song et al. | GANSim‐3D for conditional geomodeling: Theory and field application | |
CN108446516A (zh) | 一种基于力源反演的构造断裂与地面沉降形变分解法 | |
Liu et al. | 2D inverse modeling for potential fields on rugged observation surface using constrained Delaunay triangulation | |
Yang et al. | Quaternary activity of the Beihewan Fault in the southeastern Beishan Wrench Belt, Western China: Implications for crustal stability and intraplate earthquake hazards North of Tibet | |
CN105137482A (zh) | 一种沉积体古坡度的计算方法 | |
CN106295119A (zh) | 一种页岩气地层地应力计算方法 | |
CN110007365A (zh) | 一种基于信号数据稀疏空间快速计算的联合反演方法 | |
Vitousek et al. | A nonlinear, implicit one-line model to predict long-term shoreline change | |
CN102313903A (zh) | 基于波动方程外推算子的vti介质中叠前时间偏移方法 | |
CN104598672A (zh) | 一种计算发射源姿态改变引起的电磁响应误差的方法 | |
Fu et al. | Effects of spatial distribution of fault slip on calculating co‐seismic displacement: Case studies of the Chi‐Chi earthquake (Mw7. 6) and the Kunlun earthquake (Mw7. 8) | |
Jiang et al. | Numerical simulation of the propagation of electromagnetic waves in ionospheric irregularities | |
Xue et al. | Full waveform inversion of transient electromagnetic data in the time domain | |
CN105242317B (zh) | 一种纵波速度的确定方法及装置 | |
CN107288631A (zh) | 一种井点处原始地层倾角校正方法 | |
Wang et al. | Natural Source Electromagnetic Component Exploration of Coalbed Methane Reservoirs | |
CN102455440A (zh) | Vti介质中地震波的走时计算方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150128 |