【摘 要】
：
交错网格有限差分已经被广泛地应用到地震波数值模拟与偏移成像中.传统的交错网格有限差分空间导数的高阶差分系数一般是通过Taylor级数展开法求取的,而用该方法确定的差分系数来求解空间导数时,一般只是在一个较小的频率范围内才能取得比较高的精度.本文针对交错网格有限差分的空间一阶导数的求解,从其频散关系出发,采用最小二乘方法推导一种在更大频率范围内都能满足该导数数值解精度要求的任意偶数阶精度差分系数,然
【机 构】
：
中国科学院地质与地球物理研究所,中国科学院油气资源研究重点实验室,北京,100029;中国科学院大学,北京,100049
【出 处】
：
SPG/SEG 2014年国际地球物理会议
论文部分内容阅读
交错网格有限差分已经被广泛地应用到地震波数值模拟与偏移成像中.传统的交错网格有限差分空间导数的高阶差分系数一般是通过Taylor级数展开法求取的,而用该方法确定的差分系数来求解空间导数时,一般只是在一个较小的频率范围内才能取得比较高的精度.本文针对交错网格有限差分的空间一阶导数的求解,从其频散关系出发,采用最小二乘方法推导一种在更大频率范围内都能满足该导数数值解精度要求的任意偶数阶精度差分系数,然后分析该差分系数求解空间一阶导数的数值频散,最后采用该最小二乘交错网格有限差分法对弹性波场进行数值模拟.同时与传统的Taylor级数展开交错网格有限差分的数值频散及波场模拟的精度和计算效率进行比较.数值频散和弹性波场模拟分析表明:在相同的离散参数和几乎不改变计算效率的情况下,最小二乘交错网格有限差分能有效地压制数值频散,并比传统的Taylor级数展开交错网格有限差分具有更高的数值模拟精度.
其他文献
沥青组成的复杂性、形成机制与影响因素的多变性是固体沥青地质研究不得不面临的难题.晚奥陶世柴达木盆地北部属于活动大陆边缘,作者在联合沟滩间山群发现碳沥青,通过野外地质调查、有机地球化学分析与资源潜力计算,对碳沥青空间分布特征、地球化学特征、成因及地质意义进行了探讨.碳沥青呈北西-南东向分布在滩间山群a段的泥质岩与灰岩中.
白杨河铀铍矿床位于新疆雪米斯坦火山岩带西段,是产于我国西北地区的典型与岩浆作用关系密切的矿床.该矿床分布于区域近东西向查干陶勒盖一巴音布拉克大断裂(F1)上,其铀成矿与杨庄花岗斑岩侵入及后期岩脉活动和构造叠加关系密切,岩体侵入接触带、岩体次生断裂及走滑断裂为其主要控矿构造,其中主接触带附近为矿体赋存部位,NW向构造控制了深部矿体主要走向,后期近东西向走滑断裂系控制了地表次生铀矿富集,发育在主接触带
构造运动通过对元素扩散的制约而控制元素分布,从而影响成矿作用与矿床的分布.我国新疆西北部西准噶尔走滑断裂构造体系(简称"西准系")是中亚造山带巴尔喀什马蹄形构造的向东延伸部分,由于中生代成吉思—准噶尔断裂的右行走滑断裂作用而被分成了两个部分.西准系是一个多米诺式的走滑断裂构造体系,主要由达拉布特断裂、玛依勒断裂、巴尔鲁克断裂等三条NE走向的左行走滑断裂及其夹在它们之间的地块所组成,可能是晚古生代与
结合现代大陆动力学与地球动力学最新理论与研究和应用进展,作者重新审视和分析了构造-岩浆活化(地洼学说)的动力学机制和多因复成矿床的形成机理.综合分析认为,构造-岩浆活化的动力学机制与板块运动是相一致的,均来源于地幔(柱)对流和上涌或地幔柱-岩石圈板块的相互作用.由于板块构造具有多旋回性和地幔柱活动具有多期性或多幕式特点,据此,提出了多旋回构造-岩浆活化的概念,认为陈国达先生所建议的构造-岩浆活化区
实验样品辉长岩的斜长石含量约60％,辉石约30％,磁铁矿和钛铁矿约10％.电子探针分析表明,多数斜长石An=50-60,核部比边缘An略低,少量斜长石An=35-40,以及An=65-76;辉石中普通辉石(En35-42Fs12-19W043-46)约占70％,紫苏辉石(En58-64Fs32-39W02-3)约占25％,易变辉石(En49Fs35W016)约占5％,成分环带不明显.
不同于常规基于X-Y-Z坐标系从地表向地下进行偏移成像的射线偏移成像技术,全方位地下角度域偏移成像,通过射线追踪技术将地面的地震信息影射到地下局部角度域,每个成像点有四个极坐标分量(半开角、半开角方位角、地层倾角、地层倾角方位角),然后在地下局部角度域进行成像,可以得到两个全方位、3D角度域道集:共反射、共方向道集.通过塔里木H区块高密度、宽方位资料处理,可以看到,全方位地下角度域各向异性层析成像
在上世纪80年代提出了一系列新的识辨油气方法，其中“亮点”技术和AVO技术可用来直接识辨油气。AVO技术(或AVA)是利用振幅随偏移距变化特征分析和识辨岩性及油气藏的地震勘探技术，角度道集是AVA技术应用与分析的关键,已经广泛应用与目标成像和反演,但常规的直接转换法和直射线PSTM法角度道集都无法满足复杂地区的AVA分析与反演.本文提出了一种基于双平方根算子保幅波动方程叠前深度偏移方法,提取角度域
DB1井区属于塔里木盆地的高陡构造,针对DB1井全方位地震资料,存在地表激发条件多变引起能量差异较大、高程落差较大和低降速带变化剧烈引起近地表静校正问题严重、信噪比较差、VTI和HTI各向异性等问题,采取回转波层析反演技术很好解决近地表静校正问题,在不同域使用不同技术对不同噪声进行高保真压制,并进行地表一致性处理已达到振幅相位的一致性,考虑宽方位数据进行共偏移距矢量体上的数据规则化和插值,利用高保
相对于纵波勘探，多波多分量地震勘探能够得到更多的地下介质信息，从而减少地震勘探的多解性，因而受到了广泛的关注。但是当前多波多分量资料的处理采用的仍然是常规的纵波资料处理方法结合纵、横波分离来实现。它忽略了地震波场的矢量特征，且容易在波场分离中引入噪音，影响了多波多分量地震勘探的整体效果。因此，需要发展基于矢量波动方程的偏移成像方法，以充分利用多波多分量地震资料中的信息，发挥多波多分量勘探的优势。本
在地震勘探中，波阻抗反演是地震资料定量解释最为有效的方法。波阻抗参数作为联系地震、测井及地质信息的纽带，是进行储层区域评价不可缺少的内容。常规线性优化算法求解地震波阻抗反演时容易陷入局部极值并严重依赖初始模型,本文针对这一问题从求解目标函数的最优化问题入手,将共轭梯度法和遗传算法结合,提出了非线性混合优化算法(GA-CG).首先由多智能体遗传算法经邻域竞争、邻域交叉、变异得到优化种群,再用共轭梯度