[发明专利]基于VSP旅行时双加权层析反演速度场的方法及其装置

权利要求书

1.一种基于VSP旅行时双加权层析反演速度场的方法，其特征在于，它包括：

在区域内采集垂直地震剖面VSP地震数据从而得到VSP地震记录的初至时间，根据所述VSP地震记录的初至时间得到所述区域的初始速度模型；

在所述区域内的各个位置根据其距井的距离以及其距地表的深度建立所述区域的可信度模型；

基于采集VSP地震数据时的观测参数、所述区域的初始速度模型和所述可信度模型得到所述区域内网格慢度扰动值；

通过所述区域内网格慢度扰动值对所述区域的初始速度模型中网格的速度进行修正得到所述区域的新的速度模型；

基于所述新的速度模型得到初至理论走时，判断所述初至理论走时与实际走时之差是否满足预定的阀值，若满足，则输出层析反演速度场。

2.根据权利要求1所述的基于VSP旅行时双加权层析反演速度场的方法，其特征在于，基于所述新的速度模型得到初至理论走时，判断所述初至理论走时与实际走时之差是否满足预定的阀值，若不满足，则基于观测参数、所述区域的新的速度模型和所述可信度模型重新得到区域内的网格慢度扰动值，进而进行循环迭代，直至所述初至理论走时与实际走时之差满足规定的阀值。

3.根据权利要求1所述的基于VSP旅行时双加权层析反演速度场的方法，其特征在于，在步骤基于采集VSP地震数据时的观测参数、所述区域的初始速度模型和所述可信度模型得到所述区域内网格慢度扰动值中，具体包括以下步骤：

根据观测参数和所述区域的初始速度模型得到每一道VSP合成地震记录的理论初至时间、每一道射线的传播路径和射线在穿过的单元网格里的长度；

基于每一道VSP合成地震记录的理论初至时间、每一道射线的传播路径、射线在穿过的单元网格里的长度和所述可信度模型得到所述区域内网格慢度扰动值。

4.根据权利要求2所述的基于VSP旅行时双加权层析反演速度场的方法，其特征在于，在步骤基于每一道VSP合成地震记录的理论初至时间、每一道射线的传播路径、射线在穿过的单元网格里的长度和所述可信度模型得到模型的网格慢度扰动值中，具体为对于有J个网格的所述区域，I条射线信息的情况时，其所述区域内网格慢度扰动值的计算公式如下：

ΔSj=ηλjΣi=1IdTi·(LijSLi)·(1Nij·SNi)SLNi]]>

λj=Σi=1ILijα,(j=1,2...J)]]>

ρi=Σj=1JLij2-α,(i=1,2...I)]]>

SLi=Σj=1JLij,(i=1,2...I)]]>

SNi=Σj=1JNij,(i=1,2...I)]]>

SLN=Σj=1j(LijSLi)·(1Nij·SNi),(i=1,2...I)]]>

其中，L_ij表示第i条射线在第j个网格上的射线长度，N_ij表示第i条射线在第j个网格上的可信度，SL_i表示第i条射线的总长度，SN_i表示第i条射线的总的可信度，η和α表示控制收敛速度与稳定性的松弛因子，I表示总射线数，J表示总网格数。