一种基于L0半范数的绕射波成像方法及装置

摘要

本发明公开了一种基于L0半范数的绕射波成像方法，获取一个去除反射波的地震数据，作为输入数据；对地下成像空间进行离散化，并选取任意一个成像点作为当前绕射成像点，并依据输入地震数据中的炮点与检波点关系及给定的速度模型，计算上述任意绕射成像点的格林函数；循环执行未选取的其它成像点，并计算其相应的格林函数，直至得出地下成像空间所有成像点格林函数；构建基于L0半范数的绕射波成像模型；由同伦分析迭代算法求解该模型，得出绕射波成像结果。本发明还公开了一种基于L0半范数的绕射波成像装置。通过本发明可以提高地震资料成像分辨率，并增强绕射波成像信噪比，从而更容易识别与储层空间连通性有关的小尺度地质单元。

主权项

一种基于L0半范数的绕射波成像方法，其特征在于，包括以下步骤：步骤101：获取一个去除反射波的地震数据，作为输入数据；步骤102：对地下成像空间进行离散化，并选取任意一个成像点作为当前绕射成像点，并依据输入地震数据中的炮点与检波点关系及给定的速度模型，计算上述任意绕射成像点的格林函数；步骤103：循环执行上述步骤102中地下成像空间中未选取的其它成像点，并计算其相应的格林函数，直至得出地下成像空间所有成像点格林函数；步骤104：根据地下成像空间所有成像点的格林函数与初始成像模型，构建基于L0半范数的绕射波成像模型，设计如下：$\min J_{\alpha }\left(m\right):=\frac{1}{2}\left|\right|Gm-d|{|}_{l_2}^2+\alpha |\left|m\right|{|}_{l_0}$其中，min表示最小化，J_α(m)为目标函数，m为求解的绕射模型，数学符号:＝表示定义为，G为格林函数，d为去除反射波的地震数据，α为正则化因子，表示l₂范数，为l₀半范数；步骤105：由同伦分析迭代算法求解该模型，得出绕射波成像结果。