CN112099083B - 一种基于双谱谱比对数的品质因子估计方法及系统 - Google Patents
一种基于双谱谱比对数的品质因子估计方法及系统 Download PDFInfo
- Publication number
- CN112099083B CN112099083B CN202010870741.0A CN202010870741A CN112099083B CN 112099083 B CN112099083 B CN 112099083B CN 202010870741 A CN202010870741 A CN 202010870741A CN 112099083 B CN112099083 B CN 112099083B
- Authority
- CN
- China
- Prior art keywords
- quality factor
- spectrum
- bispectral
- calculating
- downlink
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 82
- 238000000034 method Methods 0.000 title claims abstract description 35
- 230000009977 dual effect Effects 0.000 claims abstract description 16
- 239000011229 interlayer Substances 0.000 claims abstract description 13
- 238000000926 separation method Methods 0.000 claims abstract description 11
- 238000004364 calculation method Methods 0.000 claims description 7
- 230000000694 effects Effects 0.000 claims description 7
- 238000009825 accumulation Methods 0.000 claims description 6
- 230000001186 cumulative effect Effects 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 5
- 230000003595 spectral effect Effects 0.000 claims description 5
- 238000010521 absorption reaction Methods 0.000 claims description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 4
- 238000007405 data analysis Methods 0.000 claims description 4
- 238000009792 diffusion process Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 238000007619 statistical method Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 5
- 238000004891 communication Methods 0.000 description 5
- 239000010410 layer Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000002301 combined effect Effects 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 238000011545 laboratory measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000010587 phase diagram Methods 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明实施例提供一种基于双谱谱比对数的品质因子估计方法及系统，该方法包括：获取零偏垂直地震剖面全波场数据，并对所述零偏垂直地震剖面全波场数据进行波场分离，获取下行波场；在所述下行波场中求取下行直达波，并计算含噪地震波的双谱；计算相邻道所述双谱谱比的对数，并基于相邻道所述双谱谱比的对数，估计层间品质因子。本发明实施例通过引入高阶统计分析实现有效信号地双谱估计，进而通过双谱谱比对数地平面拟合实现品质因子地求取，从而削弱随机噪声地影响，提高品质因子估计的稳定性及精度。
Description
技术领域
本发明涉及地球物理勘探技术领域，尤其涉及一种基于双谱谱比对数的品质因子估计方法及系统。
背景技术
品质因子Q描述了地下介质的粘滞性特征，即由于介质颗粒摩擦、不可逆热传导、晶体不完整性以及孔隙流体流动等作用地综合影响，地震波在传播过程中发生衰减，主要体现在：1、耗散效应，即能量衰减特性，表现为地震波的振幅减小，具体到单频成分，高频分量的振幅衰减程度要强于低频分量；2、频散效应，即相变特性，表现为高频分量的传播速度快于低频分量，导致地震波的相位特征发生变化。这两种效应综合在一起，在地震资料中表现为随着时间地增加，有效信号的能量及主频逐渐减小，从而影响了深部资料的信噪比及分辨率。
为了消除品质因子的影响，提高深部资料有效信号的质量，在品质因子已知的条件下，反Q滤波处理是常用的技术手段。因此，地层品质因子地准确获取是问题解决的关键。目前品质因子地求解主要有两种手段：1、实验室测量，针对地下岩心，在不同的温度、压力、流体等条件下，采用物理手段测量品质因子，这种方式得到的品质因子精度较高，但由于岩心数量地限制，往往只能得到地层品质因子的有限离散点，难以获取较宽深度范围内品质因子地连续变化；2、基于实际地震资料地品质因子估计，采用这种方式，可以获取满足资料处理及解释要求的品质因子序列，但精度容易受到地震资料信噪比地影响。通常可用的资料包括地面地震资料、垂直地震剖面(Vertical Seismic Profile，简称VSP)资料以及微测井资料，其中，VSP资料由于具有较高的信噪比以及较深的采集深度，成为求取品质因子最常用的地震资料。
到目前为止，发展了多种品质因子估计方法，包括振幅衰减、上升时间、子波模拟、解析信号、波形匹配、谱模拟、振幅谱比、质心频移和峰值频移等方法，分别从时间域和Fourier频率域中力求实现品质因子地可靠求取，其中后三种方法在实际中得到了广泛应用，但求解精度及稳定性有待进一步提高。因此，现在亟需一种基于双谱谱比对数的品质因子估计方法及系统来解决上述问题。
发明内容
针对现有技术存在的问题，本发明实施例提供一种基于双谱谱比对数的品质因子估计方法及系统。
第一方面，本发明实施例提供了一种基于双谱谱比对数的品质因子估计方法，包括：
获取零偏垂直地震剖面全波场数据，并对所述零偏垂直地震剖面全波场数据进行波场分离，获取下行波场；
在所述下行波场中求取下行直达波，并计算含噪地震波的双谱；
计算相邻道所述双谱谱比的对数，并基于相邻道所述双谱谱比的对数，估计层间品质因子。
进一步地，所述方法还包括：
对所述零偏垂直地震剖面全波场数据进行去噪和几何扩散补偿。
进一步地，所述在所述下行波场中求取下行直达波，并计算含噪地震波的双谱，包括：
基于高斯分布随机信号的高阶累积量为零，消除高斯噪声对下行波场的影响，得到含噪地震波的双谱。
进一步地，所述计算含噪地震波的双谱为：
其中，ω表示角频率，ωh表示与截止频率相关的调节参数，γ＝(πQr)-1，Qr表示参考频率处的品质因子，为无噪波形的三阶累积量。
进一步地，所述计算相邻道所述双谱谱比的对数，包括：
基于所述双谱和Fourier谱之间的关系，计算t+Δt时刻和t时刻之间的双谱振幅谱关系为：
其中，ω表示角频率，ωh表示与截止频率相关的调节参数，γ＝(πQr)-1，Qr表示参考频率处的品质因子，t和Δt表示时间。
进一步地，所述计算相邻道所述双谱谱比的对数，并基于相邻道所述双谱谱比的对数，估计层间品质因子为：
其中，ω表示角频率，ωh表示与截止频率相关的调节参数，γ＝(πQr)-1，Qr表示参考频率处的品质因子，t和Δt表示时间。
第二方面，本发明实施例提供了一种基于双谱谱比对数的品质因子估计系统，包括：
数据解析模块，用于获取零偏垂直地震剖面全波场数据，并对所述零偏垂直地震剖面全波场数据进行波场分离，获取下行波场；
双谱计算模块，用于在所述下行波场中求取下行直达波，并计算含噪地震波的双谱；
品质因子估计模块，用于计算相邻道所述双谱谱比的对数，并基于相邻道所述双谱谱比的对数，估计层间品质因子。
第三方面，本发明实施例提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如第一方面所提供的方法的步骤。
第四方面，本发明实施例提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现如第一方面所提供的方法的步骤。
本发明实施例提供的一种基于双谱谱比对数的品质因子估计方法及系统，通过引入高阶统计分析实现有效信号地双谱估计，进而通过双谱谱比对数地平面拟合实现品质因子地求取，从而削弱随机噪声地影响，提高品质因子估计的稳定性及精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的基于双谱谱比对数的品质因子估计方法的流程示意图；
图2为本发明实施例提供的粘弹性零偏VSP下行直达波示意图；
图3为本发明实施例提供的无噪资料品质因子以及含噪资料品质因子估计结果；
图4为本发明实施例提供的基于双谱谱比对数的品质因子估计系统的结构示意图；
图5为本发明实施例提供的电子设备结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
图1为本发明实施例提供的基于双谱谱比对数的品质因子估计方法的流程示意图，如图1所示，本发明实施例提供了一种基于双谱谱比对数的品质因子估计方法，包括：
步骤101，获取零偏垂直地震剖面全波场数据，并对所述零偏垂直地震剖面全波场数据进行波场分离，获取下行波场；
步骤102，在所述下行波场中求取下行直达波，并计算含噪地震波的双谱；
步骤103，计算相邻道所述双谱谱比的对数，并基于相邻道所述双谱谱比的对数，估计层间品质因子。
在本发明实施例中，首先会获取零偏VSP全波场资料，从而对其进行波场分离得到下行波场；然后，在下行波场中求取其中的下行直达波，从而根据下行直达波计算三阶谱，也就是本发明实施例中的双谱；最后，计算相邻道双谱谱比的对数，从而估计出层间品质因子。
本发明实施例提供的基于双谱谱比对数的品质因子估计方法，通过引入高阶统计分析实现有效信号地双谱估计，进而通过双谱谱比对数地平面拟合实现品质因子地求取，从而削弱随机噪声地影响，提高品质因子估计的稳定性及精度。
在上述实施例的基础上，所述方法还包括：
对所述零偏垂直地震剖面全波场数据进行去噪和几何扩散补偿。
在本发明实施例中，需要对得到的初始数据进行数据预处理，具体包括去噪和几何扩散补偿等操作。
在上述实施例的基础上，所述在所述下行波场中求取下行直达波，并计算含噪地震波的双谱，包括：
基于高斯分布随机信号的高阶累积量为零，消除高斯噪声对下行波场的影响，得到含噪地震波的双谱。
在上述实施例的基础上，所述计算含噪地震波的双谱为：
其中，ω表示角频率，ωh表示与截止频率相关的调节参数，γ＝(πQr)-1，Qr表示参考频率处的品质因子，为无噪波形的三阶累积量。
具体地，在地层吸收衰减效应地影响下，传播过程中的地震波可以表示为：
其中，t和Δt表示时间，ω表示角频率，ωh表示与截止频率相关的调节参数，γ＝(πQr)-1，Qr表示参考频率处的品质因子，U(t,ω)表示位于t时刻、频率为ω的地震波的单频分量，u(t)表示t时刻地震波的时间域波形，n(Δt)表示噪声，i表示虚数单位。
进一步地，在本发明实施例中，公式中的第一个指数项表示振幅衰减，第二个指数项表示相位变化。由于ω/ωh<1，γ<<1，因此考虑到问题求解的方便，可以忽略振幅衰减项中的(ω/ωh)-γ。
当公式中的噪声呈高斯分布时，为了消除其影响，计算地震信号的三阶累积量：
其中，τ1和τ2分别表示延迟量，c3u、和c3n分别表示含噪波形、无噪波形和高斯噪声的三阶累积量。由于高斯分布随机信号的高阶累积量为零，即中c3n(τ1,τ2)＝0，因此：
从而得到含噪地震波的三阶谱，即双谱：
在上述实施例的基础上，所述计算相邻道所述双谱谱比的对数，包括：
基于所述双谱和Fourier谱之间的关系，计算t+Δt时刻和t时刻之间的双谱振幅谱关系为：
其中，ω表示角频率，ωh表示与截止频率相关的调节参数，γ＝(πQr)-1，Qr表示参考频率处的品质因子，t和Δt表示时间。
具体地，考虑到双谱与Fourier谱之间的关系：
其中，N表示信号的长度，*表示共轭运算。通过对公式两侧取振幅谱，并忽略系数项，则得到：
Bu,amp(ω1,ω2)＝Uamp(ω1)Uamp(ω2)Uamp(ω1+ω2)；
其中，下标amp表示振幅谱。由公式可知，t+Δt时刻和t时刻之间的Fourier振幅谱关系为：
因此，t+Δt时刻和t时刻之间的双谱振幅谱关系为：
在上述实施例的基础上，所述计算相邻道所述双谱谱比的对数，并基于相邻道所述双谱的谱比对数，估计层间品质因子为：
其中，ω表示角频率，ωh表示与截止频率相关的调节参数，γ＝(πQr)-1，Qr表示参考频率处的品质因子，t和Δt表示时间。
在本发明实施例中，即为基于双谱谱比对数的品质因子求解表达式，其表示关于频率ω1和ω2的平面，通过最小二乘拟合即可求取平面斜率-Δt/Q，进而获得地震波在传播时间Δt内所穿过地层的平均Q值。
具体地，本发明实施例以Ricker子波为例，设定其初始时刻为零相位、主频为30Hz，通过对Ricker子波波形图以及对应的双谱振幅谱图进行分析可知，该子波在品质因子为50的均匀地层中传播了50ms后，经过了地层的吸收衰减效应，波形发生了变化，具体表现为振幅变小，并且时间域波形由对称变成了非对称，即相位特征由零相位变成了非零相位；进一步地，对Ricker子波对应的非零相位图以及对应的双谱振幅谱图进行分析可知，将Ricker子波的衰减后与衰减前的双谱振幅谱按照双谱谱比对数，即可得到双谱谱比对数与频率ω1和ω2之间的关系。具体地，在本发明实施例中，对双谱谱比对数与频率之间的关系进行分析，双频率域中的具有一定斜率的平面，通过最小二乘拟合得到平面的斜率为-Δt/Q＝-0.001，由于Δt＝0.05s，即可求得地层的品质因子Q＝50。
进一步地，本发明实施例提供了一种用于粘弹性零偏VSP正演模拟的速度模型和品质因子模型，速度模型和品质因子模型均设定为6层，每层厚度200m，层内为均匀介质，其中，速度模型的变化范围为500m/s～2000m/s，品质因子模型相应的变化范围为20～120。炮点坐标为(25m,0m)，检波点横坐标为25m，纵坐标从100m变化到1200m，间距为5m。
进一步地，基于上述实施例中的模型参数，进行粘弹性零偏VSP正演模拟，分别得到全波场正演模拟结果和对应的下行波场，全波场对应的上行波场以及下行直达波初至，图2为本发明实施例提供的粘弹性零偏VSP下行直达波，其中，下行波场中的直达波是品质因子求解所要利用的数据。通过在下行波场中拾取初至时间，获取各道直达波，在此基础上计算地层品质因子。图3为本发明实施例提供的无噪资料品质因子以及含噪资料品质因子估计结果，其中，在图3中，(a)为无噪资料的计算结果，可以看到估计值与真实值具有较好的吻合度，(b)为含噪资料的计算结果，其中噪声能量占有效信号能量的0.3％，可以看到，当品质因子较小时，由于地震波衰减较为剧烈，波形变化占主要成分，品质因子估计结果的精度和稳定性较高，而随着品质因子的增大，地震波衰减程度变小，此时噪声成分逐渐显现出来，从而造成品质因子估计结果与真实值有所偏离，但整体来看，品质因子估计结果与真实值的误差在可控范围之内。
图4为本发明实施例提供的基于双谱谱比对数的品质因子估计系统的结构示意图，如图4所示，本发明实施例提供了一种基于双谱谱比对数的品质因子估计系统，包括数据解析模块401、双谱计算模块402和品质因子估计模块403，其中，数据解析模块401用于获取零偏垂直地震剖面全波场数据，并对所述零偏垂直地震剖面全波场数据进行波场分离，获取下行波场；双谱计算模块402用于在所述下行波场中求取下行直达波，并计算含噪地震波的双谱；品质因子估计模块403用于计算相邻道所述双谱谱比的对数，并基于相邻道所述双谱谱比的对数，估计层间品质因子。
本发明实施例提供的基于双谱谱比对数的品质因子估计系统，通过引入高阶统计分析实现有效信号地双谱估计，进而通过双谱谱比对数地平面拟合实现品质因子地求取，从而削弱随机噪声地影响，提高品质因子估计的稳定性及精度。
本发明实施例提供的系统是用于执行上述各方法实施例的，具体流程和详细内容请参照上述实施例，此处不再赘述。
图5为本发明实施例提供的电子设备结构示意图，参照图5，该电子设备可以包括：处理器(processor)501、通信接口(Communications Interface)502、存储器(memory)503和通信总线504，其中，处理器501，通信接口502，存储器503通过通信总线504完成相互间的通信。处理器501可以调用存储器503中的逻辑指令，以执行如下方法：获取零偏垂直地震剖面全波场数据，并对所述零偏垂直地震剖面全波场数据进行波场分离，获取下行波场；在所述下行波场中求取下行直达波，并计算含噪地震波的双谱；计算相邻道所述双谱谱比的对数，并基于相邻道所述双谱谱比的对数，估计层间品质因子。
此外，上述的存储器503中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
另一方面，本发明实施例还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现以执行上述各实施例提供的基于双谱谱比对数的品质因子估计方法，例如包括：获取零偏垂直地震剖面全波场数据，并对所述零偏垂直地震剖面全波场数据进行波场分离，获取下行波场；在所述下行波场中求取下行直达波，并计算含噪地震波的双谱；计算相邻道所述双谱谱比的对数，并基于相邻道所述双谱谱比的对数，估计层间品质因子。
以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (5)
1.一种基于双谱谱比对数的品质因子估计方法，其特征在于，包括：
获取零偏垂直地震剖面全波场数据，并对所述零偏垂直地震剖面全波场数据进行波场分离，获取下行波场；
在所述下行波场中求取下行直达波，并计算含噪地震波的双谱；
计算相邻道所述双谱谱比的对数，并基于相邻道所述双谱谱比的对数，估计层间品质因子；
所述在所述下行波场中求取下行直达波，并计算含噪地震波的双谱，包括：
基于高斯分布随机信号的高阶累积量为零，消除高斯噪声对下行波场的影响，得到含噪地震波的双谱；
所述计算含噪地震波的双谱为：
在地层吸收衰减效应地影响下，传播过程中的地震波表示为：
ω/ωh<1，γ<<1；
其中，ω表示角频率，ωh表示与截止频率相关的调节参数，γ＝(πQr)-1，Qr表示参考频率处的品质因子，Q表示地层品质因子，为无噪波形的三阶累积量；t和Δt表示时间；U(t,ω)表示位于t时刻、频率为ω的地震波的单频分量，u(t)表示t时刻地震波的时间域波形，n(Δt)表示高斯噪声，i表示虚数单位；τ1和τ2分别表示延迟量，c3u、/>分别表示含噪波形、无噪波形的三阶累积量；
所述计算相邻道所述双谱谱比的对数，包括：
基于所述双谱和Fourier谱之间的关系，计算t+Δt时刻和t时刻之间的双谱振幅谱关系为：
其中，ω表示角频率，Bu,amp表示双谱振幅谱，Uamp表示Fourier振幅谱，t和Δt表示时间；
所述计算相邻道所述双谱谱比的对数，并基于相邻道所述双谱谱比的对数，估计层间品质因子为：
其中，ω表示角频率，t和Δt表示时间。
2.根据权利要求1所述的基于双谱谱比对数的品质因子估计方法，其特征在于，所述方法还包括：
对所述零偏垂直地震剖面全波场数据进行去噪和几何扩散补偿。
3.一种基于双谱谱比对数的品质因子估计系统，其特征在于，包括：
数据解析模块，用于获取零偏垂直地震剖面全波场数据，并对所述零偏垂直地震剖面全波场数据进行波场分离，获取下行波场；
双谱计算模块，用于在所述下行波场中求取下行直达波，并计算含噪地震波的双谱；
品质因子估计模块，用于计算相邻道所述双谱谱比的对数，并基于相邻道所述双谱谱比的对数，估计层间品质因子；
所述双谱计算模块具体用于：
基于高斯分布随机信号的高阶累积量为零，消除高斯噪声对下行波场的影响，得到含噪地震波的双谱；
所述计算含噪地震波的双谱为：
在地层吸收衰减效应地影响下，传播过程中的地震波表示为：
ω/ωh<1，γ<<1；
其中，ω表示角频率，ωh表示与截止频率相关的调节参数，γ＝(πQr)-1，Qr表示参考频率处的品质因子，Q表示地层品质因子，为无噪波形的三阶累积量；t和Δt表示时间；U(t,ω)表示位于t时刻、频率为ω的地震波的单频分量，u(t)表示t时刻地震波的时间域波形，n(Δt)表示高斯噪声，i表示虚数单位；τ1和τ2分别表示延迟量，c3u、/>分别表示含噪波形、无噪波形的三阶累积量；
所述计算相邻道所述双谱谱比的对数，包括：
基于所述双谱和Fourier谱之间的关系，计算t+Δt时刻和t时刻之间的双谱振幅谱关系为：
其中，ω表示角频率，Bu,amp表示双谱振幅谱，Uamp表示Fourier振幅谱，t和Δt表示时间；
所述品质因子估计模块通过以下表达式进行计算：
其中，ω表示角频率，t和Δt表示时间。
4.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现如权利要求1至2任一项所述基于双谱谱比对数的品质因子估计方法的步骤。
5.一种非暂态计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现如权利要求1至2任一项所述基于双谱谱比对数的品质因子估计方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010870741.0A CN112099083B (zh) | 2020-08-26 | 2020-08-26 | 一种基于双谱谱比对数的品质因子估计方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010870741.0A CN112099083B (zh) | 2020-08-26 | 2020-08-26 | 一种基于双谱谱比对数的品质因子估计方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112099083A CN112099083A (zh) | 2020-12-18 |
| CN112099083B true CN112099083B (zh) | 2023-10-13 |
Family
ID=73756629
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010870741.0A Active CN112099083B (zh) | 2020-08-26 | 2020-08-26 | 一种基于双谱谱比对数的品质因子估计方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112099083B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113625341A (zh) * | 2021-08-04 | 2021-11-09 | 中化地质矿山总局地质研究院 | 一种基于复赛谱分析的品质因子估计方法、装置及系统 |
| CN117055113B (zh) * | 2023-08-18 | 2024-07-05 | 中国矿业大学(北京) | 品质因子提取的方法、装置及电子设备 |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103376464A (zh) * | 2012-04-13 | 2013-10-30 | 中国石油天然气集团公司 | 一种地层品质因子反演方法 |
| CN103744114A (zh) * | 2013-12-30 | 2014-04-23 | 中国石油天然气股份有限公司 | 基于零偏垂直地震剖面数据估计品质因子的方法和装置 |
| CN103995289A (zh) * | 2014-05-19 | 2014-08-20 | 中国石油大学(华东) | 基于时频谱模拟的时变混合相位地震子波提取方法 |
| CN104199095A (zh) * | 2014-09-11 | 2014-12-10 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 提高地震记录分辨率的反褶积方法 |
| CN105005073A (zh) * | 2015-06-11 | 2015-10-28 | 中国石油大学(华东) | 基于局部相似度和评价反馈的时变子波提取方法 |
| CN106324663A (zh) * | 2015-06-17 | 2017-01-11 | 中国石油化工股份有限公司 | 一种品质因子的获取方法 |
| CN108387933A (zh) * | 2018-02-01 | 2018-08-10 | 中国石油天然气股份有限公司 | 一种确定地层品质因子的方法、装置及系统 |
| CN108845357A (zh) * | 2018-06-13 | 2018-11-20 | 成都信息工程大学 | 一种基于同步挤压小波变换估计地层等效品质因子的方法 |
| CN109765615A (zh) * | 2019-01-10 | 2019-05-17 | 中国石油天然气股份有限公司 | 一种地层品质因子反演方法及装置 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102010041920A1 (de) * | 2010-10-04 | 2012-04-05 | Siemens Aktiengesellschaft | Verfahren zur Darstellung einer Konzentration eines Kontrastmittels in einem vorbestimmten Volumenabschnitt mittels Tomosynthese und entsprechendes Tomosynthesegerät |
-
2020
- 2020-08-26 CN CN202010870741.0A patent/CN112099083B/zh active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103376464A (zh) * | 2012-04-13 | 2013-10-30 | 中国石油天然气集团公司 | 一种地层品质因子反演方法 |
| CN103744114A (zh) * | 2013-12-30 | 2014-04-23 | 中国石油天然气股份有限公司 | 基于零偏垂直地震剖面数据估计品质因子的方法和装置 |
| CN103995289A (zh) * | 2014-05-19 | 2014-08-20 | 中国石油大学(华东) | 基于时频谱模拟的时变混合相位地震子波提取方法 |
| CN104199095A (zh) * | 2014-09-11 | 2014-12-10 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 提高地震记录分辨率的反褶积方法 |
| CN105005073A (zh) * | 2015-06-11 | 2015-10-28 | 中国石油大学(华东) | 基于局部相似度和评价反馈的时变子波提取方法 |
| CN106324663A (zh) * | 2015-06-17 | 2017-01-11 | 中国石油化工股份有限公司 | 一种品质因子的获取方法 |
| CN108387933A (zh) * | 2018-02-01 | 2018-08-10 | 中国石油天然气股份有限公司 | 一种确定地层品质因子的方法、装置及系统 |
| CN108845357A (zh) * | 2018-06-13 | 2018-11-20 | 成都信息工程大学 | 一种基于同步挤压小波变换估计地层等效品质因子的方法 |
| CN109765615A (zh) * | 2019-01-10 | 2019-05-17 | 中国石油天然气股份有限公司 | 一种地层品质因子反演方法及装置 |
Non-Patent Citations (2)
| Title |
|---|
| 提高地震分辨率处理效果定量评价方法研究;刁瑞;;物探与化探(第02期);全文 * |
| 高阶统计量方法在地震勘探中的应用;尹成,伍志明,邓怀群;地球物理学进展(第03期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112099083A (zh) | 2020-12-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20150168573A1 (en) | Geologic quality factor inversion method | |
| Wang et al. | Q factor estimation based on the method of logarithmic spectral area difference | |
| CN112099083B (zh) | 一种基于双谱谱比对数的品质因子估计方法及系统 | |
| CN109669212B (zh) | 地震数据处理方法、地层品质因子估算方法与装置 | |
| CN107132577B (zh) | 一种基于谱面积变化的地震衰减估计方法 | |
| CN106019376B (zh) | 一种频率驱动空变q值模型构建的地震波补偿方法 | |
| CN108646289B (zh) | 一种估计地震品质因子的方法 | |
| CN112285778B (zh) | 一种粘声TTI介质中纯qP波的逆时偏移成像方法 | |
| Tu et al. | Improve Q estimates with spectrum correction based on seismic wavelet estimation | |
| CN111427088A (zh) | 用于识别薄互储层的地震数据低频补偿方法 | |
| CN112904415B (zh) | 基于重构源谱一致性的q层析反演方法、系统及处理设备 | |
| CN107238864A (zh) | 基于优势频带约束的地震波能量时频补偿方法 | |
| CN105093291A (zh) | 一种恢复油气储层地震反射特征的方法 | |
| CN109143345A (zh) | 基于模拟退火的品质因子q非线性反演方法及系统 | |
| Wang et al. | Application of a new wavelet threshold method in unconventional oil and gas reservoir seismic data denoising | |
| CN110568491B (zh) | 一种品质因子q的估算方法 | |
| Ben‐Feng et al. | A stable and efficient attenuation compensation method based on inversion | |
| CN111427080B (zh) | 地震数据空变广义子波的提取方法 | |
| CN114371505A (zh) | 一种基于地震分频技术的多子波反演方法及系统 | |
| CN107272060B (zh) | 一种地层品质因子的提取方法及系统 | |
| CN110794455B (zh) | 一种地震波传播能量衰减补偿方法 | |
| CN112415601A (zh) | 表层品质因子q值的确定方法及装置 | |
| CN113009569A (zh) | 一种地震偏移成像方法及装置 | |
| CN117724165B (zh) | 一种基于时变子波的品质因子估计方法 | |
| CN108303739A (zh) | 一种时频域真振幅恢复的方法和纵向能量补偿的方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |