复杂构造油气藏成为当前勘探的重点，也是进行高分辨率勘探瓶颈。为了获得精确的复杂构造成像，采用具有高角度、频散小、适应纵横向速度变化的傅里叶有限差分法实现叠前深度偏移，但该方法对陡倾角成像存在明显的误差。本发明专利技术利用全局优化的非线性多元最小二乘法对傅里叶有限差分算子中的系数进行优化，从而在保持计算效率的前提下，提高了成像角度及陡倾角的成像精度。
【技术实现步骤摘要】
本专利技术涉及地震资料处理领域，特别是涉及傅里叶有限差分法叠前深度偏移。
技术介绍
当前复杂构造或速度剧烈变化的地区成为油气勘探的重点，而对于该地区采用常规的偏移处理将得不到精确的地下构造成像。为此在处理中采用叠前深度偏移技术，叠前深度偏移能获得复杂构造及纵横向速度剧烈变化地区的精确成像，是复杂构造地区精确成像有效方法。目前常用的叠前深度偏移主要是基于两种方法kirChh0fT积分法和波动方程法。这两种方法中，波动方程叠前深度偏移是复杂构造成像最有效的手段。它没有像 kirchhoff积分法一样对方程做高频近似，而是用描述波在复杂介质中的传播算子作波场外推算子，故其更适合复杂构造的精确成像，同时还具有保幅的特性。在波动叠前深度偏移中，其核心是波场外推算子，常用的波场外推算子有频率一波数域相移法、空间一频率域有限差分法、分步Rnirier法、傅里叶有限差分法及广义屏法。频率一波数域相移法具有 90°倾角成像、频散性能好和计算效率高的优点，但它不能适应强横向速度变化；空间一频率域有限差分法能适应任意横向速度变化，但它存在着频散和陡倾角成像的问题；分步 Fourier法是在频率一空间域和频率一波数域实现波场递推，可反映速度变化的一级近似， 计算效率高，但其只适应弱的横向速度变化的介质；广义屏法是在频率一波数域和频率一空间域交替延拓波场，其实质是参考波场的内插，通过加权自动调节并计算实际波场，计算效率高，是较理想的二、三维算法，但该方法存在着陡倾角成像精度低的问题；傅里叶有限差分法是在分步Rmrier基础上通过引入频率一空间域的剩余有限差分项，对二阶速度变化进行补偿，提高对陡倾角及强横向速度变化的适应性，该方法具有相移法高角度、频散小及有限差分法适应强横向速度变化的特点，但该方法存在着计算效率及在三维中应用困难的问题，同时对陡倾角的精确成像存在着一定的误差。为了减小傅里叶有限差分陡倾角成像的误差，可通过增加有限差分校正项的阶数，但该做法会增加计算量，降低计算效率，为了在保持计算效率的同时提高算子的精度， 可对其傅里叶有限差分的校正项的系数进行优化。本专利技术采用非线性多元最小二乘法进行参数优化，在优化中考虑了速度的对比度变化、角度及频率变化等参量的影响，而常规的优化系数只考虑了速度的对比度及角度的影响，故其精度不足。采用本文的优化系数方法可明显提高其精度，从而提高对陡倾角的成像能力，使得傅里叶有限差分在保持计算效率的同时提高了陡倾角的成像精度。
技术实现思路
为更有效地克服傅里叶有限差分法叠前深度偏移中存在的上述缺陷，本专利技术目的是通过对傅里叶有限差分的校正项的系数进行优化，从而提高对陡倾角的成像能力。本专利技术的目的是提供一种在保持傅里叶有限差分法计算效率的同时提高陡倾角的成像精度。该方法步骤如下首先，在每个深度上，频率一空间域的上行波方程可分解为两部分，即与背景速度有关的背景波场部分及与变化速度有关的扰动波场部分，即权利要求1.，该方法包括步骤1对扰动项的差异算子进行非线性多元最小二乘法逼近求解优化系数，在求解优化系数的过程中，考虑了速度对比度、角度、频率对系数的影响，利用非线性多元最小二乘法进行求解；步骤2利用各自系数求偏导数，求解的系数，该系数即为优化后的傅里叶有限差分算子的系数；步骤3读入SEG/EGEA模型数据，对数据进行傅里叶分析；步骤4对频率范围内的每一个波场，采用非线性多元最小二乘法优化系数傅里叶有限差分算子对震源波场向下延拓；步骤5延拓过程中下一个深度的波场为上一个深度波场延拓后的结果，对延拓后的结果进行叠加成像，输出成像结果。2.根据权利要求1所述的方法，其特征在于利用非线性多元最小二乘法优化扰动项差异算子的系数；在优化系数中考虑了速度对比度P、角度β、频〗ω对优化系数的影响，使得优化后的系数更符合波在地下介质的传播特征。全文摘要复杂构造油气藏成为当前勘探的重点，也是进行高分辨率勘探瓶颈。为了获得精确的复杂构造成像，采用具有高角度、频散小、适应纵横向速度变化的傅里叶有限差分法实现叠前深度偏移，但该方法对陡倾角成像存在明显的误差。本专利技术利用全局优化的非线性多元最小二乘法对傅里叶有限差分算子中的系数进行优化，从而在保持计算效率的前提下，提高了成像角度及陡倾角的成像精度。文档编号G01V1/36GK102323614SQ20111014577公开日2012年1月18日 申请日期2011年6月1日 优先权日2011年6月1日专利技术者罗仁泽, 黄元溢 申请人:西南石油大学本文档来自技高网...
【技术保护点】
１．一种基于最小二乘法优化系数的傅里叶有限差分偏移方法，该方法包括：步骤１ 对扰动项的差异算子进行非线性多元最小二乘法逼近求解优化系数，在求解优化系数的过程中，考虑了速度对比度、角度、频率对系数的影响，利用非线性多元最小二乘法进行求解；步骤２ 利用各自系数求偏导数，求解的系数，该系数即为优化后的傅里叶有限差分算子的系数；步骤３ 读入ＳＥＧ／ＥＧＥＡ模型数据，对数据进行傅里叶分析；步骤４ 对频率范围内的每一个波场，采用非线性多元最小二乘法优化系数傅里叶有限差分算子对震源波场向下延拓；步骤５ 延拓过程中下一个深度的波场为上一个深度波场延拓后的结果，对延拓后的结果进行叠加成像，输出成像结果。
【技术特征摘要】
【专利技术属性】
技术研发人员：罗仁泽，黄元溢，
申请(专利权)人：西南石油大学，
类型：发明
国别省市：90
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。