Probing mid-mantle heterogeneity using PKP coda waves

We investigate the utility of PKP coda waves for studying weak scattering from small-scale heterogeneity in the mid-mantle. Coda waves are potentially a useful probe of heterogeneity in the mid-mantle because they are not preferentially scattered near the CMB, as PKP precursors are, but are sensitive to scattering at all depths. PKP coda waves have not been used for this purpose historically because of interference with other late-arriving energy due to near-surface resonance and scattering. Any study of deep mantle scattering using coda waves requires the removal of near-surface effects from the data. We have analyzed 3624 recordings of PKP precursors and coda made by stations in the Incorporated Research Institutions for Seismology (IRIS) Global Seismographic Network (GSN). To study the range and time dependence of the scattered waves, we binned and stacked envelopes of the recordings. We have considered precursors that arrive within a 20 s window before PKP and coda waves in a 60 s window after PKP. The PKP scattered waves increase in amplitude rapidly with range as predicted by scattering theory. At ranges below 125°, we predict and observe essentially no scattered energy preceding PKP. Coda amplitudes at these ranges are independent of range and provide an estimate of energy due to near-surface effects that we can expect at all ranges. We use the average coda amplitude at ranges from 120 to 125° to correct coda amplitudes at other ranges. PKP coda waves show a strong dependence on time and range and are clearly influenced by scattering in the lower mantle. PKP coda waves, however, do not provide a tighter constraint on the vertical distribution of mantle heterogeneity than is provided by precursors. This is due, in part, to relatively large scatter in coda amplitudes as revealed by a resampling analysis. Modeling using Rayleigh–Born scattering theory and an exponential autocorrelation function shows that PKP coda amplitudes are not highly sensitive to the vertical distribution of heterogeneity in the mantle. To illustrate this we consider single-scattering in two extreme models of mantle heterogeneity. One allows heterogeneity just at the CMB; the other includes heterogeneity throughout the mantle. The amplitudes of precursors are tightly constrained by our stack and support our earlier conclusion that small-scale heterogeneity is uniformly distributed throughout the lower mantle. The best-fit model includes 8 km scale length heterogeneity with an rms velocity contrast throughout the mantle of 1%.     

非线性模型线性近似的容许曲率

由于不同的非线性模型具有不同的非线性强度,使得一些非线性模型可以线性近似,而另一些则不能。本文介绍度量非线性强度的方法,提出判断非线性模型能否线性近似的数值标准——容许曲率。     

坐标相似变换模型的探讨

常用的坐标转换模型,讨论了大地高的不确定性对应用Bursa模型进行坐标转换的影响。结合工程应用的实际情况,选择平面坐标转换模型对北京地方坐标与西安80坐标进行转换,转换精度满足实际需求。     

A new approach for predicting the extreme response of structures

In this article a method to estimate the mean upcrossing intensity, μ(u), of a stochastic process is proposed. It is assumed that the stochastic process is a sum of a Gaussian process and a second-order correction term. The method is based on the two-dimensional Saddlepoint approximation. The accuracy of the method is tested on processes having analytical solutions for μ(u). Numerical examples are given where the stochastic process represents (i) the horizontal response of a floating offshore structure in a Gaussian sea, and (ii) the response of a structure subjected to a Gaussian wind velocity process. In addition, the estimates are compared to empirical upcrossing intensities of simulated responses. For case (ii), the obtained μ(u) estimates are compared to estimates obtained by numerical integration.     

起伏地表条件下的声波散射数值模拟的积分方程法

从散射理论的角度来看,起伏地表可以看作是一种特殊的扰动介质,因此应用散射积分方程求解起伏地表条件下的散射场在理论上是可行的。从三维频率域声波方程出发,由格林函数定理,得出起伏地表条件下的散射积分方程。散射积分方程为关于起伏地表的面积分和与速度扰动体有关的体积分之和,同时给出了格林函数在奇异点的积分方法。由于数值离散求解积分方程存在着计算时间太长和存储内存不足的问题,采用电磁散射积分方程的拟解析近似的方法。在假设反射函数为缓变函数的基础上,最终得到其近似表达式,因此散射场的数值求解不必再借助于代数方程组,只要进行数值积分即可。这种方法避开了传统数值计算方法存在的问题,为地震散射波场快速正演模拟打下了基础。理论分析表明,这种方法适用于小扰动的问题。当扰动较大时,拟解析近似会产生较大的误差。     

二元函数逼近在斜坡可靠度分析中的应用研究

极限平衡分析方法是斜坡稳定性评价中的常用方法,在长期的工程实践中积累了丰富的经验,但其不能考虑斜坡岩土体中实际存在的不确定性,在应用中具有一定的局限性。可靠度分析方法可有效地考虑斜坡系统内的不确定性和相关性,但因状态函数偏导数的求解比较困难,使可靠度分析方法在实际中应用不便。为解决上述问题,根据二元函数插值逼近原理,在矩形区域上构造拉格朗日不完全双二次多项式逼近状态函数,从而近似地计算状态函数的偏导数,求得状态函数的均值和方差,并利用精度较高的一次二阶矩方法来计算斜坡的可靠指标和破坏概率。据鄂西恩施地区马堡营滑坡实例分析表明,引入二元函数插值逼近的一次二阶矩方法计算结果与剩余推力法及Monte-carlo模拟方法结果一致,其精度可满足工程需求。     

Time-domain identification of dynamic properties of layered soil by using extended Kalman filter and recorded seismic data

A novel time-domain identification technique is developed for the seismic response analysis of soil-structure interaction. A two-degree-of-freedom (2DOF) model with eight lumped parameters is adopted to model the frequency-dependent behavior of soils. For layered soil, the equivalent eight parameters of the 2DOF model arc identified by the extended Kalman filter (EKF) method using recorded seismic data. The polynomial approximations for derivation of state estimators are applied in the EKF procedure. A realistic identification example is given for the layered-soil of a building site in Anchorage, Alaska in the United States. Results of the example demonstrate the feasibility and practicality of the proposed identification technique. The 2DOF soil model and the identification technique can be used for nonlinear response analysis of soil-structure interaction in the time-domain for layered of complex soil conditions. The identified parameters can be stored in a database for use in other similar soil conditions. If a universal database that covers information related to most soil conditions is developed in the future, engineers could conveniently perform time history analyses of soil-structural interaction.     

Including Geological Uncertainty and Economic Analysis in a Rapid Simulation of Hydrocarbon Exploration

Manly's approximation method has been applied to hydrocarbon discovery process modeling in order to approximate the expected value and the standard deviation of the total discovery volume for a given exploratory effort. A major benefit of this method is that it allows the model to run much faster than the regular simulation method, yet it gives accurate results. This paper extends the benefit of Manly's approach by allowing the approximation method to incorporate the uncertainties in geological parameters that drive Manly's approximation, in order to provide an approximation of the complete distribution of total discovery volume that can result from exploration activity. In addition, it allows the model to include economic parameters into an evaluation of the economic worth of the results of exploration activity, producing distributions of net present value within a short period of time. The offshore Nova Scotia Shelf Basin is selected for demonstrating the methodology.     

Point- and curve-based geometric conflation

Geometric conflation is the process undertaken to modify the coordinates of features in dataset A in order to match corresponding ones in dataset B. The overwhelming majority of the literature considers the use of points as features to define the transformation. In this article we present a procedure to consider one-dimensional curves also, which are commonly available as Global Navigation Satellite System (GNSS) tracks, routes, coastlines, and so on, in order to define the estimate of the displacements to be applied to each object in A. The procedure involves three steps, including the partial matching of corresponding curves, the computation of some analytical expression, and the addition of a correction term in order to satisfy basic cartographic rules. A numerical example is presented.