基于井数据约束的高精度层析速度反演

测井数据是一种准确可靠的先验信息。在层析反演方程组中加入井数据信息一定程度上提高了反演精度和效率,但仍无法精确描述微幅构造地质体。针对这一特点,笔者提出了基于井数据约束的高精度层析速度反演方法。首先进行常规的井约束层析速度反演,假定测井处的网格速度不更新,建立井约束反演方程组,更新速度场,得到初步的偏移成像结果;然后在此基础上利用测井速度信息约束偏移深度,对更新后的速度场进行速度分析及精细化建模,提高速度建模精度,同时减少了反演迭代更新次数,提高了反演效率。模型试算及实际资料处理结果表明,该方法具有更高的反演精度和计算效率。     

川中古隆起超压分布与形成的地温场因素

温度和压力是沉积盆地两个重要的物理场,温度影响着超压的形成和分布.本文根据钻孔实测温度和压力数据分析了川中古隆起现今压力与温度的关系;在实验室对封闭流体进行了多组温-压关系实验;利用等效镜质体反射率和包裹体测温数据恢复了川中古隆起不同井区在白垩纪抬升之前的最大古地温,并在此基础上分析了温度降低对研究区超压的影响;最后探讨了生烃增压和欠压实超压形成过程中温度的作用.研究结果表明,川中古隆起现今超压层的压力系数与温度呈正相关关系;在绝对密封的条件下,当压力大于15 MPa时,温度每变化1 ℃,压力变化1.076 MPa.川中地区不同井区自晚白垩世以来的差异性降温是现今同一超压层系超压强度不同的主要因素,此外超压层还应发生了流体的横向压力传递和泄漏.下古生界原油裂解形成超压的时间是180~110 Ma;气态烃伴生的盐水包裹体均一温度暗示了在90 Ma超压发生调整.盆地模拟结果显示温度对上三叠统须家河组的欠压实增压影响微弱.     

基于随温度变化的热系数模拟板块俯冲动力学过程

热传导系数和热膨胀系数是影响板块俯冲动力学过程的两个重要参数. 由于地球介质的不均匀性,热系数也会随深度发生变化.然而,这种变化在地球动力学模拟研究中往往被忽略.本文针对随温度变化的热传导系数和热膨胀系数, 模拟板块俯冲的动力学过程,分析热系数、黏度对板块俯冲形态的影响及其对应的地幔对流特征.结果表明,依温度变化的热传导系数和热膨胀系数会影响地幔温度及黏度分布,进而改变板块的俯冲角度;黏度是控制板块俯冲动力学演化过程的重要因素;地幔对流受黏度结构的影响,呈现分层对流及局部多个对流环等多种不同形态的对流场特征.     

沙漠地区地震勘探随机噪声建模及其在噪声压制中的应用

随机噪声是影响地震勘探有效信号的主要因素,其存在大大降低了地震记录的信噪比.在噪声压制方法不断被改进的同时,对随机噪声特性进行研究,了解噪声的产生机制是对其进行压制的先决条件,目前对噪声的研究主要是特性研究以寻找规律性,对其进行定性定量的分析还比较少.本文根据塔里木沙漠地区实际采集环境,考虑到噪声的连续性给计算带来的不便,假设各类噪声源以点源的形式分布在检波器周围,依据相应理论确定各类噪声源的源函数,其激发的噪声经由波动方程传播,将随机噪声作为各类噪声源共同作用的综合波场,建立随机噪声的理论模型.通过分析不同种噪声对地震记录的影响,选取合适的滤波方法对其进行压制,实验结果表明,通过建立沙漠地区随机噪声的理论模型,为选择有效的滤波方法,提高地震记录信噪比起到理论指导作用.     

考虑卷积完全匹配层数值色散的井间电磁三维正演

本文将以卷积完全匹配层为吸收边界条件的时域有限差分法应用到井间电磁的三维正演模拟中. 证明了卷积完全匹配层中的数值色散会因为有效延伸因子而产生,列举常规有效延伸因子和网格间距对电磁波相速度各向异性的影响,并通过波场快照验证卷积完全匹配层中数值色散的存在;进而推广三维卷积完全匹配层中最大有效延伸因子、最大网格间距与激励源主频之间的约束,以此完善卷积完全匹配层的最优参数设置方案. 在此基础上,通过正演结果中二次场的垂直磁场分量和走时来展示灵敏度分布,以此划定井间电磁勘探的优势区域,提出井间电磁正反演所需要的约束测井系列和最佳取井方案.     

适用于声波方程数值模拟的时间-空间域隐式有限差分算子优化方法

声波方程数值模拟已广泛应用于理论地震计算,同时构成了地震逆时偏移成像技术的基础.对于有限差分法而言,在满足一定的稳定性条件时,普遍存在着因网格化而形成的数值频散效应.如何有效地缓解或压制数值频散是有限差分方法研究的关键所在.有限差分格式分为显式有限差分和隐式有限差分.隐式有限差分能够进一步压制数值频散效应.因此本文提出了给定频率范围满足时间-空间域隐式有限差分频散关系的方法,并根据震源频率、波速和网格间距确定波数范围,在此基础上建立方程确定了相应的隐式有限差分系数,使得差分系数能在更大频率范围符合波场传播规律.通过频散分析和正演模拟,验证了本文方法的有效性.     

浅层气对下伏地层深度预测影响定量研究及在渤海B油田中的应用

浅层气的存在造成下伏地层地震反射出现假象,使得深度预测不准确,增大了油田开发钻井风险。为提高浅层气下伏地层深度预测精度,利用理论模型分析量化描述浅层气对下伏地层的影响范围和影响程度,同时通过正演模拟手段构建时差定量校正量版,指导储层高精度深度预测。根据研究成果和认识,成功解释了渤海B油田两口已钻井深度误差出现的原因,并进一步利用时差校正量版提高了该油田区目标储层深度预测的精度,成功指导了后续多口开发井的准确钻探。     

台风激发的第二类地脉动特征及激发模式分析

0.003~1 Hz频段的地脉动主要来源于海浪运动与固体地球的耦合作用,台风引起的强烈海浪运动往往可使地脉动能量显著增强.由于涉及大气-海洋-固体地球三个圈层之间的复杂动量传递与耦合过程,迄今为止,关于台风激发地脉动的具体源区位置及激发机制尚存在争议.本文选取日本、中国东南沿海及台湾地区的地震台站波形连续记录,研究了2008年台风“森拉克”和“黑格比”激发地脉动的时频特征,开展相应数值模拟,并与观测数据进行了对比分析研究.结果表明台风激发第二类地脉动存在两种主要模式：（1）近岸源区激发,即台风引起波浪入射至海岸反射并与后续来波相互作用形成驻波作用于海底而激发;（2）台风中心附近源区激发,即台风中心移动过程中不同时期激发的同频率波浪相向传播、相互作用产生驻波作用于海底而激发,源区位置主要集中于台风中心左后方.此外,结合波浪再分析数据、台风风场特征,我们进一步对第二类地脉动激发过程中的影响因素进行了分析,发现：第一种模式激发的地脉动与近岸源区波浪场强度、观测点至源区距离及台风中心至海岸线距离等因素相关;而第二种模式激发的地脉动则主要受台风中心附近波浪场的频率成分与传播方向影响.

     

衰减夹层GPR模拟的时频域全局反射误差

全局反射误差分析是深入研究探地雷达（GPR）吸收边界条件吸收效率的有力工具.基于常规完全匹配层（PML）的标准交错网格有限差分算法必须满足严格的CFL条件限制,即在单位时间步长内,不容许电磁波传播的距离超过单元网格尺寸.为了提高主区域所有网格节点的计算效率,并有效地吸收传播后期出现的低频隐失波,提出基于非分裂递归卷积完全匹配层（UCPML）的旋转交错网格（RSG）GPR正演算法.该算法采用不同的网格交错策略,并在边界条件中引入了吸收低频隐失波的自由可变因子,使得该算法允许选取较大的时间步长,提高了计算效率,并且实现了对低频隐失波的高效吸收.本文首先给出了RSG差分格式,推导了基于UCPML的RSG差分更新方程,探讨了数值色散的稳定性条件,然后以绕射现象严重的衰减夹层数值模拟为例,分别从波场快照、单道波记录、时间域反射误差（TDRE）、频率域反射误差（FDRE）四个方面分析了UCPML与常规PML的全局反射误差,说明了基于UCPML和RSG的GPR正演算法能更有效地吸收低频隐失波.

     

基于Low-rank分解的复杂TI介质纯qP波正演模拟与逆时偏移

近年来,面向实际应用的TI介质准P波正演模拟与逆时偏移成像技术受到空前的关注.基于常规耦合型传播方程的正演模拟方法不仅存在伪横波及频散假象干扰,而且还遭受模型参数限制（η>0）和不稳定影响;而纯qP波方程的推导繁琐,且由于方程中包含拟微分算子造成求解难度大且精度有限.为此,本文首先构建了一种适用于任意TI介质的纯qP波传播算子,然后借助Low-rank分解求取该算子中的空间-波数域矩阵,同时引入Cerjan衰减边界条件来压制边界反射干扰,最终实现了一种间接的纯qP波波场外推方案,并将其成功应用于复杂TI介质正演模拟与逆时偏移成像中.通过开展数值模拟,并与其他方法对比表明：①该方法既避免了纯qP波方程的繁琐推导,又克服了耦合型方程对模型参数的限制;②还彻底消除了残余伪横波噪音及数值频散;③且能适应较大时间或空间步长及高频震源,是一种相对准确且稳定的各向异性纵波正演与成像策略.