地层品质因子的确定方法、装置、电子设备及存储介质与流程

1.本发明涉及油气勘探技术领域，尤其涉及地层品质因子的确定方法、装置、电子设备及存储介质。


背景技术：

2.地震波在地下介质中传播过程中，由于地层的非完全弹性地震波会产生吸收衰减，导致其能量减弱和频带变窄。地震波吸收衰减通常选用品质因子q值来描述。品质因子q值可以用于提高地震资料分辨率、储层刻画、流体检测等，因此稳健可靠地估计出品质因子q值具有重要意义。
3.当前估计地层品质因子q值的方法有谱比法、质心频移法、对数谱面积差法、峰值频率法等，其中谱比法应用最为广泛。谱比法：在一定的频带范围内获取接收子波(衰减后的地震波)和参考子波(衰减前的地震波)振幅谱比的对数，基于此进行最小二乘线性拟合，得到拟合直线的斜率，根据斜率与q值的关系式估计出q值。该方法理论上估计精度高，主要缺点是对噪声非常敏感、稳定性差。
4.受地震数据采集条件的影响，获取的实际地震数据中通常包含噪声干扰，导致求取的谱比对数与真实值存在偏差。传统的谱比法求取品质因子q值是采用最小二乘法线性拟合，在进行最小二乘运算时对拟合数据取平方会放大较大偏差数据的影响，导致其拟合结果不稳定，也就是对异常噪声非常敏感。


技术实现要素：

5.本发明提供了一种地层品质因子的确定方法、装置、电子设备及存储介质，以实现对目标地层的品质因子的确定。
6.根据本发明的一方面，提供了一种地层品质因子的确定方法，包括：
7.向目标地层激发地震波，确定参考子波、相对应的接收子波和目标频带；
8.基于所述参考子波的振幅谱和所述接收子波的振幅谱确定出谱比对数和频率值，并基于多个所述谱比对数确定谱比对数平均差以及基于多个所述频率值确定频率平均差；
9.基于所述谱比对数平均差和所述频率平均差确定所述目标地层的品质因子，其中，所述品质因子表示所述地震波在所述目标地层中吸收衰减的程度。
10.根据本发明的另一方面，提供了一种地层品质因子的确定装置，包括：
11.子波接收模块，用于向目标地层激发地震波，确定参考子波和接收子波；
12.平均差计算模块，用于基于所述参考子波的振幅谱和所述接收子波的振幅谱确定出谱比对数和频率值，并基于多个所述谱比对数确定谱比对数平均差以及基于多个所述频率值确定频率平均差；
13.因子计算模块，用于基于所述谱比对数平均差和所述频率平均差确定所述目标地层的品质因子，其中，所述品质因子表示所述地震波在所述目标地层中吸收衰减的程度。
14.根据本发明的另一方面，提供了一种电子设备，所述电子设备包括：
15.至少一个处理器；以及
16.与所述至少一个处理器通信连接的存储器；其中，
17.所述存储器存储有可被所述至少一个处理器执行的计算机程序，所述计算机程序被所述至少一个处理器执行，以使所述至少一个处理器能够执行本发明任一实施例所述的地层品质因子的确定方法。
18.根据本发明的另一方面，提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于使处理器执行时实现本发明任一实施例所述的地层品质因子的确定方法。
19.本发明实施例的技术方案，通过向目标地层激发地震波，确定参考子波、相对应的接收子波和目标频带；基于参考子波的振幅谱和接收子波的振幅谱确定出谱比对数和频率值，并基于多个谱比对数确定谱比对数平均差以及基于多个频率值确定频率平均差；基于谱比对数平均差和频率平均差确定目标地层的品质因子，其中，品质因子表示地震波在目标地层中吸收衰减的程度。本发明实施例的技术方案解决了当前对品质因子的确定方法容易受到噪声影响，导致品质因子的值不准确的问题，降低了噪声干扰，提高品质因子的准确度。
20.应当理解，本部分所描述的内容并非旨在标识本发明的实施例的关键或重要特征，也不用于限制本发明的范围。本发明的其它特征将通过以下的说明书而变得容易理解。
附图说明
21.为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
22.图1是根据本发明实施例一提供的一种地层品质因子的确定方法的流程图；
23.图2是根据本发明实施例二提供的一种地层品质因子的确定方法的流程图；
24.图3是零偏vsp直达波合成记录图；
25.图4是无噪音vsp直达波合成记录谱比法和平均差法q值估计结果图；
26.图5是无噪音vsp直达波合成记录谱比法和平均差法q值估计相对误差图；
27.图6是含噪声vsp直达波合成记录谱比法和平均差法q值估计结果图；
28.图7实际零偏vsp地震资料图；
29.图8基于实际零篇vsp地震资料谱比法和平均差法q值估计结果图；
30.图9是根据本发明实施例三提供的一种地层品质因子的确定装置的结构示意图
31.图10是实现本发明实施例四的一种地层品质因子的确定方法的电子设备的结构示意图。
具体实施方式
32.为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人
员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
33.需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
34.实施例一
35.图1是根据本发明实施例一提供的一种地层品质因子的确定方法的流程图，本实施例可适用于对地层品质因子确定的情况，该方法可以由地层品质因子的确定装置来执行，该地层品质因子的确定装置可以采用硬件和/或软件的形式实现，该地层品质因子的确定装置可配置于计算机设备或者油气勘探设备中。如图1所示，该方法包括：
36.s110、向目标地层激发地震波，确定参考子波、相对应的接收子波和目标频带。
37.在本实施例中，目标地层指的是需要对其进行品质因子确定的地层，例如，在油气勘探中，在地面(或海面)激发地震波，利用检波器在沿地层垂直方向接收地震波，即可进一步对其地层的品质因子进行确定；目标频带可以是预先设置好的频带，示例性的，目标频带为5-100hz，参考子波指的是吸收衰减前的地震波的子波，接收子波指的是与通过一些探测装置探测到的与参考子波相匹配的吸收衰减后的子波。
38.在上述方案的基础上，所述参考子波和所述接收子波的频率值处于目标频带之间。
39.在本实施例中，勘探人员可以预先设置频带的范围，并沿地层垂直方向设置检波器，通过检波器检测不同的参考子波和接收子波，示例性的，沿地层垂直方向从深度0到1000m依次布置51个检波器，相邻检波器间距为20m；需要说明的是，参考子波与接收子波的频率均是处于目标频段之间的，例如当目标频带为5-100hz，相应的，参考子波和接收子波的频率均是处于5-100hz之间的。
40.在上述方案的基础上，在所述基于所述参考子波的振幅谱和所述接收子波的振幅谱确定出谱比对数和频率值之前，还包括：对所述参考子波和所述接收子波进行傅里叶变换，得到与所述参考子波对应的振幅谱以及与所述接收子波对应的振幅谱。
41.在本实施的应用场景中，当接收到参考子波和接收子波后可以对参考子波和接收子波进行傅里叶变换，得到参考子波对应的振幅谱以及与接收子波对应的振幅谱，其中，子波吸收衰减过程的振幅谱关系可以通过式子(1)表示：
[0042][0043]
其中，a0(f)为参考子波振幅谱，a(f)为接收子波振幅谱，f为频率，δt为旅行时，q为品质因子。c表示散射衰减系数，在整个地震频带内与频率无关，包括几何扩散、反射、透射等。
[0044]
s120、基于所述参考子波的振幅谱和所述接收子波的振幅谱确定出谱比对数和频率值，并基于多个所述谱比对数确定谱比对数平均差以及基于多个所述频率值确定频率平均差。
[0045]
在本实施例中，谱比对数指的是接收子波的振幅谱与参考子波的振幅谱的比值的对数，例如，比值为a，相应的谱比对数为lna；谱比对数平均差指的是基于各个谱比对数与谱比对数平均值的差值的平均值，频率平均差可以理解为将各个频率值与频率值的算术平均数的差值的平均值。
[0046]
在上述实施例的基础上，所述基于所述参考子波的振幅谱和所述接收子波的振幅谱确定出谱比对数，包括：将所述参考子波的振幅谱作为参考子波振幅谱，并将所述接收子波的振幅谱作为接收子波振幅谱；确定所述参考子波振幅谱与所述接收子波振幅谱的比值，基于所述比值的对数确定所述谱比对数。
[0047]
具体而言，针对目标频段中，相同频率对应的参考子波和接收子波，可以将参考子波的振幅谱作为参考子波振幅谱，并将接收子波的振幅谱作为接收子波振幅谱。进一步，确定接收子波振幅谱与参考子波振幅谱的比值，将比值的对数作为谱比对数。
[0048]
示例性的，将a(f)与a0(f)相除并取对数，根据公式(1)可以得到
[0049][0050]
公式(2)是关于频率的线性函数，可以改写为
[0051]
l(f)＝kf+b
ꢀꢀꢀꢀꢀꢀ
(3)
[0052]
其中，l(f)＝ln[a(f)/a0(f)]为谱比对数；k为斜率，b为截距，b＝ln(c)。
[0053]
还需要说明的是，在目标频段内，有多个频率值对应的参考子波振幅谱和接收子波振幅谱，相应的，谱比对数也是多个，每个谱比对数均按照上述公式进行计算。
[0054]
在上述实施例的基础上，所述基于多个所述谱比对数确定谱比对数平均差，包括：基于多个所述谱比对数和所述谱比对数的总数量，计算所述谱比对数的算术平均数；基于所述谱比对数的算术平均数与每个所述谱比对数之间的差值，确定所述谱比对数平均差。
[0055]
在本实施例中，谱比对数l(f)的谱比对数平均差d可以表示为：
[0056][0057]
其中，μ表示谱比对数l(f)的算术平均数，表达式如下
[0058][0059]
将公式(3)代入公式(5)可得
[0060][0061]
其中μf表示频率的算术平均数，如下表示
[0062][0063]
在上述实施例的基础上，所述基于多个所述频率值确定频率平均差，包括：基于多个所述参考子波和所述接收子波对应的频率值以及所述频率值的总数量，确定频率值的算术平均数；基于所述频率值的算术平均数与每个频率值之间的差值，确定所述频率平均差。
[0064]
在本实施例中，上述n表示频率值的总数量，μf表示频率值的算术平均数，进一步，将公式(3)和公式(6)代入公式(4)可得：
[0065][0066]
其中df表示频率的平均差，表达式如下：
[0067][0068]
进一步，由公式(8)可得
[0069][0070]
s130、基于所述谱比对数平均差和所述频率平均差确定所述目标地层的品质因子，其中，所述品质因子表示所述地震波在所述目标地层中吸收衰减的程度。
[0071]
在上述方案的基础上，所述基于所述谱比对数平均差和所述频率平均差确定所述目标地层的品质因子，包括：基于所述谱比对数平均差与所述频率平均差的比值确定目标系数，并基于所述目标系数与旅行时，确定所述目标地层的品质因子。其中，所述旅行时为所述地震波的传输时间，目标系数指的是谱比对数平均差与频率平均差的比值，即d/df，根据目标系数与q值关系，推导出新的q值估计公式如下：
[0072]
[0073]
其中，q指的是品质因子。
[0074]
本发明实施例的技术方案，通过向目标地层激发地震波，确定参考子波、相对应的接收子波和目标频带；基于参考子波的振幅谱和接收子波的振幅谱确定出谱比对数和频率值，并基于多个谱比对数确定谱比对数平均差以及基于多个频率值确定频率平均差；基于谱比对数平均差和频率平均差确定目标地层的品质因子，其中，品质因子表示地震波在目标地层中吸收衰减的程度。本发明实施例的技术方案解决了当前对品质因子的确定方法容易受到噪声影响，导致品质因子的值不准确的问题，降低了噪声干扰，提高品质因子的准确度。
[0075]
实施例二
[0076]
图2是根据本发明实施例二提供的一种地层品质因子的确定方法的流程图，本实施例为上述实施例的优选实施例。在介绍本实施例之前，先对谱比法进行介绍：
[0077]
谱比法
[0078]
当地震波在地层中传播时，由于地层的非完全弹性，地震波会产生衰减，该过程可以用模型表示如下：
[0079][0080]
其中，a0(f)为参考子波振幅谱，a(f)为接收子波振幅谱，f为频率，δt为旅行时，q为吸收衰减参数。c表示散射衰减系数，在整个地震频带内与频率无关，包括几何扩散、反射、透射等。将a(f)与a0(f)相除并取对数，根据公式(12)可以得到
[0081][0082]
公式(13)是关于频率的线性函数，可以改写为
[0083]
l(f)＝kf+b
ꢀꢀꢀꢀꢀꢀꢀꢀ
(14)
[0084]
其中l(f)＝ln[a(f)/a0(f)]为谱比对数；k为斜率，b为截距，b＝ln(c)。
[0085]
基于公式(14)，在一定频带范围内进行最小二乘线性拟合，目标函数如下
[0086][0087]
可以得到一条直线及其斜率k。根据直线斜率k与q值关系，推导出q值估计公式
[0088][0089]
上述估计q值方法称为谱比法。
[0090]
平均差法
[0091]
受采集接收条件的影响，获取的实际地震数据中通常包含噪声干扰，导致求取的谱比对数与真实值存在偏差。传统的谱比法求取q值是采用最小二乘法线性拟合，在进行最小二乘运算时对拟合数据取平方会放大较大偏差数据的影响，导致其拟合结果不稳定，也就是对噪声非常敏感。为了提高q值估计的稳定性，提出了平均差法估计q值。平均差是各变
量与其算术平均数的离差绝对值的算术平均数。在一定频带范围内对数谱比l(f)的平均差d可以表示为
[0092][0093]
其中μ表示对数谱比l(f)的算术平均数，表达式如下
[0094][0095]
将公式(14)代入公式(18)可得
[0096][0097]
其中μf表示频率的算术平均数，如下表示
[0098][0099]
将公式(14)和公式(19)代入公式(17)可得
[0100][0101]
其中df表示频率的平均差，表达式如下
[0102][0103]
由公式(21)可得
[0104][0105]
根据斜率k与q值关系，推导出新的q值估计公式如下
[0106][0107]
上述估计q值方法称为平均差法。在平均差法求取q值的过程中，各变量与其算术平均数的离差绝对值，不会放大较大偏差数据的影响，理论上其对噪声的敏感性弱于谱比法。
[0108]
在本发明实施例的应用场景中，首先选用正演合成记录检验新方法的可靠性。设定一个水平层状介质，其模型参数如表1所示。震源为雷克子波，主频45hz；沿地层垂直方向从深度0到1000m依次布置51个检波器，相邻检波器间距为20m。在该模型下模拟得到只含有初至波的零偏vsp地震记录，如附图3所示。垂直地震剖面(vertical seismic profile，vsp)。对于附图3中合成的不含噪音零偏vsp数据，分别用谱比法和平均差法两种方法估计品质因子q值，其中平均差法估计品质因子q值的流程图如附图2所示。估计过程中频带统一选取5-100hz，估计结果如附图4所示。结果表明两种方法都能较精确地估计出q值，其相对误差如附图5所示，可见在无噪音的情况下谱比法估计结果相对误差在0.5％-0.6％之间，平均差法估计结果相对误差在0.3％-0.4％之间，平均差法的相对误差小于谱比法、q值估计精度高于谱比法。
[0109]
表1
[0110][0111]
由于vsp资料道间距很小，用相邻检波器之间的直达波来估计q值时，对随机噪音非常敏感。在附图3中的合成零偏vsp资料每道中加入均值为0、标准差为0.0005(相对于参考子波最大值的0.05％)的弱随机噪音，同样采用谱比法和平均差法两种方法估计q值，选取5-100hz频带，估计结果如附图6所示。从估计结果中可以看出，受噪声影响谱比法扰动明显，与理论值偏差较大；平均差法同样存在扰动，相对于谱比法其结果与理论值偏差明显减小。该测试结果表明在噪声的干扰下平均差法估计q值的稳定性明显高于谱比法，抗噪性增强。
[0112]
接下来选取实际资料验证新方法的可靠性，附图7是某口井的零偏vsp资料。通过波场分离获取vsp资料的下行波，然后从中拾取直达初至波，在获取的初至波中估计地层品质因子q值。附图8是选取计算频带为5-100hz时谱比法和平均差法估计的q值结果。从图中可看出，常规谱比法和平均差法估计的q值整体变化趋势比较一致，但是谱比法局部存在明显扰动，稳定性弱于平均差法。实际vsp资料应用效果进一步表明平均差法可以稳定地估计地层品质因子q值，相较于谱比法估计结果更可靠，具有很好的实用性。本发明提出的基于平均差的地层品质因子估计方法。在选取的频带范围内计算谱比对数和频率的平均差，推导出品质因子q值与平均差的定量关系式。在平均差法求取品质因子q值的过程中，各变量与其算术平均数的离差绝对值，不会放大较大偏差数据的影响，这样可以有效降低噪声干扰，提高地震波吸收衰减参数估计的稳定性。
[0113]
本发明实施例的技术方案，通过向目标地层激发地震波，确定目标频带与参考子
波和接收子波；基于参考子波的振幅谱和接收子波的振幅谱确定出谱比对数和频率值，并基于多个谱比对数确定谱比对数平均差以及基于多个频率值确定频率平均差；基于谱比对数平均差和频率平均差确定目标地层的品质因子，其中，品质因子表示地震波在目标地层中吸收衰减的程度。本发明实施例的技术方案解决了当前对品质因子的确定方法容易受到噪声影响，导致品质因子的值不准确的问题，降低了噪声干扰，提高品质因子的准确度。
[0114]
实施例三
[0115]
图9是根据本发明实施例三提供的一种地层品质因子的确定装置的结构示意图。如图9所示，该装置包括：
[0116]
子波接收模块310，用于向目标地层激发地震波，确定参考子波、相对应的接收子波和目标频带；
[0117]
平均差计算模块320，用于基于所述参考子波的振幅谱和所述接收子波的振幅谱确定出谱比对数和频率值，并基于多个所述谱比对数确定谱比对数平均差以及基于多个所述频率值确定频率平均差；
[0118]
因子计算模块330，用于基于所述谱比对数平均差和所述频率平均差确定所述目标地层的品质因子，其中，所述品质因子表示所述地震波在所述目标地层中吸收衰减的程度。
[0119]
本发明实施例的技术方案，通过向目标地层激发地震波，确定参考子波、相对应的接收子波和目标频带；基于参考子波的振幅谱和接收子波的振幅谱确定出谱比对数和频率值，并基于多个谱比对数确定谱比对数平均差以及基于多个频率值确定频率平均差；基于谱比对数平均差和频率平均差确定目标地层的品质因子，其中，品质因子表示地震波在目标地层中吸收衰减的程度。本发明实施例的技术方案解决了当前对品质因子的确定方法容易受到噪声影响，导致品质因子的值不准确的问题，降低了噪声干扰，提高品质因子的准确度。
[0120]
可选的，所述平均差计算模块320，包括：
[0121]
振幅谱确定模块，用于将所述参考子波的振幅谱作为参考子波振幅谱，并将所述接收子波的振幅谱作为接收子波振幅谱；
[0122]
谱比对数确定模块，确定所述接收子波振幅谱与所述参考子波振幅谱的比值，基于所述比值的对数确定所述谱比对数。
[0123]
可选的，所述平均差计算模块320，包括：
[0124]
谱比对数算术平均数确定模块，用于基于多个所述谱比对数和所述谱比对数的总数量，计算所述谱比对数的算术平均数；
[0125]
谱比对数平均差确定模块，用于基于所述谱比对数的算术平均数与每个所述谱比对数之间的差值，确定所述谱比对数平均差。
[0126]
可选的，所述平均差计算模块320，包括：
[0127]
频率值算术平均数确定模块，用于基于多个所述参考子波和所述接收子波对应的频率值以及所述频率值的总数量，确定频率值的算术平均数；
[0128]
频率平均差确定模块，用于基于所述频率值的算术平均数与每个频率值之间的差值，确定所述频率平均差。
[0129]
可选的，所述因子计算模块330，具体用于：
[0130]
基于所述谱比对数平均差与所述频率平均差的比值确定目标系数，并基于所述目标系数与旅行时，确定所述目标地层的品质因子；
[0131]
其中，所述旅行时为所述地震波的传输时间。
[0132]
可选的，所述子波接收模块310，包括：
[0133]
子波接收单元，用于确定参考子波和接收子波；
[0134]
可选的，还包括：
[0135]
变换模块，用于对所述参考子波和所述接收子波进行傅里叶变换，得到与所述参考子波对应的振幅谱以及与所述接收子波对应的振幅谱。
[0136]
本发明实施例所提供的地层品质因子的确定装置可执行本发明任意实施例所提供的地层品质因子的确定方法，具备执行方法相应的功能模块和有益效果。
[0137]
实施例四
[0138]
图10是实现本发明实施例四的一种地层品质因子的确定方法的电子设备的结构示意图。电子设备旨在表示各种形式的数字计算机，诸如，膝上型计算机、台式计算机、工作台、个人数字助理、服务器、刀片式服务器、大型计算机、和其它适合的计算机。电子设备还可以表示各种形式的移动装置，诸如，个人数字处理、蜂窝电话、智能电话、可穿戴设备(如头盔、眼镜、手表等)和其它类似的计算装置。本文所示的部件、它们的连接和关系、以及它们的功能仅仅作为示例，并且不意在限制本文中描述的和/或者要求的本发明的实现。
[0139]
如图10所示，电子设备10包括至少一个处理器11，以及与至少一个处理器11通信连接的存储器，如只读存储器(rom)12、随机访问存储器(ram)13等，其中，存储器存储有可被至少一个处理器执行的计算机程序，处理器11可以根据存储在只读存储器(rom)12中的计算机程序或者从存储单元18加载到随机访问存储器(ram)13中的计算机程序，来执行各种适当的动作和处理。在ram 13中，还可存储电子设备10操作所需的各种程序和数据。处理器11、rom 12以及ram 13通过总线14彼此相连。输入/输出(i/o)接口15也连接至总线14。
[0140]
电子设备10中的多个部件连接至i/o接口15，包括：输入单元16，例如键盘、鼠标等；输出单元17，例如各种类型的显示器、扬声器等；存储单元18，例如磁盘、光盘等；以及通信单元19，例如网卡、调制解调器、无线通信收发机等。通信单元19允许电子设备10通过诸如因特网的计算机网络和/或各种电信网络与其他设备交换信息/数据。
[0141]
处理器11可以是各种具有处理和计算能力的通用和/或专用处理组件。处理器11的一些示例包括但不限于中央处理单元(cpu)、图形处理单元(gpu)、各种专用的人工智能(ai)计算芯片、各种运行机器学习模型算法的处理器、数字信号处理器(dsp)、以及任何适当的处理器、控制器、微控制器等。处理器11执行上文所描述的各个方法和处理，例如地层品质因子的确定方法。
[0142]
在一些实施例中，地层品质因子的确定方法可被实现为计算机程序，其被有形地包含于计算机可读存储介质，例如存储单元18。在一些实施例中，计算机程序的部分或者全部可以经由rom 12和/或通信单元19而被载入和/或安装到电子设备10上。当计算机程序加载到ram 13并由处理器11执行时，可以执行上文描述的地层品质因子的确定方法的一个或多个步骤。备选地，在其他实施例中，处理器11可以通过其他任何适当的方式(例如，借助于固件)而被配置为执行地层品质因子的确定方法。
[0143]
本文中以上描述的系统和技术的各种实施方式可以在数字电子电路系统、集成电
路系统、场可编程门阵列(fpga)、专用集成电路(asic)、专用标准产品(assp)、芯片上系统的系统(soc)、负载可编程逻辑设备(cpld)、计算机硬件、固件、软件、和/或它们的组合中实现。这些各种实施方式可以包括：实施在一个或者多个计算机程序中，该一个或者多个计算机程序可在包括至少一个可编程处理器的可编程系统上执行和/或解释，该可编程处理器可以是专用或者通用可编程处理器，可以从存储系统、至少一个输入装置、和至少一个输出装置接收数据和指令，并且将数据和指令传输至该存储系统、该至少一个输入装置、和该至少一个输出装置。
[0144]
用于实施本发明的方法的计算机程序可以采用一个或多个编程语言的任何组合来编写。这些计算机程序可以提供给通用计算机、专用计算机或其他可编程数据处理装置的处理器，使得计算机程序当由处理器执行时使流程图和/或框图中所规定的功能/操作被实施。计算机程序可以完全在机器上执行、部分地在机器上执行，作为独立软件包部分地在机器上执行且部分地在远程机器上执行或完全在远程机器或服务器上执行。
[0145]
在本发明的上下文中，计算机可读存储介质可以是有形的介质，其可以包含或存储以供指令执行系统、装置或设备使用或与指令执行系统、装置或设备结合地使用的计算机程序。计算机可读存储介质可以包括但不限于电子的、磁性的、光学的、电磁的、红外的、或半导体系统、装置或设备，或者上述内容的任何合适组合。备选地，计算机可读存储介质可以是机器可读信号介质。机器可读存储介质的更具体示例会包括基于一个或多个线的电气连接、便携式计算机盘、硬盘、随机存取存储器(ram)、只读存储器(rom)、可擦除可编程只读存储器(eprom或快闪存储器)、光纤、便捷式紧凑盘只读存储器(cd-rom)、光学储存设备、磁储存设备、或上述内容的任何合适组合。
[0146]
为了提供与用户的交互，可以在电子设备上实施此处描述的系统和技术，该电子设备具有：用于向用户显示信息的显示装置(例如，crt(阴极射线管)或者lcd(液晶显示器)监视器)；以及键盘和指向装置(例如，鼠标或者轨迹球)，用户可以通过该键盘和该指向装置来将输入提供给电子设备。其它种类的装置还可以用于提供与用户的交互；例如，提供给用户的反馈可以是任何形式的传感反馈(例如，视觉反馈、听觉反馈、或者触觉反馈)；并且可以用任何形式(包括声输入、语音输入或者、触觉输入)来接收来自用户的输入。
[0147]
可以将此处描述的系统和技术实施在包括后台部件的计算系统(例如，作为数据服务器)、或者包括中间件部件的计算系统(例如，应用服务器)、或者包括前端部件的计算系统(例如，具有图形用户界面或者网络浏览器的用户计算机，用户可以通过该图形用户界面或者该网络浏览器来与此处描述的系统和技术的实施方式交互)、或者包括这种后台部件、中间件部件、或者前端部件的任何组合的计算系统中。可以通过任何形式或者介质的数字数据通信(例如，通信网络)来将系统的部件相互连接。通信网络的示例包括：局域网(lan)、广域网(wan)、区块链网络和互联网。
[0148]
计算系统可以包括客户端和服务器。客户端和服务器一般远离彼此并且通常通过通信网络进行交互。通过在相应的计算机上运行并且彼此具有客户端-服务器关系的计算机程序来产生客户端和服务器的关系。服务器可以是云服务器，又称为云计算服务器或云主机，是云计算服务体系中的一项主机产品，以解决了传统物理主机与vps服务中，存在的管理难度大，业务扩展性弱的缺陷。
[0149]
应该理解，可以使用上面所示的各种形式的流程，重新排序、增加或删除步骤。例
如，本发明中记载的各步骤可以并行地执行也可以顺序地执行也可以不同的次序执行，只要能够实现本发明的技术方案所期望的结果，本文在此不进行限制。
[0150]
上述具体实施方式，并不构成对本发明保护范围的限制。本领域技术人员应该明白的是，根据设计要求和其他因素，可以进行各种修改、组合、子组合和替代。任何在本发明的精神和原则之内所作的修改、等同替换和改进等，均应包含在本发明保护范围之内。

技术特征：
1.一种地层品质因子的确定方法，其特征在于，包括：向目标地层激发地震波，确定参考子波、相对应的接收子波和目标频带；基于所述参考子波的振幅谱和所述接收子波的振幅谱确定出谱比对数和频率值，并基于多个所述谱比对数确定谱比对数平均差以及基于多个所述频率值确定频率平均差；基于所述谱比对数平均差和所述频率平均差确定所述目标地层的品质因子，其中，所述品质因子表示所述地震波在所述目标地层中吸收衰减的程度。2.根据权利要求1所述的方法，其特征在于，所述基于所述参考子波的振幅谱和所述接收子波的振幅谱确定出谱比对数，包括：将所述参考子波的振幅谱作为参考子波振幅谱，并将所述接收子波的振幅谱作为接收子波振幅谱；确定所述接收子波振幅谱与所述参考子波振幅谱的比值，基于所述比值的对数确定所述谱比对数。3.根据权利要求1所述的方法，其特征在于，所述基于多个所述谱比对数确定谱比对数平均差，包括：基于多个所述谱比对数和所述谱比对数的总数量，计算所述谱比对数的算术平均数；基于所述谱比对数的算术平均数与每个所述谱比对数之间的差值，确定所述谱比对数平均差。4.根据权利要求1所述的方法，其特征在于，所述基于多个所述频率值确定频率平均差，包括：基于多个所述参考子波和所述接收子波对应的频率值以及所述频率值的总数量，确定频率值的算术平均数；基于所述频率值的算术平均数与每个频率值之间的差值，确定所述频率平均差。5.根据权利要求1所述的方法，其特征在于，所述基于所述谱比对数平均差和所述频率平均差确定所述目标地层的品质因子，包括：基于所述谱比对数平均差与所述频率平均差的比值确定目标系数，并基于所述目标系数与旅行时，确定所述目标地层的品质因子；其中，所述旅行时为所述地震波的传输时间。6.根据权利要求1所述的方法，其特征在于，所述参考子波和所述接收子波的频率值处于目标频带之间。7.根据权利要求1所述的方法，其特征在于，在所述基于所述参考子波的振幅谱和所述接收子波的振幅谱确定出谱比对数和频率值之前，还包括：对所述参考子波和所述接收子波进行傅里叶变换，得到与所述参考子波对应的振幅谱以及与所述接收子波对应的振幅谱。8.一种地层品质因子的确定装置，其特征在于，包括：子波接收模块，用于向目标地层激发地震波，确定确定参考子波、相对应的接收子波和目标频带；平均差计算模块，用于基于所述参考子波的振幅谱和所述接收子波的振幅谱确定出谱比对数和频率值，并基于多个所述谱比对数确定谱比对数平均差以及基于多个所述频率值确定频率平均差；
因子计算模块，用于基于所述谱比对数平均差和所述频率平均差确定所述目标地层的品质因子，其中，所述品质因子表示所述地震波在所述目标地层中吸收衰减的程度。9.一种电子设备，其特征在于，所述电子设备包括：至少一个处理器；以及与所述至少一个处理器通信连接的存储器；其中，所述存储器存储有可被所述至少一个处理器执行的计算机程序，所述计算机程序被所述至少一个处理器执行，以使所述至少一个处理器能够执行权利要求1-7中任一项所述的地层品质因子的确定方法。10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于使处理器执行时实现权利要求1-7中任一项所述的地层品质因子的确定方法。

技术总结
本发明实施例公开了一种地层品质因子的确定方法、装置、电子设备及存储介质。该方法包括：向目标地层激发地震波，确定参考子波、相对应的接收子波和目标频带；基于参考子波的振幅谱和接收子波的振幅谱确定出谱比对数和频率值，并基于多个谱比对数确定谱比对数平均差以及基于多个频率值确定频率平均差；基于谱比对数平均差和频率平均差确定目标地层的品质因子，其中，品质因子表示地震波在目标地层中吸收衰减的程度。本发明实施例的技术方案解决了当前对品质因子的确定方法容易受到噪声影响，导致品质因子的值不准确的问题，降低了噪声干扰，提高品质因子的准确度。提高品质因子的准确度。提高品质因子的准确度。


技术研发人员：刘军 彭光荣 杨登锋 李振升 郑金云
受保护的技术使用者：中海石油深海开发有限公司
技术研发日：2023.07.13
技术公布日：2023/10/15