常Q滞弹性介质地震波动数值模拟——时域本构优化逼近

在地震波动模拟中计入常Q滞弹性阻尼,可有效降低模拟波形的误差.就时域有限差分和有限元模拟而言,常基于广义标准线性体建立阻尼介质的时域本构逼近.广义标准线性体由若干标准线性体并联得到,增加标准线性体个数能有效提高模拟精度,但计算量及计算存储将成倍增长.目前尚未有普适的标准线性体个数优化取值方案.本文基于广义标准线性体参数的非线性优化拟合方法,详细分析了时域本构逼近误差的影响因素,清楚揭示了逼近误差仅取决于频带宽度,与频带上下限取值无关这一特性,阐明了构建具有普适性标准线性体个数优化取值方案的可行性.论证了波形模拟精度主要取决于波传播距离与模拟波长的比值以及标准线性体的个数取值.综合考虑上述两个控制因素,结合在波动正反演问题中广为采纳的波形时频误差衡量准则,对不同Q值介质给出了标准线性体个数优化取值表.进一步,本文提出采用不同个数标准线性体以近似不同Q值的阻尼介质时域本构,解决了以往波动数值模拟中统一采用相同个数标准线性体而导致的计算量及计算存储浪费或模拟精度低下等问题,并基于数值实验验证了这一方法的精度.本文工作对推进滞弹性介质波动数值模拟及其在全波形反演问题中的应用具有理论价值和实践意义.     

基于广义模型约束的时间域航空电磁反演研究

由于航空电磁具有海量数据,因此快速有效的成像和反演手段至关重要.本文针对层状介质模型推导与实现了广义模型约束条件下时间域航空电磁一维反演.从正则化反演的目标函数出发,通过改变模型约束项构造Lp范数反演和聚焦反演,进而通过改变模型求解域构造出基于小波变换的稀疏约束反演.针对不同反演方法目标函数的构建方式,本文进一步从数学原理上分析不同反演方法的预期效果,并通过理论模型和实测数据进行验证.结果表明L0.8范数反演、聚焦反演和基于小波变换的稀疏约束反演可以得到更符合地下层状介质陡变界面的反演结果.     

Research note: deblended-data reconstruction using generalized blending and deblending models

We introduce a concept of generalized blending and deblending, develop its models and accordingly establish a method of deblended-data reconstruction using these models. The generalized models can handle real situations by including random encoding into the generalized operators both in the space and time domain, and both at the source and receiver side. We consider an iterative optimization scheme using a closed-loop approach with the generalized blending and deblending models, in which the former works for the forward modelling and the latter for the inverse modelling in the closed loop. We applied our method to existing real data acquired in Abu Dhabi. The results show that our method succeeded to fully reconstruct deblended data even from the fully generalized, thus quite complicated blended data. We discuss the complexity of blending properties on the deblending performance. In addition, we discuss the applicability to time-lapse seismic monitoring as it ensures high repeatability of the surveys. Conclusively, we should acquire blended data and reconstruct deblended data without serious problems but with the benefit of blended acquisition.     

一种地方与区域地震震源机制反演技术:广义极性振幅技术(二)——对实际震例的应用

为了检验广义极性振幅技术(GPAT)的实用性,我们利用GPAT反演确定了49次实际地震的震源机制、矩震级和震源深度.为了检验GPAT对地方地震、区域地震以及远震的实用性,我们选用了震级范围约为ML0.2~MS7.0之间的地震与震中距范围约在5~8000km之间的观测资料.对反演结果的分析表明,利用GPAT获取的震源机制结果在合理误差范围内是正确的,利用GPAT获取的矩震级结果是可靠的,而利用GPAT获取的震源深度总体上似乎比常规定位深度深约0.6km.总体而言,GPAT在震源机制、矩震级和震源深度的反演方面表现出良好的实用性.     

一种地方与区域地震震源机制反演技术：广义极性振幅技术（一）——原理与数值实验

基于数十年来已有的研究进展,提出了一种地方和区域地震震源机制反演技术——广义极性振幅技术（GPAT）,并通过一系列数值实验检验了这种技术的可行性和抗干扰能力.首先,从地震波场概念出发,利用P波初动极性与广义震相振幅构建矢量,建立反演系统,并给出求解技术;然后,考虑影响反演结果的各种因素,包括台站布局、台站数目、随机噪声、震中位置误差、震源深度误差和速度模型误差,分别进行了单一因素影响测试;最后,同时考虑各种因素进行了综合测试.实验结果表明,GPAT是可行的,具有良好的抗干扰能力.需要强调的是,在众多影响因素中,速度模型误差对反演结果的影响最大.     

基于GPD分布的黑河流域极端降水频率特征分析

采用广义帕雷托分布(GPD)对黑河流域极端降水的频率特征进行了分析。采用百分位数法、Hill图法、年交叉率法选取了极端降水阈值,借助L矩法对GPD分布的参数进行了估计;采用Mann-Kendall和Pettitt方法对超阈值日极端降水序列的平稳性进行了检验。结果表明：① 采用百分位数法和Hill图法得到的日极端降水阈值差别很大,在综合考虑以上两种方法的基础上,结合年交叉率法,最终确定了流域各站点日极端降水的阈值;② 根据Mann-Kendall和Pettitt方法,超阈值日极端降水序列基本满足平稳性假定;通过Kolmogorov-Smirnov检验以及理论频率曲线和经验频率曲线的拟合程度可以看出,GPD分布能够很好地拟合研究区极端降水的分布特征;③ 通过分析理论累积频率达到90%以上的极端降水可以发现,黑河流域20世纪60年代发生极端降水的次数最多,其次是90年代以后,70年代、80年代发生极端降水次数较少;在过去51年间,下游发生极端降水的次数最多,其次是中游和上游。     

均匀弹性分层介质模型中的同震位移、应变以及应力

地震发生时产生的同震位移、应变以及应力变化,特别是同震应力的变化,在地震触发等问题的研究中有着重要的意义.本文进一步发展了基于均匀弹性水平层状介质,利用广义反射透射系数矩阵和离散波数计算同震位移的方法,使之可以计算相应的应变、应力以及同震库仑应力变化.可适用于多种情况,接收点可以位于地表以及地表以下,震源类型可以是剪切位错源以及拉张位错源.通过与半无限介质的解析解相比较,结果一致,验证了方法的可靠性以及计算精度,可以用于计算地震之后库仑应力变化,为判断余震分布提供参考.在计算同震位移时,使用了梯形积分与Filon积分相结合的积分方式,即提高了同震位移计算的速度,又保证了计算精度,有利于反演问题研究.     

Multifractal characteristic analysis of near-fault earthquake ground motions

This study aims to reveal the multi-scaling behavior and quantify the irregularity of near-fault earthquake ground motions from a new perspective of multifractal theory. Based on multifractal detrended fluctuation analysis, the multifractal characteristic parameters of acceleration time series for typical near-fault ground motions are calculated, and their correlations with two period parameters (i.e., mean period T_m and characteristic period T_c) and box-counting fractal dimensions are analyzed. Numerical results of strong nonlinear dependence of generalized Hurst exponents h(q) upon the fluctuation orders q indicate that near-fault ground motions present the multifractal properties and long-range correlation obviously. Furthermore, the scaling exponent h(2) of near-fault records has a strong correlation with their periods T_m and T_c, and strongly negative correlation with their box dimension. Moreover, h(2) can be regarded as a measure of frequency content and irregularity degree of strong earthquake ground motions. Finally, it is pointed out that the long-range correlation of small and large fluctuation is the major source of multifractality of near-fault ground motions.     

Estimation of Site Effects on Stations in the Capital Circle Region Using the Generalized Inversion Method

Generalized Inversion Method has been used to estimate the spatial variation of site effects,using the digital data of SH-waves recorded by 63 stations in the Capital Circle Region of China from 2001 to 2006.We gained the site effects of all stations participating in the calculation.We found that the site effect of rock was stabile and about 1.0 from 1.0Hz to 10.0Hz,while the site effect of deposit was high in low frequencies,about 3 ～ 7 from 1.0Hz to 8.0Hz,and the site effect was protuberant at about 5.0Hz,then fell as the frequency increased.The result shows the shape and intensity of station site effects are mainly influenced by the lithology below the station,and possibly also by the local geological structure.