[发明专利]一种地震成果数据吸收补偿定量分析评价方法

权利要求书

1.一种地震成果数据吸收补偿定量分析评价方法，其特征在于：

步骤一、计算地震成果数据的振幅谱系数a；

对地震成果数据吸收补偿前后的同一时窗作归一化振幅频谱分析，基于频谱形态、高低频展布特征、主频与带宽信息，从频谱的角度去判定吸收补偿的合理性，a为主频移动值f与带宽band的乘积，即a＝f×band，f是地震成果数据吸收补偿前后主频移动值，band是吸收补偿后的带宽，band为振幅谱纵坐标为0.2时的高频与低频差值；

步骤二、计算地震成果数据的能量衰减斜率b；

对地震道从浅至深采用广义S变换作时频分析，时频谱横向上分析高频能量的补偿情况，纵向上分析中深层能量的恢复情况，b定义为能量衰减斜率，定义地震道长度为T秒，采样间隔为sam,则采样点数n＝T/sam，对于某地震道第u个采样点,u∈[0,n),它的时间深度表示为t_u＝u×sam，它对应的能量衰减斜率表示为该采样点对应的时频谱主频F_u与后一个相邻采样点对应的时频谱主频F_u+1之差，记为b_u＝F_u+1-F_u，则该地震道的能量衰减斜率表示为地震数据无吸收衰减时b＝1，完全吸收衰减时b＝0，欠补偿时为负值，其他情况取正值；

步骤三、计算地震成果数据的相关度系数c；

统计地震成果数据浅层的主频为f_main，利用主频为f_main的雷克子波生成人工合成地震记录，将吸收补偿后的地震成果数据与该人工合成地震记录作互相关计算，求取相关度系数，记为c；

步骤四、计算地震成果数据的信噪比SNR；

假设信噪比计算时窗表示为[w_ij]_M×N,在该时窗内共有N道地震数据，每道地震数据有M个地震采样点，i∈M,j∈N,w_ij表示第j道上的第i点对应的振幅数值，因为地震数据具有同相性，且噪音是随机的，设有效信号表示为signal，噪音表示为noise，则有

w_ij＝signal_ij+noise_ij

有效信号signal_i是N道地震道水平叠加的结果，于是N道有效信号的总能量表示为：

N道信号的总能量表示为

并有E＝E_signal+E_noise

由此得到信噪比SNR为

步骤五、计算地震成果数据吸收补偿评价值ε，通过评价值ε的大小，评价吸收补偿的合理性，若ε0时，表示吸收补偿为欠补偿，存在补偿不足情况，当0ε0.5时，表示吸收补偿为过补偿，存在补偿过度情况，当ε0.5或等于0.5时，表示吸收补偿为合理的；

ε值的计算方法：将步骤一到步骤四计算得到的a,b,c与SNR值代入下面的吸收补偿评价值计算公式

ε＝SNR×(0.3a+0.4abs(b)+0.3c)/max(a,b,c)×sign(b)

其中，abs(b)是求取b的绝对值运算，max(a,b,c)是求取a,b,c中最大的绝对值，sign(b)是取符号运算，如果b是负数，则sign(b)＝-1,否则sign(b)＝1。