CN113311485B - 地震沉积特征加强滤波方法及装置 - Google Patents
地震沉积特征加强滤波方法及装置 Download PDFInfo
- Publication number
- CN113311485B CN113311485B CN202010122692.2A CN202010122692A CN113311485B CN 113311485 B CN113311485 B CN 113311485B CN 202010122692 A CN202010122692 A CN 202010122692A CN 113311485 B CN113311485 B CN 113311485B
- Authority
- CN
- China
- Prior art keywords
- seismic
- dimensional
- smoothing
- filtering
- slice information
- 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
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000001914 filtration Methods 0.000 title claims abstract description 37
- 239000013049 sediment Substances 0.000 title claims abstract description 27
- 230000009466 transformation Effects 0.000 claims abstract description 39
- 230000008021 deposition Effects 0.000 claims abstract description 19
- 230000002708 enhancing effect Effects 0.000 claims abstract description 8
- 238000005728 strengthening Methods 0.000 claims abstract description 8
- 238000009499 grossing Methods 0.000 claims description 45
- 230000006870 function Effects 0.000 claims description 29
- 238000004422 calculation algorithm Methods 0.000 claims description 17
- 238000004590 computer program Methods 0.000 claims description 15
- 238000012545 processing Methods 0.000 claims description 14
- 238000012512 characterization method Methods 0.000 claims description 8
- 238000003860 storage Methods 0.000 claims description 8
- 238000006073 displacement reaction Methods 0.000 claims description 5
- 238000000605 extraction Methods 0.000 claims description 5
- 238000011160 research Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 3
- 230000002787 reinforcement Effects 0.000 claims description 3
- 230000003014 reinforcing effect Effects 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 13
- 230000004044 response Effects 0.000 abstract description 7
- 238000000354 decomposition reaction Methods 0.000 abstract description 4
- 230000007480 spreading Effects 0.000 abstract description 4
- 238000003892 spreading Methods 0.000 abstract description 4
- 238000004458 analytical method Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 7
- 238000004891 communication Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000000047 product Substances 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 108010001267 Protein Subunits Proteins 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000013050 geological sediment Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000013519 translation 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. analysis, for interpretation, for correction
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
-
- 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. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/148—Wavelet transforms
Abstract
本发明提供的一种地震沉积特征加强滤波方法及装置，基于二维小波变换获取地震切片的水平和垂向分量信息，其在不同方向上所蕴含和突出加强的地震沉积特征不尽相同：水平方向细节分量H.s(x,y)对于地震切片中地质体的横向展布特征具有明显地加强突出作用；与此对应，垂直方向细节分量V.s(x,y)对于地震切片中地质体的垂向展布特征则有明显地加强突出作用，此外开展地震切片的二维小波变换提取其不同方向的细节子图，从而分析各方位分解信息中所突出的地质体的地震沉积特征，进一步在结合工区地质概况的基础上，实现基于地震资料的地质体沉积特征加强和储层精细表征的目的。
Description
技术领域
本发明涉及地震勘探技术领域，更具体的，涉及一种地震沉积特征加强滤波方法及装置。
背景技术
在油气田勘探开发领域，地震资料通常具有覆盖范围广、横向连续性强等特点，其一直是储层沉积特征分析的重要指导信息；特别是在少井或无井工区，基于地震资料的储层沉积表征分析是寻找地下油气资源分布的重要方法手段。然而由于地下地层对地震波具有一定的吸收和衰减作用，通常地震勘探资料面临低信噪比、低主频、窄频宽等一系列制约因素，基于地震反射资料开展储层沉积特征分析通常识别的精度不足，严重影响了地震储层精细表征的准确性。
发明内容
为了解决上述问题的至少一个，本发明提供一种基于二维小波变换的地震沉积特征加强滤波方法，包括：
提取研究工区叠后地震资料的二维切片信息；
对所述二维切片信息进行平滑处理，得到平滑后的二维切片信息；
基于预设的二维小波变换滤波器对平滑后的二维切片信息进行滤波处理；
结合滤波后的所述二维切片信息的水平分量和垂向分量进行地震切片沉积特征解释。
在某些实施例中，还包括：
对所述二维切片信息采用设定的小波函数构建二维小波变换滤波器。
在某些实施例中，所述二维切片信息包括振幅、相位、相干、瞬时频率、曲率。
在某些实施例中，所述二维平滑算法包括：线性平滑算法、最临近平滑、双三次平滑以及三次样条平滑。
本发明还提供一种基于二维小波变换的地震沉积特征加强滤波装置，包括：
提取模块，提取研究工区叠后地震资料的二维切片信息；
平滑模块，对所述二维切片信息进行平滑处理，得到平滑后的二维切片信息；
滤波模块，基于预设的二维小波变换滤波器对平滑后的二维切片信息进行滤波处理；
特征解释模块，结合滤波后的所述二维切片信息的水平分量和垂向分量进行地震切片沉积特征解释。
在某些实施例中，还包括：
滤波器构建模块，对所述二维切片信息采用设定的小波函数构建二维小波变换滤波器。
在某些实施例中，所述二维切片信息包括振幅、相位、相干、瞬时频率、曲率。
在某些实施例中，所述二维平滑算法包括：线性平滑算法、最临近平滑、双三次平滑以及三次样条平滑。
在某些实施例中，一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现权利要求如上所述的地震沉积特征加强滤波方法的步骤。
在某些实施例中，一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现如上所述的地震沉积特征加强滤波方法的步骤。
本发明的有益效果如下：
本发明提供的一种地震沉积特征加强滤波方法及装置，基于二维小波变换获取地震切片的水平和垂向分量信息，其在不同方向上所蕴含和突出加强的地震沉积特征不尽相同：水平方向细节分量H.s(x,y)对于地震切片中地质体的横向展布特征具有明显地加强突出作用；与此对应，垂直方向细节分量V.s(x,y)对于地震切片中地质体的垂向展布特征则有明显地加强突出作用，此外开展地震切片的二维小波变换提取其不同方向的细节子图，从而分析各方位分解信息中所突出的地质体的地震沉积特征，进一步在结合工区地质概况的基础上，实现基于地震资料的地质体沉积特征加强和储层精细表征的目的。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种基于二维小波变换的地震沉积特征加强滤波方法流程示意图。
图2是本发明实施例输入的网格数值模型。
图3a至图3d是本实施例以图2二维网格模型为输入获取的二维小波变换结果，其中图3a是模型的低频分量信息；图3b是模型的水平分量信息；图3c是模型的垂直分量信息；图3d是模型的对角分量信息。
图4是本发明实施例输入的河道数值模型。
图5a是以图4中河道模型为输入，基于二维小波变换提取的模型水平分量信息。
图5b是以图4中河道模型为输入，基于二维小波变换提取的模型垂向分量信息。
图6是国外某工区的实际叠后地震切片。
图7a是以图6中叠后地震切片为输入，基于二维小波变换提取的该地震切片的水平分量信息。
图7b是以图6中叠后地震切片为输入，基于二维小波变换提取的该地震切片的垂向分量信息。
图8示出本发明实施例中基于二维小波变换的地震沉积特征加强滤波装置的结构示意图。
图9示出适于用来实现本申请实施例的电子设备的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
针对传统地震沉积分析方法所存在的问题，本发明的目的是提供一种新的基于二维小波变换的地震沉积特征加强滤波方法与装置。本发明是在研究了以下问题的基础之上提出的：(1)由于地层的吸收衰减作用，地震资料分辨率通常较低，地震切片上临层干涉现象严重，一定程度上影响了地下储层的地震沉积分析精度；(2)地震切片信息通常是地下具有不同展布方向的多个地质体的地震响应的综合表征，因此有必要对地震切片开展二维小波变换以分方位突出其地质体在特定方向上的地震沉积特征；(3)二维小波变换已广泛应用于图像处理等领域，其能够对图像开展有效滤波并提取低频近似分量以及对应的水平、垂直和对角三个方向上的高频分量。二维小波变换后，地震切片的水平和垂直高频分量对该方向上的沉积特征有明显加强突出效果，因此可基于二维小波变换对地震切片进行处理以达到突出加强其沉积特征的目的。
本发明具体采取以下工作步骤来实现上述技术方案：1)提取工区叠后地震资料的二维切片信息→2)对上述叠后地震切片采用平滑算法开展二维平滑处理→3)基于上述平滑后的叠后地震切片数据确定合适的小波函数类型→4)基于上述确定的小波函数构建二维小波变换滤波器→5)对平滑后的叠后地震切片开展二维小波变换→6)提取基于上述二维小波变换获取的叠后地震切片的水平和垂向分量信息→8)联合上述获取的叠后地震切片水平和垂向分量信息开展地震切片沉积特征精细解释。
地震切片技术，其通过联合地震沉积学分析手段判定储层地貌和岩性信息，从而达到分析储层成藏和判断沉积环境的目的。由于地震资料分辨率、信噪比的限制以及临层干涉等因素，基于地震切片开展储层沉积特征表征等仍存在很多问题：首先，由于不同垂向深度上临近的地质体之间的地震反射会产生干涉效应，从而导致地震垂向分辨率降低；其次，水平临近的多个地质体的地震反射波亦存在相互干涉的问题，从而降低地震资料的横向分辨能力；此外，地震波在地下传播过程中由于数据噪声、多次波等引起的杂乱反射特征以及地震偏移归位不精确等对基于地震切片识别地下地质体也具有一定的破坏作用。因此，有必要对地震切片开展进一步处理，以提高基于地震切片开展储层沉积分析的精度。
作为一种数学工具，小波变换在现今得到广泛应用，其与传统傅里叶变换或短时傅里叶变换相比对信号的处理能力更加强大。傅里叶变换缺乏刻画时域信号局部特征的能力，此外其对非平稳信号的处理效果也不太理想。作为一种时间-频率域的局部变换，小波变换将原始信号滤波为一系列具有不同时间-频率特征的小波，其不仅能够检测信号的高、低频成分，而且能够将高低频成分发生的时间与原始信号相对应，有效克服了傅里叶变换的不足。因此小波变换也被称为数学显微镜，可用于分析信号的局部特征。对于地震切片，可将其可视为一幅二维图像，从而可以借助小波变换开展图像处理以突出加强地震切片上的储层沉积特征。
地震切片自身特征决定了将小波变换应用到地震切片图像处理时必须升级为二维小波变换，其亦是在传统小波变换的基础上发展而来。通过对地震切片的二维图像开展二维小波变换，可以进一步获取其在水平、垂直、对角等方向上的高频细节信息以及对应的低频近似信息。由于地震切片可认为是地下具有不同方向特征的多个地质体的地震响应的综合表征，因此借助二维小波变换提取地震切片在水平方向和垂直方向的细节子图，由于其水平和垂向分量信息对于地下地质体在特定方向上的地震沉积特征有明显地加强突出作用，从而可以进一步联合上述水平和垂直方向子图加强突出地下地质体的地震响应特征，以提高地震储层沉积精细解释的精度。
图1示出了本发明实施例中一种基于二维小波变换的地震沉积特征加强滤波方法，包括：
S1：提取研究工区叠后地震资料的二维切片信息。
具体而言，基于叠后三维地震数据体提取特定地震层位约束和特定地震切片提取算法计算的地震切片信息，二维切片的网格模型如图2所示，地震切片类型包括但不限于：时间切片(切片垂向时间位置一致，也称为等T0切片)，水平切片(即平行于某一地层层位提取地震切片信息)，地层切片(在多个层位约束下的地震属性的等比提取切片)等。提取的地震切片可记为s(x,y,t)，其中s表示地震切片的属性信息，x,y为切片对应的平面位置坐标，t表示地震波双程旅行时，当s(x,y,t)中t为某一常数t0时，s(x,y,t)即为地震时间切片，可记为s(x,y,t＝t0)或s(x,y,t0)；地震切片的属性信息s包括但不限于：振幅、相位、相干、瞬时频率、曲率等属性。
S2：对所述二维切片信息进行平滑处理，得到平滑后的二维切片信息。
针对地震切片s(x,y,t)开展二维平滑处理，首先忽略其时间维度t从而将地震切片视为平面图像信息，记为s(x,y)。在此基础上，通过二维平滑算法过滤地震切片数据的部分高频噪声，已便提高二维小波变换后获取的地震切片的水平和垂直分量的精度和稳定性。二维平滑后的地震切片记为S(x,y)。
本发明可用的二维平滑算法类型包括但不限于：线性平滑(linear)、最临近平滑(nearest)、双三次平滑(cubic)以及三次样条平滑(spline)等。根据地震切片数据特征选择合适的二维平滑类型，本发明实施例中选择的平滑算法为基于三次样条(spline)的二维平滑算法。
S3：基于预设的二维小波变换滤波器对所述二维切片信息进行二维小波变换。
对于二维小波变换，小波基的选择一般应考虑小波的线性相位、紧支性、衰减性和正交性等：对于小波相位，如果其满足或近似满足线性特征，则可保证重构数据的精确性；小波函数的紧支宽度越窄则其衰减速度越快，其局部化特征越好，计算效率也越高；此外小波正交性越强，计算的低通和高通子带数据相关性越小。总之，所谓小波即是一种特殊的长度有限、平均值为零的波形，其具有支流分量为零的特点。
二维小波变换允许使用多种类型的正交和双正交小波，常用的小波族包括Daubechies、Coiflets、Biorthogonal以及Reverse Biorthogonal等。其中对于地震切片数据的二维小波变换，本发明实施例主要应用Daubechies小波族中的haar小波。
对于能量有限空间L2(R)中的一维信号f(t)，记f(t)∈L2(R)，其信号能量E满足：
其中，t表示信号时间，满足上式所有条件的函数集合就形成了能量有限空间L2(R)。记ψ(t)为小波基函数，对于上述一维信号f(t)开展连续小波变换可表述为：
其中，FWT为信号f(t)小波变换结果，表示发生了位移和伸缩的小波函数，σ为尺度因子表示小波函数的伸缩量，满足σ＞0且σ越大ψ(t/σ)越宽，τ表征小波的位移量，＜f(t),ψσ,τ(t)＞表示一维信号f(t)和发生位移和伸缩后的小波函数ψσ,τ(t)的内积。由上式可知，信号f(t)的小波变换相当于将信号f(t)投影到时间-尺度(τ-σ)二维空间内以便于提取其相关特征。
将上述一维小波变换推广到二维情形可得：
其中，FSWT(.)表示二维信号S(x,y)的连续二维小波变换结果，y表征第二维度；w(x,y)表示给定的二维基本小波，μx,μy分别表示小波函数在x,y方向上的位移量，保证了变换前后小波能量的一致性。
对于一维小波变换，设有L2(R)空间的子空间序列：其中设Vj的正交基函数是由一个尺度函数φ(x)经过伸缩和位移得到的：
设小波函数ψ(x)＝φ(2x)-φ(2x-1)且Wj是Vj相对于Vj+1的正交补空间，则Wj的正交基函数是ψ(x)由经过伸缩和平移得到：
构成Vj+1的正交基，φ(x)且ψ(x)和满足下列关系式：
其中，l(n),h(n)分别表示对应的低通和高通滤波器且h(n)＝(-1)nl(1-n)。
二维小波变换是在一维小波变换的基础上发展而来，其尺度函数和小波函数可由一维尺度函数与小波函数的张量积分得到：
在此基础上，对二维信号s(x,y)开展二维小波变换可表征为：
其中，A.s(x,y)为低频平滑逼近，H.s(x,y)为对应的水平方向细节子图，V.s(x,y)分别为垂直分量细节子图，D.s(x,y)为对角方向细节子图。
图3a至图3d是本实施例以图2二维网格模型为输入获取的二维小波变换结果：图3a是模型的低频分量信息；图3b是模型的水平分量信息；图3c是模型的垂直分量信息；图3d是模型的对角分量信息。如图所示，通过二维小波变换，模型中网格的水平和垂直边界分别在模型二维小波变换结果中的水平和垂直分量中得到加强突出。因此，可以通过对地震切片信息开展二维小波变换以突出加强地质体的地震沉积边界特征，从而提高地震沉积分析精度，降低地震资料解释的不确定性。
S4：结合滤波后的所述二维切片信息的水平分量和垂向分量进行地震切片沉积特征解释。
对于地震切片开展二维小波变换，获取其水平方向细节子图H.s(x,y)和垂直方向细节子图V.s(x,y)后，联合两个分量信息开展地震切片储层沉积特征的精细解释。基于二维小波变换获取地震切片的水平和垂向分量信息，其在不同方向上所蕴含和突出加强的地震沉积特征不尽相同：水平方向细节分量H.s(x,y)对于地震切片中地质体的横向展布特征具有明显地加强突出作用；与此对应，垂直方向细节分量V.s(x,y)对于地震切片中地质体的垂向展布特征则有明显地加强突出作用。为此，开展地震切片的二维小波变换提取其不同方向的细节子图，从而分析各方位分解信息中所突出的地质体的地震沉积特征，进一步在结合工区地质概况的基础上，实现基于地震资料的地质体沉积特征加强和储层精细表征的目的。
作为示例，图4是本发明实施例输入的河道数值模型。
图5a和图5b是以图4中河道模型为输入，基于二维小波变换提取的模型水平分量信息(图5a)和模型垂向分量信息(图5b)。如图所示，通过对上述河道模型开展二维小波变换，其提取的模型水平和垂直分量信息均对模型中的河道特征有明显地加强突出效果。与模型数据相比，河边边界在水平和垂直分量中均更加清晰，极大地提高了地质体的地震识别解释精度。图6是国外某工区的实际叠后地震切片，图中白色椭圆指示地震切片中的河道信息。图7a和图7b是以图6中叠后地震切片为输入，基于二维小波变换提取的该地震切片的水平分量信息(图7a)和垂向分量信息(图7b)。如图所示，通过对上述地震切片开展二维小波变换，其提取的水平和垂直分量信息均对地震切片中的河道特征有明显加强突出作用。根据图中椭圆指示可知，与原始地震切片相比，其水平分量对图中横向展布的河道特征有明显地加强突出作用；与此同时，垂向分布的河道展布特征在垂直分量上得到明显地加强突出。
从上述描述可以知晓，本发明提供的一种地震沉积特征加强滤波方法，通过二维小波变换，实现了地震切片上地层沉积信息的分层级提取，有利于加强突出目标地质体的地震沉积信息；同时提取的地震切片低频和对应残差(高频成分)均包含丰富的地震沉积信息，高低频信息的联合应用有利于提高地震沉积分析精度；进一步的，可基于地震切片信息清晰突出地质体边界，为储层横向展布精细刻画提供重要指导信息；最后本发明算法简单易行，运行效率高，无需人机交互，易于形成软件功能模块推广应用。
基于相同的发明构思，图8示出了本发明实施例提供一种基于二维小波变换的地震沉积特征加强滤波装置，包括：
提取模块1，提取研究工区叠后地震资料的二维切片信息；
平滑模块2，对所述二维切片信息进行平滑处理，得到平滑后的二维切片信息；
滤波模块3，基于预设的二维小波变换滤波器对平滑后的二维切片信息进行滤波处理；
特征解释模块4，结合滤波后的所述二维切片信息的水平分量和垂向分量进行地震切片沉积特征解释。
本发明提供的一种地震沉积特征加强滤波装置，基于二维小波变换获取地震切片的水平和垂向分量信息，其在不同方向上所蕴含和突出加强的地震沉积特征不尽相同：水平方向细节分量H.s(x,y)对于地震切片中地质体的横向展布特征具有明显地加强突出作用；与此对应，垂直方向细节分量V.s(x,y)对于地震切片中地质体的垂向展布特征则有明显地加强突出作用，此外开展地震切片的二维小波变换提取其不同方向的细节子图，从而分析各方位分解信息中所突出的地质体的地震沉积特征，进一步在结合工区地质概况的基础上，实现基于地震资料的地质体沉积特征加强和储层精细表征的目的。
在某些实施例中，还包括：
滤波器构建模块，对所述二维切片信息采用设定的小波函数构建二维小波变换滤波器。
在某些实施例中，所述二维切片信息包括振幅、相位、相干、瞬时频率、曲率。
在某些实施例中，所述二维平滑算法包括：线性平滑算法、最临近平滑、双三次平滑以及三次样条平滑。
上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为电子设备，具体的，电子设备例如可以为个人计算机、膝上型计算机、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
在一个典型的实例中电子设备具体包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如上所述的由客户端执行的方法，或者，所述处理器执行所述程序时实现如上所述的由服务器执行的方法。
下面参考图9，其示出了适于用来实现本申请实施例的电子设备的结构示意图。
如图9所示，电子设备包括中央处理单元(CPU)601，其可以根据存储在只读存储器(ROM)602中的程序或者从存储部分608加载到随机访问存储器(RAM))603中的程序而执行各种适当的工作和处理。在RAM603中，还存储有系统操作所需的各种程序和数据。CPU601、ROM602、以及RAM603通过总线604彼此相连。输入/输出(I/O)接口605也连接至总线604。
以下部件连接至I/O接口605：包括键盘、鼠标等的输入部分606；包括诸如阴极射线管(CRT)、液晶显示器(LCD)等以及扬声器等的输出部分607；包括硬盘等的存储部分608；以及包括诸如LAN卡，调制解调器等的网络接口卡的通信部分609。通信部分609经由诸如因特网的网络执行通信处理。驱动器610也根据需要连接至I/O接口606。可拆卸介质611，诸如磁盘、光盘、磁光盘、半导体存储器等等，根据需要安装在驱动器610上，以便于从其上读出的计算机程序根据需要被安装如存储部分608。
特别地，根据本发明的实施例，上文参考流程图描述的过程可以被实现为计算机软件程序。例如，本发明的实施例包括一种计算机程序产品，其包括有形地包含在机器可读介质上的计算机程序，所述计算机程序包括用于执行流程图所示的方法的程序代码。在这样的实施例中，该计算机程序可以通过通信部分609从网络上被下载和安装，和/或从可拆卸介质611被安装。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于硬件+程序类实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。虽然本说明书实施例提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或终端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本说明书实施例时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。本领域技术人员应明白，本说明书的实施例可提供为方法、系统或计算机程序产品。因此，本说明书实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本说明书实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书实施例的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。以上所述仅为本说明书实施例的实施例而已，并不用于限制本说明书实施例。对于本领域技术人员来说，本说明书实施例可以有各种更改和变化。凡在本说明书实施例的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本说明书实施例的权利要求范围之内。
Claims (10)
1.一种基于二维小波变换的地震沉积特征加强滤波方法，其特征在于，包括：
提取研究工区叠后地震资料的二维切片信息；
对所述二维切片信息进行平滑处理，得到平滑后的二维切片信息；
基于预设的二维小波变换滤波器对平滑后的二维切片信息进行滤波处理，二维小波变换通过以下公式实现：
其中，FSWT(.)表示二维信号S(x,y)的连续二维小波变换结果，y表征第二维度；w(x,y)表示给定的二维基本小波，μx,μy分别表示小波函数在x,y方向上的位移量，保证了变换前后小波能量的一致性；
结合滤波后的所述二维切片信息的水平分量和垂向分量进行地震切片沉积特征解释。
2.根据权利要求1所述的地震沉积特征加强滤波方法，其特征在于，还包括：
对所述二维切片信息采用设定的小波函数构建二维小波变换滤波器。
3.根据权利要求2所述的地震沉积特征加强滤波方法，其特征在于，所述平滑处理的算法包括：线性平滑算法、最临近平滑、双三次平滑以及三次样条平滑。
4.根据权利要求1所述的地震沉积特征加强滤波方法，其特征在于，所述二维切片信息包括振幅、相位、相干、瞬时频率、曲率。
5.一种基于二维小波变换的地震沉积特征加强滤波装置，其特征在于，包括：
提取模块，提取研究工区叠后地震资料的二维切片信息；
平滑模块，对所述二维切片信息进行平滑处理，得到平滑后的二维切片信息；
滤波模块，基于预设的二维小波变换滤波器对平滑后的二维切片信息进行滤波处理，二维小波变换通过以下公式实现：
其中，FSWT(.)表示二维信号S(x,y)的连续二维小波变换结果，y表征第二维度；w(x,y)表示给定的二维基本小波，μx,μy分别表示小波函数在x,y方向上的位移量，保证了变换前后小波能量的一致性；
特征解释模块，结合滤波后的所述二维切片信息的水平分量和垂向分量进行地震切片沉积特征解释。
6.根据权利要求5所述的地震沉积特征加强滤波装置，其特征在于，还包括：
滤波器构建模块，对所述二维切片信息采用设定的小波函数构建二维小波变换滤波器。
7.根据权利要求6所述的地震沉积特征加强滤波装置，其特征在于，所述二维切片信息包括振幅、相位、相干、瞬时频率、曲率。
8.根据权利要求6所述的地震沉积特征加强滤波装置，其特征在于，所述平滑处理的算法包括：线性平滑算法、最临近平滑、双三次平滑以及三次样条平滑。
9.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1至4任一项所述的地震沉积特征加强滤波方法。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现权利要求1至4任一项所述的地震沉积特征加强滤波方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010122692.2A CN113311485B (zh) | 2020-02-27 | 2020-02-27 | 地震沉积特征加强滤波方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010122692.2A CN113311485B (zh) | 2020-02-27 | 2020-02-27 | 地震沉积特征加强滤波方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113311485A CN113311485A (zh) | 2021-08-27 |
| CN113311485B true CN113311485B (zh) | 2024-01-30 |
Family
ID=77370103
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010122692.2A Active CN113311485B (zh) | 2020-02-27 | 2020-02-27 | 地震沉积特征加强滤波方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113311485B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116256801B (zh) | 2023-05-16 | 2023-07-21 | 中国科学院地质与地球物理研究所 | 基于图像融合的深地油气精准导航断层表征方法与系统 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6745129B1 (en) * | 2002-10-29 | 2004-06-01 | The University Of Tulsa | Wavelet-based analysis of singularities in seismic data |
| CN103454686A (zh) * | 2013-08-12 | 2013-12-18 | 中国石油天然气股份有限公司 | 基于地层切片的小尺度沉积相进行储层预测的方法及系统 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2035864B1 (en) * | 2006-06-21 | 2015-07-29 | CGG Jason (Netherlands) B.V. | Interpretation of geologic depositional systems |
-
2020
- 2020-02-27 CN CN202010122692.2A patent/CN113311485B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6745129B1 (en) * | 2002-10-29 | 2004-06-01 | The University Of Tulsa | Wavelet-based analysis of singularities in seismic data |
| CN103454686A (zh) * | 2013-08-12 | 2013-12-18 | 中国石油天然气股份有限公司 | 基于地层切片的小尺度沉积相进行储层预测的方法及系统 |
Non-Patent Citations (7)
| Title |
|---|
| 三维时频分析方法在地震层序分析中的应用;吴国忱,康仁华,印兴耀;石油大学学报(自然科学版)(第01期);第81-85页 * |
| 东营凹陷民丰北带沙四上亚段砂砾岩体层序划分;邵绪鹏;张立强;靳久强;张立宽;;科学技术与工程(第04期);第175-180页 * |
| 小波分析方法在提高地震分辨率及地震层序分析中的应用;吴国忱,印兴耀,张广智,张繁昌;石油物探(S1);第45-51页 * |
| 小波变换在某气田地震剖面滤波中的应用;刘佳宾;王四巍;王艳艳;丁聪;陈钧龙;;中国煤炭地质(第10期);第73-79页 * |
| 小波多尺度分析技术在地震信号处理中的应用;杨立强;;科技创新与应用(第08期);第287页 * |
| 小波多尺度边缘检测及其在裂缝预测中的应用;苟量,彭真明;石油地球物理勘探(第03期);第309-313页 * |
| 测井信号小波变换在层序地层划分中的应用――以南海北部莺歌海盆地D14井黄流组为例;王亚辉;赵鹏肖;左倩媚;钟泽红;李俊良;;沉积与特提斯地质(第04期);第58-62页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113311485A (zh) | 2021-08-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Li et al. | A method for low-frequency noise suppression based on mathematical morphology in microseismic monitoring | |
| Neelamani et al. | Coherent and random noise attenuation using the curvelet transform | |
| Panda et al. | Application of wavelet transforms to reservoir-data analysis and scaling | |
| Chandrasekhar et al. | Wavelet analysis of geophysical well-log data of Bombay offshore basin, India | |
| US20150120199A1 (en) | Multiple Domain Processing For Combining Reservoir Models and Seismic Data | |
| Chatterjee et al. | Multi-scale stochastic simulation with a wavelet-based approach | |
| CN102879822A (zh) | 一种基于Contourlet变换的地震多属性融合方法 | |
| Cao et al. | 3D seismic denoising based on a low-redundancy curvelet transform | |
| Karbalaali et al. | Seismic channel edge detection using 3D shearlets—a study on synthetic and real channelised 3D seismic data | |
| CN109669213A (zh) | 基于优化Morlet小波的分频扩散滤波断层强化方法 | |
| Liu et al. | Irregularly sampled seismic data reconstruction using multiscale multidirectional adaptive prediction-error filter | |
| CN113311485B (zh) | 地震沉积特征加强滤波方法及装置 | |
| Yuan et al. | Inverse spectral decomposition using an lp-norm constraint for the detection of close geological anomalies | |
| Duan et al. | High-efficiency and precision VSP wavefield separation method via DCT | |
| Chu et al. | Application of wavelet analysis to upscaling of rock properties | |
| Boustani et al. | Mapping channel edges in seismic data using curvelet transform and morphological filter | |
| CN113311486B (zh) | 地震沉积特征加强滤波方法及装置 | |
| Huang et al. | Frequency–Space-Dependent Smoothing Regularized Nonstationary Predictive Filtering | |
| Liu et al. | Improving the resolution of seismic data based on S-transform and modified variational mode decomposition, an application to Songliao Basin, Northeast China | |
| Xu et al. | Random noise attenuation using a structure-oriented adaptive singular value decomposition | |
| Likkason | Spectral analysis of geophysical data | |
| Gong et al. | Automatic time picking of microseismic data based on shearlet-AIC algorithm | |
| Zhang et al. | Separation of potential field data using 3-D principal component analysis and textural analysis | |
| Pleshkevich et al. | Parallel algorithm for one-way wave equation based migration for seismic imaging | |
| Li et al. | Adaptive linear TFPF for seismic random noise attenuation |
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 |