地震属性提取与分析-配合软件介绍

相干体属性（Coherence）•GeoEast开发了三种相干算法：相关算法Corr（C1）相似算法Semb（C2）特征值算法EIg（C3）。算法特点：相关算法Corr：运算速度最快，但是其抗噪能力最差，建议用于信噪比较高，以识别大断层为主的地震资料中。相似算法Semb：比Corr的运算速度稍慢，但其抗噪能力较强，推荐多数情况下使用该方法。特征值算法Eig：运算速度最慢，但是该方法在保证抗噪能力的同时能提高断层的横向分辨能力，在选择3点及5点提取相干时，Eig与Semb方法相当，但是在选择9点提取相干时，Eig比Semb算法对断层的横向分辨率更高，对小断层更敏感。参数说明：•MaxDip为最大倾角扫描时窗（0-15ms）该值过小会使层位信息带入相干属性中，过大会丢失断层信息，最佳值可根据实际地震资料的层位最大倾角（最陡同相轴的相邻道时间延迟量）来确定。•Window为相干时窗（20-200ms）：过小使相干体压制噪音能力下降，过大使断层分辨能力降低，合适的时窗为地震波形视周期的1.5~2倍，即地震信号的全波形长度的1.5~2倍。22相干体属性（Coherence）参数说明：3点5点9点多窗口倾角扫描（ScanDipWithMultiwindow）1.在三维地震资料解释中，地震反射层位的倾角和方位角是非常重要的信息，目前主要应用是用来对局部反射面的不连续性的估计。2.通常，利用垂直窗口进行倾角和方位角估算比在拾取的层位进行估算能提供更为稳定的估算结果。倾角和方位角体是构造导向滤波、体曲率和相干能量梯度的基础。多窗口倾角扫描（ScanDipWithMultiwindow）(Marfurt,2006)（1）最相干的反射层倾角离散扫描(瞬时倾角=具有最大相干的倾角)AnalysisPointMinimumdiptested(-200)Maximumdiptested(+200)Dipwithmaximumcoherence(+50)（2）包含分析点的最佳相干窗口搜索④①②③⑤⑥⑦⑧⑨5X5的网格crosslineinline①②③④⑤⑥⑦⑧⑨通过多窗口扫描，求取相似程度最大的窗口作为分析点的倾角、方位角估算窗口，提高了分析点的倾角和方位角精度。多窗口倾角扫描（ScanDipWithMultiwindow）(Marfurt,2006)单窗口扫描多窗口扫描Time(s)CLR多窗口倾角扫描（ScanDipWithMultiwindow）多窗口扫描界面输出数据体类型①DipAlongInline:inline方向的视倾角②DipAlongCrossline:crossline方向的视倾角③KuwaharaWindow:kuwahara扫描窗口④DipMagnitude:倾角标量⑤DipAzimuth:倾向方位角⑥DipConfidence:视倾角可信度⑦DipofAnyAzimuth:任意方位角视倾角参数说明①WindowLength(ms):扫描时窗长度(ms)，大于等于Samplinginterval的5倍。②WindowsNumber:扫描窗口数目，有1、5、9三个可选项，缺省值9。③MaxScanDip(o):最大扫描倾角(o)，范围[0,90]，缺省值30。④DipIncrement(o):倾角扫描角度增量(o)，范围[1,10]，缺省值5。⑤Ref.Velocity(m/s):反射层速度(m/s)，范围[2000,3500]，缺省值3000。⑥BeginAzimuth(o):起始方位角(o)，范围[0,180)，缺省值0。⑦AzimuthIncrement(o):方位角增量(o)，范围[0,90)，缺省值30。⑧AzimuthNumber:方位角数目，范围[1,6]。多窗口倾角扫描（ScanDipWithMultiwindow）多窗口倾角扫描计算倾角方位角在四川剑阁的应用倾角、方位角与相干比较，很好地刻画了多边形断层、河道的特征。河道相干方位角倾角倾角构造导向滤波（StructureOrientedFiltering）SeismicData：地震数据InlineDip：Inline方向的倾角数据体CrosslineDip：Crossline方向的倾角数据体KuwaharaWindow：Kuwahara数据体（记录信息）输入数据：Inline方向的倾角数据体（多窗口倾角扫描的输出结果数据体）构造导向滤波（StructureOrientedFiltering）滤波方法：①PCFilter：主分量滤波②MeanFilter：均值滤波③MedianFilter：中值滤波命名方式：前缀(用户)+固定形式后缀①ABC.SOF_PCFilter②ABC.SOF_MeanFilter③ABC.SOF_MedianFilter构造导向滤波（StructureOrientedFiltering）突出断裂和地层连续性通常相干体的算法是计算地震数据的总振幅的变化,相干能量梯度计算地震数据相干分量的振幅变化,这有利于突出一些细小的地质特征的变化，从而提高解释精度。它与相干属性在解释中的应用主要是：识别断层、裂缝带、河道和砂体边界等。利用相干能量梯度和相干属性，都很好地刻画了河道分布，但还是有细微差别。相干能量梯度（EnergyGradientCalculating）相干能量梯度（EnergyGradientCalculating）InlineEnergyGradient：Inline方向的相干能量梯度CrosslineEnergyGradient：Crossline方向的相干能量梯度GradientMagnitude：相干能量梯度标量GradientAzimuth：相干能量梯度方位EnergyRatioSimilarity：基于能量比率的相似相干OuterProductSimilarity：基于道相关的相干SobelFilterSimilarity：基于sobel滤波的特征值相干AnyAzimuthEnergyGradient：任意方位相干能量梯度CoherentEnergy：相干能量TotalEnergy：总能量相干和能量梯度属性的输出数据：可输出10个数据体体曲率属性（VolumetricCurvature）体曲率属性的输出数据：可输出构造和振幅两种曲率构造类曲率的输入数据是：沿Inline和Crossline的多窗口倾角扫描的倾角数据体振幅类曲率的输入数据是：沿Inline和Crossline的相干能量梯度的能量梯度数据体短波长(0.75)长波长(0.25)中波长(0.5)不同尺度的曲率刻画了不同尺度的地质特征多尺度体曲率属性（VolumetricCurvature）多尺度参数：SpaceFilter滤波参数FractionDerivative分数导数滤波FractionalDerivativePower分数导数指数（取值范围：0.25—2）注意：0.25代表长波长即大尺度曲率；0.50代表中波长即中尺度曲率；0.75代表短波长即小尺度曲率；值越大代表曲率尺度越小。长波长：是大尺度曲率，适用于寻找大断层、大裂缝等；中波长：是中尺度曲率，适用于寻找中等断层、裂缝等；短波长：是小尺度曲率，适用于寻找小断层、裂缝等另外：LambdaMin：Lambda最小值（缺省值：2*dx）dx：工区线道间隔ReferenceVelo：反射层参考速度值（缺省值：3000m/s）Operator_clip：滤波算子振幅截断比例（缺省值：0.01）体曲率属性（VolumetricCurvature）PrincipalCurvatures(k1andk2)：主曲率MaxandMinCurvature：最大、最小曲率Most-PositiveandMost-Negative：最大正曲率、最大负曲率ReflectorConvergenceinInline：沿inline的反射面收敛数ReflectorConvergenceinCrossline：沿crossline的反射面收敛数AngularUnconformity：角度不整合GaussianCurvature：高斯曲率MeanCurvature：平均曲率DipCurvature：倾角曲率StrikeCurvature：走向曲率体曲率属性的输出数据：构造类曲率可输出22个数据体ReflectorRotation：反射面旋转度ShapeMeasurements：形态类曲率属性Curvedness：弯曲度ShapeIndex：形态指数AzimofMin：最小曲率方位角Dome：圆顶型Bowl：碗型Saddle：马鞍型Ridge：山脊型Valley：山谷型体曲率属性的输出数据：构造类曲率可输出22个数据体注：最小曲率方位角（AzimofMin）、山脊型（Ridge）和山谷型（Valley）是玫瑰图的输入数据。体曲率属性（VolumetricCurvature）体曲率属性（VolumetricCurvature）Most-PositiveandMost-Negative：最大正曲率、最大负曲率GaussianCurvature：高斯曲率MeanCurvature：平均曲率DipCurvature：倾角曲率StrikeCurvature：走向曲率ReflectorRotation：反射面旋转度体曲率属性的输出数据：振幅类曲率可输出15个数据体ShapeMeasurements：形态类曲率属性8种Curvedness：弯曲度ShapeIndex：形态指数AzimofMin：最小曲率方位角Dome：圆顶型Bowl：碗型Saddle：马鞍型Ridge：山脊型Valley：山谷型体曲率属性的输出数据：振幅类曲率可输出15个数据体体曲率属性（VolumetricCurvature）玫瑰图（RoseDiagram）玫瑰图上某个方位的花瓣长度依赖于在该方位裂缝出现的频率，即代表裂缝分布的密度；花瓣的方向代表裂缝的方位。玫瑰图最终成为了玫瑰体，其大小与原始数据一致，每个玫瑰分布值均位于原始数据分析时窗的中心位置，这样也就相当于形成了一个裂缝分布的属性体。通过玫瑰体的时间切片可以分析裂缝分布的密度和方位。玫瑰体时间切片玫瑰图（RoseDiagram）Input输入数据：可以选择山脊类曲率分量或山谷类曲率分量和最小曲率方位角LineamentMagnitude：线性构造特征数据体（体曲率的山脊或山谷类数据体）LineamentAzimuth：线性构造特征方位角数据体（体曲率的最小曲率方位角数据体）ProcessingParameters处理参数MaximumRoseRadius：最大玫瑰图直径（单位：米）NumberofPetalsPer90：每90度的花瓣个数（取值范围：1－18，一般为6，即每个花瓣为15度）MagnitudeUpperThreshold：线性构造特征数值的上门槛值MagnitudeLowerThreshold：线性构造特征数值的下门槛值CalculateMagnitudeRange线性构造数据体范围计算，不推荐使用该功能。玫瑰图参数说明：属性比例融合示意图属性数据体比例融合实现把两个地震属性数据体按比例关系融合在一起。一般情况下把背景数据如地震数据作为A，反映地质特征如断层的属性数据（相干等）为B。取值时地震数据A的范围：全选属性B范围：选择反应断层特征的值域范围（剖面上查找）分界值Cm：选择地震数据范围内的一个值：一般取数据体A的75%左右。属性体比例融合（AttributeScaleFuse）输入背景数据A输出数据C输入特征数据B方差（Variance）方差属性体能够准确识别断层及地层不连续变化，作用与相干属性类似。TimeDelayMethod时移方法：相似系数法：SimilarCoefficientPointType所用数据所在线道位置排列方式。共有10种排列方式。WindowLength(ms)滑动时窗长度，单位：ms。DipScan(ms)倾角扫描时窗长度，单位：ms。2pl