简介：针对地面场源在地层介质中所产生瞬变电磁场的数值计算方法和响应分析问题，本文展开研究．解法方面，传统离散镜像法多采用复数运算、数字滤波等方法计算量大，针对这些问题，提出一种改进的离散镜像方法：基于Gaver-Stehfest概率变换算法将电磁场解式实数化，选用Prony方法对目标核函数进行指数级数逼近，根据离散镜像原理和近似系数闭合式求解瞬变电磁场．通过试算均质模型中瞬变电磁场并对比汉克尔变换的数字滤波法所得结果，证明该方法有效且具有较好的精度和适用性．继而基于该方法计算地表磁偶源在典型地电模型中产生的瞬变电磁场，对“地面激发一地层中测量”方式所得感应磁场水平分量响应进行分析并得出结论瞬变场水平分量响应与地电结构、观测时间、空间位置等因素有关，感应磁场水平分量响应反映出涡流场分布及其垂向梯度变化情况，在探测异常体的工作中应尽量选取零偏移距、较大偏移距位置钻孔或用较大观测延时以减小背景场对勘测结果的影响。文中所用离散镜像方法与正演计算结论可为相关研究工作提供参考依据。

简介：南海的构造与演化与资源环境等关系密切是本文研究重点。本文针对南海的东北部构造及其块体构造方向,利用所采集的区域地震剖面,通过解析地震相与构造及其演化的关系,提出以下观点:(1)构造分区特点明晰,可划分为五个不同构造单元,构造单元之间既有联系又相互独立;(2)南海沉积盆地无论表现为拉张-弱挤压-强挤压的何种构造格局,其区域构造应力场是统一的;(3)首次发现反射地震剖面上显示出两个浅俯冲点。每个块体构造层呈手风琴风箱式折曲并向东聚敛,体现沉积盆地从发育、成长、结束、消亡不同阶段在南海的表现,其块体俯冲方向以及块体包络区域性倾伏方向均与区域应力场方向一致。更多还原

简介：本文是三分量地震在陆相薄互层地震勘探中应用成功的首例。大庆长垣喇嘛甸油田处于高含水开发后期,储层预测的主要难题是密井网条件下厚度2m以上砂体的边界识别问题,常规单一纵波地震解释存在多解性,引入横波信息可增加地震解释的可靠性。本文根据纵、横波对气藏的响应特征,利用纵、横波联合振幅属性、分频属性对气藏进行了分析,能够很好地反应气藏的边界,体现了横波受流体影响较小的优势。通过测井曲线岩石物理定量的分析了该地区主要是密度、λρ和μρ三种物性参数能够很好地反映储层岩性变化。在以地震数据分频为核心的高精度地震参考标准层解释基础上利用纵、横波联合反演得出的纵、横波速度比、λρ、μρ等属性与测井解释的储层具有很好的一致性,多波属性λρ沿层切片描述了开发小层的砂体平面展布,与井砂岩图具有较高的一致性,刻画了井间砂体边界变化的细节,并指示了剩余油挖潜的有利区域。

简介：希尔伯特—黄变换（Hilbert-Huangtransform,HHT）是一种新的适合非平稳和非线性信号的分析方法,由于地震信号一般呈现出非平稳与非线性特性,因此HHT非常适合地震信号的分析。本文首先介绍了HHT中关于经验模态分解（EmpiricalModeDecomposition,EMD）的实现过程,在此基础上分析了几种基于EMD获得本征模态函数（IntrinsicModeFunctions,IMF）来计算瞬时频率的算法,其中利用了两个采样间隔瞬时频率的平均来计算瞬时频率,较好地反映了地震信号频率成分随时间变化的特征。将该方法应用于四川东北部某地区海相碳酸盐岩地层三维地震叠后偏移数据处理,提取＂三瞬＂地震属性,与传统的希尔伯特变换提取的＂三瞬＂地震属性进行对比,结果表明基于HHT的＂三瞬＂地震属性结果具有更高的分辨率,IMF2的瞬时相位能够较好地刻画台地边缘生物礁相,IMF2的瞬时频率亦具有较好的分带性。将IMF2的＂三瞬＂地震属性与钻井等资料结合分析,能够更好地识别沉积相的分布。

简介：在多个区块的测井评价工作发现,孔隙结构直接影响储集层的品质和油气层的电阻率,是测井准确评价流体性质的关键。岩石物理资料表明不同区块内影响储层孔隙结构的因素不同,但效果是一致的,即孔隙结构的复杂程度控制着储层的储集能力和渗透能力。孔隙结构的复杂程度影响了储层中导电流体的分布和含量,从而控制了储层的电阻率。储层出现低阻油气层的内因均为复杂的孔隙结构（骨架导电及工程原因除外）。测井储层评价在分析控制储层孔隙结构复杂程度的地质因素及储层分类的基础上,针对不同类型储层采用不同的模型、参数和标准,可以有效的认识储层品质和识别不同类型储层的流体性质。

简介：如何快速、精确地利用叠前深度偏移进行偏移速度分析是勘探地震学的一项重要研究内容,针对该问题,本文提出一种二阶精度广义非线性全局最优的偏移速度反演方法。我们将首先去掉速度模型修正量与成象深度差呈线性关系的假设,推导出具有二阶精度的速度模型修正量计算公式,使每一次迭代得到的速度模型尽可能地接近实际模型;然后采用广义非线性反演方法反演获得对所有道集的全局最优的速度模型修正量,不仅极大地加快了收敛速度,而且反演过程中陷入局部极小的可能性也减小了。理论模型和Marmousi模型的处理结果表明:本方法精度高、处理速度快,提高了偏移速度分析方法的实用性和对复杂构造成像的准确性。

简介：致密储层与围岩波阻抗差异小，常规储层反演方法难以有效预测出有利区。中国西部致密储层非均质性强，常规储层反演预测难度大，本文提出了基于Xu-White模型的叠前储层预测方法。首先修正Xu-White模型的试验参数，直至该模型计算出来的纵、横波速度及密度与实际测井得到的数值基本吻合，然后用修正的模型计算不同储层物性、岩性、含油性下的纵、横波速度、密度及泊松比，进而建立不同的砂泥岩互层地质模型。其次，采用Zoeppritz方程近似式进行叠前正演，分析不同储层物性、岩性、含油性时的叠前道集响应特征，在对中国西部Z区实际叠前道集优化处理的基础上，选择其中1种对岩性、含油性、物性最敏感的弹性参数，对其有利储层分布区进行预测，其预测结果具有较高的精度。

简介：角度域共成像点道集（ADCIGS）是偏移速度分析和振幅随角度变化分析（AVA）的基础数据。传统Kirchhoff叠前深度偏移（KPSDM）按偏移距组织数据,能方便的输出偏移距域共成像点道集（ODCIGS）,其高效的角度道集输出是有挑战的。本文提出基于旅行时梯度场的KPSDM角道集输出方法。其核心步骤为：（1）利用任意介质中的动态规划法旅行时计算方法提供炮点和检波点的旅行时场;（2）根据旅行时场的梯度方向计算反射张角;（3）在偏移过程中抽取ADCIGS。由于本文旅行时计算方法没有射线阴影区,也没有对速度光滑性的要求,其角度道集输出在阴影区比传统射线追踪更有优势。基于该角度道集输出方法,本文发展了一种适合大规模三维地震数据的KPSDM及角道集输出的并行实现方案。其基本思想是：（1）按照炮数据来组织输入数据;（2）旅行时场的输入与单炮覆盖范围相联系以节省内存;（3）多炮数据间采用MPI并行处理,单炮深度切片之间采用OpenMp并行处理,可进一步提高内存利用率和并行力度。数值试验结果证明本文角度道集生成方法的优越性和本文实现方案的有效性。

简介：应用多分量地震资料进行成像时通常需要先做波场分离,然后再对分离的波型进行成像。其中,波场分离可以在空间域或波数域实现。然而,由于用交错网格有限差分进行弹性波场数值模拟时,用来进行波数域波场分离的质点振动速度分量定义在不同网格节点上,本文提出了利用波数域插值方法来估算同一网格节点所需质点振动速度值;进而给出了先进行波数域插值后进行波场分离的波数域保幅波场分离方案。数值实验结果表明波数域插值方法具有较高的插值精度且保幅波场分离方法具有较好的保幅性,将本文方法进一步应用于弹性波逆时偏移可以获得保幅性较好的成像结果且对存在一定程度速度误差情况具有较好的适应性。

简介：基于斑块饱和模型计算渗透率变化的地震反射特征，为流体流动性的地震描述提供依据。基于传播矩阵理论设计反射系数与合成记录算法，实现了频率域岩石物理模型与地震响应计算的无缝连接。斑块饱和储层地震响应包含如下动力学信息：分界面处波阻抗差异、储层内部波的频散与衰减，以及顶底界面波的调谐与干涉。模拟结果表明，渗透率的增加显著降低纵波速度，使其在高、低频弹性极限之间发生频散。储层速度频散与层状构造共同导致反射系数的频变现象。在储层与围岩波阻抗接近的情况下，地震响应对渗透率变化具有敏感性，对于不同储层厚度，当围岩为高速页岩时，反射波叠加振幅随渗透率增加而增加；当围岩为低速页岩时，叠加振幅随渗透率增加而降低。

简介：本文从两种类型的地震检波器(10hz动圈式速度检波器和压电加速度检波器)及其检测介质的两种运动动参数(速度和加速度)入手,测试和对比分析了它们的频率响应函数,指出了二者的差异。又对两种检波器进行了冲击振动试验和结果对比,分析了其响应信号的特征和携带信息的能力。并于某地区在可比条件下进行了地震数据采集对比试验,对用两种检波器采集得到的单炮资料和叠加时间剖面进行了对比分析。结果表明,加速度信号更能满足当前和今后地震勘探对地震信号的高信噪比、高精度、高分辨率和大信息量的要求。

简介：标量CSAMT只适合一维及测量方向与构造方向垂直的二维情况,对于复杂的三维地电结构,CSAMT需采用张量测量。本文试图采用矢量有限元法实现三维张量CSAMT的正演模拟。为了验证算法的正确性,本文在层状介质中计算了三维CSAMT远区的电场,磁场及阻抗张量,并且与层状介质中的理论解进行了比较,接着还模拟了均匀半空间中含有三维异常体的模型,并且分析了四个阻抗张量、视电阻率及阻抗相位的响应特征。得出如下结论：采用矢量有限元法来模拟三维张量CSAMT,其电磁场及阻抗张量的实虚部计算精度都比较高,并且该方法本身满足电场法向不连续,不用进行散度校正。

简介：在直流电法有限元数值模拟中，针对常规结构化网格源点附近网格节点数值精度低的问题，并考虑计算效率和反演成像的需求，本文在二维点源问题中提出一种新的网格加密一收缩方法。其核心思想是在结构化网格计算区域中先后引入网格节点加密和网格节点收缩两个环节。通过在计算区域水平方向上均匀加密网格节点密度，降低源点奇异性的影响范围，并提升对地形的模拟能力；通过在计算区域垂向方向上将多个网格收缩到一个网格中，降低网格节点的规模，进而提高数值计算的效率。理论模型检验表明，网格加密-收缩方法能有效地解决源的奇异性问题，与均匀加密网格相比，网格节点总数降低了约80％。

简介：基于曾新吾和MacBeth提出的时域内双源累积旋转方法,建立了频率相关介质中分析多分量VSP数据的横波分裂参数提取算法（DCTF）。该算法可以在频域中针对单个频率提取横波各向异性参数（快横波的极化方向以及快、慢横波间的时间差）,从而避免了目前常用方法中使用窄带通滤波可能带来的误差。通过对地震合成记录的数值分析,确定了该算法的可行性和正确性,并与目前常用方法的应用结果进行了比较。结果显示,频率相关横波分裂参数可以利用DCTF从地震四分量数据中直接提取。在地震频率范围内,含较大尺度裂缝时各向异性参数将表现出频率相关性,这意味着在地震频率范围会出现频散。随着频率的增加,各向异性有降低的趋势。

简介：传统上,时间域航空电磁数据通过拟合迭代反演计算得到大地模型,然而,由于航空电磁数据道间的较强相关性,导致病态反演,并引起超定问题;同时电磁数据的相关性使其与模型参数的映射关系复杂,增加了反演的复杂度。采用主成分分析法将航空电磁数据变换为正交的较少数量的主成分,不仅降低了数据道间的相关性,减小了数据量,同时压制了数据的不相关噪声。本文利用人工神经网络(ANN)逼近主成分与大地模型参数间的映射关系,避免了传统反演算法中雅克比矩阵的复杂计算。层状模型的主成分神经网络与数据神经网络的反演结果对比显示,主成分神经网络反演方法网络结构简单,训练步数少,反演结果好,特别是对于含噪数据。准二维模型的主成分ANN、数据ANN以及Zhody方法的反演结果显示了主成分神经网络具有更接近真实模型的反演效果,进一步证明了主成分神经网络反演方法适合海量航空电磁探测数据反演。

简介：常规欧拉反褶积法中构造指数的选取以及分散解存在较多的问题,提出了基于联立垂向一阶导数与解析信号的欧拉齐次方程的RDAS-Euler反演方法。该方法可以更为精确的估计场源的范围及埋深,且不需考虑构造指数N的影响,避免了因构造指数不当而引起的反演误差。通过对单一地质体及组合地质体模型的实验证明本文方法能有效地完成目标体的反演工作,反演结果与理论值之间的误差小于10%,且相对于常规欧拉反褶积法更加稳定准确,能够更好的得到地质体边界及深度信息。将RDAS-Euler法应用于黑龙江省虎林盆地实测布格重力异常数据,获得了丰富的断裂信息,说明RDASEuler法增强了对断裂平面位置的识别能力。

简介：具有较高自然频率的高灵敏度检波器在采集数据时可以通过压制低频信号来相对提高高频能量，但这也造成了其低频响应差的问题，如果将高灵敏度数据与常规数据进行优势组合，就可以达到拓宽频带的目的。为此，本文提出分频段匹配滤波的方法，即在保持常规数据低频优势的前提下，对其高频端进行匹配滤波，实现不同频带范围内的优势互补，从而改善地震记录。通过引入不同主频的雷克子波模拟得到具有常规数据和高灵敏度数据特点的理论模型，论证了分频匹配滤波方法的可行性。在对野外单炮地震记录处理中发现，分频匹配滤波方法拓宽了地震记录的有效频带宽度，提高了地震记录的分辨率。

简介：地震波场数值模拟中不可避免地会出现边界反射,一般采用吸收边界条件以压制人工边界反射。目前常用的分裂式完全匹配层（PML）边界条件需要在边界处进行特殊处理,尤其是在三维情况下需要将变量分裂为三个分量,增加了数值模拟的计算时间和内存占用量。与分裂式PML吸收边界条件相比,混合吸收边界条件（HABC）具有易于实现、计算量小和吸收效果好等优点,可以提高三维波动方程数值模拟的计算效率。本文将基于一阶Higdon单程波方程的混合吸收边界条件从二维计算域发展到三维,提出了适用于三维弹性波数值模拟的混合吸收边界条件。均匀模型以及复杂模型的三维数值模拟结果表明,混合吸收边界条件与传统的完全匹配层边界条件相比,具有效率高、吸收效果好的优势。

简介：角度域共成像点道集(ADCIGS)是偏移速度分析和振幅随角度变化分析(AVA)的基础数据。传统Kirchhoff叠前深度偏移(KPSDM)按偏移距组织数据,能方便的输出偏移距域共成像点道集(ODCIGS),其高效的角度道集输出是有挑战的。本文提出基于旅行时梯度场的KPSDM角道集输出方法。其核心步骤为:(1)利用任意介质中的动态规划法旅行时计算方法提供炮点和检波点的旅行时场;(2)根据旅行时场的梯度方向计算反射张角;(3)在偏移过程中抽取ADCIGS。由于本文旅行时计算方法没有射线阴影区,也没有对速度光滑性的要求,其角度道集输出在阴影区比传统射线追踪更有优势。基于该角度道集输出方法,本文发展了一种适合大规模三维地震数据的KPSDM及角道集输出的并行实现方案。其基本思想是:(1)按照炮数据来组织输入数据;(2)旅行时场的输入与单炮覆盖范围相联系以节省内存;(3)多炮数据间采用MPI并行处理,单炮深度切片之间采用OpenMp并行处理,可进一步提高内存利用率和并行力度。数值试验结果证明本文角度道集生成方法的优越性和本文实现方案的有效性。

简介：本文提出了一种基于Kirchhoff积分偏移和逆时偏移的联合速度分析及成像方法,采用剩余曲率分析及层剥离策略进行偏移速度建模。本文方法改善了Kirchhoff积分偏移在复杂构造时计算精度不高和逆时偏移计算效率慢的缺点,兼有计算效率高、成像精度高的优点,并将其应用在反射波法隧道超前预报中。通过模型试算,发现隧道中使用逆时偏移的成像结果在多方面优于Kirchhoff积分偏移的成像结果;通过对实测数据的处理,验证本方法计算效率较高,建立的速度模型合理,成像剖面清晰,结合地质调绘资料可对隧道开挖前方地质构造做出较准确预报。