CN111751875A - 一种变偏移距vsp带限角度积分叠前时间偏移方法和装置 - Google Patents

Abstract

本发明公开了一种变偏移距VSP带限角度积分叠前时间偏移方法和装置，所述方法包括：包括以下步骤：S1.输入2维深度速度模型、变偏移距VSP共检波点道集、变偏移距VSP共炮点走时表、变偏移距VSP共检波点走时表、偏移参数；S2.建立炮点、检波点在2维地质模型中的坐标；S3.设计带限角度滤波器的带限角度加权系数；S4.实现变偏移距VSP共检波点带限角度积分叠前时间偏移成像；S5.将得到的带限角度积分叠前时间偏移重排成共成像道集，共成像道集叠加，得到变偏移VSP带限角度积分叠前时间偏移叠加成像。本发明能够直接生成变偏移距VSP带限角度的积分叠前时间偏移成像，无需生成叠前角度道集，有效提高了数据处理效率。

Description

一种变偏移距VSP带限角度积分叠前时间偏移方法和装置

技术领域

本发明涉及地震数据偏移成像方法，特别是涉及一种变偏移距VSP带限角度积分叠前时间偏移方法和装置。

背景技术

变偏移距VSP(Walkaway VSP)是指多个检波器组成的阵列陈放在观测井中，地表布设一条人工激发的多个炮点组成的炮线，是一种VSP的二维观测方式，在炮点人工激发地震数据的条件下，实现变偏移距VSP井中数据采集，采集系统立体图如图1所示；变偏移距VSP采集，其实就是在每一个炮点激发时，由各个检波点的检波器进行数据采集；变偏移距VSP可以观测到反射纵波、反射转换横波，利用偏移成像研究井旁构造、储层预测。

王楠等研究了《表驱Kirchhoff叠前时间偏移角度域成像方法》，文献利用两点射线追踪计算走时和反射张角，实现地面地震数据Kirchhoff叠前时间偏移角度域成像。王珺等研究了《克希霍夫法VSP多波联合成像》，文献实现了偏移距VSP反射P波、反射S波、透射P波、透射S波克希霍夫法联合成像。

现有的积分叠前时间偏移方法多数是针对地面地震观测方式的，变偏移距VSP的积分叠前时间偏移文献很少，数据处理较为不便。

发明内容

本发明的目的在于克服现有技术的不足，提供一种变偏移距VSP带限角度积分叠前时间偏移方法和装置，能够直接生成变偏移距VSP带限角度的积分叠前时间偏移成像，无需生成叠前角度道集，有效提高了数据处理效率。

本发明的目的是通过以下技术方案来实现的：一种变偏移距VSP带限角度积分叠前时间偏移方法，包括以下步骤：

S1.输入2维深度速度模型、变偏移距VSP共检波点道集、变偏移距VSP共炮点走时表、变偏移距VSP共检波点走时表、偏移参数；

S2.从输入的数据中获取炮点坐标、炮点深度、检波点坐标、检波点深度、井口大地坐标、井口模型坐标，建立炮点、检波点在2维地质模型中的坐标；

S3.在步骤S2的2维地质模型坐标系中，采用射线追踪计算炮点、检波点在共散射点处的射线角度，并从步骤S1输入的数据中获取角度下限、角度上限，设计带限角度滤波器的带限角度加权系数；

S4.在步骤S2的2维地质模型坐标系中，利用步骤S1的模型计算时深表，利用步骤S1中的走时表计算共散射走时，用步骤S3得到的带限角度加权系数孔径内积分求和，实现变偏移距VSP共检波点带限角度积分叠前时间偏移成像；

S5.将步骤S4得到的带限角度积分叠前时间偏移重排成共成像道集，共成像道集叠加，得到变偏移VSP带限角度积分叠前时间偏移叠加成像。

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

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

S102.输入2维网格深度速度模型{ModVel}，读取网格节点数MM×MN、网格大小dg、井口模型坐标MWMX信息；

S103.输入变偏移距VSP共炮点走时表{RayTS}；

S104.输入变偏移距VSP共检波点走时表{RayTR}；

S105.输入偏移孔径aper、成像起始坐标xmig1、成像结束坐标xmig2、成像网格dxmig参数、角度下限ang1、角度上限ang2；其中，偏移孔径是参与积分求和的范围，如1000m、2000m。角度是射线在散射点处的角度，上下限，如10°到20°成像，下限是10，上限是20。

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

S201.建立炮点在2维地质模型中的坐标：

MSZ_i＝SZ_i

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

S202.建议检波点在2维地质模型中的坐标：

MRZ_j＝RZ_j

其中，MRX_j是第j个检波点在2维地质模型中坐标，RX_j、RY_j是第j个检波点的大地坐标，RZ_j是第j个检波点的深度，WMX、WMY是井口的大地坐标，MWMX是井口模型坐标，sign是符号函数。

进一步地，所述步骤S3中，设计的带限角度加权系数如下：

WAngBF_ktr,(i,j)＝cos(θ_ktr,i)·cos(θ_ktr,j) ang1≤θ_ktr,i≤ang2且ang1≤θ_ktr,j≤ang2

其中，θ_ktr,i是第ktr道第i个炮点的各个散射点的射线角度，θ_ktr,j是第ktr道第j个检波点的各个散射点的射线角度，ang1是角度下限、ang2是角度上限，WAngBF_ktr,(i,j)是第ktr道第i个炮点第j个检波点的带限角度加权系数。

射线追踪计算炮点在共散射点处的射线角度的方式如下：

其中，MSX_i、MSZ_i是第i个炮点的模型坐标，x、z是散射点的坐标；

射线追踪计算检波点在共散射点处的射线角度的方式如下：

其中，MRX_i、MRZ_i是第j个检波点的模型坐标，x、z是散射点的坐标。

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

S401.计算步骤S1中2维深度速度模型的时深表：

ModT_k,l＝ModT_k-1,l+2·dg/ModVel_k-1,l

ModH_k,l＝ModH_k-1,l+dg

其中，ModT是模型双程垂直时，ModH是模型深度，k∈[1,MM]是模型行数，l∈[1,MN]是模型列数；

S402.用步骤S1输入的偏移参数定义成像坐标，进行如下处理：

xmig_ktr＝xmig1+(ktr-1)·dxmig ktr＝1,2,…,(xmig2-xmig1)/dxmig+1

其中，xmig_ktr是第ktr道的成像坐标，xmig1是成像起始坐标，xmig2是成像结束坐标；

MigData_ktr,j＝{0}

其中，MigData_ktr,j第j个检波点第ktr道的成像数据，先置为零；

S403.利用步骤S1输入的数据中获取偏移孔径、结合步骤S2中得到的的炮点模型坐标、步骤S402的成像坐标确定成像输入地震道范围：

xmig_ktr-aper≤MSX_kaper,j≤xmig_ktr+aper

其中，xmig_ktr是第ktr道的成像坐标，aper是偏移孔径，MSX_kaper,j是第j个检波点第ktr个成像道孔径内第kaper个炮点的模型坐标。

S404.选择第kaper个炮点xmig_ktr处的走时加上第个j检波点xmig_ktr处的走时，得到xmig_ktr处kaper、j炮检对的共散射走时TSR；

其中，

是第kaper个炮点xmig_ktr处的走时，

第j个检波点xmig_ktr处的走时；

S405.选择第j个检波点第kaper个炮点的波场数据VSPData_kaper,j及带限角度加权系数，把TSR对应的振幅带限角度加权映射为时深表对应的双程垂直，与xmig_ktr处的成像道MigData_ktr,j积分求和；

MigData_ktr,j＝MigData_ktr,j+VSPData_kaper,j·WAngBF_{ktr,(kaper,j)}

VSPData_kaper,j是第j个检波点第kaper个炮点的波场数据，WAngBF_{ktr,(kaper,j)}是第j个检波点第kaper个炮点的xmig_ktr处的带限角度加权系数；

S406.循环执行步骤S403～步骤S405，实现第j个检波点第ktr个成像道孔径内全部地震道的带限角度积分叠前时间偏移成像，得到MigData_ktr,j；

S407.循环执行步骤S402～S406，实现第j个检波点全部成像道的带限角度积分叠前时间偏移成像得到MigData_j；

S408.并行循环执行步骤S402～S407，实现全部检波点全部成像道的带限角度积分叠前时间偏移成像得到MigData。

一种变偏移距VSP带限角度积分叠前时间偏移装置，包括：

数据输入模块，用于输入2维深度速度模型、变偏移距VSP共检波点道集、变偏移距VSP共炮点走时表、变偏移距VSP共检波点走时表、偏移参数；

坐标建立模块，用于从输入的数据中获取炮点坐标、炮点深度、检波点坐标、检波点深度、井口大地坐标、井口模型坐标，建立炮点、检波点在2维地质模型中的坐标；

带限角度加权系数设计模块，用于在2维地质模型坐标系中，采用射线追踪计算炮点、检波点在共散射点处的射线角度，并从输入的数据中获取角度下限、角度上限，设计带限角度滤波器的带限角度加权系数；

积分叠前时间偏移成像模块，用于在2维地质模型中，计算时深表和共散射走时，利用得到的带限角度加权系数孔径内积分求和，实现变偏移距VSP共检波点带限角度积分叠前时间偏移成像；

共成像道集叠加模块，用于将得到的带限角度积分叠前时间偏移重排成共成像道集，共成像道集叠加，得到变偏移VSP带限角度积分叠前时间偏移叠加成像。

本发明的有益效果是：本发明通过射线追踪设计了一个带限角度滤波器；调用炮检点走时表得到共散射点走时；将共散射点走时对应的地震数据映射为双程垂直时，偏移孔径内带限角度加权并积分求和，实现共检波点波场积分叠前时间偏移成像；多个检波点的共散射成像叠加，实现变偏移距VSP积分叠前时间偏移成像叠加，有效提高了数据处理效率。

附图说明

图1为变偏移距VSP采集系统的立体图；

图2为本发明的方法流程图；

图3为实施例中输入的2维深度速度模型示意图；

图4为实施例中输入的变偏移距VSP共检波点道集示意图；

图5为实施例中输入的变偏移距VSP共炮点走时表示意图；

图6为实施例中输入的变偏移距VSP共检波点走时表示意图；

图7为实施例中同一检波点、不同散射点射线追踪的带限角度加权系数示意图；

图8为实施例中变偏移距VSP带限角度积分叠前时间偏移成像示意图；

图9为实施例中变偏移距VSP带限角度积分叠前时间偏移成像叠加成像示意图；

图10为本发明的系统原理框图。

具体实施方式

下面结合附图进一步详细描述本发明的技术方案，但本发明的保护范围不局限于以下所述。

如图2所示，一种变偏移距VSP带限角度积分叠前时间偏移方法，包括以下步骤：

S1.输入2维深度速度模型、变偏移距VSP共检波点道集、变偏移距VSP共炮点走时表、变偏移距VSP共检波点走时表、偏移参数；

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

S102.输入2维网格深度速度模型{ModVel}，读取网格节点数MM×MN、网格大小dg、井口模型坐标MWMX信息；

S103.输入变偏移距VSP共炮点走时表{RayTS}；

S104.输入变偏移距VSP共检波点走时表{RayTR}；

S105.输入偏移孔径aper、成像起始坐标xmig1、成像结束坐标xmig2、成像网格dxmig参数、角度下限ang1、角度上限ang2。

在本申请的实施例中，输入的2维深度速度模型如图3所示；输入的变偏移距VSP共检波点道集如图4所示；输入的变偏移距VSP共炮点走时表如图5所示；输入的变偏移距VSP共检波点走时表如图6所示；图3中，横坐标为长度(单位：米)，纵坐标为深度(单位：米)；图4中，横坐标为道号，纵坐标为时间(单位：毫秒)；图5～图6中，横坐标为序号(单位：无)，纵坐标为深度(单位：米)；

S2.从输入的数据中获取炮点坐标、炮点深度、检波点坐标、检波点深度、井口大地坐标、井口模型坐标，建立炮点、检波点在2维地质模型中的坐标；

S201.建立炮点在2维地质模型中的坐标：

MSZ_i＝SZ_i

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

S202.建议检波点在2维地质模型中的坐标：

MRZ_j＝RZ_j

其中，MRX_j是第j个检波点在2维地质模型中坐标，RX_j、RY_j是第j个检波点的大地坐标，RZ_j是第j个检波点的深度，WMX、WMY是井口的大地坐标，MWMX是井口模型坐标，sign是符号函数。

S3.在步骤S2的2维地质模型坐标系中，采用射线追踪计算炮点、检波点在共散射点处的射线角度，并从步骤S1输入的数据中获取角度下限、角度上限，设计带限角度滤波器的带限角度加权系数：

WAngBF_ktr,(i,j)＝cos(θ_ktr,i)·cos(θ_ktr,j) ang1≤θ_ktr,i≤ang2且ang1≤θ_ktr,j≤ang2

其中，θ_ktr,i是第ktr道第i个炮点的各个散射点的射线角度，θ_ktr,j是第ktr道第j个检波点的各个散射点的射线角度，ang1是角度下限、ang2是角度上限，WAngBF_ktr,(i,j)是第ktr道第i个炮点第j个检波点的带限角度加权系数。

射线追踪计算炮点在共散射点处的射线角度的方式如下：

其中，MSX_i、MSZ_i是第i个炮点的模型坐标，x、z是散射点的坐标；

射线追踪计算检波点在共散射点处的射线角度的方式如下：

其中，MRX_i、MRZ_i是第j个检波点的模型坐标，x、z是散射点的坐标。

在本申请的实施例中，对于同一检波点、不同散射点，射线追踪的带限角度加权系数如图7所示，图7中横坐标为序号，纵坐标为时间(单位：毫秒)。

S4.在步骤S2的2维地质模型坐标系中，利用步骤S1的模型计算时深表，利用步骤S1中的走时表计算共散射走时，用步骤S3得到的带限角度加权系数孔径内积分求和，实现变偏移距VSP共检波点带限角度积分叠前时间偏移成像：

S401.计算步骤S1中2维深度速度模型的时深表：

ModT_k,l＝ModT_k-1,l+2·dg/ModVel_k-1,l

ModH_k,l＝ModH_k-1,l+dg

其中，ModT是模型双程垂直时，ModH是模型深度，k∈[1,MM]是模型行数，l∈[1,MN]是模型列数；

S402.用步骤S1输入的偏移参数定义成像坐标，进行如下处理：

xmig_ktr＝xmig1+(ktr-1)·dxmig ktr＝1,2,…,(xmig2-xmig1)/dxmig+1

其中，xmig_ktr是第ktr道的成像坐标，xmig1是成像起始坐标，xmig2是成像结束坐标；

MigData_ktr,j＝{0}

其中，MigData_ktr,j第j个检波点第ktr道的成像数据，先置为零；

S403.利用步骤S1输入的数据中获取偏移孔径、结合步骤S2中得到的的炮点模型坐标、步骤S402的成像坐标确定成像输入地震道范围：

xmig_ktr-aper≤MSX_kaper,j≤xmig_ktr+aper

其中，xmig_ktr是第ktr道的成像坐标，aper是偏移孔径，MSX_kaper,j是第j个检波点第ktr个成像道孔径内第kaper个炮点的模型坐标。

S404.选择第kaper个炮点xmig_ktr处的走时加上第个j检波点xmig_ktr处的走时，得到xmig_ktr处kaper、j炮检对的共散射走时TSR；

其中，

是第kaper个炮点xmig_ktr处的走时，

第j个检波点xmig_ktr处的走时；

S405.选择第j个检波点第kaper个炮点的波场数据VSPData_kaper,j及带限角度加权系数，把TSR对应的振幅带限角度加权映射为时深表对应的双程垂直，与xmig_ktr处的成像道MigData_ktr,j积分求和；

MigData_ktr,j＝MigData_ktr,j+VSPData_kaper,j·WAngBF_{ktr,(kaper,j)}

VSPData_kaper,j是第j个检波点第kaper个炮点的波场数据，WAngBF_{ktr,(kaper,j)}是第j个检波点第kaper个炮点的xmig_ktr处的带限角度加权系数；

S406.循环执行步骤S403～步骤S405，实现第j个检波点第ktr个成像道孔径内全部地震道的带限角度积分叠前时间偏移成像，得到MigData_ktr,j；

S407.循环执行步骤S402～S406，实现第j个检波点全部成像道的带限角度积分叠前时间偏移成像得到MigData_j；

S408.并行循环执行步骤S402～S407，实现全部检波点全部成像道的带限角度积分叠前时间偏移成像得到MigData。

在本申请的实施例中，变偏移距VSP带限角度积分叠前时间偏移成像结果(角度10-25°，与图7检波点对应)如图8所示，图8中，横坐标为道号；纵坐标为时间(单位：毫秒)。

S5.将步骤S4得到的带限角度积分叠前时间偏移重排成共成像道集，共成像道集叠加，得到变偏移VSP带限角度积分叠前时间偏移叠加成像。在本申请的实施例中，变偏移距VSP带限角度积分叠前时间偏移成像叠加成像结果(角度10-25°)如图9所示，图9中，横坐标为道号，纵坐标为时间(单位：毫秒)。

如图10所示，一种变偏移距VSP带限角度积分叠前时间偏移装置，包括：

数据输入模块，用于输入2维深度速度模型、变偏移距VSP共检波点道集、变偏移距VSP共炮点走时表、变偏移距VSP共检波点走时表、偏移参数；

坐标建立模块，用于从输入的数据中获取炮点坐标、炮点深度、检波点坐标、检波点深度、井口大地坐标、井口模型坐标，建立炮点、检波点在2维地质模型中的坐标；

带限角度加权系数设计模块，用于在2维地质模型坐标系中，采用射线追踪计算炮点、检波点在共散射点处的射线角度，并从输入的数据中获取角度下限、角度上限，设计带限角度滤波器的带限角度加权系数；

积分叠前时间偏移成像模块，用于在2维地质模型中，计算时深表和共散射走时，利用得到的带限角度加权系数孔径内积分求和，实现变偏移距VSP共检波点带限角度积分叠前时间偏移成像；

共成像道集叠加模块，用于将得到的带限角度积分叠前时间偏移重排成共成像道集，共成像道集叠加，得到变偏移VSP带限角度积分叠前时间偏移叠加成像。

以上所述是本发明的优选实施方式，应当理解本发明并非局限于本文所披露的形式，不应该看作是对其他实施例的排除，而可用于其他组合、修改和环境，并能够在本文所述构想范围内，通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围，则都应在本发明所附权利要求的保护范围内。

Claims (6)

1.一种变偏移距VSP带限角度积分叠前时间偏移方法，其特征在于：包括以下步骤：

S1.输入2维深度速度模型、变偏移距VSP共检波点道集、变偏移距VSP共炮点走时表、变偏移距VSP共检波点走时表、偏移参数；

S2.从输入的数据中获取炮点坐标、炮点深度、检波点坐标、检波点深度、井口大地坐标、井口模型坐标，建立炮点、检波点在2维地质模型中的坐标；

S3.在步骤S2的2维地质模型坐标系中，采用射线追踪计算炮点、检波点在共散射点处的射线角度，并从步骤S1输入的数据中获取角度下限、角度上限，设计带限角度滤波器的带限角度加权系数；

S4.在步骤S2的2维地质模型坐标系中，利用步骤S1的模型计算时深表，利用步骤S1中的走时表计算共散射走时，用步骤S3得到的带限角度加权系数孔径内积分求和，实现变偏移距VSP共检波点带限角度积分叠前时间偏移成像；

S5.将步骤S4得到的带限角度积分叠前时间偏移重排成共成像道集，共成像道集叠加，得到变偏移VSP带限角度积分叠前时间偏移叠加成像。

2.根据权利要求1所述的一种变偏移距VSP带限角度积分叠前时间偏移方法，其特征在于：所述步骤S1包括：

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

S102.输入2维网格深度速度模型{ModVel}，读取网格节点数MM×MN、网格大小dg、井口模型坐标MWMX信息；

S103.输入变偏移距VSP共炮点走时表{RayTS}；

S104.输入变偏移距VSP共检波点走时表{RayTR}；

S105.输入偏移孔径aper、成像起始坐标xmig1、成像结束坐标xmig2、成像网格dxmig参数、角度下限ang1、角度上限ang2。

3.根据权利要求1所述的一种变偏移距VSP带限角度积分叠前时间偏移方法，其特征在于：所述步骤S2包括：

S201.建立炮点在2维地质模型中的坐标：

MSZ_i＝SZ_i

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

S202.建议检波点在2维地质模型中的坐标：

MRZ_j＝RZ_j

其中，MRX_j是第j个检波点在2维地质模型中坐标，RX_j、RY_j是第j个检波点的大地坐标，RZ_j是第j个检波点的深度，WMX、WMY是井口的大地坐标，MWMX是井口模型坐标，sign是符号函数。

4.根据权利要求1所述的一种变偏移距VSP带限角度积分叠前时间偏移方法，其特征在于：所述步骤S3中，设计的带限角度加权系数如下：

WAngBF_ktr,(i,j)＝cos(θ_ktr,i)·cos(θ_ktr,j)ang1≤θ_ktr,i≤ang2且ang1≤θ_ktr,j≤ang2

其中，θ_ktr,i是第ktr道第i个炮点的各个散射点的射线角度，θ_ktr,j是第ktr道第j个检波点的各个散射点的射线角度，ang1是角度下限、ang2是角度上限，WAngBF_ktr,(i,j)是第ktr道第i个炮点第j个检波点的带限角度加权系数；

射线追踪计算炮点在共散射点处的射线角度的方式如下：

其中，MSX_i、MSZ_i是第i个炮点的模型坐标，x、z是散射点的坐标；

射线追踪计算检波点在共散射点处的射线角度的方式如下：

其中，MRX_i、MRZ_i是第j个检波点的模型坐标，x、z是散射点的坐标。

5.根据权利要求1所述的一种变偏移距VSP带限角度积分叠前时间偏移方法，其特征在于：所述步骤S4包括：

S401.计算步骤S1中2维深度速度模型的时深表：

ModT_k,l＝ModT_k-1,l+2·dg/ModVel_k-1,l

ModH_k,l＝ModH_k-1,l+dg

其中，ModT是模型双程垂直时，ModH是模型深度，k∈[1,MM]是模型行数，l∈[1,MN]是模型列数；

S402.用步骤S1输入的偏移参数定义成像坐标，进行如下处理：

xmig_ktr＝xmig1+(ktr-1)·dxmig ktr＝1,2,…,(xmig2-xmig1)/dxmig+1

其中，xmig_ktr是第ktr道的成像坐标，xmig1是成像起始坐标，xmig2是成像结束坐标；

MigData_ktr,j＝{0}

其中，MigData_ktr,j第j个检波点第ktr道的成像数据，先置为零；

S403.利用步骤S1输入的数据中获取偏移孔径、结合步骤S2中得到的的炮点模型坐标、步骤S402的成像坐标确定成像输入地震道范围：

xmig_ktr-aper≤MSX_kaper,j≤xmig_ktr+aper

其中，xmig_ktr是第ktr道的成像坐标，aper是偏移孔径，MSX_kaper,j是第j个检波点第ktr个成像道孔径内第kaper个炮点的模型坐标；

S404.选择第kaper个炮点xmig_ktr处的走时加上第个j检波点xmig_ktr处的走时，得到xmig_ktr处kaper、j炮检对的共散射走时TSR；

其中，

是第kaper个炮点xmig_ktr处的走时，

第j个检波点xmig_ktr处的走时；

S405.选择第j个检波点第kaper个炮点的波场数据VSPData_kaper,j及带限角度加权系数，把TSR对应的振幅带限角度加权映射为时深表对应的双程垂直，与xmig_ktr处的成像道MigData_ktr,j积分求和；

MigData_ktr,j＝MigData_ktr,j+VSPData_kaper,j·WAngBF_{ktr,(kaper,j)}

VSPData_kaper,j是第j个检波点第kaper个炮点的波场数据，WAngBF_{ktr,(kaper,j)}是第j个检波点第kaper个炮点的xmig_ktr处的带限角度加权系数；

S406.循环执行步骤S403～步骤S405，实现第j个检波点第ktr个成像道孔径内全部地震道的带限角度积分叠前时间偏移成像，得到MigData_ktr,j；

S407.循环执行步骤S402～S406，实现第j个检波点全部成像道的带限角度积分叠前时间偏移成像得到MigData_j；

S408.并行循环执行步骤S402～S407，实现全部检波点全部成像道的带限角度积分叠前时间偏移成像得到MigData。

6.一种变偏移距VSP带限角度积分叠前时间偏移装置，采用权利要求1～5中任意一项所述的方法，其特征在于：包括：

数据输入模块，用于输入2维深度速度模型、变偏移距VSP共检波点道集、变偏移距VSP共炮点走时表、变偏移距VSP共检波点走时表、偏移参数；

坐标建立模块，用于从输入的数据中获取炮点坐标、炮点深度、检波点坐标、检波点深度、井口大地坐标、井口模型坐标，建立炮点、检波点在2维地质模型中的坐标；

带限角度加权系数设计模块，用于在2维地质模型坐标系中，采用射线追踪计算炮点、检波点在共散射点处的射线角度，并从输入的数据中获取角度下限、角度上限，设计带限角度滤波器的带限角度加权系数；

积分叠前时间偏移成像模块，用于在2维地质模型中，计算时深表和共散射走时，利用得到的带限角度加权系数孔径内积分求和，实现变偏移距VSP共检波点带限角度积分叠前时间偏移成像；

共成像道集叠加模块，用于将得到的带限角度积分叠前时间偏移重排成共成像道集，共成像道集叠加，得到变偏移VSP带限角度积分叠前时间偏移叠加成像。

Images