[发明专利]基于Synchrosqueezing变换的地震资料时频分析和衰减估计方法

摘要

本发明公开了一种基于Synchrosqueezing变换的地震资料时频分析和衰减估计方法，首次将新的时频分析工具Synchrosqueezing变换用于地震资料时频分析，该变换通过对变换域系数的重排，获得一个更加聚集的时频表示，时频分辨率大大提高，将其用于实际地震资料分析和致密砂岩模型含气性检测，能够准确界定储层的位置，指示河道与断层等地质结构，进而有利于进一步的资料解释和井位确定；提出基于Synchrosqueezing变换的地震衰减估计方法，并给出具体实现流程，对某油田致密砂岩储层三维数据体的衰减估计结果和钻井结果有着较好的一致性，该方法可以帮助地质人员指示含气储层，确定钻井位置。

主权项

基于Synchrosqueezing变换的地震资料时频分析方法，其特征在于：Synchrosqueezing变换具体包括以下步骤：步骤a：连续小波变换信号的小波变换时域和频域分别表示为：Wf(a,b)=1a∫-∞+∞f(t)ψ*(t-ba)dt=a2π∫-∞+∞F(ω)Ψ*(aω)eibωdω,---(1)]]>其中ψ(t)为基本小波，a和b分别为尺度因子和平移因子，F(ω)和Ψ(ω)为f(t)和ψ(t)的Fourier变换；假设小波函数几乎没有负频率分量，即当ω<0时，Ψ(ω)≈0；步骤b：FM解调频率的计算ωf(a,b)=-i∂bWf(a,b)Wf(a,b)|Wf(a,b)|>0,∞|Wf(a,b)|=0.---(2);]]>步骤c：时间‑尺度域到时间‑频率域的映射(1)连续形式Tf(ω,b)=∫{a:a>0,ωf(a,b)=ω,Wf(a,b)≠0}Wf(a,b)a-3/2da---(3)]]>通过公式(3)，在时间‑尺度域所有和频率ω对应的小波系数进行组合，在时间‑频率域重新将能量“挤压”到频率ω所在的位置；(2)离散形式进行数值计算时，需要对公式(3)中的尺度a和频率ω进行离散；离散化后的尺度记为{ak}，ak表示第k个尺度，其中ak>ak‑1，尺度间隔为ak‑ak‑1＝(Δa)k；对频率进行划分，记为{ωl}，ωl表示第l个频率，其中ωl>ωl‑1，频率间隔为ωl‑ωl‑1＝Δω；Synchrosqueezing变换的离散形式表示为：Tf(ωl,b)=Σak:|ωf(ak,b)-ωl|≤Δω/2Wf(ak,b)ak-3/2(Δa)k,---(4);]]>地震资料时频分析方法具体包括以下步骤：步骤1：在三维数据体地震剖面中选取典型道，对该道数据进行Synchrosqueezing变换，找出异常区域对应频率；步骤2：对整个地震剖面进行Synchrosqueezing变换，提取异常区域对应频率切片；步骤3：对整个三维数据体进行Synchrosqueezing变换，得到频率数据体，然后提取一个沿层切片，供地质人员进行地震资料解释。