基于径向道变换的svd自适应面波压制方法
【技术领域】
[0001] 本发明属于油气勘探地震资料处理技术领域，特别是一种以保护低频信息为基础 的基于径向道变换的SVD (SVD-奇异值分解）自适应面波压制方法。 现有技术
[0002] 面波是地震勘探中广泛存在的一种规则干扰波，具有低频、低速、高能量、横向衰 减较慢的特点。目前在资料处理过程中大多采用高通滤波的方法来压制它，但这同时也对 和面波处于同一频带的信号成分造成损害。一般消除面波的方法是利用面波和反射波的频 率差异及视速度差异，如一维滤波和F-K滤波。但因面波的频带往往和反射信号的频带有 重叠，即便是F-K域滤波也不能只对与有效信号相同频带的噪声进行压制，其实在滤波的 同时也对有效波的低频成分造成了损害，而且斜坡长（陡度小）的滤波响应在反褶积或谱 白化处理后,面波往往又被突出出来。
[0003] 近年来，不断有人尝试用新的方法来消除面波，如线性调频理论、tau-p变换法、 外科手术式的切除法等，这些都有一定程度的改进，其中以切除法效果较好。该方法是先 作FFT，自动分离出其低频分量，再作IFFT，建立内切除时窗，并用加权减去法剔除面波的 主要能量，其必要条件为原始单炮上面波的振幅必须比有效波强得多。但该方法在切除时 仍然会把属于该频带的一些有效信号一起切去。
[0004] 另外，也有学者提出在野外数据采集阶段采用不等间距接收方式，以便使室内处 理时可以不采用普通的二维滤波，而采用高保真云面波技术，这样就在消除面波的同时保 留信号的低频成分，达到使地震资料频谱向低频端展宽的目的，从而提高了记录的分辨 率；近几年，利用小波变换与特征值分析结合去面波也取得一定效果，该方法以较低频率处 面波的能量强于反射波、在小频率范围和小空间范围内面波能量变化缓慢为假设条件，利 用小波变换进行频带分解，用面波的视速度对面波作线性时移，再使面波逐道相干，最后 利用KL分解的方法来提取面波。但是对面波做人为时移时需要人为干涉，可操作性差，并 且面波发生频散，很难做线性时移，同时部分有效波相干性也得到提高，用KL分解时会造 成有效息损失。
[0005] 由于面波的频率特征，很多面波衰减方法在处理过程中极容易误伤有效信息，造 成与面波频带重合的有效低频成分的丢失。而随着油气勘探工作的不断深化，地震资料处 理中对地震低频信息重要性的认识越来越深入。低频信息在改善深部地层成像质量、速度 分析、叠前反演以及"低频伴影"直接油气检测等方面具有重要意义，在数据处理过程中如 何最大限度地保护好低频信息显得越来越重要。
[0006] 发明目的
[0007] 本发明目的就是针对上述问题，提供一种最大限度保护低频有效信息的基于径向 道变换的SVD自适应面波压制方法。
【发明内容】
[0008] 本发明为了解决上述问题而提出的技术方案包括：
[0009] (1)该方法所采用的技术路线
[0010] 基于径向道变换的SVD自适应面波压制方法，其实质就是先利用径向道变换使面 波频带向更低频率范围移动，然后利用低通滤波在径向道域得到低频面波信号（可以将低 通滤波频率设置得相对小点），再利用SVD将上一步得到的面波信号中的有效低频信号分 离出来，最后从原始数据中将面波减去。为了尽可能消除径向道变换插值带来的误差，可以 对整个数据先进行一次径向道变换，该方法可以根据面波的视速度合理选择作用范围，减 小对面波范围之外信号的影响。利用该方法可以在不损失有效低频信号的前提下减去面波 干扰，为后期叠前反演等特殊处理打下基础。
[0011] 基于上述技术路线的基于径向道变换的SVD自适应面波压制方法，包括：
[0012] 步骤101 :选择含有典型面波的原始单炮记录作为输入数据；
[0013] 步骤102 :对输入的含有典型面波的原始单炮记录针对面波区域作频谱分析，以 确定面波的频带范围；
[0014] 步骤103 :选择合理的速度扫描范围以及扫描步长对原始单炮记录作一次径向道 正变换；
[0015] 步骤104 :对步骤103得到的数据做一次径向道反变换，消除径向道变换过程中插 值造成的影响；
[0016] 步骤105 :设置低通滤波器，在径向道域进行低通滤波，得到含有少量低频有效信 号的面波噪声；利用最大面波视速度作为径向道变换的速度扫描约束，即径向道正变换后 径向道域只存在最大面波视速度范围内的信号，这样就只在面波存在的范围进行径向道变 换；
[0017] 步骤106 :对步骤105得到的数据作径向道反变换，还原到时间域信号；
[0018] 步骤107 :对步骤106得到的时间域信号作SVD分解，利用面波和有效信号相干性 不同将面波中的有效低频信号分离出来；
[0019] 步骤108 :利用奇异值拐点法确定奇异值范围，利用确定出的奇异值范围对分离 后的面波进行剩余干扰SVD重构，得到不含有效低频信号的面波成分；
[0020] 步骤109 :用步骤104得到的数据与步骤108分离出来的面波信号相减，得到面波 压制后的数据；
[0021] 步骤110 :输出最终去噪结果。
[0022] 上述方法中，步骤102确定的面波频带范围是I-IOHz。
[0023] 步骤104中所述径向道反变换主要是消除在步骤109相减过程中造成的误差影 响，因为步骤103中的径向道正变换和步骤106中的径向道反变换需要沿X方向进行两点 线性插值才能将地震道振幅从炮检距-双程旅行时坐标系变换到视速度-双程旅行时坐标 系，所以步骤103、106引入插值误差。
[0024] 当步骤104对整个单炮数据做一次径向道反变换后，再利用步骤109做数据相减 时就可以抵消插值带来的影响。
[0025] 径向道正变换、径向道反变换公式如式1 :
[0027]
[0028] 式中，R为径向道正变换标志，R 1为径向道反变换标志，S为原始数据，
[0029] S'为变换后数据,X为炮检距,t为双程旅行时,V为视速度,t'为平移后的双程旅 行时，（x〇，t。）为坐标原点。
[0030] 发明效果
[0031] 基于径向道变换的SVD自适应面波压制方法最大程度地压制了面波干扰，有着其 他技术不具备的优势，其具体优势和特点表现在以下几个方面：
[0032] 第一、技术效果的可靠性。该方法综合考虑面波与有效信号在频率及相干性两方 面的差异，在面波视速度范围内对面波进行压制，并且利用数据转换将面波频率变小，更利 于面波的预测，效果明显。
[0033] 第二、操作简单易实现。该方法流程及参数设置简单，运算速度快，不依据任何假 设，不受常规面波压制方法应用的限制。
[0034] 第三、低频信息的保护。利用SVD进行低频保护，在压制面波的同时最大限度地保 留了地震信号的有效低频成分，保持了有效信号的频宽，为改善深部地层成像质量、速度分 析、叠前反演以及"低频伴影"直接油气检测等提供了数据保障。
【附图说明】
[0035] 图1为基于径向道变换的SVD自适应面波压制方法流程图。
[0036] 图2为原始含面波单炮记录及频谱图，a)为含面波原始单炮记录，b)为对应单炮 记录的频谱图。
[0037] 图3为径向道域低频滤波压制面波单炮效果，a)为含面波原始单炮，b)为径向道 域低频滤波压制面波单炮，c)为压制的面波区。
[0038] 图4为径向道域低频滤波前后对比图，a)为低频滤波前，b)为低频滤波后。
[0039] 图5为基于径向道变换的SVD自适应面波压制方法压制面波单炮效果，a)为含面 波原始单炮，b)为径向道变换的SVD自适应面波压制方法压制面波单炮，c)为压制的面波 区。
[0040] 图6为区域滤波压制面波效果及压制后的频谱，a)为含面波原始单炮，b)为低截 频滤波后单炮，c)为区域滤波后的频谱。
[0041] 图7为不同方法面波压制后频谱对比，a)为基于径向道变换SVD面波压制方法， b)为径向道域低频滤波，c)为区域滤波法。
[0042] 图8为奇异值选取对比图，a)为去除第一个奇异值的结果，b)为去除前三个奇异 值的结果，c)为去除前五个奇异值的结果。
[0043] 图9为奇异值分布图。
【具体实施方式】
[0044] 本实施例以XX油田XX地区三维地震资料中的一个单炮记录，说明实际应用效果， 具体流程图见图1。该单炮记录时间长度3000ms，时间采样间隔为4ms，采样点数751，有71 道数据。采用上述