[发明专利]一种海上自适应压制鬼波的宽频逆时偏移成像方法

摘要

本发明涉及一种海上自适应压制鬼波的宽频逆时偏移成像方法，其特征在于包括以下内容：1)去噪后的叠前炮集地震记录；2)将叠前炮集地震记录进行时间域速度分析，并进行深度域速度建模，得到深度域速度体；对叠前炮集地震记录进行数据解析，拆分为数据和道头两部分信息，并利用道头信息生成炮集数据索引文件；3)读取炮集数据，将读取的炮集数据作为宽频逆时偏移输入数据体，并将深度域速度体作为宽频逆时偏移输入参数，进行三维鬼波自适应压制，完成单炮宽频逆时偏移成像得到单炮集偏移结果；4)对单炮集偏移结果进行噪音压制处理，并根据单炮偏移结果的实际空间位置坐标，对所有噪音压制后的单炮偏移结果进行叠加，获得最终的宽频成像结果。

主权项

1.一种海上自适应压制鬼波的宽频逆时偏移成像方法，其特征在于包括以下内容：1)对地震记录进行叠前噪音压制处理，获得去噪后的叠前炮集地震记录；2)将叠前炮集地震记录进行时间域速度分析，并进行深度域速度建模，得到深度域速度体；对叠前炮集地震记录进行数据解析，拆分为数据和道头两部分信息，并利用道头信息生成炮集数据索引文件；3)根据炮集数据索引文件记录的数据索引关系，并行读取炮集数据，将读取的炮集数据作为宽频逆时偏移输入数据体，并将深度域速度体作为宽频逆时偏移输入参数，进行三维鬼波的自适应压制，完成单炮宽频逆时偏移成像得到单炮集偏移结果，具体过程为：根据地震波传播理论，鬼波表达式为：式中，θ为鬼波初射角，ω为角频率，z为震源或检波点的沉放深度，v是海水速度，f是鬼波陷频点，将含鬼波新边值条件应用于逆时偏移波场延拓所用波动方程中，以实现成像过程中压制鬼波，具体形式如下：a)实际检测点为某个高度，根据逆时偏移震源波场正演所用新的修改边值条件公式(1)，利用时间方向为四阶差分精度的伪谱法计算正演波场PF：式中，P_F为正演波场，为P_F傅里叶到频率域波场，v是海水速度，t是时间，ω为角频率，x，y，z分别为空间坐标，x_s为炮点坐标位置，为频率域的震源点鬼波波场，为频率域震源函数，为空间导数算子；b)实际检波点下沉一段距离，根据逆时偏移检波点波场逆推所用新的修改边值条件公式(2)，利用时间方向为四阶差分精度的伪谱法计算逆推波场PB：式中，P_B为检波点逆推波场，为P_B傅里叶到频率域波场，v是海水速度，t是时间，ω为角频率，x，y，z分别为空间坐标，x_s为炮点坐标位置，为频率域的检波点鬼波波场，为频率域震源函数；c)将正演波场PF和逆推波场PB应用互相关成像条件公式(3)进行成像，获得单炮宽频逆时偏移成像结果：式中，PF为正演波场，PB为检波点逆推波场，R(x)偏移成像振幅，t是时间，x为空间坐标；d)对单炮宽频逆时偏移成像结果采用时间方向为四阶差分精度的伪谱法数值解法进行抗频，具体公式如下：式中，F是傅里叶算子，F‑1是逆傅里叶算子，kx、ky和kz分别是x、y和z三个方向的波数，p为波场值，c为介质速度，Δt为时间间隔，n为时间方向差分片；4)对每个单炮集偏移结果进行噪音压制处理，并根据每个单炮偏移结果的实际空间位置坐标，对所有噪音压制后的单炮偏移结果进行叠加，获得最终的宽频成像结果。