地震子波处理的二步法反褶积方法研究

针对玛湖斜坡区三块三维地震资料和赛汉塔拉凹陷二块三维地震资料连片处理中的特点,结合地质任务和处理目标要求,提出了地震数据连片处理中的地震子波处理的方法.该方法主要体现了两次反褶积,一次是采用地表一致性反褶积,将不同震源的频带拓宽到一个标准上;再一次采用相位校正反褶积,将不同震源的数据校正到相同相位上.为了保证提取的相位校正反褶积算子稳定,采用叠后地震道提取（主要考虑到叠后地震道信噪比高,算子稳定性强）,然后将该算子应用到叠前地震道,进行相位校正.     

井间地震逆高斯束共反射点叠加成像（英文）

为了解决地下复杂微小构造的精细勘探问题,使用基于高斯射线束理论的逆高斯束叠加成像方法对井间地震反射波进行成像研究。井间地震常规射线类叠加成像方法由于覆盖次数受到限制,构造复杂区精细成像质量欠佳,而以波动方程为基础的井间地震偏移方法的成像时效性不高。本文方法借鉴高斯束合成地震记录的思想,将共炮集地震数据逆高斯束分解成共反射点道集数据,然后选取合适的面元进行共反射面元数据叠加,实现了井间地震逆高斯射线束共反射点叠加。与传统的VSP-CDP叠加成像方法相比,成像范围更加广泛,且适应复杂地质构造。该方法不仅能够对二维井间地震勘探资料进行成像,针对三维井间地震资料采用基于宽线处理思路的逆高斯束叠加成像方法仍可以处理复杂构造及斜井成像问题。理论模型及实际资料试算验证了本文研究方法的有效性与稳健性。     

井间地震技术在松南油气开发中的初步应用研究

为了研究松南地区陆相砂泥岩薄互层储层岩石物性横向变化、微构造等对油气聚集的影响和探索利用高分辨率井间地震技术直接指示油气分布的可能性,开展了井间地震及地面小三维地震、VSP、测井等立体地震观测．使用了自行研制的井中炸药震源和常规地震记录仪器．采用层析成像方法、反射资料叠前偏移成像方法等对采集到的资料进行研究,得到了不同分辨率的地下构造和岩石物性信息,展示出不同方法技术在分辨薄储层能力上的差异．高分辨率的井间地震成果较好地实现了井间地层对比,查清了两井间储层横向变化和油气分布情况．层析速度与构造信息等的综合应用解释了两井油气产能差异的原因,建立了油气地质模型．所得成果显示出井间地震技术在油气开发中的潜力．     

井间地震资料精细解释方法研究与应用

井间地震资料分辨率很高,但是缺乏成熟的解释方法,本文针对井间地震资料的特点,借鉴地面地震解释技术对井间地震资料开展了精细解释方法研究.首先在深度域资料特点分析的基础上,利用区域速度和层析成像速度将深度域井间地震资料转换到时间域进行研究,并通过粗化降频建立起井间地震资料与地面地震资料之间的相似性联系,得到对井间地震低频反射特征的认识,结合井旁道合成地震记录标定建立起井间地震资料的反射特征,准确的解释了井间地震的反射层位.采用地震瞬时属性分析和时频属性分析得到了对沉积环境的认识,进一步通过波阻抗反演更好的认识井间储层特征.在实际应用中综合多种方法进行解释,并结合测井资料对比落实了井间砂体展布特征,得到对井间储层的精细描述结果,解决了两口井之间储层分布范围、横向连通性及厚度变化等问题,取得了较好的应用效果.     

大庆油田M1井3D VSP成像方法与应用

3D VSP地震数据具有高分辨率和高信噪比的特点,常规的VSP成像方法在大庆油田资料的处理中一直未见到理想的效果.依据大庆油田VSP资料的特点,设计了3D VSP资料的处理流程,主要包括三分量波场逐点定向方法和Q补偿的子波反褶积等预处理方法,所设计的流程具有高分辨率、高保真性的特点.大庆油田M1井的实际资料处理结果显示,所形成的方法可将VSP资料的频带展宽约20 Hz,同时对地面地震也进行了反褶积处理,反褶积后,地面地震的分辨率得到有效提高,波组特征更加清楚.实际数据计算表明,所形成的方法具有很好的实用性,可实现地下介质的精细刻画,可在具有同类地质条件的研究区推广应用.     

三维斜井XSP-CDP叠加成像

为了解决三维斜井井间地震数据成像问题,借鉴二维地面地震弯线及宽线处理思路,将三维井间地震成像问题转换成二维成像问题进行研究.利用三维逐段迭代方法对三维斜井井间地质模型进行正演模拟,计算上、下行反射波走时及反射点位置,通过对比正演记录与实际资料的特征,调整地质模型参数获取准确速度场.然后对波场分离后的共炮集数据进行反射波归位,得到共反射点道集数据.最后根据反射点的分布选择一条反射点分布最集中的剖面线,以一定的宽度将剖面线左右的反射点包裹起来;根据叠加宽度划分成大小相等的矩形面元,将面元内的共反射点数据投影到剖面线上进行二维叠加成像.通过模型测试,三维XSP-CDP叠加方法能有效反演地下结构特征,指导三维井间地震采集观测系统的设计并解决不共面斜井叠加成像问题.应用该方法对M地区实际资料进行成像研究,结果表明该方法可有效反演地下三维地质结构特征,为进一步数据处理提供依据.     

井间地震应用效果分析

介绍了井间地震方法在胜利油田的应用情况，列举了永安镇和罗家两种不同地质情况的成功实侧，涉及到井间地震的资料采集、处理及解释等方面。分析了井间地震与地面地震、VSP等方法相比的优劣所在，指出了井间地震的应用前景。     

多尺度地震资料联合反演方法研究

常规三维地面地震反演不可避免的存在多解性和分辨率不高的缺陷,而油藏地球物理阶段丰富的多尺度地震资料为减小多解性、提高分辨率提供了可能.基于贝叶斯反演理论,通过联合概率分布建立新的似然函数,将三维地面地震、VSP和井间地震三种多尺度资料有机地融合在一起,完善了多尺度地震资料联合反演框架及反演流程.模型测试及实际资料处理表明,联合反演算法有效地引入了小尺度地震资料中的高频信息对大尺度资料进行约束,反演结果在保留大尺度地震资料特征的基础上提高了分辨率,降低了多解性,同时促进了多种地震资料之间的相互匹配.     

空间反射结构正则化多道反褶积（英文）

在高分辨率地震资料处理中通常采用稀疏脉冲反褶积方法来展宽地震数据的有效频带,增强地震数据刻画薄层结构的能力。但是,现有的稀疏脉冲反褶积方法采用逐道运行模式,未考虑地震信号在相邻道之间的空间约束关系。受高频噪声的影响,反褶积结果不仅具有较强的多解性,且降低了地震信号的空间连续性。为此,本文提出了一种空间反射结构正则化多道稀疏脉冲反褶积方法。该方法从地震数据本身提取空间反射结构表征算子,将表征算子引入到多道稀疏脉冲反演的正则化条件,增强了稀疏脉冲反褶积方法的抗噪性和稳定性。理论模型实验和实际资料处理证实了该方法的正确性和可行性。     

测井、井间地震与地面地震数据联合约束下的地质统计学随机建模方法研究

利用能够整合测井信息与井间地震信息的地质统计学随机模拟方法,结合传统的地质统计学反演思路,得到了一种能够同时整合测井、井间地震与地面地震三种先验信息的地质统计学反演与储层建模方法.由于井间射线信息、测井信息与地面地震数据在随机反演与建模过程当中都得到了尊重,因此与传统地质统计学反演仅利用了测井与地面地震数据相比,本文的地质统计学反演与建模方法更充分地利用了先验信息,有效提高了反演的精度,降低了随机建模中的多解性.基于理论数据的测试证实了上述观点.     

地震盲反褶积综述

从动态系统角度对地震勘探中盲反褶的有关问题进行分析,将线性与非线形地震盲反褶积纳入盲系统辨识的统一框架．综述了几种重要地震盲反褶积方法的原理与效果．指明起源于瞬时混合盲反褶积(或盲源分离)的独立分量分析(ICA)方法在解决一般化地震盲反褶积问题上的重要应用前景,并给出两种可能的解决方案．     

Imaging weak zones in the foundation using frequency domain attenuation tomography

Cross-hole imaging method using Time Domain (TD) and Frequency Domain (FD) parts of cross-hole radar tomography data acquired using Step Frequency Ground Penetrating Radar (SFGPR) was implemented. This method was adopted for imaging foundation of a dam to check if the foundation was free of geological weak zones. The dam site is characterised by massive and jointed-phyllites associated with major and minor shears. The cross-hole radar tomography data was acquired in the frequency bandwidth of 250 MHz, from the deepest level gallery up to a depth of 40 m in the foundation. In TD, first arrival time and amplitudes of radio waves were inverted using Simultaneous Iterative Reconstruction Technique (SIRT) resulting in velocity and attenuation tomograms. The tomograms showed nearly uniform velocity or attenuation structure in the respective tomographic plane. Subsequently, cross-hole radar tomography data was analysed in FD for a variation of spectrum-amplitude at different frequencies. Amplitudes picked at each single frequency were then inverted using SIRT for obtaining frequency domain attenuation tomogram (FDAT). The FDAT clearly showed presence of anomalous high attenuation zones in the depth range of 23–33 m of the tomographic plane. The anomalous zones in the attenuation tomogram are weak zones in the foundation. To validate the above observations, cross-hole seismic tomography was also done in the same boreholes. Cross-hole seismic tomography results showed low velocity (p-wave) zones around the same location corresponding to the high attenuation zone in FDAT, bringing the dormant weak zone to light. This enabled fine-tuning of the reinforcement design and strengthening the weak zone. This paper discusses the cross-hole radar tomography imaging method, the results of its application in imaging weak zones in the foundation and the comparison of cross-hole radar tomography results (in TD and FD) with the cross-hole seismic tomography results.     

预条件共轭梯度反褶积及其应用

预条件共轭梯度反褶积方法是结合稀疏反褶积的实现,运用优化的预条件共轭梯度法,完成反射系数的反演。用该方法处理地震资料时可提高资料频率,展宽有效频率宽度。但由于地震信号具有时变性,因此本文将该反褶积过程中的子波用多尺度时变子波代替。由数值算例可以看出,该方法可取得较好的实用效果。     

基于广义高斯分布的地震盲反褶积方法研究

Bussgang算法是针对褶积盲源分离问题提出的,本文将其用于地震盲反褶积处理.由于广义高斯概率密度函数具有逼近任意概率密度函数的能力,从反射系数序列的统计特征出发,引入广义高斯分布来体现反射系数序列超高斯分布特征.依据反射系数序列的统计特征和Bussgang算法原理,建立以Kullback-Leibler距离为非高斯性度量的目标函数,并导出算法中涉及到的无记忆非线性函数,最终实现了地震盲反褶积.模型试算和实际资料处理结果表明,该方法能较好地适应非最小相位系统,能够同时实现地震子波和反射系数估计,有效地提高地震资料分辨率.     

Explicit expressions and numerical calculations for the Fréchet and second derivatives in 2.5D Helmholtz equation inversion

In order to perform resistivity imaging, seismic waveform tomography or sensitivity analysis of geophysical data, the Fréchet derivatives, and even the second derivatives of the data with respect to the model parameters, may be required. We develop a practical method to compute the relevant derivatives for 2.5D resistivity and 2.5D frequency-domain acoustic velocity inversion. Both geophysical inversions entail the solution of a 2.5D Helmholtz equation. First, using differential calculus and the Green's functions of the 2.5D Helmholtz equation, we strictly formulate the explicit expressions for the Fréchet and second derivatives, then apply the finite-element method to approximate the Green's functions of an arbitrary medium. Finally, we calculate the derivatives using the expressions and the numerical solutions of the Green's functions. Two model parametrization approaches, constant-point and constant-block, are suggested and the computational efficiencies are compared. Numerical examples of the derivatives for various electrode arrays in cross-hole resistivity imaging and for cross-hole seismic surveying are demonstrated. Two synthetic experiments of resistivity and acoustic velocity imaging are used to illustrate the method.     

DIRECTIONAL DECONVOLUTION OF MARINE SEISMIC REFLECTION DATA: NORTH SEA EXAMPLE1

Directional deconvolution of the signature from a marine seismic source array may be achieved in combination with prestack migration or dip moveout (DMO) processing. The benefit is demonstrated using an example profile from the southern North Sea. In particular, shallow, dipping reflectors have improved continuity and frequency content. The method could be extended to 3D data to remove both in-line and cross-line directivity effects.     

地震资料自适应时变卡尔曼反褶积

应用卡尔曼滤波的一步预测方法,并根据卡尔曼滤波方程适用于时变系统的特点,给出了利用卡尔曼滤波进行地震资料自适应时变反褶积的方法,文中给出了理论和实际资料处理的例子.     

地震资料自适应时变卡尔曼反褶积

应用卡尔曼滤波的一步预测方法,并根据卡尔曼滤波方程适用于时变系统的特点,给出了利用卡尔曼滤波进行地震资料自适应时变反褶积的方法,文中给出了理论和实际资料处理的例子.