一种地震波场能量均衡方法及处理装置

摘要

本发明提供一种地震波场能量均衡方法，包括：S1：处理装置读入经过预处理的地震数据；S2：处理装置对输入的地震数据，用异常振幅消除技术逐炮进行单炮能量均衡；S3：处理装置对单炮能量均衡处理后的地震数据，进行多炮集记录的炮间能量均衡处理。上述能量均衡方法，从能量角度来消除这两个环节的差异，没有涉及到地震波场的相对关系，在进行地震数据能量均衡后，得到的结果是浅层能量强于深层能量、近道能量强于远道能量，保持了波场的相对关系。此外本发明还提供一种地震波场能量均衡处理装置实施例，可以实现保持地震波场特征的能量均衡处理。

主权项

一种地震波场能量均衡方法，其特征在于，包括：S1：处理装置读入经过预处理的地震数据；S2：处理装置对读入的地震数据，用异常振幅消除技术逐炮进行单炮能量均衡处理；S3：处理装置对单炮能量均衡处理后的地震数据，进行多炮集记录的炮间能量均衡处理；其中，所述S2包括：S21：输入炮集中一组地震道x(i,j)，计算复数地震道的虚部y(i,j)；S22：根据地震道x(i,j)和复数地震道虚部y(i,j)，计算地震道的包络；S23：计算地震道包络A(i,j)的中值M(i)；S24：根据中值M(i)，设计检测异常振幅的准则；S25：对经过S24检测的地震道数据应用权值，衰减异常振幅的能量；S26：依次移动道窗，并重复执行S23至S25，直至处理完炮集中的所有地震道数据；所述S3包括：S31：计算每个炮集的总能量；S32：根据S31计算出的炮集总能量，对每个炮集的能量进行归一化；所述S25中衰减异常振幅的能量，具体通过下式进行：$\hat{x}(i,j)=x(i,j)/l(i,j)$(i＝0,1,2...,nt‑1,j＝tracew/2)其中，x(i,j)是地震道；i是时间采样点；nt是道长，nt为定值；j是空间采样点；tracew是窗口内道数，tracew取奇数；是异常振幅消除后的地震道；l(i,j)是权值，权值l(i,j)计算方法如下：$l(i,j)=\left\{\begin{array}{cc}\frac{A(i,j)}{M\left(i\right)·\alpha }& A(i,j)>threshold·M\left(i\right)\\ 1& A(i,j)\le threshold·M\left(i\right)\end{array}\right.$其中，threshold是门槛值，threshold>1；α是衰减系数，0＜α≤1。