CN108845357B - 一种基于同步挤压小波变换估计地层等效品质因子的方法 - Google Patents
一种基于同步挤压小波变换估计地层等效品质因子的方法 Download PDFInfo
- Publication number
- CN108845357B CN108845357B CN201810605738.9A CN201810605738A CN108845357B CN 108845357 B CN108845357 B CN 108845357B CN 201810605738 A CN201810605738 A CN 201810605738A CN 108845357 B CN108845357 B CN 108845357B
- Authority
- CN
- China
- Prior art keywords
- seismic
- quality factor
- wavelet transform
- time
- frequency
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 57
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 31
- 238000001125 extrusion Methods 0.000 title claims abstract description 19
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 17
- 238000001228 spectrum Methods 0.000 claims abstract description 24
- 238000001914 filtration Methods 0.000 claims abstract description 8
- 239000006185 dispersion Substances 0.000 claims description 4
- 239000000126 substance Substances 0.000 claims description 2
- 230000003595 spectral effect Effects 0.000 abstract description 9
- 230000002776 aggregation Effects 0.000 abstract description 3
- 238000004220 aggregation Methods 0.000 abstract description 3
- 238000000605 extraction Methods 0.000 abstract description 2
- 230000001131 transforming effect Effects 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 230000005012 migration Effects 0.000 description 4
- 238000013508 migration Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000002708 enhancing effect Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000021715 photosynthesis, light harvesting Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000008707 rearrangement Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/20—Trace signal pre-filtering to select, remove or transform specific events or signal components, i.e. trace-in/trace-out
- G01V2210/24—Multi-trace filtering
- G01V2210/242—F-k filtering, e.g. ground roll
-
- 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
本发明公开了一种基于同步挤压小波变换估计地层等效品质因子的方法,该方法是一种基于Q补偿原理的品质因子提取方法。通过采用具有极高时频分辨率和能量聚集性的同步挤压小波变换首先将相邻层位(叠后地震数据)或者相邻地震记录(垂直地震剖面)的地震数据变换到时频域,然后在一定Q值预设范围内进行反Q滤波,通过计算时频域Q补偿后的结果与参考波形频谱之间的均方误差,寻找满足最小均方误差时对应的品质因子作为最佳等效品质因子,实现品质因子的最优化估计。本发明避免了常规谱比法估计品质因子频段选择的难点,也不存在类似质心频率法等关于震源频谱为高斯谱的假设等,拓展了品质因子估计的方法,提高了品质因子的估算精度。
Description
技术领域
本发明涉及油气勘探地球物理处理方法技术领域,具体的说是一种基于同步挤压小波变换估计地层等效品质因子的方法。
背景技术
当地震波在地下传播时由于介质的非均匀性和滞弹性会遭受大地衰减,从而导致高频分量的能量耗散和速度色散。引起地震子波失真和振幅减弱的衰减的这种固有介质性质通常用品质因子Q来表示。Q对岩性、孔隙度、饱和度和渗透率等地质信息很敏感,Q估算有助于提高地震资料的分辨率,增强地层特征的精细细节,并用作直接烃类检测指示等。
目前地层品质因子提取方法包括振幅衰减法、上升时间法、谱比法、质心频率偏移法、峰值频率偏移法和解析信号法等。当没有真实可用的振幅值的时候,谱比法、质心频率偏移法和峰值频率偏移法等效果最佳。但是谱比法在计算的时候需要选择合适的频带范围,不同的频带范围结果差异很大,同时对噪声很敏感;峰值频率偏移法计算前提是假设地震子波是雷克子波,同时存在识别峰值频率时易受噪声影响的问题;质心频率偏移法是在假设震源频谱为高斯谱的前提下得出的,而实际地震资料子波谱通常是非高斯分布,从而会导致计算误差,这种方法同样在计算质心频率时候容易受到噪声影响。
发明内容
针对现有技术中存在的上述不足之处,本发明要解决的技术问题是提供一种基于同步挤压小波变换估计地层等效品质因子的方法。新的鲁棒性的品质因子高精度估计方法,采用的是基于同步挤压小波变换Q补偿的品质因子Q估计算法,进而提高品质因子Q的估算精度
本发明为实现上述目的所采用的技术方案是:一种基于同步挤压小波变换估计地层等效品质因子的方法,包括以下步骤:
对于叠后地震数据,综合利用地质信息、测井以及合成地震记录标定目标层;对于垂直地震剖面记录,利用F-K域滤波法分离下行波,在波峰点拾取每道地震记录的初至时间;
计算叠后地震数据的相邻层位处,或者垂直地震剖面记录的相邻地震道的地震数据在Q值预设范围内进行Q补偿后的结果;
计算时频域Q补偿后的结果与参考波形频谱之间的均方误差,以均方误差最小值对应的Q值作为对应层位/地震道处的最佳品质因子Q。
所述计算叠后地震数据的相邻层位处,或者垂直地震剖面记录的相邻地震道的地震数据在设定范围内进行Q补偿后的结果,包括以下步骤:
对地震数据进行同步挤压小波变换;
在时频域完成对各个相邻层位处/相邻地震道的振幅补偿;
分别计算每个Q值对应的Q补偿后的矫正信号。
所述对地震数据进行同步挤压小波变换,具体为:
其中,△ω=ωl-ωl-1.△ak=ak-1-ak.ωl是[ωl-△ω/2,ωl+△ω/2]范围内的中心频率,Ws(ak,b)是地震数据的连续小波变换,ak是第k个尺度因子,a是尺度因子,b是时移因子,Ts(ωl,b)是地震数据的同步挤压小波变换;
尺度因子a与角频率ω具有下列关系:
a=ω0/ω (4)
其中,ω0是Morlet小波的中心角频率。
所述连续小波变换采用Morlet小波,所述连续小波变换Ws(a,b)可表示为:
其中,ψ*是Morlet小波的复共轭函数,a表示尺度因子,b是时移因子,t是时间,s(t)是地震数据。
所述在时频域完成对各个相邻层位处/相邻地震道的振幅补偿,计算每个Q值对应的Q补偿后的矫正信号,具体为:
其中,ωl是[ωl-△ω/2,ωl+△ω/2]范围内的中心频率,b是时移因子,Q是品质因子,ωr是参考角频率。
所述最佳品质因子Q为:
其中,min(·)表示取结果的最小值,N为Q的预设搜索范围,m表示计算的地震记录的长度;Sr是参考层位点处的频谱或参考地震记录处的频谱,且ωl是[ωl-△ω/2,ωl+△ω/2]范围内的中心频率,br是参考频谱的时移因子,Sout是时频域Q补偿后的频谱结果且
本发明具有以下优点及有益效果:
1.本发明使用了具有极高时频分辨率和能量聚集性的同步挤压小波变换,更适合非线性非平稳地震信号的处理,保证了计算结果的更精确。
2.本发明推导了同步挤压小波变换域振幅补偿公式。发展了一种高精度同步挤压小波变换域振幅补偿方法,提高了计算结果的精确性。
3.本发明基于同步挤压小波变换Q补偿的品质因子估计方法,避免了传统品质因子估计方法对噪声敏感,需要选择频段以及估计方法的各种假设前提如子波为雷克子波等问题。
4.本发明算法运行速度快,适合大批量地震信号处理。
附图说明
图1为利用地层品质因子求取的VSP合成记录;
图2为利用本发明方法及传统谱比法提取的VSP合成记录的地层等效品质因子Q与真实值的对比图;
图3为加噪的VSP合成记录;
图4为利用本发明方法及传统谱比法提取的加噪的VSP合成记录的地层等效品质因子Q与真实值的对比图;
图5为利用本发明方法计算的叠后过井地震道的地层等效品质因子Q;
图6为本发明方法流程图。
具体实施方式
下面结合附图及实施例对本发明做进一步的详细说明。
本发明的基于同步挤压小波变换的Q补偿方法估计地层等效品质因子技术,其核心问题是进行高精度时频域Q补偿,通过利用具有高时频分辨率和能量聚集性的同步挤压小波变换,实现时频域的Q补偿(同步挤压小波变换是一种将频率重排方案与连续小波变换相结合的高分辨率时频分析方法)。然后在一定预设范围内计算Q补偿后的信号与参考信号的均方误差,在时频域通过在一定Q值范围内搜索经过反Q滤波重构的地震信号与参考地震信号均方误差最小点,实现品质因子Q估计的最优化输出。
如图6所示,基于同步挤压小波变换的地层等效品质因子估算方法,包括以下步骤:
1.对于叠后地震数据,综合利用地质信息、测井以及合成地震记录等资料准确标定层位;对于垂直地震剖面记录,利用F-K域滤波等方法准确分离下行波,拾取每道地震记录的初至时间。
2.利用基于同步挤压小波变换的Q补偿方法计算相邻层位(叠后地震数据)或者相邻地震记录(垂直地震剖面数据)在一定设定范围内进行Q补偿后的结果。在时频域完成对各个层位点(叠后地震数据)或相邻地震记录(垂直地震剖面数据)处的振幅补偿。
2.1对地震数据(叠后地震数据或垂直地震剖面数据)进行同步挤压小波变换
利用下式对地震道s(t)进行同步挤压小波变换:
其中,△ω=ωl-ωl-1.△ak=ak-1-ak.ωl是[ωl-△ω/2,ωl+△ω/2]范围内的中心频率。Ws(a,b)是地震数据s(t)的连续小波变换,ak是第k个尺度因子,a是尺度因子,b是时移因子,Ts(ωl,b)是地震数据的同步挤压小波变换。这里,我们选择Morlet小波计算连续小波变换。Morlet小波具有以下形式:
其中,ω0是Morlet小波的中心角频率。从而连续小波变换Ws(a,b)可表示为:
其中,ψ*是Morlet小波的复共轭函数,a表示尺度因子,t是时间,b是时移因子。尺度因子a与角频率ω具有下列关系:
a=ω0/ω (4)
其中,ω0是Morlet小波的中心角频率。
2.2振幅补偿
其中,i是虚数单位,△r是传播距离。波数k(ω)与品质因子Q具有下列关系:
其中v是相速度,可表示为:
ωr是参考角频率。
综合以上(5)-(7)式,将传播距离△r用传播时间t=△r/vr(vr是参考相速度)代替,可以得到下式:
由式(8),我们可以得到一个反Q滤波算子G:
在同步挤压小波变换域,令
则Λ(ωl,b)是与能量衰减有关的幅度补偿。θ(ωl,b)是与速度色散有关的相位校正,ωr是参考角频率。如式(10)-(11)所示,作为时间和频率的指数函数的幅度补偿是一个不稳定过程,其在增强带宽的同时放大了环境噪声的幅度。而相位校正是时间和频率的振荡函数,这是一个无条件稳定的过程。
2.3在预设Q范围内,利用公式(12)分别计算每个Q值对应的Q补偿后的矫正信号。
3.计算时频域Q补偿后的结果与参考波形频谱之间的均方误差,其中,均方误差最小值对应的Q值为对应层位点(叠后地震数据)或记录(垂直地震剖面数据)处的品质因子Q。
计算时频域Q补偿后的频谱结果与参考点处频谱之间的均方误差。设参考层位点(叠后地震数据)处的频谱或参考地震记录(垂直地震剖面)处的频谱为时频域Q补偿后的频谱结果为其中N为Q的预设搜索范围。则对应计算层位点(叠后地震数据)或记录处的品质因子Q(垂直地震剖面)的最佳品质因子Q为
其中,min(·)表示取结果的最小值,m表示计算的地震记录的长度,m表示计算的地震记录的长度,br是参考频谱的时移因子,b是计算幅度补偿处地震记录频谱的时移因子。
本发明具有以下实施例:
(1)如图1所示,利用雷克子波生成30道VSP合成记录。该合成记录不含噪声。
(2)图2为利用本技术及传统谱比法提取的VSP合成记录的地层等效品质因子Q与真实值对比。从图中可以看出,对于无噪声的VSP合成记录,本技术和传统谱比法提取到的品质因子Q与真实理论值都相近,误差很小。相对而言,本技术提取到的品质因子Q更接近理论值。
(3)图3为加噪的VSP合成记录。对图1所得的VSP合成记录,添加了最大幅度为0.001的随机噪声。
(4)图4为利用本技术及传统谱比法提取的加噪VSP记录的地层等效品质因子Q与真实值对比。从图中可以看到,本发明技术提取到的品质因子Q更接近理论真实值,误差较小,而传统谱比法对噪声比较敏感,误差较大。
(5)利用本技术计算的叠后过井地震道的地层等效品质因子Q。这里,利用川西地区一条叠后过井地震道数据进行分析。图5(a)是该过井地震道。图5(b)是该过井地震道对应的基于同步挤压小波变换的时频图。图5(c)是利用本发明技术提取的地层品质因子Q。在已知含气储层处(图中虚线矩形框所示),可以看到对应的时频能量较强,提取到的品质因子相对最小。符合实际情况。
Claims (4)
1.一种基于同步挤压小波变换估计地层等效品质因子的方法,其特征在于,包括以下步骤:
对于叠后地震数据,综合利用地质信息、测井以及合成地震记录标定目标层;对于垂直地震剖面记录,利用F-K域滤波法分离下行波,在波峰点拾取每道地震记录的初至时间;
计算叠后地震数据的相邻层位处,或者垂直地震剖面记录的相邻地震道的地震数据在Q值预设范围内进行Q补偿后的结果;
计算时频域Q补偿后的结果与参考波形频谱之间的均方误差,以均方误差最小值对应的Q值作为对应层位/地震道处的最佳品质因子Q;
所述计算叠后地震数据的相邻层位处,或者垂直地震剖面记录的相邻地震道的地震数据在Q值预设范围内进行Q补偿后的结果,包括以下步骤:
对地震数据进行同步挤压小波变换;
在时频域完成对各个相邻层位处/相邻地震道的振幅补偿;
分别计算每个Q值对应的Q补偿后的矫正信号;
所述在时频域完成对各个相邻层位处/相邻地震道的振幅补偿,计算每个Q值对应的Q补偿后的矫正信号,具体为:
其中,ωl是[ωl-Δω/2,ωl+Δω/2]范围内的中心频率,b是时移因子,Q是品质因子,ωr是参考角频率,Δω=ωl-ωl-1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810605738.9A CN108845357B (zh) | 2018-06-13 | 2018-06-13 | 一种基于同步挤压小波变换估计地层等效品质因子的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810605738.9A CN108845357B (zh) | 2018-06-13 | 2018-06-13 | 一种基于同步挤压小波变换估计地层等效品质因子的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108845357A CN108845357A (zh) | 2018-11-20 |
CN108845357B true CN108845357B (zh) | 2020-12-22 |
Family
ID=64211125
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810605738.9A Expired - Fee Related CN108845357B (zh) | 2018-06-13 | 2018-06-13 | 一种基于同步挤压小波变换估计地层等效品质因子的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108845357B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111337977B (zh) * | 2018-12-19 | 2022-03-29 | 中国石油天然气股份有限公司 | 品质因子的确定方法及装置 |
CN110187388B (zh) * | 2019-06-06 | 2021-01-05 | 成都信息工程大学 | 一种基于变分模态分解的稳定地震品质因子q估计方法 |
CN112394402A (zh) * | 2019-08-19 | 2021-02-23 | 中国石油化工股份有限公司 | 基于同步挤压小波变换检测微地震信号的方法和系统 |
CN112099083B (zh) * | 2020-08-26 | 2023-10-13 | 中化地质矿山总局地质研究院 | 一种基于双谱谱比对数的品质因子估计方法及系统 |
CN112904415B (zh) * | 2021-01-20 | 2022-03-15 | 中国石油大学(北京) | 基于重构源谱一致性的q层析反演方法、系统及处理设备 |
CN114035225B (zh) * | 2021-10-19 | 2023-08-15 | 成都信息工程大学 | 一种地震资料高分辨率处理bpfe方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103376464A (zh) * | 2012-04-13 | 2013-10-30 | 中国石油天然气集团公司 | 一种地层品质因子反演方法 |
CN107390267A (zh) * | 2017-07-27 | 2017-11-24 | 西安交通大学 | 一种同步挤压变换域的地震资料衰减补偿方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7376517B2 (en) * | 2005-05-13 | 2008-05-20 | Chevron U.S.A. Inc. | Method for estimation of interval seismic quality factor |
US7830745B2 (en) * | 2007-12-27 | 2010-11-09 | Schlumberger Technology Corporation | Identifying the Q-factor using microseismic event generated S-coda waves |
WO2014191011A1 (en) * | 2013-05-27 | 2014-12-04 | Statoil Petroleum As | High resolution estimation of attenuation from vertical seismic profiles |
-
2018
- 2018-06-13 CN CN201810605738.9A patent/CN108845357B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103376464A (zh) * | 2012-04-13 | 2013-10-30 | 中国石油天然气集团公司 | 一种地层品质因子反演方法 |
CN107390267A (zh) * | 2017-07-27 | 2017-11-24 | 西安交通大学 | 一种同步挤压变换域的地震资料衰减补偿方法 |
Non-Patent Citations (2)
Title |
---|
"利用同步挤压变换检测微地震信号";刘晗 等;《中国科技论文》;20151130;第10卷(第21期);第2472-2476页 * |
"基于等效Q值的时域反Q滤波方法";吕喜滨;《中国优秀硕士学位论文全文数据库-基础科学辑》;20120115(第1期);正文第11-19页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108845357A (zh) | 2018-11-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108845357B (zh) | 一种基于同步挤压小波变换估计地层等效品质因子的方法 | |
CN106597532B (zh) | 一种结合井资料与层位资料的叠前地震数据频带拓展方法 | |
RU2579164C1 (ru) | Способ обращения для определения добротности геологической среды | |
CN104849756B (zh) | 一种提高地震数据分辨率增强有效弱信号能量的方法 | |
CN109061764B (zh) | 一种分频融合波阻抗反演方法 | |
CN109669212B (zh) | 地震数据处理方法、地层品质因子估算方法与装置 | |
CN110596758B (zh) | 一种地震信号低频能量补偿方法 | |
CN110187388B (zh) | 一种基于变分模态分解的稳定地震品质因子q估计方法 | |
CN106019376B (zh) | 一种频率驱动空变q值模型构建的地震波补偿方法 | |
CN110703331A (zh) | 一种基于常q粘滞声波方程的衰减补偿逆时偏移实现方法 | |
CN112305612B (zh) | 高分辨率复谱分解时频空间域振幅随偏移距变化校正方法 | |
CN110879410A (zh) | 一种多分量地震面波勘探方法 | |
CN109143345B (zh) | 基于模拟退火的品质因子q非线性反演方法及系统 | |
CN112578438B (zh) | 一种地震子波提取方法及系统 | |
CN110749923A (zh) | 一种基于范数方程提高分辨率的反褶积方法 | |
CN114371505A (zh) | 一种基于地震分频技术的多子波反演方法及系统 | |
CN110673211B (zh) | 一种基于测井与地震数据的品质因子建模方法 | |
CN110568491B (zh) | 一种品质因子q的估算方法 | |
CN111665536B (zh) | 基于微测井子波定量化约束的井深设计方法 | |
CN110941028B (zh) | 一种对碳酸盐岩溶蚀孔洞型地热能储层位置进行定位的方法和系统 | |
CN106772571B (zh) | 一种提高相同地区不同震源叠前地震数据精度的方法 | |
CN102778692B (zh) | 一种利用地震数据确定品质因数的方法 | |
Shen et al. | Partial frequency band match filtering based on high-sensitivity data: method and applications | |
WO2021147841A1 (zh) | 利用垂直地震数据提取下行子波和衰减参数的方法、装置 | |
CN114137608B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20201222 |