向读者推荐

编者按：《地学前缘》编辑部将不定期推出“向读者推荐”栏目。根据一些论文评价指标和专家评审意见,重点推出一些文章的评价意见、主要内容和创新点提示,以期获得读者更多的关注,提升学术论文的影响力,更好地发挥学术导向作用,促进论文成果为科技创新服务。     

岩土体渗流自由面问题的重心插值无网格方法

岩土体的渗透破坏、地下工程的防渗设计等无不与渗流计算有关。针对渗流自由面问题,提出一种重心拉格朗日插值的配点型无网格方法。由于渗流自由面问题的求解区域是不规则区域,该方法通过将不规则求解区域嵌入一个正则矩形区域,在正则区域上采用重心拉格朗日插值近似未知函数,利用配点法离散渗流问题的控制方程,将重心拉格朗日插值的微分矩阵离散成代数方程表达的矩阵形式。将自由面上的边界条件通过重心拉格朗日插值离散,通过置换方程法和附加方程法施加边界条件,利用正则区域上的重心插值配点法,通过迭代确定最终自由面的位置。数值算例表明所提出的无网格方法对于求解渗流自由面问题的正确性和高精度。     

改进的均衡滤波器在位场数据边界识别中的应用

边界识别是进行位场数据解释时一项必不可少的任务.现有的边界识别滤波器存在识别边界发散和深部地质体边界模糊的缺点.本文提出增强型均衡滤波器,可有效地改善上述缺点.该滤波器是利用不同阶导数之间的组合来进行地质体边界的识别,并在运算中引入一种计算高阶垂直导数的稳定算法.通过理论模型试验证明增强型均衡滤波器能使浅部与深部地质体的边界同时清晰地显示,且相对于其它边界识别滤波器能更加准确和清晰地识别出地质体的边界.最后将增强型均衡滤波器应用于实测位场数据的解释,根据其识别结果可容易地划分出断裂的水平位置及不同地层之间的界线,并能发现更多的细节信息.     

黏弹TTI介质中旋转交错网格高阶有限差分数值模拟

以Carcione黏弹各向异性理论为基础,给出了适用于黏弹性具有任意倾斜对称轴横向各向同性介质(黏弹TTI介质)的二维三分量一阶速度-应力方程,采用旋转交错网格任意偶数阶精度有限差分格式求解该方程,并推导出了二维黏弹TTI介质完全匹配层(PML)吸收边界条件公式和相应的旋转交错网格任意偶数阶精度有限差分格式,实现了该类介质的地震波场数值模拟.数值模拟结果表明:该方法模拟精度高,边界吸收效果好,可以得到高精度的波场快照和合成记录;并且波场快照和合成记录能较好地反映地下介质的各向异性特征和黏弹性特征.     

选择性滤波同位网格有限差分法在地震波数值模拟中的应用.

引入计算空气声学领域的选择性滤波同位网格有限差分算法（SFFD法）用于二维地震波数值模拟．SFFD法使用经过优化的11点DRP同位网格差分格式,对空间一阶导数进行离散近似,同时采用选择性滤波方法来消除同位网格差分所产生的格点高频振荡,它既提高了数值模拟的精度, 又保证了求解过程的稳定性．数值实验结果表明,SFFD法能够达到O(Delta;x8, Delta;t4)阶交错网格算法同样的精度,同时该方法还具有很强的适应性,能够应用于存在着强泊松比差异的介质模型中,完整地模拟地震波传播过程中各类型的波场,并且对复杂非均匀介质的适应能力也很好．此外,由于避免了交错网格算法在曲线坐标系和一般各向异性介质的数值模拟时所需进行的复杂的插值运算, SFFD法在这些问题上也有着很好的应用前景．      

用于地震观测的瞬态滤波研究

目前,地震数据采集器使用的最后两级数字滤波器相位有最小相位和线性相位两种,相位特性直接影响着记录的地震数据波形质量.最小相位滤波的优点是:地震波形没有“前缀”波形出现,不会干扰分析人员对震相的判读;其缺点是:通带边沿频率附近不同频率分量群时延差异较大,导致最小相位滤波会带来波形失真,对初动半周期的波形影响显著.     

强震人员损失回归预测方法

在回顾近年对地震人员损失预测的研究工作的基础上,根据USGS的PAGER系统中地震人员损失回归模型,利用发生在1970-2008年间的128条地震现场调查灾害记录建立了我国东、西部区域适用的地震人员损失预测模型;并利用发生在1980-2007年间的234条地震损失记录建立了地震人员重伤数与人员死亡数的回归关系;最后利用模型对2008年四川汶川地震的人员损失进行评估计算,得到了与实际震害损失相接近的评估结果。     

地震活动性预测并划分未来地震发生的潜在区域方法研究

采用网格地震方法,选取面积百分比(33%),确定圈闭区域的半径,统计后期地震发生在圈闭区域内的概率。在我国4个不同构造环境和地震活动水平的地区,分别统计了后期地震发生在圈闭区域内的概率。结果表明,4个地区后期地震发生在圈闭区域的概率远远大于面积百分比。说明未来地震是高度地聚集在过去地震发生的区域内,采用地震活动性预测并划分未来地震发生的潜在区域的方法是可行的。     

Steerable Pyramid分解在储层断裂检测中的应用

在隐蔽油气藏勘探中,正确刻画储层断裂对于油气勘探和开发都有着至关重要的作用.储层断裂在地震数据中表现为边缘特征.但现今多尺度边缘检测的基本理论和方法都有其本身的局限性.本文详细阐述了Steerable Pyramid分解用于储层断裂检测的基本原理和方法.输入的三维地震切片经过Steerable Pyramid分解后,不仅可以分析不同尺度不同方向上断裂的特征及走向,还可以通过重构得到经S函数控制的断裂信息,该断裂信息不仅仅是断裂的边缘检测,而是增强的断裂本身.将该算法应用于实际地震切片的断裂检测,得到了清晰的断裂信息,完整细致地展现了这些地质类型的特征,为构造解释、地震相表征及储层预测提供了有用信息.     

提高规则网格最短路径方法反射波走时计算精度的走时校正技术

最短路径射线追踪方法是计算地震波走时的主要方法之一,该方法基于惠更斯原理和费玛原理,具有稳健、适于复杂介质模型的优点.为处理方便,最短路径方法中的介质模型通常以规则网格进行剖分,界面节点(界面与网格的交点)以其邻近的模型单元节点(即边界单元节点)近似表示.界面近似将导致计算误差,对于反射波尤为严重.反射波的走时精度可通过减小网格的尺寸提高,但这样会大大增加计算时间,为高精度和高效率地计算地震反射波走时,我们提出了一种基于规则网格的走时校正技术.地震波传播至或起始于边界单元节点的走时校正为地震波传播至或起始于该边界单元节点所对应的界面节点的走时.数值模型计算结果表明,走时校正方法可使反射波的走时精度提高约1～2个数量级,而其计算时间则和常规算法基本上在相同量级.