近地表勘探中炮检距对瑞利波频散成像效果的影响

在油田近地表勘探中,可以利用现有的地震大炮记录中的瑞利波信息进行频散成像,提取瑞利波频散曲线反演近地表的横波速度结构。由于大炮地震记录的采集排列很长,且不同炮检距范围内瑞利波数据的频散成像结果存在差异,从而会对频散曲线的提取和反演结果产生影响。笔者利用胜利油田LJ地区的盒子波地震数据,从频散成像的分辨率和连续性两个方面,分析了不同排列长度和不同偏移距下的频散成像结果,发现长排列可以有效提高频散成像的分辨率,且在较近偏移距下的频散成像效果更好。同时,验证了Zhang2004年提出的最优偏移距计算公式在石油近地表勘探中的有效性及适用性,并可以作为近地表瑞利波勘探时最优偏移距的选择标准。     

地质条件复杂地区瑞雷波勘探正演模拟

为研究复杂地质条件对瑞雷面波频散曲线的影响,采用有限差分地震数值模拟方法,对3种模型进行正演模拟,得出了某一时刻的波场快照与单炮记录。对模拟记录进行频散分析的结果表明:单层均匀介质没有频散现象;起伏地表两层均匀介质的频散曲线相速度只和介质本身的参数有关,不受地形的影响;水平地表含软夹层介质中,瑞雷波频散现象非常明显。      

基于频散曲线合成面波地震记录的方法

准确提取频散曲线是瑞雷波勘探的重要环节,检验各种频散曲线求取方法的正确性和稳定性至关重要。基于频散曲线,选择抽样脉冲信号作为子波,推导出了合成单炮面波地震记录的理论公式,并利用该公式,针对不同弹性层状介质模型的频散曲线合成了面波地震记录。通过对其波场特征对比和频谱分析,同时采用频率-波数域法反求其频散曲线,结果与模型频散曲线几乎相同,从而充分验证了该面波合成方法的正确性。      

Matlab环境下瑞利波有限差分正演与曲线绘制

以瑞利波频散方程为出发点的Abo-Zena传递矩阵等方法,只能研究层状介质中瑞利波的传播特性,对于非层状介质,传递矩阵无能为力。因此,不得不考虑有限元、有限差分等方法。根据弹性动力学方程,采用交错网格有限差分方法对均匀弹性半空间介质进行全波场模拟,并在Matlab环境下实现编程计算,再现了瑞利波在近地表的传播状态,实现了地震剖面曲线的正振幅充填问题,从而获得了类似野外地震勘探的剖面记录。模拟结果表明,采用一阶差分格式所得到的地震记录存在较严重的数值频散,在采用有限差分法研究瑞利波"之"字形成因时,须尽量压制数值频散现象,才能获得更接近实际情况的频散曲线。      

F-K域三层对称模型洛夫型导波频散分析

导波(分瑞利型和洛夫型2种)在波导层中传播时,其能量持续时间长,而且传播过程中携带着波导层的结构信息。基于多道分析的F-K方法(频率波数域二维傅立叶变换),分析研究了3层对称模型中的洛夫型导波的频散特性,计算了理论频散曲线,并据此合成了特定频率的导波记录。通过对合成记录的F-K域分析,获得了理论模型的速度参数。从分析结果来看,速度参数与理论合成参数相一致,说明F-K方法可作为导波的频散分析。      

Love型槽波透射记录的频域合成算法

利用射线追踪对网格化的二维模型进行走时计算,在频率域对Love型透射槽波记录进行了合成。在考虑了频散、吸收衰减等对波场振幅、相位影响的基础上,对含薄煤带模型的透射槽波理论记录进行了频域合成,并进行了波场特征分析。结果表明:频域合成的槽波透射记录在时间、频散特征、能量等属性上与理论分析一致,且与三维弹性波波动方程模拟的基阶Love型槽波波场特征相似,但频域合成算法的计算速度快1~2个数量级。     

基于NAD算法的声波方程时间四阶差分解法

波动方程数值模拟是研究地震波传播机理的重要工具,有限差分求解波动方程是当前地震波数值模拟的主要方法之一。当地下介质中的地震波速度较低或地震波高频成分丰富时,常规有限差分技术常常产生严重的数值频散误差,这种误差会降低数值模拟的精度,影响对地震波传播机理的分析。为压制地震波数值模拟时产生的数值频散误差,提高波场模拟精度,提出了基于NAD算子的时间四阶精度波动方程差分格式。根据对应的差分格式,分析了该差分格式的数值频散关系。与常规四阶精度差分算法的频散曲线相比,基于NAD时间四阶精度差分方法不但能够实现时间频散的有效压制,同时其基于更多网格点的位移分量和位移梯度分量空间微分求解方法还能够实现空间频散的有效压制。另外在相同模型条件下,基于NAD算法的声波方程时间四阶差分解法可采用大网格对模拟空间进行差分离散,减少网格数,提高计算效率。     

优化交错网格差分算子的井间地震数值模拟及其频散分析

在井间地震有限差分数值模拟中,用离散化的高阶差分方程近似连续导数的波动方程时,不可避免地会产生数值频散,而数值频散程度则直接影响到地震波数值模拟精度,因此为了得到清晰准确的地震波场记录,必须尽可能地压制数值频散。这里在一阶速度应力弹性波方程的基础上,利用两个约束条件构造拉格朗日函数获取优化差分系数,与泰勒展开差分系数下的交错网格高阶差分模拟结果比较,发现改进的优化交错网格差分算子的高阶差分数值模拟能更有效地压制数值频散,进一步提高交错网格高阶差分数值模拟的精度,为高精度井间地震数据的波场成像、纵横波联合解释等提供可靠依据。     

声波方程隐式高阶有限差分数值模拟

推导了声波方程空间二阶导数的隐式求解公式及差分系数的求解方法,讨论了该方法的数值频散特征。利用该方法分别对均匀介质及Marmousi模型进行了数值模拟,将其结果与传统的显式差分格式的模拟结果进行了对比分析。结果表明:该方法较传统的显式求解方法具有更低的数值频散、更高的计算精度。      

面向小尺度地质目标体的地震信号谐波预测拓频技术及应用

针对我国东部油田勘探开发中小断层、薄层、薄互层等小尺度地质体的探测需求,开展了地震数据拓频处理方法研究,并在胜利油田实际地震数据提高分辨率处理中进行了应用。在分析信号谐波理论的基础上,在实际地震数据提高分辨率处理中,提出了进行谐波与次谐波预测的地震数据拓频处理技术及流程。合成地震记录实例验证了利用基波预测谐波方法的可行性,在胜利油田面向小尺度地质目标体地震数据处理的应用中,验证了拓频处理技术的正确性,处理后的地震数据的频带宽度得到明显地拓展,达到了提高地震数据分辨率的目的。同时,实际数据经过拓频处理后,地震资料的信噪比也得到很好地保持,研究目标区的小断层的断点清晰。     

可控震源地震勘探的数值模拟

用可控震源激发的线性扫频（Chirp）信号模拟了二维地震勘探,采用相关处理技术把反射Chirp地震信号转化为脉冲地震信号形成了相关地震剖面,把相关地震剖面与脉冲震源产生的合成地震记录剖面进行了对比。相关地震剖面同相轴波形的最大绝对值出现的位置反映了界面的真实位置信息。相关地震剖面比脉冲震源产生的合成地震记录剖面具有更高的分辨率。对震源扫频信号进行适当地锥化处理可以部分地压制相关地震剖面中反射地震信号的旁瓣。     

分时窗多子波方法在制作高精度合成记录中地应用

人工合成地震记录,在地震资料的解释和储层预测中起着重要的作用.通过对地震子波时变和空变特点的剖析,阐述了划分时窗的必要性,进而结合合成地震记录的制作原理和方法,提出了综合利用雷克子波和井旁地震道提取子波,并进行分时窗制作高精度合成地震记录的一般方法.最后将该方法应用于实例,得到的合成记录不仅与实际地震剖面吻合较好,而且...     

高精度合成地震记录制作方法及应用

合成地震记录的精度直接影响到地震地质层位的准确标定。在总结前人工作的基础上，提出一套高精度合成地震记录的制作方法：（1）利用VSP资料能准确地反映地震波在地震中传播速度的变化规律的优点，进行漂移校正，以提高深－时转换的精度；（2）利用密度测井资料与声波时差曲线联合计算反射系数序列，以提高反射系数序列的精度；（3）利用分时窗提取子波技术，解决长期以来合成地震记录中存在的子波问题。将此方法用于合成地震记录的制作，在提高合成地震记录与地震剖面的吻合度和分辨率方面，都取得了较好的效果。     

小波变换与地震资料谱白化

从地震记录的褶积模型出发，讨论了地震记录小波变换的物理意义。经深入分析发现，恰当地选择小波函数及尺度因子α，对地震记录作了小波变换，相当于对地震记录，作谱白化，其特点是运算速度快，参数易于选择，可使不同目的层的信噪比与分辨率达到最佳折衷，模型及实际算法结果，证明了该方法的有效性。     

子波整形提高合成地震记录质量

阐述合成地震记录的制作原理,利用测井数据计算出反射系数。根据地震波传播机制,通过三维地震数据统计得到子波,由此计算出原始合成记录。再根据过井剖面建立褶积模型,对子波进行整形,得到一个随时间变化的地震子波,最终得到与地震数据的波组特征相吻合的合成地震记录。为精细储层描述打下了坚实的基础。     

Random seismic noise attenuation data using the discrete and the continuous wavelet transforms

In this paper, we propose a technique of random noise attenuation from seismic data using discrete and continuous wavelet transforms. Firstly, the discrete wavelet transform (DWT) is applied to denoise seismic data using the threshold method. After, we calculate the continuous wavelet transform of the denoised seismic seismogram, the final denoised seismic seismogram is the continuous wavelet transform coefficients at the lower scale. Application to a synthetic seismogram shows the robustness of the proposed tool for random noise attenuation. Application to real vertical seismic profile recorded in Algeria clearly shows the efficiency of the proposed tool for random noise attenuation.     

粒子群退火算法在转换波剩余静校正中的应用

求取转换波剩余静校量是非线性问题。为准确求取转换波剩余静校正量,将两种非线性全局最优化方法: 粒子群算法和模拟退火算法融合成为粒子群退火算法( PSOSA) 。为验证方法的有效性,建立了单层模型,合成带有剩余静校正量的转换波地震记录,对该地震记录利用PSOSA 求取剩余静校正量并进行静校正。结果表明: 静校正后的叠加剖面和单炮记录同相轴连续性较好,求得的剩余静校正量误差较小,证明了PSOSA 能较准确地求取转换波的剩余静校正量。     

A wavelet-based multifractal analysis of seismic data for facies identification. Application on the pilot KTB borehole

The main goal of this paper is to identify reflection from seismic data using the multifractal formalism. Firstly, the wavelet transform modulus maxima lines method is applied with a sliding window of 128 samples at the seismic seismogram data. After that, we estimate the generalized fractal dimension. Application at the noisy synthetic seismic seismogram of the pilot Kontinentales Tiefbohrprogramm de Bundesreplik Deutschland borehole shows that the three fractal dimensions related to the three first moments are a good tool that can separate between reflection dues to the change of facies and random noise.     

宽频地震处理技术在南海某工区的应用

在海洋地震勘探中产生的鬼波会抑制低频分量,造成频域的陷波点,降低地震记录的分辨率,使地震记录的同相轴发生混叠甚至产生假的同相轴。笔者采用FK域水平电缆压制鬼波法,将地震数据进行FK变换,对维纳滤波求取反褶积算子投影到FK域,并对二者进行褶积,将结果 FK反变换回时间域。鬼波压制后的模型数据频带得到了有效拓宽,低频、高频信息均得到了加强。鬼波压制后的实际数据目的层频带从12~60 Hz拓宽到8~75 Hz,三亚组能量与合成记录匹配较原始数据更好。     

溱潼凹陷边城-叶甸地区构造解释及储层预测

通过使用Landmark软件包的合成地震记录、相干技术对溱潼凹陷边城—叶甸地区地震资料进行精细构造解释，应用Gridstatpro油藏描述技术对储层进行预测。分析了该区虽多口钻井油气显示活跃，但并未发现工业油流的原因，在此基础上，进一步落实了局部构造，优选出可供上钻的有利圈闭3个。