复杂界面条件下的地震全波形反演

来源 :中国石油大学(北京) | 被引量 : 0次 | 上传用户:winningking
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文针对复杂地表起伏(高差达到千米)的地质模型,引用间断有限元方法,它利用非结构化的三角形网格剖分地质模型进行波动方程模拟。传统的有限差分方法因为使用矩形网格剖分地址模型进行波动方程模拟,仅能够适应简单起伏(高差在百米附近)的地质模型,当面对起伏剧烈的地质模型时,会因为矩形网格对起伏曲线的阶梯状近似而带来人为绕射的问题,造成模拟结果出现干扰波,带来误差。传统的有限元方法是适用于起伏地表模型波场模拟的一种方法,该方法利用非结构化的三角形网格可以很好的拟合起伏曲线,从而更好的剖分模型。但是因为该方法带来了巨大的计算量及需要极大的计算存储量的弊端,因此模拟中难以真正使用。相较而言,间断有限元法可以解决以上两个缺点,该方法既可以很好的拟合起伏曲线,也因为其间断性可以降低计算存储量,且利用并行计算可以更好的提高计算效率,从而使得该方法可以更适合于波动方程模拟。本文面对复杂地表起伏模型中波场模拟的边界吸收问题,提出PML吸收边界条件的优化方法。地表起伏模型会造成吸收边界厚度不均匀的问题,在起伏剧烈的情况下显得尤为明显。传统PML吸收边界条件,面对上述问题会因为吸收效果差而出现反射能量接近有效波场能量的边界反射,从而影响模拟结果。优化后的PML边界条件具有两个优势:一是可以根据需要设置起伏地表的边界区域,边界厚度可以调整;二是面对不同厚度的吸收边界,依然可以有效吸收,边界反射能量极弱。本文选择全波形反演方法进行波动方程反演,相较于仅仅利用地震数据中单一信息的反演方法,如相位反演、走时反演及振幅反演等,该方法利用了更多的动力学信息,得到的反演结果分辨率更高,但同时也带来反演耗时的问题。在学习前人研究的基础上,本文将选择采用网格转换的方式,把三角形网格转换到矩形网格,使得可以把间断有限元法正演程序应用到全波形反演程序当中。此外,将采用网格分区生成数据,同时在程序中使用MPI与Open MP的并行模式,以提高全波形反演的运行效率。
其他文献
相比于基于Zoeppritz方程及其近似式的叠前反演方法,基于反射率法的叠前反演在对薄互层储层以及宽入射角的地震记录的反演分析中,具有更大的优势。但是,反射率法反演方法在对实际地震资料反演时,仍然具有不足之处。首先在贝叶斯反演理论中,引入的正则化项要涉及到弹性三参数的统计相关性,在噪音比较强时,这种相关性会使得反演结果偏离地震记录的约束,使得结果分辨率下降。本论文引入了去相关算法,通过去除初始模型
地震波形是振幅、频率、相位等地震属性的综合反映。地震波形分类是常用的储层预测方法。地震波形聚类是地震波形分类中一种,它只依赖地震数据本身而不依靠先验信息,使得相同分类相的地震数据聚在同一类簇下的解释工作。K-Means(K均值)是一种基于欧式距离的常用的地震波形聚类方法,K表示需要分类的数目。本文提出了一种基于DTW(Dynamic time warping动态时间规整)算法的KMedoids(K
随着油气勘探程度的日益加大,地震勘探任务也面临着更加复杂的情况,薄储层油气藏的预测问题变得越来越重要,对于厚度小于地震数据传统分辨率的薄储层的准确预测成为了地震勘探领域的重点之一。基于这一实际情况,论文详细研究了地震数据提高分辨率处理的谱反演方法。谱反演是近几年发展起来的一种新的提高地震数据分辨率处理技术。其基本原理是通过在频率域选取合理的频率区间,建立计算与实际反射系数序列频谱之间的目标函数,反
地震波在储层岩石传播的过程中,孔隙流体的流动会导致地震波速度的频散和衰减,从而影响储层岩石的地震响应特征。而地震岩石物理建立了储层性质和地震属性之间的定量关系。基于不同流体流动机制建立的理论模型描述了流体性质、饱和度、渗透率等因素对速度频散和衰减的影响。另一方面,频率域的正演方法可以利用频散速度将储层性质与地震响应联系起来。介观尺度的流体流动通常是由于存在部分饱和流体性质差异或岩石骨架性质差异所产
苏里格气田地质储量丰富,属于典型的低孔、低渗致密砂岩储层。这类储层具有孔隙类型多样、气-水关系复杂、储层厚度变化大、含气层与围岩阻抗差小的特点,使得这类储层地球物理特征不明显,为油气识别带来了巨大的挑战。本文结合该区实际的储层及地震参数,采用地震物理模拟和岩石物理相结合的手段,设计连通与不连通孔隙两类储层模型,考虑储层与围岩波阻抗差小的特点,用添加空心微珠和硅胶粒的方式成功地降低了连通储层的速度和
油气勘探开发的不断深入凸显了储层参数精准预测的重要性。地震属性包含丰富的储层地质信息,能够对地下构造、岩性及物性进行综合地反映。因此,地震属性技术是一种有效的储层预测方法。实施地震属性技术时,通常需要提取多个地震属性,在获得大量储层信息的同时也面临着属性冗余等问题。因此,有必要对众多地震属性进行优化,以简化相关计算过程和提高储层预测精度。本文采用线性的主成分分析法与非线性的局部线性嵌入法对地震属性
水力压裂是非常规油气开发和油气田增产的重要技术手段,而微地震监测目前是评价水力压裂效果最主要的方法,其中震源定位是核心环节。本文研究了微地震监测正演模拟技术以及微地震震源定位反演方法,并通过多个正演算例对定位方法进行了测试和分析。正演部分首先研究了交错网格有限差分模拟技术,并详细阐述了PML吸收边界和震源加载方式。由于地表近似自由固体界面,地面监测时接收到的是入射波和反射波叠加的信号,而微地震监测
压制多次波的方法主要可归纳为三大类:基于有效波与多次波可分离性的直接滤波方法、基于波动理论的预测相减方法和基于全波场反演的方法。直接滤波方法多次波压制通常适用于地下构造相对简单的地震资料,而对于地下构造复杂的地区,由于不满足其假设条件,滤波法会失效,需要使用预测相减方法。最典型的预测相减方法为表层相关多次波压制方法SRME(Surface-Related Multiple Elimination)
越来越多的研究者发现,对于不同的学生背景和学科内容,教师在有效教学行为方面存在差异。为了对教师区分性效能进行检验和分析,运用增值性评价方法对北京市M区24名高中语文教师的效能进行研究,结果表明:(1)语文教师在教授积累和阅读两个内容时存在区分性效能;(2)针对不同的学生群体,同一个教师的效能也存在不同,教师职称、评价反馈能力和师生互动三个变量对不同能力的学生群体存在差异性影响,教师职称和评价反馈能
Q值(品质因子)是描述黏弹性介质对地震波吸收程度的表征指标之一,对地下异常体、构造或者含流体区域的反映往往比其他参数的敏感性更高,可用于岩石性质、断层和含流体分布的识别与解释,因此精确估计地下介质的Q值对寻找隐蔽油气藏、预测岩性、研究地壳结构等方面具有重要意义。Q层析成像在天然地震中应用相当广泛,本文将频率域Q值层析成像方法应用于主动源海洋宽角地震数据中,获取衰减模型结果,结合已知的速度信息以解释