一种浸入边界起伏地表地震波场正演模拟方法技术

本发明专利技术提供一种浸入边界起伏地表地震波场正演模拟方法，其特征在于，包括以下步骤：S1)将待模拟地层建立数值模型，确定其隐层数、鬼点坐标、镜像点坐标及镜像点物理参数；S2)将整个计算区域依时刻进行声波方程求解，得到计算区域各时刻各个点的波场值，其中，在每一时刻有限差分计算时，使用修正后的鬼点波场值代替常规有限差分计算中的鬼点波场值；S3)重复步骤S2)至完成全部时刻的波场模拟，并保存，得到地震波正演模拟结果。本发明专利技术在较大的网格采样情况下仍然能够得到与常规方法精细网格非常相似的波场模拟结果，由此表明本发明专利技术的。有效减少计算量大带来的不便。

【技术实现步骤摘要】
一种浸入边界起伏地表地震波场正演模拟方法
本专利技术涉及地质领域，具体的说，涉及一种浸入边界起伏地表地震波场正演模拟方法。
技术介绍
地震波正演模拟对分析地球深部的结构以及寻找油气等国家战略资源具有重大的意义。由于地质结构的复杂性(盐丘、风化带、山麓)给陆地地震勘探勘探带来巨大困难。其中一个重要因素就是起伏地表使得地震波场变得更为复杂，给油气等国家战略资源的勘探带来了极大的挑战。在地震波正演模拟计算中，起伏地表(自由表面)的几何形态对地震波传播具有显著影响。由于人工矩形网格离散，正演模拟建模过程中真实的处于分数网格上的地表界面，会被截断在整数网格处，因而在边界处形成“阶梯状”的数值建模结果。地震波场正演模拟过程中会在这些人工截断的“阶梯状”端点处产生很多虚假的绕射波，严重干扰波场正演模拟的结果，并影响地震波反演和成像(如全波形反演)的结果。研究起伏地表波场模拟方法得到了全球学者的关注，目前主要的起伏地表波场模拟方法从大的方面可以分为：有限差分法、伪谱法、有限元法、谱元法等。有限差分法和伪谱法在矩形网格上进行计算，因而在处本文档来自技高网...

【技术保护点】
1.一种浸入边界起伏地表地震波场正演模拟方法，其特征在于，包括以下步骤：/nS1)将待模拟地层建立数值模型，确定其隐层数、鬼点坐标、镜像点坐标及镜像点物理参数；/nS2)将整个计算区域依时刻进行声波方程求解，得到计算区域各时刻各个点的波场值，其中，在每一时刻有限差分计算时，使用修正后的鬼点波场值代替常规有限差分计算中的鬼点波场值；/nS3)重复步骤S2)至完成全部时刻的波场模拟，并保存，得到地震波正演模拟结果。/n

【技术特征摘要】
1.一种浸入边界起伏地表地震波场正演模拟方法，其特征在于，包括以下步骤：
S1)将待模拟地层建立数值模型，确定其隐层数、鬼点坐标、镜像点坐标及镜像点物理参数；
S2)将整个计算区域依时刻进行声波方程求解，得到计算区域各时刻各个点的波场值，其中，在每一时刻有限差分计算时，使用修正后的鬼点波场值代替常规有限差分计算中的鬼点波场值；
S3)重复步骤S2)至完成全部时刻的波场模拟，并保存，得到地震波正演模拟结果。


2.如权利要求1所述的一种浸入边界起伏地表地震波场正演模拟方法，其特征在于，所述步骤S1)包括：
S11)：根据有限差分空间差分阶数确定隐层数；
S12)：基于输入的不规则地表边界获得各隐层鬼点以及对应地表上的法向投影点的空间坐标；
S13)：在获得鬼点和法向投影点的空间位置的基础上，找到地表以下镜像点的坐标使法向投影点位于鬼点和镜像点连线的中点；
S14)：将镜像点处的物性参数设置于相对应的鬼点处。


3.如权利要求2所述的一种浸入边界起伏地表地震波场正演模拟方法，其特征在于，所述步骤S14)中的物性参数为速度或密度。


4.如权利要求1所述的一种浸入边界起伏地表地震波场正演模拟方法，其特征在于，某一时刻所述修正后的鬼点波场值通过如下步骤获得：
S21)：将该时刻的鬼点波场赋初值设为0；
S22)：在模型区域内进行波场空间对称插值，获得位于分数网格点上的各镜像点处...

【专利技术属性】
技术研发人员：吴迪，姚刚，李翔，
申请(专利权)人：中国石油大学北京，
类型：发明
国别省市：北京;11