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

考虑滑脱效应和温度场的煤层气渗流数值模拟

由于煤层气的解吸热效应,煤层气的运移过程是一个非等温过程。因此,温度场对煤层气渗流有着重要的影响。前人在研究煤层气渗流规律的研究中并没有同时考虑滑脱效应和温度场因素的影响,而在实际的深部开采中综合考虑滑脱效应和温度场的影响对研究深部煤层气运移规律有着重要的意义。因此本文建立了考虑滑脱效应及温度场的煤层气渗流数学模型,利用有限元数值方法研究了考虑滑脱效应和温度场耦合的煤层气运移规律;研究了考虑滑脱效应和温度场对压力场分布的影响;对考虑滑脱效应及温度场因素的煤层气产量进行了预测。得出了随着温度的升高煤储层的压力在降低,温度的升高对煤层气的产量有着负面的影响这一重要结论。     

吸水剖面综合解释方法及应用

吉林油田同位素吸水剖面测井采用~(113)In同位素追踪技术,根据同位素追踪过程解释基本可以判断配水层段吸水情况。但是在注水井存在沾污、大孔道以及串槽等影响时,吸水面积确定比较困难。在同位素测井中增加井温、流量参数,通过多参数综合解释,不仅可以对沾污影响进行合理校正,确定准确的小层吸水量,而且能够正确判断各级封隔器、配水器的工作情况,在地层存在大孔道的情况下,确定地层的吸水面积。吸水剖面综合解释方法在吉林油田已普遍应用,见到明显的效果。     

长江中下游夏季高温灾害机理及预测

利用我国1961-2003年夏季（6—8月）高温资料，建立长江中下游地区主要城市强高温及高温过程较完整的时间序列，探讨了该地区主要城市高温气候特征。分析该地区南京、杭州、南昌等城市夏季高温灾害机理，东亚副热带高压是造成长江中下游地区城市夏季高温的主要影响系统。在此基础上用均生函数-最佳子回归集构造预测模型，预测夏季月高温出现日数，通过42a高温资料预报检验，有较好的预测效果，值得在业务中应用。     

非加温型四连栋塑料温室内外温湿度关系研究

根据1999年12月～2002年8月典型天气（晴天、多云、阴天）下非加温型四连栋塑料温室中间1．5m高度气温、相对湿度的观测数据，应用逐步回归分析方法建立了冬季、春秋季、夏季典型天气下温室内日平均气温、日最高气温、日最低气温、日平均相对湿度与气象站气象要素间的关系式，为温室蔬菜品种筛选、蔬菜标准化栽培、无公害蔬菜生产和病虫害防治提供小气候方面的技术数据。     

WAFS产品中GRIB资料中国区产品评估

为了给使用WAFS产品的用户提供量化参考依据，选取WAFS产品中常用区域和层次的GRIB数据，利用由国家气象中心提供的风、温客观分析网格点资料，与WAFS中同时刻的预报场产品（风、温网格点资料），用均方根误差进行数字化形式分析比较。结果表明：WAFS提供的风、温预报，通常是短时效的风、温预报比长时效的风、温预报更接近客观分析场；低层的预报比高层的预报更接近客观分析场；风的预报以v矢量的预报优于u矢量的预报；风的误差主要来源于u矢量的误差。