基于混合倾角扫描振幅变化率的潜山储层地震识别方法与流程

技术特征：

1.基于混合倾角扫描振幅变化率的潜山储层地震识别方法，其特征在于，该基于混合倾角扫描振幅变化率的潜山储层地震识别方法包括：

步骤1，在时间域应用多道相干算法进行倾角搜索估计地层倾角值；

步骤2，在频率域应用地层倾角在地震记录上表现的时间延迟，对其进行傅里叶变换的时延特性来估计地层倾角；

步骤3，基于时间-频率域混合倾角扫描，沿地层倾角在局部层拉平的基础上计算时窗内平均绝对振幅的空间变化率；

步骤4，根据平均绝对振幅空间变化率数据体，通过连井剖面和沿层切片进行潜山内幕孔缝型储层的地震识别方法研究。

2.根据权利要求1所述的基于混合倾角扫描振幅变化率的潜山储层地震识别方法，其特征在于，该基于混合倾角扫描振幅变化率的潜山储层地震识别方法还包括，在步骤1之前，采集原始三维地震数据。

3.根据权利要求1所述的基于混合倾角扫描振幅变化率的潜山储层地震识别方法，其特征在于，在步骤2中，在道地震数据中，假设一个t时刻的地震平面波，令p和q为x轴和y轴方向上的视倾角，将(x_j,y_j,t)为中心点，则上述时窗内的平均道的能量u和全部道能量U比为相似系数σ(p,q,t)由公式给出：

$\mathrm{\σ}(p,q,t)=\frac{{\[\frac{1}{J}\underset{j=1}{\overset{J}{\Σ}}u(x_j,y_j,t-{\mathrm{px}}_j-{\mathrm{qy}}_j)\]}^2}{\frac{1}{J}\underset{j=1}{\overset{J}{\Σ}}{\[u(x_j,y_j,t-{\mathrm{px}}_j-{\mathrm{qy}}_j)\]}^2}$

j表示道数，(x_j,y_j)表示该点在坐标轴上的坐标，利用上述公式通过倾角扫描方式计算出视倾角。

4.根据权利要求1所述的基于混合倾角扫描振幅变化率的潜山储层地震识别方法，其特征在于，在步骤3中，假设在地面上接收到的某一道地震记录为x(t),与它相邻的地震记录的时间延迟为t₀，则相邻地震记录道y(t)＝x(t-to),对x(t)的傅里叶变换由公式：

$X\left(\mathrm{\ω}\right)={\∫}_{-\mathrm{\∞}}^{+\mathrm{\∞}}x\left(t\right)e^{-i\mathrm{\ω}t}dt$

再依据傅里叶变换的时延特性公式：

$Y\left(\mathrm{\ω}\right)={\∫}_{-\mathrm{\∞}}^{+\mathrm{\∞}}x(t-t_0)e^{-i\mathrm{\ω}t}dt=e^{-{\mathrm{i\ωt}}_0}X\left(\mathrm{\ω}\right)$

求对应于不同角频率的时间延迟公式：

$t_{\mathrm{\ω}}=-\frac{1}{i\mathrm{\ω}}\ln \[\frac{Y\left(\mathrm{\ω}\right)}{X\left(\mathrm{\ω}\right)}\]$

上式中，X(ω)、Y(ω)为相邻地震记录的傅里叶变换，ω为角频率，i为虚数单位；

先求取地震波有效频段对应的时间延迟值的算术平均值，并将算术平均值与空间采样间隔的比值作为估计的地层倾角值，设有效频段的起始值、终止值分别为ω_start，ω_end，选择地震记录主频附近的频段，其中包含N个角频率值，则这两道地震记录之间的时间延迟t₀公式：

$t_0=\underset{\mathrm{\ω}={\mathrm{\ω}}_{start}}{\overset{{\mathrm{\ω}}_{end}}{\Σ}}{t}_{\mathrm{\ω}}/N$

即上式表明，最终得到的时间延迟值是各个不同频率下对应的时延值的平均值。

5.根据权利要求1所述的基于混合倾角扫描振幅变化率的潜山储层地震识别方法，其特征在于，在步骤4中，平均绝对振幅空间变化率的数学表达式为：

$AVR=\sqrt{{\left(\underset{t}{\Σ}\right|\frac{dA(x,y,t-{\mathrm{\τ}}_x-{\mathrm{\τ}}_y)}{dx}\left|\right)}^2+{\left(\underset{t}{\Σ}\right|\frac{dA(x,y,t-{\mathrm{\τ}}_x-{\mathrm{\τ}}_y)}{dy}\left|\right)}^2}$

其中AVR表示平均绝对振幅空间变化率，dA/dx表示振幅在x方向的变化率，dA/dy表示振幅在y方向的变化率，τ_x、τ_y分别表示相邻地震道沿地层倾角在x和y方向上的时差，对于任意x和y，通过提取t-τ_x-τ_y时刻的振幅值达到局部层拉平的效果；

计算得到的AVR所表示的就是中心点振幅相对于其周边振幅的变化强度，从而去除了地层的背景值，使得孔缝型储层与围岩之间的反射强度差异能够更加直观地反映出来。