CN105334542A - 任意密度分布复杂地质体重力场快速、高精度正演方法 - Google Patents
任意密度分布复杂地质体重力场快速、高精度正演方法 Download PDFInfo
- Publication number
- CN105334542A CN105334542A CN201510698214.5A CN201510698214A CN105334542A CN 105334542 A CN105334542 A CN 105334542A CN 201510698214 A CN201510698214 A CN 201510698214A CN 105334542 A CN105334542 A CN 105334542A
- Authority
- CN
- China
- Prior art keywords
- prism
- matrix
- formula
- gravity field
- represent
- 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
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)
- Complex Calculations (AREA)
Abstract
本发明公开了任意密度分布复杂地质体重力场快速、高精度正演方法。本发明通过复杂地质体模型表示、棱柱体组合模型重力场计算(包括加权系数计算、二维离散卷积计算、重力场值合成)的步骤,实现了重力场正演计算在效率和精度上的统一。本发明解决了现有重力场正演方法不能同时保证计算效率和计算精度,无法满足大规模重力场三维密度反演、人机交互建模和解释的需求的问题。
Description
技术领域:
本发明涉及一种重力场正演方法,特别是任意密度分布复杂地质体重力场快速、高精度正演方法。
背景技术:
重力场正演是指根据密度分布计算重力场,反演是指根据观测重力值计算密度分布。正演是反演的基础,正演计算的效率直接影响反演计算的效率,而正演计算精度直接影响反演结果的质量。在重力勘探领域,伴随测绘技术和仪器的发展,重力测量在测量精度、空间分辨率和测量范围上都有了显著提高,为重力反演提供了大规模高精度、高分辨率数据,重力反演发展到三维密度反演阶段,成为国内外学者研究的热点。
随着计算机软硬件水平的不断提高,人机交互建模、解释也日益得到人们的重视,在地球物理勘探中发挥着越来越重要作用。人机交互建模能够使人们通过直观的方式对地质体进行建模,更容易融合地下结构的先验信息。反演方法与人机交互建模、解释方法相辅相成,将极大提高人们对地球内部结构的认识。交互建模过程中,首先对地下结构进行剖分,根据先验信息设计地质体分布。在勾勒出地质体分布后,进行正演计算,将正演结果与观测数据进行比对,再对模型进行调整,如此反复,实现建模。
现阶段,由于重力观测数据的增加,正演计算量急剧增大,导致在一般计算机上难以实现,成为制约三维密度反演发展的计算瓶颈。同时,人机交互建模对实时性有很高的要求,而正演计算是交互建模中最大计算量所在,正演计算的效率直接影响人机交互建模的效果。
针对正演计算,众多国内外学者进行了研究。正演计算首先对地质体剖分,然后根据剖分方式,采用某种方法计算重力场。文献(Zhdanov,M.S.,R.Ellis,S.Mukherjee.Three-dimensionalregularizedfocusinginversionofgravitygradienttensorcomponentdata.Geophysics,2004.69(4):925-937.)公开了一种结构化剖分方式,利用小棱柱体逼近复杂地质模型,在计算过程中,采用等体积的球体重力场公式近似计算小棱柱体重力场,提高了正演计算效率,但正演计算精度有所降低;文献(姚长利,郝天珧,管志宁,张聿文.重磁遗传算法三维反演中高速计算及有效存储方法技术.地球物理学报,2003.46(2):252-258.)采用结构化剖分方式,根据离散化后数学问题的特点,提出了“格架分离”技术和“格架等效计算方案”,较好解决了计算效率和计算精度问题,但对于大规模剖分情形,该文献所给出的正演方法的计算效率仍然比较低;文献(Tontini,F.C.,L.Cocchi,C.Carmisciano.Rapid3-DforwardmodelofpotentialfieldswithapplicationtothePalinuroSeamountmagneticanomaly(southernTyrrhenianSea,Italy).JournalofGeophysicalResearch,2009.114.)采用结构化剖分方法,采用三维傅里叶变换,给出了任意密度分布情形下重力异常正演的波数域表达式,借助三维快速傅里叶变换算法,实现了快速正演计算,该方法效率极高,但为克服截断效应,使用该方法前需要对剖分区域进行扩边,影响了正演计算精度;此外,还有学者采用非结构化剖分方式,采用有限元方法计算重力场,这种方法能够较精确刻画复杂地质体,计算精度高但计算效率很低。
剖分方式和计算方法共同决定了正演计算的效率和精度。正演计算的效率和精度是一对矛盾体,目前已有的正演方法存在的最大问题是不能同时保证计算效率和精度,无法满足大规模重力三维密度反演、人机交互建模和解释的需求。因此,寻找一种计算效率高、同时能保证计算精度的正演计算方法具有重要的现实意义。
发明内容:
本发明针对目前现有重力场正演方法不能同时保证计算效率和计算精度,无法满足大规模重力三维密度反演、人机交互建模和解释的需求问题,提出了任意密度分布复杂地质体重力场快速、高精度正演方法。
为解决上述技术问题,本发明采用以下技术方案:
任意密度分布复杂地质体重力场快速、高精度正演方法,其步骤为:
步骤一:复杂地质模型表示:建立包含所有目标区域的规则棱柱体模型,使得目标区域(包含起伏地形)完全嵌入在该棱柱体模型中;将该棱柱体划分成许多小棱柱体,每个小棱柱体密度为常值,不同棱柱体密度取值不同,以此刻画任意密度分布情形下复杂地质体;将位于空气部分的小棱柱体的密度值设为零,以此刻画起伏地形;
步骤二:棱柱体组合模型重力场计算:步骤一中给出的棱柱体组合模型重力场其计算公式为
式(1)中,(xm,yn,z0)表示观测点坐标,z0为常值;L表示z方向棱柱体剖分个数;M表示x方向棱柱体剖分个数;N表示y方向棱柱体剖分个数;(ξp,ηq,ζr)表示编号为(p,q,r)的小棱柱体几何中心坐标;ρ(ξp,ηq,ζr)表示该棱柱体的密度值;h(xm-ξp,yn-ηq,z0-ζr)表示加权系数;
实现上式的计算,分为三个环节:
首先,计算加权系数h(xm-ξp,yn-ηq,z0-ζr),其计算公式为
式(2)中,γ表示万有引力常数,Δx,Δy,Δz表示小棱柱体几何尺寸,arctan()表示反余切函数运算符,log()表示自然对数运算符;其它符号含义如下
x1=ξp-0.5Δx-xm,x2=ξp+0.5Δx-xm,y1=ηq-0.5Δy-yn,y2=ηq+0.5Δy-yn,z1=ζr-0.5Δz-z0,z2=ζr+0.5Δz-z0,μijk=(-1)i(-1)j(-1)k,i=1,2,j=1,2,k=1,2
其次,采用二维离散卷积快速计算方法来计算一层(相对z方向而言)棱柱体组合模型重力场,其计算公式为
式(3)中,表示第r层(r=1,2,…,L)棱柱体组合模型在高度面z0产生的重力场;(xm,yn,z0)表示离散观测点坐标;
最后,将各层棱柱体组合模型重力场进行累加,得到整个组合模型的重力场,即
步骤二中所述的二维离散卷积快速计算方法,其步骤为:
(1)将加权系数h(x1-ξp,y1-ηq,z0-ζr)排列成矩阵t,记为
式(5)中,矩阵元素ti,j与加权系数h(x1-ξp,y1-ηq,z0-ζr)存在关系
ti,j=h(x1-ξi+1,yj+1-ηq,z0-ζr)(6)
(2)将矩阵t补零扩展成矩阵
将矩阵分成四个块矩阵,记为
将块矩阵互换位置,得到矩阵cext
(3)将第r层密度值ρ(ξp,ηq,ζr)(p=1,2,…,M,q=1,2,…,N)排列成矩阵g,矩阵元素gi,j与密度值存在关系
gi,j=ρ(ξi,ηj,ζr)(10)
将矩阵g补零扩展成矩阵gext
式(11)中,0M×N表示M×N零矩阵;
(4)计算
式中,fft2()表示二维快速傅里叶变换;
(5)计算
式中,“.*”表示对应元素相乘运算;
(6)计算
式中,ifft2()表示二维快速傅里叶反变换;
(7)提取矩阵fext的前M行前N列,构成矩阵f,即为二维离散卷积计算结果。
本发明是一个有机整体,即在特定的模型表示方式条件下,建立棱柱体重力场叠加模型,根据一种特殊的加权系数计算公式,采用二维离散卷积快速计算方法,实现了重力场正演计算在效率和精度上的统一。
与现有技术相比,本发明具有以下优点:
(1)模型表示方法简单、灵活,很容易刻画任意密度分布复杂地质体以及起伏地形;
(2)能够实现任意密度分布情况下复杂地质体重力场的快速、高精度计算,可以满足大规模重力三维密度反演、人机交互建模和解释的需求;
(3)大规模正演计算时,算法不但计算效率和计算精度高,并且所需计算机内存小。
附图说明:
1、图1为重力场快速、高精度正演流程图;
2、图2为复杂地质模型表示;
3、图3为组合模型剖面图;
4、图4为组合模型重力场正演计算值;
5、图5为组合模型重力场理论值;
6、图6重力场理论值与计算值的差值;
图中符号说明如下:
L:表示z方向剖分小棱柱体个数;
M:表示x方向剖分小棱柱体个数;
N:表示y方向剖分小棱柱体个数;
ρ:表示密度;
具体实施方式:
下面结合附图对本发明中的方法作进一步详细描述。
1、复杂地质模型表示:
首先,建立包含所有目标区域的规则棱柱体模型,确定棱柱体在x,y,z方向的起始位置,使得目标区域(包含起伏地形)完全嵌入在该棱柱体模型中;
其次,根据实际问题需求,将棱柱体划分成许多规则小棱柱体(如图2所示),确定小棱柱体的几何尺寸Δx,Δy,Δz;
最后,根据目标区域的密度分布,对每个小棱柱体密度进行赋值,位于空气部分的小棱柱体,其密度值设为零;
2、棱柱体组合模型重力场计算:
步骤一中给出的棱柱体组合模型,其重力场计算公式为
式(12)中,(xm,yn,z0)表示观测点坐标,z0为常值;L表示z方向棱柱体剖分个数;M表示x方向棱柱体剖分个数;N表示y方向棱柱体剖分个数;(ξp,ηq,ζr)表示编号为(p,q,r)的小棱柱体几何中心坐标;ρ(ξp,ηq,ζr)表示该棱柱体的密度值;h(xm-ξp,yn-ηq,z0-ζr)表示加权系数;
实现式(12)的计算,分为三个环节:
首先,根据观测点坐标(xm,yn,z0)和小棱柱体几何中心坐标(ξp,ηq,ζr),计算加权系数h(xm-ξp,yn-ηq,z0-ζr),其计算公式为
式(13)中,γ表示万有引力常数,Δx,Δy,Δz表示小棱柱体几何尺寸,arctan()表示反余切函数运算符,log()表示自然对数运算符;其它符号含义如下
x1=ξp-0.5Δx-xm,x2=ξp+0.5Δx-xm,y1=ηq-0.5Δy-yn,y2=ηq+0.5Δy-yn,z1=ζr-0.5Δz-z0,z2=ζr+0.5Δz-z0,μijk=(-1)i(-1)j(-1)k,i=1,2,j=1,2,k=1,2
其次,采用二维离散卷积快速计算方法来计算一层(相对z方向而言)棱柱体组合模型重力场,其计算公式为
式(14)中,表示第r层(r=1,2,…,L)棱柱体组合模型在高度面z0产生的重力场;(xm,yn,z0)表示离散观测点坐标;
最后,将各层棱柱体组合模型重力场(r=1,2,…,L)进行累加,得到整个组合模型的重力场,即
步骤二中所述的二维离散卷积快速计算方法,其步骤为:
(1)将加权系数h(x1-ξp,y1-ηq,z0-ζr)排列成矩阵t,记为
式(16)中,矩阵元素ti,j与加权系数h(x1-ξp,y1-ηq,z0-ζr)存在关系
ti,j=h(x1-ξi+1,yj+1-ηq,z0-ζr)(17)
(2)将矩阵t补零扩展成矩阵
将矩阵分成四个块矩阵,记为
将块矩阵互换位置,得到矩阵cext
(3)将第r层密度值ρ(ξp,ηq,ζr)(p=1,2,…,M,q=1,2,…,N)排列成矩阵g,矩阵元素gi,j与密度值存在关系
gi,j=ρ(ξi,ηj,ζr)(21)
将矩阵g补零扩展成矩阵gext
式(22)中,0M×N表示M×N零矩阵;
(4)计算
式中,fft2()表示二维快速傅里叶变换;
(5)计算
式中,“.*”表示对应元素相乘运算;
(6)计算
式中,ifft2()表示二维快速傅里叶反变换;
(7)提取矩阵fext的前M行前N列,构成矩阵f,即为二维离散卷积计算结果。
下面对本发明方法的效果进行检验。
为了说明本发明所提出的方法用于计算任意密度分布情况下复杂地质构造重力场时的效率和精度,设计了如下组合模型(图3所示):
密度均匀的棱柱体内嵌一个密度均匀的球体,球心与棱柱体中心重合。棱柱体范围为:x方向从-10000m到10000m,y方向从-10000m到10000m,z方向从0m到3000m(z轴向下为正);球体半径为1000m。棱柱体密度为1g/cm3,球体的密度为5g/cm3。将棱柱体剖分成1000×1000×500个大小相同的小棱柱体,计算高度为-200m平面(图3中虚线所示)上的重力场,计算点个数为1000×1000。
正演算法利用Fortran语言编程实现,运行程序所用的个人台式机配置为:CPU为i7-2620,主频为2.7GHz,内存为32GB,四核八线程。运行所需时间约为60秒,由此可见正演算法效率很高。组合模型重力场正演算法计算值和理论值分别如图4和图5所示,从形态上看,两者是一致的。理论值减去计算值得到差值(图6所示),对差值进行统计,统计结果由表1给出,可知正演算法精度很高。
表1组合模型重力场理论值和正演计算误差统计量(单位:mGal)
最大值 | 最小值 | 均值 | 均方值 | |
理论值 | 145.54444 | 29.475854 | 91.824949 | 16.393163 |
误差 | 0.001254 | -0.000299 | 0.000107 | 0.000219 |
以上仅是本发明的优选实施方式,本发明的保护范围并不仅局限于上述实施例,凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理前提下的若干改进和润饰,应视为本发明的保护范围。
Claims (2)
1.任意密度分布复杂地质体重力场快速、高精度正演方法,其特征在于包含以下步骤:
步骤一、复杂地质模型表示:建立包含所有目标区域的规则棱柱体模型,使得目标区域(包含起伏地形)完全嵌入在该棱柱体模型中;将该棱柱体划分成许多小棱柱体,每个小棱柱体密度为常值,不同棱柱体密度取值不同,以此刻画任意密度分布情形下复杂地质体;将位于空气部分的小棱柱体的密度值设为零,以此刻画起伏地形;
步骤二、棱柱体组合模型重力场计算:步骤一中给出的棱柱体组合模型重力场其计算公式为
式(1)中,(xm,yn,z0)表示观测点坐标,z0为常值;L表示z方向棱柱体剖分个数;M表示x方向棱柱体剖分个数;N表示y方向棱柱体剖分个数;(ξp,ηq,ζr)表示编号为(p,q,r)的小棱柱体几何中心坐标;ρ(ξp,ηq,ζr)表示该棱柱体的密度值;h(xm-ξp,yn-ηq,z0-ζr)表示加权系数;
实现上式的计算,分为三个环节:
首先,计算加权系数h(xm-ξp,yn-ηq,z0-ζr),其计算公式为
式(2)中,γ表示万有引力常数,Δx,Δy,Δz表示小棱柱体几何尺寸,arctan()表示反余切函数运算符,log()表示自然对数运算符;其它符号含义如下
x1=ξp-0.5Δx-xm,x2=ξp+0.5Δx-xm,y1=ηq-0.5Δy-yn,y2=ηq+0.5Δy-yn,z1=ζr-0.5Δz-z0,z2=ζr+0.5Δz-z0,μijk=(-1)i(-1)j(-1)k,i=1,2,j=1,2,k=1,2
其次,采用二维离散卷积快速计算方法来计算一层(相对z方向而言)棱柱体组合模型重力场,其计算公式为
式(3)中,表示第r层(r=1,2,…,L)棱柱体组合模型在高度面z0产生的重力场;(xm,yn,z0)表示离散观测点坐标;
最后,将各层棱柱体组合模型重力场进行累加,得到整个组合模型的重力场,即
。
2.根据权利要求1所述的任意密度分布复杂地质体重力场快速、高精度正演方法,其特征在于:
步骤二中所述的二维离散卷积快速计算方法,其步骤为:
(1)将加权系数h(x1-ξp,y1-ηq,z0-ζr)排列成矩阵t,记为
式(5)中,矩阵元素ti,j与加权系数h(x1-ξp,y1-ηq,z0-ζr)存在关系
ti,j=h(x1-ξi+1,yj+1-ηq,z0-ζr)(6)
(2)将矩阵t补零扩展成矩阵
将矩阵分成四个块矩阵,记为
将块矩阵互换位置,得到矩阵cext
(3)将第r层密度值ρ(ξp,ηq,ζr)(p=1,2,…,M,q=1,2,…,N)排列成矩阵g,矩阵元素gi,j与密度值存在关系
gi,j=ρ(ξi,ηj,ζr)(10)
将矩阵g补零扩展成矩阵gext
式(11)中,0M×N表示M×N零矩阵;
(4)计算
式中,fft2()表示二维快速傅里叶变换;
(5)计算
式中,“.*”表示对应元素相乘运算;
(6)计算
式中,ifft2()表示二维快速傅里叶反变换;
(7)提取矩阵fext的前M行前N列,构成矩阵f,即为二维离散卷积计算结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510698214.5A CN105334542B (zh) | 2015-10-23 | 2015-10-23 | 任意密度分布复杂地质体重力场快速、高精度正演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510698214.5A CN105334542B (zh) | 2015-10-23 | 2015-10-23 | 任意密度分布复杂地质体重力场快速、高精度正演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105334542A true CN105334542A (zh) | 2016-02-17 |
CN105334542B CN105334542B (zh) | 2017-07-07 |
Family
ID=55285181
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510698214.5A Active CN105334542B (zh) | 2015-10-23 | 2015-10-23 | 任意密度分布复杂地质体重力场快速、高精度正演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105334542B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106646645A (zh) * | 2016-12-29 | 2017-05-10 | 中南大学 | 一种新的重力正演加速方法 |
CN106777598A (zh) * | 2016-12-02 | 2017-05-31 | 中南大学 | 任意磁化率分布复杂磁性体磁场梯度张量数值模拟方法 |
CN107402409A (zh) * | 2017-09-26 | 2017-11-28 | 西南石油大学 | 一种三维不规则地层起伏界面重力正演方法 |
CN107942399A (zh) * | 2017-11-23 | 2018-04-20 | 桂林理工大学 | 一种大距离位场向上延拓计算方法 |
CN108490496A (zh) * | 2018-03-26 | 2018-09-04 | 中国石油化工股份有限公司 | 基于拟径向基函数神经网络的重力场密度反演方法 |
CN108984939A (zh) * | 2018-07-30 | 2018-12-11 | 中南大学 | 基于3D Gauss-FFT的三维重力场正演方法 |
CN109283589A (zh) * | 2018-08-20 | 2019-01-29 | 桂林理工大学 | 一种重力场水平分量的获取方法 |
CN109313663A (zh) * | 2018-01-15 | 2019-02-05 | 深圳鲲云信息科技有限公司 | 人工智能计算辅助处理装置、方法、存储介质、及终端 |
CN109375280A (zh) * | 2018-12-10 | 2019-02-22 | 中南大学 | 一种球坐标系下重力场快速高精度正演方法 |
CN109490978A (zh) * | 2019-01-08 | 2019-03-19 | 中南大学 | 一种起伏地层的频率域快速高精度正演方法 |
CN112287534A (zh) * | 2020-10-21 | 2021-01-29 | 中南大学 | 基于nufft的二维磁异常快速正演模拟方法和装置 |
CN112800657A (zh) * | 2021-04-15 | 2021-05-14 | 中南大学 | 基于复杂地形的重力场数值模拟方法、装置和计算机设备 |
CN116911146A (zh) * | 2023-09-14 | 2023-10-20 | 中南大学 | 三维重力场全息数值模拟及cpu-gpu加速方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040201585A1 (en) * | 2003-03-31 | 2004-10-14 | Council Of Scientific And Industrial Research | Generation of three dimensional fractal subsurface structure by Voronoi Tessellation and computation of gravity response of such fractal structure |
CN101975969A (zh) * | 2010-10-19 | 2011-02-16 | 华中科技大学 | 一种基于全张量重力梯度反演的水下目标探测方法 |
CN104316972A (zh) * | 2014-10-16 | 2015-01-28 | 中国海洋石油总公司 | 一种磁源重力视强度反演成像方法 |
-
2015
- 2015-10-23 CN CN201510698214.5A patent/CN105334542B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040201585A1 (en) * | 2003-03-31 | 2004-10-14 | Council Of Scientific And Industrial Research | Generation of three dimensional fractal subsurface structure by Voronoi Tessellation and computation of gravity response of such fractal structure |
CN101975969A (zh) * | 2010-10-19 | 2011-02-16 | 华中科技大学 | 一种基于全张量重力梯度反演的水下目标探测方法 |
CN104316972A (zh) * | 2014-10-16 | 2015-01-28 | 中国海洋石油总公司 | 一种磁源重力视强度反演成像方法 |
Non-Patent Citations (2)
Title |
---|
MICHAEL S.ZHDANOV 等: "Three-dimensional regularized focusing inversion of gravity gradient tensor component data", 《GEOPHYSICS》 * |
刘银萍: "基于ExtrapolationTikhonov正则化算法的重力数据及梯度多分量数据的3D反演方法研究", 《中国博士学位论文全文数据库 基础科学辑》 * |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106777598A (zh) * | 2016-12-02 | 2017-05-31 | 中南大学 | 任意磁化率分布复杂磁性体磁场梯度张量数值模拟方法 |
CN106777598B (zh) * | 2016-12-02 | 2020-01-14 | 中南大学 | 任意磁化率分布复杂磁性体磁场梯度张量数值模拟方法 |
CN106646645A (zh) * | 2016-12-29 | 2017-05-10 | 中南大学 | 一种新的重力正演加速方法 |
CN106646645B (zh) * | 2016-12-29 | 2018-01-19 | 中南大学 | 一种重力正演加速方法 |
CN107402409A (zh) * | 2017-09-26 | 2017-11-28 | 西南石油大学 | 一种三维不规则地层起伏界面重力正演方法 |
CN107942399A (zh) * | 2017-11-23 | 2018-04-20 | 桂林理工大学 | 一种大距离位场向上延拓计算方法 |
CN109313663B (zh) * | 2018-01-15 | 2023-03-31 | 深圳鲲云信息科技有限公司 | 人工智能计算辅助处理装置、方法、存储介质、及终端 |
CN109313663A (zh) * | 2018-01-15 | 2019-02-05 | 深圳鲲云信息科技有限公司 | 人工智能计算辅助处理装置、方法、存储介质、及终端 |
CN108490496A (zh) * | 2018-03-26 | 2018-09-04 | 中国石油化工股份有限公司 | 基于拟径向基函数神经网络的重力场密度反演方法 |
CN108490496B (zh) * | 2018-03-26 | 2019-11-08 | 中国石油化工股份有限公司 | 基于拟径向基函数神经网络的重力场密度反演方法 |
CN108984939A (zh) * | 2018-07-30 | 2018-12-11 | 中南大学 | 基于3D Gauss-FFT的三维重力场正演方法 |
CN109283589A (zh) * | 2018-08-20 | 2019-01-29 | 桂林理工大学 | 一种重力场水平分量的获取方法 |
CN109375280A (zh) * | 2018-12-10 | 2019-02-22 | 中南大学 | 一种球坐标系下重力场快速高精度正演方法 |
CN109490978A (zh) * | 2019-01-08 | 2019-03-19 | 中南大学 | 一种起伏地层的频率域快速高精度正演方法 |
CN109490978B (zh) * | 2019-01-08 | 2020-04-28 | 中南大学 | 一种起伏地层的频率域快速高精度正演方法 |
CN112287534A (zh) * | 2020-10-21 | 2021-01-29 | 中南大学 | 基于nufft的二维磁异常快速正演模拟方法和装置 |
CN112287534B (zh) * | 2020-10-21 | 2022-05-13 | 中南大学 | 基于nufft的二维磁异常快速正演模拟方法和装置 |
CN112800657A (zh) * | 2021-04-15 | 2021-05-14 | 中南大学 | 基于复杂地形的重力场数值模拟方法、装置和计算机设备 |
CN116911146A (zh) * | 2023-09-14 | 2023-10-20 | 中南大学 | 三维重力场全息数值模拟及cpu-gpu加速方法 |
CN116911146B (zh) * | 2023-09-14 | 2024-01-19 | 中南大学 | 三维重力场全息数值模拟及cpu-gpu加速方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105334542B (zh) | 2017-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105334542A (zh) | 任意密度分布复杂地质体重力场快速、高精度正演方法 | |
CN102798898B (zh) | 大地电磁场非线性共轭梯度三维反演方法 | |
CN104866653B (zh) | 一种获取地下三维密度结构的方法 | |
CN106777598A (zh) | 任意磁化率分布复杂磁性体磁场梯度张量数值模拟方法 | |
CN106855904B (zh) | 一种二度体重力异常计算方法 | |
CN105549106A (zh) | 一种重力多界面反演方法 | |
CN110045432A (zh) | 基于3d-glq的球坐标系下重力场正演方法及三维反演方法 | |
CN104375195A (zh) | 时频电磁的多源多分量三维联合反演方法 | |
CN106646645A (zh) | 一种新的重力正演加速方法 | |
CN107024723B (zh) | 一种二度体磁场数值计算方法 | |
CN102636809B (zh) | 一种传播角度域共成像点道集的生成方法 | |
CN105954802A (zh) | 一种岩性数据体的转换方法及装置 | |
CN105388520A (zh) | 一种地震资料叠前逆时偏移成像方法 | |
CN112147709A (zh) | 一种基于部分光滑约束的重力梯度数据三维反演方法 | |
CN106483559A (zh) | 一种地下速度模型的构建方法 | |
CN109254327B (zh) | 三维强磁性体的勘探方法及勘探系统 | |
CN103454677A (zh) | 基于粒子群与线性加法器结合的地震数据反演方法 | |
CN114943178A (zh) | 一种三维地质模型建模方法、装置及计算机设备 | |
CN102034271A (zh) | 基于起伏地形的三维模型单元重磁异常快速处理方法 | |
CN107748834B (zh) | 一种计算起伏观测面磁场的快速、高精度数值模拟方法 | |
CN107942399A (zh) | 一种大距离位场向上延拓计算方法 | |
Tang et al. | Topographic effects on long offset transient electromagnetic response | |
CN109490978A (zh) | 一种起伏地层的频率域快速高精度正演方法 | |
CN103745118B (zh) | 一种基于磁偶极子等效源法的地磁异常数据网格化方法 | |
CN112596113A (zh) | 一种基于重力不同阶梯度特征值交点的场源位置识别方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |