白云凹陷陆架—陆坡区多次波压制技术

白云凹陷位于中国南海北部陆架—陆坡区,由于水深变化剧烈,目标区地震资料中发育的表面多次波周期变化大,较难压制。传统的预测反褶积、SRME、Radon技术及其组合,压制这种水深变化剧烈的多次波效果不佳。为此,对近年发展的一种浅水表面多次波压制技术即确定性水层多次波压制技术DWD进行了研究,提出将其与SRME和Radon传统方法联合的多次波压制思路。DWD方法是先从水速动校后的最小炮检距的自相关剖面上拾取水底双程反射时间,然后在近道外推后将t-x域地震数据变换到τ-p域,再在τ-p域通过波场外推预测出多次波模型,最后与原始数据进行自适应匹配相减压制多次波。处理结果表明,DWD+SRME+Radon方法组合能有效地压制白云凹陷目标区多次波。与原处理剖面相比,此次处理后的成果剖面品质明显提高,降低了解释风险;同时,用此次处理得到的偏移道集进行叠前波阻抗反演,结果较好地刻画了目标砂体的特征和边界。     

各向异性高分辨率Radon变换压制多次波

针对VTI介质和海洋大炮检距地震数据反射波同相轴不能严格满足双曲时距关系,一次波和多次波能量在Radon域中聚焦性较差,多次波的压制效果不理想的问题,通过在双曲Radon变换积分路径中引入各向异性参数η,描述VTI介质和海洋大炮检距地震数据的时距曲线关系,在Radon域中有效聚焦一次波和多次波能量。同时,引入基于共轭导向梯度(CGG)算法的各向异性高分辨率Radon变换,基于L₁范数对模型空间进行稀疏约束,提高Radon变换的分辨率。模型数据和海上大炮检距实际资料测试表明,该方法能够很好地压制多次波。     

基于三维抛物线Radon变换的多次波压制算法对比分析

抛物线Radon变换是地震数据多次波压制的常用算法。三维Radon变换考虑了地震波场的三维传播特性,可更精确地描述波场的传播特征,能取得更好的效果。但受离散采样和有限采集孔径的影响,传统最小二乘三维Radon变换存在着分辨率低的问题,影响一次波和多次波的分离。为此,对多次波压制中的几种高分辨率抛物线Radon变换算法进行了对比分析研究。首先详细介绍了三维Radon变换的最小二乘算法,之后分别给出了迭代重加权、低频约束以及迭代阈值收缩三种高分辨率Radon变换算法,并用模拟数据测试了每种算法的计算效果、计算效率,分析了各算法的抗噪性能及对小剩余时差多次波的压制能力。结果表明,频率域高分辨率算法只能提高变换域横向分辨率,时间频率混合域算法可以同时提高横向和纵向分辨率;各高分辨率算法均具有较好的抗噪能力,且低频约束高分辨率算法在小剩余时差多次波压制方面效果最佳。     

三维高精度保幅Radon变换多次波压制方法

三维地震勘探现已成为地震勘探的主流方式。常规二维Radon变换多次波压制方法只针对二维地震数据,并未考虑地震波场三维传播的特点,即不适用于三维地震数据处理,故亟待探寻针对三维地震数据的处理方法。文中基于对三维Radon变换多次波压制方法的深入研究,针对三维地震数据变换域分辨率低的问题,采用迭代阈值收缩的方法提高变换的分辨率;针对数据中振幅随炮检距变化的特点,引入正交多项式变换,对沿不同曲率方向地震数据振幅的变化进行拟合。模拟数据和实际数据的测试结果表明,通过三维高精度保幅Radon变换可获得高分辨率的模型域数据,能有效分离一次波与多次波,同时多项式拟合可保护有效波的振幅,高保真地实现多次波的压制。     

双曲Radon变换法多次波衰减

多次波是相干干扰,在地球物理勘探中普遍存在。作部分动校正后的共中心点道集,多次波剩余时差更接近于双曲线。研究给出时不变双曲Radon变换空间频率域衰减多次波的基本原理,同时给出了计算公式和离散采样算法。时不变双曲Radon变换在频率域是解耦的,所以该方法在计算效率上和抛物Radon变换一样,但是在分离多次波上比抛物Radon变换法更精确。水平层状介质的多次波衰减结果表明,时不变双曲Radon变换法在叠前能有效地衰减多次波,复杂SMAART模型的叠前多次波衰减实例显示了该算法的有效性和实用性。     

迭代抛物Radon变换法分离一次波与多次波

Radon变换法是进行一次波与多次波分离的常用手段,最小平方约束下的频率域抛物Radon变换将t-x域数据转换到Radon域后,因存在剪刀状发散的截断效应,用传统方法难以彻底分离一次波和多次波。针对这一缺陷,提出了迭代抛物Radon变换法,即在Radon域截取一次波聚焦点附近很小区域内的数据为初始数据,经过Radon反变换和正变换后得到新的Radon域数据,然后用初始数据覆盖对应的小区域,经过迭代,最终得到保幅效果较好的一次波,且几乎不含多次波。利用相同的方法,也可以得到几乎不含一次波的多次波。最后通过理论模型和实际资料的处理,验证了本文方法的正确性和有效性。     

分级多域迭代组合压制多次波技术

在苏丹3/7区块的地震资料处理压制多次波过程中,过去存在压制多次波不足或者损伤有效信号的问题,这些问题可用分级多域迭代组合去噪技术解决:①根据多次波的周期性先用预测反褶积压制近偏移距的周期性多次波;②用Radon变换分别在炮域、检波点域和CMP域压制远偏移距的多次波;③在CRP道集上用高精度Radon变换压制剩余的多次波。通过应用上述技术,很好地压制了研究区块较发育的多次波,从而使处理成果品质有了较大改善,获得了多次波压制好、信噪比高的地震资料,为总结该区块潜山地区油气成藏规律,识别油气聚集带提供了可靠资料。     

加权抛物Radon变换叠加速度分析

抛物Radon变换（PRT）是进行叠加速度分析较为有效的方法,当缺失原始地震记录道时,加权抛物Radon变换可改善正变换模型空间域的聚焦性。基于水平或中等倾角地层CMP道集时距曲线为双曲线的假设,给出了抛物Radon变换叠加速度求取的基本方法,同时给出加权PRT求解计算时权系数的选取原则。在数据域加权后,求解矩阵的Toeplitz结构并没有破坏,所以仍可用Levinson递推法来进行计算。加权抛物Radon变换叠加速度分析算法容易实现,且计算精度较高,能较有效地处理含多次波的地震剖面。     

致密气藏压裂水平井产能计算方法

致密气藏开发普遍采用多段压裂水平井的开发模式。为了准确评价致密气藏压裂水平井产能并确定气井的合理配产,实现气井高效开发,基于保角变换理论和气水两相渗流理论,同时将基质有效渗透率作为变量来考虑压裂施工和气井产水对储层有效渗透率的影响,建立了压裂水平井气液两相产能方程。通过实际生产数据验证,结果表明：无因次泄气边界大于0.55时,气井生产压差随配产增加呈下凹型快速增长;相同气井产能条件下,水气比越大气井所需生产压差越大;水平段方向与K_y方向平行时,渗透率各向异性程度K_y/K_x越大,相同产气量时的生产压差越小;水平段与渗透率主值方向的夹角θ < 30°时,相同产气量条件下的气井生产压差几乎不变。因此,从降低压裂水平井储层压力损失的角度来考虑,布井时必须充分考虑渗透率各向异性程度和水平井水平段方向的影响,同时注意控制气井配产和采取必要的控水排水措施,以便达到更好的开发效果。     

基于Curvelet变换的多次波去除技术

本文基于波动方程的自由表面多次波压制预测减去法,采用近期迅速发展的多尺度变换系统中的Curvelet变换替代减去法,收到了较好的效果。该变换具有最优稀疏约束条件,能使一次波在一组基函数上的投影能量尽可能小。其主要实现过程为:对地震记录数据进行Curvelet变换,以预测多次波的Curvelet域的系数为阈值,采用类似去噪的手段去除多次波。该方法能满足去除多次波过程的能量最小准则。文中在Curvelet算法数值实验中尝试了先把地震记录做Radon变换,并用Curvelet阈值法减去Radon域的多次波,再把去除了多次波的数据做反Radon变换即为多次波去除结果。数值实验结果表明,在Radon域进行Curvelet变换去除多次波的效果更好,并能更好地保留一次波的能量。     

基于共散射点道集的多次波压制方法研究

多次波压制是地震数据处理中的重要环节。一种常用的多次波压制方法是根据一次波和多次波在成像空间的曲率差异并借助Radon变换进行分离和压制。成像空间既可以是叠前时间域也可以是叠前深度域。考虑到对复杂介质的适应性以及叠前处理对计算效率的要求,提出基于共散射点(CSP)道集的高分辨率Radon变换多次波压制方法。该方法只需要一个简单的初始速度场,就可将常规共中心点道集(CMP)映射到覆盖次数更高、炮检距覆盖范围更广的共散射点道集,然后在该道集上应用高分辨率双曲Radon变换,可较好地分离一次波和多次波。相对于以水平层状介质为假设条件的CMP道集,CSP道集更适应复杂地质构造,且时距关系满足双曲规律。模型和实际资料测试结果表明,该方法可以较好地实现速度谱能量团的聚焦,有利于较复杂地质条件下的多次波压制。     

渤海海域水平拖缆数据宽频处理关键技术

近年来海上浅水区宽频地震勘探技术进展缓慢,尤其像渤海这样平均水深不超过30m的浅水海域,宽频采集难以实施,只能通过常规采集、宽频处理才有可能实现这一目的,其中鬼波及浅水多次波是制约浅水资料宽频处理的两个主要因素。针对浅水海域水平拖缆鬼波时空变剧烈难以压制的问题,研究了τ-p域自适应参数估计的水平缆鬼波压制技术;针对渤海海域海底模型难以建立、浅水多次波压制困难的问题,研究了τ-p域自适应海底模型估计的浅水多次波压制技术。渤海M工区实际数据的应用结果表明,这两项关键技术能有效拓宽地震资料的频带宽度,实现渤海海域常规水平拖缆地震资料的宽频处理。     

最优化相似加权Radon变换压制多次波

多次波压制一直是地震资料处理的难题,常规方法多种多样,但都受到使用条件的限制,在实际应用中不能达到良好的效果。Radon变换是一种压制多次波的较好方法,如何提高Radon变换的分辨率一直困扰着地球物理工作者。在此提出了最优化相似系数加权法,该方法结合了相似系数加权法的优点能够较好地控制假频和截断效应,再结合Gauss—Seidel迭代的“最早迭代能量占优”的特点,能够提高Radon域的分辨率。通过合成理论记录和实际资料计算分析对比,该方法能较好地压制多次波,具有实际应用价值。     

地空频率倾子测深法及一维正演特征分析

为了充分利用地空频率域电磁法的三分量磁场信息,参照大地电磁法倾子测量技术,提出了一种地空频率倾子测深方法。首先,从水平电偶极子在一维层状介质中的全空间磁场公式出发,推导了地空频率倾子两个分量T_x、T_y的表达式,并通过快速汉克尔变换实现地空频率域倾子响应的一维正演;然后,对典型层状模型进行一维正演模拟,分析地空频率域倾子响应特征。通过分析均匀半空间模型和一维层状模型的正演结果,得到如下结论：①倾子响应整体沿场源方向及其垂直方向对称,同时存在场源畸变和过渡带畸变效应;②随收发距增大,倾子响应形态基本保持不变,幅值快速减小;③飞行高度变化时,倾子虚部幅值保持不变,实部幅值快速减小;④通过测量不同频率的倾子响应能够有效研究地下介质电阻率的垂向变化;⑤地空频率倾子测深法应在远区观测,可尽量避开畸变效应区,并合理选择T_x和T_y分量进行解释。文中模拟分析证明了该方法的理论可行性和有效性,为地空电磁法提出一种可能的发展方向。     

浅水多次波的联合衰减技术在渤海海域LD地区的应用

针对渤海海域LD地区浅水多次波发育的特点,提出了一整套完整的浅水多次波联合衰减技术。首先使用浅水去多次波(SWD)技术来去除大部分与海底相关的短周期多次波;对于剩余的与水面相关的长周期多次波采用了二维表面多次波压制(SRME)技术进行了有效去除;三维面元均化后进行Radon域去多次波处理;最后在叠前时间偏移后的共成像道集内用高精度Radon变换对多次波进行衰减。该方法通过对SWD技术和传统的去多次波技术合理的结合,达到了衰减多次波的目的。对渤海地区实际海洋地震资料处理表明:所提出的方法能够有效地衰减浅水多次波,浅水多次波联合衰减处理流程具有快速、简介的特点,同时处理后的地震资料保幅性较好。     

琼东南深海盆地三维地震多次波组合压制方法

在琼东南深海盆地,浅层为海相沉积特征,中、深层为陆相沉积特征,原始地震资料中多次波不仅异常发育而且复杂.压制多次波的方法有很多(褶积类、速度滤波类、波动方程类),但都具有一定的局限性,不能适应复杂情况下的多次波压制.为此,对各种压制多次波的方法进行了分析,针对琼东南深海盆地三维地震资料多次波的特点,采用调频滤波、抛物线Radon变换、加权中值滤波等方法,在不同的处理阶段分别对海底多次波、层间多次波和畸变多次波进行了压制,获得了好的效果.应用这种组合方式,不仅多次波被完全压制,一次有效反射得到了清晰的显示,而且没有损害有效信号.     

Radon变换压制层间多次波技术在高石梯-磨溪地区的应用

陆上地震资料中,层间多次波与一次波的速度差异较小,一次波的速度拾取困难,利用Radon变换压制多次波往往难以取得好的应用效果。在四川盆地高石梯-磨溪地区的地震资料中,震旦系-寒武系存在多次波的干扰,致使有效反射的振幅畸变,弱的有效反射被掩盖,造成测井合成地震记录与实际地震资料匹配差。该多次波具有干扰能量强、与有效波速度差异小的特点。基于处理和解释一体化的研究思路,在速度谱上利用时间域的构造层位来约束一次波速度的拾取过程,以减小多次波干扰对速度拾取的影响,提高平面上和纵向上速度拾取的精度;利用不同权系数的叠加速度进行扫描处理和Radon变换压制多次波,确定最佳匹配权系数,实现更高精度的多次波压制。在高石梯-磨溪地区,三维地震资料经压制多次波处理后,地震剖面上多次波的能量得到有效衰减,波组特征更为清楚,测井合成地震记录和井旁地震道相关性大幅提高。应用效果表明,这一技术在压制多次波的同时,能够有效地保护一次波能量,减少压制多次波的不确定性,提高地震资料的信噪比和成像效果,是压制陆上层间多次波的有效技术。     

逆散射级数和抛物线Radon变换联合的层间多次波压制策略

层间多次波给地震成像和解释带来极大挑战,单独使用基于波动方程预测或基于滤波的方法压制层间多次波效果都不理想,前者对数据的品质要求高且运算量大,而后者则很难处理规律性较差的层间多次波。基于此,提出了一种将逆散射级数法与抛物线Radon变换相结合的层间多次波压制新策略。首先在抛物线Radon域采用自适应非线性滤波对CRP数据中的一次波进行保护处理,将未受多次波干扰的一次波和含有部分一次波的多次波能量进行分离,然后在曲波域对滤波分离出的多次波能量和逆散射级数预测的层间多次波模型数据进行匹配相减,最后将相减得到的剩余一次波和滤波分离出的一次波相融合得到层间多次波压制后的结果。实际资料测试表明该策略能高效高精度压制层间多次波,使地层结构特征更为清楚,地震成像质量显著改善,从而为油藏解释提供更可靠的地震数据。     

电场作用下W/O乳化液中水颗粒的形变模拟研究与动力学分析

采用COMSOL Multipysics^TM软件,对W/O乳化液中分散相水颗粒在电场作用下的变形特性进行数值模拟,并结合已开展的室内实验,讨论电场强度(E)、电压波形、电场频率(f)、水颗粒半径(r)和连续相黏度(μ_c)对单个水颗粒变形的影响。结果表明:增加电场强度、增大水颗粒半径、降低连续相黏度会加剧水颗粒变形,电场强度存在上临界值;当电场频率较低时其值变化对水颗粒变形影响不明显,频率较高时则影响显著,但总体来看存在一个最优电场频率(2000 Hz);交流矩形波对水颗粒变形的影响强于交流锯齿波、直流脉冲波、交流三角波、交流正弦波等波形。当E为3~4 kV/cm、r为0.25~0.5 mm、使用交流矩形波时,水颗粒的拉伸发生增速变形,从而缩短水滴间的碰撞时间和W/O乳化液的电场破乳时间。     

应用加权迭代软阈值算法的高分辨率Radon变换

Radon变换的分辨率是其进行地震数据处理的关键因素。基于Bayes反演理论的迭代加权方法改善了Radon变换的分辨率,但其收敛速度较慢;由于Radon变换空间的强相关性,常用的迭代软阈值方法（ISTA）应用于Radon变换反演时的收敛速度也较慢,且分辨率较低。为此,将迭代加权最小二乘法嵌入ISTA中,形成了加权ISTA算法。该方法引入高分辨Radon变换中加权矩阵的思想,利用Radon参数的先验信息约束反演误差函数,克服了ISTA收敛速度慢、分辨率低的缺点。合成记录和实际地震资料处理结果表明,该方法提高了Radon变换分辨率,在多次波分离、噪声压制中效果显著。