约束变密度界面反演方法

三维密度界面反演具有严重的多解性,已有方法较少利用已知信息约束反演过程,导致界面反演结果可信度降低.本文在反演过程中,利用已有地质信息作为约束条件,有效提高了三维界面反演结果的准确性.该方法具有如下特点:1)利用指数变密度模型,通过已知密度分布计算模型参数,使其更加接近实际密度分布;2)引入已知深度点约束,提高了反演结果的准确性;3)引入深度加权函数,纠正界面畸变,使其适用于界面起伏较大的情况;4)频率域正演与空间域迭代反演相结合,在提高计算速度的同时保证反演收敛.通过模型检验,证实了方法的有效性,并将该方法应用于中蒙边境地区东段莫霍面深度反演中,效果良好.     

改进的约束变密度界面反演策略及其应用

计算密度分界面的起伏变化在区域地质构造研究和石油矿产资源勘探中具有重要意义.已有密度界面反演方法更多侧重约束变密度界面反演算法,而对约束信息的准确性、研究区横向密度变化往往考虑不足,影响了最终反演结果的可信度.本文在变密度界面正反演算法基础之上,结合实际需求,提出已知深度信息约束下的变密度界面反演策略.该策略主要包括变密度约束反演算法、已知深度约束信息校验和分区变密度模型三个方面.其中反演算法提供了带已知深度约束信息的密度界面迭代反演方法;约束信息校验用于评估约束信息精度,通过调差降低约束信息的系统误差;在反演过程中引入水平密度分区以应对不同构造背景密度界面模型,提高反演结果的可信度.最后将本策略应用于南海莫霍面深度反演计算中,结果显示借助已知约束信息,利用分区密度模型能够获得更为可信的深度反演结果,验证了该策略的正确性.

     

改进的约束变密度界面反演策略及其应用

计算密度分界面的起伏变化在区域地质构造研究和石油矿产资源勘探中具有重要意义.已有密度界面反演方法更多侧重约束变密度界面反演算法,而对约束信息的准确性、研究区横向密度变化往往考虑不足,影响了最终反演结果的可信度.本文在变密度界面正反演算法基础之上,结合实际需求,提出已知深度信息约束下的变密度界面反演策略.该策略主要包括变密度约束反演算法、已知深度约束信息校验和分区变密度模型三个方面.其中反演算法提供了带已知深度约束信息的密度界面迭代反演方法;约束信息校验用于评估约束信息精度,通过调差降低约束信息的系统误差;在反演过程中引入水平密度分区以应对不同构造背景密度界面模型,提高反演结果的可信度.最后将本策略应用于南海莫霍面深度反演计算中,结果显示借助已知约束信息,利用分区密度模型能够获得更为可信的深度反演结果,验证了该策略的正确性.     

变密度界面模型重力异常反演的B样条函数法

重力异常反演是一个求解第一类非线性积分方程的问题。当三维界面起伏满足│Δｈ│〈０．４１４ｈ０和二维界面起伏满足－ｈ０〈Δｈ〈０．４１４ｈ０时，积分方程的被积函数可Δｈ的级级数形式。     

约束变密度界面反演法在渤海深部结构研究中的应用

在分析研究约束变密度界面反演基础理论及其方法应用的基础上,本文以我国渤海深部结构为研究对象,利用变密度反演法,对渤海海域主要坳陷进行沉积厚度的反演及其与断层等的综合解释工作,而后对地壳与地幔之间存在的重要分界面-莫霍面,基于已有的对莫霍面信息有较为系统、整体性反应的研究区重力数据,对其进行变密度反演计算,获得了研究区莫霍面的深度分布情况.通过与已知的地质地球物理资料,主要地震OBS测深等前人研究成果相互对比验证,取得了较为良好的效果.在此基础上,进一步证明约束变密度界面反演法在沉积界面反演、深部构造研究等方面,可以为其他勘探区域的研究提供较为准确可靠的方法参考,从而丰富发展地球物理反演工作.     

利用重力资料反演三维密度界面的质面系数法

文中给出一深度为h的水平矩形质面在地面产生的重力公式,对比了此公式算出的重力异常与三维长方体公式算出的重力异常,得出:当三维地质体的厚度与其埋深之比小于0.2—0.4时,二者的重力效应几乎相等,它们之间最大偏差在1—3％以内。用给出的重力公式进一步导出了重力反演三维密度界面的质面系数法。用理想模式检验,结果表明,本文给出的方法收敛快,反演误差较小,计算效果较好。最后用一实例,验证了方法的适用性。     

重力密度界面反演方法研究进展

重力位场的界面反演是位场处理解释中的重要问题.本文针对重力密度界面的反演问题,较全面的介绍了一些有代表性的方法,包括空间域反演方法和频率域反演方法.重点阐述了当前占主流的频率域中的Parker-Oldenburg迭代方法,在此基础上讨论了重力密度界面反演目前存在的问题和以后的发展方向.     

利用重力资料反演三维密度界面的质面系数法

文中给出一深度为h的水平矩形质面在地面产生的重力公式,对比了此公式算出的重力异常与三维长方体公式算出的重力异常,得出:当三维地质体的厚度与其埋深之比小于0.2-0.4时,二者的重力效应几乎相等,它们之间最大偏差在1-3％以内。用给出的重力公式进一步导出了重力反演三维密度界面的质面系数法。用理想模式检验,结果表明,本文给出的方法收敛快,反演误差较小,计算效果较好。最后用一实例,验证了方法的适用性。     

用四维反演方法研究地震活动的深部环境

导出了用重力和水准重复测量资料反演地壳界面密度和垂直形变随时间变化的四维反演公式，并对甘肃河西地区（包括祁连山地区）和滇西北地区的重力和水准重复测量资料进行了计算。根据计算结果，分析了2个地区莫霍西密度变化速率和垂直形变速率与地震的关系。结果表明：甘肃河西地区和滇西北地区莫霍面密度大规模减小和莫霍面大范围隆起可作为1年尺度的地震预报的深部判据。地震往往发生在莫霍面隆起区与沉降区、密度减小区与增长区的过渡地带。     

Gravity field separation, density inversion and crustal tectonics in Kang-Dian region

IntroductionTethys-Himalya tectonic belt extends from west to east across Europe, Asia and Africa, and stops at the longitude of 104(E in China. The famous North-south active tectonic belt located in China is a part of China-Mongolia central-axis strong earthquake belt and the huge crustal thickness undulation belt. The tectonic belt is a strong extrusion, tension and shear region. It is synthesized by topography gradient belt, faults, Cenozoic basin and strong earthquake belt on north-so…     

重力及重力梯度联合反演青藏高原及邻区岩石圈三维密度结构

本文利用GOCE L2观测重力梯度的五个独立分量（Txx,Tzz,Txy,Txz,Tyz）,联合EGM2008地球重力场模型计算垂直重力,反演计算了青藏高原及邻区0~120 km深度岩石圈三维密度结构.将经过低阶项改正、地形效应改正、沉积层界面起伏效应改正得到的剩余重力及重力梯度异常值作为观测值,以改正剩余量归一化权重作为观测权重,基于Tikhonov正则化理论建立反演目标函数.反演过程中,利用地震层析S波速度转换密度作为初始约束,通过非等权最小二乘迭代法计算得到最终反演密度.反演结果表明：（1）40 km深度,青藏高原内部为中地壳,表现为低密度,邻区为中下地壳,表现为高密度.青藏高原内部中地壳强低密度层主要分布在高原边界.其成因是印度板块俯冲和周围坚硬块体阻挡作用导致在高原边界形成的高应变积累闭锁区,为壳内低密度软弱物质的形成提供了条件.（2）80 km深度,青藏高原上地幔顶部显示出低密度的特征.高原内部东、中、西密度特征差异明显,低密度以95°E为中心线呈东西对称分布.以班公—怒江缝合带为中心,在拉萨块体和羌塘块体内从北向南出现了"低-高-低"的密度分布起伏特征.该特征与GRACE得到的莫霍面起伏特征一致,结合大地构造结果,这种起伏特征验证了印度、羌塘块体从南北两侧分别向喜马拉雅、拉萨地块挤入的双向俯冲模式.（3）四川盆地和鄂尔多斯盆地内,地壳高密度异常较地震波速异常明显偏低,表明古老的四川盆地和鄂尔多斯盆地比想象中更冷、更坚硬.塔里木盆地和柴达木盆地内壳、幔高密度的结构特征,对应地幔物质上涌.     

空-地-井重力异常正则化协同密度反演方法

围绕综合利用空-地-井多源重力异常联合反演提升精度的目标,提出正则化协同密度反演方法,从而有效利用多维数据的横向和纵向变化特征提高反演精度,且无需先验信息的约束.此外,还提出利用奇异值谱和深度分辨率工具评价航空数据观测高度和层数对反演分辨率的影响.通过理论模型试验本文方法在空-地、地-井等不同情况下的应用效果,结果表明空-地-井重力异常正则化协同密度反演方法能获得更高分辨率、高精度的反演结果,且证明不同的钻孔位置对反演结果有较大的影响,从而可指导实际的空-地-井联合勘探.最后,利用文顿盐丘实际勘探数据进行空地和空地井重力数据协同反演,反演结果垂向分辨率明显优于地面观测数据反演结果,盐盖的顶面及中心埋深与前人研究和地质资料解释相吻合,验证了方法正确性及实用性,为推进我国空-地-井立体重磁勘探提供了重要的技术手段.

     

青藏高原地壳三维密度结构的重力反演

在阻尼最小二乘法的基础上,提出了利用重力进行三维密度反演的选权拟合法,并利用密度分布的特点构造了约束权矩阵.为了改善单一重力反演的非唯一性,利用地震面波层析成像资料进行约束,得到了青藏高原不同深度高分辨率三维密度异常分布图.反演结果显示,青藏高原地壳内密度很不均匀,不同块体具有不同的密度异常特征,拉萨和羌塘地体为低的负密度异常区,松潘-甘孜地体东部存在较低的负密度异常区,塔里木盆地存在较高的正密度异常区.从中部到东部,上地壳物质存在差异,而下地壳物质存在一定连续性.     

基于位场分离与延拓的视密度反演

由于目前计算机内存和速度的限制,在对大面积重力资料进行三维密度反演时,已有的反演方法很难奏效.文中提出了一种基于位场分离与延拓的视密度快速反演方法：首先应用场分离的切割法对平面上的重力场进行不同深度层源的切割分离;然后运用大深度向下延拓方法将各层的场延拓至相应的深度;最后反演出各深度层的密度.反演得到的密度是各深度层密度的近似分布,称为视密度反演.本反演方法克服了传统已有方法计算时间特别长、解稳定性差的缺点,在主频1.99 GHz的微机上,反演128×128×10个密度值的计算时间小于20 s,理论模型和实际资料的应用表明反演具有较好的效果.     

垂直重力梯度反演Moho面的频谱域公式及其应用

通过求解引力相等原则下的Fredholm积分方程,可以得到不规则单一密度界面(Moho面)的起伏.本文充分参考了前人的理论研究,推导出扰动垂直重力梯度确定Moho面深度的频谱域表达式,该式具有二次项迭代精度.运用此公式进行了全球Moho面的恢复计算,并将该结果与CRUST1.0模型和GEMMA Moho模型进行了对比和验证.     

地震走时与重力数据的联合反演

对目前的地震走时与重力数据联合反演发展进行了简要的概述。首先,从两类数据单独反演的基本原理出发,对两者进行了分析比较,并由两者参数间存在的相互关系,可将两类数据归并、统一到一个系统,合理地进行联合反演。其次,根据目前国内外的相关研究,介绍了当前联合反演的两种算法:依次反演和同时反演。最后,对联合反演中存在的问题进行了讨论。