CN109242781B - 保地质边界的地震图像去噪方法及计算机可读存储介质 - Google Patents
保地质边界的地震图像去噪方法及计算机可读存储介质 Download PDFInfo
- Publication number
- CN109242781B CN109242781B CN201710557906.7A CN201710557906A CN109242781B CN 109242781 B CN109242781 B CN 109242781B CN 201710557906 A CN201710557906 A CN 201710557906A CN 109242781 B CN109242781 B CN 109242781B
- Authority
- CN
- China
- Prior art keywords
- seismic image
- representing
- denoising
- constructing
- geological
- 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 27
- 238000009792 diffusion process Methods 0.000 claims abstract description 30
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 2
- 238000013508 migration Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开一种保地质边界的地震图像去噪方法和计算机可读存储介质,该方法包括:针对原始地震图像中任一点,构建结构张量;针对原始地震图像中任一点,基于结构张量,构建扩散张量场;基于扩散张量场,构建保地质边界去噪的非线性偏微分方程;将非线性偏微分方程转化为线性方程组;对线性方程组进行求解,获得去噪后的地震图像。该方法具有较强的稳定性,能够有效保持并突出地质边界特征,具有很强的保真性。
Description
技术领域
本发明涉及油气勘探开发中的地震处理与解释领域,更具体地,涉及一种保地质边界的地震图像去噪方法。
背景技术
地震勘探技术的目标是利用地震处理与解释技术,对地下结构进行定位、识别和描述,为地下油气藏的勘探提供直观、可靠的依据。如何得到高质量的地下地震图像(地震深度偏移剖面)是地震处理与解释技术的核心内容之一。
由于地震数据在采集过程中存在噪声,直接利用偏移成像技术处理得到的地震图像往往存在信噪比低的问题,严重制约了后续的地震解释精度,因此有必要对地震图像实施去噪处理。常规偏移后去噪技术(如F-K滤波、随机噪声衰减)利用有效地震信号与噪声在频谱、斜率、幅值等方面的区别实施信号与噪声的分离,在去除噪声的同时容易损失有效信号,因此保真性不足,制约了其应用前景。
发明内容
本发明的目的是提出一种保地质边界的地震图像去噪方法,其能够克服现有方法保真性不足的缺陷。
本发明的一方面提出一种保地质边界的地震图像去噪方法,包括:
针对原始地震图像中任一点,构建结构张量;
针对所述原始地震图像中任一点,基于所述结构张量,构建扩散张量场;
基于所述扩散张量场,构建保地质边界去噪的非线性偏微分方程;
将所述非线性偏微分方程转化为线性方程组;
对所述线性方程组进行求解,获得去噪后的地震图像。
优选地,所述结构张量表示为:
其中,G表示结构张量,gx与gz分别表示地震图像沿水平方向和垂直方向的梯度,<·>表示二维高斯光滑滤波。
优选地,所述扩散张量场为:
其中,D表示扩散张量场,λ1表示结构张量G的最大特征值,λ2表示结构张量G的最小特征值,v1表示局部图像法向方向单位向量,其正交于地震图像的主结构方向,v2表示局部图像切向方向单位向量,其平行于地震图像的主结构方向。
优选地,所述保地质边界去噪的非线性偏微分方程为:
优选地,所述去噪因子α的取值范围在1至20之间。
优选地,所述线性方程组为:
优选地,利用共轭梯度算法对所述线性方程组进行求解。
本发明另一方面提供一种计算机可读存储介质,其上存储有计算机程序,其中,所述程序被处理器执行时实现以下步骤:
针对原始地震图像中任一点,构建结构张量;
针对所述原始地震图像中任一点,基于所述结构张量,构建扩散张量场;
基于所述扩散张量场,构建保地质边界去噪的非线性偏微分方程;
将所述非线性偏微分方程转化为线性方程组;
对所述线性方程组进行求解,获得去噪后的地震图像。
本发明的有益效果在于:(1)利用地震图像的结构张量,结构张量具有适应低信噪比地震数据的优势,因此本发明的保地质边界的地震图像去噪方法具有较强的稳定性;(2)能够有效保持并突出地质边界特征,具有很强的保真性。
本发明的方法和装置具有其它的特性和优点,这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的,或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述,这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施例中,相同的附图标记通常代表相同部件。
图1示出了根据本发明实施例的保地质边界的地震图像去噪方法的流程图;
图2示出了本发明实施例的原始含噪声地震图像;
图3示出了应用根据本发明实施例的保地质边界的地震图像去噪方法对图2进行处理之后的地震图像。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了使本发明更加透彻和完整,并且能够将本发明的范围完整地传达给本领域的技术人员。
图1示出了根据本发明实施例的保地质边界的地震图像去噪方法的流程图,如图1所示,根据本发明实施例的保地质边界的地震图像去噪方法包括以下步骤:
步骤1:针对原始地震图像中任一点,构建结构张量。
结构张量包含了地震图像的局部走向信息和法向信息。假设H为二维地震图像,二维地震图像H中表示空间方向信息的结构张量由图像梯度值定义,结构张量表示区域的变化方向和沿变化方向的变化量大小,地震地层纹理和断层纹理由局部各点方位信息变化关系确定。引入Gaussian函数模糊局部细节,使得结构张量突出显示区域内信号的复杂性。对二维图像,任一点的结构张量G是一个2*2的矩阵:
其中,gx与gz分别表示地震图像沿水平方向和垂直方向的梯度,<·>表示二维高斯光滑滤波。根据原始地震图像,能够确定gx与gz的值,即可以确定任一点的结构张量G。
步骤2:针对所述原始地震图像中任一点,基于所述结构张量,构建扩散张量场
结构张量G是半正定矩阵,针对半正定矩阵G,其特征值与特征向量可由求解方程|G-λI|=0得到,其中λ表示特征值,I表示单位阵,结构张量G可以表示为:
其中,
λ1表示结构张量G的最大特征值,对应于结构张量能量在第一个特征张量方向v1的能量,
λ2表示结构张量G的最小特征值,对应于结构张量能量在第二个特征张量方向v2的能量,
(λ1-λ2)/λ1表示局部线性指标,反映局部方向的一致性。
特征向量描述了地震图像局部线性结构的方向性,针对地震图像的每个点,第一个特征向量v1正交于地震图像的主结构方向,第二个特征向量v2平行于地震图像的主结构方向。
因此,根据结构张量算法的物理意义,可以计算得到地震图像中任一点的局部线性指标(λ1-λ2)/λ1、局部图像法向方向单位向量v1和局部图像切向方向单位向量v2。值得注意的是,结构张量算法能适应低信噪比地震数据,因此能用于稳健地拾取地下局部地层方向信息。
基于结构张量,可构建扩散张量场,扩散张量场D的表达式如下:
步骤3:基于扩散张量场,构建保地质边界去噪的非线性偏微分方程。
构建如下保地质边界去噪的非线性偏微分方程:
其中,X=(x,z),表示原始地震图像中的任意一点,x和z分别表示该点在原始地震图像中沿水平方向和垂直方向的坐标,f(X)表示原始地震图像,s(X)表示去噪后的地震图像,和分别表示梯度算子和散度算子,α表示去噪因子,其是正实数,用于控制去噪力度,其取值一般在1~20之间,值越大代表去噪力度越大,D(X)表示扩散张量场,其由公式(3)表示。
步骤4:将非线性偏微分方程转化为线性方程组。
利用有限差分近似,将非线性偏微分方程(3)转化为如下线性方程组:
步骤5:对线性方程组进行求解,获得去噪后的地震图像。
利用共轭梯度算法对线性方程组(5)进行求解,即可获得去噪后的地震图像s。
本发明还提供一种计算机可读存储介质,其上存储有计算机程序,其中,所述程序被处理器执行时实现以下步骤:
针对原始地震图像中任一点,构建结构张量;
针对所述原始地震图像中任一点,基于所述结构张量,构建扩散张量场;
基于所述扩散张量场,构建保地质边界去噪的非线性偏微分方程;
将所述非线性偏微分方程转化为线性方程组;
对所述线性方程组进行求解,获得去噪后的地震图像。
在一个示例中,所述结构张量表示为:
其中,G表示结构张量,gx与gz分别表示地震图像沿水平方向和垂直方向的梯度,<·>表示二维高斯光滑滤波。
在一个示例中,所述扩散张量场为:
其中,λ1表示结构张量G的最大特征值,λ2表示结构张量G的最小特征值,v1表示局部图像法向方向单位向量,其正交于地震图像的主结构方向,v2表示局部图像切向方向单位向量,其平行于地震图像的主结构方向。
在一个示例中,所述保地质边界去噪的非线性偏微分方程为:
在一个示例中,所述去噪因子α的取值范围在1至20之间。
在一个示例中,所述线性方程组为:
在一个示例中,利用共轭梯度算法对所述线性方程组进行求解。
实施例
应用根据本发明实施例的保地质边界的地震图像去噪方法对图2所示的原始含噪声地震图像进行去噪处理。从图2可以看出,在该原始含噪声地震图像中,图像地质特征被噪声淹没,难以识别地质构造特征,严重制约了地震解释的合理性与可靠性。
应用以上步骤1-5对原始地震图像进行去噪处理,所得到的地震图像如图3所示。与图2相比,该地震图像的地震层位和断层展布特征得到有效加强,且各地质边界得到有效保持。以此地震图像作为地震解释的输入,能有效提高地震解释的精度,降低勘探开发风险。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (7)
4.根据权利要求1所述的保地质边界的地震图像去噪方法,其中,所述去噪因子α的取值范围在1至20之间。
6.根据权利要求5所述的保地质边界的地震图像去噪方法,其中,利用共轭梯度算法对所述线性方程组进行求解。
7.一种计算机可读存储介质,其上存储有计算机程序,其中,所述程序被处理器执行时实现以下步骤:
针对原始地震图像中任一点,构建结构张量;
针对所述原始地震图像中任一点,基于所述结构张量,构建扩散张量场;
基于所述扩散张量场,构建保地质边界去噪的非线性偏微分方程;
将所述非线性偏微分方程转化为线性方程组;
对所述线性方程组进行求解,获得去噪后的地震图像;
所述保地质边界去噪的非线性偏微分方程为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710557906.7A CN109242781B (zh) | 2017-07-10 | 2017-07-10 | 保地质边界的地震图像去噪方法及计算机可读存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710557906.7A CN109242781B (zh) | 2017-07-10 | 2017-07-10 | 保地质边界的地震图像去噪方法及计算机可读存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109242781A CN109242781A (zh) | 2019-01-18 |
CN109242781B true CN109242781B (zh) | 2020-12-01 |
Family
ID=65083683
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710557906.7A Active CN109242781B (zh) | 2017-07-10 | 2017-07-10 | 保地质边界的地震图像去噪方法及计算机可读存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109242781B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103489159A (zh) * | 2013-09-02 | 2014-01-01 | 电子科技大学 | 基于三边结构导向滤波的三维地震数据图像降噪方法 |
CN104020492A (zh) * | 2013-07-01 | 2014-09-03 | 西安交通大学 | 一种三维地震资料的保边滤波方法 |
WO2015164388A1 (en) * | 2014-04-22 | 2015-10-29 | Westerngeco Llc | Seismic data processing |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8274859B2 (en) * | 2010-02-22 | 2012-09-25 | Landmark Graphics Corporation | Systems and methods for modeling 3D geological structures |
-
2017
- 2017-07-10 CN CN201710557906.7A patent/CN109242781B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104020492A (zh) * | 2013-07-01 | 2014-09-03 | 西安交通大学 | 一种三维地震资料的保边滤波方法 |
CN103489159A (zh) * | 2013-09-02 | 2014-01-01 | 电子科技大学 | 基于三边结构导向滤波的三维地震数据图像降噪方法 |
WO2015164388A1 (en) * | 2014-04-22 | 2015-10-29 | Westerngeco Llc | Seismic data processing |
Non-Patent Citations (3)
Title |
---|
Contrast Enhancement Based Denoising Method in Diffusion Tensor Imaging;Solwin Johnson 等;《2013 Annual IEEE India Conference》;20131213;第1-4页 * |
基于非线性扩散滤波结构信息的图像去噪方法;张建伟 等;《计算机工程与设计》;20161130;第37卷(第11期);第2972-2978页 * |
对地震图像进行保边滤波的;王绪松 等;《地球物理学进展》;20060630;第21卷(第2期);第452-457页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109242781A (zh) | 2019-01-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10234581B2 (en) | System and method for high resolution seismic imaging | |
CN109143363B (zh) | 海洋拖缆双检采集鬼波压制方法及系统 | |
CN108828670B (zh) | 一种地震数据降噪方法 | |
Florio et al. | On the estimation of the structural index from low-pass filtered magnetic data | |
Wang et al. | Robust vector median filtering with a structure-adaptive implementation | |
CN110687597B (zh) | 一种基于联合字典的波阻抗反演方法 | |
CN111781635B (zh) | 海底四分量弹性波高斯束深度偏移方法和装置 | |
Wang et al. | Denoising with weak signal preservation by group-sparsity transform learning | |
Liu et al. | A dictionary learning method with atom splitting for seismic footprint suppression | |
Wang et al. | An explicit stabilization scheme for Q-compensated reverse time migration | |
Wang et al. | Robust singular value decomposition filtering for low signal-to-noise ratio seismic data | |
Xu et al. | Ground-roll separation of seismic data based on morphological component analysis in two-dimensional domain | |
CN109242781B (zh) | 保地质边界的地震图像去噪方法及计算机可读存储介质 | |
Oboué et al. | An advanced median filter for improving the signal-to-noise ratio of seismological datasets | |
CN116719088A (zh) | 一种航空瞬变电磁数据噪声压制方法 | |
CN109655916B (zh) | 用于分离地震数据中有效波与多次波的方法及系统 | |
Tian et al. | Joint sparsity recovery for noise attenuation | |
Gómez et al. | Spectral structure-oriented filtering of seismic data with self-adaptive paths | |
CN111352158B (zh) | 地震信号增强方法及装置 | |
CN112083490B (zh) | 地震数据噪声衰减方法及装置 | |
Guo et al. | Structure-oriented filtering for seismic images using nonlocal median filter | |
CN109242770B (zh) | 图像导引的地震速度插值方法及计算机可读存储介质 | |
AU2011248910B2 (en) | System and method for compensating time and offset varying near-surface effects in seismic data background | |
Zand et al. | Least-squares reverse time migration with shifted total variation regularization | |
CN105301655A (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 |