[发明专利]基于三角剖分模型的正演方法及计算机可读存储介质

权利要求书

1.一种基于三角剖分模型的正演方法，其特征在于，该方法包括：

给定离散控制点，将速度模型进行Delaunay三角剖分；

在速度模型三角剖分的基础上，进行动态射线追踪；

利用高斯束传播算子获得理论地震记录。

2.根据权利要求1所述的基于三角剖分模型的正演方法，其中，所述动态射线追踪包括运动学射线追踪及动力学射线追踪。

3.根据权利要求2所述的基于三角剖分模型的正演方法，其中，所述运动学射线追踪通过公式(1)的一阶常微分方程组进行：

其中，x_i表示中心射线的坐标，p_i表示中心射线的慢度分量，s表示沿中心射线积分弧长，x_i(s)表示中心射线步长的坐标，p_i(s)表示中心射线步长的慢度分量，V(s)表示沿中心射线上模型的速度，τ表示射线追踪的时间步长。

4.根据权利要求2所述的基于三角剖分模型的正演方法，其中，所述动力学射线追踪通过公式(2)的一阶常微分方程组进行：

其中，P(s)和Q(s)为沿高斯束中心射线变的复数，s表示沿中心射线积分弧长，V(s)表示沿中心射线的模型速度，τ表示射线追踪时间步长，n表示射线中心坐标系垂直于射线方向的坐标。

5.根据权利要求1所述的基于三角剖分模型的正演方法，其中，所述利用高斯束传播算子获得理论地震记录包括：

步骤1，基于一条射线和一个检波器位置，利用公式(3)求取所有的单频的高斯束波场；

步骤2，对所有的检波器位置重复步骤1；

步骤3，对不同的射线，依次重复步骤1和步骤2，进而将所有射线对所有检波器位置的贡献累积求和；

步骤4，生成复数雷克子波，与步骤3得到的所述累积求和的结果相乘，然后进行逆傅里叶变换，获得所述理论地震记录；

其中，P(s)和Q(s)为沿高斯束中心射线变的复数，s为沿中心射线的积分弧长，u(s,n,ω)为单频的高斯束波场，τ(s)为中心射线弧长s处的传播时间，V(s)表示中心射线的速度，n为射线中心坐标系垂直中心射线的分量，i为虚数单位，ω为角频率。

6.根据权利要求5所述的基于三角剖分模型的正演方法，其中，所述检波器位置的获得包括：利用给定检波点位置，通过插值得到均匀分布的接收位置。

7.根据权利要求1所述的基于三角剖分模型的正演方法，其中，所述动态射线追踪包括：

步骤1，确定震源所在的三角形网格；

步骤2，确定一个出射角度，进行动态射线追踪；

步骤3，遇到所述三角形网格边界，读取所述三角形网格边界属性，确定遇到的所述三角形网格边界为反射、透射或者终止追踪，如果所述三角形网格边界为反射或者透射，计算反射或者透射对振幅的影响；

步骤4，记录中心射线路径每个步长上的变量，所述变量包括中心射线位置、慢度分量、动态射线追踪参数；

步骤5，改变出射角度，重复步骤2-步骤4，直至遍历所有出射角度。

8.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述程序被处理器执行时实现以下步骤：

步骤1，给定离散控制点，将速度模型进行Delaunay三角剖分；

步骤2，在速度模型三角剖分的基础上，进行动态射线追踪；

步骤3，利用高斯束传播算子获得理论地震记录。