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编制了实用的计算程序,既可避免繁琐的查表计算,提高计算速度,又可提高计算精度,直观生动,真正实现了人机交互,在界面的引导下,设计人员可完成全部计算,并绘出内力图形和抗滑桩截面配筋图,使设计更方便快捷,该软件可以极大地提高生产效率,降低工程造价,从而实现抗滑桩的优化设计。最后采用上述方法对某滑坡的悬臂抗滑桩进行了设计与计算。     

华南地区晚二叠和晚三叠世煤中稀土元素特征研究

在华南地区贵州六盘水、江西乐平、露庆和湖北黄石矿区晚二叠和晚三叠世煤层的稀土元素ICP—MS分析数据的基础上，系统研究了研究区煤中稀土元素的含萤分布特征、赋存特征、配分模式以及地质控制因素。含量分布表明晚二叠世稀土元素含量总体高于晚三叠世；煤中稀土元素与灰分成正相关关系，部分煤样品由于有富稀土元素的重矿物组合的存在而出现稀土元素的异常高值；同一矿区、同一时代煤中稀土元素的配分模式具有一定的相似性，而不同地区同一时代和同一地区不同时代煤中稀土元素的配分模式具有较大差异；煤的沉积环境，煤中无机矿物组成以及海水的影响是控制煤中稀土元素含量和配分模式的主要地质因素。     

淮北煤田孙庄勘探区煤层气赋存规律及影响因素分析

通过对淮北煤田孙庄勘探区地质构造、沉积特征、煤层煤质及煤层气赋存特征的分析，研究认为有利于本区煤层气赋存的主要因素是较大的埋藏深度、较高的壳质组＋镜质组、较高的变质程度、较大的煤层厚度、低透气性的盖层、较缓的地层倾角和封闭性较好的构造条件，推测煤田深部可能具有一定工业价值的煤层气资源。     

基于GIS/MIS/OA的短周期工业原料林经营模式研究

提出了短周期工业原料林经营模式的概念,探讨了短周期工业原料林经营模式的体系结构、子系统功能和主要林业数表类型,说明其与GIS/MIS/OA系统集成的必然趋势,并以岳阳纸业股份公司的纸浆原料林为例进行了实例分析。     

一种实现依比例街区自动合并的算法

针对依比例街区的图形特征，介绍了一种实现街区自动合并的有效算法，即利用数学形态学和模式识别对栅格数据进行追踪，然后提取街区轮廓线，编辑后就得到最终合并的街区。     

主砂体及其构造改造样式的盆地产铀远景评价判据

本文在介绍并分析美国和前苏联总结的砂岩型铀矿区域评价判据的基础上, 引入当今碎屑岩沉积学和盆地地球动力学方面的若干新进展, 提出以主砂体及其随后经历的构造改造样式相结合的盆地产铀远景评价判据和以分析区域构造活动史和原型盆地演化相结合的盆地产铀远景评价模式。最后以潮水盆地为例, 具体分析了盆地各个时代盖层容存后生铀矿化的可能性, 指出盆地中的主要找矿目标层位和远景成矿地段。     

泥石流汇流槽可靠度分析

汇流槽是泥石流治理工程中常用的工程措施之一。文章基于汇流槽的倾倒破坏、滑移破坏及地基破坏模式，通过对汇流槽设计影响强烈的岩土参数、几何特征、荷载因素等敏感因子的求解方法及分布特征分析，从汇流槽的抗倾稳定、抗滑稳定和基底应力要求三方面建立极限状态方程。对于每个极限状态方程，在泰勒级数对非线性极限状态方程线性化后，运用一次二阶矩的理论，先假定失效点P^＊，求解出相应的均值和方差。然后根据可靠的定义，得到相应的可靠度指标，通过迭代法求解真正失效点的可靠度指标βi。最后根据3个状态方程相互独立的假定，得出整个结构的可靠度指标β。该方法在平川泥石流防治工程中得到了成功的运用。这对于可靠度设计方法在泥石流防治工程中的运用进行了有益的探讨。