一种地震层析成像处理中的自适应加权SIRT反演方法及其系统

摘要

本发明提供了一种地震层析成像处理中的自适应加权SIRT反演方法及其系统。该反演方法及系统通过在网格慢度修正值迭代计算公式中加入了的权重修正部分，对网格的慢度修正值提供了基于射线长度的修正，解决了因为射线长度不同引入的走时误差也不同的误差问题。通过对加权因子p_i计算式中的加权幂指数h采用基于射线密度的加权幂指数h计算式，使得网格的慢度修正值c_j可以随网格射线密度的不同得到不同程度的加权修正。通过在慢度值计算单元中增加插值计算模块和平滑处理模块，对各个网格慢度值结果进行插值计算和平滑处理处理，使得该地震层析成像结果更为合理。

主权项

1.一种地震层析成像处理中的自适应加权SIRT反演方法，其特征在于，包括：获取地震层析成像处理中得到的各个炮检点的实际震动波走时数据T_im和正演理论震动波走时数据T_ic；并通过已知介质速度分布模型计算第i条射线在第j个网格内的路径长度g_ij；基于网格慢度修正值迭代计算公式对所获取的实际震动波走时数据T_im和正演理论震动波走时数据T_ic进行迭代计算，得到各个网格内的慢度修正值c_j；所述网格慢度修正值迭代计算公式为：$c_j^{\left(q\right)}=c_j^{(q-1)}+\underset{i=1}{\overset{N}{\Sigma }}\left(\frac{p_i}{L_p}\frac{g_{\mathrm{ij}}{r_i}^{(q-1)}}{\underset{j=1}{\overset{M}{\Sigma }}{g}_{\mathrm{ij}}^2}\right)$其中，j为模型网格序号；M为模型网格总数；i为射线的序号；N为射线的总条数；q为迭代计算的次数；为第j个网格的第q次迭代慢度修正值；g_ij为第i条射线在第j个网格内的长度；r_i为实际震动波走时数据T_im和正演理论震动波走时数据T_ic之差；p_i为自适应加权因子，该自适应加权因子其中，h为加权幂指数，该加权幂指数h为大于等于1的实数；L_p为所有自适应加权因子p_i的和，即$L_p=\underset{i=1}{\overset{N}{\Sigma }}{p}_i;$根据计算得到的网格的慢度修正值c_j计算该网格的慢度值s_j：