简介:针对地面场源在地层介质中所产生瞬变电磁场的数值计算方法和响应分析问题,本文展开研究.解法方面,传统离散镜像法多采用复数运算、数字滤波等方法计算量大,针对这些问题,提出一种改进的离散镜像方法:基于Gaver-Stehfest概率变换算法将电磁场解式实数化,选用Prony方法对目标核函数进行指数级数逼近,根据离散镜像原理和近似系数闭合式求解瞬变电磁场.通过试算均质模型中瞬变电磁场并对比汉克尔变换的数字滤波法所得结果,证明该方法有效且具有较好的精度和适用性.继而基于该方法计算地表磁偶源在典型地电模型中产生的瞬变电磁场,对“地面激发一地层中测量”方式所得感应磁场水平分量响应进行分析并得出结论瞬变场水平分量响应与地电结构、观测时间、空间位置等因素有关,感应磁场水平分量响应反映出涡流场分布及其垂向梯度变化情况,在探测异常体的工作中应尽量选取零偏移距、较大偏移距位置钻孔或用较大观测延时以减小背景场对勘测结果的影响。文中所用离散镜像方法与正演计算结论可为相关研究工作提供参考依据。
简介:本文给出一种既能有效衰减地震噪音又可保护地层及构造的不连续性的新方法。构造约束保边平滑技术需要已知反射局部方位和边界信息,通常这些信息由全频率地震资料估算获得,但在资料信噪比很低的情况下,噪音往往会降低估算的可靠度。对于信噪比极低的地震资料,其主频成分相对非主频成分信噪比高,所以由主频资料获取的方位和边界信息比由其它频率成分获取的更可靠。方位和边界信息通常用倾角和相干值差异来描述。由于不同频率所引起的倾角和相干值差异的变化均比地震记录的变化缓慢,所以由主频资料获取的倾角及边界信息能够近似代表所有频率成分的倾角及边界信息。Ricker子波广泛用于地震勘探,Marr小波与Ricker子波在时间和频率域均具有相同的形态,所以选用Marrl小波变换将地震数据按照倍频程分为几个分频体。扫描主频分频体,用不等权二次曲面拟合并求解极大值来获取视倾角,通过比较9个滑动窗口的相干值来确定反射边界。将这些信息用构造约束保边平滑技术可选择性地(selectively)对主频、低频、高频分频体做平滑处理,最后将平滑后的各频段地震记录合成为滤波去噪后的地震记录。理论模型和实际资料处理效果表明该方法能有效压制噪音,保护边界,保护同相轴的连续性,且灵活地保留地震记录中的有用信息。
简介:南海的构造与演化与资源环境等关系密切是本文研究重点。本文针对南海的东北部构造及其块体构造方向,利用所采集的区域地震剖面,通过解析地震相与构造及其演化的关系,提出以下观点:(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并行处理,可进一步提高内存利用率和并行力度。数值试验结果证明本文角度道集生成方法的优越性和本文实现方案的有效性。
简介:应用多分量地震资料进行成像时通常需要先做波场分离,然后再对分离的波型进行成像。其中,波场分离可以在空间域或波数域实现。然而,由于用交错网格有限差分进行弹性波场数值模拟时,用来进行波数域波场分离的质点振动速度分量定义在不同网格节点上,本文提出了利用波数域插值方法来估算同一网格节点所需质点振动速度值;进而给出了先进行波数域插值后进行波场分离的波数域保幅波场分离方案。数值实验结果表明波数域插值方法具有较高的插值精度且保幅波场分离方法具有较好的保幅性,将本文方法进一步应用于弹性波逆时偏移可以获得保幅性较好的成像结果且对存在一定程度速度误差情况具有较好的适应性。
简介:常规的时间一空间域和频率一空间域预测滤波方法假设地震记录由地震信号和随机噪声两部分构成,即所谓的加噪声模型,但是,在对随机噪声进行估算时,又假设随机噪声可以通过预测误差滤波器由地震记录中进行预测,即所谓的源噪声模型。这种前后不一致的噪声模型降低了该类方法的去噪能力和保幅性能。为此,本文提出了一种基于反演的时空域随机噪声衰减方法。它首先从地震数据中估算预测滤波算子,该算子表征了地震信号的可预测性,自适应地描述了地震信号的空间结构。在得到预测误差算子之后,将该算子作为正则化约束引入到地震信号反演系统,由含有随机噪声的地震数据直接反演地震信号。不同于常规随机噪声衰减方法,该方法将随机噪声衰减问题归结为正则化约束下的地震信号反演问题,克服了常规方法噪声模型的不一致性问题。我们采用模型数据和实际数据进行了实验分析,并与常规方法进行了效果对比。实验结果表明:与常规方法相比,本文方法在噪声压制的同时,没有对有效信号产生明显伤害,具有更好的振幅保持能力。
简介:标量CSAMT只适合一维及测量方向与构造方向垂直的二维情况,对于复杂的三维地电结构,CSAMT需采用张量测量。本文试图采用矢量有限元法实现三维张量CSAMT的正演模拟。为了验证算法的正确性,本文在层状介质中计算了三维CSAMT远区的电场,磁场及阻抗张量,并且与层状介质中的理论解进行了比较,接着还模拟了均匀半空间中含有三维异常体的模型,并且分析了四个阻抗张量、视电阻率及阻抗相位的响应特征。得出如下结论:采用矢量有限元法来模拟三维张量CSAMT,其电磁场及阻抗张量的实虚部计算精度都比较高,并且该方法本身满足电场法向不连续,不用进行散度校正。
简介:在直流电法有限元数值模拟中,针对常规结构化网格源点附近网格节点数值精度低的问题,并考虑计算效率和反演成像的需求,本文在二维点源问题中提出一种新的网格加密一收缩方法。其核心思想是在结构化网格计算区域中先后引入网格节点加密和网格节点收缩两个环节。通过在计算区域水平方向上均匀加密网格节点密度,降低源点奇异性的影响范围,并提升对地形的模拟能力;通过在计算区域垂向方向上将多个网格收缩到一个网格中,降低网格节点的规模,进而提高数值计算的效率。理论模型检验表明,网格加密-收缩方法能有效地解决源的奇异性问题,与均匀加密网格相比,网格节点总数降低了约80%。
简介:基于曾新吾和MacBeth提出的时域内双源累积旋转方法,建立了频率相关介质中分析多分量VSP数据的横波分裂参数提取算法(DCTF)。该算法可以在频域中针对单个频率提取横波各向异性参数(快横波的极化方向以及快、慢横波间的时间差),从而避免了目前常用方法中使用窄带通滤波可能带来的误差。通过对地震合成记录的数值分析,确定了该算法的可行性和正确性,并与目前常用方法的应用结果进行了比较。结果显示,频率相关横波分裂参数可以利用DCTF从地震四分量数据中直接提取。在地震频率范围内,含较大尺度裂缝时各向异性参数将表现出频率相关性,这意味着在地震频率范围会出现频散。随着频率的增加,各向异性有降低的趋势。
简介:提高地震资料的分辨率是地震数据处理流程中的重要环节,对后续的精细构造解释起到重要作用。传统的提高分辨率方法大都假设地震资料是稳态的并且噪声水平不随空间发生变化,而实际情况不满足这一假设,导致提高分辨率处理后的效果达不到预期要求。针对这一问题,本文提出了一种基于时频二次谱的提高地震资料分辨率方法。首先,文中提出了基于S变换的时频二次谱,并结合模型论述了时变子波和反射系数在时频二次谱中的特征及其可分离性;其次,依据时变子波和反射系数在时频二次谱中的特征差异,构建了二维滤波器在地震记录的时频二次谱中提取时变子波的振幅谱;再次,文中研究了噪声环境中时变提高分辨率算子设计方法,并提出了依据时频谱能量强弱相对关系自适应确定频带拓宽范围的时变提高分辨率算子设计,进行提高分辨率处理;最后,文中对该方法进行了模型和实际数据的试处理,并与传统谱模拟方法和Q补偿方法的处理结果进行了对比分析,对比结果表明:本方法不需要估计Q值,提高分辨率能力不受震源子波频带的限制,在兼顾信噪比的前提下能够充分提高不同时间局部的地震数据的分辨率。
简介:传统上,时间域航空电磁数据通过拟合迭代反演计算得到大地模型,然而,由于航空电磁数据道间的较强相关性,导致病态反演,并引起超定问题;同时电磁数据的相关性使其与模型参数的映射关系复杂,增加了反演的复杂度。采用主成分分析法将航空电磁数据变换为正交的较少数量的主成分,不仅降低了数据道间的相关性,减小了数据量,同时压制了数据的不相关噪声。本文利用人工神经网络(ANN)逼近主成分与大地模型参数间的映射关系,避免了传统反演算法中雅克比矩阵的复杂计算。层状模型的主成分神经网络与数据神经网络的反演结果对比显示,主成分神经网络反演方法网络结构简单,训练步数少,反演结果好,特别是对于含噪数据。准二维模型的主成分ANN、数据ANN以及Zhody方法的反演结果显示了主成分神经网络具有更接近真实模型的反演效果,进一步证明了主成分神经网络反演方法适合海量航空电磁探测数据反演。
简介:常规欧拉反褶积法中构造指数的选取以及分散解存在较多的问题,提出了基于联立垂向一阶导数与解析信号的欧拉齐次方程的RDAS-Euler反演方法。该方法可以更为精确的估计场源的范围及埋深,且不需考虑构造指数N的影响,避免了因构造指数不当而引起的反演误差。通过对单一地质体及组合地质体模型的实验证明本文方法能有效地完成目标体的反演工作,反演结果与理论值之间的误差小于10%,且相对于常规欧拉反褶积法更加稳定准确,能够更好的得到地质体边界及深度信息。将RDAS-Euler法应用于黑龙江省虎林盆地实测布格重力异常数据,获得了丰富的断裂信息,说明RDASEuler法增强了对断裂平面位置的识别能力。