三维地震数据离散光滑插值的共轭梯度法

针对三维地震数据插值,提出采用Laplacian算子进行光滑约束的插值方法,并借鉴Mallet研究的离散光滑插值思路,采用预条件共轭梯度法,直接生成网格节点上的值,从而回避寻求满足插值方程的函数. 为了实现其中Laplacian算子的快速求逆,文中引入Claerbout螺旋坐标系谱因式分解理论. 在螺旋坐标系下,Laplacian算子的表示矩阵具有Toeplitz结构,其快速求逆可由谱法LU分解实现. 基于二维离散光滑插值,文中还给出共轭梯度法与NMO相结合的沿时间切片逐层处理的离散光滑插值流程. 最后,应用该方法对模型数据和实际三维地震数据进行了处理.     

单斜介质中方位NMO速度Thomsen参数反演方法研究

采用启发式共轭梯度法,即随机爬山法 + 共轭梯度法,利用单斜介质中P波方位NMO速度椭圆轴向偏转角度接近于零这一特性,简化P波方位NMO速度公式,并利用多方位P波NMO速度,反演出某一初始CMP观测线与自然坐标系之间的夹角,作为进一步进行Thomsen各向异性参数反演的基础. 根据各向异性介质中方位NMO速度与Thomsen参数之间的关系,建立了利用三种波的多方位NMO速度及垂直传播速度反演单层单斜各向异性介质Thomsen各向异性参数的目标函数. 对计算的理论值添加具有一定标准差的正态分布的随机噪声,用以模拟实际观测存在的误差,通过对加噪后的数据进行多次反演的误差分析,表明了所建立的目标函数及选用的反演方法是有效可行的,而且相对稳定.     

伸展盆地转换带分析──以松辽盆地北部为例

转换带可划分为５种类型：离散型转换带，汇聚型高凸起与低凸起转换，同向型高凸起与低凸起转换带。其中当伸展增加时，高凸起转换带可能发展为低凸起转换带。转换带的形成是由于变形地壳物质组成的不均一笥或区域张应力场的不均一性，而使伸展断层呈侧列或雁列式发生。松辽地发育种类型转换带，而以同向型高凸起低凸起转换带最为发育，推测高凸起转换带是盆岭式裂谷系的典型转换带。转换带是油气有利聚集区。     

TWO－DIMENSIONAL MATHEMATICAL MODEL FOR STUDYING EROSION AND SEDIMENTATION OF A WANDERING STRETCH IN THE LOWER YELLOW RIVER

ＴＷＯ－ＤＩＭＥＮＳＩＯＮＡＬＭＡＴＨＥＭＡＴＩＣＡＬＭＯＤＥＬＦＯＲＳＴＵＤＹＩＮＧＥＲＯＳＩＯＮＡＮＤＳＥＤＩＭＥＮＴＡＴＩＯＮＯＦＡＷＡＮＤＥＲＩＮＧＳＴＲＥＴＣＨＩＮＴＨＥＬＯＷＥＲＹＥＬＬＯＷＲＩＶＥＲ￥ＺＨＡＮＧＳｈｉｑｉｒ（Ｓｅｎｉｏｒ...     

基于平面波解构滤波的速度无关动校正及速度分析方法

为了避免速度分析拾取速度参数引起的繁复费时、拾取参数不准确导致的校正效果不佳和传统动校正的拉伸畸变问题,利用速度无关动校正方法进行相关处理。通过数据解构获取局部同相轴斜率,解析斜率数据与速度的关系,由局部斜率可直接获得动校正所需时差参数,其信息足以进行动校正处理。通过模拟和实际数据进行的实验计算,证明此方法优于常规动校正处理效果。     

Coherent Chemical Variation Trends of the 55-25 Ma Magmatic Rocks in SE Tibet:N-S Direction Lithospheric Stretching of Eurasia during Early Stage of India-Eurasia Collision

The progressive indentation of India into Eurasia generated an E–W-trending orthogonal collision belt and a N–S-trending oblique collision belt. Compiling available data reveals that ～70% of the Cenozoic igneous rocks in eastern and southeastern Tibet are concentrated within an ENE-trending, ～550-km long and ～250-km wide magmatic zone(CMZ) that once separated the orthogonal and oblique collision belts. The Latitude 26°N Line is now its southern boundary. The onset timing of magmatism of the CMZ ...     

基于曲波稀疏变换的拉伸校正方法

动校正是地震数据处理中的重要步骤,但它在校正过程中会产生子波拉伸畸变效应,随着偏移距的增大,会出现主频降低、振幅扩大的现象。由于存在拉伸畸变,同相轴未被拉平,导致非同相叠加,会引起水平叠加剖面的频率失真和分辨率下降,因此,拉伸校正是提高水平叠加剖面分辨率的关键。子波拉伸畸变在曲波稀疏域中是不相干的,可以将拉伸校正视为是一个非线性优化过程。通过度量稀疏域中数据的稀疏性,使用一种快速有效的算法,来优化子波拉伸畸变生成的非线性问题,最终实现消除子波拉伸畸变的目的。曲波稀疏变换拉伸校正方法能够消除由动校正带来的子波拉伸畸变,恢复远偏移距处的高频信息,校平同相轴。综合模型数据和实际资料处理,曲波稀疏拉伸校正方法能够显著提高水平叠加剖面的分辨率。     

Stretching correction for amplitude-preserving vector wavefield reverse-time migration

The migration of multi-wave seismic data is aimed at obtaining the P- and S-wave imaging results of the amplitude preserving. But the P- and S-wave stretching effect produced by the reverse time migration of the elastic wave equation will not only reduce the vertical resolution of the migration results and the amplitude preserving of the large reflection angle. In this paper, the reverse time migration technique of amplitude preserving vector wave-field separating is used. Based on the analysis of the stretch mechanism and the influencing factors of stretch magnitude, the paper gave the stretch correcting factors. Then, realize the stretch correction method at the time that after the reverse extrapolation and before the imaging by solving the problem which is how to calculate the P-wave and Ps-wave propagation directions of imaging points at different times. The stretch correction method can improve the vertical resolution and amplitude fidelity of the imaging results and provide high fidelity input data for seismic data interpretation and inversion.© 2019 China Geology Editorial Office.     

抛物Radon变换法近偏移距波场外推

本文给出了抛物Radon变换的基本原理,以及部分动校正后的CMP道集抛物线近似有效性的证明,基于带限正反最小平方抛物Radon变换的Levinson递推算法,对缺失的近偏移距地震波场进行叠前重建和外推.给出了抛物Radon变换法地震道重建外推的基本原理和叠前地震数据规则化的处理流程,另外对于Radon域均匀采样的情形,本文给出了均匀层状介质和Marmousi模型的近偏移距外推结果,计算结果验证了算法的稳定性和适用性.