国家自然科学基金(41174119) 作品数:7 被引量:53 H指数:3 相关作者: 周辉 王颖 刘万金 袁三一 刘文岭 更多>> 相关机构: 中国石油大学(北京) 中国石油大学(北京) 中国石油 更多>> 发文基金: 国家自然科学基金 国家重点基础研究发展计划 国家科技重大专项 更多>> 相关领域: 天文地球 更多>>
基于新的差分结构的时-空域高阶有限差分波动方程数值模拟方法 被引量:10 2016年 传统的高阶有限差分波动方程数值模拟方法采用高阶差分算子近似空间偏导数,能有效抑制空间频散.然而,传统的有限差分法仅采用二阶差分算子近似时间偏导数,这使得地震波场沿时间外推的精度较低.当采用较大的时间采样间隔,传统的有限差分法模拟波场会出现明显的时间频散,甚至不稳定.本文基于新的差分结构和中心网格剖分,发展了一种空间任意偶数阶精度、时间四阶和六阶精度的时空域有限差分方法.基于对离散后的频散关系进行泰勒展开,本文推导了时空域高阶有限差分算子的差分系数.相速度分析表明时间四阶、六阶精度的差分方法能显著地减小传统时间二阶精度差分方法的时间频散.在相同的精度下与传统差分法比较,本文发展的时间四阶、六阶有限差分方法的计算效率比传统方法高.均匀和非匀均介质中的波场数值模拟实验进一步证实本文研究的时空高阶有限差分方法的优越性. 张保庆 周辉 周辉 盛善波关键词:时空域 高阶有限差分 数值模拟 Huber-Markov随机场自适应边缘保护贝叶斯地震资料叠前反演(英文) 被引量:2 2013年 由于地震数据频带范围的有限性以及噪声、正演近似等因素的影响,地震反演是一个高度不适定的反问题。为了得到稳定唯一的解,通常在反演目标函数中加入光滑化约束,常用的是Tikhonov正则化约束。这种正则化对解是全局光滑的,在构造边缘位置会产生模糊。为了解决边缘模糊的问题,同时保证抗噪的作用,本文采用Huber-Markov随机场约束边缘保护的方法求解叠前纵横波速度和密度三参数反演问题。对于待反演的参数,通过Markov邻域系统建立起纵向和横向上的约束。通过Huber边缘惩罚函数阈值的合理设置,达到在层内区域采用二次型函数抑制噪声得到平滑的结果;在边缘位置则采用线性函数,避免边缘的模糊化。通过无噪和有噪合成地震记录的反演,验证了方法的正确性,分析了正则项与反演效果的关系。 田玉昆 周辉 陈汉明 邹雅铭 关守军关键词:地震反演 马尔可夫 TIKHONOV正则化 自适应 贝叶斯 贴体坐标系二维声波方程SBP有限差分法的稳定性分析 被引量:1 2016年 有限差分方法(Finite Difference Method,FDM)是波动方程正演数值模拟领域应用最为广泛的方法之一,然而,当模拟区域不规则或者地表起伏不平时,规则网格有限差分法求解波动方程会产生阶梯状近似,影响模拟的精度。借助贴体网格技术,将不规则的物理区域转换为规则的计算域,给出了贴体坐标系下的二维声波方程及其二阶精度的分部求和(Summation by Parts,SBP)有限差分离散格式,采用Fourier谱分析方法分析了该离散格式的稳定性,得到了贴体网格二维声波方程SBP有限差分方法的稳定性条件。数值实验结果表明:1当时间采样间隔的选取满足稳定性条件时,贴体网格SBP有限差分的数值计算过程是稳定的;2与贴体网格中心差分方法相比,贴体网格SBP有限差分方法的稳定性更好。 王颖 周辉 盛善波关键词:声波方程 中心差分 提高任意广角波动方程计算效率及压制倏逝波干扰方法 被引量:2 2014年 任意广角波动方程(Arbitrarily Wide-angle Wave Equation,AWWE)是一种高精度的时空域单程波动方程,与其他时空域单程波方程相比,具有更大的成像倾角,因此是偏移成像的有力工具之一.本文改进了现有的AWWE有限差分计算方案,避免了原计算方案中对矩阵反复求逆的过程,同时减少了矩阵相乘的次数,显著地提高了AWWE数值计算的效率.数值实验证实,采用改进之后的计算方案,其计算结果与原方案完全一致.此外,本文采用了一种方法,在有倏逝波干扰的情况下,能在频率-波数(f-k)域区分出倏逝波与非倏逝波,为设计f-k滤波器提供合理的视速度门槛信息.偏移实验证实,在深度延拓时仅仅进行数次f-k滤波就能很好地抑制倏逝波干扰,提高偏移成像质量. 周辉 陈汉明 李绪宣 范廷恩关键词:有限差分 倏逝波 补偿吸收衰减的高角度叠前深度偏移方法(英文) 2012年 实际介质的吸收衰减使地震波能量减弱,主频降低,影响分辨率。在偏移过程中,可以实现吸收衰减的补偿,提高地震资料的分辨率。本文基于高角度时间-空间域单程波方程,推导出能补偿介质吸收衰减的叠前深度偏移波场外推公式。设计了一个复杂构造的地质模型,用声波方程模拟了无吸收衰减和有吸收衰减介质存在时的两种地震炮集记录。用没有补偿吸收衰减和本文讨论的补偿吸收衰减的偏移方法分别处理了两种地震资料,得到无吸收衰减、吸收衰减未得到补偿和得到补偿的三个偏移叠加剖面。补偿后的偏移剖面与没有吸收衰减的地震记录的偏移剖面几乎相同,而有吸收衰减的地震记录用无吸收衰减补偿的偏移方法处理,由于介质的吸收衰减,使吸收衰减层下方的反射信号衰减,主频降低,能量减弱。实际资料的处理结果表明,补偿吸收的偏移剖面的分辨率得到了明显改善。 周辉 林鹤 林鹤 王颖关键词:单程波方程 叠前深度偏移 谱反演在地震属性解释中的应用 被引量:28 2013年 任何一个反射系数序列都可以分解成奇、偶分量,奇分量不利于检测薄层,而少量的偶分量可以提高薄层的分辨能力,谱反演的实质就是利用偶分量在厚度趋于零时的有效干涉提高地震资料的分辨率。谱反演是在谱分解的基础上,通过反演方法使频率域目标函数达到极小而反演出反射系数、薄层厚度,在没有噪声和地震子波已知的情况下,可以识别任意薄层,并精确地反演出反射系数值。其主要特点是只采用部分频谱资料就可反演稀疏反射系数或层厚。文中利用谱反演方法对大庆油田的三维地震数据体进行了反演,得到反射系数体,并用该反射系数体与给定的宽频带子波褶积,得到提高分辨率后的地震数据体,再对其进行地震属性分析。与常规资料的地震属性分析结果相比,应用谱反演方法可提高识别沉积相和小断层的能力。 刘万金 周辉 袁三一 刘文岭关键词:地震分辨率 反射系数 地震属性分析 时间域黏滞波动方程及其数值模拟新方法 被引量:10 2016年 基于常Q模型的应力—应变关系,推导了分数阶拉普拉斯算子黏滞波动方程的一阶速度—压力形式,与现有的一阶速度—压力—应变形式相比,新推导的方程形式更简单,数值模拟时耗费内存更少。针对新推导的一阶速度—压力黏滞波动方程,采用交错网格伪谱法进行数值模拟,并利用卷积型完美匹配层(CPML)作为吸收边界,压制截断边界的反射。数值模拟实验证实,新推导的黏滞波动方程能很好地描述地震波在黏滞介质中的衰减和频散,交错网格伪谱法和CPML的组合是一种高效的数值模拟方法。 罗文山 陈汉明 王成祥 周辉 王狮虎关键词:时间域 伪谱法