CN107356964A - S变换域基于变分原理的q值估计与补偿方法 - Google Patents
S变换域基于变分原理的q值估计与补偿方法 Download PDFInfo
- Publication number
- CN107356964A CN107356964A CN201710544007.3A CN201710544007A CN107356964A CN 107356964 A CN107356964 A CN 107356964A CN 201710544007 A CN201710544007 A CN 201710544007A CN 107356964 A CN107356964 A CN 107356964A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- transformation
- msubsup
- msup
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000001228 spectrum Methods 0.000 claims description 33
- 238000012545 processing Methods 0.000 claims description 11
- 230000009466 transformation Effects 0.000 claims description 9
- 238000013461 design Methods 0.000 claims description 8
- 230000000694 effects Effects 0.000 claims description 7
- 239000013598 vector Substances 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000010521 absorption reaction Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 230000008030 elimination Effects 0.000 claims description 3
- 238000003379 elimination reaction Methods 0.000 claims description 3
- 238000011084 recovery Methods 0.000 claims description 3
- 230000001808 coupling effect Effects 0.000 abstract description 3
- 230000008569 process Effects 0.000 abstract description 3
- 230000008859 change Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000009331 sowing Methods 0.000 description 1
- 238000010183 spectrum analysis 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/301—Analysis for determining seismic cross-sections or geostructures
-
- 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
- G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
-
- 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/22—Time-domain filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/40—Transforming data representation
- G01V2210/48—Other transforms
-
- 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/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/612—Previously recorded data, e.g. time-lapse or 4D
-
- 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/63—Seismic attributes, e.g. amplitude, polarity, instant phase
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
本发明公开一种S变换域基于变分原理的Q值估计与补偿方法,将S变换引入非平稳地震记录褶积模型,首先将地震记录进行S变换,相比于时间域以及频率域方法,在时频域地震记录的非平稳特性可得更加清晰的刻画,基于变分原理进行Q值估计,不需要对地震记录进行加窗处理,可以减小窗函数对Q值估计结果的影响,此外,由于是对整个地震记录进行整体处理,可以避免在处理反射地震数据时吸收与散射的耦合效应,抗噪性能更强,反演的参数精度更高。
Description
【技术领域】
本发明属于地球物理勘探领域,涉及一种地球模型物性参数反演方法,特别涉及一种S变换域基于变分原理的Q值估计与补偿方法。
【背景技术】
地震波在地下介质中传播时会逐渐衰减,造成地震波衰减的主要因素包括几何扩散、介质的非完全弹性以及透射/反射损失等。其中由地下介质的非完全弹性引起的衰减是介质的固有特性,这种特性可以用品质因子Q来进行定量描述。非完全弹性引起的衰减会不仅导致地震波的振幅衰减,而且会造成相位畸变,导致地震资料的纵向分辨率降低。如果已知Q值,对地震数据进行反Q滤波,则可以使地震信号的高频能量得到增强,从而提高地震资料的纵向分辨率;另外,Q值也与岩性、饱和度、孔隙率等参数密切相关,可用于储层识别和烃类检测。因此,反演介质的Q值,然后进行能量补偿在地震勘探领域具有重要的实际意义。
国内外学者提出了多种Q值估计的方法。这些方法大致可以归为两类,一类为时间域方法,另一类为频率域方法。时间域方法主要包括子波模拟法、脉冲上升时间法和解析信号法等,这些方法都依赖于高保真度的地震数据,然而实际接收的地震记录由于受到几何扩散,反/透射损失等因素的影响,致使时间域Q值估计方法精度降低。频率域常用的方法包括对数谱比值法、中心频率偏移法、子波包络峰值频率偏移法、瞬时频率匹配法等,在处理反射地震数据时,这些方法首先要用一个时窗截取其中一段地震记录,然后对截取出来的的地震记录进行频谱分析,进而估算出Q值。在这类方法中,采用不同时窗对Q值估计结果影响较大,如何最佳地选择时窗的类型及长度一直是一个难题。Tonn对近十种Q值估计方法进行了比较,结果表明每种方法都有一定的适用条件。
Margrave和Lamoureux在2011年提出了一种Gabor反褶积方法,这是一种在时频域进行地震资料高分辨处理的方法,该算法将Wiener反褶积方法扩展到非平稳地震记录的情况。该方法不需要预先进行估计Q值就可以在时频域利用平滑函数进行反Q滤波,但是这种处理方法对平滑窗依赖性很强,若窗函数选取不恰当,处理效果会比较差,而且经这种方法处理后的地震记录,不能很好地反映反射系数的局部相对关系,直观地讲会产生等幅效应,类似于经自动增益控制(AGC)处理后的结果,在某些情况下不利于地震资料解释。
【发明内容】
针对现有技术的不足,本发明目的在于提供一种S变换域基于变分原理的Q值估计与补偿方法,利用反射地震资料,基于地下介质单Q值假设,在S变换域首先利用变分法进行Q值估计,然后构造一种自适应滤波器对Q值进行补偿,从而得到一种高保真度、高分辨率的地震剖面。
为达到上述目的,本发明是采取如下技术方案予以实现的:
S变换域基于变分原理的Q值估计与补偿方法,其特征在于:基于非平稳地震记录褶积模型,对叠后反射地震资料采用Q值与频率无关的假设,提出S变换域基于变分法的Q值估计与补偿方法。
进一步,具体包括以下步骤:
1)采集原始地震记录,对采集到的地震数据进行噪音压制、校正、真振幅恢复、偏移处理,处理后得到叠后反射地震记录;
2)对反射地震记录x(t)进行S变换,将地震记录从一维时间信号变到二维时频平面,S变换结果记为ST(t,f);
其中,x(t)为待分析单道信号;ST(t,f)为时频变换结果;g(t-τ,f)为随频率变化的窗函数;τ为待分析信号的时间位置;t为S变换结果的时间变量,f为频率变量;i为虚数单位;
3)在地震记录的浅层即前四分之一段,利用多项式拟合估计出参考子波;
4)在S变换域利用变分法估计Q值;
5)在S变换域设计滤波器进行Q补偿;
6)对补偿以后的S变换进行逆变换,得到补偿以后的地震记录。
进一步,所述步骤3)具体实施方案为,对浅层地震记录做傅里叶变换得到振幅谱w(f),基于地震子波振幅谱为单峰平滑函数的假设,利用以下多项式对振幅谱进行拟合,估计出参考子波;
其中K=2,N=5,将ws(f)作为参考子波的振幅谱。
进一步,所述步骤4)在S变换域利用变分法估计Q值具体如下:
经证明得到|STx(t,f)|=|ws(f)||α(t,f)||STr(t,f)| (4)
对公式(4)两边取对数
ln|STx(τ,v)|=ln|ws(f)|-πft/Q (5)
在最小二乘意义下定义目标函数
在实际应用中,由于地震数据是频带有限的,为了减少计算量和数值稳定性,需要确定积分区域,为此定义如下加权函数,
t0和t1为地震记录的起始时刻和终点时刻,f0和f1为选取的频带区间;对于Q,α(ws,Q)为关于Q的一元函数;对于ws(f),α(ws,Q)为ws(f)的泛函,为了寻找极小值,对其分别求变分可得
其中
进一步,所述步骤5)根据步骤4)计算的Q值设计滤波器进行Q补偿;
在S变换的每个点谱ST(tk,f)估计tk时刻的子波Lk(f),在每个时刻tk,设计如下形式的补偿滤波器
αk(f)=exp(βkf) (12)
其中βk为与介质粘弹性吸收衰减有关的补偿参数,补偿以后的tk时刻的点谱为
STc(tk,f)=Lk(f)αk(f) (13)
为了确定补偿参数βk,定义如下的相关系数
其中< >表示向量内积运算,|| ||表示向量的2-范数;Ck(σ)可以用来度量tk时刻补偿以后的点谱STc(tk,f)与参考子波振幅谱的相似性,Ck(σ)值越大,相似性越好,即补偿效果越好;当Ck(σ)的值达到最大时,对应的补偿参数为σk,即
σk即为tk时刻对应的补偿参数。
进一步,所述步骤6)依次对每个时刻的点谱进行Q补偿后,得到补偿以后的S变换STc(t,f),对此进行逆S变换
得到补偿以后的高分辨率地震记录。
与现有技术相比,本发明具有以下有益效果:
本发明的S变换域基于变分原理的Q值估计与补偿方法,将S变换引入非平稳地震记录褶积模型,首先将地震记录进行S变换,相比于时间域以及频率域方法,在时频域地震记录的非平稳特性可得更加清晰的刻画,基于变分原理进行Q值估计,不需要对地震记录进行加窗处理,可以减小窗函数对Q值估计结果的影响,此外,由于是对整个地震记录进行整体处理,可以避免在处理反射地震数据时吸收与散射的耦合效应,抗噪性能更强,反演的参数精度更高。
用于反射地震资料的高分辨处理,该算法首先将地震资料变换到时频域,从而可以有效刻画地震记录的非平稳特性,通过在时频域构造目标函数,进而利用变分原理对目标函数求解Q值。与常用的时间域或者频率域估计Q值算法相比,时频域计算Q可以有效克服吸收与散射的耦合效应,计算结果不受窗函数影响;在补偿过程中,针对记录各段不同的特性给出不同的补偿参数,从而使得记录的各段都能得到相对恰当地补偿。
【附图说明】
图1为合成地震记录的高分辨率处理结果;
(a)震源子波;(b)地层反射系数序列;(c)无衰减的合成地震记录;(d)含衰减(Q=50)的合成地震记录;(e)Q补偿以后的高分辨率地震记录;
图2为含衰减的合成地震记录的S变换结果;
图3为在浅层提取的参考子波的振幅谱;
图4为补偿以后的地震记录的S变换结果;
图5为实际地震数据的高分辨率处理结果
(a)某油田三维实际数据的一段剖面;(b)利用本发明处理以后的高分辨率地震剖面。
【具体实施方式】
下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,本文所描述的实施例仅仅为本发明的一部分实施例,而非全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有创造性劳动前提下所获得的所有其他实施例,均属于本发明保护范围。
下面结合附图对本发明做进一步详细描述:
本发明基于非平稳地震记录褶积模型,采用Q值与频率无关的假设,提出S变换域基于变分法的Q值估计与补偿方法。
本发明的物质基础是叠后反射地震资料。
本发明S变换域基于变分法的Q值估计与补偿方法具体步骤分别为:
1)采集地震数据,然后对采集到的地震数据做常规预处理,如噪音压制、校正、真振幅恢复、偏移等。处理后得到叠后反射地震记录。
2)对反射地震记录x(t)进行S变换,将地震记录从一维时间信号变到二维时频平面,S变换结果记为ST(t,f);
其中,x(t)为待分析单道信号;ST(t,f)为时频变换结果;g(t-τ,f)为随频率变化的窗函数;τ为待分析信号的时间位置;t为S变换结果的时间变量,f为频率变量;i为虚数单位。
ST(t,f)在频率域的等价形式为
其中,X(α)表示信号x(t)的傅里叶变换。通过公式(2),ST(t,f)可以在频率域利用FFT进行快速计算。
3)把在地震记录的浅层(地震记录的前四分之一段)估计出的子波作为参考子波。
具体实施方案为,对浅层地震记录做傅里叶变换得到振幅谱w(f),基于地震子波振幅谱为单峰平滑函数的假设,利用以下多项式对振幅谱进行拟合
其中K=2,N=5,将ws(f)作为参考子波的振幅谱。
4)S变换域利用变分法计算Q值,首先我们证明了以下结论
|STx(t,f)|=|ws(f)||α(t,f)||STr(t,f)| (4)
其中STx(t,f)为地震记录的S变换结果,ws(f)为步骤3)中提取的参考子波,为与Q值相关的衰减函数,STx(t,f)为反射系数的S变换结果。假设反射系数具有随机白噪的统计特性,则可以认为地震道的时频谱值的细节部分是由反射系数引起,主要趋势由时变子波引起。如果对具有白噪特性的反射系数频谱做平滑,其结果为单位矢量。因此,对公式(4)两边取对数
ln|STx(τ,v)|=ln|ws(f)|-πft/Q (5)
在最小二乘意义下定义目标函数
在实际应用中,由于地震数据是频带有限的,为了减少计算量和数值稳定性,需要确定积分区域,为此定义如下加权函数,
t0和t1为地震记录的起始时刻和终点时刻,f0和f1为选取的频带区间。对于Q,α(ws,Q)为关于Q的一元函数;对于ws(f),α(ws,Q)为ws(f)的泛函,为了寻找极小值,对其分别求变分可得
其中
5)S变换域根据步骤4)计算的Q值设计滤波器进行Q补。在S变换的每个点谱ST(tk,f),采用同样的方法估计tk时刻的子波Lk(f)。在每个时刻tk,设计如下形式的补偿滤波器
αk(f)=exp(βkf) (12)
其中βk为与介质粘弹性吸收衰减有关的补偿参数,补偿以后的tk时刻的点谱为
STc(tk,f)=Lk(f)αk(f)。 (13)
为了确定补偿参数βk,定义如下的相关系数
其中< >表示向量内积运算,|| ||表示向量的2-范数。Ck(σ)可以用来度量tk时刻补偿以后的点谱STc(tk,f)与参考子波振幅谱的相似性,Ck(σ)值越大,相似性越好,即补偿效果越好。当Ck(σ)的值达到最大时,对应的补偿参数为σk,即
σk即为tk时刻对应的补偿参数。
6)依次对每个时刻的点谱进行补偿后,得到补偿以后的S变换STc(t,f)。对此进行逆S变换
得到补偿以后的高分辨率地震记录。
效果分析
一、数值仿真实验
震源选择70Hz Ricker子波,如图1(a)所示。反射系数如图1(b)所示,采用如下非平稳褶积模型合成地震记录
s(t)=∫∫w(f)α(u,f)r(u)exp[2πif(t-u)]dfdu (17)
其中r(u)为反射系数序列,w(f)为震源子波的傅里叶变换,α(u,f)为与Q有关的衰减函数
H为希尔伯特变换。图1(c)为无衰减地震记录(Q=10000),图1(d)为Q=50时的地震记录,图2为地震记录的S变换结果,可以明显看出,由于衰减效应,地震记录的主频随传播时间不断向低频移动,导致地震记录在深层处的分辨率降低。为了对地震记录补偿,首先在S变换域利用变分法估计Q值,本方法Q值的估计结果为Q=52.1。然后在时频域利用滤波器进行补偿,补偿之前,首先需要浅层提取一个参考子波,如图3所示,黑色星号为提取的参考子波的振幅谱。以参考子波为依据,在每个点谱寻找合适的补偿参数进行补偿,图4为补偿以后的S变换结果。最后进行反变换,得到补偿以后的地震记录,如图1(e)所示。从补偿结果来看,地震记录的深层能量得到有效增强,分辨率也显著提高。
二、实际地震资料
为了进一步验证本发明提出的方法的有效性,将该方法用于某油田提供的三维实际地震资料的一段剖面。该剖面包含214道,每道的采样点为351,采样间隔为2ms。图5(a)为原始地震剖面,图5(b)为利用本发明处理后的高分辨率地震剖面。对比图5(a)和图5(b)可以看出用本发明处理后的地震剖面的纵向分辨率明显提高,对地层结构的刻画更加精细,对薄层的识别能力明显增强(如黑色箭头指示区域),并且没有破坏剖面原有的结构特征。
以上实施例仅用于说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细说明,所属领域的普通技术人员依然可以对本发明的具体实施方案进行修改或者等同替换,而这些并未脱离本发明精神和范围的任何修改或者等同替换,其均在本发明的权利要求保护范围之内。
Claims (6)
1.S变换域基于变分原理的Q值估计与补偿方法,其特征在于:基于非平稳地震记录褶积模型,对叠后反射地震资料采用Q值与频率无关的假设,提出S变换域基于变分法的Q值估计与补偿方法。
2.根据权利要求1所述的S变换域基于变分原理的Q值估计与补偿方法,其特征在于具体包括以下步骤:
1)采集原始地震记录,对采集到的地震数据进行噪音压制、校正、真振幅恢复、偏移处理,处理后得到叠后反射地震记录;
2)对反射地震记录x(t)进行S变换,将地震记录从一维时间信号变到二维时频平面,S变换结果记为ST(t,f);
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>S</mi>
<mi>T</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>f</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>&infin;</mi>
</mrow>
<mi>&infin;</mi>
</msubsup>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>&tau;</mi>
<mo>-</mo>
<mi>t</mi>
<mo>,</mo>
<mi>f</mi>
</mrow>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>i</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>f</mi>
<mi>&tau;</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>&tau;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>&infin;</mi>
</mrow>
<mi>&infin;</mi>
</msubsup>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mo>|</mo>
<mi>f</mi>
<mo>|</mo>
</mrow>
<msqrt>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</msqrt>
</mfrac>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mi>&tau;</mi>
<mo>-</mo>
<mi>t</mi>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msup>
<mi>f</mi>
<mn>2</mn>
</msup>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>i</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>f</mi>
<mi>&tau;</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>&tau;</mi>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,x(t)为待分析单道信号;ST(t,f)为时频变换结果;g(t-τ,f)为随频率变化的窗函数;τ为待分析信号的时间位置;t为S变换结果的时间变量,f为频率变量;i为虚数单位;
3)在地震记录的浅层即前四分之一段,利用多项式拟合估计出参考子波;
4)在S变换域利用变分法估计Q值;
5)在S变换域设计滤波器进行Q补偿;
6)对补偿以后的S变换进行逆变换,得到补偿以后的地震记录。
3.根据权利要求2所述的S变换域基于变分原理的Q值估计与补偿方法,其特征在于:所述步骤3)具体实施方案为,对浅层地震记录做傅里叶变换得到振幅谱w(f),基于地震子波振幅谱为单峰平滑函数的假设,利用以下多项式对振幅谱进行拟合,估计出参考子波;
<mrow>
<msub>
<mi>w</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mi>f</mi>
<mi>K</mi>
</msup>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
<msup>
<mi>f</mi>
<mi>n</mi>
</msup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中K=2,N=5,将ws(f)作为参考子波的振幅谱。
4.根据权利要求3所述的S变换域基于变分原理的Q值估计与补偿方法,其特征在于:所述步骤4)在S变换域利用变分法估计Q值具体如下:
经证明得到|STx(t,f)|=|ws(f)||α(t,f)||STr(t,f)| (4)
对公式(4)两边取对数
ln|STx(τ,v)|=ln|ws(f)|-πft/Q (5)
在最小二乘意义下定义目标函数
<mrow>
<mi>min</mi>
<mi>&alpha;</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>w</mi>
<mi>s</mi>
</msub>
<mo>,</mo>
<mi>Q</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>&infin;</mi>
</mrow>
<mrow>
<mo>+</mo>
<mi>&infin;</mi>
</mrow>
</msubsup>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>&infin;</mi>
</mrow>
<mrow>
<mo>+</mo>
<mi>&infin;</mi>
</mrow>
</msubsup>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mi>ln</mi>
<mfrac>
<mrow>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>f</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>w</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>&pi;</mi>
<mi>t</mi>
<mi>f</mi>
</mrow>
<mi>Q</mi>
</mfrac>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
<mi>d</mi>
<mi>t</mi>
<mi>d</mi>
<mi>f</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
在实际应用中,由于地震数据是频带有限的,为了减少计算量和数值稳定性,需要确定积分区域,为此定义如下加权函数,
<mrow>
<msub>
<mi>&chi;</mi>
<mi>&Omega;</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>f</mi>
<mo>)</mo>
<mo>&Element;</mo>
<mi>&Omega;</mi>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>f</mi>
<mo>)</mo>
<mo>&NotElement;</mo>
<mi>&Omega;</mi>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>&Omega;</mi>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>&rsqb;</mo>
<mo>&times;</mo>
<mo>&lsqb;</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
t0和t1为地震记录的起始时刻和终点时刻,f0和f1为选取的频带区间;对于Q,α(ws,Q)为关于Q的一元函数;对于ws(f),α(ws,Q)为ws(f)的泛函,为了寻找极小值,对其分别求变分可得
<mrow>
<msub>
<mi>w</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>exp</mi>
<mo>{</mo>
<mi>l</mi>
<mover>
<mi>n</mi>
<mo>&OverBar;</mo>
</mover>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mrow>
<mi>&pi;</mi>
<mi>f</mi>
</mrow>
<mi>Q</mi>
</mfrac>
<mover>
<mi>&tau;</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>}</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>Q</mi>
<mo>=</mo>
<mi>&pi;</mi>
<mfrac>
<mrow>
<munder>
<mo>&Integral;</mo>
<mi>&Omega;</mi>
</munder>
<msup>
<mi>tf</mi>
<mn>2</mn>
</msup>
<mo>&lsqb;</mo>
<mi>t</mi>
<mo>-</mo>
<mover>
<mi>&tau;</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mi>d</mi>
<mi>t</mi>
<mi>d</mi>
<mi>f</mi>
</mrow>
<mrow>
<munder>
<mo>&Integral;</mo>
<mi>&Omega;</mi>
</munder>
<mi>t</mi>
<mi>f</mi>
<mo>&lsqb;</mo>
<mi>l</mi>
<mover>
<mi>n</mi>
<mo>&OverBar;</mo>
</mover>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>ln</mi>
<mi> </mi>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mi>d</mi>
<mi>t</mi>
<mi>d</mi>
<mi>f</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
其中
<mrow>
<mi>l</mi>
<mover>
<mi>n</mi>
<mo>&OverBar;</mo>
</mover>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
</msubsup>
<mi>l</mi>
<mi>n</mi>
<mo>&lsqb;</mo>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<msub>
<mi>&chi;</mi>
<mi>&Omega;</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
</msubsup>
<msub>
<mi>&chi;</mi>
<mi>&Omega;</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mover>
<mi>&tau;</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
</msubsup>
<msub>
<mi>t&chi;</mi>
<mi>&Omega;</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
</msubsup>
<msub>
<mi>&chi;</mi>
<mi>&Omega;</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>.</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
5.根据权利要求4所述的S变换域基于变分原理的Q值估计与补偿方法,其特征在于:所述步骤5)根据步骤4)计算的Q值设计滤波器进行Q补偿;
在S变换的每个点谱ST(tk,f)估计tk时刻的子波Lk(f),在每个时刻tk,设计如下形式的补偿滤波器
αk(f)=exp(βkf) (12)
其中βk为与介质粘弹性吸收衰减有关的补偿参数,补偿以后的tk时刻的点谱为
STc(tk,f)=Lk(f)αk(f) (13)
为了确定补偿参数βk,定义如下的相关系数
<mrow>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&sigma;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo><</mo>
<msub>
<mi>ST</mi>
<mi>c</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>w</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>></mo>
</mrow>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>ST</mi>
<mi>c</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>w</mi>
<mi>s</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
其中<>表示向量内积运算,|| ||表示向量的2-范数;Ck(σ)可以用来度量tk时刻补偿以后的点谱STc(tk,f)与参考子波振幅谱的相似性,Ck(σ)值越大,相似性越好,即补偿效果越好;当Ck(σ)的值达到最大时,对应的补偿参数为σk,即
<mrow>
<msub>
<mi>&sigma;</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mi>arg</mi>
<munder>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mi>&sigma;</mi>
</munder>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&sigma;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mrow>
σk即为tk时刻对应的补偿参数。
6.根据权利要求5所述的S变换域基于变分原理的Q值估计与补偿方法,其特征在于:所述步骤6)依次对每个时刻的点谱进行Q补偿后,得到补偿以后的S变换STc(t,f),对此进行逆S变换
<mrow>
<msub>
<mi>x</mi>
<mi>c</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>&infin;</mi>
</mrow>
<mrow>
<mo>+</mo>
<mi>&infin;</mi>
</mrow>
</msubsup>
<mo>{</mo>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>&infin;</mi>
</mrow>
<mrow>
<mo>+</mo>
<mi>&infin;</mi>
</mrow>
</msubsup>
<msub>
<mi>ST</mi>
<mi>c</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>&tau;</mi>
<mo>}</mo>
<msup>
<mi>e</mi>
<mrow>
<mi>i</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>f</mi>
<mi>t</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>f</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
得到补偿以后的高分辨率地震记录。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710544007.3A CN107356964B (zh) | 2017-07-05 | 2017-07-05 | S变换域基于变分原理的q值估计与补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710544007.3A CN107356964B (zh) | 2017-07-05 | 2017-07-05 | S变换域基于变分原理的q值估计与补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107356964A true CN107356964A (zh) | 2017-11-17 |
CN107356964B CN107356964B (zh) | 2018-10-30 |
Family
ID=60292925
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710544007.3A Active CN107356964B (zh) | 2017-07-05 | 2017-07-05 | S变换域基于变分原理的q值估计与补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107356964B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108693555A (zh) * | 2018-05-16 | 2018-10-23 | 中国石油大学(北京) | 智能化时变盲反褶积宽频处理方法及装置 |
CN108919349A (zh) * | 2018-09-25 | 2018-11-30 | 中国矿业大学(北京) | 高精度反射系数求取方法和系统 |
CN108957540A (zh) * | 2018-07-05 | 2018-12-07 | 河海大学 | 一种高效提取复杂储层中衰减品质因子的方法 |
WO2019164606A1 (en) * | 2018-02-23 | 2019-08-29 | Exxonmobil Research And Engineering Company | Method of extracting intrinsic attenuation from seismic data |
CN112666609A (zh) * | 2020-11-30 | 2021-04-16 | 中国石油大学(华东) | 基于稀疏包络的地震数据低频补偿方法及装置 |
CN113740909A (zh) * | 2021-07-28 | 2021-12-03 | 西安交通大学 | 一种基于稀疏s变换和自适应对数谱比法的地震衰减估计方法、系统、设备及存储介质 |
EP4095566A4 (en) * | 2020-01-20 | 2023-06-21 | China National Petroleum Corporation | METHOD AND APPARATUS FOR EXTRACTING DOWNWARD WAVELET AND ATTENUATION PARAMETERS USING VERTICAL SEISMIC DATA |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150168573A1 (en) * | 2012-04-13 | 2015-06-18 | China National Petroleum Corporation | Geologic quality factor inversion method |
CN104749621A (zh) * | 2013-12-26 | 2015-07-01 | 中国石油化工股份有限公司 | 基于改进s变换的相对保幅点谱模拟高分辨率处理方法 |
CN106291693A (zh) * | 2015-05-21 | 2017-01-04 | 中国石油化工股份有限公司 | 一种基于广义s变换的叠前q值反演方法及系统 |
CN106353795A (zh) * | 2015-07-17 | 2017-01-25 | 中国石油化工股份有限公司 | 一种高精度的地层衰减参数计算方法 |
CN106814393A (zh) * | 2015-11-27 | 2017-06-09 | 中国石油化工股份有限公司 | 一种地层品质因子q的估算方法 |
-
2017
- 2017-07-05 CN CN201710544007.3A patent/CN107356964B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150168573A1 (en) * | 2012-04-13 | 2015-06-18 | China National Petroleum Corporation | Geologic quality factor inversion method |
CN104749621A (zh) * | 2013-12-26 | 2015-07-01 | 中国石油化工股份有限公司 | 基于改进s变换的相对保幅点谱模拟高分辨率处理方法 |
CN106291693A (zh) * | 2015-05-21 | 2017-01-04 | 中国石油化工股份有限公司 | 一种基于广义s变换的叠前q值反演方法及系统 |
CN106353795A (zh) * | 2015-07-17 | 2017-01-25 | 中国石油化工股份有限公司 | 一种高精度的地层衰减参数计算方法 |
CN106814393A (zh) * | 2015-11-27 | 2017-06-09 | 中国石油化工股份有限公司 | 一种地层品质因子q的估算方法 |
Non-Patent Citations (4)
Title |
---|
刘喜武等: "基于广义S变换的吸收衰减补偿方法", 《石油物探》 * |
叶秋焱等: "一种基于叠前数据的综合Q值求取及应用", 《石油地球物理勘探》 * |
王元君等: "时频域动态反褶积方法研究", 《西南石油大学学报(自然科学版)》 * |
魏文等: "中心频率法估算地层吸收参数", 《石油地球物理勘探》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019164606A1 (en) * | 2018-02-23 | 2019-08-29 | Exxonmobil Research And Engineering Company | Method of extracting intrinsic attenuation from seismic data |
CN108693555A (zh) * | 2018-05-16 | 2018-10-23 | 中国石油大学(北京) | 智能化时变盲反褶积宽频处理方法及装置 |
CN108693555B (zh) * | 2018-05-16 | 2019-08-02 | 中国石油大学(北京) | 智能化时变盲反褶积宽频处理方法及装置 |
CN108957540A (zh) * | 2018-07-05 | 2018-12-07 | 河海大学 | 一种高效提取复杂储层中衰减品质因子的方法 |
CN108957540B (zh) * | 2018-07-05 | 2020-02-21 | 河海大学 | 一种高效提取复杂储层中衰减品质因子的方法 |
CN108919349A (zh) * | 2018-09-25 | 2018-11-30 | 中国矿业大学(北京) | 高精度反射系数求取方法和系统 |
CN108919349B (zh) * | 2018-09-25 | 2019-10-18 | 中国矿业大学(北京) | 高精度反射系数求取方法和系统 |
EP4095566A4 (en) * | 2020-01-20 | 2023-06-21 | China National Petroleum Corporation | METHOD AND APPARATUS FOR EXTRACTING DOWNWARD WAVELET AND ATTENUATION PARAMETERS USING VERTICAL SEISMIC DATA |
US11880012B2 (en) | 2020-01-20 | 2024-01-23 | China National Petroleum Corporation | Method and apparatus for extracting downgoing wavelet and attenuation parameters by using vertical seismic data |
CN112666609A (zh) * | 2020-11-30 | 2021-04-16 | 中国石油大学(华东) | 基于稀疏包络的地震数据低频补偿方法及装置 |
CN112666609B (zh) * | 2020-11-30 | 2022-06-17 | 中国石油大学(华东) | 基于稀疏包络的地震数据低频补偿方法及装置 |
CN113740909A (zh) * | 2021-07-28 | 2021-12-03 | 西安交通大学 | 一种基于稀疏s变换和自适应对数谱比法的地震衰减估计方法、系统、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN107356964B (zh) | 2018-10-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107356964B (zh) | S变换域基于变分原理的q值估计与补偿方法 | |
CN105116442B (zh) | 岩性油气藏弱反射地震信号的重构方法 | |
CN107272062B (zh) | 一种数据驱动的地下介质q场估计方法 | |
CN105093294B (zh) | 基于可变模态分解的地震波衰减梯度估计方法 | |
Tian et al. | Variable-eccentricity hyperbolic-trace TFPF for seismic random noise attenuation | |
CN107247290B (zh) | 一种基于时空分数阶滤波的地震资料噪声压制方法 | |
CN103364832A (zh) | 一种基于自适应最优核时频分布的地震衰减定性估计方法 | |
CN103364835A (zh) | 一种地层结构自适应中值滤波方法 | |
CN104181587B (zh) | 一种地震数据振幅谱的相干值获取方法及系统 | |
CN107436450A (zh) | 一种基于连续小波变换的地震信号带宽拓展方法 | |
CN105974468A (zh) | 一种能够同时进行五维地震数据重建和噪声压制的方法 | |
CN107132579A (zh) | 一种保地层结构的地震波衰减补偿方法 | |
CN105182417B (zh) | 一种基于形态成分分析的面波分离方法及系统 | |
CN106707342B (zh) | 共炮点道集多级面波压制方法和装置 | |
Wang et al. | Seismic data denoising for complex structure using BM3D and local similarity | |
CN105652322A (zh) | 多分量地震数据的t-f-k域极化滤波方法 | |
CN106707334A (zh) | 一种提高地震资料分辨率的方法 | |
CN106526678A (zh) | 一种反射声波测井的波场分离方法及装置 | |
CN107678064A (zh) | 一种声波时差实时提取方法 | |
CN108646289B (zh) | 一种估计地震品质因子的方法 | |
CN103744114B (zh) | 基于零偏垂直地震剖面数据估计品质因子的方法和装置 | |
CN107703546A (zh) | 一种基于小波变换的新阈值函数地震资料去噪方法 | |
CN104391324A (zh) | 依赖频率的avo反演前的地震道集动校拉伸校正预处理技术 | |
CN111708087A (zh) | 一种基于DnCNN神经网络对地震数据噪声压制的方法 | |
CN104635264B (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 |