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

权利要求书

1.一种表层地震波品质因子的提取方法，其特征在于，包括：

获取待测区域内每个测点的Walkaway-VSP数据；

基于谱比法，根据所述每个测点的Walkaway-VSP数据，计算每个测点的表层地震波品质因子的值；

根据所述每个测点的表层地震波品质因子的值，以及每个测点的表层厚度，建立待测区域的表层地震波品质因子场；所述待测区域的表层地震波品质因子场，用于表征待测区域各处的表层地震波品质因子的值。

2.如权利要求1所述的方法，其特征在于，还包括：对所述每个测点的Walkaway-VSP数据预处理，具体处理包括：

对所述每个测点的Walkaway-VSP数据中的共激发点道集数据进行去噪处理；

根据去噪处理后的数据，以检波点序号重新排序，变换得到共检波点道集数据。

3.如权利要求2所述的方法，其特征在于，基于谱比法，根据所述每个测点的Walkaway-VSP数据，计算每个测点的表层地震波品质因子的值，包括：

根据每个测点的Walkaway-VSP数据变换得到的所述共检波点道集数据，拾取初至时间及续至时间，得到计算时间；

根据每个测点的Walkaway-VSP数据变换得到的所述共检波点道集数据，获得初至频谱及续至频谱，得到计算振幅谱；

根据所述计算振幅谱的对数比值计算结果，进行直线拟合得到谱比计算斜率；

根据所述计算时间及所述谱比计算斜率，按照如下公式计算得到各个续至波处的地震波品质因子的值：

Q_n＝-π(t_n+1-t₁)/k_n

其中，Q_n代表第n个续至波处的地震波品质因子值；n≥1代表可有效识别的续至波的个数；

t₁代表第一计算时间；t_n+1代表第n+1计算时间；

k_n代表第n个续至波处的谱比计算斜率；

将每个测点对应的各个续至波处的地震波品质因子的值取平均，得到每个测点的表层地震波品质因子的值。

4.如权利要求1所述的方法，其特征在于，根据所述每个测点的表层地震波品质因子的值，以及每个测点的表层厚度，建立待测区域的表层地震波品质因子场，包括：

根据所述每个测点的表层地震波品质因子的值，以及每个测点的表层厚度，拟合得到表层地震波品质因子与表层厚度的关系函数；所述关系函数用于表征表层地震波品质因子与表层厚度之间的对应关系；

将待测区域的表层厚度模型，代入所述关系函数，得到所述待测区域的表层地震波品质因子场。

5.一种表层地震波品质因子的提取装置，其特征在于，包括：

数据获取模块，用于获取待测区域内每个测点的Walkaway-VSP数据；

数据计算模块，用于基于谱比法，根据所述每个测点的Walkaway-VSP数据，计算每个测点的表层地震波品质因子的值；

表层Q场建立模块，用于根据所述每个测点的表层地震波品质因子的值以及每个测点的表层厚度，建立待测区域的表层地震波品质因子场；所述待测区域的表层地震波品质因子场，用于表征待测区域内各处的表层地震波品质因子的值。

6.如权利要求5所述的装置，其特征在于，还包括：

预处理模块，对所述每个测点的Walkaway-VSP数据预处理；

具体用于：

对所述每个测点的Walkaway-VSP数据中的共激发点道集数据进行去噪处理；

根据去噪处理后的数据，以检波点序号重新排序，变换得到共检波点道集数据。