[发明专利]基于折射波分析地质体结构的方法及存储介质

权利要求书

1.一种基于折射波分析地质体结构的方法，其特征在于，包括：

S100，将目标储层的地质体剖面图区隔成多个地质网格，以形成多个地质网格点；

S200，设置各地质网格点处地震折射波的地下传播速度的初始值，以及构建地震折射波从震源到各地质网格点的射线路径，根据地震折射波从震源到各地质网格点的实际初至旅行时与地震折射波沿射线路径传播的理论初至旅行时之间的差值对地下传播速度的初始值进行修正，以获得各地质网格点处地震折射波的地下传播速度的修正值；

S300，基于地质品质因子反演模型，根据各地质网格点处地震折射波的地下传播速度的修正值和各射线路径的地质吸收特征时间，确定各地质网格点的地质品质因子；

S400，根据各地质网格点的地质品质因子的分布分析目标储层的地质体结构。

2.根据权利要求1所述的方法，其特征在于，所述步骤S200包括以下步骤：

S210，设置各地质网格点处地震折射波的地下传播速度的初始值；

S220，构建地震折射波从震源到各地质网格点的射线路径，根据地下传播速度的初始值确定地震折射波沿射线路径传播的理论初至旅行时；

S230，基于地下速度反演模型，根据地震折射波从震源到各地质网格点的实际初至旅行时与地震折射波沿射线路径传播的理论初至旅行时之间的差值对地下传播速度的初始值进行修正，以获得各地质网格点处地震折射波的地下传播速度的修正值。

3.根据权利要求1所述的方法，其特征在于，所述步骤S200包括以下步骤：

S210，设置各地质网格点处地震折射波的地下传播速度的初始值；

S220，构建地震折射波从震源到各地质网格点的射线路径，根据地下传播速度的初始值确定地震折射波沿射线路径传播的理论初至旅行时；

S230，基于地下速度反演模型，根据地震折射波从震源到各地质网格点的实际初至旅行时与地震折射波沿射线路径传播的理论初至旅行时之间的差值对地下传播速度的初始值进行修正，以获得各地质网格点处地震折射波的地下传播速度的修正值；

S240，判断对应于当前地下传播速度的修正值的理论初至旅行时与实际初至旅行时之间的差值是否小于预设的阈值：

若是，执行步骤S300；

若否，令地下传播速度的初始值等于当前地下传播速度的修正值，返回步骤S220，以对各地质网格点处地震折射波的地下传播速度再次修正。

4.根据权利要求1至3中任一项所述的方法，其特征在于，所述步骤S200中，设置各地质网格点处地震折射波的地下传播速度的初始值，使其符合均匀变化梯度的二维速度分布；通过射线追踪法构建地震折射波从震源到各地质网格点的射线路径。

5.根据权利要求1至3中任一项所述的方法，其特征在于，所述步骤S200中，地震折射波从震源到各地质网格点的实际初至旅行时为在未叠加的单炮地震记录剖面上拾取地质网格点处地震折射波的实际初至旅行时。

6.根据权利要求1至3中任一项所述的方法，其特征在于，所述步骤S400包括以下步骤：

S410，将各地质网格点的地质品质因子对应填入地质体剖面图中，以得到地质品质因子的结构分布；

S420，根据地质体剖面图中地质品质因子的分布分析目标储层的地质体结构。

7.根据权利要求2或3所述的方法，其特征在于，所述地下速度反演模型如下：

Bδv＝δt

其中，B为地下速度反演系数矩阵，δv为各地质网格点地震折射波的地下传播速度的修正值与初始值之间的差值，δt为各地质网格点的地震折射波的实际初至旅行时与理论初至旅行时之间的差值。

8.根据权利要求1所述的方法，其特征在于，所述步骤S300中，通过以下表达式获得各射线路径的地质吸收特征时间：

其中，t^*为地质吸收特征时间，C(f)为地层响应频谱，f为地震折射子波频率。