[发明专利]基于被动源地震波形及其逆时成像的联合速度反演方法

摘要

本发明涉及一种基于被动源地震波形及其逆时成像的联合速度反演方法，属于被动源地震探测的速度建模领域。主要包括如下步骤：1)对地震数据预处理并提取震源子波；2)计算逆时延拓波场，进行逆时成像和震源定位；3)计算理论地震记录与观测地震记录之间的差；4)建立反演速度的目标函数；5)基于伴随状态法计算目标函数对模型参数的梯度以及更新方向和步长；6)更新纵、横波速度模型；7)判断是否满足迭代终止条件，如果满足则输出定位结果及速度模型，否则返回步骤2。本发明的联合速度建模方法在震源位置未知情况下，可以获得高分辨率速度模型。

主权项

1.一种基于被动源地震波形及其逆时成像的联合速度反演方法，其特征在于，包括如下步骤：/n步骤一：在地面观测被动源地震三分量数据；/n步骤二：对观测的被动源地震数据进行去噪、纵横波分离等预处理工作；/n步骤三：提取震源子波；/n步骤四：计算目标函数中的权重系数；/n步骤五：计算纵波和横波的逆时延拓波场；/n步骤六：通过逆时成像条件进行震源成像和定位；/n步骤七：构建各震源成像结果对应的惩罚函数；/n步骤八：由逆时成像获得的震源位置，计算当前模型下的正传波场及其对应的地震记录；/n步骤九：计算观测地震记录和计算地震记录之间的差值；/n步骤十：利用纵波自相关、横波自相关、纵波和横波互相关成像条件的成像结果差异和观测与计算的地震记录差值计算联合的目标函数；/n步骤十一：判断是否满足迭代终止条件，如果满足终止条件，则结束反演流程，输出计算的速度模型和震源定位结果，否则，继续进行如下步骤；/n步骤十二：根据目标函数计算伴随震源函数；/n步骤十三：由伴随震源函数计算伴随波场；/n步骤十四：基于伴随状态法计算目标函数对模型速度参数的梯度；/n步骤十五：计算目标函数对模型速度参数梯度的共轭方向；/n步骤十六：使用线性搜索法计算更新步长；/n步骤十七：计算模型修正量并更新速度模型；/n步骤十八：将更新后的速度模型作为新的初始速度模型，返回步骤四。/n