[发明专利]一种井间多尺度SIRT层析成像的方法

权利要求书

1一种井间多尺度SIRT层析成像的方法，其特征在于，在地震资料及初始模型生成后，利用一组多尺度的网格剖分探测区域；之后，先直接采用单尺度的SIRT层析成像方法，对每一种尺度的模型分别进行层析成像，后综合各种尺度的层析成像结果，生成下一次迭代的模型，直至满足精度后，将速度模型输出。

2根据权利要求1所述的一种井间多尺度SIRT层析成像的方法，其特征在于，所述方法用于井间实际地震资料的处理，包括以下步骤，

(1)地震资料预处理步骤，用于对地震资料进行预处理；

(2)拾取信号步骤，用于拾取走时地震波信号参数；

(3)初始速度模型生成步骤，用于生成初始速度模型；

(4)建立多尺度网格步骤，用于将一组不同尺度的网格剖分探测区域；

(5)井间多尺度SIRT层析成像步骤，包括，

(5-1)单尺度SIRT层析成像过程，即对每一种尺度分别进行单尺度的SIRT层析成像，

其中，射线追踪采用基于最短路径的弯曲射线追踪方法；

SIRT层析成像公式如下：

δSj(k)=1MjΣi=1n{[ti-Σj=1maijSj(k)]aij/Σj=1maij2}---(2)]]>

式中，M_j表示穿过第j个网格的射线数，n表示射线总数，m表示模型离散化后的网格数；为了加快计算速度，并突出短射线的作用，可变为下式的平均化算法，即：

δSj(k)=Σi=1nδti(k)Ni4Li/Σi=1nNi-4---(3)]]>

式中，

δti(k)=ti-Σj=1mLiNi·Sj(k)]]>

Li=Σj=1m|aij|]]>

分别表示走时残差δt_i^(k)和第i条射线的总长度L_i，N_i表示第i条射线穿经模型的网格数；

为了让解更加平稳，对上述每次求取的慢度修正量δS_j^(k)进行光滑，其平滑公式如下：

δSi,j(k)=(δSi-1,j-1(k)+δSi-1,j(k)+δSi-1,j+1(k)+δSi,j-1(k)+δSi,j(k)+δSi,j+1(k)+δSi+1,j-1(k)+δSi+1,j(k)+δSi+1,j+1(k))/9]]>

式中，i和j分别表示模型网格划分的第i行和第j列；

(5-2)结合各种尺度的层析成像结果，采用求和平均或加权求和的方法将各个尺度的SIRT层析成像结果结合，获得综合速度模型；

(5-3)迭代修正速度模型步骤，重复(5-1)～(5-2)，直到满足精度要求；

(6)输出步骤，用于将满足精度要求的速度模型输出。