首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
多层等效源曲面磁异常转换方法   总被引:1,自引:0,他引:1       下载免费PDF全文
李端  陈超  杜劲松  梁青 《地球物理学报》2018,61(7):3055-3073
在磁异常数据处理中,利用等效源技术重构磁异常场具有较好的稳定性和较高的计算精度,因而被广泛应用.传统方法是采用设置在近地表的单层等效源拟合实测磁异常数据,尽管拟合精度很高,但向上延拓之后往往会出现较大的拟合误差,即可能存在磁异常信号的"泄漏",尤其在原始数据中存在背景场时更容易出现此种误差.本文提出一种多层等效源技术方案,应用分布于不同深度范围内的等效源模拟实测数据,减少了等效源参数设置的盲目性.理论模型试验表明,采用多层等效源方法重构的磁异常及其梯度与分量,较单层等效源方法具有更高精度,可以吸收更完整的实测磁异常信息.论文详细地讨论了如何优化多层等效源设置、等效源参数选择以及计算方法,通过二维和三维理论模型试验,验证了在复杂条件下多层等效源方法的可行性和适应性,并且将该方法应用于广西某地的实测磁异常数据转换之中,取得了较好的应用效果.  相似文献   

2.
布格重力异常数据是研究地壳结构和变形的重要依据,由于观测手段、处理方法等存在差异,不同手段获得的布格重力异常可能存在数据噪声水平及起算基准不一致,从而引起多源布格重力异常数据难以有效融合的难题.为此,本文在传统等效源重力数据融合方法的基础上,通过引入贝叶斯参数优化准则,提出了一种新的多源布格重力异常融合算法.随后对川滇地区两条实测重力剖面布格异常和WGM2012模型计算布格异常开展了融合实验,结果表明:新算法不仅能显著降低不同异常中的非相干噪声,而且能将多源异常的系统偏差有效纠正.与已有布格重力异常相比,融合所得布格重力异常特征与川滇地区主要线性构造分布形态具有更好的一致性.本文给出的川滇地区融合布格重力异常同时兼顾了地表实测重力的高精度和卫星重力场全球模型的高分辨率,可为研究川滇地区地壳结构和构造区划提供更加有价值的地球物理基础资料.  相似文献   

3.
为了提高重力异常导数换算的精度,真实有效地反映地质体的异常特征,提出用余弦变换计算异常导数的新方法. 给出并证明了两个定理,利用它们推导出重力位余弦谱一般表达式以及重力异常各阶导数计算公式,建立了位场余弦谱分析理论. 模型实验中发现,用Fourier变换计算的一阶导数与理论导数偏差很大,而余弦变换计算的导数与理论异常导数拟合效果非常好,除边界几个数据因重力异常的有限截断产生的吉布斯效应残留使误差较大外,数据的计算精度均很高,误差为-009%~5%.  相似文献   

4.
利用宜川—泰安剖面高精度重力观测及GPS测量结果,计算得到该剖面较为精细的垂向一阶导数梯度和布格重力异常归一化总梯度图像.结合前人相关研究结果,本文认为垂向一阶导数梯度主要反映浅部异常场源的位置和低密度地质体的特性.利用该梯度图像,对宜川—泰安剖面经过的7条主要断裂及其空间分布特征进行分析解释,并推测各断裂的可能埋深和性质.   相似文献   

5.
本文导出了利用重力与磁异常及其高阶导数的功率谱计算密度体和磁性体上顶与下底埋深的近似公式,还提供了利用垂向或水平导数谱计算深度的方法。理论模型的试验结果表明,由近似公式计算的精度可以满足实际需要。本文最后根据一条重力剖面资料,计算了密度层上、下界面。  相似文献   

6.
区域重力异常蕴含丰富的断裂构造特征信息.采用小波多尺度分解、归一化总水平导数垂向导数(NVDR_THDR)和剖面2.5D重力异常反演方法对北京平原区高精度区域重力资料进行了处理,获得了主要断裂构造平面位置、长度、规模、汇交关系及深浅延伸特征,分析了利用重力异常识别断裂的效果.结果表明:(1)北京平原区NE向断裂平面延展长度大、连续性强,NW向断裂分段性明显,连续性弱,但NW向断裂对NE向断裂有切割改造迹象.近EW或SN向断裂分布较为局限.(2)主控断裂平面和垂向延伸特征差异明显,顺义断裂、孙河断裂、永定河断裂等为盖层断裂;南口断裂、孙河断裂西段、二十里长山断裂和张喜庄断裂为基底断裂;黄庄—高丽营断裂南段、南苑—通县断裂、礼贤断裂、夏垫断裂和皮各庄断裂东段为地壳断裂.(3)利用NVDR_THDR峰值异常带的连续性、幅值、宽度及错切关系可有效识别断裂的平面展布特征,通过小波变换获得的不同深度等效层异常结合2.5D剖面反演可有效研究断裂在基岩内部的深浅延伸情况,但重力异常的垂向分辨能力弱,断裂在新生界内部延伸特征需结合其他地质资料进一步分析.  相似文献   

7.
位场总水平导数极值位置空间变化规律研究   总被引:5,自引:2,他引:3       下载免费PDF全文
王万银 《地球物理学报》2010,53(9):2257-2270
通过对位场总水平导数函数性质的研究表明,位场总水平导数不是位函数,因而利用位场总水平导数构造新的边缘识别方法时会出现"奇点",使得计算结果的稳定性下降.对单一边界、双边界、多边界模型重力异常总水平导数和重力异常垂向导数总水平导数极值位置的空间变化规律研究表明,重力异常垂向导数总水平导数和化极磁力异常总水平导数的极值位置相同,与重力异常总水平导数的极值位置空间变化规律相似.利用位场总水平导数极大值位置能够准确识别单一直立边界地质体的边缘位置,但不能准确识别其它任何形体的边缘位置,其识别结果的偏移量大小随地质体的埋深、水平尺寸以及倾斜程度等变化,但能收敛于某一固定值;重力异常垂向导数总水平导数比重力异常总水平导数的峰值更加尖锐、横向识别能力更强,其极大值位置更靠近地质体上顶面边缘位置,但存在"次极大值"的影响.  相似文献   

8.
被评论的论文(简称原文)中应用离散余弦变换(DCT法)计算磁异常导数,并与离散Fourier变换法(DFT法)计算结果进行对比,得到“DCT法明显优于DFT法”的结论.我们对原文进行了验证.原文磁异常导数波数响应计算的错误导致DFT法计算的错误,最终导致结论的错误.原文相关两篇论文中应用DCT法分别进行重力异常的导数换算和重力异常的向上延拓,也存在同样的错误.  相似文献   

9.
不同数据用于欧拉方程的模型计算   总被引:3,自引:3,他引:0       下载免费PDF全文
在欧拉反褶积方法的研究中,过去常直接用原始位场异常值及其导数进行反演,由于数据的局限性,有时导致反演结果的发散.为改善反演结果的收敛性,本文中采用了垂直导数、水平导数及解析信号作为源数据进行欧拉反演.(1)推导了重、磁异常不同数据类型的欧拉方程;(2)找出了重、磁异常不同数据类型的欧拉方程之间的联系和差异;(3)通过模型计算比较,证明用导数或解析信号进行欧拉反演计算,可改善反演效果,提高反演结果的收敛性.  相似文献   

10.
利用卫星重力数据计算地球内部密度异常   总被引:5,自引:0,他引:5  
地球外部重力场是由地球内部物质分布决定的。本文利用RAPP81180阶卫星重力位系数根据异常源埋深与重力位系数阶数的关系和重力场分离方法计算了经度为90°E的地球密度异常剖面,得到了地球内部不同深度的密度异常结果,并初步作了地球物理解释.计算结果得到的深度为65km处的低密度异常带与地球模型的LID层和C层较为接近。  相似文献   

11.
2014年8月跨越六盘山东麓断裂带,在长度约为200公里的剖面上展开了流动重力/GPS联合剖面观测.观测结果表明,测线东端的布格重力异常约为-190 m Gal(10~(-5)ms~(-2)),西端则为-250 m Gal左右.在假设地壳均一的前提下,基于Airy均衡模型,利用布格重力异常和GPS观测数据,分别计算了测线所在剖面的莫霍面深度与均衡面深度,发现六盘山地区处于正均衡异常状态.使用布格重力异常数据反演六盘山地区的地壳密度结构,并据此地壳分层结构,计算了六盘山地区均衡面与莫霍面深度,对比显示该区域亦处于均衡正异常状态.为了确定青藏高原东北缘岩石圈有效弹性厚度和六盘山隆升机制,我们利用EMG2008自由空气异常和SRTM V18.1 DEM数据,使用自由空气重力异常导纳方法,研究了以六盘山地区为中心的青藏高原东北缘岩石圈有效弹性厚度(Te)和加载机制,发现六盘山地区的Te为5 km,岩石圈加载主要来自于地表,占总加载的95%.最后,对比六盘山地区Airy均衡异常与弹性板均衡异常,发现六盘山东麓断层处地壳承载梯度值最大,表明该断裂带吸收了较多的应变能.  相似文献   

12.
低磁纬度地区受斜磁化的影响,用常规方法很难准确确定磁性体的平面分布特征.为了研究位于低磁纬度地区普图马约盆地的磁性体分布特征,本文根据场的散度原理,假定磁△T异常为具有一定方向的矢量场,其方向与磁化方向一致,导出了磁△T异常视散度的计算方法;根据磁位与引力位的关系,在频率域中通过磁△T异常求取了磁源重力异常,并尝试利用拉普拉斯方程计算磁源重力异常垂向二阶导数.本文设计理论模型讨论了磁△T异常视散度、磁源重力异常垂向二阶导数的特征与磁性体平面分布特征的关系,证明了上述方法的有效性.进而利用上述方法推测了普图马约盆地磁性体的平面分布特征,结果表明:应用磁视散度及磁源重力异常确定的普图马约盆地磁性体分布与实际地质特征吻合较好,取得的成果对普图马约盆地相关地质研究及对低磁纬度地区的磁性体的确定有一定的参考意义.  相似文献   

13.
随着在石油勘查、矿产勘查、构造研究等地质工作中提出的应用位场异常研究规模相对较小的地质体及其细节的要求不断提高,传统重力方法在精度和分辨率方面已不能满足要求.重力梯度张量作为重力位二次导数,较之重力异常,能够更直观地反映地下密度体的赋存状态,基于重力梯度张量数据的反演能够提高地质体的定量模拟质量.但是由于受国内仪器的发展水平限制,直接获取实测的重力张量数据难度较大,目前对于重力张量数据的获取主要是通过数值计算的方法,把重力异常测量值变换成张量信息.本文是在样条理论的基础上,提出新的重力梯度张量计算方法,利用三次样条插值方法拟合重力异常,进一步地通过积分导数关系求取样条函数系数,实现了由重力异常求取重力梯度张量的数据计算.在此基础上,利用张量数据进行欧拉反演计算,对比常规欧拉反演结果,进一步分析张量欧拉反演精度.球体及水平板状体模型试验结果显示,张量欧拉反演更能有效地完成目标体的反演,其结果更加准确.最后,将张量欧拉反演应用到梅山铁矿区重力实测数据解释中,准确地反演出梅山主矿体的水平投影范围及矿体中心埋深.同时,相比常规欧拉反演,张量欧拉反演在刘家村附近、梅山村南部及黄林库—吴家洼地段显示出有更加收敛的欧拉解,这为梅山地区进一步找矿工作提供了重要的理论支撑.  相似文献   

14.
本文提出一种基于重力/GPS联合观测数据计算垂向构造应力的新方法.计算步骤如下:(1)通过重力/GPS联合观测数据计算布格重力异常;(2)依据布格重力异常数据推算莫霍面深度;(3)依据GPS观测数据,通过均衡理论计算均衡面深度;(4)依据莫霍面与均衡面之间剩余物质(壳幔物质密度差)所承受的附加浮力,计算地壳承载的垂向构造应力.本文利用上述构造应力新算法,计算了巴颜喀拉块体东边界及周边地区垂向构造应力分布,发现龙泉山断裂带以东地区垂向构造应力基本为零,龙泉山断裂带与龙门山断裂带之间地区垂向构造应力为正值,巴颜喀拉地块东部垂向构造应力为负值.鲜水河断裂带东南段周边蓄积了-40~-50 MPa的垂向构造应力,且梯度变化剧烈;松潘高原蓄积的垂向构造应力大约为-10~-20 MPa,相对较小.  相似文献   

15.
位场垂向导数零值位置空间变化规律研究(英文)   总被引:7,自引:4,他引:3  
本文推导了简单规则形体(单一边界、双边界和多边界)重力异常垂向一阶导数和重力异常垂向二阶导数零值位置的解析表达式,并研究了其空间变化规律;对难以得到其零值位置解析表达式的简单规则形体,利用其剖面图和断面图上的零值位置来研究其空间变化规律。研究结果表明,重力异常垂向二阶导数与重力异常垂向一阶导数零值位置的空间变化规律基本一致,只是重力异常垂向二阶导数零值位置比重力异常垂向一阶导数的零值位置更靠近地质体上顶面边缘位置,且其分辨能力更强。对于单边界模型,随着埋深的增加,重力异常垂向导数的零值位置均自形体上顶边缘位置向外侧偏移,但最终均能收敛于某一固定值;对于双边界模型,随着埋深的增加,重力异常垂向导数的零值位置自形体上顶边缘位置一直向外侧偏移;对于多边界模型,随着埋深的增加,重力异常垂向导数零值位置自形体上顶边缘位置向外侧收敛,部分边界的零值位置重合并消失。最后通过实际资料检验了方法的有效性和可靠性。  相似文献   

16.
我国在海域开展了大规模的航空重力勘探,这些资料对构建高精度大地水准面具有重要价值.基于此,本文提出一种利用海域航空重力测量数据快速构建大地水准面的方法.该方法基于移去-恢复法思想,利用位场最小曲率方法对航空重力数据进行高精度向下延拓并获取相应的扰动位,实现航空重力测量快速构建海域大地水准面.与斯托克斯积分计算相比,采用了处理效率更高的频率域位场转换,解决了向下延拓及垂向积分时航空重力异常数据空白及扩边问题,具有较高的位场转换精度.本文应用EGM2008模拟航空重力数据进行模型验证,计算结果与其给出的水准面的精度相当;同时,也选取GRAV-D计划的航空重力数据进行实际验证,计算结果与xGEOID18B水准面模型精度基本一致.模型验证和实际应用验证了本方法的实用性.  相似文献   

17.
李红蕾  韩建成  陈石  侍文  张贝 《地震》2021,41(1):129-140
在综合分析前人研究成果的基础上, 采用基于贝叶斯准则的等效源模型滤波增强方法来识别地质体边界。 通过理论模型实验讨论了在相同的噪声水平情况下, 不同算法对于计算结果的影响, 模型实验结果表明, 与直接计算和传统高斯滤波增强技术相比, 新算法可以有效消除传统算法中普遍存在的对于噪声的过度敏感问题, 所得到的地质体边界形态失真更小。 随后, 利用此方法对川滇地区的重力数据进行了处理, 并对该地区的重力水平总梯度(THD)、 重力垂直梯度模(AVG)、 重力解析信号模量(AAS)、 重力水平梯度垂向导数(VGHD)应用效果进行了对比分析。 研究结果提供了川滇地区的构造分区信息, 对该地区构造演化、 边界定位及动力学过程分析具有重要的指导意义。  相似文献   

18.
电性可控源导线长度是可控源勘探中的重要参数.本文借助有限差分法对三维异常体进行正演模拟,讨论了可控源导线长度对三维可控源异常的影响.计算了深部低阻异常体和高阻异常体的异常变化情况.从数值结果可以看出:(1)Ex分量和Hy分量场值大小与源长度的变化近似成正比;(2)当源长度增加到收发距的二分之一时,源中心点处的Ex分量和Hy分量场值增加倍数明显小于远离中心点处的场值增加倍数,进而对异常形态产生影响;(3)源长度的变化对于阻抗视电阻率没有影响.从而总结出电性可控源长度对于深部三维异常体响应的影响规律.  相似文献   

19.
基于布格重力异常相对于地形起伏光滑分布的约束条件,从一维自由空气重力异常数据出发,采用贝叶斯方法估算近地表岩石密度,同时采用三次B样条函数拟合布格重力异常,获取光滑分布的布格重力异常.数据拟合和光滑约束之间的权重采用Akaike贝叶斯准则(ABIC准则)自动确定.均匀剖分模型和不均匀剖分模型数据试验都验证了该方法的有效性.相关参数评价表明,足够多的样条系数可以提高估计结果的准确性,样条系数的个数接近测点数时可获得较稳定的估计结果.增大异常的噪声水平时,ABIC准则可有效地自动增大先验光滑约束的权重.云南地区两条重力剖面应用结果表明,剖面沿线的近地表密度值起伏变化明显(达2.45~2.8g·cm^-3),前寒武纪和古生代地层密度相对较高(主要为2.53~2.75g·cm^-3),而中生代密度较低(2.45~2.73g·cm^-3);本文估计的近地表密度结果与区域物性资料及地表地质特征较吻合;估计的剖面布格重力异常具有光滑性;红河断裂两侧近地表密度差异较大,可达0.4g·cm^-3.本文获得的两条剖面近地表密度结构和布格重力异常为该区深部结构与构造研究提供更可靠的重力基础数据.  相似文献   

20.
实用化的重力数据三维反演过程中需要大量的计算时间和存储空间,针对此问题本文提出一种组合优化算法:(1)利用空间域的精确性和波数域的快速性,提出一种高精度快速三维正演算法;(2)利用反演矩阵中的对称性,将重力共轭梯度反演(CG)中的主要计算分解成两次正演计算,从而优化了三维重力反演的计算效率。通过不同网格剖分数量模型的数值实验,验证了本文优化算法的计算精度和计算效率。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号