CN108919355A - 基于结构张量导引的高维s变换方法 - Google Patents
基于结构张量导引的高维s变换方法 Download PDFInfo
- Publication number
- CN108919355A CN108919355A CN201810455272.9A CN201810455272A CN108919355A CN 108919355 A CN108919355 A CN 108919355A CN 201810455272 A CN201810455272 A CN 201810455272A CN 108919355 A CN108919355 A CN 108919355A
- Authority
- CN
- China
- Prior art keywords
- structure tensor
- seismic data
- dimension
- 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.)
- Granted
Links
- 238000011426 transformation method Methods 0.000 title claims abstract description 12
- 239000011159 matrix material Substances 0.000 claims abstract description 26
- 238000000034 method Methods 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 4
- 239000000284 extract Substances 0.000 claims description 3
- 238000000605 extraction Methods 0.000 abstract description 6
- 229910052500 inorganic mineral Inorganic materials 0.000 abstract description 2
- 239000011707 mineral Substances 0.000 abstract description 2
- 230000009466 transformation Effects 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling 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/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 (3)
1.一种基于结构张量导引的高维S变换方法,其特征在于包括以下内容:
步骤1):构建结构张量矩阵,对地震数据通过结构张量矩阵的特征值与特征向量进行各向异性参数提取;
步骤2):构建基于结构张量导引的高维S变换算子;
步骤3):将地震数据采用基于结构张量导引的高维S变换算子进行高维S变换,得到符合地震数据特征的稀疏表征数据。
2.根据权利要求1所述的基于结构张量导引的高维S变换方法,其特征在于,所述步骤1)的具体过程为;
假设由地震振幅构成的标量场为u,定义结构张量矩阵:
其中,x,z为地震数据的空间位置;
结构张量矩阵S的特征矩阵为:
(S-λE)x=0
其中,λ为矩阵S对应的特征值,E为单位矩阵,x则为矩阵S的特征向量;
结构张量的特征值为:
根据特征值得到特征向量:
通过结构张量的特征值与特征向量,得到地震数据的各向异性参数:
其中,p,α,β为各向异性参数,ε1和ε2为平衡系数,取值在0~1之间。
3.根据权利要求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变换前后的地震数据空间位置。
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 true CN108919355A (zh) | 2018-11-30 |
CN108919355B 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) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110261910A (zh) * | 2019-06-27 | 2019-09-20 | 中国石油化工股份有限公司 | 基于自适应稀疏s变换的地震数据面波去除方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010111077A1 (en) * | 2009-03-27 | 2010-09-30 | Schlumberger Canada Limited | Methods to estimate subsurface deviatoric stress characteristics from borehole sonic log anisotropy directions and image log failure directions |
CN102288994A (zh) * | 2011-07-01 | 2011-12-21 | 中国海洋石油总公司 | Radon谱约束下高维地震数据规则化方法 |
CN103926616A (zh) * | 2014-04-11 | 2014-07-16 | 中国海洋石油总公司 | 一种基于叠前crp道集的多尺度各向异性扩散滤波方法 |
US20160061917A1 (en) * | 2014-09-01 | 2016-03-03 | 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
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010111077A1 (en) * | 2009-03-27 | 2010-09-30 | Schlumberger Canada Limited | Methods to estimate subsurface deviatoric stress characteristics from borehole sonic log anisotropy directions and image log failure directions |
CN102288994A (zh) * | 2011-07-01 | 2011-12-21 | 中国海洋石油总公司 | Radon谱约束下高维地震数据规则化方法 |
CN103926616A (zh) * | 2014-04-11 | 2014-07-16 | 中国海洋石油总公司 | 一种基于叠前crp道集的多尺度各向异性扩散滤波方法 |
US20160061917A1 (en) * | 2014-09-01 | 2016-03-03 | 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 | 中国石油化工股份有限公司 | 碳酸盐岩断溶体轮廓检测方法 |
Non-Patent Citations (9)
Title |
---|
BO FENG ET AL.: "Data domain automatic reflection traveltime inversion via sparse decomposition of seismic data", 《2017 SEG INTERNATIONAL EXPOSITION AND 87TH ANNUAL MEETING》 * |
TRYGVE RANDEN ET AL.: "Three-Dimensional Texture Attributes for Seismic Data Analysis", 《SEG 2000 EXPANDED ABSTRACTS》 * |
YUBANG ZHOU ET AL.: "A study of GST method in the seismic attributes", 《2018 CPS/SEG INTERNATIONAL GEOPHYSICAL CONFERENCE》 * |
李一等: "梯度结构张量分析法在三维地震资料河道砂体预测中的应用", 《天然气工业》 * |
王德营等: "二维S变换域自适应面波压制", 《中国地球物理2013》 * |
王晓凯等: "基于高维连续小波变换的地震资料不连续性检测方法研究", 《地球物理学报》 * |
王清振等: "基于高维小波变换的高抗噪性边缘检测技术", 《石油地球物理勘探》 * |
陈强等: "《沁水盆地高瓦斯矿井三维地震技术应用实践》", 30 September 2012 * |
陈海燕等: "广义S变换及其在大地电磁测深数据处理中的应用", 《地球物理学进展》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110261910A (zh) * | 2019-06-27 | 2019-09-20 | 中国石油化工股份有限公司 | 基于自适应稀疏s变换的地震数据面波去除方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108919355B (zh) | 2020-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104122540B (zh) | 基于匹配字典和压缩感知的雷达一维距离像目标识别方法 | |
CN110361778B (zh) | 一种基于生成对抗网络的地震数据重建方法 | |
CN107688201B (zh) | 基于rbm地震叠前信号聚类方法 | |
CN102879822A (zh) | 一种基于Contourlet变换的地震多属性融合方法 | |
CN105607122B (zh) | 一种基于全变分地震数据分解模型的地震纹理提取与增强方法 | |
NO345726B1 (no) | Fremgangsmåte for seismisk tolkning ved bruk av seismiske teksturattributter | |
Gao et al. | A new approach for extracting the amplitude spectrum of the seismic wavelet from the seismic traces | |
Li et al. | A deep learning method for denoising based on a fast and flexible convolutional neural network | |
CN104849757B (zh) | 消除地震信号中随机噪声系统及方法 | |
CN109669213A (zh) | 基于优化Morlet小波的分频扩散滤波断层强化方法 | |
CN104122584B (zh) | 根据地震数据确定方向性的方法及装置 | |
Xue et al. | Airborne electromagnetic data denoising based on dictionary learning | |
CN110146923A (zh) | 一种高效的高精度深度域地震子波提取方法 | |
CN108919355A (zh) | 基于结构张量导引的高维s变换方法 | |
Wu et al. | Attenuating seismic noise via incoherent dictionary learning | |
Yao-Jun et al. | Unsupervised seismic facies analysis using sparse representation spectral clustering | |
CN109738852A (zh) | 基于低秩矩阵重建的分布式源二维空间谱估计方法 | |
CN107356963B (zh) | 一种数据驱动的自适应的地震信号相干体属性分析方法 | |
CN109471164A (zh) | 基于Ho-RPCA的地震断层增强方法 | |
CN112213782B (zh) | 分相位地震数据的处理方法、装置和服务器 | |
CN109856672B (zh) | 基于深度波数谱的瞬变波包提取方法、存储介质与终端 | |
CN110489480A (zh) | 一种测井数据多属性快速切换可视化方法 | |
Xin et al. | Seismic high-resolution processing method based on spectral simulation and total variation regularization constraints | |
Liu et al. | Enhanced coherence using principal component analysis | |
CN108490414A (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 |