[发明专利]一种基于谱分解和相干属性融合流程的水道复合体沉积要素精细识别方法

权利要求书

1.一种基于谱分解和相干属性融合流程的水道复合体沉积要素精细识别方法，它包括以下步骤：

(1)数据体频率分析：根据目标水道复合体发育的位置分析地震资料频带，得到所述地震资料的最低频率、最高频率和主频，并计算得到能识别的沉积要素最大和最小厚度；

(2)形成调谐数据体：结合步骤(1)得到的数据对目标水道复合体进行S变换，S变换的时窗通过沿层开取时窗或在水道顶底层之间进行；通过计算，在垂向上生成频率连续变化的振幅和相位数据体，得到目标水道复合体内沉积要素的调谐频率；然后通过观察不同频率下对应的沉积要素的成像特征，确定不同沉积要素对应的有效频率；

(3)建立离散时频率体：将步骤(1)中所述地震资料分解成一系列时间域离散频率振幅数据体，该数据体在垂向上为时间，每个数据体只包含需要表征的沉积要素对应的频率成分，反映了特定频率振幅的空间变化特征；

(4)进行相干体运算：通过地震上波形和振幅侧向上的变化与断层或者沉积要素的边界进行相干体运算，得到沉积要素的边界信息；

(5)体融合运算：将步骤(3)中得到的不同沉积要素的对应的离散时频体分别与步骤(4)运算的相干体叠加、融合显示，得到不同沉积要素对应下的相干时频体；

(6)应用等时切片或地层切片技术对所述相干视频体进行可视化显示，重构所述目标水道复合体单一沉积要素垂向演化过程和平面分布特征，以精细识别沉积要素。

2.根据权利要求1所述的方法，其特征在于：步骤(1)中，所述沉积要素包括天然堤、点坝、单期水道、决口扇和沉积波；

计算得到能识别的沉积要素最大和最小厚度的公式如式Ⅰ和式Ⅰ-1所示：

H_max＝V/4f_max 式Ⅰ

H_min＝V/4f_min 式Ⅰ-1

式Ⅰ和式Ⅰ-1中，f_max为步骤1中目标水道复合体的最大频率；f_min为步骤1中目标水道复合体的最小频率；V为地震资料速度；H_max和H_min分别为可分辨沉积要素的最大和最小厚度。

3.根据权利要求1或2所述的方法，其特征在于：所述步骤(2)中的S变换是时间和频率的联合函数，以时间和频率为变量来描述信号的能量密度或信号强度，其根据如下式Ⅱ所示公式进行计算：

式Ⅱ中，t表示时间，f表示频率，τ为高斯窗在时间轴上的位置，时窗宽度随频率的变化而改变，频率的倒数决定了时窗的大小，s(t)表示时间域输入信号，exp为自然对数函数，π为圆周率，i表示虚数单位，dt表示积分函数。

4.根据权利要求1-3中任一项所述的方法，其特征在于：步骤(3)中根据步骤(2)中确定的目标水道复合体的沉积要素对应频率进行地震资料分解；

所述步骤3中的时频率体，包括不同厚度对应下的沉积要素频率15HZ、20HZ、25Hz、30Hz、40Hz和50Hz频率体。

5.根据权利要求1-4中任一项所述的方法，其特征在于：所述步骤(4)中，所示相干体的运算采取BG Fast算法，该算法在所示地震数据中倾角定义为水平时间梯度与垂直时间梯度平方和的均方根值，而方位角定义为垂直时间梯度与水平时间梯度的比值的反正切值。

6.根据权利要求1-5中任一项所述的方法，其特征在于：所述步骤(5)中，所述体融合运算通过兰德马克公司的Geoprobe软件阴影叠置功能实现。