[发明专利]一种大斜井炮检移动VSP自由表面多次波成像方法和装置

说明书

本发明公开了一种大斜井炮检移动VSP自由表面多次波成像方法和装置，所述方法包括以下步骤：S1.输入大斜井炮检移动VSP全波场共检波点道集、2维地质模型、偏移参数；S2.根据输入数据，建立炮点、检波点在2维地质模型中的坐标；S3.对输入的共检波点道集做快速傅里叶变换到频率波数域，设置频率域震源；S4.单检波点的自由表面多次波单程波叠前深度偏移成像；S5.完成所有检波点的炮检移动VSP自由表面多次波叠前深度偏移成像；S6.共成像道集叠加，得到炮检移动VSP自由表面多次波叠前深度偏移叠加成像。本发明有效利用自由表面多次波所携带的丰富信息，实现了炮检移动VSP自由表面多次波叠前深度偏移成像，为地质研究、钻井导向提供了准确的数据基础。

技术领域

本发明涉及地球物理勘探中地震数据偏移成像，特别是涉及一种大斜井炮检移动VSP自由表面多次波成像方法和装置。

背景技术

炮检移动VSP是针对大斜井设计的观测系统，即炮线布设于地表、位于井中检波器的正上方，炮线随着检波器移动。炮检移动VSP的优势是覆盖次更加均匀。这种观测系统主要用于斜井钻进过程中的导向入靶。通常，自由表面多次波作为干扰波被衰减了。与反射相比，自由表面多次波具有易识别、照明宽、反射角小、传播路径长的特点。

吴世萍等研究了《Walkaway VSP多次波成像技术研究》，文献将地震相干成像用于Walkaway VSP多次波成像。李建国等研究了《Walkaway VSP自由表面多次波叠前深度偏移成像》，文献采用单程波延拓实现了变偏移距VSP的自由表面多次波成像；但是由于针对炮检移动VSP的自由表面多次波成像的研究还不成熟，在地址研究、钻进导向等领域还存在着诸多不便。

发明内容

本发明的目的在于克服现有技术的不足，提供一种大斜井炮检移动VSP自由表面多次波成像方法和装置，有效利用自由表面多次波所携带的丰富信息，实现了炮检移动VSP自由表面多次波叠前深度偏移成像，为地质研究、钻井导向提供了准确的数据基础。

本发明的目的是通过以下技术方案来实现的：一种大斜井炮检移动VSP自由表面多次波成像方法，包括以下步骤：

S1.输入大斜井炮检移动VSP全波场共检波点道集、2维地质模型、偏移参数；

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

S3.对输入的共检波点道集做快速傅里叶变换到频率波数域，设置频率域震源；

S4.将频率域震源从检波点深度向上延拓至自由表面，频率波数域炮检移动VSP共检波点全波场、震源波场从自由表面向下延拓，实现单检波点的自由表面多次波单程波叠前深度偏移成像；

S5.循环执行步骤S4，完成所有检波点的炮检移动VSP自由表面多次波叠前深度偏移成像；

S6.将步骤S5的成像重排成共成像道集，共成像道集叠加，得到炮检移动VSP自由表面多次波叠前深度偏移叠加成像。

进一步地，所述步骤S1包括：

S101.输入大斜井炮检移动VSP全波场共检波点道集{VSPData}，读取炮点坐标{SX,SY}、炮点深度SZ、检波点坐标{RX,RY}、检波点深度RZ、井口大地坐标{WMX,WMY}；

S102.输入纵波2维网格层速度地质模型{V}，模型横坐标MX、模型横纵坐标MZ、井口模型坐标MWMX信息；

S103.输入成像算子最大倾角dip、成像终止深度zmig2、成像深度步长dzmig、波场{VSPData}最小频率fl、波场{VSPData}最大频率参数fh。

所述步骤S2包括：

S201.利用步骤S1中的炮点坐标、炮点深度、检波点坐标、检波点深度、井口大地坐标、井口模型坐标，建立炮点、检波点在2维地质模型中的坐标：