[发明专利]一种陆地地震数据的全波形反演方法

权利要求书

1.一种陆地地震数据的全波形反演方法，其特征在于，该方法具体包括以下步骤：

步骤一，建立陆地地震的光滑的初始化各向异性参数模型；

步骤二，从单炮地震记录中获取分频带的地震观测记录；

步骤三，根据参数模型和实际的地震观测系统，获取地震正演模拟合成数据集；

步骤四，计算合成地震数据与实际观测地震数据的残差δd_s,r；

步骤五，对于每一地震道，以体波第一个波包的中心作为高斯函数的中心点，计算高斯加权函数W_s,r(t)；其中，t代表时间变量；

步骤六，基于上述步骤四求取的数据残差δd_s,r，计算基于KR-OT目标函数的优化函数其中，x_r代表检波器的空间位置坐标；

步骤七，利用上述步骤五求取的高斯加权函数W_s,r(t)，计算优化后的KR-OT伴随场震源D_d；

步骤八，设定基于优化的KR-OT伴随场震源的目标函数计算目标函数参数模型更新的方向和大小，并利用目标函数参数模型代替步骤三中使用的参数模型；所述的步骤八中的基于优化的KR-OT伴随场震源的目标函数其中，max代表求取其中最大值，Δx_r代表检波器间距，代表地震观测数据集，代表地震合成数据集，下标s代表震源角标，下标r代表检波器角标；

步骤九，重复上述步骤三至步骤八，直至最终的数据残差满足精度要求，此时的参数模型即为全波形反演在当前频带所获得的陆地地震数据的各项异性参数模型。

2.根据权利要求1所述的一种陆地地震数据的全波形反演方法，其特征在于：所述的步骤二的具体步骤如下：从单炮地震记录中提取每炮记录的震源子波，并从低频到高频，对观测数据进行带通滤波，获取分频带的地震观测记录。

3.根据权利要求2所述的一种陆地地震数据的全波形反演方法，其特征在于：所述的步骤三的具体步骤如下：从最低频带开始，采用交错网格，开展复杂地质情形下的地震波正演模拟计算，空间精度8阶，时间精度2阶，获取地震正演模拟合成数据集。

4.根据权利要求3所述的一种陆地地震数据的全波形反演方法，其特征在于：所述的步骤四中的合成地震数据与实际观测地震数据的残差

5.根据权利要求4所述的一种陆地地震数据的全波形反演方法，其特征在于：所述的步骤五中的高斯加权函数其中，代表体波第一个波包的中心位置，σ代表方差。

6.根据权利要求5所述的一种陆地地震数据的全波形反演方法，其特征在于：所述的步骤七中的伴随场震源其中，为Blip1，1-Lipschitz函数。