[发明专利]用于近地表品质因子Q值反演的衰减曲线层析剥离方法

权利要求书

1.一种用于近地表品质因子Q值反演的衰减曲线层析剥离方法，包括：

基于微测井资料，建立近地表速度模型；

基于近地表速度模型，选取共检波点道集中的未受到噪声干扰的地震信号；

对所选取的地震信号进行频谱分析，以获得包含近地表结构中的多个层的品质因子Q值的等效衰减曲线；以及

基于所选取的地震信号的传播旅行时间，利用多炮组合消元方法，对等效衰减曲线进行层析剥离，以分别得到近地表结构中每层的衰减曲线。

2.根据权利要求1所述的方法，还包括：

利用谱比法，估算近地表结构中每层的品质因子Q值。

3.根据权利要求2所述的方法，其中，基于微测井资料，建立近地表速度模型，包括：

通过微测井观测系统的检波器采集经由不同激发深度处的多个炮激发而产生的地震信号，形成共检波点道集；以及

拾取共检波点道集中的地震信号到达检波器的初至时间，并且由此建立近地表速度模型。

4.根据权利要求3所述的方法，其中，建立近地表速度模型还包括：获得近地表结构的层位数量、厚度以及每层的速度，从而利用射线追踪方法计算出地震信号在每层的传播旅行时间。

5.根据权利要求4所述的方法，其中，选取共检波点道集中的未受到噪声干扰的地震信号，包括：

选取共检波点道集中的未受到噪声干扰的三个地震信号。

6.根据权利要求5所述的方法，其中，对所选取的地震信号进行频谱分析，以获得包含近地表结构中的多个层的品质因子Q值的等效衰减曲线，包括：

利用所选取的两个地震信号，获得包含近地表结构中的两层品质因子Q值的等效衰减曲线：

其中，上划线代表各变量的对数值，即表示x_i1(f)的对数值，i＝2,3,...,m，m表示微测井观测系统中激发点的个数，f表示频率，a₁代表与频率无关的常数，其中第i炮激发、第j道检波器接收的地震信号的振幅谱可表示为：

其中k表示第k层，g_ij表示与频率无关的常数项，t_ijk表示第i炮激发、第j道检波器接收的地震信号在第k层的传播旅行时间，其中j＝1。

7.根据权利要求6所述的方法，还包括：

选取两个地震信号，计算参考衰减曲线：

8.根据权利要求7所述的方法，其中，基于所选取的地震信号的传播旅行时间，利用多炮组合消元方法，对等效衰减曲线进行层析剥离，以分别得到近地表结构中每层的衰减曲线，包括：

将所述参考衰减曲线与所述等效衰减曲线的公式联立建立方程组：

并且通过求解方程组，得到仅包含一层品质因子Q的衰减曲线：

其中R₁(f)和R₂(f)分别表示求解后仅包含一层品质因子Q值的衰减曲线，c₁和c₂表示两个频率无关项。

9.根据权利要求8所述的方法，其中，利用谱比法，估算近地表结构中每层的品质因子Q值，包括：

通过分别对R₁(f)和R₂(f)线性拟合，求其斜率p₁＝-πQ₁^-1Δt₁和p₂＝-πQ₂^-1Δt₂，即可反演得到每层的品质因子Q₁、Q₂：

10.一种存储介质，其中存储有计算机可执行的程序，所述程序在被执行时适于实施用于近地表品质因子Q值反演的衰减曲线层析剥离方法，所述方法包括：

基于微测井资料，建立近地表速度模型；

基于近地表速度模型，选取共检波点道集中的未受到噪声干扰的地震信号；

对所选取的地震信号进行频谱分析，以获得包含近地表结构中的多个层的品质因子Q值的等效衰减曲线；以及

基于所选取的地震信号的传播旅行时间，利用多炮组合消元方法，对等效衰减曲线进行层析剥离，以分别得到近地表结构中每层的衰减曲线。