一种多次波压制方法及装置制造方法及图纸

本申请实施例公开了一种多次波压制方法及装置。所述方法提供有目的工区的测井数据以及所述目的工区的时间域的叠后地震数据；所述方法包括：根据测井数据中的声波时差测井曲线数据和密度测井曲线数据，确定目的工区的钻井位置处的地层反射系数曲线数据；基于地层反射系数曲线数据，确定待处理地震数据，以及确定时间域的叠后地震数据中目标反射源界面对应的地震子数据；其中，待处理地震数据用于表征时间域的叠后地震数据中包含多次波的地震数据；目标反射源界面表示产生多次波的反射界面；根据地震子数据，对待处理地震数据进行多次波压制处理。本申请实施例提供的技术方案，可以实现对层间多次波的有效压制。

【技术实现步骤摘要】
一种多次波压制方法及装置
本申请涉及石油地球物理勘探地震数据的处理和解释
，特别涉及一种多次波压制方法及装置。
技术介绍
当反射波向上传播遇到某个具有较大波阻抗差的界面时，反射波有可能从这个界面反射，向下传播。当反射波遇到反射界面时，再次发生反射返回地(海)面，于是形成了多次波。按照反射层的位置，多次波可分为以下几种类型：全程多次波、层间多次波、微屈多次波多次波以及鸣震。其中，层间多次波混杂在有效波(即一次波)之间，与一次波相互叠加和干涉，使一次波的波场特征畸变，模糊不清或难以识别。多次波严重时，导致地震剖面出现假的地质现象，降低了数据处理和解释工作的准确性。因此，研究压制多次波的方法具有一定的实际意义和应用价值。多次波问题一直是地震勘探领域中的难点。在陆上勘探中，由于地层浅部结构的不稳定性，使其表面多次波不像海上数据那样明显，而陆上层间多次波则更为复杂，没有较为明显的特征，往往不易正确地判断地震剖面上产生层间多次波的界面位置，因此在数据处理中很难识别和去除。目前，层间多次波压制方法以对叠前地震数据中的层间多次波进行压制为主，但由于层间多次波具有来源可能性多、周期性差以及横向分布范围变化快的特征，导致多次波波场的特征比较复杂，采用现有的方法并不能完全有效地衰减层间多次波。因此，需要一种新的压制多次波方法，来有效地衰减层间多次波。
技术实现思路
本申请实施例的目的是提供一种多次波压制方法及装置，以实现对层间多次波的有效压制。为解决上述技术问题，本申请实施例提供一种多次波压制方法及装置是这样实现的：一种多次波压制方法，提供有目的工区的测井数据以及所述目的工区的时间域的叠后地震数据；所述方法包括：根据所述测井数据中的声波时差测井曲线数据和密度测井曲线数据，确定所述目的工区的钻井位置处的地层反射系数曲线数据；基于所述地层反射系数曲线数据，确定待处理地震数据，以及确定所述时间域的叠后地震数据中目标反射源界面对应的地震子数据；其中，所述待处理地震数据用于表征所述时间域的叠后地震数据中包含多次波的地震数据；所述目标反射源界面表示产生多次波的反射界面；根据所述地震子数据，对所述待处理地震数据进行多次波压制处理。优选方案中，所述根据所述声波时差测井曲线数据和所述密度测井曲线数据，确定所述目的工区的钻井位置处的地层反射系数曲线数据，包括：根据所述声波时差测井曲线数据和所述密度测井曲线数据，确定所述目的工区的钻井位置处的波阻抗曲线数据；根据所述波阻抗曲线数据，确定所述目的工区的钻井位置处的地层反射系数曲线数据。优选方案中，采用下述公式确定所述波阻抗曲线数据：其中，PI表示所述波阻抗曲线数据中指定连续排列的N个采样点对应的波阻抗平均值，其中，所述波阻抗平均值用于分配给所述N个采样点中的任意一个采样点；N表示指定采样点数量，ρi表示所述密度测井曲线数据中指定连续排列的N个采样点中第i个采样点的密度，dti表示所述声波时差测井曲线数据中指定连续排列的N个采样点中第i个采样点的声波时差。优选方案中，采用下述公式确定所述地层反射系数曲线数据：其中，rj表示所述地层反射系数曲线数据中第j个采样点的地层反射系数，PIj和PIj-1分别表示所述波阻抗曲线数据中第j个采样点和第j-1个采样点的波阻抗。优选方案中，所述基于所述地层反射系数曲线数据，确定待处理地震数据，包括：根据指定地震子波数据和所述地层反射系数曲线数据，确定所述目的工区的钻井位置处的深度域的测井合成地震道数据；从所述时间域的叠后地震数据中获取所述目的工区的钻井位置处对应的时间域的井旁地震道数据；基于所述深度域的测井合成地震道数据和所述时间域的井旁地震道数据，确定所述时间域的叠后地震数据中的目标多次波压制时窗范围，并将所述时间域的叠后地震数据中目标多次波压制时窗范围内的地震数据作为所述待处理地震数据。优选方案中，所述基于所述深度域的测井合成地震道数据和所述时间域的井旁地震道数据，确定所述时间域的叠后地震数据中的目标多次波压制时窗范围，包括：将所述深度域的测井合成地震道数据转换为时间域的测井合成地震道数据；在所述时间域的测井合成地震道数据和所述时间域的井旁地震道数据的时窗范围内，移动预设窗长的时窗，并计算所述时间域的测井合成地震道数据和所述时间域的井旁地震道数据分别在当前时窗位置处的测井合成地震道子数据和井旁地震道子数据之间的相关性系数；确定小于预设相关性系数阈值的相关性系数对应的目标时窗位置，并基于所述目标时窗位置构成所述目标多次波压制时窗范围。优选方案中，所述根据所述地层反射系数曲线数据，确定所述时间域的叠后地震数据中目标反射源界面对应的地震子数据，包括：确定所述地层反射系数曲线数据中地层反射系数对应的采样时间；对所述地层反射系数曲线数据中的地层反射系数按照由大到小的顺序进行排序，并将排序后的地层反射系数中前指定数量的地层反射系数作为所述指定数量的目标地层反射系数，并将所述目标地层反射系数对应的采样时间作为所述目标反射源界面在所述钻井位置处的界面时间点；其中，所述地层反射系数与所述目标反射源界面一一对应；将所述指定数量的目标反射源界面在所述钻井位置处的界面时间点分别标定于所述叠后地震数据上，并对标定后的叠后地震数据进行层位追踪和解释，得到所述目标反射源界面的时间层位；将所述时间域的叠后地震数据中与所述目标反射源界面的时间层位相邻的指定时窗范围内的地震数据作为所述目标反射源界面对应的地震子数据。优选方案中，所述根据所述地震子数据，对所述待处理地震数据进行多次波压制处理，包括：根据指定地震子波数据和所述地层反射系数曲线数据，确定所述目的工区的钻井位置处的深度域的测井合成地震道数据；并将所述深度域的测井合成地震道数据转换为时间域的测井合成地震道数据；按照所述指定数量的目标地层反射系数由大到小的排列顺序，依次迭代根据与所述目标地层反射系数相关联的目标反射源界面对应的地震子数据，并对所述待处理地震数据进行多次波压制处理的步骤，直至多次波压制后的待处理地震数据中的井旁地震道子数据与所述测井合成地震道数据在所述目标多次波压制时窗范围内的测井合成地震道子数据之间的相关系数大于或等于预设相关系数阈值；其中，所述多次波压制后的待处理地震数据中的井旁地震子数据表示所述多次波压制后的待处理地震数据中在所述目的工区的钻井位置处对应的井旁地震道子数据。优选方案中，根据与所述目标地层反射系数相关联的目标反射源界面对应的地震子数据，对所述待处理地震数据进行多次波压制处理，包括：将所述地震子数据和所述待处理地震数据分别转换为频率域的地震子数据和频率域的待处理地震数据；根据所述频率域的地震子数据，确定谱矩阵；确定所述谱矩阵对应的特征向量；基于所述特征向量，对所述频率域的待处理地震数据进行多次波压制处理，得到频率域的多次波压制后的待处理地震数据；将所述频率域的多次波压制后的待处理地震数据转换为时间域的多次波压制后的待处理地震数据。优选方案中，采用下述公式确定所述谱矩阵：其中，Gi(f)表示与第i个目标反射系数相关联的第i个目标反射源界面对应的谱矩阵，Xi(f)表示与第i个目标反射系数相关联的第i个目标反射源界面对应的频率域的地震子数据，表示Xi(f)的共轭矩阵。优选方案中，采用下述公式对所述频率域本文档来自技高网...

【技术保护点】
1.一种多次波压制方法，其特征在于，提供有目的工区的测井数据以及所述目的工区的时间域的叠后地震数据；所述方法包括：根据所述测井数据中的声波时差测井曲线数据和密度测井曲线数据，确定所述目的工区的钻井位置处的地层反射系数曲线数据；基于所述地层反射系数曲线数据，确定待处理地震数据，以及确定所述时间域的叠后地震数据中目标反射源界面对应的地震子数据；其中，所述待处理地震数据用于表征所述时间域的叠后地震数据中包含多次波的地震数据；所述目标反射源界面表示产生多次波的反射界面；根据所述地震子数据，对所述待处理地震数据进行多次波压制处理。

【技术特征摘要】
1.一种多次波压制方法，其特征在于，提供有目的工区的测井数据以及所述目的工区的时间域的叠后地震数据；所述方法包括：根据所述测井数据中的声波时差测井曲线数据和密度测井曲线数据，确定所述目的工区的钻井位置处的地层反射系数曲线数据；基于所述地层反射系数曲线数据，确定待处理地震数据，以及确定所述时间域的叠后地震数据中目标反射源界面对应的地震子数据；其中，所述待处理地震数据用于表征所述时间域的叠后地震数据中包含多次波的地震数据；所述目标反射源界面表示产生多次波的反射界面；根据所述地震子数据，对所述待处理地震数据进行多次波压制处理。2.根据权利要求1所述的方法，其特征在于，所述根据所述声波时差测井曲线数据和所述密度测井曲线数据，确定所述目的工区的钻井位置处的地层反射系数曲线数据，包括：根据所述声波时差测井曲线数据和所述密度测井曲线数据，确定所述目的工区的钻井位置处的波阻抗曲线数据；根据所述波阻抗曲线数据，确定所述目的工区的钻井位置处的地层反射系数曲线数据。3.根据权利要求2所述的方法，其特征在于，采用下述公式确定所述波阻抗曲线数据：其中，PI表示所述波阻抗曲线数据中指定连续排列的N个采样点对应的波阻抗平均值，其中，所述波阻抗平均值用于分配给所述N个采样点中的任意一个采样点；N表示指定采样点数量，ρi表示所述密度测井曲线数据中指定连续排列的N个采样点中第i个采样点的密度，dti表示所述声波时差测井曲线数据中指定连续排列的N个采样点中第i个采样点的声波时差。4.根据权利要求3所述的方法，其特征在于，采用下述公式确定所述地层反射系数曲线数据：其中，rj表示所述地层反射系数曲线数据中第j个采样点的地层反射系数，PIj和PIj-1分别表示所述波阻抗曲线数据中第j个采样点和第j-1个采样点的波阻抗。5.根据权利要求1所述的方法，其特征在于，所述基于所述地层反射系数曲线数据，确定待处理地震数据，包括：根据指定地震子波数据和所述地层反射系数曲线数据，确定所述目的工区的钻井位置处的深度域的测井合成地震道数据；从所述时间域的叠后地震数据中获取所述目的工区的钻井位置处对应的时间域的井旁地震道数据；基于所述深度域的测井合成地震道数据和所述时间域的井旁地震道数据，确定所述时间域的叠后地震数据中的目标多次波压制时窗范围，并将所述时间域的叠后地震数据中目标多次波压制时窗范围内的地震数据作为所述待处理地震数据。6.根据权利要求5所述的方法，其特征在于，所述基于所述深度域的测井合成地震道数据和所述时间域的井旁地震道数据，确定所述时间域的叠后地震数据中的目标多次波压制时窗范围，包括：将所述深度域的测井合成地震道数据转换为时间域的测井合成地震道数据；在所述时间域的测井合成地震道数据和所述时间域的井旁地震道数据的时窗范围内，移动预设窗长的时窗，并计算所述时间域的测井合成地震道数据和所述时间域的井旁地震道数据分别在当前时窗位置处的测井合成地震道子数据和井旁地震道子数据之间的相关性系数；确定小于预设相关性系数阈值的相关性系数对应的目标时窗位置，并基于所述目标时窗位置构成所述目标多次波压制时窗范围。7.根据权利要求1所述的方法，其特征在于，所述根据所述地层反射系数曲线数据，确定所述时间域的叠后地震数据中目标反射源界面对应的地震子数据，包括：确定所述地层反射系数曲线数据中地层反射系数对应的采样时间；对所述地层反射系数曲线数据中的地层反射系数按照由大到小的顺序进行排序，并将排序后的地层反射系数中前指定数量的地层反射系数作为所述指定数量的目标地层反射系数，并将所述目标地层反射系数对应的采样时间作为所述目标反射源界面在所述钻井位置处的界面时间点；其中，所述地层反射系数与所述目标反射源界面一一对应；将所述指定数量的目标反射源界面在所述钻井位置处的界面时间点分别标定于所述叠后地震数据上，并对标定后的叠后地震数据进行层位追踪和解释，得到所述目标反射源界面的时间层位；将所述时间域的叠后地震数据中与所述目标反射源界面的时间层位相邻的指定时窗范围内的地震数据作为所述目标反射源界面对应的地震子数据。8.根据权利要...

【专利技术属性】
技术研发人员：戴晓峰，徐右平，孙夕平，张明，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：北京,11