[发明专利]地震波近地表地层品质因子的补偿方法及装置

说明书

技术领域

本发明涉及石油地震勘探技术领域，尤其涉及一种地震波近地表地层品质因子的补偿方法及装置。

背景技术

地震波在传播过程中能量不断损耗，高频能量衰减快，低频能量衰减慢。对地震波进行检测，能够得到地震数据。由于地震波能量的损耗，地震数据的分辨率随着传播时间的延长逐渐降低。为了提高地震数据的分辨率，需对地震波的高频能量进行补偿。

地震波的衰减程度与地层品质因子Q(quality，品质)有关，Q值越小，地震波衰减程度越高。Q的大小与地层岩性和地层速度有关，地层越松散，地层速度越低，Q值越小。地层结构中，近地表地层介质的松散程度较高，Q值较小，地震波的大部分高频能量在近地表被衰减，因此近地表的Q补偿十分重要。

相关技术主要通过确定性的Q补偿方法进行地震波的补偿，该方法需要计算地下介质的Q模型，而地下介质的Q模型是未知的。相关技术中的基于地震数据和VSP(Vertical Seismic profiling，垂直地震剖面)数据计算Q模型的方法都是针对中深层地层，通过相关技术无法获取近地表的Q模型，从而无法对近地表的地震波进行补偿。

发明内容

有鉴于此，本发明提供了一种地震波近地表地层品质因子的补偿方法及装置，能够对近地表的地震波进行补偿，缓解现有技术中无法对近地表的地震波进行补偿的问题。

第一方面，本发明实施例提供了一种地震波近地表地层品质因子的补偿方法，包括：在近地表布设微测井采集点，通过所述微测井采集点采集微测井数据；利用采集的所述微测井数据计算近地表的深度值与地层品质因子Q值之间的等效Q模型；根据近地表不同深度值的地震波传播速度和所述等效Q模型生成近地表的地震波传播速度与Q值之间的实测Q模型；获取待补偿地震波的振幅值、旅行时间、地震波主频和振幅调整系数，根据所述振幅值、所述旅行时间、所述地震波主频和所述振幅调整系数进行计算，得到所述待补偿地震波的相对Q模型；通过所述实测Q模型的低频分量约束所述相对Q模型的低频分量，得到最终的近地表Q模型，利用所述最终的近地表Q模型对所述待补偿地震波进行近地表Q补偿。

结合第一方面，本发明实施例提供了第一方面第一种可能的实施方式，其中，利用采集的所述微测井数据计算近地表的深度值与Q值之间的等效Q模型，包括：在采集的所述微测井数据中按照预设要求选取两个共检波点道集，并拾取每个所述共检波点道集的初至波时间，计算两个所述初至波时间的时间差T；利用S变换计算每个所述共检波点道集的初至波时刻的瞬时谱，将两个所述瞬时谱相除取对数得到谱比值；在所述谱比值中选取频率区间，利用积分法在选取的所述频率区间内分频段拟合得到谱比斜率K；根据两个所述共检波点道集对应的深度值以及所述谱比斜率K，采用统计法计算各深度值的谱比斜率K；通过公式拟合K-T曲线反演得到近地表的深度值与Q值之间的等效Q模型，其中T表示所述时间差，π表示圆周率。

结合第一方面，本发明实施例提供了第一方面第二种可能的实施方式，其中，根据近地表不同深度值的地震波传播速度和所述等效Q模型生成近地表的地震波传播速度与Q值之间的实测Q模型，包括：通过公式Q＝aV^b生成近地表的地震波传播速度与Q值之间的实测Q模型，其中Q表示所述Q值，V表示近地表不同深度值的地震波传播速度，a表示第一系数，b表示第二系数。

结合第一方面，本发明实施例提供了第一方面第三种可能的实施方式，其中，根据所述振幅值、所述旅行时间、所述地震波主频和所述振幅调整系数进行计算，得到所述待补偿地震波的相对Q模型，包括：通过公式Q＝-πft/ln(R*scale)计算所述待补偿地震波的相对Q模型，其中，Q表示所述相对Q模型，f表示所述地震波主频，t表示所述旅行时间，R表示所述振幅值，scale表示所述振幅调整系数。