[发明专利]一种起伏地表变偏移距VSP积分叠前深度偏移方法

说明书

本发明提供一种起伏地表变偏移距VSP积分叠前深度偏移方法，包括输入变偏移距VSP共检波点数据、二维地质模型和偏移参数；根据输入数据，建立炮点、检波点在二维地质模型中的坐标；采用两点射线追踪计算散射点到炮点的走时，保存各个散射点的走时表；采用两点射线追踪计算散射点到检波点的走时，保存各个检波点的走时表；实现起伏地表变偏移距VSP共检波点积分叠前深度偏移成像；积分叠前深度偏移重排成共成像道集，共成像道集叠加，得到变偏移VSP积分叠前深度偏移叠加成像。本发明避免了走时表重复计算，算效率高；直接从起伏地表偏移，避免了静校正对波场的影响；深度域成像，有利于井旁构造研究、储层描述，有利于指导油气开发。

技术领域

本发明涉及地球物理勘探中地震数据成像方法，属于垂直地震处理技术范畴，是一种起伏地表变偏移距VSP积分叠前深度偏移方法和装置。

背景技术

近年来，变偏移距VSP(即Walkaway VSP)已实现工业化应用。主要是通过变偏移距VSP成像，研究井旁构造、储层描述。

刘少勇等研究了《起伏地表Kirchhoff积分法叠前深度偏移方法研究与应用》，文献研究了三维地面地震的最短路径走时计算方法，实现了三维地面地震起伏地表Kirchhoff积分法叠前深度偏移。

现有的变偏移距VSP积分成像，多为时间域、平地表的，起伏地表的变偏移距VSP波场需要静校正处理。本发明提出了一种起伏地表变偏移距VSP积分叠前深度偏移方法。无需静校正处理；采用两点射线追踪，计算散射点到起伏地表炮点的走时存为走时表，计算检波点到散射点走时存为走时表；将孔径内的散射点到炮点的走时加上检波点到散射点的走时得到炮检对的共散射走时，将共散射走时对应的波场映射为散射点深度对应的波场，孔径内积分求和得到积分叠前深度偏移成像。

发明内容

本发明目的在于实现起伏地表变偏移距VSP积分叠前深度偏移成像，为深度域井旁构造研究、储层描述提供数据。

具体技术方案为：

一种起伏地表变偏移距VSP积分叠前深度偏移方法，包括以下步骤：

(1)输入变偏移距VSP共检波点数据、二维地质模型和偏移参数；

(1.1)输入变偏移距VSP共检波点波场数据{VSPData}，读取炮点坐标{SX,SY}、炮点深度SZ、检波点坐标{RX,RY}、检波点深度RZ、井口大地坐标{WMX,WMY}；

(1.2)输入二维网格层速度地质模型{ModVel}，模型横坐标MX、模型横纵坐标MZ、井口模型坐标MWMX信息；

(1.3)输入偏移孔径aper、成像起始坐标xmig1、成像结束坐标xmig2、成像网格dxmig、成像终止深度zmig2、成像深度步长dzmig参数。

(2)根据输入数据，建立炮点、检波点在二维地质模型中的坐标；

用步骤(1)中的炮点坐标、炮点深度、检波点坐标、检波点深度、井口大地坐标、井口模型坐标，建立炮点、检波点在二维地质模型中的坐标：

炮点在二维地质模型中的坐标为：

MSZ_i＝SZ_i

其中，MSX_i是第i个炮点在二维地质模型中的坐标，SX_i、SY_i是第i个炮点的大地坐标，SZ_i是第i个炮点的相对于二维地质模型的深度、表示起伏地表，WMX、WMY是井口的大地坐标，MWMX是井口模型坐标，sign是符号函数；

检波点在二维地质模型中的坐标为：