首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
The source process of February 3, 1996 Lijiang earthquake in Yunnan was studied by body waveform inversion using teleseismic data from IRIS. Two normal double-couple subevents with different strikes were obtained. The difference of the onset time between these two subevents, which are 15 km apart in space, is 7 s. The total seismic moment is 3.81 × 1018 Nm (M w=6.3). The total fault area S is about 720 km2 from the aftershock data and the average dislocation is about ū=0.18 m. Considering both the result of inversion and tectonic environment around the source, the first rupture might result from the extension along the NNW directed Zhongdian-Yongsheng fault belt where an earthquake of M=6.4 occurred in 1966. Then, the second started along the NE directed the eastern foot of Snow Mountain fault where rupture seemed to be able to propagate more easily.  相似文献   

2.

高精度的成像需要准确的速度场信息,波形反演被认为是目前具有最高分辨率的速度反演方法之一.计算效率是目前全波形反演需要考虑的一个主要问题.为了很好地解决计算效率的问题,本文引入了一种高效的无串扰同时源反演方法,并详细介绍了其原理与计算流程.同时,基于此多震源同时反演方法,本文拓展出实波场反演、虚波场反演及复波场反演的三种反演策略,进而分析复波场的实部和虚部对全波形反演的影响.实验表明,相比于复波场全波形反演,无论是实波场反演还是虚波场反演,反演分辨率有所降低;相比于实波场反演,虚波场反演对初始模型的依赖性较小,目标函数的非线性较弱;最后,通过使用组合反演策略,即初始阶段采用虚波场反演,中后期阶段采用复波场反演,不仅可以降低反演的非线性,而且能够保证高精度建模.

  相似文献   

3.
震源参数反演及精度评定的Bootstrap方法   总被引:1,自引:0,他引:1  

在震源参数反演理论研究中,地表形变与震源参数之间为复杂多维的非线性关系,针对传统泰勒级数展开的精度评定方法可能无法适用于震源参数的精度评定问题,本文将Bootstrap方法引入到震源参数非线性反演及精度评定研究中.通过对GPS地表形变观测数据实施Bootstrap重采样获取自助样本,使用遗传算法(Genetic Algorithm,GA)搜索震源参数,设计并给出了震源参数精度评定的Bootstrap方法计算流程.将本文方法用于6个模拟地震、Amatrice地震及Visso地震实验中,通过反演震源参数、获取参数的置信区间及中误差,并与Jackknife方法、Monte Carlo方法进行对比分析.实验结果表明,通过执行本文精度评定方法能够获取比Jackknife方法更加可靠的震源参数置信区间以及更加精确的精度信息.实验验证了将Bootstrap方法用于震源参数精度评定的有效性和可靠性,为研究震源参数精度评定理论研究提供了一种新的采样思路.

  相似文献   

4.
2010年高雄地震震源参数的近远震波形联合反演   总被引:4,自引:5,他引:4       下载免费PDF全文
本文改进了传统基于近震波形数据的点源震源参数反演的Cut And Paste(CAP)方法,实现了近震Pnl波、面波和远震P波、SH波的联合反演的CAPjoint算法.对2010年3月高雄地震,分别进行单独反演以及联合反演,获得各自的震源机制解及深度,其中联合反演所得的最佳双力偶机制解参数为,节面1:走向317°,倾角36°,滑移角52°,节面2:走向181°,倾角62°,滑移角114°,深度为21 km.并对不同震中距波形对本次地震以及几种典型机制解断层几何参数的敏感性进行测试.为验证联合反演方法的可靠性,本文采用重抽样思想发展而来的Bootstrap方法,对近震数据的子集及其与远震数据的联合反演所得的参数进行统计,验证了在稀疏近台条件下联合反演中添加远震数据对地震震源参数约束的作用.  相似文献   

5.
本文改进了传统基于近震波形数据的点源震源参数反演的Cut And Paste(CAP)方法,实现了近震Pnl波、面波和远震P波、SH波的联合反演的CAPjoint算法.对2010年3月高雄地震,分别进行单独反演以及联合反演,获得各自的震源机制解及深度,其中联合反演所得的最佳双力偶机制解参数为,节面1:走向317°,倾角36°,滑移角52°,节面2:走向181°,倾角62°,滑移角114°,深度为21 km.并对不同震中距波形对本次地震以及几种典型机制解断层几何参数的敏感性进行测试.为验证联合反演方法的可靠性,本文采用重抽样思想发展而来的Bootstrap方法,对近震数据的子集及其与远震数据的联合反演所得的参数进行统计,验证了在稀疏近台条件下联合反演中添加远震数据对地震震源参数约束的作用.  相似文献   

6.
用地震波波形拟合方法研究中小地震的震源机制   总被引:2,自引:0,他引:2  
采用网格搜索方法,利用区域地震三分量全波形资料,对中小地震震源机制的反演方法进行了研究。反演采用初动约束、分段求相关系数,以及幅值加权的混合方法计算误差函数,通过格点搜索法确定震源深度。理论模型试验表明,反演结果具有较好的稳定性和可靠性。云南地区区域数字地震台网的观测资料反演表明,该方法可以合理地拟合观测波形资料。将2001年7月9日(24.93°N,101.47°E)发生的MS 5.3地震的震源机制解与哈佛大学CMT进行了对比,结果表明,Z分量波形拟合略比CMT拟合的好,R分量和T分量的波形拟合明显好于CMT的对应分量。 我们对2001年4月5日在(25.77°N,102.30°E)发生的mb 4.8地震进行了震源机制反演。结果表明,该地震主要表现为走滑断层,震源机制解为(180°,80°,11°)。根据这两个地震震源机制推断的主应力轴,与该地区现代构造应力场特征基本一致。  相似文献   

7.
川滇地区速度结构的区域地震波形反演研究   总被引:22,自引:6,他引:22       下载免费PDF全文
利用云南数字地震台网的区域地震波形资料,对川滇地区的地壳上地幔速度结构进行了初步研究. 结果表明,川滇地区上地幔顶部P波速度较小,约7.8 km/s,P波速度在上地幔表现为较小的正速度梯度,S波在100~160 km深度范围内表现为弱低速层. 对于较短的观测路径,不同路径的平均P波和S波速度存在明显的横向变化. 与川滇菱形块体内部的速度结构不同,在块体边界附近可以观测到比较明显的上地壳低速层,我们认为它可能与块体边界的断裂带有关;川滇菱形块体内部存在的下地壳低速层,有利于块体向南滑动,而中上地壳没有明显低速结构,可能表明川滇菱形块体向南滑动的解耦深度至少在下地壳. 根据不同路径的反演结果,给出了云南中部地区地壳内部的平均速度结构.  相似文献   

8.
A new source location method using wave-equation based traveltime inversion is proposed to locate microseismic events accurately. With a sourceindependent strategy, microseismic events can be located independently regardless of the accuracy of the source signature and the origin time. The traveltime-residuals-based misfit function has robust performance when the velocity model is inaccurate. The new Fréchet derivatives of the misfit function with respect to source location are derived directly based on the acoustic wave equation, accounting for the influence of geometrical perturbation and spatial velocity variation. Unlike the mostly used traveltime inversion methods, no traveltime picking or ray tracing is needed.Additionally, the improved scattering-integral method is applied to reduce the computational cost. Numerical tests show the validity of the proposed method.  相似文献   

9.
We analyzed the broadband body waves of the 1992 Nicaragua earthquake to determine the nature of rupture. The rupture propagation was represented by the distribution of point sources with moment-rate functions at 9 grid points with uniform spacing of 20 km along the fault strike. The moment-rate functions were then parameterized, and the parameters were determined with the least squares method with some constraints. The centroid times of the individual moment-rate functions indicate slow and smooth rupture propagation at a velocity of 1.5 km/s toward NW and 1.0 km/s toward SE. Including a small initial break which precedes the main rupture by about 10 s, we obtained a total source duration of 110 s. The total seismic moment isM o =3.4×1020 Nm, which is consistent with the value determined from long-period surface waves,M o =3.7×1020 Nm. The average rise time of dislocation is determined to be 10 s. The major moment release occurred along a fault length of 160 km. With the assumption of a fault widthW=50 km, we obtained the dislocationD=1.3 m. From andD the dislocation velocity isD=D/0.1 m/s, significantly smaller than the typical value for ordinary earthquakes. The stress drop =1.1 MPa is also less than the typical value for subduction zone earthquakes by a factor of 2–3. On the other hand, the apparent stress defined by 2E s /M o , where andE s are respectively the rigidity and the seismic wave energy, is 0.037 MPa, more than an order of magnitude smaller than . The Nicaragua tsunami earthquake is characterized by the following three properties: 1) slow rupture propagation; 2) smooth rupture; 3) slow dislocation motion.  相似文献   

10.
《Geofísica Internacional》2014,53(2):211-220
We apply a single-step, finite-fault analysis procedure to derive a coseismic slip model for the large MW 7.4 Ometepec-Pinotepa Nacional, Mexico earthquake of 20 March 2012, using teleseismic P waveforms recorded by the Global Seismographic Network. The inversion is conducted in near-realtime using source parameters available from the USGS/NEIC and the Global Centroid Moment Tensor (gCMT) project. The fault orientation and slip angle are obtained from the gCMT mechanism assuming that the fault coincides with the shallow-dipping nodal plane. The fault dimensions and maximum rise time are based on the magnitude reported for the event. Teleseismic data from the USGS/NEIC Continuous Waveform Buffer database are used in the inversion with record start times set to the P-wave arrivals used to compute the earthquake hypocenter. The inversion is stabilized by requiring a smooth transition of slip across the fault while minimizing the seismic moment. These constraints are applied using a smoothing weight that is estimated from the inverse problem, allowing the recovery of the least-complicated rupture history in a single step. Inversion of the deconvolved, ground-displacement waveforms reveals a simple, circular rupture similar in extent to the source identified by the USGS/NEIC using body-and surface-wave data, indicating that the teleseismic P waves can provide a first-order source model for the event in near-realtime. Additional inversions conducted using velocity records identify a more-detailed rupture model characterized by an elliptical 2500 km2 source region extending updip and downdip from the hypocenter. This elliptical source preserves the orientation and overall dimensions of a dual-source slip model obtained recently by other investigators using local strong motions and global seismic waveforms. The results indicate that velocity waveforms could provide additional details of the earthquake rupture in near-realtime, finite-fault inversions using teleseismic P waves.  相似文献   

11.
井地有限线源三维电阻率反演研究   总被引:2,自引:1,他引:2       下载免费PDF全文
有限线源的三维电阻率反演,是地学探测研究的重要课题之一.本文主要讨论了以下几个方面的内容:(1)结合反演方程式改进的共轭梯度算法;(2)基于扰动法的线源Jacobi矩阵的近似计算方法及光滑系数矩阵在各个方向上的光滑因子的计算公式;(3)通过分析阻尼系数与修改量校正量对反演结果的影响,提出了利用校正系数对模型修正量进行校正;(4)较为系统的讨论了阻尼系数λ对反演分辨能力的影响,较好的改善了三维电阻率反演中的电性异常体重心的"上漂"现象.数值实验表明,改进后的共轭梯度法反演对初始模型的依赖程度较小,反演能够稳定收敛,对模型的空间位置分辨率较高.  相似文献   

12.
The new inversion algorithm developed based on the recent progress in the nonlinear programming study by us is used to invert the earthquake source process of Chi Chi earthquake M w7.6, 20 Semptember, 1999, Taiwan. A curve fault model is constructed in our inversion to make the fault model close to the real rupturing fault to reduce the influence from the discrepancy between the constructed fault model and the real rupturing fault. The results show that (1) the rupture process of the Chi Chi earthquake source lasted about 32 seconds and the main faulting occurred between 6th to 21st second after the start of the ruptures and the high slip area were mainly located at the northern segment of the fault. (2) The slip was dominated by thrust faulting. The average rake angle was 64.5°, which was very consistent with those inverted by USGS, Harvard and CWB (Central Weather Bureau of Taiwan). The amount of the moment inverted in this paper was 7.76×1020 NM, which was a slightly bigger than those inverted by USGS and Harvard. (3) A clear nucleation step existed in the source faulting process and it lasted about 6 seconds. The moment release rate accelerated obviously at the end of the nucleation step. The faulting started from the southern segment and mainly occurred at the northern segment after 10 seconds. At the end of this paper, we analyzed the reliability of the inversion result via comparing with the GPS observations and discussed its scientific signification.  相似文献   

13.
川滇地区中小震重新定位与速度结构的联合反演研究   总被引:8,自引:0,他引:8  
通过震源与速度结构联合反演, 利用2000年4月至2006年3月云南和四川区域地震台网给出的P波初至走时资料, 确定了川滇地区的三维速度结构, 同时获得了川滇地区6642次中小地震的重新定位结果。 结果表明: ① 川滇地区地震震源平均深度随震级增大而加深的特征明显, 地震震级越大, 震源深度越深, 但震源下界不超过25 km; ② 在瑞丽-龙陵、 丽江-小金河以及龙门山等断裂带以西地区, 震源深度偏浅, 大多在15 km以上, 15 km深度以下地震稀少; ③ 川滇地区中小地震分布具有与强震相同的地壳深部介质背景, 震源大多分布于正、 负异常过渡区的速度相对较高一侧, 而其下方主要为低速异常分布。  相似文献   

14.
本研究利用InSAR技术与ALOS PALSAR雷达数据,获取了2008年3月20日于田MS7.1地震视线向同震形变场,并基于该数据集和限制性最小二乘算法反演了此次地震的断层滑动分布;通过构造四大类反演方案,详细分析了InSAR观测系统中的入射角与方位角对反演结果的影响.结果表明:入射角随点位变化对反演结果有较大影响,使用其平均值将对破裂细节产生一定影响;而方位角对反演结果的影响不大,使用其平均值是一种较为理想的选择;引入入射角与方位角变化后,反演获得了较佳的于田地震同震滑动,主要集中分布于0~14 km深度附近,最大滑动量达3.2 m,矩张量为3.3×1019 N·m,相当于矩震级MW7.0.  相似文献   

15.
刘超  张勇  许力生  陈运泰 《地震学报》2008,30(4):329-339
本文利用全球地震台网(GSN)的宽频带与长周期地震波形资料,采用矩张量反演的新方法,反演得到了2008年5月12日四川省汶川县MS8.0地震及其7个较大余震(MS5.0-6.0)的矩张量解与震源时间函数等震源参数.文章首先简要叙述矩张量反演新方法的理论背景和技术途径,并以汶川大地震的一个余震为例阐述了具体的实现过程;然后给出包括主震在内的8次地震的矩张量解和震源时间函数;最后分析探讨这些结果的构造意义.本文提出的矩张量反演新方法,不但与全球矩心矩张量(GCMT)一样,可以给出点源矩张量解,而且还可以给出点源的震源时间函数. 反演得到的汶川大地震的7个较大余震的震源时间函数表明,即使是中等强度的地震也可能有复杂的震源过程;汶川大地震的多数余震发生在以逆冲为主、兼具小量走滑分量的龙门山断裂带的主断裂上,但很可能有些余震则发生在主断裂附近的次级走滑断裂上.   相似文献   

16.
全波形反演是一种高精度的地震成像方法,可以对地下介质物性参数模型进行准确的重构。然而在实际应用中,尤其是在三维复杂介质反演中,计算成本太大是该方法的一个重要缺陷。将混叠震源技术引入到频率域全波形反演中可以大幅度地降低计算成本,提高反演效率。但是使用震源编码技术也带来了两个问题:一方面,参与编码的各个震源之间会产生"串扰噪声",导致反演结果中出现假象;另一方面,基于震源编码的频率域全波形反演方法周围噪声较为敏感,使该方法对含噪数据反演质量较差。本文引入一种频率组编码方法来压制"串扰噪声",并基于震源编码技术提出一种频率域自适应全波形反演方法,通过一个与频率相关的自适应选择机制,将常规频率域全波形反演方法和基于震源编码的全波形反演方法联合起来,在保证反演质量的同时也最大程度地提高了反演效率。  相似文献   

17.
A MS8.0 earthquake occurred in Wenchuan County, Sichuan Province, China, on May 12, 2008, and subsequently, numerous aftershocks followed. We obtained the moment tensor solutions and source time functions (STFs) for the Wenchuan earthquake and its seven larger aftershocks (MS5.0~6.0) by a new technique of moment tensor inversion using the broadband and long-period seismic waveform data from the Global Seismic Network (GSN). Firstly, the theoretical background and technical flow of the new technique was brie...  相似文献   

18.
In seismological study, most of the earthquake source rupture processes are inverted via matching the waveforms. The hypocenter location and fault parameters (such as strike direction, dip angle, etc.) are assumed firstly, and the fault is divided into a …  相似文献   

19.
Abstract

The uncertainty associated with a rainfall–runoff and non-point source loading (NPS) model can be attributed to both the parameterization and model structure. An interesting implication of the areal nature of NPS models is the direct relationship between model structure (i.e. sub-watershed size) and sample size for the parameterization of spatial data. The approach of this research is to find structural limitations in scale for the use of the conceptual NPS model, then examine the scales at which suitable stochastic depictions of key parameter sets can be generated. The overlapping regions are optimal (and possibly the only suitable regions) for conducting meaningful stochastic analysis with a given NPS model. Previous work has sought to find optimal scales for deterministic analysis (where, in fact, calibration can be adjusted to compensate for sub-optimal scale selection); however, analysis of stochastic suitability and uncertainty associated with both the conceptual model and the parameter set, as presented here, is novel; as is the strategy of delineating a watershed based on the uncertainty distribution. The results of this paper demonstrate a narrow range of acceptable model structure for stochastic analysis in the chosen NPS model. In the case examined, the uncertainties associated with parameterization and parameter sensitivity are shown to be outweighed in significance by those resulting from structural and conceptual decisions.

Citation Parker, G. T. Rennie, C. D. & Droste, R. L. (2011) Model structure and uncertainty for stochastic non-point source modelling applications. Hydrol. Sci. J. 56(5), 870–882.  相似文献   

20.
A method for simultaneous determination of mixed model parameters,which have different physical dimensions or different responses to data,is presented.Mixed parameter estimation from observed data within a single model space shows instabilities and trade-offs of the solutions. We separate the model space into N-subspaces based on their physical properties or computational convenience and solve the N-subspaces systems by damped least-squares and singular-value decomposition. Since the condition number of each subsystem is smaller than that of the single global system,the approach can greatly increase the stability of the inversion. We also introduce different damping factors into the subsystems to reduce the tradeoffs between the different parameters. The damping factors depend on the conditioning of the subsystems and may be adequately chosen in a range from 0.1 % to 10 % of the largest singular value. We illustrate the method with an example of simultaneous determination of source history,source geometry,and hypocentral location from regional seismograms,although it is applicable to any geophysical inversion.  相似文献   

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

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