CN108919355B - 基于结构张量导引的高维s变换方法 - Google Patents
基于结构张量导引的高维s变换方法 Download PDFInfo
- Publication number
- CN108919355B CN108919355B CN201810455272.9A CN201810455272A CN108919355B CN 108919355 B CN108919355 B CN 108919355B CN 201810455272 A CN201810455272 A CN 201810455272A CN 108919355 B CN108919355 B CN 108919355B
- Authority
- CN
- China
- Prior art keywords
- seismic data
- dimensional
- structure tensor
- transformation
- matrix
- 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
- 238000011426 transformation method Methods 0.000 title claims abstract description 13
- 230000009466 transformation Effects 0.000 claims abstract description 34
- 239000011159 matrix material Substances 0.000 claims abstract description 31
- 238000000034 method Methods 0.000 claims abstract description 9
- 238000000605 extraction Methods 0.000 claims abstract description 7
- 230000008569 process Effects 0.000 claims description 4
- 230000000694 effects Effects 0.000 abstract description 6
- 229910052500 inorganic mineral Inorganic materials 0.000 abstract description 2
- 239000011707 mineral Substances 0.000 abstract description 2
- 238000004458 analytical method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000005070 sampling Methods 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/32—Transforming one recording into another or one representation into another
-
- 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
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变换方法,包括以下内容:构建结构张量矩阵,对地震数据通过结构张量矩阵的特征值与特征向量进行各向异性参数提取;构建基于结构张量导引的高维S变换算子;将地震数据采用基于结构张量导引的高维S变换算子进行高维S变换,得到符合地震数据特征的稀疏表征数据。本发明可以根据结构张量分析出的地震数据的特征信息,对地震数据采用改进的高维S变换,从而使地震数据特征更加突显、细化以及稀疏,提高了地震数据在处理解释时的信息利用率,有助于地震数据属性的提取,具有重要的理论和应用价值,对油气、矿产资源勘探的发展有一定的推动作用。
Description
技术领域
本发明是关于一种基于结构张量导引的高维S变换方法,涉及地震勘探领域,特别涉及高维、冗余地震数据的特征稀疏表达领域。
背景技术
时频分析作为新兴的非平稳信号处理方法,近年来受到越来越多的关注,时频分析或称时频分布,是描述信号频率随时间变化的信号处理方法。采用时间-频率联合分布表示信号,将一维的时间域信号映射到一个二维的时频平面,在时频域内对信号进行分析,全面反映观测信号的时间-频率联合特征,使研究人员同时掌握信号的时域及频域信息,而且还可以了解信号的频率是如何随时间变化的。
传统的时频变换方法有短时傅里叶变换、Gabor变换、小波变换、S变换等。其中,S变换结合了短时傅立叶变换和小波变换的优点而避免了它们的不足:频率的倒数决定了S变换中的高斯窗的尺度大小,因此具有小波变换的多分辨率分析;而且S变换含有相位因子,保留了每个频率的绝对相位特征,这是小波变换所不具有的特性。另外,S变换对于多分量信号具有可加性,是一种线性变换,具有无损可逆性。但是,一维的时频变换缺乏空间分辨率。此外,传统时频变换只是在基函数种类固定的局限下,在处理非平稳信号时,其窗口调节的灵活性略显不足,物理意义也不明确。因此,有必要将一维S变换推广到高维情况。
发明内容
针对上述问题,本发明的目的是提供一种能够突显地震数据稀疏特征,增加地震信息的利用率,有助于地震数据属性提取的基于结构张量导引的高维S变换方法。
为实现上述目的,本发明采取以下技术方案:一种基于结构张量导引的高维S变换方法,包括以下内容:
步骤1):构建结构张量矩阵,对地震数据通过结构张量矩阵的特征值与特征向量进行各向异性参数提取;
步骤2):构建基于结构张量导引的高维S变换算子;
步骤3):将地震数据采用基于结构张量导引的高维S变换算子进行高维S变换,得到符合地震数据特征的稀疏表征数据。
进一步地,所述步骤1)的具体过程为;假设由地震振幅构成的标量场为u,定义结构张量矩阵:
其中,x,z为地震数据的空间坐标;
结构张量矩阵S的特征矩阵为:
(S-λE)x=0
其中,λ为矩阵S对应的特征值,E为单位矩阵,x则为矩阵S的特征向量;结构张量的特征值为:
根据特征值得到特征向量:
通过结构张量的特征值与特征向量,得到地震数据的各向异性参数:
其中,p,α,β为各向异性参数,ε1和ε2为平衡系数,取值在0~1之间。进一步地,高维S变换算子为:
(x-x′)p=(x-x′)·cos p+(z-z′)·sin p
(z-z′)p=-(x-x′)·sin p+(z-z′)·cos p
式中,f为频率,h(x,z)为地震数据,(x,z)与(x′,z′)分别为S变换前后的地震数据空间坐标。
本发明由于采取以上技术方案,其具有以下优点:1、本发明引入了结构张量的概念,增加基函数灵活性的同时,更加稀疏的突显高维、冗余地震数据的特征,从而将物理意义明确、稀疏性更强的地震数据用于地震资料处理解释中,增加地震信息的利用率,提高处理解释的效果。2、本发明将结构张量矩阵与高维S变换相结合,凸显地震数据特征,得到更加符合地震数据特征的稀疏表征数据。3、本发明可以根据结构张量分析出的地震数据的特征信息,对地震数据采用改进的高维S变换,从而使地震数据特征更加突显、细化以及稀疏,提高了地震数据在处理解释时的信息利用率,有助于地震数据属性(方向、尺度、走时等)的提取,具有重要的理论和应用价值,对油气、矿产资源勘探的发展有一定的推动作用。
附图说明
图1是本发明的基于结构张量导引的高维S变换方法流程图;
图2是本发明实施例包含3个地震同相轴的地震剖面效果图;
图3是本发明实施例基于结构张量导引的高维S变换稀疏特征表达的地震数据效果图,图(a)、(b)和(c)分别代表同一频率下,不同方向的特征表达的地震数据效果图。
具体实施方式
以下结合附图来对本发明进行详细的描绘。然而应当理解,附图的提供仅为了更好地理解本发明,它们不应该理解成对本发明的限制。
如图1所示,本发明的基于结构张量导引的高维S变换方法,包括以下内容:
1、构建结构张量矩阵,对地震数据通过结构张量矩阵的特征值与特征向量进行各向异性参数提取。
最初,由Di Zenzo和Forstner先后于1986年和1987年提出结构张量的定义,并将其应用于特征检测。假设由地震振幅构成的标量场为u,定义结构张量矩阵(梯度平方矩阵):
其中,x,z为地震数据的空间坐标。
结构张量矩阵S的特征矩阵为:
(S-λE)x=0 (2)
其中,λ为结构张量矩阵S对应的特征值,E为单位矩阵,x则为结构张量矩阵S的特征向量。
有非零解的充分必要条件是(2)对应的特征多项式为零,即
利用求根公式:
结构张量的特征值描述了图像在一定邻域内在特征向量方向上的数值变化平均量。当λ1>>λ2≈0时,表示该点附近的图像呈线性结构;当λ1>>0且λ2>>0时,表示该点附近至少存在两个线性结构;当λ1≈0且λ2≈0时,表示该点附近图像缓变,没有线性结构。
根据特征值(4)可以得到特征向量:
较大特征值对应的特征向量x1的方向为图像中梯度变化最大的方向(orthogonalto linear features),较小特征值对应的特征向量x2与x1垂直(parallel to linearfeatures)。
通过结构张量给出的特征值与特征向量,便能根据公式(6)得到地震数据的各向异性参数:
其中,p,α,β为各向异性参数,ε1和ε2为平衡系数,调节各向异性参数的幅值,一般在0~1之间。
2、构建物理意义更加明确的基于结构张量导引的高维S变换算子。
传统的高维S变换(2D为例)是Stockwell等人1997年在一维S变换的基础上提出的,二维S变换继承了一维S变换的优点,在空域与频域都有很高的分辨率。这使得它在二维图像处理领域具有广阔的应用前景。
传统二维S变换可写成:
其中,h(x,z)为地震数据,(x,z)与(x′,z′)分别为S变换前后的地震数据空间坐标,kx和kz为x方向和z方向上的波数,公式(7)可以看作是两次一维S变换。
在传统二维S变换基础上,将传统的时频分析扩展到时间-空间频率-各向异性空间,力图通过扩展空间,进一步凸显细化数据的某些特征。然后,引入了结构张量矩阵,在高维、冗余地震信号中稀疏提取特征明确的数据。
改进后的高维S变换算子为:
(x-x′)p=(x-x′)·cos p+(z-z′)·sin p (9)
(z-z′)p=-(x-x′)·sin p+(z-z′)·cos p (10)
其中,f表示频率,与公式(7)对比,在传统S变换的基础上,引进了多个物理意义明确的参数,可以反映地震数据的空间方向,尺度信息。
通过公式(8),便可以扩展时频空间到更高维更冗余的数据空间,从而将地震数据的特征细化突显。
3、采用基于结构张量导引的高维S变换对地震数据进行稀疏特征表达。
改进后的高维S变换公式(8)中的各向异性参数(p,α,β)可以通过公式(6)得到。将地震数据采用基于结构张量导引的高维S变换算子进行高维S变换,便能得到更加符合地震数据特征的稀疏表征数据。
如图2~3所示,下面通过具体实施例详细说明本发明的基于结构张量导引的高维S变换方法的具体过程。本实施例对包含3个地震同相轴的地震剖面进行测试,具体过程:
(1)输入地震剖面数据,根据公式(4)和(5)计算得到对应地震剖面的结构张量矩阵的特征值与特征向量,从而得到相应的地震数据的各向异性参数(p,α,β);
(2)根据给出频率值f,一般满足香农采样定理即可。然后,结合地震剖面结构张量矩阵提供的各向异性参数与公式(8),便能够对地震剖面进行高维S变换,提取特征更加细化的、稀疏的、不同频率、不同方向、不同尺度的地震数据信息,有助于进一步的地震资料处理解释工作。
综上所述,与传统的二维S变换方法相比,本发明的基于结构张量导引的高维S变换方法,不仅含有地震数据的时频-空间信息,还突显了地震数据的方向、尺度信息。这样,使得地震数据得到了细化,针对地震资料处理解释的需求,可以更加稀疏的特征表达地震数据,提高了地震数据信息的利用率。
上述各实施例仅用于说明本发明,其中方法的各实施步骤等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。
Claims (1)
1.一种基于结构张量导引的高维S变换方法,其特征在于包括以下内容:
步骤1):构建结构张量矩阵,对地震数据通过结构张量矩阵的特征值与特征向量进行各向异性参数提取,具体过程为:
假设由地震振幅构成的标量场为u,定义结构张量矩阵:
其中,x,z为地震数据的空间位置;
结构张量矩阵S的特征矩阵为:
(S-λE)w=0
其中,λ为矩阵S对应的特征值,E为单位矩阵,w则为矩阵S的特征向量;
结构张量的特征值为:
根据特征值得到特征向量:
通过结构张量的特征值与特征向量,得到地震数据的各向异性参数:
其中,p,α,β为各向异性参数,ε1和ε2为平衡系数,取值在0~1之间;步骤2):构建基于结构张量导引的高维S变换算子,其中,高维S变换算子为:
(x-x′)p=(x-x′)·cosp+(z-z′)·sinp
(z-z′)p=-(x-x′)·sinp+(z-z′)·cosp
式中,f为频率,h(x,z)为地震数据,(x,z)与(x′,z′)分别为S变换前后的地震数据空间位置;
步骤3):将地震数据采用基于结构张量导引的高维S变换算子进行高维S变换,得到符合地震数据特征的稀疏表征数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810455272.9A CN108919355B (zh) | 2018-05-14 | 2018-05-14 | 基于结构张量导引的高维s变换方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810455272.9A CN108919355B (zh) | 2018-05-14 | 2018-05-14 | 基于结构张量导引的高维s变换方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108919355A CN108919355A (zh) | 2018-11-30 |
CN108919355B true CN108919355B (zh) | 2020-04-07 |
Family
ID=64402388
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810455272.9A Active CN108919355B (zh) | 2018-05-14 | 2018-05-14 | 基于结构张量导引的高维s变换方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108919355B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110261910A (zh) * | 2019-06-27 | 2019-09-20 | 中国石油化工股份有限公司 | 基于自适应稀疏s变换的地震数据面波去除方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8117014B2 (en) * | 2009-03-27 | 2012-02-14 | Schlumberger Technology Corporation | Methods to estimate subsurface deviatoric stress characteristics from borehole sonic log anisotropy directions and image log failure directions |
CN102288994B (zh) * | 2011-07-01 | 2013-04-03 | 中国海洋石油总公司 | Radon谱约束下高维地震数据规则化方法 |
CN103926616B (zh) * | 2014-04-11 | 2016-07-13 | 中国海洋石油总公司 | 一种基于叠前crp道集的多尺度各向异性扩散滤波方法 |
EP3189345A4 (en) * | 2014-09-01 | 2018-05-09 | Bioprotonics LLC | Selective sampling magnetic resonance-based method for assessing structural spatial frequencies |
WO2017160273A1 (en) * | 2016-03-14 | 2017-09-21 | Schlumberger Technology Corporation | Dip estimation via a modified structure tensor |
CN107272064A (zh) * | 2017-07-18 | 2017-10-20 | 中国石油化工股份有限公司 | 碳酸盐岩缝洞体内部结构的刻画方法 |
CN107272065A (zh) * | 2017-08-04 | 2017-10-20 | 中国石油化工股份有限公司 | 碳酸盐岩断溶体轮廓检测方法 |
-
2018
- 2018-05-14 CN CN201810455272.9A patent/CN108919355B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108919355A (zh) | 2018-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kushnir et al. | Anisotropic diffusion on sub-manifolds with application to earth structure classification | |
CN107688201B (zh) | 基于rbm地震叠前信号聚类方法 | |
CN109669213B (zh) | 基于优化Morlet小波的分频扩散滤波断层强化方法 | |
Sijtsma | Acoustic beamforming for the ranking of aircraft noise | |
Gao et al. | A new approach for extracting the amplitude spectrum of the seismic wavelet from the seismic traces | |
CN105093315B (zh) | 一种去除煤层强反射信号的方法 | |
Akinshin et al. | Accuracy of spike-train Fourier reconstruction for colliding nodes | |
CN103927761A (zh) | 一种基于稀疏表示的故障微弱信号特征提取方法 | |
Qi et al. | Seismic time–frequency spectrum analysis based on local polynomial Fourier transform | |
CN108919355B (zh) | 基于结构张量导引的高维s变换方法 | |
CN110413939B (zh) | 一种基于原子范数的到达角估计方法 | |
Yao-Jun et al. | Unsupervised seismic facies analysis using sparse representation spectral clustering | |
CN111880222A (zh) | 一种基于Shearlet变换的地震映像小断层识别增强方法 | |
Xin et al. | Seismic high-resolution processing method based on spectral simulation and total variation regularization constraints | |
Pan et al. | A Bregman adaptive sparse-spike deconvolution method in the frequency domain | |
CN110927781A (zh) | 一种面向高精度地震相干属性的地震数据优化方法及系统 | |
Liu et al. | Enhanced coherence using principal component analysis | |
Pan et al. | Automatic Seismic Lithology Interpretation via Multi-attribute Integrated Deep Learning | |
CN111766624B (zh) | 一种地震数据拓频处理方法、装置、储存介质及电子设备 | |
Ma et al. | Random noise suppression algorithm for seismic signals based on principal component analysis | |
CN112363217A (zh) | 一种地震数据随机噪声压制方法及系统 | |
Jianhua et al. | Filtering of nuclear magnetic resonance logging signal based on the generalized S transform and singular value decomposition | |
CN110542921B (zh) | 一种基于主分量分析的振幅横向变化刻画方法 | |
Bo et al. | 3-D structure-oriented adaptive Gaussian pyramid for seismic multiscale fracture detection | |
CN114428282B (zh) | 一种基于去尺度s变换的地震信号时频变换方法 |
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 |