起伏地表条件下地震波场数值模拟有限积分变换有限差分方法

针对起伏地表条件下的地震波数值模拟问题,提出将起伏的地表映射到一个规则的长方形网格坐标系中,并在此基础上推导出变换域中的波动方程。根据导出来的波动方程应用有限余弦变换有限差分方法进行地震波场的数值模拟,以解决起伏地表条件下的数值模拟。     

Modeling errors in upward continuation for INS gravity compensation

Accurate upward continuation of gravity anomalies supports future precision, free-inertial navigation systems, since the latter cannot by themselves sense the gravitational field and thus require appropriate gravity compensation. This compensation is in the form of horizontal gravity components. An analysis of the model errors in upward continuation using derivatives of the standard Pizzetti integral solution (spherical approximation) shows that discretization of the data and truncation of the integral are the major sources of error in the predicted horizontal components of the gravity disturbance. The irregular shape of the data boundary, even the relatively rough topography of a simulated mountainous region, has only secondary effect, except when the data resolution is very high (small discretization error). Other errors due to spherical approximation are even less important. The analysis excluded all measurement errors in the gravity anomaly data in order to quantify just the model errors. Based on a consistent gravity field/topographic surface simulation, upward continuation errors in the derivatives of the Pizzetti integral to mean altitudes of about 3,000 and 1,500 m above the mean surface ranged from less than 1 mGal (standard deviation) to less than 2 mGal (standard deviation), respectively, in the case of 2 arcmin data resolution. Least-squares collocation performs better than this, but may require significantly greater computational resources.     

河南省地热资源开发利用现状及展望

河南省地热资源丰富,在开封、周口、舞阳、济源、南阳盆地以及濮阳凹陷、汤阴地堑有望找到很好的热储层。现已开发的热储层为埋深８００～１３００ｍ的第三系中细砂层,水温４０～５０℃,水质优良。如果加深钻井深度至２０００ｍ,则有望找到７０℃以上的高温地热水,届时省内地热资源开发利用将上一个新台阶,达到供热采暖水平。文中对河南省地热资源的开发利用现状作一简要介绍,并对如何提高开发利用水平和资源保护提出建议。     

Closed Form Solutions of the Two-Dimensional Turning Bands Equation

The turning bands method generates realizations of isotropic Gaussian random fields by means of appropriately summed line processes. For two-dimensional simulations the relation between the isotropic correlation function of the random field and the correlation function to be simulated along the lines is given by an integral equation of Abel type. We present closed form solutions of this integral equation for almost all two-dimensional correlation models encountered in practice and discuss their numerical implementation. As an additional benefit, our tables and illustrations serve as a concise guide to correlation models useful in geostatistics.     

Electromagnetic 2.5-dimensional forward modelling with a boundary integral formulation

An effective and accurate technique for the numerical solution of 2-D electromagnetic scattering problems with 3-D sources is presented. This solution introduces a set of the usual boundary integral equations and uses a scalar Green's function. In this scalar version, the unknowns of the problem are the boundary values of the longitudinal fields and their normal derivatives in the Fourier domain. A generalization of the usual boundary integral formulation enables us to handle a large class of models composed of piecewise homogeneous domains, including contiguous domains, multiply-connected domains and unbounded domains. This formulation involves the solution of a system of linear equations, and results in a significant saving in computation time in comparison with other rigorous methods.
  The requirements for the numerical implementation of this solution are described in detail. Numerical tests were carried out using the important example of electromagnetic tomography. The specific symmetry properties of the response function in this case are illustrated. Numerical accuracy is verified over a large frequency range, up to 1  MHz.     

基于三轴加速度的波浪测量技术研究

以浮标为载体,以TMS320C6416型DSP处理器、三轴加速度传感器、电子罗盘为硬件电路核心,利用加速度频域二次积分的波高计算方法以及三轴加速度倾角补偿电子罗盘的波向测量方式,进行海洋波浪要素测量技术研究,获取了某一海域不同时段的波浪特征值。     

海底管道整体式止屈器性能研究与穿越压力预测

整体式止屈器可有效抑制海底管道屈曲传播并避免管道的大规模屈曲失效。为探究整体式止屈器的止屈特性,预测其穿越压力,利用ABAQUS软件建立有限元模型,分析了整体式止屈器厚度、长度、径厚比等参数对穿越压力和穿越形式的影响机制,揭示了整体式止屈器的止屈特性。研究结果表明:在平行穿越阶段,增加止屈器厚度、有效长度可明显提升穿越压力;而在垂直穿越阶段,穿越压力增长速率明显放缓。结合有限元模拟结果,拟合得到了整体式止屈器在不同穿越形式下的穿越压力预测公式,其相较于已有经验公式准确性更高,对整体式止屈器的设计和建造具有重要的指导意义。     

Geomorphological analysis of the morphometric characteristics that determine the volume of sediment yield of Wadi Al-Arja,South Jordan

This study addresses the morphometric variables that determine the sediment yield in Wadi Al-Arja through the analysis of the impact of different morphometric characteristics along the course of the valley on its sediment yield, as well as the analysis of spatial and formal dimensions and morphologies of the basin and its relationship to the sediment yield. The study also addresses the size of variation in the volume of sediment yield of the river tributaries that make up the water network of the valley under the differences of its morphometric and hydrological characteristics. The study found several results, most notably: The classification of Wadi Al-Arja basin according to the hypsometric integral value (72.1%) within an uneven topography, which increases the force of the erosive activity and the size of the sediment yield if the basin was exposed to moist climatic periods. The study also showed the presence of convergence in the intra-spaces between the river tributaries in the water network. This increases the volume of water flows when these tributaries meet with each other and thus increasing their erosive ability and sediment yield.The study also showed the presence of marked variation in the sediment yield of the river tributaries depending on the differences in its morphometric characteristics. The results of the step-wise regression analysis confirmed the importance of the morphometric and hydrological variables, and plant coverage in interpreting the variation in the size of the sediment yield of the river tributaries of different stream order in Wadi Al-Arja basin, where these variables interpreted 43% of the total variation, with statistical significance less than 0.05.     

GPU accelerated marine data visualization method

The study of marine data visualization is of great value. Marine data, due to its large scale, random variation and multiresolution in nature, are hard to be visualized and analyzed. Nowadays, constructing an ocean model and visualizing model results have become some of the most important research topics of ‘Digital Ocean'. In this paper, a spherical ray casting method is developed to improve the traditional ray-casting algorithm and to make efficient use of GPUs. Aiming at the ocean current data, a 3D view-dependent line integral convolution method is used, in which the spatial frequency is adapted according to the distance from a camera. The study is based on a 3D virtual reality and visualization engine, namely the VV-Ocean. Some interactive operations are also provided to highlight the interesting structures and the characteristics of volumetric data. Finally, the marine data gathered in the East China Sea are displayed and analyzed. The results show that the method meets the requirements of real-time and interactive rendering.     

A new model for the vibration isolation via pile rows consisting of infinite number of piles

In this study, on the basis of the Floquet transform method, a numerical model for the simulation of the vibration isolation via multiple periodic pile rows with infinite number of piles is established. By means of the fictitious pile method due to Muki and Sternberg, the second kind of Fredholm integral equations for the pile rows are developed by using the fundamental solutions for the half‐space and the compatibility conditions between the piles and half‐space. Employing the Floquet transform method, integral equations for the pile rows in the wavenumber domain are then derived. Solution of the integral equations yields the wavenumber domain solution for the pile rows. The space domain solution can then be retrieved by inversion of the Floquet transform. Numerical results show that the proposed model with the Floquet transform method is in a good agreement with those of the conventional direct superposition method. On the basis of the new model, influences of the spacing between neighboring piles, the Young's modulus of the piles, and the pile length on the vibration isolation effect of the pile rows are investigated. Numerical simulations conducted in this study show that compared with the direct superposition method, the efficiency of the proposed model for simulation of the vibration isolation via pile rows is very high. Copyright © 2012 John Wiley & Sons, Ltd.