利用波动方程求解地震波阻抗的方法和装置制造方法及图纸

本申请公开了利用波动方程求解地震波阻抗的方法和装置。该方法包括：根据波动方程建立目标函数φ(m)，其中，m为由傅里叶级数展开形式表示的波阻抗模型的傅里叶系数构成的矩阵；基于所述目标函数确定矩阵m中各个元素的值；将矩阵m中各个元素的值代入所述波阻抗模型。根据本申请，可以实现超深层少井或未钻进至目标层以及强非均质性储层高精度宽带波阻抗反演及应用，使得有井段与常规反演效果在波动方程反演的基础上对于局部边缘特征刻画分辨率更高，使得无井段在保持局部特征的同时提高深层的全局连续性，可为特深层断溶型储层的刻画及早期储量评估提供数据支撑。刻画及早期储量评估提供数据支撑。刻画及早期储量评估提供数据支撑。

【技术实现步骤摘要】
利用波动方程求解地震波阻抗的方法和装置


[0001]本专利技术属于地震勘探领域，更具体地，涉及一种利用波动方程求解地震波阻抗的方法和一种利用波动方程求解地震波阻抗的装置。

技术介绍

[0002]波阻抗信息作为联系地质和地球物理的一座桥梁，在地震解释中有重要作用。地震反演得到波阻抗，无论是叠前波阻抗反演还是叠后波阻抗反演，反演方法一般以基于褶积模型的方法为主。常用算法大体可分为以反褶积法为基础的直接反演(如递归反演，道积分等)、以模型为基础的反演和约束反演等，但褶积模型在平层假设的前提下，对于横向分辨率的提高存在局限性；其次常规方法缺乏超深层无井及少井情况下储层特征描述能力。对于波形反演类的方法在反演弹性参数时的主要问题在于二维及以上波动方程受多参数解耦及反演策略的限制，同时反演多个参数才求解多参数的稳定性和计算效率收到限制，函数的非凸性容易使求解过程陷入局部极小。
[0003]现有反演波阻抗的方法多是基于褶积模型，反演结果依赖于钻井记录构建的初始波阻抗模型。在数据采集不完整、少井或者钻井深度不到目标层的地区初始波阻抗模型不易构建，难以得本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种利用波动方程求解地震波阻抗的方法，其特征在于，所述方法包括：根据波动方程建立目标函数φ(m)：其中，m为由傅里叶级数展开形式表示的波阻抗模型的傅里叶系数构成的矩阵，是坐标(pΔx,qΔy)处零偏移距的观测地震数据，是所述波阻抗模型上坐标(pΔx,qΔy)处的合成地震数据，Δx、Δy和Δt分别表示在三个坐标维度下的采样步长，P、Q和H对应表示在三个坐标维度下的最大采样点数；基于所述目标函数确定矩阵m中各个元素的值；将矩阵m中各个元素的值代入所述波阻抗模型。2.根据权利要求1所述的方法，其特征在于，由傅里叶级数展开形式表示的波阻抗模型如下所示：其中，I(x,y,t)是波阻抗参数，{x,y,t}分别是在三个坐标维度下的取值，{L,M,N}分别是三个坐标维度上谐波项的最大数，{Δk,Δω}是空间域的基本波数，(a
lmn
,b
lmn
,...,h
lmn
)构成矩阵m。3.根据权利要求1所述的方法，其特征在于，基于所述目标函数确定矩阵m中各个元素的值，包括：S21，计算目标函数φ(m)对每个傅里叶系数的一阶偏导coef
i
为矩阵m中第i个元素，i＝1,2,...,R，R为矩阵m中的元素个数；S22，基于下式更新矩阵m
k
，m
k
为开始第k次迭代时矩阵m的当前取值：m
k+1
＝m
k
+αΔm
k
，其中，Δm
k
为负梯度矩阵，Δm
k,i
为负梯度矩阵Δm
k
中的第i个元素，α为步长，m1为已给出的初始傅里叶系数矩阵；S23，判断是否满足迭代终止条件，如果满足，则令矩阵m＝m
k+1
，并结束迭代；如果不满足，则令k＝k+1，并回到步骤S21，进入下次迭代。4.根据权利要求3所述的方法，其特征在于，根据下式计算一阶偏导
其中，[Δu
pq
]是所述观测地震数据与所述合成地震数据的数据残差，是波场对傅里叶系数coef
i
的Frechet导数。5.根据权利要求3所述的方法，其特征在于，所述迭代终止条件包括：目标函数的值小于预设阈值；或者迭代次数达到最大迭代次数。6.一种利用波动方程求解地震波阻抗的装置，其特征在于，所述装置包括：目标函数建立单...

【专利技术属性】
技术研发人员：朱童，张如一，李京南，
申请(专利权)人：中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：