基于欧拉反褶积方法估算华北地区重力变化的等效源参数

利用华北地区2009—2014年绝对重力与相对重力多期重复观测资料, 得到不同时空尺度的华北区域重力场动态变化图像。 采用欧拉反褶积方法, 通过对理论模型试算, 获得最优反演参数, 对引起华北地区重力场变化的场源深度和空间分布规律进行了反演计算和分析。 结果表明: 当构造指数为1时, 适合对流动重力场变化数据进行反演。 实际资料的反演结果在2009—2014年间的场源位置集中于河套断裂带。 本文研究方法可以用于重力场反演和以场求源的定量研究。 本工作为流动重力变化信号的反演和解释提供一条新思路, 也为构建地震重力预报指标体系提供了定量依据。     

辽宁地区时变重力场源变化特征分析

地震重力分析通过研究时变重力场变化获取地球内部介质物性变化信息.采用贝叶斯重力网平差方法对辽宁地区2011-2014年共计四年7期的流动重力观测资料进行处理,对研究区重力观测网总体监测能力做出分析,选用欧拉反褶积对研究区的重力变化场源深度及空间分布特征进行反演和解释.通过反演计算发现,在2013年灯塔Ms5.1地震前沈...     

鄂尔多斯南缘地区重力变化场源特征

以鄂尔多斯南缘地区布格重力异常数据及2014—2017年相对重力观测数据为基础,采用欧拉反褶积方法,对引起重力变化的场源深度进行反演,并对空间分布规律予以解释。通过构建理论模型,反演得到最优模型参数,并对实际数据进行计算和分析。为了减弱和消除欧拉解的发散性,利用水平梯度滤波法优化反演结果。结果表明:①构造指数为1时,适合对相对重力数据进行反演,当滑动窗口选择4—8倍测点间距时,可以获得较为可靠的场源参数;②在主要断裂附近,不同时间尺度的场源呈较好的一致性特征;③引起重力变化的场源深度集中在10—30 km,与鄂尔多斯南缘地区的震源深度基本一致;④不同时间尺度的场源位置相对分散,表明引起重力变化的物质流动具有随机性和波动性。     

基于欧拉反褶积方法计算辽宁地区重力变化场源特征

以2012—2018年辽宁地区的流动重力观测资料为基础,结合实际地下地质体构造特征构建理论模型,应用欧拉反褶积法对重力变化场源深度及空间分布特征进行反演和解释。结果表明:在辽宁西部辽蒙交界和海城—岫岩断裂带附近场源位置集中,地下介质活动增强;欧拉反褶积方法适用于对流动重力观测资料的重力场反演及以场至源的定量研究应用。     

时变重力场球面模型反演算法和模拟实验

时变重力场是研究地球内部介质物性变化的重要手段.本文提出了一种适用于地面流动重力测量获得的时变重力信号的场源反演方法,该方法采用球坐标系下的六面体单元来模拟场源介质,适合大尺度地震流动重力测量数据的等效源模型构建.通过引入重力时变信号的一阶光滑先验条件,压制了时变重力信号中的短周期高频分量,可用于提取与地震孕育相关的长...     

云南鲁甸MS6.5地震震源区和周边三维密度结构及重力场变化

2014年8月3日云南鲁甸发生M_S6.5地震,震源机制显示该地震为左旋走滑型,断层破裂面以北西向为主,与昭通断裂带近乎垂直.本文基于三维重力反演技术,采用Crust1.0地壳结构模型和最新Moho面深度模型联合约束,反演得到了该地震震源区及周边的三维密度结构特征.反演结果表明：鲁甸地震震源区南北两侧块体运动受到昭通断裂不同程度阻隔,在地质时间尺度上,昭通断裂北段地壳物质运移的阻挡能力相比南段更强.通过进一步调查距离震中较近的流动重力测点的点值变化发现,从2014年3月至2014年6月大海村（P1）重力点值减小,减小幅度达到19 μGal,而位于昭通断裂东北侧的江底（P5）重力点值同期变化显示为增加,其他距离较远测点在震前重力变化不明显.综上结果对鲁甸地震的成因机制给出了初步解释.本文结果对于认识孕震区深部结构及地球物理场时空变化特征具有一定的意义.对识别与局部重力变化有关的地震前兆也同样具有参考意义.     

不同高度数据联合欧拉反褶积法研究

基于不同测量高度重力场及其梯度数据可同时对应同一场源并用于反演场源位置的分析原理,拓展不同高度场数据在欧拉反褶积法中的应用范围.首先,立足于对欧拉反褶积方法的理论研究基础,提出不同高度数据融合联合欧拉反演公式.其次,在理论模型上对多种高度数据联合反演做了测试分析计算,验证了不同高度场数据融合联合欧拉反褶积法能够改善位场解释中单一观测面数据计算带来的解的发散问题,收敛过程由此改善.最后,将本文方法应用于龙门山地区实际重力数据的解释,获得了研究区断裂分布特征.     

解释位场全张量数据的张量局部波数法及其与常规局部波数法的比较

全张量探测技术以其信息量大、精度高、干扰小等优点在地球物理领域中得到广泛应用.本文提出采用张量局部波数法来进行位场全张量数据的解释,首先给出了张量局部波数的定义,然后推导出利用张量局部波数法进行反演的基本公式.本文方法在进行张量数据反演时无需事先知道场源体的类型(构造指数)即可获得场源体的位置信息,且可根据位置参数对场源体的类型进行估计.通过理论模型证明张量局部波数法可以很好地完成位场全张量数据的反演工作,并将其与常规局部波数法进行对比,证明全张量局部波数法的反演结果更加准确,即使在测点分布不合理的情况下,张量局部波数法仍可以获得准确的结果.最后应用张量局部波数法对美国得克萨斯州实测重力数据进行了反演,其反演结果与已有的研究成果相一致.     

Vinton dome地区全张量重力异常的解释

Vinton dome地区位于美国墨西哥湾,发育着大量的盐丘构造,为油气的储藏创造了有利条件.地震勘探结果在盐丘区域存在一定的发散从而难以准确地确定盐丘的准确范围,难以为钻探工作提供可靠的信息.为了获得准确的盐丘形态,现今大多采用重力与地震数据的联合解释.重力梯度异常相对原始重力异常具有更高的分辨率,Bell Geospace公司于2011-2012年在Vinton地区进行了全张量重力测量.本文推导出重力张量异常进行界面正反演的公式,并通过理论模型试验证明了重力与张量异常联合反演可有有效地提高反演结果的准确性.根据先前地震结果利用实测重力张量异常反演获得了Vinton dome地区盐丘界面的变化,此外还利用密度反演方法圈定出盐丘的具体分布,结果显示两种方法所获得的盐丘顶面位置相一致,从而为下一步的钻探工作提供了有力的保障.     

2014年8月3日云南鲁甸M6.5地震破裂方向性研究

2014年8月3日云南鲁甸发生了M6.5地震,造成了严重的人员伤亡与经济损失.基于近震波形CAP方法反演得到的震源机制解表明此次地震主要为走滑型地震,两个断层节面的走向分别为70°及160°左右.由于震源区构造复杂,分布着不同尺度及走向的断层,仅根据地震断层面解难以确认发生破裂的实际断层.利用主震和参考地震之间的P波到时差及CAP反演过程中得到的时移信息,测量得到了鲁甸地震质心位置与破裂起始位置的差异,确定了地震断层面解中走向为160°的节面为真实破裂面.本文得到的鲁甸地震破裂方向性可为进一步研究该震发震成因提供参考.     

2014年云南鲁甸6.5级地震前重力变化特征与3维反演

利用西昌测网2012-2014年6期流动重力观测资料,分析了区域重力场差分动态变化和累积动态变化特征,以及近震源区测点的时序变化曲线,认为鲁甸6.5级地震前2年并没有出现明显的趋势性累积异常,而可能的前兆信号与断层滑动面两侧的测点2014-03-2014-06时段的重力点值差异性变化相关.针对此异常采用等效源模型反演至震源深度范围内的密度值变化在10-5g/cm3量级,推测该质量源在短期内的变化可能与地壳内部介质孔隙流体的充填或运移相关,在短期构造运动等动力条件下,孕震区可能由于流体物质运移、充填并诱发断层滑动,导致地震的发生.     

川滇地区重力场动态变化及其强震危险含义

基于川滇地区2011—2014年的重力复测资料,系统分析了区域重力场时-空动态变化及其与2012年云南彝良MS5.7、2013年四川芦山MS7.0、2014年云南鲁甸MS6.5和四川康定MS6.3地震发生的关系.结合GPS、水准观测成果和区域地质构造动力环境,进一步研究了区域重力场变化的时空分布特征及其机理,讨论了近期区域重力场动态变化的强震危险含义.结果表明:1重力变化与川滇地区断裂构造活动存在密切空间联系,重力变化较好地反映了伴随活动断层的物质迁移和构造变形引起的地表重力变化效应.2重力资料对测区内2012年以来发生的4次MS5.7以上强震均有较好反映,地震前震中区及其附近观测到明显的区域性重力异常及重力变化高梯度带,可能是地震孕育过程中观测到的重力前兆信息.3区域重力场动态演化大体反映了青藏高原物质东流的动态效应,龙门山断裂带地壳受挤压隆起、面压缩率和重力上升变化的特征最为显著.4重力场的空间分布及其随时间变化与地壳垂直与水平运动及地质构造活动等观测结果有一定的对应关系,强震易发生在重力变化四象限分布中心地带或正、负异常区过渡的高梯度带上,研究区的一些重力异常部位仍存在中-长期大震危险背景.     

2017年九寨沟地震前重力场异常特征提取

2017年8月8日九寨沟发生7.0级地震,此次地震发生在塔藏断裂、岷江断裂和虎牙断裂附近,震源机制为走滑型。重力场静态异常和动态异常变化可为研究深地壳结构性质提供重要的物理场信息。本文利用WGM计算的布格重力异常数据及实测流动重力数据进行小波多尺度分解,分离得到不同深度场源特征的地壳介质横向不均匀性。结果表明,WGM布格重力异常数据和实测重力数据的小波多尺度结果较为一致,通过对比分析发现,九寨沟地震发生在区域重力场的四象限分布中心位置,不同时空尺度的重力场变化对于深入认识潜在地震危险性趋势具有一定的科学意义。     

基于分频波形拟合方法的鲁甸M_S6.5地震序列震源机制解研究

2014年鲁甸M_S6.5地震造成了极大的人员伤亡和经济损失,其孕震机理和动力学成因为国内外学者广泛关注.由于目前的反演主要采用单一的反演方法进行处理,导致结果的不确定性较大,使得该地震的发震构造及孕震机理等问题仍存在较大争议.针对目前研究存在的不足,本文基于近震波形求解震源机制解的gCAP(general Cut And Paste)方法和全波形拟合方法,分别反演了鲁甸地震主震及其9次余震的震源机制解,以此评价近震数据反演震源机制解结果的可靠性与稳定性.结果表明数据方位角覆盖对gCAP解的稳定性有较大影响.针对全波形拟合方法反演过程中,低频信号稳定、高频信号解析度高的特点,采用0.01～0.05Hz和0.01～0.2Hz分频段波形拟合思路:在低频段剔除拟合差的数据,进一步在高频段进行高解析度波形拟合,从而获得主震可靠稳定的震源机制解.研究结果表明鲁甸地震主震倾角为76°～83°、滑动角为-157°～-164°,为一次高倾角走滑型地震.两种方法获得的余震震源机制解比较一致,验证了结果的可靠性.研究结果显示,鲁甸地震的高倾角破裂特征,会导致其能量释放快速且完整,这可能是导致鲁甸地震地面破坏加重及缺乏较大震级余震的主要原因.     

欧拉反褶积法在航空重力勘探中的应用

在论述欧拉反褶积基础理论与方法的基础上,以我国某海域航空重力资料为例,对该海域内地质异常体(侵入岩、断鼻构造)及断裂构造等进行了欧拉反褶积反演计算,获得了场源深度等信息,确定了场源位置及边界,并与航空布格重力异常、垂向导数、上延、水平梯度模、剩余异常以及重磁联合反演结果等资料相互验证,取得了很好的应用效果.基于该地区的反演效果,证明欧拉反褶积法可应用于航空重力勘探中,可为其他航空重力勘探区计算场源位置及深度信息提供可靠并且准确的研究方法,从而达到辅助航空重力勘探测量的效果.     

滇东北地区三维电性构造研究

滇东北地区位处不同构造分区的交界地带,发育有多样的断裂构造,其中昭通—鲁甸断裂和小江断裂带的部分区段均具备积累高应力的条件,具有较强的孕震潜力.本文基于覆盖巧家、东川、鲁甸等地区采集的区域大地电磁阵列观测数据,通过大地电磁三维反演,获得了该地区的三维电性结构的水平和垂直切片,完整地揭示了滇东北地区的深部电性结构特征.反演结果表明,昭通—鲁甸断裂在10 km深度以内发育有多个局部高阻区,小江断裂带的巧家—东川段的东侧、东西两支交汇处在上地壳均表现为高阻,这些区域均展示出更强的孕震能力.根据低阻异常的分布,推测昭通—鲁甸断裂和小江断裂下方的低阻构造在中下地壳的多个位置发生连通.本文将研究区的地表热流、岩石圈有效弹性厚度、地震与电性结构结合分析,推测在滇东北深部低阻层可能存在下地壳流,并且下地壳流的塑性流动为鲁甸的地震活动提供了深部动力来源.     

基于时变重力数据的川滇地区场源模型反演和解释

时变重力信号能反映地球内部介质迁移引起的不均衡变化,本文利用三维欧拉反褶积方法,通过建立合适的场源模型获得最佳反演参数,并采用水平梯度滤波法消除部分虚假发散解,优化反演结果.在此基础上,利用川滇地区2015-2019年流动重力重复观测资料,分析了不同时空尺度的区域重力场动态变化特征,开展对川滇地区重力变化场源特征的定性...     

2014年云南鲁甸6.5级、永善5.0级地震前岩石圈磁场局部异常特征分析

利用2012 ～2014年南北地震带3期流动地磁矢量测量的准确可靠资料,计算得到了2012 ～2013年和2013～ 2014年相邻期岩石圈磁场的分布.分析结果表明,2014年8月3日鲁甸6.5级地震、8月17日永善5.0级地震前,川滇交界东部均存在岩石圈磁场的异常变化.同时,异常持续时间、空间范围与震级大小有一定关系,鲁甸地震前岩石圈磁场异常持续时间长、空间范围广,而永善地震则反之.     

2013年芦山Ms7.0地震前重力场源变化分析

陆地重力观测获取的时变重力信号是研究地壳内部孕震区介质变化的重要依据。利用贝叶斯平差方法处理了川西地区2010—2012年的流动重力观测资料,发现相对重力仪的格值系数在两年内变化可达2×10^-4,而且仪器存在较小的非线性漂移率;重力场累积动态变化图像显示,2013年芦山MS7.0地震发生在重力梯度带零值线与龙门山断裂带交会部位附近。进一步采用等效源方法划定了测网内以雅安、康定和石棉为中心的场源分辨能力较高的区域,并反演了地震前的地壳视密度变化。场源视密度动态变化结果表明:2013年芦山MS7.0地震前2年,在龙门山断裂带的南侧和鲜水河断裂南段的东侧区域形成了范围较大的等效视密度增加区,而在小金至康定一带和都江堰周边地区也存在小规模的视密度增加。本文基于陆地时变重力测量结果反演得到芦山地震前的地壳内部等效视密度变化,可为研究壳内深部孕震环境变化特征提供参考。     

青藏高原东南缘多尺度重力场变化特征及孕震机理分析

文中基于青藏高原东南缘2013—2016年的流动重力观测资料,获得了鲁甸和景谷地震前后不同时间尺度下区域重力场的动态演化趋势和异常特征,利用小波多尺度分解分析重力场变化与构造运动、物质密度变迁和强震孕育的关系和规律,对青藏高原东南缘深部物质变化、动力学过程及地震孕育机理进行了探讨。结果表明:1)鲁甸、景谷震前震源区重力场变化呈现与断裂构造走向一致的剧烈正、负异常过渡带和梯度带特征,显示了孕震期间震源区强烈的地壳运动和深部物质变化过程; 2) 2013年9月—2014年4月重力场变化的小波多尺度分解结果显示,青藏高原东南缘地壳、上地幔不同深度和尺度的重力场变化与地震分布及活动断裂带位置相关性显著,表明研究区地震孕育与断裂运动以及地壳、上地幔内物质密度分布变化关系密切,可能受地壳、上地幔物质运移的复杂深部动力学过程影响; 3)基于青藏高原东南缘地壳、上地幔物质运移的动力学过程,对强震多发生在重力场变化的正、负异常过渡带和梯度带附近的特征进行了初步的解释。文中的研究结果对地球构造运动、地震机理研究具有一定的参考价值。