[发明专利]考虑衰减频散的深度域地震记录合成方法

说明书

本发明实现了一种考虑衰减和频散的深度域地震记录合成方法。本发明提供的方法首先基于深度域常Q或变Q模型构建衰减频散矩阵，然后利用反射系数、深度域地震子波矩阵的傅里叶变换结果以及衰减频散矩阵直接在深度域合成包含衰减和频散效应的地震记录。

技术领域

本发明属于石油地震勘探领域，具体是涉及一种考虑衰减和频散的深度域地震记录合成方法。

背景技术

地震波在地下介质的传播过程中，其传播特征既受地层速度的影响，还会受到地下介质固有的衰减和频散效应的影响。尤其是对于深度域地震波场数据，即使在不考虑地层衰减和频散效应的情况下，地震波的波形是依赖层速度而变化的，层速度越大，地震反射波形越宽，层速度减小，则地震反射波形变窄；当地震反射界面上下的地层速度不同时，该界面处的深度域地震波形不对称，且层速度越高，波形拉伸畸变就越明显，深度域地震波的这种现象，与时间域中地震波形因地层固有的衰减和频散效应所致的能量衰减和相位畸变的现象相似，而时间域地震记录的波形是不依赖层速度的。目前，已有较多的方法可实现对包含衰减和频散效应的时间域地震记录的合成。由于地下介质的固有衰减和频散效应，它对深度域地震波场传播特征的影响也是非常值得关注和研究的问题，深度域地震波的波形特征同时受到地层速度及其衰减频散的共同影响，这与时间域地震波场表现的特性显著不同，开发包含衰减和频散效应的深度域地震记录合成方法，有助于明确地下介质的频散衰减在深度域地震波场上的表现特征及规律，对于从实际的深度域地震资料中提取用于度量地层衰减频散的品质因子Q值，提高储层预测的精度、实现地震数据的高分辨率处理等都具有重要意义。

目前，合成深度域地震记录的方法有5种：①在时域用褶积模型合成地震记录，然后时深转换到深度域；②将地层速度视为匀速，因匀速介质中子波保持不变，可直接应用褶积模型制作深度域合成记录；③变换域速度函数制作深度域合成地震记录，此方法需先建立转换函数，通过变换使速度在“新深度域”中为常数、“子波”具有时不变性，然后应用褶积模型制作合成记录；④建立深度域的速度模型，波动方程正演获得时域合成记录，再通过深度域偏移得到深度域。⑤基于时间域地震子波和深度域速度信息，在各深度域位置构建点扩散函数，与反射系数褶积后获得深度域合成地震记录。但是，这几种方法在合成地震记录时均未考虑衰减和频散效应，用这些方法均无法在深度域地震记录中体现地下介质固有的衰减和频散的影响。

发明内容

本发明的目的在于克服现有技术的不足，提供了一种考虑衰减和频散的深度域地震记录合成方法，所述方法包括以下主要步骤：

⑴构建包含衰减和频散效应的衰减频散矩阵X，矩阵X中元素的表达式如下：

式中，χ_m,n是矩阵X中索引号为(m,n)时的元素，d_m是索引号为m时的深度，k_n是索引号为n时的波数，Q_m是索引号为m时的品质因子，k_r为参考波数，exp[·]表示指数运算，i为虚数，N为深度域数据包含的采样点数量；

⑵基于反射系数和矩阵X按照下式构建矩阵A：

A＝r⊙X

式中，r为反射系数行向量，⊙表示基本积运算；

⑶基于深度域地震子波矩阵的傅里叶变换结果Φ和矩阵A，构建矩阵B，矩阵B中元素的表达式如下：

β_m,n＝exp[id_mk_n](α_m,nφ_m,n)

式中，α_m,n，β_m,n和φ_m,n分别是矩阵A，B和Φ中索引号为(m,n)时的元素；

⑷基于矩阵B，按照下式计算包含衰减和频散效应的深度域合成地震记录y：