CN104199091B - 基于方位地震数据的纹理属性方法 - Google Patents
基于方位地震数据的纹理属性方法 Download PDFInfo
- Publication number
- CN104199091B CN104199091B CN201410424869.9A CN201410424869A CN104199091B CN 104199091 B CN104199091 B CN 104199091B CN 201410424869 A CN201410424869 A CN 201410424869A CN 104199091 B CN104199091 B CN 104199091B
- Authority
- CN
- China
- Prior art keywords
- azimuth
- wide
- data
- represent
- azimuthal
- 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.)
- Expired - Fee Related
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明的基于方位地震数据的纹理属性方法,首先对宽方位角三维地震数据体进行梯度计算,得到每个方位角数据体的三个方向的梯度数据体;其次,根据不同方位角数据体所反应的地层下的地质结构的差异性,利用上一步得到的各个方位角数据体的三个方向的梯度数据体来计算目标点的相关性矩阵;最后对相关性矩阵进行特征值分解,并由得到的特征值计算基于宽方位角数据体的纹理属性体。
Description
技术领域
本发明属于地质领域,具体涉及一种方位纹理属性。
背景技术
纹理是按一定的规则或者一定的形式变化而产生的图案。在自然界中,许多植物、动物、矿物都有独特的纹理特征。利用这些纹理特征,人们可以方便的识别和区分事物。地层由于受到构造运动的影响,也会产生断层、裂缝等地质现象,从而形成地底下的纹理。在进行地质勘探中,我们可以得到一系列的由地震数据形成的图像,例如地震数据剖面图。对于这些地震图像,我们可以借助图像处理的方法、手段来计算三维地震数据的纹理属性,以便凸显地底下的内部特征,从而有效地识别断层或裂缝,为找油、找气提供更好的依据。
在如今的地质领域中,纹理属性的提取都是针对常规的叠后数据来处理的,它们以叠后地震数据体为基础输入,能够通过纹理属性的计算来预测地层中尺度较大、数量较少的断裂,但对较小尺度、数量更庞大的小型或微断裂***难于预测,其根本原因在于叠后地震资料的信息量比较小,缺乏偏移距信息和方位角信息,这样就失去了以各向异性来检测微断裂的理论依据和基础。
发明内容
针对上述问题,本发明提出一种基于方位地震数据的纹理属性方法。
本发明的技术方案为:基于方位地震数据的纹理属性方法,具体包括以下步骤:
S1:计算宽方位角地震数据的梯度;
S2:根据步骤S1得到的梯度,获得各个方位角的三个方向的三维梯度数据体,在得到三维梯度数据体的基础上,获得各个方位角的三维地震数据体在x、y和z方向上的相关性矩阵T;
S3:对相关性矩阵T进行特征值分解:Tv=λv;
其中,λ表示对角矩阵,对角线上的值为特征值,分别为λ1,λ2,λ3,它们是由大到小排列的,即λi≥λi+1,v表示特征向量矩阵,列向量表示所述特征值λ1,λ2,λ3对应的特征向量,分别为v1,v2,v3;
S4:计算基于宽方位角的方位纹理属性体;
所述宽方位角的纹理属性包括:
(1)宽方位角连续属性Ccon:
宽方位角连续属性Ccon的取值范围为0到1;
(2)宽方位角平面属性Cplane:
(3)宽方位角线条属性Cline:
所述宽方位角线条属性Cline的取值范围为0到1;
(4)宽方位角断裂***属性Cfault:
所述宽方位角断裂***属性Cfault的取值范围为0到1;
(5)宽方位角混沌纹理属性Cchaos:
(6)宽方位角倾角属性dip:
其中,atan(·)表示求反正切值,v1(1)表示特征向量v1的第一个值,v1(2)表示特征向量v1的第二个值,v1(3)表示特征向量v1的第三个值;
(7)宽方位角方位属性azimuth:
azimuth=atan2(v1(2),v1(1));
(8)宽方位角相干属性coherence:
其中,Ts(i,j)表示特征分解后的各个方位角的三维地震数据在各个方向上的相关系数。
本发明从宽方位角地震资料出发,以宽方位角三维地震资料的波场特征为基础,利用振幅随方位角变化的特征,充分利用宽方位角三维地震资料的充足信息量(如偏移距信息和方位角信息等)来预测尺度较小、数量更庞大的小型或者微断裂***,从而得到更加稳定、鲁棒性更好的方位纹理属性体。
附图说明
图1是基于宽方位地震数据的纹理属性体的计算流程图。
图2是宽方位地震数据的梯度计算示意图。
图3是宽方位地震数据的梯度计算流程图。
图4是宽方位地震数据的相关性矩阵的计算流程图。
图5是宽方位角数据的Coherence属性。
图6是常规叠后数据的Coherence属性。
图7是常规叠后数据的Dip属性。
图8是宽方位角数据的Dip属性。
图9是宽方位角数据的Dip属性与地质专家提供的溶孔、溶洞分布预测叠合图。
图10是产气井与断层关系图。
图11是常规叠后数据的Continuety属性。
图12是宽方位角数据的Continuety属性。
图13是宽方位角数据的Continuety属性与构造解释断层叠合图。
具体实施方式
总流程图如图1所示,具体方案如下:
S1:首先需要计算输入宽方位角地震数据的梯度,其目的是需要用梯度信息来反映随着方向的不同而检测到的振幅的变化,进而来获得宽方位角的地震纹理属性。其具体思路如下:对每个宽方位角的三维地震数据体,分别在x、y和z方向进行求导运算,从而获得三个方向新的数据体。具体运算采用地震数据与三维差分高斯滤波器作卷积的方式来实现,其目的是降低噪声对求导过程的影响。经典的三维高斯滤波器的表达式为:
其中,xj,yj和zj分别表示沿着x,y和z轴的以待求导的分析点中心,σ为方差或者这里称为尺度因子。这里用Sθ(x,y,z)表示各个方位角的三维原始振幅数据体,其中θ∈{0°,30°,60°,90°,120°,150°}。每个角度的三个方向的梯度数据体的计算表达式如下:
其中,*为卷积运算符号,分别为x、y和z方向的三维差分高斯滤波器,G(x,θ)、G(y,θ)和G(z,θ)即是所需要求得的x,y和z方向的三维梯度数据体,它们可以反映在x,y和z方向上原始的三维方位角振幅数据体中振幅值的变化情况。其具体示意图如图2所示:
计算宽方位角地震数据的梯度的具体流程图如图3所示:
S2:本发明在获得各个方位角的三个方向的三维梯度数据体后,紧接着在梯度结构的基础上,获得各个方位角的三维地震数据体在x、y和z方向上的相关性矩阵,该相关性矩阵即可反映出各个方位角数据体在x、y和z方向的差异,以便通过方位角数据体的纹理属性提的计算来预测地层下的地质结构。相关性矩阵是是一个3*3的对称矩阵,其定义式表示如下:
T代表在计算单点时所得到的相关性矩阵。它的矩阵因子具体表达式如公式2-2-6,其中θ∈{0°,30°,60°,90°,120°,150°}。Tij是通过在时间上在各个角度的三个方向梯度数据体内取窗口长度为Nt方式来计算获得,分别代表各个方向上的相关系数,反映的是各个方位角的三维地震数据在各个方向上的相似程度。在计算目标点的相关性矩阵的过程中,既利用了单个数据体的三个方向的数据信息,又利用了不同方位角数据体的信息,成功的利用了整个四维数据体给出的所有信息,以便后面计算出的方位角的纹理属性体更加的准确、稳定。
求相关矩阵的具体流程图如图4所示:
S3:对2-2-5中的相关性矩阵T如下的特征值分解:
Tv=λv (2-2-7)
其中,λ表示对角矩阵,对角线上的值为特征值,分别为λ1,λ2,λ3,它们是由大到小排列的,即λi≥λi+1。v表示特征向量矩阵,列向量表示以上特征值λ1,λ2,λ3对应的特征向量,分别为v1,v2,v3。
对相关性矩阵做特征值分解的目的是得到宽方位角资料中从各个角度的三维地震数据体中所反应的目标点周围的地质结构。
当λ1=λ2=λ3=0时,说明从各个角度的三维地震数据中测得的目标点周围地质结构为没有差异,振幅强度一致。
当λ1>λ2=λ3=0时,说明从各个角度的三维地震数据中测得的目标点周围地质结构为平面线性结构,在某一个方向上有变化。
当λ1>λ2>λ3=0时,说明从各个角度的三维地震数据中测得的目标点周围地质结构为线性结构。但是当λ1=λ2时,说明该线性结构是各向同性的。
当λ1>λ2>λ3>0时,说明从各个角度的三维地震数据中测得的目标点周围地质结构为非线性的。但是当λ1=λ2=λ3时,说明在目标点周围的结构为各向同性的。
S4:计算基于宽方位角的方位纹理属性体。宽方位角的纹理属性有很多,下面给出主要的几种:
(1)宽方位角连续属性
宽方位角连续性属性的取值范围为0到1,值越大,说明在目标点周围各个角度的数据体所获得的地震数据体的连续性越好,而值越小的地方则可能是断层、裂缝所在的地方。
(2)宽方位角平面属性
(3)宽方位角线条属性
宽方位角线条属性的取值范围为0到1,值越大,说明在目标点周围从各个角度的数据体中所获得的地震数据体的线性越强,值越小,则说明目标点周围并不呈现出线条状,则可能与周围的点相比变化不大。
(4)宽方位角断裂***属性
宽方位角断裂***属性的取值范围为0到1,值越大,说明在目标点周围从各个角度的数据体中所获得的地震数据体的断裂性越强,断层和裂缝越明显,值越小,则目标点周围的地质结构为连续结构。
(5)宽方位角混沌纹理属性
(6)宽方位角倾角属性
(7)宽方位角方位属性
(8)宽方位角相干属性
通过以上几种主要的宽方位角纹理属性体的计算,地层下不同的地质结构可以通过不同的属性体的计算来预测。例如溶孔、溶洞等地质构造可以在宽方位角的纹理属性中的Dip属性体中得到体现,如图9。而在宽方位角的纹理属性中的Continuety属性体中则能清晰明显的看到断层等地质结构,如图13。
下面是基于宽方位角的三维地震数据体的纹理属性的提取在四川某工区的实际应用,并与该工区的传统叠后数据体的地震纹理属性体做了对比,从结果可以看到,与传统叠后数据的地震纹理属性相比,方位地震纹理属性对裂缝和断层的预测更加明显清晰。
图9是由地质专家所绘制的纹理属性中的Dip属性图与溶孔、溶洞分布预测叠合图,可以看到,我们所计算的基于宽方位角地震资料的Dip属性与专家所预测的溶孔、溶洞位置几乎完全吻合,断层也非常清晰明显。
图10是由地质专家所绘制的产气井与断层的关系示意图,其中的线条代表断层,是由基于宽方位角地震资料的纹理属性中的Dip属性得到的。
从上述实施例可知,本发明提供的基于方位地震数据的纹理属性方法,从宽方位角地震资料出发,以宽方位角三维地震资料的波场特征为基础,利用振幅随方位角变化的特征,充分利用宽方位角三维地震资料的充足信息量(如偏移距信息和方位角信息等)来预测尺度较小、数量更庞大的小型或者微断裂***,从而得到更加稳定、鲁棒性更好的方位纹理属性体。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。
Claims (5)
1.基于方位地震数据的纹理属性方法,其特征在于,包括:
S1:计算宽方位角地震数据的梯度;
S2:根据步骤S1得到的梯度,获得各个方位角的三个方向的三维梯度数据体,在得到三维梯度数据体的基础上,获得各个方位角的三维地震数据体在x、y和z方向上的相关性矩阵T;
S3:对相关性矩阵T进行特征值分解:Tv=λv;
其中,λ表示对角矩阵,对角线上的值为特征值,分别为λ1,λ2,λ3,它们是由大到小排列的,即λi≥λi+1,v表示特征向量矩阵,列向量表示所述特征值λ1,λ2,λ3对应的特征向量,分别为v1,v2,v3;
S4:计算基于宽方位角的方位纹理属性体;
所述宽方位角的纹理属性包括:
(1)宽方位角连续属性Ccon:
宽方位角连续属性Ccon的取值范围为0到1;
(2)宽方位角平面属性Cplane:
(3)宽方位角线条属性Cline:
所述宽方位角线条属性Cline的取值范围为0到1;
(4)宽方位角断裂***属性Cfault:
所述宽方位角断裂***属性Cfault的取值范围为0到1;
(5)宽方位角混沌纹理属性Cchaos:
(6)宽方位角倾角属性dip:
其中,atan(·)表示求反正切值,v1(1)表示特征向量v1的第一个值,v1(2)表示特征向量v1的第二个值,v1(3)表示特征向量v1的第三个值;
(7)宽方位角方位属性azimuth:
azimuth=atan2(v1(2),v1(1));
(8)宽方位角相干属性coherence:
其中,Ts(i,j)表示特征分解后的各个方位角的三维地震数据在各个方向上的相关系数。
2.根据权利要求1基于方位地震数据的纹理属性方法,其特征在于,所述S1:计算宽方位角地震数据的梯度,包括:对每个宽方位角的三维地震数据体,分别在x、y和z方向进行求导运算,获得三个方向新的数据体,具体运算公式为:
其中,xj,yj和zj分别表示沿着x,y和z轴的以待求导的分析点中心,σ为方差即尺度因子;
每个角度的三个方向的梯度数据体的计算表达式如下:
其中,Sθ(x,y,z)表示各个方位角的三维原始振幅数据体,θ∈{0°,30°,60°,90°,120°,150°},*为卷积运算符号,和分别为x、y和z方向的三维差分高斯滤波器,G(x,θ)、G(y,θ)和G(z,θ)即是所需要求得的x,y和z方向的三维梯度数据体。
3.根据权利要求1基于方位地震数据的纹理属性方法,其特征在于,所述相关性矩阵T它的矩阵因子具体表达式为,
其中,θ∈{0°,30°,60°,90°,120°,150°},G(x,θ)(x,y,t)、G(y,θ)(x,y,t)、G(z,θ)(x,y,t)表示梯度数据体,Tij分别代表各个方向上的相关系数。
4.根据权利要求1基于方位地震数据的纹理属性方法,其特征在于,所述Tij通过在时间上在各个角度的三个方向梯度数据体内取窗口长度为Nt方式计算获得。
5.根据权利要求1基于方位地震数据的纹理属性方法,其特征在于,所述特征值λ1,λ2,λ3的大小关系,反应的目标点周围的地质结构;
当λ1=λ2=λ3=0时,表示目标点周围地质结构为没有差异,振幅强度一致;
当λ1>λ2=λ3=0时,表示目标点周围地质结构为平面线性结构,在一个方向上有变化;
当λ1>λ2>λ3=0时,表示目标点周围地质结构为线性结构,同时当λ1=λ2>λ3=0时,表示该线性结构是各向同性的;
当λ1>λ2>λ3>0时,表示目标点周围地质结构为非线性的,同时当λ1=λ2=λ3>0时,表示在目标点周围的结构为各向同性的。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410424869.9A CN104199091B (zh) | 2014-08-25 | 2014-08-25 | 基于方位地震数据的纹理属性方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410424869.9A CN104199091B (zh) | 2014-08-25 | 2014-08-25 | 基于方位地震数据的纹理属性方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104199091A CN104199091A (zh) | 2014-12-10 |
CN104199091B true CN104199091B (zh) | 2017-03-01 |
Family
ID=52084403
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410424869.9A Expired - Fee Related CN104199091B (zh) | 2014-08-25 | 2014-08-25 | 基于方位地震数据的纹理属性方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104199091B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021216100A1 (en) * | 2020-04-23 | 2021-10-28 | Saudi Arabian Oil Company | Method and system for sub-seismic fracture prediction |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104820237B (zh) * | 2015-04-30 | 2017-04-05 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 计算地层的方差体的方法 |
CN105242308B (zh) * | 2015-10-14 | 2017-07-07 | 中国石油天然气集团公司 | 一种宽方位地震数据的叠加方法和装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1926558A (zh) * | 2003-08-06 | 2007-03-07 | 兰德马克绘图公司 | 对任意成三角形的表面应用精确三维体纹理的***和方法 |
CN102903141A (zh) * | 2012-08-31 | 2013-01-30 | 电子科技大学 | 基于不透明度加权的多地震属性融合纹理映射体绘制方法 |
CN103364833A (zh) * | 2013-07-01 | 2013-10-23 | 西安交通大学 | 一种高精度地层倾角估计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100669251B1 (ko) * | 2005-11-25 | 2007-01-16 | 한국전자통신연구원 | 디지털 영상 품질 자동 분석 장치 및 방법 |
-
2014
- 2014-08-25 CN CN201410424869.9A patent/CN104199091B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1926558A (zh) * | 2003-08-06 | 2007-03-07 | 兰德马克绘图公司 | 对任意成三角形的表面应用精确三维体纹理的***和方法 |
CN102903141A (zh) * | 2012-08-31 | 2013-01-30 | 电子科技大学 | 基于不透明度加权的多地震属性融合纹理映射体绘制方法 |
CN103364833A (zh) * | 2013-07-01 | 2013-10-23 | 西安交通大学 | 一种高精度地层倾角估计方法 |
Non-Patent Citations (1)
Title |
---|
基于体平滑的体绘制纹理伪边界消除算法;程尧;《网络新媒体技术》;20120731;第1卷(第4期);全文 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021216100A1 (en) * | 2020-04-23 | 2021-10-28 | Saudi Arabian Oil Company | Method and system for sub-seismic fracture prediction |
Also Published As
Publication number | Publication date |
---|---|
CN104199091A (zh) | 2014-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Liu et al. | Common-azimuth seismic data fault analysis using residual UNet | |
US10145227B2 (en) | Method for estimating permeability of fractured rock formations from induced slow fluid pressure waves | |
US20180203144A1 (en) | Interferometric Microseismic Imaging Methods and Apparatus | |
CN104614782B (zh) | 一种盐湖卤水探测方法 | |
US20140129479A1 (en) | Method to aid in the exploration, mine design, evaluation and/or extraction of metalliferous mineral and/or diamond deposits | |
CN107991711B (zh) | 航空时域电磁三维条状随机断裂带模型建立及判别方法 | |
Gao et al. | Extreme curvature and extreme flexure analysis for fracture characterization from 3D seismic data: New analytical algorithms and geologic implications | |
CN104199091B (zh) | 基于方位地震数据的纹理属性方法 | |
CN104199088B (zh) | 一种提取入射角道集的方法及*** | |
Bruno et al. | Groundwater characterization in arid regions using seismic and gravity attributes: Al Jaww Plain, UAE | |
CN104216013A (zh) | 基于宽方位角资料的c3相干体的方法 | |
CN113779811B (zh) | 一种构造变形对现今地应力方位扰动的定量分析方法 | |
Chen et al. | The algorithm of 3D multi-scale volumetric curvature and its application | |
CN104316960A (zh) | 一种基于vsp的储层裂缝预测方法及*** | |
CN104749618B (zh) | 泥页岩缓倾角裂缝叠后概率定量表征方法 | |
CN104820237B (zh) | 计算地层的方差体的方法 | |
Wang et al. | A Three-Dimensional Geological Structure Modeling Framework and Its Application in Machine Learning | |
CN107797148B (zh) | 一种基于三维地质建模的航磁异常场分离方法及*** | |
CN113204855B (zh) | 地层缩短速率确定方法、***、设备及可读介质 | |
Zhdanov et al. | Principles of tensor induction well logging in a deviated well in an anisotropic medium | |
CN105093330B (zh) | 多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法 | |
CN105093309A (zh) | 基于构造应力分析进行裂缝检测及有效性验证的方法 | |
Cheng et al. | Multiscale fracture prediction technique via deep learning, seismic gradient disorder, and aberrance: Applied to tight sandstone reservoirs in the Hutubi block, southern Junggar Basin | |
Wang et al. | Fracture extraction from logging image using a dual encoder-decoder architecture with Swin transformer | |
Chopra et al. | Using automatically generated 3D rose diagrams for correlation of seismic fracture lineaments with similar lineaments from attributes and well log data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170301 Termination date: 20190825 |
|
CF01 | Termination of patent right due to non-payment of annual fee |