一种可控震源次谐波噪声的压制方法及系统.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN110888165A(43)申请公布日2020.03.17(21)申请号201811050486.4(22)申请日2018.09.10(71)申请人中国石油化工股份有限公司地址100728北京市朝阳区朝阳门北大街22号申请人中国石油化工股份有限公司石油物探技术研究院(72)发明人张洁(74)专利代理机构北京知舟专利事务所(普通合伙)11550代理人郭韫(51)Int.Cl.G01V1/36(2006.01)权利要求书3页说明书6页附图8页(54)发明名称一种可控震源次谐波噪声的压制方法及系统(57)摘要本发明提供了一种可控震源次谐波噪声的压制方法及系统，属于地震资料处理领域。该方法将含次谐波噪声的原始地震记录从空间-时间域变换到次谐波域，使得次谐波噪声与震源基波信号分离开，对次谐波噪声进行切除，再从次谐波域反变换到空间-时间域，得到压制次谐波噪声的相关前记录。利用本发明方法能够使得谐波信号更容易被检测到并得到压制，在有效保留原始地震道有效信号的情况下，消除了震源信号中的次谐波干扰噪声。CN110888165ACN110888165A权利要求书1/3页1.一种可控震源次谐波噪声的压制方法，其特征在于：所述方法将含次谐波噪声的原始地震记录从空间-时间域变换到次谐波域，使得次谐波噪声与震源基波信号分离开，对次谐波噪声进行切除，再从次谐波域反变换到空间-时间域，得到压制次谐波噪声的相关前记录。2.根据权利要求1所述的可控震源次谐波噪声的压制方法，其特征在于：所述方法包括：步骤S1：输入含次谐波噪声的原始地震记录sc(t)，并根据扫描信号参数生成次谐波信号qc(t)；步骤S2：将次谐波信号qc(t)与含次谐波噪声的原始地震记录sc(t)进行互相关处理，得到变换到次谐波域的相关记录；步骤S3：利用面波时差对步骤S2得到的变换到次谐波域的相关记录进行线性时差校正，经过线性时差校正后的次谐波噪声集中在一个水平时窗内；步骤S4：对步骤S3得到的集中在水平时窗内的次谐波噪声开取时窗，对时窗内的次谐波噪声进行置零，得到压制噪声后的地震记录；步骤S5：利用面波时差将步骤S4得到的压制噪声后的地震记录进行反线性时差校正得到反线性时差校正后的地震记录，然后利用次谐波信号qc(t)对所述反线性时差校正后的地震记录进行反相关运算，将其从次谐波域变换到空间-时间域，得到压制次谐波噪声的相关前记录。3.根据权利要求2所述的可控震源次谐波噪声的压制方法，其特征在于：所述步骤S1中输入的含次谐波噪声的原始地震记录sc(t)是实际采集到的地震记录，对其假设如下：sc(t)＝r(t)*(q(t)+qc(t))，其中，q(t)为震源基波信号，qc(t)为次谐波信号，r(t)为反射系数；所述步骤S1中的根据扫描信号参数生成次谐波信号qc(t)的操作包括：所述扫描信号参数包括：起始频率、终止频率和扫描长度；将所述起始频率的一半作为次谐波信号的起始频率，将所述终止频率的一半作为次谐波信号的终止频率；将所述扫描长度作为次谐波信号的扫描长度；根据所述次谐波信号的起始频率、次谐波信号的终止频率、扫描长度生成次谐波信号。4.根据权利要求3所述的可控震源次谐波噪声的压制方法，其特征在于：所述步骤S2的操作包括：利用下式进行互相关处理：sc(t)*qc(-t)＝(r(t)*q(t)+r(t)*qc(t))*qc(-t)＝r(t)*q(t)*qc(-t)+r(t)*qc(t)*qc(-t)＝r(t)*ccf{q(t)，qc(t)}+r(t)*acf{qc(t)}其中，qc(-t)为低次谐波，ccf{q(t)，qc(t)}为震源基波信号和次谐波信号的互相关，acf{qc(t)}为次谐波信号的自相关。5.根据权利要求2所述的可控震源次谐波噪声的压制方法，其特征在于：所述步骤S3中的利用面波时差对步骤S2得到的变换到次谐波域的相关记录进行线性时差校正的操作包括：2CN110888165A权利要求书2/3页根据不同接收道的偏移距和面波视速度，计算各道面波的到达时间，利用各道面波的到达时间计算得到面波时差，利用面波时差对步骤S2得到的变换到次谐波域的相关记录做线性时差校正。6.根据权利要求2所述的可控震源次谐波噪声的压制方法，其特征在于：所述步骤S4中的时窗的宽度为次谐波噪声的一个波长的宽度。7.一种实现权利要求1-6任一所述的可控震源次谐波噪声的压制方法的系统，其特征在于：所述系统包括：输入单元：用于输入含次谐波噪声的原始地震记录sc(t)，并根据扫描信号参数生成次谐波信号qc(t)；互相关单元：与输入单元连接，用于将次谐波信号qc(t)与含次谐波噪声的原始地震记录sc(t)进行互相关处理，得到变换到次谐波域的相关记录，并