地震记录小波域高阶相关叠加技术

水平叠加技术能够有效地压制随机噪声和多次波,但是对于由于相邻道的噪声表现出相关性和一致性而产生的某些相干干扰,常规的叠加技术容易产生假的同相轴,难以取得较高的信噪比.基于小波理论和高阶统计理论本文提出了小波域高阶相关叠加技术(HOCS),其基本原理是将经动静校正的CMP道集变换到小波域,然后在小波域求取高阶相关系数,进而加权叠加.数值实验和实际资料计算的结果表明,此方法比常规叠加更能够有效地提高信噪比,剔除常规加权叠加无法抑制的某些相关噪声.     

基于斜率分析的倾角域绕射波分离及成像研究（英文）

碳酸盐岩缝洞型储层所含缝、洞尺度差异大、岩性和速度变化剧烈,造成复杂的绕射响应特征。部分缝洞体尺度小产生的绕射能量较弱,在成像剖面上被连续反射层掩盖,而难以识别。基于此,本文发展了一种倾角域道集的绕射波成像方法,提高对成像剖面上小尺度绕射构造的分辨能力。借助高斯束偏移技术抽取倾角域成像点道集,根据在倾角域道集上绕射同相轴拟线性、反射能量拟抛物的几何形态差异,基于斜率分析方法实现波场分离,并采用汉宁窗函数提高绕射波分离效果,将分离后的绕射波倾角域道集水平叠加得到绕射波成像结果。理论分析和数值模型测试结果表明:在窗函数的基础上,应用基于斜率分析的波场分离方法得到的绕射波倾角域道集能够更好的压制噪声,叠加后的绕射波成像剖面中绕射构造更加清晰。     

TI介质局部角度域射线追踪与叠前深度偏移成像

研究与实践表明,对于长偏移距、宽方位地震数据,忽略各向异性会明显降低成像质量,影响储层预测与描述的精度.针对典型的横向各向同性(TI)介质,本文面向深度域构造成像与偏移速度分析的需要,研究基于射线理论的局部角度域叠前深度偏移成像方法.它除了像传统Kirchhoff叠前深度偏移那样输出成像剖面和炮检距域的共成像点道集,还遵循地震波在成像点处的局部方向特征、基于扩展的脉冲响应叠加原理获得入射角度域和照明角度域的成像结果.为了方便快捷地实现TI介质射线走时与局部角度信息的计算,文中讨论和对比了两种改进的射线追踪方法:一种采用从经典各向异性介质射线方程演变而来的由相速度表征的简便形式;另一种采用由对称轴垂直的TI(即VTI)介质声学近似qP波波动方程推导出来的射线方程.文中通过坐标旋转将其扩展到了对称轴倾斜的TI(即TTI)介质.国际上通用的理论模型合成数据偏移试验表明,本文方法既适用于复杂构造成像,又可为TI介质深度域偏移速度分析与模型建立提供高效的偏移引擎.     

基于delta波包叠加的时间域深度偏移

delta波包可由高斯波束经傅里叶逆变换得到,是高斯波束在时空域的表达.它最早出现在合成理论地震图的研究中,本文将其应用于偏移领域.通过delta波包叠加表达时间域格林函数,可将高斯波束偏移由频率域转换到时间域,再结合Rayleigh积分和激励时间成像条件,本文给出了基于delta波包叠加的深度偏移算法.该偏移算法可在时间域直接计算,但因包含褶积运算,成像时将耗费大量的计算时间.针对这一问题,本文提出了将褶积简化为乘积的近似公式.近似后的偏移算法,不仅保留了高斯波束偏移的优点,而且计算效率得到显著提升.文中通过两个数值算例验证了上述结论.     

频率域全波形反演方法研究进展

全波形反演方法利用叠前地震波场的运动学和动力学信息重建地下速度结构,具有揭示复杂地质背景下构造与岩性细节信息的潜力.根据研究需要,全波形反演既可在时间域也可在频率域实现.频率域相对于时间域反演具有计算高效、数据选择灵活等优势.近十几年来频率域全波形反演理论在波场模拟方法、反演频率选择策略、目标函数设置方式、震源子波处理方式、梯度预处理方法等方面取得了进展.目标函数存在大量局部极值的特性是影响反射地震全波形反演效果的重要内在因素之一.如果将Laplace域波形反演、频率域阻尼波场反演、频率域波形反演三种方法有机结合,可以降低反演的非线性程度.     

角度域叠前深度偏移叠加算子优化

叠前深度偏移是解决精细复杂构造成像的一种最佳方法,尤其是基于角度域叠前深度偏移成像方法的研究,由于射线路径的唯一性,被认为是一种能够消除假象,解决射线多路径问题的最优化方法．文章从单程波波动方程出发,推导了中点 偏移距域双平方根波动方程波场延拓算子;利用局部偏移距与反射角的关系,推导了角度域共成像点道集的构建方法．同时,为有效解决叠前深度偏移计算量大、耗时长等计算瓶颈,引入OpenMpi并行策略,利用共享存储环境下的多线程协同作业,提高了叠前深度偏移算子的计算效率;在角度域,通过调整改进叠加算子,一定程度上减少了计算量,提高了计算效率．最后,通过理论模型验证了角度域叠前偏移成像算子的正确性,评价了并行算法的计算效率及有效性;并将其应用于SEG/EAGE岩丘模型,利用“沉降观测”理论得到了不同CMP的中点 偏移距域共成像点道集和角度域共成像点道集,对比分析了两种共成像点道集的优缺点;并通过限制成像叠加角度的方法,改善了复杂模型的成像结果,验证了角度域叠前深度偏移成像算法对复杂模型的适应性．     

表驱三维角度域Kirchhoff叠前时间偏移成像方法

针对传统Kirchhoff叠前时间偏移方法的一些不足,以及振幅随入射角、方位角变化(AVA/AVAZ)分析的需要,本文提出一种基于射线理论的三维叠前时间偏移角度域成像方法.它通过横向均匀介质中稳健的射线追踪建立单程波走时和传播角度的数值表,然后在此基础上估算反射波双程走时以及在界面处传播的方位角和入射角,最后基于脉冲响应叠加原理获得三维构造图像和方位\|角度域共成像点道集.与传统方法不同之处在于,上述过程均考虑了地震波在垂向变速介质中的射线弯曲效应和三维传播特征,有利于准确提取随入射角和方位角变化的振幅和时差信息.理论模型合成数据和实际地震资料测试结果展示了方法的优越性与实用性.     

数据域初至波走时与成像域反射波走时联合层析速度建模方法

复杂地表探区,尤其是盆山过渡区的油气勘探是我国也是世界上油气勘探的重点区域,但是此类区域油气地震勘探中满足精确地震成像的速度建模一直是个没有很好解决的问题.本文提出了一种综合性的数据域初至波走时与成像域反射波走时联合层析复杂地表浅中深层速度建模方法,并针对联合层析速度反演解的非唯一性问题,深入地分析了层析反演中正则化的本质意义,指出了建立构造特征正则化方法的具体技术路线,提出了联合层析的实现流程及策略.理论和实际数据试验表明,本文提出的数据域初至走时与成像域反射走时联合层析浅中深层速度建模技术避免了常规建模方法中浅层速度模型与中深层速度模型的融合问题,较好地解决了传统成像域反射层析对近地表模型的不可控更新问题,整体提升了深度域浅中深层速度模型的建模精度,进而提高了复杂地表、复杂构造区的地震成像质量.     

基于重力数据特征的东北地区深部动力机制解析

中国东北地区处于古亚洲构造域、蒙古—鄂霍茨克构造域和环太平洋构造域叠加作用最为显著的地区,是地学研究的热点区域.为了探析欧亚大陆下西太平洋板片的俯冲形态以及揭示该区域深部地球动力学机制,利用卫星重力数据通过预处理共轭梯度快速密度反演算法获得了包含东北地区、华北部分地区及日本海海域在内的研究区域上地幔三维密度结构,结合天然地震三维层析成像结果共同揭示太平洋板片的俯冲形态和深部动力机制.俯冲的太平洋板片在日本海沟处呈高密度异常,向西横向持续扩张,深度方向上有逐渐增加趋势.不连续的高密度体俯冲至地幔转换带(410～660km)后继续水平西向俯冲,部分滞留板片可能进入下地幔;在大兴安岭断裂带下面转换带中同样发现水平分布的高密度体,推断是大兴安岭断裂带下方地幔岩石圈拆沉的结果,横向不均匀分布的俯冲板片边缘已抵至大兴安岭造山带附近,这对于研究东北地区深部动力学机制具有重要的意义.     

大兴安岭西盆地群域构造与地球物理场综述

中国东北地区大兴安岭西侧盆地群包括漠河盆地、根河盆地、拉布达林盆地、海拉尔盆地和二连盆地等,蕴藏着丰富的中、新生代油气资源.为研究该盆地群域古生代、中新生代构造演化,综合建立盆地群域地球动力学模型,补充东北亚构造演化理论,本文综述该盆地群域受控的区域构造与深部构造背景、盆地群构造特征与性质、主要控盆断裂特征、盆地群油气条件比较以及盆地群域已完成并取得重要结果的地球物理工作.归纳已有主要认识和研究结果:(1)对大兴安岭西侧的盆地群起构造控制作用的构造带包括蒙古—鄂霍茨克洋缝合带、西拉木伦河缝合带、黑河—贺根山缝合带、塔原—喜桂图缝合带、西太平洋板块俯冲带,以及额尔古纳—呼伦断裂和得尔布干断裂.(2)二连盆地、海拉尔盆地和漠河盆地的盆地构造轴向与蒙古—鄂霍茨克洋缝合带走向相关;而且三个盆地内的一级构造单元走向(隆起、坳陷和推覆带)也具有这类特点.(3)几个地学断面的综合地球物理研究表明,大兴安岭西侧盆地群岩石圈地幔厚度自北向南变厚,南部盆地基底与华北地台基底表现类似;盆地群基底电性结构因受到软流圈热物质作用可能在继续演化.(4)在盆地沉积地层方面,漠河盆地的下部是侏罗系陆相煤系地层,上部是白垩系火山岩地层;海拉尔盆地由下侏罗统的铜钵庙组、南屯组,上侏罗统的大磨拐河组和下白垩统的伊敏组共同组成扎赉诺尔群,厚约3000m;二连盆地中生代地层中,中下侏罗统主要为含煤建造,上侏罗统为火山岩建造,下白垩统主要为含油建造和含煤建造,上白垩统为砂砾岩建造.(5)盆地群整体勘探程度较低.基于上述研究结果,需要进一步研究的科学问题包括:由本研究区的地球物理、构造地质、石油地质等多学科的综合研究,解决研究区受控的区域构造应力场所包括的因素及其作用,以及在岩石圈尺度上三维空间的地球物理场表征;深部构造对盆地群域构造的作用;从晚古生代到中新生代研究区构造演化特点及其依据;从北至南约1650km长的盆地群域构造差异与依据;盆地群(域)油气条件与毗邻的松辽盆地在构造成因上的差异.     

柯西波场及在角度域成像条件中的应用

角度域成像道集是叠前深度偏移的重要输出结果,它是偏移速度分析、各向异性分析和AVA分析等研究工作的基础.目前存在的角度域成像道集的生成方法受计算效率或角度分辨率的影响,仍然满足不了实际生产的要求.角度域成像道集的生成方法可以大致分为直接法和间接法两大类.在直接方法中,波矢量方向计算和局部平面波分解是两个最重要的内容,它们共同决定角度域成像条件的实现效率和角度域成像道集的质量.为了完善现有的角度域成像道集生成方法,本文提出一种新的波矢量计算方法和局部平面波分解方法.本文先用波动方程任意时刻的柯西条件构造一个只含原波场负频率成分的柯西波场,然后根据柯西波场在时间波数域的振幅谱计算波场的波矢量方向.该方法仅在需要计算波矢量方向的时刻合成柯西波场,不需要增加额外的数据读写操作,是一种高效的波矢量计算方法.本文还以柯西波场为基础提出一种高效的局部平面波分解方法,保证角度域成像条件的实现效率.结合柯西波场和局部平面波分解方法,本文最后给出一种新的角度域成像方法.文中最后的数值实验证明该方法得到的角度域成像道集具有理想的角度分辨率,可以反映地下构造的角度照明情况.     

Wheeler域储层地震预测技术研究

利用常规地震数据进行层序划分、沉积相识别、沉积演化分析及储层预测,因受现今构造趋势的影响,具有较强的多解性.Wheeler域地震资料因具有等时性,地层旋回性清楚、展布范围明了,井震标定后便于层序及体系域划分,易于识别砂体,在进行地层的沉积演化过程分析,可以在相控约束下分析砂体的分布特征、成因,评价有利储层,减少多解性.本文通过正演探讨了同一沉积模式,不同地震频率下的地震成像,研究了不同沉积模式下的Wheeler转换方法,实现了地震资料从常规时间域转换到Wheeler域.通过Wheeler域井震标定,建立起地震地质的联系;拓展性的提出了在地震最小研究单元内,在双域(常规时间域和Wheeler域)进行解释和综合研究的技术流程;建立了基于地震沉积学的储层地震预测技术.并在准噶尔盆地西缘白垩系储层研究进行了应用.经钻井证实,这种基于地震沉积学的Wheeler域储层地震预测技术比较准确,效果良好,值得完善和推广.     

时间域航空电磁扩散特征和成像深度研究

为了研究时间域航空电磁系统扩散特征,基于连续性边界条件将电磁场向下延拓得到地下介质中各点处的频率域响应,通过傅里叶变换将其变换到时间域.计算阶跃波形发射时垂直磁偶极子和水平磁偶极子在均匀半空间介质和两层介质中产生的电磁场,通过电流密度矢量图和随时间变化的电流密度等值线展示电磁场在地下介质中的扩散过程.电磁场在地下的扩散受电阻率结构的影响.在良导地区电磁场扩散慢、衰减快,而在高阻地区电磁场扩散快、衰减慢.垂直磁偶极子在地下产生的感应电流形成一个电流环随时间向下、向外扩散.水平磁偶极子在地下产生的感应电流形成两个互相叠加的电流环.时间域电磁场扩散代表着真正物理意义上的电磁扩散,因为它展示了地下电磁场空间分布随时间的变化和场强的衰减.基于对感应电流环的研究,我们探究了扩散深度和时间域航空电磁数据成像深度的关系.时间域电磁场扩散规律的研究不仅加深了对电磁扩散特征的理解,更有助于提高航空电磁数据解释水平.     

祁连山西段疏勒河流域地貌特征及其构造意义

疏勒河流域盆地位于祁连山西段,跨越了该区多条不同方向和不同性质的活动断裂带,其所呈现的地貌特征反映了该区最新构造活动的信息。文中基于GIS空间分析技术,利用SRTM-3数字高程模型(DEM)数据,系统提取了疏勒河流域及其4个亚流域盆地的面积-高程积分曲线、疏勒河水系的Hack剖面以及河流坡降指标(SL),并对整个流域地形做了坡谱分析,获得了疏勒河流域的地貌特征。研究表明,疏勒河的流域地貌发育受到该区的阿尔金断裂、昌马断裂、托勒南山断裂以及疏勒南山断裂等强烈构造活动和区域岩性差异的影响,坡度以0°~35°的区段为主;整个流域盆地处于河流发育的"壮年期";构造活动是造成河流纵剖面发生改变的最主要因素,局部河段同时还受到岩性因素的控制和影响。这表明在活跃的造山带内部,河流地貌的发育过程中,活动断裂的构造作用是重要的控制因素。     

环域大地逆边值问题

从确定大地水准面实际出发,提出环域大地逆边值问题．文中首先建立环域大地道边值问题的数理模型．由于环域内边界待定,属自由边界,本质上是非线性问题．循传统给出环域逆边值问题的线性化形式．重点讨论并构造了线性化问题的解式,包括谱域内的解．     

基于改进NAD方法的频率域声波逆时偏移

逆时偏移成像具有高分辨率的特点,在频率域实施能提升偏移成像效率.数值求解波动方程是逆时偏移成像的关键环节,数值算法的优劣直接决定着偏移过程的计算效率与成像结果的质量.本文在频率域近似解析离散化(NAD)方法的基础上通过优化网格差分模板系数构造改进NAD方法,保证数值离散精度的同时可缓解阻抗矩阵的病态程度,以此提高频率域...     

TI介质局部角度域高斯束叠前深度偏移成像

各向异性射线理论基础上的局部角度域叠前深度偏移方法能够为深度域构造成像与基于角道集的层析反演提供有力支撑,但是对于复杂地质构造而言,高斯度叠前深度偏移在不失高效、灵活等特点的情况下,具有明显的精度优势.为此,本文研究局部角度域理论框架下的高斯束叠前深度偏移方法.为提高算法效率与实用性,文中讨论了一种从经典弹性参数表征的各向异性介质运动学和动力学射线方程演变而来的由相速度表征的简便形式,并提出了一种比较经济的各向异性高斯束近似合成方案.结合地震波局部角度域成像原理,讨论一种适合高斯束偏移的角度参数计算方法.国际上通用的理论模型合成数据试验表明：相比局部角度域Kirchhoff叠前深度偏移成像方法,本文方法具有更高的成像精度与抗噪能力,既适用于复杂构造成像,也可为TI介质深度域偏移速度分析与模型建立提供高效的偏移引擎.     

基于广义模型约束的时间域航空电磁反演研究

由于航空电磁具有海量数据,因此快速有效的成像和反演手段至关重要.本文针对层状介质模型推导与实现了广义模型约束条件下时间域航空电磁一维反演.从正则化反演的目标函数出发,通过改变模型约束项构造Lp范数反演和聚焦反演,进而通过改变模型求解域构造出基于小波变换的稀疏约束反演.针对不同反演方法目标函数的构建方式,本文进一步从数学原理上分析不同反演方法的预期效果,并通过理论模型和实测数据进行验证.结果表明L0.8范数反演、聚焦反演和基于小波变换的稀疏约束反演可以得到更符合地下层状介质陡变界面的反演结果.     

基于Gabor-Daubechies小波束叠前深度偏移的角度域共成像道集

传统炮检距域共像集（CIG）在复杂介质中因波传播的多路径而存在反射体位置不确定的问题. 角度域CIG由于克服了这一缺陷而逐步成为速度分析、AVA以及振幅保真偏移成像等研究的主要手段. 以波动理论为基础的地震偏移成像方法的发展为获得高质量的角度域CIG提供了可靠的实现途径. 其中,基于波场局域化分解和传播的小波束域波场延拓和偏移成像方法,因其波场分解基本函数和传播算子在空间和方向上的双重局域特性,而成为角度相关分析研究的有效工具. 本文在采用Gabor Daubechies框架分解的小波束叠前角度域偏移成像基础上,利用不同的叠加方法由局部角度域像矩阵得到了反射角域CIG（CRAIG）和倾角域CIG（CDAIG）. 以SEG EAGE二维盐体模型为例,通过对CRAIG和CDAIG的对比,探讨了这两种角度域CIG的特点及其在地震偏移成像中的潜在应用.     

非均匀变形域中先存构造活动性的判定

Anderson模式只适用于均匀弹性变形域,无法解释自然界中复杂的断层作用现象,实际应用中有很大的局限性.本文结合砂箱实验成果和库仑-摩尔准则,分析基底先存构造的活动性,提出了判定先存构造活动性的力学参数--活动性系数(faS),它由先存构造的产状和力学性质以及应力状态(主应力的大小和方向)等方面决定.该参数考虑到地质...