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.     

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%.     

Spectral and Multiscale Signal-to-Noise Thresholding of Spherical Vector Fields

The basic concepts of spectral and multiscale selective reconstruction of (geophysically relevant) vector fields on the sphere from error-affected data is outlined in detail. The reconstruction mechanism is formulated under the assumption that spectral as well as multiscale approximation is well-representable in terms of only a certain number of expansion coefficients at the various resolution levels. It is shown that spectral denoising by means of orthogonal expansions in terms of vector spherical harmonics reflects global a priori information of the noise (e.g., in form of a covariance tensor field), whereas multiscale signal-to-noise thresholding can be performed under locally dependent noise information within a multiresolution analysis in terms of spherical vector wavelets. An application of the multiscale formalism to Earth's magnetic field determination is presented.     

我国月降水和气温网格点资料的处理和分析

用逐步搜索逼近法的客观分析技术计算出１９５１－１９９２年中国月降水量和月平均气温的格点资料，并利用图象识别原理对格点值和站点观测值进行了比较。     

凸度量空间中渐近拟非扩张映象新的带误差的Ishikawa迭代逼近

在凸度量空间中，利用新的解析法，征明了由集合序列{On}产生的在多值渐近拟非扩张映象下的新的带误差的lshikawa迭代序列的收敛性。     

Reduced Hamiltonian formalism for general-relativistic adiabatic fluids and applications

For the dynamics of general-relativistic adiabatic fluids a reduced Hamiltonian formalism is presented. For slowly moving fluids an approximation scheme is defined. The first steps of the scheme are worked out explicitly. Radiation damping and emission, and backscattering are discussed. The usefulness of the formalism for numerical relativity is pointed out.     

微分算子展开与耦合势三维有限体积法快速计算随钻超深前视电阻率测井响应非线性Born逼近

本文将非均质各向异性地层中电磁场耦合势Helmhotz方程分解为各向异性背景场方程与散射场方程,并应用算子展开技术推导出耦合势Born级数解,建立一套三维非均质各向异性地层中随钻超深前视电阻率测井响应各阶散射电磁场与非线性Born逼近的有效计算方法.首先,利用传输线法与二维插值技术确定背景电磁场空间分布,得到零阶Born逼近解,然后根据电磁场耦合势的Born级数解推导出各阶散射电磁场方程与递推关系,通过三维有限体积法实现各阶散射场耦合势方程的离散,在此基础上结合MKL PARDISO并行技术,通过递推方式逐步计算各阶散射电磁场,并根据Born级数部分和确定不同阶次的非线性Born逼近.最后利用数值模拟结果对正演算法加以检验,并深入研究分析各阶散射场的变化特征与Born级数的收敛性质.

     

数字高程模型误差及其评价的问题综述

从数字高程模型(DEM)传递误差、基于中误差的DEM误差模型及其主要问题、DEM误差分布实验和DEM内插误差新认识几个方面分析了当前DEM误差研究的主要进展,用"中误差"讨论DEM传递误差是建立在测量误差传递理论基础之上的,但沿用"中误差"来讨论DEM内插模型逼近误差和DEM整体误差却缺乏理论依据。DEM误差分布的空间相关性实验对DEM中误差评价法所应具备的随机误差性提出质疑,却可以用基于逼近理论的DEM内插模型来解释,说明用逼近误差理论研究DEM内插误差的途径是正确、可行的。     

Hierarchical decomposition of variance with applications in environmental mapping based on satellite images

A quadtree-based image segmentation procedure (HQ) is presented to map complex environmental conditions. It applies a hierarchical nested analysis of variance within the framework of multiresolution wavelet approximation. The procedure leads to an optimal solution for determining mapping units based on spatial variability with constraints on the arrangement and shape of the units. Linkages to geostatisiics are pointed out, but the HQ decomposition algorithm does not require any homogeneity criteria. The computer implementation can be parameterized by either the number of required mapping units or the maximum within-unit variance, or it can provide a spectrum of significances of nested ANOVA. The detailed mathematical background and methodology is illustrated by a salt-affected grassland mapping study (Hortobágy, Hungary), where heterogeneous environmental characteristics have been sampled and predicted based on remotely sensed images using these principles.