Traction image method for irregular free surface boundaries in finite difference seismic wave simulation   总被引：16，自引：1，他引：16  

Wei Zhang  Xiaofei Chen 《Geophysical Journal International》2006,167(1):337-353

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.  相似文献   

Verification of the motion fields over Southeast Asia by the Florida State University (FSU) superensemble forecast     

Mastura Mahmud  Robert S. Ross 《Singapore journal of tropical geography》2006,27(2):176-190

The Florida State University (FSU) multimodel superensemble forecast is evaluated against several other operational weather models for the Southeast Asia region. The superensemble technique has demonstrated its exceptional skills in forecasting precipitation, motion and mass fields compared to either individual global operational or ensemble mean forecasts. The motion field investigation for the season of 2001 reveals that the superensemble forecasts are closer to the observed data compared to the other global member operational models through its low systematic errors at the 850 hPa level. The FSU multimodel superensemble forecasts exhibit the lowest root mean square errors (RSMEs), the highest correlation against the best observed data and the lowest systematic errors compared to the other operational model members. These forecasts have the potential to provide better daily weather predictions over the Southeast Asia region, particularly during the early northeast monsoon that often causes heavy rainfall in the equatorial part of the Southeast Asia region.  相似文献   

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

Ari Tryggvason  Björn Bergman 《Geophysical Journal International》2006,165(2):432-435

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.  相似文献   

基于有限差分法的抗滑桩计算机辅助设计   总被引：3，自引：0，他引：3  

杨佑发  黄小鹏 《中国地质灾害与防治学报》2006,17(2):95-98

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

基于预测的边缘检测方法   总被引：3，自引：0，他引：3  

王晖  卢健 《武汉大学学报(信息科学版)》2005,30(1):23-26

提出了一种基于预测的、有一定自适应性的边缘检测方法。介绍了其原理和算法实现过程。通过计算两个相邻像素的灰度值的均值、均方差和梯度来预测下一个像素的灰度值，比较预测值和真实值来判断下一个像素是否是边界点。通过与现有同类算法的比较，证明了该算法的优越性。  相似文献   

地震属性分析中水平切片的应用   总被引：6，自引：0，他引：6  

国洪伟  陈勃  刘和芝  刘伟  张秀丽 《物探化探计算技术》2005,27(1):10-15

利用地震数据体中含有的丰富地质信息,可以有效地进行储层预测。层拉平水平切片地层学解释方法是先将地震剖面目的层段做层拉平处理,然后对切片顺地层层面追踪,勾绘出各种地层现象,并以地层学的观点加以解释。该方法在吉林探区扶新南部地区的应用中,对于识别河道的展布趋势取得了较好的效果,而且发现了多期叠置的扇型沉积体。另外,还可指导反演数据体的解释,做到对岩性体从定性到定量的描述。  相似文献   

一种基于核学习的储集层渗透率预测新方法   总被引：2，自引：1，他引：2  

杨斌  匡立春  施泽进  孙中春 《物探化探计算技术》2005,27(2):119-123

基于核学习的支持向量机,是一种采用结构风险最小化原则代替传统经验风险最小化原则的新型统计学习方法,具有完备的理论基础。这里提出了核学习技术在储集层非均质特性描述中渗透率参数预测的新用途。在复杂地层中,基于支持向量机的智能和自适应模式识别能力而建立了常规测井多参数信息输入的渗透率预测模型,然后对实际油田储集层渗透率进行了预测。与常规线性回归模型预测结果相对比,所提出的方法更易于使用,很少受不确定因素的影响,并具有较强的信息整合能力以及更高的预测准确性和可信度。  相似文献   

碳酸盐岩溶洞型储层反演方法     

王晓梅  唐文榜  李学著  杨勇 《物探与化探》2005,29(1):44-46

常规声波反演的方法原理和反演技术以层状介质为基础,其研究目标多是层状储层。碳酸盐岩溶洞型储层具有非规则形态、非均匀散布的特征,常规声波反演技术有其不适应之处。笔者研究的测井约束多重反演技术,解决了非层状、非均匀储层的地震反演问题,得出反映碳酸盐岩溶洞储层信息的波阻抗数据体,提取了突出溶洞型储集体低速特征的差异波阻抗,为寻找碳酸盐岩溶洞型油气藏提供了可靠的依据。  相似文献   

现阶段地表油气化探技术难题与发展对策   总被引：2，自引：1，他引：2  

姚俊梅 《物探与化探》2005,29(6):477-480,486

虽然地表油气化探方法技术已取得显著的应用效果, 但是由于成功率计算标准不同及方法技术本身的缺陷, 导致人们对方法的认同度较低。当前, 地表油气化探技术的发展, 一方面要以可靠的实验依据、严谨的数理推导及理论模型应对来自方法技术本身的挑战, 另一方面也要利用其快速、廉价和直接的优势, 加强与地质、遥感、非震物探等方法综合应用。  相似文献   

低频泥石流特征及防治 --以四川汶川县茶园沟为例   总被引：2，自引：0，他引：2  

张金山  沈兴菊  谢洪 《中国地质灾害与防治学报》2005,16(3):37-41

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