首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
利用非线性规划研究的最新成果, 设计了一种全新的震源破裂过程的反演方法, 同时反演震源破裂的时间图像和滑动在破裂面上的分布. 该方法与目前震源过程破裂反演计算所用的两类主要方法(线性矩阵方法和全局搜索寻优方法)相比, 在线性矩阵方法对初始模型依赖大、解不稳定和全局搜索寻优方法在高维解空间寻优效率不高、易出现伪解等方面有较显著改善. 数字试验的测试表明, 本方法计算效率高, 反演结果稳定、可靠, 可有效应用于实际地震震源过程的反演研究中.  相似文献   

2.
白超英  赵瑞  李忠生 《地震学报》2009,31(4):385-395
讨论了全局选择震源初始位置下的矩阵反演求取全局解的问题.与流行的全局优化搜寻直接定位方法相比,该反演算法是一种行之有效的定位方法,具有定位精度高、计算时间少以及对噪声数据不敏感等特点.其突出的优点是在不增加计算难度和计算时间的前提下,用矩阵反演的方法确保得到全局最小值解,可适应于地震早期预警,海啸早期预警,以及大震速报等实际工作.   相似文献   

3.
微震监测定位技术的精度及效率是微震监测中的核心技术问题,震源扫描算法作为一种稳定、可靠的方法被广泛用于微震监测定位中,但是在多事件、多参数反演中,其反演效率严重制约其应用.本文将震源扫描算法与DIRECT算法思想结合,提出基于DIRECT算法的微震快速网格震源搜索定位方法,即在求解震源扫描算法目标函数最优解时,采用DIRECT搜索策略,无需划分网格大小,直接搜索目标函数最优解.通过地面及井中微震监测模型测试表明,基于DIRECT算法的微震快速网格震源搜索定位方法比传统的震源扫描算法搜索次数少、计算耗时短、计算精度高,尤其适合大空间内多事件、多参数的震源反演,为水力压裂微震监测过程中大量数据的实时、快速、准确定位提供了可能.  相似文献   

4.
介绍了震源破裂过程反演方法研究的必要性及意义,简要回顾了震源破裂过程反演研究的历史,详细阐述了现在几种主要的震源破裂过程反演方法,分析了它们各自的优势和不足。展望了震源破裂反演的发展方向,指出地震波数据同其他多类数据尤其是大地测量数据(InSAR和GPS)联合进行反演是未来震源破裂过程反演的发展趋势。  相似文献   

5.
王文萍  王庆良 《地震学报》1999,21(3):285-290
中国地震局第二地形变监测中心已监测到与1990年4月26日青海共和7.0级地震有关的大幅度隆起垂直形变资料.基于这些形变资料,本研究首先利用全局搜索性能优良的遗传算法和弹性半空间单一矩形均匀滑动模型,对该地震的震源位错参数进行了初步反演;然后以遗传算法的反演结果为初始值,利用拟牛顿阻尼最小二乘方法对该地震的震源位错参数进行了更精细的反演.结果表明,共和7.0级地震的发震断层为一走向NW60左右、倾向SW、倾角37左右的右旋走滑逆冲断层,同震断层破裂长度约37 km,宽度约6 km,错动幅度约2.7 m.反演结果还同时表明,遗传算法与最小二乘方法相结合是一种行之有效的联合反演方法,它既可以保证实现全局最优解,防止最小二乘方法的局部极值解问题,又可以有效地解决遗传算法达到全局最优解区域后收敛速度过慢的问题.   相似文献   

6.
田宵  汪明军  张雄  张伟  周立 《中国地震》2021,37(2):452-462
微地震事件的空间分布可以用来监测水力压裂过程中裂缝的发育情况。因此,震源定位是微震监测中重要的环节。震源定位依赖准确的速度模型,而震源位置和速度模型的耦合易导致线性迭代的同时反演方法陷入局部极小值。邻近算法作为一种非线性全局优化算法,能够最大程度地避免陷入局部最优解。本文将邻近算法应用于单井监测的微震定位和一维速度模型同时反演,首先利用邻近算法搜索一维速度模型,再使用网格搜索方法进行震源定位,并根据定位的走时残差产生新的速度模型,最后通过若干次迭代使其收敛到最优解。理论和实际数据结果均表明该方法能够避免局部最优解,得到较为可靠的震源位置和一维速度模型。  相似文献   

7.
主要讨论超松弛迭代法反演震源破裂过程解的稳定性。在精确反演的基础上,分别对数据扰动和模型扰动进行分析,用MATLAB作图进行对比并得出最后的结论:超松弛迭代法反演震源破裂过程的结果是可靠的。  相似文献   

8.
众所周知,破裂过程的反演依赖于震源位置和断层面参数.然而,在一些台站稀疏地区,由于震源位置的精度不高,在进行破裂过程反演时,断层面上的初始破裂点的位置难以得到有效确定.  相似文献   

9.
时间域内有限地震断层的反演问题   总被引:9,自引:18,他引:9       下载免费PDF全文
为了解决有限地震断层破裂过程研究中的正演和反演问题,本文提出:(1)用广义反射、透射系数矩阵和广义射线理论计算上地幔及远场范围体波理论地震图的混合方法.同时采用格林函数的一种新的组合方式,可节约1/3计算机内存.(2)采用基于模拟退火的热浴算法进行全局搜索反演,不仅提高了在计算过程中新模型被接收的可能性,而且使得在时间域内目标函数的计算速度提高了近二个数量级.对1994年台湾海峡地震震源破裂过程的研究结果表明,这次地震的破裂过程相对比较简单.在前4.0s内地震断层的破裂相当剧烈,错动距离的平均值随着破裂时间增加而迅速减少.开始时破裂速度较高,前4.0s的平均值为3.1km/s,之后逐渐降至2.2km/s.错动方向随破裂的变化不大.  相似文献   

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

11.
通过开展野外调查、震源深度校核、震源机制解反演等多项工作,综合研究了2011年7月威海地震序列.利用P波初动方法反演震群震源机制解,结果表明本次震群的破裂方式均为走滑型,节面Ⅰ(走向北东向)是最可能的破裂面.结合已有的研究成果,初步推断震群的发生受北东向断裂的控制.通过识别地震波形的Pn震相和sPn震相,计算得到地震震源深度,推测本次地震发生在上地壳.由于震源深度较浅,地震烈度调查工作表明震中最大烈度能达到Ⅴ度.  相似文献   

12.
在全面学习和了解震源破裂过程研究历史的基础上,从地震位移表示定理出发,梳理推导了地震震源时空破裂过程反演研究的基本原理、总结了地震破裂过程反演研究中惯用的约束条件、典型的视震源时间函数提取方法和已有的多种震源破裂过程反演方法,并在已有研究的基础上,尝试性地提出了两种新方法:等矩反投影方法和子事件成像方法。首先,为了提高区域(中强)地震震源过程的时空分辨能力,借助于2013年芦山M_W6.6地震和2014年康定M_W6.3以及M_W5.8地震等实际震例,我们探讨了利用经验格林函数技术从区域地震记录中提取视震源时间函数、进而利用不同视震源时间函数分析和反演中强地震震源过程的效果。其次,为了快速获取震源破裂过程信息,用于震后快速响应以及早期预警,在充分吸收了传统高频反投影方法的优点的基础上,提出了等矩反投影方法,用于对震源破裂过程激发的低频能量源的直接成像。通过多组数值实验测试了新发展的等矩反投影方法的可行性,并将该方法应用于2014年11月在我国四川省康定发生的两次(M_W6.3和M_W5.8)中强地震的震源破裂过程研究以检验新方法的实用性。结果表明,等矩反投影方法对于提取低频能量迁移的主要特征是有效的,是对高频能量反投影方法的必要补充。然后,基于等矩反投影方法的基本思想,并借鉴于传统的震源破裂过程反演方法,新发展了用于对震源破裂时间过程和同震滑动量分布进行直接成像的子事件成像方法。相较于震源破裂过程研究中传统的反演方法而言,新方法不需要事先经验性地给定最大破裂速度以及子断层滑动上升时间,故而在对震源破裂过程进行成像时,允许破裂速度发生变化,并且允许子断层发生多次破裂。同样,利用数值实验测试了新方法的可靠性,且为了检验该方法的实用性,将子事件成像方法用于2014年11月在我们四川省康定县发生的两次(M_W6.3和M_W5.8)中强地震的震源破裂过程研究。结果表明,子事件成像方法能够获取震源过程的主要破裂信息,但与视震源时间函数反演得到的结果相比,破裂尺度略小。最后,将多种方法应用于2015年4月发生的尼泊尔M_W7.8地震及其强余震,获得了这些地震的震源破裂过程图像,对比讨论了这些方法在揭示震源时空破裂过程复杂性信息方面各自的优缺点。具体讲,首先利用远场P波视震源时间函数分析了主震及其强余震的破裂方向,并利用P波波形资料和InSAR资料联合反演了这些事件的震源破裂过程;然后利用等矩反投影方法和子事件成像方法分析了这些事件的低频能量的时空迁移过程和震源破裂过程;最后对利用不同方法得到的结果进行了对比讨论。对于M_W7.8主震而言,联合反演结果表明其以单侧破裂为主,向南东东方向拓展约150km,并沿断层倾向方向向深部拓展约70km,这与利用等矩反投影方法和子事件成像方法直接成像结果得到的主要破裂特征较为一致,但破裂空间拓展尺度要略大,这可能是由于联合反演所给定的3km/s的最大破裂速度要大于直接成像获取的主震破裂速度所致;对于M_W6.7强余震而言,直接波形反演方法与等矩反投影方法和子事件成像方法得到的结果一致表明,该事件破裂主要沿断层倾向方向向深部拓展,并且利用直接波形反演方法搜索到的此次事件最大破裂速度为4km/s,与等矩反投影方法直接成像结果中获取的3.9km/s的破裂速度较为一致,均表明此次强余震可能是超剪破裂事件;对于M_W7.3强余震而言,由于融合了同震InSAR资料的缘故,联合反演结果表明,此次事件以朝向南—东南方向破裂为主,这与基于等矩反投影方法和子事件成像方法的直接成像结果中破裂以朝向南南西—南方向为主的主要破裂特征略有差别,并且直接成像结果表明此次事件破裂速度较慢。概而言之,本论文总结了震源破裂过程研究的发展历史,梳理推导了已有的方法;在此基础上提出了两种新方法,并通过必要的数值实验和对实际震例的应用,阐述了新方法的实用效果和优缺点。  相似文献   

13.
根据非线性规划研究的最新成果所设计的一种全新的震源破裂过程的反演方法, 用近场地震波观测资料反演了1999年9月21日发生在中国台湾省集集Mw7.6地震震源破裂过程.为了使反演中设置的断层模型与集集地震实际破裂面尽可能一致, 以尽可能减小由于断层模型设置的不确定性对震源破裂过程反演结果的影响, 设立的断层模型为与集集地震造成的主要地表破裂尽可能拟合的弯曲面模型.反演结果显示: (1) 集集地震震源的破裂大体持续了32 s, 其中主要破裂发生在第6~27 s间, 破裂主要集中发生在断层北段向东拐弯处.(2) 震源破裂以逆冲为主, 平均滑动角为64.5°, 与USGS, Harvard及CWB(台湾中央气象局)的结果相当.标量地震矩为7.76×1020牛顿米, 稍大于USGS和Harvard反演的标量地震矩.(3) 集集地震震源的破裂存在清晰的成核过程, 成核过程经历6 s后, 地震矩释放明显加速.起始破裂从断层南段开始, 10 s后破裂主要集中在断层北段发生.最后将反演结果与震后GPS观测结果进行了对照分析, 并对反演结果的科学意义进行了讨论.  相似文献   

14.
马尔科夫链蒙特卡洛方法(MCMC)是一种启发式的全局寻优算法,可以用来解决概率反演的问题.基于MCMC方法的反演不依赖于准确的初始模型,可以引入任意复杂的先验信息,通过对先验概率密度函数的采样来获得大量的后验概率分布样本,在寻找最优解的过程中可以跳出局部最优得到全局最优解.MCMC方法由于计算量巨大,应用难度较高,在地...  相似文献   

15.
根据非线性规划研究的最新成果所设计的一种全新的震源破裂过程的反演方法, 用近场地震波观测资料反演了1999年9月21日发生在中国台湾省集集Mw7.6地震震源破裂过程. 为了使反演中设置的断层模型与集集地震实际破裂面尽可能一致, 以尽可能减小由于断层模型设置的不确定性对震源破裂过程反演结果的影响, 设立的断层模型为与集集地震造成的主要地表破裂尽可能拟合的弯曲面模型. 反演结果显示: (1) 集集地震震源的破裂大体持续了32 s, 其中主要破裂发生在第6~27 s间, 破裂主要集中发生在断层北段向东拐弯处. (2) 震源破裂以逆冲为主, 平均滑动角为64.5°, 与USGS, Harvard及CWB(台湾中央气象局)的结果相当. 标量地震矩为7.76×1020牛顿米, 稍大于USGS和Harvard反演的标量地震矩. (3) 集集地震震源的破裂存在清晰的成核过程, 成核过程经历6 s后, 地震矩释放明显加速. 起始破裂从断层南段开始, 10 s后破裂主要集中在断层北段发生. 最后将反演结果与震后GPS观测结果进行了对照分析, 并对反演结果的科学意义进行了讨论.  相似文献   

16.
解析大地震的破裂模型是研究地震成因、破裂动力学机制和探究活动断层结构等研究方向的基础.在地震防灾、应急响应和危险性评估等实际问题中发挥着重要的作用.文章对近年来在大地震破裂模型反演的方法及应用上取得的进展进行综述和分析,结合北京大学地球物理专业相关研究人员在大地震破裂模型的快速反演、震源反投影、大地测量、地震波与形变观测数据的联合反演、复杂地球模型中破裂过程的反演及强地面运动模拟等方面的工作,分析不同研究方法的特长以及不足.提出从短期快速响应到长期科学问题的系统化研究思路,并对从震源模型反演到强地面运动预测的一体化自动应急系统,以及全动力学反演方法进行了分析和展望.  相似文献   

17.
震源机制解是地震学研究的重要组成部分, 它从一定程度上描述了震源的性质及其破裂过程, 对其正确、 全面地认识可深化理解地震的孕育、 发生和演化。 而地震波形自动反演系统的建立, 不仅实现了利用地震波形信息进行震源机制解的自动反演, 更重要的是加快了反演的速度, 充分体现了提取震源信息的实效性, 从而有力地帮助震情分析、 余震趋势预测、 震后救援等工作的展开。 本文系统地阐述了震源机制解反演方法, 总结了自动反演系统构建思路以及波形自动反演系统的实现, 为国内构建自动反演系统提供了可鉴的技术思路。  相似文献   

18.
用复合型法反演震源破裂过程的研究   总被引:5,自引:0,他引:5       下载免费PDF全文
本文探讨了用P波频谱研究不对称双侧破裂方式与正交破裂方式的方法,研究了用复合型法反演震源参数的可行性。 用这种方法反演了1966年3月8日邢台马兰6.8级地震及1966年3月22日东旺6.7级地震的震源参数。通过与用网格尝试法反演的结果及余震分布等进行比较,认为这种反演方法除了与过去的网格尝试法类似,可以得到最优拟合参数解外,还可以大大减少计算量,实现了网格尝试法难以实现的更复杂的震源模式的反演。 本文还用此方法研究了1962年3月19日广东河源6.1级地震及1969年7月26日广东阳江6.4级地震的破裂过程。研究结果表明:河源地震的破裂方式是正交三侧破裂。北西方向为主破裂,破裂长度8km,破裂速度0.6km/s. 南东方向破裂长度较短,约2km,破裂速度2.0km/s.正交侧指向南西,破裂长度3km,破裂速度0.5km/s,这个结果与余震分布图所显示的图象很相象。阳江地震的破裂方式是不对称双侧破裂。主破裂方向为北东东,破裂长度14km,破裂速度0.8km/s;朝南西西向的破裂长度5km,破裂速度0.6km/s。   相似文献   

19.
岩石裂纹扩展微观机制声发射定量反演   总被引:1,自引:0,他引:1       下载免费PDF全文
岩石受载内部微裂纹扩展及其震源机制反演有助于认识宏观裂纹扩展过程的非线性断裂力学行为.借助声发射监测手段,本文建立了仅涉及微裂纹张开/闭合和剪切滑移的位移不连续震源模型,通过各位置处传感器耦合质量标定及点源远场P波矩张量反演获得了含预制裂纹砂岩受载过程的震源机制解及时变响应特征,在全局坐标系下分析了微裂纹的三种断裂力学行为.结果表明:在位移不连续模型中,震源矩张量特征值与试样泊松比之间必须满足特定约束条件,该约束条件下的优化问题可采用拉格朗日乘子和Levenberg-Marquardt迭代法求解;受载砂岩裂纹扩展过程中,声发射震源以剪切滑移机制占优,微裂纹空间取向及运动方向与试样宏观主裂纹夹角平均值分别为40.9°和17.7°;对微裂纹体积分解表明岩体微观破裂机制以沿X方向I型张开为主,而沿Y方向的II型断裂滑移方向与试样全局变形方向相一致,由于试样内部晶粒分布非均质性造成了少量沿Z方向的III型平面外微裂纹滑移行为;受载砂岩裂纹扩展过程中微裂纹模式角与震源极性值变化趋势一致、利用震源震级评估的局部应力降值与实验观测结果相吻合,这两者均表明了位移不连续模型在震源机制定量反演中的适用性.  相似文献   

20.
应力场在地球动力学研究方面起着至关重要的作用。本研究给出了由震源机制解数据反演应力张量及其置信范围的网格搜索法。我们使用人工合成震源机制解数据对该算法进行了测试,人工合成的震源机制解数据由拉张型、挤压型和走滑型3种类型应力场生成,并对生成的震源机制解数据添加不同水平的噪声。测试结果表明,相对于传统方法,本研究给出的方法反演出的应力场精度得到显著提高。本文算法主要有以下3个优点:(1)由精细的网格搜索获得应力张量的全局最优解,搜索网格为1°×1°×1°×0.01,避免了搜索结果落入局部最小值;(2)在反演应力张量的过程中,本方法可以考虑震源机制解数据精度的不同,并赋予不同的权重;(3)由F检验给出应力张量结果的置信范围。我们将该方法应用于中国、越南和老挝的边界地区,获得该地区应力场的主压应力方向为SSE—NNW,主张应力方向为NEE—SWW。应力比R为0.6,这意味着应力张量的特征值接近等差数列。该地区应力场与奠边—琅勃拉邦弧形断裂的左旋走滑性质一致,所得结果对该地区的地质动力学过程研究有一定的帮助。  相似文献   

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

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