[发明专利]一种地层品质因子的计算方法及装置

权利要求书

1.一种地层品质因子的计算方法，其特征在于，所述方法包括：

根据地震数据分别计算第一时刻对应的第一振幅谱和第二时刻对应的第二振幅谱；

对所述第一振幅谱、所述第二振幅谱分别进行低通滤波处理，分别获得对应的地震子波振幅谱的第一估计值、第二估计值；

根据所述第一估计值和所述第二估计值，采用谱比法计算地层品质因子；

其中，所述对所述第一振幅谱、所述第二振幅谱分别进行低通滤波处理，包括：

计算所述第一振幅谱对应的第一频宽和所述第二振幅谱对应的第二频宽；

根据所述第一频宽设置所述第一振幅谱对应的第一低通滤波器的截止频率，根据所述第二频宽设置所述第二振幅谱对应的第二低通滤波器的截止频率；

所述根据所述第一频宽设置所述第一振幅谱对应的第一低通滤波器的截止频率，根据所述第二频宽设置所述第二振幅谱对应的第二低通滤波器的截止频率，包括：

分别设置频宽阈值一对应的第一截止频率，和频宽阈值二对应的第二截止频率，所述频宽阈值一小于所述频宽阈值二，利用如下公式计算所述第一低通滤波器的截止频率，和所述第二低通滤波器的截止频率：

式中，FEX1表示所述第一低通滤波器的截止频率，FEX2表示所述第二低通滤波器的截止频率，f1表示所述第一截止频率，f2表示所述第二截止频率，BW1表示所述第一频宽，BW2表示所述第二频宽，B1表示所述频宽阈值一，B2表示所述频宽阈值二。

2.如权利要求1所述的一种地层品质因子的计算方法，其特征在于，所述计算所述第一振幅谱对应的第一频宽和所述第二振幅谱对应的第二频宽，包括：

利用第三低通滤波器对所述第一振幅谱进行平滑处理，获得第一光滑振幅谱；

根据所述第一光滑振幅谱计算所述第一频宽；

利用第四低通滤波器对所述第二振幅谱进行平滑处理，获得第二光滑振幅谱；

根据所述第二光滑振幅谱计算所述第二频宽。

3.如权利要求1所述的一种地层品质因子的计算方法，其特征在于，所述根据所述第一频宽设置所述第一振幅谱对应的第一低通滤波器的截止频率，根据所述第二频宽设置所述第二振幅谱对应的第二低通滤波器的截止频率，包括：

若所述第一频宽小于所述频宽阈值一，则将所述第一截止频率作为所述第一低通滤波器的截止频率；

若所述第二频宽小于所述频宽阈值一，则将所述第一截止频率作为所述第二低通滤波器的截止频率；

若所述第一频宽大于所述频宽阈值二，则将所述第二截止频率作为所述第一低通滤波器的截止频率；

若所述第二频宽大于所述频宽阈值二，则将所述第二截止频率作为所述第一低通滤波器的截止频率。

4.如权利要求1所述的一种地层品质因子的计算方法，其特征在于，所述根据所述第一估计值和所述第二估计值，采用谱比法计算地层的品质因子，包括：

根据所述第一估计值和所述第二估计值，计算所述第一时刻和所述第二时刻对应的地震子波的振幅谱比；

根据所述振幅谱比获得地震子波的振幅谱比的自然对数曲线；

根据所述第一估计值或所述第二估计值，确定所述自然对数曲线中的有效范围；

在所述有效范围内对所述自然对数曲线进行最小平方直线拟合，获得拟合后直线的斜率；

根据所述斜率，计算获得所述地层品质因子。

5.如权利要求4所述的一种地层品质因子的计算方法，其特征在于，所述根据所述第一估计值或所述第二估计值，确定所述自然对数曲线中的有效范围，包括：

获取所述第一估计值对应的最大振幅值作为第一最大值，或者获取所述第二估计值对应的最大振幅值作为第二最大值；

在所述自然对数曲线中以所述第一最大值对应的第一频率或以所述第二最大值对应的第二频率作为起始点，分别向低频端和高频端搜索，确定低频起始点和高频终止点；

将所述低频起始点和所述高频终止点之间的范围作为所述有效范围。