识别断裂构造的DCT法

为了提高欧拉反褶积(Euler法)的反演精度,提出基于DCT的Euler法确定断层断点位置的反演方法,从理论上给出了基于DCT的Euler齐次方程,指出异常导数的计算精度是影响Euler法反演精度的核心问题.台阶模型实验证实:用DCT法计算的重力异常垂向一阶导数的最大误差为0.460×10-9/s2,均方差为0.189×10-9/s2,水平导数最大误差为0.182×10-9/s2,均方差为0.028×10-9/s2,与DFT法相比,其计算的异常导数具有很高的精度.通过对正演模型的计算,确定了基于DCT的Euler法模型反演的位置为台阶断面的中心点位,与理论点位相对应;在反演模型实验中,证实了DCT法具有较高的反演精度,同时在实际应用中,进一步验证了该法的可行性、实用性和准确性.     

鄂尔多斯南缘地区重力变化场源特征

以鄂尔多斯南缘地区布格重力异常数据及2014—2017年相对重力观测数据为基础,采用欧拉反褶积方法,对引起重力变化的场源深度进行反演,并对空间分布规律予以解释。通过构建理论模型,反演得到最优模型参数,并对实际数据进行计算和分析。为了减弱和消除欧拉解的发散性,利用水平梯度滤波法优化反演结果。结果表明:①构造指数为1时,适合对相对重力数据进行反演,当滑动窗口选择4—8倍测点间距时,可以获得较为可靠的场源参数;②在主要断裂附近,不同时间尺度的场源呈较好的一致性特征;③引起重力变化的场源深度集中在10—30 km,与鄂尔多斯南缘地区的震源深度基本一致;④不同时间尺度的场源位置相对分散,表明引起重力变化的物质流动具有随机性和波动性。     

反褶积试井解释方法在东海油气田中的应用

海上油气井由于测试时间短,利用常规试井解释方法往往难以较准确反映地层边界特征,而应用反褶积算法进行试井解释,基于叠加原理截取更长的压力生产史参与计算,能够得到更长时间的压力响应,从而获取更多的储层信息,为储层边界的识别提供可靠的依据,对不稳定试井解释技术是一个很好的补充和完善.尝试利用东海有限的油气井测试资料利用反褶积方法进行解释和分析,并对其应用进行归纳和总结.这对东海的勘探开发、对反褶积方法本身发展均有较重要的参考价值.     

Prestack nonstationary deconvolution based on variable-step sampling in the radial trace domain

The conventional nonstationary convolutional model assumes that the seismic signal is recorded at normal incidence. Raw shot gathers are far from this assumption because of the effects of offsets. Because of such problems, we propose a novel prestack nonstationary deconvolution approach. We introduce the radial trace （RT） transform to the nonstationary deconvolution, we estimate the nonstationary deconvolution factor with hyperbolic smoothing based on variable-step sampling （VSS） in the RT domain, and we obtain the high-resolution prestack nonstationary deconvolution data. The RT transform maps the shot record from the offset and traveltime coordinates to those of apparent velocity and traveltime. The ray paths of the traces in the RT better satisfy the assumptions of the convolutional model. The proposed method combines the advantages of stationary deconvolution and inverse Q filtering, without prior information for Q. The nonstationary deconvolution in the RT domain is more suitable than that in the space-time （XT） domain for prestack data because it is the generalized extension of normal incidence. Tests with synthetic and real data demonstrate that the proposed method is more effective in compensating for large-offset and deep data.     

Research on 3D marine electromagnetic interferometry with synthetic sources for suppressing the airwave interference

In order to suppress the airwave noise in marine controlled-source electromagnetic （CSEM） data, we propose a 3D deconvolution （3DD） interferometry method with a synthetic aperture source and obtain the relative anomaly coefficient （RAC） of the EM field reflection responses to show the degree for suppressing the airwave. We analyze the potential of the proposed method for suppressing the airwave, and compare the proposed method with traditional methods in their effectiveness. A method to select synthetic source length is derived and the effect of the water depth on RAC is examined via numerical simulations. The results suggest that 3DD interferometry method with a synthetic source can effectively suppress the airwave and enhance the potential of marine CSEM to hydrocarbon exploration.     

利用多次波提高地震数据分辨率方法

目前对多次波的有效利用仅围绕多次波成像技术展开,通过成像多次波试图获取更丰富的地下构造信息.不同于该思路,本文另辟蹊径,从利用多次波提高地震数据分辨率角度出发,对多次波的有效利用进行了深入挖掘.首先基于聚焦变换思想在聚焦域内实现多次波的降阶,通过理论推导得出聚焦域内多次波表现为原始数据的多维子波反褶积这一重要结论,从理论上证明了本文方法提高地震数据分辨率的可行性;然后采用引入整形正则化的非稳态回归自适应匹配滤波方法将聚焦域内由多次波构建的高分辨率数据分离出来,实现原始数据的高分辨率转换.与常规反褶积模型不同,该方法基于波动理论推导得出,可以适用于任意复杂情况;每一道输出结果中所有炮记录都参与了运算,从空间上加以约束,在提高纵向分辨率的同时可以改善数据的横向分辨率.最后通过模型试算和实际资料处理对本文方法的有效性、适应性和实用性进行了验证.     

潮汕坳陷多次波组合压制技术研究

潮汕坳陷中生界发育,是我国海上重要的油气勘探后备区。但该区多次波非常发育,多次波识别和压制一直是该区地震资料处理的关键和难点。合理地选取有效的多次波压制方法是获得高品质地震剖面的关键所在。针对原始资料中多次波异常发育的情况,文章分析了该区多种类型多次波的特点,尝试了零偏移距截距时间-射线参数(τ-p)域预测反褶积、自由界面多次波衰减法(Surface-Related Multiple Elimination,SRME)和双曲线Radon变换,对多次波进行组合压制,克服了仅依靠一种手段难以压制多种类型多次波的局限。从资料处理效果看,τ-p域预测反褶积方法对短周期的鸣震和微曲多次波的压制效果较好,SRME有效地压制自由表面多次波,且在近炮检距范围内压制多次波效果很好,而双曲线Radon变换是衰减长周期的层间多次波和未被充分压制的自由表面多次波很好的选择。组合压制多次波后的剖面质量明显得到改善。     

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

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