共查询到16条相似文献，搜索用时 171 毫秒
1.
针对重力与地震联合反演存在的问题,结合已有的研究成果,本文研究实现了速度和密度随机分布共网格单元模型的建模技术,以适应密度和速度剧烈变化的复杂模型及联合反演的计算要求.重力正演利用了该网格的二度半体模型,并进一步改进了地震走时的二维射线追踪计算方法,以适用于速度随机分布的网格介质.结合改进的模拟退火算法,实现了这种共网格条件下的重力与地震资料的同步联合反演.模型试验证明了重力与地震联合反演可以准确确定复杂物性界面的密度和速度结构,适用于物性界面不完全一致和物性变化剧烈的复杂模型,并且联合反演结果要优于单独的重力反演.带先验信息约束下的实际资料的联合反演,进一步证明了该方法的适用性和效果,可提高反演精度并减少多解性. 相似文献
2.
重力反演是恢复地下密度空间分布的有效工具,而选择合理的密度模型约束方法是提升重力反演分辨率和可靠性的关键.常规约束方法大多是从剖分网格空间中的密度模型出发,通过调整光滑或稀疏约束权重来匹配反演目标,但当地质体类型多样、异常分离不准确及网格剖分方案不合理时,模型约束的合理性与灵活性难以得到有效保证.为此,本文提出了一种基于密度模型稀疏表征的重力反演方法.首先假设待反演的密度模型表征为模型特征矩阵和稀疏分解系数的线性组合,之后重新推导了重力反演目标函数,并给出了分解系数的稀疏求解过程.相比现有重力反演方法,用于构建模型特征矩阵的特征模型可包含不同类型地质体的先验几何信息,分解系数的稀疏性保证了待反演目标来自于最典型的地质模式组合.最后,通过模型试验及实际资料验证了基于密度模型稀疏表征的重力反演方法的有效性. 相似文献
3.
利用重力异常揭示地壳三维密度结构是重力资料地质解释的目标和任务,密度反演的好坏至关重要.本文对华北克拉通进行了重力异常多尺度密度反演研究,首先利用小波变换对重力异常进行多尺度分解,接着利用功率谱分析方法估算各层场源的平均深度,然后利用广义密度反演方法进行各层密度反演,对华北地区进行了多层密度反演,得到其密度结构,并进行地质解释和油气藏分析.结果显示了重力场多尺度密度反演方法的有效性,对华北地区密度结构的地质含义进行了初步分析,位于造山带中的低密度异常主要反映沿造山带展布的花岗岩体.对比华北东部区浅层密度扰动与油气坳陷位置,发现油气坳陷都表现为低幅度的低密度异常区.说明利用小波多尺度反演提取的密度信息对油气勘探的部署有一定指导意义. 相似文献
4.
常规协克里金方法反演重力或重力梯度数据具有抗噪性好、加入先验信息容易等优点,其反演的地下密度分布能够识别异常体中心位置,还原异常体基本形态,但反演图像光滑,分辨率低,这是由于常规方法估计的密度协方差矩阵全局发散、平稳.为了通过协克里金方法获得聚焦的密度分布需要改善密度协方差矩阵的性质.首先,本文推导了理论密度协方差公式,其性质表明,当理论模型聚焦分布时,其密度协方差矩阵是非平稳且聚焦分布的.为了打破常规协方差矩阵全局平稳、发散的特征,本文设置密度阈值处理协方差矩阵,通过不断更新协方差矩阵来迭代实现协克里金反演,最终得到相对聚焦的反演结果.用本文方法处理重力与重力梯度数据恢复两种密度模型,均得到了与正演模型匹配的反演结果;再将方法运用于文顿盐丘的实际测量重力与重力梯度数据,反演结果与已知的地质情况匹配较好. 相似文献
5.
《地震地质》2017,(4)
深层地震资料由于受到火山岩体或风化壳盖层的影响,质量品质往往不高,难以识别盖层之下的小规模的潜山构造。而重力勘探技术恰能利用潜山构造与上覆地层存在显著密度差异的特点来识别这种特殊的地质储层。文中以华北油田虎8北潜山构造区为例,在该地区地层密度模型的基础上进行了三维重力正演。正演结果表明,深部密度异常体在地表形成的低缓重力异常,与区域或局部背景重力异常场叠加在一起,在布格重力异常曲线上无法识别目标体异常。通过利用高阶导数、滑动滤波等处理方法把异常场从区域场中剥离出来,可发现虎8北潜山目标区存在小规模局部剩余重力高。在地震资料的约束下,根据剩余重力异常的幅值及地层密度差进行密度建模,实现了三维密度反演。反演的结果确定了虎8北潜山的存在,并给出了潜山可能的埋深与规模。证明了三维高精度重力方法在确定深部潜山构造上具有较明显的地质效果,这对于探测深部低速层覆盖之下的潜山构造油气藏具有重要的指导意义。 相似文献
6.
三维重力反演是地质工作者了解地球深部构造,认知地下结构的重要手段.按照反演单元划分,三维重力反演有离散多面体(Discrete)反演和网格节点(Voxels)反演两种方式.离散多面体反演由于易于吸收先验地质信息得到的理论场能够很好地拟合观测场,因此,在实际重力反演中更受欢迎.目前离散多面体重力反演中初始模型的建立方法繁杂不一,实际应用受到很大的限制.本文本着充分挖掘利用先验信息和重力观测数据得到丰富可靠的反演结果这一原则,以离散多面体反演技术为基础,改进建模过程.在初始模型的建立中,吸收贝叶斯算法优势,采用隐马尔科夫链改善朴素贝叶斯方法的分类效果,通过最大似然函数算法求解,再采取模型降阶技术,固定所建模型中几何体的形态或密度,达到在几何体形态(x,y,z)、密度(σ)和重力值(g)五个参数中降低维数目的,从而减小高维不确定性和正演的计算量,由此反演计算的地质体密度和分布范围相对更准确,更利于重现重力模型结构.通过单位球体和任意形态几何体模拟实验,以及安徽省泥河矿区三维重力反演实践,得到非常接近实际的密度或重力值,大幅提高了三维重力反演的精度和效率,说明该方法是有效、实用的. 相似文献
7.
深部大尺度单一密度界面重力异常迭代反演 总被引：4，自引：3，他引：4
对于深部大尺度单一密度界面的重力异常反演,需顾及球形观测面和重力铅垂线变化的影响。为此,本文通过球谐展开,得到重力异常的级数展开式,并推导出积分形式的反演迭代解。该迭代方法要求设法分离出单一密度界面的重力异常,已知界面上、下层的密度差和尽可能多的深度控制点。模型和实例显示其效果良好,可用于莫霍界面、岩石层底界面等深部密度界面重力异常的反演。 相似文献
8.
V字型密度界面是一类常见的密度界面,如海沟、半地堑以及俯冲带之下的莫霍面,利用重力数据刻画此类密度界面形态对于区域构造研究、油气勘探以及物理海洋学等都具有重要意义.本文首先建立了Lp-范数形式的模型约束函数,并利用正则化原理将其与重力数据误差函数和已知深度约束函数结合形成V字型密度界面反演的目标函数,推导了目标函数的梯度表达式,并以非线性共轭梯度法为核心给出了反演流程.二维简单模型试算结果表明p=5时该方法能准确地刻画V字型密度界面起伏特征,且亦能准确地应用于二维复杂密度界面和三维界面的反演.最后将反演方法应用于挑战者深渊及邻区的实际资料处理之中,利用研究区海底地形数据和沉积层厚度数据对自由空间重力异常逐层剥离而得到莫霍面引起的重力异常,用本文方法对此重力异常进行反演,结果呈现了板块俯冲作用引起的V字型莫霍面起伏特征. 相似文献
9.
地质体的重力建模是正确解释和应用重力资料的关键问题之一.针对非规则形状变密度的三度体,本文提出了基于3D Delaunay剖分算法的重力建模方法.采用3D Delaunay剖分算法将三维目标地质体分解为若干变密度四面体体元,推导了基于四面体体元的重力正演公式,建立了剩余密度值与重力异常值的线性方程组;以变密度的长方体和倾斜台阶组合体为例,比较分析了常规块体算法和3D Delaunay剖分算法应用于重力正演的有效性,并采用共轭梯度法加密度约束条件对非规则形状变密度的倾斜台阶组合体进行了密度反演.计算结果验证了本文方法的正确性和有效性.基于3D Delaunay剖分算法的重力建模可应用于存在褶皱、断层、裂缝等复杂地质体的重力正反演计算. 相似文献
10.
直接迭代法是利用重力异常反演密度界面的一种重要方法,对迭代公式分析的结果表明,可在反演公式中添加加速收敛因子以减少反演的迭代次数从而提高效率.对于平缓的密度界面采用较小的加速收敛因子、而对起伏较大的界面采用较大的加速收敛因子可有效减少反演迭代次数.为此,本文利用重力异常总水平导数判断待反演界面起伏平缓和起伏较大的区域,并以此为基础,将总水平导数计算结果归一化至1.4～2.2之间而形成全区变化的加速收敛因子.假设待反演密度界面上下的密度差沿横向变化,将加速收敛因子加入迭代公式中形成最终的反演迭代公式.模型试算结果证实了本文方法可有效减少直接迭代法的迭代次数,并且也适用于含噪重力数据反演.最后将该方法用于渭河盆地西安凹陷的实际资料处理之中,反演结果证实了该方法的实用性. 相似文献
11.
将研究区域划分成具有固定宽度的矩形网格，以网格密度和厚度作为模型参数，在此基础上形成重力异常的反演目标函数，计算出对模型参数的偏导数矩阵，然后采用脊回归法对重力异常进行反演而同时得到密度及其界面。以此方法对理论模型进行了反演试验。 相似文献
12.
联合反演是综合地球物理研究的重要定量解释手段.本文在总结和分析重力与地震资料联合反演的研究现状基础上,利用改进的全局寻优的快速模拟退火算法,实现了重力和地震资料的约束同步联合反演.针对性地设计了密度和速度界面不完全一致的模型,理论模型的试验说明了方法的效果和适用性.结合最近完成的广东徐闻地区实际资料的处理和解释,表明该方法可准确确定复杂构造物性界面的密度和速度结构,在该地区的油气勘探中发挥了作用.在先验信息约束下,该联合反演方法要明显优于单独的重力反演. 相似文献
13.
14.
本文用压缩质面来近似二维地质体。先由地面上测得的重力数据用矩阵方法反演此压缩质面各单元的面密度,然后从面密度与体密度差的关系求得各二维质体单元的厚度,进而得到各单元的近似地壳厚度。经正演校验和反复调整结果,使计算的重力异常值与实际测量值之残差小到满足要求,从而得到较准确的地壳底部界面。 由上地幔,玄武岩层和花岗岩层的密度差异及已算出的地壳厚度,从重力异常中分解出莫霍界面和康腊界面起伏所分别引起的重力异常。将后者同样用反演地壳厚度的压缩质面法进行计算,得到康腊界面。 文中以三种假想的地壳模型和一个实测剖面为例来检验本方法,并用其他确定地壳界面方法所得到的结果相比较,表明本文提出的压缩质面法结果较好。 相似文献
15.
We present an inversion technique based on the Marquardt algorithm to estimate the depth of a 2.5-D sedimentary basin in addition
to the regional gravity anomaly that is associated with the residual gravity anomaly, wherein the density contrast varies
parabolically with depth. Forward modeling is carried out through a derived analytical gravity expression of a 2.5-D vertical
prism. Inversion of a theoretical gravity anomaly with and without a regional gravity anomaly illustrates the procedure that
it is found to be insensitive to the regional gravity effect. Furthermore, the algorithm is exemplified with the gravity anomalies
of the derived density-depth model of the Godavari subbasin, India with a parabolic density profile resulting in a more consistent
geological model rather than a constant density profile. The main advantage of this method is that it works well even when
the profile of interpretation does not bisect the strike length of the sedimentary basin. 相似文献
16.