基于矩形网格的有限差分走时计算方法

对于大多数速度场，地震波沿射线传播的初至波走时，可以用有限差分外推的方法在二维或三维数值网格上计算出来. 在保证精度的条件下，为提高计算效率和适应性，本文推导了基于任意矩形网格和局部平面波前近似的有限差分初至波走时计算方法. 另外，该方法对首波和散射波做了合适的处理，而且不会碰到传统射线法存在的阴影区和焦散区等问题. 简单模型和复杂的Marmousi模型试算的结果表明，该方法精度较高并适用于强纵、横向变速的复杂介质. 基于该方法的Kirchhoff叠前深度偏移, 在主要构造和目的层位置的成像效果上基本达到了波动方程法叠前深度偏移的位置成像效果. 由于未考虑续至波等有效能量，在成像的保幅性上不如波动方程法叠前深度偏移的效果，但其计算效率则明显高于全格林函数法和波动方程法.      

TTI介质拟声波方程数值模拟

TTI介质qP波数值模拟方法因为考虑了倾角因素,可以比VTI介质qP波数值模拟方法更加准确地描述各向异性介质中地震波场的传播规律。文中用拟声波方程对TTI介质中的地震波场进行了高阶有限差分数值模拟,在改进衰减函数分布方式后,通过坐标变换,利用改进的完全匹配层(perfectly matched layer,PML)边界控制方程对波场边界进行吸收处理,取得了良好的效果;然后分析了拟声波方程数值模拟中的稳定性问题,并对波场中的伪横波进行压制。通过对不同模型的数值模拟,验证了文中使用的TTI介质拟声波波动方程的稳定性以及所采用的PML边界控制方程的可靠性和适用性。     

对共反射面叠加的进一步探讨

共反射面元（CRS）叠加考虑了反射层的局部特征和第一菲涅耳带内的全部反射，更充分挖掘了多次覆盖数据的潜力.但在地下介质复杂并存在倾斜层时存在反射点分散的情况，从而影响了CRS叠加效果.本文从射线理论出发，在考虑反射层局部特征的情况下，推导了水平地表和起伏地表情况下计算真实反射点分散程度的公式，最终将反射点分散程度定量表达出来.通过对反射点分散程度的控制，从CRS道集中抽取出共反射点（CRP）道集，在CRP道集中而不是在CRS道集中实现叠加，其效果应比传统的CRS叠加效果要好.利用水平地表和起伏地表的模型验证了本文所推导的公式的正确性和有效性.该公式在实际资料处理中的运用尚待进一步研究.     

双复杂条件下的波动方程叠前深度偏移

复杂地表条件下的地震勘探越来越被人们所关注.双复杂条件下的叠前深度偏移方法是解决复杂地表条件和复杂地质构造成像的有效手段."波场上延"法能实现由非水平观测界面开始的偏移过程,解决复杂地表对地下构造成像的影响.复杂理论模型的试算以及实际资料处理表明,"波场上延"方法较好地克服了起伏地形对地下构造成像的影响,取得了令人满意的效果,实现了波动方程基准面校正和深度成像的有机结合.     

安达曼海东部凹陷渐新世以来断裂–构造演化特征及其成因探讨

本文研究内容为印度洋东北部边缘海安达曼海的构造演化。利用安达曼海域东部大范围二维地震数据资料及钻井数据,结合区域地质概况以及前人研究成果,选取8条具有代表性的断层并将其划分为一级和二级断裂,运用生长指数法和古落差法对断层进行定量分析,再通过计算安达曼海东部凹陷4条主测线的构造沉降量,探讨构造演化过程。结果表明:选取的生长断层中3条属于一级断裂,跨度大,几乎切穿整个地层,属于控制安达曼海域地区沉降的大断裂;另外5条属于二级断裂,控制构造带的展布情况,属于构造带的分界线。渐新世时期,印度–澳大利亚板块与欧亚板块之间处于软碰撞阶段,断层发育缓慢,上下盘落差较小,生长指数与构造沉降量也处于低值;中新世时期,板块之间的耦合效应不断增强,断层发育速度加快,此时上下盘厚度最大,是形成多处断裂带以及多种断裂样式的关键时期,各地层生长指数和构造沉降量也达到峰值。上新世至今,安达曼海沟–弧–盆体系逐渐稳定,断层活动减弱,断裂上下盘厚度差基本一致,生长指数差异较小,构造沉降量基本稳定在 1 km 左右。     

基于有效邻域波场近似的起伏地表保幅高斯束偏移

随着我国陆上地震勘探向复杂地表探区的转移,高精度、适应性强的地震成像方法在地震资料的处理、解释及后续属性分析、储层预测中具有重要意义.本文基于有效邻域波场近似理论发展了一种成像精度更高且适用于复杂起伏地表条件的叠前保幅高斯束偏移方法.在传统水平地表高斯束偏移的基础上,本文根据中心射线附近有效邻域内高斯束表征的近似波场,导出了起伏地表条件下具有相对振幅保持的高斯束偏移公式,并给出了一种精度更高的旁轴射线传播角度计算方法.同现有的高斯束偏移方法相比,本文方法不仅考虑了起伏地表对高斯束走时的线性影响,而且首次引入了由地表高程差异和近地表速度变化引起的二次时差校正项和振幅校正项,使得成像结果更加准确可靠.两个典型模型算例验证了本文方法的正确性和有效性.     

基于CSP道集的速度分析与建模

对于低信噪比地震资料,现行的叠前偏移方法应用效果欠佳的一个重要原因就是,目前的各种速度分析方法在信噪比太低时都不能建立准确的速度模型.共散射点(CSP)道集聚集了所有可能的散射能量,具有更高的覆盖次数和更大的偏移距范围.因此,对于低信噪比资料,基于CSP道集的速度分析具有一定的优越性.本文首先简单介绍了CSP道集的形成过程和基于CSP道集的速度分析方法的基本原理,提出了一种从CMP道集到CSP道集的快速映射方法,并用一个加权系数消除了映射噪音,最后用模型和实际资料验证了基于CSP道集的速度分析与建模方法的可行性和实用性.     

基于奇异谱分析的联合去噪及规则化方法

在地震数据处理中,噪音压制和数据规则化研究直接影响后期的地震处理及解释效果.本文提出一种基于奇异谱分析的联合去噪及规则化方法,在迭代时自动在地震道缺失的位置进行插值,在含有地震数据的位置压制噪声,并通过分步插值改善数据相干性.通过合成地震记录和实际资料的联合去噪及规则化处理结果表明:联合方法能够在补全地震道的同时有效地压制噪声.     

基于含气泡液体声波方程的海底冷泉数值模拟

海底冷泉羽状流与海底天然气水合物的分布密切相关,对水合物稳定带的边界具有指示作用,是未来能源勘探的重要领域.研究海底冷泉羽状流的地震响应特征,对确定天然气水合物的储集区域及成藏环境等均有重要意义.当前获得海底冷泉羽状流的地震响应主要通过数值模拟进行,然而该过程所依据的含气泡介质声速模型及随机介质理论不能完整描述海底冷泉的物理性质,采用的声波方程也不适用于高频地震波数值模拟.为了准确地实现海底冷泉羽状流地震波数值模拟,精确分析其地震响应特征,本文提出利用Keller-Miksis气泡振动模型来描述气泡在声波作用下的运动状态,同时考虑气泡间的相互作用,建立海底冷泉气泡模型.在此基础上,本文创新性地采用含气泡液体声波方程进行海底冷泉高频地震波数值模拟.数值模拟结果表明,本文提出的方法能够实现海底冷泉羽状流地震响应的高精度数值模拟.