单台地震自动定位网格搜索法及其MATLAB试验

提出了一种利用单个地震台站三分向数字地震波形记录进行单台地震自动定位的网格搜索法。该方法以台站为中心,建立辐射状网格,按可变震中距步长和方位角步长进行搜索。对一个检测到的地震事件记录波形,对每个震中距,计算预测波形特征函数与记录波形特征函数的互相关函数值,选其最大值对应的震中距;对每个方位角,合成沿径向和沿横向的波形分量,并计算P波段的径向绝对积分值与横向绝对积分值之差,选其最大值对应的方位。对于选定的方位,如果合成的径向分量与记录垂直分量乘积在P波段的积分值为负,则取该方位角为震中相对于台站的方位角,否则该角要加上180°。在MATLAB环境中对该方法的基本步骤进行了试验,初步结果表明该方法是可行的。     

计算断层自发破裂过程的离散网格划分研究

离散化网格的空间步长选取在各种数值算法中都是一个很受关注的问题,在全空间均匀介质模型和简单离散化方案下,利用边界积分方程方法研究在自发破裂求解过程中动力学参数组合Dc和Te (Dc为临界滑动弱化位移,Te为有效的断层破裂面上的初始应力)对计算网格划分的影响,初步得到了Dc和Te参数空间中的有效计算网格的选取规律,对合理、有效地运用边界积分方程方法计算地震震源的破裂过程具有重要的指导意义.      

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

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

小尺度网格层析速度建模方法在低幅构造区时深转换中的应用

渤海M油田区浅层多发育低幅构造,低幅构造区的地震数据对常规层析成像不够敏感,导致常规速度反演方法很难保证低幅构造速度建模的准确性,从而给时深转换工作带来较大误差.本文提出一种小尺度网格层析速度建模方法,该方法基于低幅构造的薄层信息,在有薄层信息的地层加密网格分布,从而实现根据地质条件划分网格的策略,进而通过解析矩阵方程实现速度模型的迭代更新.通过实际数据应用证明,对于渤海M油田利用本文提出的高精度网格层析速度建模方法是有效可行的,该方法能够有效提高低幅构造区速度建模的精度,利用该方法得到的速度模型能够大幅提高时深转换的精度,从而指导构造解释及成图工作.     

基于时空双变网格的起伏地表变坐标系正演模拟方法（英文）

坐标变换法通过将物理空间的曲网格映射为计算空间的矩形网格,将起伏地表转化为水平地表,同时将物理空间的波动方程转化为计算空间的波动方程,在计算空间完成数值模拟,坐标变换的方法对处理起伏自由边界具有较好的适应性和应用效果。本文在传统坐标变换方法的基础上,根据计算区域速度差异采用不同的网格大小和采样时间步长,提出了一种基于时空双变网格的起伏地表坐标变换正演模拟方法。在编程实现算法的基础上,通过典型模型波场模拟试算结果分析可知:(1)变网格方法与常规方法波场模拟误差在0.5%左右;(2)变网格方法计算效率视不同的变网格区域面积及变网格大小可提高几倍量级,在本文模型和计算参数下提高约5倍。(3)在满足模拟精度及频散条件要求下,变网格方法较全局细网格算法能显著节约计算内存。为此,针对起伏地表数值模拟,本文方法具有较高的模拟计算精度和一定的适应性。     

网格剖分及其精度和计算量分析

网格剖分程度直接影响着地震波正演数值模拟的计算精度及其计算量.以均匀介质模型为例,分析不同网格大小对波场模拟精度和计算量的影响,得出精细化网格剖分是实现高精度地震波正演模拟的有效方法,然而其计算量较大.以均匀倾斜介质模型为例,探讨倾斜地层网格剖分问题,数值实例分析不同震源频率对不同网格剖分方案引起的波场传播精度的影响,...     

基于CUDA的高阶旋转交错网格有限差分法的弹性波正演模拟

应用高阶旋转交错网格的有限差分法,求解了基于一阶的速度-应力弹性波方程,推导了时间二阶精度空间2M阶精度的高阶有限差分离散格式.为了提高弹性波正演模拟效率,论文探讨了常规的CPU串行计算、基于GPU平台的节点并行计算、基于GPU平台的分量、节点同时并行的3种正演模拟计算模式,给出了不同算法的流程示意图.通过对比3种算法的耗时及后两种算法的加速比,可清晰判定,基于GPU平台的节点并行算法较CPU平台的串行算法有明显加速,最大加速比可达到60倍,而基于分量、节点同时并行算法的最大加速比最大可到180左右,加速效果更为明显.最后应用Marmousi模型的15000次时间步迭代的正演计算,基于分量、节点GPU并行算法大约是节点并行计算耗时的1/3,说明了该并行设计能有效提高弹性波正演模拟的效率,节约计算耗时.     

等效交错网格高阶有限差分法标量波波场模拟

密度是岩石物性参数的重要组成部分,对储层评价、岩性解释和油藏描述等具有重要作用,因此实际资料应用中含有密度信息的正演建模是必不可少的.作为波动理论应用最广泛的实际资料处理技术——逆时偏移和全波形反演,其正演建模一般采用经典的二阶标量波方程,且不考虑密度的空间变化.而一阶速度-应力控制方程通过交错网格有限差分法正演建模,其精度高于二阶方程但计算成本过高,很难在三维实际资料中应用.鉴于非均质正演在实际生产中所面临问题,本文借用交错网格的思想,充分考虑空间变密度对波传播的影响,提出基于等效交错网格的高阶有限差分建模算法并应用于非均质标量波正演.文中从数学上给出一阶方程和二阶方程在数值模拟中的等价性证明,并分析了震源、边界条件和稳定性.最后通过简单层状模型进行数值测试,对比不同建模方法以验证新方法的准确性;利用高阶精度方法对Sigsbee-2a模型做正演测试,验证本文提出的方法的稳定性.     

多极子随钻声波测井波场模拟与采集模式分析

应用三维非均匀交错网格有限差分程序模拟了多极子随钻声波测井的波场,考察了硬地层中单极子随钻声波测井响应及软地层中声源频率对多极子随钻声波测井的影响;计算了接收阵列波形的时域相干谱及频散特征,结果表明随钻测井仪器居中时,在硬地层中可以直接利用单极子声系得到地层横波速度,但地层纵波受到钻铤波的极大干扰.而在软地层中可以利用单极子声系直接测得地层纵波速度,但地层横波速度则需要使用工作在较低频率下的四极子随钻声波仪来测得,此时测井频段又会落入钻井噪声频带而受到影响.为此本文提出数值模拟手段结合物理实验的方法得到钻铤波的影响,在实际测井中通过信号处理方法而不是用在仪器本体上刻槽的方式来更完全地消除钻铤的影响.为了避免钻井噪声的影响,本文提出六极子随钻声波仪更适合软地层的横波测量.本文用数值实验结果来说明这些方法的可行性.     

波前构建法中的网格点定位(英文)

针对波前构建法中非规则的波前四边形网格与规则的矩形网格节点之间的相对定位问题,对比分析了在计算机图形学中常用的四种网格定位方法(叉乘判断法、夹角之和检验法、射线法(交点计数检验法)以及符号比较法)在波前构建法中的具体实现问题,通过对定位方法稳健性分析及计算实例表明,矢量叉乘判断法在实现规则的矩形网格节点与非规则的波前四边形网格相对定位研究中,其稳健性和计算效率足最高的,从而为下一步做偏移或者反偏移提供了精确的网格点属性值.