Comparison of methods to model the gravitational gradients from topographic data bases

A number of methods have been developed over the last few decades to model the gravitational gradients using digital elevation data. All methods are based on second-order derivatives of the Newtonian mass integral for the gravitational potential. Foremost are algorithms that divide the topographic masses into prisms or more general polyhedra and sum the corresponding gradient contributions. Other methods are designed for computational speed and make use of the fast Fourier transform (FFT), require a regular rectangular grid of data, and yield gradients on the entire grid, but only at constant altitude. We add to these the ordinary numerical integration (in horizontal coordinates) of the gradient integrals. In total we compare two prism, two FFT and two ordinary numerical integration methods using 1" elevation data in two topographic regimes (rough and moderate terrain). Prism methods depend on the type of finite elements that are generated with the elevation data; in particular, alternative triangulations can yield significant differences in the gradients (up to tens of Eötvös). The FFT methods depend on a series development of the topographic heights, requiring terms up to 14th order in rough terrain; and, one popular method has significant bias errors (e.g. 13 Eötvös in the vertical–vertical gradient) embedded in its practical realization. The straightforward numerical integrations, whether on a rectangular or triangulated grid, yield sub-Eötvös differences in the gradients when compared to the other methods (except near the edges of the integration area) and they are as efficient computationally as the finite element methods.     

Traction image method for irregular free surface boundaries in finite difference seismic wave simulation

In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem.     

A traveltime reciprocity discrepancy in the Podvin & Lecomte time3d finite difference algorithm

We have found that the extensively used finite difference scheme time3d produces time fields which are 'asymmetric' in the sense that traveltimes computed to the right of the source are faster than traveltimes computed to the left. All finite difference schemes will, as they are approximations to the wave equation, to some extent fail to obey reciprocity perfectly. We show, however, that the errors in time3d may be significant—and unnecessarily large. An asymmetry in the point source initialization has been identified, and after correction time3d produces time fields with an improved reciprocity.     

基于有限差分法的抗滑桩计算机辅助设计

基于地基系数“m-m”法、“m-k”法、“k-k”法的原理,考虑桩顶和桩底边界条件以及桩在滑动面处位移、转角、弯矩和剪力的连续条件,可解得桩身各节点的位移和内力,提出了进行抗滑桩全桩内力计算的有限差分法。根据差分方程并用VB6.0编制了实用的计算程序,既可避免繁琐的查表计算,提高计算速度,又可提高计算精度,直观生动,真正实现了人机交互,在界面的引导下,设计人员可完成全部计算,并绘出内力图形和抗滑桩截面配筋图,使设计更方便快捷,该软件可以极大地提高生产效率,降低工程造价,从而实现抗滑桩的优化设计。最后采用上述方法对某滑坡的悬臂抗滑桩进行了设计与计算。     

VSP资料约束地层反滤波研究

浅地表的低速带和地层滤波作用导致地震信号的能量衰减,地震子波的高频快速衰减,速度频散的子波变形,对地震底层分析带来了严重的影响。地层反滤波算法就是为了消除地层对子波的影响,利用VSP中的直达下行波推算出每个地层滤波算子,计算出各个地层滤波算子的反算子,然后利用这些反算子恢复出地表记录中损失掉的高频成分,减弱地层滤波对波形的影响,增加地表记录的频率成分和波形的准确度。分析实际地震资料表明,经过地层反滤波处理后的记录频率成分明显提升,地质构造更加准确、清晰。     

深层地震勘探激发井深探讨

深层地震勘探野外资料采集的关键是在保证原始资料一定信噪比的基础上,同时兼顾分辨率,激发井深是决定资料品质的重要因素之一。通过对虚反射滤波作用分析,提出了充分利用虚反射能量设计最佳激发井深,增强下传地震波能量的采集方法。理论分析和试验资料表明,在高速层下4～7m的最佳激发岩性深度是深层勘探的理想井深。     

直升机载频率域电磁系统在我国环境和工程勘查中的应用前景

本文根据国际上采用直升机载频率域电磁系统进行环境和工程勘查的应用效果,结合我国引进的航空电磁系统的性能和特点,探讨了采用直升机载频率域电磁系统在我国进行环境和工程勘查的应用前景。     

低频泥石流特征及防治 --以四川汶川县茶园沟为例

低频泥石流指30a以上发生1次的泥石流。尽管发生频率很低，但由于自身的特征及人们认识上的原因，其造l戒的人员伤亡和财产损失却居各类泥石流之首。对低频泥石流的研究在预防和减轻泥石流灾害方面有着重要意义。文章以2003年8月9日发生在四川汶川县茶园沟的泥石流为蒯，采用形态调查法和称重法等方法确定了流量、流速及重度。分析了其物质来源及形成的自然和人为原因，探讨了低频泥石流的特征。针对低频泥石流的低频性、破坏性和隐蔽性特点，同时考虑到低频泥石流沟的防治受保护对象的重要程度及资金投入等因素的限制，提出了进行泥石流危险性分区、搞好小流域水土深持、普及泥石流防灾减灾知识及加强监测和预警预报等措施，这对减轻低频泥石流区的灾害有一定的借鉴意义。