计算视复电阻率的新的近似公式

本文推导并评价了一个新的计算视复电阻率的近似公式,其视振幅表达式为|ρ_s（ω）|≈ρ_s（ω）|_{ρs（ω）=ρs（ω）|},视相位表达式为φ_s（ω）≈sum from i=1 to n（B_i（ω）φ_i（ω））,这里B_i（ω）是按各介质复电阻率振幅计算的稀释系数。此近似公式比迄今为止已发表过的同类近似公式都精确,可在频率域激电法数值模拟计算视复电阻率时应用。     

Neural network-based GPS GDOP approximation and classification

In this paper, the neural network (NN)-based navigation satellite subset selection is presented. The approach is based on approximation or classification of the satellite geometry dilution of precision (GDOP) factors utilizing the NN approach. Without matrix inversion required, the NN-based approach is capable of evaluating all subsets of satellites and hence reduces the computational burden. This would enable the use of a high-integrity navigation solution without the delay required for many matrix inversions. For overcoming the problem of slow learning in the BPNN, three other NNs that feature very fast learning speed, including the optimal interpolative (OI) Net, probabilistic neural network (PNN) and general regression neural network (GRNN), are employed. The network performance and computational expense on NN-based GDOP approximation and classification are explored. All the networks are able to provide sufficiently good accuracy, given enough time (for BPNN) or enough training data (for the other three networks).     

Splitting phenomenon of a higher-order shallow water theory associated with a longitudinal planetary waves

The Cauchy–Poisson free boundary problem associated with a nonstationary motion of a perfect incompressible fluid circulating around the equatorial plane of a planet is considered. It is shown that the corresponding theory of a higher-order shallow approximation admits two functionally independent systems, while the classical problem for the flat bottom admits only one system.     

基于小波去噪的地铁沉降监测分析

地铁等大型工程施工过程引起的地表及其周边建筑物的变形及变形趋势分析,对于确保城市及工程安全极为重要。研究有关沉降监测方法及可靠性分析,估计地表及其周边建筑物可能产生的变形,对周边建筑物的安全,以及工程本身顺利实施有着重要的现实意义。本文将沉降观测数据序列作为带有白噪声数字信号,采用小波多孔算法进行了监测数据的去噪处理,取得了较好的效果。同时论文提出采用分段拟合的方法进行数据拟合,用对数函数曲线来描述土体沉降较快的第二阶段变形趋势,变形预报较为准确。     

Viscoelastic representation of surface waves in patchy saturated poroelastic media

Wave-induced flow is observed as the dominated factor for P wave propagation at seismic frequencies. This mechanism has a mesoscopic scale nature. The inhomogeneous unsaturated patches are regarded larger than the pore size, but smaller than the wavelength. Surface wave, e.g., Rayleigh wave, which propagates along the free surface, generated by the interfering of body waves is also affected by the mesoscopic loss mechanisms. Recent studies have reported that the effect of the wave-induced flow in wave propagation shows a relaxation behavior. Viscoelastic equivalent relaxation function associated with the wave mode can describe the kinetic nature of the attenuation. In this paper, the equivalent viscoelastic relaxation functions are extended to take into account the free surface for the Rayleigh surface wave propagation in patchy saturated poroelastic media. Numerical results for the frequency-dependent velocity and attenuation and the time-dependent dynamical responses for the equivalent Rayleigh surface wave propagation along an interface between vacuum and patchy saturated porous media are reported in the low-frequency range (0.1–1,000 Hz). The results show that the dispersion and attenuation and kinetic characteristics of the mesoscopic loss effect for the surface wave can be effectively represented in the equivalent viscoelastic media. The simulation of surface wave propagation within mesoscopic patches requires solving Biot’s differential equations in very small grid spaces, involving the conversion of the fast P wave energy diffusion into the Biot slow wave. This procedure requires a very large amount of computer consumption. An efficient equivalent approach for this patchy saturated poroelastic media shows a more convenient way to solve the single phase viscoelastic differential equations.     

Hybrid three-dimensional finite-difference and finite-element analysis of seismic wave induced fluid–structure interaction of a vertical cylinder

The two-dimensional finite-difference scheme has been extended to three dimensions to solve nonlinear hydrodynamic pressures and structural responses of a deformable, vertical and circular surface-piercing offshore cylinder during earthquakes. A complete three-dimensional analysis has been made with both the three-dimensional equations of motion and the simultaneous action of three components of ground acceleration included in the analysis. Not only the magnitude but also the direction of the acting ground motion can be varied with time. The dynamic response of a cylinder is approximated by the displacements in the fundamental modes of vibration. A comparison of the dynamic displacement of the cylinder with and without surrounding sea water has been made. The flexibility of the offshore cylinder can significantly increase the hydrodynamic pressures acting on cylinder faces, that is, the fluid-structure interaction is necessary in offshore cylinder analysis. Although the hydrodynamic pressure induced by the vertical ground acceleration of the El Centro 1979 earthquake is significant, the calculated structural dynamic response of a cylinder is very small and the corresponding resultant hydrodynamic force is almost nil. The hydrodynamic force induced by two-horizontal ground acceleration is about the same as that by three simultaneous components of ground acceleration. For a solid and stubbier circular cylinder, the vertical component of ground acceleration may be neglected.     

A STUDY ON THE ALGORITHM FOR EXTRACTING EARTHQUAKE THERMAL INFRARED ANOMALIES BASED ON THE YEARLY TREND OF LST

There are thermal infrared anomalies(TIA)before earthquake, and TIA has become one of the important parameters for assessing regional earthquake risk. However, not all of the surface infrared anomalies are related to tectonic activities or earthquakes. How to eliminate the influence of non-structural factors and extract the weak signals from strong disturbances is the key and difficult point for tectonic activities studies based on the thermal infrared remote sensing techniques. Land surface temperature(LST)background field is the basis for thermal infrared anomalies extraction. However, the established background fields in previous researches cannot eliminate the influence of climate changes, so the accuracy of thermal anomaly extraction is limited. Now an improved method is proposed. Combined with the periodic character of LST time series, harmonic analysis is lead into the process of LST background field establishment. Specifically, the yearly trend of LST is fitted based on Fourier Approximation method. As a new background field, the yearly trend is dynamic, includes the local and the yearly information. Then, based on the rule of "kσ", the earthquake anomalies, calculated by RST with the yearly trend of LST, can be extracted. At last, the effectiveness of the algorithm can be tested by the quantitative analysis of anomalies with anomaly area statistics, anomaly intensity statistics and distance index statistics. The Wenchuan earthquake was discussed again based on the proposed algorithm with MODIS land temperature products in 2008. The results show that, there were obvious pre-earthquake thermal anomalies along the Longmen Mountains faults with a longer time; but there were no anomalies when the earthquake happened; and the post-earthquake thermal anomalies occurred with much smaller amplitudes and scopes. Compared with the results derived from the traditional RST which is based on the spatial average of LST values, the TIA extracted by the new RST, which is based on the yearly trend of LST, is more fit with the active faults, and the process of the anomalies occurring and removing can be described in more detail. Therefore, as the background field to extract earthquake anomalies, the yearly trend of LST is more reliable.     

基于FFT的精确高程异常估计问题

本文讨论了FFT估计高程异常及Romberg算法解求截断系数的有关问题;并利用雪林-拉普阶方差估计FFT法的误差值,指出ψ_0=1°.0的球帽是最佳选择。     

Constructive approximation and numerical methods in geodetic research today – an attempt at a categorization based on an uncertainty principle

Current activities and recent progress on constructive approximation and numerical analysis in physical geodesy are reported upon. Two major topics of interest are focused upon, namely trial systems for purposes of global and local approximation and methods for adequate geodetic application. A fundamental tool is an uncertainty principle, which gives appropriate bounds for the quantification of space and momentum localization of trial functions. The essential outcome is a better understanding of constructive approximation in terms of radial basis functions such as splines and wavelets. Received: 4 May 1999 / Accepted: 21 May 1999     

Probability distribution of random wave forces in weakly nonlinear random waves

Based on the second-order random wave theory, the joint statistical distribution of the horizontal velocity and acceleration is derived using the characteristic function expansion method. From the joint distribution and the Morison equation, the theoretical distributions of drag forces, inertia forces and total random wave forces are determined. The distribution of inertia forces is Gaussian as that derived using the linear wave model, whereas the distributions of drag forces and total random forces deviate slightly from those derived utilizing the linear wave model. It is found that the distribution of wave forces depends solely on the frequency spectrum of sea waves associated with the first order approximation and the second order wave–wave interaction.